水动力学模型卡尔曼滤波实时校正技术.pdfVIP

水 利 学 报 2005 年6 月 SHUILI XUEBAO 第36 卷 第6 期 文章编号: 0559-9350(2005)06-0687-07 水动力学模型卡尔曼滤波实时校正技术 葛守西，程海云，李玉荣 (长江水利委员会 水文局，湖北 武汉 430010) 摘要：本文采用卡尔曼滤波技术实现圣维南方程实时校正以提高水动力学模型应用于洪水预报的精度水平。首先 利用有限差分方法将圣维南方程进行离散化并线性化，选择不同的状态变量组建了三种卡尔曼滤波状态空间方程， 对预报进行实时校正。选择长江干流清溪场至万县河段进行试算，证实此方法对提高预报精度的效果显著。 关键词：水动力学模型；卡尔曼滤波；实时校正；状态向量 中图分类号：TV131.2 文献标识码：A 水动力学模型应用于洪水预报，主要解决水文学方法无法解决或解决得不甚理想的一些课题。例如： 复杂水情、分汊水系及分洪、溃口等预报对象，以及河道断面特征发生了明显变化后的水力要素计算；尤 其当需要计算河道内任何断面任何时刻的水位、流量等要素时，水动力学模型通常为首选方法。多年来， 国内外学者对水动力学模型应用于洪水预报开展了大量研究，然而，由于其预报精度不甚理想，长期以来 未能在实时作业预报中得以广泛应用。如何通过实时校正来提高预报精度成为一个重要的课题。实时校正 技术首推卡尔曼滤波，国内外在流域洪水预报方面应用卡尔曼滤波的研究不少，而将其与水动力学模型联 合应用尚属空白，本文对此课题作初步探讨。 1 圣维南方程的离散化和线性化 扩展的圣维南方程由连续方程和动量方程组成 ∂Z 1 ∂Q 1 + Lq (1) ∂t B ∂x B  2  ∂ ∂Q ∂ Q   Z  ∂t +∂x  A +gA∂x +Sf  0 (2) 2 n Q Q Q Q 式中：S 为水力比降，S = 或S = ,其中R为水力半径，n为曼宁阻力系数，K为流量模数。 f f f 2 4 / 3 2 A R K 采用Preissmann四点隐式差分格式对连续方程和动量方程进行时间(i)和空间(j)的离散化。并令 i+1 i f =f + Δf(f代表因变量)，对差分方程进行幂级数展开并线性化，则对于j河段的连续方程和动量方程分 别简化为 收稿日期：2004-08-22 基金项目：国家自然科学基金重大项目 作者简介：葛守西(1941-)，男，浙江嘉兴人，教授级高级工程师，从事水文预报研究。 E_mail:gesx@ 1 水 利 学 报 2005 年6 月 SHUILI XUEBAO 第36 卷 第6 期