CN117148443B - 基于鬼波提取与转换的浅剖数据信噪比增强方法 - Google Patents
基于鬼波提取与转换的浅剖数据信噪比增强方法 Download PDFInfo
- Publication number
- CN117148443B CN117148443B CN202311403468.0A CN202311403468A CN117148443B CN 117148443 B CN117148443 B CN 117148443B CN 202311403468 A CN202311403468 A CN 202311403468A CN 117148443 B CN117148443 B CN 117148443B
- Authority
- CN
- China
- Prior art keywords
- ghost
- profile
- wave
- shallow
- profile data
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 19
- 238000000605 extraction Methods 0.000 title claims abstract description 19
- 230000002238 attenuated effect Effects 0.000 claims abstract description 10
- 230000001629 suppression Effects 0.000 claims abstract description 6
- 238000001228 spectrum Methods 0.000 claims description 17
- 238000001914 filtration Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000010183 spectrum analysis Methods 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 238000003786 synthesis reaction Methods 0.000 claims description 5
- 239000013535 sea water Substances 0.000 claims description 4
- 239000013598 vector Substances 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000000737 periodic effect Effects 0.000 claims description 3
- 230000002708 enhancing effect Effects 0.000 claims 5
- 230000009466 transformation Effects 0.000 claims 2
- 230000002194 synthesizing effect Effects 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000013049 sediment Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- 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
技术领域
本发明涉及地震资料处理与分析领域，是一种基于鬼波提取与转换的浅剖数据信噪比增强方法。
背景技术
浅层剖面勘测方法是一种基于声学原理的连续走航式探测水下浅地层结构和构造的地球物理方法,是进行海洋地球物理调查的常用手段之一,其利用声波在海水和海底沉积物中的传播和反射特性及规律对海底沉积物结构和构造进行连续探测，从而获得直观的海底浅部地层结构剖面。浅地层剖面方法测量结果连续性好，能快速地探测水下地层的地质特征及其分布，并且纵向分辨率较高,因此其在海洋调查中得到了广泛的应用。
采用某些类型的设备进行浅地层剖面测量时，若无法有效遮挡由地下反射回来的地震波，其到达海面后发生下行反射，被采集装置接收后形成具有一定周期、振幅较强的鬼波干扰。一般来说，鬼波与同时段的有效信号相交或重合，打断了海底以下地质层位的分布趋势，甚至会掩盖有效的层位信息。此外，在缺乏有效认识的情况下，鬼波会形成虚假构造，将严重影响资料的真实性与可靠性，进而误导后续的地质解释。
发明内容
本发明要解决的技术问题是提供一种基于鬼波提取与转换的浅剖数据信噪比增强方法，其目的是根据鬼波与一次波的周期性对其进行预测与压制，并利用原始与压制剖面的差提取出前者中的鬼波信号。在已知鬼波周期的情况下，基于鬼波预测的反过程将提取的鬼波剖面转换为一次波剖面。随后，将衰减鬼波后的剖面与鬼波转换的一次波剖面相加，从而得到信噪比显著增强的浅剖数据。
本发明是通过如下技术方案来实现的:
基于鬼波提取与转换的浅剖数据信噪比增强方法，所述方法包括(1)基于自相关谱分析的鬼波周期或延迟时提取，(2)基于时间延迟的鬼波预测，(3)基于最小平方滤波的鬼波剔除，(4)鬼波的提取与转换，(5)高信噪比剖面的合成。
进一步，所述的基于自相关谱分析的鬼波周期或延迟时提取为对于浅剖面数据中具有周期性特征的鬼波干扰，使用自相关谱分析获取其准确延迟时。具体为分析鬼波与一次波二者间的周期性，令最大的相关延迟时为m，要求m的值大于鬼波周期，创建浅剖数据s(x,t)的自相关谱，可得
其中，s(x,t)表示浅剖数据，x表示剖面中各道数据的坐标，t代表旅行时；r(x,τ)表示浅剖数据的自相关谱，τ为相关延迟时；
在基于式(1)所创建的自相关谱r(x,τ)中，通过在谱中拾取极小值的方式获得鬼波周期或延迟时，其表示为一条随坐标x变化的旅行时曲线Δτ(x)，即
Δτ(x)＝min[r(x,τ)] (2)
其中，Δτ表示鬼波的周期或延迟时；min表示通过人工拾取方式获得自相关谱中鬼波周期的过程。
所述的基于时间延迟的鬼波预测；设海面的反射系数为-1，在已知鬼波周期或延迟时的情况下，利用原始数据来预测鬼波数据，时间域的鬼波预测方程为
其中，A(x)表示鬼波传播的几何扩散效应，v为海水速度，通常取1500m/s；g(x,t)为预测的鬼波剖面。
进一步，所述的基于最小平方滤波的鬼波剔除，采用最小平方滤波技术进行鬼波压制；对于基于L2范数的自适应相减处理，通过使下式中的误差能量平方和e值为最小来确定滤波因子
e＝||s-g*a||2 (4)
其中，a表示待求取的滤波因子；向量s与g分别代表浅剖数据s(x,t)以及预测的鬼波剖面g(x,t)；*表示褶积运算，||…||2基于L2范数的最小平方约束过程；
通过求解式(4)可获得滤波因子a，当滤波因子a的长度为1时，直接给出滤波因子的计算公式
式中，上标“T”表示矩阵的转置运算；
将滤波因子a代入下式实现鬼波的压制，相应的计算公式为
p(x,t)＝s(x,t)-a*g(x,t) (6)
式中，p(x,t)表示剔除鬼波的浅剖数据。
进一步，所述的鬼波的提取与转换，通过对衰减鬼波前、后的数据作差获得原始数据中的鬼波剖面，其计算过程可表示为
g′(x,t)＝s(x,t)-p(x,t) (7)
其中，g′(x,t)表示由原始中提取的鬼波剖面；
在已知鬼波周期或延迟时的情况下，将公式(7)所提取的鬼波剖面转换为一次波剖面，该过程为鬼波预测的反过程；设海面的反射系数为-1，将鬼波转换为一次波的过程为
其中，p′(x,t)为表示利用鬼波记录g′(x,t)生成的一次波剖面。
进一步，所述的高信噪比剖面的合成；将衰减鬼波后的剖面与鬼波转换的一次波剖面相叠加，从而得到信噪比增强的浅剖数据，即
d(x,t)＝p(x,t)+p′(x,t) (9)
式中，d(x,t)表示信噪比得到显著增强的浅剖数据。
本发明与现有技术相比的有益效果：本发明方法实现了鬼波的预测与压制，并将提取的鬼波剖面转换为一次波剖面。通过衰减鬼波后的剖面与鬼波转换的一次波剖面相加，能够得到信噪比显著增强的浅剖数据。处理后的浅剖数据能够精确反映实际的地质构造与层位信息，并能为后续的多次波压制过程提供高精度的数据。
附图说明
图1为H海域Z1测线原始浅剖数据示例；
图2为浅剖数据自相关与拾取的鬼波周期；
图3为预测的鬼波剖面示例；
图4为衰减鬼波后的浅剖数据示例；
图5为原始数据中提取的鬼波剖面示例；
图6为鬼波剖面转换的一次波剖面示例；
图7为衰减鬼波后的剖面与鬼波转换的一次波剖面相加结果。
具体实施方式
下面通过实施例结合附图来对本发明的技术方案做进一步解释，但本发明的保护范围不受实施例任何形式上的限制。
本实施例是基于鬼波提取与转换的浅剖数据信噪比增强方法，具体的实施过程主要分为以下四步：(1)基于自相关谱分析的鬼波周期(或延迟时)提取；(2)基于时间延迟的鬼波预测；(3)基于最小平方滤波的鬼波剔除；(4)鬼波的提取与转换；(5)高信噪比剖面的合成。
实施例1
H海域为硬海底区域，海底起伏剧烈，平均水深约45m，Z测线为炮间距为1m，接收道数为1，道间距为1m；记录的采样间隔为0.1ms的浅地层剖面记录。进行浅地层剖面测量时，由地下反射回来的地震波到达海面后发生下行反射，被采集装置接收后形成周期约为5毫秒、振幅较强的鬼波干扰。鬼波同相轴与同时段的有效信号相交或重合，打断了海底以下地质层位的分布趋势，掩盖了有效的层位信息，因此需要进行鬼波压制与有效信号增强处理。
以下结合附图详细阐述本发明的具体实施过程：
(1)基于自相关谱分析的鬼波周期(或延迟时)提取。对于图1所示的浅剖数据，其中含有周期性特征的鬼波干扰，自相关谱分析是获取其准确延迟时的有效办法。为分析鬼波与一次波二者间的周期性，令最大的相关延迟时为m(要求m的值大于鬼波周期)，创建浅剖数据s(x,t)的自相关谱，可得
其中，s(x,t)表示浅剖数据，x表示剖面中各道数据的坐标，t代表旅行时；r(x,τ)表示浅剖数据的自相关谱，τ为相关延迟时。
由于海面的反射系数接近-1，在基于式(1)所创建的自相关谱r(x,τ)中，一次波与鬼波因极性相反其相关值表现为负值。因此，可通过在谱中拾取极小值的方式获得鬼波周期(或延迟时)，其可以表示为一条随坐标x变化的旅行时曲线Δτ(x)，即
Δτ(x)＝min[r(x,τ)] (2)
其中，Δτ表示鬼波的周期(或延迟时)；min表示通过人工拾取方式获得自相关谱中鬼波周期的过程。
令最大的鬼波相关延迟时为30毫秒，创建如图1所示浅剖数据的自相关谱，所得结果如图2所示。经过分析可知，谱中旅行时约为5毫秒的位置存在振幅将强的水平同相轴，其为鬼波与一次波的相关能量。随后，利用交互拾取的方式获得了鬼波周期曲线(见图2中白色曲线)。
(2)基于时间延迟的鬼波预测。与一次波相比，鬼波增加了海面位置上行与下行的传播过程，其旅行时大于与之相关的一次波。设海面的反射系数为-1，在已知鬼波周期(或延迟时)的情况下，则可利用原始数据来预测鬼波数据，时间域的鬼波预测方程为
其中，A(x)表示鬼波传播的几何扩散效应，v为海水速度，通常取1500m/s；g(x,t)为预测的鬼波剖面。
输入图1所示的浅剖数据，利用拾取的鬼波周期基于公式(3)进行鬼波预测，获得了如图3所示的鬼波记录，其中仅含有鬼波干扰，预测旅行时与原始记录中的基本一致。
(3)基于最小平方滤波的鬼波剔除。由于预测的鬼波(见图3)与原始数据(见图1)中的鬼波存在波形与振幅差异，需要采用最小平方滤波技术进行鬼波压制。对于基于L2范数的自适应相减处理，通过使下式中的误差能量平方和e值为最小来确定滤波因子
e＝||s-g*a||2 (4)
其中，a表示待求取的滤波因子；向量s与g分别代表浅剖数据s(x,t)以及预测的鬼波剖面g(x,t)；*表示褶积运算，||…||2基于L2范数的最小平方约束过程。
通过求解式(4)可获得滤波因子a，当滤波因子a的长度为1时，可避免线性方程组的求解，而直接给出滤波因子的计算公式
式中，a表示长度为1的滤波因子；向量s与g分别代表浅剖数据s(x,t)以及预测的鬼波剖面g(x,t)；上标“T”表示矩阵的转置运算。
将滤波因子a代入下式即可实现鬼波的压制，相应的计算公式为
p(x,t)＝s(x,t)-a*g(x,t) (6)
式中，p(x,t)表示剔除鬼波的浅剖数据。
输入图1所示的原始剖面与图3所示的鬼波记录，基于公式(5)求取滤波因子，然后根据公式(6)进行鬼波压制，所得结果见图4，原始剖面中的鬼波干扰得到了有效压制。
(4)鬼波的提取与转换。通过对衰减鬼波前、后的数据作差可以获得原始数据中的鬼波剖面，其计算过程可表示为
g′(x,t)＝s(x,t)-p(x,t) (7)
其中，g′(x,t)表示由原始中提取的鬼波剖面。
在已知鬼波周期(或延迟时)的情况下，可将公式(7)所提取的鬼波剖面转换为一次波剖面，该过程为鬼波预测的反过程。设海面的反射系数为-1，将鬼波转换为一次波的过程为
其中，p′(x,t)为表示利用鬼波记录g′(x,t)生成的一次波剖面。
通过对图1与图4所示衰减鬼波前、后的数据作差，提取原始剖面中的鬼波剖面(见图5)，然后根据公式(8)将其转换为一次波信号，所得结果如图6所示，其中的一次波信号与原始记录中的基本一致。
(5)高信噪比剖面的合成。为了压制浅剖数据s(x,t)中的随机噪音，提高数据的信噪比，改善剖面中同相轴的连续性，可将衰减鬼波后的剖面与鬼波转换的一次波剖面相叠加，从而得到信噪比增强的浅剖数据，即
d(x,t)＝p(x,t)+p′(x,t) (9)
式中，d(x,t)表示信噪比得到显著增强的浅剖数据。
为了改善图4所示衰减鬼波记录的信噪比，将其与鬼波转换的一次波剖面(见图5)相叠加，从而得到了如图7所示的浅剖数据，其信噪比得到明显提升，能够精确反映实际的地质构造与层位信息，并为后续的多次波压制过程提供了高精度的输入数据。
Claims (6)
1.一种基于鬼波提取与转换的浅剖数据信噪比增强方法，其特征在于，所述方法包括(1)基于自相关谱分析的鬼波周期或延迟时提取，(2)基于时间延迟的鬼波预测，(3)基于最小平方滤波的鬼波剔除，(4)鬼波的提取与转换，(5)高信噪比剖面的合成；所述的鬼波的提取与转换，通过对衰减鬼波前、后的数据作差获得原始数据中的鬼波成分，其计算过程可表示为
g′(x,t)＝s(x,t)-p(x,t) (7)
其中，g′(x,t)表示由原始中提取的鬼波剖面；
在已知鬼波周期或延迟时的情况下，将公式(7)所提取的鬼波剖面转换为一次波剖面，该过程为鬼波预测的反过程；设海面的反射系数为-1，将鬼波转换为一次波的过程为
其中，p′(x,t)为表示利用鬼波记录g′(x,t)生成的一次波剖面。
2.根据权利要求1所述的一种基于鬼波提取与转换的浅剖数据信噪比增强方法，其特征在于，所述的基于自相关谱分析的鬼波周期或延迟时提取为对于浅剖面数据中具有周期性特征的鬼波干扰，使用自相关谱分析获取其准确延迟时。
3.根据权利要求2所述的一种基于鬼波提取与转换的浅剖数据信噪比增强方法，其特征在于，分析鬼波与一次波二者间的周期性，令最大的相关延迟时为m，要求m的值大于鬼波周期，创建浅剖数据s(x,t)的自相关谱，得
其中，s(x,t)表示浅剖数据，x表示剖面中各道数据的坐标，t代表旅行时；r(x,τ)表示浅剖数据的自相关谱，τ为相关延迟时；
在基于式(1)所创建的自相关谱r(x,τ)中，通过在谱中拾取极小值的方式获得鬼波周期或延迟时，其表示为一条随坐标x变化的旅行时曲线Δτ(x)，即
Δτ(x)＝min[r(x,τ)] (2)
其中，Δτ表示鬼波的周期或延迟时；min表示通过人工拾取方式获得相关谱中鬼波周期的过程。
4.根据权利要求3所述的一种基于鬼波提取与转换的浅剖数据信噪比增强方法，其特征在于，所述的基于时间延迟的鬼波预测；设海面的反射系数为-1，在已知鬼波周期或延迟时的情况下，利用原始数据来预测鬼波数据，时间域的鬼波预测方程为
其中，A(x)表示鬼波传播的几何扩散效应，v为海水速度；g(x,t)为预测的鬼波剖面。
5.根据权利要求4所述的一种基于鬼波提取与转换的浅剖数据信噪比增强方法，其特征在于，所述的基于最小平方滤波的鬼波剔除，采用最小平方滤波技术进行鬼波压制；对于基于L2范数的自适应相减处理，通过使下式中的误差能量平方和e值为最小来确定滤波因子
e＝||s-g*a||2 (4)
其中，a为滤波因子；向量s与g分别代表浅剖数据s(x,t)以及预测的鬼波剖面g(x,t)；*表示褶积运算，||…||2基于L2范数的最小平方约束过程；
通过求解式(4)可获得滤波因子a，当滤波因子a的长度为1时，直接给出滤波因子的计算公式
式中，上标“T”表示矩阵的转置运算；
将滤波因子a代入下式实现鬼波的压制，相应的计算公式为
p(x,t)＝s(x,t)-a*g(x,t) (6)
式中，p(x,t)表示剔除鬼波的浅剖数据。
6.根据权利要求5所述的一种基于鬼波提取与转换的浅剖数据信噪比增强方法，其特征在于，所述的高信噪比剖面的合成，将衰减鬼波后的剖面与鬼波转换的一次波剖面相叠加，从而得到信噪比增强的浅剖数据，即
d(x,t)＝p(x,t)+p′(x,t) (9)
式中，d(x,t)表示信噪比得到显著增强的浅剖数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311403468.0A CN117148443B (zh) | 2023-10-27 | 2023-10-27 | 基于鬼波提取与转换的浅剖数据信噪比增强方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311403468.0A CN117148443B (zh) | 2023-10-27 | 2023-10-27 | 基于鬼波提取与转换的浅剖数据信噪比增强方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117148443A CN117148443A (zh) | 2023-12-01 |
| CN117148443B true CN117148443B (zh) | 2024-03-19 |
Family
ID=88902962
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311403468.0A Active CN117148443B (zh) | 2023-10-27 | 2023-10-27 | 基于鬼波提取与转换的浅剖数据信噪比增强方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117148443B (zh) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4794573A (en) * | 1988-02-11 | 1988-12-27 | Conoco Inc. | Process for separating upgoing and downgoing events on vertical seismic profiles |
| US6754591B1 (en) * | 2002-12-27 | 2004-06-22 | Chevron U.S.A. | Method for processing seismic traces to provide seismic data with enhanced signal-to-noise ratio |
| CN103308944A (zh) * | 2013-05-20 | 2013-09-18 | 中国石油天然气股份有限公司 | 一种海洋地震资料处理中鬼波衰减的方法及装置 |
| CN107678062A (zh) * | 2017-09-15 | 2018-02-09 | 上海海洋大学 | 双曲Radon域综合预测反褶积和反馈循环方法压制多次波模型构建方法 |
| CN108983284A (zh) * | 2018-06-22 | 2018-12-11 | 中国石油大学(华东) | 一种适用于海上斜缆数据的f-p域鬼波压制方法 |
| CN109143363A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 海洋拖缆双检采集鬼波压制方法及系统 |
| CN109975868A (zh) * | 2019-03-26 | 2019-07-05 | 中国科学院地质与地球物理研究所 | 一种基于泰勒展开的波动方程鬼波压制方法 |
| CN110967735A (zh) * | 2018-09-28 | 2020-04-07 | 中国石油化工股份有限公司 | 自适应的鬼波压制方法及系统 |
| CN112817047A (zh) * | 2020-12-31 | 2021-05-18 | 北京东方联创地球物理技术有限公司 | 海洋地震自适应去鬼波方法、装置、电子设备及介质 |
| CN114609671A (zh) * | 2020-12-08 | 2022-06-10 | 中国石油天然气股份有限公司 | 鬼波衰减方法、装置、计算机设备及可读存储介质 |
| CN115061197A (zh) * | 2022-05-27 | 2022-09-16 | 中国海洋大学 | 二维海面鬼波水体成像测量方法、系统、终端及测流设备 |
| CN115712146A (zh) * | 2022-12-02 | 2023-02-24 | 中国地质大学（武汉） | 基于频率慢度域延拓的鬼波参数最优化拖缆鬼波压制方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7961551B2 (en) * | 2008-03-21 | 2011-06-14 | Westerngeco L.L.C. | Determining directional propagation attributes of a seismic event |
| US8949030B2 (en) * | 2011-07-29 | 2015-02-03 | Westerngeco L.L.C. | Attenuating sea-surface ghost wave effects in seismic data |
| US20130028045A1 (en) * | 2011-07-29 | 2013-01-31 | Ralf Ferber | Seismic survey designs for attenuating sea-surface ghost wave effects in seismic data |
-
2023
- 2023-10-27 CN CN202311403468.0A patent/CN117148443B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4794573A (en) * | 1988-02-11 | 1988-12-27 | Conoco Inc. | Process for separating upgoing and downgoing events on vertical seismic profiles |
| US6754591B1 (en) * | 2002-12-27 | 2004-06-22 | Chevron U.S.A. | Method for processing seismic traces to provide seismic data with enhanced signal-to-noise ratio |
| CN103308944A (zh) * | 2013-05-20 | 2013-09-18 | 中国石油天然气股份有限公司 | 一种海洋地震资料处理中鬼波衰减的方法及装置 |
| CN109143363A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 海洋拖缆双检采集鬼波压制方法及系统 |
| CN107678062A (zh) * | 2017-09-15 | 2018-02-09 | 上海海洋大学 | 双曲Radon域综合预测反褶积和反馈循环方法压制多次波模型构建方法 |
| CN108983284A (zh) * | 2018-06-22 | 2018-12-11 | 中国石油大学(华东) | 一种适用于海上斜缆数据的f-p域鬼波压制方法 |
| CN110967735A (zh) * | 2018-09-28 | 2020-04-07 | 中国石油化工股份有限公司 | 自适应的鬼波压制方法及系统 |
| CN109975868A (zh) * | 2019-03-26 | 2019-07-05 | 中国科学院地质与地球物理研究所 | 一种基于泰勒展开的波动方程鬼波压制方法 |
| CN114609671A (zh) * | 2020-12-08 | 2022-06-10 | 中国石油天然气股份有限公司 | 鬼波衰减方法、装置、计算机设备及可读存储介质 |
| CN112817047A (zh) * | 2020-12-31 | 2021-05-18 | 北京东方联创地球物理技术有限公司 | 海洋地震自适应去鬼波方法、装置、电子设备及介质 |
| CN115061197A (zh) * | 2022-05-27 | 2022-09-16 | 中国海洋大学 | 二维海面鬼波水体成像测量方法、系统、终端及测流设备 |
| CN115712146A (zh) * | 2022-12-02 | 2023-02-24 | 中国地质大学（武汉） | 基于频率慢度域延拓的鬼波参数最优化拖缆鬼波压制方法 |
Non-Patent Citations (6)
| Title |
|---|
| f-p域海底电缆鬼波压制方法研究;马安;《中国优秀硕士学位论文全文数据库基础科学辑》(2023年第02期);全文 * |
| 中浅地层剖面质量改进方法及应用;陶华;李彦杰;;海洋技术学报(第05期);全文 * |
| 子波法去鬼波在墨西哥湾的应用;姜丹;蒲晓东;麻志国;石孟常;李君君;;物探与化探(第05期);全文 * |
| 拉东域基于最小二乘反演的斜缆数据鬼波压制方法;张威;韩立国;李洪建;叶林;张齐;;世界地质(第02期);全文 * |
| 斜缆数据tau-p域鬼波压制算法研究;郭毓;《中国优秀硕士学位论文全文数据库基础科学辑》(2020年第07期);全文 * |
| 鬼波衰减在深水多次波衰减中的作用;李慧龙;王慕雨;王子秋;张广智;石孟常;麻志国;;地球物理学进展(第06期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117148443A (zh) | 2023-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108983284B (zh) | 一种适用于海上斜缆数据的f-p域鬼波压制方法 | |
| CN113625337B (zh) | 一种极浅水高精度地震资料快速成像方法 | |
| CA2357604A1 (en) | Method for identifying and removing multiples from seismic reflection data | |
| CN112946732B (zh) | 海上立体观测系统联合压制多次波单缆处理方法及系统 | |
| CN112327362B (zh) | 速度域的海底多次波预测与追踪衰减方法 | |
| MX2014001334A (es) | Sistemas y metodos para detectar ruido del oleaje en un recolector sismico. | |
| CN108985304A (zh) | 一种基于浅剖数据的沉积层结构自动提取方法 | |
| US10488514B2 (en) | Signal processing method of raw data obtained by using chirp sub-bottom profiler and readable storage medium storing the method | |
| CN115905805A (zh) | 一种基于全局信息判别gan的das数据多尺度噪声消减方法 | |
| CN107450103B (zh) | 一种基于边界积分反算子的鬼波压制方法 | |
| Vesnaver et al. | Chirp data processing for fluid flow detection at the Gulf of Trieste (Northern Adriatic Sea) | |
| CN110967735A (zh) | 自适应的鬼波压制方法及系统 | |
| US10310115B2 (en) | Prediction and subtraction of multiple diffractions | |
| CN117148443B (zh) | 基于鬼波提取与转换的浅剖数据信噪比增强方法 | |
| Huo et al. | Suppressing reflected guided waves from contaminated borehole radar data | |
| CN110261899B (zh) | 地震数据z字形干扰波去除方法 | |
| EP0689064B1 (en) | Methods of adaptively reducing noise transients in co-sensor seismic signals | |
| CN114442152B (zh) | 地震数据的海底多次波预测方法及装置 | |
| CN114185095B (zh) | 一种三维平面波域地震数据多次波压制的方法 | |
| KR101864307B1 (ko) | 그래디언트 분석을 통한 탄성파 너울보정 및 지하지층구조 탐사 자료 처리방법 및 그 기록매체 | |
| CN112578453B (zh) | 多次波压制方法及系统 | |
| CN117148432B (zh) | 基于多尺度分量提取的浅剖数据空间插值方法 | |
| CN118655622B (zh) | 浅水地震数据叠前时间偏移处理方法、装置、设备及介质 | |
| CN113514889A (zh) | 一种提升海洋深反射地震数据中低频信号能量的处理方法 | |
| Nasif | Processing and joint interpretation of multi-resolution marine seismic datasets |
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 |