[发明专利]自适应去鬼波与宽频准零相位反褶积联合处理方法及系统

权利要求书

1.一种自适应去鬼波与宽频准零相位反褶积联合处理方法，其特征在于，包括：

输入原始地震数据；

针对所述原始地震数据进行自适应去鬼波处理，获得鬼波压制后的地震数据；

针对所述鬼波压制后的地震数据进行宽频准零相位反褶积处理，获得宽频处理后的地震数据；

其中，所述自适应去鬼波处理包括：

针对所述原始地震数据进行一维傅里叶变换，获取原始频率域地震数据；

分别计算炮检点相对于一次反射波的延迟时间，进而计算频率域去鬼波算子；

计算所述原始频率域地震数据与所述频率域去鬼波算子的乘积，获得去鬼波地震数据；

将所述去鬼波地震数据进行一维傅里叶反变换至t-x域，得到所述鬼波压制后的地震数据；

其中，通过公式(3)计算频率域去鬼波算子：

其中，A(f)为频率域去鬼波算子，R₀为海面的反射系数，ω为角频率，ω＝2πf，i表示复指数，t_s为炮点相对于一次反射波的延迟时间，t_r为检波点相对于一次反射波的延迟时间；

其中，所述宽频准零相位反褶积处理包括：

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

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

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

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

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

其中，f_n为根据数据采样率得到的折叠频率。

2.根据权利要求1所述的自适应去鬼波与宽频准零相位反褶积联合处理方法，其中，通过公式(1)计算炮点相对于一次反射波的延迟时间：

通过公式(2)计算检波点相对于一次反射波的延迟时间：

其中，x为地震数据道的偏移距，d_w表示海底深度，d_s表示激发点深度，d_r表示检波点深度，v为海水速度。

3.根据权利要求1所述的自适应去鬼波与宽频准零相位反褶积联合处理方法，其中，通过能量集中法计算反射信号的频带宽度包括：

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

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

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

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

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

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

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

6.一种自适应去鬼波与宽频准零相位反褶积联合处理系统，其特征在于，该系统包括：

存储器，存储有计算机可执行指令；

处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：

输入原始地震数据；

针对所述原始地震数据进行自适应去鬼波处理，获得鬼波压制后的地震数据；

针对所述鬼波压制后的地震数据进行宽频准零相位反褶积处理，获得宽频处理后的地震数据；

其中，所述自适应去鬼波处理包括：

针对所述原始地震数据进行一维傅里叶变换，获取原始频率域地震数据；

分别计算炮检点相对于一次反射波的延迟时间，进而计算频率域去鬼波算子；

计算所述原始频率域地震数据与所述频率域去鬼波算子的乘积，获得去鬼波地震数据；

将所述去鬼波地震数据进行一维傅里叶反变换至t-x域，得到所述鬼波压制后的地震数据；

其中，通过公式(3)计算频率域去鬼波算子：

其中，A(f)为频率域去鬼波算子，R₀为海面的反射系数，ω为角频率，ω＝2πf，i表示复指数，t_s为炮点相对于一次反射波的延迟时间，t_r为检波点相对于一次反射波的延迟时间；

其中，所述宽频准零相位反褶积处理包括：

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

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

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

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

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

其中，f_n为根据数据采样率得到的折叠频率。