降雨微分响应修正在港口湾水库洪水预报中的适用性分析.pdf

水 力 发 电第 卷第 期 年 月 .降雨微分响应修正在港口湾水库洪水预报中的适用性分析王玉丽 李京兵 张 峰 张锦堂(安徽省水文局 安徽 合肥 安徽水利水电职业技术学院资源与环境工程学院 安徽 合肥)摘 要:港口湾水库位于长江南岸一级支流水阳江上游支流西津河上 对水阳江流域防洪起着重要作用 降雨微分响应修正是基于系统和微分角度的模型误差修正方法 将该方法应用于港口湾水库流域 对该流域 场历史洪水分别进行雨量单站微分响应与流域平均微分响应 种模式的修正计算 计算结果表明 种模式的降雨微分响应修正方法在不损失预见期的情况下均能够有效地提高港口湾水库洪水预报的精度 值得进一步研究与推广关键词:洪水预报 降雨微分响应 适用性分析 修正 港口湾水库 .中图分类号:文献标识码:文章编号:()收稿日期:基金项目:安徽省自然科学基金“水科学”联合基金项目()作者简介:王玉丽()女 安徽阜阳人 硕士 主要从事水文预报与水情信息管理.水文模型的研究对象主要是时不变的离线系统引起洪水预报误差的因素有很多 主要分为预报模型初始状态值误差、模型输入误差、模型结构误差以及模型参数非优化等 以上误差均会影响预报结果的精度 有效的误差修正技术会减少模型在预报时产生的误差 进一步提高预报精度 系统响应修正方法是一种结构简单、具有物理成因机理且不增加模型参数的实时洪水修正新方法 该方法的有效性、收敛性均已得到了证明降雨资料作为洪水预报模型最重要的输入参数雨量资料的误差是洪水预报模型误差的重要来源其精度直接影响着预报结果的好坏 本文在文献、研究的基础上 将降雨微分响应的误差修正方法应用于港口湾水库流域 从修正降雨量角度出发 将模型计算概化为一个系统 通过时段降雨引起的模型微分响应来修正降雨 进而修正出口断面的流量过程 本文将雨量单站微分响应修正与流域平均微分响应修正应用于港口湾水库流域 并比较 种模式修正效果的差异 探讨分析降雨微分响应修正方法在港口湾水库流域的适用性第 卷第 期王玉丽等:降雨微分响应修正在港口湾水库洪水预报中的适用性分析 .1流域简介港口湾水库是水阳江上游控制性工程 是一座以防洪为主 结合发电、灌溉、城市供水、水产养殖、旅游等综合利用的大()型水库 港口湾水库流域系皖南山区湿润地区 四季气候温和 雨量丰沛 多年平均降水量为 降水年际变化大年内分配不均 多年平均气温为 多年平均蒸发量为 其多年平均年径流量为 亿 流域植被覆盖率为 由于流域坡降大、汇流快 所以洪水过程陡涨陡落 从流域的气象条件、下垫面条件和流量过程分析可知 该流域降雨径流关系具有蓄满产流的特点 可以按蓄满产流的理论和方法建立产流量预报方案 故选用降雨微分响应修正与三水源新安江模型结合进行产汇流计算2降雨误差微分响应估计方法介绍基于“水文系统响应理论”将模型计算与实测之间的流量误差归结为降雨误差 只考虑模型输入变量中降雨因素的变化而引起的系统响应时 利用流域出口断面的流量信息 通过建立流量增量与降雨增量间的微分响应关系 计算时段面平均降雨引起的模型微分响应来对面平均降雨误差进行反演估计 将估计得到的面平均降雨过程重新代入到模型计 算 进 而 达 到 修 正 出 口 断 面 流 量 过 程 的目的()()()()式中 为时间 为流量 为降雨为降雨对模型计算系统的微分响应假设样本长度为 表示降雨系列 流域出口断面的实测流量系列为 代入新安江模型计算得 计算流量系列为 ()()令流量模拟误差为()()()流量观测随机误差为 则流量模拟误差的矩阵表达形式为()式中 为需要通过微分响应估计求解面平均雨量误差 假设需要估计的面平均雨量长度为 则 为模型计算系统的微分响应结合式()和式()可知 的表达式为 ()()()()()()()()()()()()()式中 为列向量 表示第 时刻的面平均降雨对应模型计算系统的微分响应 式()中的每一项均可以用差分近似求解()()()()由于降雨量系列是随时间实时变化的 并且由式()可知 微分响应系列受其它时段的降雨量值影响 因此 降雨量 所对应的微分响应系列也是动态变化的根据降雨径流形成过程可知 时段出口断面的流量不会受到 时段降雨的影响 因此 时段的降雨 产生的微分响应不会对 时段之前的流量值产生任何影响 即()()()根据最小二乘估计原理 可得到第 步的降雨误差估计()()()因为 为白噪声向量 符合零均值分布的特点所以()()因此 降雨误差微分响应最小二乘估计值为()()则新的降雨系列为水 力 发 电 年 月 .()其中流域平均微分响应与雨量单站微分响应的区别在于:计算微分响应 时 是否将流域内的所有雨量站当做一个整体或是个体考虑 若为前者则为流域平均微分响应 若为后者 即为雨量单站微分响应 计算步骤如下:()是原始降雨系列 将其代入到新安江模型中计算 得到出口断面流量()计算流量模拟偏差系列()()()计算微分响应()确定需要估计误差的降雨量时段数为 对每个时段降雨()()依次增加一个单位的降雨量 其他时段降雨量不变 得到新降雨系列 将其代入模型计算出新流量过程 减去原降雨量系列计算的流量过程 得到时段降雨()对应的微分响应列()由此求得微分响应()()根据式()计算得到第 步微分响应降雨误差估计为()()()()()()根据式()计算第 步的降雨估计 即()()()()将新的降雨估计代入新安江模型计算 得到第 步的模拟流量系列()判断第 步的降雨估计是否为最优值 若是 则结束寻优 否则 转步骤()继续循环3港口湾水库流域应用 评价指标选取港口湾水库流域 年 年的 场次洪水资料 通过新安江模型计算与降雨微分响应修正 本次应用检验具体通过以下 项指标来衡量修正效果()相对径流深误差 为 ()式中 为实测径流深误差 计算径流深误差()相对洪峰误差 为 ()式中 为实测洪峰流量 为计算洪峰流量()确定性系数 为 ()式中 为实测流量系列 为实测流量均值为计算流量系列()洪峰滞时 为 ()式中 为实测洪水峰现时间 为计算洪水峰现时间()修正有效性系数 为 ()式中 为修正后的流量系列 表示修正的相对有效性 其值越大 说明修正方法的效果越明显其最大值是 计算结果本文选用港口湾水库流域 年 年间 场洪水资料 分别进行新安江模型计算以及雨量单站降雨微分响应修正与流域平均降雨微分响应修正其中新安江模型参数通过人机交互率定得到 且所有参数在计算和修正过程中保持不变 流域内雨量站主要有合庄、和阳、黄土坎、俞村、胡乐司、刘村、方塘、港口湾等 个雨量站 流域水系及雨量站点分布 如图 所示图 港口湾水库流域水系及雨量站分布港口湾流域 场洪水的新安江模型计算结果及单雨量站降雨微分响应修正与流域平均降雨微分响应修正结果 如表 所示基于以上 种模式降雨微分修正方法的计算结果 综合统计如表 所示由于港口湾水库的出口断面入流量是由实测水位、出库流量和库容曲线反推而得到的 大多呈现第 卷第 期王玉丽等:降雨微分响应修正在港口湾水库洪水预报中的适用性分析 .表 港口湾水库流域计算及降雨微分响应修正计算洪号新安江模型/雨量单站微分响应修正/流域平均微分响应修正/注:为修正计算的循环次数表 港口湾水库流域不同计算方法的评价指标统计计算方法/新安江模型 流域平均微分响应修正 雨量单站微分响应修正 锯齿波动过程 甚至出现降雨与流量过程不对应的现象 存在着异常误差 致使新安江模型计算结果精度评定整体不高 由表 可知 其中新安江模型洪量相对误差绝对值的平均值为 洪峰相对误差绝对值的平均值为 洪峰滞时绝对值的平均值为 确定性系数为 经过雨量单站微分响应修正后 洪量相对误差绝对值的平均值为 洪峰相对误差绝对值的平均值为 洪峰滞时绝对值的平均值为 确定性系数为 经过流域平均微分响应修正后 洪量相对误差绝对值的平均值为 洪峰相对误差绝对值的平均值为 洪峰滞时绝对值的平均值为 确定性系数为 由表、可知 种模式的降雨微分响应修正方法对于改进模型计算精度均是有效的 雨量单站微分响应修正的平均循环次数为 流域平均微分响应修正的循环次数为 雨量单站微分响应修正结果的精度较流域平均微分响应修正结果的精度高、效率快 经过降雨微分响应修正后的确定性系数整体大于模型计算的结果 其中有 场洪水联合修正后的确定系数小于模型的结果由表 可知 有 场次洪水的修正有效系数小于 修正无效 其余 场次修正均有效 例如由表 中 场次洪水计算及修正结果可知该场次洪量、洪峰的计算结果与实测偏差大 相对误差均大于 经过流域平均微分响应修正和雨量单站微分响应修正后 确定性系数由 分别降到、修正结果偏差较大 修水 力 发 电 年 月 .正无效 以港口湾水库流域、场次洪水修正前后流量过程为例 如图 所示 由图 可知 场次反推的入库流量过程无锯齿现象 而 场次反推的入库流量过程中锯齿现象较为明显 场次的修正有效系数为 场次的修正有效系数为 修正均有效 但 场次的修正效果较 场次明显 综上可知 实测资料的准确性及新安江模型计算结果精度对于降雨微分响应修正的效果存在一定的影响 但从整体而言 降雨微分响应修正适用于港口湾水库流域 并且效果明显图 港口湾水库实际修正效果4结语()本文基于系统响应理论 单独考虑降雨误差所采用的降雨微分响应修正法 物理意义明确在不损预见期的情况下 能够有效提升模型预报精度 由计算结果表明 降雨微分响应在港口湾水库流域洪水预报中取得了较好的修正效果 其中洪峰、洪量相对误差及洪峰滞时得到不同幅度的减小 修正后精度提高 且 种模式的降雨微分响应修正的循环次数均较为合理()对比分析雨量单站微分响应与流域平均微分响应修正方法 在港口湾水库流域洪水预报中的应用效果有差别 但总体而言可行有效 可以应用于港口湾水库的洪水预报中 并值得进一步研究与推广参考文献:包为民.水文预报.北京:中国水利水电出版社.司伟 包为民 瞿思敏.洪水预报产流误差的动态系统响应曲线修正方法.水科学进展 ():.包为民 司伟 沈国华 等.基于单位线反演的产流误差修正.水科学进展 ():.包为民 瞿思敏 李清生 等.遥测系统降雨观测误差估计方法研究.水利学报():.瞿思敏 包为民.实时洪水预报综合修正方法初探.水科学进展 ():.胡余忠 朱先武.港口湾水库入库洪水预报模型研究.水文():.李京兵.安徽大洪水.南京:河海大学出版社.司伟 包为民 瞿思敏 等.基于面平均雨量误差修正的实时洪水预报修正方法.湖泊科学 ():.杨姗姗 包为民 杨小强 等.微分响应在降雨误差修正中的应用.中国农村水利水电():.(责任编辑 安莉莉)(上接第 页)刘艳艳.内蒙古低碳经济发展研究.北京:中央民族大学.温鑫 姚文莹 何圣川 等.双碳目标下新能源消纳能力提升方法研究.机电信息():.刘同良.中国可再生能源产业区域布局战略研究.武汉:武汉大学.姚德全 陈卫中 郑宁敏 等.新能源跨省消纳互济研究.中国电力企业管理():.李铭 陆文升 梁建斌 等.面向新能源消纳的电网规划方法研究.光源与照明():.严灏.新能源电力消纳与我国新型电力系统建设.工程经济 ():.张鑫 刘飞 王世斌 等.提升特高压直流送端电网新能源消纳水平措施 .电力系统及其自动 化学报 ():.周泽昕 于溯 李勇 等.新能源经柔直送出场景下功角变化导致比相式距离保护不正确动作机制分析.中国电机工程学报 ():.(责任编辑 焦雪梅)