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

1.本发明涉及电力系统技术领域，尤其是涉及一种基于混合模型的短期电力负荷预测方法及系统。


背景技术：

2.短期负荷预测(stlf)是预测一小时至一周的电力负荷，对电力系统的规划、管理和稳定运行有着重要作用。研究表明，10gw电力公司的预测误差减少1％，每年可节省160万美元。因此，准确的短期负荷预测对电力公司减少电力浪费，保持供需平衡，降低生产成本，提高经济收益，以及调度管理和用电需求规划具有重要作用。
3.多年来许多国内外学者在负荷预测领域里不断探索研究，主要采用如下几种方法：统计学方法、机器学习方法和混合方法。由于电力负荷数据的主要特点是高度非线性和非平稳性，受天气、电价、节假日等诸多因素影响，预测难度也会随之加大。
4.其中，统计学方法内部结构简单，收敛速度快，能够很好的处理线性预测问题，但处理非线性问题精度不足；支持向量机(svm)、极限学习机(elm)的机器学习方法能够很好拟合非线性负荷数据，但不能提取时序数据中的时态特征；长短期记忆神经网络(lstm)的机器学习方法能够很好处理时序数据，但存在单一模型预测精度局限性问题；混合方法则在负荷预测方面相对于单一方法有一定的优越性，但现有的混合方法对电力负荷数据经过分解后的子序列使用相同的模型，而不同的子序列往往拥有不同的特征信息，不利于特征信息提取和提升负荷预测性能。
5.故而亟需提出一种短期电力负荷预测方法来解决所提出的问题。


技术实现要素：

6.基于此，有必要针对现有技术的不足，提供一种基于混合模型的短期电力负荷预测方法及系统，针对不同子序列设计出不同输入特征和模型，可以极大的丰富输入特征信息维度，降低模型复杂度，提高预测精度。
7.为解决上述技术问题，本发明所采用的技术方案是：
8.第一方面，提供了一种基于混合模型的短期电力负荷预测方法，其包括如下步骤：
9.采集电力负荷数据，组成电力负荷数据集，并对电力负荷数据进行归一化处理，获取归一化处理后的电力负荷数据，记为原始电力负荷数据；
10.通过ceemdan分解算法对原始电力负荷数据进行分解处理，得到多个不同频率的固有模态函数分量和一个剩余残差分量；
11.将多个不同频率的固有模态函数分量分类为高频分量子序列及低频分量子序列；
12.通过lstm预测模型提取高频分量子序列的时序特征；
13.将lstm预测模型提取的高频分量子序列的时序特征序列、原始电力负荷数据中预测点前k个历史电力负荷数据、高频分量子序列中预测点前k个历史电力负荷数据输入到elm-catboost混合预测模型中，获得高频分量子序列的预测结果；其中，所述elm-catboost
混合预测模型包括第一elm预测模型及catboost预测模型；
14.通过第二elm预测模型获取低频分量子序列的预测结果；
15.将高频分量子序列的预测结果和低频分量子序列的预测结果进行叠加，得到叠加结果；
16.将所述叠加结果进行反归一化处理后，得到基于混合模型的短期电力负荷预测结果。
17.第二方面，提供了一种基于混合模型的短期电力负荷预测系统，其包括
18.采集模块，用于采集电力负荷数据，组成电力负荷数据集，并对电力负荷数据进行归一化处理，获取归一化处理后的电力负荷数据，记为原始电力负荷数据；
19.分解模块，用于通过ceemdan分解算法对原始电力负荷数据进行分解处理，得到多个不同频率的固有模态函数分量和一个剩余残差分量；
20.分量分类模块，用于将多个不同频率的固有模态函数分量分类为高频分量子序列及低频分量子序列；
21.时序特征提取模块，用于通过lstm预测模型提取高频分量子序列的时序特征；
22.高频分量子序列预测模块，用于将lstm预测模型提取的高频分量子序列的时序特征序列、原始电力负荷数据中预测点前k个历史电力负荷数据、高频分量子序列中预测点前k个历史电力负荷数据输入到elm-catboost混合预测模型中，获得高频分量子序列的预测结果；其中，所述elm-catboost混合预测模型包括第一elm预测模型及catboost预测模型；
23.低频分量子序列预测模块，用于通过第二elm预测模型获取低频分量子序列的预测结果；
24.叠加模块，用于将高频分量子序列的预测结果和低频分量子序列的预测结果进行叠加，得到叠加结果；
25.反归一化模块，用于将所述叠加结果进行反归一化处理后，得到基于混合模型的短期电力负荷预测结果。
26.综上所述，本发明基于混合模型的短期电力负荷预测方法及系统通过lstm 预测模型提取高频分量子序列的时序特征，配合catboost预测模型和第一 elm预测模型组成的elm-catboost混合预测模型对短期电力负荷进行预测，利用ceemdan分解算法把非平稳非线性的原始电力负荷数据分解成若干个不同频率且相对平稳的固有模态函数分量，降低了模型预测难度，提高了预测准确度；另外，利用lstm预测模型提取高频分量子序列的时序特征，结合高频分量子序列的历史电力负荷数据和原始电力负荷数据共同作为 elm-catboost混合预测模型的输入特征极大的丰富了输入特征维度信息，通过使用elm-catboost混合预测模型集成了单一模型优势，拥有更高的鲁棒性和准确性，针对高、低频分量子序列采取不同的输入特征和预测模型能够降低模型复杂度，节省计算开销和提高模型预测性能。
附图说明
27.图1是本发明实施例提供的第一种基于混合模型的短期电力负荷预测方法的流程示意图；
28.图2是本发明实施例提供的第二种基于混合模型的短期电力负荷预测方法的流程
示意图；
29.图3是本发明实施例提供的第一种基于混合模型的短期电力负荷预测系统的结构框图；
30.图4是本发明实施例提供的第二种基于混合模型的短期电力负荷预测系统的结构框图；
31.图5是本发明实施例提供的一种电子设备的内部结构框图；
32.图6是本发明实施例提供的elm-catboost混合预测模型的工作原理图。
具体实施方式
33.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
34.需要说明的是，本发明的说明书和权利要求书中的术语“第一”、“第二”、“第三”、“第四”等是用于区别不同的对象，而不是用于描述特定顺序。本发明实施例的术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，示例性地，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
35.如图1是本发明实施例提供的第一种基于混合模型的短期电力负荷预测方法的流程示意图，如图1所示，该基于混合模型的短期电力负荷预测方法，包括步骤s110-步骤s160，具体如下：
36.步骤s110、采集电力负荷数据，组成电力负荷数据集，并对电力负荷数据进行归一化处理，获取归一化处理后的电力负荷数据，记为原始电力负荷数据；其中，采集的电力负荷数据为短期电力负荷历史数据。
37.所述步骤s110、采集电力负荷数据，组成电力负荷数据集，对电力负荷数据进行归一化处理，获取归一化处理后的电力负荷数据的方法，具体操作为：
38.通过公式将电力负荷数据集中的各个电力负荷数据分别映射到 [0，1]，获取归一化处理后的电力负荷数据x
*
；其中，x
*
为归一化处理后的电力负荷数据，即原始电力负荷数据；x
min
、x
max
分别为电力负荷数据集中对应的最小值、最大值的电力负荷数据，x为电力负荷数据集中的电力负荷数据。
[0039]
步骤s120、通过ceemdan分解算法对原始电力负荷数据进行分解处理，得到多个不同频率的固有模态函数分量和一个剩余残差分量。
[0040]
所述步骤s120的方法，具体包括如下步骤：
[0041]
步骤s121、在原始电力负荷数据中依次添加服从正态分布的不同幅值白噪声，构建待处理负荷序列xi(t)＝x(t)+β0ωi(t)；
[0042]
其中，xi(t)为原始电力负荷数据加入第i次白噪声后的待处理电力负荷数据，记为待处理负荷序列，x(t)为原始电力负荷数据中第t个样本数据，β0为噪声系数，ωi(t)为第i次处理对应的白噪声的幅值，1≤i≤p，p为原始电力负荷数据加入白噪声后对应的待处理
负荷序列的总数；
[0043]
步骤s122、通过emd分解方法分解待处理负荷序列，根据p个待处理负荷序列获取的imf1′
1-imf1′
p
取均值，得到第一个均值作为第一个固有模态函数及得到第一个残差序列
[0044]
步骤s123、定义ej(
·
)为待处理电力负荷数据通过emd分解方法分解后获得的第j个固有模态函数imf分量，通过emd分解方法对新的负荷序列 r1(t)+ε1e1(ωi(t))进行分解，得到第二个均值作为第二个固有模态函数及得到第二个残差序列
[0045]
步骤s124、重复上述步骤，获取剩余的固有模态函数imf
3-imfn分量，以及最后的剩余残差res分量：
[0046][0047][0048]
其中，3≤k≤n，rn(t)为剩余残差res分量。
[0049]
本实施例中，所述步骤s122中通过emd分解方法分解待处理电力负荷数据，将p个待处理负荷序列获取的imf1′
1-imf1′
p
取均值的方法，其可以包括以下步骤：
[0050]
第一步、找出待处理负荷序列xi(t)的所有极值点，并利用三次样条插值函数拟合出待处理负荷序列xi(t)的上、下包络线，计算上、下包络线的平均值将待处理负荷序列xi(t)和平均值ei(t)作差，得到差值ci(t)＝xi(t)-ei(t)；其中，e
1i
(t)、e
2i
(t)分别为第i次处理对应的待处理负荷序列的上、下包络线，ei(t)为第i次处理对应的待处理负荷序列的上、下包络线的平均值；
[0051]
第二步、如果差值ci(t)小于或等于预设基准值，或者迭代次数达到预设迭代次数，则将差值ci(t)设定为imf1′i；反之，则将ci(t)作为新的待处理负荷序列，重复执行第一步，直至满足得到的差值ci(t)小于或等于预设基准值，或者迭代次数达到预设迭代次数。
[0052]
步骤s130、将多个不同频率的固有模态函数分量分类为高频分量子序列及低频分量子序列。
[0053]
具体地，所述步骤s130的方法，具体操作为：
[0054]
通过公式将多个不同频率的固有模态函数分量分类为高频分量子序列及低频分量子序列；其中，zn表示过零率，n
zeros
表示过零次数，n表示固有模态函数分量对应的信号总长度，当固有模态函数分量的过零率＜0.01时，该固有模态函数分量对应的曲线相对平稳，频率相对较低，该固有模态函数分量记为低频分量子序列；当固有模态函数分量的过零率≥0.01时，该固有模态函数分量对应的曲线开始出现一定的波动，频率逐渐提高，该固有模态函数分量记为高频分量子序列；在离散信号中，如果相邻采样点具有不同代数
符号就表示发生了过零，过零次数与信号总长度的比例即定义为过零率。
[0055]
步骤s140、通过lstm预测模型提取高频分量子序列的时序特征；其中， lstm预测模型由第一层lstm层、第二层lstm层、一个dropout层和一个 dense层组成，其中，第一层lstm层为128维度，第二层lstm层为32维度，每层都采用relu激活函数，一个遗忘率为0.3的dropout层，一个大小为32 的dense层，批大小为64，最大迭代次数为1000，具有早期停止函数；具体地，确定lstm预测模型的输入特征，通过lstm预测模型提取高频分量子序列的时序特征。
[0056]
所述多个不同频率的固有模态函数分量分类为m个高频分量子序列及 n+1-m个的低频分量子序列，所述步骤s140的方法，具体操作包括：
[0057]
每个高频分量子序列对应一个lstm预测模型，即m个高频分量子序列对应m个lstm预测模型，lstm预测模型的第一层lstm层的输入特征为对应高频分量子序列中预测点前168个历史电力负荷数据，将高频分量子序列中预测点前168个历史电力负荷数据作为输入特征输入到lstm预测模型的第一层 lstm层中，得到第一层lstm层的输出结果；
[0058]
根据第一层lstm层的输出结果获取第二层lstm层的输出结果，作为 lstm预测模型提取高频分量子序列的时序特征；lstm预测模型中第二层 lstm层的输出结果维度为32，每个lstm预测模型中第二层lstm层的输出结果分别对应该lstm预测模型提取的对应高频分量子序列的时序特征。
[0059]
步骤s150、将lstm预测模型提取的高频分量子序列的时序特征序列、原始电力负荷数据中预测点前k个历史电力负荷数据、高频分量子序列中预测点前k个历史电力负荷数据输入到elm-catboost混合预测模型中，获得高频分量子序列的预测结果；其中，所述elm-catboost混合预测模型包括第一 elm预测模型及catboost预测模型，第一elm预测模型的隐藏层节点个数由网格搜索确定，catboost预测模型学习率为0.03，树的深度为8，树的最大数量为4000。
[0060]
如图6所示，所述步骤s150的方法，具体操作包括：
[0061]
每个高频分量子序列对应一个第一elm预测模型及catboost预测模型，即m个高频分量子序列分别对应m个第一elm预测模型及m个catboost 预测模型，将lstm预测模型提取的高频分量子序列的时序特征序列、原始电力负荷数据中预测点前k个历史电力负荷数据、对应高频分量子序列中预测点前k个历史电力负荷数据作为输入特征分别输入到第一elm预测模型及 catboost预测模型中，得到高频分量子序列的预测结果；其中，针对m个高频分量子序列得到m个高频分量子序列的预测结果，第t个高频分量子序列预测结果c
t
＝(a
t
+b
t
)/2，a
t
为第t个高频分量子序列的第一elm预测模型的输出结果，b
t
为第t个高频分量子序列的catboost预测模型的输出结果；本实施例中，k取值为24。
[0062]
本实施例中，一部分固有模态函数分量频率相对较高，不平稳，记为高频分量子序列，其拥有更多的特征信息，因此需要构造lstm预测模型提取更丰富的输入特征信息，以方便elm-catboost混合预测模型提高预测精度；另一部分固有模态函数分量频率相对较低，且相对平稳，记为低频分量子序列，适合单一模型预测，因此针对低频分量子序列利用elm预测模型进行预测。
[0063]
步骤s160、通过第二elm预测模型获取低频分量子序列的预测结果，第二 elm预测模型的隐藏层节点个数由网格搜索确定。
[0064]
所述步骤s160的方法，具体操作包括：
[0065]
每个低频分量子序列对应一个第二elm预测模型，把低频分量子序列中预测点前168个历史负荷数据作为输入特征输入到对应的第二elm预测模型中，得到低频分量子序列的预测结果；针对n+1-m个低频分量子序列分量得到n+1-m 个低频分量子序列的预测结果。
[0066]
步骤s170、将高频分量子序列的预测结果和低频分量子序列的预测结果进行叠加，得到叠加结果，所述叠加结果为高频分量子序列的预测结果和低频分量子序列的预测结果之和。
[0067]
步骤s180、将所述叠加结果进行反归一化处理后，得到基于混合模型的短期电力负荷预测结果。
[0068]
如图2所示，在一个实施例中，所述基于混合模型的短期电力负荷预测方法，还包括
[0069]
步骤s100、分别创建lstm预测模型、elm-catboot混合预测模型及第二elm预测模型，并对其进行训练；其中，所述elm-catboot混合预测模型包括第一elm预测模型及catboot预测模型。
[0070]
在本实施例中，步骤s100可描述成位于步骤s120之前，还可根据需要描述成位于步骤s110或步骤s130之前。
[0071]
所述步骤s100、分别创建lstm预测模型、elm-catboot混合预测模型及第二elm预测模型，并对其进行训练的方法，具体包括以下步骤：
[0072]
步骤s101、采集电力负荷数据，组成电力负荷数据集，并对电力负荷数据进行归一化处理，获取归一化处理后的电力负荷数据，记为原始电力负荷数据；其中，采集的电力负荷数据为短期电力负荷历史数据；本实施例中，以马来西亚某地2009年1月1日到2011年1月1日的电力负荷数据集为例，该电力负荷数据集包含17519个电力负荷数据，采样频率为1小时，选取整个电力负荷数据集的前80％的电力负荷数据作为训练集，另外20％的电力负荷数据则作为测试集。
[0073]
其中，对电力负荷数据进行归一化处理的过程与上述步骤s110类似，这里不再赘述。
[0074]
步骤s102、通过ceemdan分解算法对原始电力负荷数据进行分解处理，得到多个不同频率的固有模态函数分量和一个剩余残差分量。
[0075]
该过程与上述步骤s120类似，这里不再赘述。
[0076]
步骤s103、将多个不同频率的固有模态函数分量分类为高频分量子序列及低频分量子序列。
[0077]
该过程与上述步骤s130类似，这里不再赘述。
[0078]
步骤s104、创建lstm预测模型、elm-catboost混合预测模型及第二 elm预测模型；针对每个高频分量子序列分别创建一个lstm预测模型，所述elm-catboost混合预测模型包括第一elm预测模型及catboost预测模型，针对每个高频分量子序列分别创建一个第一elm预测模型及一个 catboost预测模型，针对每个低频分量子序列分别创建一个第二elm预测模型。
[0079]
步骤s105、通过lstm预测模型提取高频分量子序列的时序特征；确定 lstm预测模型的输入特征，通过lstm预测模型提取高频分量子序列的时序特征。
[0080]
该过程与上述步骤s140类似，这里不再赘述。
[0081]
步骤s106、将lstm预测模型提取的高频分量子序列的时序特征序列、原始电力负荷数据中预测点前k个历史电力负荷数据、高频分量子序列中预测点前k个历史电力负荷数据输入到elm-catboost混合预测模型中，获得高频分量子序列的预测结果。
[0082]
该过程与上述步骤s150类似，这里不再赘述。
[0083]
步骤s107、通过第二elm预测模型获取低频分量子序列的预测结果。
[0084]
该过程与上述步骤s160类似，这里不再赘述。
[0085]
步骤s108、将高频分量子序列的预测结果和低频分量子序列的预测结果进行叠加，得到训练结果，训练结果为高频分量子序列的预测结果和低频分量子序列的预测结果之和。
[0086]
步骤s109、将所述训练结果进行反归一化处理后，采用一种或多种的误差评价指标评估模型比对训练结果和真实结果，直至训练结果和真实结果的误差小于预设阈值，完成对lstm预测模型、elm-catboot混合预测模型及第二 elm预测模型的训练。
[0087]
常用的误差评价指标评估模型可以是平均绝对百分比误差mape，平均绝对误差mae、均方根误差rmse等，这里不做限定，具体如下：
[0088][0089][0090][0091]
其中，表示预测模型第i个样本的训练结果对应的预测值，yi表示预测模型第i个样本的训练结果对应的真实值，s为预测模型中样本的总数，预测模型可为lstm预测模型、elm-catboot混合预测模型或第二elm预测模型。
[0092]
如图2所示，为了更加清晰本发明的技术方案，下面再阐述优选实施例。
[0093]
步骤s100、分别创建lestm预测模型、elm-catboot混合预测模型及第二elm预测模型，并对其进行训练；
[0094]
步骤s110、采集电力负荷数据，组成电力负荷数据集，并对电力负荷数据进行归一化处理，获取归一化处理后的电力负荷数据，记为原始电力负荷数据；其中，采集的电力负荷数据为短期电力负荷历史数据；
[0095]
步骤s120、通过ceemdan分解算法对原始电力负荷数据进行分解处理，得到多个不同频率的固有模态函数分量和一个剩余残差分量；
[0096]
步骤s130、将多个不同频率的固有模态函数分量分类为高频分量子序列及低频分量子序列；
[0097]
步骤s140、通过lstm预测模型提取高频分量子序列的时序特征；
[0098]
步骤s150、将lstm预测模型提取的高频分量子序列的时序特征序列、原始电力负荷数据中预测点前k个历史电力负荷数据、高频分量子序列中预测点前k个历史电力负荷数据输入到预先训练得到的elm-catboost混合预测模型中，获得高频分量子序列的预测结果；
[0099]
步骤s160、通过第二elm预测模型获取低频分量子序列的预测结果；
[0100]
步骤s170、将高频分量子序列的预测结果和低频分量子序列的预测结果进行叠加，得到叠加结果；
[0101]
步骤s180、将所述叠加结果进行反归一化处理后，得到基于混合模型的短期电力负荷预测结果。
[0102]
首先，我们通过案例1验证ceemdan分解预处理的可行性。在这个案例中，lstm、elm、catboost等方法表示单一预测模型。ceemdan-lstm,ceemdan-elm,ceemdan-catboost等方法表示原始电力负荷数据经过ceemdan分解预处理技术分解处理的混合预测模型。预测结果如表1所示。
[0103]
表1带有数据预处理技术的混合模型与单一模型的性能比较
[0104][0105][0106]
通过表1，可以看出单一模型elm、catboost、lstm的rmse分别为 1023.11mw、986.73mw、1157.08mw，而对应混合模型的rmse分别为 830.42mw、901.72mw、825.10mw，分别降低了192.69mw、85.01mw、 331.98mw，并且mae、mape的分数也都优于单一预测模型。说明使用 ceemdan分解预处理技术把非线性、非平稳电力负荷序列分解成若干个不同频率大小的固有模态函数分量，能够丰富特征信息维度，有效提高了混合模型预测性能。
[0107]
通过案例2验证了针对不同子序列设计不同输入特征序列和不同预测模型能够有效提升模型预测性能。在案例2中，ceemdan-elm表示所有的子序列都采用对应子序列前24个历史负荷数据作为输入特征序列， ceemdan-lstm-elm表示所有子序列的输入特征序列由对应子序列前24个历史电力负荷数据、原始电力负荷数据前24个历史电力负荷数据和lstm预测模型提取的时序特征序列组成。ceemdan-d-lstm-elm表示不同子序列的输入特征序列不同。高频分量子序列的输入特征序列由高频分量子序列中预测点前24个历史电力负荷数据、原始电力负荷数据中预测点前24个历史电力负荷数据和lstm预测模型提取的时序特征序列组成。低频分量子序列的输入特征为对应低频分量子序列中预测点前24个历史电
力负荷数据。预测结果如表2所示。
[0108]
表2高频分量子序列混合特征输入与单一特征输入的性能比较
[0109][0110][0111]
通过表2，可以看出混合特征输入模型的rmse、mape、mae都比单一特征输入模型低。说明lstm预测模型能够很好的提取负荷序列中的时序特征，原始电力负荷数据中也包含了一些影响子序列预测性能的隐含特征。此外，我们也可以看出ceemdan-d-lstm-elm预测模型的rmse、mape、mae均最低，说明了不同的子序列拥有不同的特征信息，需要设计不同的输入特征序列。
[0112]
通过案例3验证了针对高低频分量子序列使用不同预测模型能够显著降低模型计算时间和提高模型预测准确性。在案例3中，采用马来西亚电力负荷数据集进行试验比较。ceemdan-lstm-elm-catboost表示针对所有子序列均采用相同预测模型进行预测。ceemdan-lstm-elm-catboost-elm表示针对高频分量子序列采用lstm-elm-catboost混合预测模型进行预测，针对低频分量子序列采用elm预测模型进行预测。预测结果如表3所示。
[0113]
表3高低频分量子序列不同预测模型与相同预测模型的性能比较
[0114][0115]
通过表3，可以看出高低频分量子序列不同预测模型的计算时间明显少于相同预测模型的计算时间。ceemdan-lstm-elm-catboost的计算时间为 6388.15秒，而ceemdan-lstm-elm-catboost-elm的计算时间为2852.40 秒，减少了一半多的时间，此外rmse、mae、mape也都有不同程度的降低，说明了针对不同子序列设计不同的预测模型能够有效降低模型计算时间和提高模型预测性能。因为高频分量子序列包含了更多的特征信息，组合模型能够学习到更加复杂的隐含特征，而低频分量子序列更平缓，单一模型就能有很好的预测性能。
[0116]
由上可以看出，本发明公开了一种基于混合模型的短期电力负荷预测方法，通过lstm预测模型提取高频分量子序列的时序特征，配合catboost预测模型和第一elm预测模型组成的elm-catboost混合预测模型对短期电力负荷进行预测，利用ceemdan分解算法把非平稳非线性的原始电力负荷数据分解成若干个不同频率且相对平稳的固有模态函数分
量，降低了模型预测难度，提高了预测准确度；另外，利用lstm预测模型提取高频分量子序列的时序特征，结合高频分量子序列的历史电力负荷数据和原始电力负荷数据共同作为 elm-catboost混合预测模型的输入特征极大的丰富了输入特征维度信息，通过使用elm-catboost混合预测模型集成了单一模型优势，拥有更高的鲁棒性和准确性，针对高、低频分量子序列采取不同的输入特征和预测模型能够降低模型复杂度，节省计算开销和提高模型预测性能。
[0117]
图3是本发明实施例提供的一种基于混合模型的短期电力负荷预测系统的结构框图，如图3所示，对应于上述基于混合模型的短期电力负荷预测方法，本发明还提供一种基于混合模型的短期电力负荷预测系统，该基于混合模型的短期电力负荷预测系统包括用于执行上述基于混合模型的短期电力负荷预测方法的模块，该系统可以被配置于电子设备等终端，应用本发明基于混合模型的短期电力负荷预测系统，通过lstm预测模型提取高频分量子序列的时序特征，配合catboost预测模型和第一elm预测模型组成的elm-catboost混合预测模型对短期电力负荷进行预测，利用ceemdan分解算法把非平稳非线性的原始电力负荷数据分解成若干个不同频率且相对平稳的固有模态函数分量，降低了模型预测难度，提高了预测准确度；另外，利用lstm预测模型提取高频分量子序列的时序特征，结合高频分量子序列的历史电力负荷数据和原始电力负荷数据共同作为elm-catboost混合预测模型的输入特征极大的丰富了输入特征维度信息，通过使用elm-catboost混合预测模型集成了单一模型优势，拥有更高的鲁棒性和准确性，针对高、低频分量子序列采取不同的输入特征和预测模型能够降低模型复杂度，节省计算开销和提高模型预测性能。
[0118]
具体地，如图3所示，该基于混合模型的短期电力负荷预测系统包括采集模块110、分解模块120、分量分类模块130、时序特征提取模块140、高频分量子序列预测模块150、低频分量子序列预测模块160、叠加模块170及反归一化模块180。
[0119]
采集模块110，用于采集电力负荷数据，组成电力负荷数据集，并对电力负荷数据进行归一化处理，获取归一化处理后的电力负荷数据，记为原始电力负荷数据；
[0120]
分解模块120，用于通过ceemdan分解算法对原始电力负荷数据进行分解处理，得到多个不同频率的固有模态函数分量和一个剩余残差分量；
[0121]
分量分类模块130，用于将多个不同频率的固有模态函数分量分类为高频分量子序列及低频分量子序列；
[0122]
时序特征提取模块140，用于通过lstm预测模型提取高频分量子序列的时序特征；
[0123]
高频分量子序列预测模块150，用于将lstm预测模型提取的高频分量子序列的时序特征序列、原始电力负荷数据中预测点前k个历史电力负荷数据、高频分量子序列中预测点前k个历史电力负荷数据输入到elm-catboost混合预测模型中，获得高频分量子序列的预测结果；其中，所述elm-catboost 混合预测模型包括第一elm预测模型及catboost预测模型；
[0124]
低频分量子序列预测模块160，用于通过第二elm预测模型获取低频分量子序列的预测结果；
[0125]
叠加模块170，用于将高频分量子序列的预测结果和低频分量子序列的预测结果进行叠加，得到叠加结果；
[0126]
反归一化模块180，用于将所述叠加结果进行反归一化处理后，得到基于混合模型
的短期电力负荷预测结果。
[0127]
如图4所示，在一个实施例中，该基于混合模型的短期电力负荷预测系统还包括模型创建模块190，用于分别创建lstm预测模块、elm-catboot混合预测模型及第二elm预测模型，并对其进行训练。
[0128]
图5是本发明实施例提供的一种电子设备的内部结构框图，如图5所示，本发明提供的电子设备包括通过系统总线连接的通过系统总线连接的存储器、处理器；所述存储器上存储有计算机程序，所述处理器用于提供计算和控制能力，以支撑整个电子设备的运行，所述处理器执行所述计算机程序时实现上述的一种基于混合模型的短期电力负荷预测方法。
[0129]
存储器可以包括非易失性存储介质和内存储器，该非易失性存储介质存储有操作系统，还可存储有计算机程序，该计算机程序被处理器执行时，可使得处理器实现基于混合模型的短期电力负荷预测方法。
[0130]
该内存储器中也可储存有计算机程序，该计算机程序被处理器执行时，可使得处理器执行基于混合模型的短期电力负荷预测方法。本领域技术人员可以理解，图5中示出的结构，仅仅是与本技术方案相关的部分结构的框图，并不构成对本技术方案所应用于其他的电子设备的限定，具体的电子设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
[0131]
在一个实施例中，本技术提供的基于混合模型的短期电力负荷预测方法可实现为一种计算机程序的方式，计算机程序可以在如图5所示的电子设备上运行。电子设备的存储器中可存储组成该基于混合模型的短期电力负荷预测系统的各个程序模块，比如，图3所示的采集模块110、分解模块120、分量分类模块130、时序特征提取模块140、高频分量子序列预测模块150、低频分量子序列预测模块160、叠加模块170及反归一化模块180。各个程序模块构成的计算机程序使得处理器执行本说明中描述的本技术各个实施例的基于混合模型的短期电力负荷预测系统的步骤。例如，图5所示的电子设备可以通过如图3所示的基于混合模型的短期电力负荷预测系统中的采集模块110采集电力负荷数据，组成电力负荷数据集，并对电力负荷数据进行归一化处理，获取归一化处理后的电力负荷数据，记为原始电力负荷数据；通过分解模块通过ceemdan分解算法对原始电力负荷数据进行分解处理，得到多个不同频率的固有模态函数分量和一个剩余残差分量；通过分量分类模块130将多个不同频率的固有模态函数分量分类为高频分量子序列及低频分量子序列；通过时序特征提取模块140 通过lstm预测模型提取高频分量子序列的时序特征；通过高频分量子序列预测模块150将lstm预测模型提取的高频分量子序列的时序特征序列、原始电力负荷数据中预测点前k个历史电力负荷数据、高频分量子序列中预测点前k 个历史电力负荷数据输入到elm-catboost混合预测模型中，获得高频分量子序列的预测结果；其中，所述elm-catboost混合预测模型包括第一elm 预测模型及catboost预测模型；通过低频分量子序列预测模块160通过第二 elm预测模型获取低频分量子序列的预测结果；通过叠加模块170将高频分量子序列的预测结果和低频分量子序列的预测结果进行叠加，得到叠加结果；通过反归一化模块180将所述叠加结果进行反归一化处理后，得到基于混合模型的短期电力负荷预测结果。
[0132]
在一个实施例中，提出了一种电子设备，包括存储器和处理器，所述存储器和处理器，所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时，使得所述处理器
执行以下步骤：步骤s110、采集电力负荷数据，组成电力负荷数据集，并对电力负荷数据进行归一化处理，获取归一化处理后的电力负荷数据，记为原始电力负荷数据；其中，采集的电力负荷数据为短期电力负荷历史数据；步骤s120、通过ceemdan分解算法对原始电力负荷数据进行分解处理，得到多个不同频率的固有模态函数分量和一个剩余残差分量；步骤 s130、将多个不同频率的固有模态函数分量分类为高频分量子序列及低频分量子序列；步骤s140、通过lstm预测模型提取高频分量子序列的时序特征；步骤s150、将lstm预测模型提取的高频分量子序列的时序特征序列、原始电力负荷数据中预测点前k个历史电力负荷数据、高频分量子序列中预测点前k个历史电力负荷数据输入到预先训练得到的elm-catboost混合预测模型中，获得高频分量子序列的预测结果；步骤s160、通过第二elm预测模型获取低频分量子序列的预测结果；步骤s170、将高频分量子序列的预测结果和低频分量子序列的预测结果进行叠加，得到叠加结果；步骤s180、将所述叠加结果进行反归一化处理后，得到基于混合模型的短期电力负荷预测结果。
[0133]
应当理解，在本技术实施例中，处理器可以是中央处理单元(centralprocessing unit，cpu)，该处理器还可以是其他通用处理器、数字信号处理器 (digital signal processor，dsp)、专用集成电路(application specific integratedcircuit，asic)、现成可编程门阵列(field-programmable gate array，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。其中，通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
[0134]
本领域普通技术人员可以理解的是实现上述实施例的方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成。该计算机程序包括程序指令，计算机程序可存储于一存储介质中，该存储介质为计算机可读存储介质。该程序指令被该电子设备中的至少一个处理器执行，以实现上述方法的实施例的流程步骤。
[0135]
因此，本发明还提供一种存储介质。该存储介质可以为计算机可读存储介质。该存储介质存储有计算机程序，其中计算机程序包括程序指令。该程序指令被处理器执行时使处理器执行如下步骤：步骤s110、采集电力负荷数据，组成电力负荷数据集，并对电力负荷数据进行归一化处理，获取归一化处理后的电力负荷数据，记为原始电力负荷数据；其中，采集的电力负荷数据为短期电力负荷历史数据；步骤s120、通过ceemdan分解算法对原始电力负荷数据进行分解处理，得到多个不同频率的固有模态函数分量和一个剩余残差分量；步骤s130、将多个不同频率的固有模态函数分量分类为高频分量子序列及低频分量子序列；步骤s140、通过lstm预测模型提取高频分量子序列的时序特征；步骤s150、将lstm预测模型提取的高频分量子序列的时序特征序列、原始电力负荷数据中预测点前k个历史电力负荷数据、高频分量子序列中预测点前k 个历史电力负荷数据输入到预先训练得到的elm-catboost混合预测模型中，获得高频分量子序列的预测结果；步骤s160、通过第二elm预测模型获取低频分量子序列的预测结果；步骤s170、将高频分量子序列的预测结果和低频分量子序列的预测结果进行叠加，得到叠加结果；步骤s180、将所述叠加结果进行反归一化处理后，得到基于混合模型的短期电力负荷预测结果。
[0136]
所述存储介质可以是u盘、移动硬盘、只读存储器(read-only memory，rom)、磁碟或者光盘等各种可以存储程序代码的计算机可读存储介质。
[0137]
综上所述，本发明一种基于混合模型的短期电力负荷预测方法及系统通过 lstm
预测模型提取高频分量子序列的时序特征，配合catboost预测模型和第一elm预测模型组成的elm-catboost混合预测模型对短期电力负荷进行预测，利用ceemdan分解算法把非平稳非线性的原始电力负荷数据分解成若干个不同频率且相对平稳的固有模态函数分量，降低了模型预测难度，提高了预测准确度；另外，利用lstm预测模型提取高频分量子序列的时序特征，结合高频分量子序列的历史电力负荷数据和原始电力负荷数据共同作为 elm-catboost混合预测模型的输入特征极大的丰富了输入特征维度信息，通过使用elm-catboost混合预测模型集成了单一模型优势，拥有更高的鲁棒性和准确性，针对高、低频分量子序列采取不同的输入特征和预测模型能够降低模型复杂度，节省计算开销和提高模型预测性能。
[0138]
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
[0139]
在本发明所提供的几个实施例中，应该理解到，所揭露的系统和方法，可以通过其它的方式实现。例如，以上所描述的系统实施例仅仅是示意性的。例如，各个单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式。例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。
[0140]
本发明实施例方法中的步骤可以根据实际需要进行顺序调整、合并和删减。本发明实施例装置中的单元可以根据实际需要进行合并、划分和删减。另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。该集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，终端，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。
[0141]
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明的保护范围应以所附权利要求为准。