天然气液化流程模拟中立方型状态方程的计算
天然气引射式液化流程的计算和分析
初 始 运行 时 ，原 料气 冷 却所 需冷 量 由外加 制冷 剂 提 供 。 当 第 三 换 热 器 出 口 的 天 然 气 温 度 达 到 ２０３．１５ Ｋ 时 ，天然 气 从拉 发尔 喷 管喷 出并进 入 引射 器 ，由于 没有被 引 射气 体 ，故此刻 仅 以天 然气 组 分 从 引 射 器 出 口 流 出 ， 控 制 其 出 口 温 度 使 之 在 １．２ ＭＰａ压 力下 的气 化率 为 ０．６７。本 文选 用 ＳＲＫ 方 程 计算 天 然气 气液 相 的相平 衡 ，通 过 比较 计算 结 果 和实验结果，可得出用 ＳＲＫ 方程计算烃类混合物 的气液 相 平衡 是可 靠 的 。利用 ＡＳＰＥＮ．ＨＹＳＹＳ计算 得 前 述 天 然 气 组 成 条 件 下 ， 引 射 器 出 口温 度 为 １５２．４２６ ３ Ｋ，其 气液 相 组成 见表 ２：
１ 液 化流 程描 述
天然气引射式液化流程 的引射气体 由天然气 、 回流 复热气 两 种气 体组 成 。１ ｋｍｏｌ／ｈ、４．０ ＭＰａ、３００Ｋ 的原料 天 然气经 加 压 、水冷 后 ，压力 为 ２０ ＭＰａ，温 度 为 ３１０ Ｋ。回流气 为 ０．６７ ｋｍｏｌ／ｈ，１．２ ＭＰａ，经 多 次 引射 、 回流 、复 热后 ，温 度为 ３０７ Ｋ，然 后经 加 压 、水 冷 至压 力为 ２０ ＭＰａ，温 度 为 ３１Ｏ Ｋ。天然气 和 回流 复热 气 以相 同压 力 ，相 同温度 汇 合进 入 高压 引射 气体 管道 ，此 时总摩 尔 流量 为 １．６７ ｋｍｏｌ／ｈ。汇 合气 经第 一 换热器 被 冷却 至 ２８６．２２ Ｋ，接 着进 入 外 加制 冷剂 的第 二换 热 器 ，被 冷 却至 ２１８．２７ Ｋ，然后 进入 第三 换热 器 ，最 终被冷 却 至 ２０３．１５ Ｋ，进 入 高 压拉发尔喷管 ，以高压、低温 的超 临界状态进入拉 发尔喷管 ，之后以超临界速度喷入引射器，并从中 间罐 引射 出气相 组 分 ，在 混合 室 中组成 引射 及被 引 射气 体 的 混合气 体 ，再经 引射 器 的扩 压管 减速 ，恢 复静 压 至 １．２ １Ｖ［Ｐａ。此 过 程 中 ， 由于 急骤 膨 胀促 使 混 合气 体 出现 降温 ，并在 临界 状 态下 出现 两相 流 。 在分 离器 １中分 出 的气相 即 为 回流 、复热气 体 ，此
《石油与天然气化工》2007年1～6期题目索引
（ ：７ ４） ２ ５
其 甲烷 燃 烧 活 性 的研 究 … … … … … … （ ）：１ １ １ 秋 林 轻 烃 回 收 装 置 冰 堵 分 析 及 改 进 措 施 … （ ）：４ １ １
天 然 气 中含 蜡 成 分 对 处 理 装 置 运 行 的 影 响 分 析 及对 策 … …… ……… ……… ……… （ ）２２ ４ ：８
（ ６）：５ ４２
（ ）：８ ３ １５
大 涝 坝 站 导 热 油 系 统 运 行 工 况 分 析 及 防护 对 策
… … … … … … … … － … … … … … … 一
基 于 Ｃ ｅ ｓｐ ｈ ｍ ｅ Ｌ的 甲 醇预 精 馏 塔 模 拟 与 分 析
… ． … … ・ … … … … … ・ … … … ・ … … － －
（ ５）：６ ３６
（ ）：７ ２ ９
气体水合 物 结晶驱 动力定 义辨析 Ｌ Ｇ 管路 保 冷 厚 度 的 计 算 Ｎ
Ｌ Ｆ Ａ１： ａ ｅ
一
… … … （ ） ３０ ５ ：７
复 合 掺 杂钙 钛 矿 氧 化 物 催 化 剂 的 制 备 方 法 …
… … … … … … … … … … … … … … …
ＨＹ 分 子 筛 合 成 Ｅ Ｂ 的 研 究 … … … … … （ ）：５ ＴＥ ５ ３３
（ ） ２ １ ：２
甲 醇 乙 醇一 步 合 成 异 丁 醛 反 应 用 Ｖ Ｏ 。 催 化 剂
的改性研 究 ……… …… ……… …… （ ５）：５ ３６
反 相 微 乳 法 制 备 二 甲 醚 燃 烧 催 化 剂 Ｌ ＭＡ１ Ｏ．一 Ａ ． 。 ．
天然气的计算方法
天然气计算方法体积计算•标准立方米 (Sm³)：天然气在 15°C（59°F）和 1 个大气压 (101325 Pa) 下的体积。
•立方英尺 (ft³)：天然气在指定温度和压力下的体积。
•百万标准立方英尺 (MMscf)：100 万标准立方英尺。
能量计算•百万英热单位 (MMBtu)：天然气中包含的能量，等于 100 万英热单位。
•标准立方英尺当量 (Mcfe)：天然气的体积，其能量含量等于 1000 立方英尺甲烷的能量含量。
换算公式•Sm³转换为 MMBtu：MMBtu = Sm³ x 热值（MJ/m³）/ 1055.06•ft³转换为 Sm³：Sm³ = ft³ x 温度系数 / 压力系数•Mcfe 转换为 MMBtu：MMBtu = Mcfe x 1温度和压力系数•温度系数：将天然气从标准条件转换为特定温度时的因子。
•压力系数：将天然气从标准条件转换为特定压力时的因子。
温度和压力系数可以通过查阅燃气协会或其他权威机构提供的表格来获得。
热值天然气的热值因其成分而异，通常以百万焦耳每立方米 (MJ/m³) 为单位。
热值通常由气体供应商提供，或可以通过实验室分析确定。
示例计算•将 100 Sm³天然气转换为 MMBtu：热值为 38 MJ/m³，因此 MMBtu = 100 x38 / 1055.06 = 3.6 MMBtu•将 1,000 ft³天然气在 60°F 和 5 psia 下转换为 Sm³：温度系数为 0.964，压力系数为 1.075，因此 Sm³ = 1,000 x 0.964 / 1.075 = 901.2 Sm³。
天然气液化流程模拟中立方型状态方程的计算
天然气 的液 化 … 流 程 主 要 分 为 三 种 ： 阶式 （ 级 联式 ） 液化 流程 、 丙烷 预 冷混 合 制 冷剂 液 化流 程 、 带 膨胀 机 的液 化流程 。 目前应用 最广泛 的是 丙烷 预冷 混合 制冷 剂 液 化 流 程 。在 各 种 流程 中 降 低 液 化 天然 气 的生 产成本 ， 都是 重 点要 解 决 生 产过 程 中 的 节 能降耗 。对于 天然 气 液 化 的整 个 过 程 系统 来 说 ， 不 能够 只是 考虑单 个 设 备 的节 能 ， 应 该考 虑 整 个 而 过程 系统节 能 ， 因此要 使 天然 气 液 化工 艺 设 计 的 能耗最 小 、 费用最小 ， 就必 须将整 个系 统集成 起来作 为一个 有机结 合 的整 体来 看 待 ， 到 整体 设 计 最优 达 化 。由 此 ， 要 对 液 化 流 程 进 行 模 拟 计 算 及 优 需 化 。以下便针对 天 然气 液 化循 环 流程 中最 常用
为 后 续 计 算 焓 、 等 参 数 提 供 相 应 的解 决 办 法 ， 判 断 选择 的 状 态 方 程 是 否 符 合要 求 。 熵 并
关键 词 天 然气
液化 流程
状 态方程
相 平衡ห้องสมุดไป่ตู้ 算
随着我 国不断 引进液化 天然气 以及大力 发展 国 内液化 天然气 产业 ， 化 天 然气 在 我 国能 源 结 构 中 液 占据越来 越重 要 的地 位 。因 此 ， 化天 然 气 技术 的 液 发展显得 至关 重要 。在 国外 ， 化 天然 气 技 术 的研 液 究 非常成 熟 ， 是 由于知识产 权及 经济 利益 的限制 ， 但
摘 要 天 然 气 液 化 流 程 广 泛 采 用 的 是 丙 烷 预 冷 混 合 制 冷 剂 液 化 流 程 。 为 了 进 一 步 优 化 流
天然气液化中气液相平衡计算研究详解
态方程
杂，参数不广谱，在大分子和缔合性流体的描叙发
挥越来越大作用。
混合法则
1
维里型
混合 法则
2
范德华型
3
与组成相关的交互作
用参数的混合法则
范德华单流体混合法则
bm
xi x jbij am
xi x j aij
ij
ij
bij
(bi
bj ) 2
am
aiaj (1 kij )
综合选择
• ①状态方程选择： • SRK方程：形式简单，改善了气液相逸度计
平衡常数
SRK
算 法 框 图
编程思路与方法
• 1.方程简化： • 2.内嵌算法：压缩因子方程（试算法） • 3.循环节控制：迭代精度，初值选取
1.方程简化
ln i
bi b
(Z-1) ln(Z B)
A B
2 a
n j 1
z
j
(ai
a
j
)0.5
(1
ki,
j
)
bi b
ln
1
B Z
• double S[m_size],X1[m_size],AA[m_size]; • for (int i=0;i<m_size;++i) • double XK=0.0; • S[i]=(bi[i]/b)*(m_Zl-1)-log(m_Zl-B); • for(int j=0;j<m_size;++j) • XK+=pow(ai[i]*ai[j],0.5)*Zi[j];
结论
• （1）通过绪论和各种方程，模型的对比研究，我 们可以看出，SRK方程在不失立方型方程形式简 单的情况下，大大的改善了计算气液相平衡的效 果，并且对方程所需的基础数据少，能够很好的 获得纯组分及混合物的气液相组成摩尔分数，气 化率等。
天然气液化流程模拟及其工艺计算
4 NACE Standard TM0284 - 96. Evaluation of pipeline and pressure vessel steels for resistance to hydrogen induced cracking.
·98 ·
图 2 透平示意图
对重烃分离器及节流的计算 , 可参见文献〔4〕。 透平膨胀机的计算〔3〕,透平示意图如图 2 所示 ,等熵 效率取为 0. 75 ,计算方法如下 。第一步 :由透平膨胀 机入口参数 pi1 、Ti1 、计算入口焓 Hi1 和入口熵 S i1 ;
© 1995-2005 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved.
图 1 工艺流程图
L —K 方程的形式为 :
Z
=
prV r Tr
=
1
+
B Vr
+
C
V
2 r
+
D
V
5 r
cr
9+7
c4
T
3r V
2 r
β+
γ
V
2 r
exp
-
γ
V
2 r
uc
式中 : B
= b1 -
b2/ T r -
b3/
T
2 r
-
b4/
T
3 r
C = c1 -
立方型方程,常熟状态方程,理论型状态方程额的
立方型方程,常熟状态方程,理论型状态方程额的立方型方程，即该方程中必须有一个三次多项式（aX³ +bX² +cX + d=0），其中a≠0。
这种方程的解可以由Viete公式求得。
Viete公式表示：如果有aX³ +bX² + cX + d=0，用p,q,r代表它的特征多项式的根（即X³+pX²+qX+r=0），则有：p = -a/3q = (3a c- bc²)/ (6a²)r = -(27a²d-9acb+2b³c)/(54a³)常熟状态方程是一种计算空气动力学物理量的抽象方程。
它与流体动力学的经典方程组相对应，实际问题的几个参数的变化本质上可以表示为这些方程的变形，使这些方程变成一种特定状态的常熟方程。
常熟状态方程可以与恒定温度的完全确定流体的经典动力学方程组等价转换，从而研究不同温度和压力状态下复杂流体的运动。
理论型状态方程是一种描述一般性热力学系统和可能出现在一般性物理系统中的多种各种状态变化的理论。
它以一种通用的、系统的方式来描述各种状态变化。
该方法可以定量地表达热力学系统的状态变化，并可以推断出影响因素。
它没有考虑系统中存在的具体物质，只考虑系统的扩散，聚合，分裂等基本变化的总的影响，以及系统本身的学习，应用和储存能力。
因此，理论型状态方程可以极大地概括和提高系统中状态变化的模拟仿真能力。
总之，立方型方程、常熟状态方程和理想状态方程都是用来研究一般性热力学系统和可能出现在一般性物理系统中的多种各种状态变化的常见方程类型，它们的用途各有不同，都能够提供有效的模拟仿真能力。
天然气液化流程
• l)能耗较高，比级联式液化流程高10%·20%左右;
• 2)混合制冷剂的合理配比较为困难;
• 3)流程计算须提供各组分可靠的平衡数据与物性参数，计 算困难。
• 带预冷的混合冷剂制冷循环，简称c3/MRC工艺， 是在MRc工艺基础上开发出来的新一代液化工艺， 也可以视其为对传统的阶式循环的改进。C3似 RC循环采用丙烷预冷(或者氨制冷预冷)与混合制 冷剂(NZ+Cl~c4)联合作用方式。
• 在混合制冷剂液化流程中，天然气被预冷换热器、液化换 热器、过冷换热器逐步冷却，最后经过闪蒸于储罐中储存。
与级联式液化流程相比，其优点是:
• 1)机组设备少，只需一台循环压缩机;流程简单;投资省， 投资费用比经典级联式液化流程约低15%一20%
• 2)管理方便;
• 3)混合制冷剂组分可以部分或全部从天然气本身提取与补 充。
• C3/MRC混合冷剂循环天然气冷却装置见图2一3。 • 图2一3是丙烷预冷混合制冷剂循环液化天然气流
程图。流程由三部分组成:1，混合制冷剂循环;2， 丙烷预冷循环;3，天然气液化回路。在此液化流 程中，丙烷预冷循环用于预冷混合制冷剂和天然 气，而混合制冷剂循环用于深冷和液化天然气。
• 混合制冷剂经两级压缩压缩至高压，首先用水冷却，带走一部分热量， 然后通过丙烷预冷循环预冷。预冷后进入气液分离器分离成气相和液 相，液相经第一换热器冷却后，节流、降温、降压，与返流的混合制 冷剂混合后，为第一个换热器提供冷量，冷却天然气和从分离器出来 的气相和液相两股混合制冷剂。气相制冷剂经第一换热器冷却后，进 入气液分离器分离成气相和液相，液相经第二个换热器冷却后节流、 降温、降压，与返流的混合制冷剂混合后为第二个换热器提供冷量， 冷却天然气和从分离器出来的气相和液相两股混合制冷剂。从第二换 热器出来的气相制冷剂，经第三换热器冷却后，节流、降温后进入第 三换热器，冷却天然气和气相混合制冷剂。
天然气液化流程模拟中立方型状态方程的计算
天然气液化流程模拟中立方型状态方程的计算
李佩铭;焦文玲
【期刊名称】《石油与天然气化工》
【年(卷),期】2007(036)006
【摘要】天然气液化流程广泛采用的是丙烷预冷混合制冷剂液化流程.为了进一步优化流程,减少能源的消耗,需要对整个流程进行模拟,而模拟过程中热力学参数的计算便是整个流程计算的基础.本文讨论使用两种立方型状态方程(SRK和PR方程)对热力学参数进行相平衡计算,为后续计算焓、熵等参数提供相应的解决办法,并判断选择的状态方程是否符合要求.
【总页数】4页(P444-447)
【作者】李佩铭;焦文玲
【作者单位】深圳市燃气集团股份有限公司;深圳市燃气集团股份有限公司
【正文语种】中文
【中图分类】TE6
【相关文献】
1.丙烷预冷混合制冷剂天然气液化流程模拟计算 [J], 沈立龙;赖建波
2.立方型状态方程的空气焦-汤效应数值模拟 [J], 程阳;李东升;胡佳成;余施美
3.立方型状态方程在 R134a 热力性质计算中的研究 [J], 刘晖;王昱凯
4.超临界水的热力学模拟:一种可用到4273K、2GPa的立方型状态方程 [J], 郭涛;胡家文;王向辉;靳丽花;王艳哲
5.天然气液化流程模拟及其工艺计算 [J], 刘新伟;李海国;刘芙蓉
因版权原因，仅展示原文概要，查看原文内容请购买。
液化天然气换算公式
液化天然气换算公式
液化天然气（LNG）是将天然气冷却至极低温下（约
162°C），使其转变成液态的过程。
液化天然气的体积较小，
便于储藏和运输。
在液化天然气的换算过程中，常用的单位是以体积为基础的，通常使用立方米（m³）或者标准立方英尺（scf）。
要进行液化天然气的换算，需要考虑两个主要因素：温度和
压力。
温度换算：
液化天然气的标准温度通常是162°C，也可以转换为开尔文（K）为111.15K。
压力换算：
液化天然气的压力通常会有不同的单位，常用的单位有帕斯
卡（Pa）和巴（ba）等。
常见的LNG压力可以使用常压
（1atm或1.01325×10^5Pa）或者标准压力（101.3kPa）
进行换算。
换算公式：
液化天然气的体积换算公式可以根据温度和压力进行计算。
1.从立方米（m³）到标准立方英尺（scf）的换算公式为：
1m³=35.3147scf
2.从标准立方英尺（scf）到立方米（m³）的换算公式为：
1scf=0.0283168m³
3.从帕斯卡（Pa）到巴（ba）的换算公式为：
1Pa=0.00001ba
4.从巴（ba）到帕斯卡（Pa）的换算公式为：
1ba=100000Pa
使用这些换算公式，可以根据给定的温度、压力和体积单位，在液化天然气的换算过程中进行计算。
立方型状态方程
立方型状态方程立方型状态方程是一种用于描述气体状态的数学模型。
它是基于理想气体状态方程和Van der Waals方程的推导而来，通过引入一个修正因子来更加准确地描述气体的行为。
下面将从立方型状态方程的定义、推导、应用以及优缺点等多个方面进行详细阐述。
一、定义立方型状态方程是指一种用于描述气体状态的数学公式，其形式为：P = (RT)/(V-b) - a/(V^2)其中，P表示气体压强，V表示气体体积，T表示气体温度，R为普适气体常量，a和b为修正因子。
这个公式也被称为德拜-胡克斯（DHB）方程。
二、推导立方型状态方程是基于理想气体状态方程和Van der Waals方程的推导而来。
理想气体状态方程可以表达为：PV=nRT其中n表示物质量。
但是，在高压和低温下，由于分子间相互作用力的存在，理想气体模型并不能完全适用于实际情况。
因此，在这些情况下需要引入修正因子。
Van der Waals 方程可以表示为：(P+a(n/V)^2)(V-nb)=nRT其中a和b为修正因子，分别表示分子间吸引力和排斥力。
然而，Van der Waals方程仍然存在一些缺陷，比如在高压下会出现负压等问题。
为了解决这些问题，科学家们提出了立方型状态方程。
该方程引入了一个修正项-a/(V^2)，用于更准确地描述气体的行为。
同时，它也可以通过改变修正因子a和b来适应不同的气体。
三、应用立方型状态方程在热力学和化学工程中有广泛的应用。
其中一些典型的应用包括：1. 气体储存和输送：立方型状态方程可以帮助工程师计算气体在不同温度、压力和体积下的密度，从而更好地设计气体储存和输送系统。
2. 化学反应：立方型状态方程可以作为化学反应过程中气态物质的基础模型，在计算反应热、平衡常数等参数时起到重要作用。
3. 石油工业：立方型状态方程可以帮助石油工业从原油中分离出不同组分，并计算其物理性质。
4. 环境保护：立方型状态方程可以用于计算大气污染物的扩散和传输，从而帮助环境保护部门制定有效的控制措施。
物理天然气体积计算公式
物理天然气体积计算公式天然气是一种重要的能源资源，广泛应用于工业生产、家庭供暖、发电等领域。
在天然气的生产、储存和运输过程中，准确计算天然气的体积是非常重要的。
本文将介绍物理天然气体积计算的公式及其应用。
物理天然气体积计算公式是通过理想气体状态方程来计算的。
理想气体状态方程描述了气体的状态与其压力、体积和温度之间的关系，其数学表达式为PV=nRT，其中P为气体的压力，V为气体的体积，n为气体的摩尔数，R为气体常数，T为气体的温度。
根据理想气体状态方程，可以推导出物理天然气体积计算的公式如下：V = (nRT)/P。
其中，V为天然气的体积，n为天然气的摩尔数，R为气体常数，T为天然气的温度，P为天然气的压力。
在使用物理天然气体积计算公式时，需要注意以下几点：1. 温度单位的转换，在计算天然气体积时，温度的单位通常使用开尔文（K），而常见的温度单位是摄氏度（℃）。
因此，在使用公式计算天然气体积时，需要将温度转换为开尔文单位。
转换公式为T(K) = T(℃) + 273.15。
2. 压力单位的转换，在不同的国家和地区，压力的单位也有所不同，常见的压力单位有帕斯卡（Pa）、标准大气压（atm）、毫米汞柱（mmHg）等。
在使用物理天然气体积计算公式时，需要将压力转换为统一的单位。
常见的压力单位换算关系为1 atm = 101325 Pa = 760 mmHg。
3. 摩尔数的计算，摩尔数是指一种物质中含有的摩尔数量，通常用符号n表示。
在计算天然气体积时，需要知道天然气的摩尔数。
摩尔数的计算可以通过天然气的质量和摩尔质量来进行。
摩尔质量是指1摩尔物质的质量，常用符号M表示，其单位为g/mol。
摩尔数的计算公式为n = m/M，其中m为天然气的质量，M为天然气的摩尔质量。
通过以上几点注意事项，可以准确地使用物理天然气体积计算公式来计算天然气的体积。
下面将通过一个实际的例子来说明如何使用物理天然气体积计算公式。
假设有一定质量的天然气，其质量为1000克，温度为25摄氏度，压力为1 atm。
天然气液化中采用pr方程的气液相平衡计算
天然气液化中采用pr方程的气液相平衡计算随着工业的发展和技术的进步，液化天然气已经成为一种重要的能源来源，特别是在需要能量的地方。
为了更好地利用天然气液化的能源，人们必须正确地计算气液平衡。
而采用Pr方程计算气液相平衡是计算过程中最重要的一步，因此，有必要重点研究。
Pr方程是指压力与温度的比值， p = T / Tr其中Tr是绝对温度，T是温度，p是压力。
Pr方程与理想气体呈现出一定的关系，因此，在天然气液化过程中，可以用Pr方程来计算压力。
计算过程中，以Pr方程为基础，可以计算出压力和温度的依赖关系，从而准确地计算出气液相平衡条件及它们之间的关系。
由于天然气的性质有所不同，采用Pr方程进行计算时，需要选择合适的参数，即定义温度和压力之间的关系，这取决于天然气特性。
Pr方程所表示的关系可以根据不同类型的天然气来更改，以获得更准确的计算结果。
除了Pr方程外，在天然气液化研究中还应该考虑一系列的因素，如地压，温度等，以及天然气的性质，如沸点，沸程，挥发性等。
这些因素都会影响到最终计算出来的结果，因此，在计算气液平衡时，必须考虑这些因素。
由于Pr方程模型本身的复杂度，以及外部因素的影响，计算气液相平衡条件是一个很复杂的问题，因此，在实际应用中，需要进行大量的实验研究来验证和调整计算结果。
尽管Pr方程计算天然气液化中气液平衡条件的过程比较复杂，但它是重要的一步，它可以提供精确的计算结果，因此，使用Pr方程来计算气液相平衡仍然是必要的。
总之，采用Pr方程计算天然气液化中气液相平衡是一项重要的任务，它对于改善天然气液化工艺及最优化能源使用具有重要意义。
Pr方程可以提供精确的计算结果，但计算过程仍然复杂，必须考虑外部因素的影响，同时需要大量的实验研究。
天然气液化计算技术详解
1.3火用 损失 1.3.1 压缩机火用 损失 • △Ex=
• • • • • • • •
1 G H .2 H 1 1 T0 S .2 S1 / 3600
KW Kg/h KJ/kg KJ/kg KJ/kg. ℃ KJ/kg. ℃ K
2氮膨胀液化工艺计算说 明
2.1 液化工艺流程
氮膨胀工艺流程国内中小型 天然气液化厂大部分用此流 程。 优点： 流程简单，运行管理方便 安全性好，无爆炸危险 技术成熟 缺点： 能耗高，比混合制冷剂液化 流程能耗高40%
2.1.1 天然气流程 经过预处理（脱硫、脱碳、脱水、脱汞）的天然气由1点进入冷 箱。经过1号换热器预冷后脱除重烃，经2、3、4号换热器液化、 过冷，经6点节流降压到24点。24点即为成品LNG，由管路余 压输送到LNG储罐。 2.1.2 制冷剂流程 氮气经压缩机压缩、水冷后，依次进入低压增压机、中压增压 机增压，然后由7点进入冷箱，经1号换热器预冷后进入中压膨 胀机，膨胀后温度、压力均降低，同时输出膨胀功带动中压增 压机。膨胀后的氮气由9点进入3号换热器继续冷却，由10点进 入低压膨胀机，膨胀后温度、压力进一步降低，同时输出膨胀 功带动低压增压机。经第二级膨胀后的制冷剂由换热器的最冷 端11点逐级进入各换热器，为天然气和制冷剂提供冷量。
KW Kg/h KJ/kg. ℃ KJ/kg. ℃ K
式中 △Ex——多股流换热器火用损失 Gi——换热器中第i股流的流量 Si1——换热器中第i股流入口熵 Si2——换热器中第i股流出口熵 T0——环境温度
1.3.5 水冷器火用 损失
Q t2 n / 3600 T0 S 2 S1 G T0 t t1
天然气混合制冷液化流程模拟
资料范本本资料为word版本，可以直接编辑和打印，感谢您的下载天然气混合制冷液化流程模拟地点：__________________时间：__________________说明：本资料适用于约定双方经过谈判，协商而共同承认，共同遵守的责任与义务，仅供参考，文档可直接下载或修改，不需要的部分可直接删除，使用时请详细阅读内容天然气混合制冷液化流程模拟摘要混合制冷剂天然气液化工艺是目前应用最广泛的液化工艺。
本文在分析天然气液化装置中常用的混合制冷剂液化循环的几种基本工艺的基础上，根据天然气和混合制冷剂热物性的特点，选择了PR（Peng-Robinson）方程来计算这两种混合物的相平衡特性。
利用HYSYS软件研究了混合制冷剂流程的冷箱制冷部分，建立了冷箱模拟计算模型，研究了混合制冷剂组分对液化过程的影响。
关键词: 混合制冷液化循环；流程模拟；HYSYS；冷箱Simulation of Mixed Refrigerant Cycle for Natural Gas LiquefactionAbstractThe Mixed-Refrigerant Cycle (MRC) is the most widely used liquefaction process nowadays. Several MRC cycles for Liquefied Natural Gas (LNG) production were analyzed in this paper,--Based on the thermodynamic properties of natural gas and mixed-refrigerant, the Peng-Robinson (PR) equation was selected to calculate the phase equilibrium. The cold box in the liquefaction cycle was simulated by used HYSYS software, the model was established and the components of the mixed refrigerant influence on liquefaction process were studied.Key words: MRC; simulation; HYSYS; cold box目录TOC \o "1-3" \h \z \u HYPERLINK \l "_Toc201665123" 第1章前言 PAGEREF _Toc201665123 \h 1HYPERLINK \l "_Toc201665124" 1.1 工业背景和研究意义 PAGEREF _Toc201665124 \h 1HYPERLINK \l "_Toc201665125" 1.1.1 世界液化天然气工业的发展 PAGEREF _Toc201665125 \h 2HYPERLINK \l "_Toc201665126" 1.1.2 中国液化天然气工业的发展 PAGEREF _Toc201665126 \h 2HYPERLINK \l "_Toc201665127" 1.2 国内外研究现状 PAGEREF _Toc201665127 \h 4HYPERLINK \l "_Toc201665128" 1.3 研究内容 PAGEREF_Toc201665128 \h 5HYPERLINK \l "_Toc201665129" 第2章混合制冷液化流程 PAGEREF _Toc201665129 \h 6HYPERLINK \l "_Toc201665130" 2.1 混合制冷液化流程 PAGEREF _Toc201665130 \h 6HYPERLINK \l "_Toc201665131" 2.2 混合制冷剂液化流程的分类PAGEREF _Toc201665131 \h 6HYPERLINK \l "_Toc201665132" 2.2.1 闭式混合制冷剂液化流程PAGEREF _Toc201665132 \h 6HYPERLINK \l "_Toc201665133" 2.2.2 开式混合制冷剂液化流程PAGEREF _Toc201665133 \h 8HYPERLINK \l "_Toc201665134" 2.2.3 丙烷预冷混合制冷剂液化流程 PAGEREF _Toc201665134 \h 9HYPERLINK \l "_Toc201665135" 2.2.4 CII液化流程 PAGEREF_Toc201665135 \h 12HYPERLINK \l "_Toc201665136" 2.2.5 新型两级混合制冷剂液化流程 PAGEREF _Toc201665136 \h 14HYPERLINK \l "_Toc201665137" 第3章冷箱 PAGEREF_Toc201665137 \h 16HYPERLINK \l "_Toc201665138" 3.1 冷箱简介 PAGEREF_Toc201665138 \h 16HYPERLINK \l "_Toc201665139" 3.2 冷箱的技术关键 PAGEREF _Toc201665139 \h 16HYPERLINK \l "_Toc201665140" 3.2 液化天然气领域冷箱的应用PAGEREF _Toc201665140 \h 16HYPERLINK \l "_Toc201665141" 第4章天然气液化流程模拟软件PAGEREF _Toc201665141 \h 19HYPERLINK \l "_Toc201665142" 4.1 HYSYS简介 PAGEREF_Toc201665142 \h 19HYPERLINK \l "_Toc201665143" 4.2 HYSYS中各个模块的性质与原理 PAGEREF _Toc201665143 \h 20HYPERLINK \l "_Toc201665144" 4.2.1 气液分离器 PAGEREF_Toc201665144 \h 20HYPERLINK \l "_Toc201665145" 4.2.2 壳管式换热器 PAGEREF _Toc201665145 \h 21HYPERLINK \l "_Toc201665146" 4.2.3 LNG换热器 PAGEREF_Toc201665146 \h 24HYPERLINK \l "_Toc201665147" 4.2.4 阀门 PAGEREF_Toc201665147 \h 26HYPERLINK \l "_Toc201665148" 4.3 HYSYS的实际应用 PAGEREF _Toc201665148 \h 26HYPERLINK \l "_Toc201665149" 第5章天然气液化流程模拟PAGEREF _Toc201665149 \h 28HYPERLINK \l "_Toc201665150" 5.1 概述 PAGEREF_Toc201665150 \h 28HYPERLINK \l "_Toc201665151" 5.2 液化流程模拟步骤 PAGEREF _Toc201665151 \h 28HYPERLINK \l "_Toc201665152" 5.2.1 输入条件 PAGEREF_Toc201665152 \h 28HYPERLINK \l "_Toc201665153" 5.2.2 流程搭建 PAGEREF_Toc201665153 \h 30HYPERLINK \l "_Toc201665154" 5.3 流程模拟计算 PAGEREF_Toc201665154 \h 33HYPERLINK \l "_Toc201665155" 5.3.1 收敛计算 PAGEREF_Toc201665155 \h 33HYPERLINK \l "_Toc201665156" 5.3.2制冷剂组分对换热的影响PAGEREF _Toc201665156 \h 34HYPERLINK \l "_Toc201665157" 5.3.3 结果分析 PAGEREF_Toc201665157 \h 35HYPERLINK \l "_Toc201665158" 第6章结论与展望 PAGEREF_Toc201665158 \h 37HYPERLINK \l "_Toc201665159" 6.1 结论 PAGEREF_Toc201665159 \h 37HYPERLINK \l "_Toc201665160" 6.2 展望 PAGEREF_Toc201665160 \h 37HYPERLINK \l "_Toc201665161" 参考文献 PAGEREF_Toc201665161 \h 38HYPERLINK \l "_Toc201665162" 致谢 PAGEREF_Toc201665162 \h 39第1章前言1.1 工业背景和研究意义天然气作为一种清洁优质燃料，是当今世界能源消耗中的重要组成部分，其开发和利用已在全球受到普遍关注[1]。
2.3 立方型状态方程
立方型状态方程的求解方法
（1）利用立方型状态方程的基本形式求解。 已知V、T，求解P 解题思路： ① 查表找出物质的Tc,Pc,ω 等参数，计算出 a,b,α (T)等参数。 ② 计算过程中，统一单位。 ③ 列出相应的状态方程，带入相应的参数。 ④ 根据状态方程所涉及到的未知参数，采用相应的 公式求解出。 ⑤ 带入数据求解。 举例：例2-2。（陈忠秀教材P9例2－1）
举例：例2-4（陈新志教材P18例2－2）。
定义参数h：
b B h V Z
PR方程可以表示成另外一种形式：
1 A h Z 2 1 h B 1 2h h
PV Z RT
2.3.5 立方型状态方程的根及其求解方法
给定T和V，由立方型状 态方程可直接求得P 。但大 多数情况是由T和P求 V 。 当T > Tc 时，立方型状态 方程有一个实根，它是气体 容积。 当T<Tc时，高压下立方 型状态方程有一个实根，它 是液体容积。低压存在三个 不同实根，最大的V值是蒸 气容积，最小的V值是液体 容积，中间的根无物理意义。
231vanderwaalsvdw方程第一个同时计算汽液两相表达临界点的方程其它立方型方程的基础形式简单ab是常数准确度低实际应用少计算常数采用了临界等温线在临界点的条件prtvbav???2关于vdw常数和临界压缩因子zc临界等温线在c点的斜率和曲率等于零??0222?????????????cccctvabvrtvp??0624322??????????????cccctvabvrtvp解方程组得方程常数389cccvbvrta??ccccccccccccvrtvvrtvvrtvabvrtp8389322???????可得到375
（3） 利用立方型状态方程的三次展开式求解 已知T、P，求解V。 ◆在临界点，方程有三重实根，所求实根即为Vc； ◆当T<Tc，压力为相应温度下的饱和蒸气压时，方 程有三个实根，最大根是气相摩尔体积VV （ZV ）， 最小根是液相摩尔体积VL （ZL ），中间根无物理 意义。 ◆其他情况下，方程有一实根和两个虚根，其实根 为 液 相 摩 尔 体 积 VL （ ZL ） 或 气 相 摩 尔 体 积 VV （ZV）。
立方型状态方程
立方型状态方程一、什么是立方型状态方程 1.1 状态方程的定义 1.2 立方型状态方程的特点二、立方型状态方程的数学表达 2.1 立方型状态方程的一般形式 2.2 理想气体状态方程与立方型状态方程的关系三、应用 3.1 立方型状态方程的意义 3.2 立方型状态方程在热力学中的应用 3.3 立方型状态方程在化学反应中的应用四、实例分析 4.1 实例一：计算气体的压强 4.2 实例二：计算气体的体积 4.3实例三：计算气体的温度五、优缺点比较 5.1 立方型状态方程的优点 5.2 立方型状态方程的缺点 5.3 立方型状态方程与其他状态方程的比较六、总结一、什么是立方型状态方程 1.1 状态方程的定义状态方程是用来描述气体或其他物质在特定条件下的状态的方程。
通过状态方程，我们可以获得物质的压强、体积、温度等参数之间的关系。
1.2 立方型状态方程的特点立方型状态方程是一种常见的状态方程形式，可以表示为P = aV^3 + bV^2 + cV + d，其中P为压强，V为体积，a、b、c、d为常数。
立方型状态方程的特点是方程中包含了V的三次方项，因此称为立方型。
二、立方型状态方程的数学表达 2.1 立方型状态方程的一般形式立方型状态方程可以表示为P = aV^3 + bV^2 + cV + d。
其中，a、b、c、d为常数，代表了物质的特性和条件。
2.2 理想气体状态方程与立方型状态方程的关系立方型状态方程是理想气体状态方程的一种推广形式。
理想气体状态方程可以表示为PV = nRT，其中P为压强，V为体积，n为物质的物质量，R为气体常数，T为温度。
通过对理想气体状态方程进行适当的变换，可以得到立方型状态方程。
三、应用 3.1 立方型状态方程的意义立方型状态方程可以用来描述物质的状态和性质，对于研究物质的性质和过程具有重要意义。
通过立方型状态方程，可以计算和预测物质的压强、体积和温度等参数之间的关系。
3.2 立方型状态方程在热力学中的应用在热力学中，我们经常需要计算气体在不同条件下的状态参数。
LKP方程天然气液化中混合物实际焓熵的计算
c4 2 3 2 6 3 ρ r 3 β + γρ r （ 6 － 3 β） － 3 γ ρ r exp (－ γρ r ) T3 r 牛顿迭代式为： ρr = ρr － 式中 f（ ρ r ） f' （ ρ r ）
[
]
}
（ 16 ）
天然气和混合制冷剂的余焓和余熵表达式见式 （ 21 ） 、 （ 22 ） ， 其中上标 （ 0 ） 代表简单流体， 上标 （ r ） 代表参考流体。 H m － H m， H m － H m， 0 0 = RT c RT c
(
（ 21 ）
（ 22 ）
。当计算简单流体的参数时， I = 1； 当
4 式
天然气液化中混合物理想焓熵的计算
［8 ］
I = 2。 计算参考流体的参数时， T r ≤1 时， 当 p r ≤1 、 气相对比密度初值 ρ r， 0
0． 8 = 0． 5 （ 3 － T r ） p r ； 液相对比密度初值 ρ r， 0 = 15 － 8 T r
LKP 方程［8 － 9］
LKP 方程的对比态表达式见式 （ 4 ） ～ （ 8 ） 。 对
于混合物， 还要考虑混合物的临界性质， 其表达式见 式（ 9 ） ～ （ 14 ） 。 Z = Z （ 0） + Zi = ω （ r） － Z （ 0） ） （ r） （ Z ω （ 4）
c4 γ γ β + 2 exp － 2 3 2 v v Tr vr r r B = b1 － C = c1 － D = d1 + p pr = pc b2 b3 b4 － － 3 T r T2 Tr r c2 c3 + T r T3 r d2 Tr T Tr = Tc
b3 2 b4 + 3 T2 Tr r － Tr vr
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
参考文献(14条)
1.李佩铭;焦文玲 液化天然气工业链与国际贸易[期刊论文]-煤气与热力 2006(12)
2.郭东海 混合制冷剂液化天然气循环(MRC)参数选定及优化初探 1992 3.朱开宏 化工过程流程模拟 1993
4.顾安忠;鲁雪生;汪荣顺 液化天然气应用的基础研究[期刊论文]-石油与天然气化工 2001(02)
关键词 天然气液化流程状态方程相平衡计算
随着我国不断引进液化天然气以及大力发展国 内液化天然气产业，液化天然气在我国能源结构中 占据越来越重要的地位。因此，液化天然气技术的 发展显得至关重要。在国外，液化天然气技术的研 究非常成熟，但是由于知识产权及经济利益的限制， 我国对天然气液化方面的研究还只是起步阶段。
已知Ｐ、ｒ、Ｆ和ｚ四个参数，求解ｙ、Ｌ、戈、），，是相 平衡计算的主要任务。在计算ｙ、￡、戈、ｙ和日、Ｓ时， 压缩因子、逸度系数、理想焓、理想熵、余焓、余熵可 以用状态方程或者该状态方程推导得到的相应表达 式求取，日、ｓ的计算则以相平衡的计算为基础。
２ 天然气气液相平衡计算难点
在应用最广泛的混合制冷剂液化流程１９１ （ＭＲｃ）中，该循环以ｃ。至ｃ，的碳氢化合物及Ｎ： 等六种以上的多组分混合制冷剂为工质，进行逐级 冷凝、蒸发、节流膨胀得到不同温度水平的制冷量， 以达到逐步冷却和液化天然气的目的。
式中，Ｐ为压力；Ｒ为摩尔气体常数，Ｒ＝８．１３１４
Ｊ／（ｍｏｌ／Ｋ）；，为温度（Ｋ）；ｙ。为摩尔体积（ｍ３／
ｍ０１）；ｂ、６；、ａ、ｏｉ、ｏｉ、ｍ。为与气体种类有关的常数；ｚ。
为组分ｉ的摩尔分数；ｒ…为组分ｉ的临界温度（Ｋ）；
Ｐ。为组分ｉ的临界压力（Ｐａ）；ｚ，为组分ｊ的摩尔分
数；Ｋｉ为二元交互作用系数；丁。为组分ｉ的对比温 度；叫．为组分ｉ的偏心因子。
式Ⅲ１如下：
ｐ
㈩∽’
其中 ６＝勋。６。６。＝０．０７７８Ｒ丁。／Ｐ。
。仅＝强。乃（Ⅱｄ）ｄ （Ⅱｄ）盯＝（１一Ｋ。）￣／Ⅱ。ｄｉⅡＪａ』
石油与天然气化工
ＣＨＥＭｌＣＡＬ ＥＮＧｌＮＥＥＲｌＮＧ ＯＦ ＯｌＬ＆ＧＡＳ
２００７
ｏｉ＝（ｏ．４５７２４尺２《／Ｐ。）
仅。＝［１＋（ｏ．３７４６４＋１．５４２２６伽。一ｏ．２６９９２叫；）
天然气的液化¨。流程主要分为三种：阶式（级 联式）液化流程、丙烷预冷混合制冷剂液化流程、带 膨胀机的液化流程。目前应用最广泛的是丙烷预冷 混合制冷剂液化流程‘２。。在各种流程中降低液化 天然气的生产成本，都是重点要解决生产过程中的 节能降耗。对于天然气液化的整个过程系统来说， 不能够只是考虑单个设备的节能，而应该考虑整个 过程系统节能‘３。，因此要使天然气液化工艺设计的 能耗最小、费用最小，就必须将整个系统集成起来作 为一个有机结合的整体来看待，达到整体设计最优 化。由此，需要对液化流程进行模拟计算及优 化Ｈ’“。以下便针对天然气液化循环流程中最常用 的混合制冷剂液化流程的热力学相平衡计算∞’７‘进 行分析和讨论。
万方数据
时的精度。ｓＲＫ方程是比较成功的一种改进，其形
ｐ＝器一南 式【１３】为：
ｐ ２可］一瓦万了万
㈩ｏ １’
其中 ６＝昆。６。 ６。＝０．０８６６４Ｒ丁．。／Ｐ．。
。２罾￣／ｎ。ｏ，ｚ。弓（１一Ｋ“）
ｏ。＝（ｏ．４２７４７Ｒ２ Ｚ．／Ｐ。．。）／［１＋ｍ：（１一
￡？）］２
ｍ。＝ｏ．４８＋１．５７４伽。一ｏ．１７６训；
黧组分。豢黼嬲铘衡比 ＭＰ８““ 分数
摩尔分数
摩尔分数
”“ｌ““
紫
撇麓撇蕊辫蕊
ＣＨ。 ０．８４５０ ０．７８４０ Ｏ．７７６８ ０．９０６１ ０．９０７２ １．１５６０ １．１６７９
：，‘．。ｊ ｃ２Ｈ６
Ｃ３Ｈ８
ｏ．１４７６ ０．００７４
ｏ．２０３７ Ｏ．０１２３
ｏ．２１０５ ０．０１２７
ｏ．０９１４ ｏ．０９０２ ｏ．４４８０ ｏ．４２８５ ０．００２５ ０．００２５ Ｏ．２０３０ ０．２００４
ＣＨ。 Ｏ．６１９９ ０．３７２５ Ｏ．３７３６ ０．８６７３ ０．８６５２ ２．３２８０ ２．３１５５
二’。。：ｃ２Ｈ６ ｏ．３２３ｌ ｏ．５１９３ ｏ．５１７４ ｏ．１２７０ ｏ．１２９６ ｏ．２４４０ ｏ．２５０５
Ｃ３Ｈ８ Ｏ．０５７０ ０．１０８２ Ｏ．１０９０ ０．００５７ ０．００５２ Ｏ．０５２０ ０．０４８１
｝
式中，ｍ为ⅣＦ；ｓ为收敛条件判别值；Ⅳ为组 分数。
（２）根据ＳＲＫ方程、ＰＲ方程和以上的平衡方 程，可以求出气液相组分的摩尔分数、气液相平衡常 数等，流程框图见图１。
（３）本文给出使用ＰＲ方程的两组算例。已知 压力、温度、流量和各组分的摩尔分数，统一将流量 设定为１ ｍｏｌ，并将计算结果与参考文献［１］给出的 国外实验数据进行对比，见表１。由表ｌ可以看出， 本文模拟的结果满足一定的精度，符合实际。
5.Vink K J Comparison of baseload liquefaction process 1998 6.Mukhopadhyay M;Awasthi R K-value predictions for the methane -ethane-propane System 1981(06) 7.Soave G Equilibrium constants from a modified Redich-Kwong equation of state 1972 8.冯霄 化工节能原理与技术 2004 9.石玉美;顾安忠 天然气液化流程中的气液相平衡计算 1997(03) 10.童景山 化工热力学 1995 11.Evans L B Advances in process flowsheeting systems 1985 12.Finn A J;Johnson G L;Tomlinson T R Development in natural gas liquefaction 1999(04) 13.Terry Lavin Comparison of liquefaction process 1998(03) 14.顾安忠 液化天然气技术 2004
（２）求解状态方程中的系数：在ｓＲＫ和ＰＲ方 程中，利用已知量求出状态方程中的系数，如ｂ、ａ、
Ａ、Ｂ。
（３）计算压缩因子：由步骤（２）计算出的结果， 然后根据式（２）和式（５）求出压缩因子ｚ，计算方法 可选用迭代法。
（４）计算逸度系数：状态方程计算气液相平衡 的关键在于求出两相的分逸度系数，将上述步骤求 出的结果代人式（３）和式（６）求出逸度系数。
万方数据
上述步骤完成后，得到状态方程求出的结果，该 结果用于求解在已知温度ｒ，压力尸的给定条件下 各相的摩尔分数和摩尔流量。
４ 给定条件下的ｓＲＫ、ＰＲ方程的计算
已知温度丁，压力Ｐ条件下的相平衡计算即
（ＪＰ，丁）闪蒸计算¨’”１，用于确定各相的摩尔分数、摩
尔流量，从而为焓和熵的计算提供已知参数。求解
石油与天然气化工
ＣＨＥＭｌＣＡＬ ＥＮＧＩＮＥＥＲＩＮＧ ＯＦ ＯｌＬ＆ＧＡＳ
天然气液化流程模拟中立方型状态方程的计算８
李佩铭 焦文玲 （深圳市燃气集团股份有限公司）
摘 要 天然气液化流程广泛采用的是丙烷预冷混合制冷剂液化流程。为了进一步优化流 程，减少能源的消耗，需要对整个流程进行模拟，而模拟过程中热力学参数的计算便是整个流程计 算的基础。本文讨论使用两种立方型状态方程（ｓＲＫ和ＰＲ方程）对热力学参数进行相平衡计算， 为后续计算焓、熵等参数提供相应的解决办法，并判断选择的状态方程是否符合要求。
度系数咖｛时，ｚ。为戈：；计算气相逸度系数咖；时，ｚ。为
ｙ。。
３．２ ＰＲ方程‘１２３
由于ＲＫ方程和ｓＲＫ方程预测液体密度时精
度较差，Ｐｅｎｇ和Ｒｏｂｉｎｓｏｎ提出了改进形式的状态方
程。由于ｓＲＫ和ＰＲ方程有很好的温度函数，在预
测蒸汽压时有明显优点，而在预测稠密区的摩尔体
ｐ２＝可嚣：］一一万万南砰 积方面，ＰＲ方程比ｓＲＫ方程更优越。ＰＲ方程形
（１一砟５）］２
式中，ｏｄ、ｏ。、ｄ。均为与气体种类有关的常数。 压缩因子大小直接影响到各逸度的大小，从而
也影响气液相平衡计算结果值，也即气相流量ｙ、液
相流量￡、液相摩尔分率ｚ蒯、气相摩尔分率ｙ删，进 而影响到焓熵值的正确性，所以正确计算压缩因子
是系统计算是否可靠的关键因素之一。 ＰＲ方程用压缩因子ｚ表示为：
１ 天然气和混合制冷剂的热物性参数。１＇钊
在天然气液化循环流程的热力学计算中，天然 气及制冷剂的热物性计算是整个流程计算的基础。
因此，是否得到准确的天然气及混合制冷剂的热物 性参数，成为衡量流程模拟计算过程准确度的标准 之一。
由于天然气和制冷剂组分类似，而且状态接近， 求解各自的热物性参数时，采用相同方法进行。在 天然气液化流程中，涉及到的物性参数¨。有压力Ｐ、 温度丁、总流量Ｆ、总流量中的摩尔组分ｚ、气相流量 ｙ、液相流量￡、液相摩尔分率戈、气相摩尔分率ｙ、焓 日，熵｜ｓ。除了以上的物性参数外，还会涉及压缩因 子、逸度系数、理想焓、理想熵、余焓、余熵等参数。
立方型方程的形式并不是其计算精度的主要原 因，方程中参数的选取才是提高方程的精度和广泛 适用性的关键。因此选用常用的ＳＲＫ方程和ＰＲ方 程来分别计算天然气和混合制冷剂的相平衡特性并 进行比较分析。
３．１ ｓＲＫ方程㈦１３１ ＳＲＫ方程是Ｓｏａｖｅ对ＲＫ方程的改进，ＲＫ方程
于１９４９年提出，得到了广泛的应用。该方程是由德 瓦尔斯方程衍生而来。由于ＲＫ方程应用方便，许 多研究工作者都想在原方程基础上进一步修正，以 便改进其精度，特别希望改进用它计算气液相平衡
Ｚ３一（１一日）Ｚ２＋（Ａ一３日２—２Ｂ）Ｚ一
（ＡＢ—Ｂ２一曰３）＝０
（５）
式中，压缩因子ｚ＝ＰＶ／（ＲＴ）；Ａ＝ａ仅Ｐ／（ＲＴ）２；
Ｂ＝ｂＰ／（ＲＴ）。
ＰＲ方程的逸度系数的表达为：