河道崩岸、护岸、清淤、生态需水量等相关论文 河道糙率反问题在实时洪水预报中的应用.pdfVIP

第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 河道糙率反问题在实时洪水预报中的应用 吴晓玲，王船海 (河海大学，水文水资源与水利工程科学国家重点实验室，南京 210098) 摘要： 针对洪水预报中糙率参数的实时修正问题，采用卡尔曼滤波技术，实现动量方程 实时校正糙率参数，以提高水动力学模型的洪水预报精度。对复杂地形所引起的局部动量损失 进行处理，使反演的糙率参数更显客观。 选择长江干流寸滩至万县河段进行试验，应用2004 年该河段洪水资料作验证。结果表明， 反演的糙率参数对提高实时洪水预报精度有效。 关键词： 实时洪水预报；卡尔曼滤波；糙率参数；反问题 1 前言 对于水动力学计算中糙率参数的选择问题，国内外的学者倾注了大量的研究精力。杨世孝 等[1]依据反问题方法研究河网的非恒定流计算，利用各时段各河段的糙率平均值动态调整河段 的平均糙率；金忠青等[2]根据河道规模及水力特性，把河道划分为若干等级，采用最优化中的 复合形法对不同等级的河道糙率分别求解；董文军等[3]利用最小二乘逼近思想建立曼宁糙率的 辨识模型；李光炽[4] 以曼宁公式及动量方程为基础，依据实测水位资料借助卡尔曼滤波技术反 推河道糙率。葛守西等[5, 6]针对应用卡尔曼滤波技术解决河段糙率问题，提出以水位、流量的 涨差作为观测量，建立以糙率为状态量的线性滤波器方程，通过实测水位、流量的变化对糙率 进行反推计算。 Fread 等人[7]通过流量分级，分段拟和平均流量与糙率的关系函数，运用牛顿־拉斐森算法 得到当水位误差函数取最小时糙率与平均流量的关系；Khatibi 等人[8, 9]应用最小二乘法，依据 实测水位解决潮汐河道糙率的最优估计问题，通过比较三种以水位作为“验证项” 的最小二乘构 造，得出相对高潮水位而言，低潮位对糙率的反推估计更加有效的结论。 受上述文献的启发，本研究采用四点线性隐格式[10]求解圣维南方程组，分离动量方程中 的水位及糙率项，依据卡尔曼滤波量测方程形式，建立水位与糙率的对应关系，运用滤波器的 最优估计原理，反推河段糙率；同时针对河底高程变化剧烈的河段考虑局部动量损失的影响， 把反推得到的糙率应用到实时洪水预报中，通过实时动态调整河道糙率参数提高洪水预报精 度。 国家支撑计划项目(2006BAC05B02)和国家自然科学基金资助项目资助. - 902 - 第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 2 模型的建立 2.1 河道糙率反问题-离散连续方程和动量方程 描述河道水流运动的圣维南方程组为 ∂ ∂  Z Q B + q  ∂ ∂  t x Q αQ2 Z u Q ∂ ∂ ∂ 2  + ( ) +gA +gn 4/3 0 (1) t x A x R ∂ ∂ ∂ f Q Z t 上边界条件：up ( , , ) 0  下边界条件 f Q Z t  ：down ( , , ) 0 式中：q 为旁侧入流；Q、A 、B 、Z 为河道断面流量、过水面积、河宽和水位；n 为河道糙率； R 为水力半径；α为动量校正系数。设河道被（L+1 ）个断面分成L 个子河段，对任一子河段 （j ，j +1 ）应用四点线性隐格式离散式(1)，得到如下差分方程，  k +1 k +1 k +1 k +1 k