CN115222098A - 一种基于分布式在线学习的短期电力负荷预测方法 - Google Patents

一种基于分布式在线学习的短期电力负荷预测方法 Download PDF

Info

Publication number
CN115222098A
CN115222098A CN202210712764.8A CN202210712764A CN115222098A CN 115222098 A CN115222098 A CN 115222098A CN 202210712764 A CN202210712764 A CN 202210712764A CN 115222098 A CN115222098 A CN 115222098A
Authority
CN
China
Prior art keywords
power load
data
meteorological
load prediction
training
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
Application number
CN202210712764.8A
Other languages
English (en)
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.)
Information and Telecommunication Branch of State Grid Qinghai Electric Power Co Ltd
Original Assignee
Information and Telecommunication Branch of State Grid Qinghai Electric Power Co Ltd
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 Information and Telecommunication Branch of State Grid Qinghai Electric Power Co Ltd filed Critical Information and Telecommunication Branch of State Grid Qinghai Electric Power Co Ltd
Priority to CN202210712764.8A priority Critical patent/CN115222098A/zh
Publication of CN115222098A publication Critical patent/CN115222098A/zh
Pending legal-status Critical Current

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
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R22/00Arrangements for measuring time integral of electric power or current, e.g. electricity meters
    • G01R22/06Arrangements for measuring time integral of electric power or current, e.g. electricity meters by electronic methods
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology
    • 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
    • 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Business, Economics & Management (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Primary Health Care (AREA)
  • Artificial Intelligence (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • Operations Research (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Power Engineering (AREA)
  • Quality & Reliability (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种基于分布式在线学习的短期电力负荷预测方法,包括离线和在线两个阶段,离线阶段,在每一个气象传感器端,将气象测量数据和预处理后,结合对应的时间和负荷值,形成该传感器下的训练数据。然后利用极限学习机进行离线回归学习,得到在该气象数据下电力负荷预测模型和对应的模型权重系数。在线阶段,在每一个传感器端,将气象数据预处理后,形成测量数据指纹带入对应的电力负荷预测模型,利用加权求和的方法,得到电力负荷最终估计值。同时,利用在线数据对电力负荷预测模型进行参数更新。该方法具有结构简单,实现方便,预测精度高的优点。

Description

一种基于分布式在线学习的短期电力负荷预测方法
技术领域
本发明属于智能电网和人工智能技术领域,具体涉及一种基于分布式在线学习的短期电力负荷预测方法。
背景技术
电力系统负荷预测是电力系统安全、稳定、经济运行的保障,也是编排调度计划、交易计划的重要基础。随着智能电网的建设,对电网规划、营销与调度提出了更高要求,要求其朝着信息化、互联化和智能化的方向发展。输电量数据是电力市场中一项重要的经济指标,准确地对未来输电量进行预测可对未来的生产规划进行合理安排,进而对电力资源进行有效调度,既能保障人们用电可靠和电网安全稳定运行,也能使电力行业的效益极大地提升。因此,如何精确地对输电量进行预测是目前亟待解决的重要问题。随着电网运行电量数据逐步向繁荣庞杂、多维度的趋势发展,传统的预测方法逐步淘汰。
电力负荷预测大致可以分为传统方法和人工智能方法两种。传统方法又可分为线性回归、自回归移动平均值、灰色模型等,这些方法模型简单、计算速度快,但是由于电力负荷具有非线性的特点,这些传统方法并不能很好地解决非线性问题。近年来人工智能方法快速发展,突破了传统方法不能很好解决非线性问题的不足。人工智能方法主要可分为人工神经网络、支持向量机、深度学习等。随着人工智能理论和技术的快速发展,智能预测方法日益受到重视。通过借助机器学习算法,如BP神经网络、支持向量机、随机森林算法、极限学习机等,可以有效处理高维负荷数据。
随着机器学习算法研究的深入,极限学习机(ELM)由于其结构简单、训练速度快且具有较高的泛化能力等特点,已经在模式识别、信号处理、短期预测等领域有了很多应用成果。相比于传统的基于梯度下降的BP算法来训练单隐层前馈神经网络,ELM具有更好的泛化能力和更快的训练速度。其次,目前大多数负荷预测中都利用气象传感器数据采集气象数据,由于各传感器处理速度和通信时延的不同,从而导致负荷预测实时性不高。
检索发现公开号为CN114066072A和CN113962458A中国专利申请均是利用原始的极限学习机作为学习框架进行离线训练,同时利用人工智能算法对ELM参数进行优化,得到ELM算法参数的初始值,其参数更新都是在离线阶段,并没有利用后续测量数据。
发明内容
本发明的目的在于,针对现有技术存在的缺陷,提出一种基于分布式在线学习短期电力负荷预测方法,该方法利用分布式学习和在线学习技术实现电力物联网下负荷的预测,能够降低实现复杂度,提高网络预测模型的实时性和精度。
为了达到以上目的,本发明提供一种基于分布式在线学习的短期电力负荷预测方法,包括离线和在线两个阶段,具体步骤如下:
离线阶段
步骤1、构建基于分布式的训练数据
步骤1-1、在每一个气象传感器测量端获取气象数据,并结合时间信息,构建时间和气象数据的训练样本指纹,同时记录训练样本指纹所对应的电力负荷实际值,以电力负荷实际值为标签构建训练数据库;
步骤1-2、对训练样本指纹进行预处理,得到每个气象传感器测量端的训练数据集合;
步骤2、分布式离线学习
步骤2-1、在每一个气象传感器测量端,利用极限学习机(ELM)对各自的训练样本数据集合进行回归学习,得到每一个气象传感器测量端的电力负荷预测模型;
步骤2-2、利用训练误差,计算每一个电力负荷预测模型的权重系数;
在线阶段
步骤1、电力负荷预测
步骤1-1、在每一个气象传感器测量端,对得到的时间和气象数据进行数据预处理;
步骤1-2、采用预处理后的时间和气象数据,利用离线阶段步骤2-1得到的电力负荷预测模型,分别计算电力负荷预测值;
步骤1-3、结合离线阶段步骤2-2计算得到的模型的权重系数,利用线性权重求和的方法,得到电力负荷预测的最终值;
步骤2、电力负荷预测模型更新
利用在线阶段收集的训练数据,对电力负荷预测模型的参数进行更新。
本发明将分布式学习和在线学习用于电力负荷短期预测,分布式学习能够克服传感器所测量的气象数据不同步的问题,而在线学习不需要将全部的训练样本重新代入模型训练,只需要利用新的训练样本对负荷预测模型进行更新。本发明采用在线顺序极限学习机(OS-ELM)进行学习,OS-ELM与ELM在学习形式上完全不一样。同时本发明利用后续的测量值数据作为训练样本,对已经得到的ELM模型进行参数更新,减少学习时间。本发明的参数更新从本质上讲是把后续的样本作为训练集训练,并且参数更新位于在线阶段。
本发明进一步的采用如下技术方案:
所述离线阶段的步骤1-2中训练样本指纹预处理具体步骤如下:
步骤1-2-1、利用中值滤波技术对传感器测量的气象数据进行滤波,以去除传感器的测量噪声;
步骤1-2-2、对滤波后的气象数据进行归一化处理;
步骤1-2-3、以处理后的气象数据和时间作为训练样本指纹,构成训练数据集合。
所述步骤1-2-2中,采用Min-Max归一化对气象数据进行归一化处理,具体计算公式如下:
Figure BDA0003707485920000041
式中,xn为待处理数据在归一化前的值,x′n为待处理数据在归一化后的值,xmax为待处理数据的最大值,xmin为待处理数据的最小值。
所述离线阶段的步骤2-1中,利用极限学习机训练样本数据集合进行回归学习到电力负荷预测模型的具体操作如下:
步骤2-1-1、假设有N个独立样本(xj,tj),j=1,2,…N,那么单隐层前馈神经网络可以描述为:
Figure BDA0003707485920000042
式中,L表示隐层节点的数目,βi表示第i个隐层节点与输出层的连接权值,h(x)表示激活函数,ai,表示第i个隐层节点与输入层的连接权值且i=1,2,…L,xj表示第j个独立样本,bi表示第i个隐层节点的偏差值且i=1,2,…L,oj表示输出层的输出;
步骤2-1-2、采用下式使得ELM算法的误差尽可能小:
Figure BDA0003707485920000043
式中,tj表示第j个训练样本标签;
步骤2-1-3、由于(ai,bi)是随机生成的,在学习过程中不需要调整,即存在βi使得:
Figure BDA0003707485920000051
公式矩阵形式可以表示为Hβ=T,其中
Figure BDA0003707485920000052
Figure BDA0003707485920000053
式中,H是隐层输出矩阵,β为输出权重且β=[β1,....βL]T
步骤2-1-4、由于N>>L,即训练样本数远大于隐层节点数,为了使误差最小,即:
Figure BDA0003707485920000054
式中,
Figure BDA0003707485920000055
为输出权重β的估计值;
那么线性方程Hβ=T的最小二乘解为:
Figure BDA0003707485920000056
最后,通过求解β,得到每个传感器测量端的电力负荷预测模型:
Figure BDA0003707485920000057
其中
Figure BDA0003707485920000058
表示电力负荷预测值,
Figure BDA0003707485920000059
表示气象测量值,
Figure BDA00037074859200000510
为βi的估计值。
所述离线阶段步骤2-2中计算电力负荷预测模型权重系数的具体步骤如下:
定义分布式学习中第i个气象传感器的离线训练误差为ei,则该气象传感器所测量的气象数据对应的电力负荷预测模型的权重系数为wi,其计算公式为:
Figure BDA0003707485920000061
所述在线阶段步骤1-3中,采用下式计算电力负荷预测的最终值λ:
Figure BDA0003707485920000062
式中,λi为第i个气象传感器测量端的电力负荷预测值。
所述在线阶段步骤2中电力负荷预测模型更新的具体操作如下:
步骤2-1、根据后续电力负荷和气象信息的测量值,构造在线数据更新集
在每一个气象传感器测量端,将对应的气象测量数据和测量时间经过数据预处理后,得到在线更新样本的指纹信息;
步骤2-2、在线负荷预测模型参数更新
利用在线极限学习机(OS-ELM)中的参数更新处理技术,对电力负荷预测模型中的参数进行更新。
所述在线阶段步骤2-2中,在线负荷预测模型参数更新的具体操作如下:
步骤2-2-1、对于离线阶段的训练数据块Γ0={chunk1},数据块大小为N0,离线训练考虑最小化目标函数为||H0β-T0||,
初始隐层输出矩阵H0为:
Figure BDA0003707485920000063
其中,x1,
Figure BDA0003707485920000064
分别表示训练数据块Γ0中第1个和第N0个训练数据的气象测量值,T0表示训练数据块Γ0中训练数据的负荷值;
则初始输出权重β(0)可以表示为:
Figure BDA0003707485920000065
Figure BDA0003707485920000066
步骤2-2-2、在线阶段增加一个训练数据块Γ1={chunk2},数据块大小为N1,那么最小化目标函数为
Figure BDA0003707485920000071
隐层输出矩阵H1
Figure BDA0003707485920000072
输出权重β(1)可以表示为:
Figure BDA0003707485920000073
其中,
Figure BDA0003707485920000074
推导可得:
Figure BDA0003707485920000075
其中,
Figure BDA0003707485920000076
分别表示训练数据块Γ1中第1个和第N1个训练数据的气象测量值,T1表示训练数据块Γ1中训练数据的负荷值,且
Figure BDA0003707485920000077
步骤2-2-3、利用递推法,当增加数据块Γk+1={chunkk+2}时,数据块的大小为Nk+1,对应的隐层输出矩阵Hk+1可以表示为:
Figure BDA0003707485920000078
那么输出权重β(k+1)为:
Figure BDA0003707485920000079
其中,
Figure BDA0003707485920000081
分别表示训练数据块Γk+1中第1个和第Nk+1个训练数据的气象测量值,且
Figure BDA0003707485920000082
Tk+1表示训练数据块Γk+1中训练数据的负荷值。
所述步骤2-2-3中,在计算β(k+1)的过程中用到了
Figure BDA0003707485920000083
为了避免在递归过程中计算矩阵的逆,使用Woodbury公式推导出方程式:
Figure BDA0003707485920000084
Figure BDA0003707485920000085
那么最终更新过程可以表示为:
Figure BDA0003707485920000086
Figure BDA0003707485920000087
式中,I表示单位阵。
本发明的方法中离线阶段,在每一个气象传感器端,将气象测量数据和预处理后,结合对应的时间和负荷值,形成该传感器下的训练数据。然后利用极限学习机进行离线回归学习,得到在该气象数据下电力负荷预测模型和对应的模型权重系数。在线阶段,在每一个传感器端,将气象数据预处理后,形成测量数据指纹带入对应的电力负荷预测模型,利用加权求和的方法,得到电力负荷最终估计值。同时,利用在线数据对电力负荷预测模型进行参数更新。该方法具有结构简单,实现方便,预测精度高的优点。
本发明的有益效果如下:
(1)本发明利用极限学习机来学习传感器数据和电力负荷之间的非线性关系。由于极限学习机具有更好的泛化能力和更快的训练速度,因此能够得到更好的预测性能。
(2)本发明利用分布式学习框架训练负荷预测模型,每一个传感器测量端分别利用各自数据测量值进行离线学习。因此,能克服传感器测量数据不同步的问题,提高预测实时性。
(3)本发明利用在线学习方法来更新预测模型参数。在线学习中不需要将全部的训练样本重新代入模型训练,而只需要利用新的训练样本对负荷预测模型进行更新。因此能够利用后续测量数据得到更精确的预测模型,从而提高预测精度。
附图说明
下面结合附图对本发明作进一步的说明。
图1为本发明的流程图。
图2为本发明中ELM学习算法的结构框图。
图3为本发明中在线电力负荷预测的流程图。
图4为本发明中试验结果的RMSE图。
具体实施方式
实施例一
一种基于分布式在线学习的短期电力负荷预测方法,是基于分布式学习和在线学习理论对电力负荷进行预测。如图1所示,本发明的方法包括离线和在线两个阶段,具体步骤如下:
离线阶段
步骤1、构建基于分布式的训练数据
步骤1-1、在每一个气象传感器测量端获取气象数据,并结合时间信息,构建时间和气象数据的训练样本指纹,同时记录训练样本指纹所对应的电力负荷实际值,以电力负荷实际值为标签构建训练数据库。气象数据包括温度、压强、风速、湿度。
步骤1-2、对时间和气象数据的训练样本指纹进行预处理,得到每个气象传感器测量端的训练数据集合。本发明采用的预处理方法为中值滤波和归一化处理。训练样本指纹预处理的具体步骤如下:
步骤1-2-1、利用中值滤波技术对传感器测量的气象数据进行滤波,以去除传感器的测量噪声。中值滤波是一种非线性平滑技术,它将每一点的值设置为该点某邻域窗口内的所有点值的中值,从而消除孤立的噪声点。具体的操作过程如下:用一个奇数点的移动窗口,将窗口中心点的值用窗口内各点的中值代替。假设窗口内有5个点,其值为1,2,3,4和5,那么此窗口内各点的中值即为3,也就是用3来代替中心点的像素值。中值滤波在滤除噪声的同时,能够保护信号的边缘,并且算法简单,易于实现。
步骤1-2-2、对滤波后的气象数据进行归一化处理。在电力负荷预测的传感器测量数值中,不同测量参数具有不同的量纲和量纲单位,数值间的差别很大,从而影响离线学习结果。为消除测量值量纲和取值范围差异影响,本发明采用Min-Max归一化,对数据进行规范化即归一化处理。原始数据经过数据标准化处理后,各指标处于同一数量级,落入一个特定的区域,便于进行综合对比评价。Min-Max归一化又称为最大最小标准化,它是通过将不同量级的数据统一转化为同一个量级,利用线性函数将原始数据线性化的方法转换到[0,1]的范围,计算结果为归一化后的数据,能够适用在数值比较集中的场景下。
采用Min-Max归一化对气象数据进行归一化处理,具体计算公式如下:
Figure BDA0003707485920000101
式中,xn为待处理数据在归一化前的值,x′n为待处理数据在归一化后的值,xmax为待处理数据的最大值,xmin为待处理数据的最小值。最大最小标准化最大的优点是简单,容易计算,能够应用于数值型的数据,并且不受数据量级的影响。
经过数据预处理后,每一个传感器的训练样本指纹信息为(时间,测量值),大小为1x2。
步骤1-2-3、以处理后的气象数据和时间作为训练样本指纹,构成训练数据集合。
步骤2、分布式离线学习
步骤2-1、在每一个气象传感器测量端,利用极限学习机(ELM)对各自的训练样本数据集合进行离线回归学习,得到每一个气象传感器测量端的电力负荷预测模型。离线回归学习利用分布式机器学习框架,分别在每个传感器测量端进行离线学习。
在介绍ELM算法之前,需要先了解Moore-Penrose广义逆以及线性方程的最小范数最小二乘解的定义。Moore-Penrose广义逆的定义为:如果B(m×n的矩阵)是A(n×m的矩阵)的Moore-Penrose广义逆,那么:
ABA=A,BAB=B,
(AB)T=AB,(BA)T=BA
并且B可以记为A+
线性方程的最小范数最小二乘解的定义为:如果x1是线性方程Hx=y的最小范数最小二乘解,x1∈Rn,y∈Rm,那么:
Figure BDA0003707485920000111
其中,||·||代表的是欧式空间的范数。
定理1:对于矩阵G,Gy是线性方程Hx=y的最小范数最小二乘解的充要条件为:G=H+,即矩阵G是H的Moore-Penrose广义逆。
如图2所示,ELM算法的结构框中最左边是输入层,n表示样本的特征维度,xj表示第j个独立样本,可以表示为xj=[xj1,xj2,…xjn]T;中间的是隐层,L表示隐层节点的数目,L的值越大,网络的表达能力越强,第i个隐层节点与输入层的连接权值为ai,i=1,2,…L,可以表示为ai=[ai1,ai2,…,ain]T,是一个n维向量,第i个隐层节点的偏差值为bi,是一个数值,(ai,bi)是随机生成的,学习过程中不需要调整,h(x)表示的是激活函数,第i个隐层节点与输出层的连接权值为βi;最右边的是输出层,输出为oj
本发明中利用极限学习机训练样本数据集合进行回归学习到电力负荷预测模型的具体操作如下:
假设有N个独立样本(xj,tj),j=1,2,…N,那么这个单隐层前馈神经网络(SLFN)可以描述为:
Figure BDA0003707485920000121
式中,L表示隐层节点的数目,βi表示第i个隐层节点与输出层的连接权值,h(x)表示激活函数,ai,表示第i个隐层节点与输入层的连接权值且i=1,2,…L,xj表示第j个独立样本,bi表示第i个隐层节点的偏差值且i=1,2,…L,oj表示输出层的输出;
ELM算法的目的是让误差尽可能小,即:
Figure BDA0003707485920000122
式中,tj表示第j个训练样本标签;
由于(ai,bi)是随机生成的,在学习过程中不需要调整,即存在βi使得:
Figure BDA0003707485920000123
公式矩阵形式可以表示为Hβ=T,其中
Figure BDA0003707485920000124
Figure BDA0003707485920000131
式中,H是隐层输出矩阵,β为输出权重且β=[β1,....βL]T,β1为第1个输出权重,βL为第L个输出权重;
对于公式Hβ=T,当N=L时,即训练样本数等于隐层节点数,H为可逆方阵,那么β=H-1T,单隐层前馈神经网络(SLFN)的误差可以接近于零;但是在实际应用中,由于N>>L,即训练样本数远大于隐层节点数,为了使误差最小,即:
Figure BDA0003707485920000132
式中,
Figure BDA0003707485920000133
为输出权重β的估计值;
根据定理1,线性方程Hβ=T的最小二乘解为:
Figure BDA0003707485920000134
因此,通过求解β,就能得到每个传感器测量端的电力负荷预测模型,输出权重参数。传感器测量端的电力负荷预测模型如下所示:
Figure BDA0003707485920000135
其中,
Figure BDA0003707485920000136
表示电力负荷预测值,
Figure BDA0003707485920000137
表示气象测量值,
Figure BDA0003707485920000138
为βi的估计值。
步骤2-2、利用训练误差,计算每一个电力负荷预测模型的权重系数。计算电力负荷预测模型权重系数的具体步骤如下:
定义分布式学习中第i个气象传感器的离线训练误差为ei,则该气象传感器所测量的气象数据对应的电力负荷预测模型的权重系数为wi,其计算公式为:
Figure BDA0003707485920000139
离线阶段利用分布式的学习框架,在每一个传感器测量端,对气象测量数据和时间进行预处理,形成训练数据库。然后利用极限学习机进行离线回归学习,得到每个测量值的电力负荷预测模型和其所对应的权重系数。
在线阶段
步骤1、电力负荷预测,其流程如图3所示
步骤1-1、在第i个气象传感器测量端,对得到的时间和气象数据进行与离线阶段的步骤1-2相同的数据预处理。
步骤1-2、将预处理后的数据带入离线阶段步骤2-1得到的电力负荷预测模型中,分别计算电力负荷预测值,得到电力负荷预测值λi
步骤1-3、结合离线阶段步骤2-2计算得到的模型的权重系数,利用线性权重求和的方法,得到电力负荷预测的最终值。采用下式计算电力负荷预测的最终值λ:
Figure BDA0003707485920000141
式中,λi为第i个气象传感器测量端的电力负荷预测值。
步骤2、电力负荷预测模型更新
利用在线阶段收集的训练数据,采用在线极限学习机(OS-ELM)技术中参数更新算法对每个电力负荷预测模型中的参数进行更新。电力负荷预测模型更新的具体操作如下:
步骤2-1、根据后续电力负荷和气象信息的测量值,构造在线数据更新集
在第i个气象传感器测量端,将对应的气象测量数据和测量时间经过数据预处理后,得到在线更新样本的指纹信息。
步骤2-2、在线负荷预测模型参数更新
利用在线极限学习机(OS-ELM)中的参数更新处理技术,对电力负荷预测模型中的参数进行更新。在线负荷预测模型参数更新的具体操作如下:
步骤2-2-1、对于离线阶段的训练数据块Γ0={chunk1},该数据块大小为N0,离线训练考虑最小化目标函数为||H0β-T0||,初始隐层输出矩阵H0为:
Figure BDA0003707485920000151
a1表示第1个隐层节点与输入层的连接权值,aL表示第L个隐层节点与输入层的连接权值,b1表示第1个隐层节点的偏差值,bL表示第L个隐层节点的偏差值,x1,
Figure BDA0003707485920000152
分别表示训练数据块Γ0中第1个和第N0个训练数据的气象测量值,T0表示训练数据块Γ0中训练数据的负荷值;初始输出权重β(0)可以表示为:
Figure BDA0003707485920000153
其中K0为系数且
Figure BDA0003707485920000154
步骤2-2-2、在线阶段增加一个训练数据块Γ1={chunk2},该数据块大小为N1,那么最小化目标函数为
Figure BDA0003707485920000155
隐层输出矩阵H1为:
Figure BDA0003707485920000156
输出权重β(1)可以表示为:
Figure BDA0003707485920000157
其中,
Figure BDA0003707485920000158
推导可得:
Figure BDA0003707485920000159
其中,
Figure BDA00037074859200001510
分别表示训练数据块Γ1中第1个和第N1个训练数据的气象测量值,T1表示训练数据块Γ1中训练数据的负荷值,K1为系数且
Figure BDA0003707485920000161
步骤2-2-3、利用递推法,当增加数据块Γk+1={chunkk+2}时,数据块的大小为Nk+1,对应的隐层输出矩阵Hk+1可以表示为:
Figure BDA0003707485920000162
那么输出权重β(k+1)为:
Figure BDA0003707485920000163
其中,
Figure BDA0003707485920000164
分别表示训练数据块Γk+1中第1个和第Nk+1个训练数据的气象测量值,Kk+1为系数且
Figure BDA0003707485920000165
Tk+1表示训练数据块Γk+1中训练数据的负荷值,
Figure BDA0003707485920000166
Figure BDA0003707485920000167
分别表示训练数据块Γk+1中第1个和第Nk+1个电力负荷值。
在计算β(k+1)的过程中用到了
Figure BDA0003707485920000168
为了避免在递归过程中计算矩阵的逆,使用Woodbury公式推导出方程式:
Figure BDA0003707485920000169
Figure BDA00037074859200001610
那么最终更新过程可以表示为:
Figure BDA00037074859200001611
Figure BDA00037074859200001612
式中,
Figure BDA00037074859200001613
I表示单位阵。
在线阶段,在各个传感器测量端,利用电力负荷预测模型得到电力负荷预测中间值,然后利用各个模型权重系数进行线性加权,得到电力负荷最终预测值。此外,为了充分利用后续训练数据,利用在线学习的方式,对电力负荷预测模型进行参数更新,获得更好的预测模型。
为了验证算法性能,本发明对算法进行了电力负荷预测实验。当训练样本数目为9000,图4描述了不同测试样本下的估计误差的均方根误差(RMSE)。由图4可知,随着测试数目的增加,仿真误差RMSE呈递减趋势。因此通过对预测误差的统计分析,本发明的方法完全能够满足实际的电力负荷预测场景。
除上述实施例外,本发明还可以有其他实施方式。凡采用等同替换或等效变换形成的技术方案,均落在本发明要求的保护范围。

Claims (9)

1.一种基于分布式在线学习的短期电力负荷预测方法,其特征在于,包括离线和在线两个阶段,具体步骤如下:
离线阶段
步骤1、构建基于分布式的训练数据
步骤1-1、在每一个气象传感器测量端获取气象数据,并结合时间信息,构建时间和气象数据的训练样本指纹;
步骤1-2、对训练样本指纹进行预处理,得到每个气象传感器测量端的训练数据集合;
步骤2、分布式离线学习
步骤2-1、在每一个气象传感器测量端,利用极限学习机对各自的训练样本数据集合进行回归学习,得到每一个气象传感器测量端的电力负荷预测模型;
步骤2-2、利用训练误差,计算每一个电力负荷预测模型的权重系数;
在线阶段
步骤1、电力负荷预测
步骤1-1、在每一个气象传感器测量端,对得到的时间和气象数据进行数据预处理;
步骤1-2、采用预处理后的时间和气象数据,利用离线阶段步骤2-1得到的电力负荷预测模型,分别计算电力负荷预测值;
步骤1-3、结合离线阶段步骤2-2计算得到的模型的权重系数,利用线性权重求和的方法,得到电力负荷预测的最终值;
步骤2、电力负荷预测模型更新
利用在线阶段收集的训练数据,对电力负荷预测模型的参数进行更新。
2.根据权利要求1所述一种基于分布式在线学习的短期电力负荷预测方法,其特征在于,所述离线阶段的步骤1-2中训练样本指纹预处理具体步骤如下:
步骤1-2-1、利用中值滤波技术对传感器测量的气象数据进行滤波,以去除传感器的测量噪声;
步骤1-2-2、对滤波后的气象数据进行归一化处理;
步骤1-2-3、以处理后的气象数据和时间作为训练样本指纹,构成训练数据集合。
3.根据权利要求2所述一种基于分布式在线学习的短期电力负荷预测方法,其特征在于,所述步骤1-2-2中,采用Min-Max归一化对气象数据进行归一化处理,具体计算公式如下:
Figure FDA0003707485910000021
式中,xn为待处理数据在归一化前的值,x’n为待处理数据在归一化后的值,xmax为待处理数据的最大值,xmin为待处理数据的最小值。
4.根据权利要求1所述一种基于分布式在线学习的短期电力负荷预测方法,其特征在于,所述离线阶段的步骤2-1中,利用极限学习机训练样本数据集合进行回归学习到电力负荷预测模型的具体操作如下:
步骤2-1-1、假设有N个独立样本(xj,tj),j=1,2,…N,那么单隐层前馈神经网络可以描述为:
Figure FDA0003707485910000022
式中,L表示隐层节点的数目,βi表示第i个隐层节点与输出层的连接权值,h(x)表示激活函数,ai,表示第i个隐层节点与输入层的连接权值且i=1,2,…L,xj表示第j个独立样本,bi表示第i个隐层节点的偏差值且i=1,2,…L,oj表示输出层的输出;
步骤2-1-2、采用下式使得ELM算法的误差尽可能小:
Figure FDA0003707485910000031
式中,tj表示第j个训练样本标签;
步骤2-1-3、由于(ai,bi)是随机生成的,在学习过程中不需要调整,即存在βi使得:
Figure FDA0003707485910000032
公式矩阵形式可以表示为Hβ=T,其中
Figure FDA0003707485910000033
Figure FDA0003707485910000034
式中,H是隐层输出矩阵,β为输出权重且β=[β1,....βL]T
步骤2-1-4、由于N>>L,即训练样本数远大于隐层节点数,为了使误差最小,即:
Figure FDA0003707485910000035
式中,
Figure FDA0003707485910000036
为输出权重β的估计值;
那么线性方程Hβ=T的最小二乘解为:
Figure FDA0003707485910000037
最后,通过求解β,得到每个传感器测量端的电力负荷预测模型:
Figure FDA0003707485910000038
其中
Figure FDA0003707485910000041
表示电力负荷预测值,
Figure FDA0003707485910000042
表示气象测量值,
Figure FDA0003707485910000043
为βi的估计值。
5.根据权利要求1所述一种基于分布式在线学习的短期电力负荷预测方法,其特征在于,所述离线阶段步骤2-2中计算电力负荷预测模型权重系数的具体步骤如下:
定义分布式学习中第i个气象传感器的离线训练误差为ei,则该气象传感器所测量的气象数据对应的电力负荷预测模型的权重系数为wi,其计算公式为:
Figure FDA0003707485910000044
6.根据权利要求1所述一种基于分布式在线学习的短期电力负荷预测方法,其特征在于,所述在线阶段步骤1-3中,采用下式计算电力负荷预测的最终值λ:
Figure FDA0003707485910000045
式中,λi为第i个气象传感器测量端的电力负荷预测值。
7.根据权利要求1所述一种基于分布式在线学习的短期电力负荷预测方法,其特征在于,所述在线阶段步骤2中电力负荷预测模型更新的具体操作如下:
步骤2-1、根据后续电力负荷和气象信息的测量值,构造在线数据更新集
在每一个气象传感器测量端,将对应的气象测量数据和测量时间经过数据预处理后,得到在线更新样本的指纹信息;
步骤2-2、在线负荷预测模型参数更新
利用在线极限学习机中的参数更新处理技术,对电力负荷预测模型中的参数进行更新。
8.根据权利要求6所述一种基于分布式在线学习的短期电力负荷预测方法,其特征在于,所述在线阶段步骤2-2中,在线负荷预测模型参数更新的具体操作如下:
步骤2-2-1、对于离线阶段的训练数据块Γ0={chunk1},数据块大小为N0,离线训练考虑最小化目标函数为||H0β-T0||,
初始隐层输出矩阵H0为:
Figure FDA0003707485910000051
其中,x1,
Figure FDA0003707485910000052
分别表示训练数据块Γ0中第1个和第N0个训练数据的气象测量值,T0表示训练数据块Γ0中训练数据的负荷值;
则初始输出权重β(0)可以表示为:
Figure FDA0003707485910000053
Figure FDA0003707485910000054
步骤2-2-2、在线阶段增加一个训练数据块Γ1={chunk2},数据块大小为N1,那么最小化目标函数为
Figure FDA0003707485910000055
隐层输出矩阵H1为:
Figure FDA0003707485910000056
输出权重β(1)可以表示为:
Figure FDA0003707485910000057
其中,
Figure FDA0003707485910000058
推导可得:
Figure FDA0003707485910000059
其中,
Figure FDA00037074859100000510
分别表示训练数据块Γ1中第1个和第N1个训练数据的气象测量值,T1表示训练数据块Γ1中训练数据的负荷值,且
Figure FDA00037074859100000511
步骤2-2-3、利用递推法,当增加数据块Γk+1={chunkk+2}时,数据块的大小为Nk+1,对应的隐层输出矩阵Hk+1可以表示为:
Figure FDA0003707485910000061
那么输出权重β(k+1)为:
Figure FDA0003707485910000062
其中,
Figure FDA0003707485910000063
分别表示训练数据块Γk+1中第1个和第Nk+1个训练数据的气象测量值,且
Figure FDA0003707485910000064
Tk+1表示训练数据块Γk+1中训练数据的负荷值。
9.根据权利要求8所述一种基于分布式在线学习的短期电力负荷预测方法,其特征在于,所述步骤2-2-3中,在计算β(k+1)的过程中用到了
Figure FDA0003707485910000065
为了避免在递归过程中计算矩阵的逆,使用Woodbury公式推导出方程式:
Figure FDA0003707485910000066
Figure FDA0003707485910000067
那么最终更新过程可以表示为:
Figure FDA0003707485910000068
Figure FDA0003707485910000069
式中,I表示单位阵。
CN202210712764.8A 2022-06-22 2022-06-22 一种基于分布式在线学习的短期电力负荷预测方法 Pending CN115222098A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210712764.8A CN115222098A (zh) 2022-06-22 2022-06-22 一种基于分布式在线学习的短期电力负荷预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210712764.8A CN115222098A (zh) 2022-06-22 2022-06-22 一种基于分布式在线学习的短期电力负荷预测方法

Publications (1)

Publication Number Publication Date
CN115222098A true CN115222098A (zh) 2022-10-21

Family

ID=83607639

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210712764.8A Pending CN115222098A (zh) 2022-06-22 2022-06-22 一种基于分布式在线学习的短期电力负荷预测方法

Country Status (1)

Country Link
CN (1) CN115222098A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115937355A (zh) * 2023-02-14 2023-04-07 国网湖北省电力有限公司电力科学研究院 电力物联网负荷预测指纹图像构造方法、系统及介质

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115937355A (zh) * 2023-02-14 2023-04-07 国网湖北省电力有限公司电力科学研究院 电力物联网负荷预测指纹图像构造方法、系统及介质

Similar Documents

Publication Publication Date Title
Ke et al. Short-term electrical load forecasting method based on stacked auto-encoding and GRU neural network
Xiaoyun et al. Short-term prediction of wind power based on deep long short-term memory
CN110414788B (zh) 一种基于相似日和改进lstm的电能质量预测方法
CN111461463B (zh) 一种基于tcn-bp的短期负荷预测方法、系统及设备
CN113554466B (zh) 一种短期用电量预测模型构建方法、预测方法和装置
CN111178616B (zh) 基于负相关学习和正则化极限学习机集成的风速预测方法
CN106503867A (zh) 一种遗传算法最小二乘风电功率预测方法
CN111241755A (zh) 电力负荷预测方法
CN113095598A (zh) 一种多能负荷预测方法、系统、设备和介质
CN114462718A (zh) 基于时间滑动窗口的cnn-gru风电功率预测方法
CN111898825A (zh) 一种光伏发电功率短期预测方法及装置
CN115456287A (zh) 一种基于长短时记忆网络的综合能源系统多元负荷预测方法
CN110738363B (zh) 一种光伏发电功率预测方法
CN115358437A (zh) 基于卷积神经网络的供电负荷预测方法
CN115169742A (zh) 一种短期风力发电功率预测方法
CN113128666A (zh) 基于Mo-S-LSTMs模型的时间序列多步预测方法
CN117151770A (zh) 一种基于注意力机制的lstm碳价预测方法和系统
CN116975645A (zh) 一种基于vae-mrcnn的工业过程软测量建模方法
CN114091647A (zh) 一种基于bp神经网络的太阳10.7厘米射电流量预报方法
CN115222098A (zh) 一种基于分布式在线学习的短期电力负荷预测方法
Zhao et al. A dendritic neuron based very short-term prediction model for photovoltaic power
CN117407704A (zh) 可再生能源发电功率预测方法及其计算机设备、储存介质
CN114254828B (zh) 一种基于混合卷积特征提取器和gru的电力负荷预测方法
CN115759343A (zh) 一种基于e-lstm的用户电量预测方法和装置
CN116129212A (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