步兆彬,夏友剛,孔凡新,江廣旭,楊繼沖,何鎖盈,高 明,張政清

(1.山東華聚能源股份有限公司 趙樓綜合利用電廠,山東 菏澤 274700;2.山東大學(xué) 能源與動力工程學(xué)院/高效節(jié)能及儲能技術(shù)與裝備山東省工程實驗室,山東 濟南 250061;3.山東交通學(xué)院 汽車工程學(xué)院,山東 濟南 250357)

0 引 言

自然通風(fēng)逆流式濕式冷卻塔因結(jié)構(gòu)簡單、冷卻效率高,常作為火力發(fā)電系統(tǒng)的冷卻設(shè)備[1-2],對其強化增效的研究也一直受到廣泛關(guān)注[3-4]。濕式冷卻塔的填料區(qū)占全塔換熱量的60%以上[5];而配水區(qū)的配水方式直接影響塔內(nèi)循環(huán)水的分布形態(tài),且對全塔氣水比的分布起到?jīng)Q定性作用[6]。目前,對填料區(qū)而言,非均勻布置是主要的強化措施,包括非等高[7]和非等片距布置[8]。相比于填料的非等高布置,填料的非等片距布置可以避免因填料高度不同而帶來的局部配水不均的問題,是更合理的一種填料布置方式[9]。SHAHALI等研究發(fā)現(xiàn)片距較小的填料熱力性能較好但通風(fēng)性能較差[10]。ZHANG等[11]、CHEN等[8,12]研究發(fā)現(xiàn)非等片距布置可有效強化冷卻塔的換熱,但需對填料的半徑分割點進行優(yōu)化設(shè)計才能達到較好的效果。

目前,分區(qū)配水方式是配水區(qū)最主要的性能強化方式[13-14],且二分區(qū)和三分區(qū)是最常用的分區(qū)方式。ZHOU等研究發(fā)現(xiàn),適當(dāng)減少內(nèi)區(qū)配水量同時增加外區(qū)配水量,可有效增強全塔換熱效果[15]。ZHANG等將分區(qū)配水和填料的非等片距布置進行耦合,形成了兩區(qū)協(xié)同增效方案[11]。結(jié)果表明,填料和配水兩區(qū)協(xié)同增效方案比單一方案具有更明顯的增效效果。

此外,冷卻塔雨區(qū)的水滴下落過程中沖擊水面形成較大的淋水噪聲,對周圍環(huán)境產(chǎn)生嚴(yán)重影響。對冷卻塔降噪的研究也一直受到研究人員的關(guān)注。周揚建立了淋水噪聲計算模型,研究了冷卻塔運行工況對淋水噪聲的影響規(guī)律[16]。黃平等對淋水噪聲進行了測試,發(fā)現(xiàn)濕式冷卻塔的淋水噪聲多為中高頻且穩(wěn)定的噪聲[17]。目前,冷卻塔的降噪措施主要是在冷卻塔周圍修建聲屏障,阻礙噪聲的傳播[18]。聲屏障的高度需要超過冷卻塔的進風(fēng)口,這勢必對冷卻塔的熱力和通風(fēng)性能造成影響[19]。無風(fēng)工況下聲屏障會影響冷卻塔的正常進風(fēng),而側(cè)風(fēng)工況下又一定程度上起到防風(fēng)墻的作用[20]。因此,聲屏障對冷卻塔性能的影響要視具體工況和運行條件而定,需深入研究。

隨著技術(shù)的進步,涉及填料和配水的兩區(qū)協(xié)同增效方案以及增設(shè)聲屏障用以降低冷卻塔噪聲的方式在冷卻塔的設(shè)計施工中已經(jīng)被廣泛應(yīng)用。但是,聲屏障對填料和配水協(xié)同優(yōu)化的濕式冷卻塔熱力及通風(fēng)性能影響的研究還未見報道。本文針對帶有聲屏障的兩區(qū)協(xié)同增效冷卻塔展開研究,基于數(shù)值計算的方法研究聲屏障對填料和配水協(xié)同優(yōu)化的濕式冷卻塔熱力及通風(fēng)性能影響規(guī)律,為帶有聲屏障的兩區(qū)協(xié)同增效冷卻塔的設(shè)計提供理論基礎(chǔ)和指導(dǎo)。

1 物理模型

1.1 幾何模型及物性參數(shù)

本文以魯西南地區(qū)某300 MW火電機組配備的濕式冷卻塔為研究對象。該塔總高度141.7 m,喉部高度113.02 m,出口直徑52.23 m,喉部直徑47.5 m,進風(fēng)口高度7.728 m,塔底直徑92.842 m。該塔的名義淋水面積為5 500 m2,設(shè)計填料為高度1 m的S波填料,設(shè)計配水系統(tǒng)采用均勻配水方式。該塔設(shè)計的運行參數(shù)及氣象數(shù)據(jù)如下:進塔水溫45 ℃,循環(huán)水流量39 186 m3/h,環(huán)境干球溫度28.3 ℃,相對濕度89%,大氣壓強99.6 kPa,累年平均風(fēng)速2.5 m/s。

為了降低雨區(qū)淋水噪聲對周圍辦公環(huán)境的影響,后續(xù)改造中,在辦公區(qū)方向增設(shè)2部分相互垂直的聲屏障(長度分別為105 m和54 m),聲屏障的高度為12.4 m。

為了增強該塔的熱力和通風(fēng)性能,對冷卻塔進行了涉及填料和配水的兩區(qū)協(xié)同增效改造。如圖1所示,冷卻塔內(nèi)區(qū)和外區(qū)分別采用30 mm和26 mm的填料,形成填料的非等片距布置方式,同時在總配水量一定的前提下增加外區(qū)配水、減少內(nèi)區(qū)配水,形成冷卻塔的分區(qū)配水模式。表1、2分別為所涉及填料的熱力和阻力參數(shù),其中A、m、A0和n0為性能實驗擬合參數(shù)。

圖 1 非等片距和分區(qū)配水協(xié)同增效Fig.1 Synergistic enhancement scheme of non-equidistant fillings and non- uniform water distribution

表 1 填料的熱力性能參數(shù)

表 2 填料的阻力性能參數(shù)

1.2 簡化假設(shè)

對冷卻塔的幾何結(jié)構(gòu)及其內(nèi)部氣水流動做如下簡化假設(shè),在此基礎(chǔ)上,對氣水流動過程進行建模計算。

1) 氣體為理想氣體,滿足氣體狀態(tài)方程;

2) 忽略塔內(nèi)細節(jié)結(jié)構(gòu),以阻力源項對其產(chǎn)生的氣動阻力進行計算;

3) 氣水流動皆處于穩(wěn)定狀態(tài);

4) 塔內(nèi)水滴為剛性球,只在z向運動。

2 數(shù)學(xué)模型及其離散化

2.1 控制方程

對濕空氣運動描述的控制方程為通用形式,并已經(jīng)有廣泛的應(yīng)用[5, 21],此處不再贅述。

根據(jù)守恒原理推導(dǎo)得出的循環(huán)水運動的控制方程如式(1)～(3)所示[5]。

(1)

(2)

(3)

式中:g為當(dāng)?shù)刂亓铀俣?m/s2;ρw為循環(huán)水的密度,kg/m3;ρ為濕空氣的密度,kg/m3;vw,z為循環(huán)水水滴在豎直方向上的速度,m/s;mw為單個水滴的質(zhì)量,kg;fz為單個水滴收到的來自氣流的阻力,N。

2.2 計算域及邊界條件

為了充分考慮外界環(huán)境對冷卻塔內(nèi)部氣水流動狀態(tài)的影響,本文計算過程中采取的計算域為高600 m、直徑800 m的大圓柱空間區(qū)域,冷卻塔置于圓柱空間區(qū)域的底部中心。計算域及對應(yīng)的邊界條件類型如圖2所示。

圖 2 計算域及邊界條件類型Fig.2 Calculation domain and boundary conditions

在無風(fēng)工況下,整個計算域的出口條件為壓力出口,入口邊界條件為壓力入口;在側(cè)風(fēng)環(huán)境下,計算域出口為壓力出口,入口為速度入口。外部空間區(qū)域與冷卻塔內(nèi)部區(qū)域以交界面邊界類型連接。假設(shè)聲屏障為實體墻,其邊界條件為壁面。為了模擬實際狀態(tài)下的側(cè)風(fēng)環(huán)境,在側(cè)風(fēng)工況下進口風(fēng)速大小與高度呈風(fēng)速廓線關(guān)系,如式(4)所示[22]。

(4)

式中:vref為參考位置的環(huán)境風(fēng)速,m/s;vz為高度z處的環(huán)境風(fēng)速,m/s;z表示位置高度,m。

2.3 計算方法

以Fluent軟件中集成的氣體流動模型計算濕空氣運動過程,編寫UDF程序計算循環(huán)水控制方程。以用戶自定義標(biāo)量UDS計算循環(huán)水控制方程中的參數(shù)。氣水運動的相互作用,包括傳熱、傳質(zhì)和動量交換等過程,以氣水控制方程的源項實現(xiàn)[23-25]。

以殘差和出塔水溫、出塔氣溫等關(guān)鍵參數(shù)的變化為依據(jù)判斷控制方程計算的收斂性。當(dāng)能量方程殘差小于10-6,其余殘差小于10-4,且塔內(nèi)關(guān)鍵參數(shù)不發(fā)生變化時,判斷計算過程收斂。

2.4 模型驗證

首先驗證網(wǎng)格的獨立性,其次驗證模型計算的準(zhǔn)確性?；贗CEM CFD進行網(wǎng)格劃分。分別創(chuàng)建了網(wǎng)格數(shù)量為115萬、140萬、188萬和196萬的4套網(wǎng)格系統(tǒng),以設(shè)計工況為計算工況進行計算,用以驗證網(wǎng)格的獨立性。計算結(jié)果如表3所示。綜合考慮計算精度和計算資源,本文選擇188萬的網(wǎng)格系統(tǒng)進行后續(xù)計算研究。

表 3 網(wǎng)格獨立性驗證

出塔水溫的設(shè)計值是基于能量守恒和質(zhì)量守恒原理,并結(jié)合修正系數(shù)綜合計算得出,在工程實踐中被證明和初建冷卻塔的實測值較為一致,常作為模型驗證的數(shù)據(jù)[3,26]。另外,本文研究對象已運行多年,填料面臨老化等問題,出塔水溫的實測值與最初的數(shù)據(jù)出現(xiàn)較大偏差,已不適合用于模型的驗證,因此本文將計算值和設(shè)計值進行比較,驗證模型的準(zhǔn)確性。在網(wǎng)格獨立性基礎(chǔ)上,對不同運行工況進行計算,比較計算值和設(shè)計值以驗證模型準(zhǔn)確性。如表4所示,本文模型最大誤差為3.53%,符合計算精度要求。

表 4 模型準(zhǔn)確性驗證

3 結(jié)果與討論

側(cè)風(fēng)是影響冷卻塔運行的重要變量。本文重點研究側(cè)風(fēng)工況下聲屏障對填料和配水協(xié)同優(yōu)化的濕式冷卻塔熱力性能的影響。首先分析聲屏障和兩區(qū)協(xié)同增效方案對塔內(nèi)溫度場、速度場和水蒸氣組分分布的影響,在此基礎(chǔ)上定量分析聲屏障和兩區(qū)協(xié)同增效方案對熱力和通風(fēng)性能的影響。

3.1 塔內(nèi)物理場演變規(guī)律

本文所涉及冷卻塔的累年平均風(fēng)速為2.5 m/s,在此側(cè)風(fēng)工況下研究冷卻塔內(nèi)各物理場的演變規(guī)律才具有實際意義,包括溫度場、空氣速度場和水蒸氣質(zhì)量分?jǐn)?shù)分布。為了充分對比,本文給出3種結(jié)構(gòu)下冷卻塔的計算結(jié)果,即未經(jīng)任何改造的原始塔、增加聲屏障的原始塔以及增加聲屏障及兩區(qū)協(xié)同增效改造的冷卻塔。

圖3(a)為塔內(nèi)溫度場云圖。由于側(cè)風(fēng)會誘導(dǎo)進風(fēng)口上沿氣流紊亂,故3種結(jié)構(gòu)下在進風(fēng)口的迎風(fēng)側(cè)都會出現(xiàn)小范圍的高溫區(qū)。增加聲屏障會在一定程度上縮減迎風(fēng)側(cè)進風(fēng)口上沿的高溫區(qū)范圍,此時聲屏障起到擋風(fēng)墻的作用。另外,由于外界冷風(fēng)不易進入冷卻塔中心,導(dǎo)致原始塔塔心出現(xiàn)明顯的高溫區(qū)。對比發(fā)現(xiàn),增加聲屏障會減小塔心高溫區(qū)的范圍。在增加聲屏障的基礎(chǔ)上,實施填料和配水的協(xié)同優(yōu)化措施會進一步減小迎風(fēng)側(cè)進風(fēng)口上沿和塔心區(qū)域的高溫區(qū)范圍。從溫度云圖的對比可知,在本文研究工況下,聲屏障和兩區(qū)協(xié)同增效方案都對塔內(nèi)溫度場的均勻分布產(chǎn)生有利影響。

圖3(b)為冷卻塔內(nèi)空氣流速分布云圖。在2.5 m/s側(cè)風(fēng)風(fēng)速下冷卻塔內(nèi)部空氣流速呈非對稱性分布特征,在進風(fēng)口迎風(fēng)側(cè)附近以及背風(fēng)側(cè)出風(fēng)口附近都出現(xiàn)小范圍的高流速區(qū)。另外,空氣繞流導(dǎo)致的旋渦現(xiàn)象致使進風(fēng)口上沿出現(xiàn)低風(fēng)速區(qū),但迎風(fēng)側(cè)的低風(fēng)速區(qū)范圍明顯高于背風(fēng)側(cè)。增加聲屏障和實施兩區(qū)協(xié)同增效方案會使迎風(fēng)側(cè)進風(fēng)口附近的高風(fēng)速區(qū)在一定程度上減小,但對塔內(nèi)整體的空氣流速分布影響并不明顯。

圖3(c)為水蒸氣質(zhì)量分?jǐn)?shù)分布云圖。與塔內(nèi)溫度場分布類似,由于空氣的繞流作用在進風(fēng)口上沿的迎風(fēng)側(cè)附近空氣流速變慢,出現(xiàn)水蒸氣質(zhì)量分?jǐn)?shù)較高的區(qū)域。另外,塔心區(qū)域空氣流速較小,氣水傳質(zhì)充分,出現(xiàn)水蒸氣質(zhì)量分?jǐn)?shù)較高的區(qū)域。對比可知,增加聲屏障后塔心區(qū)域的水蒸氣質(zhì)量分?jǐn)?shù)明顯減小,而實施兩區(qū)協(xié)同增效方案后塔心區(qū)域的水蒸氣質(zhì)量分?jǐn)?shù)進一步減小。

(a) 冷卻塔內(nèi)溫度場云圖

(b) 冷卻塔內(nèi)速度場云圖

(c) 冷卻塔內(nèi)水蒸氣質(zhì)量分?jǐn)?shù)云圖

3.2 熱力性能變化規(guī)律

循環(huán)水溫降和冷卻數(shù)是評價冷卻塔熱力性能最常用的2個參數(shù)[27-28],其中循環(huán)水溫降是冷卻塔熱力性能最直接的體現(xiàn)。水溫降隨環(huán)境風(fēng)速的變化如圖4所示?？梢钥闯?水溫降的變化趨勢在6 m/s風(fēng)速處存在臨界點。對比可知,在僅增加聲屏障這一改造下,冷卻塔的水溫將會略有減低。這可能是由于聲屏障在一定程度上對冷卻塔的進風(fēng)產(chǎn)生了削弱作用。但是,增加兩區(qū)協(xié)同增效方案后冷卻塔的循環(huán)水溫降又有明顯的提高。這說明兩區(qū)協(xié)同增效方案可以有效抵消聲屏障的不利影響。定量計算可知,聲屏障結(jié)合兩區(qū)協(xié)同增效方案后水溫降平均增加了0.06 ℃,在降低噪聲的同時保證了冷卻塔的正常運行。由于冷卻塔的循環(huán)水量巨大,較小的溫降即可引起較大的全塔散熱量變化。以本文冷卻塔39 186 m3/h循環(huán)水量計,0.06 ℃的溫降增量可引起每小時約9 875 MJ散熱量的變化。

圖 4 循環(huán)水溫降隨側(cè)風(fēng)風(fēng)速的變化Fig.4 The variation of circulating water temperature drop with crosswind

圖5為冷卻數(shù)隨側(cè)風(fēng)風(fēng)速的變化。可以看出:與循環(huán)水溫降的變化規(guī)律類似,冷卻數(shù)隨側(cè)風(fēng)風(fēng)速的增加也是先減小后增大。在無風(fēng)工況下增加聲屏障前后冷卻數(shù)相差不大,在側(cè)風(fēng)工況下聲屏障使冷卻數(shù)減小。但增加兩區(qū)協(xié)同增效方案后,冷卻數(shù)明顯增加。這說明,增加兩區(qū)協(xié)同增效方案可抵消聲屏障對冷卻數(shù)的不利影響。定量計算可知,聲屏障結(jié)合兩區(qū)協(xié)同增效方案可使冷卻數(shù)平均增加3.8%。

圖 5 冷卻數(shù)隨側(cè)風(fēng)風(fēng)速的變化Fig.5 The variation of Merkel number with crosswind

3.3 全塔通風(fēng)性能

冷卻塔通風(fēng)量隨環(huán)境風(fēng)速的變化趨勢如圖6所示?？梢钥闯?隨著環(huán)境風(fēng)速的不斷變大,通風(fēng)量變化趨勢也存在臨界點,同樣出現(xiàn)在6 m/s風(fēng)速處。通過對比發(fā)現(xiàn),在各個側(cè)風(fēng)風(fēng)速環(huán)境下原始塔的全塔通風(fēng)量都是最大的。這是由于采用聲屏障會在一定程度上影響冷卻塔的正常進風(fēng);另一方面在冷卻塔內(nèi)增加小片距填料也會形成更大的通風(fēng)阻力,最終導(dǎo)致全塔通風(fēng)量的減小。雖然在經(jīng)過改造后,冷卻塔的全塔通風(fēng)量有所減小,但是一方面填料的非均勻布置在塔內(nèi)形成了更為合理的空氣動力場分布,另一方面增加了整個填料區(qū)的熱力性能。

圖 6 全塔通風(fēng)量隨側(cè)風(fēng)風(fēng)速的變化Fig.6 The variation of ventilation with crosswind

4 結(jié) 論

1) 在本文研究工況條件下聲屏障和兩區(qū)協(xié)同增效方案可均勻塔內(nèi)溫度場、空氣流速分布和水蒸氣質(zhì)量分?jǐn)?shù)分布。

2) 兩區(qū)協(xié)同增效方案可在一定程度上削弱聲屏障的不利影響,使循環(huán)水溫降和冷卻數(shù)比原始塔分別平均增加了0.06 ℃和3.8%。

3) 雖然通過改造后全塔通風(fēng)量有所減小,但空氣動力場分布更為合理,最終使冷卻塔熱力性能有所增加。