周秀全，黃海波，鄭 寧，施曉旺

浙江泰達(dá)安全技術(shù)有限公司，浙江 杭州 310000

0 引言

現(xiàn)階段我國(guó)部分水庫(kù)在蓄水后發(fā)生了一體系列次生地質(zhì)災(zāi)害，對(duì)當(dāng)?shù)丨h(huán)境和生產(chǎn)生活安全造成較大影響.根據(jù)現(xiàn)有文獻(xiàn)研究［1-3］，由于水庫(kù)蓄水造成的滑坡、崩塌地質(zhì)災(zāi)害已經(jīng)成為庫(kù)區(qū)地質(zhì)災(zāi)害的主要來(lái)源，并且這類地質(zhì)災(zāi)害受當(dāng)?shù)亟涤晏卣骱蛶?kù)水位變動(dòng)影響，具有復(fù)發(fā)性和周期性的特點(diǎn).另外，由于滑坡、崩塌等地質(zhì)災(zāi)害發(fā)生機(jī)理復(fù)雜，傳統(tǒng)的勘察手段和分析方法無(wú)法精確預(yù)測(cè)滑坡破壞機(jī)理，主流的滑坡防治手段仍然以預(yù)警和預(yù)防為主.滑坡一旦發(fā)生，其量能巨大，無(wú)法通過(guò)人工干預(yù)阻止其繼續(xù)發(fā)展，因此提前規(guī)避滑坡災(zāi)害及其影響區(qū)是最為經(jīng)濟(jì)的方案［4-8］.

針對(duì)庫(kù)水位滑坡地質(zhì)災(zāi)害的發(fā)生特點(diǎn)，本文依托于層次分析法，以某水庫(kù)滑坡體為研究對(duì)象，在對(duì)其地形地質(zhì)和水文地質(zhì)條件分析的基礎(chǔ)上，開(kāi)展滑坡至滑因子研究并建立滑坡易發(fā)性定量化評(píng)價(jià)方法.

1 研究區(qū)概況

研究區(qū)水庫(kù)位于山西省呂梁山脈西麓的南川河上，是一座以防洪灌溉為主，兼顧發(fā)電、養(yǎng)殖的年調(diào)節(jié)中型水庫(kù)，控制流域面積309 km2.隨著社會(huì)經(jīng)濟(jì)的發(fā)展，水庫(kù)下游鄉(xiāng)鎮(zhèn)、工礦企業(yè)的數(shù)量、規(guī)模迅速擴(kuò)大，城市生活、工業(yè)用水量驟增，而且需水方式也發(fā)生了明顯的變化，要求水庫(kù)的調(diào)節(jié)方式必須適應(yīng)發(fā)展的需要.由于庫(kù)區(qū)水位變動(dòng)幅度增大，現(xiàn)狀條件下水庫(kù)周邊巖體穩(wěn)定性受到一定程度影響，部分區(qū)域已經(jīng)發(fā)生明顯的崩塌滑坡跡象，一方面可能影響水庫(kù)正常運(yùn)行，另一方面威脅到附近居民的安全.

1.1 區(qū)域地質(zhì)

在區(qū)域地質(zhì)構(gòu)造上，研究區(qū)跨越呂梁-太行斷塊之次級(jí)構(gòu)造單元呂梁塊隆的離石-中陽(yáng)復(fù)向斜.該復(fù)向斜大致呈菱形，南北向展布，位于呂梁塊隆中西部，西鄰興縣-石樓南北向褶皺帶，東接關(guān)帝山穹狀隆起.該菱形復(fù)向斜大部分被第四系覆蓋，基巖一般只出露于溝谷，向斜軸部地層為石炭系—二疊系，翼部地層為寒武系—奧陶系.

研究區(qū)出露的地層主要有太古宇呂梁群、寒武系中統(tǒng)徐莊組及新生界第四系.巖性包括變粒巖、混合花崗巖、黑云斜長(zhǎng)片麻巖和角閃斜長(zhǎng)片麻巖，總厚度大于1000 m，在區(qū)內(nèi)大面積出露，是構(gòu)成基巖山體的主要地層.第四系地層是庫(kù)周及庫(kù)底的主要上部巖土層，分布于河道及溝谷，主要有棕紅色低液限黏土、灰黃色低液限粉土、低液限粉土、級(jí)配不良砂及卵石混合土層，一般具二元結(jié)構(gòu)，厚度5～30 m.

1.2 庫(kù)區(qū)地質(zhì)

水庫(kù)壩址位于南川河上，壩址以上河流長(zhǎng)30.2 km，河道縱坡15.4%，壩址附近河道縱坡12.2‰，河道穩(wěn)定.南川河庫(kù)區(qū)段整體走向334°，河床寬80～350 m.兩岸發(fā)育有漫灘及Ⅰ級(jí)階地，在近壩區(qū)寬30～60 m.庫(kù)區(qū)沿河兩岸除局部覆蓋有第四系中、上更新統(tǒng)黃土外，基本是由呂梁群混合花崗巖、變粒巖、混合巖化黑云斜長(zhǎng)片麻巖、角閃斜長(zhǎng)片麻巖構(gòu)成基巖山體，山頂高程均在1200 m以上.

庫(kù)區(qū)發(fā)育的地層主要為太古宇呂梁群、寒武系徐莊組及第四系.庫(kù)區(qū)主要淤積物為級(jí)配不良砂、低液限黏土及低液限粉土，厚1～17.5 m，主要分布于近壩庫(kù)岸.壩體堆積物為低液限粉土、低液限黏土，偶含小礫石，結(jié)構(gòu)松散，厚0～36.5 m.

研究區(qū)地質(zhì)構(gòu)造較為簡(jiǎn)單，主要發(fā)育兩條主要斷層：F1逆斷層位于庫(kù)區(qū)中廟附近，走向北北東并橫跨河床，傾向南東東，地表處傾角約75°，東岸廟后溝口可見(jiàn)斷層上下盤均為太古宇變質(zhì)巖，垂直斷距大于50 m；F2逆斷層位于壩軸線下游600 m處，走向北西并橫跨河床，傾向西南，地表處傾角62°，上盤為太古宇變質(zhì)巖，下盤為寒武系頁(yè)巖，接觸處有角礫巖，斷層帶寬約1 m，垂直斷距大于50 m.

庫(kù)區(qū)內(nèi)地下水可分為變質(zhì)巖類裂隙水及松散巖類孔隙水兩種類型.變質(zhì)巖類裂隙水貯存于巖層風(fēng)化帶及裂隙中，補(bǔ)給主要為大氣降水.排泄為向下游的徑流或向溝谷的側(cè)向徑流.松散巖類孔隙水含水層巖性主要為全新統(tǒng)低液限粉土及卵石混合土，大氣降水及兩岸山體地下水是主要的補(bǔ)給來(lái)源，地下水徑流方向與地形傾斜方向基本一致，向下游河谷低凹地帶徑流.

1.3 庫(kù)區(qū)滑坡地質(zhì)災(zāi)害

庫(kù)區(qū)兩岸支流和沖溝眾多，庫(kù)盆呈樹(shù)枝狀.

上游庫(kù)段總體為較開(kāi)闊的峽谷型庫(kù)盆，正常蓄水位高程以上岸坡多較陡，覆蓋層厚度變化比較大，基巖出露普遍，但風(fēng)化強(qiáng)烈.庫(kù)岸以巖質(zhì)岸坡為主，局部為土質(zhì)岸坡.兩岸沖溝較多，且比較寬長(zhǎng).庫(kù)尾（江邊寨）、河流急轉(zhuǎn)彎凸岸、支流和沖溝的出口交匯處等部位河谷谷底一般比較寬闊，形成比較開(kāi)闊的庫(kù)盆，沉積有較厚土層.該庫(kù)段岸邊村寨眾多，居民房屋較多；同時(shí)，有多條高速公路及鐵路通過(guò)庫(kù)區(qū)，庫(kù)岸失穩(wěn)的危害性較大.

下游庫(kù)段為比較狹窄的峽谷型庫(kù)盆，正常蓄水位高程910 m以上基本為陡峻基巖岸坡，覆蓋層較薄，基本為巖質(zhì)岸坡.兩岸沖溝少且短，右岸發(fā)育有谷頂河.河流急轉(zhuǎn)彎凸岸、支流和沖溝的出口交匯處等部位不存在寬闊谷底.該庫(kù)段的村寨多分布在兩岸山頂，岸坡上居民房屋極少，也沒(méi)有高速公路、鐵路等通過(guò)，庫(kù)岸失穩(wěn)的危害性很小.

綜合判斷庫(kù)區(qū)范圍內(nèi)產(chǎn)生的庫(kù)岸失穩(wěn)形式主要為塌岸、塌滑體，局部基巖陡壁存在小規(guī)模崩塌的可能.

2 基于層次分析法的評(píng)價(jià)思路

層次分析法（Analytic Hierarchy Process，AHP）首先應(yīng)用于安全與系統(tǒng)工程中的風(fēng)險(xiǎn)分析，其主要思路是將決策某事件發(fā)生的因素進(jìn)行劃分并量化，在此基礎(chǔ)上進(jìn)行決策分析［9-10］.層次分析法每一級(jí)矩陣代表一個(gè)影響事件發(fā)生的因子，通過(guò)實(shí)地調(diào)研分析對(duì)每個(gè)因子賦予權(quán)重值，然后進(jìn)行打分從而得到事件發(fā)生的概率.層次分析法的基本工作流程如圖1.

圖1 層次分析法工作流程Fig.1 Workflow of analytic hierarchy process

通過(guò)對(duì)該水庫(kù)基本地形地質(zhì)條件分析認(rèn)為，影響滑坡穩(wěn)定的災(zāi)害因素包括地形地貌、地質(zhì)、氣象、水文、構(gòu)造以及水位變化等，其中庫(kù)水位變化過(guò)程是影響研究區(qū)滑坡危險(xiǎn)性的重要考量因素［11-15］.綜合分析，庫(kù)區(qū)主要的滑坡致害因子分為災(zāi)害發(fā)育程度、地形地質(zhì)情況和誘發(fā)因子3個(gè)大類（分別用B1、B2、B3表示）.其中，災(zāi)害發(fā)育程度包括災(zāi)害點(diǎn)密度和災(zāi)害發(fā)展面積；地形地質(zhì)情況包括庫(kù)區(qū)地層巖性、分布高程、坡向、坡度；誘發(fā)因子包括庫(kù)水消落高度、地震、降雨強(qiáng)度（由C1—C9表示，見(jiàn)圖2）.

圖2 基于層次分析法的滑坡災(zāi)害因子指標(biāo)Fig.2 Landslide hazard factor index based on AHP

3 基于層次分析法的庫(kù)內(nèi)滑坡穩(wěn)定性評(píng)價(jià)

3.1 評(píng)價(jià)指標(biāo)的確定

合理的評(píng)價(jià)指標(biāo)是確保層次分析法精確度的重要前提.通過(guò)研究大量文獻(xiàn)中有關(guān)庫(kù)水位變動(dòng)過(guò)程中的滑坡穩(wěn)定性問(wèn)題［16］，綜合確定了影響庫(kù)水位變動(dòng)下滑坡災(zāi)害發(fā)生的9種主要因子并對(duì)其等級(jí)進(jìn)行劃分（表1）.

表1 致滑因子及其量化指標(biāo)綜合分類表Table 1 Comprehensive classification of landslide factors and their quantitative indexes

3.2 評(píng)價(jià)矩陣構(gòu)建

本研究采用兩層判斷矩陣法，即構(gòu)造A-B和BC.每個(gè)判斷矩陣之間的相互關(guān)系通過(guò)以往水庫(kù)區(qū)滑坡災(zāi)害的研究認(rèn)識(shí)和遙感解譯的統(tǒng)計(jì)分析成果來(lái)綜合判斷.為了直觀反映判斷各個(gè)子因素C1—C9之間的關(guān)聯(lián)過(guò)程，采用運(yùn)籌學(xué)的標(biāo)度法［17-20］（表2）.

表2 運(yùn)籌標(biāo)度指標(biāo)Table 2 Scale indexes of logistics operation research

通過(guò)對(duì)災(zāi)害發(fā)育程度、地形地質(zhì)和誘發(fā)因素3個(gè)B矩陣與危險(xiǎn)性A之間的關(guān)聯(lián)度分析，建立了A-B矩陣如下：

然后，分別將B1—B3矩陣與C之間的關(guān)系采用同樣方法建立B-C矩陣：

其中，B1-C表示災(zāi)害發(fā)育程度與災(zāi)害點(diǎn)分布、災(zāi)害點(diǎn)密度之間的相關(guān)性；B2-C表示地形地質(zhì)因素與地貌、巖性、坡度和坡向之間的相關(guān)性；B3-C表示誘發(fā)因素與地震、降雨和庫(kù)水變動(dòng)之間的相關(guān)性.

為了證明該矩陣的合理性，通過(guò)單層排序法計(jì)算每個(gè)判斷矩陣所產(chǎn)生的特征根及其特征向量，然后利用特征向量建立由C到B、由B到A的層次單排序計(jì)算.最后，依據(jù)計(jì)算結(jié)果對(duì)每個(gè)判斷矩陣進(jìn)行一致性檢驗(yàn)，若結(jié)果小于0.1，表明該判斷矩陣是合理可信的，若不滿足小于0.1則重新調(diào)整矩陣結(jié)構(gòu)和判斷方式，直至完成一致性檢驗(yàn).具體步驟如下.

1）對(duì)矩陣C采用和積法對(duì)每列向量進(jìn)行歸一化處理：

2）對(duì)歸一化后得到的矩陣Cg相加求和：

3）將求和的Ch矩陣再次歸一化得到特征向量：

4）利用特征向量W獲取特征根：

5）進(jìn)行一致性檢驗(yàn)，其判別公式如下：

其中，Rc表示一致性檢驗(yàn)計(jì)算比率，R表示隨機(jī)一致性指標(biāo)，其取值參考表3.

表3 隨機(jī)一致性指標(biāo)R取值Table 3 Random consistency index R value

通過(guò)以上分析步驟得到10個(gè)判斷矩陣之間的特征向量和特征根如下：

A-B矩陣中，W1＝（0.122，0.503，0.347），λ＝2.97，計(jì)算求得其Rc值為0.01，滿足一致性檢驗(yàn)條件；

B1-C矩陣中，W2＝（0.45，0.45），λ＝1.99，計(jì)算得Rc＝0.02，滿足一致性檢驗(yàn)標(biāo)準(zhǔn)；

B1-C矩陣中，W3＝（1.66，0.99，0.99，0.50），λ＝3.99，計(jì)算的Rc＝0.004，滿足一致性檢驗(yàn)標(biāo)準(zhǔn)；

B1-C矩陣中，W1＝（1.82，0.66，0.35），λ＝3.97，計(jì)算求得Rc＝0.03，同樣滿足一致性檢驗(yàn)標(biāo)準(zhǔn).

一致性檢驗(yàn)合格后綜合層次分析成果，獲取在危險(xiǎn)性分析過(guò)程中B1—B3和C1—C9的權(quán)重值，見(jiàn)表4.

從表4可以看出，對(duì)庫(kù)區(qū)滑坡影響最大的權(quán)重因子分別為C1（地災(zāi)點(diǎn)密度）、C3（地層巖性）和C8（庫(kù)水變動(dòng)幅度），其權(quán)重值分別為0.27、0.21和0.33.這一規(guī)律在滑坡地災(zāi)分析的宏觀文獻(xiàn)上能夠得到驗(yàn)證.

表4 水庫(kù)滑坡區(qū)至滑因子權(quán)重評(píng)估成果Table 4 Weight evaluation results of landslide factors in reservoir landslide area

3.3 滑坡危險(xiǎn)性分區(qū)

為了分析庫(kù)區(qū)蓄水后不同滑坡的危險(xiǎn)性，建立基于層次分析權(quán)重的危險(xiǎn)性定量評(píng)估模型：

其中，Mi表示庫(kù)區(qū)不同部位的滑坡危險(xiǎn)性評(píng)估值；Ci表示不同因子的權(quán)重值；λi表示不同影響因子的災(zāi)害危險(xiǎn)性賦值.

為了進(jìn)一步提高庫(kù)區(qū)滑坡危險(xiǎn)性評(píng)估成果的直觀可視化程度，建立庫(kù)區(qū)基于滑坡危險(xiǎn)性等級(jí)的指標(biāo)，庫(kù)區(qū)滑坡發(fā)生概率分為極高、高、中等、低共4個(gè)等級(jí)，然后利用MapGIS疊加功能生成滑坡區(qū)危險(xiǎn)性分區(qū)圖（圖3）.從圖3可以發(fā)現(xiàn)，庫(kù)區(qū)左岸中游與右岸中游為滑坡地質(zhì)災(zāi)害極高危險(xiǎn)地區(qū)；水庫(kù)中上游和庫(kù)尾段為滑坡易發(fā)高危區(qū)，庫(kù)區(qū)上游兩側(cè)為滑坡地質(zhì)災(zāi)害低發(fā)區(qū).

圖3 基于層次分析法的庫(kù)區(qū)滑坡危險(xiǎn)性分區(qū)Fig.3 Landslide risk zoning in reservoir area based on AHP

為了驗(yàn)證計(jì)算模型設(shè)計(jì)的合理性，調(diào)取研究區(qū)衛(wèi)星遙感圖像的地質(zhì)災(zāi)害點(diǎn)解譯成果（見(jiàn)圖4），可以看出，研究區(qū)滑坡地質(zhì)災(zāi)害極高危地段（圖中紅色部位）均存在地質(zhì)災(zāi)害發(fā)生點(diǎn)，說(shuō)明滑坡危險(xiǎn)性區(qū)劃圖能夠較好地反映現(xiàn)狀條件下庫(kù)區(qū)潛在滑坡體可能的分布范圍.建議重點(diǎn)關(guān)注紅色和黃色區(qū)域的地質(zhì)災(zāi)害問(wèn)題，加強(qiáng)巡查并落實(shí)合理的滑坡治理措施.

圖4 研究區(qū)遙感解譯地質(zhì)災(zāi)害發(fā)生點(diǎn)分布圖Fig.4 Distribution map of geohazard sites in the study area interpreted by remote sensing

4 結(jié)論

為研究和分析水庫(kù)蓄水后的庫(kù)內(nèi)滑坡地質(zhì)災(zāi)害潛在危險(xiǎn)性和易發(fā)程度，開(kāi)展了基于層次分析法的庫(kù)區(qū)滑坡易發(fā)性評(píng)價(jià)研究，結(jié)論如下.

1）庫(kù)區(qū)上游庫(kù)段總體為較開(kāi)闊的峽谷型庫(kù)盆，正常蓄水位高程以上岸坡多較陡，覆蓋層厚度變化比較大，基巖出露普遍，但風(fēng)化強(qiáng)烈.下游庫(kù)段為比較狹窄的峽谷型庫(kù)盆，村寨多分布在兩岸山頂，岸坡上居民房屋極少，也沒(méi)有高速公路、鐵路等通過(guò)，庫(kù)岸失穩(wěn)的危害性很小.綜合判斷庫(kù)區(qū)范圍內(nèi)產(chǎn)生的庫(kù)岸失穩(wěn)形式主要為塌岸、塌滑體.

2）庫(kù)區(qū)左岸中游與右岸中游為滑坡地質(zhì)災(zāi)害極高危險(xiǎn)地區(qū)，水庫(kù)中上游和庫(kù)尾段為滑坡易發(fā)高危區(qū)，庫(kù)區(qū)上游兩側(cè)為滑坡地質(zhì)災(zāi)害低發(fā)區(qū).研究區(qū)衛(wèi)星遙感圖像的地質(zhì)災(zāi)害點(diǎn)解譯成果表明滑坡地質(zhì)災(zāi)害極高危地段均存在地質(zhì)災(zāi)害發(fā)生點(diǎn)，說(shuō)明本研究獲得的滑坡危險(xiǎn)性區(qū)劃圖能夠較好地反映現(xiàn)狀條件下庫(kù)區(qū)潛在滑坡體可能的分布范圍.

3）本研究為庫(kù)區(qū)蓄水條件下的滑坡地質(zhì)災(zāi)害易發(fā)性半定量化研究提供了一定參考和借鑒，今后的研究中將進(jìn)一步細(xì)化滑坡因子數(shù)量和權(quán)重值的計(jì)算合理性，為滑坡易發(fā)性評(píng)價(jià)提供更為豐富的計(jì)算參考.