CN112835103B - 自适应去鬼波与宽频准零相位反褶积联合处理方法及系统 - Google Patents

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

Info

Publication number
CN112835103B
CN112835103B CN202011632555.XA CN202011632555A CN112835103B CN 112835103 B CN112835103 B CN 112835103B CN 202011632555 A CN202011632555 A CN 202011632555A CN 112835103 B CN112835103 B CN 112835103B
Authority
CN
China
Prior art keywords
seismic data
frequency
ghost
domain
zero phase
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011632555.XA
Other languages
English (en)
Other versions
CN112835103A (zh
Inventor
徐洪斌
周云和
何跃明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Dongfang Lianchuang Geophysical Technology Co ltd
Original Assignee
Beijing Dongfang Lianchuang Geophysical Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Dongfang Lianchuang Geophysical Technology Co ltd filed Critical Beijing Dongfang Lianchuang Geophysical Technology Co ltd
Priority to CN202011632555.XA priority Critical patent/CN112835103B/zh
Publication of CN112835103A publication Critical patent/CN112835103A/zh
Application granted granted Critical
Publication of CN112835103B publication Critical patent/CN112835103B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/364—Seismic filtering

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

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

Description

自适应去鬼波与宽频准零相位反褶积联合处理方法及系统
技术领域
本发明涉及地球物理勘探地震资料的处理领域，更具体地，涉及一种自适应去鬼波与宽频准零相位反褶积联合处理方法及系统。
背景技术
鬼波对海洋地震勘探资料处理品质带来了严重的不利影响，一直以来，去鬼波是海洋地震勘探资料处理的主要重点之一。但由于鬼波周期复杂多变，它会随震源和检波器沉放深度、反射深度、偏移距、入射角、出射角、地层倾角等多种因素变化，鬼波压制难度大，当鬼波去除不彻底时，提取子波与真实的地震子波相去甚远，不满足反褶积所需要的子波最小相位假设，后续进行预测反褶积处理会破坏有效地震波的波动特征，脉冲反褶积处理会造成处理品质明显降低，反褶积处理很难发挥其应有的作用，难以实现高分辨率处理的目的。
因此，有必要开发一种自适应去鬼波与宽频准零相位反褶积联合处理方法及系统，解决海洋地震勘探的宽频处理问题。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种自适应去鬼波与宽频准零相位反褶积联合处理方法及系统，其能够解决海洋地震勘探资料子波复杂、宽频处理中可靠拓频及准零相位化的技术问题，提高海洋地震资料处理成果的分辨率和品质，为资料解释提供可靠的基础数据，进而提高勘探开发的成功率，有很好的应用前景。
第一方面，本公开实施例提供了一种自适应去鬼波与宽频准零相位反褶积联合处理方法，包括：
输入原始地震数据；
针对所述原始地震数据进行自适应去鬼波处理，获得鬼波压制后的地震数据；
针对所述鬼波压制后的地震数据进行宽频准零相位反褶积处理，获得宽频处理后的地震数据。
优选地，所述自适应去鬼波处理包括：
针对所述原始地震数据进行一维傅里叶变换，获取原始频率域地震数据；
分别计算炮检点相对于一次反射波的延迟时间，进而计算频率域去鬼波算子；
计算所述原始频率域地震数据与所述频率域去鬼波算子的乘积，获得去鬼波地震数据；
将所述去鬼波地震数据进行一维傅里叶反变换至t-x域，得到所述鬼波压制后的地震数据。
优选地，通过公式(1)计算炮点相对于一次反射波的延迟时间：
Figure BDA0002875276210000021
通过公式(2)计算检波点相对于一次反射波的延迟时间：
Figure BDA0002875276210000022
其中，ts为炮点相对于一次反射波的延迟时间，tr为检波点相对于一次反射波的延迟时间，x为地震数据道的偏移距，dw表示海底深度，ds表示激发点深度，dr表示检波点深度，v为海水速度。
优选地，通过公式(3)计算频率域去鬼波算子：
Figure BDA0002875276210000031
其中，A(f)为频率域去鬼波算子，R0为海面的反射系数，ω为角频率，ω＝2πf，i表示复指数。
优选地，所述宽频准零相位反褶积处理包括：
针对所述鬼波压制后的地震数据的炮集内每道数据每个时窗的地震子波，在t-x域炮集内完成每道各时窗内的脉冲反褶积运算，进而获得f-x域的炮集数据；
通过能量集中法计算反射信号的频带宽度；
根据所述频带宽度，确定宽频准零相位滤波算子；
根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频处理后的地震数据。
优选地，通过能量集中法计算反射信号的频带宽度包括：
计算炮记录的平均功率谱，进而在所述平均功率谱上计算极值点；
所述极值点对应的频率为中心频率，确定反褶积期望频率；
以所述反褶积期望频率为中心向两边分别进行功率谱累加，设定能量门限值，计算所述频带宽度。
优选地，通过公式(4)计算每个时窗的有效信号的低频端频率与高频端频率，确定所述频带宽度：
Figure BDA0002875276210000032
其中，σ为高斯形功率谱的频率均方差，f0为中心频率，P(0)为零频率处对应的功率谱，P(f)为频率f处对应的功率谱。
优选地，所述宽频准零相位滤波算子为：
Figure BDA0002875276210000041
其中，fn为根据数据采样率得到的折叠频率。
优选地，根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频处理后的地震数据包括：
计算频率域将所述宽频准零相位滤波算子与所述f-x域的炮集数据的乘积，将计算结果进行一维傅里叶反变换至t-x域，获得所述宽频处理后的地震数据。
第二方面，本公开实施例还提供了一种自适应去鬼波与宽频准零相位反褶积联合处理系统，包括：
输入原始地震数据；
针对所述原始地震数据进行自适应去鬼波处理，获得鬼波压制后的地震数据；
针对所述鬼波压制后的地震数据进行宽频准零相位反褶积处理，获得宽频处理后的地震数据。
优选地，所述自适应去鬼波处理包括：
针对所述原始地震数据进行一维傅里叶变换，获取原始频率域地震数据；
分别计算炮检点相对于一次反射波的延迟时间，进而计算频率域去鬼波算子；
计算所述原始频率域地震数据与所述频率域去鬼波算子的乘积，获得去鬼波地震数据；
将所述去鬼波地震数据进行一维傅里叶反变换至t-x域，得到所述鬼波压制后的地震数据。
优选地，通过公式(1)计算炮点相对于一次反射波的延迟时间：
Figure BDA0002875276210000051
通过公式(2)计算检波点相对于一次反射波的延迟时间：
Figure BDA0002875276210000052
其中，ts为炮点相对于一次反射波的延迟时间，tr为检波点相对于一次反射波的延迟时间，x为地震数据道的偏移距，dw表示海底深度，ds表示激发点深度，dr表示检波点深度，v为海水速度。
优选地，通过公式(3)计算频率域去鬼波算子：
Figure BDA0002875276210000053
其中，A(f)为频率域去鬼波算子，R0为海面的反射系数，ω为角频率，ω＝2πf，i表示复指数。
优选地，所述宽频准零相位反褶积处理包括：
针对所述鬼波压制后的地震数据的炮集内每道数据每个时窗的地震子波，在t-x域炮集内完成每道各时窗内的脉冲反褶积运算，进而获得f-x域的炮集数据；
通过能量集中法计算反射信号的频带宽度；
根据所述频带宽度，确定宽频准零相位滤波算子；
根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频处理后的地震数据。
优选地，通过能量集中法计算反射信号的频带宽度包括：
计算炮记录的平均功率谱，进而在所述平均功率谱上计算极值点；
所述极值点对应的频率为中心频率，确定反褶积期望频率；
以所述反褶积期望频率为中心向两边分别进行功率谱累加，设定能量门限值，计算所述频带宽度。
优选地，通过公式(4)计算每个时窗的有效信号的低频端频率与高频端频率，确定所述频带宽度：
Figure BDA0002875276210000061
其中，σ为高斯形功率谱的频率均方差，f0为中心频率，P(0)为零频率处对应的功率谱，P(f)为频率f处对应的功率谱。
优选地，所述宽频准零相位滤波算子为：
Figure BDA0002875276210000062
其中，fn为根据数据采样率得到的折叠频率。
优选地，根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频处理后的地震数据包括：
计算频率域将所述宽频准零相位滤波算子与所述f-x域的炮集数据的乘积，将计算结果进行一维傅里叶反变换至t-x域，获得所述宽频处理后的地震数据。
其有益效果在于：
针对要解决的海洋子波拓频处理技术难点，依据自适应去鬼波技术具有良好的去鬼波和宽频准零相位反褶积具有显著拓频的特点，发明了基于自适应去鬼波与宽频准零相位反褶积联合的海洋地震资料宽频处理方法。本发明可较好地解决海洋地震勘探资料子波复杂，宽频处理中可靠拓频及准零相位化的技术问题，可提高海洋地震资料处理成果的分辨率和品质，为资料解释提供可靠的基础数据，进而提高勘探开发的成功率，有很好的应用前景。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的一个实施例的自适应去鬼波与宽频准零相位反褶积联合处理方法的步骤的流程图。
图2a和图2b分别示出了根据本发明的一个实施例的原始地震数据的子波及频谱的示意图。
图3a和图3b分别示出了根据本发明的一个实施例的鬼波压制后的地震数据的子波及频谱的示意图。
图4a和图4b分别示出了根据本发明的一个实施例的宽频处理后的地震数据的子波及频谱的示意图。
图5a、图5b、图5c分别示出了根据本发明的一个实施例的处理前、仅去鬼波后以及去鬼波与宽频准零相位反褶积处理后的单炮记录的对比示意图。
图6a、图6b、图6c分别示出了根据本发明的一个实施例的处理前、仅去鬼波后以及去鬼波与宽频准零相位反褶积处理后的PSTM剖面的对比示意图。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
图1示出了根据本发明的一个实施例的自适应去鬼波与宽频准零相位反褶积联合处理方法的步骤的流程图。
本发明提供一种自适应去鬼波与宽频准零相位反褶积联合处理方法，包括：
步骤101，输入原始地震数据；
具体地，输入加载过观测系统信息道头的地震数据，原始地震数据中包含鬼波，其子波与频谱分别如图2a、图2b所示。
步骤102，针对原始地震数据进行自适应去鬼波处理，获得鬼波压制后的地震数据；在一个示例中，自适应去鬼波处理包括：
针对原始地震数据进行一维傅里叶变换，获取原始频率域地震数据；
分别计算炮检点相对于一次反射波的延迟时间，进而计算频率域去鬼波算子；
计算原始频率域地震数据与频率域去鬼波算子的乘积，获得去鬼波地震数据；
将去鬼波地震数据进行一维傅里叶反变换至t-x域，得到鬼波压制后的地震数据。
在一个示例中，通过公式(1)计算炮点相对于一次反射波的延迟时间：
Figure BDA0002875276210000081
通过公式(2)计算检波点相对于一次反射波的延迟时间：
Figure BDA0002875276210000082
其中，ts为炮点相对于一次反射波的延迟时间，tr为检波点相对于一次反射波的延迟时间，x为地震数据道的偏移距，dw表示海底深度，ds表示激发点深度，dr表示检波点深度，v为海水速度。
在一个示例中，通过公式(3)计算频率域去鬼波算子：
Figure BDA0002875276210000091
其中，A(f)为频率域去鬼波算子，R0为海面的反射系数，ω为角频率，ω＝2πf，i表示复指数。
具体地，假设地震记录x(t)由一次波和各种类型的鬼波叠加而成，设s(t)为地震记录中的一次反射波，则地震记录可以表示为：
Figure BDA0002875276210000092
其中，R0为海面的反射系数，近似为-1，ts和tr为炮检点相对于一次反射波的延迟时间，对公式(6)做一维傅里叶变换，得到原始频率域地震数据为：
Figure BDA0002875276210000093
分别通过公式(1)、(2)计算炮检点相对于一次反射波的延迟时间。
根据公式(7)可以得到鬼波为：
Figure BDA0002875276210000094
则公式(7)可以写成：
X(f)＝S(f)G(f) (9)
从公式(8)可见，存在鬼波干扰的原始频率域地震数据可以看作是一次反射波与一个滤波算子为G(f)的乘积结果。因此，根据(9)式设计一个反滤波算子A(f)为公式(3)，即为频率域去鬼波算子，可以去除原始频率域地震数据中的鬼波干扰。
图3a和图3b分别示出了根据本发明的一个实施例的鬼波压制后的地震数据的子波及频谱的示意图。
在频率域计算原始频率域地震数据与频率域去鬼波算子的乘积，获得去鬼波地震数据。将经过消除鬼波陷波影响的去鬼波地震数据进行一维傅里叶反变换至t-x域，得到鬼波压制后的地震数据，其子波及频谱分别如图3a、图3b所示。
步骤103，针对鬼波压制后的地震数据进行宽频准零相位反褶积处理，获得宽频处理后的地震数据。在一个示例中，宽频准零相位反褶积处理包括：
针对鬼波压制后的地震数据的炮集内每道数据每个时窗的地震子波，在t-x域炮集内完成每道各时窗内的脉冲反褶积运算，进而获得f-x域的炮集数据；
通过能量集中法计算反射信号的频带宽度；
根据频带宽度，确定宽频准零相位滤波算子；
根据f-x域的炮集数据与宽频准零相位滤波算子，计算宽频处理后的地震数据。
在一个示例中，通过能量集中法计算反射信号的频带宽度包括：
计算炮记录的平均功率谱，进而在平均功率谱上计算极值点；
极值点对应的频率为中心频率，确定反褶积期望频率；
以反褶积期望频率为中心向两边分别进行功率谱累加，设定能量门限值，计算频带宽度。
在一个示例中，通过公式(4)计算每个时窗的有效信号的低频端频率与高频端频率，确定频带宽度：
Figure BDA0002875276210000101
其中，σ为高斯形功率谱的频率均方差，f0为中心频率，P(0)为零频率处对应的功率谱，P(f)为频率f处对应的功率谱。
在一个示例中，宽频准零相位滤波算子为：
Figure BDA0002875276210000111
其中，fn为根据数据采样率得到的折叠频率。
在一个示例中，根据f-x域的炮集数据与宽频准零相位滤波算子，计算宽频处理后的地震数据包括：
计算频率域将宽频准零相位滤波算子与f-x域的炮集数据的乘积，将计算结果进行一维傅里叶反变换至t-x域，获得宽频处理后的地震数据。
具体地，若将地震子波作为反滤波的输入，期望输出则为d(t)尖脉冲。脉冲反褶积的基本思想在于设计一个反子波a(t)算子，用它把已知的输入地震信号变换成给定的期望输出尖脉冲信号，这就是脉冲反褶积。脉冲反褶积的基本方程是：
Figure BDA0002875276210000112
一般情况下，地震子波为未知的，为在未知子波的情况下求出反滤波因子，必须对地震子波及反射系数序列加上一定的假设条件，包括：
1、假设反射系数序列R(t)是随机的白噪序列，即其自相关为：
Figure BDA0002875276210000113
2、假设地震子波是最小相位的。
根据假设1，地震子波的自相关Rrr可以用地震记录的自相关Rxx来代替。根据假设2，可知地震子波的Z变换B(z)的零点全部在单位圆外，也即反滤因子a(t)的Z变换A(z)＝1/B(z)的分母多项式的零点全在单位圆外，故a(t)是稳定及物理可实现的。这脉冲反褶积的的方程变为：
Figure BDA0002875276210000121
这就是脉冲反褶积的基本方程，其系数矩阵中各元素可直接由地震记录求得。当求取了反褶积因子a(t)后，令其与地震记录x(t)进行褶积运算，即：
s(t)＝a(t)*x(t) (13)
则s(t)即为经过脉冲反褶积之后输出地新的地震记录。如前所述，脉冲反褶积在压缩地震子波，提高分辨率的同时，也增强的有效信号之外的低频端和高频部分的噪声能量，降低有效反射信号的信噪比。
输入加载过观测系统信息道头的地震炮集数据；在输入的炮集记录上，分时窗利用自相关函数计算炮集内每道数据本时窗的地震子波；针对炮集内每道数据每个时窗的地震子波，利用脉冲反褶积的方法，期望输出为尖脉冲，解托布利兹矩阵求得反褶积反子波算子，在t-x域炮集内完成每道各时窗内的反褶积运算，对运算结果做一维傅里叶变换，获得f-x域的炮集数据。
为了克服脉冲反褶积的缺陷，根据地震资料有效信号频带范围，在频率域设计一个衰减函数，使其在有效频带f1及f2之间接近于1，而在f1和f2处迅速衰减6分贝，在有效频带之外迅速衰减到极小值，对有效信号低频端和高频端的噪声信号做压制衰减处理，保证在提高有效信号分辨率的同时，不降低有效信号的信噪比。其效果就相当于在时间域设计一个宽频零相位滤波器m(t)，对经过脉冲反褶积处理后的地震信号做宽频准零相位褶积运算。
应用能量集中法对反射信号的频带宽度进行估算，其基本原理是针对高斯形有效信号的频谱，自适应估算信号的中心频率后设定一个限定范围，然后将信号能量集中到这个范围来估算出频带宽度，即f1和f2。
高斯形功率谱为公式(4)，按照半功率点的定义，幅度衰减3分贝的变化量△f3db和σ之间的关系如下：
△f3db＝2.335σ (14)
对于任何正态分布，期望为μ，标准差是σ，期望μ的每一边的值都占50％，μ-σ和μ+σ之间有68％的值，μ-2σ和μ+2σ之间有95％的值。如果变量有n个样本，则μ和σ由下式定义：
Figure BDA0002875276210000131
Figure BDA0002875276210000132
在估算频带宽度的过程中，首先求出炮记录的平均功率谱为：
Figure BDA0002875276210000133
其中，E(f)为炮记录中每道频率为f的振幅，N为每炮的道数，然后在平均功率谱上计算极值点，极值点对应的频率为中心频率f0，极值点两边的次极值点的平均定义为反褶积期望频率fm。以反褶积期望输出主频率fm为中心向两边分别进行功率谱累加，设定能量门限值为70％，根据公式(4)可以自动算出每个处理时窗内的f1和f2。
图4a和图4b分别示出了根据本发明的一个实施例的宽频处理后的地震数据的子波及频谱的示意图。
根据炮记录上有效波的频带范围，在频率域低频端f1及高频端f2设计频域衰减算子为公式(5)。计算频率域将宽频准零相位滤波算子与f-x域的炮集数据的乘积；将计算结果进行一维傅里叶反变换至t-x域，获得宽频准零相位反褶积处理后的地震数据，其子波及频谱分别如图4a、图4b所示。
实施例1
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
图5a、图5b、图5c分别示出了根据本发明的一个实施例的处理前、仅去鬼波后以及去鬼波与宽频准零相位反褶积处理后的单炮记录的对比示意图。
经本发明方法处理前后的单炮记录对比如图5a、图5b、图5c所示，可见鬼波压制好、波组特征好。
图6a、图6b、图6c分别示出了根据本发明的一个实施例的处理前、仅去鬼波后以及去鬼波与宽频准零相位反褶积处理后的PSTM剖面的对比示意图。
经本发明方法处理前后的PSTM剖面对比如图6a、图6b、图6c所示，可见分辨率高、波组特征好；有利于解释层位追踪与地震反演，质量得到显著提高。
本发明提供了一种自适应去鬼波与宽频准零相位反褶积联合处理系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
输入原始地震数据；
针对原始地震数据进行自适应去鬼波处理，获得鬼波压制后的地震数据；
针对鬼波压制后的地震数据进行宽频准零相位反褶积处理，获得宽频处理后的地震数据。
在一个示例中，自适应去鬼波处理包括：
针对原始地震数据进行一维傅里叶变换，获取原始频率域地震数据；
分别计算炮检点相对于一次反射波的延迟时间，进而计算频率域去鬼波算子；
计算原始频率域地震数据与频率域去鬼波算子的乘积，获得去鬼波地震数据；
将去鬼波地震数据进行一维傅里叶反变换至t-x域，得到鬼波压制后的地震数据。
在一个示例中，通过公式(1)计算炮点相对于一次反射波的延迟时间：
Figure BDA0002875276210000151
通过公式(2)计算检波点相对于一次反射波的延迟时间：
Figure BDA0002875276210000152
其中，ts为炮点相对于一次反射波的延迟时间，tr为检波点相对于一次反射波的延迟时间，x为地震数据道的偏移距，dw表示海底深度，ds表示激发点深度，dr表示检波点深度，v为海水速度。
在一个示例中，通过公式(3)计算频率域去鬼波算子：
Figure BDA0002875276210000153
其中，A(f)为频率域去鬼波算子，R0为海面的反射系数，ω为角频率，ω＝2πf，i表示复指数。
在一个示例中，宽频准零相位反褶积处理包括：
针对鬼波压制后的地震数据的炮集内每道数据每个时窗的地震子波，在t-x域炮集内完成每道各时窗内的脉冲反褶积运算，进而获得f-x域的炮集数据；
通过能量集中法计算反射信号的频带宽度；
根据频带宽度，确定宽频准零相位滤波算子；
根据f-x域的炮集数据与宽频准零相位滤波算子，计算宽频处理后的地震数据。
在一个示例中，通过能量集中法计算反射信号的频带宽度包括：
计算炮记录的平均功率谱，进而在平均功率谱上计算极值点；
极值点对应的频率为中心频率，确定反褶积期望频率；
以反褶积期望频率为中心向两边分别进行功率谱累加，设定能量门限值，计算频带宽度。
在一个示例中，通过公式(4)计算每个时窗的有效信号的低频端频率与高频端频率，确定频带宽度：
Figure BDA0002875276210000161
其中，σ为高斯形功率谱的频率均方差，f0为中心频率，P(0)为零频率处对应的功率谱，P(f)为频率f处对应的功率谱。
在一个示例中，宽频准零相位滤波算子为：
Figure BDA0002875276210000162
其中，fn为根据数据采样率得到的折叠频率。
在一个示例中，根据f-x域的炮集数据与宽频准零相位滤波算子，计算宽频处理后的地震数据包括：
计算频率域将宽频准零相位滤波算子与f-x域的炮集数据的乘积，将计算结果进行一维傅里叶反变换至t-x域，获得宽频处理后的地震数据。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。

Claims (6)

1.一种自适应去鬼波与宽频准零相位反褶积联合处理方法，其特征在于，包括：
输入原始地震数据；
针对所述原始地震数据进行自适应去鬼波处理，获得鬼波压制后的地震数据；
针对所述鬼波压制后的地震数据进行宽频准零相位反褶积处理，获得宽频处理后的地震数据；
其中，所述自适应去鬼波处理包括：
针对所述原始地震数据进行一维傅里叶变换，获取原始频率域地震数据；
分别计算炮检点相对于一次反射波的延迟时间，进而计算频率域去鬼波算子；
计算所述原始频率域地震数据与所述频率域去鬼波算子的乘积，获得去鬼波地震数据；
将所述去鬼波地震数据进行一维傅里叶反变换至t-x域，得到所述鬼波压制后的地震数据；
其中，通过公式(3)计算频率域去鬼波算子：
Figure FDA0003177621000000011
其中，A(f)为频率域去鬼波算子，R0为海面的反射系数，ω为角频率，ω＝2πf，i表示复指数，ts为炮点相对于一次反射波的延迟时间，tr为检波点相对于一次反射波的延迟时间；
其中，所述宽频准零相位反褶积处理包括：
针对所述鬼波压制后的地震数据的炮集内每道数据每个时窗的地震子波，在t-x域炮集内完成每道各时窗内的脉冲反褶积运算，进而获得f-x域的炮集数据；
通过能量集中法计算反射信号的频带宽度；
根据所述频带宽度，确定宽频准零相位滤波算子；
根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频处理后的地震数据；
其中，所述宽频准零相位滤波算子为：
Figure FDA0003177621000000021
其中，fn为根据数据采样率得到的折叠频率。
2.根据权利要求1所述的自适应去鬼波与宽频准零相位反褶积联合处理方法，其中，通过公式(1)计算炮点相对于一次反射波的延迟时间：
Figure FDA0003177621000000022
通过公式(2)计算检波点相对于一次反射波的延迟时间：
Figure FDA0003177621000000023
其中，x为地震数据道的偏移距，dw表示海底深度，ds表示激发点深度，dr表示检波点深度，v为海水速度。
3.根据权利要求1所述的自适应去鬼波与宽频准零相位反褶积联合处理方法，其中，通过能量集中法计算反射信号的频带宽度包括：
计算炮记录的平均功率谱，进而在所述平均功率谱上计算极值点；
所述极值点对应的频率为中心频率，确定反褶积期望频率；
以所述反褶积期望频率为中心向两边分别进行功率谱累加，设定能量门限值，计算所述频带宽度。
4.根据权利要求1所述的自适应去鬼波与宽频准零相位反褶积联合处理方法，其中，通过公式(4)计算每个时窗的有效信号的低频端频率与高频端频率，确定所述频带宽度：
Figure FDA0003177621000000031
其中，σ为高斯形功率谱的频率均方差，f0为中心频率，P(0)为零频率处对应的功率谱，P(f)为频率f处对应的功率谱。
5.根据权利要求1所述的自适应去鬼波与宽频准零相位反褶积联合处理方法，其中，根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频处理后的地震数据包括：
计算频率域将所述宽频准零相位滤波算子与所述f-x域的炮集数据的乘积，将计算结果进行一维傅里叶反变换至t-x域，获得所述宽频处理后的地震数据。
6.一种自适应去鬼波与宽频准零相位反褶积联合处理系统，其特征在于，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
输入原始地震数据；
针对所述原始地震数据进行自适应去鬼波处理，获得鬼波压制后的地震数据；
针对所述鬼波压制后的地震数据进行宽频准零相位反褶积处理，获得宽频处理后的地震数据；
其中，所述自适应去鬼波处理包括：
针对所述原始地震数据进行一维傅里叶变换，获取原始频率域地震数据；
分别计算炮检点相对于一次反射波的延迟时间，进而计算频率域去鬼波算子；
计算所述原始频率域地震数据与所述频率域去鬼波算子的乘积，获得去鬼波地震数据；
将所述去鬼波地震数据进行一维傅里叶反变换至t-x域，得到所述鬼波压制后的地震数据；
其中，通过公式(3)计算频率域去鬼波算子：
Figure FDA0003177621000000041
其中，A(f)为频率域去鬼波算子，R0为海面的反射系数，ω为角频率，ω＝2πf，i表示复指数，ts为炮点相对于一次反射波的延迟时间，tr为检波点相对于一次反射波的延迟时间；
其中，所述宽频准零相位反褶积处理包括：
针对所述鬼波压制后的地震数据的炮集内每道数据每个时窗的地震子波，在t-x域炮集内完成每道各时窗内的脉冲反褶积运算，进而获得f-x域的炮集数据；
通过能量集中法计算反射信号的频带宽度；
根据所述频带宽度，确定宽频准零相位滤波算子；
根据所述f-x域的炮集数据与所述宽频准零相位滤波算子，计算宽频处理后的地震数据；
其中，所述宽频准零相位滤波算子为：
Figure FDA0003177621000000042
其中，fn为根据数据采样率得到的折叠频率。
CN202011632555.XA 2020-12-31 2020-12-31 自适应去鬼波与宽频准零相位反褶积联合处理方法及系统 Active CN112835103B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011632555.XA CN112835103B (zh) 2020-12-31 2020-12-31 自适应去鬼波与宽频准零相位反褶积联合处理方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011632555.XA CN112835103B (zh) 2020-12-31 2020-12-31 自适应去鬼波与宽频准零相位反褶积联合处理方法及系统

Publications (2)

Publication Number Publication Date
CN112835103A CN112835103A (zh) 2021-05-25
CN112835103B true CN112835103B (zh) 2021-10-08

Family

ID=75926103

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011632555.XA Active CN112835103B (zh) 2020-12-31 2020-12-31 自适应去鬼波与宽频准零相位反褶积联合处理方法及系统

Country Status (1)

Country Link
CN (1) CN112835103B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113341462A (zh) * 2021-06-10 2021-09-03 广州海洋地质调查局 一种海洋地震数据零相位处理方法
CN113514889B (zh) * 2021-07-13 2022-06-21 中山大学 一种提升海洋深反射地震数据中低频信号能量的处理方法
CN114755723A (zh) * 2022-03-29 2022-07-15 北京东方联创地球物理技术有限公司 海上地震数据确定极性零相位化方法、装置、设备及介质
CN114859409B (zh) * 2022-04-11 2023-03-31 中山大学 一种大洋岩石圈间断面信息的获取方法及装置
CN116559940A (zh) * 2023-03-27 2023-08-08 广州海洋地质调查局 一种潜山精细成像的地震数据处理方法及装置
CN116774280B (zh) * 2023-06-25 2024-05-24 中海石油(中国)有限公司深圳分公司 一种鬼波压制定量质控方法、装置、电子设备及存储介质

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10345469B2 (en) * 2016-01-29 2019-07-09 Cgg Services Sas Device and method for correcting seismic data for variable air-water interface
CN106526677B (zh) * 2016-10-26 2018-12-21 中海石油(中国)有限公司 一种海上自适应压制鬼波的宽频逆时偏移成像方法
CN106896409B (zh) * 2017-03-14 2018-12-07 中国海洋石油集团有限公司 一种基于波动方程边值反演的变深度缆鬼波压制方法
CN107179550B (zh) * 2017-07-05 2018-12-07 西安交通大学 一种数据驱动的地震信号零相位反褶积方法
CN108919357B (zh) * 2018-05-16 2019-10-11 中国海洋石油集团有限公司 一种基于频谱重构的鬼波压制方法
CN110749923A (zh) * 2018-07-24 2020-02-04 中国石油化工股份有限公司 一种基于范数方程提高分辨率的反褶积方法
CN110967735A (zh) * 2018-09-28 2020-04-07 中国石油化工股份有限公司 自适应的鬼波压制方法及系统

Also Published As

Publication number Publication date
CN112835103A (zh) 2021-05-25

Similar Documents

Publication Publication Date Title
CN112835103B (zh) 自适应去鬼波与宽频准零相位反褶积联合处理方法及系统
CN112817047B (zh) 海洋地震自适应去鬼波方法、装置、电子设备及介质
Xu et al. Monochromatic noise removal via sparsity-enabled signal decomposition method
CN110967735A (zh) 自适应的鬼波压制方法及系统
CN112817040B (zh) 宽频准零相位反褶积处理方法、装置、电子设备及介质
CN110646841B (zh) 时变稀疏反褶积方法及系统
Tian et al. Statistical analysis of split spectrum processing for multiple target detection
CN111060879A (zh) 一种基于二维匹配滤波结果的联合副瓣抑制方法
CN112213773B (zh) 一种地震分辨率提高方法及电子设备
Santoso et al. Performance of various speckle reduction filters on Synthetic Aperture Radar image
CN110703332B (zh) 一种鬼波压制方法
CN113156514B (zh) 基于主频波数域均值滤波的地震数据去噪方法及系统
CN115712146A (zh) 基于频率慢度域延拓的鬼波参数最优化拖缆鬼波压制方法
CN116047504A (zh) 一种改进反褶积抑制探地雷达多次波的方法
Weiss et al. Wavelet-based denoising of underwater acoustic signals
CN110837119B (zh) 一种增强反q值补偿稳定性的地震资料处理方法及系统
CN113009464A (zh) 基于线性约束最小方差准则的稳健自适应脉冲压缩方法
CN109459788B (zh) 地层品质因子计算方法及系统
CN117269928B (zh) 基于动目标检测雷达的多普勒过采样投影杂波抑制方法
Li et al. Ultrasound image enhancement using dynamic filtering
CN114740530B (zh) 基于双曲时窗约束的中高频拟线性噪声压制方法及装置
CN113358927B (zh) 一种基于区域核函数的多分量线性调频信号时频分析方法
CN111694057B (zh) 一种压制地震资料涌浪噪声的方法、存储介质及设备
Yang et al. Sidelobe suppression of SAR images by spectrum shaping
CN112526604B (zh) 基于目的层谱分析的自适应低频补偿方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant