CN109931055B - 盆地深层复合成因的流体压力预测方法 - Google Patents
盆地深层复合成因的流体压力预测方法 Download PDFInfo
- Publication number
- CN109931055B CN109931055B CN201910095928.5A CN201910095928A CN109931055B CN 109931055 B CN109931055 B CN 109931055B CN 201910095928 A CN201910095928 A CN 201910095928A CN 109931055 B CN109931055 B CN 109931055B
- Authority
- CN
- China
- Prior art keywords
- compaction
- effective stress
- vertical effective
- fluid pressure
- region
- 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.)
- Expired - Fee Related
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000002131 composite material Substances 0.000 title claims abstract description 17
- 238000005056 compaction Methods 0.000 claims abstract description 39
- 238000009933 burial Methods 0.000 claims abstract description 23
- 230000008021 deposition Effects 0.000 claims abstract description 23
- 238000013178 mathematical model Methods 0.000 claims abstract description 15
- 238000010276 construction Methods 0.000 claims abstract description 10
- 230000009471 action Effects 0.000 claims abstract description 8
- 210000003194 forelimb Anatomy 0.000 claims description 6
- 238000002679 ablation Methods 0.000 claims 2
- 230000015572 biosynthetic process Effects 0.000 abstract description 2
- 238000001125 extrusion Methods 0.000 description 5
- 239000011435 rock Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000002706 hydrostatic effect Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供盆地深层复合成因的流体压力预测方法,调查研究区的构造、沉积背景,遴选出研究区内构造稳定、沉积同物源的地区;对构造稳定和沉积同物源地区进行泥岩压实规律研究,得到最大埋深期泥岩压实曲线斜率的数学模型;求取构造复杂地区最大埋深期的压实与加载曲线,通过声速与垂向有效应力作图得到加载曲线;求取构造作用前后垂向有效应力差值的数学模型,最大埋深期地层的垂向有效应力在受到挤压、抬升或断裂作用时,垂向有效应力将减小到现今大小,对于该有效应力差值的影响因素,构建地质模型,在此基础上,研究其与距断层或背斜前肢距离及剥蚀厚度的关系,得到构造作用前后垂向有效应力差值的数学模型;求取复合成因的流体压力。
Description
技术领域
本发明属于油气勘探流体压力预测技术领域,具体涉及盆地深层复合成因的流体压力预测方法。
技术背景
精确的流体压力钻前预测在油气勘探中至关重要,历来受到油气地质学家和油气工程师的高度重视。
现有的盆地流体压力预测方法较多,如等效深度法,压力=深度下的岩石平均密度*g*深度-(等效深度下的岩石平均密度-等效深度下的流体平均密度)*g*深度,正常压实趋势法,压力=静水压力+﹛[(岩石平均密度-地层水平均密度)*g*深度差]/[log(h2测井值/h1测井值)+0.435*测井曲线上任一点的温度系数*两高度间的低温梯度*高度差]﹜*log(估算具有正常压实趋势地区的特定地球物理参数/估算具有异常孔隙压力地区的特定地球物理参数),压缩曲线法(Dobrynin et al.,1989)、Eaton幂指数法(Eaton,1972)、Dutta法(Dutta,1998,2002)、Fillippone法(Fillippone,1982)等,但这些方法大多只能计算由孔隙或温度异常等单一因素引起的流体压力,当存在复合成因的流体压力(如欠压实、构造挤压及超压传递等)时,使用这些“单打一”的方法将无法进行准确的压力预测。
发明内容
本发明提供盆地深层复合成因的流体压力预测方法,解决复杂地质背景下流体压力预测的难题。
本发明的技术方案是:盆地深层复合成因的流体压力预测方法,包括
(1)调查清楚研究区的构造、沉积背景,遴选出研究区内构造稳定、沉积同物源的地区;
(2)对构造稳定和沉积同物源区进行泥岩压实规律研究:
首先,利用地球物理测井资料,读取纯泥岩段的测井序列值;其次,分别绘制所述测井序列值与埋深的关系图,制作综合压实曲线,得到压实斜率;最后,研究影响压实斜率的地质因素本身及其与压实斜率的经验关系,所述地质因素包括沉积速率、地层密度,最终得到压实斜率与沉积速率的关系,亦即最大埋深期泥岩压实曲线斜率的数学模型,c=F(V),其中,V为沉积速率,ρ为地层密度,c为泥岩压实曲线斜率;
(3)求取构造复杂地区最大埋深期的压实与加载曲线,将步骤(2)中构造稳定、沉积同物源地区求取的最大埋深期泥岩压实规律数学模型,应用于构造复杂地区,求取该地区的最大埋深期的压实曲线,进而通过声速与垂向有效应力作图得到加载曲线VES=G(VP),其中VP为声速,VES为垂向有效应力;
(4)求取构造作用前后垂向有效应力差值的数学模型,最大埋深期地层的垂向有效应力在受到挤压、抬升或断裂作用时,垂向有效应力将减小到现今大小,对于该有效应力差值的影响因素,首先构建地质模型,在此基础上,研究其与距断层或背斜前肢距离及剥蚀厚度的关系,最终得到构造作用前后垂向有效应力差值的数学模型△VES=D(H/FDIST),其中△VES是有效应力差值,FDIST为距断层或背斜前肢距离,H为剥蚀厚度;
(5)复合成因的流体压力的求取,P=S-(VES-△VES),VES=G(VP),△VES=D(H/FDIST),S为上覆负荷,VP为声速,p为流体压力。
方案进一步地,上述步骤(1)中所述构造稳定、沉积同物源的地区是无明显水平挤压、断裂作用的地区。
方案进一步地,上述步骤(2)中所述测井序列值包括声波、密度、电阻率、自然伽马。
方案进一步地,上述步骤(3)中所述构造复杂地区为存在构造挤压、抬升剥蚀和断裂作用的地区。
本发明的优点是:本发明提出了复杂地质背景下,盆地深层复合成因的流体压力预测方法,具有考虑因素多、预测准、实用性强等优点,解决了深层复合成因流体压力预测的难题。
附图说明
图1是本发明流程示意图;
具体实施方式
下面结合附图对本发明做清楚完整的描述,以使本领域的技术人员在不需要作出创造性劳动的条件下,能够充分实施本发明。
本发明的具体实施方式是:如图1所示,盆地深层复合成因的流体压力预测方法,包括
(1)调查清楚研究区的构造、沉积背景,遴选出研究区内构造稳定、沉积同物源的地区,所述构造稳定、沉积同物源的地区是无明显水平挤压、断裂作用的地区;
(2)对构造稳定和沉积同物源区进行泥岩压实规律研究:
首先,利用地球物理测井资料,读取纯泥岩段的测井序列值,所述测井序列值包括声波、密度、电阻率、自然伽马;其次,分别绘制所述测井序列值与埋深的关系图,制作综合压实曲线,得到压实斜率;最后,研究影响压实斜率的地质因素本身及其与压实斜率的经验关系,所述地质因素包括沉积速率、地层密度,最终得到压实斜率与沉积速率的关系,亦即最大埋深期泥岩压实曲线斜率的数学模型,c=F(V),其中,V为沉积速率,ρ为地层密度,c为泥岩压实曲线斜率;
(3)求取构造复杂地区最大埋深期的压实与加载曲线,所述构造复杂地区为存在构造挤压、抬升剥蚀和断裂作用的地区,将步骤(2)中构造稳定、沉积同物源地区求取的最大埋深期泥岩压实规律数学模型,应用于构造复杂地区,求取该地区的最大埋深期的压实曲线,进而通过声速与垂向有效应力作图得到加载曲线VES=G(VP),其中VP为声速,VES为垂向有效应力;
(4)求取构造作用前后垂向有效应力差值的数学模型,最大埋深期地层的垂向有效应力在受到挤压、抬升或断裂作用时,垂向有效应力将减小到现今大小,对于该有效应力差值的影响因素,首先构建地质模型,在此基础上,研究其与距断层或背斜前肢距离及剥蚀厚度的关系,最终得到构造作用前后垂向有效应力差值的数学模型△VES=D(H/FDIST),其中△VES是有效应力差值,FDIST为距断层或背斜前肢距离,H为剥蚀厚度;
(5)复合成因的流体压力的求取,P=S-(VES-△VES),VES=G(VP),△VES=D(H/FDIST),S为上覆负荷,VP为声速,p为流体压力。
以上对本发明的较佳实施例进行了描述,需要指出的是,本发明并不局限于上述特定实施方式,其中未尽详细描述的设备和结构应该理解为用本领域中的普通方式予以实施;任何熟悉本领域的技术人员,在不脱离本发明技术方案范围情况下,依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化及修饰,均仍属于本发明技术方案保护的范围内。
Claims (4)
1.盆地深层复合成因的流体压力预测方法,其特征在于,包括
(1)调查清楚研究区的构造、沉积背景,遴选出研究区内构造稳定、沉积同物源的地区;
(2)对构造稳定和沉积同物源区进行泥岩压实规律研究:
首先,利用地球物理测井资料,读取纯泥岩段的测井序列值;其次,分别绘制所述测井序列值与埋深的关系图,制作综合压实曲线,得到压实斜率;最后,研究影响压实斜率的地质因素本身及其与压实斜率的经验关系,所述地质因素包括沉积速率、地层密度,最终得到压实斜率与沉积速率的关系,亦即最大埋深期泥岩压实曲线斜率的数学模型,c=F(V),其中,V为沉积速率,ρ为地层密度,c为泥岩压实曲线斜率;
(3)求取构造复杂地区最大埋深期的压实与加载曲线,将步骤(2)中构造稳定、沉积同物源地区求取的最大埋深期泥岩压实规律数学模型,应用于构造复杂地区,求取该地区的最大埋深期的压实曲线,进而通过声速与垂向有效应力作图得到加载曲线VES=G(VP),其中VP为声速,VES为最大埋深期垂向有效应力;
(4)求取构造作用前后垂向有效应力差值的数学模型,最大埋深期地层的垂向有效应力在受到挤压、抬升或断裂作用时,垂向有效应力将减小到现今大小,对于该垂向有效应力差值的影响因素,首先构建地质模型,在此基础上,研究其与距断层或背斜前肢距离及剥蚀厚度的关系,最终得到构造作用前后垂向有效应力差值的数学模型△VES=D(H/FDIST),其中△VES是最大埋深期垂向有效应力与现今垂向有效应力差值,FDIST为距断层或背斜前肢距离,H为剥蚀厚度;
(5)复合成因的流体压力的求取,P=S-(VES-△VES),VES=G(VP),△VES=D(H/FDIST),S为上覆负荷,VP为声速,p为流体压力。
2.根据权利要求1所述的盆地深层复合成因的流体压力预测方法,其特征在于,上述步骤(1)中所述构造稳定、沉积同物源的地区是无明显水平挤压、断裂作用的地区。
3.根据权利要求1所述的盆地深层复合成因的流体压力预测方法,其特征在于,上述步骤(2)中所述测井序列值包括声波、密度、电阻率、自然伽马。
4.根据权利要求1所述的盆地深层复合成因的流体压力预测方法,其特征在于,上述步骤(3)中所述构造复杂地区为存在构造挤压、抬升剥蚀和断裂作用的地区。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910095928.5A CN109931055B (zh) | 2019-01-31 | 2019-01-31 | 盆地深层复合成因的流体压力预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910095928.5A CN109931055B (zh) | 2019-01-31 | 2019-01-31 | 盆地深层复合成因的流体压力预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109931055A CN109931055A (zh) | 2019-06-25 |
CN109931055B true CN109931055B (zh) | 2020-09-25 |
Family
ID=66985340
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910095928.5A Expired - Fee Related CN109931055B (zh) | 2019-01-31 | 2019-01-31 | 盆地深层复合成因的流体压力预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109931055B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112664188B (zh) * | 2020-12-31 | 2023-09-22 | 中国海洋石油集团有限公司 | 一种基于砂体传递模式的地层孔隙压力预测方法 |
CN113187463B (zh) * | 2021-04-14 | 2022-08-02 | 中海石油(中国)有限公司海南分公司 | 一种基于地层超压单成因贡献率的随钻孔隙压力预测方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4981037A (en) * | 1986-05-28 | 1991-01-01 | Baroid Technology, Inc. | Method for determining pore pressure and horizontal effective stress from overburden and effective vertical stresses |
US5200929A (en) * | 1992-03-31 | 1993-04-06 | Exxon Production Research Company | Method for estimating pore fluid pressure |
US5282384A (en) * | 1992-10-05 | 1994-02-01 | Baroid Technology, Inc. | Method for calculating sedimentary rock pore pressure |
US6473696B1 (en) * | 2001-03-13 | 2002-10-29 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
US7299132B2 (en) * | 2005-08-08 | 2007-11-20 | Schlumberger Technology Corp. | Method and system for pre-drill pore pressure prediction |
CN104453879B (zh) * | 2014-11-14 | 2017-04-05 | 中国海洋石油总公司 | 钻前压力的预测方法 |
CN107015290B (zh) * | 2017-03-13 | 2018-05-08 | 西北大学 | 一种改造型断陷盆地原始沉积面貌恢复的方法 |
CN106979006A (zh) * | 2017-05-17 | 2017-07-25 | 中国神华能源股份有限公司 | 地层压力的确定方法和装置 |
CN108150160B (zh) * | 2017-12-21 | 2021-07-06 | 西安石油大学 | 一种地层中欠压实作用超压的求取方法 |
-
2019
- 2019-01-31 CN CN201910095928.5A patent/CN109931055B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN109931055A (zh) | 2019-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1966934B (zh) | 一种随钻预测钻头底下地层坍塌压力和破裂压力的方法 | |
CN105651676B (zh) | 一种水平井规则开发井网下的储层非均质性表征方法 | |
CN101025084B (zh) | 一种随钻预测钻头底下地层孔隙压力的方法 | |
CN108089227B (zh) | 新的基于三维地震资料的地层孔隙压力预测方法 | |
CN109931055B (zh) | 盆地深层复合成因的流体压力预测方法 | |
CN111706317B (zh) | 一种确定加密调整区低渗储层剩余油分布状况的方法 | |
CN112502701A (zh) | 一种低渗透储层综合地质-工程的分类评价方法 | |
CN113790044B (zh) | 一种多机制复合构成超压地层孔隙压力的评价方法 | |
KR20160056359A (ko) | Gis 기반 실시간 지진피해 예측 방법 | |
CN104122602A (zh) | 一种获取沉积盆地剥蚀厚度的方法及装置 | |
CN110187391B (zh) | 一种利用声波时差确定地层层序界面的方法 | |
CN114660269B (zh) | 一种古潜山天然气动态成藏过程恢复方法 | |
Ugwu et al. | Integrated approach to geopressure detection in the X-field, Onshore Niger Delta | |
Christeson et al. | Physical properties and seismic structure of Izu‐B onin‐M ariana fore‐arc crust: Results from IODP E xpedition 352 and comparison with oceanic crust | |
CN111025391B (zh) | 一种断层活动性的定量评价方法 | |
CN107478768B (zh) | 一种泥页岩古抗剪强度确定方法 | |
Dutta et al. | Estimation of formation fluid pressure using high-resolution velocity from inversion of seismic data and a rock physics model based on compaction and burial diagenesis of shales | |
Sun et al. | Characteristics and prediction of weathered volcanic rock reservoirs: A case study of Carboniferous rocks in Zhongguai paleouplift of Junggar Basin, China | |
CN115680638B (zh) | 一种利用压力衰减梯度识别超压顶封层的方法 | |
CN112241025B (zh) | 一种井震联合地层压力确定方法及系统 | |
CN110297280B (zh) | 一种碳酸盐岩超压空间分布特征的预测方法 | |
Woillez et al. | Insights in fault flow behaviour from onshore Nigeria petroleum system modelling | |
Yang et al. | Adit deformation prediction based on ICEEMDAN dispersion entropy and LSTM-BP | |
Chen et al. | Resonance oscillations of the Soufrière Hills Volcano (Montserrat, WI) magmatic system induced by forced magma flow from the reservoir into the upper plumbing dike | |
Wilson | Seismic evaluation of differential tectonic subsidence, compaction, and loading in an interior basin |
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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Wang Gang Inventor after: Fan Changyu Inventor after: Chen Jie Inventor after: Han Xiaojie Inventor before: Fan Changyu Inventor before: Wang Gang Inventor before: Chen Jie Inventor before: Han Xiaojie |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200925 |