水力学模型在淮河洪水预报系统中集成与应用.pdfVIP

x?褡麓蒸嚣霪慧冀熏冀雾夔鬻黧鞠豳蘑霾霾鬻阂 第四部分 水利信息学的新进展 水力学模型在淮河洪水预报系统’ 中的集成与应用 彭期冬 禹雪中 张建立骆辉煌 (中国水利水电科学研究院 北京 100038 pengqd@iwhr,com) 摘要 针对传统洪水预报方法中马斯京根法的局限，本文采用水力学模型模拟洪水在河道 中的演进过程。根据淮河洪水预报的需要，对糙率、边界条件、初始条件、行蓄洪区和分洪道等 进行了技术处理。将水力学模型进行了集成到淮河洪水预报系统中，达到操作简单，计算快捷， 结果准确的目标。该模型对淮河2002～2005年洪水模拟情况表明，模型模拟精度高，完全满足 洪水预报要求。本文工作对建立我国自主知识产权的洪水预报有重要意义。 关键词 洪水预报；水力学模型；集成；应用 1前言 传统洪水预报中，一般采用马斯京根流量演算法对河道汇流进行模拟。马斯京根法的基本原理是根 据河道人流和起始条件，通过逐时段求解的水量平衡方程和蓄泄方程，计算出流过程。马斯京根法计算 得到流量后，一般需要通过水位流量关系计算断面水位；然而，淮河由于洪水顶托、水利工程影响等因素 的作用，河道水位流量关系十分不稳定，由流量推求水位十分困难。而运用水力学模型来模拟河道中的 洪水演进过程，则可以克服马法的上述不足。但是一般的水力学模型要应用到洪水预报系统中，必须对 其中一些关键技术问题进行处理，同时还要通过一定方式集成到洪水预报系统的平台上，方便用户使用。 本文将一维水力学模型[1]集成到淮河洪水预报系统中[2]，在实际运用中取得了较好的结果。 2淮河洪水预报中主要技术问题的处理 水力学模型中最重要的参数是河道的糙率，此外模型计算需要以已知的流量、水位过程作为边界条 件，因此必须与水文预报模型相结合，在统一平台上进行集成，才能满足洪水预报的要求。 2．1河道糙率 ‘。 模型中将淮河河道糙率分为沿程阻力和局部阻力进行计算，沿程阻力主要反映河道边界粗糙程度对 水流运动产生的阻力，局部阻力主要反映河道断面形态变化产生的阻力。对于沿程阻力咒，，拟采用多组 次实测资料率定计算的方法，确定不同流量级下的沿程阻力系数。模型采用的综合糙率为反映沿程阻力 的糙率和局部阻力糙率的叠加。 2．2边界条件 水力学计算中，上边界一般采用流量过程，下边界一般采用水位过程或者水位流量关系作为定解 条件。 上边界预报期间的流量过程可以由水文模型得到。 下游边界如果采用实测水位过程作为边界条件，则洪水预报就没有预见期，并且这种方法在实际应 *本文受中国水利水电科学研究院科研专项(2004KY01)资助。 ·591· 水力学与水利信息学进展2007鞠黼黼黼麟嬲灞黼黼獭嬲黼黝黝燃㈣粼燃黝㈨㈣瓣％#。∥．， 用中也是不现实的。因此制定满足洪水预报的合理的下边界，在应用中显得尤为重要。模型中根据下边 界为湖泊、闸门或固定断面，分别采用不同的方法进行处理。 当下边界为湖泊时，一般都有较好的水位一容积关系，由湖泊的容积变化过程可得到相应的水位过 程，作为模型计算的下边界。模型中设置专门模块用于计算湖泊水位，由预报开始时刻的湖泊水位，得到 其初始容积，根据预报期间的出入湖流量过程及湖面降雨情况，计算其容积随时间的变化过程。 当下边界为闸门时，可建立闸门的流量一水位关系，或者闸门水位与特征站位的洪水要素之间的关 系，以数学函数的方式作为计算的边界条件。还可以人为设定闸门的控制方案，进行各种可能方案的计 算和分析。 当下边界为某个固定断面时，根据断面的水位一流量关系，确定下游边界条件，进行计算和分析。通 常情况下，一个断面的水位一流量关系不是固定的，有多个曲线，在实际应用时要根据当时的情况选用最 为适用的曲线关系。当下边界断面的水位一流量关系不理想，如变化过大时，还可以采用一种“水位+流 量”的下边界条件，即已有实测水位过程加上预报时段的流量过程，这个流量过程可以通过水文学的方法 得到。在实测水位