青岛市大沽河流域防汛信息系统洪水预报数学模型的研究.pdfVIP

青岛市大沽河流域防汛信息系统 洪水预报数学模型研究 张彬 郑毅 方红卫 (清华大学水沙科学教育部重点实验室，北京100084) 摘要 防汛信息系统涉及降雨过程、产流规律、流域内各类水工建筑物的调蓄、洪水在河 道中的演进以及对可能造成灾害的评估等各个方面，其中流域内各类水工建筑物的调蓄模型 和洪水在河道中的演进模拟是整个防汛信息系统的核心，它的精度直接影响到整个系统的水 平。本文研究的重点是建立青岛市大沽河流域防汛信息系统中的两个模型，一是各类水库的 调蓄模型，该模型通过多年实测的水文资料和水库调度方式反映出降雨量和水库泄流的关系 曲线，为河道洪水的演进提供边界条件；二是增加内边界条件(橡胶坝)以后的洪水演进模型， 用橡胶坝的调度方式代替坝前的水位流量关系，再分段采用圣维南方程组和四点偏心格式进 行计算。用大沽河流域7个站的实测水位资料对模型进行了率定和验证，计算与实测吻合较 好，表明该研究所反映的降雨泄流关系模型和一维内边界洪水演进模型在大沽河防汛系统中 有较高的计算精度，可以推广至其他防汛信息系统的研究之中。 关键词 防汛信患系统：洪水预报，水库泄流曲线，圣维南方程，内边界 1前言 随着信息技术的快速发展，防汛信息系统的研制越来越受到防汛调度部门的重视，也为科 学高效快捷的制定正确的防洪决策提供了平台。 村庄，总人口179．3万人，耕地251．89万亩。大沽河主河道虽经多年治理，但下游段仍达不到 0年一遇的防洪标准。域内两座大型，五座中型水库也都存在着不同程度的病险因素。青岛 市水利局在1999年制定大沽河流域防汛规划时，计划对大沽河进行进一步整治，拓宽主河槽， 加高加厚两岸堤防，狭窄段予以退堤，在现有险工防护工程的基础上，对全部险工段进行防护 的安全行洪标准‘11。 但要完成不同阶段的防汛目标，必须编制正确的防汛预案，并在防汛决策的信息平台上予 以实现。目前青岛市这方面的工作正在完善中，本文的研究是在大沽河防汛规划[1]和大洁河 防汛工作手册[2]的基础之上，采用一维水力学数学模型的方法研究大沽河干流在不同洪水频 率下沿程流量与水位，主要工作为：①在CRESTID模型的基础上添加橡胶坝内边界子程序模 ·182· 块，避免因橡胶坝运行造成流量不恒定以及水位变化不连续导致的模型计算误差，可以正确计 算大沽河干流9座橡胶坝运行时的水位变化；②推求降雨——水库泄流工作曲线，为数学模型 计算提供边界条件。模型完成后验证数学模型中的糙率，在此基础上对大沽河干流规划前后 不同频率洪水进行计算分析，并提出各种防汛预案。 2非恒定流数学模型中的数值方法 连续方程： aQ Y xaAa_．4 8t=q 动量方程： adQ：dz署一adzz+g涨=象(2) (m)。 对于上述非线性偏微分方程组，采用四点偏心Preissmann格式差分方法求解[3--51。 数学模型中有关问题的处理 3．1大沽河主河道糙率的确定 糙率确定是否合理直接影响洪水在某一流量过程中的水面变化的计算精度，也是一个非 常复杂的问题。由于动床阻力的存在，边壁阻力在河流运动中很难量化，所以完全用理论方法 得到河道阻力难度较大，本模型采用的方法是根据不同流量级的水面线来率定河道糙率，再将 合理的糙率值用于河道计算，先统一将河道糙率定为0．035，在得到初步结果后将河道分段， 逐渐调整各河段糙率使结果与已知的现状资料(大沽河干流共设有7个监测点的水位资料)相 符合，并将此结果作为这段河道的糙率。 3．2各段河道的进口流量 大沽河自产芝水库以下共有9座大中型水库向大沽河主河道汇流，研究区域的人口流量 由产芝水库的下泄流量决定，各段河道的流量分别由各水库的汇人流量加上上游的来流量合 成。由于大沽河流域水库的运行方式与当地的降雨量直接相关，所以可以由每个水库的降雨 ——泄洪关系曲线得到各水库下泄的汇人流量。 首先可以通过资料中各水库的运行方式即每个降雨量对应着不同的泄流量，点绘一张降 雨——泄流关系点图，再用数据点拟合出以下方程作为水库的降雨—