天然气液化混合冷剂配方优化研究

孟毅明 王东军 陈 博

1. 中国石油大学(北京) 2. 中国石油天然气管道局天津设计院

天然气液化混合冷剂配方优化研究

孟毅明1王东军2陈 博2

1. 中国石油大学(北京) 2. 中国石油天然气管道局天津设计院

混合冷剂制冷是目前天然气低温液化的重要手段，由于操作工况复杂、冷剂介质多相共存，制冷剂的选择及配方优化已成为天然气混合冷剂液化技术面临的主要问题。本文采用均匀设计实验方法对混合冷剂配方进行研究，首先依据不同制冷剂组分在不同温区制冷的原理，初步选定基本冷剂组分；然后通过HYSYS模拟不同组成混合冷剂的制冷效果，通过实验结果分析混合冷剂中各组分的作用；在保持操作压力及处理量不变的条件下，以冷剂的最小循环量为目标参数进行优化分析，最终形成混合冷剂优化配方，从而避免了采用复杂原理探讨和试算的传统方法。

天然气液化 混合冷剂 流程模拟 配方 均匀设计 优化

随着环境污染问题的日益加剧，天然气作为一种优质清洁燃料，在能源、交通、化工等领域的应用越来越广泛。液化天然气便于运输和使用，因此，开展天然气液化技术研究对扩大天然气的应用范围具有十分重要的意义。

目前，天然气液化主要通过制冷实现，大致分为级联式、膨胀式和混合冷剂3种方式。与前两者相比，混合冷剂技术具有流程简单、机组设备少、投资少、能耗低等特点。目前世界上80%的基本负荷型天然气液化装置都采用了混合冷剂液化流程[1]。

由于操作工况复杂、冷剂介质多相共存，天然气混合冷剂液化技术面临的主要挑战是制冷剂的组成确定及配方优化[2]。鉴于混合冷剂制冷过程中各单纯制冷工质间复杂的热力学交互作用，采用高效的均匀设计确定初步试验方案，通过HYSYS流程模拟软件对混合制冷过程进行模拟研究，归纳预测混合冷剂优化组成并进一步模拟验证。

1 混合冷剂天然气液化流程

混合冷剂天然气液化流程包括制冷循环和天然气液化，具体工艺流程见图1[3]。

1.1 制冷循环

低压混合冷剂气体首先进入一级压缩机C101压缩，经出口冷却器冷却降温，使部分高沸点组分凝结后进入气液分离器V-l。分离器底部液相经预冷换热器换热并通过节流阀VLV-1节流后，与来自主换热器LNG-2的逆流低压气体汇合，作为预冷换热器LNG-l的冷流；分离出的气相混合冷剂，经高压压缩机C102压缩，由换热器E-2降温后进入分离器V-2。分离器底部液相经预冷换热器后进入主换热器LNG-2中冷却，经节流阀VLV-3节流后与来自过冷换热器LNG-3的低压冷剂汇合，并返回主换热器LNG-2提供冷量；气相先后进入主换热器LNG-2和过冷换热器LNG-3完成冷凝和过冷，混合冷剂经节流阀VLV-4节流降温后为天然气和冷剂自身提供冷量。

1.2 天然气液化

净化后的天然气先后经预冷换热器LNG-1、主换热器LNG-2逐步冷却，然后进入气液分离器V-4脱除气体中少量重组分。重组分返回换热器E-3为冷剂预冷提供部分冷量，天然气中的轻组分则继续进入过冷换热器LNG-3冷凝、液化，液化天然气经节流阀VLV-5节流降压分离出不凝气后注人LNG储罐。

2 天然气液化混合冷剂配方研究

混合冷剂组成的确定目前主要采用优化模型计算法[4-5]和原理优化法[6]。计算法是通过选取实验影响因素编写优化方程，利用热力学公式优化冷剂组成，该方法需要对制冷系统建立复杂的热力学模型，优化计算困难;原理优化法是结合流程模拟和原理判断，通过提高或降低混合物中某一种组分的含量来最终确定混合物组成，该方法需要大量的试算且不易寻求最优结果。鉴于混合冷剂制冷过程是复杂的热力学效应交互作用的结果，本文采用原理初选结合均匀实验设计的混合冷剂配方优化方法，通过HYSYS软件对混合制冷过程进行流程模拟，利用少量具有代表性的实验结果建立考察指标与冷剂组成间定量关系，进而预测混合冷剂的优化组成，并对预测结果进行模拟验证。

2.1 混合冷剂初选

混合冷剂制冷主要是通过蒸发吸热及复杂的换热过程实现的，这与阶式制冷的原理相类似[7]，即高沸点冷剂组分主要在高温区发挥制冷效果，低沸点冷剂组分主要在低温区发挥制冷效果。不同的是，混合冷剂中各组分可共同在单压缩循环下逐级蒸发、冷凝，实现不同沸点组分在不同温区制冷。

各个组分都有其特定的制冷温区范围，制冷温区位于-180～200 ℃的多种冷剂组分在0.3 MPa及4.0 MPa压力下的温焓特性如图2所示。

由图2可知，各制冷剂组分均存在特定的制冷温区，且随着冷剂组分的相对分子质量的增加，冷剂组分的沸点逐渐升高，相应的制冷温区也向较高温区移动，气化潜热增大，液化所需的功耗较大。相邻组分间的制冷温区有部分重叠现象。相比而言，在同一制冷温区，为达到相同制冷效果，相对分子质量高的组分功耗小但制冷效果不足；相对分子质量低的组分制冷效果较好，但功耗大。因此，需均衡选择适宜的冷剂组分，温区过多重叠的组分制冷效果是相似的，在冷剂筛选中应尽量避免温区重叠。

为满足常规天然气液化要求，需选择-180～200 ℃的温区范围。对于N2及C1～C5制冷体系，首先选择N2、CH4组分作为低温区制冷组分。鉴于C2H4与CH4及C3H6的制冷温区重叠较小，因此保留C2H4组分，剔除同温区的C2H6组分。鉴于C3H6及C5H12能够连续覆盖中高温区，因此，剔除中间温区重叠的C4H10组分。综上所述，最终选择N2、CH4、C2H4、C3H8、C5H125种组分作为混合冷剂的基本组成。

2.2 混合冷剂配方均匀实验研究

均匀实验设计的实验点具有代表性强、试验次数少的特点。通过数学拟合，可建立考察指标与冷剂组成的定量关系，进而实现定量分析[8]。

本文首先确定实验影响因素及优化指标，再结合流程模拟对不同组成的配方进行均匀实验，并对实验结果进行优化分析及验证，最终确定适宜的混合冷剂配方。

流程模拟参考图1中的混合冷剂液化流程，固定操作压力、处理气量等条件，考察不同混合冷剂配比对制冷效果的影响，并以最低冷剂循环量作为考察不同配方混合冷剂效果的指标，进行均匀实验。

根据某工程实际，选取原料气的组成及操作条件分别如表1和表2所示。

表1 原料天然气组成Table1Componentsoffeednaturalgas组分HeN2CO2CH4C2H6C3H8i-C4n-C4i-C5n-C5n-C6y/%0.020.640.001996.052.460.520.10.10.050.020.03

表2 混合冷剂天然气液化工艺操作参数Table2 Operatingparameterofnaturalgasliquefactionprocessbymixedrefrigerant天然气处理量/(108m3·a-1)冷剂低压操作压力/MPa冷剂高压操作压力/MPa最低液化温度/℃1.810.34.0-160

现利用初步选择的5种冷剂组分进行配方评选实验。参考APCI[9]提出的混合制冷剂构成的一般原则，确定其摩尔分数及变化范围如表3所示。

表3 实验因素及变化范围Table3 Experimentalfactorsandtheirvariationranges冷剂组分y(冷剂组分)/%y(变化范围)/%N2X10.15～0.21CH4X20.28～0.33C2H4X30.12～0.20C3H8X40.11～0.20C5H12X50.06～0.10

通过对组成如表4所示的混合冷剂进行模拟实验，调节冷剂循环量以满足天然气液化所需冷量，同时保证HYSYS工艺模型运行正常。将不同组成混合冷剂实验得到的最低冷剂循环量分别记录，实验结果如图3所示。

表4 均匀实验设计表Table4 Uniformexperimentaldesigntable实验序号y(混合冷剂组成)/%X1X2X3X4X5He10.1560.2950.1520.1550.0960.14620.1620.3100.1840.2000.0880.05630.1680.3250.1280.1460.0800.15340.1740.2850.1600.1910.0720.11850.1800.3000.1920.1370.0640.12760.1860.3150.1360.1820.1000.08170.1920.3300.1680.1280.0920.09080.1980.2900.2000.1730.0840.05590.2040.3050.1440.1190.0760.152100.2100.3200.1760.1640.0680.062

由图3可知，不同配方混合冷剂的制冷效率不同，其中10#配方混合冷剂的循环量最低，为2 197 kmol/h。为便于优化分析，采用式(1)的计算值作为优化指标，均匀实验结果及对应的优化指标如表5所示。

表5 液化天然气混合冷剂配方均匀设计实验结果Table5 Uniformdesignexperimentalresultsbymixedrefrigerantformulation实验序号冷剂最小循环量/(kmol·h-1)优化指标(Y)1#255739.10832#222944.86323#280235.68884#250139.98405#253239.49456#247440.42047#276836.12728#221545.14679#326530.627910#219745.5166

(1)

式中，Y为优化指标；M为转换因子，取100 000 kmol/h；L为循环量，kmol/h。

2.3 实验结果优化

对表5中的实验数据进行二次多项式回归处理，得到的回归方程如式(2)[12]：

Y=-57.501 423 2-412.103 619 5X2+

1 653.773 145 5X4+212.604 388 86X3X3-

4 370.907 584X4X4+1 580.076 157 8X1X2-

2 926.197 693 9X1X4+1 190.255 918 7X2X4+235.380 408 94X3X4

(2)

对该方程进行显著性检验得出：其相关系数R为0.999 9，显著性水平p为0.006 5，剩余标准差S为0.042。因此，该回归方程在实验范围内，回归方程能较好地描述混合冷剂组成与最小冷剂循环量的关系。同时，X5未引入方程，表明在约定范围内，相应的C5H12含量对结果影响相对不明显，可用于组分归一化调整。

以各组分合理范围为约束条件，以优化指标(Y)最大化为优化目标，对式(2)进行最优化求解，得出冷剂组分的优化配方见表6。优化计算得到的各组分摩尔分数加和小于1，为保证N2、CH4、C2H4、C3H8的组成为优化值，需在0.6%～1%范围内适当调整X5含量以使得组成归一化。如表6所示，优化计算得到的最小冷剂循环量为2 014.37 kmol/h。对优化配方进行模拟实验验证，得到最小冷剂循环量为1 976 kmol/h。结果表明，实验值与优化计算值较为接近，且低于均匀设计的全部实验点，从而得到约束条件下天然气液化混合冷剂的优化配方。通过拟合方程(2)可对该制冷系统不同冷剂配方下的最小冷剂循环量进行预测，同时也可对不同条件的混合冷剂配比进行优化。

表6 混合冷剂制冷液化天然气均匀设计实验优化结果Table6 Uniformdesignexperimentaloptimizationresultsbymixedrefrigerantformulation项目冷剂循环量/(kmol·h-1)y(冷剂组分)/%X1X2X3X4X5优化计算2014.370.210.330.20.16960.0667实验验证19760.210.330.20.16960.0904

2.4 因素分析

根据得到的拟合方程(2)，将其他组分固定为最优配比的条件下，分析各冷剂组分含量的变化对制冷效果的影响，结果如图4所示。

由图4可知，C3H8因处于过渡温区，适度增加C3H8含量能够降低冷剂循环量，提高制冷效率;但当其增加至一定值后(最高点对应的冷剂组成)，过高的C3H8含量则会降低总体的制冷效率。在目前冷剂组成范围内，N2、CH4、C2H4含量的增加均有利于降低冷剂循环量，增强冷剂制冷效果。

表7中的显著性检验结果也佐证了上述结论。

表7 实验因素显著性检验Table7 Significancetestofexperimentalfactors因素偏相关系数t检验值显著性水平pX2-0.999741.04020.0006X40.999740.38810.0006X3X30.996712.33930.0065X4X4-0.999978.35420.0002X1X20.999532.87730.0009X1X4-0.999531.51700.0010X2X40.999428.36340.0012X3X40.98916.71160.0215

由显著性检验结果可知，X5项未被选入方程。因此，在加量范围内C5H12的加入对天然气液化的影响并不显著。根据显著性水平p的大小可知，各因素及其交互作用对冷剂循环量的影响由大到小的顺序为：X2gt;X4gt;X1X2gt;X1X4gt;X2X4gt;X3gt;X3X4，其中的乘积项表明，N2、CH4、C2H4、C3H8等冷剂组分间存在一定的交互作用。

3 结 论

本文通过对天然气液化混合冷剂进行理论和模拟实验研究，形成了原理初选、模拟实验与均匀设计优化分析相结合的天然气液化混合冷剂配方高效评选方法，主要研究结论如下：

(1) 不同制冷剂组分在不同温区制冷，应根据天然气组成确定制冷温区，从而选择对应温区的混合冷剂组分。

(2) 在同一制冷温区，为达到相同制冷效果，分子量高的组分功耗小但制冷效果不足；分子量低的组分制冷效果较好，但功耗大。因此，需均衡选择适宜的冷剂组分，同时在冷剂组分筛选中应尽量避免冷剂组分的温区重叠。

(3) 提出采用原理初选、HYSYS模拟实验与均匀设计优化分析相结合的天然气液化混合冷剂配方高效评选方法，即通过均匀设计选取少量具有代表性的模拟实验，回归分析混合冷剂中各组分的作用，并通过优化计算确定混合冷剂的优化配方，对优化配方模拟实验验证表明，此方法是可行的。

[1] 赵敏, 厉彦忠. C3/MRC液化流程中原料气成分及制冷剂组分匹配[J]. 化工学报, 2009, 60(增刊1)： 50-57.

[2] 曹文胜, 吴集迎, 鲁雪生， 等. 撬装型混合制冷剂液化天然气流程的热力学分析[J]. 化工学报, 2008, 59(增刊2)： 53-59.

[3] 王勇, 张玉玺, 白剑锋. LNG制冷HYSYS计算模型研究[J]. 天然气与石油, 2012, 30(4)：30-32.

[4] 尹全森. 混合制冷剂循环优化设计和动态特性研究[D]. 哈尔滨： 哈尔滨工业大学, 2010.

[5] GEOFFREY G H，CHEN J. University of Leeds, Leeds LS2 9JT, UK:A computer simulation program for mixed-refrigerant air conditioning[J]. Energy, 2013: 46-155.

[6] 蒋浩, 高永和, 王文军， 等.MRC流程混合制冷剂的选择[J]. 煤气与热力, 2012, 32(2)： 4-6.

[7] 张镨, 鹿来运, 郭开华. 一种混合工质遴选、配比方法探讨及实验验证[J]. 低温工程, 2013(6)： 26-30.

[8] 毛东森. 均匀设计及其在化学化工中的应用[J]. 化工进展, 1997(5): 16-19.

[9] 牛亚楠. 多元混合制冷剂小型天然气液化装置的模拟研究[D]. 上海： 同济大学, 2007.

[10] 王柱, 方开泰. 带有定性因素均匀设计的均匀性度量准则[J]. 数据统计与管理, 2000, 19(3): 28-32.

[11] 高辉, 胡良平, 郭晋， 等. 如何正确处理正交设计和均匀设计定量资料[J]. 中西医结合学报, 2008, 6(8): 873-881.

[12] 姚钟尧. 回归分析法在均匀设计数据分析中的地位[J]. 特种橡胶制品, 1998, 19(2): 35-41.

Studyonmixedrefrigerantformulationoptimizationfornaturalgasliquefaction

MengYiming1,WangDongjun2,ChenBo2

(1.ChinaUniversityofPetroleum(Beijing),Beijing102249,China)(2.ChinaPetroleumPipelineBureauTianjinDesignInstitute,Tianjin120107,China)

Mixed refrigerant method is an important means to liquefy natural gas at low temperature. Due to the complex operating condition and cooling medium multiphase coexistence,selection and formulation optimization of mixed refrigerant become the major challenge for the natural gas liquefaction technique. This article used uniform design experiments to study mixed refrigerant formulation. That is, preliminarily select the basic components by the principle that different refrigerant has an effect in its particular refrigerating range; then simulate the refrigerating effect of different components by using process simulation software ASPEN HYSYS, and analyze the function of each component with less experiments; under the constant pressure and treatment capacity, taking the minimum circulation quantity of mixed refrigerant as target parameter to make optimal analysis, finally obtain the optimization formulation of mixed refrigerant, which avoid adopting the conventional approach with complex principle analysis and trial calculation.

natural gas liquefaction, mixed refrigerant, process simulation, formulation, uniform design, optimization

孟毅明(1991-)，男，天津市人，中国石油大学(北京)化学工程学院，硕士研究生。E-mailmengym7@163.com

TE624.9

A

10.3969/j.issn.1007-3426.2015.03.014

2015-04-24；编辑康莉