CN107589447A - 基于曲率的地层压力校正方法 - Google Patents
基于曲率的地层压力校正方法 Download PDFInfo
- Publication number
- CN107589447A CN107589447A CN201610533977.9A CN201610533977A CN107589447A CN 107589447 A CN107589447 A CN 107589447A CN 201610533977 A CN201610533977 A CN 201610533977A CN 107589447 A CN107589447 A CN 107589447A
- Authority
- CN
- China
- Prior art keywords
- strata pressure
- curvature
- data volume
- pressure
- compressional structure
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
一种基于曲率的地层压力校正方法,包括以下步骤:步骤1:获得地层压力三维数据体P;步骤2:针对与所述三维数据体P相同范围的地震解释层位,提取曲率属性,将所述曲率属性归一化,作为挤压构造对地层压力的贡献率a;步骤3:基于地层资料拟合挤压构造应力函数σT;步骤4:根据以下公式计算校正后的地层压力三维数据体PM:PM=P+a·σT。该方法在地层压力校正过程中引入曲率属性,在空间上可以有效地反映挤压构造对地层压力的贡献,能够获得更准确的地层压力预测结果。
Description
技术领域
本发明涉及地球物理技术,特别涉及一种基于曲率的地层压力校正方法。
背景技术
获得准确的地层压力分布是水平井压裂技术的关键,对于页岩气的开发也是至关重要。目前用于地层压力预测、监测和检测的方法主要有地质分析法、地震法、dc指数法和测井法等。地质分析法与地震法都属于钻前预测,在普查新区的深部地层、无钻井和测井资料时,可以用来宏观划定异常压力带,但精度较低。dc指数法属于随钻预测,对岩性的判断不够准确,计算压力值也受到一定的影响。此外,目前基于三维地震数据的地层压力预测技术中并未考虑断层因素的影响,因此在断层附近的地层压力结果与实际地质情况存在着一定的误差。业内期待一种更准确的地层压力预测方法,以获得更符合实际地质情况的三维地层压力数据。
发明内容
本发明的目的是提供一种基于曲率的地层压力校正方法,其通过分析断层附近地层压力变化等因素,能够获得更准确的地层压力预测结果。
本发明采用以下解决方案:
一种基于曲率的地层压力校正方法,包括以下步骤:
步骤1:获得地层压力三维数据体P;
步骤2:针对与所述三维数据体P相同范围的地震解释层位,提取曲率属性,将所述曲率属性归一化,作为挤压构造对地层压力的贡献率a;
步骤3:基于地层资料拟合挤压构造应力函数σT;
步骤4:根据以下公式(4)计算校正后的地层压力三维数据体PM:
PM=P+α·σT (4)。
优选地,所述步骤3是通过以下步骤实现的:
收集地层资料,获得挤压构造应力函数σT的离散值;
基于所述离散值,根据以下公式(3)求取拟合参数c和d:
σT=(H-c)ρ0g/d (3)
其中,ρ0为地层的平均密度;g为重力加速度;H为泥岩声波时差平直段的开始深度;
通过所述公式(3)拟合挤压构造应力函数σT。
与现有技术相比,本发明的有益效果在于:
(1)本发明在地层压力校正过程中引入曲率属性,在空间上可以有效地反映挤压构造对地层压力的贡献,并且提取曲率属性是成熟的技术,在实际应用过程中易于实现;
(2)经过本发明修正的地层压力预测结果更符合地质规律,尤其是在地质构造变化剧烈的地区,更符合实际地质情况。
附图说明
通过结合附图对本公开示例性实施例进行更详细的描述,本公开的上述以及其它目的、特征和优势将变得更加明显。
图1示出了根据示例性实施例的基于曲率的地层压力校正方法的流程图。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施例。虽然附图中显示了本公开的优选实施例,然而应该理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本公开更加透彻和完整,并且能够将本公开的范围完整地传达给本领域的技术人员。
本发明的示例性实施例基于常规的三维地层压力预测数据,以曲率属性和测井数据为基础,在断层附近进行地层压力校正,最终获得更符合实际地质情况的三维地层压力数据,预测结果可以用于页岩气等热点领域。
以下参考图1详细描述本发明的基于曲率的地层压力校正方法的原理和步骤。
杨进等(2009)针对挤压构造提出地层孔隙压力预测的修正模型,如以下公式(1)所示:
PM=K·P+a·σT (1)
其中,PM为校正后的地层压力预测值,单位为MPa;P为常规模型计算的地层压力值,单位为MPa;a为挤压构造对地层压力的贡献率,单位为%;K为构造压力因子;σT为相对挤压构造应力。
测井资料受人为因素影响较少,且能随井深连续变化,资料连续性好,纵向分辨率高、可靠性好,地层的声波速度、密度和电阻率等参数都与地层孔隙压力存在一定的关系,能够表现出一定的规律性。因此,目前通常利用测井资料预测地层压力。公式(1)中提到的常规模型利用测井资料中对地层孔隙压力敏感的曲线对地层压力进行拟合,是一个非线性回归的过程。在常规模型中,只考虑测井资料与地层孔隙压力的关系,未考虑其它地质因素的影响。应用常规模型预测地层压力是现有技术,本领域技术人员可根据实际需要选择适当的模型进行地层压力预测。
挤压构造应力对地层压力的作用可视为侧向的压力作用,其主要表现在两个方面:压力作用使孔隙体积降低,岩石的渗透能力因孔隙度的减小而降低。因此可以将构造压力因子K定义为孔隙度的函数,如以下公式(2)所示:
其中,K为构造压力因子,为随深度H变化的地层孔隙度,b为待定系数。
根据地层资料拟合本地区的挤压构造应力函数,如以下公式(3)所示:
σT=(H-c)ρ0g/d (3)
其中,PM为相对构造挤压应力,单位为MPa;ρ0为地层的平均密度,单位为g/cm3;g为重力加速度,单位为m/s2;H为泥岩声波时差平直段的开始深度,单位为m;c和d为拟合参数。
如果忽略构造压力因子K在地层孔隙压力修正模型的作用,即K=1,则以上公式(1)可简化为公式(4):
PM=P+a·σT (4)
其中,PM为校正后的地层压力预测值,单位为MPa;P为常规模型计算的地层压力值,单位为MPa;σT为相对挤压构造应力;a为挤压构造对地层压力的贡献率,单位为%。
依据以下两点原因,可以在以上公式(4)中忽略参数K:(1)构造压力因子K主要体现为侧向的压力作用,而公式(4)中通过曲率属性获得的参数a体现的是挤压构造对地层压力的影响,因此在一定程度上抵消了忽略构造压力因子K对校正结果的影响;(2)计算构造压力因子K需要获得三维的孔隙度数据,原公式(1)主要用于二维地层孔隙压力校正,而根据示例性实施例的方法主要针对三维数据进行地层孔隙压力校正,就需要考虑三维的孔隙度数据的可靠性,众所周知三维孔隙度数据的获取难度非常大,因此可忽略参数K。
依据以下两点原因,在示例性实施例中可以采用曲率属性近似获得挤压构造对地层压力的贡献率a:(1)曲率属性用于描述地质体的几何变化,与地震反射体的弯曲程度相对应,对岩层的弯曲、褶皱和裂缝、断层等反应敏感,是寻找地质体构造特征的有效手段,此外其能很好地反映由沉积作用或低级序断层引起的波形变化;(2)曲率属性基于构造解释层位提取,容易获取且精度较高,曲率属性绝对值较大时地层弯曲程度大,地层受到的挤压应力也较大。
基于以上描述的基本原理,根据示例性实施例的基于曲率的地层压力校正方法包括以下步骤:
步骤1:获得地层压力三维数据体P
地层压力三维数据体P是根据现有的常规模型计算的。如前所述,常规模型利用测井资料中对地层孔隙压力敏感的曲线对地层压力进行拟合,是一个非线性回归的过程。在常规模型中,只考虑测井资料与地层孔隙压力的关系,未考虑其它地质因素的影响。应用常规模型预测地层压力是现有技术,本领域技术人员可根据实际需要选择适当的模型进行地层压力预测,获得地层压力三维数据体P。
步骤2:针对与所述三维数据体P相同范围的地震解释层位,提取曲率属性,将所述曲率属性归一化,作为挤压构造对地层压力的贡献率a
曲率属性反应的是地震反射体的弯曲程度,即曲率属性越大地震反射体弯曲程度越大,所产生的挤压构造的应力越大,即挤压构造对地层压力的贡献率a越大。因此,在本步骤中,将提取的曲率属性进行归一化,作为挤压构造对地层压力的贡献率a。针对某一地震解释层位提取曲率属性是非常成熟的技术,在实际应用过程中易于通过曲率属性获得挤压构造对地层压力的贡献率a。
步骤3:基于地层资料拟合挤压构造应力函数σT
收集地层资料可以获得挤压构造应力函数σT的离散值,基于这些离散值,可以通过以下公式(3)求取拟合参数c和d:
σT=(H-c)ρ0g/d (3)
其中,σT为挤压构造应力函数,单位为MPa;ρ0为地层的平均密度,单位为g/cm3;g为重力加速度,单位为m/s2;H为泥岩声波时差平直段的开始深度,单位为m;c和d为拟合参数。
在求取拟合参数c和d之后,即可根据公式(3)拟合挤压构造应力函数σT,以便计算任意一点的挤压构造应力。
步骤4:根据以下公式(4)计算校正后的地层压力三维数据体PM:
PM=P+a·σT (4)。
应用示例
在本示例中,按照以下步骤进行基于曲率的地层压力校正:
步骤1:根据已知常规模型计算地层压力三维数据体P,并根据工区实际情况,获得工区范围和目的层位信息。
步骤2:针对与三维数据体P相同范围的地震解释层位,提取曲率属性,对曲率属性进行归一化,将其取值范围限制在[0,1]区间内,将归一化后的曲率属性作为挤压构造对地层压力的贡献率a。由于贡献率a为沿层属性,与三维数据体P在纵向上(即z方向)范围不同,只有x、y方向范围相同,为了获得与三维数据体P在纵向上范围一致的贡献率a,将a的值在纵向上进行复制,即贡献率a的纵向上为常值,只有x、y方向的值是变化的。
步骤3:通过工区地层资料获得靶区的挤压构造应力函数σT的离散值,基于该离散值,根据公式(3)求取拟合参数c=500,d=1000。因此,拟合的挤压构造应力函数σT为:
σT=(H-500)ρ0g/1000。
步骤4:根据公式(4)计算校正后的地层压力三维数据体PM。
获得的地层压力三维数据体PM可以用于页岩气开发评价等方面,在页岩气有利预测预测中起着至关重要的作用。
上述技术方案只是本发明的一种实施例,对于本领域内的技术人员而言,在本发明公开的原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施例的描述,因此前面的描述只是优选的,而并不具有限制性的意义。
Claims (2)
1.一种基于曲率的地层压力校正方法,包括以下步骤:
步骤1:获得地层压力三维数据体P;
步骤2:针对与所述三维数据体P相同范围的地震解释层位,提取曲率属性,将所述曲率属性归一化,作为挤压构造对地层压力的贡献率a;
步骤3:基于地层资料拟合挤压构造应力函数σT;
步骤4:根据以下公式(4)计算校正后的地层压力三维数据体PM:
PM=P+a·σT (4)。
2.根据权利要求1所述的基于曲率的地层压力校正方法,其中所述步骤3是通过以下步骤实现的:
收集地层资料,获得挤压构造应力函数σT的离散值;
基于所述离散值,根据以下公式(3)求取拟合参数c和d:
σT=(H-c)ρ0g/d (3)
其中,ρ0为地层的平均密度;g为重力加速度;H为泥岩声波时差平直段的开始深度;
通过所述公式(3)拟合挤压构造应力函数σT。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610533977.9A CN107589447A (zh) | 2016-07-07 | 2016-07-07 | 基于曲率的地层压力校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610533977.9A CN107589447A (zh) | 2016-07-07 | 2016-07-07 | 基于曲率的地层压力校正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107589447A true CN107589447A (zh) | 2018-01-16 |
Family
ID=61045944
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610533977.9A Pending CN107589447A (zh) | 2016-07-07 | 2016-07-07 | 基于曲率的地层压力校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107589447A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113187463A (zh) * | 2021-04-14 | 2021-07-30 | 中海石油(中国)有限公司海南分公司 | 一种基于地层超压单成因贡献率的随钻孔隙压力预测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011005353A1 (en) * | 2009-07-06 | 2011-01-13 | Exxonmobil Upstream Research Company | Method for seismic interpretation using seismic texture attributes |
CN102200588A (zh) * | 2011-03-22 | 2011-09-28 | 成都理工大学 | 地震资料的波形相似性体曲率分析方法 |
CN103982179A (zh) * | 2014-05-26 | 2014-08-13 | 中国地质大学(北京) | 一种油藏储层的古压力定量反演探测方法 |
CN104267429A (zh) * | 2014-09-30 | 2015-01-07 | 中国石油天然气股份有限公司 | 确定地层压力的方法及装置 |
CN105445791A (zh) * | 2015-11-25 | 2016-03-30 | 成都理工大学 | 一种基于多种地震属性的地层孔隙压力预测方法 |
CN105468886A (zh) * | 2014-09-04 | 2016-04-06 | 中国石油化工股份有限公司 | 基于岩石物理学参数计算地层压力的方法 |
-
2016
- 2016-07-07 CN CN201610533977.9A patent/CN107589447A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011005353A1 (en) * | 2009-07-06 | 2011-01-13 | Exxonmobil Upstream Research Company | Method for seismic interpretation using seismic texture attributes |
CN102200588A (zh) * | 2011-03-22 | 2011-09-28 | 成都理工大学 | 地震资料的波形相似性体曲率分析方法 |
CN103982179A (zh) * | 2014-05-26 | 2014-08-13 | 中国地质大学(北京) | 一种油藏储层的古压力定量反演探测方法 |
CN105468886A (zh) * | 2014-09-04 | 2016-04-06 | 中国石油化工股份有限公司 | 基于岩石物理学参数计算地层压力的方法 |
CN104267429A (zh) * | 2014-09-30 | 2015-01-07 | 中国石油天然气股份有限公司 | 确定地层压力的方法及装置 |
CN105445791A (zh) * | 2015-11-25 | 2016-03-30 | 成都理工大学 | 一种基于多种地震属性的地层孔隙压力预测方法 |
Non-Patent Citations (1)
Title |
---|
杨进 等: "挤压构造地层压力预测模型研究", 《石油学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113187463A (zh) * | 2021-04-14 | 2021-07-30 | 中海石油(中国)有限公司海南分公司 | 一种基于地层超压单成因贡献率的随钻孔隙压力预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106368691B (zh) | 基于岩石物理地震信息三维异常孔隙压力预测方法 | |
Zhang et al. | Integration of discrete fracture reconstruction and dual porosity/dual permeability models for gas production analysis in a deformable fractured shale reservoir | |
CN105277982B (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
KR101415199B1 (ko) | 전기비저항 및 밀도 검층 데이터를 이용한 셰일가스전 스윗스팟 수평시추구간의 슬로우니스값, 영률, 포아송비 및 취성 추정방법 | |
CN104453873A (zh) | 页岩油气经济有效层段的评价方法 | |
CN107907910B (zh) | 一种不同岩性油藏横波测井确定方法 | |
CN105545299B (zh) | 一种地下不同层位煤岩层中氡气析出浓度原位测定方法 | |
US11879316B2 (en) | Geostatistical analysis of microseismic data in fracture modeling | |
CN106295042B (zh) | 一种煤层顶板稳定性测井定量评价方法 | |
CN104502997A (zh) | 一种利用裂缝密度曲线预测裂缝密度体的方法 | |
KR101415196B1 (ko) | 전기비저항 및 밀도 검층 데이터를 이용한 셰일가스전 스윗스팟 수평시추구간의 파쇄구간 선정방법 | |
CN105697002A (zh) | 一种用于识别煤系地层岩性的方法 | |
CN106599377A (zh) | 基于测井资料和煤体结构参量定量划分煤体结构的方法 | |
CN105467438A (zh) | 一种基于三模量的泥页岩地应力三维地震表征方法 | |
KR101415197B1 (ko) | 전기비저항 및 중성자 검층 데이터를 이용한 셰일가스전 스윗스팟 수평시추구간의 파쇄구간 선정방법 | |
CN105931125A (zh) | 一种致密油分段多簇体积压裂水平井产量预测方法 | |
Chen et al. | Prediction of sweet spots in shale reservoir based on geophysical well logging and 3D seismic data: A case study of Lower Silurian Longmaxi Formation in W4 block, Sichuan Basin, China | |
CN105370270B (zh) | 由偶极声波纵、横波时差确定页岩气储层含气饱和度的方法 | |
CN104483706B (zh) | 一种基于煤岩力学参数的煤体结构测井定量识别方法 | |
CN104153768A (zh) | 一种评价花岗岩储层储集性能的方法 | |
CN107728205A (zh) | 一种地层压力预测方法 | |
WO2018231594A1 (en) | Estimation of mechanical properties of transversely isotropic media | |
US10400591B2 (en) | Systems and methods for acoustic testing of laminated rock to determine total organic carbon content | |
CN108661632B (zh) | 一种高精度碳氧比测井方法 | |
Wang et al. | Effects of in-situ stress and joint on permeability of the coal bed in Linfen block, southeastern Ordos Basin, China |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180116 |