[发明专利]一种自适应的多阶频散面波压制方法、装置及设备

权利要求书

1.一种自适应的多阶频散面波压制方法，其特征在于，包括：

获取原始地震炮集数据，对所述原始地震炮集数据进行数据处理后根据预设条件对其进行选取得到频散谱分析源数据；

对所述频散谱分析源数据进行非线性信号相似性分析，得到各控制点的高精度频散谱；

对所述高精度频散谱进行分阶拾取，得到所述各控制点的多阶频散曲线；

根据所有控制点在工区内的相对位置和网格设置情况，对所述多阶频散曲线中的每一阶频散曲线进行插值处理，得到每一阶频散曲线在工区内的各个网格位置处的局部频散曲线；

对所述原始地震炮集数据中的每一道地震数据，根据用户设置的参数选取所述每一道地震数据的相邻道地震数据，根据所述局部频散曲线以及各相邻道的偏移距信息对所述相邻道地震数据分别做相移处理，并对相移处理后的各相邻道地震数据做叠加，得到所述每一道地震数据的面波噪声预测结果；

根据所述每一道地震数据的面波噪声预测结果，将其与原始地震炮集数据中对应的地震数据进行匹配相减，得到消除各阶面波噪声后的地震炮集数据。

2.根据权利要求1所述的方法，其特征在于，所述对所述频散谱分析源数据进行非线性信号相似性分析，具体包括：

使用以下公式计算非线性信号的相似度表征：

其中，S_NLSC(ω,V_ph；σ)为规则化后的非线性信号相似度表征；S_π为背景值；S_NL(ω,V_ph)为一个控制点中的两个接收点处的地震数据对某一频率ω和数据频率在ω的相速度V_ph的非线性信号相关；σ为用于调整频散谱分辨率的系数；ω为圆频率，V_ph为相速度；

根据计算出的相似度表征对所述频散谱分析源数据进行相似性分析。

3.根据权利要求1所述的方法，其特征在于，所述对所述多阶频散曲线中的每一阶频散曲线进行插值处理，得到每一阶频散曲线在工区内的各个网格位置处的局部频散曲线，具体包括：

根据工区网格以及控制点的相对位置对所述多阶频散曲线中的每一阶频散曲线进行插值，得到整个工区范围内各个网格位置处每一阶的频散曲线。

4.根据权利要求1所述的方法，其特征在于，所述根据所述局部频散曲线以及各相邻道的偏移距信息对所述相邻道地震数据分别做相移处理，具体包括：

根据所述局部频散曲线以及各相邻道的偏移距信息通过以下公式对所述相邻道地震数据分别做相移处理：

其中，U_S(ω)为重建后的面波数据，U(ω)为原始数据，X_win为预测道面波时窗的最大范围，h为预测到面波时窗，i为虚部，ω为圆频率，v_p为相速度。

5.根据权利要求1所述的方法，其特征在于，所述将其与原始地震炮集数据中对应的地震数据进行匹配相减，得到消除各阶面波噪声后的地震炮集数据，具体包括：

根据以下公式将其与原始地震炮集数据中对应的地震数据进行匹配相减：

其中，x(t)为原始地震炮集数据，N为计算次数，g_i(t)为自适应滤波器，m_i(t)为面波数据，E为误差能量。

6.根据权利要求1所述的方法，其特征在于，所述对所述原始地震炮集数据进行数据处理后根据预设条件对其进行选取得到频散谱分析源数据，具体包括：

对获取的原始地震炮集数据进行预处理，得到第一地震炮集数据；

基于用户给定的参数，对所述第一地震炮集数据进行滤波和归一化处理，得到第二地震炮集数据；

在工区内不同位置确定控制点，根据所述控制点处的面波发育情况在所述控制点处选取符合预设条件的面波数据，并从所述第二地震炮集数据中选取符合所述预设条件的面波数据，将选出的符合所述预设条件的面波数据作为对应各控制点的频散谱分析源数据。