彭景平，葛云征，陈凤云，刘蕾，吴浩宇，刘伟民*

( 1. 自然资源部第一海洋研究所 海洋工程环境研究中心，山东 青岛 266061)

1 引言

海洋温差能转换的实质是将储存在海水中的太阳能转换成为电能[1-2]，海洋覆盖地球约71%的面积，是巨大太阳能接收器。海洋是地球上巨大的可再生能源载体，而温差能是诸多海洋能中储量最大的可再生能源[3]。随着常规化石能源的紧缺和造成的环境污染问题日益严重，而海洋温差能又具备受季节及昼夜的影响小、储量大和稳定性好的优势，海洋温差能开发利用倍受国际研究者的青睐[4]。但是由于表层海水与深层海水仅有约为15～25℃的温差，可利用理论循环热效率通常为3%～5%[5-6]。提高海洋温差能利用热效率成为国内外研究人员的研究重点和方向。

国内外学者对提升海洋温差能发电系统利用效率方面做了一些研究，包括工质选取和热力循环形式等方面。Chen 等[7]对朗肯循环进行了理论和试验研究，分析了R717、R22、R125、R134a、R152a 和R245fa6中工质在不同透平进口压力下的循环热效率，研究结果表明R717 是海洋温差能最适合的工质。黄靖伦等[8]采用R245fa/R123 非共沸混合工质对有机朗肯循环余热回收系统的热力性能进行了研究，分析了分别采用R245fa、R123 单工质和R245fa/R123 混合工质循环系统的热力性能，结果显示混合工质可提高系统的㶲效率。Aydin 等[9]开展了采用太阳能再热提高朗肯循环效率的研究，研究结果表明在透平进口设置太阳能再热可以使循环效率从1.9%提高到3%。Kusuda等[10]对双级朗肯循环进行了研究，研究结果显示双级朗肯循环系统可以减小换热器内的不可逆损失，有利于提高循环系统的输出功。Kalina[11]在1984 年提出了一种新的热力循环系统，该循环系统热效率是朗肯循环的1.6～1.9 倍。Uehara 等[12]提出了一种新型高效海洋温差能热力循环，采用氨-水非共沸混合工质，采用了两个透平，并采用了中间抽气回热的方法，系统循环热效率高于Kalina 循环[13-14]。Yuan 等[15-16]通过在循环系统中安装引射器对Uehara 循环加以改进，采用氨水再热-引射吸收的方法提高循环效率，研究结果表明温差为24℃时，循环效率达到5.27%，循环热效率高于Uehara 循环。Yoon等[17]提出了一种带有两级透平和回热器的单工质海洋热能转换热力循环，通过对透平效率、换热器压降、透平出口压力等参数分析系统性能，研究结果表明系统热效率略高于Uehara 循环。

目前，有关蒸发压力、工质质量分数等因素对海洋温差能热力循环效率的研究较少，本文提出一种新的采用非共沸工质海洋温差能回热循环系统，选取蒸发压力、工质质量分数、蒸发冷凝温度参数对循环系统进行热力分析。

2 热力循环系统介绍

本文提出的高效海洋温差能回热循环系统由一个主循环和两个热回收支路构成，选取R134a/R123混合物作为系统工质，其原理图和T-S-w图如图1、图2 所示。循环系统由蒸发器、气液分离器、两级透平发电机、工质泵、回热器、吸收器和冷凝器组成。其工作原理为：蒸发器中表层热海水将混合工质加热成为气液两相混合溶液，气液两相混合工质在气液分离器中分离为气相和液相，气相工质蒸气进入透平做功，液相工质溶液进入热回收支路。工质蒸气在第一级透平做功后的乏气，有少部分抽出进入抽气热回收支路，剩余的乏气进入第二级透平做功。在吸收器中两条热回收支路的工质与第二级透平排出的乏气混合进入冷凝器，深层冷海水将混合工质冷却为液态后，由工质泵输送到回热支路回热器进行预热，然后进入蒸发器进入下一次循环。两条热回收支路通过回热器1、2、3 回收利用分离器分离出的液态工质和双级透平中间抽气的热量，可有效提高循环系统的热效率。

图1 新型海洋温差能回热循环系统原理示意图Fig. 1 Schematic diagram of new ocean thermal energy conversion regenerative cycle system

图2 新型海洋温差能回热循环系统T-S-w 图Fig. 2 T-S-w diagram of new ocean thermal energy conversion regenerative cycle

3 仿真模型建立

3.1 数值模型

在分析提出的循环热力过程中的各节点位置如图1 所示。在数值计算过程中，为了在理论研究中可以明确循环热力过程中各基本热力参数关系，定量分析提出的热力循环系统，拟进行以下简化假设：

（1）循环系统中各分离、交汇处的蒸气压力相同；

（2）冷凝器出口的工质为饱和液体；

（3）分离器分离出的液体和气体均为饱和状态；

（4）两级透平内的气态工质和工质泵内的液态工质均为等熵过程；

（5）换热器的窄点端差选取2℃。

根据质量守恒和能量守恒方程，系统中主要设备的数值模型建立如下：

蒸发器和冷凝器的能量守恒方程可以表示为

式中，Qe为蒸发器换热量，单位：kW；cw为温海水比热，单位：kJ/(kg·℃)；mw为温海水质量流量，单位：kg/s；Twi为蒸发器进口温海水温度， 单位：℃；Two为蒸发器出口温海水温度， 单位：℃；m1为蒸发器出口（1 点）工质质量流量，单位：kg/s；h1为蒸发器出口（1 点）工质焓值，单位：kJ/kg；h11为蒸发器进口（11 点）工质焓值，单位：kJ/kg；Qc为冷凝器换热量，单位：kW；cc为冷海水比热，单位：kJ/(kg· ℃)；mc为冷海水质量流量，单位：kg/s；Tci为冷凝器进口温海水温度， 单位：℃；Tco为冷凝器出口冷海水温度， 单位：℃；m5为冷凝器出口（5 点）工质质量流量，单位：kg/s；h5为冷凝器出口（5 点）工质焓值，单位：kJ/kg；h6为冷凝器进口（6 点）工质焓值，单位：kJ/kg。

系统中3 个回热器的能量守恒方程可以表示为式中，m8为回热器1 进口（8 点）混合工质溶液的质量流量，单位：kg/s；h8为回热器1 进口（8 点）混合工质溶液的焓值，单位：kJ/kg；h9为回热器1 出口（9 点）混合工质溶液的焓值，单位：kJ/kg；m13为回热器1 进口（13 点）低浓度工质溶液的质量流量，单位：kg/s；h13为回热器1 进口（13 点）低浓度工质溶液的焓值，单位：kJ/kg；h14为回热器1 出口（14 点）低浓度工质溶液的焓值，单位：kJ/kg；m9为回热器2 进口（回热器1 出口，9 点）混合工质溶液的质量流量，单位：kg/s；h10为回热器2 出口（10 点）混合工质的焓值，单位：kJ/kg；h16为回热器2 出口（16 点）抽气工质的焓值，单位：kJ/kg；m10为回热器3 进口（10 点）混合工质溶液的质量流量，单位：kg/s；h11为回热器3 出口（11 点）混合工质溶液的焓值，单位：kJ/kg；h12为回热器3 进口（12 点）低浓度工质溶液的焓值，单位：kJ/kg。

气液分离器的质量守恒方程和能量守恒方程可表示为

式中，m2为一级透平进口（2 点）气相工质质量流量，单位：kg/s；m12为分离器液相出口（12 点）工质质量流量， 单位：kg/s；h2为一级透平进口（2 点）气相工质焓值，单位：kJ/kg；h12为分离器液相出口（12 点）工质焓值， 单位：kJ/kg。

吸收器的能量守恒方程可表示为

式中，m5为冷凝器进口混合工质的质量流量，单位：kg/s；h5为冷凝器进口混合工质的焓值，单位：kJ/kg；m4为二级透平出口高浓度R134a 蒸气的质量流量，单位：kg/s；h4为二级透平出口高浓度R134a 蒸气的焓值，单位：kJ/kg；m15为膨胀阀1 出口低浓度R134a 溶液的质量流量，kg/s；h15为膨胀阀1 出口低浓度R134a 溶液的焓值，单位：kJ/kg；m17为膨胀阀2 出口低浓度R134a 溶液的质量流量，单位：kg/s；h17为膨胀阀2 出口低浓度R134a 溶液的焓值，单位：kJ/kg。

双级透平：在提出的循环系统中采用的透平形式是双级透平，在一级透平出口处进行抽气，抽出进入抽气热回收支路。中间抽气点的压力按下式选取

透平内工质等熵膨胀做功，透平的能量守恒方程可表示为

式中，Wt1为一级透平输出功，单位：kW；h3为透平中间（3 点）抽气工质焓值，单位：kJ/kg；Wt2为二级透平输出功，单位：kW；m3为透平中间（3 点）抽气工质质量流量，单位：kg/s；h4为二级透平出口（4 点）工质焓值，单位：kJ/kg。

工质泵的能量守恒方程可表示为

循环热效率和系统效率可以表示为

3.2 计算流程

依据可利用海洋温差能的海域实际海水温度情况，表层温海水温度选取26℃，深层冷海水温度选取5℃，新型海洋温差能回热循环系统的理论分析初始条件组合如表1 所示。循环系统中各节点的状态参数及工质的热力学性质使用自主研发的热物性计算软件进行计算。依据质量守恒和能量守恒定律，循环系统计算流程如图3 所示。

图3 新型海洋温差能回热循环系统热力计算流程图Fig. 3 Flowchart for thermodynamic calculation of new ocean thermal energy conversion regenerative cycle

表1 循环系统的理论分析初始条件组合Table 1 Initial conditions for theoretical analysis of circulatory system

4 结果讨论

在海洋温差能利用系统中影响系统循环热效率的因素主要有工质质量分数、蒸发压力、蒸发和冷凝温度。本文通过数值计算研究了工质质量分数和蒸发压力参数对循环热效率和系统净输出功的影响。在循环热效率的计算过程中采用单位工质进行计算，在系统净输出功率计算过程中采用装机容量为10 kW的模型进行量算。

4.1 工质质量分数对循环热效率和净输出功率的影响

在给定的温差能循环系统初始条件：表层温海水温度为26℃，冷海水温度为5℃，蒸发器和冷凝器窄点端差均为2℃，蒸发压力为0.6 MPa。循环热效率随着混合工质质量分数的增大先增大后减小，如图4所示。循环热效率在混合工质质量分数小于0.9 时，循环热效率随工质质量分数的增大快速升高，在混合工质质量分数为0.9～0.93 之间时，系统热效率变化幅度不大，并在混合工质质量分数为0.91 时，系统热效率取得极大值5.28%，当混合工质质量分数超过0.93时，随着工质质量分数的增加，循环热效率迅速减小。出现上述变化规律的原因是：在蒸发温度不变的情况下，混合工质的气化率随着工质质量分数增加而增大，透平输出功率随之增大。同时，随着混合工质气化率的增加，蒸发器内的换热量增大。蒸发器内换热量的增大比透平输出功率增大的趋势更快，如图5所示。因而出现了循环热效率随混合工质质量分数的增大先增大后减小的变化规律。

图4 循环热效率和净输出功率随工质质量分数的变化Fig. 4 Variation of cycle thermal efficiency and net output with the change of working fluid mass fraction

图5 换热量和透平输出功率随工质质量分数的变化Fig. 5 Variation of heat transfer rate and turbine output with the change of working fluid mass fraction

随着混合工质质量分数的增大，系统净输出功率先增大后减小，在工质质量分数为0.96 时取得极大值3.83 kW，如图4 所示。在装机容量不变的情况下，蒸发器冷凝器内的换热量基本不变。随着混合工质气化率的增加，系统内混合工质的质量流量逐渐减小，随之减小的是工质泵的功率，系统的净输出功率增大。当工质质量分数继续增大，分离器内液相溶液过少导致液相热回收支路失去了回热功能，循环热效率减小迅速导致净输出功率减小。

4.2 蒸发压力对循环热效率和净输出功率的影响

在给定温差能循环系统初始条件，表层温海水温度为26℃，深层冷海水温度为5℃，混合工质质量分数为0.9。循环热效率随蒸发压力的增大先增大后减小，蒸发压力由0.55 MPa 增加到0.595 MPa 时，循环热效率从4.60%升高到5.26%，在蒸发进口压力继续增加到0.604 MPa 时，循环热效率减小到5.00%，如图6所示。在蒸发温度和工质质量分数不变的情况下，随着蒸发压力的增加，混合工质的气化率减小导致透平输出功率减小。此时，随换热器内混合工质气化率减少，换热器内换热量减小，相比于蒸发器换热量透平输出功率的变化幅度更大，如图7 所示。因此出现了循环热效率随蒸发压力的增大先增大后减小的变化规律。

图6 循环热效率和净输出功率随蒸发压力的变化Fig. 6 Variation of cycle thermal efficiency and net output with the change of evaporation pressure

图7 透平做功和蒸发器换热量随蒸发压力的变化Fig. 7 Variation of turbine output and heat transfer rate with the change of evaporation pressure

系统净输出功率随着蒸发压力的增加先缓慢增大后快速减小，在蒸发压力为0.58 MPa 时取得极大值3.57 kW，如图6 所示。随着蒸发压力的升高蒸发器内工质气化率逐渐减小，单位工质的透平输出功率降低，在装机功率不变的情况下，工质质量流量随蒸发压力的升高而增加，工质泵耗功随之增大。在循环热效率的影响下，随着蒸发压力的升高海水泵的耗功先减小后增大。在海水泵耗功、透平输出功率和工质泵耗功耦合作用下，系统净输出功率随蒸发压力的增加呈现出先增大后减小的变化规律。

4.3 验证对比

通过在一定的运行工况下与Uehara 循环和Yoon循环的计算结果对比，以验证文中提出的新型海洋温差能回热循环具有较高的循环热效率。文献[12]中运行工况为：工质采用氨的质量分数为0.9 的氨水混合工质，蒸发压力为0.8 MPa，蒸发温度和冷凝温度分别为299.15 K 和279.15 K。在上述工况下，Uehara 循环和提出新型海洋温差能回热循环的循环热效率分别为5.3%和5.6%，新循环的热效率比Uehara 循环高了5.66%，提出的新循环系统各节点计算参数见表2，节点编号见图1。在文献[17]所述的运行工况下，将提出的新型海洋温差能回热循环与Uehara 循环、Yoon循环计算结果进行对比，如表3 所示。新热力循环的系统效率为2.532%，高于Uehara 循环和Yoon 循环的2.379%和2.401%。

表2 文献[12]所述工况下提出的循环各节点参数Table 2 Parameters of the proposed cycle under the conditions described in reference [12]

表3 文献[17]所述工况下Uehara 循环、Yoon 循环和新循环的参数比较Table 3 Comparison of parameters of Uehara cycle, Yoon cycle and proposed cycle under the conditions described in reference [17]

5 结论

本文依据热力学定律，选取混合工质质量分数和蒸发压力性能参数对提出的新型海洋温差能回热循环进行热力分析研究。得到了提出的新热力循环系统循环热效率和净输出功率的变化规律，并与其他热力循环形式进行了对比，得到了如下结论。

（1）在给定的运行工况下，提出的新型海洋温差能回热循环系统净输出功率和热效率均随混合工质质量分数的增加先增大后减小，随蒸发压力的增加先增大后减小。

（2）新型海洋温差能回热循环的热效率最佳值点和净输出功率最大值点并不是同一个参数值点。以工质质量分数为变量时，循环热效率在工质质量分数为0.91 时取得极值5.28%，净输出功率在工质质量分数为0.96 时取得极值3.83 kW；以蒸发压力为变量时，工质质量分数为0.9 时，循环热效率在蒸发压力为0.595 MPa 时取得极值5.26%，净输出功率在蒸发压力为0.58 MPa 时取得极值3.57 kW。

（3）在相同的温冷海水和工质质量分数等运行参数下，与Uehara 循环和Yoon 循环进行对比，提出的新型海洋温差能回热循环性能最佳。在文献[12]的运行参数下，提出的新循环热效率比Uehara 循环高5.66%；在文献[17]的运行参数下，提出的新循环系统效率为2.532%高于Uehara 循环的2.379%和Yoon 循环的2.401%。