CN103439741B - 一种零值法单砂体逐层剥离预测方法 - Google Patents
一种零值法单砂体逐层剥离预测方法 Download PDFInfo
- Publication number
- CN103439741B CN103439741B CN201310417935.5A CN201310417935A CN103439741B CN 103439741 B CN103439741 B CN 103439741B CN 201310417935 A CN201310417935 A CN 201310417935A CN 103439741 B CN103439741 B CN 103439741B
- Authority
- CN
- China
- Prior art keywords
- sand body
- zero point
- seismic
- interest
- zone
- 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.)
- Active
Links
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所述的零值法单砂体逐层剥离预测方法,其特征在于在步骤(1)中,通过对测井曲线求取反射系数,并根据阀值做弱反射系数过滤后,计算出主要砂层的反射系数界面,并依此计算各单砂体的离散反射记录及多砂体复合反射合成记录,并给出各单砂体的合成记录零值点时间。
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 CN103439741A (zh) | 2013-12-11 |
CN103439741B true 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) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103969685B (zh) * | 2014-04-23 | 2016-07-13 | 长江大学 | 一种薄互层地震信号的处理方法 |
CN104991275B (zh) * | 2015-07-09 | 2018-03-09 | 北京诺克斯达石油科技有限公司 | 一种特征切片薄互层分析法 |
CN106405637B (zh) * | 2015-07-28 | 2018-06-12 | 中国石油化工股份有限公司 | 一种用于砂泥岩薄互层砂体厚度预测的方法及系统 |
CN106526670B (zh) * | 2016-09-21 | 2019-04-19 | 中石化石油工程技术服务有限公司 | 一种碎屑岩储层中地震属性砂体空间分布描述及评价的方法 |
CN107025360A (zh) * | 2017-04-21 | 2017-08-08 | 中国石油大港油田勘探开发研究院 | 一种利用测井资料恢复古地层压力方法 |
CN111090117B (zh) * | 2018-10-24 | 2021-11-19 | 中国石油化工股份有限公司 | 一种相控正演约束下的有效储层预测方法及系统 |
CN111580156B (zh) * | 2019-02-18 | 2022-12-02 | 中国石油天然气股份有限公司 | 一种地震零值切片自动拾取方法及系统 |
CN111812743B (zh) * | 2019-04-10 | 2022-02-11 | 中国石油化工股份有限公司 | 一种储层单砂体的识别方法 |
Citations (2)
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 |
CN102879822A (zh) * | 2012-09-28 | 2013-01-16 | 电子科技大学 | 一种基于Contourlet变换的地震多属性融合方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9134457B2 (en) * | 2009-04-08 | 2015-09-15 | Schlumberger Technology Corporation | Multiscale digital rock modeling for reservoir simulation |
-
2013
- 2013-09-16 CN CN201310417935.5A patent/CN103439741B/zh active Active
Patent Citations (2)
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 |
CN102879822A (zh) * | 2012-09-28 | 2013-01-16 | 电子科技大学 | 一种基于Contourlet变换的地震多属性融合方法 |
Non-Patent Citations (4)
Title |
---|
High-frequency Miocene sequence stratigraphy, offshore Louisiana: Cycle framework and influence on production distribution in a mature shelf province;Tucker F et al.;《AAPG Bulletin》;20030228;第87卷(第2期);第197-230页 * |
冯延状,等.利用地震资料进行古地形恢复方法研究及应用.《地层学杂志》.2007,第31卷第527-530页. * |
岳英,等.微地震相三维模型正演储层预测方法与应用.《中国石油和化工标准与质量》.2013,第170-171页. * |
李国发,等.基于模型的薄互层地震属性分析及其应用.《石油物探》.2011,第50卷(第2期),第144-149页. * |
Also Published As
Publication number | Publication date |
---|---|
CN103439741A (zh) | 2013-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103439741B (zh) | 一种零值法单砂体逐层剥离预测方法 | |
Alkhalifah | Velocity analysis using nonhyperbolic moveout in transversely isotropic media | |
Abdel-Fattah et al. | Integrated seismic attributes and stochastic inversion for reservoir characterization: insights from Wadi field (NE Abu-Gharadig Basin, Egypt) | |
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) | |
Bujakowski et al. | Integrated seismic and magnetotelluric exploration of the Skierniewice, Poland, geothermal test site | |
CN106405637A (zh) | 一种用于砂泥岩薄互层砂体厚度预测的方法及系统 | |
Ebong et al. | 3D structural modelling and fluid identification in parts of Niger Delta Basin, southern Nigeria | |
Yang et al. | Moving source profile data processing, modelling and comparison with 3D surface seismic data at the CO2SINK project site, Ketzin, Germany | |
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 | |
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 | |
Zeng et al. | Multistage progradational clinoform-set characterisation and evolution analysis of the Early Oligocene in the Baiyun Sag, Pearl River Mouth Basin, South China Sea | |
CN105954800A (zh) | 利用地震纹理特性识别小断层的方法 | |
Deng et al. | Hydrocarbon accumulation conditions and key exploration and development technologies for PL 19–3 oilfield | |
Zhenwu et al. | Current situation and trend of geophysical technology in CNPC | |
Murty et al. | Mapping the sub-trappean Mesozoic sediments in the western part of Narmada–Tapti region of Deccan Volcanic Province, India | |
Zhang et al. | The velocity modeling techniques of igneous rock—A case study of Ha-La-Ha-Tang area | |
Shillington et al. | P‐and S‐wave velocities of consolidated sediments from a seafloor seismic survey in the North Celtic Sea Basin, offshore Ireland | |
Hu et al. | Horizontal well path design with seismic inversion for a shale oil reservoir in the Ordos Basin, China | |
Gritto et al. | Estimating subsurface topography from surface-to-borehole seismic studies at the Rye Patch geothermal reservoir, Nevada, USA | |
Masinde et al. | Spatial and temporal control of the volcanic activities over the distribution of sediments in the Magadi basin, Kenya Rift | |
Laurenzano et al. | Numerical modelling of earthquake strong ground motion in the area of Vittorio Veneto (NE Italy) | |
Guan et al. | Application of Facies Controlled Prestack Inversion Technology in Prediction of Deep Tight Gas in Continental Faulted Basin | |
Shelokhov et al. | Shallow velocity model from the transient electromagnetic method data: results of application in Eastern and Western Siberia | |
Ata et al. | The improvements of three-dimensional seismic interpretation in comparison with the two-dimensional seismic interpretation in Al-Amal oil field, Gulf of Suez, Egypt |
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 |
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 |
|
TR01 | Transfer of patent right |