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

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