[发明专利]全波形反演中去除波动方程模拟直达波的方法及装置

权利要求书

1.一种全波形反演中去除波动方程模拟直达波的方法，其特征在于，包括以下步骤：

输入速度模型；

基于声波波动方程对所述速度模型进行波动方程正演，得到模拟的地震记录；

计算直达波到达地表接收器的时间T₀(x)，x为检波点位置到震源位置的距离；

根据公式计算输入时间窗口T(x)，α为控制所述输入时间窗口的拉伸系数；

删除所述地表接收器接收到的其接收时间小于所述T(x)时刻的地震记录。

2.根据权利要求1所述的全波形反演中去除波动方程模拟直达波的方法，其特征在于，所述速度模型为二维速度模型。

3.根据权利要求1所述的全波形反演中去除波动方程模拟直达波的方法，其特征在于，所述声波波动方程为双程波动方程∂vx∂t=∂P∂x∂vz∂t=∂P∂z∂P∂t=υ2∂vx∂x+υ2∂vz∂z+s(x0,z0,t),]]>其中，s(x₀,z₀,t)为震源，P为应力向量，v_x和v_z分别为模拟地震剖面P_rec(t,x)、时间-空间波场P(t,x,z)，其采用改进的非分裂PML吸收边界条件，即

∂vx∂t=∂P∂x-Ωxx∂vz∂t=∂P∂z-Ωzz∂P∂t=υ2∂vx∂x+υ2∂vz∂z-υ2Ψxx-υ2Ψzz,---(1)]]>

其中Ω_xx、Ω_zz和Ψ_xx、Ψ_zz为引入的中间辅助变量，它们对应的辅助方程为：

∂Ωxx∂t+σxΩxx=σx∂Txx∂x∂Ωzz∂t+σzΩzz=σz∂Tzz∂z∂Ψxx∂t+σxΨxx=σx∂Vx∂x∂Ψzz∂t+σzΨzz=σz∂Vz∂z,---(2)]]>

联立（1）式和（2）式可得到PML波动方程问题的求解。