一种基于VMD-IPSO-GRU的电力负荷预测方法技术

一种基于VMD

【技术实现步骤摘要】
一种基于VMD
‑
IPSO
‑
GRU的电力负荷预测方法


[0001]本专利技术涉及电力负荷预测
，尤其涉及一种电力负荷预测方法。

技术介绍

[0002]能源是我们人类赖以生存的基础与保障，其中电力能源对我国来说是最重要的能源之一。因此，保障电力的稳定供应，对我国的经济及发展、居民的工作与生活以及社会的稳定有至关重要的作用。目前无论是人们生活或是生产的工具，如通讯工具、生产设备等大多数都是由电力驱动的。一旦电力系统无法保证稳定，进而导致停电事故的发生，那么该地区将直接处于瘫痪状态，人们的通勤和休息被打乱，公司的运作被干扰，工厂的生产也将无法进行。因此对电力负荷进行预测，保证电力系统的正确运行一直是相关领域的研究热点。
[0003]传统的电力负荷预测方法包括时间序列法、回归分析法等，这些模型结构简单，但对电力负荷的预测精度也较低。而随机森林、人工神经网络、支持向量回归等机器学习方法没有考虑到负荷数据在时序上的相关性，RNN循环神经网络因其具有时序性，被广泛运用在时间序列的预测上。GRU和LSTM都是一种特殊的RNN，相对于LSTM，GRU结构更简单、收敛速度更快，但利用GRU进行预测时需要依靠经验来对部分参数进行确定，具有不确定性。因此本专利技术使用了改进的粒子群算法对这些参数进行寻优来进行电力负荷的预测。

技术实现思路

[0004]针对现有电力负荷预测技术中存在的上述不足，本专利技术提供一种基于VMD
‑
IPSO
‑
GRU的电力负荷预测方法。
[0005]本专利技术采用VMD变分模态分解将原始数据分解为多个分量，使用改进的粒子群算法IPSO对GRU网络的参数进行优化，该模型具有预测速度快、准确性高、泛化能力强的特点。
[0006]本专利技术的目的是提出一种基于VMD
‑
IPSO
‑
GRU的电力负荷预测方法，包括如下步骤：
[0007]步骤1：获取某地区的电力负荷数据及该地区的天气数据，保留其中的日期、电力消耗数据和温度数据，构成原始数据集；
[0008]步骤2：对原始电力负荷时间序列及气象数据、时间数据进行数据清洗，补全缺失值，剔除偏差过大的异常值；
[0009](2.1)缺失数据处理：检查原始数据集中是否存在缺失值，由于电力负荷具有连续性，因此对缺失数据使用均值法进行填充，即对该数据缺失点的前后两条数据之和求平均值；
[0010](2.2)偏差数据处理：由于电力负荷呈现周期性，且短期内的变化幅度较小，因此对识别到的异常数据中的偏差数据使用一周前同一采样位置的数据替代；
[0011](2.3)归一化处理：根据公式(1)对数据进行归一化处理
[0012][0013]其中y
i
为归一化后的值，x
i
为数据原始值值，x
max
和x
min
分别为数据的最大值和最小值。
[0014]步骤3：获得最佳分解模态数K，将电力负荷时间序列f(t)进行VMD分解，得到K个不同中心频率的模态u
k
＝(u1,u2,
…
,u
k
)，将VMD分解后的子序列与温度因素T结合作为输入向量I
k
＝(u
k
,q1,q2,
…
,q
n
)，同时划分训练集和测试集。
[0015](3.1)首先是VMD的构建，通过对每一个模态进行希尔伯特变换得到其解析信号和单边谱，根据梯度平方准则和高斯平滑度计算确定其带宽，计算过程如公式(2)所示
[0016][0017]式中，k为模态的分解个数，为对时刻t的微分，δ(t)为狄拉克函数，u
k
和ω
k
为分解后的第k个模态分量的数学表达式和中心频率，|| ||2为二范式函数，f(t)为待处理的信号。最终分解得到的所有模态之和必须与f(t)一致。
[0018](3.2)其次是VMD的求解。引入拉格朗日算子，将每个模态确定带宽的约束问题转变为无约束问题进行求解。处理过程如公式(3)所示：
[0019][0020]其中，α为二次惩罚因子，λ为拉格朗日乘法算子。
[0021]使用交替方向乘子法更新各个模态分量u
k
(t)和对应的中心频率ω
k
，求得公式(2)的最优解，通过傅里叶变换，计算分解得到的k个模态分量的最优值，如公式(4)所示:
[0022][0023]按照相同的过程，得到中心频率的更新公式，如公式(5)所示：
[0024][0025]由此，通过VMD分解出电力负荷数据的k个模态分量u
k
和中心频率ω
k
；
[0026]步骤4：建立IPSO
‑
GRU预测模型。
[0027](4.1)设置IPSO的维度空间d＝2，粒子群规模为20,最大迭代次数为100，惯性权重ω
max
＝0.9，ω
min
＝0.4，学习因子c1＝c2＝1.5；
[0028](4.2)设置隐藏层节点数n0、n1，取值范围为[1,200]；
[0029](4.3)初始化粒子群中的每个粒子i的位置速度个体最优解群体最
优解其中粒子i的位置代表节点数为n0,n1；
[0030](4.4)将每个粒子的位置参数作为神经网络节点训练GRU模型，计算适应度值，并比较选择适应度小的位置更新
[0031](4.5)对粒子群中的每个粒子的惯性权重、速度与位置进行更新，具体如公式(6)、(7)、(8)所示：
[0032][0033][0034]其中，分别代表粒子在第k+1、k次迭代中第d维的速度、位置；为第k次迭代时第d维的个体最优解、群体最优解；r1、C2为加速度常数，分别代表粒子的个体学习因子和群体学习因子，r1r2为[0,1]内的随机数，用于满足粒子寻优过程的随机性；
[0035][0036]其中，ω为惯性权重，ω
max
和ω
min
分别表示惯性权重的最大值和最小值，k为当前迭代次数，k
max
为最大迭代次数；
[0037](4.6)重复步骤(4.3)、(4.4)，直到满足最大迭代次数时停止，将寻找到的全局最优解代入GRU进行训练；
[0038]步骤5：将所有子序列的预测值序列重构并进行反归一化处理，得到最终预测值并根据平均绝对百分比误差MAPE来判断，其计算公式如公式(9)所示：
[0039][0040]其中，y
i
表示的中真实值，表示神经网络模型的预测值。
[0041]将本专利技术提出的方法与其他经典预测方法在同一数据集上进行准确度对比，使用单一GRU模型、PSO
‑
GRU模型、IPSO
‑
GRU模型、VMD本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于VMD
‑
IPSO
‑
GRU的电力负荷预测方法，包括如下步骤：步骤1：获取某地区的电力负荷数据及该地区的天气数据，保留其中的日期、电力消耗数据和温度数据，构成原始数据集；步骤2：对原始电力负荷时间序列及气象数据、时间数据进行数据清洗，补全缺失值，剔除偏差过大的异常值；(2.1)缺失数据处理：检查原始数据集中是否存在缺失值，由于电力负荷具有连续性，因此对缺失数据使用均值法进行填充，即对该数据缺失点的前后两条数据之和求平均值；(2.2)偏差数据处理：由于电力负荷呈现周期性，且短期内的变化幅度较小，因此对识别到的异常数据中的偏差数据使用一周前同一采样位置的数据替代；(2.3)归一化处理：根据公式(1)对数据进行归一化处理其中y
i
为归一化后的值，x
i
为数据原始值值，x
max
和x
min
分别为数据的最大值和最小值；步骤3：获得最佳分解模态数K，将电力负荷时间序列f(t)进行VMD分解，得到K个不同中心频率的模态u
k
＝(u1,u2,
…
,u
k
)，将VMD分解后的子序列与温度因素T结合作为输入向量I
k
＝(u
k
,q1,q2,
…
,q
n
)，同时划分训练集和测试集；(3.1)首先是VMD的构建，通过对每一个模态进行希尔伯特变换得到其解析信号和单边谱，根据梯度平方准则和高斯平滑度计算确定其带宽，计算过程如公式(2)所示式中，k为模态的分解个数，为对时刻t的微分，δ(t)为狄拉克函数，u
k
和ω
k
为分解后的第k个模态分量的数学表达式和中心频率，|| ||2为二范式函数，f(t)为待处理的信号；最终分解得到的所有模态之和必须与f(t)一致；(3.2)其次是VMD的求解；引入拉格朗日算子，将每个模态确定带宽的约束问题转变为无约束问题进行求解；处理过程如公式(3)所示：其中，α为二次惩罚因子，λ为拉格朗日乘法算子；使用交替方向乘子法更新各个模态分量u
k
(t)和对应的中心频率ω
k
，求得公式(2)的最优解，通过傅里叶变换，计算分解得到的k个模态分量的最优值，如公式(4)所示:
按照相同的过程，得到中心频率的更新公式，如公式(5)所示：由此，通过VMD分解出电力负荷数据的k个模态分量u
k
和中心频率ω
k
；步骤4：建立IPSO
‑
GRU预测模型；(4.1)设置IPSO的维度空...

【专利技术属性】
技术研发人员：蒋一波，朱航锟，周泽宝，周轲，
申请(专利权)人：浙江工业大学，
类型：发明
国别省市：