泥页岩井壁流-固-化耦合安全密度窗口计算方法
张世锋1, 汪海阁2, 邱正松3, 曹文科1, 黄洪春2, 陈志学2
1. 常州大学,江苏常州 213164
2. 中国石油集团工程技术研究院有限公司,北京 102206
3. 中国石油大学(华东),山东青岛 266580

第一作者简介:张世锋(1987-),男,山东济南人,博士,常州大学石油工程学院讲师,主要从事井壁稳定性、钻井液及岩石力学方面的研究工作。地址:江苏省常州市武进区常州大学(武进校区)石油工程学院,邮政编码:213164。E-mail: shifeng@cczu.edu.cn

摘要

针对页岩气开发钻井工程中亟待解决的泥页岩钻井井壁安全密度窗口计算难题,通过对Heidug-Wong的广义Biot有效应力原理进行修正,使用Weibull统计模型描述水化应变相关的强度损伤,考虑钻井液封堵层,结合流动-扩散耦合模型建立了钻井液抑制-封堵-渗透作用下的泥页岩安全密度窗口计算方法。以此为基础,分析了钻井液封堵性、抑制性等关键参数对安全密度窗口的影响:提高钻井液封堵性,能够阻缓泥页岩井壁孔隙压力传递,减少溶质扩散;钻井液的抑制性,尤其是对强度损伤的抑制,对水化显著的泥页岩井壁坍塌压力具有决定性的作用,钻井液封堵-抑制性能的提高可以降低坍塌压力、提高破裂压力,从而拓宽安全密度窗口,延长井壁的坍塌周期。若泥页岩内存在渗透流动,钻井液渗透作用能够拓宽安全密度窗口,且钻井液封堵性越强,其影响幅度越大。利用安全密度窗口计算方法可以有效分析坍塌压力、破裂压力与钻井液防塌性能的关系,可用于安全密度窗口的准确计算及钻井液性能的优化。图11表3参29

关键词: 泥页岩; 钻井液; 坍塌压力; 破裂压力; 井壁稳定; 安全密度窗口; 页岩气开发
中图分类号:TE345 文献标志码:A 文章编号:1000-0747(2019)06-1197-09
Calculation of safe drilling mud density window for shale formation by considering chemo-poro-mechanical coupling effect
ZHANG Shifeng1, WANG Haige2, QIU Zhengsong3, CAO Wenke1, HUANG Hongchun2, CHEN Zhixue2
1. Changzhou University, Changzhou 213164, China
2. CNPC Engineering Technology Research Institute Co., Ltd., Beijing 102206, China
3. China Petroleum University (East China), Qingdao 266580, China
Abstract

It is difficult to define safe drilling mud density window for shale sections. To solve this problem, the general Biot effective stress principle developed by Heidug and Wong was modified. The Weibull statistical model was used to characterize the hydration strain-related strength damage. Considering drilling fluid sealing barrier on shale, a calculation method of safe drilling mud density for shale formation under drilling fluid sealing-inhibition-penetration effect has been established combined with a flow-diffusion coupling model. The effects of drilling fluid sealing and inhibiting parameters on safe drilling fluid window were analyzed. The study shows that enhancing drilling fluid sealing performance can slow the pore pressure transmission and reduce solute diffusion; the inhibiting property of drilling fluid, especially inhibition to strength damage, is crucial for the wellbore collapse pressure of shale section with strong hydration property. The improvement of drilling fluid sealing and inhibition performance can lower collapse pressure and enhance fracturing pressure, and thus making the safe drilling fluid density window wider and the collapse cycle of wellbore longer. If there is osmosis flow in shale, induced osmosis flow can make the gap between collapse pressure and fracturing pressure wider, and the stronger the sealing ability of drilling fluid, the wider the gap will be. The safe drilling fluid density calculation method can analyze the relationships between collapse pressure, fracturing pressure and drilling fluid anti collapse performance, and can be used to optimize drilling fluid performance.

Keyword: shale; drilling fluid; collapse pressure; fracturing pressure; wellbore stability; safe mud density window; shale gas development
0 引言

随着页岩油气勘探与开发的发展, 泥页岩层段突出的井壁稳定问题引发了越来越多的关注[1, 2]。“ 多元协同” 防塌钻井液技术在泥页岩钻井中得到了较为广泛的应用, 即通过“ 封堵-抑制-有限活度平衡” 方法增加井壁稳定性[3], 但如何准确预测钻井液作用下泥页岩的水化应力及水化损伤, 进而计算泥页岩段的安全密度窗口仍是国内外学者持续研究的问题, 也是目前泥页岩段钻井亟待解决的关键技术。与采用从微观角度统计分析然后拓展到宏观尺度估算泥页岩水化变形的方法相比[4, 5, 6], 直接进行宏观本构模型的研究更为简单[7, 8, 9]。Heidug、Ghassemi与Zhou等提出并应用了化学线弹性的泥页岩本构方程, 认为泥页岩有效应力是由泥页岩变形、孔隙压力及孔隙流体化学势变化所控制的[7, 8, 9]; Zhang等[10, 11]使用CT扫描观察了水化裂缝的扩展规律, 发现水化应力的变化过程与Heidug等[7]提出的考虑水化学应力的广义Biot有效应力方程所阐述的规律有所不同。另外, 泥页岩中黏土矿物水化引起水化裂缝的扩展并伴随着泥页岩的强度损伤[10, 11, 12, 13, 14], 现有的主要流-固-化耦合模型并未考虑。

本文基于水化应力的变化规律, 对Heidug-Wong广义Biot有效应力原理[7]进行修正, 使用Weibull统计模型描述水化应变相关的强度损伤, 同时考虑钻井液封堵层, 基于流动-扩散耦合模型建立井壁孔隙压力、溶质浓度的求解方法, 进而建立钻井液作用下泥页岩安全密度窗口计算方法, 探索准确计算泥页岩井段坍塌压力、破裂压力及其随时间变化规律的方法, 为钻井过程中的钻井液防塌性能设计提供依据。

1 有效应力原理修正与损伤统计模型
1.1 泥页岩水化应力变化特性及Heidug-Wong有效应力原理的修正

1.1.1 细观水化裂缝宽度表征及水化应力演化特性

如图1所示, 使用Griffith裂缝表示水化张应力及围压作用下的泥页岩黏土矿物水化膨胀裂缝[11]

图1 承受净张应力的Griffith水化裂缝模型[9]

裂缝内水化应力(${{\sigma }_{\text{hf}}}$)为张应力, 克服裂缝闭合压力(${{\sigma }_{\text{c}}}-\alpha p$)引起水化裂缝的扩展。根据Griffith准则, 水化过程中的裂缝缝宽分布模型可以表述为:

$\Delta y\left( x \right)=\frac{2\left( 1-\nu \right)}{{{G}_{\text{s}}}}\left[ {{\sigma }_{\text{hf}}}-\left( {{\sigma }_{\text{c}}}-\alpha p \right) \right]\sqrt{{{L}_{\text{h}}}^{2}-{{x}^{2}}}$ (1)

忽略孔隙压力的影响, 则

${{\sigma }_{\text{hf}}}-{{\sigma }_{\text{c}}}=\frac{{{G}_{\text{s}}}}{2(1-\nu )\sqrt{L{{_{\text{h}}^{{}}}^{2}}-{{x}^{2}}}}\Delta y\left( x \right)$ (2)

Zhang等[10, 11]观察了0.03 MPa孔隙压力、0.10 MPa和20.00 MPa围压条件下, 页岩吸水过程中裂缝的扩展, 共计存在2种裂缝扩展趋势:①$\Delta y\left( x \right)$持续降低, 则${{\sigma }_{\text{hf}}}$< ${{\sigma }_{\text{c}}}$且${{\sigma }_{\text{hf}}}$持续减小; ②$\Delta y\left( x \right)$先增加后降低, 则初始阶段${{\sigma }_{\text{hf}}}$> ${{\sigma }_{\text{c}}}$, 随着时间的延长, ${{\sigma }_{\text{hf}}}$< ${{\sigma }_{\text{c}}}$且${{\sigma }_{\text{hf}}}$持续减小, 由于围压${{\sigma }_{\text{c}}}$保持恒定, 可以推断${{\sigma }_{\text{hf}}}$初始阶段较大, 随后随着时间的延长逐渐降低。Heidug等[7]认为黏土矿物层间化学势差的变化是黏土矿物膨胀的驱动力。随着溶质的运移、水分子的扩散, 孔隙流体与黏土矿物层间的化学势差降低, 则黏土矿物层间化学势差的变化逐渐减小, 造成水化膨胀压降低。因此, 可以推断:

${{\sigma }_{\text{hf}}}\propto \frac{\text{d}\left[ \sum\limits_{k}{\mu \left( {{c}_{k}} \right)} \right]}{\text{d}t}$ (3)

1.1.2 泥页岩溶质扩散及残余水化膨胀压

泥页岩中含有大量的黏土矿物(见图2a), 水和溶质可以在多尺度的孔隙中进行传递(见图2b), 这些孔隙包括纳米尺度黏土矿物层间孔及大孔隙[15]。大孔中的渗透流动与膜效率相关, 当膜效率为零时, 则不存在渗透流动。对于黏土矿物层间孔, 化学势差驱动下的溶质或水分子运移在层间产生水化膨胀斥力, 引起黏土矿物层间距的增加。随着化学势差的降低, 层间的水化膨胀斥力逐渐降低, 直至与外部应力平衡(见图2c), 此时溶质运移被阻止, 化学势差保持恒定, 层间存在残余的水化膨胀斥力, 即为残余膨胀压[15, 16]

图2 泥页岩及其中溶质扩散、水化应力分布示意图

使用Fick定律描述泥页岩中溶质的扩散[17, 18]

${{J}_{\text{da}}}=-\phi {{D}_{\text{e}}}\frac{\text{d}c}{\text{d}x}$ (4)

考虑外部应力对黏土矿物层间孔内水、溶质运移的阻碍作用时, ${{D}_{\text{e}}}$可以表示为[17, 18, 19]

$\left\{ \begin{align} & {{D}_{\text{e}}}={{\alpha }_{\text{ }\!\!\pi\!\!\text{ }}}{{D}_{\text{a}}} \\ & {{\alpha }_{\text{ }\!\!\pi\!\!\text{ }}}\text{=}1-\frac{{{\omega }^{2}}{{C}_{\text{ac}}}{{K}_{\text{h}}}RT}{{{D}_{\text{a}}}\phi {{\rho }_{\text{w}}}g} \\ & {{D}_{\text{a}}}=\phi \tau {{D}_{\text{s}}} \\ \end{align} \right.$ (5)

1.1.3 泥页岩广义Biot有效应力原理的修正

Heidug-Wong建立的与化学势相关的泥页岩有效应力随时间的变化为[7]

$\frac{\partial {{\sigma }_{\text{e, }i\text{, }j}}}{\partial t}=\frac{\partial {{\sigma }_{i\text{, }j}}}{\partial t}-\left( \alpha \frac{\partial p}{\partial t}-\sum\limits_{k}{{{\omega }_{k}}\frac{\partial {{\mu }_{k}}}{\partial t}} \right){{\delta }_{i\text{, }j}}$ (6)

对(6)式积分可得:

${{\sigma }_{\text{e, }i\text{, }j}}={{\sigma }_{i\text{, }j}}-\left( \alpha p-\int\limits_{0}^{t}{\sum\limits_{k}{{{\omega }_{k}}\frac{\partial {{\mu }_{k}}}{\partial {{c}_{k}}}\frac{\partial {{c}_{k}}}{\partial t}\text{d}t}} \right){{\delta }_{i, j}}$ (7)

假设孔隙流体中仅有一种与膨胀相关的溶质, 根据水化应力与化学势差正相关的变化规律, 考虑残余水化膨胀压, 使用下述方程描述水化应力为:

$\int\limits_{0}^{t}{\sum\limits_{k}{{{\omega }_{k}}\frac{\partial {{\mu }_{k}}}{\partial {{c}_{k}}}\frac{\partial {{c}_{k}}}{\partial t}\text{d}t}}={{\omega }_{0}}\frac{\partial c}{\partial t}+{{\pi }_{\text{d}}}$ (8)

Zhang等[10, 11]的观察结果中, 泥页岩仅有少量水化裂缝在残余水化膨胀压的作用下保持张开, 因此黏土矿物中残余膨胀压对泥页岩有效应力的影响可使用系数${{\alpha }_{\text{d}}}$来校正, 修正后的有效应力模型可表示为:

${{\sigma }_{\text{e, }i\text{, }j}}={{\sigma }_{i\text{, }j}}-\left[ \alpha p-\left( {{\omega }_{\text{0}}}\frac{\partial c}{\partial t}+{{\alpha }_{\text{d}}}{{\pi }_{\text{d}}} \right) \right]{{\delta }_{i\text{, }j}}$ (9)

1.2 含水化应变损伤的泥页岩Weibull统计损伤模型

泥页岩水化过程中体积膨胀引起水化裂缝扩展、强度损伤[10, 11, 12, 13, 14], 国内外学者建立了不同的理论模型来描述岩石细观损伤与宏观破坏的关系[20, 21, 22, 23, 24, 25], 其中Weibull统计损伤模型是广泛采用的模型之一[23, 24, 25]。Zhang等[11]的实验结果表明, 随着泥页岩水化膨胀应变由7.160× 10-4增加至1.662× 10-3、3.170× 10-3, 水化前强度为62.2 MPa的Mancos泥页岩, 水化后强度由43.1 MPa降至27.9 MPa、0 MPa(水化破碎)。因此, 选用水化膨胀应变作为强度基元参数, 泥页岩水化膨胀Weibull统计损伤可以表示为:

$D\text{=}1-\text{exp}\left[ -{{\left( \frac{{{\varepsilon }_{\text{v, chem}}}}{{{\varepsilon }_{\text{v0, chem}}}} \right)}^{m}} \right]$ (10)

地层应力作用下, 钻井液中泥页岩破坏可以看作水化膨胀应力与剪切应力的两级加载, 其总损伤效应可以表示为[25]

${{D}_{\text{m}}}\text{=}{{D}_{\text{ms}}}+D-{{D}_{\text{ms}}}D$ (11)

假设未水化泥页岩强度损伤为零, 则

${{D}_{\text{m}}}\text{=}D$ (12)

根据Heidug-Wong假设[7], 泥页岩水化过程中存在局部的物理-化学平衡, 水化作用速率远远高于应变速率, 即可假设水化膨胀过程中压力及外部约束恒定, 则水化膨胀体积应变随时间的变化可以表示为:

$\frac{\partial {{\varepsilon }_{\text{v, chem}}}}{\partial t}=\frac{\alpha -1}{{{G}_{\text{K}}}}\sum\limits_{k}{{{\omega }_{k}}\frac{\partial {{\mu }_{k}}}{\partial {{c}_{k}}}\frac{\partial {{c}_{k}}}{\partial t}}$ (13)

假设孔隙流体中仅有一种与膨胀相关的溶质, 令${{\omega }_{1}}\text{=}\frac{\alpha -1}{{{G}_{\text{K}}}}{{\omega }_{\text{k0}}}\frac{\partial \mu }{\partial c}$, 水化应变、溶质浓度连续变化, 对(13)式积分, 并联立(10)式可得:

$\left\{ \begin{align} & D\text{=}1-{{\text{e}}^{-{{\left( \gamma c \right)}^{m}}}} \\ & \gamma \text{=}\frac{{{\omega }_{1}}}{{{\varepsilon }_{\text{v0, chem}}}} \\ \end{align} \right.\ \ \ $ (14)

从(14)式可以看出, 随着黏土矿物层间孔内水分子、溶质含量的增加, 损伤程度增加。

2 安全密度窗口计算模型
2.1 考虑钻井液封堵层的井筒模型

沥青、高分子聚合物等在泥页岩表面黏附, 另外, 聚合醇、硅酸盐、铝酸盐等可以在泥页岩孔隙内滤液前缘形成沉淀封堵孔隙[1]。因此, 钻井液封堵层由黏附层及封堵剂充填的表层泥页岩构成, 具有有限半透膜的性质, 且较空白泥页岩更加致密。

据此, 可建立如图3所示的井筒物理模型, 井壁内侧存在封堵层, 能够影响钻井液内水、溶质向井壁的运移, 进而影响井壁处的孔隙压力和溶质浓度。

图3 1/2井筒物理模型

2.2 封堵层、泥页岩渗流模型

泥页岩孔隙流体的质量守恒方程可以表示为[2]

$\frac{1-{{G}_{\text{K}}}/{{G}_{\text{Km}}}-\phi }{{{G}_{\text{Km}}}}\frac{\partial p}{\partial t}-\frac{K}{\eta }{{\nabla }^{2}}p=-{{10}^{-6}}\frac{{{\sigma }_{\text{m}}}K}{\eta }\frac{{{\rho }_{\text{f}}}RT}{{{M}_{\text{s}}}{{C}_{\text{sD}}}}{{\nabla }^{2}}c$ (15)

基于连续性方程, 忽略化学势、固体变形、溶质吸附对流体流动的影响, 封堵层孔隙流体质量守恒方程可以表示为:

$\frac{{{\phi }_{\text{s}}}}{{{G}_{\text{Kf}}}}\frac{\partial p}{\partial t}-\frac{{{K}_{\text{s}}}}{\eta }{{\nabla }^{2}}p=-{{10}^{-6}}\frac{{{\sigma }_{\text{ms}}}{{K}_{\text{s}}}}{\eta }\frac{{{\rho }_{\text{f}}}RT}{{{M}_{\text{s}}}{{C}_{\text{sD}}}}{{\nabla }^{2}}c$ (16)

假设封堵层渗透率与孔隙度间的关系与泥页岩一致, 则封堵层孔隙度[14]可表示为:

${{\phi }_{\text{s}}}=\phi \frac{\lg {{K}_{\text{s}}}-a}{\lg K-a}$ (17)

如图3所示, 渗流边界条件为:

$\left\{ \begin{align} & {{\left. p \right|}_{r={{r}_{\text{w}}}-{{d}_{\text{s}}}}}={{p}_{\text{df}}} \\ & {{\left. \frac{{{K}_{\text{s}}}}{\eta }\nabla {{p}_{\text{s}}} \right|}_{r={{r}_{\text{w}}}}}=-\frac{K}{\eta }\nabla p \\ & {{\left. p \right|}_{r=\infty }}={{p}_{\text{0}}} \\ & {{\left. p \right|}_{t=0}}={{p}_{\text{0}}} \\ \end{align} \right.$ (18)

2.3 泥页岩及封堵层溶质扩散模型

泥页岩溶质的质量守恒方程为:

$\phi \frac{\partial c}{\partial t}-{{D}_{\text{e}}}{{\nabla }^{2}}c=0$ (19)

封堵层内溶质的质量守恒方程为:

${{\phi }_{\text{s}}}\frac{\partial c}{\partial t}-{{D}_{\text{es}}}{{\nabla }^{2}}c=0$ (20)

根据(4)式和(5)式, 假设封堵层和泥页岩的扩散校正系数、几何特征参数一致, 则封堵层的扩散系数可表示为:

${{D}_{\text{es}}}={{D}_{\text{e}}}\frac{{{\phi }_{\text{s}}}}{\phi }$ (21)

无限远处, 由于泥页岩孔隙流体中K+或其他抑制剂含量极少, 可忽略。因此, 扩散边界条件为:

$\left\{ \begin{align} & {{\left. c \right|}_{r={{r}_{\text{w}}}-{{d}_{\text{s}}}}}={{c}_{\text{f}}} \\ & {{\left. {{D}_{\text{es}}}\nabla {{c}_{\text{s}}} \right|}_{r={{r}_{\text{w}}}}}=-{{D}_{\text{e}}}\nabla c \\ & {{\left. c \right|}_{r=\infty }}={{c}_{\text{0}}} \\ & {{\left. c \right|}_{t=0}}={{c}_{\text{0}}} \\ \end{align} \right.$ (22)

2.4 钻井液作用下考虑水化应力及强度损伤的泥页岩安全密度窗口计算

基于Lemaitre J[26]提出的应变等价性假说, 修正后考虑损伤的有效应力原理为:

$\sigma _{i\text{, }j}^{\text{* }}=\frac{{{\sigma }_{i\text{, }j}}-\left[ \alpha p-\left( {{\omega }_{\text{0}}}\frac{\partial c}{\partial t}+{{\alpha }_{\text{d}}}{{\pi }_{\text{d}}} \right) \right]{{\delta }_{i\text{, }j}}}{1-D}$ $\left( i, j=1, 2, 3 \right)$ (23)

库伦摩尔破坏准则是目前井壁稳定分析中广泛应用的剪切破坏强度准则之一[23], 可以表示为:

$\left( \sigma _{1}^{* }-\sigma _{3}^{* } \right)-\left( \sigma _{1}^{* }+\sigma _{3}^{* } \right)\sin \varphi =2C\cos \varphi $ (24)

拉伸破坏强度准则可以表示为:

$\sigma _{3}^{* }=-{{\sigma }_{\text{t}}}$ (25)

使用Comsol有限元软件求解(15)、(16)式可得井壁孔隙压力, 求解(19)、(20)式可得溶质浓度, 并根据(14)式求取与溶质浓度相关的损伤。使用MATLAB软件根据地应力求取井周应力分布、井壁主应力, 根据修正后考虑损伤的有效应力及剪切和拉伸破坏准则, 求解坍塌压力、破裂压力, 具体求解流程如图4所示。

图4 泥页岩安全密度窗口计算流程图

3 泥页岩井壁安全密度窗口计算

泥页岩防塌钻井液主要利用其封堵、抑制性能提高井壁稳定性, 若泥页岩存在渗透现象, 也可通过增加钻井液无机盐浓度以诱导孔隙流体流出, 降低孔隙压力, 提高井壁稳定性。采用建立的安全密度窗口计算方法, 根据文献[2, 15-18, 23-25, 27-29]中泥页岩封堵、扩散、膨胀等参数的测试结果, 计算安全密度窗口(基础参数见表1), 探讨钻井液不同防塌性能对泥页岩井壁孔隙压力、溶质浓度、水化应力、强度损伤及安全密度窗口的影响。

表1 泥页岩安全密度窗口计算基础参数表
3.1 钻井液封堵性能对井壁孔隙压力及安全密度窗口的影响

3.1.1 钻井液封堵性能对井壁泥页岩孔隙压力传递的影响

泥页岩及封堵层膜效率为零时, 由(15)、(16)式可知, 井壁孔隙压力的影响因素为钻井液封堵层厚度及渗透率。取表1中的基础参数, 计算不同封堵层渗透率、厚度条件下井壁处的孔隙压力(见图5):①井壁孔隙压力与时间正相关, 表现为初期快速上升, 1.5 d内曲线均出现拐点, 随后缓慢上升; ②封堵层厚度一定, 封堵层渗透率越低, 井壁孔隙压力越低; ③封堵层厚度对井壁孔隙压力的影响程度随封堵层渗透率的增大而大幅降低, 如当Ks=1× 10-10 μ m2时, 封堵层厚度的变化对井壁孔隙压力的影响变得很小, 10 d后3条线接近重合; ④封堵层渗透率一定, 封堵层厚度越小, 井壁孔隙压力越高。可以看出, 封堵层渗透率与厚度是控制孔隙压力传递的关键因素。

图5 渗透率及封堵层厚度对井壁处孔隙压力的影响

3.1.2 钻井液渗透作用对井壁泥页岩孔隙压力的影响

根据渗流-扩散耦合模型可知, 化学势差作用下, 若钻井液活度(无机盐浓度)保持恒定, 则泥页岩井壁处化学势差诱导渗透压力的影响因素包括封堵层厚度、渗透率、膜效率。设定p0为30 MPa, Ks为(0.01~1.00)× 10-10 μ m2, ds为1, 2 mm; 封堵层膜效率取0.3、0.5, 采用表1中相关参数计算p组合1— 12条件下(见表2)井壁处孔隙压力(见图6)。

表2 泥页岩井壁孔隙压力计算设定参数组合表

图6 不同参数组合情况下井壁处孔隙压力变化

分析可知, 在钻井液的渗透作用下:①井壁孔隙压力随时间逐渐降低至最小值, 然后逐渐增加并趋于平稳; ②封堵层厚度一定, 封堵层渗透率越低, 井壁孔隙压力降低幅度越大, 降至最低值后增加的速率越缓慢, 如当ds=1 mm, σ ms=0.3时, 渗透率由1.00× 10-10 μ m2降低为0.10× 10-10 μ m2, 0.01× 10-10 μ m2, 孔隙压力最低值由20.39 MPa降低为16.86 MPa、14.44 MPa, 10 d后孔隙压力由21.87 MPa降低为20.76 MPa、15.36 MPa。③封堵层厚度增加时, 井壁孔隙压力降低速率及降低至最低值后增加的速率均减缓, 孔隙压力最低值小幅降低, 如当Ks=0.10× 10-10 μ m2时, ds由1 mm(p组合5)增至2 mm(p组合7), 孔隙压力最低值由16.86 MPa降低为16.84 MPa, 10 d后孔隙压力由20.76 MPa降低为19.67 MPa。④封堵层膜效率的提高对井壁孔隙压力的影响程度随封堵层渗透率的减小而降低。如当Ks=1.00× 10-10 μ m2时, σ ms由0.3增至0.5, ds为1, 2 mm条件下, 孔隙压力最低值降幅分别为Δ pD1=0.03 MPa, Δ pD2=0.16 MPa; 当Ks=0.10× 10-10 μ m2, 0.01× 10-10 μ m2时, 膜效率的变化对井壁孔隙压力影响甚微, 曲线几乎重合。可以看出, 封堵层渗透率与厚度是钻井液渗透作用影响孔隙压力的关键因素。

3.2 钻井液封堵-抑制作用对井壁泥页岩有效应力及强度损伤的影响

根据与化学势相关的有效应力原理, 井壁泥页岩水化膨胀压${{\sigma }_{\text{he}}}$可以表示为:

$\left\{ \begin{align} & {{\sigma }_{\text{e, }i\text{, }j}}={{\sigma }_{i\text{, }j}}-{{\sigma }_{\text{he}}}{{\delta }_{i\text{, }j}} \\ & {{\sigma }_{\text{he}}}\text{=}{{\omega }_{\text{0}}}\frac{\partial c}{\partial t}+{{\alpha }_{\text{d}}}{{\pi }_{\text{d}}} \\ \end{align} \right.\ $ (26)

根据渗流-扩散耦合模型可知, 若孔隙压力保持恒定, 则井壁处泥页岩水化应力的影响因素包括溶质浓度、膨胀系数、残余膨胀压; 强度损伤D(14(式))的影响因素包括γ m。如m取值1.5, γ 取值0.005~0.030, 则D的取值为0~0.994 5。

设定溶质在泥页岩中的扩散系数为1× 10-10 m2/s, 在封堵层中扩散系数为(0.01~1.00)× 10-11 m2/s, 封堵层厚度为1~3 mm, 泥页岩中抑制剂浓度为零, 钻井液中抑制剂(KCl)浓度为100 g/L, 计算不同扩散系数下井壁处溶质浓度的变化规律(见图7)。分析表明:①井壁处溶质浓度与时间正相关, 表现为初期快速上升, 约1 d后出现拐点, 随后缓慢上升; ②扩散系数越低, 井壁处溶质浓度越低, 且随时间延长, 扩散系数对井壁处溶质浓度的影响程度初期差异较大, 后期趋于稳定; ③封堵层越厚, 井壁处溶质浓度越低, 同样随时间延长, 对井壁处溶质浓度的影响程度初期差异较大, 后期趋于稳定。可见, 降低扩散系数可以减少溶质扩散。根据(5)式, 溶质的有效扩散系数与外部应力约束相关, 随着外部应力的增加, 溶质的扩散减少。

图7 扩散系数对井壁处溶质浓度的影响

图8为扩散系数对泥页岩中水化膨胀压的影响(水化膨胀系数ω 0取2.0 MPa· s· m3/mol), 初始时刻, 未发生溶质/水的扩散, 虽然扩散系数不同, 但由于泥页岩性质及两侧溶质浓度相同, 即化学势差相同, 水化膨胀压相同; 随着时间的延长, 水化膨胀压逐渐降低, 并趋于恒定; 扩散系数的降低对水化膨胀压的影响较大, 时间为10 d时, Des=1× 10-12 m2/s与Des= 1× 10-13 m2/s的水化膨胀压差为0.31 MPa, 而Des= 1× 10-11 m2/s与Des=1× 10-12 m2/s的水化膨胀差为2.58 MPa, 约为前者的8倍。

图8 扩散系数对水化膨胀压的影响

图9为扩散系数、损伤系数对岩石强度损伤的影响计算结果, 可以看出:①岩石强度损伤与时间正相关, 表现为初期快速上升, 约1 d左右出现拐点, 随后缓慢上升; ②损伤系数越大, 岩石强度损伤越大, 且随时间延长对岩石强度损伤的影响程度逐步加大; ③扩散系数降低, 由于溶质运移/水扩散变缓, 引起的水化膨胀应变减小, 岩石的强度损伤降低。

图9 扩散系数、损伤系数对岩石强度损伤的影响

3.3 钻井液封堵-抑制-渗透作用对井壁安全密度窗口的影响

根据上述讨论, 钻井液封堵性及渗透性主要影响井壁孔隙压力传递及溶质扩散, 可采用封堵层渗透率、扩散系数、厚度进行表征。封堵层厚度可调整幅度较小, 且厚度较低的封堵层更加致密, 渗透率更低[2], 综合考虑(17)式和(21)式, 选取封堵层渗透率作为主要表征参数。钻井液抑制性主要影响井壁泥页岩的有效应力及强度, 可采用膨胀系数、损伤系数及残余水化膨胀压进行表征。因此, 钻井液封堵层渗透率、膨胀系数、损伤系数、残余水化膨胀压是影响井壁安全密度窗口的主要参数。

采用表1中的基础参数, 设定ds=0.001 m, σ ms= 1.5σ m, m=1.5, 同时假设钻井液抑制性提高, ω 0γ α dπ d同时降低, 且降低幅度一致, 根据钻井液抑制性的不同, 设计如表3所示的参数组合计算不同封堵性(Ks)、抑制性(ω 0, γ , α dπ d)及反渗透作用(σ m)条件下泥页岩井壁的安全密度窗口变化(见图10)。

表3 泥页岩井壁安全密度窗口计算设定参数组合表

图10 封堵层渗透率、膨胀系数、损伤系数对泥页岩井壁破裂压力和坍塌压力当量密度的影响

分析可知, ①井壁坍塌压力当量密度与时间正相关, 表现为初期快速上升, 约1 d后出现拐点, 随后缓慢上升, 井壁破裂压力当量密度与时间负相关, 表现为初期快速降低, 约1 d后出现拐点, 随后缓慢降低; ②钻井液抑制性提高, ω 0γ α dπ d降低, 坍塌压力当量密度大幅降低, 破裂压力当量密度略有提高, 如Ks= 1.0× 10-10 μ m2, 无渗透作用时, 随着参数组合3变为参数组合2和参数组合1, 10 d后安全密度窗口由无安全密度窗口(坍塌压力当量密度2.07 g/cm3, 破裂压力当量密度1.98 g/cm3)分别变为1.56~2.01 g/cm3, 1.10~2.03 g/cm3; ③封堵层渗透率降低, 坍塌压力当量密度降低、破裂压力当量密度提高, 安全密度窗口拓宽; 封堵层渗透率对坍塌压力当量密度的影响幅度随着坍塌压力当量密度的增加而提高, 如钻井液封堵层渗透率由1.0× 10-10 μ m2降至1.0× 10-11 μ m2时, 无反渗透作用情况下, 10 d后参数组合1的安全密度窗口1.10~2.03 g/cm3拓宽为参数组合7的1.10~2.08 g/cm3; 参数组合3由无安全密度窗口(坍塌压力当量密度2.07 g/cm3, 破裂压力当量密度1.98 g/cm3)拓宽为参数组合9的1.98~2.03 g/cm3; ④若泥页岩内存在渗透流动, 钻井液渗透作用能够降低坍塌压力当量密度, 提高破裂压力当量密度, 拓宽安全密度窗口, 且渗透率越低, 其影响程度越大, 如在钻井液的渗透作用下, 封堵层渗透率为1.0× 10-10 μ m2, 10 d后参数组合1的安全密度窗口由1.10~2.03 g/cm3拓宽为参数组合4的安全密度窗口1.09~2.04 g/cm3; 封堵层渗透率为0.1× 10-10 μ m2时, 10 d后参数组合7的安全密度窗口由1.10~2.08 g/cm3拓宽为参数组合10的安全密度窗口1.02~2.15 g/cm3

由此可知, 对于水化膨胀强度损伤显著的泥页岩, 钻井液抑制性的提高能够显著降低井壁坍塌压力当量密度; 钻井液封堵性提高, 封堵层渗透率降低, 能够阻缓压力的传递及溶质的扩散, 延长坍塌周期, 拓宽安全密度窗口, 且封堵层渗透率对坍塌压力当量密度的影响幅度随着坍塌压力当量密度的增加而提高; 若泥页岩内存在渗透流动, 钻井液的渗透作用可以降低孔隙压力, 提高有效应力, 从而实现坍塌压力当量密度的降低及破裂压力当量密度的增加, 拓宽安全密度窗口, 且渗透率越低, 其影响程度越大。

ω 0γ α dπ d是与钻井液抑制性相关的3个参数, 保持γ 不变, ω 0取值0, 1, 2 MPa· s· m3/mol; α dπ d取值0, 0.5, 1.0 MPa; Ks取值1.0× 10-11 μ m2, 若封堵层无渗透作用, 安全密度窗口计算结果见图11, 分析可知:①ω 0的提高引起坍塌压力当量密度的增加和破裂压力当量密度的降低; 随着时间的延长, ω 0的变化对坍塌压力和破裂压力影响幅度逐渐降低, 0.25 d后影响可以忽略, 如随着ω 0增加, 初始时刻安全密度窗口由参数组合13的1.09~2.90 g/cm3变化为参数组合14时的1.19~2.71 g/cm3和参数组合9时的1.28~2.53 g/cm3; 0.25 d时参数组合13、14、6的安全密度窗口均为1.98~2.03 g/cm3。②α dπ d的增加引起坍塌压力当量密度的增加和破裂压力当量密度的降低, 如10 d后参数组合15的安全密度窗口1.94~2.05 g/cm3变化为参数组合16的1.96~2.04 g/cm3和参数组合9的1.98~2.03 g/cm3。与图10中不同钻井液抑制性(ω 0, γ , α dπ d)对坍塌压力当量密度的影响幅度相比, ω 0α dπ d的变化对安全密度窗口的影响较小, 即钻井液抑制性提高引起的强度损伤(γ )降低是坍塌压力当量密度变化的决定性因素。

图11 封堵层渗透率、膨胀系数、损伤系数对泥页岩井壁当量钻井液密度的影响

4 结论

提高钻井液封堵性, 能够阻缓泥页岩井壁孔隙压力传递, 减少溶质扩散; 钻井液的抑制性, 尤其是对泥页岩强度损伤的抑制, 对水化显著的泥页岩井壁坍塌压力具有决定性作用。钻井液封堵-抑制性能提高, 能够降低坍塌压力、提高破裂压力, 从而拓宽安全密度窗口, 延长井壁的坍塌周期; 坍塌压力越高, 钻井液封堵性对其影响程度越大。

若泥页岩内存在渗透流动, 钻井液渗透作用能够拓宽安全密度窗口, 且钻井液封堵性越强, 其影响幅度越大。

通过修正广义Biot有效应力原理、使用Weibull统计模型描述水化应变相关的强度损伤建立的安全密度窗口计算方法可以有效分析坍塌压力、破裂压力与钻井液防塌性能的关系, 可用于安全密度窗口的准确计算及钻井液性能的优化。

符号注释:

a— — 常数, -24.6; c— — 溶质的摩尔浓度, mol/m3; c0— — 泥页岩内抑制剂初始浓度, mol/m3; cdf— — 钻井液中抑制剂浓度, mol/m3; cf— — 钻井液内溶质的摩尔浓度, mol/m3; cid— — 黏土矿物层间溶质含量, mol/m3; cip— — 流体进入孔隙的溶质浓度, mol/m3; cis— — 流动源孔隙溶质浓度, mol/m3; ck— — k组分溶质摩尔浓度, mol/m3; cs— — 封堵层内溶质的摩尔浓度, mol/m3; C— — 泥页岩内聚力, MPa; Cac— — 溶质离子平均浓度, mol/m3; CsD— — 溶质平均浓度, mol/m3; d— — 黏土矿物层间距, 0.1 nm; ds— — 封堵层的厚度, mm; D— — 水化膨胀损伤变量, 无因次; Da— — 溶质在孔隙中的自由扩散系数, m2/s; De— — 有效扩散系数, m2/s; Des— — 封堵层内溶质的扩散系数, m2/s; Dm— — 水化与荷载作用下岩石的总损伤变量, 无因次; Dms— — 荷载损伤变量, 无因次; Ds— — 溶质在水中的自由扩散系数, m2/s; g— — 重力加速度, 9.8 m/s2; GK— — 泥页岩的体积模量, MPa; GKf— — 液体体积模量, MPa; GKm— — 泥页岩基质体积模量, MPa; Gs— — 剪切模量, MPa; Jda— — 溶质扩散通量, mol/(m2· s); i, j— — 分量编号; k— — 溶质组分编号; K— — 渗透率, 103 μ m2; Kh— — 流动效率, m/s; Ks— — 封堵层的渗透率, 103 μ m2; L— — 距井壁的距离, m; Lh— — 裂缝的半长, m; m— — 损伤系数, 无因次; Ms— — 溶质的摩尔质量, kg/mol; p— — 孔隙压力, MPa; p0— — 泥页岩初始孔隙压力, MPa; pdf— — 钻井液液柱压力, MPa; ps— — 钻井液封堵层内孔隙压力, MPa; r— — 距井筒中心的距离, m; rw— — 井筒半径, m; R— — 通用气体常数, 8.314 J/(mol· K); t— — 时间, s; T— — 绝对温度, K; x— — 距裂缝中心点处x方向的距离, m; α — — Biot系数, 无因次; α d— — 残余膨胀压校正系数, 无因次; α π — — 扩散校正系数, 无因次; γ — — 损伤系数, m3/mol; Δ y(x)— — x处的裂缝宽度, m; δ i, j— — 系数, i=j时为1, ij时为0, 无因次; ${{\varepsilon }_{\text{v0, chem}}}$— — 平均水化体积应变量, 无因次; ${{\varepsilon }_{\text{v, chem}}}$— — 水化体积应变量, 无因次; η — — 钻井液黏度, mPa· s; μ — — 化学势, J/mol; μ k— — 孔隙流体中k组分的化学势, J/mol; ν — — 泊松比, 无因次; π d— — 黏土矿物层间水化膨胀压, MPa; ${{\rho }_{\text{f}}}$— — 钻井液的密度, kg/m3; ${{\rho }_{\text{w}}}$— — 水的密度, kg/m3; ${{\sigma }_{\text{c}}}$— — 围压, MPa; ${{\sigma }_{\text{e, }i\text{, }j}}$— — 有效应力分量, MPa; ${{\sigma }_{\text{h}}}$— — 最小水平地应力, MPa; ${{\sigma }_{\text{he}}}$— — 水化应力对有效应力的贡献, MPa; ${{\sigma }_{\text{hf}}}$— — 裂缝内水化膨胀压力, MPa; ${{\sigma }_{\text{H}}}$— — 最大水平地应力, MPa; ${{\sigma }_{i\text{, }j}}$— — 地应力张量的分量, MPa; ${{\sigma }_{\text{m}}}$— — 泥页岩的膜效率, 无因次; ${{\sigma }_{\text{ms}}}$— — 泥页岩封堵层的膜效率, 无因次; ${{\sigma }_{\text{t}}}$— — 抗拉强度, MPa; ${{\sigma }_{\text{v}}}$— — 上覆岩层压力, MPa; $\sigma _{1}^{* }$— — 最大有效主应力, MPa; $\sigma _{3}^{* }$— — 最小有效主应力, MPa; $\sigma _{i\text{, }j}^{* }$— — 考虑损伤的有效应力, MPa; τ — — 孔隙的几何特征参数, 无因次; ϕ — — 孔隙度, %; ${{\phi }_{\text{s}}}$— — 封堵层的孔隙度, %; φ — — 泥页岩内摩擦角, (° ); ω — — 黏土微孔/层间渗透效率, 无因次; ${{\omega }_{\text{0}}}$— — 与溶质浓度变化率相关的水化膨胀系数, MPa· s· m3/mol; ${{\omega }_{\text{1}}}$— — 浓度相关损伤系数, m3/mol; ${{\omega }_{k}}$— — k组分化学势相关的水化膨胀系数, MPa· mol/J; ${{\omega }_{\text{k0}}}$— — 单一组分化学势相关的水化膨胀系数, MPa· mol/J。

参考文献
[1] MA T, CHEN P. A wellbore stability analysis model with chemical-mechanical coupling for shale gas reservoirs[J]. Journal of Natural Gas Science and Engineering, 2015, 26: 72-98. [本文引用:2]
[2] 张世锋, 邱正松, 汪海阁, . 泥页岩流-固-化耦合模型参数的蒙特卡洛估计[J]. 石油学报, 2016, 37(7): 921-929.
ZHANG Shifeng, QIU Zhengsong, WANG Haige, et al. Monte Carlo estimation on the parameters of shale chemo-poro-mechanical coupling model[J]. Acta Petrolei Sinica, 2016, 37(7): 921-929. [本文引用:3]
[3] 邱正松, 徐加放, 吕开河, . “多元协同”稳定井壁新理论[J]. 石油学报, 2007, 28(2): 117-119.
QIU Zhengsong, XU Jiafang, LYU Kaihe, et al. A multivariate cooperation principle for well-bore stabilization[J]. Acta Petrolei Sinica, 2007, 28(2): 117-119. [本文引用:1]
[4] SUN L, HIRVI J T, SCHATZ T, et al. Estimation of montmorillonite swelling pressure: A molecular dynamics approach[J]. The Journal of Physical Chemistry C, 2015, 119(34): 19863-19868. [本文引用:1]
[5] 康毅力, 杨斌, 李相臣, . 页岩水化微观作用力定量表征及工程应用[J]. 石油勘探与开发, 2017, 44(2): 301-308.
KANG Yili, YANG Bin, LI Xiangchen, et al. Quantitative characterization of micro forces in shale hydration and field applications[J]. Petroleum Exploration and Development, 2017, 44(2): 301-308. [本文引用:1]
[6] 薛华庆, 周尚文, 蒋雅丽, . 水化作用对页岩微观结构与物性的影响[J]. 石油勘探与开发, 2018, 45(6): 1075-1081.
XUE Huaqing, ZHOU Shangwen, JIANG Yali, et al. Effects of hydration on the microstructure and physical properties of shale[J]. Petroleum Exploration and Development, 2018, 45(6): 1075-1081. [本文引用:1]
[7] HEIDUG W K, WONG S W. Hydration swelling of water-absorbing rocks: A constitutive model[J]. International Journal for Numerical and Analytical Methods, 1996, 20(6): 403-430. [本文引用:7]
[8] GHASSEMI A, TAO Q, DIEK A. Influence of coupled chemo-poro-thermoelastic processes on pore pressure and stress distributions around a wellbore in swelling shale[J]. Journal of Petroleum Science and Engineering, 2009, 67: 57-64. [本文引用:2]
[9] ZHOU X X, GHASSEMI A. Finite element analysis of coupled chemo-poro-thermo-mechanical effects around a wellbore in swelling shale[J]. International Journal of Rock Mechanics & Mining Sciences, 2009, 46: 769-778. [本文引用:2]
[10] ZHANG S F, SHENG J. Effects of salinity and confining pressure on hydration-induced fracture propagation and permeability of Mancos shale[J]. Rock Mechanics and Rock Engineering, 2017, 50(11): 2955-2972. [本文引用:5]
[11] ZHANG S F, SHENG J. Effect of water imbibition on fracture generation in mancos shale under isotropic and anisotropic stress conditions[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2018, 144(2): 1-10. [本文引用:7]
[12] AL-BAZALI T. A novel experimental technique to monitor the time-dependent water and ions uptake when shale interacts with aqueous solutions[J]. Rock Mechanics and Rock Engineering, 2013, 46(5): 1145-1156. [本文引用:2]
[13] CHENG J Y, WAN Z J, ZHANG Y D, et al. Experimental study on anisotropic strength and deformation behavior of a coal measure shale under room dried and water saturated conditions[J]. Shock and Vibration, 2015: 1-13. [本文引用:2]
[14] WONG R. Swelling and softening behavior of La Biche shale[J]. Canadian Geotechnical Journal, 1998, 35(2): 206-221. [本文引用:3]
[15] BOURG I C, SPOSITO G. Connecting the molecular scale to the continuum scale for diffusion processes in smectite-rich porous media[J]. Environmental Science & Technology, 2010, 44(6): 2085-2091. [本文引用:2]
[16] TRIPATHY S, SCHANZ T. Compressibility behavior of clays at large pressures[J]. Canadian Geotechnical Journal, 2007, 44(3): 355-362. [本文引用:1]
[17] ZOU W L, YE J B, HAN Z, et al. Effect of montmorillonite content and sodium chloride solution on the residual swelling pressure of an expansive clay[J]. Environmental Earth Sciences, 2018, 77(19): 1-12. [本文引用:2]
[18] ALTMANN S, TOURNASSAT C, GOUTELARD F, et al. Diffusion-driven transport in clayrock formations[J]. Applied Geochemistry, 2012, 27(2): 463-478. [本文引用:2]
[19] MALUSIS M A, SHACKELFORD C D. Coupling effects during steady-state solute diffusion through a semipermeable clay membrane[J]. Environmental Science & Technology, 2002, 36(6): 1312-1319. [本文引用:1]
[20] CAO W G, LI P, ZHAO M H. On statistical damage constitutive model and its parameters for rock based on normal distribution[J]. Chinese Hydrogeology & Engineering Geology, 2008, 32(3): 11-14. [本文引用:1]
[21] LI S, XU J, TAO Y, et al. Study on damage constitutive model of rocks based on lognormal distribution[J]. Journal of Coal Science & Engineering, 2007, 13(4): 430-433. [本文引用:1]
[22] LIU C, YANG L, CAO W. A statistical damage softening constitutive model for rock and back analysis of its parameters[J]. Chinese Journal of Underground Space and Engineering, 2007, 3(3): 453-457. [本文引用:1]
[23] WANG Z L, LI Y C, WANG J G. A damage-softening statistical constitutive model considering rock residual strength[J]. Computers and Geosciences, 2007, 33(1): 1-9. [本文引用:3]
[24] WEIBULL W. A statistical distribution function of wide applicability[J]. Journal of Applied Mechanics, 1951, 18(3): 293-297. [本文引用:2]
[25] 张慧梅, 雷利娜, 杨更社. 温度与荷载作用下岩石损伤模型[J]. 岩石力学与工程学报, 2014, 33(2): 3391-3396.
ZHANG Huimei, LEI Lina, YANG Gengshe. Damage model of rock under temperature and load[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(2): 3391-3396. [本文引用:3]
[26] LEMAITRE J. A course on damage mechanics[M]. 2nd ed. Berlin: Springer Science & Business Media, 2012. [本文引用:1]
[27] YANG Y, APLIN A C. A permeability-porosity relationship for mudstones[J]. Marine and Petroleum Geology, 2010, 27(8): 1692-1697. [本文引用:1]
[28] LOON L R V, BAEYENS B, BRADBURY M H. Diffusion and retention of sodium and strontium in opalinus clay: Comparison of sorption data from diffusion and batch sorption measurements, and geochemical calculations[J]. Applied Geochemistry, 2005, 20(12): 2351-2363. [本文引用:1]
[29] GARCÍAGUTIÉRREZ M, CORMENZANA J L, MISSANA T, et al. Diffusion experiments in Callovo-Oxfordian clay from the Meuse/ Haute-Marne URL, France. Experimental setup and data analyses[J]. Physics & Chemistry of the Earth Parts B, 2008, 33: 125-130. [本文引用:1]