CN103375161B - 监测孔隙型地层压力梯度的方法 - Google Patents

监测孔隙型地层压力梯度的方法 Download PDF

Info

Publication number
CN103375161B
CN103375161B CN201210114213.8A CN201210114213A CN103375161B CN 103375161 B CN103375161 B CN 103375161B CN 201210114213 A CN201210114213 A CN 201210114213A CN 103375161 B CN103375161 B CN 103375161B
Authority
CN
China
Prior art keywords
gas
formation
gradient
water layer
stratum
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.)
Active
Application number
CN201210114213.8A
Other languages
English (en)
Other versions
CN103375161A (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 Petroleum and Chemical Corp
Sinopec Research Institute of Petroleum Engineering
Original Assignee
China Petroleum and Chemical Corp
Sinopec Research Institute of Petroleum Engineering
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 Petroleum and Chemical Corp, Sinopec Research Institute of Petroleum Engineering filed Critical China Petroleum and Chemical Corp
Priority to CN201210114213.8A priority Critical patent/CN103375161B/zh
Publication of CN103375161A publication Critical patent/CN103375161A/zh
Application granted granted Critical
Publication of CN103375161B publication Critical patent/CN103375161B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明涉及一种监测孔隙型碳酸盐岩地层压力梯度方法。其包括以下步骤:步骤1:建立孔隙型碳酸盐岩地层的dc指数曲线和其趋势线dcn;步骤2:在钻井过程中,判断地层为气层还是水层;步骤3:当地层为气层时,建立出口流量O曲线和其趋势线On,并计算气层的压力梯度;当地层为水层时,则建立钻井液的出口电导率和入口电导率的比值(Cout/Cin)曲线及其趋势线(Cout/Cin)n,并计算水层的压力梯度,其中,Gn为静水压力梯度,Gb为流体浮力作用产生的压力梯度。

Description

监测孔隙型地层压力梯度的方法
技术领域
本发明涉及一种监测孔隙型地层压力梯度的方法,特别是一种监测孔隙型碳酸盐岩地层压力梯度的方法。
背景技术
在钻井过程中,异常压力特别是异常高压对钻井安全极具危害作用,因此受到石油工程界的高度重视。如不能准确监测异常压力梯度,则在钻井过程中有可能发生井涌、井喷等复杂事故。
目前,异常压力梯度的监测模型包括钻前地震预测、随钻录井监测、钻后测井评价等方法,这些方法基本上都是基于碎屑岩地层的欠压实原理,即随着地层埋深的增大,孔隙度呈指数式减少。通常,地震层速度、录井dc指数、测井声波时差等曲线的趋势线表现为一条随深度增加而增大的直线。当遇到欠压实地层时,由于孔隙度变大,相应的曲线也会产生异常变化,其变化幅度与欠压实程度即压力异常幅度呈函数关系。通常对于欠压实地层可以采用等效深度法、伊顿法、比率法等进行地层压力的定量计算。
对于碳酸盐岩地层,在录井上基本普遍采用Sigma指数法进行异常压力梯度随钻监测。碳酸盐岩地层的沉积以化学胶结作用为主,所以正常压力地层的dc指数、Sigma指数随埋深的增大并不增大,其“趋势线”为一垂线(见图1),这说明上覆应力并未作用于地层流体上,因此等效深度法和伊顿法并不适用。对于孔隙型碳酸盐岩地层,dc指数和Sigma指数会降低,偏离“趋势线”。某些地区(如川东北地区)的孔隙型碳酸盐岩地层的孔隙是由于硫酸盐热化学还原反应(TSR)而产生的增容作用,这是一种降压机制,所以比率法不完全适用。另外,孔隙型碳酸盐岩地层的增压机制有烃类裂解、流体充注和构造挤压作用。因此,通常所使用的任何一种单一的监测方法均不足以反应孔隙型碳酸盐岩地层的压力梯度的变化。
由此可见,孔隙型碳酸盐岩地层以化学胶结作用为主,其地层压力既有增压机制又有降压机制、既有自源成因例如烃类裂解、TSR反应,又有它源成因例如流体充注、构造挤压。由于孔隙型碳酸盐岩地层的这种复杂性,难以对其进行异常压力的随钻准确监测。
发明内容
针对现有技术中所存在的上述技术问题,本发明提出了一种监测孔隙型碳酸盐岩地层压力梯度的方法,能够较为准确地监测孔隙型碳酸盐岩地层。
根据本发明,提出了一种监测孔隙型碳酸盐岩地层压力梯度的方法,其包括以下步骤:
步骤1:建立孔隙型碳酸盐岩地层的dc指数曲线和其趋势线dcn
步骤2:在钻井过程中,判断地层为气层还是水层;
步骤3:当为气层时,建立钻井液出口流量O的曲线和其趋势线On,并通过式(A)来计算压力梯度
G p g = G n × d c dc n × O O n + G b 式(A),
取压力梯度Gpg的最大值作为该地层段的压力梯度;
当为水层时,建立钻井液出口电导率和入口电导率的比值(Cout/Cin)曲线和其趋势线(Cout/Cin)n,并通过式(B)来计算压力梯度
G p w = G n × d c dc n × ( C o u t / C i n ) n ( C o u t / C i n ) + G b 式(B),
取压力梯度Gpw的最大值作为该地层段的压力梯度;
其中,Gn为静水压力梯度,Gb为流体浮力作用产生的压力梯度。
根据本发明的方法,在计算孔隙型碳酸盐岩气层的压力梯度Gpg时,充分考虑了各种因素对地层压力梯度的影响,因此与传统的Sigma指数相比能够较为准确地监测气层压力梯度。例如,在本发明的方法中,Gn代表了静水压力梯度的作用,Gb代表了气体浮力的作用,O/On代表了烃类裂解/充注增压机制的作用,dc/dcn代表了由于TSR反应产生的降压机制的作用,整个比值法计算方法则考虑了化学胶结作用的影响。因此根据本发明的方法能较为精确地监测孔隙型碳酸盐岩气层的压力梯度。
根据本发明的方法,在计算孔隙型碳酸盐岩水层的压力梯度Gpw时,也充分考虑了各种因素对地层压力梯度的影响,因此与传统的Sigma指数相比也能够较为准确地监测水层压力梯度。例如,在本发明的方法中,Gn代表了静水压力梯度的作用,Gb代表了地层水浮力的作用,dc/dcn代表了由于TSR反应产生的降压机制的作用,代表了地层水矿化度变化的作用,整个比值法计算过程则代表了化学胶结作用的影响。因此本发明的方法能较为精确地监测孔隙型碳酸盐岩水层的压力梯度。
根据本发明的方法,还包括在步骤1之前的预备步骤:检测碳酸盐岩样品的储集空间类型。所述碳酸盐岩样品取自正钻井的地层,以决定在随后的计算过程中,本方法是否适用。
根据本发明的方法,在钻井过程中,需要判断地层为气层还是水层,以便于使用相应的地层压力梯度计算公式来计算地层压力梯度。判断某一地层为气层还是水层对于本领域的技术人员是已知的。例如,当气测见良好显示时,为气层;当钻井液出口电导率与入口电导率的比值明显下降时,为水层。
根据本发明的方法,对于同一钻井地区,Gn和Gb为固定值。在一个实施例中,对于气层,Gn=1.02g/cm3,Gb=0.05g/cm3。对于水层,Gn=1.02g/cm3,Gb=0.2g/cm3
根据本发明的方法,在钻井过程中,dc指数、O和(Cout/Cin)为每一米或半米记录一次值。
根据本发明的方法,dcn根据下面公式建立:dcn=a×H+b,其中,a为dc指数回归方程的斜率,b为dc指数回归方程的截距,H为井深。对于不同的地层,由于dc指数回归方程不同,因此a和b也会不同。在一个气层实施例中,a=0,b=1.13,因此dcn=1.13,即dcn表现为一条不随钻井深度变化的直线;在一个水层实施例中,dcn=0.9。
根据本发明的方法,(Cout/Cin)n和On使用与dcn的建立方法相同的方法来建立。在一个气层实施例中,On=36%;在一个水层实施例中,(Cout/Cin)n=1.0。
与目前广泛使用的Sigma指数监测孔隙型碳酸盐岩地层压力梯度的方法相比,本发明的方法综合考虑了对地层压力梯度的各种影响因素。对于气层,这些影响因素包括静水压力梯度的作用、化学胶结作用、TSR反应降压机制的作用、烃类裂解/充注增压机制的作用以及地层流体的浮力作用;对于水层,这些影响因素包括除了静水压力梯度的作用、化学胶结作用、地层流体的浮力作用和TSR反应降压机制的作用外,还包括地层水的矿化度变化的作用,因而本发明的方法对于孔隙型碳酸盐岩气层和水层均能够较为准确地监测其压力梯度。
附图说明
下面结合附图来对本发明作进一步详细说明,其中:
图1是正常地层的dc指数、Sigma指数曲线图;
图2是实施本发明的方法的流程图;
图3是根据第一实施例的碳酸盐岩样品的核磁共振T2弛豫谱;
图4是根据第一实施例的测量参数的曲线图;
图5是根据第二实施例的碳酸盐岩样品的核磁共振T2弛豫谱;
图6是根据第二实施例的测量参数的曲线图。
具体实施方式
下面将结合附图对本发明作进一步说明。
图2是实施本发明的方法的流程,其用于测定孔隙型碳酸盐岩地层的压力梯度。所述方法包括以下步骤:
步骤11:建立孔隙型碳酸盐岩地层的dc指数曲线和其趋势线dcn
步骤12:判断地层为气层/水层;
步骤13:当为气层时,建立钻井液出口流量O的曲线和其趋势线On
步骤14:计算气层的压力梯度,并取其最大值作为该地层段的压力梯度;
步骤15:当为水层时,则建立钻井液出口电导率和入口电导率的比值(Cout/Cin)曲线和其趋势线(Cout/Cin)n
步骤16:计算水层的压力梯度,并取其最大值作为该地层段的压力梯度。
可以理解地是,在步骤11之前,还可具有预备步骤:检测碳酸盐岩样品的储集空间类型。在一个实施例中,通过岩石的核磁共振T2弛豫谱来判断岩石的储集空间类型。碳酸盐岩的储集空间通常分为孔隙和裂缝,其在核磁共振T2弛豫谱分别表示为具有不同弛豫时间和不同强度的峰,例如T2<100ms的峰为孔隙峰,T2>100ms的峰为裂缝峰。当孔隙峰的强度或面积占优时,则认为地层为孔隙型碳酸盐岩地层;而裂缝峰的强度或面积占优时,则认为地层为裂缝型碳酸盐岩地层。
实施例1:
选取M1井5505m深度的碳酸盐岩样品,使用核磁共振T2弛豫谱来检测其储集空间的类型,如图3所示。从图3中可以看出,碳酸盐岩样品的储集空间为孔隙型,因此该地层为孔隙型碳酸盐岩地层。
根据本发明的方法,每一米记录一次录井参数的数值,在5505-5592m深度的地层中,气测见良好显示,因此该地层为气层,如表1所示(在表1中,仅选择了计算气层压力梯度所需要的部分录井参数),由此M1井在5505-5592m深度的地层为孔隙型碳酸盐岩气层。图4以曲线方式形象地显示了录井参数的变换情况。
通过将dc指数回归计算,得出dc趋势线dcn取值为1.13。同样钻井液出口流量O的趋势线On=36%,如在表1和图4所示。由于该地层为气层,所以Gn=1.02g/cm3,Gb=0.05g/cm3
使用本发明的方法来计算该地层的压力梯度Gpg其中Gpg最大值为1.06g/cm3,如表1所示,因此使用本发明的方法判断该气层的压力梯度Gpg为1.06g/cm3。使用Sigma指数的伊顿法计算的该地层的压力梯度(FPSigma)为2.31g/cm3。而该气层的实测压力梯度为1.06g/cm3。由此可见,与使用Sigma指数的伊顿法相比,本发明的方法能更精确地计算气层的压力梯度。
表1
井深 dc dcn dc/dcn O On O/On Gpg FP Sigma
5505 0.95 1.13 0.84 36.41 36.00 1.01 0.92 1.38
5506 0.73 1.13 0.65 39.51 36.00 1.10 0.77 1.30
5507 0.95 1.13 0.84 39.55 36.00 1.10 0.99 2.18
5508 0.89 1.13 0.79 39.47 36.00 1.10 0.93 1.53
5509 0.92 1.13 0.81 39.54 36.00 1.10 0.96 2.28
5510 0.94 1.13 0.83 39.29 36.00 1.09 0.98 2.07
5511 0.92 1.13 0.81 39.17 36.00 1.09 0.95 2.29
5512 0.96 1.13 0.85 39.10 36.00 1.09 0.99 2.11
5513 0.98 1.13 0.87 38.66 36.00 1.07 1.00 2.09
5514 0.91 1.13 0.81 38.25 36.00 1.06 0.92 2.31
5515 0.96 1.13 0.85 38.98 36.00 1.08 0.99 2.19
5516 1.03 1.13 0.91 39.02 36.00 1.08 1.06 1.89
5517 1.01 1.13 0.89 39.04 36.00 1.08 1.04 1.90
5518 1.02 1.13 0.90 39.20 36.00 1.09 1.05 1.68
5519 0.99 1.13 0.88 39.17 36.00 1.09 1.02 1.66
5520 0.97 1.13 0.86 39.27 36.00 1.09 1.01 1.65
5521 1.03 1.13 0.91 39.20 36.00 1.09 1.06 1.73
5522 1.00 1.13 0.88 39.15 36.00 1.09 1.03 1.62
5523 1.02 1.13 0.90 39.11 36.00 1.09 1.05 1.75
5524 1.02 1.13 0.90 39.14 36.00 1.09 1.05 1.76
5525 0.98 1.13 0.87 39.63 36.00 1.10 1.02 1.73
5526 1.01 1.13 0.89 39.88 36.00 1.11 1.06 1.70
5527 1.04 1.13 0.92 38.74 36.00 1.08 1.06 1.68
实施例2:
选取M2井6750m深度的碳酸盐岩样品,使用核磁共振T2弛豫谱来检测其储集空间的类型,如图5所示。从图5中可以看出,碳酸盐岩样品的储集空间为孔隙型,因此该地层为孔隙型碳酸盐岩地层。
根据本发明的方法,每一米记录一次录井参数的数值,在6750-7000m深度的地层中,钻井液出口电导率与入口电导率的比值明显下降,左偏离趋势线,因此该地层为水层,如表2所示(在表2中,仅选择了计算水层压力梯度所需要的部分录井参数),由此M2井在6750-7000m深度的地层为孔隙型碳酸盐岩水层。图6以曲线方式形象地显示了录井参数的变换情况。
通过将dc指数回归计算,得出dc趋势线dcn取值为0.9。同样钻井液出口电导率与入口电导率的比值(Cout/Cin)的趋势线(Cout/Cin)n=1.0,如在表2和图6所示。由于该地层为水层,所以Gn=1.02g/cm3,Gb=0.2g/cm3
使用本发明的方法计算该地层的压力梯度Gpw其中Gpw最大值为1.12g/cm3,如表2所示,因此使用本发明的方法判断该水层的压力梯度Gpw为1.12g/cm3。使用Sigma指数的伊顿法计算的该地层的压力梯度(FPSigma)为2.78g/cm3。而实测该水层的压力梯度为1.15g/cm3。由此可见,与使用Sigma指数的伊顿法相比,本发明的方法能更精确地计算水层的压力梯度。
表2
根据本发明的方法能较为精确地计算孔隙型碳酸盐岩地层的压力梯度,为保证钻井安全提供有力的技术支持。
虽然已经参考优选实施例对本发明进行了描述,但在不脱离本发明的范围的情况下,可以对其进行各种改进。本发明并不局限于文中公开的实施例,而是包括落入权利要求的范围内的所有技术方案。

Claims (10)

1.一种监测孔隙型碳酸盐岩地层压力梯度的方法,包括以下步骤:
步骤1:建立孔隙型碳酸盐岩地层的dc指数曲线和其趋势线dcn
步骤2:在钻井过程中,判断地层为气层还是水层;
步骤3:当为气层时,建立钻井液出口流量O的曲线和其趋势线On,并通过式(A)来计算压力梯度
G p g = G n &times; d c dc n &times; O O n + G b 式(A),
取压力梯度Gpg的最大值作为该地层段的压力梯度;
当为水层时,则建立钻井液出口电导率和入口电导率的比值(Cout/Cin)曲线和其趋势线(Cout/Cin)n,并通过式(B)来计算压力梯度
G p w = G n &times; d c dc n &times; ( C o u t / C i n ) n ( C o u t / C i n ) + G b 式(B),
取压力梯度Gpw的最大值作为该地层段的压力梯度;
其中,Gn为静水压力梯度,Gb为流体浮力作用产生的压力梯度。
2.根据权利要求1所述的方法,其特征在于,还包括在步骤1之前的预备步骤:检测碳酸盐岩样品的储集空间类型。
3.根据上述权利要求1或2所述的方法,其特征在于,对于同一钻井地区,所述Gn和Gb为固定值。
4.根据上述权利要求3所述的方法,其特征在于,对于气层,Gn=1.02g/cm3,Gb=0.05g/cm3
5.根据上述权利要求3所述的方法,其特征在于,对于水层,Gn=1.02g/cm3,Gb=0.2g/cm3
6.根据权利要求1或2所述的方法,其特征在于,所述dc指数、O和(Cout/Cin)为每一米或半米记录一次值。
7.根据权利要求6所述的方法,其特征在于,对于气层,dcn=1.13,On=36%。
8.根据权利要求6所述的方法,其特征在于,对于水层,dcn=0.9;(Cout/Cin)n=1.0。
9.根据权利要求1或2所述的方法,其特征在于,在钻井过程中,当气测见良好显示时,地层为气层。
10.根据权利要求1或2所述的方法,其特征在于,在钻井过程中,当钻井液出口电导率与入口电导率的比值明显下降时,地层为水层。
CN201210114213.8A 2012-04-18 2012-04-18 监测孔隙型地层压力梯度的方法 Active CN103375161B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210114213.8A CN103375161B (zh) 2012-04-18 2012-04-18 监测孔隙型地层压力梯度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210114213.8A CN103375161B (zh) 2012-04-18 2012-04-18 监测孔隙型地层压力梯度的方法

Publications (2)

Publication Number Publication Date
CN103375161A CN103375161A (zh) 2013-10-30
CN103375161B true CN103375161B (zh) 2015-12-16

Family

ID=49460940

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210114213.8A Active CN103375161B (zh) 2012-04-18 2012-04-18 监测孔隙型地层压力梯度的方法

Country Status (1)

Country Link
CN (1) CN103375161B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103821505B (zh) * 2014-02-28 2016-10-05 中国石油大学(北京) 砂岩油气输导层地球物理-地质-地球化学检测方法及装置
CN106202772B (zh) * 2016-07-18 2019-06-14 三一海洋重工有限公司 半潜平台压载水舱的设计水压头计算方法、装置及设备

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1991004500A1 (en) * 1989-09-20 1991-04-04 Chevron Research And Technology Company Pore pressure prediction method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2482457A1 (en) * 2002-04-10 2004-03-18 Schlumberger Technology Corporation Method, apparatus and system for pore pressure prediction in presence of dipping formations

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1991004500A1 (en) * 1989-09-20 1991-04-04 Chevron Research And Technology Company Pore pressure prediction method
CN1052530A (zh) * 1989-09-20 1991-06-26 切夫里昂研究和技术公司 空隙压力预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
异常高压随钻预测理论与方法;王志战,路黄生;《录井工程》;20110930;第22卷(第3期);全文 *

Also Published As

Publication number Publication date
CN103375161A (zh) 2013-10-30

Similar Documents

Publication Publication Date Title
CN106285646B (zh) 基于多信息融合的钻井漏失层位识别方法
CN104863577B (zh) 利用地震纵波传播时间预测地层孔隙压力的方法
CN105607146B (zh) 一种曲流河砂体规模的定量表征方法
CN106837297A (zh) 一种识别井间连通性及油水动态预测的方法
CN104695950A (zh) 火山岩油藏产能预测方法
CN106761677A (zh) 页岩气水平井单井产能的测井预测方法
CN106401574B (zh) 一种钻前高温地热井地层压力的预测方法
CN102707333A (zh) 页岩气资源量/储量的测量方法
CN104500054B (zh) 地层孔隙压力的确定方法及装置
CN104114811A (zh) 定位诱导的地下岩层裂缝中的示踪支撑剂的与岩性和井眼状况无关的方法
CN104379870A (zh) 用于确定储层的径向流响应的渗透率或迁移率的方法
CN105652329A (zh) 一种评估煤层顶板视涌水量的方法和装置
CN108301825A (zh) 一种高温高压储层孔隙压力分析方法
CN109254138A (zh) 一种基于砂岩微观特征的富水性评价方法
CN112083515B (zh) 致密砂岩低阻储层挖掘效应的定量表征及含气性评价方法
Zhou et al. Sensitivity analysis in permeability estimation using logging and injection-falloff test data for an anthracite coalbed methane reservoir in Southeast Qinshui Basin, China
CN106368686A (zh) 基于岩石泊松比计算地层压力的方法
Odling et al. Fractional flow in fractured chalk; a flow and tracer test revisited
CN112363242A (zh) 基于测录井融合的储层流体识别方法与装置
CN107515430B (zh) 一种地震法探测盐湖卤水的方法
CN101936157B (zh) 一种利用测井资料检测高压盐水层孔隙压力的方法
Izgec Understanding waterflood performance with modern analytical techniques
CN103375161B (zh) 监测孔隙型地层压力梯度的方法
CN104790943A (zh) 一种油气储层含油性与孔隙性综合指数的计算方法
CN103375162B (zh) 监测裂缝型地层压力梯度的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant