許 波, 謝謨文, 何興東, 胡 嫚

(1.北京科技大學(xué) 土木與環(huán)境工程學(xué)院, 北京 100083; 2.中國電建集團(tuán)昆明勘察設(shè)計研究院有限公司, 云南 昆明 650051; 3.西南大學(xué) 工程技術(shù)學(xué)院， 重慶 400715)

?

基于邊坡單元的水庫庫岸穩(wěn)定性層次分析與評價

許 波1, 謝謨文1, 何興東2, 胡 嫚3

(1.北京科技大學(xué) 土木與環(huán)境工程學(xué)院, 北京 100083; 2.中國電建集團(tuán)昆明勘察設(shè)計研究院有限公司, 云南 昆明 650051; 3.西南大學(xué) 工程技術(shù)學(xué)院， 重慶 400715)

摘要：[目的] 分析影響孟底溝水庫庫岸穩(wěn)定的影響因素,為水庫庫岸穩(wěn)定性評價提供新思路。[方法] 在工程地質(zhì)條件分析及不良地質(zhì)現(xiàn)象解譯的基礎(chǔ)上,運(yùn)用邊坡單元將各岸坡段劃分為獨(dú)立的研究對象,分析影響庫岸穩(wěn)定的諸多因素,把邊坡單元的三維安全系數(shù)作為穩(wěn)定性評價的基本指標(biāo),對影響邊坡穩(wěn)定的各因素采用層次分析法,建立基于邊坡單元的庫岸穩(wěn)定性分段層次分析與評價方法。[結(jié)果] 將該方法運(yùn)用在工程實(shí)例中,結(jié)果表明其具有一定的可靠性和實(shí)用性。[結(jié)論] 以邊坡單元為研究對象，運(yùn)用層次分析法可以有效地分析水庫庫岸的穩(wěn)定性，為工程決策提供支持。

關(guān)鍵詞：庫岸穩(wěn)定; 遙感解譯; 邊坡單元; 層次分析法; 三維分析

文獻(xiàn)參數(shù)： 許波, 謝謨文,何興東.基于邊坡單元的水庫庫岸穩(wěn)定性層次分析與評價[J].水土保持通報，2016,36(3)：290-296.DOI:10.13961/j.cnki.stbctb.2016.03.050

水庫蓄水過程及運(yùn)行階段，改變河谷岸坡的自然平衡條件后，庫岸巖土體的穩(wěn)定性稱為庫岸穩(wěn)定[1]。在此過程中，庫岸所處的地質(zhì)環(huán)境將發(fā)生顯著改變，自然平衡條件遭到破壞，引起岸坡沖刷、剝蝕、坍塌，庫岸逐漸后退，直到達(dá)到新的平衡狀態(tài)為止。由于大部分水庫都會有庫岸不穩(wěn)定問題[2]，國內(nèi)外很多專家學(xué)者都對庫岸穩(wěn)定性做了一些研究。

在庫水對岸坡的水巖作用研究方面：Fabio Luino[3]以降水誘發(fā)滑坡為對象，從不同角度對水巖(土)體系中地下水在邊坡失穩(wěn)破壞中的作用機(jī)制進(jìn)行了分析研究，對失穩(wěn)提出了不同的解釋；仵彥卿[4]總結(jié)了地下水對巖土體產(chǎn)生物理、化學(xué)以及力學(xué)的作用；馬水山[5]等從水巖作用的材料力學(xué)效應(yīng)、水力學(xué)效應(yīng)、化學(xué)效應(yīng)和地震效應(yīng)4個方面進(jìn)行了分析；王思敬等[6]論述了水庫地區(qū)水環(huán)境系統(tǒng)的變化及水巖作用類型和機(jī)制，在此基礎(chǔ)上討論了庫區(qū)地質(zhì)環(huán)境的變化和環(huán)境災(zāi)害的發(fā)生，提出了以水環(huán)境變化和水巖作用為中心的水庫地質(zhì)環(huán)境評價和災(zāi)害預(yù)測研究方向。在庫水作用下岸坡穩(wěn)定性的分析方法研究方面：Morgenstem在不考慮孔壓消散的假定基礎(chǔ)上，利用極限平衡法探討了庫水位變化對均質(zhì)岸坡安全系數(shù)的影響；梁燕、李同錄、李萍、陶福平等通過黃金分割方法優(yōu)化計算分析了庫水位下降時岸坡內(nèi)孔隙水壓力的變化及由此引起的岸坡穩(wěn)定性變化；劉新喜等[7]，朱冬林等[8]，Lane等[9]，利用強(qiáng)度折減法分析了水位變化對邊坡安全系數(shù)的影響。

在這些分析方法中，每一種方法所考慮的因素都是有限且固定的，而現(xiàn)場的地質(zhì)環(huán)境是復(fù)雜多變的，所以，這些方法在實(shí)際工程中的應(yīng)用還需要進(jìn)一步的完善。而層次分析法AHP[10]根據(jù)問題的性質(zhì)和要達(dá)到的總目標(biāo)，將問題分解為不同的組成因素，并按照因素間的相互關(guān)聯(lián)影響以及隸屬關(guān)系將因素按不同層次聚集組合，形成一個多層次的分析結(jié)構(gòu)模型，從而最終使問題歸結(jié)為最低層(供決策的方案、措施等)相對于最高層(總目標(biāo))的相對重要權(quán)值的確定或相對優(yōu)劣次序的排定。所以應(yīng)用層次分析法對水庫區(qū)庫岸穩(wěn)定性進(jìn)行評價是一種比較理想的方法。

本文以孟底溝水電站庫區(qū)為對象進(jìn)行庫岸穩(wěn)定研究。該水電站為四川省雅礱江中游第4個梯級水電站，位于甘孜州九龍縣與涼山州木里縣的交界的雅礱江干流上。研究區(qū)域位于青藏高原的東緣沙魯里山脈與大雪山之間，地處松潘—甘孜造山帶南部祝桑弧形構(gòu)造帶以南，位于“川滇巨型菱形斷塊”內(nèi)。

本文擬運(yùn)用邊坡單元的理論方法，并采用層次分析法來確定各評價指標(biāo)的權(quán)系數(shù)。將邊坡單元理論和層次分析法結(jié)合在一起，利用層次分析法確定影響邊坡單元穩(wěn)定的各因素的權(quán)系數(shù)，利用基于GIS研發(fā)的3D SlopeGIS模塊計算邊坡單元的三維安全系數(shù)，以其作為庫岸穩(wěn)定性評價的基本指標(biāo)，并參照邊坡單元實(shí)際的地形地貌、巖性構(gòu)造、不良地質(zhì)現(xiàn)象發(fā)育情況等因素予以修正，最終實(shí)現(xiàn)對整個庫區(qū)岸坡穩(wěn)定性評價。

1研究方法

1.1邊坡單元基本理論

大范圍研究區(qū)域內(nèi)的滑坡災(zāi)害分區(qū)方法一般是基于分區(qū)單元分析其危險度值或安全系數(shù)。分區(qū)單元可以是規(guī)則的網(wǎng)格單元或復(fù)雜的邊坡單元。規(guī)則網(wǎng)格單元的算法相對簡單，但這種均勻分布的網(wǎng)格單元與實(shí)際的邊坡破壞單元及力學(xué)機(jī)制、地質(zhì)、地形及水文沒有任何聯(lián)系。邊坡單元定義為一塊與鄰近區(qū)域具有明顯不同地形特征的區(qū)域[11]，而地形的形成本身可以反映地質(zhì)及水文地質(zhì)條件的長期作用效果，因此將邊坡單元作為分區(qū)單元的適用性是極強(qiáng)的，對基于三維力學(xué)模型的滑坡災(zāi)害分區(qū)及常用的多因素統(tǒng)計方法都是適用的。

一個邊坡單元實(shí)際上是一個曲面，可以用主傾斜方向、主傾斜角、邊坡寬度、邊坡長度、邊坡高度及中心點(diǎn)定義其屬性。

從地形的視覺上看，邊坡單元可以用特征地形線來進(jìn)行劃分。本文采用谷底線和山脊線來劃分邊坡單元，即一個匯水區(qū)的左右兩側(cè)各為一個邊坡單元。邊坡單元劃分的具體方法采用GIS水理分析功能，其步驟如下：

(1) 數(shù)字高程DEM柵格數(shù)據(jù)中凹點(diǎn)檢查和處理：對DEM柵格數(shù)據(jù)數(shù)值插值誤差引起的標(biāo)高極低(與周圍點(diǎn)相比)的孤立柵格進(jìn)行標(biāo)高充填處理。

(2) 水流方向計算：根據(jù)相鄰柵格的標(biāo)高比較得到水流方向。

(3) 水流匯集計算：從周圍流入到該單元的單元數(shù)。

(4) 匯水區(qū)，即山脊線的抽出：根據(jù)水流方向和水流匯集，即可抽出匯水區(qū)。這個匯水區(qū)的邊界即是要求的山脊線。但由于對一個大范圍的研究區(qū)域可以有不同規(guī)模的匯水區(qū)，即要求確定一個閥值來決定匯水區(qū)的規(guī)模。

為了得到谷底線的分布，我們采用以下的公式(1)對DEM柵格數(shù)據(jù)進(jìn)行處理，即可得到一個與原地形反向的地形柵格數(shù)據(jù)。在此反向的地形柵格數(shù)據(jù)中，原山脊線的位置變?yōu)楣鹊拙€，而原谷底線的位置變?yōu)樯郊咕€。這樣，基于這個反向的地形柵格數(shù)據(jù)采用同樣的水理分析流程，即可得到谷底線的分布。對谷底線和山脊線兩個多邊形數(shù)據(jù)集進(jìn)行聯(lián)合運(yùn)算，即可得到邊坡單元的數(shù)據(jù)集。

R-DEM=DEM×(-1)+DEMmax

(1)式中：R-DEM——反向的地形柵格數(shù)據(jù); DEM——原地形柵格數(shù)據(jù); DEMmax——原地形柵格數(shù)據(jù)的最大標(biāo)高值。

邊坡單元的大小選定取決于研究對象區(qū)域內(nèi)歷史滑坡的平均尺寸，即對于同一區(qū)域內(nèi)具有相同或類似的地形、地質(zhì)及水文條件的研究對象，可以認(rèn)為其將來的可能滑坡破壞與過去的滑坡破壞在規(guī)模上是一致的。由于在邊坡單元的劃分過程中，匯水集閥值決定了邊坡單元的大小，因此在實(shí)際應(yīng)用時，需要進(jìn)行多個閥值的對比試算以確定合適的閥值水平。

1.2層次分析法

層次分析法AHP是一種多指標(biāo)分析評價方法，具有精度高，使用方便的特點(diǎn)。作為一種多因素分析方法，雖然AHP有一定的主觀因素在里面，但是基于專家知識，逐一對比每兩個影響因子之間的關(guān)系，可以有效地將人為誤差降至最低。同時在確定各層次各因素之間的權(quán)重時，如果只是定性的結(jié)果，則常常不容易被別人接受，因而采用一致矩陣法，即：① 不把所有因素放在一起比較，而是兩兩相互比較；② 對比時采用相對尺度，以盡可能減少性質(zhì)不同的諸因素相互比較的困難，以提高準(zhǔn)確度。

1.3邊坡穩(wěn)定的控制參數(shù)

邊坡的穩(wěn)定性受諸多因素的影響，這些因素包括：地形地貌〔包括坡度坡向、巖體完整度、蓄水位(庫水位)等〕、巖性構(gòu)造(包括巖性、河谷結(jié)構(gòu)、河谷類型、岸坡結(jié)構(gòu)、斷層等)和不良地質(zhì)現(xiàn)象(包括滑坡、崩塌、堆積體、泥石流等)等。對于這些影響因素，通過層次分析法分析其對穩(wěn)定性的影響程度，結(jié)合三維安全系數(shù)評價其綜合穩(wěn)定性，建立的邊坡單元穩(wěn)定性綜合評估指標(biāo)如式(2)所示，其修正系數(shù)如表所示(由于在計算邊坡單元的三維安全系數(shù)時已考慮坡度、坡向等因素，故在地形地貌項中不再考慮)。

ISU=SF3D+RII

(2)

式中：ISU——邊坡單元穩(wěn)定性指數(shù)； SF3D——邊坡單元三維安全系數(shù)；RII——邊坡單元穩(wěn)定性修正指數(shù)。具體修正系數(shù)如表1所示。

表1　邊坡單元穩(wěn)定性綜合評估指標(biāo)修正系數(shù)

建立上述庫岸穩(wěn)定性評價指標(biāo)后，運(yùn)用層次分析法的原理，計算邊坡單元穩(wěn)定性指標(biāo)修正系數(shù)各層的指標(biāo)權(quán)重。

對比各指標(biāo)的相互重要程度，構(gòu)造判斷矩陣如表2所示。

表2　準(zhǔn)則層判斷矩陣Z

通過計算得出矩陣Z的最大特征值為：

λmax=3.053 6

對應(yīng)的權(quán)向量為：

W=(0.091 4，0.690 9，0.217 7)T

一致性指標(biāo)：CI=0.026 8

隨機(jī)一致性指標(biāo)：RI=0.58

一致性比率：CR=CI/RI=0.0462<0.1CR=CI/RI

通過一致性檢驗(yàn)。

1.3.1地形地貌對比各指標(biāo)的相互重要程度，構(gòu)造判斷矩陣如表3所示。

表3　地形地貌判斷矩陣A

通過計算得出矩陣A的最大特征值為：

λmax=3.094

對應(yīng)的權(quán)向量為：

W=(0.186 5，0.687 0，0.126 5)T

一致性指標(biāo)：CI=0.047 0

隨機(jī)一致性指標(biāo)：RI=0.58

一致性比率：CR=CI/RI=0.081 0<0.1

通過一致性檢驗(yàn)。

相對于目標(biāo)層，其權(quán)向量為：

W=(0.017 0，0.062 8，0.011 6)T

1.3.2巖性構(gòu)造對比各指標(biāo)的相互重要程度，構(gòu)造判斷矩陣如表4所示。

表4　巖性構(gòu)造判斷矩陣B

通過計算得出矩陣B的最大特征值為：

λmax=4.165 1

對應(yīng)的權(quán)向量為：

W=(0.252 8，0.059 7，0.084 3，0.603 2)T

一致性指標(biāo)：CI=0.069 4

隨機(jī)一致性指標(biāo)：RI=0.90

一致性比率：CR=CI/RI=0.077 1<0.1

通過一致性檢驗(yàn)。

相對于目標(biāo)層，其權(quán)向量為：

W=(0.174 7，0.041 3，0.005 82，0.416 7)T

1.3.3不良地質(zhì)現(xiàn)象對比各指標(biāo)的相互重要程度，構(gòu)造判斷矩陣如表5所示。

表5　不良地質(zhì)現(xiàn)象判斷矩陣C

通過計算得出矩陣C的最大特征值為：

λmax=4.173 8

對應(yīng)的權(quán)向量為：

W=(0.566 9，0.150 7，0.049 7，0.232 8)T

一致性指標(biāo)：CI=0.057 9

隨機(jī)一致性指標(biāo)：RI=0.90

一致性比率：CR=CI/RI=0.064 4<0.1

通過一致性檢驗(yàn)。

相對于目標(biāo)層，其權(quán)向量為：

W=(0.123 4，0.032 8，0.010 8，0.050 7)T

通過以上計算，得出各指標(biāo)的權(quán)重如表6所示，對于指標(biāo)層各指標(biāo)的取值如表7所示。

表6　邊坡單元穩(wěn)定性綜合評估指標(biāo)修正系數(shù)權(quán)重

表7　邊坡單元穩(wěn)定性修正指數(shù)指標(biāo)層各項目取值

綜上，可定義邊坡單元穩(wěn)定性修正指數(shù)RII如式(3)所示，可得邊坡單元穩(wěn)定性指數(shù)如式(4)所示。由于邊坡單元的穩(wěn)定指標(biāo)對庫岸段的貢獻(xiàn)與其規(guī)模正相關(guān)，而其規(guī)模與面積直接相關(guān)，故在得出邊坡單元的穩(wěn)定指標(biāo)后，可對庫岸段的邊坡單元的穩(wěn)定指標(biāo)進(jìn)行對面積的加權(quán)平均，作為庫段的穩(wěn)定性指標(biāo)，即岸段穩(wěn)定性指數(shù)如式(5)所示，為方便各庫段的比較評判，可把庫段穩(wěn)定性指數(shù)標(biāo)準(zhǔn)百分制化，如式(6)所示：

RIIi=∑xiyi

(3)

式中：xi——項目指標(biāo)值；yi——標(biāo)值對應(yīng)權(quán)重。

ISUi=SF3Di+RIIi

(4)

(5)

式中：Ai——邊坡單元的面積(m2)；ISUi——邊坡單元的穩(wěn)定性指數(shù)。

(6)

式中：Ii——庫岸段穩(wěn)定性指數(shù)；ISi——標(biāo)準(zhǔn)百分制化穩(wěn)定性指數(shù)。

1.4邊坡穩(wěn)定的分析

3D SlopeGIS是基于GIS開發(fā)的一個系統(tǒng)插件，可以對某一具體邊坡進(jìn)行綜合穩(wěn)定分析及加固分析，而且可以采用不同的模型計算邊坡三維安全系數(shù)，同時可以進(jìn)行最危險滑動面的搜索，此外，在考慮地震、地下水及外荷載等的情況下，研究邊坡在各種荷載組合下的三維穩(wěn)定性。

2結(jié)果與分析

2.1研究區(qū)庫岸分段

由于庫岸邊坡的穩(wěn)定性情況與該段的地形地貌、地層巖性、構(gòu)造和岸坡結(jié)構(gòu)等直接相關(guān)，其破壞形式也不一樣，如“V”型峽谷與“U”寬谷、沉積巖與侵入巖河段區(qū)分明顯；順向坡可能更多的產(chǎn)生滑坡破壞，而反向坡更多的可能是產(chǎn)生崩塌破壞；有斷層通過的河段穩(wěn)定性較差。故可根據(jù)以上對地形地貌、地層巖性、構(gòu)造和岸坡結(jié)構(gòu)等基本地質(zhì)條件的分析，把研究區(qū)分成12段(如附圖16所示)。

2.2研究區(qū)邊坡單元劃分

從地形學(xué)上看，邊坡單元可以用地形上的山脊線和谷底線來劃分。采用水文模型工具，經(jīng)過統(tǒng)合臨近要素，得到研究區(qū)域邊坡單元合理的劃分，最終結(jié)果如圖1所示。結(jié)合庫岸分段結(jié)果，可以得到每一個庫段的邊坡單元分布情況，本文以第一段的邊坡單元(圖2)為例詳細(xì)說明。同時對每個邊坡單元進(jìn)行編號。

2.3研究區(qū)邊坡單元三維安全系數(shù)

得到研究區(qū)的邊坡單元后，針對每一段的每一個邊坡單元進(jìn)行三維安全系數(shù)計算。此處運(yùn)用3D SlopeGIS邊坡三維分析軟件中的Janbu[12]模型計算每個邊坡單元的三維安全系數(shù)。由于研究區(qū)域受地下水影響不明顯，故在運(yùn)用3D SlopeGIS計算邊坡單

元穩(wěn)定性時將地下水直接予以考慮，而在運(yùn)用層次分析方法時不再考慮。此處依然以第一段的邊坡單元為例(如表8所示)。

圖1　最終的研究區(qū)邊坡單元(附山體陰影)

圖2　研究區(qū)RBS01段邊坡單元

編號SF3D編號SF3D編號SF3D編號SF3D11.49781.328151.335221.35521.23491.496161.521231.12231.339101.312171.485241.21141.428111.375181.401251.41951.217121.482191.37761.488131.492201.22171.478141.341211.203

2.4研究區(qū)庫岸穩(wěn)定綜合性評價

根據(jù)上文定義的邊坡穩(wěn)定性指數(shù)，可以計算得到各邊坡單元的穩(wěn)定性指數(shù)，同樣以第一段的邊坡單元為例(如表9所示)。結(jié)合第一段各邊坡單元的面積(表10)，可以得到第一段的穩(wěn)定性指數(shù)。按照同樣的步驟，可以得到其余各段的穩(wěn)定性指數(shù)(如表11所示)。

為方便研究區(qū)各庫段之間的比較評判，把庫段穩(wěn)定性指數(shù)標(biāo)準(zhǔn)百分制化后，結(jié)果如表12所示。

表9　RBS01段各邊坡單元的穩(wěn)定性指數(shù)

表10　RBS01段各邊坡單元面積

表11　各岸坡段的穩(wěn)定性指數(shù)

表12　庫岸段標(biāo)準(zhǔn)百分制化指數(shù)

按照工程上常用的對庫岸穩(wěn)定性劃分的四分法，把庫岸的穩(wěn)定性分成穩(wěn)定、基本穩(wěn)定、潛在不穩(wěn)定和不穩(wěn)定，定義90(含)分以上為穩(wěn)定段，80(含)～90分為基本穩(wěn)定段，70(含)～80分為潛在不穩(wěn)定段，70分以下為不穩(wěn)定段。

從表12可以看出，90分以上的庫段有RBS05，RBS08；80～90分的庫段有RBS01，RBS03，RBS04，RBS06，RBS07，RBS09，RBS12；70～80分的庫段為RBS02，RBS10，RBS11。

由以上劃分結(jié)果可得出庫區(qū)的穩(wěn)定段為RBS05，RBS08；基本穩(wěn)定段為RBS01，RBS03，RBS04，RBS06，RBS07，RBS09，RBS12；潛在不穩(wěn)定段為RBS02，RBS10，RBS11，無不穩(wěn)定岸坡，結(jié)果如圖3所示。

圖3　研究區(qū)庫岸穩(wěn)定性分段示意圖

3討論與結(jié)論

(1) 通過流程化的GIS數(shù)據(jù)處理方法，可以高效的把多源數(shù)據(jù)加入地理信息系統(tǒng)數(shù)據(jù)庫，不僅實(shí)現(xiàn)了對庫區(qū)基本地質(zhì)條件和不良地質(zhì)現(xiàn)象的解譯分析，而且便于與其他平臺進(jìn)行數(shù)據(jù)的交換。

(2) 通過對庫岸穩(wěn)定性的影響因素分析，將研究區(qū)劃分為獨(dú)立的邊坡單元，以此作為研究對象，計算其三維安全系數(shù)，作為庫岸穩(wěn)定評價的基本指標(biāo)。

(3) 采用層次分析法，對庫岸穩(wěn)定評價的基本指標(biāo)予以修正，以此建立了基于邊坡單元的定量三維穩(wěn)定性指標(biāo)與定性指標(biāo)相結(jié)合的庫岸穩(wěn)定性層次分析評價方法。

穩(wěn)定性評價指標(biāo)體系的建立包括評價指標(biāo)的選擇及指標(biāo)量化兩方面的內(nèi)容，這兩方面的選取均需要具有豐富的專業(yè)知識和現(xiàn)場經(jīng)驗(yàn)，本文部分指標(biāo)的選取和取值有待進(jìn)一步的商榷。結(jié)合現(xiàn)場的情況來看，最后的穩(wěn)定程度評價結(jié)果基本符合實(shí)際情況。

[參考文獻(xiàn)]

[1]中華人民共和國國家發(fā)展和改革委員會.DL/T5336—2006水電水利工程水庫區(qū)工程地質(zhì)勘察技術(shù)規(guī)程[S]. 2006.

[2]地質(zhì)礦產(chǎn)部編寫組.長江三峽工程庫岸穩(wěn)定性研究[M].北京：地質(zhì)出版社,1988.

[3]Luino F. Sequence of instability processes triggered by heavy rainfall in the northern Italy [J]. Geomorphology, 2005,66(1):13-39.

[4]仵彥卿.地下水與地質(zhì)災(zāi)害[J].地下空間與工程學(xué)報,1999,19(4):303-310.

[5]馬水山,雷俊榮,張保軍,等.滑坡體水巖作用機(jī)制與變形機(jī)理研究[J].長江科學(xué)院院報,2005,22(5):37-39.

[6]王思敬,馬鳳山,杜永廉.水庫地區(qū)的水巖作用及其地質(zhì)環(huán)境影響[J].工程地質(zhì)學(xué)報,1996,4(3):1-9.

[7]劉新喜,夏元友,張顯書,等.庫水位下降對滑坡穩(wěn)定性的影響[J].巖石力學(xué)與工程學(xué)報,2005,24(8):1439-1444.

[8]朱冬林,任光明,聶德新,等.庫水位變化下對水庫滑坡穩(wěn)定性影響的預(yù)測[J].水文地質(zhì)工程地質(zhì),2002,29(3):6-9.

[9]Lane P A, Griffiths D V. Assessment of Stability of Slopes under Drawdown Conditions[J]. Journal of Geotechnical & Geoenvironmental Engineering, 2000,126(5):443-450.

[10]孫宏才,田平,王蓮芬.網(wǎng)絡(luò)層次分析法與決策科學(xué)[M].北京:國防工業(yè)出版社,2011.

[11]謝謨文,蔡美峰.信息邊坡工程學(xué)的理論與實(shí)踐[M].北京:科學(xué)出版社,2005.

[12]Janbu N. Slope stability computations[C]∥Embankment-Dam Eng., Publication of Wiley & Sons Incorporated, New York, 1973:47-86.

收稿日期：2015-05-13修回日期：2015-07-13

通訊作者：胡嫚(1988—),女(漢族),貴州省畢節(jié)市人,博士,講師。主要從事地質(zhì)災(zāi)害調(diào)查研究、巖土工程數(shù)值模擬、防災(zāi)減災(zāi)的研究及教學(xué)工。E-mail:humanyyyes@126.com。

文獻(xiàn)標(biāo)識碼：A

文章編號：1000-288X(2016)03-0290-07

中圖分類號：P642.2

Slope Unit-based Application of Analytic Hierarchy Process in Stability Evaluation of Reservoir Bank

XU Bo1, XIE Mowen1, HE Xingdong2, HU Man3

(1.CollegeofGivilandEnvironmentalEngineering,UniversityofScienceandTechnologyBeijing,Beijing100083,China; 2.PowerchinaKunmingEngineeringCorporationLinited,Kunming,Yunnan650051,China; 3.CollegeofEngineeringandTechnology,SouthwestUniversity,Chongqing400715,China)

Abstract：[Objective] To analyze the influence factors of the Mengdigou reservoir bank stability in order to provide new train of thought to evaluate the stability of reservoir bank. [Methods] On the basis of the analysis of the engineering geological conditions and the interpretation of unfavorable geological phenomenon, this paper classified the bank into different segments. Then using the principle of the slope unit, we then taken these segments as independent research objects and analyzed the factors affecting the bank stability. And then we used three-dimension factor safety of slope units as an index of the stability evaluation and used analytic hierarchy process(AHP) to analyze the factors affecting slope stability. Finally we established hierarchy analysis and evaluation methods for the stability of reservoir bank based on the slope unit. [Results] The method was tested in the engineering practice, showing that the method is reliable and practical. [Conclusion] Tanking the slope unit as the research object, using the AHP can effectively analyze the stability of reservoir bank, which can support the decision of project.

Keywords:stability of reservoir bank; remote sensing; slope unit; analytic hierarchy process; three dimensional analysis

資助項目：國家自然科學(xué)基金項目“高山峽谷水庫區(qū)區(qū)域滑坡早期動態(tài)識別研究”(41372370)

第一作者：許波(1986—),女(漢族),安徽省淮南市人,博士研究生，主要研究方向?yàn)檫吰路€(wěn)定性、滑坡安全監(jiān)測。E-mail:xuboustb@163.com。