CN113673775A - 基于cnn-lstm及深度学习的时空组合预测方法 - Google Patents

基于cnn-lstm及深度学习的时空组合预测方法 Download PDF

Info

Publication number
CN113673775A
CN113673775A CN202110981226.4A CN202110981226A CN113673775A CN 113673775 A CN113673775 A CN 113673775A CN 202110981226 A CN202110981226 A CN 202110981226A CN 113673775 A CN113673775 A CN 113673775A
Authority
CN
China
Prior art keywords
lstm
cnn
layer
data
model
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.)
Pending
Application number
CN202110981226.4A
Other languages
English (en)
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.)
Shenzhen Wanzhida Enterprise Management Co ltd
Original Assignee
Wuhan University of Science and Engineering WUSE
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 Wuhan University of Science and Engineering WUSE filed Critical Wuhan University of Science and Engineering WUSE
Priority to CN202110981226.4A priority Critical patent/CN113673775A/zh
Publication of CN113673775A publication Critical patent/CN113673775A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION 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/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • 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/044—Recurrent networks, e.g. Hopfield 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/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/08—Learning methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/004—Generation forecast, e.g. methods or systems for forecasting future energy generation
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/76—Power conversion electric or electronic aspects

Landscapes

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

Abstract

本发明公开了基于CNN‑LSTM及深度学习的时空组合预测方法，包括：S1、通过MI互信息算法对原始数据集的时间和空间相关性进行分析，消除数据间的冗余信息，对原始数据进行降维；S2、建立MI‑CNN‑LSTM模型，对降维后的数据进行预测，CNN用于提取各站点空间信息，LSTM用于获取时间序列数据间的依赖信息；S3、在LSTM网络的LSTM层与输出层之间增加AT层，得到MI‑CNN‑ALSTM模型，并对训练集数据进行模型学习；S4、引入PSO算法对训练参数寻优，得到MI‑CNN‑ALSTM‑PSO模型；S5、对测试集进行预测得到最终的预测值。本发明提出的MI‑CNN‑ALSTM‑PSO时空组合预测模型进一步提高风电功率预测精度和稳定性。

Description

基于CNN-LSTM及深度学习的时空组合预测方法
技术领域
本发明涉及风电预测技术领域，具体涉及一种基于CNN-LSTM及深度学习的时空组合预测方法。
背景技术
在过去的几十年里，全球风力发电装机容量增长迅速。目前，风能作为一种清洁、绿色的可再生能源，被广泛的应用于实际生活中。所以风电场的风电预测为城市建设输送电能与电路建设能提供有效参考依据。
目前风电预测模型可分为三类：时间序列模型，机器学***均模型(Moving Average Model,MA)、自回归移动平均模型(Autoregressive Moving Average Model,ARMA)、自回归整合移动平均模型(Autoregressive Integrated Moving Average Model,ARIMA)等。时间序列模型在线性的、平稳的时间序列预测中具有较好的预测效果。为了改善在非线性和非平稳时间序列的预测效果，一些机器学***稳多步风速预测；Kejun Wang用K-means聚类算法处理数值天气预报，选择对预测精度影响巨大的天气数值作为DBN(深度信念网络)模型的输入，预测结果优于BP神经网络。Yan Jiang提出了基于增强变分模态分解、深度特征选择和多误差修正的概率风速预测框架，结果表明该方法具有较好的预测能力，可能更适合于具有非平稳性和非高斯性的数据；Peng Lu提出了基于时空(ST)分析，多输出支持向量机(MSVM)模型，来预测多个风电场的风电功率；Li Han提出空间重构(PSR)来重构风动力***的相空间，其次使用主成分分析(PCA)最大程度的减少延迟时间和相位尺寸不当所带来的影响，进而使用资源分配网络(RAN)建立预测模型，分析结果表明，在48个小时中的预测误差约为3％，明显低于其他预测方法所获得的误差，并且该模型可以预测不同容量、不同海拔高度的涡轮机风电功率；LiuHongChi提出基于VMD-BPNN组合模型对风电功率预测精度提升。Shengdong Du提出了一种基于双向长短期(Bi-LSTM)的时间注意编码解码模型处理时间序列问题，能自适应学***稳性，提高预测精度。Ping Jiang提出了V-SVM(V支持向量机)混合风速预测模型，选择相邻站点的风速特征作为输入，并采用CS(杜鹃搜索)算法调整模型参数，该模型具有良好的非线性拟合能力和更高的精度。
但是以上方法主要在风速或风电预测过程中，使用信号分解、权重分配、深度学习、回归模型和参数优化等方法进行研究，主要根据时间序列的前后依赖关系进行数据预测，而且大多数风电风速预测方法仅使用单个站点的特征数据，没有考虑风电场相邻站点特征的潜在时空联系，忽略其在空间上的特征信息，使得预测准确性有一定局限性。
发明内容
为了更好利用多站点时空数据预测风电功率，本文利用两个风电场的多站点风电数据的时间和空间的相关性以及深度学习方法进行预测研究，提供了一种基于CNN-LSTM及深度学习的时空组合预测方法，本发明的深度学习算法包括MI算法、AT算法和PSO算法。
本发明克服其技术问题所采用的技术方案是：
一种基于CNN-LSTM及深度学习的时空组合预测方法，包括如下步骤：
S1、通过MI互信息算法对原始数据集的时间和空间相关性进行分析，消除数据间的冗余信息，对原始数据进行降维；
S2、建立MI-CNN-LSTM模型，对步骤S1降维后的数据进行预测，其中，CNN网络用于提取各站点空间信息，LSTM网络用于获取时间序列数据间的依赖信息；
S3、在LSTM网络的LSTM层与输出层之间增加AT层，得到MI-CNN-ALSTM模型，并对训练集数据进行模型学习；
S4、引入PSO算法对训练参数寻优，得到MI-CNN-ALSTM-PSO模型；
S5、通过MI-CNN-ALSTM-PSO模型对测试集进行预测，得到最终的预测值。
进一步地，步骤S1中，原始数据集是采集的过去t小时内n个站点的m维特征的风电数据，表示为：
X＝(X1，X2,...,Xm)n×t＝(X1，X2,...,Xt)∈Rn×m×t
上式中，n为观测站点编号，t为观测时间序列的窗口长度，m为数据特征维度；
则t时刻特征数据描述为：
Figure BDA0003229182060000031
上式中，WP为风电功率，WS为风速，WD为风向，TE为温度，PA为气压，DE为空气密度，i为站点编号。
进一步地，步骤S2中，将n个站点的m维特征数据组成n行m列矩阵，并将t个时刻的矩阵构造为时间分布层；
CNN网络的输入层维度为n×m×t，卷积层特征图的输出
Figure BDA0003229182060000041
为：
Figure BDA0003229182060000042
上式中，Xi+n,j+m是输入矩阵的第n行m列的值，fcov(.)是选择激活函数，wn,m为卷积核n行m列的权重，bn,m是卷积核偏差，k为滑动窗口大小，输入矩阵使用多个卷积核进行卷积操作；
CNN池化层通过滤波器和滑动窗口步长进行降维下采样，再通过展平层将数据展平为一维阵列
Figure BDA0003229182060000043
其中
Figure BDA0003229182060000044
作为时间序列输入到LSTM网络的LSTM层。
进一步地，步骤S2中，LSTM网络内部包括三个门结构和一个用于存储记忆的状态模块；
设Ct为本LSTM单元存储的状态信息，xt为输入层的输入，ht为本单元隐含层的输出，ft为遗忘门，it为输入门，
Figure BDA0003229182060000045
为当前时刻信息，ot为输出门，“×”表示矩阵元素相乘，“+”表示相加运算，σ为sigmoid函数；
遗忘门：用于控制上一单元状态Ct-1被遗忘的程度，其表达式如下：
ft＝σ(Wf*[ht-1,xt]+bf)
输入门：用于控制哪些信息被加入到本单元中，其表达式如下：
it＝σ(Wi*[ht-1,xt]+bi)
单元存储的状态信息：用于根据ft和it将新信息有选择的记录到Ct中，其表达式如下：
Figure BDA0003229182060000046
Figure BDA0003229182060000047
输出门：用于将Ct激活，并控制Ct被过滤的程度，其表达式如下：
ot＝σ(Wo*[ht-1,xt]+bo)
ht＝ot*tanh(Ct)
其中，Wf、Wi、
Figure BDA0003229182060000051
Wo分别为ft、it、
Figure BDA0003229182060000052
ot对应的权重矩阵，bf、bi、
Figure BDA0003229182060000053
bo分别为ft、it、
Figure BDA0003229182060000054
ot对应的偏置项，tanh为双曲正切激活函数，定义如下：
σ(x)＝1/(1+e-x)
tanh(x)＝(ex-e-x)/(ex+e-x)
传统的LSTM网络结构包括输入层、LSTM层和输出层，输出层则依据下式将ht经过一个全连接层得到最终预测值yt：
yt＝σ(Wy*ht+by)
上式中，Wy和by分别为权重矩阵和偏置项。
进一步地，步骤S3中，在LSTM网络的LSTM层与输出层之间增加AT层之后，通过下式得到注意力机制矢量Gt：
Figure BDA0003229182060000055
上式中，
Figure BDA0003229182060000056
为softmax激活函数，Wh为权重向量，bh为偏置项；
依据下式将Gt与ht作元素乘，得到AT层的输出
Figure BDA0003229182060000057
Figure BDA0003229182060000058
根据步骤S2中的yt＝σ(Wy*ht+by)，将
Figure BDA0003229182060000059
作为后续LSTM输出层的输入，并得到预测值yt：
Figure BDA00032291820600000510
进一步地，步骤S4具体包括如下步骤：
S41、初始化：随机初始化粒子参数Pξ，并定义初始种群的粒子，赋予其初始速度和位置；
S42、迭代：每个粒子在探索空间中单独的搜寻最优解，将其记录为当前个体极值，迭代函数如下：
vξ+1＝vξ+c1*rand()*(pbestξ-xξ)+c2*rand()*gbestξ-xξ)
xξ+1＝xξ+vξ+1
上式中，vξ是第ξ次迭代的速度，c1和c2是学习因子，rand()是介于(0,1)之间的随机数，粒子通过追寻局部最优解pbest(Pξ)和全局最优解gbest(Pξ)来确定后续位置，xξ是第ξ次迭代的位置；
S43、寻优：每个粒子单独进行寻优，并将个体极值与整个粒子群里的其他粒子共享，找到最优的那个个体极值作为整个粒子群的当前全局最优解；
S44、评估：每次迭代过程均对组合模型进行训练并计算出预测值及评价指标RMSE；
S45、停止迭代：MI-CNN-ALSTM-PSO模型在搜索过程中达到最大迭代次数且满足RMSE最优则停止，否则返回步骤S42继续更新粒子。
本发明的有益效果是：
1、由于原始数据集为多站点多维特征风电数据，本发明通过MI算法对数据集的时间和空间相关性进行分析，消除数据间的冗余信息，对原始数据集进行降维。
2、本发明考虑各站点的数据相关性，对降维后的数据利用CNN网络提取各站点空间信息，通过LSTM网络获取时间依赖信息，提出CNN-LSTM时空组合预测模型，该模型比常规组合模型有更好的预测性能。
3、本发明设计了CNN-LSTM时空预测模型，并采用MI、AT、PSO等深度学习算法，对其进行改进，提出了MI-CNN-ALSTM-PSO的时空组合预测模型，进一步提高风电功率预测精度和稳定性。
附图说明
图1为本发明实施例的基于CNN-LSTM及深度学习的时空组合预测方法的流程示意图。
图2为本发明实施例的LSTM网络的内部结构示意图。
图3为本发明实施例的CNN-LSTM模型的结构示意图。
图4为本发明实施例的MI-CNN-LSTM模型的结构示意图。
图5为本发明实施例的ALSTM模型的结构示意图。
图6为本发明实施例的MI-CNN-ALSTM模型的结构示意图。
图7为本发明实施例的PSO算法的流程示意图。
图8为本发明实施例的A区站A8风电数据曲线图。
图9为本发明实施例的B区站B8风电数据曲线图。
图10为本发明实施例的A8站点实际风电风速曲线图。
图11为本发明实施例的A区风电和风速数据ACF和PACF图。
图12为本发明实施例的各个模型预测误差对比图，其中，图12(a)为A区预测模型对比图，图12(b)为B区预测模型对比图，图12(a)和(b)中示出的是模型M1、M2、M3、M4、M5，图12(c)和(d)中示出的是模型M6、M7、M8、M9、M10。
图13为本发明实施例的A区预测精度与稳定性对比图。
图14为本发明实施例的B区预测精度与稳定性对比图。
具体实施方式
为了便于本领域人员更好的理解本发明，下面结合附图和具体实施例对本发明做进一步详细说明，下述仅是示例性的不限定本发明的保护范围。
如图1所示，本发明的实施例公开了一种基于CNN-LSTM及深度学习的时空组合预测方法，包括如下步骤：
S1、通过MI互信息算法对原始数据集的时间和空间相关性进行分析，消除数据间的冗余信息，对原始数据进行降维。
S2、建立MI-CNN-LSTM模型，对步骤S1降维后的数据进行预测，其中，CNN网络用于提取各站点空间信息，LSTM网络用于获取时间序列数据间的依赖信息。
S3、在LSTM网络的LSTM层与输出层之间增加AT层，得到MI-CNN-ALSTM模型，并对训练集数据进行模型学习。
S4、引入PSO算法对训练参数寻优，得到MI-CNN-ALSTM-PSO模型。
S5、通过MI-CNN-ALSTM-PSO模型对测试集进行预测，得到最终的预测值。
一、关于MI互信息算法
在风电场预测过程中，各风电站点多维特征变量之间存在一定的相关性。互信息法利用计算互信息函数来确定相关因子，是一种非线性的分析方法，在时间序列预测和特征筛选等领域有着广泛的应用。两个随机变量x和y的互信息可以定义为：
Figure BDA0003229182060000081
其中，p(x,y)是x和y的联合概率密度函数，p(x)和p(y)分别是x和y的边缘概率密度函数。
本实施例中，步骤S1通过MI互信息算法对原始数据集的时间和空间相关性进行分析，消除数据间的冗余信息，对原始数据进行降维。其中，原始数据集是采集的过去t小时内n个站点的m维特征的风电数据，表示为：
X＝(X1，X2,...,Xm)n×t＝(X1，X2,...,Xt)∈Rn×m×t (2)
上式中，n为观测站点编号，t为观测时间序列的窗口长度，m为数据特征维度；
则t时刻特征数据描述为：
Figure BDA0003229182060000082
上式中，WP为风电功率，WS为风速，WD为风向，TE为温度，PA为气压，DE为空气密度，i为站点编号。
本文多维特征风电功率预测问题，是对给定观测序列历史数据X＝(X1，X2,...,Xt-1)，Xt∈Rn×m×t，寻找观测特征变量X和目标预测变量Y(t)之间的非线性映射关系，即找到一个非线性映射函数使得Y(t)＝F(X1，X2,...,Xt-1)。
二、关于CNN卷积神经网络
卷积神经网络(CNN)是一种前馈多层神经网络，从输入到输出只有一个方向，同时具有强大的网格数据处理能力。
CNN卷积步骤主要包括：
卷积层，通过滑动窗口数据与卷积核做内积，提取不同位置数据的权重，得到一个新的矩阵，即特征图。
池化层，是一种降维采样操作，通过类似卷积层的滑动窗口方式得到窗口数据的平均数或最大值，其作用是逐渐降低数据体的空间尺寸，减少网络参数数量，降低计算资源耗费，同时也能够有效控制过拟合。通过这种下采样的方式，可最大限度激活数值大的特征，去除干扰和噪声信息。
全连接层，将池化后提取特征图的所有神经元与本层单元完全连接，再经过输出层的几个隐藏单元得到最后的输出结果。
三、LSTM网络
LSTM可解决循环神经网络(Recurrent Neural Network，RNN)长时间序列训练过程中存在的梯度***、梯度消失与无法长久保存历史信息的问题。相比传统RNN网络，LSTM网络能够在更长的时间序列中有更好的表现。
LSTM网络内部包括三个门结构和一个用于存储记忆的状态模块，如图2所示；
设Ct为本LSTM单元存储的状态信息，xt为输入层的输入，ht为本单元隐含层的输出，ft为遗忘门，it为输入门，
Figure BDA0003229182060000091
为当前时刻信息，ot为输出门，“×”表示矩阵元素相乘，“+”表示相加运算，σ为sigmoid函数。
遗忘门：用于控制上一单元状态Ct-1被遗忘的程度，其表达式如下：
ft＝σ(Wf*[ht-1,xt]+bf) (4)
输入门：用于控制哪些信息被加入到本单元中，其表达式如下：
it＝σ(Wi*[ht-1,xt]+bi) (5)
单元存储的状态信息：用于根据ft和it将新信息有选择的记录到Ct中，其表达式如下：
Figure BDA0003229182060000101
Figure BDA0003229182060000102
输出门：用于将Ct激活，并控制Ct被过滤的程度，其表达式如下：
ot＝σ(Wo*[ht-1,xt]+bo) (8)
ht＝ot*tanh(Ct) (9)
式(4)-(9)中，Wf、Wi、
Figure BDA0003229182060000103
Wo分别为ft、it、
Figure BDA0003229182060000104
ot对应的权重矩阵，bf、bi、
Figure BDA0003229182060000105
bo分别为ft、it、
Figure BDA0003229182060000106
ot对应的偏置项，tanh为双曲正切激活函数，定义如下：
σ(x)＝1/(1+e-x) (10)
tanh(x)＝(ex-e-x)/(ex+e-x) (11)
传统的LSTM网络结构包括输入层、LSTM层和输出层，输出层则依据下式将ht经过一个全连接层得到最终预测值yt：
yt＝σ(Wy*ht+by) (12)
上式中，Wy和by分别为权重矩阵和偏置项。
四、CNN-LSTM模型
CNN具有很好的空间网格数据特征提取能力，LSTM则可以较好的分析时间序列数据间的依赖关系，因此本发明将CNN与LSTM结合。
具体地，步骤S2中，针对本发明的多站点多维特征数据，将n个站点m个特征的数据组成n行m列矩阵，并将t个时刻的矩阵数据构造为时间分布层，以获取更多时间序列数据的步长信息。
CNN模型输入层维度为n×m×t(6×6×30)，CNN卷积层特征图的输出
Figure BDA0003229182060000107
为：
Figure BDA0003229182060000111
上式中，Xi+n,j+m是输入矩阵的第n行m列的值，fcov(.)是选择激活函数，wn,m为卷积核n行m列的权重，bn,m是卷积核偏差，k为滑动窗口大小，输入矩阵使用多个卷积核进行卷积操作。
CNN池化层则通过2×2或3×3滤波器和滑动窗口步长(stride＝2)进行降维下采样，最大限度激活权重大的特征，去除干扰和噪声信息，再通过展平层将数据展平为一维阵列
Figure BDA0003229182060000112
其中
Figure BDA0003229182060000113
作为时间序列输入到LSTM网络的LSTM层。最后，通过完全连接层和输出层获得目标站点的预测结果，如图3所示。
五、关于MI-CNN-LSTM模型
MI序列选择方法，可以实现对称图的特征处理，有效去除约50％的基于模型的特征和95％的对称图像素，且没有明显的损失，这说明MI可以剔除不相关和冗余的信息，识别最重要的特征。本实施例中，多维特征原始数据噪声较多，直接输入CNN-LSTM模型，会对模型预测性能造成干扰，因此基于CNN-LSTM模型，通过MI互信息算法对原始数据集X进行特征选择后再作为CNN的输入，从而减少CNN卷积的数据维度，在保持模型预测精度的前提下减少算法消耗，构造如图4所示的MI-CNN-LSTM模型。
六、MI-CNN-ALSTM模型
传统的LSTM模型，LSTM层直接连接到输出层。本实施例则将LSTM与注意力机制相结合，在LSTM隐含层的最后一个时刻加入注意力机制；现有技术中的LSTM-Attention方法能有效分析输入输出的相关性，以获得更多重要时间的整体特征，通过学习特征权重，提高网络表示能力。而本实施例要解决的是风速剧烈变化导致风电预测结果不稳定的问题，为了让模型更好的作出权重判断，提高模型的预测稳定性，本实施例提出了ALSTM模型，即，在原LSTM层与输出层之间增加AT层(即，Attention层)，其结构如图5所示；利用注意力机制将LSTM层t时刻的输出与MI-CNN-ALSTM模型最终预测结果相关联。AT层使用的激活函数为softmax函数，用于计算隐含层向量的权重；而输出层采用的激活函数为sigmoid函数，用于把向量映射到输出。
具体地，步骤S3中，在LSTM网络的LSTM层与输出层之间增加AT层之后，通过式(14)得到注意力机制矢量Gt：
Figure BDA0003229182060000121
式(14)中，
Figure BDA0003229182060000122
为softmax激活函数，Wh为权重向量，bh为偏置项；
依据下式将Gt与ht作元素乘，得到AT层的输出
Figure BDA0003229182060000123
Figure BDA0003229182060000124
根据步骤S2中的yt＝σ(Wy*ht+by)，将
Figure BDA0003229182060000125
作为后续LSTM输出层的输入，并得到预测值yt：
Figure BDA0003229182060000126
在步骤S2的MI-CNN-LSTM模型的基础上，通过对LSTM网络添加注意力机制，得到如图6所示的MI-CNN-ALSTM模型。
七、PSO参数优化
在使用CNN和LSTM对多维特征进行预测过程中，CNN滤波器个数(n_filters)、CNN卷积核(n_kernel_size)、LSTM隐含层(n_lstm)、LSTM迭代次数(n_batch)等参数对模型预测性能影响较大，需对以上超参数进行优化。而PSO算法是基于鸟群捕食的行为研究的一种进化计算技术，通过群体中个体间的协作和信息共享来寻找最优解。因此，本实施例利用PSO对MI-CNN-ALSTM模型进行的以上参数进行寻优，采用不同参数对测试集进行迭代预测，并通过评价指标RMSE(Root Mean Square Error，均方根差)作为粒子最优解的评估标准，找到最优参数，从而得到最优的风电预测结果。
具体地，步骤S4所述的PSO算法具体步骤如下，如图7所示：
S41、初始化：随机初始化粒子参数Pξ＝(n_filters,n_kernel_size,n_lstm,n_batch)，并定义初始种群的粒子，赋予其初始速度和位置；
S42、迭代：每个粒子在探索空间中单独的搜寻最优解，将其记录为当前个体极值，迭代函数如下：
vξ+1＝vξ+c1*rand()*(pbestξ-xξ)+c2*rand()*gbestξ-xξ) (17)
xξ+1＝xξ+vξ+1 (18)
式(17)和(18)中，vξ是第ξ次迭代的速度，c1和c2是学习因子，rand()是介于(0,1)之间的随机数，粒子通过追寻局部最优解pbest(Pξ)和全局最优解gbest(Pξ)来确定后续位置，xξ是第ξ次迭代的位置；
S43、寻优：每个粒子单独进行寻优，并将个体极值与整个粒子群里的其他粒子共享，找到最优的那个个体极值作为整个粒子群的当前全局最优解，满足公式(17)的计算过程；
S44、评估：每次迭代过程均对组合模型进行训练并计算出预测值及RMSE；
S45、停止迭代：MI-CNN-ALSTM-PSO模型在搜索过程中达到最大迭代次数且满足RMSE最优则停止，否则返回步骤S42继续更新粒子。
八、MI-CNN-ALSTM-PSO模型
通过上述步骤S1-S4，得到MI-CNN-ALSTM-PSO模型，如图1所示，所述MI-CNN-ALSTM-PSO模型主要包括以下部分：
(1)MI互信息：将过去t小时内n个站点的m维特征作为输入数据(30×6×6＝1080)，利用MI对原始数据(1080维)进行时间维度的初步特征选择，得到其中权重最高的前360维数据，组成30×6×2的网格数据。
(2)CNN卷积：对MI特征选择后的数据，进行卷积与池化，充分提取数据空间特征，数据展平成具有时序特征与空间特征的一维阵列，作为LSTM的输入。
(3)CNN-LSTM模型：CNN提取多站点空间信息，LSTM获时间序列记忆信息，CNN-LSTM模型充分考虑风电场各站点的时间和空间信息，进一步提高多站点预测能力。
(4)AT注意力机制：通过结合了局部AT机制的LSTM-AT模型对训练集数据进行模型学习。进一步注意LSTM在最后时刻的输出权重，提高预测能力。
(5)PSO参数选优：采用PSO算法对模型训练参数进行寻优。
(6)数据预测：采用PSO寻优之后的网络模型对测试集进行预测，得到其预测值。
九、实验及分析
本实施例对所述的基于CNN-LSTM及深度学习的时空组合预测方法进行实验及分析。
(一)实验数据集
实验数据来源美国国家可再生能源实验室(NERL)，取自圣克鲁斯岛(约西经120°，北纬34°)风电场和圣米格尔岛(约西经120°，北纬34°)风电场。从2019年7月28日到10月5日(10000条数据)，数据集采样时间间隔为10分钟。本数据集可在美国国家可再生能源实验室(NERL)官网获得。
所选风电场位于太平洋西南岸，地势西北低，主导风向为西北方向，风电站的研究区域如表1所示。
表1站点位置信息
Figure BDA0003229182060000141
Figure BDA0003229182060000151
选取A8和B8站点作为目标站点，其风电数据如图8和图9所示，实验采用前8000条数据(实曲线部分)作为训练集，后2000条数据(虚曲线部分)作为测试集；这里要说明的是，图8和9中的水平实线均代表平均值。该站点风电数据有如下特点：
1)风电数据具有高随机性与波动性；
2)该区域地理位置比较特殊，位于海上岛屿，风速变化强烈，风向难以捉摸，大气环境也变化无常；这些外部条件会导致风电预测更加困难，也可以进一步验证本实发明提出的组合预测模型的有效性。
(二)数据处理与分析
1)归一化
将有量纲的原始数据集变换为0～1范围之内无量纲的标量，这样不仅可以使后续模型处理更加快速，而且便于把不同单位或量级的指标能够进行比较和加权。
y＝(x-Min)/(Max-Min) (19)
其中，x为原始数据，Min为数据集x中的最大值，Min为数据集x中的最小值，y为归一化后的输出值。
2)数据相关性
图10为A8站点的风速和风电功率的趋势对比，从变化趋势可知，风电和风速数据具有一致的同步变化趋势。另外，A区各站点实际风电数据具有高度同步性。
表2 A区各个站点风电与预测点风电的相关性
Figure BDA0003229182060000161
表2为A区站点与目标站点距离和目标站点WP的相关性分析。从中可看出，临近A8站点的特征，与其相关性较大。反之，则相关性较小，风电数据相关性与距离成正相关。通过A区各站点的风电数据相关性分析可知，A5、A6、A8、A10、A11、A12站点距离和相关性具有良好的参考性，因此用这些站点的相关数据进行预测。
B区站点相关性分析如表3所示，同理，B区选用B8作为待预测中心站点，而选用B6、B7、B8、B9、B10、B13站点的相关数据进行预测。
表3 B区各个站点风电与预测点风电的相关性
Figure BDA0003229182060000171
3)风电自相关性分析
为了确定后续预测模型的时间步长，采用ACF和PACF准确分析风电数据的自相关性。从图11可以看出，在风电数据区间内，ACF图特征是拖尾，而PACF图特征是截尾，因此该风电数据满足自回归(AR)模型的特性，而从ACF图中可以看出风电数据在滞后30左右的时候进入置信区间之内，因此预测模型的时间步初步确定为30。
4)特征维度
表4为预测站点多维特征相关性分析。从中可以看出风电和风速具有高度相关性，风向也是影响风电一个重要因素，而密度、气压等，影响风电的程度微乎其微。
表4预测站点多维特征相关性对比
Figure BDA0003229182060000172
(三)评价指标
在本次实验中使用如下三种评估指标。其中：N为风电预测样本总个数，y(t)为t时刻原始观测站记录风电数据，
Figure BDA0003229182060000181
为t时刻风电预测值。MAE、RMSE和MASE值越小，表示模型预测精度越高。
平均绝对误差：
Figure BDA0003229182060000182
平均方根误差：
Figure BDA0003229182060000183
平均绝对标度误差：
Figure BDA0003229182060000184
(四)对比实验基准模型
表5实验对比基准模型
Figure BDA0003229182060000185
各个模型主要参数设置如表6所示。
表6各个模型主要参数
Figure BDA0003229182060000186
Figure BDA0003229182060000191
(五)预测模型实验数据分析
本节主要对组合预测模型的MI互信息、CNN卷积、LSTM-AT、PSO粒子群优化算法几个部分算法和预测模型，进行实验验证和对比分析。
1)CNN-LSTM卷积模型
各种模型预测指标对比如表9所示，CNN-LSTM的10组平均误差指标MAE、RMSE和MASE分别为0.196，0.414，0.356，比单一的LSTM、VMD-LSTM、EEMD-LSTM、TCN模型误差有明显降低，说明CNN-LSTM有更高的预测精度。
2)MI-CNN-LSTM模型
表7为MI特征选择并采用MI-CNN-LSTM模型预测的结果对比。当MI选择特征维度为900时，其预测指标RMSE为0.398，此时其中WP、WS、WD、TE为180维、DE为19维、PA为161维，相比CNN-LSTM(1080维)模型预测误差降低4％。而选择360维数据时，WP、WS均为180，预测结果最好(RMSE＝0.265)。
表7 MI特征维度选择分析
Figure BDA0003229182060000192
Figure BDA0003229182060000201
说明其他特征维度输入到CNN-LSTM模型中存在数据冗余，影响预测结果。这也与前面相关性分析得到的结果相似。
表8给出了A区MI特征选择360维数据的MI权重前60维数据，其中前20维特征都是A5、A8、A10、A11站点的风电和风速，这也说明MI特征选择的结果与表10各站点距离和特征相关性分析得到的结果相似。
表8 MI特征选择结果
Figure BDA0003229182060000202
从表9的预测结果误差指标可以看出，MI-CNN-LSTM的最终误差指标MAE、RMSE和MASE分别为0.060，0.265，0.164，其误差指标比CNN-LSTM模型降低约70％，36％和50％，说明了MI-CNN-LSTM的预测效果更好。
3)MI-CNN-ALSTM模型
MI-CNN-LSTM-AT组合模型对网格数据进行卷积，并对LSTM结构进行改进，增加局部AT。从表9可以看出，MI-CNN-ALSTM的误差指标MAE、RMSE、MASE分别为0.058、0.247、0152，比MI-CNN-LSTM模型的误差指标分别低8％、10％、7％，说明加入AT后的MI-CNN-ALSTM模型的预测精度明显优于MI-CNN-LSTM，从图12的预测误差对比柱状图中也可以看出，MI-CNN-ALSTM所有评价指标在两个区域的预测误差均为最优。
表9 A区站、B区站预测模型误差对比(MAE(kw),RMSE(kw),MASE(kw))
Figure BDA0003229182060000211
4)MI-CNN-ALSTM-PSO参数寻优
表10为PSO参数寻优结果，当参数n_filters、n_kernel_size、n_lstm、n_batch分别为46、2、37、64时，模型预测精度最高，误差最小，故选择该组参数作为模型训练参数。
表10 MI-CNN-LSTM-AT模型PSO粒子群参数寻优
Figure BDA0003229182060000212
Figure BDA0003229182060000221
(六)组合模型预测精度与稳定性分析
利用模型M3至M10对A区和B区站点风速数据分别进行20组预测实验，并对实验结果的误差指标RMSE进行对比，图13和图14为20组预测结果的箱线图对比。
从图中可以看出：1)CNN-LSTM的预测误差RMSE平均数为0.48，优于EEMD-LSTM、VMD-LSTM和TCN，可以说明CNN卷积在处理多站点相关性的时候明显优于信号分解。2)增加AT和MI的MI-CNN-ALSTM模型平均误差RMSE为0.24，预测效果也有明显提升。3)从稳定性的角度来看，本实施例所提出的模型经过PSO粒子群算法优化之后，其标准差为0.021，比MI-CNN-ALSTM标准差0.025稳定，预测误差进一步降低，同时预测模型的误差范围也更加平稳，更容易训练出更优的模型。
(七)实验总结
综合上述A区和B区的对比实验数据，可以得到如下结论：
1)MI作用：经过MI特征提取，消除冗余数据，如DE、PA特征几乎可以忽略不计。对提取后数据进行预测，实验效果提升明显。
2)常规预测模型：ARIMA模型能预测风电的主要趋势，也即是在风速发生突变与风电出现峰值时不能准确预测；而SVR与LSTM单个模型虽然能预测风电功率的大致走向，但是在预测时存在预测滞后的问题，并不能准确预测下一时刻的风电。TCN、VMD-LSTM和EEMD-LSTM等模型在多站点建立相关模型预测，也难以捕捉其空间相关性，并且算法时间开销较大，预测效果也不好。
3)CNN-LSTM模型：通过将各站点组成的网格卷积进行卷积，CNN卷积可提取多站点的空间信息，建立多站点之间的联系，同时与LSTM结合，又能保证时间特征上的完整性。CNN-LSTM充分结合两个模型的算法优势，提高风电预测性能。
4)MI-CNN-ALSTM模型：从表3可以看出，模型经过MI特征选择后，以CNN-LSTM为基础，在其输出单元最后时刻增加注意力机制，使输出时更加关注上一时刻的权重特征，从而进一步提高预测能力。
5)MI-CNN-ALSTM-PSO模型：主要作用包括：提升模型在风电急剧变化时的适应能力；能够更加准确的预测风电发生突变出现峰值与低谷的情况；基于PSO优化的MI-CNN-ALSTM模型使得预测精度更加稳定。
以上仅描述了本发明的基本原理和优选实施方式，本领域人员可以根据上述描述做出许多变化和改进，这些变化和改进应该属于本发明的保护范围。

Claims (6)

1.一种基于CNN-LSTM及深度学习的时空组合预测方法，其特征在于，包括如下步骤：
S1、通过MI互信息算法对原始数据集的时间和空间相关性进行分析，消除数据间的冗余信息，对原始数据进行降维；
S2、建立MI-CNN-LSTM模型，对步骤S1降维后的数据进行预测，其中，CNN网络用于提取各站点空间信息，LSTM网络用于获取时间序列数据间的依赖信息；
S3、在LSTM网络的LSTM层与输出层之间增加AT层，得到MI-CNN-ALSTM模型，并对训练集数据进行模型学习；
S4、引入PSO算法对训练参数寻优，得到MI-CNN-ALSTM-PSO模型；
S5、通过MI-CNN-ALSTM-PSO模型对测试集进行预测，得到最终的预测值。
2.根据权利要求1所述的方法，其特征在于，步骤S1中，原始数据集是采集的过去t小时内n个站点的m维特征的风电数据，表示为：
X＝(X1，X2,...,Xm)n×t＝(X1，X2,...,Xt)∈Rn×m×t
上式中，n为观测站点编号，t为观测时间序列的窗口长度，m为数据特征维度；
则t时刻特征数据描述为：
Figure FDA0003229182050000011
上式中，WP为风电功率，WS为风速，WD为风向，TE为温度，PA为气压，DE为空气密度，i为站点编号。
3.根据权利要求2所述的方法，其特征在于，步骤S2中，将n个站点的m维特征数据组成n行m列矩阵，并将t个时刻的矩阵构造为时间分布层；
CNN网络的输入层维度为n×m×t，卷积层特征图的输出
Figure FDA0003229182050000012
为：
Figure FDA0003229182050000013
上式中，Xi+n,j+m是输入矩阵的第n行m列的值，fcov(.)是选择激活函数，wn,m为卷积核n行m列的权重，bn,m是卷积核偏差，k为滑动窗口大小，输入矩阵使用多个卷积核进行卷积操作；
CNN池化层通过滤波器和滑动窗口步长进行降维下采样，再通过展平层将数据展平为一维阵列
Figure FDA0003229182050000021
其中
Figure FDA0003229182050000022
作为时间序列输入到LSTM网络的LSTM层。
4.根据权利要求3所述的方法，其特征在于，步骤S2中，LSTM网络内部包括三个门结构和一个用于存储记忆的状态模块；
设Ct为本LSTM单元存储的状态信息，xt为输入层的输入，ht为本单元隐含层的输出，ft为遗忘门，it为输入门，
Figure FDA0003229182050000023
为当前时刻信息，ot为输出门，“×”表示矩阵元素相乘，“+”表示相加运算，σ为sigmoid函数；
遗忘门：用于控制上一单元状态Ct-1被遗忘的程度，其表达式如下：
ft＝σ(Wf*[ht-1,xt]+bf)
输入门：用于控制哪些信息被加入到本单元中，其表达式如下：
it＝σ(Wi*[ht-1,xt]+bi)
单元存储的状态信息：用于根据ft和it将新信息有选择的记录到Ct中，其表达式如下：
Figure FDA0003229182050000024
Figure FDA0003229182050000025
输出门：用于将Ct激活，并控制Ct被过滤的程度，其表达式如下：
ot＝σ(Wo*[ht-1,xt]+bo)
ht＝ot*tanh(Ct)
其中，Wf、Wi、
Figure FDA0003229182050000026
Wo分别为ft、it、
Figure FDA0003229182050000027
ot对应的权重矩阵，bf、bi、
Figure FDA0003229182050000028
bo分别为ft、it、
Figure FDA0003229182050000029
ot对应的偏置项，tanh为双曲正切激活函数，定义如下：
σ(x)＝1/(1+e-x)
tanh(x)＝(ex-e-x)/(ex+e-x)
传统的LSTM网络结构包括输入层、LSTM层和输出层，输出层则依据下式将ht经过一个全连接层得到最终预测值yt：
yt＝σ(Wy*ht+by)
上式中，Wy和by分别为权重矩阵和偏置项。
5.根据权利要求4所述的方法，其特征在于，步骤S3中，在LSTM网络的LSTM层与输出层之间增加AT层之后，通过下式得到注意力机制矢量Gt：
Figure FDA0003229182050000031
上式中，
Figure FDA0003229182050000032
为softmax激活函数，Wh为权重向量，bh为偏置项；
依据下式将Gt与ht作元素乘，得到AT层的输出
Figure FDA0003229182050000033
Figure FDA0003229182050000034
根据步骤S2中的yt＝σ(Wy*ht+by)，将
Figure FDA0003229182050000035
作为后续LSTM输出层的输入，并得到预测值yt：
Figure FDA0003229182050000036
6.根据权利要求1所述的方法，其特征在于，步骤S4具体包括如下步骤：
S41、初始化：随机初始化粒子参数Pξ，并定义初始种群的粒子，赋予其初始速度和位置；
S42、迭代：每个粒子在探索空间中单独的搜寻最优解，将其记录为当前个体极值，迭代函数如下：
vξ+1＝vξ+c1*rand()*(pbestξ-xξ)+c2*rand()*gbestξ-xξ)
xξ+1＝xξ+vξ+1
上式中，vξ是第ξ次迭代的速度，c1和c2是学习因子，rand()是介于(0,1)之间的随机数，粒子通过追寻局部最优解pbest(Pξ)和全局最优解gbest(Pξ)来确定后续位置，xξ是第ξ次迭代的位置；
S43、寻优：每个粒子单独进行寻优，并将个体极值与整个粒子群里的其他粒子共享，找到最优的那个个体极值作为整个粒子群的当前全局最优解；
S44、评估：每次迭代过程均对组合模型进行训练并计算出预测值及评价指标RMSE；
S45、停止迭代：MI-CNN-ALSTM-PSO模型在搜索过程中达到最大迭代次数且满足RMSE最优则停止，否则返回步骤S42继续更新粒子。
CN202110981226.4A 2021-08-25 2021-08-25 基于cnn-lstm及深度学习的时空组合预测方法 Pending CN113673775A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110981226.4A CN113673775A (zh) 2021-08-25 2021-08-25 基于cnn-lstm及深度学习的时空组合预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110981226.4A CN113673775A (zh) 2021-08-25 2021-08-25 基于cnn-lstm及深度学习的时空组合预测方法

Publications (1)

Publication Number Publication Date
CN113673775A true CN113673775A (zh) 2021-11-19

Family

ID=78546098

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110981226.4A Pending CN113673775A (zh) 2021-08-25 2021-08-25 基于cnn-lstm及深度学习的时空组合预测方法

Country Status (1)

Country Link
CN (1) CN113673775A (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114913296A (zh) * 2022-05-07 2022-08-16 中国石油大学(华东) 一种modis地表温度数据产品重建方法
CN114936669A (zh) * 2022-04-06 2022-08-23 武汉大学 一种基于数据融合的混合船舶横摇预测方法
CN114997543A (zh) * 2022-08-03 2022-09-02 通号通信信息集团有限公司 人流量预测方法、模型训练方法、电子设备、可读介质
CN115099497A (zh) * 2022-06-28 2022-09-23 中国水利水电科学研究院 一种基于cnn-lstm的实时洪水预报智能方法
CN115587617A (zh) * 2022-10-17 2023-01-10 国网黑龙江省电力有限公司 长时间序列的序列信息捕捉方法及装置、短期用电量预测方法及装置
CN116087042A (zh) * 2023-02-16 2023-05-09 中国矿业大学 一种粉尘沉积量空间分布实时监测***及方法
CN116110234A (zh) * 2023-04-11 2023-05-12 城云科技（中国）有限公司 基于人工智能的车流量预测方法、装置及其应用
CN116861467A (zh) * 2023-08-31 2023-10-10 中国科学院软件研究所 一种基于上下文特征的数据库异常查询访问控制方法
CN117034784A (zh) * 2023-09-21 2023-11-10 河北工程大学 一种非稳态轧制下热轧变形区动态力学参数的建模方法
CN117674098A (zh) * 2023-11-29 2024-03-08 国网浙江省电力有限公司丽水供电公司 面向不同渗透率的多元负荷时空概率分布预测方法及***
CN117851922A (zh) * 2024-03-08 2024-04-09 山东新泽仪器有限公司 一种空气浓度监测仪器运行监测方法
CN118036789A (zh) * 2023-06-21 2024-05-14 中国科学院广州能源研究所 一种混合电动汽车负荷分类预测方法及相关设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110443725A (zh) * 2019-08-14 2019-11-12 国网甘肃省电力公司信息通信公司 基于cnn+lstm的风电功率预测模型
CN110991700A (zh) * 2019-11-08 2020-04-10 北京博望华科科技有限公司 基于深度学习改进的天气与用电相关性预测方法及装置
CN112598180A (zh) * 2020-12-21 2021-04-02 北京华能新锐控制技术有限公司 一种分布式区域风电功率预测方法
CN112686464A (zh) * 2021-01-07 2021-04-20 云南电力技术有限责任公司 短期风电功率预测方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110443725A (zh) * 2019-08-14 2019-11-12 国网甘肃省电力公司信息通信公司 基于cnn+lstm的风电功率预测模型
CN110991700A (zh) * 2019-11-08 2020-04-10 北京博望华科科技有限公司 基于深度学习改进的天气与用电相关性预测方法及装置
CN112598180A (zh) * 2020-12-21 2021-04-02 北京华能新锐控制技术有限公司 一种分布式区域风电功率预测方法
CN112686464A (zh) * 2021-01-07 2021-04-20 云南电力技术有限责任公司 短期风电功率预测方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
周恒俊 ET AL.: "基于MIPCA与GRU网络的光伏出力短期预测方法", 电力***及其自动化学报, vol. 32, no. 9, pages 55 - 62 *

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114936669A (zh) * 2022-04-06 2022-08-23 武汉大学 一种基于数据融合的混合船舶横摇预测方法
CN114936669B (zh) * 2022-04-06 2024-05-17 武汉大学 一种基于数据融合的混合船舶横摇预测方法
CN114913296B (zh) * 2022-05-07 2023-08-11 中国石油大学(华东) 一种modis地表温度数据产品重建方法
CN114913296A (zh) * 2022-05-07 2022-08-16 中国石油大学(华东) 一种modis地表温度数据产品重建方法
CN115099497A (zh) * 2022-06-28 2022-09-23 中国水利水电科学研究院 一种基于cnn-lstm的实时洪水预报智能方法
CN115099497B (zh) * 2022-06-28 2023-02-17 中国水利水电科学研究院 一种基于cnn-lstm的实时洪水预报智能方法
CN114997543A (zh) * 2022-08-03 2022-09-02 通号通信信息集团有限公司 人流量预测方法、模型训练方法、电子设备、可读介质
CN114997543B (zh) * 2022-08-03 2023-01-13 通号通信信息集团有限公司 人流量预测方法、电子设备、可读介质
CN115587617A (zh) * 2022-10-17 2023-01-10 国网黑龙江省电力有限公司 长时间序列的序列信息捕捉方法及装置、短期用电量预测方法及装置
CN116087042B (zh) * 2023-02-16 2023-09-08 中国矿业大学 一种粉尘沉积量空间分布实时监测***及方法
CN116087042A (zh) * 2023-02-16 2023-05-09 中国矿业大学 一种粉尘沉积量空间分布实时监测***及方法
CN116110234A (zh) * 2023-04-11 2023-05-12 城云科技（中国）有限公司 基于人工智能的车流量预测方法、装置及其应用
CN118036789A (zh) * 2023-06-21 2024-05-14 中国科学院广州能源研究所 一种混合电动汽车负荷分类预测方法及相关设备
CN116861467A (zh) * 2023-08-31 2023-10-10 中国科学院软件研究所 一种基于上下文特征的数据库异常查询访问控制方法
CN116861467B (zh) * 2023-08-31 2023-12-26 中国科学院软件研究所 一种基于上下文特征的数据库异常查询访问控制方法
CN117034784A (zh) * 2023-09-21 2023-11-10 河北工程大学 一种非稳态轧制下热轧变形区动态力学参数的建模方法
CN117674098A (zh) * 2023-11-29 2024-03-08 国网浙江省电力有限公司丽水供电公司 面向不同渗透率的多元负荷时空概率分布预测方法及***
CN117674098B (zh) * 2023-11-29 2024-06-07 国网浙江省电力有限公司丽水供电公司 面向不同渗透率的多元负荷时空概率分布预测方法及***
CN117851922A (zh) * 2024-03-08 2024-04-09 山东新泽仪器有限公司 一种空气浓度监测仪器运行监测方法
CN117851922B (zh) * 2024-03-08 2024-06-07 山东新泽仪器有限公司 一种空气浓度监测仪器运行监测方法

Similar Documents

Publication Publication Date Title
CN113673775A (zh) 基于cnn-lstm及深度学习的时空组合预测方法
Ma et al. A hybrid attention-based deep learning approach for wind power prediction
Huang et al. Hybrid deep neural model for hourly solar irradiance forecasting
Lipu et al. Artificial intelligence based hybrid forecasting approaches for wind power generation: Progress, challenges and prospects
Zhang et al. Multi-head attention-based probabilistic CNN-BiLSTM for day-ahead wind speed forecasting
Zhang et al. Short-term power load forecasting using integrated methods based on long short-term memory
Liang et al. Clusterfomer: clustering as a universal visual learner
Li et al. A novel hybrid model for multi-step ahead photovoltaic power prediction based on conditional time series generative adversarial networks
CN117117859B (zh) 基于神经网络的光伏发电功率预测方法及***
CN114693064B (zh) 一种建筑群方案生成性能评价方法
Chen et al. Solar energy forecasting with numerical weather predictions on a grid and convolutional networks
Shang et al. Decomposition-based wind speed forecasting model using causal convolutional network and attention mechanism
CN114792158A (zh) 基于时空融合图神经网络的多风电场短期功率预测方法
CN111222689A (zh) 基于多尺度时间特征的lstm负荷预测方法、介质及电子装置
CN115099461A (zh) 基于双分支特征提取的太阳辐射预测方法及***
Xu et al. AM-ConvGRU: A spatio-temporal model for typhoon path prediction
Phiri et al. Deep learning techniques for solar tracking systems: A systematic literature review, research challenges, and open research directions
Xuan et al. A comprehensive evaluation of statistical, machine learning and deep learning models for time series prediction
Dong et al. Tropical cyclone track prediction with an encoding-to-forecasting deep learning model
CN113837434A (zh) 太阳能光伏发电预测方法、装置、电子设备和存储介质
CN113935458A (zh) 基于卷积自编码深度学习的空气污染多站点联合预测方法
CN117194926A (zh) 陆上风电基地吊装窗口期预测方法及***
CN112597890A (zh) 基于多维泰勒网的人脸识别方法
CN116894524A (zh) 一种基于组合网络模型的大气温度预测方法
Xu et al. The hidden-layers topology analysis of deep learning models in survey for forecasting and generation of the wind power and photovoltaic energy

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20240613

Address after: 1003, Building A, Zhiyun Industrial Park, No. 13 Huaxing Road, Tongsheng Community, Dalang Street, Longhua District, Shenzhen City, Guangdong Province, 518000

Applicant after: Shenzhen Wanzhida Enterprise Management Co.,Ltd.

Country or region after: China

Address before: 430081 No. 947 Heping Avenue, Qingshan District, Hubei, Wuhan

Applicant before: WUHAN University OF SCIENCE AND TECHNOLOGY

Country or region before: China