三峽庫首區滑坡空間分布特征分析及危險性評價

李　雪　李井岡　劉小利　王秋良

1　中國地震局地震研究所(地震大地測量重點實驗室)，武漢市洪山側路40號，430071

三峽庫首區滑坡空間分布特征分析及危險性評價

李雪1李井岡1劉小利1王秋良1

1中國地震局地震研究所(地震大地測量重點實驗室)，武漢市洪山側路40號，430071

摘要：采用遙感和GIS技術從宏觀上對三峽庫首區三峽工程二期、三期治理滑坡和新增滑坡的空間分布與地形地貌、地層巖性、斷裂構造、水系流域和人為影響等因素的相關程度進行分析，獲取孕災環境特征規則，并對其危險性進行評價。結果表明，地形因子中的高程、坡向和高程變異系數與滑坡的相關性較高；研究區內斷裂構造對滑坡發育的控制作用不明顯；研究區滑坡體的空間分布受到不同因子組合的約束，可根據不同因子組合約束的權重對研究區滑坡危險性進行評價。

關鍵詞：三峽水庫；遙感；空間分布；相關性；滑坡危險性分析

三峽工程建設以來，地質災害隱患點明顯增多，地質災害險情發生頻繁。雖然通過二期、三期地質災害防治監測預警工程[1-3]，已治理滑坡的危害性大為降低，但滑坡體所處位置仍具備孕災特征，必須對二期、三期治理滑坡和新增滑坡進行長期監測，并根據情況開展防治預警工作。

當前滑坡災害防治監測預警的主要方法有GPS監測、地下水監測、滑坡深部位移監測和應力應變監測等[4]。盡管專業監測的精度更高、更準確，但成本高昂，難以大范圍開展。而結合遙感與GIS空間分析方法可以有效提取孕育滑坡的環境因子，并對滑坡的危險性進行評價，為滑坡監測預警提供有力支持[5]。

1方法技術

1.1滑坡空間分布分析的技術流程

利用遙感和GIS中的特征提取和空間分析技術，從宏觀上對滑坡體的空間分布特征進行分析，并對滑坡體與地形地貌、地層巖性、斷裂構造、水系流域和人為影響等因素的相關性進行分析。通過對孕育滑坡的環境影響因素的定量分析，提取特征因子并確定因子與滑坡之間的相關性。再對相關因子進行模式分析，獲取孕育滑坡的環境規則，最后根據規則對不同區域的滑坡易發性進行分級，評價研究區滑坡危險性。滑坡空間分布分析的技術流程如圖1所示。

圖1　滑坡空間分布分析技術流程Fig.1　Flow chart of landslide spatial distribution analysis

1.2研究區域概況

研究區域(110.4°～110.9°E，30.8°～31.1°N)位于三峽庫首區秭歸至巴東段，面積約1 600 km2，地貌類型以侵蝕溶蝕與侵蝕剝蝕型為主，區內地層出露齊全，發育良好，除缺失泥盆系下統、石炭系上統、白堊系的一部分和第三系外，自前震旦系至第四系均有出露[6]。研究區以長江為骨干，組成樹枝狀水文網，江河水量充沛[7]。

搜集研究區內三峽工程二期搬遷避讓、監測預警和治理滑坡41處，三期搬遷避讓、監測預警和治理滑坡20處，新增滑坡16處，共計77處。其中47處滑坡體前緣低于145 m，19處滑坡體前緣介于145～175 m，11處滑坡體前緣高于175 m，所有滑坡體后緣均高于175 m。

用于提取滑坡體特征因子的輔助數據包括具有30 m多光譜和15 m全色分辨率的Landsat ETM遙感影像，具有20 m垂直精度、30 m水平精度的ASTER GDEM數字高程模型，以及研究區基礎地質和基礎地理信息數據。結合遙感影像和湖北省地質礦產局測制的1∶20萬地質圖(巴東幅)，獲得研究區地質構造略圖(圖2)。

圖2　研究區地質構造略圖Fig.2　Geological sketch map of the studied area

2滑坡分布因子分析

2.1滑坡因子概率分析

地形因素是孕育滑坡的重要影響因素。利用數字高程模型(DEM)可以提取滑坡體的高程、坡度、坡向等地形因子[8]，通過對地形因子的統計分析，可以獲取滑坡發生的地形特征。然而，并不是所有地形因子對滑坡孕育的影響都同等重要。王秋良[9]對北川地區地震觸發滑坡災害的地形因子分析后認為，坡度、坡度變率和地形起伏度等地形因子對該區域影響較大。夏金梧[10]對長江上游地區滑坡分布特征分析后認為，坡度和地表切割深度等地形因子對滑坡孕育的影響最大。

本文根據滑坡目錄提取了77處滑坡點的12種地形因子，分別為高程、高程變異系數、坡度、坡度變率、坡向、坡向變率、坡形、地表切割深度、地形起伏度、地面粗糙度、平面曲率、剖面曲率。為檢驗這些地形因子與滑坡的相關程度，在研究區隨機選取100處非滑坡點，并提取其地形因子，分別比較滑坡點和隨機點12種地形因子直方圖曲線的相似度發現，除高程、坡向和高程變異系數3項地形因子差別較大外，其他9項地形因子在滑坡點與隨機點的直方圖曲線十分相似。這表明，高程、坡向和高程變異系數相比于其他9項地形因子，與滑坡發生的相關性更高。因此，本文選擇高程、坡向和高程變異系數作為后續滑坡因子分析的地形因子。

2.2滑坡因子貢獻分析

滑坡的孕育與發生除跟地形因子有關外，還與地層巖性、地質構造、水系流域等自然環境和人為影響等社會環境有關[11]，不同類型影響因子對滑坡分布的貢獻各不相同[12-14]。因此，我們選擇的滑坡因子除地形因子外，還包括地層巖性因子、斷裂構造因子、水系流域因子和人為影響因子。其中地層巖性因子為地質構造略圖(圖2)中的12種主要地層時代和巖性；斷裂構造因子為滑坡點距最近斷裂構造的距離；水系流域因子為滑坡點距最近河流、水系的距離；人為影響因子為滑坡點距最近居民點的距離。

本文采用決策樹分析中的信息增益率計算不同影響因子對滑坡分布分析的貢獻:

(1)

其中，A為樣本的特征空間；S為樣本隨機變量；H(S)與H(A)為信息熵函數，用于衡量隨機變量S和A中的信息量；H(S|A)為條件熵函數，表示當A取值已知時隨機變量S的不確定性程度(信息量)。采用決策樹算法進行分析時，每次計算子節點劃分閾值時都取A0，使得當前狀況下GainRatio(S,A0)的值為最大。假設樣本的某一特征在決策樹分析中被作為節點計算了m次，則該特征的貢獻度為：

Contribution(Ai)=

(2)

本文利用77處滑坡點的特征因子作為正例，100處隨機點的特征因子作為反例，代入C5.0決策樹算法，對不同特征因子在滑坡分析中的貢獻度進行分析，獲得如圖3所示的因子貢獻分析決策樹。生成的決策樹具有16個節點和19個葉子節點，正例判別正確率97.4%，反例判別正確率100%，總體精度98.9%。

3區域滑坡危險性分析

根據滑坡因子貢獻分析獲得的決策樹，可以提取出滑坡分布規則。從決策樹的根節點出發，沿子節點到達判別為滑坡的葉子節點，路徑上所有節點判別條件的集合為一條規則。從圖3可知，決策樹共存在11個判別為滑坡的葉子節點，故存在11條滑坡分布規則。對每條規則路徑上的條件進行整理、歸并后獲得如表1所示的滑坡分布規則。利用滑坡分布規則對研究區范圍內的滑坡因子進行判別分析，可以得到滿足滑坡分布規則的二值化圖像。圖4為11條滑坡分布規則及其疊加后的判別結果，圖中黑色為滑坡危險區，白色為非危險區。由圖可見，不同規則判別出的危險區范圍各不相同，這是由于決策樹算法具有判別的唯一性，即一個樣本對應唯一的葉子節點。將11條規則判別出的危險區疊加后，即獲得規則疊加分析結果。

圖3　滑坡因子貢獻分析決策樹Fig.3　Decision tree of landslide factors contribution analysis

規則編號規則條件GCDCBYXSWDPXRD1(515,613]2<=0.0086>752(515,613]2<=753<=364{3,4,8,9,10,12}4<=3642<=0.04765<=3642>0.0476<=41956<=3642>0.0476>208>41957(364,515]{8,11}<=11748(364,515]3<=1174>1549(497,515]3<=1174<=15410(364,515]2<=1174>7411(464,515]2<=1174<=74

圖4　滑坡規則危險性判別結果Fig.4　Results of landslide risk identification

由于判別規則生成的結果為0/1影像，沒有滑坡危險程度區分。本文對11條規則判別的滑坡危險區進行統計分析，計算了不同判別規則下滑坡數量比值(R1)和滑坡面積比值(R2)。根據兩種比值計算不同規則判別出的危險區的權重(R)：

(3)

其中，函數N(x)=(x-min(x))/(max(x)-min(x))為歸一化函數，歸一區間為[0,1]。11條規則判別的滑坡危險區統計數據和其相應權重如表2所示。

表2　滑坡危險區統計

根據滑坡規則判別危險區的權重值。將滑坡規則的危險等級由高到低劃分為5級：規則4判別的滑坡區危險性等級最高，面積約為89.1 km2，主要分布于香溪河右岸，涼臺河、童莊河兩側，鑼鼓洞河兩河口鎮至長江干流段兩側以及長江干流北岸零星分布；規則3判別的滑坡區危險性等級次之，面積約為130.2 km2，主要分布于板橋河、青干河兩側，香溪河左岸，以及長江干流兩側等區域；規則{5,9,10,11}和規則{1,2,6,7,8}判別的滑坡區危險性等級相對較小，面積分別為89.1 km2和61.4 km2；規則判別以外的區域滑坡危險性極低，面積約1 257.2 km2。研究區不同等級滑坡危險區分布如圖5所示。

圖5　研究區滑坡危險性分析結果Fig.5　Result of landslide risk analysis

4結語

1)研究區地形因子中，高程、坡向和高程變異系數與滑坡體分布的相關性較高，而坡度、地形起伏度等其余9項因子的概率分布與研究區隨機點的概率分布一致，說明其與滑坡體分布的相關性較小；

2)經數據挖掘模式分析，研究區內滑坡的空間分布具有11條規則，其中規則4對應的區域滑坡危險性最高，主要分布于香溪河右岸、涼臺河及童莊河兩側，鑼鼓洞河兩河口鎮至長江干流段兩側以及長江干流北岸零星分布。

本文方法給出的滑坡危險性評價為基于孕育滑坡環境因素的靜態評價，如引入地震動對滑坡體影響的地震因子和消落區水位變化因子，則可對滑坡危險性進行動態評價。

參考文獻

[1]季偉峰. 三峽庫區三期地質災害防治監測預警的主要方法[J]. 探礦工程(巖土鉆掘工程), 2008(7):14-17(Ji Weifeng. Main Method of Monitoring and Early-Warning for Control of the Third Period of Geological Hazard in Three Gorges Reservoir Region[J]. Exploration Engineering(Rock & Soil Drilling and Tunneling), 2008(7):14-17)

[2]吳永鋒. 三峽庫區地質災害治理若干問題的思考[J]. 人民長江, 2005, 36(3):10-13(Wu Yongfeng. Consideration on Some Issues of Geologic Disaster Control in the Three Gorges Reservoir Area[J]. Yangtze River, 2005, 36(3):10-13)

[3]梁學戰. 三峽庫區水位升降作用下岸坡破壞機制研究[D]. 重慶:重慶交通大學, 2013(Liang Xuezhan. Research on Slope Failure Mechanism under the Action of the Three Gorges Reservoir Water Level Fluctuation[D]. Chongqing:Chongqing Jiaotong University, 2013)

[4]易武，孟召平，易慶林. 三峽庫區滑坡預測理論與方法[M]. 北京: 科學出版社，2011(Yi Wu, Meng Zhaoping, Yi Qinglin. Theory and Method of Landslide Stability Prediction in the Three Gorges Reservoir Area[M]. Beijing:Science Press, 2011)

[5]王志旺. 基于GIS技術的區域滑坡分形特征分析與危險性評價[D].武漢：中國地質大學(武漢)，2010(Wang Zhiwang. GIS-Based Methods for Fractal Analysis and Hazard Estimation of Regional Landslides[D]. Wuhan: China University of Geosciences (Wuhan), 2010)

[6]田陵君.長江三峽工程庫岸穩定性[M].北京:中國科學技術出版社，1992(Tian Lingjun. Bank Stability of Three Gorges Project Reservoir Yangtze River[M]. Beijing: China Science and Technology Press, 1992)

[7]長江水利委員會.三峽工程地質研究[M].武漢:湖北科技出版社，1997(Changjiang Water Resources Commission. Geological Study of the Three Gorges Project[M]. Wuhan: Hubei Science and Technology Press, 1997)

[8]王治華. 三峽水庫區城鎮滑坡分布及發育規律[J]. 中國地質災害與防治學報, 2007, 18(3):33-38(Wang Zhihua. Landslide Distribution and Development in the Towns of the Three Gorges Reservior Area[J]. The Chinese Journal of Geological Hazard and Control, 2007, 18(3):33-38)

[9]王秋良. 北川地區地震災害的地形因子分析[J]. 大地測量與地球動力學, 2008, 28(6):68-72(Wang Qiuliang. Analysis on Terrain Factors of Earthquake Disaster in Beichuan Area[J]. Journal of Geodesy and Geodynamics, 2008, 28(6):68-72)

[10]夏金梧, 李長安. 三峽庫區新構造運動與滑坡耦合關系探討[J]. 人民長江, 2005, 36(3):16-18(Xia Jinwu, Li Changan. Discussion on Coupling Relation of New Tectonic Movement and Landslide in the Three Gorges Reservoir Area[J]. Yangtze River, 2005, 36(3):16-18)

[11]梁學戰, 唐紅梅. 三峽庫區及鄰近地區滑坡發育宏觀地學背景分析[J]. 重慶交通大學學報:自然科學版, 2009, 28(1):100-104(Liang Xuezhan, Tang Hongmei. Analysis on Macro Geography Background of Landslide Development at Three Gorges Reservoir Area and Nearby[J]. Journal of Chongqing Jiaotong University:Natural Science, 2009, 28(1):100-104)

[12]喬建平, 吳彩燕, 田宏嶺. 三峽庫區云陽-巫山段地層因素對滑坡發育的貢獻率研究[J]. 巖石力學與工程學報, 2004, 23(17):2 920-2 924(Qiao Jianping, Wu Caiyan, Tian Hongling. Contribution Rate Research of Stratum to Landslide Growth of Yunyang-Wushan Segment in Three Gorges Reservoir Region[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(17):2 920-2 924)

[13]喬建平, 吳彩燕. 三峽水庫區云陽-巫山段斜坡高差因素對滑坡發育的貢獻率研究[J]. 中國地質災害與防治學報, 2005, 16(4):16-19(Qiao Jianping, Wu Caiyan. The Influence of Relative Elevation on Landslide Occurrence within Yunyang-Wushan Section in the Three Gorges Reservoir Region[J]. The Chinese Journal of Geological Hazard and Control, 2005, 16(4):16-19)

[14]樊曉一, 喬建平, 陳永波. 層次分析法在典型滑坡危險度評價中的應用[J]. 自然災害學報, 2004, 13(1):72-76(Fan Xiaoyi, Qiao Jianping, Chen Yongbo. Application of Analytic Hierarchy Process in Assessment of Typical Landslide Danger Degree [J]. Journal of Natural Disasters, 2004, 13(1):72-76)

Foundation support:Director Fund of the Institute of Seismology, CEA, No. IS201116002, IS201056076；National Natural Science Foundation of China,No.41572354.

About the first author:LI Xue, associate professor, majors in remote sensing information extraction and application of remote sensing for disaster reduction, E-mail：leexue1211@126.com.

The Spatial Distribution Characteristics and Risk Assessment of Landslides in the Head Region of the Three Gorges Reservoir

LIXue1LIJinggang1LIUXiaoli1WANGQiuliang1

1Key Laboratory of Earthquake Geodesy, Institute of Seismology, CEA, 40 Hongshance Road, Wuhan 430071, China

Abstract:The risk of landslides, which have been treated in phase Ⅱ and Ⅲ of the Three Gorges project, are increasing because of the cyclical loading and unloading of the Three Gorges reservoir. The reactivated landslides and new landslides would threaten people’s lives and the safety of property. In order to mitigate the risk, remote sensing (RS) and geographic information systems (GIS) technology are used to study the correlation between the distribution of landslides and other factors, such as topography, lithology, faults, water basins and anthropogenic influences. The results show that the elevation, aspect and elevation variation coefficient in topography have higher correlation with the landslides. Moreover, faults have no obvious effect on the quantity of landslides. The spatial distribution of the landslides in the study area is restricted by the combinations of factors. According to the weights of these restrictions, a risk assessment of landslides is carried out.

Key words:Three Gorges reservoir; remote sensing; correlation; spatial distribution; landslide risk analysis

收稿日期：2015-05-07

第一作者簡介：李雪，副研究員，主要從事遙感信息提取、遙感減災應用研究，E-mail：leexue1211@126.com。

DOI：10.14075/j.jgg.2016.07.016

文章編號：1671-5942(2016)07-0630-05

中圖分類號：P315

文獻標識碼：A

項目來源：中國地震局地震研究所所長基金(IS201116002,IS201056076)；國家自然科學基金(41572354)。