CN107085239B - 海底电缆双检资料陆检微分合并方法 - Google Patents
海底电缆双检资料陆检微分合并方法 Download PDFInfo
- Publication number
- CN107085239B CN107085239B CN201710315194.8A CN201710315194A CN107085239B CN 107085239 B CN107085239 B CN 107085239B CN 201710315194 A CN201710315194 A CN 201710315194A CN 107085239 B CN107085239 B CN 107085239B
- Authority
- CN
- China
- Prior art keywords
- land
- data
- inspection
- differential
- seismic channel
- 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 39
- 238000001514 detection method Methods 0.000 title claims abstract description 22
- 238000007689 inspection Methods 0.000 claims abstract description 94
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 49
- 230000008859 change Effects 0.000 claims abstract description 15
- 230000005012 migration Effects 0.000 claims abstract description 15
- 238000013508 migration Methods 0.000 claims abstract description 15
- 238000009795 derivation Methods 0.000 claims abstract description 3
- 230000008569 process Effects 0.000 claims abstract description 3
- 239000011435 rock Substances 0.000 claims description 6
- 238000009792 diffusion process Methods 0.000 claims description 3
- 239000013535 sea water Substances 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 9
- 238000001228 spectrum Methods 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 3
- 230000004069 differentiation Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000009466 transformation Effects 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
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/14—Signal detection
- G01V2210/142—Receiver location
- G01V2210/1425—Land surface
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/14—Signal detection
- G01V2210/142—Receiver location
- G01V2210/1427—Sea bed
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/52—Move-out correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/59—Other corrections
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
本发明提供一种海底电缆双检资料陆检微分合并方法，包括：对海底电缆双检采集地震资料的水检地震道集和陆检地震道集进行预处理；根据工区的实际地质情况，利用地质经验和认识估计所在工区的海底反射系数r；将陆检资料Z的每个地震道对时间t进行求导，即得到陆检微分资料；扫描水检资料和陆检微分资料近偏移距地震道的振幅信息；将所有参与运算的近偏移距地震道振幅扫描信息进行平均，获得尺度变换因子k；将水检资料与陆检微分资料根据尺度变换因子k与海底反射系数r进行合并叠加。该海底电缆双检资料陆检微分合并方法能够在保证有效波不受损害的前提下，更好地压制海底鸣震干扰，能够达到拓宽频带、提高分辨率的目的。
Description
技术领域
本发明涉及油气勘探地震数据处理领域，特别是涉及到一种海底电缆双检资料陆检微分合并方法。
背景技术
在海洋环境中，海上拖缆是高效率获得地震数据的一种方法，然而这种方法在障碍物附近、浅水处的施工能力受到限制，同时也限制了观测系统的多变性。海底电缆双检接收技术是针对滩浅海地区勘探的地震采集新技术，它将压电检波器和速度检波器结合起来进行地震信号采集，在障碍物附近保持炮检距和方位角分布的一致性，且有多变的采集观测系统的能力，并能在水深从很浅到100米的范围内施工，大大提高了地震资料的采集质量，对海上地震勘探常见的鸣震具有较好的压制效果，从而提高了地震资料解释精度。目前，海底电缆双检采集数据在滩浅海地震勘探中显示出极大优越性，取得了很好的勘探效果。近年来，由于对这种双检数据的相位、振幅以及频率等差异认识不足，一般根据水检和陆检资料对有效波和鸣震的不同响应，采用直接合并或者调90度相位合并的方式进行双检合并以压制海底鸣震，直接合并的方法，没有考虑双检资料的相位和频率差异，使得有效波不能同相叠加，受到一定程度的损害，鸣震不能彻底消除；调90度相位合并的方法，虽然考虑了双检资料的相位差异，但是没有考虑其频率差异，造成叠加剖面分辨率和频率成分不一致，影响后续的偏移成像以及综合研究。为此我们发明了一种新的海底电缆双检资料陆检微分合并方法，解决了以上技术问题。
发明内容
本发明的目的是提供一种适用于海底电缆双检地震采集数据用于实现水检和陆检资料合并以压制海底鸣震干扰的海底电缆双检资料陆检微分合并方法。
本发明的目的可通过如下技术措施来实现：海底电缆双检资料陆检微分合并方法，该海底电缆双检资料陆检微分合并方法包括：步骤1，对海底电缆双检采集地震资料的水检地震道集和陆检地震道集进行预处理；步骤2，根据工区的实际地质情况，利用地质经验和认识估计所在工区的海底反射系数r；步骤3，将陆检资料Z的每个地震道对时间t进行求导，即得到陆检微分资料；步骤4，扫描水检资料和陆检微分资料近偏移距地震道的振幅信息；步骤5，将所有参与运算的近偏移距地震道振幅扫描信息进行平均，获得尺度变换因子k；步骤6，将水检资料与陆检微分资料根据尺度变换因子k与海底反射系数r进行合并叠加。
本发明的目的还可通过如下技术措施来实现：
在步骤1，进行水检和陆检地震道集的分选、线性动校正、二次定位、球面扩散补偿、去噪以及地表一致性补偿这些处理，得到预处理后的水检和陆检的地震单炮数据。
在步骤2，根据工区的实际地质情况，利用地质经验和认识，利用岩石密度和地震波速度信息估计所在工区的海底反射系数r，其估算公式为：
其中，ρ0和v0分别为海水密度和地震波速度，ρ1和v1分别为海底岩层的密度和地震波速度。
在步骤4，扫描水检单炮资料和陆检微分单炮资料近偏移距地震道的振幅信息，要求参与扫描的地震单炮资料偏移距小于500米，获得一系列的水检单炮振幅信息Ai和陆检微分单炮振幅信息Bi，i＝1,2,…,n，n为偏移距小于500米的地震道道数。
在步骤5，将小于500米的近偏移距地震道水检和陆检微分单炮振幅扫描信息Ai和Bi进行平均得到Aavg和Bavg，获得尺度变换因子
在步骤6，将水检资料P与陆检微分资料根据尺度变换因子k与海底反射系数r进行合并叠加，合并公式为：
本发明中的海底电缆双检资料陆检微分合并方法，依据水检和陆检资料对有效波和鸣震的不同响应以及所记录的地球物理量(水检记录的是加速度分量，陆检记录的是速度分量)的差异，首先通过微分运算将陆检资料转换到加速度域，此时水检资料和陆检微分结果在同一域中，代表同样的地球物理量，然后通过调整尺度变换因子，将水检资料同陆检微分结果进行合并得到叠加剖面，这种方法充分考虑了双检资料中有效波与鸣震在振幅、相位以及频率成分的差异，能够在保证有效波不受损害的前提下，更好地压制海底鸣震干扰，得到拓频的叠加结果。该发明充分考虑了双检资料中有效波与鸣震在振幅、相位以及频率成分方面的差异，能够在保证有效波不受损害的前提下，更好地压制海底鸣震干扰；该发明利用双检接收地震资料中的水检和陆检微分进行合并，能够达到拓宽频带、提高分辨率的目的。
附图说明
图1为本发明的海底电缆双检资料陆检微分合并方法的一具体实施例的流程图；
图2为本发明的一具体实施例中预处理后的水检和陆检地震单炮数据的示意图；
图3为本发明的一具体实施例中陆检微分资料的示意图；
图4为本发明的一具体实施例中陆检微分合并叠加剖面的示意图；
图5为本发明的一具体实施例中直接合并叠加剖面的示意图；
图6为本发明的一具体实施例中调90度相位合并叠加剖面的示意图；
图7为本发明的一具体实施例中陆检微分合并方法中水检与合并剖面残差的示意图；
图8为本发明的一具体实施例中陆检微分合并方法中水检与合并剖面残差频谱的示意图；
图9为本发明的一具体实施例中直接合并方法中水检与合并剖面残差的示意图；
图10为本发明的一具体实施例中直接合并方法中水检与合并剖面残差频谱的示意图；
图11为本发明的一具体实施例中调90度相位合并方法中水检与合并剖面残差的示意图；
图12为本发明的一具体实施例中调90度相位合并方法中水检与合并剖面残差频谱的示意图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合附图所示，作详细说明如下。
如图1所示，图1为本发明的海底电缆双检资料陆检微分合并方法的流程图。本发明的海底电缆双检资料陆检微分合并方法，首先通过微分运算将陆检资料转换到加速度域，此时水检资料和陆检微分结果在同一域中，代表同样的地球物理量，然后通过调整尺度变换因子，将水检资料同陆检微分结果进行合并得到叠加剖面。包括以下步骤：
在步骤101，对海底电缆双检采集地震资料的水检地震道集和陆检地震道集进行预处理，主要进行水检和陆检地震道集的分选、线性动校正、二次定位、球面扩散补偿、去噪以及地表一致性补偿等处理，得到预处理后的水检和陆检的地震单炮数据。
在步骤102，根据工区的实际地质情况，利用地质经验和认识，利用岩石密度和地震波速度信息估计所在工区的海底反射系数r，其估算公式为：
其中，ρ0和v0分别为海水密度和地震波速度，ρ1和v1分别为海底岩层的密度和地震波速度。
在步骤103，将预处理后的陆检单炮数据Z的每个地震道对时间t进行求导运算，即：得到陆检微分单炮资料。
在步骤104，扫描水检单炮资料和陆检微分单炮资料近偏移距地震道的振幅信息，要求参与扫描的地震单炮资料偏移距小于500米，获得一系列的水检单炮振幅信息Ai和陆检微分单炮振幅信息Bi，i＝1,2,…,n，n为偏移距小于500米的地震道道数；
在步骤105，将小于500米的近偏移距地震道水检和陆检微分单炮振幅扫描信息Ai和Bi进行平均得到Aavg和Bavg，获得尺度变换因子
在步骤106，将水检资料与陆检微分资料根据尺度变换因子k与海底反射系数r进行合并叠加，合并公式为：
在应用本发明的一具体实施例中，包括以下步骤：
步骤1，输入预处理后的水检和陆检地震单炮数据(图2)；
步骤2，根据该探区的地质认识估算海底反射系数为0.45；
步骤3，将陆检资料的每个地震道对时间进行求导运算，得到陆检微分资料(图3)；
步骤4，扫描水检资料和陆检微分资料近偏移距地震道的振幅信息，此次参与扫描的地震资料偏移距小于500米；
步骤5，将所有参与运算的近偏移距地震道振幅扫描信息进行平均，获得该资料的尺度变换因子为0.00002；
步骤6，将水检资料与陆检微分资料根据以上步骤求取的尺度变换因子与海底反射系数利用公式进行合并叠加，得到陆检微分合并叠加结果(图4)。将陆检微分合并方法的叠加剖面(图4)与直接合并方法(图5)以及调90度相位合并方法的叠加剖面(图6)进行相比，可以看到陆检微分合并法鸣震去除效果更好。将三种方法的水检剖面与合并剖面相减可以得到残差剖面，对残差剖面进行频谱分析，效果显示在图7到图12中，可以看到陆检微分合并方法的残差频谱只要集中在高频部分(鸣震)，而其他两种方法的残差频谱还包含低频成分(有效波)，这也证明了本发明的陆检微分合并方法在压制鸣震的时候，能够保证有效波不受损害，优于其他两种合并方法。
Claims (2)
1.海底电缆双检资料陆检微分合并方法，其特征在于，该海底电缆双检资料陆检微分合并方法包括：
步骤1，对海底电缆双检采集地震资料的水检地震道集和陆检地震道集进行预处理；
步骤2，根据工区的实际地质情况，利用地质经验和认识估计所在工区的海底反射系数r；
步骤3，将陆检资料Z的每个地震道对时间t进行求导，即得到陆检微分资料；
步骤4，扫描水检资料和陆检微分资料近偏移距地震道的振幅信息；
步骤5，将所有参与运算的近偏移距地震道振幅扫描信息进行平均，获得尺度变换因子k；
步骤6，将水检资料与陆检微分资料根据尺度变换因子k与海底反射系数r进行合并叠加；
在步骤2，根据工区的实际地质情况，利用地质经验和认识，利用岩石密度和地震波速度信息估计所在工区的海底反射系数r，其估算公式为：
其中，ρ0和v0分别为海水密度和地震波速度，ρ1和v1分别为海底岩层的密度和地震波速度；
在步骤5，将小于500米的近偏移距地震道水检和陆检微分单炮振幅扫描信息Ai和Bi进行平均得到Aavg和Bavg，获得尺度变换因子i＝1,2,…,n，n为偏移距小于500米的地震道道数；
在步骤6，将水检资料P与陆检微分资料根据尺度变换因子k与海底反射系数r进行合并叠加，合并公式为：
在步骤1，进行水检和陆检地震道集的分选、线性动校正、二次定位、球面扩散补偿、去噪以及地表一致性补偿这些处理，得到预处理后的水检和陆检的地震单炮数据。
2.根据权利要求1所述的海底电缆双检资料陆检微分合并方法，其特征在于，在步骤4，扫描水检单炮资料和陆检微分单炮资料近偏移距地震道的振幅信息，要求参与扫描的地震单炮资料偏移距小于500米，获得一系列的水检单炮振幅信息Ai和陆检微分单炮振幅信息Bi，i＝1,2,…,n，n为偏移距小于500米的地震道道数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710315194.8A CN107085239B (zh) | 2017-05-05 | 2017-05-05 | 海底电缆双检资料陆检微分合并方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710315194.8A CN107085239B (zh) | 2017-05-05 | 2017-05-05 | 海底电缆双检资料陆检微分合并方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107085239A CN107085239A (zh) | 2017-08-22 |
| CN107085239B true CN107085239B (zh) | 2019-08-09 |
Family
ID=59612715
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710315194.8A Active CN107085239B (zh) | 2017-05-05 | 2017-05-05 | 海底电缆双检资料陆检微分合并方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107085239B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108303741A (zh) * | 2018-02-11 | 2018-07-20 | 中国石油化工股份有限公司 | 基于水检积分的地震双检合并鸣震压制方法 |
| CN110208852B (zh) * | 2019-05-15 | 2020-06-02 | 中国地质大学（北京） | 陆地双检原始地震数据评价及数据处理质控方法 |
| CN112147716A (zh) * | 2020-09-15 | 2020-12-29 | 中海油田服务股份有限公司 | 一种水陆检信号的合并方法和装置 |
| CN112327360B (zh) * | 2020-10-29 | 2022-12-09 | 中海油田服务股份有限公司 | 一种实现噪声处理的方法、装置、计算机存储介质及终端 |
| CN112415595B (zh) * | 2020-10-30 | 2024-05-28 | 中国石油天然气集团有限公司 | 海底节点二次定位的方法及装置 |
| CN113075733B (zh) * | 2021-03-26 | 2022-03-29 | 成都理工大学 | 一种水陆双检的道序质控方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101576621A (zh) * | 2008-05-07 | 2009-11-11 | 王振华 | 海底电缆双检地震勘探的数据处理方法及数据处理装置 |
-
2017
- 2017-05-05 CN CN201710315194.8A patent/CN107085239B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101576621A (zh) * | 2008-05-07 | 2009-11-11 | 王振华 | 海底电缆双检地震勘探的数据处理方法及数据处理装置 |
Non-Patent Citations (4)
| Title |
|---|
| Numerical simulation and feature analysis of OBC dual-sensor data;Ning Qin et al.;《2017 CGS/SEG International Geophysical Conference》;20170420;第874-877页 * |
| 双检资料处理方法研究;崔荣富等;《内蒙古石油化工》;20121231(第11期);第7-8页 * |
| 海底电缆双检资料消除鸣震处理技术;熊小娟;《中国石油和化工标准与质量》;20110115(第1期);第53页 * |
| 速度检波器与压电检波器相位差异分析及应用;任立刚等;《地球物理学进展》;20150228;第30卷(第1期);第455-457页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107085239A (zh) | 2017-08-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107085239B (zh) | 海底电缆双检资料陆检微分合并方法 | |
| Schwenk et al. | Surface-wave methods for anomaly detection | |
| She et al. | Shallow crustal structure of the middle‐lower Yangtze River region in eastern China from surface‐wave tomography of a large volume airgun‐shot experiment | |
| Kimman et al. | Characteristics of seismic noise: fundamental and higher mode energy observed in the northeast of the Netherlands | |
| CN106896409B (zh) | 一种基于波动方程边值反演的变深度缆鬼波压制方法 | |
| NO20141128A1 (no) | Konstruksjon og anvendelse av vinkelsamlere fra tredimensjonal avbildning av multippelbølgefelter | |
| CN107884829A (zh) | 一种联合压制浅海obc地震资料多次波的方法 | |
| Alaei et al. | Seismic imaging of fault damaged zone and its scaling relation with displacement | |
| CN107367761B (zh) | 一种大沙漠区深层碳酸盐岩储层广角反射波成像方法 | |
| Jervis et al. | Making time-lapse seismic work in a complex desert environment for CO2 EOR monitoring—Design and acquisition | |
| Colombo et al. | Near-surface full-waveform inversion in a transmission surface-consistent scheme | |
| Lee et al. | Resolution analysis of shallow marine seismic data acquired using an airgun and an 8-channel streamer cable | |
| Luo et al. | Joint Inversion of Active Sources and Ambient Noise for Near‐Surface Structures: A Case Study in the Balikun Basin, China | |
| CN108303741A (zh) | 基于水检积分的地震双检合并鸣震压制方法 | |
| Reilly et al. | The case for separate sensor processing: Meeting the imaging challenge in a producing carbonate field in the Middle East | |
| Lee et al. | High-resolution shallow marine seismic surveys off Busan and Pohang, Korea, using a small-scale multichannel system | |
| Benfield et al. | Extracting geologic information directly from high-resolution full-waveform inversion velocity models—A case study from offshore Trinidad | |
| Schapper et al. | Anisotropic velocities and offset vector tile prestack-migration processing of the Durham Ranch 3D, Northwest Colorado | |
| Fu et al. | Imaging the ice sheet and uppermost crustal structures with a dense linear seismic array in the Larsemann Hills, Prydz Bay, East Antarctica | |
| CN105676278B (zh) | 一种处理地震波数据的方法和装置 | |
| Faggetter et al. | Time-lapse imaging using 3D ultra-high-frequency marine seismic reflection data | |
| Dorn et al. | Automatic fault extraction in hard and soft rock environments | |
| Tyiasning et al. | Anisotropy signatures in the Cooper Basin of Australia: Stress versus fractures | |
| Wen et al. | Highlighting display of geologic bodies based on directivity filtering | |
| Wang et al. | Research on vertical cable seismic multiples imaging processing |
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 |