淮河中游大型人工新河分洪效果分析

I【科技论坛
淮河中游大型人工新河分洪效果分析
郁玉锁赵凯
(安徽省淮河河道管理局蚌埠233000中水淮河规划设计研究有限公司合肥230000)
【摘要】茨淮新河及怀洪新河是1970年后在淮北平原开挖的两条大型人工河道，工程主要目的是分泄中游洪水，减轻淮河干流的防洪压力。

本文借助一、二维耦合水动力数学模型，模拟了淮河中游干支流河网的洪水演进情况，研究了茨淮新河以及怀洪新河2003年运用对淮河干流洪水传播的影响，从水力学的角度评估了典型洪水条件下两条人工新河的实际分洪效益，并与水文学方法计算的成果进行了对比。

【关键词】淮河中游大型人工河道分洪效果水动力数学模型
1概述
茨淮新河及怀洪新河是横亘在淮北平原上的两条大型 人工河道，工程始建于20世纪70年代初，分别于1991年 和2002年建成。

其中茨淮新河全长134km，最大分洪流量 2000m3/S，主要作用是遇沙颍河或淮河发生特大洪水时，分 泄沙颍河洪水，尽量减少沙颍河来水对淮河干流的影响;怀 洪新河是茨淮新河的接力河道，全长121km，最大分洪流量 2000m3/,，其主要作用是分泄淮河干流洪水人洪泽湖，尽量 减少上游来水对蚌埠以下河段的防洪压力。

茨淮新河及怀洪新河建成后，历经多次分洪，效益显著。

为充分发挥两条人工新河的作用，进一步完善工程调度运行 方案，本文利用淮河中游一、二维耦合水动力数学模型，分析 模拟人工新河运用与干流水位的响应关系，系统评估典型年 洪水条件下茨淮新河及怀洪新河实际运用的效果。

2数学模型简介
淮河中游地势平缓、河网交错、行蓄洪区众多，相互之 间的水力联系复杂，开展河道洪水演进的模拟研究是解析 淮河干流洪水变化特性、合理安排治淮工程布局、实现洪水 科学调度的关键所在。

淮河中游的洪水演进方法可分为水 文学和水力学两种，水文学方法通过马斯京根方法演算各 断面的流量，然后用扩散波非线性水位法推求水面线;水力 学方法的核心是求解圣维南方程组或浅水方程组，相对而 言计算精度较高，但需要详细的河道地形资料。

虞邦义等基于对淮河中游大型河网水流运动的分析，利用MIKE水流模拟系统建立了淮河中游水动力数学模 型。

该模型采用一、二维嵌套的模型架构，对淮河中游河道、湖泊、行蓄洪区、水利工程枢纽等河网典型要素进行分别概 化，将蓄洪区视为零维调蓄节点，将干支流河道水流运动处 理成一维模型，将洪泽湖、沿淮行洪区、重要生产圩等宽浅 水域处理成二维模型，通过堰闸等水工建筑物过流计算实 现不同模型和不同区域的耦合。

目前，该模型已初步应用于 淮河中游河道整治及行洪区调整工程、沿淮重要区域洪水风险图绘制以及行蓄洪区优化调度等方面，为淮河进一步 治理和保护提供了重要的技术支持手段。

鉴于该模型能客观反映淮河中游洪水运动规律，精确模 拟各种调度方案对河网水动力所产生的影响，故本文采用其 作为研究淮河中游大型人工新河分洪效果的计算工具。

3实际分洪过程及效果分析
3.1实际分洪过程
在2003年6月下旬~7月中旬期间，茨淮新河茨河铺 闸共有3次运用过程，共分泄沙颍河洪水4.73亿m3,最大 分洪流量1580 m3/,。

具体分洪情况见表1。

2003年7月怀洪新河何巷闸共3次开闸分洪，总分洪水 量17.18亿m3,最大分洪流量1670m3/,。

具体分洪过程见表2。

3.2分洪效果分析
图1为茨淮新河分洪与不分洪时正阳关的水位过程。

从 图中可以看出，茨淮新河若不分洪，正阳关第一次的洪峰水位 将达26.78m，超保证水位0.38m，相应鲁台子洪峰流量为 8532m3/,;第二次洪峰将达26.87m，超保证水位0.47m，相应鲁 台子洪峰流量7889m3/,。

两次洪峰保证水位以上的洪水历时 将达137h。

茨淮新河分洪后，正阳关洪峰水位及鲁台子洪峰 流量均有所降低，其中第一次和第二次正阳关洪峰水位的降 幅较明显，分别为0.33m和0.17m，相应鲁台子洪峰流量分别 削减了 642 m3/,和270 m3/,。

上述计算表明，为保证淮干正阳 关以下防洪安全，在2003年洪水期实时开启茨河铺闸分泄沙 颍河洪水是十分必要的，效果也是明显的。

图2为怀洪新河分洪与不分洪时吴家渡的水位过程。

从图中可以看出，怀洪新河若不分洪，吴家渡第一次洪峰水 位将达22.52m，超保证水位(22.48m)0.04m，相应洪峰流量 9422m3/,，第二次洪峰水位22.37m，相应洪峰流量8688m3/,。

怀洪新河分洪后，吴家渡站洪峰水位将降低0.37~0.63m，洪峰 流量也将减少350~850m3/,，且三次洪峰的水位均在保证水位 以下。

上述计算表明，2003年怀洪新河的适时分洪，为缓解蚌 埠城市圈堤及淮北大堤的防汛压力起到了关键的作用。

■【科技论坛
注：茨淮新河对应正阳关站(水位）及鲁台子站（流量），怀 洪新河对应吴家渡站。

4水文学与水力学方法成果对比分析
《2003年淮河暴雨洪水》(以下简称文献1)采用水文学
方法,分析了茨淮新河和怀洪新河2003年分洪的效果。

表3 给出了两者主要成果的对比情况，由表3可知，水文学及水 动力模型在削减洪峰流量、降低洪峰水位以及缩短警戒水 位以上历时等方面还存在一定的差异，初步分析原因主要 有以下几点：
4.1洪水演算方法
在分析人工新河不分洪时干流的洪峰流量时，文献1 采用马斯京根法将茨河铺的实测分洪过程演算到鲁台子 站,并与原鲁台子实测流量叠加。

由于马斯京根法是在出流 与槽蓄量单一函数关系的假定下导出的，是运动波的差分 解，这种假定在河底比降较大的山区河道与实际基本符合， 有足够的精度。

而对于淮河中游平原河网地区，河底比降较 小,洪水演进受上游来流、下游水位以及闸坝调控的联合作 用,上述假定与实际情况有偏差。

此种情况较宜采用动力波
，即直接求解圣维南方程组。

4.2经验相关关系
在分析人工新河不分洪时干流的洪峰水位时，文献1 普遍采用了水位流量的相关关系及上下游水位的相关关 系。

受洪水涨落及回水顶托的影响，吴家渡的水位流量呈逆 时针的绳套曲线，特别是遇到复式洪峰时，水位流量的绳套 关系更加复杂。

此外，正阳关站与鲁台子站的水位相关关系 受鲁台子流量及其南北区间洪水组成的影响，十分复杂。

简 单的经验在相关处理上述复杂的情况是有难度的。

4.3分洪之间的互相影响
在评估人工新河的分洪效果时，分洪过程效果应作为 一个整体进行分析。

这是由于一次分洪在削减本次洪峰流 量，降低本次洪峰水位的同时，也改变了下一次分洪时河道 的起涨水位等条件，因此，每次分洪的效果互相影响，很难 截然分开。

以怀洪新河2003年分洪效果为例，吴家渡第二 次洪峰的削峰效果是何巷闸第一次和第二次分洪的综合作 用，而吴家渡第三次洪峰的削峰效果则包含何巷闸三次分 洪的综合作用。

5结论与讨论
本文应用淮河中游一、二维耦合水动力数学模型重演
了 2003年典型洪水，研究了茨淮新河、怀洪新河的分洪效 果，并与水文学方法的成果进行比较，初步讨论了两者差异 的成因。

研究表明，水力学模型能够更为精细地模拟淮河中 游的水流运动，能有效地分析人工新河分洪的效果。

今后可 以将水文学及水力学方法相结合，补充分析闸坝调度对区 域汇流条件的影响，以进一步提高调度效果评估的精度■
：
7/8
7/18
7/28
8/7
B 期
图1 2003年正阳关水位过程线图
日期
图2 2003年吴家渡水位过程线图
表1 2003年6月下旬& 7月中旬茨河铺闸分洪情况统计表
表2 2003年7月上旬（7月下旬何巷闸分洪情况统计表
分洪起讫时间
最大
分洪流量
分洪水量
分洪期间淮河上、下游
控制站水情
淮河
田家庵
淮河吴家渡
洪峰流量洪峰水位洪峰流量
(月一日时:分）(m3/s )(K m 3)(m3/s )(m )(m3/s )7—4 10:00~7—7 20:001590 3.40524.2621.9486207-9 12:0〇3 7-18 14:30
1670
8.48324.0721.7479207—22 02:00 4 7—27 21:301510
5.29423.20
21.53
7430
合计
17.182
表3
结果对比表
分洪起讫时间
最大
分洪流量
分洪水量
分洪期间下游控制站水情
沙颍河阜阳闸淮河正阳关
洪峰流量洪峰水位洪峰流量
(月一日时:分）(m3/s )(K m 3)(m3/s )(m )(m3/s )7—2 20:3(5 7—6 09:451240 2.11248026.4578907—12 21:5767—17 08:151580 2.104关闭26.70
7620
7—19 13:15 7 7—20 17:18
607
0.511关闭
持续回落
合计
4.725
人工
河道年份
运用次序降低干流洪峰 水位（m )削减干流洪峰
流量(m3/s )文献1本文文献1本文茨淮新河2003
第一次0.190.33640642第二次0.210.17390270怀洪新河
2003
第一次
0.40.581440802第二次0.450.631580788第三次
0.36
0.37
920
370
(I I
o o o o o o o o o o o o o o o o o
o
o
50505050 o 5443322115
5000500000500000500000500000
5
5
2.1.1.0.
2 2 2
2
(3
匀
爷
M M n l K。