1 000 MW機組冷卻塔加裝導風通道的方法研究

周蘭欣，郝 穎

（華北電力大學能源動力與機械工程學院，河北保定071003）

1 000 MW機組冷卻塔加裝導風通道的方法研究

周蘭欣，郝 穎

（華北電力大學能源動力與機械工程學院，河北保定071003）

以冷卻塔相關傳熱傳質理論為基礎，應用計算機模擬軟件，建立了1 000 MW機組自然通風濕式冷卻塔的三維模型。針對冷卻塔中心區域換熱效率低的問題，提出了在雨區加裝導風通道的方案，通過導風通道，將塔外冷空氣直接導入塔內中心區域進行換熱，并對不同導風通道的結構尺寸及個數進行了計算和模擬。計算結果表明，加裝優化后的導風通道，可降低冷卻塔的出塔水溫。

機組；冷卻塔；導風；通道；模擬；計算；加裝；方案

0　概述

自然通風濕式冷卻塔是發電廠冷端系統中重要的熱力設備［1］。評價冷卻塔工作性能的主要指標就是出塔水溫，經驗表明：水溫每降低1℃，凝汽器真空度提高1%［2］。

目前，許多冷卻水塔的工作流程，是讓環境冷空氣進入冷卻塔雨區后，沿徑向向塔中心區域流動，在流動過程中，由于與水滴進行熱質交換，使空氣溫度升高，濕度增大，流速降低。現提出在雨區加裝導風通道的方法，簡單易行，可使環境冷空氣經導風通道，直接到達雨區的中心區域，改善中心區域的換熱狀況。通過對1 000 MW機組自然通風濕式冷卻塔進行模擬計算，結果表明，出塔水溫可降低0.66℃。

1　計算模型

對自然通風濕式冷卻塔進行數值模擬，通常根據冷卻水傳熱傳質特點，對噴淋區、填料區、雨區三個區域進行模擬計算。在雨區和噴淋區采用離散相模型模擬氣水間的相互作用，其中空氣為連續相，水滴為離散相。填料區的氣水間傳熱傳質以及阻力計算，采用自定義函數添加源項的方法進行求解計算［3—4］。

1.1 連續相（濕空氣）控制方程

描述濕空氣的運動狀態控制方程為質量、動量、能量及組分守恒方程。方程的通用形式［5］：

式（1）中：ρ為空氣密度；ui為速度矢量；φ為通用變量，分別表示各向速度分量（u、v、w）、水蒸氣組分Y、溫度T、湍動能k和湍流耗散率ε；Γφ為廣義擴散系數；Sφ為廣義擴散源項。

1.2 離散相（水滴）控制方程

冷卻塔內的冷卻水以水滴形式下落，與空氣進行傳熱傳質，使用離散項模型，模擬和追蹤水滴運動軌跡，通過水滴與空氣間耦合計算熱質傳遞［5］。水滴的溫度變化為：

式（2）中：mp為液滴質量；Tp為液滴溫度；T∞為濕空氣溫度；hfg為液滴的汽化潛熱，h為液滴表面傳熱系數。

水滴的蒸發速率為：

式（3）中：Ap為氣液接觸面積；hd為表面傳質系數；Cs為液滴表面蒸汽摩爾溶度；C∞為濕空氣中蒸汽摩爾溶度；Mw為液滴摩爾質量。

1.3 填料區傳熱傳質計算

填料區是冷卻塔內進行熱質交換的主要區域，作用是增加水滴和空氣的接觸時間和換熱表面積。填料區結構緊湊且數量眾多，精確模擬較難實現。因此，對填料區采用用戶自定義函數添加源項的方法進行模擬［6］。

根據填料區結構特點，對該區域進行離散化，如圖1所示。填料區徑向離散為m個柱體，沿垂直方向離散為n層，整個區由此離散為m×n×1個網格，每個網格的循環水流動和換熱情況，可用溫度和質量流率進行描述。計算所得參數進入下一網格進行引導計算，以此計算后，得出填料區內的全部換熱量。

計算填料區域內的傳熱傳質情況時，以Poppe模型作為數學基礎進行分析，單位高度的傳質系數可由淋水密度和空氣流率得出［7-8］。通過實驗數據得出1.2m高度填料的傳質系數表達式為：

圖1　填料內熱質傳遞簡化模型

式（4）中：Mep為Merkel數；Lfill為填料高度；A為汽水接觸面積；mw為淋水密度；ma為空氣質量流率。傳熱系數與傳質系數的關系由劉易斯關系式確定［］：

劉易斯數Lef可由Bosjnakovics公式得到［7］。

基于以上理論，得出第n層出口循環水質量流率為：

第n層的循環水出口溫度為：

當空氣達到飽和狀態時，在第n層填料控制單元中，單位體積內的附加源項可通過式（8）描述，它通過用戶自定義函數添加源項至控制方程來實現［9］。

式（8）、式（9）中：ΔLfi為當前填料層的高度，Cpw和Cpv分別為水的比熱和飽和水蒸氣的比熱，ΔTnw、分別為第n層填料內水溫變化量、水蒸汽冷凝量、質量源項和能量源項。

1.4 阻力計算

在冷卻塔中包括了氣水兩相間的作用力，可通過離散項模型進行模擬。此外，還包括了結構阻力，該部分阻力，可通過自定義函數添加動量源項，利用控制方程的形式進行求解計算。

動量方程的一般形式為［10］：

式（10）中：Vp為空氣通過幾何邊界面的垂直速度分量；K為壓損系數，它是一個通過實驗數據擬合得出的經驗關系式［5，11］，表達式為：

式（11）中：mw、ma分別為水滴、濕空氣的質量流量；L為填料高度。進風口、配水管路和收水器的壓損系數，根據實驗數據確定［12］，分別為Kjfk=0.5、Kgw=0.5和Kssq=3.5。

2　冷卻塔熱力特性模擬

2.1 幾何模型及邊界條件

以某1 000 MW機組為例，塔高165 m，頂部直徑為80.08 m，進風口高度11.64 m，進風口上緣直徑為123.942 m，進風口下緣直徑為133.33 m，填料層直徑為123.62 m，喉部直徑為75.21 m，填料厚度為1.2 m，淋水面積為12 000 m2。考慮到周圍環境對冷卻塔換熱的影響，建立了直徑為500 m、高度為500 m的圓柱體計算區域，如圖2所示。利用GAMBIT軟件生成幾何模型，采用四面體和楔形網格混合的方式，劃分計算區域網格。對該模型進行網格無關性驗證，最終確定網格總數為130萬。

設置了物理邊界條件，將環境區域外圍側面設為壓力入口，環境區域頂面設為壓力出口。離散相邊界條件為：將水池面、地面設為逃逸面，冷卻塔內壁設為反射面。

圖2　計算邊界

2.2 模擬結果分析

參考工況參數：進塔水量為21 503 kg/s，空氣干球溫度為29.96℃，相對濕度為70.3%，進塔水溫為42.3℃，水滴當量直徑為3 mm。

經計算得到冷卻塔X—Y截面濕空氣溫度等值線分布圖以及氣流速度等值線分布圖，如圖3、圖4所示。

圖3　冷卻塔X-Y截面濕空氣溫度分布

圖4　冷卻塔X-Y截面濕空氣速度分布

從圖3、圖4可知，冷卻塔內的氣水流場為軸對稱情況，在雨區進風口處的空氣流速較高，溫度最低。環境冷空氣從冷卻塔外圍向中心區域流動，在此過程中，冷空氣與水滴進行熱質交換，流速降低，隨著空氣冷卻能力的下降，使得空氣溫度沿徑向逐漸升高。空氣濕度增大，并在中心區域處達到最高值。通過對離散項水滴進行跟蹤計算，最終得到出塔水溫為30.62℃。

3　加裝導風通道

3.1 模型概述

針對冷卻塔中心區域換熱效率較低的問題，提出加裝導風通道的方案，管道一端布置在進風口處，另一端沿徑向通向塔中心區域，通道內的新鮮空氣未經過換熱，溫度和濕度都較低，具有較好的冷卻性能。導風通道將空氣直接引入塔中心區域進行換熱，可提高中心區域的換熱效果。初步設置了6個導風通道，長度為50 m，管道截面高5 m、寬5 m。導風通道的結構布置，如圖5所示。

圖5　冷卻塔加裝導風通道

3.2 模擬結果分析

加裝導風通道后，X—Y截面空氣溫度分布和氣流速度分布，如圖6、圖7所示。相比于原型塔的X—Y截面分布圖，塔中心區域的空氣溫度降低約2℃，塔內空氣流速更為均勻。

圖6　加裝導風通道后X-Y截面濕空氣溫度分布

圖8、圖9分別為原型塔與加入導風通道后的填料底部橫截面溫度分布。

從圖8、圖9可知，加裝導風通道后，塔內中心區域的溫度有所下降，靠近通道出口附近的溫度下降較為明顯，溫度分布更為均勻，中心區域降溫的最大值約4.5℃，冷卻效果明顯。

圖7　加裝導風通道后X-Y截面濕空氣速度分布

圖8　未裝導風通道填料底面溫度分布

圖9　加裝導風通道后填料底面溫度分布

3.3 導風通道結構優化

3.3.1 導風通道結構尺寸優化

進入塔中心區域的空氣量與導風通道的橫截面積有關，同時，也會影響冷卻塔的冷卻效果及出塔水溫。因此，選擇合理的導風通道尺寸，可使加裝導風通道冷卻塔的換熱效果達到最優。設置管道個數為6個，確定管道高度為5 m、分別取管道寬為2.5 m、5 m、7.5 m、10 m、12.5 m，管道長為30 m、35 m、40 m、45 m、50 m進行模擬。模擬計算的結果，如表1所示。

表1　不同尺寸導風通道模擬結果比較℃

由表1可知，當導風通道的寬度為10 m、長度為45 m時，出塔水溫最低，比不加導風通道時下降了約0.48℃。

3.3.2 導風通道個數優化

導風通道的數量同樣影響進入塔中心區域的空氣量，數量過少，從導風通道進入的空氣量就較少，使換熱效果不明顯，數量過多，將影響周圍雨區的換熱。因此，確定管道高為5 m，取管道寬為2.5 m、5 m、7.5 m、10 m、12.5 m，管道長為30 m、35 m、40 m、45 m、50 m，管道個數分別設為8個、10個、12個，分別進行模擬。模擬的結果，如表2、表3、表4所示。

表2　8個導風管道的模擬結果℃

表3　10個導風管道的模擬結果℃

表4　12個導風管道的模擬結果℃

由表2～表4可知，當導風通道的長度為45 m、通道個數為8個時，出塔水溫最低，比不加導風通道時下降了約0.53℃。

3.3.3 導風通道長度配置優化

為了使換熱更加均勻，得到更好的換熱效果，可選取不同長度的管道進行混合布置，使塔內不同半徑上的空氣質量均有所改善，混合布置方案，如表5所示。

表5　導風通道長度混合布置方案

當導風管的長度取35 m、45 m進行混合布置，且管道數為10個時，出塔水溫較低。加裝優化后的導風通道的模擬結果，如圖10、圖11所示。

圖10　X-Y截面濕空氣溫度分布

圖11　填料底面溫度分布

通過模擬可知，選取導風通道的管長為35 m、45 m進行混合布置，出塔水溫為最低，比長度相同導風管降低了0.13℃。因此，采用不同長度的導風通道布置，可獲得較低的出塔水溫。

4　結語

建立了比較完整的自然通風濕式冷卻塔模型，分析計算了1 000 MW機組的熱力特性，提出了在雨區加裝導風通道的方案，以提高中心區域換熱效率。同時，對導風通道的主要參數進行了優化。

（1）通過加裝導風通道，直接將低溫空氣引進冷卻塔中心區域，模擬分析了改進后的冷卻塔熱力特性，得到了最優結果，使出塔水溫下降了0.53℃。

（2）采用不同長度的導風通道布置，以改善冷卻塔不同半徑處的換熱效果，使換熱更均勻。計算結果表明，選用通風管道長度為35 m、45 m進行混合布置，可使冷卻效果最佳，進一步降低了出塔水溫約0.13℃。

［1］趙振國.冷卻塔［M］.北京：中國水利水電出版社，1997.

［2］李秀云，林萬超.冷卻塔的節能潛力分析［J］.中國電力，1997，30 （10）：34-36.

［3］趙元賓.側風對于自然通風逆流濕式冷卻塔傳熱傳質影響機制的研究［D］.山東大學.2009.

［4］Al-waked R，Behina M.CFD simulation of wet cooling tower ［J］.Applied Thermal Engineering，2006，（26）：382-395.

［5］N.Williamson，S.Armfield，M.Behnia.Numerical simulation of flow in a natural draft wet cooling tower-The effect of radial thermo fluid fields［J］.Applied Thermal Engineering，2008，28 （2）：178-189.

［6］FLUENT：User's Guide［M］.Fluent Inc，2003.

［7］N.Williamson，S.Armfield，M.Behnia.Numerical simulation of flow in a natural draft wet cooling tower-The effect of radial thermo fluid fields［J］.Applied Thermal Engineering，2008，28（2）：178-189.

［8］J.C.Kloppers.A Critical Evaluation and Refinement of the Performance of Wet-cooling Towers［D］.South Africa：University of Stellenbosch，2003.

［9］Y.Fourneir，V.Boyer.Improvements to the N3S-AERO heat exchanger and cooling tower simulation code［J］.The 12th IARH Symposium in Cooling Tower and Heat Exchangers，Sydney，Australia，2001（1）：31-38.

［10］趙順安.海水冷卻塔［M］.北京：中國水利出版社，2007.

［11］J.C.Kloppers，D.G.Kroger.The Lewis factor and its influence on the performance prediction of wet-cooling tower ［J］.International Journal of Thermal Sciences，2005，44（9）：879-884.

［12］周蘭欣，蔣波，陳素敏，自然通風濕式冷卻塔熱力特性數值模擬［J］.水利學報，2009，40（2）：208-213.

Method Study on the Installation of Catheter Channels in the Cooling Tower of a 1 000 MW Power Unit

ZHOU Lan-xin，HAO Ying
（School of Energy Power and Mechanical Engineering，North China Electric Power University，Baoding 071003，Hebei，China）

Based on the relative theories of heat and mass transfer in the cooling tower，a three-dimension model of the natural draft wet cooling tower of a 1 000 MW power unit has been built with computer simulation software. Considering the low heat transfer efficiency of the central region，this paper proposes the solution of installing several catheter channels in the rain zone.The outer cold wet air can transfer heat in the central region directly through the catheter channels.This paper also simulates different structure size and number of catheter channels.Calculation shows that outlet water temperature drops with the installation of catheter channels.

unit；cooling tower；catheter；channel；simulation；calculation；installation；plan

TK264.1

A

1672-0210（2016）01-0015-06

2015-12-08

周蘭欣（1956-），男，教授，主要從事電廠冷端設備的節能與系統優化方面的研究工作。