PML边界条件下二维粘弹性介质波场模拟

在波场模拟中较多使用的有限差分法和有限元法都或多或少地存在一些缺陷,如为了提高计算精度而使运算效率降低,在泊松比变化大的界面存在稳定性问题等。伪谱法是一种计算精度高且计算效率也较高的方法,它是对空间坐标通过快速傅里叶变换在时间域作差分运算,避免了求偏导数;不存在有限差分法和有限元法对高频成分限制的问题,可以实现全频带地震波模拟;对内存容量的需求也远远低于有限差分法和有限元法。为此,采用了伪谱法进行粘弹性介质波场模拟。首先给出了开尔芬粘弹性介质中的波动方程,推导出二维粘弹性介质最佳匹配层(PML)吸收边界条件及其相应的伪谱法计算公式;然后对均匀介质模型进行了模拟,结果表明,该方法的边界吸收效果很好;最后通过数值模拟分析了具有不同粘滞系数介质对地震波的吸收和衰减,结果表明,随着粘滞系数的增大,粘弹性反射波的主频向低频方向移动,高频吸收明显,有效频带变窄,振幅降低。此外,在伪谱法中引入匹配层边界条件后,只需对空间做一维傅里叶变换,大大提高了伪谱法的计算效率。