基于耦合反射/透射系数单程波传播算子的地震波模拟研究

单程波近似实际上是一种多次前向散射和单次后向散射近似.利用单程波近似来描述波传播可以极大地节省地震数值模拟的计算时间和内存，实现地震波长距离传播模拟和三维地震模拟快速计算.本文基于单程波近似和波动积分方程的分离变量逼近,从广义Lippmann-Schwinger波动积分方程推导出耦合反射/透射系数的单程波传播算子.该算子由两部分构成：分离变量Fourier单程波传播算子和薄板间的反射/透射系数表达.前者将常规的Fourier分裂步单程波传播算子(SSF)推广适应横向强速度变化介质和大角度传播波场.后者是利用垂直波数来表示反射/透射系数，自然耦合到波场传播的计算过程中，其为地质界面倾角的隐式表达，精确描述振幅随入射角的变化，能适应任意复杂的模型.通过两个数值算例和一个实际地质模型的计算，本文将该方法和边界元法进行了比较，结果表明：在算例给出的介质横向速度变化情况下，本文提出的方法在相位和振幅方面与全波数值方法基本吻合.