CN113958315B - 一种基于自吸-本构模型的页岩地层坍塌压力预测方法 - Google Patents
一种基于自吸-本构模型的页岩地层坍塌压力预测方法 Download PDFInfo
- Publication number
- CN113958315B CN113958315B CN202110668186.8A CN202110668186A CN113958315B CN 113958315 B CN113958315 B CN 113958315B CN 202110668186 A CN202110668186 A CN 202110668186A CN 113958315 B CN113958315 B CN 113958315B
- Authority
- CN
- China
- Prior art keywords
- shale
- self
- absorption
- stress
- rock
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000010521 absorption reaction Methods 0.000 claims abstract description 87
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 58
- 238000005553 drilling Methods 0.000 claims abstract description 48
- 230000006378 damage Effects 0.000 claims abstract description 47
- 239000011435 rock Substances 0.000 claims abstract description 40
- 230000008569 process Effects 0.000 claims abstract description 30
- 238000002474 experimental method Methods 0.000 claims abstract description 28
- 239000012530 fluid Substances 0.000 claims abstract description 23
- 230000036571 hydration Effects 0.000 claims abstract description 20
- 238000006703 hydration reaction Methods 0.000 claims abstract description 20
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 18
- 239000011148 porous material Substances 0.000 claims description 24
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 7
- 229910052753 mercury Inorganic materials 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 239000007788 liquid Substances 0.000 claims 1
- 230000008859 change Effects 0.000 description 4
- 239000000243 solution Substances 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241000695274 Processa Species 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing 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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- Animal Husbandry (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Marine Sciences & Fisheries (AREA)
- General Business, Economics & Management (AREA)
- Agronomy & Crop Science (AREA)
- Geophysics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种基于自吸‑本构模型的页岩地层坍塌压力预测方法,包括:基于损伤力学理论,结合页岩动态自吸方程、损伤下的页岩本构方程,构建水化损伤下的页岩自吸‑本构模型;对页岩岩心开展三轴力学实验获得原状页岩力学参数;对页岩岩心开展自吸实验,基于页岩动态自吸方程,计算出页岩平均吸水长度;获取井壁岩石主应力分布;基于井壁岩石主应力分布、原状页岩力学参数、页岩平均吸水长度,采用所述页岩自吸‑本构模型,实现钻井过程中的页岩地层坍塌压力预测。本发明结合页岩自吸特征与本构特征,构建了一种新的坍塌压力预测方法,可为现场确定安全钻井液密度提供科学依据。
Description
技术领域
本发明涉及页岩地层稳定井壁地层能力评价方法领域,尤其涉及一种基于自吸-本构模型的页岩地层坍塌压力预测方法。
背景技术
页岩气是目前石油工业的热点。由于页岩气具有巨大的潜能,世界上多个国家都逐渐增加了页岩气勘探开发的投入。页岩地层钻井过程中,井壁垮塌失稳现象严重,从而严重制约了页岩气的高效开发。实际钻井过程中,钻井液密度需要高于地层坍塌压力,从而保持井壁不垮塌。因此,需要一种准确的坍塌压力预测方法,指导钻井过程中的钻井液密度优选。
页岩的特点之一是具有强水化特征,在钻井过程中,随着页岩与钻井液接触,页岩发生水化作用,力学特征发生改变,从而影响地层坍塌压力。目前,页岩地层坍塌压力研究中,已经形成了明确认识:当页岩与钻井液接触,随水化时间增加,页岩吸水量越大,水化损伤越明显,坍塌压力增大,井壁稳定性降低。但是,尚未形成有效的理论方法建立“吸水特征—损伤程度—坍塌压力”的定量关系,从而导致实际钻井过程中,坍塌压力预测不准确,井壁失稳现象频发。
发明内容
本申请为了解决上述技术问题提供一种基于自吸-本构模型的页岩地层坍塌压力预测方法。
本申请通过下述技术方案实现:
一种基于自吸-本构模型的页岩地层坍塌压力预测方法,包括以下步骤:
基于Weibull强度分布理论和损伤力学定义,构建损伤下的页岩本构方程;
基于损伤力学理论,结合页岩动态自吸方程、损伤下的页岩本构方程,构建水化损伤下的页岩自吸-本构模型;
制取页岩岩心,对页岩岩心开展三轴力学实验获得原状页岩力学参数;
对页岩岩心开展自吸实验,基于页岩动态自吸方程,计算出页岩平均吸水长度;
获取井壁岩石主应力分布;
基于井壁岩石主应力分布、原状页岩力学参数、页岩平均吸水长度,采用所述页岩自吸-本构模型,实现钻井过程中的页岩地层坍塌压力预测。
其中,所述原状页岩力学参数包括页岩弹性模量、泊松比、内聚力、抗压强度-围压的相关性;抗压强度-围压的相关性可通过以下方式获得:通过开展不同围压下的三轴力学实验,获得不同围压下的抗压强度,从而进行线性拟合,确定抗压强度-围压的相关性方程。
进一步的,所述损伤下的页岩本构方程为:
进一步的,所述损伤系数D的表达式为:
上式中,L为页岩岩样长度,Ls为页岩平均吸水长度。
进一步的,所述的自吸实验具体包括:
S1,对页岩岩心开展自吸实验,获取不同自吸时间下的页岩岩心;
S2,对不同自吸时间下的页岩岩心开展压汞实验,获取自吸过程中的页岩平均孔径。
进一步的,所述的自吸实验采用端部自吸方式,水沿岩心一个端面进入页岩内部,吸水方向平行页岩岩心轴向方向。
进一步的,所述页岩动态自吸方程为:
上式中,Ls为页岩平均吸水长度,σ为表面张力,θ为接触角,μw为水黏度,p为水密度,λa(t)为自吸过程中的页岩平均孔径。
进一步的,所述页岩自吸-本构模型为:
Ls为页岩平均吸水长度,σ为表面张力,θ为接触角,μw为水黏度,p为水密度,λa(t)为自吸过程中的页岩平均孔径。
进一步的,井壁岩石主应力分布的获取方法具体包括:
基于地层三向地应力,采用井眼坐标转换方程,得到井周应力分布:
上式中,σr、σθ、σz为柱坐标系下径向、周向和轴向正应力;τθz、τrθ、τrz分别为柱坐标下θz、rθ、rz平面切应力;σxx、σyy、σzz分别为直角坐标下X、Y、Z轴方向上地应力分量;σxy、σxz、σyz分别为直角坐标系下XY、XZ、YZ平面上地应力分量;θc井周角;pi为钻井液液柱压力;
基于井周应力分布,得到井壁任意位置的主应力分布,如下所示:
上式中,σ1为井壁岩石最大主应力;σ3井壁岩石最小主应力。
与现有技术相比,本申请具有以下有益效果:
本发明的方法考虑了页岩自吸过程中,水化作用导致的页岩孔隙结构变化,形成了动态的页岩自吸特征表征方法,进而利用动态自吸特征建立水化损伤表征参数,最终结合页岩自吸特征与本构特征,构建了一种新的坍塌压力预测方法,可为现场确定安全钻井液密度提供科学依据。
附图说明
此处所说明的附图用来提供对本申请实施方式的进一步理解,构成本申请的一部分,并不构成对本发明实施方式的限定。
图1为本发明的流程框图;
图2是不同围压下的原状页岩本构关系图;
图3是不同自吸时间下的页岩本构关系图;
图4是页岩自吸示意图;
图5是自吸过程中的页岩平均吸水长度;
图6a是不同钻井时间下的地层坍塌压力分布图;
图6b是钻井时间为24小时的页岩应力—应变曲线图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,下面将结合实施方式中的附图,对本发明实施方式中的技术方案进行清楚、完整地描述。显然,所描述的实施方式是本发明一部分实施方式,而不是全部的实施方式。通常在此处附图中描述和示出的本发明实施方式的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施方式。基于本发明中的实施方式,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式,都属于本发明保护的范围。
需要说明的是,在不冲突的情况下,本发明中的实施方式及实施方式中的特征可以相互组合。需要说明的是,本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。
如图1所示,本发明公开的基于自吸-本构模型的页岩地层坍塌压力预测方法,包括以下步骤:
步骤A、构建构建水化损伤下的页岩自吸-本构模型。
A.1、基于Weibull强度分布理论和损伤力学定义,构建损伤下的页岩本构方程,具体如下:
依据Weibull强度分布理论,本构方程表达式为:
上式中,σi、εi分别为页岩应力、应变;E、u分别为页岩弹性模量与泊松比;为内聚力;σ3为围压或最小主应力,σc为不同围压下的抗压强度;F、m是Weibull强度统计参数,基于页岩三轴力学参数,查阅系数分布版图,确定具体数值。
由于钻井过程中,钻井液与井壁岩石接触,钻井液内部的水侵入页岩,造成水化损伤。当岩体发生损伤时,基于损伤力学定义,本申请将损伤下的页岩本构方程改写为:
通过对比可知:原状页岩本构方程与损伤下的页岩本构方程的主要差异在于损伤系数D,损伤系数用于表征损伤程度。
A.2、基于Lucas-Washburn吸水理论,考虑自吸过程中的页岩平均孔径变化,页岩动态自吸方程为:
上式中,Ls为页岩平均吸水长度,σ为表面张力,θ为接触角,μw为水黏度,p为水密度,λa(t)为自吸过程中的页岩平均孔径;其中,表面张力、接触角、水黏度、水密度为水-页岩的界面参数,可直接查阅资料获取。采用上式自吸方程,可以得到不同自吸时间t下的页岩平均吸水长度Ls。
A.3、基于损伤力学理论,结合页岩动态自吸方程、损伤下的页岩本构方程,构建水化损伤下的页岩自吸-本构模型,具体如下:
自吸方程与本构方程的结合依据损伤力学理论。本申请中涉及的损伤指水化损伤,即水与页岩相互作用产生的损伤。根据现有研究认识,当页岩吸水越显著,造成的水化损伤越强。因此,利用页岩自吸方程求取平均吸水长度,以平均吸水长度代表页岩水化损伤,损伤系数D的表达式如下:
上式中,L为页岩岩样长度,Ls为吸水长度,采用自吸方程进行计算。显然,吸水长度越大,水化损伤越强;当吸水长度等于岩样长度时,代表岩样完全吸水,损伤系数达到最大,此时D=1。
因此,最终形成考虑水化损伤的页岩自吸-本构模型,如下所示:
根据原状页岩力学参数、Weibull强度统计系数、损伤系数D,采用考虑水化损伤的页岩自吸-本构模型,不同钻井时间下,造成不同的损伤,从而形成不同的本构关系,如图3所示。
步骤B、制取页岩岩心,通过野外采集得到的页岩岩样露头,通过钻机钻取标准圆柱样岩心。
步骤C、对页岩岩心展开三轴力学实验和页岩吸水实验。
C1、三轴力学实验
对页岩岩心开展三轴力学实验,获取原状页岩力学参数。
原状页岩力学参数表示未受到外部损伤的页岩力学参数。力学参数主要指弹性模量、泊松比、内聚力、抗压强度-围压的相关性。其中,抗压强度-围压的相关性通过开展不同围压下的三轴力学实验,获得不同围压下的抗压强度,进行线性拟合,从而确定抗压强度-围压的相关性方程。
C2、页岩吸水实验,具体如下:
S1,对页岩岩心开展自吸实验,获取不同自吸时间下的页岩岩心。其中,自吸实验采用端部自吸方式,水沿岩心一个端面进入页岩内部,吸水方向平行页岩岩心轴向方向。
S2,对不同自吸时间下的页岩岩心开展压汞实验,获取自吸过程中页岩动态孔隙结构参数。压汞实验是获取岩石孔隙结构的方法,其中,孔隙结构参数指页岩平均孔径λa。
由于水化作用,页岩自吸过程中,孔隙结构发生变化。基于此,本申请中对不同自吸时间下的页岩岩样开展压汞实验,从而获取自吸过程中的页岩平均孔径λa(t)。
S3、基于自吸过程中页岩动态孔隙结构特征,采用以上页岩动态自吸方程,计算得到不同自吸时间t下的页岩平均吸水长度Ls。
步骤D、通过井周应力分布方程求取井壁主应力。
基于地层三向地应力(垂向地应力、水平最大地应力、水平最小地应力),采用井眼坐标转换方程,得到井周应力分布:
式中:σr、σθ、σz为柱坐标系下径向、周向和轴向正应力;τθz、τrθ、τrz为柱坐标下θz、rθ、rz平面切应力;σxx、σyy、σzz为直角坐标下X、Y、Z轴方向上地应力分量;σxy、σxz、σyz为直角坐标系下XY、XZ、YZ平面上地应力分量;θc井周角;pi为钻井液液柱压力。
基于井周应力分布,可以得到井壁任意位置(井周角θc从0°到360°)的主应力分布,如下所示:
式中,σ1为最大主应力;σ3为最小主应力。
其中,钻井过程中,三向地应力(垂向地应力、水平最大地应力、水平最小地应力)是已知参数,地应力分量σxx、σyy、σzz、σxy、σxz、σyz可以通过三向地应力求取。因此,井壁任意位置(井周角从0°到360°)上的最大和最小主应力主要与钻井液柱压力pi相关。
E、基于井壁岩石主应力,采用页岩自吸-本构模型,实现钻井过程中的页岩地层坍塌压力预测,具体如下:
在不同钻井时间下,首先计算井壁最大主应力σ1和最小主应力σ3。在此基础上,借助页岩自吸-本构模型,计算得到当前最小主应力(围压)条件下的页岩应力—应变曲线(σi—εi)。基于应力—应变曲线(σi—εi),确定当前最小主应力σ3下的页岩应力σi的峰值σmax。当σmax等于井壁岩石最大主应力σ1时,井壁发生坍塌破坏。此时,井筒的液柱压力为地层坍塌压力。
基于上述基于自吸-本构模型的页岩地层坍塌压力预测方法,本申请公开一应用例。
应用例1
本应用例为中国某油田黑色页岩层位实际钻井情况,以直井为例,地质力学特征为:深度2650m,垂向地应力2.32g/cm3,水平最大地应力2.20g/cm3,水平最小地应力1.82g/cm3,孔隙压力1.10g/cm3。表面张力为86.5mN·m-1,接触角为40.1°,流体粘度默认为1mPa·s,流体密度为1.03g/cm3。
基于如图1所示的方法流程图,本应用例的步骤如下:
步骤1、制取页岩标准圆柱岩心。
表1原状页岩力学参数
力学参数 | 数值 |
弹性模量 | 32GPa |
泊松比 | 0.23 |
内聚力 | 15.4MPa |
抗压强度-围压相关性 | σ<sub>c</sub>=1.2σ<sub>3</sub>+72 |
步骤3、基于页岩三轴力学实验结果,查阅系数分布版图,可确定Weibull强度统计系数F的数值为2.3,Weibull强度统计系数m的数值为3.2。
步骤4、对页岩开展自吸实验,获取不同自吸时间下的页岩岩心。本应用例自吸实验采用端部自吸方式,水沿岩心一个端面进入页岩内部,吸水方向平行页岩岩心轴向方向,如图4所示。
步骤5、对不同自吸时间下的页岩岩心开展压汞实验,获取自吸过程中页岩动态孔隙结构参数。其中,孔隙结构参数主要指页岩平均孔径λa。
由于水化作用,页岩自吸过程中,孔隙结构发生变化。基于此,本应用例对自吸时间下的页岩岩样开展压汞实验,从而可以得到自吸过程中的页岩平均孔径,如下所示:
步骤6、基于自吸过程中页岩动态孔隙结构特征,采用页岩动态自吸方程,计算页岩内部平均吸水长度。
因页岩动态自吸方程为:
根据表面张力σ、接触角θ、水黏度μw、水密度p、自吸过程中的页岩平均孔径λa(t),采用上述的页岩动态自吸方程,可以得到不同自吸时间t下的页岩平均吸水长度Ls,如图5所示。
步骤7、根据步骤6获得的页岩平均吸水长度Ls,采用以下公式表达式如下:可求得损伤系数D:
上式中,L为页岩岩样长度,Ls为吸水长度。
步骤8、通过井周应力分布方程求取主应力。
基于地层三向地应力(垂向地应力、水平最大地应力、水平最小地应力),采用井眼坐标转换方程,得到井周应力分布:
式中:σr、σθ、σz为柱坐标系下径向、周向和轴向正应力;τθz、τrθ、τrz为柱坐标下θz、rθ、rz平面切应力;σxx、σyy、σzz为直角坐标下X、Y、Z轴方向上地应力分量;σxy、σxz、σyz为直角坐标系下XY、XZ、YZ平面上地应力分量;θc井周角;pi为钻井液液柱压力。
基于井周应力分布,可以得到井壁任意位置(井周角从0°到360°)的主应力分布,如下所示:
式中,σ1为最大主应力;σ3为最小主应力。
其中,钻井过程中,三向地应力(垂向地应力、水平最大地应力、水平最小地应力)是已知参数,地应力分量σxx、σyy、σzz、σxy、σxz、σyz可以通过三向地应力求取。因此,井壁任意位置(井周角从0°到360°)上的最大和最小主应力主要与钻井液柱压力pi相关。
步骤9、基于井壁岩石主应力,采用页岩自吸-本构模型,实现钻井过程中的页岩地层坍塌压力预测,具体如下:
在不同钻井时间下,首先计算井壁最大主应力σ1和最小主应力σ3。在此基础上,结合步骤2获得的力学参数、步骤3获得的Weibull强度统计系数、步骤7获得的损伤系数,采用页岩自吸-本构模型,计算得到当前最小主应力σ3条件下的页岩应力—应变曲线(σi—εi),。
基于应力—应变曲线(σi—εi),确定当前最小主应力σ3下的页岩应力σi的峰值σmax。当σmax等于井壁岩石最大主应力σ1时,井壁发生坍塌破坏。此时,井筒的液柱压力为地层坍塌压力。
考虑不同钻井时间下,井壁页岩水化损伤不同,从而计算得到坍塌压力随钻井时间的变化,如图6a所示。
为进一步说明计算过程,以钻井时间24小时为例,如图6b所示,基于页岩自吸-本构模型与井周主应力分布,计算得到:当钻井液柱压力pi=34.28MPa时(即是说最小主应力=34.28MPa),最大主应力为98.91MPa,与应力-应变曲线上的应力峰值(σmax=98.91MPa)相等。由此说明,此时井壁岩石发生坍塌破坏,地层坍塌压力为34.28MPa,转换为坍塌压力当量密度为1.32g/cm3。其中,当量密度转换方法为:
上式中,Pcd为坍塌压力当量密度;Pc为坍塌压力;g为重力加速度,常数;H为深度。
根据工区钻井资料,实际钻井过程中,钻井初期(24h以内)采用1.3g/cm3钻井液密度,稳定井壁良好。在钻井后期(24h以后),1.3g/cm3钻井液密度无法稳定地层,增加至1.36g/cm3钻井液密度,从而实现井壁稳定。实际钻井情况与预测坍塌压力变化趋势相吻合,从而证明了预测得到的坍塌压力具有适用性。
本发明的坍塌压力预测方法耦合了页岩吸水特征与本构关系,对指导现场钻井液密度优选具有重要意义。
以上所述的具体实施方式,对本申请的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种基于自吸-本构模型的页岩地层坍塌压力预测方法,其特征在于:包括以下步骤:
基于Weibull强度分布理论和损伤力学定义,构建损伤下的页岩本构方程;
基于损伤力学理论,结合页岩动态自吸方程、损伤下的页岩本构方程,构建水化损伤下的页岩自吸-本构模型;
制取页岩岩心,对页岩岩心开展三轴力学实验获得原状页岩力学参数;
对页岩岩心开展自吸实验,基于页岩动态自吸方程,计算出页岩平均吸水长度;
获取井壁岩石主应力分布;
基于井壁岩石主应力分布、原状页岩力学参数、页岩平均吸水长度,采用所述页岩自吸-本构模型,实现钻井过程中的页岩地层坍塌压力预测;具体的,根据井壁岩石最小主应力、原状页岩力学参数、页岩平均吸水长度,采用页岩自吸-本构模型,计算得到当前最小主应力条件下的页岩应力—应变曲线;基于应力—应变曲线,确定当前最小主应力下的页岩应力σi的峰值σmax;当σmax等于井壁岩石最大主应力时,井壁发生坍塌破坏;此时,井筒的液柱压力为地层坍塌压力;
所述损伤下的页岩本构方程为:
所述页岩动态自吸方程为:
上式中,Ls为页岩平均吸水长度,σ为表面张力,θ为接触角,μw为水黏度,p为水密度,λa(t)为自吸过程中的页岩平均孔径。
2.根据权利要求1所述的一种基于自吸-本构模型的页岩地层坍塌压力预测方法,其特征在于:所述原状页岩力学参数包括页岩弹性模量、泊松比、内聚力、抗压强度-围压的相关性;
抗压强度-围压的相关性通过以下方式获得:通过开展不同围压下的三轴力学实验,获得不同围压下的抗压强度,从而进行线性拟合,确定抗压强度-围压的相关性方程。
4.根据权利要求1所述的一种基于自吸-本构模型的页岩地层坍塌压力预测方法,其特征在于:所述的自吸实验具体包括:
S1,对页岩岩心开展自吸实验,获取不同自吸时间下的页岩岩心;
S2,对不同自吸时间下的页岩岩心开展压汞实验,获取自吸过程中的页岩平均孔径。
5.根据权利要求1或4所述的一种基于自吸-本构模型的页岩地层坍塌压力预测方法,其特征在于:所述的自吸实验采用端部自吸方式,水沿岩心一个端面进入页岩内部,吸水方向平行页岩岩心轴向方向。
7.根据权利要求1所述的一种基于自吸-本构模型的页岩地层坍塌压力预测方法,其特征在于:井壁岩石主应力分布的获取方法具体包括:
基于地层三向地应力,采用井眼坐标转换方程,得到井周应力分布:
上式中,σr、σθ、σz为柱坐标系下径向、周向和轴向正应力;τθz、τrθ、τrz分别为柱坐标下θz、rθ、rz平面切应力;σxx、σyy、σzz分别为直角坐标下X、Y、Z轴方向上地应力分量;σxy、σxz、σyz分别为直角坐标系下XY、XZ、YZ平面上地应力分量;θc井周角;pi为钻井液液柱压力;
基于井周应力分布,得到井壁任意位置的主应力分布,如下所示:
上式中,σ1为井壁岩石最大主应力;σ3井壁岩石最小主应力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110668186.8A CN113958315B (zh) | 2021-06-16 | 2021-06-16 | 一种基于自吸-本构模型的页岩地层坍塌压力预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110668186.8A CN113958315B (zh) | 2021-06-16 | 2021-06-16 | 一种基于自吸-本构模型的页岩地层坍塌压力预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113958315A CN113958315A (zh) | 2022-01-21 |
CN113958315B true CN113958315B (zh) | 2022-05-17 |
Family
ID=79460227
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110668186.8A Active CN113958315B (zh) | 2021-06-16 | 2021-06-16 | 一种基于自吸-本构模型的页岩地层坍塌压力预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113958315B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116066092B (zh) * | 2023-02-03 | 2024-04-26 | 西南石油大学 | 一种基于硬度实验的页岩地层坍塌压力增量预测方法 |
CN115961950B (zh) * | 2023-02-03 | 2024-06-25 | 西南石油大学 | 基于岩石声学信息的水岩作用下钻井液密度下限预测方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104483334A (zh) * | 2014-12-22 | 2015-04-01 | 中国矿业大学 | 一种可视化表征水泥基材料水分空间分布和含量的方法 |
CN108071389A (zh) * | 2017-12-10 | 2018-05-25 | 西南石油大学 | 一种页岩强度准则及页岩气井井眼坍塌压力预测方法 |
CN108593771A (zh) * | 2018-04-27 | 2018-09-28 | 西南石油大学 | 损伤强度计算方法和损伤强度计算装置 |
CN109187228A (zh) * | 2018-09-30 | 2019-01-11 | 西南石油大学 | 一种页岩地层钻井液稳定井壁能力的室内评价方法 |
CN109632510A (zh) * | 2019-01-15 | 2019-04-16 | 中石化重庆涪陵页岩气勘探开发有限公司 | 一种预测水化损伤页岩强度的方法 |
CN110208487A (zh) * | 2019-05-24 | 2019-09-06 | 西南石油大学 | 一种基于ct扫描的页岩水化损伤测试方法 |
WO2019237050A1 (en) * | 2018-06-07 | 2019-12-12 | Texas Tech University System | Method to generate microfractures by chemical reaction in low carbonate mineral content shale reservoirs |
CN111927446A (zh) * | 2020-09-11 | 2020-11-13 | 西南石油大学 | 一种水化页岩地层井壁坍塌失稳预测方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105003247A (zh) * | 2014-04-23 | 2015-10-28 | 中国海洋石油总公司 | 低密度钻井可行性快速评价方法 |
CN105334107B (zh) * | 2015-10-28 | 2017-09-29 | 西安石油大学 | 基于地层粘弹性的泥页岩井壁稳定钻井液密度确定方法 |
CN105863626A (zh) * | 2016-05-10 | 2016-08-17 | 西南石油大学 | 一种钻井液与泥页岩地层理化作用的评价方法 |
CN109555515A (zh) * | 2018-10-08 | 2019-04-02 | 中国石油大学(北京) | 地层坍塌压力确定方法和装置 |
CN111060420B (zh) * | 2019-12-31 | 2020-11-03 | 成都理工大学 | 一种描述页岩裂缝-孔隙流体自吸的方法 |
-
2021
- 2021-06-16 CN CN202110668186.8A patent/CN113958315B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104483334A (zh) * | 2014-12-22 | 2015-04-01 | 中国矿业大学 | 一种可视化表征水泥基材料水分空间分布和含量的方法 |
CN108071389A (zh) * | 2017-12-10 | 2018-05-25 | 西南石油大学 | 一种页岩强度准则及页岩气井井眼坍塌压力预测方法 |
CN108593771A (zh) * | 2018-04-27 | 2018-09-28 | 西南石油大学 | 损伤强度计算方法和损伤强度计算装置 |
WO2019237050A1 (en) * | 2018-06-07 | 2019-12-12 | Texas Tech University System | Method to generate microfractures by chemical reaction in low carbonate mineral content shale reservoirs |
CN109187228A (zh) * | 2018-09-30 | 2019-01-11 | 西南石油大学 | 一种页岩地层钻井液稳定井壁能力的室内评价方法 |
CN109632510A (zh) * | 2019-01-15 | 2019-04-16 | 中石化重庆涪陵页岩气勘探开发有限公司 | 一种预测水化损伤页岩强度的方法 |
CN110208487A (zh) * | 2019-05-24 | 2019-09-06 | 西南石油大学 | 一种基于ct扫描的页岩水化损伤测试方法 |
CN111927446A (zh) * | 2020-09-11 | 2020-11-13 | 西南石油大学 | 一种水化页岩地层井壁坍塌失稳预测方法 |
Non-Patent Citations (8)
Title |
---|
"Experimental study on crack propagation in shale formations considering hydration and wettability";liang lixi;《OURNAL OF NATURAL GAS SCIENCE AND ENGINEERING 》;20150331;全文 * |
"Experimental study on hydration damage mechanism of shale from the Longmaxi Formation in southern Sichuan Basin";Liu, X;《Petroleum 》;20160331;全文 * |
"Hydration index and hydrated constitutive model of clay shale using acoustic frequency spectrum";Liu, XJ;《ENERGY SCIENCE & ENGINEERING》;20190831;全文 * |
"多孔介质自发渗吸研究进展";蔡建超,等;《力学进展》;20121125;第42卷(第6期);全文 * |
"水相自吸诱发页岩裂缝起裂扩展行为研究";杨斌;《中国优秀博士论文全文库工程科技I辑》;20210115(第1期);全文 * |
"钻井液侵入泥页岩地层规律及其对井壁稳定性影响研究";邓富元;《中国优秀硕士论文全文库工程科技I辑》;20200615(第6期);全文 * |
"页岩地层钻完井岩石动力学形变研究";李泽;《中国优秀硕士论文全文库工程科技I辑》;20200615(第6期);全文 * |
李瑶."基于核磁共振和驱替试验联测的岩心自吸行为研究".《 第32届全国天然气学术年会(2020)论文集》.2020,全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN113958315A (zh) | 2022-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108868748B (zh) | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 | |
US11326449B2 (en) | Method for determining three-dimensional in-situ stress based on displacement measurement of borehole wall | |
CN113958315B (zh) | 一种基于自吸-本构模型的页岩地层坍塌压力预测方法 | |
US10928282B2 (en) | Method and device for determining elasticity of cement stone utilized in well cementing of oil-gas well | |
CN104778303B (zh) | 一种裂缝型地层防塌钻井液性能参数的设计方法 | |
CN107038290A (zh) | 一种考虑增压效应力的井壁坍塌压力的计算方法 | |
CN111927446B (zh) | 一种水化页岩地层井壁坍塌失稳预测方法 | |
CN113625364B (zh) | 一种基于双重校正的泥页岩地层孔隙压力的计算方法 | |
CN102128025B (zh) | 一种盐膏层造斜钻井液密度确定方法 | |
CN106644736A (zh) | 一种测量储层岩石的脆性指数的方法及装置 | |
CN105334107A (zh) | 基于地层粘弹性的泥页岩井壁稳定钻井液密度确定方法 | |
CN110984970A (zh) | 一种利用地层测试确定启动压力梯度的方法 | |
CN102661894A (zh) | 裂隙岩体试件与裂隙岩体水岩耦合试验方法 | |
CN114594530A (zh) | 一种计算地层中高角度天然裂缝内摩擦系数的方法 | |
CN107290799B (zh) | 一种岩石可压性的确定方法 | |
CN109555515A (zh) | 地层坍塌压力确定方法和装置 | |
CN112412434B (zh) | 一种改进的疏松砂岩地应力计算方法 | |
CN112304838B (zh) | 基于岩心三轴压缩的油页岩储层可改造性室内评价方法 | |
CN106772674B (zh) | 一种裂缝性砂岩储层含油饱和度的计算方法 | |
CN112836944A (zh) | 一种深井、超深井井壁稳定评价模型的建立方法 | |
CN104712299B (zh) | 适合气井控水增气压裂的设计方法 | |
CN113449353B (zh) | 考虑地层流体渗流作用的水平井井壁稳定分析方法 | |
CN113609702B (zh) | 一种页岩地层井壁天然裂缝扩展压力计算方法 | |
CN114776286A (zh) | 井壁稳定评价方法、装置、设备和钻井液处理剂优选方法 | |
CN115270411A (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 |