2021年地質災害區(qū)劃與分區(qū)評價

**歐陽光明*創(chuàng)編2021.03.07*歐陽光明*創(chuàng)編2021.03.07第五章地質災害區(qū)劃與分區(qū)評價歐陽光明(2021.03.07)

第一節(jié)地質災害易發(fā)區(qū)劃分及分區(qū)評價一、評價思路與方法地質災害發(fā)育現(xiàn)狀是其易發(fā)性的客觀反映，要想準確的進行地質災害易發(fā)性分區(qū)，必須依賴遙感解譯和野外實際調查工作，本著這一思路，此次調查十分重視對基礎地質元素的搜集與分析，野外工作結束時已形成本區(qū)地質災害易發(fā)性基本輪廓，即初步的定性分區(qū)結果；同時考慮到地質環(huán)境條件的復雜性，通過對影響地質災害發(fā)育的諸多因素分析，采用半定量方法進行分區(qū)計算，作為對定性評價的補充，最后綜合兩種結果，形成本區(qū)地質災害的易發(fā)性分區(qū)圖。由于地質災害易發(fā)性的評價結果受到多種因素的影響，而這些因素本身存在著不確定性、模糊性以及各因素之間相互作用的復雜性；如何將復雜的地質因素盡可能的定量化，使分析和評價結果最大限度的符合客觀實際情況，是地質工作者廣為探討的問題。層次分析法(簡稱AHP)是一種定性和定量相結合的系統(tǒng)化、層次化的分析方法，它不僅適用于存在不確定性和主觀信息的情況，還允許以合乎邏輯的方式運用經驗、洞察力和直覺。由于它在處理復雜的決策問題上的實用性和有效性，因此在各學科模糊綜合評判中得以廣泛應用。**歐陽光明*創(chuàng)編2021.03.07**歐陽光明*創(chuàng)編2021.03.07地理信息系統(tǒng)在最近的30多年內得到了驚人的發(fā)展，廣泛應用于資源調查、環(huán)境評估、災害預測等眾多領域，借助GIS系統(tǒng)可以完成制圖、數(shù)字地形分析、空間決策支持、空間分析統(tǒng)計等任務，在GIS平臺上進行易發(fā)性區(qū)劃可以在一定程度上避免傳統(tǒng)區(qū)劃工作量大、工作強度大、工作精度不高以及主觀影響大的不足。本次工作擬采用基于層次分析法和GIS空間分析統(tǒng)計方法相結合的工作方法對區(qū)內地質災害易發(fā)性進行評價和區(qū)劃，主要技術路線和方法如下：1、確定評價單元和評價因子，利用層次分析法確定各因子和各要素的權值。2、對各評價因子指標進行量化，并采用歸一化數(shù)值變換方法統(tǒng)一量綱。3、在評價指標權值確定和數(shù)據歸一化的基礎上，利用GIS系統(tǒng)的空間分析功能進行數(shù)據的空間疊加與統(tǒng)計。4、經統(tǒng)計分析確定易發(fā)性區(qū)劃的分界點，將評價結果分成不同等級。5、在GIS分析成圖的基礎上綜合考慮各種因素，進行修改完善，最終編制工作區(qū)地質災害易發(fā)性區(qū)劃圖。二、理論與方法層次分析法是美國運籌學家T.L.Saaty于20世紀70年代中期提出的一種定性與定量相結合的多準則決策的系統(tǒng)分析方法，其基本原理是把復雜系統(tǒng)分解成目標、準則、方案等層次，在此基礎上進行定性和定量分析的決策。它把人的決策思維過程層次化、數(shù)量化、模型化，并用數(shù)學手段為分析、決策提供定量的依據，是一種對非定量事件進行定量分析的有效方法，特別是在目標因素結構復雜且缺少必要的數(shù)據情況下，需要將決策者的經驗判斷定量化時該法非常實用。該方法適用于多準則、多目標或無結構特征的復雜問題的決策分析，它按照各因子相互之間的內在支配關系，建立層次結構模型，通過因子的兩兩比較，建立判斷矩陣，進行層析排序，確定各因子的相對重要性。（一）層次分析法的主要步驟運用層次分析法建模，大體上可按下面4個步驟進行：1、建立層次結構模型：綜合分統(tǒng)所析系涉及的目標、范圍、準則，約束條件，確定綜合評價體系中各因索之間的關系，并根據指標的隸屬關系進行上下分層排列，形成危險性綜合評價體系層次結構。2、將問題中的各個要素劃分為不同的層次結構,以框架結構說明各層次之間的從屬關系。3、構造判斷矩陣：層次分析法的基礎主要靠對各層次兩兩元素相互重要程度差異給出判斷，并將這些判斷用數(shù)值表示出來，形成判斷矩陣，即兩兩比較矩陣。該矩陣是層次分析法的出發(fā)點，也是整個權值確定過程中的關鍵。在構造判斷矩陣的過程中，可以通過向專家發(fā)放問卷進行因子間的相對重要程度打分,并用矩陣表示打分4、檢驗層次分析結果如有誤差,需對判斷矩陣的元素取值通過專2021.03.072021.03.07**歐陽光明*創(chuàng)編2021.03.07*歐陽光明*創(chuàng)編家意見進行調整,從新運算。（二）評價指標判斷矩陣的構建信息是系統(tǒng)分析最基礎的數(shù)據。任何系統(tǒng)分析都要掌握一定的信息才能進行。層次分析法也需要有相應的信息作為分析的基礎，其信息主要來源于人們對不同層次中各個因素之間的相對重要性所做出的判斷。通過引入適當?shù)呐袛鄻硕葘⑦@些判斷用數(shù)值的形式表示出來構成判斷知陣。以便比較本層次各因素與某一因素之間的相對重要性。設B層次中的元素Bi,Bz,Bs,……從與上一層次A中的元素A有關系，則可以通過判斷知陣表示出來，如表511所示。表511判斷矩陣示意表AB1B2B3BnB1b11b12b13b1nB2b21b22b23b2nBnbn1Bn2Bn3bnn表中bij表示對危險性綜合評價而言Bi和Bj相對重要程度的數(shù)值表示。bij的取值是根據表512的T.Satty19標度含義來確定。表512判斷矩陣元需Aij的19度標度法標度含義1兩個因子相比較，兩者冋樣重要3兩個因子相比較，其中一個比另一個稍微重要5兩個因子相比較，其中一個相對另一個來說比較重要7兩個因子相比較，其中一個相對另一個來說非常重要9兩個因襲相比較，其中一個相對另一個來說極其重要2,4,6,8介于上面兩個相鄰判斷值的中間倒數(shù)若i與j相比較的判斷值為bij,則j與i比較的判斷值為1/bij三）特征值與特征向量計算根據判斷矩陣，利用線性代數(shù)知識，精確的求出T的最大特征根所對應的特征向量。所求特征向量即為各評價因素的重要性排序，經歸一化后即為同一層次相應因素對于上一層次某因素相對重**歐陽光明*創(chuàng)編2021.03.07**歐陽光明*創(chuàng)編2021.03.07**歐陽光明*創(chuàng)編2021.03.07**歐陽光明*創(chuàng)編2O21.O3.O7要性的排序權值。本次評價采用和積法進行求解，具體步驟如下：1、將判斷矩陣每一列歸一化：2、每一列經正規(guī)化的判斷矩陣按行相加：3、對向量，■A.■..：■■'做正規(guī)化處理，依次所得到的I「即為所求特征向量。4、計算判斷矩陣的最大特征根入max,式中：（TA）i表示向量TA的第i個元素。（四）一致性檢驗為避免其他因素對判斷矩陣的干擾,在實際應用中要求判斷矩陣滿足大體上的一致性,需進行一致性檢驗。只有通過檢驗,才能說明判斷矩陣在邏輯上是合理的,才能繼續(xù)對結果進行分析。對判斷矩陣進行一致性檢驗,計算公式：CR=CI/RI（1）式中，CR（consistencyratio）為一致性比例。當CRvO.lO時，認為判斷矩陣的一致性是可以接受的,否則應對判斷矩陣作適當修正。CI（consistencyindex）為一致性指標，按下式計算：CI=（入maxn）/（n1）（2）式中：max-判斷矩陣的最大特征根；n-成對比較因子的個數(shù)；RI（randomindex）-隨機一致性指標，可查表確定，如表513所示。表513平均隨機一致性指數(shù)RI階數(shù)n123456789RI000.580.91.121.241.321.411.45當CRvO.l時，就認為判斷矩陣具有滿意的一致性，否則就需要重新調整，直到具有滿意的一致性為止。三、評價因子選取及敏感性分析影響地質災害形成的自然因素眾多,歷史地震地質災害發(fā)生的數(shù)量、分布范圍、活動規(guī)模都直接反映了地層巖性、地形地貌、現(xiàn)存新老滑坡以及有關地震動力環(huán)境對地震誘發(fā)地質災害的控制作用;此外,土地利用、地下水、植物條件等因素也對震后地質災害形成起到一定程度的影響。本文采用統(tǒng)計學方法，對研究區(qū)山地災害點與各因子的每個屬性進行相對頻率組合的定量計算方法，綜合天水市秦州區(qū)震后地質災害發(fā)育情況，本次危險性區(qū)劃分析中選用了10個影響因子，主要包括：地質構造、地形坡度、海拔高程及水系發(fā)育情況等。（一）地質構造地質構造因素對地質災害點的發(fā)育控制作用十分明顯，在區(qū)域地質構造比較復雜，褶皺比較強烈，新構造運動比較活動的地區(qū)，地質災害比較發(fā)育。其影響主要表現(xiàn)在：①地質構造決定了地貌形態(tài)的分布，對地質災害發(fā)育的臨空條件起到間接的控制作用；②地質構造帶巖石破碎、風化嚴重，使得邊坡的連續(xù)性和完整性受到破壞，是地下水最豐富和活動的地區(qū)，降低了巖體的抗剪強度；③在構造應力作用下，巖體內節(jié)理、裂隙發(fā)育，為崩塌發(fā)育提供了條件；④活動斷層造成地表破裂，巖層結構發(fā)生破壞，非活動斷層作為地震波的反射界面，可能導致巖體的拉力破壞；⑥斷裂構造控制著水系的發(fā)育和人類工程活動的分布，對地質災害的威脅對象起到間接的控制作用研究中，通過GIS軟件緩沖區(qū)分析和數(shù)據統(tǒng)計功能，對研究區(qū)內災害點與斷裂距離分布關系做了統(tǒng)計：首先，對研究區(qū)內的斷裂做距離緩沖處理，分別得到0500，5001000，10002000及大于2000米四個緩沖區(qū)；然后利用GIS統(tǒng)計功能，對每個緩沖區(qū)內的災害數(shù)量、緩沖區(qū)面積進行統(tǒng)計，計算每個緩沖區(qū)內災害點密度。詳細數(shù)據如表514所示，災害點與斷裂的分布關系和敏感性關系，如圖511所示。表514研究區(qū)地質災害點與斷裂距離統(tǒng)計關系斷裂緩沖距離(m)<500500100010002000>2000緩沖區(qū)面積(km2)291.83248.35528.141285.43面積百分比(％)12.4010.5522.4454.61災害點數(shù)量342842100災害點密度(個/km2)0.120.110.080.08圖511斷裂構造與地質災害分布關系圖(二)地質災害頻率對天水市秦州區(qū)按1km*1km網格進行離散，形成2534個空間離散網格，依據災害點在網格上的分布，進行基于GIS的統(tǒng)計。該計算包括單元面積上災害發(fā)生的頻率及地質災害面積模數(shù)比。地質災害頻率比：設第(i,j)單元內災害頻率為f(i,j)，單元面積為S(ij)，單元內災害的頻率密度為pf(i,j),整個研究區(qū)面積為S,災害總數(shù)為f,總頻率密度為pf,則：第(i,j)單元格災害頻率比為：Rf(i,j)=pf(i,j)/pf其中,pf(i,j)=f(i,j)/S(i,j)；pf=f/S。地質災害面積模數(shù)比：設第(i,j)單元內災害體分布面積為Ss(i,j),單元面積為S(i,j)，單元內災害的面積模數(shù)為ps(i,j),整個研究區(qū)面積為S,災害點總面積為s,總面積模數(shù)為ps,則：第(i,j)單元格災害面積模數(shù)為：RS(i,j)=ps(i,j)/ps其中,ps(i,j)=s(i,j)/S(i,j)；ps=s/S。經計算可得,單個像元上最大出現(xiàn)地質災害的頻率為8。對空間災害點的頻率分布進行歸一,可形成圖512災害發(fā)生頻率歸一化分布圖及513地質災害面積模數(shù)比歸一化圖。圖512地質災害頻率比歸一化圖

圖513地質災害面積模數(shù)比歸一化圖坡度及坡度變率利用工作區(qū)1：5萬DEM數(shù)據提取坡度數(shù)據。根據前文中的分析，由于工作區(qū)內滑坡、崩塌災害主要分布于10?！?0。之間的斜坡,10°以下斜坡基本不發(fā)生滑坡、崩塌等災害,因此本次評價將60。以上斜坡的易發(fā)程度定義為1,10。以下易發(fā)程度定義為0,將坡度數(shù)據進行0—1之間的線性歸一化,得到坡度歸一化結果圖。坡度變化率是對地形基本因子——坡度變化情況進行量化的指標,由于斜坡拉張應力區(qū)的分布與斜坡坡度呈正相關聯(lián)系,因此隨著斜坡坡度變化率增大的斜坡坡腳地帶形成的最大剪應力也不斷增大,斜坡也就愈容易產生變形破壞。本次通過DEM對全區(qū)坡度變化率數(shù)據進行提取,然后進行0—1之間歸一化處理之后參與評價。圖514坡度歸一化圖圖515坡度變率歸一化（四）坡向及坡形變率坡形可以利用地表的曲率進行描述和量化，直線形和凸型斜坡在曲率上的體現(xiàn)是曲率＞0凹型坡和階梯型坡的曲率V0，因此，可利用ArcGIS平臺從DEM數(shù)據中提取調查區(qū)地表曲率信息，（平面曲率I、激活坡向數(shù)據。II、從【Surface】菜單中選擇【DeriveSlope】命令。III、生成平面曲率層面SlopeofAspect然后進行斜坡坡形的歸一化。由于滑坡和崩塌主要發(fā)育在直線型斜坡和凸型斜坡上，因此，當曲率V0時，坡面為凹型或階梯型，易發(fā)程度最低；當曲率〉0時，坡面為直線型和凸型，易發(fā)程度較高，按照曲率的大小進行0一1之間的線性歸一化，得到斜坡坡形指標歸一化結果（圖516，517）。圖516坡向歸一化圖圖517坡形曲率歸一化圖（五）海拔高程海拔高程對地質災害的控制作用主要表現(xiàn)在，一方面，海拔高程影響了地下水的分布，特別是潛水層的分布，松散巖土體構成的斜坡體內的地下水多為潛水，高程越高，潛水分布越少，對斜坡的影響越小。另一方面，海拔高程對人類活動范圍起控制作用，人類大多居住在海拔較低的河流沿岸，也多在海拔較低的地方進行生產活動，如開墾耕地。這些因素影響著地質災害的發(fā)育。因此，海拔高程也是地質災害危險性評價考慮的因素之一。利用DEM的高程信息進行求解，最后進行柵格化和歸一化處理（圖518）。圖518海拔高程與災害點關系（六）溝壑密度前已述及，溝壑密度是地形發(fā)育階段和地表抗蝕能力的重要特征值，對地質災害的發(fā)育有重要的影響作用。本次工作主要利用ARCGIS平臺中的Hydrology工具集，基于工作區(qū)1：5萬柵格DEM提取各流域單元的溝壑密度（圖519），主要步驟如下：1、對工作區(qū)dem數(shù)據進行洼地填平。2、利用GIS水文分析，得到提取區(qū)域的水流方向矩陣、水流累計矩陣。3、給定不同集水閥值，將水流方向累計矩陣中高于此閥值的格網連接起來得到矢量的溝壑網絡。4、對上一步提取的不同集水閥值下的溝谷網絡依據與實際形態(tài)的擬合程度進行對比分析，確定提取水文網和溝壑流域網絡最終的集水閥值。5、利用上一步確定的集水閥值分別提取水文網和流域溝壑網絡，并計算各流域的溝谷總長度和面積。6、依據得到的溝壑總長度和面積求得各流域的溝壑密度值。7、將各流域的溝壑密度進行歸一化處理并轉換為柵格數(shù)據參與評價。圖519溝壑密度圖（七）植被指數(shù)通過天水市秦州區(qū)1999年8月ETM+遙感數(shù)據，選擇近紅外波段4和可見光紅波段3,進行計算求取植被指數(shù)NDVI,之后將計算結果進行歸一化處理參與評價（圖5110）。圖5110植被指數(shù)圖（NDVI）（八）水系發(fā)育水系是誘發(fā)因子中對地質災害影響較大的一個因素,基于DEM提取區(qū)內水系發(fā)育情況,參與易發(fā)性評價,500m緩沖區(qū)內,近30%災害點落入區(qū)內，如緩沖區(qū)達到1000m,則災害點個數(shù)達到120余個,占60%以上。（圖5111）。圖5111水系發(fā)育與災害點關系圖（九）人類工程活動人類工程活動對地質環(huán)境的影響是極為復雜的,區(qū)內對地質環(huán)境改造較為強烈的人類活動即為公路、鐵路等線狀工程的修建,本次評價將工作區(qū)內的公路（包括國道、省道及縣主要干道）、鐵路做為基準線,間隔500m做緩沖區(qū)分析,分別向兩邊做三個緩沖區(qū),再經柵格化和歸一化處理后參與評價（圖5113）。圖5113人類活動歸一化圖四、因子權重確定（一）建立易發(fā)性評價層次結構地質災害易發(fā)區(qū)系指容易產生地質災害的區(qū)域,因此在選取評判因子時要依據工作區(qū)內地質災害發(fā)育的特點來選取,所選取的評判因子,應能全面反映區(qū)內地質災害的發(fā)育特點和孕災條件。本次**歐陽光明*創(chuàng)編2021.03.07**歐陽光明*創(chuàng)編202l.03.07評價以地質災害易發(fā)性做為目標層，選擇了發(fā)育因子、基礎因子和誘發(fā)因子構成準則層即二級評判因子，并選取了對地質災害易發(fā)性影響較為明顯的13個因子構成。對地質災害易發(fā)性評價指標體系進行分析，確定評價體系的層次結構，如圖5114所示：圖5114評價體系結構圖(二)構造判斷矩陣1、在層次結構中，對于從屬于(或影響)上一層的每個因素的同一層諸因素進行兩兩比較，比較其對于準則的重要程度，并按事前規(guī)定的標度定量化，構成矩陣形式，即判斷矩陣。判斷矩陣中各元素的數(shù)值由多名經驗豐富的專家集中群體智慧對各因素的相對重要性進行評估打分確定。2、根據層次分析法的基本原理，用10個指標建立起判斷矩陣，將各個指標間的相對重要程度表示出來。其中10個指標為地質構造(Bl);地質災害頻率(B2);坡度及坡度變率(B3);坡向及坡形變率(B4);海拔高程(B5);溝壑密度(B6);植被指數(shù)(B7);降雨量(B8);地震(B9);人類工程活動(B10)。通過向專家咨詢，進行指標因子間相互重要性進行打分，并結合研究區(qū)的敏感性評價分析結果，最終構建判斷矩陣如表所示：表5l5判斷矩陣標度層次單排序結果一覽表B722111212230.0554B83211/2231/21220.0611B91/231/21/2221/21/2130.0900B101/321/31/51/31/21/31/21/310.2190注：Bl—地質構造；B2—地質災害頻率；B3—坡度及坡度變率；B4—坡向及坡形變率；B5—海拔高程；B6—溝壑密度；B7—植被指數(shù)；B8—降雨量；B9—地震；B10—人類工程活動。3、求特征根通過計算可得特征向量，歸一化處理后即權重數(shù)W及最大特征根為：W=[0.ll67，0.l236，0.05ll，0.047l，0.0964，0.l396，0.0554，0.06ll，0.