声波方程频率域高精度正演的17点格式及数值实现

频率域正演计算是频率域全波形反演的基础.传统的最优9点格式只具有二阶精度,不能满足高精度地震成像的需要.本文考虑两个四阶精度的格式,即经典的四阶9点格式和优化的17点格式.17点格式可将最小波长内所需网格点数减小到2.56.通过在简单模型和Overthrust模型上的数值实验,比较分析了三种格式的正演效果;简单模型数值实验显示了17点格式克服频散误差的能力优于四阶9点格式和最优9点格式;复杂模型数值实验则进一步承认了算法的可行性.

基于MCPML边界条件的频率域声波方程高精度模拟

由于地下介质结构复杂,混合网格二阶精度的9点差分已不能满足频率域声波方程数值模拟的需要。在四阶精度9点差分格式的基础上,推导了四阶精度17点差分格式,并引入非分裂的多轴卷积完全匹配层(MCPML)边界条件。相比传统完全匹配层(PML)边界条件,MCPML结合了卷积完全匹配层(CPML)边界条件和多轴完全匹配层(MPML)边界条件的优点,在不增加计算时间的同时改进了对大入射角波场的吸收效果,而且非分裂的方法易于编程实现。应用Marmousi模型进行试算,验证了方法的有效性。

基于LSQR算法的二维声波方程频率域正演模拟与数值实现

频率域正演是频率域全波形反演的基础,高效准确地完成频率域正演计算是目前该领域的一个热门问题。这里利用优化17点差分格式和PML完全匹配层构建波场,并使用LSQR算法求解波场,通过在简单层状模型和Overthrust模型上计算验证可知,LSQR算法可大幅提高频率域正演的计算速度和模拟精度,从而为此方向的研究提供一个可行的技术手段。