复杂边界及实际地形上溃坝洪水流动过程模拟

建立了基于无结构三角网格下采用有限体积法求解的二维水动力学模型,用于模拟溃坝洪水在复杂边界及实际地形上的流动过程。该模型采用Roe格式的近似Riemann解计算界面水流通量,结合空间方向的TVD-MUSCL格式及时间方向的预测-校正格式,可使模型在时空方向具有二阶计算精度。模型中引入最小水深概念,提出了有效的干湿界面处理方法。模拟了理想条件下溃坝水流过程,研究不同最小水深取值对干河床上洪水演进的影响,并用两组简单溃坝水流的水槽试验资料对模型进行验证。采用该模型模拟了实际溃坝洪水的流动过程,所得计算结果与实测资料及已有模型计算结果较为符合。     

基于结构网格的溃坝水流数值模拟

针对溃坝水流数值模拟面临的复杂地形和不规则边界等问题,基于结构网格建立了适应复杂地形和不规则边界的溃坝水流数值模拟有限体积模型（HydroM2D）。模型基于具有守恒特性的二维浅水方程,利用HLLC格式的近似Riemann解计算网格界面通量,利用MUSCL-Hancock法不断向前积分,使模型在时空上具有二阶精度;对源项进行离散处理确保模型的稳定性;模型引入有效干湿边界和不规则地形边界处理方法,准确模拟了干湿单元的动态交替和复杂边界上的水流特性。最后分别利用水槽试验、物理模型和实际算例对模型进行验证。结果表明,该模型对不同情景下的溃坝洪水模拟结果和实测资料以及现有模型模拟结果具有较高的一致性,模拟精度较高,稳定性较好,具有推广应用价值。     

具有复杂计算域和地形的潮汐流动数值模拟

采用无结构网格上的Roe型二阶精度迎风型FVM格式的有限体积方法对具有复杂计算域和地形的胜利油田海域的潮汐流动进行数值模拟。为保证Roe格式能够应用于复杂地形条件下的计算,采用Ropers格式来解决通量梯度项与源项的平衡问题,对摩擦力源项采用分步法求解以增加格式的稳定性。应用此方法对地形复杂的胜利油田海域进行了实际模拟,数值计算结果和实测结果吻合较好。     

基于改进SPH模型的溃坝洪水演进模拟方法

溃坝洪水演进模拟的准确性是制约水库洪水预演有效性的关键。基于光滑粒子流体动力学(Smoothed Particle Hydrodynamics, SPH)方法提出了适用于溃坝洪水演进分析的数值模拟方法。通过设置溃口粒子与粒子库,基于黎曼不变量对SPH粒子状态进行修正,构建施加边界条件的改进SPH溃坝洪水演进模型,将SPH瞬时全溃整体模型转换为考虑溃口水流变化的入流边界模型,实现SPH方法与溃口计算模型的耦合。以Malpasset溃坝事件为例,检验了该模型计算溃坝洪水的精度,结果表明该模型精度相对较高,与实测值吻合较好;应用该模型模拟了某水库溃坝洪水演进预演过程,评估其对下游输水干渠及交叉建筑物排水倒虹吸的洪水冲击风险,结果表明在上游水库遭遇超标准洪水漫顶溃坝工况下,洪水演进至排水倒虹吸处的最大洪水位未超过校核洪水位。改进SPH模型精度高,可靠性强,与溃口计算模型耦合性好,可作为溃坝洪水演进模拟的通用手段之一。     

基于非结构混合网格求解二维浅水方程的一种数值方法

复杂边界下的流场数值模拟常基于非结构化网格进行求解,建立一种在非结构混合网格上求解水深平均的二维浅水方程的模型,以便精确模拟复杂边界、提高计算效率。该模型时间项离散采用隐格式使得模型具有较好的稳定性,对流项和扩散项分别采用总变差减小(Total Variation Diminishing,TVD)格式和构造辅助点的方法来离散,同时采用水深平均的标准k-ε模型来封闭湍流模型。选用两个经典验证算例检验模型,计算结果表明,基于非结构混合网格开发的模型具有较高的精度,且收敛性能较好。     

复杂边界条件多洪源防洪保护区洪水风险分析

为了定量获取防洪保护区在多洪源和复杂边界条件下的溃堤洪水风险信息,以非恒定流控制方程为理论基础,建立了多洪源一维河网水动力学模型和防洪保护区二维洪水演进模型,利用溃坝模型实现河道与保护区的耦联,并采用局部网格加密和相似建筑物模拟等方法处理保护区内道路等复杂边界的导阻水作用。利用所建模型模拟了长江、汉江和东荆河3种不同洪水来源, 在4种不同位置溃堤情况下汉南至白庙长江干堤防洪保护区的洪水淹没情景,采用基于淹没水深的损失率关系法对比分析了4种计算方案的淹没面积、经济损失和受灾人口。结果表明:模型构建合理、稳定性和适应性好,复杂边界对洪水演进过程影响明显,不同洪源溃堤情形的风险信息差异较大;在计算条件下,以长江发生1954年型300年一遇洪水向新溃口情形下的淹没损失最严重,其淹没面积达3 790 km2,受灾人口为196.8万人,经济损失约802亿元。研究成果可为洪水风险管理与避洪转移决策提供有力的技术支撑。     

漫溃堤洪水联算全二维水动力模型及应用

为应对河道洪峰流量增大和漫溃堤长历时相伴发生的洪灾现象,借鉴全二维气相色谱理论提出全二维水动力模型概念,建立了模拟河道和灌区洪水演进的漫溃堤洪水联算全二维水动力模型,并采用Roe格式的近似Riemann解对界面通量进行数值求解。模型内通过漫溃堤堰流公式成功实现河道与灌区的耦联,考虑溃口展宽变化,加密处理河道网格,采用热启动与干湿水深理论对模型进行优化,并利用加大糙率法对村庄较为密集的地形进行优化处理,尽可能反应地面真实情况。将该模型应用于黄河宁蒙段河道与左右岸灌区的漫溃堤洪水演进模拟,计算结果合理可靠,流场分布均匀光滑,初步验证了模型的精度及可靠性,研究成果对河道溃决洪水的精细仿真模拟和该地区洪水风险分析决策具有重要意义。     

基于Godunov格式的溃坝水流数学模型

为了更好地把握溃坝洪水风险,减小因溃坝洪水而造成的人员生命和财产损失,建立了基于Godunov格式的一维、二维溃坝水流耦合数学模型。一维溃坝水流模型采用HLL格式的有限体积法求解,二维溃坝水流模型采用基于非结构网格的Roe格式离散求解,在一维、二维模型的链接处采用重叠计算区域的方法实现一维模型和二维模型之间的水力要素信息交换。经弯道溃坝算例和断面突变溃坝算例验证,该耦合模型具有良好的可靠性和适用性,验证后的耦合模型为大尺度的溃坝水流数值模拟打下了基础。     

溃坝洪水分析及影响研究——以陕北地区瑶镇水库为例

为研究溃坝洪水带来的影响,以陕北地区的瑶镇水库为例,采用洪水演进模拟和GIS技术等先进的溃坝分析技术,对水库遭遇五十年一遇洪水下右岸与非溢流坝连接的副坝段溃坝洪水演进过程、水深及淹没范围、流速分布及流速矢量分布展开评估,重点对溃坝洪水对水库下游产生的影响进行分析,从而得出相应结论。研究结果认为:受溃坝洪水影响最大的是瑶镇地区,水流最大流速均发生在主河道区域,五百年一遇洪水淹没范围与五十年一遇洪水淹没面积相比增加明显。瑶镇水库溃坝不会对下游的采兔沟水库安全产生很大影响。研究结果对本区域溃坝洪水及未来研究提供了参考依据。     

易贡滑坡堰塞湖溃坝洪水分析

滑坡堰塞坝体主要由块石、碎石土等松散材料组成,随着上游水位的不断上升,极易失稳,一旦决口将对给下游人民的生命财产安全造成极大的威胁。因此,研究堰塞坝溃坝问题具有重要的学术意义和应用价值。2000 年 4 月 9 日,西藏林芝地区波密县易贡藏布河扎木弄沟发生大规模山体滑坡堵塞易贡藏布江,形成坝高60m,长约2500m,库容可达288108m3,体积约28108～30108m3的滑坡堰塞湖, 2000年6月10日堰塞坝溃决。本文以易贡堰塞湖溃坝为例,从连续性方程及Navier Stokes方程出发,结合标准型湍流模型,并采用VOF方法进行自由面处理,基于流体计算软件Fluent模拟分析了溃坝洪水在下游弯曲河道的演进过程及不同位置的流速变化。数值模拟结果与实测资料记录基本一致,表明该模型能够模拟溃坝洪水在地形复杂弯曲河道中的演进过程。     

Three-dimensional flood routing of a dam break based on a high-precision digital model of a dense urban area

Dense buildings are the major factor affecting urban flood routing. Currently, the study of urban dam-break flood routing primarily focuses on a simplified terrain model and 2D shallow water equations, which ignore the effects of dense urban buildings. Furthermore, the complex interactions between the dam-break wave and the wall surfaces of buildings are not reflected in the results. To tackle these problems, three-dimensional flood routing with a high-precision digital model of an urban area is studied in this paper. Firstly, the vector data of various land types is extracted from a remote sensing image, and the NURBS algorithm based on the TIN algorithm is introduced to construct a three-dimensional terrain model. Coupled with the vector data and the terrain model, a three-dimensional digital model of the urban area is established. Next, a three-dimensional \(k{-}\varepsilon\) turbulence model is proposed for the flood routing simulation. A polyhedral grid with a second-order accuracy and a discrete format is used to divide the digital model of the urban area, and the governing equations are solved using the PISO algorithm. Finally, the superiority of the 3D mathematical model and the computational efficiency of the polyhedral mesh model are validated according to the urban flood routing experiments of Testa and Soares-Frazão. An urban reservoir, located in SZ City, China, is modelled to show that a dam-break flood in an urban area exhibits significant three-dimensional characteristics. Moreover, due to the surrounding buildings, the flood exhibits complex three-dimensional turbulence phenomena, including collision, reflection and vortices.     

Study of the comprehensive risk analysis of dam-break flooding based on the numerical simulation of flood routing. Part I: model development

Dam-break floods have been of increasing concern to safety engineers and decision makers. The presence of complex terrain in inundation areas multiplies the simulation difficulty of flood routing. In previous studies, representing the flood routing parameters empirically does not reflect the characteristics of flood routing, which strongly influences the accurate assessment of the dam-break consequences. The basis for carrying out dangerous reservoir reinforcement is just engineering safety, without considering the actual situation of downstream areas. In this study, a comprehensive risk analysis of the dam-break flood was implemented based on the numerical simulation of flood routing. First, coupled with the volume of fluid method, a three-dimensional k–? turbulence mathematical model was developed for flood routing in complex inundation areas. Then, based on the flow parameters obtained through computational fluid dynamics modeling, the attribute measure methodology was used for the evaluation of consequences combined with the calculation of the dam-break consequences (loss of life, economic loss, social and environmental influence). Furthermore, a methodology containing the combined weight method and the technique for order performance by similarity to ideal solution method was proposed for risk ranking of dangerous reservoirs due to its logical consideration of scalar values that simultaneously account for both the best and worst alternatives. Finally, a sensitivity analysis was performed to provide information about the stability of risk ranking. The aforementioned model and methodology are applied to a case involving five reservoirs in the Haihe River Basin in China for Part II of this study.     

云南山区水库溃坝洪水及其演进分析

受地形影响,西南山区水源以水库为主,由于大部分水库修建时间较早,在西南季风气候降水集中影响下,存在较大的溃坝风险。溃坝洪水突发性与破坏性极强,进行水库溃坝洪水计算及洪水演进分析,是在水库大坝发生突发性安全事故时科学应对的基础。结合西南山区实际,优选溃坝洪水计算与洪水演进模型,并以云南省昌宁县河西水库为例,分析确定相关参数,分析指出了西南山区水库溃坝形态以全溃为主,具有溃坝洪水量极大、洪水演进迅速的特点。研究不仅可为河西水库制定大坝安全管理应急预案提供技术支撑,也能为该地区水库溃坝洪水及其演进分析提供参考借鉴。     

基于力学过程的蓄滞洪区洪水风险评估模型及应用

为定量评估分蓄洪工程启用过程中蓄滞洪区的洪水风险等级,创建了基于力学过程的蓄滞洪区洪水风险评估模型。该模型采用二维水动力学模块计算蓄滞洪区的洪水演进过程,利用洪水中人体跌倒失稳公式及洪水中房屋、农作物损失的计算关系式,评估各类受淹对象的洪水风险等级。然后将二维水动力学模块计算的洪水要素与两个物理模型试验值进行对比,表明二维水动力学模块的计算精度良好。最后计算了荆江分洪工程启用时分洪区内洪水的演进过程,并评估洪灾中群众的危险等级和财产损失。计算结果表明:洪水演进至140 h时,蓄滞洪区群众、房屋、水稻和棉花的平均损失率分别为85%、59%、63%和72%。模型中提出的采用基于受淹对象失稳机制的洪水风险分析方法,比以往经验水深法划分风险等级的适用性更好,不仅能为洪水风险管理及蓄滞洪区启用标准制定提供参考,也能推广应用于溃坝或堰塞湖溃决等极端洪水风险评估。     

Modelling flash flood propagation in urban areas using a two-dimensional numerical model

This paper reports on the numerical modelling of flash flood propagation in urban areas after an excessive rainfall event or dam/dyke break wave. A two-dimensional (2-D) depth-averaged shallow-water model is used, with a refined grid of quadrilaterals and triangles for representing the urban area topography. The 2-D shallow-water equations are solved using the explicit second-order scheme that is adapted from MUSCL approach. Four applications are described to demonstrate the potential benefits and limits of 2-D modelling: (i) laboratory experimental dam-break wave in the presence of an isolated building; (ii) flash flood over a physical model of the urbanized Toce river valley in Italy; (iii) flash flood in October 1988 at the city of Nîmes (France) and (iv) dam-break flood in October 1982 at the town of Sumacárcel (Spain). Computed flow depths and velocities compare well with recorded data, although for the experimental study on dam-break wave some discrepancies are observed around buildings, where the flow is strongly 3-D in character. The numerical simulations show that the flow depths and flood wave celerity are significantly affected by the presence of buildings in comparison with the original floodplain. Further, this study confirms the importance of topography and roughness coefficient for flood propagation simulation.     

非粘性沙坝溃决数学模型及求解

考虑一维陡坡浅水方程与Exner方程耦合构建漫顶溃坝数学模型,利用Roe格式黎曼近似解法求解,并用预测-校正法处理源项保证水流、泥沙侵蚀输运和底床演变过程的计算稳定.模型实例计算结果表明,该方法能较好地计算非粘性模型沙坝漫顶溃决问题,可为溃坝数值计算提供新的思路.     

瞬时溃坝最大流量计算新通式推导及验证

为了计算瞬时溃坝最大流量,基于堰流与波流量相等原理,建立了瞬时溃坝最大流量与堰流关系,理论推导得到同一公式可以计算大坝全溃、横向局部溃坝、垂向局部溃坝及横垂向局部溃坝的瞬时最大流量的新通式,给出了该通式的基本和上限、下限及一般等通式。研究发现:瞬时溃坝流量与溃口堰坎类型有关,可以把瞬时溃坝最大流量、逐渐溃坝洪水或漫坝洪水的计算公式和相应系数统一到不同堰坎类型的堰流量计算中;而这些新通式相应系数的取值是被实验确认和经典著作肯定的,保证了系数取值的可靠和准确。这些新通式容易理解、形式简单、计算简便,参数取值可靠,把溃坝洪水计算的复杂问题转变为处理经典水力学的简单堰流计算。同时通过对比分析,充分论证了一般通式的合理性,并用国内外溃坝实例验证了计算方法和成果较科学、可靠、合理。     

Study of the comprehensive risk analysis of dam-break flooding based on the numerical simulation of flood routing. Part II: Model application and results

The present model and methodology described in Part I of this work are applied to perform a comprehensive risk analysis of the dam-break flood of five reservoirs in the Haihe River Basin, China. The results indicate that the three-dimensional numerical simulation considering complex terrain can reflect the characteristics of flood routing and the three-dimensional phenomena. Based on the simulation results, it can be concluded that the risk grades of the consequences induced by a gradual or instantaneous dam break of the Dongwushi reservoir are extremely serious, as determined through the attribute synthetic approach. The results obtained from ranking the risk by the technique for order performance by similarity to ideal solution method are that the Dongwushi reservoir has the most serious consequences when the dam breaks followed by the Lincheng reservoir, the Miaogong reservoir and the Yunzhou reservoir, and the Youyi reservoir has the least severe consequences. Though the ranking of the relatively comprehensive risk coincides with that of the consequences, the dam safety measured by the dam failure probability plays an important role in ranking the risk. A sensitivity analysis is performed by individually increasing the weight of each criterion by 20 %, and the ranking order is not changed, suggesting that the evaluation model is reasonable, stable and reliable.