CN109840587A - 基于深度学习的水库入库流量预测方法 - Google Patents

Abstract

本发明涉及流域水情预测技术领域，公开了一种基于深度学习的水库入库流量预测方法，包括如下步骤：取得历史数据，利用DBN模型进行学习，取得各控制站及水库的历史流量数据与水库的入库流量之间的对应关系，进而取得无雨情况下预测的水库入库流量、有雨情况下预测的水库入库流量和有雨情况下的差值delta，通过LSTM训练学习，取得有雨情况下预测的入库流量差值，进而取得最终预测的水库入库流量。本发明基于深度学习的水库入库流量预测方法，将深度置信网络与长短期记忆网络算法相融合应用于入库流量的预测，提高了对入库流量的预报精度，并且提升了模型的可靠性和可扩展性。

Description

基于深度学习的水库入库流量预测方法

技术领域

本发明涉及流域水情预测技术领域，具体涉及一种基于深度学习的水库入库流量预测方法。

背景技术

目前，对于水库入库流量预测的方法主要分为两类：基于传统的水文学研究方法以及基于人工神经网络等传统机器学习的方法。在基于传统的水文学研究方法中，一般都是通过计算土壤含水量、降雨量、上下游水库蓄水泄水情况等多种物理量值，并通过模拟产汇流机制来预报水情。传统的水库入库流量预测模型一般具有明确的物理关系，但相对来说，也存在一定问题：物理模型比较复杂，很难收集到建立模型所需的所有按时间空间分布的细节数据，例如地形数据、河道特性数据、流域土壤特性分布和降雨径流数据；此外，作为非线性过程的产汇流过于复杂，难以通过物理模型准确模拟；再者，不同流域的情况千差万别，而所用的水文模型一般又很难将各种复杂的情况都充分考虑。

在基于人工神经网络等传统机器学习的方法中，代表性的有决策树、各种聚类算法、支持向量机、神经网络等，这些方法一般从历史水文数据入手，并从数据之间的关系出发寻找内在联系的模式，通过对历史数据分析和学习从而对未来某时段的入库流量进行预测。这些方法虽然相对于传统方法在预测精度上有所提升，但是这些方法通常都需要对数据进行提前分析，再把分析结果作为模型的输入。并且大多数都是利用浅层模型在学习，容易陷入局部最优、计算时间过长等等问题，而且不能模拟大规模复杂的数学计算，所以在入库流量预测方面也难以取得理想效果。此外，这些网络模型都较为单一、结构较为简单，缺乏对具体问题的分析和设计，扩展性较差。

发明内容

本发明的目的就是针对上述技术的不足，提供一种基于深度学习的水库入库流量预测方法，将深度置信网络与长短期记忆网络算法相融合应用于入库流量的预测，提高了对入库流量的预测精度，并且提升了模型的可靠性和可扩展性。

为实现上述目的，本发明所设计基于深度学习的水库入库流量预测方法，包括如下步骤：

A)将水库周围流域范围划分为设置有若干控制站的第一监测区和设置有若干雨量站的第二监测区，按照时间点取得水库的历史流量数据、控制站的历史流量数据以及雨量站的历史雨量数据，流量数据是每个时间点测量一次的瞬时流量，雨量数据是相邻两个时间点之间的累计雨量，每个时间点上所有雨量站的雨量数据之和大于10mm为有雨情况，小于等于10mm为无雨情况；

B)利用深度置信网络(DBN)模型对无雨情况下各控制站的历史流量数据进行学习，选取各控制站以及水库在无雨情况下的历史流量数据，输入DBN模型，取得各控制站及水库的历史流量数据与水库的入库流量之间的对应关系，进而取得无雨情况下预测的水库入库流量F1；

C)选取各控制站以及水库在有雨情况下的历史流量数据，输入DBN模型，通过所述步骤B)取得的无雨情况下各控制站及水库的历史流量数据与水库的入库流量之间的对应关系，取得有雨情况下预测的水库入库流量F2，然后让有雨情况下预测的水库入库流量F2与有雨情况下的真实水库入库流量作差，取得有雨情况下的差值delta；

D)通过长短期记忆网络(LSTM)以流域范围内雨量站的历史雨量数据作为LSTM的输入，并以所述步骤C)中取得的有雨情况下的流量差值delta作为标签，经过训练学习取得有雨情况下预测的入库流量差值F3；

E)将所述步骤B)中取得的无雨情况下预测的水库入库流量F1、所述步骤C)中取得的有雨情况下预测的水库入库流量F2和所述步骤D)中取得的有雨情况下预测的入库流量差值F3融合取得最终预测的水库入库流量F＝F1+F2+F3。

优选地，所述步骤B)中，选取各控制站以及水库在无雨情况下的历史流量数据，转化为一个一维矩阵输入DBN模型，DBN模型是一种概率生成模型，其组成元件是受限玻尔兹曼机(Restricted Boltzmann Machine，RBM)，该DBN模型由1个RBM和1个人工神经网络(ANN)组成，其中RBM由一个显层和一个隐藏层构成，ANN由一个隐藏层和一个输出层构成，DBN模型中各全连接层的输出定义为如下公式：

y^l＝f(w^lx^l-1+b^l)

式中，y^l为第l层的输出，w^l为第l层的全连接的权值，x^l-1为第l-1层的输出，b^l为第l层的偏置，f为激活函数，其中RBM和ANN中隐藏层的激活函数是sigmoid激活函数，定义为如下公式：

ANN中输出层没有激活函数，为线性输出。

优选地，DBN模型学***均平方误差作为损失函数来计算地真值与估计结果的差值，并通过随机梯度下降法进行优化，假设给定N个训练样本，y′和y分别代表预测的速度向量和真值，则损失函数定义为如下公式：

通过DBN不断对各控制站以及水库历史流量数据的学习和特征提取，找到各控制站以及水库历史流量数据与入库流量之间的对应关系，最后取得各时间点无雨情况下预测的入库流量值F1，其中，模型的输出是1×1的矩阵。

优选地，所述步骤D)中，根据雨量站的位置关系，从空间位置角度对雨量站进行粗粒度化聚类，分成若干个界限明显的区域，然后在每个时间点上对每个区域雨量站的历史雨量数据进行求和，再将这若干个求和的值作为LSTM的输入，并以所述步骤C)中取得的有雨情况下的流量差值delta作为LSTM网络训练的标签，LSTM由1个cell组成，并且进行时间步展开，时序持续长度为四倍历史天数，每做一次预测，先输入四倍历史天数个数时间点的历史雨量数据，形成一维矩阵，训练LSTM时，将无雨情况时的历史雨量数据和标签都置为0，并且将值小于0的delta值也置为0，以去除不合理数据，LSTM优化的损失函数为如下公式：

式中，为N个训练样本，y′和y分别代表预测的速度向量和真值，最后经过训练学习取得有雨情况下预测的入库流量差值F3。

优选地，时间点选取一天中的2点、8点、14点、20点，即相邻时间点间隔6小时，为反映流量和降雨量的变化，历史流量数据和历史雨量数据为预测时间点前9天至前1天的流量数据和雨量数据。

本发明与现有技术相比，具有以下优点：

1、通过结合两种深度学习模型以及对数据的分类，充分有效地对相关水文数据进行了分析与挖掘，建立了最适合的水库流量预测模型；

2、在数据层面，按有雨、无雨情况对流量数据进行分类、组合，充分考虑了降雨量对入库流量的影响；在模型层面，首先，利用深度置信网络优秀的性能，克服了一般神经网络因随机初始化权值参数而容易陷入局部最优和训练时间长的缺点，使其能够在有雨、无雨两种情况下精准地挖掘流量数据中隐含的规律和特征，从而得到有雨、无雨两种情况下入库流量的预测值，其次利用长短期记忆网络在拟合时序数据上优秀的性能来处理在有雨情况下雨量与流量差delta的拟合问题，得到有雨情况下入库量量的预测差值；

3、通过有机地将上述三个部分得到的结果融合起来得到最后预测的入库流量值，最大程度的发挥各模型的优势并且充分利用了已有的数据资源，使得模型的预测的效果更加精确，也使模型具有更好的鲁棒性和适用性。

附图说明

图1为本发明基于深度学习的水库入库流量预测方法的流程图；

图2为深度置信网络模型结构图；

图3为长短期记忆网络模型结构图；

图4为实例预测值与真值的拟合曲线图。

具体实施方式

下面结合附图和具体实施例对本发明作进一步的详细说明。

一种基于深度学习的水库入库流量预测方法，其特征在于：包括如下步骤：

A)将水库周围流域范围划分为设置有三个控制站的第一监测区和设置有六百个雨量站的第二监测区，按照时间点取得水库的历史流量数据、控制站的历史流量数据以及雨量站的历史雨量数据，时间点选取一天中的2点、8点、14点、20点，即相邻时间点间隔6小时，流量数据是每隔6小时测量一次的瞬时流量，雨量数据是6小时的雨量数据，历史流量数据和历史雨量数据为预测时间点前9天至前1天的流量数据和雨量数据，定义每个时间点上所有雨量站的雨量数据之和大于10mm为有雨情况，小于等于10mm为无雨情况；

B)利用深度置信网络(DBN)模型对无雨情况下各控制站的历史流量数据进行学习，选取各控制站以及水库在无雨情况下的历史流量数据，输入DBN模型，取得各控制站及水库的历史流量数据与水库的入库流量之间的对应关系，进而取得无雨情况下预测的水库入库流量F1，

其中，如图2所示，选取各控制站以及水库在无雨情况下的历史流量数据，转化为一个一维矩阵输入DBN模型，DBN模型是一种概率生成模型，其组成元件是受限玻尔兹曼机(Restricted Boltzmann Machine，RBM)，该DBN模型由1个RBM和1个人工神经网络(ANN)组成，其中RBM由一个显层和一个隐藏层构成，ANN由一个隐藏层和一个输出层构成，DBN模型中各全连接层的输出定义为如下公式：

y^l＝f(w^lx^l-1+b^l)

式中，y^l为第l层的输出，w^l为第l层的全连接的权值，x^l-1为第l-1层的输出，b^l为第l层的偏置，f为激活函数，其中RBM和ANN中隐藏层的激活函数是sigmoid激活函数，定义为如下公式：

ANN中输出层没有激活函数，为线性输出，DBN模型学***均平方误差作为损失函数来计算地真值与估计结果的差值，并通过随机梯度下降法进行优化，假设给定N个训练样本，y′和y分别代表预测的速度向量和真值，则损失函数定义为如下公式：

通过DBN不断对各控制站以及水库历史流量数据的学习和特征提取，找到各控制站以及水库历史流量数据与入库流量之间的对应关系，最后取得各时间点无雨情况下预测的入库流量值F1，其中，模型的输出是1×1的矩阵；

C)选取各控制站以及水库在有雨情况下的历史流量数据，输入DBN模型，通过所述步骤B)取得的无雨情况下各控制站及水库的历史流量数据与水库的入库流量之间的对应关系，取得有雨情况下预测的水库入库流量F2，然后让有雨情况下预测的水库入库流量F2与有雨情况下的真实水库入库流量作差，取得有雨情况下的差值delta；

D)通过长短期记忆网络(LSTM)以流域范围内雨量站的历史雨量数据作为LSTM的输入，并以所述步骤C)中取得的有雨情况下的流量差值delta作为标签，经过训练学习取得有雨情况下预测的入库流量差值F3，其中，根据雨量站的位置关系，从空间位置角度对雨量站进行粗粒度化聚类，分成若干个界限明显的区域，然后在每个时间点上对每个区域雨量站的历史雨量数据进行求和，再将这若干个求和的值作为LSTM的输入，并以所述步骤C)中取得的有雨情况下的流量差值delta作为LSTM网络训练的标签，LSTM是一种特殊的循环神经网络，其关键组成结构是细胞cell，一个LSTM cell是由3个门限结构和1个状态向量传输线组成的，门限分别是遗忘门、传入门、输出门，其中状态向量传输线负责长程记忆，3个门限负责短期记忆的选择，因此门限设置可以对输入向量进行删除或者添加操作，从而选择性地让信息通过，因此LSTM能够学习长期依赖，故非常适合用于处理雨量这种与时间序列高度相关的数据，并且通过遗忘门和输出门忘记部分信息可以解决梯度消失的问题，如图3所示，LSTM由1个cell组成，并且进行时间步展开，时序持续长度为四倍历史天数，每做一次预测，先输入四倍历史天数个数时间点的历史雨量数据，形成一维矩阵，训练LSTM时，将无雨情况时的历史雨量数据和标签都置为0，并且将值小于0的delta值也置为0，以去除不合理数据，LSTM优化的损失函数为如下公式：

式中，为N个训练样本，y′和y分别代表预测的速度向量和真值，最后经过训练学习取得有雨情况下预测的入库流量差值F3；

E)将所述步骤B)中取得的无雨情况下预测的水库入库流量F1、所述步骤C)中取得的有雨情况下预测的水库入库流量F2和所述步骤D)中取得的有雨情况下预测的入库流量差值F3融合取得最终预测的水库入库流量F＝F1+F2+F3。

如图4是用该模型对三峡水库入库流量进行预测得到的真实流量和预测流量的拟合曲线图，其中纵轴表示入库流量，横轴表示时间序列，这里选择的是从2017年5月10日的2点到10月31日20点的数据，图中虚线表示模型的预测入库流量值，实线表示对应的真实流量值，可以看到预测入库流量值曲线和真实流量值曲线之间拟合程度是非常高的，从而也可以反映出模型预测的高准确率。

Claims (5)

1.一种基于深度学习的水库入库流量预测方法，其特征在于：包括如下步骤：

A)将水库周围流域范围划分为设置有若干控制站的第一监测区和设置有若干雨量站的第二监测区，按照时间点取得水库的历史流量数据、控制站的历史流量数据以及雨量站的历史雨量数据，流量数据是每个时间点测量一次的瞬时流量，雨量数据是相邻两个时间点之间的累计雨量，每个时间点上所有雨量站的雨量数据之和大于10mm为有雨情况，小于等于10mm为无雨情况；

B)利用深度置信网络(DBN)模型对无雨情况下各控制站的历史流量数据进行学习，选取各控制站以及水库在无雨情况下的历史流量数据，输入DBN模型，取得各控制站及水库的历史流量数据与水库的入库流量之间的对应关系，进而取得无雨情况下预测的水库入库流量F1；

C)选取各控制站以及水库在有雨情况下的历史流量数据，输入DBN模型，通过所述步骤B)取得的无雨情况下各控制站及水库的历史流量数据与水库的入库流量之间的对应关系，取得有雨情况下预测的水库入库流量F2，然后让有雨情况下预测的水库入库流量F2与有雨情况下的真实水库入库流量作差，取得有雨情况下的差值delta；

D)通过长短期记忆网络(LSTM)以流域范围内雨量站的历史雨量数据作为LSTM的输入，并以所述步骤C)中取得的有雨情况下的流量差值delta作为标签，经过训练学习取得有雨情况下预测的入库流量差值F3；

E)将所述步骤B)中取得的无雨情况下预测的水库入库流量F1、所述步骤C)中取得的有雨情况下预测的水库入库流量F2和所述步骤D)中取得的有雨情况下预测的入库流量差值F3融合取得最终预测的水库入库流量F＝F1+F2+F3。

2.根据权利要求1所述基于深度学习的水库入库流量预测方法，其特征在于：所述步骤B)中，选取各控制站以及水库在无雨情况下的历史流量数据，转化为一个一维矩阵输入DBN模型，DBN模型是一种概率生成模型，其组成元件是受限玻尔兹曼机(RestrictedBoltzmann Machine，RBM)，该DBN模型由1个RBM和1个人工神经网络(ANN)组成，其中RBM由一个显层和一个隐藏层构成，ANN由一个隐藏层和一个输出层构成，DBN模型中各全连接层的输出定义为如下公式：

y^l＝f(w^lx^l-1+b^l)

式中，y^l为第l层的输出，w^l为第l层的全连接的权值，x^l-1为第l-1层的输出，b^l为第l层的偏置，f为激活函数，其中RBM和ANN中隐藏层的激活函数是sigmoid激活函数，定义为如下公式：

ANN中输出层没有激活函数，为线性输出。

3.根据权利要求2所述基于深度学***均平方误差作为损失函数来计算地真值与估计结果的差值，并通过随机梯度下降法进行优化，假设给定N个训练样本，y′和y分别代表预测的速度向量和真值，则损失函数定义为如下公式：

通过DBN不断对各控制站以及水库历史流量数据的学习和特征提取，找到各控制站以及水库历史流量数据与入库流量之间的对应关系，最后取得各时间点无雨情况下预测的入库流量值F1，其中，模型的输出是1×1的矩阵。

4.根据权利要求1所述基于深度学习的水库入库流量预测方法，其特征在于：所述步骤D)中，根据雨量站的位置关系，从空间位置角度对雨量站进行粗粒度化聚类，分成若干个界限明显的区域，然后在每个时间点上对每个区域雨量站的历史雨量数据进行求和，再将这若干个求和的值作为LSTM的输入，并以所述步骤C)中取得的有雨情况下的流量差值delta作为LSTM网络训练的标签，LSTM由1个cell组成，并且进行时间步展开，时序持续长度为四倍历史天数，每做一次预测，先输入四倍历史天数个数时间点的历史雨量数据，形成一维矩阵，训练LSTM时，将无雨情况时的历史雨量数据和标签都置为0，并且将值小于0的delta值也置为0，以去除不合理数据，LSTM优化的损失函数为如下公式：

式中，为N个训练样本，y′和y分别代表预测的速度向量和真值，最后经过训练学习取得有雨情况下预测的入库流量差值F3。

5.根据权利要求1所述基于深度学习的水库入库流量预测方法，其特征在于：时间点选取一天中的2点、8点、14点、20点，即相邻时间点间隔6小时，历史流量数据和历史雨量数据为预测时间点前9天至前1天的流量数据和雨量数据。