时间域广义2.5D一阶波动方程曲网格有限差分法数值模拟

2.5D地震波场数值模拟通过在2D地质模型中施加点源从而计算得到3D地震波场。首先通过傅里叶变换得到了一种适用于混合(声波、弹性各向同性、弹性各向异性)介质和各种边界条件(声波自由地表、固体自由地表和固—液边界)的2.5D时域广义一阶波动方程,并采用曲线网格有限差分法求解该波动方程。在各种均匀介质模型(声波、弹性各向同性和弹性各向异性)中,通过对比2.5D数值解与3D解析解和3D数值解,不仅验证了推导的方程和数值求解方法的正确性,而且验证了2.5D数值方法相比3D数值方法在计算效率和内存占用方面有很大的优势。2D数值方法由于线源假设,其解与2.5D数值解相比存在较大的振幅误差和相移,难以直接应用。数值实验结果表明,该2.5D数值模拟方法适用于含各种边界(声波自由地表、固体自由地表和固—液边界)的地质模型。不同于2D波场数值模拟方法,2.5D波场数值模拟方法可直接应用于实际的点源观测数据处理,如2.5D逆时偏移成像。

TTI介质声波方程分裂式PML吸收边界条件研究

针对倾斜横向各向同性(TTI)介质声波波动方程的特点,研究了TTI介质分裂式完全匹配层(Split perfectly matched layer,SPML)吸收边界条件。首先对常见的几种TTI介质声波波动方程进行了归纳,并从TTI介质波场传播稳定性的角度进行对比分析,结果表明,引入横波分量的TTI介质纵横波耦合方程适用于TTI介质。然后从TTI介质纵横波耦合二阶波动方程出发,推导得到其一阶波动方程的形式,进而推导出一阶波动方程形式的SPML波动方程,并给出了高阶交错网格有限差分算法的具体实现过程。数值模拟结果表明,SPML吸收边界条件能达到很好的人工边界反射吸收效果,相比优化海绵吸收边界条件,其人工边界反射吸收效果更好。

非均匀介质地震波传播交错网格高阶有限差分法模拟

采用常规的二阶声波方程有限差分方法对于非均匀介质进行了数值模拟时 ,其数值模拟精度较低。而采用一阶双曲型标量波动方程 ,则无须对介质的弹性常数进行空间求导。根据Taylor级数展开式 ,推导出了交错网格一阶空间导数的任意偶数阶精度展开式和相应差分系数计算式以及一阶双曲型标量波动方程交错网格任意偶数阶精度差分格式 ,并给出了该差分算法的稳定性条件。用该差分算法对均匀介质模型、非均匀介质模型和Marmousi模型进行了数值模拟试验 ,并与伪谱法进行了对比。结果表明 ,一阶双曲型标量波动方程交错网格高阶差分法的模拟精度与伪谱法的精度非常接近 ,计算效率高 ,且适合于模拟非均匀介质。

TI介质地震波场数值模拟边界条件处理

地震波场数值模拟研究中,边界条件处理非常重要,它直接影响模拟地震数据的精度.针对TI(Transversely Isotropic)介质一阶速度—应力弹性波动方程组,利用时间域分裂法给出了PML(Perfectly matched layer,完全匹配层)吸收边界条件,其基本思想是在所研究的区域外引入吸收层,地震波从有效计算区域传播到吸收层时不发生任何反射,且在吸收层内按传播距离呈指数衰减,不发生反射.利用交错网格高阶差分技术对PML吸收边界条件求解,对比不同边界条件数值模拟结果表明,PML吸收边界条件较应用广泛的Clayton吸收边界有更好的吸收效果,能更好地满足地震数值模拟的需要.

Numerical modeling of seismic wavefields in transversely isotropic media with a compact staggered-grid finite difference scheme

To deal with the numerical dispersion problem, by combining the staggeredgrid technology with the compact finite difference scheme, we derive a compact staggered- grid finite difference scheme from the first-order velocity-stress wave equations for the transversely isotropic media. Comparing the principal truncation error terms of the compact staggered-grid finite difference scheme, the staggered-grid finite difference scheme, and the compact finite difference scheme, we analyze the approximation accuracy of these three schemes using Fourier analysis. Finally, seismic wave numerical simulation in transversely isotropic （VTI） media is performed using the three schemes. The results indicate that the compact staggered-grid finite difference scheme has the smallest truncation error, the highest accuracy, and the weakest numerical dispersion among the three schemes. In summary, the numerical modeling shows the validity of the compact staggered-grid finite difference scheme.