2. 天然气物性与压缩基础

天然气压缩机站的流程设计与能耗分析，最终都要落到一组可计算的物性参数上。无论是估算压缩机的轴功率、校核管道的温降与压降，还是判断级间冷却器的负荷，输入端的天然气组成、状态方程、压缩因子和热值决定了计算结果是否可信。本节的任务不是罗列物性公式，而是把"现场气质数据"转化为"可用于流程与能耗计算的数据体系"，并说明每个参数在后续计算中扮演什么角色。

一、从组成到摩尔分数：物性计算的起点

天然气是典型的多组分混合物，主要成分为甲烷，通常还含有乙烷、丙烷、丁烷及更重的烃类，以及氮气、二氧化碳、硫化氢等非烃组分。工程上最常用的输入形式是摩尔分数（或体积分数，在理想气体近似下二者数值接近）。例如某站场气质分析给出的典型干气组成（示例假设值）为：甲烷 92.0%、乙烷 4.0%、丙烷 1.5%、正丁烷 0.5%、氮气 1.5%、二氧化碳 0.5%。这组数据是后续所有计算的源头，一旦组成偏差较大，压缩因子、热值和功率都会随之偏移。

由组成可直接计算两个基础量：视相对分子质量与视临界参数。视相对分子质量按各组分摩尔分数加权求和；视临界压力和视临界温度同样按 Kay 规则加权。以甲烷（M=16.04，Pc=4.60 MPa，Tc=190.6 K）与二氧化碳（M=44.01，Pc=7.38 MPa，Tc=304.2 K）为例，含 0.5% CO₂ 的混合气视临界压力会比纯甲烷略高、视临界温度略升，这正是后续对比压力、对比温度计算的依据。

二、状态方程与压缩因子：真实气体的核心修正

理想气体方程 PV=nRT 在低压常温下尚可近似，但压缩机站常见出口压力在 4～10 MPa 量级，此时气体分子间作用力和分子体积不可忽略，必须引入压缩因子 Z 修正：PV = ZnRT。Z 的物理含义是"实际气体体积与同温同压理想气体体积之比"，Z 偏离 1 越远，说明气体越偏离理想行为。

工程上最常用的状态方程是Standing-Katz 图版法及其拟合式，以对比压力 Pr=P/Pc 和对比温度 Tr=T/Tc 为自变量查取 Z。对于含酸性气体的天然气，需采用 Wichert-Aziz 修正调整视临界参数后再查图。更精确的场合可使用 PR（Peng-Robinson）或 SRK 状态方程，通过迭代求解 Z 的三次方程根。下面给出一个可复现的示例计算思路：

已知：P = 6.0 MPa，T = 313 K（40℃）
混合气视临界参数（示例）：Pc = 4.62 MPa，Tc = 198 K
对比参数：Pr = 6.0 / 4.62 ≈ 1.30
          Tr = 313 / 198 ≈ 1.58
查 Standing-Katz 图或拟合式 → Z ≈ 0.86（示例值）
结论：该工况下实际气体体积比理想气体小约 14%，
      质量流量相同时，实际体积流量与压缩功均需按 Z 修正。

需要强调：Z 不是常数，它随压力和温度变化。压缩机从入口到出口，压力升高、温度升高，Z 沿压缩路径持续变化。因此能耗计算中不能用一个入口 Z 值"一算了之"，而应在每一级进出口分别取 Z，或沿压缩过程积分。

三、热值：把能耗换算成"单位气量成本"

热值分为高位热值（Hh，包含烟气中水蒸气凝结潜热）和低位热值（Hl，不包含）。压缩机站能耗分析中，热值的作用是把"压缩功"与"输送气量"联系起来：同样压缩 1 标准立方米天然气，热值越高，单位能量对应的气量越小，比能耗指标越好看。热值同样由组成加权计算，甲烷约 35.9 MJ/m³（高位），乙烷约 65.3 MJ/m³，丙烷约 93.6 MJ/m³，因此重烃含量上升会显著抬高热值，也会改变压缩因子。

实际工程中，热值数据常来自在线色谱分析仪，需注意其给出的基准（0℃或 20℃、干基或湿基）与计算所用标准状态一致，否则会出现百分之几的系统偏差。

四、压缩过程热力学：多变指数与功率估算

压缩机站能耗分析最核心的公式是多变压缩功。理想等温压缩功最小，绝热压缩功最大，实际过程介于两者之间，用多变指数 n 描述。对于天然气，n 通常在 1.25～1.35 之间，取决于冷却效果和气体性质。单级多变压缩功可写为：

W = (n/(n−1)) · Z · R · T₁ · [(P₂/P₁)^((n−1)/n) − 1]
其中 R 为气体常数（按视相对分子质量换算），T₁ 为入口绝对温度，
Z 取该级进出口压缩因子的平均值。

以示例工况计算：入口 P₁=2.0 MPa、T₁=303 K，出口 P₂=6.0 MPa，n=1.30，Z 平均取 0.90，视相对分子质量 17.5（R≈475 J/(kg·K)），则单位质量压缩功约为 1.30/0.30 × 0.90 × 475 × 303 × [(3.0)^(0.2308) − 1] ≈ 1.30/0.30 × 0.90 × 475 × 303 × 0.286 ≈ 1.60×10⁵ J/kg 量级。再乘以质量流量并除以压缩机效率（通常 0.80～0.88），即得轴功率。这一串计算中，Z、n、R 三个参数都来自本节建立的物性体系，任何一个取值不当都会直接放大到功率结果上。

五、把物性数据组织成可复用的计算表

为了让后续流程计算与能耗分析可追溯，建议为每个计算工况建立一张物性参数表，至少包含：组成（摩尔分数）、视相对分子质量、视临界压力与温度、工况压力与温度、对比参数、压缩因子、热值、多变指数及取值依据。下表为示例结构（数值为教学假设）：

这张表的价值在于：当气质发生变化或工况调整时，只需更新组成与工况两列，其余参数按既定方法重算，能耗结果即可快速复现与对比。物性数据体系不是一次性工作，而是贯穿流程设计与能耗分析全过程的"基准尺"。