上海煤气 2019年第4期〈〈 1
LNG 热物性参数的计算
同济大学机械与能源工程学院 罗 洋 周伟国 贾云飞
摘要：液化天然气(LNG)是多元混合物，其热物性参数的计算不能按照理想状态处理。引入各参数相应混合规则，分别介绍了LNG 导热系数、动力黏度、密度及定压比热容的计算方法，并且以某气源为例，用Matlab 编程计算出不同温度下LNG 的这四种基本热物性参数，并和HYSYS 的计算结果进行了比较，结果表明该计算方法可行。
关键词：LNG 热物性参数 计算
天然气液化后的体积约为同质量气态天然气体积的1/625，大大节省了储存空间及运输成本，具有明显的经济优势和发展潜力。天然气在其液态状态输送过程中，其热物性参数会随着温度和压力的变化而不断变化。较为准确的热物性参数是天然气在液化、换热和运输等环节流程模拟及动态分析的基础，同时也是提高流程模拟分析准确性的关键。本文主要针对LNG 导热系数λ、动力黏度µ、密度ρ和定压比热容Cp 四个热物性参数的计算方法进行总结和验证。
1 混合物热物性计算方法
天然气是多元混合物，分子之间的尺寸、形状和极性等存在较大差异，分子间相互作用与纯物质中同种分子间的相互作用存在本质上的差别，如果用纯组分热物性、摩尔分数以及分子量等参数建立函数关系计算混合物热物性，则忽略了各组分之间分子作用力，导致计算出现偏差；并且方程中一些参数值或参数的关联式仅适合于纯物质，用状态方程处理混合物体系时，离不开方程中有关参数的混合规则；其次，对于某些组分的热物性参数数据缺乏时，只能采用估算值进行计算，使得计算值和实测值偏差较大。天然气作为混合物，需要综合考虑各组分之间的相互作用对于整体效应的影响。
由此需要引进混合规则，根据对应态原理，混
合物可以看作具有一套按一定规则求出的假临界参数、性质均一的虚拟的纯物质，其对应的物性参
数需要通过混合规则求出。由于这种假想参数强烈依赖于混合物的成分，因此完全由实验确定是非常困难的。现在普遍的做法是，利用各种混合规则而由纯物质的参数求混合物的假想参数。
混合规则即为混合物的虚拟参数与其组成的各纯物质参数之间的关联式，在此采用LNG 各组成的由统计力学得到的二次型混合规则求得，即：
(1) 式中：
——
混合物的假想参数； ——组分i 的摩尔分数； ——组分j 的摩尔分数；
——表征组分i 和j 的相互作用参数。
2 LNG 热物性参数计算方法
2.1 导热系数
导热系数直接影响着传热传质的计算，所以对于混合物选择适当的计算方法尤为重要。液体混合物的导热系数一般由单组分热导率通过混合规则导出，目前较为成熟的混合物热导率模型多针对二组分混合物，多组分液体混合物的热导率公式相对较少，以Li 模型(参看C. C. Li 发表于AIChE Journal 1976年第5期的文章Thermal Conductivity of Liquid Mixtures)介绍一下较为方便、准确。先对混合液体做以下假设：
(1)液相中，能量输送依靠的是分子之间的碰撞，如果液体混合物与它的临界状态相差甚远，即可假设碰撞的频率与相邻分子的数量和尺寸成正
2 〉〉2019年第4期 上海煤气
比，则对于热力学状态而言，碰撞过程的近似加权参数应该是体积分数，而不是摩尔分数；
(2)相互导热系数可以通过纯组分调和平均值来近似求得。
基于以上假设，本文采用热导率经验关联式与Li 模型结合使用的方法计算LNG 的导热系数。
LNG 各组分的热导率可以由以下经验公式计算，
有机物采用公式(2)，
无机物采用公式(3)，
公式中的相关参数的取值见表
1。
(2) (3)
Li 模型如下：
(4) (5) (6)
式中：——
组分i 的体积分数；
——组分i 纯液体的摩尔体积；
——LNG 混合物的导热系数；
——组分i 的导热系数，W/(m·K)。
表1 导热系数计算公式参数值
2.2 动力黏度计算
黏度是工程设计中必不可少的基础数据，目前黏度计算的方法主要采用经验关联式法。
目前较好的研究方法是在甲烷迁移参数计算的基础上，通过引入混合规则进行虚拟纯物质的修正。
目前常见的LNG 动力黏度计算方法主要有Jamieson-Lohrenz 经验关联式、Teja-Rice 法、
Eyring 动力黏度模型和统一动力黏度法。本文采用的Teja-Rice 法对于在完全非极性到高极性水-有机物混合物范围内的许多混合物，均可得到良好的结果，计算误差平均为
2%。
Teja-Rice 法中采用的混合规则规定如下：
(7)
(8)
(9) (10)
(11) (12) (13)
式中：
——混合物的虚拟临界摩尔体积； ——混合物的虚拟临界温度；
——混合物的虚拟临界压力；
——混合物的虚拟临界相对分子质量；
——临界压缩因子； R ——摩尔气体常数，R=8.314 5 J/(mol K)；
——重量平均摩尔分子质量；
——重量平均相对分子质量。
Teja-Rice 法计算LNG 黏度的公式如下：
(14)
(15)
(16)
(17)
式中：——LNG 混合物黏度，Pa·s ；
——LNG 混合物的特性参数； ——组分i 的特性参数。
在公式(16)中，上标r 1、r 2代表两种参考流体，可选天然气中摩尔组分最大的两种组分。黏度值r1是参考流体r 1在温度T[(T c )r 1/T c,m ]下求得的，r2是参考流体r 2在温度T[(T c )r 2/T c,m ]下求得的。本文取r 1为甲烷，r 2为乙烷。液体黏度的经验公式参数见表2，天然气组分物性数据见表3。
表2 液体黏度计算公式参数值
