基于误差修正的负荷日前预测修正方法、设备和存储介质

1.本发明属于电力负荷预测领域，具体涉及一种基于误差修正的负荷日前预测修正方法、设备和存储介质。


背景技术：

2.电力负荷预测是保证电力系统稳定且经济运行的重要基础，不同时间跨度的负荷预测对电网有着不同的应用目的。其中，精确的短期负荷预测结果能帮助电力系统工作人员制定合理生产计划，维持供需平衡和确保电网安全，同时减少资源浪费和用电成本。负荷水平受日类型、天气、气候、特殊活动等因素的影响。负荷预测技术可分为统计学方法和人工智能方法。人工神经网络是一种软计算技术，它不需要预测人员对底层物理系统进行显式建模，其通过对历史数据的简单学习建模，构建输入变量与电力需求之间的映射关系，进行预测。换句话说，预测者不必指定输入和输出变量之间的函数形式，而在构建多元线性回归模型时必须指定。因此，常通过人工神经网络建模进行电力负荷预测，但由于电网负荷预测具有一定的随机性，稳定性相对较低，影响因素也相对较多，导致误差较大，人工神经网络的应用具有一定的难度，特别是对于规模越小的电网，其存在预测偏差的可能性就越大。
3.园区级短期负荷预测的目的是规划一个最优的电力调度计划，以满足负荷消耗并进行微网能量管理。微网能量管理是根据系统的电负荷、热负荷需求、大气环境、电网电价、燃气价格等信息，协调微网内的可控型分布式电源和储能等设备，保证微网安全稳定运行，实现微网最优经济运行。
4.由于园区级微网负荷规模较小，波动性和不可预测性较强，且受天气、体感参数等可测量和预知的因素较小，而受到随机启动负荷的影响较大。因此尽管短期负荷预测对于微电网很重要，但这方面的研究依然较少。已有一些工作集中在微电网建筑负荷消耗的日前预测上，尚无当日及时根据实际负荷做出反馈和调整的修正算法。这导致之前提出的日前负荷预测算法要么所需信息输入量较大且难以获得，不具备实用性；要么预测精度欠佳且无法应对突发负荷增大的情况。


技术实现要素：

5.本发明的目的在于克服现有技术的缺陷，提供一种基于误差修正的负荷日前预测修正方法、设备和存储介质。
6.为实现上述目的，本发明采用了如下技术方案：
7.一种基于误差修正的的负荷日前预测修正方法，包括以下步骤：
8.步骤1：基于历史负荷和历史气象数据，建立日前预测模型；
9.步骤2：基于日前预测模型得到待预测日预测结果，保存待预测日预测结果和待预测日的实际负荷数据，累积多组的待预测日预测结果及对应的实际负荷数据；
10.步骤3：利用多组的待预测日预测结果及对应的实际负荷数据，采用粒子群算法进行修正算法的关键参数的整定；
11.步骤4：在待预测日当日的修正时段内，定期基于日前预测负荷值、实际负荷值和修正启动偏差门槛值door检查是否满足修正条件，当满足修正条件时启动修正，通过修正算法修正该日的日前预测结果；
12.所述修正算法为：
13.ci＝real
i-predicti，i∈(n，n x*m]
[0014][0015][0016]
其中，reali为第i个采样点的实际负荷值，predicti为第i个采样点的日前预测负荷值，t
now
表示启动修正的时刻，t
end
为修正时段结束时间点，n为t
now
前x个小时采样点数据的编号，m为每个小时内采样点的数目，k为衰减系数，predict
i*
为第i个采样点修正后的日前预测负荷值；第i个采样点时间为ti，采样点i∈(n，n x*m]对应的采样时间段为ti∈(t
now-xh，t
now
]，采样点i∈[n x*m，24*m)对应的修正时间段为ti∈[t
now
，t
end
]。
[0017]
优选的，步骤4中所述的修正条件为，启动修正的时刻前一个时间段内所有采样点的预测负荷值和实际负荷值的偏差的绝对值均大于修正启动偏差门槛值door，而且负荷预测值和实际值的偏差的符号相同。
[0018]
优选的，步骤3中所述的修正算法的关键参数为修正时段结束时间点t
end
、修正启动偏差门槛值door和衰减系数k中的一个或多个。
[0019]
优选的，还包括步骤5，每隔一定时间重复步骤3。
[0020]
优选的，在步骤1中根据日类型分别建立工作日和休息日的日前预测模型。
[0021]
优选的，步骤1中日前预测模型采用卷积神经网络模型，所述历史气象数据输入至卷积神经网络模块的全连接层中；所述历史负荷数据输入至卷积神经网络模块的卷积层中。
[0022]
优选的，所述历史气象数据包括天气数据和天气类型，所述天气数据输入至卷积神经网络模型的两个全连接层，所述天气类型输入至卷积神经网络模型的一个全连接层，所述历史负荷数据输入至卷积神经网络模型的卷积层和池化层，处理后的所有数据合并经过卷积神经网络模块的三个全连接层得到最终输出。
[0023]
优选的，基于粒子群算法整定修正参数的具体步骤如下：
[0024]
步骤31，将累积的多组的待预测日预测结果及对应的实际负荷数据组成的数据集投入粒子群算法；
[0025]
步骤32，初始化参数，设置粒子群算法的参数，包括惯性因子、局部速度因子、全局速度因子和粒子群数量；
[0026]
步骤33，随机初始化若干个粒子，粒子的三维坐标为修正时段结束时间点t
end
、修正启动偏差门槛值door和衰减系数k，设置修正时段结束时间点t
end
、修正启动偏差门槛值door和衰减系数k的变化范围；
[0027]
步骤34，对每个粒子按数据集的平均修正次数和平均精度计算适应度，寻找全局最优解；
[0028]
步骤35，检查是否满足终止条件，如果不满足终止条件则进入步骤36，满足终止条件则进入步骤37；
[0029]
步骤36，基于以下公式进行速度和位置更新，更新后进入步骤34；
[0030]vi，j
(t 1)＝ωv
i，j
(t) c1*rand(0，1)*[pbest
i，j-x
i，j
(t)] c2*rand(0，1)*[gbest
i，j-x
i，j
(t)]
[0031]
x
i，j
(t 1)＝x
i，j
(t) v
i，j
(t 1)
[0032]vi，j
是粒子的速度；pbest和gbest分别定义为粒子最佳位置和种群最佳位置；rand(0，1)是介于(0，1)之间的随机数；x
i，j
(t)是粒子的当前位置；c1和c2是学习因子，ω为压缩因子；
[0033]
步骤37，根据每次迭代得到的最优解拟合曲线得到帕累托最优解，输出最优解的修正时段结束时间点t
end
和修正启动偏差门槛值door和修正函数的衰减系数k。
[0034]
本发明还提供一种微网能量管理设备，包括一个或多个处理器以及一个或多个存储计算机程序的存储器，所述一个或多个处理器，用于执行所述计算机程序，以执行本发明上述的基于误差修正的负荷日前预测修正方法的步骤。
[0035]
本发明还一种存储有计算机程序的存储介质，当所述计算机程序被一个或多个处理器执行时，致使所述一个或多个处理器执行本发明上述的基于误差修正的负荷日前预测修正方法的步骤。
[0036]
本发明的基于误差修正的负荷日前预测修正方法，对园区级微网建立负荷预测模型，基于负荷预测模型对待预测日进行日前负荷预测；在积累多日的预测负荷值和实际负荷值后，通过粒子群算法进行修正算法的关键参数的整定，得到最优参数；将最优参数内置于修正算法，在待预测日当日对负荷预测结果进行即时修正，为微网能量管理系统提供更加准确、灵活的预测结果，从而保证微网安全稳定运行，实现微网最优经济运行。
[0037]
此外，本发明定期的通过粒子群算法优化修正算法中的三个关键参数，修正时段结束时间点t
end
、修正启动偏差门槛值door和衰减系数k，基于本发明的修正算法对负荷预测结果进行即时修正，随着微网系统的不断运行学习，能够得到越来越准确的预测结果。
[0038]
此外，本发明的日前预测模型根据日类型分别建立工作日和休息日的日前预测模，而且针对历史负荷和历史气象数据的不同输入特性，经过不同的输入通道进行处理，再进入全连接层合并，从而在前期提取出更有效的信息进行训练，实现对卷积神经网络的改进，且未增加卷积神经网络概率层的网络深度。
附图说明
[0039]
图1是本发明实施例的负荷日前预测修正方法的步骤示意图；
[0040]
图2是本发明实施例的负荷预测模型示意图；
[0041]
图3是本发明实施例的修正算法参数整定和更新逻辑流程图；
[0042]
图4是本发明实施例的判断是否满足修正条件的逻辑流程图；
[0043]
图5是本发明实施例的原日前预测值、修正后预测值及实际负荷值的对比图；
[0044]
图6是本发明实施例的日前预测与修正后预测对比图。
具体实施方式
[0045]
以下结合附图1至6给出的实施例，进一步说明本发明的基于误差修正的负荷日前预测修正方法的具体实施方式。本发明的基于误差修正的负荷日前预测修正方法不限于以
下实施例的描述。
[0046]
本发明的一种基于误差修正的负荷日前预测修正方法，包括以下步骤：
[0047]
步骤1：基于历史负荷和历史气象数据，建立日前预测模型；
[0048]
步骤2：基于日前预测模型得到待预测日预测结果，保存待预测日预测结果和待预测日的实际负荷数据，累积多组的待预测日预测结果及对应的实际负荷数据；
[0049]
步骤3：利用累积的多组的待预测日预测结果及对应的实际负荷数据，采用粒子群算法进行修正算法的关键参数的整定；优选的，所述的修正算法的关键参数为修正时段结束时间点t
end
、修正启动偏差门槛值door和衰减系数k中的一个或多个。
[0050]
步骤4：在待预测日当日的修正时段内，定期基于日前预测负荷值、实际负荷值和修正启动偏差门槛值door检查是否满足修正条件，当满足修正条件时启动修正，通过修正算法修正该日的日前预测结果；
[0051]
所述修正算法为：
[0052]ci
＝real
i-predicti，i∈(n，n x*m]
[0053][0054][0055]
其中，reali为第i个采样点的实际负荷值，predicti为第i个采样点的日前预测负荷值，t
now
表示启动修正的时刻，t
end
为修正时段结束时间点，n为t
now
前x个小时采样点数据的编号，m为每个小时内采样点的数目，k为衰减系数，predict
i*
为第i个采样点修正后的日前预测负荷值；第i个采样点时间为ti，采样点i∈(n，n x*m]对应的采样时间段为ti∈(t
now-xh，t
now
]，采样点i∈[n x*m，24*m)对应的修正时间段为ti∈[t
now
，t
end
]，-xh表示前x个小时。
[0056]
本发明的基于误差修正的负荷日前预测修正方法，对微网建立负荷预测模型，基于负荷预测模型对待预测日进行日前负荷预测；在积累多日的预测负荷值和实际负荷值后，通过粒子群算法进行修正算法的关键参数的整定，得到最优参数；将最优参数内置于修正算法，在待预测日当日对负荷预测结果进行即时修正，为微网能量管理系统提供更加准确、灵活的预测结果，从而保证微网安全稳定运行，实现微网最优经济运行。
[0057]
本发明的基于误差修正的负荷日前预测修正方法，适用于园区级微电网系统，以下结合附图，对本发明的一个某园区级微电网系统的实施例作详细说明，本发明的流程图如图1所示：
[0058]
考虑某园区级微电网系统，根据上文所述方法，对其进行负荷日前预测修正获得最终预测结果。
[0059]
本实施例在python3.0下编程实现，并调用mysql数据库、keras包实现日前预测部分，将日前预测结果写入数据库；pandas包实现粒子群算法部分，用以得到修正算法三个关键参数的帕累托最优解；在待预测日当日再次通过读写数据库，完成预测结果的实时修正。
[0060]
步骤1，基于历史负荷和历史气象数据，建立日前预测模型。
[0061]
首先累计园区微网一定的历史负荷数据，历史负荷数据可以存入历史负荷数据表；获取历史气象数据，可以通过气象爬虫从网络上抓取该园区位置的历史气象数据，写入气象数据库表。然后基于获取的历史负荷数据和历史气象数据训练日前预测模型，日前预
测模型可以采用人工神经网络，此为本领域的现有技术不再赘述。
[0062]
优选的，步骤1中根据日类型分别建立工作日和休息日的日前预测模型，分别建立工作日日前预测模型和休息日日前预测模型。
[0063]
优选的，步骤1中日前预测模型采用卷积神经网络模型。
[0064]
如图2所示，本实施例中日前预测模型，在输入端引入不同的输入数据处理方式。历史气象数据具有复杂非线性关系，输入至卷积神经网络模块的全连接层中；而历史负荷数据具有较强时序性，输入至卷积神经网络模块的卷积层中，采用了两种输入数据处理方式；处理后的数据再进行合并输入三层全连接层。本实施例将各类不同特性的输入特征能经过不同的输入通道进行处理，再进入全连接层合并，从而在前期提取出更有效的信息进行训练，实现对卷积神经网络的改进，且未增加卷积神经网络概率层的网络深度。当然，根据输入的不同数据，还可以采用三种或更多的数据处理方式。
[0065]
本实施例采用多输入数据处理方式的卷积神经网络模型来建立日前预测模型，输入积累的历史负荷数据和历史气象数据。由于历史气象数据和历史负荷数据均包含待预测日负荷信息，因此每个输入量必须得到符合其时序特征的处理。且由于不同日类型下负荷曲线形状不同，例如工作日和休息日的负荷曲线明显不同，工作日情况下不同的天气类型其负荷曲线肯定也不同，因此需要分类训练和预测，根据日类型分别建立工作日和休息日的日前预测模型。
[0066]
同时，建立日前预测模型时卷积层、池化层的顺序和参数都需要精心设计，才能更好的提取时间序列的数据特征。一般来说，深度越深的神经网络能够在更复杂的任务上获得更高的精度，但也会导致更长的训练时间、优化收敛和过拟合问题。因此，模型设计的本质是在模型的复杂性和准确性之间进行权衡。通常情况下，卷积层和激活函数结合在一起提取输入输出的非线性关系。
[0067]
优选的，所述历史气象数据包括天气数据和天气类型，本实施例的日前预测模型选择了历史气象数据中的天气数据(如气温数据，还可以包括湿度、降水、日照、风力等数据)经过卷积神经网络模块的两个全连接层、历史气象数据中的天气类型经过卷积神经网络模块的一个全连接层，历史负荷数据经过卷积神经网络模块的一个卷积层和一个池化层，通过加入池化层，在不显著降低精度的前提下降低训练参数，优选在卷积层中结合激活函数一起提取输入输出的非线性关系，最终所有数据再合并经过卷积神经网络模块的三个全连接层得到最终输出，即得到日前预测结果。
[0068]
步骤2：基于日前预测模型得到待预测日预测结果，保存待预测日预测结果和待预测日的实际负荷数据，累积多组的待预测日预测结果及对应的实际负荷数据。
[0069]
由待预测日对应输入得到待预测日预测结果，即修正前的待预测日一次预测结果。例如，输入待预测日期：2019年7月11日；日类型：工作日(编码1)；温度：最高温度(25)、最低温度(12)；天气类型：晴(编码例如2)；历史负荷值。历史负荷值如可以输入一周前、两日前、一日前三项历史负荷数据，根据需要也可以输入更多次作为参考，如多周前的同类型日历史负荷值。根据待预测日园区微网的实际运行负荷得到实际负荷数据，均保存到数据库中。对每日的日前预测模型分多个时段进行预测，保存每个时段的预测结果，以及每个时段的实际负荷数据。例如一天划分为24个时段，每个小时一个时段；或者一天划分为96个时段，即15分钟一个时段。
[0070]
参见图1的日前预测主程序部分，本实施例划分为96个时段，进行待预测日预测时首先判断是工作日还是休息日，如果是工作日则采用工作日日前预测模型，如果是休息日则采用休息日预测模型；然后获得96点预测负荷值，写入数据库预测负荷数据表。园区的微网能量管理系统，在运行时，定期将实际负荷数据写入数据库实际负荷数据表，通用可以按96个时段进行保存记录。经过多日的预测和运行，累积多组的待预测日预测结果及对应的实际负荷数据。
[0071]
步骤3：利用多组的待预测日预测结果及对应的实际负荷数据，采用粒子群算法进行修正方法的关键参数的整定。
[0072]
步骤3中所述的修正算法的关键参数为修正时段结束时间点t
end
、修正启动偏差门槛值door和衰减系数k中的一个或多个。本实施例中需要进行整定的修正参数包括修正时段结束时间点t
end
和修正启动偏差门槛值door和修正函数的衰减系数k。根据累积的多日的一次预测结果和实际负荷数据，以平均修正次数最少和平均最终精度最高为目标，进行参数整定，确定修正启动偏差门槛值和修正时段，得到修正方案。
[0073]
粒子群算法(pso)模拟鸟集群飞行觅食的行为,通过鸟之间的集体协作使群体达到目的。用位置、速度和适应度值来表征粒子，粒子在迭代中不断更新自己的位置和速度。通过为目标参数设置合适的变化范围，用粒子群算法解决多目标优化问题得到帕累托最优解，即在修正次数最少的情况下获得最高的预测精度，并将结果作为修正算法的内置参数。
[0074]
由于影响最终精度的关键问题在修正启动条件和修正方法，即修正问题主要集中在三个参数：修正时段结束时间点t
end
和修正启动偏差门槛值door和修正函数的衰减系数k。初始化一群随机粒子(随机解)，然后通过迭代找到最优解。在每一次迭代中，粒子通过跟踪两个“极值”来更新自己。在找到这两个最优值pbest和gbest后，粒子通过下面的公式来更新自己的速度和位置。
[0075]vi，j
(t 1)＝ωv
i，j
(t) c1*rand(0，1)*[pbest
i，j-x
i，j
(t)] c2*rand(0，1)*[gbest
i，j-x
i，j
(t)]
[0076]
x
i，j
(t 1)＝x
i，j
(t) v
i，j
(t 1)
[0077]vi，j
是粒子的速度；pbest和gbest分别定义为粒子最佳位置和种群最佳位置；rand(0，1)是介于(0，1)之间的随机数；x
i，j
(t)是粒子的当前位置；c1和c2是学习因子，ω为压缩因子。学习因子c1和c2，压缩因子ω为设置的常数，可根据需要进行调整，一般有常用的取值范围。
[0078]
本实施例中粒子群算法的参数设置为惯性因子0.5，局部速度因子0.1，全局速度因子0.5，粒子群数量50。
[0079]
本实施例中等待训练的三个参数，修正时段结束时间点t
end
、修正启动偏差门槛值door和衰减系数k的变化范围分别设置为[64，72]、[20，40]、[0.005，3]。这个范围主要是考虑到该园区为工业园区，每日负荷变化规律具有一定周期性，且观察多日预测结果可以发现，较大偏差大多集中在96个时段点的第64到第72点之间，也即下午16点到18点之间，这个时间段接近于下班时间，负荷会出现大幅下降，因此出现误差的可能性也会增大。由于该园区负荷峰值在200千瓦到600千瓦之间，因此比较合适的修正启动门槛在10％-20％之间，因此设置为20到40。
[0080]
参数设定完毕之后，将多个工作日和休息日的预测结果和实际负荷的数据集作为
训练库。每轮计算的适应值即为整个数据集的平均修正次数和平均精度，在完成设定轮数计算后，在所有适应值集合中找到帕累托最优解，其对应的输入变量组即为最优参数组。该逻辑如图3所示。
[0081]
基于粒子群算法ps0整定三个修正参数的具体步骤如下：
[0082]
步骤31，将累积的多组的待预测日预测结果及对应的实际负荷数据组成的数据集投入粒子群算法；该数据集即为按一定时间间隔的(例如，根据负荷预测相关标准定义，数据间隔为15分钟，则每日有96个点数据)每日日前预测结果数据和实际负荷数据组成的多日数据集。例如将2019年7-8月日前预测结果和实际负荷数据投入粒子群算法，参数优化得到的修正算法用于2019年9-10月的修正参数，或者用于2020年7-8月的修正参数。
[0083]
步骤32，初始化参数，设置粒子群算法的参数设置为惯性因子0.5，局部速度因子0.1，全局速度因子0.5，粒子群数量50；
[0084]
步骤33，随机初始化若干个粒子，粒子的三维坐标为修正时段结束时间点t
end
、修正启动偏差门槛值door和衰减系数k，设置修正时段结束时间点t
end
、修正启动偏差门槛值door和衰减系数k的变化范围；本实施例中初始化50个粒子，修正时段结束时间点t
end
、修正启动偏差门槛值door和衰减系数k的变化范围分别设置为[64，72]、[20，40]、[0.005，3]；粒子的个数可根据需求设置，修正时间段[64，72]可根据采样点个数具体设置，如每日设置96个采样点，选取[64，72]采样点对应的负荷数据；启动门槛值也可根据具体的负荷大小进行调整，如[20，40]指的是门槛值为20-40kw之间，衰减系数[0.005，3]可根据经验和实验进行调节；
[0085]
步骤34，对每个粒子按数据集的平均修正次数和平均精度计算适应度，寻找全局最优解；
[0086]
步骤35，检查是否满足终止条件，例如达到一定迭代次数则终止，如果不满足终止条件则进入步骤36，满足终止条件则进入步骤37；
[0087]
步骤36，基于以下公式进行速度和位置更新，更新后进入步骤34；
[0088]vi，j
(t 1)＝ωv
i，j
(t) c1*rand(0，1)*[pbest
i，j-x
i，j
(t)] c2*rand(0，1)*[gbest
i，j-x
i，j
(t)]
[0089]
x
i，j
(t 1)＝x
i，j
(t) v
i，j
(t 1)
[0090]vi，j
是粒子的速度；pbest和gbest分别定义为粒子最佳位置和种群最佳位置；rand(0，1)是介于(0，1)之间的随机数；x
i，j
(t)是粒子的当前位置；c1和c2是学习因子，ω为压缩因子；
[0091]
参数组的三个参数，修正时段结束时间点t
end
、修正启动偏差门槛值door和衰减系数k，可以理解为三维空间中的点的三个坐标，寻找粒子帕累托最优解的最佳位置和种群最佳位置得到最优解，点的最优解的坐标即为最优参数组。
[0092]
步骤37，根据每次迭代得到的最优拟合曲线得到帕累托最优解，输出最优解的修正时段结束时间点t
end
和修正启动偏差门槛值door和修正函数的衰减系数k。每次迭代的最优解拟合成曲线的形式，拐点的位置认为是帕累托最优解。
[0093]
步骤4：在待预测日当日的修正时段内，定期基于日前预测负荷值、实际负荷值和修正启动偏差门槛值door检查是否满足修正条件，当满足修正条件时启动修正，通过修正算法修正该日的日前预测结果。
[0094]
将上述步骤3中参数的整定结果，即修正时段结束时间点t
end
和修正启动偏差门槛值door和修正函数的衰减系数k设置在修正算法中。在待预测日修正时段内进行误差的动态检测，每间隔单位时间间隔就对过去一小时的负荷序列进行检测，若满足修正条件，如当误差高于修正启动偏差门槛值door时，则启动修正算法，更新该日的日前预测结果。不断重复上述动态检测步骤判断是否启动修正，直至修正时段结束，即可得到实时更新的负荷预测值供能量管理系统使用。修正算法与该时段误差大小有关，即在原日前负荷预测结果叠加一随时间衰减项，衰减系数k由前述粒子群算法整定得到。修正算法始终在线监测，直至到达整定的每日时间结束点。
[0095]
所述修正启动偏差门槛值door为一正数，步骤4中所述修正条件基于日前预测负荷值、实际负荷值和修正启动偏差门槛值door进行检查。例如预测负荷值与实际负荷值差的绝对值大于修正启动偏差门槛值door时启动修正。例如，启动修正的时刻前一个时间段内所有采样点的预测负荷值和实际负荷值的偏差的绝对值均大于修正启动偏差门槛值door启动修正。
[0096]
本实施中步骤4中所述修正条件为，当处于修正值时段的预测负荷序列和实际负荷序列满足下面两个条件即可启动修正：
[0097]
abs(real
i-predicti)＞door，i∈(n，n x*m]，ti∈(t
now-xh，t
now
]，
[0098]
(real
i-predicti)(real
i 1-predict
i 1
)＞0，i∈(n，n x*m-1]，ti∈(t
now-xh，t
now
]，
[0099]
其中，door为修正程序参数之一修正启动偏差门槛值，reali为第i个采样点的实际负荷值，predicti为第i个采样点的日前预测负荷值，t
now
为启动修正的时刻，n为t
now
前x个小时采样点数据的编号，m为每个小时内采样点的数目，第i个采样点时间为ti。
[0100]
即启动修正的时刻前一个时间段x小时内所有采样点的预测负荷值和实际负荷值的偏差的绝对值均大于修正启动偏差门槛值door，而且负荷预测值和实际值的偏差的符号相同，则符合修正条件。所述前一个时间段x小时可以根据需要进行调整，本实施例优选设定为1个小时，根据需要也可以设置为半小时或者2小时等。
[0101]
参见图4的实施例，以15分钟数据间隔为例，判断t之前过去一个小时内4个数据点(t，t-1/4h，t-1/2h，t-3/4h)的预测负荷值和实际负荷值的绝对值是否大于修正启动偏差门槛值door，即m为4，如果全部大于修正启动偏差门槛值door、并且(yn-xn)即(预测负荷值-实际负荷值)的符号相同，则根据修正方法进行修正。
[0102]
修正算法为：
[0103]ci
＝real
i-predicti，i∈(n，n x*m]
[0104][0105][0106]
其中，reali为第i个采样点的实际负荷值，predicti为第i个采样点的日前预测负荷值，t
now
表示启动修正的时刻，t
end
为修正时段结束时间点，每日修正时段可以为从凌晨1点到t
end
，n为t
now
前x个小时采样点数据的编号，本实施例中x取值为1，m为每个小时内采样点的数目，本实施例中m取值为4，k为衰减系数，用以控制修正量的衰减速度，predict
i*
为第i个采样点修正后的日前预测负荷值；第i个采样点时间为ti，采样点i∈(n，n x*m]对应的采样时间段为ti∈(t
now-xh，t
now
]，采样点i∈[n x*m，24*m)对应的修正时间段为ti∈[t
now
，
t
end
]。
[0107]
假设全天15分钟采样一次，则全天有96个采样点，当前时刻为下午1点，判断前一个小时的数据是否满足修正条件，即x＝1，m＝4，则此时采样点时间为ti(12，13]，采样点i属于(48，48 1x4]，即采样点编号i＝49，50，51，52，判断这4个采样点的数据是否满足修正启动条件；如满足修正启动条件，则获取t
end
，即从当前时刻开始到修正的具体时间，比如计算获取t
end
为下午2点，即修正时刻从(13，14]，此时对应的采样点[53，54，55，56]，即依据修正公式修正采样点[53，54，55，56]的预测负荷值。以2019年7月11日为例，基于上述修正方法修正后的最终预测曲线与原预测曲线对比如图5所示。
[0108]
本例中精度评估标准采用日平均预测准确率，表达式如下：
[0109][0110]
其中，y
′i表示当日第i点预测值，yi表示当日第i点实际值。
[0111]
如图6所示，连续预测16日的日前预测精度和经过修正算法后的最终精度对比。可见经过修正后的预测结果更加准备，有利于保证微网安全稳定运行，实现微网最优经济运行。
[0112]
步骤5：每隔一定时间重复步骤3，进行修正算法的关键参数的整定。
[0113]
随着微网系统的运行，累积的待预测日预测结果及对应的实际负荷数据越来越多，每隔一定时间重复步骤3，对修正算法的关键参数重新进行整定，更新修正模型。
[0114]
本发明还提供一种微网能量管理设备，包括一个或多个处理器以及一个或多个存储计算机程序的存储器，所述一个或多个处理器，用于执行所述计算机程序，以执行本发明的基于误差修正的负荷日前预测修正方法的步骤。
[0115]
本发明还提供一种存储有计算机程序的存储介质，当所述计算机程序被一个或多个处理器执行时，致使所述一个或多个处理器执行本发明的基于误差修正的负荷日前预测修正方法的步骤。
[0116]
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。