許文柱，周文杰*，郝 帥，邵佳棟，張 進，趙 陽

（1.杭州電子科技大學能量利用與自動化研究所，杭州 310018；2.浙江大學制冷與低溫研究所，杭州 310027）

0 引言

隨著全球溫室效應的加劇，國際社會對減少溫室氣體排放的重視程度日漸提高，載冷劑在工業(yè)和商業(yè)制冷中的作用越來越被重視。在此背景下，冰漿作為一種優(yōu)良的載冷劑在制冷中的作用越發(fā)凸顯。冰漿的工業(yè)應用約有45年的歷史[1]。在冰漿系統(tǒng)的應用中，冰漿輸送的能耗是須重點考慮的問題，其中泵功更是重中之重。

一個典型的冰漿系統(tǒng)，如日本CAPCOM公司本部大樓[2]的冰漿制冷空調系統(tǒng)，其能耗占建筑總能耗的4%，其中，由于樓層很高造成的輸送損耗尤為巨大。因此，研究不同冰濃度冰漿的輸送泵功，從而提高冰漿系統(tǒng)的泵送效率可以顯著地節(jié)省冰漿的輸送成本。Frei等[3]研究發(fā)現(xiàn)，泵的性能隨著冰濃度增大而下降。張敬齋等[4]的研究也發(fā)現(xiàn)，相變顆粒會在一定程度上損壞泵。Hansen等[5]的研究認為，冰的蓄冷能力隨冰濃度增大而提高。N?rgaarda等[6]通過檢測泵的進出口壓差，得到了泵運行中的實際揚程，再通過檢測泵流量，得到了不同冰濃度下泵的性能曲線。徐寶成等[7]的研究發(fā)現(xiàn)，根據(jù)流體的不同流量調節(jié)泵的輸入功率可以減小能耗?；菸臉s[8]發(fā)現(xiàn)，在濕球溫度為16℃左右，負載37%左右，近5 h時間段內，與舊的控制邏輯相比，變頻調節(jié)的冷卻泵節(jié)能率為13%。劉超等[9]根據(jù)油的性質設置變頻調節(jié)技術的思路，按照冰漿對泵影響的不同設置變頻調節(jié)技術即可實現(xiàn)工程上的節(jié)能。冰漿對泵功的影響不是簡單的梯度問題。梁坤峰等[10]的研究發(fā)現(xiàn)，泵中冰濃度的分布是不均勻的。本文參考N?rgaarda的做法，使用進化策略算法，將不同流態(tài)冰濃度下的最優(yōu)泵功視為黑箱優(yōu)化問題進行研究。從經(jīng)濟效益的角度研究一個典型的冰漿冷卻系統(tǒng)，尋找其最優(yōu)泵功，為冰漿系統(tǒng)的運行提供建議。

1 冰漿制冷系統(tǒng)案例

CAPCOM大樓位于日本大阪中央?yún)^(qū)，是CAPCOM的研發(fā)部門，共20層。公司辦公室位于2～16層，每層面積570 m2。據(jù)計算，由于辦公室內的計算機產(chǎn)生熱量，使得辦公室的熱負荷為151 W/m2，高于大樓的平均水平。因此，必須采用氣體驅動的吸收式制冷系統(tǒng)為周邊區(qū)域提供冷卻。樓宇辦公間內部區(qū)域的冷卻由一個冰漿冷卻系統(tǒng)提供。該系統(tǒng)由兩個制冰機組組成，每個機組的冷卻能力為272 kW。位于17層的冰漿機產(chǎn)生的冰漿儲存在兩個儲罐中。冰漿的低溫使得設計者選用12℃而不是15℃的冷空氣，在日本這是正常設計。這將對空氣流量的要求從41 m3/h降到了32 m3/h，使得設備尺寸更小，空氣分配的電力需求更低，節(jié)省了資本和運營成本。冰漿冷卻系統(tǒng)的工作過程為：不需要或少需要冰漿時，制冰機工作產(chǎn)生的冰漿被泵送到儲藏室中；當有冷卻需求時，冰漿被從17層的儲藏室中取出，經(jīng)單管運輸系統(tǒng)被泵送到2～16層中每一層的空氣處理機組中，經(jīng)過冰漿盤管后，沒有用完的冰漿又進入回流通道，從本層泵送回17層的儲藏室。通過對上述系統(tǒng)的分析可知，長達16層的輸送管道使得輸送冰漿的耗能十分巨大，除了產(chǎn)生冰漿的制冰機耗能以及冰漿分離器耗能外，輸送冰漿的耗能即泵功是系統(tǒng)主要耗能，因此研究如何合理地減小冰漿系統(tǒng)的泵功具有重要意義。

2 冰漿泵送模型與算法設計

冰漿系統(tǒng)的單個泵簡化模型如圖1所示，模型是由儲罐、運輸泵、測試泵和控制閥組成的簡單冰漿輸送管道。冰漿儲存在儲罐中，運輸泵指被研究泵以外的所有泵，測試泵即為被研究泵。圖中的兩個流量計和兩個壓力傳感器用于監(jiān)控。實驗所用的泵為Grundfos CR2-50離心泵，溶液為16%質量分數(shù)的丙二醇水溶液，冰顆粒直徑為0.1～0.3 mm，泵轉速為2 800 r/min。從N?rgaarda等[6]的實驗圖中取點，所得數(shù)據(jù)分別如表1和表2所列。使用MATLAB2018a軟件進行數(shù)據(jù)擬合，得到與冰濃度相關的經(jīng)驗公式。

圖1 冰漿系統(tǒng)的單個泵簡化模型圖[6]Fig.1 Simplified drawing of ice slurry pumping model for single pump

表1 不同冰濃度下的泵流量－效率（Q-η）數(shù)據(jù)Tab.1 Pump efficiency data under different ice concentrations

表2 不同冰濃度下的泵流量－揚程（Q-H）數(shù)據(jù)Tab.2 Q-H data at different ice concentrations

冰濃度的定義公式為：

式中：C為冰漿中的冰濃度（質量分數(shù)），%；mice為冰顆粒的質量，kg；m為冰漿流體的質量，kg。

泵效率計算公式：

式中：η為泵效率，%；Pe為有效功率，W；P為泵的軸功率（即耗電量），W。有效功率也可以表達為：

將式（3）代入式（2）可得[11]：

式中：H為泵的總揚程，m；ρ為冰漿流體的密度，kg/m3；g為重力加速度，m/s2；V˙為實測體積流量，m3/s。冰漿流體的密度可以使用多相流密度公式計算[12]：

式中：ρice為冰密度，kg/m3；ρcs為多相流體密度，kg/m3。查詢ASHRAE手冊得到：16%丙二醇溶液的密度約為1 020 kg/m3，冰的密度約為917 kg/m3。

使用MATLAB的cftool工具對表2數(shù)據(jù)進行擬合，結果為：

式中：Q為流量（以m3/s為單位時，數(shù)據(jù)過小，為方便計算設為m3/h）。該經(jīng)驗公式只適用于0～30%的冰濃度。用MATLAB對式（6）繪圖，冰濃度為0～30%質量分數(shù)，流量為0～5 m3/h，得到圖2。

圖2 冰濃度與揚程的經(jīng)驗函數(shù)圖Fig.2 Empirical function diagram of ice concentration and head

用同樣的方式擬合處理表1冰濃度與效率數(shù)據(jù)，可得：

該經(jīng)驗公式適用的效率范圍為0～1，超過此范圍沒有意義。受限于泵的性能，超出泵可運送的流量與冰濃度的取值也無意義。用MATLAB對式（7）繪圖，冰濃度為0～30%，流量為0～5 m3/h，泵效率為0～100%，得到圖3。

圖3 冰濃度與泵效率經(jīng)驗函數(shù)圖Fig.3 Empirical function diagram of ice concentration and pump efficiency

將上面的經(jīng)驗揚程與效率式（6）與式（7）代入理論泵效率式（4），可得泵功如式（8）。用MATLAB對式（8）繪圖，冰濃度為0～30%，流量為0～5 m3/h，泵功為0～750 W，得到圖4。

圖4 冰濃度與泵功經(jīng)驗函數(shù)圖Fig.4 Empirical function diagram of ice concentration and pump power

張強等[13]使用典型高維復雜函數(shù)進行的測試表明，進化策略算法具有很好的收斂精度和計算速度。張鑫等[14]的思路是使用μ、λ進化策略算法進行優(yōu)化。算法流程如圖5所示。主要步驟如下：

圖5 算法流程圖Fig.5 Algorithm flow char

（1）首先，將求解最優(yōu)冰濃度的問題轉換為數(shù)學問題，即求解函數(shù)式（8）的最小值。將式（8）作為MATLAB算法程序的求解函數(shù)。

（2）其次，設置一個初始種群，這是第一代用于進化的親本（可隨機生成）。親本是由冰濃度與流量組成的矩陣，第一行為冰濃度，第二行為流量，一列為代表函數(shù)P的一組解。綜合考慮運行時間與結果的精確度，將初始種群規(guī)模設為10，即初始種群矩陣中有10組P函數(shù)的解。

（3）接著開始種群的進化。第一步，交叉親本。例如，將親本中的隨機兩組解“冰濃度10%、流量2 m3/h”與“冰濃度20%、流量1 m3/h”中的冰濃度或流量交換，變成“冰濃度10%、流量1 m3/h”與“冰濃度20%、流量2 m3/h”，使數(shù)據(jù)解的范圍更加廣泛，避免出現(xiàn)僅有低冰濃度對小流量，高冰濃度對大流量的情況發(fā)生。第二步，開始變異，給每一組解的兩個量都加上一個變量。該變量滿足某一零均值和某一方差的高斯分布，其中的方差在進化策略算法中被稱為變異程度。加上該變量等于給每組解增加了一個變異程度，可以表征這組解變異了多少。當種群進化到收斂時，變異程度也會慢慢減小，整個種群更加容易收斂?？偟膩碚f，交叉和變異都是為了使冰濃度與流量的值在允許的范圍內更加“隨機”，使得泵功最小的一組解出現(xiàn)的概率更高。

（4）最后，對完成變異的數(shù)據(jù)即子代進行選擇。一般來說，在μ、λ進化策略算法中，親本與子代的比例最好為1∶7。鑒于前文設定親本為10，因此子代應取70。將70個子代代入泵功公式計算，即可得到本次子代解對應的泵功矩陣。比較70組解的泵功矩陣，選出其中泵功最小的10組作為新的親本即第二代親本進行下一次進化。最終，經(jīng)過一代一代的進化，種群必然收斂于最優(yōu)解。

本質上，進化策略算法是一種隨機取值后計算比較適應度的概率性算法。其優(yōu)點是：冰漿最優(yōu)冰濃度是實際問題，其數(shù)學模型必然很復雜，用平常的方法難以分析。使用進化策略算法可以忽略其函數(shù)性質，很容易求出其最優(yōu)值，且避免了局部最優(yōu)的問題。隨著迭代次數(shù)的增長，最優(yōu)值更加精確，但運行時間也大幅加長。為應用最優(yōu)冰濃度算法，首先須定義“最優(yōu)”，并將“最優(yōu)”的標準作為進化策略算法里的適應度。求出最佳適應度的過程即為最優(yōu)化冰濃度的過程。

從工程實際出發(fā)，最優(yōu)的冰濃度即為經(jīng)濟效益最高的冰濃度，即在滿足同等蓄冷能力下，使得輸送泵功最小的冰濃度。表現(xiàn)在算法中即泵功為適應度，泵功越小，適應度越好（這是確定種群選擇的標準）。

張雪等[15]研究表明，冰顆粒相變過程能在很小的溫度（冰漿溫度）變化范圍內吸收大量的能量。設定必須滿足h能量的蓄冷能力，根據(jù)冰融化相變吸收的能量以及熱量傳遞守恒定律得：

式中：t為時間，s；hls為冰（冰顆粒）的潛熱，參考馮卓宏等[16]的實驗數(shù)據(jù)為335 kJ/kg；h為焓變，在等質量情況下，即為制冷量。焓變h越大，相變時所吸收的能量越多，冰漿的蓄冷能力越強。

即滿足式（10）時：

判斷為滿足供應需求。可以將式（10）作為算法的選擇條件。使用算法后得到最優(yōu)冰濃度。取h=1.388 9×104W（5×104kJ/h），泵的最高運行功率為750 W，通過算法運算，多次迭代得到，冰濃度為19.73%、流量為0.756 9 m3/h（2.102 5×10-4m3/s）時，泵功最小，為417.5 W，如圖6所示。

圖6 最優(yōu)算法運行結果曲線Fig.6 Running results of optimal algorithm

3 結果分析

由圖2可以看出，當冰濃度增加時，等流量下泵的揚程變化不大，當冰濃度達到30%的系統(tǒng)極限時，泵的流量僅在3.1 m3/h左右。這種情況出現(xiàn)的原因是泵的功率有極限，在一定轉速下，隨著冰濃度的增加，輸送單位流量的冰漿耗能更高。為了使冰漿能達到所需的揚程，只能犧牲流量，先輸送小部分冰漿。這里也表現(xiàn)出冰漿與普通水溶液的性質不同，顯然輸送清水的能耗遠比輸送冰漿的要低。由于泵本身的Q-H特性，流量必然隨著揚程的上升而下降，但在高冰濃度下，流量下降得更快。因此，總的來說，隨著冰濃度的增加，泵所能輸送的流量更小。

由圖3可以看出，當冰濃度增加時，泵的效率是一直下降的，也就是說，冰濃度越高，泵送越困難。流量上升，使泵效率先上升達到最大值后下降，當流量上升到泵所能承受的極值時，效率為0，此時的泵功無限大，無意義。

圖2、圖3的結果符合二相流對泵的影響規(guī)律[17]：當輸送的冰濃度增加時，揚程在小流量時下降較小，在大流量時下降幅度較大。隨著冰濃度的增加，泵的最大流量逐漸減小，為達到所需的揚程，功率就會增加，效率降低。這也說明該經(jīng)驗公式在一定范圍是可信的。

圖4左上角異常數(shù)據(jù)的出現(xiàn)是由于模型都是基于實驗的，僅考慮了其數(shù)學特性。實際上，當所需輸送的冰漿的冰濃度超過泵的輸送極限時，就無法泵送，表現(xiàn)在數(shù)學模型上就是泵功無限大。顯然，流量過大和冰濃度過大都會使泵功無限大。接近泵送極限時，泵功斷崖式下降，這表示模型在這部分的預測不準確。綜合流量與冰濃度的影響，本系統(tǒng)輸送高冰濃度的冰漿時流量不宜超過3.5 m3/h。泵的實際運行范圍如圖7所示，理論運行范圍如圖8所示。

圖7 泵的實際運行范圍圖Fig.7 Actual operation range of pump

圖8 泵的理論運行范圍圖Fig.8 Theoretical operating range of pump

出現(xiàn)泵的實際運行范圍與理論不符的原因是，在一定的流量以及冰濃度下（此時流量沒有超過泵的最大容量），泵效率盡管不為0（泵仍在以額定功率運行或超額定功率運行），但已經(jīng)很低了，達到泵送所需流量的泵功盡管沒有達到無窮大，但已經(jīng)超過了泵的極限功率（即斷崖式下降部分）。這種情況下，泵容易堵塞、損壞，甚至導致系統(tǒng)崩潰。因此，在算法設計上，必須避免這種情況的發(fā)生。

從圖6算法運行案例知道了制冷量為1.388 9×104W下的最優(yōu)冰濃度、流量配置與泵功。取不同制冷量運行算法，可以得到更多狀態(tài)下的最優(yōu)冰濃度、流量配置與泵功。利用變頻調節(jié)技術設置泵功輸入，用冰漿分離器（或制冰機）控制冰濃度，用閥門控制流量，即可實現(xiàn)最節(jié)能的冰漿泵送。應當注意的是，由于泵的功率存在上限，因此系統(tǒng)一段時間內可以供應的蓄冷能力h是有極限的，應用算法時，不應設置超越極限的蓄冷能力h值，否則會出現(xiàn)運行錯誤。

4 結論

本文以單個泵為研究對象，使用進化策略算法給出了一組較優(yōu)的冰漿泵送布置方案：Grundfos CR2-50離心泵額定轉速為2 800 r/min、輸送1.388 9×104W的制冷量時，冰濃度為19.73%、流量為2.102 5×10-4m3/s的冰漿流體可使運行泵功最小，節(jié)省運輸成本。本文的研究方法和所得結論可為國內建設類似冰漿系統(tǒng)提供借鑒。