CN115330085A - 基于深度神经网络且无未来信息泄露的风速预测方法 - Google Patents

基于深度神经网络且无未来信息泄露的风速预测方法 Download PDF

Info

Publication number
CN115330085A
CN115330085A CN202211136465.0A CN202211136465A CN115330085A CN 115330085 A CN115330085 A CN 115330085A CN 202211136465 A CN202211136465 A CN 202211136465A CN 115330085 A CN115330085 A CN 115330085A
Authority
CN
China
Prior art keywords
wind speed
prediction
data
sequence data
prediction 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
CN202211136465.0A
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.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN202211136465.0A priority Critical patent/CN115330085A/zh
Publication of CN115330085A publication Critical patent/CN115330085A/zh
Pending legal-status Critical Current

Links

Images

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
    • 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/06Energy or water supply

Landscapes

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

Abstract

本发明公开了一种基于深度神经网络且无未来信息泄露的风速预测方法,包括如下步骤:步骤1:分别风速序列数据进行数据筛选和数据预处理;根据子模态能量占比大小筛选有效成分Tc并作为预测模型的输出值;采用实时滚动分解策略对风速序列数据进行预处理以得到无信息泄露的预测模型的输入值;以输入值和输出值构建数据集;步骤2:确定预测步长:按照实际的预测需求,确定需要预测的未来时间步数;步骤3:构建预测模型:构建结合注意力机制的双向长短期记忆网络作为预测模型,利用数据集对预测模型进行训练和测试;步骤4:风速预测:利用训练得到的预测模型预测风速,得到风速预测结果。本发明方法规避了数据泄露的风险,提高了实用性。

Description

基于深度神经网络且无未来信息泄露的风速预测方法
技术领域
本发明属于风速预测技术领域,具体的为一种基于深度神经网络且无未来信息泄露的风速预测方法。
背景技术
准确的风速预测在天气预报精度的提高、气象灾害的预警以及清洁风电的发展等方面至关重要。在风电领域,准确的短期风速预测可以让运营商对风电场的整体运行进行合理调控,维持稳定的供能并降低对电网的冲击。随着国内陆上风电场向低风速区逐步发展,越来越多的风电场选址在山地区域。由于地形变化等因素的影响,山地区域容易产生尾流效应并出现流动分离再附、剪切流等复杂的流动现象。因此,为了保障山区风电场的安全平稳运行,亟需发展准确预测短期风速的方法。
目前常用的预测方法分为如下几类:物理方法、统计方法、人工智能方法以及混合方法。物理方法需消耗大量的计算资源与时间,在中、长期预测中表现良好,但却难以满足短期预测的时效性要求。统计方法由于其固有的线性假设,导致难以对非线性风速序列进行合理的建模及预测。人工智能方法仅使用了浅层神经网络,缺乏对风速序列深层次非线性特征的提取能力。
随着计算机能力的飞速发展,采用深度神经网络进行风速预测也变得切实可行。虽然基于深度学习的短期风速预测模型取得了良好的预测结果,但不可否认的是,单一的神经网络模型可能难以适用于不同的风速序列。并且,实测风速序列中存在着大量的噪音信息,这严重影响了风速预测的精度。因此,结合数据预处理方法的混合预测模型得到了广泛的发展,其中模态分解方法应用地最为广泛。数据预处理操作确实能显著降低风速的随机性与高频噪声,提高风速的可预测性以及预测模型的性能,极大地提高了风速预测的精度。但是这类混合预测模型在数据预处理时对全部的风速序列数据进行了分解,在分解后的子模态上划分训练集和测试集并建立预测模型。这意味着本应未知的测试集数据被视为已知,不可避免的造成了未来信息的泄露。而当有新数据被添加进原始风速序列数据中时,需要重新进行分解,然后再代入子模态预测模型中进行预测。但即使只在原始数据的末端增加部分新数据,分解后的子模态在序列的末端也会产生明显的变化。这部分末端数据是时序预测中最为关键的信息。这意味着在原始数据的子模态上建立的预测模型可能难以应用于新的数据,导致此类混合预测模型缺少实用性。
发明内容
有鉴于此,本发明的目的在于提供一种基于深度神经网络且无未来信息泄露的风速预测方法,采用事后评估的方式,此方法规避了数据泄露的风险,提高了实用性。
为达到上述目的,本发明提供如下技术方案:
一种基于深度神经网络且无未来信息泄露的风速预测方法,包括如下步骤:
步骤1:分别风速序列数据进行数据筛选和数据预处理;根据子模态能量占比大小筛选有效成分Tc并作为预测模型的输出值;采用实时滚动分解策略对风速序列数据进行预处理以得到无信息泄露的预测模型的输入值;以输入值和输出值构建数据集;
步骤2:确定预测步长:按照实际的预测需求,确定需要预测的未来时间步数;
步骤3:构建预测模型:构建结合注意力机制的双向长短期记忆网络作为预测模型,利用数据集对预测模型进行训练和测试;
步骤4:风速预测:利用训练得到的预测模型预测风速,得到风速预测结果。
进一步,所述步骤一中,从风速序列数据中筛选出有效成分Tc的方法为:
S1:将风速序列数据划分训练集、验证集和测试集,以划分后的数据组建三个数据组,第一个数据组包括训练集,第二个数据组包括验证集,第三个数据组包括测试集;
S2:对三个数据组中的风速序列数据分别进行模态分解;
S3:基于能量原理判断指标,选择合适的子模块得到能代表风速序列数据绝大部分能量的有效分成Tci,i=1,2或3,分别表示从三个数据组的风速序列数据中筛选得到的有效分成;
S4:按顺序将有效成分Tci进行拼接,得到预测模型的目标输出值Tc。
进一步,所述步骤S3中,当前m个子模态能量之和大于等于原始风速序列99%的能量时,Tci为前m个子模态之和,子模态的ER计算方式为:
Figure BDA0003852310430000021
其中,ER(k)表示第k个子模态的能量占比;
Figure BDA0003852310430000022
表示第k个子模态序列;X表示风速序列数据;N表示风速序列数据X的长度。
进一步,所述步骤一中,采用实时滚动分解策略对风速序列数据进行预处理的方法为:
(1)确定窗口长度L:分析风速序列数据X的频谱特征,寻找频率峰值点,得到对应的周期,选择周期数作为RTRD所需的窗口长度;
(2)确定滑动步幅s;
(3)重构风速信号:根据窗口长度L和滑动步幅s将风速序列数据X重构为二维矩阵XR,矩阵形状为(N-L+s,L),其中,N为风速序列数据X的长度;
(4)执行RTRD:逐行分解XR,得到无信息泄露的分解结果
Figure BDA0003852310430000031
进一步,所述预测模型包括依次设置的一层输入层、两层Bi-LSTM网络、一层注意力层和一层全连接层。
进一步,所述LSTM网络的原理为:
ft=σ(Wxfxt+Whfht-1+bf)
gt=tanh(Wxgxt+Whght-1+bg)
it=σ(Wxixt+Whiht-1+bi)
ct=ft⊙ct-1+it⊙gt
ot=σ(Wxoxt+Whoht-1+bo)
yt=ht=ot⊙tanh(ct)
其中,xt为输入向量,ct为长期状态,ht为短期状态,yt为输出向量;Wxf,Wxg,Wxi,Wxo分别为与xt连接的权重矩阵,Whf,Whg,Whi,Who分别为与ht-1连接的权重矩阵,bf,bg,bi,bo分别为四层偏置项;ft为遗忘门的控制器,由ft决定ct-1的哪些部分应被删除;gt为LSTM网络的中间输出,其作用为分析xt与ht-1;it为输入门的控制器,判断gt的哪些重要部分应被添加进ct;ot为输出门的控制器,判断应读取ct中的哪部分并应用到ht和yt中;σ表示激活函数。
进一步,所述注意力层采用Luong-Attention模型,其计算方法为:
score计算方法为:
Figure BDA0003852310430000032
注意力权重αts计算方法为:
Figure BDA0003852310430000033
其中,
Figure BDA0003852310430000034
和ht分别为编码器输出的全部隐藏状态以及最后一个隐藏状态;W表示权重矩阵;
上下文向量
Figure BDA0003852310430000035
计算方法为:
Figure BDA0003852310430000036
其中,αts是注意力权重,
Figure BDA0003852310430000041
是全部隐藏状态;
注意力向量at计算方法为:
Figure BDA0003852310430000042
其中,
Figure BDA0003852310430000043
为上下文向量;ht为最后一个隐藏状态;Wc表示将要学习的模型参数。
进一步,将注意力向量at输入到全连接层,得到风速预测结果。
本发明的有益效果在于:
现有的单点风速组合预测模型普遍采用了全分解的数据预处理策略,这种数据预处理策略存在数据泄露的风险,可能导致此类组合预测模型缺少工程上的实用性,因此,准确预测风速的有效成分显得更为关键和可行。考虑到实测风速中存在大量的高频噪声,这部分高频噪声是风速随机性与波动性的重要来源,而风速的低频有效成分占据了风速的绝大部分能量,本发明通过无信息泄露的数据预处理方法提取风速序列的有效成分,建立未来信息无泄漏的表达方法,是得到可靠实用的风场预测结果的前提。
考虑到传统的组合预测模型普遍存在数据泄露问题,风速实测数据中也存在一定程度的噪音成分,本发明基于深度神经网络且无未来信息泄露的风速预测方法,采用无信息泄露的数据预处理策略的短期风速预测框架,旨在提供一种准确预测风速序列有效成分的方法。该预测方法包括数据筛选、基于实时滚动分解的数据预处理以及融合Attention机制的Bi-LSTM神经网络(Bi-LSTM-Attention)三部分内容,避免了数据泄露的风险,并对噪音成分进行了有效处理,且本发明方法采用事后评估,所以通过模态分解得到的真值是延迟获得的。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行说明:
图1为本发明基于深度神经网络且无未来信息泄露的风速预测方法的原理图;
图2为模型训练阶段、验证阶段和测试阶段(即使用阶段)数据的流动情况说明;
图3为所有模型的单步预测曲线图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明,以使本领域的技术人员可以更好的理解本发明并能予以实施,但所举实施例不作为对本发明的限定。
如图1所示,本发明基于深度神经网络且无未来信息泄露的风速预测方法,包括如下步骤。
步骤1:分别风速序列数据进行数据筛选和数据预处理;根据子模态能量占比大小筛选有效成分Tc并作为预测模型的输出值;采用实时滚动分解(Real-Time RollingDecomposition,RTRD)策略对风速序列数据进行预处理以得到无信息泄露的预测模型的输入值;以输入值和输出值构建数据集。
具体的,从风速序列数据中筛选出有效成分Tc的方法为:
S1:将风速序列数据划分训练集、验证集和测试集,本实施例中,风速序列数据按照6:2:2划分训练集(1~0.6N)、验证集(0.6N~0.8N)和测试集(0.8N~N),以划分后的数据组建三个数据组,第一个数据组包括训练集,第二个数据组包括训练集,第三个数据组包括测试集。
S2:对三个数据组中的风速序列数据分别进行模态分解。采用模态分解方法对第一个数据组(1~0.6N)、第二个数据组(0.6N~0.8N)和第三个数据组(0.8N~N)范围内的风速数据进行分解,其中VMD、SSA算法需指定分解层数,而EMD与ICEEMDAN算法能自适应确定分解层数。
S3:基于能量原理判断指标,选择合适的子模块得到能代表风速序列数据绝大部分能量的有效分成Tci,i=1,2或3,分别表示从三个数据组的风速序列数据中筛选得到的有效分成。具体的,根据子模态能量占比(Energy ratio,ER)大小筛选有效成分Tc需要的子模态数量(子模态按低频至高频的顺序排列)为:当前m个子模态能量之和大于等于原始风速序列99%的能量时,可认为Tci即为前m个子模态之和。子模态的ER计算方式为:
Figure BDA0003852310430000051
其中,ER(k)表示第k个子模态的能量占比;
Figure BDA0003852310430000052
表示第k个子模态序列;X表示风速序列数据;N表示风速序列数据X的长度。
S4:按顺序将有效成分Tci进行拼接,得到预测模型的目标输出值Tc,Tc是后续构建模型目标输出的数据来源,即真值。
具体的,采用实时滚动分解策略对风速序列数据进行预处理的方法为:
(1)确定窗口长度L:分析风速序列数据X的频谱特征,寻找频率峰值点,得到对应的周期,选择合适的周期数作为RTRD所需的窗口长度;
(2)确定滑动步幅s,当采取逐步预测的方式进行短期风速预测时,滑动步幅应设为1;
(3)重构风速信号:根据窗口长度L和滑动步幅s将风速序列数据X重构为二维矩阵XR,矩阵形状为(N-L+s,L),其中,N为风速序列数据X的长度;
(4)执行RTRD:逐行分解XR,得到无信息泄露的分解结果
Figure BDA0003852310430000053
Figure BDA0003852310430000054
每一行数据的前m个子模态相加,可以构成输入矩阵
Figure BDA0003852310430000055
计算
Figure BDA0003852310430000056
与Tc的相关系数;利用相关系数准则判断不同分解算法下构成
Figure BDA0003852310430000061
的最优子模态数量。
步骤2:确定预测步长:按照实际的预测需求,确定需要预测的未来时间步数。
步骤3:构建预测模型:构建结合注意力机制的双向长短期记忆网络(Bi-LSTM-Attention)作为预测模型,利用数据集对预测模型进行训练和测试。具体的,本实施例中,预测模型包括依次设置的一层输入层、两层Bi-LSTM网络、一层注意力层和一层全连接层。
本实施例的LSTM网络采用如下方式进行t时刻隐藏状态的更新:
ft=σ(Wxfxt+Whfht-1+bf)
gt=tanh(Wxgxt+Whght-1+bg)
it=σ(Wxixt+Whiht-1+bi)
ct=ft⊙ct-1+it⊙gt
ot=σ(Wxoxt+Whoht-1+bo)
yt=ht=ot⊙tanh(ct)
其中,xt为输入向量,ct为长期状态,ht为短期状态,yt为输出向量;Wxf,Wxg,Wxi,Wxo分别为与xt连接的权重矩阵,Whf,Whg,Whi,Who分别为与ht-1连接的权重矩阵,bf,bg,bi,bo分别为四层偏置项;ft为遗忘门的控制器,由ft决定ct-1的哪些部分应被删除;gt为LSTM网络的中间输出,其作用为分析xt与ht-1;it为输入门的控制器,判断gt的哪些重要部分应被添加进ct;ot为输出门的控制器,判断应读取ct中的哪部分并应用到ht和yt中;σ表示激活函数。
本实施例的注意力层采用Luong-Attention模型,其计算方法为:
score计算方法为:
Figure BDA0003852310430000062
注意力权重αts计算方法为:
Figure BDA0003852310430000063
其中,
Figure BDA0003852310430000064
和ht分别为编码器输出的全部隐藏状态以及最后一个隐藏状态;W表示权重矩阵;
上下文向量
Figure BDA0003852310430000065
计算方法为:
Figure BDA0003852310430000066
其中,αts是注意力权重,
Figure BDA0003852310430000071
是全部隐藏状态;
注意力向量at计算方法为:
Figure BDA0003852310430000072
其中,
Figure BDA0003852310430000073
为上下文向量;ht为最后一个隐藏状态;Wc表示将要学习到的模型参数。
最后,将注意力向量at输入到全连接层,得到风速预测结果。
步骤4:风速预测:利用训练得到的预测模型预测风速,得到风速预测结果。
图2说明了数据在神经网络的训练阶段、验证和测试阶段是如何运行的。在图2(a)中,当构造模型的输入时,模态分解遵循实时滚动分解策略,该策略在蓝色数据窗格上进行。这种策略保证了模型的输入始终使用预测值之前的信息,即未来信息的保密性得到保证,这意味着没有信息泄露。而预测值主要是指预测风速的有效分量,用Tc表示,遵循“真值”的含义。在图2(b)中,在模型训练过程中,由于数据已在完整的时间轴上收集,因此真值Tc1可以在训练神经网络时通过模态分解获得。因此,可以直接计算误差并用于更新神经网络的参数。在图2(b)中,在验证和测试过程中,神经网络的参数已经确定,预测结果与真实值之间的误差Tc2或Tc3以延迟方式计算,而不是实时计算。这意味着误差可以在一段时间的操作后计算,并且它们仅用于证明预测的准确性。总之,无论是在训练过程中还是在验证和测试过程中,都始终坚持实时滚动分解策略,避免信息泄露是我们算法的首要原则。
实验验证
下面结合具体实例验证本发明基于深度神经网络且无未来信息泄露的风速预测方法。
选取中国云南某测风塔70m高度处连续91天的10min间隔风速数据X作为输入信号,选取日期为2010.4.1~2010.6.30,共有13104组。
一、预测结果误差评估标准
本实施例选取了三种常用的性能指标作为衡量预测模型误差的判断标准,包括MAPE、MAE和RMSE,其定义分别如下。
Figure BDA0003852310430000074
Figure BDA0003852310430000075
Figure BDA0003852310430000076
其中:xn代表目标输出值(又称真值),
Figure BDA0003852310430000077
代表预测值。MAPE、MAE和RMSE的值越小,预测精度越高。
此外,定义RMSE、MAPE、MAE三个误差指标的相对提升比Pmetric,以进一步定量评价不同模型的性能表现。Pmetric的定义如下
Figure BDA0003852310430000081
其中,metric代表MAPE、MAE和RMSE三种指标,Ea,Eb代表模型a和模型b的预测误差。
二、预测模型与不同模型的对比
本实施例选择的深度学习平台为基于GPU的TensorFlow 2.3,Python 3.7版本,构建Bi-LSTM-Attention模型。训练过程采用“Adam”优化器,学习率lr设为0.001,损失函数选为“Mean Squared Error”,迭代次数固定为100个epoch。模型的各层参数由人工选择与随机搜索共同确定。
表1显示了持续模型(Persistence,PR)、BP模型、GRU、LSTM和本发明提出的Bi-LSTM-Attention进行了对比。表1和表2显示了不同预测模型单步预测和1-6步预测的RMSE、PMSE、MAPE、PMAPE、MAE和PMAE误差结果。
表1基于VMD方法的各模型预测性能—单步预测
Figure BDA0003852310430000082
表2基于VMD方法的各模型预测性能—1-6步预测
Figure BDA0003852310430000083
从以上所有神经网络模型的对比中可以看出,本发明方法能够准确的预测占有风速序列绝大部分能量的有效成分Tc,模型的单步和1-6步预测的精度均优于PR模型,且本发明所提出的Bi-LSTM-Attention预测模型取得了最佳的预测精度,这充分证明了本发明提出的RTRD-Bi-LSTM-Attention预测方法的有效性。具体来说,单步预测下,Bi-LSTM-Attention预测模块相对于PR模型的精度提升比Pmetric在23.79%~32.77%之间;1-6步预测下,Bi-LSTM-Attention预测模块相对于PR模型的精度提升比Pmetric在9.90%~18.09%之间。而当预测步长增大时,所有神经网络模型的Pmetric均有所下降,说明神经网络模型在多步风速预测时性能有所下降。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例,本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换,均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。

Claims (8)

1.一种基于深度神经网络且无未来信息泄露的风速预测方法,其特征在于:包括如下步骤:
步骤1:分别风速序列数据进行数据筛选和数据预处理;根据子模态能量占比大小筛选有效成分Tc并作为预测模型的输出值;采用实时滚动分解策略对风速序列数据进行预处理以得到无信息泄露的预测模型的输入值;以输入值和输出值构建数据集;
步骤2:确定预测步长:按照实际的预测需求,确定需要预测的未来时间步数;
步骤3:构建预测模型:构建结合注意力机制的双向长短期记忆网络作为预测模型,利用数据集对预测模型进行训练和测试;
步骤4:风速预测:利用训练得到的预测模型预测风速,得到风速预测结果。
2.根据权利要求1所述基于深度神经网络且无未来信息泄露的风速预测方法,其特征在于:所述步骤一中,从风速序列数据中筛选出有效成分Tc的方法为:
S1:将风速序列数据划分训练集、验证集和测试集,以划分后的数据组建三个数据组,第一个数据组包括训练集,第二个数据组包括验证集,第三个数据组包括测试集;
S2:对三个数据组中的风速序列数据分别进行模态分解;
S3:基于能量原理判断指标,选择合适的子模块得到能代表风速序列数据绝大部分能量的有效分成Tci,i=1,2或3,分别表示从三个数据组的风速序列数据中筛选得到的有效分成;
S4:按顺序将有效成分Tci进行拼接,得到预测模型的目标输出值Tc。
3.根据权利要求2所述基于深度神经网络且无未来信息泄露的风速预测方法,其特征在于:所述步骤S3中,当前m个子模态能量之和大于等于原始风速序列99%的能量时,Tci为前m个子模态之和,子模态的ER计算方式为:
Figure FDA0003852310420000011
其中,ER(k)表示第k个子模态的能量占比;
Figure FDA0003852310420000012
表示第k个子模态序列;X表示风速序列数据;N表示风速序列数据X的长度。
4.根据权利要求1所述基于深度神经网络且无未来信息泄露的风速预测方法,其特征在于:所述步骤一中,采用实时滚动分解策略对风速序列数据进行预处理的方法为:
(1)确定窗口长度L:分析风速序列数据X的频谱特征,寻找频率峰值点,得到对应的周期,选择周期数作为RTRD所需的窗口长度;
(2)确定滑动步幅s;
(3)重构风速信号:根据窗口长度L和滑动步幅s将风速序列数据X重构为二维矩阵XR,矩阵形状为(N-L+s,L),其中,N为风速序列数据X的长度;
(4)执行RTRD:逐行分解XR,得到无信息泄露的分解结果
Figure FDA0003852310420000021
5.根据权利要求1-4任一项所述基于深度神经网络且无未来信息泄露的风速预测方法,其特征在于:所述预测模型包括依次设置的一层输入层、两层Bi-LSTM网络、一层注意力层和一层全连接层。
6.根据权利要求5所述基于深度神经网络且无未来信息泄露的风速预测方法,其特征在于:所述LSTM网络的原理为:
ft=σ(Wxfxt+Whfht-1+bf)
gt=tanh(Wxgxt+Whght-1+bg)
it=σ(Wxixt+Whiht-1+bi)
ct=ft⊙ct-1+it⊙gt
ot=σ(Wxoxt+Whoht-1+bo)
yt=ht=ot⊙tanh(ct)
其中,xt为输入向量,ct为长期状态,ht为短期状态,yt为输出向量;Wxf,Wxg,Wxi,Wxo分别为与xt连接的权重矩阵,Whf,Whg,Whi,Who分别为与ht-1连接的权重矩阵,bf,bg,bi,bo分别为四层偏置项;ft为遗忘门的控制器,由ft决定ct-1的哪些部分应被删除;gt为LSTM网络的中间输出,其作用为分析xt与ht-1;it为输入门的控制器,判断gt的哪些重要部分应被添加进ct;ot为输出门的控制器,判断应读取ct中的哪部分并应用到ht和yt中;σ表示激活函数。
7.根据权利要求5所述基于深度神经网络且无未来信息泄露的风速预测方法,其特征在于:所述注意力层采用Luong-Attention模型,其计算方法为:
score计算方法为:
Figure FDA0003852310420000022
注意力权重αts计算方法为:
Figure FDA0003852310420000023
其中,
Figure FDA0003852310420000024
和ht分别为编码器输出的全部隐藏状态以及最后一个隐藏状态;W表示权重矩阵;
上下文向量
Figure FDA0003852310420000031
计算方法为:
Figure FDA0003852310420000032
其中,αts是注意力权重,
Figure FDA0003852310420000033
是全部隐藏状态;
注意力向量at计算方法为:
Figure FDA0003852310420000034
其中,
Figure FDA0003852310420000035
为上下文向量;ht为最后一个隐藏状态;Wc表示将要学到的模型参数。
8.根据权利要求6所述基于深度神经网络且无未来信息泄露的风速预测方法,其特征在于:将注意力向量at输入到全连接层,得到风速预测结果。
CN202211136465.0A 2022-09-19 2022-09-19 基于深度神经网络且无未来信息泄露的风速预测方法 Pending CN115330085A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211136465.0A CN115330085A (zh) 2022-09-19 2022-09-19 基于深度神经网络且无未来信息泄露的风速预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211136465.0A CN115330085A (zh) 2022-09-19 2022-09-19 基于深度神经网络且无未来信息泄露的风速预测方法

Publications (1)

Publication Number Publication Date
CN115330085A true CN115330085A (zh) 2022-11-11

Family

ID=83915084

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211136465.0A Pending CN115330085A (zh) 2022-09-19 2022-09-19 基于深度神经网络且无未来信息泄露的风速预测方法

Country Status (1)

Country Link
CN (1) CN115330085A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116662766A (zh) * 2023-08-01 2023-08-29 浙江大学 基于数据二维重构的风速预测方法及装置、电子设备

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116662766A (zh) * 2023-08-01 2023-08-29 浙江大学 基于数据二维重构的风速预测方法及装置、电子设备
CN116662766B (zh) * 2023-08-01 2023-10-03 浙江大学 基于数据二维重构的风速预测方法及装置、电子设备

Similar Documents

Publication Publication Date Title
CN109272146B (zh) 一种基于深度学习模型和bp神经网络校正的洪水预测方法
CN109272156B (zh) 一种超短期风电功率概率预测方法
CN111079998B (zh) 基于长短时序相关性注意力机制模型的流量预测方法
CN111144644B (zh) 基于变分异方差高斯过程回归的短期风速预测方法
CN109583588B (zh) 一种短期风速预测方法及***
CN113554466A (zh) 一种短期用电量预测模型构建方法、预测方法和装置
CN112016839B (zh) 一种基于qr-bc-elm的洪涝灾害预测预警方法
CN113887787A (zh) 一种基于长短时记忆网络和nsga-ii算法的洪水预报模型参数多目标优化方法
CN114169251A (zh) 一种超短期风电功率预测方法
CN112183877A (zh) 一种基于迁移学习的光伏电站故障智能诊断方法
CN115239029A (zh) 考虑功率时序及气象相依特性的风电功率预测方法及***
CN115330085A (zh) 基于深度神经网络且无未来信息泄露的风速预测方法
CN112767692A (zh) 基于SARIMA-GA-Elman组合模型的短时交通流量预测***
CN115796327A (zh) 一种基于vmd和iwoa-f-gru模型的风电功率区间预测方法
CN112561203B (zh) 一种基于聚类和gru实现水位预警的方法及***
CN115217152A (zh) 一种沉管隧道管节接头张合变形的预测方法及装置
CN115204362A (zh) 一种机械设备剩余寿命区间预测方法
CN114372615A (zh) 一种考虑空间相关性的短期光伏功率预测方法及***
Shen et al. An interval analysis scheme based on empirical error and MCMC to quantify uncertainty of wind speed
Zhang et al. A Novel Combined Model Based on Hybrid Data Decomposition, MSWOA and ENN for Short-Term Wind Speed Forecasting
Huang et al. Probabilistic prediction intervals of wind speed based on explainable neural network
CN117709488B (zh) 基于RUN-XGBoost的大坝渗流预测方法
Zheng et al. Research on predicting remaining useful life of equipment based on health index
LU504816B1 (en) Interpretable deep learning model for extreme climate prediction
CN117829822B (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