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

详细信息   

• 责任者：李振春^c中国石油大学(华东)^d1963^f教授 ; 田坤 ; 黄建平 ; 曹晓莉 ; 李娜 ; 李庆洋
• 刊名：地球物理学进展
• 出版年：2013
• 卷期：28^b6总122
• 页码：2984～2992
• 索书号：P260.6/208-1
• 图表：^c参38插图
• 文摘期号：2014/06

摘要

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