韓永勝,趙文峰,高國(guó)紅,駱 勇,4
(1.云南召瀘高速公路有限責(zé)任公司,云南 曲靖 655099;2.四川省煤礦瓦斯(煤層氣)工程研究中心,四川 成都 610072;3.云南交投集團(tuán)投資有限公司,云南 昆明 650228;4.四川省煤田地質(zhì)工程勘察設(shè)計(jì)研究院,四川 成都 610072)
遇水搭橋,逢山鑿隧,天塹變通途。隨著交通建設(shè)的快速發(fā)展,隧道工程日益增多[1~5]。在公路、鐵路隧道選線時(shí),雖然盡可能避免穿越含煤地層或煤礦采空區(qū),但往往改道建設(shè)會(huì)增加成倍的投入,不得不采用隧道穿越,相應(yīng)地增加了工程建設(shè)的風(fēng)險(xiǎn)[6~8]。目前,國(guó)內(nèi)外學(xué)者在采空區(qū)穩(wěn)定性評(píng)價(jià)方面進(jìn)行了大量的研究工作,而針對(duì)煤礦采空區(qū)瓦斯危害性評(píng)價(jià)與治理方法少有關(guān)注[5,9~12]。為有效防治公路隧道因穿越采空區(qū)而導(dǎo)致瓦斯事故的發(fā)生,建立隧道采空區(qū)瓦斯危害性評(píng)價(jià)體系和方法,評(píng)價(jià)其危害程度,對(duì)隧道的安全、經(jīng)濟(jì)、高效施工具有現(xiàn)實(shí)意義。
根據(jù)公路隧道施工特點(diǎn),參考以往鐵路、公路瓦斯隧道以及煤礦采空區(qū)風(fēng)險(xiǎn)評(píng)價(jià)分析的資料及成果[7~15],總結(jié)歸納出影響公路隧道采空區(qū)瓦斯危害性的3個(gè)一級(jí)指標(biāo),共10個(gè)二級(jí)評(píng)價(jià)指標(biāo),以此建立危險(xiǎn)性評(píng)價(jià)指標(biāo)體系(圖1)。
圖1 公路隧道采空區(qū)瓦斯危害性評(píng)價(jià)指標(biāo)體系
在采用層次分析法(Analytic Hierarchy Process,AHP)和關(guān)系矩陣法分別計(jì)算權(quán)重系數(shù)的基礎(chǔ)上,建立組合權(quán)重模型最終確定評(píng)價(jià)指標(biāo)權(quán)重系數(shù)。
層次分析法是一種廣泛應(yīng)用于不同領(lǐng)域的一種較好的權(quán)重確定方法[16]。為了能定量確定各指標(biāo)的權(quán)重,首先需要建立同層次指標(biāo)間的判別矩陣。對(duì)于同層次指標(biāo)的兩兩比較,重要性標(biāo)度按照表1取值[17]。利用數(shù)學(xué)方法計(jì)算出判別矩陣的最大特征根λmax及其對(duì)應(yīng)的特征向量W,得到各指標(biāo)的權(quán)重系數(shù)。采用隨機(jī)一致性比率CR來(lái)判別矩陣的一致性。若CR<0.1,則認(rèn)為判別矩陣具有可接受的不一致性;否則,重新賦值和修正計(jì)算[18]。
表1 標(biāo)度方法
(1)
式中:n為矩陣的階數(shù);RI為同階平均隨機(jī)一致性指標(biāo),其值為0.00(1,2階矩陣)、0.52(3階矩陣)、0.89(4階矩陣)和1.12(5階矩陣)。
在構(gòu)建一級(jí)評(píng)價(jià)指標(biāo)和二級(jí)評(píng)價(jià)指標(biāo)判斷矩陣的基礎(chǔ)上,計(jì)算各評(píng)價(jià)指標(biāo)在該級(jí)中的權(quán)重系數(shù),并利用CR值進(jìn)行一致性檢驗(yàn)。公路隧道采空區(qū)瓦斯危害性一、二級(jí)評(píng)價(jià)指標(biāo)判斷矩陣、權(quán)重值及CR值計(jì)算結(jié)果見表2~6。
表2 一級(jí)評(píng)價(jià)指標(biāo)判斷矩陣、權(quán)重值及CR值
表3 瓦斯地質(zhì)因素判斷矩陣、權(quán)重值及CR值
表4 工程因素判斷矩陣、權(quán)重值及CR值
表5 人為因素判斷矩陣、權(quán)重值及CR值
根據(jù)關(guān)系矩陣相互作用原理,主對(duì)角線元素表示系統(tǒng)的主要影響因素,主對(duì)角線以外元素表示因素之間的相互作用關(guān)系。主對(duì)角線上元素取值為空,表示各因素不影響其自身,只能通過(guò)與其他因素相互作用來(lái)影響系統(tǒng);主對(duì)角線以外元素取值一般采用專家半定量取值法[19],根據(jù)相互作用程度取為0,1,2,3,4(表7)。圖2表示N個(gè)參數(shù)相互作用的關(guān)系矩陣,Iij表示主要因素Pi對(duì)主要因素Pj的影響,Ci為第i行非主對(duì)角線元素值之和,Ei為第i列非主對(duì)角線元素值之和[6,20]。
表6 各層次評(píng)價(jià)指標(biāo)權(quán)重
表7 關(guān)系矩陣的專家半定量取值法
圖2 N個(gè)參數(shù)相互作用關(guān)系矩陣
通常計(jì)算所有參與評(píng)價(jià)指標(biāo)的活動(dòng)性指數(shù)ki:
(2)
該指數(shù)表示每一指標(biāo)的因果值總和占系統(tǒng)總因果值的百分比,即各指標(biāo)的權(quán)重。
根據(jù)上述方法,構(gòu)建公路隧道采空區(qū)瓦斯危害性評(píng)價(jià)關(guān)系矩陣編碼表,并計(jì)算ki值(表8)。關(guān)系矩陣編碼表中主對(duì)角線元素Pi(i=1,2,…,10)表示采空區(qū)瓦斯災(zāi)害危險(xiǎn)性評(píng)價(jià)指標(biāo):P1為瓦斯?jié)舛?;P2為瓦斯資源量;P3為連通及封閉性能;P4為采空區(qū)積水情況;P5為埋深;P6為采空區(qū)與隧道關(guān)系;P7為開挖斷面大小;P8為瓦斯管理;P9為瓦斯隧道施工經(jīng)驗(yàn);P10為作業(yè)人員素質(zhì)。
表8 評(píng)價(jià)指標(biāo)相互作用關(guān)系矩陣編碼及活動(dòng)性指數(shù)
綜合考慮層次分析法和關(guān)系矩陣法所確定的權(quán)重系數(shù),將組合權(quán)重模型定義為:
Wi=tpi+(1-t)ki
(3)
式中:pi,ki分別為層次分析法、關(guān)系矩陣法確定的權(quán)重系數(shù);t為系數(shù),取值為0.5,i=1,2,…,10。通過(guò)組合權(quán)重計(jì)算得到二級(jí)指標(biāo)權(quán)重系數(shù)(表9)。
表9 采空區(qū)瓦斯危害性評(píng)價(jià)指標(biāo)權(quán)重系數(shù)
綜合隧道工程、煤礦采空區(qū)危險(xiǎn)性評(píng)價(jià)的成熟經(jīng)驗(yàn),公路隧道采空區(qū)瓦斯危害性劃分為四個(gè)等級(jí)(表10),即較低(Ⅰ)、一般(Ⅱ)、較高(Ⅲ)、極高(Ⅳ)。劃分等級(jí)時(shí)遵循采空區(qū)瓦斯危害性越嚴(yán)重則量化值越大的原則。
表10 采空區(qū)瓦斯危害性評(píng)價(jià)等級(jí)標(biāo)準(zhǔn)
根據(jù)二級(jí)評(píng)價(jià)指標(biāo)危害性程度評(píng)分量化值和確定的權(quán)重系數(shù),按疊加法求得其綜合評(píng)價(jià)值y,對(duì)評(píng)價(jià)樣本進(jìn)行綜合評(píng)判。公式如下:
y=∑αixi
(4)
式中:xi為i類評(píng)價(jià)指標(biāo);αi為xi的權(quán)重。
令物元R=(事物,特征,量值)=(N,C,V)。事物有n個(gè)特征,即C={c1,c2,…,cn};也有相應(yīng)的量值,即V={v1,v2,…,vn}??赏貙W(xué)評(píng)價(jià)模型按如下步驟構(gòu)建[21]。
(1)確定經(jīng)典域、節(jié)域
設(shè)標(biāo)準(zhǔn)事物為N0 j,關(guān)于特征Ci的量值范圍為Vi j=[ai j,bi j](i=1,2,…,n;j=1,2,…,m),則經(jīng)典域R0 j和節(jié)域物元Rp分別為:
(5)
(6)
式中:Ci(i=1,2,…,n)為決定N0 j的因素;V0ij為N0 j關(guān)于特征Ci所確定的量值范圍;Vpi為Np關(guān)于特征Ci在P條件下所取的量值范圍。
(2)確定待評(píng)物元
對(duì)于待評(píng)采空區(qū)N,把所測(cè)得的數(shù)據(jù)或分析結(jié)果用物元表示,稱為采空區(qū)待評(píng)物元。
(7)
式中:Vi為關(guān)于Ci的量值,即待評(píng)采空區(qū)所得的具體數(shù)值。
(3)計(jì)算待評(píng)物元關(guān)于各級(jí)的關(guān)聯(lián)度
(8)
對(duì)每個(gè)特征Ci,取αi為權(quán)重。
(4)確定評(píng)價(jià)等級(jí)
云南省召夸至瀘西高速公路擬建隧道一座,即白兆隧道。該隧道為分離式隧道,右幅起點(diǎn)樁號(hào)K8+370,止點(diǎn)樁K9+633,右幅長(zhǎng)1263 m,左幅起點(diǎn)樁號(hào)ZK8+400,止點(diǎn)樁號(hào)ZK9+650,左幅長(zhǎng)1250 m。隧道開挖斷面150 m2,圍巖以玄武巖、泥頁(yè)巖、泥質(zhì)粉砂巖、砂巖為主,圍巖分級(jí)為Ⅳ2~Ⅴ2級(jí),復(fù)合式襯砌。據(jù)隧址區(qū)工程勘察成果,會(huì)遇到C10,C14,C17,C22共4層大于0.70 m的煤層,瓦斯含量2.55~3.60 m3/t。里程樁號(hào)左幅ZK9+005—ZK9+098 m、右幅K9+005—K9+090 m 為高瓦斯地段。預(yù)測(cè)了5個(gè)煤礦采空區(qū),編號(hào)為1~5號(hào)。其中,對(duì)隧道施工有直接影響的有3~5號(hào)采空區(qū)(圖3)。本文取其中3~5號(hào)采空區(qū)作為待評(píng)樣本。
圖3 百兆隧道采空區(qū)分布示意
(1)3號(hào)采空區(qū)
根據(jù)鉆孔采樣測(cè)試成果,3號(hào)采空區(qū)瓦斯?jié)舛葹?.05%。預(yù)測(cè)瓦斯資源量為778 m3,占采空區(qū)體積的5%。采空引起的裂隙已連通至風(fēng)氧化帶或地表,地表斷層F919切割采空區(qū),封閉性能差。采空區(qū)處于無(wú)水或少水狀態(tài)。采空區(qū)埋深30~60 m,該采空區(qū)與隧道部分重疊,并橫穿隧道左、右幅。
(2)4號(hào)采空區(qū)
根據(jù)鉆孔采樣測(cè)試成果,4號(hào)采空區(qū)瓦斯?jié)舛葹?.26%。預(yù)測(cè)瓦斯資源量為1249 m3,占采空區(qū)體積的8%。采空引起的裂隙已連通至風(fēng)氧化帶或地表,封閉性能差。采空區(qū)處于無(wú)水或少水狀態(tài)。采空區(qū)埋深40~60 m,位于隧道軸線及底板下方15 m,并橫穿隧道左、右幅。
(3)5號(hào)采空區(qū)
根據(jù)鉆孔采樣測(cè)試成果,5號(hào)采空區(qū)瓦斯?jié)舛葹?.13%,預(yù)測(cè)瓦斯資源量為485 m3,占采空區(qū)體積的4%。采空引起的裂隙已連通至地表以下10 m左右,封閉性能差。采空區(qū)處于無(wú)水或少水狀態(tài)。采空區(qū)埋深63~85 m,其頂板位于隧道軸線底板下方5~17 m,即隧道位于采空區(qū)垮落帶和裂隙帶內(nèi)。
3~5號(hào)采空區(qū)瓦斯地質(zhì)因素與工程因素相關(guān)指標(biāo)有以上差異,其余指標(biāo)相同。如隧道采用雙向六車道高速公路標(biāo)準(zhǔn)建設(shè),設(shè)計(jì)速度100 km/h,路基寬度33.5 m,開挖斷面150 m2。施工隊(duì)伍瓦斯管理制度完善,符合煤礦安全規(guī)程的要求,但有少數(shù)次要項(xiàng)目不落實(shí)。有瓦斯隧道施工經(jīng)驗(yàn),作業(yè)人員全部培訓(xùn),責(zé)任心較強(qiáng),素質(zhì)較好。
按照采空區(qū)瓦斯危害性等級(jí)劃分標(biāo)準(zhǔn)和百兆隧道3~5號(hào)采空區(qū)實(shí)際參數(shù),采用分段線性函數(shù)的方法確定評(píng)分值,二級(jí)評(píng)價(jià)指標(biāo)評(píng)分取值見表11。
表11 二級(jí)評(píng)價(jià)指標(biāo)評(píng)分取值
基于綜合評(píng)判法模型的和方法,可以疊加求得各采空區(qū)綜合評(píng)分值(表12)。
根據(jù)采空區(qū)危害性等級(jí)劃分標(biāo)準(zhǔn),確定經(jīng)典域?yàn)椋?/p>
節(jié)域?yàn)椋?/p>
確定待評(píng)價(jià)物元模型及指標(biāo)如下:
以3號(hào)采空區(qū)關(guān)于“較低(Ⅰ)”等級(jí)的關(guān)聯(lián)度計(jì)算為例。
首先,計(jì)算待評(píng)物元各指標(biāo)關(guān)于“較低(Ⅰ)”等級(jí)的關(guān)聯(lián)度:
K1(v1)=-ρ(0.05,〈0,0.5〉)-1=-0.95
K1(v2)=-ρ(0.05,〈0,0.1〉)-1=-0.95
K1(v3)=-ρ(8,〈0,25〉)-1=7.00
(1)次梁桁架的上弦桿長(zhǎng)度方向采用鋼管搭接,搭接長(zhǎng)度不小于500mm,每個(gè)連接點(diǎn)不少于4個(gè)直角扣件連接。
=-0.47
K1(v5)=-ρ(50,〈0,100〉)-1=49
=-0.27
=-0.47
=-0.14
=-0.14
其次,計(jì)算待評(píng)物元關(guān)于“較低(Ⅰ)”等級(jí)的關(guān)聯(lián)度:
=0.192×(-0.95)+0.163×(-0.95)+0.163×7.00+0.055×(-0.47)+0.076×
49.00+0.086×(-0.27)+0.056×
(-0.47)+0.097×(-0.14)+0.066×
(-0.14)+0.046×(-0.14)=4.42
同理,可計(jì)算3號(hào)采空區(qū)關(guān)于一般(Ⅱ)、較高(Ⅲ)、極高(Ⅳ)等級(jí)及4,5號(hào)采空區(qū)各級(jí)的關(guān)聯(lián)度,計(jì)算結(jié)果見表13,以此確定3~5號(hào)采空區(qū)均為“較低”級(jí)別。
表13 各級(jí)關(guān)聯(lián)度及可拓評(píng)價(jià)結(jié)果
(1)基于組合賦權(quán)模型確定了公路隧道采空區(qū)瓦斯危害性指標(biāo)權(quán)重系數(shù),運(yùn)用綜合評(píng)判法和可拓學(xué)評(píng)判法建立了公路隧道采空區(qū)瓦斯危害性評(píng)價(jià)體系。
(2)評(píng)價(jià)體系成功運(yùn)用到百兆隧道采空區(qū)危害性評(píng)價(jià)中,兩種評(píng)價(jià)方法的分析評(píng)價(jià)結(jié)果一致,且與現(xiàn)場(chǎng)實(shí)際情況相吻合。由此可見,這兩種評(píng)價(jià)方法對(duì)隧道采空區(qū)瓦斯災(zāi)害危險(xiǎn)性進(jìn)行評(píng)價(jià)是可行、可信的。
(3)根據(jù)采空區(qū)危害性評(píng)價(jià)結(jié)果,擬穿越的3個(gè)采空區(qū)危害性程度較低,采取優(yōu)化揭露采空區(qū)通風(fēng)方案,保證了隧道施工順利穿越采空區(qū)。