GIS和層次分析法在地震地質(zhì)災害小區(qū)劃中的應用

我國地處環(huán)太平洋地震帶及喜馬拉雅地震帶的交匯處，斷裂活動頻繁，是一個多震的國家。同時地震誘發(fā)的地質(zhì)災害形式多樣，主要類型有滑坡、崩塌、泥石流、砂土液化、震陷、地裂縫。一旦發(fā)生地震，將有可能對人民的生命財產(chǎn)和社會經(jīng)濟的發(fā)展造成重大損失［1］。地震地質(zhì)災害的巨大危害性，使得對它的研究也越來越受到人們的重視。然而，地震地質(zhì)災害影響因素眾多，形成機制復雜的特點，給研究者的工作帶來很多難題。近年來隨著GIS技術的飛速發(fā)展，使這些問題得以解決。部分研究者將GIS、層次分析法單獨運用到地質(zhì)災害領域［2－4］，部分研究者將GIS與信息量法、層次分析法等數(shù)學模型結合運用到地質(zhì)災害方面［5－8］，取得良好的效果。GIS在地質(zhì)災害中得到了廣泛的應用，而在地震地質(zhì)災害小區(qū)劃中應用較少。根據(jù) 《工程場地地震安全性評價》中的規(guī)定，在區(qū)域地震小區(qū)劃工作中必須進行不同概率水平下的地震地質(zhì)災害小區(qū)劃。為此，筆者以GIS為平臺，結合層次分析法，以某地區(qū)為研究對象，在該地區(qū)地震地質(zhì)災害調(diào)查及地震動參數(shù)研究成果的基礎上，對其在50年超越概率10%下的地震地質(zhì)災害進行小區(qū)劃。

1 GIS與層次分析法的結合

層次分析法 （A HP）與GIS的結合主要從以下3個方面考慮：①利用Arcgis采集數(shù)據(jù)并對基礎數(shù)據(jù)進行處理，建立地震地質(zhì)災害基礎數(shù)據(jù)庫。②全面分析地震地質(zhì)災害單災種的影響因子，篩選出單災種的主要影響因子，采用層次分析法建立單災種層次模型確定各評價因子的權重，為單災種的空間分析提供定量指標；③在單災種評價結果的基礎上，再次利用層次分析法建立研究區(qū)地震地質(zhì)災害層次模型，計算單災種在區(qū)域地震地質(zhì)災害中的權重，在此基礎上，以Arcgis為平臺，利用其強大的空間分析功能進行研究區(qū)地震地質(zhì)災害小區(qū)劃。

2 地震地質(zhì)災害數(shù)據(jù)庫的建立

研究區(qū)地處渭河流域，區(qū)內(nèi)有多條河流穿過，為典型的黃土地貌，地層沉積連續(xù)、厚度大。場區(qū)及近場區(qū)有多條構造斷裂，受構造斷裂帶的控制，該區(qū)可劃分為黃土丘陵、黃土梁峁、黃土臺塬、洪積扇及河流階地5個地貌單元。其中河流階地從1級到5級，部分階地被剝蝕。經(jīng)過詳細的野外地質(zhì)調(diào)查分析，研究區(qū)未來受地震的影響性較大，存在的地震地質(zhì)災害主要有滑坡、黃土震陷及砂土液化。通過室內(nèi)相關資料收集、分析，野外地質(zhì)填圖 （1∶10000），對研究區(qū)內(nèi)所有的地質(zhì)災害案例及災害點經(jīng)行詳細的調(diào)查和記錄，包括災害類型、時間、地點、經(jīng)緯度坐標、規(guī)模及影響范圍。以1∶25000地形圖為底圖，對區(qū)內(nèi)所有的資料進行數(shù)字化和相關數(shù)據(jù)的錄入，利用Arccatalog建立該區(qū)的地震地質(zhì)災害數(shù)據(jù)庫。

3 地震地質(zhì)災害小區(qū)劃

1）地震滑坡評價 ①地震峰值加速度。地震峰值加速度是指地震時地面運動的最大加速度，它是反映地震動的重要指標，根據(jù)該區(qū)50年超越概率10%下地震動參數(shù)的分析結果，以Arcgis為平臺建立地震峰值加速度分區(qū)圖層。②滑坡災害面密度。災害面密度的定性分析是根據(jù)災害點的空間分布情況，統(tǒng)計每個區(qū)域的災害面密度，災害點面密度與滑坡危險性成正比?；旅婷芏鹊奶崛》椒ㄊ且訟rcgis為平臺，將滑坡分布圖層轉化為柵格圖層，以10 m×10 m劃分單元網(wǎng)格，然后以每個網(wǎng)格為圓心，半徑為500 m，逐個網(wǎng)格進行搜索，以圓內(nèi)部滑坡面積與圓面積之比作為該網(wǎng)格的滑坡災害的面密度。利用Arcgis的重分類功能將該指標分為4級。③坡度和坡高。坡高和坡度是引起滑坡的基本影響因子，隨著坡度和坡高的增大，引起斜坡剪應力增大，坡腳處出現(xiàn)應力集中，導致新滑坡的形成及老滑坡的復活。以Arcgis為平臺，利用1∶2.5萬DEM數(shù)據(jù)，提取研究區(qū)的坡度、坡高，按照調(diào)查統(tǒng)計結果進行分級。④地層巖性。地層巖性對滑坡活動的影響是一個比較復雜的因素，綜合考慮研究區(qū)地層巖性的各種因素，按照地層巖性對滑坡的影響程度，選取一種綜合參數(shù) ，在Arcgis中以地貌單元為基礎，對不同地層巖性進行賦值，取值范圍為1～10，將其劃分為4個等級，作為地層巖性的評價指標。

滑坡各影響因子的分級標準如表1所示。按層次分析法的步驟，首先以影響滑坡的5個因子為基礎建立層次分析法模型，采用Santy引入的標度方法結合專家意見建立判斷矩陣。判斷矩陣具有較好的一致性，對特征向量歸一化處理后，得到滑坡各因子的權重，見表2。

表1 滑坡影響因子分級標準

運用Arc GIS的柵格運算功能，將研究區(qū)5種評價因子圖層按照層次分析法所確定的權重，進行空間疊加計算，從而得到50年超越概率10%下該地震滑坡區(qū)劃柵格圖件。根據(jù)Arcgis的nat ural Break工具及實際調(diào)查結果，將研究區(qū)滑坡危險性劃分為4級：高危險區(qū)、中危險區(qū)、低危險區(qū)及無危險區(qū)，如圖1所示。

表2 滑坡危險性評價各因子判斷矩陣及權重

2）黃土震陷評價 研究區(qū)黃土震陷主要分布在2級以上階地及黃土臺塬上，該地區(qū)黃土地層沉積連續(xù)，厚度大，特別是馬蘭黃土的大孔隙、遇水易濕陷等特點，為黃土震陷提供了有利的條件。由于臺塬及階地馬蘭黃土坡度相對平緩，沉積在橫向上連續(xù)穩(wěn)定，因此，按水平層狀模型考慮黃土場地，采用一維計算模型，按分層總和思想計算黃土場地的震陷量。計算步驟如下：

①根據(jù)鉆探資料及試驗資料分層，分層總數(shù)為n。

②利用地震危險性分析得到的αmax計算各層可能遭受的動剪應力：

圖1 滑坡危險性評價圖

③根據(jù)計算所得動剪應力進行動三軸試驗，繪制每一層土的震陷曲線并確定震陷系數(shù)εpi。

④求黃土場地在垂向上的總震陷量S：

在典型的黃土塬及階地上，取原狀土樣數(shù)組，根據(jù)室內(nèi)動三軸試驗所得數(shù)據(jù)，采用上述方法，計算各個點在50年超越概率10%下黃土的累積震陷量。根據(jù)震陷量將其劃分為不震陷、輕微震陷、中等震陷及嚴重震陷4級。依據(jù)黃土震陷結果，利用Arcgis軟件繪制50年超越概率10%作用下研究區(qū)黃土震陷評價圖，如圖2所示。

圖2 黃土震陷及砂土液化評價圖

3）砂土液化評價 根據(jù)勘探孔地層、標貫資料及歷史資料判定，液化的區(qū)域主要分布在渭河河漫灘，呈零星分布，其余地貌單元均不具有液化的條件，綜合考慮各個鉆孔的資料，初步判定，可能液化的鉆孔共有7個。采用規(guī)范判別法和灰色理論模型法2種評判方法對研究區(qū)砂土液化進行了液化判別，在50年超越概率10%情況下，僅有3處可能發(fā)生砂土液化現(xiàn)象，以2種評判方法所得結果為基礎，基于Arcgis平臺，繪制50年超越概率10%作用下砂土液化評價圖，按液化指數(shù)將其劃分為不液化、輕微液化、中等液化、嚴重液化，如圖2所示。

4）地震地質(zhì)災害危險性指數(shù)模型的建立 把研究區(qū)看作是一個具有一定層次性和不確定性的系統(tǒng)，根據(jù)系統(tǒng)分析理論，采用層次分析法和GIS相結合的方法來計算地震地質(zhì)災害危險性指數(shù)，根據(jù)地震地質(zhì)災害的危險性指數(shù)進行小區(qū)劃。地震地質(zhì)災害危險指數(shù)定義為某一地區(qū)的某一柵格位置上的各種地質(zhì)災害產(chǎn)生的疊加影響的綜合，其表達式為：

式中，wj為第j個柵格的地震地質(zhì)災害危險性指數(shù)；Xj為單災種第j個柵格的計算值；δi為第i種地質(zhì)災害的權重；n為地震地質(zhì)災害的種類。

表3 地震地質(zhì)災害單災種判斷矩陣及權重

5）地震地質(zhì)災害小區(qū)劃結果 根據(jù)層次分析法計算各單災種的權重，如表3所示。利用Arcgis的空間分析功能，采用柵格加權疊加，對該區(qū)內(nèi)3種地震地質(zhì)災害柵格圖層進行空間分析，得到地震地質(zhì)災害危險性柵格圖，通過統(tǒng)計分析，利用頻率和頻數(shù)分布直方圖，確定地震地質(zhì)災害危險性指數(shù)的分區(qū)閾值，使用Arcgis的重分類功能對獲得的柵格圖進行重分類，同時對所得柵格圖進行矢量化，綜合考慮各種因素，修改完善后最終得到50年超越概率10%下研究區(qū)的地震地質(zhì)災害小區(qū)劃圖，如圖3所示。

根據(jù)以上方法，將研究區(qū)地震地質(zhì)災害劃分為4個等級，即無危險區(qū)、低危險區(qū)、中危險區(qū)、高危險區(qū)。在50年超越概率10%情況下，研究區(qū)共劃分為195個小區(qū)，其中無危險區(qū)43個、低危險區(qū)67個、中危險區(qū)61個、高危險區(qū)24個，各個等級區(qū)域面積分別為133.2、89.93、36.5、21.3k m3，面積百分比分別為47.5%、32.0%、13.0%、7.5%。

圖3 地震地質(zhì)災害小區(qū)劃圖

［1］劉鳳民，張立海，劉海青，等 .中國地震次生地質(zhì)災害危險性評價 ［J］.地質(zhì)力學學報，2006，12（2）：127－131.

［2］陳曉利，葉洪，程菊紅 .GIS技術在區(qū)域地震滑坡危險性預測中的應用——以龍陵地震滑坡為例 ［J］.工程地質(zhì)學報，2006，14 （3）：333－338.

［3］姚玉增，任群智 .層次分析法在山地地質(zhì)災害危險性評價中的應用——以遼寧凌源地區(qū)為例 ［J］.水文地質(zhì)工程地質(zhì)，2010，37 （2）：130－134.

［4］莊建琦，崔鵬，葛永剛，等 .“5·12”汶川地震崩塌滑坡危險性評價——以都汶路沿線為例 ［J］.巖石力學與工程學報，2010，29 （2）：3735－3742.

［5］光磊 .基于GIS的信息量法在滑坡危險性評價中的應用 ［J］.地質(zhì)與資源，2005，14（3）：231－233.

［6］阮沈勇，黃潤秋 .基于GIS的信息量法模型在地質(zhì)災害危險性區(qū)劃中的應用 ［J］.成都理工學院院報，2001，28（1）：89－92.

［7］杜軍，楊青華 .基于GIS與AHP耦合的汶川震后次生地質(zhì)災害風險評估 ［J］.中國水土保持，2009，11（3）：14－17.

［8］張永軍，侯云龍，劉武 .層次分析法和GIS空間分析法在蘭州市區(qū)地質(zhì)災害易發(fā)性評價中的應用 ［J］.甘肅地質(zhì)，2009，18（4）：84－88.