一种基于贝叶斯网络模型的电力负荷峰值预测方法和装置与流程

本发明涉及用电预测技术领域，具体涉及一种基于贝叶斯网络模型的电力负荷峰值预测方法和装置。

背景技术

近年来，随着国民经济的快速增长，居民的生活水平不断提高，电能在终端能源消费中占比也越来越高，电力负荷整体呈现较高的增长速度，尤其负荷峰值的增长速度更是高于电力负荷的增长速度。目前中国负荷峰值不断创新高，高峰时期的供需不平衡不断加深。因此，非常有必要加强对电力负荷峰值的预测研究，以便为电网的运营和交易等提供更加科学有效的决策依据，实现系统安全、可靠及经济运行。现有技术中电力负荷峰值预测方法有多元线性回归模型、人工神经网络、支持向量机和小波变换算法等方法，其中的多元线性回归模型进行电力负荷峰值的预测，根据历史负荷变化的趋势预测电力负荷峰值，但是当外界影响因素变化较大时，预测的电力负荷峰值偏差较大；其中的人工神经网络虽然能够处理海量信息输入，具备自适应学习能力，但是神经网路模型学习速率慢，面对海量数据时模型学习耗时较长且大多结果无法解释其因果关系；其中的支持向量机有难以处理大规模训练样本的缺点；其中的小波变换算法则通常需与其他算法相结合等。综上所述，现有技术中的电力负荷峰值预测方法具有当外界影响因素变化较大时预测结果偏差较大、学习速率慢、学习结果无法解释以及无法处理大规模样本的缺点。

技术实现要素：

为了克服上述现有技术中在外界影响因素变化较大时预测结果偏差较大、学习速率慢、无法处理大规模样本的不足，本发明提供一种基于贝叶斯网络模型的电力负荷峰值预测方法和装置，先以获取的气象数据和数值化处理后的时间数据为聚类特征，对以获取的负荷数据进行聚类，然后根据预先构建的贝叶斯网络模型对电力负荷峰值进行预测，其中的贝叶斯网络模型根据聚类结果构建，本发明在外界影响因素变化较大时预测结果偏差较小，学习速率快，能够处理大规模样本。

为了实现上述发明目的，本发明采取如下技术方案：

一方面，本发明提供一种基于贝叶斯网络模型的电力负荷峰值预测方法，包括：

以获取的气象数据和数值化处理后的时间数据为聚类特征，对获取的负荷数据进行聚类；

根据预先构建的贝叶斯网络模型对电力负荷峰值进行预测；

所述贝叶斯网络模型根据聚类结果构建。

所述贝叶斯网络模型按下述过程构建：

基于聚类结果，采用智能搜索算法和贝叶斯评分准则学习贝叶斯网络结构；

基于聚类结果和贝叶斯网络结构，采用极大似然估计方法学习贝叶斯网络参数。

所述基于聚类结果，采用智能搜索算法和贝叶斯评分准则学习贝叶斯网络结构，包括：

设最佳贝叶斯网络结构与初始贝叶斯网络结构相同，根据贝叶斯评分准则计算预设步数内每一步的贝叶斯网络结构评分；

判断预设步数内第一步以外的每一步的贝叶斯网络结构评分是否均低于预设步内第一步的贝叶斯网络结构评分；若是，则停止搜索，输出预设步数内第一步的贝叶斯网络结构作为最优网络结构；否则，通过智能搜索算法继续搜索最优网络结构。

根据贝叶斯评分准则，按下式计算预设步数内每一步的贝叶斯网络结构评分：

其中，score表示预设步数内每一步的贝叶斯网络结构评分，l表示极大似然估计的对数似然函数值，m表示与概率密度函数有关的节点数量，p表示当前贝叶斯网络结构中节点数量。

所述通过智能搜索算法继续搜索最优网络结构，包括：

将预设步数内评分最高的贝叶斯网络结构作为新一轮的初始贝叶斯网络结构，通过在新一轮的初始贝叶斯网络结构上增加一条边、改变一条边的方向、或删除一条边，生成候选贝叶斯网络结构；

判断候选贝叶斯网络结构是否满足藐视准则；

若是，则将该候选贝叶斯网络结构作为最优网络结构，并用该候选贝叶斯网络结构对应的禁忌对象更新智能禁忌表；

否则，将预设步数内评分次高的贝叶斯网络结构作为新一轮的初始贝叶斯网络结构，重复上述步骤，直至找到最优网络结构。

另一方面，本发明还提供一种基于贝叶斯网络模型的电力负荷峰值预测装置，包括：

聚类模块，用于以获取的气象数据和数值化处理后的时间数据为聚类特征，对获取的负荷数据进行聚类；

预测模块，根据预先构建的贝叶斯网络模型对电力负荷峰值进行预测；

所述贝叶斯网络模型根据聚类结果构建。

所述装置还包括建模模块，所述建模模块包括：

贝叶斯网络结构学习单元，用于基于聚类结果，采用智能搜索算法和贝叶斯评分准则学习贝叶斯网络结构；

贝叶斯网络参数学习单元，用于基于聚类结果和贝叶斯网络结构，采用极大似然估计方法学习贝叶斯网络参数。

贝叶斯网络结构学习单元包括：

计算单元，用于设最佳贝叶斯网络结构与初始贝叶斯网络结构相同，根据贝叶斯评分准则计算预设步数内每一步的贝叶斯网络结构评分；

判断单元，用于判断预设步数内第一步以外的每一步的贝叶斯网络结构评分是否均低于预设步内第一步的贝叶斯网络结构评分；若是，则停止搜索，输出预设步数内第一步的贝叶斯网络结构作为最优网络结构；否则，通过智能搜索算法继续搜索最优网络结构。

所述计算单元按下式计算预设步数内每一步的贝叶斯网络结构评分：

其中，score表示预设步数内每一步的贝叶斯网络结构评分，l表示极大似然估计的对数似然函数值，m表示与概率密度函数有关的节点数量，p表示当前贝叶斯网络结构中节点数量。

所述判断单元具体用于：

将预设步数内评分最高的贝叶斯网络结构作为新一轮的初始贝叶斯网络结构，通过在新一轮的初始贝叶斯网络结构上增加一条边、改变一条边的方向、或删除一条边，生成候选贝叶斯网络结构；

判断候选贝叶斯网络结构是否满足藐视准则；

若是，则将该候选贝叶斯网络结构作为最优网络结构，并用该候选贝叶斯网络结构对应的禁忌对象更新智能禁忌表；

否则，将预设步数内评分次高的贝叶斯网络结构作为新一轮的初始贝叶斯网络结构，重复上述步骤，直至找到最优网络结构。

与最接近的现有技术相比，本发明提供的技术方案具有以下有益效果：

本发明提供的基于贝叶斯网络模型的电力负荷峰值预测方法中，先以获取的气象数据和数值化处理后的时间数据为聚类特征，对获取的负荷数据进行聚类，然后根据预先构建的贝叶斯网络模型对电力负荷峰值进行预测，其中的贝叶斯网络模型根据聚类结果构建，本发明在外界影响因素变化较大时预测结果偏差较小，学习速率快，能够处理大规模样本；

本发明提供的基于贝叶斯网络模型的电力负荷峰值预测装置包括聚类模块和预测模块，聚类模块用于以获取的气象数据和数值化处理后的时间数据为聚类特征，对获取的负荷数据进行聚类，预测模块根据预先构建的贝叶斯网络模型对电力负荷峰值进行预测，其中的贝叶斯网络模型根据聚类结果构建，本发明在外界影响因素变化较大时预测结果偏差较小，学习速率快，能够处理大规模样本；

本发明中的将贝叶斯网络模型与聚类进行结合，采用智能搜索算法和贝叶斯评分准则学习贝叶斯网络结构，并基于聚类结果和贝叶斯网络结构，采用极大似然估计方法学习贝叶斯网络参数，能够有效处理不完整数据，和其他技术结合进行因果分析，使先验知识和数据有效结合，且有效避免对数据的过渡拟合，对环境因素变化较大时依然具有较高的准确度。

附图说明

图1是本发明实施例中基于贝叶斯网络模型的电力负荷峰值预测方法流程图。

具体实施方式

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

实施例1

本发明实施例1提供一种基于贝叶斯网络模型的电力负荷峰值预测方法，具体流程图如图1所示，具体过程如下：

s101：以获取的气象数据和数值化处理后的时间数据为聚类特征，对获取的负荷数据进行聚类；

s102：根据预先构建的贝叶斯网络模型对电力负荷峰值进行预测，其中的贝叶斯网络模型根据s102得到的聚类结果构建。

可以在s101的以获取的气象数据和数值化处理后的时间数据为聚类特征，对获取的负荷数据进行聚类之前先获取气象数据、时间数据和负荷数据，并对时间数据进行数值化处理。

获取的气象数据包括温度、气压、风速和湿度；时间数据包括月份、星期、日和节假日。

上述s101中对包括月份、星期、日和节假日的时间数据具体按下式进行数值化处理：

montht＝[m1,m2,…,mi,…,m12]

weekt＝[a1,a2,…,aj,…,a6]

dayt＝[d1,d2,…,dp,…,d31]

holidayt＝[h1,h2,…,hq,…,h31]

其中，montht为历史第t日所在的月份，mi为montht中的第i天；

weekt为历史第t日所在的星期，aj为weekt中的第j星期；

dayt为历史第t日，dp为dayt中的第p天；

holidayt为历史第t日所在的节假日，hq为holidayt的第q个节假日。

上述s101中，可以采用k-means算法对负荷数据进行聚类。

上述s102中的贝叶斯网络模型具体按下述过程构建：

基于聚类结果，采用智能搜索算法和贝叶斯评分准则学习贝叶斯网络结构；

基于聚类结果和贝叶斯网络结构，采用极大似然估计方法学习贝叶斯网络参数。

上述学习贝叶斯网络结构的具体过程如下：

设最佳贝叶斯网络结构与初始贝叶斯网络结构相同，根据贝叶斯评分准则计算预设步数内每一步的贝叶斯网络结构评分；

判断预设步数内第一步以外的每一步的贝叶斯网络结构评分是否均低于预设步内第一步的贝叶斯网络结构评分；若是，则停止搜索，输出预设步数内第一步的贝叶斯网络结构作为最优网络结构；否则，通过智能搜索算法继续搜索最优网络结构。

根据贝叶斯评分准则，按下式计算预设步数内每一步的贝叶斯网络结构评分：

其中，score表示预设步数内每一步的贝叶斯网络结构评分，l表示极大似然估计的对数似然函数值，m表示与概率密度函数有关的节点数量，p表示当前贝叶斯网络结构中节点数量；

通过智能搜索算法继续搜索最优网络结构，包括：

将预设步数内评分最高的贝叶斯网络结构作为新一轮的初始贝叶斯网络结构，通过在新一轮的初始贝叶斯网络结构上增加一条边、改变一条边的方向、或删除一条边，生成候选贝叶斯网络结构；

判断候选贝叶斯网络结构是否满足藐视准则；

若是，则将该候选贝叶斯网络结构作为最优网络结构，并用该候选贝叶斯网络结构对应的禁忌对象更新智能禁忌表；

否则，将预设步数内评分次高的贝叶斯网络结构作为新一轮的初始贝叶斯网络结构，重复上述步骤，直至找到最优网络结构。

上述的极大似然估计是通过似然函数求得，似然函数可以是样本联合概率密度f(x1,x2,...,xn)的任意正常数倍，可记为l(μ,σ)。未知p元正态分布总体的x～np(μ,σ)，设x(i)＝(xi1,...,xip)'，i＝1,2,...,n，是从p元正态分布总体x抽取的n个简单随机样本，即x(i),...,x(n)独立同分布与总体x，则观测样本阵为：

其中，表示样本的均值，a表示立方差矩阵，l(μ,σ)表示样本x(i),...,x(n)的似然函数，是μ的函数，且b(μ)＞＝0，即b(μ)非负定。

求解得正态分布的极大似然估计

实施例2

基于同一发明构思，本发明实施例2还提供一种基于贝叶斯网络模型的电力负荷峰值预测装置，包括聚类模块和预测模块，下面对上述几个模块的功能进行详细说明：

其中的聚类模块，用于以获取的气象数据和数值化处理后的时间数据为聚类特征，对获取的负荷数据进行聚类；

其中的预测模块，根据预先构建的贝叶斯网络模型对电力负荷峰值进行预测，其中的贝叶斯网络模型根据聚类结果构建。

上述的气象数据包括温度、气压、风速和湿度，时间数据包括月份、星期、日和节假日。

本发明实施例2提供的基于贝叶斯网络模型的电力负荷峰值预测装置还包括建模模块，该建模模块包括：

贝叶斯网络结构学习单元，用于基于聚类结果，采用智能搜索算法和贝叶斯评分准则学习贝叶斯网络结构；

贝叶斯网络参数学习单元，用于基于聚类结果和贝叶斯网络结构，采用极大似然估计方法学习贝叶斯网络参数。

上述的贝叶斯网络结构学习单元具体包括：

计算单元，用于设最佳贝叶斯网络结构与初始贝叶斯网络结构相同，根据贝叶斯评分准则计算预设步数内每一步的贝叶斯网络结构评分；

判断单元，用于判断预设步数内第一步以外的每一步的贝叶斯网络结构评分是否均低于预设步内第一步的贝叶斯网络结构评分；若是，则停止搜索，输出预设步数内第一步的贝叶斯网络结构作为最优网络结构；否则，通过智能搜索算法继续搜索最优网络结构。

计算单元根据贝叶斯评分准则，按下式计算预设步数内每一步的贝叶斯网络结构评分：

其中，score表示预设步数内每一步的贝叶斯网络结构评分，l表示极大似然估计的对数似然函数值，m表示与概率密度函数有关的节点数量，p表示当前贝叶斯网络结构中节点数量；

上述判断单元通过智能搜索算法继续搜索最优网络结构的具体过程如下：

将预设步数内评分最高的贝叶斯网络结构作为新一轮的初始贝叶斯网络结构，通过在新一轮的初始贝叶斯网络结构上增加一条边、改变一条边的方向、或删除一条边，生成候选贝叶斯网络结构；

判断候选贝叶斯网络结构是否满足藐视准则；

若是，则将该候选贝叶斯网络结构作为最优网络结构，并用该候选贝叶斯网络结构对应的禁忌对象更新智能禁忌表；

否则，将预设步数内评分次高的贝叶斯网络结构作为新一轮的初始贝叶斯网络结构，重复上述步骤，直至找到最优网络结构。

为了描述的方便，以上所述装置的各部分以功能分为各种模块或单元分别描述。当然，在实施本申请时可以把各模块或单元的功能在同一个或多个软件或硬件中实现。

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

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

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

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

最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，所属领域的普通技术人员参照上述实施例依然可以对本发明的具体实施方式进行修改或者等同替换，这些未脱离本发明精神和范围的任何修改或者等同替换，均在申请待批的本发明的权利要求保护范围之内。