城市地下道路分(合)流匝道通風阻力特性

袁浩庭， 陳 超, *， 李 瓊, 2， 王陸瑤， 李亞茹

(1. 北京工業大學建筑工程學院， 北京 100124； 2. 華北科技學院建筑工程學院， 北京 101601)

城市地下道路分(合)流匝道通風阻力特性

袁浩庭1， 陳 超1, *， 李 瓊1, 2， 王陸瑤1， 李亞茹1

(1. 北京工業大學建筑工程學院， 北京 100124； 2. 華北科技學院建筑工程學院， 北京 101601)

為分析多點進出結構的城市地下道路空氣流動特性，以長沙市營盤路湘江隧道為原型，通過現場實測、縮尺模型試驗以及CFD軟件數值模擬方法，對分(合)匝道通風阻力特性進行了研究，考察了雷諾數、風量比、分岔角度對分(合)流匝道阻力特性的影響規律。研究結果表明： 1)主隧道與匝道風量比、主隧道與匝道夾角是影響分(合)流匝道局部阻力系數的關鍵因素； 2)城市地下道路分(合)流匝道局部阻力變化特性，不宜簡單套用通風管道的三通構件的參數； 3)結合最小二乘法和Matlab軟件對計算結果進行多因素回歸分析，給出了基于長沙營盤路湘江隧道的分(合)流匝道主隧道和匝道的局部阻力系數關聯式。研究結果可為復雜結構城市地下道路通風系統阻力特性分析及通風工程優化設計提供方法參考。

多點進出城市地下道路； 分(合)流匝道； 局部阻力系數； 縮尺模型實驗； 數值模擬

0 引言

為緩解城市用地緊張、城市中心區主干道交通擁堵、縮短城區間的道路里程，多點進出城市地下道路已成為目前我國城市地下道路建設的主要隧道形式。例如，中國深圳新經濟特區前海正在新建的一條穿越CBD區域的長約5.257 km的城市地下道路，沿線設置9個進出地面匝道的出入口、10個進出地下匝道的出入口。復雜的隧道結構特征，直接影響城市地下道路通風工程設計，特別是分(合)流匝道通風阻力特性會對多點進出城市地下道路空氣流動特性產生直接影響。

關于隧道局部構件阻力特性的研究主要集中在公路隧道。李峰[1]對公路隧道彎曲風道進行了研究，得到了彎曲風道局部阻力損失系數隨彎曲半徑與隧道當量直徑的關系； 溫玉輝[2]采用數值模擬的方法對影響短道流態的4個因素(排風道與隧道軸線夾角、送風方向與隧道軸線夾角、送排風口間距和送排風流量比)進行單獨分析； 喻映華等[3]通過數值模擬的方法對隧道緊急停車帶的局部阻力系數進行了研究。同時部分學者對隧道三通局部阻力損失做了相關研究。石平等[4]結合《雁門關公路隧道運營通風與防災技術研究》，采用理論分析、現場測試及數值模擬等方法，對三通局部阻力損失、送排風短道的優化進行了研究，得到了一系列有益于通風設計的結論； 馬佳[5]運用數值模擬方法對影響公路隧道三通局部阻力系數的因素進行了研究，發現傾斜角度、主支流量比等因素對局部阻力系數有影響； 茍紅松等[6]利用理論計算對隧道施工通風分風三通位置的選擇做了研究，發現主風管與分風管直徑比是影響風管出口風量的關鍵因素。國內外關于三通局部阻力損失相關的研究大多集中在管道三通局部構件上。COSTA等[7]對夾角為90°的三通局部阻力系數進行了研究，分析了雷諾數和流量比對局部阻力系數的影響； Abdulwahhab等[8]對夾角為90°的三通的局部阻力損失進行了模擬和試驗研究，結果表明流量比對其阻力損失影響很大； 我國的很多學者也對管道三通局部阻力特性進行了研究[9-12]，分析了分流比、雷諾數、管間夾角和管徑比等因素的影響。

目前，國內關于多點進出城市地下道路通風系統基礎數據取值多是來自參考文獻[13]，其關于支管與主管夾角對管道空氣流動的影響，僅僅給出了夾角為30°、45°和90°時的3種情況； 現行規范[14]涉及到多點進出隧道主隧道與分(合)流匝道夾角對隧道空氣流動的影響，也只是給出了夾角為90°時的情況。考慮到實際隧道工程分(合)流匝道與常見通風三通管道的不同性，以及實際隧道主隧道與分(合)流匝道角度通常都小于30°的實際情況，本文以多點進出城市地下道路的分(合)流匝道通風阻力特性為重點研究對象，通過現場實測、縮尺模型試驗和數值模擬的方法，研究并分析雷諾數Re、主隧道與分(合)流匝道風量比q、主隧道與分(合)流匝道夾角θ等因素對分(合)流匝道局部阻力系數的影響規律，以期為城市地下道路通風系統設計提供參考。

1 材料與方法

1.1 縮尺模型試驗

1.1.1 試驗臺搭建

長沙市營盤路湘江隧道(以下簡稱長沙隧道)位于長沙市中心城區，為雙洞單向隧道，由南線和北線2條獨立且平行的隧道構成，主線為雙車道、總長2.7 km，每條隧道均有2個入口和2個出口，其中東側出(入)口對接湘江東岸商業區，西側出(入)口對接湘江西岸居住區，如圖1所示。 其中，A匝道長648 m、B匝道長425 m、C匝道長632 m、D匝道長752 m，均為單車道。匝道雖為單車道，但同時考慮了1條應急車道，因此主、匝隧道斷面尺寸幾乎一樣。主、匝道斷面結構尺寸見圖2。

圖1 線路示意圖

(a) 主隧道斷面構造圖

(b) 匝道斷面構造圖

根據流體力學理論以及流動相似理論，本文以長沙隧道為原型，基于雷諾相似準則和幾何相似準側搭建1∶8縮尺模型比例試驗臺，試驗臺尺寸見表1。

表1 模型試驗臺尺寸

考慮到實際隧道通風速度通常都在紊流區，雷諾數Re>5×104，滿足自模區的要求，因此主要考慮在幾何相似和邊界條件相似的前提下，使模型的空氣流動充分發展成完全紊流，進入阻力平方區，即自動模型區，確保與原型流動保持相似。基于管內自模區臨界雷諾數大于5×104，且可適當降低的設計原則[15]，模型試驗最小風速按0.69 m/s，對應臨界雷諾數Re>3.5×104。

按照幾何相似準則，試驗場地無法滿足模型長度要求。為此，本研究采用變態模型法對模型的長度進行補償修正，且將變態率控制在4以內。變態率

k=CL/CD=4[15]。

(1)

式中：k為變態率；CD為徑向比尺，CD=DP/Dm(DP為原型當量直徑，Dm為模型當量直徑)，本試驗臺為8；CL為長度方向的比尺，CL=LP/Lm(LP為原型長度，Lm為模型長度)。

根據式(1)，模型長度雖可縮短到2 700/32=84.37 m，但試驗場地仍然難以滿足要求。因此進一步采用等效模擬法，即加阻力格柵的方式，以局部阻力代替沿程阻力，對模型長度進行補償修正。試驗模型的最終長度為30 m。

圖3為1∶8模型試驗臺示意圖，試驗臺主要由風機、軟連接、整流段、格柵段、主隧道段、變徑段和匝道段等組成。其中變徑段長5 m，風管材料采用鍍鋅板。

圖3 1∶8縮尺隧道模型試驗臺

1.1.2 測點布置及數據采集

縮尺試驗臺各功能段構成及其風速、風壓等參數采樣點位置如圖4所示，風速、風壓測試儀器以及通風機等性能參數見表2，各采樣點位置按照上游大于2倍當量直徑、下游大于3～5倍當量直徑的原則確定。分流試驗時，主隧道風機正轉； 合流試驗時，主隧道風機反轉。當主隧道風速穩定后，采集各個采樣點的風速和風壓，風速每20 s采集1次，連續采集30 min； 風壓每5 min采集1次，采集30 min，共采集6組數據。

圖4 分(合)流段通風試驗測點布置圖

Fig. 4 Layout of monitoring points for ventilation of confluence and distributary sections

1.2 數值計算分析

1.2.1 物理模型及網格劃分

根據圖5建立相應的CFD計算物理模型，采用Gambit軟件建立相應的結構型網格模型，如圖6所示。通過網格獨立性檢驗，模型長度方向(z軸方向)、寬度方向(x軸方向)和高度方向(y軸方向)的網格尺度分別為0.1、0.062 5、0.062 5 m； 總網格數為55 250個。

圖5 計算物理模型

圖6 網格模型

1.2.2 模擬工況

根據長沙隧道不同工況下的實測結果，直隧道通常風速為2～7 m/s，支隧道風速為1～4 m/s。為考察雷諾數、風量比(包括上下游主隧道的風量比q1以及上游主隧道與匝道的風量比q2)以及主隧道與匝道的夾角等因素對分(合)流匝道局部阻力特性的影響規律，本研究分3種工況進行分析討論，具體工況見表3。

表3 計算工況

1.2.3 模型計算及求解

1)模型求解條件。應用Fluent15.0軟件對分(合)流匝道局部阻力系數進行分析。采用k-epsilon(2eqn)湍流模型，方程的離散采用有限體積法，擴散項采用中心差分格式離散。為提高計算精度，對流項采用二階迎風格式離散，速度與壓力的耦合計算采用SIMPLEC 算法。

2)邊界條件。進出口條件： 設置速度入口或速度出口邊界，采用湍流強度和水力直徑定義湍流，采用壓力出口或壓力入口邊界條件； 壁面條件： 固體壁面上采用無滑移條件，當量粗糙度為1 mm。

3)計算流體條件。空氣溫度為303 K、密度為1.225 kg/m3、動力黏度為1.789 4×10-5kg/(m·s)。

1.2.4 模型有效性驗證

依據流體力學理論，風道構件局部阻力系數

ξ=Δp/p。

(2)

式中：ξ為局部阻力系數； Δp為2個測點之間的全壓差；p為主隧道動壓。

圖7示出分(合)流匝道局部阻力系數與風速的關系圖，是模擬結果的有效性驗證，分岔匝道和合流匝道(θ=0°)局部阻力系數計算值與試驗結果具有較好的一致性。

2 結果及分析

2.1 雷諾數的影響(工況1)

圖8示出風量比q1=q2=0.5、主隧道與匝道夾角θ=0°時，主隧道風速變化對分流匝道、合流匝道局部阻力系數的影響有限，局部阻力系數變化很小。說明當流動處于紊流狀態即阻力平方區時，局部阻力系數不再隨Re變化。鞏啟濤等[9]發現，對于三通局部阻力系數，進入局部阻力系數的“阻力平方區”后，三通局部阻力系數基本趨于穩定。

(a) 分流匝道

(b) 合流匝道

Fig. 7 Relationships between local resistance coefficient of confluence and distributary ramps and air speed

2.2 風量比的影響(工況2)

圖9示出主隧道與匝道夾角為0°和30°時，風量比變化對分(合)流匝道局部阻力系數的影響規律。分流匝道情況下，主隧道和匝道的局部阻力系數都隨風量比呈先下降后增大的趨勢，其最小值分別出現在q1=0.6和q2=0.5左右，該規律不同于文獻[13]關于矩形風管分流三通中直通管道局部阻力系數的變化規律； 合流匝道情況下，主隧道和匝道的局部阻力系數隨風量比的增加呈單值增大的趨勢，但q1和q2小于0.4時的主隧道、匝道的局部阻力系數均為負值，此規律與文獻[13]關于矩形風管合流三通旁通管的趨勢相同，但不同于直通管的變化規律。

(a) 分流匝道

(b) 合流匝道

Fig. 8 Relationships between local resistance coefficient of confluence and distributary ramps andRe

(a) 主隧道

(b) 匝道

Fig. 9 Relationships between local resistance coefficient of confluence and distributary ramps andq

2.3 主隧道與匝道夾角的影響(工況3)

2.3.1 分流匝道

圖10示出主隧道與匝道夾角分別為0°、10°、20°、30°時，分流匝道局部阻力系數的差異性。由圖10(a)可知，當風量比小于0.5且風量比一定時，分岔角度變化對主隧道局部阻力的影響有限； 但當風量比大于0.5后，主隧道與匝道夾角變化對直隧道局部阻力系數的影響明顯增大，特別是對匝道局部阻力系數的影響尤為顯著。

(a) 主隧道

(b) 匝道

圖10分流匝道局部阻力系數與分岔角度的關系

Fig. 10 Relationships between local resistance coefficient of distributary ramp andθ

圖11為風量比q1=q2=0.5條件下的速度場計算云圖。由圖可知，隨著主隧道與匝道夾角的增大，進入主隧道的流線分布基本沒有發生變化，這意味著在q1=q2=0.5條件下，夾角變化對主隧道局部阻力系數的影響有限； 而隨著夾角的增大，氣流偏轉角度增大，氣流能快速充滿整個隧道斷面并均勻流動，匝道斷面平均壓力減小，局部阻力系數增大。

(a) 0°

(b) 10°

(c) 20°

(d) 30°

圖11風量比q1=q2=0.5時分流匝道流速分布及流線圖

Fig. 11 Flow velocity distribution and streamlines of distributary ramp whenq1=q2=0.5

2.3.2 合流匝道

圖12示出主隧道與匝道夾角分別為0°、10°、20°、30°時，合流匝道局部阻力系數的差異性。圖12(a)示出隨著主隧道與匝道夾角的增加，主隧道局部阻力系數呈單值上升趨勢； 但隨著風量比q1逐漸增大，這種影響在減弱。圖12(b)反映了風量比的變化對匝道局部阻力系數的影響不可忽略，特別是風量比大于0.5后，這種影響更明顯。

(a) 主隧道

(b) 匝道

圖12合流匝道局部阻力系數與分岔角度的關系

Fig. 12 Relationships between local resistance coefficient of confluence ramp andθ

圖13為風量比q1=0.17、q2=0.83時的流速分布及流線圖。由圖可知，隨著主隧道與匝道夾角的增大，主隧道漩渦尺寸在逐漸增大，致使主隧道的局部損失增大。在匝道風量比q2=0.83的條件下，隨著主隧道與匝道夾角的增大渦旋減弱，但是在30°夾角時渦旋又增多，這也是合流匝道局部阻力系數在0°、10°、20°、30°之間先減小后增大的原因。

2.4 局部阻力系數計算關聯式

通過上述分析可知，風量比、主隧道與匝道夾角是影響分(合)流匝道局部阻力系數的關鍵因素。基于上述分析結果，根據最小二乘法并采用Matlab軟件對計算結果進行多因素回歸分析，可得到合流匝道主隧道局部阻力系數ξ13、匝道局部阻力系數ξ23關聯式(3)—(4)； 分流匝道主隧道局部阻力系數ξ31、匝道局部阻力系數ξ32關聯式(5)—(6)。

(3)

(4)

(5)

(6)

(a) 0°

(b) 10°

(c) 20°

(d) 30°

圖13風量比q1=0.17、q2=0.83時合流匝道流速分布及流線圖

Fig. 13 Flow velocity distribution and streamlines of confluence ramp whenq1=0.17 andq2=0.83

3 結論與建議

基于長沙營盤路湘江隧道的實測、1∶8縮尺模型試驗以及CFD數值分析結果，得到以下關于城市地下道路分(合)流匝道通風阻力特性研究結果：

1)主隧道與匝道風量比、主隧道與匝道夾角是影響分(合)流匝道局部阻力系數的關鍵因素。

2)城市地下道路分(合)流匝道局部阻力變化特性，尤其是分(合)流主隧道隨風量比的變化規律完全不同于文獻[13]，所以不宜簡單套用通風管道的三通構件的參數。

3)基于長沙營盤路湘江隧道的計算結果，根據最小二乘法以及Matlab軟件進行多因素回歸分析，得出關于風量比和分岔角度的分(合)流匝道局部阻力系數計算關聯式。

4)本研究結果基于長沙市營盤路湘江隧道的結構特征，即主隧道與匝道斷面積基本相同。然而隨著各種特長復雜城市隧道的修建，主隧道和匝道的結構特性也在發生變化，所以建議進一步研究分析多個斷面積比條件下的分(合)流匝道局部阻力特性，為城市地下道路通風工程設計計算、相關標準與規范的制定提供方法參考。

[1] 李峰.公路隧道通風彎曲風道及火災排煙風口優化研究[D]. 西安： 長安大學, 2009.

LI Feng. Study of the optimization of bending-duct in road tunnel and smoke outlets′layout in fire[D]. Xi′an: Chang′an University, 2009.

[2] 溫玉輝.公路隧道排送通風系統升壓模式分析研究[J].廣東公路交通, 2011(4): 43.

WEN Yuhui. Study of pressure increasing principle of highway tunnel with exhaust and blow ventilation system [J]. Guangdong Highway Communications, 2011(4): 43.

[3] 喻映華, 雷波, 畢海權. 隧道緊急停車帶流場和局部阻力數值模擬[J]. 現代隧道技術, 2006, 43(1): 49.

YU Yinghua, LEI Bo, BI Haiquan. Study of the airflow field and the coefficient of local resistance at the emergency parking place in a tunnel[J]. Modern Tunnelling Technology, 2006, 43(1): 49.

[4] 石平, 夏永旭. 公路隧道送排風短道流場數值分析及優化研究[J].地下空間與工程學報, 2009, 5(2): 402.

SHI Ping, XIA Yongxu. Numerical analysis and optimization research of short duct flow field of outlet and inlet in highway tunnel[J].Chinese Journal of Underground Space and Engineering, 2009, 5(2): 402.

[5] 馬佳. 長大公路隧道通風系統局部結構優化研究[D]. 西安： 長安大學, 2015.

MA Jia. Study of the optimization of the long highway tunnel ventilation system on local structure[D]. Xi′an: Chang′an University, 2015.

[6] 茍紅松, 吳元金, 李永生, 等. 隧道施工通風分風三通位置選擇研究[J]. 隧道建設, 2016, 36(4): 384.

GOU Hongsong, WU Yuanjin, LI Yongsheng, et al. Selection of location of ventilation t-junction in tunnel construction[J]. Tunnel Construction, 2016, 36(4): 384.

[7] COSTA N P, MAIA R, PROENCA M F, et al. Edge effects on the flow characteristics in a 90 degree junction[J]. Journal of Fluids Engineering, 2006, 128: 1204.

[8] ABDULWAHHAB M, INJETI N K, DAKHI S F. Numerical prediction of pressure loss of fluid in a T-junction[J]. International Journal of Energy and Environment, 2013, 4(2): 253.

[9] 鞏啟濤, 楊俊紅, 韓奎, 等. 大管徑T型三通數值模擬及局部阻力特性分析[J]. 動力工程學報, 2016, 36(9): 753.

GONG Qitao, YANG Junhong, HAN Kui, et al. Characteristic analysis of the flow and local resistance in large pipe tees[J]. Journal of Chinese Society of Power Engineering, 2016, 36(9): 753.

[10] 石喜, 呂宏興, 朱德蘭, 等. PVC三通管水流阻力與流動特征分析[J]. 農業機械學報, 2013, 44(1): 73.

SHI Xi, LYU Hongxing, ZHU Delan, et al. Flow resistance and characteristics of PVC tee pipes[J]. Transactions of the Chinese Society of Agricultural Machinery, 2013, 44(1): 73.

[11] 陳江林, 呂宏興, 石喜, 等. T型三通管水力特性的數值模擬與試驗研究[J]. 農業工程學報, 2012, 28(5): 73.

CHEN Jianglin, LYU Hongxing, SHI Xi, et al. Numerical simulation and experimental study of hydrodynamic characteristics of T-type pipes[J]. Transactions of the Chinese Society of Agricultural Engineering, 2012, 28(5): 73.

[12] 陳偉業, 呂宏興, 石喜, 等. 等徑 PVC 三通管局部水頭損失系數試驗研究[J]. 灌溉排水學報, 2013, 32 (1): 128.

CHEN Weiye, LYU Hongxing, SHI Xi, et al. Experimental study of local loss parameter of diameter PVC pipe tee[J]. Journal of Irrigation and Drainage, 2013, 32 (1): 128.

[13] 陸耀慶. 實用供熱空調設計手冊[M]. 北京： 中國建筑工業出版社， 2008.

LU Yaoqing. Practical design handbook for heating and air conditioning[M]. Beijing: China Architecture & Building Press, 2008.

[14] 公路隧道通風設計細則: JTG/T D70/2-02—2014[S].北京： 人民交通出版社， 2014.

Guidelines for design of ventilation of highway tunnels: JTG/T D70/2-02—2014[S]. Beijing: China Communications Press, 2014.

[15] 方磊. 長大公路隧道通風物理模型試驗研究[D]. 西安： 長安大學, 2005.

FANG Lei. The physical model experiment research of long highway tunnel ventilation[D]. Xi′an: Chang′an University, 2005.

VentilationResistanceCharacteristicsofConfluenceandDistributaryRampsofUrbanTunnel

YUAN Haoting1, CHEN Chao1, *, LI Qiong1, 2, WANG Luyao1, LI Yaru1

(1.CollegeofArchitectureandCivilEngineering,BeijingUniversityofTechnology,Beijing100124,China; 2.ArchitecturalEngineeringCollege,NorthChinaInstituteofScienceandTechnology,Beijing101601,China)

The ventilation resistance characteristics of confluence and distributary ramps of Xiangjiang River-crossing Tunnel on Yingpan Road in Changsha is studied by field monitoring, scale model experiment and software CFD numerical simulation; and then the influences of Reynolds number, air volume ratio and bifurcation angle on resistance characteristic of confluence and distributary ramps of the tunnel are analyzed so as to know the air flowing characteristics of urban underground passageway with bifurcated ramps. The results show that: 1) The air volume ratio and bifurcation angle are the key factors on the local resistance coefficients of confluence and distributary ramps. 2) The parameters of triplet joint component of ventilation pipe cannot be applied to characterize the variation of local resistance of confluence and distributary ramps. 3) The multiple-factor regression analysis of calculation results is carried out by least square method and software Matlab, which provides the expression equation of local resistance coefficient of confluence and distributary ramps of Xiangjiang River-crossing Tunnel on Yingpan Road in Changsha. The study results can provide reference for analysis of resistance characteristics of ventilation system and optimization of ventilation engineering design of urban underground passageway with complex structure.

urban bifurcate tunnel; confluence and distributary ramp; local resistance coefficient; scale model experiment; numerical simulation

2017-03-01;

2017-05-04

國家自然科學基金項目(51378024)； 北京市自然科學基金項目(8162006)

袁浩庭(1990—)，男，河北滄州人，北京工業大學暖通專業在讀碩士，研究方向為隧道通風。E-mail： yhtbjut1990@163.com。*通信作者： 陳超， E-mail： chenchao@bjut.edu.cn。

10.3973/j.issn.2096-4498.2017.11.008

U 45

A

2096-4498(2017)11-1409-08