CN114779334A - 一种基于统计理论模型的地表一致性振幅处理方法 - Google Patents

一种基于统计理论模型的地表一致性振幅处理方法 Download PDF

Info

Publication number
CN114779334A
CN114779334A CN202210328218.4A CN202210328218A CN114779334A CN 114779334 A CN114779334 A CN 114779334A CN 202210328218 A CN202210328218 A CN 202210328218A CN 114779334 A CN114779334 A CN 114779334A
Authority
CN
China
Prior art keywords
amplitude
time window
cmp
sampling point
gather
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.)
Granted
Application number
CN202210328218.4A
Other languages
English (en)
Other versions
CN114779334B (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.)
Tongji University
Original Assignee
Tongji University
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 Tongji University filed Critical Tongji University
Priority to CN202210328218.4A priority Critical patent/CN114779334B/zh
Publication of CN114779334A publication Critical patent/CN114779334A/zh
Application granted granted Critical
Publication of CN114779334B publication Critical patent/CN114779334B/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
    • G01V1/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • 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

本发明公开了一种基于统计理论模型的地表一致性振幅处理方法，属于地球物理勘探技术领域，包括：获取NMO(动校正)后的CMP(共中心点)道集；对CMP道集开时窗，滑动时窗，统计时窗在每个位置上的能量趋势，以作为CMP道集中与时窗中心相对应的样点经振幅处理后应达到的能量值；在约束条件下求取时窗中心采样点的振幅补偿系数，以使CMP道集中相对应的样点经振幅补偿后达到相应的所述能量值；输出CMP道集中各样点经振幅补偿后的振幅处理结果。本发明基于对数据的统计分析，自适应数据的特点，校正来自同一稳定反射层的子波的振幅畸变，无需反褶积，校正在单个CMP上进行，方法灵活高效，成本低。

Description

一种基于统计理论模型的地表一致性振幅处理方法
技术领域
本发明涉及地球物理勘探技术领域，特别涉及一种基于统计理论模型的地表一致性振幅处理方法。
背景技术
地震勘探利用弹性波估计地下岩石介质构造和弹性参数，对地震勘探的抽象描述是：地下介质构成一个系统，地震子波是输入，观测到的波场是输出。数字地震处理的物理基础是地震波在实际介质中的激发、传播和接收过程，而影响地震波波形的因素可归为与地表、反射界面、传播和采集有关的4方面因素。勘探地震学研究地震波旅行时、振幅和相位与地下介质的关系，任何其他影响因素都是希望去除的因素，地表一致性校正是用来消除地表和接收因素的技术。
反射子波的振幅理论上应仅由地下介质和传播过程决定，但由于源检端近地表因素的影响，以及源检本身的影响，导致振幅畸变。地表一致性振幅校正的本质是剥离掉近地表因素以及源检本身的影响，只留下地下介质和传播过程的影响。
褶积是一个包容性十分广泛的议题，什么样的褶积，就对应一种反褶积，由此导致了反褶积的普遍性。比如脉冲反褶积、预测反褶积、子波整形反褶积等，地表一致性反褶积也是一种典型的反褶积技术。反褶积实际上是褶积和线性反演理论的结合，本质上是在Bayes反演理论框架下，估计反子波，从而消除褶积模型中不想要的因素。把近地表和地下介质看作滤波器，TanerandKoehler(1981)引入了地表一致性褶积模型，并通过求解各个分量再反褶积，达到地表一致性校正的目的。
然而，基于褶积模型迭代地求解各个滤波冲激响应需要遍历数据集，其计算量十分巨大(CaryandLorentz,1993；Guo andZhou,2001；VanVossen etal.,2006)，计算时间长，成本高。
发明内容
针对现有技术存在的常规地表一致性振幅处理中反褶积方法的计算量大、计算时间长的问题，本发明的目的在于提供一种基于统计理论模型的地表一致性振幅处理方法。
为实现上述目的，本发明的技术方案为：
第一方面，本发明提供一种基于统计理论模型的地表一致性振幅处理方法，包括以下步骤：
获取NMO后的CMP道集；
对所述CMP道集开时窗，滑动时窗，统计时窗在每个位置上的能量趋势，以作为所述CMP道集中与时窗中心相对应的样点经振幅处理后应达到的能量值；
在约束条件下求取时窗中心采样点的振幅补偿系数，以使所述CMP道集中相对应的样点经振幅补偿后达到相应的所述能量值；
输出所述CMP道集中各样点经振幅补偿后的振幅处理结果。
优先权，所述NMO后的CMP道集获取步骤为：将经过NMO处理的CMP道集作为地表一致性振幅处理的输入数据。
优先权，所述能量趋势为时窗内各道采样点的能量和的中值。
优先权，所述时窗大小为5×20，即每道5个时间采样点，共20道。
优先权，用于求取时窗中心采样点的振幅补偿系数的目标泛函为：
Figure BDA0003572230690000021
Figure BDA0003572230690000022
其中，ci是振幅补偿系数，cr是当前道，t1和t2是时窗的时间范围，
Figure BDA0003572230690000023
是求中位数，si(t)是输入的CMP的第i道。
优先权，通过高斯-牛顿法求解非线性方程组，即得到所述CMP道集中每个样点的振幅补偿系数。
第二方面，本发明还提供一种电子设备，包括存储有可执行程序代码的存储器以及与所述存储器耦合的处理器；其中，所述处理器调用所述存储器中存储的可执行程序代码，执行如上所述的方法。
第三方面，本发明还提供一种计算机可读存储介质，存储有计算机程序，所述计算机程序被处理器运行时执行如上所述的方法。
采用上述技术方案，本发明的有益效果在于：与传统地表一致性振幅处理方法相比，本发明方法从对数据的统计分析出发，自适应CMP道集的特点，基于对能量趋势的统计，用高斯-牛顿法求解振幅补偿系数使采样点能量达到统计的能量趋势上。本方法是数据驱动的方法，对单个CMP校正，避免了常规反褶积方法遍历数据集求解每道的反褶积分量，提高了地表一致性振幅处理的效率，降低了计算成本。
附图说明
图1为本发明实施例一的方法流程图；
图2为本发明实施例一中模拟的无源检端复杂因素影响的CMP道集示意图；
图3为本发明实施例一中模拟的受源检端复杂因素影响的CMP道集示意图；
图4为本发明实施例一中输入的实际数据示意图；
图5为对图3所示数据进行振幅补偿处理后的数据示意图；
图6为图4的实际数据进行振幅补偿处理后CMP道集示意图；
图7为本发明实施例一中NMO后的含噪CMP道集示意图；
图8为对图7进行振幅补偿处理后CMP道集示意图；
图9为本发明实施例二的结构示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明。在此需要说明的是，对于这些实施方式的说明用于帮助理解本发明，但并不构成对本发明的限定。此外，下面所描述的本发明各个实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互组合。
需要说明的是，在本发明的描述中，术语“上”、“下”、“左”、“右”、“前”、“后”等指示的方位或位置关系为基于附图所示对本发明结构的说明，仅是为了便于描述本发明的简便，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
对于本技术方案中的“第一”和“第二”，仅为对相同或相似结构，或者起相似功能的对应结构的称谓区分，不是对这些结构重要性的排列，也没有排序、或比较大小、或其他含义。
另外，除非另有明确的规定和限定，术语“安装”、“连接”应做广义理解，例如，连接可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个结构内部的连通。对于本领域的普通技术人员而言，可以根据本发明的总体思路，联系本方案上下文具体情况理解上述术语在本发明中的具体含义。
实施例一
一种基于统计理论模型的地表一致性振幅处理方法，如图1所示，包括以下步骤：
S1、获取NMO(动校正)后的CMP道集。
本实施例中，如图2所示，以子波褶积反射系数为地表一致性振幅处理后的标准数据，即模拟的无源检端复杂因素影响的CMP道集示意图，该标准数据包括长乘宽为200乘100个样点。对图2的标准数据随机施加振幅衰减，例如振幅衰减道数占总道数的30％，标准数据中的振幅衰减范围[0.1:0.5]，如图3所示，将振幅衰减处理后的数据作为要处理的输入CMP道集。要处理的实际数据如图4所示。
S2、对CMP道集开时窗，滑动时窗，统计时窗在每个位置上的能量趋势，以作为CMP道集中与时窗中心相对应的样点经振幅处理后应达到的能量值。
本实施例中，配置时窗大小为5×20，即每道5个时间采样点，共20道。在时窗的每个位置上，每道采样点的能量和的中值即表示该位置上时窗的能量趋势，并以此作为位于时窗中心的采样点经振幅处理后应达到的能量值。
当时窗滑动一个位置时，其中心位置发生变化，即可求得CMP道集中另一个样点经振幅处理后应达到的能量值，依此类推，通过不断滑动时窗，即可求得CMP道集中各个样点经振幅处理后应达到的能量值。
S3、在约束条件下求取时窗中心采样点的振幅补偿系数，以使CMP道集中相对应的样点经振幅补偿后达到相应的能量值。
其中S2的计算步骤较多，因此本实施例中通过以下目标泛函求取时窗中心采样点的振幅补偿系数，从而提高计算速度。
Figure BDA0003572230690000041
Figure BDA0003572230690000042
其中，ci是振幅补偿系数，cr是当前道，t1和t2是时窗的时间范围，
Figure BDA0003572230690000043
是求中位数，si(t)是输入的CMP的第i道。通过高斯-牛顿法求解上述的非线性方程，即可求得时窗中心采样点的振幅补偿系数，通过滑动时窗连续处理，比反褶积方法大幅提高了计算速度和效率。
S4、输出CMP道集中各样点经振幅补偿后的振幅处理结果。
如图5和图6所示，在S3中已经计算出CMP道集中各个样点的振幅补偿系数后，对应的即可输出经振幅补偿后的CMP道集。
另外，本实施例提出的方法步骤对于含有噪音的NMO后的CMP道集也具有较好的地表一致性振幅处理效果。图7是含噪数据，应用本实施例公开的方法步骤进行地表一致性振幅处理后如图8，可见对含噪数据也有好的处理效果。
实施例二
一种电子设备，如图9所示，包括存储有可执行程序代码的存储器以及与该存储器耦合的处理器；其中，处理器调用存储器中存储的可执行程序代码，执行如实施例一公开的方法步骤。
实施例三
一种计算机存储介质，计算机存储介质中存储有计算机程序，计算机程序被处理器运行时执行如实施例一公开的方法步骤。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上结合附图对本发明的实施方式作了详细说明，但本发明不限于所描述的实施方式。对于本领域的技术人员而言，在不脱离本发明原理和精神的情况下，对这些实施方式进行多种变化、修改、替换和变型，仍落入本发明的保护范围内。

Claims (7)

1.一种基于统计理论模型的地表一致性振幅处理方法，其特征在于：包括以下步骤：
获取NMO后的CMP道集；
对所述CMP道集开时窗，滑动时窗，统计时窗在每个位置上的能量趋势，以作为所述CMP道集中与时窗中心相对应的样点经振幅处理后应达到的能量值；
在约束条件下求取时窗中心采样点的振幅补偿系数，以使所述CMP道集中相对应的样点经振幅补偿后达到相应的所述能量值；
输出所述CMP道集中各样点经振幅补偿后的振幅处理结果。
2.根据权利要求1所述的方法，其特征在于：所述NMO后的CMP道集获取步骤为：经过NMO(动校正)处理后的CMP(共中心点)道集为地表一致性振幅处理的输入数据。
3.根据权利要求1所述的方法，其特征在于：所述能量趋势为时窗内各道采样点的能量和的中值。
4.根据权利要求3所述的，其特征在于：用于求取时窗中心采样点的振幅补偿系数的目标泛函为：
Figure FDA0003572230680000011
Figure FDA0003572230680000012
其中，ci是振幅补偿系数，cr是当前道，t1和t2是时窗的时间范围，
Figure FDA0003572230680000013
是求中位数，si(t)是输入的CMP的第i道。
5.根据权利要求4所述的，其特征在于：通过高斯-牛顿法求解非线性方程组，即得到所述CMP道集中每个样点的振幅补偿系数。
6.一种电子设备，其特征在于：包括存储有可执行程序代码的存储器以及与所述存储器耦合的处理器；其中，所述处理器调用所述存储器中存储的可执行程序代码，执行如权利要求1-5任一项所述的方法。
7.一种计算机可读存储介质，存储有计算机程序，其特征在于：所述计算机程序被处理器运行时执行如权利要求1-5任一项所述的方法。
CN202210328218.4A 2022-03-30 2022-03-30 一种基于统计理论模型的地表一致性振幅处理方法 Active CN114779334B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210328218.4A CN114779334B (zh) 2022-03-30 2022-03-30 一种基于统计理论模型的地表一致性振幅处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210328218.4A CN114779334B (zh) 2022-03-30 2022-03-30 一种基于统计理论模型的地表一致性振幅处理方法

Publications (2)

Publication Number Publication Date
CN114779334A true CN114779334A (zh) 2022-07-22
CN114779334B CN114779334B (zh) 2023-06-06

Family

ID=82427885

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210328218.4A Active CN114779334B (zh) 2022-03-30 2022-03-30 一种基于统计理论模型的地表一致性振幅处理方法

Country Status (1)

Country Link
CN (1) CN114779334B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5764516A (en) * 1995-12-29 1998-06-09 Atlantic Richfield Company Method and system for surface-consistent phase and time lag correction of seismic data
CN102109615A (zh) * 2010-11-19 2011-06-29 恒泰艾普石油天然气技术服务股份有限公司 地震资料处理真地表条件下的振幅补偿方法和装置
CN103645499A (zh) * 2013-11-08 2014-03-19 中国石油大学(北京) 基于叠后反射波能量统计的地表一致性振幅补偿方法
WO2016141598A1 (zh) * 2015-03-12 2016-09-15 中国矿业大学(北京) 共姿态道集地表一致性振幅补偿方法及装置
CN106226815A (zh) * 2016-08-15 2016-12-14 中国石油集团川庆钻探工程有限公司地球物理勘探公司 转换波等效速度比谱的获取方法及转换波速度分析方法
US20170176617A1 (en) * 2015-12-18 2017-06-22 Saudi Arabian Oil Company Automated near surface analysis by surface-consistent refraction methods
CN109212606A (zh) * 2017-07-06 2019-01-15 中国石油化工股份有限公司 海上地震资料地表一致性振幅补偿处理方法及系统
CN109932748A (zh) * 2019-03-01 2019-06-25 中国石油天然气集团有限公司 一种地表一致性振幅补偿处理方法、装置及存储介质
CN112379429A (zh) * 2020-11-02 2021-02-19 中国石油天然气集团有限公司 地震数据的振幅补偿方法及装置
CN113109872A (zh) * 2021-04-07 2021-07-13 中国煤炭地质总局物测队 一种基于覆盖次数加权的叠前振幅补偿处理方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5764516A (en) * 1995-12-29 1998-06-09 Atlantic Richfield Company Method and system for surface-consistent phase and time lag correction of seismic data
CN102109615A (zh) * 2010-11-19 2011-06-29 恒泰艾普石油天然气技术服务股份有限公司 地震资料处理真地表条件下的振幅补偿方法和装置
CN103645499A (zh) * 2013-11-08 2014-03-19 中国石油大学(北京) 基于叠后反射波能量统计的地表一致性振幅补偿方法
WO2016141598A1 (zh) * 2015-03-12 2016-09-15 中国矿业大学(北京) 共姿态道集地表一致性振幅补偿方法及装置
US20170176617A1 (en) * 2015-12-18 2017-06-22 Saudi Arabian Oil Company Automated near surface analysis by surface-consistent refraction methods
CN106226815A (zh) * 2016-08-15 2016-12-14 中国石油集团川庆钻探工程有限公司地球物理勘探公司 转换波等效速度比谱的获取方法及转换波速度分析方法
CN109212606A (zh) * 2017-07-06 2019-01-15 中国石油化工股份有限公司 海上地震资料地表一致性振幅补偿处理方法及系统
CN109932748A (zh) * 2019-03-01 2019-06-25 中国石油天然气集团有限公司 一种地表一致性振幅补偿处理方法、装置及存储介质
CN112379429A (zh) * 2020-11-02 2021-02-19 中国石油天然气集团有限公司 地震数据的振幅补偿方法及装置
CN113109872A (zh) * 2021-04-07 2021-07-13 中国煤炭地质总局物测队 一种基于覆盖次数加权的叠前振幅补偿处理方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
TING LEI等: "Surface-Consistent Residual Statics, Phase, and Amplitude Corrections: A Statistical Way" *
周阳等: "一种基于逆时偏移的高效照明补偿成像方法" *
蔡希玲;刘学伟;王彦娟;吕英梅;聂红梅;: "地表一致性统计相关分析法及其应用" *
韦红等: "基于相控的振幅补偿法在渤海B油田储层预测中的应用" *

Also Published As

Publication number Publication date
CN114779334B (zh) 2023-06-06

Similar Documents

Publication Publication Date Title
AU2008346986B2 (en) Spectral shaping inversion and migration of seismic data
US10705238B2 (en) Method and apparatus for processing seismic data
US8855934B2 (en) Noise attenuation of seismic data
WO2008112036A1 (en) Imaging of multishot seismic data
CN110174702B (zh) 一种海上地震数据低频弱信号恢复的方法和系统
CN109765616B (zh) 一种保幅波场延拓校正方法及系统
Chen et al. Iterative deblending using shaping regularization with a combined PNMO-MF-FK coherency filter
CN102819043A (zh) 阵列信号随机噪声自适应模型去噪方法
EP2113792A1 (en) Spectral shaping inversion and migration of seismic data
CN105116443A (zh) 一种低频信号的能量补偿方法及装置
Lokshtanov Multiple suppression by single channel and multichannel deconvolution in the tau-p domain
Rickett* Successes and challenges in 3D interpolation and deghosting of single-component marine-streamer data
CN113687417A (zh) 一种三维叠前地震数据层间多次波预测和压制方法
CN114779334A (zh) 一种基于统计理论模型的地表一致性振幅处理方法
Ji et al. Noise attenuation methods for point-receiver land seismic data
AU2013216645B2 (en) Spectral shaping inversion and migration of seismic data
Mills et al. Processing ground roll for the study of near-surface Rayleigh wave dispersion
Han et al. Automatic time picking and velocity determination on full waveform sonic well logs
Henley More coherent noise attenuation in the radial trace domain
CN118011473A (zh) 一种高分辨率地震成像处理方法及装置和设备
Momo et al. Pseudo-Vz deghosting of deep-tow seismic data and implications for low-frequency signal estimation
Berkhout* et al. From removing to using ghost reflections
CN113219530A (zh) 非稳态盲反褶积方法及装置
Lee et al. Cross‐equalization for time‐lapse sparker seismic data
Chen et al. Empirical mode decomposition based smoothing in the flattened domain

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