CN115021335B - 一种考虑储能和新能源的多时段鲁棒可靠性评估方法 - Google Patents
一种考虑储能和新能源的多时段鲁棒可靠性评估方法 Download PDFInfo
- Publication number
- CN115021335B CN115021335B CN202210683294.7A CN202210683294A CN115021335B CN 115021335 B CN115021335 B CN 115021335B CN 202210683294 A CN202210683294 A CN 202210683294A CN 115021335 B CN115021335 B CN 115021335B
- Authority
- CN
- China
- Prior art keywords
- ith
- new energy
- energy storage
- storage element
- representing
- 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
- 238000004146 energy storage Methods 0.000 title claims abstract description 143
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000004088 simulation Methods 0.000 claims abstract description 39
- 230000009194 climbing Effects 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 7
- 238000012614 Monte-Carlo sampling Methods 0.000 claims description 4
- 239000013598 vector Substances 0.000 claims description 4
- 238000007599 discharging Methods 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 102000002274 Matrix Metalloproteinases Human genes 0.000 claims description 2
- 108010000684 Matrix Metalloproteinases Proteins 0.000 claims description 2
- 238000002948 stochastic simulation Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 8
- 230000006870 function Effects 0.000 description 10
- 238000011156 evaluation Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 4
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 230000001483 mobilizing effect Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
- H02J3/48—Controlling the sharing of the in-phase component
-
- 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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- 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
- H02J3/50—Controlling the sharing of the out-of-phase component
-
- 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/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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
- 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
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems 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)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种考虑储能和新能源的多时段鲁棒可靠性评估方法,涉及电网规划领域。方法具体包括:先获取包含储能元件在内的电网系统各元件工作参数和强迫停运率;然后采用随机模拟方法,抽取电网系统的故障元件,进而评估电网系统整体的鲁棒可靠性。在鲁棒可靠性计算中,采用电网系统的多时段鲁棒切负荷模型,基于考虑在新能源不确定性下的多时段火电机组限制、储能约束,计算模拟的电网系统的鲁棒切负荷量。本发明可用于准确评估在考虑新能源不确定性下的包含储能元件的电网系统的鲁棒可靠性。
Description
技术领域
本公开涉及电网规划领域,尤其涉及含有储能元件和新能源的新型电网系统的鲁棒可靠性评估。
背景技术
在能源革命和数字革命相融并进的背景下,以风电、光伏等新能源为代表的能源结构调整和能源高效变换是实现碳达峰、碳中和目标的重要途径。为实现大规模新能源场站的友好接入,促进新能源消纳,需要通过调动储能参与的方法构建有功功率相对可控的联合发电系统,实现包括新能源在内的多种调节资源的协调调控。全面快速的评估含有储能元件和新能源的新型电网系统可靠性是重大和紧迫的国家需求。
然而,现阶段对于含有储能元件和新能源的电网系统多时段最优潮流研究较少。此外,在实际应用中,由于新能源出力的不确定性带来的火电机组和储能元件的出力变化难以准确衡量,导致评估系统的多时段鲁棒可靠性的难度增大。
发明内容
针对上述现有技术,本发明所解决的技术问题在于如何克服新能源出力的不确定性带来的火电机组和储能元件的出力变化难以准确衡量,导致评估系统的多时段鲁棒可靠性的难度增大。
为了解决上述技术问题,本发明提出一种考虑储能和新能源的多时段鲁棒可靠性评估方法,所述方法包括下述步骤:
根据电网系统的拓扑结构、各元件强迫停运率和各元件在各时段的工作参数范围,建立模拟电网系统;
使模拟电网系统具有初始工作状态;
步骤S:在当前模拟电网系统下,根据各元件强迫停运率,采用随机模拟方法使元件故障,得到模拟电网系统新的拓扑结构;
基于新的拓扑结构,计算当前模拟电网系统的鲁棒可靠性指标值;
判断当前模拟电网系统的鲁棒可靠性指标值是否满足为最终结果的条件,若满足,则当前模拟电网系统的可靠性指标值为电网系统在当前初始工作状态下对应的可靠性指标;否则,返回步骤S;
所述电网系统包括储能元件、新能源组、火电机组;
所述当前模拟电网系统的储能元件不同时充放电。
在上述技术方案中,以实际电网系统中各时段为单位模拟电网系统,在当前模拟的电网系统中,同时含有储能元件和新能源的信息电网系统,且当前模拟电网系统不同时充放电,符合实际情况,从而使基于当前模拟电网系统评估的鲁棒可靠性能够真实反映实际电网系统的鲁棒可靠性,适应新型电网系统的计算需求。
在上述技术方案中,所述当前模拟电网系统满足下述电力网格约束:
其中:n表示母线总数;R代表所有新能源机组的集合;G代表所有火电机组的集合;T为模拟时长;表示不考虑新能源不确定性时第i个火电机组在t时刻的出力;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准放电功率;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准充电功率;Bij为导纳矩阵的第i行第j列的元素;/>表示不考虑新能源波动的t时刻第j个母线处的电压相角;/>为t时刻第i个母线的切负荷后的负荷大小;/>表示不考虑新能源不确定性时第i个新能源节点在t时刻的有功出力;/>表示表示不考虑新能源波动的t时刻第1个母线处的电压相角;/>为t时刻母线i的初始负荷;
所述当前模拟电网系统的一种实现方式是采用考虑新能源波动区间的直流潮流模型,以有效评估规划时段内不同时刻新能源的不确定性出力对电网系统的影响。在该模型中,第i个母线和第j个母线间的线路在t时刻流过的功率Pij,t表示为:
式中:bij为第i个母线和第j个母线间的线路电纳;为不考虑新能源波动的t时刻第i个母线处的电压相角;/>为不考虑新能源波动的t时刻第j个母线处的电压相角;/>表示第k个新能源节点在t时刻的新能源不确定出力;
为/>的第i行第k列元素;/>为/>的第j行第k列元素;
为以节点1为参考节点的直流潮流逆矩阵,如下所示:
在直流潮流逆矩阵中,/>为电网导纳矩阵去掉第一行第一列后的子矩阵;
δi,t为δt的第i个分量;δj,t为δt的第j个分量;其中:
αt为与相角有关的参与因子向量,其元素分量由与储能元件节点相关的参与因子和与火电机组节点相关的参与因子构成;
在最坏情况下,线路潮流应满足下述线路容量约束:
其中,表示第i母线和第j个母线间的线路容量;Ωi表示与第j个母线相连的母线集合;zij,k,t表示第i个母线和第j个母线间的线路在t时刻与/>对应的松弛变量,/>表示第i个母线和第j个母线间的线路在t时刻与/>对应的松弛变量;/>表示第k个新能源节点在t时刻的新能源不确定出力最大值;/>表示第k个新能源节点在t时刻的新能源不确定出力最小值。
在上述技术方案中,在规划时段内不同时刻新能源的不确定性出力情况下,对所述储能元件进行建模量化,采用下式描述输出电力:
其中,表示不考虑新能源不确定性时第i个储能元件在t时刻的基准放电功率;表示不考虑新能源不确定性时第i个储能元件在t时刻的基准充电功率;/>表示t时刻第i个储能元件节点考虑新能源不确定性后的总输出功率,S表示储能元件节点集合;T为模拟时长;/>表示第i个储能元件在t时刻的参与因子;/>表示第j个新能源节点在t时刻的新能源不确定出力;R代表所有新能源机组的集合;
所述储能元件的基本状态下应满足以下三式:
βi,t∈{0,1},i∈S,t=1,...,T
其中,表示第i个储能元件的最大充电率;/>表示第i个储能元件的最大放电功率,βi,t表示第i个储能元件在t时刻的工作状态,当其值为1时代表充电状态,当其值为0时代表放电状态,通过引入0-1变量保证储能元件不能同时充放电,更符合电网系统的实际情况,从而使鲁棒可靠性指标的评估更符合实际情况。
在上述技术方案中,在规划时段内不同时刻新能源的不确定性出力情况下,所述火电机组采用下述模型描述火电出力:
其中,表示考虑不确定性前的第i个火电机组在t时刻的出力,/>表示考虑不确定性后的第i个火电机组在t时刻的出力,/>表示第i个火电机组在t时刻的参与因子,G代表所有火电机组的集合;
所述新能源组对各时刻各新能源节点的不确定出力建模,采用下述区间数表述:
其中,表示考虑不确定性前的第i个新能源节点在t时刻的出力;/>表示考虑不确定性后的第i个新能源节点在t时刻的出力;/>表示第i个新能源节点在t时刻的新能源不确定出力;/>表示第i个新能源节点在t时刻的新能源不确定出力最大值;/>表示第i个新能源节点在t时刻的新能源不确定出力最小值。
基于上述建立的考虑新能源不确定性的储能元件模型和火电机组模型,适应了新型电网系统的计算需求,提供了一种建立最优切负荷模型以计算电网系统的切负荷的实施方式,从而实现鲁棒可靠性指标计算。
在上述技术方案中,所述当前模拟电网系统的初始工作状态通过下述任一方法获得:使电网系统正常工作的各元件或组件具有随机值,或者根据电网系统各元件或组件的实际值。
在上述技术方案中,所述随机模拟方法包括蒙特卡洛抽样法、重要性抽样法或等分散抽样法。
在上述技术方案中,所述条件为下述任一种:当前模拟电网系统的鲁棒可靠性指标值是否收敛或者当前是否已经达到设定模拟次数。
作为上述技术方案的进一步改进,将鲁棒可靠性指标进一步确定为鲁棒期望切负荷概率,和/或鲁棒期望切负荷量,即:鲁棒可靠性指标可以为其中之一或者同时使用两者。其中,将鲁棒期望切负荷概率记作LOLP,则有:
其中,P(Xm)表示第m次模拟的电网系统Xm的概率,Moff表示模拟电网系统中出现切负荷事件的样本合集;
将鲁棒期望切负荷量记作EDNS,则有:
式中:为t时刻第i个母线的初始负荷,/>为t时刻第i个母线的切负荷后的负荷大小;NL为负荷母线数量,T为模拟时长。
在上述技术方案中,所述当前模拟电网系统的鲁棒切负荷基于多时段鲁棒切负荷模型获取;
所述多时段鲁棒切负荷模型包括目标函数和约束条件;
所述目标函数如下:
式中:为t时刻第i个母线的初始负荷,/>为t时刻第i个母线的切负荷后的负荷大小;NL为负荷母线数量,T为模拟时长;
所述约束条件包括考虑新能源不确定性的储能元件约束和火电机组约束,通过储能元件约束和火电机组约束求解上述目标函数,进而实现有效评估规划时段内不同时刻新能源的不确定性出力对电网系统的影响。
所述储能元件约束如下:
储能元件在新能源波动的最差情况下也应满足如下充放电速率限制:
式中:表示不考虑新能源不确定性时第i个储能元件在t时刻的基准充电功率;表示不考虑新能源不确定性时第i个储能元件在t时刻的基准放电功率;/>表示第i个储能元件节点在t时刻的参与因子;/>表示第k个新能源节点在t时刻的新能源不确定出力最大值;/>表示第k个新能源节点在t时刻的新能源不确定出力最小值;/>为第i个储能元件的最大充电功率;/>为第i个储能元件的最大放电功率;S表示储能元件节点集合;T为模拟时长;
在最差情况下应满足如下容量限制:
其中,表示不考虑新能源波动的第i个储能元件在t时刻的电量;/>表示不考虑新能源波动的第i个储能元件在t-1时刻的电量;Ei,max表示第i个储能元件的储电量最大值;Ei,min表示第i个储能元件的储电量最小值;Ei,f表示第i个储能元件的初始值;ηc表示充电效率;ηd表示放电效率;Δt表示模拟的时间间隔;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准充电功率;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准放电功率;S表示储能元件节点集合;T为模拟时长;/>表示不考虑新能源波动的第i个储能元件在T时刻的电量;/>表示考虑新能源波动后第i个储能元件t时刻的实际最大电量;/>表示考虑新能源波动后第i个储能元件t时刻的实际最小电量;/>表示第j个新能源节点在t时刻的新能源不确定出力最小值;/>表示第j个新能源节点在t时刻的新能源不确定出力最大值;
所述火电机组约束如下:
火电机组在新能源出力的最差条件下也能满足出力限制如下:
其中,表示考虑新能源不确定性前的第i个火电机组在t时刻的出力;/>表示第i个火电机组在t时刻的参与因子;/>表示第i台火电机组最大出力上限,/>表示第i台火电机组最小出力下限;
火电机组的爬坡在最差条件下也应满足爬坡速率限制:
其中,表示考虑新能源不确定性前的第i个火电机组在t时刻的出力;/>表示考虑新能源不确定性前的第i个火电机组在t-1时刻的出力;/>表示第i个火电机组在t-1时刻的参与因子;/>表示第j个新能源节点在t时刻的新能源不确定出力最小值;/>表示第j个新能源节点在t时刻的新能源不确定出力最大值;/>表示第i台火电机组的最大爬坡速率;/>表示第i台火电机组的最小爬坡速率。
在上述约束条件中,考虑新能源不确定性出力下的最坏情况,降低了计算难度,同时又符合工程需求,能够为包含储能和新能源的新型电网系统的规划提供指导。
在上述技术方案中,参与因子一种方式是根据储能元件和火电机组的重要程度设定值,也可以根据有效运行时长所占比重设定值。作为上述技术方案的改进,所述参与因子满足下述约束,能够保证评估的鲁棒可靠性可能收敛,从而使评估值具有:
式中:S表示储能元件节点集合;T为模拟时长;G代表所有火电机组的集合;表示第i个火电机组在t时刻的参与因子;/>表示第j个储能元件节点在t时刻的参与因子。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1、为一种实施方式下的多时段鲁棒可靠性评估流程图;
图2、为一种实施方式下多时段鲁棒切负荷模型框架图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
本发明针对评估新型电网系统的鲁棒可靠性评估方法,在鲁棒可靠性评估中需要考虑电网系统在最差情况下满足的条件。本发明方法的一种实施方式下的多时段鲁棒可靠性评估流程图如图1所示,在图1中,系统指含有储能元件和新能源的电网系统。流程图的操作可以不按顺序实现。相反,操作可以以反转顺序或同时实现。此外,可以向流程图添加一个或多个其他操作。可以从流程图中移除一个或多个操作。
如图1所示,本发明方法的一种实施方式包括下述步骤:
S1:获取实际电网系统的拓扑结构、电网系统中各个元件的运行参数和强迫停运率,建立模拟电网系统;各个元件的运行参数包括新能源在各时刻的波动区间;
S2:初始化系统状态,设置用于计算当前模拟电网系统鲁棒可靠性指标值的迭代次数m=1;
S3:采用蒙特卡洛抽样法,根据各元件强迫停运率抽取系统故障元件;
S4:根据步骤S3的故障元件,修正当前模拟电网系统的拓扑结构,并计算当前模拟电网系统的鲁棒切负荷;
S5:计算系统的鲁棒期望切负荷量、鲁棒期望切负荷概率作为鲁棒可靠性指标值;
S6:判断迭代次数是否达到设定值或鲁棒可靠性指标是否收敛,若是,转入步骤S7;否则系统迭代次数m=m+1,转入步骤S3;
S7:鲁棒可靠性评估结束,将当前模拟电网系统的鲁棒可靠性指标值作为被模拟的实际电网系统的鲁棒可靠性。
上述评估电网系统鲁棒可靠性的方法,若把模拟的电网系统限定到某个时段,可以获得某个时段电网系统的鲁棒可靠性。
上述实际电网系统中包括储能元件、新能源组、火电机组;在模拟电网系统中,储能元件不同时放电。
在步骤3中,蒙特卡洛抽样法可用其它随机抽样方法替换,比如重要性抽样法或等分散抽样法等。
在步骤S3中,判断鲁棒可靠性指标是否收敛,可以将鲁棒期望切负荷概率是否收敛、或者鲁棒期望切负荷量是否收敛,或者两者均是否收敛作为判断条件。
在步骤S4中,当元件故障时,则该元件退出运行,对于发电机而言是出力为0;对于线路而言是线路开断,流过的功率为0;对于储能元件而言是该储能元件退出运行,与电网分离等。因此,修正当前模拟电网系统的拓扑结构,使除故障元件之外的元件重新构成可工作的模拟电网系统,该模拟电网系统可以功率不平衡,以使其更符合真实电网系统运行状态,从而使评估的鲁棒可靠性更准确地反映真实的电网系统性能。
在步骤S4中当前模拟电网系统的鲁棒切负荷的计算,本发明提供的一种实施方式是建立下述目标函数和约束条件,通过使目标函数取得最小值,从而获得当前模拟电网系统的鲁棒切负荷。
目标函数:
式中:为t时刻第i个母线的初始负荷,/>为t时刻第i个母线的切负荷后的负荷大小;NL为负荷母线数量,T为模拟时长;
约束条件参照图2所示:
对各时刻各新能源节点的不确定出力进行建模,采用区间数表示:
其中,R代表所有新能源机组的集合;表示考虑不确定性前的第i个新能源节点在t时刻的出力,/>表示考虑不确定性后的第i个新能源节点在t时刻的出力,/>表示第i个新能源节点在t时刻的新能源不确定出力最大值;/>表示第i个新能源节点在t时刻的新能源不确定出力最小值;/>表示第i个新能源节点在t时刻的新能源不确定出力。
①火电机组约束
考虑新能源出力不确定性后,火电出力可以表示为:
其中,分别表示考虑不确定性前后的第i个火电机组在t时刻的出力,/>表示第i个火电机组在t时刻的参与因子,G代表所有火电机组的集合。
考虑火电机组在新能源出力的最差条件下也能满足出力限制:
其中,表示考虑新能源不确定性前的第i个火电机组在t时刻的出力;/>表示第i个火电机组在t时刻的参与因子;/>分别表示第i台火电机组最大出力上限和最小出力下限。
类似的,火电机组的爬坡在最差条件下也应满足爬坡速率限制:
其中,表示考虑新能源不确定性前的第i个火电机组在t时刻的出力;/>表示考虑新能源不确定性前的第i个火电机组在t-1时刻的出力;/>表示第i个火电机组在t-1时刻的参与因子;/>表示第j个新能源节点在t时刻的新能源不确定出力最小值;/>表示第j个新能源节点在t时刻的新能源不确定出力最大值;/>表示第i台火电机组的最大爬坡速率;/>表示第i台火电机组的最小爬坡速率。
②储能元件约束
类似的,储能的输出电力可以表示为:
其中,表示不考虑新能源不确定性时第i个储能元件在t时刻的基准放电功率;表示不考虑新能源不确定性时第i个储能元件在t时刻的基准充电功率;/>表示t时刻第i个储能元件节点考虑新能源不确定性后的总输出功率,S表示储能元件节点集合;T为模拟时长;/>表示第i个储能元件在t时刻的参与因子;/>表示第j个新能源节点在t时刻的新能源不确定出力;R代表所有新能源机组的集合。
储能元件的基本状态下应满足以下三式:
βi,t∈{0,1},i∈S,t=1,...,T
其中,表示第i个储能元件的最大充电率;/>表示第i个储能元件的最大放电功率,βi,t表示第i个储能元件在t时刻的工作状态,当其值为1时代表充电状态,当其值为0时代表放电状态。
储能在新能源波动的最差情况下也应满足充放电速率限制,即
式中:表示不考虑新能源不确定性时第i个储能元件在t时刻的基准充电功率;表示不考虑新能源不确定性时第i个储能元件在t时刻的基准放电功率;/>表示第i个储能元件节点在t时刻的参与因子;/>表示第k个新能源节点在t时刻的新能源不确定出力最大值;/>表示第k个新能源节点在t时刻的新能源不确定出力最小值;/>为第i个储能元件的最大充电功率;/>为第i个储能元件的放电功率;S表示储能元件节点集合;T为模拟时长;
类似的,在最差情况下还应满足容量限制:
其中,表示不考虑新能源波动的第i个储能元件在t时刻的电量;/>表示不考虑新能源波动的第i个储能元件在t-1时刻的电量;Ei,max表示第i个储能元件的储电量最大值;Ei,min表示第i个储能元件的储电量最小值;Ei,f表示第i个储能元件的初始值;ηc表示充电效率;ηd表示放电效率;Δt表示模拟的时间间隔;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准充电功率;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准放电功率;S表示储能元件节点集合;T为模拟时长;/>表示不考虑新能源波动的第i个储能元件在T时刻的电量;/>表示考虑新能源波动后第i个储能元件t时刻的实际最大电量;/>表示考虑新能源波动后第i个储能元件t时刻的实际最小电量;/>表示第j个新能源节点在t时刻的新能源不确定出力最小值;/>表示第j个新能源节点在t时刻的新能源不确定出力最大值。
上述涉及的与储能元件相关的参与因子和与火电机组相关的参与因子可以为设定值,也可以是满足下述约束的值:
③参与因子的约束:
④电力网络的约束:
其中,其中:n表示母线总数;R代表所有新能源机组的集合;G代表所有火电机组的集合;T为模拟时长;表示不考虑新能源不确定性时第i个火电机组在t时刻的出力;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准放电功率;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准充电功率;Bij为导纳矩阵的第i行第j列的元素;/>表示不考虑新能源波动的t时刻第j个母线处的电压相角;/>为t时刻母线i的切负荷后的负荷大小;/>表示不考虑新能源不确定性时第i个新能源节点在t时刻的有功出力;/>表示表示不考虑新能源波动的t时刻第1个母线处的电压相角;/>为t时刻母线i的初始负荷。
δ1,t=0,i=2,...,n,t=1,...,T
/>
其中,为电网导纳矩阵去掉第一行第一列后的子矩阵;/>为第j个新能源节点在t时刻的不确定出力;Δθi,t为母线i在t时刻考虑新能源不确定性后电压相角的变化量;αt为与相角有关的参与因子向量,其元素分量由与储能元件节点相关的参与因子和与火电机组相关的参与因子构成;δj,t为δt的第j个分量,αi,t为αt的第i个分量;δt为与相角有关的辅助变量,αt为与相角有关的参与因子向量,其元素分量由与储能元件节点相关的参与因子和与火电机组节点相关的参与因子构成。
因此,采用直流潮流模型,考虑新能源波动区间后,母线i和母线j间的线路在t时刻流过的功率可以表示为:
式中:bij为第i个母线和第j个母线间的线路电纳;为不考虑新能源波动的t时刻第i个母线处的电压相角;/>为不考虑新能源波动的t时刻第j个母线处的电压相角;/>为表示第k个新能源节点在t时刻的新能源不确定出力。
在最坏情况下,线路潮流应满足线路容量约束,即:
其中,表示第i母线和第j个母线间的线路容量;Ωi表示与第j个母线相连的母线集合;zij,k,t表示第i个母线和第j个母线之间的线路在t时刻与/>对应的松弛变量,表示第i个母线和第j个母线之间的线路在t时刻与/>对应的松弛变量;/>表示第k个新能源节点在t时刻的新能源不确定出力最大值;/>表示第k个新能源节点在t时刻的新能源不确定出力最小值。
综上,通过上述目标函数在相关约束下取得最小值计算,能够获得步骤S4中的鲁棒切负荷,进而可计算鲁棒期望切负荷量EDNS、鲁棒期望切负荷概率LOLP如下:
式中:为t时刻第i个母线的初始负荷,/>为t时刻第i个母线的切负荷后的负荷大小;NL为负荷母线数量,T为模拟时长,Moff表示模拟电网系统中出现切负荷事件的样本合集。/>
式中:Xm为第m次模拟的电网系统,P(Xm)表示第m次模拟的电网系统Xm的概率。
在本发明方法的上述实施过程中,是对同时含有储能元件和新能源的新型电力系统进行鲁棒可靠性评估,引入0-1变量保证储能元件不能同时充放电,符合实际情况。在鲁棒可靠性评估中,相对可靠性评估,其需要考虑电网系统的最坏情况。最坏情况是指:考虑到电网系统中火电机组、储能元件和传输线路等各元件必须工作在安全范围内,且火电机组出力和储能元件的充放电功率、储电量变化速率不能过大的约束条件,当新能源出力取波动区间内的某一值时,某些约束刚好到达将要越限的临界情况。在本发明中,通过获取鲁棒切负荷进行鲁棒可靠性评估。在获取鲁棒切负荷时,包含考虑新能源不确定性的储能元件约束和火电机组约束,适应了新型电力系统的计算需求,能够有效评估规划时段内不同时刻新能源的不确定性出力对电力系统的影响。在储能元件约束和火电机组约束均考虑了在新能源不确定性出力下的最坏情况,降低了计算难度,同时又符合工程需求,能够为包含储能和新能源的新型电力系统的规划提供指导。
通过以上的实施方式的描述,所属领域的技术人员可以清楚地了解到本公开可借助软件加必需的通用硬件的方式来实现,当然也可以通过专用硬件包括专用集成电路、专用CPU、专用存储器、专用元器件等来实现。一般情况下,凡由计算机程序完成的功能都可以很容易地用相应的硬件来实现,而且,用来实现同一功能的具体硬件结构也可以是多种多样的,例如模拟电路、数字电路或专用电路等。但是,对本公开而言更多情况下,软件程序实现是更佳的实施方式。
尽管以上结合附图对本发明的实施方案进行了描述,但本发明并不局限于上述的具体实施方案和应用领域,上述的具体实施方案仅仅是示意性的、指导性的,而不是限制性的。本领域的普通技术人员在本说明书的启示下和在不脱离本发明权利要求所保护的范围的情况下,还可以做出很多种的形式,这些均属于本发明保护之列。
Claims (8)
1.一种考虑储能和新能源的多时段鲁棒可靠性评估方法,其特征在于,所述方法包括下述步骤:
根据电网系统的拓扑结构、各元件强迫停运率和各元件在各时段的工作参数范围,建立模拟电网系统;
使模拟电网系统具有初始工作状态;
步骤S:在当前模拟电网系统下,根据各元件强迫停运率,采用随机模拟方法使元件故障,得到模拟电网系统新的拓扑结构,并计算当前模拟电网系统的鲁棒切负荷;
基于新的拓扑结构,计算当前模拟电网系统的鲁棒可靠性指标值,鲁棒可靠性指标值为鲁棒期望切负荷概率和/或鲁棒期望切负荷量;
判断当前模拟电网系统的鲁棒可靠性指标值是否满足为最终结果的条件,若满足,则当前模拟电网系统的可靠性指标值为电网系统在当前初始工作状态下对应的可靠性指标;否则,返回步骤S;
所述电网系统包括储能元件、新能源组、火电机组;
所述当前模拟电网系统的储能元件不同时充放电;
所述当前模拟电网系统的鲁棒切负荷基于多时段鲁棒切负荷模型获取;
所述多时段鲁棒切负荷模型包括目标函数和约束条件;
所述目标函数如下:
式中:为t时刻第i个母线的初始负荷,/>为t时刻第i个母线的切负荷后的负荷大小;NL为负荷母线数量,T为模拟时长;
所述约束条件包括考虑新能源不确定性的储能元件约束和火电机组约束;
所述储能元件约束如下:
储能元件在新能源波动的最差情况下也应满足如下充放电速率限制:
式中:表示不考虑新能源不确定性时第i个储能元件在t时刻的基准充电功率;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准放电功率;/>表示第i个储能元件节点在t时刻的参与因子;/>表示第k个新能源节点在t时刻的新能源不确定出力最大值;/>表示第k个新能源节点在t时刻的新能源不确定出力最小值;/>为第i个储能元件的最大充电功率;/>为第i个储能元件的最大放电功率;S表示储能元件节点集合;T为模拟时长;
在最差情况下应满足如下容量限制:
其中,表示不考虑新能源波动的第i个储能元件在t时刻的电量;/>表示不考虑新能源波动的第i个储能元件在t-1时刻的电量;Ei,max表示第i个储能元件的储电量最大值;Ei,min表示第i个储能元件的储电量最小值;Ei,f表示第i个储能元件的初始值;ηc表示充电效率;ηd表示放电效率;Δt表示模拟的时间间隔;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准充电功率;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准放电功率;S表示储能元件节点集合;T为模拟时长;/>表示不考虑新能源波动的第i个储能元件在T时刻的电量;/>表示考虑新能源波动后第i个储能元件t时刻的实际最大电量;/>表示考虑新能源波动后第i个储能元件t时刻的实际最小电量;/>表示第j个新能源节点在t时刻的新能源不确定出力最小值;/>表示第j个新能源节点在t时刻的新能源不确定出力最大值;
所述火电机组约束如下:
火电机组在新能源出力的最差条件下也能满足出力限制如下:
其中,表示考虑新能源不确定性前的第i个火电机组在t时刻的出力;/>表示第i个火电机组在t时刻的参与因子;/>表示第i台火电机组最大出力上限,/>表示第i台火电机组最小出力下限;
火电机组的爬坡在最差条件下也应满足爬坡速率限制:
其中,表示考虑新能源不确定性前的第i个火电机组在t时刻的出力;/>表示考虑新能源不确定性前的第i个火电机组在t-1时刻的出力;/>表示第i个火电机组在t-1时刻的参与因子;/>表示第j个新能源节点在t时刻的新能源不确定出力最小值;/>表示第j个新能源节点在t时刻的新能源不确定出力最大值;/>表示第i台火电机组的最大爬坡速率;/>表示第i台火电机组的最小爬坡速率;G代表所有火电机组的集合;R代表所有新能源节点的集合。
2.根据权利要求1所述的方法,其特征在于,所述当前模拟电网系统满足下述电力网格约束:
其中:n表示母线总数;R代表所有新能源机组的集合;G代表所有火电机组的集合;T为模拟时长;表示不考虑新能源不确定性时第i个火电机组在t时刻的出力;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准放电功率;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准充电功率;Bij为导纳矩阵的第i行第j列的元素;/>表示不考虑新能源波动的t时刻第j个母线处的电压相角;/>为t时刻第i个母线的切负荷后的负荷大小;/>表示不考虑新能源不确定性时第i个新能源节点在t时刻的有功出力;/>表示表示不考虑新能源波动的t时刻第1个母线处的电压相角;/>为t时刻母线i的初始负荷;
所述当前模拟电网系统采用考虑新能源波动区间的直流潮流模型,则第i个母线和第j个母线间的线路在t时刻流过的功率Pij,t表示为:
式中:bij为第i个母线和第j个母线间的线路电纳;为不考虑新能源波动的t时刻第i个母线处的电压相角;/>为不考虑新能源波动的t时刻第j个母线处的电压相角;/>表示第k个新能源节点在t时刻的新能源不确定出力;
为/>的第i行第k列元素;/>为/>的第j行第k列元素;
为以节点1为参考节点的直流潮流逆矩阵,如下所示:
在直流潮流逆矩阵中,/>为电网导纳矩阵去掉第一行第一列后的子矩阵;
δi,t为δt的第i个分量;δj,t为δt的第j个分量;
其中:
αt为与相角有关的参与因子向量,其元素分量由与储能元件节点相关的参与因子和与火电机组节点相关的参与因子构成;
在最坏情况下,线路潮流应满足下述线路容量约束:
其中,表示第i母线和第j个母线间的线路容量;Ωi表示与第j个母线相连的母线集合;zij,k,t表示第i个母线和第j个母线间的线路在t时刻与/>对应的松弛变量,/>表示第i个母线和第j个母线间的线路在t时刻与/>对应的松弛变量;/>表示第k个新能源节点在t时刻的新能源不确定出力最大值;/>表示第k个新能源节点在t时刻的新能源不确定出力最小值。
3.根据权利要求1所述的方法,其特征在于,所述储能元件采用下式描述输出电力:
其中,表示不考虑新能源不确定性时第i个储能元件在t时刻的基准放电功率;/>表示不考虑新能源不确定性时第i个储能元件在t时刻的基准充电功率;/>表示t时刻第i个储能元件节点考虑新能源不确定性后的总输出功率,S表示储能元件节点集合;T为模拟时长;/>表示第i个储能元件在t时刻的参与因子;/>表示第j个新能源节点在t时刻的新能源不确定出力;R代表所有新能源机组的集合;
所述储能元件的基本状态下应满足以下三式:
βi,t∈{0,1},i∈S,t=1,...,T
其中,表示第i个储能元件的最大充电率;/>表示第i个储能元件的最大放电功率,βi,t表示第i个储能元件在t时刻的工作状态,当其值为1时代表充电状态,当其值为0时代表放电状态;
所述火电机组采用下述模型描述火电出力:
其中,表示考虑不确定性前的第i个火电机组在t时刻的出力,/>表示考虑不确定性后的第i个火电机组在t时刻的出力,/>表示第i个火电机组在t时刻的参与因子,G代表所有火电机组的集合;
所述新能源组对各时刻各新能源节点的不确定出力建模,采用下述区间数表述:
其中,表示考虑不确定性前的第i个新能源节点在t时刻的出力;/>表示考虑不确定性后的第i个新能源节点在t时刻的出力;/>表示第i个新能源节点在t时刻的新能源不确定出力;/>表示第i个新能源节点在t时刻的新能源不确定出力最大值;/>表示第i个新能源节点在t时刻的新能源不确定出力最小值。
4.根据权利要求1所述的方法,其特征在于,所述初始工作状态通过下述任一方法获得:使电网系统正常工作的各元件或组件具有随机值,或者根据电网系统各元件或组件的实际值。
5.根据权利要求1所述的方法,其特征在于,所述随机模拟方法为蒙特卡洛抽样法、重要性抽样法或等分散抽样法。
6.根据权利要求1所述的方法,其特征在于,所述条件为下述任一种:当前模拟电网系统的鲁棒可靠性指标值是否收敛或者当前是否已经达到设定模拟次数。
7.根据权利要求1所述的方法,其特征在于,将第m次模拟的电网系统记作Xm,使用P(Xm)表示第m次模拟的电网系统Xm的概率,使用Moff表示模拟电网系统中出现切负荷事件的样本合集;
若将鲁棒期望切负荷概率记作LOLP,则有:
若将鲁棒期望切负荷量记作EDNS,则有:
式中:为t时刻第i个母线的初始负荷,/>为t时刻第i个母线的切负荷后的负荷大小;NL为负荷母线数量,T为模拟时长。
8.根据权利要求1或3所述的方法,其特征在于,所述参与因子满足下述约束:
式中:S表示储能元件节点集合;T为模拟时长;G代表所有火电机组的集合;表示第i个火电机组在t时刻的参与因子;/>表示第j个储能元件节点在t时刻的参与因子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210683294.7A CN115021335B (zh) | 2022-06-16 | 2022-06-16 | 一种考虑储能和新能源的多时段鲁棒可靠性评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210683294.7A CN115021335B (zh) | 2022-06-16 | 2022-06-16 | 一种考虑储能和新能源的多时段鲁棒可靠性评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115021335A CN115021335A (zh) | 2022-09-06 |
CN115021335B true CN115021335B (zh) | 2024-03-26 |
Family
ID=83074101
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210683294.7A Active CN115021335B (zh) | 2022-06-16 | 2022-06-16 | 一种考虑储能和新能源的多时段鲁棒可靠性评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115021335B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005866A (zh) * | 2015-07-30 | 2015-10-28 | 国家电网公司 | 一种电力系统无功优化控制策略的鲁棒性评估方法 |
CN108183512A (zh) * | 2018-02-23 | 2018-06-19 | 南方电网科学研究院有限责任公司 | 一种接入新能源的电力系统的可靠性评估方法 |
CN109217294A (zh) * | 2018-09-14 | 2019-01-15 | 大唐新疆清洁能源有限公司 | 风电场大型电网安全性评估方法 |
CN109636027A (zh) * | 2018-12-07 | 2019-04-16 | 华电电力科学研究院有限公司 | 一种基于蒙特卡罗法的多能互补系统供能可靠性评估方法 |
CN111030094A (zh) * | 2019-12-05 | 2020-04-17 | 广东电网有限责任公司 | 一种多能互补园区供能可靠性评估方法 |
AU2020101660A4 (en) * | 2020-08-04 | 2020-09-10 | North China Electric Power University | Reliability calculation method of distribution system |
CN113901679A (zh) * | 2021-12-13 | 2022-01-07 | 中国南方电网有限责任公司超高压输电公司广州局 | 电力系统的可靠性分析方法、装置和计算机设备 |
-
2022
- 2022-06-16 CN CN202210683294.7A patent/CN115021335B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005866A (zh) * | 2015-07-30 | 2015-10-28 | 国家电网公司 | 一种电力系统无功优化控制策略的鲁棒性评估方法 |
CN108183512A (zh) * | 2018-02-23 | 2018-06-19 | 南方电网科学研究院有限责任公司 | 一种接入新能源的电力系统的可靠性评估方法 |
CN109217294A (zh) * | 2018-09-14 | 2019-01-15 | 大唐新疆清洁能源有限公司 | 风电场大型电网安全性评估方法 |
CN109636027A (zh) * | 2018-12-07 | 2019-04-16 | 华电电力科学研究院有限公司 | 一种基于蒙特卡罗法的多能互补系统供能可靠性评估方法 |
CN111030094A (zh) * | 2019-12-05 | 2020-04-17 | 广东电网有限责任公司 | 一种多能互补园区供能可靠性评估方法 |
AU2020101660A4 (en) * | 2020-08-04 | 2020-09-10 | North China Electric Power University | Reliability calculation method of distribution system |
CN113901679A (zh) * | 2021-12-13 | 2022-01-07 | 中国南方电网有限责任公司超高压输电公司广州局 | 电力系统的可靠性分析方法、装置和计算机设备 |
Non-Patent Citations (2)
Title |
---|
Alvaro Lorca et al..Multistage Robust Unit Commitment With Dynamic Uncertainty Sets and Energy Storage.《IEEE TRANSACTIONS ON POWER SYSTEMS》.2017,第32卷(第3期),第1678-1688页. * |
杨毅等.考虑安全性与可靠性的微电网电能优化调度 .《中国电机工程学报》.2014,第34卷(第19期),第3080-3088页. * |
Also Published As
Publication number | Publication date |
---|---|
CN115021335A (zh) | 2022-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Huo et al. | Chance-constrained optimization for multienergy hub systems in a smart city | |
Karaki et al. | Probabilistic performance assessment of autonomous solar-wind energy conversion systems | |
Li et al. | A multi-state model for the reliability assessment of a distributed generation system via universal generating function | |
Li et al. | Optimal planning of energy storage system in active distribution system based on fuzzy multi-objective bi-level optimization | |
Aien et al. | On possibilistic and probabilistic uncertainty assessment of power flow problem: A review and a new approach | |
CN104376504B (zh) | 一种基于解析法的配电系统概率可靠性评估方法 | |
CN108183512B (zh) | 一种接入新能源的电力系统的可靠性评估方法 | |
CN108400595B (zh) | 一种考虑新能源出力相关性的电压暂降随机预估方法 | |
Wang et al. | Optimal planning of stand-alone microgrids incorporating reliability | |
CN110350518B (zh) | 一种用于调峰的电网储能容量需求评估方法及系统 | |
CN115456304A (zh) | 考虑台风影响的海上风电场可靠性指标计算方法及装置 | |
CN111641205A (zh) | 一种基于随机优化的主动配电网故障管理方法 | |
CN113256152A (zh) | 基于贝叶斯网络的气电耦合系统可靠性评估方法及系统 | |
CN112257274B (zh) | 配电系统运行柔性的量化评估方法及系统 | |
Rouhani et al. | Soft clustering based probabilistic power flow with correlated inter temporal events | |
CN107179688B (zh) | 考虑蒙特卡罗状态抽样截尾的电力系统可靠性分析方法 | |
CN104283208B (zh) | 大型电网概率性可用输电能力的分解协调计算方法 | |
CN115021335B (zh) | 一种考虑储能和新能源的多时段鲁棒可靠性评估方法 | |
Li et al. | Identifying critical nodes in power grids containing renewable energy based on electrical spreading probability | |
CN109617100B (zh) | 一种数据驱动的风电场储能容量规划方法 | |
Bai et al. | Event-triggered forecasting-aided state estimation for active distribution system with distributed generations | |
CN115719967A (zh) | 一种用于提高供电可靠性的有源配电网储能装置优化配置方法 | |
CN113780722B (zh) | 配电网的联合规划方法、装置、计算机设备和存储介质 | |
Yundai et al. | Analysis of aggregated wind power dependence based on optimal vine copula | |
Niu et al. | Power supply reliability evaluation method for DC distribution network with power energy router |
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 |