3. 管道水力计算基础方程与适用条件

原油长输管道的水力计算，本质上是在回答三个相互关联的问题：油品在管内的流量与流速如何对应、沿程能量如何转换与损失、以及压降与流量之间服从什么定量规律。连续性方程、伯努利方程和达西-魏斯巴赫公式分别承担这三项任务，它们不是彼此独立的公式，而是一条从"质量守恒"到"能量守恒"再到"摩擦损失定量化"的递进链条。理解这条链条的物理意义、适用前提和单位制统一方法，是后续摩阻系数选取、泵站扬程计算和工况校核的基础。

连续性方程：质量守恒的数学表达

对管道内稳定流动的原油，取两个过流断面1和2，连续性方程写作：

ρ₁·A₁·v₁ = ρ₂·A₂·v₂ = q_m
对于不可压缩流体（ρ₁ ≈ ρ₂ = ρ）：
A₁·v₁ = A₂·v₂ = q_V

式中 q_m 为质量流量（kg/s），q_V 为体积流量（m³/s），A 为过流面积（m²），v 为断面平均流速（m/s）。物理意义十分明确：稳定流动时，单位时间内通过任一断面的质量必然相等，流体不会在管内"堆积"或"凭空产生"。

对长输管道而言，原油虽有一定压缩性，但在常规压力范围内密度变化通常小于1%，工程上按不可压缩处理是合理的。由此得到一个非常实用的推论：管径不变则流速不变，管径变化则流速与面积成反比。例如某段管道由 Φ720×8 变为 Φ630×8，内径由704 mm变为614 mm，面积比约为1.315，若上游流速为1.2 m/s，则下游流速约为1.58 m/s。这一关系在变径、分输、支线接入的计算中反复使用。

连续性方程的适用前提是：稳定流（各断面参数不随时间变化）、单相连续介质、断面间无质量注入或抽出。若存在分输阀室或注入点，需按节点分段列方程，把分输量作为边界条件处理。

伯努利方程：能量守恒与压头转换

伯努利方程描述断面间机械能的转换与损失。以单位重量流体为基准（能量头形式），取断面1为上游、断面2为下游，沿流动方向列式：

z₁ + p₁/(ρg) + v₁²/(2g) = z₂ + p₂/(ρg) + v₂²/(2g) + h_f + h_j

三项依次为位置水头、压力水头、速度水头，h_f 为沿程水头损失，h_j 为局部水头损失。物理意义是：上游断面的总机械能等于下游断面总机械能加上两断面间被摩擦和局部阻力耗散掉的能量。若两断面间设有泵站，则在左侧（上游侧）加上泵扬程 H_p，这正是泵站选型与水力计算衔接的接口。

长输管道计算中，位置水头 z 的变化相对压力水头往往很小，但在地形起伏段不可忽略；速度水头在等径管段中因 v₁ = v₂ 而相互抵消，这也是等径长输管道常用"压降即水头损失"简化计算的原因。需要特别强调参照约定：本文及后续计算一律以流动方向为正方向，高程以统一基准面起算，泵扬程计入上游侧，损失项一律取正值置于下游侧。切换基准面或改变正方向时，必须在计算书中明确标注，否则极易出现符号错误。

伯努利方程的适用前提是：稳定流、不可压缩流体、断面间无外功输入输出（有泵或阀门时需补充相应项）、所取断面为渐变流断面。阀门、弯头、三通等局部阻力集中的位置不宜直接取作计算断面。

达西-魏斯巴赫公式：沿程摩阻的定量计算

沿程水头损失由达西-魏斯巴赫公式给出：

h_f = λ · (L/d) · v²/(2g)
或以压降形式表示：
Δp_f = λ · (L/d) · (ρv²/2)

式中 λ 为沿程阻力系数（无量纲），L 为管段长度（m），d 为管道内径（m），v 为平均流速（m/s），ρ 为密度（kg/m³）。公式结构说明：损失与管长成正比、与内径成反比、与速度平方成正比。这解释了长输管道为什么"增大管径降压降"效果显著——管径增大一倍，同等流量下流速降为四分之一，摩阻约降为十六分之一（尚未计入 λ 随雷诺数变化的影响）。

λ 不是常数，它取决于流态与管壁粗糙度：层流区 λ = 64/Re，与粗糙度无关；紊流光滑区、过渡区、粗糙区分别适用不同的经验或半经验关系（如柯列勃洛克公式、尼古拉兹公式等），这是后续章节专门讨论的内容。本节只需明确：达西-魏斯巴赫公式本身是定义式，λ 的取值才是计算精度的关键。

该公式的适用前提是：圆管、满管、稳定流、单相牛顿流体。对含蜡、含气或非牛顿特性的原油，需修正或改用相应流变模型。

单位制统一：工程计算中最常见的错误来源

三个方程联用时，单位不一致是最容易出错也最容易被忽视的环节。工程上常见混用：流量用 m³/h、直径用 mm、长度用 km、压力用 MPa、密度用 kg/m³。若直接代入公式而不换算，结果可能相差数个数量级。推荐做法是全部换算为 SI 基本单位后再计算：流量 m³/s、内径 m、长度 m、流速 m/s、密度 kg/m³、压力 Pa、水头 m。

常用换算关系如下表，建议在计算书首页固定列出，形成习惯：

一个典型示例：某管段长 50 km，内径 704 mm，输量 1200 m³/h，原油密度 860 kg/m³，运动粘度 25 cSt，λ 取 0.022。换算后 q_V = 0.3333 m³/s，d = 0.704 m，A = 0.3893 m²，v = 0.856 m/s，L = 50000 m。代入达西-魏斯巴赫公式得 h_f ≈ 0.022 × (50000/0.704) × 0.856²/(2×9.81) ≈ 58.4 m 油柱，对应压降约 0.49 MPa。此结果属于教学示例的合理估算，实际工程须以实测粘度、实际内径和校核后的 λ 重新计算。

把三个方程放在一起看，计算路径就清晰了：连续性方程由输量和管径确定流速，达西-魏斯巴赫公式由流速、管径、管长和 λ 算出沿程损失，伯努利方程再把损失、高程差和泵扬程统一到能量平衡中，求出各断面压力或所需扬程。任何一步的单位或符号错误，都会沿链条传递到最终结果。因此在实际计算中，建议按"统一单位—算流速—算损失—列能量方程—校核结果量级"的顺序推进，并对每一步结果做量级合理性判断，例如长输管道每公里压降通常在 0.01～0.03 MPa 量级，若算出每公里数兆帕，首先应检查单位换算而非怀疑公式本身。