CN109767043A - 一种电力负荷时间序列大数据智能建模与预测方法 - Google Patents

一种电力负荷时间序列大数据智能建模与预测方法 Download PDF

Info

Publication number
CN109767043A
CN109767043A CN201910043089.2A CN201910043089A CN109767043A CN 109767043 A CN109767043 A CN 109767043A CN 201910043089 A CN201910043089 A CN 201910043089A CN 109767043 A CN109767043 A CN 109767043A
Authority
CN
China
Prior art keywords
time series
frequency
cluster
electric load
component amount
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910043089.2A
Other languages
English (en)
Other versions
CN109767043B (zh
Inventor
刘辉
于程名
吴海平
尹恒鑫
刘泽宇
陈浩林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Xindu Energy Co ltd
Dragon Totem Technology Hefei Co ltd
Original Assignee
Central South University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Central South University filed Critical Central South University
Priority to CN201910043089.2A priority Critical patent/CN109767043B/zh
Publication of CN109767043A publication Critical patent/CN109767043A/zh
Application granted granted Critical
Publication of CN109767043B publication Critical patent/CN109767043B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种电力负荷时间序列大数据智能建模与预测方法,对电力负荷时间序列历史数据进行小波分解,将电力负荷时间序列历史数据分解为电力负荷高频时间序列历史数据和低频时间序列历史数据,然后对所有时间序列进行集成化处理,对时间序列进行聚类,对每种聚类类别的时间序列基于Elman神经网络建立负荷预测模型,最后对预测的分解后的电力负荷进行重构,完成对电力负荷时间序列的智能建模,从而实现对不同特征的电力负荷进行高效、智能预测。

Description

一种电力负荷时间序列大数据智能建模与预测方法
技术领域
本发明属于电力系统领域,具体涉及一种电力负荷时间序列大数据智能建模与预测方法。
背景技术
电力系统在人们生活中发挥中重要的作用,近年来,电力系统如何规划,对电力系统负荷如何调度,来确保人们获得安全、经济、优质的电能受到越来越多的关注。电力负荷大小作为衡量电力系统性能的重要指标之一,对电力负荷进行预测能为优化电力系统提供选择方案,为电力系统进行优质管理起到有效的作用。
传统对电力负荷时间序列进行预测的主要方法是先将非平稳时间序列转换成平稳时间序列,然后使用ARMA模型、TAR模型和ARIMA模型等模型对平稳时间序列进行预测。传统方法建立预测模型较为简单,对电力负荷时间序列的平稳性要求较高,且对非平稳时间序列进行处理时未能寻求关键特征,电力负荷的预测精度和准确度不高。
发明内容
针对上述技术问题,本发明提供一种电力负荷时间序列大数据智能建模与预测方法,实现对不同特征的电力负荷进行高效、智能、精确预测。
为实现上述技术目的,本发明采用如下技术方案:
一种电力负荷时间序列大数据智能建模与预测方法,包括以下步骤:
步骤1,获取历史时间段内的电力负荷值,组成电力负荷时间序列历史数据;
步骤2,对电力负荷时间序列历史数据进行n层小波分解,得到n层电力负荷高频时间序列历史数据和1层电力负荷低频时间序列历史数据;
步骤3,按照时间窗,将n层高频时间序列历史数据和1层低频时间序列历史数据分别截断得到n*Q个高频时间序列分量{xh,i}和Q个低频时间序列分量{xl,i},其中时间窗长度为L1;
其中,xh,i表示第i个高频时间序列分量,xl,i表示第i个低频时间序列分量;
步骤4,按照时间子窗,将所有的高频时间序列分量xh,i和低频时间序列分量xl,i都分为L1/L2个时间子窗,对每个时间子窗内的电力负荷值进行集成化,得到维度均为L1/L2的高频时间序列集成分量{x'h,i}和低频时间序列集成分量{x'l,i};其中,时间子窗长度为L2,L1为L2的整数倍;
步骤5,设聚类簇数为k,设置相似性度量函数,分别对高频时间序列集成分量{x'h,i}和低频时间序列集成分量{x'l,i}进行标准化处理,然后按相似性度量函数进行层次聚类,得到k个高频时间序列聚类类别和相应的k个高频聚类中心,以及k个低频时间序列聚类类别和k个低频聚类中心;
步骤6,分别针对2*k个聚类类别,构建各负荷预测模型的训练样本;
针对每个聚类类别,将其中m个元素作为与聚类类别对应的负荷预测模型的训练样本Y={yj,j=1,2,…,m};其中,yj=[xj,1,xj,2,…xj,l,…,xj,L1/L2]表示训练样本Y中的第j组训练向量,xj,l表示第j组训练向量中的第l个电力负荷集成值,所述元素是指聚类类别下的高频时间序列集成分量或低频时间序列集成分量;
步骤7,分别对2*k个聚类类别,建立2*k个负荷预测模型;
对每个聚类类别建立对应的负荷预测模型的过程如下:
步骤7.1,设局部变量l=1;
步骤7.2,对于每个负荷预测模型,将对应聚类类别的负荷预测模型训练样本Y={yj,j=1,2,…,m}中的每个训练向量yj内的xj,l,xj,l+1,xj,l+2,xj,l+3,xj,l+4为输入、xj,l+5为输出,训练Elman神经网络,判断l=L1/L2-5是否成立:若成立则进入步骤7.3,否则l=l+1,重复步骤7.2;
步骤7.3,针对2*k个聚类类别,每个聚类类别利用各自的训练样本训练Elman神经网络,分别得到与聚类类别对应的2*k个负荷预测模型;
步骤8,使用步骤7.3得到的负荷预测模型对目标时间T的电力负荷进行实时预测,包括以下步骤:
步骤8.1,实时采集当前时刻及近期的电力负荷值,组成电力负荷时间序列实时数据;
步骤8.2,利用步骤2中的小波分解方法,将电力负荷时间序列实时数据分解为n层电力负荷高频时间序列实时数据和1层电力负荷低频时间序列实时数据;
步骤8.3,利用步骤4中的方法,按照时间子窗,将n层电力负荷高频时间序列实时数据和1层电力负荷低频时间序列实时数据进行集成化,得到n层高频时间序列集成分量和1层低频时间序列集成分量;
步骤8.4,利用步骤5设置的相似性度量函数,分别计算n层高频时间序列集成分量与k个高频聚类中心中同时间子窗的部分序列之间的相似性、1层低频时间序列集成分量与k个低频聚类中心中同时间子窗的部分序列之间的相似性,为n个高频时间序列集成分量和1个低频时间序列集成分量分别找到聚类类别,将n个高频时间序列集成分量以及1个低频时间序列集成分量中与目标时间T最近的各自5个电力负荷集成值进行标准化处理,并分别输入到与聚类类别对应的负荷预测模型中,得到n个高频时间序列集成分量以及1个低频时间序列集成分量中各自的下一时刻的电力负荷标准化值,并分别进行反标准化后插入到n个高频时间序列集成分量以及1个低频时间序列集成分量的首位,更新n个高频时间序列集成分量以及1个低频时间序列集成分量;
步骤8.5,判断n个高频时间序列集成分量以及1个低频时间序列集成分量的首位数据所在时刻是否超过目标时间T,若是,进入步骤8.6,否则,返回步骤8.4;
步骤8.6,采用小波重构方法,将n个高频时间序列集成分量以及1个低频时间序列集成分量进行重构,得到电力负荷时间序列预测向量。
进一步地,将小波分解层次n的临界分解层数设为N,聚类簇数k的临界聚类簇数设为K,且首次进行步骤2小波分解的小波分解层次n设为初始值n0,首次进行步骤5聚类簇数k设置的初始值为k0
在步骤6还根据各个聚类类别下的元素构建参数寻优样本:针对每个聚类类别,将各个聚类类别下除去相应的训练样本所剩余的p个元素,作为与聚类类别对应的负荷预测模型的参数调优样本R={re,e=1,2,…,p};其中,re=[re,1,re,2,…,re,l,…,re,L1/L2]表示参数调优样本R中的第e组调优向量,re,l表示第e组调优向量中的第l个电力负荷值;
在步骤7后还包括:步骤C,寻优小波分解层次n和聚类簇数k,得到最优小波分解层次nbest和最优聚类簇数kbest,具体为:
步骤C1,定义局部变量l=1,进入步骤C2;
步骤C2,依次将每个聚类类别的参数调优样本R={re,e=1,2,…,p}中的每个调优向量内的re,l,re,l+1,re,l+2,re,l+3,re,l+4,分别输入到对应聚类类别的负荷预测模型中,判断l=L1/L2-5是否成立,若成立,则进入步骤C3,否则,l=l+1,重复步骤C2;
步骤C3,依据调优向量所属的时间窗,将所有聚类类别的负荷预测模型的输出值组成所属时间窗的长度为l=L1/L2-5的n个高频时间序列预测标准化分量和1个低频时间序列预测标准化分量,分别进行反标准化处理,得到n个高频时间序列预测分量xh,j”和1个低频时间序列预测分量xl”,利用小波重构公式得到电力负荷时间序列预测向量;
步骤C4,利用步骤4中的方法,按照时间子窗,将电力负荷时间序列历史数据进行集成化,得到电力负荷时间序列历史数据集成序列,计算所有的电力负荷时间序列预测向量与电力负荷时间序列历史数据集成序列中对应时间子窗组成的电力负荷时间序列向量之间的均方误差之和,判断均方误差之和是否小于设定误差极限值,若是,则将当前小波分解层次作为最优小波分解层次,将当前聚类簇数作为最优聚类簇数,并得到对应聚类类别的最优聚类中心,建模完成,否则,进入步骤C5;
步骤C5,判断小波分解层次n达到临界分解层数N且聚类簇数k达到临界聚类簇数K是否同时成立,若是,则将历次执行步骤C4所计算的均方误差之和中的最小值所对应的小波分解层次和聚类簇数,分别作为最优小波分解层次和最优聚类簇数,并得到对应聚类类别的最优聚类中心,建模完成,否则,进入步骤C6;
步骤C6,判断聚类簇数k是否达到临界聚类簇数K,若是,令k=k0,n=n+1,然后返回步骤2,否则,令k=k+1,然后返回步骤2。
进一步地,所述步骤4中对每个时间子窗内的电力负荷值进行集成化的具体过程为:
步骤A1,定义局部变量x'0=0;
步骤A2,针对每个时间子窗,分别计算电力负荷集成值x'l其中,wl为第l个时间子窗内的所有电力负荷值,max(wl)为第l个时间子窗内最大的电力负荷值,min(wl)为第l个时间子窗内最小的电力负荷值,α为集成化控制参数;
步骤A3,将x'1,x'2,…,x'L1/L2共L1/L2个时间子窗对应的电力负荷集成值按时间次序组成与高频时间序列分量{xh,i,i=1,2,3,…,n*Q}对应的高频时间序列集成分量{x'h,i,i=1,2,3,…,n*Q}或与低频时间序列分量{xl,i,i=1,2,3,…,Q}对应的低频时间序列集成分量{x'l,i,i=1,2,3,…,Q}。
进一步地,进行小波分解所采用的小波函数为Haar小波函数、dbn小波函数、Mexican Hat小波函数或Meyer小波函数。
进一步地,所述标准化处理的方法包括min-max标准化方法或Z-score标准化方法。
进一步地,所述相似性度量函数选取欧氏距离、马氏距离或明氏距离。
进一步地,聚类中心的获得方法为:分别针对每个聚类类别,将聚类类别内的所有的电力负荷时间序列的平均值作为该聚类类别的聚类中心。
进一步地,L1取值为24小时,L2取值为30分钟。
进一步地,所述步骤8.1中的电力负荷时间序列的时间跨度不小于5*L2,所述目标时间T与步骤8.1中的开始实施预测的当前时刻的时间差小于L1。
有益效果
本发明提供了一种电力负荷时间序列大数据智能建模与预测方法,对电力负荷时间序列历史数据进行小波分解,将电力负荷时间序列历史数据分解为电力负荷高频时间序列历史数据和低频时间序列历史数据,然后对所有时间序列进行集成化处理,对时间序列进行聚类,对每种聚类类别的时间序列基于Elman神经网络建立负荷预测模型,最后对预测的分解后的电力负荷进行重构,完成对电力负荷时间序列的智能建模,从而实现对不同特征的电力负荷进行高效、智能预测。采用小波分解和聚类方法,可以将非平稳的电力负荷时间序列按不同频率进行分类,从而可采用Elman神经网络建立的负荷预测模型对电力负荷时间序列这种非平稳的序列进行准确预测。
同时,在进行小波分解时设置分解层次初始值,在对时间序列进行聚类时设置初始聚类簇数,以建立初始的负荷预测模型;然后利用电力负荷时间序列历史数据构建的参数调优样本寻找最优小波分解层次和最优聚类簇数,使负荷预测模型输出的预测向量在预设误差范围之内,实现负荷预测模型对不同特征的电力负荷进行精确预测。
附图说明
图1为本发明所述的方法流程图。
具体实施方式
下面对本发明的实施例作详细说明,本实施例以本发明的技术方案为依据开展,给出了详细的实施方式和具体的操作过程,对本发明的技术方案作进一步解释说明。
如图1所示,一种电力负荷时间序列大数据智能建模与预测方法,包括以下步骤:
步骤1,获取历史时间段内的电力负荷值,组成电力负荷时间序列历史数据;
所述电力负荷时间序列历史数据包括有功功率时间序列历史数据、无功功率时间序列历史数据;其采样时间间隔设为Δt,数据时间跨度至少为2个月。
步骤2,对电力负荷时间序列历史数据进行n层小波分解,设置小波函数,设小波分解层次为n,初始值n=n0,得到n层电力负荷高频时间序列历史数据和1层电力负荷低频时间序列历史数据;
其中,进行小波分解所采用的小波函数为Haar小波函数、dbn小波函数、MexicanHat小波函数和Meyer小波函数中的一种。
步骤3,设置时间窗长度L1为24小时,以每天时间凌晨零点为时间起点,将n层高频时间序列历史数据和1层低频时间序列历史数据按时间窗分别截断,得到n*Q个高频时间序列分量{xh,i,i=1,2,3,…,n*Q}和Q个低频时间序列分量{xl,i,i=1,2,3,…,Q};其中xh,i表示第i个高频时间序列分量,xl,i表示第i个低频时间序列分量;
步骤4,设置时间子窗长度L2为30分钟,将所有的高频时间序列分量xh,i和低频时间序列分量xl,i都分为48个时间子窗,对每个时间子窗内的电力负荷值进行集成化,得到维度均为48的高频时间序列集成分量{x'h,i,i=1,2,3,…,n*Q}和低频时间序列集成分量{x'l,i,i=1,2,3,…,Q};
其中,对每个时间子窗内的电力负荷值进行集成化的具体过程为:
步骤A1,定义局部变量x'0=0;
步骤A2,针对每个时间子窗,分别计算电力负荷集成值x'l其中,wl为第l个时间子窗内的所有电力负荷值,max(wl)为第l个时间子窗内最大的电力负荷值,min(wl)为第l个时间子窗内最小的电力负荷值,α为集成化控制参数;
步骤A3,将x'1,x'2,…,x'48共48个时间子窗对应的电力负荷集成值按时间次序组成与高频时间序列分量{xh,i,i=1,2,3,…,n*Q}对应的高频时间序列集成分量{x'h,i,i=1,2,3,…,n*Q}或与低频时间序列分量{xl,i,i=1,2,3,…,Q}对应的低频时间序列集成分量{x'l,i,i=1,2,3,…,Q}。
步骤5,设聚类簇数为k,设置相似性度量函数,分别对高频时间序列集成分量{x'h,i,i=1,2,3,…,n*Q}和低频时间序列集成分量{x'l,i,i=1,2,3,…,Q}进行标准化处理,然后进行层次聚类,得到k个高频时间序列聚类类别和相应的k个高频聚类中心,以及k个低频时间序列聚类类别和k个低频聚类中心;
其中,相似性度量函数优先选取欧氏距离、马氏距离或明氏距离;标准化处理的方法为min-max标准化方法或Z-score标准化方法;每个聚类中心的获得方法为:分别针对每个聚类类别,将聚类类别内所有的电力负荷时间序列的平均值作为对应聚类类别的聚类中心。
步骤6,分别针对2*k个聚类类别,构建各负荷预测模型的训练样本和参数寻优样本;
针对每个聚类类别,将其中m个元素作为与聚类类别对应的负荷预测模型的训练样本Y={yj,j=1,2,…,m};其中,yj=[xj,1,xj,2,…xj,l,…,xj,48]表示训练样本Y中的第j组长度为48的训练向量,xj,l表示第j组训练向量中的第l个电力负荷集成值,所述元素是指聚类类别下的高频时间序列集成分量或低频时间序列集成分量;
针对每个聚类类别,将各个聚类类别下除去相应的训练样本所剩余的p个元素,作为与聚类类别对应的负荷预测模型的参数调优样本R={re,e=1,2,…,p};其中,re=[re,1,re,2,…,re,l,…,re,48]表示参数调优样本R中的第e组长度为48的调优向量,re,l表示第e组调优向量中的第l个电力负荷集成值。
步骤7,分别对2*k个聚类类别,建立2*k个负荷预测模型;
对每个聚类类别建立对应的负荷预测模型的过程包括:
步骤7.1,设局部变量l=1,进入步骤7.2;
步骤7.2,对于每个负荷预测模型,将对应聚类类别的负荷预测模型训练样本Y={yj,j=1,2,…,m}中的每个训练向量yj内的xj,l,xj,l+1,xj,l+2,xj,l+3,xj,l+4为输入、xj,l+5为输出,训练Elman神经网络,判断l=L1/L2-5=43是否成立:若成立则进入步骤7.3,否则l=l+1,重复步骤7.2;
步骤7.3,针对2*k个聚类类别,每个聚类类别利用各自的训练样本训练Elman神经网络,分别得到与聚类类别对应的2*k个负荷预测模型;
为获得最优的负荷预测模型,实现对不同特征电力负荷进行精确预测,本发明还包括步骤C,寻优小波分解层次n和聚类簇数k,得到最优小波分解层次nbest和最优聚类簇数kbest,具体过程为:
步骤C1,定义局部变量l=1,进入步骤C2;
步骤C2,依次将每个聚类类别的参数调优样本R={re,e=1,2,…,p}中的每个调优向量内的re,l,re,l+1,re,l+2,re,l+3,re,l+4,分别输入到对应聚类类别的负荷预测模型中,判断l=L1/L2-5=43
是否成立:若成立则进入步骤C3,否则l=l+1,重复步骤C2;
步骤C3,依据调优向量所属的时间窗,将所有聚类类别的负荷预测模型的输出值组成所属时间窗的长度为l=L1/L2-5=43的n个高频时间序列预测标准化分量和1个低频时间序列预测标准化分量,并分别进行反标准化处理,得到n个高频时间序列预测分量xh,j”和1个低频时间序列预测分量xl”,再利用小波重构公式得到电力负荷时间序列预测向量;
步骤C4,利用步骤4中的方法,按照长度为L2为30分钟的时间子窗,将电力负荷时间序列历史数据进行集成化,得到电力负荷时间序列历史数据集成序列,计算所有的电力负荷时间序列预测向量与电力负荷时间序列历史数据集成序列中对应时间子窗组成的电力负荷时间序列向量之间的均方误差之和,判断均方误差之和是否小于设定误差极限值,若是,则将当前小波分解层次n作为最优小波分解层次nbest,将当前聚类簇数k作为最优聚类簇数kbest,并得到kbest个聚类中心,建模完成,否则,进入步骤C5;
步骤C5,判断小波分解层次n达到临界分解层数N且聚类簇数k达到临界聚类簇数K是否同时成立,若是,则将历次执行步骤C4所计算的均方误差之和中的最小值所对应的小波分解层次n和聚类簇数k作为最优小波分解层次nbest和最优聚类簇数kbest,并得到kbest个聚类中心,建模完成,否则,进入步骤C6;
步骤C6,判断聚类簇数k是否达到临界聚类簇数K,若是,令k=k0,n=n+1,然后返回步骤2,否则,令k=k+1,然后返回步骤2。
步骤8,使用步骤7.3得到的负荷预测模型对目标时间T的电力负荷进行实时预测,包括以下步骤:
步骤8.1,实时采集当前时刻及近期的电力负荷值,组成时间跨度至少为2.5个小时的电力负荷时间序列实时数据;其中,目标时间T与实施预测的起始时刻的时间差小于24个小时;
步骤8.2,利用步骤2中的小波分解方法,将电力负荷时间序列分解为nbest层电力负荷高频时间序列实时数据和1层电力负荷低频时间序列实时数据;
步骤8.3,利用步骤4中的方法,以30分钟为一个时间子窗,将nbest层电力负荷高频时间序列实时数据和1层电力负荷低频时间序列实时数据进行集成化,得到nbest层高频时间序列集成分量和1层低频时间序列集成分量;
步骤8.4,利用步骤5设置的相似性度量函数,分别计算nbest层高频时间序列集成分量与kbest个高频聚类中心中同时间子窗的部分序列之间的相似性、1层低频时间序列集成分量与kbest个低频聚类中心中同时间子窗的部分序列之间的相似性,为n个高频时间序列集成分量和1个低频时间序列集成分量分别找到聚类类别,将nbest个高频时间序列集成分量以及1个低频时间序列集成分量中与目标时间T最近的各自5个电力负荷集成值进行标准化处理,并分别输入到与聚类类别对应的负荷预测模型中,得到nbest个高频时间序列集成分量以及1个低频时间序列集成分量中各自的下一时刻的电力负荷标准化值,并分别进行反标准化后插入到nbest个高频时间序列集成分量以及1个低频时间序列集成分量的首位,更新nbest个高频时间序列集成分量以及1个低频时间序列集成分量;
步骤8.5,判断nbest个高频时间序列集成分量以及1个低频时间序列集成分量的首位数据所在时刻是否超过目标时间T,若是,进入步骤8.6,否则,返回步骤8.4;
步骤8.6,采用小波重构方法,将nbest个高频时间序列集成分量以及1个低频时间序列集成分量进行重构,得到电力负荷时间序列预测向量。
以上实施例为本申请的优选实施例,本领域的普通技术人员还可以在此基础上进行各种变换或改进,在不脱离本申请总的构思的前提下,这些变换或改进都应当属于本申请要求保护的范围之内。

Claims (9)

1.一种电力负荷时间序列大数据智能建模与预测方法,其特征在于,包括以下步骤:
步骤1,获取历史时间段内的电力负荷值,组成电力负荷时间序列历史数据;
步骤2,对电力负荷时间序列历史数据进行n层小波分解,得到n层电力负荷高频时间序列历史数据和1层电力负荷低频时间序列历史数据;
步骤3,按照时间窗,将n层高频时间序列历史数据和1层低频时间序列历史数据分别截断得到n*Q个高频时间序列分量{xh,i}和Q个低频时间序列分量{xl,i},其中时间窗长度为L1;
其中,xh,i表示第i个高频时间序列分量,xl,i表示第i个低频时间序列分量;
步骤4,按照时间子窗,将所有的高频时间序列分量xh,i和低频时间序列分量xl,i都分为L1/L2个时间子窗,对每个时间子窗内的电力负荷值进行集成化,得到维度均为L1/L2的高频时间序列集成分量{x'h,i}和低频时间序列集成分量{x'l,i};其中,时间子窗长度为L2,L1为L2的整数倍;
步骤5,设聚类簇数为k,设置相似性度量函数,分别对高频时间序列集成分量{x'h,i}和低频时间序列集成分量{x'l,i}进行标准化处理,然后按相似性度量函数进行层次聚类,得到k个高频时间序列聚类类别和相应的k个高频聚类中心,以及k个低频时间序列聚类类别和k个低频聚类中心;
步骤6,分别针对2*k个聚类类别,构建各负荷预测模型的训练样本;
针对每个聚类类别,将其中m个元素作为与聚类类别对应的负荷预测模型的训练样本Y={yj,j=1,2,…,m};其中,yj=[xj,1,xj,2,…xj,l,…,xj,L1/L2]表示训练样本Y中的第j组训练向量,xj,l表示第j组训练向量中的第l个电力负荷集成值,所述元素是指聚类类别下的高频时间序列集成分量或低频时间序列集成分量;
步骤7,分别对2*k个聚类类别,建立2*k个负荷预测模型;
对每个聚类类别建立对应的负荷预测模型的过程如下:
步骤7.1,设局部变量l=1;
步骤7.2,对于每个负荷预测模型,将对应聚类类别的负荷预测模型训练样本Y={yj,j=1,2,…,m}中的每个训练向量yj内的xj,l,xj,l+1,xj,l+2,xj,l+3,xj,l+4为输入、xj,l+5为输出,训练Elman神经网络,判断l=L1/L2-5是否成立:若成立则进入步骤7.3,否则l=l+1,重复步骤7.2;
步骤7.3,针对2*k个聚类类别,每个聚类类别利用各自的训练样本训练Elman神经网络,分别得到与聚类类别对应的2*k个负荷预测模型;
步骤8,使用步骤7.3得到的负荷预测模型对目标时间T的电力负荷进行实时预测,包括以下步骤:
步骤8.1,实时采集当前时刻及近期的电力负荷值,组成电力负荷时间序列实时数据;
步骤8.2,利用步骤2中的小波分解方法,将电力负荷时间序列实时数据分解为n层电力负荷高频时间序列实时数据和1层电力负荷低频时间序列实时数据;
步骤8.3,利用步骤4中的方法,按照时间子窗,将n层电力负荷高频时间序列实时数据和1层电力负荷低频时间序列实时数据进行集成化,得到n层高频时间序列集成分量和1层低频时间序列集成分量;
步骤8.4,利用步骤5设置的相似性度量函数,分别计算n层高频时间序列集成分量与k个高频聚类中心中同时间子窗的部分序列之间的相似性、1层低频时间序列集成分量与k个低频聚类中心中同时间子窗的部分序列之间的相似性,为n个高频时间序列集成分量和1个低频时间序列集成分量分别找到聚类类别,将n个高频时间序列集成分量以及1个低频时间序列集成分量中与目标时间T最近的各自5个电力负荷集成值进行标准化处理,并分别输入到与聚类类别对应的负荷预测模型中,得到n个高频时间序列集成分量以及1个低频时间序列集成分量中各自的下一时刻的电力负荷标准化值,并分别进行反标准化后插入到n个高频时间序列集成分量以及1个低频时间序列集成分量的首位,更新n个高频时间序列集成分量以及1个低频时间序列集成分量;
步骤8.5,判断n个高频时间序列集成分量以及1个低频时间序列集成分量的首位数据所在时刻是否超过目标时间T,若是,进入步骤8.6,否则,返回步骤8.4;
步骤8.6,采用小波重构方法,将n个高频时间序列集成分量以及1个低频时间序列集成分量进行重构,得到电力负荷时间序列预测向量。
2.根据权利要求1所述的方法,其特征在于,将小波分解层次n的临界分解层数设为N,聚类簇数k的临界聚类簇数设为K,且首次进行步骤2小波分解的小波分解层次n设为初始值n0,首次进行步骤5聚类簇数k设置的初始值为k0
在步骤6还根据各个聚类类别下的元素构建参数寻优样本:针对每个聚类类别,将各个聚类类别下除去相应的训练样本所剩余的p个元素,作为与聚类类别对应的负荷预测模型的参数调优样本R={re,e=1,2,…,p};其中,re=[re,1,re,2,…,re,l,…,re,L1/L2]表示参数调优样本R中的第e组调优向量,re,l表示第e组调优向量中的第l个电力负荷值;
在步骤7后还包括:步骤C,寻优小波分解层次n和聚类簇数k,得到最优小波分解层次nbest和最优聚类簇数kbest,具体为:
步骤C1,定义局部变量l=1,进入步骤C2;
步骤C2,依次将每个聚类类别的参数调优样本R={re,e=1,2,…,p}中的每个调优向量内的re,l,re,l+1,re,l+2,re,l+3,re,l+4,分别输入到对应聚类类别的负荷预测模型中,判断l=L1/L2-5是否成立,若成立,则进入步骤C3,否则,l=l+1,重复步骤C2;
步骤C3,依据调优向量所属的时间窗,将所有聚类类别的负荷预测模型的输出值组成所属时间窗的长度为l=L1/L2-5的n个高频时间序列预测标准化分量和1个低频时间序列预测标准化分量,分别进行反标准化处理,得到n个高频时间序列预测分量xh,j”和1个低频时间序列预测分量xl”,利用小波重构公式得到电力负荷时间序列预测向量;
步骤C4,利用步骤4中的方法,按照时间子窗,将电力负荷时间序列历史数据进行集成化,得到电力负荷时间序列历史数据集成序列,计算所有的电力负荷时间序列预测向量与电力负荷时间序列历史数据集成序列中对应时间子窗组成的电力负荷时间序列向量之间的均方误差之和,判断均方误差之和是否小于设定误差极限值,若是,则将当前小波分解层次作为最优小波分解层次,将当前聚类簇数作为最优聚类簇数,并得到对应聚类类别的最优聚类中心,建模完成,否则,进入步骤C5;
步骤C5,判断小波分解层次n达到临界分解层数N且聚类簇数k达到临界聚类簇数K是否同时成立,若是,则将历次执行步骤C4所计算的均方误差之和中的最小值所对应的小波分解层次和聚类簇数,分别作为最优小波分解层次和最优聚类簇数,并得到对应聚类类别的最优聚类中心,建模完成,否则,进入步骤C6;
步骤C6,判断聚类簇数k是否达到临界聚类簇数K,若是,令k=k0,n=n+1,然后返回步骤2,否则,令k=k+1,然后返回步骤2。
3.根据权利要求1所述的方法,其特征在于,所述步骤4中对每个时间子窗内的电力负荷值进行集成化的具体过程为:
步骤A1,定义局部变量x'0=0;
步骤A2,针对每个时间子窗,分别计算电力负荷集成值x'l其中,wl为第l个时间子窗内的所有电力负荷值,max(wl)为第l个时间子窗内最大的电力负荷值,min(wl)为第l个时间子窗内最小的电力负荷值,α为集成化控制参数;
步骤A3,将x'1,x'2,…,x'L1/L2共L1/L2个时间子窗对应的电力负荷集成值按时间次序组成与高频时间序列分量{xh,i,i=1,2,3,…,n*Q}对应的高频时间序列集成分量{x'h,i,i=1,2,3,…,n*Q}或与低频时间序列分量{xl,i,i=1,2,3,…,Q}对应的低频时间序列集成分量{x'l,i,i=1,2,3,…,Q}。
4.根据权利要求1所述的方法,其特征在于,进行小波分解所采用的小波函数为Haar小波函数、dbn小波函数、Mexican Hat小波函数或Meyer小波函数。
5.根据权利要求1所述的方法,其特征在于,所述标准化处理的方法包括min-max标准化方法或Z-score标准化方法。
6.根据权利要求1所述的方法,其特征在于,所述相似性度量函数选取欧氏距离、马氏距离或明氏距离。
7.根据权利要求1所述的方法,其特征在于,聚类中心的获得方法为:分别针对每个聚类类别,将聚类类别内的所有的电力负荷时间序列的平均值作为该聚类类别的聚类中心。
8.根据权利要求1所述的方法,其特征在于,L1取值为24小时,L2取值为30分钟。
9.根据权利要求1所述的方法,其特征在于,所述步骤8.1中的电力负荷时间序列的时间跨度不小于5*L2,所述目标时间T与步骤8.1中的开始实施预测的当前时刻的时间差小于L1。
CN201910043089.2A 2019-01-17 2019-01-17 一种电力负荷时间序列大数据智能建模与预测方法 Active CN109767043B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910043089.2A CN109767043B (zh) 2019-01-17 2019-01-17 一种电力负荷时间序列大数据智能建模与预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910043089.2A CN109767043B (zh) 2019-01-17 2019-01-17 一种电力负荷时间序列大数据智能建模与预测方法

Publications (2)

Publication Number Publication Date
CN109767043A true CN109767043A (zh) 2019-05-17
CN109767043B CN109767043B (zh) 2020-11-24

Family

ID=66454126

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910043089.2A Active CN109767043B (zh) 2019-01-17 2019-01-17 一种电力负荷时间序列大数据智能建模与预测方法

Country Status (1)

Country Link
CN (1) CN109767043B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110222882A (zh) * 2019-05-21 2019-09-10 国家电网公司西南分部 一种电力系统中长期负荷的预测方法和装置
CN111009898A (zh) * 2019-12-13 2020-04-14 深圳供电局有限公司 一种智慧园区多能协同供电方法及其系统、终端设备
CN111222688A (zh) * 2019-12-06 2020-06-02 北京国电通网络技术有限公司 一种商业楼宇的日负荷预测方法
CN113298300A (zh) * 2021-05-14 2021-08-24 上海电机学院 一种基于hp滤波器及重标极差法的负荷序列聚类方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023147877A1 (en) * 2022-02-04 2023-08-10 Telefonaktiebolaget Lm Ericsson (Publ) Adaptive clustering of time series from geographic locations in a communication network

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101592673A (zh) * 2009-02-18 2009-12-02 中南大学 铁路沿线风速预测的方法
RU2009122419A (ru) * 2009-06-15 2010-12-20 Министерство обороны Российской Федерации Государственное образовательное учреждение профессионального высшего образования Военно- Способ и устройство прогнозирования нестационарного временного ряда
CN102338808A (zh) * 2011-08-26 2012-02-01 天津理工大学 一种在线的风电场短期风速混合预测方法
JP5570409B2 (ja) * 2010-12-22 2014-08-13 三菱電機株式会社 データ圧縮装置及びデータ圧縮方法及びデータ圧縮プログラム
CN104573371A (zh) * 2015-01-16 2015-04-29 上海大学 基于最优小波包变换和ar-svm的非平稳风速预测方法
CN106295798A (zh) * 2016-08-29 2017-01-04 江苏省电力试验研究院有限公司 经验模式分解与Elman神经网络组合风电预测方法
CN107102969A (zh) * 2017-04-28 2017-08-29 湘潭大学 一种时间序列数据的预测方法和系统
CN107358294A (zh) * 2017-07-21 2017-11-17 河北工程大学 一种基于Elman神经网络的需水预测方法
CN107563565A (zh) * 2017-09-14 2018-01-09 广西大学 一种考虑气象因素变化的短期光伏分解预测方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101592673A (zh) * 2009-02-18 2009-12-02 中南大学 铁路沿线风速预测的方法
RU2009122419A (ru) * 2009-06-15 2010-12-20 Министерство обороны Российской Федерации Государственное образовательное учреждение профессионального высшего образования Военно- Способ и устройство прогнозирования нестационарного временного ряда
JP5570409B2 (ja) * 2010-12-22 2014-08-13 三菱電機株式会社 データ圧縮装置及びデータ圧縮方法及びデータ圧縮プログラム
CN102338808A (zh) * 2011-08-26 2012-02-01 天津理工大学 一种在线的风电场短期风速混合预测方法
CN104573371A (zh) * 2015-01-16 2015-04-29 上海大学 基于最优小波包变换和ar-svm的非平稳风速预测方法
CN106295798A (zh) * 2016-08-29 2017-01-04 江苏省电力试验研究院有限公司 经验模式分解与Elman神经网络组合风电预测方法
CN107102969A (zh) * 2017-04-28 2017-08-29 湘潭大学 一种时间序列数据的预测方法和系统
CN107358294A (zh) * 2017-07-21 2017-11-17 河北工程大学 一种基于Elman神经网络的需水预测方法
CN107563565A (zh) * 2017-09-14 2018-01-09 广西大学 一种考虑气象因素变化的短期光伏分解预测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘辉等: "《基于小波分析法与滚动式时间序列法的风电场风速短期预测优化算法》", 《中南大学学报(自然科学版)》 *
刘辉等: "《基于小波分析法与神经网络的非平稳风速信号短期预测优化算法》", 《中南大学学报(自然科学版)》 *
胡蓓蓓: "《HILBERT-HUANG变换在高频金融时间序列中的应用》", 《中国优秀硕士学位论文全文数据库 经济与管理科学辑》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110222882A (zh) * 2019-05-21 2019-09-10 国家电网公司西南分部 一种电力系统中长期负荷的预测方法和装置
CN111222688A (zh) * 2019-12-06 2020-06-02 北京国电通网络技术有限公司 一种商业楼宇的日负荷预测方法
CN111222688B (zh) * 2019-12-06 2024-01-12 北京国电通网络技术有限公司 一种商业楼宇的日负荷预测方法
CN111009898A (zh) * 2019-12-13 2020-04-14 深圳供电局有限公司 一种智慧园区多能协同供电方法及其系统、终端设备
CN113298300A (zh) * 2021-05-14 2021-08-24 上海电机学院 一种基于hp滤波器及重标极差法的负荷序列聚类方法

Also Published As

Publication number Publication date
CN109767043B (zh) 2020-11-24

Similar Documents

Publication Publication Date Title
CN109767043A (zh) 一种电力负荷时间序列大数据智能建模与预测方法
CN107766990B (zh) 一种光伏电站发电功率的预测方法
CN111929748B (zh) 一种气象要素预报方法和系统
CN102274022B (zh) 一种基于脑电信号的睡眠状态监测方法
CN106371610A (zh) 一种基于脑电信号的驾驶疲劳的检测方法
CN101178703B (zh) 基于网络分割的故障诊断谱聚类方法
CN102085095A (zh) 一种检测心电图中st段的方法、系统和心电检测仪器
CN109558975A (zh) 一种电力负荷概率密度多种预测结果的集成方法
CN111882114B (zh) 一种短时交通流量预测模型构建方法及预测方法
CN109117992A (zh) 基于wd-la-wrf模型的超短期风功率预测方法
CN111242355A (zh) 一种基于贝叶斯神经网络的光伏概率预测方法及系统
CN111292124A (zh) 一种基于优化组合神经网络的需水预测方法
CN112186761A (zh) 一种基于概率分布的风电功率场景生成方法及系统
CN113673768A (zh) 基于vmd和cnn-lstm的风功率概率预测模型
CN105159443A (zh) 基于PCA和Granger因果的脑网络特征提取方法
CN114298132A (zh) 风电功率预测方法及装置、电子设备
CN106056167A (zh) 一种基于高斯核混合人工蜂群算法的归一化可能性模糊熵聚类方法
Zhang et al. Hybrid model for renewable energy and load forecasting based on data mining and EWT
CN113449920A (zh) 一种风电功率预测方法、系统及计算机可读介质
CN112200346A (zh) 一种天气波动过程划分与匹配的短期风电功率预测方法
CN116646927A (zh) 一种基于分段滤波与纵横聚类的风电功率预测方法
CN116865232A (zh) 一种基于风速误差修正的中长期风电功率预测方法及系统
CN110321389A (zh) 一种基于时间序列聚类的电能表设备类别识别方法
Qin Software reliability prediction model based on PSO and SVM
CN111222688A (zh) 一种商业楼宇的日负荷预测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230828

Address after: 2201, 22nd Floor, Building 8, Yard 15, Jiajia Garden, Fengtai District, Beijing, 100000

Patentee after: Beijing Xindu Energy Co.,Ltd.

Address before: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province

Patentee before: Dragon totem Technology (Hefei) Co.,Ltd.

Effective date of registration: 20230828

Address after: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province

Patentee after: Dragon totem Technology (Hefei) Co.,Ltd.

Address before: Yuelu District City, Hunan province 410083 Changsha Lushan Road No. 932

Patentee before: CENTRAL SOUTH University