[发明专利]基于最优双曲线积分路径叠加的多次波压制方法

权利要求书

1.一种基于最优双曲线积分路径叠加的多次波压制方法，其特征在于，包括以下步骤：

a、将三维海洋地震勘探采集的数据，通过互易定理补全数据，得到炮道相同的地震数据，并按照道头信息中的炮号抽道集，得到共炮点道集d(x_k,y_k,t,x_s,y_s)；按照检波号抽道集，得到共检波点道集r(x_r,y_r,t,x_k,y_k)，其中，(x_k,y_k)为任一点位置炮点或检波点，(x_s,y_s)为炮点坐标，(x_r,y_r)为检波点坐标，t为时间；

b、对共炮点道集d(x_k,y_k,t,x_s,y_s)与共检波点道集r(x_r,y_r,t,x_k,y_k)进行傅里叶正变换得到频率域的共炮点道集D(x_k,y_k,w,x_s,y_s)与共检波点道集R(x_r,y_r,w,x_k,y_k)，其中，w为角频率；

c、将共炮点道集D(x_k,y_k,w,x_s,y_s)与共检波点道集R(x_r,y_r,w,x_k,y_k)在频率域乘积得到对应不同位置(x_k,y_k)的纵测线多次波贡献道集M_xy(x_r,y_r,x_s,y_s,x_k,y_k,w)；

M_xy(x_r,y_r,x_s,y_s,x_k,y_k,w)＝r₀R(x_r,y_r,w；x_k,y_k)D(x_k,y_k,w；x_s,y_s)

其中，r₀表示地层反射系数，为-1；

d、对不同位置的纵测线多次波贡献道集M_xy(x_r,y_r,x_s,y_s,x_k,y_k,w)沿纵测线方向求和得到联络测线多次波贡献道集M_y(x_r,y_r,x_s,y_s,y_k,w)，并通过傅里叶逆变换变换回时间域，得到时间域的联络测线多次波贡献道集m_y(x_r,y_r,x_s,y_s,y_k,t)，

e、对联络测线多次波贡献道集作抛物Radon变换，将抛物Radon变换的结果与未作抛物Radon变换的联络测线多次波贡献道集内积，

m_yradon(x_r,y_r,x_s,y_s,y_k,t)＝Radon(m_y(x_r,y_r,x_s,y_s,y_k,t))

m_yradon(x_r,y_r,x_s,y_s,y_k,t)为Radon变换后的结果，m_d(x_r,y_r,x_s,y_s,y_k,t)为内积的结果；

f、通过阈值限制对比振幅值拾取局部极大值，获取双曲线顶点坐标；

[x_m,t_m]＝findpeaks(m_d(x_r,y_r,x_s,y_s,y_k,t))

[x_m,t_m]为同相轴顶点的横向坐标组成的矩阵；

g、对所有同相轴顶点坐标进行统计，选取出现频率最高的多个作为路径叠加的顶点，

为路径叠加的顶点，i＝1，2，3...；

h、分别选取路径叠加的多个顶点的x_max(i)作为Radon变换的偏移距，进行Radon变换，将Radon变换的结果叠加得到若干个预测的多次波m_(i)(x_r,y_r,x_s,y_s,t)，对比t_max(i)的大小找到时间方向的不同区间

i、在不同区间中，选取上述若干个预测的多次波中振幅较大的部分组成最终的多次波预测结果；

j、采用时间域匹配减去预测的多次波，得到不含多次波的地震数据。

2.根据权利要求1所述的一种基于最优双曲线积分路径叠加的多次波压制方法，其特征在于：步骤g，所述顶点为3—5个。