基于VMD-LSTM的水库中长期径流预测模型.pdf

水 力 发 电第 卷第 期 年 月 .基于 的水库中长期径流预测模型朱非林 侯添甜 梁一帆 任瑞杰(河海大学水文水资源学院 江苏 南京)摘 要:针对径流时间序列非平稳、非线性、变异性强等复杂特点及传统的直接预测方法干扰性强等问题 采取“分解预测重构”的策略 耦合变分模态分解()方法和长短期记忆神经网络()构建基于 的中长期径流预测模型 首先采用 方法对历史中长期径流序列进行信号分解和预处理 得到一系列相对平稳的子序列 然后利用 模型对每个子序列进行建模和预测 最后加和各个子序列结果得到重构的径流预测 池潭水库的月径流预测结果表明 与单一 模型比 其均方根误差、平均绝对误差、平均绝对百分比误差分别下降、和 时间尺度的改变会影响径流预报的准确性 预测步长从 个月增至 个月 确定性系数由 降至 关键词:中长期径流预测 分解预测重构 变分模态分解 长短期记忆神经网络 池潭水库 .中图分类号:文献标识码:文章编号:()收稿日期:基金项目:国家重点研发计划项目()国家自然科学基金资助项目()作者简介:朱非林()男 江苏南通人 教授 博士主要从事流域水库群防洪预报调度、水风光储清洁能源系统互补调度方面的研究.0引言中长期径流预报是流域水资源配置、水旱灾害防御、水利工程运行调度等诸多工作的重要决策依据 影响中长期径流预测的多方面因子与预测要素之间具有十分繁杂的联系 中长期径流预报一直是水文水资源和大气科学领域内的一个难点方向近十几年来 中长期径流预报得到了国内外学者的普遍关注 相关预报理论方法得到了快速发展与短期径流预报相比 因为没有可靠的气象预报且预见期较长、大气环流等影响因素异常复杂 中长期径流预报的方法原理普遍缺乏机理 现有的中长期径流预测方法可分为时间序列统计模型和机器学习模型 时间序列预测方法是指通过类比推导或第 卷第 期朱非林等:基于 VMDLSTM 的水库中长期径流预测模型 .时间序列线性延伸所反映出的发展方向 以预测未来情况 该方法一般适用于线性或平稳的数据序列 难以刻画序列中的非平稳特征 机器学习通过研究计算机如何模拟或实现人的学习行为来获得提高自身性能的新知识 它具有从大数据集中自动总结归纳信息的能力 可以捕获径流序列中的非平稳和非线性特征 例如 卢敏等将支持向量机()应用于径流预测中 可以较好地处理高度非线性问题 赵铜铁钢为提高径流预报精度 将随机森林模型应用于长江中上游枯水期径流预报中 等利用人工神经网络和相空间重构方法对河流动态进行了预测 尽管机器学习方法已在水文预报领域取得大量应用 仍难以满足复杂情形下的应用需求 近年来 计算机技术和人工智能已经发展至新阶段 以深度学习为代表的机器学习方法受到广泛关注 作为更有效的深度学习方法 长短期记忆神经网络()可以从原始数据中辨析出更深层次的特征数据 解决现有径流预测方法存在的长期记忆能力欠缺和梯度消失、梯度爆炸等问题在进一步提升中长期径流预测的精度和效率方面具有明显优势受气候、下垫面及人类活动等综合作用 由多种频率组成的径流序列具有随机性强、平稳性低且线性特征不明显的特点 采用单一预报模型直接预测径流的精度受到一定限制 信号分解技术可将径流序列分解为若干相对稳定的分量 有效去除序列中的干扰信息 提高信噪比 将其与机器学习模型结合 可提高预测精度 近年来 信号分解技术在诸多领域取得了成功应用 等将经验模态分解与支持向量机相结合 实现渭河流域的月径流预测 周婷等构建了基于小波分解的 模型 并对响洪甸水库径流过程进行了预测 李继清等采用极点对称模态分解()方法处理原始径流序列 提出了 耦合模型 现有研究表明 组合预测模型是进一步提升原始模型预测效果的有效手段 但传统分解方法仍面临最佳筛选次数难以确定、抗噪能力弱、分量波形混叠、趋势项粗略等问题 变分模态分解()可以有效克服现有分解方法存在的端点发散效应和分量波形混叠的缺点 具有良好的抗噪性和更坚实的数学理论基础鉴于此 本文采用“分解预测重构”的思路结合变分模态分解方法与长短期记忆神经网络 构建中长期 组合预测模型 先利用 方法将径流分解为若干相对平稳的子序列 再构建 模型对子序列进行单独预测和二次重构 从而探讨不同预测步长对径流预测精度的影响 为水库中长期径流预测提供借鉴和参考1研究方法 变分模态分解变分模式分解()具有自动调整的能力 是一种完全非递归的处理信号和分解模态的方法 该方法利用非递归方式 可以降低非平稳性和非线性强的时间序列复杂程度 分解得到多个不同频率的分量构建、分析和求解变分问题是 的核心 分解原始序列为 个子序列 其具体步骤为:为获取各子分量的单向频谱 利用希尔伯特变换对其解析信号进行逐个计算 针对各子分量 根据其中心频率 把频谱调至相应基带 依据高斯平滑度对信号进行解调 估算分解模态的频率范围 构建以模态估计带宽之和最小为目标 以所有模态之和与原始信号相等为约束条件的变分问题 数学表达式为()()()式中 和 分别为第 个模态分量和其中心频率()为单位脉冲函数利用 算法进行径流序列分解的计算流程见图 图 算法用于径流序列分解的计算流程 长短期记忆神经网络长短期记忆神经网络()是一种特殊的循环神经网络()对于有长程依赖关系的时间序列预水 力 发 电 年 月 .测问题 神经网络可以解决传统神经网络存在的弊端 如梯度消失和梯度爆炸等 在非线性径流序列预测方面具有更明显的优势如图 所示 单元由遗忘门、输入门、细胞状态 和输出门 构成:图 径流预测原理结构()遗忘门 决定应遗忘的信息 即()()()输入门 选择记录到细胞状态的信息 即()()()()()更新细胞状态 即()()输出门 选择携带到下一个神经元的信息 即()()()()图 基于 的中长期径流组合预测模型计算流程式中、均为 权 值 矩 阵、为偏置向量、分别为前一模块和当前模块的输入和输出、分别为前一模块和当前模块的细胞状态 为输入门的细胞状态()为 函数 径流预测模型本文采用“分解预测重构”的思路 将 与 两种方法相耦合 发挥两种方式各自的长处构建了基于 的中长期径流组合预测模型模型计算过程如图 所示 主要有以下几个步骤:()调试 参数 分解径流序列 得到一系列从高频到低频的平稳模态分量()()对于 分解得到的每个分量 单独构建相应的 模型 将一系列平稳分量作归一化处理 设置预测步长为、个月 通过已分解的实际径流序列训练和预测 模型 得到各成分的预测值()加和重构各分量预测值 得到组合模型预测结果()选取一些误差指标 如均方根误差等 评价模型性能 误差评价指标为了评估径流预测效果 本文以均方根误差()、平均绝对误差()和平均绝对百分比误差()作为评价指标 计算公式如下 ()()第 卷第 期朱非林等:基于 VMDLSTM 的水库中长期径流预测模型 .()()式中 为测试集的样本数 为实际径流值为模型预测值2研究区及其数据概况金溪作为福建省闽江支流富屯溪的最大一级支流 其流域属于亚热带湿润季风型山地气候 降雨具有量大但年内时程分配不均的特点 池潭水库作为金溪干流水电梯级的第一级龙头水库 设计以发电为主 兼顾防洪等综合利用 其地理位置如图 所示 研究所采取的数据为池潭水库 年年共 逐月径流资料 其中前 年数据用于模型训练 后 年的数据用于模型验证图 池潭水库的地理位置示意表 不同分解层数下各分量的中心频率分解层数 3研究结果与分析 月径流序列 分解结果 分解效果主要受分解层数 影响 若 取值过大 相邻模态分量的中心频率则会过于相近引起分量波形混叠问题 若 取值偏小 部分原始信号中的信息容易被忽略 影响后续预测精度 不同的模态区别主要在于中心频率的差异 当出现相似频率 即选取此模态数 为确定合适的模态数值 本文观察不同 值下中心频率的分布 表 为不同分解层数下各分量的中心频率 为 与 为 时中心频率趋于稳定 因此选取 作为最终的分解层数将实际径流序列分解为 个分量 分解结果见图 与原始单一径流序列相比 分解后的各分量均具有较为明显的变化趋势 其中 分量 的波形起伏最为平缓 振荡及频率大幅降低 后续各分量的变化幅度逐渐减小 但振荡频率则逐渐增大例如、波动十分剧烈 但由于其变幅相对较小 故 在径流重构中 仍占主导地位 对 预测结果的影响为探究 方法对中长期径流预测模型效果的影响 分别采用 单一模型与 组合模型对径流进行预测 预测结果见图 由图 可知 在相同预测步长下 两组模型预测的径流变化与实际径流过程较为吻合 表明 方法的非线性拟合能力较强 与单一模型相比 基于分解重构策略的 组合预测模型的预测值与实测值的接近程度明显较高 尤其在径流极值处 组合模型预测的效果更优 精度更高 结合误差分析 组合模型的预报误差与单一模型的直接预报误差相比显著降低 表明 分解可以有效提升 单一模型的径流预测精度 预测步长对 组合模型预测结果的影响 种不同预测步长下 组合模型对池潭水库径流序列预测结果对比分析如图 所示 当预测步长为 时 点线吻合程度最高 随着预测步长增加 吻合程度逐渐降低 为得到更直观的结果将不同预测步长下实测径流与组合模型预测结果做线性拟合分析 得到右列散点图 其中 为确定性水 力 发 电 年 月 .系数 表示因变量 的变异中可由自变量 解释的部分所占的百分比 即拟合程度 计算公式为 ()()()图 月径流序列的 分解结果图 组合模型与 单一模型径流预测结果对比式中 为径流系列平均值可以看出 随着预测步长增加 逐渐减小拟合效果降低 同时 拟合线斜率分别为 、图像偏离 直线的程度逐渐加大 预测精度降低 表明预测步长的增加会导致模型预测性能降低 模型效果评价为了进一步对 组合模型的有效性进行定量评估 本文选取、三种误差指标对经过变分模态分解和未经分解的模型预测结果进行了对比(见表)由表 和图 可知 当预测步长一定时 与 单一模型相比 分解能显著降低 模型的预测误差 当预测步长为第 卷第 期朱非林等:基于 VMDLSTM 的水库中长期径流预测模型 .图 不同预测步长下的 组合模型预测结果对比表 单一 模型和 组合模型评价指标计算结果预测步长模型/时 组合模型的、指标值降幅分别 为 、确 定 性 系 数 为 当预测步长为 时 组合模型的指标值分别降低 了 、确 定 性 系 数 为 当预测步长为 时 组合模型的指标值分别降低 了 、确 定 性 系 数 为 由此可见 组合模型在预测精度上明显提高图、分别为 单一模型、模型多步预测的径流误差 由表 以及图、可知不管是 单一模型还是 模型 预报误差均随预报步长增加而增大 随着预测步长的增加 单一预测模型的误差指标增幅为 左右 组水 力 发 电 年 月 .图 单一模型多步预测的径流误差图 组合模型多步预测的径流误差合预测模型的指标值增幅为 左右综上所述 本文构建的 中长期径流组合预测模型可以充分发挥 与 方法各自的优势 表现出明显优于 单一模型的预测精度和稳定性4结论中长期径流预测是水资源开发与调度、水旱灾害防治、水库运行与管理的重要支撑 按照“分解预测重构”的模型架构 本文结合信号分解与神经网络 构建了基于 的中长期径流组合预测模型 以金溪流域池潭水库的月径流预测为实例对模型预测性能进行了分析论证 主要结论如下:()基于“分解预测重构”策略的 组合预测模型结合了 方法平稳化处理技术和 神经网络模型深度学习的优势 可有效降低径流序列非平稳性对预测精度的影响 具有预测准确率高和稳定性好等优点()耦合 分解方法可以显著提高 单一模型的预测精度()时间尺度的选择影响径流预报的准确性对于 单一模型和 组合模型 预测步长增加 预测误差随之增大()模型径流预测效果良好 可以为水库中长期调度计划编制与水资源规划管理提供决策依据 也提供了相关时间序列预测的研究思路参考文献:张美英 何杰.时间序列预测模型研究综述.数学的实践与认识 ():.廖明 詹总谦 呙维 等.动态数据驱动模式下的湖泊流域降雨径流模拟.遥感学报 ():.().贺圣平 王会军 李华 等.机器学习的原理及其在气候预测中的潜在应用.大气科学学报 ():.卢敏 张展羽 冯宝平.支持向量机在径流预报中的应用探讨.人民长江 ():.赵铜铁钢 杨大文 蔡喜明 等.基于随机森林模型的长江上游枯水期径流预报研究.水力发电学报 ():.():.常新雨 周建中 方威 等.黄龙滩水库中长期径流预报方法研究.水力发电 ():.顾冬 郁怀光 吴宇 等.基于相似测点残差信息挖掘的大坝变形异常值识别方法研究/.水力发电:.:./.:.周婷 金菊良 李荣波 等.基于小波支持向量机的径流预测性能优化分析.水力发电学报 ():.李继清 王爽 段志鹏 等.基于 神经网络组合模型的中长期径流预报.应用基础与工程科学学报 ():.闫国辉 乔长录 陈伏龙.基于 模型的玛纳斯河年径流预测.中国农村水利水电():.():.殷兆凯 廖卫红 王若佳 等.基于长短时记忆神经网络()的降雨径流模拟及预报.南水北调与水利科技():.陶凯 吴定会.基于 的短期风电功率预测.控制工程 ():.(责任编辑 陈 萍)