漳卫南运河洪水资源实时调配数值模型研究.pdfVIP

《水资源研 究》第 32卷 第 1期 总第 118期 2011年 3月 漳卫南运河洪水资源实时调配数值模型研究 雷 四 华，吴 永 祥 ，王 高 旭 ，汤 绍 青 南京水利科学研 究院，水文水资源与水利工程科学国家重点实验室，江苏 南京 210029 摘 要 ：研 究洪水资源在汛期 内如何实时配置 ，有利于洪水资源的蓄存 ，通过洪水资源实时调配数字模 型计算 ， 可模拟各种洪水资源利用方案 的合理可行性 。系统介绍 了漳卫南运河洪水资源数字模型所采用的数学模型方 法及模型河 网各类节点处理方法 ，并给 出了模 型验证及应用结果 。结果表 明：该模型适合计算 多种情景洪水及 引水利用过程 ，计算效率 、计算精度满足要求。 关键词 ：洪水资源 ；实时调配模型 ；计算水力学 ；漳卫南运河 漳卫南运河水 系是海河 流域五 大水 系之一 ，位 于太 岳山 以 1．2 方程离散及数值 方法 东，滏 阳河 、子牙河 以南 ，黄河 、马颊河 以北，由漳河 、卫河 、卫运 对流体运动方程采用简化 的四点线性隐格式方法，可得数 河 、漳卫新河和南运河组成 ，流经 山西 、河南 、河北 、山东 4省及 值差分方程 J。水位边界条件计算递推公式见下式 3 。在树 天津市人渤海 ，流域面积 37700km 。流域地处温带半 干旱 、半 状 河网中，对 内节 点上边界均为水位 过程 的下游 河段递 推公式 湿润季风气候 区，多年平均年 降雨量为 549IllIll，水资源总量 系数可 由上游河段得 出，如河道 1、河道 2为上游河段 ，河道 3 102．1亿 m ，全流域人均 212m ，低于海河流域 350m 的人均 为下游河段 ，河道 1的第 Ⅳ个断面 、河道2的第 个 断面 、河道3 值 ，不足全 国的 1／9，是我 国水资源严重短缺 的地 区之一 ，而且 的第 +1个 断面为汇合节点断面。 时空分布不均 ，水生态环境状况恶化程度高… 。针对 日趋严重 河 道 1的递 推公式为 ： 的水资源供需矛 盾 ，如果能尽可 能利用漳卫南运河洪水 资源 ，将 z P 一 Q 3 洪水资源化 ，对缓解海河流域水资源供需矛盾 ，促进社会经济持 河道 2的递推公式为 ： 续发展具有重要 意义。 z： P 一 Q： 4 相容关系方程有：z ：z z：+和 Q Q Q+ l 数学模型 将河道 1、河道 2的递推关系代人相容关系的流量方程得 ： 1．1 流体运动方程 Q 甓+甓～c+甓， 河 网水位流量过程求解的基本方程采用 明渠非恒定流连续 再用相容关系方程 的水位相等条件整理上式得 ：Q+ 方 程 见下式 1 和运动方程 见下式 2 联 合求解 。 P + 一 + +1。 警a+Ot qL 1 其中系数关系为：P 甓+和+t 亩+亩 式 中， 可理解为流量随流程 的变化率 ；oA为断面面积随时间 将河道 3的递推关系转为与河道 1、河道 2同式得 ：z+ 0 P：+一 +，Q+。。 变化率 ；qL为单位长度上的旁侧汇流流量 ，汇人流量取正号 ，出 流流量 为负号。 其 中系数关 系为 ：P+ + 一 … + ～ 。 at + 杀 +去+譬 0 z 1．3 堰 闸过流计算 式 中，流量模式 K AC ，C为谢才系数 ，R为水力半径 ，A为过 对于 自由出流 孑L流 、堰流 有计算公式 ： mbh 流 面 积 5 收稿 日期 ：2009—09—25 资助项 目：“十一五 ”国家科技 支撑计 划课题 2006BAB14B07 研 究 内容 作者简介 ：雷四华，男，高级工程师，主要从事水文水资源、水利信息化研究。 15 式 中，m为综合流量系数 ；6为 闸孔 净宽 ；h为计算过 流水深 闸 对 内节点根据河 网关系分别处理 ，漳河 、卫河与卫运河汇交 前水位 一堰高程 ，Z—d ；孔流 h o，。为闸前水头或总水头 n 节点 ，由递推公式及相容关 系处理 ；卫运河与岔河 、减河相交节 Ⅱ 点 由水 闸控制 ，采用水 闸调度 公式处理 ；岔河 、减河与漳卫新河 ho+— ；堰流 h h0。 g 汇交节点 ，由递推公式及相容关系处理 。 由于推导公式 Z P一Q 水位为上边界 的P、值 已知 ，联 3 洪水资源利用原则 合求解可得 Q 0函数 。 ，丫r】、 洪水资源利用是指在不成灾 的情况下 ，尽量利用水库、拦河 利用迭代公式 Q Q；一 ，可求的 Z、Q。 ＼Y i， 闸坝 、自然洼地 、人工湖泊 、地下水库等蓄水工程拦 蓄洪水 ；以及 淹没 出流 闸后水位大于堰底高程 有计算公式 p Obh · 延长洪水在河道 、蓄滞洪区等滞 留的时间；恢复河流及湖泊 、洼 ／