摘要

通常地震波声波方程正演模拟所采用的方程是描述波在均匀介质中传播的波动方程,该方程用于非均质介质时,其模拟结果存在较大误差,因此需要在界面处引入过渡层,即认为弹性参数是连续变化的。采用有限元方法求解非均质声波方程,将单元内的弹性参数插值,使介质上任意一点的位移和应力始终是一个连续函数,这样对于突变界面就能通过过渡层来近似模拟,当该层单元格足够小时可以有效控制误差,提高正演模拟精度。最终通过数值模型算例验证了该方法的有效性。