CN114648147A - 一种基于ipso-lstm的风电功率预测方法 - Google Patents

Abstract

本发明涉及一种基于IPSO‑LSTM的风电功率预测方法。该方法首先对风电功率数据进行清洗以及归一化处理，然后确定LSTM初始网络结构，根据IPSO算法优化更新LSTM的参数，最后确定最佳LSTM网络结构，最后使用训练好的模型进行风电功率短期预测，最后将模型预测结果进行反归一化并输出。该方法充分融合了LSTM强大的预测功能以及IPSO的优化能力，提高了模型的准确率。本发明通过具体实验验证了方法的有效性。

Description

一种基于IPSO-LSTM的风电功率预测方法

技术领域

本发明属于电力系统时序数据预测技术领域，具体地说是一种基于IPSO-LSTM的风电功率预测方法。

背景技术

随着能源与环境问题的日益突出，风能作为清洁能源中最重要的可生能源之一，已受到世界各国的高度重视。国内风力发电水平也保持着稳步增长，2018年风力发电水平达到了3660亿kW·h，占全部发电量的5.2％。大规模的风电接入给电力系统的稳定运行和调度带来影响，若能对风电功率进行可靠和准确的预测，则有利于电网的稳定运行。

针对风电功率的预测方法，风电功率的预测模型主要包括时间序列模型和基于数值天气预报的预测模型，后者又包括物理模型和统计模型。在统计模型中，有基于LSTM网络的短期风力功率预测模型。还有一些作者提出基于主成分分析和遗传算法来优化BP神经网络，进行风电场短期功率的预测。还有采用CNN和GRU相结合的网络进行超短期风电预测。同时，还有一些采用人工智能算法和深度学习神经网络相结合的方法。例如，采用一种基于改进的萤火虫算法的短期风电功率预测和采用一种使用遗传算法改进LSTM的短期风电功率预测方法。

针对粒子群算法的应用，主要用于对单目标或多目标的优化，可以利用粒子群算法对深度学习网络进行参数优化。LSTM作为一种改进的时间递归神经网络(RNN)，可以学习时间序列长短期信息，可以用来处理和预测时间序列中的间隔和延迟事件。同时，也存在一些缺点，例如神经网络隐藏层数目、时间步长、批处理大小难以确定。在实际应用中，这些参数都是依靠经验来确定，具有很大的随机性。可以使用IPSO来优化LSTM的参数，进而来更好地进行风电功率预测。

发明内容

针对现有风电功率预测方法中存在的不足，本发明解决的技术问题是根据风电功率历史数据来预测未来一段时间内的风电功率。采用改进粒子群算法优化的LSTM模型来进行风电功率预测，提高了模型的准确率和鲁棒性。

本发明为实现上述目的所采用的技术方案是：一种基于IPSO-LSTM的风电功率预测方法，包括以下步骤：

1)从电力系统监测站采集风电相关数据；

2)对于采集的风电功率数据进行数据清洗，并将数据划分为训练集和测试集；

3)分别对训练集和测试集数据进行归一化处理；

4)用训练集数据输入IPSO-LSTM模型用于优化网络参数，用测试集数据进行测试，从而获取优化后的IPSO-LSTM模型；

5)采集实际现场一段时间内的风电功率相关数据序列，并输入优化后的IPSO-LSTM模型，获取风电功率数据的预测结果；

6)对风电功率数据的预测结果再进行反归一化处理，结合时间点、风速得到预测的未来设定时间内风电功率的数据。

所述风电相关数据为时间点、风速、风电功率的数据。

所述数据清洗具体为：使用孤立森林算法进行异常值检测，对于异常值的处理，采用附近的正常值进行修正，对于数据集中的缺失值，采用均值的方法进行填充。

所述数据归一化处理具体为：采用最小-最大标准化对样本数据进行线性变换，将结果值映射到[0,1]之间；

其中，样本数据为风电功率数据，转换表达式为：

x'＝(x-x_min)/(x_max-x_min)

上式中，x、x'分别为样本数据的实际值和归一化后的值，x_min、x_max分别为样本数据集中的最小值和最大值。

所述网络参数需要预先设置其范围，包括设置粒子群粒子个数，最大迭代次数，神经网络隐藏层单元个数取值范围[a,b]，时间步长取值范围[1,p],批处理大小取值范围[1,d]；其中，a,b,d,p均为自然数。

是通过改进的粒子群算法IPSO不断的迭代更新上述参数，最终确定出LSTM网络隐藏层单元个数，时间步长，批处理大小的最优值，从而获取优化后的IPSO-LSTM模型。

所述改进的粒子群算法为：从惯性权重和加速因子两个方面解决粒子容易陷入局部最优解和收敛速度慢的问题，惯性权重表达式和自调整加速因子表达式分别表示为：

其中，ω_max，ω_min分别是惯性权重的最大值、最小值，c_1s和c_2s分别为c₁和c₂的初值，c_1c和c_2c分别是c₁和c₂的终值，t为当前迭代的次数，T为总迭代的次数。

所述LSTM网络结构包括：遗忘门、输入门、输出门；LSTM网络在t时刻保持单元状态的记忆，通过遗忘门f_t和输入门i_t进行调节；遗忘门的作用是让细胞记住或忘记之前的状态C_t-1，输入门的作用是允许或阻止传入信号更新单元状态；输出门的作用是控制单元状态C输出和传输到下一个单元；将归一化后的训练集样本数据依次经过遗忘门、输入门、输出门处理后输出当前数据的预测数据。

所述LSTM网络模型是由多个单元结构组成，该结构能通过更新内部状态来长时间存储信息，函数如下：

f_t＝σ(w_f·[h_t-1,x_t]+b_f)

i_t＝σ(w_i·[h_t-1,x_t]+b_i)

o_t＝σ(w_o·[h_t-1,x_t]+b_o)

h_t＝o_t·tanh(C_t)

式中，x为LSTM单元的输入向量；h为LSTM单元的输出向量，f、i、o分别表示遗忘门、输入门和输出门；C表示单元状态；下标t表示时刻；σ、tanh分别为sigmoid、tanh激活函数；w和b分别表示权重和偏差矩阵。

本发明具有以下有益效果及优点：

1)本发明为风电功率预测提供一种切实有效的方法，能够较为准确的预测出未来四天内风电功率的变化趋势，为电力系统的调度，维护提供有效的数据支持。

2)本发明以风电功率的历史数据作为模型的输入，充分发挥LSTM对时间序列数据预测的高效性，并采用粒子群算法对参数进行优化，提高了模型的准确度。

附图说明

图1为本发明一种基于IPSO-LSTM风电功率预测方法结构图；

图2为LSTM单元结构；

图3为不同模型预测结果图；

具体实施方式

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明的具体实施方法做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但本发明能够以很多不同于在此描述的其他方式来实施，本领域技术人员可以在不违背发明内涵的情况下做类似改进，因此本发明不受下面公开的具体实施的限制。

除非另有定义，本文所使用的所有技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本发明。

本发明提出一种基于IPSO-LSTM的风电功率预测模型，LSTM模型具有很强的预测能力，能针对时间序列数据，具有很好的效果。粒子群方法易于实现，收敛速度快，使用IPSO来优化LSTM神经网络参数的选择，从而提高模型的预测精度。

如图1所示，基于IPSO-LSTM的风电功率预测方法，包括以下步骤：

步骤1：电力系统监测站有风电功率监控装置，从监测站采集风电功率数据以及相应的时间点，风速，等相关数据；

步骤2：对采集到的风电功率数据进行数据清洗，包括对缺失值和异常值处理；

步骤3：将数据划分为训练集和测试集，然后将归一化后的数据作为LSTM模型的输入。

步骤4：根据IPSO算法不断迭代更新出LSTM的最优参数。

步骤5：对训练好的模型进行测试和预测，并将模型预测结果进行反归一化处理并输出，得到在未来一段时间的风电功率预测值。

所述数据清洗具体为：使用孤立森林算法进行异常值检测，对于异常值的处理，采用附近的正常值进行修正。对于数据集中的缺失值，采用均值的方法进行填充。

所述数据归一化处理具体为：采用最小-最大标准化对样本数据进行线性变换，将结果值映射到[0,1]之间，样本数据为风电功率数据，转换表达式为：

x'＝(x-x_min)/(x_max-x_min)

上式中，x、x'分别为按照步骤2处理后的样本数据实际值和归一化后的值，x_min、x_max分别为数据集中的最小值和最大值。

传统粒子群算法描述：在一个n维的搜索空间中，由m个粒子构成种群X，X＝{x₁,...,x_i,...,x_m}，其中第i个粒子的位置可描述为x_i＝(x_i1,x_i2,...,x_in)^T，第i个粒子的速度可描述为v_i＝(v_i1,v_i2,...,v_in)^T。粒子的局部最优解为p_i＝(p_i1,p_i2,...p_in)^T，种群的全局最优解为p_g＝(p_g1,p_g2,...,p_gn)^T。在迭代寻优的过程中，每个粒子以其经历过的最优适应值处的解作为局部最优解，以所有粒子经历过最优适应值处的解作为全局最优解，并根据这两个最优解来更新粒子位置

和速度

其表达式分别为：

式中，d＝1,2,...n；i＝1,2,...m；m为种群规模；t为当前的进化代数；r₁和r₂为0到1之间的随机数，c₁和c₂为加速因子，一般取值为c₁＝c₂＝2，ω为惯性权重。

所述改进的粒子群算法(IPSO)：由于传统粒子群算法存在寻优时间长、容易陷入局部最优的问题，该方法主要从惯性权重和加速因子两个方面，解决粒子容易陷入局部最优解和收敛速度慢的问题。惯性权重表达式和自调整加速因子表达式分别表示为：

其中，ω_max，ω_min分别是惯性权重的最大值、最小值，c_1s和c_2s分别为c₁和c₂的初值，c_1c和c_2c分别是c₁和c₂的终值，t为当前迭代的次数，T为总迭代的次数。

所述LSTM网络：

f_t＝σ(w_f·[h_t-1,x_t]+b_f)

i_t＝σ(w_i·[h_t-1,x_t]+b_i)

o_t＝σ(w_o·[h_t-1,x_t]+b_o)

h_t＝o_t·tanh(C_t)

如图2所示，LSTM是由多个单元结构组成，该结构能通过更新内部状态来长时间存储信息，式中，x为LSTM单元的输入向量；h为LSTM单元的输出向量，f、i、o分别表示遗忘门、输入门和输出门；C表示单元状态；下标t表示时刻；σ、tanh分别为sigmoid、tanh激活函数；w和b分别表示权重和偏差矩阵。

LSTM的关键是单元状态C,它在t时刻保持单元状态的记忆，通过遗忘门f_t和输入门i_t进行调节。遗忘门的作用是让细胞记住或忘记之前的状态C_t-1，输入门的作用是允许或阻止传入信号更新单元状态。输出门的作用是控制单元状态C输出和传输到下一个单元。

IPSO-LSTM模型初始化：设置粒子群粒子个数为40，最大迭代次数为600，神经网络隐藏层单元个数取值范围[10,50]，时间步长取值范围[1,40],批处理大小取值范围[1,50]。

通过不断的迭代更新，最终确定出LSTM网络隐藏层单元个数为30，时间步长为5，批处理大小为10。

实例分析1：

步骤1数据采集：采集某地区的电力系统监测点的风电功率数据，每一小时进行一次数据的采集，共采集12天数据，即从2019年3月1日至3月12日，共计288小时。

步骤2数据清洗：使用孤立森林算法进行异常值检测，对于异常值的处理，采用附近的正常值进行修正。对于数据集中的缺失值，采用均值的方法进行填充。同时进行归一化处理，并将处理后的数据划分为训练集和测试集。以2019年3月1日至2019年3月8日的数据作为模型的训练集，2019年3月9日至2019年3月12日的数据作为测试集。

步骤3IPSO-LSTM模型初始化：设置粒子群粒子个数为40，最大迭代次数为600，神经网络隐藏层单元个数取值范围[10,50]，时间步长取值范围[1,40],批处理大小取值范围[1,50]。

步骤4通过不断的迭代更新，最终确定出LSTM网络隐藏层单元个数为30，时间步长为5，批处理大小为10。

步骤5模型预测：使用训练好的模型根据历史数据并进行未来四天的风电功率。

如图3所示，为充分表明模型预测精度的提高，本实验分别采用LSTM模型，PSO-LSTM模型，GA-LSTM模型以及IPSO-LSTM模型进行训练和预测，并计算模型预测的均方根误差(RMSE)，平均绝对误差(MAE)，平均绝对百分比误差(MAPE)。

表1不同模型的误差对比

由表1可知，本发明方法IPSO-LSTM模型比其他三种模型的精度都要高。以平均绝对百分比误差为例，IPSO-LSTM模型比LSTM模型提高了1.79％，比PSO-LSTM模型提高了1.51％，比GA-LSTM模型提高了0.79％。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种基于IPSO-LSTM的风电功率预测方法，其特征在于，包括以下步骤：

S1、从电力系统监测站采集风电相关数据；

S2、对于采集的风电功率数据进行数据清洗，并将数据划分为训练集和测试集；

S3、分别对训练集和测试集数据进行归一化处理；

S4、用训练集数据输入IPSO-LSTM模型用于优化网络参数，用测试集数据进行测试，从而获取优化后的IPSO-LSTM模型；

S5、采集实际现场一段时间内的风电功率相关数据序列，并输入优化后的IPSO-LSTM模型，获取风电功率数据的预测结果；

S6、对风电功率数据的预测结果再进行反归一化处理，结合时间点、风速得到预测的未来设定时间内风电功率的数据。

2.根据权利要求1所述的一种基于IPSO-LSTM的风电功率预测方法，其特征在于，所述风电相关数据为时间点、风速、风电功率的数据。

3.根据权利要求1所述的一种基于IPSO-LSTM的风电功率预测方法，其特征在于，所述数据清洗具体为：使用孤立森林算法进行异常值检测，对于异常值的处理，采用附近的正常值进行修正，对于数据集中的缺失值，采用均值的方法进行填充。

4.根据权利要求1所述的一种基于IPSO-LSTM的风电功率预测方法，其特征在于，所述数据归一化处理具体为：采用最小-最大标准化对样本数据进行线性变换，将结果值映射到[0,1]之间；

其中，样本数据为风电功率数据，转换表达式为：

x'＝(x-x_min)/(x_max-x_min)

上式中，x、x'分别为样本数据的实际值和归一化后的值，x_min、x_max分别为样本数据集中的最小值和最大值。

5.根据权利要求1所述的一种基于IPSO-LSTM的风电功率预测方法，其特征在于，所述网络参数需要预先设置其范围，包括设置粒子群粒子个数，最大迭代次数，神经网络隐藏层单元个数取值范围[a,b]，时间步长取值范围[1,p],批处理大小取值范围[1,d]；其中，a,b,d,p均为自然数。

6.根据权利要求5所述的一种基于IPSO-LSTM的风电功率预测方法，其特征在于，是通过改进的粒子群算法IPSO不断的迭代更新上述参数，最终确定出LSTM网络隐藏层单元个数，时间步长，批处理大小的最优值，从而获取优化后的IPSO-LSTM模型。

7.根据权利要求6所述的一种基于IPSO-LSTM的风电功率预测方法，其特征在于，所述改进的粒子群算法为：从惯性权重和加速因子两个方面解决粒子容易陷入局部最优解和收敛速度慢的问题，惯性权重表达式和自调整加速因子表达式分别表示为：

其中，ω_max，ω_min分别是惯性权重的最大值、最小值，c_1s和c_2s分别为c₁和c₂的初值，c_1c和c_2c分别是c₁和c₂的终值，t为当前迭代的次数，T为总迭代的次数。

8.根据权利要求1所述的一种基于IPSO-LSTM的风电功率预测方法，其特征在于，所述LSTM网络结构包括：遗忘门、输入门、输出门；LSTM网络在t时刻保持单元状态的记忆，通过遗忘门f_t和输入门i_t进行调节；遗忘门的作用是让细胞记住或忘记之前的状态C_t-1，输入门的作用是允许或阻止传入信号更新单元状态；输出门的作用是控制单元状态C输出和传输到下一个单元；将归一化后的训练集样本数据依次经过遗忘门、输入门、输出门处理后输出当前数据的预测数据。

9.根据权利要求8所述的一种基于IPSO-LSTM的风电功率预测方法，其特征在于，所述LSTM网络模型是由多个单元结构组成，该结构能通过更新内部状态来长时间存储信息，函数如下：

f_t＝σ(w_f·[h_t-1,x_t]+b_f)

i_t＝σ(w_i·[h_t-1,x_t]+b_i)

o_t＝σ(w_o·[h_t-1,x_t]+b_o)

h_t＝o_t·tanh(C_t)

式中，x为LSTM单元的输入向量；h为LSTM单元的输出向量，f、i、o分别表示遗忘门、输入门和输出门；C表示单元状态；下标t表示时刻；σ、tanh分别为sigmoid、tanh激活函数；w和b分别表示权重和偏差矩阵。

Images