王敬林，李秀地，王曉波

（1后勤工程學(xué)院軍事建筑工程系 重慶 400041 2重慶博其樂(lè)土木工程技術(shù)咨詢有限公司 重慶 400069）

現(xiàn)代城市的快速發(fā)展和天然資源的不可重復(fù)始終是一個(gè)矛盾。壓縮天然氣（CNG）是一種最理想的車用替代能源，它具有成本低、效益高、無(wú)污染、使用安全便捷等特點(diǎn)，正日益顯示出強(qiáng)大的發(fā)展?jié)摿Α．?dāng)前，由于加氣站數(shù)量不足等原因，許多城市均出現(xiàn)車輛“排長(zhǎng)隊(duì)”加氣的現(xiàn)象，因此都在積極規(guī)劃新建加氣站項(xiàng)目。同時(shí)，CNG加氣站一般修建在城市交通主干道邊上，有時(shí)與輕軌或者地鐵之間相距很近，如果處理不好會(huì)給城市帶來(lái)一定的安全隱患。近年來(lái)，CNG加氣站發(fā)生爆炸的事件屢有發(fā)生，給人們帶來(lái)了深刻的教訓(xùn)，如2005年10月8日，重慶市九龍坡區(qū)石橋鋪白馬凼一加氣站氣井漏氣，并發(fā)生爆炸，一人被炸重傷，附近上千人疏散，消防隊(duì)出動(dòng)5輛消防車，耗用5小時(shí)才排險(xiǎn)成功。為預(yù)防事故的發(fā)生或者減少事故帶來(lái)的損失，2009年8月31日，成都市在三環(huán)路外側(cè)的通能清江加氣站內(nèi)進(jìn)行了最大規(guī)模的CNG站事故應(yīng)急演練。

目前有關(guān)CNG加氣站設(shè)計(jì)、風(fēng)險(xiǎn)安全分析與預(yù)防的文獻(xiàn)較多[1][2][3][4]，但針對(duì)CNG加氣站的爆炸對(duì)鄰近工程的影響的文章卻很少，本文針對(duì)重慶市某新建CNG加氣站地處毗鄰的輕軌三號(hào)線隧洞安全控制保護(hù)線范圍內(nèi)，利用有限元數(shù)值模擬，探討了CNG加氣站在不利工況下發(fā)生爆炸對(duì)鄰近輕軌隧洞影響的數(shù)值分析與評(píng)價(jià)，為類似工程處理提供借鑒。

1 工程概況

為緩解重慶市公交車加氣難問(wèn)題，重慶市經(jīng)濟(jì)委員會(huì)立項(xiàng)批復(fù)同意由重慶燃?xì)猓瘓F(tuán)）有限責(zé)任公司和重慶市公交站場(chǎng)集團(tuán)合作建設(shè)江北龍頭寺公交站場(chǎng)CNG加氣站。該項(xiàng)目地處龍頭寺公交站場(chǎng)內(nèi)，距輕軌三號(hào)線隧洞邊線31.2m，處于輕軌交通保護(hù)控制線內(nèi)[1]。

龍頭寺公交站場(chǎng)CNG加氣站設(shè)計(jì)參數(shù)如下：進(jìn)氣壓力：0.2MPa；天然氣日處理量：（8.0～10.0）×104m3/d；儲(chǔ)氣能力：采用儲(chǔ)氣井儲(chǔ)氣方式，分為高、中、低三個(gè)壓力段分別儲(chǔ)存天然氣（高壓井1口，中壓井2口，低壓井3口），每口井深150m，儲(chǔ)氣量4500m3，儲(chǔ)氣井共占地9m2。加氣站與輕軌隧洞平面位置關(guān)系如圖1所示。

圖1 加氣站與輕軌隧洞平面位置關(guān)系示意圖

CNG加氣站方案按照國(guó)家標(biāo)準(zhǔn)GB50156《汽車加氣站設(shè)計(jì)與施工規(guī)范》、GB50028《城鎮(zhèn)燃?xì)庠O(shè)計(jì)規(guī)范》設(shè)計(jì)，并采用了一定的技術(shù)安全防護(hù)措施，如將工藝區(qū)改為半地下式，工藝區(qū)屋頂上側(cè)開(kāi)設(shè)10m×25m寬的框洞，泄壓面積滿足《建筑設(shè)計(jì)防火規(guī)范》的泄壓要求，同時(shí)采用500mm厚的鋼筋混凝土防爆隔離墻，上至地表，下至車庫(kù)地坪下1m，設(shè)備控制線纜采用埋地方式引入等，剖面關(guān)系如圖2所示。

圖2 加氣站與輕軌隧洞剖面位置關(guān)系示意圖

2 場(chǎng)地巖土工程特性

擬建場(chǎng)地位于重慶市龍頭寺公交站場(chǎng)內(nèi)，屬構(gòu)造剝蝕殘丘地貌，地形開(kāi)闊、平坦。場(chǎng)地地質(zhì)構(gòu)造位于龍王洞背斜南西翼，巖層產(chǎn)狀為245°∠17°。巖體中主要發(fā)育兩組裂隙，綜合分析基巖裂隙不發(fā)育。

地層巖土特征：擬建場(chǎng)地地層主要有第四系素填土、南側(cè)分布有粉質(zhì)粘土、下伏基巖為侏羅系中統(tǒng)沙溪廟組砂巖和砂質(zhì)泥巖。中風(fēng)化巖體較完整，基本質(zhì)量等級(jí)為Ⅳ級(jí)，巖體為整體塊狀結(jié)構(gòu)。

根據(jù)提供的地質(zhì)報(bào)告，計(jì)算采用的巖土和材料物理力學(xué)參數(shù)見(jiàn)表1。

3 CNG等效TNT換算

CNG泄漏出的氣體燃料和空氣混合形成可燃混合物，在無(wú)限大氣中擴(kuò)展，形成很大面積的可燃蒸氣云，一旦遇到點(diǎn)火源，此蒸氣云即可發(fā)生大面積的無(wú)約束可燃蒸氣云爆炸。要精確計(jì)算壓縮天然氣的爆炸能量是比較困難的。雖然容器內(nèi)可燃?xì)怏w的量已知，而且在容器爆炸時(shí)又幾乎全部流出，但由于這些氣體一般以球狀或其它形態(tài)在空間擴(kuò)散，只有外圍一部分可燃?xì)怏w與大氣中的氧混合形成爆炸性氣體，所以并不是全部可燃?xì)怏w都參與反應(yīng)。參與反應(yīng)的可燃?xì)怏w量的多少與許多因素有關(guān)，如容器周圍的氣流情況、氣體的爆炸極限范圍和出現(xiàn)火源的時(shí)間等。因此，一般只能是估算，即假定參與爆炸反應(yīng)氣體所占的百分比，然后按這些可燃?xì)怏w的燃燒熱計(jì)算其爆炸能量。

表1 巖土物理力學(xué)參數(shù)計(jì)算取值

計(jì)算這類蒸氣云爆炸效應(yīng)有以下幾類方法：第一類是經(jīng)驗(yàn)方法，其中最重要的是TNT當(dāng)量法；第二類是理論方法，其中包括自模擬解及作為其簡(jiǎn)化近似的解析方法；第三類是球形和半球形的數(shù)值模型；第四類是數(shù)值方法，如AUTOREAGAS商用軟件。

經(jīng)驗(yàn)方法TNT當(dāng)量法的爆炸能量是用TNT當(dāng)量來(lái)表示的，其最大優(yōu)點(diǎn)是簡(jiǎn)便易行，因而長(zhǎng)期以來(lái)在爆炸波效應(yīng)預(yù)報(bào)以及爆炸事故調(diào)查中一直被廣泛應(yīng)用。雖然存在一定的誤差，但在稍遠(yuǎn)的距離上，相同能量的壓縮天然氣爆炸與TNT爆炸所產(chǎn)生的超壓還是相近似的。其基本概念是將無(wú)約束蒸氣云所產(chǎn)生的爆炸波等同于相同能量的TNT爆轟所產(chǎn)生的爆炸波，如式（1）所示：

式中：WTNT為可燃?xì)怏w蒸氣云的TNT當(dāng)量，kg；

α為可燃?xì)怏w蒸氣云的TNT當(dāng)量系數(shù)，一般取0.04；

Wf為蒸氣云爆炸中可燃?xì)怏w總質(zhì)量，kg；

Qf為可燃?xì)怏w的燃燒熱，本文中取值為55.164MJ/kg；

QTNT為TNT的燃燒熱，一般取4.52MJ/kg。

由于確定安全距離應(yīng)盡量保守，現(xiàn)假設(shè)所儲(chǔ)存的氣體（Vf=4500m3）全部泄漏到敞開(kāi)空間，與空氣混合形成爆炸混合物，發(fā)生蒸氣云地面半球爆炸。

CNG的主要成分是甲烷 （CH4），比重為0.55，密度為0.6875kg/m3。將數(shù)據(jù)代入上式，則CNG無(wú)約束可燃蒸氣云的等效TNT量為式（2）所示：

WTNT=0.04X4500X0.6875X55.164/4.52=1510kg （2）

4 數(shù)值分析

采用程序：目前應(yīng)用較廣的是美國(guó)勞倫斯利弗莫爾國(guó)家實(shí)驗(yàn)室（Lawrence Livermore National Lab.）的LS-DYNA計(jì)算分析程序，LS-DYNA是一個(gè)大型通用顯式動(dòng)力有限元程序，計(jì)算功能更為強(qiáng)大，本文即采用該軟件進(jìn)行數(shù)值模擬計(jì)算。

4.1 計(jì)算條件

計(jì)算模型：為建模方便，將等效TNT炸藥看做立方體裝藥，并假設(shè)研究對(duì)象以處理區(qū)中心截面對(duì)稱。為保守起見(jiàn)，考慮炸藥在處理區(qū)地表面爆炸，并按平面應(yīng)變問(wèn)題進(jìn)行數(shù)值計(jì)算，計(jì)算模型如圖3所示。

圖3 計(jì)算模型

計(jì)算所采用的單位為g-cm-μs-Mbar（注：1Mbar=105 MPa）。模型尺寸如圖2所示。

單元?jiǎng)澐郑簡(jiǎn)卧愋筒捎肧OLID164，單元形狀為規(guī)則四邊形。單元尺寸為10cm，共計(jì)300155個(gè)單元。由于爆炸數(shù)值模擬過(guò)程中涉及炸藥爆轟產(chǎn)物、空氣、鋼筋混凝土及巖土等介質(zhì)的大變形響應(yīng)，為避免炸藥爆轟流體網(wǎng)格在計(jì)算過(guò)程中的畸變，對(duì)炸藥和空氣單元采用歐拉算法，每個(gè)時(shí)步執(zhí)行1次輸運(yùn)計(jì)算。采用Van Leer+HIS（二階精度）方法進(jìn)行輸運(yùn)計(jì)算，該方法是目前LS-DYNA中最為精確的方法。對(duì)鋼筋混凝土和巖土采用拉格朗日算法，對(duì)歐拉介質(zhì)和拉格朗日介質(zhì)采用流固耦合算法。鋼筋混凝土單元、巖土單元之間的接觸采用共節(jié)點(diǎn)耦合在一起。

邊界條件：在對(duì)稱邊界上采用對(duì)稱邊界條件。為模擬半無(wú)限巖土介質(zhì)，在巖土的右邊界和下邊界采用非反射邊界條件，地表采用自由邊界條件。空氣邊界的范圍用這樣的方法確定，即從炸藥起爆開(kāi)始到?jīng)_擊波到達(dá)邊界的持續(xù)時(shí)間，應(yīng)與所研究點(diǎn)的壓力-時(shí)間相比足夠地長(zhǎng)。

材料及其參數(shù)：本次爆炸數(shù)值模擬過(guò)程中涉及炸藥、空氣、鋼筋混凝土及巖土等材料。

空氣模型簡(jiǎn)化為非粘性理想氣體，沖擊波的膨脹假設(shè)為絕熱過(guò)程。則空氣的狀態(tài)方程為

式中，ρ為密度；γ為絕熱指數(shù)；E為內(nèi)能密度。

用LS-Dyna中的MAT_NULL材料和EOS_LINEAR_POLYNOMIAL狀態(tài)方程分別表示空氣的本構(gòu)關(guān)系。空氣材料的模型參數(shù)，見(jiàn)表2所示。

表2 空氣材料的模型參數(shù)

用MAT_HIGH_EXPLOSIVE_BURN材料和EOS_JWL狀態(tài)方程表示TNT炸藥的本構(gòu)關(guān)系。其中JWL狀態(tài)方程模擬炸藥爆轟過(guò)程中壓力、內(nèi)能和比容的關(guān)系，表達(dá)式為

式中，P為材料壓力，V為相對(duì)體積，E為內(nèi)能密度，A、B、R1、R2、ω為材料常數(shù)。

炸藥的爆轟過(guò)程采用“Programmed+beta burn”技術(shù)模擬，每個(gè)炸藥單元的點(diǎn)火時(shí)間由該單元形心至起爆點(diǎn)的距離和爆速確定。TNT炸藥材料及JWL狀態(tài)方程參數(shù)，見(jiàn)表3所示。

表3 TNT炸藥材料及JWL方程參數(shù)

鋼筋混凝土、表層覆土及下面的中風(fēng)化砂巖都采用彈塑性模型。用MAT_PLASTIC_KINEMATIC表示其本構(gòu)關(guān)系。材料物理力學(xué)參數(shù)，見(jiàn)表1所示。

因?yàn)樗淼揽v向較長(zhǎng)，橫向尺寸相對(duì)較小，可以按照平面應(yīng)變問(wèn)題處理，選取隧道圍巖最不利斷面進(jìn)行計(jì)算。

單元建立和網(wǎng)格劃分，巖（土）體用平面實(shí)體單元PLANE2來(lái)模擬；錨桿擋墻（肋柱）面板、隧道襯砌結(jié)構(gòu)和基坑邊墻均用梁?jiǎn)卧狟EAM3來(lái)模擬；錨桿采用桿單元LINK1單元模擬。采用三角形節(jié)點(diǎn)單元進(jìn)行有限元網(wǎng)格劃分。

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

4.2.1 處理區(qū)防爆墻的防爆作用模擬

炸藥起爆后，鋼筋混凝土防爆墻及鄰近巖土介質(zhì)的變形與破壞，見(jiàn)圖4。

圖4 不同時(shí)刻的破壞結(jié)果

可見(jiàn)，用LS-DYNA軟件可以形象地模擬炸藥爆炸時(shí)爆炸波的傳播及介質(zhì)的破壞過(guò)程。數(shù)值模擬表明，在本次模擬情況下，炸藥爆炸時(shí)處理區(qū)的鋼筋混凝土防爆墻被完全炸壞。防爆墻與下穿道之間的覆蓋土被完全破壞，處理區(qū)底部在沖擊波的直接沖擊下形成了較深的炸坑。防爆墻及鄰近巖土介質(zhì)的破壞，會(huì)消耗一部分爆炸能量，因此，處理區(qū)的鋼筋混凝土防爆墻具有一定的耗能保護(hù)作用。

4.2.2 輕軌隧洞的響應(yīng)分析

現(xiàn)選取隧洞迎爆面的墻面中間作為研究對(duì)象，其x向應(yīng)力、y向應(yīng)力、有效應(yīng)力（米塞斯應(yīng)力）、最大剪應(yīng)力、最終位移及振速等數(shù)值模擬結(jié)果，見(jiàn)圖5所示。為節(jié)省篇幅，迎爆面的拱腳及墻角處的數(shù)值模擬結(jié)果略。

圖5 隧道迎爆面墻壁中間處的響應(yīng)

在LS-DYNA中規(guī)定，拉應(yīng)力為正，壓應(yīng)力為負(fù)。從數(shù)值模擬結(jié)果可以看出，隧道迎爆面中間處的應(yīng)力，x向應(yīng)力的最大壓應(yīng)力為0.019MPa，最大拉應(yīng)力為0.015 MPa；y向應(yīng)力最大壓應(yīng)力為0.11MPa，最大拉應(yīng)力為0.22MPa；最大米塞斯應(yīng)力為0.2MPa，最大剪應(yīng)力為0.11MPa。最終位移為0.098cm，x向振速（最大振速）為12cm/s。

可見(jiàn)，本次爆炸情況下，隧道迎爆面墻壁中間處的壓應(yīng)力、拉應(yīng)力、最大剪應(yīng)力及最終位移等都較小，材料處于彈性變形狀態(tài)。最大振速值（12cm/s）也未超出國(guó)家標(biāo)準(zhǔn)《爆破安全規(guī)程》（GB6722-2003）中關(guān)于交通隧道允許振速15cm/s的規(guī)定。因此，在隧道墻壁中間不會(huì)產(chǎn)生破壞。

5 結(jié)論與建議

通過(guò)顯式有限元數(shù)值計(jì)算分析，得到如下結(jié)論：

（1）數(shù)值計(jì)算結(jié)果表明，在充分考慮CNG混合氣體極端爆炸情況下，在鄰近輕軌三號(hào)線隧洞周邊產(chǎn)生的應(yīng)力、位移較小，材料處于彈性變形狀態(tài)。隧洞周邊振動(dòng)速度也小于國(guó)家標(biāo)準(zhǔn)《爆破安全規(guī)程》（GB6722-2003）規(guī)定的交通隧洞的允許振速。因此，本項(xiàng)目CNG加氣站的建設(shè)對(duì)輕軌三號(hào)線隧洞的安全性影響較小，各項(xiàng)指標(biāo)符合相關(guān)規(guī)定。

(2)本次爆炸模擬情況下，防爆墻及鄰近巖土介質(zhì)的破壞，會(huì)消耗一部分爆炸能量，因此，處理區(qū)的鋼筋混凝土防爆墻具有一定的耗能保護(hù)作用。

（3）隧洞迎爆面墻壁中間的最大振速值為12cm/s（允許振速為15cm/s），為提高該隧洞的安全系數(shù)，建議將鋼筋混凝土防爆墻適當(dāng)加厚10～20cm，以進(jìn)一步提高防爆效果。

（4）用LS-DYNA軟件模擬炸藥爆炸時(shí)爆炸波的傳播及介質(zhì)的破壞過(guò)程是可行的。

[1]黃海波，楊建軍，李開(kāi)國(guó)等.CNG加氣站設(shè)備失效與爆炸燃燒風(fēng)險(xiǎn)評(píng)價(jià)[J].西華大學(xué)學(xué)報(bào)（自然科學(xué)版），2005，24(4).

[2]姚東，楊明，楊豪.CNG汽車加氣站火災(zāi)危險(xiǎn)性分析與預(yù)防[J].工業(yè)安全與環(huán)保，2007，33(12):42-44.

[3]譚金會(huì)，何太碧，楊菡，林秀蘭.CNG加氣站設(shè)備安全風(fēng)險(xiǎn)評(píng)價(jià)的關(guān)鍵問(wèn)題[J].天然氣工業(yè)，2008，28(11).

[4]黃鄭華，李建華，張良.壓縮天然氣加氣站的安全設(shè)計(jì)[J].油氣儲(chǔ)運(yùn)，2005，24(6)．