CN105305439A - 一种考虑输入变量相关性的概率动态潮流计算方法及*** - Google Patents
一种考虑输入变量相关性的概率动态潮流计算方法及*** Download PDFInfo
- Publication number
- CN105305439A CN105305439A CN201510824118.0A CN201510824118A CN105305439A CN 105305439 A CN105305439 A CN 105305439A CN 201510824118 A CN201510824118 A CN 201510824118A CN 105305439 A CN105305439 A CN 105305439A
- Authority
- CN
- China
- Prior art keywords
- variable
- node
- power
- probability
- prime
- 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
Links
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种考虑输入变量相关性的概率动态潮流计算方法及***。考虑电力***的功频静特性,在常规潮流计算模型中增加频率待求变量,以建立概率动态潮流模型,可计算得到***频率的概率分布。此外,本发明考虑了输入变量相关性对概率潮流计算的影响。该模型更符合实际的电力***运行情况,从而为电力***的安全稳定分析提供更加完整和准确的综合评估。
Description
技术领域
本发明属于电力***稳态分析技术领域,更具体地,涉及一种考虑输入变量相关性的概率动态潮流计算方法及***。
背景技术
潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件,确定电力***各部分稳态运行状态参数的计算,是一种最基本、最重要、应用最广泛的电气计算,是电力***分析的基础。
在常规的潮流模型中,一般都忽略了发电机组和负荷的功频静态特性,将除平衡机以外的发电机有功出力设置为固定值,导致将***所有的不平衡功率分配给平衡节点,但这并不符合电力***的实际运行情况。为了能更准确地反映具有调频能力的发电机和负荷的有功功率随***频率变化而变化的特性,动态潮流的概念应运而生。传统的动态潮流是基于将不平衡功率按照一定的比例分配给具有调频能力的发电机,并在此基础上对电压再进行一次修正,但是这种算法收敛性差、求解速度慢。
此外,现代电力***中负荷和电源的不确定性也在逐步增大:用户侧负荷具有自主性和聚集效应,且负荷预测存在误差;随着近年来风电的大力发展,风电的波动性和不确定性同样需要重点关注。如何在动态潮流计算中考虑这些不确定性因素,并得到支路潮流、节点电压和***频率的概率分布结果也成为一个亟待解决的问题。
目前较为常用的概率潮流计算方法是蒙特卡洛模拟法、点估计法、半不变量法,其中,半不变量法由于计算速度快,且计算精度较高,在电力***实时分析评估中得到了广泛的应用。但该方法要求输入随机变量相互独立,这在实际电网中并不总能满足,比如相距较近的风电场出力往往具有较强的相关性。若在计算过程当中忽略了输入变量之间的相关性,则会影响所得支路潮流、节点电压和***频率的概率分布结果的准确性。因此研究计及输入变量相关性的概率动态潮流算法并探究相关性对计算结果的影响也十分必要。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种考虑输入变量相关性的概率动态潮流计算方法及***,建立了考虑电力***功频静态特性概率潮流模型,并计及输入变量的相关性,采用半不变量法进行概率潮流计算。该方法除了能得到节点电压和支路功率的概率分布,还能够得到***频率的概率分布,可以为电力***分析提供更加完整和准确的综合评估。
为实现上述目的,按照本发明的一个方面,提供了一种考虑输入变量相关性的概率动态潮流计算方法,包括如下步骤:(1)根据电网结构参数、发电机出力和输入随机变量的概率特性,建立常规概率潮流模型;(2)考虑发电机和负荷的功频静态特性,在常规概率潮流模型中增设频率待求变量,并增加平衡节点的节点有功功率方程,由此得到概率动态潮流模型中的修正方程;(3)根据概率动态潮流模型中的修正方程,计算电力***的基础运行点;基础运行点包括状态变量的基础运行点X0、输出变量的基础运行点Z0、状态变量对节点注入功率的灵敏度系数矩阵S0和输出变量对节点注入功率的灵敏度系数矩阵T0;(4)根据输入随机变量的概率特性,利用状态变量对节点注入功率的灵敏度系数矩阵S0和输出变量对节点注入功率的灵敏度系数矩阵T0,计算状态变量X和输出变量Z的各阶半不变量;(5)根据状态变量的基础运行点X0和输出变量的基础运行点Z0,以及状态变量X和输出变量Z的各阶半不变量,求取状态变量X和输出变量Z的累积概率分布函数和概率密度分布函数。
优选地,所述步骤(4)进一步包括如下步骤:(4-1)将输入随机变量W=[w1,w2,…,wl]T进行分块,得到W=[W′W″]T,其中,W′为相互独立的输入随机变量,W″为具有相关性的输入随机变量,其相关系数矩阵为CW;(4-2)对CW进行cholesky分解,得到CW=GGT,进而得到具有相关性的输入随机变量W″=AGY+μ,其中,G为下三角矩阵,A为对角矩阵,对角元素为W″中对应变量的标准差,Y为一组相互独立的服从标准正态分布的随机变量,μ为各个随机变量的期望;(4-3)将状态变量对节点注入功率的灵敏度系数矩阵S0和输出变量对节点注入功率的灵敏度系数矩阵T0进行分块,得到 其中,S′0为S0中与相互独立的输入随机变量对应的元素构成的子矩阵,S″0为S0中与具有相关性的输入随机变量对应的元素构成的子矩阵,T′0为T0中与相互独立的输入随机变量对应的元素构成的子矩阵,T″0为T0中与具有相关性的输入随机变量对应的元素构成的子矩阵;(4-4)计算状态变量X的各阶半不变量ΔX(k)和输出变量Z的各阶半不变量ΔZ(k)分别为: 其中,ΔW′(k)为W′的k阶半不变量,ΔY(k)为Y的k阶半不变量,S0 ′(k)、T0 ′(k)、(S0″AG)(k)和(T0″AG)(k)分别为S0′、T0′、T0″AG和S0″AG中元素的k次幂构成的矩阵。
优选地,修正方程具体为:
其中,ΔP为除平衡节点以外节点的有功功率修正量,ΔQ为pq节点的无功功率修正量,ΔPn为平衡节点的有功功率修正量,Δθ为除平衡节点以外节点的电压相角修正量,ΔV为pq节点的电压幅值修正量,Δf为***频率修正量, 为概率动态潮流模型的雅克比矩阵,H、N、K和L分别为非平衡节点有功功率不平衡量对节点电压相角的偏导数、非平衡节点有功功率不平衡量对节点电压幅值的偏导数、非平衡节点无功功率不平衡量对节点电压相角的偏导数和非平衡节点无功功率不平衡量对节点电压幅值的偏导数;C为发电机和负荷的功频静态特性向量,对C中的任一元素Ci,当其对应除平衡节点以外的发电机节点时,Ci=-KGi-KDi,当其对应负荷节点时,Ci=-KDi,KGi为除平衡节点以外的发电机节点连接的发电机的功频静态特性系数,KDi为非平衡节点连接的负荷的频率调节效应系数;Cn=-KGn-KDn为平衡节点的功频静态特性参数,KGn为平衡节点连接的发电机的功频静态特性系数,KDn为平衡节点连接的负荷的频率调节效应系数;Hn为平衡节点有功功率不平衡量对节点电压相角的偏导数,Nn为平衡节点有功功率不平衡量对节点电压幅值的偏导数。
优选地,状态变量对节点注入功率的灵敏度系数矩阵S0=J-1,输出变量对节点注入功率的灵敏度系数矩阵T0=G0S0,其中,G0为潮流基准点处支路功率对状态变量的偏导数矩阵。
按照本发明的另一方面,提供了一种考虑输入变量相关性的概率动态潮流计算***,其特征在于,包括:第一模块,用于根据电网结构参数、发电机出力和输入随机变量的概率特性,建立常规概率潮流模型;第二模块,用于考虑发电机和负荷的功频静态特性,在常规概率潮流模型中增设频率待求变量,并增加平衡节点的节点有功功率方程,由此得到概率动态潮流模型中的修正方程;第三模块,用于根据概率动态潮流模型中的修正方程,计算电力***的基础运行点;所述基础运行点包括状态变量的基础运行点X0、输出变量的基础运行点Z0、状态变量对节点注入功率的灵敏度系数矩阵S0和输出变量对节点注入功率的灵敏度系数矩阵T0;第四模块,用于根据输入随机变量的概率特性,利用状态变量对节点注入功率的灵敏度系数矩阵S0和输出变量对节点注入功率的灵敏度系数矩阵T0,计算状态变量X和输出变量Z的各阶半不变量;第五模块,用于根据状态变量的基础运行点X0和输出变量的基础运行点Z0,以及状态变量X和输出变量Z的各阶半不变量,求取状态变量X和输出变量Z的累积概率分布函数和概率密度分布函数。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下有益效果:考虑电力***的功频静特性,在常规潮流计算模型中增加频率待求变量,以建立概率动态潮流模型,可计算得到***频率的概率分布。此外,本发明考虑了输入变量相关性对概率潮流计算的影响。该模型更符合实际的电力***运行情况,从而为电力***的安全稳定分析提供更加完整和准确的综合评估。
附图说明
图1是本发明实施例的考虑输入变量相关性的概率动态潮流计算方法流程图;
图2是计及和忽略相关性的概率动态潮流所计算出的典型支路有功功率概率密度分布与蒙特卡洛模拟法得结果的对比图;
图3是计及和忽略相关性的概率动态潮流所计算出的***频率概率密度分布与蒙特卡洛模拟法得到结果的对比图;
图4是计及和忽略相关性的概率动态潮流所计算出的典型支路有功功率累积概率分布与蒙特卡洛模拟法得结果的对比图;
图5是计及和忽略相关性的概率动态潮流所计算出的***频率累积概率分布与蒙特卡洛模拟法得到结果的对比图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示,本发明实施例的考虑输入变量相关性的概率动态潮流计算方法包括如下步骤:
(1)根据电网结构参数、发电机出力和输入随机变量(包括间歇式能源出力和随机负荷)的概率特性,建立常规概率潮流模型;
其中,输入随机变量一般认为服从正态分布。这是由于输入随机变量通常为随机负荷和间歇式能源发电,随机负荷通常都采用正态分布模型,而在间歇式能源发电当中,风电占有很大的比例。而根据实际风电场统计数据表明,在风电功率的超短期预测中,用正态分布来描述风电功率的预测误差也是适用的。
(2)考虑发电机和负荷的功频静态特性,在常规概率潮流模型中增设频率待求变量,并增加平衡节点的节点有功功率方程,由此得到概率动态潮流模型中的修正方程;
在潮流模型中增加了一个方程和一个未知数,总方程数和总未知数数目相等,仍然可以求解。
修正方程具体为:
其中,ΔP为除平衡节点以外节点的有功功率修正量,ΔQ为pq节点的无功功率修正量,ΔPn为平衡节点的有功功率修正量,Δθ为除平衡节点以外节点的电压相角修正量,ΔV为pq节点的电压幅值修正量,Δf为***频率修正量, 为概率动态潮流模型的雅克比矩阵,H、N、K、L为常规概率潮流模型的雅克比矩阵中的分块子矩阵,分别为非平衡节点有功功率不平衡量对节点电压相角的偏导数、非平衡节点有功功率不平衡量对节点电压幅值的偏导数、非平衡节点无功功率不平衡量对节点电压相角的偏导数、非平衡节点无功功率不平衡量对节点电压幅值的偏导数;C为发电机和负荷的功频静态特性向量,对C中的任一元素Ci,当其对应除平衡节点以外的发电机节点时,Ci=-KGi-KDi,当其对应负荷节点时,Ci=-KDi,KGi为除平衡节点以外的发电机节点连接的发电机的功频静态特性系数,KDi为非平衡节点连接的负荷的频率调节效应系数;Cn=-KGn-KDn为平衡节点的功频静态特性参数,KGn为平衡节点连接的发电机的功频静态特性系数,KDn为平衡节点连接的负荷的频率调节效应系数;Hn为平衡节点有功功率不平衡量对节点电压相角的偏导数,Nn为平衡节点有功功率不平衡量对节点电压幅值的偏导数。
(3)利用输入随机变量的期望值,采用牛顿法计算电力***的基础运行点,包括状态变量的基础运行点X0(状态变量包括节点电压和***频率)、输出变量的基础运行点Z0(输出变量包括支路功率)、状态变量对节点注入功率的灵敏度系数矩阵S0=J-1和输出变量对节点注入功率的灵敏度系数矩阵T0=G0S0,其中,G0为潮流基准点处支路功率对状态变量的偏导数矩阵。
(4)采用半不变量法计算概率动态潮流:根据输入随机变量的概率特性求取其各阶半不变量,并根据历史统计数据求取输入随机变量的相关系数矩阵CW,通过对相关系数矩阵CW进行cholesky分解来修正灵敏度系数矩阵S0和T0,从而将具有相关性的输入随机变量转换为不相关的随机变量,再计算状态变量X和输出变量Z的各阶半不变量;
进一步包括如下步骤:
(4-1)将输入随机变量W=[w1,w2,…,wl]T进行分块,得到W=[W′W″]T,其中,W′为相互独立的输入随机变量,W″为具有相关性的输入随机变量,其相关系数矩阵为CW;
(4-2)对CW进行cholesky分解,得到CW=GGT,G为下三角矩阵,其元素可由 得到,从而具有相关性的输入随机变量W″就可表示为W″=AGY+μ,其中,A为对角矩阵,对角元素为W″中对应变量的标准差,Y为一组相互独立的服从标准正态分布的随机变量,μ为各个随机变量的期望;
Y=[y1,y2,…,yl]T为一组不相关的的服从标准正态分布的随机变量。不相关的变量Y的相关系数矩阵CY为单位矩阵I,故
CY=ρ(Y,YT)=ρ(BW,WTBT)=Bρ(W,WT)BT
=BCWBT=BGGTBT=(BG)(BG)T=I
由以上式子可得,取B=G-1,即Y=G-1W″,从而具有相关性的输入随机变量W″就可表示为一组不相关的服从标准正态分布的随机变量Y的线性表达式,即W″=AGY+μ,其中A为对角矩阵,对角元为相应变量的标准差,μ为相应变量的期望。
(4-3)将S0和T0进行分块,得到 其中,S′0为S0中与相互独立的输入随机变量对应的元素构成的子矩阵,S″0为S0中与具有相关性的输入随机变量对应的元素构成的子矩阵,T′0为T0中与相互独立的输入随机变量对应的元素构成的子矩阵,T″0为T0中与具有相关性的输入随机变量对应的元素构成的子矩阵;
(4-4)由式 计算状态变量X的各阶半不变量ΔX(k)和输出变量Z的各阶半不变量ΔZ(k),其中,ΔW′(k)为W′的k阶半不变量,ΔY(k)为Y的k阶半不变量,S0 ′(k)、T0 ′(k)、(S0″AG)(k)和(T0″AG)(k)分别为S0′、T0′、T0″AG和S0″AG中元素的k次幂构成的矩阵。
为能应用半不变量的可加性这一性质,要求输入变量之间相互独立。由于输入随机变量均服从正态分布,因此只需将具有相关性的随机变量W″转换为不相关的随机变量Y,即等价于Y为一组相互独立的随机变量。
(5)根据状态变量的基础运行点X0和输出变量的基础运行点Z0,以及状态变量X的各阶半不变量ΔX(k)和输出变量Z的各阶半不变量ΔZ(k),求取状态变量X和输出变量Z的累积概率分布函数和概率密度分布函数。
由于随机负荷和风电功率的预测误差在超短期预测范围内服从正态分布,由此求出的状态变量X(节点电压、***频率)和输出变量Z(支路功率)也近似服从正态分布,其期望即为基准状态X0、Z0,方差即为X、Z的二阶半不变量。
以下结合具体实施例对本发明的技术方案和技术效果作进一步说明。
在本发明一个实施例中,为了验证所述考虑输入变量相关性的概率动态潮流计算方法的有效性和准确性,在修改后的IEEE30节点测试***仿真。假定风电接入节点22和27,接入风机为双馈风机,采用恒功率因数控制(默认为1),因此风机出力不对频率波动做出响应。假定这两个节点的风电出力具有相关性,且相关系数矩阵为
设定***中常规发电机组均为火电机组,各节点发电机参数如表1所示。***有功负荷的频率调节效应系数设为KD=1.5。风电节点参数如表2所示。
表1
表2
1)与蒙特卡洛模拟算法进行对比
为了适应在线实时应用,本发明采用的概率潮流算法为半不变量法,以获得较快的计算速度。为验证该方法的准确性,将采用蒙特卡洛模拟法的计算结果作为对比。表3~7分别为本发明所述的半不变量法和采用蒙特卡洛模拟法所计算出的支路有功功率、支路无功功率、电压幅值、电压相角和频率的对比情况,其中,支路有功功率的期望值相对误差大部分都小于0.1%,最大为0.93%,标准差相对误差大部分小于0.5%,最大为0.94%。支路无功功率的期望值相对误差最大为1.17%,标准差相对误差最大为0.8%,电压幅值的期望值相对误差最大为0.009%,标准差相对误差最大为0.86%;电压相角的期望值相对误差最大为0.15%,标准差相对误差最大为0.33%;***频率的期望值相对误差为0.0022%,标准差相对误差为0.29%。图2给出了考虑相关性和不考虑相关性的概率动态潮流所计算出的典型支路有功功率的概率密度分布与蒙特卡洛模拟法得结果的对比图;图3给出了考虑相关性和不考虑相关性的概率动态潮流所计算出的***频率的概率密度分布与蒙特卡洛模拟法得结果的对比图;图4给出了考虑相关性和不考虑相关性的概率动态潮流所计算出的典型支路有功功率的累积概率分布与蒙特卡洛模拟法得结果的对比图;图5给出了考虑相关性和不考虑相关性的概率动态潮流所计算出的***频率的累积概率分布与蒙特卡洛模拟法得结果的对比图;图2至图5中的考虑相关性所对应的曲线与蒙特卡洛模拟给出的结果拟合情况良好。可见,本发明所提方法与蒙特卡洛模拟算法相比,计算误差均在可接受范围内,且计算时间比蒙特卡洛模拟法短,具有实时在线应用的意义。
表3
表4
表5
表6
表7
2)忽略输入变量相关性的影响
在以往的概率动态潮流计算中,往往没有考虑输入变量相关性的影响,但是根据仿真结果表明,若忽略输入变量相关性,将会对支路功率、电压甚至频率的计算结果产生影响。表8、表9分别给出了忽略相关性对支路有功和***频率的影响。可见,若忽略了输入变量间的相关性,对期望值的影响不大,但是对标准差影响较大。其中,支路功率的标准差最大相对误差超过了100%,频率的标准差最大相对误差也达到了25%,此外,图2至图5中的不考虑相关性所对应的曲线与蒙特卡洛模拟给出的结果曲线差别较大,也说明忽略输入变量相关性会给电力***的安全稳定分析带来较大偏差,因此在计算***概率动态潮流时,需要计及输入变量相关性。
表8
表9
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种考虑输入变量相关性的概率动态潮流计算方法,其特征在于,包括如下步骤:
(1)根据电网结构参数、发电机出力和输入随机变量的概率特性,建立常规概率潮流模型;
(2)考虑发电机和负荷的功频静态特性,在常规概率潮流模型中增设频率待求变量,并增加平衡节点的节点有功功率方程,由此得到概率动态潮流模型中的修正方程;
(3)根据概率动态潮流模型中的修正方程,计算电力***的基础运行点;所述基础运行点包括状态变量的基础运行点X0、输出变量的基础运行点Z0、状态变量对节点注入功率的灵敏度系数矩阵S0和输出变量对节点注入功率的灵敏度系数矩阵T0;
(4)根据输入随机变量的概率特性,利用状态变量对节点注入功率的灵敏度系数矩阵S0和输出变量对节点注入功率的灵敏度系数矩阵T0,计算状态变量X和输出变量Z的各阶半不变量;
(5)根据状态变量的基础运行点X0和输出变量的基础运行点Z0,以及状态变量X和输出变量Z的各阶半不变量,求取状态变量X和输出变量Z的累积概率分布函数和概率密度分布函数。
2.如权利要求1所述的考虑输入变量相关性的概率动态潮流计算方法,其特征在于,所述步骤(4)进一步包括如下步骤:
(4-1)将输入随机变量W=[w1,w2,…,wl]T进行分块,得到W=[W′W″]T,其中,W′为相互独立的输入随机变量,W″为具有相关性的输入随机变量,其相关系数矩阵为CW;
(4-2)对CW进行cholesky分解,得到CW=GGT,进而得到具有相关性的输入随机变量W″=AGY+μ,其中,G为下三角矩阵,A为对角矩阵,对角元素为W″中对应变量的标准差,Y为一组相互独立的服从标准正态分布的随机变量,μ为各个随机变量的期望;
(4-3)将状态变量对节点注入功率的灵敏度系数矩阵S0和输出变量对节点注入功率的灵敏度系数矩阵T0进行分块,得到 其中,S′0为S0中与相互独立的输入随机变量对应的元素构成的子矩阵,S″0为S0中与具有相关性的输入随机变量对应的元素构成的子矩阵,T′0为T0中与相互独立的输入随机变量对应的元素构成的子矩阵,T0″为T0中与具有相关性的输入随机变量对应的元素构成的子矩阵;
(4-4)计算状态变量X的各阶半不变量ΔX(k)和输出变量Z的各阶半不变量ΔZ(k)分别为:
其中,ΔW′(k)为W′的k阶半不变量,ΔY(k)为Y的k阶半不变量,S0′(k)、T0′(k)、(S0″AG)(k)和(T0″AG)(k)分别为S0′、T0′、T0″AG和S0″AG中元素的k次幂构成的矩阵。
3.如权利要求1或2所述的考虑输入变量相关性的概率动态潮流计算方法,其特征在于,修正方程具体为:
其中,ΔP为除平衡节点以外节点的有功功率修正量,ΔQ为pq节点的无功功率修正量,ΔPn为平衡节点的有功功率修正量,Δθ为除平衡节点以外节点的电压相角修正量,ΔV为pq节点的电压幅值修正量,Δf为***频率修正量, 为概率动态潮流模型的雅克比矩阵,H、N、K和L分别为非平衡节点有功功率不平衡量对节点电压相角的偏导数、非平衡节点有功功率不平衡量对节点电压幅值的偏导数、非平衡节点无功功率不平衡量对节点电压相角的偏导数和非平衡节点无功功率不平衡量对节点电压幅值的偏导数;C为发电机和负荷的功频静态特性向量,对C中的任一元素Ci,当其对应除平衡节点以外的发电机节点时,Ci=-KGi-KDi,当其对应负荷节点时,Ci=-KDi,KGi为除平衡节点以外的发电机节点连接的发电机的功频静态特性系数,KDi为非平衡节点连接的负荷的频率调节效应系数;Cn=-KGn-KDn为平衡节点的功频静态特性参数,KGn为平衡节点连接的发电机的功频静态特性系数,KDn为平衡节点连接的负荷的频率调节效应系数;Hn为平衡节点有功功率不平衡量对节点电压相角的偏导数,Nn为平衡节点有功功率不平衡量对节点电压幅值的偏导数。
4.如权利要求3所述的考虑输入变量相关性的概率动态潮流计算方法,其特征在于,状态变量对节点注入功率的灵敏度系数矩阵S0=J-1,输出变量对节点注入功率的灵敏度系数矩阵T0=G0S0,其中,G0为潮流基准点处支路功率对状态变量的偏导数矩阵。
5.一种考虑输入变量相关性的概率动态潮流计算***,其特征在于,包括:
第一模块,用于根据电网结构参数、发电机出力和输入随机变量的概率特性,建立常规概率潮流模型;
第二模块,用于考虑发电机和负荷的功频静态特性,在常规概率潮流模型中增设频率待求变量,并增加平衡节点的节点有功功率方程,由此得到概率动态潮流模型中的修正方程;
第三模块,用于根据概率动态潮流模型中的修正方程,计算电力***的基础运行点;所述基础运行点包括状态变量的基础运行点X0、输出变量的基础运行点Z0、状态变量对节点注入功率的灵敏度系数矩阵S0和输出变量对节点注入功率的灵敏度系数矩阵T0;
第四模块,用于根据输入随机变量的概率特性,利用状态变量对节点注入功率的灵敏度系数矩阵S0和输出变量对节点注入功率的灵敏度系数矩阵T0,计算状态变量X和输出变量Z的各阶半不变量;
第五模块,用于根据状态变量的基础运行点X0和输出变量的基础运行点Z0,以及状态变量X和输出变量Z的各阶半不变量,求取状态变量X和输出变量Z的累积概率分布函数和概率密度分布函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510824118.0A CN105305439B (zh) | 2015-11-24 | 2015-11-24 | 一种考虑输入变量相关性的概率动态潮流计算方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510824118.0A CN105305439B (zh) | 2015-11-24 | 2015-11-24 | 一种考虑输入变量相关性的概率动态潮流计算方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105305439A true CN105305439A (zh) | 2016-02-03 |
CN105305439B CN105305439B (zh) | 2017-12-12 |
Family
ID=55202370
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510824118.0A Active CN105305439B (zh) | 2015-11-24 | 2015-11-24 | 一种考虑输入变量相关性的概率动态潮流计算方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105305439B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105846437A (zh) * | 2016-05-16 | 2016-08-10 | 中国电力科学研究院 | 一种计及互动相关性的cholesky分解半不变量潮流计算方法 |
CN106295862A (zh) * | 2016-08-01 | 2017-01-04 | 中国电力科学研究院 | 一种配电网无功优化数据处理方法 |
CN106451455A (zh) * | 2016-08-29 | 2017-02-22 | 甘肃省电力公司风电技术中心 | 基于节点电压相关性的含分布式电源***随机潮流方法 |
CN107684423A (zh) * | 2016-08-05 | 2018-02-13 | 深圳先进技术研究院 | 一种胎儿心电分离方法及装置 |
CN107834562A (zh) * | 2017-11-16 | 2018-03-23 | 大连海事大学 | 基于Matlab矩阵运算的快速分解法系数矩阵计算法 |
CN108336740A (zh) * | 2018-02-06 | 2018-07-27 | 重庆大学 | 一种考虑外网不确定性和静态频率特性的等值概率潮流方法 |
CN109412145A (zh) * | 2018-10-16 | 2019-03-01 | 国网上海市电力公司 | 一种基于同步测量数据的主动配电网动态特性评估方法 |
CN114936454A (zh) * | 2022-05-18 | 2022-08-23 | 东北电力大学 | 一种基于半不变量法的交直流混联***概率潮流计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103825269A (zh) * | 2014-02-24 | 2014-05-28 | 华中科技大学 | 一种考虑电力***功频静特性的快速概率潮流计算方法 |
CN104410069A (zh) * | 2014-12-05 | 2015-03-11 | 国家电网公司 | 一种计及响应相关性的动态概率潮流计算方法 |
CN104485665A (zh) * | 2014-12-17 | 2015-04-01 | 河海大学 | 计及风速预测误差时空相关性的动态概率潮流计算方法 |
CN104638644A (zh) * | 2015-02-12 | 2015-05-20 | 河海大学 | 一种含风电场电力***动态随机最优潮流获取方法 |
-
2015
- 2015-11-24 CN CN201510824118.0A patent/CN105305439B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103825269A (zh) * | 2014-02-24 | 2014-05-28 | 华中科技大学 | 一种考虑电力***功频静特性的快速概率潮流计算方法 |
CN104410069A (zh) * | 2014-12-05 | 2015-03-11 | 国家电网公司 | 一种计及响应相关性的动态概率潮流计算方法 |
CN104485665A (zh) * | 2014-12-17 | 2015-04-01 | 河海大学 | 计及风速预测误差时空相关性的动态概率潮流计算方法 |
CN104638644A (zh) * | 2015-02-12 | 2015-05-20 | 河海大学 | 一种含风电场电力***动态随机最优潮流获取方法 |
Non-Patent Citations (3)
Title |
---|
徐彩凤等: "基于半不变量的电力***静态稳定性概率分析", 《浙江电力》 * |
石东源等: "计及输入变量相关性的半不变量法概率潮流计算", 《中国电机工程学院》 * |
陈雁等: "考虑输入变量相关性的概率潮流计算方法", 《中国电机工程学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105846437A (zh) * | 2016-05-16 | 2016-08-10 | 中国电力科学研究院 | 一种计及互动相关性的cholesky分解半不变量潮流计算方法 |
CN106295862A (zh) * | 2016-08-01 | 2017-01-04 | 中国电力科学研究院 | 一种配电网无功优化数据处理方法 |
CN106295862B (zh) * | 2016-08-01 | 2022-03-18 | 中国电力科学研究院 | 一种配电网无功优化数据处理方法 |
CN107684423A (zh) * | 2016-08-05 | 2018-02-13 | 深圳先进技术研究院 | 一种胎儿心电分离方法及装置 |
CN106451455A (zh) * | 2016-08-29 | 2017-02-22 | 甘肃省电力公司风电技术中心 | 基于节点电压相关性的含分布式电源***随机潮流方法 |
CN107834562A (zh) * | 2017-11-16 | 2018-03-23 | 大连海事大学 | 基于Matlab矩阵运算的快速分解法系数矩阵计算法 |
CN107834562B (zh) * | 2017-11-16 | 2019-12-03 | 大连海事大学 | 基于Matlab矩阵运算的快速分解法系数矩阵计算法 |
CN108336740A (zh) * | 2018-02-06 | 2018-07-27 | 重庆大学 | 一种考虑外网不确定性和静态频率特性的等值概率潮流方法 |
CN108336740B (zh) * | 2018-02-06 | 2021-04-20 | 重庆大学 | 一种考虑外网不确定性和静态频率特性的等值概率潮流方法 |
CN109412145A (zh) * | 2018-10-16 | 2019-03-01 | 国网上海市电力公司 | 一种基于同步测量数据的主动配电网动态特性评估方法 |
CN109412145B (zh) * | 2018-10-16 | 2022-03-29 | 国网上海市电力公司 | 一种基于同步测量数据的主动配电网动态特性评估方法 |
CN114936454A (zh) * | 2022-05-18 | 2022-08-23 | 东北电力大学 | 一种基于半不变量法的交直流混联***概率潮流计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105305439B (zh) | 2017-12-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105305439A (zh) | 一种考虑输入变量相关性的概率动态潮流计算方法及*** | |
CN103825269B (zh) | 一种考虑电力***功频静特性的快速概率潮流计算方法 | |
WO2019091356A1 (zh) | 一种概率潮流计算方法及*** | |
CN107133406A (zh) | 一种电力***静态电压稳定域边界的快速搜索方法 | |
CN109494724B (zh) | 基于lu分解的大电网戴维南等值参数在线辨识方法 | |
CN105048468B (zh) | 基于分布式计算的输配电网一体化电压稳定评估方法 | |
CN103810646B (zh) | 一种基于改进投影积分算法的有源配电***动态仿真方法 | |
CN106849078B (zh) | 一种计及死区非线性的电力***稳态频率预测方法 | |
CN105207204B (zh) | 一种计及一次调频不确定性的概率潮流分析方法 | |
CN104156542B (zh) | 一种基于隐式投影的有源配电***稳定性仿真方法 | |
CN105514971A (zh) | 一种适用于各种运行模式微电网的潮流计算方法 | |
CN110429648B (zh) | 考虑风速随机波动的小干扰稳定裕度概率评估方法 | |
CN103985058B (zh) | 一种基于改进多中心校正内点法的可用输电能力计算方法 | |
CN105391059A (zh) | 一种基于电流量测变换的分布式发电***状态估计方法 | |
CN104269867A (zh) | 一种节点扰动功率转移分布均衡度分析方法 | |
CN105634018B (zh) | 一种基于随机响应面法和内点法的含风电随机最优潮流求解方法 | |
CN111884259A (zh) | 一种考虑***小干扰稳定特性的场站级风电机组自适应等值方法 | |
CN105184027A (zh) | 一种基于交互式多模型算法的电力负荷建模方法 | |
Kontis et al. | Dynamic equivalencing of active distribution grids | |
Weng et al. | A geometric analysis of power system loadability regions | |
Hu et al. | Frequency prediction model combining ISFR model and LSTM network | |
Chen et al. | Quantitative stability analysis of complex nonlinear hydraulic turbine regulation system based on accurate calculation | |
Lin et al. | Uncertainty quantification in dynamic simulations of large-scale power system models using the high-order probabilistic collocation method on sparse grids | |
Wang et al. | A three-machine equivalent system frequency response model and its closed-form solution | |
CN110474323B (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 |