CN114021847A - 基于isgu混合模型的短期电力负荷预测方法及系统 - Google Patents

Abstract

本发明涉及一种基于ISGU混合模型的短期电力负荷预测方法及系统，该方法包括以下步骤：步骤S1：输入电力负荷数据及相关因素，并进行数据预处理，包括对星期、节假日因素进行数字化处理以及对电力负荷特征值进行筛选，然后划分得到训练集和预测集；步骤S2：构建结合GRU网络和ISSA算法的ISGU混合模型，并通过训练集和预测集对ISGU混合模型进行训练和测试，得到训练好的ISGU混合模型；步骤S3：基于训练好的ISGU混合模型进行短期电力负荷预测，获得电力负荷预测结果。该方法及系统有利于提高短期电力负荷预测精度，进而提高电力系统的稳定性和安全性。

Description

基于ISGU混合模型的短期电力负荷预测方法及系统

技术领域

本发明属于电力负荷预测技术领域，具体涉及一种基于ISGU混合模型的短期电力负荷预测方法及系统。

背景技术

随着我国社会经济的迅猛发展，市场逐渐对电能供应的安全性、可靠性、经济性等方面要求提高。因此，电力公司在运行中面临着诸多经济技术挑战，在这些挑战中，电力系统的调度、潮流分析、规划和控制最为突出。面对电能需求随时间不断变化的特点，电力部门应尽可能保证电能供需处于动态平衡状态，以此追求最优的发电设备利用率与调度经济性，避免产能过剩或供应不足。精确的电力负荷预测不仅有助于规划供电设备的运行检修计划，保证供电可靠性，并且还能为新机组投入计划与电网增容改建计划给予一定参考。现如今，多样化预测模型、智能化电网改造给负荷预测提供技术以及数据支持。

近年来，多种预测方法在短期电力负荷预测中得到广泛应用与发展。1、利用时间序列法进行短期负荷的预测，其优势是模型简单，计算迅速，需要的数据容量不大，对较为平稳的时间序列进行预测，其效果较为良好。然而，电力负荷存在波动性，电力负荷预测存在不确定性、条件性，故使得时间序列预测精确较低；2、利用BP神经网络进行短期负荷的预测，然而反向传播(BP)神经网络利用反向传播的算法来反复更新神经网络中的权重与阈值的方式，使其容易陷入局部最优；3、利用极限学习机(ELM)算法进行短期负荷的预测，ELM算法的单隐含层的结构以及利用最小二乘法更新输出层权重的特点，使得其学习速度快、泛化性能较好，然而输入层初始化权重与阈值的随机生成仍会影响预测结果；4、利用支持向量机(SVM)算法进行短期负荷的预测，SVM鲁棒性好、泛化能力强，但其模型训练耗时较长，面对大规模样本数据较为乏力。5、利用循环神经网络(RNN)进行短期负荷的预测，RNN擅于挖掘时间序列数据的特点，但其存在梯度爆炸与梯度消散等问题。6、利用引入参数寻优的混合模型，例如GA-KELM、PSO-SVM、CS-SVM等等。参数寻优算法的引入可以客观性的配置预测模型算法的参数，使得预测模型得精度提高。但某些参数寻优算法存在初始种群多样性不足、容易陷入局部最优等缺点。

发明内容

本发明的目的在于提供一种基于ISGU混合模型的短期电力负荷预测方法及系统，该方法及系统有利于提高短期电力负荷预测精度，进而提高电力系统的稳定性和安全性。

为实现上述目的，本发明采用的技术方案是：一种基于ISGU混合模型的短期电力负荷预测方法，包括以下步骤：

步骤S1：输入电力负荷数据及相关因素，并进行数据预处理，包括对星期、节假日因素进行数字化处理以及对电力负荷特征值进行筛选，然后划分得到训练集和预测集；

步骤S2：构建结合GRU网络和ISSA算法的ISGU混合模型，并通过训练集和预测集对ISGU混合模型进行训练和测试，得到训练好的ISGU混合模型；

步骤S3：基于训练好的ISGU混合模型进行短期电力负荷预测，获得电力负荷预测结果。

进一步地，步骤S1中，对星期、节假日因素进行数字化处理的具体方法为：

设各星期类型平均负荷所组成的矩阵为W_k＝[w_mon,w_tue,w_wed,w_thu,w_fri,w_sat,w_sun]；则星期因素映射方式如公式(1)所示：

其中W_map为星期因素映射值，缩放倍数b∈(0,∞),d∈(0,∞)为常数；

节假日映射与星期映射同理，即利用节假日平均负荷值与其对应的星期类型的平均负荷值进行处理；当两者映射值分别求出后，将两者映射值整合为一体，如公式(2)所示：

WH_map＝W_map+H_map (2)

其中WH_map为星期-节假日映射值，H_map为节假日映射值。

进一步地，步骤S1中，采用PCCs与MI算法对电力负荷特征值进行筛选，PCCs用于筛选气候、星期、节假日等特征值，MI用于筛选历史电力负荷特征值，筛选出的特征值记为C_v，其原理如公式(3)(4)所示：

其中，x为电力负荷信号，F为相关因素，I为样本总个数，i为第i个样本，F_i为第i个样本的相关因素，

为所有样本相关因素平均值，x为第i个样本的电力负荷信号值，

为所有样本电力负荷信号平均值，J为样本总个数，j为第j个样本，这里I＝J；p(*,·)为联合分布函数，q(·)为边缘分布函数。

进一步地，步骤S2中，构建GRU网络进行电力负荷预测，其输入参数为筛选出的电力负荷特征值C_v；同时将ISSA算法引入GRU网络的构建，利用ISSA算法对GRU网络中的隐含层节点数与迭代次数进行寻优，以得到结合GRU网络和ISSA算法的ISGU混合模型。

进一步地，GRU网络由重置门、更新门与输出门组成；GRU的重置门与更新门分别确定前一时刻状态记忆变量h_λ-1的遗忘程度与保留程度，其数学模型分别如公式(5)(6)所示：

其中C_v为输入向量，σ(·)为sigmoid激活函数，W_r为重置门权重，r为重置门状态，W_u为更新门权重，z为更新门状态；

根据重置门状态r_λ与C_v、h_λ-1得出当前时刻候选状态记忆变量

其数学模型如公式(7)所示：

其中tanh(·)为双曲正切激活函数，W_o为计算

部分的权重；

最后根据更新门状态z_λ对h_λ-1、

进行遗忘与保留，由此获得GRU网络输出，其公式如(8)所示：

将h_λ输入至全连接层中得到负荷预测结果Y^pred。

进一步地，ISSA算法的实现方法为：

首先采用对立学习算法生成初始化种群，生成种群规模为M的初始化随机种群，记为P’_M×N，则对立种群的生成如公式(9)所示：

其中lb_n表示第n个待寻优参数的下限，ub_n表示第n个待寻优参数的上限；

根据2M个种群个体所对应适应度的排序，取前M个适应度最优的种群个体作为最终的初始化种群，记为P_M×N；

发现者为种群中适应度较好的个体，其进行食物搜索并指导其余个体的运动；发现者位置更新公式如(10)所示：

其中R₂∈[0,1]和ST∈[0.5,1]分别代表警报值和安全阈值，rand与Q是[0,1]之间的随机数；L表示1×N的矩阵，其元素均为1；

采用鸟群算法对发现者位置更新思路进行替代；

在麻雀种群的加入者位置更新中引入随机游走策略，其公式如(11)～(13)所示：

A⁺＝A^T(AA^T)^-1 (13)

其中P_pd为发现者当前最佳位置；A为1×N矩阵，其元素随机赋值1或-1；a(t)为控制因子，c(t)为满足柯西分布的随机数；t为当前迭代次数，t_max为最大迭代次数；

为避免适应度低的加入者往搜索边界进行更新，引入随机游走位置更新策略；

麻雀群中的警戒者是随机产生，其位置更新公式如(14)所示：

其中P_best为当前全局最优位置，P_worst为当前全局最差位置，β为[0,1]之间随机数，K为为[-1,1]之间随机数，f_m为当前麻雀个体适应度，f_g与f_w为当前全局最优与最差适应度；ε为趋近于0的数，其用于避免分母为0；

最后采用贪心算法决定各麻雀个体位置更新结果；

ISSA中的随机数rand、Q、β、K均采用PWMLCM混沌映射产生，其公式如(15)所示：

其中x_r∈[0,1]，控制参数s∈[0,0.5]，x^l _r为第l代随机数。

本发明还提供了一种基于ISGU混合模型的短期电力负荷预测系统，包括存储器、处理器以及存储于存储器上并能够被处理器运行的计算机程序指令，当处理器运行该计算机程序指令时，能够实现上述的方法步骤。

与现有技术相比，本发明具有以下有益效果：提供了一种基于ISGU混合模型的短期电力负荷预测方法及系统，该方法及系统可以克服负荷波动性与非线性、电力负荷预测的条件性以及算法参数设置无客观依据等问题给电力负荷预测所带来的困难，从而提高短期电力负荷预测精度，保证电能供需处于动态平衡，减少资源的浪费与产能成本，提高电力系统的稳定性和安全性。

附图说明

图1为本发明实施例的方法实现流程图；

图2为本发明实施例中考虑星期-节假日因素的节假日预测结果图；

图3为本发明实施例中GRU结构图；

图4为本发明实施例中ISSA参数寻优过程图。

具体实施方式

下面结合附图及实施例对本发明做进一步说明。

应该指出，以下详细说明都是示例性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。

需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。

如图1所示，本实施例提供了一种基于ISGU混合模型的短期电力负荷预测方法，包括以下步骤：

步骤S1：输入电力负荷数据及相关因素，并进行数据预处理，包括对星期、节假日因素进行数字化处理以及对电力负荷特征值进行筛选，然后划分得到训练集和预测集。

步骤S2：构建结合GRU网络和ISSA算法的ISGU混合模型，并通过训练集和预测集对ISGU混合模型进行训练和测试，得到训练好的ISGU混合模型。

步骤S3：基于训练好的ISGU混合模型进行短期电力负荷预测，获得电力负荷预测结果。

1、对星期、节假日因素进行数字化处理

星期与节假日因素与电力负荷存在一定相关性，然而星期与节假日为非数字化因素，因此本发明提出一种基于日平均负荷值的星期-节假日因素映射方法，其具体方法为：

设各星期类型平均负荷所组成的矩阵为W_k＝[w_mon,w_tue,w_wed,w_thu,w_fri,w_sat,w_sun]；则星期因素映射方式如公式(1)所示：

其中W_map为星期因素映射值，缩放倍数b∈(0,∞),d∈(0,∞)为常数。

节假日映射与星期映射同理，即利用节假日平均负荷值与其对应的星期类型的平均负荷值进行处理；当两者映射值分别求出后，将两者映射值整合为一体，如公式(2)所示：

WH_map＝W_map+H_map (2)

其中WH_map为星期-节假日映射值，H_map为节假日映射值。

引入考虑星期-节假日因素，并采用基于日平均负荷值的星期-节假日因素映射方法后进行负荷预测的结果如图2所示。

2、特征值筛选

由于短期电力负荷预测的条件性特点，使影响电力负荷的相关因素众多。因此，本发明采用PCCs与MI算法对电力负荷特征值进行筛选，PCCs用于筛选气候、星期、节假日等特征值，MI用于筛选历史电力负荷特征值，筛选出的特征值记为C_v，其原理如公式(3)(4)所示：

其中，x为电力负荷信号，F为相关因素，I为样本总个数，i为第i个样本，F_i为第i个样本的相关因素，

为所有样本相关因素平均值，x为第i个样本的电力负荷信号值，

为所有样本电力负荷信号平均值，J为样本总个数，j为第j个样本，这里I＝J；p(*,·)为联合分布函数，q(·)为边缘分布函数。

3、基于ISGU混合模型进行短期电力负荷预测

构建GRU网络进行电力负荷预测，其输入参数为筛选出的电力负荷特征值C_v；同时将ISSA算法引入GRU网络的构建，利用ISSA算法对GRU网络中的隐含层节点数与迭代次数进行寻优，以得到结合GRU网络和ISSA算法的ISGU混合模型。

GRU网络由重置门、更新门与输出门组成，其结构如图3所示。GRU对LSTM网络的简化改进，不仅减少网络所需要训练的参数，缩短模型训练时间，并且也解决RNN的梯度消散问题。GRU的重置门与更新门分别确定前一时刻状态记忆变量h_λ-1的遗忘程度与保留程度，其数学模型分别如公式(5)(6)所示：

其中C_v为输入向量，σ(·)为sigmoid激活函数，W_r为重置门权重，r为重置门状态，W_u为更新门权重，z为更新门状态。

根据重置门状态r_λ与C_v、h_λ-1得出当前时刻候选状态记忆变量

其数学模型如公式(7)所示：

其中tanh(·)为双曲正切激活函数，W_o为计算

部分的权重。

最后根据更新门状态z_λ对h_λ-1、

进行遗忘与保留，由此获得GRU网络输出，其公式如(8)所示：

将h_λ输入至全连接层中得到负荷预测结果Y^pred。

GRU网络参数设置对预测精度影响较大，然而其网络参数设置未有客观标准，因此本发明提出ISSA算法对GRU的隐含层节点数与迭代次数进行寻优。通过ISSA在对预测模型参数寻优的过程如图4所示。

Jiankai Xue等人在2020年根据麻雀觅食和抗掠夺行为的启发提出麻雀算法。但SSA的初始种群缺少多样性，个别位置更新公式对搜索产生一定限制，导致SSA算法并不一定能搜索到较优秀结果。因此，本发明对SSA的初始种群生成与位置更新公式进行改进，ISSA算法的实现方法如下所示。

首先采用对立学习算法生成初始化种群，以此克服SSA算法初始化种群缺少多样性问题。现生成种群规模为M的初始化随机种群，记为P’_M×N，则对立种群的生成如公式(9)所示：

其中lb_n表示第n个待寻优参数的下限，ub_n表示第n个待寻优参数的上限。

根据2M个种群个体所对应适应度的排序，取前M个适应度最优的种群个体作为最终的初始化种群，记为P_M×N。

发现者为种群中适应度较好的个体，其进行食物搜索并指导其余个体的运动；发现者位置更新公式如(10)所示：

其中R₂∈[0,1]和ST∈[0.5,1]分别代表警报值和安全阈值，rand与Q是[0,1]之间的随机数；L表示1×N的矩阵，其元素均为1。

SSA中发现者在R₂<ST时的位置更新公式使得发现者每个维度都在缩小，则发现者搜索食物范围受到限制。因此，本发明采用鸟群算法对发现者位置更新思路进行替代。

在麻雀种群的加入者位置更新中引入随机游走策略，其公式如(11)～(13)所示：

A⁺＝A^T(AA^T)^-1 (13)

其中P_pd为发现者当前最佳位置；A为1×N矩阵，其元素随机赋值1或-1；a(t)为控制因子，c(t)为满足柯西分布的随机数；t为当前迭代次数，t_max为最大迭代次数。

适应度低的加入者为获得更多能量需飞往其他地方进行觅食，而非跟随发现者。在SSA中该类加入者的位置更新方式容易使其往搜索边界进行更新。为避免适应度低的加入者往搜索边界进行更新，本发明引入随机游走位置更新策略。

麻雀群中的警戒者是随机产生，其位置更新公式如(14)所示：

其中P_best为当前全局最优位置，P_worst为当前全局最差位置，β为[0,1]之间随机数，K为为[-1,1]之间随机数，f_m为当前麻雀个体适应度，f_g与f_w为当前全局最优与最差适应度；ε为趋近于0的数，其用于避免分母为0。

最后采用贪心算法决定各麻雀个体位置更新结果。

ISSA中的随机数rand、Q、β、K均采用PWMLCM混沌映射产生，其公式如(15)所示：

其中x_r∈[0,1]，控制参数s∈[0,0.5]，x^l _r为第l代随机数。混沌随机数的产生：先随机生成初始随机数，然后根据公式(15)与第一代随机数生成第二代，由此类推。

本实施例还提供了一种基于ISGU混合模型的短期电力负荷预测系统，包括存储器、处理器以及存储于存储器上并能够被处理器运行的计算机程序指令，当处理器运行该计算机程序指令时，能够实现上述的方法步骤。

以上所述，仅是本发明的较佳实施例而已，并非是对本发明作其它形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型，仍属于本发明技术方案的保护范围。

Claims (7)

1.一种基于ISGU混合模型的短期电力负荷预测方法，其特征在于，包括以下步骤：

步骤S1：输入电力负荷数据及相关因素，并进行数据预处理，包括对星期、节假日因素进行数字化处理以及对电力负荷特征值进行筛选，然后划分得到训练集和预测集；

步骤S2：构建结合GRU网络和ISSA算法的ISGU混合模型，并通过训练集和预测集对ISGU混合模型进行训练和测试，得到训练好的ISGU混合模型；

步骤S3：基于训练好的ISGU混合模型进行短期电力负荷预测，获得电力负荷预测结果。

2.根据权利要求1所述的基于ISGU混合模型的短期电力负荷预测方法，其特征在于，步骤S1中，对星期、节假日因素进行数字化处理的具体方法为：

设各星期类型平均负荷所组成的矩阵为W_k＝[w_mon,w_tue,w_wed,w_thu,w_fri,w_sat,w_sun]；则星期因素映射方式如公式(1)所示：

其中W_map为星期因素映射值，缩放倍数b∈(0,∞),d∈(0,∞)为常数；

节假日映射与星期映射同理，即利用节假日平均负荷值与其对应的星期类型的平均负荷值进行处理；当两者映射值分别求出后，将两者映射值整合为一体，如公式(2)所示：

WH_map＝W_map+H_map (2)

其中WH_map为星期-节假日映射值，H_map为节假日映射值。

3.根据权利要求1所述的基于ISGU混合模型的短期电力负荷预测方法，其特征在于，步骤S1中，采用PCCs与MI算法对电力负荷特征值进行筛选，PCCs用于筛选气候、星期、节假日等特征值，MI用于筛选历史电力负荷特征值，筛选出的特征值记为C_v，其原理如公式(3)(4)所示：

其中，x为电力负荷信号，F为相关因素，I为样本总个数，i为第i个样本，F_i为第i个样本的相关因素，

为所有样本相关因素平均值，x为第i个样本的电力负荷信号值，

为所有样本电力负荷信号平均值，J为样本总个数，j为第j个样本，这里I＝J；p(*,·)为联合分布函数，q(·)为边缘分布函数。

4.根据权利要求1所述的基于ISGU混合模型的短期电力负荷预测方法，其特征在于，步骤S2中，构建GRU网络进行电力负荷预测，其输入参数为筛选出的电力负荷特征值C_v；同时将ISSA算法引入GRU网络的构建，利用ISSA算法对GRU网络中的隐含层节点数与迭代次数进行寻优，以得到结合GRU网络和ISSA算法的ISGU混合模型。

5.根据权利要求4所述的基于ISGU混合模型的短期电力负荷预测方法，其特征在于，GRU网络由重置门、更新门与输出门组成；GRU的重置门与更新门分别确定前一时刻状态记忆变量h_λ-1的遗忘程度与保留程度，其数学模型分别如公式(5)(6)所示：

其中C_v为输入向量，σ(·)为sigmoid激活函数，W_r为重置门权重，r为重置门状态，W_u为更新门权重，z为更新门状态；

根据重置门状态r_λ与C_v、h_λ-1得出当前时刻候选状态记忆变量

其数学模型如公式(7)所示：

其中tanh(·)为双曲正切激活函数，W_o为计算

部分的权重；

最后根据更新门状态z_λ对h_λ-1、

进行遗忘与保留，由此获得GRU网络输出，其公式如(8)所示：

将h_λ输入至全连接层中得到负荷预测结果Y^pred。

6.根据权利要求4所述的基于ISGU混合模型的短期电力负荷预测方法，其特征在于，ISSA算法的实现方法为：

首先采用对立学习算法生成初始化种群，生成种群规模为M的初始化随机种群，记为P’_M×N，则对立种群的生成如公式(9)所示：

其中lb_n表示第n个待寻优参数的下限，ub_n表示第n个待寻优参数的上限；

根据2M个种群个体所对应适应度的排序，取前M个适应度最优的种群个体作为最终的初始化种群，记为P_M×N；

发现者为种群中适应度较好的个体，其进行食物搜索并指导其余个体的运动；发现者位置更新公式如(10)所示：

其中R₂∈[0,1]和ST∈[0.5,1]分别代表警报值和安全阈值，rand与Q是[0,1]之间的随机数；L表示1×N的矩阵，其元素均为1；

采用鸟群算法对发现者位置更新思路进行替代；

在麻雀种群的加入者位置更新中引入随机游走策略，其公式如(11)～(13)所示：

A⁺＝A^T(AA^T)^-1 (13)

其中P_pd为发现者当前最佳位置；A为1×N矩阵，其元素随机赋值1或-1；a(t)为控制因子，c(t)为满足柯西分布的随机数；t为当前迭代次数，t_max为最大迭代次数；

为避免适应度低的加入者往搜索边界进行更新，引入随机游走位置更新策略；

麻雀群中的警戒者是随机产生，其位置更新公式如(14)所示：

其中P_best为当前全局最优位置，P_worst为当前全局最差位置，β为[0,1]之间随机数，K为为[-1,1]之间随机数，f_m为当前麻雀个体适应度，f_g与f_w为当前全局最优与最差适应度；ε为趋近于0的数，其用于避免分母为0；

最后采用贪心算法决定各麻雀个体位置更新结果；

ISSA中的随机数rand、Q、β、K均采用PWMLCM混沌映射产生，其公式如(15)所示：

其中x_r∈[0,1]，控制参数s∈[0,0.5]，x^l _r为第l代随机数。

7.一种基于ISGU混合模型的短期电力负荷预测系统，其特征在于，包括存储器、处理器以及存储于存储器上并能够被处理器运行的计算机程序指令，当处理器运行该计算机程序指令时，能够实现如权利要求1-6所述的方法步骤。

Images