CN109343119A - 低信噪比地震资料有效信号主能量保幅处理方法 - Google Patents

低信噪比地震资料有效信号主能量保幅处理方法 Download PDF

Info

Publication number
CN109343119A
CN109343119A CN201811335341.9A CN201811335341A CN109343119A CN 109343119 A CN109343119 A CN 109343119A CN 201811335341 A CN201811335341 A CN 201811335341A CN 109343119 A CN109343119 A CN 109343119A
Authority
CN
China
Prior art keywords
main energy
seismic data
noise ratio
processing method
relative 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.)
Pending
Application number
CN201811335341.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.)
Beijing Co-Sail Petroleum Technology CoLtd
Original Assignee
Beijing Co-Sail Petroleum Technology CoLtd
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 Beijing Co-Sail Petroleum Technology CoLtd filed Critical Beijing Co-Sail Petroleum Technology CoLtd
Priority to CN201811335341.9A priority Critical patent/CN109343119A/zh
Publication of CN109343119A publication Critical patent/CN109343119A/zh
Pending legal-status Critical Current

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
    • 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/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
    • G01V1/362Effecting static or dynamic corrections; Stacking

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、利用已有或者初步成果建立解释数据库;
2、根据已有钻井和测井信息进行层位标定;
3、选取地质层位明显的典型井的测井声波曲线进行初始合成记录制作,作为后续处理的标准依据;
4、加载原始地震数据和几何库加载;
5、野外静校正处理;
6、信号及噪音频率分析;
7、根据步骤6分析进行主能量段噪音衰减;
8、叠加速度分析;
9、扩展主能量反褶积及初叠加;
10、主能量剩余静校正处理;
11、叠前时间偏移速度分析;
12、扩展主能量叠前时间偏移;
13、对应于新成果制作新合成记录;
14、效果对比,满足地质要求则输出最终结果,否则回到步骤7调整去噪参数、反褶积参数和叠前时间偏移参数重新进行处理工作,直到满足地质要求。
15、输出叠前成像道集和叠加结果。
作为优选,所述的步骤1中所述成果是指收集到的原处理的地震资料结果或者是初步处理后的地震资料结果,该结果可在处理过程中进行跟新。
作为优选,所述的步骤2中所述层位标定是来自钻井分层信息,测井数据由现场第三方测量提供。
作为优选,所述的步骤3中所述初始合成记录,主要是利用步骤1和步骤2数据,选取合适主频的地震子波,进行合成记录制作,以备后续参考。
作为优选,所述的步骤4中由第三方提供原始数据和野外测量成果,加载后形成含有几何关系的内部地震数据库。
作为优选,所述的步骤5主要是利用第三方或者自己计算得到的静校正量对地震数据进行校正。
作为优选,所述的步骤6中的信号及噪音分析是指分析不同频率段信号和噪音的能量分布,确定有效信号的主能量分布范围,后续处理以此能量范围作为约束。
作为优选,所述的步骤7的主能量段噪音衰减为根据不同噪音类型选取针对性处理技术在步骤7确定的主能量段内进行针对性处理。
作为优选,所述的步骤9中扩展是指以步骤6分析得到的主能量所在最小、最大频率为基础,分别向外扩展10Hz,对此频率范围内的地震资料进行处理。
作为优选,所述的步骤13新合成记录是指根据步骤12处理得到的地震资料成果,重新利用步骤2的数据合成新的井旁道合成记录。
作为优选,所述的步骤14根据合成记录和地震数据相关系数最终确定处理参数和处理流程。
本发明的有益效果:在本发明中,以原始记录浅、中、深频谱和分频扫描处理为基础,分析地震记录原始有效信号、噪音的主能量分布范围,参考初始合成记录中地震子波主频率范围,确定原始地震资料信号的主能量分布频段,利用井震匹配作为约束和验证的手段,以扩展主能量频率段地震数据为处理对象,进行保幅叠前去噪、扩展主能量反褶积、主能量剩余静校正、扩展叠前时间偏移等处理。由于在处理全过程中限制了频率范围,更容易确定各个处理步骤中的参数,可有效提高20%以上的工作效率。本发明中的基于井震匹配的低信噪比地震资料有效信号主能量保幅处理方法在地震资料处理中有广泛的应用前景。
附图说明
下面结合附图和具体实施方式来详细说明本发明;
图1为本发明的流程图。
具体实施方式
为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体实施方式,进一步阐述本发明。
实施例1:低信噪比地震资料有效信号主能量保幅处理方法,包含以下步骤:
在步骤101,利用老成果建立解释库,老成果为原先处理的结果,如果没有,则可以利用本次处理过的初步结果。在本步骤中建立的解释数据库主要是为了后续的制作合成记录服务。
在步骤102中,把收集到的井位坐标信息、测井数据和层位信息标定输入到步骤101所建立的解释数据库中,确定井位位置,通过分析过井剖面的构造特征和信噪比,选取过典型井的纵横向地震剖面。
在步骤103中,计算在步骤102中确定的典型井的合成记录,并且根据与井旁地震道相关性选择合适主频的地震子波进行计算。在此步骤中的地震子波主频作为后续处理有效频带的参考之一。
在步骤104中,录入原始地震数据和几何库,并把两者合并形成内部格式,用于后续的资料处理。
在步骤105中,把由第三方计算的野外静校正值(单位:毫秒),或者根据大炮初值计算的静校正值输入步骤104形成的数据,并进行校正。
在步骤106中,分析地震数据的噪音类型和有效信号的频率范围,方法如下:首先按照工区平面位置以2公里X 2公里的间隔选择单炮记录;然后对这些单炮进行噪音分析,包括异常振幅分析,随机噪音分析,浅层折射波分析,多次波分析等;最后进行有效信号主能量分析,这里进行两种频率分析,一是分别分析资料0-1500毫秒、1500-3000毫秒、3000-8000毫秒的浅、中、深层的振幅谱分析,以-20分贝作为阈值,确定资料的频率范围;二是对这些单炮进行频率扫描,如果资料的采样率为2毫秒,则具体扫描频率为0-10赫兹、10-20赫兹,20-40赫兹,30-60赫兹,40-80赫兹,50-100赫兹,60-120赫兹,70-140赫兹,80-160赫兹,90-180赫兹,100-200赫兹,110-220赫兹,120-150赫兹。在这里根据有效信号强弱,结合第一步中结果,确定地震资料有效信号在浅、中、深层的主能量频带分布范围。
在步骤107中,根据步骤106中确定的噪音类型,和有效信号频率分布范围结合步骤103中所用子波主频范围,各向频率的高端、低端扩展10赫兹,在这个扩展的频率范围内进行各种噪音的衰减。
在步骤108中,进行叠加速度分析,这是一迭代过程,初始速度分析网格为1000米X1000米,在反褶积和剩余静校正后的速度分析网格为500米X 500米,并且后一次的速度分析值是在前一次速度分析值的基础上进行。
在步骤109中,以宽展主能量频带内去噪资料作为输入,以地表一致性反褶积作为主要反褶积技术进行处理,在这一步骤中主要参数是反褶积的预测步长,一般按照预测步长8毫秒、12毫秒、16毫秒、20毫秒、24毫秒和48毫秒来测试,对所得结果进行叠加,按照步骤106的手段分析频率,选择合适的反褶积步长。
在步骤110中,输入步骤109中的道集和叠加结果,在扩展频带内对资料进行剩余静校正处理,该处理和步骤108交替进行处理,直到速度值趋于稳定和静校正量趋于一个采样点为止,然后把静校正值加到步骤109数据的道头上。
在步骤111中,输入步骤110中产生的数据,在扩展主能量频段内进行叠前时间偏移速度分析,该步骤和步骤112配合进行,初次速度分析网格为1000米X 1000米,后续逐步加密到500米 X 500米,个别构造复杂地区速度分析网格为250米 X 250米,最终形成偏移速度场。
在步骤112中,输入步骤110的地震数据和步骤111的偏移速度场,在扩展主能量频段内对资料进行叠前时间偏移处理。
在步骤113中,把步骤112中的叠前时间偏移结果输入到步骤101建立的解释库中,重新根据新数据确定地震子波的频率范围并制作新的合成记录。
在步骤114中,分析步骤113中的结果,重点分析两方面的内容:首先分析新地震子波频带宽度和新成像结果频带宽度,正常情况下地震子波的频带宽度应该高于地震子波的频带宽度,如果地震资料的频带宽度低于地震子波的频带宽度,则需要在步骤107、109和步骤112中进一步把频率向两端扩展,直到两者的频带宽度基本一致。其次计算合成记录和井旁地震道的相关系数,如果相关系数平均值在0.85以上,则说明所选频率范围和偏移速度场合适,输出偏移结果用于后续研究;如果相关系数平均值在0.85以下,则重点调整步骤109的反褶积步长和步骤111中的偏移速度场,直到相关系数平均值在0.85以上。
在步骤115中输出最终结果用于后续研究。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (11)

1.低信噪比地震资料有效信号主能量保幅处理方法,其特征在于,包括以下步骤:
(1)、利用已有或者初步成果建立解释数据库;
(2)、根据已有钻井和测井信息进行层位标定;
(3)、选取地质层位明显的典型井的测井声波曲线进行初始合成记录制作,作为后续处理的标准依据;
(4)、加载原始地震数据和几何库加载;
(5)、野外静校正处理;
(6)、信号及噪音频率分析;
(7)、根据步骤(6)分析进行主能量段噪音衰减;
(8)、叠加速度分析;
(9)、扩展主能量反褶积及初叠加;
(10)、主能量剩余静校正处理;
(11)、叠前时间偏移速度分析;
(12)、扩展主能量叠前时间偏移;
(13)、对应于新成果制作新合成记录;
(14)、效果对比,满足地质要求则输出最终结果,否则回到步骤(7)调整去噪参数、反褶积参数和叠前时间偏移参数重新进行处理工作,直到满足地质要求;
(15)、输出叠前成像道集和叠加结果。
2.根据权利要求1所述的低信噪比地震资料有效信号主能量保幅处理方法,其特征在于,所述的步骤(1)中所述成果是指收集到的原处理的地震资料结果或者是初步处理后的地震资料结果,该结果可在处理过程中进行跟新。
3.根据权利要求1所述的低信噪比地震资料有效信号主能量保幅处理方法,其特征在于,所述的步骤(2)中所述层位标定是来自钻井分层信息,测井数据由现场第三方测量提供。
4.根据权利要求1所述的低信噪比地震资料有效信号主能量保幅处理方法,其特征在于,所述的步骤(3)中所述初始合成记录,主要是利用步骤(1)和步骤(2)数据,选取合适主频的地震子波,进行合成记录制作,以备后续参考。
5.根据权利要求1所述的低信噪比地震资料有效信号主能量保幅处理方法,其特征在于,所述的步骤(4)中由第三方提供原始数据和野外测量成果,加载后形成含有几何关系的内部地震数据库。
6.根据权利要求1所述的低信噪比地震资料有效信号主能量保幅处理方法,其特征在于,所述的步骤(5)主要是利用第三方或者自己计算得到的静校正量对地震数据进行校正。
7.根据权利要求1所述的低信噪比地震资料有效信号主能量保幅处理方法,其特征在于,所述的步骤(6)中的信号及噪音分析是指分析不同频率段信号和噪音的能量分布,确定有效信号的主能量分布范围,后续处理以此能量范围作为约束。
8.根据权利要求1所述的低信噪比地震资料有效信号主能量保幅处理方法,其特征在于,所述的步骤(7)的主能量段噪音衰减为根据不同噪音类型选取针对性处理技术在步骤(7)确定的主能量段内进行针对性处理。
9.根据权利要求1所述的低信噪比地震资料有效信号主能量保幅处理方法,其特征在于,所述的步骤(9)中扩展是指以步骤(6)分析得到的主能量所在最小、最大频率为基础,分别向外扩展10Hz,对此频率范围内的地震资料进行处理。
10.根据权利要求1所述的低信噪比地震资料有效信号主能量保幅处理方法,其特征在于,所述的步骤(13)新合成记录是指根据步骤(12)处理得到的地震资料成果,重新利用步骤(2)的数据合成新的井旁道合成记录。
11.根据权利要求1所述的低信噪比地震资料有效信号主能量保幅处理方法,其特征在于,所述的步骤(14)根据合成记录和地震数据相关系数最终确定处理参数和处理流程。
CN201811335341.9A 2018-11-10 2018-11-10 低信噪比地震资料有效信号主能量保幅处理方法 Pending CN109343119A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811335341.9A CN109343119A (zh) 2018-11-10 2018-11-10 低信噪比地震资料有效信号主能量保幅处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811335341.9A CN109343119A (zh) 2018-11-10 2018-11-10 低信噪比地震资料有效信号主能量保幅处理方法

Publications (1)

Publication Number Publication Date
CN109343119A true CN109343119A (zh) 2019-02-15

Family

ID=65314673

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811335341.9A Pending CN109343119A (zh) 2018-11-10 2018-11-10 低信噪比地震资料有效信号主能量保幅处理方法

Country Status (1)

Country Link
CN (1) CN109343119A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113835123A (zh) * 2020-06-24 2021-12-24 中国石油化工集团有限公司 基于地质目标叠前偏移成像的地震采集参数分析方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120243372A1 (en) * 2011-03-25 2012-09-27 Saudi Arabian Oil Company Simultaneous Wavelet Extraction and Deconvolution in the Time Domain
CN103217713A (zh) * 2012-01-19 2013-07-24 中国石油化工股份有限公司 油气勘探地震速度分析数据优化方法
CN104216014A (zh) * 2014-09-25 2014-12-17 郭平 一种地震信号分频处理方法
WO2014207060A1 (en) * 2013-06-28 2014-12-31 Cgg Services Sa Parametervariation improvement for seismic data using sensitivity kernels

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120243372A1 (en) * 2011-03-25 2012-09-27 Saudi Arabian Oil Company Simultaneous Wavelet Extraction and Deconvolution in the Time Domain
CN103217713A (zh) * 2012-01-19 2013-07-24 中国石油化工股份有限公司 油气勘探地震速度分析数据优化方法
WO2014207060A1 (en) * 2013-06-28 2014-12-31 Cgg Services Sa Parametervariation improvement for seismic data using sensitivity kernels
CN104216014A (zh) * 2014-09-25 2014-12-17 郭平 一种地震信号分频处理方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
刘志江: "精细三维地震勘探在大民屯凹陷的应用研究", 《大庆石油地质与开发》 *
刘治凡: "利用VSP资料标定层位", 《新疆石油地质》 *
李振江,等: "松辽盆地南部深层地震资料处理方法研究", 《世界地质》 *
李金山,等: "主能量脉冲反褶积", 《中国海洋大学学报》 *
王世瑞,等: "地震数据目标处理技术在松辽南部的应用", 《石油地球物理勘探》 *
符力耘,等: "库车坳陷复杂高陡构造地震成像研究", 《地球物理学报》 *
解凤英: "含噪声地震资料自动编辑的新算法", 《美国勘探地球物理学家学会第62届年会论文集》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113835123A (zh) * 2020-06-24 2021-12-24 中国石油化工集团有限公司 基于地质目标叠前偏移成像的地震采集参数分析方法
CN113835123B (zh) * 2020-06-24 2024-02-09 中国石油化工集团有限公司 基于地质目标叠前偏移成像的地震采集参数分析方法

Similar Documents

Publication Publication Date Title
CN102323617B (zh) 一种复杂地表的二维地震资料连片处理方法
CN109669212B (zh) 地震数据处理方法、地层品质因子估算方法与装置
CN104237945B (zh) 一种地震资料自适应高分辨处理方法
CN1873441A (zh) 地震数据处理中压制与激发源无关的背景噪声的方法
CN102099705A (zh) 处理地震数据以供解读的系统和方法
EA199900747A1 (ru) Способ обработки сигналов сейсмических данных
CN113640881B (zh) 多偏移距二维横向高分辨率瞬态面波探测方法
CN103645499B (zh) 基于叠后反射波能量统计的地表一致性振幅补偿方法
CN103630932A (zh) 一种地震数据分形保幅方法
CN101852864B (zh) 一种利用地表一致性统计频谱分析技术处理海量地震数据的方法
CN104133244B (zh) 一种震源横波勘探的近地表静校正处理方法
CA2053113C (en) Method for high-resolution seismic recording using detectors planted at shallow depths
CN103217709B (zh) 一种提高地震数据信噪比和分辨率的面波衰减方法
CN109343119A (zh) 低信噪比地震资料有效信号主能量保幅处理方法
Kuehnert et al. Rockfall localization based on inter-station ratios of seismic energy
CN107678065B (zh) 提高地震分辨率的保构造井控空间反褶积方法和装置
CN102778692B (zh) 一种利用地震数据确定品质因数的方法
CN114137606A (zh) 一种稳健的谱模拟反褶积方法
CN112526611A (zh) 表层地震波品质因子的提取方法及装置
CN113064204B (zh) 一种井中地震吸收参数提取及地层岩性检测方法
CN112799132B (zh) 微局部线性噪声压制方法及装置
Han et al. Automatic time picking and velocity determination on full waveform sonic well logs
CN113835123B (zh) 基于地质目标叠前偏移成像的地震采集参数分析方法
CN113589366B (zh) 基于全波形反演的宽频融合建模方法
Urosevic et al. Acquisition of the Nested 3D seismic survey at Harvey

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20190215

RJ01 Rejection of invention patent application after publication