CN108646294A - 一种复杂断裂系统下深水浊积砂岩储层高精度反演方法 - Google Patents

一种复杂断裂系统下深水浊积砂岩储层高精度反演方法 Download PDF

Info

Publication number
CN108646294A
CN108646294A CN201810694506.5A CN201810694506A CN108646294A CN 108646294 A CN108646294 A CN 108646294A CN 201810694506 A CN201810694506 A CN 201810694506A CN 108646294 A CN108646294 A CN 108646294A
Authority
CN
China
Prior art keywords
cut
layer
inversion
tomography
model
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
CN201810694506.5A
Other languages
English (en)
Other versions
CN108646294B (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.)
Beijing Research Center of CNOOC China Ltd
CNOOC China Ltd
Original Assignee
Beijing Research Center of CNOOC China Ltd
CNOOC China Ltd
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 Beijing Research Center of CNOOC China Ltd, CNOOC China Ltd filed Critical Beijing Research Center of CNOOC China Ltd
Priority to CN201810694506.5A priority Critical patent/CN108646294B/zh
Publication of CN108646294A publication Critical patent/CN108646294A/zh
Application granted granted Critical
Publication of CN108646294B publication Critical patent/CN108646294B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/51Migration
    • G01V2210/512Pre-stack
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/624Reservoir parameters

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种复杂断裂系统下深水浊积砂岩储层高精度反演方法,包括以下步骤:1)根据断层解释数据建立三维断层框架模型;2)根据地震层位解释数据建立三维层位框架模型;3)根据构造演化特征及地质先验信息判断断层与层位的发育先后,以此为依据定义断层与层位的空间削截关系,将步骤1)得到的三维断层框架模型与步骤2)得到的三维层位框架模型组合得到高精度的三维地震地层格架;4)采用“多体相控迭代法”构建反演低频模型;5)迭代优化叠前反演过程,直至反演结果满足空间误差和损失函数最小时为止,输出最终反演结果。

Description

一种复杂断裂系统下深水浊积砂岩储层高精度反演方法
技术领域
本发明涉及一种高精度储层反演方法,尤其涉及一种复杂断裂系统下深水浊积砂岩储层高精度反演方法。
背景技术
深水浊流沉积储层油田已成为当今世界石油储量和产量的增长来源,该类储层为典型的重力流沉积,储层横向突变,纵向多期叠置,迁移摆动频繁,储层精细刻画难度大。再加上构造上断裂系统发育,多期断层相互切割,组合关系复杂,如何对这种复杂断裂系统下的深水浊积砂岩储层进行精细描述是保证后续油田开发生产的关键环节。
由于重力流沉积的突发性及与此相伴生的盐底辟、泥底辟、逆冲背斜等构造响应,造成这类储集体的空间展布、岩性、物性及含油气性等特征表现得较为复杂,叠后反演往往难以刻画此类储层展布特征,在生产实践中,叠前反演逐步成为储层精细描述的关键技术。深水浊流储层的这种复杂性对油藏描述的精细程度要求较高,对叠前反演技术的要求也相对提高。目前,有部分学者开展了针对深水浊积储层的叠前反演研究。例如,闵小刚等基于AVO正演模型对近道多次波进行压制,进而提取PG属性,结合岩石物理分析,对工区储层的含油性、横向连通性、厚度、孔渗等特征取得了较好的认识;黄饶等首先应用叠前AVA同时反演技术提取纵横波阻抗、纵横波速度比等弹性参数,然后利用属性转换技术求取泊松比和烃类指示属性等转换参数,较好地预测了储层岩性及含油气性分布规律;解吉高等用纵波/横波速度比属性预测具有浊积水道展布形态的地震异常体的岩性,在砂岩储集体刻画基础上,用叠前重构的泊松阻抗技术有效识别砂岩孔隙内的流体类型等。但对于复杂断裂系统下的深水浊积砂岩储层,以上技术没有提供完善的解决方案。
发明内容
针对上述问题,本发明的目的是提供一种复杂断裂系统下深水浊积砂岩储层高精度反演方法。
为实现上述目的,本发明采取以下技术方案:一种复杂断裂系统下深水浊积砂岩储层高精度反演方法,包括以下步骤:1)根据断层解释数据建立三维断层框架模型;2)根据地震层位解释数据建立三维层位框架模型;3)根据构造演化特征及地质先验信息判断断层与层位的发育先后,以此为依据定义断层与层位的空间削截关系,将步骤1)得到的三维断层框架模型与步骤2)得到的三维层位框架模型组合得到高精度的三维地震地层格架;4)采用“多体相控迭代法”构建反演低频模型;5)迭代优化叠前反演过程,直至反演结果满足空间误差和损失函数最小时为止,输出最终反演结果。
所述步骤1)中,对共轭断层采用削截关系拆分-组合的定义方法,具体过程如下:①根据共轭断层的接触范围及空间接触关系,选定晚期发育断层为被削截断层,早期发育断层为削截断层;②将被削截断层拆分为三个断层,即削截断层与被削截断层完全未接触部分定义为一条断层,削截断层与被削截断层接触部分定义为两条断层;③三条拆分后断层分别与早期发育的削截断层按照不同的断层削截方向定义接触关系;④将定义好接触关系的断层重新融合为晚期发育的被削截断层,完成共轭断层削截关系定义。
所述步骤4)中的反演低频模型包括纵横波速度模型和密度模型。
所述步骤4)中,采用“多体相控迭代法”构建纵横波速度、密度等反演低频模型的过程如下:①根据密度测井曲线的低频趋势线采用马尔科夫链-蒙特卡洛算法计算得到初始的密度体数据;②在地震偏移速度体数据基础上,利用井点位置处声波测井曲线并采用马尔科夫链-蒙特卡洛算法得到校正后的偏移速度数据;③利用Gassmann方程计算横波速度数据体;④将三维地震地层格架内充填上述步骤计算的速度、密度,得到用于叠前反演的初始低频模型;⑤利用初始低频模型进行叠前反演得到初次反演结果,对初次反演结果提取能反映储层展布的敏感地震体属性,并结合先期地质认识修正地震属性对储层展布不符合的区域,通过合理添加虚拟井的方式使得地震属性反映的储层展布能够符合地质认识,得到地震相数据体;⑥利用地震相数据体修改参与反演的低频模型。
本发明由于采取以上技术方案,其具有以下优点:1、本发明适用于复杂断裂系统下的深水浊积砂岩储层,能够用于预测储层岩性及含油气性分布规律。2、本发明通过创新“共轭断层削截关系拆分-组合定义法”,能够准确构建共轭断层、花状断层等复杂断层大量发育时的三维断层框架模型,提高断层附近储层的预测精度,为断层附近井位设计提供可靠依据。3、通过创新“多体相控迭代法”构建反演低频模型,相比传统叠前反演方法能够提高反演结果对储层横向边界的刻画精度,为浊流储层突变边界的识别与刻画提供重要方法。
附图说明
图1是本发明方法的流程示意图;
图2是本发明创新发明的共轭断层削截关系拆分-组合定义法示意图;
图3是本发明建立的三维断层框架模型图;
图4是本发明建立的三维层位框架模型;
图5是本发明添加虚拟井得到的符合储层展布地质认识的地震属性图;
图6(a)是实际目的层段的原始地震剖面图;
图6(b)是采用本发明方法得到的叠前反演剖面图;
图7(a)是本发明得到的反演结果提取的均方根振幅属性;
图7(b)是常规叠前反演结果提取的均方根振幅属性。
具体实施方式
下面结合附图和实施案例对本发明进行详细的描述。
如图1所示,本发明提出了一种复杂断裂系统下深水浊积砂岩储层高精度反演方法,包括以下步骤:
1)根据断层解释数据建立三维断层框架模型。
断层组合样式多样,包括平行、亚平行、Y字形、反Y字形、入字形及共轭断层,其中,共轭断层削截关系定义较为困难,本发明对共轭断层采用削截关系拆分-组合的定义方法,具体过程如下(如图2所示):
①根据共轭断层的接触范围及空间接触关系,选定晚期发育断层为被削截断层,早期发育断层为削截断层;
②将被削截断层拆分为三个断层,即削截断层与被削截断层完全未接触部分定义为一条断层,削截断层与被削截断层接触部分定义为两条断层(相同断层数据,不同断层名称标记);
③三条拆分后断层分别与早期发育的削截断层按照不同的断层削截方向定义接触关系;
④将定义好接触关系的断层重新融合为晚期发育的被削截断层,完成共轭断层削截关系定义。
图3为建立的三维断层框架模型图。
2)根据地震层位解释数据建立三维层位框架模型(如图4所示)。
3)根据构造演化特征及地质先验信息判断断层与层位的发育先后,以此为依据定义断层与层位的空间削截关系,将步骤1)得到的三维断层框架模型与步骤2)得到的三维层位框架模型组合得到高精度的三维地震地层格架。
4)采用“多体相控迭代法”构建纵横波速度、密度等反演低频模型,其具体步骤为:
①根据密度测井曲线的低频趋势线,采用马尔科夫链-蒙特卡洛算法计算得到初始的密度体数据;
②在地震偏移速度体数据基础上,利用井点位置处声波测井曲线并采用马尔科夫链-蒙特卡洛算法得到校正后的偏移速度数据;
③利用Gassmann方程计算横波速度数据体;
④将三维地震地层格架内充填上述步骤计算的速度、密度,得到用于叠前反演的初始低频模型;
⑤利用初始低频模型进行叠前反演得到初次反演结果,对初次反演结果提取能反映储层展布的敏感地震体属性,并结合先期地质认识修正地震属性对储层展布不符合的区域,通过合理添加虚拟井的方式使得地震属性反映的储层展布能够符合地质认识(如图5所示),得到地震相数据体;
⑥利用地震相数据体修改参与反演的低频模型。
5)迭代优化叠前反演过程,直至反演结果满足空间误差和损失函数最小时为止,输出最终反演结果。
图6、图7可以说明本发明所述的反演方法的技术效果:图6(a)为实际目的层段的原始地震剖面图;如图6(b)为采用本发明方法得到的叠前反演剖面图,可以看出储层边界特征清晰,与原始地震剖面特征吻合。图7(a)所示为本发明得到的反演结果提取的均方根振幅属性,图7(b)所示为常规叠前反演结果提取的均方根振幅属性。对比可以看到,本发明得到的反演结果浊积水道特征明显,较常规叠前反演结果而言,对储层的边界刻画程度大大提升。
上述各实施例仅用于对本发明的目的、技术方案和有益效果进行了进一步详细说明,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种复杂断裂系统下深水浊积砂岩储层高精度反演方法,包括以下步骤:
1)根据断层解释数据建立三维断层框架模型;
2)根据地震层位解释数据建立三维层位框架模型;
3)根据构造演化特征及地质先验信息判断断层与层位的发育先后,以此为依据定义断层与层位的空间削截关系,将步骤1)得到的三维断层框架模型与步骤2)得到的三维层位框架模型组合得到高精度的三维地震地层格架;
4)采用“多体相控迭代法”构建反演低频模型;
5)迭代优化叠前反演过程,直至反演结果满足空间误差和损失函数最小时为止,输出最终反演结果。
2.如权利要求1所述的一种复杂断裂系统下深水浊积砂岩储层高精度反演方法,其特征在于:所述步骤1)中,对共轭断层采用削截关系拆分-组合的定义方法,具体过程如下:
①根据共轭断层的接触范围及空间接触关系,选定晚期发育断层为被削截断层,早期发育断层为削截断层;
②将被削截断层拆分为三个断层,即削截断层与被削截断层完全未接触部分定义为一条断层,削截断层与被削截断层接触部分定义为两条断层;
③三条拆分后断层分别与早期发育的削截断层按照不同的断层削截方向定义接触关系;
④将定义好接触关系的断层重新融合为晚期发育的被削截断层,完成共轭断层削截关系定义。
3.如权利要求1所述的一种复杂断裂系统下深水浊积砂岩储层高精度反演方法,其特征在于:所述步骤4)中的反演低频模型包括纵横波速度模型和密度模型。
4.如权利要求3所述的一种复杂断裂系统下深水浊积砂岩储层高精度反演方法,其特征在于:所述步骤4)中,采用“多体相控迭代法”构建纵横波速度、密度等反演低频模型的过程如下:
①根据密度测井曲线的低频趋势线采用马尔科夫链-蒙特卡洛算法计算得到初始的密度体数据;
②在地震偏移速度体数据基础上,利用井点位置处声波测井曲线并采用马尔科夫链-蒙特卡洛算法得到校正后的偏移速度数据;
③利用Gassmann方程计算横波速度数据体;
④将三维地震地层格架内充填上述步骤计算的速度、密度,得到用于叠前反演的初始低频模型;
⑤利用初始低频模型进行叠前反演得到初次反演结果,对初次反演结果提取能反映储层展布的敏感地震体属性,并结合先期地质认识修正地震属性对储层展布不符合的区域,通过合理添加虚拟井的方式使得地震属性反映的储层展布能够符合地质认识,得到地震相数据体;
⑥利用地震相数据体修改参与反演的低频模型。
CN201810694506.5A 2018-06-29 2018-06-29 一种复杂断裂系统下深水浊积砂岩储层高精度反演方法 Active CN108646294B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810694506.5A CN108646294B (zh) 2018-06-29 2018-06-29 一种复杂断裂系统下深水浊积砂岩储层高精度反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810694506.5A CN108646294B (zh) 2018-06-29 2018-06-29 一种复杂断裂系统下深水浊积砂岩储层高精度反演方法

Publications (2)

Publication Number Publication Date
CN108646294A true CN108646294A (zh) 2018-10-12
CN108646294B CN108646294B (zh) 2019-09-06

Family

ID=63749975

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810694506.5A Active CN108646294B (zh) 2018-06-29 2018-06-29 一种复杂断裂系统下深水浊积砂岩储层高精度反演方法

Country Status (1)

Country Link
CN (1) CN108646294B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112230301A (zh) * 2020-09-18 2021-01-15 西南石油大学 一种深水水道成因类型划分方法
CN112327365A (zh) * 2020-11-06 2021-02-05 大庆油田有限责任公司 一种复杂断块构造建模方法
CN113933897A (zh) * 2020-07-13 2022-01-14 中国石油天然气股份有限公司 一种基于二维叠前资料的膏岩尖灭线预测方法及系统

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101980053A (zh) * 2010-10-18 2011-02-23 中国石油化工股份有限公司 一种复杂礁滩储层预测方法
CN102096107A (zh) * 2009-12-09 2011-06-15 中国石油天然气股份有限公司 一种根据声波时差和密度反演孔隙扁度进行储层渗透性评价的方法
CN102353989A (zh) * 2011-08-24 2012-02-15 成都理工大学 基于自适应基质矿物等效弹性模量反演的横波速度估算方法
WO2014028030A1 (en) * 2012-08-17 2014-02-20 Landmark Graphics Corporation Systems and methods for imaging seismic data
CN104200115A (zh) * 2014-09-12 2014-12-10 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种基于地质统计学模拟的全地层速度建模方法
CN104297800A (zh) * 2014-10-09 2015-01-21 中国海洋石油总公司 一种自相控叠前反演方法
CN104360385A (zh) * 2014-11-26 2015-02-18 中国石油集团川庆钻探工程有限公司地球物理勘探公司 结合地震地质认识综合建立初始深度层速度模型的方法
CN105182444A (zh) * 2015-09-09 2015-12-23 中国海洋石油总公司 一种高分辨率层序地层格架约束的地质统计学反演方法
CN105652323A (zh) * 2015-08-11 2016-06-08 中国石油化工股份有限公司 一种储层预测方法
CN106094032A (zh) * 2016-08-30 2016-11-09 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种构建地层速度模型的方法
CN106932819A (zh) * 2017-02-23 2017-07-07 河海大学 基于各向异性马尔科夫随机域的叠前地震参数反演方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102096107A (zh) * 2009-12-09 2011-06-15 中国石油天然气股份有限公司 一种根据声波时差和密度反演孔隙扁度进行储层渗透性评价的方法
CN101980053A (zh) * 2010-10-18 2011-02-23 中国石油化工股份有限公司 一种复杂礁滩储层预测方法
CN102353989A (zh) * 2011-08-24 2012-02-15 成都理工大学 基于自适应基质矿物等效弹性模量反演的横波速度估算方法
WO2014028030A1 (en) * 2012-08-17 2014-02-20 Landmark Graphics Corporation Systems and methods for imaging seismic data
CN104200115A (zh) * 2014-09-12 2014-12-10 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种基于地质统计学模拟的全地层速度建模方法
CN104297800A (zh) * 2014-10-09 2015-01-21 中国海洋石油总公司 一种自相控叠前反演方法
CN104360385A (zh) * 2014-11-26 2015-02-18 中国石油集团川庆钻探工程有限公司地球物理勘探公司 结合地震地质认识综合建立初始深度层速度模型的方法
CN105652323A (zh) * 2015-08-11 2016-06-08 中国石油化工股份有限公司 一种储层预测方法
CN105182444A (zh) * 2015-09-09 2015-12-23 中国海洋石油总公司 一种高分辨率层序地层格架约束的地质统计学反演方法
CN106094032A (zh) * 2016-08-30 2016-11-09 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种构建地层速度模型的方法
CN106932819A (zh) * 2017-02-23 2017-07-07 河海大学 基于各向异性马尔科夫随机域的叠前地震参数反演方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
曹彤 等: ""高分辨率层序格架约束的地震反演及应用"", 《科技导报》 *
马良涛 等: ""不确定性反演关键参数的地质含义及正演模型反演研究——以渤海海域W油田为例"", 《地球物理学进展》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113933897A (zh) * 2020-07-13 2022-01-14 中国石油天然气股份有限公司 一种基于二维叠前资料的膏岩尖灭线预测方法及系统
CN113933897B (zh) * 2020-07-13 2024-03-01 中国石油天然气股份有限公司 一种基于二维叠前资料的膏岩尖灭线预测方法及系统
CN112230301A (zh) * 2020-09-18 2021-01-15 西南石油大学 一种深水水道成因类型划分方法
CN112327365A (zh) * 2020-11-06 2021-02-05 大庆油田有限责任公司 一种复杂断块构造建模方法
CN112327365B (zh) * 2020-11-06 2022-07-26 大庆油田有限责任公司 一种复杂断块构造建模方法

Also Published As

Publication number Publication date
CN108646294B (zh) 2019-09-06

Similar Documents

Publication Publication Date Title
CN107817535B (zh) 薄夹层的确定方法和装置
CN102867330B (zh) 基于区域划分的空间复杂层位重构方法
CN105652323B (zh) 一种储层预测方法
CN108646294B (zh) 一种复杂断裂系统下深水浊积砂岩储层高精度反演方法
Toro et al. Thermal and kinematic evolution of the Eastern Cordillera fold and thrust belt, Colombia
CN110579802B (zh) 一种天然气水合物储层物性参数的高精度反演方法
CN108363100A (zh) 基于封存单元与粗糙集的煤层气甜点区地震地质识别方法
CN103901478A (zh) 一种井震信息联合确定储层沉积特征和分布的方法
Tavani et al. Transpressive inversion of a Mesozoic extensional forced fold system with an intermediate décollement level in the Basque‐Cantabrian Basin (Spain)
CN104635261B (zh) 山前带重磁电震联合分体建模方法
CN109541685B (zh) 一种河道砂体识别方法
CN112394404B (zh) 一种渐进式储层精细表征方法
CN104749617A (zh) 一种多尺度裂缝储层正演模型建立方法
CN105607120A (zh) 一种基于时移测井的地震相约束构建初始模型的方法
CN105089652A (zh) 一种拟声波曲线重构与稀疏脉冲联合反演方法
CN105005079A (zh) 一种井曲线反演方法
CN104297800B (zh) 一种自相控叠前反演方法
CN107942379A (zh) 一种提高复杂断块速度模型精度的方法
Taixian et al. Techniques for high-efficient development of offshore fluvial oilfields
Tao et al. Geological structure and dynamic mechanism of the Termit rift basin in West African rift system
Zhang et al. Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China
CN112649893B (zh) 一种面向薄储层的多资料多参数融合建模方法及系统
Faidhllah et al. Building of a Static Model for Jeribe Formation in Jambour Oil Field
CN107991702B (zh) 正断层时序的简单精准判定方法
Shuwen et al. Multiple facies-controlled prediction of single sand body of meandering rivers

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant