[发明专利]一种基于散射波场级数法的层间多次波压制方法

权利要求书

1.一种基于散射波场级数法的层间多次波压制方法，其特征在于，包括依次进行的以下步骤：

S1、层间多次波源头分析

在叠加数据体上拾取产生层间多次波的层位，确定层间多次波源头的时窗范围；

S2、数据准备

对偏移后或动校正后的叠前道集进行一维化处理；

S3、多次波模型预测

数据驱动与层位及时窗相关的多次波模型预测；

S4、自适应相减

含多次波的数据与多次波模型自适应相减得到多次波压制后的数据；

S5、叠前道集恢复，然后进行叠加

把多次波压制后的数据重新分选形成新的叠前道集，然后进行叠加，最终形成不含多次波的地震剖面；

所述新的叠前道集为不含多次波的道集。

2.根据权利要求1所述的一种基于散射波场级数法的层间多次波压制方法，其特征在于，所述步骤S2按照以下步骤顺序进行：

S21、从全工区叠前道集中抽取单次剖面；

S22、将各个单次剖面抽取出来形成n个叠后数据，完成叠前道集一维化的处理；

步骤S21中选取n个单次剖面的过程选自方法一或方法二，其中，

方法一

炮检距域叠前时间偏移后道集，偏移距按一定增量规则递增，选取相同的偏移距组成一个单次剖面，则全工区可获得n个单次剖面，这时n取值为最大偏移距或偏移距分组增量；

方法二

OVT域叠前时间偏移后道集，一个OVT子集就是全工区的一个单次剖面，则全工区按OVT子集分选可获得n个单次剖面，这时n取值为最大覆盖次数；

所述步骤S5中，把多次波压制后的数据按不同的偏移距或者OVT子集将它们对应的单次剖面重新分选形成新的叠前道集，然后进行叠加，最终形成不含多次波的地震剖面。

3.根据权利要求2所述的一种基于散射波场级数法的层间多次波压制方法，其特征在于，所述步骤S3按照以下步骤顺序进行：

S31、假设从地表测得的散射波场称为地震数据D＝(φ_s)_ms，根据D的阶数展开扰动计算V，可得

V＝V₁+V₂+V₃+……，          式①；

S32、将式①和散射波方程结合起来，得到一阶、二阶和三阶逆散射级数：

D＝G₀V₁G₀        式②

G₀V₂G₀＝-G₀V₁G₀V₁G₀          式③

G₀V₃G₀＝-G₀V₁G₀V₂G₀-G₀V₂G₀V₁G₀-G₀V₁G₀V₁G₀V₁G₀    式④；

S33、通过式②，由散射场和背景场格林函数来推导出V₁；

再由式③和V₁，推导出V₂，以此类推，完成对逆散射级数的求解；

其中，三阶逆散射级数的第三项式结合背景介质的格林函数公式：

以及单频平面波场b₁：

且满足形成层间多次波的限制条件z₁z₂且z₂z₃，即可得到层间多次波场的预测算法。