CN109143363B - 海洋拖缆双检采集鬼波压制方法及系统 - Google Patents

海洋拖缆双检采集鬼波压制方法及系统 Download PDF

Info

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
Application number
CN201710467201.6A
Other languages
English (en)
Other versions
CN109143363A (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.)
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
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 China Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201710467201.6A priority Critical patent/CN109143363B/zh
Publication of CN109143363A publication Critical patent/CN109143363A/zh
Application granted granted Critical
Publication of CN109143363B publication Critical patent/CN109143363B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/30Noise handling
    • G01V2210/32Noise 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包括海底一次反射信号及海平面反射鬼波信号:
Figure GDA0002941801690000021
假设P0(ω,kx,0)为海面处接收到的压力波场,则无鬼波上行波场信号Pup(ω,kx,z)可以表示为:
Figure GDA0002941801690000022
建立无鬼波上行波场信号与压力波场信号的关系为:
Figure GDA0002941801690000023
其中kz为垂直波数,Z为深度。
进一步地,根据牛顿第二定律,压力波场信号p和垂直速度波场信号vz之间的关系为:
Figure GDA0002941801690000024
在频率波数域进行计算:
Figure GDA0002941801690000025
即:
Figure GDA0002941801690000026
其中,Vz(ω,kx,z)表示频率波数域垂直速度波场信号,ρ表示密度,ω表示频率,k表示波数,z表示拖缆深度。进一步地,进一步包括建立无鬼波信号与压力波场信号和垂直速度波场信号的关系,由方程(3)和(6),可得:
Figure GDA0002941801690000031
由以上可知:
Figure GDA0002941801690000032
进一步地,基于最小平方算法,得到:
Figure GDA0002941801690000033
海洋拖缆双检采集技术可以直接获得压力波场P(ω,kx,z)和垂直速度波场Vz(ω,kx,z),因此通过上式可以直接获得宽频无鬼波上行波场Pup(ω,kx,z)。
经过双检采集压力波场和垂直速度波场合并,即可消除鬼波,获得宽频地震数据。
根据本发明另一方面,提供一种海洋拖缆双检采集鬼波压制系统,包括:
存储器,存储有可执行指令;
处理器,运行所述存储器存储的所述可执行指令,执行以下步骤:
获取压力波场信号和垂直速度波场信号;
建立鬼波无鬼波上行波场信号与压力波场信号的关系;
根据牛顿第二定律,建立压力波场信号和垂直速度波场信号之间的关系;
基于最小平方算法,合并压力波场和垂直速度波场,消除鬼波信号。根据本发明的又一个方面,提供一种记录介质,所述记录介质中存储有计算机可执行指令;所述计算机可执行指令运行时,执行以下步骤:
获取压力波场信号和垂直速度波场信号;
建立鬼波无鬼波上行波场信号与压力波场信号的关系;
根据牛顿第二定律,建立压力波场信号和垂直速度波场信号之间的关系;
基于最小平方算法,合并压力波场和垂直速度波场,消除鬼波信号。
本发明采用基于最小平方的双检数据合并方法,可实现水、陆检数据合并,不依赖地下假设信息,稳定性好,受采集环境等制约少,能够适应复杂海洋地形和地质情况。
附图说明
通过结合附图对本公开示例性实施方式进行更详细的描述,本公开的上述以及其它目的、特征和优势将变得更加明显,其中,在本公开示例性实施方式中,相同的参考标号通常代表相同部件。
图1示出了海洋双检采集技术水检信号叠加剖面。
图2示出了海洋双检采集技术叠加剖面。
图3示出了频谱对比。
图4示出了本发明的海洋拖缆双检采集鬼波压制方法流程图。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式,然而应该理解,可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本公开更加透彻和完整,并且能够将本公开的范围完整地传达给本领域的技术人员。
海洋拖缆双检采集中,鬼波对一次波的波形及频带宽度产生极大的影响,鬼波压制是提高海上地震资料分辨率、保真度及宽频采集效果的重要因素。
海洋拖缆双检采集技术的核心部分为宽频原始数据的采集和水、陆检数据的合并处理,双检数据的合并方法选取不当并不能实现鬼波的有效压制,低频缺失、陷波现象等问题并无法得到有效消除。因此研究最优的水陆检数据合并算法可以发挥拖缆双检宽频采集技术的最佳效果,从而提高低频能量,压制陷波,拓宽原始资料频带。
对于海洋拖缆双检采集技术数据处理,本发明采用基于最小平方的双检数据合并方法,能够适应复杂海洋地形和地质情况。实际资料处理结果表明,该技术在很好的去除鬼波的同时极大的拓宽了地震资料的频带,尤其是低频端,有利于后续资料的处理解释,且能有效压制随机噪声,提高信噪比。
如图4所示,本公开提出了一种海洋拖缆双检采集鬼波压制方法,包括:
获取压力波场信号和垂直速度波场信号;
建立鬼波无鬼波上行波场信号与压力波场信号的关系;
根据牛顿第二定律,建立压力波场信号和垂直速度波场信号之间的关系;
基于最小平方算法,合并压力波场和垂直速度波场,消除鬼波信号。
对于海洋拖缆采集,检波器信号由两部分组成:海底一次反射及海平面反射鬼波:
Figure GDA0002941801690000051
假设P0(ω,kx,0)为海面处接收到的压力波场,根据波场延拓原理,则无鬼波上行波场信号Pup(ω,kx,z)可以表示为:
Figure GDA0002941801690000052
即:
Figure GDA0002941801690000053
根据牛顿第二定律,压力波场p和垂直速度波场vz之间的关系为:
Figure GDA0002941801690000054
在频率波数域进行计算:
Figure GDA0002941801690000055
即:
Figure GDA0002941801690000056
其中,Vz(ω,kx,z)表示频率波数域垂直速度波场信号,ρ表示密度,ω表示频率,k表示波数,z表示拖缆深度。由方程3和6,可得:
Figure GDA0002941801690000061
由以上可知:
Figure GDA0002941801690000062
基于最小平方算法,可得:
Figure GDA0002941801690000063
海洋拖缆双检采集技术可以直接获得压力波场P(ω,kx,z)和垂直速度波场Vz(ω,kx,z),因此通过上式可以直接获得宽频无鬼波上行波场Pup(ω,kx,z)。
经过双检采集压力波场和垂直速度波场合并,即可消除鬼波,获得宽频地震数据。
根据本发明另一方面,提供一种基于最小平方算法的海洋拖缆双检采集鬼波压制系统,包括:
存储器,存储有可执行指令;
处理器,运行所述存储器存储的所述可执行指令,执行以下步骤:
获取压力波场信号和垂直速度波场信号;
建立鬼波无鬼波上行波场信号与压力波场信号的关系;
根据牛顿第二定律,建立压力波场信号和垂直速度波场信号之间的关系;
基于最小平方算法,合并压力波场和垂直速度波场,消除鬼波信号。
本发明采用基于最小平方的双检数据合并方法,可实现水、陆检数据合并,不依赖地下假设信息,稳定性好,受采集环境等制约少,能够适应复杂海洋地形和地质情况。
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
为了验证基于最小平方拖缆双检鬼波压制方法的实用效果,我们选取海上双检地震资料进行处理和分析,数据为海上水检、陆检地震记录,震源深度15m,拖缆深度25m。
基于最小平方算法进行水检、陆检数据合并,图1为海洋双检采集技术水检信号叠加剖面,图2为海洋双检采集技术叠加剖面,可见,经过鬼波压制,地震记录分辨率明显提高,由鬼波产生的虚假同相轴被很好的去除。
通过频谱分析,如图3,同样可以看出,水检记录在30hz频率附近存在明显的陷波,经过水陆检合并,陷波处能量得到有效提升,且低频能量有所提高,地震信号频带得到有效拓宽。
本发明采用基于最小平方的双检数据合并方法,能够适应复杂海洋地形和地质情况,实际资料处理结果表明,基于最小平方的海洋拖缆双检采集鬼波压制技术在很好的去除鬼波的同时极大的拓宽了地震资料的频带,有利于后续资料的处理解释。
以上已经描述了本公开的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。

Claims (6)

1.一种海洋拖缆双检采集鬼波压制方法,其特征在于,包括:
获取压力波场信号和垂直速度波场信号;
建立无鬼波上行波场信号与压力波场信号的关系;
根据牛顿第二定律,建立压力波场信号和垂直速度波场信号之间的关系;
基于最小平方算法,合并压力波场和垂直速度波场,消除鬼波信号;
其中,根据牛顿第二定律,压力波场信号p和垂直速度波场信号vz之间的关系为:
Figure FDA0002941801680000011
在频率波数域进行计算:
Figure FDA0002941801680000012
即:
Figure FDA0002941801680000013
其中,Vz(ω,kx,z)表示频率波数域垂直速度波场信号,ρ表示密度,ω表示频率,k表示波数,z表示拖缆深度,P0(ω,kx,0)为海面处接收到的压力波场。
2.根据权利要求1所述的海洋拖缆双检采集鬼波压制方法,其特征在于,压力波场信号P包括海底一次反射信号及海平面反射鬼波信号:
Figure FDA0002941801680000014
假设P0(ω,kx,0)为海面处接收到的压力波场,则无鬼波上行波场信号Pup(ω,kx,z)表示为:
Figure FDA0002941801680000015
建立无鬼波上行波场信号与压力波场信号的关系为:
Figure FDA0002941801680000016
其中kz为垂直波数,z为深度。
3.根据权利要求2所述的海洋拖缆双检采集鬼波压制方法,其特征在于,进一步包括建立无鬼波信号与压力波场信号和垂直速度波场信号的关系,由方程(3)和(6)得到:
Figure FDA0002941801680000021
由以上可知:
Figure FDA0002941801680000022
4.一种海洋拖缆双检采集鬼波压制系统,其特征在于,包括:
存储器,存储有可执行指令;
处理器,运行所述存储器存储的所述可执行指令,执行以下步骤:
获取压力波场信号和垂直速度波场信号;
建立鬼波无鬼波上行波场信号与压力波场信号的关系;
根据牛顿第二定律,建立压力波场信号和垂直速度波场信号之间的关系;
基于最小平方算法,合并压力波场和垂直速度波场,消除鬼波信号;
其中,根据牛顿第二定律,压力波场信号p和垂直速度波场信号vz之间的关系为:
Figure FDA0002941801680000023
在频率波数域进行计算:
Figure FDA0002941801680000024
即:
Figure FDA0002941801680000025
其中,Vz(ω,kx,z)表示频率波数域垂直速度波场信号,ρ表示密度,ω表示频率,k表示波数,z表示拖缆深度,P0(ω,kx,0)为海面处接收到的压力波场。
5.根据权利要求4所述的海洋拖缆双检采集鬼波压制系统,其特征在于,压力波场信号P包括海底一次反射信号及海平面反射鬼波信号:
Figure FDA0002941801680000031
假设P0(ω,kx,0)为海面处接收到的压力波场,则无鬼波上行波场信号Pup(ω,kx,z)表示为:
Figure FDA0002941801680000032
建立无鬼波上行波场信号与压力波场信号的关系为:
Figure FDA0002941801680000033
其中kz为垂直波数,z为深度。
6.根据权利要求5所述的海洋拖缆双检采集鬼波压制系统,其特征在于,进一步包括建立无鬼波信号与压力波场信号和垂直速度波场信号的关系,由方程(3)和(6)得到:
Figure FDA0002941801680000034
由以上得到:
Figure FDA0002941801680000035
CN201710467201.6A 2017-06-19 2017-06-19 海洋拖缆双检采集鬼波压制方法及系统 Active CN109143363B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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