[发明专利]自适应的鬼波压制方法及系统

说明书

公开了一种自适应的鬼波压制方法及系统。该方法可以包括：确定地震波场数据中每一道鬼波延迟时间的搜索范围，计算搜索范围内每一道鬼波延迟时间对应的去鬼波结果；根据去鬼波结果，计算搜索范围内最大峭度系数对应的鬼波延迟时间为最优鬼波延迟时间；根据最优鬼波延迟时间，计算鬼波滤波算子；根据地震波场数据与鬼波滤波算子，计算无鬼波上行波场。本发明通过峭度系数最大化的方法精确估计最优鬼波延迟时间，进而利用确定性鬼波压制算法，可最大消除粗糙海水面对鬼波压制的影响，实现海洋拖缆数据鬼波压制。

技术领域

本发明涉及海洋地震资料宽频处理领域，更具体地，涉及一种自适应的鬼波压制方法及系统。

背景技术

海洋地震勘探中，检波器除接收到海底一次反射波外，还接收到与自由表面相关的鬼波。鬼波作为一种特殊的噪声，与海底一次反射波形相似，叠加在一次波的尾部，影响地震记录的分辨率，甚至产生虚假同相轴。

在海洋地震勘探中，由于受到外部采集环境的影响，海水表面随着时间和空间而不断起伏变化。在平坦海水表面假设下进行资料处理，会引入误差，降低对地下介质成像的分辨率。

Jovanovich基于经典散射理论，研究了粗糙海水表面对鬼波能量的影响；Lawsand Kragh提出利用低频信号重构粗糙海水表面形态，降低海水起伏对统计反褶积的影响；Orji等基于克西霍夫-亥姆霍兹积分，基于地震建模实现时变海水表面估算；Egorov等实现了粗糙海水表面反射系数估算，降低海水起伏对反射系数的影响。

因此对于海洋拖缆采集，鬼波的存在限制了数据的频带和精度，常规的鬼波压制技术受采集环境等制约因素明显，鬼波延时难以精确估算，数据拓频效果难以保证，不利于技术普及推广，制约海洋海洋宽频地震勘探的发展。因此，有必要开发一种自适应的鬼波压制方法及系统。

公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

发明内容

本发明提出了一种自适应的鬼波压制方法及系统，其能够峭度系数最大化的方法精确估计最优鬼波延迟时间，进而利用确定性鬼波压制算法，可最大消除粗糙海水面对鬼波压制的影响，实现海洋拖缆数据鬼波压制。

根据本发明的一方面，提出了一种自适应的鬼波压制方法。所述方法可以包括：确定地震波场数据中每一道鬼波延迟时间的搜索范围，计算所述搜索范围内每一道鬼波延迟时间对应的去鬼波结果；根据所述去鬼波结果，计算所述搜索范围内最大峭度系数对应的鬼波延迟时间为最优鬼波延迟时间；根据所述最优鬼波延迟时间，计算鬼波滤波算子；根据所述地震波场数据与所述鬼波滤波算子，计算无鬼波上行波场。

优选地，所述去鬼波结果为：

其中，U_n为去鬼波结果，Δt_n为第n道的鬼波延迟时间，其搜索范围为[Δt_min,Δt_max]，ε为白噪因子，G(Δt_n)表示鬼波滤波算子，P_n表示地震波场数据。

优选地，所述最优鬼波延迟时间为：

其中，u(Δt_n)为去鬼波结果的时间域表达式，k[u_n(Δt_n)]表示u(Δt_n)的峭度系数值，E表示期望值。

优选地，所述鬼波滤波算子为：

G(x,z,f)＝1+R·e^-i2πfΔt (3)

其中，G(x,z,f)表示鬼波滤波算子，R为海水表面反射系数，f表示频率，Δt表示鬼波延迟时间。