齊 信

(1.中國地質(zhì)調(diào)查局武漢地質(zhì)調(diào)查中心，湖北武漢 430205;2.成都理工大學(xué)地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護國家重點實驗室，四川成都 610059)

0 引言

泥石流是發(fā)育在斜坡上或溝谷中，飽含大量泥、砂、石的固、液相顆粒流體，是山區(qū)最常見的地質(zhì)災(zāi)害之一［1］。泥石流具有運動快速、暴發(fā)突然、能量巨大等特征，攜帶大量固體物質(zhì)沖出到溝口平坦地區(qū)堆積，常常淤埋鐵路、公路、橋梁、農(nóng)田等，給人類生命財產(chǎn)造成嚴(yán)重的危害。僅2010年全國發(fā)生泥石流災(zāi)害高達上百起，包括2010年8月7日甘肅舟曲特大泥石流［2-3］，2010年8月13～14日映秀鎮(zhèn)群發(fā)性泥石流［4-5］，2010 年 8 月 13 日龍池鎮(zhèn)群發(fā)性泥石流［6］，2010年8月13日綿竹清平群發(fā)性泥石流［7］等。鐵路工程屬于重大生命線工程，一旦遭受泥石流災(zāi)害的襲擊和破壞，后果極其嚴(yán)重，不幸的是類似災(zāi)害性事件時常發(fā)生，現(xiàn)在營運的鐵路沿線已發(fā)現(xiàn)和確認(rèn)的泥石流溝多達2000條，泥石流淤埋、沖毀車站和鐵路線路的嚴(yán)重災(zāi)害事件已超過300起。成都-蘭州鐵路(簡稱成蘭鐵路)位于四川省和甘肅省境內(nèi)，線路全長435.188km，運營總長700.071km，成蘭線于2011年2月26日開工，預(yù)計2017年竣工［8］。成蘭鐵路線從隧道延伸經(jīng)過水溝墩溝堆積扇，新塘關(guān)雙線大橋正好建設(shè)在扇形地的建設(shè)帶部位，過溝里程為DK223+258。水溝墩溝建設(shè)帶表面拱度大，溝槽左右擺動，淤積度大，對新塘關(guān)雙線大橋橋墩的危害大。扇形地即定帶位于扇形地的前端，寬度大，該帶越過G213國道直沖岷江對岸。野外調(diào)查表明水溝墩溝歷史上暴發(fā)過大規(guī)模泥石流，有典型泥石流堆積扇體。該溝經(jīng)過一段休整和物源累積后，具有再暴發(fā)泥石流的可能性，是一條隱蔽性極強的潛在泥石流溝。尤其是成蘭鐵路修建中及其運行后期，將長期受其泥石流災(zāi)害的威脅。因此，對水溝墩溝泥石流成因、動力學(xué)運動特征及對鐵路工程影響評價研究，有利于深入掌握水溝墩溝泥石流的形成條件、運動規(guī)律、動力學(xué)特征、危害方式、危險范圍，供鐵路部門制定相應(yīng)的泥石流防止對策，對成蘭鐵路線的長遠(yuǎn)安全運行和服務(wù)提供科學(xué)保障。

1 水溝墩溝流域特征

水溝墩溝位于松潘縣安宏鄉(xiāng)以南，岷江東側(cè)，水溝墩溝流域面積1.98km2，平均坡度約25°，最高海拔高度為3505m，溝口海拔高度為2685m，相對高差920m，總流域長1985m，主溝長1894m，平均坡降為387‰。水溝墩溝是典型的溝谷型泥石流溝，水溝墩溝泥石流溝可以劃分為3個區(qū)域，分別為形成區(qū)、流通區(qū)和堆積區(qū)，分區(qū)界限明顯清晰(圖1)，流域發(fā)育較為單一，溝道總體較為順直，整體走向以NE向為主。各分區(qū)要素見表1。

表1 水溝墩溝泥石流流域要素Table 1 Elements of the valley of the Shuigoudun debris flow

圖1 水溝墩溝流域圖Fig.1 Drainage map of the Shuigoudun gully

2 水溝墩溝泥石流成因分析

2.1 豐富的松散固體物質(zhì)

豐富的松散固體物質(zhì)是泥石流暴發(fā)的物質(zhì)基礎(chǔ)，其儲量大小決定了泥石流暴發(fā)規(guī)模和破壞程度大小。從地層巖性上看，水溝墩溝地層為三疊系新都橋組，巖性多是粉砂巖、千枚巖、變質(zhì)砂巖，這些巖石易于風(fēng)化破碎，極易形成崩坡積物等;從空間分布上看，水溝墩溝泥石流物源體主要為形成區(qū)和流通區(qū)段溝谷兩側(cè)分布的滑坡、崩塌體和溝谷內(nèi)泥石流松散堆積物。本次物源統(tǒng)計采用室內(nèi)遙感解譯和野外調(diào)查相結(jié)合方法，水溝墩溝物源共計大約10處(圖2)，結(jié)合溝道內(nèi)泥石流堆積體，總物質(zhì)體積約119.6×104m3。其中欠穩(wěn)定物源體積21.3×104m3，占物源總量的17.8%;局部欠穩(wěn)定物源體積46.8×104m3，占物源總量的39.2%;可參與近期泥石流活動量約36.9×104m3，占物源總量的30.9%。

圖2 水溝墩溝流域內(nèi)松散物源分布圖Fig.2 Distribution map of loose solid material in the Shuigoudun gully

2.2 降雨

降雨是泥石流暴發(fā)的誘發(fā)因子和動力條件。一方面，降雨改變邊坡土體的力學(xué)性能，導(dǎo)致其內(nèi)聚力下降，基質(zhì)吸力減小，抗剪強度降低，再加上水對滑動體的靜、動水壓力，都可以誘發(fā)滑坡、崩塌形成，為泥石流發(fā)生提供物質(zhì)基礎(chǔ);另一方面，降雨在溝道匯流后，使溝道受到了強烈的地表徑流或洪水沖刷侵蝕，溝道堆積物受“揭底作用”，演變?yōu)椤皾L雪球模式”，最終以“消防水管效應(yīng)”機制暴發(fā)泥石流。尤其是有前期降雨條件下，如遇短歷時強降雨，此溝暴發(fā)泥石流的可能性非常大。水溝墩溝流域地處四川盆地西北邊緣，屬于亞熱帶濕潤季風(fēng)區(qū)，降雨量豐富，多年平均降雨量838.2mm，歷年24h最大降雨量125.9mm(1984年8月3日)，2h最大降雨量76mm(1990年8月7日)，6月～9月為雨季，降雨占全年總雨量的70%～90%。特別是短歷時的強降雨對泥石流的激發(fā)起著決定性的作用，強降雨匯流入坡面沖溝，攜帶坡面松散物質(zhì)形成坡面泥石流，坡面泥石流匯流過程中沖刷、剝蝕沖溝，繼而擴大沖溝規(guī)模，如此滾動發(fā)展，沿途不斷有物質(zhì)加入，匯入主溝后行成較大規(guī)模泥石流。

2.3 陡峻地形

陡峻的地形條件有利于泥石流的暴發(fā)。首先，陡峻的地形給崩滑體物質(zhì)提供強大勢能，物質(zhì)在運動過程中將勢能轉(zhuǎn)化為動能，對溝道和坡面產(chǎn)生強烈沖刷和切割，且具有較強的搬運能力和沖擊動力，泥石流破壞強度增大;其次，陡峻的地形導(dǎo)致溝谷內(nèi)崩滑體及坡面松散堆積體更加不穩(wěn)定，有利于泥石流物質(zhì)的累積，為泥石流的暴發(fā)提供充足的物源。水溝墩溝相對高差920m，平均坡降達到387‰。如此大的高差和縱坡降加速了雨水的匯流，攜帶松散固體物質(zhì)，形成泥石流。

3 水溝墩溝泥石流靜動力學(xué)特征

3.1 泥石流容重

泥石流的容重獲取方法有直接測量法、現(xiàn)場攪拌法、經(jīng)驗公式法等。由于泥石流的暴發(fā)突發(fā)性、沖擊力大等特點，難以直接測得天然泥石流容重［9］。采用現(xiàn)場攪拌法進行泥石流容重的測定，往往因人目睹不同而異，同時因選取不同代表性堆積物而結(jié)果差距很大。在無法取得代表性堆積物樣品時，可根據(jù)泥石流溝易發(fā)程度數(shù)量化標(biāo)準(zhǔn)得到泥石流的容重，通過現(xiàn)場試驗和綜合判斷，確定水溝墩溝泥石流的容重為1.69 t/m3。

3.2 泥石流流速

泥石流流速是泥石流動力學(xué)特征的重要參數(shù)之一，也是計算其它動力學(xué)相關(guān)參數(shù)的重要因子。目前泥石流流速的計算公式多為半經(jīng)驗公式或經(jīng)驗公式法，一般分為兩類:①粘性泥石流流速計算;②稀性泥石流流速計算公式。根據(jù)泥石流容重計算結(jié)果和現(xiàn)場堆積扇沉積物沉積形態(tài)綜合判斷水溝墩溝泥石流為稀性泥石流，故采用稀性泥石流流速通用計算公式［10］，計算結(jié)果表明水溝墩溝泥石流具有較大的流速，這與水溝墩溝具有較大的溝道縱坡降有關(guān)。

式中:

VC——泥石流斷面平均流速(m/s);

φ——泥石流泥砂修正系數(shù);

γH——泥石流固體物質(zhì)重度(t/m3);

γc——泥石流重度(t/m3);

R——水力半徑(m);

I——泥石流水力坡度(‰)，用溝床平均縱坡降代替，本溝為0.387;

n——清水河床糙率系數(shù)。根據(jù)公式(1)和(2)計算，水溝墩溝泥石流的平均流速為7.17 m/s。

3.3 泥石流流量

常用的泥石流流量計算方法有雨洪法和泥痕法，泥痕法較為準(zhǔn)確，然而泥痕法要求野外能采集到新近發(fā)生泥石流的新鮮斷面，這樣計算的數(shù)據(jù)才為準(zhǔn)確。在水溝墩溝泥石流流域內(nèi)不易找到發(fā)生泥石流的新鮮斷面，因此本文采用雨洪法計算不同頻率下的泥石流的流量。雨洪法是假設(shè)泥石流與暴雨同頻率、且同步發(fā)生，按水文方法計算出斷面不同頻率下的小流域暴雨洪峰流量［9］，然后選用堵塞系數(shù)，按公式(3)計算泥石流流量［7］。

式中:

QC——頻率為P的泥石流洪峰值流量(m3/s);

Qp——頻率為P的暴雨洪水設(shè)計流量(m3/s);

φ——泥石流泥沙修正系數(shù);

Dc——泥石流堵塞系數(shù)，本文取2.00。

表2 泥石流流量計算參數(shù)及結(jié)果Table 2 Calculation results for debris flow discharge parameters

3.4 一次泥石流總量和最大沖出量

一次泥石流總量和最大沖出量是泥石流設(shè)防中最重要的參數(shù)，是攔和排等工程措施必須考慮的指標(biāo)。一般有實測法和公式法確定，本次按照公式法計算，一次泥石流總量和最大沖出量按下式計算［12］:

式中:

QH——一次泥石流固體物質(zhì)最大沖出量(m3);

Q——次泥石流總量(m3);

γc——泥石流重度(t/m3);

γW——水的重度(t/m3);

γH——泥石流固體物質(zhì)的重度(t/m3)。

其中，K值的變化隨流域面積的大小而變化，流域面積F=1.97km2，則K=0.202。本次泥石流持續(xù)時間取60min。依據(jù)公式(4)和(5)，計算結(jié)果見表3，計算結(jié)果和現(xiàn)場調(diào)查具有較好的一致性。

表3 一次泥石流總量和最大沖出量計算結(jié)果Table 3 Calculation results for debris flow transport material volume in total and contained solid material volume

3.5 泥石流最大淤積厚度

泥石流具有大沖大淤特點，尤其是鐵路工程，要考慮一次泥石流淤積厚度對工程的影響。通過不同頻率條件下泥石流的容重、屈服應(yīng)力及堆積區(qū)坡度等參數(shù)計算出水溝墩溝泥石流在不同頻率條件下可能淤積厚度［13］。本次采用泥石流體泥沙體積濃度與屈服應(yīng)力相結(jié)合而推倒的泥石流的最大淤積厚度公式進行計算［14］(略有改動)，得出水溝墩溝一次泥石流最大淤積厚度為3.70m。實際的泥石流淤積厚度可以小于最大淤積厚度，因此用最大淤積厚度作為泥石流的淤積厚度，結(jié)果會更為安全和可靠。

式中:

Hmax——一次泥石流最大淤積厚度;

γc——泥石流容重(t/m3);

θ——泥石流淤積坡度(°);

K、r——與泥石流的粘粒特性和粘粒含量有關(guān)的參數(shù);

g——重力加速度(m/s2);

c——泥沙體積濃度。

4 水溝墩溝泥石流對鐵路工程的影響評價

4.1 水溝墩溝泥石流危險性評價

4.1.1 泥石流危險度評價

泥石流的危險程度簡稱危險度，它是指遭到泥石流損害的可能性的大小［15］，危險度評價是對泥石流的發(fā)展趨勢所進行的預(yù)測評價，是對泥石流溝自身危險性大小的定量劃分。泥石流危險因子的賦值與其權(quán)重乘積和即為泥石流的危險度，采用劉希林等單溝泥石流危險度的計算方法［15］:

式中:

M——泥石流暴發(fā)規(guī)模;

F——泥石流暴發(fā)頻率;

S1——泥石流流域面積;

S2——主溝長度;

S3——流域相對高差;

S6——流域切割密度;

S9——不穩(wěn)定溝床比例。

通過將水溝墩溝的流域參數(shù)帶入公式(7)，計算得到泥石流危險度為0.46，根據(jù)危險度的判斷標(biāo)準(zhǔn)，水溝墩溝泥石流危險度為中度危險，可能造成重大災(zāi)害和危害。

4.1.2 泥石流最大危險范圍預(yù)測

泥石流危險范圍預(yù)測是對泥石流沖出最大距離進行估算方法的研究，是對泥石流災(zāi)害危險性評價的重要環(huán)節(jié)，用來計算泥石流最大危險范圍。汶川地震后，強震區(qū)豐富的松散固體物質(zhì)導(dǎo)致泥石流規(guī)模的放大效應(yīng)，即使泥石流流域面積較小，其洪峰流量比通常的要大數(shù)倍［16］。本次采用強震區(qū)一次泥石流沖出范圍危險預(yù)測模型［17］。經(jīng)過計算，水溝墩溝一次泥石流最大沖出距離0.494km，一次泥石流最大堆積寬度為0.491km。

式中:

Lf——預(yù)測的一次泥石流最大沖出距離(km);

Bf——預(yù)測的一次泥石流最大堆積寬度(km);

W——松散固體物質(zhì)儲量(×104m3);

H——流域最大相對高差(km)。

圈定泥石流最大危險范圍方法如下:①判斷泥石流堆積扇起點位置，以該點為起點，按照泥石流可能流動方向，標(biāo)一直線;②堆積扇起點為頂點，以該直線為頂點的平分線，繪出泥石流兩側(cè)邊界;③在頂點平分線上取Lf為泥石流最大堆積長度，垂直于平分線(Lf)取Bf為泥石流最大堆積寬度;④以Lf線段末端點和Bf線段兩個端點為基本點作圓弧，這樣構(gòu)成的圓弧與泥石流兩側(cè)邊界圍成范圍即為泥石流最大危險范圍(圖3)。經(jīng)過野外實地驗證，泥石流危險范圍略大于泥石流堆積扇范圍，表明預(yù)測結(jié)果在研究者可控范圍之內(nèi)，預(yù)測結(jié)果是可靠的。

4.2 水溝墩溝泥石流對鐵路工程影響評價

泥石流的危害方式可以簡化為兩種類型:一是沖，二是淤，泥石流的危害往往是既沖又淤，沖淤結(jié)合。水溝墩溝的縱坡降非常大，且溝道狹窄，泥石流

圖3 水溝墩溝泥石流最大危險范圍預(yù)測圖Fig.3 Maximal risk range of the Shuigoudun debris flow

4.3 水溝墩溝泥石流減災(zāi)防災(zāi)建議

水溝墩溝泥石流對成蘭鐵路及新塘關(guān)雙線大橋影響是潛在的，為使鐵路線正常、安全運行，必須對水溝墩溝泥石流采用相關(guān)減災(zāi)防災(zāi)措施。水溝墩溝泥石流危險度大小為0.46，屬于中度危險，水溝墩溝泥石流防治應(yīng)以治為主防為輔，實施生物和土建工程綜合治理即可防止泥石流的發(fā)生發(fā)展。

根據(jù)對水溝墩溝主溝堆積扇顆粒組成野外實地測試表明，溝床中碎石、塊石，最大塊徑為0.5m×0.6m×1.0m。通過對扇體溝口積塊碎石(粒徑＞10cm)共進行了1000塊顆粒的統(tǒng)計分析工作，而按塊數(shù)統(tǒng)計表明:10～20cm占80.47%，可見堆積區(qū)顆粒以10～20cm的碎石為主，沒有太大的石塊沖到堆積扇。鑒于以上分析，水溝墩溝泥石流雖然大塊石含量相對較少，但其沖擊破壞力仍在一定程度上對鐵路的建設(shè)構(gòu)成威脅，同時泥石流對堆積區(qū)的鐵路工程建設(shè)會造成潛在的泥沙淤積或侵蝕沖刷下切作用。因此，在設(shè)計和施工中，需要采取相應(yīng)的手段以增大橋墩的抗沖擊能力，同時還要在橋墩抗沖刷和下切能力加大的基礎(chǔ)上考慮泥石流可能淤埋的高度，以增大橋涵的凈空高度。其建議防治措施概括為三點:

(1)考慮提高線路標(biāo)高，以增大大橋的排洪凈空，同時盡可能使線路與主溝方向正交，不得己斜角角度不要大于10°～20°。

(2)需要采取相應(yīng)的手段以增大橋墩的抗沖擊力能力，尤其是橋墩迎水面要加強結(jié)構(gòu)的整體性和建材的耐磨性，同時還要在橋墩處設(shè)防以提高溝道的抗沖刷能力，防止泥石流因下切、揭底拉槽沖刷導(dǎo)致橋墩處巖土被掏空而引起橋體的垮塌。

(3)布置排導(dǎo)槽和攔擋壩工程，用以控制組成泥石流固體物質(zhì)的和雨洪徑流，削弱泥石流流量、下泄總量和能量，減少泥石流對新塘關(guān)大橋的沖刷、撞擊、淤積等危害。

5 結(jié)論

(1)水溝墩溝泥石流動靜力學(xué)參數(shù)表明該溝是一條典型的暴雨類、溝谷型、稀性、中型、低頻率泥石流溝。豐富的松散固體物質(zhì)、陡峻的地形和豐富的降水控制了水溝墩溝具有再次暴發(fā)泥石流的可能性。同時，水溝墩溝泥石流較高的平均流速和較大淤積厚度，表明水溝墩溝泥石流的直進性強，泥石流流體沖於量大，具有較強的沖擊力和侵蝕力，大沖大淤，沖淤結(jié)合。

(2)水溝墩溝泥石流危險度為0.46，屬于中度危險;一次泥石流最大堆積長度為494m，最大堆積寬度為491m，最大危險范圍為0.20 km2，成蘭鐵路新塘關(guān)雙線大橋完全處在水溝墩溝泥石流危害的范圍之內(nèi)，泥石流對鐵路工程及橋梁的危害模式為:先沖-又淤-再沖。

(3)成蘭鐵路及新塘關(guān)雙線大橋在設(shè)計和施工中，必須采取相應(yīng)的手段以增大橋墩的抗沖擊能力和下切能力，同時考慮泥石流可能淤埋的高度，以增大橋涵的凈空高度，減少泥石流對新塘關(guān)大橋橋墩的沖刷、撞擊、淤積等危害。

［1］吳積善，田連權(quán)，康志成，等.泥石流及其綜合治理［M］.北京:科學(xué)出版社，1993:40-48.WU Jishan，TIAN Liangui，KANG Zhicheng，et al.Debris flowand general prevention［M］.Beijing:Science Press，1993:40-48.

［2］Tang C，Rengers N.，van Asch Th.W.J.，et al.Triggering conditions and depositional characteristics of a disastrous debris flow event in Zhouqu city，Gansu Province，northwestern China［J］.Nat.Hazards Earth Syst.Sci，2011，11:2903-2912.

［3］胡向德，畢遠(yuǎn)宏，魏潔.甘肅舟曲三眼峪溝泥石流的形成條件與發(fā)展趨勢［J］.中國地質(zhì)災(zāi)害與防治學(xué)報，2011，22(2):55-60.HU Xiangde，BI Yuanhong，WEI Jie.Conditions leading to the destructive debris flow at the Sanyanyu gulley and its occurrence trend［J］.The Chinese Journalof Geological Hazard and Control，2011，22(2):55-60.

［4］許強.四川省8.13特大泥石流災(zāi)害特點、成因與啟示［J］.工程地質(zhì)學(xué)報，2010，18(5):596-608.XU Qiang.The 13 August 2010 catastrophic debris flows in Sichuan province:characteristics，genetic mechanism and suggestions ［J］.Journal of Engineering Geology，2010，18(5):596-608.

［5］唐川，李為樂，丁軍，等.汶川映秀鎮(zhèn)“8.14”特大泥石流災(zāi)害調(diào)查［J］.地球科學(xué)(中國地質(zhì)大學(xué)學(xué)報)，2011，36(1):172-180.TANG Chuan，LIWeile， DING Jun， etal.Field investigation and research on giant debris flow on August 14，2010 in Yingxiu town ，epicenter of Wenchuan earthquake［J］.Earth Science(JournalofChina University of Geosciences)，2011，36(1):172-180.

［6］余斌，馬煜，張建楠，等.汶川地震后四川省都江堰市龍池鎮(zhèn)群發(fā)泥石流災(zāi)害［J］.山地學(xué)報，2011，29(6):738-746.YU Bin，MA Yu，ZHANG Jiannan，et al.The group debris flow hazards after the Wenchuan earthquake in Longchi，Dujiangyan，Sichuan province［J］.Journal of Mountain Science，2011，29(6):738-746.

［7］倪化勇，鄭萬模，唐業(yè)旗，等.綿竹清平8·13群發(fā)泥石流成因、特征與發(fā)展趨勢［J］.水文地質(zhì)工程地質(zhì)，2011，38(3):129-138.NI Huayong，ZHENG Wanmo，TANG Yeqi，et al.Formation，characteristics and trend of the group debris flows occurred on August 13 in Qingping，Mianzhu county［J］.Hydrogeology＆ Engineering Geology，2011，38(3):129-138.

［8］http://baike.baidu.com/view/1410689.htm.

［9］鐵永波，唐川，蘇小琴.四川省九龍縣石頭溝泥石流泥石流動力學(xué)特征及危險性評價研究［J］.水土保持研究，2008，15(5):168-170.TIE Yongbo，TANG Chuan，SU Xiaoqin.Shitougou debris flow dynamic characters and its hazard risk assessment in Jiulong of Sichuan provience［J］.Researchof Soil and Water Conservation，2008，15(5):168-170.

［10］費祥俊，舒安平.泥石流運動機制與災(zāi)害防治［M］.北京:清華大學(xué)出版社，2003:80-95.FEI Xiangjun，SHU Anping.Movement mechanism and disaster control for debris flow［M］.Beijing:Tsinghua University Press，2003:80-95.

［11］四川省水利廳.四川省中小流域暴雨洪水計算手冊［M］.成都，1984:10-52.Si Chuan provincial water resources department.Sichuan province smalland medium-sized basin storm flood calculation manual［M］.Chengdou，1984:10-52.

［12］中華人民共和國國土資源部.泥石流災(zāi)害防治工程勘查規(guī)范［S］.2006.The ministry of land and resources of the People’s Republic of China.Specification of Geological Investigation for Debris Flow Stabilization［S］.2006.

［13］中國科學(xué)院成都山地災(zāi)害與環(huán)境研究所，中國科學(xué)院東川泥石流觀測研究站.云南蔣家溝泥石流觀測研究［M］.北京:科學(xué)出版社，1990:53-196.IMDE(Institution of Mountain Disaster and Environment，CAS)，Research and observation of debris flow in Jiangjia Gully，Yunnan，China［M］.Beijing:Science Press，1990:53-196.

［14］余斌.不同容重的泥石流淤積厚度計算方法研究［J］.防災(zāi)減災(zāi)工程學(xué)報，2010，30(2):207-211.YU Bin.Study on the method for deposition depth calculation of debris flow with different densities［J］.Journal of Disaster Prevention and Mitigation Engineering，2010，30(2):207-211.

［15］劉希林，唐川.泥石流危險性評價［M］.北京:科學(xué)出版社，1995:17-30.LIU Xilin，TANG Chuan.Hazard assessment of debris flow［M］.Beijing:Science Press，1995:17-30.

［16］Tang C，Zhu J，Li w L.Rainfall triggered debris flows after Wenchuan earthquake.Bul1.Eng.Geo1.Environ.，2009，68:187-194.

［17］朱靜，常鳴，丁軍，等.汶川震區(qū)暴雨泥石流危險范圍預(yù)測研究［J］.工程地質(zhì)學(xué)報，2012，20(1):7-14.ZHU Jing，CHANG Ming，DING Jun，et al.Prediction for hazard zones of rainstorm induced debris flows in the WenChuan earthquake epicenters［J］.Journal of Engineering Geology，2012，20(1):7-14.