CN105277984A - 一种时移地震互约束频率一致性处理方法 - Google Patents
一种时移地震互约束频率一致性处理方法 Download PDFInfo
- Publication number
- CN105277984A CN105277984A CN201410320369.0A CN201410320369A CN105277984A CN 105277984 A CN105277984 A CN 105277984A CN 201410320369 A CN201410320369 A CN 201410320369A CN 105277984 A CN105277984 A CN 105277984A
- Authority
- CN
- China
- Prior art keywords
- seismic
- auto
- data
- time
- reflection point
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于油气勘探地震资料处理领域的一种时移地震频率一致性处理方法。包括:时移地震两期资料数据分选;自相关计算;自相关结果互均化计算;求取滤波算子;用每一个滤波算子对相应共中心反射点或共反射点的两期地震道进行褶积处理,就可以达到两期地震资料互约束频率一致性。本发明优势和特点表现在:第一、方法效果的可靠性;第二、实现了两期资料的互约束处理;第三、可实现两期资料空变的频率匹配。本方法在处理过程中,以两期资料共中心反射点或共反射点为匹配标识,在频率一致性处理中实现了空变的处理,提高了频率一致性处理结果的精度。
Description
技术领域
本发明属于油气勘探地震资料处理领域,特别是是一种时移地震频率一致性处理方法。
现有技术
时移地震主要目的是以地震的手段监测油气储层在油气开发过程中储层参数发生一系列的变化。但是在实际处理过程中,由于不同时期资料采集过程中各种外部因素的影响,会使时移地震剖面上非油藏部分带有不应该有的差异。时移地震互约束频率一致性处理的目的是通过对同一区块不同时间采集的地震资料的处理,消除各种因素的影响,使两期资料的频率有合理的同一性和差异性,保证不同时期地震响应的频率差异仅来自于油气储层的变化。
目前,时移地震互约束频率一致性处理方法不是很多且各有优缺点,目前,检索到的时移地震互约束频率一致性处理技术有以下几种:
1、不同步长反褶积处理方法:首先对两期资料都进行不同反褶积步长的测试处理;对两期资料的测试结果进行分析、对比;通过分析、对比结果选取一组能够使两期资料频率最为接近的反褶积参数对资料进行反褶积处理。此方法的缺点:参数的选取更多是靠经验判断;处理中没有很好的体现两期资料的互约束。
2、匹配滤波技术:把作为参考标准的一组资料称为基础资料,要校正的资料称为监测资料,那么匹配滤波技术可以描述如下:首先,在基础资料和监测资料中不包含地下油藏属性变化的区域内取一个时间窗;然后,在此时间窗内让两种资料尽可能匹配,计算出匹配滤波器;最后,用此滤波器作用于监测资料。此方法的缺点:全区用同一个校正因子,无法满足实际资料的复杂变化情况;处理中没有很好的体现两期资料的互约束;匹配滤波器不仅仅受频率控制,还包含相位、能量等因素,一致性处理因素不单一。
3、滤波法:其做法就是设计一个滤波器,对两期资料进行滤波处理。此方法的缺点:虽然可以统一两期资料的频宽,但无法统一两期资料的频率构成;参数选取缺乏科学监控手段;易损伤储层变化引起的频率差异。
发明内容
本发明的目的是针对现有技术存在的问题,提供一种能够实现空变的互约束频率一致性处理消除非储层因素引起频率变化,从而提高时移地震互约束频率一致性处理的精度,继而提高时移地震对储层变化的监控能力的时移地震互约束频率一致性处理方法。
本发明的处理方法包括:
①输入地震数据:输入时移地震两期数据,数据是经过时差校正、能量校正处理后的共中心点道集或偏移后的共反射点道集。
②数据匹配抽取:按照共中心反射点或共反射点为标识道头,对时移地震两期资料进行分选,将两期地震资料具有相同共中心反射点或共反射点的数据放在同一个共中心反射点或共反射点域内。
③自相关计算:设计一个资料信噪比较高且不包含具有储层变化的反射信息位置的计算时窗,在设计的自相关时窗内对所有地震道进行自相关计算处理,得到每一个地震道的自相关结果。
④自相关结果互均化计算:将同一个共中心反射点或共反射点域的两期地震数据的自相关结果用n次方根的方法进行互均化叠加,使每一个共中心反射点或共反射点只有一个互均化后的自相关结果。
⑤求取滤波算子:利用步骤④计算出的自相关结果,求取其反算子,将这个反算子做为滤波算子。
⑥空变滤波:利用步骤⑤得出的滤波算子,用每一个滤波算子对相应共中心反射点或共反射点的两期地震道进行褶积处理,就可以达到两期地震资料互约束频率一致性。
⑦输出:将滤波处理后的数据输出,并对两期时移地震数据进行分选输出,完成互约束频率一致性处理。
本发明的技术方案进一步包括:
(1)时移地震两期资料数据分选
按照共中心反射点或共反射点为标识道头,对时移地震两期资料进行分选,将两期地震资料具有相同共中心反射点或共反射点的数据放在同一个共中心反射点或共反射点域内。
(2)自相关计算
自相关计算的结果取决于计算时窗内是否包含储层变化信息及时窗内地震数据信噪比高低,所以在自相关计算时,计算时窗选择资料信噪比较高且不包含具有储层变化的反射信息位置。在设计的自相关时窗内对所有地震道进行自相关计算处理,得到每一个地震道的自相关结果。
设e(t)为地震记录,其自相关函数为r(τ),τ为时间延迟。自相关是一个时间序列在不同时间点上波形相似性的量度,这是一个移动求和过程,用下式表示:
(3)自相关结果互均化计算
将同一个共中心反射点或共反射点域的两期地震数据的自相关结果进行互均化叠加,使每一个共中心反射点或共反射点只有一个互均化后的自相关结果。设m为一个共中心点或共反射点自相关结果个数,u为自相关结果,n为设定的次方根参数,均化结果为:
(4)求取滤波算子
对互均化后的每一个自相关结果求取其一个反算子,将这个反算子做为滤波算子。设存在一个滤波算子f,滤波算子与(2)式自相关结果褶积可以得到一个尖脉冲函数δ,即:
S*f=δ……………………………………………………(3)
则f可以表示为:
f为自相关结果S的反算子,就是所求取滤波算子。
(5)空变滤波
用每一个滤波算子对相应共中心反射点或共反射点的两期地震道进行褶积处理,就可以达到两期地震资料互约束频率一致性。
本发明有着其他方法不具备的优势,其具体优势和特点表现在以下几个方面:
第一、方法效果的可靠性。本方法通过对时移地震两期资料进行多道统计的方法求取两期资料共同的滤波算子,通过滤波算子的应用实现两期资料的频率匹配,通过处理实例证明效果明显且稳定。
第二、实现了两期资料的互约束处理。由于本方法在滤波算子求取过程中采用两期资料多道统计的方法,从而在频率一致性处理过程中真正的实现了互约束处理的目的。
第三、可实现两期资料空变的频率匹配。本方法在处理过程中,以两期资料共中心反射点或共反射点为匹配标识,在频率一致性处理中实现了空变的处理,提高了频率一致性处理结果的精度。
附图说明
图1是本发明的一种具体实施例的流程图。
图2是互约束频率一致性处理前时移地震两期资料的叠加剖面,左图为1期资料,右图为2期资料,线框内为质量监控比对的标志层。
图3是互约束频率一致性处理后时移地震两期资料的叠加剖面,左图为1期资料,右图为2期资料,线框内为质量监控比对的标志层。
图4是采用本发明方法处理前、后时移地震两期资料的频谱对比图,其中左图为处理前时移地震两期资料的频谱,右图为处理后时移地震两期资料的频谱,频谱分析范围为图2中线框范围。
具体实施方式
本次以XX油田HG地区两期三维地震资料为目标靶区,应用本方法对该资料进行处理,以验证本方法的效果。该实际资料分别于1992年和2010年完成采集,地震资料时间长度6000ms,时间采样间隔为2ms,采样点数3000。输入的地震数据已经完成了观测系统退化匹配、噪音衰减、能量补偿、反褶积及时差校正工作。
处理过程:(按照图1所示处理流程进行)
1、按照共中心反射点为标识道头,对时移地震两期资料进行分选,将两期地震资料具有相同共中心反射点的数据放在同一个共中心反射点域内。
2、设计一个资料信噪比较高且不包含具有储层变化反射信息的计算时窗,在设计的自相关时窗内对所有地震道进行自相关计算处理,得到每一个地震道的自相关结果。
其过程中设e(t)为地震记录,其自相关函数为r(τ),τ为时间延迟。自相关是一个时间序列在不同时间点上波形相似性的量度,这是一个移动求和过程,如式(1)。
3、将同一个共中心反射点域的两期地震数据的自相关结果用n次方根的方法进行互均化叠加,使每一个共中心反射点只有一个互均化后的自相关结果。
其过程中设m为一个共中心点或共反射点自相关结果个数,u为自相关结果,n为设定的次方根参数,均化结果为式(2)所示。
4、利用上一步骤计算出的自相关结果,求取其反算子,将这个反算子做为滤波算子。其求取过程如式(4)所示。
5、利用前一步骤得出的滤波算子,用每一个滤波算子对相应共中心反射点或共反射点的两期地震道进行褶积处理;
6、将滤波处理后的数据输出,并对两期时移地震数据进行分选。对互约束频率一致性处理前后两期数据分别进行动校正、叠加成像(图2、图4)进行质量监控。
效果分析:图2是互约束频率一致性处理前两期资料剖面对比,可以看出,两期资料频率存在一定的差异,图4中左图是两期资料互约束频率一致性处理前的频谱对比,明显可以看出两期资料频率成分存在不同。图3互约束频率一致性处理后两期资料剖面对比,可以看出,经过处理后两期资料频率非常接近,图4中右图是本方法处理后的两期资料频谱对比,此图表明,经过处理后两期资料的频率成分基本相同。
不难看出:本方法可以较好的消除时移地震两期资料存在的频率差异,利于后期的地震综合解释研究。
Claims (2)
1.一种时移地震互约束频率一致性处理方法,其特征是包括,
①输入地震数据:输入时移地震两期数据,数据是经过时差校正、能量校正处理后的共中心点道集或偏移后的共反射点道集;
②数据匹配抽取:按照共中心反射点或共反射点为标识道头,对时移地震两期资料进行分选,将两期地震资料具有相同共中心反射点或共反射点的数据放在同一个共中心反射点或共反射点域内;
③自相关计算:设计一个资料信噪比较高且不包含具有储层变化的反射信息位置的计算时窗,在设计的自相关时窗内对所有地震道进行自相关计算处理,得到每一个地震道的自相关结果;
④自相关结果互均化计算:将同一个共中心反射点或共反射点域的两期地震数据的自相关结果用n次方根的方法进行互均化叠加,使每一个共中心反射点或共反射点只有一个互均化后的自相关结果;
⑤求取滤波算子:利用步骤④计算出的自相关结果,求取其反算子,将这个反算子做为滤波算子;
⑥空变滤波:利用步骤⑤得出的滤波算子,用每一个滤波算子对相应共中心反射点或共反射点的两期地震道进行褶积处理,达到两期地震资料互约束频率一致性;
⑦输出:将滤波处理后的数据输出,并对两期时移地震数据进行分选输出,完成互约束频率一致性处理。
2.根据权利要求1所述的时移地震互约束频率一致性处理方法,其特征是包括:
在自相关计算时,设e(t)为地震记录,其自相关函数为r(τ),τ为时间延迟,自相关是一个时间序列在不同时间点上波形相似性的量度,这是一个移动求和过程,用下式表示:
在自相关结果互均化计算时,设m为一个共中心点或共反射点自相关结果个数,u为自相关结果,n为设定的次方根参数,均化结果为:
在求取滤波算子时,设存在一个滤波算子f,滤波算子与(2)式自相关结果褶积得到一个尖脉冲函数δ,即:
S*f=δ……………………………………………………(3)
f为自相关结果S的反算子,就是所求取滤波算子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410320369.0A CN105277984B (zh) | 2014-07-04 | 2014-07-04 | 一种时移地震互约束频率一致性处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410320369.0A CN105277984B (zh) | 2014-07-04 | 2014-07-04 | 一种时移地震互约束频率一致性处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105277984A true CN105277984A (zh) | 2016-01-27 |
CN105277984B CN105277984B (zh) | 2018-02-02 |
Family
ID=55147322
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410320369.0A Active CN105277984B (zh) | 2014-07-04 | 2014-07-04 | 一种时移地震互约束频率一致性处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105277984B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105700012A (zh) * | 2016-01-28 | 2016-06-22 | 中国海洋石油总公司 | 一种海上非重复性时移地震分频互均衡处理方法 |
CN105866837A (zh) * | 2016-04-06 | 2016-08-17 | 中国海洋石油总公司 | 一种海上非重复性时移地震多层系互均衡处理方法 |
CN106842296A (zh) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | 一种时移地震互约束频率一致性处理方法 |
CN110824561A (zh) * | 2019-11-14 | 2020-02-21 | 中国石油集团东方地球物理勘探有限责任公司 | 度量时移地震数据可重复性的方法 |
CN111045079A (zh) * | 2019-12-20 | 2020-04-21 | 核工业北京地质研究院 | 一种增强地震反射特征的数据处理方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1111911A (en) * | 1965-01-29 | 1968-05-01 | Atlantic Richfield Co | Seismic exploration methods and apparatus |
CN1230260A (zh) * | 1996-09-13 | 1999-09-29 | Pgs张量公司 | 时滞储层监控方法 |
CN101937101A (zh) * | 2010-08-04 | 2011-01-05 | 中国海洋石油总公司 | 一种鉴定能否实施时移地震的方法 |
CN102147478A (zh) * | 2010-12-29 | 2011-08-10 | 中国海洋大学 | 复杂油藏的叠前低频信号识别方法 |
US20120143510A1 (en) * | 2007-05-25 | 2012-06-07 | Aftab Alam | High resolution attributes for seismic data processing and interpretation |
CN102809762A (zh) * | 2012-08-13 | 2012-12-05 | 成都理工大学 | 基于全频带地震信息挖掘的储层成像技术 |
-
2014
- 2014-07-04 CN CN201410320369.0A patent/CN105277984B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1111911A (en) * | 1965-01-29 | 1968-05-01 | Atlantic Richfield Co | Seismic exploration methods and apparatus |
CN1230260A (zh) * | 1996-09-13 | 1999-09-29 | Pgs张量公司 | 时滞储层监控方法 |
US20120143510A1 (en) * | 2007-05-25 | 2012-06-07 | Aftab Alam | High resolution attributes for seismic data processing and interpretation |
CN101937101A (zh) * | 2010-08-04 | 2011-01-05 | 中国海洋石油总公司 | 一种鉴定能否实施时移地震的方法 |
CN102147478A (zh) * | 2010-12-29 | 2011-08-10 | 中国海洋大学 | 复杂油藏的叠前低频信号识别方法 |
CN102809762A (zh) * | 2012-08-13 | 2012-12-05 | 成都理工大学 | 基于全频带地震信息挖掘的储层成像技术 |
Non-Patent Citations (2)
Title |
---|
郝振江等: "基于相关参数的时移地震互均化质量监控方法研究", 《地球物理学进展》 * |
郭念民等: "非重复性时移地震数据关键处理技术应用研究", 《石油地球物理勘探》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106842296A (zh) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | 一种时移地震互约束频率一致性处理方法 |
CN105700012A (zh) * | 2016-01-28 | 2016-06-22 | 中国海洋石油总公司 | 一种海上非重复性时移地震分频互均衡处理方法 |
CN105866837A (zh) * | 2016-04-06 | 2016-08-17 | 中国海洋石油总公司 | 一种海上非重复性时移地震多层系互均衡处理方法 |
CN110824561A (zh) * | 2019-11-14 | 2020-02-21 | 中国石油集团东方地球物理勘探有限责任公司 | 度量时移地震数据可重复性的方法 |
CN110824561B (zh) * | 2019-11-14 | 2021-11-23 | 中国石油集团东方地球物理勘探有限责任公司 | 度量时移地震数据可重复性的方法 |
CN111045079A (zh) * | 2019-12-20 | 2020-04-21 | 核工业北京地质研究院 | 一种增强地震反射特征的数据处理方法 |
CN111045079B (zh) * | 2019-12-20 | 2021-11-30 | 核工业北京地质研究院 | 一种增强地震反射特征的数据处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105277984B (zh) | 2018-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103376464B (zh) | 一种地层品质因子反演方法 | |
CN105277984A (zh) | 一种时移地震互约束频率一致性处理方法 | |
CN105425289B (zh) | 确定低频波阻抗的方法和装置 | |
CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
CN104483704B (zh) | 基于avo异常类型约束的剩余相位校正方法 | |
CN107065013B (zh) | 一种地震尺度下的层速度确定方法及装置 | |
Song et al. | Microseismic event estimation and velocity analysis based on a source-focusing function | |
CN110687592B (zh) | 一种微地震震相识别初至拾取方法、设备和存储介质 | |
CN104570110A (zh) | 一种基于纵横波匹配的多分量资料联合速度分析方法 | |
CN108845357A (zh) | 一种基于同步挤压小波变换估计地层等效品质因子的方法 | |
CN102692647A (zh) | 一种高时间分辨率的地层含油气性预测方法 | |
CN112327358A (zh) | 一种粘滞性介质中声波地震数据正演模拟方法 | |
CN102928875B (zh) | 基于分数阶傅里叶域的子波提取方法 | |
CN104391324A (zh) | 依赖频率的avo反演前的地震道集动校拉伸校正预处理技术 | |
CN103675904B (zh) | 一种井震匹配解释性目标处理方法及装置 | |
CN117031539A (zh) | 一种自监督深度学习地震数据低频重建方法及系统 | |
Wang et al. | A simple approach to nonstationary source waveform estimation in seismic records | |
CN109655883A (zh) | 一种针对目标的地震分频方法及系统 | |
Xie et al. | An improved ant-tracking workflow based on divided-frequency data for fracture detection | |
CN112578475B (zh) | 基于数据挖掘的致密储层双甜点识别方法 | |
US20230072038A1 (en) | Method and apparatus for removing tube wave interference from optical fiber acoustic wave sensing seismic data | |
CN105388519A (zh) | 一种提高地震资料分辨率的方法 | |
CN106054245B (zh) | 地震资料高分辨率处理定量评价方法 | |
CN113960660B (zh) | 一种基于正演模拟的动校正畸变区域自动识别及切除方法 | |
CN102353991B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |