山地旅游城市滑坡災害風險性評價
——以重慶市為例

馬 煒，陳麗聰，冉啟香，劉增力，蔣亞芳

(1.國家林業局調查規劃設計院，北京 100714；2.中林天合森林認證中心有限公司，北京 100013；3.重慶市巫溪縣林業局，重慶 405800)

1 引言

滑坡是主要的地質災害類型[1]，2016年全國共發生地質災害9710起，其中滑坡7403起，對地區經濟與旅游造成重大的影響，與2015年同期相比，地質災害數量、造成的經濟損失、人員傷亡都有所增加[2]。旅游業是重慶市的經濟可持續發展的重要助力與依托[3]，然而重慶市60%以上的旅游資源均位于重大地質災害的多發地段[4]，旅游旺季4～10月同時也是地質災害高發的季節[5]，重慶市作為西南山區典型城市，氣候多變、地質活動強烈、地形地貌復雜，且處于長江及主干流流域，是地質災害滑坡的高發區[6]。加之重慶市城市建設日新月異，城鎮化水平不斷提升，導致生態壓力不斷增大[7]，建設中稍有不慎容易導致山地環境破壞，加劇滑坡災害的發生。科學有效的系統開展滑坡災害危險性評測有利于保障旅游地居民與游客安全，加強旅游地地質災害防范，保證旅游的可持續發展[8]。

滑坡災害相關研究隨著信息技術的發展而不斷深入[9]。滑坡災害風險評測是理清滑坡誘發因素基礎上，運用科學的方法，建立滑坡風險評價機制的過程[10]。滑坡災害的誘發因素多樣，其中降水與地震為最主要的因素，而滑坡災害發生后其危害性受其體量與所處地理位置的影響[11]。因此科學的評測滑坡災害的危險性需要建立合理的滑坡災害評測模型。國內關于滑坡災害危險性評價模型可分為定性分析模型與定量分析模型，如主觀臆斷模型、統計分析模型等[12]。隨著信息技術的發展，滑坡危險性評價方法不斷豐富，有統計分析法、信息熵法、突變級數法、AHP因子加權法、模糊綜合評價法、人工神經網絡方法等[13～15]，但定性與定量綜合使用成為主流研究方法[10]。重慶市作為滑坡高危區，已有學者運用GIS技術基于重慶市已有滑坡數據進行災害時空分布、重點區域三峽庫區的萬州區、巫山縣，北碚區、涪陵區的滑坡情況進行分析，從宏觀上得出重慶市的滑坡誘發因素主要有降水、人類活動、自然災害、水庫水的漲落，滑坡密度高發區集中在主城區及忠縣、萬州區等為核心的三峽庫區[16～20]，基于已有基礎數據的滑坡風險預測模型構建還缺乏系統的研究，同時對于重慶市各縣市的滑坡危險性分級評測未有涉及。

基于此，本文依據重慶市已有情況，從滑坡的發生發展的誘發因素及影響過程，選取10個指標，運用ArcGIS空間分析技術，采用定量與定性相結合的方法，構建由地理環境因素、地質構造因素、人類活動因素組成的滑坡災害危險性評測體系，以期運用基礎數據構建滑坡災害風險等級評價模型，對重慶市各個縣的滑坡災害危險性進行評價，為重慶市旅游地地質災害防范提供理論基礎。

2 研究區概況與數據來源

2.1 研究區概況

重慶市地處我國西南部，是我國直轄市、國家中心城市，跨東經105°11′～110°11′、北緯28°10′～32°13′，是我國著名的山城，山地面積占76%，是國家歷史文化名城，旅游資源豐富，擁有世界文化遺產1處，國家5A級景區7個，4A級景區達76個，2016年全市接待境內外游客達到4.51億人次，實現旅游總收入2645.21億元，居2016年中國旅游城市吸引力排行榜第三，僅次于上海和北京。然而，重慶市2016年地質災害發生數量排名全國前十，由地質災害造成的經濟損失、人員傷亡失蹤也排名全國前十，1950～2011年間發生滑坡災害16500多起，是西南地區受地質災害影響的典型區域。

2.2 數據來源

本文關于重慶市山體滑坡災害危險性評價的10個指標，其中高度與坡度數據是根據中國科學院計算機網絡信息中心(地理空間數據云)的GDEM數據集中的2015年重慶市30 m DEM數據提取獲得；多年平均降水量數據為中國氣象科學數據共享網上的重慶市11個區域內和14個區域外的站點的多年平均降水量數據進行克里金插值處理得到；NDVI數據源自地理空間數據云平臺的MOD13Q1數據集，預處理過程中對2015全年的NDVI 產品求年平均，再進行重分類，得到研究區重慶

市的NDVI指標數據；巖性數據參考1：2500000中國地質圖，進行矢量化處理，裁剪出重慶市巖性數據；城鎮、河流、人口密度等基礎數據是2015年遙感影像解譯數據，影像源自國際科學數據服務平臺；地震點、斷裂帶、滑坡數據來自重慶市地質災害防治相關圖件矢量化。

3 研究方法

3.1 評價體系構建

在綜合其他學者相關研究成果的基礎上，本著代表性、可操作性為原則，本文根據滑坡發生發展的誘發因素及影響過程選取山地滑坡災害相關的10個具體指標，分別為坡度、高程、巖性、多年平均降水量、歸一化植被指數(NDVI)、地震密度、距河流距離、斷裂帶距離、人口密度、距城鎮距離，將10個指標的原始數據，通過插值工具、核密度工具、緩沖區工具將原始數據轉換為圖形數據，并重分類為山地滑坡災害相關的三類因素，分別為地理環境因素、地質構造因素及人類活動因素，用來綜合評價重慶市各縣市滑坡災害危險性等級(圖1)。

圖1 重慶市滑坡危險性等級評價體系構建

3.1.1 地理環境因素

地理環境因素對滑坡的影響是直接的、深刻的，是形成山地滑坡災害的基礎和根本，在山地滑坡災害研究中具有不可代替的作用。本文選取坡度、高程、NDVI 、巖性、距河流距離及多年平均降水量6個指標作為滑坡發生的影響因子：①坡度是滑坡的前提條件之一，以往研究表明坡度與滑坡之間存在一定的關系，坡度越高，對滑坡影響越大；②高程的大小對滑坡體具有的勢能大小有著直接影響，進而影響到滑坡的危險程度，高程越大，滑坡的危險度越高；③歸一化植被指數能夠很好地反映地表覆被情況，而植被能夠保水固土，植被覆蓋及地表水的搬運、侵蝕會影響岸坡穩定性，影響滑坡的發生發展，NDVI指數越高，滑坡概率越小；④巖石結構構造決定斜坡巖牢固性、可破壞性，為滑坡發生、發展提供物質基礎，本文參考相關研究，將地質圖矢量化并數提取研究區巖石類型信息，并分為5個等級巖組(Ⅰ.砂巖、砂泥夾白云巖等巖組，Ⅱ.泥砂巖、互層巖巖組，Ⅲ.變質巖、漿巖巖組，Ⅳ.碳酸鹽巖巖組，Ⅴ.黏土、砂卵礫石巖組)，等級越高，在外力作用下越容易發生滑坡；⑤河流是滑坡體的最終歸集地，距離河流長短決定滑坡體的影響范圍，同時滑坡體進入河流，對河流產生一系列影響，進而增加滑坡災害的危險性，距離河流越近，滑坡影響越大；⑥降水是引起山地滑坡發生的主要外界因子，滑坡的發生數量、規模與持續過程降雨及暴雨量等的關系十分明顯，降水入滲會減小土體的抗剪度及土體與基巖的摩擦阻力，并增加土體重度和內部的動水壓力，從而誘發地質災害，降水量越大，滑坡發生概率越大。

3.1.2 地質構造因素

地質構造因素是山地滑坡發生的內營力，特別是地震和斷裂帶對山地滑坡產生具有重要的影響，因此選用地震點核密度和距斷裂帶距離2個指標。①地震產生的地震力直接作用于斜坡巖土導致滑坡災害的發生，且脆弱山坡中積累的損害可能會受該地區先前地震遺留的影響，地震點核密度代指地震帶分布的空間集聚狀況，地震點核密度越高，滑坡越容易發生；②中國的大陸地震主要受活動構造的控制，易發生地震災害的環境往往是地質構造極為復雜、斷裂發育、巖石破碎，從而間接影響滑坡災害的發生，滑坡災害發生受距離活動斷裂帶距離的控制，距離斷裂帶越近，滑坡發生可能性越高。

3.1.3 人類活動因素

人類工程活動和經濟建設對自然環境進行破壞一定程度上造成滑坡災害的可能性，基于可行性原則，本文選用距城鎮距離及人口密度來作為滑坡災害性的評價因素。①人口密度一定程度上反映區域人口活動的強度，人口活動強度對區域地形地貌、植被環境都具有一定影響，也就間接影響到山地滑坡災害，因此本文認為人口密度越高，滑坡災害性等級越高；②距離城鎮的距離影響到滑坡對城鎮潛在的危險程度，也能在一定程度上反映山地滑坡危險性，距離越近，災害等級越高。

3.2 指標權重確定

每個評價因子對滑坡的影響是不同的，合理的確定每個滑坡因子的貢獻率才能確保評價的準確性。本文引用已有學者相關研究成果[21]，并運用AHP法確定指標權重，如表1所示。

表1 滑坡災害危險性評價因子權重賦值表

3.3 評價模型

GIS疊加分析是對某一地區的相同比例尺的多層矢量數據進行空間邏輯上的運算，以獲得對某一地區的綜合分析，來滿足用戶的需求。在矢量柵格數據內部，像元之間的各種運算來實現疊加運算，而進行疊加運算時的各層面數據是運用數學關系建立的，通過疊加運算來發現某地新的空間規律或空間變化。設兩個圖層Xij(i=1，m；j=1，n)，Yij(i=1，m；j=1，n)，Zij=(i=1，m；j=1，n)，則矢量柵格的邏輯運算規則為

本文首先利用歐幾里德度量即歐式距離法來測算兩點間的真實距離，即測算滑坡危險性評級的距離河流距離、距離斷裂帶距離、距離城鎮距離3個指標，其二維空間表達式為

(1)

式(1)表示X與Y的歐式距離，其中X1，Y1代表X點的橫坐標和縱坐標，X2，Y2表示Y的橫坐標和縱坐標，并根據測算結果建立數據庫；同時運用經驗貝葉斯克里金插值法計算出各地的年平均降水量以及人口密度指標得到兩個因子的數據庫；運用核密度方法得到地震點核密度數據庫；并運用GIS其余空間分析方法提取高程、巖性、坡度、NDVI數據建立數據庫。根據確定好的各評價因子權重，利用已經建好的10個指標數據庫進行GIS疊加分析，文中運用Arcgis10.2軟件中的加權疊加方法進行矢量柵格疊加分析，加權疊加公式為

A=∑Bn×Wn，

(2)

式(2)中A為疊加結果，Bn代表第n個指標的矢量柵格分析圖層，Wn為第n個指標的權重值。

根據疊加法則，設y1為地理環境因素，y2為地質構造因素，y3為人類活動因素，x1～x6為y1的6個因子，x7～x8為y2的2個因子，x9～x10為y3的2個因子，f函數表示各層上屬性與用戶需求之間的關系，則最終滑坡危險性評價加權疊加用公式可表示為

C=f(y1，y2，y3)，

(3)

式(3)中y1=f(x1，…x6)，y2=f(x7，x8)，y3=f(x9，x10)，C為疊加結果。

4 結果與分析

4.1 評價因子分析

以重慶市基礎數據將滑坡災害危險性評價的10個指標數據庫進行信息量提取并按照自然斷點法進行分類，將分類后的數據歸一化，并按照單一因素滑坡災害誘發危險性概率重分類為5個等級，等級越高，對滑坡危險性貢獻越大，制作10個指標不同等級空間分異圖(圖2)，并按照自然環境因素、地質構造因素、人類活動因素進行加權疊加得到不同等級分布圖(圖3)。

由圖2可知，10個指標距河流距離與距城鎮距離等級全市較高空間分異不明顯；坡度與高程較為接近，東北部大巴山山脈、東南部與南部有巫山、大婁山等山脈，整體東部高于等級高于西部；年平均降水東南部、東北部等級最高；NDVI等級西部高于東部，最高區域為長江流域；巖性等級分布于坡度與高程相反；地震點核密度等級最高區域分布在巫山縣、纂江縣、南川縣以及榮昌縣；距斷裂帶距離整個市域等級較高，等級最高區域在大巴山、巫山等褶皺帶上；人口密度等級自西向東遞減。

如圖3所示，重慶市滑坡危險性評價地理環境因素危險性等級整體較高，達3.13，在全市中地理環境因素評價危險性較高的區域為開縣、云陽縣、萬州區、忠縣、梁平縣、豐都縣、涪陵區等長江主干流流域。地質構造因素危險性評價全市空間分異較大，平均等級達2.7，滑坡危險性等級較高的區域主要集中在斷裂帶周邊，城口縣、巫溪縣、巫山縣、奉節縣東部區域組成的大巴山斷褶帶，萬州區、石柱土家族自治縣、豐都縣、武隆縣、北碚區、璧山縣、永川區所處斷裂帶都屬于高風險區域。人類活動滑坡危險性評價全市整體等級較高，西部區域尤為明顯。

圖2 重慶市滑坡災害危險性10個指標等級空間分異

圖3 重慶市滑坡災害三類因素危險性等級空間分異

4.2 重慶市滑坡災害危險性評價

將通過加權疊加處理好的重慶市滑坡風險地理環境因素、地質構造因素及人類活動因素等級數據再一次進行GIS加權疊加運算，并按照自然斷點法進行分級，得到最終的重慶市滑坡災害危險性等級分布圖，等級越高，滑坡災害危險性越大(圖4)，同時提取各等級面積信息，得到各等級面積所占百分比分別為1級滑坡風險地區占總面積的10.24%，2級占27.2%，3級占33.3%，4級占20.7%，5級占總面積的8.55%。可知，重慶市全市滑坡風險災害等級3級、4級占比非常大，滑坡災害風險等級較高的區域主要分布在長江三峽庫區的開縣、云陽縣、萬州區、忠縣等地區，同時長江主干流流域如嘉陵江流域北碚區高滑坡風險等級較高，此外東北部沿大巴山山麓滑坡災害風險等級也較高。就整體區域而言東部滑坡災害危險性大于西部，北部地區大于南部，評

圖4 重慶市滑坡災害危險性等級空間分異

價結果空間分布與地理環境因素等級評價結果的相似度較高。

為進一步明晰重慶市滑坡災害危險性的詳細分布，提取各個縣的滑坡災害風險等級面積數據比上每個縣的總面積，得到重慶市各個縣滑坡災害危險性各等級所占百分比，同時對5個滑坡災害風險等級分別賦值為1-5，等級越高，分值越高，將各個縣每個滑坡風險等級所占面積比乘以其等級所附分值再求和，得到各個縣滑坡危險性分數如表2所示。由表3可知，位于東部的長江三峽腹心、三峽庫區淹沒縣石柱土家族自治縣5級滑坡災害風險面積占51%，居重慶市滑坡災害危險性最高，石柱土家族自治縣海拔較高，且雨水充沛，處于渝東褶皺地帶與長江上游，滑坡危險性大，但是其NDVI等級高，即森林覆蓋率較高，一定程度上降低滑坡災害的擴大；其次是位于重慶西北部的主城區之一的北碚區居滑坡災害危險性居第2位，北碚區4級以上滑坡災害危險性區域占到總區域面積的73%，該人口密度大，嘉陵江穿流而過，巖組較容易發生突變，增加了滑坡災害的危險性；位于三峽庫區的萬州區、云陽縣、忠縣、豐都縣、開縣，滑坡災害危險性得分較高，在全市危險性排名居前10，其中萬州區4級以上的滑坡災害危險性面積占總區域面積的70.9%，排名第3；位于重慶西部、大都市區內的璧山縣與渝北區滑坡災害危險性3級以上的區域占比分別超90%與80%，危險性排名居全市第5和第7；城口縣位于川、渝、陜三省交界處，處于大巴山弧形斷褶帶上，海拔較高，滑坡災害危險性3級以上占比達77%，居全市第9；位于西北角的潼南縣為盆地淺山地區，遠離地震與斷裂帶，滑坡災害危險性全市最小。

表2 重慶市各縣市滑坡災害危險性等級分布詳細表

續表2

縣/區名稱滑坡災害危險性各等級所占面積百分比1級(%)2級(%)3級(%)4級(%)5級(%)危險性得分排名巴南縣12.5331.5643.6311.690.592.56324長壽區8.1045.8738.057.150.832.46725榮昌縣0.5655.0342.511.780.122.45926綦江縣20.6531.1336.0611.230.922.40627大足縣3.3660.1131.214.710.622.39128巫山縣18.8737.1834.618.680.672.35129墊江縣9.9863.3719.925.890.832.24230彭水苗族土家族自治縣35.0340.5318.524.731.201.96531潼南縣24.7865.069.750.420.00081.85832

4.3 重慶市滑坡災害危險性評價模型驗證

將重慶市多年等級以上的地質災害圖件進行矢量化，并提取圖件中的滑坡點，與評價出的滑坡災害風險等級矢量圖進行疊加，來驗證評價結果，如圖5所示，整體擬合度較好，根據本文模型評價出的中高等級風險區，滑坡較為密集，城口縣、開縣、云陽縣、萬州區表現尤為突出。但是在重慶東北部巫山縣存在一個密集滑坡點，由于巫山縣處于長江三峽庫區腹心，跨長江巫峽兩岸，年降水量分布不均，集中在夏季且多暴雨，同時受人類開礦等因素影響形成滑坡聚集區，其滑坡多發生在長江及主要支流一帶，海拔在400 m以下[22]，在本文評價模型中，巫山縣3級以上較多，但在全市排名較為靠后，后續應加強模型實測能力。

圖5 重慶市滑坡災害危險性等級預測與實際滑坡點的擬合

5 結論與討論

(1)重慶市滑坡危險性評價地理環境因素危險性等級整體較高，達3.13，危險性較高的區域為開縣、云陽縣、萬州區、忠縣、梁平縣、豐都縣、涪陵區等長江主干流流域。地質構造因素危險性評價全市空間分異較大，平均等級達2.7，滑坡危險性等級較高的區域主要集中在斷裂帶周邊。人類活動滑坡危險性評價全市整體等級較高，西部區域尤為明顯。

(2)重慶市滑坡災害風險等級較高區域面積較大，5級占區域面積的8.55%，4級占20.7%，3級占33.3%，2級占27.2%，1級占10.24%，就整體區域而言東部滑坡災害危險性大于西部，北部地區大于南部，三峽庫區及長江干流區域及大巴山脈區域滑坡災害危險性等級較高，評價結果空間分布與地理環境因素等級評價結果的較為接近。

(3)重慶市各縣區滑坡災害危險性等級排名前3的縣市分別是三峽庫區的石柱土家族自治縣與萬州區，以及西北部主城區北碚區，其余三峽庫區云陽縣、忠縣、豐都縣、開縣，以及西北部大都市區內的璧山縣與渝北區都在前十，危險性等級較高。位于西北角的潼南縣為

本文設計的重慶市滑坡災害危險性等級評價模型在全市區域驗證性較好，三峽庫區表現最為明顯，但是同樣為三峽庫區的巫山縣受地理位置、降水因素的影響導致成為滑坡實際高危區，本模型由于坡度與高程的設置導致模型驗證性不足。導致滑坡的因素非常多，如何能夠根據區域特殊性構建合理的滑坡災害評價模型應予以進一步的思考。

同時，滑坡影響旅游基礎設施尤其是旅游公路的暢通，對旅游安全造成一定的影響，本文研究的滑坡高危區域為三峽庫區及主城區，而三峽庫區與主城區又為重慶市旅游重要區域，地質災害對旅游可持續發展的影響、旅游地如何提高自身基礎設施以提高地質災害的應對能力都值得進一步的思考。