孫 欣，林孝松，何錦峰，韓 賾

（重慶交通大學(xué)河海學(xué)院，重慶400074）

0 引 言

在中國，由于山區(qū)面積比重大、人口多，因此山洪造成的災(zāi)害損失往往十分嚴重，是國內(nèi)自然災(zāi)害造成人員傷亡的主要災(zāi)害種類之一。山洪災(zāi)害具有明顯的多發(fā)性、無序的突發(fā)性和強烈的破壞性，防御難度很大[1]。

本文以重慶市巴南區(qū)安瀾鎮(zhèn)為例，將地理信息系統(tǒng) （geographic information system，GIS）與層次分析法相結(jié)合引入山洪災(zāi)害危險性評價，充分發(fā)揮GIS強大的空間分析能力和層次分析法善于多要素綜合評價的長處，最大限度地將評價過程中的因素量化和權(quán)重判別中的主觀因素與客觀因素有效的結(jié)合在一起，大大提高了山洪災(zāi)害評價工作效率，使山洪災(zāi)害危險性評價結(jié)果更具有科學(xué)性、準確性，評價結(jié)果對山洪災(zāi)害的防治具有重要意義。

安瀾鎮(zhèn)位于巴南區(qū)山洪災(zāi)害的高易發(fā)區(qū)，屬深丘－低山過渡地帶，地勢起伏較大，地形坡度陡，一般地形坡度大于35°，多砂巖質(zhì)陡崖，下伏地層主要為侏羅系遂寧組及沙溪廟組砂巖夾泥巖，在各砂巖陡崖邊緣由于卸荷作用，危巖崩塌較發(fā)育。

安瀾鎮(zhèn)山洪主要是暴雨引發(fā)的滑坡和崩塌，規(guī)模以中小型為主，多處于欠穩(wěn)定狀態(tài)。特別是近10年來，隨著三峽水庫175m蓄、降水工程的實施，每年汛期災(zāi)害更是頻頻發(fā)生，影響著城鎮(zhèn)規(guī)劃建設(shè)及居民的正常生活、工作，威脅交通的暢通，危及農(nóng)戶和耕地，給全鎮(zhèn)經(jīng)濟建設(shè)和社會發(fā)展造成了巨大影響。

1 評價因子選取

山洪災(zāi)害的形成因素眾多，而每一種因子又包含不同的表現(xiàn)形式，使得山洪災(zāi)害系統(tǒng)內(nèi)各因子之間因果表達關(guān)系較為復(fù)雜[2]。參考目前山洪災(zāi)害危險性分析所選取指標內(nèi)容，并結(jié)合鎮(zhèn)域山洪災(zāi)害的特點，山洪危險性評價主要考慮地形、降水、水系、植被、人類活動等5個因子。

（1）地形因子（U1） 地形與山洪危險程度密切相關(guān)。一般認為，地形對形成洪水的影響主要表現(xiàn)在2個方面：地形高程及地形變化程度。山高坡陡地區(qū)，水流湍急，河道狹窄，洪水的影響范圍小。而到了平原地區(qū)，水流漸緩，河道開闊，洪水的影響范圍變大[3]。

（2）降水因子（U2） 降水是暴發(fā)山洪的前提條件，也是山洪致災(zāi)的主要因素。山洪災(zāi)害的嚴重程度受降水強度、歷時和范圍的影響較大，一般降水強度越大、歷時越長、范圍越廣，對山洪形成所起的作用也越大。

（3）水系因子（U3） 水系的空間分布在很大程度上決定了研究區(qū)域遭受洪水侵襲的難易程度，距離河道、水庫等越近，則山洪危險程度越高。河流級別越高，其影響范圍越大，同一級別的河流如果其所處地形不同，影響范圍也會不一樣，平原區(qū)的河流具有更強的影響力。

（4）植被因子（U4） 不同的植被類型對降雨的截留、蓄水能力不同，同時也影響著降雨匯流的時間和大小。不同的植被覆蓋度對山洪的發(fā)生有著重要的影響，植被稀疏的地方，匯流迅速，易形成洪峰，容易暴發(fā)山洪災(zāi)害；植被覆蓋度大的地方，植被減緩降雨匯流，延緩洪峰的形成，降低山洪暴發(fā)的機率。

（5）人類活動因子（U5） 人類不合理的活動是造成山洪災(zāi)害加劇的重要因素之一?，F(xiàn)在人類活動范圍及強度都是空前的，對山洪災(zāi)害的影響也在日益加深。不合理的人類活動對山洪災(zāi)害的影響主要是影響氣候、破壞山體穩(wěn)定、影響河道泄洪等。

結(jié)合安瀾鎮(zhèn)山洪災(zāi)害實際情況，得到安瀾鎮(zhèn)山洪災(zāi)害危險性評價體系，如圖1所示。

圖1 安瀾鎮(zhèn)山洪災(zāi)害危險性評價指標體系圖

2 指標數(shù)據(jù)獲取

根據(jù)研究內(nèi)容和上述評價指標的選取，收集大量相關(guān)資料，利用GIS對相關(guān)資料進行數(shù)字化處理，構(gòu)建評價數(shù)據(jù)庫，并基于安瀾鎮(zhèn)山洪災(zāi)害的特征、分布規(guī)律及專家系統(tǒng)法將每個指標分為4級，其分別賦值為1、2、3、4，具體的分級標準及賦值見表1。

表1 指標分等賦值表

圖2 安瀾鎮(zhèn)坡度圖

圖3 安瀾鎮(zhèn)起伏度圖

圖4 安瀾鎮(zhèn)暴雨量圖

圖5 安瀾鎮(zhèn)年均降水量圖

（1）坡度（C1）和起伏度（C2）數(shù)據(jù)獲取。在1∶1萬的CAD地形圖基礎(chǔ)上，利用MapGIS對等高線進行賦值。在Arc－GIS空間分析模塊支持下，利用Slope函數(shù)提取坡度，由于安瀾鎮(zhèn)面積較小，選取50m×50m柵格大小較為合適，結(jié)果如圖2所示。地形起伏度隨著統(tǒng)計單元半徑的增大，地形起伏度的值也隨之增大，但增大到一定程度即趨于穩(wěn)定。基于安瀾鎮(zhèn)DEM[4]數(shù)據(jù)選取150m×150m作為統(tǒng)計單元，結(jié)果如圖3所示。

（2）暴雨量（C3）和年均降水量（C4）數(shù)據(jù)獲取。收集到重慶市暴雨等值線圖和重慶市以及周邊區(qū)縣共44個氣象臺站多年實測降水量平均值，利用ArcGIS軟件空間插值功能獲取[5]，結(jié)果如圖4所示。年均降水量需要根據(jù)重慶市各個區(qū)縣以及周邊區(qū)縣共44個縣級氣象臺站的高程以及多年平均降水量數(shù)據(jù)，在SPSS中進行線性回歸模擬得到年均降水量的回歸方程[6]，采用反權(quán)重法，得到安瀾鎮(zhèn)的年均降雨量分布模擬數(shù)據(jù)，結(jié)果如圖5所示。

（3）河水匯流量（C5）數(shù)據(jù)獲取。在ArcGIS水文分析模塊支持下，利用生成的DEM數(shù)據(jù)，通過無洼地DEM生成、水流方向、匯流量積累。流量累積值亦代表各網(wǎng)格的上游集流網(wǎng)格數(shù)量，將之乘上單位網(wǎng)格面積便可得到每個網(wǎng)格點的上游集水面積，亦代表每個柵格的匯流能力，得到河水匯流量數(shù)據(jù)[7－8]，結(jié)果如圖6所示。

（4）植被覆蓋率（C6）數(shù)據(jù)獲取。利用收集到的土地利用現(xiàn)狀圖及相關(guān)植被統(tǒng)計數(shù)據(jù)，將小斑塊合并，并對斑塊進行賦值量化，有林地85、竹林70、經(jīng)濟林50、灌木林45、草地30、旱地20、水田15，結(jié)果如圖7所示。

（5）房屋（C7）和道路（C8）數(shù)據(jù)獲取。數(shù)據(jù)來源于土地利用類型，將房屋和道路要素建立緩沖區(qū)，利用距離遠近來表示對山洪災(zāi)害影響的大小，參考專家意見，房屋周圍人類活動頻繁，影響距離為300m，結(jié)果如圖8所示。道路的等級不同，人類活動影響力就有所差異。將道路分為2部分，高速路影響距離為300m，其他道路影響力為1 5 0m，結(jié)果如圖9所示。

圖6 安瀾鎮(zhèn)河水匯流量圖

圖7 安瀾鎮(zhèn)植被覆蓋率圖

圖8 安瀾鎮(zhèn)房屋緩沖區(qū)圖

圖9 安瀾鎮(zhèn)道路緩沖區(qū)圖

3 指標權(quán)重確定

層次分析法 （analytic hierarehy process，AHP）是由美國T.L.Saaty等人20世紀70年代提出的一種定性和定量相結(jié)合的多準則決策方法[9－11]。本文應(yīng)用層次分析法確定權(quán)重系數(shù)，向一些對山洪災(zāi)害有一定研究和認識的專家征詢意見，采用常用的1－9標度方法[12]，綜合構(gòu)造判斷矩陣，獲取8大指標的權(quán)重系數(shù)和最大特征值，對判斷矩陣進行一致性檢驗，得到各個指標權(quán)重 （表2）。

表2 各指標層次總排序權(quán)重表

從影響山洪災(zāi)害危險性因素的總排序來看，主要因子地形、降水、水系、植被、人類活動其權(quán)值分別是：0.308、0.308、0.154、0.154、0.076。

4 危險性評價分區(qū)及評價

4.1 評價指標重分類

利用綜合評價法構(gòu)建了山洪災(zāi)害危險性評價指標體系，層次分析法確定各個指標權(quán)重，山洪災(zāi)害危險性由加權(quán)后的指標分數(shù)決定，分數(shù)越高危險性越高[13]。需要將分級后的各指標進行重分類，重分類的值由表1決定。ArcGIS空間疊加分析功能，使各個指標按照權(quán)重綜合評定安瀾鎮(zhèn)山洪災(zāi)害危險性，得到危險性評價結(jié)果。

4.2 危險性評價分區(qū)圖

目前國內(nèi)外相關(guān)文獻關(guān)于危險度分級標準的問題很少提及，對山洪災(zāi)害危險性評價結(jié)果的危險度分級方法國際上沒有統(tǒng)一的標準[14]。多數(shù)采取的是依照專家個人經(jīng)驗對山洪災(zāi)害危險度評價結(jié)果進行分級。在GIS軟件支持下，全自動對評價結(jié)果進行defined interval（定義間隔分類）、equal interval（等間隔分級）、quantile（等量分類）、standard deviation（標準差分類）、natural breaks（自然斷點分類）、geometry interval（幾何級數(shù)分類）分級劃分對比實驗。實驗表明幾何級數(shù)分類的分級劃分方法分區(qū)效果明顯，更能體現(xiàn)危險性的分布規(guī)律。將結(jié)果分為微危險、低危險、中危險、高危險共4級，得到山洪危險度評價成果圖[15]（圖10），并進行重分類統(tǒng)計，得到山洪災(zāi)害各個等級的危險區(qū)面積 （表3）。

表3 各個等級危險區(qū)網(wǎng)格數(shù)量及面積表

圖10 安瀾鎮(zhèn)山洪災(zāi)害危險性評價分區(qū)圖

5 結(jié)論及防治建議

從最終的危險性評價分區(qū)圖可以看出：安瀾鎮(zhèn)山洪災(zāi)害危險性在空間上面積較廣，等級較高。鎮(zhèn)內(nèi)山洪災(zāi)害中危險區(qū)和高危險區(qū)主要集中在沿河地帶，從海拔高度上具有相對集中性。中危險區(qū)和高危險區(qū)主要分布在安瀾鎮(zhèn)西北方，與地形走向東南高、西北低相互吻合，呈反向關(guān)系。研究結(jié)果與實際情況吻合較好，可以將研究區(qū)山洪災(zāi)害危險性評價分區(qū)圖用于該區(qū)山洪災(zāi)害防治工作：

（1）健全和完善有關(guān)法律法規(guī)，特別要保證山洪災(zāi)害重點防治區(qū)內(nèi)退耕還林和移民搬遷生態(tài)環(huán)境保護等方面的政策、法律、法規(guī)得到嚴格執(zhí)行。

（2）根據(jù)山洪災(zāi)害危險評價圖，編制山洪災(zāi)害防治預(yù)案，建立山洪災(zāi)害預(yù)防領(lǐng)導(dǎo)、指揮及組織機構(gòu)，進行山洪災(zāi)害普查，明確山洪災(zāi)害范圍與影響程度，確定避災(zāi)預(yù)警程序和臨時轉(zhuǎn)移人口的路線和地點。

（3）根據(jù)安瀾鎮(zhèn)山洪災(zāi)害危險評價分區(qū)圖，調(diào)整山洪災(zāi)高危險區(qū)土地利用結(jié)構(gòu)。

（4）加強分區(qū)圖所示危險性高的河道管理力度，控制水土流失，嚴格禁止侵占泄洪河道行為，疏通洪水、宣泄渠道。

[1]傅春梅，徐剛.重慶市山洪災(zāi)害的危害、成因及防治對策[J].太原師范學(xué)院學(xué)報：自然科學(xué)版，2011，10（01）：145－148.

[2]劉廷祥，吳泉源，鮑文東，等.GIS技術(shù)支持的黃水河流域山洪災(zāi)害風(fēng)險性分析[J].山東師范大學(xué)學(xué)報：自然科學(xué)版，2009，24（02）：92－95.

[3]李紅霞，楊志，高宇，等.基于GIS的黃土滑坡危險性評價的指標體系的建立[J].寧夏工程技術(shù)，2007，6（01）：70－72.

[4]陸忠艷，袁子鵬，陳艷秋，等.基于DEM 的遼寧鄉(xiāng)鎮(zhèn)漬澇預(yù)報[J].中國農(nóng)業(yè)氣象，2008，29（04）：463－466.

[5]林孝松，余情.重慶市農(nóng)業(yè)氣候資源優(yōu)度評價[J].中國農(nóng)業(yè)氣象，2007，28（03）：241－245.

[6]余情，郭躍，林孝松，等.基于模糊物元的重慶市農(nóng)業(yè)氣候資源優(yōu)度評價[J].海南師范大學(xué)學(xué)報：自然科學(xué)版，2008，21（02）：208－212.

[7]莫奇京，羅蔚生，楊桂賢，等.基于ASTER GDEM的石漠化綜合治理中的小流域提取研究——以廣西石漠化綜合治理應(yīng)用為例[J].林業(yè)調(diào)查規(guī)劃，2011，36（05）：22－25，35.

[8]唐余學(xué)，廖向花，李晶，等.基于 GIS的重慶市山洪災(zāi)害區(qū)劃[J].氣象科技，2011，39（04）：423－428.

[9]Tumpal P T，Sinaga，A N，Lee Y W，et al.GIS mapping of Tsunami vulnerability：a case study of the Jembrana regency in Bali，Indonesia[J].KSCE J Civil Engin，2011，15（03）：537－543.

[10]唐紅梅，陳洪凱，唐芬，等.重慶庫區(qū)松散土體滑坡危險性分區(qū)及評價[J].災(zāi)害學(xué)，2004，19（01）：1－6.

[11]Bathrellos G D，Gaki－Papanastassiou K，Skilodimou H D，et al.Potential suitability for urban planning and industry development using natural hazard maps and geological geomorphologic parameters[J].Envir Earth Sci，2012，66（02）：537－548.

[12]李景保，喻小紅，金濤.湖南四水流域洪澇災(zāi)害特性與減災(zāi)戰(zhàn)略[J].水土保持通報，2002，22（05）：57－60.

[13]Kazuyuki O，Hiroshi S，Tadahiro K，et al.Characterization of debris flows by rainstorm condition at a torrent on the Mount Yakedake volcano，Japan[J].Geomorphology，2012，136（01）：88－94.

[14]楊志雙，韓玉龍，張浩然，等.基于可拓理論的西南典型山區(qū)泥石流的危險性分級[J].河南理工大學(xué)學(xué)報：自然科學(xué)版，2013，32（01）：35－39.

[15]李娜，王靜.山洪災(zāi)害特點分析及風(fēng)險圖制作方法探討[J].水利水電技術(shù)，2011，42（02）：57－61.