CN103439741A - 一种零值法单砂体逐层剥离预测方法 - Google Patents

一种零值法单砂体逐层剥离预测方法 Download PDF

Info

Publication number
CN103439741A
CN103439741A CN2013104179355A CN201310417935A CN103439741A CN 103439741 A CN103439741 A CN 103439741A CN 2013104179355 A CN2013104179355 A CN 2013104179355A CN 201310417935 A CN201310417935 A CN 201310417935A CN 103439741 A CN103439741 A CN 103439741A
Authority
CN
China
Prior art keywords
sand body
discrete
layer
interest
seismic
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
CN2013104179355A
Other languages
English (en)
Other versions
CN103439741B (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.)
Petrochina Co Ltd
Original Assignee
Research Institute Of Petroleum Exploration & Development Dagang Oil Field Of Cnpc
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 Research Institute Of Petroleum Exploration & Development Dagang Oil Field Of Cnpc filed Critical Research Institute Of Petroleum Exploration & Development Dagang Oil Field Of Cnpc
Priority to CN201310417935.5A priority Critical patent/CN103439741B/zh
Publication of CN103439741A publication Critical patent/CN103439741A/zh
Application granted granted Critical
Publication of CN103439741B publication Critical patent/CN103439741B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种零值法单砂体逐层剥离预测方法，所述的方法包括以下步骤：(1)离散零值点求取层位标定：根据已知井测井资料建立砂泥岩薄互层速度模型，并进行离散合成记录标定分析，计算每个单砂体反射波的所有零值点；(2)分尺度地震精细解释：对包含目的层的一个波长范围内的等效地层进行小尺度精细层位解释，追踪出复合波的零值点、波峰、波谷层位；(3)分尺度古地形恢复：针对目的层分别进行大尺度、中尺度、小尺度古地形恢复，并进行简化求解；(4)离散零值点指导下的地震属性切片：某一单砂体的离散合成记录零值点处的地震属性切片。采用本发明的预测方法可提高地震储层预测分辨能力，克服薄互层砂体地震分辨率的限制。

Description

一种零值法单砂体逐层剥离预测方法
技术领域
本发明属于石油物探储层预测的技术领域，更具体地说，本发明涉及一种基于三维地震资料的砂体分布预测方法。 
背景技术
在石油地球物理勘探储层预测领域，地震勘探是最有效的勘探方法。砂体三维地震储层预测技术很多，包括递推反演、模型反演、多参数反演等，主要是以已知地质规律和钻井、测井资料、地质模型为约束，对地下岩石空间结构和物理性质进行成像，如井约束波阻抗反演、模型控制波阻抗反演、随机模拟波阻抗反演、曲线重构波阻抗反演、有色反演、分频反演等波阻抗反演，还有各种三维地震属性反演。但是根据瑞利准则，薄层的反射在地震剖面上可以分辨的必要条件是其厚度大于四分之一波长。但根据目前地震资料的分辨率，厚度小于10米的地层是难以分辨的。 
我国东部陆相断陷盆地储层多为砂泥岩薄互层特征，是勘探和开发工作的重点和难点。这类储层厚度薄、横向变化快、沉积过程复杂，单砂层厚度仅为2-5米，砂岩厚度远小于地震分辨率。三维砂泥岩薄互层模型的地震正演结果表明：多个薄砂层反射相互干涉之后形成一个复合波，地震反射记录难以分辨每个单砂层。 
应用地震资料提取属性及地层切片是储层预测的重要手段，但基于目前的属性提取方法，各种属性都只能反映多套砂体的叠合分布，无法预测出每个单砂体的平面分布情况。在实际工作中，需要根据地球物理储层预测的单砂体平面厚度变化来指导水平井的钻探，所以急需一种突破地震分辨率限制的薄储层单砂体分布预测方法。 
发明内容
为了解决现有技术中的上述技术问题，本发明的目的在于提供一种零值法单砂体逐层剥离预测方法。采用本发明的预测方法可提高地震储层预测分辨能力，克服薄互层砂体地震分辨率的限制。 
为了实现上述目的，本发明采用了以下技术方案： 
一种零值法单砂体逐层剥离预测方法，其特征在于包括以下步骤： 
(1)离散零值点求取层位标定：根据已知井测井资料建立砂泥岩薄互层速度模型，并进行 离散合成记录标定分析，计算每个单砂体反射波的所有零值点； 
(2)分尺度地震精细解释：根据合成记录标定的目的层上下宏观范围内可等时追踪的大尺度等效地层作为约束背景；按照目的层上下宏观范围内可等时追踪的最小等效地层解释范围锁定目标；再对包含目的层的一个波长范围内的等效地层进行小尺度精细层位解释，追踪出复合波的零值点、波峰、波谷层位； 
(3)分尺度古地形恢复：针对目的层分别进行大尺度、中尺度、小尺度古地形恢复，锁定一个各尺度古地形背景较一致的单元范围进行简化求解； 
(4)离散零值点指导下的地震属性切片：某一单砂体的离散合成记录零值点处的地震属性切片，该砂体无反射记录贡献，因此反映的只是其它砂体的分布特征，即可剥离掉该砂层对地震复合波的影响； 
其中，所述大尺度是指包含目的层的三级层序单元顶底界面等时范围；中尺度是指包含目的层的两个波长范围；小尺度是指包含目的层的一个波长范围。 
其中，在步骤(2)中，通过对测井曲线求取反射系数，并根据阀值做弱反射系数过滤后，计算出主要砂层的反射系数界面，并依此计算各单砂体的离散反射记录及多砂体复合反射合成记录，并给出各单砂体的合成记录零值点时间。 
本发明的技术方案相比现有技术具有以下有益效果： 
在本发明的预测方法中分尺度古地形恢复分别以目的层上下宏观和微观范围内的等时层位为对象，遵循各尺度古地形一致的原则锁定研究目标范围，将高频古地形恢复结果与沉积微相分析紧密结合，使该项储层预测技术赋予了地震地质综合研究的意义；而且各单砂体离散合成记录零值点的远间隔地震属性切片技术可避免常规砂体标定位置地震属性切片的多砂体干涉叠合效果，能有效剥离单砂体对地震复合波的影响，达到预测单砂体分布的效果。采用本发明的逐层剥离预测方法分别在复杂三维地震地质模型以及X地区进行了试验和应用，理论结合实际取得了较好的效果。 
附图说明
图1a为实施例1所述两层砂体的砂泥岩互层三维地震地质模型，其中图1b为第一层砂体，图1c为第二层砂体。 
图2为实施例1所述三维模型用雷克子波正演形成三维数据体。 
图3为实施例1中根据模型速度计算离散合成记录图。 
图4a是对应第一层砂体复合波的零值点振幅切片；图4b是对应第二层砂体复合波的零值点振幅切片；图4c是对应第一层砂体标定位置145ms处的振幅切片。 
图5是X地区地层对比图。 
图6是X地区XXl井离散合成地震记录分析结果。 
图7a是根据离散合成记录计算的零值在XX井区进行单砂体零值点地震属性切片的结果； 
图7b是根据目的层范围内一个波长的小尺度解释层位进行的古地形恢复结果。 
具体实施方式
以下将结合具体实施例对本发明所述的零值法单砂体逐层剥离技术的发明过程和实际应用方法做进一步的阐述，以帮助本领域的技术人员对本发明的发明构思、技术方案有更完整、准确和深入的理解。 
实施例1 
如附图1所示，本实施例所述的针对X地区某层位地层地质及砂泥岩厚度特点，首先，在小尺度条件下建立一个单层厚度分别为5米和7米，背景泥岩为3700m/s，砂岩速度为4400m/s的三维地震地质模型。如图2所示，用最近该地区实际地震资料的主频为26Hz的雷克子波进行褶积形成三维正演地震数据体。对此正演数据体进行离散合成记录及零值法地震属性切片试验，可验证零值法砂体逐层剥离预测技术的可行性。此基于三维地质地质模型的储层预测方法试验与常规的基于二维模型正演方法的区别在于能有效验证地震储层预测方法的空间预测的符合率。第二，对砂泥岩薄互层速度模型或实际井目的层的测井资料按lms的间隔进行抽稀后计算出离散反射系数，再根据用户自定义的门槛值或缺省值(-0.01-0.01)进行反射系数过滤。用户自定义主频的子波计算出总合成记录及每个单砂层复合反射记录并用不同的颜色表示，同时给出每个反射复合波的所有零值点时间。该离散合成记录方法方便用户分析每个单砂体对总合成记录的贡献及单砂体的反射形式及零值时间，结果如图3和图6所示。第三，将X地区某三级地层作为大尺度研究对象开展宏观沉积体系和沉积相带划分，背景约束目的层储层横向预测范围；以X地区A、B、C油组为小尺度研究对象，开展古地形恢复研究，锁定XX井区为砂层发育区；以A产层为微尺度研究单元，精细追踪地震反射波的零值、峰值及波谷层位，开展高频古地形恢复和沉积微相研究，结果如图5和图7b所示。第四，利用X地区目的层范围上下的各等时地层解释层位分别开展大尺度、中尺度、小尺度 古地形恢复，锁定各尺度古地形背景较一致的XX井区为砂体发育区，进行相控储层预测和零值法地震属性切片应用试验。第五，以离散合成记录计算的反射记录零值点指导进行地震属性切片，在已知三维模型的正演数据体复合波主能量上能够找到相对独立地反映第一套砂体和第二套砂体的振幅切片，它们分别为单层砂体复合波的过零点位置，如图4所示，首先出现的是第二套砂体，然后出现的是第一套砂体；在X地区XX单元对A层底部砂体反射波零值点附近进行地震属性切片，剥离该砂体对复合波的影响，准确预测出A层顶部砂体的平面分布情况，结果如图7a所示。最后，综合小尺度古地形恢复及零值地震属性切片结果，对薄互层储层单砂体平面分布进行综合评价。 
对于本领域的普通技术人员而言，具体实施例只是结合附图对本发明进行了示例性描述，显然本发明具体实现并不受上述方式的限制，只要采用了本发明的方法构思和技术方案进行的各种非实质性的改进，或未经改进将本发明的构思和技术方案直接应用于其它场合的，均在本发明的保护范围之内。 

Claims (2)

1.一种零值法单砂体逐层剥离预测方法，其特征在于包括以下步骤：
(1)离散零值点求取层位标定：根据已知井测井资料建立砂泥岩薄互层速度模型，并进行离散合成记录标定分析，计算每个单砂体反射波的所有零值点；
(2)分尺度地震精细解释：根据合成记录标定的目的层上下宏观范围内可等时追踪的大尺度等效地层作为约束背景；按照目的层上下宏观范围内可等时追踪的最小等效地层解释范围锁定目标；再对包含目的层的一个波长范围内的等效地层进行小尺度精细层位解释，追踪出复合波的零值点、波峰、波谷层位；
(3)分尺度古地形恢复：针对目的层分别进行大尺度、中尺度、小尺度古地形恢复，锁定一个各尺度古地形背景较一致的单元范围进行简化求解；
(4)离散零值点指导下的地震属性切片：某一单砂体的离散合成记录零值点处的地震属性切片，该砂体无反射记录贡献，因此反映的只是其它砂体的分布特征，即可剥离掉该砂层对地震复合波的影响；
其中，所述大尺度是指包含目的层的三级层序单元顶底界面等时范围；中尺度是指包含目的层的两个波长范围；小尺度是指包含目的层的一个波长范围。
2.根据权利要求1所述的零值法单砂体逐层剥离预测方法，其特征在于在步骤(2)中，通过对测井曲线求取反射系数，并根据阀值做弱反射系数过滤后，计算出主要砂层的反射系数界面，并依此计算各单砂体的离散反射记录及多砂体复合反射合成记录，并给出各单砂体的合成记录零值点时间。
CN201310417935.5A 2013-09-16 2013-09-16 一种零值法单砂体逐层剥离预测方法 Active CN103439741B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310417935.5A CN103439741B (zh) 2013-09-16 2013-09-16 一种零值法单砂体逐层剥离预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310417935.5A CN103439741B (zh) 2013-09-16 2013-09-16 一种零值法单砂体逐层剥离预测方法

Publications (2)

Publication Number Publication Date
CN103439741A true CN103439741A (zh) 2013-12-11
CN103439741B CN103439741B (zh) 2015-10-14

Family

ID=49693445

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310417935.5A Active CN103439741B (zh) 2013-09-16 2013-09-16 一种零值法单砂体逐层剥离预测方法

Country Status (1)

Country Link
CN (1) CN103439741B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103969685A (zh) * 2014-04-23 2014-08-06 长江大学 一种薄互层地震信号的处理方法
CN104991275A (zh) * 2015-07-09 2015-10-21 北京诺克斯达石油科技有限公司 一种特征切片薄互层分析法
CN106405637A (zh) * 2015-07-28 2017-02-15 中国石油化工股份有限公司 一种用于砂泥岩薄互层砂体厚度预测的方法及系统
CN107025360A (zh) * 2017-04-21 2017-08-08 中国石油大港油田勘探开发研究院 一种利用测井资料恢复古地层压力方法
CN106526670B (zh) * 2016-09-21 2019-04-19 中石化石油工程技术服务有限公司 一种碎屑岩储层中地震属性砂体空间分布描述及评价的方法
CN111090117A (zh) * 2018-10-24 2020-05-01 中国石油化工股份有限公司 一种相控正演约束下的有效储层预测方法及系统
CN111580156A (zh) * 2019-02-18 2020-08-25 中国石油天然气股份有限公司 一种地震零值切片自动拾取方法及系统
CN111812743A (zh) * 2019-04-10 2020-10-23 中国石油化工股份有限公司 一种储层单砂体的识别方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5583825A (en) * 1994-09-02 1996-12-10 Exxon Production Research Company Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data
WO2012118864A2 (en) * 2011-02-28 2012-09-07 Schlumberger Technology Corporation Multiscale digital rock modeling for reservoir simulation
CN102879822A (zh) * 2012-09-28 2013-01-16 电子科技大学 一种基于Contourlet变换的地震多属性融合方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5583825A (en) * 1994-09-02 1996-12-10 Exxon Production Research Company Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data
WO2012118864A2 (en) * 2011-02-28 2012-09-07 Schlumberger Technology Corporation Multiscale digital rock modeling for reservoir simulation
CN102879822A (zh) * 2012-09-28 2013-01-16 电子科技大学 一种基于Contourlet变换的地震多属性融合方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
TUCKER F ET AL.: "High-frequency Miocene sequence stratigraphy, offshore Louisiana: Cycle framework and influence on production distribution in a mature shelf province", 《AAPG BULLETIN》, vol. 87, no. 2, 28 February 2003 (2003-02-28), pages 197 - 230 *
冯延状，等: "利用地震资料进行古地形恢复方法研究及应用", 《地层学杂志》, vol. 31, 31 December 2007 (2007-12-31), pages 527 - 530 *
岳英，等: "微地震相三维模型正演储层预测方法与应用", 《中国石油和化工标准与质量》, 10 July 2013 (2013-07-10), pages 170 - 171 *
李国发，等: "基于模型的薄互层地震属性分析及其应用", 《石油物探》, vol. 50, no. 2, 31 March 2011 (2011-03-31), pages 144 - 149 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103969685A (zh) * 2014-04-23 2014-08-06 长江大学 一种薄互层地震信号的处理方法
CN103969685B (zh) * 2014-04-23 2016-07-13 长江大学 一种薄互层地震信号的处理方法
CN104991275A (zh) * 2015-07-09 2015-10-21 北京诺克斯达石油科技有限公司 一种特征切片薄互层分析法
CN104991275B (zh) * 2015-07-09 2018-03-09 北京诺克斯达石油科技有限公司 一种特征切片薄互层分析法
CN106405637A (zh) * 2015-07-28 2017-02-15 中国石油化工股份有限公司 一种用于砂泥岩薄互层砂体厚度预测的方法及系统
CN106405637B (zh) * 2015-07-28 2018-06-12 中国石油化工股份有限公司 一种用于砂泥岩薄互层砂体厚度预测的方法及系统
CN106526670B (zh) * 2016-09-21 2019-04-19 中石化石油工程技术服务有限公司 一种碎屑岩储层中地震属性砂体空间分布描述及评价的方法
CN107025360A (zh) * 2017-04-21 2017-08-08 中国石油大港油田勘探开发研究院 一种利用测井资料恢复古地层压力方法
CN111090117A (zh) * 2018-10-24 2020-05-01 中国石油化工股份有限公司 一种相控正演约束下的有效储层预测方法及系统
CN111580156A (zh) * 2019-02-18 2020-08-25 中国石油天然气股份有限公司 一种地震零值切片自动拾取方法及系统
CN111812743A (zh) * 2019-04-10 2020-10-23 中国石油化工股份有限公司 一种储层单砂体的识别方法
CN111812743B (zh) * 2019-04-10 2022-02-11 中国石油化工股份有限公司 一种储层单砂体的识别方法

Also Published As

Publication number Publication date
CN103439741B (zh) 2015-10-14

Similar Documents

Publication Publication Date Title
CN103439741B (zh) 一种零值法单砂体逐层剥离预测方法
Anees et al. Channel identification using 3D seismic attributes and well logging in lower Shihezi Formation of Hangjinqi area, northern Ordos Basin, China
CA2920499C (en) Stratigraphic function
Pugin et al. Near-surface mapping using SH-wave and P-wave seismic land-streamer data acquisition in Illinois, US
Norden et al. Geological modelling of the Triassic Stuttgart Formation at the Ketzin CO2 storage site, Germany
CN102466815B (zh) 三叠系碎屑岩油气藏识别方法
Abdel-Fattah et al. 3D geometric modeling of the Abu Madi reservoirs and its implication on the gas development in Baltim area (Offshore Nile Delta, Egypt)
EP2994781A1 (en) Estimation of q-factor in time domain
Bujakowski et al. Integrated seismic and magnetotelluric exploration of the Skierniewice, Poland, geothermal test site
Ebong et al. 3D structural modelling and fluid identification in parts of Niger Delta Basin, southern Nigeria
Qin et al. An interactive integrated interpretation of GPR and Rayleigh wave data based on the genetic algorithm
Alvarez et al. Gravity, seismic reflection, and tomographic constraints on the subduction-to-strike-slip transition at the southeastern Caribbean Plate boundary zone
Kim et al. Lithofacies modeling by multipoint statistics and economic evaluation by NPV volume for the early Cretaceous Wabiskaw Member in Athabasca oilsands area, Canada
Zhou et al. Seismic diffraction imaging for improved coal structure detection by a structure-oriented moving average error filter
CN105954800A (zh) 利用地震纹理特性识别小断层的方法
Khan et al. Three-dimensional structural modeling (3D SM) and joint geophysical characterization (JGC) of hydrocarbon reservoir: a case study of the Kadanwari field in Middle Indus Basin (MIB), Southeastern Pakistan
Bakke et al. Seismic modeling of outcrop analogues: Techniques and applications
Deng et al. Hydrocarbon accumulation conditions and key exploration and development technologies for PL 19–3 oilfield
Murty et al. Mapping the sub-trappean Mesozoic sediments in the western part of Narmada–Tapti region of Deccan Volcanic Province, India
Zhao et al. Tidal sand ridges seismic identification and its application in the Xihu Depression, East China Sea Shelf Basin: enlightenment to hydrocarbon exploration
Chongwain et al. Seismic attribute analysis for reservoir description and characterization of M-field, Douala Sub-Basin, Cameroon
Werner et al. Chapter 11 Structural interpretation of the Barents and Kara Seas from gravity and magnetic data
Kempka 3D structural geological model of the Mallik Anticline, Canadian Arctic
Shillington et al. P‐and S‐wave velocities of consolidated sediments from a seafloor seismic survey in the North Celtic Sea Basin, offshore Ireland
Zhang et al. The velocity modeling techniques of igneous rock—A case study of Ha-La-Ha-Tang area

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230821

Address after: 100007 No. 9 North Main Street, Dongcheng District, Beijing, Dongzhimen

Patentee after: PETROCHINA Co.,Ltd.

Address before: 300280 Exploration and Development Research Institute of Xingfu Avenue, Dagang Oilfield, Dagang District, Tianjin City

Patentee before: DAGANG OIL FIELD OF CNPC