CN109301877B - 一种分布式电源及节点负荷典型运行场景集生成方法 - Google Patents
一种分布式电源及节点负荷典型运行场景集生成方法 Download PDFInfo
- Publication number
- CN109301877B CN109301877B CN201811067310.XA CN201811067310A CN109301877B CN 109301877 B CN109301877 B CN 109301877B CN 201811067310 A CN201811067310 A CN 201811067310A CN 109301877 B CN109301877 B CN 109301877B
- Authority
- CN
- China
- Prior art keywords
- scene
- state
- scenes
- probability
- time
- 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 38
- 230000009467 reduction Effects 0.000 claims abstract description 31
- 239000011159 matrix material Substances 0.000 claims abstract description 21
- 230000003068 static effect Effects 0.000 claims abstract description 21
- 238000005070 sampling Methods 0.000 claims abstract description 19
- 230000007704 transition Effects 0.000 claims abstract description 15
- 238000005315 distribution function Methods 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims abstract description 7
- 238000009826 distribution Methods 0.000 claims description 34
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 description 8
- 238000011161 development Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 230000002457 bidirectional effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/14—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B70/00—Technologies for an efficient end-user side electric power management and consumption
- Y02B70/30—Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
- Y02B70/3225—Demand response systems, e.g. load shedding, peak shaving
-
- Y—GENERAL 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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
- Y04S20/00—Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
- Y04S20/20—End-user application control systems
- Y04S20/222—Demand response systems, e.g. load shedding, peak shaving
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种分布式电源及节点负荷典型运行场景集生成方法,确定所要生成的分布式电源及节点负荷场景集的初始时刻为t0及该时刻的离散状态数为M0,设定最大能处理的场景序列数为Nmax,典型分布式电源及节点负荷场景数目为N0;根据分布式电源及节点负荷功率经验概率分布函数,采用等间距抽样生成M0个概率相等的t0时刻的初始静态场景;通过状态离散,初始静态场景离散为M0个状态并构造转移概率矩阵;由得到的M0个状态分布通过N1次MCMC抽样得到t1时刻的M1=M0*N1个场景序列;对包含M1个场景的原始场景集进行变速率削减;将不超过Nmax个长度为L的场景序列进行场景削减,生成最终N0个的典型场景集。本发明对分布式电源及负荷场景生成具有普适性。
Description
技术领域
本发明涉及主动配电网态势感知技术领域,特别是一种分布式电源及节点负荷典型运行场景集生成方法。
背景技术
随着能源短缺和环境污染问题日益紧迫,优化能源结构,实现清洁低碳发展,越来越受到世界各国的关注和重视。可再生能源,尤其是小容量、分散化、以就地消纳为主的分布式电源(Distributed Genenation,DG)将成为世界能源发展的主流。
对于处于系统末端、直接面向电力用户的配电网而言,分布式电源的高密度接入使配电网由“无源”向“有源”发展、潮流由“单向”向“双向”转变,对电网运行的可靠性和安全性以及电网调控能力都带来了极大挑战。若仍沿用传统配电网的调控运行模式,分布式电源密集接入地区的配电网将存在局部电压越限、线路功率越限的可能。而随着经济发展电力用户的需求由“稳定”向“优质”逐步提高,采用必要技术手段对有源配电网进行管控,以实现资源的有效利用和电网的安全经济运行是电网发展的必然要求。
主动配电网(Active Distribution Network,ADN)是一种通过使用灵活的网络拓扑结构对潮流进行管控,进而主动控制、主动管理分布式能源的配电系统。它是消纳分布式电源的主要载体,同时也连接着大电网和微电网共同构筑坚强智能电网。与类似的“有源配电网”概念强调配电网接有分布式电源这一物理特性相比,主动配电网应具有更高的可观可控水平,能够积极应对区域内分布式电源的随机性与波动性特征。它的“主动性”体现在,通过控制中心可以观测到配网的负荷和分布式电源的实时状态,进而借助态势感知技术进行风险预判,并积极制定应对策略,通过控制中心对分布式电源、储能、无功补偿装置、可控负荷等灵活性资源的调控实现主动管控。它对电网运行状态进行主动控制、对可控资源实现主动管理的愿景是依托主动配电网态势感知技术实现的,为了使主动配电网能够先于配电网运行实时状态制定管控计划,主动配电网态势预测技术亟待研究发展,而主动配电网态势预测是在态势觉察、态势理解的基础上开展的,具体内容是针对配电网中各种不确定因素如负荷需求、分布式电源出力等进行预测,根据预测结果对未来场景进行安全分析评估并针对分析场景进行预警。
而现有的主动配电网典型运行场景生成方法主要都集中在分布式电源尤其是风电功率的典型场景模拟,且大部分涉及场景生成的方法都是生成典型的静态场景集服务于后续的电网优化调度工作,对于连续时序动态场景生成的研究仍然十分匮乏;对于场景削减而言,现有的绝大部分场景削减算法都与场景规模呈现平方倍的时间复杂度,效率较为低下。
发明内容
本发明所要解决的技术问题是克服现有技术的不足而提供一种分布式电源及节点负荷典型运行场景集生成方法,针对配电网中分布式电源出力以及负荷需求的不确定性进行预测,采用场景分析法,通过抽样将连续概率分布的随机量离散为包含一系列不同出现概率的场景集合,实现随机问题到确定性问题的转化,在场景生成阶段整合基于历史功率数据经验概率分布的静态生成方法和基于马尔可夫链-蒙特卡洛方法的动态生成方法实现连续时序动态场景生成,在场景削减阶段采用基于概率距离的变速率动态削减的场景削减算法,降低算法时间复杂度,通过两个阶段的结合实现从数据到数据的分布式电源及节点负荷典型运行场景集生成,对主动配电网中各种类型的分布式电源及负荷场景生成具有普适性。
本发明为解决上述技术问题采用以下技术方案:
根据本发明提出的一种分布式电源及节点负荷典型运行场景集生成方法,包括以下步骤:
步骤10、确定所要生成的分布式电源及节点负荷场景集的初始时刻为t0及该时刻的离散状态数为M0,设定最大能处理的场景序列数为Nmax,典型分布式电源及节点负荷场景数目为N0,预设目标时刻为tL;
步骤40、对由步骤30得到的M0个状态分布通过N1次马尔可夫链-蒙特卡洛方法MCMC抽样得到t1时刻的M1=M0*N1个场景序列,每个场景序列对应的发生概率为被抽样场景概率除以抽样次数;
步骤50、对包含M1个场景的原始场景集进行变速率削减:比较M1与设定最大能处理的场景序列数Nmax的大小,若M1大于Nmax则进行场景削减,使场景数削减至Nmax,否则返回步骤40,更新状态分布数M0为M1,模拟时间轴前进至t2时刻继续动态场景生成,直至模拟时间轴得到预设时刻tL;
步骤60、将不超过Nmax个长度为L的场景序列进行场景削减,生成最终N0个的典型场景集。
作为本发明所述的一种分布式电源及节点负荷典型运行场景集生成方法进一步优化方案,所述步骤10中,对于所确定的离散状态数M0,根据初始时刻t0的对应的月份和时段的历史功率数据来设置,其大小考虑转移概率矩阵P的阶数以及所生成模拟序列对原始数据的拟合度。
其中,F(*)为根据格利文科定理由历史功率数据获得的经验分布函数。
作为本发明所述的一种分布式电源及节点负荷典型运行场景集生成方法进一步优化方案,步骤30中的状态离散过程为:将步骤20中所选历史功率数据中功率的上下限分别设为Pmax和Pmin,则出力区间为[Pmin,Pmax],根据精度需求将出力区间划分为M0个离散状态,每个状态覆盖的功率区间宽度为(Pmax-Pmin)/M0,这M0个状态分别记作X1,X2,…,XM0,Xn为第n个离散状态,n=1,2,...,M0,将初始静态场景分别与这M0个状态对应即得到离散状态
作为本发明所述的一种分布式电源及节点负荷典型运行场景集生成方法进一步优化方案,步骤30中的转移概率矩阵如下式所示:
P{Xt+1=Xj|Xt=Xi}=pij i,j∈[1,2,...,M0]
其中,Xt+1为t+1时刻的状态,Xt为t时刻的状态,Xj为M0个状态中第j个状态,Xi为M0个状态中第i个状态,P{Xt+1=Xj|Xt=Xi}表示从t时刻的状态Xt转移到t+1时刻的状态Xt+1的概率,pij为从状态Xi转移到状态Xj的概率,nij为所选时段的历史功率数据经过离散化后每两个相邻时刻由状态Xi转移到状态Xj的总次数。
作为本发明所述的一种分布式电源及节点负荷典型运行场景集生成方法进一步优化方案,步骤40中一次马尔可夫链-蒙特卡洛方法抽样具体步骤如下:
①从M0个状态中取一个状态Xx,x=1,2,...,M0;
②由标准均匀分布产生随机数u;
③根据转移概率矩阵P计算累计概率分布矩阵Pcum,矩阵Pcum中的第x行第y列的元素为Pcum,xy,其中x,y=1,2,...,M0;
④将随机数u与Pcum的第x行进行比较,若Pcum,x(y-1)<u≤Pcum,xy,则选取t1时刻为状态Xy;同时状态Xy对应的光伏输出功率Sy=Sy,min+u(Sy,max-Sy,min),Sy.min为状态Xy对应区间的最小功率,Sy.max为状态Xy对应区间的最大功率,随机数u取值范围在0到1之间。
作为本发明所述的一种分布式电源及节点负荷典型运行场景集生成方法进一步优化方案,步骤40中,
其中,Pxk为转移概率矩阵P中第x行第k列的元素。
作为本发明所述的一种分布式电源及节点负荷典型运行场景集生成方法进一步优化方案,步骤50中,对包含M1个场景的原始场景集进行变速率削减的具体步骤如下:
5.1、计算场景间二范数d(cv,cw):
5.2、确定削减速率V
由各场景最近场景组成集合D,按照概率距离从小到大排序,当前20%场景的平均概率距离小于后20%场景的平均概率距离10%时,采用大倍速削减,即V=0.01*M1;当前20%场景的平均概率距离大于后20%场景的平均概率距离10%且小于15%时,采用小倍速削减,即V=0.001*M1;当前20%场景的平均概率距离大于后20%场景的平均概率距离的15%时,采用逐个削减,即V=1;
其中,p(v)为M1个场景中第v个场景发生的概率,d(cv,cw)为场景cv与场景cw的二范数,cv、cw分别为M1个场景中第v、w个场景;
5.3、根据选定速率削减V个满足下式的场景cs*:
其中p(s)分别为场景cs*与场景cs发生的概率,cs*为当前M1个场景中被削减的第s*个场景,cs为当前M1个场景中除了场景cs*外的第s个场景,p(w)为场景cw发生的概率,cw为当前M1个场景中的第w个场景,p(z)为当前M1个场景中第z个场景发生的概率;
5.5、变更场景总数M1=M1-V*1,返回步骤5.2重复削减过程,直至削减后的场景数等于Nmax。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
本发明的分布式电源及节点负荷典型运行场景集生成方法,整合了基于历史功率数据经验概率分布的静态生成方法和基于马尔可夫链-蒙特卡洛方法的动态生成方法:采用完全基于历史功率数据的经验概率分布函数,对主动配电网中各种类型的分布式电源及负荷场景具有可推广性和实用性,采用等间距抽样来生成初始场景集,对抽样样本的数量要求不高,应用马尔可夫链-蒙特卡洛方法动态生成运行场景,较好表征了连续时间断面输出功率的相关性;基于概率距离的变速率动态削减的场景削减方法,能够快速有效削减随着模拟时间轴前进而呈指数级增长的场景序列个数,得到能充分逼近整个样本空间的典型场景集。
附图说明
图1是本发明实施例的流程图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
如图1所示,本发明的分布式电源及节点负荷典型运行场景集生成策略,包括下述步骤:
步骤10、确定所要生成的分布式电源及节点负荷场景集的初始时刻为t0及该时刻的离散状态数为M0,设定最大能处理的场景序列数为Nmax,典型分布式电源及节点负荷场景数目为N0,预设目标时刻为tL。
步骤20、根据由历史功率数据获取的分布式电源及节点负荷功率经验概率分布函数,采用等间距抽样生成M0个概率相等的t0时刻的初始静态场景为t0时刻的第n个初始静态场景,n=1,2,…,M0;如式(1)-(3)所示:
Y=F(x) (1)
x=F-1(Y) (2)
式(1)中x为对于t0时刻的功率变量,F(x)为根据格利文科定理由历史功率数据获得的经验分布函数;式(2)为功率变量x的概率分布函数逆函数;式(3)中,设生成的初始场景数为M0,将概率分布函数分为M0个等概率空间,即将概率分布函数y轴等分为M0份,并对每个区间的中点yn进行抽样,将yn代入概率分布逆函数F-1得到的x*n即为第n个样本值。
P{Xt+1=Xj|Xt=Xi}=pij i,j∈[1,2,...,M0] (4)
本步骤状态离散过程为:将步骤20中所选历史功率数据中功率的上下限分别设为Pmax和Pmin,则出力区间为[Pmin,Pmax],根据精度需求将出力区间划分为M0个离散状态,每个状态覆盖的功率区间宽度为(Pmax-Pmin)/M0,这M0个状态分别记作X1,X2,…,XM0,Xn为第n个离散状态,n=1,2,...,M0,将初始静态场景分别与这M0个状态对应即得到离散状态式(4)表示若时刻t对应的状态为Xt,则下一时刻t+1对应的Xt+1转移到Xt的概率为pij,pij为所构造转移概率矩阵P=(pij)M0*M0中的元素,其计算式(5)中,nij为所选时段的历史功率数据经过离散化后每两个相邻时刻由状态Xi转移到状态Xj的总次数;Xt+1为t+1时刻的状态,Xt为t时刻的状态,Xj为M0个状态中第j个状态,Xi为M0个状态中第i个状态,P{Xt+1=Xj|Xt=Xi}表示从t时刻的状态Xt转移到t+1时刻的状态Xt+1的概率,pij为从状态Xi转移到状态Xj的概率,nij为所选时段的历史功率数据经过离散化后每两个相邻时刻由状态Xi转移到状态Xj的总次数。
步骤40、对由步骤30得到的M0个状态分布通过N1次马尔可夫链-蒙特卡洛方法(Markov Chain Monte Carlo,MCMC)抽样得到t1时刻的M1=M0*N1个场景序列,每个场景序列对应的发生概率为被抽样场景概率除以抽样次数;本步骤中一次马尔可夫链-蒙特卡洛方法抽样具体步骤如下:
步骤40中一次马尔可夫链-蒙特卡洛方法抽样具体步骤如下:
①从M0个状态中取一个状态Xx,x=1,2,...,M0;
②由标准均匀分布产生随机数u;
③根据转移概率矩阵P计算累计概率分布矩阵Pcum,矩阵Pcum中的第x行第y列的元素为Pcum,xy,其中x,y=1,2,...,M0;Pcum如式(6):
其中,Pxk为转移概率矩阵P中第x行第k列的元素;;
④将随机数u与Pcum的第x行进行比较,若Pcum,x(y-1)<u≤Pcum,xy,则选取t1时刻为状态Xy;同时状态Xy对应的光伏输出功率Sy=Sy,min+u(Sy,max-Sy,min),Sy.min为状态Xy对应区间的最小功率,Sy.max为状态Xy对应区间的最大功率,随机数u取值范围在0到1之间。
步骤50、对包含M1个场景的原始场景集进行变速率削减:比较M1与设定最大能处理的场景序列数Nmax的大小,若M1大于Nmax则进行场景削减,使场景数削减至Nmax,否则返回步骤40,更新状态分布数M0为M1,模拟时间轴前进至t2时刻继续动态场景生成,直至模拟时间轴得到预设时刻tL。
步骤50中,对包含M1个场景的原始场景集进行变速率削减的具体步骤如下:
5.1、计算场景间二范数d(cv,cw):
5.2、确定削减速率V
由各场景最近场景组成集合D,按照概率距离从小到大排序,当前20%场景的平均概率距离小于后20%场景的平均概率距离10%时,采用大倍速削减,即V=0.01*M1;当前20%场景的平均概率距离大于后20%场景的平均概率距离10%且小于15%时,采用小倍速削减,即V=0.001*M1;当前20%场景的平均概率距离大于后20%场景的平均概率距离的15%时,采用逐个削减,即V=1;
其中,p(v)为M1个场景中第v个场景发生的概率,d(cv,cw)为场景cv与场景cw的二范数,cv、cw分别为M1个场景中第v、w个场景;
5.3、根据选定速率削减V个满足下式的场景cs*:
其中p(s)分别为场景cs*与场景cs发生的概率,cs*为当前M1个场景中被削减的第s*个场景,cs为当前M1个场景中除了场景cs*外的第s个场景,p(w)为场景cw发生的概率,cw为当前M1个场景中的第w个场景,p(z)为当前M1个场景中第z个场景发生的概率;该步骤删除了场景集中与其他场景概率距离最小的场景cs*,通过计算概率距离,可以将不具备代表性且发生概率很小的场景更快速地削减出场景集;
5.5、变更场景总数M1=M1-V*1,返回步骤5.2重复削减过程,直至削减后的场景数等于Nmax。
步骤60、将不超过Nmax个长度为L的场景序列进行场景削减,生成最终N0个的典型场景集。
以上步骤,形成了本分布式电源及节点负荷典型运行场景集生成策略,在生成初始静态场景基础上,动态生成连续时间断面运行场景,基于概率距离的变速率动态削减,快速有效削减呈指数增长的场景序列个数,得到能充分逼近整个样本空间的典型场景集,实现了配电网中分布式电源出力以及负荷需求的不确定性进行预测。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替代,都应当视为属于本发明的保护范围。
Claims (7)
1.一种分布式电源及节点负荷典型运行场景集生成方法,其特征在于,包括以下步骤:
步骤10、确定所要生成的分布式电源及节点负荷场景集的初始时刻为t0及该时刻的离散状态数为M0,设定最大能处理的场景序列数为Nmax,典型分布式电源及节点负荷场景数目为N0,预设目标时刻为tL;
步骤40、对由步骤30得到的M0个状态分布通过N1次马尔可夫链-蒙特卡洛方法MCMC抽样得到t1时刻的M1=M0*N1个场景序列,每个场景序列对应的发生概率为被抽样场景概率除以抽样次数;
步骤50、对包含M1个场景的原始场景集进行变速率削减:比较M1与设定最大能处理的场景序列数Nmax的大小,若M1大于Nmax则进行场景削减,使场景数削减至Nmax,否则返回步骤40,更新状态分布数M0为M1,模拟时间轴前进至t2时刻继续动态场景生成,直至模拟时间轴得到预设时刻tL;
步骤60、将不超过Nmax个长度为L的场景序列进行场景削减,生成最终N0个的典型场景集;
步骤50中,对包含M1个场景的原始场景集进行变速率削减的具体步骤如下:
5.1、计算场景间二范数d(cv,cw):
5.2、确定削减速率V
由各场景最近场景组成集合D,按照概率距离从小到大排序,当前20%场景的平均概率距离小于后20%场景的平均概率距离10%时,采用大倍速削减,即V=0.01*M1;当前20%场景的平均概率距离大于后20%场景的平均概率距离10%且小于15%时,采用小倍速削减,即V=0.001*M1;当前20%场景的平均概率距离大于后20%场景的平均概率距离的15%时,采用逐个削减,即V=1;
其中,p(v)为M1个场景中第v个场景发生的概率,d(cv,cw)为场景cv与场景cw的二范数,cv、cw分别为M1个场景中第v、w个场景;
5.3、根据选定速率削减V个满足下式的场景cs*;
其中p(s)分别为场景cs*与场景cs发生的概率,cs*为当前M1个场景中被削减的第s*个场景,cs为当前M1个场景中除了场景cs*外的第s个场景,p(w)为场景cw发生的概率,cw为当前M1个场景中的第w个场景,p(z)为当前M1个场景中第z个场景发生的概率;
5.5、变更场景总数M1=M1-V*1,返回步骤5.2重复削减过程,直至削减后的场景数等于Nmax。
2.根据权利要求1所述的一种分布式电源及节点负荷典型运行场景集生成方法,其特征在于,所述步骤10中,对于所确定的离散状态数M0,根据初始时刻t0的对应的月份和时段的历史功率数据来设置,其大小考虑转移概率矩阵P的阶数以及所生成模拟序列对原始数据的拟合度。
6.根据权利要求1所述的一种分布式电源及节点负荷典型运行场景集生成方法,其特征在于,步骤40中一次马尔可夫链-蒙特卡洛方法抽样具体步骤如下:
①从M0个状态中取一个状态Xx,x=1,2,...,M0;
②由标准均匀分布产生随机数u;
③根据转移概率矩阵P计算累计概率分布矩阵Pcum,矩阵Pcum中的第x行第y列的元素为Pcum,xy,其中x,y=1,2,...,M0;
④将随机数u与Pcum的第x行进行比较,若Pcum,x(y-1)<u≤Pcum,xy,则选取t1时刻为状态Xy;同时状态Xy对应的光伏输出功率Sy=Sy,min+u(Sy,max-Sy,min),Sy,min为状态Xy对应区间的最小功率,Sy,max为状态Xy对应区间的最大功率,随机数u取值范围在0到1之间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811067310.XA CN109301877B (zh) | 2018-09-13 | 2018-09-13 | 一种分布式电源及节点负荷典型运行场景集生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811067310.XA CN109301877B (zh) | 2018-09-13 | 2018-09-13 | 一种分布式电源及节点负荷典型运行场景集生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109301877A CN109301877A (zh) | 2019-02-01 |
CN109301877B true CN109301877B (zh) | 2021-08-17 |
Family
ID=65166762
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811067310.XA Active CN109301877B (zh) | 2018-09-13 | 2018-09-13 | 一种分布式电源及节点负荷典型运行场景集生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109301877B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110336333A (zh) * | 2019-07-31 | 2019-10-15 | 广东电网有限责任公司 | 一种区域综合能源系统的场景预测方法 |
CN110516944B (zh) * | 2019-08-20 | 2022-04-15 | 国网江苏省电力有限公司 | 一种配电网多阶段典型运行场景生成方法 |
CN110765582B (zh) * | 2019-09-24 | 2022-08-16 | 河海大学常州校区 | 基于马尔可夫链的自组织中心K-means微电网场景划分方法 |
CN111401755B (zh) * | 2020-03-19 | 2022-04-19 | 国电南瑞科技股份有限公司 | 基于马尔科夫链的多新能源出力场景生成方法、装置及系统 |
CN111401757B (zh) * | 2020-03-19 | 2022-04-22 | 西安交通大学 | 一种计及时空相关性的多区域电力负荷时序场景建模方法 |
CN112861365B (zh) * | 2021-02-23 | 2022-09-06 | 国网冀北电力有限公司智能配电网中心 | 低压配电网负荷状态估计方法、装置及终端设备 |
CN117293885B (zh) * | 2023-11-27 | 2024-03-22 | 华北电力大学 | 一种海岛微能网多目标调度优化方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012023816A (ja) * | 2010-07-12 | 2012-02-02 | Toshiba Corp | 情報処理装置およびそのプログラム |
CN105205502A (zh) * | 2015-10-30 | 2015-12-30 | 山东大学 | 一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法 |
CN106684889A (zh) * | 2017-03-24 | 2017-05-17 | 河海大学 | 一种基于场景法的主动配电网随机无功优化方法 |
CN106780134A (zh) * | 2016-12-13 | 2017-05-31 | 中国电力科学研究院 | 风电功率生成算法中风电状态数上限的确定方法及系统 |
CN107301472A (zh) * | 2017-06-07 | 2017-10-27 | 天津大学 | 基于场景分析法和电压调节策略的分布式光伏规划方法 |
CN107451349A (zh) * | 2017-07-25 | 2017-12-08 | 华中科技大学 | 一种用于生成风电功率序列的方法 |
CN107591800A (zh) * | 2017-09-04 | 2018-01-16 | 国网江苏省电力公司南京供电公司 | 基于场景分析的含分布式电源配电网运行状态预测方法 |
CN107947192A (zh) * | 2017-12-15 | 2018-04-20 | 西安理工大学 | 一种下垂控制型孤岛微电网的无功优化配置方法 |
-
2018
- 2018-09-13 CN CN201811067310.XA patent/CN109301877B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012023816A (ja) * | 2010-07-12 | 2012-02-02 | Toshiba Corp | 情報処理装置およびそのプログラム |
CN105205502A (zh) * | 2015-10-30 | 2015-12-30 | 山东大学 | 一种基于马尔柯夫蒙特卡罗的负荷特性综合分类方法 |
CN106780134A (zh) * | 2016-12-13 | 2017-05-31 | 中国电力科学研究院 | 风电功率生成算法中风电状态数上限的确定方法及系统 |
CN106684889A (zh) * | 2017-03-24 | 2017-05-17 | 河海大学 | 一种基于场景法的主动配电网随机无功优化方法 |
CN107301472A (zh) * | 2017-06-07 | 2017-10-27 | 天津大学 | 基于场景分析法和电压调节策略的分布式光伏规划方法 |
CN107451349A (zh) * | 2017-07-25 | 2017-12-08 | 华中科技大学 | 一种用于生成风电功率序列的方法 |
CN107591800A (zh) * | 2017-09-04 | 2018-01-16 | 国网江苏省电力公司南京供电公司 | 基于场景分析的含分布式电源配电网运行状态预测方法 |
CN107947192A (zh) * | 2017-12-15 | 2018-04-20 | 西安理工大学 | 一种下垂控制型孤岛微电网的无功优化配置方法 |
Non-Patent Citations (1)
Title |
---|
"含风电电力系统的场景分析方法及其在随机优化中的应用";马溪原;《中国博士学位论文全文数据库(电子期刊)》;20150630;第[21]-[29]页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109301877A (zh) | 2019-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109301877B (zh) | 一种分布式电源及节点负荷典型运行场景集生成方法 | |
CN108306303B (zh) | 一种考虑负荷增长和新能源出力随机的电压稳定评估方法 | |
CN111786417A (zh) | 一种面向分布式新能源消纳的主动配电网多目标区间优化调度的方法 | |
Tian et al. | Coordinated planning with predetermined renewable energy generation targets using extended two-stage robust optimization | |
WO2018177529A1 (en) | Method to accelerate the processing of multiperiod optimal power flow problems | |
Wang | A big data framework for stock price forecasting using fuzzy time series | |
Ramalakshmi | Optimal siting and sizing of distributed generation using fuzzy-EP | |
CN116388174A (zh) | 一种适应层次化调度的分布式光伏集群短期出力预测方法 | |
Sharma et al. | Distributed generation integration optimization using fuzzy logic controller. | |
CN113422371B (zh) | 一种基于图卷积神经网络的分布式电源就地电压控制方法 | |
Liu et al. | A novel learning method for traffic flow forecasting by seasonal SVR with chaotic simulated annealing algorithm | |
PADMA et al. | Application of fuzzy and ABC algorithm for DG placement for minimum loss in radial distribution system | |
CN108110756A (zh) | 考虑不确定性因素的工业园区配电网规划方法 | |
Korad et al. | Zonal do-not-exceed limits with robust corrective topology control | |
CN116722608A (zh) | 一种基于光伏逆变器的无功功率补偿系统 | |
Zakka et al. | Optimal capacitor placement using fuzzy and artificial bee colony algorithm for maximum loss reduction | |
CN116054179A (zh) | 一种基于事件触发的电力系统无功趋优控制系统及方法 | |
CN113824161B (zh) | 一种DoS攻击下的交直流微电网集群系统控制方法 | |
CN110829440B (zh) | 一种三相中低压一体化配电网电压控制方法及系统 | |
CN112348235A (zh) | 风光母线负荷自适应预测方法、装置、计算机设备 | |
Fernández-Fernández et al. | A dynamic programming approach for power curtailment decision making on pv systems | |
You | A Joint Chance Constrained Economic Dispatch Model Considering Wind Generation and Dynamic Line Rating | |
CN112801379B (zh) | 一种基于云计算与大数据的智能电网分布式能源管理系统 | |
CN116683472A (zh) | 无功功率补偿方法、装置、设备以及存储介质 | |
CN116683471A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |