張文泉，劉　毓(1.山東科技大學(xué)礦業(yè)與安全工程學(xué)院，山東青島266590; 2.礦山災(zāi)害預(yù)防與控制國家重點實驗室(培育基地)，山東青島266590; .山東科技大學(xué)地球科學(xué)與工程學(xué)院，山東青島266590)

?

巖溶隧道突涌水危險性的多級模糊綜合評價*

張文泉1，2，劉毓3
(1.山東科技大學(xué)礦業(yè)與安全工程學(xué)院，山東青島266590; 2.礦山災(zāi)害預(yù)防與控制國家重點實驗室(培育基地)，山東青島266590; 3.山東科技大學(xué)地球科學(xué)與工程學(xué)院，山東青島266590)

摘要:為降低巖溶隧道突涌水的危害，減少人員傷亡及經(jīng)濟(jì)損失，采用理論分析與實例驗證相結(jié)合的方法，通過分析各地質(zhì)因素的影響，構(gòu)造以地層因素、地質(zhì)構(gòu)造、地形地貌、水文地質(zhì)為一級評價指標(biāo)，以巖溶化程度、巖層組合等10個因素為二級評價指標(biāo)的評價體系，并依此建立了巖溶隧道突涌水危險性多級模糊綜合評價模型，對齊岳山隧道的2個區(qū)段進(jìn)行了突涌水危險性評價。評價過程中，采用了改進(jìn)的方法確定各影響因素權(quán)重，克服了傳統(tǒng)方法需要一致性檢驗的缺陷，針對各影響因素的不同特點，采用定性與定量相結(jié)合的方式分別確定其隸屬度。將得到的評價結(jié)果與實際施工結(jié)果進(jìn)行比較，結(jié)果表明，理論分析結(jié)果跟現(xiàn)場實際狀況基本一致。

關(guān)鍵詞:巖溶隧道;突涌水;多級模糊綜合評價;評價結(jié)果

0　引言

中國西部山區(qū)巖溶分布廣泛，水文地質(zhì)條件較為復(fù)雜，加之開挖不當(dāng)?shù)热藶橐蛩赜绊?，一系列不良巖溶地質(zhì)災(zāi)害接踵而至，其中巖溶隧道突涌水災(zāi)害尤為突出，造成了大量人員傷亡及經(jīng)濟(jì)損失［1－2］。

施工前對巖溶隧道突涌水進(jìn)行危險性評價，分析災(zāi)害的危險程度，提前做好防治工作，有很重要的現(xiàn)實意義。許振浩［3］等將突水突泥災(zāi)害的評價分為3個階段，并建立起一種新型防災(zāi)機(jī)制。毛邦燕［4］等通過研究其致災(zāi)機(jī)理，對巖溶突涌水進(jìn)行了評價，提出針對性的治理方案?？镄牵?］通過對巖溶地質(zhì)災(zāi)害進(jìn)行深入分析，建立了地質(zhì)災(zāi)害預(yù)警系統(tǒng)。上述研究成果為突涌水災(zāi)害的評價及防治提供了寶貴經(jīng)驗，但其在運用層析分析法計算權(quán)重時，難以保證判斷矩陣的一致性，文中對此進(jìn)行改進(jìn)，避免了一致性檢驗需要的大量計算，并將其運用于巖溶隧道突涌水的多級模糊評價中。

1　評價指標(biāo)體系的構(gòu)建

巖溶隧道突涌水受多個方面因素影響，如隧道周圍氣候條件、地層條件，外界施工影響等。這里僅考慮地質(zhì)因素方面，地質(zhì)因素主要包括地層因素、地質(zhì)構(gòu)造、地形地貌及水文地質(zhì)等主要因素，又可繼續(xù)將它們細(xì)分為巖溶化程度、巖層組合條件等指標(biāo)。

1．1地層因素

巖溶化程度的大小影響導(dǎo)水構(gòu)造的形態(tài)及分布特征，對巖溶隧道突涌水的影響主要體現(xiàn)在削弱巖層強(qiáng)度方面。巖層組合條件是指可溶性巖層與非可溶性巖層組合形態(tài)，由于性質(zhì)不同，它們之間不同的組合形態(tài)呈現(xiàn)出的實際狀態(tài)也有所不同。在可溶性巖層與非可溶巖層接觸帶，由于地下水的相互流通，對可溶巖層的侵蝕得以加重。巖層產(chǎn)狀能夠改變地下水的滲流特征，影響地下水與可溶性巖石的化學(xué)反應(yīng)程度。向斜或背斜翼部巖層傾角為25°～65°時，對巖溶的發(fā)育促進(jìn)作用最強(qiáng)［6］。

1．2地質(zhì)構(gòu)造

不良地質(zhì)構(gòu)造按其對突涌水的影響可分為多個方面，如斷裂破碎帶的寬度、斷裂影響帶的寬度及褶皺的構(gòu)造形態(tài)等。斷裂破碎帶的存在改變了地層中原有的滲透系數(shù)等參數(shù)，為地下水的滲流提供了一種良好的途徑，斷層及褶皺的富水性及導(dǎo)水性等特征的不同決定了對突涌水的貢獻(xiàn)不同。

1．3地形地貌

地表巖溶和地下巖溶相互連通，相互影響，地表巖溶形態(tài)(溶溝、溶蝕洼地、漏斗等)是地上外界環(huán)境與地下巖溶系統(tǒng)的媒介。地表巖溶的匯水面積大小決定著對地下巖溶的補(bǔ)給水量的多少，從而在一定程度上控制地下巖溶的發(fā)育程度。

1．4水文地質(zhì)

水文地質(zhì)條件的優(yōu)劣影響著巖溶隧道突涌水的危險性程度，主要通過地下水位的高低及地下水水壓力的大小等來表現(xiàn)。地下水位的選取一般以地下水水面到隧道底板的高程為準(zhǔn)，隨著高程的降低，地下水的溶蝕性有所降低，危害性也逐漸減弱。地下水水壓力是突涌水的動力因素，隨著水壓力的增大，其對隔水層的影響加深，超出隔水層的承壓界限后，突涌水就會發(fā)生。

根據(jù)模糊綜合評價原理，可將這些因素分為3個層次，具體結(jié)構(gòu)及從屬關(guān)系如圖1所示。

圖1　巖溶隧道突涌水危險性評價體系Fig．1 Risk assessment system of water inrush in karst tunnels

2　突涌水多級模糊綜合評價

多級模糊綜合評價作為一種多因素決策技術(shù)，在實際工程中被廣泛應(yīng)用。它結(jié)合了層次分析法及模糊綜合評判法的優(yōu)點，引進(jìn)隸屬度這一模糊中介，通過建立隸屬函數(shù)，與權(quán)重集合進(jìn)行相應(yīng)的矩陣運算依照最大隸屬度原則最終獲得評價結(jié)果［7］。

巖溶隧道突涌水影響因素眾多，并且每一個主要因素還由若干個子因素決定，多級模糊綜合評價系統(tǒng)性強(qiáng)、評價全面，可以很好地綜合各因素的作用。

2. 1確定因素集

針對每級評價對象的不同，參照圖1確立因素集U = { U1，U2，U3，U4} = {地層因素，地質(zhì)構(gòu)造，地形地貌，水文地質(zhì)} ; U1= { u1，u2，u3} = {巖溶化程度，巖層組合，巖層傾角} ; U2= { u4，u5，u6} ={斷裂破碎帶寬度，斷裂影響帶寬度，褶皺構(gòu)造形態(tài)} ; U3= { u7，u8} = {地表巖溶形態(tài)，地表匯水面積} ; U4= { u9，u10} = {地下水位，地下水水壓力}。

2. 2確定并求解權(quán)向量

根據(jù)層次分析法的9標(biāo)度法確定初始判斷矩陣A = (aij)n×n，(i，j =1，2，3…，n)，aij為指標(biāo)Ai對Aj的標(biāo)度。

傳統(tǒng)層次分析法求解權(quán)向量的方法是計算判斷矩陣A各行的乘積的n次方根，將得到的向量正規(guī)化，得到權(quán)向量。由于得到的權(quán)向量是否具有一致性不能確定，需繼續(xù)對其進(jìn)行一致性檢驗，計算量巨大。

設(shè)矩陣B = (bij)n×n，(i，j =1，2，3，…，n)，若對于任意i，j，k(1≤k≤n)，滿足bii=0．5，0＜bij＜1，bij+ bji=1，bij= bik－bjk+0．5，則稱B為模糊一致性矩陣，其滿足一致性要求。可依此對A進(jìn)行修正，具體步驟如下:運用轉(zhuǎn)換公式

將A轉(zhuǎn)換成矩陣A' = (yij(x) )n×n，(i，j =1，2，3…n)，yii(x) =0．5，yij(x) + yji(x) =1，為滿足0＜yij(x)＜1，得x≥18．隨著x取值增大，轉(zhuǎn)換后各標(biāo)度范圍逐漸變小，對各因素權(quán)重表述不夠清晰，故x取值18．

再通過求矩陣A'的行和，即

然后利用轉(zhuǎn)換公式

將矩陣A'轉(zhuǎn)化為A″= (aij)n×n，(i，j = 1，2，3…n)，矩陣模糊一致性矩陣［8］。

最后，通過求約束規(guī)劃問題的解［9］來確定權(quán)向量

求得由于權(quán)重為非負(fù)值，(6)式中a的取值范圍需要滿足a≥αij，隨a取值增大，各影響因素的權(quán)重之差逐漸縮小，在實際工程運用中不利于區(qū)分各指標(biāo)的影響程度，故取a = 1．

上述方法通過將初始判斷矩陣轉(zhuǎn)化成模糊一致性矩陣，改進(jìn)了傳統(tǒng)AHP法需要一致性檢驗的缺陷，減少了計算量，提供了一種求解權(quán)向量的新思路。

2. 3評價集的構(gòu)建

通過總結(jié)過去巖溶隧道突涌水危險性評價的成果，結(jié)合實際案例的現(xiàn)場狀況，綜合分析各致災(zāi)因素的影響作用，根據(jù)災(zāi)害可能造成的后果建立評價集: V = { V1，V2，V3，V4，V5} = {高風(fēng)險，較高風(fēng)險，中等風(fēng)險，較低風(fēng)險，低風(fēng)險}。

2. 4隸屬度的確定

確定評價指標(biāo)的隸屬度時，根據(jù)影響突涌水的方式，將其進(jìn)行定性或定量分析計算。對于定量指標(biāo)，根據(jù)模糊數(shù)學(xué)理論，通過構(gòu)建隸屬度函數(shù)來表達(dá)，定性指標(biāo)無法用具體數(shù)值進(jìn)行量化，但是其具有明顯的劃分特征。

根據(jù)突涌水定量指標(biāo)的性質(zhì)，隸屬度確定采用梯形隸屬度函數(shù)

定性指標(biāo)方面，可按其自身特征及影響程度進(jìn)行分級，再運用模糊統(tǒng)計的方法確定其隸屬度，對于構(gòu)造隸屬度函數(shù)較為困難的定量指標(biāo)，將其歸入定性指標(biāo)來進(jìn)行探討。各因素具體取值見表1．

2. 5多級模糊綜合評價過程

評價的過程是由評價指標(biāo)的最低層次到最高層次依次逐級進(jìn)行評價，綜合考慮各因素之間的從屬關(guān)系，最終得出總的評價結(jié)果。具體過程以評價體系自上而下為A，B，C 3層進(jìn)行介紹:

首先，求得A層權(quán)向量A = = (a1，a2，…，am)，B層權(quán)向量Ai= (ai1，ai2，…，ain)，i = (1，2，…，m)．

表1　定性指標(biāo)的隸屬度取值Tab．1 Membership of qualitative factors

表2　單因素隸屬度評價值Tab．2 Evaluation value of each factor’s membership

其次，分析計算各評價指標(biāo)的隸屬度，將所得到的結(jié)果對比表2得出C層評價矩陣

式中i = (1，2，…，m)．

對B層進(jìn)行一級模糊評價，得

式中i = (1，2，…，m)．

對B層的模糊評價的結(jié)果可以作為對A層模糊評價時的評價矩陣，即R = Bi．對A層進(jìn)行二級模糊評價

根據(jù)最大隸屬度原則，bt= max{ bk}，t項所處等級即為突涌水的危險性級別。

3　實例驗證

宜萬鐵路齊岳山隧道全長10 528 m，最大埋深670 m．隧道地處巖溶地區(qū)，附近地層多由可溶性巖石構(gòu)成，伴有多條斷層，沿斷層溶溝、巖溶洼地等發(fā)育，巖溶管道沿斷層張性裂隙溝通地表巖溶及地下巖溶。選取隧道中巖溶較發(fā)育的2個區(qū)段對其進(jìn)行評價［10－11］。

3. 1研究區(qū)域概況

第1區(qū)段(隧道里程DK362 + 700—DK364 + 380，長度1 680 m) :位于齊岳山巖溶壟脊槽谷，地表巖溶發(fā)育，巖溶形態(tài)各異且分布廣泛，地上巖溶與地下巖溶連通較好，深部巖溶管道、溶隙發(fā)育良好，存在地下暗河系統(tǒng)。地下巖溶發(fā)育到最深350 m處，巖溶系統(tǒng)中含水層具有較強(qiáng)的富水性。隧道在該區(qū)段位于地下水深部循環(huán)帶，地下水補(bǔ)給徑流區(qū)，線路附近有F5，F(xiàn)6，F(xiàn)7，F(xiàn)84條斷層。

第2區(qū)段(隧道里程DK361 + 257—DK362 + 700，長度1 443 m) :位于齊岳山背斜斜坡地帶，地表巖溶洼地、溶溝分布較多，地下巖溶主要包括大魚泉、小魚泉地下河系統(tǒng)，多處發(fā)育地下溶洞，巖層接觸帶巖溶管道發(fā)育良好，巖溶系統(tǒng)中含水層具有較強(qiáng)的富水性。此區(qū)段屬地下水匯集排泄區(qū)，接近水平徑流帶，隧道進(jìn)口略高于大魚泉地下河排泄口，該段線路經(jīng)過F1，F(xiàn)2，F(xiàn)3，F(xiàn)4斷層。

3. 2安全性評價

通過分析地質(zhì)資料，顯示該區(qū)段地層因素和地質(zhì)構(gòu)造影響程度較為接近，且明顯強(qiáng)于地形地貌及水文地質(zhì)的影響，水文地質(zhì)條件對突涌水的貢獻(xiàn)略強(qiáng)于地形地貌。其中，地層因素方面各指標(biāo)影響程度巖溶化程度＞巖層組合＞巖層傾角;地質(zhì)構(gòu)造方面斷裂破碎帶寬度＞褶皺構(gòu)造形態(tài)＞斷層影響帶寬度;地形地貌方面地表巖溶形態(tài)＞地表匯水面積;水文地質(zhì)方面地下水位＞地下水壓力。

基于上述結(jié)論，根據(jù)表1運用AHP的9標(biāo)度法確定一級指標(biāo)的初始判斷矩陣

運用公式(1)對其進(jìn)行轉(zhuǎn)化，得

再運用公式(2)，(3)繼續(xù)進(jìn)行轉(zhuǎn)化，得到模糊一致性矩陣

最后，依據(jù)公式(6)進(jìn)行求解的一級指標(biāo)權(quán)重值(0. 326，0. 270，0. 186，0. 218)，依次進(jìn)行二級權(quán)重值的求解，具體數(shù)值見表3．

表3　各評價指標(biāo)權(quán)重Tab．3 Weight of each factor

將上述所得權(quán)重值與工程資料進(jìn)行對比分析，發(fā)現(xiàn)其符合該區(qū)段各評價指標(biāo)的影響程度。

從地質(zhì)資料中獲取各指標(biāo)的相應(yīng)取值，根據(jù)2. 4節(jié)提到的方法，分析及計算得出各指標(biāo)的隸屬度，對照表2得出第1段區(qū)域C層評價矩陣

對B層進(jìn)行一級模糊評價，得

得到對A層進(jìn)行二級模糊評價單因素評價矩陣

對A層進(jìn)行二級模糊評價，得第1區(qū)段綜合評價值

B = (0．203 0．291 0．242 0．166 0．098)

對第2區(qū)段各項指標(biāo)進(jìn)行同樣的運算，得第2區(qū)段多級模糊綜合評價值

B = (0．255 0．209 0．209 0．162 0．150)

據(jù)最大隸屬度原則得知，第1區(qū)段評價結(jié)果為突涌水高危險區(qū)，第2區(qū)段評價為突涌水極危險區(qū)。

3. 3結(jié)果比較

在實際施工過程中，所選區(qū)段多次發(fā)生突涌水事故。其中，第1區(qū)段共發(fā)生突涌水25次，最大突水量超過2 500 m3/h;第2區(qū)段共發(fā)生突涌水23次，最大突水量1 300 m3/h．

將評價結(jié)果和實際施工狀況進(jìn)行對比，發(fā)現(xiàn)兩者的結(jié)果基本一致，說明該評價模型能夠較為準(zhǔn)確地表達(dá)出實際狀態(tài)下突涌水的危險程度。

4　結(jié)論

文中將多級模糊綜合評價方法運用到巖溶隧道突涌水安全性評價中，通過分析巖溶隧道突涌水影響因素，構(gòu)建了巖溶隧道突涌水危險性多級模糊評價模型，并結(jié)合實例進(jìn)行了驗證。得到如下結(jié)論

1)多級模糊綜合評價能夠充分考慮突涌水各個因素的影響，綜合考慮各影響因素之間的關(guān)系，對呈現(xiàn)模糊性的實際問題做出較為科學(xué)合理的量化;

2)將評價模型運用于工程實例，得到的結(jié)果驗證了評價模型的可行性，為巖溶隧道突涌水災(zāi)害預(yù)測及防治提供了一種有效途徑;

3)運用多級模糊綜合評價解決實際問題時，計算較為復(fù)雜，進(jìn)行權(quán)向量求解時對初始判斷矩陣的確定主觀性較強(qiáng)，隨著應(yīng)用的增多這方面有待改進(jìn)。

參考文獻(xiàn)References

［1］許振浩，李術(shù)才，李利平，等．基于風(fēng)險動態(tài)評估與控制的巖溶隧道施工許可機(jī)制［J］．巖土工程學(xué)報，2011，33(11) : 1 715－1 725. XU Zhen-hao，LI Shu-cai，LI Li-ping，et al．Construction permit mechanism of karst tunnels based on dynamic assessment and management of risk［J］．Chinese Journal of Geotechnical Engineering，2011，33 (11) : 1 715－1 725.

［2］李利平，李術(shù)才，陳軍，等．基于巖溶突涌水風(fēng)險評價的隧道施工許可機(jī)制及其應(yīng)用研究［J］．巖石力學(xué)與工程學(xué)報，2011，30(7) : 1 345－1 354. LI Li-ping，LI Shu-cai，CHEN Jun，et al．Construction license mechanism and its application based on karst water inrush risk evaluation［J］．Chinese Journal of Rock Mechanics and Engineering，2011，30 (7) : 1 345－1 354.

［3］許振浩，李術(shù)才，李利平，等．基于層次分析法的巖溶隧道突水突泥風(fēng)險評估［J］．巖土力學(xué)，2011，32(6) : 1 757－1 765. XU Zhen-hao，LI Shu-cai，LI Li-ping，et al．Risk assessment of water or mud inrush of karst tunnels based on analytic hierarchy process［J］．Rock and soil Mechanics ，2011，32(6) : 1 757－1 765.

［4］毛邦燕，許模，蔣良文．隧道巖溶突水、突泥危險性評價初探［J］．中國巖溶，2010，29(2) : 183－189. MAO Bang-yan，XU Mo，JIANG Liang-wen．Preliminary study of risk assessment of water and mud inrush in karst tunnels［J］．Carsologica Sinica，2010，29(2) : 183 －189.

［5］匡星，白明洲，王成亮，等．基于模糊評價方法的隧道巖溶突水地質(zhì)災(zāi)害綜合預(yù)警方法［J］．公路交通科技，2010，27(11) : 100－103. KUANG Xing，BAI Ming-zhou，WANG Cheng-liang，et al．Research of comprehensive warning of water inrush hazards in karst tunnel based on fuzzy evolution method ［J］．Journal of Highway and Transportation Research and Development，2010，27(11) : 100－103.

［6］周宗青，李術(shù)才，李利平，等．巖溶隧道突涌水危險性評價的屬性識別模型及其工程應(yīng)用［J］．巖土力學(xué)，2013，34(3) : 818－826. ZHOU Zong-qing，LI Shu-cai，LI Li-ping，et al．Attribute recognition model of fatalness assessment of water inrush in karst tunnels and its application［J］．Rock and Soil Mechanics，2013，34(3) : 818－826.

［7］謝季堅，劉承平．模糊數(shù)學(xué)方法及其應(yīng)用［M］．武漢:華中科技大學(xué)出版社，2006. XIE Ji-jian，LIU Cheng-ping．Fuzzy mathematics and its applications［M］．Wuhan: Huazhong University of Science and Technology Press，2006.

［8］李永，胡向紅，喬箭．改進(jìn)的模糊層次分析法［J］．西北大學(xué)學(xué)報:自然科學(xué)版，2005，35(1) : 11－12. LI Yong，HU Xiang-hong，QIAO Jian．Improved fuzzy AHP［J］．Northwest University Journal: Nature Science，2005，35(1) : 11－12.

［9］呂躍進(jìn)．基于模糊一致矩陣的模糊層次分析法的排序［J］．模糊數(shù)學(xué)與系統(tǒng)，2002，16(2) : 79－85. LV Yue-jin．Sort of fuzzy AHP based on fuzzy consistent matrix［J］．Fuzzy Systems and Marhematics，2002，16 (2) : 79－85.

［10］王成亮．巖溶區(qū)隧道施工地質(zhì)災(zāi)害風(fēng)險評價方法研究［D］．北京:北京交通大學(xué)，2010．WANG Cheng-liang．Study on risk assessment method of geological disaster during tunnel construction in karst area［D］．Beijing: Beijing Jiaotong University，2010.

［11］沈詳明，劉坡拉，汪繼鋒．基于層次分析法的鐵路巖溶隧道突水風(fēng)險評價［J］．鐵道工程學(xué)報，2010(12) : 56－63. SHEN Xiang-ming，LIU Po-la，WANG Ji-feng．Evaluation of water inrush risks of karst tunnel with analysis hierarchy process［J］．Journal of Railway Engineering Society，2010(12) : 56－63.

Multilevel fuzzy comprehensive evaluation of water inrush in karst tunnels

ZHANG Wen-quan1，2，LIU Yu3
(1. College of Mining and Safety Engineering，Shandong University of Science and Technology，Qingdao 266590，China; 2. National Key Laboratory of Mine Disaster Prevention and Control(Cultivation Base)，Qingdao 266590，China; 3. College of Earth Science and Engineering，Shandong University of Science and Technology，Qingdao 266590，China)

Abstract:In order to reduce the damage of water inrush in karst tunnels and control the casualties and economic losses，this article built a multilevel fuzzy comprehensive evaluation model to evaluate two sections of the Qiyueshan Tunnel by combining theoretical analysis and verification method．This model was built according to a risk assessment system which had two levels．Each level was made through analysis of each geological factor’s influence．First level has stratum factors，geological structures，topography and hydrogeological．Second level has karst degree，rock combination and other eight factors．In the process of evaluation，when it comes to determine the weight of each factor，this article used an improved method to overcome the traditional method’s defect which requiring consistency-test．As the difference between each factor，this article also used a combination of qualitative and quantitative manner to determine their memberships．At last，by compared with each other，the evaluation result showed a satisfying conclusion compared with real condition．It’s feasible to use multilevel fuzzy comprehensive theory to evaluate water inrush in karst tunnels．

Key words:karst tunnel; water inrush; multilevel fuzzy comprehensive evaluation; evaluation results

通訊作者:張文泉(1965－)，男，山東濰坊人，教授，E-mail: wenquanzhang@163．com

*收稿日期:2015－10－17責(zé)任編輯:劉潔

DOI:10．13800/j．cnki．xakjdxxb．2016．0207

文章編號:1672－9315(2016) 02－0187－06

中圖分類號:U 25

文獻(xiàn)標(biāo)志碼:A