工业低温余热回收用于多周期供冷的设计与优化

高 为, 周宇昊, 张凌玮, 王彧斐

(1. 华电电力科学研究院有限公司, 浙江 杭州310030; 2. 浙江省蓄能与建筑节能技术重点实验室,浙江 杭州310030; 3. 中国石油大学(北京) 化学工程与环境学院, 北京 昌平102249)

1 前 言

在我国工业领域中，常有大量低温余热资源没有得到回收利用。目前对低温余热的定义不尽相同[1-3]，一般可认为200 ℃以下的余热为低温余热。在我国石化行业中，低温余热的占比可达80%[3]。由于温度等级低，回收利用技术难度高，低温余热的浪费情况非常严重。据估计，我国北方每年损失近7.6×1014J，约合2.6 亿吨标准煤当量(tce)的低温余热[4]。与此同时，大量化石燃料燃烧产生高品质能源，以满足城市建筑物的能源需求。如果能在城市供热供冷中有效地利用工业低温余热，则可以减少化石燃料的消耗和温室气体的排放量，并提高能源的利用效率[5]。当前已有一些研究致力于将低温余热整合到区域供暖系统中，例如，Fang 等[2]提出一种估算低温余热利用潜力的方法，Tveit 等[6]为一座造纸厂与市政能源系统的整合设计换热网络。

在传统的建筑物制冷方式中，独立的小型制冷机用于满足每个房间或建筑物的冷量需求，如家用空调。近年来，利用大型制冷站集中制冷并将冷量输送至周边用户的区域供冷技术得到了重视和发展。Söderman[7]提出一种包括冷却设备在内的区域供冷优化模型。中央冷却设备的制冷方式可分为压缩式制冷和吸收式制冷2 种，相比于主要由电驱动的压缩式制冷，吸收式制冷能够利用低温余热资源来驱动制冷循环的运转[3,8]，从而提高能源的利用效率，这使得吸收式制冷技术在区域供冷中具有较强的吸引力。相比于直接利用低温余热进行区域供暖，将低温余热用于制冷供冷的研究较少，而且很少有研究能在考虑用户需求随季节波动的基础上给出各周期下具体的操作方案。

本研究旨在通过回收工业低温余热来提供驱动吸收式制冷机所需的热量，从而满足用户的冷量需求。本文的主要工作是建立相关数学模型，以制冷的年成本最小化为目标函数，运用数学规划法确定在已知冷量需求下的余热回收方案。建筑物对冷量的需求往往随季节变化，因此在模型中加入了多周期的描述，并考虑了不同周期下的系统操作优化。

2 问题表述

对于某些类型的建筑物，例如医院，往往一年四季都有冷量需求，被称为冷阱。将每个季节视为一个周期，周期的集合记为S = {T}。假设冷阱各周期的冷量需求量Qcd(T)已知，在离冷阱合适的距离范围内有一工厂，冷阱需要的冷量由该工厂提供。工厂内部有若干条余热流股i，集合为I = {i}，其入口温度θhin(i)、目标温度θhout(i)、热容流率Fh(i)、热负荷Qh(i)均已知。在不与制冷系统进行整合的情况下，这些热量由循环冷却水直接排放到环境中。通过以工艺软水为介质与这些余热流股换热，可以得到较高温度的热水，用以驱动LiBr 吸收式制冷机组的运转。完成制冷后，热水的温度降低。由于在工厂内部各余热流股与制冷机组间的距离往往不可忽略而且并不相等，因此需要用泵将热水输送至各余热流股的换热器，从而完成循环。在整个过程中，工厂的年费用 (total annual cost，TAC)可由式(1)表示。

式中：Qh0为所有余热流股原本需要的冷却负荷，Qr为回收得到的用于制冷的热量，Qcold为LiBr 机组的制冷量。故余热流股和制冷机组需要的总冷却公用工程负荷Qcu,total为

3 模型的建立

3.1 换热模型

为实现热量的梯级回收，采用分级换热网络超结构[10]描述余热回收。由于冷阱的供冷负荷随季节变化，因此有关变量需要引入周期集合 S 来表示，如图 1 所示。图中，θh(i,k,T)为周期T 下热物流i 在k 温区的热端温度，θwin(T)为周期T 下热水进行换热的入口温度，θwout(T)为周期T 下热水进行换热的出口温度，θw(k,T)为周期 T 下热水在 k 温区的热端温度，mw(T)为周期T 下热水的质量流率，q(i,k,T)为周期T 下热水在温区k 从热物流i 回收的热量，qcu(i,T)为周期T 下热物流i 需要的冷公用工程负荷。热水和余热流股需要满足进出口温度约束。

图1 分级超结构示意图Fig.1 The stage-wise superstructure

式中：zT(i,k,T)为决定周期T 下热水在温区k 与热物流 i 之间的换热器是否存在的二元变量，Δθ(i,k,T)和Δθ(i,k+1,T)分别为该换热器热端和冷端的传热温差。Δθmax为传热温差上限，取250 ℃。Δθmin为最小传热温差，取Δθmin= 10 ℃。

周期T 下热物流i 的冷却器的负荷为

式中：COP(T)为周期T 下制冷机的性能系数(coefficient of performance)，它与进入制冷机的热水温度(即θwout(T))有关。在制得的冷媒水温度为 7 ℃的情况下，不同的 θwout(T)所对应的 COP(T)由文献[9]提供。COP(T)与θwout(T)并非呈线性关系，因此为降低求解难度，不能直接对两者进行线性拟合，而应进行分段线性拟合，结果如式(26)所示。值得说明的是，在不同的冷媒水温度下，COP(T)与θwout(T)有着不同的对应关系，冷媒水温度可以作为优化变量。但目前将其作为变量纳入模型中会造成求解困难，因此在本文中，冷媒水温度设为定值7 ℃。

3.2 泵输送模型

如果热水在第k 个温区分流成多股与多条余热流股进行换热，则在稳定流动的状态下每条热水支路的压降均相等。然而，各余热流股 i 到制冷机组的距离不相等，使得单独输送时距离造成的压头损失Hf,branch(i)不相等，因此第k 温区的总压头损失Hf,total(k,T)应为各支路的Hf,branch(i)的最大值，这里采用不等式约束来表示，其上限值将由目标函数来限制，如式(27)所示。

3.3 目标函数计算

简化后的工厂年费用TAC’中只包括Cexchanger、Cpump,cap和Cpump,op三项，计算公式为

式中：AF 为年金因子(annualized factor)，α、β、γ 为换热器设备成本参数，a、b、c 为泵的设备成本参数，g 为重力加速度，η 为泵的效率，uce为单位电价，Prated为泵的额定功率，t(T)为各周期的时间。

4 案例研究与结果分析

4.1 案例数据

基于上述混合整数非线性规划(MINLP)模型，进行一项案例研究。工厂的余热流股的有关数据见表1，其中各流股的输送压头Hf,branch(i)作为参数直接给出，仅用来表示余热流股与制冷机组的距离远近。各周期的时间和各周期下医院的冷量需求量见表2。

表1 余热流股数据Table 1 The data of hot streams

表2 制冷量数据表Table 2 The data of cooling demand

各方程中涉及的有关参数取值如下[13]：cpw= 4.2 kJ⋅kg-1⋅℃-1，hw= 1.8 kW⋅m-2⋅℃-1，hh(i)= 2.0 kW⋅m-2⋅℃-1，g = 9.81 N⋅kg-1，AF= 0.264，α= 11 000 $，β= 150 $⋅m-2，γ= 1，a = 8 600 $，b = 7 310 $⋅kW-1，c = 0.2，uce= 0.1 $⋅(kW⋅h)-1，η = 0.7。

4.2 结果与讨论

本研究使用的求解工具为GAMS 软件的DICOPT 求解器，得到的各周期下的余热回收方案如图2～4所示，TAC’和各项费用见表3。

图2 春/秋季余热回收方案图Fig.2 Waste heat recovery scheme in spring/autumn

图3 冬季余热回收方案图Fig.3 Waste heat recovery scheme in winter

图4 夏季余热回收方案图Fig.4 Waste heat recovery scheme in summer

从 Cexchanger、Cpump,cap 和 Cpump,op 三项的计算方程中可以看出，影响 TAC’大小的核心变量是热水的出口温度θwout(T)。当θwout(T)增大时，LiBr 机组的 COP 值升高，需要回收的热量减小，且热水进出口温差增大，因此热水的质量流量降低，泵的费用减小。但与此同时换热器的传热温差减小，会导致所需要的换热面积增大，使换热器设备费用增大。反之，当θwout(T)减小时，泵的费用增大，换热器的费用减小。因此目标函数求解的核心就是换热器费用与泵费用之间的权衡。得到的结果表明，为使得 TAC’值最小，θwout(T)应处于 120 ～ 135 ℃。

由于各周期下医院的冷量需求差异明显，工厂的余热回收方案也随之发生变化。随着制冷量的减小，需要参与热量回收的余热流股数减少，未用于制冷的余热流股仍由原先设置的冷却器进行冷却。流股H1 的温位虽然不低，但其热量在4 个周期下均未得到回收，这是因为H1 距离制冷机组较远，若将热水输送至H1 处则需要泵提供较高的压头，使得泵的费用增大，因此在其他3 条流股的热量能够满足制冷需求的情况下，完全可以不用回收H1 的热量。由此可见，本案例中余热量大于制冷所需的热量，若制冷所需的热量增大，则可能会出现H1 的余热也被利用的情况。在制冷量最小的冬季，热水出口温度降低至120.6 ℃时，仅与距离制冷机组最近的H2 流股换热即可满足制冷机组热量需求，此时距离较远的流股H3 和H4 的换热器可以关停。这说明余热流股与制冷机组之间的距离对余热回收方案有着不可忽视的影响，应优先选择回收距制冷机组较近的余热流股的热量。

在制冷量最大的夏季，换热器需要提供较大的换热面积。在这种情况下，春季和秋季的热水可以取较高的出口温度，不仅使换热器面积的使用率较高，还能减小热水的质量流率，从而降低泵的费用。而在冬季，如果选择只回收流股 H2 的热量，则热水的出口温度不能取太高，这是因为传热温差的限制，当 θwout(T)增大时，H2 流股可以提供的热量会减小，直至热量的需求无法得到满足，从而要求再回收另一条流股的热量，造成泵的费用增大。因此，应在保证能用最少的流股满足热量需求的前提下提高θwout(T)，否则会使TAC’增大。

由于承担了医院的供冷任务，工厂需要增加制冷机组、换热器和泵的成本。但与此同时，出售冷量可以为工厂带来收入。可以求得工厂供冷的总投资TAC 为769 564 $。取供冷单价为120 $⋅(MW⋅h)-1，则供冷收入为2 589 120 $，大于TAC。由此可见，回收工业低温余热用于吸收式制冷，不仅可以实现能量的有效利用，还能为工厂带来可观的经济收益。

表3 经济性数据表Table 3 The economic data

5 结 论

本研究建立了工业余热用于区域供冷的多周期优化方法。通过对MINLP 模型的求解，可以得到各周期下的余热回收方案。结果表明，核心优化变量是得到的热水温度θwout(T)，它导致了泵的费用和换热器费用之间的权衡。当θwout(T)增大时，泵的费用降低，换热器费用增大。在研究的案例中，各周期下θwout(T)的取值均在120 ～ 135 ℃。在不同周期中，由于制冷量的变化，参与热量回收的工业流股数也会相应变化，冬季制冷量最小，因此只有一条工业余热物流参与余热制冷，其他余热物流由循环水系统冷却，反之夏季制冷量最大，则参与余热制冷的物流最多。在满足热量需求的前提下，热水应优先回收距制冷机组较近的余热流股的热量，从而减少泵的费用，使工厂的年费用TAC’最小。该方法可以优化出系统的最优设计和每个季节下系统的具体操作，经济性评估的结果表明，本案例中回收工业低温余热用于吸收式制冷的总成本为769 564 $，远小于收入2 589 120 $，经济效益显著。