2. 原油物性与温降的耦合关系

在原油加热输送的热力计算中，最容易出错的地方往往不是公式本身，而是把密度、比热容、导热系数和粘度当成常数代入。实际上，这四个物性都随温度变化，而温度又沿管线不断下降，于是形成一条“温度降低 → 物性改变 → 热力与水力结果改变”的耦合链。理解这条链，是后续温降计算和保温设计能否贴近实际的前提。

四个物性随温度的变化规律

密度随温度升高而减小，对多数原油在输送温度区间内近似线性，可用 ρ(T)=ρ₂₀[1−α(T−20)] 估算，体积膨胀系数 α 大致在 6×10⁻⁴～9×10⁻⁴ /℃ 量级。密度本身对温降方程影响不大，但它决定质量流量与体积流量的换算，也影响自然对流与散热计算中的格拉晓夫数。

比热容随温度升高略有增大，且对含蜡、含胶质沥青质的原油更敏感：在析蜡点附近，蜡结晶释放潜热，使表观比热容出现明显峰值。工程上常取输送温度区间的平均值，例如 1.9～2.2 kJ/(kg·℃)；若跨越析蜡点，用单一平均值会低估该段的“放热能力”，导致温降算得偏快。

导热系数随温度升高而略降，原油一般在 0.12～0.16 W/(m·℃)。它主要影响管壁与保温层的导热热阻，对原油本体温降的直接作用较小，但在计算总传热系数 K 时不可忽略，尤其是保温层导热系数随温度上升而增大的特性更值得关注。

粘度是四者中温度敏感性最强的，常呈指数规律，可用 ν(T)=ν₀·e^(−b(T−T₀)) 或 ASTM 双对数式拟合。以某含蜡原油示例（合理假设值）为例：50 ℃ 时粘度约 30 mPa·s，40 ℃ 时约 60 mPa·s，30 ℃ 时约 150 mPa·s，接近凝点时急剧上升。粘度的剧烈变化会通过雷诺数、摩阻系数放大到压降上，并间接影响沿程摩擦生热。

物性如何进入热力计算

稳态温降的苏霍夫公式可写成 T(x)=T₀+(T起−T₀)·e^(−KπDx/(Gc))，其中 G 为质量流量、c 为比热容、K 为总传热系数。可以看出：比热容 c 直接出现在指数分母上，c 越大，同样热损失下温降越慢；K 值由原油导热系数、管壁、保温层和环境外膜共同决定，保温层导热系数随温度升高而增大，会使高温段 K 偏大、温降偏快。粘度不直接出现在温降公式里，但它通过压降改变沿程摩阻生热，在长距离、大流量管线上这部分热量可占总热损失的百分之几到十几，忽略它会使计算的终点温度偏低。

更关键的是，上述公式默认 c 和 K 为常数。当管线跨越析蜡点或凝点附近时，物性变化剧烈，必须分段计算：把全线按 5～10 ℃ 的温差或按里程分成若干段，每段用该段平均温度下的物性重新计算 K 和 c，再逐段推进。这就是“耦合”二字的实际含义——温度决定物性，物性又决定下一段的温度。

一个可复算的示例

设某段输油管：外径 D=0.325 m，质量流量 G=40 kg/s，起点油温 55 ℃，环境温度 5 ℃，总传热系数取 K=1.6 W/(m²·℃)。若全线统一取 c=2.0 kJ/(kg·℃)，则单位长度温降参数 KπD/(Gc)=1.6×3.1416×0.325/(40×2000)=4.08×10⁻⁵ /m，即每公里温降约 4.08×10⁻² ℃ 的指数衰减系数。若该原油在 40 ℃ 以下比热容因析蜡升至 2.4 kJ/(kg·℃)，则低温段衰减系数降为约 3.4×10⁻⁵ /m，温降明显放缓。用统一 c=2.0 计算会高估低温段温降，导致终点温度算低、加热站间距设计偏保守。

反过来，若忽略粘度变化，用 55 ℃ 的低粘度算全线压降，会低估压降和泵功率；而实际低温段粘度成倍上升，摩阻增大，若泵扬程余量不足，就会出现“温降算得差不多、实际却输不动”的情况。因此工程上常把热力计算与水力计算交替迭代：先假定温度分布算粘度和压降，再用压降算摩阻生热修正温度分布，直到前后两次温度偏差小于 0.5 ℃ 为止。

归纳起来，密度影响流量换算与对流换热，比热容主导温降幅度，导热系数通过 K 值影响散热强度，粘度则通过压降与摩阻生热间接参与热平衡。四者中粘度最敏感、比热容最直接、导热系数最“隐蔽”。在后续的温降计算与保温层设计中，凡是跨越析蜡点、凝点或大温差工况，都应坚持分段取值、热水平衡迭代，这样才能让计算结果既安全又经济。