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

说明书

本发明属于地球物理勘探中的噪声压制技术领域，涉及一种基于最优双曲线积分路径叠加的多次波压制方法。本发明对常规的三维表面相关多次波压制方法进行优化改进，利用联络测线多次波贡献道集与其经抛物Raon变换得到的结果内积，通过拾取内积结果的局部极大值确定双曲线顶点坐标，并对其进行统计选取出现频次较高的少数顶点坐标作为路径叠加顶点，将顶点水平位置坐标作为联络测线多次波贡献道集的偏移距选取结果，得到若干个预测多次波子集，按照顶点位置的时间坐标确定不同区间，以区间振幅最大为原则选择筛选多次波子集，得到最终多次波预测结果，最后通过匹配减去实现多次波压制。

技术领域

本发明属于地球物理勘探中的噪声压制技术领域，具体涉及一种海洋地震勘探中的多次波压制方法，特别涉及一种基于最优双曲线积分路径叠加的多次波压制方法。

背景技术

多次波在海洋地震勘探中，因其周期性和运动学形态特征，会降低地震数据的信噪比，常被作为一种需要去除的相干噪声。目前，较为常用的多次波压制技术是2D SRME技术。然而在地层倾斜情况时，真实的传播路径为三维空间函数，2D方法无法刻画联络测线方向的变化，因此多次波预测问题转变为3D，需引入3D SRME。但3D预测受限于观测系统采样稀疏，直接叠加会引入预测假象，引起波形畸变，使得多次波预测产生了严重的空间假频。

2D SRME方法是由Berkhout和Verschuur等人(Verschuur D J,Berkhout AJ,Wapenaar C P A.Adaptive surface-related multiple elimination[J].GEOPHYSICS,1992,57(9):1166-1177.)于1992年提出，该方法通过消除数据中表面反射率的影响来消除所有与表面相关的多次波，将多次波去除设计为一种反演过程，在该过程中，可以估计源和地表的反射率特性，反演残差就是预测的多次波，属于数据驱动，仅适用于简单地形情况。随着地震勘探难度的加大，3D SRME方法由E.J.van Dedem(Dedem E J V.3D surface-related multiple elimination and interpolation[J].SEG Technical ProgramExpanded Abstracts,1999,17(1):1321.)于1999年提出，在2D SRME基础上，引入联络测线方向，适应于更加复杂的地形。然而，3D预测过程受观测系统联络测线方向采样稀疏限制，因此一系列解决办法被相继提出。数据重建得到3DSRME所需要的密集采样数据方法由Anatoly Baumstein等人(Baumstein A.3D SRME:Data reconstruction and applicationto field data[J].SEG Technical Program Expanded Abstracts,1999,23(1).)于1999年提出；之后，基于稀疏反演的多次波压制方法于2002年由Van Dedem和Verschuur(DedemE J V,Verschuur D J.3D Surface-related Multiple Prediction Using SparseInversion:Experience With Field Data[J].Seg Technical Program ExpandedAbstracts,2002,21(1):2094.)提出，对稀疏的联络测线多次波贡献道集利用稀疏反演在模型空间重建，再通过振幅和相位进行校正叠加得到预测多次波；利用curvelet变换分离一次波与多次波由Xiang Wu(Wu X,Hung B.Adaptive primary-multiple separationusing 3D curvelet transform[J].ASEG Extended Abstracts,2015,2015(1):1155.)于2015年提出；基于数据规则化和稀疏反演的三维表面多次波压制方法由方云峰等人(方云峰,聂红梅,张丽梅,et al.基于数据规则化和稀疏反演的三维表面多次波压制方法[J].地球物理学报,2016,v.59(02):269-277.)于2016年提出，通过数据规则化与反规则化技术,使得数据分布满足三维SRME的要求，之后通过稀疏反演代替联络测线多次波贡献道集求和。但是，上述这些方法采用的稀疏反演方法以及数据重构或者curvelet变换的稳定性较弱，对于不同的复杂地形影响较大。