中藥濃縮用機(jī)械蒸汽再壓縮系統(tǒng)模型與仿真*

湯添鈞 劉成剛 肖聰 李偉華

（蘇州科技大學(xué)環(huán)境科學(xué)與工程學(xué)院）（河北樂恒化工設(shè)備制造有限公司)

應(yīng)用研究

中藥濃縮用機(jī)械蒸汽再壓縮系統(tǒng)模型與仿真*

湯添鈞**劉成剛 肖聰 李偉華

（蘇州科技大學(xué)環(huán)境科學(xué)與工程學(xué)院）（河北樂恒化工設(shè)備制造有限公司)

針對(duì)中藥濃縮用機(jī)械蒸汽再壓縮（MVR）系統(tǒng)，利用能量守恒式、壓縮比與換熱溫差關(guān)系式、蒸汽壓縮機(jī)特性曲線擬合式建立了系統(tǒng)的數(shù)學(xué)模型。以某中藥制藥企業(yè)的MVR系統(tǒng)為例，將數(shù)學(xué)模型求解值與實(shí)際測(cè)量值進(jìn)行對(duì)比，發(fā)現(xiàn)模型求解精度較高，實(shí)際傳熱系數(shù)與壓縮機(jī)總效率是造成模型誤差的主要因素。對(duì)模型進(jìn)行仿真，得到蒸發(fā)器傳熱系數(shù)下降后蒸發(fā)量與壓縮比的預(yù)測(cè)值，對(duì)指導(dǎo)生產(chǎn)具有一定的意義。

中藥濃縮 機(jī)械蒸汽再壓縮 建模 仿真 蒸發(fā)器 傳熱系數(shù) 壓縮比

隨著經(jīng)濟(jì)的快速增長，我國正面臨著嚴(yán)重的能源短缺問題。2013年醫(yī)藥制造業(yè)共消耗2 179.11萬t標(biāo)準(zhǔn)煤[1]。作為中藥制備中一項(xiàng)重要的耗能單元，蒸發(fā)濃縮工藝存在著蒸汽用量大、利用效率低等問題，具有較大的節(jié)能潛力。為提高蒸發(fā)濃縮工藝蒸汽的利用效率，多效蒸發(fā)（MED）技術(shù)已得到較為廣泛的應(yīng)用[2]。但與MED技術(shù)相比，一種設(shè)計(jì)合理的機(jī)械蒸汽再壓縮（MVR）系統(tǒng)，具有蒸汽利用效率高、結(jié)構(gòu)更為緊湊的特點(diǎn)[3]，目前已經(jīng)進(jìn)入國家重點(diǎn)節(jié)能技術(shù)推廣目錄。

在MVR系統(tǒng)中，利用蒸汽壓縮機(jī)將蒸發(fā)出的二次蒸汽進(jìn)行再壓縮，使之重新成為熱源加熱物料，從而實(shí)現(xiàn)了二次蒸汽的全部回收[4]。在試驗(yàn)研究方面，文獻(xiàn) [5-6]測(cè)試了離心式與單螺桿式壓縮機(jī)驅(qū)動(dòng)的MVR試驗(yàn)臺(tái)的蒸發(fā)能力；在實(shí)際工程應(yīng)用方面，文獻(xiàn) [7-8]探究了MVR技術(shù)在工業(yè)含鹽廢水處理方面的節(jié)能效果。目前，MVR系統(tǒng)在中藥濃縮領(lǐng)域的理論與應(yīng)用研究，正逐步受到學(xué)者與制藥企業(yè)的重視。

本文首先對(duì)中藥濃縮用MVR系統(tǒng)建立數(shù)學(xué)模型，并實(shí)測(cè)了某中藥制藥企業(yè)提取車間MVR系統(tǒng)，然后將模型求解值與實(shí)際測(cè)量值進(jìn)行對(duì)比，最后對(duì)數(shù)學(xué)模型進(jìn)行仿真，得到傳熱系數(shù)下降后蒸發(fā)量與壓縮比的預(yù)測(cè)值。

1 基于MVR的中藥濃縮系統(tǒng)及數(shù)學(xué)模型

1.1 系統(tǒng)描述

適應(yīng)于中藥蒸發(fā)濃縮的特點(diǎn)，中藥濃縮用MVR系統(tǒng)基本流程如圖1所示，系統(tǒng)主要包含蒸汽壓縮機(jī)、蒸發(fā)器、氣液分離器、真空泵及其它輔助設(shè)備。

圖1 MVR系統(tǒng)流程

蒸發(fā)流程為：物料罐中的藥物原液由溫度t0經(jīng)兩級(jí)預(yù)熱至飽和溫度te后進(jìn)入蒸發(fā)器，在蒸發(fā)器管程內(nèi)沸騰蒸發(fā)；氣液兩相流通過蒸發(fā)器下部的平衡管進(jìn)入分離器中，氣相被壓縮機(jī)吸入，液相則通過分離器底部的管道返回蒸發(fā)器中。壓縮后的過熱蒸汽經(jīng)噴水處理，變成同壓力下的飽和蒸汽進(jìn)入蒸發(fā)器殼程，放出潛熱，最終進(jìn)入冷凝水罐。由于冷凝水有較高的溫度，使其進(jìn)入預(yù)熱器與原液換熱，將熱量傳遞給原液，從而完成蒸發(fā)流程。

系統(tǒng)中蒸汽熱力過程見圖2。當(dāng)二次蒸汽從點(diǎn)A經(jīng)過一多變過程壓縮至點(diǎn)S時(shí)，其絕對(duì)壓力上升，對(duì)應(yīng)的飽和溫度也相應(yīng)上升。經(jīng)過噴水處理消除過熱，至飽和狀態(tài)點(diǎn)B后進(jìn)入蒸發(fā)器，熱交換放出潛熱后到達(dá)點(diǎn)C，回收的潛熱量為圖中陰影部分。蒸汽在壓縮機(jī)進(jìn)出口壓力下對(duì)應(yīng)的飽和溫差Δt，即為MVR系統(tǒng)的換熱溫差。

圖2 蒸汽熱力過程

1.2 系統(tǒng)數(shù)學(xué)模型

建立MVR系統(tǒng)數(shù)學(xué)模型遵從以下幾點(diǎn)假設(shè)：

（1）忽略系統(tǒng)內(nèi)的能量損失；

（2）二次蒸汽在蒸發(fā)器內(nèi)完全冷凝；

（3）忽略物料濃度升高對(duì)沸點(diǎn)的影響。

1.2.1 系統(tǒng)能量守恒

在MVR系統(tǒng)中，蒸發(fā)器的常見形式為管殼式，將其簡化成一個(gè)基本數(shù)學(xué)模型，其傳熱過程的計(jì)算式如下：

式中Φ——換熱器換熱量，kW;

A——換熱面積，m2；

K——傳熱系數(shù)，kW/(m2·K)；

tc——蒸發(fā)器殼程飽和蒸汽溫度，℃；

te——原液蒸發(fā)溫度，℃；

Δt——換熱溫差，℃。

換熱量的大小對(duì)壓縮機(jī)的工作狀態(tài)至關(guān)重要。壓縮機(jī)處理的二次蒸汽量有限，蒸發(fā)出質(zhì)量為me的二次蒸汽所需要的換熱量為：

式中Φ——蒸發(fā)出二次蒸汽所需要的換熱量,kW；

me——二次蒸汽蒸發(fā)量，kg/s；

Δhv——二次蒸汽潛熱，kJ/kg。

在蒸發(fā)器中，傳熱系數(shù)的選擇和計(jì)算較為復(fù)雜，影響因素較多，如換熱器結(jié)垢、管壁磨損等，因此傳熱系數(shù)的大小對(duì)系統(tǒng)設(shè)計(jì)有較大影響。聯(lián)列式（1）與式（2），可得到壓縮機(jī)與蒸發(fā)器聯(lián)合運(yùn)行的能量方程：

式（3）表明，蒸發(fā)量與換熱器的換熱面積以及傳熱系數(shù)均成正比，與換熱溫差也呈線性關(guān)系。

1.2.2 壓縮比與換熱溫差的關(guān)系

以水提取的藥液為研究對(duì)象，在40～80℃溫度范圍內(nèi)，水蒸氣的飽和壓力溫度擬合曲線方程如下：

式中p——水蒸氣飽和壓力，kPa；

t——水蒸氣飽和溫度，℃。由于模型中忽略系統(tǒng)內(nèi)的能量損失，因此原液蒸發(fā)溫度與壓縮機(jī)入口水蒸氣溫度近似相等（te=tin）。將水蒸氣從pA壓縮到pS，經(jīng)過噴水處理后，壓縮比與換熱溫差關(guān)系如下：

式中ε——壓縮比。

若蒸發(fā)器所需的換熱溫差越低，壓縮機(jī)提供的壓縮比也就越小，同時(shí)需要的壓縮功也就越低，對(duì)系統(tǒng)節(jié)能越有利。但是結(jié)合式（3）可知，換熱溫差減小后，為了提供足夠的換熱量，需增加換熱面積或提高傳熱系數(shù)，二者一定程度上受到技術(shù)和成本的限制，因此，換熱溫差的取值對(duì)MVR系統(tǒng)的運(yùn)行影響重大。

1.2.3 蒸汽壓縮機(jī)特性

系統(tǒng)所需的壓縮比是由蒸汽壓縮機(jī)提供的，壓縮機(jī)也是整個(gè)系統(tǒng)的主要耗能設(shè)備，壓縮機(jī)的運(yùn)行效果直接決定了系統(tǒng)的蒸發(fā)能力。某型號(hào)離心式蒸汽壓縮機(jī)在12 800、14 400、16 000 r/min三個(gè)轉(zhuǎn)速下的特性曲線如圖3所示，式（6）為壓縮比與流量的擬合關(guān)系式。

式中α、β、χ、δ——分別為擬合方程的系數(shù)；

ms、mc——分別為喘振點(diǎn)流量和堵塞點(diǎn)流量。

根據(jù)壓縮機(jī)的多變效率曲線，可以計(jì)算出蒸汽壓縮機(jī)的實(shí)際輸入功率：

式中Wλ——壓縮機(jī)輸入功率，kW；

hout——壓縮后水蒸氣焓值，kJ/kg；

hin——壓縮機(jī)入口水蒸氣焓值，kJ/kg；

ηpol——壓縮機(jī)多變效率；

ηm——傳動(dòng)效率，取0.75；

圖3 離心式蒸汽壓縮機(jī)特性曲線

ηs——電機(jī)效率，取0.9。

在MVR系統(tǒng)中，根據(jù)系統(tǒng)能量守恒、壓縮比與換熱溫差關(guān)系、蒸汽壓縮機(jī)特性曲線，可建立描述系統(tǒng)運(yùn)行的數(shù)學(xué)模型：

對(duì)于MVR系統(tǒng)而言，在確定蒸發(fā)溫度、壓縮機(jī)與蒸發(fā)器特性參數(shù)后，方程組封閉，可求解出系統(tǒng)在不同工況下的蒸發(fā)量、壓縮比、蒸發(fā)器換熱量、壓縮機(jī)輸入功率等核心參數(shù)。

2 實(shí)例分析

以某中藥制藥企業(yè)蒸發(fā)濃縮工藝節(jié)能改造項(xiàng)目為例，對(duì)新增設(shè)MVR系統(tǒng)進(jìn)行建模求解與測(cè)試。該系統(tǒng)的設(shè)計(jì)蒸發(fā)量為2880 kg/h，蒸發(fā)溫度為70℃，系統(tǒng)選用離心式蒸汽壓縮機(jī)（裝機(jī)功率110 kW，特性曲線已確定），同時(shí)選擇降膜式蒸發(fā)器，換熱面積為210 m2，設(shè)計(jì)傳熱系數(shù)為1.50 kW/(m2·K)。

2.1 模型校核

MVR系統(tǒng)在原液預(yù)熱、建立真空后，蒸發(fā)器中存在一定量的閃蒸蒸汽。啟動(dòng)壓縮機(jī)，隨著壓縮機(jī)轉(zhuǎn)速的提升，系統(tǒng)的蒸發(fā)量、壓縮比將不斷增加。利用數(shù)學(xué)模型求解出12 800、14 400、16 000 r/min三個(gè)轉(zhuǎn)速下的蒸發(fā)量、壓縮比、輸入功率等參數(shù)的理論值，并與實(shí)際測(cè)量值進(jìn)行對(duì)比。

圖4中的曲線交點(diǎn)即為壓縮機(jī)與蒸發(fā)器共同工作狀態(tài)點(diǎn)，是三個(gè)轉(zhuǎn)速下壓縮機(jī)所需換熱量曲線與蒸發(fā)器提供的換熱量曲線的交點(diǎn)，也是數(shù)學(xué)模型的圖形解。若圖4中無曲線交點(diǎn)，則說明系統(tǒng)設(shè)計(jì)出現(xiàn)偏差，需重新選定蒸發(fā)器參數(shù)。求解出的計(jì)算值與實(shí)際測(cè)量值見表1。

從表1中的數(shù)據(jù)可以看出，由模型求解出的蒸發(fā)量、壓縮比的理論值與實(shí)測(cè)值誤差較小，而壓縮機(jī)輸入功率的誤差則較大。系統(tǒng)測(cè)試運(yùn)行中，壓縮機(jī)轉(zhuǎn)速提升至16 000 r/min，工況點(diǎn)A滿足設(shè)計(jì)要求，此時(shí)實(shí)測(cè)蒸發(fā)量為2900 kg/h，壓縮比為1.330，換熱量為1 876.9 kW，壓縮機(jī)電機(jī)輸入功率為106 kW。

圖4 壓縮機(jī)轉(zhuǎn)速對(duì)系統(tǒng)工況的影響

表1 各工況點(diǎn)計(jì)算值與測(cè)量值

造成誤差的主要原因有：（1）中藥提取液在循環(huán)泵的驅(qū)動(dòng)下形成泡沫，增加了傳熱熱阻，使傳熱系數(shù)低于設(shè)計(jì)值。（2）壓縮機(jī)總效率實(shí)際值低于理論值。

2.2 模型仿真

系統(tǒng)經(jīng)過一段時(shí)間運(yùn)行后，蒸發(fā)器的壁面污垢厚度增加，傳熱系數(shù)降低。壓縮機(jī)轉(zhuǎn)速恒定于16 000 r/min，若蒸發(fā)器換熱系數(shù)從1.50 kW/(m2·K)下降至0.74 kW/(m2·K)，應(yīng)用數(shù)學(xué)模型仿真出傳熱系數(shù)下降后的情況如圖5所示。

圖5 傳熱系數(shù)對(duì)蒸發(fā)量、壓縮比的影響

由圖5中曲線可以看出，當(dāng)傳熱系數(shù)下降后，蒸發(fā)量急劇下降，而壓縮比卻迅速上升，運(yùn)行工況點(diǎn)會(huì)從點(diǎn)A沿壓縮機(jī)特性曲線向右下方移動(dòng)，并接近壓縮機(jī)喘振點(diǎn)，系統(tǒng)運(yùn)行將極不穩(wěn)定。

由模型仿真結(jié)果可知，當(dāng)傳熱系數(shù)下降50.7%后，蒸發(fā)量理論值從3008kg/h下降至1 796.4 kg/h，壓縮比從1.318上升至1.395。在實(shí)際生產(chǎn)中，當(dāng)轉(zhuǎn)速達(dá)到設(shè)定的16 000 r/min，穩(wěn)定蒸發(fā)量小于2000 kg/h或壓縮比高于1.38后，應(yīng)考慮對(duì)蒸發(fā)器列管進(jìn)行清洗，以提高蒸發(fā)器傳熱系數(shù)，保證系統(tǒng)正常運(yùn)行。通過對(duì)模型的仿真，預(yù)測(cè)出蒸發(fā)器傳熱惡化后系統(tǒng)的蒸發(fā)量與壓縮比，對(duì)系統(tǒng)實(shí)際運(yùn)行具有一定指導(dǎo)意義。

3 結(jié)論

（1）針對(duì)中藥濃縮用機(jī)械蒸汽再壓縮系統(tǒng)，利用系統(tǒng)能量守恒式、壓縮比與換熱溫差關(guān)系式、蒸汽壓縮機(jī)特性曲線擬合式建立了系統(tǒng)的數(shù)學(xué)模型。

（2）以某中藥制藥企業(yè)的MVR系統(tǒng)為例，將數(shù)學(xué)模型求解值與實(shí)測(cè)值進(jìn)行對(duì)比，結(jié)果表明，該模型求解精度較高。傳熱系數(shù)與壓縮機(jī)總效率是造成模型誤差的主要因素。

（3）MVR系統(tǒng)測(cè)試運(yùn)行中，壓縮機(jī)轉(zhuǎn)速提升至16 000 r/min后，實(shí)測(cè)蒸發(fā)量為2 900 kg/h，壓縮比為1.330，換熱量為1 876.9 kW，壓縮機(jī)輸入功率為106 kW，滿足了設(shè)計(jì)要求。

（4）對(duì)數(shù)學(xué)模型進(jìn)行仿真，得到傳熱系數(shù)下降后蒸發(fā)量與壓縮比的預(yù)測(cè)值，對(duì)系統(tǒng)實(shí)際運(yùn)行具有一定指導(dǎo)意義。

[1]中華人民共和國國家統(tǒng)計(jì)局.中國統(tǒng)計(jì)年鑒 [M].北京：中國統(tǒng)計(jì)出版社，2015.

[2]田榮海，陳國禮，吳俊生，等.五效蒸發(fā)裝置模型與仿真 [J].華東理工大學(xué)學(xué)報(bào)，2000，26（5）：568-572.

[3]修允，陳嵐，劉瑜鳳，等.藥用蒸餾水機(jī)基本傳熱結(jié)構(gòu)性能分析 [J].中國醫(yī)藥工業(yè)雜志，2014，45（6）：577-581.

[4]顧承真，洪厚勝，張志強(qiáng)，等.羅茨壓縮機(jī)驅(qū)動(dòng)MVR熱泵系統(tǒng)的實(shí)驗(yàn)研究 [J].化工進(jìn)展，2015，34（6）：1602-1606.

[5]龐衛(wèi)科，林文舉，潘麒麟，等.離心風(fēng)機(jī)驅(qū)動(dòng)機(jī)械蒸汽再壓縮熱泵系統(tǒng)的性能分析 [J].機(jī)械工程學(xué)報(bào)，2013，49（12）：142-146.

[6]王力威，莊景發(fā)，楊魯偉，等.單螺桿水蒸氣壓縮機(jī)驅(qū)動(dòng)的MVR系統(tǒng)性能實(shí)驗(yàn)研究 [J].中國科學(xué)院大學(xué)學(xué)報(bào)，2015，32（1）：38-45.

[7]李清方，劉中良，龐會(huì)中，等.基于機(jī)械蒸汽壓縮蒸發(fā)的油田污水脫鹽系統(tǒng)及分析 [J].化工學(xué)報(bào)，2011，62（7）：1963-1969.

[8]唐文騫，楊同蓮，宋冬寶.機(jī)械蒸汽再壓縮技術(shù)在鈦白黑鈦液濃縮中應(yīng)用[J].化學(xué)工程，2015，43（9）：21-24.

Model and Simulation of Mechanical Vapor Recompression System for Concentrated Traditional Chinese Medicine

Tang Tianjun Liu Chenggang Xiao Cong Li Weihua

Aiming at the mechanical vapor recompression(MVR) system for concentrating traditional Chinese medicine,the mathematical model of the system is established by means of the formula of energy conservation,compression ratio and heat transfer temperature difference,and the fitting of the characteristic curve of the steam compressor.Taking an MVR system of a Chinese pharmaceutical enterprise as an example,the mathematical model is compared with the actual measured value,and it is found that the accuracy of the model is higher.The actual heat transfer coefficient and the total compressor efficiency are the main factors causing model errors.The model is simulated,and then the predicted value of evaporation and compression ratio is obtained after the evaporator heat transfer coefficient is decreased,which has certain significance for guiding production.

Chinese traditionalmedicine concentrate;Mechanicalvaporrecompression;Modeling;Simulation;Evaporator;Heat transfer coefficient;Compression ratio

TK 11+5

10.16759/j.cnki.issn.1007-7251.2017.06.006

2016-10-10)

*鎮(zhèn)江市經(jīng)信類專項(xiàng)資金扶持節(jié)能與循環(huán)經(jīng)濟(jì)項(xiàng)目（鎮(zhèn)經(jīng)信[2015]103號(hào)）。

**湯添鈞，男，1992年生，碩士。蘇州市，215009。