CN109633744B - 地震子波的提取方法、装置、设备及存储介质 - Google Patents

地震子波的提取方法、装置、设备及存储介质 Download PDF

Info

Publication number
CN109633744B
CN109633744B CN201910049493.0A CN201910049493A CN109633744B CN 109633744 B CN109633744 B CN 109633744B CN 201910049493 A CN201910049493 A CN 201910049493A CN 109633744 B CN109633744 B CN 109633744B
Authority
CN
China
Prior art keywords
seismic
wavelets
record information
phase
extracting
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
CN201910049493.0A
Other languages
English (en)
Other versions
CN109633744A (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.)
Ocean University of China
CNOOC China Ltd Shanghai Branch
Original Assignee
Ocean University of China
China Oilfield Services Ltd Shanghai Branch
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 Ocean University of China, China Oilfield Services Ltd Shanghai Branch filed Critical Ocean University of China
Priority to CN201910049493.0A priority Critical patent/CN109633744B/zh
Publication of CN109633744A publication Critical patent/CN109633744A/zh
Application granted granted Critical
Publication of CN109633744B publication Critical patent/CN109633744B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

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

地震子波的提取方法、装置、设备及存储介质
技术领域
本发明实施例涉及地震勘探技术领域，尤其涉及一种地震子波的提取方法、装置、设备及存储介质。
背景技术
地震子波是一段具有确定的起始时间、能量有限且有一定延续长度的信号。在地震勘探中，地震子波可直接反映出地震波的传播特征，其在地震资料的处理和解释环节有着重要的作用，尤其在地震资料的处理中，地震子波的正确与否对处理结果的精度有着不可忽视的影响。在地震正演时，必须通过褶积结合地震子波来获取地震记录；在地震反演时，地震子波的准确性又会涉及波形匹配问题，进而影响到最终的反演精度。因此，如何从地震记录中提取准确的地震子波一直以来是地震勘探领域的研究热点。
地震子波提取方法总体上可分为确定性子波提取方法、统计性子波提取方法两大类。确定性子波提取方法是一种利用声波测井和密度测井资料的子波提取方法，其首先需要计算出反射系数序列，然后结合井旁地震道的数据通过反褶积技术求出地震子波。这种方法的优点是不需要对反射系数序列的分布做任何假设就可以得到较为准确的地震子波。然而传统的确定性子波提取方法在提取子波的过程中均需要利用测井资料，而在实际地震勘探中，测井资料往往难以获取，因此该类方法在实际地震子波提取中的应用受到了较大的限制。
Robinson于1975年提出了统计性子波提取方法，这类方法在不需要测井信息的情况下也可以得到地震子波的估计。统计性子波提取方法通常基于如下假设：地震子波是时不变的，地下反射率是具有白噪谱的随机序列，且地震道不含噪声，因此地震记录的自相关就可给定出地震子波自相关的一个估计，即可获取地震子波的振幅谱。而对于子波的相位谱，这类方法通常假设子波相位为最小相位、零相位或最大相位的一种。由于统计性子波提取方法不依赖于测井资料，其较好地解决了测井资料难以获取的问题，然而该类方法均是基于地震子波时不变的假设，即假设浅、深层均为同一子波，而实际上地震子波具有时变的特性，经过地层衰减后，其在浅层表现为高频，而在深层表现为低频，因此传统的统计性子波提取方法无法实现时变子波的提取，其无法反映真实的子波变化。
发明内容
本发明实施例提供一种地震子波的提取方法、装置、设备及存储介质，以实现对地震子波的提取，提高地震子波的准确性。
第一方面，本发明实施例提供了一种地震子波的提取方法，该方法包括：
获取地震记录信息；
根据所述地震记录信息确定多个同相轴；
根据所述多个同相轴进行子波提取，获得所述地震记录信息对应的子波。
进一步地，根据所述地震记录信息确定多个同相轴，包括：
根据所述地震记录信息获得同相加权的叠加速度谱；
获取所述叠加速度谱在不同时刻分别对应的速度极值；
根据所述不同时刻和分别对应的速度极值确定多个同相轴。
进一步地，根据所述地震记录信息获得同相加权的叠加速度谱按照如下公式计算：
Figure GDA0002811609150000031
其中，n为道号，1≤n≤N；xn为第n道的偏移距；τ表示位于第n道的传感器初次检测到地震波的时刻；b(v,τ)为同相加权因子，其计算公式为：
Figure GDA0002811609150000032
其中，λ表示阶数且λ≥2；L表示时窗的样点数；C为常数，取地震波平均振幅的0.01-0.001。
进一步地，获取所述叠加速度谱在不同时刻分别对应的速度极值，包括：
采用等值线追踪法获取所述叠加速度谱在不同时刻分别对应的速度极值。
进一步地，根据所述不同时刻和分别对应的速度极值确定多个同相轴按照如下公式计算：
Figure GDA0002811609150000033
其中，τn表示表示位于第n道的传感器初次检测到地震波的时刻，vn表示在τn时刻叠加速度谱的速度极值。
进一步地，根据所述多个同相轴进行子波提取，获得所述地震记录信息对应的子波，包括：
对所述多个同相轴分别进行傅里叶变换，获得多个振幅谱；
计算所述多个振幅谱的平均值，获得平均振幅谱；
多所述平均振幅谱进行傅里叶反变换，获得所述地震记录信息对应的子波。
第二方面，本发明实施例还提供了一种地震子波的提取装置，该装置包括：
地震记录信息获取模块，用于获取地震记录信息；
多个同相轴确定模块，用于根据所述地震记录信息确定多个同相轴；
子波提取模块，用于根据所述多个同相轴进行子波提取，获得所述地震记录信息对应的子波。
进一步地，所述多个同相轴确定模块，还用于：
根据所述地震记录信息获得同相加权的叠加速度谱；
获取所述叠加速度谱在不同时刻分别对应的速度极值；
根据所述不同时刻和分别对应的速度极值确定多个同相轴。
第三方面，本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如本发明实施例所述的地震子波的提取方法。
第四方面，本发明实施例提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如本发明实施例所述的地震子波的提取方法。
本发明实施例，首先获取地震记录信息，然后根据地震记录信息确定多个同相轴，最后根据多个同相轴进行子波提取，获得地震记录信息对应的子波。本发明实施例提供的地震子波的提取方法，根据地震记录信息对应的多个同相轴进行子波提取，可以提高提取的子波的准确性，从而为后续地震勘探的处理及反演解释工作奠定良好的基础。
附图说明
图1是本发明实施例一中一种地震子波的提取方法的流程图；
图2是本发明实施例一中的一种地震记录信息的示意图；
图3a是本发明实施例一中的浅层同相轴的振幅谱；
图3b是本发明实施例一中的中层同相轴的振幅谱；
图3c是本发明实施例一中的深层同相轴的振幅谱；
图4是本发明实施例二中的一种地震子波的提取装置的结构示意图；
图5是本发明实施例三中的一种计算机设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
实施例一
图1为本发明实施例一提供的一种地震子波的提取方法的流程图，本实施例可适用于对地震记录信息中的子波进行提取的情况，该方法可以由地震子波的提取装置来执行，该装置可由硬件和/或软件组成，并一般可集成在具有地震子波的提取功能的设备中，该设备可以是服务器、移动终端或服务器集群等电子设备。如图1所示，具体包括如下步骤：
步骤110，获取地震记录信息。
其中，地震波记录信息可以是根据埋在不同深度的传感器检测到地震强度组成的。地震记录信息由距离和时间两个变量表征，例如，可以利用d(x,t)表征地震记录信息，其中，x表示与震源之间的距离，t表示时间。示例性的，图2为本发明实施例中的一种地震记录信息的示意图，如图2所示，纵轴表示时间，横轴表示触感器的序号，每个传感器对应与震源的距离，地震记录信息表征的是在某一时刻距离震源一定长度的位置处的振动强度。
步骤120，根据地震记录信息确定多个同相轴。
其中，同相轴是地震记录上各道振动相位相同的极值(俗称波峰或波谷)的连线。可选的，根据地震记录信息确定多个同相轴，可通过下述方式实施：根据地震记录信息获得同相加权的叠加速度谱；获取叠加速度谱在不同时刻分别对应的速度极值；根据不同时刻和分别对应的速度极值确定多个同相轴。
其中，根据所述地震记录信息获得同相加权的叠加速度谱按照如下公式计算：
Figure GDA0002811609150000061
其中，n为道号，1≤n≤N；xn为第n道的偏移距；τ表示位于第n道的传感器初次检测到地震波的时刻；b(v,τ)为同相加权因子，其计算公式为：
Figure GDA0002811609150000062
其中，λ表示阶数且λ≥2；L表示时窗的样点数；C为常数，取地震波平均振幅的0.01-0.001。本实施例中，λ越大同相加权因子的分辨率越高。
获取叠加速度谱在不同时刻分别对应的速度极值的方式可以是：采用等值线追踪法获取叠加速度谱在不同时刻分别对应的速度极值。
其中，根据不同时刻和分别对应的速度极值确定多个同相轴按照如下公式计算：
Figure GDA0002811609150000063
其中，τn表示表示位于第n道的传感器初次检测到地震波的时刻，vn表示在τn时刻叠加速度谱的速度极值。
步骤130，根据多个同相轴进行子波提取，获得地震记录信息对应的子波。
具体的，根据多个同相轴进行子波提取，获得地震记录信息对应的子波，可通过下述方式实施：对多个同相轴分别进行傅里叶变换，获得多个振幅谱；计算多个振幅谱的平均值，获得平均振幅谱；多平均振幅谱进行傅里叶反变换，获得地震记录信息对应的子波。
示例性的，图3a-图3c为浅、中、深层三个同相轴对应的振幅谱。
本实施例的技术方案，首先获取地震记录信息，然后根据地震记录信息确定多个同相轴，最后根据多个同相轴进行子波提取，获得地震记录信息对应的子波。本发明实施例提供的地震子波的提取方法，根据地震记录信息对应的多个同相轴进行子波提取，可以提高提取的子波的准确性，从而为后续地震勘探的处理及反演解释工作奠定良好的基础。
实施例二
图4为本发明实施例二中的一种地震子波的提取装置的结构示意图。如图4所示，该装置包括：地震记录信息获取模块210，多个同相轴确定模块220和子波提取模块230。
地震记录信息获取模块210，用于获取地震记录信息；
多个同相轴确定模块220，用于根据地震记录信息确定多个同相轴；
子波提取模块230，用于根据多个同相轴进行子波提取，获得地震记录信息对应的子波。
可选的，多个同相轴确定模块220，还用于：
根据地震记录信息获得同相加权的叠加速度谱；
获取叠加速度谱在不同时刻分别对应的速度极值；
根据不同时刻和分别对应的速度极值确定多个同相轴。
可选的，根据地震记录信息获得同相加权的叠加速度谱按照如下公式计算：
Figure GDA0002811609150000081
其中，n为道号，1≤n≤N；xn为第n道的偏移距；τ表示位于第n道的传感器初次检测到地震波的时刻；b(v,τ)为同相加权因子，其计算公式为：
Figure GDA0002811609150000082
其中，λ表示阶数且λ≥2；L表示时窗的样点数；C为常数，取地震波平均振幅的0.01-0.001。
可选的，多个同相轴确定模块220，还用于：
采用等值线追踪法获取叠加速度谱在不同时刻分别对应的速度极值。
可选的，根据不同时刻和分别对应的速度极值确定多个同相轴按照如下公式计算：
Figure GDA0002811609150000083
其中，τn表示表示位于第n道的传感器初次检测到地震波的时刻，vn表示在τn时刻叠加速度谱的速度极值。
可选的，子波提取模块230，还用于：
对多个同相轴分别进行傅里叶变换，获得多个振幅谱；
计算多个振幅谱的平均值，获得平均振幅谱；
多平均振幅谱进行傅里叶反变换，获得地震记录信息对应的子波。
上述装置可执行本发明前述所有实施例所提供的方法，具备执行上述方法相应的功能模块和有益效果。未在本实施例中详尽描述的技术细节，可参见本发明前述所有实施例所提供的方法。
实施例三
图5为本发明实施例三提供的一种计算机设备的结构示意图。图5示出了适于用来实现本发明实施方式的计算机设备312的框图。图5显示的计算机设备312仅仅是一个示例，不应对本发明实施例的功能和使用范围带来任何限制。设备312典型的是承担地震子波的提取功能的计算设备。
如图5所示，计算机设备312以通用计算设备的形式表现。计算机设备312的组件可以包括但不限于：一个或者多个处理器316，存储装置328，连接不同系统组件(包括存储装置328和处理器316)的总线318。
总线318表示几类总线结构中的一种或多种，包括存储器总线或者存储器控制器，外围总线，图形加速端口，处理器或者使用多种总线结构中的任意总线结构的局域总线。举例来说，这些体系结构包括但不限于工业标准体系结构(Industry StandardArchitecture，ISA)总线，微通道体系结构(Micro Channel Architecture，MCA)总线，增强型ISA总线、视频电子标准协会(Video Electronics Standards Association，VESA)局域总线以及外围组件互连(Peripheral Component Interconnect，PCI)总线。
计算机设备312典型地包括多种计算机系统可读介质。这些介质可以是任何能够被计算机设备312访问的可用介质，包括易失性和非易失性介质，可移动的和不可移动的介质。
存储装置328可以包括易失性存储器形式的计算机系统可读介质，例如随机存取存储器(Random Access Memory，RAM)330和/或高速缓存存储器332。计算机设备312可以进一步包括其它可移动/不可移动的、易失性/非易失性计算机系统存储介质。仅作为举例，存储系统334可以用于读写不可移动的、非易失性磁介质(图5未显示，通常称为“硬盘驱动器”)。尽管图5中未示出，可以提供用于对可移动非易失性磁盘(例如“软盘”)读写的磁盘驱动器，以及对可移动非易失性光盘(例如只读光盘(Compact Disc-Read Only Memory，CD-ROM)、数字视盘(Digital Video Disc-Read Only Memory，DVD-ROM)或者其它光介质)读写的光盘驱动器。在这些情况下，每个驱动器可以通过一个或者多个数据介质接口与总线318相连。存储装置328可以包括至少一个程序产品，该程序产品具有一组(例如至少一个)程序模块，这些程序模块被配置以执行本发明各实施例的功能。
具有一组(至少一个)程序模块326的程序336，可以存储在例如存储装置328中，这样的程序模块326包括但不限于操作系统、一个或者多个应用程序、其它程序模块以及程序数据，这些示例中的每一个或某种组合中可能包括网络环境的实现。程序模块326通常执行本发明所描述的实施例中的功能和/或方法。
计算机设备312也可以与一个或多个外部设备314(例如键盘、指向设备、摄像头、显示器324等)通信，还可与一个或者多个使得用户能与该计算机设备312交互的设备通信，和/或与使得该计算机设备312能与一个或多个其它计算设备进行通信的任何设备(例如网卡，调制解调器等等)通信。这种通信可以通过输入/输出(I/O)接口322进行。并且，计算机设备312还可以通过网络适配器320与一个或者多个网络(例如局域网(Local AreaNetwork，LAN)，广域网Wide Area Network，WAN)和/或公共网络，例如因特网)通信。如图5所示，网络适配器320通过总线318与计算机设备312的其它模块通信。应当明白，尽管图中未示出，可以结合计算机设备312使用其它硬件和/或软件模块，包括但不限于：微代码、设备驱动器、冗余处理单元、外部磁盘驱动阵列、磁盘阵列(Redundant Arrays ofIndependent Disks，RAID)系统、磁带驱动器以及数据备份存储系统等。
处理器316通过运行存储在存储装置328中的程序，从而执行各种功能应用以及数据处理，例如实现本发明上述实施例所提供的地震子波的提取方法。
实施例四
本发明实施例四还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如本发明实施例所提供的地震子波的提取方法。
当然，本发明实施例所提供的一种计算机可读存储介质，其上存储的计算机程序不限于如上所述的方法操作，还可以执行本发明任意实施例所提供的地震子波的提取方法中的相关操作。
本发明实施例的计算机存储介质，可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括——但不限于无线、电线、光缆、RF等等，或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本发明操作的计算机程序代码，所述程序设计语言包括面向对象的程序设计语言—诸如Java、Smalltalk、C++，还包括常规的过程式程序设计语言—诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。

Claims (9)

1.一种地震子波的提取方法，其特征在于，包括：
获取地震记录信息；
根据所述地震记录信息确定多个同相轴；
根据所述多个同相轴进行子波提取，获得所述地震记录信息对应的子波；
其中，根据所述多个同相轴进行子波提取，获得所述地震记录信息对应的子波，包括：
对所述多个同相轴分别进行傅里叶变换，获得多个振幅谱；
计算所述多个振幅谱的平均值，获得平均振幅谱；
对所述平均振幅谱进行傅里叶反变换，获得所述地震记录信息对应的子波。
2.根据权利要求1所述的方法，其特征在于，根据所述地震记录信息确定多个同相轴，包括：
根据所述地震记录信息获得同相加权的叠加速度谱；
获取所述叠加速度谱在不同时刻分别对应的速度极值；
根据所述不同时刻和分别对应的速度极值确定多个同相轴。
3.根据权利要求2所述的方法，其特征在于，根据所述地震记录信息获得同相加权的叠加速度谱按照如下公式计算：
Figure FDA0003072907740000011
其中，n为道号，1≤n≤N；xn为第n道的偏移距；τ表示位于第n道的传感器初次检测到地震波的时刻，d(xn,t)为地震记录信息，E(v,τ)为同相加权的叠加速度谱，b(v,τ)为同相加权因子，其计算公式为：
Figure FDA0003072907740000012
其中，λ表示阶数且λ≥2；L表示时窗的样点数；C为常数，取地震波平均振幅的0.01-0.001。
4.根据权利要求3所述的方法，其特征在于，获取所述叠加速度谱在不同时刻分别对应的速度极值，包括：
采用等值线追踪法获取所述叠加速度谱在不同时刻分别对应的速度极值。
5.根据权利要求3所述的方法，其特征在于，根据所述不同时刻和分别对应的速度极值确定多个同相轴按照如下公式计算：
Figure FDA0003072907740000021
其中，τn表示表示位于第n道的传感器初次检测到地震波的时刻，vn表示在τn时刻叠加速度谱的速度极值。
6.一种地震子波的提取装置，其特征在于，包括：
地震记录信息获取模块，用于获取地震记录信息；
多个同相轴确定模块，用于根据所述地震记录信息确定多个同相轴；
子波提取模块，用于根据所述多个同相轴进行子波提取，获得所述地震记录信息对应的子波；
其中，子波提取模块，还用于：
对多个同相轴分别进行傅里叶变换，获得多个振幅谱；
计算多个振幅谱的平均值，获得平均振幅谱；
对平均振幅谱进行傅里叶反变换，获得地震记录信息对应的子波。
7.根据权利要求6所述的装置，其特征在于，所述多个同相轴确定模块，还用于：
根据所述地震记录信息获得同相加权的叠加速度谱；
获取所述叠加速度谱在不同时刻分别对应的速度极值；
根据所述不同时刻和分别对应的速度极值确定多个同相轴。
8.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1-5中任一所述的地震子波的提取方法。
9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1-5中任一所述的地震子波的提取方法。
CN201910049493.0A 2019-01-18 2019-01-18 地震子波的提取方法、装置、设备及存储介质 Active CN109633744B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910049493.0A CN109633744B (zh) 2019-01-18 2019-01-18 地震子波的提取方法、装置、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910049493.0A CN109633744B (zh) 2019-01-18 2019-01-18 地震子波的提取方法、装置、设备及存储介质

Publications (2)

Publication Number Publication Date
CN109633744A CN109633744A (zh) 2019-04-16
CN109633744B true CN109633744B (zh) 2021-07-13

Family

ID=66062164

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910049493.0A Active CN109633744B (zh) 2019-01-18 2019-01-18 地震子波的提取方法、装置、设备及存储介质

Country Status (1)

Country Link
CN (1) CN109633744B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110988986B (zh) * 2019-12-25 2021-01-01 成都理工大学 改善深层碳酸盐岩储层刻画精度的地震资料低频增强方法
CN117687095B (zh) * 2023-12-14 2024-08-27 中海石油(中国)有限公司上海分公司 一种基于方向分量数据提取与子波谱拓展的时变反褶积方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102928875B (zh) * 2012-11-07 2016-05-11 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于分数阶傅里叶域的子波提取方法
CN103675915B (zh) * 2013-12-20 2016-07-06 中国石油天然气股份有限公司 基于地震资料估计地层横向相对品质因子的方法和装置
US20150316666A1 (en) * 2014-05-05 2015-11-05 The Board Of Trustees Of The Leland Stanford Junior University Efficient Similarity Search of Seismic Waveforms
CN107132573B (zh) * 2017-02-28 2019-04-26 中国石油化工股份有限公司 一种应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的方法
CN109738951B (zh) * 2019-01-03 2019-12-20 国家深海基地管理中心 一种基于地震同相轴子波谱的时变反褶积方法

Also Published As

Publication number Publication date
CN109633744A (zh) 2019-04-16

Similar Documents

Publication Publication Date Title
Stork et al. Application of machine learning to microseismic event detection in distributed acoustic sensing data
Grigoli et al. Automated seismic event location by travel‐time stacking: An application to mining induced seismicity
CN108279436B (zh) 一种拾取地震波初至的方法及装置
US10295688B2 (en) Attenuating internal multiples from seismic data
CN106646611B (zh) 一种纵波和转换波地震数据时间匹配的方法及装置
CN107065013B (zh) 一种地震尺度下的层速度确定方法及装置
CN109633744B (zh) 地震子波的提取方法、装置、设备及存储介质
Zhou et al. An improved automatic picking method for arrival time of acoustic emission signals
CN110687592A (zh) 一种微地震震相识别初至拾取方法、设备和存储介质
CN108508481B (zh) 一种纵波转换波地震数据时间匹配的方法、装置及系统
CN113222029B (zh) 基于数据融合的地震成像方法、装置、设备及存储介质
CN117518233A (zh) 一种剩余时差校正的微地震定位方法及相关设备
CN114200522B (zh) 深度域地震子波提取方法和装置、存储介质及电子设备
CN115079257A (zh) 基于融合网络的q值估计及地震衰减补偿方法
CN110646846B (zh) Vti介质各向异性参数确定方法、装置和设备
CN113721296A (zh) 远震数据处理方法及装置
CN118091749B (zh) 基于动态时间规整的微地震事件初至拾取方法及系统
US20170336523A1 (en) Seismic signal processing method, apparatus and system
CN112649876A (zh) 建立地震偏移速度模型的方法和装置
CN113219526A (zh) 一种地震波初至自动拾取与筛选方法
CN115128675A (zh) 层间多次波的压制方法、装置、设备和存储介质
CN110824555B (zh) 地震能量均衡方法及装置、计算机可读存储介质
US11194068B1 (en) Systems and methods for object location detection such as detecting airplane crash location
CN110488348B (zh) 深度域地震数据的拓频处理方法、装置及存储介质
CN117849875B (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
TA01 Transfer of patent application right

Effective date of registration: 20190620

Address after: Room A901, Area A, 9th floor, 388 Tongxie Road, Changning District, Shanghai, 2003

Applicant after: Shanghai Branch of China National Offshore Oil Corporation

Applicant after: Ocean University of China

Address before: 200335 Zhonghai Oil Building, 388 Tongxie Road, Changning District, Shanghai (Shanghai)

Applicant before: Shanghai Branch of China National Offshore Oil Corporation

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant