錢喜玲，閆小燕，趙江平

(1.西安建筑科技大學(xué) 材料與礦資學(xué)院，陜西 西安，710055；2.陜西省西安市公安消防支隊(duì)，陜西 西安，710055)

0　引言

地下綜合管廊作為21世紀(jì)新型城市市政基礎(chǔ)設(shè)施建設(shè)的重要標(biāo)志之一，現(xiàn)已被越來(lái)越多的城市采用。地下綜合管廊(Utility tunnel)是指建于城市地下用于容納兩類及以上城市工程管線的構(gòu)筑物及附屬設(shè)施[1]，地下綜合管廊內(nèi)敷設(shè)有自來(lái)水、污水、熱力、電力、電信、天然氣等各種管線。使用綜合管廊帶來(lái)便利的同時(shí)也帶來(lái)了一定的危險(xiǎn)性，特別是電纜管線和天然氣管線具有巨大的火災(zāi)隱患。地下綜合管廊一旦發(fā)生火災(zāi)，煙氣會(huì)沿著管艙迅速蔓延，產(chǎn)生的煙氣不能及時(shí)排出，熱量大量聚集，使得管艙內(nèi)溫度快速升高，危險(xiǎn)性極大，同時(shí)該類設(shè)施的建設(shè)費(fèi)用遠(yuǎn)高于傳統(tǒng)的直埋式鋪設(shè)的成本[2]，因此一旦在管艙內(nèi)出現(xiàn)異常情況，勢(shì)必造成巨大的經(jīng)濟(jì)損失。

為了預(yù)防地下綜合管廊事故，許多學(xué)者已進(jìn)行了相關(guān)研究。林俊等[3]采用FDS分析軟件，對(duì)熱源恒定情況下，不同通風(fēng)風(fēng)速和不同防火分區(qū)長(zhǎng)度地下綜合管廊火災(zāi)時(shí)的煙氣流動(dòng)特性進(jìn)行研究，最終對(duì)城市綜合管廊建設(shè)中的防火分區(qū)設(shè)置和通風(fēng)設(shè)備的安裝提出建議；Kim HS[4]等運(yùn)用CFD模擬軟件，建立以體積熱源為火源的圓形和矩形綜合管廊尺度，對(duì)空氣流速及煙氣流量進(jìn)行對(duì)比分析，得知圓形綜合管廊模型煙氣溫度高于矩形；李文婷[5]根據(jù)測(cè)量電纜質(zhì)量損失速率和熱釋放速率參數(shù)，運(yùn)用FDS分析軟件，對(duì)地下綜合管廊電纜火災(zāi)煙氣蔓延過(guò)程和熱量傳遞過(guò)程進(jìn)行模擬分析，得到電纜管艙內(nèi)部分溫度高達(dá) 1 000℃時(shí)，可能引起纜線轟燃和管艙結(jié)構(gòu)變形；方自虎[6]以深圳地下綜合管廊為研究對(duì)象，用甲烷代替天然氣，采用向上和水平2種管道泄漏口開(kāi)口方向，利用Fluent軟件模擬了天然氣管艙內(nèi)甲烷泄漏的擴(kuò)散過(guò)程，得到距離泄漏口達(dá)到一定的距離后，甲烷濃度的變化與泄漏口的開(kāi)口方向無(wú)關(guān)；趙永昌等[7]在隧道模型基礎(chǔ)上，搭建了1∶3∶6的小尺寸地下綜合管廊實(shí)體模型，對(duì)電力艙室火災(zāi)通過(guò)油池火試驗(yàn)研究發(fā)現(xiàn)，不同火源功率下，煙氣溫度均呈現(xiàn)冪函數(shù)衰減，火源功率較大時(shí)，溫度衰減梯度也較大，電力艙內(nèi)火災(zāi)煙氣蔓延過(guò)程中存在煙氣分層現(xiàn)象；李朝陽(yáng)等[8]利用Fluent軟件對(duì)埋地含硫天然氣管道穿孔泄漏后的甲烷、硫化氫氣體的擴(kuò)散進(jìn)行了數(shù)值模擬。

通過(guò)閱讀大量文獻(xiàn)發(fā)現(xiàn)，國(guó)內(nèi)外學(xué)者對(duì)地下綜合管廊的研究主要集中于電纜火災(zāi)，對(duì)天然氣泄漏擴(kuò)散規(guī)律研究甚少。因此，本文采用Fluent模擬軟件對(duì)某地下綜合管廊天然氣泄漏后在管艙內(nèi)的擴(kuò)散規(guī)律進(jìn)行數(shù)值仿真模擬，根據(jù)天然氣爆炸下限濃度，將質(zhì)量濃度轉(zhuǎn)化為探頭報(bào)警響應(yīng)時(shí)間，為預(yù)防天然氣泄漏后發(fā)生火災(zāi)提供一定的參考依據(jù)。

1　物理模型

以西安市某地下綜合管廊天然氣管艙為研究對(duì)象，其簡(jiǎn)化的物理模型如圖1所示。為體現(xiàn)約束條件對(duì)泄漏擴(kuò)散規(guī)律的影響，模型中保留X方向和Y方向的約束。在機(jī)械通風(fēng)條件下，圓形泄漏孔徑設(shè)置在模型下方中間位置，直徑為20 mm。由于城鎮(zhèn)天然氣中甲烷的組分一般大于95%，所以模擬時(shí)以甲烷為介質(zhì)，天然氣管道內(nèi)氣體平均溫度T1=300 K[9]，模擬時(shí)忽略泄漏氣體與空氣溫度的差異，管艙內(nèi)溫度Ta=300 K壓力，Pa=101 325 Pa。根據(jù)《石油化工可燃?xì)怏w和有毒氣體檢測(cè)報(bào)警設(shè)計(jì)規(guī)范》(GB/50493-2009)中的規(guī)定，可燃?xì)怏w處于封閉或局部通風(fēng)不良的半敞開(kāi)廠房?jī)?nèi)，每隔15 m可設(shè)置一臺(tái)探測(cè)器，且探測(cè)器距其所覆蓋范圍內(nèi)的任意釋放源不宜大于7.5 m[10]，故選擇天然氣管艙頂部每15 m設(shè)置一個(gè)天然氣報(bào)警探測(cè)器，以監(jiān)測(cè)探測(cè)器報(bào)警響應(yīng)時(shí)間。

圖1　天然氣管艙模型Fig.1　Natural gas pipe model

1.1　模型簡(jiǎn)化

為了簡(jiǎn)化計(jì)算過(guò)程，根據(jù)管艙內(nèi)的實(shí)際布置情況，在滿足工程精度的條件下可忽略一些次要因素，具體假設(shè)如下：

1)將計(jì)算域進(jìn)行簡(jiǎn)化，忽略天然氣管道本身對(duì)泄漏氣體擴(kuò)散的影響；

2)天然氣泄漏口處的泄漏壓力恒定[11]；

3)在數(shù)值計(jì)算中忽略空氣粘性的影響[12]；

4)假定泄漏氣體為理想氣體，氣體之間不發(fā)生化學(xué)反應(yīng)；

5)泄漏過(guò)程中，溫度恒定且與外界無(wú)熱量交換[13]。

1.2　邊界條件

入口條件：空氣入口邊界選用速度入口，地下綜合管廊天然氣管艙內(nèi)平時(shí)通風(fēng)量為每小時(shí)6次，事故通風(fēng)排風(fēng)量為每小時(shí)12次[14]，換算得到進(jìn)口風(fēng)速Vair=3.5 m/s，天然氣泄漏口選用壓力入口，泄漏壓力為P1=103.3 kPa和P2=200.0 kPa。

出口條件：天然氣出口選用壓力出口，出口壓力為101 325 Pa。

壁面條件：固體壁面為無(wú)滑移條件。

Fluent模擬采用非耦合隱式求解器，由于Realizablek-ε湍流模型已被有效地用于各種不同類型的流動(dòng)模擬[15]，其中包含有射流和混合流等，所以計(jì)算選用Realizablek-ε湍流模型，天然氣在大氣中的擴(kuò)散過(guò)程實(shí)際上是天然氣在空氣中的運(yùn)動(dòng)，涉及到化學(xué)組分的傳輸[16]，故選擇組分傳輸模型(Species transport)，并定義組分為甲烷-空氣。

2　模擬結(jié)果與分析

2.1　氣體泄漏初期擴(kuò)散過(guò)程

根據(jù)《城市綜合管廊工程技術(shù)規(guī)范》(GB/T50086-2015)中的規(guī)定，地下綜合管廊中天然氣報(bào)警濃度設(shè)定值(上限值)不應(yīng)大于其爆炸下限值(體積分?jǐn)?shù))的20%[14]。天然氣的爆炸極限為5%～15%(體積分?jǐn)?shù))，換算成質(zhì)量分?jǐn)?shù)爆炸極限為2.7%～8.6%，選擇爆炸下限濃度的20%作為天然氣報(bào)警的臨界值，換算成質(zhì)量分?jǐn)?shù)天然氣報(bào)警上限約為0.6%。以下針對(duì)天然氣泄漏后最先的濃度探測(cè)器報(bào)警時(shí)，不同泄漏壓力下天然氣在管艙中的運(yùn)移規(guī)律，泄漏壓力分別為103.3 kPa和200.0 kPa時(shí)，天然氣泄漏后濃度分布情況如圖2所示。

圖2　壓力不同時(shí)管道泄漏天然氣擴(kuò)散規(guī)律Fig.2　Affected area under different operating pressure of the pipeline

由圖2分析可以得出，泄漏初期，天然氣管道內(nèi)的壓力大于管艙內(nèi)的壓力，氣體向管艙內(nèi)的泄漏呈射流態(tài)，氣體進(jìn)入管艙后氣流快速?zèng)_擊管艙頂棚，受到壁面的約束，氣流將在沖擊頂棚后沿水平方向向泄漏口左右兩側(cè)運(yùn)動(dòng)。從天然氣濃度分布等值線圖中看出，管艙內(nèi)氣體質(zhì)量分?jǐn)?shù)從0.6%至8.6%之間的區(qū)域均已在爆炸極限內(nèi)，在此區(qū)域內(nèi)若氧氣含量充足，遇火源點(diǎn)燃將會(huì)發(fā)生爆炸事故；泄漏口處氣體質(zhì)量分?jǐn)?shù)均高于32.6%，若火源在泄漏口處，氣體高速噴射擴(kuò)散遇到火源點(diǎn)燃將發(fā)生燃燒事故，這是因?yàn)樘烊粴鈬娚渌俣容^快，流動(dòng)性較大，不足以發(fā)生爆炸。由于在管艙最左端的頂部設(shè)有機(jī)械進(jìn)風(fēng)口，天然氣偏向管艙右端的排風(fēng)口側(cè)進(jìn)行擴(kuò)散，泄漏壓力為103.3 kPa時(shí)，最先的濃度探測(cè)器報(bào)警時(shí)間是2.15 s，泄漏壓力為200.0 kPa時(shí)，最先的濃度探測(cè)器報(bào)警時(shí)間是0.45 s，泄漏壓力對(duì)探測(cè)器的響應(yīng)時(shí)間有顯著影響，響應(yīng)時(shí)間隨泄漏壓力的增大而減小。

泄漏壓力分別為103.3 kPa和200.0 kPa時(shí)，天然氣泄漏后從泄漏口處向管艙右端排風(fēng)口側(cè)的擴(kuò)散距離與時(shí)間的變化規(guī)律如圖3所示。

圖3　不同泄漏壓力下擴(kuò)散距離與時(shí)間變化關(guān)系Fig.3　Spread distance vs．spread time under different operating pressure

圖3中的結(jié)果表明：泄漏時(shí)間相同時(shí)，泄漏壓力對(duì)氣體擴(kuò)散距離有顯著影響，且氣體擴(kuò)散距離隨泄漏壓力的增大而增大；對(duì)所得數(shù)據(jù)進(jìn)行擬合可知，泄漏壓力分別為103.3 kPa和200 kPa時(shí)，泄漏氣體擴(kuò)散距離與時(shí)間的函數(shù)關(guān)系分別為：y=2.12 158+2.03 699x(x≥1)和y=9.75 579+2.87 278x(x≥1)，故泄漏壓力相同時(shí)，氣體擴(kuò)散距離與時(shí)間呈一次函數(shù)，即y=a+bx(x≥1)當(dāng)其曲線斜率表示氣體在管艙中的擴(kuò)散速度，隨著泄漏壓力的增大氣體擴(kuò)散速度增大。

2.2　探測(cè)器報(bào)警響應(yīng)時(shí)間

天然氣泄漏后，天然氣濃度探測(cè)器的報(bào)警響應(yīng)時(shí)間是安全控制的重要數(shù)據(jù)，因此，從地下綜合管廊安全控制來(lái)說(shuō)，對(duì)天然氣泄漏后探測(cè)器報(bào)警響應(yīng)時(shí)間的研究至關(guān)重要。管艙中每相隔15 m設(shè)置一個(gè)濃度探測(cè)器，且選擇爆炸下限濃度的20%作為探測(cè)器響應(yīng)臨界點(diǎn)，泄漏壓力分別為103.3 kPa和200.0 kPa時(shí)，天然氣泄漏后管艙中的探測(cè)器報(bào)警響應(yīng)時(shí)間如圖4所示。

圖4　不同泄漏壓力與探測(cè)器響應(yīng)時(shí)間關(guān)系Fig.4　The different operating pressure vs. the detector response time

圖4中的結(jié)果表明：泄漏壓力為103.3 kPa時(shí)，最先的濃度探測(cè)器在天然氣泄漏2.15 s時(shí)響應(yīng)，且60 m處的探測(cè)器沒(méi)有響應(yīng)，泄漏壓力為200.0 kPa時(shí)，最先的濃度探測(cè)器響應(yīng)時(shí)間是0.45 s，且泄漏口上游距離泄漏口60 m處的探測(cè)器在天然氣泄漏12.1 s時(shí)響應(yīng)，而下游距離泄漏口90 m處的探測(cè)器在2.75 s時(shí)已經(jīng)響應(yīng)，可見(jiàn)受機(jī)械進(jìn)風(fēng)的影響，泄漏口下游的探測(cè)器比上游的探測(cè)器先響應(yīng)，且泄漏口下游的濃度探測(cè)器在一定的時(shí)間范圍內(nèi)均會(huì)響應(yīng)。

圖5　報(bào)警濃度與探測(cè)器響應(yīng)時(shí)間關(guān)系Fig.5　The alarm concentration vs. the response time of detector

選擇爆炸下限濃度的10%和20%作為探測(cè)器探頭響應(yīng)臨界點(diǎn)，泄漏壓力分別為103.3 kPa和200.0 kPa時(shí)，天然氣泄漏后報(bào)警濃度與探測(cè)器響應(yīng)時(shí)間的關(guān)系如圖5所示。圖5中結(jié)果表明：受機(jī)械進(jìn)風(fēng)影響，擴(kuò)散氣體沿泄漏口下游方向遷移，泄漏口上游方向45 m處以前所設(shè)置的濃度探測(cè)均未響應(yīng)；當(dāng)泄漏壓力相同時(shí)，報(bào)警濃度設(shè)置對(duì)探測(cè)器響應(yīng)時(shí)間的影響極小，同一位置探測(cè)器的報(bào)警響應(yīng)時(shí)間不超過(guò)1 s，故管艙中每15 m設(shè)置一個(gè)探測(cè)器和選擇爆炸下限的20%作為探測(cè)器響應(yīng)的臨界濃度值是合理的。

3　結(jié)論

1)以探測(cè)器保護(hù)半徑為7.5 m核算，當(dāng)泄漏壓力為103.3 kPa時(shí)，濃度探測(cè)器報(bào)警時(shí)間是2.15 s，當(dāng)泄漏壓力為200.0 kPa時(shí)，濃度探測(cè)器報(bào)警時(shí)間是0.45 s，報(bào)警響應(yīng)時(shí)間隨著泄漏壓力的增大而減小，在常規(guī)中壓輸出壓力下，響應(yīng)時(shí)間最大值為2.15 s。

2)對(duì)所得數(shù)據(jù)進(jìn)行擬合可知，泄漏壓力相同時(shí)，氣體擴(kuò)散距離與時(shí)間呈正相關(guān)關(guān)系，其曲線斜率表示氣體在管艙中的擴(kuò)散速度。

3)當(dāng)泄漏壓力相同時(shí)，報(bào)警濃度設(shè)置對(duì)探測(cè)器響應(yīng)時(shí)間的影響極小，同一位置探測(cè)器的報(bào)警響應(yīng)時(shí)間不超過(guò)1 s，故管艙中選擇爆炸下限的20%作為探測(cè)器響應(yīng)的臨界濃度值是合理的。

[1]Xing-hui Zhang, Yu-xi Guan, Yu-fan Liao. Fire risk analysis and prevention of urban comprehensive pipeline corridor[J]. Procedia Engineering 2016, 7(6): 463 -468．

[2]錢七虎，陳曉強(qiáng)．國(guó)內(nèi)外地下綜合管線廊道發(fā)展的現(xiàn)狀、問(wèn)題及對(duì)策[J]．地下空間與工程學(xué)報(bào)，2007(2)：191-194.

QIAN Qihu, CHEN Xiaoqiang. Situation problems and counter measures of utility tunnel 'development in China and Abroad[J]. Chinese Journal of Underground Space andEngineering, 2007(2): 191-194.

[3]林俊，叢北華，韓新，等．基于模擬CFD分析的城市綜合管廊火災(zāi)特性研究[A]．上海防災(zāi)救災(zāi)研究所20周年慶典會(huì)議研究短文集[C]，2009．

LIN Jun, CONG Beihua, HAN Xin, et al. Study on fire characteristics of urban integrated corridor based on simulated CFD analysis[A].Study on the 20th Anniversary of the Shanghai Disaster Relief Research Institute[C], 2009.

[4]Kim HS, KIMYJ. Characteristics of smoke con-centration profiles with underground utility tunnel fire[J]. 2005, 10(1): 94-98．

[5]李文婷．綜合管溝電纜火災(zāi)數(shù)值模擬研究[D]．北京：首都經(jīng)濟(jì)貿(mào)易大學(xué)，2012．

[6]方自虎，王家遠(yuǎn)．共同溝內(nèi)燃?xì)鈹U(kuò)散規(guī)律的數(shù)值分析[J]．武漢大學(xué)學(xué)報(bào)，2009，42(2)：215-218．

FANG Zihu, WANG Jiayuan. Mumerical analysis of gas diffusion properties in utility tunnel[J].Engineering Journal of Wuhan University, 2009, 42(2): 215-218．

[7]趙永昌，朱國(guó)慶，高云驥. 城市地下綜合管廊火災(zāi)煙氣溫度場(chǎng)研究[J]．消防科學(xué)與技術(shù)，2017，36(1)：37-40．

ZHAO Yongchang, ZHU Guoqing, GAO Yunji. Study on flue gas temperature field of urban underground comprehensive corridor[J]. Fire Science and Technology, 2017,36(1): 37-40.

[8]李朝陽(yáng)，馬貴陽(yáng)，徐柳. 架空及埋地天然氣管道泄漏擴(kuò)散數(shù)值研究[J]. 中國(guó)安全科學(xué)生產(chǎn)技術(shù)，2011，7(7):66-69.

LI Zhaoyang,MA Guiyang,Xu Liu.Numerical study on the diffusion concentration of overhead and buried nataral gas pipeline[J]. Journal of Safety Science and Technology, 2011, 7(7): 66-69.

[9]于菲菲，王世杰．受限空間泄漏孔徑對(duì)氣體擴(kuò)散影響的模擬研究[J]．安全與環(huán)境學(xué)報(bào),2015,15(4):159-162．

YU Feifei, WANG Shijie. Simulation study on the effect of leakage aperture on gas diffusion in confined spaces[J].Journal of Safety and Environment, 2015,15(4): 159-162.

[10]中華人民共和國(guó)住房和城鄉(xiāng)建設(shè)部. 石油化工可燃?xì)怏w和有毒氣體檢測(cè)報(bào)警設(shè)計(jì)規(guī)范：GB/50493-2009[S]. 北京：中國(guó)計(jì)劃出版社，2009.

[11]Veynante D, Vervish L. Turbulent combustion modeling[J]. Progress in Energy and Combustion Science,2002(28): 193-266.

[12]孟超，趙晶．居室天然氣泄漏擴(kuò)散過(guò)程仿真研究[J]．中國(guó)安全生產(chǎn)科學(xué)技術(shù)，2011，7(5)：153-158．

MENG Chao，ZHAO Jing. Study on simulation of leakage and diffusion process for natural gas in residence[J]. Journal of Safety Science and Technology, 2011,7(5):153-158.

[13]孫恩吉，李紅果，王敏. 基于Realizablek-ε湍流模型的氨氣泄漏數(shù)值模擬研究[J]. 中國(guó)安全生產(chǎn)科學(xué)技術(shù)，2017，13(2)：114-118.

SUN Enji, LI Hongguo, WANG Min. Study on numerical simulation of ammonia leakage based on realizablek-εturbulence model.[J]. Journal of Safety Science and Technology, 2017, 13(2):114-118.

[14] 中華人民共和國(guó)住房和城鄉(xiāng)建設(shè)部. 城市綜合管廊工程技術(shù)規(guī)范：GB/T50838-2015[S]. 北京：中國(guó)計(jì)劃出版社，2015.

[15]王福軍．計(jì)算流體動(dòng)力學(xué)分析—CFD軟件原理與應(yīng)用[M]．北京：清華大學(xué)出版社，2004．

[16]黃小美，張婧，周陽(yáng),等．街道天然氣管道泄漏擴(kuò)散的數(shù)值模擬分析[J]．煤氣與熱力，2016，36(2)：32-38．

HUANG Xiaomei, ZHANG Jing, ZHOU Yang, et al. Numerical simulation analysis of leakage and diffusion of natural gas pipeline in street [J].Gas and Heat, 2016，32(2):32-38．