CN109931055B - 盆地深层复合成因的流体压力预测方法 - Google Patents

盆地深层复合成因的流体压力预测方法 Download PDF

Info

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
ves
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
Application number
CN201910095928.5A
Other languages
English (en)
Other versions
CN109931055A (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.)
Northwestern University
Original Assignee
Northwestern University
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 Northwestern University filed Critical Northwestern University
Priority to CN201910095928.5A priority Critical patent/CN109931055B/zh
Publication of CN109931055A publication Critical patent/CN109931055A/zh
Application granted granted Critical
Publication of CN109931055B publication Critical patent/CN109931055B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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)中所述构造复杂地区为存在构造挤压、抬升剥蚀和断裂作用的地区。
CN201910095928.5A 2019-01-31 2019-01-31 盆地深层复合成因的流体压力预测方法 Expired - Fee Related CN109931055B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 西安石油大学 一种地层中欠压实作用超压的求取方法

Also Published As

Publication number Publication date
CN109931055A (zh) 2019-06-25

Similar Documents

Publication Publication Date Title
He et al. The shale gas “sweet window”:“The cracked and unbroken” state of shale and its depth range
CN105468886B (zh) 基于岩石物理学参数计算地层压力的方法
Rojstaczer et al. The influence of formation material properties on the response of water levels in wells to Earth tides and atmospheric loading
Ren et al. Acoustic velocity and electrical resistance of hydrate bearing sediments
WO2017024523A1 (zh) 一种射线弹性参数的反演方法
CN107784150B (zh) 一种裂缝封闭性预测方法
CN108089227B (zh) 新的基于三维地震资料的地层孔隙压力预测方法
CN109931054A (zh) 致密砂岩储层压力的预测方法
CN109138975B (zh) 一种基于时移测井数据的求解相渗特征曲线的新方法
WO2012035036A1 (en) Method of predicting the pressure sensitivity of seismic velocity within reservoir rocks
CN109143358B (zh) 一种获得深层、超深层碎屑岩地层压力结构的方法及装置
CN109931055B (zh) 盆地深层复合成因的流体压力预测方法
CN106289964B (zh) 一种泥页岩古抗压强度确定方法
CN106979006A (zh) 地层压力的确定方法和装置
Dai et al. Real-time ultrasonic features and damage characterization of deep shale
CN110187391B (zh) 一种利用声波时差确定地层层序界面的方法
AU2015275302B2 (en) Method and apparatus of determining stiffness coefficients of formation
Hansont et al. Fracture orientation analysis by the solid earth tidal strain method
Sun et al. Characteristics and prediction of weathered volcanic rock reservoirs: A case study of Carboniferous rocks in Zhongguai paleouplift of Junggar Basin, China
Olson et al. Estimating natural fracture producibility in tight gas sandstones: Coupling diagenesis with geomechanical modeling
CN108708715A (zh) 快速预测页岩气层地层破裂压力梯度的测录井方法
CN110671102B (zh) 一种气井临界出砂压差的确定方法及系统
Fischer et al. The creation of an asymmetric hydraulic fracture as a result of driving stress gradients
Wangen Generation of overpressure by cementation of pore space in sedimentary rocks
CN115099014B (zh) 一种基于随钻录井的天然气井地质探明储量计算方法

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

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

CB03 Change of inventor or designer information
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200925

CF01 Termination of patent right due to non-payment of annual fee