趙永輝,萬志杰,涂國祥,鄧 輝
1.西藏自治區水利電力規劃勘測設計研究院,西藏 拉薩 850000;2.成都理工大學 地質災害防治與地質環境保護國家重點實驗室,四川 成都 610059
近些年,受全球氣候變化影響,中國西南地區堰塞湖事件頻發,給當地經濟造成重大損失,對其破壞力談虎色變。事實上,堰塞湖事件在全世界屢見不鮮,如1911年,受地震影響,塔吉克斯坦帕米爾高原穆爾加布河谷右岸產生山體滑坡,形成壩高約500多米的薩列茲堰塞湖,威脅下游塔吉克斯坦、阿富汗、烏茲別克斯坦等國家[1-2];1959年8月17日,美國蒙大拿州強震引發美國地震史上有記載的最大山體滑坡,致使麥迪遜河堵塞,并在不到一個月的時間內形成了水深約60 m、長約10 km的堰塞湖[3];1999年10月6日,新西蘭亞當斯山體滑坡,形成120 m高的堰塞體,阻塞波伊魯阿河,堰塞湖庫容約500×104~700×104m3[4];2008年5月12日,中國四川汶川地震觸發山體滑坡,滑坡堆積體阻塞河道,形成壩高約105 m,庫容約3.2×108m3的唐家山堰塞湖[5];2014年8月3日,中國云南昭通魯甸地震引發紅石巖水電站上游河段山體滑坡,造成曲靖市會澤縣牛欄江堵塞,形成庫容約2.6×108m3的堰塞湖,堰塞體高約96 m、方量約1 200×104m3[6]。特別是2018年,中國西藏金沙江堰塞湖和雅魯藏布江堰塞湖事件后,堰塞湖瞬間成為當時地質災害的代名詞,給中國西藏高原防治減災工作增加了難度和壓力。
截至目前,眾多專家學者對堰塞湖相關專題做了大量研究,推動了堰塞湖科研工作的不斷發展。鑒于西藏高原地區氣候、水文、地質等多方面的特殊性,針對中國西藏高原堰塞湖演化模式、易發區劃等方面的研究相對較少,而其又是高海拔地區防災減災研究工作的重點。文章選擇6處西藏典型堰塞湖事件作為主要分析和論述對象,以期通過收集西藏典型堰塞湖相關資料,明晰西藏堰塞湖演化模式和易發性區劃,旨在為今后西藏防災減災和堰塞湖相關研究提供支撐與參考。
西藏典型堰塞湖主要位于西藏東部和南部,其中昌都市分布2處,即玉曲堰塞湖和金沙江堰塞湖;林芝市分布4處,分別為帕隆藏布堰塞湖、易貢藏布堰塞湖、雅魯藏布江堰塞湖和嘎貢溝堰塞湖,見圖1。

①玉曲堰塞湖;②金沙江堰塞湖;③帕隆藏布堰塞湖;④易貢藏布堰塞湖;⑤雅魯藏布江堰塞湖;⑥嘎貢溝堰塞湖
玉曲位于中國西藏昌都市左貢縣東南部,阿左隆巴溝為玉曲一級支溝,東經98°28′4.29″、北緯28°51′54.88″。玉曲堰塞湖為阿左隆巴溝溝谷型泥石流引發的次生地質災害。
玉曲河谷坡度約為40°~50°,屬于高山峽谷區。阿左隆巴溝呈EW向,縱長約為7.6 km,橫寬約為3.3 km,匯水面積約為1.3×107m2。溝谷后緣及兩側基巖山體主要為三疊系上統察瓦龍群變質砂巖、板巖和礫巖(T3C),巖層產狀C為N14°E/60°SW,后緣溝壁坡向及坡度P為N10°E/40°~50°SW,溝道底部堆積結構松散的冰磧物和崩積碎石土(Qcol+gl),溝口主要為松散泥石流堆積扇。經遙感影像解譯,溝口泥石流堆積扇分兩期形成,其中第一期堆積扇(Qdf-1)分布面積較廣,結構稍密,表面植被長勢較好,第二期堆積扇(Q1df-2)面積相對較小。
2014年8月之前阿左隆巴溝就有泥石流爆發歷史,溝口分布大量的第一期泥石流堆積扇(Qdf-1)。在此基礎上,2014年8月,以溝道底部冰磧物和崩積碎石土,以及第一期泥石流堆積扇為物源,冰川融水和降雨作為泥石流主要誘發因素,阿左隆巴溝爆發二次泥石流,泥石流堆積體堵塞玉曲河道,形成寬約265 m、高約25 m的堰塞體,上游湖區水位急劇上漲,庫容達500×104m3[7]。由于堰塞壩為泥石流堆積體,其主要由結構松散、架空、抗沖刷性能差的碎礫石土組成,加之,河流凹岸侵蝕、凸岸沉積作用,堰塞體與玉曲左岸連接部位產生潰決,形成自然排導槽。后期,在堰塞湖湖水下蝕和側蝕作用下,下泄排導槽不斷擴大,直至堰塞體侵蝕殆盡,見圖2。
2014年8月,泥石流已攜帶大量溝道松散物源沖出阿左隆巴溝。據遙感解譯,溝內仍殘留大量松散物源,為再次爆發泥石流提供了物源基礎。鑒于此,玉曲阿左隆巴溝段仍存在形成堰塞湖的可能。

a-堰塞湖DEM;b-堰塞湖形成前遙感影像(據Googel Earth);c-堰塞湖形成后遙感影像(據Googel Earth)
金沙江堰塞湖位于中國西藏昌都市江達縣東南部,東經98°42′0.00″、北緯31°4′48.00″,其為金沙江右岸山體滑坡引發的次生地質災害。
金沙江堰塞湖地處高山峽谷區,區域型斷裂較發育,以逆斷層為主,區域穩定性較差。金沙江右岸滑源區斜坡坡度和坡向P約為S15°W/32°SE,出露地層主要為二疊系-三疊系崗托巖組絹云石英片巖(PT1g),產狀C為N148°W/67°SW,淺表層呈強風化,深部主要為強風化-弱風化巖體。斜坡山體主要發育兩組結構面:即J1為S155°E/56°NE,J2約為S30°W/37°SE[8-9]。
2018年10月11日,受強降雨影響,滑源區巖土體沿SW105°整體失穩,產生山體滑坡,堵塞金沙江河谷,形成縱長約220 m、橫寬約130 m、方量約2.2×107m3的滑坡堰塞體及庫容約為2.9×108m3的金沙江堰塞湖。隨著堰塞湖湖水位不斷上漲,于10月12日產生漫壩,并開始自然泄流,截至10月13日,堰塞湖完成泄流。由于10月11日滑坡為后緣滑坡邊界以外斜坡巖土體提供了更為有利的臨空環境,2018年11月3日,滑坡后緣不穩定性巖土體產生二次滑移,隨之而來的是金沙江二次堵江,堰塞湖庫容約為5.24×108m3。與第一次相比,金沙江第二次滑坡堵江堰塞湖規模更大,更具威脅。經人工排導,第二次堰塞湖于11月12日開始泄洪,13日堰塞湖水位與下游水位基本持平,至此,金沙江堰塞湖險情解除[10-14],見圖3。
特別是,由于金沙江右岸坡度較陡,且滑坡源區距離金沙江水位線約為600~1 300 m,斜坡整體失穩后極速沖入金沙江,形成極具能量的滑坡涌浪,致使滑坡對岸(左岸)坡表植被沖蝕殆盡,該區域為滑坡涌浪影響區。因此,與金沙江堰塞體相連接的左右兩岸岸坡處變形破壞跡象截然不同,其中左岸為涌浪影響區,表層殘留物為干燥泥漿,右岸為剪切-滑移區,分布多組順坡向滑痕。

a-滑坡源三維影像;b-堰塞湖遙感影像(據李俊峰,有修改);c-堰塞湖現場照片(據李俊峰,有修改)
帕隆藏布堰塞湖位于西藏林芝市波密縣東南部古鄉溝溝口處,古鄉溝為帕隆藏布一級支溝,東經95°26′30.80″、北緯29°54′47.46″。帕隆藏布堰塞湖為古鄉溝溝谷型泥石流引發的次生地質災害。
古鄉溝縱長約6.0 km,橫寬約為5.3 km,溝口寬度約為0.3 km,匯水面積約為為2.6×107m2,總體呈SN向發育,與帕隆藏布近垂直。古鄉溝山體主要為黑云母花崗閃長巖(Yj1),溝道底部堆積結構松散的冰磧物和崩積碎石土(Qcol+gl)。
古鄉溝地處西藏低緯高山峽谷區,降雨充沛,雪線約為5 000 m,以上為古鄉海洋性冰川,古鄉溝為冰川泥石流溝。在冰川融水、降雨等影響下,古鄉溝內地表徑流流量不斷增大,水力條件也隨之不斷改善。當水力條件達到一定峰值時,溝道內出現侵蝕破壞,主要包括下切侵蝕、溯源侵蝕和側向侵蝕。其中,下切侵蝕和溯源侵蝕作用導致物源區溝道松散堆積體迅速形成泥石流,側向侵蝕致使溝谷岸坡部位松散堆積體失穩破壞,使泥石流物源不斷增加。最終,泥石流以摧古拉朽之勢沖出古鄉溝,堵塞帕隆藏布,形成帕隆藏布堰塞湖。迄今為止,古鄉溝已暴發幾十余次泥石流,為典型高頻泥石流溝。其中1953年9月29日泥石流規模最大,屬于特大型冰川泥石流,沖出松散碎屑顆粒約為1.1×107m3,形成長約5 km、寬約1 km的帕隆藏布堰塞湖[15],見圖4。
易貢藏布堰塞湖位于西藏林芝市波密縣西部扎本龍溝溝口處,扎本龍溝為易貢藏布一級支溝,東經94°56′15.13″、北緯30°10′41.28″。易貢藏布堰塞湖為扎本龍溝“崩滑流”地質災害引發的次生地質災害。
扎本龍溝縱長約7.3 km,橫寬約為3.0 km,匯水面積約為為1.6×107m2,總體呈SW向發育,與易貢藏布呈斜交關系。扎本龍溝后緣坡度和坡向P為N40°W/40°~50°SW,出露以花崗巖、大理巖和板巖為主的混雜巖AnZgd,發育兩組優勢裂隙:J1為S32°E/32°NE,J2為S42°E/59°NE[16]。
2000年4月9日前,氣溫回暖,扎木弄溝后緣海拔5 000m附近冰雪不斷融化,加之,期間持續性降雨,大量水流灌入后緣山體“X”型裂隙中,致使崩滑源區山體飽和。2000年4月9日,扎木弄溝后緣山體產生崩滑地質災害。由于崩滑源距離溝底垂直高差較大,其整體沿“X”型裂隙失穩的同時,以強大的沖擊力撞擊溝底,并迅速解體為高速碎屑流,以銳不可當之勢刮鏟溝道底部松散碎屑物質,最終碎屑流攜帶大量松散體沖出扎木弄溝,堵塞易貢藏布,形成長約2.5 km、寬約2.5 km、方量約2.8×108~3.0×108m3的堰塞壩,其上游即為易貢堰塞湖。2000年之前,易貢藏布就有堵江史,易貢湖即為碎屑流堵江所致[17-18],見圖5。

a-堰塞湖DEM;b-堰塞湖三維模型(據Googel Earth);c-易貢藏布扎本龍溝“崩滑流”照片[19]
雅魯藏布江堰塞湖位于西藏林芝市波密縣西部色東普溝溝口,色東普溝為雅魯藏布江一級支溝,東經94°56′24.23″、北緯29°44′54.07″。雅魯藏布江堰塞湖為色東普溝“崩滑流”地質災害引發的次生地質災害。
色東普溝縱長約1.1 km,橫寬約為0.5 km~0.8 km,匯水面積約為7.0×107m2,總體近SN向發育,與雅魯藏布江斜交。色東普溝后緣山體坡度和坡向P約為N70°E/40°~50°SE,地層巖性為派鄉巖組片巖和片麻巖(AnP),巖層產狀C為N60°E/75°~85°NW,后緣溝道堆積大量松散冰磧物和崩積碎石土(Qcol+gl)[20]。
色東普溝冰川屬于海洋型冰川,在近些年全球變暖的大背景下,其消融速度加快,冰崩災害頻繁。同時,色東普溝處于高寒地區,溝谷兩側山體在凍融循環破壞作用下導致巖體強度較低,碎裂松動巖體分布較廣,致使溝內時常發生巖崩。加之,2018年,研究區降雨量偏多,溝道松散體容重不斷增大,抗剪強度隨之不斷降低。2018年10月11日,色東普溝后緣溝道松散體產生剪切-滑移變形,即色東普溝滑坡。由于色東普溝匯水面積較大,且溝道底部主要為松散碎屑堆積體,滑坡“啟動”后,松散架空的滑體瞬間崩解為碎屑流,并以“破竹之勢”下切和側向侵蝕溝道巖土體,最終于色東普溝口形成堰塞壩,完全堵塞雅魯藏布江[21],見圖6。

a-“崩滑”區照片;b-碎屑流區照片;c-碎屑流堆積體(堰塞體)照片
嘎貢溝堰塞湖位于西藏林芝市朗縣北部雅魯藏布江一級支溝嘎貢溝內,東經93°6′34.33″、北緯29°11′1.86″。嘎貢溝堰塞湖為溝谷左岸山體滑坡引發的次生地質災害。
嘎貢溝地處岡底斯巖漿弧構造單元,位于雅魯藏布江結合帶內,走滑為研究區主要構造特征,屬于區域穩定性較差區。嘎貢溝左岸山體滑源區坡度P為40°~50°,主要為花崗片麻巖(Pt2-3b),發育一組優勢裂隙J1,產狀為N30°W/5°~30°SW,少量發育裂隙J2(N60°W/60°~80°NE)和J3(N50°W/75°SW)[22]。
嘎貢溝左岸為凹岸,在河流凹岸侵蝕作用下,左岸斜坡花崗片麻巖(Pt2-3b)順坡內J1結構面整體向嘎貢溝溝谷剪切-滑移,堆積堵塞河道,形成縱長約1 600 m、橫寬約600~900 m、高約130 m、估算方量約2.2×108m3的滑坡塊碎石土堰塞壩。隨著湖水不斷匯集,上游形成面積約1.2×105m2的嘎貢堰塞湖。后期堰塞湖湖水不斷上漲直到漫壩,在河水凹岸侵蝕作用下,堰塞體于左岸連接處產生潰決,形成自然下泄通道,經歷漫長的地質歷史后,進而演化為現今的嘎貢溝瀑布[22]。此外,滑坡區凍融風化較強,后緣滑壁處碎裂松動巖體較發育,現今堰塞體與溝谷岸坡連接處均殘留大量崩坡積體。

a-堰塞體潰決口上端;b-堰塞體(滑坡堆積體);c-堰塞體潰決口下端
西藏典型堰塞湖實為次生地質災害,其原生地質災害往往是滑坡、崩塌、泥石流等原生地質災害。筆者將西藏堰塞湖演化過程歸納為三類,即泥石流堵江型、滑坡堵江型和“崩滑流”堵江型。
泥石流堵江型堰塞湖的演化過程伴隨有明顯的階段性。具體為:(1)松散堆積體形成階段。在西藏高寒高海拔山區凍融、風化等強烈作用下,溝谷后緣及溝道底部滑坡堆積體、崩塌堆積體和冰水堆積體較發育,該區域作為泥石流爆發的形成區,提供充足的松散物源,見圖8(a);(2)泥石流形成階段。在強降雨和冰川融水作用下,松散的滑坡堆積體、崩塌堆積體和冰水堆積體物源瞬間解體為流體,形成溝谷型泥石流,見圖8(b);(3)在溝道流水作用下,泥石流迅速經過流通區,以迅雷不及掩耳之勢沖出溝口,堆積于溝口河道,堵塞河道,形成堰塞湖,泥石流堆積體構成堰塞體Qdf,見圖8(c)。前述西藏玉曲堰塞湖和帕隆藏布堰塞湖就屬該類堰塞湖。

a-松散堆積體;b-泥石流形成階段;c-堰塞體形成階段
滑坡堵江型堰塞湖是指河谷斜坡產生滑坡,并堵塞河道形成堰塞湖的過程,由滑坡堆積體構成堰塞壩。具體為:西藏地質構造活動強烈,河谷下切侵蝕作用強烈,但是河谷下切初始階段屬于寬谷時期,河谷斜坡穩定性較好,見圖9(a)。隨著河谷下切侵蝕不斷加劇,岸坡坡度不斷變陡,臨空條件不斷改善,形成西藏高山峽谷地貌。當河谷下切至一定深度時,斜坡巖土體在降雨、地震等其它誘發因素影響下整體失穩,形成滑坡地質災害,見圖9(b)?;w以摧枯拉朽之勢剪切破壞斜坡巖土體,最終堆積于河道,形成堰塞湖,堰塞體物質組成主要為滑坡堆積體Qdel,見圖9(c)。這類高原堰塞湖常常依托于大型或特大型滑坡,一般小規模滑坡很難產生堵江,如西藏嘎貢溝堰塞湖和金沙江堰塞湖。

a-原始斜坡;b-滑坡階段;c-堰塞體形成階段
“崩滑流”堵江型堰塞湖的形成需要崩塌、滑坡和碎屑流的共同“參與”。在全球氣候變暖影響下,西藏高海拔溝谷后緣時常發育崩塌(巖崩、冰崩)或滑坡地質災害。在大量冰川融水、大氣降水等水力條件下,崩滑體迅速解體,形成高速遠程碎屑流,沿途刮鏟溝道松散堆積體,最終沖出溝口堵塞江河,形成堰塞湖,見圖10。該類型堰塞湖壩體物質組成主要為碎屑流堆積體Qcol+del,西藏易貢堰塞湖和雅魯藏布江堰塞湖就屬該類型。

a-崩塌(冰崩、巖崩)、滑坡;b-碎屑流階段;c-堰塞體形成階段
綜上所述,三種堰塞湖演化模式既有相似之處,又有其不同之點,具體見表1。
此外,堰塞壩形成后,上游堰塞湖淹沒區周邊岸坡將孕育崩塌、滑坡等地質災害,如雅魯藏布江堰塞湖淹沒道路,湖區岸坡發育滑坡,見圖11(a)、圖11(b)。堰塞壩潰決或人工泄洪后,下游洪水區岸坡同樣易誘發崩塌、滑坡等地質災害,如金沙江堰塞湖洪水下泄后,致使下游金沙江索橋和公路被沖毀,見圖11(c)、圖11(d)。由此可見,堰塞湖本身為崩塌、滑坡、泥石流等原生地質災害的次生災害,其又可延伸出下一級次生地質災害。至此,形成了以“‘崩滑流’—堵江(堰塞湖)—淹沒—洪水—‘崩滑’”為主線的地質災害鏈[23]。

表1 西藏典型堰塞湖演化模式分類

a-雅魯藏布江堰塞湖淹沒公路;b-雅魯藏布江堰塞湖區滑坡;c-金沙江堰塞湖洪水沖毀下游索橋[24];d-金沙江堰塞湖下游洪水區滑坡[24]
西藏典型堰塞湖世紀形成時間不算久遠,除嘎貢溝堰塞湖不詳外,其它典型堰塞湖最近一次堵江事件均產生于21世紀,主要分布于西藏東部和東南部。其主要原因在于西藏西部主要為隔壁荒漠地貌,地形較緩,干旱少雨,河流分布較少。相比之下,藏東屬于高原溫帶半干旱季風氣候,中國藏南地區屬于亞熱帶、熱帶季風氣候,大江大河密布,高山峽谷地貌發育。隨著全球變暖帶來的溫室效應,西藏東南高山峽谷區冰川大面積消失,雪線快速后退,冰川融水和大氣降水不斷增多。加之,在凍融、風化等作用下,西藏高海拔河谷區斜坡山體架空松散,強度較低。同時,冰川退卻過程中,高海拔區凍土加速消融,形成大量松散堆積物,為原生地質災害(崩塌、滑坡、泥石流)的發育提供了充足的松散物源。當藏東南河谷斜坡及溝谷產生原生地質災害后,極易產生堵江事件,形成堰塞湖。不難發現,西藏典型堰塞湖形成必須滿足兩個條件:一是大江大河區;二是原生地質災害高易發區。即當西藏某一地區同時滿足上述兩個條件時,該區域可初步判定為堰塞湖高易發區。因此,“兩區”,即:大江大河區,原生地質災害高易發區,這是形成西藏典型堰塞湖事件的充分條件,并以此作為西藏典型堰塞湖易發性的初步判據。

a-西藏西部阿里地區普蘭縣寬谷地貌;b-西藏中部日喀則薩迦縣寬谷地貌;c-西藏中部日喀則市康馬縣寬谷地貌;d-西藏南部聶拉木縣喜馬拉雅山南坡典型深切河谷地貌;e-西藏東部林芝市波密縣帕隆藏布深切河谷地貌
據不完全統計,西藏大江大河主要分布于藏東、藏中、藏南和藏西,其中藏東、藏南發育金沙江、瀾滄江、怒江、雅魯藏布江,及其主要支流,該地區山高坡陡,多為高山峽谷地貌,原生地質災害(崩塌、滑坡、泥石流)較發育;藏西和藏中主要發育森格藏布(獅泉河)、朗欽藏布(象泉河)、當卻藏布(馬泉河)等,但該地區多為戈壁荒漠,河谷寬廣,原生地質災害相對不發育。在此基礎上,基于西藏典型堰塞湖易發性的初步判據,西藏堰塞湖易發區主要集中于藏東和藏南,其中藏東以巴青縣—工布江達縣—朗縣—隆子縣為界線,以西為低易發區,以東為高易發區;藏南高易發區主要位于喜馬拉雅山南坡吉隆縣—聶拉木縣—亞東縣一帶。此外,雅魯藏布江日喀則仁布縣—曲水縣段為峽谷地帶,且該段具有發育原生地質災害的潛力,因此,該段亦為西藏堰塞湖易發區,見圖12。
(1)近些年,受全球溫室效應影響,西藏冰川不斷消融,雪線快速后退,降雨量偏高,地質災害頻發頻發。西藏典型堰塞湖實為次生地質災害,其原生地質災害往往是滑坡、崩塌、泥石流等地質災害,以此將西藏堰塞湖演化過程歸納為三類,即泥石流堵江型、滑坡堵江型和“崩滑流”堵江型。堰塞體形成后,上游堰塞湖湖區岸坡易產生崩塌、滑坡等地質災害,堰塞壩潰決或人工泄洪后,下游洪水區岸坡同樣易誘發崩塌、滑坡等地質災害。由此,形成一條以“‘崩滑流’—堵江(堰塞湖)—淹沒—洪水—‘崩滑’”為主線的地質災害鏈。
(2)西藏典型堰塞湖分布于西藏東部和東南部大江大河原生地質災害易發區,由此提出西藏堰塞湖易發區初步判據——大江大河區、原生地質災害高易發區。即當滿足充分條件“兩區”時,可初步判定為堰塞湖高易發區,未來西藏堰塞湖預測應以高山河谷區的崩塌、滑坡、泥石流等原生地質災害為主。依據該判據,初步建立西藏堰塞湖易發性區劃,具體以巴青縣—工布江達縣—朗縣—隆子縣—亞東縣—聶拉木縣—吉隆縣為界線,即藏東、藏南及雅魯藏布江日喀則仁布縣—曲水縣段為堰塞湖高易發區,藏西和藏中為堰塞湖低易發區,以上地區也是未來西藏堰塞湖預防的重點區域。
研究區地處西藏高原高山峽谷區,山高溝深,交通不便。特別的,喜馬拉雅山脈北部阿里地區、薩迦縣、康馬縣局部地區海拔5 000m以上,高寒缺氧,部分現場一手資料實屬來之不易,十分具有研究價值。鑒于中國西藏堰塞湖數量之多、規模之大、成因之復雜,而筆者學術水平有限,西藏堰塞湖事件研究中尚有諸多不足之處,觀點的繆誤和資料的遺漏在所難免,敬請各位專家學者批評指正。
【致謝:感謝西藏自治區水利電力規劃勘測設計研究院萬志杰高級工程師,成都理工大學地質災害防治與地質環境保護國家重點實驗室涂國祥教授、鄧輝教授在西藏地質災害研究過程中給予的幫助。特別感謝西藏自治區水利廳、西藏自治區水利電力規劃勘測設計研究院各級領導,在西藏雅魯藏布江堰塞湖、金沙江堰塞湖等地質災害發生后第一時間派駐專業技術人員進駐現場,以及對后期西藏自治區地質災害研究工作的支持。】