[发明专利]宽频准零相位反褶积处理方法、装置、电子设备及介质

权利要求书

1.一种宽频准零相位反褶积处理方法，其特征在于，包括：

针对炮集内每道数据每个时窗的地震子波，在t-x域炮集内完成每道各时窗内的反褶积运算，进而获得f-x域的炮集数据；

通过能量集中法计算反射信号的频带宽度；

根据所述频带宽度，确定宽频准零相位滤波算子；

根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频准零相位反褶积处理后的地震数据；

其中，通过能量集中法计算反射信号的频带宽度包括：

计算炮记录的平均功率谱，进而在所述平均功率谱上计算极值点；

所述极值点对应的频率为中心频率，确定反褶积期望频率；

以所述反褶积期望频率为中心向两边分别进行功率谱累加，设定能量门限值，计算所述频带宽度；

其中，所述宽频准零相位滤波算子为：

其中，f_n为根据数据采样率得到的折叠频率，f₁为有效频带的低频端点，f₂为有效频带的高频端点。

2.根据权利要求1所述的宽频准零相位反褶积处理方法，其中，所述反褶积期望频率为所述极值点两边的次极值点的平均值。

3.根据权利要求1所述的宽频准零相位反褶积处理方法，其中，通过公式(1)计算每个时窗的有效信号的低频端频率与高频端频率，确定所述频带宽度：

其中，σ为高斯形功率谱的频率均方差，f₀为中心频率，P(0)为零频率处对应的功率谱，P(f)为频率f处对应的功率谱。

4.根据权利要求1所述的宽频准零相位反褶积处理方法，其中，根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频准零相位反褶积处理后的地震数据包括：

在频率域将所述宽频准零相位滤波算子与所述f-x域的炮集数据进行宽频准零相位反褶积计算；

将计算结果进行一维傅里叶反变换至t-x域，获得所述宽频准零相位反褶积处理后的地震数据。

5.根据权利要求4所述的宽频准零相位反褶积处理方法，其中，通过公式(3)进行宽频准零相位反褶积计算：

S(t)＝a(t)*x(t)*m(t,f) (3)

其中，a(t)为反褶积因子，x(t)为原始地震记录。

6.一种宽频准零相位反褶积处理装置，其特征在于，包括：

反褶积模块，针对炮集内每道数据每个时窗的地震子波，在t-x域炮集内完成每道各时窗内的反褶积运算，进而获得f-x域的炮集数据；

频带宽度计算模块，通过能量集中法计算反射信号的频带宽度；

滤波算子设置模块，根据所述频带宽度，确定宽频准零相位滤波算子；

宽频准零相位滤波模块，根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频准零相位反褶积处理后的地震数据；

其中，通过能量集中法计算反射信号的频带宽度包括：

计算炮记录的平均功率谱，进而在所述平均功率谱上计算极值点；

所述极值点对应的频率为中心频率，确定反褶积期望频率；

以所述反褶积期望频率为中心向两边分别进行功率谱累加，设定能量门限值，计算所述频带宽度；

其中，所述宽频准零相位滤波算子为：

其中，f_n为根据数据采样率得到的折叠频率，f₁为有效频带的低频端点，f₂为有效频带的高频端点。

7.一种电子设备，其特征在于，所述电子设备包括：

存储器，存储有可执行指令；

处理器，所述处理器运行所述存储器中的所述可执行指令，以实现权利要求1-5中任一项所述的宽频准零相位反褶积处理方法。

8.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1-5中任一项所述的宽频准零相位反褶积处理方法。