CN107748934B - 基于改进bp神经网络算法的配网电采暖负荷预测方法 - Google Patents

基于改进bp神经网络算法的配网电采暖负荷预测方法 Download PDF

Info

Publication number
CN107748934B
CN107748934B CN201711029996.9A CN201711029996A CN107748934B CN 107748934 B CN107748934 B CN 107748934B CN 201711029996 A CN201711029996 A CN 201711029996A CN 107748934 B CN107748934 B CN 107748934B
Authority
CN
China
Prior art keywords
electric heating
layer
heating load
neuron
output
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.)
Active
Application number
CN201711029996.9A
Other languages
English (en)
Other versions
CN107748934A (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.)
State Grid Xinjiang Electric Power Co Economic And Technological Research Institute
State Grid Corp of China SGCC
North China Electric Power University
Original Assignee
State Grid Xinjiang Electric Power Co Economic And Technological Research Institute
State Grid Corp of China SGCC
North China Electric Power 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 State Grid Xinjiang Electric Power Co Economic And Technological Research Institute, State Grid Corp of China SGCC, North China Electric Power University filed Critical State Grid Xinjiang Electric Power Co Economic And Technological Research Institute
Priority to CN201711029996.9A priority Critical patent/CN107748934B/zh
Publication of CN107748934A publication Critical patent/CN107748934A/zh
Application granted granted Critical
Publication of CN107748934B publication Critical patent/CN107748934B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315Needs-based resource requirements planning or analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06Electricity, gas or water supply
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

一种基于改进BP神经网络算法的配网电采暖负荷预测方法,包括获取历史数据以及参数设置步骤S1,相关系数计算步骤S2,权系数赋值步骤S3,各层神经元计算输出步骤S4,结果输出判断步骤S5,上下限判断步骤S6,神经元学习误差计算步骤S7,基于学习误差的权系数修正步骤S8,随机修正权系数步骤S9,随机修正偏倚和学习率步骤S10,判断是否到最后一组数据步骤S11和最后的预测步骤S12。本发明能够配网电采暖负荷预测方法提高收敛速度,避免输出平坦区,提高权值改变幅度,考虑了弃风供暖和新建绿色建筑对电采暖推广的影响,使电采暖负荷预测结果更接近实际值。

Description

基于改进BP神经网络算法的配网电采暖负荷预测方法
技术领域
本发明涉及一种负荷预测方法,具体的,涉及一种基于改进BP神经网络算法的配网电采暖负荷预测方法。
背景技术
配网负荷预测,是配电网规划的基础,然而随着大量电采暖负荷的接入,采暖内大幅度增加了配电网的负荷,配电网的规划需及时调整适应大规模电采暖负荷的接入,因此,准确的预测电采暖负荷是当前配电网规划前提和基础。
电采暖负荷预测与历史数据的关系较小,主要受以下两个因素的影响:1)我国推广电采暖的主要目的之一是弃风供暖,大规模消纳多余风电,提高新能源利用效率,电采暖的推广受当地年风力发电量利用率的影响较大;2)随着我国新建建筑节能标准的提高,大量新建建筑采用节能75%设计标准,相对于传统采暖,推广电采暖有助于用户节约采暖费用,提高了用户采用电采暖的积极性。目前,还没有考虑以上两个因素的电采暖负荷预测方法。虽然,BP神经网络法可考虑以上影响因素对电采暖负荷预测的影响,但是目前还没有基于BP神经网络法的电采暖负荷预测方法研究。
目前,BP神经网络方法大量应用于配网负荷预测,步骤为:
1)随机赋值权系数;
2)计算各层输出;
3)判断输出与期望值是否满足要求,若满足转向步骤5;
4)计算各层学习误差,修正权系数,转向步骤2);
5)基于计算的权,预测负荷。
但是,该技术应用到电采暖负荷预测后,存在的问题是:权是随机给定的,大量的反复实验导致收敛过程慢;对于给定的权修正量很小,这样优化是通过沿局部改善的方向逐渐进行调整的,很难达到满意的预测结果;在输出接近边界值0或1的条件下,易出现平坦区,权改变极小,容易使训练过程停顿。
因此,如何优化神经网络算法,提高收敛速度,避免输出平坦区,改善训练效果成为现有技术亟需解决的技术问题。
发明内容
本发明的目的在于提出一种基于改进BP神经网络算法的配网电采暖负荷预测方法,能够克服现有BP神经网络算法的缺点,提高电采暖负荷预测的精度,提高收敛速度,避免输出平坦区,提高权值改变幅度。
为达此目的,本发明采用以下技术方案:
一种基于改进BP神经网络算法的配网电采暖负荷预测方法,包括如下步骤:
获取历史数据以及参数设置步骤S1:获取表示输入的历年规划的新建建筑的建筑面积和历年风力发电量利用率,表示输出的历年电采暖负荷值,所述历年包括同时有风力发电和电采暖负荷的时间段,设有m层的神经网络,所述神经网络共有n个神经元,并从1到n顺序编号,所述神经网络计算的数据组数为Vmax=year-1,year是历史数据的年数,令迭代次数it=1,给定最大迭代次数itmax,设训练组数变量v=1,给定输出误差εy,给定神经元接近上限值的误差ε和接近下限值的误差ε,给定学习率η和各神经元的偏倚θi,i=3,4,...,n,电采暖负荷期望输出为y;
相关系数计算步骤S2,分别利用历年的风力发电量利用率、历年新建建筑面积与历年电采暖负荷值计算相关系数;
权系数赋值步骤S3:将计算得到的相关系数分别赋值给各层神经元对应相应变量的权系数ωij,ωij为位于第k-1层总第j个神经元到位于第k层总第i个神经元的权系数,其中i的序号范围对应从第2层到第m层的神经元,j的序号范围对应从第1层到第m-1层的神经元;
各层神经元计算输出步骤S4,对于第k层总第i个神经元输出
Figure BDA0001446114420000031
利用如下公式计算:/>
Figure BDA0001446114420000032
其中j为第k-1层的总第j个神经元,k取值从第2层到第m层,式中θi是第i个神经元的偏倚;
结果输出判断步骤S5,计算输出层输出值与电采暖负荷期望输出y之差的绝对值,如果该绝对值小于输出误差εy,或迭代次数超过itmax,则转向步骤S11,否则进入步骤S6;
上下限判断步骤S6,通过将神经元
Figure BDA0001446114420000033
输出值与神经元接近上限值的误差ε或接近下限值的误差ε进行比较判断神经元/>
Figure BDA0001446114420000034
输出值是否接近其上下限值,如果是则进入步骤S9,否则进入步骤S7;
神经元学习误差计算步骤S7,计算隐藏和输出层学习误差步骤,若该层为输出层,即k=m,则第i个神经元学习误差
Figure BDA0001446114420000035
y是电采暖负荷期望输出;若该层为隐藏层,即k≠m,则学习误差/>
Figure BDA0001446114420000036
l为第k+1层的神经元的序号;
基于学习误差的权系数修正步骤S8,根据学习误差修正权系数,令it=it+1,
Figure BDA0001446114420000037
然后转向步骤S10;
随机修正权系数步骤S9,令it=it+1,随机产生一个很小的正数ε,令ωij(it)=ωij(it-1)+ε
随机修正偏倚和学习率步骤S10,基于[-0.2,0.2]的均匀分布函数,随机产生一个数,修正偏倚和学习率,然后转向步骤S4;
判断是否到最后一组数据步骤S11,令v=v+1,再将v与Vmax比较,判断是否到最后一组数据,如果若v≤Vmax转向步骤S3,否则进入步骤S12;
预测步骤S12,依据计算的各层权系数、神经元偏倚、学习率,给定的预测年的规划新建建筑的建筑面积和风力发电量利用率,预测电采暖负荷值。
可选的,在步骤S2中,首先对历年的风力发电量利用率、历年新建建筑面积与历年电采暖负荷值进行归一化计算,再进行相关系数的计算。
可选的,在步骤S2中,归一化的公式为:
p=(q-0.7qmin)/(1.3qmax-0.7qmin)
其中,p是归一化后的值,q是归一化前的值,qmax,qmin是归一化变量q的最大值、最小值。
可选的,给定输出误差εy为0.01,给定神经元接近上限值的误差ε和接近下限值的误差ε均为0.01。
可选的,步骤S1的各用户能够按采暖时间分成5类,包括学校、行政事业单位、公司企业、商业宾馆和居民。
可选的,在对每一组数据进行训练时的电采暖负荷期望输出y为该组数据中的电采暖负荷值。即采用历年的真实的电采暖负荷值,作为神经网络训练的期望输出值。
可选的,学习率η在[0,1]区间内按均匀分布随机产生。
可选的,在获取历史数据以及参数设置步骤S1中,itmax=500。
可选的,在步骤S9中,ε在[0,0.1]区间按均匀分布条件随机产生。
本发明的配电网电采暖负荷预测方法,能够提高收敛速度,避免输出平坦区,提高权值改变幅度,考虑了弃风供暖和新建绿色建筑对电采暖推广的影响,使电采暖负荷预测结果更接近实际值,有助于电网公司确定针对性的规划方案,不仅保证采暖期内电采暖用户可靠供电,而且能够降低投资,提高电网公司运营的经济效益。
附图说明
图1是根据本发明的基于改进BP神经网络算法的配网电采暖负荷预测方法的流程图;
图2是根据本发明具体实施例的BP神经网络的示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部结构。
参见图1,示出了根据本发明的基于改进BP神经网络算法的配网电采暖负荷预测方法的流程图,该方法包括如下步骤:
获取历史数据以及参数设置步骤S1:获取表示输入的历年规划的新建建筑的建筑面积和历年风力发电量利用率,表示输出的历年电采暖负荷值,所述历年包括同时有风力发电和电采暖负荷的时间段,设有m层的神经网络,所述神经网络共有n个神经元,并从1到n顺序编号,所述神经网络计算的数据组数为Vmax=year-1,year是历史数据的年数,令迭代次数it=1,给定最大迭代次数itmax,设训练组数变量v=1,给定输出误差εy,给定神经元接近上限值的误差ε和接近下限值的误差ε,给定学习率η和各神经元的偏倚θi,i=3,4,...,n,电采暖负荷期望输出为y;
在该步骤中,由于输入变量为两个,即历年规划的新建建筑的建筑面积和历年风力发电量利用率,因此,该神经网络每层的神经元个数为2,输入神经元x1、x2无需具有偏倚值,因此神经元的偏倚从3到n。
相关系数计算步骤S2,分别利用历年的风力发电量利用率、历年新建建筑面积与历年电采暖负荷值计算相关系数。
进一步的,在步骤S2中,首先对历年的风力发电量利用率、历年新建建筑面积与历年电采暖负荷值进行归一化计算。
具体的,归一化公式为:
p=(q-0.7qmin)/(1.3qmax-0.7qmin)
其中,p是归一化后的值,q是归一化前的值,qmax,qmin是归一化变量q的最大值、最小值。
归一化计算为现有技术中常见的归一化计算方法,可以利用如下公式进行归一化计算。
公式为
Figure BDA0001446114420000061
其中,在计算历年风力发电量和历年电采暖负荷的相关系数时,Xi是历年风力发电量数据,Xj是历年电采暖负荷值,Cov表示协方差,Var表示方差;同理,将历年风力发电量数据替换为历年新建建筑面积,可计算历年新建建筑面积和电采暖负荷的相关系数。
权系数赋值步骤S3:将计算得到的相关系数分别赋值给各层神经元对应相应变量的权系数ωij,ωij为位于第k-1层总第j个神经元到位于第k层总第i个神经元的权系数,其中i的序号范围对应从第2层到第m层的神经元,j的序号范围对应从第1层到第m-1层的神经元。
例如ω31表示位于第1层的总第1个神经元到位于第2层的总第3个神经元的权系数,ω53表示位于第2层的总第3个神经元到位于第3层的总第5个神经元的权系数,其中如步骤S1所述,整个神经网络的神经元是从1到n顺序编号的。
各层神经元计算输出步骤S4,对于第k层总第i个神经元输出
Figure BDA0001446114420000071
利用如下公式计算:/>
Figure BDA0001446114420000072
其中j为第k-1层的总第j个神经元,k取值从第2层到第m层,式中θi是第i个神经元的偏倚。
结果输出判断步骤S5,计算输出层输出值与电采暖负荷期望输出y之差的绝对值,如果该绝对值小于输出误差εy,或迭代次数超过itmax,则转向步骤S11,否则进入步骤S6。
可选的,在对每一组数据进行训练时的电采暖负荷期望输出y为该组数据中的电采暖负荷值。即采用历年的真实的电采暖负荷值,作为神经网络训练的期望输出值。
上下限判断步骤S6,通过将神经元
Figure BDA0001446114420000073
输出值与神经元接近上限值的误差ε或接近下限值的误差ε进行比较判断神经元/>
Figure BDA0001446114420000074
输出值是否接近其上下限值,如果是则进入步骤S9,否则进入步骤S7。
例如,若神经元x3或x4或x5的输出值接近其上下限值1和0,即1-xj≤ε1或xj-0≤ε0j=3,4,5,ε1,ε0分别是上下限的误差限值,则转向步骤S9。
神经元学习误差计算步骤S7,计算隐藏和输出层学习误差步骤,若该层为输出层,即k=m,则第i个神经元学习误差
Figure BDA0001446114420000075
y是电采暖负荷期望输出;若该层为隐藏层,即k≠m,则学习误差/>
Figure BDA0001446114420000076
l为第k+1层的神经元的序号;
基于学习误差的权系数修正步骤S8,根据学习误差修正权系数,令it=it+1,
Figure BDA0001446114420000081
然后转向步骤S10;
可选的,对于学习率η能够在[0,1]区间内按均匀分布随机产生。
随机修正权系数步骤S9,令it=it+1,随机产生一个很小的正数ε,令ωij(it)=ωij(it-1)+ε
随机修正偏倚和学习率步骤S10,基于[-0.2,0.2]的均匀分布函数,随机产生一个数,修正偏倚和学习率,然后转向步骤S4;
在该步骤中,该修正公式类似于步骤S9的修正公式。
判断是否到最后一组数据步骤S11,令v=v+1,再将v与Vmax比较,判断是否到最后一组数据,如果若v≤Vmax转向步骤S3,否则进入步骤S12;
预测步骤S12,基于各层的权系数和给定的输入,预测电采暖负荷值。
下面以一具体示例,结合图2的神经网络对本发明进行进一步的说明。
实施例1:
获取历史数据以及参数设置步骤S1:获取供电区域内历年规划的新建建筑的建筑面积,历年的风力发电量利用率,历年的电采暖负荷;设有m=3层神经网络,即一个输入层,一个隐藏层,一个输出层,神经元总数是n=5,因此输入层神经元为x1,x2,隐藏层神经元为x3,x4,输出层神经元为x5,电采暖负荷期望输出为y,神经网络计算的数据组数为Vmax=year-1,year是历史数据的年数,令迭代次数it=1,训练组数变量v=1,给定最大迭代次数itmax,给定输出误差限值εy,给定神经元接近上限值的误差ε和接近下限值的误差ε,给定学习率η和各神经元的偏倚θi,i=3,4,...,n。
相关系数计算步骤S2,分别计算历年的风力发电量利用率与电采暖负荷值的相关系数ξ1,历年的新建建筑面积与电采暖负荷值的相关系数ξ2
权系数赋值步骤S3,将步骤S2计算的相关系数赋值给隐藏层和输出层的权系数,即ω31=ω41=ω53=ξ1,ω32=ω42=ω52=ξ2
各层神经元计算输出步骤S4,对于第k层总第i个神经元输出
Figure BDA0001446114420000093
利用如下公式计算:/>
Figure BDA0001446114420000091
其中j为第k-1层的总第j个神经元,k取值从第2层到第m层,式中θi是第i个神经元的偏倚。
进一步的,能够以隐藏层的神经元x3为例,有
Figure BDA0001446114420000092
神经元x4,x5计算方式与神经元x3的计算方法类似。
结果输出判断步骤S5,若|x5-y|<εy εy为输出误差限值,或者it>ITmax,ITmax为最大迭代次数,则转向步骤S11,否则进入步骤S6。
上下限判断步骤S6,若神经元x3或x4或x5的输出值接近其上下限值1和0,即1-xj≤ε1或xj-0≤ε0 j=3,4,5,ε1,ε0分别是上下限的误差限值,如果是则进入步骤S9,否则进入步骤S7。
神经元学习误差计算步骤S7,计算隐藏和输出层学习误差,对于输出层神经元x5的学习误差d5=x5(1-x5)(y-x5),隐藏层神经元x3,x4的学习误差d3=x3(1-x3)d5ω53,d4=x4(1-x4)d5ω54
基于学习误差的权系数修正步骤S8,根据学习误差修正权系数,令it=it+1,
Figure BDA0001446114420000101
η为学习率,在[0,1]区间内按均匀分布随机产生,然后转向步骤S10;
随机修正权系数步骤S9,令it=it+1,随机产生一个很小的正数ε,令ωij(it)=ωij(it-1)+ε
其中,ε在[0,0.1]区间按均匀分布条件随机产生。
随机修正偏倚和学习率步骤S10,基于[-0.2,0.2]的均匀分布函数随机产生偏倚Δθ3,Δθ4,Δθ5,修正偏倚θ3(it)=θ3(it-1)+Δθ3,θ4(it)=θ4(it-1)+Δθ4,θ5(it)=θ5(it-1)+Δθ5,基于[-0.2,0.2]均匀分布函数随机产生学习率修正值Δη,修正学习率η(it)=η(it-1)+Δη,然后转向步骤S4;
判断是否到最后一组数据步骤S11,v=v+1,判断是否到最后一组数据,如果若v≤Vmax,转向步骤S3,否则进入步骤S12;
预测步骤S12,基于各层的权系数和给定的输入,依据计算的权系数、神经元偏倚、学习率,给定的预测年的规划新建建筑的建筑面积和风力发电量利用率,预测电采暖负荷值。
其中,步骤S1的各用户能够按采暖时间分成5类,包括学校、行政事业单位、公司企业、商业宾馆和居民,以分别训练神经网络并进行计算。
实施例2:
步骤S1:获取步骤
以新疆某市为例,下表给出了2013-2016年该城市的规划新建建筑面积、风力发电量利用率、电采暖负荷值;
Figure BDA0001446114420000102
Figure BDA0001446114420000111
步骤S2:计算相关系数
归一化后的数值如下表所示:
年份 规划新建建筑面积 风力发电量利用率 电采暖负荷
2012 0.0000 0.0000 0.0000
2013 0.1331 0.2736 0.3870
2014 0.2314 0.6700 0.3913
2015 0.4560 0.3841 0.6383
2016 1.0000 1.0000 1.0000
第一组数据为2012年规划新建建筑面积和风力发电量利用率预测2013年的电采暖负荷值,令迭代次数it=1,itmax=500,v=1,输出误差限值εy=0.01,ε1=0.01,ε0=0.01,θ3=0.1756,θ4=0.0472,θ5=0.0953,η=0.542;
计算的新建建筑面积和电采暖负荷的相关系数为ξ1=0.92,风力发电量利用率和电采暖负荷的相关系数为ξ2=0.95;
步骤S3:赋值权系数步骤
ω31=ω4`=ω53=0.92,ω32=ω42=ω54=0.95;
步骤S4:计算各层神经元输出
经计算隐藏层神经元x3=0.544,x4=0.512,输出层神经元x5=0.746;
步骤S5:结果输出判断步骤,若|0.746-0.387|=0.368>0.01,it=1<itmax=500,转向步骤S6;
步骤S6:神经元输出上下限判断判断步骤
神经元x3,x4,x5的计算值与上下限0和1之间的误差均大于0.01,转向步骤S7;
步骤S7:计算隐藏和输出层神经元的学习误差步骤
神经元x5的学习误差d5=-0.068,隐藏层神经元x3,x4的学习误差d3=-0.0155,d4=-0.0161;
步骤S8:修正权系数步骤
计算的迭代步骤it=2,计算的权系数为ω53(it)=0.95,ω54(it)=0.985,ω31(it)=0.927,ω32(it)=0.958,ω41(it)=0.928,ω42(it)=0.959
步骤S9:随机修正权系数步骤
第1组,第1次迭代,不满足要求,算法不经过此步骤;
步骤S10:随机修正学习率和偏倚步骤
随机修正后的偏倚值为θ3=0.0214,θ4=-0.1021,θ5=0.1310,学习率为η=0.249;
步骤S11:将训练组数变量v加1,再判断是否到最后一组数据步骤;
训练组数不到最后一组数据,转向步骤S3;
步骤S12:预测电采暖负荷值步骤
经363次迭代,确定了权系数、学习率和各神经元的偏倚如下:ω53(it)=0.872,ω54(it)=0.854,ω31(it)=0.831,ω32(it)=0.818,ω41(it)=0.814,ω42(it)=0.823,θ3=0.254,θ4=0.153,θ5=0.142,η=0.357。
在给出预测的2017年的规划新建建筑面积167.57万平方,风力发电量利用率为75.36%条件下,下表给出了2013年到2017年,改进的BP设计网络方法和传统的BP神经网络方法电采暖负荷的预测结果对比。
年份 传统BP神经网络法预测值 本文方法预测值 实际电采暖负荷(MW)
2012 - - 7.48
2013 19.58 20.82 21.99
2014 20.11 23.58 22.15
2015 29.41 30.24 31.41
2016 40.23 42.35 44.97
2017 25.32 28.91 -
可见相对于现有技术,本发明的配电网电采暖负荷预测,提高收敛速度,避免输出平坦区,提高权值改变幅度,考虑了弃风供暖和新建绿色建筑对电采暖推广的影响,使电采暖负荷预测结果更接近实际值,有助于电网公司确定针对性的规划方案,不仅保证采暖期内电采暖用户可靠供电,而且能够降低投资,提高电网公司运营的经济效益。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施方式仅限于此,对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单的推演或替换,都应当视为属于本发明由所提交的权利要求书确定保护范围。

Claims (8)

1.一种基于改进BP神经网络算法的配网电采暖负荷预测方法,包括如下步骤:
获取历史数据以及参数设置步骤S1:获取表示输入的历年规划的新建建筑的建筑面积和历年风力发电量利用率,表示输出的历年电采暖负荷值,所述历年包括同时有风力发电和电采暖负荷的时间段,设有m层的神经网络,所述神经网络共有n个神经元,并从1到n顺序编号,所述神经网络计算的数据组数为Vmax=year-1,year是历史数据的年数,令迭代次数it=1,给定最大迭代次数itmax,设训练组数变量v=1,给定输出误差εy,给定神经元接近上限值的误差ε和接近下限值的误差ε,给定学习率η和各神经元的偏倚θi,i=3,4,...,n,电采暖负荷期望输出为y;
相关系数计算步骤S2,分别利用历年的风力发电量利用率、历年新建建筑面积与历年电采暖负荷值计算相关系数;
权系数赋值步骤S3:将计算得到的相关系数分别赋值给各层神经元对应相应变量的权系数ωij,ωij为位于第k-1层总第j个神经元到位于第k层总第i个神经元的权系数,其中i的序号范围对应从第2层到第m层的神经元,j的序号范围对应从第1层到第m-1层的神经元;
各层神经元计算输出步骤S4,对于第k层总第i个神经元输出
Figure FDA0003996815700000011
利用如下公式计算:
Figure FDA0003996815700000012
其中j为第k-1层的总第j个神经元,k取值从第2层到第m层,式中θi是第i个神经元的偏倚;
结果输出判断步骤S5,计算输出层输出值与电采暖负荷期望输出y之差的绝对值,如果该绝对值小于输出误差εy,或迭代次数超过itmax,则转向步骤S11,否则进入步骤S6;
上下限判断步骤S6,通过将神经元
Figure FDA0003996815700000021
输出值与神经元接近上限值的误差ε或接近下限值的误差ε进行比较判断神经元/>
Figure FDA0003996815700000022
输出值是否接近其上下限值,如果是则进入步骤S9,否则进入步骤S7;
神经元学习误差计算步骤S7,计算隐藏和输出层学习误差步骤,若该层为输出层,即k=m,则第i个神经元学习误差
Figure FDA0003996815700000023
y是电采暖负荷期望输出;若该层为隐藏层,即k≠m,则学习误差/>
Figure FDA0003996815700000024
l为第k+1层的神经元的序号;
基于学习误差的权系数修正步骤S8,根据学习误差修正权系数,令it=it+1,
Figure FDA0003996815700000025
然后转向步骤S10;
随机修正权系数步骤S9,令it=it+1,随机产生正数ε,令ωij(it)=ωij(it-1)+ε,在步骤S9中,ε在[0,0.1]区间按均匀分布条件随机产生;
随机修正偏倚和学习率步骤S10,基于[-0.2,0.2]的均匀分布函数,随机产生一个数,修正偏倚和学习率,然后转向步骤S4;
判断是否到最后一组数据步骤S11,令v=v+1,再将v与Vmax比较,判断是否到最后一组数据,如果若v≤Vmax转向步骤S3,否则进入步骤S12;
预测步骤S12,依据计算的各层权系数、神经元偏倚、学习率,给定的预测年的规划新建建筑的建筑面积和风力发电量利用率,预测电采暖负荷值。
2.根据权利要求1所述的配网电采暖负荷预测方法,其特征在于:
在步骤S2中,首先对历年的风力发电量利用率、历年新建建筑面积与历年电采暖负荷值进行归一化计算,再进行相关系数的计算。
3.根据权利要求2所述的配网电采暖负荷预测方法,其特征在于:
在步骤S2中,归一化的公式为:
p=(q-0.7qmin)/(1.3qmax-0.7qmin)
其中,p是归一化后的值,q是归一化前的值,qmax,qmin是归一化变量q的最大值、最小值。
4.根据权利要求1所述的配网电采暖负荷预测方法,其特征在于:
给定输出误差εy为0.01,给定神经元接近上限值的误差ε和接近下限值的误差ε均为0.01。
5.根据权利要求1所述的配网电采暖负荷预测方法,其特征在于:
步骤S1的各用户能够按采暖时间分成5类,包括学校、行政事业单位、公司企业、商业宾馆和居民。
6.根据权利要求1所述的配网电采暖负荷预测方法,其特征在于:
在对每一组数据进行训练时的电采暖负荷期望输出y为该组数据中的电采暖负荷值,即采用历年的真实的电采暖负荷值,作为神经网络训练的期望输出值。
7.根据权利要求1所述的配网电采暖负荷预测方法,其特征在于:
学习率η在[0,1]区间内按均匀分布随机产生。
8.根据权利要求1所述的配网电采暖负荷预测方法,其特征在于:
在获取历史数据以及参数设置步骤S1中,itmax=500。
CN201711029996.9A 2017-10-26 2017-10-26 基于改进bp神经网络算法的配网电采暖负荷预测方法 Active CN107748934B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711029996.9A CN107748934B (zh) 2017-10-26 2017-10-26 基于改进bp神经网络算法的配网电采暖负荷预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711029996.9A CN107748934B (zh) 2017-10-26 2017-10-26 基于改进bp神经网络算法的配网电采暖负荷预测方法

Publications (2)

Publication Number Publication Date
CN107748934A CN107748934A (zh) 2018-03-02
CN107748934B true CN107748934B (zh) 2023-06-23

Family

ID=61252728

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711029996.9A Active CN107748934B (zh) 2017-10-26 2017-10-26 基于改进bp神经网络算法的配网电采暖负荷预测方法

Country Status (1)

Country Link
CN (1) CN107748934B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108253519B (zh) * 2017-12-27 2021-06-22 国网北京市电力公司 基于电采暖的用电控制方法
CN108510118A (zh) * 2018-04-02 2018-09-07 张龙 一种基于物联网的楼宇供暖节能预测分析终端
CN109816155A (zh) * 2019-01-04 2019-05-28 国网重庆市电力公司经济技术研究院 一种考虑温度因素的神经网络机器学习负荷预测方法
CN110135619A (zh) * 2019-04-02 2019-08-16 国网能源研究院有限公司 一种中长期电供暖需求的预测方法及系统
CN110083951B (zh) * 2019-04-30 2023-06-02 贵州电网有限责任公司 一种基于变压器相关运行数据的固体绝缘寿命预测方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103124072A (zh) * 2012-12-21 2013-05-29 辽宁省电力有限公司电力科学研究院 考虑负荷特性的电网动态无功优化系统及方法
CN103559655A (zh) * 2013-11-15 2014-02-05 哈尔滨工业大学 基于数据挖掘的微网新型馈线负荷的预测方法
CN103580063A (zh) * 2013-11-13 2014-02-12 国家电网公司 一种基于需方响应消纳大规模并网风电的方法
CN104008427A (zh) * 2014-05-16 2014-08-27 华南理工大学 基于bp神经网络的中央空调冷负荷的预测方法
CN104376389A (zh) * 2014-12-10 2015-02-25 国电南京自动化股份有限公司 基于负载均衡的主从式微电网功率负荷预测系统及其方法
CN106447122A (zh) * 2016-10-12 2017-02-22 国网上海市电力公司 一种区域型能源互联网及其一体化优化规划方法
CN106845701A (zh) * 2017-01-11 2017-06-13 东南大学 一种基于热网和房屋热惯性的综合能源系统优化方法
CN106877409A (zh) * 2017-04-13 2017-06-20 国网山东省电力公司菏泽供电公司 通过电转气储能技术提高电力系统对风电消纳能力的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7502768B2 (en) * 2004-02-27 2009-03-10 Siemens Building Technologies, Inc. System and method for predicting building thermal loads
US8843416B2 (en) * 2009-09-11 2014-09-23 NetESCO LLC Determining energy consumption in a structure

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103124072A (zh) * 2012-12-21 2013-05-29 辽宁省电力有限公司电力科学研究院 考虑负荷特性的电网动态无功优化系统及方法
CN103580063A (zh) * 2013-11-13 2014-02-12 国家电网公司 一种基于需方响应消纳大规模并网风电的方法
CN103559655A (zh) * 2013-11-15 2014-02-05 哈尔滨工业大学 基于数据挖掘的微网新型馈线负荷的预测方法
CN104008427A (zh) * 2014-05-16 2014-08-27 华南理工大学 基于bp神经网络的中央空调冷负荷的预测方法
CN104376389A (zh) * 2014-12-10 2015-02-25 国电南京自动化股份有限公司 基于负载均衡的主从式微电网功率负荷预测系统及其方法
CN106447122A (zh) * 2016-10-12 2017-02-22 国网上海市电力公司 一种区域型能源互联网及其一体化优化规划方法
CN106845701A (zh) * 2017-01-11 2017-06-13 东南大学 一种基于热网和房屋热惯性的综合能源系统优化方法
CN106877409A (zh) * 2017-04-13 2017-06-20 国网山东省电力公司菏泽供电公司 通过电转气储能技术提高电力系统对风电消纳能力的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于云理论和元胞自动机理论的城市配电网空间负荷预测;刘自发 等;《中国电机工程学报》;20130405;第33卷(第10期);98-105 *

Also Published As

Publication number Publication date
CN107748934A (zh) 2018-03-02

Similar Documents

Publication Publication Date Title
CN107748934B (zh) 基于改进bp神经网络算法的配网电采暖负荷预测方法
US11581740B2 (en) Method, system and storage medium for load dispatch optimization for residential microgrid
CN104636822B (zh) 一种基于elman神经网络的居民负荷预测方法
CN103413174B (zh) 基于深度学习方法的短期风速多步预测方法
CN105354620A (zh) 一种风机发电功率的预测方法
CN103106544B (zh) 一种基于t-s型模糊神经网络的光伏发电预测系统
Yu et al. The forecast of the electrical energy generated by photovoltaic systems using neural network method
CN104050517A (zh) 基于grnn神经网络的光伏发电预测方法
CN104636985A (zh) 一种改进bp神经网络的输电线路无线电干扰预测方法
CN109002915A (zh) 基于Kmeans-GRA-Elman模型的光伏电站短期功率预测方法
CN110009135B (zh) 一种基于宽度学习的风电功率预测方法
CN104598765A (zh) 一种基于弹性自适应神经网络的建筑物能耗预测方法
CN105574615A (zh) 一种基于空间相关性与ga的小波-bp神经网络风电功率预测方法
CN109636054A (zh) 基于分类和误差组合预测的太阳能光伏发电量预测方法
CN112149905A (zh) 一种基于小波变换和小波神经网络的光伏电站短期功率预测方法
CN113379142B (zh) 一种基于风速修正和融合模型的短期风电功率预测方法
CN103793887A (zh) 基于自适应增强算法的短期电力负荷在线预测方法
CN116167531A (zh) 一种基于数字孪生的光伏发电预测方法
CN113158573A (zh) 一种基于深度学习的小水电群区域配网最优潮流获取方法
Xu et al. Short-term wind power forecasting based on Elman neural network with particle swarm optimization
CN104134103A (zh) 利用修正的bp神经网络模型预测热油管道能耗的方法
CN104036356A (zh) 一种利用分形算法对电网未来运行状态进行预测的方法
CN110555784A (zh) 基于门控循环神经网络风速修正的风电功率预测方法
CN105205572A (zh) 基于EMD分解和Elman算法的光伏出力预测方法
CN113757781A (zh) 基于bp神经网络预测模型的供热负荷预测节能控制方法及系统

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