分布式能源站冷卻塔冷卻性能影響因素研究

路亮亮 劉文韜 谷菁 趙大周

（1 石家莊華電供熱集團有限公司 河北石家莊 050000 2 華電電力科學(xué)研究院有限公司 浙江杭州 310030 3 浙江省蓄能與建筑節(jié)能技術(shù)重點實驗室 浙江杭州 310030）

1 概述

分布式能源站是指分布在用戶端的能源綜合利用系統(tǒng)，以天然氣、可再生能源等綠色資源，為用戶端提供熱、電、冷等多種能源［1］，實現(xiàn)能源梯級利用，有效地提高了能源利用的安全性和靈活性。分布式能源站多采用電制冷供能，冷卻塔［2］是影響電制冷效率的關(guān)鍵因素。冷卻塔是利用水與空氣流動接觸后進行冷熱交換產(chǎn)生蒸汽，蒸汽揮發(fā)帶走熱量達到蒸發(fā)散熱，來散去工作系統(tǒng)產(chǎn)生的余熱降低水溫的蒸發(fā)散熱裝置［3］。

冷卻塔是常用的散熱設(shè)備，按氣水接觸方式，可分為干式、濕式和干濕混合式［4］。濕式冷卻塔因具有換熱效率高、運行穩(wěn)定、維護費用低等優(yōu)點被廣泛應(yīng)用于電廠、制冷、空調(diào)等領(lǐng)域。濕冷塔中空氣充當(dāng)冷卻介質(zhì)，所以環(huán)境對濕冷塔的運行性能有著重要的影響，換熱水溫、空氣溫度、空氣濕度等因素［5］均會對換熱效果產(chǎn)生影響。除以上因素外，由于分布式能源站靠近用戶，在其設(shè)計布置時還要考慮其噪聲對周圍環(huán)境的影響，多采用設(shè)置隔音墻減少噪音，圍墻也成為換熱效果重要影響因素之一。

相比于實驗研究［5］，采用數(shù)值計算及計算機流體動力學(xué)模擬［6-7］具有高效、低成本和可實現(xiàn)多工況分析的特點。計算流體動力學(xué)CFD 目前已被廣泛應(yīng)用。FLUENT 是1 種較為典型的CFD 商用軟件，可模擬流體、熱傳遞和化學(xué)反應(yīng)等，具有豐富的物理模型、先進的數(shù)值方法和強大的前后處理功能。

本文采用數(shù)值計算研究空氣含濕量對冷卻塔冷卻效果的影響，采用Fluent 軟件模擬冷卻塔周圍流場，以確定圍墻對冷卻塔的影響，為分布式能源站工程應(yīng)用奠定理論基礎(chǔ)。

2 數(shù)值計算

分布式能源站冷卻塔為水和空氣直接接觸換熱，很多因素均會對換熱效果產(chǎn)生影響。換熱后，空氣含濕量會發(fā)生變化。按《冷卻塔驗收測試規(guī)程》（T/CECS 118—2017）［8］，冷卻塔的冷卻效果可用冷卻數(shù)來衡量，冷卻數(shù)是表征冷卻塔內(nèi)包括噴淋區(qū)、淋水填料區(qū)和填料下的雨區(qū)在內(nèi)的臨水裝置的熱、質(zhì)交換能力的特征數(shù)［6］。本文采用相同冷卻數(shù)條件下計算不同空氣含濕量得到的出塔水溫，以研究空氣含濕量冷卻塔冷卻效果的影響，具體見公式（1）。

式中：Ωh為冷卻數(shù)；CW為水的比熱；Δt 為溫差；K 為蒸發(fā)冷卻帶走熱量系數(shù)；Δhmh為平均焓差。

蒸發(fā)冷卻帶走熱量系數(shù)計算式見式（2）。

式中：t2為出塔水溫。

平均焓差計算式見式（3）。

式中：Xs根據(jù)η 和ξ 查Xs計算圖求得，η、ξ 和δh″根據(jù)式（4）～（6）計算。

式中：h1為進塔空氣焓；h1″為空氣溫度等于進塔水溫時飽和空氣焓；h2″為空氣溫度等于出塔水溫時飽和空氣焓；λ 為氣水比。

干球溫度一定條件下，通過焓濕圖可得到不同空氣含濕量下的比焓，進而可計算得到不同出塔溫度。本文研究空氣含濕量這一影響因素對冷卻效果的影響，取空氣含濕量在0.5～1.5 之間的相等間隔數(shù)據(jù)，得到比焓，結(jié)果見表1，通過公式（1）～（7）計算得到出塔溫度，結(jié)果見圖1。

圖1 不同比焓下出塔溫度和冷卻數(shù)的關(guān)系

表1 不同空氣組分?jǐn)?shù)下比焓值

可以看出，空氣含濕量越大，比焓越大；相同入口水溫條件下，出塔溫度隨空氣含濕量增大而增大，說明冷卻效果變差。這是由于在環(huán)境空氣中，入口空氣濕度越大，水分含量越多，在熱交換的過程中，所吸收的水分和熱量就會越少，因此才會降低冷卻的效率。

3 流場模擬

3.1 模型構(gòu)建

3.1.1 幾何模型

本文對石家莊某分布式能源站冷卻塔開展數(shù)值模擬，簡化后建立的計算模型如圖2 所示，網(wǎng)格分布如圖3 所示。冷卻塔周圍空氣為湍流流動，模型充分考慮空氣擴散域，其中最內(nèi)側(cè)的黃色矩形為電制冷冷卻塔，外側(cè)的黃色為圍墻墻體，最外側(cè)為氣流流動的流體域。冷卻塔部分尺寸為35 m×15 m×6 m，圍墻內(nèi)流體域尺寸為66 m×46 m×8m，外圍流體域尺寸為166 m×96 m×38 m。

圖2 計算模型

圖3 網(wǎng)格分布

3.1.2 邊界條件

入口條件采用速度入口，速度分別設(shè)為2.78 m/s，出口條件采用壓力出口。壁面設(shè)置為無滑移邊界條件。在定義所有邊界條件時必須給定溫度和水蒸汽質(zhì)量分?jǐn)?shù)。

3.1.3 計算模型

流體運動以雷諾數(shù)Re 的大小決定，分為湍流與層流。本文流體為湍流運動，在Fluent 中有很多湍流計算方法，本文采用標(biāo)準(zhǔn)k-epsilon 模型進行N-S 方程求解。模擬采用穩(wěn)態(tài)壓力基處理器，采用組分模型，設(shè)定為濕空氣。求解方法采用Simple 方法。

3.2 存在圍墻模擬分析

3.2.1 南北向風(fēng)力模擬

研究夏季南北向風(fēng)力2.78 m/s，氣溫32 ℃時的氣流運行情況，模擬結(jié)果如圖4 和圖5 所示。

圖4 南側(cè)進風(fēng)速度分布圖

圖5 南側(cè)進風(fēng)流線分布圖

選取典型截面進行研究，可以看出，當(dāng)南北風(fēng)吹來時，遇到墻體氣流向上運動，墻后側(cè)形成低速區(qū)（動力陰影區(qū)［9］），并在冷卻塔的左上方形成旋渦，導(dǎo)致冷卻塔上方的高溫氣體不易擴散，影響散熱效果。

動力陰影區(qū)影響高度計算公式［10］，見式（7）。

式中：Hk為動力陰影區(qū)的影響高度；A 為迎風(fēng)面的面積。

計算可得Hk=5.75 m，略小于圍墻高度6 m，對冷卻塔周圍流體擴散形成不良效果。

3.2.2 東西向風(fēng)力模擬

研究東西向風(fēng)力2.78 m/s，氣溫32 ℃時，冷卻塔周圍的氣流分布情況。模擬結(jié)果如圖6 和圖7 所示。

圖6 東側(cè)進風(fēng)速度分布圖

圖7 東側(cè)進風(fēng)流線分布圖

可以看出，當(dāng)風(fēng)從東側(cè)吹來，在圍墻與冷卻塔東側(cè)、冷卻塔頂部、冷卻塔與圍墻西側(cè)形成了空氣旋流，西側(cè)的旋流由冷卻塔頂部形成，一定程度上影響了冷卻塔的冷卻效果。

東西側(cè)進風(fēng)時，圍墻內(nèi)仍然為低速區(qū)，同時冷卻塔頂部存在“動力陰影區(qū)”，根據(jù)動力陰影區(qū)影響高度計算公式計算得到Hk=6.89 m。即東西向進風(fēng)由于迎風(fēng)面積大于南北方向進風(fēng)，因此動力陰影區(qū)的影響范圍大于南北向。

同時冷卻塔為東西兩側(cè)橫向進風(fēng)，冷卻塔頂部排出的高濕空氣受動力陰影區(qū)影響會再次從下方進入冷卻塔，空氣濕度聚集增加，冷卻效果進一步降低。

3.3 去掉圍墻后模擬

當(dāng)去掉圍墻后，研究夏季南北向風(fēng)力2.78 m/s、氣溫32 ℃時的氣流運行情況，模擬結(jié)果如圖8 和圖9 所示。

圖8 去掉圍墻后南側(cè)進風(fēng)速度分布圖

可以看出，當(dāng)去掉圍墻后，空氣碰到冷卻塔后仍會在上方產(chǎn)生渦流，但渦流影響范圍小了很多，僅在冷卻塔下風(fēng)向形成“動力陰影區(qū)”。根據(jù)動力陰影區(qū)影響高度計算公式計算得到Hk=2.85 m，遠(yuǎn)小于有圍墻時的動力陰影區(qū)大小，因此較有圍墻時，散熱效果得到了加強。

4 結(jié)論

本文分別采用數(shù)值計算和Fluent 軟件模擬研究空氣濕度和設(shè)置圍墻對分布式能源站冷卻塔冷卻性能的影響，得出如下2 個結(jié)論：

（1）相同冷卻數(shù)下，出塔溫度隨空氣含濕量增大而增大，空氣含濕量的增加將不利于冷卻塔散熱，冷卻效果更差。

（2）由于“去工業(yè)化”增加的圍墻影響了氣流正常運動，圍墻后及冷卻塔上方產(chǎn)生渦流，動力陰影區(qū)是未設(shè)置圍墻的2 倍以上，且東西側(cè)方向的動力陰影區(qū)大于南北側(cè)，動力陰影區(qū)越大，冷卻效果越差，熱量無法有效擴散出去。