CPML条件下黏弹介质地震波传播的有限差分模拟与特征分析

在波动方程有限差分波场数值模拟中,为了使计算得到的边界波场值更接近于真实的边界波场值,需要使用边界条件以减少来自计算区域边界的人为反射能量。传统的分裂式完全匹配层(SPML)吸收边界条件不能有效吸收掠射波,后来发展的不分裂卷积完全匹配层(CPML)能够较好地处理这个问题,并且CPML在处理边界问题时无需对波场进行非物理分裂。本文针对Kelvin-Voigt模型的黏弹性介质,采用高阶交错网格有限差分方法进行数值模拟,同时利用CPML吸收边界条件进行边界处理。数值模拟结果表明:与弹性波相比,大角度入射产生的低频掠射波对黏弹性波的影响更大;CPML吸收边界条件对低频掠射波的吸收效果比传统的SPML吸收边界条件更好。