4. 热平衡与温降基本方程
要定量回答“原油从首站到末站会降多少度”，必须把沿管道轴向取一个微元段，对这段油流做能量守恒分析。轴向温降微分方程正是这一分析的结果，它把散热损失、摩擦生热、焦耳-汤姆逊效应和流动携带的热量统一在一个式子里。理解这个方程的建立过程和每一项的物理意义，是后续解析求解、数值计算和保温设计的基础。
一、微元段热平衡的建立
取管道轴向坐标 x 自起点（首站出站）指向终点，管外土壤温度记为 T₀，油流温度记为 T(x)。在距起点 x 处截取长度 dx 的微元段，设质量流量为 G（kg/s），比热容为 c（J/(kg·℃)），管径为 D，总传热系数为 K（W/(m²·℃)）。
在稳态条件下，微元段内油流的内能不变，进入与离开微元段的热量必须平衡。单位时间内：
- 流入带进的热量：G·c·T(x)；
- 流出带走的热量：G·c·T(x+dx)；
- 经管壁向环境散失的热量：K·πD·dx·(T−T₀)；
- 摩擦生热：单位长度摩擦压降为 dp/dx（Pa/m），微元段内摩擦功率为 G·(dp/dx)·dx/ρ，折算为温升项；
- 焦耳-汤姆逊效应：原油压缩性很小，该项通常可忽略，但在高压天然气或含气原油中需保留。
将上述各项代入能量守恒并整理，得到轴向温降微分方程的标准形式：
G·c·(dT/dx) = −K·πD·(T − T₀) + G·(dp/dx)/ρ
即：dT/dx = −[K·πD/(G·c)]·(T − T₀) + (1/(ρ·c))·(dp/dx)若忽略摩擦生热（长距离、低流速或高黏度工况下摩擦项相对较小），方程退化为最常见的简化形式：
dT/dx = −a·(T − T₀)，其中 a = K·πD/(G·c)参数 a 的量纲为 1/m，其倒数 1/a 具有长度量纲，称为“温降特征长度”，反映油温向环境温度逼近的快慢：a 越大（K 大、流量小、管径大），温降越快。
二、各项物理意义
方程左边 dT/dx 是油温沿轴向的变化率，负值表示降温。右边第一项 −a(T−T₀) 是散热驱动项：油温高于环境温度时持续放热，温差越大散热越快，且随温差减小而自动减弱，这正是温降曲线呈指数衰减的根本原因。第二项 (1/(ρc))·(dp/dx) 是摩擦生热项：原油沿管道流动克服摩阻，机械能转化为内能，使油温略有回升，压降越大、流速越高，该项越显著。对于高含蜡、高黏原油的热力计算，摩擦项在站间后半段可能不可忽略，工程上常以“摩擦生热修正系数”或直接保留该项参与数值积分。
三、稳态与非稳态的适用条件
上述方程是在稳态假设下建立的，即流量、油温、环境温度、传热系数均不随时间变化。实际运行中，这一假设是否成立决定了能否直接使用解析解。
| 工况类型 | 适用条件 | 方程形式与处理方式 |
|---|---|---|
| 稳态运行 | 流量稳定、出站温度恒定、地温与总传热系数基本不变，连续运行时间远大于热惯性时间 | 可直接用 dT/dx = −a(T−T₀) 解析求解，得到指数温降公式 |
| 启动/停输 | 油温随时间明显变化，管壁与土壤蓄热尚未建立稳定温度场 | 需引入时间项，采用非稳态方程 ∂T/∂t + v·∂T/∂x = −a(T−T₀)，配合数值方法求解 |
| 变流量/变出站温度 | 输量调整、加热炉切换、间歇输送 | 按准稳态分段处理或直接数值积分，注意 a 随流量变化 |
| 停输再启动 | 停输期间油温按径向和轴向同时下降，存在析蜡风险 | 必须用非稳态模型，重点校核最低油温与再启动压力 |
判断能否按稳态处理，可用“热惯性时间”作粗略估计：土壤温度场建立通常需要数天到数周，若工况变化持续时间远小于该时间尺度，则应按非稳态处理；若管道已连续稳定运行较长时间，且流量、出站温度波动在允许范围内，稳态方程足够用于日常温降核算与保温设计。
四、一个教学示例
假设某段管道：G = 50 kg/s，c = 2000 J/(kg·℃)，D = 0.4 m，K = 1.5 W/(m²·℃)，T₀ = 5 ℃，出站油温 T₁ = 60 ℃。计算 a = K·πD/(G·c) = 1.5×3.1416×0.4/(50×2000) ≈ 1.885×10⁻⁵ 1/m，特征长度 1/a ≈ 53 km。若忽略摩擦项，则任意距离 x 处油温为 T(x) = T₀ + (T₁−T₀)·e^(−a·x)。在 x = 30 km 处，T ≈ 5 + 55×e^(−0.5655) ≈ 5 + 55×0.568 ≈ 36.2 ℃。这一结果说明：在给定参数下，30 km 内油温下降约 24 ℃，温降曲线前段陡、后段缓，符合指数规律。上述数值为教学假设，实际工程应以现场物性、地温和传热系数测试值为准。
需要强调的是，稳态方程给出的是“运行稳定后”的轴向温度分布，它不能回答“停输 8 小时后油温降到多少”这类问题。工程上做保温设计时，通常先用稳态方程确定正常输送所需的保温层厚度和出站温度，再用非稳态模型校核停输、再启动等极端工况下的最低油温是否高于凝点或析蜡点。两者配合，才能既满足正常运行的经济性，又保证特殊工况的安全性。
综上，轴向温降微分方程的本质是“流动携带的热量 = 散热损失 + 摩擦生热”的微元表达。掌握每一项的来源和量级，明确稳态与非稳态的边界条件，才能在温降计算与保温设计中选对模型、用对参数、判对结果。