CN112182709B - 大型水库叠梁门分层取水设施的下泄水温快速预测方法 - Google Patents
大型水库叠梁门分层取水设施的下泄水温快速预测方法 Download PDFInfo
- Publication number
- CN112182709B CN112182709B CN202011040789.5A CN202011040789A CN112182709B CN 112182709 B CN112182709 B CN 112182709B CN 202011040789 A CN202011040789 A CN 202011040789A CN 112182709 B CN112182709 B CN 112182709B
- Authority
- CN
- China
- Prior art keywords
- model
- reservoir
- water
- hidden layer
- training
- 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.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 125
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012549 training Methods 0.000 claims abstract description 56
- 238000012360 testing method Methods 0.000 claims abstract description 44
- 238000012216 screening Methods 0.000 claims abstract description 12
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 9
- 230000007787 long-term memory Effects 0.000 claims abstract description 7
- 230000006403 short-term memory Effects 0.000 claims abstract description 7
- 238000013135 deep learning Methods 0.000 claims abstract description 6
- 238000009960 carding Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 26
- 230000006870 function Effects 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 230000015654 memory Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 238000010206 sensitivity analysis Methods 0.000 claims description 7
- 230000004913 activation Effects 0.000 claims description 6
- 239000003657 drainage water Substances 0.000 claims description 6
- 238000012821 model calculation Methods 0.000 claims description 6
- 230000006872 improvement Effects 0.000 claims description 5
- 230000007246 mechanism Effects 0.000 claims description 4
- 230000008020 evaporation Effects 0.000 claims description 3
- 238000001704 evaporation Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000001556 precipitation Methods 0.000 claims description 3
- 238000007599 discharging Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 6
- 238000010276 construction Methods 0.000 description 4
- 230000002411 adverse Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005842 biochemical reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013136 deep learning model Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/048—Activation functions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/049—Temporal neural networks, e.g. delay elements, oscillating neurons or pulsed inputs
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Mathematical Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- Data Mining & Analysis (AREA)
- Computer Hardware Design (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了大型水库叠梁门分层取水设施的下泄水温快速预测方法,包括以下步骤:S1、梳理影响水库下泄水温的主要因素;S2、初步搭建水库下泄水温预测模型框架;S3、形成数据集S1;S4、选出模型的输入因子集X2;S5、对数据集S2做标准化处理;S6、构建基于长短期记忆网络深度学习算法的大型水库叠梁门分层取水下泄水温预测模型;S7、将数据集S2分为训练样本集和测试样本集;S8、采用长短期记忆网络算法和Holdout训练方法训练并生成水库下泄水温预测模型;S9、筛选出最优隐层层数、训练批量、最大迭代次数和隐层节点数;S10、确定水库下泄水温预测模型,并得到下泄水温预测结果。
Description
技术领域
本发明属于水库水温预测领域,具体涉及一种大型水库叠梁门分层取水设施的下泄水温快速预测方法。
背景技术
高坝大库建成蓄水后,坝上河道水位抬升形成水库,改变了原始河流水体的水动力和热动力特征,造成库区水体水温态势的改变,不仅影响水体的物理化学特性和生化反应速率,而且关系到水生生物的生长繁殖和水生态***的结构。为减缓水库建成后水温变化对下游生态环境的影响,需要采取相应的改善措施,分层取水是减缓水温分层不利生态环境影响的重要工程措施。
2000年以来,我国有近20座大型水电工程采取叠梁门分层取水措施,以减缓高坝大库下泄低温水对水生生物和农作物的不利影响。目前这些工程主要依托叠梁门分层取水措施,开展水库生态调度,通过改变取水口位置和水库的径流过程,改变大坝下游水体的水动力和热动力特性,调节下泄水体的温度,从而减缓水电工程的不利生态环境影响。
但是,目前水库分层取水设施在实际应用中尚缺乏科学有效的指导。现阶段分层取水设施的运行多是基于水库水位条件或是电站设计阶段、环评阶段提出的运行方案来开展调度工作,尚存在以下问题:
(1)仅根据水位条件来开展分层取水调度,难以引取到目标水温层水体;
(2)电站设计阶段或者是环评阶段提出的运行方案,多是利用数值模拟模型针对典型年工况设计而成,与电站面临的实际工况场景差异较大;
(3)现有数值模拟模型计算耗时较长,构建模型构建专业性强,构建过程复杂,对使用人员的专业性要求比较高,因此难以根据水库面临的工况场景,快速给出分层取水设施的调度运行方案,不适合直接指导分层取水设施的运行管理。
发明内容
发明目的:为解决现有技术的不足,本发明提供了一种大型水库叠梁门分层取水设施的下泄水温快速预测方法。
技术方案:大型水库叠梁门分层取水设施的下泄水温快速预测方法,包括以下步骤:
S1、基于对水库水温结构形成物理机理的分析,梳理影响水库下泄水温的主要因素;
S2、将步骤S1得到的影响水库下泄水温的主要因素纳入模型的输入因子集X1,以所述下泄水温为模型输出因子Y,初步搭建水库下泄水温预测模型框架,经初步筛选,模型的输入因子集X1包括当前时段所处的月份M、入流量Qi t、出流量QO t、入流水温Ti t、坝上水位Zt、取水口深度Dt、叠梁门运行层数Lt,以及水库当前时段气温Ta t、气压Pt、湿度Ht、太阳辐照度At、风速Wt、云量Ct、蒸发量Et、降水量Rt,预测变量即输出因子Y为水库当前时段的下泄水温To t;
S3、采集所述输入因子集X1和输出因子Y的长时间序列资料,整理形成数据集S1,则所述数据集S1={X1,Y}={M、t、Qi t、Qo t、Ti t、Zt、Dt、Lt、Ta t、Pt、Ht、At、Wt、Ct、Et、Rt、To t},其中t为当前时刻在日内所处时段,t∈[0,24];
S4、基于所述数据集S1,利用敏感性分析及逐步测试,对输入因子进行二次筛选,最终选出模型的输入因子集X2;
S5、对数据集S2做标准化处理;
S6、以经过步骤S5标准化处理后的影响因子数据X2为模型输入向量,以所述标准化处理后的预测变量Y为模型输出向量,构建基于长短期记忆网络(Long and Short-TermMemory network,LSTM)深度学习算法的大型水库叠梁门分层取水下泄水温预测模型,并根据所述标准化处理后的数据集确定训练样本集和测试样本集,其中,所述水库叠梁门分层取水下泄水温预测模型的函数表达式为To t=f(X2),f表示LSTM计算模块;
S7、将所述数据集S2划分为训练样本集和测试样本集;
S8、设定所述水温预测模型的隐层层数、训练批量、隐层节点数和最大迭代次数,结合训练样本集,采用长短期记忆网络算法和Holdout训练方法训练并生成水库下泄水温预测模型;
S9、利用测试样本集对所述水库下泄水温预测模型能进行测试,测试不同隐层层数、训练批量、隐层节点数和最大迭代次数的结构参数组合对所述水库下泄水温预测模型精度、收敛情况和计算耗时的影响,根据测试结果,在保证水库下泄水温预测模型收敛的基础上,重点考虑模型精度,结合模型计算耗时,筛选出最优隐层层数、训练批量、最大迭代次数和隐层节点数;
S10、根据所述最优隐层层数、最优训练批量、最优迭代次数和最优隐层节点数,确定水库下泄水温预测模型,并得到下泄水温预测结果。
进一步地,步骤S6中的LSTM计算模块包括一个输入门,一个记忆细胞,一个遗忘门和一个输出门,计算过程如下:
ft=σ(Wxfxt+Whfht-1+bf)
it=σ(Wxixt+Whiht-1+bi)
Ct=ft·Ct-1+it·tanh(Wxcxt+Whcht-1+bc)
ot=σ(Wxoxt+Whoht-1+bo)
ht=ot·tanh(Ct)
ft、it、Ct、ot和ht分别表示遗忘门、输入门、细胞状态、输出门和隐含层的记忆状态或输出;
Whi表示权重系数矩阵;
bi表示偏置项,
σ为Sigmoid激活函数;
tanh为双曲正切激活函数。
进一步地,步骤S1中影响水库下泄水温的主要因素包括水库调度方式、叠梁门运行方式、入流水温、气象条件。
进一步地,步骤S4包括:以下步骤:
S41、利用敏感性分析,计算各输入因子对输出因子Y的影响程度,将输入因子集分为高影响度输入因子集X1h和低影响度输入因子集X1s;
S42、高影响度输入因子集X1h直接选定为模型输入因子,低影响度输入因子集X1s通过逐步测试,逐一加入到模型输入因子,若加入后对模型精度的提高度5%,则加入模型输入因子,反之则不加入;最终筛选出的模型输入因子集记为X2;,最终整理形成的数据集记为S2,则S2={X2,Y}。
进一步地,步骤S7中将所述数据集S2确定训练样本和测试样本包括以下步骤:
基于二八定律,将数据集划分为训练样本和测试样本,其中,所述训练样本约占所述数据集总年份的80%,所述测试样本约占所述数据集中总年份的20%。
进一步地,步骤S9中测试结构参数组合对模型性能影响测试进一步包含以下步骤:
在统一设定其他参数不同的基础上,依次测试隐层数、训练批量、迭代次数和隐层节点数对模型精度、收敛情况和计算速度的影响,在保证模型收敛的基础上,重点考虑模型精度,结合模型计算耗时,筛选出最优隐层层数、训练批量、最大迭代次数和隐层节点数设定组合。
本发明公开的大型水库叠梁门分层取水设施的下泄水温快速预测方法,具有以下有益效果:
(1)该方法从水库水温结构形成机理出发,初步筛选了下泄水温影响因素,而后通过敏感性分析和逐步筛选,最终选定了模型的输入变量,通过两次优化筛选,在赋予AI模型以物理意义的同时,又从数理统计的角度对模型的输入变量进行把控,为模型预测精度的优化提供了双重保证;
(2)该方法基于深度学习领域中的LSTM计算模块从水库历史数据中学习水库下泄水温影响因子与下泄水温之间的映射关系,相比于传统反向传播神经网络(BP)、支持向量回归(SVR)等机器学习模型,深度学习模型具备更强的学习能力,能够基于对水库历史数据的学习,给出更为精确的下泄水温预测结果,改善预测效果;
(3)本发明提供的基于深度学习的大型水库叠梁门分层取水下泄水温预测模型,模型计算耗时较短,通过参数调整,一般能在几分钟内完成训练过程,在几秒钟之内完成预测过程,远远低于传统数值模拟模型,能够结合水库面临场景,快速给出水库下泄水温预测结果,保证了模型的实用性。
附图说明
图1为本发明公开的大型水库叠梁门分层取水设施的下泄水温快速预测方法的流程示意图;
图2a和图2b为本发明实施例提供的水库下泄水温预测模型的隐层层数与预测精度平均绝对误差MAE、计算耗时之间关系的示意图;
图3a和图3b为本发明实施例提供的水库下泄水温预测模型的训练批量与预测精度平均绝对误差MAE、计算耗时之间关系的示意图;
图4a和图4b为本发明实施例提供的水库下泄水温预测模型的隐层节点数与预测精度平均绝对误差MAE、计算耗时之间关系的示意图;
图5a和图5b为本发明实施例提供的水库下泄水温预测模型的迭代次数与预测精度平均绝对误差MAE、计算耗时之间关系的示意图;
图6为本发明实施例提供的水库下泄水温预测模型的预测结果与实测结果对比示意图。
具体实施方式:
为使本申请的目的、技术方案和优点更加清楚,下面将结合本申请具体实施例及相应的附图对本申请技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。以下结合附图和具体实施例对本发明作具体的介绍。
如图1所示,大型水库叠梁门分层取水设施的下泄水温快速预测方法,包括以下步骤:
S1、基于对水库水温结构形成物理机理的分析,梳理影响水库下泄水温的主要因素;
S2、将步骤S1得到的影响水库下泄水温的主要因素纳入模型的输入因子集X1,以所述下泄水温为模型输出因子Y,初步搭建水库下泄水温预测模型框架,经初步筛选,模型的输入因子集X1包括当前时段所处的月份M、入流量Qi t、出流量QO t、入流水温Ti t、坝上水位Zt、取水口深度Dt、叠梁门运行层数Lt,以及水库当前时段气温Ta t、气压Pt、湿度Ht、太阳辐照度At、风速Wt、云量Ct、蒸发量Et、降水量Rt,预测变量即输出因子Y为水库当前时段的下泄水温To t;
S3、采集所述输入因子集X1和输出因子Y的长时间序列资料,整理形成数据集S1,则所述数据集S1={X1,Y}={M、t、Qi t、Qo t、Ti t、Zt、Dt、Lt、Ta t、Pt、Ht、At、Wt、Ct、Et、Rt、To t},其中t为当前时刻在日内所处时段,t∈[0,24];
S4、基于所述数据集S1,利用敏感性分析及逐步测试,对输入因子进行二次筛选,最终选出模型的输入因子集X2;
S5、对数据集S2做标准化处理;
S6、以经过步骤S5标准化处理后的影响因子数据X2为模型输入向量,以所述标准化处理后的预测变量Y为模型输出向量,构建基于长短期记忆网络(Long and Short-TermMemory network,LSTM)深度学习算法的大型水库叠梁门分层取水下泄水温预测模型,并根据所述标准化处理后的数据集确定训练样本集和测试样本集,其中,所述水库叠梁门分层取水下泄水温预测模型的函数表达式为To t=f(X2),f表示LSTM计算模块;
S7、将所述数据集S2划分为训练样本集和测试样本集;
S8、设定所述水温预测模型的隐层层数、训练批量、隐层节点数和最大迭代次数,结合训练样本集,采用长短期记忆网络算法和Holdout训练方法训练并生成水库下泄水温预测模型;
S9、利用测试样本集对所述水库下泄水温预测模型能进行测试,测试不同隐层层数、训练批量、隐层节点数和最大迭代次数的结构参数组合对所述水库下泄水温预测模型精度、收敛情况和计算耗时的影响,根据测试结果,在保证水库下泄水温预测模型收敛的基础上,重点考虑模型精度,结合模型计算耗时,筛选出最优隐层层数、训练批量、最大迭代次数和隐层节点数;
S10、根据所述最优隐层层数、最优训练批量、最优迭代次数和最优隐层节点数,确定水库下泄水温预测模型,并得到下泄水温预测结果。
进一步地,步骤S6中的LSTM计算模块包括一个输入门,一个记忆细胞,一个遗忘门和一个输出门,计算过程如下:
ft=σ(Wxfxt+Whfht-1+bf)
it=σ(Wxixt+Whiht-1+bi)
Ct=ft·Ct-1+it·tanh(Wxcxt+Whcht-1+bc)
ot=σ(Wxoxt+Whoht-1+bo)
ht=ot·tanh(Ct)
ft、it、Ct、ot和ht分别表示遗忘门、输入门、细胞状态、输出门和隐含层的记忆状态或输出;
Whi表示权重系数矩阵;
bi表示偏置项,
σ为Sigmoid激活函数;
tanh为双曲正切激活函数。
进一步地,步骤S1中影响水库下泄水温的主要因素包括水库调度方式、叠梁门运行方式、入流水温、气象条件。
进一步地,步骤S4包括:以下步骤:
S41、利用敏感性分析,计算各输入因子对输出因子Y的影响程度,将输入因子集分为高影响度输入因子集X1h和低影响度输入因子集X1s;
S42、高影响度输入因子集X1h直接选定为模型输入因子,低影响度输入因子集X1s通过逐步测试,逐一加入到模型输入因子,若加入后对模型精度的提高度5%,则加入模型输入因子,反之则不加入;最终筛选出的模型输入因子集记为X2;,最终整理形成的数据集记为S2,则S2={X2,Y}。
进一步地,步骤S7中将所述数据集S2确定训练样本和测试样本包括以下步骤:
基于二八定律,将数据集划分为训练样本和测试样本,其中,所述训练样本约占所述数据集总年份的80%,所述测试样本约占所述数据集中总年份的20%。
进一步地,步骤S9中测试结构参数组合对模型性能影响测试进一步包含以下步骤:
在统一设定其他参数不同的基础上,依次测试隐层数、训练批量、迭代次数和隐层节点数对模型精度、收敛情况和计算速度的影响,在保证模型收敛的基础上,重点考虑模型精度,结合模型计算耗时,筛选出最优隐层层数、训练批量、最大迭代次数和隐层节点数设定组合。
水库下泄水温预测模型的计算精度的评价指标采用平均绝对误差(MAE),计算公式按式(1)计算得到:
其中,隐层层数对模型精度和计算耗时的影响如附图2a和图2b所示,随着隐层层数的增加,模型精度呈现出先增加后减小的趋势,在隐层层数等于5时,模型平均精度达到最高;模型的训练耗时随着隐层层数的增加呈逐渐上升趋势;
训练批量对模型精度和计算耗时的影响如附图3a和图3b所示,训练批量较小时,模型容易发散,同时随着训练批量值的增大,模型收敛所需的迭代次数增多;结合训练耗时测试结果,随着批量值的增加,模型的训练耗时逐渐降低,综合考虑两项因素,批量值为64,时模型易于收敛,且计算耗时相对较低;
隐层节点数对模型精度和计算耗时的影响如附图4a和图4b所示,随着隐层节点数的增加模型的平均误差呈现出逐渐下降趋势,根据本次测试结果隐层节点数为38时,模型的平均误差相对最低;模型的训练耗时随着隐层节点数的增加呈逐渐上升趋势;
迭代次数对模型精度和计算耗时的影响如附图5a和图5b所示,随着迭代次数的增加,模型评价误差逐渐降低,迭代次数高于300后,模型精度逐渐趋于稳定,根据本次测试结果,迭代次数为500时,模型的预测精度相对最高;模型的计算耗时随迭代次数增加呈逐渐上升趋势;
根据测试结果,最终确定最优隐层层数为5、最优训练批量为64、最大迭代次数为500、最优隐层节点数为38;
步骤S10中,预测结果如图6所示,模型的预测误差0.385℃,训练耗时为5.4min,预测耗时为4.2s。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的***,装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在此提供的算法和显示不与任何特定计算机、虚拟***或者其它设备固有相关。各种通用***也可以与基于在此的示教一起使用。根据上面的描述,构造这类***所要求的结构是显而易见的。此外,本发明也不针对任何特定编程语言。应当明白,可以利用各种编程语言实现在此描述的本发明的内容,并且上面对特定语言所做的描述是为了披露本发明的最佳实施方式。
以上仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (5)
1.大型水库叠梁门分层取水设施的下泄水温快速预测方法,其特征在于,包括以下步骤:
S1、基于对水库水温结构形成物理机理的分析,梳理影响水库下泄水温的主要因素;
S2、将步骤S1得到的影响水库下泄水温的主要因素纳入模型的输入因子集X1,以所述下泄水温为模型输出因子Y,初步搭建水库下泄水温预测模型框架,经初步筛选,模型的输入因子集X1包括当前时段所处的月份M、入流量Qi t、出流量QO t、入流水温Ti t、坝上水位Zt、取水口深度Dt、叠梁门运行层数Lt,以及水库当前时段气温Ta t、气压Pt、湿度Ht、太阳辐照度At、风速Wt、云量Ct、蒸发量Et、降水量Rt,预测变量即输出因子Y为水库当前时段的下泄水温To t;
S3、采集所述输入因子集X1和输出因子Y的长时间序列资料,整理形成数据集S1,则所述数据集S1={X1,Y}={M、t、Qi t、Qo t、Ti t、Zt、Dt、Lt、Ta t、Pt、Ht、At、Wt、Ct、Et、Rt、To t},其中t为当前时刻在日内所处时段,t∈[0,24];
S4、基于所述数据集S1,利用敏感性分析及逐步测试,对输入因子进行二次筛选,最终选出模型的输入因子集X2;
S5、对数据集S2做标准化处理;
S6、以经过步骤S5标准化处理后的影响因子数据X2为模型输入向量,以所述标准化处理后的预测变量Y为模型输出向量,构建基于长短期记忆网络深度学习算法的大型水库叠梁门分层取水下泄水温预测模型,并根据所述标准化处理后的数据集确定训练样本集和测试样本集,其中,所述水库叠梁门分层取水下泄水温预测模型的函数表达式为To t=f(X2),f表示LSTM计算模块;
S7、将所述数据集S2划分为训练样本集和测试样本集;
S8、设定所述水温预测模型的隐层层数、训练批量、隐层节点数和最大迭代次数,结合训练样本集,采用长短期记忆网络算法和Holdout训练方法训练并生成水库下泄水温预测模型;
S9、利用测试样本集对所述水库下泄水温预测模型能进行测试,测试不同隐层层数、训练批量、隐层节点数和最大迭代次数的结构参数组合对所述水库下泄水温预测模型精度、收敛情况和计算耗时的影响,根据测试结果,在保证水库下泄水温预测模型收敛的基础上,重点考虑模型精度,结合模型计算耗时,筛选出最优隐层层数、训练批量、最大迭代次数和隐层节点数;
S10、根据所述最优隐层层数、最优训练批量、最优迭代次数和最优隐层节点数,确定水库下泄水温预测模型,并得到下泄水温预测结果,其中:
步骤S4包括以下步骤:
S41、利用敏感性分析,计算各输入因子对输出因子Y的影响程度,将输入因子集分为高影响度输入因子集X1h和低影响度输入因子集X1s;
S42、高影响度输入因子集X1h直接选定为模型输入因子,低影响度输入因子集X1s通过逐步测试,逐一加入到模型输入因子,若加入后对模型精度的提高度5%,则加入模型输入因子,反之则不加入;最终筛选出的模型输入因子集记为X2,最终整理形成的数据集记为S2,则S2={X2,Y}。
2.如权利要求1所述的大型水库叠梁门分层取水设施的下泄水温快速预测方法,其特征在于,步骤S6中的LSTM计算模块包括一个输入门,一个记忆细胞,一个遗忘门和一个输出门,计算过程如下:
ft=σ(Wxfxt+Whfht-1+bf)
it=σ(Wxixt+Whiht-1+bi)
Ct=ft·Ct-1+it·tanh(Wxcxt+Whcht-1+bc)
ot=σ(Wxoxt+Whoht-1+bo)
ht=ot·tanh(Ct)
ft、it、Ct、ot和ht分别表示遗忘门、输入门、细胞状态、输出门和隐含层的记忆状态或输出;
Whi表示权重系数矩阵;
bi表示偏置项,
σ为Sigmoid激活函数;
tanh为双曲正切激活函数。
3.如权利要求1所述的大型水库叠梁门分层取水设施的下泄水温快速预测方法,其特征在于,步骤S1中影响水库下泄水温的主要因素包括水库调度方式、叠梁门运行方式、入流水温、气象条件。
4.如权利要求1所述的大型水库叠梁门分层取水设施的下泄水温快速预测方法,其特征在于,步骤S7中将所述数据集S2确定训练样本和测试样本包括以下步骤:
基于二八定律,将数据集划分为训练样本和测试样本,其中,所述训练样本占所述数据集总年份的80%,所述测试样本占所述数据集中总年份的20%。
5.如权利要求1所述的大型水库叠梁门分层取水设施的下泄水温快速预测方法,其特征在于,步骤S9中测试结构参数组合对模型性能影响测试进一步包含以下步骤:
在统一设定其他参数不同的基础上,依次测试隐层数、训练批量、迭代次数和隐层节点数对模型精度、收敛情况和计算速度的影响,在保证模型收敛的基础上,重点考虑模型精度,结合模型计算耗时,筛选出最优隐层层数、训练批量、最大迭代次数和隐层节点数设定组合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011040789.5A CN112182709B (zh) | 2020-09-28 | 2020-09-28 | 大型水库叠梁门分层取水设施的下泄水温快速预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011040789.5A CN112182709B (zh) | 2020-09-28 | 2020-09-28 | 大型水库叠梁门分层取水设施的下泄水温快速预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112182709A CN112182709A (zh) | 2021-01-05 |
CN112182709B true CN112182709B (zh) | 2024-01-16 |
Family
ID=73944746
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011040789.5A Active CN112182709B (zh) | 2020-09-28 | 2020-09-28 | 大型水库叠梁门分层取水设施的下泄水温快速预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112182709B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113033618B (zh) * | 2021-03-03 | 2022-03-15 | 四川大学 | 基于支持向量回归的分层水库取水下泄水温预测模型及预测方法 |
CN113283191A (zh) * | 2021-05-26 | 2021-08-20 | 北京华可实工程技术有限公司 | 一种引水渠道水温迭代预测方法 |
CN113704696B (zh) * | 2021-08-31 | 2022-05-27 | 珠江水资源保护科学研究所 | 一种水库水温结构判别方法及判别设备 |
CN116595381B (zh) * | 2023-06-13 | 2024-04-02 | 大连理工大学 | 一种水库分层水温模拟方法及*** |
CN116502569B (zh) * | 2023-06-28 | 2024-03-08 | 长江三峡集团实业发展(北京)有限公司 | 一种基于数值模拟的河道型水库下泄水温调控方法及装置 |
CN117421558B (zh) * | 2023-10-26 | 2024-06-21 | 华中科技大学 | 一种梯级水库运行规则提取及其模型训练方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105046069A (zh) * | 2015-07-06 | 2015-11-11 | 中国水利水电科学研究院 | 一种大型水电站库水温度数值预测方法 |
CN105300558A (zh) * | 2015-11-20 | 2016-02-03 | 中国电建集团贵阳勘测设计研究院有限公司 | 一种水库叠梁门分层取水水温观测装置及其观测方法 |
CN105976101A (zh) * | 2016-04-29 | 2016-09-28 | 武汉大学 | 一种基于svm与dpy的预报-决策耦合水库调度方法 |
CN106407671A (zh) * | 2016-09-08 | 2017-02-15 | 河海大学 | 面向产粘沉性卵鱼类繁殖需求的梯级水库调控***及方法 |
CN109543912A (zh) * | 2018-11-29 | 2019-03-29 | 中国水利水电科学研究院 | 基于深度学习的水库最优调度决策模型生成方法 |
CN110009193A (zh) * | 2019-03-07 | 2019-07-12 | 华北水利水电大学 | 考虑下游农业作物适宜生长的水库生态调度方法 |
CN110059443A (zh) * | 2019-05-08 | 2019-07-26 | 河海大学 | 一种分层水库取水下泄水温的快速预测方法 |
CN110175416A (zh) * | 2019-05-30 | 2019-08-27 | 三峡大学 | 基于主成分分析和神经网络的三峡水库水温预测方法 |
WO2020098256A1 (zh) * | 2018-11-14 | 2020-05-22 | 平安科技(深圳)有限公司 | 基于全卷积神经网络的语音增强方法、装置及存储介质 |
CN111400973A (zh) * | 2020-04-21 | 2020-07-10 | 中国水利水电科学研究院 | 一种基于水文监测数据构建流量-水面宽的关系曲线的方法 |
CN211030535U (zh) * | 2019-08-09 | 2020-07-17 | 中国水利水电科学研究院 | 一种通水调温*** |
CN111563706A (zh) * | 2020-03-05 | 2020-08-21 | 河海大学 | 一种基于lstm网络的多变量物流货运量预测方法 |
CN112633584A (zh) * | 2020-12-29 | 2021-04-09 | 中国地质大学(武汉) | 基于改进的LSTM-seq2seq模型的河流突发水污染事故水质预测方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011159919A2 (en) * | 2010-06-16 | 2011-12-22 | Conocophillips Company | In situ methanogenesis modeling and risk analysis |
US8523425B2 (en) * | 2010-07-01 | 2013-09-03 | Mark Malmquist | Reservoir tank water mixing system |
US10604814B2 (en) * | 2017-09-27 | 2020-03-31 | International Business Machines Coporation | Manufacturing process control with deep learning-based predictive model for hot metal temperature of blast furnace |
-
2020
- 2020-09-28 CN CN202011040789.5A patent/CN112182709B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105046069A (zh) * | 2015-07-06 | 2015-11-11 | 中国水利水电科学研究院 | 一种大型水电站库水温度数值预测方法 |
CN105300558A (zh) * | 2015-11-20 | 2016-02-03 | 中国电建集团贵阳勘测设计研究院有限公司 | 一种水库叠梁门分层取水水温观测装置及其观测方法 |
CN105976101A (zh) * | 2016-04-29 | 2016-09-28 | 武汉大学 | 一种基于svm与dpy的预报-决策耦合水库调度方法 |
CN106407671A (zh) * | 2016-09-08 | 2017-02-15 | 河海大学 | 面向产粘沉性卵鱼类繁殖需求的梯级水库调控***及方法 |
WO2020098256A1 (zh) * | 2018-11-14 | 2020-05-22 | 平安科技(深圳)有限公司 | 基于全卷积神经网络的语音增强方法、装置及存储介质 |
CN109543912A (zh) * | 2018-11-29 | 2019-03-29 | 中国水利水电科学研究院 | 基于深度学习的水库最优调度决策模型生成方法 |
CN110009193A (zh) * | 2019-03-07 | 2019-07-12 | 华北水利水电大学 | 考虑下游农业作物适宜生长的水库生态调度方法 |
CN110059443A (zh) * | 2019-05-08 | 2019-07-26 | 河海大学 | 一种分层水库取水下泄水温的快速预测方法 |
CN110175416A (zh) * | 2019-05-30 | 2019-08-27 | 三峡大学 | 基于主成分分析和神经网络的三峡水库水温预测方法 |
CN211030535U (zh) * | 2019-08-09 | 2020-07-17 | 中国水利水电科学研究院 | 一种通水调温*** |
CN111563706A (zh) * | 2020-03-05 | 2020-08-21 | 河海大学 | 一种基于lstm网络的多变量物流货运量预测方法 |
CN111400973A (zh) * | 2020-04-21 | 2020-07-10 | 中国水利水电科学研究院 | 一种基于水文监测数据构建流量-水面宽的关系曲线的方法 |
CN112633584A (zh) * | 2020-12-29 | 2021-04-09 | 中国地质大学(武汉) | 基于改进的LSTM-seq2seq模型的河流突发水污染事故水质预测方法 |
Non-Patent Citations (5)
Title |
---|
Research on Resource Management System of Multidimensional Reservoir Simulation;Ruihe Wang;2020 IEEE International Conference on Power, Intelligent Computing and Systems (ICPICS);100-102 * |
基于多源数据融合的金沙江下游- 三峡梯级水库调度***设计与应用;杨旭;水利水电技术;154-165 * |
水库水温数值预测方法;胡平;水利学报;1045-1053 * |
水电站叠梁门多层取水下泄水温公式;高学平;陈弘;宋慧芳;;中国工程科学(第12期);65-69 * |
糯扎渡水电站叠梁门试运行期实测水温与数值模拟水温对比分析;李坤;曹晓红;温静雅;简云忠;;水利水电技术(第11期);159-165+189 * |
Also Published As
Publication number | Publication date |
---|---|
CN112182709A (zh) | 2021-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112182709B (zh) | 大型水库叠梁门分层取水设施的下泄水温快速预测方法 | |
CN110298501B (zh) | 基于长短时记忆神经网络的电负荷预测方法 | |
CN108764539B (zh) | 一种梯级电站的上下游水位预测方法 | |
CN111860982A (zh) | 一种基于vmd-fcm-gru的风电场短期风电功率预测方法 | |
CN112116147A (zh) | 一种基于lstm深度学习的河流水温预测方法 | |
CN109146162B (zh) | 一种基于集成循环神经网络的概率风速预测方法 | |
CN113705877B (zh) | 基于深度学习模型的实时月径流预报方法 | |
CN113033618B (zh) | 基于支持向量回归的分层水库取水下泄水温预测模型及预测方法 | |
CN112036658B (zh) | 一种基于随差遗忘长短期记忆的风电功率实时预测方法 | |
CN111144644B (zh) | 基于变分异方差高斯过程回归的短期风速预测方法 | |
CN112434848A (zh) | 基于深度信念网络的非线性加权组合风电功率预测方法 | |
CN113139329B (zh) | 一种基于水文相似性和人工神经网络的新安江模型参数率定方法 | |
CN112288164A (zh) | 一种计及空间相关性和修正数值天气预报的风功率组合预测方法 | |
CN105512832A (zh) | 基于时变权最小方差的城市需水量组合预测方法 | |
CN106778838A (zh) | 一种预测空气质量的方法 | |
CN114139305A (zh) | 基于汽轮机调节级压力预测模型的单阀门流量特性优化方法 | |
CN112801416A (zh) | 基于多维水文信息的lstm流域径流量预测方法 | |
CN113344288A (zh) | 梯级水电站群水位预测方法、装置及计算机可读存储介质 | |
CN114219131A (zh) | 一种基于lstm的流域径流预测方法 | |
CN116757057A (zh) | 基于pso-ga-lstm模型的空气质量预测方法 | |
CN114169639A (zh) | 一种企业用电量预测方法 | |
CN115239029A (zh) | 考虑功率时序及气象相依特性的风电功率预测方法及*** | |
CN110909492A (zh) | 一种基于极端梯度提升算法的污水处理过程软测量方法 | |
CN116611702A (zh) | 一种面向楼宇综合能源管理的集成学习光伏发电预测方法 | |
CN114861555A (zh) | 一种基于Copula理论的区域综合能源***短期负荷预测方法 |
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 |