CN109143363B - 海洋拖缆双检采集鬼波压制方法及系统 - Google Patents
海洋拖缆双检采集鬼波压制方法及系统 Download PDFInfo
- Publication number
- CN109143363B CN109143363B CN201710467201.6A CN201710467201A CN109143363B CN 109143363 B CN109143363 B CN 109143363B CN 201710467201 A CN201710467201 A CN 201710467201A CN 109143363 B CN109143363 B CN 109143363B
- Authority
- CN
- China
- Prior art keywords
- wave field
- ghost
- field signal
- pressure
- wave
- 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 22
- 230000001629 suppression Effects 0.000 title claims description 11
- 238000001514 detection method Methods 0.000 title description 26
- 238000007689 inspection Methods 0.000 claims abstract description 35
- 230000006835 compression Effects 0.000 claims abstract description 9
- 238000007906 compression Methods 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 14
- 238000012876 topography Methods 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 description 17
- 230000000694 effects Effects 0.000 description 7
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
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/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提出了一种海洋拖缆双检采集鬼波压制方法及系统,所述方法包括:获取压力波场信号和垂直速度波场信号;建立无鬼波上行波场信号与压力波场信号的关系;根据牛顿第二定律,建立压力波场信号和垂直速度波场信号之间的关系;基于最小平方算法,合并压力波场和垂直速度波场,消除鬼波信号。本发明采用基于最小平方的双检数据合并方法,可实现水、陆检数据合并,不依赖地下假设信息,稳定性好,受采集环境等制约少,能够适应复杂海洋地形和地质情况。
Description
技术领域
本发明属于海洋地震宽频采集领域,具体涉及基于最小平方算法的海洋拖缆双检采集鬼波压制方法。
背景技术
区别于海洋常规采集,拖缆双检采集技术在电缆的同一位置,同时拖带水检和陆检两套检波器,分别用于接收压力场和垂直速度场信号,双检采集作为唯一一种同时接收压力波场和垂直速度波场的采集方式,数据具有很好说服力。
水、陆检波器由于对地震信号具有不同的滤波,信号在振幅、相位、频率等方面存在差异,数据合并算法主要部分为振幅匹配,常规算法为陆检数据乘以匹配算子,再与水检数据求和。但由于两种检波器自身性能不同和外部条件影响,再加上海底反射系数也是在不断变化的,直接匹配求和对复杂地下结构适应性较差,无法取得理想合并效果,需要过多的人工干预。
因此对于海洋拖缆双检采集技术,其合并算法的精度和适用性决定了采集技术的应用范围和效果,常规的双检合并技术对地下假设信息依赖性强,稳定性差,受采集环境等制约因素明显,宽频采集效果难以保证,不利于技术普及推广,制约海洋拖缆双检采集技术的实际应用。优化海洋拖缆双检采集数据合并算法,利于消除鬼波影响,促进该技术的实际应用。
发明内容
本发明针对海洋拖缆双检采集原始数据,基于最小平方算法,合并处理水陆、检信号,压制鬼波,以获取海洋宽频地震原始资料,用于解决海洋地震采集中鬼波影响严重、低频缺失、陷波现象等问题。
根据本发明的一个方面,提供一种海洋拖缆双检采集鬼波压制方法,包括:
获取压力波场信号和垂直速度波场信号;
建立鬼波无鬼波上行波场信号与压力波场信号的关系;
根据牛顿第二定律,建立压力波场信号和垂直速度波场信号之间的关系;
基于最小平方算法,合并压力波场和垂直速度波场,消除鬼波信号。进一步地,压力波场信号P包括海底一次反射信号及海平面反射鬼波信号:
假设P0(ω,kx,0)为海面处接收到的压力波场,则无鬼波上行波场信号Pup(ω,kx,z)可以表示为:
建立无鬼波上行波场信号与压力波场信号的关系为:
其中kz为垂直波数,Z为深度。
进一步地,根据牛顿第二定律,压力波场信号p和垂直速度波场信号vz之间的关系为:
在频率波数域进行计算:
其中,Vz(ω,kx,z)表示频率波数域垂直速度波场信号,ρ表示密度,ω表示频率,k表示波数,z表示拖缆深度。进一步地,进一步包括建立无鬼波信号与压力波场信号和垂直速度波场信号的关系,由方程(3)和(6),可得:
由以上可知:
进一步地,基于最小平方算法,得到:
海洋拖缆双检采集技术可以直接获得压力波场P(ω,kx,z)和垂直速度波场Vz(ω,kx,z),因此通过上式可以直接获得宽频无鬼波上行波场Pup(ω,kx,z)。
经过双检采集压力波场和垂直速度波场合并,即可消除鬼波,获得宽频地震数据。
根据本发明另一方面,提供一种海洋拖缆双检采集鬼波压制系统,包括:
存储器,存储有可执行指令;
处理器,运行所述存储器存储的所述可执行指令,执行以下步骤:
获取压力波场信号和垂直速度波场信号;
建立鬼波无鬼波上行波场信号与压力波场信号的关系;
根据牛顿第二定律,建立压力波场信号和垂直速度波场信号之间的关系;
基于最小平方算法,合并压力波场和垂直速度波场,消除鬼波信号。根据本发明的又一个方面,提供一种记录介质,所述记录介质中存储有计算机可执行指令;所述计算机可执行指令运行时,执行以下步骤:
获取压力波场信号和垂直速度波场信号;
建立鬼波无鬼波上行波场信号与压力波场信号的关系;
根据牛顿第二定律,建立压力波场信号和垂直速度波场信号之间的关系;
基于最小平方算法,合并压力波场和垂直速度波场,消除鬼波信号。
本发明采用基于最小平方的双检数据合并方法,可实现水、陆检数据合并,不依赖地下假设信息,稳定性好,受采集环境等制约少,能够适应复杂海洋地形和地质情况。
附图说明
通过结合附图对本公开示例性实施方式进行更详细的描述,本公开的上述以及其它目的、特征和优势将变得更加明显,其中,在本公开示例性实施方式中,相同的参考标号通常代表相同部件。
图1示出了海洋双检采集技术水检信号叠加剖面。
图2示出了海洋双检采集技术叠加剖面。
图3示出了频谱对比。
图4示出了本发明的海洋拖缆双检采集鬼波压制方法流程图。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式,然而应该理解,可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本公开更加透彻和完整,并且能够将本公开的范围完整地传达给本领域的技术人员。
海洋拖缆双检采集中,鬼波对一次波的波形及频带宽度产生极大的影响,鬼波压制是提高海上地震资料分辨率、保真度及宽频采集效果的重要因素。
海洋拖缆双检采集技术的核心部分为宽频原始数据的采集和水、陆检数据的合并处理,双检数据的合并方法选取不当并不能实现鬼波的有效压制,低频缺失、陷波现象等问题并无法得到有效消除。因此研究最优的水陆检数据合并算法可以发挥拖缆双检宽频采集技术的最佳效果,从而提高低频能量,压制陷波,拓宽原始资料频带。
对于海洋拖缆双检采集技术数据处理,本发明采用基于最小平方的双检数据合并方法,能够适应复杂海洋地形和地质情况。实际资料处理结果表明,该技术在很好的去除鬼波的同时极大的拓宽了地震资料的频带,尤其是低频端,有利于后续资料的处理解释,且能有效压制随机噪声,提高信噪比。
如图4所示,本公开提出了一种海洋拖缆双检采集鬼波压制方法,包括:
获取压力波场信号和垂直速度波场信号;
建立鬼波无鬼波上行波场信号与压力波场信号的关系;
根据牛顿第二定律,建立压力波场信号和垂直速度波场信号之间的关系;
基于最小平方算法,合并压力波场和垂直速度波场,消除鬼波信号。
对于海洋拖缆采集,检波器信号由两部分组成:海底一次反射及海平面反射鬼波:
假设P0(ω,kx,0)为海面处接收到的压力波场,根据波场延拓原理,则无鬼波上行波场信号Pup(ω,kx,z)可以表示为:
即:
根据牛顿第二定律,压力波场p和垂直速度波场vz之间的关系为:
在频率波数域进行计算:
即:
其中,Vz(ω,kx,z)表示频率波数域垂直速度波场信号,ρ表示密度,ω表示频率,k表示波数,z表示拖缆深度。由方程3和6,可得:
由以上可知:
基于最小平方算法,可得:
海洋拖缆双检采集技术可以直接获得压力波场P(ω,kx,z)和垂直速度波场Vz(ω,kx,z),因此通过上式可以直接获得宽频无鬼波上行波场Pup(ω,kx,z)。
经过双检采集压力波场和垂直速度波场合并,即可消除鬼波,获得宽频地震数据。
根据本发明另一方面,提供一种基于最小平方算法的海洋拖缆双检采集鬼波压制系统,包括:
存储器,存储有可执行指令;
处理器,运行所述存储器存储的所述可执行指令,执行以下步骤:
获取压力波场信号和垂直速度波场信号;
建立鬼波无鬼波上行波场信号与压力波场信号的关系;
根据牛顿第二定律,建立压力波场信号和垂直速度波场信号之间的关系;
基于最小平方算法,合并压力波场和垂直速度波场,消除鬼波信号。
本发明采用基于最小平方的双检数据合并方法,可实现水、陆检数据合并,不依赖地下假设信息,稳定性好,受采集环境等制约少,能够适应复杂海洋地形和地质情况。
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
为了验证基于最小平方拖缆双检鬼波压制方法的实用效果,我们选取海上双检地震资料进行处理和分析,数据为海上水检、陆检地震记录,震源深度15m,拖缆深度25m。
基于最小平方算法进行水检、陆检数据合并,图1为海洋双检采集技术水检信号叠加剖面,图2为海洋双检采集技术叠加剖面,可见,经过鬼波压制,地震记录分辨率明显提高,由鬼波产生的虚假同相轴被很好的去除。
通过频谱分析,如图3,同样可以看出,水检记录在30hz频率附近存在明显的陷波,经过水陆检合并,陷波处能量得到有效提升,且低频能量有所提高,地震信号频带得到有效拓宽。
本发明采用基于最小平方的双检数据合并方法,能够适应复杂海洋地形和地质情况,实际资料处理结果表明,基于最小平方的海洋拖缆双检采集鬼波压制技术在很好的去除鬼波的同时极大的拓宽了地震资料的频带,有利于后续资料的处理解释。
以上已经描述了本公开的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (6)
4.一种海洋拖缆双检采集鬼波压制系统,其特征在于,包括:
存储器,存储有可执行指令;
处理器,运行所述存储器存储的所述可执行指令,执行以下步骤:
获取压力波场信号和垂直速度波场信号;
建立鬼波无鬼波上行波场信号与压力波场信号的关系;
根据牛顿第二定律,建立压力波场信号和垂直速度波场信号之间的关系;
基于最小平方算法,合并压力波场和垂直速度波场,消除鬼波信号;
其中,根据牛顿第二定律,压力波场信号p和垂直速度波场信号vz之间的关系为:
在频率波数域进行计算:
其中,Vz(ω,kx,z)表示频率波数域垂直速度波场信号,ρ表示密度,ω表示频率,k表示波数,z表示拖缆深度,P0(ω,kx,0)为海面处接收到的压力波场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710467201.6A CN109143363B (zh) | 2017-06-19 | 2017-06-19 | 海洋拖缆双检采集鬼波压制方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710467201.6A CN109143363B (zh) | 2017-06-19 | 2017-06-19 | 海洋拖缆双检采集鬼波压制方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109143363A CN109143363A (zh) | 2019-01-04 |
CN109143363B true CN109143363B (zh) | 2021-05-25 |
Family
ID=64804595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710467201.6A Active CN109143363B (zh) | 2017-06-19 | 2017-06-19 | 海洋拖缆双检采集鬼波压制方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109143363B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110208852B (zh) * | 2019-05-15 | 2020-06-02 | 中国地质大学(北京) | 陆地双检原始地震数据评价及数据处理质控方法 |
CN112444881B (zh) * | 2019-08-30 | 2024-04-09 | 中国石油化工股份有限公司 | 一种鬼波压制方法 |
CN112578454B (zh) * | 2019-09-29 | 2024-03-22 | 中国石油化工股份有限公司 | 一种任意弯曲拖缆鬼波压制方法及系统 |
CN111538088B (zh) * | 2020-05-07 | 2021-07-27 | 中国科学院地质与地球物理研究所 | 一种海上斜缆波场校正方法 |
CN112130208B (zh) * | 2020-09-03 | 2024-02-06 | 中海油田服务股份有限公司 | 一种处理电缆中鬼波的方法和装置 |
CN112147716A (zh) * | 2020-09-15 | 2020-12-29 | 中海油田服务股份有限公司 | 一种水陆检信号的合并方法和装置 |
CN117148443B (zh) * | 2023-10-27 | 2024-03-19 | 胜利信科(山东)勘察测绘有限公司 | 基于鬼波提取与转换的浅剖数据信噪比增强方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8693282B2 (en) * | 2010-05-25 | 2014-04-08 | Westerngeco L.L.C. | Deghosting seismic data |
CA2890630A1 (en) * | 2012-11-09 | 2014-05-15 | Shell Internationale Research Maatschappij B.V. | Method for deghosting seismic data acquired by a marine seismic source and receiver assembly |
EP3014307A4 (en) * | 2013-06-25 | 2017-04-05 | Westerngeco LLC | Seismic wavefield deghosting and noise attenuation |
CN103901473B (zh) * | 2014-04-14 | 2016-10-05 | 中国海洋石油总公司 | 一种基于非高斯性最大化的双检信号上下行波场分离方法 |
US10107929B2 (en) * | 2014-12-18 | 2018-10-23 | Pgs Geophysical As | Methods and systems to determine ghost operators from marine seismic data |
US10156651B2 (en) * | 2015-04-21 | 2018-12-18 | Pgs Geophysical As | Methods and systems to source deghost and source designature marine seismic data |
-
2017
- 2017-06-19 CN CN201710467201.6A patent/CN109143363B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109143363A (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109143363B (zh) | 海洋拖缆双检采集鬼波压制方法及系统 | |
Yu et al. | Attenuation of noise and simultaneous source interference using wavelet denoising | |
US10224900B2 (en) | Systems and methods for de-noising seismic data | |
US9581709B2 (en) | Suppressing 4D-noise by weighted stacking of simultaneously acquired wave-fields | |
AU2014280832B2 (en) | Seismic data spectrum restoring and broadening | |
WO2010082126A2 (en) | Processing seismic data | |
GB2449564A (en) | Suppressing ghosts in marine seismic data acquired from over/under streamers at independent x,y positions | |
Tognarelli et al. | High-resolution coherency functionals for velocity analysis: An application for subbasalt seismic exploration | |
NO337134B1 (no) | Fremgangsmåte for å prosessere seismiske data svarende til innsamlinger oppnådd for samme sone ved hjelp av seismiske mottakere plassert på havbunnen og ved hjelp av mottakere plassert ved havoverflaten | |
CN109975873B (zh) | 一种逆时偏移成像去除低频噪音的方法及系统 | |
CN110703332B (zh) | 一种鬼波压制方法 | |
CN108919345B (zh) | 一种海底电缆陆检噪声的衰减方法 | |
Bunting et al. | Marine broadband case study offshore China | |
CN111352158B (zh) | 地震信号增强方法及装置 | |
Shu et al. | Simultaneous sources deblending via frequency-varying filter in cross-spread azimuth-offset domain | |
AU2011248910A1 (en) | System and method for compensating time and offset varying near-surface effects in seismic data background | |
D'Agosto et al. | Modeling and removal of ground roll from horizontal component of C-waves | |
CN109471170B (zh) | 一种地震数据处理方法和装置 | |
Basu et al. | Pre-conditioning of data before PZ summation in OBC survey–a case study | |
Ogagarue et al. | An improved dual sensor summation method with application to four-component (4-C) seafloor seismic data from the niger delta | |
CN116088055A (zh) | 用于海洋obn数据的上下行波场分离的方法 | |
Soudani et al. | 3D methodology for OBC pre-processing | |
Geng et al. | Some Factors Impact on Deghosting by Over/Under Towed-streamer Acquisition-Synthetic Data Analysis | |
Orji et al. | Imaging 3D sea surfaces from 3D dual-sensor towed streamer data | |
CN117388920A (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 |