[发明专利]一种基于稀疏反演的层间多次波压制方法

权利要求书

1.一种基于稀疏反演的层间多次波压制方法，其包括以下步骤：

1）在叠加剖面上选取介于强反射界面之间的基准面，在炮集上分离来自基准面上地层的波场数据δP和下地层的波场数据

2）拾取δP初至一次波振幅峰值对应的脉冲估计基准面以上地层脉冲响应的初始值，并通过匹配滤波估计震源子波初始值，并设定反演脉冲响应与震源子波的迭代次数；

3）迭代更新基准面以上地层的脉冲响应、基准面以下地层的脉冲响应以及震源子波；

4）匹配层间多次波预测结果与原始地震数据得到振幅校正因子，利用振幅校正因子更新脉冲响应和震源子波；

5）以步骤4）的参数更新结果作为新的初始值，继续重复步骤3）～4）直至达到给定迭代次数，停止迭代；

6）从地震数据P中减去预测得到的层间多次波M得到一次波P₀。

2.如权利要求1所述的一种基于稀疏反演的层间多次波压制方法，其特征在于：所述步骤3）迭代更新基准面以上地层的脉冲响应、基准面以下地层的脉冲响应以及震源子波的具体计算过程为：

3.1）更新基准面以上地层的脉冲响应δX：

利用最速下降法得到δX在i+1次迭代时的最速下降方向：

ΔδX=VSiH-P‾VHX‾i]]>

变换ΔδX到时间域，拾取每道振幅峰值对应的脉冲后变换至频率域，得到脉冲响应的更新方向ΔδX，乘以不依赖于频率的步长α，更新脉冲响应与拟合残差：

δX_i+1＝δX_i+αΔδX

V＝V+αΔV₁

其中，

α=-Σωvec(ΔV1)Hvec(V)Σωvec(ΔV1)Hvec(ΔV1)]]>

ΔV1=-ΔδX‾Si+X‾iΔδX‾HP‾;]]>

3.2）更新基准面以下地层的脉冲响应

利用最速下降法计算在i+1次迭代时的最速下降方向：

ΔX‾=VSiH-VP‾HδXi+1]]>

同样的，变换到时间域，拾取每道最大振幅值对应的脉冲后变回频率域得到的更新方向乘以不依赖于频率的更新步长β，更新脉冲响应与拟合残差：

X‾i+1=X‾i+βΔX‾‾]]>

V＝V+βΔV₂

其中，

β=-Σωvec(ΔV2)Hvec(V)Σωvec(ΔV2)Hvec(ΔV2)]]>

ΔV2=-ΔX‾‾Si+ΔX‾‾δXi+1HP‾;]]>

3.3）匹配脉冲响应与更新震源子波。