CN113193556A - 基于概率预测模型的短期风电功率预测方法 - Google Patents
基于概率预测模型的短期风电功率预测方法 Download PDFInfo
- Publication number
- CN113193556A CN113193556A CN202110517888.6A CN202110517888A CN113193556A CN 113193556 A CN113193556 A CN 113193556A CN 202110517888 A CN202110517888 A CN 202110517888A CN 113193556 A CN113193556 A CN 113193556A
- Authority
- CN
- China
- Prior art keywords
- prediction
- residual
- point
- function
- wind power
- 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
Links
Images
Classifications
-
- 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
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出一种基于概率预测模型的短期风电功率预测方法,即:1)由变分模态分解将风电功率历史数据分解为多个分量,针对每一维分量分别建立泄露积分型回声状态网络模型进行训练预测,将各个预测结果进行重构,得到风电功率点预测值;2)应用回声状态分位数回归网络对点预测的残差进行建模,得到不同分位数条件之下的残差预测值;3)将点预测值与残差预测值进行整合,在点预测的基础上由残差预测进一步提高预测精度,得到风电功率的概率预测值。本发明针对风电功率具有随机性、波动性特点,结合点预测模型与残差预测模型,得到概率预测模型,对风电功率进行精准的预测,对确保电力系统安全、经济、稳定的运行具有重大意义。
Description
技术领域
本发明涉及一种短期风电功率预测方法,具体地说,涉及一种基于概率预测模型的短期风电功率预测方法。本发明属于新能源发电过程中风电功率预测技术领域。
背景技术
自20世纪70年代世界性能源危机开始以来,许多国家都更加重视可再生能源的研究、开发与利用。我国也一直致力于优化电源结构,包括水能、风能、太阳能在内的可再生能源发电所占的比重越来越大,其中,风能具有蕴量巨大、分布广泛、利用技术成熟等优点,是国际公认的最具有大规模开发潜力的可再生能源之一,风能的开发利用已经成为我国可持续发展战略的重要组成部分,合理地开发利用风能在延缓全球气候变暖、保护生态环境和推进可持续发展等方面具有十分重要的意义。
然而,风电具有随机性、波动性,在并网时可能会对大电网带来难以调控等诸多问题,不仅会使系统的调度运行成本大大增高,还有可能给系统造成隐性危害,破坏系统的安全性和稳定性。因此,对风电功率进行准确有效地预测是电力系统重要的、不可缺少的环节,逐步提高风电功率预测精度成为当今的研究热点。
发明内容
针对风电功率具有随机性、波动性的特点,本发明的目的是提供一种基于概率预测模型的短期风电功率预测方法。
为实现上述目的,本发明采用以下技术方法:一种基于概率预测模型的短期风电功率预测方法,它包括以下步骤:
S1、由变分模态分解将风电功率历史数据分解为多个分量,针对每一维分量分别建立泄露积分型回声状态网络模型进行训练预测,将各个预测结果进行重构,得到风电功率点预测值;
S2、应用回声状态分位数回归网络对点预测的残差进行建模,得到不同分位数条件之下的残差预测值;
S3、将点预测值与残差预测值进行整合,在点预测的基础上由残差预测进一步提高预测精度,得到风电功率的概率预测值。
优选地,步骤S1中,所述变分模态分解将原始时间序列分解为K个固有模态分量,每个模态的带宽有限且具有不同的中心频率,使用交替方向乘子法和维纳滤噪法更新各模态分量的中心频率,同时将各模态函数解调到对应的基频带,最终实现所有模态分量重构之后近似等于原始时间序列和各模态分量带宽之和最小的目的;
原始时间序列被分解为K个分量,保证分解序列是具有中心频率的有限带宽的模态分量,同时各模态的估计带宽之和最小,约束条件为所有模态之和与原始信号相等,则相应约束表达式如下式所示:
式中,k为需要分解的模态个数;uk为分解出的第k个模函数分量;ωk为分解出的第k阶模函数的中心频率;δ(t)为狄拉克函数;*为卷积运算符;f(t)为原始信号;求解目标函数即可得到uk,ωk的解;
为求解目标函数,引入二次惩罚项α和拉格朗日算子λ将约束变分问题转化为非约束变分问题,其中α控制时间序列有噪声时的重构精度,λ用于调节约束的严格性;增广拉格朗日表达式如下:
式中,k为需要分解的模态个数;uk为分解出的第k个模函数分量;ωk为分解出的第k阶模函数的中心频率;δ(t)为狄拉克函数;*为卷积运算符;f(t)为原始信号;使用交替方向乘子法交替更新uk,ωk和λ即可求得式(1)的最优解。
优选地,步骤S1中,变分模态分解先将风电功率历史数据分解得到K个模态,然后对每一个模态都建立Leaky ESN进行训练预测,将各模态预测结果进行求和重构,得到风电功率的初步预测值;
所述Leaky ESN的储备池是由泄露积分型神经元组成,其神经元具有独立的状态动力学信息,利用各种方式来适应网络学习任务的时序特征;
假定输入层输入单元为u(k)=[u1(k),u2(k),…,uK(k)]T,储备池的神经状态为x(k)=[x1(k),x2(k),…,xN(k)]T,输出层输出单元为y(k)=[y1(k),y2(k),…,yL(k)]T;LeakyESN网络的状态更新方程为:
x(k+1)=(1-α)x(k)+αf(Winu(k+1)+Wx(k)+Wbacky(k))
式中,f(·)表示神经元激活函数;u(k+1)表示系统在k+1时刻的输入单元;x(k)表示储备池状态向量在k时刻的取值;y(k)表示系统在k时刻的输出单元;Win、W、Wback分别表示输入、储备池和反馈连接权矩阵;α是泄漏率;
ESN的输出方程为:
y(k+1)=fout(Wout[u(k+1),x(k+1),y(k)])
其中,y(k)和y(k+1)分别表示系统在k时刻和k+1时刻的输出单元;x(k+1)表示储备池状态向量在k+1时刻的取值;u(k+1)表示系统在k+1时刻的输入单元;fout(·)表示输出函数,根据问题的不同,fout(·)可以取线性函数或者S型函数;Wout表示输出连接权矩阵。
优选地,步骤S2中,应用回声状态分位数回归网络对点预测的残差进行建模,得到不同分位数条件之下的残差预测值的方法:
建立从输入层到隐含层的连接:
在上式基础上建立从隐含层到输出层的连接:
式中,为t时刻预测变量y的条件分位数估计;gj,t为t时刻输入层第j个神经元的输出;b(a)分别为隐含层到时输出层的连接权重、偏置;f(·)表示输出函数;τ为预测变量的第τ分位点;w(τ)、b(τ)分别为τ分位点下的连接权重、偏置。
优选地,步骤S2中,将分位数回归应用于ESN网络的训练过程,得到回声状态分位数回归网络;
所述ESN前向输出公式为:
y(t+1)=fout(Wout(u(k+1),x(k+1)))
式中,y(k+1)为系统在k+1时刻的输出单元;u(k+1)表示系统在k+1时刻的输入单元;x(k+1)表示储备池状态向量在k+1时刻的取值;fout(·)表示激活函数;Wout表示输出连接权矩阵;
将ESN输出层的输入向量定义为S(t):
S(t)=fconcate(u(k+1),x(k+1))
式中,fconcate为将两向量纵向拼接函数;u(k+1)表示系统在k+1时刻的输入单元;x(k+1)表示储备池状态向量在k+1时刻的取值;
则回声状态分位数回归网络为:
QY(τ|S)=fout(S,Wout(τ),b(τ))
式中,QY(τ|S)为预测变量Y关于输入向量S的τ分位点估计;fout为激活函数,通常选恒等函数;Wout(τ)为输出层的权重矩阵;b(τ)为输出层的偏置;
定义模型的损失函数为:
式中,ρτ为τ分位数对应的概率;Yi为与网络输出对应的实际值;S为网络输出层的输入向量;
为了进一步提高模型的泛化能力,在损失函数中添加正则项:
式中,f'cost为正则化后的损失函数;||·||2为求矩阵的第二范数;
训练后得到τ分位点的权重矩阵和偏置向量,通过模型的前向传播预测变量Y在τ分位点的条件分布预测:
优选地,步骤S3中,将两个子模型得到的点预测值与残差预测值进行整合,用残差预测值对点预测值进行修正,具体方法为:
式中,W和Xt表示点预测模型的参数和输入;f(·)表示输出函数;
S3.2.利用上述结果计算出预测残差:
S3.3将点预测值作为残差预测模型的输入特征,对残差进行预测:
优选地,每个模态分量都能用一个调幅-调频的模态函数uk(t)来表示,如:
优选地,输出连接权矩阵Wout,采用岭回归的方法获得,即:
Wout=YXT(XXT+θI)-1
式中,θ为正则化系数;I为单位矩阵;假设uteach(k)是训练样本的输入信号,yteach(k)是训练样本的期望输出信号;把所有的uteach(k)和x(k)储存起来并用矩阵X∈R(K+N)*T的行来表示,对应的目标输出值y(k)表示成Y∈RL*T的行向量。
本发明将点预测值与残差预测值进行整合,用残差预测值对点预测值进行修正,得到概率预测模型,该概率预测模型包括训练阶段和测试阶段,数据集包括点预测模型训练数据集、残差预测模型训练数据集、概率预测模型测试数据集;
在训练阶段,首先使用点预测模型训练数据集对点预测模型进行训练;应用此模型对残差预测模型训练数据集进行虚拟预测,进而结合点预测值和真实数据得到残差;接下来,将点预测结果与残差预测模型训练数据集中的历史数据和相关特征结合,训练残差预测模型;
在测试阶段,首先使用概率预测模型测试数据集中的历史数据和相关特征进行点预测;然后使用残差预测模型对残差数据进行建模,其中将点预测值与历史数据和其他相关特征一起用作输入特征,最后将点预测值与残差预测值合并以产生概率预测值。
本发明与现有技术相比,具有如下优点:
(1)结合了变分模态分解和泄露积分性回声状态网络得到点预测模型,可有效克服风电功率具有随机性、波动性的问题,对风电功率进行较为准确的初步预测。
(2)本发明针对残差数据具有波动大的特点,将分位数回归应用于回声状态网络的训练过程中,得到回声状态分位数回归网络,作为残差预测模型对残差进行预测,可以得到在不同分位数下的多个预测结果,形成一个预测区间,充分的反应残差的波动趋势,可以良好的解决上述问题。
(3)本发明将点预测值和残差预测值整合,得到最终的短期风电功率概率预测值,由残差预测值对点预测值进行修正,在点预测的基础上进一步提高预测精度。
附图说明
图1是本发明基于概率预测模型预测短期风电功率的方法流程图;
图2是本发明点预测模型流程图;
图3是本发明泄露积分型回声状态网络结构图;
图4是本发明概率预测模型框架图。
具体实施方式
以下结合附图对本发明作进一步说明。
如图1所示,本发明提出了一种基于概率预测模型的短期风电功率预测方法,它包括以下步骤:
S1、由变分模态分解将风电功率历史数据分解为多个分量,针对每一维分量分别建立泄露积分型回声状态网络模型进行训练预测,将各个预测结果进行重构,得到风电功率点预测值。
为了克服风电功率具有随机性、波动性的问题,对风电功率进行较为准确的初步预测。本发明选取变分模态分解和泄露积分性回声状态网络相结合的方法构建点预测模型,得到风电功率的点预测值。如图2所示,风电功率历史数据经过变分模态分解处理得到K个模态,对每一个模态都建立Leaky ESN进行训练预测,将各模态预测结果进行求和重构,得到风电功率的初步预测值。
S11、变分模态分解
由于风电功率数据往往具有较大的噪音,在使用一些数据前需要适当的处理算法来降噪。降噪的方法通常有小波变换、局部均值分解、模态分解等方法。近年来,随着变分模态分解的提出和研究,越来越多的风电功率预测研究开始使用变分模态分解来进行风电预测数据的预处理。变分模态分解采用非递归的求解变分的方法得到模态分量,有效的避免了模态混叠、过包络和欠包络等问题,对于噪声较大的信号具有很高的适应性。
变分模态分解能把原始时间序列分解为K个固有模态分量,每个模态的带宽有限且具有不同的中心频率,使用交替方向乘子法和维纳滤噪法更新各模态分量的中心频率,同时将各模态函数解调到对应的基频带,最终实现所有模态分量重构之后近似等于原始时间序列和各模态分量带宽之和最小的目的。
每个模态分量都能用一个调幅-调频的模态函数uk(t)来表示,如:
原始信号被分解为K个分量,保证分解序列是具有中心频率的有限带宽的模态分量,同时各模态的估计带宽之和最小,约束条件为所有模态之和与原始信号相等,则相应约束表达式如下式所示:
式中,k为需要分解的模态个数;uk为分解出的第k个模函数分量;ωk为分解出的第k阶模函数的中心频率;δ(t)为狄拉克函数;*为卷积运算符;f(t)为原始信号。求解目标函数即可得到uk,ωk的解。
为求解目标函数,引入二次惩罚项α和拉格朗日算子λ将约束变分问题转化为非约束变分问题,其中α控制时间序列有噪声时的重构精度,λ用于调节约束的严格性。增广拉格朗日表达式如下:
式中,k为需要分解的模态个数;uk为分解出的第k个模函数分量;ωk为分解出的第k阶模函数的中心频率;δ(t)为狄拉克函数;*为卷积运算符;f(t)为原始信号。使用交替方向乘子法交替更新uk,ωk和λ即可求得式(1)的最优解。
S12、泄露积分型回声状态网络
Leaky ESN网络是ESN网络的一种改进模型,其储备池是由泄露积分型神经元组成。这种类型的神经元具有独立的状态动力学信息,可以利用各种方式来适应网络学习任务的时序特征。Leaky ESN与ESN的拓扑结构相同,如图3所示。
左侧为K个输入节点;中间是储备池网络,由N个内部节点以及稀疏的节点连接权值构成;右侧是L个输出节点。实线表示网络的必要连接,而虚线表示不同情况下可能存在的连接。假定输入层输入单元为u(k)=[u1(k),u2(k),…,uK(k)]T,储备池的神经状态为x(k)=[x1(k),x2(k),…,xN(k)]T,输出层输出单元为y(k)=[y1(k),y2(k),…,yL(k)]T。LeakyESN网络的状态更新方程为:
x(k+1)=(1-α)x(k)+αf(Winu(k+1)+Wx(k)+Wbacky(k))
式中,f(·)表示神经元激活函数;u(k+1)表示系统在k+1时刻的输入单元;x(k)表示储备池状态向量在k时刻的取值;y(k)表示系统在k时刻的输出单元;Win、W、Wback分别表示输入、储备池和反馈连接权矩阵;α是泄漏率。
由上式可以看出,ESN是Leaky ESN在α=1时的特例,其储备池是由泄露积分型神经元组成。这种类型的神经元具有独立的状态动力学信息,可以利用各种方式来适应网络学习任务的时序特征。改变一个参数对结果的影响可以由另一个参数进行弥补,并不影响网络的回声状态特性。因此,只要参数选择合适,Leaky ESN的储备池性能优于基本ESN。
ESN的输出方程为:
y(k+1)=fout(Wout[u(k+1),x(k+1),y(k)])
其中,y(k)和y(k+1)分别表示系统在k时刻和k+1时刻的输出单元;x(k+1)表示储备池状态向量在k+1时刻的取值;u(k+1)表示系统在k+1时刻的输入单元;fout(·)表示输出函数,根据问题的不同,fout(·)可以取线性函数或者S型函数;Wout表示输出连接权矩阵,采用岭回归的方法获得Wout,即:
Wout=YXT(XXT+θI)-1
式中,θ为正则化系数;I为单位矩阵。假设uteach(k)是训练样本的输入信号,yteach(k)是训练样本的期望输出信号。把所有的uteach(k)和x(k)储存起来并用矩阵X∈R(K+N)*T的行来表示,对应的目标输出值y(k)表示成Y∈RL*T的行向量。
Leaky ESN的数学模型可以看作是作用在储备池状态神经元上的低通滤波器。泄漏率α控制前一时刻神经元状态的保持程度,其截止频率由参数α决定。较小的α取值导致内部神经元状态x(k)的较慢变化,从而明显增强了ESN的短时记忆能力。特别值得一提的是,Leaky ESN在训练阶段只改变输出连接权值,而其他权值固定不变。
S2、应用回声状态分位数回归网络对点预测的残差进行建模,得到不同分位数条件之下的残差预测值。
残差数据具有波动大的特点,普通的预测模型难以反应残差数据的变化规律。而分位数回归模型对数据的分布情况没有严格限制,主要是描述响应变量在不同分位数下的变化规律,突出局部间的关联性。它可以得到在不同分位数下的多个预测结果,形成一个预测区间,以此对数据的波动进行反应,可以良好的解决上述问题。同时ESN具有较好的时间序列拟合能力,同时兼具快速的训练速度。因此,本发明将分位数回归应用于回声状态网络的训练过程中,得到回声状态分位数回归网络,构成残差预测模型,得到残差预测值。
S21、分位数回归神经网络
分位数回归神经网络(QRNN)模型是一种非参数的分位数回归方法。对于简单的单隐层前馈神经网络,建立从输入层到隐含层的连接:
在上式基础上建立从隐含层到输出层的连接:
式中,为t时刻预测变量y的条件分位数估计;gj,t为t时刻输入层第j个神经元的输出;b(a)分别为隐含层到时输出层的连接权重、偏置;f(·)表示输出函数;τ为预测变量的第τ分位点;w(τ)、b(τ)分别为τ分位点下的连接权重、偏置。
QRNN通过使用分位数回归使得神经网络具有概率性预测能力,训练出的参数具有大样本理论下的渐进优良性。
S22、回声状态分位数回归网络
与全连接神经网络相比,ESN具有较好的拟合能力,具有快速的训练速度。因此,本发明提出将分位数回归应用于ESN网络的训练过程,得到回声状态分位数回归网络(ESQRN)。ESQRN的数学推导公式如下。
ESN前向输出公式为:
y(t+1)=fout(Wout(u(k+1),x(k+1)))
式中,y(k+1)为系统在k+1时刻的输出单元;u(k+1)表示系统在k+1时刻的输入单元;x(k+1)表示储备池状态向量在k+1时刻的取值;fout(·)表示激活函数;Wout表示输出连接权矩阵。
将ESN输出层的输入向量定义为S(t):
S(t)=fconcate(u(k+1),x(k+1))
式中,fconcate为将两向量纵向拼接函数;u(k+1)表示系统在k+1时刻的输入单元;x(k+1)表示储备池状态向量在k+1时刻的取值。
则回声状态分位数回归网络为:
QY(τ|S)=fout(S,Wout(τ),b(τ))
式中,QY(τ|S)为预测变量Y关于输入向量S的τ分位点估计;fout为激活函数,通常选恒等函数;Wout(τ)为输出层的权重矩阵;b(τ)为输出层的偏置。
定义模型的损失函数为:
式中,ρτ为τ分位数对应的概率;Yi为与网络输出对应的实际值;S为网络输出层的输入向量。
为了进一步提高模型的泛化能力,在损失函数中添加正则项:
式中,fc'ost为正则化后的损失函数;||·||2为求矩阵的第二范数。
训练后得到τ分位点的权重矩阵和偏置向量,通过模型的前向传播预测变量Y在τ分位点的条件分布预测:
由于利用了ESN的较好的拟合能力和分位数回归的优点,ESQRN能在较少计算资源的开销下完成准确的残差概率预测,同时模型对异常值和异方差分布具有较强的鲁棒性。
S3、将点预测值与残差预测值进行整合,在点预测的基础上由残差预测进一步提高预测精度,得到风电功率的概率预测值。
将上述两个子模型得到的点预测值与残差预测值进行整合,用残差预测值对点预测值进行修正,得到概率预测模型,具体预测框架如图4所示,该概率预测模型包括训练阶段和测试阶段,数据集包括点预测模型训练数据集T1、残差预测模型训练数据集T2、概率预测模型测试数据集T3;
在训练阶段,首先使用点预测模型训练数据集T1对点预测模型进行训练;应用此模型对残差预测模型训练数据集T2进行虚拟预测,进而结合点预测值和真实数据得到残差εt;接下来,将点预测结果与残差预测模型训练数据集T2中的历史数据和相关特征结合,训练残差预测模型;
在测试阶段,首先使用概率预测模型测试数据集T3中的历史数据和相关特征进行点预测;然后使用残差预测模型对残差数据进行建模,其中将点预测值与历史数据和其他相关特征一起用作输入特征,最后将点预测值与残差预测值合并以产生概率预测值。
对于某一个时刻的概率预测模型预测结果是一个区间范围,会更有效的反应预测风电功率的波动范围,提供未来数据更多的不确定性信息。
式中,W和Xt表示点预测模型的参数和输入;f(·)表示输出函数。
利用上述结果可以计算出预测残差:
综上所述,本发明用残差预测值对点预测值进行修正,将点预测值与残差预测值进行整合,得到风电功率的概率预测值,在点预测的基础上由残差预测进一步提高预测精度。这对规划电力系统的运行调度计划,确保系统的安全、经济、稳定运行具有重大意义。
以上所述的实施例仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (9)
1.一种基于概率预测模型的短期风电功率预测方法,其特征在于,包括以下步骤:
S1、由变分模态分解将风电功率历史数据分解为多个分量,针对每一维分量分别建立泄露积分型回声状态网络模型进行训练预测,将各个预测结果进行重构,得到风电功率点预测值;
S2、应用回声状态分位数回归网络对点预测的残差进行建模,得到不同分位数条件之下的残差预测值;
S3、将点预测值与残差预测值进行整合,在点预测的基础上由残差预测进一步提高预测精度,得到风电功率的概率预测值。
2.如权利要求1所述的一种基于概率预测模型的短期风电功率预测方法,其特征在于,步骤S1中,所述变分模态分解将原始时间序列分解为K个固有模态分量,每个模态的带宽有限且具有不同的中心频率,使用交替方向乘子法和维纳滤噪法更新各模态分量的中心频率,同时将各模态函数解调到对应的基频带,最终实现所有模态分量重构之后近似等于原始时间序列和各模态分量带宽之和最小的目的;
原始时间序列被分解为K个分量,保证分解序列是具有中心频率的有限带宽的模态分量,同时各模态的估计带宽之和最小,约束条件为所有模态之和与原始信号相等,则相应约束表达式如下式所示:
式中,k为需要分解的模态个数;uk为分解出的第k个模函数分量;ωk为分解出的第k阶模函数的中心频率;δ(t)为狄拉克函数;*为卷积运算符;f(t)为原始信号;求解目标函数即可得到uk,ωk的解;
为求解目标函数,引入二次惩罚项α和拉格朗日算子λ将约束变分问题转化为非约束变分问题,其中α控制时间序列有噪声时的重构精度,λ用于调节约束的严格性;增广拉格朗日表达式如下:
式中,k为需要分解的模态个数;uk为分解出的第k个模函数分量;ωk为分解出的第k阶模函数的中心频率;δ(t)为狄拉克函数;*为卷积运算符;f(t)为原始信号;使用交替方向乘子法交替更新uk,ωk和λ即可求得式(1)的最优解。
3.如权利要求1所述的一种基于概率预测模型的短期风电功率预测方法,其特征在于,步骤S1中,变分模态分解先将风电功率历史数据分解得到K个模态,然后对每一个模态都建立Leaky ESN进行训练预测,将各模态预测结果进行求和重构,得到风电功率的初步预测值;
所述Leaky ESN的储备池是由泄露积分型神经元组成,其神经元具有独立的状态动力学信息,利用各种方式来适应网络学习任务的时序特征;
假定输入层输入单元为u(k)=[u1(k),u2(k),…,uK(k)]T,储备池的神经状态为x(k)=[x1(k),x2(k),…,xN(k)]T,输出层输出单元为y(k)=[y1(k),y2(k),…,yL(k)]T;Leaky ESN网络的状态更新方程为:
x(k+1)=(1-α)x(k)+αf(Winu(k+1)+Wx(k)+Wbacky(k))
式中,f(·)表示神经元激活函数;u(k+1)表示系统在k+1时刻的输入单元;x(k)表示储备池状态向量在k时刻的取值;y(k)表示系统在k时刻的输出单元;Win、W、Wback分别表示输入、储备池和反馈连接权矩阵;α是泄漏率;
ESN的输出方程为:
y(k+1)=fout(Wout[u(k+1),x(k+1),y(k)])
其中,y(k)和y(k+1)分别表示系统在k时刻和k+1时刻的输出单元;x(k+1)表示储备池状态向量在k+1时刻的取值;u(k+1)表示系统在k+1时刻的输入单元;fout(·)表示输出函数,根据问题的不同,fout(·)可以取线性函数或者S型函数;Wout表示输出连接权矩阵。
4.如权利要求1所述的一种基于概率预测模型的短期风电功率预测方法,其特征在于,步骤S2中,应用回声状态分位数回归网络对点预测的残差进行建模,得到不同分位数条件之下的残差预测值的方法:
建立从输入层到隐含层的连接:
在上式基础上建立从隐含层到输出层的连接:
5.如权利要求1所述的一种基于概率预测模型的短期风电功率预测方法,其特征在于,步骤S2中,将分位数回归应用于ESN网络的训练过程,得到回声状态分位数回归网络;
所述ESN前向输出公式为:
y(t+1)=fout(Wout(u(k+1),x(k+1)))
式中,y(k+1)为系统在k+1时刻的输出单元;u(k+1)表示系统在k+1时刻的输入单元;x(k+1)表示储备池状态向量在k+1时刻的取值;fout(·)表示激活函数;Wout表示输出连接权矩阵;
将ESN输出层的输入向量定义为S(t):
S(t)=fconcate(u(k+1),x(k+1))
式中,fconcate为将两向量纵向拼接函数;u(k+1)表示系统在k+1时刻的输入单元;x(k+1)表示储备池状态向量在k+1时刻的取值;
则回声状态分位数回归网络为:
QY(τ|S)=fout(S,Wout(τ),b(τ))
式中,QY(τ|S)为预测变量Y关于输入向量S的τ分位点估计;fout为激活函数,通常选恒等函数;Wout(τ)为输出层的权重矩阵;b(τ)为输出层的偏置;
定义模型的损失函数为:
式中,ρτ为τ分位数对应的概率;Yi为与网络输出对应的实际值;S为网络输出层的输入向量;
为了进一步提高模型的泛化能力,在损失函数中添加正则项:
式中,f′cost为正则化后的损失函数;||·||2为求矩阵的第二范数;
训练后得到τ分位点的权重矩阵和偏置向量,通过模型的前向传播预测变量Y在τ分位点的条件分布预测:
6.如权利要求1所述的一种基于概率预测模型的短期风电功率预测方法,其特征在于,步骤S3中,将两个子模型得到的点预测值与残差预测值进行整合,用残差预测值对点预测值进行修正,具体方法为:
式中,W和Xt表示点预测模型的参数和输入;f(·)表示输出函数;
S3.2.利用上述结果计算出预测残差:
S3.3将点预测值作为残差预测模型的输入特征,对残差进行预测:
8.如权利要求3所述的一种基于概率预测模型的短期风电功率预测方法,其特征在于,输出连接权矩阵Wout,采用岭回归的方法获得,即:
Wout=YXT(XXT+θI)-1
式中,θ为正则化系数;I为单位矩阵;假设uteach(k)是训练样本的输入信号,yteach(k)是训练样本的期望输出信号;把所有的uteach(k)和x(k)储存起来并用矩阵X∈R(K+N)*T的行来表示,对应的目标输出值y(k)表示成Y∈RL*T的行向量。
9.如权利要求6所述的一种基于概率预测模型的短期风电功率预测方法,其特征在于,将点预测值与残差预测值进行整合,用残差预测值对点预测值进行修正,得到概率预测模型,该概率预测模型包括训练阶段和测试阶段,数据集包括点预测模型训练数据集、残差预测模型训练数据集、概率预测模型测试数据集;
在训练阶段,首先使用点预测模型训练数据集对点预测模型进行训练;应用此模型对残差预测模型训练数据集进行虚拟预测,进而结合点预测值和真实数据得到残差;接下来,将点预测结果与残差预测模型训练数据集中的历史数据和相关特征结合,训练残差预测模型;
在测试阶段,首先使用概率预测模型测试数据集中的历史数据和相关特征进行点预测;然后使用残差预测模型对残差数据进行建模,其中将点预测值与历史数据和其他相关特征一起用作输入特征,最后将点预测值与残差预测值合并以产生概率预测值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110517888.6A CN113193556B (zh) | 2021-05-12 | 2021-05-12 | 基于概率预测模型的短期风电功率预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110517888.6A CN113193556B (zh) | 2021-05-12 | 2021-05-12 | 基于概率预测模型的短期风电功率预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113193556A true CN113193556A (zh) | 2021-07-30 |
CN113193556B CN113193556B (zh) | 2023-02-03 |
Family
ID=76981256
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110517888.6A Active CN113193556B (zh) | 2021-05-12 | 2021-05-12 | 基于概率预测模型的短期风电功率预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113193556B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113610665A (zh) * | 2021-08-12 | 2021-11-05 | 东北电力大学 | 一种基于多延迟输出回声状态网的风力发电功率预测方法 |
CN113673768A (zh) * | 2021-08-24 | 2021-11-19 | 河海大学 | 基于vmd和cnn-lstm的风功率概率预测模型 |
CN113988481A (zh) * | 2021-12-23 | 2022-01-28 | 南京鼐威欣信息技术有限公司 | 一种基于动态矩阵预测控制的风功率预测方法 |
CN114021297A (zh) * | 2021-11-18 | 2022-02-08 | 吉林建筑科技学院 | 基于回声状态网络的复杂管网泄漏定位方法 |
CN115456312A (zh) * | 2022-11-09 | 2022-12-09 | 广东电网有限责任公司肇庆供电局 | 基于辛几何模态分解的短期电力负荷预测方法及系统 |
CN115498660A (zh) * | 2022-11-17 | 2022-12-20 | 华北电力大学 | 风电场调频动态建模方法、装置及电子设备 |
CN116663865A (zh) * | 2023-07-31 | 2023-08-29 | 山东建筑大学 | 考虑源荷不确定性的综合能源系统日前调度方法及系统 |
CN117892113A (zh) * | 2024-03-13 | 2024-04-16 | 广东工业大学 | 一种自适应vmd和双重降维注意力机制的风电功率预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103117546A (zh) * | 2013-02-28 | 2013-05-22 | 武汉大学 | 一种超短期风电功率滑动预测方法 |
CN109242212A (zh) * | 2018-10-16 | 2019-01-18 | 中国矿业大学 | 一种基于变模式分解与长短记忆网络的风电预测方法 |
CN109272156A (zh) * | 2018-09-12 | 2019-01-25 | 河海大学 | 一种超短期风电功率概率预测方法 |
CN111612275A (zh) * | 2020-05-29 | 2020-09-01 | 云南电网有限责任公司 | 一种预测区域用户负荷量的方法及装置 |
CN112270454A (zh) * | 2020-11-19 | 2021-01-26 | 国网北京市电力公司 | 极端因素影响下的电力系统短期负荷预测方法与预测装置 |
-
2021
- 2021-05-12 CN CN202110517888.6A patent/CN113193556B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103117546A (zh) * | 2013-02-28 | 2013-05-22 | 武汉大学 | 一种超短期风电功率滑动预测方法 |
CN109272156A (zh) * | 2018-09-12 | 2019-01-25 | 河海大学 | 一种超短期风电功率概率预测方法 |
CN109242212A (zh) * | 2018-10-16 | 2019-01-18 | 中国矿业大学 | 一种基于变模式分解与长短记忆网络的风电预测方法 |
CN111612275A (zh) * | 2020-05-29 | 2020-09-01 | 云南电网有限责任公司 | 一种预测区域用户负荷量的方法及装置 |
CN112270454A (zh) * | 2020-11-19 | 2021-01-26 | 国网北京市电力公司 | 极端因素影响下的电力系统短期负荷预测方法与预测装置 |
Non-Patent Citations (4)
Title |
---|
严雪晴: "基于灰色预测模型的广东省货运总量预测研究", 《数学的实践与认识》 * |
牛安敏 等: "基于VMD-Leaky-ESN的电力系统短期负荷多步预测方法", 《新疆大学学报》 * |
王新友等: "基于改进型ESN的短期风电功率预测", 《自动化与仪器仪表》 * |
范锦谕: "基于深度学习算法的含风电场区域的等效负荷短期预测", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113610665A (zh) * | 2021-08-12 | 2021-11-05 | 东北电力大学 | 一种基于多延迟输出回声状态网的风力发电功率预测方法 |
CN113610665B (zh) * | 2021-08-12 | 2024-01-30 | 东北电力大学 | 一种基于多延迟输出回声状态网的风力发电功率预测方法 |
CN113673768A (zh) * | 2021-08-24 | 2021-11-19 | 河海大学 | 基于vmd和cnn-lstm的风功率概率预测模型 |
CN114021297B (zh) * | 2021-11-18 | 2023-12-01 | 吉林建筑科技学院 | 基于回声状态网络的复杂管网泄漏定位方法 |
CN114021297A (zh) * | 2021-11-18 | 2022-02-08 | 吉林建筑科技学院 | 基于回声状态网络的复杂管网泄漏定位方法 |
CN113988481B (zh) * | 2021-12-23 | 2022-05-03 | 南京鼐威欣信息技术有限公司 | 一种基于动态矩阵预测控制的风功率预测方法 |
CN113988481A (zh) * | 2021-12-23 | 2022-01-28 | 南京鼐威欣信息技术有限公司 | 一种基于动态矩阵预测控制的风功率预测方法 |
CN115456312A (zh) * | 2022-11-09 | 2022-12-09 | 广东电网有限责任公司肇庆供电局 | 基于辛几何模态分解的短期电力负荷预测方法及系统 |
CN115498660A (zh) * | 2022-11-17 | 2022-12-20 | 华北电力大学 | 风电场调频动态建模方法、装置及电子设备 |
US11847530B1 (en) | 2022-11-17 | 2023-12-19 | North China Electric Power University | Frequency modulation dynamic modeling method and device for wind farm, and electronic device |
CN116663865A (zh) * | 2023-07-31 | 2023-08-29 | 山东建筑大学 | 考虑源荷不确定性的综合能源系统日前调度方法及系统 |
CN116663865B (zh) * | 2023-07-31 | 2023-11-14 | 山东建筑大学 | 考虑源荷不确定性的综合能源系统日前调度方法及系统 |
CN117892113A (zh) * | 2024-03-13 | 2024-04-16 | 广东工业大学 | 一种自适应vmd和双重降维注意力机制的风电功率预测方法 |
CN117892113B (zh) * | 2024-03-13 | 2024-06-11 | 广东工业大学 | 一种自适应vmd和双重降维注意力机制的风电功率预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113193556B (zh) | 2023-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113193556B (zh) | 基于概率预测模型的短期风电功率预测方法 | |
Mi et al. | Wind speed prediction model using singular spectrum analysis, empirical mode decomposition and convolutional support vector machine | |
He et al. | A hybrid system for short-term wind speed forecasting | |
Ma et al. | A hybrid attention-based deep learning approach for wind power prediction | |
Li et al. | A ship motion forecasting approach based on empirical mode decomposition method hybrid deep learning network and quantum butterfly optimization algorithm | |
CN111461463B (zh) | 一种基于tcn-bp的短期负荷预测方法、系统及设备 | |
CN111832825B (zh) | 融合长短期记忆网络和极限学习机的风电预测方法及系统 | |
Chen et al. | Model-free primal-dual methods for network optimization with application to real-time optimal power flow | |
Li et al. | A hybrid deep interval prediction model for wind speed forecasting | |
Zou et al. | A hybrid model based on multi-stage principal component extraction, GRU network and KELM for multi-step short-term wind speed forecasting | |
Wang et al. | An improved Wavenet network for multi-step-ahead wind energy forecasting | |
Alamaniotis | Synergism of deep neural network and elm for smart very-short-term load forecasting | |
CN117114160A (zh) | 一种短期光伏功率预测方法 | |
Duong-Ngoc et al. | Short term load forcast using deep learning | |
CN112508279A (zh) | 基于空间相关性的区域分布式光伏预测方法及系统 | |
CN116843083A (zh) | 基于混合神经网络模型的碳排放预测系统及方法 | |
Wu et al. | Explainable temporal dependence in multi-step wind power forecast via decomposition based chain echo state networks | |
Alharbi et al. | Short-term solar irradiance forecasting model based on bidirectional long short-term memory deep learning | |
Huang et al. | Short-term PV power forecasting based on CEEMDAN and ensemble DeepTCN | |
Qi et al. | Short-term prediction of PV power based on fusions of power series and ramp series | |
Liu et al. | A hybrid deep learning model based on parallel architecture TCN-LSTM with Savitzky-Golay filter for wind power prediction | |
Zhu et al. | Wind Speed Short-Term Prediction Based on Empirical Wavelet Transform, Recurrent Neural Network and Error Correction | |
Gao et al. | TgDLF2. 0: Theory-guided deep-learning for electrical load forecasting via transformer and transfer learning | |
Qader | Power management in a hydrothermal system considering maintenance using Lagrangian relaxation and augmented Lagrangian methods | |
CN116822722A (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 |