基于VMD-PSO-BPNN的短期电力负荷预测模型建立方法技术

本发明专利技术公开了一种基于VMD‑PSO‑BPNN的短期电力负荷预测模型建立方法，该短期电力负荷预测模型应用于造纸企业的电力负荷预测，包括下述步骤：首先获取造纸企业总有效负荷数据；利用VMD分解算法对预处理后的数据序列进行分解；利用滞后自相关选择模型的输入变量；利用PSO‑BPNN对分解序列进行建模；利用训练样本对PSO‑BPNN模型进行训练，建立预测模型，并进行造纸企业用电负荷预测，并对预测效果进行分析。本发明专利技术方法建立的基于VMD‑PSO‑BPNN的短期电力负荷预测模型，具有收敛快，预测结果无滞后的特点。

【技术实现步骤摘要】
基于VMD-PSO-BPNN的短期电力负荷预测模型建立方法
本专利技术涉及造纸企业智能用电
，具体涉及一种基于VMD-PSO-BPNN的短期电力负荷预测模型建立方法。
技术介绍
制浆造纸过程存在大量间歇用电设备，通过优化这些间歇用电设备的开机、停机计划，可有效实现用电削峰填谷，降低成产成本。目前制浆造纸企业对电力需求侧的预测主要通过历史经验购电，这导致所购的常出现过剩或者不足的问题。因此，如何对企业用电负荷进行预测，进而实现优化排产和合理购电是提高造纸企业经济效益增强稳定生产的一个重要途径。近三年来，主要的研究方向之一是采用拆分重构的方法对电力负荷进行预测。由于造纸企业的用电无周期性，且负荷波动频率高，利用VMD可以将高频率的波形拆分成几个不同且频率较低的波形的特点。而且相比于EMD(经典模态分解)，VMD算法是一种更为精确的数学算法，而且对噪声和采样频率都很敏感，因此十分适合于造纸企业用电负荷的拆分。神经网络模型根据不同的实际情况来选择适合解决相应问题的结构参数。但是当需要解决的问题相对比较复杂时，用基本的人工神经网络(ANN)一般达不到实际应用的需求，基于优化算法进行神经网络优化的结合算法是解决上述问题的方法之一。BPNN(BP神经网络)是一种针对非线性、非周期、无规律、无结构性或半结构性数据建模最常用、效果最佳的模型，结合数据挖掘建立具有时间序列特征的BPNN来建立造纸企业用电负荷预测模型是十分实用的。PSO(粒子群优化算法)是一类概率型的全局优化算法。非确定算法的优点在于算法能有更多机会求解全局最优解。由于BPNN预测存在过拟合以及陷入局部最优等问题，因此采用PSO算法优化BPNN的权值和阈值，防止以上问题的发生是十分合适的。
技术实现思路
本专利技术的目的是为了解决现有技术中的上述缺陷，提供一种基于VMD-PSO-BPNN的短期电力负荷预测模型建立方法，利用VMD分解算法把造纸企业用电负荷进行分解，利用滞后自相关法选择预测模型的输入，最后利用PSO对BPNN的权值和阈值进行优化，建立电力负荷预测模型，预测并对预测效果进行评价，从而能够精确预测造纸企业用电负荷。本专利技术的目的可以通过采取如下技术方案达到：一种基于VMD-PSO-BPNN的短期电力负荷预测模型建立方法，该短期电力负荷预测模型应用于造纸企业的电力负荷预测，所述的建立方法包括以下步骤：S1、获取造纸企业数据的用电数据；S2、利用VMD分解算法，对预处理后的负荷序列进行序列分解；S3、利用滞后自相关法对每个分解序列选取输入变量；S4、设置BPNN网络的隐藏层神经元数，以及BPNN网络的权值和阈值，把分解序列的训练集输入初始的BPNN网络中，把拟合结果和实际结果之间残差作为适应度值，利用PSO算法更新权值和阈值的大小，寻找最优的拟合结果，把最优拟合结果对应BPNN网络进行输出，利用训练好的BPNN网络对分解序列进行预测，把所有分解序列的预测结果进行叠加，得到短期电力负荷预测模型。进一步地，所述的步骤S2中通过搜寻约束分变模型最优解来实现信号自适应分解，将原始负荷序列分解成一系列具有稀疏特性的模态分量，即将原始序列分解为不同频率的序列，具体包括：S201、对于每个模态u(t)，通过希尔伯特变换计算与之相关的解析信号，计算公式如下：式中，H(t)为模态解析信号，δ(t)是狄拉克分布，t是采样时间点，j为虚数，*表示卷积；S202、对各模态解析信号预估计的中心频率ωk进行混合，将每个模态的频谱调制到相应的基频带，其公式如下：S203、计算式(2)中基频带的梯度平方L2范数，估计出各模态分量的带宽，对应的约束变分模型为：式中，f(t)＝∑ku(t)；S204、采用二次惩罚函数项和拉格朗日乘子算子得到一个无约束问题，最后求解该问题的公式为：式中，{uk}＝{u1,u2,···,uk}代表分解得到的k个IMF分量，{ωk}＝{ω1,ω2,···,ωk}表示各分量的中心频率，∑K表示各模态分量求和，λ(t)为拉格朗日乘数，α是数据保真约束的平衡参数，f(t)为原始信号。进一步地，所述的步骤S3中，通过滞后自相关方法找出过去用电负荷对当前用电负荷的影响，使用自相关函数作为选择信息特征子集的指导，即通过自相关的滞后阶次来选取输入变量，当滞后自相关系数的绝对值大于0.8时，用这一滞后时刻对应的有效功率作为模型的输入，其表达式为：式中，是给定时间序列中所有X的平均值，X＝{Xt:t∈T}，为时间序列数据集，rk为测量时间t和t-k时间序列的线性相关性。进一步地，所述的步骤S4包括：S401、把数据处理后的输入和输出数据集分为训练集和测试集；S402、初始化BPNN的权值和阈值，并设置BPNN参数，其中BPNN参数包括：输入层层数、隐藏层层数、输出层的层数、训练次数、训练目标和学习速率；S403、初始化PSO优化算法参数，PSO优化算法参数包括：学习因子c1和c2、粒子群规模N、惯性权重wmax和wmin、最大迭代次数和最大速度Vmax；S404、设置当前迭代次数i＝1；S405、训练BPNN，并计算各个粒子的适应度，对适应度的大小进行排序，把每个粒子作为当前种群的局部最优，记为Ppbest，把粒子群中适应度最小的粒子作为全局最优，记为Pgbest；S406、更新粒子的速度和位置；S407、判断是否满足结束的条件，如果不满足，迭代次数i＝i+1，并重新执行步骤S405-步骤S407，直到满足结束条件；S408、输出最终的权值和阈值给BPNN，得到PSO-BPNN模型；S409、输入处理好的测试集的输入变量数据，通过PSO-BPNN模型进行预测，输出预测负荷曲线。进一步地，所述的每一层的状态只影响下一层神经元的状态，其中隐含层和输出层的激励函数选取双曲正弦Sigmoid函数，学习算法选择动量梯度下降算法函数traingdm。进一步地，所述的步骤S406中对第k次迭代粒子i的第n维的速度和位置的更新变化公式如下，其中，1≤n≤N：式中：是第k次迭代粒子i位置矢量的第n维分量，xi＝(xi1，xi2，···，xiN)，范围限定在[Xmin,n，Xmax,n]内，是第k次迭代粒子i飞行速度矢量的第n维分量，Vi＝(vi1,vi2,…,viN)，范围限定在[-Vmax,n，Vmax,n]内，pbestin代表粒子i位置矢量的第n维分量自身经历过的最好位置，粒子i自身经历最好的位置记为，pbesti＝(pbesti1,pbesti2,…,pbestin)，gbestn代表群体的第n维分量的所有微粒经历过的最好位置，群体所有粒子经过的最好位置记为，gbest＝(gbest1，gbest2，···，gbestn)，c1、c2是加速度常数，调节学习最大步长，rand()是随机函数，取值范围[0,1]，以增加搜索随机性，ω是惯性因子，非负数。进一步地，所述的步骤S405中选择误差作为适应度。进一步地，所述的建立方法还包括模型评价步骤：S5、对预测模型的预测结果进行预测效果分析，过程如下：S501、将预测模型得出的用电负荷的预测值和实际值进行比较；S502、根据预测模型预测效果的评价指标进行模型预测效果的分析；其中，模型预测效果的评价指标包括均方根误差平方和RMSE、平均绝对本文档来自技高网...

【技术保护点】
1.一种基于VMD‑PSO‑BPNN的短期电力负荷预测模型建立方法，该短期电力负荷预测模型应用于造纸企业的电力负荷预测，其特征在于，所述的建立方法包括以下步骤：S1、获取造纸企业数据的用电数据；S2、利用VMD分解算法，对预处理后的负荷序列进行序列分解；S3、利用滞后自相关法对每个分解序列选取输入变量；S4、设置BPNN网络的隐藏层神经元数，以及BPNN网络的权值和阈值，把分解序列的训练集输入初始的BPNN网络中，把拟合结果和实际结果之间残差作为适应度值，利用PSO算法更新权值和阈值的大小，寻找最优的拟合结果，把最优拟合结果对应BPNN网络进行输出，利用训练好的BPNN网络对分解序列进行预测，把所有分解序列的预测结果进行叠加，得到短期电力负荷预测模型。

【技术特征摘要】
1.一种基于VMD-PSO-BPNN的短期电力负荷预测模型建立方法，该短期电力负荷预测模型应用于造纸企业的电力负荷预测，其特征在于，所述的建立方法包括以下步骤：S1、获取造纸企业数据的用电数据；S2、利用VMD分解算法，对预处理后的负荷序列进行序列分解；S3、利用滞后自相关法对每个分解序列选取输入变量；S4、设置BPNN网络的隐藏层神经元数，以及BPNN网络的权值和阈值，把分解序列的训练集输入初始的BPNN网络中，把拟合结果和实际结果之间残差作为适应度值，利用PSO算法更新权值和阈值的大小，寻找最优的拟合结果，把最优拟合结果对应BPNN网络进行输出，利用训练好的BPNN网络对分解序列进行预测，把所有分解序列的预测结果进行叠加，得到短期电力负荷预测模型。2.根据权利要求1所述的基于VMD-PSO-BPNN的短期电力负荷预测模型建立方法，其特征在于，所述的步骤S2中通过搜寻约束分变模型最优解来实现信号自适应分解，将原始负荷序列分解成一系列具有稀疏特性的模态分量，即将原始序列分解为不同频率的序列，具体包括：S201、对于每个模态u(t)，通过希尔伯特变换计算与之相关的解析信号，计算公式如下：式中，H(t)为模态解析信号，δ(t)是狄拉克分布，t是采样时间点，j为虚数，*表示卷积；S202、对各模态解析信号预估计的中心频率ωk进行混合，将每个模态的频谱调制到相应的基频带，其公式如下：S203、计算式(2)中基频带的梯度平方L2范数，估计出各模态分量的带宽，对应的约束变分模型为：式中，f(t)＝∑ku(t)；S204、采用二次惩罚函数项和拉格朗日乘子算子得到一个无约束问题，最后求解该问题的公式为：式中，{uk}＝{u1,u2,···,uk}代表分解得到的k个IMF分量，{ωk}＝{ω1,ω2,···,ωk}表示各分量的中心频率，∑K表示各模态分量求和，λ(t)为拉格朗日乘数，α是数据保真约束的平衡参数，f(t)为原始信号。3.根据权利要求1所述的基于VMD-PSO-BPNN的短期电力负荷预测模型建立方法，其特征在于，所述的步骤S3中，通过滞后自相关方法找出过去用电负荷对当前用电负荷的影响，使用自相关函数作为选择信息特征子集的指导，即通过自相关的滞后阶次来选取输入变量，当滞后自相关系数的绝对值大于0.8时，用这一滞后时刻对应的有效功率作为模型的输入，其表达式为：式中，是给定时间序列中所有X的平均值，X＝{Xt:t∈T}，为时间序列数据集，rk为测量时间t和t-k时间序列的线性相关性。4.根据权利要求1所述的基于VMD-PSO-BPNN的短期电力负荷预测模型建立方法，其特征在于，所述的步骤S4包括：S401、把数据处理后的输入和输出数据集分为训练集和测试集；S402、初始化BPNN的权值和阈值，并设置BPNN参数，其中BPNN参数包括：输入层层数、隐藏层层数、输出层的层数、训练次数、训练目标和学习速率；S403、初始化PSO优化算法参数，PSO优化算法参数包括：学习因子c1和c2、粒子群规模N、惯性权...

【专利技术属性】
技术研发人员：洪蒙纳，李继庚，满奕，胡雨沙，
申请(专利权)人：华南理工大学，
类型：发明
国别省市：广东,44