1. 引言
泥石流是山区常见的地质灾害之一,具有暴发突然、高流速及破坏力巨大等特点 [1] ,给当地人民群众的生命安全造成极大的威胁,同时也给当地的生态环境造成破坏,严重阻碍当地的经济发展和生态建设 [2] 。而风景区内的泥石流威胁更大,造成的后果也更严重。2005年8月11日,我国四川省境内磨西河流域内群发泥石流,泥石流不仅冲毁房屋、公路和水电站,更造成1200余名游客被困在海螺沟景区内,造成数千万元的经济损失 [3] ;泥石流发生之后,整个景区基础设施陷于瘫痪,景观和生态环境遭受严重破坏,游客大幅减少,旅游收入蒙受重大损失 [4] 。对海螺沟内泥石流开展研究并提出防治建议,对其防灾减灾具有重要意义。海螺沟是我国著名的风景区,风景区内泥石流的防治不同于威胁城镇、交通泥石流的防治,不仅要有效控制泥石流,还要达到保护景观、保护游客的目的,治理工程要与景观相协调 [5] 。谭万沛在20世纪90年代提出海螺沟地质灾害的防治对策主要为加强宣传、树立警示标记、加强工作人员防灾培训、合理选择游览时间、建立监测网络等措施 [6] ;何德伟等提出制定景区防灾规划、建立泥石流预测预报和预警系统、设置警示标识、系统开展治理和生态保护 [4] ;倪化勇提出了加强游客和工作人员防灾避灾宣传、树立警示牌、沟口跨越工程、热水沟和黄崩溜小沟修建排导槽或堡坎 [7] ;晏鄂川等对海螺沟3个营地附近的地质灾害问题分别提出了防治对策,其中热水沟的治理方案是疏导 [8] ;倪化勇研究了磨西河流域泥石流爆发的临界雨量值,为该区域泥石流的降雨监测提供了降雨阈值 [9] 。根据实地调研发现,海螺沟内群测群防网络健全、工作人员防灾意识强,有部分治理工程和专业监测,但对游客的防灾减灾宣传有待加强;热水沟有专业监测设备但未进行治理,二号营地处于半营业状态。本文通过对2010年热水沟泥石流动力学参数计算,结合2013年泥石流调查成果,发现热水沟泥石流以小规模为主,但沟内滑坡堵溃会加大泥石流规模和破坏力,综合前人的防治经验,提出专业监测、工程治理、加强宣传相结合的防治模式。
2. 流域背景概况
热水沟流域内海拔最高点为5888 m,最低点为2425 m,相对高差为3463 m,主沟床平均纵比降为343‰;山坡坡度一般在30˚~36˚,最大达53˚;主沟以“V”型谷为主,沟宽3~8 m;是一条典型的极高山陡坡大比降小流域沟谷。根据其地貌及泥石流物源补给特征,将热水沟划分为冻融区、形成区、流通区和堆积区。冻融区主要为海拔3500 m以上的区域,沟道比降为589‰,植被稀疏,岩体破碎,崩落现象突出。形成区高程范围2625~3500 m,沟道比降为350‰,沟道两岸冰碛物大面积分布,受沟谷下切影响,冰碛体发生较多的崩塌、滑坡。流通区高程范围2440~2625 m,沟道比降为170‰,沟道较宽阔,相对平缓,略曲折,海螺沟二号营地及景区公路位于该区域。堆积区高程范围2425~2440 m,堆积扇中等大小,堆积体坡度平缓,约为10˚ (如图1)。
流域内岩浆岩集中在中上游,为斜长花岗岩混合岩(γo52),岩体风化严重,节理发育,表层破碎,冰雪崩积锥分布较多。变质岩集中在中下游,为下二叠统中段(P12)灰岩、白云岩、石英片岩,下二叠统上段(P13)大理岩夹云母石英片岩,上二叠统(P2)绿片岩、石英岩夹板岩、大理岩,基岩大面积出露,残坡积物厚度不大。泥石流堆积物(Qhapl+df)集中在沟口堆积扇和下游沟道中。冰碛物(Qhgl)集中分布在距沟口约2 km的沟道两岸和沟口处;冰碛物较松散,发育5个大滑坡,均处于活动状态,是热水沟泥石流固体物质的主要来源 [8] 。目前热水沟松散堆积物的总储量估算约为195 × 104 m3,强活动物源储量约为9 × 104 m3,具备再发泥石流的物源条件。
3. 泥石流动力学特征
2010年8月热水沟爆发小规模泥石流,泥石流冲入海螺沟,冲毁部分温泉泳池。根据现场调查结果,对该次泥石流的流速、流量、撞击力、一次泥石流总量和固体物质总量进行计算。
3.1. 流速
热水沟泥石流爆发时,水石混杂,似山洪爆发,浆体较稀,为稀性泥石流,采用西南地区的稀性泥石流流速计算公式 [10] 计算其流速:
(1)
热水沟取1/n清水河床糙率系数为10,水力半径R用平均水深1.5 m代替,水力坡度I用流通区沟床纵坡169.7‰代替,泥石流固体物质比重rH取2.65 t/m3,泥石流泥沙修正系数
为0.65。则热水沟泥石流流速为6.52 m/s。
Figure 1. Layout of debris flow monitoring equipment in Reshuigou
图1. 热水沟概貌和泥石流监测设备布置图
3.2. 流量
泥石流流量计算方法有两种,一是雨洪法;二是形态调查法。将形态调查法的结果与雨洪法结果进行比较,可以确定一次泥石流的发生频率。
1) 雨洪法
假设泥石流与暴雨同频率、且同步发生,按下式 [10] [11] 计算泥石流流量Qc:
(2)
其中:Qp为暴雨洪峰流量,DC堵塞系数,Φ为泥砂修正系数。
计算结果如表1所示。
2) 形态调查法
调查泥石流泥位和过流断面,根据公式 [10] 求泥石流断面峰值流量Qc:
(3)
其中:Vc为泥石流流速,WC为过流断面面积。
用上式计算的流通区上段和下段泥石流峰值流量Qc分别为124 m3/s和193 m3/s。
3.3. 撞击力
泥石流撞击力是泥石流防治工程设计的重要参数。
泥石流体整体冲压力计算公式 [10] :
(4)
泥石流中石块的冲击力的计算参照以下公式 [10] :
(5)
其中:λ为建筑物形状系数,γc为泥石流重度,α为受力面与冲击方向夹角,VC为泥石流流速,g为重力加速度,r为动能折减系数,Q为石块质量,C1 + C2为弹性变形系数。
通过计算,热水沟泥石流的整体冲压力
约为92.5 kPa,泥石流中石块的冲击力Pd为9.8 kPa。
3.4. 一次泥石流总量及输砂量
计算法根据泥石流历时T(s)和最大流量Qc (m3/s),按泥石流暴涨暴落的特点,将其过程线概化成五角形计算。按下式 [10] 计算一次泥石流总量Q(m3)。
(6)
F = 10~100 km2,K = 0.0378。
一次泥石流冲出的固体物质总量QH (m3):
(7)
根据计算,得到形态调查与不同频率下的一次泥石流过程总量与冲出固体物质总量,见表2。
Table 1. Debris flow discharge based on storm flood peak
表1. 雨洪法计算泥石流流量
Table 2. Volume and total solid matter of debris flow
表2. 一次泥石流过程与固体物质总量表
4. 泥石流防治对策
4.1. 监测预警
通过实地调查,海螺沟二号营地和景区公路位于泥石流危险区,热水沟泥石流已对二号营地的正常运营产生影响。2010年在中国地质调查局项目的依托下对该沟进行监测预警,热水沟泥石流的监测方案为雨量监测、地声监测和视频监测相结合的综合监测方案(图1),雨量计、视频设备和地声警报设备置于二号营地,无线地声遥测仪安装于二号营地上游约1.5 km处沟道内,并于2010年底完成设备安装。
雨量监测主要是预报特定降雨强度下可能会发生泥石流灾害 [12] ;视频监测了解泥石流在发生过程中的运动特征 [13] 。地声监测主要是对泥石流进行提前预警 [14] 。
经监测,2013年7月12日23点30分热水沟内滑坡堵溃形成中等规模泥石流,方量约12万m3,无人员伤亡,二号营地内的部分建筑被破坏和於埋。无线地声监测仪提前约2分钟发出警报,成功预报了该次泥石流,获取的监测资料为该地区泥石流的科学研究提供了丰富的资料 [15] 。
4.2. 防治建议
1) 加强监测:群专结合的监测是预防泥石流的有效措施,我国地质灾害的群测群防具有许多成功的案例和经验。热水沟内群测群防网络健全,应该继续予以加强;而专业的监测更是必不可少,目前的监测方案取得了较理想的效果;由于沟内滑坡堵溃易形成较大规模的泥石流,可在原有监测方案的基础上,增加滑坡监测。
2) 加强对游客防灾避灾知识的宣传:海螺沟是著名的旅游景区,每年接待游客30万人次,旅游旺季同时是泥石流频发的季节,而游客大多数防灾避灾意识薄弱,而二号营地正处于热水沟内,受泥石流威胁大,因此需加强对游客防灾避灾的宣传。除树立标识标牌外,还应印发宣传册,给每位入住二号营地的游客发放一份;导游为游客讲解自然风光的同时还应告知泥石流的防灾减灾知识。
3) 采取一定的工程措施:热水沟内河道起伏大,局部堵塞,应该进行疏通与清理。由于河道较窄、二号营地又紧邻河道,滑坡堵溃后泥石流规模较大,仅靠排导槽难以达到理想的治理效果,因此建议同时修建排导槽、谷坊坝等工程措施进行治理,并修建防护堤保护二号营地。
5. 结论
1) 热水沟位于海螺沟风景区内,泥石流对二号营地和景区公路造成较大的威胁,沟内物源丰富,近年来泥石流频发。该沟泥石流一般规模较小,流速较大,冲击力较大,沟内滑坡堵溃会加大泥石流规模和破坏力,研究其动力学参数对该沟的防灾减灾具有重要的科学意义。
2) 由于热水沟处于海螺沟风景区内,威胁二号营地和公路,建议在现有的专业监测基础上增加沟内滑坡监测,工程治理应以排导为主辅以稳、拦,同时还应采用多种方式加强对游客防灾避灾知识的宣传。
3) 目前,在热水沟开展了雨量、视频和地声监测预警手段,积累了大量的监测数据,为开展热水沟泥石流形成机理的研究提供科学依据。
基金项目
国家自然科学基金项目“地震扰动的泥石流源区砾石土强度恢复动态过程实验研究”(项目编号:41502337),中国地质调查局地质调查项目“四川西昌等城镇和金沙江下游水电开发区工程地质调查”(项目编号:DD20160722),四川省国土资源厅联合科技攻关项目“四川地质灾害综合防治关键技术与方法研究”(项目编号:0773-1641GNSC01693)。