› 实在气体压缩因子 Z 模拟器 返回
热力学

实在气体压缩因子 Z 模拟器

理想气体公式 PV=nRT 虽然好用，但在高压或低温下与实在气体的偏差会很大。本工具使用 Lee-Kesler/Pitzer 对应状态相关，根据压力、温度、临界常数和偏心因子，实时计算压缩因子 Z。可一眼看出与理想气体的误差，培养压缩机、管道和储存罐设计的直觉。

参数设置
压力 P
MPa
系统的绝对压力。用于计算 Pitzer 的 P_r = P / P_c
温度 T
K
系统的绝对温度。用于计算换算温度 T_r = T / T_c
临界压力 Pc
MPa
物质特有的临界压力。甲烷 4.6、N₂ 3.40、CO₂ 7.38、水 22.06 MPa
临界温度 Tc
K
物质特有的临界温度。甲烷 190.6、N₂ 126.2、CO₂ 304.1、水 647.1 K
偏心因子 ω
分子形状修正项。Ar 0.000、CH₄ 0.011、C₂H₆ 0.099、CO₂ 0.225、H₂O 0.344
计算结果
—
换算温度 Tr
—
换算压力 Pr
—
Z⁽⁰⁾ 单纯流体项
—
Z⁽¹⁾ 偏心补正项
—
压缩因子 Z
—
实摩尔体积 (cm³/mol)
—
分子运动 — 引力与斥力的支配区域

在中等压力下 Z<1 时，分子间引力占优，分子会有所聚集。在超高压下 Z>1 时，分子大小的斥力占优，分子相互排斥。颜色用绿色（≈1）→橙色→红色表示当前的 Z。

Z vs Pr — 换算压力对压缩因子的影响
Z vs Tr — 换算温度对压缩因子的影响
理论与主要公式

$$Z=\frac{PV}{nRT}=Z^{(0)}+\omega\,Z^{(1)}$$

Z=1 为理想气体。Z<1 为饱和线附近的中压区域（分子间引力主导），Z>1 为超高压区域（分子大小的斥力主导）。ω 是分子形状修正的偏心因子。

$$Z^{(0)}=1+\left(0.083-\frac{0.422}{T_r^{1.6}}\right)\frac{P_r}{T_r},\quad Z^{(1)}=\left(0.139-\frac{0.172}{T_r^{4.2}}\right)\frac{P_r}{T_r}$$

Vogel-Schlosser 对 Pitzer 二项相关的简化形式。使用换算坐标 T_r=T/T_c、P_r=P/P_c。在低～中压（P_r<2）的单相气体中精度良好。

$$V_{\text{ideal}}=\frac{RT}{P},\qquad V_{\text{real}}=Z\cdot V_{\text{ideal}}$$

理想摩尔体积与实摩尔体积（R=8.314 J/(mol·K)）。Z 正是表示"实气体的体积是理想气体的多少倍"。

压缩因子 Z 的基本概念

🙋
所谓的"实在气体"，就是高中学到的 PV=nRT 不成立的气体吧？但为什么要用一个"Z"这样的数字来表示呢？
🎓
很好的问题。实际上 PV=nRT 是基于"气体分子大小为零""分子间无作用力"这样的假设推导出来的。在常温、低压下，这个假设的偏差不到 0.5%，但压力超过数 MPa，或者接近临界点时，偏差就会变大。压缩因子 Z = PV/(nRT) 就是把这个偏差用一个数字来表示。Z=1 表示理想气体，Z≠1 表示"相对于理想气体偏离了多少"。这样用一个数字就能表示，在设计和教科书中，首先求出这个 Z 是业界的标准做法。
🙋
那么，Z 通常是大于 1 还是小于 1 呢？
🎓
有趣的是，两种情况都可能出现。在中等压力（数 MPa～十几 MPa）下，分子之间有"轻微相吸"的效果占优，Z<1，通常在 0.6～0.95 之间。如果继续升高压力，分子被挤在一起，分子的体积效应开始显现，Z>1。极端的例子，水氢气在 100 MPa 时，Z 可以上升到 1.7 左右。按照默认值（甲烷 5 MPa、350 K），Z ≈ 0.955，也就是实气体比理想气体体积小 4.5%。
🙋
对每种物质都要画一张图表会很麻烦吧，你们怎么计算的？
🎓
这里引入了"对应状态原理"和"Pitzer 偏心因子"的概念。如果用 T_c 和 P_c 来归一化温度和压力，得到换算坐标 (T_r, P_r)，那么许多物质的 Z 曲线竟然能重合在一起！这就是对应状态原理。但如果分子不是完全球形，会有轻微的偏差，用补正参数 ω（偏心因子）来表示，Z = Z⁽⁰⁾ + ω·Z⁽¹⁾，这样从甲烷到 n-辛烷都能用同一个式子处理。本工具使用的是这个式子的简化版（Vogel-Schlosser），在低～中压时精度在 2% 以内。
🙋
在实际应用中，Z 算错了会导致什么样的失败呢？能举个具体的例子吗？
🎓
最常见的是天然气管道设计。天然气主要是甲烷，在 5 MPa 下的实际流量用理想气体公式算的话，会高估 4～5%。有人甚至因此少选了一级压缩机，结果吃了大亏。还有更极端的例子，CCS（CO₂ 固碳）中，CO₂ 在 10 MPa、常温下输送时，Z 会降到 0.4～0.6，用理想气体公式的话容器体积会算得大 2 倍左右。在冷冻机的 R-410A 或氢气加氢站的 70 MPa 氢气罐容量计算中，Z 修正也是必不可少的。"选压缩机时最早要出现的数字"，可以这么说。
🙋
那么，Z=1 的假设在什么范围内是可以接受的？
🎓
一个经验法则是 T_r > 2 且 P_r < 0.5。比如空气在常温、常压（T_r ≈ 2.3、P_r ≈ 0.027）时，Z = 0.9997，基本上可以认为是 1。反过来，如果 T_r < 1.5 或 P_r > 1，就必须把 Z 修正考虑进去。进一步，如果 T_r < 1，就进入了饱和蒸汽的区域，可能出现两相，超出了这个简化 Pitzer 相关的适用范围。处理临界点附近或液相时，要切换到 SRK、Peng-Robinson、Lee-Kesler 完整相关这样的状态方程。

常见问题

压缩因子 Z 是表示实在气体与理想气体偏差程度的无量纲数，定义为 Z = PV/(nRT)。对于理想气体，Z=1；但实在气体具有分子大小和分子间引力，因此 Z≠1。在中等压力下，分子间引力为主导，Z<1（约 0.6～0.95）；在超高压下，分子大小的排斥力为主导，Z>1。在压缩机选型、管道设计和冷冻循环计算中，使用理想气体近似会产生数十%的误差，因此 Z 的评估在实际应用中极其重要。
Pitzer 的二项相关是在对应状态原理中引入偏心因子 ω 的修正，形式为 Z = Z⁽⁰⁾(T_r, P_r) + ω·Z⁽¹⁾(T_r, P_r)。本工具采用 Vogel-Schlosser 近似 Z⁽⁰⁾ = 1 + (0.083 − 0.422/T_r^1.6)·P_r/T_r 和 Z⁽¹⁾ = (0.139 − 0.172/T_r^4.2)·P_r/T_r，在低～中压（约 P_r < 2 且 T_r > 0.8）的单相气相中精度可达 2% 以内。在临界点附近、两相区域和液相中误差会增大，需要使用 SRK/PR/Lee-Kesler 完整相关等状态方程。
偏心因子 ω 是 Pitzer 引入的分子形状修正参数，从 T_r = 0.7 时的换算饱和蒸汽压 P_r,sat 定义为 ω = -log10(P_r,sat) - 1。完全球形、无极性的单纯流体（Ar、Kr、Xe）的 ω ≈ 0；甲烷为 0.011、乙烷为 0.099、n-庚烷为 0.349、水为 0.344。标准化学工程手册（Reid-Prausnitz-Poling 等）中列有许多物质的数值。本工具的默认值 0.011 对应甲烷。
误差取决于具体条件。在低压、高温条件下（T_r > 2、P_r < 0.5），误差小于 1%，在实际应用中无问题。但在常温 350 K、5 MPa 下输送天然气（甲烷主体）时，Z ≈ 0.955，误差约 5%；在 10 MPa 下输送二氧化碳用于 CCS 管道时，Z ≈ 0.4～0.6，误差达 40～60%。这直接影响压缩机段数选择和储存罐尺寸，因此在超临界或高压区域必须应用 Z 修正。本工具右侧卡片显示"理想摩尔体积 vs 实摩尔体积"，其比值 = Z，可一目了然地看出误差规模。

实际应用场景

天然气管道与压缩机设计：甲烷主体的天然气在 5～10 MPa 下长距离输送时，Z 在 0.85～0.95 之间。不考虑 Z 的话，气体容量、流速和摩擦损失的计算会有数%的误差，导致压缩机段数选择失误或达不到合同流量。北美和欧洲的长距离管道中，AGA-8 或 GERG-2008 状态方程是标准，本工具的 Pitzer 相关相当于这些的简化预检版本。

CO₂ 回收与地质储存（CCS）和超临界输送：从火电厂或钢铁厂回收的 CO₂ 液化或超临界态地下储存，通常在 7～15 MPa、常温下输送。CO₂ 的 T_c=304 K、P_c=7.38 MPa，正好在临界点附近，Z 在 0.2～0.6 之间大幅变动。理想气体的话会把管道口径和泵的规格算错 1 倍左右，Z 修正是 CCS 设计的第一步。

高压氢气储存罐与加氢站：燃料电池汽车（FCV）用的高压氢气在 35 MPa 或 70 MPa 下储存。氢气的 T_c=33 K 很低，常温下 T_r 很大，Z>1，在 70 MPa 时 Z ≈ 1.4～1.5。也就是说，同样体积的罐子里，实际能装的氢气比理想气体少 30～40%。FCV 航程计算如果用 Z=1，会过高估计，对加氢站设计和车辆开发造成致命错误。

冷冻循环与制冷剂选择：HFC（R-134a、R-410A）和新型 HFO 制冷剂（R-1234yf）在蒸发器、冷凝器的工作温度范围内，T_r 在 0.7～0.95，Z 在 0.7～0.9 的区域。COP（性能系数）和冷凝器尺寸的计算，需要用 Pitzer 或 PR/SRK 状态方程来评估 Z，这是空调、热泵效率设计的核心。

常见误解与注意事项

最大的陷阱是"常温常压下 Z=1 就行"这个观念被带到高压区域。用本工具的滑块试试，把甲烷压到 10 MPa 时，Z 会跌到 0.85 以下。理想气体的话管道流量会高估 15%，可能导致压缩机能力不足。经验法则是"P_r 超过 0.3 就要加上 Z 修正，P_r 超过 1 就一定要加"。常温常压的 Z ≈ 1 只是在很低压下才成立，设备设计的前提条件中不能包含这个假设。

其次，"Pitzer 相关并非万能"。本工具的 Vogel-Schlosser 简化 Pitzer 在 T_r > 0.8、P_r < 2 的单相气体中精度在 2% 以内，但在临界点附近（T_r ≈ 1）时，分子间关联会以爆炸性方式变动，此时需要用 Lee-Kesler 完整相关或 Peng-Robinson 状态方程。此外两相区域（饱和蒸汽和饱和液共存）中，Z 本身有两个值，本工具的公式不适用。极性较强的物质如水（ω=0.344）或氨，也应该用更高级的相关（SRK、PRSV、SAFT）。

最后，"Z 不仅是气体的参数"。Z 也可以对液相定义，液相中 Z 通常很小，在 0.01～0.3 之间（因为摩尔体积比气相小 2～3 个数量级）。本工具针对气相，但在化学工程中，通常用同一个状态方程（PR 或 SRK）同时计算气相和液相的 Z，从两者的差异推算饱和蒸汽压和气液平衡。Z 的概念本身是一个强大的工具，能统一处理气相、液相和超临界相。

使用指南

  1. 输入压力（0.1～1000 bar）、温度（250～500 K）以及物质的临界压力、临界温度和偏心因子ω。可以选择 CO₂（Pc=73.8 bar、Tc=304.13 K、ω=0.228）或 CH₄（Pc=46.0 bar、Tc=190.6 K、ω=0.011）等实在气体的参数。
  2. 模拟器自动计算换算温度 Tr=T/Tc 和换算压力 Pr=P/Pc，使用 Lee-Kesler 相关公式计算 Z⁽⁰⁾（单纯流体项）和 Z⁽¹⁾（偏心补正项）。
  3. 得到压缩因子 Z=Z⁽⁰⁾+ωZ⁽¹⁾，实摩尔体积 V=ZRT/P 以 cm³/mol 为单位实时显示。可直观看出与理想气体（Z=1）的偏差。

具体计算示例

氮气（N₂：Pc=33.9 bar、Tc=126.2 K、ω=0.039）在压力 50 bar、温度 300 K 下压缩的情况：Tr=2.38、Pr=1.47 为换算坐标，经过 Z⁽⁰⁾≈0.920、Z⁽¹⁾≈0.015，得到 Z≈0.926。实摩尔体积相对于理想值 480 cm³/mol 会压缩到约 445 cm³/mol，可用于工业气体精馏装置的流量设计。

实务中的注意点