基于MIKE FLOOD的城区溃坝洪水模拟研究
王欣;王玮琦;黄国如
【摘 要】The safety of dam not only affects the benefit of the project,but also relates cloly to the safety of people's life and property.Dam-break flood simulation can evaluate the impacts of dam-break,and it is of great significance to the making of emergency plans so as to control flood and reduce disasters.Rearches are carried out on the dam-break flood developing process in the downstream of the Minzhi rervoir bad on the MIKE FLOOD model,which couples MIKE 11 and MIKE 21.By adopting two methods,a sudden dam-break method (for sudden partial dam-break as well as sudden full dam-break) and a gradual dam-break method,simulations of the flow graphs at the breach dam under four operating situations and the flood developing process in the downstream are respectively made in this study.The rearch results show that the peak flow of the sudden dam-break is larger at the beginning of the dam break,and the gradual dam-break peak flow is relatively small when the deformation of epage failure develops to the upper part of the dam.And th
en,as the water level in the rervoir falls gradually,the flow discharge becomes smaller,until the rervoir is emptied.The dam-break flood has a greater destructive power against the upstream area of the Hengling village,where the submerged depth is larger.The flood velocity is up to 5 m/s in the residential and commercial areas of the middle reach of the Minzhi River.The flood has a certain damage to the buildings.The left side of the Xiangnan village is low-lying,the flooding is the most rious,and it still has a waterlogging depth of 3 m after the flood subsidence.%大坝安全不仅影响工程效益,还影响人民的生命和财产安全,溃坝洪水模拟可以对水库大坝的失事影响做出评估,对制定应急预案和防洪减灾具有重要意义.以深圳市龙华新区民治水库及下游片区为研究对象,基于MIKE FLOOD将MIKE 11模型和MIKE21模型进行动态耦合,对溃坝洪水在下游的演进过程进行仿真模拟.模型采用瞬间溃(瞬间部分溃和瞬间全溃)以及逐渐溃两种溃决方式,分别模拟4种工况下的溃口流量过程线以及下游洪水演进过程.结果表明:瞬间溃的洪峰流量较大,出现在溃坝开始时刻,而逐渐溃的洪峰流量相对较小,出现在渗透破坏变形发展至上部坝体坍塌时刻,之后均随库区水位逐渐降低,下泄流量变小,直至库区水体排空.溃坝洪水对上游地区横岭村附近破坏较大,淹没水深较深.民治河中游段居民和商业区附近洪水流速接近5 m/s,对建筑物
有一定破坏力,左侧向南村地势较低,淹没情况最为严重,并且在洪水消退后仍有3 m左右积水.民治河下游地区在洪水消退后也有少量积水.
【期刊名称】《水利水运工程学报》
【年(卷),期】2017(000)005
【总页数】7页(P67-73)
【关键词】溃坝洪水;MIKE FLOOD;数值模拟;淹没水深;淹没范围
【作 者】王欣;王玮琦;黄国如
【作者单位】华南理工大学土木与交通学院,广东广州510640;华南理工大学土木与交通学院,广东广州510640;华南理工大学土木与交通学院,广东广州510640;华南理工大学亚热带建筑科学国家重点实验室,广东广州510640
【正文语种】中 文
【中图分类】TV122+.4
水库大坝是人类历史上重要的水工建筑物，在兴利和调蓄洪水等方面发挥着不可替代的作用，同时对区域内国民经济发展也承担着重要角色[1]。大坝一旦发生溃决，将会给城市带来灾难性破坏，造成不可估量的损失[2-3]。近些年随着计算机和数值计算方法的迅速发展，溃坝洪水演进在复杂计算域和地形条件下的高性能数值模拟研究成为国内外学术界和工程界关注的前沿研究领域之一[4]。如何使用计算机仿真技术科学准确地预测、模拟和显示洪水淹没范围，对于防洪救灾和损失评估具有十分重要的意义。不少学者在这方面进行了深入研究并取得了一定成果[5-9]。一维模型对溃坝洪水演进进行模拟研究时，虽计算效率高、适应性好，但在模拟洪水漫过河堤进入洪泛区演进时具有很大难度。二维模型则在处理复杂水流在无固定路径的地表演进时，能够描述洪水淹没的时间、范围、流速、水深分布等信息，且具有良好精度；但在模拟洪水在河道内演进时，计算效率及适应性远不如一维模型且需要详尽的河道地形资料。鉴于此，本文以深圳市龙华新区上游民治水库为研究对象，利用MIKE FLOOD将一、二维模型进行动态耦合，发挥一、二维模型模拟洪水在河道中及地表演进时各自的优势，同时回避其计算效率、网格精度及适用性等问题，更为合理真实地模拟上游民治水库大坝溃决后溃坝洪水的演进过程。
MIKE FLOOD是一个耦合的水力模型，能够完整模拟一维地下排水管网系统水流过程和二
维地表漫流过程。MIKE FLOOD集成了MIKE Urban CS，MIKE11及MIKE21三个独立的软件模块。根据不同的应用情境可将其中的MIKE Urban CS或者MIKE11与MIKE21进行动态耦合，耦合后的MIKE FLOOD具有以下优点：能够在暴雨洪水过程中，模拟一维河道洪水超出河堤后，从决堤口向淹没区域的二维洪水演进过程；能够模拟一维河网与二维模拟区域连接处的动量传输；能够与外部软件如ArcGIS等进行复杂的数据交互，通过外部专业软件对其输入数据进行预处理等[10]。
MIKE11用于模拟一维河道水体的流态，集成了水动力模块、降雨径流模块、构筑物模块以及溃坝模块等，几乎涵盖了河流模拟的各个方面，在模拟过程中采用六点Abbott-Ionescu有限差分格式对圣维南方程组求解[11]。MIKE21属于平面二维自由表面流模型，广泛运用于海洋、湖泊、河道及蓄滞洪区的流场、流速、水位等方面的模拟，能够获得不同水文要素的时空分布及洪水淹没信息[10]。
本文利用MIKE11与MIKE 21模型在MIKE FLOOD平台上进行动态耦合，采用侧向连接方式[12](见图1)。侧向连接即河道通过两岸与二维模型计算区域进行水流交换，它允许MIKE21网格单元从侧面连接到MIKE11的部分河段甚至整个河段。
侧向连接中一维河道和二维网格单元间的水量交换采用堰流公式近似计算[13]，计算方法如下：
式中：q为交换水量(m3/s)；W为宽度(m)，一般取单元格和河道相连的边长；C为堰流系数；k为堰指数；Hus和Hds分别为堰上游和下游水位(m)；Hw为堰顶高程(m)。本次研究选择MIKE11中的河岸标记作为堰的位置。
民治水库位于深圳市龙华新区民治街道民治社区东南方约2.5 km，在2个小型水库民乐水库及雅宝水库(原名油柑水库)的下游，包含2个小水库在内的总集雨面积为4.5 km2。民治水库正常库容271 万m3，总库容400.2 万m3，正常蓄水位79.58 m，校核水位82.36 m，死水位69 m，坝顶长320 m，其设计标准为50年一遇设计，100年一遇校核。民治街道位于民治水库下游，为龙华新区与南山、福田、罗湖、龙岗五区交汇处。下辖9个社区工作站，21个社区居委会，总面积26 km2，人口42.7万。根据地形条件以及可能的淹没区域，确定研究范围扩展至民治水库下游4.25 km2。为分析受灾情况选取4个特征点进行分析比较(图2)，分别为民治水库和民治河上游间重灾居民区横岭村(特征点1)，民治河道中游西侧向南村(特征点2)，民治河下游水尾新村附近(特征点3)，以及民治河道中游段居民区和商业区附近(特征点4)。
利用MIKE11构建一维模型，从ArcGIS中提取河道的相关地理信息，按MIKE11河网文件的格式导入MIKE11中，并设定河道名称、里程数、地理标识等信息，制成河网文件。从河道横断面CAD图中提取断面资料，并根据里程数标志在河道上，对河道赋予断面信息。由于民治水库通过一条暗涵与下游民治河相连，其过流能力相对于洪水在地表演进时的流量很小，因此在进行本次溃坝洪水模拟时，根据其过流能力，在河道上游边界给定3 m3/s的恒定入流，河道水位作为下游开边界。初始水深设为1 m，初始流量为0.1 m3/s，民治河断面较为规整，水流较为通畅，两侧岸壁为土砂或石质，由于没有实测资料加以验证，本文根据天然河道糙率取值经验，取糙率为0.025。将上述所有文件导入模拟文件编辑器，设定模拟时间步长、结果输出文件名等，运行MIKE11，至此一维模型建立完毕。
利用MIKE21构建二维模型，根据计算范围提取计算边界以及河道两岸边界，制作成Mesh文件所需要的格式，导入非结构化网格生成器中。本文采用非结构三角形网格，网格划分完成后，节点数为9 314个，网格总数为17 911个，最大网格为400 m2，最小网格为79.2 m2，平均网格大小237 m2。导入研究区域的高程散点，根据网格划分进行插值，得到模拟区域的地形文件。研究区域边界取无滑移闭边界，闭边界不与外界进行水量交换，即法向流速为0，上边界以溃口流量过程作为流量边界条件。
将独立的MIKE11模型文件和MIKE21模型文件建立完成并可正常运行之后，将两者在MIKE FLOOD平台上进行耦合连接。首先导入需要连接的两个模型文件，然后选择侧向连接方式。非结构化网格与MIKE11进行侧向连接时，需要设定连接处的耦合线，以确定MIKE21与MIKE11的连接位置。本文将二维地形中一维河道的边界坐标导出，制作成耦合线所需要的格式，保证水量交换点即为河道两侧堤岸，同时，沿程河道的所有计算点都将与MIKE21相连。通过耦合线进行连接后，一维河道模型左侧岸线与72个MIKE21网格单元进行连接，右侧岸线与69个MIKE21网格单元进行连接。
鉴于本次模拟的民治水库主坝为土石坝，发生逐渐溃的概率较大，同时由于潜在的地震、滑坡等因素引发的瞬间溃的危害更大，因此，本文采用瞬间溃和逐渐溃两种溃决方式进行模拟。对于土石坝瞬间溃的溃口流量计算，根据相关资料[14]，本文采用下式进行溃坝洪峰流量计算：
式中：Qmax为溃决最大流量(m3/s)；g为重力加速度(m/s2)；B为坝长(m)；bm为最终溃口宽度(m)；H0为溃决时水深(m)。
对于逐渐溃的溃口流量计算，假设大坝溃决时，坝体先渗透变形破坏发展到一定程度后再
瞬间局部破坏。坝体渗透变形先从小孔开始，逐渐扩展，此时溃坝洪水流量根据文献按下式计算[15]：
式中：H为水库水位(m)；A为水流流过小孔的断面面积(m2)；Hp为溃口小孔处的高程(m)；f为达西摩擦系数，根据Moody曲线由d50粒径计算；L为小孔沿水流方向的长度(m)；D为小孔宽度(m)。假设小孔初始形状为矩形，宽度逐渐线性扩展到30 m后再瞬间部分溃决，溃决后的流量再根据式(3)进行计算。
上述两种溃决方式的溃口流量过程均采用下式计算：
式中：Q为溃口流量(m3)；V为库容(m3)；t为时间(s)。
分别模拟计算4种工况下溃口流量过程线以及下游洪水演进过程，具体工况设定见表1。
根据表1中的4种不同工况，采用式(2)～(4)，分别推算出各种工况下溃口洪水流量过程线，见图3。由图3可知，对于瞬间溃，当溃口宽度为100和320 m时，最大下泄流量分别为4 220, 10 098和13 912 m3/s，库容水量排空历时约52, 22和23 min。可见在溃坝初始时刻流量即达到最大，随着溃坝洪水下泄，水库水位不断下降，洪水流量也快速下降，并且溃口
宽度越大，溃坝初始时刻的最大洪水流量也越大，其下泄时长越短。而对于逐渐溃，当溃坝开始后，随着渗透变形管道小孔不断扩大，下泄流量也不断增加，溃口发展到一定程度后管道上部坝体开始垮塌，类似于瞬间溃，溃口流量迅速增大达到最大值，接着流量变小，直至库区水量排空，其中渗透变形破坏阶段最大流量为1 206 m3/s，之后迅速垮塌，最大流量达3 632 m3/s，总历时约68 min。各工况下洪水流量参数见表2。
