云南省大關(guān)縣地質(zhì)災(zāi)害風(fēng)險性評價及分區(qū)研究

和子祺,李波,蔡映坤,劉志明,馬云燕

(昆明理工大學(xué)國土資源工程學(xué)院,昆明 650093)

大關(guān)縣作為云南省省級劃定的水富-大關(guān)-巧家-東川地質(zhì)災(zāi)害危險區(qū)、地質(zhì)活動斷裂帶、滇東北地震活躍區(qū),生態(tài)環(huán)境十分脆弱。受區(qū)域地形地貌的影響,人類活動主要集中在深切河谷及山前堆積區(qū),特殊的高原山地地區(qū)加之不合理的人類活動及降雨等誘發(fā)因素的影響,致使區(qū)內(nèi)地質(zhì)災(zāi)害類型多樣且各類型地質(zhì)災(zāi)害密集發(fā)育,是云南省遭受地質(zhì)災(zāi)害最為嚴(yán)重的縣市之一[1-6]。

而為有效地預(yù)防地質(zhì)災(zāi)害,需要首先完成地質(zhì)災(zāi)害的易發(fā)性評價,還需要從區(qū)域上進(jìn)行地質(zhì)災(zāi)害危險性評價,對區(qū)域的社會屬性進(jìn)行深入了解,進(jìn)行研究區(qū)的易損性研究。在易發(fā)性評價、危險性評價、易損性評價基礎(chǔ)上完成地質(zhì)災(zāi)害風(fēng)險評價,反應(yīng)研究區(qū)地質(zhì)災(zāi)害總體風(fēng)險水平。

自2015年大關(guān)縣完成1∶5萬地質(zhì)災(zāi)害詳細(xì)調(diào)查以來,指導(dǎo)地質(zhì)災(zāi)害防治的區(qū)劃以及相關(guān)建議嚴(yán)重滯后,與現(xiàn)狀嚴(yán)重不符,影響當(dāng)?shù)厝嗣竦纳敭a(chǎn)安全及經(jīng)濟(jì)的發(fā)展。因此,亟待進(jìn)行新一次的地質(zhì)災(zāi)害調(diào)查與風(fēng)險評價工作,在以往已有資料的基礎(chǔ)上進(jìn)一步提高精度,充分利用新理論調(diào)查、研判、評價地質(zhì)災(zāi)害風(fēng)險,為地質(zhì)災(zāi)害防治工作提供依據(jù)。本文結(jié)合現(xiàn)有的地質(zhì)災(zāi)害數(shù)據(jù)資料,以ArcGIS為平臺,對大關(guān)縣進(jìn)行地質(zhì)災(zāi)害風(fēng)險性評價,編制地質(zhì)災(zāi)害風(fēng)險區(qū)劃圖,對該區(qū)地質(zhì)災(zāi)害防治和地質(zhì)災(zāi)害風(fēng)險管控提供依據(jù),并以此研究成果為大關(guān)縣地質(zhì)災(zāi)害防治“十四五”規(guī)劃提供支撐。

1 研究區(qū)概況以及地質(zhì)災(zāi)害情況

大關(guān)縣位于云南省東北部,昭通地區(qū)的腹心地帶?？h境氣溫平面差異不大,垂直差異突出,為山區(qū)立體氣候,全年平均氣溫15.0℃, 最高40.3℃,最低-6.4℃?？h境內(nèi)降水分布不均衡,年平均降水量992.9 mm[7]。整個縣域地貌以峽谷山岳地形地貌為主,山高坡陡,溝谷縱橫,地勢南北高,中間低,相對高差大。大關(guān)縣境內(nèi)除太古界、元古界和古生界、中生界的地層個別缺失外,其余地層均有出露。地層出露巖性主要為淺海相碳酸鹽巖,濱海相及海陸交互相的砂巖、粉砂巖,巖漿巖主要為二疊系玄武巖。溝谷和斜坡地帶廣布第四系松散堆積層。大關(guān)處于揚(yáng)子準(zhǔn)地臺邊緣,川滇經(jīng)向構(gòu)造帶北段。區(qū)內(nèi)以褶皺和斷裂構(gòu)造為主,總體為北東向。境內(nèi)構(gòu)造形跡北部和南部主要以南北向或北西向構(gòu)造為主,中部構(gòu)造形跡主要為北東向。地層褶皺寬闊,多數(shù)以短軸背斜出現(xiàn),按展布方向可分為南北、北東、北西和東西向構(gòu)造?？h域內(nèi)新構(gòu)造運(yùn)動強(qiáng)烈,近代地震活動強(qiáng)烈,河谷下切,河流深切和側(cè)蝕作用加劇[8]。根據(jù)最新云南省的地質(zhì)災(zāi)害調(diào)查工作的開展,區(qū)內(nèi)共發(fā)育地質(zhì)災(zāi)害點(diǎn)567處,包括滑坡299處、崩塌210處、泥石流56條及地面塌陷2處,多沿區(qū)內(nèi)各主要河道及其支流谷坡區(qū)集中發(fā)育。其中滑坡以土質(zhì)滑坡為主,且多為淺層、推移式、小型滑坡;崩塌以巖質(zhì)崩塌為主,且多為墜落式、小型崩塌;泥石流以溝谷型、低-中頻、發(fā)育期、中-小型的泥石流為主,水源類型均為暴雨型,按物質(zhì)組成多為泥石流型;2處地面塌陷特征不明顯,僅局部地表出現(xiàn)開裂等變形跡象。

2 大關(guān)縣地質(zhì)災(zāi)害風(fēng)險性評價及分區(qū)

地質(zhì)災(zāi)害風(fēng)險性評價及分區(qū)是一項系統(tǒng)的、綜合的評價手段,包括了地質(zhì)災(zāi)害易發(fā)性評價、地質(zhì)災(zāi)害危險性評價以及承災(zāi)體易損性評價的內(nèi)容,對地質(zhì)災(zāi)害的防治具有重要的理論和實際意義。通過對云南省大關(guān)縣的野外地質(zhì)災(zāi)害調(diào)查和相關(guān)資料的收集,建立了該區(qū)的風(fēng)險性評價模型,即風(fēng)險性=危險性×易損性。評價步驟首先對研究區(qū)進(jìn)行了易發(fā)性評價,采用柵格單元作為評價單元,由于進(jìn)行全縣域地質(zhì)災(zāi)害風(fēng)險評價影響因子較多,不同影響因子間相互聯(lián)系、相互制約,通過對各常用評價方法進(jìn)行比較,本文采用基于層次分析法-信息量法模型進(jìn)行;其次,以10 a平均降雨量為誘發(fā)因素,再利用信息量法對其進(jìn)行信息量值計算,結(jié)合地質(zhì)災(zāi)害易發(fā)性指數(shù)對研究區(qū)進(jìn)行了地質(zhì)災(zāi)害危險性評價;第三,以收集到的最新大關(guān)縣年鑒及第三次全國土地普查數(shù)據(jù)為基礎(chǔ),基于承災(zāi)體易損性各屬性特點(diǎn),選擇人口密度、道路密度、耕地密度及GDP密度作為本次地質(zhì)災(zāi)害易損性評價指標(biāo)對研究區(qū)進(jìn)行地質(zhì)災(zāi)害易損性評價;最后依據(jù)完成的危險性評價、易損性評價,利用ARCGIS軟件的疊加分析功能,完成大關(guān)縣的地質(zhì)災(zāi)害風(fēng)險性區(qū)劃。

2.1 大關(guān)縣地質(zhì)災(zāi)害易發(fā)性評價及分區(qū)

2.1.1 評價方法

(1) 信息量法

主要通過各指標(biāo)區(qū)間已知災(zāi)點(diǎn)信息量值的計算,表征各指標(biāo)區(qū)間對發(fā)生該類型地質(zhì)災(zāi)害的貢獻(xiàn)值,最終進(jìn)行區(qū)域地質(zhì)災(zāi)害易發(fā)性的評價。對某種因素特定狀態(tài)下地質(zhì)災(zāi)害信息量計算公式為:

(1)

式中,Ixi→A為對應(yīng)因素x在i狀態(tài)或區(qū)間條件下,地質(zhì)災(zāi)害A發(fā)生的信息量值;Ni為對應(yīng)因素x在i狀態(tài)或區(qū)間條件下,地質(zhì)災(zāi)害A面積或點(diǎn)數(shù);N為調(diào)查區(qū)地質(zhì)災(zāi)害A總面積或總點(diǎn)數(shù);Si為對應(yīng)因素x在i狀態(tài)或區(qū)間條件下的分布面積;S為調(diào)查區(qū)總面積。

當(dāng)Ixi→A>0時,反映了地質(zhì)災(zāi)害A在對應(yīng)因素x在i狀態(tài)或區(qū)間下發(fā)生地質(zhì)災(zāi)害的可能性較大,或者說該狀態(tài)或區(qū)間有利于地質(zhì)災(zāi)害A的發(fā)生;當(dāng)Ixi→A<0時,則反之說明不利于地質(zhì)災(zāi)害A的發(fā)生;當(dāng)Ixi→A=0時,反映了不提供有關(guān)地質(zhì)災(zāi)害A發(fā)生與否的任何信息[9]。

每個評價單元地質(zhì)災(zāi)害的易發(fā)性是各因素綜合影響的結(jié)果,各影響因素又存在不同的狀態(tài)或區(qū)間,各因素各狀態(tài)區(qū)間地質(zhì)災(zāi)害發(fā)生的總信息量值計算公式為:

(2)

式中,I為每個評價單元地質(zhì)災(zāi)害發(fā)生的總信息量值;Ni為對應(yīng)因素在i狀態(tài)或區(qū)間條件下地質(zhì)災(zāi)害的面積或地質(zhì)災(zāi)害的數(shù)量;Si為對應(yīng)因素在i狀態(tài)或區(qū)間分布的總面積;N為調(diào)查區(qū)地質(zhì)災(zāi)害的總面積或總地質(zhì)災(zāi)害數(shù)量;S為調(diào)查區(qū)總面積。

(2) 層次分析法

層次分析法是一種將復(fù)雜的問題轉(zhuǎn)化為有序的階梯層次結(jié)構(gòu)的多準(zhǔn)則決策方法,具有實用性、系統(tǒng)性、簡潔性。在對各影響因素與地質(zhì)災(zāi)害及各影響因素之間的相關(guān)性進(jìn)行深入分析的基礎(chǔ)上,利用層次分析法能有效地將各影響因素的權(quán)重進(jìn)行數(shù)字化[10]。

(3) 易發(fā)性指數(shù)

信息量模型未對各評價因子之間的關(guān)系及其權(quán)重進(jìn)行分析考慮,本次研究為了提高評價的真實性、客觀性,將各次級因子信息量值(Ii)與其對應(yīng)因子的權(quán)重值(ωi)計算出易發(fā)性指數(shù)(Yi)(見式3),最終通過易發(fā)性指數(shù)進(jìn)行地質(zhì)災(zāi)害易發(fā)性的分析評價。

Yi=Ii×ωi

(3)

2.1.2 評價流程

(1) 評價單元的確定

常用的評價單元類型主要包括行政單元、柵格單元、均一條件單元及斜坡單元等?？紤]評價精度問題,選擇了柵格單元作為本次研究的評價單元?；谇捌趯W(xué)者的經(jīng)驗公式[11],1∶50 000的調(diào)查評價工作柵格單元的大小計算結(jié)果為32.85 m。

Gs=7.49+0.000 6s-2.0×10-9s2

+2.9×10-15s3

(4)

式中,Gs為經(jīng)驗網(wǎng)格大小;s為原始等高線數(shù)據(jù)精度的分母(50 000)。

利用ARCGIS的矢量柵格轉(zhuǎn)化功能以DEM數(shù)據(jù)為基礎(chǔ),進(jìn)行柵格單元的劃分,共劃分出24 319 122個柵格單元,按柵格單元進(jìn)行各影響因素數(shù)據(jù)的提取。

(2) 評價因子的確定

① 因子選取

本次評價基于各孕災(zāi)地質(zhì)條件與各類地質(zhì)災(zāi)害的相關(guān)性及其各影響因素之間的相關(guān)性的分析,最終選取坡度、高程、斷層距離、工程地質(zhì)巖組、水系距離以及道路距離6個影響因素作為滑坡、崩塌的評價因子;溝坡坡度、相對高差、溝壑密度、構(gòu)造密度、工程地質(zhì)巖組、道路密度及滑坡崩塌災(zāi)點(diǎn)密度7個影響因素作為泥石流的評價因子。

② 次級因子劃分

在分析各類型地質(zhì)災(zāi)害的發(fā)育與其各因子區(qū)間分布特征的基礎(chǔ)上,將各評價因子進(jìn)行次級評價因子劃分,滑坡崩塌共劃分出30個次級因子,泥石流劃分出38個次級因子。

(3) 評價因子的信息量計算

評價因子信息量的計算,是在統(tǒng)計各類地質(zhì)災(zāi)害在各因子區(qū)間內(nèi)分布的基礎(chǔ)上進(jìn)行的,基于各類地質(zhì)災(zāi)害的性質(zhì),滑坡崩塌進(jìn)行各區(qū)間內(nèi)災(zāi)點(diǎn)分布比率統(tǒng)計,而泥石流進(jìn)行各區(qū)間內(nèi)泥石流流域面積比率的統(tǒng)計。

① 滑坡

各評價因子滑坡總信息量值的大小順序為:工程地質(zhì)巖組>斷層距離>水系距離>高程>坡度>道路距離。根據(jù)各次級評價因子信息量值的正負(fù)判斷,滑坡在10°～45°坡度區(qū)間,<1 500 m高程區(qū)間,<500 m、1 500～2 500 m斷層距離區(qū)間,層狀較堅硬工程巖組、層狀軟硬相間工程巖組、層狀軟弱工程巖組、松散結(jié)構(gòu)軟弱土層4個工程地質(zhì)巖組,<500 m水系距離區(qū)間及<200 m道路距離區(qū)間內(nèi)易發(fā)生滑坡。

② 崩塌

各評價因子崩塌總信息量值的大小順序為:坡度>斷層距離>道路距離>水系距離>工程地質(zhì)巖組>高程。根據(jù)各次級評價因子信息量值的正負(fù)判斷,崩塌易發(fā)區(qū)間分別為:>25°坡度區(qū)間,<1 500 m高程區(qū)間,<1 500 m斷層距離區(qū)間,層狀較堅硬工程巖組、層狀軟硬相間工程巖組、層狀軟弱工程巖組3個工程地質(zhì)巖組,<500 m水系距離區(qū)間及<400 m道路距離區(qū)間。

③ 泥石流

各評價因子泥石流總信息量值的大小順序為:構(gòu)造密度>滑坡崩塌災(zāi)點(diǎn)密度>道路密度>坡度>工程地質(zhì)巖組>溝壑密度>相對高差。根據(jù)各次級評價因子信息量值的正負(fù)判斷,泥石流易發(fā)區(qū)間分別為:>25°坡度區(qū)間,100～1 600 m相對高差區(qū)間,2～3 km/km2溝壑密度區(qū)間,>0.2 km/km2構(gòu)造密度區(qū)間,層狀較堅硬工程巖組、層狀軟硬相間工程巖組兩個工程地質(zhì)巖組,1～5 km/km2道路密度區(qū)間及<0.5、2～3.5個/km2滑坡崩塌災(zāi)點(diǎn)密度區(qū)間。

表1 泥石流各次級評價因子易發(fā)性指數(shù)計算結(jié)果統(tǒng)計表

(4)評價因子權(quán)重計算

①構(gòu)建層次模型

由于本次研究不考慮次級評價因子的權(quán)重,只計算各評價因子的權(quán)重,因此只構(gòu)建由目標(biāo)圖層和指標(biāo)圖層組成的層次模型,基于滑坡崩塌和泥石流所選取的評價因子不同,分別建立滑坡崩塌及泥石流兩個層次模型。

②構(gòu)造判斷矩陣

在孕災(zāi)地質(zhì)條件分析的基礎(chǔ)上,結(jié)合各評價因子總信息量值的大小及專家打分,用Santy1-9標(biāo)度方法進(jìn)行調(diào)查區(qū)判斷矩陣的構(gòu)造,由于各評價因子分別與滑坡、崩塌、泥石流的相關(guān)性有所不同,分別進(jìn)行判斷矩陣的構(gòu)造。

③權(quán)重值計算

采用和積法對以上各判斷矩陣根據(jù)公式進(jìn)行特征向量的計算,也就是各評價因子的權(quán)重值。

④一致性檢驗

根據(jù)公式計算各判斷矩陣的λmax、CI及CR構(gòu)建的各判斷矩陣CR均<0.1,因此構(gòu)建的各判斷矩陣均具有滿意的一致性。

⑤易發(fā)性指數(shù)計算

表2 滑坡崩塌各次級評價因子易發(fā)性指數(shù)計算結(jié)果統(tǒng)計表

利用公式進(jìn)行各災(zāi)害類型各次級評價因子易發(fā)性指數(shù)計算,計算結(jié)果見表1、表2。

2.1.3 易發(fā)性評價結(jié)果

將滑坡、崩塌及泥石流易發(fā)性指數(shù)按等級重分類后從低至高依次賦值為1、2、3、4,再進(jìn)行疊加,形成縣域地質(zhì)災(zāi)害易發(fā)性指數(shù)分布圖,再將調(diào)查區(qū)劃分為極高、高、中、低4個易發(fā)性等級(圖1)。區(qū)內(nèi)地質(zhì)災(zāi)害易發(fā)性等級劃分及其地質(zhì)災(zāi)害分布情況見表3。

表3 地質(zhì)災(zāi)害易發(fā)性分級統(tǒng)計表

2.2 大關(guān)縣地質(zhì)災(zāi)害危險性評價及分區(qū)

地質(zhì)災(zāi)害危險性是指在某種誘發(fā)條件下,一定范圍內(nèi)在某一時間段發(fā)生地質(zhì)災(zāi)害可能性的大小。地質(zhì)災(zāi)害危險性評價除了要解決地質(zhì)災(zāi)害會在什么地方發(fā)生、以何種方式發(fā)生、發(fā)生的類型及其破壞方式等,還要了解發(fā)生的條件以及概率大小等問題。

地質(zhì)災(zāi)害危險性的評價,就是利用地質(zhì)學(xué)原理、數(shù)學(xué)方法等相結(jié)合,充分考慮地質(zhì)災(zāi)害系統(tǒng)的各種特征,對系統(tǒng)內(nèi)部各影響因素及其相互關(guān)系進(jìn)行綜合的、動態(tài)地分析,建立區(qū)域地質(zhì)災(zāi)害發(fā)展模式,進(jìn)而分析確定區(qū)域內(nèi)現(xiàn)今或是將來一段時間內(nèi)、特定條件下,地質(zhì)災(zāi)害發(fā)生的可能性。

圖1 地質(zhì)災(zāi)害易發(fā)性指數(shù)分布圖

2.2.1 評價指標(biāo)體系

本次研究在地質(zhì)災(zāi)害易發(fā)性評價基礎(chǔ)上,選擇10 a平均降雨量為誘發(fā)因素進(jìn)行地質(zhì)災(zāi)害危險性分析評價。以國家地球系統(tǒng)科學(xué)數(shù)據(jù)中心官方網(wǎng)站下載的2011～2020年逐月降雨數(shù)據(jù)為基礎(chǔ),通過提取與疊加分析,得出調(diào)查區(qū)內(nèi)10 a平均降雨柵格數(shù)據(jù)。結(jié)合滑坡崩塌點(diǎn)在降雨區(qū)間的統(tǒng)計分析,計算出各降雨區(qū)間的信息量值及降雨強(qiáng)度指數(shù),再疊加地質(zhì)災(zāi)害易發(fā)性指數(shù)得到地質(zhì)災(zāi)害危險性指數(shù)。

2.2.2 危險性評價

將地質(zhì)災(zāi)害易發(fā)性指數(shù)按極高、高、中、低4個等級依次賦值為4、3、2、1,再對應(yīng)乘以易發(fā)性指數(shù)權(quán)重值0.7對易發(fā)性指數(shù)柵格進(jìn)行重新賦值,最終與降雨量強(qiáng)度指數(shù)進(jìn)行疊加得到調(diào)查區(qū)的地質(zhì)災(zāi)害危險性指數(shù)圖。利用ARCGIS自然間斷點(diǎn)分級法將危險性指數(shù)劃分為極高、高、中、低4個危險等級區(qū)間,如圖2所示。區(qū)內(nèi)地質(zhì)災(zāi)害危險性等級劃分及其地質(zhì)災(zāi)害分布情況見表4。

圖2 地質(zhì)災(zāi)害危險性指數(shù)分布圖

表4 地質(zhì)災(zāi)害危險性分級統(tǒng)計表

2.3 大關(guān)縣地質(zhì)災(zāi)害易損性評價及分區(qū)

地質(zhì)災(zāi)害易損性的評價主體是承災(zāi)體,而承災(zāi)體易損性是一個包括物質(zhì)易損性、環(huán)境易損性、經(jīng)濟(jì)易損性及社會易損性等多重屬性的集合體。要客觀、合理地進(jìn)行地質(zhì)災(zāi)害易損性評價,就要選擇合理的指標(biāo),真實客觀地反映調(diào)查區(qū)承災(zāi)體的基本特征,而且要考慮信息的可提取性。

以收集到的最新大關(guān)縣年鑒及第三次全國土地普查數(shù)據(jù)為基礎(chǔ),基于承災(zāi)體易損性各屬性特點(diǎn),選擇人口密度、道路密度、耕地密度及GDP密度作為本次地質(zhì)災(zāi)害易損性評價指標(biāo)。根據(jù)各易損性評價指標(biāo)的易損性指數(shù)進(jìn)行賦值,再進(jìn)行疊加得到最終的地質(zhì)災(zāi)害易損性指數(shù)分布圖。將易損性指數(shù)分為極高、高、中、低4個等級區(qū)間,如圖3所示。區(qū)內(nèi)地質(zhì)災(zāi)害易損性等級劃分及其分布情況見表5。

表5 地質(zhì)災(zāi)害易損性分級統(tǒng)計表

2.4 大關(guān)縣地質(zhì)災(zāi)害風(fēng)險性評價及分區(qū)

地質(zhì)災(zāi)害風(fēng)險性評價是對地質(zhì)災(zāi)害的發(fā)生及其造成人員、經(jīng)濟(jì)等損失可能性的預(yù)評價,是在充分調(diào)查分析地質(zhì)災(zāi)害各孕災(zāi)地質(zhì)條件及其各誘發(fā)因素的基礎(chǔ)上,結(jié)合范圍內(nèi)人口、經(jīng)濟(jì)等可能造成損失的大小、程度及其抗災(zāi)強(qiáng)度等的分析評價,最終進(jìn)行的地質(zhì)災(zāi)害綜合性評價。

通過地質(zhì)災(zāi)害風(fēng)險性指數(shù)=地質(zhì)災(zāi)害危險性指數(shù)×易損性指數(shù)的計算,得到縣域地質(zhì)災(zāi)害風(fēng)險性指數(shù)分布圖(圖4)。根據(jù)風(fēng)險性指數(shù)大小,利用自然間斷點(diǎn)法將調(diào)查區(qū)劃分為極高、高、中、低4個風(fēng)險等級(圖5)。

圖4 地質(zhì)災(zāi)害風(fēng)險性指數(shù)分布圖

圖5 地質(zhì)災(zāi)害風(fēng)險性分區(qū)圖

表6 地質(zhì)災(zāi)害風(fēng)險性分區(qū)統(tǒng)計表

根據(jù)大關(guān)縣地質(zhì)災(zāi)害風(fēng)險性分區(qū)圖(參見圖5)綜合分析野外實際調(diào)查情況包括災(zāi)點(diǎn)分布及其危險性、風(fēng)險性等進(jìn)行判斷,結(jié)合孕災(zāi)地質(zhì)條件的分析,對調(diào)查區(qū)進(jìn)行地質(zhì)災(zāi)害風(fēng)險區(qū)的圈定,最終劃分出3個極高風(fēng)險區(qū)、5個高風(fēng)險區(qū)、2個中風(fēng)險區(qū)及4個低風(fēng)險區(qū),見表6?？芍獦O高風(fēng)險區(qū)占縣域面積的1.33%,高風(fēng)險區(qū)占縣域面積的10.02%,中風(fēng)險區(qū)占縣域面積的19.5%,低風(fēng)險區(qū)占縣域面積的69.15%。

根據(jù)風(fēng)險性分區(qū)圖以及野外實地調(diào)查情況,將這4個風(fēng)險性等級區(qū)域劃分為3級11個亞區(qū),包括5個重點(diǎn)防治區(qū)、2個次重點(diǎn)防治區(qū)和4個一般防治區(qū)。其中風(fēng)險性極高或高定為重點(diǎn)防治區(qū),風(fēng)險性中定為次重點(diǎn)防治區(qū),風(fēng)險性低定位一般防治區(qū)。重點(diǎn)防治區(qū)包括①木桿重點(diǎn)防治區(qū),分布于木桿向斜山的地貌區(qū)木桿河右岸斜坡下部,呈不規(guī)則片狀分布。②高橋重點(diǎn)防治區(qū),分布于高橋河背斜山地地貌區(qū)新開河匯入高橋河地帶,呈片狀分布。③吉利重點(diǎn)防治區(qū),分布于大關(guān)河北東向斷褶河谷地貌區(qū)中部,呈不規(guī)則環(huán)帶狀分布。④大關(guān)縣中部重點(diǎn)防治區(qū),分布于大關(guān)河、灑漁河、洛澤河及其主要支流谷坡區(qū)中下部,呈不規(guī)則帶狀、環(huán)帶狀分布。⑤天星重點(diǎn)防治區(qū),分布于洛澤河斷褶山地地貌區(qū)洛澤河支流陳家溝中下游流域,呈不規(guī)則片狀分布。次重點(diǎn)防治區(qū)包括①木桿次重點(diǎn)防治區(qū),分布于木桿向斜山地地貌區(qū)內(nèi)木桿河近東西段及主要支流下游谷坡區(qū),呈不規(guī)則片狀分布。②大關(guān)縣中部次重點(diǎn)防治區(qū),分布于大關(guān)河、灑漁河、洛澤河、高橋河下游及其主要支流谷坡區(qū)中上部,呈不規(guī)則帶狀、環(huán)帶狀分布。一般防治區(qū)包括①大關(guān)縣北西部一般防治區(qū),分布于調(diào)查區(qū)北西部,主要分布于木桿、高橋、壽山及吉利西部、悅樂北部,呈片狀分布。②大關(guān)縣中部一般防治區(qū),主要分布于吉利、壽山、天星3個鄉(xiāng)鎮(zhèn)交界地帶,呈不規(guī)則片狀分布。③大關(guān)縣東部一般防治區(qū),主要分布于天星鎮(zhèn)羅漢壩景區(qū)范圍內(nèi),呈不規(guī)則片狀分布。④大關(guān)縣南部一般防治區(qū),主要分布于上高橋、悅樂中南部、翠華西南部、玉碗南部及天星和翠華交界地帶,呈不規(guī)則片狀分布。

3 結(jié)論

根據(jù)地質(zhì)災(zāi)害風(fēng)險性評價體系,分別開展了研究區(qū)地質(zhì)災(zāi)害易發(fā)性評價、危險性評價和易損性評價。根據(jù)研究區(qū)特征,因地制宜,對于滑坡和崩塌選取了6個評價因子,泥石流選取了7個評價因子作為地質(zhì)災(zāi)害易發(fā)性評價指標(biāo),再選取10 a平均降雨量作為地質(zhì)災(zāi)害危險性評價指標(biāo),最后選擇了人口密度、道路密度、耕地密度及GDP密度作為易損性評價指標(biāo)。

采用信息量法模型和層次分析法進(jìn)行地質(zhì)災(zāi)害危險性評價其考慮了各因子的權(quán)重也考慮了各評級因子內(nèi)部各等級狀態(tài)的信息量值,使得對地質(zhì)災(zāi)害風(fēng)險性評價更加準(zhǔn)確、科學(xué)。

使用ArcGIS強(qiáng)大的數(shù)據(jù)處理、空間分析、統(tǒng)計功能完成了研究區(qū)地質(zhì)災(zāi)害易發(fā)性、危險性、易損性、風(fēng)險性評價和區(qū)劃,評價結(jié)果符合大關(guān)縣地質(zhì)災(zāi)害現(xiàn)狀,具有較高的準(zhǔn)確性,為大關(guān)縣地質(zhì)災(zāi)害防治以及地質(zhì)災(zāi)害風(fēng)險管控提供依據(jù)。

根據(jù)大關(guān)縣地質(zhì)災(zāi)害風(fēng)險性區(qū)劃,重點(diǎn)防治區(qū)建議防治措施以工程措施、群測群防及普適型監(jiān)測為主,輔以少量搬遷避讓和生物工程,在未防治前加強(qiáng)監(jiān)測預(yù)警。次重點(diǎn)防治區(qū)建議防治措施以群測群防及普適型監(jiān)測為主,工程措施、搬遷避讓為輔,在工程治理和搬遷避讓完工前應(yīng)加強(qiáng)監(jiān)測。一般防治區(qū)建議防治措施以加強(qiáng)地質(zhì)災(zāi)害防治培訓(xùn)宣傳,對管理人員提供防治專業(yè)知識培訓(xùn),對群眾提供不定期的基本常識的普及宣傳教育以及定期組織相關(guān)人員對地質(zhì)災(zāi)害隱患點(diǎn)進(jìn)行動態(tài)巡查、排查、核查。