王述紅， 朱寶強(qiáng)， 張 澤

(東北大學(xué) 資源與土木工程學(xué)院， 遼寧 沈陽(yáng) 110819)

巖體由各種各樣雜亂無(wú)章的結(jié)構(gòu)面和塊體組成，結(jié)構(gòu)面按成因可分為原生結(jié)構(gòu)面、次生結(jié)構(gòu)面和構(gòu)造結(jié)構(gòu)面.大量的工程實(shí)踐表明，巖體的破壞常常是沿著控穩(wěn)結(jié)構(gòu)面的破壞，在工程擾動(dòng)情況下，結(jié)構(gòu)面的力學(xué)效應(yīng)及巖坡和隧道的穩(wěn)定性主要取決于結(jié)構(gòu)面規(guī)模和工程地質(zhì)特性[1-2]，因此對(duì)結(jié)構(gòu)面進(jìn)行“分級(jí)”對(duì)于進(jìn)一步深入研究結(jié)構(gòu)面特性[3]及分析巖體穩(wěn)定性有重要意義.

目前，對(duì)于巖體分級(jí)[4-6]和圍巖分級(jí)[7-9]的研究較多，而針對(duì)結(jié)構(gòu)面的分級(jí)方法，主要為宏觀分析，如谷德振[10]在1979年依據(jù)結(jié)構(gòu)面的規(guī)模將其分為5級(jí)，然而該分級(jí)體系的宏觀性導(dǎo)致了其更適合于大范圍的區(qū)域地質(zhì)穩(wěn)定性研究，不便于指導(dǎo)工程實(shí)踐.張倬元等[11]又根據(jù)結(jié)構(gòu)面對(duì)巖體力學(xué)行為的影響程度，將其劃分為貫通性宏觀結(jié)構(gòu)面、顯現(xiàn)結(jié)構(gòu)面和隱微結(jié)構(gòu)面，該分級(jí)同樣較為宏觀.后來(lái)，向波等[12]在總結(jié)了前人研究的基礎(chǔ)上，加入了部分定量指標(biāo)，并結(jié)合相關(guān)國(guó)家標(biāo)準(zhǔn)將結(jié)構(gòu)面分為更加細(xì)致的5種級(jí)別.Sun等[13]結(jié)合四川省峨眉山某水電站工程，將現(xiàn)場(chǎng)結(jié)構(gòu)面細(xì)分為4組9個(gè)亞類(lèi)，為類(lèi)似地區(qū)結(jié)構(gòu)面的分級(jí)提供了參考.

上述研究均不同程度地推進(jìn)了結(jié)構(gòu)面分級(jí)體系的發(fā)展，但研究多為宏觀上的定性判斷，不便于指導(dǎo)工程設(shè)計(jì)和施工.鑒于此，本文將多指標(biāo)評(píng)價(jià)的理想點(diǎn)模型引入結(jié)構(gòu)面分級(jí)中，通過(guò)對(duì)模型權(quán)重求解進(jìn)行改進(jìn)，提出了一種基于改進(jìn)理想點(diǎn)模型的結(jié)構(gòu)面定量分級(jí)方法，并將該方法應(yīng)用于重慶市興隆隧道現(xiàn)場(chǎng)實(shí)測(cè)的典型結(jié)構(gòu)面分級(jí)中，將分級(jí)結(jié)果與實(shí)際工程進(jìn)行對(duì)比分析，驗(yàn)證了所提方法的可靠性和實(shí)用性.

1 結(jié)構(gòu)面分級(jí)的多指標(biāo)評(píng)價(jià)體系

影響結(jié)構(gòu)面分級(jí)的因素有許多，本文對(duì)結(jié)構(gòu)面的定量分級(jí)主要是結(jié)合向波等[12]提出的定性及半定量分級(jí)體系(見(jiàn)表1)，該分級(jí)體系中結(jié)構(gòu)面分級(jí)的主要影響因素有：①結(jié)構(gòu)面兩側(cè)巖壁巖性、強(qiáng)度及風(fēng)化程度；②結(jié)構(gòu)面的粗糙度；③結(jié)構(gòu)面充填物性質(zhì)；④結(jié)構(gòu)面富水程度.

表1 結(jié)構(gòu)面分級(jí)體系

根據(jù)以上所述各影響因素，為便于計(jì)算，指標(biāo)應(yīng)選取定量或半定量指標(biāo)，對(duì)于定性指標(biāo)，則應(yīng)進(jìn)行定量化處理，類(lèi)比“巖體分級(jí)”[4-6]及“圍巖分級(jí)”[7-9]中指標(biāo)選取方法和指標(biāo)獲取可行性，確定以下7個(gè)指標(biāo)對(duì)結(jié)構(gòu)面分級(jí)進(jìn)行定量表征，主要包括：泉流量(X1)、巖石基本質(zhì)量指標(biāo)(X2)、單軸飽和抗壓強(qiáng)度(X3)、巖體完整性系數(shù)(X4)、風(fēng)化系數(shù)(X5)、充填系數(shù)(X6)、粗糙系數(shù)(X7)，并以此為依據(jù)將結(jié)構(gòu)面分為5級(jí)，多指標(biāo)分級(jí)體系如圖1所示.

圖1 結(jié)構(gòu)面多指標(biāo)分級(jí)體系

2 結(jié)構(gòu)面分級(jí)的改進(jìn)理想點(diǎn)模型

理想點(diǎn)模型是由Hwang和Yoon[14]于1981年提出的基于多指標(biāo)的多方案分析方法.該方法操作簡(jiǎn)單，不需要復(fù)雜的計(jì)算過(guò)程，結(jié)果十分客觀[15].本文將其引入結(jié)構(gòu)面的分級(jí)中，充分利用原始結(jié)構(gòu)面數(shù)據(jù)信息，對(duì)結(jié)構(gòu)面進(jìn)行定量分級(jí)，并針對(duì)該方法權(quán)重確定較為主觀的缺點(diǎn)，提出了對(duì)該方法進(jìn)行改進(jìn)，利用層次分析法求得主觀權(quán)重；然后結(jié)合熵權(quán)法確定客觀權(quán)重，并利用最小熵原理優(yōu)化主、客觀權(quán)重，由此確定分級(jí)指標(biāo)的最優(yōu)組合權(quán)重.建立結(jié)構(gòu)面分級(jí)的改進(jìn)理想點(diǎn)模型，從而對(duì)結(jié)構(gòu)面進(jìn)行合理分級(jí)，為后續(xù)抗剪強(qiáng)度確定及巖體穩(wěn)定性分析奠定基礎(chǔ).其主要流程如圖2所示，為便于工程應(yīng)用，筆者基于Matlab平臺(tái)編制了該部分程序.

圖2 改進(jìn)理想點(diǎn)模型分級(jí)流程

2.1 結(jié)構(gòu)面分級(jí)指標(biāo)量化

對(duì)于泉流量(X1)，結(jié)合實(shí)際工程確定其定量評(píng)價(jià)指標(biāo)見(jiàn)表2.對(duì)于巖石基本質(zhì)量指標(biāo)(X2)、單軸飽和抗壓強(qiáng)度(X3)、巖體完整性系數(shù)(X4)，其指標(biāo)量化采用國(guó)標(biāo)分級(jí)標(biāo)準(zhǔn)，見(jiàn)表3.對(duì)于風(fēng)化系數(shù)(X5)、充填系數(shù)(X6)、粗糙系數(shù)(X7)，結(jié)合GSI巖體質(zhì)量分級(jí)法中的SCR部分確定其定量化評(píng)價(jià)指標(biāo)，見(jiàn)表4.

綜上，可得到結(jié)構(gòu)面分級(jí)各指標(biāo)的定量化匯總結(jié)果，見(jiàn)表5.

表2 結(jié)構(gòu)面富水程度評(píng)價(jià)指標(biāo)量化

表3 結(jié)構(gòu)面兩側(cè)巖壁質(zhì)量評(píng)價(jià)指標(biāo)量化

表4 結(jié)構(gòu)面表面特征評(píng)價(jià)指標(biāo)量化

表5 結(jié)構(gòu)面定量分級(jí)指標(biāo)

2.2 建立分級(jí)指標(biāo)矩陣

設(shè)有m個(gè)待分級(jí)結(jié)構(gòu)面樣本，含有n個(gè)分級(jí)指標(biāo)，對(duì)應(yīng)的評(píng)估值xij(i=1,2,…,m;j=1,2,…,n)，則分級(jí)指標(biāo)矩陣為

(1)

2.3 分級(jí)指標(biāo)矩陣標(biāo)準(zhǔn)化

由于結(jié)構(gòu)面各分級(jí)指標(biāo)具有不同量綱及量綱單位，為了消除該影響，需要對(duì)指標(biāo)進(jìn)行無(wú)量綱化，經(jīng)過(guò)標(biāo)準(zhǔn)化后的分級(jí)指標(biāo)矩陣為

(2)

對(duì)于收益型指標(biāo)：

對(duì)于消耗型指標(biāo)：

2.4 加權(quán)標(biāo)準(zhǔn)化分級(jí)指標(biāo)矩陣

2.4.1 層次分析法確定主觀權(quán)重

2.4.2 熵權(quán)法確定客觀權(quán)重

熵權(quán)法是一種客觀賦權(quán)法，它僅依賴(lài)于數(shù)據(jù)本身的離散性，目前已經(jīng)廣泛應(yīng)用于工程技術(shù)領(lǐng)域，采用該方法確定分級(jí)指標(biāo)客觀權(quán)重的具體步驟如下：

1)建立原始結(jié)構(gòu)面數(shù)據(jù)評(píng)價(jià)矩陣：

P=(pij)m×n.

(3)

式中，pij為指標(biāo)初始評(píng)價(jià)值.

2)歸一化原始矩陣，得到標(biāo)準(zhǔn)矩陣：

P*=(pij*)m×n，

(4)

(5)

3)設(shè)ωj′為第j個(gè)評(píng)價(jià)指標(biāo)的熵權(quán)，即所求客觀權(quán)重，則

(6)

2.4.3 確定最優(yōu)組合權(quán)重

為彌補(bǔ)單一賦權(quán)帶來(lái)的缺陷，本文利用層次分析法求得各分級(jí)指標(biāo)的主觀權(quán)重ω*，利用熵權(quán)法確定客觀權(quán)重ω′，ω′=(ω1′，ω2′，…，ωn′)，以此來(lái)修正主觀權(quán)重，最終得到最優(yōu)組合權(quán)重ω，為使三者盡可能接近，依據(jù)最小熵原理對(duì)主、客觀權(quán)重進(jìn)行組合，即

(7)

然后利用拉格朗日乘子法得到最優(yōu)組合權(quán)重為

(8)

2.4.4 加權(quán)標(biāo)準(zhǔn)化分級(jí)指標(biāo)矩陣

將最優(yōu)組合權(quán)重與標(biāo)準(zhǔn)化分級(jí)指標(biāo)矩陣相乘，得到加權(quán)標(biāo)準(zhǔn)化矩陣F=(fij)m×n=(ωjyij)m×n，即

(9)

2.5 構(gòu)建評(píng)價(jià)指標(biāo)體系

2.5.1 理想點(diǎn)和反理想點(diǎn)的確定

評(píng)價(jià)指標(biāo)一般分為收益型和消耗型[17-18]，收益型指標(biāo)越大越好，消耗型指標(biāo)越小越好，此時(shí)結(jié)構(gòu)面所屬等級(jí)也就越低.本文中泉流量屬于消耗型指標(biāo)，其余指標(biāo)均為收益型.若指標(biāo)單調(diào)變化，則可定義如下理想點(diǎn)和反理想點(diǎn).

當(dāng)指標(biāo)為收益型時(shí)，定義

(10)

當(dāng)指標(biāo)為消耗型時(shí)，定義

(11)

2.5.2 構(gòu)建理想點(diǎn)評(píng)價(jià)函數(shù)

各分級(jí)方案與理想點(diǎn)的距離采用歐氏距離來(lái)度量，具體含義為距離理想點(diǎn)越近，反理想點(diǎn)越遠(yuǎn)，分級(jí)結(jié)果就越好，具體公式為

(12)

2.6 貼近度分析

對(duì)于一組結(jié)構(gòu)面數(shù)據(jù)，當(dāng)分級(jí)指標(biāo)距理想點(diǎn)越近，距反理想點(diǎn)越遠(yuǎn)時(shí)，表明該結(jié)構(gòu)面等級(jí)越靠近理想點(diǎn)(即，等級(jí)越高)，該程度可以用貼近度β來(lái)表示，即

(13)

基于上述原理，可根據(jù)貼近度判斷各結(jié)構(gòu)面所屬等級(jí)，判斷規(guī)則見(jiàn)表6.

3 工程實(shí)例

重慶市三環(huán)高速公路合川至長(zhǎng)壽段興隆隧道位于重慶市渝北區(qū)木耳鎮(zhèn)良橋村，隧道在構(gòu)造上位于重慶-沙坪向斜北西翼與龍王洞背斜之間，隧址區(qū)圍巖巖性主要為侏羅系上沙溪廟組中風(fēng)化砂巖和泥巖，巖體結(jié)構(gòu)呈大塊狀，富水程度較低.

本文以現(xiàn)場(chǎng)實(shí)測(cè)的15條具有代表性的結(jié)構(gòu)面數(shù)據(jù)為例，進(jìn)一步驗(yàn)證改進(jìn)理想點(diǎn)模型應(yīng)用于結(jié)構(gòu)面分級(jí)中的可靠性和實(shí)用性，結(jié)合現(xiàn)場(chǎng)各工點(diǎn)原位試驗(yàn)結(jié)果，采用文獻(xiàn)[12]所述方法確定各結(jié)構(gòu)面所屬等級(jí).

3.1 確定無(wú)量綱化分級(jí)指標(biāo)矩陣

對(duì)實(shí)測(cè)的15條結(jié)構(gòu)面的各分級(jí)指標(biāo)按式(2)進(jìn)行無(wú)量綱化，得到無(wú)量綱化的分級(jí)指標(biāo)見(jiàn)表7.

3.2 確定分級(jí)指標(biāo)的組合權(quán)重

利用層次分析法，結(jié)合結(jié)構(gòu)面特征參數(shù)及專(zhuān)家意見(jiàn)，確定各分級(jí)指標(biāo)的主觀權(quán)重為

表7 實(shí)測(cè)結(jié)構(gòu)面分級(jí)指標(biāo)(無(wú)量綱化)

ω*=(0.108 7,0.110 7,0.142 5,
0.085 7,0.079 5,0.308 2,0.164 7).

然后利用式(3)～式(6)中熵權(quán)法計(jì)算得到各分級(jí)指標(biāo)客觀權(quán)重為

ω′=(0.031 7,0.085 0,0.280 8,
0.053 4,0.069 3,0.302 1,0.177 7).

然后由式(7)，式(8)確定最優(yōu)組合權(quán)重為

ω=(0.074 4,0.069 4,0.210 7,
0.026 5,0.049 5,0.389 6,0.189 9).

3.3 加權(quán)標(biāo)準(zhǔn)化分級(jí)指標(biāo)矩陣

確定最優(yōu)組合權(quán)重后，將表5中分級(jí)指標(biāo)無(wú)量綱化，并結(jié)合表7中無(wú)量綱化后的15條實(shí)測(cè)結(jié)構(gòu)面，代入式(9)可求得加權(quán)標(biāo)準(zhǔn)化分級(jí)指標(biāo)矩陣為

3.4 確定理想點(diǎn)和反理想點(diǎn)

將加權(quán)后的標(biāo)準(zhǔn)化矩陣代入式(10)，式(11),可求得理想點(diǎn)與反理想點(diǎn)為

3.5 確定貼近度

將3.3節(jié)及3.4節(jié)計(jì)算結(jié)果代入式(12)，式(13)，確定實(shí)測(cè)結(jié)構(gòu)面及各等級(jí)結(jié)構(gòu)面貼近度分別見(jiàn)表8和表9.

3.6 結(jié)果分析

結(jié)合表6中分級(jí)標(biāo)準(zhǔn)及表8和表9中貼近度計(jì)算結(jié)果，確定各條結(jié)構(gòu)面所屬分級(jí)結(jié)果見(jiàn)表10，各方法與實(shí)際分級(jí)數(shù)量的誤差對(duì)比結(jié)果如圖3所示.

表8 實(shí)測(cè)結(jié)構(gòu)面貼近度計(jì)算結(jié)果

表9 各等級(jí)結(jié)構(gòu)面貼近度計(jì)算結(jié)果

表10 結(jié)構(gòu)面定量分級(jí)結(jié)果

圖3 結(jié)構(gòu)面分級(jí)結(jié)果誤差對(duì)比

由表10及圖3中的對(duì)比結(jié)果可知，本文所提改進(jìn)理想點(diǎn)模型在應(yīng)用于結(jié)構(gòu)面分級(jí)時(shí)，僅第5條結(jié)構(gòu)面等級(jí)劃分結(jié)果(Ⅳ級(jí))略高于實(shí)際分級(jí)(Ⅲ級(jí))，整體分級(jí)結(jié)果與實(shí)際較為接近.進(jìn)一步分析出現(xiàn)分級(jí)偏差的結(jié)構(gòu)面可知，該結(jié)構(gòu)面貼近度計(jì)算結(jié)果為0.521 2，與分級(jí)標(biāo)準(zhǔn)中的0.508 3較為接近，分級(jí)時(shí)將介于兩級(jí)之間的結(jié)構(gòu)面判別為其中一級(jí)，事實(shí)上，由于結(jié)構(gòu)面等級(jí)越高，其性質(zhì)也越差，因此按高等級(jí)結(jié)構(gòu)面分析實(shí)際工程時(shí)更有利于保證施工安全，即判別結(jié)果偏于保守，也說(shuō)明了本文方法具有一定的可行性.另外，改進(jìn)理想點(diǎn)模型分級(jí)的準(zhǔn)確率為93.3%(1條分級(jí)出現(xiàn)偏差，14條分級(jí)準(zhǔn)確)，遠(yuǎn)優(yōu)于改進(jìn)前理想點(diǎn)模型(準(zhǔn)確率80%，3條分級(jí)出現(xiàn)偏差，12條分級(jí)準(zhǔn)確)，表明了本文方法應(yīng)用于結(jié)構(gòu)面分級(jí)中具有較好的工程實(shí)用性.

4 結(jié) 論

1)本文考慮結(jié)構(gòu)面分級(jí)的多種影響因素，利用7個(gè)常用指標(biāo)進(jìn)行定量化表征，建立了多指標(biāo)的結(jié)構(gòu)面定量分級(jí)體系，并依據(jù)最小熵原理對(duì)層次分析法及熵權(quán)法所求權(quán)重進(jìn)行優(yōu)化，避免了權(quán)重求解的主觀性，將所求最優(yōu)組合權(quán)重引入到理想點(diǎn)模型中，構(gòu)建了較為客觀的結(jié)構(gòu)面定量分級(jí)模型.

2)將改進(jìn)的理想點(diǎn)模型應(yīng)用于重慶市興隆隧道現(xiàn)場(chǎng)實(shí)測(cè)的結(jié)構(gòu)面分級(jí)中，其分級(jí)準(zhǔn)確率達(dá)93.3%，分級(jí)結(jié)果與實(shí)際情況更為吻合，且遠(yuǎn)優(yōu)于改進(jìn)前的準(zhǔn)確率(80%)，說(shuō)明所提方法實(shí)際可行，為巖體結(jié)構(gòu)面定量分級(jí)提供了一種新的方法.

3)理想點(diǎn)模型作為一種多目標(biāo)決策評(píng)價(jià)方法，在應(yīng)用其進(jìn)行結(jié)構(gòu)面分級(jí)時(shí)，評(píng)價(jià)指標(biāo)體系仍有待完善，嘗試用更科學(xué)的方法確定各指標(biāo)的動(dòng)態(tài)組合權(quán)重，將會(huì)使結(jié)構(gòu)面分級(jí)結(jié)果更為準(zhǔn)確且符合實(shí)際.