一种多能源系统多类型负荷联合预测方法及系统

1.本发明涉及综合能源系统技术领域，更具体地，涉及一种多能源系统多类型负荷联合预测方法及系统。


背景技术：

2.近年来，多能源系统由于其能源灵活高效利用的优点成为研究热点。与传统的独立能源系统相比，多能源系统通过能量转换和存储实现能源供需平衡，能源联系更加紧密。因此需要提前准确预测负荷，才能确保多能源系统的高效运行和经济调度。然而，多能源系统的冷热电负荷预测不同于传统的独立能源系统中的负荷预测，由于多能源系统的能量转换，负荷间存在复杂的耦合关系，此外冷热电负荷在不同季节下存在波动，给多能源系统的负荷预测带来挑战。
3.在本发明技术之前，传统的多种能源系统中，负荷影响因素主要包括历史负荷、天气和日期类型。因此，传统的多能源系统的预测技术存在诸多缺陷，具体如下：1)由于多种能源之间的相互转换，导致多能源系统中不同负荷间存在耦合关系，多能源系统负荷受季节变化的影响较大，但现有研究尚无法计及不同季节负荷的耦合关系变化；2)对多能源系统负荷耦合特点的负荷预测较少，现有预测模型的预测精度不足的问题；3)现有预测方法结合历史数据和长短期记忆神经网络算法，利用模型的非线性映射和自学习能力，传递有效的时间信息，取得较好的短期负荷预测效果，但在训练过程中会忽略全局的历史数据信息，存在部分局限性，导致预测模型泛化能力较差。这些原因，导致传统的负荷预测性能较差，预测准确性不足。


技术实现要素：

4.鉴于上述问题，本发明提出了一种多能源系统多类型负荷联合预测方法及系统，该方案通过考虑不同季节的负荷耦合关联，充分挖掘基于多类信息的负荷预测，提高多能源系统负荷预测性能。
5.根据本发明实施例第一方面，提供一种多能源系统多类型负荷联合预测方法。
6.在一个或多个实施例中，优选地，所述一种多能源系统多类型负荷联合预测方法包括：
7.获得历史气象数据和负荷测量数据，在剔除异常时间序列数据和填补缺失数据后，进行对每个数据单独归一化处理，生成时间序列样本集；
8.获取季节的时间划分节点，将所述时间序列样本集按照四季划分，计算每个季节的耦合相关性，生成多能预测模型的输入数据；
9.搭建每个季节的双向长短期记忆神经网络负荷预测模型，输入训练样本进行预测模型学习训练并调整优化模型参数；
10.将测试样本输入所述双向长短期记忆神经网络负荷预测模型，生成冷、热、电负荷预测值。
11.在一个或多个实施例中，优选地，所述获得历史气象数据和负荷测量数据，在剔除异常时间序列数据和补全缺失数据后，对每个数据单独归一化处理，生成时间序列样本集，具体包括：
12.获取历史气象数据和负荷测量数据，其中，所述历史气象数据包括温度、湿度、气压、风速、露点、降水量；所述负荷测量数据包括冷负荷、热负荷、电负荷数据，所述历史气象数据和所述负荷测量数据均为时间序列数据；
13.对所述历史气象数据和所述负荷测量数据采用四分位法剔除异常的时间序列数据，采用最小二乘回归法填充缺失时间序列数据，生成高质量气象数据和负荷数据；
14.对所述高质量气象数据与负荷测量数据进行归一化处理，生成时间序列样本集，其中，所述时间序列样本集的时间跨度为3年。
15.在一个或多个实施例中，优选地，所述获取季节的时间划分节点，将所述时间序列样本集按照四季划分，计算每个季节的耦合相关性，生成多能预测模型的输入数据，具体包括：
16.获取季节的时间划分节点，将所述时间序列样本集按照四季时间节点划分为四个季节样本集；
17.计算四个季节样本中冷负荷、热负荷、电负荷的最大互信系数，衡量每个季节冷热电负荷的耦合相关性；
18.根据负荷间互信息系数的统计值，建立联合特征选择方法；
19.选择所述耦合相关性最大的负荷作为多能预测模型的输入数据；
20.将所述多能预测模型的输入数据划分为训练样本集、验证样本集、测试样本集，其中，所述训练样本集占比80％，所述验证样本集占比10％，所述测试样本集占比10％；
21.其中，所述最大互信系数的计算公式如下：
[0022][0023]
其中，x为比较对象表示所分析的负荷特征数据，y为参考对象表示所参照的负荷特征数据，n为数据样本包含的数量，p为联合概率计算，a表示表示在二维空间中x方向上划分的网格数量，b表示在二维空间中y方向上划分的网格数量，mic(x,y)表示比较对象数据x和参考数据y之间的最大相互信息系数，取值范围在0到1之间。
[0024]
在一个或多个实施例中，优选地，所述搭建每个季节的双向长短期记忆神经网络负荷预测模型，输入训练样本进行预测模型学习训练并调整优化模型参数，具体包括：
[0025]
搭建每个季节对应的双向长短期记忆神经网络负荷预测模型，所述双向长短期记忆神经网络负荷预测模型的超参数包括双向神经网络隐层数量、隐层神经元数量、学习率、l2正则化参数；
[0026]
将所述训练样本集作为输入数据，将对应每类负荷预测样本数据作为输出数据；
[0027]
采用网格法确定预测模型的最优超参数，输入到所述双向长短期记忆神经网络负荷预测模型。
[0028]
在一个或多个实施例中，优选地，所述将测试样本输入所述双向长短期记忆神经网络负荷预测模型，生成冷、热、电负荷预测值，具体包括：
[0029]
将所述测试样本集中的数据加入到所述双向长短期记忆神经网络负荷预测模型中，获得归一化预测值；
[0030]
对所述归一化预测值进行反归一化处理，获得多能源系统的冷、热、电负荷预测值。
[0031]
在一个或多个实施例中，优选地，所述根据负荷间互信息系数的统计值，建立联合特征选择方法，具体包括：
[0032]
计算各负荷特征间的最大互信息系数，统计满足最大互信息系数大于阈值0.3以上的负荷特征数量；
[0033]
当所述负荷特征数量大于等于2，说明冷热电负荷特征间的耦合关联性强，联合选择冷热电负荷共同作为预测模型输入特征；
[0034]
当所述负荷特征数量为1，说明冷热电负荷中有两种负荷特征的耦合关联性强，联合选择耦合关联性强的两种负荷特征作为其预测模型的输入特征，另外一种负荷则单独进行预测；
[0035]
当所述负荷特征数量为0，说明冷热电负荷间耦合关联相互影响较弱，将冷热电负荷分别单独进行预测。
[0036]
在一个或多个实施例中，优选地，所述采用网格法确定预测模型的最优超参数，具体包括：
[0037]
拟定双向长短期记忆神经网络的隐层数量、每层网络的神经元数量、学习率和l2正则化参数的取值范围；
[0038]
通过在取值范围内的循环设定不同的超参数进行负荷预测和交叉验证；
[0039]
选择负荷平均预测误差最小的超参数值作为多能源系统负荷预测模型的最优超参数。
[0040]
根据本发明实施例第二方面，提供一种多能源系统多类型负荷联合预测系统。
[0041]
在一个或多个实施例中，优选地，所述一种多能源系统多类型负荷联合预测系统包括：
[0042]
样本生成模块，用于获得历史气象数据和负荷测量数据，在剔除异常时间序列数据和补全缺失数据后，对每个数据单独归一化处理，生成时间序列样本集；
[0043]
季节划分模块，用于获取季节的时间划分节点，将所述时间序列样本集按照四季划分，计算每个季节的耦合相关性，生成多能预测模型的输入数据；
[0044]
参数设置模块，用于搭建每个季节的双向长短期记忆神经网络负荷预测模型，输入训练样本进行预测模型学习训练并调整优化模型参数；
[0045]
参数预测模块，用于将测试样本输入所述双向长短期记忆神经网络负荷预测模型，生成冷、热、电负荷预测值。
[0046]
根据本发明实施例第三方面，提供一种计算机可读存储介质，其上存储计算机程序指令，所述计算机程序指令在被处理器执行时实现如本发明实施例第一方面中任一项所述的方法。
[0047]
根据本发明实施例第四方面，提供一种电子设备，包括存储器和处理器，所述存储器用于存储一条或多条计算机程序指令，其中，所述一条或多条计算机程序指令被所述处理器执行以实现本发明实施例第一方面中任一项所述的步骤。
[0048]
本发明的实施例提供的技术方案可以包括以下有益效果：
[0049]
1)本发明实施例中，对多能源系统数据进行清洗，剔除负荷异常值，填补负荷缺失值，提高多能源系统的数据质量，为负荷预测提供基础数据。
[0050]
2)本发明实施例中，选用时间跨度包含为包含3年较长的历史实际数据，训练样本包含较长时期的多能源负荷数据，并且包含同期历史多能源系统负荷数据，可适用于不同季节的负荷预测；由于负荷是时间序列数据具有时序性变化规律，可根据历年同期和最近时期历史数据的不同时间尺度的综合数据来综合挖掘未来多能源系统负荷变化规律。
[0051]
3)本发明实施例中，采用最大互信息系数分析方法，分析挖掘多能源系统负荷间的耦合相关性，通过联合特征选择方法，选择负荷预测模型的输入特征，进而提高负荷预测的效率和精度。
[0052]
4)本发明实施例中，基于双向长短期记忆神经网络建立负荷预测模型，通过神经元双向传递有效的时间序列数据信息，根据模型自学习能力充分挖掘多能源系统负荷的全局数据规律，预测精度高且不容易过拟合。
[0053]
本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
[0054]
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
[0055]
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0056]
图1是本发明一个实施例的一种多能源系统多类型负荷联合预测方法的流程图。
[0057]
图2是本发明一个实施例的一种多能源系统多类型负荷联合预测方法中的获得历史气象数据和负荷测量数据，在剔除异常时间序列数据和补全缺失数据后，对每个数据单独归一化处理，生成时间序列样本集的流程图。
[0058]
图3是本发明一个实施例的一种多能源系统多类型负荷联合预测方法中的获取季节的时间划分节点，将所述时间序列样本集按照四季划分，计算每个季节的耦合相关性，生成多能预测模型的输入数据的流程图。
[0059]
图4为多能源系统结能源交互结构图。
[0060]
图5是本发明一个实施例的一种多能源系统多类型负荷联合预测方法中的搭建每个季节的双向长短期记忆神经网络负荷预测模型，输入训练样本进行预测模型学习训练并调整优化模型参数的流程图。
[0061]
图6为双向长短期记忆神经网络结构图。
[0062]
图7是本发明一个实施例的一种多能源系统多类型负荷联合预测方法中的将测试样本输入所述双向长短期记忆神经网络负荷预测模型，生成冷、热、电负荷预测值的流程图。
[0063]
图8是本发明一个实施例的一种多能源系统多类型负荷联合预测方法中的根据负
荷间互信息系数的统计值，建立联合特征选择方法的流程图。
[0064]
图9是本发明一个实施例的一种多能源系统多类型负荷联合预测方法中的采用网格法确定预测模型的最优超参数的流程图。
[0065]
图10是本发明一个实施例的一种多能源系统多类型负荷联合预测系统的结构图。
[0066]
图11是本发明一个实施例的四分位法异常数据剔除结果。
[0067]
图12是本发明一个实施例中的最小二乘回归填补缺失数据结果。
[0068]
图13是本发明一个实施例中的分季节负荷间最大互信息系数图。
[0069]
图14是本发明一个实施例中的四个季节多能源系统负荷预测结果图。
[0070]
图15是本发明一个实施例中一种电子设备的结构图。
具体实施方式
[0071]
在本发明的说明书和权利要求书及上述附图中的描述的一些流程中，包含了按照特定顺序出现的多个操作，但是应该清楚了解，这些操作可以不按照其在本文中出现的顺序来执行或并行执行，操作的序号如101、102等，仅仅是用于区分开各个不同的操作，序号本身不代表任何的执行顺序。另外，这些流程可以包括更多或更少的操作，并且这些操作可以按顺序执行或并行执行。需要说明的是，本文中的“第一”、“第二”等描述，是用于区分不同的消息、设备、模块等，不代表先后顺序，也不限定“第一”和“第二”是不同的类型。
[0072]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0073]
近年来，多能源系统由于其能源灵活高效利用的优点成为研究热点。与传统的独立能源系统相比，多能源系统通过能量转换和存储实现能源供需平衡，能源联系更加紧密。因此需要提前准确预测负荷，才能确保多能源系统的高效运行和经济调度。然而，多能源系统的冷热电负荷预测不同于传统的独立能源系统中的负荷预测，由于多能源系统的能量转换，负荷间存在复杂的耦合关系，此外冷热电负荷在不同季节下存在波动，给多能源系统的负荷预测带来挑战。
[0074]
在本发明技术之前，传统的多种能源系统中，负荷影响因素主要包括历史负荷、天气和日期类型。因此，传统的多能源系统的预测技术存在诸多缺陷，具体如下：1)由于多种能源之间的相互转换，导致多能源系统中不同负荷间存在耦合关系，多能源系统负荷受季节变化的影响较大，但现有研究尚无法计及不同季节负荷的耦合关系变化；2)对多能源系统负荷耦合特点的负荷预测较少，现有预测模型的预测精度不足的问题；3)现有预测方法结合历史数据和长短期记忆神经网络算法，利用模型的非线性映射和自学习能力，传递有效的时间信息，取得较好的短期负荷预测效果，但在训练过程中会忽略全局的历史数据信息，存在部分局限性，导致预测模型泛化能力较差。这些原因，导致传统的负荷预测性能较差，预测准确性不足。
[0075]
本发明实施例中，提供了一种多能源系统多类型负荷联合预测方法及系统。该方案通过考虑不同季节的负荷耦合关联，充分挖掘基于多类信息的负荷预测，提高多能源系统负荷预测性能。
[0076]
根据本发明实施例第一方面，提供一种多能源系统多类型负荷联合预测方法。
[0077]
图1为本发明一个实施例的一种多能源系统多类型负荷联合预测方法的流程图。
[0078]
如图1所述，在一个或多个实施例中，优选地，所述一种多能源系统多类型负荷联合预测方法包括：
[0079]
s101、获得历史气象数据和负荷测量数据，在剔除异常时间序列数据和补全缺失数据后，对每个数据单独归一化处理，生成时间序列样本集；
[0080]
s102、获取季节的时间划分节点，将所述时间序列样本集按照四季划分，计算每个季节的耦合相关性，生成多能预测模型的输入数据；
[0081]
s103、搭建每个季节的双向长短期记忆神经网络负荷预测模型，输入训练样本进行预测模型学习训练并调整优化模型参数；
[0082]
s104、将测试样本输入所述双向长短期记忆神经网络负荷预测模型，生成冷、热、电负荷预测值。
[0083]
图2是本发明一个实施例的一种多能源系统多类型负荷联合预测方法中的获得历史气象数据和负荷测量数据，在剔除异常时间序列数据和补全缺失数据后，对每个数据单独归一化处理，生成时间序列样本集的流程图。
[0084]
如图2所示，在一个或多个实施例中，优选地，所述获得历史气象数据和负荷测量数据，在剔除异常时间序列数据和补全缺失数据后，对每个数据单独归一化处理，生成时间序列样本集，具体包括：
[0085]
s201、获取历史气象数据和负荷测量数据，其中，所述历史气象数据包括温度、湿度、气压、风速、露点、降水量；所述负荷测量数据包括冷负荷、热负荷、电负荷数据，所述历史气象数据和所述负荷测量数据均为时间序列数据；
[0086]
s202、对所述历史气象数据和所述负荷测量数据采用四分位法剔除异常的时间序列数据，采用最小二乘回归法填充缺失时间序列数据，生成高质量气象数据和负荷数据；
[0087]
s203、对所述高质量气象数据与负荷测量数据进行归一化处理，生成时间序列样本集，其中，所述时间序列样本集的时间跨度为3年。
[0088]
在本发明实施例中，具体进行了对于数据的初步的预处理，这个预处理过程中，主要是对于不正常的数据的剔除和对于缺少的数据的补全，但是在此过程中由于不同的数据的基准不同，因此，在最终进行了归一化处理，生成的数据最终是一种时间序列的样本数据。
[0089]
具体的，考虑到不同季节的近期以及去年同期历史数据包含了较为全面的多能源系统负荷变化规律，本实施例选用时间跨度3年较长时期的历史实际多能源负荷数据。除了历史负荷数据，多能源系统负荷的主要影响因素还包含气象和日期信息，因此选用时间和气象等非电气特征的时间序列数据，以增强负荷预测模型在不同季节时期与不同环境条件下的适用性。
[0090]
具体的，通过四分位法对多能源系统数据异常值进行检测的过程如下：依据数据偏离平均水平程度来检测多能源系统的异常数据，统计出样本的上四分位数q1，下四分位数q3，计算四分位位距qr，公式如下：
[0091]
qr＝q
1-q3[0092]
其中，在区间[q3—1.5qr,q1 1.5qr]范围外的数据判定为异常值并剔除，然后采
用最小二乘回归方法填充多能源系统缺失数据，将缺失区间邻近的数据样本作为训练数据，建立回归模型并训练，最后对局部缺失部分以模型回归结果进行填补。
[0093]
通过数据清洗，解决多能源系统因为测量设备出现故障，噪声干扰以及数据传输错误所导致的采集值异常和缺失现象，为多能源系统负荷预测模型提供高质量的基础数据。
[0094]
具体的，进行时间序列样本集生成的过程如下：考虑到负荷预测输入数据量纲不统一以及不同特征数据间的量级差异较大，为消除特征单位和尺度差异的影响，同等看待每种特征，对特征进行归一化处理，其中，数据样本集归一化公式具体如下：
[0095][0096]
式中，x为每种特征样本数据，x
max
和x
min
分别为该特征样本的最大值和最小值，x
′
为归一化后的特征数值。
[0097]
如图3所示，在一个或多个实施例中，优选地，所述获取季节的时间划分节点，将所述时间序列样本集按照四季划分，计算每个季节的耦合相关性，生成多能预测模型的输入数据，具体包括：
[0098]
s301、获取季节的时间划分节点，将所述时间序列样本集按照四季时间节点划分为四个季节样本集；
[0099]
举例说明，考虑到冷热电负荷在不同季节下的变化和波动规律不相同，因此根据季节对时间序列数据进行划分，其中春季为3至5月、夏季为6至8月、秋季为9至11月、冬季为12月至次年2月，按年份顺序排列相同季节数据，获得四个季节的时间序列数据。
[0100]
s302、计算四个季节样本中冷负荷、热负荷、电负荷的最大互信系数，衡量每个季节冷热电负荷的耦合相关性；
[0101]
s303、根据负荷间互信息系数的统计值，建立联合特征选择方法；
[0102]
s304、选择所述耦合相关性最大的负荷作为多能预测模型的输入数据；
[0103]
s305、将所述多能预测模型的输入数据划分为训练样本集、验证样本集、测试样本集，其中，所述训练样本集占比80％，所述验证样本集占比10％，所述测试样本集占比10％；
[0104]
本发明实施例中，考虑到双向长短期记忆神经网络需要进行模型训练和参数调整，对预测总样本集进行划分，输入80％的训练样本集进行负荷预测模型学习训练，根据10％的验证样本集调整预测模型超参数，最后输入10％的测试样本集验证模型的负荷预测效果。
[0105]
其中，所述最大互信系数的计算公式如下：
[0106][0107]
其中，x为比较对象表示所分析的负荷特征数据，y为参考对象表示所参照的负荷特征数据，n为数据样本包含的数量，p为联合概率计算，a表示表示在二维空间中x方向上划分的网格数量，b表示在二维空间中y方向上划分的网格数量，mic(x,y)表示比较对象数据x和参考数据y之间的最大相互信息系数，取值范围在0到1之间。
[0108]
如图4所示为多能源系统结能源交互结构图，分别提供了输入层、中间层和输出层之间的能源交互关系。在本发明实施例中，对于不同的季节分别进行了不同的层级的数据分析，这种分析后，不同类型的负荷信息分别完成了占比对于冷、热、电负荷的耦合相关性的计算，计算后的数据将会分别分配为样本集、验证集和测试集。
[0109]
图5是本发明一个实施例的一种多能源系统多类型负荷联合预测方法中的搭建每个季节的双向长短期记忆神经网络负荷预测模型，输入训练样本进行预测模型学习训练并调整优化模型参数的流程图。
[0110]
如图5所示，在一个或多个实施例中，优选地，所述搭建每个季节的双向长短期记忆神经网络负荷预测模型，输入训练样本进行预测模型学习训练并调整优化模型参数，具体包括：
[0111]
s501、搭建每个季节对应的双向长短期记忆神经网络负荷预测模型，所述双向长短期记忆神经网络负荷预测模型的超参数包括双向神经网络隐层数量、隐层神经元数量、学习率、l2正则化参数；
[0112]
s502、将所述训练样本集作为输入数据，将对应每类负荷预测样本数据作为输出数据；
[0113]
s503、采用网格法确定预测模型的最优超参数，输入到所述双向长短期记忆神经网络负荷预测模型。
[0114]
在本发明实施例中，为了能够分别对各个季节单独预测，搭建每个季节的双向长短期记忆神经网络负荷预测模型，输入训练样本进行预测模型学习训练并调整优化模型参数，其中模型参数是通过优化获得的。具体的，双向长短期记忆神经网络结构图如图6所示。具体的，双向长短期记忆神经网络负荷预测模型，步骤如下：
[0115]
首先，输入门输入前后向当前时刻的状态信息x
(t)
和上一时刻隐藏层的状态信息y
(t-1)
，经过非线性函数变换后，通过遗忘门对状态信息进行筛选，确定需要被遗忘和清楚的冗余状态信息，同时确定有用信息进入新的神经元状态。最后，输出门利用前后向新的神经元状态，建立时间序列前后之间的联系，输出负荷预测最终结果。其中，各变量之间的具体计算公式如下：
[0116]f(t)
＝σ(w
fyy(t-1)
w
fx
x
(t)
bf)
[0117]i(t)
＝σ(w
iyy(t-1)
w
ix
x
(t)
bi)
[0118]u(t)
＝tanh(w
uyy(t-1)
w
ux
x
(t)
bu)
[0119]c(t)
＝i
(t)
gu
(t)
f
t
·c(t-1)
[0120]o(t)
＝σ(w
oyy(t-1)
w
ox
x
(t)
bo)
[0121]y(t)
＝o
(t)
gtanh(c
(t)
)
[0122]
y'
(t)
＝g(w
y'fyf(t)
w
y'byb(t)
b
y'
)
[0123]
其中，wfy,wfx,wix,wiy,wuy,wux,woy,wox分别对应遗忘门f、输入门i和u，输出门o，隐藏层状态y由输出门输出和当前神经元状态c共同决定，预测结果y’由前向yf和后向yb预测结果共同决定，此外公式中的σ、g和tanh为激活函数。
[0124]
为了比较负荷预测效果，建立预测结果评估指标，采用平均绝对百分误差mape与均方根误差rmse两种指标衡量预测值与实际值之间的偏差。具体误差公式如下式所示：
[0125][0126][0127]
其中，n为样本数量，y
′i表示模型的预测值，yi表示实际值。
[0128]
如图7所示，在一个或多个实施例中，优选地，所述将测试样本输入所述双向长短期记忆神经网络负荷预测模型，生成冷、热、电负荷预测值，具体包括：
[0129]
s701、将所述测试样本集中的数据加入到所述双向长短期记忆神经网络负荷预测模型中，获得归一化预测值；
[0130]
具体的，考虑到将测试集样本输入训练好的双向长短期记忆神经网络模型，通过模型所保存的权重和偏置参数进行计算，可以获得模型的预测结果。
[0131]
s702、对所述归一化预测值进行反归一化处理，获得多能源系统的冷、热、电负荷预测值。
[0132]
在本发明实施例中，由于模型是通过归一化处理后的数据训练获得的，因此在执行实际的样本数据的运算过程中，同样需要以归一化的预测值进行输出，但是输出后的数据还需要反馈为各类负荷的实际值，因此在此基础上进行了反归一化处理，获得各类负荷的实际预测值。
[0133]
如图8所示，在一个或多个实施例中，优选地，所述根据负荷间互信息系数的统计值，建立联合特征选择方法，具体包括：
[0134]
s801、计算各负荷特征间的最大互信息系数，统计满足最大互信息系数大于阈值0.3以上的负荷特征数量；
[0135]
s802、当所述负荷特征数量大于等于2，说明冷热电负荷特征间的耦合关联性强，联合选择冷热电负荷共同作为预测模型输入特征；
[0136]
s803、当所述负荷特征数量为1，说明冷热电负荷中有两种负荷特征的耦合关联性强，联合选择耦合关联性强的两种负荷特征作为其预测模型的输入特征，另外一种负荷则单独进行预测；
[0137]
s804、当所述负荷特征数量为0，说明冷热电负荷间耦合关联相互影响较弱，将冷热电负荷分别单独进行预测。
[0138]
在本发明实施例中，对于根据负荷之间互信系数统计值进行了详细的联合特征选择，分别在不同的负荷特征数量下选择不同的预测方式。
[0139]
如图9所示，在一个或多个实施例中，优选地，所述采用网格法确定预测模型的最优超参数，具体包括：
[0140]
s901、拟定双向长短期记忆神经网络的隐层数量、每层网络的神经元数量、学习率和l2正则化参数的取值范围；
[0141]
s902、通过在取值范围内的循环设定不同的超参数进行负荷预测和交叉验证；
[0142]
s903、选择负荷平均预测误差最小的超参数值作为多能源系统负荷预测模型的最优超参数。
[0143]
根据本发明实施例第二方面，提供一种多能源系统多类型负荷联合预测系统。
[0144]
在一个或多个实施例中，优选地，所述一种多能源系统多类型负荷联合预测系统包括：
[0145]
样本生成模块1001，用于获得历史气象数据和负荷测量数据，在剔除异常时间序列数据和补全缺失数据后，对每个数据单独归一化处理，生成时间序列样本集；
[0146]
季节划分模块1002，用于获取季节的时间划分节点，将所述时间序列样本集按照四季划分，计算每个季节的耦合相关性，生成多能预测模型的输入数据；
[0147]
参数设置模块1003，用于搭建每个季节的双向长短期记忆神经网络负荷预测模型，输入训练样本进行预测模型学习训练并调整优化模型参数；
[0148]
参数预测模块1004，用于将测试样本输入所述双向长短期记忆神经网络负荷预测模型，生成冷、热、电负荷预测值。
[0149]
以下提供具体实施例对本发明技术方案作进一步实施解释说明：
[0150]
根据某地区实际多能源系统2017至2019年四年共1095天数据，包含冷热电负荷每30分钟一个点数据，同时还包含多能源系统所处位置主要气象站的温度、湿度、降水量、风速、露点和气压每30分钟一个点的数据，针对该多能源系统数据进行提前半小时的短期负荷预测，预测的具体流程如下：
[0151]
s1101、提升多能源系统的数据质量，本发明实施例对多能源系统数据进行清洗，四分位法异常数据剔除结果如图11所示，最小二乘回归填补缺失数据结果如图12所示。
[0152]
s1102、考虑到去年同期数据，根据季节划分多能源系统数据集，春季276天、夏季276天、秋季273天和冬季270天。
[0153]
s1103、分析多能源系统负荷间的季节耦合关系进行联合特征选择，本发明实施例计算多能源系统冷热电最大互信息系数，分季节负荷间最大互信息系数如图13所示，设置负荷特征选择阈值为0.3，得到多能源系统负荷联合特征选择结果如表1所示。
[0154]
表1
[0155][0156]
其中，春季和秋季冷热电负荷间最大互信息系数大于阈值0.3的统计个数都大于
等于2，说明春秋多能源系统的冷热电负荷间具有较强的耦合关联性，因此联合选择冷热电负荷共同作为负荷预测模型的输入特征。夏季，冷电负荷的最大互信息系数为0.44大于0.3，说明冷电负荷耦合关联性较强，因此冷负荷和电负荷共同选择作为冷电负荷的预测模型输入进行联合预测，而热负荷与其它两类负荷的最大互信息系数都小于0.3，说明热负荷与其它两种负荷的耦合关联性较弱，因此进行单独预测。冬季，冷热负荷的最大互信息系数为0.47，说明冷热负荷耦合关联性较强，因此冷负荷和热负荷共同选作冷热负荷预测模型的输入特征进行联合预测，而电负荷与其它负荷的最大互信息系数较小，耦合关联程度较弱，因此单独进行预测。
[0157]
s1104、本发明实施例采用python代码编写，基于kears框架搭建双向长短期记忆神经网络，设置了双向长短期记忆神经网络激活函数为relu，损失函数为均方误差mse，batch_size为24，其余参数由表2所示，在取值范围内由网格法循环搜索调优确定。
[0158]
表2
[0159][0160]
s1105、通过训练好的双向长短期记忆神经网络负荷预测模型对不同季节的多能源系统负荷进行预测，为了对比本文联合特征选择方法和双向长短期记忆神经网络在预测精度方面的提升，设置4组实验进行比较。实验1：不考虑负荷之间的相关关系，采用长短期记忆神经网络(lstm)建立冷热电负荷单独预测模型。实验2：根据冷热电负荷耦合关系进行输入特征联合选择，采用长短期记忆神经网络(lstm)建立负荷联合预测模型。实验3：不考虑负荷之间的相关关系，采用双向长短期记忆神经网络(bilstm)建立冷热电负荷单独预测模型。实验4：根据冷热电负荷耦合关系进行输入特征联合选择，采用双向长短期记忆神经网络(bilstm)建立负荷联合预测模型。通过上述实验，四个季节多能源系统负荷预测结果如图14所示，每个季节的预测误差如表3至表6所示。
[0161]
本发明实施例从联合特征选择负荷预测结果可知，把耦合关联程度高的负荷联合选择作为预测模型输入预测效果更好。实验2的平均rmse和mape分别为88.52kw和比实验1小0.4％。实验4的平均rmse和mape分别为111.82kw和比实验3小0.5％。特别在秋季，负荷预测精度大幅提高，说明多能源系统秋季负荷之间存在紧密耦合影响的关系。因此验证了组合特征选择在负荷预测中的优越性。为了体现双向长短期记忆神经网络预测模型的优势，对比bilstm和lstm负荷预测模型性能。实验3和实验4每个季节的预测误差明显小于实验1和实验2。此外，从图14中可以看出bilstm的负荷预测曲线更符合实际负荷曲线，这表明bilstm可以更有效地学习时间序列的波动规律，从而提高多能源系统冷热电负荷预测的准确性。
[0162]
表3
[0163][0164]
表4
[0165][0166][0167]
表5
[0168][0169]
表6
[0170][0171]
本发明实施例从多能源系统不同季节的负荷预测结果可知，春、夏、冬三季负荷波动较为规律，预测效果较好。夏季负荷的预测误差最小，而秋季和冬季冷负荷波动较大，导致mape小幅增加，在5％～7％之间。从负荷种类的预测结果可知，电负荷和热负荷预测误差较小，变化相对稳定在5％以内。
[0172]
根据本发明实施例第三方面，提供一种计算机可读存储介质，其上存储计算机程序指令，所述计算机程序指令在被处理器执行时实现如本发明实施例第一方面中任一项所述的方法。
[0173]
根据本发明实施例第四方面，提供一种电子设备。图15是本发明一个实施例中一种电子设备的结构图。图15所示的电子设备为通用负荷联合预测装置，其包括通用的计算
机硬件结构，其至少包括处理器1501和存储器1502。处理器1501和存储器1502通过总线1503连接。存储器1502适于存储处理器1501可执行的指令或程序。处理器1501可以是独立的微处理器，也可以是一个或者多个微处理器集合。由此，处理器1501通过执行存储器1502所存储的指令，从而执行如上所述的本发明实施例的方法流程实现对于数据的处理和对于其它装置的控制。总线1503将上述多个组件连接在一起，同时将上述组件连接到显示控制器1504和显示装置以及输入/输出(i/o)装置1505。输入/输出(i/o)装置1505可以是鼠标、键盘、调制解调器、网络接口、触控输入装置、体感输入装置、打印机以及本领域公知的其他装置。典型地，输入/输出装置1505通过输入/输出(i/o)控制器1506与系统相连。
[0174]
本发明的实施例提供的技术方案可以包括以下有益效果：
[0175]
1)本发明实施例中，对多能源系统数据进行清洗，剔除负荷异常值，填补负荷缺失值，提高多能源系统的数据质量，为负荷预测提供基础数据。
[0176]
2)本发明实施例中，选用时间跨度包含为包含3年较长的历史实际数据，训练样本包含较长时期的多能源负荷数据，并且包含同期历史多能源系统负荷数据，可适用于不同季节的负荷预测；由于负荷是时间序列数据具有时序性变化规律，可根据历年同期和最近时期历史数据的不同时间尺度的综合数据来综合挖掘未来多能源系统负荷变化规律。
[0177]
3)本发明实施例中，采用最大互信息系数分析方法，分析挖掘多能源系统负荷间的耦合相关性，通过联合特征选择方法，选择负荷预测模型的输入特征，进而提高负荷预测的效率和精度。
[0178]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。
[0179]
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0180]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0181]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0182]
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围
之内，则本发明也意图包含这些改动和变型在内。