3. 管柱与井壁接触力学基础

管柱在井眼中的摩阻与扭矩,本质上来源于管柱与井壁之间的接触摩擦。要把后续的摩阻扭矩计算做准,必须先厘清四个基本量的物理含义及其相互关系:接触正压力、摩擦系数、接触面积和摩阻力。本节从库仑摩擦的基本模型出发,说明各量的定义、影响因素、适用条件,并给出工程上常用的简化处理方式。

一、库仑摩擦模型与四个基本量

在钻井工程中,管柱与井壁的摩擦通常按干摩擦或边界润滑条件下的库仑摩擦处理,其基本关系为:

F_f = μ · N
T_f = μ · N · r

其中:
F_f —— 切向摩阻力(沿管柱轴向或滑动方向),N
T_f —— 摩擦扭矩(绕管柱轴线),N·m
μ   —— 摩擦系数(无量纲)
N   —— 接触处的法向正压力(合力),N
r   —— 摩擦力作用点到管柱轴线的力臂,m

这里需要强调一个容易混淆的点:库仑摩擦中真正进入公式的是接触正压力 N,而不是接触面积。这是干摩擦定律的经典结论——摩擦力与法向载荷成正比,与表观接触面积无关。但这一结论有明确的适用前提,后文会专门讨论。

二、接触正压力的来源与计算

接触正压力是管柱因自重、弯曲、张力等因素压向井壁时产生的法向反力。工程上最常用的估算方法是单位长度接触力(线接触载荷)乘以接触长度:

N = q_n · L

q_n ≈ √( (w·sinα)² + (2T·sin(θ/2)/L_c)² )

其中:
q_n —— 单位长度法向接触力,N/m
w   —— 管柱单位长度浮重,N/m
α   —— 井斜角
T   —— 该段管柱轴向张力,N
θ   —— 井眼狗腿角(该段两端井斜方位变化合成角),rad
L_c —— 狗腿段长度,m
L   —— 接触段长度,m

第一项 w·sinα 反映重力在垂直井壁方向的分量,在水平段最大、直井段为零;第二项反映张力在狗腿处的“绷弦”效应,张力越大、狗腿越急,正压力越大。这解释了为什么大位移井、水平井的摩阻问题远比直井突出,也解释了为什么解卡时盲目增大上提张力反而可能加剧摩阻。

工程提示:在直井段,若管柱居中且无狗腿,q_n 近似为零,理论上不产生接触摩阻;实际直井中的摩阻多来自井眼曲率、管柱弯曲和井壁不规则,需按实测摩阻反演修正。

三、摩擦系数的取值与影响因素

摩擦系数 μ 是摩阻扭矩计算中不确定性最大的参数。它并非材料常数,而是受接触面状态、钻井液润滑性、井壁泥饼、岩性、相对滑动速度、温度压力等共同影响的综合参数。工程上通常按工况分类取值:

工况典型 μ 范围说明
套管—套管(油基泥浆)0.15~0.20润滑性好,取值偏低
套管—套管(水基泥浆)0.20~0.30常见旋转下入工况
钻杆—裸眼(砂岩)0.25~0.35泥饼质量影响显著
钻杆—裸眼(页岩/泥岩)0.30~0.45易缩径、泥饼厚,偏高
卡瓦/键槽工况0.40 以上局部嵌入,需专门评估

需要说明:上表为工程经验参考区间,实际项目应通过摩阻反演(用起下钻实测大钩载荷反算 μ)进行标定,而不是直接套用。同一口井不同井段、不同工况的 μ 可能差异明显。

四、接触面积的作用与适用条件辨析

既然库仑摩擦与面积无关,为什么工程中还要关注接触面积?原因在于适用条件:

  • 干摩擦/边界润滑:真实接触面积由微凸体决定,与表观面积无关,F=μN 成立,这是钻井摩阻计算的主流假设。
  • 流体动压润滑或泥饼承载:当接触面间存在连续液膜或软泥饼时,法向载荷部分由流体/泥饼承担,摩擦行为偏离纯库仑模型,此时需引入等效摩擦系数或承载力修正。
  • 局部嵌入与塑性接触:卡瓦、键槽、缩径处接触应力超过岩石或管材屈服强度时,接触面积与压入深度相关,摩阻可能随载荷非线性增长。

因此,接触面积在接触应力校核(判断是否压溃、是否嵌入)和偏离库仑模型时的修正中才是关键量,而在常规摩阻扭矩累加计算中不直接进入公式。理解这一点,可以避免把“接触面积大所以摩阻大”这类错误直觉带入计算。

接触受力示意(井眼截面) 井壁 管柱 接触点 N(正压力) F_f = μN w·sinα 摩阻扭矩计算关系链 井眼轨迹参数 α、θ、L_c 管柱浮重 w、张力 T 接触状态参数 摩擦系数 μ 接触长度 L 单位长度接触力 q_n q_n = √[(w·sinα)² + (2T·sin(θ/2)/L_c)²] 接触正压力 N = q_n · L (沿接触段累加) 摩阻力 / 扭矩 F_f = μN T_f = μN·r 适用条件(库仑模型成立前提) ① 干摩擦或边界润滑,无连续动压液膜 ② 接触应力低于管材/岩石屈服,无塑性嵌入 ③ 摩擦系数按工况标定,非材料常数 ④ 接触面积不直接进入 F=μN,仅用于应力校核

五、从基本关系到工程计算的组织方式

把上述关系串起来,一段管柱的摩阻扭矩计算可按以下顺序组织:先由井眼轨迹(井斜角、狗腿角、段长)和管柱参数(浮重、张力)算出各段的单位长度接触力 q_n;再乘以该段实际接触长度得到正压力 N;然后按工况选取摩擦系数 μ;最后累加得到总摩阻和总扭矩。这个链条中,q_n 和 μ 是误差的主要来源,而 N 的累加方式(分段长度、是否考虑屈曲)决定了计算精度。

以一个水平段示例说明(以下为教学假设参数,非实测):某水平段长 100 m,井斜角 90°,管柱单位长度浮重 w = 300 N/m,该段张力 T = 200 kN,狗腿角 θ = 3°(0.0524 rad),狗腿段长 L_c = 30 m。则重力项 w·sinα = 300 N/m;绷弦项 2T·sin(θ/2)/L_c ≈ 2×200000×0.0262/30 ≈ 349 N/m;合成 q_n ≈ √(300² + 349²) ≈ 460 N/m;接触正压力 N ≈ 460×100 = 46 kN;取 μ = 0.25,则摩阻力 F_f ≈ 0.25×46 ≈ 11.5 kN。若管柱半径 r = 0.0635 m,对应扭矩约 0.73 kN·m。可见在水平段,绷弦项甚至可能超过重力项,张力控制对摩阻影响很大。

注意:上述算例仅用于说明计算流程与量级关系,实际工程必须以本井实测摩阻反演标定 μ,并考虑屈曲、井径变化、泥饼厚度等因素,不可直接引用示例数值作为设计依据。

六、小结与判断要点

掌握本节内容,关键是把四个量的角色分清:正压力 N 是“因”,摩擦系数 μ 是“状态”,摩阻力 F_f 是“果”,接触面积则是判断模型是否适用的“校核量”而非直接计算量。当现场摩阻异常偏大时,应优先排查 q_n 是否被低估(狗腿、张力、屈曲)和 μ 是否取值偏低(泥饼、岩性、润滑失效),而不是简单归因于“接触面积大”。只有在接触应力接近或超过屈服、或存在明显液膜承载时,才需要跳出库仑模型,引入接触面积相关的修正。