慢度法全波场模拟及其水平慢度积分的实现

叠前地震数据的反演需要以某种正演算法为基础。本文从弹性动力学方程出发,采用慢度法计算层状介质的叠前地震记录。同反射率法一样,慢度法也能够提供地震波场的完全解,包括地震波在地层中传播时形成的一次波、多次波和转换波信息,而且其中间结果容易解释。具体实现步骤为:①通过傅立叶变换和汉克尔交换将本构方程和弹性动力学方程变换到频率-慢度域;②在频率-慢度域得到反射、透射系数矩阵递推公式;③分别对频率和慢度进行积分变换得到时-空域的地震道集。文中还讨论了慢度积分计算中解的稳定性问题。