📚 课程目录第4章支撑剂在压裂液中的沉降速度理论

4. 支撑剂在压裂液中的沉降速度理论

支撑剂在裂缝中能否被携送到远端，很大程度上取决于它与压裂液之间的相对运动——也就是沉降速度。沉降速度是支撑剂输送评价的“基准量”：它决定了单颗粒在给定液体中下落的快慢，也决定了工程上需要多高的排量、多黏的液体、多细的砂比才能把支撑剂送到目标位置。本节从最基础的斯托克斯定律出发，逐步引入阻力系数修正与非牛顿流体模型，并说明各自的适用边界。

一、斯托克斯定律：低雷诺数下的解析解

设支撑剂颗粒为球形，直径 d_p、密度 ρ_p，压裂液密度 ρ_f、黏度 μ。颗粒在静止液体中自由沉降，重力与浮力之差被流体阻力平衡时达到匀速，该速度即沉降速度 v_s。当颗粒周围流动处于层流（蠕变流）状态、颗粒雷诺数 Re_p = ρ_f v_s d_p/μ 远小于 1 时，阻力服从斯托克斯阻力 F_D = 3πμd_pv_s，由此得到斯托克斯沉降速度：

v_s = (ρ_p − ρ_f) · g · d_p² / (18 μ)

其中：g = 9.81 m/s²；d_p 单位 m；μ 单位 Pa·s；v_s 单位 m/s

该式的物理含义非常清晰：沉降速度与粒径的平方成正比、与黏度成反比、与密度差成正比。举一个教学示例（假设参数，非实测）：取 40/70 目石英砂中值粒径 d_p = 0.212 mm = 2.12×10⁻⁴ m，ρ_p = 2650 kg/m³，压裂液 ρ_f = 1000 kg/m³、μ = 0.010 Pa·s（近似线性瓜胶稀溶液的量级），则

v_s = (2650 − 1000) × 9.81 × (2.12e-4)² / (18 × 0.010)
    ≈ 1650 × 9.81 × 4.49e-8 / 0.18
    ≈ 4.04e-3 m/s ≈ 4.0 mm/s

代入雷诺数校验：Re_p = 1000 × 4.04e-3 × 2.12e-4 / 0.010 ≈ 0.086，远小于 1，说明本例用斯托克斯式是自洽的。若把黏度降到 0.003 Pa·s，v_s 升至约 13.5 mm/s，Re_p 约 0.95，已接近斯托克斯区上限，误差开始显现——这正是需要阻力系数修正的原因。

二、阻力系数修正：从层流走向过渡区与紊流区

更一般的表达是把阻力写成阻力系数形式：F_D = C_D·(πd_p²/4)·(ρ_fv_s²/2)，平衡后得通用沉降速度式：

v_s = sqrt( 4 g d_p (ρ_p − ρ_f) / (3 C_D ρ_f) )

关键在 C_D 随 Re_p 的变化。层流区（Re_p < 1 左右）C_D = 24/Re_p，代回即还原为斯托克斯式；过渡区（约 1 < Re_p < 1000）常用 Schiller–Naumann 关联式 C_D = (24/Re_p)(1 + 0.15Re_p^0.687)，或工程上更简便的 C_D = 24/Re_p + 3/√Re_p + 0.34；牛顿区（Re_p > 约 1000）C_D ≈ 0.44，此时 v_s 与 d_p^0.5 成正比，不再依赖黏度。由于 C_D 依赖 v_s，而 v_s 又待求，实际计算需迭代：先按斯托克斯式估初值 → 算 Re_p 与 C_D → 用通用式更新 v_s → 重复至收敛（通常 3～5 次即可）。

此外，真实支撑剂并非光滑圆球。工程上用球度与圆度描述形状，并以形状修正系数（或等效沉降直径）修正 C_D：棱角状颗粒阻力更大、沉降更慢，若直接套用球形公式会高估沉降速度。对 40/70 目石英砂，球度常在 0.7～0.9 之间，属于必须考虑形状修正的范畴。

三、非牛顿流体中的沉降速度模型

压裂液多为非牛顿流体，最常用的是幂律模型 τ = K·γ̇ⁿ，其中 K 为稠度系数（Pa·sⁿ），n 为流性指数。n = 1 时退化为牛顿流体；n < 1 为剪切变稀（瓜胶、滑溜水体系常见）。对幂律流体中球形颗粒的斯托克斯型沉降，常用模型为：

v_s = [ (ρ_p − ρ_f) g d_p^(n+1) / (18 K X) ]^(1/n)

X 为与 n 相关的修正常数（如取 X ≈ (2^(n+1) · 3^((1-n)/2)) / ... 的近似形式，
工程上常按 n 查表或采用经验关联式确定）

不必死记常数，重要的是理解趋势：n 越小（剪切变稀越强），沉降速度对粒径和密度差的敏感度下降，而对 K 的依赖上升；也就是说，增稠（提高 K）比单纯改变密度差更能有效抑制沉降。对屈服应力型流体（如交联冻胶，可用 Herschel–Bulkley 模型描述），还存在“临界屈服”现象：当颗粒所受净重力小于屈服应力产生的支撑力时，颗粒根本不沉降，v_s = 0。这解释了为什么交联冻胶的携砂能力远强于同黏度的线性胶。

需要强调，上述幂律沉降式同样建立在“颗粒周围为稳态剪切、忽略壁面与颗粒群干扰”的假设上。实际裂缝中支撑剂浓度高，颗粒间相互干扰使群沉降速度显著低于单颗粒速度，工程上常用 Richardson–Zaki 型关系 v_群 = v_s(1 − φ)^m（φ 为固相体积分数）做浓度修正。

四、模型选择与适用条件对照

选错模型是沉降评价最常见的系统性误差来源。判断顺序应为：先算 Re_p 定流域，再看流体是否非牛顿，最后看颗粒形状与浓度是否需要修正。

模型	核心表达式	适用条件	典型误用
斯托克斯	v_s = Δρ·g·d_p²/(18μ)	牛顿流体、球形、Re_p ≪ 1	粗粒径或低黏液中高估精度
阻力系数通用式	v_s = √(4gd_pΔρ/(3C_Dρ_f))	牛顿流体、全流域，需迭代	忽略形状修正
幂律沉降式	v_s ∝ (Δρ·g·d_pⁿ⁺¹/K)^1/n	幂律流体、Re_p 较小、单颗粒	用于交联冻胶（有屈服应力）
屈服型判据	净重力 < 屈服支撑力时 v_s=0	Herschel–Bulkley 类冻胶	仍按幂律式算得非零速度

实操建议：现场快速估算时，先用斯托克斯式算一个下限参考值，再用 Re_p 判断是否需要切换到阻力系数迭代；对交联体系务必先判断是否存在屈服应力，否则会把“不沉降”误判为“缓慢沉降”。

注意单位一致性：d_p 用米、μ 用 Pa·s、K 用 Pa·sⁿ，若混用 mm、mPa·s 会带来数量级错误；幂律式中 n 出现在指数上，单位换算不能简单按线性比例处理。

把上述判断落到一次具体的携砂设计上：若目标是把 40/70 目砂在 1000 kg/m³、10 mPa·s 的线性胶中输送，按斯托克斯式得 v_s ≈ 4 mm/s，Re_p ≈ 0.09，模型自洽；若改用交联冻胶，则必须先判断屈服应力是否足以托住颗粒，此时“沉降速度”可能直接为零，评价重点应从“算速度”转为“算临界支撑条件”。模型不是越复杂越好，而是要与实际流变参数、粒径范围和浓度条件匹配，并在报告中写明所选模型与 Re_p 校验结果，才能保证沉降评价可复现、可比较。