[发明专利]基于ROE格式的一维溃坝洪水模拟方法

说明书

本发明公开的基于ROE格式的一维溃坝洪水模拟方法，在计算机中读取计算区域的基础数据并对数据进行处理；采用Godunov格式的有限体积法离散圣维南方程，采用预估‑校正法对圣维南方程进行时间积分，计算在半时间水平上计算出网格平均水深，流速，流量等水力要素；采用ROE格式的数值通量函数求出各个网格左右界面上的流量和动量；采用校正器将网格平均水力要素推进到下一个完整的时间步长；最后输出模型模拟的结果文件。本发明的基于ROE格式的一维溃坝洪水模拟方法，计算准确，用时较短，可以模拟复杂的溃坝水流，并实现大尺度的溃坝水流模拟。

技术领域

本发明属于数值模拟技术领域，具体涉及基于ROE格式的一维溃坝洪水模拟方法。

背景技术

溃坝洪水是指堤坝或其他挡水建筑物瞬时溃决，发生水体突泄所形成的洪水。其形成原因一般可分自然和人为两大因素。如超标准洪水、冰凌、地震等导致大坝溃决属于自然因素；设计不周、施工不良、管理不善、战争破坏等导致大坝溃决属于人为因素。近年来，随着全球气候的不断变化，极端降雨等事件发生的频率不断增加，由此引发的超标洪水对大坝的威胁也在相应增加。准确、及时的溃坝洪水计算是对水库和堤防的失事影响做出定量估算，并合理确定水库或堤坝防洪设计标准以及避险措施的有效手段。

随着计算机技术的发展，使用数值模型手段模拟溃坝洪水过程逐渐进入人们的视野。通过数值模型可以在短时间构建并便捷计算不同初始溃坝条件下的结果，且可更加精细地描述水流在下游演进的物理过程，使得数值模型在溃坝洪水的模拟中得到了更广泛的应用。溃坝洪水计算的主要内容是根据溃坝坝址处的流量和水位过程线,求得下游洪水演进过程中沿程各处的流量、水深、流速、波前和洪峰到达的时间等。但由于溃坝水流的构成十分复杂，通常包含激波，亚临界流，超临界流等水流流态，导致对其进行数值模拟时容易出现数值振荡和不稳定的现象。此外，由于溃坝洪水在下游的演进过程尺度较大，采用二维模型计算耗时较长，难以满足及时评估洪水风险的需求。

发明内容

本发明的目的在于提供基于ROE格式的一维溃坝洪水模拟方法，以实现对溃坝洪水演进过程的准确、高效模拟。

本发明所采用的技术方案是：基于ROE格式的一维溃坝洪水模拟方法，具体按照如下步骤实施：

步骤1、在计算机中读取历史数据，所述历史数据包括河道断面数据以及曼宁系数，溃坝流量过程，计算时长，并根据河道断面间距将计算区域划分为不同大小的网格，然后对河道断面数据进行处理，获取各个网格的水深-过流面积、水深-湿周关系；

步骤2、采用Godunov格式的有限体积法离散圣维南方程，得到通量项、底坡源项、摩阻源项和时间步进；

步骤3、采用预估-校正法对离散后的圣维南方程进行时间积分，计算出步骤1中每个网格上的水力要素，所述水力要素包括平均水深和流速，再根据步骤1中所求水深-过流面积关系获得对应的过流面积；使用Superbee限制器对水深、流速、过流面积进行限制，避免出现异常的梯度，然后采用MUSCL型格式外推构造每个网格左右界面上相应的水深、流速、过流面积，使得求解具有二阶空间精度；

步骤4、根据步骤3中所求得的每个网格上的平均水深、流速、过流面积以及网格左右界面上相应的水深、流速、过流面积，采用ROE格式的数值通量函数求出各个网格左右界面上的流量和动量；

步骤5、对步骤3求得的每个网格上的平均水深、流速、过流面积采用半隐式计算方法进行更新校正，并将其推进到下一个时间步长；

步骤6、将步骤4和步骤5中计算的流量、动量、平均水深、流速、过流面积根据时间步长进行更新并输出，获得各网格单元在各个时刻的水力要素值。

本发明的特点还在于，

步骤1中，对河道断面数据进行处理的具体方法为：