CN110212524A - 一种区域电力负荷预测方法 - Google Patents
一种区域电力负荷预测方法 Download PDFInfo
- Publication number
- CN110212524A CN110212524A CN201910510630.6A CN201910510630A CN110212524A CN 110212524 A CN110212524 A CN 110212524A CN 201910510630 A CN201910510630 A CN 201910510630A CN 110212524 A CN110212524 A CN 110212524A
- Authority
- CN
- China
- Prior art keywords
- load
- time series
- group
- prediction model
- shot
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000002596 correlated effect Effects 0.000 claims abstract description 36
- 230000007787 long-term memory Effects 0.000 claims abstract description 24
- 238000012360 testing method Methods 0.000 claims abstract description 14
- 230000005611 electricity Effects 0.000 claims abstract description 11
- 238000007689 inspection Methods 0.000 claims abstract description 10
- 230000000875 corresponding effect Effects 0.000 claims description 15
- 238000012353 t test Methods 0.000 claims description 13
- 238000005457 optimization Methods 0.000 claims description 10
- 230000006403 short-term memory Effects 0.000 claims description 7
- 230000007774 longterm Effects 0.000 claims description 4
- 230000015654 memory Effects 0.000 claims 3
- 230000006870 function Effects 0.000 description 11
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 8
- 230000008569 process Effects 0.000 description 5
- 230000003442 weekly effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000013277 forecasting method Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 229940068517 fruit extracts Drugs 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011068 loading method Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 125000006850 spacer group Chemical group 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请公开一种区域电力负荷预测方法,包括:根据用电区域的历史负荷数据,构造N组不同时间间隔步长的时间序列特征数据;分别对N组时间序列特征数据进行相关性检验,得到M组通过检验的相关序列;所述相关序列具有序列自相关性;将M组所述相关序列输入多序列长短期记忆网络预测模型中,集成所述多序列长短期记忆网络预测模型的输出结果,预测目标时间点的负荷数据。本申请提供的多序列长短期记忆网络预测模型,基于多个相关序列并捕获负荷数据的长距离依赖性,捕捉如相邻时间点相关性、日相关性、周相关性、月相关性之间的联系,建立特征与负荷之间的非线性关系,提高了负荷预测的准确性,同时还可防止负荷数据样本过度拟合。
Description
技术领域
本申请涉及电力系统技术领域,尤其涉及一种区域电力负荷预测方法。
背景技术
作为使电网更智能化的关键环节之一,区域电力负荷预测在电力系统的规划和运行中起着重要作用。诸如专家系统、灰色系统理论和人工神经网络等许多方法被用于负荷预测,然而由于其负荷预测是一种极其不稳定的时间序列模型,这些方法建立的预测模型很难达到理想的效果,因此需探索一种更优化的负荷预测方法。
相关技术中提出由粒子群优化算法优化的最小二乘支持向量机模型、基于支持向量机的分销网络方法,或者使用随机的基于森林的集成系统等负荷预测方法,这些方法虽然已在传统时间序列模型的基础上,提高了函数的拟合能力并提高了样本的训练效率,但没有考虑到时间序列数据之间的相关性,导致负荷预测的准确性仍旧较低。此外,上述预测方法仅建立特征与负荷之间的非线性映射关系,忽略连续负荷数据样本之间的相关关系,并且存在负荷数据样本过度拟合的问题。
发明内容
本申请提供一种区域电力负荷预测方法,以解决负荷预测准确性低的问题。
本申请提供一种区域电力负荷预测方法,包括:
根据用电区域的历史负荷数据,构造N组不同时间间隔步长的时间序列特征数据;
分别对N组时间序列特征数据进行相关性检验,得到M组通过检验的相关序列,M≦N;所述相关序列具有序列自相关性;
将M组所述相关序列输入多序列长短期记忆网络预测模型中,集成所述多序列长短期记忆网络预测模型的输出结果,预测目标时间点的负荷数据;其中,所述多序列长短期记忆网络预测模型包括若干组长短期记忆网络预测模型。
可选地,所述时间序列特征数据为:
Ei={et-Ki+1,et-2×Ki+1,...,et-j×Ki+1,...,et-n×Ki+1}
式中,Ei表示时间序列特征数据,Ki为时间间隔步长,n为时间序列特征数据的维度,et-j×Ki+1为Ei中第j个时间点对应的负荷变量;i=1,2…,N;j=1,2…,n。
可选地,所述分别对N组时间序列特征数据进行相关性检验,包括:
根据所述时间序列特征数据中任意两个负荷变量之间的Spearman相关系数,以及所述时间序列特征数据的维度n,计算t检验统计量T;
判断t检验统计量T是否大于或等于阈值tα,v,所述阈值tα,v是与显著性水平α和自由度v有关的参数;
如果t检验统计量T大于或等于阈值tα,v,则所述时间序列特征数据中任意两个负荷变量之间具备相关性;反之,则所述时间序列特征数据中任意两个负荷变量不相关。
可选地,所述得到M组通过检验的相关序列,包括:
相关性检验完成后,从所述时间序列特征数据中筛选出非相关负荷变量,所述非相关负荷变量与所述时间序列特征数据中的其他负荷变量都不相关;
将所述非相关负荷变量从所述时间序列特征数据中移除,将移除处理后的时间序列特征数据作为所述相关序列;
丢弃N-M组不具有序列自相关性的时间序列特征数据,从而得到M组通过检验的相关序列。
可选地,所述时间序列特征数据中任意两个负荷变量(Xj,Yj)之间的Spearman相关系数rs为:
式中,Rj为负荷变量Xj的秩,Qj为负荷变量Yj的秩。
可选地,所述t检验统计量T通过如下公式计算:
式中,rs为所述时间序列特征数据中任意两个负荷变量(Xj,Yj)之间的Spearman相关系数。
可选地,在将M组所述相关序列输入多序列长短期记忆网络预测模型之前,所述方法还包括:
对所述相关序列进行如下归一化处理:
式中,xnorm为负荷变量的归一化值,x为所述负荷变量的当前值,xmax为所述负荷变量的最大值,xmin为所述负荷变量的最大值。
可选地,按照如下公式预测目标时间点的负荷数据S:
式中,Sk为第k组长短期记忆网络预测模型输出的目标时间点负荷值,Wk为第k组长短期记忆网络预测模型对应的权重系数,L为所述多序列长短期记忆网络预测模型中所包括的长短期记忆网络预测模型的数量。
可选地,所述方法还包括:
对所述长短期记忆网络预测模型的超参数进行贝叶斯优化,所述超参数包括迭代次数、批尺寸和隐含层节点个数。
可选地,所述长短期记忆网络预测模型采用单隐含层结构。
本申请提供一种基于多序列周期性异构特征交互的区域电力负荷预测方法,根据用电区域的历史负荷数据,来构造N组不同时间间隔步长的时间序列特征数据,在进行相关性检验后,得到M组相关序列,所述相关序列是指通过校验的具有序列自相关性的时间序列特征数据,M组相关序列作为多序列长短期记忆网络(Long Short-Term Memory,LSTM)模型的输入,来预测目标时间点的区域用户的用电负荷。本申请是基于Spearman相关系数和t检验,来衡量电力消费者的时间序列相关性,将若干组并行的LSTM预测模型进行集成,构建多序列LSTM预测模型,一组LSTM预测模型对应一种时间序列类型,通过M组相关序列作为输入并捕获负荷数据中的长距离依赖性,以提高负荷预测的准确性。本申请提供的多序列LSTM预测模型并行集成架构,基于M组相关序列并捕获负荷数据的长距离依赖性,捕捉如相邻时间点相关性、日相关性、周相关性、月相关性之间的联系,建立特征与负荷之间的非线性关系,提高了负荷预测的准确性,同时还可防止负荷数据样本过度拟合。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例示出的区域电力负荷预测方法流程图;
图2为本申请实施例示出的基于多序列LSTM的短期负荷预测框架示意图;
图3a为本申请实施例示出的传统LSTM模型的预测值与真实值的对比图;
图3b为本申请实施例示出的由图2所示的负荷预测架构输出的预测值与真实值的对比图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,为本申请提供的一种区域电力负荷预测方法流程图,所述方法包括如下步骤:
步骤S100,根据用电区域的历史负荷数据,构造N组不同时间间隔步长的时间序列特征数据。
在执行步骤S100之前,需要提取电力系统中记录的用户区域的历史负荷(即耗电量)数据,可选择建立用电区域的负荷数据库,并对负荷数据库中存储的历史负荷数据进行预处理,这里所述的预处理包括但不限于异常值处理和空缺值填充等。对数据进行预处理的过程可以参照现有相关技术,本实施例中不再赘述。
步骤S200,分别对N组时间序列特征数据进行相关性检验,得到M组通过检验的相关序列;所述相关序列具有序列自相关性,自相关性度量的是同一序列在两个不同时期之间的相关程度,形象的讲就是度量自己过去的行为对现在的影响。在此步骤中,由于可能丢弃了未通过相关性检验的时间序列特征数据,即丢弃了不具有序列自相关性的时间序列特征数据,因此M≦N。对于给定的用电区域,其电力生产和管理通常具有较为明显的周期性和规律性,若无重大异常事件发生,一般假设各区域的负荷曲线应具有以下属性:
(1)短时间内的负荷曲线波动应该很小,即相邻时间点的负荷具有相关性。
(2)对于不同日期的同一时间点,例如6月12号PM18:00和6月13号PM18:00,其负荷曲线是相似的,即以日为周期的负荷具有相关性。
(3)具有相同周类型的每日负荷曲线是相似的,其中周类型包括周一、周二…周日,例如对于第一周的周三PM18:00和第二周的周三PM18:00,两者的负荷曲线是相似的,即以周为周期的负荷具有相关性。
(4)以月(或者季节)为周期的负荷曲线是相似的,比如6月12号PM18:00和7月12号PM18:00,即以月(或者季节)为周期的负荷具有相关性。
在历史负荷数据中存在多种时间序列相关性,提取相应的多种时间跨度的序列数据并存储在相应数据结构中,从而构成N组时间序列特征数据,其中N为自然数。本实施例中,分别以相邻时间点、日、周和月作为周期为例,构造四种类型的时间序列特征数据。在实际应用中,还可选择其他类型的周期单位作为时间间隔,比如季节、年等,不限于本实施例所述。
具体地,所述时间序列特征数据可表示为:
Ei={et-Ki+1,et-2×Ki+1,...,et-j×Ki+1,...,et-n×Ki+1}
式中,Ei表示时间序列特征数据,Ki为时间间隔步长,n为时间序列特征数据的维度,et-j×Ki+1为Ei中第j个时间点对应的负荷变量;i=1,2…,N;j=1,2…,n。
对于相邻时间点序列,其时间间隔步长为K1,加载预测时间点之前连续间隔K1的n个时间点的负荷变量数据,构成相邻时间点的时间序列特征数据E1,则有:
E1={et-K1+1,et-2×K1+1,...,et-j×K1+1,...,et-n×K1+1}
对于日间隔时间序列,其时间间隔步长为K2,加载预测时间点之前间隔K2的n个日期内相同时间点的负荷变量数据,构成与日有关的时间序列特征数据E2,则有:
E2={et-K2+1,et-2×K2+1,...,et-j×K2+1,...,et-n×K2+1}
对于周间隔时间序列,其时间间隔步长为K3,加载预测时间点之前间隔K3的n个周内相同时间点的负荷变量数据,构成与周有关的时间序列特征数据E3,则有:
E3={et-K3+1,et-2×K3+1,...,et-j×K3+1,...,et-n×K3+1}
对于月间隔时间序列,其时间间隔步长为K4,加载预测时间点之前间隔K4的n个月内相同时间点的负荷变量数据,构成与月有关的时间序列特征数据E4,则有:
E4={et-K4+1,et-2×K4+1,...,et-j×K4+1,...,et-n×K4+1}
进一步地,在步骤S200中,所述分别对N组时间序列特征数据进行相关性检验,包括:
步骤a,根据所述时间序列特征数据中任意两个负荷变量之间的Spearman相关系数,以及所述时间序列特征数据的维度n,计算t检验统计量T。
对于前述四种类型的时间序列特征数据Ei中的任一种,设Xj和Yj为时间序列特征数据中任意两个负荷变量,令Rj为负荷变量Xj的秩,Qj为负荷变量Yj的秩。如果Xj与Yj同步,则Rj和Qj也是同步的,为了定义秩之间的一致性,Spearman相关系数rs定义如下:
其中,
则有:
t检验用于参数统计中的相关性测试,在零假设下,t检验统计量T通过如下公式计算:
步骤b,判断t检验统计量T是否大于或等于阈值tα,v,所述阈值tα,v是与显著性水平α和自由度v有关的参数。显著性水平α和自由度v可以根据实际对相关性的要求进行选取,当显著性水平α和自由度v确定时,即可确定阈值tα,v,这一过程可具体参照现有t检验的相关方法,从而判定两个负荷变量之间是否相关。
步骤c,如果t检验统计量T大于或等于阈值tα,v,则所述时间序列特征数据中任意两个负荷变量之间具备相关性;反之,如果t检验统计量T小于阈值tα,v,则所述时间序列特征数据中任意两个负荷变量不相关。
由于Spearman相关系数消除了影响整体数据的一些干扰因素,适用于大数据量的相关性判定,因此Spearman相关系数适用于计算负荷预测中的时间序列相关性。
如图2所示,为本实施例提供的基于多序列LSTM的短期负荷预测框架,当提取数据(即历史负荷数据)后,构造N组时间序列特征数据,分别对应相邻时间序列、日间隔时间序列、周间隔时间序列以及月间隔时间序列这四种类型,然后按照上述方式进行相关性检验,并从时间序列特征数据中筛选出非相关负荷变量,所述非相关负荷变量与时间序列特征数据中的其他负荷变量都不相关,然后将筛选出的非相关负荷变量从所述时间序列特征数据中移除,将移除处理后的时间序列特征数据作为所述相关序列,以保证相关序列具有自相关性;丢弃N-M组不具有序列自相关性的未通过检验的时间序列特征数据,从而得到M组通过检验的相关序列。
受限于现有负荷预测模型和方法,一般而言,预测下一个目标时间点的负荷值,主要依赖于前几步的数据,即距离目标时间点较近距离的数据,而不考虑更长步数(长距离)前的数据的影响。一般预测问题也仅需要近距离就能取得不错的效果,也比较合理。然而,对比这种周期性的负荷数据,长距离的关系可能更有用处,日、周、月间隔的时间序列相当于人为构造出长距离关系,例如假设每个时间点间隔单位为1小时,那么一日24小时,就间隔了24个时间点距离,一周则间隔的时间点更多,距离更长。
图2所示的多序列负荷预测架构则充分利用了不同的序列相关性,并捕获负荷数据中的长距离依赖性,以提高负荷预测的准确性,长距离依赖是指与离当前时刻间隔较远的历史负荷数据对当前预测具有相关性。该多序列负荷预测架构能根据序列相关性检验结果,提取多个周期的相关序列,自适应选择多个相关序列作为输入。该多序列负荷预测架构中,具体地,包括多个并行的LSTM预测模型,每个LSTM模型对应一种时间序列类型,比如相邻时间序列对应1号LSTM模型,日间隔时间序列对应2号LSTM模型,周间隔时间序列对应3号LSTM模型,月间隔时间序列对应4号LSTM模型,利用多个LSTM模型的递归神经网络,并行捕获多组相关序列中的信息数据,每个单独的LSTM预测模型都采用单隐含层结构,图2所示的预测架构能够防止负荷数据样本过度拟合。
步骤S300,将M组所述相关序列输入多序列长短期记忆网络预测模型中,集成所述多序列长短期记忆网络预测模型的输出结果,预测目标时间点的负荷数据。本实施例中所述的目标时间点为下一时刻,优选为1小时间隔,当然也可根据具体情况进行设定和修改,不过不推荐多步预测。
在本实施例可选的方案中,在步骤S200之后,并且在步骤S300之前,所述方法还包括:
对所述相关序列进行如下归一化处理:
式中,xnorm为负荷变量的归一化值,x为所述负荷变量的当前值,xmax为所述负荷变量的最大值,xmin为所述负荷变量的最大值,通过归一化处理即可将相关序列中的各个负荷变量的值放缩到[0,1]区间内,从而便于LSTM预测模型的运算和处理。
在各个LSTM模型运算完成后,集成方法是利用预设的权重,对各个LSTM模型的输出结果进行线性加权计算,更具体地,按照如下公式预测目标时间点的负荷数据S:
式中,Sk为第k组长短期记忆网络预测模型输出的目标时间点负荷值,L为所述多序列长短期记忆网络预测模型中所包括的长短期记忆网络预测模型的数量,k可表示为如图2中LSTM模型的序号,i=1、2、3或4,即L=4;Wk为第k组长短期记忆网络预测模型对应的权重系数,权重系数Wk可以根据实际应用进行设定和选取。
LSTM模型的超参数对负荷预测准确性的影响很大,这是因为深度学习对于这些超参数很敏感,不合理的超参数会使得LSTM模型预测不准,甚至严重偏离,因此需要选择一个优秀的超参数优化方式,此外,由于本申请时间序列众多,还需要考虑降低预测运算的复杂度,因此在本实施例其他可能的实现方式中,选择对长短期记忆网络预测模型的超参数进行贝叶斯优化,所述的超参数包括迭代次数、批尺寸和隐含层节点个数等。
贝叶斯优化是一种逼近思想,当计算非常复杂、迭代次数较高时能起到很好的效果,多用于超参数确定,它是基于数据使用贝叶斯定理估计目标函数的后验分布,然后再根据分布选择下一个采样的超参数组合。它充分利用了前一个采样点的信息,其优化的工作方式是通过对目标函数形状的学习,并找到使结果向全局最大提升的参数,从而提高负荷预测的准确性和运算效率。基于贝叶斯优化进行超参数优化,与传统网格搜索相比,参迭代次数少,速度快,即便在非凸问题下也已然能够稳定地执行超参数优化。
本实施例贝叶斯优化选取的Acquisition Function如下:
其中,f(X)为超参数的目标函数值,f(X+)为到目前为止最优超参数的目标函数值,μ(x)为高斯过程所得到的目标函数的均值,σ(x)为高斯过程所得到的目标函数的方差,即f(X)的后验分布。ξ为trade-off系数,如果没有该系数,POI函数会倾向于取在X+周围的点,而加入trade-off系数ξ进行权衡,不断选取新的X,使得POI(X)最大,从而达到优化参数的目的。
本实施例给出对EKPC电站的分时负荷数据进行仿真预测的示例,首先提取EKPC电站5年左右的历史负荷数据,检测异常值和缺失值,并未发现缺失值和异常值。将负荷序列数据转化为多个时间跨度的序列数据并存储在相应数据结构中,以便每次提取。由于此数据是分时数据,那么K1=1,K2=24,K3=7×24,K4=30×24,需要注意的是,K4的取值应与真实月间隔对应。
在每次预测下一时刻的负荷数据时,提取相关的4种类型时间序列特征数据,然后进行基于Spearman相关系数和t检验的相关性检验过程,进而确定是否作为多序列LSTM模型框架的输入,以及丢弃不具有自相关性的未通过检验的时间序列特征数据,每个单独的LSTM预测模型都采用单隐含层结构。仿真实验数据表明,绝大部分情况下,检验结果均为有相关性。然后对需要输入多序列LSTM预测模型的相关序列进行归一化处理,最终集成按照W={W1,W2,W3,W4}={0.4,0.2,0.2,0.2}的权重进行线性加权运算,得到下一时刻的负荷数据S=0.4×S1+0.2×S2+0.2×S3+0.2×S4,其中,S1、S2、S3、S4分别为如图2中1号LSTM模型、2号LSTM模型、3号LSTM模型、4号LSTM模型的输出结果。
如图3a和图3b所示,由于本申请采用多个时间跨度的序列特征以及捕捉负荷数据中的长距离依懒性,与传统仅使用单一相邻时间序列(近距离)的LSTM模型相比,传统LSTM模型中预测值与真实值之间的均方根误差为159.76,而本申请预测架构中预测值与真实值之间的均方根误差为121.46,预测准确性具有较为明显的提升,并且还在一定程度上缓解了常见LSTM模型的滞后错位现象。
本示例进一步针对隐含层节点个数、批尺寸(batch_size)和迭代次数(epoch)进行贝叶斯优化超参数优化,假定参数域服从离散均匀分布,则最终确定LSTM超参数如下:
相邻时间序列对应的LSTM模型:隐含层节点个数为64,batch_size为8192,epoch为64;
日间隔时间序列对应的LSTM模型:隐含层节点个数为64,batch_size为32,epoch为32;
周间隔时间序列对应的LSTM模型:隐含层节点个数为32,batch_size为16,epoch为64;
月间隔时间序列对应的LSTM模型:隐含层节点个数为16,batch_size为2,epoch为128。
由以上技术方案可知,根据用电区域的历史负荷数据,来构造N组不同时间间隔步长的时间序列特征数据,在进行相关性检验后,得到M组相关序列作为多序列长短期记忆网模型的输入,来预测目标时间点的区域用户的用电负荷。本申请是基于Spearman相关系数和t检验,来衡量电力消费者的时间序列相关性,将若干组(本实施例中示出的为四组)并行的LSTM预测模型进行集成,构建多序列LSTM预测模型,一组LSTM预测模型对应一种时间序列类型,通过M组相关序列作为输入并捕获负荷数据中的长距离依赖性,以提高负荷预测的准确性。本申请提供多序列LSTM预测模型并行集成架构,基于M组相关序列并捕获负荷数据的长距离依赖性,捕捉如相邻时间点相关性、日相关性、周相关性、月相关性之间的联系,建立特征与负荷之间的非线性关系,提高了负荷预测的准确性,同时还可防止负荷数据样本过度拟合。
本领域技术人员在考虑说明书及实践这里公开的发明后,将容易想到本发明的其它实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本发明的真正范围和精神由所附的权利要求指出。
应当理解的是,本发明并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。本发明的范围仅由所附的权利要求来限制。
Claims (10)
1.一种区域电力负荷预测方法,其特征在于,包括:
根据用电区域的历史负荷数据,构造N组不同时间间隔步长的时间序列特征数据;
分别对N组时间序列特征数据进行相关性检验,得到M组通过检验的相关序列,M≦N;所述相关序列具有序列自相关性;
将M组所述相关序列输入多序列长短期记忆网络预测模型中,集成所述多序列长短期记忆网络预测模型的输出结果,预测目标时间点的负荷数据;其中,所述多序列长短期记忆网络预测模型包括若干组长短期记忆网络预测模型。
2.根据权利要求1所述的方法,其特征在于,所述时间序列特征数据为:
Ei={et-Ki+1,et-2×Ki+1,...,et-j×Ki+1,...,et-n×Ki+1}
式中,Ei表示时间序列特征数据,Ki为时间间隔步长,n为时间序列特征数据的维度,et-j×Ki+1为Ei中第j个时间点对应的负荷变量;i=1,2…,N;j=1,2…,n。
3.根据权利要求2所述的方法,其特征在于,所述分别对N组时间序列特征数据进行相关性检验,包括:
根据所述时间序列特征数据中任意两个负荷变量之间的Spearman相关系数,以及所述时间序列特征数据的维度n,计算t检验统计量T;
判断t检验统计量T是否大于或等于阈值tα,v,所述阈值tα,v是与显著性水平α和自由度v有关的参数;
如果t检验统计量T大于或等于阈值tα,v,则所述时间序列特征数据中任意两个负荷变量之间具备相关性;反之,则所述时间序列特征数据中任意两个负荷变量不相关。
4.根据权利要求1-3任一项所述的方法,其特征在于,所述得到M组通过检验的相关序列,包括:
相关性检验完成后,从所述时间序列特征数据中筛选出非相关负荷变量,所述非相关负荷变量与所述时间序列特征数据中的其他负荷变量都不相关;
将所述非相关负荷变量从所述时间序列特征数据中移除,将移除处理后的时间序列特征数据作为所述相关序列;
丢弃N-M组不具有序列自相关性的时间序列特征数据,从而得到M组通过检验的相关序列。
5.根据权利要求3所述的方法,其特征在于,所述时间序列特征数据中任意两个负荷变量(Xj,Yj)之间的Spearman相关系数rs为:
式中,Rj为负荷变量Xj的秩,Qj为负荷变量Yj的秩。
6.根据权利要求5所述的方法,其特征在于,所述t检验统计量T通过如下公式计算:
式中,rs为所述时间序列特征数据中任意两个负荷变量(Xj,Yj)之间的Spearman相关系数。
7.根据权利要求1所述的方法,其特征在于,在将M组所述相关序列输入多序列长短期记忆网络预测模型之前,所述方法还包括:
对所述相关序列进行如下归一化处理:
式中,xnorm为负荷变量的归一化值,x为所述负荷变量的当前值,xmax为所述负荷变量的最大值,xmin为所述负荷变量的最大值。
8.根据权利要求1所述的方法,其特征在于,按照如下公式预测目标时间点的负荷数据S:
式中,Sk为第k组长短期记忆网络预测模型输出的目标时间点负荷值,Wk为第k组长短期记忆网络预测模型对应的权重系数,L为所述多序列长短期记忆网络预测模型中所包括的长短期记忆网络预测模型的数量。
9.根据权利要求1所述的方法,其特征在于,所述方法还包括:
对所述长短期记忆网络预测模型的超参数进行贝叶斯优化,所述超参数包括迭代次数、批尺寸和隐含层节点个数。
10.根据权利要求1或9所述的方法,其特征在于,所述长短期记忆网络预测模型采用单隐含层结构。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910510630.6A CN110212524A (zh) | 2019-06-13 | 2019-06-13 | 一种区域电力负荷预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910510630.6A CN110212524A (zh) | 2019-06-13 | 2019-06-13 | 一种区域电力负荷预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110212524A true CN110212524A (zh) | 2019-09-06 |
Family
ID=67792521
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910510630.6A Pending CN110212524A (zh) | 2019-06-13 | 2019-06-13 | 一种区域电力负荷预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110212524A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110765363A (zh) * | 2019-09-27 | 2020-02-07 | 复旦大学 | 一种基于高斯分布表示的深度推荐模型 |
CN110796282A (zh) * | 2019-09-18 | 2020-02-14 | 国网江苏省电力有限公司 | 一种基于相关性系数的行业波动周期监测方法 |
CN111260206A (zh) * | 2020-01-14 | 2020-06-09 | 中国计量大学 | 一种光伏发电影响因素评价模型、构建方法及应用 |
CN112418623A (zh) * | 2020-11-12 | 2021-02-26 | 国网河南省电力公司郑州供电公司 | 基于双向长短时记忆网络和滑动窗输入的反窃电识别方法 |
CN112651537A (zh) * | 2019-10-10 | 2021-04-13 | 国网河北省电力有限公司 | 一种光伏发电超短期功率预测方法及系统 |
CN113158572A (zh) * | 2021-04-27 | 2021-07-23 | 上海交通大学 | 一种短期负荷的预测方法及装置 |
CN113162118A (zh) * | 2021-05-18 | 2021-07-23 | 山东建筑大学 | 一种风力发电机组海上低压穿越检测方法 |
CN113837503A (zh) * | 2021-11-29 | 2021-12-24 | 新风光电子科技股份有限公司 | 基于人工智能的电网负荷预测方法及系统 |
CN114069617A (zh) * | 2021-11-11 | 2022-02-18 | 广东电网有限责任公司 | 一种馈线负荷的预测方法及装置 |
CN117688505A (zh) * | 2024-02-04 | 2024-03-12 | 河海大学 | 一种植被大范围区域化负异常的预测方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106127360A (zh) * | 2016-06-06 | 2016-11-16 | 国网天津市电力公司 | 一种基于用户特性分析的多模型负荷预测方法 |
CN107392368A (zh) * | 2017-07-17 | 2017-11-24 | 天津大学 | 一种基于气象预报的办公建筑动态热负荷组合预测方法 |
CN109034504A (zh) * | 2018-09-14 | 2018-12-18 | 云南电网有限责任公司 | 一种短期负荷预测模型建立方法及装置 |
CN109376953A (zh) * | 2018-11-21 | 2019-02-22 | 国网冀北电力有限公司 | 一种中长期电量预测方法及系统 |
CN109523057A (zh) * | 2018-10-18 | 2019-03-26 | 国网山东省电力公司经济技术研究院 | 一种考虑经济转型背景的区域电网电力负荷预测方法 |
CN109858700A (zh) * | 2019-02-01 | 2019-06-07 | 华北水利水电大学 | 基于相似样本筛选的bp神经网络供暖系统能耗预测方法 |
-
2019
- 2019-06-13 CN CN201910510630.6A patent/CN110212524A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106127360A (zh) * | 2016-06-06 | 2016-11-16 | 国网天津市电力公司 | 一种基于用户特性分析的多模型负荷预测方法 |
CN107392368A (zh) * | 2017-07-17 | 2017-11-24 | 天津大学 | 一种基于气象预报的办公建筑动态热负荷组合预测方法 |
CN109034504A (zh) * | 2018-09-14 | 2018-12-18 | 云南电网有限责任公司 | 一种短期负荷预测模型建立方法及装置 |
CN109523057A (zh) * | 2018-10-18 | 2019-03-26 | 国网山东省电力公司经济技术研究院 | 一种考虑经济转型背景的区域电网电力负荷预测方法 |
CN109376953A (zh) * | 2018-11-21 | 2019-02-22 | 国网冀北电力有限公司 | 一种中长期电量预测方法及系统 |
CN109858700A (zh) * | 2019-02-01 | 2019-06-07 | 华北水利水电大学 | 基于相似样本筛选的bp神经网络供暖系统能耗预测方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110796282A (zh) * | 2019-09-18 | 2020-02-14 | 国网江苏省电力有限公司 | 一种基于相关性系数的行业波动周期监测方法 |
CN110765363A (zh) * | 2019-09-27 | 2020-02-07 | 复旦大学 | 一种基于高斯分布表示的深度推荐模型 |
CN112651537A (zh) * | 2019-10-10 | 2021-04-13 | 国网河北省电力有限公司 | 一种光伏发电超短期功率预测方法及系统 |
CN111260206A (zh) * | 2020-01-14 | 2020-06-09 | 中国计量大学 | 一种光伏发电影响因素评价模型、构建方法及应用 |
CN112418623A (zh) * | 2020-11-12 | 2021-02-26 | 国网河南省电力公司郑州供电公司 | 基于双向长短时记忆网络和滑动窗输入的反窃电识别方法 |
CN113158572A (zh) * | 2021-04-27 | 2021-07-23 | 上海交通大学 | 一种短期负荷的预测方法及装置 |
CN113162118A (zh) * | 2021-05-18 | 2021-07-23 | 山东建筑大学 | 一种风力发电机组海上低压穿越检测方法 |
CN114069617A (zh) * | 2021-11-11 | 2022-02-18 | 广东电网有限责任公司 | 一种馈线负荷的预测方法及装置 |
CN113837503A (zh) * | 2021-11-29 | 2021-12-24 | 新风光电子科技股份有限公司 | 基于人工智能的电网负荷预测方法及系统 |
CN113837503B (zh) * | 2021-11-29 | 2022-03-29 | 新风光电子科技股份有限公司 | 基于人工智能的电网负荷预测方法及系统 |
CN117688505A (zh) * | 2024-02-04 | 2024-03-12 | 河海大学 | 一种植被大范围区域化负异常的预测方法及系统 |
CN117688505B (zh) * | 2024-02-04 | 2024-04-19 | 河海大学 | 一种植被大范围区域化负异常的预测方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110212524A (zh) | 一种区域电力负荷预测方法 | |
Wang et al. | A seasonal GM (1, 1) model for forecasting the electricity consumption of the primary economic sectors | |
Guo et al. | A deep learning model for short-term power load and probability density forecasting | |
CN113962364B (zh) | 一种基于深度学习的多因素用电负荷预测方法 | |
Abuella et al. | Solar power forecasting using artificial neural networks | |
Zhu et al. | Using a self-adaptive grey fractional weighted model to forecast Jiangsu’s electricity consumption in China | |
Yang et al. | A hybrid forecasting approach applied in wind speed forecasting based on a data processing strategy and an optimized artificial intelligence algorithm | |
Zhao et al. | Multi-step wind speed and power forecasts based on a WRF simulation and an optimized association method | |
Wang et al. | Analysis and application of forecasting models in wind power integration: A review of multi-step-ahead wind speed forecasting models | |
CN106960252A (zh) | 基于长短时记忆神经网络的电力负荷预测方法 | |
CN107622329A (zh) | 基于多时间尺度长短时记忆神经网络的电力负荷预测方法 | |
De Giorgi et al. | Error analysis of short term wind power prediction models | |
Khotanzad et al. | A neuro-fuzzy approach to short-term load forecasting in a price-sensitive environment | |
CN106952181A (zh) | 基于长短时记忆神经网络的电力负荷预测系统 | |
CN102930155B (zh) | 获取电力需求的预警参数的方法及装置 | |
Muzumdar et al. | Designing a robust and accurate model for consumer-centric short-term load forecasting in microgrid environment | |
Liu et al. | Heating load forecasting for combined heat and power plants via strand-based LSTM | |
CN113554466A (zh) | 一种短期用电量预测模型构建方法、预测方法和装置 | |
Moustris et al. | Medium, short and very short-term prognosis of load demand for the Greek Island of Tilos using artificial neural networks and human thermal comfort-discomfort biometeorological data | |
Moon et al. | Interpretable Short‐Term Electrical Load Forecasting Scheme Using Cubist | |
Colak et al. | Multi-time series and-time scale modeling for wind speed and wind power forecasting part I: Statistical methods, very short-term and short-term applications | |
CN111859054B (zh) | 气象卫星数据的处理方法及装置 | |
Alhmoud et al. | Short-term load forecasting for Jordan power system based on NARX-ELMAN neural network and ARMA model | |
Heydari et al. | Mid-term load power forecasting considering environment emission using a hybrid intelligent approach | |
Høverstad et al. | Effects of data cleansing on load prediction algorithms |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190906 |