2. 原油物性参数获取与处理

水力计算的所有结果——沿程摩阻、所需扬程、泵轴功率、温降分布——最终都建立在原油物性参数之上。物性取错,后面的计算再精细也没有意义。本节围绕密度、粘度、凝点、比热容四类关键物性,说明它们随温度变化的规律、在水力计算中的具体用法,以及工程上常见的数据来源与处理方式。

一、四类关键物性及其温度依赖

密度。原油密度随温度升高而下降,常用线性近似:ρ(t)=ρ₂₀[1−α(t−20)],其中ρ₂₀为20℃标准密度,α为体积膨胀系数,轻质原油约0.0007~0.0009 ℃⁻¹,重质原油可达0.0009~0.0010 ℃⁻¹。密度直接进入伯努利方程的重位压头项,也用于质量流量与体积流量的换算。长输管道沿线温度变化大时,若全程用20℃密度,重油管道的高差压降会产生可观偏差。

粘度。粘度是水力计算中影响最大的物性,摩阻与粘度呈近似一次方关系(层流)或0.25次方关系(紊流光滑区)。原油粘度随温度变化剧烈,常用两种模型:ASTM D341 粘温方程 lg lg(ν+0.7)=A−B·lg(T+273.15),适合宽温域拟合;以及 Arrhenius 型指数式 ν=ν₀e^(−bT),适合窄温域快速估算。含蜡原油在析蜡点附近粘度会陡增,出现非牛顿特性,此时不能简单套用牛顿流体公式。

凝点。凝点决定管输温度下限与是否需加热输送。含蜡原油凝点通常为0~35℃,管输时一般要求油温高于凝点3~5℃以上,否则蜡晶析出、粘度骤升、甚至出现凝管风险。凝点数据用于确定加热站出口温度、允许最低进站温度以及停输允许时间。

比热容。比热容用于热力计算:加热站热负荷 Q=G·c·ΔT,以及沿程温降的苏霍夫公式。原油比热容一般为1.8~2.2 kJ/(kg·℃),随温度升高略增,含蜡原油在析蜡区间因潜热释放会出现表观比热容峰值,工程上常按温度分段取值。

原油物性—温度关系与水力计算数据链 物性随温度变化趋势(示意) 温度 t / ℃ 物性值 粘度 ν 密度 ρ 比热容 c 凝点 T凝 析蜡点附近粘度陡增 物性数据在水力计算中的使用链 数据来源 化验报告 · 标准曲线 · 经验关联式 · 在线仪表 按管段温度分段取值 出站温度 → 沿程温降 → 进站温度 水力/热力计算 ρ → 重位压头、流量换算 ν → 雷诺数、摩阻系数 λ c → 热负荷、温降;T凝 → 温度下限 泵站扬程 · 加热负荷 · 允许停输时间 校核:油温 ≥ 凝点 + 3~5℃

二、物性数据在水力计算中的使用方法

关键在于“分段取值”而非“全程一个值”。长输管道沿线油温从出站70℃左右逐步降至进站30℃上下,粘度可能相差数倍。正确做法是沿管线按里程划分若干计算段,每段取平均油温对应的物性,逐段计算摩阻后累加。以某Φ529×8 mm、年输量800万吨的含蜡原油管道为例(示例参数):出站温度65℃,进站温度32℃,若全程按出站温度取粘度,摩阻将被低估约20%~30%,泵站扬程选型偏小;反之按进站温度取值则偏保守、浪费投资。

具体操作上,先由热力计算得到沿线温度分布 t(x),再对每一计算段取 t̄=(t进+t出)/2,代入粘温方程求 ν(t̄),计算该段雷诺数 Re=4Q/(πdν),判断流态后选取摩阻系数公式,得到该段水力坡降 i,最后 Σi·L 即为全线摩阻。密度同理按段取平均温度对应值,用于重位压头与流量换算。比热容则直接进入温降计算,与水力计算交替迭代——因为温降依赖流量,而流量对应的摩阻又依赖温度,工程上通常先假定温度分布、算水力、再修正温度,迭代2~3次即可收敛。

含蜡原油在析蜡点以下呈非牛顿特性(如宾汉姆流体),表观粘度随剪切率变化。此时应使用流变仪测得的流变曲线,按管段实际剪切率取值,不能直接套用牛顿流体的粘温方程。

三、常见数据来源与处理要点

工程上物性数据主要有四个来源:一是原油化验报告,包括20℃密度、凝点、比热容及若干温度点的粘度,这是最权威的依据;二是标准粘温曲线,由化验数据按ASTM D341拟合得到,可外推到未测温度点;三是经验关联式,在缺乏实测数据时按原油类别估算,精度有限,仅用于方案阶段;四是管道在线密度计、温度变送器等仪表数据,用于运行期校核。

处理时需注意三点:第一,统一温度基准,化验报告常给20℃密度,计算中要换算到实际油温;第二,粘度单位统一,化验常用 mm²/s(cSt),计算时须与流量、管径单位匹配,避免量纲错误;第三,对异常点做甄别,化验数据中个别温度点粘度偏离拟合曲线过大时,应复查化验过程而非直接采用。建议建立一张“温度—密度—粘度—比热容”对照表,作为全线分段计算的输入底表,既便于复核,也便于运行期与实测数据比对。

切勿用不同原油、不同批次的物性数据混算同一条管道。原油物性随产地和批次差异显著,混用会导致摩阻与温降计算结果失真,进而影响泵站选型与加热负荷配置。