基于成因分析法和周期叠加法的中长期天然来水预测分析.pdf

水文水资源水利规划与设计 年第 期:/基于成因分析法和周期叠加法的中长期天然来水预测分析徐 峥 左茂盛 任 腾(国家电网公司东北分部 辽宁 沈阳 沈阳工业大学机械工程学院 辽宁 沈阳)摘要:由于地理位置原因 东北地区年际降水量变化大 科学有效地预报该地区中长期来水是合理进行水库调度工作的关键 文章讨论了二松流域和鸭绿江流域的气候特点 基于其间云峰、水丰、白山、丰满 座水库 年的年天然来水资料 运用成因分析法和周期叠加法实现了水库天然来水的预测 并验证了两方法预报中长期来水的可行性 研究数据表明:成因分析法模型与实际值拟合度超过 周期叠加法预报合格率超过 年来水预测成功率达到 其中周期叠加法预报效果更好 两种方法可有效地实现东北地区中长期来水的预报 为水文要素预报研究提供了实践参考关键词:中长期来水预报 水文要素 成因分析法 周期叠加法 方差分析中图分类号:文献标识码:文章编号:()收稿日期:基金项目:辽宁省重点研发计划(/)作者简介:徐 峥(年)男 工程师:引言云峰、水丰、白山、丰满 座水库为东北地区重要的控制性水利工程 兼顾发电和防洪等功能 其对东北地区农业、养殖业和旅游业等经济产业的发展及人民安全具有重大意义 与国内其它地区相比 东北地区水库存在径流量少、库容系数大、来水量年际变化大的特点 因此水库及水系的极端来水预报对东北地区流域的水库调度工作显得尤为重要 同时 水库流域来水成因复杂、影响因素众多 受天文、气象、水文、地理、地质、人类社会活动等多方面因素的共同影响为此 学者们运用各种方法及数学模型对此进行了大量研究 传统预测方法是基于概率论的数理统计方法 此类方法受限于基本假定和样本容量随着数学理论的发展 大量数值方法涌现 为水文预报工作开辟了新的思路和方法 如模糊数学方法、人工神经网络方法、灰色系统理论、小波分析等方法 各种方法的不足在于:模糊数学法因其信息模糊化的特点使得客观性降低 人工神经网络模型外延泛化能力较差 灰色系统理论适用范围 对于水文数据变化趋势 拟合灰度难以控制 因此精度难以提高 小波分解预报精度受限于其拓延方式 预报效果不稳定 所以每种方法都有其各自的适用条件 或存在有待深入研究的问题 仍需进一步的探索 成因分析法与周期叠加法将抽象的复杂问题转变为具体的数学模型 在中长期径流定量预报工作中应用广泛 是该领域的重要研究课题本文采用成因分析法与周期叠加法分别研究上述区域水文要素的特点 并对研究区域中长期来水情况进行分析和预报 以探究两种方法中长期来水预报在东北地区的准确性和可行性 研究区概况预报对象选取云峰、水丰、白山、丰满 座水库 其 集 水 面 积 分 别 为、其中前两座水库位于鸭绿江干流 后两座水库位于第二松花江干流 可汇集相应江段上游来水 是东北地区重要的控制性水利工程鸭绿江流域地处温带与亚寒带的交界处 其华夏系构造及山体走向与海岸线一致 构成了天然的气候屏障 阻挡了冬季盛行的西北气流和夏 年第 期水利规划与设计水文水资源季盛行的东南、西南气流 对该地区大气环流及天气系统具有重要影响作用 加上山体高大地形、地势复杂 形成山地气候 以上都使得该研究区内的气候具有特殊性 第二松花江流域地处温带季风气候区 气候特征显著 气候状况受季风强度影响巨大 另外 由于地理位置及地形原因 影响流域大暴雨的地面天气北上台风、蒙古气旋、华北气旋、东北气旋、冷涡和冷峰等多种天气系统都会对流域水文产生影响 造就了复杂多样的天气系统 综上 研究区流域来水状况相对复杂 因此 中长期来水预测具有较大难度鸭绿江、二松流域均位于我国东北地区 来水特性存在一定关联性 同时又具有独特性 年鸭绿江流域云峰水库、二松流域白山水库来水距平对比状况如图 所示 由图 可见 大部分年份两流域来水状况相似 但也存在差异(如 年鸭绿江丰水、二松平水 年鸭绿江偏丰、二松偏枯的情形)因此 分别对两流域进行来水预测 可证明预测方法在东北地区流域具有普适性图 年鸭绿江流域云峰水库、二松流域白山水库来水距平对比曲线 研究方法 成因分析法 方法概况普遍观点 大气环流、海洋温度及太阳活动变化规律等等与中长期水文要素状况密切相关当某种水文现象与其影响因素之间存在明确的因果关系时 可通过观测资料或实验数据分析 建立水文要素与其影响因素之间的定量关系 这种从水文现象的成因即确定性规律出发解决水文问题的方法 称为成因分析法 其在水文预报、降雨径流分析中应用较为广泛 预报因子的确定鸭绿江与第二松花江的水量主要是靠降雨补给 洪枯水量相差悬殊 地球表面约/被海洋覆盖 大气环流中水汽很大部分来自海洋 海洋温度是大气系统变化的重要因素 研究区靠近日本海面向太平洋 前一年西太平洋平均海温(日本海、黄海、东海、南海)对研究区本年降水影响很大东亚太平洋副高压直接控制和影响着太平洋上的台风活动 对中国天气系统具有很大影响 极涡是一个大尺度气旋性环流系统 它的变化对中高纬度乃至整个半球天气系统的活动和天气气候有着直接影响 太阳辐射是地球能源的主要来源 太阳活动(例如太阳黑子、耀斑等)会直接或间接的影响地球气候 太阳黑子是基本的太阳活动且易于观测其数量能够有效地反映太阳活动变化 对预测气候变化及自然灾害有着重要作用 南方涛动指数()反映了太平洋东西两侧气压增强和减弱的演变情况 变化与厄尔尼诺的发生密切相关(为负数 对应厄尔尼诺事件 反之 为正数则对应拉尼娜事件)因此 上述影响因子能够作为水文预报的依据来源 且常用于中国北方地区流域的中长期水文要素预报 故可选择上述影响因子作为此次水文预报的主因子 即:前一年 月至本年 月降水量()前一年西太平洋海表温度距平指数()前一年东亚太平洋副高压面积指数()前一年亚洲极涡强度指数()前一年年平均太阳黑子数()当年 月平均太阳黑子数()前一年年平均南方涛动值()模型建立根据研究地区历史年降水量和预报因子数据进行线性拟合 建立数学模型为:()式中 年降水量 预报因子数据 拟合系数 周期叠加法 方法概况周期叠加法是水文要素预报常用的有效方法 将随时间变化的水文要素看作有限个周期波的叠加 公式为:()()()()水文水资源水利规划与设计 年第 期式中()水文要素序列()第 个周期波序列()误差值通过周期分析方法可以将等时距的水文要素观测样本分离成若干个周期波 将分离出的周期波外延后再进行线性叠加 便可获得预报结果 周期分离过程如图 所示周期分析的方法有简单分波法、傅里叶分析法、小波分析法、方差分析法等 本文采用方差分析法识别周期设某水文要素随时间变化的等时矩样本序列为()()将其分为 组 ()的取值如下:(为偶数时)(为奇数时)()将样本序列排成表 的形式 其中 为组别 为每组的项数 为组平均值表 试验周期分组表每项组数 试验周期分组 对于不同的 可计算得到相应的方差比 为:/()()()()()()当 分别取 时 可计算得 个 值 由、和选定的信度 可以查出相应的 个 挑选最大的 值 与对应的 值比较 如果 则表明在这一信度上不存在周期 要重新选择信度 如果 则表明存在周期长度为 的周期 各组均值即为第一周期波各年的振幅 按年份排列起来就构成了第一周期波序列 然后从样本中减去第一周期波序列 形成新序列 再重复上述过程 寻找新周期 直到不能识别或不想识别为止 然后对所识别的周期波进行外延及线性叠加即可进行预测 数据处理以云峰水库为例 对 年的年天然来水量进行周期识别 各分组单因素方差分析结果见表 可知 在分组数 时 统计量取最大值 其相伴概率 满足显著性水平 的 检验 可以相信存在周期长度为 年的第一周期 各组均值即为周期波隔年的振幅 同样 根据“新序列”识别出周期长度为 的第二周期 继续分析未能识别出更多周期 故水库天然来水量存在、年共两个周期然后将各周期波进行叠加 便可得到预测结果 将 年各水库预报拟合值与实测序列进行对比 对比结果如图 所示 结果与讨论 成因分析法分别利用白山、丰满、云峰、水丰水库长序列图 周期分离流程图 年第 期水利规划与设计水文水资源表 各分组下单因素方差分析结果分组数 原序列新序列 统计量 值(相伴概率)统计量 值(相伴概率)资料 进行多元线性拟合分析 以丰满水库为例 所得模型拟合公式如下:()式中 年天然年来水量 预报因子数据预测模型拟合分析情况见表 根据表 可知 统计量对应 通过显著性指标 的 检验 可以认定该模型满足可靠性要求以丰满水库为例 其模型预测结果与实际值对比情况见表 年天然来水预测结果为偏枯实测天然来水量与成因分析模型预测量拟合对比如图 所示 丰满、白山、云峰、水丰水库年天然来水量拟合曲线相关系数分别为 、其中水丰水库拟合效果最好 由图 可见 成因分析法整体拟合效果较好 但对于大波动数据拟合结果偏差表 预测模型拟合分析表模型平方和自由度均方回归 残差 总计 周期叠加法以云峰水库为例 其水文要素拟合具体情况见表 对 年云峰水库天然来水量预测为偏丰根据水文情报预报规范 以实测入库径流量的作为允许误差 绝对误差小于允许误差则为合格预报 云峰水库预测不合格年份有、年共 次 模型合格率为 其余研究区域模型合格率均可达到 以上 实例分析 年 二松流域白山、丰满水库实际天然图 周期叠加法拟合效果图水文水资源水利规划与设计 年第 期图 成因分析法拟合效果图表 丰满水库成因分析法预测与实际来水对比分析表年份相对误差/年份相对误差/来水量较历年均值分别偏少 和 来水性质分别为平水和偏枯 鸭绿江流域云峰、水丰水库实际天然来水量较历年均值分别偏多 和 来水均偏丰 成因分析、周期叠加两种方法预测结果见表 由表 可知 绝大多数预测结果相对误差均在 以内 可以满足预报要求 仅运用成因分析法预测云峰水库来水情况偏差较大 综上可以认为上述两种方法可以有效预报东北地区中长期天然来水情况表 周期叠加法预报数据年份相对误差/年份相对误差/注:绝对误差 实际来水量 预测来水量 相对误差 绝对误差/实际来水量 年第 期水利规划与设计水文水资源表 年天然来水预测结果分析白山丰满云峰水丰实际天然来水较历年平均/年来水性质平水偏枯偏丰偏丰成因分析法预测来水性质平水偏枯偏枯偏丰与实际来水相对误差/预测结论正确正确错误正确周期叠加法预测来水性质平水偏枯偏丰偏丰与实际来水相对误差/预测结论正确正确正确正确 结论与展望本文以中长期气候预报为理论基础 基于我国东北地区二松和鸭绿江流域 年的来水情况 运用成因分析与周期叠加两种方法建立了相应预报模型 并应用于研究地区 年来水预报主要得到以下结论()成因分析法选取影响研究区域气候的关键因素为预报因子 建立数学模型在 年的预报中 座水库来水有 座定性预报正确 在一定程度上能够起到预报作用()基于长序列历史资料 采用方差分析方法实现分离周期 从历史来水情况及预报检验来看 周期叠加法预报效果良好 从 年实际预报结果来看 周期叠加法可以正确完成定性预报为了提高 座水库中长期来水预报的准确性还应该综合考虑以下方面:()本文采用的两种方法都需要大量的历史数据 且数据越多 模型越准确 故应尽可能多的选取历史数据()模型方面 进一步研究建立关于预报误差的新模型 前一年的数据资料可以及时作为下一年预报的数据依据 以提高预报准确性参考文献 徐宗学 刘琳 杨晓静.极端气候事件与旱涝灾害研究回顾与展望.中国防汛抗旱():.王彩翼.考虑中长期水文预报的北方水库综合运用经济效益分析.东北水利水电():.李敏.基于可变模糊集理论的径流预报方法研究.大连理工大学.():.():.祖佳.基于优化算法的组合模型在水库年径流预测中的应用.水利规划与设计():.庞园 曾慧 张明珠.改进的灰色关联法在流溪河流域地下水环境质量评价中的应用研究.广东水利水电():.等.:():.赵永龙 丁晶 邓育仁.相空间小波网络模型及其在水文中长期预测中的应用.水科学进展():.王子臣 冯天琼.鸭绿江流域的水文概况.水文().孙虹 李鸿雁 郭道华 等.基于物理成因识别的第二松花江汛期径流预报.水利水电技术().孟悦.松花江流域气候年代际变化分析及预测.兰州大学.()():.():.().王姣妍.气候变化影响下的新疆降水量演变规律.水利规划与设计():.喇承芳 喇承龙 刘根生.周期均值叠加法在黄河上游中长期径流预报中的应用.甘肃水利水电技术():.李俊伟.周期叠加外延最近邻抽样回归组合模型在韩江流域年径流预测中的应用.广东水利水电 ():.彭贵芬 刘盈曦.云南年降水多周期叠加与连续特大气象干旱关系及趋势预估研究.云南大学学报(自然科学版)():.周建康.实用水文统计 .北京:中国水利水电出版社.