自适应去鬼波与宽频准零相位反褶积联合处理方法及系统技术方案

本申请公开了一种自适应去鬼波与宽频准零相位反褶积联合处理方法及系统。该方法可以包括：输入原始地震数据；针对原始地震数据进行自适应去鬼波处理，获得鬼波压制后的地震数据；针对鬼波压制后的地震数据进行宽频准零相位反褶积处理，获得宽频处理后的地震数据。本发明专利技术解决海洋地震勘探资料子波复杂、宽频处理中可靠拓频及准零相位化的技术问题，提高海洋地震资料处理成果的分辨率和品质，为资料解释提供可靠的基础数据，进而提高勘探开发的成功率，有很好的应用前景。

【技术实现步骤摘要】
自适应去鬼波与宽频准零相位反褶积联合处理方法及系统
本专利技术涉及地球物理勘探地震资料的处理领域，更具体地，涉及一种自适应去鬼波与宽频准零相位反褶积联合处理方法及系统。
技术介绍
鬼波对海洋地震勘探资料处理品质带来了严重的不利影响，一直以来，去鬼波是海洋地震勘探资料处理的主要重点之一。但由于鬼波周期复杂多变，它会随震源和检波器沉放深度、反射深度、偏移距、入射角、出射角、地层倾角等多种因素变化，鬼波压制难度大，当鬼波去除不彻底时，提取子波与真实的地震子波相去甚远，不满足反褶积所需要的子波最小相位假设，后续进行预测反褶积处理会破坏有效地震波的波动特征，脉冲反褶积处理会造成处理品质明显降低，反褶积处理很难发挥其应有的作用，难以实现高分辨率处理的目的。因此，有必要开发一种自适应去鬼波与宽频准零相位反褶积联合处理方法及系统，解决海洋地震勘探的宽频处理问题。公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
本专利技术提出了一种自适应去鬼波与宽频准零相位反褶积联合处理方法及系统，其能够解决海洋地震勘探资料子波复杂、宽频处理中可靠拓频及准零相位化的技术问题，提高海洋地震资料处理成果的分辨率和品质，为资料解释提供可靠的基础数据，进而提高勘探开发的成功率，有很好的应用前景。第一方面，本公开实施例提供了一种自适应去鬼波与宽频准零相位反褶积联合处理方法，包括：输入原始地震数据；针对所述原始地震数据进行自适应去鬼波处理，获得鬼波压制后的地震数据；针对所述鬼波压制后的地震数据进行宽频准零相位反褶积处理，获得宽频处理后的地震数据。优选地，所述自适应去鬼波处理包括：针对所述原始地震数据进行一维傅里叶变换，获取原始频率域地震数据；分别计算炮检点相对于一次反射波的延迟时间，进而计算频率域去鬼波算子；计算所述原始频率域地震数据与所述频率域去鬼波算子的乘积，获得去鬼波地震数据；将所述去鬼波地震数据进行一维傅里叶反变换至t-x域，得到所述鬼波压制后的地震数据。优选地，通过公式(1)计算炮点相对于一次反射波的延迟时间：通过公式(2)计算检波点相对于一次反射波的延迟时间：其中，ts为炮点相对于一次反射波的延迟时间，tr为检波点相对于一次反射波的延迟时间，x为地震数据道的偏移距，dw表示海底深度，ds表示激发点深度，dr表示检波点深度，v为海水速度。优选地，通过公式(3)计算频率域去鬼波算子：其中，A(f)为频率域去鬼波算子，R0为海面的反射系数，ω为角频率，ω＝2πf，i表示复指数。优选地，所述宽频准零相位反褶积处理包括：针对所述鬼波压制后的地震数据的炮集内每道数据每个时窗的地震子波，在t-x域炮集内完成每道各时窗内的脉冲反褶积运算，进而获得f-x域的炮集数据；通过能量集中法计算反射信号的频带宽度；根据所述频带宽度，确定宽频准零相位滤波算子；根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频处理后的地震数据。优选地，通过能量集中法计算反射信号的频带宽度包括：计算炮记录的平均功率谱，进而在所述平均功率谱上计算极值点；所述极值点对应的频率为中心频率，确定反褶积期望频率；以所述反褶积期望频率为中心向两边分别进行功率谱累加，设定能量门限值，计算所述频带宽度。优选地，通过公式(4)计算每个时窗的有效信号的低频端频率与高频端频率，确定所述频带宽度：其中，σ为高斯形功率谱的频率均方差，f0为中心频率，P(0)为零频率处对应的功率谱，P(f)为频率f处对应的功率谱。优选地，所述宽频准零相位滤波算子为：其中，fn为根据数据采样率得到的折叠频率。优选地，根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频处理后的地震数据包括：计算频率域将所述宽频准零相位滤波算子与所述f-x域的炮集数据的乘积，将计算结果进行一维傅里叶反变换至t-x域，获得所述宽频处理后的地震数据。第二方面，本公开实施例还提供了一种自适应去鬼波与宽频准零相位反褶积联合处理系统，包括：输入原始地震数据；针对所述原始地震数据进行自适应去鬼波处理，获得鬼波压制后的地震数据；针对所述鬼波压制后的地震数据进行宽频准零相位反褶积处理，获得宽频处理后的地震数据。优选地，所述自适应去鬼波处理包括：针对所述原始地震数据进行一维傅里叶变换，获取原始频率域地震数据；分别计算炮检点相对于一次反射波的延迟时间，进而计算频率域去鬼波算子；计算所述原始频率域地震数据与所述频率域去鬼波算子的乘积，获得去鬼波地震数据；将所述去鬼波地震数据进行一维傅里叶反变换至t-x域，得到所述鬼波压制后的地震数据。优选地，通过公式(1)计算炮点相对于一次反射波的延迟时间：通过公式(2)计算检波点相对于一次反射波的延迟时间：其中，ts为炮点相对于一次反射波的延迟时间，tr为检波点相对于一次反射波的延迟时间，x为地震数据道的偏移距，dw表示海底深度，ds表示激发点深度，dr表示检波点深度，v为海水速度。优选地，通过公式(3)计算频率域去鬼波算子：其中，A(f)为频率域去鬼波算子，R0为海面的反射系数，ω为角频率，ω＝2πf，i表示复指数。优选地，所述宽频准零相位反褶积处理包括：针对所述鬼波压制后的地震数据的炮集内每道数据每个时窗的地震子波，在t-x域炮集内完成每道各时窗内的脉冲反褶积运算，进而获得f-x域的炮集数据；通过能量集中法计算反射信号的频带宽度；根据所述频带宽度，确定宽频准零相位滤波算子；根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频处理后的地震数据。优选地，通过能量集中法计算反射信号的频带宽度包括：计算炮记录的平均功率谱，进而在所述平均功率谱上计算极值点；所述极值点对应的频率为中心频率，确定反褶积期望频率；以所述反褶积期望频率为中心向两边分别进行功率谱累加，设定能量门限值，计算所述频带宽度。优选地，通过公式(4)计算每个时窗的有效信号的低频端频率与高频端频率，确定所述频带宽度：其中，σ为高斯形功率谱的频率均方差，f0为中心频率，P(0)为零频率处对应的功率谱，P(f)为频率f处对应的功率谱。优选地，所述宽频准零相位滤波算子为：其中，fn为根据数据采样率得到的折叠频率。优选地，根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频处理后的地震数据包括：计算频率域将所述宽频准零相位滤波算子与所述f-x域的炮集数据的乘积，将计算结果进行一维傅里叶反变换至t-x域，获得所述宽频本文档来自技高网...

【技术保护点】
1.一种自适应去鬼波与宽频准零相位反褶积联合处理方法，其特征在于，包括：/n输入原始地震数据；/n针对所述原始地震数据进行自适应去鬼波处理，获得鬼波压制后的地震数据；/n针对所述鬼波压制后的地震数据进行宽频准零相位反褶积处理，获得宽频处理后的地震数据。/n

【技术特征摘要】
1.一种自适应去鬼波与宽频准零相位反褶积联合处理方法，其特征在于，包括：
输入原始地震数据；
针对所述原始地震数据进行自适应去鬼波处理，获得鬼波压制后的地震数据；
针对所述鬼波压制后的地震数据进行宽频准零相位反褶积处理，获得宽频处理后的地震数据。


2.根据权利要求1所述的自适应去鬼波与宽频准零相位反褶积联合处理方法，其中，所述自适应去鬼波处理包括：
针对所述原始地震数据进行一维傅里叶变换，获取原始频率域地震数据；
分别计算炮检点相对于一次反射波的延迟时间，进而计算频率域去鬼波算子；
计算所述原始频率域地震数据与所述频率域去鬼波算子的乘积，获得去鬼波地震数据；
将所述去鬼波地震数据进行一维傅里叶反变换至t-x域，得到所述鬼波压制后的地震数据。


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



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



其中，ts为炮点相对于一次反射波的延迟时间，tr为检波点相对于一次反射波的延迟时间，x为地震数据道的偏移距，dw表示海底深度，ds表示激发点深度，dr表示检波点深度，v为海水速度。


4.根据权利要求2所述的自适应去鬼波与宽频准零相位反褶积联合处理方法，其中，通过公式(3)计算频率域去鬼波算子：



其中，A(f)为频率域去鬼波算子，R0为海面的反射系数，ω为角频率，ω＝2πf，i表示复指数。


5.根据权利要求1所述的自适应去鬼波与宽频准零相位反褶积联合处理方法，其中，所述宽频准零相位反褶积处理包括：
针对所述鬼波压制后的地震数据的炮集内每道数据每个时窗的地震子波，在t-x域炮集内完成每道各时窗内的脉冲反褶积运算，进而获得f-x域的炮集数据；
通过能量集中法计算反射信号的频带宽度；
根据所...

【专利技术属性】
技术研发人员：徐洪斌，周云和，何跃明，
申请(专利权)人：北京东方联创地球物理技术有限公司，
类型：发明
国别省市：北京;11