熱泵錯(cuò)流多效蒸發(fā)系統(tǒng)的模型與模擬

， ，，

(福州大學(xué) 石油化工學(xué)院，福建 福州 350108)

?開發(fā)與研究?

熱泵錯(cuò)流多效蒸發(fā)系統(tǒng)的模型與模擬

鄭曉桂，阮奇，周江建，張龍

(福州大學(xué) 石油化工學(xué)院，福建 福州 350108)

為了降低錯(cuò)流多效蒸發(fā)過程的能耗，將蒸汽噴射式熱泵技術(shù)運(yùn)用到錯(cuò)流多效蒸發(fā)系統(tǒng)中，并將系統(tǒng)中的物料衡算、熱量衡算方程組用矩陣方程的形式描述，使所建立的熱泵錯(cuò)流多效蒸發(fā)模型具有較強(qiáng)的通用性。用收斂速度快、穩(wěn)定性好的迭代法結(jié)合矩陣法模擬求解熱泵錯(cuò)流四效蒸發(fā)模型，結(jié)果表明。蒸汽噴射式熱泵技術(shù)具有較明顯的節(jié)能效果，熱泵的噴射系數(shù)越大，或熱泵的抽汽位置越往后效，節(jié)省的生蒸汽消耗量越大，節(jié)能效果越好，但同時(shí)蒸發(fā)器的傳熱面積也越大，存在最佳的熱泵噴射系數(shù)和最佳的熱泵抽汽位置。

熱泵 ； 錯(cuò)流多效蒸發(fā) ； 節(jié)能 ； 數(shù)學(xué)模型

Keywords:heat pump;cross flow multi-effect evaporation;energy saving;mathematical model

0 引言

蒸發(fā)是濃縮溶液(通常是水溶液)最重要、最常用的化工單元操作，但蒸發(fā)過程要除去大量的水分，需要消耗大量的加熱蒸汽(生蒸汽)，而采用多效蒸發(fā)能節(jié)省大量的生蒸汽，故多效蒸發(fā)在化工、輕工、食品工業(yè)、海水淡化及高鹽廢水處理等行業(yè)中得到了廣泛應(yīng)用。已報(bào)道的有關(guān)多效蒸發(fā)系統(tǒng)的模型、算法和節(jié)能措施的研究主要集在并流、平流和逆流[1-3]多效蒸發(fā)系統(tǒng)，而有關(guān)錯(cuò)流多效蒸發(fā)系統(tǒng)的研究報(bào)道較少。錯(cuò)流多效蒸發(fā)的特點(diǎn)是既融合了并流多效蒸發(fā)與逆流多效蒸發(fā)各自的優(yōu)點(diǎn)又克服了它們各自的缺點(diǎn)[4]，比較有代表性，故本文選擇錯(cuò)流多效蒸發(fā)系統(tǒng)為研究對(duì)象，并將節(jié)能效果較顯著的蒸汽噴射式熱泵技術(shù)運(yùn)用于該系統(tǒng)，建立其數(shù)學(xué)模型，研究模型的求解方法，模擬并分析熱泵技術(shù)對(duì)錯(cuò)流多效蒸發(fā)過程能耗的影響，這對(duì)降低該過程的能耗有重要的意義。

1 數(shù)學(xué)模型

1.1工藝流程

含蒸汽噴射式熱泵技術(shù)的錯(cuò)流多效蒸發(fā)系統(tǒng)的工藝流程如圖1所示。系統(tǒng)共有n效，其中逆流部分共q效，任意效用i表示，則i=1，2，…，q-1，q；并流部分共n-q效，任意效用q+j表j=1，2，…，n-q-1，n-q。用cq代表抽汽位置，cq可在如圖1所示的逆流部分的任意第i效，也可以在并流部分除末效(第n效)以外的任意第q+j效。來自鍋爐的高壓生蒸汽進(jìn)入蒸汽噴射熱泵，將第cq效低壓的部分二次蒸汽抽吸進(jìn)熱泵，兩者混合后作為加熱蒸汽從第1效進(jìn)入，第1效蒸發(fā)產(chǎn)生的二次蒸汽作為第2效的加熱蒸汽進(jìn)入第2效，以此類推直至末效，末效產(chǎn)生的二次蒸汽則全部進(jìn)入冷凝器(圖1中未畫出)冷凝除去。

原料液進(jìn)入第q+1效即并流部分第1效進(jìn)行蒸發(fā)，溶液依次在并流部分各效蒸發(fā)濃縮后，離開第n效即并流部分最后一效的溶液則進(jìn)入逆流部分第1效即第q效繼續(xù)蒸發(fā)，溶液依次在逆流部分各效蒸發(fā)濃縮后，離開第1效即逆流部分最后一效的濃縮液(完成液)作為產(chǎn)品。

圖1 熱泵錯(cuò)流多效蒸發(fā)系統(tǒng)流程圖

1.2系統(tǒng)物料衡算

假設(shè)溶質(zhì)不揮發(fā)，對(duì)逆流部分任意第i效和并流部分任意第q+j效蒸發(fā)器的溶質(zhì)進(jìn)行衡算有

F0x0=Fixi=Fq+jxq+j

(1)

逆流部分從第q效到第i效蒸發(fā)器料液中的溶質(zhì)組分的衡算式為：

(2)

并流部分從第q+1效到第q+j效蒸發(fā)器料液中的溶質(zhì)組分的衡算式為：

(3)

系統(tǒng)總的蒸發(fā)水量W0為逆流部分和并流部分各效蒸發(fā)水分量之和，即：

(4)

1.3系統(tǒng)熱量衡算

對(duì)各效蒸發(fā)器進(jìn)行熱量衡算可以求出各效蒸發(fā)水分量。假定各效加熱蒸汽的冷凝水在相應(yīng)的飽和溫度下排出，對(duì)圖1所示的并流部分任意第q+j效(j=1，2，…，n-q-1，n-q)蒸發(fā)器進(jìn)行熱量衡算，將二次蒸汽近似視為飽和蒸汽[5]，將溶液的比熱容近似認(rèn)為是由溶質(zhì)和水的比熱容成線性加和關(guān)系構(gòu)成，利用式(1)和式(3)并考慮熱損失及濃縮熱，采用與文獻(xiàn)[1，5]類似的方法推導(dǎo)、整理，可得：

Wq+j={αq+jDq+j+[F0c0-c*(Wq+1+…+Wq+j-1)]βq+j}ηq+j

(5)

對(duì)圖1所示的逆流部分任意第i效(i=1，2，…，q-1，q)蒸發(fā)器進(jìn)行熱量衡算，利用式(1)和式(2)并考慮熱損失及濃縮熱，采用與文獻(xiàn)[3]類似的方法推導(dǎo)、整理，可得：

Wi={αiDi+[F0c0-c*(Wi+1+…+Wq+

(6)

式(5)和式(6)中：ηi、ηq+j、αi、αq+j、βi、βq+j計(jì)算可參考文獻(xiàn)[]；Di為逆流各效蒸發(fā)器的加熱蒸汽量，定義為

(7)

式(7)中：當(dāng)i-1≠cq時(shí)，ui-1=0；當(dāng)i-1=cq時(shí)，ui-1=ucq；Dq+j為并流各效蒸發(fā)器的加熱蒸汽量，定義為：

Dq+j=Wq+j-1-Dsuq+j-11≤j≤n-q

(8)

式(8)中：當(dāng)q+j-1≠cq時(shí)，uq+j-1=0；

當(dāng)q+j-1=cq時(shí)，uq+j-1=ucq。

1.4傳熱面積的計(jì)算

逆流部分任意第i效蒸發(fā)器的傳熱面積按下式[7]計(jì)算：

Ai=Qi/KiΔti=Diri/KiΔti

(9)

式(9)中第i效的加熱蒸汽汽化潛熱ri和有效傳熱溫度差Δti根據(jù)文獻(xiàn)[1]的方法計(jì)算。若是并流部分任意第q+j效蒸發(fā)器，將式(9)中的變量下標(biāo)改為q+j即可。

1.5熱泵噴射系數(shù)的計(jì)算

在多效蒸發(fā)系統(tǒng)中采用熱泵技術(shù)的工作原理是，將高溫高壓的生蒸汽從蒸汽噴射式熱泵的噴嘴高速噴出，以造成噴嘴處的壓力遠(yuǎn)低于蒸發(fā)器的蒸發(fā)壓力，兩者形成一定的壓差，在此壓差的作用下，蒸發(fā)產(chǎn)生的部分低壓二次蒸汽被吸入噴射泵的混合室中與高壓生蒸汽混合，混合蒸汽在噴射泵的擴(kuò)散室中降速升壓升溫，而后進(jìn)入第1效蒸發(fā)器作為加熱蒸汽。蒸汽噴射式熱泵技術(shù)是將二次蒸汽壓縮后利用其潛熱從而達(dá)到節(jié)能(節(jié)省生蒸汽消耗量)的目的。噴射系數(shù)不僅是蒸汽噴射式熱泵的重要技術(shù)指標(biāo)和設(shè)計(jì)依據(jù)[8]，也是影響蒸發(fā)過程的重要參數(shù)?？刹捎梦墨I(xiàn)[9]導(dǎo)出的公式計(jì)算蒸汽噴射熱泵對(duì)第cq效(被抽汽效)二次蒸汽的噴射系數(shù)ucq，即

ucq=φ1φ2φ3ξ[1+(Hs-H0)/(H0-Hcq)]1/2-1

(10)

式(10)中Hs、H0和Hcq分別為生蒸汽、混合蒸汽和第cq效蒸發(fā)器二次蒸汽的焓值，J/kg；φ1、φ2和φ3分別為工作噴嘴、混合室和擴(kuò)散室的速度系數(shù)，可分別取值為0.95、0.975和0.9；ξ為修正系數(shù)，計(jì)算時(shí)其值可取1.1。

2 模型求解

多效蒸發(fā)系統(tǒng)計(jì)算就是聯(lián)立求解系統(tǒng)物料衡算式、熱量衡算式和傳熱速率方程式等，從而求出加熱蒸汽消耗量、各效蒸發(fā)水分量和各效傳熱面積等。本文建立的是含蒸汽噴射式熱泵技術(shù)的錯(cuò)流多效蒸發(fā)系統(tǒng)的數(shù)學(xué)模型，模型中的眾多計(jì)算式實(shí)質(zhì)上是一個(gè)復(fù)雜的非線性方程組，求解相當(dāng)困難。迭代法結(jié)合矩陣法[1，3]是求解復(fù)雜多效蒸發(fā)系統(tǒng)數(shù)學(xué)模型的高效算法，決定采用該法求解本文的數(shù)學(xué)模型，為此須將模型中的非線性方程組寫成矩陣方程的形式。

2.1系統(tǒng)物料衡算和熱量衡算方程組的矩陣形式

由式(5)、式(6)可得各效蒸發(fā)水分量的計(jì)算式n個(gè)，由式(4)可得總蒸發(fā)水分量計(jì)算式1個(gè)，上述計(jì)算式共n+1個(gè)，待求的未知量有n+1個(gè)(Ds，W1，W2，…，Wn)。將上述各式寫成矩陣方程，該矩陣是一個(gè)稀疏矩陣，將其寫成分塊矩陣會(huì)使原矩陣結(jié)構(gòu)簡(jiǎn)單，運(yùn)算更容易。分析該矩陣的結(jié)構(gòu)，將其分成8塊處理比較適宜，其具體結(jié)構(gòu)如下：

(11)

A1、A3矩陣如下：

其它塊矩陣參見文獻(xiàn)[6]。

2.2矩陣方程的特點(diǎn)及通用性

矩陣方程式(11)是描述含熱泵技術(shù)的進(jìn)料位置在第q+1效的錯(cuò)流(先并流后逆流)n效蒸發(fā)系統(tǒng)的數(shù)學(xué)模型，具有結(jié)構(gòu)清晰、意義明確、易于編程求解和通用性強(qiáng)等特點(diǎn)。當(dāng)熱泵在任意第cq(cq=1，2，…，n-1)效抽汽時(shí)，矩陣方程式(11)中的分塊矩陣A1、A3中的ucq≠0，而令u1，u2，…，ucq-1，ucq+1，…，un-1=0，即可實(shí)現(xiàn)熱泵在除第n效以外的任意效抽汽的功能；當(dāng)系統(tǒng)不設(shè)置熱泵時(shí)，令分塊矩陣A1和A3中的ucq=0，其它參數(shù)的取值同上所述，則該矩陣方程簡(jiǎn)化為沒有熱泵技術(shù)的錯(cuò)流多效蒸發(fā)模型。若將矩陣方程式(11)中的參數(shù)q設(shè)置為0，將矩陣方程中有關(guān)逆流部分的參數(shù)(i=1，2，…，q-1，q)均令為0，則剩下的j=1，2，…，n-1，n部分均為并流部分，進(jìn)料位置也由第q+1效變?yōu)榈?效，該矩陣方程簡(jiǎn)化為并流多效蒸發(fā)模型。同理，矩陣方程式(11)亦可簡(jiǎn)化為逆流多效蒸發(fā)模型，上述情況充分說明該模型具有很強(qiáng)的通用性。

2.3算例及討論

擬設(shè)計(jì)一個(gè)先并流后逆流的錯(cuò)流四效(并流和逆流各兩效)蒸發(fā)系統(tǒng)，用于將流量為6.67 kg/s的氫氧化鈉溶液從初始濃度0.15(質(zhì)量分?jǐn)?shù)，下同)濃縮到0.4。已知年操作時(shí)間為7 200 h，進(jìn)料溫度為25 ℃，生蒸汽溫度為180 ℃，末效冷凝器中的二次蒸汽溫度為45 ℃，各效蒸發(fā)器的傳熱系數(shù)分別為1 100、1 180、1 250、800 W/(m2·℃)，采用各效蒸發(fā)器傳熱面積相等的設(shè)計(jì)原則，按以下六種情況進(jìn)行設(shè)計(jì)：a.無節(jié)能措施的錯(cuò)流四效蒸發(fā)，編號(hào)①；b.在①的基礎(chǔ)上加蒸汽噴射式熱泵節(jié)能措施，且抽氣位置在第1效，噴射系數(shù)為0.1、0.15、0.20，編號(hào)分別為②、③、④；c.在③噴射系數(shù)為0.15的基礎(chǔ)上，抽氣位置在第3效，編號(hào)⑤；d.在④噴射系數(shù)為0.20的基礎(chǔ)上，抽氣位置在第3效，編號(hào)⑥。

先并流后逆流的錯(cuò)流多效蒸發(fā)系統(tǒng)的通用數(shù)學(xué)模型及其求解算法已用Visual Basic 6.0 語言開發(fā)成通用設(shè)計(jì)軟件，分別將上述六種設(shè)計(jì)情況的設(shè)計(jì)條件輸入并運(yùn)行軟件，設(shè)計(jì)結(jié)果如表1 所示。表1中Ai為任意第i效的傳熱面積，因按各效等面積原則設(shè)計(jì)，則四效總的傳熱面積是其四倍；表1中的生蒸汽消耗量對(duì)無熱泵的第①這種情況即為第1效的加熱蒸汽量D1，對(duì)有熱泵的第②～⑥種情況則為進(jìn)入蒸汽噴射熱泵的生蒸汽量Ds。

表1 算例模擬計(jì)算結(jié)果

對(duì)于蒸發(fā)器材質(zhì)為不銹鋼或碳鋼的多效蒸發(fā)系統(tǒng)，生蒸汽費(fèi)用占年總費(fèi)用的88%～96%[9](在后面的分析中為了論述的方便取此值為90%)，而蒸發(fā)系統(tǒng)的設(shè)備投資折舊費(fèi)用(主要取決于蒸發(fā)器的傳熱面積)僅占年總費(fèi)用的4%～12%(在后面的分析中為了論述的方便取此值為10%)。根據(jù)文獻(xiàn)[10]的觀點(diǎn)并結(jié)合表1的設(shè)計(jì)結(jié)果作如下的分析與討論：

2.3.1 第②種情況采用熱泵技術(shù)噴射系數(shù)為0.1且抽氣位置在第1效的設(shè)計(jì)與第①種無節(jié)能措施的設(shè)計(jì)相比，生蒸汽的用量可以節(jié)省11.22%，而傳熱面積僅增加44.67%，占年總費(fèi)用90%左右的生蒸汽節(jié)省11.22%帶來的經(jīng)濟(jì)效益大于僅占年總費(fèi)用10%左右的設(shè)備投資折舊費(fèi)用因傳熱面積增加44.67%而增加的費(fèi)用，這說明熱泵技術(shù)對(duì)于錯(cuò)流多效蒸發(fā)系統(tǒng)是經(jīng)濟(jì)合理的有效節(jié)能措施。

2.3.2 ②、③和④三種情況的設(shè)計(jì)結(jié)果表明，熱泵的抽汽位置相同時(shí)，熱泵的噴射系數(shù)越大，節(jié)省的生蒸汽量越多，節(jié)能效果越好，但蒸發(fā)器的傳熱面積也增加越多。

2.3.3 ⑤和③兩種情況、⑥和④兩種情況的設(shè)計(jì)結(jié)果均表明，熱泵的噴射系數(shù)相同時(shí)，熱泵的抽汽位置越往后效，節(jié)省的生蒸汽量越多，節(jié)能效果越好，但蒸發(fā)器的傳熱面積也增加越多。

2.3.4 2.3.2和2.3.3中討論的現(xiàn)象可做如下分析：一方面，熱泵的噴射系數(shù)ucq越大，被抽吸利用的二次蒸汽的量就越多，或熱泵的抽汽位置越往后效，被抽的二次蒸汽的汽化潛熱已在前面效得到越多次的利用，這些均使二次蒸汽潛熱的利用更為充分，因而生蒸汽的消耗量會(huì)越節(jié)省，達(dá)到越節(jié)能的目的；另一方面，由式(11)可知，在生蒸汽的溫度和焓一定時(shí)，ucq越大，或抽汽位置越往后效即被抽的二次蒸汽的溫度越低焓越小，均使混合蒸汽的焓越小，混合蒸汽的溫度也越低，致使蒸發(fā)系統(tǒng)的總有效傳熱溫度差越小，各效蒸發(fā)器的有效傳熱溫度差也越小，各效傳熱面積則增大；反之亦然。

綜上所述，存在最佳的噴射系數(shù)ucq和最佳的抽汽位置cq，使生蒸汽的操作費(fèi)用和傳熱面積的投資費(fèi)用之和最小，這是熱泵多效蒸發(fā)系統(tǒng)的優(yōu)化設(shè)計(jì)問題，有待今后繼續(xù)深入研究。綜合考慮，生蒸汽節(jié)省量和傳熱面積增加量的影響，表1中第④種設(shè)計(jì)結(jié)果最好，第⑥種設(shè)計(jì)結(jié)果最差。

3 結(jié)論

所建立的含蒸汽噴射式熱泵技術(shù)的錯(cuò)流多效蒸發(fā)系統(tǒng)的數(shù)學(xué)模型，用由八塊分塊矩陣構(gòu)成的矩陣方程表達(dá)，可實(shí)現(xiàn)熱泵在任意效抽汽，易于簡(jiǎn)化成無熱泵的錯(cuò)流多效蒸發(fā)模型，亦可簡(jiǎn)化成并流或逆流多效蒸發(fā)模型，具有很強(qiáng)的通用性。蒸汽噴射式熱泵技術(shù)具有較明顯的節(jié)能效果。熱泵的噴射系數(shù)越大，或熱泵的抽汽位置越往后效，節(jié)省的生蒸汽消耗量越大，節(jié)能效果越好，但同時(shí)蒸發(fā)器的傳熱面積也越大，存在最佳的熱泵噴射系數(shù)和最佳的熱泵抽汽位置，這有待今后繼續(xù)深入研究。

符號(hào)說明：

A，任意效蒸發(fā)器傳熱面積，m2；c0，原料液的比熱容，J/(kg·℃) ；c*，水的比熱容，J/(kg·℃) ；D，任意效加熱蒸汽流量，kg/s；Ds，生蒸汽的流量，kg/s；F0，原料液的流量，kg/s；F，任意效完成液流量，kg/s；H，任意效二次飽和蒸汽的焓，J/kg；K，任意效蒸發(fā)器的傳熱系數(shù)，W/(m2·℃-1) ；Q，任意效傳熱速率，W；T0，混合蒸汽溫度，℃；Ts，生蒸汽溫度，℃；T，任意效二次飽和蒸汽溫度，℃；t，任意效溶液沸點(diǎn)，℃；W0，總蒸發(fā)水分量，kg/s；W，任意效蒸發(fā)水分量，kg/s；x0，原料液濃度，質(zhì)量分?jǐn)?shù)；x，任意效完成液濃度，質(zhì)量分?jǐn)?shù)；α，蒸發(fā)系數(shù)；β，自蒸發(fā)系數(shù)；η，熱利用系數(shù)。

下角標(biāo)：i，逆流部分的任意第i效蒸發(fā)器；q+j，并流部分的任意第q+j效蒸發(fā)器。

[1]阮 奇，陳文波，黃詩煌，等. 復(fù)雜并流多效蒸發(fā)系統(tǒng)的數(shù)學(xué)模型與矩陣解法[J].中國(guó)工程科學(xué)，2001，3(4)：36-41.

[2]李 玲，阮 奇. 平流多效蒸發(fā)系統(tǒng)的數(shù)學(xué)模型與求解[J].化工學(xué)報(bào)， 2009，60(1)：104-111.

[3]阮 奇，黃詩煌，葉長(zhǎng)燊，等. 復(fù)雜逆流多效蒸發(fā)系統(tǒng)常規(guī)設(shè)計(jì)的模型與算法[J].化工學(xué)報(bào)，2001，52(7)：616-621.

[4]時(shí) 鈞，汪家鼎，余國(guó)琮，等. 化學(xué)工程手冊(cè)[M].北京：化學(xué)工業(yè)出版社，1996.

[5]上?；W(xué)院，成都科技大學(xué)，大連工學(xué)院．化學(xué)工程(第一冊(cè))[M].北京：化學(xué)工業(yè)出版社，1980：316-352.

[6]王 勇. 復(fù)雜多效蒸發(fā)過程模擬與操作優(yōu)化研究[D].福州：福州大學(xué)化學(xué)化工學(xué)院，2006.

[7]譚天恩，竇 梅，等. 化工原理(上冊(cè)).4版[M].北京：化學(xué)工業(yè)出版社，2013.

[8]索科洛E R，津格爾H M. 噴射器[M].黃秋云譯. 北京：科學(xué)出版社，1977.

[9]王 權(quán)，向雄彪. 蒸汽噴射壓縮器噴射系數(shù)計(jì)算方法研究[J].太陽能學(xué)報(bào)，1997，18(3)：314-321.

[10]梁 虎，王 黎，朱 平. 多效蒸發(fā)系統(tǒng)優(yōu)化設(shè)計(jì)研究[J].化學(xué)工程，1997，25(6)：48-51.

ModelingandSimulationofCrossFlowMulti-effectEvaporationSystemwithHeatPump

ZHENGXiao-gui,RUANQi,ZHOUJiang-jian,ZHANGLong

(School of Chemical Engineering,Fuzhou University,Fuzhou 350108,China)

In order to reduce the energy consumption of crossflow multi-effect evaporation process,steam jet heat pump technology is applied to the cross flow multi-effect evaporation system. Equations of material balance and energy balance of the system are expressed in the form of matrix equation so that the model of cross flow multi-effect evaporation with heat pump has strong versatility. Iteration method combined with matrix method is used to simulate and solve the quadruple effect cross-flow evaporation system.This algorithm has some advantages of good stability and quick convergence speed.The result indicates that steam jet heat pump technology has obvious energy saving effect.The larger injection coefficient or more backward extraction location of heat pump is set,the more fresh steam consumption is saved and the better energy saving effect is obtained.At the same time,the heat transfer area of evaporator is larger.So,there exist optimal injection coefficient and extraction location of heat pump.

2014-02-25

國(guó)家大學(xué)創(chuàng)新創(chuàng)業(yè)訓(xùn)練計(jì)劃項(xiàng)目(201210386012)；福建省大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計(jì)劃項(xiàng)目(FJ121038605)

鄭曉桂(1991-)，女，本科，研究方向?yàn)榛み^程模擬與優(yōu)化、傳質(zhì)與分離；聯(lián)系人：阮 奇，教授，E-mail:hys@fzu.edu.cn。

TQ050.2

A

1003-3467(2014)05-0023-05