LNG船用BOG再液化工艺系统效能优化

2011-02-26方江敏黎志昌

低温工程 2011年5期

谈震方江敏黎志昌

(华南理工大学机械与汽车工程学院广州 510640)

1 引言

随着LNG在全球能源市场所占份额日益扩大，LNG(liquefied natural gas)海运行业得到了空前的发展。在LNG海运过程中如何提高BOG(boil off gas，LNG蒸发气)的回收利用率，成了当前国际建造LNG运输船的前沿课题［1］。为提高LNG海运效率，降低运输成本，用BOG再液化系统工艺流程对LNG船进行优化设计有着巨大的运用潜力。

2 模拟流程

2.1 流程建模及说明

LNG船用BOG再液化装置工艺所涉及的单元设备有压缩机、海水换热器、换热器、多股流换热器等［2］。在用ASPEN PLUS化工模拟软件对LNG船用BOG再液化装置进行整个工艺流程模拟计算时，分别选用如下的模块:压缩机(Cpr)、冷却器(Cooler)、双股流换热器(HeatX)、多股流换热器(MHeatX)、节流阀(Valve)。使用上述模块对LNG船用BOG再液化流程建立模型图，如图1所示。

图1 工艺模拟流程图Fig.1 Picture of simulation process

再液化系统采用复叠式制冷系统，由BOG制冷循环、丙烯制冷循环和混合制冷剂制冷循环构成，即由一个以液货为制冷介质的开式循环和两个分别采用丙烯、混合制冷剂的闭式制冷循环组成。3个制冷循环的说明如下:

(1)BOG制冷循环:LNG液货舱中的BOG进入回热器中，为混合制冷剂的过冷提供冷量，再由BOG压缩机压缩，压缩后的高温BOG通过海水换热器1带走热量，降温后的BOG进入丙烯多股流换热器预冷，再通过BOG深冷换热器进行深冷，最后由BOG节流阀节流降温。节流后的BOG冷凝液进入气液分离器，液相直接返回液货舱，气相BOG与LNG液货舱中的BOG汇合在此进入下一次BOG制冷循环。

(2)丙烯制冷循环:丙烯制冷剂经过丙烯压缩机压缩后，进入海水换热器2降温，再由丙烯节流阀节流降温。进入丙烯预冷多股流换热器，为BOG和混合制冷剂的预冷提供冷量，换热后进入下一轮循环。

(3)混合制冷剂制冷循环:混合制冷剂由混合制冷剂压缩机压缩后，依次经过海水换热器3、预冷多股流换热器、混合制冷剂换热器、深冷多股流换热器、回热器进行换热直至被液化过冷。经回热器后的过冷混合制冷剂由混合制冷剂节流阀降压降温后，进入深冷多股流换热器提供冷量，流出的混合制冷剂再利用剩余冷量为混合制冷剂换热器提供冷量后，进入下一轮混合制冷剂循环。

2.2 流程性能参数计算公式

通过软件的模拟计算得到各节点的压力、温度、密度、质量流量等热力学参数。运用表1中所示的计算公式，计算出工艺流程分析中用到的重要性能参数。

表1 重要性能参数计算公式Table 1 Important performance parameters

3 流程优化

3.1 优化目标

在流程优化问题中，目标函数十分重要。目标函数不同，液化流程的性能以及设计变量的最终优化的结果相差将会很大［3］。追求高效、低能耗、低投资是再液化流程设计的目的。所以应以降低流程功耗和增加再液化系统中3个制冷循环的制冷效率作为优化的目标。本文以BOG压缩机出口压力、丙烯预冷温度、丙烯压缩机出口压力、丙烯制冷节流阀出口压力、混合制冷压缩机出口压力作为优化的变量，以流程中压缩机功耗最小为优化的目标函数。以此得到如下的优化目标函数:

式中:Pmin为总功耗的最小值;a为BOG压缩机出口压力;b为丙烯压缩机出口压力;c为丙烯预冷温度;d为丙烯制冷节流阀出口压力;e为混合制冷压缩机出口压力;Pcpr1为BOG压缩机功耗;Pcpr2为丙烯压缩机功耗;Pcpr3为混合制冷剂压缩机功耗。

以 a、b、c、d、e为优化变量，在优化计算中寻找其最佳取值，使得在该最佳取值状态下，工艺流程的总功耗最低。

3.2 约束条件

在进行流程优化时，为了让优化方法得到的参数进行流程计算时具有物理意义，在优化时设定如下约束条件:

(1)各制冷循环的压缩机入口工质必须为气相;

(2)丙烯制冷循环中，经海水换热器冷却后的丙烯制冷剂，至少被部分液化;

(3)海水换热器入口温度必须高于出口温度;

(4)BOG制冷循环中，BOG经深冷多股流换热器HX4冷却后，至少被部分液化;

(5)混合制冷剂制冷循环中，混合制冷剂经回热器HX1冷却后，至少被部分液化;

(6)根据热力学第二定律，在优化计算中必须设定各换热器的熵增大于零;

(7)各节流阀都产生降温，并且有气液两相产生;

(8)产生目标函数的最优值时对应生成的一组流程已知参数，以此为已知参数计算整个流程性能时，各换热器的热量要平衡［4］。

3.3 优化方法

由于优化计算中优化变量有5个，求解难度较大。为了简化求解，采用变量轮换法［5］也称为交替方向法，基本原理较为简单，它对目标函数的解析性质没有特别的要求，适用于各变量之间无本质上联系或沿坐标轴方向搜索比较容易的特殊结构。图2为使用变量轮换法进行优化计算的框图。为了使优化的结果精确，在优化计算中设定了变量的变化值，即优化变量步长Δx，如表2所示。