CN111552923B - 一种基于通用分布的负荷预测方法及负荷预测*** - Google Patents
一种基于通用分布的负荷预测方法及负荷预测*** Download PDFInfo
- Publication number
- CN111552923B CN111552923B CN202010322604.3A CN202010322604A CN111552923B CN 111552923 B CN111552923 B CN 111552923B CN 202010322604 A CN202010322604 A CN 202010322604A CN 111552923 B CN111552923 B CN 111552923B
- Authority
- CN
- China
- Prior art keywords
- temperature
- power
- load
- data
- histogram
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000005315 distribution function Methods 0.000 claims abstract description 36
- 238000012545 processing Methods 0.000 claims description 28
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000003491 array Methods 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 15
- 230000001186 cumulative effect Effects 0.000 claims description 10
- 238000000926 separation method Methods 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000009467 reduction Effects 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 12
- 230000000694 effects Effects 0.000 description 9
- 238000013528 artificial neural network Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 238000012549 training Methods 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000000611 regression analysis Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010205 computational analysis Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000013277 forecasting method Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 238000005057 refrigeration Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Operations Research (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Mathematical Physics (AREA)
- Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Evolutionary Biology (AREA)
- Development Economics (AREA)
- Primary Health Care (AREA)
- Probability & Statistics with Applications (AREA)
- Quality & Reliability (AREA)
- Algebra (AREA)
- Entrepreneurship & Innovation (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种基于通用分布的负荷预测方法及预测***,以实现电力***短期负荷功率预测。该方法将历史气象信息数据和电气数据相结合,针对不同的温度增量水平分别构建相应的统计分析模型。首先采用直方图模型来近似表征实际分布,然后选取适当的模型对实际分布进行拟合,本文采用通用分布函数模型进行拟合,并且依据通用分布CDF函数的逆函数的闭合解析表达式,得出一定置信水平下预测功率的置信区间。
Description
技术领域:
本发明属于电力负荷预测技术领域,具体涉及一种基于通用分布的气象因素对负荷用电影响的负荷预测方法及预测***。
背景技术:
随着国民经济水平的提高,对电力的需求也越来越大,要进行电力的科学合理的调度,需要考虑的因素很多,其中气温是重要的影响因素。气温的变化在民用电和工业用电上都会造成区域负荷功率的显著变化。一方面,民用大功率耗电设备(例如制冷、取暖等用电设备)的使用与温度等气象因素密切相关;另一方面,气温变化同样会影响企业生产计划从而造成工业负荷功率的变化。因此,分析负荷功率与气温的定量关系并将其应用到负荷预测过程中是极为必要的。
目前,现有的用电量预测的研究方法可以分为两大类,一类是基于概率统计的间接计算方法,较为常见的是多元回归模型。在回归分析中,回归方程的因变量往往是电力负荷,可以是短期负荷,也可以是长期负荷,自变量是影响电力负荷的各种因素,比如经济、政策、电价和气候等。而回归分析的缺点在于它的实现要求有一定的样本数量,并且样本要有较好的分布规律和确定的发展趋势,计算工作量较大。另一类是基于用电量历史数据的间接法。比如常规数学模型,例如灰色预测模型和神经网络模型。对于电力***,对其影响的供电机组、电网容量、用户情况、等信息是已知的,但是,影响负荷的其他很多因素,比如气象情况、行政与管理政策的变化、地区经济活动等是难以确切知道的,因此,电力负荷是灰色***。灰色理论将无规律的历史数据列经累加生成后,使其变为具有指数增长规律的上升形状数列,由于一阶微分方程解的形成即是指数增长形式,所以可对生成后数列建立微分方程模型。所以灰色模型实际上是生成数列所建模型。GM模型所得数据必须经过逆生成,即累减生成做还原后才能应用。其缺点在于预测过程没有考虑负荷受经济、气候、政策的影响,而是直接对用电量的时间序列数据进行处理分析,此外,原始数据离散程度越大,灰度越大,预测精度越低。神经网络理论利用神经网络的学习功能,让计算机学习包含在历史负荷数据中的映射关系,再利用这种映射关系预测未来负荷。BP神经网络能够模拟任意输入和输出之间的非线性映射关系。一个三层BP网络由输入层、隐含层和输出层组成。训练方法是,对于一组输入样本,通过BP神经网络计算实际输出。用BP网络的实际输出与输出样本之间的误差来修正网络的连接权值,直至二者的误差达到设定值。一般输入样本包含负荷、温度、天气情况和日期类型,训练时的目标样本是分时段负荷,模型应用是由前一条的负荷数据来预测第二天的负荷值。其缺点在于处理非文本类历史数据信息和调度人员的已有经验知识时存在一定的困难;对输入数据要求高,在样本数比较小的情况下,预测精度降低;对于离训练样本比较远的数据,预测效果也不理想。
现有技术中,基于数值天气预报集合预报结果的风电功率概率预测方法,具体方法步骤为:1.对数值天气预报集合预报的每个成员建立短期风电功率预测模型,分别输入数值天气预测结果后获得各个所述成员的风功率预测结果,得到天气预报集合预报的集合的风功率预测结果。2.根据所述集合中各成员的所述风电功率预测结果识别各个所述集合的误差类型。3.根据所述集合的风功率预测结果将各个所述集合的功率水平进行划分。4.计算不同误差类型和功率水平的集合的相对误差的集合。5.用核估计的方法得到所述集合中各个样本的概率密度函数。6.采用非参数拟合的方法拟合概率密度分布和拟合回归函数。7.对拟合回归结果进行回归校验。8.计算一定置信水平下的误差上限和误差下限。9.根据误差上下限计算得到一定置信度水平下功率预测的估计区间。
在现有的用电量预测的研究方法中,多元回归模型的模型简单、理论完善,能够充分考虑各种影响因素,但容易因为因素选取不当而造成预测精度不高。常规数学模型短期预测效果好,但是过度依赖于数学物理机制,在长期预测和现实应用中有较大的局限性。
发明内容:
针对现有技术存在的以上问题,本发明公开了一种基于通用分布的负荷预测方法及预测***,并应用于负荷功率预测。
本发明的上述技术问题主要通过以下技术方案得以解决:
一种基于通用分布负荷预测方法,包括以下步骤:
步骤1,采集设定时间段内待预测负荷的分时刻负荷功率及温度的历史数据:
步骤2,对步骤1所采集负荷功率以及温度历史数据的进行处理,计算得到包含温度增量和功率比值的数组,其中,所述温度增量是指每个时间点同一时刻后一天减前一天的温度增量,所述功率比值是指每个时间点同一时刻后一天的功率与前一天的功率的比值;
步骤3,对步骤2获得的包含温度增量和功率比值的数组进行分箱处理,所有数组分别存入其对应的温度增量等级数据箱中:
步骤4,基于直方图,建立各等级数据箱中负荷功率实际分布模型:
步骤5,采用通用分布函数模型对各等级数据箱内数据进行拟合,得到每一等级数据箱所对应的通用分布函数;
步骤6,根据待预测负荷的当前功率值和预测日的温度,基于步骤5拟合得到的通用分布函数,确定待预测负荷在预测日温度条件下的功率比值置信区间,实现给定气象下的功率预测。
本发明进一步包括以下优选方案。
在步骤1中,采集待预测负荷在最近三个月内的全天24小时分时刻的负荷功率以及对应时刻的温度数据;其中,在历史数据中,负荷功率的采样周期为1min,即对应每个负荷功率数据点为1min时刻的有功功率;温度数据的采样周期为5min,即对应每个温度数据点为5min时刻的温度。
在步骤2中,进一步包括以下内容:
步骤2.1,对于步骤1所采集的历史数据,分别建立负荷功率时间序列和温度时间序列,其中,负荷功率时间序列为每1min一个时段的时间序列;温度时间序列为每5min一个时段的时间序列;
步骤2.2,对负荷功率时间序列进行降频处理,每隔5min取一个点,将功率与温度统一为5min一个时间段;
步骤2.3,求出每个时间点同一时刻后一天减前一天的温度增量以及对应的后一天的功率与前一天的功率的比值;从而得到温度增量序列与功率比值序列,二者对应起来形成数组,数组第一列为温度增量,第二列为功率比值。
在步骤3中,具体包括以下内容:
步骤3.1,对温度增量进行等级划分,按照温度增量划分为R个等级,即R个数据箱,所述R个等级是在温度增量最大值和最小值之间平均划分R个等级;
步骤3.2,根据每个历史数据数组中的“温度增量”的大小,将步骤2计算得到的所有数组分别存入其对应的温度增量等级数据箱中,完成数组的分箱过程。
优选地,所述R=7。
在步骤4中,根据各等级数据箱中的“温度增量”,对各数据箱分别绘制功率比值的概率密度直方图。
在步骤4中,具体包括以下内容:
步骤4.1,确定直方图的总宽度;由不同温度增量等级数据箱样本集合中的功率比值的最大值与最小值之差求得极差,直方图的总宽度(即组距与组数的乘积)取大于极差的数;
步骤4.2,确定直方图的组数、组距和各组界限,假设直方图的组数为N;每组的宽度为1/N,即组距;
步骤4.3,确定各组频数,把样本集合中的功率比值数据按照大小划分到直方图各组,统计各组中数据的个数,即为频数;
步骤4.4,绘制概率密度直方图,横坐标为功率比值,纵坐标为概率密度,直方图每组对应一个矩形,矩形的宽度为组距,高度为各组的概率密度,概率密度与频数的换算关系:概率密度=频数/总样本数*组数。
优选N=50。
在步骤5中,通用分布定义如下:
若连续型随机变量X服从一个形状参数为α、β和γ的通用分布,则记为:
X~V(α,β,γ)
其中,形状参数α、β和γ分别满足:
α>0,β>0,-∞<γ<+∞
通用分布的概率密度函数(Probability Density Function,PDF)f(x)表达式为:
通用分布的累积分布函数(Cumulative Distribution Function,CDF)F(x)表达式为:
F(x)=(1+e-α(x-γ))-β
CDF的反函数F-1(c)表达式为:
其中,f(x)为通用分布函数的概率密度,F(x)为通用分布函数的累计分布值,x为连续型随机变量X中的变量值,c表示F(x)。
在步骤6中,将待预测日各时间点的温度预报值,与当前日对应时间点的温度值做差得到温度增量序列,根据所对应所属的不同温度增量数据箱,利用步骤5所确定的不同温度增量数据箱所拟合得到的通用分布函数形状参数值,计算通用分布CDF逆函数F-1(c)的表达式,从而获取设定置信水平下的负荷预测功率的置信区间:
其中,wl,up和wl,low分别为负荷功率置信区间的上界和下界,conf为置信水平。
优选置信水平conf为0.95。
本申请还公开了一种基于前述负荷预测方法的负荷预测***,包括历史数据库、数据采集及处理单元、数据分箱处理单元、直方图单元、拟合计算单元、负荷功率预测计算单元;其特征在于:
所述数据采集及处理单元采集待预测负荷历史数据中的分时刻负荷功率及温度,并将负荷功率和温度采样值统一为相同的采样频率,建立包含温度增量序列与功率比值序列的数组;
所述数据分箱处理单元对数据采集及处理单元所建立的数组基于温度增量进行分箱处理;
所述直方图单元针对每一分箱数据绘制功率比值概率密度直方图,并计算直方图中各分组中的功率比值概率密度;
所述拟合计算单元采用通用分布函数模型对各等级数据箱内数据进行拟合,得到每一等级数据箱所对应的通用分布函数;
所述负荷功率预测计算单元根据输入的当前功率值和预测日的天气预报温度值,基于拟合计算单元拟合得到的各对应等级数据箱的通用分布函数,计算得到待预测负荷在预测日温度条件下的功率比值置信区间,实现给定气象下的功率预测。
相对于现有技术,本申请具有以下有益的技术效果。
通过概率建模的方法,建立负荷变化量与气象信息变化量的关联概率模型,采用直方图模型来近似表征不同温度增量条件下负荷功率变化量的实际分布,选取通用分布模型对实际分布直方图进行拟合,该模型能够有效表征分布规律的尖峰特性和偏轴特性且拟合精度高。利用该分布函数的累积概率密度函数可逆的数学性质,在确定的置信水平条件下能够解析求出负荷功率在已知气象预测结果下的增量区间,实现对电力***负荷功率预测结果加入气象影响因素后的动态修正。基于在不同日期的相同时刻下,温度增加了多少度,负荷功率相应地会增加多少倍,探究温度增量与功率比值之间的定量关系,具有创新性与可实现性,所建立的负荷变化量与温度变化量的关联概率模型,能够为分析温度变化对负荷功率的影响提供定量的数据分析方法,且该方法亦适用于分析其他外部因素如湿度、降雨量等对负荷功率变化的定量影响关系。
附图说明:
图1为本申请一种通用分布负荷预测方法流程示意图;
图2为本申请实施例第4个数据箱的概率分布直方图;
图3为本申请实施例对第4个箱内的“功率比值”进行拟合后的效果图;
图4为本申请一种通用分布负荷预测***结构框图。
具体实施方式:
下面结合说明书附图以及具体实施例对本发明的技术方案做进一步详细介绍。
如附图1所示为本申请公开的一种通用分布负荷预测方法流程示意图。所述包括以下步骤:
一种基于通用分布负荷预测方法,包括以下步骤:
步骤1,采集设定时间段内待预测负荷的分时刻负荷功率及温度的历史数据:
本申请优选采集待预测负荷在最近三个月内的全天24小时分时刻的负荷功率以及对应时刻的温度数据;其中,在历史数据中,负荷功率的采样周期为1min,即对应每个负荷功率数据点为1min时刻的有功功率;温度数据的采样周期为5min,即对应每个温度数据点为5min时刻的温度。
步骤2,对步骤1所采集负荷功率以及温度历史数据的进行处理,计算得到包含温度增量和功率比值的数组,其中,所述温度增量是指每个时间点同一时刻后一天减前一天的温度增量,所述功率比值是指每个时间点同一时刻后一天的功率与前一天的功率的比值;
步骤2.1,对于步骤1所采集的历史数据,分别建立负荷功率时间序列和温度时间序列,其中,负荷功率时间序列为每1min一个时段的时间序列;温度时间序列为每5min一个时段的时间序列;
步骤2.2,对负荷功率时间序列进行降频处理,每隔5min取一个点,将功率与温度统一为5min一个时间段;
步骤2.3,求出每个时间点同一时刻后一天减前一天的温度增量以及对应的后一天的功率与前一天的功率的比值;从而得到温度增量序列与功率比值序列,二者对应起来形成数组,数组第一列为温度增量,第二列为功率比值。
步骤3,对步骤2获得的包含温度增量和功率比值的数组进行分箱处理,所有数组分别存入其对应的温度增量等级数据箱中;
步骤3具体包括以下内容:
步骤3.1,对温度增量进行等级划分,按照温度增量划分为R个等级,即R个数据箱,所述R个等级是在温度增量最大值和最小值之间平均划分R个等级;
步骤3.2,根据每个历史数据数组中的“温度增量”的大小,将步骤2计算得到的所有数组分别存入其对应的温度增量等级数据箱中,完成数组的分箱过程。
其中,在本申请的优选方案中,优选R=7。
步骤4,建立基于直方图的各等级数据箱中负荷功率实际分布模型;
在步骤4中,具体包括以下内容:
步骤4.1,确定直方图的总宽度;将温度增量序列中的最大值减去最小值,得到差值,直方图总宽度取大于这个差值的最小整数;
步骤4.2,确定直方图的组数、组距和各组界限,假设直方图的组数为N;每组的宽度为1/N,即组距;
根据大量的统计结果以及计算分析,本申请选取N=50。
步骤4.3,确定各组频数,把样本集合中的功率比值数据按照大小划分到直方图各组,统计各组中数据的个数,即为频数;
步骤4.4,绘制概率密度直方图,横坐标为功率比值,纵坐标为概率密度,直方图每组对应一个矩形,矩形的宽度为组距,高度为各组的概率密度,概率密度与频数的换算关系:概率密度=频数/总样本数*组数。
步骤5,采用通用分布函数模型对各等级数据箱内数据进行拟合,得到每一等级数据箱所对应的通用分布函数;
通用分布定义如下:
若连续型随机变量X服从一个形状参数为α、β和γ的通用分布,则记为:
X~V(α,β,γ)
其中,形状参数α、β和γ分别满足:
α>0,β>0,-∞<γ<+∞
通用分布的概率密度函数(Probability Density Function,PDF)f(x)表达式为:
通用分布的累积分布函数(Cumulative Distribution Function,CDF)F(x)表达式为:
F(x)=(1+e-α(x-γ))-β
CDF的反函数F-1(c)表达式为:
其中,f(x)为通用分布函数的概率密度,F(x)为通用分布函数的累计分布值。
步骤6,根据待预测负荷的当前功率值和预测日的温度,基于步骤5拟合得到的通用分布函数,确定待预测负荷在预测日温度条件下的功率比值置信区间,实现给定气象下的功率预测。
在步骤6中,将待预测日各时间点的温度预报值,与当前日对应时间点的温度值做差得到温度增量序列,根据所对应所属的不同温度增量数据箱,利用步骤5所确定的不同温度增量数据箱所拟合得到的通用分布函数形状参数值,计算通用分布CDF逆函数F-1(c)的表达式,从而获取设定置信水平下的负荷预测功率的置信区间:
其中,wl,up和wl,low分别为负荷功率置信区间的上界和下界。
优选置信水平conf为0.95。
如附图4所示,本申请还公开了一种基于前述负荷预测方法的负荷预测***,包括历史数据库、数据采集及处理单元、数据分箱处理单元、直方图单元、拟合计算单元、负荷功率预测计算单元。
所述数据采集及处理单元采集待预测负荷历史数据中的分时刻负荷功率及温度,并将负荷功率和温度采样值统一为相同的采样频率,建立包含温度增量序列与功率比值序列的数组;所述数据分箱处理单元对数据采集及处理单元所建立的数组基于温度增量进行分箱处理;所述直方图单元针对每一分箱数据绘制功率比值概率密度直方图,并计算直方图中各分组中的功率比值概率密度;所述拟合计算单元采用通用分布函数模型对各等级数据箱内数据进行拟合,得到每一等级数据箱所对应的通用分布函数;所述负荷功率预测计算单元根据输入的当前功率值和预测日的天气预报温度值,基于拟合计算单元拟合得到的各对应等级数据箱的通用分布函数,计算得到待预测负荷在预测日温度条件下的功率比值置信区间,实现给定气象下的功率预测。
下面以某工业园区作为实施例进行介绍:
步骤1,采集设定时间段内待预测负荷的分时刻负荷功率及温度的历史数据在本实施例中,采集2016/11/1-2016/11/30;2017/7/1-2017/8/31一共三个月的全天24小时分时刻的负荷功率以及对应时刻的温度。其中,负荷功率的采样周期为1min,即对应每个负荷功率数据点为1min时刻的有功功率;温度数据的采样周期为5min,即对应每个温度数据点为5min时刻的温度。
步骤2,数据的处理:
步骤2.1,所收集到的数据中,功率为每1min一个时段的时间序列,对应每个数据点为1min时刻的有功功率;温度为每5min一个时段的时间序列,对应每个数据点为5min时刻的温度;每一天24小时,共有144个1min时间段和288个5min时间段,因此每一功率时间序列P共有144个数据,每一温度时间序列共有288个数据。
步骤2.2,对功率时间序列进行降频采样,每隔5min取一个点,将功率与温度统一为5min一个时间段;
步骤2.3,求出每个时间点同一时刻后一天减前一天的温度增量以及对应的后一天的功率与前一天的功率的比值;从而得到温度增量序列与功率比值序列,二者对应起来形成数组,数组第一列为温度增量,第二列为功率比值。
步骤3,基于温度增量预测等级的分箱处理:
步骤3.1,对温度增量进行等级划分:
按照温度增量划分为7个等级,即7个数据箱。将对应数组分别存入温度增量等级数据箱。在本实施例中,温度增量的取值范围在[-15.5℃,18.5℃]之内。将其分为7个箱,每个箱内的温度增量取值范围为[-16,-11],[-11,-6],[-6,-1],[-1,4],[4,9],[9,14],[14,19]。
步骤3.2,经过步骤2的数据处理后,历史数据已被划分为了以历史时刻为单位的数组;根据每个历史数据数组中的“温度增量”的大小,将数据库中所有数组分别存入其对应的温度增量等级数据箱中,完成数组的分箱过程。
步骤4,基于直方图的实际分布模型:
根据各等级数据箱中的“温度增量”样本,对各数据箱分别绘制功率比值的概率密度直方图;
步骤4.1,确定直方图的总宽度;将温度增量序列中的最大值减去最小值,得到差值,直方图总宽度取大于这个差值的最小整数;
步骤4.2,确定直方图的组数、组距和各组界限,假设直方图的组数为N=50;每组的宽度为1/N,即组距;
步骤4.3,确定各组频数,把样本集合中的功率比值数据按照大小划分到直方图各组,统计各组中数据的个数,即为频数;步骤4.4,绘制概率密度直方图,横坐标为功率比值,纵坐标为概率密度,直方图每组对应一个矩形,矩形的宽度为组距,高度为各组的概率密度,概率密度与频数的换算关系:概率密度=频数/总样本数*组数。
步骤5,箱内数据拟合:
采用通用分布函数模型对箱内数据进行拟合;
通用分布定义如下:
若连续型随机变量X服从一个形状参数为α、β和γ的通用分布,则记为:
X~V(α,β,γ)其中,形状参数α、β和γ分别满足:
α>0,β>0,-∞<γ<+∞
通用分布的概率密度函数(Probability Density Function,PDF)f(x)表达式为:
通用分布的累积分布函数(Cumulative Distribution Function,CDF)F(x)表达式为:
F(x)=(1+e-α(x-γ))-β
CDF的反函数F-1(c)表达式为:
对于不同的温度增量等级分别进行拟合,拟合求得的α,β和γ值分别如下表1所示:表1
直方图是实际分布,连续曲线是通用分布模型效果。如图2所示为本申请实施例第4个数据箱的概率分布直方图,图3为本申请实施例对第4个箱内的“功率比值”进行拟合后的效果图。
步骤6,根据待预测负荷的当前功率值和预测日的温度,基于步骤5拟合得到的通用分布函数,确定待预测负荷在预测日温度条件下的功率比值置信区间,实现给定气象下的功率预测。
求取一定置信水平(以95%为例)下功率比值的最大值与最小值,以及给定气象下的功率预测。假定已知当前的功率值为10,可以计算得出95%置信水平下明天的功率最小值与最大值如下表2所示:
表2
直方图是实际分布,连续曲线是通用分布模型效果。如图2所示为本申请实施例第4个数据箱的概率分布直方图,图3为本申请实施例对第4个箱内的“功率比值”进行拟合后的效果图。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。
Claims (12)
1.一种基于通用分布的负荷预测方法,其特征在于,所述负荷预测方法包括以下步骤:
步骤1,采集设定时间段内待预测负荷的分时刻负荷功率及温度的历史数据;
步骤2,对步骤1所采集负荷功率以及温度历史数据进行处理,计算得到包含温度增量和功率比值的数组,其中,所述温度增量是指每个时间点同一时刻后一天温度减前一天温度所得的温度增量,所述功率比值是指每个时间点同一时刻后一天的功率与前一天的功率的比值;
步骤3,对步骤2获得的包含温度增量和功率比值的数组进行分箱处理,即将所有数组分别存入其对应的温度增量等级数据箱中;
步骤4,基于直方图,建立各等级数据箱中负荷功率实际分布模型;
步骤5,采用通用分布函数模型对各等级数据箱内数据进行拟合,得到每一等级数据箱所对应的通用分布函数;
步骤6,根据待预测负荷的当前功率值和预测日的温度,基于步骤5拟合得到的通用分布函数,确定待预测负荷在预测日温度条件下的功率比值置信区间,实现给定气象下的功率预测。
2.根据权利要求1所述的负荷预测方法,其特征在于:
在步骤1中,采集待预测负荷在最近三个月内的全天24小时分时刻的负荷功率以及对应时刻的温度数据;其中,在历史数据中,负荷功率的采样周期为1min,即对应每个负荷功率数据点为1min时刻的有功功率;温度数据的采样周期为5min,即对应每个温度数据点为5min时刻的温度。
3.根据权利要求1所述的负荷预测方法,其特征在于:
在步骤2中,进一步包括以下内容:
步骤2.1,对于步骤1所采集的历史数据,分别建立负荷功率时间序列和温度时间序列,其中,负荷功率时间序列为每1min一个时段的时间序列;温度时间序列为每5min一个时段的时间序列;
步骤2.2,对负荷功率时间序列进行降频处理,每隔5min取一个点,将功率与温度的时间序列统一为5min一个时段;
步骤2.3,求出每个时间点同一时刻后一天温度减前一天温度所得的温度增量以及对应的后一天的功率与前一天的功率所得的功率比值;从而得到温度增量序列与功率比值序列,二者对应起来形成数组,数组第一列为温度增量,第二列为功率比值。
4.根据权利要求1所述的负荷预测方法,其特征在于:
在步骤3中,具体包括以下内容:
步骤3.1,对温度增量进行等级划分,按照温度增量划分为R个等级,即R个数据箱,所述R个等级是指在温度增量最大值和最小值之间平均划分的R个等级;
步骤3.2,根据每个历史数据数组中的温度增量的大小,将步骤2计算得到的所有数组分别存入其对应的温度增量等级数据箱中,完成数组的分箱过程。
5.根据权利要求4所述的负荷预测方法,其特征在于:
所述R=7。
6.根据权利要求1所述的负荷预测方法,其特征在于:
在步骤4中,根据各等级数据箱中的温度增量,对各数据箱分别绘制功率比值的概率密度直方图。
7.根据权利要求1所述的负荷预测方法,其特征在于:
在步骤4中,具体包括以下内容:
步骤4.1,确定直方图的总宽度;将温度增量序列中的最大值减去最小值,得到差值,直方图的总宽度(即组距与组数的乘积)取大于这个差值的最小整数;
步骤4.2,确定直方图的组数、组距和各组界限,假设直方图的组数为N;每组的宽度为1/N,即组距;
步骤4.3,确定各组频数,把样本集合即分箱处理的所有数组中的功率比值数据按照大小划分到直方图各组,统计各组中数据的个数,即为频数;
步骤4.4,绘制概率密度直方图,横坐标为功率比值,纵坐标为概率密度,直方图每组对应一个矩形,矩形的宽度为组距,高度为各组的概率密度,概率密度与频数的换算关系:概率密度=频数/总样本数*组数。
8.根据权利要求7所述的负荷预测方法,其特征在于:
直方图的组数N=50。
11.根据权利要求10所述的负荷预测方法,其特征在于:
置信水平conf为0.95。
12.一种利用前述权利要求1-11任一项权利要求所述负荷预测方法的负荷预测***,包括历史数据库、数据采集及处理单元、数据分箱处理单元、直方图单元、拟合计算单元、负荷功率预测计算单元;其特征在于:
所述数据采集及处理单元采集待预测负荷历史数据中的分时刻负荷功率及温度,并将负荷功率和温度采样值统一为相同的采样频率,建立包含温度增量序列与功率比值序列的数组;
所述数据分箱处理单元对数据采集及处理单元所建立的数组基于温度增量进行分箱处理;
所述直方图单元针对每一分箱数据绘制功率比值概率密度直方图,并计算直方图中各分组中的功率比值概率密度;
所述拟合计算单元采用通用分布函数模型对各等级数据箱内数据进行拟合,得到每一等级数据箱所对应的通用分布函数;
所述负荷功率预测计算单元根据输入的当前功率值和预测日的天气预报温度值,基于拟合计算单元拟合得到的各对应等级数据箱的通用分布函数,计算得到待预测负荷在预测日温度条件下的功率比值置信区间,实现给定气象下的功率预测。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010322604.3A CN111552923B (zh) | 2020-04-22 | 2020-04-22 | 一种基于通用分布的负荷预测方法及负荷预测*** |
PCT/CN2021/086342 WO2021213192A1 (zh) | 2020-04-22 | 2021-04-11 | 一种基于通用分布的负荷预测方法及负荷预测*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010322604.3A CN111552923B (zh) | 2020-04-22 | 2020-04-22 | 一种基于通用分布的负荷预测方法及负荷预测*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111552923A CN111552923A (zh) | 2020-08-18 |
CN111552923B true CN111552923B (zh) | 2021-05-04 |
Family
ID=71998399
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010322604.3A Active CN111552923B (zh) | 2020-04-22 | 2020-04-22 | 一种基于通用分布的负荷预测方法及负荷预测*** |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN111552923B (zh) |
WO (1) | WO2021213192A1 (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110689169B (zh) * | 2019-09-04 | 2023-05-30 | 上海工程技术大学 | 一种基于分数Levy稳定运动模型的短期电力负荷预测方法 |
CN111552923B (zh) * | 2020-04-22 | 2021-05-04 | 国网江苏省电力有限公司苏州供电分公司 | 一种基于通用分布的负荷预测方法及负荷预测*** |
CN113744082A (zh) * | 2021-08-25 | 2021-12-03 | 国网(苏州)城市能源研究院有限责任公司 | 一种综合刻画源荷两侧不确定性的二维分箱方法 |
CN114037282B (zh) * | 2021-11-10 | 2024-06-21 | 国网江苏省电力有限公司南京供电分公司 | 一种基于用电特征的用电量标签*** |
CN114142472B (zh) * | 2021-12-06 | 2023-08-08 | 浙江华云电力工程设计咨询有限公司 | 基于混合高斯分布概率密度的风光容量配置方法及*** |
CN114254838B (zh) * | 2022-01-07 | 2024-06-28 | 深圳供电局有限公司 | 一种短期电力负荷预测影响因子的确定方法 |
CN114322321B (zh) * | 2022-01-20 | 2023-05-16 | 中国石油大学(华东) | 一种燃气热水器的灰色预测可拓恒温控制方法 |
CN115130778B (zh) * | 2022-07-19 | 2024-06-25 | 国网安徽省电力有限公司滁州供电公司 | 基于大数据的一体化电力调度自动化*** |
CN115313377B (zh) * | 2022-08-26 | 2023-08-08 | 中国长江三峡集团有限公司 | 一种电力负荷预测方法及*** |
CN116245221B (zh) * | 2023-01-09 | 2024-03-08 | 上海玫克生储能科技有限公司 | 一种负荷实时预测方法、装置和电子设备 |
CN116204805B (zh) * | 2023-04-24 | 2023-07-21 | 青岛鑫屋精密机械有限公司 | 一种微压氧舱及数据管理*** |
CN116258355B (zh) * | 2023-05-15 | 2023-08-11 | 国网浙江省电力有限公司永康市供电公司 | 适用于多点功率估计的配电台区负荷曲线分解方法及装置 |
CN116865246B (zh) * | 2023-06-27 | 2023-12-26 | 广东电网有限责任公司广州供电局 | 一种基于快速响应的工业用户负荷可行域预测方法及*** |
CN117150891B (zh) * | 2023-08-15 | 2024-04-26 | 幂光新材料科技(上海)有限公司 | 基于数据驱动的led灯珠功率智能预测方法及*** |
CN117495434B (zh) * | 2023-12-25 | 2024-04-05 | 天津大学 | 电能需求预测方法、模型训练方法、装置及电子设备 |
CN117852928B (zh) * | 2024-03-08 | 2024-06-04 | 国网北京市电力公司 | 一种近零能耗建筑负荷预测方法、装置、设备及介质 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7587330B1 (en) * | 2003-01-31 | 2009-09-08 | Hewlett-Packard Development Company, L.P. | Method and system for constructing prediction interval based on historical forecast errors |
KR101012863B1 (ko) * | 2008-09-25 | 2011-02-08 | 한국전력공사 | 고객기준부하 산출을 위한 부하 예측 비교 분석 시스템 |
CN102930356B (zh) * | 2012-11-16 | 2016-06-22 | 广东电网公司电力调度控制中心 | 基于气象因素敏感度的短期负荷预测方法 |
CN103268524B (zh) * | 2013-06-03 | 2017-03-29 | 国家电网公司 | 一种提高电网短期负荷预测准确度的方法 |
CN104463344B (zh) * | 2014-10-29 | 2018-01-16 | 广东电网有限责任公司电力调度控制中心 | 电网短期负荷预测方法和*** |
CN106408223A (zh) * | 2016-11-30 | 2017-02-15 | 华北电力大学(保定) | 基于气象相似日及误差校正的短期负荷预测方法 |
CN106971240A (zh) * | 2017-03-16 | 2017-07-21 | 河海大学 | 一种变量选择与高斯过程回归的短期负荷预测方法 |
CN111008727A (zh) * | 2019-10-28 | 2020-04-14 | 许昌许继软件技术有限公司 | 一种配电台区负荷预测方法及装置 |
CN111552923B (zh) * | 2020-04-22 | 2021-05-04 | 国网江苏省电力有限公司苏州供电分公司 | 一种基于通用分布的负荷预测方法及负荷预测*** |
-
2020
- 2020-04-22 CN CN202010322604.3A patent/CN111552923B/zh active Active
-
2021
- 2021-04-11 WO PCT/CN2021/086342 patent/WO2021213192A1/zh active Application Filing
Also Published As
Publication number | Publication date |
---|---|
CN111552923A (zh) | 2020-08-18 |
WO2021213192A1 (zh) | 2021-10-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111552923B (zh) | 一种基于通用分布的负荷预测方法及负荷预测*** | |
CN111260136A (zh) | 一种基于arima-lstm组合模型的楼宇短期负荷预测方法 | |
CN109726865A (zh) | 基于emd-qrf的用户负荷概率密度预测方法、装置和存储介质 | |
CN106600037B (zh) | 一种基于主成分分析的多参量辅助负荷预测方法 | |
CN111680841B (zh) | 基于主成分分析的短期负荷预测方法、***及终端设备 | |
CN114372360A (zh) | 用于电力负荷预测的方法、终端及存储介质 | |
CN113837488B (zh) | 能源消耗数据预测方法、***、设备 | |
WO2022021727A1 (zh) | 城市综合体用电量预测方法、装置、电子设备及存储介质 | |
CN112257962B (zh) | 一种台区线损预测方法及装置 | |
CN110751416A (zh) | 一种用水量的预测方法、装置及设备 | |
CN113177366A (zh) | 一种综合能源***规划方法、装置和终端设备 | |
CN115545333A (zh) | 一种多负荷日类型配电网负荷曲线预测方法 | |
CN115187134A (zh) | 基于网格的配电网规划方法、装置及终端设备 | |
CN115860797A (zh) | 一种适应新电价改革形势的电量需求预测方法 | |
CN115829235A (zh) | 一种基于大数据分析的用电调度方法 | |
CN113902181A (zh) | 公变重过载的短期预测方法及设备 | |
CN113139698A (zh) | 负荷预测方法、装置及设备 | |
CN110414776B (zh) | 分行业用电特性快速响应分析*** | |
CN112330089A (zh) | 装备制造类企业的综合能效监测方法及监测*** | |
CN116826745A (zh) | 一种电力***背景下的分层分区短期负荷预测方法及*** | |
CN113469440B (zh) | 一种预测月尺度生活需水量的方法 | |
CN116050636A (zh) | 光伏电站的出力预测方法、装置、设备及介质 | |
Ignatiadis et al. | Forecasting residential monthly electricity consumption using smart meter data | |
CN114282658A (zh) | 一种针对流量序列分析及预测的方法、设备和介质 | |
Kien et al. | Application of Sarima model in load forecasting in Hanoi city |
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 |