李元灵, 王军朝, 陈 龙, 刘建康, 杨东旭, 张佳佳

(1.中国地质科学院 探矿工艺研究所, 成都 611734; 2.中国地质调查局 地质灾害防治技术中心, 成都 611734)

帕隆藏布流域位于西藏东南部，是我国海洋性冰川的主要分布地区[1]。由于该区域地质构造复杂、新构造运动活跃、山高谷深、降水充沛及现代冰川活跃，导致流域内发生的冰川泥石流数量多、规模大、分布集中、暴发频率高、危害大[2]。近几十年来，帕隆藏布流域内频频暴发的泥石流灾害严重影响着当地居民的生活，并阻碍G318的正常通行。国内学者亦对该区域的泥石流分外关注，开展了大量的研究工作，分析该区域泥石流形成的气象背景、地貌条件、泥石流启动特征等。20世纪50年代为修建川藏公路，以齐树春为代表的工程技术人员和西藏自治区交通厅科学研究所的专家对帕隆藏布流域的冰川泥石流进行了探索研究，积累了大量基础资料，对冰川泥石流也形成了初步认识；60,70年代，以中科院兰州冰川冻土研究所的施雅风[3]、杜榕桓[4]等为主的研究人员进行了半定位观测研究，并对区内的冰川泥石流进行了广泛的考察研究，提出了著名的古乡冰川泥石流流速计算公式；70年代中期—90年代，中科院成都山地灾害与环境研究所和兰州冰川冻土研究所等对区内泥石流进行了较为全面的考察、勘测和研究，较为透彻地分析了冰川泥石流的形成机制及发展趋势，提出了川藏公路西藏境内沿线的泥石流流速计算公式[5]；2000年以来，成都山地灾害与环境研究所在区内持续开展观测及研究工作，得到了不少研究成果[6]。

区内泥石流一直较为活跃，尤其在目前全球变暖的气候条件下，区内冰川强烈消融，释放出大量松散冰碛物及冰雪融水，为泥石流活动提供了良好的条件。近年来，区内泥石流时有发生。2016年9月5日，帕隆藏布流域古乡段相邻几条泥石流沟集中暴发，造成较为严重的财产损失，并阻碍了G318国道的正常通行。这几条泥石流暴发时间及空间上均高度集中，分析此次群发性泥石流事件，利于掌握该区域泥石流灾害的现状，推测其发展趋势，指导该区泥石流的监测预警工作，为防灾减灾工作提供支撑，对川藏铁路及高等级公路的修建有重要参考价值。

1 区域概况

帕隆藏布位于西藏自治区东南部，属青藏高原强烈切割的边缘山地，山高谷深，地形坡度大。地跨冈底斯陆块、雅鲁藏布江结合带、喜马拉雅陆块东部南迦巴瓦构造结，为一强烈挤压、碰撞、旋扭走滑、急剧伸展隆升、地质构造极其复杂的造山带。近3 Ma以来的地壳垂直隆升速率达3 mm/a[5]。帕隆藏布流域范围内主要断裂为嘉黎断裂，沿帕隆藏布河谷展布，总体走向近东西向，呈向北凸出的弧形。该区域地震活动频繁，20世纪40年代至今发生3级以上地震近200余次[7]。地层上归属于藏滇地层大区，岩浆岩、变质岩相对发育，未变质的沉积岩零星出露。高山区基岩裸露，岩性以花岗岩、片麻岩为主，坡麓为崩坡积物，河谷为冲洪积物。帕隆藏布流域海洋性冰川发育，由于朝向不同，两岸沟谷地貌特征及冰川发育特征存在明显差异。流域北岸为阳坡，沟谷内冰川因接受更多日照而更易消融，沟谷上游多发育有冰斗和冰蚀围谷，围谷内堆积有大量冰碛物，沟谷下游沟道狭窄，岸坡陡立；南岸为阴坡，日照相较于北岸要少，因而消融量少，冰川运动更为强烈，因而南岸沟谷多为冰川U型谷。区内降雨充沛，主要受益于来自印度洋的暖湿气流。波密县城年均降雨多在800～1 000 mm。在通麦、东久一带，由于逆江而上的印度洋气流受高山阻截，形成多雨、暴雨区，年降水量可达1 127 mm。基于帕隆藏布流域上述地貌、地质、气象水文等特征，该区域泥石流形成条件好，暴发频次高、规模大[8]。

2 泥石流活动特征

2016年9月5日，位于帕隆藏布流域古乡段的比通沟、赤担隆巴、角隆沟同时暴发泥石流，其中比通沟泥石流致使电塔架设施工工棚损毁、G318短时中断，赤担隆巴泥石流致使当地居民发电设施被毁，道路被淤埋，造成一定经济损失。本次暴发泥石流的地理位置如图1所示。由于比通沟、赤担隆巴泥石流分别位于帕隆藏布南北两岸，在地形地貌、物源特征等方面均具有各自的特性，也代表着帕隆藏布该段南北两岸的泥石流的特征，且二者暴发规模较大、致灾较重，以下论述比通沟与赤担隆巴的泥石流活动特征及成因。

2.1 比通沟泥石流

2016年9月5日清晨6：30，比通沟泥石流暴发，部分泥石流冲出物跃上距沟底4.6 m的钢桥，淤埋公路，阻碍了G318的通行(图2)，并损毁了左岸的电力施工工棚，所幸未造成人员伤亡。在汇入主河前，泥石流强烈磨蚀新修筑的排导槽，在排导槽表面形成约5～10 cm深的凹槽，致使混凝土内钢筋裸露出来，对排导槽的使用寿命造成不良影响。

图1研究区位置

图2比通沟泥石流跃上钢桥淤埋公路

根据实地调查及试验，测算得泥石流的容重为γc=1.8 t/m3，为粘性泥石流。选定出山口处稳定沟道，根据粘性泥石流流速的通用计算公式(1)，计算该处泥石流流速Vc=8.77 m/s，按公式(2)计算峰值流量为447.27 m3/s，按公式(3)计算其本次过流总量为2.01万m3。

(1)

Qc=Wc·Vc

(2)

Q=K·T·Qc

(3)

式中：Vc为泥石流流速；nc为粘性泥石流的沟床糙率；Hc为平均泥深；Ic为泥位纵坡率；Qc为峰值流量；Wc为过流断面面积；Q为过流总量；T为泥石流历时；K取0.037 8。

2.2 赤担隆巴泥石流

同日中午11点，位于帕隆藏布南岸的赤担隆巴暴发泥石流，损毁多个小型发电设施，摧毁转经房，淤埋道路(图3)，致使当地居民无法出行，对其生活造成严重影响。据实地调查，沟道内泥痕3 m，沟道沿途多处宽缓处可见大量巨石、块石架空堆积，靠沟道内侧块石块径较大，往外侧逐渐减小。由于主沟道方向在出山口处以近直角的角度由向北转为向西，泥石流未沿原沟道冲出，而是翻过沟道，堆积于老扇之上，新堆积扇宽60余米。部分泥石流冲出物质直抵主河，运动过程中掏蚀河堤形成凹槽，由于主河水流速度较大，未在河岸出现大量堆积、挤压主河的现象。

图3赤担隆巴泥石流淤埋公路

根据上文所述方法，计算得出山口处泥石流流速Vc=6.56 m/s，峰值流量为Qc=267.64 m3/s，过流总量约2.2万m3。本次群发泥石流的基本信息如表1所示。

表1 泥石流基本特征

3 泥石流成因分析

3.1 地形地貌条件

地形地貌条件是泥石流形成的必要条件，控制着泥石流的启动、搬运与堆积[9]。帕隆藏布流域地貌上属深切割高山峡谷地貌，山势险峻、河流深切，为泥石流的形成提供了良好的地形地貌条件。

比通沟流域形态呈漏斗状(图4)，利于地表径流的汇流集中。流域面积23.7 km2，主沟长7.1 km。流域最高点高程为5 560 m，沟口与帕隆藏布交汇处高程为2 429 m，流域相对高差为3 131 m，巨大的流域高差为泥石流活动提供了良好的地势条件。沟道平均纵比降321.7‰，物源区和流通区平均纵坡降约为464.4‰。沟域两侧地势陡峭，平均坡度大于32°，沟道纵坡降较大，尤其是主沟上游段及各支沟纵坡降都在450‰以上，有利于降雨的汇集。比通沟物源区呈宽U型谷，利于物源的积累聚集。流通区沟道狭窄，呈V型谷，易形成堵溃，从而增加泥石流的峰值流量。根据实地调查，流通区下段约500 m长沟道内产生堵溃点5处，对比无人机航拍影像及遥感影像，推测沟道内堵溃点达数十处。这使泥石流在冲出的过程中，不断地产生“堵溃循环”，从而使得泥石流物质逐级积累，规模增大。这也是本次比通沟暴发泥石流规模不大，龙头却足以跃上沟口钢架桥的原因。

图4比通沟泥石流流域全貌

赤担隆巴流域形态近矩形(图5)，流域面积27.5 km2，主沟长8 km，主沟道上下游坡度较为一致，平均纵比降233.9‰，为泥石流发育提供良好的地势条件。沟源处发育有冰蚀围谷，利于物源、水源的聚集。沟谷上下游均呈宽U型谷，沟道整体较为宽缓，泥石流暴发后多顺沟下泄，不易产生堵溃。沟道两岸山坡坡度较为一致，约30°。流域后缘较为平缓，沟道纵剖面呈直线。因而从地貌特征来看，赤担隆巴流域发育程度较高，但仍具有形成泥石流的地形条件。

3.2 物源条件

该区域由于地质构造活动剧烈、冰川作用、寒冻风化及冻融风化作用强烈等原因，松散物质储量丰富[8]。该区域泥石流物源包括冰碛物源、崩滑物源、沟道物源、坡面物源几类。帕隆藏布流域海洋性冰川广布，使得该区域内冰碛物分布广、规模大，冰碛物厚度达数十米甚至数百米。而近代以来受全球气候变暖的影响，现代冰川剧烈消融、退缩，冰川中的固体物质不断地释放出来，补充泥石流松散物源。并且，现代冰川前缘冰碛物结构松散、稳定性差，尤在冰雪融水或大气降水条件下易失稳转化为泥石流物质。因而，冰碛物源是该区域泥石流的主要物源。

图5赤担隆巴泥石流流域全貌

通过实地调查及遥感影像调查对本次群发性泥石流的物源进行调查，查明比通沟与赤担隆巴泥石流物源类型及其总量和动储量。调查发现，比通沟泥石流的物源以沟源处冰斗内巨厚松散的冰碛物源为主(表2)。冰碛物磨圆差、无分选、未胶结，经沟道侵蚀形成高陡临空面(图6)，稳定性差，易在降雨或冰雪融水作用下失稳进入沟道，参与泥石流活动。沟道两岸可见大量崩滑物源，虽规模不大，但由于流通区沟道狭窄，崩滑物源进入沟道后造成多处堵塞(图8)，产生堵溃放大效应，对泥石流活动的影响较大。统计比通沟泥石流物源如表2所示。

赤担隆巴泥石流松散冰碛物源占其物源总量的大部分(表2)。沟源处左岸现代冰川冰碛物前缘高陡，易失稳进入沟道参与泥石流活动(图7)。赤担隆巴泥石流沟道宽缓，两岸山坡坡度相对较缓，植被茂密，因而沟道两岸的坡面物源及崩滑物源较少，少量沟道松散物质在泥石流的冲刷下参与泥石流活动(图9)。

表2 泥石流物源统计 104m3

图6比通沟沟源处冰碛物源

图7赤担隆巴沟源处冰碛物源

图8比通沟泥石崩滑物源

图9赤担隆巴泥石流崩滑物源与沟道物源

3.3 水动力条件

分析泥石流的水动力条件，应选取泥石流沟口及沟源处气象站所获的气象数据进行研究[10]。然而，研究区内区域气象站较少，仅在帕隆藏布河谷有零星分布，当前无高山气象观测站[11]。比通沟沟口建有区域气象站(站名：古乡U4027)，但由于其未完整获取泥石流暴发当天的气象数据，因此选择其附近的古乡U4023区域气象站所获取的数据进行分析(气象站位置如图1内标注所示)。

本次群发性泥石流前期的日平均温度及日降雨变化过程如图10所示。泥石流暴发前，区内温度持续下降，但接连出现强降雨，可见本次泥石流的主要诱发因素为强降雨。本次群发性泥石流暴发前发生一次降雨过程，尤其泥石流暴发前两日研究区内出现强降雨，沟口气象站收集到的累计降雨21.8 mm。研究区当前无高山气象站，参照前人根据古乡沟沟口观测站(海拔2 530 m)与沟源高山站(海拔3 750 m)降雨量观测数据建立的降雨随高程变化规律[12]，对比此次暴发泥石流的主要形成区(海拔3 500～4 000 m一带)与沟口高差，推测主要形成区前两日累计降雨约50 mm左右。充沛的前期降雨使得沟内松散物质含水量增加，趋近饱和，稳定性降低，为后续泥石流的暴发做好了水源与物源的准备。泥石流暴发当天，再次出现强降雨，沟口日降雨23.9 mm，推测形成区降雨达55 mm，短时强降雨使得沟内物源启动，暴发泥石流。

图10泥石流暴发前后日平均气温和日降雨

图11泥石流暴发当日小时降雨

图11为泥石流暴发当日小时降雨分布情况。9月5日6：30比通沟暴发泥石流，沟口气象站收集到的小时降雨仅0.5 mm，但凌晨2点沟口的小时降雨为3.3 mm，为泥石流暴发前当日最大小时降雨量。如前文所述比通沟的地形条件及物源条件特征，其沟道存在多处堵溃点，因而认为凌晨2点时，沟源处的短时强降雨导致冰斗内的冰碛物源启动，形成泥石流，但由于沟道狭窄，沟岸崩滑物源丰富，形成数十处堵点，泥石流运动过程中产生“堵溃循环”，泥石流物质逐级停滞累积，最终延迟至6：30泥石流龙头冲出沟口。因此，本次比通沟泥石流的激发小时雨量沟口观测值为3.3 mm，按照上文同样方法，推测主要形成区小时雨量约为7.6 mm。

9月5日11：30赤担隆巴暴发泥石流，沟口气象站收集到11：00的小时雨量为3.6 mm，推测主要形成区小时雨量为8.3 mm，与上文分析的比通沟泥石流激发雨量相近，认为当日11：00的短时强降雨导致赤担隆巴沟源处的冰碛物源启动，形成泥石流。由于沟道宽缓顺直，沟岸崩滑物源较少，未产生沟道堵塞，泥石流直接顺沟下泄，堆积于老扇之上。泥石流冲出沟口的时间延迟较短。

4 结 论

(1) 北岸比通沟泥石流运动特征以堵溃为主。因堵溃放大效应而形成较高的峰值流量，一方面使其龙头部分能够从深切的沟道中翻出，淤埋公路及工棚，形成灾害。南岸赤担隆巴泥石流则直接顺沟下泄，冲出出山口后淤埋发电设施及道路，形成灾害。

(2) 本次暴发的泥石流均拥有的良好地形及物源条件。比通沟流域形态呈“漏斗”状，利于水源、物源向沟道内聚集，同时也储存了大量的松散冰碛物源。主沟纵比降大，为泥石流运动提供了良好的地势条件。流通区沟道狭窄，崩滑物源发育，沟道内存在数十处堵点，使泥石流产生堵溃放大效应，增加泥石流峰值流量，提高其致灾能力。赤担隆巴沟源处发育有冰蚀围谷，同样储存有大量松散冰碛物源。冰蚀围谷坡度较陡，利于物源、水源聚集，启动泥石流。主沟纵比降较大，沟道宽缓顺直，沟岸崩滑物源较少，利于泥石流运动、下泄。

(3) 本次群发性泥石流的主要诱发因素为强降雨。泥石流暴发前两日出现强降雨，使得沟内松散物质含水量增加，趋近饱和，稳定性降低，为后续泥石流的暴发做好了水源与物源的准备。泥石流暴发当日再次出现强降雨，沟口气象站收集到的累计降雨23.9 mm，推测沟源处降雨达55 mm，集中强降雨直接激发了泥石流。经对暴发当日小时雨量进行分析，结合比通沟地形及物源特征，认为其由凌晨2：00的集中降雨启动，但由于沟道堵塞严重，龙头数小时后才到达沟口。赤担隆巴则由11：00的集中降雨启动，随后顺沟下泄，堆积于老扇之上。经分析，沟口观测的激发本次泥石流的小时雨量为3.3～3.6 mm，推测主要形成区激发泥石流的小时雨量约为7.6～8.3 mm。

致谢：特别感谢加拿大Laval大学方向博士给予本文的热情指导。