基于多層次物理力學(xué)參數(shù)的小區(qū)域地震滑坡危險(xiǎn)性評(píng)估——以長江上游石棉縣城及周邊為例

宋 志，倪化勇，周洪福，馮 偉

(1.中國地質(zhì)調(diào)查局成都地質(zhì)調(diào)查中心，成都 610081；2.西南交通大學(xué)土木工程學(xué)院，成都 610031)

?

基于多層次物理力學(xué)參數(shù)的小區(qū)域地震滑坡危險(xiǎn)性評(píng)估
——以長江上游石棉縣城及周邊為例

宋 志1，2，倪化勇1，周洪福1，馮 偉2

(1.中國地質(zhì)調(diào)查局成都地質(zhì)調(diào)查中心，成都 610081；2.西南交通大學(xué)土木工程學(xué)院，成都 610031)

通過實(shí)地調(diào)查、遙感解譯、資料收集等手段，獲取滑坡崩塌體編錄、松散堆積層、地質(zhì)單元的巖土體物理力學(xué)參數(shù)，使得滑坡編錄、地質(zhì)調(diào)查數(shù)據(jù)與區(qū)域Newmark位移模型有機(jī)結(jié)合。研究表明，在滑坡編錄等3個(gè)層次中，由第一層次到第三層次，物理力學(xué)參數(shù)精度逐漸下降，這也反映了滑坡編錄在危險(xiǎn)性評(píng)價(jià)中所占據(jù)的重要性，更能與實(shí)際相吻合。通過對(duì)長江上游石棉縣城地質(zhì)災(zāi)害潛在危險(xiǎn)性的評(píng)估，得出了不同尺度峰值加速度下危險(xiǎn)性分布區(qū)域與規(guī)律，經(jīng)與危險(xiǎn)性線性擬合，在峰值加速度a=0.3時(shí)，區(qū)域危險(xiǎn)區(qū)面積呈大規(guī)模急劇上升，為區(qū)域毀滅性災(zāi)難的臨界值。同時(shí)，石棉縣城隨著峰值加速度數(shù)值增大，危險(xiǎn)區(qū)從滑坡編錄控制逐漸過渡到坡度控制，顯示了多層次物理力學(xué)參數(shù)下危險(xiǎn)性評(píng)估的合理性。

滑坡；危險(xiǎn)性；力學(xué)參數(shù)；滑坡編錄

0 引言

地震滑坡危險(xiǎn)性評(píng)估是國內(nèi)外研究的熱點(diǎn)問題，尤其在近年來地震頻發(fā)的大環(huán)境下，該研究日益引起廣泛重視。地震滑坡危險(xiǎn)性評(píng)估的技術(shù)方法主要分為統(tǒng)計(jì)分析法、綜合指標(biāo)法、確定性評(píng)價(jià)方法。其中，確定性評(píng)價(jià)方法，由于其強(qiáng)大的GIS分析功能、斜坡動(dòng)力的適應(yīng)性，被廣泛運(yùn)用到不同尺度的區(qū)域地震滑坡評(píng)估[1-2]。

在確定性評(píng)價(jià)方法中，最具代表性的模型為Newmark位移模型。葛華等[3]將Newmark位移模型運(yùn)用到映秀震中區(qū)地震滑坡危險(xiǎn)性評(píng)價(jià)，得出了位移圖與滑坡概率分布圖，并開展了編圖示范。王濤等[4]以汶川8.0級(jí)地震為例開展了基于簡化Newmark位移模型的區(qū)域地震滑坡危險(xiǎn)性快速評(píng)估，得出了區(qū)域地震滑坡危險(xiǎn)分區(qū)。這些基于Newmark位移模型的地震滑坡危險(xiǎn)性研究從源頭到結(jié)果采用ArcGIS軟件編制了一系列圖件，取得了較好的認(rèn)識(shí)與結(jié)果。但是也存在一些不足，主要表現(xiàn)為區(qū)域地質(zhì)圖代表巖土體，適用于應(yīng)急快速評(píng)估，而精度不夠，缺乏實(shí)際應(yīng)用；滑坡編錄成果貢獻(xiàn)少，缺乏實(shí)例進(jìn)行相互驗(yàn)證[5-8]。針對(duì)以上不足，本文以巖土體物理力學(xué)參數(shù)為出發(fā)點(diǎn)，從不同精確程度、不同獲取方法，以長江上游石棉縣城及周邊為例，開展大比例尺小區(qū)域的地震地質(zhì)災(zāi)害危險(xiǎn)性評(píng)估，以期為復(fù)雜山區(qū)地震地質(zhì)災(zāi)害的危險(xiǎn)性評(píng)估和防災(zāi)減災(zāi)提供方法技術(shù)。

1 研究區(qū)概況

1.1 研究區(qū)地理位置

長江上游石棉縣城位于大渡河中游地段，為雅安市石棉縣政府所在地，區(qū)域內(nèi)人口密集，經(jīng)濟(jì)較為發(fā)達(dá)，為石棉縣及川西地區(qū)重要的經(jīng)濟(jì)、文化、政治中心之一。研究區(qū)平面上呈矩形，長5.5 km，寬4.4 km，總面積為24.2 km2，涵蓋了石棉縣城的企事業(yè)機(jī)關(guān)單位、人口建筑密集區(qū)，為地質(zhì)災(zāi)害風(fēng)險(xiǎn)的重點(diǎn)防治區(qū)[1-4]。

G108線、省道211線與2012年建成通車的京昆高速公路G5在縣城交匯?？h城經(jīng)G5高速公路向北至雅安120 km，距成都260 km，南距西昌150 km；經(jīng)省道211線至瀘定110 km。研究區(qū)作為G5高速與S211的連接點(diǎn)，交通便利，是雅安市南大門，也是川藏南線和川滇西線的交通要道。

1.2 研究區(qū)地質(zhì)環(huán)境條件

區(qū)內(nèi)氣象冬季溫暖干燥，春末夏初干旱多風(fēng)，夏季悶熱，具四季不分明的亞熱帶季風(fēng)氣候特征。區(qū)內(nèi)多年平均年降水量為803.1 mm，但在時(shí)空上分布不均。降水主要集中在5—9月份，占全年降水量的86.4%；同時(shí)，山地降雨多于河谷地帶，且多以暴雨或陣雨出現(xiàn)。研究區(qū)水系主要為三級(jí)河流——大渡河及其支流楠椏河，大渡河由西至東貫穿全區(qū)，楠椏河由南至北，于縣城處匯入大渡河。

地形地貌上，研究區(qū)屬于緩脊峽谷地貌，地勢西北、東部、西南高，中部及河谷地帶低。最高點(diǎn)位于大坪崗，海拔2030 m，最低點(diǎn)位于區(qū)內(nèi)大渡河最下游汪家崗附近，海拔850 m，相對(duì)高差1180 m，地形起伏度為487.6 m/km2。受河流強(qiáng)烈切割，研究區(qū)地勢陡峭，坡度由河流面0°至陡坡82°范圍分布。據(jù)ArcGIS坡度統(tǒng)計(jì)，區(qū)內(nèi)平均坡度30.8°，達(dá)到陡坡分級(jí)標(biāo)準(zhǔn)，切割強(qiáng)度顯著，頻繁引發(fā)各類地質(zhì)災(zāi)害。

地質(zhì)構(gòu)造上，研究區(qū)位于鮮水河、龍門山、安寧河三大著名斷裂帶的交匯部位，受地質(zhì)構(gòu)造作用影響顯著，呈三面河流深切割地形，尤其是安寧河斷裂帶分支——石棉斷裂直接從研究區(qū)近南北向通過(見圖1)，使得區(qū)內(nèi)巖體破碎，不良工程地質(zhì)現(xiàn)象發(fā)育[9-11]。

圖1 研究區(qū)石棉斷裂分布簡圖Fig.1 Distribution of Shimian fault in the research area

綜上，從孕災(zāi)背景條件上，研究區(qū)位于鮮水河、安寧河、龍門山三大著名斷裂帶的交匯部位，構(gòu)造運(yùn)動(dòng)強(qiáng)烈，歷史地震頻發(fā)。研究區(qū)地質(zhì)環(huán)境條件脆弱，主要表現(xiàn)為：構(gòu)造運(yùn)動(dòng)強(qiáng)烈且石棉斷裂直接從縣城通過，位于暴雨中心常發(fā)生群發(fā)性地質(zhì)災(zāi)害，山體深切割巖體破碎具有滑坡高易發(fā)的孕災(zāi)背景，人類工程活動(dòng)強(qiáng)烈加劇地質(zhì)災(zāi)害的發(fā)生幾率，是地震活動(dòng)頻繁和地質(zhì)災(zāi)害高易發(fā)區(qū)。

2 計(jì)算方法與過程

2.1 Newmark位移模型方法簡述

Newmark位移模型基于極限平衡理論開展累積位移分析，主要的原理是：在地震加速度等耦合作用下，滑坡沿滑動(dòng)面發(fā)生瞬時(shí)位移，并不斷累積，在施加的地震加速度大于斜坡臨界荷載加速度情況下，滑坡便產(chǎn)生，并通過二者加速度的差值進(jìn)行二次積分，得到永久位移值[12]。模型的基本原理示意見圖2。

ac—臨界加速度；a—地面加速度；D—Newmark位移量圖2 Newmark位移模型原理圖Fig.2 Schematic diagram of Newmark displacement model

先利用邊坡靜力學(xué)公式計(jì)算滑坡靜態(tài)安全系數(shù)F，再通過重力加速度g、邊坡坡度β以及安全系數(shù)F間接求出臨界加速度ac，最后通過地震加速度a與臨界加速度ac的差值部分進(jìn)行二次積分，獲取Nermark位移量D，從而進(jìn)行不同尺度下地震滑坡危險(xiǎn)性預(yù)測評(píng)估。在施加臨界加速度ac的過程中，目前研究的方法基本相似，唯一的區(qū)別在于靜力學(xué)平衡方式選取的不同。考慮到研究區(qū)地震時(shí)，處于少雨或無雨季節(jié)，可不考慮地下水的影響，將靜力學(xué)方程進(jìn)行簡化處理。各公式如下：

(1)

(2)

式中，c′—有效內(nèi)聚力，N/m2；—巖土體重度，N/m3；z—破壞面深度，m；β—邊坡坡度傾角，(°)；′—有效內(nèi)摩擦角，(°)；F—靜態(tài)安全系數(shù)，無量綱；g—重力加速度，9.8 m/s2；ac—臨界加速度，m/s2。

Newmark位移計(jì)算過程中需進(jìn)行地震加速度與臨界加速度差值的二次積分，涉及區(qū)域數(shù)據(jù)量大，且隨時(shí)程變化地震參數(shù)難以獲取，難以進(jìn)行有效的地震滑坡區(qū)劃。為此許多學(xué)者采用經(jīng)驗(yàn)公式與簡化地震參數(shù)方法開展計(jì)算。如Jibson[12]通過大量數(shù)據(jù)與實(shí)際滑坡進(jìn)行對(duì)比分析，得出顯著回歸方程，最終驗(yàn)證精確率高，滿足地震滑坡累計(jì)位移量的需要。該方程的優(yōu)點(diǎn)在于地震滑坡位移計(jì)算僅涉及峰值加速度、矩震級(jí)2個(gè)區(qū)域地震參數(shù)，計(jì)算結(jié)果滿足地震滑坡預(yù)測評(píng)估需要。其中，歷史地震的矩震級(jí)Mw利用面波震級(jí)Ms，通過經(jīng)驗(yàn)公式換算得到。具體公式如下：

(3)

Mw=0.844Ms+0.951

(4)

式中，D—Nermark位移，cm；Mw—矩震級(jí)；Ms—面波震級(jí)；a—峰值加速度(可近似代替地面加速度值)，m/s2。

以上基于回歸估計(jì)的Nermark位移經(jīng)驗(yàn)公式，優(yōu)點(diǎn)在于根據(jù)預(yù)測需要將區(qū)域面波震級(jí)、峰值加速度進(jìn)行賦值，便能進(jìn)行不同尺度下地震滑坡危險(xiǎn)性評(píng)估，簡化了流程，可操作性強(qiáng)。

2.2 多層次物理力學(xué)參數(shù)獲取

危險(xiǎn)性評(píng)價(jià)的準(zhǔn)確度取決于參數(shù)獲取的精度，其中，物理力學(xué)參數(shù)是影響滑坡危險(xiǎn)性評(píng)價(jià)的重要因子。由于研究區(qū)為1∶5000高精度地形圖，不能根據(jù)以往傳統(tǒng)小比例尺按1∶50000地質(zhì)圖進(jìn)行巖組劃分，而是采用實(shí)地調(diào)查的巖組進(jìn)行大比例尺細(xì)分。

巖土體的劃分分為3個(gè)層次，第一層次為已調(diào)查的滑坡崩塌體編錄數(shù)據(jù)，第二層次為松散堆積層遙感解譯數(shù)據(jù)，第三層次是區(qū)域地質(zhì)圖件反映的地質(zhì)單元數(shù)據(jù)。在3個(gè)層次中，由第一層次到第三層次，物理力學(xué)參數(shù)精度逐漸下降，這也反映了滑坡編錄在危險(xiǎn)性評(píng)價(jià)所占據(jù)的重要性，更能與實(shí)際相吻合。

研究區(qū)內(nèi)滑坡編錄的數(shù)量為17個(gè)，區(qū)域點(diǎn)密度為0.7個(gè)/km2。松散堆積層的物理力學(xué)指標(biāo)與地質(zhì)單元相差甚大，不以單純的區(qū)域地質(zhì)圖上的地質(zhì)單元確定，而是進(jìn)行不同類別的詳細(xì)劃分。松散堆積層主要分為沖洪積堆積、冰水堆積、殘坡堆積、崩坡堆積、殘坡堆積，區(qū)內(nèi)共劃分22組不同物理力學(xué)指標(biāo)的松散堆積層，編號(hào)為SS01—SS22；地質(zhì)單元?jiǎng)t根據(jù)區(qū)域地質(zhì)圖件直接勾繪，劃分為6組，編號(hào)為DZ01—DZ06；河流不參與參數(shù)劃分(見表1、圖3)。

表1 研究區(qū)物理力學(xué)參數(shù)取值

圖3 研究區(qū)不同物理力學(xué)參數(shù)分布圖Fig.3 The distribution diagram of different physical and mechanical parameters in research area

通過對(duì)研究區(qū)開展大比例尺的不同物理力學(xué)參數(shù)的劃分，涵蓋了滑坡編錄的重點(diǎn)參數(shù)、松散堆積層的次重點(diǎn)參數(shù)以及區(qū)域地質(zhì)單元的次要參數(shù)，將區(qū)域巖土體的劃分有機(jī)組合在一起，有針對(duì)性地開展研究區(qū)滑坡危險(xiǎn)性評(píng)價(jià)。

3 結(jié)果與分析

通過獲取的物理力學(xué)參數(shù)值和坡度值，利用ArcGIS柵格計(jì)算功能，對(duì)研究區(qū)斜坡靜態(tài)安全系數(shù)F按照公式(1)、臨界加速度ac按照公式(2)開展計(jì)算，結(jié)果見圖4、圖5。

根據(jù)上節(jié)所述，對(duì)于預(yù)測性Nermark位移危險(xiǎn)性評(píng)估，地震動(dòng)參數(shù)只需要峰值地震加速度a、面波震級(jí)Ms便能利用Jibson[12]的回歸方程進(jìn)行計(jì)算。根據(jù)《建筑抗震設(shè)計(jì)規(guī)范》[13]以及相關(guān)資料查詢，設(shè)定研究區(qū)的基本地震參數(shù)分別為：①a=0.1g、Ms=7；②a=0.2g、Ms=7；③a=0.3g、Ms=8；④a=0.4g、Ms=9。為了便于簡化，不同尺度以動(dòng)峰值加速度單因子為準(zhǔn)。

通過公式(3)、公式(4)開展Nermark位移值D計(jì)算，并將其值作為不同工況下地震危險(xiǎn)性評(píng)價(jià)依據(jù)。同時(shí)在ArcGIS中采用自然斷點(diǎn)法分為4類，分別代表高、中、低、極低危險(xiǎn)性程度，計(jì)算結(jié)果見圖6。

從圖6可以看出，高危險(xiǎn)區(qū)范圍隨動(dòng)峰值加速度的增長而增大，呈顯著正相關(guān)。利用ArcGIS軟件的柵格重分類功能，得到不同地震工況下高、中危險(xiǎn)區(qū)分布面積(見表2)。

圖4 靜態(tài)安全系數(shù)F圖Fig.4 The image of static safety coefficient

圖5 臨界加速度ac圖Fig.5 The image of critical accelaration

地震工況分布面積/km2高危險(xiǎn)區(qū)中危險(xiǎn)區(qū)a=01g278401a=02g412629a=03g6781002a=04g11621428

對(duì)包括汶川地震在內(nèi)的歷史地震誘發(fā)的滑坡進(jìn)行統(tǒng)計(jì)分析，滑坡數(shù)量規(guī)模隨震級(jí)(動(dòng)峰值加速度)呈顯著正相關(guān)，并且存在一個(gè)危險(xiǎn)區(qū)面積急劇上升的拐點(diǎn)。該拐點(diǎn)是區(qū)域大面積災(zāi)難的臨界值，一旦震級(jí)(或動(dòng)峰值加速度)逼近或達(dá)到該拐點(diǎn)，對(duì)區(qū)域滑坡災(zāi)害將是毀滅性的。

對(duì)表2不同峰值加速度對(duì)應(yīng)的高、中危險(xiǎn)區(qū)面積分布的相關(guān)性分析，得出冪函數(shù)線性擬合曲線。其中，高危險(xiǎn)區(qū)與峰值加速度的R2值為0.9950、中危險(xiǎn)區(qū)與峰值加速度的R2值為0.9966，相關(guān)性均好。高、中危險(xiǎn)區(qū)線性擬合曲線(見圖7)顯示，在峰值加速度a=0.3g時(shí)，區(qū)域危險(xiǎn)區(qū)面積呈大規(guī)模急劇上升，為區(qū)域毀滅性災(zāi)難的臨界值。

圖7 高、中危險(xiǎn)區(qū)與峰值加速度擬合曲線Fig.7 The fitted curves of medium and high risky areas and peak acceleration

3.1a=0.1g危險(xiǎn)區(qū)分布特征

當(dāng)峰值加速度為0.1g時(shí)，全區(qū)高危險(xiǎn)分布面積為2.78 km2，占全區(qū)面積11.49%，中危險(xiǎn)分布面積為4.01 km2，占全區(qū)面積16.57%。調(diào)查區(qū)危險(xiǎn)區(qū)總體受地震影響較弱，危險(xiǎn)區(qū)分布區(qū)域與常規(guī)滑坡危險(xiǎn)區(qū)編錄調(diào)查一致，主要分布于河谷到山地陡坡地段，受單體滑坡影響顯著。地震作用影響小，區(qū)域危險(xiǎn)區(qū)范圍受已有滑坡編錄控制。

3.2a=0.2g危險(xiǎn)區(qū)分布特征

當(dāng)峰值加速度為0.2g時(shí)，全區(qū)高危險(xiǎn)分布面積為4.12 km2，占全區(qū)面積17.02%，中危險(xiǎn)分布面積為6.29 km2，占全區(qū)面積25.99%。危險(xiǎn)區(qū)面積較a=0.1g時(shí)顯著增長，增長率達(dá)48.02%。高危險(xiǎn)區(qū)分布特征表現(xiàn)為：沿河谷地形急變陡坡呈帶分布，高海拔區(qū)域受地形坡度控制顯著，已有滑坡調(diào)查成片分布。

3.3a=0.3g危險(xiǎn)區(qū)分布特征

當(dāng)峰值加速度為0.3g時(shí)，全區(qū)高危險(xiǎn)分布面積為6.78 km2，占全區(qū)面積28.02%，中危險(xiǎn)分布面積為10.02 km2，占全區(qū)面積41.40%。危險(xiǎn)區(qū)面積較a=0.2g時(shí)顯著增長，增長率達(dá)64.56%。表現(xiàn)的特征為已有危險(xiǎn)區(qū)域向四周范圍擴(kuò)散且順坡向增長明顯、坡度巖性控制顯著、危險(xiǎn)區(qū)成群成片發(fā)育等。

3.4a=0.4g危險(xiǎn)區(qū)分布特征

當(dāng)峰值加速度為0.4g時(shí)，全區(qū)高危險(xiǎn)分布面積為11.62 km2，占全區(qū)面積48.02%，中危險(xiǎn)分布面積為14.28 km2，占全區(qū)面積59.01%。危險(xiǎn)區(qū)面積較a=0.3g時(shí)顯著增長，增長率達(dá)71.38%。a=0.4g工況下區(qū)內(nèi)危險(xiǎn)區(qū)面積呈幾何倍數(shù)增長，遠(yuǎn)超低峰值時(shí)增長速度，全區(qū)基本被高、中危險(xiǎn)區(qū)域包圍。區(qū)內(nèi)表現(xiàn)的特征為：危險(xiǎn)區(qū)面積分布急劇上升，環(huán)境承載力達(dá)到極限，山體危險(xiǎn)區(qū)成片發(fā)育，坡度成為危險(xiǎn)性主控因素。石棉縣城及附近區(qū)域危險(xiǎn)區(qū)廣泛發(fā)育，除平緩河谷地帶外，基本為高、中危險(xiǎn)區(qū)。

4 結(jié)論

與以往地震滑坡危險(xiǎn)性評(píng)估相比，本文基于多層次物理力學(xué)參數(shù)下小區(qū)域地質(zhì)災(zāi)害潛在危險(xiǎn)性評(píng)估方法是通過實(shí)地調(diào)查、遙感解譯、資料收集手段，獲取滑坡崩塌體編錄、松散堆積層、地質(zhì)單元的巖土體物理力學(xué)參數(shù)，使得滑坡編錄、地質(zhì)調(diào)查數(shù)據(jù)與區(qū)域Newmark位移模型有機(jī)結(jié)合。在滑坡編錄等3個(gè)層次中，由第一層次到第三層次，物理力學(xué)參數(shù)精度逐漸下降，這也反映了滑坡編錄在危險(xiǎn)性評(píng)價(jià)所占據(jù)的重要性，更能與實(shí)際相吻合。

通過對(duì)長江上游石棉縣城地質(zhì)災(zāi)害潛在危險(xiǎn)性評(píng)估，得出了峰值加速度a分別為0.1g、0.2g、0.3g、0.4g尺度下危險(xiǎn)性分布區(qū)域。當(dāng)a=0.1g時(shí)，研究區(qū)危險(xiǎn)區(qū)受單體滑坡影響顯著，受已有滑坡編錄控制；當(dāng)a=0.2g時(shí)，危險(xiǎn)區(qū)呈現(xiàn)沿河谷地形急變陡坡呈帶分布、高海拔區(qū)域受地形坡度控制顯著、已有滑坡調(diào)查成片分布；當(dāng)a=0.3g時(shí)，已有危險(xiǎn)區(qū)域向四周范圍擴(kuò)散且順坡向增長明顯、坡度巖性控制顯著、危險(xiǎn)區(qū)成群成片發(fā)育；當(dāng)a=0.4g時(shí)，危險(xiǎn)區(qū)面積分布急劇上升，環(huán)境承載力達(dá)到極限、山體危險(xiǎn)區(qū)成片發(fā)育、坡度成為危險(xiǎn)性主控因素。隨著峰值加速度數(shù)量增大，危險(xiǎn)區(qū)從滑坡編錄控制逐漸過渡到坡度控制，顯示了多層次物理力學(xué)參數(shù)下危險(xiǎn)性評(píng)估的合理性。通過對(duì)高、中危險(xiǎn)區(qū)線性擬合曲線分析，在峰值加速度a=0.3g時(shí)，區(qū)域危險(xiǎn)區(qū)面積呈大規(guī)模急劇上升，為區(qū)域毀滅性災(zāi)難的臨界值。

[1] 黃潤秋．汶川地震地質(zhì)災(zāi)害研究[M]．北京：科學(xué)出版社，2009：1～944.

HUANG Run-qiu. Geohazard research on the Wenchuan earthquake[M]. Beijing: Science Press, 2009: 944.

[2] Keefer D K. Investigating landslides caused by earthquakes: A historical review[J]. Surveys in Geophysics, 2002, 23(6): 473～510.

[3] 葛華，陳啟國，王德偉．地震滑坡危險(xiǎn)性評(píng)價(jià)及編圖——以映秀震中區(qū)為例[J]．中國地質(zhì)，2013，40(2)：644～652.

GE Hua, CHEN Qi-guo, WANG De-wei. The assessment and mapping of seismic landslide hazards: A case study of Yingxiu area, Sichuan Province[J]. Geology in China, 2013, 40(2): 644～652.

[4] 王濤，吳樹仁，石菊松，等．基于簡化Newmark 位移模型的區(qū)域地震滑坡危險(xiǎn)性快速評(píng)估——以汶川MS8.0級(jí)地震為例[J]．工程地質(zhì)學(xué)報(bào)，2013，21(1)：16～24.

WANG Tao, WU Shu-ren, SHI Ju-song, et al. Case study on rapid assessment of regional landslide hazard based on simplified Newmark displacement model: WenchuanMs8.0 earthquake[J]. Journal of Engineering Geology, 2013, 21(1): 16～24.

[5] 吳樹仁，王濤，石玲，等．2008 汶川大地震極端滑坡事件初步研究[J]．工程地質(zhì)學(xué)報(bào)，2010，18(2)：145～159.

WU Shu-ren, WANG Tao, SHI Ling, et al. Study on catastrophic landslide triggered by 2008 great Wenchuan earthquake, Sichuan, China[J]. Journal of Engineering Geology, 2010, 18(2): 145～159.

[6] Huang Runqiu, Li Weile. Analysis of the geo-hazards triggered by the 12, May 2008 Wenchuan earthquake, China[J]. Bulletin of Engineering Geology and the Environment, 2009, 68(3): 363～371.

[7] Parker R N, Densmore A L, Rosser N J, et al. Mass wasting triggered by 2008 Wenchuan earthquake is greater than orogenic growth[J]. Nature Geoscience, 2011, 4(7): 449～452.

[8] 王濤．汶川地震重災(zāi)區(qū)地質(zhì)災(zāi)害危險(xiǎn)性評(píng)估研究[D]．北京：中國地質(zhì)科學(xué)院，2010.

WANG Tao. Study on seismic landslide hazard assessment in Wenchuan earthquake severely afflicted area[D]. Beijing: Chinese Academy of Geological Sciences, 2010.

[9] 四川省華地建設(shè)工程有限責(zé)任公司．“4.20”蘆山強(qiáng)烈地震災(zāi)區(qū)雅安市石棉縣地質(zhì)災(zāi)害詳細(xì)調(diào)查報(bào)告[R]．成都：四川省華地建設(shè)工程有限責(zé)任公司，2014.

HuaDi Construction Engineering Limited Liability Company, Sichuan Province. Survey report of geohazard in ‘4.20’ Lushan strong earthquake in Shimian County, Yaan City, Sichuan Province[R]. Chengdu: HuaDi Construction Engineering Limited Liability Company, Sichuan Province, 2014.

[10] 四川省華地建設(shè)工程有限責(zé)任公司．四川省雅安市石棉縣重大地質(zhì)災(zāi)害治理工程復(fù)查復(fù)核報(bào)告[R]．成都：四川省華地建設(shè)工程有限責(zé)任公司，2013.

HuaDi Construction Engineering Limited Liability Company, Sichuan Province. The rechecks report of the regulation project of serious geo hazards in Shimian, Yaan City, Sichuan Province[R]. Chengdu: HuaDi Construction Engineering Limited Liability Company, Sichuan Province, 2013.

[11] 中國地質(zhì)調(diào)查局成都地質(zhì)調(diào)查中心．四川省石棉縣地質(zhì)災(zāi)害詳細(xì)調(diào)查[R]．成都：中國地質(zhì)調(diào)查局成都地質(zhì)調(diào)查中心，2009.

Chengdu Institute of Geology and Mineral Resources. Detailed survey of the geo hazard in Shimian, Sichuan Province[R]. Chengdu: Chengdu Institute of Geology and Mineral Resources, 2009.

[12] Jibson R W. Regression methods for estimating coseismic landslidedisplacement [J]. Engineering Geology, 2007, 90：209～218.

[13] 中華人民共和國建設(shè)部．建筑抗震設(shè)計(jì)規(guī)范(GB 50011-2010)[S]．北京：中國建筑工業(yè)出版社，2010.

Ministry of Construction of the People′s Republc of China. Code for seismic design of building (GB 50011-2010). Beijing: China Building Industry Press, 2010.

RISK ASSESSMENT OF SEISMIC LANDSLIDE WITHIN SMALL REGION BASED ON MULTI-LEVEL PHYSICAL AND MECHANICAL PARAMETERS: A CASE STUDY OF SHIMIAN AND ADJACENT AREAS IN THE UPPER REACHES OF YANGTZE RIVER

SONG Zhi1 2, NI Hua-yong1, ZHOU Hong-fu1, FENG Wei2

(1.Chengdu Institute of Geology and Mineral Resources, Chengdu 610008, China;2.SchoolofCivilEngineering,SouthwestJiaotongUniversity,Chengdu610031,China)

Compared to the risk assessments of seismic landslide before, this paper achieves the landslide and collapse inventory, loose accumulated layer and the physical and mechanical parameters of rock soil mass of geological units by field survey, remote sensing interpretation and data collection. It makes an organic integration of landslide inventory, geological survey and regional Newmark displacement model. During three degrees of the landslide inventory, it can been seen that the physical and mechanical parameters get smaller gradually from the first degree to the third degree, which also inflects the importance of the landslide inventory for risk assessment and that it accords better with reality. The dangerous areas and regular under different peak acceleration are got by assessing the potential risk of geological hazard in Shimian City on the upper Yangzi River. Through the linear fitting of risks, when peak acceleration is 0.3, the dangerous areas increase greatly on the large scale that is the critical value of the regional devastating disaster. Meanwhile, the dangerous area gradually transit from that under the control of landslide inventory to that of slope, which shows the validity of the risk assessment under multi-level physical and mechanical parameters.

landslide; risk; mechanical parameters; landslide inventory

1006-6616(2016)03-0760-11

2016-04-09

中國地質(zhì)調(diào)查局“成渝經(jīng)濟(jì)區(qū)宜賓—萬州沿江發(fā)展帶1∶5萬環(huán)境地質(zhì)調(diào)查”項(xiàng)目(121201010000150015)

宋志(1982-)，男，四川省宣漢縣人，高級(jí)工程師，主要從事巖土工程、地質(zhì)災(zāi)害研究。E-mail：35842126@qq.com

P642.22

A