[发明专利]陆地地震勘探资料的降维自适应层间多次波压制方法

权利要求书

1.一种陆地地震勘探资料的降维自适应层间多次波压制方法，其特征是，所述层间多次波压制方法在叠前CMP道集通过采用虚同相轴方法，实现陆地地震资料的层间多次波预测和压制；包括输入数据准备过程、构建参考道过程和压制多次波过程；具体步骤如下：

A.输入数据准备：

A1.对原始地震数据通过常规地震初步处理，得到准确动校正后的道集；

A2.根据数据道头中的CMP号和偏移距，将上述动校正后的道集重排至CMP-偏移距域；

B.参考道构建：对叠前CMP道集做加权叠加，得到高信噪比的参考道；

C.采用虚同相轴方法压制多次波：自上而下依次选取每一层地层w，执行步骤C1～C4：

C1.进行数据切割：记来自地层w的反射波旅行时为t_w，将步骤B所述参考道切割成仅包含来自地层w的反射波的数据R_upper(t)和旅行时大于t_w的数据R_lower(t)；对叠前CMP道集做相同切割，得到仅包含来自地层w的反射波的数据G_upper(t)和旅行时大于t_w的数据G_lower(t)；

C2.在叠前CMP道集中，利用R_upper(t)与每道的G_lower(t)进行互相关，构建虚同相轴数据G_vir(t)；

C3.利用G_vir(t)褶积G_lower(t)，重构得到与地层w有关的层间多次波，作为预测得到的多次波数据；

C4.将步骤C3所述预测得到的多次波数据与叠前CMP道集中的多次波进行匹配，将经匹配后的多次波从叠前CMP道集中减去，即完成该地层的层间多次波的压制；

C5.自上而下依次选取第1至n层地下反射层，执行C1-C4的步骤，完成全部地层的层间多次波的压制。

2.如权利要求1所述层间多次波压制方法，其特征是，步骤A2所述将上述动校正后的道集重排至CMP-偏移距域，具体是按照CMP号从小到大排列，相同CMP号的道集按偏移距从小到大排列。

3.如权利要求1所述层间多次波压制方法，其特征是，步骤C4所述匹配具体通过文献“吴静,吴志强,胡天跃等.基于构建虚同相轴压制地震层间多次波.地球物理学报.2013,56(3):985–994”记载的多道最小二乘匹配滤波方法，将步骤C3所述预测得到的多次波数据与叠前CMP道集中的多次波做振幅和相位的匹配。

4.如权利要求1所述层间多次波压制方法，其特征是，步骤B通过计算叠前道集与参考道间的相似度，将所述相似度作为权重函数用于加权叠加，构建得到参考道。

5.如权利要求4所述层间多次波压制方法，其特征是，具体通过式1进行加权叠加得到参考道：

其中，T_i(t)为叠前CMP道集中的第i道；t为旅行时；W_i为第i道记录的权重函数；S(t)为加权叠加得到的参考道；

所述权重函数W_i通过式2表示的相似度确定：

其中，ρ_XY为X和Y之间的相关系数，即其中E(.)表示求取数学期望，反映两道记录间的相似度；X为叠前CMP道集中的第i道；Y为等权重叠加得到的道集；ε为设置的阈值，当ρ_XY相似系数小于ε时，将该道的权重置为零。

6.如权利要求1所述层间多次波压制方法，其特征是，步骤C2具体通过式3将R_upper(t)与每道的G_lower(t)进行互相关，得到虚同相轴数据：

其中，t和τ都表示时间。

7.如权利要求1所述层间多次波压制方法，其特征是，步骤C3具体采用式4进行所述褶积，预测得到与地层w有关的层间多次波：

其中，G_IM(t)为预测的层间多次波；t和τ都表示时间。