【摘 要】
：
本文提出利用NAD方法数值离散频率域声波方程，在计算区域的边界处考虑完美匹配层（PML）吸收边界条件，最终导致一个复值线性方程组。接下来分析了这个线性方程组系数矩阵的稀疏结构
【机 构】
：
清华大学数学科学系 北京 100084
论文部分内容阅读
本文提出利用NAD方法数值离散频率域声波方程，在计算区域的边界处考虑完美匹配层（PML）吸收边界条件，最终导致一个复值线性方程组。接下来分析了这个线性方程组系数矩阵的稀疏结构特点以及特征值性质，将所要求的复线性方程组同解地转化为一个实值分块二乘二线性方程组，然后采用不精确旋转分块三角预处理子（IRBTP）加速Krylov子空间方法（例如：GMRES方法）来求解这个实线性方程组。数值试验表明：该方法在计算时间方面明显优于现有一些求解方法，从而也就加速了频率域正演计算过程。最后进行在各种不同介质中的波场模拟，通过与其它数值策略比较，体现了NAD方法在压制数值频散和节约计算方面具有明显的优势。
其他文献
三维地震勘探是一种高精度、多信息量的勘探技术,在实际生产中得到越来越广泛的应用.在地震资料处理中,正演模拟的精度直接影响到反演结果的准确性,复杂构造地区三维资料的处
全波形反演方法是一种经典的数据域波动方程层析方法.本文针对声介质波动方程反射波旅行时反演，指出该方法目前存在的两个问题：偏移成像结果的不准确性对反偏移生成的反射波记
全波形反演能够利用叠前地震波场的全波形信息,不仅包括旅行时等运动学信息,还包括相位、振幅等等动力学信息,因此反演结果具有很高的精度;但其作为一个典型的地球物理反问题
在地震勘探中,射线理论做为波动理论的一种高频近似,以其灵活高效的特点在地震波正演与成像研究中得到了广泛的应用.射线追踪作为射线理论的核心,通常包括运动学射线追踪和动
在地震勘探中,地震数据往往是多维信号,按单道进行处理往往会忽略剖面横向上的特征,所以在这里本文探讨一下2D-EMD的实现过程.在去噪应用中,本文将一个原始地震剖面分解成若
随着当前地震勘探形式的发展,基于各向同性介质的叠前深度偏移和速度建模已经不能满足对地下复杂构造准确描述的要求.随着勘探精度的逐步提高,各向异性的影响不能忽略,准确地
地震逆时偏移(Whitmore,1983)是目前精度最高的成像方法,是一种基于波动理论的深度域偏移方法.逆时偏移是将检波点接收到的波场进行逆时延拓和震源波场正向延拓,然后利用成像
One of long-standing problems with underwater seismic studies is the inaccurate timing because the internal clock of Ocean Bottom Seismograph(OBS)is unable to s
初至波走时层析成像是近地表速度结构建模的重要方法.传统的三维走时层析反演在应用中存在诸多问题:一是射线追踪技术固有的计算效率低、对复杂模型计算不稳定;二是对于大规
在基于波动方程的地震反演及偏移中,正演计算是最为关键的因素之一.而目前正演计算的主流方法中,差分方法凭借简单快速易于并行的优势,得到了广泛应用.然而,差分方法也同时具