CN103913767A - 多块三维地震数据叠前融合处理能量不均的解决方法 - Google Patents
多块三维地震数据叠前融合处理能量不均的解决方法 Download PDFInfo
- Publication number
- CN103913767A CN103913767A CN201310351715.7A CN201310351715A CN103913767A CN 103913767 A CN103913767 A CN 103913767A CN 201310351715 A CN201310351715 A CN 201310351715A CN 103913767 A CN103913767 A CN 103913767A
- Authority
- CN
- China
- Prior art keywords
- energy
- energy adjusting
- seismic data
- bin
- cmp
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及多块三维地震数据叠前融合处理能量不均的解决方法,首先对叠前CMP道集进行动校正和切除,并计算信噪比和统计每一个CMP面元的有效覆盖次数,根据信噪比和有效覆盖次数的相关函数计算能量调整因子,然后对能量调整因子进行平滑滤波,最后将能量调整因子应用到CMP道集中进行能量调整,最终实施叠前偏移处理。解决了现有技术不能完好解决多块三维地震数据融合处理能量不均造成的偏移划弧问题。
Description
技术领域
本发明涉及石油天然气地震勘探领域,具体涉及一种多块三维地震数据叠前融合处理能量不均的解决方法。
背景技术
叠前时间(或深度)偏移技术已成为复杂地区地震勘探的关键处理技术。当进行多块三维地震数据融合处理时,由于工区局部覆盖次数相差悬殊、能量严重不均,造成叠前偏移成像处理时出现画弧现象。这种假象对地质构造的真实性有着严重的影响。
近几年来,为了解决多块三维地震数据融合处理时由于能量不均引起的偏移划弧问题,通常采用区域振幅均衡、局部能量标定、地震道插值等方法解决能量差异问题。这些方法的缺点在于资料处理过程中人为地改变了数据的振幅相对关系,成果数据的保真度变差,这种解决办法特别是在岩性勘探和高精度勘探中属于慎重使用的方法。
发明内容
本发明的目的是提供一种多块三维地震数据叠前融合处理能量不均的解决方法,用以解决现有技术不能完好解决多块三维地震数据融合处理能量不均造成的偏移划弧问题。
为实现上述目的,本发明的方案是:一种多块三维地震数据叠前融合处理能量不均的解决方法,包括如下步骤:
(1)对完成叠前预处理的多块三维地震数据的CMP道集进行动校正和动校拉伸切除,在动校拉伸切除后的CMP道集上,计算地震数据的信噪比SNR,对经过动校拉伸切除后的CMP道集进行纯波叠加,并提取每一个CMP面元的覆盖次数Nm;
(2)根据步骤(1)中得到的覆盖次数和信噪比,计算基于覆盖次数和信噪比的实际能量调整因子Yms,其中,Yms=Gm/SNR,SNR为信噪比,Gm为基于覆盖次数的能量调整因子,Bm是第m个面元的均方根振幅,M为振幅级别5000。
(3)对步骤(2)中得到的能量调整因子Yms进行平滑,去除异常跳跃值,提取光滑的能量调整因子Gms;
(4)对平滑后得到的能量调整因子Gms进行指数化处理,得到最终的能量调整因子G,并在叠前CMP道集上施加该能量调整因子G,对地震数据叠前道集进行能量调整,其中,
(5)用经过能量调整后的多块三维地震数据进行叠前时间偏移或深度偏移。
步骤(1)中计算地震数据的信噪比SNR的公式为:As为有效信号频带范围内的均方根振幅,A为全频带信号均方根振幅。
步骤(2)中第m个面元的均方根振幅Bm的计算方法为:其中i=1,2,3,…,L,L是某一地震道分析时窗内的总样点数,Xi是该地震道分析时窗内第i个样点的振幅值,m是第m个面元序号。
步骤(3)中对能量调整因子Yms进行平滑的公式为:Gms=aYms+(1-a)G(ms-1),其中,Gms为第m个面元的平滑值,G(ms-1)为第m-1个面元的平滑值,a为平滑系数,a的取值范围为[0,1]。
本发明达到的有益效果:本发明通过动校正和动校拉伸切除,去除拉伸畸变引起的能量异常,并通过计算信噪比和覆盖次数能量调整因子,对能量调整因子进行平滑滤波,并将该能量调整因子应用到CMP道集中,对CMP道集进行了基于信噪比和覆盖次数的能量调整,使得能量分布更加均匀,最后对能量调整后的多块三维地震数据实施叠前偏移处理,消除了多块三维地震数据融合处理时由于覆盖次数不均以及偏移距不规则导致能量不均引起的偏移划弧现象,该方法极大地降低了石油地震勘探的成本,提高了工作效率。
附图说明
图1是为动校正后CMP道集切除;
图2a为能量调整前数据体能量图;
图2b为能量调整后数据体能量图;
图3a为地震数据能量调整前偏移图;
图3b为地震数据能量调整后偏移图。
具体实施方式
下面结合附图对本发明做进一步详细的说明。
本发明的具体方案如下:
(1)动校正和切除
在做完叠前去噪、静校正、反褶积等精细预处理的多块三维地震数据的CMP(共反射点)道集上,首先对道集采用准确的叠加速度进行动校正,然后对道集进行动校拉伸切除,去除拉伸畸变引起的能量异常,消除异常能量引起的统计误差。
如图1所示,图中下方的标记线为合适的切除,相对上方的标记线比较去除动校拉伸畸变频率要彻底。
(2)计算数据的覆盖次数
对叠前CMP(共反射点)道集进行纯波叠加,在纯波叠加数据体道头中提取每一个CMP(共反射点)面元的覆盖次数Nm。
(3)基于覆盖次数的能量均衡
将覆盖次数的能量因子用于叠加数据体,根据覆盖次数进行能量调整。设A0为未作能量调整前的纯波叠加数据体的某一CMP面元的能量,Am为根据覆盖次数进行能量调整后的叠加数据体上某一CMP面元的能量,按公式Am=A0/Nm对叠加数据体进行能量调整,其中Nm为覆盖次数。
如图2(a),能量调整前,能量分布差异大;如图2(b),能量调整后,能量分布均匀。
(4)计算统计时窗内的能量调整因子
对做完基于覆盖次数能量调整后的纯波叠加数据体进行均方根振幅统计。首先修改纯波叠加数据体道头信息,将炮号修改成inline号,检波点道头修改成inline+crossline的组合号,然后在定义时窗内,计算第m个面元的均方根振幅,其中i=1,2,3,…,L;L是某一地震道分析时窗内的总样点数;Xi是该地震道分析时窗内第i个样点的振幅值;m是第m个面元序号;Bm是第m个CMP面元的均方根振幅。对每个输入CMP面元地震道,以M=5000为振幅级别,给定一个以dB为单位的增益因子,计算公式为Gm是统计时窗内基于覆盖次数的能量调整因子。
(5)计算统计时窗内的地震数据的信噪比
在动校切除后的CMP(共反射点)道集上,采用能量比值法计算信噪比。在给定的时窗内,先计算有效信号频带范围内的均方根振幅AS,再计算全频带信号均方根振幅A,则信噪比可以表示为
(6)计算与信噪比有关的能量调整因子
根据得到的基于覆盖次数的能量调整因子Gm,计算基于覆盖次数和信噪比的能量调整因子Yms,计算公式为Yms=Gm/SNR,由该公式可知:信噪比高则相对能量因子要变小,信噪比低则变高。
(7)平滑能量调整因子
对能量因子进行平滑,去除异常值,得到光滑的能量调整因子。采用指数平滑方式对振幅因子进行平滑,平滑公式:Gms=aYms+(1-a)G(ms-1),其中Gms为第m个面元的平滑值,G(ms-1)为第m-1个面元的平滑值,a为平滑系数,起取值范围为[0,1],Yms为进行平滑处理前基于覆盖次数和信噪比的能量调整因子。
(8)地震数据叠前道集能量调整
对叠前道集进行能量调整,在叠前CMP道集上施加能量调整因子,其中Gms为步骤(7)中对Yms进行平滑得到的能量调整因子。
(9)叠前时间偏移或深度偏移
用经过能量调整后的多块三维地震数进行叠前时间偏移或深度偏移,解决多块三维地震数据融合处理时由于覆盖次数不均以及偏移距不规则导致能量不均引起的偏移划弧问题。
如图3(a),能量调整前,有偏移画弧现象;如图3(b),能量调整后,消除画弧现象,还原地质构造的真实性。
Claims (5)
1.一种多块三维地震数据叠前融合处理能量不均的解决方法,其特征在于,包括如下步骤:
(1)对完成叠前预处理的多块三维地震数据的CMP道集进行动校正和动校拉伸切除,计算地震数据的信噪比SNR;对经过动校拉伸切除后的CMP道集进行纯波叠加,提取每一个CMP面元的覆盖次数Nm;
(2)根据覆盖次数对经过纯波叠加的数据体进行基于覆盖次数的能量调整,根据信噪比计算实际能量调整因子Yms,其中,Yms=Gm/SNR,SNR为信噪比,Gm为基于覆盖次数的能量调整因子,Bm是第m个面元的均方根振幅,M为振幅级别5000;
(3)对步骤(2)中得到的能量调整因子Yms进行平滑,去除异常跳跃值,提取光滑的能量调整因子Gms;
(4)对平滑后得到的能量调整因子Gms进行指数化处理,得到最终的能量调整因子G,并在叠前CMP道集上施加该能量调整因子G,对地震数据叠前道集进行能量调整,其中,
(5)用经过能量调整后的多块三维地震数据进行叠前时间偏移或深度偏移。
2.根据权利要求1所述的多块三维地震数据叠前融合处理能量不均的解决方法,其特征在于,步骤(1)中计算地震数据的信噪比SNR的公式为:As为有效信号频带范围内的均方根振幅,A为全频带信号均方根振幅。
3.根据权利要求1所述的多块三维地震数据叠前融合处理能量不均的解决方法,其特征在于,步骤(2)中,根据公式Am=A0/Nm对经过纯波叠加的数据体进行能量调整,其中,Nm为一个CMP面元的覆盖次数,A0为未作能量调整前的纯波叠加数据体中一个CMP面元的能量,Am为进行能量调整后的纯波叠加数据体中一个CMP面元的能量。
4.根据权利要求1所述的多块三维地震数据叠前融合处理能量不均的解决方法,其特征在于,步骤(2)中第m个面元的均方根振幅Bm的计算方法为:其中i=1,2,3,…,L,L是某一地震道分析时窗内的总样点数,Xi是该地震道分析时窗内第i个样点的振幅值,m是第m个面元序号。
5.根据权利要求1所述的多块三维地震数据叠前融合处理能量不均的解决方法,其特征在于,步骤(3)中对能量调整因子Yms进行平滑的公式为:Gms=aYms+(1-a)G(ms-1),其中,Gms为第m个面元的平滑值,G(ms-1)为第m-1个面元的平滑值,a为平滑系数,a的取值范围为[0,1]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310351715.7A CN103913767A (zh) | 2013-08-13 | 2013-08-13 | 多块三维地震数据叠前融合处理能量不均的解决方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310351715.7A CN103913767A (zh) | 2013-08-13 | 2013-08-13 | 多块三维地震数据叠前融合处理能量不均的解决方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103913767A true CN103913767A (zh) | 2014-07-09 |
Family
ID=51039591
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310351715.7A Pending CN103913767A (zh) | 2013-08-13 | 2013-08-13 | 多块三维地震数据叠前融合处理能量不均的解决方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103913767A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105137481A (zh) * | 2015-08-24 | 2015-12-09 | 中国石油天然气集团公司 | 一种调整用于叠前偏移成像的地震数据能量一致性的方法 |
CN105425295A (zh) * | 2015-12-18 | 2016-03-23 | 中国石油集团东方地球物理勘探有限责任公司 | 地震数据分阶段叠加权值获取方法及分阶段叠加方法 |
CN109471170A (zh) * | 2018-09-19 | 2019-03-15 | 中国海洋石油集团有限公司 | 一种地震数据处理方法和装置 |
CN110824555A (zh) * | 2018-08-09 | 2020-02-21 | 中国石油天然气股份有限公司 | 地震能量均衡方法及装置、计算机可读存储介质 |
CN112558160A (zh) * | 2020-11-06 | 2021-03-26 | 中国海洋石油集团有限公司 | 一种方位差异三维地震叠前融合处理方法和*** |
CN113109872A (zh) * | 2021-04-07 | 2021-07-13 | 中国煤炭地质总局物测队 | 一种基于覆盖次数加权的叠前振幅补偿处理方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090187347A1 (en) * | 2006-06-30 | 2009-07-23 | Erik Hans Saenger | VH Signal Integration Measure for Seismic Data |
US20100265797A1 (en) * | 2009-04-17 | 2010-10-21 | Robertsson Johan O A | Interferometric seismic data processing |
CN102023310A (zh) * | 2009-09-19 | 2011-04-20 | 中国石油化工集团公司 | 三维观测***均匀性分析方法 |
CN102998704A (zh) * | 2012-12-13 | 2013-03-27 | 石颖 | 一种地球物理勘探地震资料处理方法 |
-
2013
- 2013-08-13 CN CN201310351715.7A patent/CN103913767A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090187347A1 (en) * | 2006-06-30 | 2009-07-23 | Erik Hans Saenger | VH Signal Integration Measure for Seismic Data |
US20100265797A1 (en) * | 2009-04-17 | 2010-10-21 | Robertsson Johan O A | Interferometric seismic data processing |
CN102023310A (zh) * | 2009-09-19 | 2011-04-20 | 中国石油化工集团公司 | 三维观测***均匀性分析方法 |
CN102998704A (zh) * | 2012-12-13 | 2013-03-27 | 石颖 | 一种地球物理勘探地震资料处理方法 |
Non-Patent Citations (4)
Title |
---|
吴清岭等: "基于覆盖次数的叠前振幅归一化处理在大庆油田的应用", 《大庆石油地质与开发》 * |
张学涛等: "CMP域逐点自动能量调整法在地震资料连片处理中的应用", 《油气地球物理》 * |
苏世龙等: "两种数据规则化地震处理技术应用探讨", 《勘探地球物理进展》 * |
郭树祥: "埕岛桩海地区连片地震资料的数据规则化处理", 《石油物探》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105137481A (zh) * | 2015-08-24 | 2015-12-09 | 中国石油天然气集团公司 | 一种调整用于叠前偏移成像的地震数据能量一致性的方法 |
CN105425295A (zh) * | 2015-12-18 | 2016-03-23 | 中国石油集团东方地球物理勘探有限责任公司 | 地震数据分阶段叠加权值获取方法及分阶段叠加方法 |
CN110824555A (zh) * | 2018-08-09 | 2020-02-21 | 中国石油天然气股份有限公司 | 地震能量均衡方法及装置、计算机可读存储介质 |
CN109471170A (zh) * | 2018-09-19 | 2019-03-15 | 中国海洋石油集团有限公司 | 一种地震数据处理方法和装置 |
CN112558160A (zh) * | 2020-11-06 | 2021-03-26 | 中国海洋石油集团有限公司 | 一种方位差异三维地震叠前融合处理方法和*** |
CN112558160B (zh) * | 2020-11-06 | 2024-05-28 | 中国海洋石油集团有限公司 | 一种方位差异三维地震叠前融合处理方法和*** |
CN113109872A (zh) * | 2021-04-07 | 2021-07-13 | 中国煤炭地质总局物测队 | 一种基于覆盖次数加权的叠前振幅补偿处理方法 |
CN113109872B (zh) * | 2021-04-07 | 2022-06-21 | 中国煤炭地质总局物测队 | 一种基于覆盖次数加权的叠前振幅补偿处理方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103913767A (zh) | 多块三维地震数据叠前融合处理能量不均的解决方法 | |
CN102323617B (zh) | 一种复杂地表的二维地震资料连片处理方法 | |
CN102778693B (zh) | 一种基于反射波层拉平提取并消除的绕射波分离处理方法 | |
CN102998704B (zh) | 一种地球物理勘探地震资料处理方法 | |
CN108196305B (zh) | 一种山地静校正方法 | |
CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
CN109307890A (zh) | 基于上下行波场分解的逆时偏移方法及*** | |
CN103616723A (zh) | 基于avo特征的crp道集真振幅恢复方法 | |
CN106033125B (zh) | 压制叠前大角度道集干涉的提频方法 | |
CN103713323A (zh) | 一种全方位各向异性保幅成像与抽道集方法 | |
CN104483704B (zh) | 基于avo异常类型约束的剩余相位校正方法 | |
CN104570124B (zh) | 一种适合井间地震大角度反射条件的延拓成像方法 | |
CN103984011A (zh) | 一种动态q补偿偏移方法 | |
CN103792579B (zh) | 一种压制动校拉伸的动校正方法 | |
CN109765615A (zh) | 一种地层品质因子反演方法及装置 | |
CN102914796A (zh) | 一种基于高斯束的获取纵横波偏移速度的控制方法 | |
CN104932010A (zh) | 一种基于近道镶边稀疏Radon变换的绕射波分离方法 | |
CN104820242A (zh) | 一种面向叠前反演的道集振幅分频补偿方法 | |
CN108845357A (zh) | 一种基于同步挤压小波变换估计地层等效品质因子的方法 | |
CN104330826A (zh) | 一种去除复杂地表条件下多种噪音的方法 | |
CN103954998B (zh) | 基于avo特征的剩余振幅补偿方法 | |
CN102053260B (zh) | 获得地震纵波的方位速度的方法及处理地震数据的方法 | |
CN104749623B (zh) | 一种地震资料成像处理方法 | |
CN103984016B (zh) | 转换波各向异性振幅随入射角变化道集抽取方法 | |
CN105093318B (zh) | 一种自适应波动方程波场延拓静校正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140709 |
|
RJ01 | Rejection of invention patent application after publication |