4. 摩阻系数与流态判别
在天然气输气管道压降计算中，沿程摩阻系数 λ（部分文献记作 f）是决定水力损失大小的关键参数之一。它并非固定常数，而是随雷诺数 Re 和管壁相对粗糙度变化的量。工程上若把 λ 取错，压降结果可能偏差百分之几十，进而影响管径选型、压缩机站间距和输量核算。因此，本节的核心任务是：先算准雷诺数、判明流态，再选对摩阻系数公式，并掌握紊流区公式的迭代求解方法。
一、雷诺数的计算与流态判别
雷诺数表征惯性力与黏性力的比值，是判别流态的唯一无量纲依据。对于圆管满流，其定义为：
Re = (v · d) / ν = (4 · Q) / (π · d · ν)式中：v 为管内平均流速（m/s）；d 为管道内径（m）；ν 为天然气在操作温度、压力下的运动黏度（m²/s）；Q 为工况体积流量（m³/s）。注意 Q 必须换算到管内实际压力、温度条件下的工况体积流量，而不是标准状态流量，否则 Re 会严重失真。天然气的运动黏度随压力升高而减小、随温度升高而增大，高压输气管道中 ν 通常在 1×10⁻⁶～3×10⁻⁶ m²/s 量级，具体应查取对应工况下的物性数据。
流态判别的一般界限为：Re < 2000 为层流；2000 ≤ Re ≤ 4000 为过渡区（临界区）；Re > 4000 为紊流。输气管道绝大多数运行在紊流区，Re 常在 10⁶～10⁸ 量级，因此层流公式在实际干线中很少用到，但过渡区与紊流区的区分、以及紊流内部“水力光滑—混合摩擦—完全粗糙”三个子区的划分，才是选公式的真正难点。
二、常用摩阻系数公式及其适用性
层流区摩阻系数只与 Re 有关，公式为 λ = 64/Re，与粗糙度无关，计算直接。紊流区则复杂得多：水力光滑区可用 Blasius 型公式（如 λ = 0.3164/Re^0.25，适用于 Re < 10⁵）或 Prandtl 光滑管公式；完全粗糙区 λ 仅取决于 ε/d，可用 λ = 1/(2·lg(3.7d/ε))² 直接计算；而混合摩擦区 λ 同时受 Re 与 ε/d 影响，最经典的表达就是 Colebrook-White 公式：
1/√λ = −2 · lg( ε/(3.7d) + 2.51/(Re·√λ) )该式由光滑管与粗糙管两项叠加而成，物理意义清晰、精度高，被公认为紊流过渡区的基准公式。但它的致命特点是 λ 同时出现在等式两边，属于隐式方程，无法直接解出，必须迭代或查 Moody 图。工程计算中应优先用迭代法求数值解，Moody 图仅作定性核对。
三、Colebrook-White 公式的迭代求解
迭代的基本思路是把隐式方程改写为不动点形式：先假定一个初值 λ₀，代入右端算出一个新值，再回代，直到相邻两次结果之差小于允许误差。由于右端含 √λ，直接迭代收敛较慢，工程上更常用的是把方程改写为关于 1/√λ 的迭代式：
令 x = 1/√λ
x_(k+1) = −2 · lg( ε/(3.7d) + 2.51·x_k / Re )
收敛判据：|x_(k+1) − x_k| < 1×10⁻⁶
最后 λ = 1 / x²初值可取完全粗糙区公式的结果 x₀ = −2·lg(ε/(3.7d))，或直接取 λ₀ = 0.02。下面给出一个教学示例（参数为合理假设，非实测）：内径 d = 0.6 m，绝对粗糙度 ε = 0.05 mm，工况下 Re = 2.0×10⁷。先算 ε/(3.7d) = 0.00005/(3.7×0.6) ≈ 2.25×10⁻⁵。取 x₀ = −2·lg(2.25×10⁻⁵) ≈ 9.30，则 λ₀ ≈ 0.01156。第一次迭代：2.51×9.30/2.0×10⁷ ≈ 1.17×10⁻⁶，与 2.25×10⁻⁵ 相加后取对数，得 x₁ ≈ 9.19；第二次 x₂ ≈ 9.19，已基本稳定，λ ≈ 1/9.19² ≈ 0.01184。可见初值选得好时，两三次迭代即可收敛。
四、工程取值中的几个关键判断
第一，粗糙度 ε 的取值往往比公式选择更影响结果。新无缝钢管、内涂层管可取 0.03～0.05 mm；普通碳钢新管约 0.05～0.1 mm；长期运行、有内腐蚀或积液污染的管道应取 0.2～0.5 mm 甚至更大。若沿用新管粗糙度去核算老管道，压降会被明显低估。第二，判断是否进入完全粗糙区很有价值：一旦 Re 超过 665 − 765·lg(ε/d)，λ 与 Re 无关，此时增大输量不再改善摩阻系数，压降近似与流量平方成正比，这一结论对输量调节和增输改造的判断非常直接。第三，过渡区（2000～4000）流态不稳定，工程设计中应尽量避免管道长期在此区间运行，若无法避开，可按紊流公式偏保守取值并留足压降裕量。
实际计算时建议形成固定顺序：由工况流量、内径、黏度算 Re → 结合 ε/d 判流态与子区 → 层流用 64/Re、完全粗糙区用显式公式、混合摩擦区用 Colebrook-White 迭代 → 将 λ 回代压降公式。迭代收敛后还应做一次合理性核对：λ 是否落在 0.008～0.03 的常见区间内，与同管径、同粗糙度量级的经验值是否接近。若偏差过大，优先复查 Re 的流量是否误用了标况值、ν 是否取错工况、ε 是否与实际管况不符，而不是反复调整迭代初值。
掌握本节方法后，摩阻系数就不再是一个“查表取值”，而是一条由工况参数出发、经流态判别和公式选择、最终由迭代收敛确定的计算链。这条链的每一步输入都必须与管道真实运行状态对应，压降计算结果才具备工程可信度。