过冷法海水冰浆制备系统能耗特性分析*

宋文吉，黎福超3，陈明彪，冯自平,2

过冷法海水冰浆制备系统能耗特性分析*

宋文吉1,2,†，黎福超1,2,3，陈明彪1，冯自平1,2

（1. 中国科学院广州能源研究所，广州 510640；2. 中国科学技术大学 能源科学与技术学院，广州 510640；3. 华南理工大学，广州 510641）

海水制冰浆技术可应用于海岛及海上风电的负荷转移和削峰以及冷冻保存海产品等，具有广阔应用前景。此外，海水制冰的原料近乎免费，并且利用海水作为冷却介质可以降低冷凝温度，进而提高系统能效，具有明显优势。通过建立双回路单级压缩循环制冰系统模型，分析利用过冷法制取海水冰浆的经济性和可行性，为海水制冰浆技术的实际应用提供理论指导。结果显示，该系统在计算条件下冬季循环性能系数平均值为4.76，夏季平均值为3.61，在蒸发器内较宽的海水流速范围内可维持较长时间正常运行。系统采用板式水冷冷凝器，并与空气源机组进行定量对比，结果表明循环机组采用水冷冷凝器在各指定工况点下的效益均高于采用空冷器，在设计工况下水冷机组单位投资蓄冷量为2.85 kW∙h/元，对应的空冷机组仅为1.37 kW∙h/元。

冰浆；海水；过冷法；双回路单级压缩制冷循环；板式换热器；效率

0 引 言

冰浆因具有较高的能量储存密度（334 kJ/kg）和良好的换热特性，是极佳的冷量存储介质。利用海水制冰浆具有独特优势，一方面海水制冰的原料丰富而且几乎免费，另一方面利用海水作为冷却介质可以降低冷凝温度，可以达到提高系统能效的目的。此外，储冷的成本约为储电成本的10%[1]，经济优势明显。

海水制冰浆应用前景广阔。海上风电与海岛产电不易储存，可以用多余的电量制取冰浆，将电能以冷量的形式存储起来。例如夜晚电量需求小，此时利用海水制冰浆存储起来，白天用电高峰期时又可以融化冰浆快速释放冷量，供给建筑物的空调装置。除此之外，冰浆还可用于冷冻保存海产品，其保鲜效果比传统的冰块冷冻更好，是一种良好的保存介质[2]。

国内外已经有不少学者对海水制冰浆相应技术进行了研究。王振等[3]就超声波对海水流化冰制取过程中的作用效果和作用机理进行了研究，发现在33.0 kHz的实验工况下海水过冷度可以被完全消除，制冰速度得到明显提升，此外，相对于无超声波工况，冰晶粒径更小、更绵柔，不易造成冰堵。CHEN等[4]评估了使用过冷水动态制冰法来淡化海水的可行性和经济性，结果表明该方法比现有的间接接触冷冻法淡化海水更加快速有效，且由于动态冰浆换热性能更好等原因，理论能耗仅为间接接触冷冻法的58%。袁久峰等[5]针对目前海水流化冰机长时间运行出现冰堵和能耗升高的问题，选型设计了一种海水流化冰机，可实现连续稳定的制取海水流化冰，但未进行经济性评估。HEKMATSYAR等[6]对涂有聚四氟乙烯的刮刀式海水冰浆发生器的生产速度进行了研究，发现在海水盐度较高的情况下提高泵或刮刀的转速都会降低冰浆的生产速度，且能耗升高。KONG等[7]对使用钛合金双管蒸发器的制冰系统进行实验研究，结果显示在名义条件下制冰机的性能系数（coefficient of performance, COP）为1.66，而随着海水进口温度的降低，冰浆含冰率和冰浆产率都大大增加。综上，不少文献对海水制取冰浆技术的实际应用做出了理论指导，但目前国内外研究大多关注制冷蒸发侧海水制成冰浆的循环过程，而对于冷凝侧将海水同时作为冷却介质以及相关的系统优化研究仍较少。

为此，本文利用海水制冰时既利用海水作为原材料且同时作为冷却介质的优点，采用直接蒸发式的过冷法冰浆制备系统，建立海水制冰的系统模型，分析其在不同运行条件下的能耗特性以及经济性，以期为海水制取冰浆技术的实际应用提供理论指导。

1 计算模型与数据处理方法

1.1 计算模型

拟设计的海水制冰浆装置如图1所示，采用单级压缩制冷双回路循环，制冷剂为R22，循环设计了3 ℃的过热度和2 ℃的过冷度，相较于一般的三回路循环效率更高。冷凝器和蒸发器均采用板式换热器，冷凝器板片数为25片，换热面积12 m2。蒸发器板片数为57片，换热面积28.6 m2。板片参数见表1。

图1 海水制冰浆装置简图

表1 水平波纹形式板式换热器板片参数

海水制冰浆装置设置了蓄冰池，池内维持在凝固点温度，即−1.83 ℃。为避免池内的冰晶在循环过程中进入蒸发器内而引发冰堵，海水在进入蒸发器之前经加热器加热，使其温度升高到0.1 ℃，加热可通过与高温的海水换热实现，故不计算加热电耗。在蒸发器内，海水会获得一定的过冷度，过冷度过低会降低产冰率，过高易使蒸发器内形成冰堵，影响装置正常运作。为此将过冷度控制为1.5 ℃[8]，过冷的海水在促晶器中结晶，进入蓄冰池，未形成冰晶的海水过滤后再次进入循环。由于循环过程中池内海水盐度会上升，需要及时更换海水，更换海水损失的冷量可通过回热器回收，这里不做讨论。

冷凝器侧也使用海水水冷，假设海水经过冷凝器后温升为5 ℃。

选取夏季（6、7、8月）以及冬季（12、1、2月）东海海水数据作为计算数据，如表2所示。

表2 东海夏季和冬季海水、空气数据

1.2 系统主要性能计算方法

板式换热器海水侧换热系数采用准则关系式描述[9]：

板式换热器制冷剂沸腾换热系数[10]：

板式换热器制冷剂冷凝换热系数[11]：

总的传热系数由下式计算：

沸腾换热量可通过水侧计算得到：

冷凝器散热量可由下式计算：

换热器所需换热面积计算式：

制冰量可由下式计算：

空气冷却器制冷剂水平微内肋管冷凝换热系数[12]：

空气冷却器空气侧换热系数[9]：

㶲效率计算式：

蒸发器海水侧出口段壁温计算式：

系统循环性能系数计算式：

2 计算结果分析

2.1 制冷循环效率分析

在保持过冷度为−1.5 ℃的条件下，改变蒸发器内海水流速以控制制冷量变化。低海水流速时换热系数低，冰浆产量少，不能充分利用换热设施；增大海水流速，在保持相同过冷度的情况下，则会增大产生冰堵的可能性[8]。因此计算中控制海水流速在0.05 ～ 0.6 m/s之间变化。计算出的系统整体循环COP和循环㶲效率、制冰量如图2所示。

冬季循环COP平均值为4.76，夏季平均值为3.61。低海水流速下COP和循环㶲效率最高，COP可达5.6，㶲效率可达0.5；且冬季COP比夏季高1以上，这是由于冬季海水温度较低，制冷剂冷凝温度下降，经济性提高。虽然夏季冷凝温度较冬季高，单位制冷剂压缩耗功也更高，但夏季循环㶲效率反而比冬季高0.1以上，主要原因是夏季环境温度也较冬季升高，且是主要的影响因素，此时制造同样的冷量放出的冷量㶲更大，冷量的品质更高。高海水流速时制冰量较低流速时成倍增长，制冰浆的效率大大提高，但此时经济性下降。相应的每小时冰浆蓄冷量与系统耗电量关系曲线如图3所示，很显然斜率越大代表经济性越好。实际生产时应该综合考虑，选择适中的海水流速，保证制冰量的同时提高经济性，减少冰堵的发生。

图3 夏、冬季蓄冷量与耗电量关系曲线

海水出口段的壁温如图4所示。出口段（对应制冷剂干度为0.99的部分）的壁温最低，实际生产过程中应该控制这部分的壁温，壁温过低会直接导致换热器内冰晶的生成，造成频繁冰堵进而影响设施正常工作。在本例中，认为壁面温度高于−5 ℃即可[14]，因此海水流速在0.05 ～ 0.6 m/s范围内变化设备可保持较长时间正常工作，最高流速可达1.1 m/s。若要进一步提高海水流速并降低冰堵的可能性，可在海水侧换热表面涂抹抑制冰晶生成的涂层[15]。

图4 出口段壁温随海水流速变化规律

2.2 与空气冷却器对比

为进一步探究该采用水冷板式冷凝器的海水制冰浆循环装置经济性，将其与等换热面积的空冷器进行定量对比。空冷器换热管参数如表3所示。

表3 空气冷却器水平式低翅片换热管参数

保持蒸发器侧换热器、参数不变，将冷凝器侧换为空气冷却器，且保持换热面积与板式冷凝器一致，为12 m2。此外假设空气经过空冷器后温升为10 ℃。计算制冷量变化时的系统循环COP和循环㶲效率如图5所示。

由图5可知，在每个工况点下水冷的经济性都要高于空冷。夏季水冷系统循环COP均值与风冷差值为1.30，高出56.3%，冬季时差值为1.15，高出31.9%。在低负荷运行时，空冷与水冷的COP值和循环㶲效率较为接近，但负荷越高，空冷的这两个性能指标与水冷相差就越大，这一点在冬季环境下最为明显，最高负荷时COP差值达到了2.2。这是由于水冷板式换热器整体换热系数要比空气冷却器高出许多，在高负荷运转时该趋势更加明显，经济性更高。而且在高温或高负荷下，使用板式水冷冷凝器经济性比使用空冷器高得多，特别是在夏季高负荷时，COP和㶲效率的差值分别达到了2.1、0.21，水冷经济性远高于空冷。与之相对应的每小时蓄冷量与耗电量关系曲线如图6所示，图中能更直观地看出高负荷下空冷耗电量远大于水冷，而低负荷时两者相差不大。

图6 夏冬季水冷、空冷蓄冷量与耗电量关系曲线

分析采用水冷冷凝器和空冷冷凝器时整体机组的经济性，考虑投资成本和运营成本，通过文献以及调研可知：①海上风电电价按0.8元/(kW∙h)计算；②机组设备费按照制冷量计算，空冷机组为0.8元/W，水冷机组按照0.9元/W计算[16]；③安装费用按照设备费的10%计算。

以200 kW制冷量作为机组的额定制冷量，得出空冷与水冷机组的初投资成本如表4所示。年运营成本是将夏季与冬季的功耗特性综合平均而计算得出，其中包括水泵、压缩机及空冷器风机电耗，不同制冷功率下运行年费用如表5所示。

表4 两种机组初投资成本明细

表5 不同制冷量条件下两种机组的运营成本

由表4和表5可知，水冷机组初投资成本稍高于空冷机组，但空冷机组年运营成本要高于水冷机组，且在制冷功率较大时特别明显。综合投资成本与运营成本，首年单位投资蓄冷量数据如表6所示。

表6 不同制冷量条件下两种机组的首年单位投资蓄冷量

水冷机组在设计工况即200 kW制冷功率下具有最高的单位投资蓄冷量，为2.85 kW∙h/元，对应的空冷机组仅为1.37 kW∙h/元。且空冷机组的单位投资蓄冷量随着制冷功率上升呈先增大再减小的趋势，因此在设计工况下达不到最佳经济性，再次说明高负荷采取水冷机组具有更好的经济性。

为探究海水盐度变化时装置运行的功耗情况，在夏季和冬季的气温、水温条件下，令蒸发器侧海水流速为0.25 m/s，并保持过冷度为1.5 ℃，海水盐度在5‰ ～ 50‰之间，因此海水的出口温度在−4.55 ～−1.78 ℃之间变化。计算出各盐度下的系统循环COP和功耗如图7所示，与之相对应的每小时蓄冷量与耗电量关系曲线如图8所示。

随着盐度的增大，蒸发器侧海水出口温度降低，蒸发温度也随着降低，因此COP曲线随着下降。但采用水冷的循环COP稳定地大于空冷COP，随着盐度增大，COP差值稍有增大。盐度增大时由于制冷量的增大，蒸发温度的下降，压缩功持续增加，且盐度越大增速越快，在高盐度时采用空冷循环的压缩功耗明显大于采用水冷循环时的功耗，该趋势可从图8上直观看出。

图8 夏冬季变盐度条件下水冷、空冷蓄冷量与耗电量关系曲线

为比较板式水冷冷凝器和空冷器在不同气温情况下的实用性和经济性，令蒸发温度为−5℃，冷凝器侧入口水温与入口气温相等，温度设置在10 ～ 35℃之间，换热面积相同，计算出不同环境温度下系统循环COP和循环㶲效率如图9所示。

图9 水冷、空冷的系统COP值、㶲效率随冷凝器入口温度变化规律

显而易见，随着冷凝器入口水温、气温的升高，水冷与空冷的循环COP都呈现出下降的趋势，但是水冷循环的COP总是大于空冷，并且温度越低，二者之间的差值越大，在10 ℃时差值达到了2.1。即使在35 ℃条件下，水冷COP仍比空冷大82%，差值为1.18。此外，水冷循环㶲效率也总是大于空冷的㶲效率，差值恒保持在0.1以上。与之相对应的每小时蓄冷量与耗电量关系曲线如图10所示，与图9中COP不同的是，在高温时制冷负荷更高的情况下，空冷耗电量远大于水冷耗电量。

图10 变冷凝器入口温度条件下水冷、空冷蓄冷量与耗电量关系曲线

与COP变化趋势不同的是，随着温度的升高循环㶲效率先上升后下降，该结果由式（11）分析可知。由于蒸发温度不变，影响循环㶲效率的因素可归结为环境温度和COP，前半部分㶲效率上升是环境温度上升的结果，此时环境温度是主要影响因素，后半段温度影响减弱，COP下降并且其影响占据主导地位。

3 结 论

采用海水水冷的海水制冰浆装置具有较高的经济性和可行性，整个系统的循环COP最高可达5.6，系统循环㶲效率最高可达0.5，冬季循环COP平均为4.76，夏季平均为3.61。在较宽的海水流速范围内可维持较长时间正常运行。

采用水冷冷凝器整体传热系数高，可降低制冷剂冷凝温度，使其在各个指定工况点下效益都要高于采用空冷冷凝器，夏季时水冷平均COP比空冷高56.3%，平均㶲效率高37.4%，冬季时水冷平均COP比空冷高31.9%，水冷平均㶲效率比空冷高24.2%。蓄冷量与耗电量关系曲线表明，低制冷负荷下两者耗电量相差不大，负荷升高时空冷耗电量远大于水冷耗电量。经济性分析也表明在设计工况下水冷机组单位投资蓄冷量为2.85 kW∙h/元，对应的空冷机组仅为1.37 kW∙h/元。在水温、气温相等的情况下，水冷的经济性也总是高于空冷，且温度越低水冷越显示出优越性。但海水具有一定的腐蚀性，在使用时需要注意进行防腐处理。

努塞尔数

雷诺数

普朗特数

常数，流体被加热时为0.4，被冷却时为0.3

换热系数，W/(m2·℃)

放大系数，取1

l液相导热系数，W/(m·℃)

h水力直径，m

eq等效雷诺数

质量流率，kg/(m2·s)

eq等效质量流率，kg/(m2·s)

l、g液、汽相动力黏度，Pa·s

m平均干度

l、g液、汽相密度，kg/m3

h重力驱动冷凝换热系数，W/(m2·℃)

shear剪切力驱动冷凝换热系数，W/(m2·℃)

s蒸汽饱和温度，℃

w壁面温度，℃

重力加速度，m/s2

汽化潜热，J/kg

板片长度，m

l液相对流换热系数，W/(m2·℃)

总传热系数，W/(m2·℃)

I、II两侧换热面换热系数，W/(m2·℃)

板片厚度，m

污垢热阻，取1.7 × 10−5m2·℃/W

0制冷量，W

w水流量，kg/s

c,w水定压比热容，J/(kg·℃)

w,i、w,o进、出口水温度，℃

k冷凝散热量，W

f制冷剂流量，kg/s

c,s压缩机绝热效率，取0.8

1压缩机入口的制冷剂比焓值，kJ/kg

2理想状态下压缩机出口的制冷剂比焓值，kJ/kg

4蒸发器入口的制冷剂比焓值，kJ/kg

换热面积，m2

换热量，W

Δlm对数换热温差，℃

w,s凝固点温度，℃

w,f给水温度，℃

相变潜热，kJ/kg

热流密度，W/m2

i管内径，m

NF迎风面风速，m/s

校正系数，取0.946

0环境温度，℃

1蒸发温度，℃

0单位制冷剂制冷量，kJ/kg

net单位制冷剂循环实际耗功量，kJ/kg

el压缩机电效率，取值参考文献 [17]

m出口段热流密度，W/m2

m出口段水侧换热系数，W/(m2·℃)

m出口段壁温，℃

m出口段局部总传热系数，W/(m2·℃)

l,eq两相流中仅液相部分的传热系数，W/(m2·℃)

干度，取0.99

[1] 安建才, 司爱国, 高玉国. 冰浆制取技术及其应用的研究进展[J]. 电力与能源进展, 2021, 9(2): 11-12. DOI: 10.12677/aepe.2021.92002.

[2] LOSADA V, BARROS-VELÁZQUEZ J, AUBOURG S P. Rancidity development in frozen pelagic fish: influence of slurry ice as preliminary chilling treatment[J]. LWT - food science and technology, 2007, 40(6): 991-999. DOI:10.1016/j.lwt.2006.05.011.

[3] 王振, 李保国, 罗权权. 超声波对海水流化冰制取过程的影响[J]. 包装工程, 2020, 41(17): 65-70. DOI: 10.19554/j.cnki.1001-3563.2020.17.009.

[4] CHEN D, ZHANG C S, RONG H W, et al. Experimental study on seawater desalination through supercooled water dynamic ice making[J]. Desalination, 2020, 476: 114233. DOI: 10.1016/j.desal.2019.114233.

[5] 袁久峰, 贺莹炜, 殷娟. 渔船用海水流化冰机设计探究[J]. 内燃机与配件, 2021(12): 211-212. DOI: 10.3969/ j.issn.1674-957X.2021.12.099.

[6] HEKMATSYAR D F, RAYHAN F A, PAMITRAN A S. Optimization of ice slurry generator with scraper blade coated by teflon for fishery industry[J]. AIP conference proceedings, 2019, 2062: 020030. DOI:10.1063/1. 5086577.

[7] KONG K S, KIM J D. Haracteristic analysis of slurry ice maker with a double tube according to sea water temperature[J]. Journal of power system engineering, 2019, 23(6): 88-94. DOI: 10.9726/kspse.2019.23.6.088.

[8] BÉDÉCARRATS J P, DAVID T, CASTAING-LASVIGNOTTES J. Ice slurry production using supercooling phenomenon[J]. International journal of refrigeration, 2010, 33(1): 196-204. DOI:10.1016/j. ijrefrig.2009.08.012.

[9] 史美中, 王中铮. 热交换器原理与设计[M]. 6版. 南京: 东南大学出版社, 2018.

[10] LONGO G A, MANCIN S, RIGHETTI G, et al. A new model for refrigerant boiling inside brazed plate heat exchangers (BPHEs)[J]. International journal of heat and mass transfer, 2015, 91: 144-149. DOI: 10.1016/j. ijheatmasstransfer.2015.07.078.

[11] MANCIN S, DEL COL D, ROSSETTO L. Condensation of superheated vapour of R410A and R407C inside plate heat exchangers: experimental results and simulation procedure[J]. International journal of refrigeration, 2012, 35(7): 2003-2013. DOI: 10.1016/j.ijrefrig.2012.06.001.

[12] 石文星, 田长青, 王宝龙.空气调节用制冷技术[M]. 5版. 北京: 中国建筑工业出版社, 2016.

[13] ARIMA H, KIM J H, OKAMOTO A, et al. Local boiling heat transfer characteristics of ammonia in a vertical plate evaporator[J]. International journal of refrigeration, 2010, 33(2): 359-370. DOI:10.1016/j.ijrefrig.2009.09.017.

[14] 吴锐. 过冷法动态冰浆制取的分析及实验[D]. 武汉: 华中科技大学, 2006.

[15] KAUFFELD M, GUND S. Ice slurry-history, current technologies and future developments[J]. International journal of refrigeration, 2019, 99: 264-271. DOI:10.1016/j.ijrefrig.2019.01.010.

[16] 陈明彪, 宋文吉, 王瑛滢, 等. 跨季节蓄冷特性与系统设计优化[J]. 太阳能学报, 2022, 43(6): 1-7. DOI:10.19912/j.0254-0096.tynxb.2020-1043.

[17] 马一太, 刘忠彦, 李敏霞. 容积式制冷压缩机电效率分析[J]. 制冷学报, 2013, 34(3): 1-7. DOI: 10.3969/j. issn.0253-4339.2013.03.001.

Energy Consumption Characteristics Analysis of Seawater Ice Slurry Generation System Using Supercooling Method

SONG Wenji1,2,†, LI Fuchao1,2,3, CHEN Mingbiao1, FENG Ziping1,2

(1. Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences, Guangzhou 510640, China; 2. School of Energy Science and Engineering, University of Science and Technology of China, Guangzhou 510640, China; 3. South China University of Technology, Guangzhou 510641, China)

Technology for producing ice slurry from seawater can be applied to the load shifting and peak shaving for island and offshore wind power, as well as to freeze and preserve seafood, which has broad application prospects. In addition, the source material for producing ice slurry from seawater is virtually cost-free, and its use as a cooling medium enables reduce condensation temperatures, thereby increasing system energy efficiency, so there are distinct advantages. The economy and feasibility of using the supercooling method to produce seawater ice slurry are analyzed by modeling a single-stage compression cycle ice production system with a double loop, which provides theoretical guidance for the practical application of technology for producing ice slurry from seawater. The results suggest that under the calculated conditions the average coefficient of performance of the system is 4.76 in winter and 3.61 in summer, which can maintain normal operation for a long time in a wide range of flow speeds. The system uses a plate-type water-cooled condenser, and the quantitative comparison with air source system suggests that the efficiency of water-cooled system is higher than that of air-cooled system at the concerning operating point, and the cool-storage capacity per unit investment is 2.85 kW∙h/yuan for water-cooled units at design conditions, corresponding to only 1.37 kW∙h/yuan for air-cooled units.

ice slurry; seawater; supercooling method; double-loop single-stage compression refrigeration cycle; plate heat exchanger; efficiency

2095-560X（2023）06-0491-08

TK02

A

10.3969/j.issn.2095-560X.2023.06.001

2023-01-23

2023-03-15

国家重点研发计划项目（2021YFE0112500）

宋文吉，E-mail：songwj@ms.giec.ac.cn

宋文吉, 黎福超, 陈明彪, 等. 过冷法海水冰浆制备系统能耗特性分析[J]. 新能源进展, 2023, 11(6): 491-498.

： SONG Wenji, LI Fuchao, CHEN Mingbiao, et al. Energy consumption characteristics analysis of seawater ice slurry generation system using supercooling method[J]. Advances in new and renewable energy, 2023, 11(6): 491-498.

宋文吉（1978-），男，博士，研究员，博士生导师，主要从事冰浆技术及创新应用、大规模储电及控制技术研究。