4. 沿程摩阻系数确定方法
沿程摩阻系数 λ(也称达西摩阻系数)是达西—魏斯巴赫公式中的关键参数,它直接决定沿程水头损失的计算结果。在原油长输管道的工艺计算中,λ 取值的偏差往往被放大到全线摩阻和泵站扬程上,因此确定 λ 的方法必须与流态、管壁粗糙度和介质物性相匹配。本节围绕"先判流态、再选公式、最后校核"的主线展开。
一、先判别流态:层流与湍流的界线
流态判别依据雷诺数 Re = vd/ν,其中 v 为平均流速(m/s),d 为管道内径(m),ν 为原油运动黏度(m²/s)。工程上以 Re = 2000 作为层流与湍流的分界,Re < 2000 为层流,Re > 4000 为完全湍流,2000~4000 为过渡区。层流时 λ 只与 Re 有关,用哈根—泊肃叶公式 λ = 64/Re 即可精确计算;过渡区流态不稳定,工程上通常避免在该区间长期运行,计算时按湍流偏保守取值。
湍流区还需进一步区分水力光滑区、混合摩擦区和粗糙区(阻力平方区),判据是相对粗糙度 ε/d 与 Re 的组合。常用界限为:Re < 59.5/(ε/d)^(8/7) 属光滑区;超过该值进入混合摩擦区;当 Re > 665 − 765·lg(ε/d) 时进入粗糙区,此时 λ 仅取决于 ε/d,与 Re 无关。原油管道常用钢管绝对粗糙度 ε 取 0.05~0.1 mm(新管)至 0.2 mm(结蜡、腐蚀后),计算时应按实际运行状态选取。
二、莫迪图:λ 的图形化查取
莫迪图以 Re 为横坐标(对数刻度)、λ 为纵坐标,绘出不同 ε/d 下 λ 随 Re 变化的曲线族。使用时先算 Re 和 ε/d,在图上定位交点读取 λ。莫迪图的价值在于直观展示三个湍流分区的过渡关系:左侧陡降段为层流,中部曲线随 Re 增大而下降为光滑区与混合摩擦区,右侧水平段为粗糙区。工程估算时可先查图得到量级,再用公式精确计算。
三、常用经验公式对比
湍流区 λ 的显式计算以 Colebrook-White 公式精度最高,但它是隐式方程,需迭代求解:
1/√λ = −2·lg[ ε/(3.7d) + 2.51/(Re·√λ) ]
工程上常用 Haaland 公式作为显式近似,无需迭代,精度与 Colebrook-White 相差约 1%~2%:
1/√λ = −1.8·lg[ (ε/d / 3.7)^1.11 + 6.9/Re ]
此外,Blasius 公式 λ = 0.3164/Re^0.25 仅适用于光滑区(Re < 10⁵);粗糙区可用 λ = 1/[2·lg(3.7d/ε)]² 直接计算。选用原则是:层流用 64/Re;光滑区可用 Blasius 或 Haaland;混合摩擦区优先 Colebrook-White,快速估算用 Haaland;粗糙区用粗糙区专用式。对长输管道,若全线 Re 与 ε/d 变化不大,可统一用 Colebrook-White 迭代一次全线 λ;若各站间油温、黏度差异显著,应分段计算 Re 后分别取 λ。
实际工程中,建议将 Re、ε/d、所选公式及计算结果列入计算表逐段记录,便于复核与交接。当同一管段在不同季节或不同输量下跨越流态分区时,应以最不利工况(通常是低输量、低流速、高黏度)校核 λ,确保泵站选型留有余量。