施 晨，胡代清，衣傳寶，劉玉成，許 昌，李林敏

(1.河海大學(xué)能源與電氣學(xué)院，江蘇南京211100；2.安徽金寨抽水蓄能有限公司，安徽六安237300；3.國網(wǎng)新源控股有限公司，北京100761)

0 引 言

隨著電網(wǎng)結(jié)構(gòu)不斷復(fù)雜化、擾動(dòng)的不斷多樣化，抽水蓄能電站作為電網(wǎng)電能品質(zhì)的調(diào)節(jié)器，近年建設(shè)規(guī)模不斷擴(kuò)大[1]。而隨著大型地下工程建設(shè)技術(shù)的成熟，抽水蓄能電站的設(shè)計(jì)建設(shè)逐漸傾向于采用地下廠房方案。這種布局不僅具有人防作用，同時(shí)能最大程度地保護(hù)電站建設(shè)區(qū)域的生態(tài)環(huán)境[2]。但由于電站深埋地下，電站與外界通風(fēng)換熱途徑較為復(fù)雜，電站內(nèi)易出現(xiàn)溫濕度分布不均勻和超標(biāo)等熱濕環(huán)境問題，不僅影響機(jī)電設(shè)備的安全運(yùn)行，同時(shí)也影響工作人員的身體健康與工作效率。

國內(nèi)外學(xué)者針對水電站地下廠房溫濕度環(huán)境問題進(jìn)行了大量研究工作。在理論計(jì)算方面，劉琳[3]結(jié)合射流理論對白鶴灘水電站地下廠房整體進(jìn)風(fēng)網(wǎng)絡(luò)進(jìn)行風(fēng)速等參量的計(jì)算確定了空氣處理方案；Tuve[4]采用射流理論對室內(nèi)環(huán)境中的速度、溫度等重要參量的分布進(jìn)行預(yù)測。在模型試驗(yàn)上，何喆[5]對瑯琊山水電站地下廠房建立了相似模型，通過模型試驗(yàn)的數(shù)據(jù)分析得出溫度、速度分布規(guī)律，為流場優(yōu)化提供了參考；在數(shù)值計(jì)算上，丹麥科學(xué)家Nielsen等[6]于1974年首次將基于k-ε湍流模型的CFD計(jì)算方法對室內(nèi)空氣流動(dòng)進(jìn)行了三維數(shù)值模擬，而在水電站地下廠房的數(shù)值計(jì)算上；董蕾[7]對厄瓜多爾科卡科多辛克雷水電站地下廠房發(fā)電機(jī)層的氣流組織進(jìn)行了數(shù)值模擬及分析，并運(yùn)用不同的評價(jià)指標(biāo)確定了通風(fēng)系統(tǒng)的最佳送風(fēng)速度。當(dāng)前針對地下廠房溫濕度的數(shù)值研究多集中于地下廠房單獨(dú)某層的溫濕度環(huán)境分析，而針對主廠房全廠的數(shù)值研究較為缺乏。

抽水蓄能電站的主廠房是電站運(yùn)維的關(guān)鍵區(qū)域，該區(qū)域布置復(fù)雜，機(jī)電設(shè)備密集，散熱量巨大，運(yùn)維人員往來頻繁，研究其環(huán)境溫濕度的分布規(guī)律具有重要意義，過去的研究主要采用針對部分區(qū)域的溫濕度問題進(jìn)行測量和數(shù)值建模方面的研究，解決了部分問題，但是地下廠房的整體通風(fēng)方式影響整體溫濕度、速度場分布，且相互影響，相關(guān)耦合性非常強(qiáng)。本文采用現(xiàn)場實(shí)測與數(shù)值模擬相結(jié)合的方法，對華東某抽水蓄能電站地下主廠房全廠區(qū)域進(jìn)行數(shù)值建模，結(jié)合現(xiàn)場試驗(yàn)數(shù)據(jù)對不同季節(jié)工況下的全廠熱濕環(huán)境進(jìn)行數(shù)值模擬，分析不同工況下全廠各樓層的溫濕度分布規(guī)律，提出優(yōu)化運(yùn)行方案，為抽水蓄能電站廠房溫濕度環(huán)境控制與管理提供理論依據(jù)。

1 數(shù)學(xué)模型

1.1 通用控制方程

將流體動(dòng)力學(xué)物理守恒定律中的質(zhì)量守恒、動(dòng)量守恒與能量守恒方程寫成通用形式如下[8]：

(1)

式中，φ為通用變量，可代表速度、溫度等變量；Γ為廣義擴(kuò)散系數(shù)；S為廣義源項(xiàng)。式(1)從左往右的四項(xiàng)依次是非穩(wěn)態(tài)項(xiàng)、對流項(xiàng)、擴(kuò)散項(xiàng)和源項(xiàng)。

在地下廠房熱濕環(huán)境的計(jì)算中，相對濕度是一個(gè)重要的考察對象，本文使用組分輸運(yùn)方程來計(jì)算相對濕度的輸運(yùn)規(guī)律，令φ=cs，Γ=Dsρ，S=Ss，則式(1)可寫成：

(2)

式中，cs為組分s的體積濃度；Ds為組分s的擴(kuò)散系數(shù)；Ss為組分s的產(chǎn)生率。

1.2 湍流模型

抽水蓄能電站地下廠房風(fēng)口眾多，流動(dòng)復(fù)雜，流動(dòng)過程中具有不斷變化、隨機(jī)的湍流現(xiàn)象，本文采用標(biāo)準(zhǔn)k-ε湍流模型來計(jì)算[9，10]。其中湍流動(dòng)能k與湍動(dòng)能耗散率ε分別被定義為

(3)

(4)

式中，ρ是流體密度；i、j、k分別代表各坐標(biāo)分量；ui、uj代表平均相對速度分量；μ為有效粘性系數(shù)。k和ε的輸運(yùn)方程如下：

(5)

(6)

其中，Gk是由于平均速度梯度引起的湍動(dòng)能源項(xiàng)，由下式計(jì)算：

(7)

而Gb是由于浮力引起的湍動(dòng)能源項(xiàng)，對于不可壓縮流體，其值為0，而對于可壓縮流體，有：

(8)

式中，Prl是Prandtl數(shù)，在該模型中可取Prl=0.85；gi是重力加速度在i方向的分量；β為熱膨脹系數(shù)。YM代表可壓湍流中脈動(dòng)擴(kuò)張的貢獻(xiàn)，對于不可壓流體，YM=0。模型常數(shù)C1ε、C2ε、Cμ、σk、σε的取值分別為[11]C1ε=1.44，C2ε=1.92，Cμ=0.09，σk=1.0，σε=1.22。

2 物理模型

2.1 模型與網(wǎng)格劃分

華東某抽水蓄能電站主廠房共分4層，自上而下依次為發(fā)電機(jī)層、中間層、水輪機(jī)層以及蝸殼層。發(fā)電機(jī)層在靠近進(jìn)廠交通洞區(qū)域設(shè)有安裝間，面積略大，發(fā)電機(jī)層尺寸為140.3 m×20.00 m(長×寬)，其余3層尺寸均為120.2 m×20.00 m(長×寬)，主廠房總高37.8 m，4層凈高自上而下依次為17.7、5.25、5.95、6.65 m。

主廠房進(jìn)風(fēng)來源共有3條途徑：①交通洞→9號(hào)洞→3號(hào)洞→主廠房安裝間頂部空調(diào)；②通風(fēng)兼安全洞→端副廠房頂部空調(diào)；③交通洞→9號(hào)洞→水輪機(jī)層進(jìn)風(fēng)口。主廠房頂部共有46個(gè)進(jìn)風(fēng)口，4臺(tái)組合式空調(diào)機(jī)組，層間夾墻風(fēng)口共有44個(gè)。新風(fēng)從主廠房頂部進(jìn)入，通過層間風(fēng)口聯(lián)通，最終從中間層下游側(cè)排向母線洞，主廠房空氣循環(huán)見圖1。

圖1 主廠房空氣循環(huán)示意

對主廠房的通風(fēng)網(wǎng)絡(luò)以及三維立體結(jié)構(gòu)以1∶1的比例進(jìn)行建模，并對其進(jìn)行網(wǎng)格劃分，所建立的三維幾何模型如圖2所示。

圖2 主廠房幾何模型

表1 進(jìn)口邊界條件

為方便研究，根據(jù)實(shí)際情況對模型進(jìn)行簡化：①簡化主廠房發(fā)電機(jī)層頂棚弧形結(jié)構(gòu)以及頂部橋機(jī)等繁瑣物件，主廠房內(nèi)球閥、管道等非主要物件省略或以規(guī)則實(shí)體等提及替代。②由于空調(diào)送風(fēng)管道采用雙層保溫材料，將空調(diào)機(jī)組至送風(fēng)口管道省略，近似認(rèn)為空調(diào)出口送風(fēng)參數(shù)等于送風(fēng)口送風(fēng)參數(shù)。③省略主廠房壁面中層間通風(fēng)管道以及軸流風(fēng)機(jī)，僅保留風(fēng)口，以風(fēng)口對應(yīng)的現(xiàn)場測試數(shù)據(jù)作為計(jì)算的邊界條件。

2.2 邊界條件與數(shù)值細(xì)節(jié)

電站主廠房頂端途徑①與②進(jìn)風(fēng)口處各有2臺(tái)組合式空調(diào)機(jī)組，且兩側(cè)各運(yùn)行1臺(tái)，備用1臺(tái)，途徑③處設(shè)置風(fēng)機(jī)。本文計(jì)算春、夏兩種工況，分別以2017年3月20日與7月15日現(xiàn)場測試的溫濕度、壁溫?cái)?shù)據(jù)以及文獻(xiàn)[12]中的熱源散熱量作為參考數(shù)據(jù)，以確定邊界條件。各進(jìn)風(fēng)途徑風(fēng)速分別由各自風(fēng)機(jī)進(jìn)風(fēng)量換算得出，途徑③溫濕度為直接測量數(shù)據(jù)，空調(diào)處理后的空氣焓值由下式給出[13]：

(9)

式中，Q為空調(diào)制冷量；ma為空氣質(zhì)量流量，兩者均可由空調(diào)銘牌得出；h1為初態(tài)比焓；h2為終態(tài)比焓，最終溫濕度參數(shù)由h2查詢焓濕圖得到。兩種工況下各進(jìn)口邊界條件設(shè)置如表1所示。除主廠房頂部46個(gè)風(fēng)口向主廠房發(fā)電機(jī)層的射流送風(fēng)外，發(fā)電機(jī)層通過上下游側(cè)壁面處的夾墻風(fēng)管向下面3層送風(fēng)形成循環(huán)，本文簡化風(fēng)管建模，僅保留風(fēng)口，對應(yīng)的送風(fēng)風(fēng)口分別為速度進(jìn)口(volecity-inlet)與速度出口(volecity-outlet)，具體數(shù)值根據(jù)現(xiàn)場測試數(shù)據(jù)確定。最終出口為中間層下游側(cè)的回風(fēng)口，邊界條件為壓力出口(pressure-outlet)。

壓力-速度耦合采用Simple算法，湍流模型采用k-ε雙方程模型，由于涉及組分以及傳熱計(jì)算，開啟組分方程以及能量方程，同時(shí)采用DO輻射模型[14，15]，重力加速度g=9.8 m/s2。本文計(jì)算采用FLUENT求解器求解。

3 計(jì)算結(jié)果分析與優(yōu)化方案論證

3.1 模型可靠性驗(yàn)證

在現(xiàn)場測試時(shí)，針對主廠房4層共布置了20個(gè)測量點(diǎn)，每層測點(diǎn)5個(gè)，其中發(fā)電機(jī)層熱源較少，不布置測點(diǎn)。測量數(shù)據(jù)時(shí)采用手持溫濕度監(jiān)測設(shè)備，測點(diǎn)高度均為1.7 m。各層測點(diǎn)在水平方向上的分布如圖3所示。

圖3 各層測點(diǎn)示意

以春季工況為例，將2017年3月20日的各層實(shí)測值分別取平均值，與對應(yīng)位置處1.7 m高度處的計(jì)算值對比，如圖4所示。

圖4 可靠性驗(yàn)證

由圖4可知，各層溫度和濕度的實(shí)測值與計(jì)算值總體吻合。經(jīng)過計(jì)算，各層溫度計(jì)算值與實(shí)測值平均誤差為0.5 K，最大誤差1.24 K，最小誤差0.29 K。相對濕度計(jì)算值與實(shí)測值平均誤差0.67%，最大及最小誤差分別為1.50%和0.22%。

3.2 計(jì)算結(jié)果分析

圖5 春季與夏季工況的溫濕度計(jì)算結(jié)果

根據(jù)春季、夏季兩種工況的計(jì)算結(jié)果，分別截取y=3.65 m豎直截面處2種工況下的溫度與相對濕度云圖，結(jié)果如圖5所示。根據(jù)計(jì)算結(jié)果可以看出，在春夏2種工況下，發(fā)電機(jī)層的溫度在全廠中均為最低，其余3層溫度相對較高。且由于所選取截面頂部進(jìn)風(fēng)途徑①下游側(cè)空調(diào)運(yùn)行，另一端的進(jìn)風(fēng)途徑②下游側(cè)空調(diào)處于備用狀態(tài)，未開啟運(yùn)行，在發(fā)電機(jī)層左側(cè)區(qū)域看到自上而下的射流組織，且出現(xiàn)溫度較低的低溫區(qū)，溫度分布不均。在4層空間中發(fā)電機(jī)層電氣設(shè)備散熱量最大，但發(fā)電機(jī)層空間范圍最大且直接接受頂部低溫射流送風(fēng)，使該層溫度降低。而下3層空間較小且均利用發(fā)電機(jī)層氣流進(jìn)行通風(fēng)，通風(fēng)效果不佳。發(fā)電機(jī)層氣流存在不均勻性，這主要由于發(fā)電機(jī)層的安裝間范圍較大且不存在熱源，發(fā)電機(jī)層熱源集中在機(jī)組附近，熱源分布不均。相對濕度上，2個(gè)工況的發(fā)電機(jī)層頂部射流送風(fēng)口相對濕度較高，這是由于氣流經(jīng)過空調(diào)處理后溫度降低，空氣中水分迅速達(dá)到露點(diǎn)溫度并隨著溫度降低不斷下降，從而達(dá)到降低濕度的目的。含濕量得到降低的氣流離開空調(diào)之后幾乎以100%的相對濕度進(jìn)入溫度較高的發(fā)電機(jī)層，與發(fā)電機(jī)層內(nèi)的熱空氣混合后，溫度迅速升高，相對濕度迅速降低。

兩種工況所不同的是：溫度上，主廠房整體區(qū)域夏季高于春季，但均能低于302.2 K；濕度上，水輪機(jī)層夏季明顯高于春季，且部分情況下高于75%，對機(jī)電設(shè)備安全運(yùn)行形成一定威脅，也影響運(yùn)維人員的人體舒適度。這是因?yàn)?，水輪機(jī)層進(jìn)風(fēng)途徑③引進(jìn)主廠房的新風(fēng)相對濕度接近100%，而溫度與水輪機(jī)層環(huán)境溫度幾乎沒有差異，新風(fēng)含濕量較高，造成了水輪機(jī)層相對濕度較高。除此之外，因熱空氣會(huì)因自然對流而上升，水輪機(jī)層的濕空氣通過樓梯口進(jìn)入中間層，一定程度上造成了中間層部分區(qū)域相對濕度的升高。

3.3 優(yōu)化方案計(jì)算及論證

目前電站的運(yùn)行過程中，出現(xiàn)了相對濕度過高的問題，為提升機(jī)電設(shè)備檢修環(huán)境與運(yùn)維人員的人體舒適度，提出以下優(yōu)化運(yùn)行方案：假定在水輪機(jī)層進(jìn)風(fēng)途徑③處加設(shè)組合式空調(diào)機(jī)組，根據(jù)式(9)以及焓濕圖進(jìn)行計(jì)算，假定選型為申菱40HP的LD130型號(hào)組合式空調(diào)機(jī)組，優(yōu)化前后途徑③進(jìn)風(fēng)參數(shù)如表2所示。

表2 優(yōu)化前后途徑③進(jìn)風(fēng)參數(shù)

圖6 優(yōu)化前后溫濕度環(huán)境對比

其余送風(fēng)參數(shù)與夏季工況下的參數(shù)一致，截取水輪機(jī)層相對高度zr=1.7 m處的水平截面的溫濕度云圖進(jìn)行分析，其結(jié)果如圖6所示。由計(jì)算結(jié)果可得，優(yōu)化前，由于水輪機(jī)的發(fā)熱及其與周圍空氣的換熱導(dǎo)致水輪機(jī)周圍溫度較高，溫度分布不均勻，原空調(diào)系統(tǒng)對該層空氣的冷卻效果不明顯。優(yōu)化后，水輪機(jī)層溫度有所下降，且溫度分布更加均勻。且水輪機(jī)層送風(fēng)溫度的降低，使其與上下層之間溫度差降低，優(yōu)化后樓梯口的溫度均勻性優(yōu)于優(yōu)化之前。優(yōu)化后水輪機(jī)層環(huán)境溫度下降，有利于機(jī)電設(shè)備的安全運(yùn)行。另一方面，水輪機(jī)層相對濕度在優(yōu)化前存在超標(biāo)的情況，圖6c中大部分區(qū)域相對濕度達(dá)到了80%。優(yōu)化后，相對濕度得到了有效降低，如圖6d所示，大部分區(qū)域相對濕度在74%以內(nèi)。可見，優(yōu)化措施對于控制水輪機(jī)層的空氣濕度具有較好的效果。

4 結(jié) 論

本文建立了基于CFD的抽水蓄能電站地下主廠房全廠溫濕度、速度場的數(shù)值計(jì)算模型，通過試驗(yàn)數(shù)據(jù)驗(yàn)證了模型的可靠性，分析了全廠溫濕度和速度的分布特性。針對運(yùn)行過程中，水輪機(jī)層濕度過高的問題，提出在進(jìn)風(fēng)口加設(shè)組合式空調(diào)機(jī)組除濕的方案。通過數(shù)值計(jì)算，在夏季工況下，水輪機(jī)層大部分區(qū)域相對濕度從最高80%下降到74%以內(nèi)，而且通風(fēng)效果明顯提升，優(yōu)化效果較明顯。