多轴完全匹配层的非分裂实现

作者：李振春; 田坤; 黄建平; 曹晓莉; 李娜; 李庆洋

来源：地球物理学进展, 2013, (06): 2984-2992.

• 卷积完全匹配层
• 多轴完全匹配层
• 吸收边界条件
• 数值模拟
• 不分裂 convolution PML
• multiaxial PML
• absorbing boundary condition
• Forward modeling
• unsplit

完全匹配层吸收边界(PML)已经被证明是非常有效的边界吸收技术,对体波和面波的吸收都具有非常好的效果,已经被广泛应用于弹性波的数值模拟中.但是在某些情况下传统的PML和C-PML(卷积完全匹配层)还是存在一定的问题,比如在一些介质中有固有的不稳定性等等.基于分裂的多轴完全匹配层(M-PML)通过引入多个衰减剖面可以有效地消除不稳定的现象.本文采用C-PML中不分裂的递推卷积方法,对M-PML进行了实现,这样在实现过程中不用分裂变量,与原来的分裂方法相比,直观易懂,方便简单,易于编程,而且可以减少计算量和存储量,节约计算成本.基于一阶速度-应力方程,采用交错网格高阶有限差分法对不同的介质模型进行了数值模拟,结果证明了本文的不分裂M-PML可以消除不稳定性,得到比较好的计算结果.

• 出版日期2013-12-15
• 单位中国石油大学（华东）