CN105140967A - 一种含新能源电力系统调峰需求的评估方法 - Google Patents

一种含新能源电力系统调峰需求的评估方法 Download PDF

Info

Publication number
CN105140967A
CN105140967A CN201510672598.3A CN201510672598A CN105140967A CN 105140967 A CN105140967 A CN 105140967A CN 201510672598 A CN201510672598 A CN 201510672598A CN 105140967 A CN105140967 A CN 105140967A
Authority
CN
China
Prior art keywords
net load
predicated error
sample
scene
period
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.)
Granted
Application number
CN201510672598.3A
Other languages
English (en)
Other versions
CN105140967B (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.)
Huazhong University of Science and Technology
State Grid Zhejiang Electric Power Co Ltd
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
Huazhong University of Science and Technology
State Grid Zhejiang Electric Power Co Ltd
China Electric Power Research Institute Co Ltd CEPRI
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 Huazhong University of Science and Technology, State Grid Zhejiang Electric Power Co Ltd, China Electric Power Research Institute Co Ltd CEPRI filed Critical Huazhong University of Science and Technology
Priority to CN201510672598.3A priority Critical patent/CN105140967B/zh
Publication of CN105140967A publication Critical patent/CN105140967A/zh
Application granted granted Critical
Publication of CN105140967B publication Critical patent/CN105140967B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers
    • 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Abstract

本发明公开了一种含新能源电力系统调峰需求的评估方法,包括:统计发生概率大于0.95的正净负荷预测误差边界值,发生概率大于0.95的负净负荷预测误差边界值;生成位于正净负荷预测误差边界和负净负荷预测误差边界之间的净负荷预测误差代表场景;根据所得的净负荷预测误差代表场景,对次日可能发生的净负荷场景进行预测;计算预测得到的净负荷场景中的最大峰谷差,作为所求的电力系统次日的调峰需求。本发明可以根据日前预测的新能源出力的随机特征而制定出更加准确的调峰需求容量,从而使得电力系统在保证安全运行的前提下更加经济。

Description

一种含新能源电力系统调峰需求的评估方法
技术领域
本发明属于含有风电、光伏发电等新能源电力系统运行技术领域,更具体地,涉及一种考虑新能源随机特性的电力系统调峰需求的评估方法。
背景技术
风电出力、光伏发电出力随机性强,预测精度低,运行方式复杂多变,其大规模并网后,将使电力系统调度运行面临巨大压力。调峰压力是大规模风电、光伏发电接入后电力系统调度运行面临的主要压力之一。准确地计算、确定含大规模新能源电力系统的调峰需求,是缓解电力系统调峰压力的有效手段之一。
当前,电力系统的调峰需求的确定主要有2种方法。第一种方法,从多年的电力系统监测到的负荷曲线中,选出系统运行的最小负荷率和日最大负荷,然后采用(日最大负荷)×(1-最小负荷率)的方法计算系统的调峰需求。第二种方法,将电力系统过去一年或多年的每个时刻的调峰需求从大到小进行排序,形成调峰需求持续曲线。然后,将调峰需求持续曲线进行分级,得到一个离散的调峰需求概率分布。最后,通过调峰需求的概率分布,可以得到不同置信度水平下,系统的调峰需求。第一种方法简单直接,但是由于没有考虑负荷和新能源电力的随机特性,难以适用于含大规模新能源的电力系统,且直接选取最小负荷率和日最大负荷的方法,往往使得所确定的电力系统调峰需求过于保守,降低电力系统运行的经济性。第二种方法在一定程度上可以考虑负荷和新能源的随机特性,但凭借经验划分调峰需求水平所得到概率分布,误差较大。此外,其在计算过程中,仅考虑过去的历史数据,没有利用或借助于最新的负荷和新能源电力的预测信息等,一定程度上也降低了所确定调峰需求的准确性。
综上,上述两种传统的调峰需求的确定方法,其评估的结果往往过于保守,不够准确,可用于在电力系统规划和设计阶段大致评估系统所需要的调峰容量。但是,并不适用于准确评估电力系统运行时的调峰需求,特别是难以适用于含有大规模新能源的电力系统中。
发明内容
针对现有技术的缺陷,本发明的目的在于提供一种含新能源电力系统调峰需求的评估方法,旨在解决现有技术难以跟踪新能源出力的随机变化且制定的调峰需求过于保守,导致电力系统运行不经济的问题。
为实现上述目的,本发明提供了一种含新能源电力系统调峰需求的评估方法,包括下述步骤:
S1:获取样本数据;
所述样本数据包括日负荷测量曲线、日风力发电测量曲线、日光伏发电测量曲线和日负荷预测曲线、日风力发电预测曲线、日光伏发电预测曲线;
S2:对所述样本数据进行处理,获得由负荷曲线、风电出力曲线、光伏发电出力曲线构成的净负荷曲线次日可能发生的预测净负荷曲线场景;
S3:根据所述预测净负荷曲线场景获得次日调峰需求评估值 P r = m a x j = 1 , ... n 0 { m a x t = 1 , ... , 96 ( P ‾ j , t A ) - m i n t = 1 , ... , 96 ( P ‾ j , t A ) } .
更进一步地,步骤S1中,所述样本数据中的曲线是指由每天间隔15min共96个时段的数据点组成的日曲线。
更进一步地,步骤S2具体为:
S21:根据公式 P i , t A = P i , t L , A - P i , t W , A - P i , t S , A P i , t F = P i , t L , F - P i , t W , F - P i , t S , F 分别将所述样本数据转换为实际净负荷样本和预测净负荷样本;
S22:根据所述实际净负荷样本和预测净负荷样本获得净负荷的预测误差样本;
S23:对净负荷预测误差样本的每个时刻的净负荷预测误差进行统计,获得发生概率大于0.95的正净负荷预测误差边界值和负净负荷预测误差边界值;
S24:根据所述正净负荷预测误差边界值和负净负荷预测误差边界值获得位于正净负荷预测误差边界与负净负荷预测误差边界之间的净负荷预测误差代表场景;
S25:根据n0个所述净负荷预测误差代表场景获得次日可能发生的预测净负荷曲线场景;
其中,t为时段的序号,t=1,2,…,96;i为样本的序号,i=1,2,…,N,N为样本总数,N=180~1825天;为第i个样本的第t个时段的实际的净负荷,为第i个样本的第t个时段的预测的净负荷;为第i个样本的第t个时段的实际的负荷,为第i个样本的第t个时段的预测的负荷;为第i个样本的第t个时段的实际的风电出力,为第i个样本的第t个时段的预测的风电出力;为第i个样本的第t个时段的实际的光伏发电出力,为第i个样本的第t个时段的预测的光伏发电出力;
所述净负荷的预测误差样本包括N个净负荷预测误差Ri,t为第i个净负荷预测误差样本的第t个时段的净负荷预测误差。
更进一步地,步骤S23具体为:
(1)将N个净负荷预测误差样本中每个时刻的净负荷预测误差分为大于0的正净负荷预测误差组和小于0的负净负荷预测误差组;
(2)将所述正净负荷预测误差组按照从小到大的顺序进行排序,并将负净负荷预测误差组中的数据按照从大到小的顺序进行排序;
(3)将排序之后的正净负荷预测误差组中的位于第个数据作为第t个时段的正净负荷预测误差的边界值,将排序之后的负净负荷预测误差组中的位于第个数据作为第t个时段的负净负荷预测误差的边界值。
其中,分别表示第t个时段的正和负净负荷预测误差组中所包含的数据的个数;表示取的整数。
更进一步地,步骤S24具体为:
(1)采用k均值聚类的方法将N个净负荷预测误差样本聚成n0类,获得n0个净负荷预测误差的代表场景;
其中,每类的概率等于该类所包含的样本数目除以N,聚类的中心作为该类的代表场景,n0为给定聚类数目,1<n0<N;
(2)将所述净负荷预测误差代表场景的每个时段的值分别与对应时刻的正净负荷预测误差边界值和负净负荷预测误差边界值进行大小比较,将第t时段中所有大于的值置为所有小于的值置为
(3)分别计算N个净负荷预测误差样本和n0个净负荷预测误差代表场景的每个时刻的均值和标准差;
(4)获得净负荷预测误差样本与净负荷预测误差代表场景之间的偏差 &epsiv; = m a x 1 &le; t &le; 96 { m a x { | E t 1 , N - E t 1 , n 0 | / E t 1 , N , | E t 2 , N - E t 2 , n 0 | / E t 2 , N } } ;
(5)判断ε≤ε,若是,则所生成的n0条净负荷预测误差曲线即为所求;若否,则转入步骤(1)中重新进行聚类;
其中,ε为最大偏差的阈值,0≤ε≤1,ε为偏差;分别为N个净负荷预测误差样本和n0个净负荷预测误差代表场景中第t时段的均值; 分别为N个净负荷预测误差样本和n0个净负荷预测误差代表场景中第t时段的标准差。
更进一步地,步骤S25具体为:
(1)获得次日的负荷预测曲线、风电出力预测曲线、光伏发电出力预测曲线,并根据公式计算次日预测净负荷曲线;
(2)根据所得的次日预测净负荷曲线和净负荷预测误差代表场景,利用公式计算得到次日可能发生的预测净负荷曲线场景;
其中,为次日预测净负荷曲线的第t个时段的预测净负荷;为次日负荷预测曲线中第t个时段的负荷;为次日风电出力预测曲线中第t个时段的风电出力;为次日光伏发电出力预测曲线中第t个时段的光伏发电出力;为次日可能发生的第j个预测净负荷曲线场景的第t个时段的实际净负荷;为第j个净负荷预测误差代表场景的第t个时段的净负荷预测误差,n0为预测净负荷曲线场景和预测净负荷预测误差代表场景的数目。
更进一步地,步骤S3中,所述次日调峰需求评估值其中,t为时段的序号,t=1,2…L,96;j为净负荷预测误差代表场景的序号,1≤j≤n0,n0为预测净负荷曲线代表场景数目,1<n0<N,N为样本总数,N=180~1825天;为次日可能发生的第j个预测净负荷曲线场景的第t个时段的实际净负荷;Pr为所求的次日调峰需求评估值。
本发明在制定调峰需求的过程中考虑了风电出力、光伏发电出力的随机性,考虑了负荷预测误差、风电出力预测误差以及光伏发电出力预测误差的随机性,且该调峰需求的确定基于日前预测的负荷曲线、风电出力曲线、光伏发电曲线而制定。因此,与传统的确定电力系统的调峰需求方法相比,本发明的方法可以根据日前(短期)预测的新能源出力的随机特征而制定出更加准确的调峰需求容量,从而使得电力系统在保证安全运行的前提下更加经济。
附图说明
图1是本发明实施例提供的一种含新能源电力系统调峰需求的评估方法的实现流程图;
图2本发明实施例提供的净负荷预测误差曲线样本示意图;
图3本发明实施例提供的净负荷预测误差发生概率大于0.95的正、负边界值示意图;
图4本发明实施例提供的净负荷预测误差代表场景(20个)示意图;
图5本发明实施例提供的次日可能发生的预测净负荷曲线场景(20个)示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
传统的调峰需求确定的方法往往难以跟踪新能源出力的随机变化、制定的调峰需求过于保守、导致电力系统运行不经济。本发明的目的在于克服传统调峰需求制定方法的不足,提供了一种可以考虑新能源出力随机变化特征的电力系统调峰需求的确定方法。所提的方法可以根据日前预测得到的新能源出力曲线,动态地制定电力系统的调峰需求,对新能源出力的变化就有很好的适应性,从而使得电力系统安全、经济运行。
本发明提出了一种适用于满足含有大规模新能源电力系统安全运行所需调峰容量的评估方法。该方法可以考虑负荷、风电、光伏发电的不确定性以及预测的不准确性,且在计算过程中利用了日前预测负荷、预测风电出力、预测光伏发电出力等最新信息,可以更好地跟踪和反映新能源出力的随机变化,从而更准确地计算出电力系统安全运行所需的调峰容量,使得电力系统安全、经济运行。
如图1所示,本发明提供的一种含新能源电力系统调峰需求的评估方法,包括以下步骤:
S1:从电力系统的数据采集与监视控制系统(SupervisoryControlAndDataAcquisition,以下简称SCADA系统)中,获取计算所需的数据。
计算所需的数据包括:电力系统半年以上的实际监测到的日负荷曲线、日风力发电曲线、日光伏发电曲线,和SCADA系统预测的日负荷曲线、日风力发电曲线、日光伏发电曲线。
其中,日负荷曲线、日风力发电曲线、日光伏发电曲线是指由每天间隔15min共96个时段的数据点组成的日曲线。
S2:预测由负荷曲线、风电出力曲线、光伏发电出力曲线构成的净负荷曲线次日可能发生的场景。
预测的步骤如下:
S21:将第1步得到的数据,按照式(1)转换成实际净负荷样本和预测净负荷样本。
P i , t A = P i , t L , A - P i , t W , A - P i , t S , A ...... ( 1. a )
P i , t F = P i , t L , F - P i , t W , F - P i , t S , F ...... ( 1. b )
式中,t表示时段的序号,t=1,2,…,96;i表示样本的序号,i=1,2,…,N,N为样本总数,一般取0.5~5年的样本,即N=180~1825天;分别表示第i个样本的第t个时段的实际的、预测的净负荷;分别表示第i个样本的第t个时段的实际的、预测的负荷;分别表示第i个样本的第t个时段的实际的、预测的风电出力;分别表示第i个样本的第t个时段的实际的、预测的光伏发电出力。
S22:根据上述S21得到的N个实际净负荷样本和预测净负荷样本,采用式(2)计算净负荷的预测误差,则可得到N个净负荷预测误差样本。
R i , t = ( P i , t A - P i , t F ) / P i , t A , i = 1 , 2 , ... , N , t = 1 , 2 , ... , 96 ... ... ( 2 )
式中,Ri,t为第i个净负荷预测误差样本的第t个时段的净负荷预测误差。
S23:对S22得到N个净负荷预测误差样本的每个时刻的净负荷预测误差进行统计,获得发生概率大于0.95的正净负荷预测误差边界值和负净负荷预测误差边界值,分别记为:其中,分别表示第t个时段的发生概率大于0.95的正和负净负荷预测误差边界值。
具体做法为:1)把N个净负荷预测误差样本中每个时刻的净负荷预测误差分为大于0和小于0的两组,分别称为正净负荷预测误差组和负净负荷预测误差组,假设每组数据长度分别为分别表示第t个时段的正和负净负荷预测误差组中所包含的数据的个数,t=1,2,…,96;2)将正净负荷预测误差组按照从小到大的顺序进行排序,将负净负荷预测误差组中的数据按照从大到小的顺序进行排序;3)分别取出排序之后的两组数据中位于第个和第个的数据作为第t个时段的正净负荷预测误差和负预测误差的边界值,分别记为其中,表示取的整数。
S24:生成位于正净负荷预测误差边界和负净负荷预测误差边界之间的净负荷预测误差代表场景。
具体做法为:1)给定聚类数目n0(1<n0<N),采用k均值聚类的方法将S23所得的N个净负荷预测误差样本聚成n0类,每类的概率等于该类所包含的样本数目除以N,聚类的中心作为该类的代表场景,则得到n0个净负荷预测误差的代表场景;2)将所得的净负荷预测误差代表场景的每个时段t=1,2,…,96的值分别与对应时刻的进行大小比较,将第t(t=1,2,…,96)时段中所有大于的值置为所有小于的值置为得到调整后的净负荷预测误差代表场景3)分别计算N个净负荷预测误差样本和n0个调整后的净负荷预测误差代表场景的每个时刻的均值和标准差,均值和标准差的计算方法如式(3);4)按照式(4)计算净负荷预测误差样本与净负荷预测误差代表场景两者之间的偏差;式(4)的含义是:将所有时段中的均值和标准差的最大偏差作为最终的偏差值;5)给定最大偏差的阈值ε,0≤ε≤1,若ε≤ε,则所生成的n0条净负荷预测误差曲线即为所求;否则,转1)重新进行聚类。
E t 1 , N = &Sigma; i = 1 N R i , t / N , E t 2 , N = &Sigma; i = 1 N ( R i , t ) 2 / N - ( E t 1 , N ) 2 ...... ( 3. a )
E t 1 , n 0 = &Sigma; i = 1 n 0 R &OverBar; j , t p j , E t 2 , n 0 = &Sigma; j = 1 n 0 ( R &OverBar; j , t ) 2 p j - ( E t 1 , n 0 ) 2 ...... ( 3. b )
&epsiv; = m a x 1 &le; t &le; 96 { m a x { | E t 1 , N - E t 1 , n 0 | / E t 1 , N , | E t 2 , N - E t 2 , n 0 | / E t 2 , N } } ...... ( 4 )
式中,ε为偏差;分别为N个净负荷预测误差样本和n0个净负荷预测误差代表场景中第t时段的均值;分别为N个净负荷预测误差样本和n0个净负荷预测误差代表场景中第t时段的标准差;为第j个净负荷预测误差代表场景的第t个时段的净负荷预测误差,pj为第j个的概率,j为净负荷预测误差代表场景的序号,1≤j≤n0
S25:根据S24所得的n0个净负荷预测误差代表场景,计算次日可能发生的预测净负荷曲线场景。
具体做法为:1)从SCADA系统中获取次日的负荷预测曲线、风电出力预测曲线、光伏发电出力预测曲线,采用式(1.b)计算次日预测净负荷曲线;根据步骤S24计算得到的净负荷预测误差代表场景,采用式(5)计算得到次日可能发生的预测净负荷曲线场景。
P &OverBar; j , t A = P &OverBar; t F / ( 1 - R &OverBar; j , t ) , j = 1 , 2 , ... n 0 ...... ( 5 ) ; 式中,为次日可能发生的第j个预测净负荷曲线场景的第t个时段的预测净负荷;为次日预测净负荷曲线的第t个时段的预测净负荷;步骤S24所得的第j个净负荷预测误差代表场景的第t个时段的净负荷预测误差,n0指的是预测净负荷曲线场景和预测净负荷预测误差代表场景的数目。
S3:根据第2步得到的次日可能发生的预测净负荷曲线场景,采用式(6)计算次日的调峰需求。式(6)的含义是:将预测净负荷曲线场景中的最大峰谷差,作为次日的调峰需求。其中峰谷差等于预测净负荷曲线场景中所有时段的最大值减去所有时段的最小值。
P r = m a x j = 1 , ... n 0 { m a x t = 1 , ... , 96 ( P &OverBar; j , t A ) - m i n t = 1 , ... , 96 ( P &OverBar; j , t A ) } ...... ( 6 ) ; 式中,Pr即为评估得到的次日调峰需求。
本发明的优点和积极效果是:在制定调峰需求的过程中考虑了风电出力、光伏发电出力的随机性,考虑了负荷预测误差、风电出力预测误差以及光伏发电出力预测误差的随机性,且该调峰需求的确定基于日前预测的负荷曲线、风电出力曲线、光伏发电曲线而制定。因此,与传统的确定电力系统的调峰需求方法相比,本发明的方法可以根据日前(短期)预测的新能源出力的随机特征而制定出更加准确的调峰需求容量,从而使得电力系统在保证安全运行的前提下更加经济。
以下结合附图对本发明实施案例作进一步详细说明。
实施步骤1:从SCADA系统中获取计算数据,包括从2015年1月1日到2015年6月30日共半年的实际监测日负荷曲线、日风力发电曲线、日光伏发电曲线;预测的日负荷曲线、日风力发电曲线、日光伏发电曲线。
实施步骤2:计算净负荷曲线样本和净负荷曲线样本的预测误差曲线。所得的净负荷误差曲线样本如附图2所示。
实施步骤3:计算发生概率大于0.95的正、负净负荷预测误差边界。预测得到的96个时段的净负荷预测误差边界如附图2所示。从图3可以看出,最大正净负荷预测误差略大于0.3,最小净负荷预测误差的绝对值略小于0.3。
实施步骤4:生成在正、负净负荷预测误差边界之间的净负荷预测误差样本。此实施例中,取聚类数目n0=20,均值和标准差最大偏差的阈值ε=0.3。所生成的净负荷预测误差代表场景如附图4所示。从图4可以看出,所生成的净负荷预测误差场景落在正、负净负荷预测场景的边界之内。
实施步骤5:从SCADA系统中获取次日(2015-7-3)的预测净负荷曲线,计算次日可能发生的预测净负荷曲线场景。生成的次日可能发生的预测净负荷曲线场景如附图5所示。从附图5可以看出,该日(2015-7-3)实际所发生的场景落在本发明方法所生成的场景之内,从而说明本发明方法可以较为准确预测出实际净负荷曲线可能发生的场景。
实施步骤6:取所有可能发生的预测净负荷曲线场景中的最大峰谷差作为系统的调峰需求。每个场景的峰谷差如表1所示。从表1中可见,第3个场景的峰谷差最大,等于2065.5MW。因此,系统的调峰需求为2065.5MW。在本实施例中,2015-7-3日实际的峰谷差为1537MW。可见,本发明方法在满足实际需求的同时,仍有一定冗余以应对新能源出力的随机变化。
表1代表的场景的峰谷差(MW)
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.一种含新能源电力系统调峰需求的评估方法,其特征在于,包括下述步骤:
S1:获取样本数据;
所述样本数据包括日负荷测量曲线、日风力发电测量曲线、日光伏发电测量曲线和日负荷预测曲线、日风力发电预测曲线、日光伏发电预测曲线;
S2:对所述样本数据进行处理,获得净负荷曲线次日可能发生的预测净负荷曲线场景;
S3:根据所述预测净负荷曲线场景获得次日调峰需求评估值 P r = m a x j = 1 , ... n 0 { m a x t = 1 , ... , 96 ( P &OverBar; j , t A ) - min t = 1 , ... , 96 ( P &OverBar; j , t A ) } .
2.如权利要求1所述的评估方法,其特征在于,步骤S1中,所述样本数据中的曲线是指由每天间隔15min共96个时段的数据点组成的日曲线。
3.如权利要求1或2所述的评估方法,其特征在于,步骤S2具体为:
S21:根据公式 P i , t A = P i , t L , A - P i , t W , A - P i , t S , A P i , t F = P i , t L , F - P i , t W , F - P i , t S , F 分别将所述样本数据转换为实际净负荷样本和预测净负荷样本;
S22:根据所述实际净负荷样本和预测净负荷样本获得净负荷的预测误差样本;
S23:对净负荷预测误差样本的每个时刻的净负荷预测误差进行统计,获得发生概率大于0.95的正净负荷预测误差边界值和负净负荷预测误差边界值;
S24:根据所述正净负荷预测误差边界值和负净负荷预测误差边界值获得位于正净负荷预测误差边界与负净负荷预测误差边界之间的净负荷预测误差代表场景;
S25:根据n0个所述净负荷预测误差代表场景获得次日可能发生的预测净负荷曲线场景;
其中,t为时段的序号,t=1,2,…,96;i为样本的序号,i=1,2,…,N,N为样本总数,N=180~1825天;为第i个样本的第t个时段的实际的净负荷,为第i个样本的第t个时段的预测的净负荷;为第i个样本的第t个时段的实际的负荷,为第i个样本的第t个时段的预测的负荷;为第i个样本的第t个时段的实际的风电出力,为第i个样本的第t个时段的预测的风电出力;为第i个样本的第t个时段的实际的光伏发电出力,为第i个样本的第t个时段的预测的光伏发电出力;
所述净负荷的预测误差样本包括N个净负荷预测误差Ri,t为第i个净负荷预测误差样本的第t个时段的净负荷预测误差。
4.如权利要求3所述的评估方法,其特征在于,步骤S23具体为:
(1)将N个净负荷预测误差样本中每个时刻的净负荷预测误差分为大于0的正净负荷预测误差组和小于0的负净负荷预测误差组;
(2)将所述正净负荷预测误差组按照从小到大的顺序进行排序,并将负净负荷预测误差组中的数据按照从大到小的顺序进行排序;
(3)将排序之后的正净负荷预测误差组中的位于第个数据作为第t个时段的正净负荷预测误差的边界值,将排序之后的负净负荷预测误差组中的位于第个数据作为第t个时段的负净负荷预测误差的边界值;
其中,分别表示第t个时段的正和负净负荷预测误差组中所包含的数据的个数;表示取的整数。
5.如权利要求3所述的评估方法,其特征在于,步骤S24具体为:
(1)采用k均值聚类的方法将N个净负荷预测误差样本聚成n0类,获得n0个净负荷预测误差的代表场景;
其中,每类的概率等于该类所包含的样本数目除以N,聚类的中心作为该类的代表场景,n0为给定聚类数目,1<n0<N;
(2)将所述净负荷预测误差代表场景的每个时段的值分别与对应时刻的正净负荷预测误差边界值和负净负荷预测误差边界值进行大小比较,将第t时段中所有大于的值置为所有小于的值置为
(3)分别计算N个净负荷预测误差样本和n0个净负荷预测误差代表场景的每个时刻的均值和标准差;
(4)获得净负荷预测误差样本与净负荷预测误差代表场景之间的偏差 &epsiv; = m a x 1 &le; t &le; 96 { m a x { | E t 1 , N - E t 1 , n 0 | / E t 1 , N , | E t 2 , N - E t 2 , n 0 | / E t 2 , N } } ;
(5)判断ε≤ε,若是,则所生成的n0条净负荷预测误差曲线即为所求;若否,则转入步骤(1)中重新进行聚类;
其中,ε为最大偏差的阈值,0≤ε≤1,ε为偏差;分别为N个净负荷预测误差样本和n0个净负荷预测误差代表场景中第t时段的均值;分别为N个净负荷预测误差样本和n0个净负荷预测误差代表场景中第t时段的标准差。
6.如权利要求3所述的评估方法,其特征在于,步骤S25具体为:
(1)获得次日的负荷预测曲线、风电出力预测曲线、光伏发电出力预测曲线,并根据公式计算次日预测净负荷曲线;
(2)根据所得的次日预测净负荷曲线和净负荷预测误差代表场景,利用公式计算得到次日可能发生的预测净负荷曲线场景;
其中,为次日预测净负荷曲线的第t个时段的预测净负荷;为次日负荷预测曲线中第t个时段的负荷;为次日风电出力预测曲线中第t个时段的风电出力;为次日光伏发电出力预测曲线中第t个时段的光伏发电出力;为次日可能发生的第j个预测净负荷曲线场景的第t个时段的实际净负荷;为第j个净负荷预测误差代表场景的第t个时段的净负荷预测误差,n0为预测净负荷曲线场景和预测净负荷预测误差代表场景的数目。
7.如权利要求3所述的评估方法,其特征在于,步骤S3中,所述次日调峰需求评估值其中,t为时段的序号,t=1,2,…,96;j为净负荷预测误差代表场景的序号,1≤j≤n0,n0为预测净负荷代表场景的数目,1<n0<N,N为样本总数,N=180~1825天;为次日可能发生的第j个预测净负荷曲线场景的第t个时段的实际净负荷;Pr为所求的次日调峰需求评估值。
CN201510672598.3A 2015-10-16 2015-10-16 一种含新能源电力系统调峰需求的评估方法 Expired - Fee Related CN105140967B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510672598.3A CN105140967B (zh) 2015-10-16 2015-10-16 一种含新能源电力系统调峰需求的评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510672598.3A CN105140967B (zh) 2015-10-16 2015-10-16 一种含新能源电力系统调峰需求的评估方法

Publications (2)

Publication Number Publication Date
CN105140967A true CN105140967A (zh) 2015-12-09
CN105140967B CN105140967B (zh) 2017-07-18

Family

ID=54726224

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510672598.3A Expired - Fee Related CN105140967B (zh) 2015-10-16 2015-10-16 一种含新能源电力系统调峰需求的评估方法

Country Status (1)

Country Link
CN (1) CN105140967B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108462165A (zh) * 2018-01-19 2018-08-28 长沙理工大学 一种新能源接入电力系统的负荷特性评估方法
CN109002912A (zh) * 2018-07-10 2018-12-14 国网湖南省电力有限公司电力科学研究院 一种水风光互补发电系统调峰能力评估方法
CN113139675A (zh) * 2021-03-08 2021-07-20 深圳职业技术学院 一种园区综合能源系统综合电负荷区间预测方法
CN113410846A (zh) * 2021-07-12 2021-09-17 广东电网有限责任公司 一种结合网络重构的配电网日前无功电压双层优化方法
CN113962598A (zh) * 2021-11-11 2022-01-21 国网天津市电力公司 新能源日运行调峰需求测算方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103259285A (zh) * 2013-05-03 2013-08-21 国家电网公司 含大规模风电电力系统的短期运行优化方法
CN103997039A (zh) * 2014-05-30 2014-08-20 国网甘肃省电力公司 基于概率区间预测的计及风电接纳的旋转备用区间预测方法
CN104268800A (zh) * 2014-09-30 2015-01-07 清华大学 基于场景库的风电并网调峰平衡判定方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103259285A (zh) * 2013-05-03 2013-08-21 国家电网公司 含大规模风电电力系统的短期运行优化方法
CN103997039A (zh) * 2014-05-30 2014-08-20 国网甘肃省电力公司 基于概率区间预测的计及风电接纳的旋转备用区间预测方法
CN104268800A (zh) * 2014-09-30 2015-01-07 清华大学 基于场景库的风电并网调峰平衡判定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
尤毅 等: "多时间尺度下基于主动配电网的分布式电源协调控制", 《电力系统自动化》 *
赵书强 等: "基于风电预测误差随机性的火储联合相关机会规划调度", 《中国电机工程学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108462165A (zh) * 2018-01-19 2018-08-28 长沙理工大学 一种新能源接入电力系统的负荷特性评估方法
CN108462165B (zh) * 2018-01-19 2021-04-20 长沙理工大学 一种新能源接入电力系统的负荷特性评估方法
CN109002912A (zh) * 2018-07-10 2018-12-14 国网湖南省电力有限公司电力科学研究院 一种水风光互补发电系统调峰能力评估方法
CN113139675A (zh) * 2021-03-08 2021-07-20 深圳职业技术学院 一种园区综合能源系统综合电负荷区间预测方法
CN113410846A (zh) * 2021-07-12 2021-09-17 广东电网有限责任公司 一种结合网络重构的配电网日前无功电压双层优化方法
CN113410846B (zh) * 2021-07-12 2022-03-15 广东电网有限责任公司 一种结合网络重构的配电网日前无功电压双层优化方法
CN113962598A (zh) * 2021-11-11 2022-01-21 国网天津市电力公司 新能源日运行调峰需求测算方法及装置
CN113962598B (zh) * 2021-11-11 2024-05-07 国网天津市电力公司 新能源日运行调峰需求测算方法及装置

Also Published As

Publication number Publication date
CN105140967B (zh) 2017-07-18

Similar Documents

Publication Publication Date Title
Shi et al. Hybrid forecasting model for very-short term wind power forecasting based on grey relational analysis and wind speed distribution features
CN103259285B (zh) 含大规模风电电力系统的短期运行优化方法
CN103545832B (zh) 一种基于发电预测误差的光伏系统储能容量配置方法
CN103258142B (zh) 基于天气模式的风电场短期风速预测方法
CN103475021B (zh) 一种基于统计模型的风电场弃风电量确定方法
CN105140967A (zh) 一种含新能源电力系统调峰需求的评估方法
CN102479347B (zh) 基于数据驱动的风电场短期风速预测方法和系统
CN107528350B (zh) 一种适应中长期电源规划的风电出力典型场景生成方法
CN102478584B (zh) 基于小波分析的风电场风速预测方法及系统
CN104882905A (zh) 一种考虑暂态安全约束的新能源接纳能力评估方法
CN104077632A (zh) 一种基于深度神经网络的风电场功率预测方法
CN103065202A (zh) 一种基于组合核函数的风电场超短期风速预测方法
CN103810535A (zh) 电力系统风电消纳能力评估方法
CN104463356A (zh) 一种基于多维信息人工神经网络算法的光伏发电功率预测方法
CN103473393A (zh) 一种考虑随机概率的输电裕度控制模型建模方法
CN104820868A (zh) 一种基于风电概率进行预测的调度展示方法
CN103345585A (zh) 基于支持向量机的风功率预测校正方法及系统
CN103326394B (zh) 一种计及风电波动性的多场景概率优化调度方法
CN102904248B (zh) 基于风电出力不确定性集合的电力系统调度方法
CN102749471A (zh) 一种短期风速、风电功率预测方法
CN105305485A (zh) 一种消纳大规模间歇性能源的安全约束经济调度方法
CN104239962A (zh) 一种基于相关性分析的区域电网总风功率短期预测方法
CN105956713A (zh) 一种新能源年/月电量计划制定方法
CN107330538A (zh) 一种气候变化条件下水库适应性调度规则编制的方法
CN106208136A (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170718

Termination date: 20211016