CN113433589A - 一种基于数理统计的风化壳储层底界面识别方法 - Google Patents

一种基于数理统计的风化壳储层底界面识别方法 Download PDF

Info

Publication number
CN113433589A
CN113433589A CN202110886389.4A CN202110886389A CN113433589A CN 113433589 A CN113433589 A CN 113433589A CN 202110886389 A CN202110886389 A CN 202110886389A CN 113433589 A CN113433589 A CN 113433589A
Authority
CN
China
Prior art keywords
reservoir
weathering crust
weathering
crust
ancient
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.)
Pending
Application number
CN202110886389.4A
Other languages
English (en)
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.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
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 Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN202110886389.4A priority Critical patent/CN113433589A/zh
Publication of CN113433589A publication Critical patent/CN113433589A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/624—Reservoir 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

本发明涉及储层底界面识别技术领域，具体涉及一种基于数理统计的风化壳储层底界面识别方法；包括以下步骤：S1、以井资料为依据构建岩石物理模型，S2、确定不同厚度花岗岩风化壳储层与地震响应之间的联系，S3、通过印模法研究确定花岗岩潜山古地貌，并在古地貌图上选取不同暴露时间点，统计不同暴露时间点在地震剖面上对应的风化壳储集层顶底界面相对时间差，以正演模拟为基准进行多样点拟合分析，获得相应的定量关系式，S4、以全区风化壳储层顶界面为输入获得全区风化壳储层底界面标定结果，以及风化壳储层时间厚度定量结果图；本发明的方法可行性和可靠性极佳，为花岗岩潜山风化壳的圈定工作提供了一个良好的思路和方法。

Description

一种基于数理统计的风化壳储层底界面识别方法
技术领域
本发明涉及储层底界面识别技术领域，具体涉及一种基于数理统计的风化壳储层底界面识别方法。
背景技术
古潜山油气储层具有很好的开发前景与巨大的勘探开发潜力，在我国渤海湾盆地、东海盆地、珠江口盆地、琼东南盆地、南海北部湾盆地、辽河盆地等地区都有广泛分布。琼东南盆地深水区松南低凸起及其周缘普遍发育前古近系花岗岩风化壳，钻探结果表示该区基底潜山风化壳富集大量油气，但由于基底储层物性横纵向变化较大、底界面地震响应杂乱，储层顶底标定较难，而储层顶底的标定能够直接影响储层预测的精确度，所以针对花岗岩风化壳储层顶底标定的研究格外重要。就目前国内外研究情况来看，在工程勘探方面对基底风化壳的研究较多，而相关于花岗岩风化壳储层的内容较少；胡志伟 (2017)总结出近几年国内渤海盆地蓬莱9-1及柴达木盆地东坪地区在花岗岩基岩风化壳中取得重大突破，奠定花岗岩储层研究的必要性；徐守立等(2019)综合花岗岩潜山的储层特性及控制因素，在纵向上进行了风化壳储层的识别和分带划分；姚元锋(2018)总结了我国花岗岩岩石风化带特征以及风化带划分的一些常用技术和方法，为风化壳纵向分带提供依据；马平(2019)通过对花岗岩风化程度与其所处的微地形相关性进行统计分析，获得地形与风化程度之间的关系。众多学者从纵向上划分了花岗岩风化壳储层，同时提出了一些统计划分的方法，但是少有从横向出发标定研究区内风化壳储层顶底界面方法的研究。
综上所述，研发一种基于数理统计的风化壳储层底界面识别方法，仍是储层底界面识别技术领域中急需解决的关键问题。
发明内容
针对现有技术所存在的上述缺点，本发明的目的在于提供一种基于数理统计的风化壳储层底界面识别方法，本发明从已钻探花岗岩潜山风化壳储集层物性特征出发，进行正演模拟获得相应的地震响应特征，基于古潜山暴露时间与风化程度之间存在的一定关系，通过印模法恢复潜山古地貌，进行多样点统计拟合分析，获得古地貌不同高点与风化壳储层在地震剖面上顶底时差的定量关系式，运用该关系式进行全区风化壳底界面的识别与标定工作，本发明的方法可行性和可靠性极佳，为花岗岩潜山风化壳的圈定工作提供了一个良好的思路和方法。
为实现上述目的，本发明提供了如下技术方案：
一种基于数理统计的风化壳储层底界面识别方法，包括以下步骤：
S1、以井资料为依据构建岩石物理模型。
S2、采用地震正演模拟方法研究确定不同厚度花岗岩风化壳储层与地震响应之间的联系。
S3、从古潜山暴露时间与风化程度关系分析入手，通过印模法研究确定花岗岩潜山古地貌，并在古地貌图上选取不同暴露时间点，统计不同暴露时间点在地震剖面上对应的风化壳储集层顶底界面相对时间差，以正演模拟为基准进行多样点拟合分析，获得相应的定量关系式。
S4、以全区风化壳储层顶界面为输入获得全区风化壳储层底界面标定结果，以及风化壳储层时间厚度定量结果图。
本发明进一步设置为：在所述步骤S1中，所述岩石物理模型为楔形模型，所述楔形模型以强风化带和中风化带作为可成为储集空间的风化壳储集层、以盖层为砂泥岩储集层和以致密的花岗岩为基岩建立。
本发明进一步设置为：在所述步骤S2中，所述地震正演模拟方法是指通过波动方程正演获得叠加剖面，模拟不同厚度的风化壳储集层在地震剖面上的响应。
本发明进一步设置为：在所述步骤S3中，所述印模法是指将前古近系花岗岩潜山顶面结束剥蚀开始上覆沉积时期视为一等时面，根据沉积补偿原理，利用上覆新沉积地层与残余古地貌之间存在镜像关系反映古地貌形态。
本发明进一步设置为：在所述步骤S3中，所述多点拟合分析的方法为在研究区中选取古地貌不同高度特征的10个点，并选取其进行地震剖面反射特征分析。
本发明进一步设置为：在所述步骤S3中，所述定量关系式为： y＝-0.0136x2+1.2088x+11.143，式中，x表示古地貌基准点上的时间高度，y表示该点对应地震剖面上风化壳储集层顶底界面时间差。
本发明进一步设置为：在所述步骤S4中，通过对古地貌风化零界面以上的古地貌时间高度进行提取，作为拟合公式的输入，获得全区风化壳时间厚度分布图。
本发明进一步设置为：在所述步骤S4后，验证风化壳顶底标定结果的准确性。
有益效果
采用本发明提供的技术方案，与已知的公有技术相比，具有如下有益效果：
本发明从已钻探花岗岩潜山风化壳储集层物性特征出发，进行正演模拟获得相应的地震响应特征，基于古潜山暴露时间与风化程度之间存在的一定关系，通过印模法恢复潜山古地貌，进行多样点统计拟合分析，获得古地貌不同高点与风化壳储层在地震剖面上顶底时差的定量关系式，运用该关系式进行全区风化壳底界面的识别与标定工作，本发明的方法可行性和可靠性极佳，为花岗岩潜山风化壳的圈定工作提供了一个良好的思路和方法。
附图说明
图1为一种基于数理统计的风化壳储层底界面识别方法的流程图；
图2a为一种基于数理统计的风化壳储层底界面识别方法中花岗岩风化壳正演模拟图；
图2b为一种基于数理统计的风化壳储层底界面识别方法中花岗岩风化壳正演模拟图；
图3为一种基于数理统计的风化壳储层底界面识别方法中全区风化壳底界面标定结果图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
下面结合实施例对本发明作进一步的描述。
实施例1
请参照图1所示，一种基于数理统计的风化壳储层底界面识别方法，包括以下步骤：
步骤一、以井资料为依据构建岩石物理模型。
岩石物理模型为楔形模型，所述楔形模型以强风化带和中风化带作为可成为储集空间的风化壳储集层、以盖层为砂泥岩储集层和以致密的花岗岩为基岩建立。
步骤二、采用地震正演模拟方法研究确定不同厚度花岗岩风化壳储层与地震响应之间的联系。
地震正演模拟方法是指通过波动方程正演获得叠加剖面，模拟不同厚度的风化壳储集层在地震剖面上的响应。
如图2a和图2b所示，图中a.风化壳楔形模型图，b.风化壳波动方程正演模拟结果图，由图2a和图2b可知，楔形模型能够全方位的展示不同厚度的风化壳储集层在地震剖面上的响应，通过波动方程正演获得叠加剖面，可以看出当花岗岩风化壳厚度不足25米时，风化壳砾质风化带底部至黏土带顶，波阻抗呈一直增大趋势，形成一个低频单峰反射，风化壳顶底界面难以有效识别、区分；当风化壳厚度大于25米时，风化壳砾质风化带底面和黏土带顶面可反应为两条连续的强振幅反射，对应两个波峰，可对风化壳顶底面进行有效区分。
因此认为地震剖面上潜山顶面对应的强波峰代表潜山花岗岩风化壳的构造形态，而底界面的识别与风化壳的厚度有着较为直接的关系，古潜山的高度与受风化的时间即风化壳储集层厚度有着直接的关系，通过顶界面追踪与潜山古地貌形态的恢复能够实现风化壳底界面的定量识别和追踪。
步骤三、从古潜山暴露时间与风化程度关系分析入手，通过印模法研究确定花岗岩潜山古地貌，并在古地貌图上选取不同暴露时间点，统计不同暴露时间点在地震剖面上对应的风化壳储集层顶底界面相对时间差，以正演模拟为基准进行多样点拟合分析，获得相应的定量关系式。
印模法是指将前古近系花岗岩潜山顶面结束剥蚀开始上覆沉积时期视为一等时面，根据沉积补偿原理，利用上覆新沉积地层与残余古地貌之间存在镜像关系反映古地貌形态。
多点拟合分析的方法为在研究区中选取古地貌不同高度特征的 10个点，并选取其进行地震剖面反射特征分析。
定量关系式为：y＝-0.0136x2+1.2088x+11.143，式中，x表示古地貌基准点上的时间高度，y表示该点对应地震剖面上风化壳储集层顶底界面时间差。
步骤四、以全区风化壳储层顶界面为输入获得全区风化壳储层底界面标定结果，以及风化壳储层时间厚度定量结果图。
通过对古地貌风化零界面以上的古地貌时间高度进行提取，作为拟合公式的输入，获得全区风化壳时间厚度分布图。
在步骤四后，验证风化壳顶底标定结果的准确性。
如图3所示，其中图3a为已知风化壳顶界面全区平面图，图3c 为古地貌通过拾取拟合得到的全区风化壳储集层时间厚度图，可以看出该区风化壳储集层的分布特征，与已知钻井YL-A钻探结果吻合度高，满足地质规律，图3b为顶界面通过叠加厚度得到的风化壳底界面，可以看出顶底界面在保持良好的一致性的情况下，在潜山风化程度高的地方有良好的指示性，在风化程度低的地方也有较好的准确性，通过顶底界面的追踪可以直接圈定风化壳储集层的分布情况，为油气的圈定提供直接的依据。
本发明从已钻探花岗岩潜山风化壳储集层物性特征出发，进行正演模拟获得相应的地震响应特征，基于古潜山暴露时间与风化程度之间存在的一定关系，通过印模法恢复潜山古地貌，进行多样点统计拟合分析，获得古地貌不同高点与风化壳储层在地震剖面上顶底时差的定量关系式，运用该关系式进行全区风化壳底界面的识别与标定工作，本发明的方法可行性和可靠性极佳，为花岗岩潜山风化壳的圈定工作提供了一个良好的思路和方法。
以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不会使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (8)

1.一种基于数理统计的风化壳储层底界面识别方法，其特征在于，包括以下步骤：
S1、以井资料为依据构建岩石物理模型；
S2、采用地震正演模拟方法研究确定不同厚度花岗岩风化壳储层与地震响应之间的联系；
S3、从古潜山暴露时间与风化程度关系分析入手，通过印模法研究确定花岗岩潜山古地貌，并在古地貌图上选取不同暴露时间点，统计不同暴露时间点在地震剖面上对应的风化壳储集层顶底界面相对时间差，以正演模拟为基准进行多样点拟合分析，获得相应的定量关系式；
S4、以全区风化壳储层顶界面为输入获得全区风化壳储层底界面标定结果，以及风化壳储层时间厚度定量结果图。
2.根据权利要求1所述的一种基于数理统计的风化壳储层底界面识别方法，其特征在于，在所述步骤S1中，所述岩石物理模型为楔形模型，所述楔形模型以强风化带和中风化带作为可成为储集空间的风化壳储集层、以盖层为砂泥岩储集层和以致密的花岗岩为基岩建立。
3.根据权利要求1所述的一种基于数理统计的风化壳储层底界面识别方法，其特征在于，在所述步骤S2中，所述地震正演模拟方法是指通过波动方程正演获得叠加剖面，模拟不同厚度的风化壳储集层在地震剖面上的响应。
4.根据权利要求1所述的一种基于数理统计的风化壳储层底界面识别方法，其特征在于，在所述步骤S3中，所述印模法是指将前古近系花岗岩潜山顶面结束剥蚀开始上覆沉积时期视为一等时面，根据沉积补偿原理，利用上覆新沉积地层与残余古地貌之间存在镜像关系反映古地貌形态。
5.根据权利要求1所述的一种基于数理统计的风化壳储层底界面识别方法，其特征在于，在所述步骤S3中，所述多点拟合分析的方法为在研究区中选取古地貌不同高度特征的10个点，并选取其进行地震剖面反射特征分析。
6.根据权利要求1所述的一种基于数理统计的风化壳储层底界面识别方法，其特征在于，在所述步骤S3中，所述定量关系式为：y＝-0.0136x2+1.2088x+11.143，式中，x表示古地貌基准点上的时间高度，y表示该点对应地震剖面上风化壳储集层顶底界面时间差。
7.根据权利要求1所述的一种基于数理统计的风化壳储层底界面识别方法，其特征在于，在所述步骤S4中，通过对古地貌风化零界面以上的古地貌时间高度进行提取，作为拟合公式的输入，获得全区风化壳时间厚度分布图。
8.根据权利要求1所述的一种基于数理统计的风化壳储层底界面识别方法，其特征在于，在所述步骤S4后，验证风化壳顶底标定结果的准确性。
CN202110886389.4A 2021-08-03 2021-08-03 一种基于数理统计的风化壳储层底界面识别方法 Pending CN113433589A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110886389.4A CN113433589A (zh) 2021-08-03 2021-08-03 一种基于数理统计的风化壳储层底界面识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110886389.4A CN113433589A (zh) 2021-08-03 2021-08-03 一种基于数理统计的风化壳储层底界面识别方法

Publications (1)

Publication Number Publication Date
CN113433589A true CN113433589A (zh) 2021-09-24

Family

ID=77762688

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110886389.4A Pending CN113433589A (zh) 2021-08-03 2021-08-03 一种基于数理统计的风化壳储层底界面识别方法

Country Status (1)

Country Link
CN (1) CN113433589A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114265116A (zh) * 2022-03-02 2022-04-01 中海油研究总院有限责任公司 花岗岩类潜山风化壳沟脊幅度定量统计方法及装置

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103048690A (zh) * 2012-12-11 2013-04-17 成都理工大学 基于最优地震子波提取的快速匹配投影分解的地层反射拾取技术
RU2572462C1 (ru) * 2014-09-25 2016-01-10 ФГУП Центральный научно-исследовательский институт геологии нерудных полезных ископаемых Способ прогноза и поиска хемокластогенных магнезитов в кайнозойских депрессионных структурах
CN106873032A (zh) * 2017-02-16 2017-06-20 中国海洋石油总公司 一种花岗岩潜山储层裂缝预测方法
CN106886043A (zh) * 2017-03-01 2017-06-23 成都理工大学 基于地震数据深度学习的储层检测方法
CN108196302A (zh) * 2017-11-28 2018-06-22 中国石油天然气股份有限公司 白云岩缝洞发育储层甜点区预测方法及装置
EP3399143A1 (en) * 2015-12-30 2018-11-07 PetroChina Company Limited Method and device for determining dissolution degree of karst reservoir
CN109143399A (zh) * 2017-06-28 2019-01-04 中国石油化工股份有限公司 一种识别碳酸盐岩层序界面的方法
CN111474601A (zh) * 2019-01-24 2020-07-31 中国石油天然气股份有限公司 一种逐级预测火山岩优质储层的方法
CN111751877A (zh) * 2019-03-26 2020-10-09 中国石油天然气股份有限公司 一种地震数据多重积分相干裂缝预测方法与装置

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103048690A (zh) * 2012-12-11 2013-04-17 成都理工大学 基于最优地震子波提取的快速匹配投影分解的地层反射拾取技术
RU2572462C1 (ru) * 2014-09-25 2016-01-10 ФГУП Центральный научно-исследовательский институт геологии нерудных полезных ископаемых Способ прогноза и поиска хемокластогенных магнезитов в кайнозойских депрессионных структурах
EP3399143A1 (en) * 2015-12-30 2018-11-07 PetroChina Company Limited Method and device for determining dissolution degree of karst reservoir
CN106873032A (zh) * 2017-02-16 2017-06-20 中国海洋石油总公司 一种花岗岩潜山储层裂缝预测方法
CN106886043A (zh) * 2017-03-01 2017-06-23 成都理工大学 基于地震数据深度学习的储层检测方法
CN109143399A (zh) * 2017-06-28 2019-01-04 中国石油化工股份有限公司 一种识别碳酸盐岩层序界面的方法
CN108196302A (zh) * 2017-11-28 2018-06-22 中国石油天然气股份有限公司 白云岩缝洞发育储层甜点区预测方法及装置
CN111474601A (zh) * 2019-01-24 2020-07-31 中国石油天然气股份有限公司 一种逐级预测火山岩优质储层的方法
CN111751877A (zh) * 2019-03-26 2020-10-09 中国石油天然气股份有限公司 一种地震数据多重积分相干裂缝预测方法与装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
邱宇 等: "琼东南盆地深水区松南低凸起花岗岩风化壳储层底界面识别技术", 《海洋地质前沿》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114265116A (zh) * 2022-03-02 2022-04-01 中海油研究总院有限责任公司 花岗岩类潜山风化壳沟脊幅度定量统计方法及装置

Similar Documents

Publication Publication Date Title
CN111766630B (zh) 一种盆地古地貌的恢复方法
US10365387B1 (en) Method for secondary exploration of old oil area in fault subsidence basin
US5586082A (en) Method for identifying subsurface fluid migration and drainage pathways in and among oil and gas reservoirs using 3-D and 4-D seismic imaging
O'Byrne et al. Sequence, parasequence, and intraparasequence architecture of the Grassy Member, Blackhawk Formation, Book Cliffs, Utah, USA
CN105651676A (zh) 一种水平井规则开发井网下的储层非均质性表征方法
CN109541685B (zh) 一种河道砂体识别方法
Pringle et al. Capturing stratigraphic and sedimentological complexity from submarine channel complex outcrops to digital 3D models, Karoo Basin, South Africa
Gebreyohannes et al. Regional groundwater flow modeling of the Geba basin, northern Ethiopia
CN111706317B (zh) 一种确定加密调整区低渗储层剩余油分布状况的方法
Miller et al. Back to basics of sequence stratigraphy: early Miocene and mid-Cretaceous examples from the New Jersey paleoshelf
CN112394404B (zh) 一种渐进式储层精细表征方法
CN113093286B (zh) 一种冷泉发育区储层非均质性的反演方法
CN111768503B (zh) 一种基于三维地质模型的海砂资源量估算方法
CN113433589A (zh) 一种基于数理统计的风化壳储层底界面识别方法
CN106443772B (zh) 一种去底辟原始地层厚度恢复方法
CN109584365B (zh) 一种用于海洋吹填造陆工程土方量造价核算的Voronoi图地质分层模型构建方法
Peijs-van Hilten et al. Heterogeneity modeling and geopseudo upscaling applied to waterflood performance prediction of an incised valley reservoir: Countess YY Pool, southern Alberta, Canada
Erzeybek Balan Characterization and modeling of paleokarst reservoirs using multiple-point statistics on a non-gridded basis
Kantak et al. Sediment thickness and a west-east geologic cross section in the Caracas Valley
Sweetkind et al. Construction of 3-D geologic framework and textural models for Cuyama Valley groundwater basin, California
Curkan Reservoir characterization of channel-belt strata, McMurray Formation, northeastern Alberta
Zhang et al. A novel three-dimensional mine area hydrogeological model based on groundwater modeling systems
Kreitler Lineations and faults in the Texas coastal zone
CN114076982B (zh) 一种基于波形特征差异的岩溶古地貌恢复方法及装置
Paulen et al. Bedrock topography and drift thickness models from the Timmins area, northeastern Ontario: An application of GIS to the Timmins overburden drillhole database

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20210924

RJ01 Rejection of invention patent application after publication