基于数据驱动集合卡尔曼滤波的电力负荷预测方法和装置

技术领域

本发明涉及电力系统技术领域，尤其涉及一种基于数据驱动集合卡尔曼滤波的电力负荷预测方法和装置。

背景技术

近年来随着智能电网的快速发展和基于非侵入式终端的用电信息采集系统渐渐完善，电力大数据逐渐成为人们研究的热点。在大数据背景下分析和挖掘电力负荷数据之间的联系，运用合理的方法探索电力用户负荷的未来发展变化规律，对于提升电网运行的稳定性和经济性具有重大的现实意义。

短期负荷由于受天气变化、社会活动等因素的影响，在时间序列上表现为确定性系统的不规则状态，即系统的混沌性质。但是各影响因素中大部分都具有潜在的规律性，从而为电力负荷的有效预测提供了基础条件。目前用于短期负荷预测的算法主要有神经网络法、时间序列法、回归分析法、支持向量机法、模糊预测法等。负荷预测的核心问题是如何利用历史数据建立预测模型，对未来某一时刻的负荷进行预测。随着非侵入式终端的推广应用，以及天气预测水平的不断提高，准确获取海量历史数据已经易于实现。因此，短期负荷预测的核心问题是预测方法的精度高低。另外，基于非侵入式终端采样得到的海量数据中势必含有各种噪声，这些噪声的存在会对电力负荷的预测精度产生恶劣影响。因此，如何应用先进的噪声滤除技术来实现精准的数据同化以实现高精度的电力负荷预测，这也是一个值得深入研究的方向。

发明内容

本发明的目的在于：公开了一种基于数据驱动集合卡尔曼滤波的电力负荷预测方法和装置，解决了负荷预测精度不高的问题。

为实现上述目的，本发明的技术方案如下：一种基于数据驱动集合卡尔曼滤波的电力负荷预测方法，包括步骤：

获取电力用户负荷数据，作为电力负荷训练数据集；

基于电力负荷训练数据集进行多维相空间重构，得到电力负荷的多维向量；

基于电力负荷的多维向量对电力负荷进行预测，选取预测值点集得到滤波器数据集合，通过集合卡尔曼滤波对滤波器数据集合中的数据元素进行校正，得到电力负荷预测值。

进一步的，多维相空间重构的方法为：基于电力负荷训练数据集和延迟坐标方法进行多维相空间重构。

进一步的，所述延迟坐标方法中的延迟时间确定方法为：

电力负荷训练数据集中的电力负荷时间序列为{x(i)}，i＝1,2,…,n，存在其自相关函数

进一步的，所述延迟坐标方法中的嵌入维数确定方法为：

使用虚假最临近点法进行选取，应用电力负荷训练数据集中的电力负荷时间序列，从嵌入维数2开始计算虚假最临近点的比例，逐步增加维数，直至虚假最临近点的比例小于设定比例或者其数目不再随着维数的增加而减少，将此时的维数选作嵌入维数。

进一步的，基于电力负荷的多维向量对电力负荷进行预测，方法为：

选取k个与电力负荷多维向量中最后一个向量最临近的向量，再取这k个最临近向量的平均值作为下个延迟时间的多维向量，实现下个延迟时间的电力负荷状态的预测。

进一步的，使用循环计算的方法来确定最临近向量的数目，从k＝1开始，直至电力负荷状态预测误差小于设定值以及计算时间小于设定时间为止。

进一步的，选取预测值点集得到滤波器数据集合，方法为：

采取对称采样的策略选取2m+1个Sigma点及其权重来表征电力负荷预测的状态，预测均值为Sigma点集的中心，其它所有点关于均值对称分布，且所有权重之和为1，m为电力负荷状态变量的维数；Sigma点组成预测值点集。

一种基于数据驱动集合卡尔曼滤波的电力负荷预测装置，包括：

电力负荷训练数据集获取模块，用于获取电力用户负荷数据，作为电力负荷训练数据集；

多维向量获取模块，用于基于电力负荷训练数据集进行多维相空间重构，得到电力负荷的多维向量；

预测模块，用于基于电力负荷的多维向量对电力负荷进行预测，选取预测值点集得到滤波器数据集合，通过集合卡尔曼滤波对滤波器数据集合中的数据元素进行校正，得到电力负荷预测值。

进一步的，多维相空间重构的方法为：基于电力负荷训练数据集和延迟坐标方法进行多维相空间重构；

所述延迟坐标方法中的延迟时间确定方法为：

电力负荷训练数据集中的电力负荷时间序列为{x(i)}，i＝1,2,…,n，存在其自相关函数

所述延迟坐标方法中的嵌入维数确定方法为：

使用虚假最临近点法进行选取，应用电力负荷训练数据集中的电力负荷时间序列，从嵌入维数2开始计算虚假最临近点的比例，逐步增加维数，直至虚假最临近点的比例小于设定比例或者其数目不再随着维数的增加而减少，将此时的维数选作嵌入维数。

进一步的，基于电力负荷的多维向量对电力负荷进行预测，方法为：

选取k个与电力负荷多维向量中最后一个向量最临近的向量，再取这k个最临近向量的平均值作为下个延迟时间的多维向量，实现下个延迟时间的电力负荷状态的预测；

使用循环计算的方法来确定最临近向量的数目，从k＝1开始，直至电力负荷状态预测误差小于设定值以及计算时间小于设定时间为止；

选取预测值点集得到滤波器数据集合，方法为：

采取对称采样的策略选取2m+1个Sigma点及其权重来表征电力负荷预测的状态，预测均值为Sigma点集的中心，其它所有点关于均值对称分布，且所有权重之和为1，m为电力负荷状态变量的维数；Sigma点组成预测值点集。

有益效果：本发明基于电力负荷的多维向量对电力负荷进行预测，选取预测值点集得到滤波器数据集合，通过集合卡尔曼滤波对滤波器数据集合中的数据元素进行校正，得到准确的电力负荷预测值，通过结合集成集合卡尔曼滤波和相空间重构方法能得到更好的负荷预测效果。

附图说明

图1是本发明具体实施方式中的一种负荷预测方法流程图；

图2是电热水器负荷预测值、实际值及误差图；

图3是空调负荷预测值、实际值及误差图；

图4是总负荷预测值与实际的比较图；

图5是总负荷预测误差图；

图6是本发明方法与基于传统相空间重构的负荷预测误差图。

具体实施方式

下面结合附图对本发明作进一步详细说明。

本发明所述的数据驱动，是指针对没有模型的随机系统，使用大量的采集数据对系统行为进行预测。

实施例1：

如图1所示，一种基于数据驱动集合卡尔曼滤波的电力负荷预测方法，包括步骤：

步骤1，获取非侵入式终端采集的电力用户负荷数据，作为电力负荷训练数据集；

电力负荷训练数据集，包括各种采样时间间隔的电力用户负荷数据。

非侵入式终端，可以是单相终端，也可以是三相终端。

步骤2，基于电力负荷训练数据集和延迟坐标方法进行多维相空间重构，得到电力负荷的多维向量；

延迟坐标方法是通过一维时间序列的不同延迟时间来构建多维的相空间矢量。

在相空间重构过程中，首先需要选择延时坐标方法中用到的两个参数，即延迟时间和嵌入维数。嵌入维数是指原来是一维的时间序列，通过延迟坐标的方法处理后，变成多维时间序列，最终的维数值就是嵌入维数。

使用自相关函数法来确定延迟时间。假设电力用户负荷数据{x(i)}，i＝1,2,…,n，存在其自相关函数

x(i)代表i时刻电力负荷的值，i的范围为1,2,…,n，n为电力负荷时间序列的最大长度；

下降系数可以根据具体应用进行选择，下降系数一个常用的经验值为1-1/e。

关于嵌入维数，使用虚假最临近点法进行选取。对于多维相空间中的每一个向量，都有一个欧几里德距离的最邻近点。当相空间的维数增加时，两点之间的距离会发生变化。如果距离变化超过判断门槛值，就称这两点为虚假临近点。应用电力负荷训练数据集中的电力负荷时间序列，从嵌入维数2开始计算虚假最临近点的比例，然后逐步增加维数，直至虚假最临近点的比例小于设定比例或者其数目不再随着维数的增加而减少，即可将此时的维数选作嵌入维数。

根据嵌入维数和延迟时间，得到多维向量，构成了多维相空间，从而完成了相空间的重构。

延迟时间、嵌入维数选取结束后，基于电力负荷训练数据集来重构相空间和使用该相空间进行各种时间尺度的电力负荷预测。

虚假最临近点的判断门槛值可以根据具体情况灵活选择，判断门槛值通常大于10且小于50。

步骤3，基于电力负荷的多维向量对电力负荷进行预测，选取预测值点集得到滤波器数据集合，通过集合卡尔曼滤波对滤波器数据集合中的数据元素进行校正，得到准确的电力负荷预测值。

根据连续函数的性质，如果两个向量

使用循环计算的方法来确定最临近向量的数目，从k＝1开始，直至电力负荷状态预测误差小于设定值以及计算时间小于设定时间为止。

通过无迹变换的方法选择预测值点集，直至滤波器数据集合中的数目达到预先设定值q。将滤波器数据集合的q个数据输入集合卡尔曼滤波器，对电力负荷状态预测值进行优化滤波，从而得到一个更高精度的电力负荷预测值。

通过无迹变换的方法选择预测值点集的方法为：采取对称采样的策略选取2m+1(m为电力负荷状态变量的维数)个Sigma点及其相关权重来表征电力负荷预测的状态，预测均值为Sigma点集的中心，其它所有点关于均值对称分布，且所有权重之和为1。

可以根据Sigma点组成的滤波器预测状态集合，依次产生相应的状态误差集合、观测误差集合、预测状态误差协方差、预测状态误差与观测误差之间的协方差、观测误差协方差。在产生预测状态集合的基础上，使用集合卡尔曼滤波方程对预测集合中的数据元素进行校正。

实施例2：

一种基于数据驱动集合卡尔曼滤波的电力负荷预测装置，包括：

电力负荷训练数据集获取模块，用于获取电力用户负荷数据，作为电力负荷训练数据集；

多维向量获取模块，用于基于电力负荷训练数据集进行多维相空间重构，得到电力负荷的多维向量；

预测模块，用于基于电力负荷的多维向量对电力负荷进行预测，选取预测值点集得到滤波器数据集合，通过集合卡尔曼滤波对滤波器数据集合中的数据元素进行校正，得到电力负荷预测值。

进一步的，多维相空间重构的方法为：基于电力负荷训练数据集和延迟坐标方法进行多维相空间重构；

所述延迟坐标方法中的延迟时间确定方法为：

电力负荷训练数据集中的电力负荷时间序列为{x(i)}，i＝1,2,…,n，存在其自相关函数

所述延迟坐标方法中的嵌入维数确定方法为：

使用虚假最临近点法进行选取，应用电力负荷训练数据集中的电力负荷时间序列，从嵌入维数2开始计算虚假最临近点的比例，逐步增加维数，直至虚假最临近点的比例小于设定比例或者其数目不再随着维数的增加而减少，将此时的维数选作嵌入维数。

进一步的，基于电力负荷的多维向量对电力负荷进行预测，方法为：

选取k个与电力负荷多维向量中最后一个向量最临近的向量，再取这k个最临近向量的平均值作为下个延迟时间的多维向量，实现下个延迟时间的电力负荷状态的预测；

使用循环计算的方法来确定最临近向量的数目，从k＝1开始，直至电力负荷状态预测误差小于设定值以及计算时间小于设定时间为止；

选取预测值点集得到滤波器数据集合，方法为：

采取对称采样的策略选取2m+1个Sigma点及其权重来表征电力负荷预测的状态，预测均值为Sigma点集的中心，其它所有点关于均值对称分布，且所有权重之和为1，m为电力负荷状态变量的维数；Sigma点组成预测值点集。

案例：

在电力负荷的预测过程中，由于数学模型的简化、各种影响因子的忽略、不当的参数选取、训练样本不完整等因素的存在，导致预测值与真实值之间一般总有误差。一种简单的评价指标是绝对误差，即负荷预测值与真实值之间的差值。但是，这种评价指标对误差大小的表达不够直观。本发明使用相对误差来作为评价指标，更加直观。假设y(t)和

训练样本数据：

训练样本数据来自非侵入式终端量测得到的河北省保定市五个家庭用户2019年7月10日到7月15日的用电数据，包括电热水器、空调的负荷以及家庭总负荷，采样间隔为15分钟。表1给出了所采集的训练样本数据。

表1训练样本集

对于一维电力负荷时间序列，对每一个状态选取3个Sigma点组成点集。中心点的权重W

表2主要参数

实验结果：

通过本发明方法分别得到热水器、空调和家庭总负荷的预测值。图2给出了一天中电热水器负荷的预测值、实际值和相对误差。可以看出，误差范围主要集中在+7～-10之间，预测效果令人满意。图3给出了一天中空调负荷的预测值、实际值和相对误差，误差范围主要集中在+15～-15之间，预测效果比较令人满意。由于空调的使用受个人习惯、家庭经济状况等众多因素影响，因此随机性较大，预测误差相对较大。图4为家庭总负荷的预测值、实际值的比较结果(纵轴数值应放大10倍)，相对误差如图5所示，误差范围主要集中在+9％～-8％之间，取得了较好的预测效果。同时，将单纯基于相空间重构方法的预测结果与实际值作比较，相对误差如图6所示，误差范围为+13％～-20％。对比图5和图6发现集成集合卡尔曼滤波的相空间重构方法能得到更好的负荷预测效果。

本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。