论文部分内容阅读
随着全球气候变化异常和人类社会经济活动的不断增强,极端天气和地质事件增多,尤其是汶川地震后的4年里泥石流活动非常活跃,暴发了数次特大规模泥石流,汛期的实践表明,泥石流是震后最为活跃的灾种,泥石流堵河成为一种频繁的灾害事件。本文面对极端激发条件和汶川地震震后成灾环境急剧演变状态下,我国特大泥石流灾害频发的严峻形势,面向国家山区建设中急需解决的泥石流防灾减灾问题,以岷江映秀至汶川段多处泥石流沟堵江,以及舟曲、清平特大泥石流灾害为背景,采用现场调查、理论分析、模型试验和数值模拟相结合的技术手段,研究泥石流堵塞大河的特点及机理,建立泥石流堵塞大河各种模式的判据等,为城镇和重大基础设施防御泥石流堵河引发的大范围灾害提供科学依据。主要进行了以下工作并取得了研究成果:(1)基于现场调查和查阅文献收集的资料,归纳出泥石流堵塞大河存在着潜坝壅堵、堰塞坝局部堵塞和堰塞坝全堵三种堵河模式,并提炼出了三种堵河模式的特点及成灾模式。潜坝壅堵主要发生于运动速度较高的稀性泥石流或者容重较小的粘性泥石流,其特点是泥石流始终在水下运动,泥石流停积体淤高河床,改变河道纵断面,以上游壅水形式造成淹没灾害。以堰塞坝形式造成大河堵塞主要发生于粘性(或高浓度)泥石流。一般情况的局部堵塞,主要是堰塞坝束窄水流、改变河道,影响大河行洪能力;但当堵塞系数大于0.8以后,由于泥石流严重推移大河水流,可能造成异岸冲刷灾害。泥石流完全堵断大河的事件不多,但一旦发生,堰塞湖会造成上游淹没灾害,溃坝时会造成下游超常洪水;在宽谷段甚至会迫使主河改道,形成更大范围的严重灾害。(2)针对岷江流域泥石流多半是以堰塞坝形式局部堵塞主河的现象,运用大型室内水槽试验再现了交汇角为900时粘性泥石流入汇主河的过程,揭示了泥石流入汇主河过程的特点及泥石流堵塞大河的机理。重点研究了堰塞坝局部堵塞大河的特点,试验结果表明:泥石流入汇主河的运动过程是后续泥石流逐步推移龙头运动的过程;泥石流总量对堵河有重要的影响,总量越大堵河越严重;在泥石流堵河过程中,主河流速、堰塞坝物质流失量和上游水位壅高值都随堵塞系数呈非线性规律:主河流速随R先增大后减小,当R小于0.8时,堰塞坝的物质流失量少,上游水位壅高值呈线性增加,当R大于0.8时,堰塞坝的物质流失严重,上游水位壅高值呈指数增加。(3)根据泥石流堰塞坝堵河模式的特点,基于水沙动力过程,以泥石流在坡度变化沟道内的运动方程为基础,考虑堰塞坝和主河水流的相互作用,给出了泥石流在主河中运移距离的计算公式;鉴于泥石流堵塞大河过程中物质流失严重的现象,根据水槽实验数据确定了流失系数;以周必凡提出的泥石流堵断主河总量的计算公式为基础,提出了按流失系数加以修正的泥石流阻塞大河物质总量条件。从而建立了基于动力学方程、物质总量的泥石流堵塞大河综合判据。(4)基于堵塞系数R对于预估泥石流堵河灾害效应的关键作用,以模型试验数据为依据,导出了堵塞系数R的表达式,此式可用来预估泥石流堵河的程度,也可作为泥石流局部堵塞大河的另一判据,可用来指导泥石流防治工程的设计。(5)根据潜坝壅堵模式的特点,基于异重流理论导出了泥石流龙头运动简化方程,基于有坎宽顶堰模型给出了潜坝壅堵上游水位计算式;对影响泥石流龙头运动的各个因素进行了单因素分析,分析表明泥石流粘滞系数和切应力对龙头运动速度影响不明显,绕流阻力系数对龙头运动速度影响最大;给出了泥石流龙头在主河里的运移距离随时间的变化曲线,在此基础上给出了潜坝堵河的最小固体物质总量计算式。(6)根据泥石流堰塞坝堵河模式的成灾模式,运用数值模拟方法分析主河在不同堵塞系数时的水流流场,再基于信息熵原理定量地分析主河流场有序度,得出当堵塞系数R为0.80-0.9时,主河流场熵值急剧下降,表明流场有序度不断增高,流向以偏向异岸为主,且此时水流速度大,水流偏转又集中发生在近河岸处,对河岸具有很强的冲刷能力,R=0.80可作为主河水流偏转冲刷异岸的临界点,把堵塞系数R>0.80作为泥石流局部堵河异岸受灾判据。通过岷江流域银杏至映秀段泥石流堵河实例验证,表明数值模拟结果与实际情况吻合,泥石流局部堵河异岸受灾判据可行,进而提出了泥石流灾害效应与判据框图。基于堰塞坝局部堵塞主河使上游水位壅高致灾的现象,给出了水位壅高值计算式,为计算淹没范围提供了依据。(7)针对不同堵河模式的灾害效应,提出相应的减灾对策。对于全堵模式首先考虑抬高淹没区线路标高,若受其他条件限制无法实施时,则考虑限制泥石流入汇主河的流量和总量,以减轻或消除泥石流堵塞大河。对于局部堵塞模式主要考虑异岸受灾,采用导流工程、支挡工程等工程措施保护异岸线路工程免受冲刷。对于壅堵模式,则考虑小范围抬高线路标高。