CN109272146A - 一种基于深度学习模型和bp神经网络校正的洪水预测方法 - Google Patents

一种基于深度学习模型和bp神经网络校正的洪水预测方法 Download PDF

Info

Publication number
CN109272146A
CN109272146A CN201810965624.5A CN201810965624A CN109272146A CN 109272146 A CN109272146 A CN 109272146A CN 201810965624 A CN201810965624 A CN 201810965624A CN 109272146 A CN109272146 A CN 109272146A
Authority
CN
China
Prior art keywords
model
flood
data
basin
prediction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810965624.5A
Other languages
English (en)
Other versions
CN109272146B (zh
Inventor
王继民
朱跃龙
张�成
张鹏程
朱晓晓
张玲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201810965624.5A priority Critical patent/CN109272146B/zh
Publication of CN109272146A publication Critical patent/CN109272146A/zh
Application granted granted Critical
Publication of CN109272146B publication Critical patent/CN109272146B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services
    • G06Q50/265Personal security, identity or safety
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • General Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Linguistics (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Security & Cryptography (AREA)
  • Computing Systems (AREA)
  • Molecular Biology (AREA)
  • Educational Administration (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Physics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Primary Health Care (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于深度学习模型和BP神经网络校正的洪水预测方法,步骤为:1、对历史洪水过程数据进行归一化处理;2、对归一化后的历史洪水过程数据分析得到各监测点对出口流量的影响时间;3、利用滑动窗口从历史场次洪水数据中建立预测模型的输入和输出值,建立模型训练数据集TRSet1;4、利用TRSet1训练并建立深度学习的洪水预测模型CNNFM;5、建立实时误差校正模型训练数据集TRSet2;6、利用TRSet2训练并建立基于BP神经网络的误差校正模型BPCM;7、对实时监测的流量和雨量,利用模型CNNFM对实时数据进行预测,利用实时误差校正模型BPCM进行校正,获得最终的预测值。本发明通过深度学习自动提取数据特征的特点对训练数据进行建模,并通过BP神经网络进行实时校正,以提高模型预测的准确率。

Description

一种基于深度学习模型和BP神经网络校正的洪水预测方法
技术领域
本发明属于信息处理技术领域,特别涉及一种基于深度学习模型和BP神经网络校正的洪水预测方法。
背景技术
我国国土辽阔,河流水资源众多,洪涝灾害频发,对我国经济发展、社会的进步带来了阻碍,因此开展水文预报研究十分重要。传统上大多使用基于物理过程的概念性水文模型描述水文过程,这类方法较为成熟,能够达到较好的预测效果。但是这些模型复杂,针对不同地区模型适应性较差,模型参数率定难度较大。因此基于数据驱动的水文过程预测方法日益得到发展。近年来,我国已经建立了较为完善的水文信息监测网络,采集了大量的水文数据,这些数据中蕴含了水文过程的内在规律,如何利用数据挖掘技术建立的基于数据驱动的预测模型,提高水文过程预测的准确率,是一个重要的研究方向。
发明内容
发明目的:为了克服现有技术中存在的不足,本发明提供一种通过深度学习自动提取数据特征的特点对训练数据进行建模,并通过BP神经网络进行实时校正,以提高模型预测的准确率的基于深度学习模型和BP神经网络校正的洪水预测方法。
技术方案:为实现上述目的,本发明提供一种基于深度学习模型和BP神经网络校正的洪水预测方法,包括如下步骤:
(1)对历史洪水过程数据进行归一化处理;
(2)对归一化后的历史洪水过程数据序列进行分析,分析降雨量和蒸发量影响因素对流域出口流量的影响时间范围;
(3)利用滑动窗口从历史场次洪水数据中建立预测模型的输入和输出值,建立模型训练数据集TRSet1;
(4)建立基于深度学习的洪水预测模型CNNFM,利用训练数据TRSet1训练预测模型得到模型参数,建立流域的洪水预测模型;
(5)流域的洪水预测模型是指适应某具体流域的洪水预测模型,通过模型训练过程中,从输入参数中立定的模型参数体现了当前流域的特征。
(6)建立实时误差校正模型训练数据集TRSet2,输入与TRSet1相同,流域洪水预测模型的预测误差作为输出;
(7)建立基于BP神经网络的误差校正模型BPCM,利用TRSet2训练校正模型,获得模型参数,建立流域实时误差校正模型;
(8)流域实时误差校正模型是指适应具体流域的BP校正模型。该模型中的参数体现了当前流域的特征。
(9)对实时监测的降雨量和流量,利用流域洪水预测模型CNNFM对实时数据进行预测,利用流域实时误差校正模型BPCM进行校正,获得最终的预测值。
进一步的,所述步骤(1)中对历史洪水过程数据进行归一化处理的具体步骤如下:
采用离差标准化即min-max标准化将洪水过程预测涉及到的降雨量、蒸发量、流量和水位数据进行线性变化,通过转化公式将原始值映射到[0,1]之间;转化公式如下:
其中x*为转换后的值,x为原始值,xmin=min(x),xmax=max(x);经过min-max标准化后,原始值被规约到[0,1]之间;其中,max(x)取x对应指标的历史最大值,min(x)取x对应指标的历史最小值。
进一步的,所述步骤(2)中分析降雨量和蒸发量影响因素对流域出口流量的影响时间范围;的具体步骤如下:
对归一化后的历史洪水过程数据序列采用皮尔逊相关系数分析输出量即流域出口流量与不同时间范围的各输入量即流域内各雨量站的降雨量和蒸发站的蒸发量的相关性,确定输入量对输出量影响的时间范围,各输入量的影响范围的最大值作为整体输入量对输出量影响的时间范围。
进一步的,所述步骤(3)中建立模型训练数据集TRSet1的具体步骤如下:
利用步骤(2)中得到的时间范围作为滑动窗口的宽度,按照步长为1提取历史洪水过程数据作为输入,对应输出为预测预见期为宽度的输出序列;具体的输入量对输出量影响的最大时间范围为w,即表示某个输入量的t0时刻监测值,从t1到tw时间内,逐渐影响输出量,并最终影响消失;
上式表示预测模型的输入为I的矩阵,有m+1个输入量,包括预测量自身;输出为O的向量,预见期为n。
进一步的,所述步骤(4)中建立基于深度学习的洪水预测模型CNNFM,采用卷积神经网络构建,其中预测模型CNNFM包括卷积层、池化层和激活函数。
进一步的,根据权利要求1所述的一种基于深度学习模型和KNN实时校正的洪水预测方法,其特征在于,所述步骤(5)中建立误差校正模型训练数据集TRSet2的具体步骤如下:
输入与步骤(3)中的预测模型数据集TRSet1输入相同,输出为预测模型CNNFM的预测误差,其中对输入为模型预测输出为 误差为校正模型的输入输出数据集中每个输入输出对为
其中,
其中,I表示输入,O表示输出。
进一步的,根据权利要求1所述的一种基于深度学习模型和BP实时校正的洪水预测方法,其特征在于,所述步骤(7)中训练校正模型,具体步骤如下:
建立BP神经网络校正模型BPCM,利用TRSet2训练校正模型,通过BP神经网络建立输入数据和输出误差之间的非线性关系,获得模型参数。
有益效果:本发明与现有技术相比具有以下优点:
本发明提出的基于卷积神经网络和实时校正的洪水过程预测方法,从监测的降雨数据、流量数据等中利用数据挖掘发现数据之间的关系,相对传统的基于物理过程的预测过程来说,参数较少,同时更加容易设置相关参数。同时,在预测模型的基础上,通过结合组合校正模型,对预测值进行实时调整,提高了预测的准确率,提高了峰现时间和峰值预报准确率。
附图说明
图1是本发明的流程图;
图2是实施例中不同的时间差中di与q之间的关系图;
图3为具体实施例中屯溪流域报汛站网分布图;
图4为具体实施例中1997060608号洪水不同预见期预报情况图;
图5为具体实施例中2001050108号洪水1h预见期预报校正对比情况图;
图6为具体实施例中2001050108号洪水2h预见期预报校正对比情况图;
图7为具体实施例中2001050108号洪水4h预见期预报校正对比情况图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图1所示,本发明包括如下步骤:
(1)对历史洪水过程数据进行规格化。历史洪水过程数据中可能会存在缺失值,缺失的数据可能会存在偏差估计,导致样本数据不能很好的代表总体,因此需要对缺失值进行处理,针对洪水过程中缺失值的处理参考《水文资料整编规范》(SL247-2012)。实际的洪水过程预测往往涉及到降雨量、蒸发量、流量、水位等多种数据,,为了消除这些不同的属性和指标间存在的数量级和量纲的差异,采用min-max标准化,也称为离差标准化,它将原始数据进行线性变化,通过转化公式将原始值映射到[0,1]之间。转化公式如(1):
其中x*为转换后的值,x为原始值,xmin=min(x),xmax=max(x)。经过min-max标准化后,原始值被规约到[0,1]之间。其中,max(x)取x对应指标的历史最大值,min(x)取x对应指标的历史最小值。
(2)分析流量信息与降雨量、蒸发量之间的相关性,获取预测输入数据的时间窗口宽度。对各场历史洪水过程数据序列采用皮尔逊相关系数分析因变量(流域出口流量)与各自变量(流域内各雨量站的降雨量、蒸发站的蒸发量等)的相关性,确定各自变量对因变量影响的时间范围,各自变量的影响范围的并集作为整体输入量对输出量影响的时间范围,获取预测输入数据的长度。
若某场洪水过程有流域出口流量q序列qt0,qt1,qt2,…,qtn,某雨量站/蒸发站监测的雨量/蒸发量di序列为dit0,dit1,dit2,…,ditn,分析di对q的影响时间范围,如图2所示,逐步将di前移1个时间单位,并分别分析与q之间的相关性,相关系数一般逐渐增加,到最高值,然后逐渐减小。
若当时间差=w时,相关系数任然大于等于50%,当时间差=w+1时,相关系数小于50%,则该di对q影响的时间宽度为w,即认为dt(n-k)至dtn的降雨数据都对qtn有影响。分析所有自变量对因变量的影响时间宽度,取最大值作为步骤(3)中的建立训练数据的滑动窗口宽度。
(3)建立预测模型训练数据集。利用步骤(1)中分析得到的时间范围w作为滑动窗口的宽度,按照步长为1提取历史洪水过程数据以及同步的降雨、蒸发等作为输入,对应输出为指定预测预见期的流量。如下为一个输入输出对,所有输入量对输出量影响的最大时间范围为w。
表示预测模型的输入为I的矩阵,输出为O,有m+1个输入量,其中,di表示某个测站监测的降雨量或蒸发量序列等,q表示流量序列,预见期为n。
(4)建立基于卷积神经网络的预测模型CNNFM,并利用训练数据进行训练。模型主要包括卷积层、池化层、激活函数。卷积层采用Google深度学习框架TensorFlow来实现,使用以下程序来进行卷积神经网络前向计算过程中卷积层的模块设计。
filter_weight=tf.get_varialble('weight',[patch,insize,outsize],initializer) (2)
biases=tf.get_varialble('biases',[outsize],initializer) (3)
式(3-1),(3-2)中filter_weight与biases分别为通过tf.get_varialble方法创建的卷积层中滤波器的权重向量和阈值向量。filter_weight中包含一个4维矩阵的参数,前二个维度为滤波器的尺寸(长宽),后两个维度为当前层的深度和滤波器的深度;biases中包含的参数为与权重向量个数相匹配的阈值向量,initializer为权重阈值的初始化策略。TensorFlow提供conv2d函数来实现卷积层的核心算法:
tf.nn.conv2d(input,filter_weight,strides,padding) (4)
其中input为接收上一层的输出,即本层的输入,filter_weight提供了当前卷积层的权重参数,strides为滤波器在卷积过程中的步长参数,padding为TensorFlow提供的卷积步长策略。
为了减少模型中的参数,并降低计算量,在卷积层后添加池化层,对卷积层输出得到的特征映射做抽样或聚合,减少计算量,避免过拟合。TensorFlow提供以下两种池化策略函数tf.nn.max_pool和tf.nn.avg_pool:
tf.nn.max_pool(input,ksize,strides,padding) (5)
tf.nn.avg_pool(input,ksize,strides,padding) (6)
input为上层的输出,即为池化层的输入,ksize为池化窗口大小参数,strides为池化窗口移动步长,padding为池化操作后大小保留的策略。
激活函数的选择的对于整个卷积神经网络的预测结果有着很大的影响,选择合适的激活函数可以提高网络模型的收敛速度。本发明使用矫正的线性激活函数Relu函数,当神经元输出的值小于0时恒为0,否则保持原来的值不变。这样,激活函数导数大于0时恒为1,可以将梯度很好的传递到上一层,避免了梯度消散的问题,加快了误差反向传播,极大的提高了训练的速度。
(5)流域的洪水预测模型是指适应某具体流域的洪水预测模型,通过模型训练过
程中,从输入参数中立定的模型参数体现了当前流域的特征。
(6)建立误差校正模型BPCM训练数据集。其输入与步骤(3)中的模型训练数据集输入相同,输出为预测模型CNNM的预测误差。假设对输入 模型预测输出为误差为校正模型的输入输出数据集中每个输入输出对为
其中,
(7)建立基于BP神经网络的误差校正模型BPCM,采用(5)产生的训练数据,训练BPCM,建立输入数据和预测模型预测的误差之间的非线性关系。
(8)流域实时误差校正模型是指适应具体流域的BP校正模型。该模型中的参数体
现了当前流域的特征。
(9)实时监测降雨数据是指当前预测点之前的实时监测的降雨数据和流量数据等输入数据。
(10)利用流域洪水预测模型CNNFM,以实时监测降雨数据为输入,进行预测得到初始预测值。
(11)利用流域实时误差校正模型BPCM,以实时监测数据为输入,对初始预测值进行校正。将校正误差作用于CNNFM预测模型的结果predValue,则最终预测结果为predValue-e。采用纳什效率系数(Nash-Sutcliffe efficiency coefficient,NSE)来评定预测结果的质量。
(12)评价预测结果
1)纳什效率系数
对水文过程模拟的效果检验使用纳什效率系数(Nash-Sutcliffe efficiencycoefficient,NSE)来评定,NSE一般用来验证水文过程模拟结果的好坏。NSE的取值范围为负无穷到1,NSE的值越接近1是模型预测的质量越好,模型可信度越高;NSE的值越接近0,表示模型预测的结果接近平均值水平,预测结果总体可信,但存在较大误差;NSE的值远远小于0,预测结果是不可信的。NSE计算方法如式7)。
2)均方根误差(RMSE)
实施例:
为了验证本发明的效果,选取安徽省屯溪流域,建立模型预测该流域内的洪水过程,流域集水面积2696.76平方公里,地处亚热带季风气候,气候适宜,年平均温度17℃,雨量充沛。选取从1982年到2002年间发生的33场洪水数据作为实验的研究资料。所有数据都经过水文数据整编处理,不存在缺失数据。
下面从预见期、实时校正结果两个方面,进行对比实验,分析实验结果并检测模型的可用性。
1)从实验模型的预见期出发。以预见期为出发点设计对比实验,进行流量预报,设置预见期分别为1小时、2小时、4小时的卷积神经网络流量预测实验,使用NSE、RMSE、MAE三个评价指标对模型的预测结果进行分析。
2)从对预测结果的实时校正方面出发。将预测结果结合基于得分加权的组合校正模型来进行实时校正,结合预测结果与NSE、RMSE、MAE评价指标,测试未经过实时校正的预报模型的性能,与未经过实时校正的预报模型的测试结果进行对比。
1.数据准备
本节以屯溪流域(包括11个报讯雨量站、1个蒸发站,如图3)降雨量及流量数据作为实验数据,选取从1982年到2002年间发生的33场洪水数据作为实验的研究资料。
表3.1水雨情报汛站网
屯溪流域位于我国东南沿海[56],流域集水面积2696.76平方公里,地处亚热带季风气候,气候适宜,年平均温度17℃,雨量充沛,年降水量1600毫米,年内降水分配不均,其中4-6月降水较为充沛,易发生洪涝灾害;7-9月降水较少,旱灾频繁。部分降雨及流量数据如表1。
表11989年部分降雨流量表
实验数据为屯溪流域流量雨量数据,其中1982年~1994年间23场洪水数据作为训练集,1995年~2002年间10场洪水数据作为测试集。
历史洪水数据在水文数据整编过程进行了异常和缺失处理,因此本实验无需要处理。我们首先对历史水文数据进行皮尔逊相关系数分析得到预测输入时间段,设当前时刻为t,选取t-14时刻到t时刻屯溪流量与11个雨量站的的数据进行皮尔逊相关性分析。
由上表可以看出,t时刻屯溪流量与t-8时刻流量之间,随着时间越来越接近t时刻,变量之间的相关性越高;其余11个雨量站与屯溪流量之间的相关性随着时间越来越接近t时刻呈现出变量相关性先单调递增再单调递减的趋势,综合来看,选取当前时刻t时前8个时刻的历史信息作为输入数据。
(3)采用min-max标准化,也称为离差标准化,它将原始数据进行线性变化原始值被规约到[0,1]之间。
以下文中预见期为1小时为例,使用滑动窗口技术对训练和测试数据进行滑动切分,将预见期前8小时流量及雨量数据作为输入数据,后1小时流量数据作为输出数据建立预测模型。训练样本数据和测试样本数据(归一化前)如表2。
表2屯溪雨量流量1h预见期数据集
以上训练数据为预见期为1小时时实验数据,hisData表示预见期前8小时的流量雨量数据,predLiuliang1为预见期预测的流量值。当预见期为2小时、4小时时,上述训练集格式需作相应改动。
2.CNNFM模型建立
CNNFM模型输入层-输出层参数确定:由以上分析,输入层为预见期前8h流量雨量数据,总计共12×8个神经元节点,输出层为预见期大小(假设为4h,共4个神经元)。
CNNFM模型卷积层参数确定:其中滤波器中初始权重向量初始化策略采取TensorFlow中提供的tf.truncated_normal,它从截断的正太分布中获取随机值;阈值向量采取tf.Constant(0.1),将阈值向量初始化为0.1。共设置3个卷积层,节点个数分别为128、256、512,设置卷积层参数为式(9),其中filter_weight卷积核大小为3×3,strides滑动步长为1。
tf.nn.conv2d(x_input,filter_weight,strides=[1,1,1,1]) (9)
CNNFM模型池化层参数确定:考虑到屯溪流域降雨充沛,但年内降水分配不均的特点,且数据中多存在降雨为0mm的时刻数据,为保留数据的纹理信息,在减少参数的同时,保留更多的特征,采用maxpooling最大值池化策略来设计池化层,在每一个卷积层后添加一个池化层,池化函数为式(10),其中池化窗口大小为2*2,滑动步长为2。
tf.nn.max_pool(x_input,ksize=[1,2,2,1],strides=[1,2,2,1]) (10)
CNNFM模型激活函数及其他参数确定:使用Relu激活函数对卷积池化输出结果进行保留映射,在其后添加1024个神经元节点的全连接层,使其连接到相应数目的输出节点。设置学习率时使用tf.train.AdamOptimizer自适应学习率优化算法,不需要人为调整学习率,使用默认参数就可能获得最优值。训练总轮数设置为100000。
3.实验分析
1)不带实时校正的模型在不同预见期内准确率对比
对预见期为1h、2h、4h进行基于卷积神经网络的洪峰预测,选择预报的最后一小时预报结果作为表格统计结果,测试结果如表3。
表3 1h、2h、4h预见期洪水模拟特征值
以测试数据中1997060608号洪水预报为例,流量过程预测如图4。
对预见期为1h、2h、4h进行基于卷积神经网络流量过程预测,计算得到NSE、RMSE、MAE这三个指标,结果如表2。
表2不同预见期洪水过程预测实验结果评价指标表
根据实验结果所示,对于测试数据中10场洪水过程而言,卷积神经网络模型预测洪峰在三种预见期下平均相对误差分别为6.5%、10.2%与14.14%;纳什效率系数均值分别为97.71%、94.57%与89.17%;在1998050108号和1999062215号洪水峰现时间上存在较大误差,总体预报结果良好。由应用结果可以看出,RMSE与MAE均随预见期延长而提高,总体看来积神经网络模型能够很好的模拟屯溪流域的洪水过程,并且显示出随预见期延长,预报误差逐渐增大的趋势。按照《水文情报预报规范》(GB/T 22482-2008),卷积神经网络模型在1h、2h、4h预见期下水文过程模拟结果分别达到了甲级、甲级、乙级标准。
2)卷积神经网络预报结合实时校正的分析
本节针对预测结果从带有实时校正以及不带实时校正的结果来分别展示预测结果,以预见期中最后一小时预测结果进行实验分析。
1.使用BP神经网络模型的校正模型BPCM来对屯溪流域预测数据进行实时校正,已知测试数据集为1995年~2002年间共计10场洪水数据。
以2001050108号洪水数据为例,预见期为1h时,预测结果以及校正结果如图5。
预见期为2h时,预测结果及校正结果如图6。
预见期为4h时,预测结果及校正结果如图7。
2.对预见期为1h、2h、4h进行基于卷积神经网络流量预测,计算得到NSE、RMSE、MAE这三个指标,结果如表4。
表4结合实时校正的洪水过程预测实验结果评价指标表
根据实验结果所示,随着预见期延长,洪水预报结果精度逐渐下降,并且在对预测结果进行校正之后,与不带校正的预测结果相比,预见期为1h时,RMSE、MAE有小幅度缩减,NSE有提高,经过校正模型校正效果较好,三种预见期下纳什效率系数均超过0.9,按照《水文情报预报规范》(GB/T 22482-2008),卷积神经网络模型在1h、2h、4h预见期下水文过程模拟结果均达到了甲级标准。

Claims (7)

1.一种基于深度学习模型和BP神经网络校正的洪水预测方法,其特征在于,包括如下步骤:
(1)对历史洪水过程数据进行归一化处理;
(2)对归一化后的历史洪水过程数据序列进行分析,分析降雨量和蒸发量影响因素对流域出口流量的影响时间范围;
(3)利用滑动窗口从历史场次洪水数据中建立预测模型的输入和输出值,建立模型训练数据集TRSet1;
(4)建立基于深度学习的洪水预测模型CNNFM,利用训练数据TRSet1训练预测模型得到模型参数,建立流域的洪水预测模型;
(5)流域的洪水预测模型是指适应某具体流域的洪水预测模型,通过模型训练过程中,从输入参数中立定的模型参数体现了当前流域的特征。
(6)建立实时误差校正模型训练数据集TRSet2,输入与TRSet1相同,流域洪水预测模型的预测误差作为输出;
(7)建立基于BP神经网络的误差校正模型BPCM,利用TRSet2训练校正模型,获得模型参数,建立流域实时误差校正模型;
(8)流域实时误差校正模型是指适应具体流域的BP校正模型。该模型中的参数体现了当前流域的特征。
(9)对实时监测的降雨量和流量,利用流域洪水预测模型CNNFM对实时数据进行预测,利用流域实时误差校正模型BPCM进行校正,获得最终的预测值。
2.根据权利要求1所述的一种基于深度学习模型和BP神经网络校正的洪水预测方法,其特征在于,所述步骤(1)中对历史洪水过程数据进行归一化处理的具体步骤如下:
采用离差标准化即min-max标准化将洪水过程预测涉及到的降雨量、蒸发量、流量和水位数据进行线性变化,通过转化公式将原始值映射到[0,1]之间;转化公式如下:
其中x*为转换后的值,x为原始值,xmin=min(x),xmax=max(x);经过min-max标准化后,原始值被规约到[0,1]之间;其中,max(x)取x对应指标的历史最大值,min(x)取x对应指标的历史最小值。
3.根据权利要求1所述的一种基于深度学习模型和BP神经网络校正的洪水预测方法,其特征在于,所述步骤(2)中分析降雨量和蒸发量影响因素对流域出口流量的影响时间范围;的具体步骤如下:
对归一化后的历史洪水过程数据序列采用皮尔逊相关系数分析输出量即流域出口流量与不同时间范围的各输入量即流域内各雨量站的降雨量和蒸发站的蒸发量的相关性,确定输入量对输出量影响的时间范围,各输入量的影响范围的最大值作为整体输入量对输出量影响的时间范围。
4.根据权利要求1所述的一种基于深度学习模型和BP神经网络校正的洪水预测方法,其特征在于,所述步骤(3)中建立模型训练数据集TRSet1的具体步骤如下:
利用步骤(2)中得到的时间范围作为滑动窗口的宽度,按照步长为1提取历史洪水过程数据作为输入,对应输出为预测预见期为宽度的输出序列;具体的输入量对输出量影响的最大时间范围为w,即表示某个输入量的t0时刻监测值,从t1到tw时间内,逐渐影响输出量,并最终影响消失;
上式表示预测模型的输入为I的矩阵,有m+1个输入量,包括预测量自身;输出为O的向量,预见期为n。
5.根据权利要求1所述的一种基于深度学习模型和BP神经网络校正的洪水预测方法,其特征在于,所述步骤(4)中建立基于深度学习的洪水预测模型CNNFM,采用卷积神经网络构建,其中预测模型CNNFM包括卷积层、池化层和激活函数。
6.根据权利要求1所述的一种基于深度学习模型和BP神经网络校正的洪水预测方法,其特征在于,所述步骤(5)中建立误差校正模型训练数据集TRSet2的具体步骤如下:
输入与步骤(3)中的预测模型数据集TRSet1输入相同,输出为预测模型CNNFM的预测误差,其中对输入为模型预测输出为 误差为校正模型的输入输出数据集中每个输入输出对为其中,
其中,I表示输入,O表示输出。
7.根据权利要求1所述的一种基于深度学习模型和BP神经网络校正的洪水预测方法,其特征在于,所述步骤(7)中训练校正模型,具体步骤如下:
建立BP神经网络校正模型BPCM,利用TRSet2训练校正模型,通过BP神经网络建立输入数据和输出误差之间的非线性关系,获得模型参数。
CN201810965624.5A 2018-08-23 2018-08-23 一种基于深度学习模型和bp神经网络校正的洪水预测方法 Active CN109272146B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810965624.5A CN109272146B (zh) 2018-08-23 2018-08-23 一种基于深度学习模型和bp神经网络校正的洪水预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810965624.5A CN109272146B (zh) 2018-08-23 2018-08-23 一种基于深度学习模型和bp神经网络校正的洪水预测方法

Publications (2)

Publication Number Publication Date
CN109272146A true CN109272146A (zh) 2019-01-25
CN109272146B CN109272146B (zh) 2021-10-19

Family

ID=65154150

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810965624.5A Active CN109272146B (zh) 2018-08-23 2018-08-23 一种基于深度学习模型和bp神经网络校正的洪水预测方法

Country Status (1)

Country Link
CN (1) CN109272146B (zh)

Cited By (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109978235A (zh) * 2019-03-04 2019-07-05 宁波市气象服务中心 一种基于样本学习的积涝水位预测方法
CN109993354A (zh) * 2019-03-24 2019-07-09 北京工业大学 一种针对能源消费量预测的方法
CN110163419A (zh) * 2019-04-28 2019-08-23 河海大学 一种中小河流流域洪水预报的方法
CN110414713A (zh) * 2019-06-27 2019-11-05 电子科技大学 一种基于同步数据流压缩的径流实时预测方法
CN110728353A (zh) * 2019-09-10 2020-01-24 宁波市气象台 一种基于bp人工神经网络的城市积涝水位预报方法
CN110969282A (zh) * 2019-10-17 2020-04-07 天津大学 一种基于lstm复合网络的径流量稳定预测方法
CN111027193A (zh) * 2019-12-02 2020-04-17 浙江工业大学 一种基于回归模型的短期水位预测方法
CN111160472A (zh) * 2019-12-30 2020-05-15 深圳集智数字科技有限公司 对象的目标特征的预测方法、装置、存储介质与电子设备
CN111666656A (zh) * 2020-05-09 2020-09-15 江苏微之润智能技术有限公司 基于微波雨衰的降雨估算方法及降雨监测系统
CN111832810A (zh) * 2020-06-23 2020-10-27 西安电子科技大学 洪水预报模型、信息处理方法、存储介质、计算机设备
CN112039711A (zh) * 2020-09-08 2020-12-04 中国联合网络通信集团有限公司 流量预测方法及设备
CN112084461A (zh) * 2020-08-10 2020-12-15 河海大学 一种基于历史洪水学习的knn实时校正方法
CN112132333A (zh) * 2020-09-16 2020-12-25 安徽泽众安全科技有限公司 一种基于深度学习的短期水质水量预测方法及系统
CN112149892A (zh) * 2020-09-17 2020-12-29 中国水利水电科学研究院 一种洪水重现期输出方法和装置
CN112418542A (zh) * 2020-12-03 2021-02-26 浙江知水信息技术有限公司 基于气象数据利用机器深度学习实现汛情提前预警的方法
CN112525438A (zh) * 2020-10-15 2021-03-19 国网浙江省电力有限公司杭州供电公司 一种sf6密度继电器的漏气监测方法及系统
CN112560338A (zh) * 2020-12-10 2021-03-26 东北大学 基于自适应深度学习的复杂工业系统智能预报方法、装置、设备及存储介质
CN112631215A (zh) * 2020-12-10 2021-04-09 东北大学 工业过程运行指标智能预报方法、装置、设备及存储介质
CN112950062A (zh) * 2021-03-24 2021-06-11 郑州大学 基于动态系统响应曲线和lstm的山洪预报实时校正方法
CN112990108A (zh) * 2021-04-19 2021-06-18 四川省水利科学研究院 一种基于卷积神经网络实现堤坝护坡的系统
CN113139329A (zh) * 2021-03-04 2021-07-20 河海大学 一种基于水文相似性和人工神经网络的新安江模型参数率定方法
CN113566931A (zh) * 2021-07-22 2021-10-29 水利部南京水利水文自动化研究所 基于边缘计算的闸前反射式水位计智能校准方法及系统
CN113705931A (zh) * 2021-09-17 2021-11-26 中国长江电力股份有限公司 一种利用k最邻近法预测径流要素的方法
CN114386693A (zh) * 2022-01-11 2022-04-22 福建中锐网络股份有限公司 基于PSO-LSTM-Attention模型的水库水位预测预警方法
CN116011678A (zh) * 2023-02-13 2023-04-25 天津大学 一种耦合峰值修正的洪水过程自动寻优预报方法
CN117973236A (zh) * 2024-04-01 2024-05-03 北京适创科技有限公司 一种对仿真结果的校正方法、系统、设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102880755A (zh) * 2012-09-25 2013-01-16 河海大学 一种极端降水的定量预测方法及系统
CN104239489A (zh) * 2014-09-05 2014-12-24 河海大学 利用相似性搜索和改进bp神经网络预测水位的方法
CN106650767A (zh) * 2016-09-20 2017-05-10 河海大学 基于聚类分析和实时校正的洪水预报方法
CN107203810A (zh) * 2017-05-22 2017-09-26 河海大学 一种基于深度网络的降水量预测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102880755A (zh) * 2012-09-25 2013-01-16 河海大学 一种极端降水的定量预测方法及系统
CN104239489A (zh) * 2014-09-05 2014-12-24 河海大学 利用相似性搜索和改进bp神经网络预测水位的方法
CN106650767A (zh) * 2016-09-20 2017-05-10 河海大学 基于聚类分析和实时校正的洪水预报方法
CN107203810A (zh) * 2017-05-22 2017-09-26 河海大学 一种基于深度网络的降水量预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YUN BAI等: "Daily reservoir inflow forecasting using multiscale deep feature learning with hybrid models", 《JOURNAL OF HYDROLOGY》 *
马凯凯等: "数据驱动的中小河流智能洪水预报方法对比研究", 《中国科学技术大学学报》 *

Cited By (38)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109978235A (zh) * 2019-03-04 2019-07-05 宁波市气象服务中心 一种基于样本学习的积涝水位预测方法
CN109978235B (zh) * 2019-03-04 2022-12-27 宁波市气象服务中心 一种基于样本学习的积涝水位预测方法
CN109993354A (zh) * 2019-03-24 2019-07-09 北京工业大学 一种针对能源消费量预测的方法
CN110163419A (zh) * 2019-04-28 2019-08-23 河海大学 一种中小河流流域洪水预报的方法
CN110163419B (zh) * 2019-04-28 2022-08-19 河海大学 一种中小河流流域洪水预报的方法
CN110414713A (zh) * 2019-06-27 2019-11-05 电子科技大学 一种基于同步数据流压缩的径流实时预测方法
CN110728353A (zh) * 2019-09-10 2020-01-24 宁波市气象台 一种基于bp人工神经网络的城市积涝水位预报方法
CN110969282A (zh) * 2019-10-17 2020-04-07 天津大学 一种基于lstm复合网络的径流量稳定预测方法
CN111027193A (zh) * 2019-12-02 2020-04-17 浙江工业大学 一种基于回归模型的短期水位预测方法
CN111160472A (zh) * 2019-12-30 2020-05-15 深圳集智数字科技有限公司 对象的目标特征的预测方法、装置、存储介质与电子设备
CN111666656A (zh) * 2020-05-09 2020-09-15 江苏微之润智能技术有限公司 基于微波雨衰的降雨估算方法及降雨监测系统
CN111832810A (zh) * 2020-06-23 2020-10-27 西安电子科技大学 洪水预报模型、信息处理方法、存储介质、计算机设备
CN111832810B (zh) * 2020-06-23 2023-06-27 西安电子科技大学 洪水预报模型、信息处理方法、存储介质、计算机设备
CN112084461B (zh) * 2020-08-10 2024-05-14 河海大学 一种基于历史洪水学习的knn实时校正方法
CN112084461A (zh) * 2020-08-10 2020-12-15 河海大学 一种基于历史洪水学习的knn实时校正方法
CN112039711B (zh) * 2020-09-08 2023-03-24 中国联合网络通信集团有限公司 流量预测方法及设备
CN112039711A (zh) * 2020-09-08 2020-12-04 中国联合网络通信集团有限公司 流量预测方法及设备
CN112132333B (zh) * 2020-09-16 2024-02-02 安徽泽众安全科技有限公司 一种基于深度学习的短期水质水量预测方法及系统
CN112132333A (zh) * 2020-09-16 2020-12-25 安徽泽众安全科技有限公司 一种基于深度学习的短期水质水量预测方法及系统
CN112149892A (zh) * 2020-09-17 2020-12-29 中国水利水电科学研究院 一种洪水重现期输出方法和装置
CN112525438A (zh) * 2020-10-15 2021-03-19 国网浙江省电力有限公司杭州供电公司 一种sf6密度继电器的漏气监测方法及系统
CN112418542A (zh) * 2020-12-03 2021-02-26 浙江知水信息技术有限公司 基于气象数据利用机器深度学习实现汛情提前预警的方法
CN112560338A (zh) * 2020-12-10 2021-03-26 东北大学 基于自适应深度学习的复杂工业系统智能预报方法、装置、设备及存储介质
CN112631215A (zh) * 2020-12-10 2021-04-09 东北大学 工业过程运行指标智能预报方法、装置、设备及存储介质
CN112560338B (zh) * 2020-12-10 2022-03-25 东北大学 基于自适应深度学习的复杂工业系统智能预报方法、装置、设备及存储介质
WO2022121932A1 (zh) * 2020-12-10 2022-06-16 东北大学 基于自适应深度学习的复杂工业系统智能预报方法、装置、设备及存储介质
CN113139329B (zh) * 2021-03-04 2022-10-18 河海大学 一种基于水文相似性和人工神经网络的新安江模型参数率定方法
CN113139329A (zh) * 2021-03-04 2021-07-20 河海大学 一种基于水文相似性和人工神经网络的新安江模型参数率定方法
CN112950062B (zh) * 2021-03-24 2023-02-03 郑州大学 基于动态系统响应曲线和lstm的山洪预报实时校正方法
CN112950062A (zh) * 2021-03-24 2021-06-11 郑州大学 基于动态系统响应曲线和lstm的山洪预报实时校正方法
CN112990108A (zh) * 2021-04-19 2021-06-18 四川省水利科学研究院 一种基于卷积神经网络实现堤坝护坡的系统
CN113566931A (zh) * 2021-07-22 2021-10-29 水利部南京水利水文自动化研究所 基于边缘计算的闸前反射式水位计智能校准方法及系统
CN113566931B (zh) * 2021-07-22 2023-06-09 水利部南京水利水文自动化研究所 基于边缘计算的闸前反射式水位计智能校准方法及系统
CN113705931A (zh) * 2021-09-17 2021-11-26 中国长江电力股份有限公司 一种利用k最邻近法预测径流要素的方法
CN113705931B (zh) * 2021-09-17 2023-04-18 中国长江电力股份有限公司 一种利用k最邻近法预测径流要素的方法
CN114386693A (zh) * 2022-01-11 2022-04-22 福建中锐网络股份有限公司 基于PSO-LSTM-Attention模型的水库水位预测预警方法
CN116011678A (zh) * 2023-02-13 2023-04-25 天津大学 一种耦合峰值修正的洪水过程自动寻优预报方法
CN117973236A (zh) * 2024-04-01 2024-05-03 北京适创科技有限公司 一种对仿真结果的校正方法、系统、设备及存储介质

Also Published As

Publication number Publication date
CN109272146B (zh) 2021-10-19

Similar Documents

Publication Publication Date Title
CN109272146A (zh) 一种基于深度学习模型和bp神经网络校正的洪水预测方法
CN109299812A (zh) 一种基于深度学习模型和knn实时校正的洪水预测方法
CN106650767B (zh) 基于聚类分析和实时校正的洪水预报方法
CN114254561B (zh) 一种内涝预测方法、系统及存储介质
Wilby et al. Detection of conceptual model rainfall—runoff processes inside an artificial neural network
CN104715292A (zh) 基于最小二乘支持向量机模型的城市短期用水量预测方法
CN107526083B (zh) 一种基于天气雷达数据的强对流风力等级预测方法
CN107145965A (zh) 一种基于相似度匹配和极限学习机的河流洪水预测方法
CN108898251A (zh) 考虑气象相似性和功率波动的海上风电场功率预测方法
CN108428017A (zh) 基于核极限学习机分位数回归的风电功率区间预测方法
CN104090974B (zh) 展延水库后续来水的动态数据挖掘方法及系统
CN110363349A (zh) 一种基于ascs的lstm神经网络水文预测方法及系统
CN107705556A (zh) 一种基于支持向量机和bp神经网络结合的交通流预测方法
CN114970377B (zh) 基于新安江与深度学习耦合模型的场次洪水预报方法和系统
CN106875048A (zh) 基于典型经验模式的应急洪水预报方法
CN108021773A (zh) 基于dss数据库读写的分布式水文模型多场次洪水参数率定方法
CN109492748A (zh) 一种基于卷积神经网络的电力系统的中长期负荷预测模型建立方法
CN110837915A (zh) 一种基于混合集成深度学习的电力系统低压负荷点预测及概率预测方法
CN105005708A (zh) 一种基于ap聚类算法的广义负荷特性聚类方法
CN110276477A (zh) 一种基于分层贝叶斯网络和增量学习的洪水预报方法
CN115186857A (zh) 一种基于集成学习的神经网络水库水位预测方法
CN111143498B (zh) 小河流洪水预报方法
CN116432828A (zh) 一种缺资料流域径流量智能预测方法
CN109063388A (zh) 风环境模拟的微气候建筑选址设计方法
CN106570582A (zh) 一种用于建立输电线路舞动跳闸风险预测的网络模型的方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant