4. 井壁稳定性力学与化学耦合分析基础

井壁失稳很少是单一因素造成的。钻开地层后，井眼周围原本平衡的地应力被打破，岩石骨架承担的有效应力重新分布；与此同时，钻井液滤液进入近井带，与黏土矿物发生水化、膨胀、分散或溶蚀，改变岩石的强度参数与变形特性。前者是力学问题，后者是化学问题，而真实井壁的坍塌或破裂，往往是两者叠加的结果。本节先建立井周应力与坍塌/破裂压力的基本算法，再说明化学作用如何"改写"强度参数，最后用一个简化井例把两条线索串起来。

一、井周应力分布：从地应力到井壁上的有效应力

分析前先约定符号与坐标系：取垂向主应力为 σ_v，两个水平主应力为 σ_H（最大）和 σ_h（最小），压应力为正；井眼轴线沿垂向，井壁位置用从最大水平主应力方向起算的周向角 θ 表示。在井壁处（井眼半径 r = r_w），弹性解给出的三个主应力为：

• 径向应力 σ_r = p_w（井筒液柱压力），是最小主应力方向；
• 切向应力 σ_θ = σ_H + σ_h − 2(σ_H − σ_h)cos2θ − p_w；
• 轴向应力 σ_z = σ_v − 2ν(σ_H − σ_h)cos2θ，ν 为泊松比。

由表达式可见，当 θ = 90° 或 270°（即最小水平主应力方向）时 cos2θ = −1，切向应力达到最大，该方位是压性剪切破坏（坍塌）最易发生的部位；而 θ = 0° 方向切向应力最小，是拉伸破裂的敏感方位。把总应力换算为有效应力需引入 Biot 系数 α：σ′ = σ − αp_p，p_p 为孔隙压力。井壁是否破坏，最终由有效应力与岩石强度比较决定。

二、坍塌压力与破裂压力的基本算法

坍塌压力对应井壁剪切破坏的临界液柱压力。采用 Mohr–Coulomb 准则，令最大、最小有效主应力为 σ′₁、σ′₃，强度参数为黏聚力 C 与内摩擦角 φ，破坏条件为 σ′₁ = σ′₃·tan²(45°+φ/2) + 2C·tan(45°+φ/2)。将 σ′_θ、σ′_r 代入并令其临界相等，可反解出维持井壁不坍塌所需的最小井筒压力，即坍塌压力 p_c。工程上常写成等效钻井液密度 ρ_c = p_c/(gH)（H 为垂深），便于与设计密度直接比较。

破裂压力则由拉伸破坏控制：当井壁最小有效主应力变为拉应力并达到岩石抗拉强度 T 时产生裂缝，通常发生在 θ = 0° 方位，令 σ′_θ = −T 即可解出 p_f。安全钻井液密度窗口即位于 p_c 与 p_f 之间，并需叠加孔隙压力与漏失压力约束。

三、化学作用如何改变强度参数

钻井液滤液侵入后，化学作用主要通过三条路径削弱岩石：其一，黏土水化使颗粒间胶结弱化，黏聚力 C 明显下降；其二，吸水膨胀产生附加膨胀应力，改变近井带有效应力状态；其三，离子交换与溶蚀改变颗粒接触性质，使内摩擦角 φ 略有降低。宏观表现就是 C、φ 随含水量或浸泡时间下降，坍塌压力相应上升——同一地层，用清水钻进比用抑制性钻井液更容易坍塌，原因正在于此。

因此耦合分析的关键，是把强度参数写成含水率 w 或浸泡时间 t 的函数：C = C(w)、φ = φ(w)，再代回坍塌压力公式，得到"化学损伤后的 p_c(t)"。抑制性钻井液的作用，本质上就是减缓 C、φ 随时间的衰减速率，从而把 p_c 控制在设计密度之下。

四、简单井例：把两条线索合起来算

设某页岩井垂深 H = 2500 m，上覆压力 σ_v = 55 MPa，水平主应力 σ_H = 48 MPa、σ_h = 38 MPa，孔隙压力 p_p = 25 MPa，泊松比 ν = 0.25，Biot 系数 α = 0.9。以上均为教学假设值，用于演示流程，不代表实测数据。

先看坍塌敏感方位 θ = 90°，切向总应力 σ_θ = 48 + 38 + 2×(48−38) − p_w = 106 − p_w（MPa）。取有效应力 σ′_θ = σ_θ − 0.9×25 = 83.5 − p_w，径向有效应力 σ′_r = p_w − 22.5。设新鲜岩心试验得 C = 6 MPa、φ = 28°，则 K = tan²(59°) ≈ 2.77。代入破坏条件 σ′_θ = K·σ′_r + 2C√K，即 83.5 − p_w = 2.77(p_w − 22.5) + 2×6×1.664，整理得 p_w ≈ 31.5 MPa，对应等效密度约 1.28 g/cm³。

再考虑化学损伤：若钻井液抑制性不足，浸泡 3 天后 C 降至 3.5 MPa、φ 降至 25°（K ≈ 2.46），同样代入可解出 p_w 升至约 35 MPa，等效密度约 1.43 g/cm³。这意味着原本 1.30 g/cm³ 的钻井液在新鲜状态下安全，浸泡后却已低于坍塌压力，井壁会在最小水平主应力方位出现剪切崩落。破裂压力一侧，令 θ = 0° 处 σ′_θ = −T（取 T = 2 MPa），可解出 p_f 约 60 MPa 量级，密度窗口上限约 2.4 g/cm³，说明本例中坍塌是主要矛盾，控制手段应集中在提高抑制性、适度提高密度，而非担心压裂漏失。

五、耦合分析的操作要点

实际做耦合分析时，建议按以下顺序推进：先用地应力与孔隙压力数据算出力学控制的 p_c、p_f 基线；再用岩心浸泡试验测不同时间点的 C、φ，拟合 C(t)、φ(t)；然后把时变强度代回公式，得到 p_c(t) 曲线；最后把设计钻井液密度与 p_c(t) 对比，若设计密度在浸泡后期低于 p_c(t)，就必须提高抑制剂加量或调整密度。判断依据不是单一时刻的数值，而是"设计密度是否始终高于整个浸泡周期内的 p_c(t)"。

把上述流程固化为可复用的计算表，有助于现场快速判断。表中每一行对应一个浸泡时间点，最后一列给出该时刻的安全结论，可直接用于钻井液配方比选。

这张表清楚显示：同一口井、同一密度，安全与否取决于浸泡时长与抑制剂效果。若换用抑制性更强的钻井液体系，使 3 天时 C 仅降至 5.2 MPa，则 p_c 可维持在 1.35 g/cm³ 附近，设计密度即可覆盖整个作业周期。力学给出窗口的边界，化学决定边界随时间移动的速度，二者共同决定井壁能否稳住。