張金山,謝 洪,王小丹,范繼輝
(1. 中國(guó)科學(xué)院 山地災(zāi)害與地表過(guò)程重點(diǎn)實(shí)驗(yàn)室,四川 成都 610041;2. 中國(guó)科學(xué)院·水利部 成都山地災(zāi)害與環(huán)境研究所,四川 成都 610041)
?
西藏尖姆普曲泥石流*
張金山1,2,謝 洪1,2,王小丹2,范繼輝2
(1. 中國(guó)科學(xué)院 山地災(zāi)害與地表過(guò)程重點(diǎn)實(shí)驗(yàn)室,四川 成都 610041;2. 中國(guó)科學(xué)院·水利部 成都山地災(zāi)害與環(huán)境研究所,四川 成都 610041)
尖姆普曲為雅魯藏布江支流帕隆藏布江左岸的一級(jí)支溝,是一條大型冰川泥石流溝,2007年曾暴發(fā)特大規(guī)模的泥石流,堵斷帕隆藏布江,導(dǎo)致對(duì)岸的川藏公路斷道,并造成人員傷亡。目前幾乎每年都會(huì)暴發(fā)中小規(guī)模的泥石流。尖姆普曲高差巨大、岸坡陡峻,給泥石流暴發(fā)提供了優(yōu)越的能量條件,流域所在區(qū)域隆升強(qiáng)烈,地質(zhì)構(gòu)造復(fù)雜,巖石破碎,加上地震頻發(fā),風(fēng)化強(qiáng)烈,給泥石流發(fā)育提供了大量的固體物質(zhì)條件;流域所在區(qū)位于西南季風(fēng)影響區(qū),又為水汽通道,降水豐富且多暴雨,冰川活動(dòng)也十分強(qiáng)烈,為泥石流暴發(fā)提供了大量的水源條件。根據(jù)調(diào)查訪問(wèn),尖姆普曲內(nèi)泥石流的主要激發(fā)水源為暴雨,若疊加增溫和雨水潛熱導(dǎo)致冰雪融水,則可產(chǎn)生大規(guī)模泥石流。根據(jù)實(shí)測(cè)資料進(jìn)行了泥石流流速、流量等參數(shù)的計(jì)算,結(jié)果表明暴雨和冰雪融水組合條件較好時(shí),激發(fā)的大規(guī)模泥石流可再次堵斷帕隆藏布江。
泥石流;尖姆普曲;西藏;帕隆藏布江;松繞村
尖姆普曲位于西藏波密縣松繞村境內(nèi),為雅魯藏布江支流帕隆藏布江左岸的一級(jí)支溝,溝口與東偏南方向的波密縣城直線距離約45 km,距西偏南方向的林芝地區(qū)政府所在的八一鎮(zhèn)約100 km(圖1)。
圖1 尖姆普曲地理位置示意圖
從2006年的遙感影像看,尖姆普曲溝口已存在較大的堆積扇(圖2),并存在泥石流堆積臺(tái)地(圖3),因此認(rèn)為該溝是一條老的泥石流溝。
2007年9月4日,尖姆普曲發(fā)生泥石流,泥石流越過(guò)帕隆藏布江沖毀溝對(duì)岸的川藏公路約150 m(圖4),導(dǎo)致溝口村莊及過(guò)往路人1人死亡、7人失蹤、9人受傷,2戶(hù)民房被毀、6戶(hù)民房受損,沖毀農(nóng)田2 hm2,沖走牲畜40余頭,溝口下游兩岸大片森林被毀,泥石流短暫堵斷帕隆藏布江(約1 h),1座跨江吊橋(松饒吊橋)和1輛汽車(chē)被毀,川藏公路交通中斷43 h。據(jù)調(diào)查訪問(wèn),該溝近年來(lái)小規(guī)模泥石流幾乎每年都會(huì)暴發(fā),2009年還暴發(fā)過(guò)較大規(guī)模的泥石流,但未堵斷主河。
圖2 遙感影像顯示的尖姆普曲溝口堆積(googleearth影像 2006年4月30日)
圖3 尖母普曲溝口
圖4 尖母普曲2007年9月4日泥石流沖過(guò)帕隆藏布江危害川藏公路及村莊
1.1 地貌特征
尖姆普曲發(fā)源于念青唐古拉山東段余脈,由南向北注入帕隆藏布江,為帕隆藏布江中游左側(cè)支流,流域面積為17.65 km2,主溝長(zhǎng)度7.04 km,溝床平均縱比降294‰,流域最高峰海拔5 590 m,溝口帕隆藏布江水面海拔2 439 m,相對(duì)高差3 151 m,流域內(nèi)山坡平均坡度達(dá)35°,地形十分陡峻(圖5)。在流域海拔4 000 m以上的山坡和溝谷,發(fā)育有現(xiàn)代冰川,總面積約1.45 km2,屬海洋性溫冰川,補(bǔ)給與消融均十分強(qiáng)烈(圖6)。
圖5 尖姆普曲流域圖
圖6 尖姆普曲上游的冰川
1.2 地質(zhì)條件
尖姆普曲流域在地質(zhì)構(gòu)造上處于印度板塊和歐亞板塊的縫合帶的北側(cè),發(fā)育于察隅—波密—林芝近東西向深大斷裂帶的中段。流域所在區(qū)域經(jīng)前寒武紀(jì)及古生代構(gòu)造運(yùn)動(dòng),產(chǎn)生一套中度變質(zhì)巖系,燕山構(gòu)造運(yùn)動(dòng)中形成大量斷裂。
流域出露地層以前震旦系為主,巖性主要為片麻巖、片巖、變粒巖及大理巖。由于強(qiáng)烈的變質(zhì)作用和寒凍風(fēng)化作用,以及構(gòu)造的破壞和地震的影響,巖體破碎,崩塌滑坡發(fā)育,是泥石流活動(dòng)的主要固體物質(zhì)來(lái)源。喜山運(yùn)動(dòng)以來(lái),本區(qū)及鄰近地區(qū)發(fā)生強(qiáng)烈的隆升和向東擠出走滑[1-2],區(qū)內(nèi)有NE、EW、NNW、NWW四個(gè)方向的構(gòu)造體系,它們彼此以斜接或反接形式復(fù)合于雅魯藏布江大拐彎頂端和南迦巴瓦峰附近[3-4],造成本區(qū)斷裂切割密度很大,巖層被切割的支離破碎,加之流域上游晝夜溫差大,寒凍風(fēng)化作用強(qiáng)烈,在強(qiáng)烈的雪崩作用下補(bǔ)給冰川,在消融區(qū)冰川表面形成深厚的表磧塊礫層,在全球變暖的總體趨勢(shì)下,流域冰川退縮,表磧、側(cè)磧、終磧等冰川堆積物成為泥石流的另一固體物質(zhì)來(lái)源。
尖姆普曲所處位置是西藏兩個(gè)強(qiáng)地震帶之一的藏東南強(qiáng)地震帶中部,對(duì)區(qū)域有影響的特大地震有:1897年印度阿落姆8.7級(jí)特大地震及1950年8月我國(guó)西藏察隅8.5級(jí)特大地震; 1976年以后共發(fā)生影響該區(qū)(92°~98°E, 28°~32°N)的6級(jí)以上的地震1次,5次以上的地震4次,4級(jí)以上的地震78次以上[5]。頻頻發(fā)生的中、小地震破壞了山體的穩(wěn)定性,促使崩塌滑坡形成,加大松散固體物質(zhì)的儲(chǔ)備量,增強(qiáng)松散物質(zhì)參與泥石流活動(dòng)的可能性,在加快泥石流活躍周期,以及加大泥石流的規(guī)模等方面都起到了積極的作用。
1.3 氣象水文條件
尖姆普曲地處西藏高原的東南邊緣,深受印度洋西南季風(fēng)的影響,水汽溯雅魯藏布江下游河谷逆流而上,翻越帕隆藏布江南岸山嶺的埡口后直入本區(qū),帶來(lái)充沛的降水。由于海洋性季風(fēng)的影響,雨季長(zhǎng)達(dá)6個(gè)月左右(4月下旬到9月中旬);冬季則受西風(fēng)環(huán)流控制,以固態(tài)降水為主。根據(jù)波密氣象站和古鄉(xiāng)溝零星觀測(cè)資料,年平均氣溫8.6℃,降雨884.5 mm。年內(nèi)降雨分配極不均勻,雨季降水占年降水量的70%以上(圖7),最大月份與最小月份降雨懸殊達(dá)20余倍。另外,流域降水量也受海拔高度的影響,降水量隨海拔明顯增大,據(jù)與該溝溝口位置相距不到20 km的古鄉(xiāng)溝流域的觀測(cè)資料,流域中夏季降水垂直遞減率達(dá)80 mm/100 m[6]。
圖7 尖姆普曲所在區(qū)域各月降水分配圖
尖姆普曲上游冰川面積約1.45 km2,冰川融水為流域徑流重要來(lái)源之一,但因缺乏觀測(cè)資料,故使用“冰川的年平均消融量大致等于雪線處的年消融量”的原理,計(jì)算得出年平均冰川消融徑流[7]。
先用式(1)計(jì)算雪線處的消融深A(yù)(mm):
A=1.33(9.66+ts)2.85。
(1)
式中:ts為計(jì)算年夏季流域雪線處6-8月的平均氣溫(℃),利用波密氣象站觀測(cè)資料及氣溫垂直遞減原理,算得流域雪線處夏季平均氣溫為7.5 ℃,代入式(1)得尖姆普曲2007年冰川消融深為4 370 mm。
尖姆普曲流域的冰川屬海洋性溫冰川,冰的密度宜取0.9 g/cm3,計(jì)算消融徑流深h(mm):
h=4 370×0.9=2 933.1。
(2)
最后根據(jù)冰川面積計(jì)算2007年冰川的年平均消融徑流Q(m3/s):
Q=(2 933.1×1.45×1 000)/(365×24×3 600)=0.18。
(3)
據(jù)2014年7月考察時(shí)實(shí)測(cè)流量為6.2 m3/s,而據(jù)附近培龍溝、古鄉(xiāng)溝等流域研究資料,7月流量一般占年總徑流量的26%左右,即流域年平均徑流量約為1.98 m3/s,可見(jiàn)冰雪融水僅占年徑流量的9.1%。
尖姆普曲冰川面積占流域總面積的8.2%,利用面積分割法[6],根據(jù)同緯度同為海洋性冰川的海螺溝冰川末端多年觀測(cè)徑流分割比率,算得尖姆普曲冰雨融水占全年徑流量的比率是:
(0.483×0.082)/0.38=0.104,
(4)
即在尖姆普曲流域每年徑流問(wèn)題構(gòu)成中占10.4%。
由以上分析可見(jiàn),尖姆普曲流域平時(shí)的徑流主要以降水和地下水,冰雪融水僅占年平均徑流量的10%左右。因此,尖姆普曲泥石流的激發(fā)條件以暴雨為主,如果暴雨伴隨著流域溫度增加,加上雨水潛熱,則冰川消融量會(huì)大大增加,且與冰磧層等碎屑物相互作用,從而形成大規(guī)模泥石流。
2.1 泥石流的形成機(jī)理
調(diào)研中未了解到尖姆普曲在2007年之前暴發(fā)泥石流的記錄,但調(diào)查到流域下游段兩岸及溝口保存的高出溝床數(shù)十米混雜堆積臺(tái)地,且明顯分為多級(jí),所以分析認(rèn)為該堆積物臺(tái)地為歷史泥石流堆積形成,在區(qū)域地殼持續(xù)抬升的背景下,主河帕隆藏布江下切而形成。其暴發(fā)頻率分析認(rèn)定為百年一遇(P=1%)。
據(jù)分析,該類(lèi)泥石流的形成過(guò)程為:前期降雨使坡面和溝岸碎屑物質(zhì)浸潤(rùn)飽和,致使土體平衡條件遭受破壞。當(dāng)大量的水流在坡面或溝谷運(yùn)動(dòng)時(shí),一方面,沖刷一部分碎屑物質(zhì),使其與水流一起產(chǎn)生運(yùn)動(dòng),另一方面,一部分冰雪融水滲入到冰磧層中,使其飽和失穩(wěn)。溝床流動(dòng)的水體不斷掀揭沖刷不穩(wěn)定的溝床與邊坡,引發(fā)大量物質(zhì)坍塌,從而不斷擴(kuò)大自身的規(guī)模和含沙量,若堵塞溝道還會(huì)造成潰決洪水,形成泥石流[6]。
根據(jù)殘留的泥石流堆積物中細(xì)粒物質(zhì)少(圖8),塊石含量高,長(zhǎng)徑2 m者在溝內(nèi)多見(jiàn),堆積龍頭明顯,分選性差等特點(diǎn),認(rèn)為該溝暴發(fā)的泥石流為高容重泥石流,估計(jì)其容重可能達(dá)2.2~2.3 t/m3。
圖8 尖姆普曲內(nèi)泥石流堆積物顆粒級(jí)配圖
2.2 泥石流特征參數(shù)計(jì)算
ZHANG Ming-ming, ZHAO Pei, YU Yue-qing, ZHANG Cui-gai, GAO Wei, LI Fang
在溝谷下游,泥石流翻越了堆積臺(tái)地,對(duì)臺(tái)地上的植被造成破壞,經(jīng)分析為泥石流超高堆積所致。選擇溝道內(nèi)岸邊該次泥石流的堆積物質(zhì)的高度認(rèn)定為泥石流所達(dá)到的高度,測(cè)量了此處的橫斷面(圖9)及縱比降,可據(jù)此計(jì)算這次泥石流的最大流速和流量。
圖9 尖姆普曲下游溝道斷面圖
尖姆普曲泥石流雖然容重大,但因其塊石含量高,粘粒含量少,在運(yùn)動(dòng)特征上更接近西北地區(qū)的稀性泥石流,因此流速計(jì)算選用鐵道部第一勘測(cè)設(shè)計(jì)院推薦公式(5),該公式是西北地區(qū)的經(jīng)驗(yàn)公式[8]。
(5)
式中:Vc為泥石流斷面平均流速(m/s);Rc為泥石流流體水力半徑(m),可近似取泥石流泥位深度;I為泥石流流面縱坡(‰);a為阻力系數(shù),可查表獲得。
為了確定選用公式與參數(shù)的正確性及最后計(jì)算結(jié)果的可靠性,還選用通用式(6)計(jì)算流速,該公式是根據(jù)西藏古鄉(xiāng)溝、云南蔣家溝和甘肅火燒溝199次泥石流3 000多陣次觀測(cè)資料得出,資料來(lái)源于不同地區(qū)多種泥石流,其適應(yīng)范圍寬[8],可用于本研究區(qū)泥石流的參數(shù)計(jì)算。
(6)
式中:Vc為泥石流斷面平均流速(m/s);nc為粘性泥石流的溝床糙率,可用內(nèi)插法查表獲得;Hc為計(jì)算斷面的平均泥深(m),Ic為泥石流流面縱坡(‰),一般可用溝床縱坡代替。
泥石流的總量可用式(7)計(jì)算:
Wc=19TQc/72
(7)
用式(5)和式(6)計(jì)算的尖姆普曲泥石流流速與流量,及使用式(7)計(jì)算的泥石流總量如表所示(表1)。
表1 尖姆普曲泥石流特征值
2.3 泥石流堵河計(jì)算
據(jù)臨近尖姆普曲的原102道班地段附近帕隆藏布江加馬其美水文觀測(cè)斷面資料,1953-1996年帕隆藏布江43年平均流量是420.8 m3/s。其中,7月流量最大,44年里,7月平均流量是1 288.2 m3/s[6]。因觀測(cè)數(shù)據(jù)的欠缺,沒(méi)能得到2009年9月4日泥石流暴發(fā)時(shí)的主河流量,但根據(jù)一般情況推測(cè),此時(shí)主河流量應(yīng)較平時(shí)流量大為增加,再考慮9月份帕隆藏布江的流量處于回落期,因此取7月平均流量來(lái)參考泥石流暴發(fā)時(shí)的主河流量,事后實(shí)測(cè)尖姆普曲入?yún)R段帕隆藏布江流速約6 m/s。與計(jì)算得到的泥石流特征參數(shù)相比,泥石流流速、流量均較主河稍大,具有堵斷主河的可能性。為具體分析堵河可能性大小及堵河程度,對(duì)該次泥石流的堵河特征進(jìn)行計(jì)算。
通過(guò)對(duì)比分析,張金山得出的岷江上游泥石流堵河可能性的經(jīng)驗(yàn)公式具有參數(shù)獲取相對(duì)容易的優(yōu)點(diǎn),模型表達(dá)式體現(xiàn)出了主河與入?yún)R泥石流的對(duì)比關(guān)系,具有明確的物理意義,故被本研究所采用。該公式通過(guò)實(shí)例分析及專(zhuān)家問(wèn)卷調(diào)查,結(jié)合泥石流對(duì)主河的影響度的概念和表達(dá)式,并考慮數(shù)據(jù)取得的難易程度及可靠性而確立,其表達(dá)式為:
(8)
式中:R為泥石流堵河危險(xiǎn)度,數(shù)值越大表示堵斷主河的可能性越大;ρn、Qn和Jn分別為泥石流容重、流量和泥石流溝床比降;Kz、Qz和Jz分別為入?yún)R處主河寬度、流量及比降[9]。
經(jīng)計(jì)算,尖姆普曲2007年“9·4”泥石流的堵河危險(xiǎn)度為3.84(表2),也就是說(shuō),支溝泥石流的能量、物質(zhì)量在一定寬度的入?yún)R口處表現(xiàn)出相對(duì)于主河的絕對(duì)優(yōu)勢(shì),因此堵斷主河的可能性非常大,這與實(shí)際是相符的。
表2 泥石流堵河危險(xiǎn)性計(jì)算表
(1)尖姆普曲高差巨大、岸坡陡峻,給泥石流暴發(fā)提供了優(yōu)越的能量條件,流域所在區(qū)域隆升強(qiáng)烈,地質(zhì)構(gòu)造復(fù)雜,巖石破碎,加上地震頻發(fā),風(fēng)化強(qiáng)烈,給泥石流發(fā)育提供了大量的固體物質(zhì);流域所在區(qū)位于西南季風(fēng)影響區(qū),又為水汽通道,降水豐富且多暴雨,冰川活動(dòng)也十分強(qiáng)烈,為泥石流暴發(fā)提供了豐沛的水源。
(2)根據(jù)調(diào)查訪問(wèn),尖姆普曲內(nèi)泥石流的主要激發(fā)水源為暴雨,若疊加增溫和雨水潛熱導(dǎo)致冰雪融水,則可產(chǎn)生大規(guī)模泥石流。根據(jù)實(shí)測(cè)資料進(jìn)行了泥石流流速、流量等參數(shù)的計(jì)算,結(jié)果表明暴雨和冰雪融水組合條件較好時(shí),激發(fā)的大規(guī)模泥石流可堵斷帕隆藏布江。
(3)本文在調(diào)查、計(jì)算的基礎(chǔ)上對(duì)尖姆普曲2007年9月4日的一場(chǎng)泥石流進(jìn)行了初步的研究,泥石流特征參數(shù)的計(jì)算結(jié)果符合實(shí)際,堵河可能性的計(jì)算結(jié)果也反映了真實(shí)的情況,因此本文所采用的方法可以用于該類(lèi)泥石流的科學(xué)研究和工程實(shí)踐,但在使用過(guò)程中應(yīng)注意進(jìn)行檢驗(yàn)。
[1] 鐘大賚, 吳根耀, 丁林. 藏東新生代陸內(nèi)變形特征[M]. 廣州: 廣東科技出版社, 1998: 335-364.
[2] 潘裕生. 青藏高原的形成與隆升[J]. 地學(xué)前緣, 1999, 6(3): 153-163.
[3] 張進(jìn)江, 季建清, 鐘大賚, 等. 東喜馬拉雅南迦巴瓦構(gòu)造結(jié)的構(gòu)造格局及形成過(guò)程探討[J]. 中國(guó)科學(xué):D輯, 2003, 33(4): 373-383.
[4] 章振根, 劉玉海, 王天武, 等. 南迦巴瓦峰地區(qū)地質(zhì)[M]. 北京: 科學(xué)出版社, 1992.
[5] 呂儒仁, 唐邦興, 朱平一. 西藏泥石流與環(huán)境[M]. 成都: 成都科技大學(xué)出版社, 1999.
[6] 中國(guó)科學(xué)院·水利部成都山地災(zāi)害與環(huán)境研究所, 西藏自治區(qū)交通廳科學(xué)研究所. 川藏公路典型山地災(zāi)害研究[M]. 成都: 成都科技大學(xué)出版社, 1999.
[7] 中國(guó)科學(xué)院·水利部成都山地災(zāi)害與環(huán)境研究所, 中國(guó)科學(xué)院蘭州冰川凍土研究所, 西藏自治區(qū)交通廳科學(xué)研究所. 川藏公路南線(西藏境內(nèi))山地災(zāi)害及防治對(duì)策[M]. 北京: 科學(xué)出版社, 1995.
[8] 國(guó)家防汛抗旱指揮部辦公室, 中國(guó)科學(xué)院·水利部成都山地災(zāi)害與環(huán)境研究所. 山洪泥石流滑坡災(zāi)害及防治[M]. 北京: 科學(xué)出版社, 1994: 180-181.
[9] 張金山,謝洪. 岷江上游泥石流堵河可能性的經(jīng)驗(yàn)公式判別[J]. 長(zhǎng)江流域資源與環(huán)境, 2008, 17(4): 651-655.
《災(zāi)害學(xué)》被評(píng)為“RCCSE中國(guó)核心學(xué)術(shù)期刊(A)”
《災(zāi)害學(xué)》編輯部2015年5月收到中國(guó)科教評(píng)價(jià)網(wǎng)、武漢大學(xué)中國(guó)科學(xué)評(píng)價(jià)研究中心、武漢大學(xué)圖書(shū)館的證書(shū),在第四屆《中國(guó)學(xué)術(shù)期刊評(píng)價(jià)研究報(bào)告(武大版)》(2015-2016年)中,《災(zāi)害學(xué)》雜志被評(píng)為“RCCSE中國(guó)核心學(xué)術(shù)期刊(A)”。
《災(zāi)害學(xué)》編輯部
Debris-flowofJianmupuquRavineinTibet
Zhang Jinshan1, 2, Xie Hong1, 2, Wang Xiaodan1and Fan Jihui2
(KeyLaboratoryofMountainHazardsandEarthSurfaceProcesses,InstituteofMountainHazardsandEnvironmen,CAS,Chengdu610041,China)
JianmupuquRavine,abranchravineofParlungZangboRiver,isalargeglacierdebris-flowravine.Superdebris-flowshavebrokenoutin2007,andhaveblockedParlungZangboRiverandblockedoffSichuanTibethighwaywhichlocateattheothersideoftheParlungZangboRiver,themoreseriousishavetakencasualties.Now,debris-flowswithmediumand/orsmallscalesoccuralmosteveryyear.ThesteepslopeandthehugealtitudedifferenceoftheJianmupuquRavineprovidepredominantenergyconditionsforthedebris-flow.Thedrainagebasinislocatedinthefastupliftedarea,wherethecomplicatedgeologicstructure,thecrackedrock,andthefrequentearthquakemaketherocksexperiencestrongweathering,thusplentyofgranularmaterialsareavailablefortheformationofdebris-flows.Theregionislocatedinthesouthwestmonsooninfluencearea,andwatervaporchannel,abundantandconcentratedrainfall,aswellasthestrongglacieractivityprovideswatersourcefordebris-flow.Accordingtoinvestigationandanalysis,mostdebris-flowsintheravineareinducedbytherainstorm,whenthemeltedwaterisoverlaid,thelargescaledebris-flowsmayoccur.Parametriccalculationsuchastheflowvelocityandtherunoffisconductedaccordingtothemonitoringdata.Theresultshowsthatlargedebris-flowscanbearousedwhentherainstormandthemeltedwaterarecombinedwell,maybeblockingoffParlungZangboRiveragain.
debris-flow;JianmupuquRavine;Tibet;theParlungZangboriver;Songraovillage
2014-11-25
2015-01-13
水利部公益性行業(yè)科研專(zhuān)項(xiàng)經(jīng)費(fèi)項(xiàng)目(201301037);中國(guó)科學(xué)院STS計(jì)劃項(xiàng)目(KFJ-EKFJ-EW-STS-094)
張金山(1972-),男,甘肅古浪人,副研究員,理學(xué)博士,主要研究方向?yàn)槟嗍鳒p災(zāi)理論與方法. E-mail:zjszj@163.com
P642.23;X43
A
1000-811X(2015)03-0099-05
10.3969/j.issn.1000-811X.2015.03.019
張金山,謝洪,王小丹,等. 西藏尖姆普曲泥石流[J].災(zāi)害學(xué), 2015,30(3):099-103. [Zhang Jinshan, Xie Hong, Wang Xiaodan,et al. Debris-flow of Jianmupuqu Ravine in Tibet[J].Journal of Catastrophology, 2015,30(3):099-103.]