工业水处理, 2019, 39(7): 19-23 doi: 10.11894/iwt.2019-0037

试验研究

3倍浓缩低温多效海水淡化蒸发设备传热系数研究

刘治川,

Study on heat transfer coefficient on triple-concentrated low temperature multi-effect seawater desalination evaporator

Liu Zhichuan,

收稿日期: 2019-04-19  

基金资助: 国家科技支撑计划.  2015BAB08B00
北京市重点研发计划.  Z171100000717009

Received: 2019-04-19  

Fund supported: 国家科技支撑计划.  2015BAB08B00
北京市重点研发计划.  Z171100000717009

作者简介 About authors

刘治川(1971-),高工电话:13691070468,E-mail:liuzhichuan@ceedi.cn , E-mail:liuzhichuan@ceedi.cn

摘要

直接计算3倍浓缩低温多效海水淡化蒸发设备传热系数所需的一些参数难以测量。推导出平均传热系数计算公式。借助流程模拟软件PROⅡ建立1 000 t/d装置数学模型,通过实际运行数据验证模型可靠性,再利用模型求取部分参数。在此基础上计算了2.5~3倍浓缩范围内的平均传热系数。研究表明,产品水闪蒸汽量约占总凝液量的2%,总沸点升约占总换热温差38%,3倍浓缩平均传热系数可达2 900 W/(m2·℃)。平均传热系数在3倍浓缩时数值最小,但随浓缩比的变化不大。

关键词: 3倍浓缩 ; 低温多效 ; 海水淡化 ; 传热系数

Abstract

In the triple-concentrated low temperature multi-effect seawater desalination evaporator, it is difficult to measure some parameters which are necessary for transfer coefficient of evaporator. The formula for calculating average heat transfer coefficient was derived. The mathematical model of the 1 000 t/d plant was established by using PROⅡ process simulation software. The reliability of the model was verified by the actual running data. Then the partial parameters in the formula could be obtained by using the model. On the basis, the average heat transfer coefficient was calculated in the range of 2.5-3 times concentration. The results showed that the total flash vapor of product water accounted for about 2% of the total condensate, and the total boiling point rise accounted for about 38% of the total heat transfer temperature difference. The average heat transfer coefficient could reach 2 900 W/(m2·℃) at 3 times concentration. The average heat transfer coefficient was the smallest when the concentration ratio was 3, but it did not change much with different concentration ratio.

Keywords: triple-concentrated ; low temperature multi-effect ; seawater desalination ; heat transfer coefficient

PDF (0KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

刘治川. 3倍浓缩低温多效海水淡化蒸发设备传热系数研究. 工业水处理[J], 2019, 39(7): 19-23 doi:10.11894/iwt.2019-0037

Liu Zhichuan. Study on heat transfer coefficient on triple-concentrated low temperature multi-effect seawater desalination evaporator. Industrial Water Treatment[J], 2019, 39(7): 19-23 doi:10.11894/iwt.2019-0037

3倍浓缩低温多效海水淡化技术处理后的浓盐水总含盐质量分数达10%,这在排放盐水可资源化利用的条件下是一种优选方式,有利于下游盐场和盐化工企业节省土地及能源,促进循环经济发展。为减少结硫酸钙垢的可能性,3倍浓缩低温多效海水淡化技术一般会采用顺流流程,其特点是蒸发设备高温端低浓度,低温端高浓度,各效温度浓度分布不同于目前2倍浓缩及以下技术常用的平流和逆流设备,后两者蒸发设备处理后的盐水总含盐质量分数小于7%,特点是高温端高浓度,低温端低浓度。因此前者设备的传热系数可能与后两者不同,需要测定前者的传热系数。

低温多效蒸发设备通常采用水平管降膜蒸发器。目前该类设备传热系数缺乏通用性高的理论公式,多来自于经验公式、试验测试或已运行装置的经验数据。经验公式基于有限的试验数据,如果超出了参数变化范围,或操作条件有所差异,公式将会给出与实际偏差较大的预测值1。试验测试,通常都是采用小型单效设备针对独立的工艺条件进行。工业设备效数多,每效有大量传热管,管内外传热过程互相耦合,各效间热量流动互相影响,实际运行工况与小型单一测试环境有较大差距,因此实际运行状态下的综合传热性能可能不同于小型测试结果,应当了解其实际传热性能。

目前对于水平管降膜蒸发器传热研究的文献报道多是利用小型设备在特定条件下关注相关因素对传热系数的影响2-4,较少提及工业环境下的传热系数及其测定方法。出于投资和操作的考虑,工业设备通常不会设置直接测定传热系数所需的全部仪表。笔者根据工业装置仪表设置状况,推导出测量参数和计算参数相结合的平均传热系数测定公式,借助流程模拟软件PROⅡ建立1 000 t/d 3倍浓缩顺流工业中试装置的数学模型,利用装置实际运行数据验证模型的可靠性后求取计算参数,进而计算平均传热系数,以期为大型3倍浓缩顺流装置提供传热数据参考。

1 蒸发设备概况

1.1 蒸发流程

3倍浓缩顺流蒸发工艺流程见图1

图1

图1   3倍浓缩顺流蒸发工艺流程


蒸发器共14效,分成9个效组。海水在效组内是平行进料,在效组间的流动方向相对于蒸汽是顺流。海水经冷凝器及3台预热器预热后进入首效。为充分利用高品质进料蒸汽压力能,进料蒸汽通过喷射器将第9效部分二次蒸汽抽吸至首效一起作为加热蒸汽。为减少进料蒸汽用量,利用第10、5、2效二次蒸汽逐级预热海水,各自的凝液流入本效闪蒸罐。第1效管侧蒸汽凝液包括进料蒸汽和来自第9效的二次蒸汽抽气,经冷凝水缓冲罐和冷凝水泵送出界区。第1~13效壳侧海水产生的二次蒸汽在第2~14管侧被海水冷凝,成为产品水。第14效壳侧二次蒸汽在冷凝器管侧被壳侧海水冷凝,成为产品水。为回收产品水中的热量,第2~13效产品水设有闪蒸罐,通过接收本效及来自上一级闪蒸罐(第2效没有)的产品水,再向下一效闪蒸,回收闪蒸汽。第14效和冷凝器管侧产品水与第13效闪蒸罐的产品水混合,经产品水泵排出。每效海水泵出口设有循环线,返回至本效进料入口,以满足各效换热管喷淋量要求。

在进出界区处设有海水、进料蒸汽、冷凝水、盐水、产品水流量计,首效设有压力表,冷凝器入口、预热器出口和各效蒸发器盐水侧设有温度计。

1.2 蒸发设备结构

蒸发设备总长80 m。每效蒸发器内径2.6 m,3 524根换热管。前6效换热管长3.5 m,后8效换热管长2.47 m,14效总换热面积8 118.4 m2。蒸发器布管情况见图2,铝合金换热管结构见图3

图2

图2   蒸发器布管


图3

图3   换热管结构


2 测定方法

在推导公式前先对涉及到的参数符号进行统一说明:

A——总换热面积,m2

Ai——i效换热面积,m2

B——盐水出料量,t/h;

Bi——i效盐水出料量,t/h;

Cp——进料海水定压热容,J/(g·℃);

F——海水进料量,t/h;

Fi——i效海水进料量,t/h;

M——动力蒸汽量,t/h;

Mc——去冷凝器蒸汽量,t/h;

Mi——i效凝液量,t/h;

Mai——i效去热泵蒸汽量,t/h;

Mpi——i效去海水预热器蒸汽量,t/h;

Mfi——产品水闪蒸罐进i效蒸汽量,t/h;

Ki——i效传热系数,W/(m2·℃);

${\bar K}$——平均传热系数,W/(m2·℃);

Qi——i效传热量,kJ/h;

QL——热损失,kJ/h;

Ri——i效蒸发潜热,J/g;

${\bar R}$——平均蒸发潜热,J/g;

S——盐度,%;

T1——1效加热蒸汽饱和温度,℃;

TS——饱和温度,℃;

t1——1效海水进料温度,℃;

t14——14效盐水温度,℃;

tc——冷凝器海水进料温度,℃;

Δi——i效海水沸点升,℃;

Δt——总有效换热温差,℃;

Δti——i效有效换热温差,℃。

对顺流多效蒸发流程进行物料衡算和热量衡算,流程示意见图4,衡算结果见式(1)、式(2)。

图4

图4   顺流多效蒸发物料衡算和热量衡算


(1)

(2)

由于各效凝液处于重力流动状态,用于直接计算各效传热量的Mi难以测量。

图1可知下列初始条件:

Mpi=(i=1~14,但i≠2、5、10)

根据式(1)有:

根据式(2)有:

Qi=KiAiΔti,则

Ki=${\bar K}$,tit/14,则得出式(3):

(3)

其中,则得出式(4):

(4)

由式(4)可知,MFBt1tcT1t14设有相关仪表可测,A是已知定值,Cp和${\bar R}$是物性数据,只要获得ΔiMfi,即可求取值。

根据文献〔5〕,44~72 ℃范围内,蒸发潜热与温度拟合关系式为R=-2.434 5TS+2 503.6,单位J/g。由于各效温度近似等温差分布,为简化计算,取首末效温度下的平均值。进料海水含盐质量分数为3.3%,温度在18~68 ℃间变化,Cp计算公式可简化为温度的一次关系式6,即:Cp=1.004 9-0.016 2S+3.526×10-4S2-(3.250 6-1.479 5S+0.077 7S2)×10-4t,单位J/(g·℃)。Cptct1下计算的平均值。

Δi与各效海水的浓度和温度有关。各效盐水侧设有温度计,产品水侧不设温度计,各效盐水侧也不设在线浓度仪,因此Δi无法直接测出或通过沸点升与海水浓度温度关系直接推算出。第2~13效产品水闪蒸罐的闪蒸汽量也难以直接测量。这两个参数需要通过数学模型计算获得。

3 装置建模及模型验证

流程模拟软件PROⅡ可用于描述3倍浓缩低温多效海水淡化过程,可较准确预测海水浓度在10%范围内不同浓度温度下海水的沸点升7。利用PROⅡ建立图1流程的数学模型,其热力学模型选择电解质模型。模型关键输入值与实际运行值相同,实际运行值见表1。即使最终浓缩比相同,各效蒸发器温度浓度分布不同,总沸点升也可能不同。由于模型中是通过规定压力来确定各效温度,为减少试算工作量,又保证各效的温度浓度分布与实际结果接近,仅选取模型中关键效(第1、2、5、9、10、14效)盐水温度输入值与运行值相同。其余效的盐水温度对物料和热量衡算影响小,因而对沸点升的影响也小。蒸发设备在海边闲置了较长时间并经改造后重新启用,保温层一定程度受损,且受环境温度和风速等影响,散热量难以准确确定。为模型输入热损失估算值,模型计算结果与运行结果的比较见表2

表1   装置实际运行值

工况浓缩比F/(t·h-1B/(t·h-1M/(t·h-1t1/℃T1/℃t14/℃
13.0260.0319.853.42763.072.3146.86
2.7760.4021.833.37558.070.2446.04
2.5060.8724.393.2356.067.8145.03
23.0558.2719.133.37859.070.7444.77
2.7758.1221.013.24956.568.3844.19
2.5758.4522.722.98557.568.5845.24

新窗口打开| 下载CSV


表2   模型计算值和实际运行值比较

工况浓缩比QL/%M1/(t·h-1偏差/%tc/℃偏差/%
计算值运行值计算值运行值
13.021.434.494.49022.1222.23-0.495
2.772.094.8954.890.10222.6020.838.49
2.502.074.6394.64-0.02222.1220.497.94
23.051.94.8084.81-0.04219.6418.217.84
2.772.064.6364.64-0.08620.5419.336.25
2.571.514.2934.290.07021.3822.35-4.33

新窗口打开| 下载CSV


为便于计算结果具有可比性,两工况测量数据集中取自调试运行早期前后相接的一段时间。其间,换热管表面目测洁净,没有明显污垢。为保证数据取值具有代表性和准确性,各组数据为设备稳定运行后8~25 h平均值。

蒸发设备效数多、流程长、众多影响因素互相关联,运行过程具有时变和滞后等特点,而且不能完全消除仪表的异常波动,因此数学模型计算值与实际运行数据会存在一定偏差。表2两个工况下M1计算值与测量值偏差很小时,tc的平均偏差分别为5.64%和6.14%,说明所建模型与实际值吻合较好,可以用沸点升和淡水闪蒸气量的计算值代替实际值。通过累加上一效壳侧海水和本效管侧蒸汽凝液的温度差求取总沸点升,累加产品水闪蒸气量获得总闪蒸气量。

4 结果与讨论

4.1 产品水总闪蒸汽量、总沸点升

X代表总凝液量计算中可直接测量的部分,Y是总闪蒸气量,是总凝液量计算中需要借助模型计算的部分。产品水总闪蒸气量和总沸点升等计算结果见表3

表3   淡水总闪蒸气量和总沸点升计算结果

工况浓缩比X/(t·h-1Y/(t·h-1T1-t14/℃Δt/℃
13.0239.440.8440.289.8025.4515.6572.090.385
2.7738.130.7838.919.4224.2014.7802.000.389
2.5036.040.7036.748.6922.7814.0911.910.381
23.0538.480.8639.349.9325.9716.0402.190.382
2.7736.700.7337.439.2624.1914.9301.950.383
2.5735.230.6835.918.8423.3414.5041.890.379

新窗口打开| 下载CSV


表3可知,产品水总闪蒸汽量占总凝液的比例小,平均值为2.0%;总沸点升占总温差的比例大,2.50~3.05倍浓缩比时,在37.9%~38.9%。

4.2 平均传热系数

系统运行条件复杂,不同工况、不同浓缩比下各效海水喷淋量、蒸汽量、温度条件都不相同。为便于计算结果可量化比较,对影响传热系数的主要参数用平均值描述:规定海水喷淋量、蒸汽量和蒸汽侧气速为14效的平均值,海水侧温度取6个关键效盐水侧温度的平均值。由于每效海水采用循环进料方式,同一工况中实测的各效海水进料量表现相对稳定。温度数据来自实测值,蒸汽量和蒸汽侧气速来自各组模型计算值。各参数平均值计算结果见表4

表4   影响传热系数相关参数平均值

工况浓缩比海水喷淋量/(L·m-1·s-1海水侧温度/℃蒸汽量/(kg·h-1蒸汽侧气速/(m·s-1
13.020.049 959.432 8908.72
2.770.050 759.542 8338.87
2.500.050 257.692 5578.66
23.050.046 658.292 8649.19
2.770.046 957.742 7229.34
2.570.046 857.922 5978.83

新窗口打开| 下载CSV


两个工况传热系数的计算结果见图5。实际运行过程不受控制的因素多,即使是海水进料相同的工况,也很难进行不同浓缩比下理想化的对比试验。因此,对于图5中传热系数随浓缩比先增加再减少的现象,试图通过有限的测量数据从单一的微观理论解释宏观综合表现是困难的。但总体而言,在每个工况的操作范围内,浓缩比约为3时传热系数最低,浓缩比对传热系数的影响较小。蒸汽量大、喷淋量大且蒸发温度高的工况1的传热系数整体高于工况2,其在3倍浓缩时的传热系数为2 900 W/(m2·℃)。

图5

图5   传热系数与浓缩比关系


数学模型可直接计算各效凝液量和换热量,并可根据式(3)计算传热系数。模型相关直接计算值与前述推导计算值的比较见表5。从表5可知,利用与实际运行数据吻合较好的数学模型直接计算平均传热系数同样便捷可靠。

表5   总凝液量及平均传热系数的模型值与推导值比较

工况浓缩比偏差/%K/(W·m-2·℃-1偏差/%
模型值推导值模型值推导值
13.0239.9540.28-0.822 8752 902-0.92
2.7739.1338.910.572 9902 9780.39
2.5036.9636.740.602 9682 9550.43
23.0539.5439.340.512 7852 7760.34
2.7737.6037.430.452 8492 8420.25
2.5735.8835.91-0.082 7982 805-0.25

新窗口打开| 下载CSV


5 结论

3倍浓缩顺流工业装置受到仪表设置条件限制,蒸发设备各效凝液量和沸点升无法测量时,利用经过实际运行数据验证的数学模型,可便捷可靠地完成平均传热系数计算。蒸发设备首末效操作温度47~72 ℃,各效蒸发器顶部海水平均喷淋量0.05 L/(m·s),管内平均气速8.7 m/s条件下,3倍浓缩顺流蒸发设备平均传热系数可达到2 900 W/(m2·℃)左右。2.5~3倍浓缩比范围内,淡水总闪蒸汽量占总凝液的比例小,为2%左右;总沸点升占总温差的比例大,为38%左右。研究中还发现,平均传热系数在浓缩比3时数值最小,但随浓缩比的变化不大。根据1 000 t/d 3倍浓缩顺流工业化中试装置计算的3倍浓缩平均传热系数可为同类大型装置提供传热数据参考。

参考文献

龚路远.大型水平管降膜蒸发器热力性能和参数分布特性研究[D].大连:大连理工大学, 2017.

URL     [本文引用: 1]

Shen Shengqiang , Chen Xue , Mu Xingsen , et al.

Heat transfer characteristics of horizontal tube falling film evaporation for desalination

[J]. Desalination and Water Treatment, 2015, 55 (12): 3343- 3349.

DOI:10.1080/19443994.2014.946726      [本文引用: 1]

解利昕, 周文萌, 陈飞.

水平管降膜蒸发器的传热性能

[J]. 化工进展, 2014, 33 (11): 2878- 2881.

URL    

齐春华, 李炎, 俞永江, .

水平管降膜蒸发海水淡化的传热试验研究

[J]. 化学工业与工程, 2011, 28 (5): 1- 5.

URL     [本文引用: 1]

严家騄, 俞晓福, 王永青. 水和水蒸气热力性质图表[M]. 北京: 高等教育出版社, 2004: 2.

[本文引用: 1]

Sharqawy M H , Lienhard J H , Zubair S M .

Thermophysical properties of seawater:a review of existing correlations and data

[J]. Desalination and Water Treatment, 2010, 16 (1/2/3): 354- 380.

URL     [本文引用: 1]

刘治川.

三倍浓缩顺流低温多效海水淡化工艺的模拟与优化

[J]. 中国给水排水, 2018, 34 (11): 98- 103.

URL     [本文引用: 1]

/