基于经验小波消噪和神经网络融合的中长期径流预报方法技术

本发明专利技术涉及水文预报领域，公开了一种基于经验小波消噪和神经网络融合的中长期径流预报方法，针对现有技术的不足，本发明专利技术首先采用经验小波变换(EWT)对中长期径流时间序列进行分解和重构，通过消除分解模态中的高频分量对中长期径流时间序列的冗余噪声进行处理，在此基础上，以RBF神经网络模型、ELM模型和Elman神经网络模型作为基预报模型对中长期径流时间序列进行预报，并将三种基模型的预报结果作为GRNN神经网络的预报因子，对预报径流做进一步的模拟，从而对每个基模型的预报信息进行融合，提升中长期径流时间序列的预报精度。

【技术实现步骤摘要】
基于经验小波消噪和神经网络融合的中长期径流预报方法
本专利技术涉及水文学中的水文预报领域，特别涉及一种基于经验小波消噪和神经网络融合的中长期径流预报方法。
技术介绍
中长期径流预报(遇见期3天～1年)是近几十年来水文研究的热点问题之一，提高中长期径流预报精度对水资源调配、防洪、减灾等方面具有重要意义。现有中长期径流预报方法可以分为两类：物理模型和经验模型。物理模型一般通过模拟流域水循环的物理过程对未来径流进行预报。物理模型的构建需要获取较长时间段(遇见期)内的降水、蒸发和气温等气象数据，增加了采用物理模型进行中长期预报的难度。中长期水文预报中涉及的各种经验模型主要包括时间序列模型和机器学习模型。然而，以往的研究表明，时间序列模型对非线性和非平稳性较强的径流时间序列预报能力有限。在各种机器学习方法中，人工神经网络(ANN)包括反向传播神经网络(BPNN)、径向基函数(RBF)神经网络、广义回归神经网络(GRNN)、Elman神经网络和极限学习机(ELM)等在中长期径流预报领域取得了丰硕的成果，但是单一ANN模型的预报性能仍存在很大的改进空间。其中利用数据预处理技术对径流时间序列进行预处理是一种能有效提升中长期径流预报精度的方法。此外，另一种能够提升单一ANN模性预报性能的方法是将多个ANN模型的预报信息进行融合，综合考虑每个模型的优势，进而提高中长期径流预报精度。
技术实现思路
专利技术目的：针对现有技术中存在的问题，本专利技术提供一种基于经验小波消噪和神经网络融合的中长期径流预报方法，能够有效提升中长期径流时间序列的预报精度。技术方案：本专利技术提供了一种基于经验小波消噪和神经网络融合的中长期径流预报方法，其特征在于，包括以下步骤：步骤一：获取水文站点的历史实测径流数据，建立中长期径流时间序列，并将样本数据分为训练期和检验期；步骤二：采用EWT将径流时间序列分解为几个独立的经验模态，去除分解经验模态中频率最高的模态，并对剩余的经验模态进行线性求和，得到重构径流时间序列以消除原始径流时间序列的冗余噪声；步骤三：对经EWT处理后的径流时间序列进行相空间重构，构造相空间矩阵作为基预报模型的输入因子，建立ANN基预报模型的输入、输出矩阵；步骤四：采用步骤三确定的训练期输入、输出矩阵分别对RBF、ELM和Elman三种神经网络模型进行训练；步骤五：将三种基预报模型的输出作为GRNN模型的预报因子，再次建立输入、输出矩阵，对GRNN模型进行训练，并将检验样本代入训练好的GRNN模型，得到检验期的预测值；步骤六：采用四种评价指标包括均方根误差(RMSE)、平均绝对误差(MAE)、平均绝对百分误差(MAPE)和相关系数(R)评价对预报结果进行评价。进一步地，在所述步骤二中，对于原始径流时间序列x(t)，EWT的分解过程如下：(1)根据快速傅立叶变换算法(FFT)计算原始径流时间序列x(t)的傅立叶频谱F(ω)；(2)将傅立叶频谱F(ω)的频域[0,π]自适应的分割为K个带宽不等的频带，[0,ω1],[ω1,ω2],...,[ωK-1,π]，其中，小波滤波器的边界ωn,n＝1,2,...,N-1选取为两个连续局部极大值之间的中间频率；(3)根据Meyer小波构造经验小波，确定经验小波函数ψk(ω)和经验尺度函数式中，β(x)＝x4(35-84x+70x2-20x3)；(4)重构原始径流时间序列得到不同的经验模态。优选地，在所述(4)中，原始径流时间序列通过下式进行重构：式中，*为卷积运算，Wx(0,t)为近似系数，Wx(k,t)为细节系数，近似系数和细节系数分别通过下式给出：式中，ψk(t)和ψk(t)分别表示经验小波函数和经验尺度函数，和分别表示ψk(τ-t)和的复共轭，F-1(·)表示傅里叶逆变换，和分别表示ψ(ω)和的傅里叶变换；则经验模态uk(t)的定义如下：优选地，RBF神经网络的隐含层使用高斯传递函数作为激活函数，高斯激活函数定义为：式中，x＝[x1,x2,...,xn]T为n维的输入向量，ci＝[ci1,ci2,...,cin]T为第i个隐含层神经元的中心，qi表示高斯函数的宽度，·为欧氏范数；则输出层第j个结点的响应为：式中，wij为第i个隐含层结点与第j个输出层结点之间的连接权重，h和m分别表示隐含层和输出层结点个数；极限学习机是一种单隐层前馈神经网络，给定K个训练样本{xk,yk},k＝1,2,...,K，其中X＝[x1,x2,...,xK]T∈RM为输入向量，Y＝[y1,y2,...,yK]T∈RN为输出向量；则第k个训练样本的第j个隐含层神经元的输出表示为其中，J表示隐含神经元数目，g(·)表示激活函数，βj和aj分别表示输入层和隐含层之间想连接权重和偏置；ELM模型第k个训练样本的输出可以表示为：式中，wj表示隐含层和输出层之间的权重；隐含层核映射矩阵表示如下：ELM网络的输出表示为：(11)；式中，w＝[w1,w2,...,wJ]T；ELM模型的目标是找到最合适的w，使得网络输出值和实测值Y之间的误差最小，通过求解以下优化问题得到隐含层与输出层之间的关联权重：方程(12)的解根据MoorePenrose(MP)广义逆求得：式中，是矩阵D的广义逆矩阵；则ELM模型第k个训练样本的输出表示为：Elman神经网络是一种动态递归神经网络。除了输入层、隐含层和输出层外，Elman神经网络还包括一个特殊的递归层，用于记忆隐含层神经元以前时刻的输出信息，然后将该信息用作隐藏层的输入。进一步地，在所述步骤五中，假设已知随机向量x，y和联合概率密度函数p(x,y)，若x的观测值为X，则y关于x的回归表示为：式中，E[y|X]为给定输入X下输出Y的预测值，密度函数p(X,y)未知，由Parzen非参数估计进行估算：式中，Xj和Yj分别表示第j个样本的实测输入值和实测输出值，N是样本集的个数，d是x的维数，σ是平滑参数；根据式(15)和(16)，GRNN模型的输出可以表示如下：优选地，在所述步骤六中，RMSE，MAE，MAPE及R的计算公式如下：式中，qp(i)为预测值；qo(i)为实测值；和分别表示为实测值和预测值的均值；N是样本集的个数。有益效果：与现有技术相比，通过本专利技术建立的中长期径流预报方法，能够取得下列有益效果：1)本专利技术采用经验小波变换对径流时间序列进行分解与重构，能够去除天然径流时间序列中的冗余噪声，降低径流时间序列的非线性波动对径流预报的干扰，进而提高预报模型的性能；2)针对单一ANN模型的预报精度难以达到预期水平，且每种模型各具优缺点，以三种ANN预报模型的预报结果作为输入因子，采用GRNN神经网络模型对月径流时间序列进行进一步预报，以综合不同ANN预报模型的优势，实现各基模型预报信息的融合，从而获得精度更高的径流预报信息。附图说明图1为基于经验小波消噪和神经网络融合的中长期径流预报模型流程图；图2为乌东德站径流时间序列的EWT分解以及原始与消噪径流时间序列的对比图；图3为乌东德站检验期GNE和EWT-GNE模型的月径流预测结果对比图。具体实施方式下面结合附图对本专利技术进行详细的介绍。实施方式1：本专利技术以长江上游流域乌东德水文站月径流时间序列为实施例，进行实例仿真，以验证本专利技术的效果。图1为本专利技术提供本文档来自技高网...

【技术保护点】
1.一种基于经验小波消噪和神经网络融合的中长期径流预报方法，其特征在于，包括以下步骤：步骤一：获取水文站点的历史实测径流数据，建立中长期径流时间序列，并将样本数据分为训练期和检验期；步骤二：采用EWT将径流时间序列分解为几个独立的经验模态，去除分解经验模态中频率最高的模态，并对剩余的经验模态进行线性求和，得到重构径流时间序列以消除原始径流时间序列的冗余噪声；步骤三：对经EWT处理后的径流时间序列进行相空间重构，构造相空间矩阵作为基预报模型的输入因子，建立ANN基预报模型的输入、输出矩阵；步骤四：采用步骤三确定的训练期输入、输出矩阵分别对RBF、ELM和Elman三种神经网络模型进行训练；步骤五：将三种基预报模型的输出作为GRNN模型的预报因子，再次建立输入、输出矩阵，对GRNN模型进行训练，并将检验样本代入训练好的GRNN模型，得到检验期的预测值；步骤六：采用四种评价指标包括均方根误差(RMSE)、平均绝对误差(MAE)、平均绝对百分误差(MAPE)和相关系数(R)评价对预报结果进行评价。

【技术特征摘要】
1.一种基于经验小波消噪和神经网络融合的中长期径流预报方法，其特征在于，包括以下步骤：步骤一：获取水文站点的历史实测径流数据，建立中长期径流时间序列，并将样本数据分为训练期和检验期；步骤二：采用EWT将径流时间序列分解为几个独立的经验模态，去除分解经验模态中频率最高的模态，并对剩余的经验模态进行线性求和，得到重构径流时间序列以消除原始径流时间序列的冗余噪声；步骤三：对经EWT处理后的径流时间序列进行相空间重构，构造相空间矩阵作为基预报模型的输入因子，建立ANN基预报模型的输入、输出矩阵；步骤四：采用步骤三确定的训练期输入、输出矩阵分别对RBF、ELM和Elman三种神经网络模型进行训练；步骤五：将三种基预报模型的输出作为GRNN模型的预报因子，再次建立输入、输出矩阵，对GRNN模型进行训练，并将检验样本代入训练好的GRNN模型，得到检验期的预测值；步骤六：采用四种评价指标包括均方根误差(RMSE)、平均绝对误差(MAE)、平均绝对百分误差(MAPE)和相关系数(R)评价对预报结果进行评价。2.根据权利要求1所述的基于经验小波消噪和神经网络融合的中长期径流预报方法，其特征在于，在所述步骤二中，对于原始径流时间序列x(t)，EWT的分解过程如下：(1)根据快速傅立叶变换算法(FFT)计算原始径流时间序列x(t)的傅立叶频谱F(ω)；(2)将傅立叶频谱F(ω)的频域[0,π]自适应的分割为K个带宽不等的频带，[0,ω1],[ω1,ω2],...,[ωK-1,π]，其中，小波滤波器的边界ωn,n＝1,2,...,N-1选取为两个连续局部极大值之间的中间频率；(3)根据Meyer小波构造经验小波，确定经验小波函数ψk(ω)和经验尺度函数式中，(4)重构原始径流时间序列得到不同的经验模态。3.根据权利要求2所述的基于经验小波消噪和神经网络融合的中长期径流预报方法，其特征在于，在所述(4)中，原始径流时间序列通过下式进行重构：式中，*为卷积运算，Wx(0,t)为近似系数，Wx(k,t)为细节系数，近似系数和细节系数分别通过下式给出：式中，ψk(t)和ψk(t)分别表示经验小波函数和经验尺度函数，和分别表示ψk(τ-t)和的复共轭，F-1(·)表示傅里叶逆变换，和分别表示ψ(ω)和的傅里叶变换；则经验模态uk(t)的定义如下：4.根据权利要求1至3中任一项所述的基于经验小波消噪和神经网络融合的中长期径流预报方法，其特征在于，在所述步骤四中，RBF神经...

【专利技术属性】
技术研发人员：彭甜，倪伟，张楚，夏鑫，纪捷，薛小明，
申请(专利权)人：淮阴工学院，
类型：发明
国别省市：江苏,32