CN111927445A - 基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法 - Google Patents

基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法 Download PDF

Info

Publication number
CN111927445A
CN111927445A CN202010925184.8A CN202010925184A CN111927445A CN 111927445 A CN111927445 A CN 111927445A CN 202010925184 A CN202010925184 A CN 202010925184A CN 111927445 A CN111927445 A CN 111927445A
Authority
CN
China
Prior art keywords
drilling
illite
montmorillonite
pressure
formation pressure
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
CN202010925184.8A
Other languages
English (en)
Other versions
CN111927445B (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.)
China National Petroleum Corp
CNPC Xibu Drilling Engineering Co Ltd
Original Assignee
China National Petroleum Corp
CNPC Xibu Drilling Engineering Co 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 China National Petroleum Corp, CNPC Xibu Drilling Engineering Co Ltd filed Critical China National Petroleum Corp
Priority to CN202010925184.8A priority Critical patent/CN111927445B/zh
Publication of CN111927445A publication Critical patent/CN111927445A/zh
Application granted granted Critical
Publication of CN111927445B publication Critical patent/CN111927445B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/06Measuring temperature or pressure
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geophysics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Earth Drilling (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及随钻地层压力计算方法技术领域,是一种基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法,其先根据随钻岩矿数据计算随钻地层的岩石密度、泥质含量、垂向有效应力、声波速度和孔隙度,再利用计算得到的各项参数计算随钻地层压力。本发明所述获取随钻地层压力的方法利用连续、及时的获取的随钻岩矿分析数据,结合测井不受欠压实成因限制的解释方法,建立了多种压力成因机制下的地层压力计算方法,且本发明所述获取随钻地层压力的方法计算的随钻地层压力与测井解释地层压力相关性极高,满足了现场对地层压力的评价,进而可以指导钻井液密度的使用,开展地层压力风险预警。

Description

基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法
技术领域
本发明涉及随钻地层压力计算方法技术领域,是一种基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法。
背景技术
国内有关随钻地层压力的解释方法,主要基于测井数据或综合录井数据,尤其是地层压力成因复杂地区,还没有使用随钻岩矿数据计算随钻地层压力的方法或成果。
准噶尔盆地南缘地区地质条件复杂,经历了多期构造运动。钻探难度大,且随钻测井费用较高,近年来部署的井,均未开展随钻测井。目前应用测井资料对地层压力进行的评价通常采用的是钻后电缆测井资料,但是钻后评价无法实时指导钻井工作的进行、保证施工安全。
申请号为2018109242273的中国专利申请文献提供了一种基于综合录井参数的随钻压力预测方法,该方法通过获取目标井的第一综合录井参数,根据所述第一综合录井参数和压力预测模型,预测所述目标井中对应深度的地层的压力系数。其中,压力预测模型是基于已钻井的压力系数及第二综合录井参数,结合DNN模型训练得到,已钻井与目标井的地质情况相当。由于综合录井参数数据量大的特点,根据已钻井的压力系数及第二综合录井参数训练得到的压力预测模型的预测准确度较高;又由于已钻井与目标井的地质情况相当,因此,可根据该压力预测模型和目标井对应的第一综合录井参数,精确预测目标井中对应深度的地层的压力系数。
而上述基于综合录井参数的随钻压力预测方法,无法从地层压力的成因出发对地层压力进行针对性计算。
樊洪海(《异常地层压力分析方法与应用》,201603)在一书中指出声波在地层中的传播速度与孔隙度及垂直有效应力有明显的关系,尤其是泥岩地层更是如此。目前绝大多数测井资料评价地层压力的方法都是声波测井。而且补偿声波测井减少了井径变化的影响,提高了声波测井资料的精度。
邢培俊等人(邢培俊等.利用测井资料确定粘土矿物的方法对比.中国石油大学学报(自然科学版)2008,第32卷第2期,53-57)利用蒙脱石、伊利石等粘土矿物在声波等测井资料上的差异,建立了声波等测井资料与粘土矿物的关系。
发明内容
本发明提供了一种基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法,克服了上述现有技术之不足,其能有效解决现有随钻压力预测方法无法从地层压力的成因机制出发以获取随钻地层压力的问题。
本发明的技术方案是通过以下措施来实现的:一种基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法,包括:
采用公式(1)计算声波时差,
Figure BDA0002668227320000021
公式(1)中:
FAC为声波时差,μs/ft,
RIS为随钻伊/蒙混层相对含量,
b0、b1、b2、b3为与伊/蒙混层相关的系数;
采用公式(2)计算声波速度,
vp=1/FAC (2)
公式(2)中:
vp为声波速度,km/s;
采用公式(3)计算岩石密度,采用公式(4)计算泥质含量,采用公式(5)计算垂向有效应力,
Figure BDA0002668227320000022
Vsh=RIS+RM+RK+RI (4)
Figure BDA0002668227320000023
公式(3)至(5)中:
ρ为岩石密度,kg/m3,
vp为声波速度,km/s,
Vsh为泥质含量,%,
σev -为垂向有效应力,MPa,
pe为地层压力,MPa,
Figure BDA0002668227320000024
为高泉背斜已钻井对应地层的平均孔隙度,%,
h为井深,
g为重力加速度,
RQ、RIS、RA、RM、RK、RI、RD、RC、RAN、RR、RB、RCU、RN分别为石英、伊/蒙混层、斜长石、蒙脱石、高岭土、伊利石、白云石、方解石、硬石膏、石盐、重晶石、方沸石及其他矿物的相对含量,
DQ、DIS、DA、DM、DK、DI、DD、DC、DAN、DR、DB、DCU、DO分别为石英、伊/蒙混层、斜长石、蒙脱石、高岭土、伊利石、白云石、方解石、硬石膏、石盐、重晶石、方沸石及其他矿物的标准密度;
采用公式(6)计算地层压力,
Figure BDA0002668227320000031
公式(6)中:
ρ为岩石密度,kg/m3,
vp为声波速度,km/s
Vsh为泥质含量,%,
σev -为垂向有效应力,MPa,
pe为地层压力,MPa,
Figure BDA0002668227320000032
为高泉背斜已钻井对应地层的平均孔隙度,%。
下面是对上述发明技术方案的进一步优化或/和改进:
上述公式(3)中,石英、伊/蒙混层、斜长石、蒙脱石、高岭土、伊利石、白云石、方解石、硬石膏、石盐、重晶石、方沸石及其他矿物的相对含量采用现有X-衍射全岩分析对随钻采集的岩屑进行矿物含量分析而获得。
本发明所述获取随钻地层压力的方法利用连续、及时的获取的随钻岩矿分析数据,结合测井不受欠压实成因限制的解释方法,建立了多种压力成因机制下的地层压力计算方法,且本发明所述获取随钻地层压力的方法计算的随钻地层压力与测井解释地层压力相关性高达93.5%,满足了现场对地层压力的评价,进而可以指导钻井液密度的使用,开展地层压力风险预警。
附图说明
附图1为本发明所述方法的主要流程。
具体实施方式
本发明不受下述实施例的限制,可根据本发明的技术方案与实际情况来确定具体的实施方式。
下面结合实施例对本发明作进一步描述:
实施例1:如附图1所示,该基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法,包括:
采用公式(1)计算声波时差,
Figure BDA0002668227320000033
公式(1)中:
FAC为声波时差,μs/ft,
RIS为随钻伊/蒙混层相对含量,
b0、b1、b2、b3为与伊/蒙混层相关的系数;
采用公式(2)计算声波速度,
vp=1/FAC (2)
公式(2)中:
vp为声波速度,km/s;
采用公式(3)计算岩石密度,采用公式(4)计算泥质含量,采用公式(5)计算垂向有效应力,
Figure BDA0002668227320000041
Vsh=RIS+RM+RK+RI (4)
Figure BDA0002668227320000042
公式(3)至(5)中:
ρ为岩石密度,kg/m3,
vp为声波速度,km/s,
Vsh为泥质含量,%,
σev -为垂向有效应力,MPa,
pe为地层压力,MPa,
Figure BDA0002668227320000043
为高泉背斜已钻井对应地层的平均孔隙度(所述地层的泥岩孔隙度、渗透率较低,因此取平均值),%,
h为井深,
g为重力加速度,
RQ、RIS、RA、RM、RK、RI、RD、RC、RAN、RR、RB、RCU、RN分别为石英、伊/蒙混层、斜长石、蒙脱石、高岭土、伊利石、白云石、方解石、硬石膏、石盐、重晶石、方沸石及其他矿物的相对含量,
DQ、DIS、DA、DM、DK、DI、DD、DC、DAN、DR、DB、DCU、DO分别为石英、伊/蒙混层、斜长石、蒙脱石、高岭土、伊利石、白云石、方解石、硬石膏、石盐、重晶石、方沸石及其他矿物的标准密度;
采用公式(6)计算地层压力,
Figure BDA0002668227320000044
公式(6)中:
ρ为岩石密度,kg/m3,
vp为声波速度,km/s,
Vsh为泥质含量,%,
σev -为垂向有效应力,MPa,
pe为地层压力,MPa,
Figure BDA0002668227320000051
为高泉背斜已钻井对应地层的平均孔隙度,%。
蒙脱石的晶体结构是以水合阳离子及水分子作为层间物的2:1型黏土矿物。在埋藏后,随压力的增加,层间水将逐渐被释放出来,造成层间塌陷,形成伊利石/蒙脱石混层(伊/蒙混层),进而逐渐向伊利石转化。伊利石/蒙脱石混层是蒙脱石向伊利石转化的中间产物。在压力等多种因素的影响下,由层间水移出而导致的层间塌陷,使得蒙脱石晶体结构中形成少量伊利石的结构层,它们杂乱或有规律分布于蒙脱石结构层之间。在转化作用早期,伊/蒙混层中的伊利石层含量较少,在蒙脱石层之间无规律分布;随转化作用的增强,伊利石层的相对含量逐渐升高,分布逐步有规律,最终蒙脱石层消失,伊/蒙混层转化为伊利石。但这个反应过程受压力的影响,甚至会因超压而停止。因此可以依据伊蒙混层的相对含量变化,建立岩矿数据与地层压力的关系。
本发明所述获取随钻地层压力的方法中,有两个关键点:
关键点一:通过公式(1)获取随钻声波时差数据,将多口水平井的测井声波时差数据与通过公式(1)获取的随钻声波时差数据验证对比,表明采用公式(1)计算声波时差的精度达到90%以上。
关键点二:通过公式(2)至(5)分别求取各项参数(vp、ρ、Vsh、σev -
Figure BDA0002668227320000052
),再利用公式(6)计算随钻地层压力(pe)数据,通过与测井数据计算的随钻地层压力数据的验证对比,表明采用本发明所述获取随钻地层压力的方法的计算精度达到90%以上。
采用本发明所述获取随钻地层压力的方法计算的随钻地层压力与测井解释地层压力相关性高达93.5%,满足了现场对地层压力的评价,进而可以指导钻井液密度的使用,开展地层压力风险预警。
本发明所述获取随钻地层压力的方法利用连续、及时的获取的随钻岩矿分析数据,结合测井不受欠压实成因限制的解释方法(公式6),建立了多种压力成因机制下的地层压力计算方法,即本发明所述获取随钻地层压力的方法计是一种不受欠压实成因限制的随钻地层压力获取方法。
公式6为适合准噶尔盆地南缘高泉背斜的地层压力测井多参数模型。
公式(3)中,石英、伊/蒙混层、斜长石、蒙脱石、高岭土、伊利石、白云石、方解石、硬石膏、石盐、重晶石、方沸石及其他矿物的相对含量可采用现有X-衍射全岩分析(XRD)技术对随钻采集的岩屑进行矿物含量分析而获得。
这里所述的其他矿物为随钻采集的岩屑中,除了石英、伊/蒙混层、斜长石、蒙脱石、高岭土、伊利石、白云石、方解石、硬石膏、石盐、重晶石、方沸石以外的其他矿物。
以上技术特征构成了本发明的实施例,其具有较强的适应性和实施效果,可根据实际需要增减非必要的技术特征,来满足不同情况的需求。

Claims (2)

1.一种基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法,其特征在于包括:
采用公式(1)计算声波时差,
Figure FDA0002668227310000011
公式(1)中:
FAC为声波时差,μs/ft,
RIS为随钻伊/蒙混层相对含量,
b0、b1、b2、b3为与伊/蒙混层相关的系数;
采用公式(2)计算声波速度,
vp=1/FAC (2)
公式(2)中:
vp为声波速度,km/s;
采用公式(3)计算岩石密度,采用公式(4)计算泥质含量,采用公式(5)计算垂向有
Figure FDA0002668227310000012
Vsh=RIS+RM+RK+RI (4)
Figure FDA0002668227310000013
公式(3)至(5)中:
ρ为岩石密度,kg/m3,
vp为声波速度,km/s,
Vsh为泥质含量,%,
σev为垂向有效应力,MPa,
pe为地层压力,MPa,
Figure FDA0002668227310000014
为高泉背斜已钻井对应地层的平均孔隙度,%,
h为井深,
g为重力加速度,
RQ、RIS、RA、RM、RK、RI、RD、RC、RAN、RR、RB、RCU、RN分别为石英、伊/蒙混层、斜长石、蒙脱石、高岭土、伊利石、白云石、方解石、硬石膏、石盐、重晶石、方沸石及其他矿物的相对含量,
DQ、DIS、DA、DM、DK、DI、DD、DC、DAN、DR、DB、DCU、DO分别为石英、伊/蒙混层、斜长石、蒙脱石、高岭土、伊利石、白云石、方解石、硬石膏、石盐、重晶石、方沸石及其他矿物的标准密度;
采用公式(6)计算地层压力,
Figure FDA0002668227310000015
公式(6)中:
ρ为岩石密度,kg/m3,
vp为声波速度,km/s
Vsh为泥质含量,%,
σev为垂向有效应力,MPa,
pe为地层压力,MPa,
Figure FDA0002668227310000021
为高泉背斜已钻井对应地层的平均孔隙度,%。
2.根据权利要求1所述基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法,其特征在于公式(3)中,石英、伊/蒙混层、斜长石、蒙脱石、高岭土、伊利石、白云石、方解石、硬石膏、石盐、重晶石、方沸石及其他矿物的相对含量采用现有X-衍射全岩分析对随钻采集的岩屑进行矿物含量分析而获得。
CN202010925184.8A 2020-09-06 2020-09-06 基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法 Active CN111927445B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010925184.8A CN111927445B (zh) 2020-09-06 2020-09-06 基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010925184.8A CN111927445B (zh) 2020-09-06 2020-09-06 基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法

Publications (2)

Publication Number Publication Date
CN111927445A true CN111927445A (zh) 2020-11-13
CN111927445B CN111927445B (zh) 2023-10-27

Family

ID=73308610

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010925184.8A Active CN111927445B (zh) 2020-09-06 2020-09-06 基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法

Country Status (1)

Country Link
CN (1) CN111927445B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112664186A (zh) * 2020-12-25 2021-04-16 中法渤海地质服务有限公司 一种利用甲烷碳同位素资料预警和识别他源高压的方法
CN114033362A (zh) * 2021-10-28 2022-02-11 成都亿页图蓝科技有限公司 一种基于录井和测井数据的页岩气声波曲线重建方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150267527A1 (en) * 2014-03-21 2015-09-24 Conocophillips Company Method for analysing pore pressure in shale formations
CN106979006A (zh) * 2017-05-17 2017-07-25 中国神华能源股份有限公司 地层压力的确定方法和装置
CN109236286A (zh) * 2018-10-23 2019-01-18 西南石油大学 一种新型地层孔隙压力随钻测量的方法
CN111506863A (zh) * 2020-05-29 2020-08-07 中国石油天然气股份有限公司西南油气田分公司川中油气矿 一种自动计算气藏地层压力的方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150267527A1 (en) * 2014-03-21 2015-09-24 Conocophillips Company Method for analysing pore pressure in shale formations
CN106979006A (zh) * 2017-05-17 2017-07-25 中国神华能源股份有限公司 地层压力的确定方法和装置
CN109236286A (zh) * 2018-10-23 2019-01-18 西南石油大学 一种新型地层孔隙压力随钻测量的方法
CN111506863A (zh) * 2020-05-29 2020-08-07 中国石油天然气股份有限公司西南油气田分公司川中油气矿 一种自动计算气藏地层压力的方法及系统

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112664186A (zh) * 2020-12-25 2021-04-16 中法渤海地质服务有限公司 一种利用甲烷碳同位素资料预警和识别他源高压的方法
CN114033362A (zh) * 2021-10-28 2022-02-11 成都亿页图蓝科技有限公司 一种基于录井和测井数据的页岩气声波曲线重建方法
CN114033362B (zh) * 2021-10-28 2024-04-26 成都亿页图蓝科技有限公司 一种基于录井和测井数据的页岩气声波曲线重建方法

Also Published As

Publication number Publication date
CN111927445B (zh) 2023-10-27

Similar Documents

Publication Publication Date Title
CN107766662B (zh) 一种页岩气水平井测试分段评价方法
CN106501856B (zh) 亚地震断层定量预测方法
CN108301825B (zh) 一种高温高压储层孔隙压力分析方法
CN109632459B (zh) 一种页岩可压裂性评价方法
EP0568643A4 (en) MANUFACTURE OF FINE FILMS FROM POLYESTER.
CN111927445A (zh) 基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法
Abbas et al. Stability analysis of highly deviated boreholes to minimize drilling risks and nonproductive time
CN112505784B (zh) 一种煤矿井下盾构机掘进路线选择方法
CN104153768A (zh) 一种评价花岗岩储层储集性能的方法
CN112560246B (zh) 一种目标井散点地层压力系数的预测方法
Lorenz et al. Subsurface fracture spacing: Comparison of inferences from slant/horizontal and vertical cores
CN114462237A (zh) 一种泥页岩油气甜点预测方法
Neeamy et al. Building 1D mechanical earth model for Zubair oilfield in Iraq
Mohammed et al. Building 1D mechanical earth model for Zubair oilfield in Iraq
CN110764161B (zh) 一种通过钻探资料判识碳酸盐岩断裂破碎带的综合方法
Amanipoor et al. The application of géostatistical methods to prepare the 3D petrophysical model of oil reservoir
Walker et al. Straining at the leash: Understanding the full potential of the deep-water, subsalt Mad Dog field, from appraisal through early production
Al-Sudani et al. Geomechanical study to predict the onset of sand production formation
CN105259576A (zh) 一种利用地震统计特征的油气藏识别方法
CN114910953B (zh) 一种复杂断裂带不规则断面的解释方法
Phipps et al. The Babbage Field, Block 48/2a, UK North Sea
Huang et al. Development characteristics and prediction of fractures in deep-ultradeep tight sandstone reservoirs of Kuqa Depression: a case study of the Bashijiqike and Baxigai formations in Dabei 9 Block
CN116148929A (zh) 井震结合上覆地层等效密度估算方法
Dou et al. Seismic prediction technology for shale gas “double sweet spots”: A case study of the Ordovician Ulalik Formation in the western margin of the Ordos Basin, China
CN117892480A (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