对唐家山堰塞湖溃坝洪水计算的主要认识_徐照明

合集下载

唐家山堰塞湖排险的探索与思考

唐家山堰塞湖成功排险的探索与思考摘要:在交通、通讯完全中断的情况下，在地质情况不完全了解、施工环境极其恶劣的条件下，在时间紧迫、任务繁重的前提下，成功排除极高危险堰塞湖险情，创造了世界奇迹。

本文从堰塞湖排险的设计组织、统筹规划、施工管理等方面进行了简要介绍，并进一步对堰塞湖应急排险进行探索与思考。

主题词:堰塞湖排险施工探索与思考1、概况1.1堰塞湖的定义堰塞湖是指地震后引起的大规模山体滑坡，河水冲击泥土、山石而造成堆积，堵截河谷或河床后贮水而形成的湖泊。

堰塞湖形成过程有四个过程，一是原有的水系；二是原有水系被堵塞物堵住。

堵塞物可能是火山熔岩流，可能是地震活动等原因引起的山崩滑坡体，可能是泥石流，亦可能是其他物质。

三是河谷、河床被堵塞后，流水聚集并且往四周漫溢。

四是储水到一定程度便形成堰塞湖。

1.2唐家山堰塞湖成因唐家山堰塞湖形成过程是汶川大地震造成了龙门山300KM的地震断裂带，山体坠落，封闭河道，原有水系被堵塞物堵住。

河谷、河床被堵塞后，流水聚集并且往四周漫溢，水位迅速积蓄抬高，造成洪峰，威胁下游。

1.3唐家山堰塞湖规模及险情唐家山堰塞湖位于北川县上游通口河上，堰体长803.4m，宽611.8m，面积30.72万m2，体积2037万m3，堰体最高点高程791.9m，堰体低高程667.5m，最大坝高124.4m，垭口处坝高82.65m。

总库容3.15亿方，回水长度约23km。

若产生溃坝将影响下游北川县，江油市，绵阳市的涪城区、科学城、农科区以及三台县130多万人口，甚至影响到遂宁市，属此次特大地震中形成的极高危级堰塞湖。

2、设计方案与施工规划2.1设计方案及方案调整2.1.1最初抢险方案。

利用地震形成的垭口，采用工程机械开挖明渠，降低水位。

有“高、中、低”三个泄流渠方案，低标准方案，泄流渠进口高程747m、渠底宽28m，泄流漕总长670m，上游平缓段纵坡0.6%，下游陡坡段纵坡分别为24%和16%；中标准方案，泄流渠进口高程745m、渠底宽22m，泄流漕总长680m，上游平缓段纵坡0.6%，下游陡坡段纵坡分别为24%和16%；高标准方案，泄流渠进口高程742m、渠底宽13m，泄流漕总长695m，上游平缓段纵坡0.6%，下游陡坡段纵坡分别为24%和16%。

堰塞湖逐渐溃决洪水模拟及溃口变化影响分析

{bm
=
W3Eρ，ρ
=
KH 100
（ 1）
K = φW －0． 577
式中，W 为总溃坝水量，万 m3 ； E 为每米长度坝的方量，
m2； K 为冲刷系数； H 为坝体高度 m； φ 为与坝身材料及
密实度有关的系数，由表1 查出。其中溃口口门宽 bm 还
受坝址断面的河谷宽度限制。
表1 φ值
Tab． 1 Values of φ
土类
材料及密实度
φ
1 较松的均质土坝，压实不好的辗压、冲填、水中倒土坝 12． 500
2
质量较好的均质土坝
6． 700
3
坚硬的均质土坝，土料较多的土石混合坝
3． 650
4
土料较少的土石混合坝，质量较差的堆石坝
1． 680
5
较密实பைடு நூலகம்堆石坝
0． 495
1． 1． 2 峰前流量过程
峰前流量过程的界定首先是分析峰现时间。谢任
纵观历次堰塞湖溃决事件，堰塞坝以土石结构居多，溃决形式多数以渐溃方式存在。渐溃方式下溃坝洪水计算的关键问题是其溃口参数及演变过程的确定。但迫于坝体土质特征差异，溃口形态差别较大。当前不少学者提出应用泥沙冲刷理论进行溃口变化过程动态仿真，已得到实际运用。如张建国等人采用侵蚀崩塌的侧向扩展溃决理论，基于坝体土质参数建立了溃口侵蚀扩展模型，溃口变化过程与实测较为吻合［2］。郭非凡在土石坝溃坝 BＲEACH 模型基础上，引入基于能量平衡原理和刚体极限平衡法的边坡稳定公式，计算溃口间歇性横向扩展，总体效果也较良好［3］。徐照明等还通过构建 MIKE11 溃坝洪水模型［4］，从泥
沙输移特性考虑，对唐家山堰塞湖不同溃决历时、溃口形状及溃口发展过程情况下的溃坝洪水进行了模拟，总体结果令人满意。与此同时，还有不少学者利用统计学理论整理出溃坝溃口形状、溃决时间的经验公式，如黄河水利委员会经验公式［5］、铁道科学研究院经验公式［6］、谢任之经验公式等［5］，这些经验公式结构相对简单，可较容易计算出溃口最终形态。如黄毅等人在 DAMBＲK 溃口理论的基础上6］，应用黄河水利委员会及铁道科学研究院的经验公式确定了溃口参数，并以板桥水库溃坝洪水为例进行验正，认为计算效率较高。谢任之经实际数据分析所提出的溃口计算公式及洪水计算方法［8］，同样在板桥水库溃坝洪水模拟上具有较高精度。虽然基于泥沙冲刷理论的溃口计算方法可较精细化地模拟溃口演变过程，但需要的土质参数较多，过程也较烦琐，在资料稀缺及应急条件下往往略显乏力。而经验公式相对而言计算简便，模型参数少，但不可避免的是多数参数为经验取值，适用性仍有待进一步验证。

唐家山堰塞湖溃坝洪水分析及泄流冲刷模拟

第３９卷第２２期２００８年１１月人民长江ＹａｎｇｔｚｅＲｉｖｅｒＶ０１．３９．Ｎｏ．２２ＮＯＶ．，２００８文章编号：１００１—４１７９（２００８）２２—００７９—０４唐家山堰塞湖溃坝洪水分析及泄流冲刷模拟朱勇辉范北林卢金友张细兵杨文俊渠庚（长江水利委员会长江科学院，湖北武汉４３００１０）摘要：土石坝渍坝数学模型ＢＲＥＳＺＨＵ建立在渍坝试验及原型渍坝案例中所观察到的溃坝机理基础之上。

模型先后用不同国家的多组溃坝试验资料进行率定和验证，并被成功应用于原型堤坝溃决案例的模拟，结果良好。

“５·１２”地震唐家山堰塞湖险情发生后，其不断升高的水位和不断增大的湖容给下游百万群众的生命财产安全造成了巨大威胁。

运用ＢＲＥＳＺＨＵ模型并结合坝下游溃坝洪水演进模型针对堰塞湖上游可能出现的不同频率洪水、坝体的不同溃决方案和不同溃决过程等对数１０种工况下唐家山堰塞湖的调洪、溃坝及洪水传播过程进行了计算与分析，为抢险方案和应急预案的制定提供了有力的技术支持。

险情结束后运用ＢＲＥＳＺＨＵ模型及时对湖水下泄过程中泄流渠断面发展及坝址处洪水过程等进行了模拟，结果表明模型计算的下泄洪水过程（水位、流量）及泄流渠断面发展等与实测情况符合较好。

关键词：渍坝；洪水；泄流；冲刷；模型；唐家山堰塞湖中图分类号：ｐ３１５．９文献标识码：Ａ１唐家山堰塞湖形成及除险概况２００８年５月１２１３四川汶川发生特大地震，距＝ｊＥ川县城约４ｋｍ的涪江支流通口河右岸唐家山发生大规模滑坡，滑坡堵塞通口河形成堰塞坝，通口河水壅高形成唐家山堰塞湖。

通口河系涪江一级支流，主河长１７３ｋｍ，控制集水面积４５２０ｋｍ２。

，流域内多年平均年降水量为１３５５．４ｍｍ，降雨主要集中在５～９月，占全年降水的８６．３％，历年１ｄ最大降水量为３２３．４ｍｍ。

，多年平均流量１１８ｍ３／ｓ，年径流量３７．１亿ｍ。

径流年内分配不均，其中７～９月多年平均流量２４９ｍ３／ｓ，相应径流量１９．８亿ｍ３，占全年径流量的５３．５％。

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

数学建模竞赛题目唐家山堰塞湖泄洪问题研究摘要：本文研究了唐家山堰塞湖及其泄洪规律，并对其蓄水、溃坝、泛洪等三个过程进行了模型分析。

问题一，是堰塞湖的蓄水过程。

本文采取了两种方法得到蓄水量和水位高程关系曲线，一是采用曲线拟合的方法；二是通过数字高程模型（DEM），使用连通性分析算法。

对比发现，前者在709至742米之间的精度较高，而后者普适性好。

而后，本文建立调洪演算模型得到了降雨量与水位高程的关系曲线，并得到不同降雨量影响水位高程达到泄洪槽水位和警戒水位的具体时间。

问题二，是堰塞湖的溃坝过程。

通过对坝体物质组成及其上下游地理情况的分析，可以得到溃坝的方式应属于逐渐溃决。

本文建立了由溃口冲刷扩展速率，溃口动态出流计算公式联合求解的溃口流量非恒定过程计算模型，并进行了灵敏度分析。

得出结论：溃口在1小时达到稳定且溃口最大流量约为140003/m s。

问题三，是堰塞湖的泛洪过程。

本文通过种子蔓延算法求得淹没区域，其总面积约为1412.64万平方米，其中绵阳市的淹没面积为481.10万平方米，人员撤离18万左右。

考虑到圣维南方程组过于复杂，且求解困难，本文采用李斯特万公式和黄河水利委员会水利科学研究经验公式模拟了洪水演进过程。

结论指出洪水的影响地区有北川县城、通口镇、含增镇、青莲镇、绵阳县城，到达时间分别是0.121h，2.125h，3.047h，5.241h，7.856h。

问题四，本文讨论了降低唐家山堰塞湖水位的方案，并提出了使用虹吸管的建议；进行了导流槽效用分析，提出如何以最合理的方式构建导流槽；提出了将唐家山堰塞湖改造为大型水库的想法；对苦竹坝下游的堰塞湖的形势进行了分析，建议在上游未来水之前，及时对苦竹坝下游的堰塞湖采取人工开挖或爆破拆除，以免危害北川县城。

参赛队号随机A8目录目录 (2)一、问题重述 (3)二、问题分析 (3)三、模型假设 (4)四、问题1模型建立与求解 (4)4.1水位高程和蓄水量 (4)4.1.1曲线拟合模型 (4)4.1.2基于DEM的蓄水量计算 (6)4.1.3连通性分析算法 (7)4.1.4唐家山堰塞湖蓄水量计算 (10)4.2降雨量和水位高程 (11)4.2.1调洪演算模型 (11)4.2.2 堰塞湖的蓄水特征 (12)4.2.3由多年平均降雨推求水位高程 (13)4.2.4由天气预报降雨推求水位高程 (14)五、问题2模型建立与求解 (17)5.1溃坝模型 (17)5.1.1溃口流量 (17)5.1.2溃口参数选择 (18)5.2 对唐家山堰塞湖的分析及假设 (19)5.2.1溃坝方式分析 (19)5.2.2残留坝高的确定 (19)5.2.3溃口冲刷速率分析 (20)5.3溃口冲刷稳定时间的确定 (20)5.4溃坝模型的灵敏度分析 (21)5.5进一步讨论 (23)六问题3模型建立与求解 (23)6.1洪水淹没区域的计算 (23)6.1.1淹没区域计算模型 (23)6.1.2数据计算 (25)6.1.3淹没面积以及人员撤离方案 (26)6.2洪水演进计算 (28)七、问题4模型建立与求解 (29)7.1 降低堰塞湖水位方案讨论 (29)7.2 关于导流槽方案效用分析 (30)7.3、唐家山堰塞湖处置分析 (31)7.4、应对唐家山下游的堰塞湖采取应急处理措施 (32)参考文献 (33)一、问题重述汶川地震给我国的给人民生命财产和国民经济造成了极大的损失。

浅谈唐家山堰塞湖抢险安全管理

浅谈唐家山堰塞湖应急排险的安全管理摘要：在方圆不足一平方公里的堰塞体上，在人员数量多、机械设备多、安全隐患多、施工强度大、条件艰苦的情况下，创造了施工作业零伤亡的良好业绩。

本文从安全制度建立、安全作业监督管控、安全观测与警戒、应急预案建立与完善等方面进行了阐述。

关键词：排险安全管理1 前言唐家山堰塞湖位于四川省北川县上游通口河上，是因“5.12”汶川大地震造成了龙门山300KM的地震断裂带，唐家山山体大面积滑坡，封闭河道，原有水系被堵塞物堵住。

河谷、河床被堵塞后，水位迅速积蓄抬高，严重威胁下游。

由于排险时间要求紧、施工强度高、作业面狭窄、人员和机械设备多、相互干扰大，所以本堰塞湖应急排险工程的特点决定了安全管理的重要性和难度性，稍有不慎，便会发生群死群伤的特大安全事故。

其安全隐患分析情况如下：1.1堰塞体本身是由山体滑坡后形成的堆积体,结构松散,加之湖水不断上涨,压力增大,易形成较大的渗漏通道，造成管涌等甚至坝体溃决的危险。

1.2由于地震,堰塞体两边山坡已形成较大裂缝，而所有排险人员驻地就设在堰塞体上，施工也在峡谷中的堰塞体上，若发生新的滑坡塌方，势必造成特大安全事故。

1.3根据排险方案，施工采取以“开挖渠道引流、分段同时实施”为主、铅丝石笼护坡为辅的方式，作业面狭小、施工人员多、相互干扰大，易发生渠道边坡垮塌掩埋、人员砸伤等安全事故。

1.4气候条件差，高温炎热，能使用的水源缺乏，卫生条件恶劣，易发生职业健康卫生事故。

2 建立安全制度2.1安全组织机构与职责在唐家山堰塞湖应急排险前线指挥部的领导下，建立了专门的安全管理机构，专人专责。

在现场从各施工单位抽调安全知识丰富、责任心强的干部战士组成专业安全分队负责所有安全管理工作，主要有渠道开挖施工安全分队，主要负责渠道开挖施工作业安全；左岸山坡裂缝发展情况监测分队，负责观察边坡裂缝变化情况；右岸山坡滑坡监测分队，主要对山体滑坡、塌方等情况发出预警预报，防止发生塌方对渠道施工人员造成伤害；上游来水量监测分队，主要监测湖水上涨情况；下游渗漏变化情况监测分队，主要对堰塞体下游渗水点、出水量、浑浊度等情况进行巡视。

唐家山堰塞湖可能起溃水位分析研究

第39卷第22期2008年11月　　人　民　长　江Y angtze River 　Vol.39,No.22Nov.,2008收稿日期作者简介李书飞,男,长江水利委员会设计院规划处,工程师,硕士。

文章编号:1001-4179(2008)22-0073-03唐家山堰塞湖可能起溃水位分析研究李书飞　胡维忠(长江水利委员会设计院,湖北武汉430010)摘要:确定起溃水位在溃坝洪水分析计算中十分重要。

在唐家山堰塞湖可能起溃水位分析研究中,采用堰流公式法、渠道水面线法和MI KE 11模型等3种方法进行了泄流渠过流能力计算,分析计算可能起溃水位的影响。

在此基础上,分析了开挖泄流渠、扩宽泄流渠渠底以及降低泄流渠渠底高程对降低可能起溃水位的效果。

计算分析表明从符合实际情况和偏安全的角度出发,在可能起溃水位分析采用渠道水面线法较合适。

降低泄流渠的底高程扩宽渠底与宽度比,更能降低起溃水位。

在唐家山堰塞湖开挖泄流方案制定中,应尽可能降低泄流渠渠底高程,从而降低起溃水位和可泄水量,减轻溃坝洪水对下游人民生命财产的威胁。

关　键　词:溃坝洪水;溃坝洪峰流量;可能起溃水位;渠道水面线法;唐家山堰塞湖中图分类号:P315.9 文献标识码:A1　概述唐家山滑坡形成的原因是强烈度地震的作用,其基岩主要堆积于坝体左岸及右岸下部,上部坡残积碎石土层主要堆积右岸上部及堰塞坝坝前。

堰塞坝最可能的破坏形式是在大流量水流冲刷作用下的滑塌破坏。

在此种情况下,分析溃坝洪水对下游人民生命财产的威胁、提供安全避险策略以及采取何种工程措施降低溃坝洪灾损失至关重要。

溃坝洪水分析的主要内容之一是坝址溃坝流量及其过程线。

坝址溃坝洪峰流量取决于起溃水位、溃坝历时、溃口形状及其发展过程等因素。

起溃水位与堰塞坝坝体物质组成和结构特征密切相关,这里仅从洪水调节的角度界定可能起溃水位的概念,以下的分析也是基于这一基本概念。

可能起溃水位是指在坝址发生一定频率洪水,库水位上升至最高溃坝情况下的水位。

汶川地震中唐家山堰塞湖泄洪问题

第18组刘慧妍何栓张丰宇汶川地震中唐家山堰塞湖泄洪问题摘要2008年5月12日14：28在我国四川汶川地区发生了8.0级强烈地震，给人民生命财产和国民经济造成了极大的损失。

地震引发的次生灾害也相当严重，特别是地震的造地运动形成了三十多个高悬于灾区人民头上的堰塞湖，对下游人民的生命财产和国家建设构成巨大的威胁，其中以唐家山堰塞湖尤为严重。

本文就汶川地震中唐家山堰塞湖泄洪问题展开了讨论针对问题一：本文首先利用DEM图确定堰塞湖不同高程下的湖面面积，然后采用微元的思想对堰塞湖的体积进行积分，以坝前水位为自变量，总蓄水量为因变量建立了唐家山堰塞湖的水位高程与总蓄水量的模型Ⅰ，再在模型Ⅰ的基础上，通过分析堰塞湖每天的新增水量来源，将新增水量与降雨量之间的关系转化为坝前水位与降雨量的关系，之后，根据实际降雨量，进行数据拟合，最终建立降雨模型Ⅱ。

针对问题二：文中采用了溃坝洪水预报DAMBRK模型。

建立了溃坝模型Ⅲ。

首先，利用为微元法推导出溃坝模式下的流量公式，并采用迭代法求出坝址处流量过程线和坝前水位变化趋势。

然后，利用迭代法求解圣维南方程组，得到洪水在河道中的演进过程。

即洪峰沿河道到达下游各城镇的时间、流量和水位。

针对问题三：首先，根据题目要求，本文建立了堰塞湖1/3溃决下，洪水的水流速度、淹没区域以及到达时间等模型，即三分之一溃决计算模型Ⅳ。

利用DEM图和MATLAB软件模拟了汶川地震堰塞湖1/3溃决后的洪水演进情况，然后，设计了“人员撤离紧迫性”的量化方法，并在此基础上建立了人员撤离模型Ⅴ，拟定了下游人员撤离方案，为政府相关部门提供了数据支持，最后对模型的精确性进行了评价。

针对问题四：根据本文所建立的水位高程与总蓄水量的数学模型Ⅰ，降雨模型Ⅱ，溃坝模型Ⅲ，三分之一溃决计算模型Ⅳ和人员撤离模型Ⅴ。

在前对前三个问题进行求解的基础之上，分析当时所采取对策的正确性和改进的可能性，并且对地震后次生山地灾害展开了讨论，如：对于山地灾害而言，降雨是重要的诱发因素。

唐家山堰塞湖泄洪对通口水电站设计的全面检验及通口水库的滞洪作用

唐家山堰塞湖泄洪对通口水电站设计的全面检验及通口水库的滞洪作用杨光伟;李西瑶【摘要】2008年5月12甘四川汶川发生8.0级地震,导致通口河口上游山体垮塌,在庸家山形成特大型堰塞湖,交通中断.建于典型岩溶地段的通口水电站大坝在堰塞湖下游约23km处,为碾压混凝土高坝,泄洪消能系统采用溢流表孔+宽尾墩+戽式消力池的方式.工程建在高烈度地区,距离震中较近,在汶川特大地震中震损轻微,防渗体系完好,在唐家山堰塞湖泄洪过程中的滞洪削峰为保护下游沿岸和城市发挥了很大的作用.本文对此进行了介绍.【期刊名称】《水电站设计》【年(卷),期】2010(026)001【总页数】7页(P15-21)【关键词】水电站;岩溶;枢纽布置;泄流消能;基础防渗;地震;滞洪削峰;通口水电站【作者】杨光伟;李西瑶【作者单位】中国水电顾问集团成都勘测设计研究院,四川,成都,610072;中国水电顾问集团成都勘测设计研究院,四川,成都,610072【正文语种】中文【中图分类】TV7421 工程概况1.1 工程简介2008年 5月 12日 14时 28分四川汶川发生8.0级地震,导致通口河中上游山体垮塌,在唐家山形成特大型堰塞湖,交通中断。

建于典型的岩溶地段的通口水电站大坝在堰塞湖下游约 23km处,电站位于涪江右岸一级支流通口河中游末段的北川羌族自治县通口镇上游 2km,左岸属江油市,右岸属北川县,与北川县城及江油市均有公路相通。

坝、厂区河道平均比降为5‰,坝址多年平均流量为 93.4m3/s。

枢纽工程主要由拦河坝、右岸引水系统、地面主(副)厂房和升压站等组成,电站主要任务是发电。

工程等别为Ⅲ等,主要永久性建筑物按 3级设计,次要建筑物按 4级设计。

水库正常蓄水位 598.00m,总库容 3610万 m3,调节库容 97万m3,具日调节性能。

电站设计水头 48.0m,引用流量 106.504m3/s,电站装机 45MW,保证出力9.60MW,多年平均发电量 2.31亿kW·h,年利用小时数 5133h。

唐家山堰塞湖成功排险的探索与思考

难重重，５月２５日下午进行首次方案调整，调整为
高程６７５ｎ，堰体最大高度１４４ｎ，垭口处堰高６．ｌ２．ｌ８．５ｎ。总库容３１亿ｎ，回水长度约２ｍ。若２６ｌ．５ｌ３ｋ发生溃坝将影响下游北川县，江油市，绵阳市的涪城
堰体长８３４ｎ，宽６１８ｎ，面积３．２万ｎ，体０．ｌ１．ｌ０７ｌ
积２０７万１，堰体最高点高程７１９１，堰体底部３ＴＩ９．ＩＴ
因自５月２日至２日，唐家山堰塞湖上空气３５象条件极差，能见度非常低，直升飞机吊运设备困
区、科学城、农科区以及三台县人口１０多万，甚至３影响到遂宁市，属 “ ．２汶川特大地震中形成的极５１” 高危级堰塞湖。
机械施工与人工爆破 “ 管齐下 ” 工爆破为主的施双人工方案。后因机械施工情况良好，爆破施工未实
泄流渠进口高程７５０ｎ、渠底宽２ｌ４．ｌ２ｎ，泄流漕总
理；立足早处理，防止因水量增大增加处理难度；立
即实施工程排险方案并确保施工人员安全的原则进行
施工规划。根据设计方案，按照 “ 员轮班作业，设人备持续作业 ” 的原则进行方案实施。
３施工组织
３１组织机构．
长６０ｎ，上游平缓段纵坡０６，下游陡坡段纵８ｌ．％坡分别为２％和１％；高标准方案，泄流槽进口４６高程７２０ｎ、渠底宽１ｌ４．ｌ３ｎ，泄流槽总长６５ｎ，９ｌ上游平缓段纵坡０６，下游陡坡段纵坡分别为．％

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

数学建模竞赛题目唐家山堰塞湖泄洪问题研究摘要：本文以新闻报道中的数据和数字高程图中的地貌信息等为依据，研究堰塞湖的泄洪问题。

对唐家山堰塞湖未来18天的水位高程进行了预测，模拟了堰塞湖的溃坝，并对洪水的蔓延进行了分析，计算了1/3溃坝时下游地区的淹没情况。

第一问，将湖体简化成横截面为梯形的几何体，建立水面高程与蓄水量的关系模型，根据数据，回归分析确定参数。

并对该模型进行了改进，以堰塞湖湖底的等高线为底，等高线之间高度差为高，借用积分思想，将湖体看成由底面积不同的柱体叠加而成，建立水位高程与蓄水量的关系模型。

将堰塞湖的集雨地区网格化，运用附近县城的降雨量插值计算每个网格点的降雨量。

引入延迟因子、入渗因子、蒸发因子，根据预报降雨量的50%，80%，100%，150%计算每日堰塞湖的蓄水量，通过水位高程与蓄水量的关系，求出堰塞湖水位每日上升的高度。

第二问，把唐家山堰塞湖泄洪过程的数据作为研究溃坝过程的重要依据。

设定溃口形状参数，建立溃口形状及洪水流量在溃坝过程中的连续性方程组。

为了计算机求解的方便，转换成差分方程模型，迭代求解。

并对差分方程模型的结果进行了分析，对模型的合理性进行了评价。

第三问，把洪水分成很多具有一定体积的“分子”，应用面向对象的思想，规定“分子”的属性（体积和能量）和动作（跃迁）建立“分子”扩散模型。

把地表网格化，取每个格子所在区域的平均高度作为该格子的高度。

洪水下泻过程即为分子在网格上的跃迁过程。

分子每跃迁一步都会使自身能量变化，如果跃迁后能量小于0，则该分子结束它的生命周期。

通过记录“分子”到达的各格子的时间等数据，确定洪水到达各地区的时间，淹没范围及水量。

并给出淹没地区的危险等级，制定人员转移安排方案，见7.3.2。

第四问，根据上述模型，分析了政府采取政策的合理性。

建立最佳泄洪槽深度的模型，对其进行讨论。

并对地质灾害后的次生山地灾害，给出了预防建议。

最后对模型的不足提出了改进方向和建议。

简析唐家山堰塞体溃口宽度的模拟计算论文

简析唐家山堰塞体溃口宽度的模拟计算论文简析唐家山堰塞体溃口宽度的模拟计算论文唐家山堰塞体为地震诱发山体滑坡而堆积形成的天然坝,相对堆石坝具有以下特点:泥沙颗粒不均匀,组成物质松散,渗流严重,岸坡不稳定,坝体稳定性更低,其力学稳定性、抗冲性较差.受二次滑坡的影响,坝体呈现明显的双层结构,上层为含水性较好的粘性土层,下层为较易崩塌的砂土层.唐家山溃口的形成发展过程也较为复杂,为人工开挖、溯源冲刷和溃口崩决的混合过程,在溃决过程中存在一定的漫流,同时双层土体结构导致水头差较大,渗流明显.匡尚富提出,天然坝的溃决过程可以分为溢流侵蚀溃决型、滑动崩决型、逆流渐进溃决型等.对比唐家山堰塞体的实际情况,其侵蚀方式为溢流侵蚀,而堰塞体则与多发生滑动崩决型溃决的坝体接近.1唐家山堰塞湖溃口唐家山堰塞湖库水位持续上升,严重威胁水库下游人民群众的生命和财产安全.因此在坝顶位置沿河开挖了溢洪道宣泄洪水以减轻溃坝洪水灾害.溢洪道底部尺寸为695m×4m,两侧边坡1∶1.5,入口处高程为742m.由于底部及边坡均由碎石堆砌而成,溢洪道可蚀性强.2008年6月7日7:06时,溢洪道开始排水,入口处的初始水位为740.37m,初始水面宽为58m,初始水深0.18m.到6月10日1∶30时,水位攀升至最大值743.10m,随后逐渐降落,在6月10日8:00时降落至742.18m.该水位持续了1h,随后迅速降落,从9:00时的742.17m降落至11:30时的735.81m,对应的最大洪峰流量为6900m3/s.随后,在水流的持续冲刷及重力作用下,溢洪道迅速垮塌.在6月10日15:00时,溢洪道达到临时稳定状态,此时水面宽145m,河底高程721.2m.最后相对稳定的溃口宽为131.2m,溃口深为20.8m.2计算结果及分析为了寻求合适的模型以模拟唐家山堰塞体的溃决,采用以Osman 公式为代表的河岸展宽模型及以Froehlich方法和垦务局模型为代表的土石坝溃决模型,比较分析其应用于计算唐家山溃口模拟的结果.2.1河岸展宽模型Osman和Thorne提出计算河岸横向冲刷距离的公式,在Δt(s)时间内,粘性土河岸被水流横向冲刷后退的`距离为ΔB=Cl×Δt×τf-τce-1.3τcγbk(1)式中,ΔB为Δt(s)时间内河岸因水流侧向冲刷而后退的距离,m;γbk 为河岸土体容重,kN/m3;τf为作用在河岸上的水流切应力,N/m2;τc为河岸土体的启动切应力,N/m2;Cl为河岸抗冲系数.2.2土石坝溃坝模型在Fu等中,采用Macchione提出的预测土坝溃决洪水的方法模拟预报了唐家山溃口的洪水过程,但是Fu等的模拟计算表明,需要采用不同的特征速度ve来分别模拟洪峰的峰现时间和洪峰流量,特征流速ve 分别取值0.33m/s和1.9m/s可以模拟得到洪峰流量和洪峰流量的发生时间,而Macchione特征流速范围为0.05—0.10m/s.即Fu等的模拟中特征速度ve取值0.33m/s才能较为合理地模拟唐家山堰塞湖溃口的洪峰流量,相比原文献中的区间范围取值偏大约3倍.因此,本节不再对比分析Macchione方法的适用性.对比Froehlich和Macchione公式的适用范围可发现:Macchione 为土石坝的溃坝模拟,用于验证Macchione模型的溃坝实例多为土坝,仅少量坝体为细沙、土和卵石的混合物,而溃决原因多为管涌,少部分为漫流和渗流导致的溃堤,此外,Macchione假定溃决断面为三角形,在溃口的发展过程中与实际相差较大.而Froehlich和垦物局的方法则是较多地依赖水深,与坝体结构等关系不大,与实际情况不符,导致计算结果与实测不符.3结论从计算与实测结果的对比可以发现,Osman方法、Froehlich和垦务局的模拟结果均不是很理想,不能满足唐家山堰塞体的溃口宽度计算,需要针对堰塞体的具体特点改进模型.由于Osman公式为模拟非粘性河岸展宽的水动力学-土力学方法,是基于物理机理的河岸坍塌模型.接下来仅讨论Osman公式中尚需要改进的地方.对于将Osman公式应用于唐家山堰塞体,除夏军强等指出的两处不足(不能考虑孔隙水压力和静水压力对河岸稳定性的影响,和假定河岸崩塌时的平面滑动面必须限于通过坡脚),本文认为还存在以下困难:一是不能考虑土体具有明显的双层结构情况下的大水压差;二是导流渠弯道中二次环流的影响难以估计.有学者针对各项不足提出改进方法.如Darby和Thorne针对不考虑孔隙水压力和静水压力对稳定性的影响及滑动面必须通过坡脚的假定这两个不足改进了Osman模型.在模拟顺治河段的河床变形中,Darby和Thorne考虑河岸可能发生平面滑动和圆弧滑动两种破坏情况,夏军强等指出,虽然计算值与实测值仍存在一定的误差,但同时考虑两种破坏情况更符合实际.弯道对河岸展宽的作用则可根据Thorne和Osman中的方法得到改进,弯道剪切力与平均剪切力的关系式为τbτ=frcW,其中,τb为弯道外测的剪切力,τ为平均剪切力,τ=rwRS,rc为弯道曲率半径,W为河宽,将τb代替式(1)中的τf即可考虑弯道作用,该关系根据文献中的图即可获.但是,双层土体导致的大水压差仍需要深入的研究.。

5A-唐家山堰塞湖泄洪问题研究(1)

第39卷第16期2009年8月数学的实践与认识M A TH EM A T I CS I N PRA CT I CE AND TH EO R Y V o l 139　N o 116　A ugu st,2009　唐家山堰塞湖泄洪问题研究曾向荣,　郝红星,　孙博良(国防科技大学信息系统与管理学院,湖南长沙　410073)摘要:　研究了唐家山堰塞湖的泄洪规律,并对其蓄水、溃坝、泛洪等三个过程进行了模型分析.关键词:　唐家山;堰塞湖;泄洪;蓄水;溃坝;泛洪1　问题描述和分析收稿日期:2009204230基金项目:国防科技大学数学建模创新实践基地(S M 0802) 地震后形成的唐家山堰塞湖对下游人民的生命财产和国家建设构成了巨大的威胁.因此对唐家山堰塞湖水情的关注是非常必要的.分析堰塞湖的各种水情,采取相应的对策,最大程度的减少人民生命财产和国民经济的损失.为了研究唐家山堰塞湖的泄洪规律,需要分析堰塞湖三个位置、三个阶段的水情,如图1所示,三个位置分别是坝前、坝体、坝后,三个阶段则为蓄水阶段、溃坝阶段、泛洪阶段.图1　堰塞湖三个位置和三个阶段示意图对于蓄水阶段,可以通过数字高程模型(D E M ),使用连通性分析算法得到蓄水量和水位高程关系曲线;又由于存在蓄水量和水位高程变化的历史资料,也可以采用曲线拟合的方法.降水量对水位的影响则可以使用调洪演算模型得出.对于溃坝阶段,溃口流量过程可以使用溃口冲刷扩展速率,溃口动态出流计算公式联合求解的溃口流量非恒定过程计算模型进行分析.对于泛洪阶段,基于D E M ,可以采用搜索算法进行淹没面积的计算.而洪水到达时间等值的计算则需要在经验公式和基于圣维南方程的数值模拟之间进行选择.2　蓄水阶段的模型建立与求解2.1　曲线拟合模型根据地震期间新闻提供的数据[127],以水位高程为自变量,蓄水量为纵坐标得到分布图2所示.分析图中点的分布特征可以得到蓄水量与水位高程的关系并不是简单的线性关系,而根据实际情况容易得到,由于水位高程的增大会使蓄水面积增大,因此随着水位高程的增大,需水量的增大速度应当逐步增大.考虑到散点的分布特征,本论文引入二次曲线、指数函数和双曲线等几种预测模型,通过已知点用M A TLAB 对上面三种模型进行拟合,可以得到各模型参数以及误差如表1所示.表1　曲线拟合参数计算结果模型参数值误差二次曲线模型(y =ax 2+bx +c (a ≠0))a =6.2945b =-8608c =29459001203.7指数函数模型(y =y 0e Αt +b )y 0=658Α=-3b =1754831417双曲线模型y =a +bc +t a =59100b =-3362400c =-6009367从表1可以得出,运用二次曲线进行拟合得到的表达式的拟合误差比较小,因此本文选择水位高程h (m )与蓄水量x (104m 3)的模型为:h =6.2945x 2-8608x +2.9459×106(x >683)拟合曲线如图3所示.图2　水位高程与蓄水量数据点示意图图3　水位高程与蓄水量拟合曲线从图3可以看出,用二次曲线对水位高程为709～742米时的蓄水量拟合程度比较好,但是二次曲线只能拟合当水位高程大于683米的情况,而且当水位高程很大时,该曲线得到的蓄水量将会和实际差别很大,因为该曲线不能反映实际的地形条件.2.2　基于D E M 的蓄水量计算[829]如果考虑更精确地计算唐家山堰塞湖的蓄水量,那么就应该立足于唐家山堰塞湖的具体地形,因此,考虑基于D E M 来计算.本文将数字地图用规则网格覆盖,然后进行网格抽取,由于规则网格将区域空间切为规则的格网单元,而每个单元网格的高程本文用该网格中心点的高程作为该网格内的平均高程,这样就可以得到该地图的高程矩阵.规则格网的高程矩阵,很容易进行计算机处理,例如可以进行等高线计算,坡度坡向以及自动提取流域地形等.如果数据量不大,可以运用手工83数　学　的　实　践　与　认　识39卷填表的方法采集的网格化高程数据,但是假如数据量过大,有必要研究D E M 文件的存储格式,进行程序化处理.在后面的洪水演进中,根据网格化高程数据计算地形的连通性是非常重要的,此部分将进行连通性分析.当选定一个边缘网格后,本文按照低于此网格高程数据h 的标准寻找其周围的网格并进行标注.当选定一个网格后,下一步要对其相邻的其他三个网格的高程进行判断,连通性分析如图4所示.连通算法的表述为:选定一个需要计算联通的网格,其高程值记为h ,将其高程值与周围网格高程值进行比较,按照低于h 的标准判断其可能的方向,并对可能的方向数量进行累加,并将其判断过的网格进行标注,然后取其中一个方向继续判断,把方向数量减1,循环判断.但是在后来的判断过程中可能遇到方向指示已经标注过格网,这样为了不重复标记,对于此方向不进行计数.每次累加方向后判断方向数量是否为零,不为零时继续判断,直到方向数量为零的时候,连通范围也就判断出.整体算法流程图可以表述如图5.在该流程图中,m 表示该网格周围被判断过的网格的数目,流程图比较重要的是:1)判断该网格与周围的连通关系;2)判断周围网格是否标记;3)继续进行连通性判断.图4　连通性分析算法示意图图5　连通性分析算法流程图2.3　唐家山堰塞湖蓄水量计算在提取连通性网格的时候,同时对于网格数量记数为p ,对于提取网格表示为H (i ,j ),这样堰塞湖湖面面积、蓄水量按如下公式进行计算S y =p ×d x2V y =d x 2×∑h (i ,j )其中S y 表示堰塞湖湖面面积,d x 表示正方形格网边长,V y 表示堰塞湖蓄水量,h (i ,j )表示格网高程数值.对于格网淹没的高度按如下公式计算H y =h -h (i ,j )其中H y 表示格网水面高度,h (i ,j )表示格网水位高程数值.D E M 计算得到的蓄水量与水9316期曾向荣,等:唐家山堰塞湖泄洪问题研究位高程曲线与拟合曲线对比如图6.图6　拟合曲线与D E M 计算曲线对比图前面已经提到,随着水位高程的增加,每增加单位高程的蓄水量的增加逐步变大,而从上图中的两条曲线分析,曲线的斜率逐渐变大,符合实际情况.另一方面,由于二次曲线拟合是一条抛物线,而抛物线的拐点在水位高程为683米处,由于随着水位高程的降低,蓄水量增加显然不合理,因此在抛物线对称轴左侧不适宜采用拟合公式,而运用调洪演算得到的曲线较好地克服了这种拟合的缺点.虽然拟合曲线在媒体报道点附近比调洪演算曲线更好地符合了实际数据值(3亿立方米),但是曲线拟合并不能很好地反应拟合点界限以外的数据.而运用调洪演算得到的曲线由于根据出入水量的关系进行演算,能够在较大的范围内很好地反应水位高程与蓄水量的关系.2.4　调洪演算模型1)水量平衡方程在堰塞湖的蓄水过程中遵循水量平衡,即在某一时段内(∃t =t n -t n -1),入湖流量与泄流量之差应等于该时段内堰塞湖蓄水量的变化值,由于蒸发水量很少,因此忽略堰塞湖蒸发水量,表示为:Q n +Q n -12∃t -q n +q n -12∃t =W n -W n -1(1)式中:Q n -1,q n -1分别为时段初入、出湖流量;Q n ,q n 分别为时段末入、出湖流量;Wn -1,W n 分别为时段初、末湖蓄水量;∃t 为时段长.由于不考虑泄洪,因此q n ,q n -1,W n =0,得:Q n +Q n -12∃t =W n -W n -1(2) 2)数值解析法调洪计算原理数值解析法采用堰塞湖水量平衡的微分方程.得:d h d t =Q (t )F (h )h (t 0)=h 0(起调水位)(3) 上式中,F (h )为堰塞湖水面面积,湖水面的面积由上面讨论的D E M 算法求得;Z (t )为水位高程;Q (t )为入湖流量;h 0,t 0分别为起调水位和起调时间.04数　学　的　实　践　与　认　识39卷图7　降水汇入堰塞湖流图2.5　堰塞湖的蓄水特征由于要考虑降雨量与水位高程的关系,入湖流量为堰塞湖上游流量、区间入流和直接降到湖面的水量,后两者都与降雨量有关,如图7所示.因此,得到Q =f (R )+M(4) 式中:Q 为入湖流量,R 为降雨量(mm ),f (R )为关于R 的函数,M 为堰塞湖上游流量.则第t 日的入湖流量为Q (t )=f (R (t ))+Md t(5) 假设相邻两日的堰塞湖上游流量相等,得Q (t )=d f (R (t ))d t(6) 联立(3)和(6)即可推算降雨和水位高程的关系2.6　由多年平均降雨推求水位高程如果要在5月25日前对接下来17天唐家山堰塞湖水位高程进行预测,那就只能依赖历年的这些天的平均降雨情况了,因为天气预报很难准确预测17天的降雨情况.根据历史日降雨量系列,分析计算多年平均情况下5月25日至6月12日降水过程,按起调水位723m (5月25日8点实测值),进行唐家山堰塞湖未来堰前水位预估.多年平均来水情况下堰塞湖堰前水位过程如图8.图8　多年平均来水情况下唐家山堰塞湖堰前水位变化过程根据上述分析计算结果可知,若唐家山堰塞湖以上流域未来一段时间降雨按多年平均考虑,到6月4日,唐家山堰塞湖蓄水将超过泄洪槽740m ,如果没挖泄洪槽,则6月11日后将会发生漫顶.2.7　由天气预报降雨推求水位高程根据天气预报结果,按其预报的上限、中间值和下限整理,得到唐家山堰塞湖以上流域5月25日至6月12日的降雨预报成果,根据多年平均日降水过程,按起调水位723m (5月25日8点实测值),进行唐家山堰塞湖未来堰前水位预估多年平均来水情况下堰塞湖堰前水位过程见图9.1416期曾向荣,等:唐家山堰塞湖泄洪问题研究图9　不同预报来水情况下堰塞湖堰前水位变化过程对比图根据上述分析结果,若唐家山堰塞湖以上流域未来一段时间降雨按天气预报结果考虑,到6月5日,唐家山堰塞湖蓄水将超过泄洪槽740m ,如果没挖泄洪槽,则6月12日后将会发生漫顶.预测当降雨量为天气预报降雨量中值的50%、80%和150%时,水位的结果如图10:图10　水位高程预测曲线从上面图表可以看出,随着降雨量的增加,水位高程增加的越快,天气预报降雨量的50%、80%、100%和150%时,堰塞湖湖水超过泄洪槽的时间依次为:6月9日、6月7日、6月5日和6月5日.当降雨量为天气预报降雨量的150%时,6月11日将发生漫顶,而为天气预报降雨量的50%、80%、100%时,6月12日以后水位高程达到警戒水位.3　溃坝阶段的模型建立与求解3.1　溃口流量溃口流量利用溃口动态出流公式[10]计算Q b =c v k s [3.1b i (h -h b )1.5+2.45z (h -h b )2.5]式中,c v 为对行近流速的修正,b i 为溃口瞬时底宽(m ),h 为堰塞湖水位高程(m ),h b 为溃口底高程(m ),z 为溃口边坡(1:z ),k s 为尾水影响出流的淹没修正.24数　学　的　实　践　与　认　识39卷当(h l -h b ) (h -h b )>0.67时,k s =1.0-27.8h l -h b h -h b-0.673,否则,k s =110.上式中h l 为尾水位高程(m ).c v 由下式计算c v =1.0+0.023Q 2b [B 2d(h -h b m )2(h -h b )]式中B d 为堰塞湖坝前宽度(m ),h b m 为最终溃口底高程(m ).如果溃口因管涌引起,那么可以应用孔流方程.Q b 可用以下孔流方程替代:Q b =4.8A Θ(h -h θ)1 2其中A Θ=2b i (h p -h b )式中h p 为管涌中心高程线(m ),h θ=h p ,如果尾水位高程h l >h p ,则h θ=h l .当坝前水位下降或管涌扩大到满足如下条件h <3h p -2h b则溃口流量由孔流转变为宽顶堰流.3.2　溃口参数选择早期由于对溃坝机理不甚了解,一般假设溃坝形式为瞬时全溃或顺势部分溃.对混凝土拱坝可以适用,但对土石坝和混凝土重力坝则不太适用.溃口形成如图11所示图11　溃口形成示意图[10]对逐渐溃来说,主要是确定两个参数:最终溃口底宽b 和形状z (即溃口边坡).最终溃口底宽b 与溃口平均宽度b λ关系如下b =b λ-z h d式中h d 为溃口底部以上水高,一般以坝高近似代替.假设溃口底宽从一点开始,然后在溃决时间范围Σ内逐渐以线性或非线性速度增长,直到形成最终的溃口底宽b 和最终溃口底高程h b m .如果溃决时间Σ小于一分钟,那溃口底宽开始于限定的b 值,而不是某一点.与冲刷溃决相比,这更像是崩塌溃决.溃口底高程h b 是时间Σ的函数h b =h d -(h d -h b m )t b ΣΘ0,　0<t b <Σ 最终溃口底高程h b m 通常是水库或泄水河道段底部高程,但不是必须的.t b 为溃口形成开始算起的时间.Θ0为非线性程度参数,1ΦΘ0Φ4.通常假设为线性速度增长,因而Θ0=1.溃口瞬时底宽则由下式确定b i =b t b ΣΘ0,　0<t b <Σ 溃坝模拟时,当库水位高程h 超过某一指定高程h f 时,溃口开始形成.对于土石坝,溃3416期曾向荣,等:唐家山堰塞湖泄洪问题研究口平均宽度b λ为0.5h Φb λΦ8h d ,如果是管涌引起的溃决,模拟管涌初始溃口Θ0Ε2是合适的.3.3　溃口冲刷速率分析溃口冲刷速率假定溃口冲刷形成的总时间为T ,总的溃坝深度为Z ,按照下述方法确定,可得到任意时间内的溃口发展数度V d .由以上分析,本文可以假定溃口发展模式符合以下曲线(如图12).由溃口发展速度基本模式可以得到,当溃口发展时间t ΦT 2时,冲刷速率V d =2(V m ax T )×t 当溃口发展时间T 2Φt ΦT 时,冲刷速率V d =-2(V m ax T )×t +2V m ax 当溃口发展时间t >T ,冲刷速率为Vd =0,可见V m ax =2Z m ax T .上述各式中,V m ax 为溃口冲刷过程中的最大冲刷速率,t 为是时间,Z m ax 为给定可能冲刷溃口高度.3.4　溃口冲刷稳定时间的确定就目前的研究水平,对溃口流量与冲刷速度的关系仍然不是十分清楚,存在较大的不确定性,很难定量描述.对于如何确定溃口冲刷稳定时间,文章采用与宽顶堰公式和肖克列奇公式进行比较的方式进行确定.将溃口冲刷速率公式带入溃口动态出流公式,在残留坝底为720m 的条件下,本文可以得到在溃口形成时间变化时,溃口流量的变化曲线,如图13.图12　溃口形成速率示意图图13　溃口流量随时间变化图将其于宽顶堰公式和肖克列奇公式得出的结果进行比较.宽顶堰公式Q m ax =m b 2g H +Αv 22g32式中,m 为系数;Α为水头损失系数;H 为坝前水深(m );v 为流速(m s );g 为重力加速度(m s 2).肖克列奇公式Q m ax =8g 14bH 3244数　学　的　实　践　与　认　识39卷式中,B 为坝体宽度;b 为决口宽度(m );H 为坝前水深(m ).通过计算,可以得到在相应溃口特征尺寸的最大流量分别为15500m 3 s 和13500m 3 s .与上述曲线结果比较,曲线中1小时左右溃口达到稳定的结果与这两个公式的计算结果最为接近.因此,可以大致确定在在堰塞湖蓄水漫顶1小时后溃口形态基本达到稳定.同时,这一小时内的最大溃口流量14000m 3 s 可以作为下游风险评估的基本条件.3.5　溃坝模型的灵敏度分析按照上述讨论的非恒定流数学模型,由图3的水位蓄水量关系,假定不同的溃坝水位,不同的残留坝体高程,计算了各种情况下的溃坝最大流量及其溃口洪水过程,其中,各状态下的溃坝最大流量情况如图14～图16.图14　747水位、不同残留坝高、不同溃口形成时间的溃口最大流量图15　745水位、不同残留坝高、不同溃口形成时间的溃口最大流量据图14、图15、图16各种情况的计算结果可见,如果溃口冲刷至弱风化顶部720m 高程后停止冲刷,以堰顶高程747m 作为起始溃坝水位,得到溃坝最大流量约为16500m 3 s ;以堰顶高程745m 作为起始溃坝水位,得到溃坝最大流量约为14000m 3 s ;以堰顶高程742m 作为起始溃坝水位,得到溃坝最大流量约为11000m 3 s .又当固定水位而改变残留坝高时,溃口最大流量的变化同样非常明显.因此,可以看出由于残留坝高和水位的变化对溃口最大流量的影响还是相当剧烈的.5416期曾向荣,等:唐家山堰塞湖泄洪问题研究图16　742水位、不同残留坝高、不同溃口形成时间的溃口最大流量4　泄洪阶段的模型建立与求解4.1　淹没区域计算模型基于数字高程模型(D E M )[11]求取给定水位条件下的淹没区,应当区分两种情形.第一种:凡是高程值低于给定水位的点,皆计入淹没区;第二种则需考虑“流通”淹没的情形,即洪水只淹没它能流到的地方.例如对于环形山(一种中间低洼、四周环形隆起的地形),第一种淹没计算可能导致环形山内外都生成淹没区,而在第二种淹没中,外来的洪水若未及山顶,只能在山环外形成淹没区.这两种情形都具有实际意义.第一种相当于整个地区大面积均匀降水的情形,所有低洼处都可能积水成灾;第二种相当于高发洪水向邻域泛滥,例如洪水决堤,或局部暴雨引起的暴涨洪水向四周扩散.通常,将这两种情形分别称为“无源淹没”和“有源淹没”.在有源淹没中,洪流不仅可能被阻于环形山外,也可能为任何高地所阻隔,一条流域的洪水一般不会越山而过在相邻流域形成淹没区.可以看出,唐家山堰塞湖溃坝洪水问题是有源淹没问题.本文采取种子蔓延算法处理唐家山溃坝洪水有源淹没问题.种子蔓延算法与前面讨论的连通性分析算法相似,是一种基于种子空间特征的扩散探测算法,其核心思想是将给定的种子点作为一个对象,赋予特定的属性,在某一平面区域上沿4个(或8个)方向游动扩散,求取满足给定条件、符合数据采集分析精度、且具有连通关联分布的点的集合.利用种子蔓延算法计算淹没区,就是按给定水位条件,求取满足精度、连通性要求的点的集合,该集合给出的连续平面就是所要求算的淹没区范围.满足水位条件但与种子点不具备连通关联性的其他连续平面,将不能进入集合区内.淹没区计算的准确性在很大程度上依赖D E M 的分辨率.每一个象元都代表着地面上一个区域,且每一个象元都拥有自己的高程值和精度.依据象元的高程值精度,选择恰当的阈值作为判断象元归属的条件,该值称为种子蔓延探测分辨率.设空间数据精度为m ,探测分辨率为w ,经分析和计算实验,其关系式可以表达为:w =1nm +K (K <m )式中n 为倍频系数,通常取2;K 为平衡参数,用以平衡探测分辨率和运算效率间的矛盾.种子蔓延算法计算淹没面积的流程图17所示.种子点的起始位置一般选在水库堤坝、河岸边界等特征点处,将满足所有条件的连通关联淹没点,存入缓存区,并不断地进行累加,从而使64数　学　的　实　践　与　认　识39卷淹没区域不断扩大.淹没区面积计算模型如图18,有源淹没区的淹没范围计算比较复杂,所有像元与种子点必须连通,程序沿4个方向寻找满足条件的像元.有源淹没所获得的区域,必须与地形图上的行政区划层进行叠加计算,才能最后得出淹没区的实地位置.事先应将矢量行政区划层转换成网格数据.图17　淹没区域计算算法流程图图18　淹没区面积计算模型由于预先得知坝体发生溃塌造成堰塞湖内1 3的蓄水突然下泻,不考虑下渗等作用的影响,洪水淹没容积是可以计算得到的.由资料可以得到唐家山堰塞湖的最大蓄水容量为3亿立方米,发生溃坝时堰塞湖蓄水量按3亿立方米计算.又由于下流河道以及堰塞湖会储存部分洪水,约为2270万立方米,因此淹没容积为:W ′=30000 3-2270=7730(万立方米) 在种子点的起始位置和水面的水位高程h 的情况下,使用种子蔓延算法可以得到洪水淹没面积、容积和洪水淹没水深分布.因此,可以通过改变水面高程h ,使洪水淹没容积为7730万立方米,从而得到洪水淹没区域.4.2　淹没面积以及人员撤离方案计算机实现得到1 3溃坝影响到的地区有北川县城、通口镇、含增镇和青莲镇、绵阳县城(见图19),淹没总面积为1412.64万平方米,其中绵阳市的淹没面积为481.10万平方米,人员撤离18万左右.4.3　洪水演进计算要求得洪水到达时间,必须对洪水演进进行推算,由于圣维南方程组的拟线性及水流的不连续性,溃坝问题的分析求解难度非常大,因此采取近似方法和经验公式[12214].1)溃坝下游流量计算.坝址下游沿程各河段控制断面溃坝最大流量采用李斯特万公式计算:Q LM =WW Q M +L V K7416期曾向荣,等:唐家山堰塞湖泄洪问题研究图19　各地区淹没范围式中,W 为堰塞湖蓄水量,Q LM 为距溃坝L (m )控制断面溃坝最大流量(m 3 s ),W 为堰塞湖蓄水量(m 3),Q M 为坝址最大流量(m 3 s ),L 为控制断面距堰塞湖坝址的距离(m ),V 为特大洪水的最大流速(山区取3～5m s ,丘陵区取2～3m s ,平原区取1～2m s ),K 为经验系数(山区取1.1-1.5,丘陵区取1.0,平原区取0.8～0.9).考虑唐家山以下河道形态,取V K =5.0.2)溃坝洪水传播时间和流量过程线.黄河水利委员会水利科学研究所根据实验得出了溃坝洪水传播时间及概化流量过程线.洪水到达时间计算t 1=K 1L 1.75(10-h 0)1.3W 0.2H 0.350式中,h 0为溃坝洪水到达前下游计算断面的平均水深,单位是m ,可以根据河道横切面粗略计算;K 1为系数,取为0.70×10-3;t 1是洪水起涨时间,单位是s .最大流量到达时间计算84数　学　的　实　践　与　认　识39卷t 2=K 2L 1.4W 0.2H 0.50h0.25M 式中,K 2为系数,取值0.8～1.2;h M 是最大流量时的平均水深,首先求出最大流量,再由水位流量关系曲线求得,单位是m ;t 2是最大流量到达时间,单位是s ,计算结果如表2所示.表2　洪水演进计算结果城镇名称所属县市与溃坝距离(km )洪水到达时间(h )最大流量到达时间(h )北川县城北川县2.10.1210.364通口镇北川县22.52.1254.225含增镇江油市28.13.0475.613青莲镇江油市48.75.2416.082绵阳县城绵阳市65.97.8568.510:[1]　唐家山堰塞湖抢险全记录[EB OL ].h ttp : www .new ssc .o rg ,200829228.[2]　地震灾区气象服务专报(第41期)[N ],200825225.[3]　四川省气象台.四川汶川地震灾区堰塞湖气象服务专报(第60期)[N ],20082622.[4]　中央气象台.四川省汶川天气预报[EB OL ].h ttp : www .sina .com .cn ,2008205217.[5]　四川汶川地震灾区气象服务专报(70期)[N ],2008206206.[6]　新华网直播车实时报道.唐家山堰塞湖抢险进展和人员撤离[EB OL ].www .xinhuanet .com .2008205231.[7]　长江委组织专家对唐家山堰塞滑坡体进行应急研究分析[EB OL ].h ttp : www .hw cc .com .cn ,200825221.[8]　陈凯.洪水淹没计算与评估系统[D ].硕士学位,2004.13218.[9]　刘仁义,刘南.基于G IS 的复杂地形洪水淹没区计算方法[J ].地理学报,2001,56(1):126.[10]　彭雪辉.风险分析在我国大坝安全上的应用[D ].硕士学位,2003.52254.[11]　吕杰堂,王治华,周成虎.西藏易贡滑坡堰塞湖的卫星遥感监测方法初探[J ].地球学报,2002,23:3632368.[12]　张振,张光科,刘超,杨中华.小型水库溃坝洪水演算探讨[J ].东北水利水电,2007,25(11):325.[13]　付磊,张洪明,姚激.水库调洪演算的数值解析法[J ].中国水利水电科学研究院《科研简报》特刊,2006,4(4):75277.[14]　中国水科院.四川汶川抗震救灾[R ].2008,7.The Research of Tangj i ashan Barr ier Lake Sp illwayZEN G X iang 2rong ,　HAO Hong 2x ing ,　SU N Bo 2liang(D epartm en t of Info rm ati on System and M anagem en t ,N ati onal U n iversityof D efen se T echno logy ,Changsha 410073,Ch ina )Abstract :　T he paper studied the sp ill w ay law of T angjiashan barrier lake ,and did model analysis of its w ater sto rage ,dam b reak and flooding .Keywords :　T angjiashan ;barrier lake ;sp ill w ay ;w ater sto rage ;dam b reak ;flooding 9416期曾向荣,等:唐家山堰塞湖泄洪问题研究。

堰塞湖泄洪问题

全国第五届研究生数学建模竞赛题目建立唐家山堰塞湖泄洪问题的相关数学模型合理性很好。

漫顶）。

溃坝风险大，最后，关键字一、问题的提出四川省汶川县5月12日发生里氏8.0级强烈地震，给人民生命财产和国民经济造成了极大的损失。

地震引发的次生灾害也相当严重，特别是地震导致大量山体滑坡、崩塌，形成三十多个高悬于灾区人民头上的堰塞湖。

如果堰塞湖水位不断上涨，溃决以后会造成下游大量沿江地区洪水灾害，波及范围更广。

与此同时，堰塞湖不仅阻断交通，还淹没上游道路，阻碍救援人员进入灾区，迫使救援人员只能徒步翻山越岭。

最令人担心的是由松散土组成的规模比较大的堰塞湖，这种堰塞湖相对比较容易溃决。

堰塞湖最好的处理办法就是逐步降低水位，以人工挖掘、爆破、拦截等方式来引流或疏通湖道，使其汇入主流流域或分散到水库，以免造成洪灾。

但由于此次的堰塞湖规模非常大，小的50万立方米以上，一般都在100万立方米、200万立方米，大的还有上千万方甚至上亿方。

如此大体积短期处理难度非常大，加上路不通，重型机械上不去。

因此需要做好风险预案，一边进行监测，随时观察变化情况，一边预估下游危险，避免一旦发生溃决造成太大损失。

此次的三十多个堰塞湖规模大小不等，影响范围不同，危险程度各异。

其中以唐家山堰塞湖尤为严重，被水利部抗震救灾指挥部前方专家列为1号风险。

唐家山堰塞湖坝底高710米，回水尾部高750米，最大水深41米，坝体体积约2037万立方米，集水面积3550立方千米，最大可蓄水量3.12 亿立方米。

坝区两岸基岩为寒武系下统清平组薄层硅质岩、砂岩、泥灰岩、泥岩组成，岩层软硬相间。

两岸坡分布有残坡积的碎石土，碎石土由粉质壤土、岩屑和块石组成。

堰塞湖坝体由基岩挤压或解体形成的碎裂岩、残坡积碎石土和苦竹坝库区沉积的含泥粉细砂组成。

该湖的堰塞体沿河流方向达800多米，从最终的实际情况看，从坝顶溢出而溃坝的可能性比其它原因溃坝的可能性大得多。

鉴于唐家山堰塞湖的威胁严重程度，加强对震后灾害规律的研究为国家抗震救灾提供更有力的科学支撑显得尤为重要。

堰塞湖防汛演练心得体会

堰塞湖防汛演练心得体会堰塞湖防汛演练心得体会堰塞湖防汛演练是为了提高应急准备和应对突发事件的能力，保障人民生命财产安全。

通过参与这次演练，我深刻体会到了防汛工作的重要性和团队合作的力量。

首先，充分认识到了堰塞湖防汛工作的紧迫性和重要性。

堰塞湖是一种极具威胁性的自然灾害，一旦出现堰塞湖崩溃泄洪，将给下游地区带来巨大的洪涝灾害。

在演练中，我们通过模拟情景感受到了水流的威力和冲击力，意识到了堰塞湖防汛工作的难度和风险。

只有提前做好准备，确保各项防汛设施的完好和作业人员的熟练，才能有效应对堰塞湖崩溃的危险。

其次，加深了我对团队合作的认识和理解。

堰塞湖防汛工作需要多个部门和不同职能的人员齐心协力完成。

在演练中，我们分工明确，各司其职，密切配合，形成了高效的工作模式。

特别是沟通和协调能力的重要性得到了充分体现。

只有大家共同努力，密切配合，才能应对突发事件，保障人民生命财产安全。

此外，还强调了科学决策和技术支撑的重要性。

在演练中，我们利用先进的技术手段进行了水流模拟和演练过程的模拟分析，为指挥决策提供了科学依据。

技术人员的专业能力得到了充分发挥，为防灾减灾工作提供了有力的支持。

只有依靠科学技术，才能更好地应对突发事件和提高防灾减灾能力。

最后，还要强调宣传教育的重要性。

当灾害发生时，应对充满不确定性和变数。

通过宣传教育，可以增加公众的风险意识，提高自防自救的能力，减轻灾害损失。

而在演练中，我们还可以利用宣传教育的机会，向公众展示堰塞湖防汛工作的重要性和团队合作的力量，进一步营造社会关注和支持防汛工作的氛围。

通过参与这次演练，我深刻体会到了堰塞湖防汛工作的复杂性和困难性，以及团队合作和科学技术的重要性。

我相信，只要我们充分认识到了这些问题，提高应急准备和应对突发事件的能力，我们就能更好地保障人民生命财产安全，为社会的稳定和发展作出更大的贡献。

唐家山堰塞湖溃坝洪水演进及下泄过程复演

唐家山堰塞湖溃坝洪水演进及下泄过程复演
张细兵;卢金友;范北林;朱永辉
【期刊名称】《人民长江》
【年(卷),期】2008(039)022
【摘要】将溃坝模型和传统一维非恒定流模型相结合,建立了溃坝洪水演进数学模型.采用DEM数字模型插值技术,快速获取了模型计算所需的断面地形,适合于资料缺乏情况下的应急模拟计算.采用该模型,进行了唐家山堰塞湖典型溃坝模式下的洪水演进计算和湖水实际下泄过程的复演计算.实测资料的验证结果表明,计算与实测的水位、流量过程符合较好,模型可靠.建立的模型及相关处理方法值得在堰塞湖应急除险及其它溃坝问题模拟计算中推广应用.
【总页数】3页(P76-78)
【作者】张细兵;卢金友;范北林;朱永辉
【作者单位】长江水利委员会,长江科学院,湖北,武汉,430010;长江水利委员会,长江科学院,湖北,武汉,430010;长江水利委员会,长江科学院,湖北,武汉,430010;长江水利委员会,长江科学院,湖北,武汉,430010
【正文语种】中文
【中图分类】P315.9
【相关文献】
1.尾砂库溃坝洪水演进过程模拟 [J], 周扬;张强
2.水库溃坝洪水与洪水演进过程分析计算 [J], 秦晨
3.坝体溃决过程与溃坝洪水演进耦合数值模拟 [J], 姜治兵;崔丹;程子兵
4.唐家山堰塞湖应急除险工程堰塞坝下游洪水演进复演验证 [J], 张细兵;卢少为;范北林;朱勇辉
5.基于HEC-RAS的百花滩水电站溃坝洪水演进过程及影响分析 [J], 蒋林杰;付成华;程馨玉;刘健
因版权原因，仅展示原文概要，查看原文内容请购买。

唐家山堰塞湖成功排险的探索与思考

唐家山堰塞湖成功排险的探索与思考
涂云;詹登民;方盛;龚军仓
【期刊名称】《水利水电技术》
【年(卷),期】2008(039)008
【摘要】@@ 1 概况rn唐家山堰塞湖位于四川省北川县上游通口河上,堰体长803.4 m,宽611.8 m,面积30.72万m2,体积2 037万m3,堰体最高点高程791.9 m,堰体底部高程667.5 m,堰体最大高度124.4 m,垭口处堰高82.65 m.总库容3.15亿m3,回水长度约23 km.若发生溃坝将影响下游北川县,江油市,绵阳市的涪城区、科学城、农科区以及三台县人口130多万,甚至影响到遂宁市,属"5.12"汶川特大地震中形成的极高危级堰塞湖.
【总页数】4页(P21-23,26)
【作者】涂云;詹登民;方盛;龚军仓
【作者单位】中国人民武装警察部队,水电第三总队,四川,成都,610036;中国人民武装警察部队,水电第三总队,四川,成都,610036;中国人民武装警察部队,水电第三总队,四川,成都,610036;中国人民武装警察部队,水电第三总队,四川,成都,610036
【正文语种】中文
【中图分类】P642.2(271)
【相关文献】
1.唐家山堰塞湖应急排险设计及综合整治的设想 [J], 陈五一;郑家祥;谭建;唐朝阳
2.唐家山堰塞湖排险设备保障方案 [J], 陈同俭;范双柱
3.浅谈唐家山堰塞湖应急排险的安全管理 [J], 涂云
4.唐家山堰塞湖成功排险综述 [J], 李虎章;刘松林
5.惊心动魄30天——唐家山堰塞湖排险镜头回放 [J], 苏竹青;罗勤（图）;王永（图）
因版权原因，仅展示原文概要，查看原文内容请购买。