实时校正中的旁侧入流反演方法 Inverse analysis of lateral inflow in real-time correction.pdf

第加卷第i期 水科学进展 V01．20．No．1 2009年1月 ADVANCEsINWATERscIENCE Jan．．2009 实时校正中的旁侧入流反演方法 吴晓玲，王船海，向小华 (河海大学水文水资源与水利工程科学国家蘑点实验室，江苏南京210098) 摘要：在水文与水动力学相结合的预报模型中，对于水文模型计算有误差，造成水动力学区间旁侧入流与实际不吻 合情况，依据水量平衡原理，对旁侧人流进行反演校正。采用四点线性隐格式求解圣维南方程组，结合卡尔曼滤波 技术建立以水位、流量作为综合状态量的卡尔曼滤波状态空间方程。采用反问题的思想把校正后的水位、流量不平 衡量反馈到旁侧人流中，修正原有的区间来水量，进而修正预报期的旁侧人流。选择长江干流寸滩至万县河段，以 及长江下游大通至镇江段进行演算，分别采用三峡水库135m蓄水阶段的水文资料和1998年水文资料进行6h洪水 预报，结果表明旁侧人流反演校正方法是可行的。 关键词：实时校正；反问题；卡尔曼滤波；旁侧人流 中图分类号：P338．1 文献标识码：A 文章编号：|001-6791(200910boos2-06 水文与水动力学相结合的洪水预报技术受到越来越多的广泛关注…，由于在构建水文预报模型过程中会有 很多误差导致模型预报不准确㈨2。因此对于洪水预报模型来说，误差校正是必不可少的，它使模型预报沿着正 确的轨迹进行引3。水力计算的误差主要受河道糙率，旁侧人流，上、下游边界条件，对于非恒定流计算还包括 初始条件【4J等因素影响。一般情况下，初始条件随时间推移对流态的影响减弱；用充分信任的实测过程作为 上、下游边界条件；通过数场洪水进行参数率定，尽量降低边界条件和糙率不准确对模型计算的影响L5J。此时 的水动力学模型计算结果与控制站的实测过程相比，仍存在水位、流量的非随机性有色噪声，为了白化该有色 噪声，认为此水位、流量的不平衡量是由水文、水动力学模型相关联时，水文模型出流计算不准确导致的。针 对该问题，采取反问题【6,7]的思路，把控制点处的水位、流量的不平衡量，通过卡尔曼滤波技术反馈到水文模 型计算，从而得到对旁侧人流的修正，达到校正水文、水动力学耦合预报模型的效果。 本文以离散并线性化的圣维南方程组为基础，建立以断面水位、流量为状态量的卡尔曼滤波方程，运用卡 尔曼滤波反问题的分析方法把水位、流量不平衡量反推到旁侧人流中，对其进行修正。从而为准确计算区间旁 侧人流提供了一条新途径，并选择两种典型流域验证本文方法的普遍性。 1 计算方法 1．1 离散线性化圣维南方程组 描述河道水流运动的圣维南方程组为 『曰箐+篑=叮 I篑+袅(譬)州瓦3ZⅧ：廿：。 。’ 式中口为单位河长I-分布的旁侧人流量；Q、A、B、Z分别为河道断面流量、过水面积、河宽和水位：n 收稿日期：2007．1I．16 基金项目：“十一五”国家科技支撑计划课题(2006BAc015Ⅸ】2) 作者简介：吴晓玲(1981一)，女(满族)，黑龙江牡丹江人，博士研究生，主要从事水动力学洪水预报方面的研究。 E-mail：freebir7237@hhu．edu．cn 万方数据 第l期 吴晓玲。等：实时校正中的旁侧人流反演方法 53 为河道糙率；R为水力半径；口为动量校正系数。 如下差分方程： f饼：；一谚“+甲cj7川k+I+喇“=毋+q(破-+才) ，¨ 怕7 【喇“+州t-flk川+I+垆FTk川+,一嘲“=q(醵l+钟)一‘ 式中 G：酉Bj+{Ax／；2：竺弓全兰一100-(Q，+。一钟)；弓：会惫一(口Ⅱ妒+_