CN112630822B - 一种ovt域叠前地震数据处理方法及系统 - Google Patents

一种ovt域叠前地震数据处理方法及系统 Download PDF

Info

Publication number
CN112630822B
CN112630822B CN201910950889.2A CN201910950889A CN112630822B CN 112630822 B CN112630822 B CN 112630822B CN 201910950889 A CN201910950889 A CN 201910950889A CN 112630822 B CN112630822 B CN 112630822B
Authority
CN
China
Prior art keywords
data
ovt
seismic data
domain
amplitude
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
CN201910950889.2A
Other languages
English (en)
Other versions
CN112630822A (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 CN201910950889.2A priority Critical patent/CN112630822B/zh
Publication of CN112630822A publication Critical patent/CN112630822A/zh
Application granted granted Critical
Publication of CN112630822B publication Critical patent/CN112630822B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/51Migration
    • G01V2210/512Pre-stack

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

本发明提出了一种OVT域叠前地震数据处理方法及系统。该方法包括:将原始地震数据根据不同的方位角与偏移距进行分选,得到OVT域道集;在OVT域对数据进行瞬时振幅和瞬时相位的分解;选取处理中心面元,将中心面元一定的空间范围内的道集基于叠加原理进行叠加;对所有中心面元重复以上处理步骤,得到所有处理中心面元的结果,对输出的相位数据与振幅数据进行重构。通过将原始数据分选为不同偏移距、不同方位角的OVT道集,利用OVT道集中信号相关性,在保相位的基础上增强被噪音掩盖的有效弱反射信号,整体提高叠前资料信噪比,使得叠前道集、叠加和成像剖面效果均得到明显改善。

Description

一种OVT域叠前地震数据处理方法及系统
技术领域
本发明属于地球物理勘探领域,涉及地震资料叠前预处理技术,更具体地,涉及一种OVT域叠前地震数据处理方法及系统。
背景技术
随着地震勘探技术的发展,探勘深度逐步加深。超深层复杂油气藏高精度地震波成像成为了新的勘探方向,其核心问题是实现深层目标地质体反射波的同相位叠加。宽频带的、无道间静校时差和道间子波特征一致的、规则无假频采样、无(层间)多次波、信噪比高的叠前数据体是实现超深层油藏目标高精度成像的基础。依据目标探区的数据观测方式、复杂近地表因素引起的噪音、反射同相轴的道间时差、反射子波特征的道间不一致、数据不规则、高频成分损失、目标层反射信号弱等情况,制定合适的数据预处理流程、QC方法以及合适的关键技术组合是目前处理此类问题的主要出发点。
由于勘探深度大,导致有效信号能量弱、信噪比低是预处理中的难点问题。目前,常规地震数据叠前去噪技术主要存在两方面的问题,一是在去除噪音的同时,会影响到有效弱信号;二是叠前道集上有效果,但叠加及成像剖面上无明显效果。
发明内容
针对深层复杂地下构造形成的弱反射和低信噪比问题,本发明提出了将地震数据分选到OVT(OVT全称为Offset Vector Tile,即矢量偏移距面元)域,由于OVT道集具有相同的偏移距和方位角信息,增强了数据间的相关性,在OVT域对数据进行瞬时振幅和瞬时相位的分解,在复数域实现弱信号的增强,将弱信号振幅增强后数据与相位增强后数据进行重构,在保数据频率、相位的基础上增强被噪音掩盖的有效弱反射信号,整体提高叠前资料信噪比。
根据本发明的一个方面,提供一种OVT域叠前地震数据处理方法,该方法包括:
将原始地震数据根据不同的方位角与偏移距进行分选,得到OVT域道集;
在OVT域对数据进行瞬时振幅和瞬时相位的分解;
选取处理中心面元,将中心面元一定的空间范围内的道集基于叠加原理进行叠加;
对所有中心面元重复以上处理步骤,得到所有处理中心面元的结果,对输出的相位数据与振幅数据进行重构。
具体地,利用地震数据玫瑰属性图,对原始数据进行不同偏移距与方位角的分组;
依据每组方位角与偏移距,对原始地震数据进行分选,得到多个OVT域道集。
进一步地,利用希尔伯特变换在复数域,实现OVT数据相位与振幅的分解。
优选地,希尔伯特数学表达式为:
其中xr为时间域地震数据,τ为辅助变量,xh为希尔伯特变换后地震数据。
优选地,瞬时振幅表达式为:
瞬时相位表达式为:
进一步地,对中心面元进行相关叠加,求取归一化后的叠加结果作为振幅或相位数据的处理结果。
进一步地,叠加后的相位数据与振幅数据分别为和/>则重构表达式为:
根据本发明的另一方面,提供一种OVT域叠前地震数据处理系统,该系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
将原始地震数据根据不同的方位角与偏移距进行分选,得到OVT域道集;
在OVT域对数据进行瞬时振幅和瞬时相位的分解;
选取处理中心面元,将中心面元一定的空间范围内的道集基于叠加原理进行叠加;
对所有中心面元重复以上处理步骤,得到所有处理中心面元的结果,对输出的相位数据与振幅数据进行重构。
进一步地,利用希尔伯特变换在复数域,实现OVT数据相位与振幅的分解,希尔伯特数学表达式为:
其中xr为时间域地震数据,τ为辅助变量,xh为希尔伯特变换后地震数据;
瞬时振幅表达式为:
瞬时相位表达式为:
进一步地,叠加后的相位数据与振幅数据分别为和/>则重构表达式为:
本发明提出了基于OVT域的深层弱信号增强方法,通过将原始数据分选为不同偏移距、不同方位角的OVT道集,由于OVT道集具有相同的偏移距和方位角信息,增强了数据间的相关性,在OVT域对数据进行瞬时振幅和瞬时相位的分解,在复数域实现弱信号的增强,将弱信号振幅增强后数据与相位增强后数据进行重构,在数据相对保真的基础上增强被噪音掩盖的有效弱反射信号,整体提高叠前资料信噪比,使得叠前道集、叠加和成像剖面效果均得到明显改善。
附图说明
通过结合附图对本公开示例性实施方式进行更详细的描述,本公开的上述以及其它目的、特征和优势将变得更加明显,其中,在本公开示例性实施方式中,相同的参考标号通常代表相同部件。
图1示出了根据本发明实施例的OVT域叠前地震数据处理处理方法流程图。
图2示出了根据本发明实施例的玫瑰图和偏移距分布图。
图3a-3c示出了根据本发明实施例的OVT域数据分选示意图。
图4示出了根据本发明实施例的处理前(左)后(右)CMP道集及速度谱。
图5示出了根据本发明实施例的处理前(左)后(右)叠加剖面及差剖面。
图6示出了根据本发明实施例的处理前(左)后(右)叠加剖面。
图7示出了根据本发明实施例的处理前(左)后(右)叠加剖面。
图8示出了现有商业软件增益处理后显示的差剖面。
图9示出了根据本发明实施例的处理前(左)后(右)信噪比属性剖面。
图10示出了根据本发明实施例的处理前(左)后(右)偏移剖面。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式,然而应该理解,可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本公开更加透彻和完整,并且能够将本公开的范围完整地传达给本领域的技术人员。
本发明属于地震资料叠前预处理技术领域。针对深层复杂地下构造形成的弱反射和低信噪比问题,本发明提出了基于OVT域的深层弱信号增强方法,通过将原始数据分选为不同偏移距、不同方位角的OVT道集,由于OVT道集具有相同的偏移距和方位角信息,增强了数据间的相关性,在OVT域对数据进行瞬时振幅和瞬时相位的分解,在复数域实现弱信号的增强,将弱信号振幅增强后数据与相位增强后数据进行重构,在数据相对保真的基础上增强被噪音掩盖的有效弱反射信号,整体提高叠前资料信噪比,使得叠前道集、叠加和成像剖面效果均得到明显改善。
基于OVT(Offset Vector Tile,OVT)域的数据处理为弱信号的增强提供了较好的实现途径,其原因是OVT域道集是共中心点域的一个子集,是一个考虑了方位角与偏移距的CMP单次覆盖剖面道集,OVT道集的分选方式是将玫瑰图划分为同心面元,每个面元的中心对应着一个方位角和一个偏移距,即为OVT道集,该道集中的相邻道,具有很强的相关性,有利于弱信号增强算法的实现。
按照本发明方法的原理,深层地震数据弱信号增强技术的实现过程为:将原始地震数据根据不同的方位角与偏移距进行分选,得到OVT域道集。由于OVT道集具有相同的偏移距和方位角信息,增强了数据间的相关性,在OVT域对数据进行瞬时振幅和瞬时相位的分解,对相位数据与振幅数据分别选取不同的处理面元作为中心面元,将中心面元一定的空间范围内道集基于叠加原理进行叠加,可以实现地震数据有效信号的增强。
首先,本发明提出了一种OVT域叠前地震数据处理方法,该方法包括:
将原始地震数据根据不同的方位角与偏移距进行分选,得到OVT域道集;
在OVT域对数据进行瞬时振幅和瞬时相位的分解;
选取处理中心面元,将中心面元一定的空间范围内的道集基于叠加原理进行叠加;
对所有中心面元重复以上处理步骤,得到所有处理中心面元的结果,对输出的相位数据与振幅数据进行重构。
如图1所示,根据本发明实施例的OVT域叠前地震数据处理方法,具体实现步骤为:
首先:输入CMP道集,然后将原始地震数据根据不同的方位角与偏移距进行分选,得到OVT域道集。
具体地,可以利用地震数据玫瑰属性图,对原始数据进行不同偏移距与方位角的分组;依据每组方位角与偏移距,对原始地震数据进行分选,得到多个OVT道集。
图2是一个实测三维探区的玫瑰图(上)和相应的偏移距分布图(下)。把玫瑰图的一个圆环投影到偏移距分布图,就得到对应的偏移距分布值。将玫瑰图划分为同心面元,每个面元的中心对应着一个方位角和一个偏移距,将每个面元投射到CDP域,则某些CDP会有一道或者多道,根据偏移距的差、方位角的差和偏离CDP中心位置的差等综合信息选出几何意义上的最优道,形成一个CDP域的单次剖面,即为分选出的OVT道集,该道集中的相邻道,具有很强的相关性。
图3是OVT域数据提取示意图。图3a)是一个OVT的原始数据,将其划分成网格数据(如图3b所示)。在多于一道的网格中,根据偏移距的差、方位角的差和偏离CDP中心位置的差,选出几何意义上的最优道,形成一个CDP域的单次剖面道集(图3c)。该道集即为具有相同偏移距和方位角的OVT道集,将此数据作为深层弱信号增强的处理数据,能有效提高信噪比的同时避免假同相轴的产生。
接下来,利用希尔伯特变换在复数域,实现OVT数据相位与振幅的分解。
具体地,在OVT道集域利用希尔伯特变换将数据变换到复数域,实现相位与振幅的分离。如图1所示,可分为余弦相位道集、正弦相位道集和瞬时振幅道集。
将数据在复数域进行分离,利用的希尔伯特数学表达式为:
其中xr为时间域地震数据,τ为辅助变量,xh为希尔伯特变换后地震数据。
则瞬时振幅表达式为:
瞬时相位表达式为:
余弦相位道集和正弦相位道集表达式的分子是不同的,但是两者求出的角度是同一个。因此可以选择余弦相位道集和正弦相位道集之一进行角度的求解。
接下来,选取处理中心面元,在数据的一定空间范围(如线点号3*3)内,对中心面元进行相关叠加,并求取归一化后的叠加结果作为振幅或相位数据的处理结果。具体地,将余弦相位道集和振幅道集视为同一个中心面元的两个不同特征表达的道集。两个道集的处理方式可以一样:例如在3*3数据体内,将数据体中所有CMP中相同位置的道集进行同向叠加并记录叠加次数,然后除以叠加次数作为中心面元每道的输出结果,最终得到该中心面元处理后的余弦相位道集和振幅道集。
最后,对所有中心面元重复以上处理步骤,得到所有处理中心面元的结果,对输出的相位数据与振幅数据进行重构,输出CMP道集,得到最终弱信号增强后结果。
针对每个面元数据在一定的局部空间范围内,对分离后的振幅数据与相位数据进行同相相关归一化叠加,最后对数据振幅和相位进行重构,得到弱信号增强后的地震数据。对相位数据与振幅数据分别选取不同的处理面元作为中心面元,由于同一个OVT道集中地震数据具有很强的相关性,因此将中心面元一定的空间范围内道集基于叠加原理可以实现地震数据有效信号的增强。弱信号增强处理后得到相位数据与振幅数据分别为则重构表达式为:
根据本发明的另一方面,还提供一种OVT域叠前地震数据处理系统,该系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
将原始地震数据根据不同的方位角与偏移距进行分选,得到OVT域道集;
在OVT域对数据进行瞬时振幅和瞬时相位的分解;
选取处理中心面元,将中心面元一定的空间范围内的道集基于叠加原理进行叠加;
对所有中心面元重复以上处理步骤,得到所有处理中心面元的结果,对输出的相位数据与振幅数据进行重构。
具体地,利用希尔伯特变换在复数域,实现OVT数据相位与振幅的分解,希尔伯特数学表达式为:
其中xr为时间域地震数据,τ为辅助变量,xh为希尔伯特变换后地震数据;
瞬时振幅表达式为:
瞬时相位表达式为:
叠加后的相位数据与振幅数据分别为和/>则重构表达式为:
本发明面向复杂地下构造形成的弱反射和低信噪比问题,提出了基于OVT域的深层弱信号增强方法,通过将原始数据分选为不同偏移距、不同方位角的OVT道集,利用OVT道集中信号相关性,在保相位的基础上增强被噪音掩盖的有效弱反射信号,整体提高叠前资料信噪比,使得叠前道集、叠加和成像剖面效果均得到明显改善。
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
为了验证算法的有效性及对实际数据的处理效果,进行了实际数据的应用处理。
如图4左为原始处理前CMP数据及对应的速度谱,从其速度谱可知在浅层速度谱能量聚焦性较好。在CMP道集上其深层由于信息能量弱,信噪比较低,其速度谱在深层速度谱能量聚焦性较差,这不利用速度的拾取和后续处理的进行;图4右为经过本发明处理后的CMP数据及对应的速度谱,速度谱在深层位置能量聚焦性明显变好,CMP道集中数据同相轴连续性和信噪比都得到提升,弱的双曲同相轴显露出来。
图5左为工区所有原始CMP数据叠加后得到的叠加剖面,中间为本发明处理后得到的新的叠加剖面,右边为处理前后两个剖面的差剖面,对比处理前后的叠加剖面可知,剖面在中深层信噪比和同相轴连续性明显变好并且差剖面中不含有明显有效信号和同相轴,表明本发明处理后没有损伤有效信号也没有引入假的同相轴信息,说明本发明具有良好的相对保真性。
图6为叠加剖面加了相同增益后用主流商业软件显示的结果,加了相同增益后数据能量级别一致,该显示更能体现处理前后信噪比差别,本发明处理后的叠加剖面信噪比提高明显,特别是在中深层,图7也可以得出类似结论。图8为商业软件显示的差剖面,同时也做了增益处理,差剖面中无有效信号,反映的基本为噪音,体现了该发明的先进性。
为了进一步说明处理前后的效果,对处理前后叠加面进行了信噪比的定量分析,信噪比属性图9说明本发明处理后整体信噪比提高,被噪音淹没的弱有效信号得到增强,同相轴更为连续性。图10为偏移对比剖面,作为一种叠前处理方法,对叠前道集的处理在偏移成像剖面上也得到了良好的应用效果,偏移剖面从浅至深信噪比都得到不同程度的提高,中深层目标区同相轴更为清晰连续。
实际数据的应用结果表明,本发明无论是在叠前道集、速度谱、叠加剖面、信噪比属性、偏移剖面上均取得了较好的应用效果,克服了常规方法只对叠前有效,偏移后效果不明显的缺点,通过多种途径的质量监控也说明了该发明具有良好的相对保真性,这对于保真处理以及解释工作的开展具有重要意义,保护了地震数据中的微弱有效信号不受损失,有效改善了偏移成像的质量。
以上已经描述了本公开的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。

Claims (8)

1.一种OVT域叠前地震数据处理方法,其特征在于,该方法包括:
将原始地震数据根据不同的方位角与偏移距进行分选,得到OVT域道集;
在OVT域对数据进行瞬时振幅和瞬时相位的分解;
选取处理中心面元,将中心面元一定的空间范围内的道集基于叠加原理进行叠加;
对所有中心面元重复以上处理步骤,得到所有处理中心面元的结果,对输出的相位数据与振幅数据进行重构;
其中,利用希尔伯特变换在复数域,实现OVT数据相位与振幅的分解;
其中,希尔伯特数学表达式为:
其中xr为时间域地震数据,τ为辅助变量,xh为希尔伯特变换后地震数据。
2.根据权利要求1所述的OVT域叠前地震数据处理方法,其特征在于,利用地震数据玫瑰属性图,对原始数据进行不同偏移距与方位角的分组;
依据每组方位角与偏移距,对原始地震数据进行分选,得到多个OVT域道集。
3.根据权利要求1所述的OVT域叠前地震数据处理方法,其特征在于,瞬时振幅表达式为:
瞬时相位表达式为:
4.根据权利要求3所述的OVT域叠前地震数据处理方法,其特征在于,对中心面元进行相关叠加,求取归一化后的叠加结果作为振幅或相位数据的处理结果。
5.根据权利要求4所述的OVT域叠前地震数据处理方法,其特征在于,
叠加后的相位数据与振幅数据分别为和/>则重构表达式为:
6.一种OVT域叠前地震数据处理系统,其特征在于,该系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
将原始地震数据根据不同的方位角与偏移距进行分选,得到OVT域道集;
在OVT域对数据进行瞬时振幅和瞬时相位的分解;
选取处理中心面元,将中心面元一定的空间范围内的道集基于叠加原理进行叠加;
对所有中心面元重复以上处理步骤,得到所有处理中心面元的结果,对输出的相位数据与振幅数据进行重构;
其中,利用希尔伯特变换在复数域,实现OVT数据相位与振幅的分解;
其中,希尔伯特数学表达式为:
其中xr为时间域地震数据,τ为辅助变量,xh为希尔伯特变换后地震数据。
7.根据权利要求6所述的OVT域叠前地震数据处理系统,其特征在于,
瞬时振幅表达式为:
瞬时相位表达式为:
8.根据权利要求7所述的OVT域叠前地震数据处理系统,其特征在于,
叠加后的相位数据与振幅数据分别为和/>则重构表达式为:
CN201910950889.2A 2019-10-08 2019-10-08 一种ovt域叠前地震数据处理方法及系统 Active CN112630822B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910950889.2A CN112630822B (zh) 2019-10-08 2019-10-08 一种ovt域叠前地震数据处理方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910950889.2A CN112630822B (zh) 2019-10-08 2019-10-08 一种ovt域叠前地震数据处理方法及系统

Publications (2)

Publication Number Publication Date
CN112630822A CN112630822A (zh) 2021-04-09
CN112630822B true CN112630822B (zh) 2024-06-18

Family

ID=75283427

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910950889.2A Active CN112630822B (zh) 2019-10-08 2019-10-08 一种ovt域叠前地震数据处理方法及系统

Country Status (1)

Country Link
CN (1) CN112630822B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115184994B (zh) * 2022-06-08 2023-04-18 北京东方联创地球物理技术有限公司 一种多方位角采集的三维地震数据融合处理方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105277985A (zh) * 2014-07-24 2016-01-27 中国石油化工股份有限公司 一种基于图像处理的ovt域地震数据规则化方法
CN105629306A (zh) * 2014-10-27 2016-06-01 中国石油化工股份有限公司 一种信噪比模型建立方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6446007B1 (en) * 1998-11-05 2002-09-03 Exxonmobil Upstream Research Company Method for controlled-amplitude prestack time migration of seismic data
AU2013401761B2 (en) * 2013-09-25 2020-10-15 China Petroleum & Chemical Corporation Method and Device for Increasing Frequency of Seismic Digital Signal
CN104459772B (zh) * 2013-09-25 2017-08-18 中国石油化工股份有限公司 一种地震数字信号的提频方法及装置
CN105044769B (zh) * 2015-06-10 2017-11-10 中国石油集团川庆钻探工程有限公司地球物理勘探公司 提高地震信号的分辨率的方法
CN106054253A (zh) * 2016-08-01 2016-10-26 中国石油天然气集团公司 一种道集同相轴校正的方法及装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105277985A (zh) * 2014-07-24 2016-01-27 中国石油化工股份有限公司 一种基于图像处理的ovt域地震数据规则化方法
CN105629306A (zh) * 2014-10-27 2016-06-01 中国石油化工股份有限公司 一种信噪比模型建立方法

Also Published As

Publication number Publication date
CN112630822A (zh) 2021-04-09

Similar Documents

Publication Publication Date Title
CN109143363B (zh) 海洋拖缆双检采集鬼波压制方法及系统
CN102073064B (zh) 一种利用相位信息提高速度谱分辨率的方法
CN105510975B (zh) 提高地震数据信噪比的方法及装置
CN116520419B (zh) 一种热流体裂缝通道识别方法
WO2021127382A1 (en) Full waveform inversion in the midpoint-offset domain
CN112630822B (zh) 一种ovt域叠前地震数据处理方法及系统
Chu et al. Inversion of source parameters for moderate earthquakes using short-period teleseismic P waves
Li et al. Surface‐wave tomography of eastern and central Tibet from two‐plane‐wave inversion: Rayleigh‐wave and Love‐wave phase velocity maps
CN103048641A (zh) 具有联合约束优化形式的矢量阵稳健聚焦处理方法
CN110780346A (zh) 一种隧道超前探测复杂地震波场的分离方法
Wang et al. Iterative deblending using unsupervised learning with double-deep neural networks
CN115993641A (zh) 一种提取被动源面波频散曲线的方法
CN114545495A (zh) 一种ps波地震道数据的ovt道集处理方法、装置及电子设备
CN114076985B (zh) 可控震源的噪音压制方法、装置、设备和存储介质
CN112083490B (zh) 地震数据噪声衰减方法及装置
Luo et al. Slice-Relation-Clustering Framework via Horizontal Angle Information for 3-D Tree Roots Reconstruction
CN112415592B (zh) 地震资料五维谱解析噪声压制方法、存储介质及计算设备
CN113391348B (zh) 用于叠前反演的共反射点道集构建方法及装置
Shi et al. 3D high-resolution radon transform based on a strong sparse Lp–1 norm and its applications
CN117826248B (zh) 基于多尺度观测背景噪声聚束的面波频散提取方法
Jia et al. A bipolar-bisection piecewise encoding scheme for multi-source reverse time migration
JP7295979B2 (ja) 海底に埋在する1つ又は複数の物体を検出するための方法、装置、及びコンピュータプログラム
CN111781645B (zh) 一种利用海底地层和初至波联合求取海洋地震子波方法
Liu et al. Research on a matching detection method for magnetic anomaly of underwater targets
Li et al. N-th root slant stack for enhancing weak seismic signals

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