CN112180445A - 地震勘探资料噪音去噪方法、存储介质和计算机系统 - Google Patents

地震勘探资料噪音去噪方法、存储介质和计算机系统 Download PDF

Info

Publication number
CN112180445A
CN112180445A CN201910596349.9A CN201910596349A CN112180445A CN 112180445 A CN112180445 A CN 112180445A CN 201910596349 A CN201910596349 A CN 201910596349A CN 112180445 A CN112180445 A CN 112180445A
Authority
CN
China
Prior art keywords
data
seismic
denoising
noise
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.)
Withdrawn
Application number
CN201910596349.9A
Other languages
English (en)
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 CN201910596349.9A priority Critical patent/CN112180445A/zh
Publication of CN112180445A publication Critical patent/CN112180445A/zh
Withdrawn legal-status Critical Current

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
    • 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/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

地震勘探资料噪音去噪方法、存储介质和计算机系统
技术领域
本发明涉及地球物理勘探技术领域，尤其涉及一种基于偏移距校正的地震勘探资料噪音去噪方法，以及用于存储能够实现该方法的计算机程序的存储介质和用于实施该方法的计算机系统。
背景技术
在地震勘探三维(或宽线)资料处理中往往存在相干噪音十分严重的情况。随着技术的不断发展，地震勘探资料处理对保幅的要求越来越高，而保幅处理的效果很大程度依赖于资料噪音压制的质量。相干噪音往往与有效波混叠在一起，尤其是在三维(或宽线)的非纵测线观测情况下，在共检波点道集(或共炮点道集)记录中，线性噪音的视速度会随炮检点相对位置关系变化而变化。
如图1所示，炮点为S，相邻两个检波点分别为Ri，Ri+1，如果炮点和检波点之间有非纵偏移距dy，第i道的纵向偏移距为dxi，则任意道德炮检偏移距Fi为：
Figure BDA0002117786450000011
随着i的变化，该方程为典型的双曲线方程，相邻两道之间的道间距离
Figure BDA0002117786450000012
它不是常数。这种情况会导致在地震资料中出现以下两个问题：
①相干噪音呈现双曲线化(如图4所示)；
②道间距不是常数。
对于后一问题，道间距的不一致会给傅里叶变换引入误差，此外，相干噪音的双曲线化会进一步降低目前各种基于固定视速度的线性噪音压制技术的处理效果。
发明内容
针对上述目前三维(或宽线)地震资料处理中普遍存在的线性噪音压制效果不理想的技术问题，本发明提供了一种基于偏移距校正的地震勘探资料线性噪音压制方法，以提高地震勘探资料的精度，从而为进行解释分析岩石物理特征分布规律奠定基础。
具体地，本发明提供了一种地震勘探资料噪音去噪方法，包括以下步骤：
将探区的地震勘探资料的地震数据分选成共检波点道集；
从共检波点道集的各道地震数据的道头信息中提取各道的炮检偏移距和炮检方位角，根据各道的炮检偏移距和炮检方位角以及探区的线性噪音的视速度，确定各道的偏移距校正量；
基于各道的偏移距校正量对共检波点道集的各道地震数据进行偏移距校正，使得其中的非线性相干噪音变换为线性相干噪音；
对经过校正的共检波点道集进行线性噪音去噪处理；
基于各道的偏移距校正量对经过去噪处理的共检波点道集的各道地震数据进行偏移距逆校正，以变换回原道集。
根据本发明的实施例，在上述方法中，在将探区的地震勘探资料的地震数据分选成共检波点道集之前，先对探区的地震勘探资料进行预处理。
根据本发明的实施例，在上述方法中，所述预处理包括对地震勘探资料赋上空间坐标信息，以及对地震勘探资料进行静校正和振幅恢复处理。
根据本发明的实施例，在上述方法中，具体地，根据观测系统的炮检坐标，将探区的地震勘探资料的地震数据分选成共检波点道集。
根据本发明的实施例，在上述方法中，所述探区的线性噪音包括探区的折射波和/或探区的面波。
根据本发明的实施例，在上述方法中，所述偏移距校正量为时差校正量。
根据本发明的实施例，在上述方法中，所述时差校正量通过下式计算：
Figure BDA0002117786450000021
式中，v为探区的线性噪音的视速度，Fi为第i道的炮检偏移距和θi为第i道的炮检方位角，Δti为第i道的时差校正量。
此外，本发明还提供了一种存储介质，其中存储有用于实现上述地震勘探资料噪音去噪方法的程序。
另外，本发明还提供了一种计算机系统，包括处理器和存储介质，其中所述处理器用于执行存储在所述存储介质中的用于实现上述地震勘探资料噪音去噪方法的程序。
与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：
本发明依据勘探工区的相干噪音表现最强的数据域的炮检坐标，对地震记录进行道集分选，形成共检波点道集，采用偏移距校正的手段，将其中非线性相干噪音变换成线性相干噪音，以使地震记录符合线性噪音压制方法的假设条件，提高去噪处理的质量，改善资料成像和保真度，达到保幅去噪、不损失有效波成分的处理目的。采用该发明提供的去噪方法可为获得高品质的处理成果以及后续的储层地质分析奠定良好的基础，具有很好的应用前景。
附图说明
通过结合附图阅读下文示例性实施例的详细描述可更好地理解本公开的范围。其中所包括的附图是：
图1是本发明实的基于偏移距校正的线性噪音压制方法原理的示意图；
图2是本发明实施例一的基于偏移距校正的线性噪音压制处理的流程图；
图3是本发明实施例二的基于偏移距校正的线性噪音压制处理的流程图；
图4是本发明实施例二的进行偏移距校正前的共检波点道集记录的示意图；
图5是本发明实施例二的进行偏移距校正后的共检波点道集记录的示意图；
图6是本发明实施例二的进行偏移距校正后的共检波点道集在完成线性噪音压制后的示意图；
图7是本发明实施例二的在进行线性噪音压制后完成偏移距逆校正的共检波点道集记录的示意图；
图8是本发明实施例三的未应用基于偏移距校正的线性噪音压制处理方法的叠加剖面的示意图；
图9为本发明实施例三的已应用基于偏移距校正的线性噪音压制处理方法的叠加剖面的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，以下将结合附图及实施例来详细说明本发明的实施方法，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。
实施例一
本发明实施例一主体构思是采用如图2所示的方法步骤对地震勘探资料进行噪音压制处理，以获得高品质的地震资料。
如图2所示，该方法主要包括以下几个步骤：
S110，将探区的地震勘探资料的地震数据分选成共检波点道集；
S120，从共检波点道集的各道地震数据的道头信息中提取各道的炮检偏移距和炮检方位角，根据各道的炮检偏移距和炮检方位角以及探区的线性噪音的视速度，确定各道的偏移距校正量；
S130，基于各道的偏移距校正量对共检波点道集的各道地震数据进行偏移距校正，使得其中的非线性相干噪音变换为线性相干噪音；
S140，对经过校正的共检波点道集进行线性噪音去噪处理；
S150，基于各道的偏移距校正量对经过去噪处理的共检波点道集的各道地震数据进行偏移距逆校正，以变换回原道集。
其中，步骤S120中所提及的偏移距校正量为时差校正量，也称为正常时差校正量。
具体地，所述时差校正量可以通过下式确定：
Figure BDA0002117786450000041
式中，v为探区的线性噪音的视速度，Fi为第i道的炮检偏移距和θi为第i道的炮检方位角，Δti为第i道的时差校正量。
实施例二
实际上，在实际应用时不能直接对探区的地震勘探资料实施上述方法。这是因为现场获得的地震勘探资料还需要进行进一步的处理，在此也称为预处理。
具体地，如图3至图7所示，在本实施例二中，本发明提供的地震勘探资料噪音去噪方法包括如下步骤：
S210，对探区的地震勘探资料进行预处理，所述预处理包括对地震勘探资料加载空间属性(即赋于空间坐标)，以及对地震勘探资料进行静校正和振幅恢复处理。
S220，根据观测系统的炮检坐标，将经过预处理的探区的地震勘探资料的地震数据分选成共检波点道集。
S230，从共检波点道集的各道地震数据的道头信息中提取各道的炮检偏移距和炮检方位角，根据各道的炮检偏移距和炮检方位角以及探区的线性噪音的视速度，确定各道的偏移距校正量。其中，因为偏移距是因炮检距所引起的正常时差，因此采用通过下式获得的正常时差校正量作为偏移距校正量：
Figure BDA0002117786450000051
式中，v为探区的线性噪音的视速度，Fi为第i道的炮检偏移距和θi为第i道的炮检方位角，Δti为第i道的正常时差校正量。
S240，基于各道的偏移距校正量对共检波点道集的各道地震数据进行偏移距校正，使得其中的非线性相干噪音变换为线性相干噪音。
具体地，将图1中Ri投影到Ri'，原来非纵偏移距dy的排列ly,变换到零纵测线l0。通过这一对非纵测线的排列进行纵测线化的变换过程,实现了非线性相干噪音的线性化，同时也将非常量的偏移距增量(道间距)变换成常量增量(道距)。经过这样的变换，实现了非线性相干噪音的线性化，使得地震资料的数据能够严格满足基于固定速度的线性噪音压制方法的假设条件。
S250，对经过校正的共检波点道集进行基于固定速度的线性噪音去噪处理；
S260，基于各道的偏移距校正量对经过去噪处理的共检波点道集的各道地震数据进行偏移距逆校正，以变换回原道集。
实施例三
图8为未应用本发明的基于偏移距校正的线性噪音压制处理方法的叠加剖面，图9为本实施例应用了本发明的基于偏移距校正的线性噪音压制方法的叠加剖面。通过对比可知，后者的剖面品质得到了明显改善。
该实施例充分证明了本发明在地震资料的处理中能提高去噪的质量，改善资料成像和保真度，为后续进行储层地质分析奠定了良好的基础，具有很好的应用前景。
应该理解的是，本发明所公开的实施例不限于这里所公开的特定处理步骤或材料，而应当延伸到相关领域的普通技术人员所理解的这些特征的等同替代。还应当理解的是，在此使用的术语仅用于描述特定实施例的目的，而并不意味着限制。
说明书中提到的“实施例”意指结合实施例描述的特定特征、或特性包括在本发明的至少一个实施例中。因此，说明书通篇各个地方出现的短语“实施例”并不一定均指同一个实施例。
本领域的技术人员应该明白，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。本领域的技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储介质(RAM)、内存、只读存储介质(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM或技术领域内所公知的任意其它形式的存储介质中。
虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所公开的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的保护范围，仍须以所附的权利要求书所界定的范围为准。

Claims (9)

1.一种地震勘探资料噪音去噪方法，其特征在于，包括以下步骤：
将探区的地震勘探资料的地震数据分选成共检波点道集；
从共检波点道集的各道地震数据的道头信息中提取各道的炮检偏移距和炮检方位角，根据各道的炮检偏移距和炮检方位角以及探区的线性噪音的视速度，确定各道的偏移距校正量；
基于各道的偏移距校正量对共检波点道集的各道地震数据进行偏移距校正，使得其中的非线性相干噪音变换为线性相干噪音；
对经过校正的共检波点道集进行线性噪音去噪处理；
基于各道的偏移距校正量对经过去噪处理的共检波点道集的各道地震数据进行偏移距逆校正，以变换回原道集。
2.如权利要求1所述的地震勘探资料噪音去噪方法，其特征在于：
在将探区的地震勘探资料的地震数据分选成共检波点道集之前，所述方法还包括：对探区的地震勘探资料进行预处理。
3.如权利要求2所述的地震勘探资料噪音去噪方法，其特征在于：
所述预处理包括对地震勘探资料赋上空间坐标信息，以及对地震勘探资料进行静校正和振幅恢复。
4.如权利要求1所述的地震勘探资料噪音去噪方法，其特征在于，根据观测系统的炮检坐标，将探区的地震勘探资料的地震数据分选成共检波点道集。
5.如权利要求1所述的地震勘探资料噪音去噪方法，其特征在于，所述探区的线性噪音包括探区的折射波和/或探区的面波。
6.如权利要求1所述的地震勘探资料噪音去噪方法，其特征在于，所述偏移距校正量为时差校正量。
7.如权利要求6所述的地震勘探资料噪音去噪方法，其特征在于，所述时差校正量通过下式计算：
Figure FDA0002117786440000021
式中，v为探区的线性噪音的视速度，Fi为第i道的炮检偏移距，θi为第i道的炮检方位角，Δti为第i道的时差校正量。
8.一种存储介质，其中存储有用于实现如权利要求1至7中任意一项所述的地震勘探资料噪音去噪方法的程序。
9.一种计算机系统，包括处理器和存储介质，其中所述处理器用于执行存储在所述存储介质中的用于实现如权利要求1至7中任意一项所述的地震勘探资料噪音去噪方法的程序。
CN201910596349.9A 2019-07-03 2019-07-03 地震勘探资料噪音去噪方法、存储介质和计算机系统 Withdrawn CN112180445A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910596349.9A CN112180445A (zh) 2019-07-03 2019-07-03 地震勘探资料噪音去噪方法、存储介质和计算机系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910596349.9A CN112180445A (zh) 2019-07-03 2019-07-03 地震勘探资料噪音去噪方法、存储介质和计算机系统

Publications (1)

Publication Number Publication Date
CN112180445A true CN112180445A (zh) 2021-01-05

Family

ID=73915820

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910596349.9A Withdrawn CN112180445A (zh) 2019-07-03 2019-07-03 地震勘探资料噪音去噪方法、存储介质和计算机系统

Country Status (1)

Country Link
CN (1) CN112180445A (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030055569A1 (en) * 2001-09-14 2003-03-20 Andrew P. Shatilo Method for attenuating noise in seismic data
US20050024990A1 (en) * 2001-12-22 2005-02-03 Andreas Laake Method of seismic surveying and a seismic surveying arrangement
US20050128874A1 (en) * 2003-12-15 2005-06-16 Chevron U.S.A. Inc. Methods for acquiring and processing seismic data from quasi-simultaneously activated translating energy sources
US20070076525A1 (en) * 2005-10-04 2007-04-05 Craft Kenneth L Coherent wave energy removal from seismic data
GB201013228D0 (en) * 2010-08-06 2010-09-22 Geco Technology Bv Seismic acquisition and filtering
CN101893719A (zh) * 2010-04-16 2010-11-24 西安石油大学 一种地震纵横波波场分离与去噪方法
CN102193107A (zh) * 2010-03-05 2011-09-21 西安石油大学 一种地震波场分离与去噪方法
FR2976088A1 (fr) * 2011-05-30 2012-12-07 Cggveritas Services Sa Procede de tomographie non lineaire pour un axe de symetrie principal d'un modele de vitesse anisotrope et dispositif
CN105301653A (zh) * 2014-07-04 2016-02-03 中国石油化工股份有限公司 一种多因素联合识别的钻机噪音压制方法
CN106324677A (zh) * 2015-07-08 2017-01-11 中国石油化工股份有限公司 一种位移场剩余动校正方法及装置
CN107843928A (zh) * 2016-09-21 2018-03-27 中国石油化工股份有限公司 一种基于网格搜索算法的炮点纠偏方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030055569A1 (en) * 2001-09-14 2003-03-20 Andrew P. Shatilo Method for attenuating noise in seismic data
US20050024990A1 (en) * 2001-12-22 2005-02-03 Andreas Laake Method of seismic surveying and a seismic surveying arrangement
US20050128874A1 (en) * 2003-12-15 2005-06-16 Chevron U.S.A. Inc. Methods for acquiring and processing seismic data from quasi-simultaneously activated translating energy sources
US20070076525A1 (en) * 2005-10-04 2007-04-05 Craft Kenneth L Coherent wave energy removal from seismic data
CN102193107A (zh) * 2010-03-05 2011-09-21 西安石油大学 一种地震波场分离与去噪方法
CN101893719A (zh) * 2010-04-16 2010-11-24 西安石油大学 一种地震纵横波波场分离与去噪方法
GB201013228D0 (en) * 2010-08-06 2010-09-22 Geco Technology Bv Seismic acquisition and filtering
GB2482538A (en) * 2010-08-06 2012-02-08 Geco Technology Bv Filtering of seismic data to remove coherent noise from reflection data
FR2976088A1 (fr) * 2011-05-30 2012-12-07 Cggveritas Services Sa Procede de tomographie non lineaire pour un axe de symetrie principal d'un modele de vitesse anisotrope et dispositif
CN105301653A (zh) * 2014-07-04 2016-02-03 中国石油化工股份有限公司 一种多因素联合识别的钻机噪音压制方法
CN106324677A (zh) * 2015-07-08 2017-01-11 中国石油化工股份有限公司 一种位移场剩余动校正方法及装置
CN107843928A (zh) * 2016-09-21 2018-03-27 中国石油化工股份有限公司 一种基于网格搜索算法的炮点纠偏方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘明乾等: "《非纵测线地震资料相干噪声压制》", 石油物探, vol. 50, no. 3, pages 307 - 309 *
崔永福等: "《不规则观测系统数据规则化及在相干噪声 压制中的应用》", 石油物探, vol. 55, no. 4 *

Similar Documents

Publication Publication Date Title
Tian et al. Parabolic-trace time-frequency peak filtering for seismic random noise attenuation
CN110907995B (zh) 井中vsp地震数据的逆时偏移方法及装置
CN113625337B (zh) 一种极浅水高精度地震资料快速成像方法
MX2013012969A (es) Procesamiento de datos sismicos de expansion espacial, metodo y aparato.
CN113805234B (zh) 在被动源地震数据中增强面波的处理方法
CN111736224B (zh) 一种压制叠前地震资料线性干扰方法、存储介质及设备
CN113064203A (zh) 共轭梯度归一化lsrtm方法、系统、存储介质及应用
CN116068619A (zh) 一种自适应的多阶频散面波压制方法、装置及设备
CN111239814B (zh) 基于同相轴分频追踪平滑的浅剖数据机械干扰压制方法
CN112180445A (zh) 地震勘探资料噪音去噪方法、存储介质和计算机系统
CN114428343A (zh) 基于归一化互相关的Marchenko成像方法及系统
US11467302B1 (en) Seismic first breaks onset times determination by interface tracking using level-sets method
CN112213784B (zh) 复杂地表地震数据一次处理快速静校正方法
CN111538088B (zh) 一种海上斜缆波场校正方法
CN111352157B (zh) 一种横波静校正方法及系统
CN109884701B (zh) 地质体散射角度导向深度成像方法
CN116381790A (zh) 一种三维叠前地震记录相干噪声衰减方法
CN110737024B (zh) 地震波信号的处理方法、装置、设备及存储介质
Zhang et al. The ridgelet transform with non-linear threshold for seismic noise attenuation in marine carbonates
CN113075734B (zh) 一种基于信噪比约束的剩余曲率谱计算方法及装置
CN112711074B (zh) 一种地震初至波的去噪方法及装置
CN112444878A (zh) 地震数据分频初至层析静校正方法、存储介质和计算设备
CN117930357A (zh) 三维地震数据噪音压制方法、装置、电子设备及储存介质
Yoneshima et al. Migration‐Based Local Event‐Location Workflow for Ocean‐Bottom Seismometer (OBS) Records in Subduction Zones: A Practical Approach for Addressing a Large Number of Events
CN113917539A (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
WW01 Invention patent application withdrawn after publication

Application publication date: 20210105

WW01 Invention patent application withdrawn after publication