張家坪泥石流物源活動(dòng)性演化及沖出范圍模擬研究

常鳴　楊濤　周宇　姚成

摘要：為了研究泥石流內(nèi)物源活動(dòng)性演化及沖出范圍，利用搜集的四期遙感影像和DEM，結(jié)合GIS與RS，對(duì)不同時(shí)期張家坪泥石流物源做出詳細(xì)解譯，得出不同活動(dòng)率條件下的張家坪泥石流物源演化趨勢(shì)。根據(jù)變化趨勢(shì)預(yù)測(cè)該泥石流未來(lái)不會(huì)產(chǎn)生較多新的物源，但老崩滑體有擴(kuò)大的可能性，同時(shí)利用FLO-2D軟件模擬三種降雨頻率下的張家坪泥石流沖出范圍及堵江程度，研究結(jié)果能夠?yàn)閺埣移耗嗍魑镌捶乐未胧┨峁﹨⒖純r(jià)值，同時(shí)為沖出范圍內(nèi)的災(zāi)后重建工作提供科學(xué)指導(dǎo)依據(jù)。

關(guān)鍵詞：張家坪;泥石流;遙感解譯;物源演化;沖出量

中圖分類號(hào)：P642文獻(xiàn)標(biāo)志碼：A

Abstract：The detailed descriptions of the Zhangjiaping debris-flow source at different periods and its evolution process under different activity rates were investigated in this study with four remote sensing images，DEM，and GIS and RS techniques.It is predicted that the debris-flow will not generate more new sources in the future，but the existing collapses has a possibility of expansion.FLO-2D software was also used to simulate the flushing ranges and the degree of river blocking under three rainfall frequencies for the Zhangjiaping debris-flow.The research results can provide reference for the source prevention and control measures at Zhangjiaping debris-flow，and provide scientific guidance for post-disaster reconstruction within the runout zone.

Key words：Zhangjiaping;debris-flow;remote interpretation;material evolution;run-out volume

汶川地震后，陳曉清等[1]基于衛(wèi)星影像數(shù)據(jù)在重災(zāi)區(qū)對(duì)自然災(zāi)害體進(jìn)行解譯、統(tǒng)計(jì)，得出崩塌、滑坡災(zāi)害面積達(dá)2 264.53 km2，初步估算汶川震區(qū)地質(zhì)災(zāi)害的水土流失量達(dá)55.86億t。新增松散固體物源產(chǎn)生了大量的次生地質(zhì)災(zāi)害，嚴(yán)重威脅災(zāi)區(qū)居民的生命財(cái)產(chǎn)安全，震后次生災(zāi)害將進(jìn)入活躍期，泥石流的活躍期將持續(xù)10～20 a[2-4]。目前汶川地震已過(guò)去10年，其物源已發(fā)生不同強(qiáng)度的變化，但地質(zhì)災(zāi)害的危險(xiǎn)性仍然存在。唐川等[5]以汶川地震高烈度區(qū)的北川縣城及湔江河谷為研究區(qū)，對(duì)“9·24”暴雨前后利用高精度影像進(jìn)行遙感解譯結(jié)果進(jìn)行對(duì)比研究，得出暴雨前后的滑坡和泥石流變化情況。常鳴等[6]對(duì)綿遠(yuǎn)河流域的20條泥石流溝進(jìn)行解譯，利用GIS統(tǒng)計(jì)出崩塌、滑坡在坡度、坡向、高程、溝道縱坡降、地層巖性、溝壑密度因子上的分布特點(diǎn)。李凌婧[7]利用遙感技術(shù)對(duì)綿遠(yuǎn)河地區(qū)的泥石流進(jìn)行解譯，并對(duì)綿遠(yuǎn)河地質(zhì)災(zāi)害發(fā)育規(guī)律的空間分析。周縱橫等[8]對(duì)綿遠(yuǎn)河地區(qū)的降雨時(shí)間和空間分布情況以及泥石流的形態(tài)、溝道縱坡降、物源儲(chǔ)量、物源密度等方面進(jìn)行調(diào)查。后來(lái)蔣志林[9]、常鳴等[10]利用遙感與GIS技術(shù)先后對(duì)映秀地區(qū)和龍池地區(qū)泥石流流域內(nèi)的崩滑體物源做出面積變化分析。喬建平等[11]對(duì)汶川災(zāi)區(qū)泥石流物源的主要類型進(jìn)行分類，并根據(jù)各自的特點(diǎn)建立了泥石流物源動(dòng)儲(chǔ)量的啟動(dòng)地質(zhì)模式。劉洋等[12]選取8個(gè)影響物源分布的因子利用層次分析-信息量法對(duì)都江堰龍池地區(qū)的泥石流進(jìn)行敏感性分析。龍溪河流域發(fā)生多處泥石流災(zāi)害后，馬煜等[13]對(duì)都江堰龍溪河沿岸泥石流災(zāi)害分布規(guī)律、流域面積、高程、距斷層距離、坡度、巖性等進(jìn)行研究。程霄[14]利用遙感和GIS技術(shù)對(duì)汶川縣映秀地區(qū)的崩滑體變化量敏感性進(jìn)行評(píng)價(jià)。雖然前人已對(duì)強(qiáng)震區(qū)的物源變化已做出研究，但在研究中并未考慮其崩滑體的活動(dòng)性，古崩滑體仍然具有增強(qiáng)危險(xiǎn)的可能。唐晨曉等[15]利用強(qiáng)震區(qū)都江堰龍池鎮(zhèn)多期影像遙感解譯（2008年、2009年、2011年、2013年、2015年），觀察崩滑體面積演化過(guò)程的同時(shí)，根據(jù)每期解譯結(jié)果研究滑體的活動(dòng)狀態(tài)。汶川8.0級(jí)地震誘發(fā)了大量的崩滑體，在余震及后續(xù)的強(qiáng)降雨作用下進(jìn)一步增強(qiáng)其活動(dòng)性[16-17]。在此基礎(chǔ)上，本文以張家坪溝泥石流流域?yàn)檠芯繀^(qū)，研究在基于不同活動(dòng)率條件下的崩滑體演化結(jié)果并對(duì)張家坪溝的物源活動(dòng)性做出預(yù)測(cè)分析，使強(qiáng)震區(qū)的地質(zhì)災(zāi)害研究更為具體化。此外，對(duì)于泥石流的沖出量研究有理論模型和數(shù)值模擬等手段，本文利用FLO-2D軟件模擬在不同降雨頻率條件下的泥石流暴發(fā)情況，得出泥石流的沖出規(guī)模和危險(xiǎn)性分區(qū)。本文的研究成果可為現(xiàn)階段張家坪的防災(zāi)減災(zāi)措施提供一定的參考價(jià)值，同時(shí)為震區(qū)多年的物源變化研究提供科學(xué)參考。

1 研究區(qū)概況

張家坪泥石流位于岷江右岸，溝口坐標(biāo)為103°28′31.61″E， 31°02′59.19″N，隸屬于阿壩州汶川縣映秀鎮(zhèn)張家坪村。張家坪溝流域平面形態(tài)呈不規(guī)則的菱形，最大縱向長(zhǎng)度2.2 km，平均寬度0.9 km，總面積為1.71 km2。整個(gè)流域?qū)儆谏钋蓄~割構(gòu)造侵蝕的低-中山地貌，地勢(shì)總體為西北高，東南低，流域最高點(diǎn)高程為1 995 m，最低點(diǎn)在溝口處高程為880 m，最大相對(duì)高差達(dá)1 115 m。根據(jù)地形數(shù)據(jù)，張家坪溝上游段的縱坡降達(dá)750‰，在溝口附近為185‰，全流域平均縱坡降為386‰。整個(gè)流域內(nèi)有5條具有溝槽的次級(jí)支溝，縱坡降均較大，且表現(xiàn)出上游陡下游漸緩的特點(diǎn)，這種地形條件下，流域內(nèi)的輸沙能力較強(qiáng)，為松散固體物質(zhì)的搬運(yùn)和參與泥石流活動(dòng)均提供了有利的地形條件。

張家坪溝位于強(qiáng)震區(qū)，距離汶川地震震中心較近，在地震后曾發(fā)生兩次規(guī)模較大的泥石流災(zāi)害事件：2010年“8·13”泥石流和2013 年“7·10”泥石流。根據(jù)映秀氣象站發(fā)布的降雨資料，2010年8月13日的一整天，共計(jì)降雨量達(dá)126.8 mm。從8月12日下午17時(shí)到8月14日凌晨2時(shí)的33個(gè)小時(shí)內(nèi)，該地區(qū)累計(jì)降雨量為162.1 mm。在2013年“7·10” 降雨期間，映秀鎮(zhèn)的降雨總量達(dá)576 mm。強(qiáng)降雨過(guò)程導(dǎo)致流域內(nèi)松散物源處于飽和狀態(tài)，在不利條件下極易啟動(dòng)暴發(fā)泥石流災(zāi)害事件，泥石流易造成岷江被堵塞，形成堰塞湖，造成河水雍高，進(jìn)而淹沒(méi)道路及房屋，威脅居民及附近人口的生命財(cái)產(chǎn)安全。

根據(jù)搜集到的張家坪泥石流勘察資料，張家坪溝存在固體物源總量達(dá)92.57萬(wàn)m2，可能參與泥石流活動(dòng)的動(dòng)儲(chǔ)量約31.95萬(wàn)m2，物源十分豐富，本文利用搜集的各期高精度的遙感影像綜合分析張家坪溝流域的物源分布及活動(dòng)率變化情況。

2 遙感解譯及活動(dòng)率賦值

2.1 數(shù)據(jù)來(lái)源

本文搜集到研究區(qū)的四期高精度遙感影像，主要包括2009年2月10日、2011年4月26日、2014年12月19日及2017年5月1日分辨率為0.5 m的衛(wèi)星影像，對(duì)搜集的遙感影像進(jìn)行圖像增強(qiáng)、幾何糾正等處理工作后，這四期影像清晰、色彩飽和，能清晰識(shí)別地物特征、判定地質(zhì)災(zāi)害特性。進(jìn)行精確配準(zhǔn)后利用ARCGIS軟件針對(duì)不同時(shí)期的崩滑體災(zāi)害進(jìn)行詳細(xì)解譯、分析。為簡(jiǎn)化文本，文中數(shù)據(jù)分析時(shí)僅用年代表影像時(shí)期進(jìn)行研究描述。

2.2 遙感解譯及活動(dòng)率分級(jí)

根據(jù)四期高精度遙感影像，由地表植被情況及崩滑體的滑動(dòng)痕跡及裸露和堆積狀態(tài)，開(kāi)展崩滑體的遙感解譯及活動(dòng)率分析。首先對(duì)地震后2009年的遙感影像進(jìn)行解譯，其次在2009年解譯的崩滑體基礎(chǔ)上對(duì)2011年的影像進(jìn)行解譯，根據(jù)其活動(dòng)率對(duì)崩滑體的屬性賦值，依次均在前一期影像的基礎(chǔ)上遙感解譯并進(jìn)行活動(dòng)率賦值。本文將崩滑體活動(dòng)率劃分為如下4級(jí)標(biāo)準(zhǔn)。

（1）極高活動(dòng)率。因缺乏地震前后的數(shù)據(jù)，第一期（2009年）影像則根據(jù)現(xiàn)實(shí)際解譯情況，將其歸類為極高活動(dòng)率;此外，在之后多期遙感影像中由降雨引起擴(kuò)大的崩滑體及新增的崩滑體也定義為極高活動(dòng)率。

（2）高活動(dòng)率。在上一期解譯的崩滑體在本期影像中活動(dòng)面積大于崩滑體面積40%以上。

（3）中活動(dòng)率。在上一期解譯的崩滑體在本期影像中活動(dòng)面積占崩滑體面積10%～40%。

（4）低活動(dòng)率。在上一期解譯的崩滑體在本期影像中活動(dòng)面積小于崩滑體面積10%，或者崩滑體在本期影像中沒(méi)有活動(dòng)跡象。

3 解譯結(jié)果統(tǒng)計(jì)分析

3.1 張家坪溝2009年遙感影像解譯結(jié)果

張家坪溝屬于汶川8.0級(jí)地震影響較為嚴(yán)重的區(qū)域，產(chǎn)生的自然災(zāi)害也較多，根據(jù)統(tǒng)計(jì)及遙感解譯2009年的崩滑體，得出張家坪溝共產(chǎn)生79個(gè)崩滑體，崩滑體面積共15.64萬(wàn)m2。張家坪溝流域面積共1.71 km2，則產(chǎn)生的崩滑體面積占流域面積的10.93%。物源量較大，極有暴發(fā)大規(guī)模泥石流的可能，見(jiàn)圖1。

3.2 張家坪溝2011年遙感影像解譯結(jié)果

基于2009年影像解譯的崩滑體，在2011年影像上再進(jìn)一步解譯。在此間，2010年8月13日該區(qū)域暴發(fā)泥石流災(zāi)害時(shí)間。經(jīng)解譯及統(tǒng)計(jì)，共有150個(gè)崩滑體，較上一期新增有71處，崩滑體總面積共19.26萬(wàn)m2，占流域面積的11.26%，與2009年相比面積增加23.12%。流域內(nèi)新增崩滑體面積3.62萬(wàn)m2，主要是由余震或地震使斜坡達(dá)到極限狀態(tài)后，在強(qiáng)降雨作用下產(chǎn)生較多的新崩滑體或者在原崩滑體基礎(chǔ)上進(jìn)一步擴(kuò)大，根據(jù)崩滑體活動(dòng)率劃分標(biāo)準(zhǔn)將其全部歸為極高。同時(shí)，由2009年影像解譯的79個(gè)古崩滑體中，活動(dòng)率為高的面積為14.4萬(wàn)m2，占總面積的74.76%;活動(dòng)率為中度的面積為0.97萬(wàn)m2，占總面積的5.05%;活動(dòng)率為低的面積為0.37萬(wàn)m2，占總面積的1.62%。由統(tǒng)計(jì)可知在2009-2011年間，張家坪溝流域內(nèi)的崩滑體活動(dòng)率較高，主要原因是在2010年8月間映秀鎮(zhèn)內(nèi)暴發(fā)持續(xù)性降雨，使大部分崩滑體達(dá)到飽和程度，在降雨條件下，崩滑體向不同方向擴(kuò)大，部分發(fā)生垮塌，崩滑體活動(dòng)率為低的面積所占比例最小。

3.3 張家坪溝2014年遙感影像解譯結(jié)果

基于2009年、2011年解譯的崩滑體，在2014年影像上再進(jìn)一步解譯。在此間，2013年7月10日該流域內(nèi)也暴發(fā)泥石流災(zāi)害事件。經(jīng)遙感解譯后統(tǒng)計(jì)，共有163個(gè)崩滑體，較上一期新增有13處，崩滑體總面積共19.96萬(wàn)m2，占流域面積的11.70%，較上一期增加0.44%。流域內(nèi)新增崩滑體面積為0.67萬(wàn)m2，占總面積的3.37%，主要是由降雨產(chǎn)生的新崩滑體或原崩滑體擴(kuò)大的區(qū)域，根據(jù)崩滑體活動(dòng)率劃分標(biāo)準(zhǔn)將其全部歸為極高;此外，由2009年、2011年影像解譯的150個(gè)古崩滑體中，活動(dòng)率為高的面積為10.29萬(wàn)m2，占總面積的51.52%;活動(dòng)率為中度的面積為4.27萬(wàn)m2，占總面積的21.40%;活動(dòng)率為低的面積為4.73萬(wàn)m2，占總面積的23.71%。由此統(tǒng)計(jì)可知在2011-2014年間，流域內(nèi)崩滑體活動(dòng)率為極高的面積所占的比例最小，活動(dòng)率為高的面積所占比例仍是最大。主要原因是在2013年7月間，該流域內(nèi)也暴發(fā)持續(xù)性強(qiáng)降雨，流域內(nèi)已經(jīng)產(chǎn)生大量的崩滑體，古崩滑體再次被激活，作為補(bǔ)充物源形成泥石流災(zāi)害事件，見(jiàn)圖1。

3.4 張家坪溝2017年遙感影像解譯結(jié)果

基于2009年、2011年、2014年影像解譯的崩滑體，在2017年影像上進(jìn)一步解譯。經(jīng)統(tǒng)計(jì)，共有167個(gè)崩滑體，較上一期新增有4處，崩滑體總面積共20.26萬(wàn)m2，占流域面積的11.85%。流域內(nèi)新增崩滑體面積為0.26萬(wàn)m2，占總面積的1.3%，在此階段新增的崩滑體面積最少;此外，由2010年、2011年及2014年影像解譯的163個(gè)古崩滑體中，活動(dòng)率為高的面積為7.47萬(wàn)m2，占總面積的36.9%;活動(dòng)率為中度的面積為6.90萬(wàn)m2，占總面積的34.1%;活動(dòng)率為低的面積為5.59萬(wàn)m2，占總面積的27.6%。由此統(tǒng)計(jì)可知在2014-2017年間，流域內(nèi)在沒(méi)有強(qiáng)降雨等外動(dòng)力下，新產(chǎn)生的崩滑體較少，主要是古崩滑體在持續(xù)活動(dòng)，并逐漸呈現(xiàn)衰減趨勢(shì)。主要在經(jīng)過(guò)2014-2017幾年間山體的自然恢復(fù)作用及未發(fā)生持續(xù)性暴雨等外動(dòng)力，植被恢復(fù)率逐漸增加，崩滑體活動(dòng)率逐漸降低，面積相對(duì)較小的崩滑體擾動(dòng)小，恢復(fù)穩(wěn)定的速度更快，見(jiàn)表1。

3.5 演化規(guī)律分析經(jīng)過(guò)上統(tǒng)計(jì)結(jié)果，得出基于多期遙感影像的不同活動(dòng)率崩滑體的演化圖，詳見(jiàn)圖2。以黃潤(rùn)秋[18]收集震前、震后的大量地質(zhì)災(zāi)害作為基礎(chǔ)數(shù)據(jù)，對(duì)地質(zhì)災(zāi)害的特點(diǎn)和規(guī)律進(jìn)行分析，得出地質(zhì)災(zāi)害將以4～5 a為一個(gè)高峰周期，呈震蕩式的下降，并最終恢復(fù)到震前的水平。結(jié)合張家坪溝整理的震后基礎(chǔ)數(shù)據(jù)，在不考慮持續(xù)性暴雨及地震等強(qiáng)大外動(dòng)力條件下，對(duì)張家坪流域內(nèi)的崩滑體活動(dòng)率從2009年到2017年的變化及趨勢(shì)做出分析，見(jiàn)圖2。

極高活動(dòng)率的崩滑體在2009-2017年間的面積逐漸減小，由2014年及2017年的數(shù)據(jù)可得張家坪溝的新增面積已經(jīng)降至最低，在未來(lái)將不會(huì)產(chǎn)生較多的新崩滑體，但古崩滑體的活動(dòng)性有再次增強(qiáng)的可能。高活動(dòng)率的崩滑體在2011-2017年間面積逐漸減小，預(yù)測(cè)至少需要5 a面積將會(huì)恢復(fù)至中度及以下?tīng)顟B(tài)。中活動(dòng)率的崩滑體面積在2011-2017年間面積逐漸增加，未來(lái)將進(jìn)行震蕩式的波動(dòng)下降，及活動(dòng)率為高的崩滑體恢復(fù)至中度狀態(tài)，預(yù)測(cè)還需要7～10 a恢復(fù)至穩(wěn)定狀態(tài)。低活動(dòng)率的崩滑體在2011-2017年間的面積逐年增加，但之間也會(huì)有震蕩式的平穩(wěn)恢復(fù)期，根據(jù)變化趨勢(shì)預(yù)測(cè)未來(lái)將會(huì)增加更多的面積，但也會(huì)有部分崩滑體會(huì)有震蕩式的活動(dòng)，直至最后的植被逐漸恢復(fù)，至少還需要10 a時(shí)間山體才能達(dá)到穩(wěn)定狀態(tài)。但現(xiàn)階段仍有滑體處于活動(dòng)狀態(tài)，而且流域的地形條件在經(jīng)過(guò)洪水、泥石流的沖刷、下切后更有利于水動(dòng)力條件的形成和物源的起動(dòng)，可導(dǎo)致激發(fā)泥石流的臨界雨強(qiáng)降低，再次發(fā)生泥石流的可能性仍然較大，仍然需要加強(qiáng)泥石流的工程防治、預(yù)警預(yù)報(bào)、危險(xiǎn)范圍分區(qū)等工作。

4 不同降雨頻率下張家坪泥石流沖出范圍

張家坪流域面積為1.71 km2，主溝長(zhǎng)度1.8 km，平均縱坡降為317‰，且張家坪溝流域內(nèi)的物源量巨大，在這樣的地形條件下，其活動(dòng)性很強(qiáng)，極易暴發(fā)泥石流災(zāi)害事件。前文通過(guò)遙感影像結(jié)合GIS技術(shù)分析其物源的活動(dòng)性，在此條件下，利用DEM數(shù)據(jù)，運(yùn)用FLO-2D軟件模擬不同頻率條件下的泥石流沖出規(guī)模，達(dá)到高效的防災(zāi)減災(zāi)技術(shù)。

4.1 FLO-2D的基本原理及參數(shù)選取

FLO-2D是1988年O′Brien提出的基于非牛頓流體模式及中央有限差分的數(shù)值模擬軟件。該軟件在洪水災(zāi)害管理、城市淹沒(méi)分析、泥石流災(zāi)害危險(xiǎn)性劃分等方面有廣泛應(yīng)用[19-20]。賈濤[21]等結(jié)合FLO-2D軟件和GIS技術(shù)以泥石流沖量進(jìn)行分區(qū)，建立泥石流堆積扇的危險(xiǎn)性分區(qū)模型。黃勛[22]等利用FLO-2D軟件構(gòu)建了一套適用于我國(guó)西南山地城鎮(zhèn)的泥石流定量風(fēng)險(xiǎn)評(píng)價(jià)的理論體系和技術(shù)流程。梁鴻熙[23]利用FLO-2D分析泥石流流動(dòng)及堆積特性與黏性系數(shù)和屈服應(yīng)力的關(guān)系。龔柯[24]等利用FLO-2D軟件模擬汶川縣綿虒鎮(zhèn)地區(qū)在P=2%條件下發(fā)生的泥石流，并構(gòu)建泥石流危險(xiǎn)性評(píng)價(jià)模型。常鳴等[25]已運(yùn)用FLO-2D軟件模擬汶川震區(qū)多條泥石流溝，其技術(shù)已非常成熟。開(kāi)展數(shù)值模擬的主要參數(shù)有：泥石流流量、體積濃度（CV）、層流阻滯系數(shù)（K）、曼寧系數(shù)（n）等。采用雨洪法計(jì)算出不同降雨頻率下的泥石流流量，其它參數(shù)可結(jié)合張家坪泥石流的實(shí)際情況進(jìn)行取值，見(jiàn)表2。

4.2 張家坪泥石流模擬沖出范圍

根據(jù)參數(shù)得出不同降雨頻率下的模擬結(jié)果見(jiàn)圖3至圖5，不同降雨頻率下的沖出范圍見(jiàn)圖6。由模擬結(jié)果得：在降雨頻率為P=5%（20年一遇）的條件下，泥石流堆積扇面積0.73 萬(wàn)m2。泥石流沖出溝口進(jìn)入岷江，堆積扇平均堆積厚度1m左右，溝道堆積物大多處于溝道中下游。在降雨頻率為P=2%（50年一遇）的條件下，泥石流堆積扇面積1.52

萬(wàn)m2。泥石流沖入岷江，堆積扇平均堆積厚度在3 m左右，溝道堆積物大多

處于溝道中下游，溝道平均堆積厚度為2 m左右。在降雨頻率為P=1%（100年一遇）的條件下，大量松散物質(zhì)沖入岷江，擠壓河道，迫使水位抬高。堆積扇面積為2.29 萬(wàn)m2，平均堆積厚度4～5 m。

5 結(jié)論

本文通過(guò)多期遙感影像和DEM數(shù)據(jù)，開(kāi)展物源解譯及活動(dòng)性評(píng)價(jià)工作，統(tǒng)計(jì)分析出張家坪溝泥石流基于不同活動(dòng)率條件下的物源演化結(jié)果，同時(shí)運(yùn)用FLO-2D軟件模擬在不同降雨頻率下的泥石流沖出范圍。

（1）在不考慮持續(xù)性暴雨及地震等強(qiáng)大外動(dòng)力條件下，極高活動(dòng)率的崩滑體在2009-2017年間的面積逐漸減小，崩滑體的新增面積已經(jīng)降至最低，在未來(lái)將不會(huì)產(chǎn)生較多的新崩滑體，但老崩滑體的擴(kuò)大的可能性很高;高活動(dòng)率的崩滑體在2011-2017年間面積逐漸減小，預(yù)測(cè)至少需要5 a會(huì)恢復(fù)至中度及以下?tīng)顟B(tài);中等活動(dòng)率和低活動(dòng)率的崩滑體面積在2011-2017年間面積逐漸增加;根據(jù)變化趨勢(shì)預(yù)測(cè)至少還需要10 a時(shí)間山體才能達(dá)到穩(wěn)定狀態(tài)。現(xiàn)階段仍有大部分崩滑體處于活動(dòng)狀態(tài)，再次發(fā)生泥石流的可能性仍然存在，需要加強(qiáng)張家坪泥石流的監(jiān)測(cè)預(yù)警工作。

（2）運(yùn)用FLO-2D模擬張家坪泥石流不同降雨頻率下的沖出范圍及形成不同堵江程度。在降雨頻率為P=5%（20年一遇）、2%（50年一遇）、1%（100年一遇）的條件下，泥石流沖出范圍分別為0.73萬(wàn)m2、1.52萬(wàn)m2、2.29萬(wàn)m2，平均堆積厚度分別為1 m、3 m、4～5 m，存在不同程度的堵江現(xiàn)象，仍需對(duì)張家坪泥石流開(kāi)展長(zhǎng)期監(jiān)測(cè)工作。

參考文獻(xiàn)（References）：

[1] 陳曉清，李智廣，崔鵬，等.5.12汶川地震重災(zāi)區(qū)水土流失初步估算[J].山地學(xué)報(bào)，2009，27（1）：122-127.（CHEN X Q，LI Z G，CUI P，et al.Estimation of soil erosion caused by the 5·12 Wenchuan Earthquake[J].Journal of Mountain Science，2009.27（1）：122-127.（in Chinese）） DOI：10.16089/j.cnki.1008-2786.2009.01.018.

[2] CUI P，CHEN X Q，ZHU Y Y，et al.The Wenchuan Earthquake （May 12，2008），Sichuan Province，China，and resulting geohazards[J].Natural Hazards，2011，56（1）：19-36.DOI：10.1007/s11069-009-9392-1.

[3] TANG C，ZHU J，LI W L，et al.Rainfall-triggered debris flows following the Wenchuan earthquake[J].Bulletin of Engineering Geology and the Environment，2009，68：187-194.DOI：10.1007/s10064-009-0201-6.

[4] 唐川.汶川地震區(qū)暴雨滑坡泥石流活動(dòng)趨勢(shì)預(yù)測(cè)[J].山地學(xué)報(bào)，2010，28（3）：341-349.（TANG C.Activity tendency prediction of rainfall induced landslides and debris flows in the Wenchuan Earthquake Areas[J].Journal of Mountain Science，2010，28（3）：341-349.（in Chinese）） DOI：10.16089/j.cnki.1008-2786.2010.03.009.

[5] 唐川，丁軍，齊信，等.汶川地震高烈度區(qū)暴雨滑坡活動(dòng)的遙感動(dòng)態(tài)分析[J].地球科學(xué)-中國(guó)地質(zhì)大學(xué)學(xué)報(bào)，2010，35（2）：317-323.（TANG C，DING J，QI X，et al.Dynamic analysis on rainfall-induced landslide activity in high seismic intensity areas of the Wenchuan Earthquake using remote sensing image[J].Earth Science-Journal of China University of Geosciences，2010，35（2）：317-323.（in Chinese）） DOI：10.3799/dqkx.2010.033.

[6] 常鳴，唐川，李為樂(lè)，等.汶川地震區(qū)綿遠(yuǎn)河流域泥石流形成區(qū)的崩塌滑坡特征[J].山地學(xué)報(bào)，2012，30（5）：561-569.（CHANG M，TANG C，LI W L，et al.The characteristic of collapse and landslide by Wenchuan Earthquake in debris flow for mative region along the Mianyuan River Basin，China[J].Journal of Mountain Science，2012，30（5）：561-569.（in Chinese）） DOI：10.16089/j.cnki.1008-2786.2012.05.017.

[7] 李凌婧，姚鑫，張永雙，等.汶川地震綿遠(yuǎn)河流域地質(zhì)災(zāi)害遙感提取與分布特征研究[J].工程地質(zhì)學(xué)報(bào)，2014，22（1）：46-55.（LI L J，YAO X，ZHANG Y S，et al.Rs-based extraction and distribution characteristics of geo-hazards triggered by Wenchuan Earthquake in Mianyuan River[J].Journal of Engineering Geology，2014，22（1）：46-55.（in Chinese）） DOI：10.13544/j.cnki.jeg.2014.01.015.

[8] 周縱橫，任光明，許英杰，等.汶川震區(qū)清平鄉(xiāng)綿遠(yuǎn)河流域泥石流發(fā)育特征分析[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào)，2017，28（1）：30-35.（ZHOU Z H，REN G M，XU Y J，et al.Characteristics of the debris flow in Qingping country of Mianyuan river basin，epicenter Wenchuan earthquake[J].Chinese Journal of Geological Hazard and Control，2017，28（1）：30-35.（in Chinese）） DOI：10.16031/j.cnki.issn.1003-8035.2017.01.05.

[9] 蔣志林.基于RS的震區(qū)泥石流物源演化特征分析——以四川省汶川縣映秀地區(qū)為例[D].成都：成都理工大學(xué)，2014.（JIANG Z L.Research on the materials evolution characteristics of debris flows based on RS in earthquake area–Case study in Yingxiu Area of Wenchuan in Sichuan province[D].Chengdu：Chengdu University of Technology，2014.（in Chinese））

[10] [ZK（#]常鳴，唐川，蔣志林，等.強(qiáng)震區(qū)都江堰市龍池鎮(zhèn)泥石流物源的遙感動(dòng)態(tài)演化[J].山地學(xué)報(bào)，2014，32（1）：89-97.（CHANG M，TANG C，JIANG Z L，et al.Dynamic evolution process of sediment supply for debris flow occurrence in Longchi of Dujiangyan，Wenchuan Earthquake Area[J].Journal of Mountain Science，2014，32（1）：89-97.（in Chinese）） DOI：10.16089/j.cnki.1008-2786.2014.01.012.

[11] 喬建平，黃棟，楊宗佶，等.汶川地震極震區(qū)泥石流物源動(dòng)儲(chǔ)量統(tǒng)計(jì)方法討論[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào)，2012，23（2）：1-6.（QIAO J P，HUANG D，YANG Z Y，et al.Statistical method on dynamic reserve of debris flow′s source materials in meizoseismal area of Wenchuan earthquake region[J].Chinese Journal of Geological Hazard and Control，2012，23（2）：1-6.（in Chinese）） DOI：10.16031/j.cnki.issn.1003-8035.2012.02.011.

[12] 劉洋，唐川，李為樂(lè)，等.四川省都江堰市龍池地區(qū)群發(fā)性泥石流物源敏感性分析[J].災(zāi)害學(xué)，2013，28（2）：107-113.（LIU Y，TANG C，LI W L，et al.Sensitivity analysis of massive debris flow sources in Longchi Area，Dujiangyan City，Sichuan Province[J].Journal of Catastrophology，2013，28（2）：107-113.（in Chinese）） DOI：10.3969/j.issn.1000-811X.2013.02.022.

[13] 馬煜，李彩俠.汶川強(qiáng)震區(qū)龍溪河流域泥石流災(zāi)害分布規(guī)律研究——以四川省都江堰龍溪河8四川省泥石流為例[J].科學(xué)技術(shù)與工程，2015，15（10）：27-31.（MA Y，LI C X.Research on distribution of earthquake induced debris flow in Longxi River Area：a case study in the Longxi River of Dujiangyan，Sichuan，China[J].Science and technology and engineering，2015，15（10）：27-31.（in Chinese）） DOI：10.3969/j.issn.1671-1815.2015.10.005.

[14] 程霄.汶川映秀地區(qū)震后泥石流敏感性研究[D].成都：成都理工大學(xué)，2016.（CHENG X.Susceptibility assessment of debris flow after the earthquake in Yingxiu Area，Wenchuan County[D].Chengdu：Chengdu University of Technology，2016.（in Chinese））

[15] TANG C X，VAN WESTEN C J，TANYAS H，et al.Analysing post-earthquake landslide activity using multi-temporal landslide inventories near the epicentral area of the 2008 Wenchuan earthquake[J].Natural Hazards & Earth System Sciences，2016，16（12）：1-26.DOI：10.5194/nhess-16-2641-2016.

[16] 崔鵬，韋方強(qiáng)，何思明，等.5·12汶川地震誘發(fā)的山地災(zāi)害及減災(zāi)措施[J].山地學(xué)報(bào)，2008，26（3）：280-282.（CUI P，WEI F Q，HE S M，et al.Mountain disasters induced by the earthquake of May 12 in Wenchuan and the disasters mitigation[J].Journal of Mountain Science，2008，26（3）：280-282.（in Chinese）） DOI：10.16089/j.cnki.1008-2786.2008.03.007.

[17] 程霄，朱靜，王帥永，等.汶川縣漁子溪流域泥石流源區(qū)的崩滑體空間分布特征[J].長(zhǎng)江科學(xué)院院，2016，33（8）：52-58.（CHENG X，ZHU J，WANG S Y，et al.Spatial distribution characteristics of landslides and collapse masses in debris flow formation region along the Yuzixi River Basin of Wenchuan County，Sichuan Province[J].Journal of Yangtze River Scientific Research Institute，2016，33（8）：52-58.（in Chinese）） DOI：10.11988/ckyyb.20150497.

[18] 黃潤(rùn)秋，李為樂(lè).汶川大地震觸發(fā)地質(zhì)災(zāi)害的斷層效應(yīng)分析[J].工程地質(zhì)學(xué)報(bào)，2009，17（1）：19-28.（HUANG R Q，LI W L.Fault effect analysis of geo-hazard trjggered by Wenchuan Earthquake[J].Journal of Engineering Geology，2009，17（1）：19-28.（in Chinese）） DOI：10.3969/j.issn.1004-9665.2009.01.003.

[19] O′BRIEN J S.FLO-2D reference manual version 2009[R].2009.

[20] O′BRIEN J S，JULIEN P Y，F(xiàn)ULLERTON W T.Two-dimensional water flood and mudflow simulation[J].Journal of Hydraulic Engineering，1993，119（2）：244-261.DOI：10.1061/（ASCE）0733-9429（1993）119：2（244）.

[21] 賈濤，唐川，王納納.基于FLO-2D與沖量模型的泥石流危險(xiǎn)度分區(qū)方法及應(yīng)用[J].水電能源科學(xué)，2015，33（2）：152-155.（JIA T，TANG C，WANG N N.Method and application of debris flow hazard zoning based on FLO-2D and impulse model[J].International Journal Hydroelectric Energy，2015，33（2）：152-155.（in Chinese））

[22] 黃勛，唐川.基于數(shù)值模擬的泥石流災(zāi)害定量風(fēng)險(xiǎn)評(píng)價(jià)[J].地球科學(xué)進(jìn)展，2016，31（10）：1047-1055.（HUANG X，TANG C.Quantitative risk assessment of catastrophic debris flows through numerical simulation[J].Advances in Earth Science，2016，31（10）：1047-1055.（in Chinese）） DOI：10.11867/j.issn.1001-8166.2016.10.1047.

[23] 梁鴻熙，尚敏，徐鑫.基于FLO-2D數(shù)值模擬的泥石流流動(dòng)與堆積影響因素研究[J].工程地質(zhì)學(xué)報(bào)，2016，24（2）：228-234.（LIANG H X，SHANG M，XU X.Research on the influence factors of flow and deposition of debris flow based on the FLO-2D simulation[J].Journal of Engineering Geology，2016，24（2）：228-234.（in Chinese） DOI：10.13544/j.cnki.jeg.2016.02.008.

[24] 龔柯，楊濤，夏晨皓，等.基于FLO-2D的泥石流危險(xiǎn)性評(píng)價(jià)-以四川省汶川縣綿虒鎮(zhèn)簇頭溝為例[J].水資源與水工程學(xué)報(bào)，2017，28（6）：134-138.（GONG K，YANG T，XIA C H，et al.Assessment on the hazard of debris flow based on FLO-2D：A case study of debris flow in Cutou Gully，Wenchuan，Sichuan[J].Journal of Water Resources and Water Engineering，2017，28（6）：134-138.（in Chinese） DOI：10.11705/j.issn.1672-643X.2017.06.23.

[25] 常鳴.基于遙感及數(shù)值模擬的強(qiáng)震區(qū)泥石流定量風(fēng)險(xiǎn)評(píng)價(jià)研究[D].成都：成都理工大學(xué)，2014.（CHANG M.Quantitative risk assessment of debris flow in coseismic area based on remote sensing and numerical simulation[D].Chengdu：Chengdu University of Technology，2014.（in Chinese））