陳金華，馬 晴，夏 磊，彭運林，郭建平

（1.重慶大學(xué) 國家級低碳綠色建筑國際聯(lián)合研究中心，三峽庫區(qū)生態(tài)環(huán)境教育部重點實驗室，重慶 400045；2.重慶大學(xué) 城市建設(shè)與環(huán)境工程學(xué)院，重慶 400045；3.重慶博建建筑設(shè)計有限公司，重慶 400020；4.云南省水利水電勘測設(shè)計研究院，云南 昆明 650021）

正在興建的干河泵站是目前亞洲最大的地下抽水泵站，該泵站主體廠房位于地下，如果通風(fēng)空調(diào)不良，長期在此環(huán)境下工作的人會感到疲倦、頭暈等，而且還會對機(jī)電設(shè)備造成危害.所以合理設(shè)計地下廠房的通風(fēng)空調(diào)系統(tǒng)，對維護(hù)地下廠房合理的溫、濕度環(huán)境，保障機(jī)電設(shè)備正常運轉(zhuǎn)和人員身體健康具有重要意義[1－2].對大型地下廠房通風(fēng)空調(diào)系統(tǒng)的研究，目前主要有3種方法：根據(jù)相似理論進(jìn)行模型試驗[3]，CFD計算流體模擬分析[4－8]和網(wǎng)絡(luò)模型數(shù)值模擬[9].本次模型試驗是在數(shù)值模擬計算基礎(chǔ)上進(jìn)行的，試驗?zāi)Ｐ偷脑O(shè)計采用了模擬計算的優(yōu)化方案.模型試驗的目的是通過實驗得出在夏季和過渡季節(jié)兩個季節(jié)，水泵不同運行模式及不同送風(fēng)溫度下，泵站各層工作區(qū)的溫度分布狀況，以此來驗證主廠房優(yōu)化送風(fēng)方案的合理性和可靠性.

1 地下廠房通風(fēng)系統(tǒng)方案

泵站地下廠房通風(fēng)系統(tǒng)主要由主廠房通風(fēng)系統(tǒng)、檢修閥室通風(fēng)系統(tǒng)、電纜豎井通風(fēng)系統(tǒng)等組成.

泵站主交通洞設(shè)送風(fēng)機(jī)室，將通過一層灌漿廊道自然冷卻或加熱的室外空氣沿主交通洞頂拱送風(fēng)管送至主廠房上部，在主廠房吊頂上從端部引2支送風(fēng)管利用噴口頂送風(fēng)方式將空氣送至電機(jī)層的各機(jī)組段.同時利用主廠房四周的防濕隔墻做通風(fēng)道，通過設(shè)置軸流風(fēng)機(jī)將電動機(jī)層的新鮮空氣通過夾墻風(fēng)管側(cè)送至廠房的中間層、水泵層和閥層.

主廠房的排風(fēng)設(shè)2個系統(tǒng)，一個系統(tǒng)通過工作廊道、母線廊道排至工作豎井，并通過設(shè)置在工作豎井頂部的風(fēng)機(jī)排至室外；另一系統(tǒng)由球閥室頂部排至二層灌漿廊道，通過設(shè)置在端部的風(fēng)機(jī)排至室外.整個地下廠房的通風(fēng)示意圖見圖1.

2 氣流組織分析

主廠房共4層，每層之間通過夾墻風(fēng)管和樓梯間等相連，各層之間氣流相互影響、貫通，構(gòu)成一個復(fù)雜的洞室通風(fēng)群，不同的氣流組織形式對電站的設(shè)備運行、工藝要求都會產(chǎn)生重要的影響，如何確定正確、合理的氣流組織是整個電站設(shè)計中的重點.

本試驗利用CFD軟件進(jìn)行數(shù)值計算，確定合適的送風(fēng)方式，為熱態(tài)模型試驗奠定基礎(chǔ).

圖1 整個廠房的通風(fēng)示意圖Fig.1 The entire plant ventilation schemes

2.1 氣流組織參數(shù)確定

1）送風(fēng)量.確定送風(fēng)量的基本原則為：送至廠房內(nèi)總空氣量，能夠排出廠內(nèi)的余熱和余濕，使廠內(nèi)的熱濕環(huán)境達(dá)到機(jī)組運行、工藝生產(chǎn)的要求.根據(jù)電站提供的采暖與空氣調(diào)節(jié)計算書可知，主廠房內(nèi)排除余熱所需總送風(fēng)量為53 809m3／h，排除余濕所需總送風(fēng)量為19 425m3／h，故確定主廠房送風(fēng)機(jī)送風(fēng)量54 000m3／h.

2）風(fēng)口形式.地下廠房為高大空間，需采用大風(fēng)量、長射流型的風(fēng)口形式，參考同類廠房采用圓形噴口進(jìn)行頂送風(fēng).

3）風(fēng)口數(shù)量、風(fēng)口尺寸及風(fēng)口風(fēng)速.本次計算中，采用假定參數(shù)法確定風(fēng)口的相關(guān)參數(shù)，即根據(jù)泵站的初設(shè)情況考慮風(fēng)口數(shù)量，假定送風(fēng)口風(fēng)速確定送風(fēng)口尺寸，并對假定的風(fēng)口布置形式進(jìn)行氣流組織模擬計算，對其通風(fēng)效果進(jìn)行檢驗.

根據(jù)泵站初設(shè)，風(fēng)口布置形式分為2種模式.模式1：采用8個風(fēng)口（只在電機(jī)層，不包括旁邊的安裝間），8個風(fēng)口均勻布置，每個電機(jī)上平均布置2個.模式2：采用18個風(fēng)口（其中安裝間上布置4個），每個電機(jī)上布置2個，相鄰電機(jī)上空布置2個.

當(dāng)采用模式1時，假定風(fēng)口流速分別取為8，10和12m／s 3種，對應(yīng)的噴口直徑分別為550，500和450mm.當(dāng)采用模式2時，假定風(fēng)口流速取為9和10m／s，對應(yīng)的噴口直徑為350mm，詳細(xì)的工況設(shè)置見表1，其中括號內(nèi)為安裝間噴口速度.

表1 噴口參數(shù)設(shè)置Tab.1 Nozzle parameter settings

2.2 計算過程及結(jié)果

計算過程中，采用三維模型，將動量方程與能量方程進(jìn)行分離迭代求解，送風(fēng)噴口設(shè)置為速度入口，排風(fēng)口設(shè)置為壓力出口，各層的夾墻風(fēng)機(jī)設(shè)置為壓力提升邊界.計算根據(jù)表1進(jìn)行，分為多種工況，每種工況均對整個廠房內(nèi)的速度場進(jìn)行分析，并按《水力發(fā)電廠廠房采暖通風(fēng)與空氣調(diào)節(jié)設(shè)計規(guī)程》的規(guī)定對計算結(jié)果進(jìn)行驗證，選出最優(yōu)氣流組織模式.

在工況1～工況3中，豎直切面速度云圖見圖2～圖4.在工況4～工況6中，工作區(qū)內(nèi)的速度云圖見圖5～圖7.

圖2 工況1豎直切面速度云圖Fig.2 The velocity nephogram of vertical section under condition one

圖3 工況2豎直切面速度云圖Fig.3 The velocity nephogram of vertical section under condition two

圖4 工況3豎直切面速度云圖Fig.4 The velocity nephogram of vertical section under condition three

圖5 工況4工作區(qū)速度云圖Fig.5 The velocity nephogram of work area under condition four

圖6 工況5工作區(qū)速度云圖Fig.6 The velocity nephogram of work area under condition five

圖7 工況6工作區(qū)速度云圖Fig.7 The velocity nephogram of work area under condition six

根據(jù)計算結(jié)果可知，不同風(fēng)口數(shù)量和風(fēng)口速度對電機(jī)層氣流組織影響較大，其余各層，由于采用側(cè)送風(fēng)，受噴口參數(shù)影響較??；所有計算工況中，當(dāng)噴口速度為12m／s時，電機(jī)層工作區(qū)平均流速為1.53m／s，不滿足規(guī)范要求，工況4～工況6的工作區(qū)平均流速雖然小于規(guī)范要求的0.8m／s，但從圖5～圖7可以明顯看出，這3種工況在局部區(qū)域有流速超標(biāo)現(xiàn)象.

通過在6種風(fēng)口數(shù)量和尺寸組合工況下的模擬計算，最終確定最佳的風(fēng)口布置方案為：電機(jī)層機(jī)組段布置直徑為500mm的頂送風(fēng)口8個，送風(fēng)速度應(yīng)控制在8～10m／s.

3 模型建立

3.1 各種相似比例尺的確定[10]

本項目以相似理論為依據(jù)，指導(dǎo)模型設(shè)計以及整個模型試驗研究工作.模型需要確定的比例尺主要有幾何比例尺、溫度比例尺、速度比例尺、熱量比例尺和風(fēng)量比例尺.

由于廠房內(nèi)有熱源散發(fā)熱量，廠房壁面吸熱量小于熱源散發(fā)的熱量，因而存在需通風(fēng)帶走的余熱.送風(fēng)溫度低于廠房內(nèi)的氣溫，送風(fēng)氣流屬于非等溫射流.射流所受的重力大于浮力，重力和浮力之差稱為有效重力，射流向下彎曲的程度取決于有效重力的大小，所以應(yīng)采用阿基米德模型律.

模型的幾何比例尺Cl是模型設(shè)計必須首先確定的參數(shù).為方便模型制作和保證試驗結(jié)果的準(zhǔn)確，同時考慮試驗場地的實際情況，本實驗選定的幾何比例尺為1∶10.確定了模型試驗的幾何比例尺Cl，進(jìn)而就要確定模型送風(fēng)溫度、送風(fēng)量、送風(fēng)速度，進(jìn)行熱量模擬.根據(jù)阿基米德模型律，有Arm＝Arp.由此可得：

其中重力加速度可以認(rèn)為是常數(shù)，即Cg＝1.模型與原型氣流的送風(fēng)絕對溫度之比CT0≈1.為了對比和方便計算，一般取溫差比例尺CΔT0＝1.模型地重慶和原型地云南大氣壓比例不等于1.當(dāng)模型雷諾數(shù)處于自模區(qū)時，可計算出模型的各種比例尺，結(jié)果見表2.

3.2 模型整體設(shè)計思路

在確定了各種相似比例尺后，按照實際工程圖紙，確定模型所需的幾大部分，具體分以下4部分.

模型主體結(jié)構(gòu)：模型骨架采用角鋼搭建，外圍壁面采用新型復(fù)合保溫板，層與層之間鋪設(shè)木工板，設(shè)備采用鍍鋅鋼板制作.發(fā)熱系統(tǒng)：電纜采用額定發(fā)熱量為15W／m的燈帶，其他設(shè)備、照明散熱都使用額定功率為15W／m的白熾燈.送排風(fēng)系統(tǒng)：利用空調(diào)實現(xiàn)不同溫度的送風(fēng).檢測系統(tǒng)：按照測試目標(biāo)要求，共布置了125個溫度測點.

表2 干河泵站通風(fēng)模型相似比例尺Tab.2 The ventilation model similar scale of Ganhe pumping station

4 模型試驗

本文模型試驗采用了模擬計算優(yōu)化的送風(fēng)方案，在電機(jī)層機(jī)組段頂部設(shè)置8個風(fēng)口.模型電機(jī)層送風(fēng)口速度3.62m／s，對應(yīng)的原型送風(fēng)口速度為10.3m／s.

在總風(fēng)量為54 000m3／h時，主廠房中間層、水泵層、球閥層及檢修閥層的計算溫度為限值溫度，從安全角度來考慮，以上各層的送風(fēng)量在前面計算的基礎(chǔ)上乘以富裕系數(shù)1.2.即總風(fēng)量調(diào)整為54 000 m3／h×1.2＝64 800m3／h，由于主廠房閥層相對濕度在計算時達(dá)到92.3%，因此加大送風(fēng)量，將總風(fēng)量取為66 000m3／h.

根據(jù)泵站提供夏季設(shè)計送風(fēng)溫度為23℃來確定模型夏季送風(fēng)溫度，應(yīng)電站要求，本試驗增加一個比實際送風(fēng)溫度略高的工況用以對比.過渡季節(jié)實際送風(fēng)溫度為16℃.當(dāng)?shù)囟就L(fēng)室外干球溫度為4℃，直接送入室內(nèi)，完全可以保證冬季室內(nèi)空氣設(shè)計最低溫度10℃.所以本項目試驗的重點在夏季和過渡季節(jié).因此，在以上3種溫度情況下，分別開啟1臺，2臺以及3臺機(jī)組安排試驗.

具體的試驗工況安排見表3，每種工況進(jìn)行2～3次重復(fù)試驗.

表3 試驗工況方案表Tab.3 Test conditions

5 模型試驗結(jié)果與分析

5.1 試驗數(shù)據(jù)無因次化

模型試驗要確定溫度場，溫度屬于待定參數(shù).為了由模型試驗結(jié)果預(yù)測原型，必需將待定參數(shù)無因次化.

各測點無因次溫度為：

式中：to為送風(fēng)溫度；te為排風(fēng)溫度；ti為測點測量溫度.

5.2 模型試驗結(jié)果轉(zhuǎn)換為原型值

原型電機(jī)層熱源總發(fā)熱量為24.6kW，外壁面無法做到絕對保溫，故冬夏兩季皆有向外傳熱可能性，取10%被“壁面吸收”，則由送風(fēng)帶走的發(fā)電機(jī)層的送風(fēng)余熱量為：

原型送排風(fēng)溫差Δt為：

令為模型射流空間平均無因次溫差，to為送風(fēng)溫度，則原型射流空間任意高度平均溫度為：

5.3 電機(jī)層溫度場模型試驗數(shù)據(jù)處理

電機(jī)層夏季工況溫度場模型試驗數(shù)據(jù)處理見圖8～圖10.

圖8 夏季工況1臺機(jī)組不同送風(fēng)溫度典型測點無因次溫度對比圖Fig.8 The measuring point dimensionless temperature contrast for a typical units of different air supply temperature

圖9 夏季工況2臺機(jī)組不同送風(fēng)溫度典型測點無因次溫度對比圖Fig.9 The measuring point dimensionless temperature contrast for two typical units of different air supply temperature

各試驗工況下，在廠房中部2臺機(jī)組（即2＃，3＃機(jī)組）測點的無因次溫度低于兩端的2臺機(jī)組（即1＃，4＃機(jī)組）測點的無因次溫度，這主要是由于兩端發(fā)熱量大于中部發(fā)熱量所致.各測試方案工作區(qū)無因次溫差分布的特征為：最高溫度出現(xiàn)在1＃和4＃機(jī)組處，且排風(fēng)溫度高于工作區(qū)平均溫度.各測點無因次溫差標(biāo)準(zhǔn)偏差不超過0.2.電機(jī)層工作區(qū)的溫度分布均勻性較好.

圖10 夏季工況3臺機(jī)組不同送風(fēng)溫度典型測點無因次溫度對比圖Fig.10 The measuring point dimensionless temperature contrast for three typical units of different air supply temperature

5.4 模型試驗結(jié)果

根據(jù)前述電機(jī)層夏季工況試驗數(shù)據(jù)處理的方法，我們可以得到夏季和過渡季節(jié)共9種試驗工況下各層的實驗結(jié)果，見圖11～圖12.

圖11 夏季各工況下各層平均溫度分布Fig.11 The average temperature distribution of each layer in summer

圖12 過渡季節(jié)各工況下各層平均溫度分布Fig.12 The average temperature distribution of each layer in transition season

從圖11～圖12可以看出，各工況下，廠房內(nèi)球閥2層溫度最高.夏季當(dāng)廠房送風(fēng)溫度為24.5℃時，球閥2層溫度達(dá)29.46℃，廠房平均溫度27.72℃，仍小于有關(guān)規(guī)范和設(shè)計要求的溫度（30℃），則方案可行.過渡季節(jié)時，廠房內(nèi)的最高溫度為25.36℃，遠(yuǎn)低于30℃，可以采用變風(fēng)量調(diào)節(jié)，減小總的送風(fēng)量，降低輸送能耗，實現(xiàn)節(jié)能.

6 結(jié) 論

干河泵站是目前亞洲修建的最大的地下抽水泵站.同地下水電站廠房一樣，也需要依靠一定的通風(fēng)空調(diào)方式來調(diào)節(jié)廠房內(nèi)的溫濕度，使室內(nèi)空氣參數(shù)達(dá)到規(guī)范的要求，保證運行檢修人員的身心健康及保障設(shè)備的正常運轉(zhuǎn).本文通過氣流組織模擬試驗以及通風(fēng)系統(tǒng)模型試驗對干河泵站地下廠房通風(fēng)空調(diào)系統(tǒng)進(jìn)行研究，得出結(jié)論如下：

1）不同頂送風(fēng)參數(shù)對電機(jī)層氣流組織影響較大，但對其余各層影響較小，對于頂送風(fēng)，在電機(jī)層滿足空氣射流原理、總送風(fēng)量54 000m3／h不變的情況下，通過在6種風(fēng)口數(shù)量和尺寸組合工況下的模擬計算，最終確定最佳的風(fēng)口布置方案為：電機(jī)層機(jī)組段布置直徑為500mm的頂送風(fēng)口8個，送風(fēng)速度控制在8～10m／s.

2）在夏季最高送風(fēng)溫度24.5℃時，廠房內(nèi)最高溫度為29.46℃，平均溫度為27.72℃，滿足相關(guān)設(shè)計規(guī)范對于溫度的要求，泵站通風(fēng)系統(tǒng)方案可行.電機(jī)層各測點無因次溫差的標(biāo)準(zhǔn)偏差最大不超過0.2，這充分反映了電機(jī)層工作區(qū)氣流分布的均勻性，驗證了數(shù)值模擬優(yōu)化的氣流組織方案.

3）在部分機(jī)組開啟及過渡季節(jié)運行時，室內(nèi)溫度較低，主廠房負(fù)荷減少，可減少送風(fēng)量使主廠房的熱濕環(huán)境達(dá)到規(guī)程要求，同時通風(fēng)系統(tǒng)采用變風(fēng)量運行可進(jìn)一步實現(xiàn)節(jié)能.

[1] 謝龍祥.水電站地下廠房通風(fēng)空調(diào)設(shè)計探討[J].浙江水利科技，1999，3：45－46.XIE Long－xiang.Ventilation and air conditioning design of hydropower station underground factory building[J].Zhejiang Hydrotechnics，1999，3 ：45－46.（In Chinese）

[2] MYREFELT S.The reliability and availability of heating ventilation and air conditioning systems[J].Energy and Buildings，2004，6：17－20.

[3] 李曉豐.瑯琊山水電站地下廠房通風(fēng)模型裝置設(shè)計與實驗研究[D].重慶：重慶大學(xué)城市建設(shè)與環(huán)境工程學(xué)院，2003：8－10.LI Xiao－feng.Model installation design and study on ventilation experiment of Lang Ya Shan hydroelectric power station’s underground houses[D].Chongqing：School of Urban Construction ＆ Environmental Engineering，Chongqing University，2003：8－10.（In Chinese）

[4] 龍?zhí)煊?，采增基，董宏明，?水電站地下廠房通風(fēng)氣流組織數(shù)值模擬[J].建筑熱能通風(fēng)空調(diào)，2000，4：15－18.LONG Tian－yu，CAI Zeng－ji，DONG Hong－ming，etal.The numerical simulation for the ventilation of a underground main workshop in hydropower station[J].Building Energy ＆ Environment，2000，4：15－18.（In Chinese）

[5] NIELSEN P V.The selection of turbulence models for prediction of room airflow[J].ASHRAE Transactions，1998，104：1119－1126.

[6] NIELSEN P V.Description of supply openings in numerical models for room air distribution[J].ASHRAE Transactions，1992，98（1）：963－971.

[7] CHEN Qing－yan，XU Wei－ran.A zero equation turbulence model for indoor air flow simulation[J].Energy and Building，1998，28：137－144.

[8] 周良，張國強，彭建國，等.建筑節(jié)能及暖通空調(diào)仿真軟件的現(xiàn)狀和發(fā)展[J].湖南大學(xué)學(xué)報：自然科學(xué)版，2000，27（6）：103－107.ZHOU Liang，ZHANG Guo－qiang，PENG Jiang－guo，etal.Literature review of building energy simulation software[J].Journal of Hunan University：Natural Sciences，2000，27（6）：103－107.（In Chinese）

[9] 肖益民.水電站地下洞室群自然通風(fēng)網(wǎng)絡(luò)模擬及應(yīng)用研究[D].重慶：重慶大學(xué)城市建設(shè)與環(huán)境工程學(xué)院，2005：4－9.XIAO Yi－min.Studies on natural ventilation network simulation and application of the underground excavation cluster of hydropower station[D].Chongqing：School of Urban Construction ＆ Environmental Engineering，Chongqing University，2005：4－9.（In Chinese）

[10] 李安桂，唐懷林，張建飛，等.HHHT抽水蓄能電站地下主廠房母線層通風(fēng)模型試驗研究[J].制冷與空調(diào)，2009，23（5）：1－6.LI An－gui，TANG Huai－lin，ZHANG Jian－fei，etal.Ventilation model tests on HHHT pumped－storage powerhouse busbar layer[J].Refrigeration and Air Conditioning，2009，23（5）：1－6.（In Chinese）