对于超临界进料气体，汽化因子 (Kᵢ) 无法从蒸气压表中查得——它是使用非对称热力学框架确定的。控制关系式为 Kᵢ = (Hᵢ · γᵢ) / (φᵢ · P)*，其中 Hᵢ 是气体在溶剂中的亨利常数，γᵢ* 采用非对称规定（当 xᵢ → 0 时，γᵢ* → 1），φᵢ 是气相逸度系数，P 是总压。该表达式取代了传统的拉乌尔定律 K 值，因为纯气体蒸气压在临界温度以上是没有定义的。
确定氢气或氧气等超临界进料的 Kᵢ 完全取决于带有非对称活度系数的亨利定律。Kᵢ = yᵢ/xᵢ 的值变为 (Hᵢ · γᵢ*) / (φᵢ · P)，将溶解度、分子相互作用和气相非理想性联系在一起——这是将溶解气体视为永远不会冷凝成纯液体的溶质的直接结果。
超临界气体的热力学基础
为何传统的气液平衡会失效
经典的 K 因子，yᵢ/xᵢ = Pᵢˢᵃᵗ/P，依赖于纯组分蒸气压。
对于温度超过其临界温度的组分（例如 298 K 下的氢气），Pᵢˢᵃᵗ 不存在。气体无法仅通过压力液化，因此标准对称规定变得毫无意义。
因此，设计中试规模气液反应器的工程师必须转向一个将超临界组分视为溶解溶质而非可冷凝蒸气的框架。平衡条件 fᵢⱽ = fᵢᴸ 仍然成立，但液相逸度必须参考亨利定律而不是纯液体逸度。
非对称规定与亨利定律
超临界溶质的液相逸度表示为 fᵢᴸ = γᵢ · Hᵢ · xᵢ。*
此处 Hᵢ 是亨利常数（溶质-溶剂对在系统温度下的性质），γᵢ* 是一个活度系数，经归一化处理使得当摩尔分数 xᵢ 趋近于零时 γᵢ* → 1。这种非对称规定将标准态锚定在无限稀释处，这正是亨利定律适用的区域。
Hᵢ 捕捉了极限溶解行为——气体与液体之间的相互作用越强，Hᵢ 越小，溶解浓度越大。当 xᵢ 并非极小时，γᵢ* 随后会对偏离理想稀释行为的情况进行修正。
Kᵢ 的完整表达式
使气相和液相逸度相等，得到 φᵢ · yᵢ · P = γᵢ · Hᵢ · xᵢ。*
重新排列可直接得出汽化因子：
[ K_i = \frac{y_i}{x_i} = \frac{H_i \cdot \gamma_i^*}{\phi_i \cdot P} ]
每一项都是可获取的：Hᵢ 来自文献关联式或有限的一组中试溶解度测量数据；γᵢ* 来自针对非对称规定调整的过量吉布斯能模型（例如 NRTL、UNIQUAC）；φᵢ 来自立方型状态方程（Peng-Robinson、Soave-Redlich-Kwong）；P 仅仅是反应器操作压力。每当进料流在其临界温度以上进入时，该表达式是拉乌尔定律类比法的逻辑且必要的替代方案。
在中试装置环境中应用该方程
获取亨利常数
Hᵢ 是计算的基石。
它对温度有很强的依赖性，并且在中等压力到高压下，必须通过 Krichevsky-Kasarnovsky 方程 [ \ln H_i(P) = \ln H_i(P_{\text{sat}}) + \frac{\bar{v}i^\infty}{RT}(P - P{\text{sat}}) ] 对压力进行修正，其中 (\bar{v}_i^\infty) 是溶质在无限稀释下的偏摩尔体积。
对于常见有机溶剂中的许多中试气体（H₂、O₂、CO），可靠的 Hᵢ 值已列成表格，或可以通过基团贡献法预测。在缺乏数据的情况下，几次台架规模的溶解度实验可以迅速建立反应器温度条件下的亨利常数。
估算活度系数 (γᵢ*)
即使在相对稀的溶液中，如果溶剂复杂或存在多种溶质相互作用，γᵢ 也可能偏离 1。*
当超临界气体的液相摩尔分数保持在 0.01 以下时，假设 γᵢ* = 1 通常引入的误差小于 5%——这对于早期中试装置的尺寸估算是可以接受的。
对于浓度较高的系统，UNIQUAC 或 NRTL 等模型可以针对非对称规定重新参数化。必要的二元相互作用参数通常通过气体溶解度数据回归得出，或通过分子相似性估算。因此，中试装置工程师的决策树分为两支：为了快速范围界定使用 γᵢ* = 1，或者当准确的质量传递预测至关重要时，投入精力使用完整的活度模型。
计算气相逸度系数 (φᵢ)
在许多中试反应器常见的中等压力（5–20 bar）下，φᵢ 与 1 的偏差很少超过几个百分点。
然而，高压加氢或氧化反应通常在 50–200 bar 下运行，此时气相非理想性变得显著。与气体混合物和温度相匹配的立方型状态方程可以简便地提供 φᵢ。
例如，在 Peng-Robinson 框架中，φᵢ 是对比温度、压力和偏心因子的函数。对于纯氢气，在 100 bar 下，φᵢ 修正可使 Kᵢ 比理想气体假设增加 10–30%，因此忽略它会导致低估气相摩尔分数并扭曲质量传递推动力的计算。
压力的作用
压力显式地出现在 Kᵢ 表达式的分母中，使得当所有其他项保持不变时，Kᵢ 与 P 成反比。
这意味着随着反应器压力升高，给定液相浓度下的平衡气相摩尔分数会下降——这直接反映了亨利定律：更多的气体溶解，因此保留在顶部空间或气泡相中的气体减少。
在中试装置设计中，这种压力效应必须与经常被忽视的 Hᵢ 本身的压力依赖性相结合。在高压下，Hᵢ 会增加，部分抵消 P 的影响。因此，完整的 Kᵢ 模型将 Krichevsky-Kasarnovsky 修正与分母耦合，确保汽化因子准确反映实际操作条件下的溶解度增强。
理解权衡与局限性
稀溶液假设
假设 γᵢ* = 1 很诱人，但并不总是安全的。在高浓度的加氢混合物中，或当溶解气体分数超过几个摩尔百分比时，活度系数可能偏离 1 达 10–20%。 由此产生的 Kᵢ 误差会传播到液相浓度估算中，并最终影响反应器体积预测。对于早期可行性研究，这种误差通常是可以容忍的；但对于最终放大，则不行。
温度和压力敏感性
轻质气体的亨利常数通常在 50 K 的温度变化下翻倍或三倍，且逸度系数随压力变化。如果反应器的温度分布发生变化——这在启动或催化剂失活期间很常见，那么在一个操作点有效的 Kᵢ 值可能会变得不准确。中试团队在对反应器在整个操作包络内进行建模时，必须将 Kᵢ 视为动态变量，而不是静态数字。
多组分混合物
当气体进料包含多种超临界组分（例如合成气中的 H₂ 和 CO）时，逸度系数 φᵢ 变为混合物属性，状态方程中的交叉相互作用参数就变得重要。使用纯组分 φᵢ 估算值可能会引入系统性偏差。 只要有可能，中试装置工程师应根据少量的总压和气相组成测量来验证 Kᵢ 预测，并根据需要调整状态方程混合规则。
为中试装置目标做出正确选择
确定 Kᵢ 时适当的详细程度完全取决于您需要该数字用于何种目的。
- 如果您的主要关注点是快速可行性筛选： 使用简化形式 Kᵢ ≈ Hᵢ/P，配合文献 Hᵢ 值，并假设 γᵢ* = φᵢ = 1。这可以在一个下午的时间内给出气液分配和最小溶剂流率的数量级估算。
- 如果您的主要关注点是准确的质量传递建模： 采用完整表达式，使用经压力修正的 Hᵢ（通过 Krichevsky-Kasarnovsky 关系式），用于 φᵢ 的最先进状态方程，以及至少一阶活度模型用于 γᵢ*。您的鼓泡塔或搅拌釜模拟将提供经得起放大审查的 KLa 值。
- 如果您的主要关注点是高压操作（30 bar 以上）： 绝不要忽略逸度系数修正。立方型状态方程——拟合至已发表的二元数据或测量的露点——是避免系统性低估气相浓度的最低要求。
- 如果您的主要关注点是溶剂筛选或多组分进料： 投资进行实验性溶解度运行，以回归 Hᵢ 和 γᵢ* 参数。这项努力通过大幅减少反应器尺寸估算的不确定性并避免下游分离装置的意外情况，物有所值。
准确确定超临界进料的 Kᵢ 并非教科书式的练习——它是将溶解度、质量传递和反应器尺寸连接成可靠中试设计的关键环节。
汇总表：
| 参数 | 符号 | 在 $K_i$ 方程中的作用 | 关键估算方法 |
|---|---|---|---|
| 汽化因子 | $K_i$ | 衡量气液平衡 ($y_i/x_i$) | 计算为 $(H_i \cdot \gamma_i^*) / (\phi_i \cdot P)$ |
| 亨利常数 | $H_i$ | 定义极限溶解行为 | Krichevsky-Kasarnovsky 方程修正 |
| 活度系数 | $\gamma_i^*$ | 修正液相非理想性 | NRTL / UNIQUAC（非对称规定） |
| 逸度系数 | $\phi_i$ | 修正气相非理想性 | 立方型状态方程（例如 Peng-Robinson） |
| 系统压力 | $P$ | 代表反应器操作压力 | 显式分母；与 $K_i$ 成反比 |
优化气液反应器热力学模型需要精确的中试规模测试。LABPARK 为大学、科研院所和企业提供优质的化工、生物工艺与生物技术以及环境与水处理领域的教育和职业培训单元操作中试装置。通过可靠、可放大的系统赋能您的研究和培训——立即联系 LABPARK 讨论您的项目！
