CN105005080A - 一种利用振幅比率属性识别地层圈闭尖灭线的方法 - Google Patents

一种利用振幅比率属性识别地层圈闭尖灭线的方法 Download PDF

Info

Publication number
CN105005080A
CN105005080A CN201510333999.6A CN201510333999A CN105005080A CN 105005080 A CN105005080 A CN 105005080A CN 201510333999 A CN201510333999 A CN 201510333999A CN 105005080 A CN105005080 A CN 105005080A
Authority
CN
China
Prior art keywords
reservoir
pinch
amplitude ratio
trap
ratio attribute
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
CN201510333999.6A
Other languages
English (en)
Other versions
CN105005080B (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.)
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
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, Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co filed Critical China Petroleum and Chemical Corp
Priority to CN201510333999.6A priority Critical patent/CN105005080B/zh
Publication of CN105005080A publication Critical patent/CN105005080A/zh
Application granted granted Critical
Publication of CN105005080B publication Critical patent/CN105005080B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种利用振幅比率属性识别地层圈闭尖灭线的方法,属于圈闭识别技术领域。本发明通过获取相关测井曲线及地球物理参数,建立以测井曲线为基础的正演模型;对所建立的正演模型进行不同储层、不同砂泥结构的一系列去砂模拟实验,确定储层和不同的砂泥岩互层结构对地震响应的影响及贡献度;在精细解释层位的基础上,根据储层和不同的沙泥岩互层结构对地震响应的影响规律及贡献度,优选参数,提取振幅比率属性,突出储层在地震中的变化点,从而识别地层圈闭的尖灭线。本发明通过振幅比率识别地层圈闭尖灭线,能够确定圈闭形态及储层展布特征,为勘探开发的整体评价与部署提供依据。

Description

一种利用振幅比率属性识别地层圈闭尖灭线的方法
技术领域
本发明涉及一种利用振幅比率属性识别地层圈闭尖灭线的方法,属于圈闭识别技术领域。
背景技术
隐蔽性油气藏已经成为当前勘探开发的重点,而在隐蔽性油气藏中地层及与其相关的圈闭类型占剧了较大的比例。在这种情况下,如何定量的刻画地层分布范围(即地层尖灭线)成为了重要的课题。大部分情况下受地震资料分辨率和信噪比等诸多因素的限制,合理、准确的解释地层尖灭线是很困难的,而地质体的各项异性,更直接影响了精细刻画地层尖灭线。现在所使用的识别尖灭线的技术大都建立在精细解释的基础之上,而地震解释受到多方面因素的限制,使剖面上视尖灭点与实际尖灭点差异较大,难以精细描述尖灭线。而其他的识别方法(相位、正演延伸及反射层夹角外推公式等)也更多的是在地震剖面上的刻画尖灭线,无法展现平面上尖灭线的形态。
发明内容
本发明的目的是提供一种利用振幅比率属性识别地层圈闭尖灭线的方法,以解决现有尖灭线识别过程中存在无法展现平面上尖灭线形态的问题。
本发明为解决上述技术问题而提供一种利用振幅比率属性识别地层圈闭尖灭线的方法,该方法包括以下步骤:
1)获取相关测井曲线及地球物理参数;
2)以测井曲线为基础建立正演模型;
3)对所建立的正演模型进行不同储层、不同砂泥结构的一系列去砂模拟实验,确定储层和不同的砂泥岩互层结构对地震响应的影响及贡献度;
4)根据所确定的储层和不同的砂泥岩互层结构对地震响应的影响规律及贡献度优选参数,提取振幅比率属性,确定储层在地震中的变化点,识别地层圈闭的尖灭线。
所述步骤4)中优选出的参数为包含一个完整地震波形的时窗。
所述步骤4)中振幅比率属性是通过计算一个时窗内地震波形中波峰的平均振幅值与波谷的平均幅值的比值得到。
所述时窗是通过储层影响的地震波形时间厚度优选得到。
所述步骤3)中储层和不同的砂泥岩互层结构对地震响应的影响为有砂体时波峰振幅大于无砂体时的波峰振幅。
所述储层对波峰振幅值的贡献度要大于对波谷值的贡献度。
本发明的有益效果是:本发明通过获取相关测井曲线及地球物理参数,建立以测井曲线为基础的正演模型;对所建立的正演模型进行不同储层、不同砂泥结构的一系列去砂模拟实验,确定储层和不同的砂泥岩互层结构对地震响应的影响及贡献度;在精细解释层位的基础上,根据储层和不同的沙泥岩互层结构对地震响应的影响规律及贡献度,优选参数,提取振幅比率属性,突出储层在地震中的变化点,从而识别地层圈闭的尖灭线。本发明通过振幅比率识别地层圈闭尖灭线,能够确定圈闭形态及储层展布特征,为勘探开发的整体评价与部署提供依据。
附图说明
图1是利用振幅比率属性识别地层圈闭尖灭线方法的流程图;
图2是本发明实施例中探区的地震资料的频谱分析图;
图3是本发明实施例中探区南北向过井剖面示意图;
图4是本发明实施例中探区相关资料建立的理论正演模型示意图;
图5是本发明实施例中探区目的层系去砂模拟分析示意图;
图6是本发明实施例中探区目的层系普通振幅属性图;
图7是本发明实施例中探区目的层系振幅比率属性图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的说明。
振幅比率属性,是分别统计时窗内波峰、波谷值,在求取比值的一种属性,大部分情况用于反应储层产生的强振幅。但在分析、总结区域地震响应特征的基础上,它能够敏感的识别储层对地震响应的影响,突出储层形态及尖灭线。正演技术分析地震响应特征,就是利用已知资料(测井、录井、地震、地质)建立符合区域沉积模式的正演模型,根据正演模型和相应的地层物性参数,按照弹性波动理论模拟地震响应,将模拟地震道与实际剖面进行比较分析,根据需要对模型进行迭代,分析、总结各变量参数对地震响应的影响及贡献度。
如图1所示,本发明通过获取相关测井曲线及地球物理参数,建立以测井曲线为基础的正演模型;对所建立的正演模型进行不同储层、不同砂泥结构的一系列去砂模拟实验,确定储层和不同的砂泥岩互层结构对地震响应的影响及贡献度;在精细解释层位的基础上,根据储层和不同的沙泥岩互层结构对地震响应的影响规律及贡献度,优选参数,提取振幅比率属性,突出储层在地震中的变化点,从而识别地层圈闭的尖灭线。可以检验解释结果的正确性及优选对地震响应影响最大的参数,以此为依据选择适合的圈闭识别方法。地震正演将区域沉积模式和地震响应有机的联系起来,使地震反射特征具有双重意义,既具有地球物理意义,又具有明确的地质含义。该方法的具体实施过程如下:
1.获取地球物理模型及参数
对地区的地震资料进行频谱分析,获取地震资料的频谱范围及主频,区域的砂体的岩石物性参数及相关测井曲线(DT.SP等)。
本实施例中获取地球物理模型及参数如图2所示,对探区的地震资料进行频谱分析,可知地震数据有效频带为10-120Hz,主频65Hz,储层速度为2500米/秒、密度1.6克/立方厘米,砂体厚度约为4米,同时获取相关测井曲线(DT,SP等)数据。
2.根据步骤1获取的地震频域资料,包括带宽、主频、储层速度、子波类型等,以及区域砂体的岩石物性参数及相关测井曲线为基础建立准确的正演模型。
本实施例地质模式的指导下,使用相关资料建立的理论正演模型(正极性)如图4所示,从模型中可以看出在目标砂体存在时,虽然砂体反射被下方强反射掩盖,砂体所在位置的波形无法反应砂体的真实形态,但是由于波形叠加,使有砂体时的波谷振幅要大于无砂体时的波谷振幅(黑圈为砂体存在的范围),验证了本发明的可行性。
3.针对步骤2中所建立的正演模型选择重点井位进行以测井为基础的不同储层、不同沙泥结构的一系列去砂模拟实验(正演的一种方法),进行模拟地震道的响应特征与实际地震响应特征的对比、分析,总结储层和不同的沙泥岩互层结构对地震响应的影响规律及贡献度。通过一些列去砂实验得到的地震响应的影响规律因地形的不同而不同,本实施例中的去砂实验模拟的结果如图5所示(左边为去砂模拟地震道,右边为实际地震道),从中可以总结影响规律和贡献度:有砂体时的波峰振幅要大于无砂体时的波峰振幅;砂体储层存在时,目的层位的波峰、波谷振幅值均较强,反之,均减弱,但波峰减弱的数值要远大于波谷,也就是说储层对地震波峰值的影响要远大于波谷值,即储层对波峰振幅值的贡献度要大于对波谷值的贡献度。所以波峰与波谷的振幅比率能够较好刻画地层及部分特殊储层结构(如强反射界面上的薄砂体等)的岩性类圈闭。
4.应用步骤3中得到地震响应规律,在精细解释层位的基础上,根据储层和不同的沙泥岩互层结构对地震响应的影响规律及贡献度,在地震模拟道上读取储层所影响的完整地震波形(即一个波峰与一个波谷)的时间厚度(本实施例中得到的时间厚度为20ms),时窗参数为受储层影响的地震波形时间厚度20ms(振幅比率属性的时窗需要包含完整的地震波形)。然后通过计算地震波形中波峰及波谷各自的平均振幅值(根据采样点运算)及两者比值,提取振幅比率属性。突出储层在地震中的变化点,刻画地层圈闭的尖灭线。图6中所示的为目的层的振幅属性图,可以看出并无储层形态及圈闭边界的轮廓,但图7中的振幅比率属性图能够清晰看出圈闭形态及储层展布特征,分为3个条带,这也与井1、井2的实钻结果较为一致。
上述方法和结果表明,建立符合实际的地质模型和应用符合实际的参数,开展地震正演模拟分析,根据储层和不同的沙泥岩互层结构对地震响应的影响规律及贡献度,优选参数,提取振幅比率属性,突出储层在地震中的变化点,刻画地层圈闭的尖灭线,为勘探开发的整体评价与部署提供依据。

Claims (6)

1.一种利用振幅比率属性识别地层圈闭尖灭线的方法,其特征在于,该方法包括以下步骤:
1)获取相关测井曲线及地球物理参数;
2)以测井曲线为基础建立正演模型;
3)对所建立的正演模型进行不同储层、不同砂泥结构的一系列去砂模拟实验,确定储层和不同的砂泥岩互层结构对地震响应的影响及贡献度;
4)根据所确定的储层和不同的砂泥岩互层结构对地震响应的影响规律及贡献度优选参数,提取振幅比率属性,确定储层在地震中的变化点,识别地层圈闭的尖灭线。
2.根据权利要求1所述的利用振幅比率属性识别地层圈闭尖灭线的方法,其特征在于,所述步骤4)中优选出的参数为包含一个完整地震波形的时窗。
3.根据权利要求2所述的利用振幅比率属性识别地层圈闭尖灭线的方法,其特征在于,所述步骤4)中振幅比率属性是通过计算一个时窗内地震波形中波峰的平均振幅值与波谷的平均幅值的比值得到。
4.根据权利要求3所述的利用振幅比率属性识别地层圈闭尖灭线的方法,其特征在于,所述时窗是通过储层影响的地震波形时间厚度优选得到。
5.根据权利要求1所述的利用振幅比率属性识别地层圈闭尖灭线的方法,其特征在于,所述步骤3)中储层和不同的砂泥岩互层结构对地震响应的影响为有砂体时波峰振幅大于无砂体时的波峰振幅。
6.根据权利要求5所述的利用振幅比率属性识别地层圈闭尖灭线的方法,其特征在于,所述储层对波峰振幅值的贡献度要大于对波谷值的贡献度。
CN201510333999.6A 2015-06-16 2015-06-16 一种利用振幅比率属性识别地层圈闭尖灭线的方法 Active CN105005080B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510333999.6A CN105005080B (zh) 2015-06-16 2015-06-16 一种利用振幅比率属性识别地层圈闭尖灭线的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510333999.6A CN105005080B (zh) 2015-06-16 2015-06-16 一种利用振幅比率属性识别地层圈闭尖灭线的方法

Publications (2)

Publication Number Publication Date
CN105005080A true CN105005080A (zh) 2015-10-28
CN105005080B CN105005080B (zh) 2018-04-13

Family

ID=54377814

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510333999.6A Active CN105005080B (zh) 2015-06-16 2015-06-16 一种利用振幅比率属性识别地层圈闭尖灭线的方法

Country Status (1)

Country Link
CN (1) CN105005080B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107728230A (zh) * 2017-09-11 2018-02-23 中国石油大港油田勘探开发研究院 一种基于地球物理三联技术的岩性油藏预测方法
CN109917455A (zh) * 2019-02-19 2019-06-21 中国石油天然气集团有限公司 一种地层圈闭的刻画方法及系统
CN110794477A (zh) * 2018-08-03 2020-02-14 中国石油天然气股份有限公司 一种井震结合识别超覆地层尖灭线的方法
CN110888164A (zh) * 2019-12-10 2020-03-17 中国石油大学(华东) 一种确定三角洲前缘席状砂储层砂体尖灭点的方法
CN111025393A (zh) * 2019-12-28 2020-04-17 中海石油(中国)有限公司上海分公司 针对含薄煤层地层的储层预测方法、装置、设备及介质
CN111856569A (zh) * 2019-04-30 2020-10-30 中国石油天然气集团有限公司 地层砂体预测方法及装置
CN112394398A (zh) * 2019-08-13 2021-02-23 中国石油化工股份有限公司 一种基于岩相与地震相的地层尖灭线识别方法
CN112922589A (zh) * 2021-02-03 2021-06-08 中国石油天然气股份有限公司 尖灭线确定方法、装置、终端及存储介质
US11175422B2 (en) 2017-07-05 2021-11-16 Petrochina Company Limited Method and apparatus for identifying low permeable conglomerate diagenetic trap
CN114578452A (zh) * 2022-05-07 2022-06-03 东营市京川石油技术有限公司 一种定量计算地下反射系数振幅占比的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6745129B1 (en) * 2002-10-29 2004-06-01 The University Of Tulsa Wavelet-based analysis of singularities in seismic data
CN103513278A (zh) * 2012-06-19 2014-01-15 中国石油化工股份有限公司 利用地震波组厚度进行储层预测的方法
US20140177390A1 (en) * 2012-12-20 2014-06-26 General Electric Company Modeling of parallel seismic textures

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6745129B1 (en) * 2002-10-29 2004-06-01 The University Of Tulsa Wavelet-based analysis of singularities in seismic data
CN103513278A (zh) * 2012-06-19 2014-01-15 中国石油化工股份有限公司 利用地震波组厚度进行储层预测的方法
US20140177390A1 (en) * 2012-12-20 2014-06-26 General Electric Company Modeling of parallel seismic textures

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
孙志华 等: "层序地层学技术方法应用初探", 《石油地球物理勘探》 *
牟中海 等: "《地震地层学》", 31 March 2013 *
邵雨 等: "岩性圈闭识别与描述方法研究", 《CPS/SEG2004国际地球物理会议论文集》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11175422B2 (en) 2017-07-05 2021-11-16 Petrochina Company Limited Method and apparatus for identifying low permeable conglomerate diagenetic trap
CN107728230A (zh) * 2017-09-11 2018-02-23 中国石油大港油田勘探开发研究院 一种基于地球物理三联技术的岩性油藏预测方法
CN107728230B (zh) * 2017-09-11 2019-06-25 中国石油大港油田勘探开发研究院 一种基于地球物理三联技术的岩性油藏预测方法
CN110794477A (zh) * 2018-08-03 2020-02-14 中国石油天然气股份有限公司 一种井震结合识别超覆地层尖灭线的方法
CN109917455A (zh) * 2019-02-19 2019-06-21 中国石油天然气集团有限公司 一种地层圈闭的刻画方法及系统
CN111856569A (zh) * 2019-04-30 2020-10-30 中国石油天然气集团有限公司 地层砂体预测方法及装置
CN112394398A (zh) * 2019-08-13 2021-02-23 中国石油化工股份有限公司 一种基于岩相与地震相的地层尖灭线识别方法
CN110888164A (zh) * 2019-12-10 2020-03-17 中国石油大学(华东) 一种确定三角洲前缘席状砂储层砂体尖灭点的方法
CN111025393A (zh) * 2019-12-28 2020-04-17 中海石油(中国)有限公司上海分公司 针对含薄煤层地层的储层预测方法、装置、设备及介质
CN112922589A (zh) * 2021-02-03 2021-06-08 中国石油天然气股份有限公司 尖灭线确定方法、装置、终端及存储介质
CN112922589B (zh) * 2021-02-03 2023-08-22 中国石油天然气股份有限公司 尖灭线确定方法、装置、终端及存储介质
CN114578452A (zh) * 2022-05-07 2022-06-03 东营市京川石油技术有限公司 一种定量计算地下反射系数振幅占比的方法

Also Published As

Publication number Publication date
CN105005080B (zh) 2018-04-13

Similar Documents

Publication Publication Date Title
CN105005080A (zh) 一种利用振幅比率属性识别地层圈闭尖灭线的方法
CN105445800B (zh) 一种厚层砂体顶部分异岩性油藏的识别方法
CN110412649B (zh) 一种单期水下分流河道的识别方法
CN104516018B (zh) 一种地球物理勘探中岩性约束下的孔隙度反演方法
CN108802812A (zh) 一种井震融合的地层岩性反演方法
CN109425896A (zh) 白云岩油气储层分布预测方法及装置
CN104155701B (zh) 一种利用叠前地震资料及井信息的多尺度裂缝预测方法
CN103573250A (zh) 一种计算水平井井眼到地层上下界面距离的方法
CN103149589A (zh) 火成岩油气勘探方法及装置
CN107132573B (zh) 一种应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的方法
CN103733089B (zh) 用于包括不确定性估计的地下表征的系统和方法
CN104142516B (zh) 一种薄单砂层厚度预测方法
CN103792571A (zh) 点约束贝叶斯稀疏脉冲反演方法
CN103410504A (zh) 一种水平井/大斜度井确定地层真电阻率的方法及装置
CN104049275A (zh) 复杂盆缘超剥带隐蔽圈闭边界的识别方法及识别系统
CN104155687A (zh) 一种相控叠后声波阻抗反演方法
CN109655896A (zh) 古河道识别及刻画的方法及系统
CN106291677A (zh) 一种基于匹配追踪方法的叠后声波阻抗反演方法
CN106443774A (zh) 一种提高起伏地表叠前深度偏移成像精度的方法
CN109655894B (zh) 碳酸盐岩古河道地震反演低频模型构建方法及系统
CN106597543A (zh) 一种地层沉积相划分方法
CN105319585A (zh) 一种利用薄层干涉振幅恢复识别油气藏的方法
CN104678432A (zh) 砂砾岩裂缝识别的方法
CN104991270A (zh) 一种利用地震正演模拟的叠置砂体识别方法
CN105607120A (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