CN104156783B - 计及气象累积效应的电力系统最大日负荷预测系统及方法 - Google Patents
计及气象累积效应的电力系统最大日负荷预测系统及方法 Download PDFInfo
- Publication number
- CN104156783B CN104156783B CN201410367888.2A CN201410367888A CN104156783B CN 104156783 B CN104156783 B CN 104156783B CN 201410367888 A CN201410367888 A CN 201410367888A CN 104156783 B CN104156783 B CN 104156783B
- Authority
- CN
- China
- Prior art keywords
- load
- temperature
- model
- meteorological
- day
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000000694 effects Effects 0.000 title claims abstract description 26
- 230000001186 cumulative effect Effects 0.000 claims abstract description 53
- 238000009825 accumulation Methods 0.000 claims abstract description 17
- 238000012937 correction Methods 0.000 claims abstract description 17
- 238000012216 screening Methods 0.000 claims abstract description 12
- 230000008569 process Effects 0.000 claims abstract description 9
- 230000035945 sensitivity Effects 0.000 claims description 26
- 238000004458 analytical method Methods 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 11
- 240000002853 Nelumbo nucifera Species 0.000 claims description 8
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 8
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 8
- 238000000611 regression analysis Methods 0.000 claims description 7
- 230000005611 electricity Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000035772 mutation Effects 0.000 claims description 4
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 238000002715 modification method Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 2
- 238000012417 linear regression Methods 0.000 claims description 2
- 230000001932 seasonal effect Effects 0.000 claims description 2
- 238000000205 computational method Methods 0.000 claims 1
- 238000000605 extraction Methods 0.000 claims 1
- 239000008246 gaseous mixture Substances 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 102000003712 Complement factor B Human genes 0.000 description 1
- 108090000056 Complement factor B Proteins 0.000 description 1
- 241000406668 Loxodonta cyclotis Species 0.000 description 1
- 238000004378 air conditioning Methods 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 238000013277 forecasting method Methods 0.000 description 1
- 238000009472 formulation 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
- 238000003012 network analysis Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种计及气象累积效应的电力系统最大日负荷预测系统及方法,由数据采集模块将气象、负荷数据读入系统;通过数据筛分处理模块将负荷数据分为不含累积效应部分及累积效应显著部分,并确定影响电网负荷的关键气象指数;在负荷预测建模模块中,首先建立非累积日负荷基础预测模型,在此基础上建立累积日负荷修正模型,由此生成综合预测模型;最后在负荷预测模块中,将建模模块结果与气象及负荷数据库相联结获得负荷预测值输出到软件图形界面,提供给用户。本发明能够反映气象条件对电网负荷的时间性和累积性影响,帮助电网运行人员准确地掌握负荷的变动规律,及时根据天气预报的气象情况预测最大日负荷变化,提高短期负荷预测的精度。
Description
技术领域
本发明涉及电力系统考虑气象因素的中短期负荷预测方法技术,是一种考虑气象因素综合影响的电力系统最大日负荷预测建模优化系统及其方法。
背景技术
电力系统最大日负荷预测短期负荷预测是电网调度的重要环节,是电力系统中短期负荷预测领域的重要组成部分,对电力系统节能减排等都有重要影响。准确的负荷预测对于保持电力系统的安全稳定运行、保障人们生产活动和生活的有序进行具有重要的意义。随着居民生活条件的逐步改善,空调、暖气高频率地出现在日常生活工作中,电网统调负荷快速增长,调温负荷所占的比例逐年增大,从而构成用电峰荷,拉大电网负荷峰谷差。电网气象敏感负荷不断上升,气候对负荷特性的影响越来越大,还存在着复杂的累积和滞后等效应。这些因素使电网负荷规律性变得很复杂,给负荷建模、负荷管理、电网规划带来许多困难。
发明内容
本发明将根据电网实际的负荷、经济、各种气象因素等历史数据,系统分析负荷特性,对电网统调日最大负荷及气象因素进行分析研究,提出了一种计及气象累积效应的电力系统最大日负荷预测模型及优化方法,包括不同气象指数对最大日负荷的交叉影响分析方法、基于地市电网负荷形成省级多因素混合气象指数分析方法、气象指数与负荷关系的灵敏度实用化建模、温度累积效应的建模和优化方法等。通过该方法能够反映气象条件对电网负荷的时间性和累积性影响,帮助电网运行人员准确地掌握负荷的变动规律,及时根据天气预报的气象情况预测最大日负荷变化,提高短期负荷预测的精度,完成日调度计划的制定。
本发明解决上述问题所采用的技术方案是:一种计及气象累积效应的电力系统最大日负荷预测系统,包括数据采集模块、数据筛分处理模块、负荷预测建模模块、负荷预测模块和软件图形界面模块。
数据采集模块是将外部气象、负荷数据传递到系统中进行解析;数据筛分处理模块是将采集来的数据进行季节性筛选,并且将负荷分解为累积效应显著负荷和累积效应不显著负荷两大类分别存储;负荷预测建模模块,包括非累积日负荷基础预测模型建模模块及累积日负荷修正模型建模模块,非累积日负荷基础预测模型建模模块是提取累积效应不显著的负荷进行正常工作日模型的建立和预测,累积日负荷修正模型建模模块是在非累积日负荷基础预测模型基础上针对累积日进行建模和预测修正,最后建立综合预测模型;负荷预测模块是将负荷预测建模模块结果与气象及负荷数据库相联接获取负荷预测值;软件图形界面模块是对以上预测结果的输出。
所述的计及气象累积效应的电力系统最大日负荷预测系统的预测方法,包括以下步骤:
首先采用相关度分析、关联误差分析等方法进行特征选择,确定影响电网负荷的关键气象指数;通过回归分析和灵敏度计算方法得到具体电网负荷对关键气象指数的灵敏度,并通过最小二乘法拟合逼近得到气象指数与负荷关系的灵敏度实用化表达方式;通过数据筛选分离出累积效应不显著负荷并建立基础预测模型,在此模型上对累积效应显著负荷建模,采用不同气象因素预测负荷的误差交叉修正方法,利用最优化方法求解最优修正因子,得到考虑了温度累积效应的实际负荷预测工程化模型,并借助灰色关联度、预测误差率等评估修正结果。
本发明的工作原理及过程如下:
(1)利用负荷、各种气象因素等历史数据,计算多因素混合气象指数;
(2)分析不同气象因素包括温度、湿度、风速、降雨量以及四个多因素混合气象指数(温湿指数、实感温度、舒适度、寒湿指数)对负荷的交叉作用,采用相关性分析方法、关联误差分析等方法,分析影响最大日负荷的各种气象因素及其相关度,通过多种特征选择方法确定影响最大日负荷的主要气象因素;
(3)根据地市多种气象指数和地市负荷的相关性,建立地市多种气象指数和地市负荷推导模型,基于地市电力负荷和地市气象指数计算省级多因素混合气象指数;
(4)筛分累积效应显著与非显著日,获取这两种情况下的气象数据和负荷数据,以方便对正常负荷和含累积效应负荷建模;
(5)根据步骤(2)得到的影响最大日负荷的气象因素,构建气象指数与负荷关系的灵敏度实用化建模,得到不含温度累积效应时的正常负荷预测工程化模型;
(6)将实际负荷与上述正常负荷预测值之差当作由累积效应引起的负荷偏差,利用不同于正常负荷预测的气象指数对该负荷偏差建立优化模型;
(7)求解由累积效应引起负荷偏差建立优化模型,获取修正模型系数;
(8)根据实际负荷值与含累积效应的负荷预测值,计算预测误差率,并用灰色关联度评估修正结果,获取最优的气象因素交叉修正方案。
与现有技术方案相比,本发明的有益效果是:
建立了考虑温度累积效应影响的最大日负荷预测模型,根据历年负荷和气象数据,分析了电网最大日负荷与气象之间的关系及影响规律,建立的模型以气象指数作为输入量,根据气象部门的天气预报情况判断未来负荷变化趋势。所建模型考虑了气象因素对负荷影响的复杂性,在建立负荷预测模型时应该将正常的负荷与受累积效应影响显著的负荷分开来考虑,这样在数学层面上避免了“重复”运算。考虑了单一气象因素和多因素混合气象因素,利用不同的气象因素及多因素混合气象指数交叉修正预测负荷,体现了不同气象条件对负荷的交叉影响。提出基于地级电用电情况的省级地区多因素混合气象指数计算方法,处理各地市间由于地域原因产生的气象差异。提出的灵敏度简化表达方式简化了模型的复杂程度,使其在工程上更容易实现,同时减轻了预测人员的繁杂度,提高了预测工作的效率。由于累积效应最直接影响的电量,因此该模型处理累积效应时采用直接修正负荷的方法,而不采用修正气象指数的方法,避免了误差的多级放大。建立优化模型,采用优化方法,有效减小了预测误差,提高了电网短期负荷预测的准确度,为合理安排发电计划提供依据。
附图说明
图1为本发明所述的计及气象累积效应的电力系统最大日负荷预测系统的结构连接图。
图2为本发明所述的多气象因素交叉修正方法。
图3为本发明所述的非累积日负荷基础预测模型建模模块的整体实现过程。
图4为本发明所述的累积日负荷修正模型建模模块的整体实现过程。
图5为本发明实施例所述的夏季最大负荷与综合最高温度的曲线图。
图6表示本发明实施例所述的最高温度与负荷关系灵敏度实用化表达及相关性。
图7表示本发明实施例所述的负荷偏差(ΔL)与温差(ΔT1)的关系。
图8表示本发明实施例所述的负荷偏差(△L)与温差(△T2)的关系。
图9表示本发明实施例所述的实际负荷值与最高温度初步预测值之间的关系。
图10表示本发明实施例所述的实际负荷值与平均实感温度二次拟合修正预测的负荷值的相关关系。
具体实施方式
下面参照附图,详细叙述本发明的具体实施方案。
如图1所示,本发明所述的计及气象累积效应的电力系统最大日负荷预测系统,包括数据采集模块、数据筛分处理模块、负荷预测建模模块、负荷预测模块和软件图形界面模块。
本发明所述的电力系统最大日负荷预测建模中考虑气象累积效应影响的建模方法,包括影响负荷变动的关键气象因素选择,建立不含累积效应的基础预测模型,对累积效应显著的负荷进行预测修正建模。
如图2、图3和图4所示,建立考虑气象因素综合影响的电力系统最大日负荷预测模型,包括如下步骤:
1、在影响负荷变动的关键气象因素选择模块中,分析气象对负荷的交叉作用,采用计算负荷与各气象指数的相关度以及方差分析等方法,对各影响因子进行比较,相关度大并且方差小的气象指数即确认是关键的影响因子。
在大多数研究中认为,各个气象因素直接作用于电力负荷。实际上,气象因素对于电力负荷的影响规律是非常复杂的,而且往往存在着不同气象因素的交互影响,这就需要进一步分析多个气象因素产生的耦合效果(多因素混合气象指数)及其对电力系统的影响规律。
在影响短期负荷的气象因素中,温度很明显是影响最大的一个,低温和高温闷热可以引起相对较大的高峰负荷。事实上,气象因素对负荷的影响是通过改变人体对环境的舒适度感觉而实现的。为此,本发明引入生物气象学中的温湿指数、寒湿指数、实感温度、舒适度指数4个指标来综合衡量气象因素对电力负荷的影响,并用多因素混合气象指数与温度、湿度等单一因素的相关关系进行比较,找出对负荷特征的变化有更好的跟随性和描述效果的气象指数。
相关性分析是指对两个或多个具备相关性的变量元素进行分析,从而衡量两个变量因素的相关密切程度。相关性的元素之间需要存在一定的联系或者概率才可以进行相关性分析。本方法计算各气象因素以及各多因素混合气象指数与电网还原负荷的相关系数,找出对负荷的影响大且影响较为稳定的因素,并用回归分析法建立关系方程。利用方差分析等方法,对各影响因子进行比较,相关度大并且方差小的气象指数即认为是关键的影响因子。将不同的影响因子分别用在正常负荷建模和累积效应修正建模上,生成交叉修正方案。
以某地区实际电网为例,对2008年到2013年夏季最大日负荷和气象数据进行分析,负荷与各气象之间的关联程度如表1所示,从表内可以看出最大温度及平均实感温度的相关度最高且各年相关度的方差较小,说明这两个因素在该电网区域,对夏季最大日负荷的影响相对较大且影响较为稳定。
表1 某电网夏季日最大日负荷与加权多因素混合气象的相关度
因此,夏季最大日负荷预测计算时,选择的关键的影响因子可为最大温度和平均实感温度。
2、基于地市电力负荷计算省级多因素混合气象指数。
由于地域的原因,各市的气象条件存在一定的差别,在分析统调总负荷与气象的关系时,若只用某市的气象数据来分析必定会对准确性造成一定影响。用各市的实际气象,结合各市的实际用电情况形成地区加权多因素混合气象因素,例如加权综合温度,在这个温度基础上进行负荷与温度的关系研究将会更具科学性。
以某地区综合温度为例,本发明将采取如下形式的公式计算某地区综合温度:
其中,N表示地区个数,Pi表示各市的负荷,P∑表示整个地区的总负荷,ωi是各地负荷占总调负荷的权重,T是形成的地区综合温度。
下面我们以2011年为例:得到2011年夏季8月1日-8月31日综合温度与电网最大负荷的曲线,如图5所示。
3、筛分累积效应显著与非显著数据
根据累积效应的显著程度,需要对现有数据进行分离,筛选出累积效应显著的数据和不显著的数据。首先我们留下夏季累积效应不显著的数据,其定义如下:
1)留下4-9月数据;
2)留下28-38度以外的数据;
3)去掉节假日及周末数据;
4)去掉拉闸限电日数据;
5)去掉降雨量超过10mm的数据;
6)去掉气象指数突升或突降日开始2日(T-2,T-1,T0,T1,T2)的数据;
其次,留下夏季累积效应显著的数据,即原始数据与不显著数据的差集。
4、建立不含累积效应的基础预测模型模块,建模方案如下:
以温度为例,对夏季正常负荷进行预测建模。将负荷分为经济负荷和气象负荷,其中气象负荷与温度存在关联性,通过历史样本拟合温度与负荷的灵敏度;经济负荷则通过历史样本,推算线性增长关系。两者叠加求得未来某时间段的负荷。
(1)根据每年的历史样本,对夏季累积效应不显著统调负荷需求与最高温度进行一元二次回归分析,做出拟合曲线,求得每年的拟合函数及拟合度;
(2)根据拟合得到的拟合函数f(x),对f(x)求一阶导数可得各年统调负荷需求对最高温度的灵敏度。将N个温度区间分别定义为[Ti1-Ti2],i=1,2…N,则第i个温度区间的灵敏度计算公式为:
Sj=f(Ti1-Ti2) (4-1)
(3)计算这几年相同温度区间的灵敏度平均值,取这几年每个温度区间的灵敏度平均值作为新一年夏季统调负荷需求对温度的灵敏度,再对拟合温度与该灵敏度平均值进行一元一次回归分析,做出拟合曲线,求出各温度对应的灵敏度关系式;
(4)当温度变化ΔT(℃)时,负荷的变化量为:
其中,ΔT=T0-Tr,T0表示预测日的温度,Tr表示参考日的温度。
(5)随着经济的发展,每年夏季的经济负荷均有一定的增长量,此时算出每月的负荷平均增长量ΔLtime。
(6)同时考虑负荷随时间的变化量和负荷随温度的变化量,最终正常负荷的预测公式为:
L′=Lr+ΔL+ΔLtime (4-3)
其中,Lr为参考日的地区最大日负荷。
以某地区2008年-2013年夏季数据为例,表2是地区最大日负荷对最高温度的灵敏度,图6是最高温度与负荷关系的灵敏度实用化表达。
表2 统调最大日负荷对最高温度的灵敏度(单位:MW/℃)
温度区间 | 2008 | 2009 | 2010 | 2011 | 2012 | 2013 | 平均灵敏度 |
[20,21] | -18.48 | 5.11 | -33.81 | -70.17 | 73.26 | -91.53 | -22.60 |
[21,22] | -8.31 | 34.51 | -19.06 | -33.00 | 77.23 | -60.30 | -1.49 |
[22,23] | 1.86 | 63.92 | -4.30 | 4.16 | 81.20 | -29.07 | 19.63 |
[23,24] | 12.02 | 93.32 | 10.46 | 41.32 | 85.17 | 2.16 | 40.74 |
[24,25] | 22.19 | 122.72 | 25.22 | 78.49 | 89.15 | 33.39 | 61.86 |
[25,26] | 32.36 | 152.12 | 39.98 | 115.65 | 93.12 | 64.61 | 82.97 |
[26,27] | 42.53 | 181.52 | 54.74 | 152.82 | 97.09 | 95.84 | 104.09 |
[27,28] | 52.69 | 210.93 | 69.50 | 189.98 | 101.06 | 127.07 | 125.20 |
[28,29] | 62.86 | 240.33 | 84.26 | 227.14 | 105.03 | 158.30 | 146.32 |
[29,30] | 73.03 | 269.73 | 99.01 | 264.31 | 109.01 | 189.53 | 167.44 |
[30,31] | 83.20 | 299.13 | 113.77 | 301.47 | 112.98 | 220.75 | 188.55 |
[31,32] | 93.37 | 328.53 | 128.53 | 338.64 | 116.95 | 251.98 | 209.67 |
[32,33] | 103.53 | 357.94 | 143.29 | 375.80 | 120.92 | 283.21 | 230.78 |
[33,34] | 113.70 | 387.34 | 158.05 | 412.96 | 124.89 | 314.44 | 251.90 |
[34,35] | 123.87 | 416.74 | 172.81 | 450.13 | 128.86 | 345.67 | 273.01 |
[35,36] | 134.04 | 446.14 | 187.57 | 487.29 | 132.84 | 376.89 | 294.13 |
[36,37] | 144.20 | 475.54 | 202.33 | 524.46 | 136.81 | 408.12 | 315.24 |
[37,38] | 154.37 | 504.95 | 217.09 | 561.62 | 140.78 | 439.35 | 336.36 |
[38,39] | 164.54 | 534.35 | 231.84 | 598.78 | 144.75 | 470.58 | 357.47 |
[39,40] | 174.71 | 563.75 | 246.60 | 635.95 | 148.72 | 501.81 | 378.59 |
5、对累积效应显著的负荷进行预测修正建模,建模方案如下:
设定正常负荷预测模型对正常负荷预测是精准的,那么当累积效应显著时,用该模型计算出的预测负荷与实际负荷的偏差,可以认为主要是由累积效应引起的。因此,对累积效应进行修正负荷建模,实际上是对以上负荷偏差进行建模,则实际负荷即正常负荷与由累积效应引起的负荷偏差的叠加。
以温度为例,建立夏季累积效应修正负荷模型:
计算所有温度差值:
ΔT1=T0-T-1 (5-1)
及
ΔT2=T0-T-2 (5-2)
其中,T0是气温突变日的气温,T-1、T-2分别是突变日前一天和两天的气温,该公式的结果作为累积效应修正公式的输入变量,反映了温度变化与人们感知惯性造成的偏差,与累积效应的物理现象相符合。
以L-2为基准,按正常负荷预测模型算出L0′、L1′、L2′,并计算所有负荷偏差值ΔL:
ΔL=L-L′ (5-3)
利用残差最小化的思想,通过最小二乘法进行多元回归分析,得到累积效应修正量ΔL′。
二元线性回归的表达式为:
ΔL′=f(ΔT1,ΔT2)=k1ΔT1+k2ΔT2+k3 (5-4)
二元二次回归的表达式为:
ΔL′=f1(ΔT1,ΔT2)=k1ΔT1 2+k2ΔT2 2+k3ΔT1ΔT2+k4ΔT1+k5ΔT2+k6 (5-5)
其中,k1,k2,k3,k4,k5,k6为待求的优化系数。综合考虑正常负荷模型与累积效应修正负荷模型,得到最终的实际负荷模型为:
L″=Lr+ΔLT+ΔLtime+f(ΔT1,ΔT2)=L″=L′+ΔL′ (5-6)
可建立优化模型:
若选取两项对电网短期负荷影响大且影响较为稳定的气象因素A和B,一项用于预测正常负荷L′,一项用于累积效应修正量ΔL′,则正常负荷模型及累积效应修正负荷模型中可用变量的组合有四种,如表3所示,本发明按这四种方式进行建模分析,找出其中最佳的方案。
表3 两种气象因素建模的组合方法
L′ | f | |
方案一 | A | A |
方案二 | A | B |
方案三 | B | A |
方案四 | B | B |
气象因素A和B的选择可根据步骤1所求出的关键的影响因子确定。
1)方案一模型:
(1)f(ΔA1,ΔA2)为线性方程:
(2)f(ΔA1,ΔA2)为二次方程:
2)方案二模型:
(1)f(ΔB1,ΔB2)为线性方程:
(2)f(ΔB1,ΔB2)为二次方程:
3)方案三模型:
(1)f(ΔA1,ΔA2)为线性方程:
(2)f(ΔA1,ΔA2)为二次方程:
4)方案四模型:
(1)f(ΔB1,ΔB2)为线性方程:
(2)f(ΔB1,ΔB2)为二次方程:
其中,a、b、k1、k2、k3、k4、k5、k6∈R为待优化系数,可利用现代内点法对该优化模型进行求解,以求出优化系数。
下面我们以方案一为例,关键气象因素A选择最高气温。利用4-3的初步预测模型对筛选出来的累积日进行预测,得到初步预测的负荷为P和预测偏差值ΔL。负荷偏差(ΔL)与ΔT1的关系以及负荷偏差(ΔL)与ΔT2的关系如图7、图8所示。由图7、图8可以看出,按原不含累积效应的基础预测模型算出的负荷与实际负荷的差值,在一定温度区间内与温度差呈现明显的负相关性,体现了负荷随气象因素变化的累积效应,即如果温度由持续高温向低温突变,则负荷的减少程度不明显;反之,如果温度由持续低温向高温突变,则负荷的上升程度不明显。
本发明对预测结果进行分析,采用方案一和方案三分别某电网2008年-2013年夏季,气象累积日进行预测,预测结果比较如表4所示,其中关键气象因素A选择最高气温,关键气象因素B选择平均实感温度。
表4 方案一和方案三的预测结果
由表4可以直观地看出加了累积效应修正后模型的整体预测精度较未考虑累积效应的基础预测模型有较大提高,同时方案三在整体上比方案一所得的模型预测结果更好,这说明考虑综合的平均实感温度要比单一考虑最高温度的预测模型更为精准,并且二次拟合修正模型预测结果比线性拟合修正模型的预测结果更好。图9为实际负荷值与最高温度初步预测值之间的关系。图10为实际负荷值与平均实感温度二次拟合修正预测的负荷值的相关关系。
综上所述,最终确定的累积效应修正模型为:
式中:P0为基准日负荷、ΔTe为平均实感温度变化量,且ΔTe=Te1-Te0、Te1为预测日平均实感温度、Te0为基准日平均实感温度、ΔT1=T0-T-1、ΔT2=T0-T-2、T0为基准日最高温度、T-1、T-2分别是累积日前一天和两天的最高温度。
而
其中,T0:基准日最高摄氏气温(℃);
T1:预测日最高摄氏气温(℃);
Rh:相对湿度(100%)
V:风速(m/s)。
Claims (1)
1.一种基于计及气象累积效应的电力系统最大日负荷预测系统的预测方法,其特征在于,所述的电力系统最大日负荷预测系统包括数据采集模块、数据筛分处理模块、负荷预测建模模块、负荷预测模块;
数据采集模块是将外部气象、负荷数据传递到系统中进行解析;数据筛分处理模块是将采集来的数据进行季节性筛选,并且将负荷分解为累积效应显著负荷和累积效应不显著负荷两大类分别存储;负荷预测建模模块,包括非累积日负荷基础预测模型建模模块及累积日负荷修正模型建模模块,非累积日负荷基础预测模型建模模块是提取累积效应不显著的负荷进行正常工作日模型的建立和预测,累积日负荷修正模型建模模块是在非累积日负荷基础预测模型基础上针对累积日进行建模和预测修正,最后建立综合预测模型;负荷预测模块是将负荷预测建模模块结果与气象及负荷数据库相联接获取负荷预测值;
所述的预测方法包括如下步骤:
(1)采用相关性分析方法,分析历史数据中影响最大日负荷的各种气象因素,包括温度、湿度、风速、降雨量以及多因素混合气象指数,多因素混合气象指数包括温湿指数、寒湿指数、实感温度指数和舒适度指数,计算各气象因素以及各多因素混合气象指数与电网还原负荷的相关系数,找出对负荷的影响大且影响稳定的因素,并用回归分析法建立关系方程;利用方差分析方法,对各影响因子进行比较,相关度大并且方差小的气象指数即认为是关键的影响因子;将不同的影响因子分别用在正常负荷建模和累积效应修正建模上,生成交叉修正方案;
(2)根据地级市多种气象指数和地级市负荷的相关性,基于地级市电力负荷和地级市气象指数计算省级多因素混合气象指数;计算方法为用各市的实际气象指数,结合各市的实际用电情况形成地区加权多因素混合气象因素;
(3)筛分气象累积效应显著与非显著日,提出一种气象指数与负荷关系的灵敏度实用化建模,建立不含温度累积效应时的正常负荷预测工程化模型;建模方法为:
将负荷分为经济负荷和气象负荷,其中气象负荷与温度存在关联性,通过历史样本拟合温度与负荷的灵敏度;经济负荷则通过历史样本,推算线性增长关系;两者叠加求得未来某时间段的负荷;
具体的建模步骤如下:
1)根据每年的历史样本,对夏季累积效应不显著统调负荷需求与最高温度进行一元二次回归分析,做出拟合曲线,求得每年的拟合函数及拟合度;
2)根据拟合得到的拟合函数f(x),对f(x)求一阶导数可得各年统调负荷需求对最高温度的灵敏度;将N个温度区间分别定义为[Ti1-Ti2],i=1,2...N,则第i个温度区间的灵敏度计算公式为:
Sj=f(Ti1-Ti2)
3)计算多年内相同温度区间的灵敏度平均值,取每个温度区间的灵敏度平均值作为新一年夏季统调负荷需求对温度的灵敏度,再对拟合温度与该灵敏度平均值进行一元一次回归分析,做出拟合曲线,求出各温度对应的灵敏度关系式;
4)当温度变化ΔT(℃)时,负荷的变化量为:
其中,ΔT=T0-Tr,T0表示预测日的温度,Tr表示参考日的温度;ΔLT为在正常气温条件下得到的负荷增量;
5)随着经济的发展,每年夏季的经济负荷均有一定的增长量,此时算出每月的负荷平均增长量ΔLtime;
6)同时考虑负荷随时间的变化量和负荷随温度的变化量,最终正常负荷的预测公式为:
L′=Lr+ΔLT+ΔLtime
其中,Lr为参考日的地区最大日负荷;
(4)在正常负荷预测工程化模型上,提出不同气象因素预测误差交叉修正方法,利用优化理论求解最优修正因子,建立考虑了气象累积效应的最大日负荷预测工程化模型;建模方案为:设定正常负荷预测模型对正常负荷预测是精准的,那么当累积效应显著时,用该模型计算出的预测负荷与实际负荷的偏差,认为主要是由累积效应引起的;因此,对累积效应进行修正负荷建模,实际上是对以上负荷偏差进行建模,则实际负荷即正常负荷与由累积效应引起的负荷偏差的叠加;
建立夏季温度累积效应修正负荷模型的具体步骤如下:
计算所有温度差值:
ΔT1=T0-T-1
及
ΔT2=T0-T-2
其中,T0是气温突变日的气温,T-1、T-2分别是突变日前一天和两天的气温,ΔT1、TΔ2作为累积效应修正公式的输入变量,反映了温度变化与人们感知惯性造成的偏差,与累积效应的物理现象相符合;
计算所有突变日的负荷偏差值ΔL:
ΔL=L-L′
其中,L为突变日的实际负荷值利用残差最小化的思想,通过最小二乘法进行多元回归分析,得到累积效应修正量ΔL′;
二元线性回归的表达式为:
ΔL′=f(ΔT1,ΔT2)=k1ΔT1+k2ΔT2+k3
二元二次回归的表达式为:
ΔL′=f1(ΔT1,ΔT2)=k1ΔT1 2+k2ΔT2 2+k3ΔT1ΔT2+k4ΔT1+k5ΔT2+k6
其中,k1,k2,k3,k4,k5,k6为待求的优化系数;综合考虑正常负荷模型与累积效应修正负荷模型,得到最终的实际负荷模型为:
L″=Lr+ΔLT+ΔLtime+f(ΔT1,ΔT2)=L″=L′+ΔL′
建立优化模型判断得出:
其中,n为历史累积日样本个数,ΔLi为第i个累积日的负荷偏差值,ΔL′i为第i个累积日计及累积效应的修正量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410367888.2A CN104156783B (zh) | 2014-07-29 | 2014-07-29 | 计及气象累积效应的电力系统最大日负荷预测系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410367888.2A CN104156783B (zh) | 2014-07-29 | 2014-07-29 | 计及气象累积效应的电力系统最大日负荷预测系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104156783A CN104156783A (zh) | 2014-11-19 |
CN104156783B true CN104156783B (zh) | 2017-07-25 |
Family
ID=51882279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410367888.2A Expired - Fee Related CN104156783B (zh) | 2014-07-29 | 2014-07-29 | 计及气象累积效应的电力系统最大日负荷预测系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104156783B (zh) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104376384B (zh) * | 2014-11-27 | 2017-07-25 | 广西大学 | 一种基于电力大数据分析的台风日最大日负荷预测系统 |
CN105160588A (zh) * | 2015-05-27 | 2015-12-16 | 南京国云电力有限公司 | 一种用电负荷特征分析的方法 |
CN105069525B (zh) * | 2015-07-30 | 2019-06-04 | 广西大学 | 全天候96点日负荷曲线预测及优化修正系统 |
CN105787259A (zh) * | 2016-02-17 | 2016-07-20 | 国网甘肃省电力公司武威供电公司 | 一种多元气象因素对负荷变化影响相关性的分析方法 |
CN106485068B (zh) * | 2016-10-08 | 2018-12-28 | 东南大学 | 一种空调负荷的聚合建模方法 |
CN106849057B (zh) * | 2016-12-30 | 2019-06-07 | 西安理工大学 | 基于现代内点法和灵敏度分析法的分布式风电源优化方法 |
CN107017624A (zh) * | 2017-04-13 | 2017-08-04 | 云南电网有限责任公司电力科学研究院 | 电力负荷周周期特性建模方法 |
CN107767037B (zh) * | 2017-09-30 | 2021-04-02 | 广东电网有限责任公司信息中心 | 一种用户用电构成解析方法 |
CN108229742B (zh) * | 2018-01-04 | 2021-10-22 | 国网浙江省电力公司电力科学研究院 | 一种基于气象数据以及数据趋势的负荷预测方法 |
CN108054779B (zh) * | 2018-01-04 | 2020-05-22 | 阳光电源股份有限公司 | 一种光储充换电站及其云储能系统优化控制方法 |
CN108053082B (zh) * | 2018-01-05 | 2021-06-01 | 国网湖南省电力有限公司 | 基于温度区间分解的电网中长期负荷预测方法 |
CN108932557A (zh) * | 2018-04-28 | 2018-12-04 | 云南电网有限责任公司临沧供电局 | 一种基于气温累积效应和灰色关联度的短期负荷预测模型 |
CN109886455B (zh) * | 2019-01-14 | 2020-12-08 | 国网江苏省电力有限公司苏州供电分公司 | 雷电气候下超短期负荷预测方法 |
CN109978261A (zh) * | 2019-03-27 | 2019-07-05 | 新奥数能科技有限公司 | 确定负荷预测模型的方法、装置、可读介质及电子设备 |
CN110445136A (zh) * | 2019-08-12 | 2019-11-12 | 安徽立卓智能电网科技有限公司 | 一种基于功率预测的新能源场站无功电压控制方法和平台 |
CN111210060B (zh) * | 2019-12-30 | 2023-04-28 | 国网宁夏电力有限公司信息通信公司 | 一种工作日期间机房温度预测方法 |
CN111191936B (zh) * | 2019-12-31 | 2020-09-15 | 深圳市国家气候观象台 | 一种台风风雨综合影响指数计算方法及存储装置 |
CN112113316B (zh) * | 2020-09-18 | 2022-03-11 | 国网辽宁省电力有限公司电力科学研究院 | 一种空调负荷提取的方法 |
CN112149309B (zh) * | 2020-09-30 | 2023-07-21 | 贵州电网有限责任公司 | 一种考虑温度-负荷相关性的调温负荷曲线拟合方法 |
CN113819581B (zh) * | 2021-08-23 | 2022-11-18 | 广西电网有限责任公司电力科学研究院 | 一种夏季空调负荷容量及其气温灵敏性分析方法 |
CN113705932B (zh) * | 2021-09-22 | 2023-02-14 | 华中科技大学 | 一种短期负荷预测方法及装置 |
CN114677019B (zh) * | 2022-03-30 | 2024-06-28 | 山东大学 | 基于随机源矩阵和层次分析的用电特征辨识方法及系统 |
CN116090620B (zh) * | 2022-12-30 | 2023-10-13 | 厦门大学 | 一种基于气温惯性效应的电量预测方法、终端设备及介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930356A (zh) * | 2012-11-16 | 2013-02-13 | 广东电网公司电力调度控制中心 | 基于气象因素敏感度的短期负荷预测方法 |
CN103268524A (zh) * | 2013-06-03 | 2013-08-28 | 国家电网公司 | 一种提高电网短期负荷预测准确度的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100830095B1 (ko) * | 2007-11-12 | 2008-05-20 | 충남대학교산학협력단 | 냉방부하 예측방법 |
-
2014
- 2014-07-29 CN CN201410367888.2A patent/CN104156783B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930356A (zh) * | 2012-11-16 | 2013-02-13 | 广东电网公司电力调度控制中心 | 基于气象因素敏感度的短期负荷预测方法 |
CN103268524A (zh) * | 2013-06-03 | 2013-08-28 | 国家电网公司 | 一种提高电网短期负荷预测准确度的方法 |
Non-Patent Citations (1)
Title |
---|
基于气象累积和ACA-GRNN短期电力负荷预测;李如琦等;《继电器》;20080216;第36卷(第4期);第58-61页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104156783A (zh) | 2014-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104156783B (zh) | 计及气象累积效应的电力系统最大日负荷预测系统及方法 | |
CN104156786B (zh) | 一种考虑气象多因素影响的非工作日最大日负荷预测系统 | |
Chen et al. | Evaluation of economic transformation and upgrading of resource-based cities in Shaanxi province based on an improved TOPSIS method | |
CN105260803B (zh) | 一种系统用电量预测方法 | |
CN105069525B (zh) | 全天候96点日负荷曲线预测及优化修正系统 | |
CN106712061B (zh) | 一种基于电动汽车可调度能力的日内优先调度方法 | |
CN105096004B (zh) | 一种水库群供调水系统实时调度方法 | |
CN108491969A (zh) | 基于大数据的空间负荷预测模型构建方法 | |
CN109858728B (zh) | 基于分行业用电特性分析的负荷预测方法 | |
CN103854068A (zh) | 一种居民小区短期负荷预测方法 | |
CN104376384B (zh) | 一种基于电力大数据分析的台风日最大日负荷预测系统 | |
CN105046953A (zh) | 一种短时交通流组合预测方法 | |
CN104008432A (zh) | 基于emd-kelm-ekf的微网短期负荷预测方法 | |
CN107122856A (zh) | 新态势下的空间饱和负荷预测方法 | |
CN101290713A (zh) | 一种结合周相似特性的分形交通流预测方法 | |
CN104751249A (zh) | 一种空间负荷预测方法 | |
CN105846425A (zh) | 一种基于通用风电预测误差模型的经济调度方法 | |
CN104657791A (zh) | 一种基于相关性分析的风电场群风速分布预测方法 | |
CN106548272A (zh) | 一种电动汽车快充设施综合性能的评价方法 | |
CN103985000A (zh) | 基于函数型非参数回归的中长期典型日负荷曲线预测方法 | |
CN103218668A (zh) | 一种基于地理加权泊松回归的县级交通事故的预测方法 | |
CN110852496A (zh) | 一种基于lstm循环神经网络的天然气负荷预测方法 | |
CN103530700A (zh) | 城区配电网饱和负荷综合预测方法 | |
CN106779139A (zh) | 基于小波分解和二阶灰色神经网络的短期风速预测方法 | |
CN105954821B (zh) | 一种用于数值大气模式的典型降雨事件选取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170725 |
|
CF01 | Termination of patent right due to non-payment of annual fee |