CN103066593A - 含多类型分布式电源的弱环配电网三相潮流计算方法 - Google Patents
含多类型分布式电源的弱环配电网三相潮流计算方法 Download PDFInfo
- Publication number
- CN103066593A CN103066593A CN2012105590580A CN201210559058A CN103066593A CN 103066593 A CN103066593 A CN 103066593A CN 2012105590580 A CN2012105590580 A CN 2012105590580A CN 201210559058 A CN201210559058 A CN 201210559058A CN 103066593 A CN103066593 A CN 103066593A
- Authority
- CN
- China
- Prior art keywords
- node
- type
- sequence
- delta
- power
- 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
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种含多类型分布式电源的弱环配电网三相潮流计算方法,首先,基于序分量计算方法,结合弱环配电网的特点和道路的回路分析法,在配电序网中提出一种有效的三相不平衡弱环配电网改进潮流算法。接着,对PQ、PQ(V)、PV和PI节点类型DG的迭代计算模型进行了详细的公式推导,并可非常简单地将其引入到所提出的潮流算法中。本发明能高效处理多类型DG和弱环网,其计算过程清晰,编程简单,容易实现,保留了面向支路的前推回代法计算速度快、收敛性稳定的优点。最后,通过IEEE-37母线测试算例验证了本发明的有效性和良好的收敛性,以及具有较强的处理多类型DG和环网的能力。
Description
技术领域
本发明涉及一种含多类型分布式电源的弱环配电网三相潮流计算方法,属于电力***分析与计算技术领域。
背景技术
配电网潮流计算是配电***分析的一项重要内容,它是对配电***规划设计和运行方式的合理性、可靠性及经济性进行定量分析的重要依据。一方面,配电网一般是环形设计,开环运行,但在配电网进行重构和优化分析时,需要计算闭环运行时的弱环配电网三相潮流;另一方面,随着分布式发电技术迅速发展,越来越多分布式电源(distributed generation,DG)接入配电网,给配电网的网络结构、功率损耗、电压分布和潮流计算产生了巨大的影响;首先,配电网从传统单电源***变成了多电源***,潮流的流向由单向变成了不定向;其次,传统配电网中一般仅包含2种节点类型:Vθ节点(平衡节点)和PQ节点;而随着各种分布式电源接入配电网,***中增加了新的节点类型:PQ(V)节点、PV节点和PI节点。因此,研究含多类型分布式电源的弱环配电网三相潮流计算方法是非常必要的。
发明内容
发明目的:本发明针对现有配电网三相潮流计算方法无法适用于含多类型分布式电源的三相不平衡弱环配电网的三相潮流计算,提出了一种含多类型分布式电源的弱环配电网三相潮流计算方法。
技术方案:一种含多类型分布式电源的弱环配电网三相潮流计算方法,包括如下步骤:
A. 建立各类型分布式电源(distributed generation,DG)的潮流计算模型
1)PQ节点类型DG
常规潮流计算中,PQ节点的注入有功和无功功率为给定值,但这一处理方式不能简单推广到处理PQ型DG。现有含分布式电源的弱环配电网潮流算法中均把PQ型DG简单的处理成***的负负荷。但三相对称DG在三相不平衡弱环配电网中,由于三相电压不对称,输出的三相功率并不相等。而且,考虑到发电机运行特性,负序和零序有功功率都可以认为等于0。因此,对于三相不平衡的弱环配电***,这样的处理方式就不够合理。故本发明认为PQ型DG输出恒定有功和无功功率为DG节点的正序有功和正序无功功率,即
式中:PDG和QDG分别为PQ型DG的有功和无功输出。
2)PQ(V)节点类型DG
该类型DG的处理方法类似于PQ型DG,不同之处在于,迭代过程中,需要根据最新DG节点正序电压迭代值不断更新DG节点正序无功功率,然后求出DG节点新的正序注入电流,开始下一次迭代。其计算模型为
式中取值有以下2种情况:
a. 采用无励磁调节能力同步发电机作为接口时,DG发出的无功功率为
式中:PDG、EDGq、Xd、I1,DG,i分别为DG机组的有功输出、空载电势、同步电抗、端电压。
b. 采用异步发电机的风机作为接口时,DG吸收的无功功率为
式中:PDG、U1,DG,i分别为DG的有功输出、机端电压;x为异步电机定子漏抗与转子漏抗之和;xp为异步电机励磁电抗与机端并联电容等效电抗。
3)PV节点类型DG
考虑到发电机运行特性和三相对称DG存在不对称运行状态情况,本发明认为PV节点类型DG输出恒定有功功率为DG节点正序有功功率、输出恒定电压幅值为DG节点正序电压幅值,但是其输出无功功率却是未知的。因此,问题的关键就是求出满足DG节点正序电压幅值与PV型DG电压额定值相等的DG节点正序无功功率。则有
针对PV节点类型DG,可采用开环阻抗矩阵(戴维南等值阻抗矩阵)来处理PV型DG节点,在一个含有nDG,PV个PV型DG的三相弱环配电网的正序网络中,若在每个PV型DG节点处开环后出现nDG,PV个开环点,则存在
ΔU1,DG=Z1,DGΔI1,DG
式中:ΔU1,DG、ΔI1,DG为开环点校正的正序电压、正序电流矩阵(nDG,PV×1阶);Z1,DG为从开环点看进去的戴维南等值阻抗矩阵(nDG,PV×nDG,PV阶)。
针对任一弱环配电网的正序网络中,从道路矩阵T1中把各PV型DG节点所对应行向量提取出来组成一个新的矩阵T1,DG,则有
把ΔU1,DG、ΔI1,DG、Z1,DG表示为
ΔU1,DG=Δe1,DG+jΔf1,DG
ΔI1,DG=Δc1,DG+jΔd1,DG
Z1,DG=R1,DG+jX1,DG
第k次迭代后,设第i个PV型DG节点正序电流改变量为,则其正序视在功率的改变量为
则该PV型DG节点正序有功功率增加量为
因为PV型DG节点正序有功功率为常数,所以,根据上式可以求出
而该PV型DG节点正序无功功率增量为
则可得出
而第k+1次迭代时,第i个PV型DG节点的正序无功功率为
然后求出DG节点新的注入正序电流,开始下一次迭代。当|ΔU1,DG|满足收敛精度时,停止迭代。
4)PI节点类型DG
考虑到类似的情况,本发明认为PI节点类型DG输出恒定有功功率为DG节点正序有功功率、输出恒定电流幅值为DG节点正序电流幅值。则有
相应的DG输出无功功率可按下式计算得出:
因此,第k+1次迭代时相应DG节点正序无功功率为
然后求出DG节点新的注入正序电流,开始下一次迭代。
针对第i个DG节点,其节点注入的正序电流可用下式计算:
在潮流迭代过程中,若PQ(V)、PV和PI型DG节点出现无功功率越界,则将其转换成PQ型DG节点处理,且QDG取各节点类型DG的无功上界或下界,然后重新计算。
B. 含多类型分布式电源的弱环配电网三相潮流算法流程如下
1)确定弱环配电网络的树,给节点编号,规定树中的根节点的编号最小,设为“0”,其余节点按其离根节点的远近来编号,离根节点越远的节点编号越大。而树支的编号则规定为取两端节点编号中的大者,连支的编号依次从树支的最大编号后加1。
2)确定弱环配电网拓扑结构参数,包括节点数,支路数,回路数。不妨设三相弱环配电网有N个节点、b条支路和l个独立回路,假设首节点(根节点)是电源且作为参考节点,则独立节点个数为n=N-1,独立支路条数b=n+l。其中,树支n条,连支l条。
3)确定分布式发电弱环配电网中DG接入情况,包括DG接入的节点类型和对应的位置(节点号),接入DG的总数量,各类型DG并网参数,各类型并网DG的各自数量。
4)设首节点是电源且作为参考节点,首节点三相电压相量矩阵为Uabc,0(3×1阶) ,各节点三相电压相量矩阵为Uabc,n(3n×1阶),在弱环配电***三序网络中,可以得出首节点的三序电压矩阵为U012,0=AUabc,0(3×1阶),各节点三序电压矩阵为U012,n (3n×1阶)。其中,令a=ej2π/3, , 。
5)计算各序网络参数Zs,b。Zs,b为基于支路i的序阻抗Zs,bi形成的对角阵(n×n阶),其中,下标s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型;支路i的三相阻抗为Zabc,bi,则有Z012,bi=AZabc,biA-1,其中, , 。
6)计算各序网的道路矩阵Ts和回-支关联矩阵Bs。其中,下标s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型。
7)计算各序网中回路序阻抗矩阵,以及其逆矩阵Ys,1=(Zs,1)-1。其中,下标s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型。
11)针对不同类型DG依据步骤A中对应的潮流计算模型分别计算出第i个DG节点在第k次迭代时的和,然后求出第k次迭代时第i个DG节点的注入正序电流,i=1,2,…,nDG。其中和分别为第i个DG节点在第k次迭代时的有功功率和无功功率,为第i个DG节点在第k-1次迭代时的节点电压,nDG为接入***的DG个数。
13)计算第k次迭代时的。其中,下标s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型。
有益效果:含多类型分布式电源的弱环配电网三相潮流计算方法,首先基于序分量计算方法,结合弱环配电网的特点和道路的回路分析法,在配电序网中提出一种有效的三相不平衡弱环配电网改进潮流计算方法,充分利用序分量法在处理三相不平衡***接入对称DG和PV节点时的优势;然后对PQ、PQ(V)、PV和PI节点类型DG的迭代计算模型进行了详细的公式推导,并可非常简单地引入到所提潮流计算程序中实现。本发明能高效处理多类型DG和弱环网,其计算过程清晰,编程简单,容易实现,保留了面向支路的前推回代法计算速度快、收敛性稳定的优点。最后,通过IEEE-37母线测试算例分析并验证了本发明的有效性和良好的收敛性,以及具有较强的处理多类型DG和环网的能力。
附图说明
图1为本发明实施例的流程图;
图2为IEEE-37母线含多类型分布式电源的三相不平衡弱环配电网***原理图。
具体实施方式
下面结合具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
含多类型分布式电源的弱环配电网三相潮流计算方法,具体包括如下步骤:
1.建立各类型分布式电源(distributed generation,DG)的潮流计算模型
(1)PQ节点类型DG
常规潮流计算中,PQ节点的注入有功和无功功率为给定值,但这一处理方式不能简单推广到处理PQ型DG。现有含分布式电源的弱环配电网潮流算法中均把PQ型DG简单的处理成***的负负荷。但三相对称DG在三相不平衡弱环配电网中,由于三相电压不对称,输出的三相功率并不相等。而且,考虑到发电机运行特性,负序和零序有功功率都可以认为等于0。因此,对于三相不平衡的弱环配电***,这样的处理方式就不够合理。故本发明认为PQ型DG输出恒定有功和无功功率为DG节点的正序有功和正序无功功率,即
式中:PDG和QDG分别为PQ型DG的有功和无功输出。
针对第i个PQ型DG节点,节点注入的正序电流可用下式计算:
(2)PQ(V)节点类型DG
该类型DG的处理方法类似于PQ型DG,不同之处在于,迭代过程中,需要根据最新DG节点正序电压迭代值不断更新DG节点正序无功功率,然后将其代入式(2)求出DG节点新的正序注入电流,开始下一次迭代。其计算模型为
式中取值有以下2种情况:
a. 采用无励磁调节能力同步发电机作为接口时,DG发出的无功功率为
式中:PDG、EDGq、Xd、U1,DG,i分别为DG机组的有功输出、空载电势、同步电抗、端电压。
b. 采用异步发电机的风机作为接口时,DG吸收的无功功率为
式中:PDG、U1,DG,i分别为DG的有功输出、机端电压;x为异步电机定子漏抗与转子漏抗之和;xp为异步电机励磁电抗与机端并联电容等效电抗。
(3)PV节点类型DG
考虑到发电机运行特性和三相对称DG存在不对称运行状态情况,本发明认为PV节点类型DG输出恒定有功功率为DG节点正序有功功率、输出恒定电压幅值为DG节点正序电压幅值,但是其输出无功功率却是未知的。因此,问题的关键就是求出满足DG节点正序电压幅值与PV型DG电压额定值相等的DG节点正序无功功率。
针对PV节点类型DG,可采用开环阻抗矩阵(戴维南等值阻抗矩阵)来处理PV型DG节点,在一个含有nDG,PV个PV型DG的三相弱环配电网的正序网络中,若在每个PV型DG节点处开环后出现nDG,PV个开环点,则存在
ΔU1,DG=Z1,DGΔI1,DG(6)
式中:ΔU1,DG、ΔI1,DG为开环点校正的正序电压、正序电流矩阵(nDG,PV×1阶);Z1,DG为从开环点看进去的戴维南等值阻抗矩阵(nDG,PV×nDG,PV阶)。
针对任一弱环配电网的正序网络中,从道路矩阵T1中把各PV型DG节点所对应行向量提取出来组成一个新的矩阵T1,DG,则有
把ΔU1,DG、ΔI1,DG、Z1,DG表示为
ΔU1,DG=Δe1,DG+jΔf1,DG(8)
ΔI1,DG=Δc1,DG+jΔd1,DG(9)
Z1,DG=R1,DG+jX1,DG(10)
则该PV型DG节点正序有功功率增加量为
因为PV型DG节点正序有功功率为常数,所以,代入式(14)得
而该PV型DG节点正序无功功率增量为
把式(15)代入式(17)可得
而第k+1次迭代时,第i个PV型DG节点的正序无功功率为
然后将其代入式(2)求出DG节点新的注入正序电流,开始下一次迭代。当|ΔU1,DG|满足收敛精度时,停止迭代。
(4)PI节点类型DG
考虑到类似的情况,本发明认为PI节点类型DG输出恒定有功功率为DG节点正序有功功率、输出恒定电流幅值为DG节点正序电流幅值。相应的DG输出无功功率可按下式计算得出:
因此,第k+1次迭代时相应DG节点正序无功功率为
并将其代入式(2)求出DG节点新的注入正序电流,开始下一次迭代。
在潮流迭代过程中,若PQ(V)、PV和PI型DG节点出现无功功率越界,则将其转换成PQ型DG节点处理,且QDG取各节点类型DG的无功上界或下界,然后重新计算。
2.含多类型分布式电源的弱环配电网三相潮流算法
针对具有N个节点、b条支路和l个独立回路的三相弱环配电网,假设首节点(根节点)是电源且作为参考节点,则独立节点个数为n=N-1,支路条数b=n+l。其中,树支n条,连支l条。
对于连通图中一棵选定的树,由于基本回路中仅包含一条连支,基本回路数等于连支数,基本回路-支路(下面简称“回-支”)关联矩阵B描述基本回路、树支、连支之间的联系。其中回-支关联矩阵B是一个l×b阶矩阵,假定连支支路的正方向都是从大号节点指向小号节点,基本回路的正方向与连支支路的正方向相同,如果支路j在回路j内,且二者方向相同,则B(j,k)=1,如果支路k在回路j内,且二者方向相反,则B(j,k)=-1,如果支路k不在回路j内,则B(j,k)=0。
一个节点的道路是指节点沿树到根所经过的路径上的支路集合,节点的道路强调的是路径上的支路,对于一个给定的树,节点的道路是唯一的、只由树支组成,可用道路-支路关联矩阵(简称道路矩阵)T描述。其中道路矩阵T是一个n×b阶矩阵,假定道路的正方向都是从电源点(即根节点)指向各节点,各树支支路正方向与道路正方向相同,如果支路k在道路i上,则T(i,k)=1,反之T(i,k)=0。道路矩阵T是一个稀疏矩阵,利用稀疏技术可以降低内存需求。
在弱环配电序网络中,定义节点注入序电流向量矩阵为Is,n(n×1阶),支路序电流向量矩阵为Is,b(b×1阶),连支序电流(也即回路序电流)向量矩阵为Is,l(l×1阶)。在解耦的各序网模型电路中,可以获得各序网络的道路矩阵和回-支关联矩阵分别为T0、T1、T2和B0、B1、B2,并依据KCL电流定律,支路序电流Is,b与节点注入序电流Is,n、回路序电流Is,l满足如下等式:
其中,s=0,1,2,分别表示序网络模型中的零序、正序和负序网络。
对任一弱环配电***序分量电路模型中,基于欧姆定律,支路特性约束为
Us,b=Zs,bIs,b(23)
依据KVL电压定律,基本回路约束为
BsUs,b=0(24)
其中,Us,b为配电网支路序电压矩阵(n×1阶);Zs,b为基于支路i的序阻抗Zs,bi形成的对角阵(n×n阶),s=0,1,2,分别表示序网络模型中的零序、正序和负序网络模型。
由式(22)和式(23)代入式(24)可得:
设电源节点三相电压相量矩阵为Uabc,0(3×1阶) ,各节点三相电压相量矩阵为Uabc,n(3n×1阶),则可以得出电源节点的三序电压相量矩阵为U012,0=A Uabc,0 (3×1阶),各节点三序电压相量矩阵为U012,n(3n×1阶),那么,在各序网络模型中,可知任一节点与电源节点的序电压差等于从此节点开始沿着该节点的道路到达电源节点所经支路的支路序电压之和,即(设1n=[1,1,…,1]T,共n个1;a=ej2π/3, , ):
Us,n=1nUs,0-ΔUs,n(30)
式(29)是含多类型分布式电源的弱环配电网三相潮流算法计算的核心,如图1所述,含多类型分布式电源的弱环配电网三相潮流计算流程具体步骤如下(k为迭代次数):
1)确定弱环配电网络的树,给节点编号,规定树中的根节点的编号最小,设为“0”,其余节点按其离根节点的远近来编号,离根节点越远的节点编号越大。而树支的编号则规定为取两端节点编号中的大者,连支的编号依次从树支的最大编号后加1。
2)确定弱环配电网拓扑结构参数,包括节点数,支路数,回路数。不妨设三相弱环配电网有N个节点、b条支路和l个独立回路,假设首节点(根节点)是电源且作为参考节点,则独立节点个数为n=N-1,独立支路条数b=n+l。其中,树支n条,连支l条。
3)确定分布式发电弱环配电网中DG接入情况,包括DG接入的节点类型和对应的位置(节点号),接入DG的总数量,各类型DG并网参数,各类型并网DG的各自数量。
4)设首节点是电源且作为参考节点,首节点三相电压相量矩阵为Uabc,0(3×1阶) ,各节点三相电压相量矩阵为Uabc,n(3n×1阶),在弱环配电***三序网络中,可以得出首节点的三序电压矩阵为U012,0=AUabc,0(3×1阶),各节点三序电压矩阵为U012,n (3n×1阶)。其中,令a=ej2π/3, , 。
5)计算各序网络参数Zs,b。Zs,b为基于支路i的序阻抗Zs,bi形成的对角阵(n×n阶),其中,下标s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型;支路i的三相阻抗为Zabc,bi,则有,其中, , 。
6)计算各序网的道路矩阵Ts和回-支关联矩阵Bs。其中,下标s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型。
8)给弱环配电网各节点三相电压赋初始值。其中En[E,E,…,E]T,共n个E,E为3×3单位矩阵。
11)针对不同类型DG依据步骤A中对应的潮流计算模型分别计算出第i个DG节点在第k次迭代时的和,然后求出第k次迭代时第i个DG节点的注入正序电流,i=1,2,…,nDG。其中和分别为第i个DG节点在第k次迭代时的有功功率和无功功率,为第i个DG节点在第k-1次迭代时的节点电压,nDG为接入***的DG个数。
算例分析
图2为IEEE-37母线三相不平衡配电网,变压器Δ-Δ接线方式,增加5条环路,在母线12、25、29、30和35接入5个DG***分别为DG-1、DG-2、DG-3、DG-4和DG-5,其单线图如图2所示。
在IEEE-37母线含多类型DG的弱环配电***中,5个DG***DG-1、DG-2、DG-3、DG-4和DG-5依次分别接入PQ、PQ(V)-1、PQ(V)-2、PV和PI节点类型DG***,各节点类型DG并网参数如表1所示,以及其无功界限设置如表2所示。
表1 算例中各节点类型DG并网参数
表2 算例中各节点类型DG无功界限
选择了表3所示的7种运行方式进行分析讨论,7种运行方式下潮流收敛情况如表3所示,以及7种运行方式下接入DG的无功越界及无功输出情况如表6所示。
表3 算例的7种运行方式及潮流收敛迭代次数
表4 算例中DG的无功越界及无功输出情况
由表3可见,与其他方式相比,方式6和方式7中出现了迭代次数增加近似一倍,这是因为:方式6和方式7中出现接入DG无功输出越界,潮流程序自动将无功越界DG转换成PQ型DG后重新进行潮流计算,故迭代次数会增加近似一倍。另外,从表4中也可以看出在方式6和方式7中接入的PV型DG出现无功输出越界。
从表3和表4中可知,环路闭合后收敛次数减少,且随着闭合环路的增加潮流程序收敛性稳定;***在接入DG及增加DG时,潮流的收敛性变化很小,但在DG出现无功越界时,迭代次数会有所增加,但并不影响潮流收敛;显示了本发明具有良好的收敛性和处理DG及环网能力。在方式7和方式6都同时接入了4种不同节点类型DG的情况下,潮流稳定收敛,表明了本发明具有较强的处理多种不同节点类型DG同时并网的能力和处理DG无功越界的能力。
Claims (1)
1.一种含多类型分布式电源的弱环配电网三相潮流计算方法,其特征在于,包括如下步骤:
步骤一,建立各类型分布式电源的潮流计算模型
1)PQ节点类型DG
将PQ节点类型DG输出恒定有功和无功功率确定为DG节点的正序有功和正序无功功率,即
式中:PDG和QDG分别为PQ节点类型DG的有功和无功输出;
2)PQ(V)节点类型DG
PQ(V)节点类型DG的处理方法类似于PQ节点类型DG,不同之处在于,迭代过程中,需要根据最新DG节点正序电压迭代值不断更新DG节点正序无功功率,然后求出DG节点新的正序注入电流,开始下一次迭代,其计算模型为
式中取值有以下2种情况:
a. 采用无励磁调节能力同步发电机作为接口时,DG发出的无功功率为
式中:PDG、EDGq、Xd、U1,DG,i分别为DG机组的有功输出、空载电势、同步电抗、端电压;
b. 采用异步发电机的风机作为接口时,DG吸收的无功功率为
式中:PDG、U1,DG,i分别为DG的有功输出、机端电压;x为异步电机定子漏抗与转子漏抗之和;xp为异步电机励磁电抗与机端并联电容等效电抗;
3)PV节点类型DG
考虑到发电机运行特性和三相对称DG存在不对称运行状态情况,将PV节点类型DG输出恒定有功功率确定为DG节点正序有功功率,将输出恒定电压幅值确定为DG节点正序电压幅值,求出满足DG节点正序电压幅值与PV型DG电压额定值相等的DG节点正序无功功率;则有
针对PV节点类型DG,采用开环阻抗矩阵来处理PV型DG节点,在一个含有nDG,PV个PV型DG的三相配电网的正序网络中,若在每个PV型DG节点处开环后出现nDG,PV个开环点,则存在
ΔU1,DG=Z1,DGΔI1,DG
式中:ΔU1,DG、ΔI1,DG为开环点校正的正序电压、正序电流矩阵;Z1,DG为从开环点看进去的戴维南等值阻抗矩阵;
针对任一弱环配电网的正序网络中,从道路矩阵T1中把各PV型DG节点所对应行向量提取出来组成一个新的矩阵T1,DG,则有
把ΔU1,DG、ΔI1,DG、Z1,DG表示为
ΔU1,DG=Δe1,DGjΔf1,DG
ΔI1,DG=Δc1,DGjΔd1,DG
Z1,DG=R1,DGjX1,DG
式中:Δe1,DG和Δf1,DG分别表示ΔU1,DG的实部和虚部矩阵;Δc1,DG和Δd1,DG分别表示ΔI1,DG的实部和虚部矩阵;R1,DG和X1,DG分别表示Z1,DG的实部和虚部矩阵;
第k次迭代后,设第i个PV型DG节点正序电流改变量为,则其正序视在功率的改变量为
则该PV型DG节点正序有功功率增加量为
因为PV型DG节点正序有功功率为常数,所以,根据上式可以求出
而该PV型DG节点正序无功功率增量为
则可得出
而第k+1次迭代时,第i个PV型DG节点的正序无功功率为
然后求出DG节点新的注入正序电流,开始下一次迭代;当|ΔU1,DG|满足收敛精度时,停止迭代;
4)PI节点类型DG
考虑到类似的情况,本发明认为PI节点类型DG输出恒定有功功率为DG节点正序有功功率、输出恒定电流幅值为DG节点正序电流幅值;则有
相应的DG输出无功功率可按下式计算得出:
因此,第k+1次迭代时相应DG节点正序无功功率为
然后求出DG节点新的注入正序电流,开始下一次迭代;
针对第i个DG节点,其节点注入的正序电流可用下式计算:
在潮流迭代过程中,若PQ(V)、PV和PI型DG节点出现无功功率越界,则将其转换成PQ型DG节点处理,且QDG取各节点类型DG的无功上界或下界,然后重新计算;
步骤二,含多类型分布式电源的弱环配电网三相潮流算法流程如下
1)确定弱环配电网络的树,给节点编号,规定树中的根节点的编号最小,设为“0”,其余节点按其离根节点的远近来编号,离根节点越远的节点编号越大;而树支的编号则规定为取两端节点编号中的大者,连支的编号依次从树支的最大编号后加1;
2)确定弱环配电网拓扑结构参数,包括节点数,支路数,回路数;不妨设三相弱环配电网有N个节点、b条支路和l个独立回路,假设首节点是电源且作为参考节点,则独立节点个数为n=N-1,独立支路条数b=n+l;其中,树支n条,连支l条;
3)确定分布式发电弱环配电网中DG接入情况,包括DG接入的节点类型和对应的位置,接入DG的总数量,各类型DG并网参数,各类型并网DG的各自数量;
4)设首节点是电源且作为参考节点,首节点三相电压相量矩阵为,各节点三相电压相量矩阵为Uabc,n,在弱环配电***三序网络中,可以得出首节点的三序电压矩阵为U012,0=AUabc,0,各节点三序电压矩阵为U012,n;其中,令a=ej2π/3, , ;
5)计算各序网络参数Zs,b;Zs,b为基于支路i的序阻抗Zs,bi形成的对角阵(n×n阶),其中,下标s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型,支路i的三相阻抗为Zabc,bi,则有Z012,bi=AZabc,biA-1,其中, , ;
6)计算各序网的道路矩阵Ts和回-支关联矩阵Bs;其中,下标s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型;
11)针对不同类型DG依据步骤A中对应的潮流计算模型分别计算出第i个DG节点在第k次迭代时的和,然后求出第k次迭代时第i个DG节点的注入正序电流,i=1,2,…,nDG;其中和分别为第i个DG节点在第k次迭代时的有功功率和无功功率,为第i个DG节点在第k-1次迭代时的节点电压,nDG为接入***的DG个数;
15)计算k次迭代时节点i三相电压相量,i=1,2,…,n;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210559058.0A CN103066593B (zh) | 2012-12-20 | 2012-12-20 | 含多类型分布式电源的弱环配电网三相潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210559058.0A CN103066593B (zh) | 2012-12-20 | 2012-12-20 | 含多类型分布式电源的弱环配电网三相潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103066593A true CN103066593A (zh) | 2013-04-24 |
CN103066593B CN103066593B (zh) | 2014-08-13 |
Family
ID=48109105
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210559058.0A Expired - Fee Related CN103066593B (zh) | 2012-12-20 | 2012-12-20 | 含多类型分布式电源的弱环配电网三相潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103066593B (zh) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103683284A (zh) * | 2013-12-26 | 2014-03-26 | 国家电网公司 | 一种含分布式电源的配电网三相不平衡潮流计算方法 |
CN103971026A (zh) * | 2014-05-28 | 2014-08-06 | 南京工程学院 | 主动配电网潮流计算通用方法 |
CN104158198A (zh) * | 2013-05-15 | 2014-11-19 | 株式会社日立制作所 | 配电网优化潮流控制装置和方法 |
WO2014187047A1 (zh) * | 2013-05-24 | 2014-11-27 | Ma Shanjuan | 一种遵守二维守恒定律能量方向的判别方法 |
CN104466956A (zh) * | 2014-12-12 | 2015-03-25 | 国家电网公司 | 一种分布式电源的潮流计算方法及装置 |
CN104578159A (zh) * | 2015-01-13 | 2015-04-29 | 国家电网公司 | 一种含分布式电源的配电网三相潮流修正方法 |
CN104917197A (zh) * | 2015-07-01 | 2015-09-16 | 中国电力科学研究院 | 一种并行计算主动配电网三相不平衡潮流的方法 |
CN108387818A (zh) * | 2018-01-23 | 2018-08-10 | 中国石油大学(华东) | 一种适用于含分布式电源树状配电网的故障测距方法 |
CN108471119A (zh) * | 2018-04-24 | 2018-08-31 | 长沙理工大学 | 含智慧社区的配电网三相不平衡动态潮流模型预测控制方法 |
CN108964057A (zh) * | 2017-05-18 | 2018-12-07 | 广东电网有限责任公司佛山供电局 | 一种基于改进前推回代法的配电网潮流计算方法 |
CN109149583A (zh) * | 2018-09-27 | 2019-01-04 | 西安石油大学 | 有源配电网简洁前推回代潮流计算方法 |
CN109390971A (zh) * | 2018-11-07 | 2019-02-26 | 国家电网有限公司 | 一种基于门当户对遗传算法的配电网多目标主动重构方法 |
CN110518603A (zh) * | 2019-08-22 | 2019-11-29 | 浙江大学 | 含环网的三相不平衡的中低压完整配电网潮流计算方法 |
CN110808597A (zh) * | 2019-11-06 | 2020-02-18 | 山东电力工程咨询院有限公司 | 主动配电网中考虑三相不平衡的分布式电源规划方法 |
CN110880760A (zh) * | 2019-11-29 | 2020-03-13 | 国网辽宁省电力有限公司葫芦岛供电公司 | 一种低压配电网线损与无功补偿协同优化方法 |
CN111725812A (zh) * | 2020-06-04 | 2020-09-29 | 海仿(上海)科技有限公司 | 大规模配电***的潮流计算方法、装置、设备和存储介质 |
CN112039076A (zh) * | 2020-05-21 | 2020-12-04 | 国网电力科学研究院有限公司 | 综合分布式电源与负荷的配电网潮流动态等值方法及*** |
CN113009276A (zh) * | 2021-03-04 | 2021-06-22 | 国网宁夏电力有限公司电力科学研究院 | 一种基于阻抗矩阵的智能配电网故障定位方法 |
CN117154735A (zh) * | 2023-08-30 | 2023-12-01 | 国网江苏省电力有限公司扬州供电分公司 | 一种适用于新能源发电台区的潮流计算方法、设备及介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012005210A (ja) * | 2010-06-15 | 2012-01-05 | Hitachi Ltd | 配電系統潮流シミュレーション装置、配電系統潮流シミュレーション方法およびそのプログラム |
CN102709907A (zh) * | 2012-05-29 | 2012-10-03 | 上海电力学院 | 一种含分布式电源的弱环配电网潮流计算方法 |
-
2012
- 2012-12-20 CN CN201210559058.0A patent/CN103066593B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012005210A (ja) * | 2010-06-15 | 2012-01-05 | Hitachi Ltd | 配電系統潮流シミュレーション装置、配電系統潮流シミュレーション方法およびそのプログラム |
CN102709907A (zh) * | 2012-05-29 | 2012-10-03 | 上海电力学院 | 一种含分布式电源的弱环配电网潮流计算方法 |
Non-Patent Citations (2)
Title |
---|
李红伟等: "一种基于回路分析法的弱环配电网潮流直接算法", 《电力***保护与控制》, vol. 39, no. 22, 16 November 2011 (2011-11-16), pages 6 - 10 * |
王良缘等: "基于序分量的辐射网三相潮流求解", 《电力自动化设备》, vol. 23, no. 1, 31 January 2003 (2003-01-31), pages 21 - 23 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104158198A (zh) * | 2013-05-15 | 2014-11-19 | 株式会社日立制作所 | 配电网优化潮流控制装置和方法 |
WO2014187047A1 (zh) * | 2013-05-24 | 2014-11-27 | Ma Shanjuan | 一种遵守二维守恒定律能量方向的判别方法 |
CN103683284A (zh) * | 2013-12-26 | 2014-03-26 | 国家电网公司 | 一种含分布式电源的配电网三相不平衡潮流计算方法 |
CN103683284B (zh) * | 2013-12-26 | 2016-08-17 | 国家电网公司 | 一种含分布式电源的配电网三相不平衡潮流计算方法 |
CN103971026A (zh) * | 2014-05-28 | 2014-08-06 | 南京工程学院 | 主动配电网潮流计算通用方法 |
CN103971026B (zh) * | 2014-05-28 | 2017-02-15 | 南京工程学院 | 主动配电网潮流计算通用方法 |
CN104466956B (zh) * | 2014-12-12 | 2018-06-01 | 国家电网公司 | 一种分布式电源的潮流计算方法及装置 |
CN104466956A (zh) * | 2014-12-12 | 2015-03-25 | 国家电网公司 | 一种分布式电源的潮流计算方法及装置 |
CN104578159A (zh) * | 2015-01-13 | 2015-04-29 | 国家电网公司 | 一种含分布式电源的配电网三相潮流修正方法 |
CN104917197A (zh) * | 2015-07-01 | 2015-09-16 | 中国电力科学研究院 | 一种并行计算主动配电网三相不平衡潮流的方法 |
CN104917197B (zh) * | 2015-07-01 | 2018-01-19 | 中国电力科学研究院 | 一种并行计算主动配电网三相不平衡潮流的方法 |
CN108964057A (zh) * | 2017-05-18 | 2018-12-07 | 广东电网有限责任公司佛山供电局 | 一种基于改进前推回代法的配电网潮流计算方法 |
CN108387818A (zh) * | 2018-01-23 | 2018-08-10 | 中国石油大学(华东) | 一种适用于含分布式电源树状配电网的故障测距方法 |
CN108471119A (zh) * | 2018-04-24 | 2018-08-31 | 长沙理工大学 | 含智慧社区的配电网三相不平衡动态潮流模型预测控制方法 |
CN108471119B (zh) * | 2018-04-24 | 2021-07-06 | 长沙理工大学 | 含智慧社区的配电网三相不平衡动态潮流模型预测控制方法 |
CN109149583B (zh) * | 2018-09-27 | 2022-04-05 | 西安石油大学 | 有源配电网简洁前推回代潮流计算方法 |
CN109149583A (zh) * | 2018-09-27 | 2019-01-04 | 西安石油大学 | 有源配电网简洁前推回代潮流计算方法 |
CN109390971A (zh) * | 2018-11-07 | 2019-02-26 | 国家电网有限公司 | 一种基于门当户对遗传算法的配电网多目标主动重构方法 |
CN109390971B (zh) * | 2018-11-07 | 2021-02-26 | 国家电网有限公司 | 一种基于门当户对遗传算法的配电网多目标主动重构方法 |
CN110518603A (zh) * | 2019-08-22 | 2019-11-29 | 浙江大学 | 含环网的三相不平衡的中低压完整配电网潮流计算方法 |
CN110518603B (zh) * | 2019-08-22 | 2020-12-08 | 浙江大学 | 含环网的三相不平衡的中低压完整配电网潮流计算方法 |
CN110808597A (zh) * | 2019-11-06 | 2020-02-18 | 山东电力工程咨询院有限公司 | 主动配电网中考虑三相不平衡的分布式电源规划方法 |
CN110880760A (zh) * | 2019-11-29 | 2020-03-13 | 国网辽宁省电力有限公司葫芦岛供电公司 | 一种低压配电网线损与无功补偿协同优化方法 |
CN112039076A (zh) * | 2020-05-21 | 2020-12-04 | 国网电力科学研究院有限公司 | 综合分布式电源与负荷的配电网潮流动态等值方法及*** |
CN111725812A (zh) * | 2020-06-04 | 2020-09-29 | 海仿(上海)科技有限公司 | 大规模配电***的潮流计算方法、装置、设备和存储介质 |
CN113009276A (zh) * | 2021-03-04 | 2021-06-22 | 国网宁夏电力有限公司电力科学研究院 | 一种基于阻抗矩阵的智能配电网故障定位方法 |
CN117154735A (zh) * | 2023-08-30 | 2023-12-01 | 国网江苏省电力有限公司扬州供电分公司 | 一种适用于新能源发电台区的潮流计算方法、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN103066593B (zh) | 2014-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103066593B (zh) | 含多类型分布式电源的弱环配电网三相潮流计算方法 | |
CN103094900B (zh) | 考虑相序混合法的分布式发电配电网三相潮流计算方法 | |
Ju et al. | An extension of FBS three-phase power flow for handling PV nodes in active distribution networks | |
CN106208102B (zh) | 一种基于辅助问题原理的主动配电网分布式无功优化方法 | |
CN102842908B (zh) | 含多变压器支路的配电网三相解耦潮流计算方法 | |
CN106446458B (zh) | 考虑分布式电源的弱环配电网潮流计算方法 | |
CN109802392B (zh) | 大规模配电网潮流计算方法及装置 | |
CN103353918A (zh) | 一种配电网三相牛顿拉夫逊法潮流计算初值选取方法 | |
CN103984822A (zh) | 一种基于全数字实时仿真装置的三相潮流实现方法 | |
CN102842907B (zh) | 基于道路矩阵的配电网三相解耦潮流计算方法 | |
CN111668843B (zh) | 基于相序混合法的低压配电网三相潮流计算方法 | |
CN106130025B (zh) | 一种基于回路分析法计算三相配电网潮流的线性方法 | |
CN103956735B (zh) | 一种分布式发电***的谐波潮流分析方法 | |
CN107121604A (zh) | 一种不对称短路电流直流分量衰减时间常数获取方法 | |
CN102891485A (zh) | 基于序分量法的弱环配电网三相解耦潮流计算方法 | |
CN110323751A (zh) | 一种含分布式电源的弱环配电网潮流计算方法 | |
CN109149583A (zh) | 有源配电网简洁前推回代潮流计算方法 | |
CN112688332B (zh) | 基于序分量坐标的混合交直流微电网潮流分析方法及*** | |
Olamaei et al. | An efficient method for load flow analysis of distribution networks including PV nodes | |
Yang et al. | Power flow calculation for unbalanced three-phase distribution network with DGs based on phase-sequence hybrid modeling | |
Liao et al. | Measurement-based dynamic equivalent modeling for small and medium hydropower generator group | |
Zhou et al. | Dynamic Model Validation and Harmonic Stability Analysis of Offshore Wind Power Plants | |
CN103928941A (zh) | 一种发电机聚合的解析求解方法 | |
Chittora et al. | Coherency based dynamic equivalencing of electric power system | |
Rusilawati | Determination of generator steady state stability limit using losses concept and RBFNN |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140813 Termination date: 20191220 |
|
CF01 | Termination of patent right due to non-payment of annual fee |