小流域地震地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)

李佩佩，沈軍輝，燕俊松，陳 亮

(地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國家重點(diǎn)實(shí)驗(yàn)室(成都理工大學(xué))，四川 成都 610059)

小流域地震地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)

李佩佩，沈軍輝，燕俊松，陳 亮

(地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國家重點(diǎn)實(shí)驗(yàn)室(成都理工大學(xué))，四川 成都 610059)

本文以北川縣楊家溝為小流域地震地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)的典型實(shí)例，根據(jù)地震地質(zhì)災(zāi)害特征，選取坡度、坡向、地層巖性、距斷裂帶距離、斜坡結(jié)構(gòu)、高程和坡形7個(gè)評(píng)價(jià)因子，在各評(píng)價(jià)因子對(duì)地震地質(zhì)災(zāi)害影響程度分析基礎(chǔ)上，采用層次分析法確定評(píng)價(jià)因子的權(quán)重，并在GIS平臺(tái)下對(duì)各評(píng)價(jià)因子進(jìn)行綜合分析處理，得到小流域地震地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)圖。研究表明，楊家溝地震地質(zhì)災(zāi)害最大影響因子為距斷裂帶距離，其次為斜坡結(jié)構(gòu)、坡度、地層巖性、高程、坡形及坡向；地震地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)預(yù)測(cè)表明，楊家溝流域內(nèi)高危險(xiǎn)區(qū)面積8.98 km2，占流域面積的40.43%。對(duì)比楊家溝實(shí)際地震地質(zhì)災(zāi)害分布情況，占面積70.74%或占總數(shù)66.38%的地震地質(zhì)災(zāi)害位于高危險(xiǎn)區(qū)內(nèi)，表明危險(xiǎn)性評(píng)價(jià)成果可信度較高。研究成果對(duì)小流域地震地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)方法研究具有一定的意義。

危險(xiǎn)性；地震地質(zhì)災(zāi)害；層次分析法；GIS

0 引言

近年來的地震活動(dòng)特征表明地殼活動(dòng)進(jìn)入相對(duì)活躍期，地震誘發(fā)的一系列地質(zhì)災(zāi)害給災(zāi)區(qū)人民帶來了沉重災(zāi)難?，F(xiàn)階段，相關(guān)部門組織實(shí)施縣(市)地質(zhì)災(zāi)害調(diào)查與區(qū)劃工作取得了一定成果，但總體上評(píng)價(jià)工作精度較低，評(píng)價(jià)成果還不能有效地應(yīng)用于小流域(流域面積在5～30 km2的自然閉合集水區(qū)[1])的防災(zāi)減災(zāi)之中。小流域流域面積小，人口分散，地質(zhì)災(zāi)害防災(zāi)減災(zāi)工作難度大。為了使評(píng)價(jià)成果更好地應(yīng)用于小流域防災(zāi)減災(zāi)的具體實(shí)施，為流域建設(shè)規(guī)劃及地質(zhì)災(zāi)害高危險(xiǎn)區(qū)的綜合治理，開展高精度小流域地震地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)是十分必要的。

本文借鑒現(xiàn)有地震地質(zhì)災(zāi)害研究成果[2-17]，結(jié)合野外調(diào)研，并考慮原始資料獲得的難易程度，選取影響地震地質(zhì)災(zāi)害評(píng)價(jià)因子，以北川縣楊家溝為典型實(shí)例，采用層次分析法，基于GIS 平臺(tái)，進(jìn)行小流域地震地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)。

1 地震地質(zhì)災(zāi)害評(píng)價(jià)因子選取及量化賦值

地震地質(zhì)災(zāi)害發(fā)生主要受坡度、坡向、地層巖性、斜坡結(jié)構(gòu)、高程、坡形、地面粗糙度等斜坡自身性質(zhì)的控制，也受到地震動(dòng)峰值加速度、地震烈度等觸發(fā)因素的影響。本文借鑒現(xiàn)有的地震地質(zhì)災(zāi)害形成條件、發(fā)育特征等研究成果，結(jié)合野外調(diào)研及小流域自身特征，選取了坡度、坡向、地層巖性、距斷裂帶距離、斜坡結(jié)構(gòu)、高程、坡形七個(gè)影響因素為地震地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)的評(píng)價(jià)因子。

1.1 坡度

坡度對(duì)地震地質(zhì)災(zāi)害發(fā)生的影響主要有兩方面。一方面反映在該斜坡地段是否具備有效臨空面，而有效臨空面又與地形坡度有很大關(guān)系[9]，一般坡度越大地震地質(zhì)災(zāi)害越發(fā)育。另一方面，隨著坡度的增加，坡面有效受雨面積不斷減小，導(dǎo)致坡面徑流量隨坡度增大而減少，沖刷能力相應(yīng)減弱，有效受雨面積大的斜坡在地震作用下更容易發(fā)生崩滑。受上述兩方面共同影響，斜坡的坡度越大，地震地質(zhì)災(zāi)害愈發(fā)育，但當(dāng)坡度大于某個(gè)值時(shí)，穩(wěn)定性又會(huì)增加[10]。結(jié)合汶川地震地質(zhì)災(zāi)害的相關(guān)研究成果，當(dāng)斜坡的坡度在20°～45°之間時(shí)，地質(zhì)災(zāi)害最發(fā)育，大于45°次之，小于20°災(zāi)害發(fā)育最少[9-11]。因此，將坡度按危險(xiǎn)性程度進(jìn)行賦值，危險(xiǎn)性越大，賦值越高，反之越低，最低賦值為1。

1.2 坡向

無論緯度、高程或坡度，其不同坡向上的日照、太陽輻射、溫度、濕度和降水等級(jí)都有很大的差異[13]，它間接地影響斜坡穩(wěn)定性。根據(jù)應(yīng)力波的理論[14]，巖土體普遍存在裂隙，當(dāng)壓縮的地震波傳播時(shí)遇到自由面，將反射成一個(gè)大小相等、方向相反的拉伸波，壓縮波和拉伸波的疊加使巖土體產(chǎn)生拉裂破壞。唐春安等[8]通過模擬動(dòng)態(tài)沖擊載荷模型試驗(yàn)，應(yīng)力波從左側(cè)向右側(cè)傳播，在右側(cè)首先產(chǎn)生拉裂破壞，隨著時(shí)間增加，最終斜坡破壞，這與實(shí)地調(diào)查中背靠震源方向滑坡明顯比面向震源方向滑坡發(fā)育相吻合。表明應(yīng)力波方向不同(與震源位置有關(guān))，邊坡的破壞位置不同。因此，按照地震波傳播方向的關(guān)系進(jìn)行分類賦值，背靠震源方向賦值最高，面向震源方向賦值最低為1，其他的方向介于1與最高值。

1.3 地層巖性

地震地質(zhì)災(zāi)害發(fā)育程度與巖性關(guān)系密切，不同巖石地層單元地質(zhì)災(zāi)害分布數(shù)量差別較大。一般，硬巖和較硬巖地層中常發(fā)生崩塌災(zāi)害，較堅(jiān)硬巖類中通?；伦罴邪l(fā)育，這與斜坡在非地震破壞的巖性特征有一定差別[16]。地震地質(zhì)災(zāi)害發(fā)育最多的巖性為砂巖、粉砂巖、灰?guī)r、閃長(zhǎng)巖、花崗巖等，其次為頁巖、砂礫巖、黏土巖為主、偶夾板巖、千枚巖等。根據(jù)上述不同巖性對(duì)地震地質(zhì)災(zāi)害的貢獻(xiàn)，依次賦值。

1.4 距斷裂帶距離

地震地質(zhì)災(zāi)害的空間分布主要受發(fā)震斷裂帶的控制。地震波對(duì)坡體的強(qiáng)烈沖擊是觸發(fā)地質(zhì)災(zāi)害的決定性因素[17]。距發(fā)震斷層越近，地震波對(duì)坡體的作用越強(qiáng)烈，就越容易破壞。根據(jù)相關(guān)研究成果統(tǒng)計(jì)[17-19]，80%以上的地震地質(zhì)災(zāi)害分布在斷層兩側(cè)10 km以內(nèi)，強(qiáng)發(fā)育區(qū)在斷層7 km范圍內(nèi)。按距斷裂帶距離由遠(yuǎn)到近依次由小到大賦值。

1.5 斜坡結(jié)構(gòu)

斜坡結(jié)構(gòu)綜合體現(xiàn)了地層產(chǎn)狀與斜坡坡度和坡向的空間狀況及組合形式[20]，在很大程度上決定了斜坡巖土體變形的方式和強(qiáng)度[21]，對(duì)地質(zhì)災(zāi)害的分布起著重要的作用。根據(jù)斜坡坡度、坡向和巖層傾向、傾角四者在空間上的相互組合[22]，將斜坡結(jié)構(gòu)類型劃分為：順向飄傾坡、順向?qū)用嫫?、順向伏傾坡、順斜坡、橫向坡、逆斜坡、和逆向坡7類。在地震作用下，順向、斜向—順向坡內(nèi)有利于層面與節(jié)理面的貫通，更容易發(fā)生破壞。將順向坡賦值最高，將逆向坡賦值為1。其中順向坡中的飄傾坡，傾角小于坡度，最危險(xiǎn)；其次為坡度大于10°，傾角與坡度相等的層面坡；最后為傾角大于坡度的伏傾坡。

1.6 高程

根據(jù)斜坡地震動(dòng)放大效應(yīng)[23]，山頂?shù)牡卣鸺铀俣缺壬侥_放大效應(yīng)明顯。主要由于強(qiáng)震作用下，地震波豐富的波長(zhǎng)成分與地形尺寸耦合作用，產(chǎn)生的動(dòng)力效應(yīng)超過或遠(yuǎn)遠(yuǎn)超過斜坡巖土體的強(qiáng)度進(jìn)而產(chǎn)生破壞。在汶川地震中74.8%的崩滑地質(zhì)災(zāi)害主要分布在800～2 000 m[24]，其中1 000～1 500 m密度最大，高于2 000 m地質(zhì)災(zāi)害明顯減少。故將高程800～2 000 m賦值最高，高于2 000 m賦值次之，低于800 m的賦值為1。

1.7 坡形

坡面形態(tài)主要有凹形坡、斜線坡和凸形坡[13]三類。坡形影響地震波的傳播方向和能量的釋放,從而影響崩塌、滑坡的發(fā)育。一般用地面曲率來衡量坡形，地面曲率是對(duì)地形表面扭曲程度的定量化度量因子[25]，地面曲率決定地表及土壤中物質(zhì)移動(dòng)的相對(duì)速度，以及地表徑流對(duì)坡面的侵蝕程度。地面曲率大于0為凸形坡，表明移動(dòng)加速，侵蝕強(qiáng)度大；小于0為凹形坡，表明移動(dòng)減速，侵蝕強(qiáng)度小。凸形坡，從基部到頂部，坡度隨坡形而急劇復(fù)雜變化，導(dǎo)致凸形坡上部坡緣的高應(yīng)力集中區(qū),在其它條件一致的情況下，斜線坡和凹形坡對(duì)應(yīng)力在剖面上的分布影響較小[22]。因此，凸形坡比其他坡形穩(wěn)定性要差。根據(jù)斜坡曲率值由正到負(fù)依次由大到小賦值，最小賦值為1。

2 層次分析法原理

將地震地質(zhì)災(zāi)害影響的因素視為系統(tǒng)，按照分解、比較、判斷、綜合的思維方法進(jìn)行決策，并按各因素間的相互關(guān)系分層次，形成一個(gè)多層次的結(jié)構(gòu)體系。建立判別因子比較矩陣，最后求得各因子的相對(duì)權(quán)重。

根據(jù)前人的研究數(shù)據(jù)資料、專家意見和本文作者對(duì)北川縣城野外考察的認(rèn)識(shí)，加以平衡后給出的[26]。一般引用Saaty等學(xué)者提出1～9標(biāo)度對(duì)重要性判斷結(jié)果進(jìn)行量化,依次對(duì)兩兩指標(biāo)進(jìn)行比對(duì)賦值，將賦值結(jié)果矩陣化，得n×n階矩陣。本文用方根法對(duì)判斷矩陣求解特征向量和最大特征值。

由于層次分析法是解決一些現(xiàn)象測(cè)度的一種方法(表1)，結(jié)果的不一致性是客觀存在的，但這種不一致性需要一個(gè)度，因此一致性檢驗(yàn)就非常有必要。一致性指標(biāo)為：

C.R.=C.I./R.I.

其中 C.I.=(λmax-n)/(n-1)，λmax為判斷矩陣A的最大特征值，R.I.為隨機(jī)一致性指標(biāo)[28]。對(duì)于一致性比例，一般在C.R.=0時(shí)，可以稱A是完全一致性矩陣，當(dāng)C.R.<0.1時(shí)，該判斷矩陣具有一致性，當(dāng)C.R.>0.1時(shí)，則該判斷矩陣不具有一致性。

表1 層次分析法評(píng)價(jià)尺度

3 實(shí)例應(yīng)用

3.1 楊家溝概況

楊家溝位于汶川地震重災(zāi)區(qū)北川縣陳家壩鄉(xiāng)，為侵蝕構(gòu)造中山地形，地形坡度一般為20°～40°，流域最高高程2 200 m，溝口處高程800 m，相對(duì)高差500～1 000 m。溝道全長(zhǎng)約9 km，流域面積約20.69 km2。楊家溝溝谷上游呈“V”型，下游呈“U”型，谷寬50～200 m。地層巖性主要為志留系韓家店組(Sh)頁巖、絹云母千枚巖；志留系茂縣組(Sm)絹云母千枚巖夾細(xì)粒粉砂巖；奧陶系寶塔組(Ob)塊狀粒晶灰?guī)r、生物碎屑灰?guī)r；寒武系油房組(∈y)變凝灰質(zhì)砂巖；寒武系邱家河組(∈q)深灰-黑灰色碳質(zhì)板巖，變碳質(zhì)粉砂巖,局部有四宗坪侵入體(Ts)為輝綠玢巖。NE向北川逆沖斷層從楊家溝溝口通過，斷層傾向北西，傾角60°～70°。受構(gòu)造走向控制，巖層走向以北東走向?yàn)橹鳌蟽?nèi)斜坡結(jié)構(gòu)主要以順斜坡、橫向坡、逆斜坡、和逆向坡為主。根據(jù)中國地震局對(duì)汶川地震圈定的地震烈度圖，楊家溝地震烈度達(dá)Ⅺ度。

受上述地形地貌、地層巖性及構(gòu)造等的控制，汶川地震時(shí)楊家溝共發(fā)育大小地質(zhì)災(zāi)害119處(圖1)，類型以中小型崩滑為主。距離斷裂帶4 km范圍內(nèi)發(fā)育最為密集，造成大量人員傷亡。而且耕地、林地及農(nóng)作物大面積受損，水、電、通訊基本被摧毀，鄉(xiāng)村道路60%被破壞。

圖1 楊家溝流域特征Fig.1 The characteristics of Yangjiagou watershed

3.2 判斷矩陣的構(gòu)造及一致性檢驗(yàn)

對(duì)影響地震地質(zhì)災(zāi)害發(fā)育的7個(gè)評(píng)價(jià)因子, 按照(表2) 規(guī)則進(jìn)行兩兩比較，矩陣結(jié)果見表3。由于小流域內(nèi)降雨差異基本可以忽略，但是有些因子對(duì)降雨敏感性強(qiáng)，故在因子對(duì)比時(shí)將降雨影響明顯的因子在對(duì)比時(shí)賦值有所調(diào)整。計(jì)算得到矩陣的最大特征值為7.428 8，各個(gè)評(píng)價(jià)因子的權(quán)重為{0.129 8，0.023 4，0.190 8，0.459 9，0.094 2，0.051 3，0.051 3}(表2)。一致性指標(biāo)為C.R.=0.054<0.1，表明各評(píng)價(jià)因子相互比較在合理范圍內(nèi)。評(píng)價(jià)因子權(quán)重計(jì)算結(jié)果表明，距斷裂帶距離是對(duì)地質(zhì)災(zāi)害影響最嚴(yán)重的因子，其他依次是坡度、斜坡結(jié)構(gòu)、地層巖性、高程、坡形、坡向。

表2 評(píng)價(jià)因子權(quán)重的成對(duì)比較矩陣

3.3 數(shù)據(jù)處理及危險(xiǎn)性評(píng)價(jià)

在ARCGIS 10.1平臺(tái)上，采用Transverse_Mercator投影系統(tǒng)。本文用地震前1∶50 000地形圖上的等高線、高程點(diǎn)得到數(shù)字高程模型(DEM)。坡度、坡向、坡形3個(gè)因子圖層均用ARCGIS10.1表面分析獲得。巖性與斷裂均來自于1∶50 000陳家壩幅地質(zhì)圖。斜坡結(jié)構(gòu)根據(jù)現(xiàn)場(chǎng)調(diào)查得到的產(chǎn)狀點(diǎn)數(shù)據(jù)，通過克里金插值法網(wǎng)格化處理，得到空間上連續(xù)的巖層傾向、傾角網(wǎng)格化數(shù)據(jù)，該數(shù)據(jù)柵格單元要與坡度、坡向保持一致。當(dāng)作圖區(qū)域地質(zhì)構(gòu)造比較復(fù)雜時(shí)，不能簡(jiǎn)單地進(jìn)行空間插值，應(yīng)當(dāng)進(jìn)行分區(qū)分片處理——如斷層、褶皺等進(jìn)行巖層產(chǎn)狀分塊處理[21]。然后利用ARCGIS10.1柵格計(jì)算器用巖層傾向減去坡向，巖層傾角減去坡度，將得到的兩個(gè)圖層進(jìn)行疊加并重分類賦值，得到最終的斜坡結(jié)構(gòu)圖層。由于斷裂帶在圖上為線性表示，首先利用ArcGIS 10.1鄰域分析生成半徑差為1 km的多環(huán)緩沖區(qū)，轉(zhuǎn)化為與其他因子?xùn)鸥駟卧嗤臇鸥駡D層。按照第2節(jié)中各評(píng)價(jià)因子中各個(gè)級(jí)別對(duì)地震地質(zhì)災(zāi)害的影響程度重分類標(biāo)準(zhǔn)，將楊家溝各因子進(jìn)行分類賦值(表3)，得到各因子分區(qū)圖。最后應(yīng)用ARCGIS10.1平臺(tái)空間分析模塊中的柵格計(jì)算工具對(duì)各因子加權(quán)疊加，得到最終的研究區(qū)危險(xiǎn)性評(píng)價(jià)分區(qū)圖。

表3 楊家溝地震地質(zhì)災(zāi)害評(píng)價(jià)因子及其分類賦值

最后參照地質(zhì)災(zāi)害風(fēng)險(xiǎn)度分級(jí)標(biāo)準(zhǔn)，結(jié)合研究區(qū)實(shí)際情況，使用Arcgis10.1重分類將研究區(qū)劃分為高度危險(xiǎn)區(qū)、中度危險(xiǎn)區(qū)、低度危險(xiǎn)區(qū)、極低度危險(xiǎn)區(qū)4類，得到研究區(qū)地震地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)結(jié)果(表4、圖2)。

表4 地質(zhì)災(zāi)害危險(xiǎn)性統(tǒng)計(jì)

圖2 楊家溝危險(xiǎn)性分區(qū)成果圖Fig.2 Risk zoning map of Yangjiagou

根據(jù)圖2對(duì)楊家溝地震地質(zhì)災(zāi)害進(jìn)行初步分析，地質(zhì)災(zāi)害高危險(xiǎn)區(qū)集中在距斷裂帶4 km的區(qū)域及其一些零散區(qū)域內(nèi)。高危險(xiǎn)區(qū)約占楊家溝總面積的43.4%，地質(zhì)災(zāi)害70%以上分布在該區(qū)，且面積大于2×104m2的地質(zhì)災(zāi)害幾乎全部在該區(qū)域內(nèi)，與實(shí)際比較一致。在高危險(xiǎn)區(qū)無論是地質(zhì)災(zāi)害數(shù)量還是面積遠(yuǎn)遠(yuǎn)大于其他區(qū)，其中汶川地震發(fā)生的沙壩里滑坡(圖1中的A、圖3)形成堰塞湖，對(duì)下游造成嚴(yán)重威脅。

流域內(nèi)一老滑坡體(圖1中的B、圖4)，在本次地震中除局部滑塌外，整體未有明顯變形復(fù)活，其原因是由于老滑坡體早期的高位勢(shì)能得到充分耗散，并且深厚覆蓋層堆積體吸收地震波的能量，使得老滑坡體在汶川強(qiáng)烈的地震荷載作用下只出現(xiàn)局部垮塌,沒有出現(xiàn)整體失穩(wěn)。在圖2成果圖中為低危險(xiǎn)區(qū)。

圖3 沙壩里滑坡Fig.3 Shabali landslide

圖4 老滑坡體Fig.4 Old landslide

4 結(jié)論

(1)本文以楊家溝為小流域地震地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)的典型實(shí)例，選取坡度、坡向、地層巖性、距斷裂帶距離、斜坡結(jié)構(gòu)、高程、坡形7個(gè)評(píng)價(jià)因子。用層次分析法構(gòu)建各評(píng)價(jià)因子比較矩陣，得到各評(píng)價(jià)因子的權(quán)重，得出對(duì)地震地質(zhì)災(zāi)害影響最大的因子為距發(fā)震斷裂帶距離，其次為斜坡結(jié)構(gòu)、坡度、地層巖性、高程、坡形及坡向。

(2)基于GIS 平臺(tái)分析處理得到的楊家溝危險(xiǎn)性評(píng)價(jià)成果圖，將楊家溝按危險(xiǎn)性程度分為高危險(xiǎn)區(qū)、中危險(xiǎn)區(qū)、低危險(xiǎn)區(qū)與極低危險(xiǎn)區(qū)4 類，其中高度危險(xiǎn)區(qū)面積為8.99 km2，占流域總面積的43.4%，主要集中于距NE向北川斷裂4 km范圍內(nèi)，其他區(qū)域僅零散分布。

(3)評(píng)價(jià)結(jié)果與汶川地震地質(zhì)災(zāi)害實(shí)際情況進(jìn)行對(duì)比，占面積70.74%或占總數(shù)66.38%的地震地質(zhì)災(zāi)害位于高危險(xiǎn)區(qū)內(nèi)，危險(xiǎn)性評(píng)價(jià)成果可信度較高。表明在地震地質(zhì)災(zāi)害影響因素分析及合理選用影響因子基礎(chǔ)上，采用層次分析法應(yīng)用于小流域地震地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)是可行的。

[1] 趙愛軍.小流域綜合治理模式研究[D].湖北，華中農(nóng)業(yè)大學(xué)，2005:19-23. ZHAO Aijun.The Research on comprehensive control model in small watershed[D].Huazhong Agricultural University，2005:19-23

[2] 唐川,朱靜,張翔瑞.GIS支持下的地震誘發(fā)滑坡危險(xiǎn)區(qū)預(yù)測(cè)研究[J].地震研究，2001, 24(1):73-81. TANG Chuan,ZHU Jing,ZHANG Qiangrui. GIS Based earthquake triggered landslide hazard prediction[J].Journal of Seismological Research,2001,24(1):73-81

[3] 黃潤秋,許強(qiáng).中國典型災(zāi)難性滑坡[M].北京，科學(xué)出版社,2008:20-40. HUANG Runqiu,XU Qiang. Typical catastrophic landslides in China[M].Beijing,Science Press,2008:20-40.

[4] 黃潤秋,李為樂.汶川地震觸發(fā)崩塌滑坡數(shù)量及其密度特征分析[J].地質(zhì)災(zāi)害與環(huán)境保護(hù), 2009,20(3):1-7. HUANG Runqiu, LI Weile.Analysis on the number and density of landslides triggered by the 2008 Wenchuan earthquanke, China[J].Journal of Geological Hazards and Environment Preservation,2009,20(3):1-7.

[5] 姚鑫,許沖,戴福初,等.四川汶川Ms8級(jí)地震引發(fā)的滑坡與地層巖性、坡度的相關(guān)性.地質(zhì)通報(bào), 2009,28(8):1156-1162. YAO Xin, XU Chong, DAI Fuchu,et al.Contribution of strata lithology and slope gradient to landslides triggered by Wenchuan Ms 8 earthquake, Sichuan, China[J].Geological Bulletin of China, 2009, 28(8):1156-1162

[6] 殷躍平,張永雙,馬寅生,等.青海玉樹MS 7. 1級(jí)地震地質(zhì)災(zāi)害主要特征[J].工程地質(zhì)學(xué)報(bào),2010,18(3):289-296. YIN Yueping,ZHANG Yongshuang,MA Yinsheng,et al.Research on major characteristics of geohazards in duced by the Yushu Ms7.1 earthquake[J].Journal of Engineering Geology,2010,18(3):289-296.

[7] 劉應(yīng)輝.汶川地震區(qū)都汶公路沿線崩塌滑坡災(zāi)害特征與評(píng)價(jià)[D].蘭州，蘭州大學(xué),2009:18-22. LIU Yinghui.The charaeteristie and evaluation of collapse and landslide disaster along Du-wen highway in wenchuan earthquakeregion[D].Lanzhou,Lanzhou University,2009:18-22.

[8] 楊德宏，范文.基于ArcGIS的地質(zhì)災(zāi)害易發(fā)性分區(qū)評(píng)價(jià)—以旬陽縣為例[J].中國地質(zhì)災(zāi)害與防治學(xué)報(bào)，2015,27(1):82-86. YANG Dehong，F(xiàn)AN Wen.Zoning of probable occurrence level of geological disastersbased on ArcGIS—a case of Xunyang[J].The Chinese Journal of Geological Hazard and Control，2015.27(1):82-86.

[9] 梁京濤.遙感和GIS 在汶川地震災(zāi)區(qū)地質(zhì)災(zāi)害調(diào)查與評(píng)價(jià)中的應(yīng)用研究———以青川縣為例[D].成都，成都理工大學(xué),2009:30-35. LIANG Jingtao.Research on survey of earthquake triggering geohazards and risk assessiment in Wenchuan earthquake stricken areas using RS and GIS-a case study of Qingchuan county[D].Chengdu，Chengdu University of Technology，2009:30-35.

[10] 祁生文,許強(qiáng),劉春玲,等.汶川地震極重災(zāi)區(qū)地質(zhì)背景及次生斜坡災(zāi)害空間發(fā)育規(guī)律[J].工程地質(zhì)學(xué)報(bào),2009,17(1):39-49. QI Shengwen，XU Qiang，LIU Chunlin，et al.Slope instabilities in the severest disaster areas of 5.12Wenchuan earthquake[J]. Journal of Engineering Geology,2009,17(1):39-49.

[11] 王秀英,聶高眾.汶川MS 8.0 級(jí)地震誘發(fā)崩滑特點(diǎn)及其與地震動(dòng)參數(shù)對(duì)應(yīng)關(guān)系初析[J].巖土工程學(xué)報(bào), 2009,31(9):1378-1383. WANG Xiuying, NIE Gaozhong.Characteristics of landslides induced by Wenchuan Ms 8.0 earthquake and preliminary analysis of their relations with ground motion parameters[J].Chinese Journal of Geotechnical Engineering,2009,31(9):1378-1383.

[12] 莊建琦,崔鵬,葛永剛,等.5.12汶川地震崩塌滑坡分布特征及評(píng)價(jià)因子評(píng)價(jià)———以都江堰至汶川公路沿線為例[J].地質(zhì)科技情報(bào), 2009,28(2):16-22. ZHUANG Jianqi，CUI Peng, GE Yonggang,et al.Distribution characteristics and impact factors assessment of collapses and landslides caused by 5.12 Wenchuan earthquake:taking Dujiangyan-Wenchuan highway as a sample[J].Geological Science and Technology Information, 2009,28(2):16-22.

[13] 徐茂其.地貌要素對(duì)自然邊坡穩(wěn)定性的影響[C].自然邊坡穩(wěn)定性分析暨華鎣山邊坡變形趨勢(shì)研討會(huì),1991. XU Maoqi.Effect of landscape elements on natural slope stability[C].Stability analysis and deformation trends seminar in Huaying Mountain Slope of the natural slope,1991.

[14] Cho S H,Ogata y,Kaneko K. Strain rate dependency of the dynamic fracture process in rock. International journal of Rock Mechanics and Minning Sciences, 2003,40:763-777.

[15] 唐春安,左宇軍,秦泗鳳.汶川地震中邊坡淺表層破裂與拋射模式及其動(dòng)力學(xué)解釋[C].汶川大地震工程震害調(diào)查分析與研究,2009. TANG Chunan,ZUO Yujun,QIN Sifeng.Slope shallow chipped and aerosol propellant mode and its dynamics explanation in Wenchuan earthquake[C].Wenchuan Earthquake Engineering Disaster Investigation Analysis and Research,2009.

[16] 童立強(qiáng),祁生文,劉春玲,等.喜馬拉雅山東南地區(qū)地質(zhì)災(zāi)害發(fā)育規(guī)律初步研究[J].工程地質(zhì)學(xué)報(bào),2007,15(6):721-729. TONG Liqiang,QI Shengwen,LIU Chunlin.Preliminary study of geohazard development patternsin the southeast region of himalaya mountains[J].Journal of Engineering Geology,2007,15(6):721-729.

[17] 黃潤秋,李為樂.5. 12汶川大地震觸發(fā)地質(zhì)災(zāi)害的斷層效應(yīng)研究[J].工程地質(zhì)學(xué)報(bào)，2009, 17 (1):19-28. HUANG Runqiu,LI Weile.Analys is of geo-hazard triggered by wenchaun earthquake[J].Journal of Engineering Geology,2009,17(1):19-28.

[18] 黃潤秋,李為樂.汶川地震觸發(fā)地質(zhì)災(zāi)害的發(fā)育分布規(guī)律研究[J].巖石力學(xué)與工程學(xué)報(bào)，2008, 27 (12)：2585-2592. HUANG Runqiu,LI Weile.Research on development and distribution rules of geohazards induced by Wenchuan earthquake on 12th may,2008[J].Chinese Journal of Rock Mechanics and Engineering,2008, 27 (12)：2585-2592.

[19] 殷躍平.汶川八級(jí)地震滑坡特征分析[J].工程地質(zhì)學(xué)報(bào),2009, 17(1):29-38. YIN Yueping.Features of landslide striggered by the wenchuan earthquake[J].Journal of Engineering Geology,2009, 17(1):29-38.

[20] 喬彥肖,馬中社,呂鳳軍.汶川地震地質(zhì)災(zāi)害發(fā)育特點(diǎn)及動(dòng)因機(jī)制分析[J].中國地質(zhì),2009,36(3):736-741. QIAO Yanxiao, MA Zhongshe, LYU Fengjun Characteristics and dynamic cause mechanism of the Wenchuan earthquake geological hazards[J].Geology in China,2009,36(3):736-741.

[21] 李景富,牛瑞卿.基于GIS的斜坡結(jié)構(gòu)圖自動(dòng)化制圖方法研究[J],人民長(zhǎng)江,2009,40(19):38-40. LI Jingfu, NIU Ruiqing. Research on automatic drawing method of slope structure map based on GIS[J],Yangtze River,2009,40(19):38-40.

[22] 邱丹丹,牛瑞卿,趙艷南,等.斜坡單元支持下地震滑坡危險(xiǎn)性區(qū)劃——以蘆山地震為例[J].吉林大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,45(5):1470-1478. Qiu Dandan，NIU Ruiqing，ZHAO Yannan，et al.Risk zoning of earthquake-induced landslides based on slope units：acase study on lushan earthquake[J].Journal of Jilin University (Earth Science Edition)，2015,45(5)：1470-1478.

[23] 張倬元,王士天,王蘭生,等.工程地質(zhì)分析原理[M].北京，地質(zhì)出版社,2008:112-150. ZHANG Zhuoyuan, WANG Shitian, WANG Lansheng, et al. Principles of engineering geological analysis [M]. Beijing, Geological Publishing House, 2008:112-150.

[24] 祁生文.強(qiáng)震觸發(fā)斜坡失穩(wěn)機(jī)理及動(dòng)力過程研究結(jié)題報(bào)告[R].成都，成都理工大學(xué)，2009:10-17. QI Shengwen. Earthquake triggered slope instability mechanism and dynamic process research report [R].Chengdu, Chengdu University of Technology, 2009:10-17.

[25] 范紅艷,戚鵬程,黃天勇.關(guān)于GIS 中三種地面曲率的探討[J].南陽師范學(xué)院學(xué)報(bào),2011,10(6):66-69. FAN Hongyan, QI Pengcheng, HUANG Tianyong. Discussion on the three kinds of ground curvature in GIS[J]. Journal of Nanyang Normal University, 2011,10 (6): 66-69.

[26] 眭封云.庫區(qū)公路邊坡穩(wěn)定性分析[D].陜西：長(zhǎng)安大學(xué)，2006：32-40. SUI Fengyun.The highway slope stability analysis of reservoir area[D].Shanxi:Chang’an University,2006:32-40.

Risk assessment of seismic geological hazards in small watersheds

LI Peipei,SHEN Junhui,YAN Junsong,CHEN Liang

(StateKeyLaboratoryofGeohazardPreventionandGeoenvironmentProtection(ChengduUniversityofTechnology)，Chengdu,Sichuan610059,China)

This paper takes Yangjia Gully, locating in Beichuan county of Sichuan province, as a typical case of seismic geological hazard risk assessment in small watershed, and selects 7 evaluation factors including slope, exposure, lithology, distance from the fault, slope structure, elevation and slope shape, by analyzing the characteristics of seismic geological hazards. Weights of these factors are also determined by using hierarchy analysis method according to their influential extents. Consequently, a seismic geological hazard assessment map of small watershed is obtained through a comprehensive analysis in GIS. The outcomes indicate that the most influential factor to the forming of seismic geological hazards in Yangjia Gully is the distance from the fault, and slope structure, slope, lithology, elevation, slope shape and exposure, in sequence, are relatively less influential. A prediction of seismic geological hazards is also conducted and demonstrates that the area of high risk in Yangjia Gully comes to 8.98 km2, accounting for 40.43% of the watershed. Besides, comparing with the actual state, the area and quantity of seismic geological hazards in the high risk area account for 70.74% and 66.38% respectively of the total amount, which indicates that the prediction is highly reliable. Thus, this study has a certain significance to promote the research on seismic geological hazard risk assessment in small watershed.

risk assessment; seismic geological hazard; analytical hierarchy process(AHP)；geographical information system(GIS)

10.16031/j.cnki.issn.1003-8035.2017.01.20

2016-07-20;

2016-09-03

高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金(博導(dǎo)類) 資助課題( 20115122110012)；中國地質(zhì)調(diào)查局地質(zhì)調(diào)查子項(xiàng)目(12120113011000)

李佩佩(1991-)，男，甘肅平?jīng)鋈?，在讀碩士研究生，主要從事地質(zhì)災(zāi)害評(píng)價(jià)與防治工作。E-mail:2675715414@qq.com

沈軍輝(1964-)，男，浙江奉化人，博士，教授，博士生導(dǎo)師，主要從事工程地質(zhì)、環(huán)境地質(zhì)研究。E-mail:jhshen2001@sina.com

P694

A

1003-8035(2017)01-0128-07