唐家山堰塞体渗流稳定及溃决模式分析

在准确获取唐家山堰塞体地质结构(从上到下依次为①层碎石土、②层块碎石和③层似层状碎裂岩3层结构)及相关渗透参数基础上,采用Visual Modflow可视化三维软件,模拟4种堰塞湖水位(710,720,730和740m)条件下不同土层的渗流速度和渗透坡降,表明由于②,③两层颗粒粗大、渗透性好,总体表现出稳定流的渗水特点,整体坡降变化稳定,中间不会出现如低渗透黏性土的紊流状态、渗透坡降出现拐点以及管涌渗透特点。而下游侧①层碎石土由于沿第③层形成的贯通性渗流在下游侧穿越该层渗出形成的最大坡降会大于允许坡降,将导致堰塞体下游发生零散或局部渗流破坏,且临界堰塞湖水位为726m。随堰塞湖水位抬升,整个堰塞坝体出现的溃决模式为:下游侧表层碎石土层因渗透破坏和溯源侵蚀,同时因进口段地表水流漫顶淘刷,最终导致上部第①层碎石土被侵蚀、淘刷带走。随着第①层被淘刷、水流速度加大又进而会带动第②层块碎石被逐渐冲刷下切,但不会发生整体溃决,而第③层似层状碎裂岩将保持稳定,侵蚀和淘刷的下限深度就是第③层似层状碎裂岩顶部。

唐家山堰塞湖应急除险技术实践

简要介绍了唐家山堰塞湖应急除险工程概况,归纳了堰塞湖有关地形地质、水情等资料的快速获取与应用技术,从形成条件、形成过程等方面对堰塞坝形成机制进行了分析,对堰塞坝做了安全性分析评估,并对溃坝洪水进行了分析计算;概要论述了唐家山堰塞湖应急除险总体方案、开漕引流方案及除险效果。

堰塞坝溃坝数学模型研究与应用

针对堰塞坝坝体土石料的宽级配特性,引入与水流方向垂直的附加作用力来考虑粗颗粒对细颗粒的阻拦、遮蔽作用以及细颗粒对粗颗粒的包围、填实作用,提出了一个可模拟堰塞坝漫顶溃决过程溃口发展规律与流量过程的数值模型和相应的计算方法,利用该模型对唐家山堰塞坝泄流过程进行了模拟,得出的泄流槽发展规律与洪水流量过程与实测资料接近,验证了该模型和计算方法的合理性。进一步,利用笔者建议的数学模型及数值计算方法,比较分析了唐家山堰塞坝除险过程中泄流槽断面型式对堰塞坝泄流过程的影响,发现堰塞湖在采用泄流槽引流除险时,泄流槽深度与断面型式对其泄流过程具有重要影响,增加泄流槽深度,可明显提高泄流效率,但堰塞湖下游将承受更大的风险。对于同样深度的梯形泄流槽,如果将槽底部断面减小,形成复合梯形泄流槽,不仅可减少开挖工作量,而且没有明显降低泄流效率,同时后者的泄流过程更为平缓,最大洪峰流量减小,出现的时间滞后,堰塞湖下游承受的风险也将降低。

土石坝漫顶溃决模型的分析比较

为解决由现场溃坝过程记录资料匮乏导致的大量土石坝漫顶溃坝模型得不到工程验证的问题,在分析总结17种土石坝漫顶溃决机理模型基础上,以唐家山堰塞坝1/3溃坝事故为工程背景,选取BREACH模型、DB模型和IWHR模型对其溃坝过程进行数值模拟,研究这3种模型计算得到的筑坝材料冲蚀速率、流量过程线和水位变化曲线与实测值的差异。应用BREACH模型研究了筑坝材料的黏聚力、内摩擦角、不均匀系数和孔隙率对流量过程线的影响。结果表明:(a)IWHR模型预测的洪峰流量最大,BREACH模型次之,DB模型最小;(b)与实测数据相比,DB模型和IWHR模型预测的峰现时间提前,BREACH模型预测的峰现时间推迟;(c)随着黏聚力和内摩擦角的增大,计算峰值流量减小,峰现时间延后;(d)随着孔隙率和不均匀系数的增大,计算峰值流量增大,峰现时间提前;(e)不均匀系数和孔隙率对溃坝计算结果的敏感性尤为显著。

几种溃坝模型在溃决洪水模拟中的适用性比较

使用一维溃坝模型对唐家山堰塞湖溃决过程进行模拟,研究了推移质公式和边岸侵蚀模型对模拟结果的影响。溃决洪水流量过程线、溃口展宽过程和床面冲刷过程的对比表明,Meyer-Peter&Muller公式在描述漫顶溃坝初始阶段时存在计算冲刷偏小的现象,其可能原因是输沙率与Shields数的高阶(1.5次)非线性关系。本文的比较结果,有助于提高溃坝条件下推移质公式和边岸模型适用性的认识。

唐家山堰塞湖泄洪问题研究

研究了唐家山堰塞湖的泄洪规律,并对其蓄水、溃坝、泛洪等三个过程进行了模型分析.