[发明专利]一种多源混合地震数据表层相关多次波压制方法

权利要求书

1.一种多源混合地震数据表层相关多次波压制方法，其特征在于包括以下步骤：

1)对于给定的时空域地震数据，按炮点、检波点和偏移距信息进行排列；

2)对排列后地震数据每一道进行傅里叶变换，得到频率域的切片数据P_bl；

3)利用最小二乘方法，求取混合算子的逆；

4)取频率域的切片数据P_bl与混合算子的逆相乘，得到预测的多次波M_bl；

5)对频率域预测得到的多次波进行傅里叶逆变换，得到时空域预测的多次波记录；

6)利用最小能量准则，将得到的多次波记录从原始给定的多源混合地震数据中自适应减去；

7)根据深层多次波压制的效果，判断是否仍有多次波存在，若仍有，则返回步骤2)；反之，则输出多次波压制后的多源混合地震记录。

2.如权利要求1所述方法，其特征在于：所述步骤3)中，混合算子的逆(Γ_bl)^-1求取方法为：

设多源混合地震数据混合算子为Γ_bl，采用最小二乘理论对其进行求逆处理，该算子的最小二乘逆为：

式中，上标H代表矩阵的共轭转置，μ为稳定求逆过程的阻尼因子，I为单位矩阵。

3.如权利要求1所述方法，其特征在于：所述步骤4)中，预测的多次波M_bl为：

M_bl＝P_0,bl(Γ_bl)^-1P_bl

其中，P_0,bl为多源混合地震数据一次波记录。

4.如权利要求3所述方法，其特征在于：所述步骤7)中，若仍有多次波存在，则将预测的多次波M_bl中的多源混合地震数据一次波记录P_0,bl返回步骤2)中，作为得到频率域的切片数据P_bl，进行下一次多次波的预测和相减。

5.如权利要求1所述方法，其特征在于：所述步骤6)中，将得到的多次波从原始给定的多源混合地震数据中自适应减去的方法为：

P_0,bl⁽ⁿ⁺¹⁾＝P_bl-AM_bl

其中，P_0,bl⁽ⁿ⁺¹⁾表示第n次迭代后的多次波压制结果，n表示迭代次数，迭代次数n小于3；A为根据最小能量准则估计出的子波。

6.如权利要求5所述方法，其特征在于：所述子波A满足：

E(A)＝||P_bl-AM_bl||²＝min，

其中，E(A)为相减后数据的整体能量。