本发明涉及电力系统负荷预测技术领域，具体而言涉及一种基于改进shapley值模型的短期负荷预测方法。
背景技术：
在电力系统运行控制和计划管理中，负荷预测决定了发电、输电和配电的合理安排，是电力系统规划的重要组成部分。其中，短期符合预测最主要的应用是为发电计划程序提供数据，用来确定满足安全要求、运行约束、以及自然环境和设备的运行方案。如何提高预测精度是目前研究短期负荷预测理论与方法的中心和重点。
有部分研究人员针对单一预测方法预测精度不高的问题，提出了组合预测模型进行预测的理念，针对单一预测方法同时对呈现不同平稳特征的三种分量进行预测将产生较大误差的问题，组合预测方法可以有效结合各个预测方法的优势，提高预测精度，难点是权值的分配问题。
例如专利号为cn201811334912.7的发明专利中提出一种基于改进shapley值的风电汇聚趋势性分状态量化方法，针对仅利用测风数据及气象数据对汇聚效应预测精度较低、传统shapley值法在单一模型预测结果偏差过大时仍参与组合的问题，本发明基于改进shapley值法对风电汇聚效应的趋势性进行量化分析，在对不同风电输出状态量化分析的基础上，得到各风电输出状态下的风电持续出力曲线，进而构建基于汇聚特性分析的风电持续出力曲线分状态组合预测模型，并建立预测精度评价体系。相对于单一的预测模型，风电持续出力曲线的分状态组合预测方法能更准确地描述风电汇聚的趋势，为风电基地扩建后的外送输电容量规划提供一定的理论依据。但是这类方法比较适合稳定运行环境下的电力系统负荷预测，并不适用于实际运行过程中出现不稳定特征的复杂负荷数据的预测能力，其中包含的权值分配方法也不适用于短期负荷预测方法。
而随着国家电网的进一步推进，需求侧用户终端的智能化水平不断提高，用户负荷时序特征显著改变，短期负荷数据具有非平稳性的特点，单一的负荷预测模型和常规的组合预测模型忽略负荷数据的时序性特点，难以达到满意的预测准确度。新形势下如何更加精准科学地分析和挖掘负荷的内部特征，提高短期负荷预测的精确性亟待解决。
技术实现要素：
本发明目的在于提供一种基于改进shapley值模型的短期负荷预测方法，通过希尔伯特黄(hht)变换对非平稳负荷进行重构得到随机、周期、趋势分量；通过改进shapley值模型确定组合预测各个预测方法的权重分配，并分别应用于随机、周期、趋势分量的预测，将得到的各个预测分量进行叠加得到最终预测值。相对于现有预测方法，本发明提出的预测方法具有更高的精确度和稳定性。
为达成上述目的，结合图1，本发明提出一种基于改进shapley值模型的短期负荷预测方法，所述短期负荷预测方法包括：
s1：通过经验模态分解对呈现不稳定特征的负荷数据进行分解，得到有限个本征模态函数；
s2：通过希尔伯特变换得到各个本征模态函数对应的时间-频率曲线；
s3：按各个时间-频率曲线交叠最少为原则将各个本征模态函数重构为随机、周期、趋势分量；
s4：针对随机、周期、趋势分量分别采用改进shapley值组合预测方法进行预测，得到随机、周期、趋势分量对应的预测分量；
s5：将随机、周期、趋势分量对应的预测分量进行叠加后得到最终预测分量。
进一步的实施例中，步骤s1中，所述通过经验模态分解对呈现不稳定特征的负荷数据进行分解，得到有限个本征模态函数的过程包括以下步骤：
s11：设l＝1；
s12：采用三次样条函数对负荷数据进行插值，得到其上包络线emax(t)以及下包络线emin(t)，计算二者均值ml(t)；
s13：通过负荷功率s(t)分解得到候选模态函数序列cl(t)：
cl(t)＝s(t)-ml(t)
s14：判断cl(t)是否满足由局部最大值点及最小值点分别拟合生成的上包络线、下包络线的均值等于零，如果不满足，则将cl(t)视为新的s(t)，重复步骤s12至s13，直至满足条件得到模态函数序列cl(t)；
s15：从负荷曲线s(t)中去除cl(t)分量得到余项r(t)，计算公式如下：
r(t)＝s(t)-cl(t)；
s16：将r(t)视为新的s(t)，l加1，重复步骤s12至s16，计算得到其他候选模态函数序列cl(t)，直至r(t)无法再分离模态函数序列；
此时，各固有模态函数cl(t)和余项r(t)满足以下公式：
其中，num是候选模态函数序列的总数，s(t)是初始负荷功率。
进一步的实施例中，步骤s2中，所述通过希尔伯特变换得到各个本征模态函数对应的时间-频率曲线的过程包括以下步骤：
s21：将本征模态函数cl(t)与
s22：以本征模态函数cl(t)与全通移相网络
其中，al(t)是振幅函数，
s23：计算振幅函数al(t)，计算公式如下：
s24：计算瞬时频率fl(t)，计算公式如下：
进一步的实施例中，步骤s4中，所述针对随机、周期、趋势分量分别采用改进shapley值组合预测方法进行预测，得到随机、周期、趋势分量对应的预测分量的过程包括以下步骤：
s41：采用多种预测方法针对分量k进行预测，获取预测误差集合fk；
s42：通过针对分量k的预测误差集合fk得到其所有子集，即组合误差值ek(*)；
s43：计算各个预测方法针对分量k的shapley值eik；
s44：计算得到针对分量k的预测权值集合wk；
s45：根据针对分量k的预测权值集合wk，计算得到分量k的预测结果。
进一步的实施例中，步骤s44中，所述计算得到针对分量k的预测权值集合wk还包括以下步骤：
判断预测权值集合wk中是否存在小于0的元素，如果存在，去除权值为负值的预测方法，返回步骤s42，重新计算针对分量k的预测权值集合wk，否则，进入步骤s45，计算得到分量k的预测结果。
进一步的实施例中，步骤s4中，所述针对随机、周期、趋势分量分别采用改进shapley值组合预测方法进行预测，得到随机、周期、趋势分量对应的预测分量的过程包括以下步骤：
采用灰色预测、elman神经网络和svm进行组合预测，通过改进shapley值法确定各预测方法的权重分配，组合预测模型如下：
yk＝ωakyak ωbkybk ωckyck
式中：yk为状态k下的组合预测值；ωak、ωbk、ωck分别为灰色预测、elman神经网络和svm三种在分量k下的权重；yak、ybk、yck分别为三种预测方法在分量k下的预测值。
进一步的实施例中，所述采用灰色预测、elman神经网络和svm进行组合预测，通过改进shapley值法确定各预测方法的权重分配的过程包括以下步骤：
s401：设集合i＝{1，2，3}，集合i包含的任何子集p、q，其表示为3种预测方法中的任一组合，e(p)、e(q)表示各内部缺陷预测误差；
对于集合i的任何子集p、q，都有：
e(p) e(q)≥e(p∪q)；
s402：设
e(3)表示3种预测模型进行组合预测由于内部缺陷而产生的总误差，即e(3)＝∑i∈ixi；
s403：定义ei为第i种预测模型所产生误差绝对值的平均值，e为组合预测的误差总和，则有：
式中：样本个数表示为m，第i种预测模型的第σ个样本产生的误差绝对值表示为eiσ，n为预测模型总数；
s404：根据下述公式基于shapley值法分配权重：
式中：ω|p|表示加权因子，p-{i}表示组合中去除模型i′，i′表示参与组合预测的某个负荷预测模型；e′t为第i种预测模型的误差量，即shapley值，|p|为组合预测模型内预测方法的个数；
通过下述公式计算组合预测模型中各预测方法权重分配：
以上本发明的技术方案，与现有相比，其显著的有益效果在于：
(1)在具有较大波动特征的负荷预测中，通过hht变换对负荷进行重构，得到波动特征集中且变化频率平稳的分量，对重构分量分别进行预测，有效降低了负荷波动性对负荷预测的影响。
(2)通过改进shapley值法确定组合预测内各预测方法的权值，叠加重构后三种分量的预测值得到最终预测结果，相比单一预测模型和未改进的shapley值组合预测模型，预测准确度和算法稳定性得到很大的提高。
(3)在对负荷数据进行重构后，针对各个分量设置与之对应的权值，分别采用组合预测模型进行预测，相比于现有预测方法，针对性更强，预测准确度更高。
应当理解，前述构思以及在下面更加详细地描述的额外构思的所有组合只要在这样的构思不相互矛盾的情况下都可以被视为本公开的发明主题的一部分。另外，所要求保护的主题的所有组合都被视为本公开的发明主题的一部分。
结合附图从下面的描述中可以更加全面地理解本发明教导的前述和其他方面、实施例和特征。本发明的其他附加方面例如示例性实施方式的特征和/或有益效果将在下面的描述中显见，或通过根据本发明教导的具体实施方式的实践中得知。
附图说明
附图不意在按比例绘制。在附图中，在各个图中示出的每个相同或近似相同的组成部分可以用相同的标号表示。为了清晰起见，在每个图中，并非每个组成部分均被标记。现在，将通过例子并参考附图来描述本发明的各个方面的实施例，其中：
图1为基于改进shapley值模型的短期负荷预测方法的流程图。
图2为用于训练和测试的历史负荷数据的示意图。
图3为历史负荷数据中30日数据经emd分解后的本征模态分量的示意图。
图4为经希尔伯特变换后各个本征模态分量的时间-频率曲线的示意图。
图5为随机、周期、趋势分量情况的示意图。
图6为本方法的预测结果与其他方法预测结果的对比示意图。
图7为采用shapley值模型得到的权值分配结果(表1)。
图8为采用改进shapley值模型得到的权值分配结果(表2)。
图9为各个方法预测结果的对比(表3)。
具体实施方式
为了更了解本发明的技术内容，特举具体实施例并配合所附图式说明如下。
如图1所示，本发明提出的基于改进shapley值模型的短期负荷预测方法包括如下步骤：
步骤1，通过经验模态分解(emd)对呈现不稳定特征的负荷数据进行分解，得到有限个本征模态函数(imf)。
希尔伯特黄(hht)变换可以根据负荷内部不平稳的时序性特征，通过emd分解将呈现不平稳特点的负荷分解为几种具有不同特征时间尺度的固有模态分量，通过希尔伯特变换可以得到各模态分量的瞬时频率-时间曲线，为负荷曲线的重构提供条件。
本方法采用图2所示的负荷数据，提取其中18日数据(每日采集96次)作为算例进行分析，将18*96个数据进行emd分解，得到各个本征模态分量，如图3所示，具体步骤为：
1)采用三次样条函数对初始负荷数据进行插值，得到其上包络线emax(t)以及下包络线emin(t)，计算二者均值m1(t)。
2)通过负荷功率p(t)分解得到候选模态函数序列c1(t)：
c1(t)＝s(t)-m1(t)(1)
3)判断c1(t)是否满足由局部最大值点及最小值点分别拟合生成的上包络线、下包络线的均值等于零，如果不满足，则将c1(t)视为新的s(t)，重复步骤1)、2)，直至满足条件得到模态函数序列c1(t)。
4)从负荷曲线s(t)中去除c1(t)分量得到余项r(t)，计算公式如下：
r(t)＝s(t)-c1(t)(2)
5)将r(t)视为新的s(t)，重复步骤1)、2)，3)，4)，计算得到其他候选模态函数序列cl(t)，直至r(t)无法再分离模态函数序列。
此时，各固有模态函数cl(t)和余项r(t)满足以下公式：
步骤2，通过希尔伯特变换得到各个本征模态函数对应的时间-频率曲线。
对经过emd分解之后的本征模态分量cl做希伯尔特变换，得各个模态分量的瞬时频率，结果如图4所示，主要步骤为：
1)将本征模态函数cl(t)与
2)以本征模态函数cl(t)与全通移相网络
3)计算振幅函数，计算公式如下：
4)计算瞬时频率，计算公式如下：
步骤3，按各个时间-频率曲线交叠最少为原则将各个本征模态函数重构为随机、周期、趋势分量，结果如图5所示。
步骤4，针对随机、周期、趋势分量分别采用改进shapley值组合预测方法进行预测，得到随机、周期、趋势分量对应的预测分量。
组合预测模型如下：
yk＝ωakyak ωbkybk ωckyck(8)
式中：yk为状态k下的组合预测值；ωak,ωbk,ωck分别为灰色预测、elman神经网络和svm三种在分量k下的权重；yak，ybk，yck分别为三种预测方法在分量k下的预测值。
为得到各个预测方法的权值，本专利构造shapley值权值模型，并做以下定义：
1)对于i集合(i＝{1，2，3})包含的任何子集p、q，其表示为3种方法中的任一组合，e(p)、e(q)表示各内部缺陷预测误差。
2)对于i的任何子集p、q，都有：
e(p) e(q)≥e(p∪q)(9)
3)设
4)e(3)表示3种预测模型进行组合预测由于内部缺陷而产生的总误差，即e(3)＝∑i∈ixi。
定义ei为第i种预测模型所产生误差绝对值的平均值，e为组合预测的误差总和，则有：
式中：样本个数表示为m；第i种预测模型的第σ个样本产生的误差绝对值表示为eiσ，n为预测模型总数。
式(12)、(13)表示基于shapley值法分配权重的公式：
式中：ω|p|表示加权因子，p-{i}表示组合中去除模型i′，i′表示参与组合预测的某个负荷预测模型；e′t为i模型的误差量，即shapley值，|p|为组合预测模型内预测方法的个数。
组合预测模型中各预测方法权重分配可以通过式(14)计算：
针对某一预测方法误差过大影响预测精度，对shapley值模型进行改进，本方法舍弃负权值系数代表的预测方法，对剩余的组合预测剩余方法重新进行权重分配，避免预测结果较差的预测方法干扰组合预测结果。采用shapley值模型得到的权值分配结果结果如表1所示。采用改进shapley值模型得到的权值分配结果如表2所示。
步骤5，针对随机、周期、趋势分量对应的预测分量进行叠加得到最终预测分量。
为了验证本发明提出的预测方法的性能，本发明采用单一预测模型、未改进的shapley组合预测模型和改进后shapley值的组合预测模型三种方案对非平稳负荷进行短期预测，从模型精确度和稳定性两个角进行对比分析。比对结果如图6、表3所示，结果表明，本文提出的预测方法具有更高的精确度和稳定性。
在本公开中参照附图来描述本发明的各方面，附图中示出了许多说明的实施例。本公开的实施例不必定义在包括本发明的所有方面。应当理解，上面介绍的多种构思和实施例，以及下面更加详细地描述的那些构思和实施方式可以以很多方式中任意一种来实施，这是因为本发明所公开的构思和实施例并不限于任何实施方式。另外，本发明公开的一些方面可以单独使用，或者与本发明公开的其他方面的任何适当组合来使用。
虽然本发明已以较佳实施例揭露如上，然其并非用以限定本发明。本发明所属技术领域中具有通常知识者，在不脱离本发明的精神和范围内，当可作各种的更动与润饰。因此，本发明的保护范围当视权利要求书所界定者为准。
1.一种基于改进shapley值模型的短期负荷预测方法，其特征在于，所述短期负荷预测方法包括：
s1：通过经验模态分解对呈现不稳定特征的负荷数据进行分解，得到有限个本征模态函数；
s2：通过希尔伯特变换得到各个本征模态函数对应的时间-频率曲线；
s3：按各个时间-频率曲线交叠最少为原则将各个本征模态函数重构为随机、周期、趋势分量；
s4：针对随机、周期、趋势分量分别采用改进shapley值组合预测方法进行预测，得到随机、周期、趋势分量对应的预测分量；
s5：将随机、周期、趋势分量对应的预测分量进行叠加后得到最终预测分量。
2.根据权利要求1所述的基于改进shapley值模型的短期负荷预测方法，其特征在于，步骤s1中，所述通过经验模态分解对呈现不稳定特征的负荷数据进行分解，得到有限个本征模态函数的过程包括以下步骤：
s11：设l＝1；
s12：采用三次样条函数对负荷数据进行插值，得到其上包络线emax(t)以及下包络线emin(t)，计算二者均值ml(t)；
s13：通过负荷功率s(t)分解得到候选模态函数序列cl(t)：
cl(t)＝s(t)-ml(t)
s14：判断cl(t)是否满足由局部最大值点及最小值点分别拟合生成的上包络线、下包络线的均值等于零，如果不满足，则将cl(t)视为新的s(t)，重复步骤s12至s13，直至满足条件得到模态函数序列cl(t)；
s15：从负荷曲线s(t)中去除cl(t)分量得到余项r(t)，计算公式如下：
r(t)＝s(t)-cl(t)；
s16：将r(t)视为新的s(t)，l加1，重复步骤s12至s16，计算得到其他候选模态函数序列cl(t)，直至r(t)无法再分离模态函数序列；
此时，各固有模态函数cl(t)和余项r(t)满足以下公式：
其中，num是候选模态函数序列的总数。
3.根据权利要求1所述的基于改进shapley值模型的短期负荷预测方法，其特征在于，步骤s2中，所述通过希尔伯特变换得到各个本征模态函数对应的时间-频率曲线的过程包括以下步骤：
s21：将本征模态函数cl(t)与
s22：以本征模态函数cl(t)与全通移相网络
其中，al(t)是振幅函数，
s23：计算振幅函数al(t)，计算公式如下：
s24：计算瞬时频率fl(t)，计算公式如下：
4.根据权利要求1所述的基于改进shapley值模型的短期负荷预测方法，其特征在于，步骤s4中，所述针对随机、周期、趋势分量分别采用改进shapley值组合预测方法进行预测，得到随机、周期、趋势分量对应的预测分量的过程包括以下步骤：
s41：采用多种预测方法针对分量k进行预测，获取预测误差集合fk；
s42：通过针对分量k的预测误差集合fk得到其所有子集，即组合误差值ek(*)；
s43：计算各个预测方法针对分量k的shapley值eik；
s44：计算得到针对分量k的预测权值集合wk；
s45：根据针对分量k的预测权值集合wk，计算得到分量k的预测结果。
5.根据权利要求4所述的基于改进shapley值模型的短期负荷预测方法，其特征在于，步骤s44中，所述计算得到针对分量k的预测权值集合wk还包括以下步骤：
判断预测权值集合wk中是否存在小于0的元素，如果存在，去除权值为负值的预测方法，返回步骤s42，重新计算针对分量k的预测权值集合wk，否则，进入步骤s45，计算得到分量k的预测结果。
6.根据权利要求1所述的基于改进shapley值模型的短期负荷预测方法，其特征在于，步骤s4中，所述针对随机、周期、趋势分量分别采用改进shapley值组合预测方法进行预测，得到随机、周期、趋势分量对应的预测分量的过程包括以下步骤：
采用灰色预测、elman神经网络和svm进行组合预测，通过改进shapley值法确定各预测方法的权重分配，组合预测模型如下：
yk＝ωakyak ωbkybk ωckyck
式中：yk为状态k下的组合预测值；ωak、ωbk、ωck分别为灰色预测、elman神经网络和svm三种在分量k下的权重；yak、ybk、yck分别为三种预测方法在分量k下的预测值。
7.根据权利要求6所述的基于改进shapley值模型的短期负荷预测方法，其特征在于，所述采用灰色预测、elman神经网络和svm进行组合预测，通过改进shapley值法确定各预测方法的权重分配的过程包括以下步骤：
s401：设集合i＝{1，2，3}，集合i包含的任何子集p、q，其表示为3种预测方法中的任一组合，e(p)、e(q)表示各内部缺陷预测误差；
对于集合i的任何子集p、q，都有：
e(p) e(q)≥e(p∪q)；
s402：设
e(3)表示3种预测模型进行组合预测由于内部缺陷而产生的总误差，即e(3)＝∑i∈ixi；
s403：定义ei为第i种预测模型所产生误差绝对值的平均值，e为组合预测的误差总和，则有：
式中：样本个数表示为m，第i种预测模型的第σ个样本产生的误差绝对值表示为eiσ，n为预测模型总数；
s404：根据下述公式基于shapley值法分配权重：
式中：ω|p|表示加权因子，p-{i}表示组合中去除模型i′，i′表示参与组合预测的某个负荷预测模型；e′t为第i种预测模型的误差量，即shapley值，|p|为组合预测模型内预测方法的个数；
通过下述公式计算组合预测模型中各预测方法权重分配：
