九寨沟县勿角镇甲勿沟泥石流运动参数特征与堵河分析

2023-12-18蒋清明郝红兵杨学之蒙明辉

地质灾害与环境保护 2023年4期

蒋清明,郝红兵,杨学之,蒙明辉

(1.四川省华地建设工程有限责任公司,成都 610081;2.四川省地质灾害防治工程技术研究中心,成都 610081)

泥石流是一种分布广泛,具有突发性、偶然性且破坏力极强的地质灾害。泥石流的形成有三大条件,即地形条件、物源条件、水源条件,在西部山区,地形主要以高山峡谷为主,为泥石流的形成创造了良好的地形条件;“5.12”地震对震区内山体扰动较大,新增了大量的松散固体物源,临沟临坎处稳定性较差的第四系松散堆积体易启动参与泥石流过程;随着西部大开发的进一步加强加深,在西部山区具有较强的人类工程活动,人类工程活动也新增了部分泥石流固体物源,进一步增强了泥石流的形成规模;近年来极端天气增多,强降雨天气持续增加,为泥石流的形成提供了较好的降雨条件。

由于泥石流易发性逐渐增加、破坏性极强等特点,前人对泥石流开展了较为深入的研究。许多学者主要从地形地貌条件、松散固体物质条件、植被发育条件及降水情况等多方面深入分析了泥石流的形成条件[1-3];也有部分学者从泥石流的运动特征方面进行了较为深入的研究[4];丰富的理论研究成果为泥石流的防治提供了理论依据及经验。

2020年8月15～17日,九寨沟勿角镇甲勿沟流域经历持续强降雨过程,根据附近勿角镇甲勿气象观测站数据可知,8月15日降雨量为21.9 mm,8月16日降雨量为38.2 mm,8月17日降雨量为37.3 mm,3 d降雨量总和达到97.4 mm,单日降雨量接近5 a一遇降雨强度,在连续经历强降雨天气后,甲勿沟于8月17日凌晨发生大规模泥石流过程,泥石流持续时间近3 h,沿途冲毁旅游道路,淤埋村民房屋,泥石流沟口处绵九高速公路项目部几乎被完全淤埋,直接经济损失超过4 000万元,沟口处冲出的固体物质方量约20.0×104m3。

受九寨沟“8.8”地震影响[5],再加上甲勿沟本身有利的地形条件、物源条件等[6],甲勿沟泥石流本身属于易发性较强的泥石流沟[7,8],一旦甲勿沟再次发生大规模泥石流过程,将直接威胁沟域内的熊猫园、甲勿池景区及游客中心、旅游道路等旅游设施,除此之外还威胁沿沟的居民区、沟口的勿角场镇以及绵九高速公路等,其危害性较大[9]。因此,对泥石流的基本特征、运动参数特征及堵河可能性等进行研究,对泥石流的防治具有重要的参考意义。

1 流域概况

甲勿沟泥石流位于九寨沟县勿角镇甲勿村,属汤珠河左岸一级支流,流域面积43.9 km2,主沟全长10.7 km,最高点位于流域西南侧山脊部位,高程3 335 m,最低点位于主沟与汤珠河交汇处,高程1 780 m,相对高差1 555 m,沟道平均纵坡降145.9‰。流域属强烈切割区,总体地势西南高、东北低,河谷深切,地形起伏大。山坡坡度一般在20°以上,部分地区坡度在50°以上,总体上具有岸坡陡峻,切割深度大,呈现上陡下缓的空间特征。甲勿沟流域地形呈叶脉状,共发育有9条主要支沟,其中二道坪沟、6#支沟为泥石流沟,其余7条支沟均为清水沟(图1)。

图1 甲勿沟流域概况图

2 泥石流分区特征

根据泥石流形成条件和松散固体物源的分布情况,将甲勿沟沟域划分为清水区、形成流通区、堆积区等三大片区[10]。

2.1 清水区特征

清水区主要分布于主沟2 460 m高程以上至沟源头部位。面积约32.04 km2,主沟2 460 m以上至沟源段长2.3 km,最大高程3 335 m,相对高差875 m,平均纵坡降380.4‰。区内沟谷形态多呈“V”型槽谷;沟道岸坡陡峻,岸坡坡角35°～60°。清水区森林植被发育较好,保水性好,不良地质现象较弱,主要为泥石流的形成汇集水源和提供水动力条件。

2.2 形成流通区特征

形成流通区主要分布于沟域中下游,流域面积约11.73 km2,分布高程1 815～2 460 m,相对高差645 m,该段沟道长度约5.2 km,平均纵坡降为124.9‰。该区沟谷形态多呈“U”型槽谷,沟道上游局部地段呈“V”形峡谷(见图2);岸坡坡角在20°～60°之间,沟底为第四系冲洪积层及前期泥石流堆积层,岸坡局部基岩裸露,出露地层主要为生物碎屑灰岩及石英岩状砂岩等。该区发育有较多崩塌和滑坡等不良地质现象,为泥石流提供了大量松散固体物源;形成流通区沟道较长,其冲淤特征视不同沟段的弯曲程度以及纵坡降差异而表现出不同的特点,纵坡较大处表现为以冲为主,局部稍缓部位以冲淤平衡为主[11]。

图2 形成流通区地形地貌

2.3 堆积区特征

堆积区主要分布于主沟高程1 815 m以下至主河汤珠河段,该段主沟沟长约0.61 km,沟道高程1 770～1 815 m,沟宽50～100 m,平均纵坡降73.8‰,流通堆积区面积约0.13 km2。该区段沟道平均纵坡比降变缓,沟床变宽,水动力条件减弱,沟道冲淤特征以淤积为主[11]。

3 泥石流运动特征参数计算

泥石流运动特征参数计算主要从形态调查法以及雨洪法两个方面入手,既可验证泥石流运动特征参数计算的准确性,也可为泥石流的堵河分析提供依据。形态调查法根据现场试验、现场实测泥痕、泥位及沟道断面等进行计算,雨洪法计算参数主要根据《四川省中小流域暴雨洪水计算手册》及相关规范[12]中阐述的方法进行计算得到。

3.1 泥石流容重

(1) 配浆法重度取值

甲勿沟泥石流近期发生过泥石流过程,根据泥石流流体性质,在沟道上中下不同部位共进行了12次重度实验,配浆法重度计算公式见泥石流勘查规范[12],此处不再进行赘述,其结果见表1。

(2) 雨洪法重度取值

雨洪法重度取值是结合配浆法以及查表法进行综合取值。查表法是在现状调查基础上,综合了流域内地质环境条件、沟道特征及物源发育分布等特性,带预测性,更能反映泥石流特征及实际情况。泥石流过后沟域内局部物源发生了启动,物源特征、启动条件及堵塞程度在一定程度上发生了变化,采用配浆法及查表法确定泥石流重度参数均存在误差。因此,雨洪法重度取值考虑不同断面现场实际及查表法等综合分析取值确定,结果见表2。

3.2 泥石流流速

由配浆法及查表法结果可知,泥石流属于稀性泥石流,对流速的计算采用原铁道部第二勘察设计院推荐的西南地区经验公式(式1)进行计算[12]。

(1)

式中,Vc为泥石流断面平均流速;1/n为清水河床糙率系数,查水文手册;R为水力半径,一般用平均泥深代替;I为泥石流水力坡度,一般按沟床纵坡代替;φ为泥沙修正系数;γH为泥石流固体物质容重。

根据公式及各断面部位沟道特征,形态调查法及雨洪法10 a一遇降雨强度(P=10%)各断面部位泥石流流速计算结果见表3。

通过形态调查法及雨洪法泥石流流速计算结果对比可知,甲勿沟“8.17”泥石流流速接近雨洪法10 a一遇降雨强度下的流速,说明甲勿沟“8.17”泥石流的强度接近10 a一遇降雨条件下的泥石流强度。

表2 泥石流重度综合取值统计表

表3 泥石流流速计算表(形态调查法及雨洪法)

3.3 泥石流峰值流量

(1) 形态调查法

采用形态调查法对泥石流流量进行计算[12],计算公式为:

Qc=WcVc

(2)

式中,Qc为泥石流断面峰值流量(m3/s);Wc为泥石流过流断面面积(m2),采用调查测绘获取的沟道形态参数;Vc为泥石流断面平均流速(m/s)。

据此求得形态调查法各断面位置泥石流峰值流量,计算结果见表4。

表4 甲勿沟“8.17”泥石流流量计算统计表(形态调查法)

(2) 雨洪法

暴雨洪峰流量计算根据相关规范[12]及计算手册[13]公式计算地表水汇水流量,计算公式为:

(3)

式中,Q为最大洪峰流量(m3/s);ψ为洪峰径流系数;i为最大平均暴雨强度;s为暴雨雨力;n为暴雨公式指数;F为集水面积(km2);τ为流域汇流时间(h)。

暴雨洪峰流量计算结果见表5。

在计算求得暴雨洪峰流量的基础上,根据规范[12],采用下式(式4)进行泥石流峰值流量的计算:

Qc=(1+φ)QpDc

(4)

式中,Qc为泥石流断面峰值流量(m3/s);φ为泥沙修正系数;Qp为暴雨洪峰流量;Dc为堵塞系数,按规范查表确定。

各断面泥石流峰值流量计算结果详见表5。

表5 泥石流暴雨洪峰流量与峰值流量计算结果(雨洪法)

通过形态调查法及雨洪法泥石流峰值流量计算结果对比可知,甲勿沟“8.17”泥石流峰值流量比雨洪法10 a一遇降雨强度下的峰值流量偏小较多,如主沟服务中心上游形态调查法流量为80.261 m3/s,雨洪法流量为115.55 m3/s,又如主沟沟口形态调查法流量为127.92 m3/s,雨洪法流量为161.79 m3/s,所有计算断面形态调查法计算结果均比雨洪法计算结果小。推断主要原因是因为现场调查时能准确测量最高泥位线,但由于泥石流后期强度逐渐变弱,泥石流的携砂能力减弱,沟道内逐渐淤积有固体物质,导致现场测量的过流断面偏小,进而导致形态调查法计算得到的泥石流峰值结果偏小。

4 泥石流堵河分析

甲勿沟沟口为汤珠河,“8.17”泥石流过程中,大量固体物源进入主河,造成主河的堵塞以及河床的抬升,主河河水以及甲勿沟洪流在交汇处形成漫流之势,大量河水越过河堤,河水携带大量泥砂及碎石,淤埋了阳山村安置点以及勿角场镇,淤埋深度平均达0.5～1.0 m,最大深度处达到2.0 m(图3)。

针对泥石流堵塞主河问题,陈德明[14]、唐川[15]、张金山[16]等国内许多学者对其进行了相应研究,提出了泥石流堵河的经验判别式。

参照陈德明等提出的泥石流堵河可能性的经验判别式,对可能堵河阈值进行计算。

(5)

式中,γ支为泥石流流体重度(t/m3);γ主为主河河水重度(t/m3);Q支为泥石流峰值流量(m3/s);Q主为泥石流发生时主河河水流量,取汤珠河洪期平均流量(m3/s);V支为泥石流入汇时的流速(m/s);V主为泥石流入汇处主河河水流速(m/s);α为泥石流沟入汇角(°);Cγ为阈值,为1.44。

前文第3节对泥石流运动特征参数进行了计算,根据泥石流流量、流速等计算结果,采用判别式进行堵河阈值计算,计算结果见表6。