黄河溃堤过程数学模型及其模拟方法.pdfVIP

中目料茔(E辑) 第32卷第5期 SCIENCEINCHINA(SeriesE) 2002年I(1H 黄河溃堤过程数学模型及其模拟方法4 梁林倪晋仁 『vg A 北京大学环境工程系水壮科学救育部重点实验室．北京10(1871’ A．GL．Borthwick B．D．Rogers (牛津大学工程科学系．OXl3PI，英国。 摘要 溃堤卸溃坝都属于非匣定含间断的浅水问题．但两者既有类似又有不霹之处． 渍坝p算侧重于出流过程，溃堤计算则还需要考察口门水流特征和形态变化过程．鉴 于将现行渍坝p算方法用于溃堤计算时会遇到许多困难，基于四叉树周格，采用 Godunov量有限体积方法，联台浅水方程、泥沙输穆方程、河床变形方程眦及口门横向 展宽计算方法，建立了适用于溃堤过程模拟的数学模型，并采用黄河典型河段的实际 地形和水流、泥沙条件进行了应用研究，得到了黄河溃堤过程中可能出现的口门水流特 征和形态的典型变化过程． 关键词 黄河溃堤 间断浅水方程Godunov有限体口门扩展 黄河因泥沙淤积严重而导致河床不断抬升并形成“悬河”，也使得堤防工程难度加大历 控导工程，黄河堤防的抗洪能力得到提高．但是，在高水位浸泡下仍会有发生溃决的可能”】． 历史上黄河口门堵复多在汛后进行，但是如今黄河下游两岸平原目前已成为人口、社会和经济 发展非常集中的地区，决口事件势必会对国计民生造成巨大的影响．因此，需要充分考虑黄河 下游可能出现的险情，以便增强在出险后的应变能力，为堵口技术的改进和完善提供科学依据 黄河堤防多由不耐冲刷的松散沙土为主建造，决口后口门会迅速扩展，口门附近很快出 现深坑，因而不利于堵日工作的进行溃堤过程和低水头土石坝的溃坝过程在许多方面类似， 其水流特征都可简化为非恒定含间断的泼水问题．但其不同之处也非常明显．溃坝水流就是 河道的主流．溃堤水流则是河道水流的一个侧支流．溃坝问题多关心坝体溃决后的流量过程、 洪水演进速度和淹没范围，溃堤问题则更关注如何在最短的时间内将决口堵复．这就需要对 决口口门区域的流场、地形变化及口门展宽速度等有比较深入的把握 针对堤坝决口和洪水演进的数值模拟研究已有不少成果例如．美国国家气象局的溃坝 到了广泛应用i31近年来，基于无结构网格的有限体积法已被用于溃坝数值模拟研究”…，该 方法综合了有限差分和有限元方法的优点，既可通过灵活的生成网格系统来满足处理复杂边 200i—iI一05临稿，2002—06ti0收修改稿 ·国家自然科学基金和水利部联合资助项目(批准号 万方数据 第5期 梁 林等黄河溃堤过程数学模型及其模拟h鞋、 界和地形问题的要求，义可避免有限元法计算量和守恒误差大的缺点．并且能自动捕捉阳J断 和消除虚假振动．然而．将溃坝过程模拟方法用于溃堤过程模拟时．日前还存在一些问题例 如．『垌格系统需要能够协调河道的大尺度网格和堤防的小尺度网格，这意味着需要解决局郭 加密的问题溃堤r1门展宽过程则需要在口门处动态地改变边界、增加同格．一般的有结构网 格不适于处理大小尺度变化的问题，而无结构网格在动态增加网格方面也存在困难． 本文基于浅水方程、泥沙输移方程、河床变形方程和口门展宽方程提出了一个综合的数 学模型，在四叉树网格系统的基础上采用Godunov型有限体积法求解双曲刊质量和动量守恒 方程当在复杂地形条件下使用某些Godunov型格式时，由于底坡的影响会产生守恒误差…“ 为了避免这种误差．可以对浅水方程某些项进行适当调整以平衡源项和通量梯度项的影响”J 1溃堤数学模型 1．i二维水动力学模型 一般河流的平面尺度远大于水深，因此河道水流适合采用沿深度平均的二维水动力学模 型来描述．在假设静压分布条件下．浅水方程可表示为… 堕+堂盟+塑：o． … 3t dr 3y