CN104917197B - 一种并行计算主动配电网三相不平衡潮流的方法 - Google Patents
一种并行计算主动配电网三相不平衡潮流的方法 Download PDFInfo
- Publication number
- CN104917197B CN104917197B CN201510379552.2A CN201510379552A CN104917197B CN 104917197 B CN104917197 B CN 104917197B CN 201510379552 A CN201510379552 A CN 201510379552A CN 104917197 B CN104917197 B CN 104917197B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- mtd
- sequence
- mtr
- 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 27
- 239000011159 matrix material Substances 0.000 claims abstract description 33
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 13
- 238000002347 injection Methods 0.000 claims description 16
- 239000007924 injection Substances 0.000 claims description 16
- 230000005540 biological transmission Effects 0.000 claims description 14
- 230000001360 synchronised effect Effects 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000016507 interphase Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 3
- 238000000819 phase cycle Methods 0.000 abstract 2
- 238000010586 diagram Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/50—Arrangements for eliminating or reducing asymmetry in polyphase networks
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供一种并行计算主动配电网三相不平衡潮流的方法,该方法包括:(1)基于补偿法建立主动配电网中各元件序分量模型;(2)由所述各元件序分量模型建立网络的正序、负序和零序拓扑图,生成正序、负序和零序节点导纳矩阵;(3)根据序分量潮流算法请求全网络的不平衡潮流。本发明中在建模中增加了虚拟线路,简化了非全相线路的处理,消除了变压器不同接线方式的影响,建立了实际网络中不同种类、不同接线方式负荷的序分量模型,精确度很高,正序网络可以根据网络的特点灵活选择算法,并且三序网络间可并行计算。
Description
技术领域
本发明涉及一种计算三相不平衡潮流的方法,具体涉及一种并行计算主动配电网三相不平衡潮流的方法。
背景技术
随着分布式电源(DG)在配网侧渗透率逐渐提高,具有主动调节与较强控制能力的主动配电网(ADN)应运而生,成为了未来智能配电网的主要发展模式。ADN的研究和实践涉及了众多技术领域,其中配电网快速仿真与模拟(DFSM)是其基础与核心功能之一,作为DFSM的基本仿真工具的潮流计算,要求能够快速、精确计算配网三相不平衡潮流,具备超实时仿真功能。
配电网具有不同于输电网的显著特点:辐射状结构或接近辐射状结构运行;阻抗比较大;三相参数不平衡;存在非全相线路;三相负荷不平衡,负荷类型和接线方式多样。所以配电网的潮流计算与输电网差别较大,尤其是加入分布式电源后,配网潮流就显得更加复杂多变。配网三相潮流问题是一组非线性方程的求解问题。
目前,对配网不平衡潮流计算有三种方法:相分量法、序分量法以及两者相结合的混合算法。相分量法比较直观,三相之间不能解耦,无法并行计算,潮流计算量较大,计算效率不高,在处理变压器支路时能力较弱。基于补偿的序分量法可以解耦,便于并行计算。但是,现有序分量法在处理非全相线路时较为繁杂,而且没有计及实际配网中不同种类三相负荷的接线方式,不能精确、快速的计算配网不平衡潮流。
发明内容
为了克服上述现有技术的不足,本发明提供一种并行计算主动配电网三相不平衡潮流的方法,该方法中增加了虚拟线路,简化了非全相线路的处理,并消除了变压器不同接线方式的影响。
为了实现上述发明目的,本发明采取如下技术方案:
一种并行计算主动配电网三相不平衡潮流的方法,所述方法包括如下步骤:
(1)基于补偿法建立主动配电网中各元件序分量模型;
(2)由所述各元件序分量模型建立网络的正序、负序和零序拓扑图,生成正序、负序和零序节点导纳矩阵;
(3)根据序分量潮流算法请求全网络的不平衡潮流。
优选的,所述步骤(1)中,所述各元件序分量模型包括非全相线路序分量模型和负荷序分量模型。
优选的,建立所述非全相线路序分量模型包括如下步骤:
步骤3-1、增加虚拟线路,满足增加虚拟线路后与原线路等值;
步骤3-2、增加虚拟线路后,配电网中线路都为三相线路,一段三相不平衡配电网线路km,采用π型等值电路表示,序分量模型满足:
式中,1,2,0分别表示正序、负序、零序,表示线路并联序导纳矩阵,表示并联自导纳,表示并联互导纳,其中i,j=0,1,2;i≠j;表示线路串联序导纳矩阵,表示串联自导纳,表示串联互导纳,其中i,j=0,1,2;i≠j,为k节点序电压,为m节点序电压,表示线路km上注入序电流;
步骤3-3、把矩阵非对角元分离出来,只保留矩阵对角元时,将后面分离出的两部分看为补偿量,此时,序分量之间已无耦合关系,
式中,分别为线路km上的正序、负序、零序注入电流, 分别为k(m)点正序、负序、零序电压,分别为m点正序、负序、零序电压,表示线路补偿注入序电流。
优选的,所述负荷序分量模型分为恒功率、恒电流和恒阻抗三种基本类型,配电负荷的连接方式包括接地星形或不接地三角形接法的三相平衡或不平衡负荷,建立所述负荷序分量模型包括如下步骤:
步骤4-1、恒功率负荷/恒电流负荷:
假设k点接有恒功率/恒电流负载,
当负荷为星型接线时,k节点三相注入电流为,
式中,Sx为k节点x相负荷,为k节点x相电压;
当负荷为三角型接线时,k节点三相注入电流为,
Ia,b,c=CIab,bc,ca
式中,Iab,bc,ca为k节点ab,bc,ca相间电流;
步骤4-2、恒阻抗负荷:
当负荷采用三角形接线方式时,通过星三角变换将三角形连接的阻抗变为星型连接的阻抗,恒阻抗负荷的相分量模型看为3×3阶的对角阵;
三相负荷导纳为
式中:yxx为x相的负荷导纳,为x相的负荷功率因数角;
序分量电流和电压存在以下关系:
式中,为k节点三序注入电流,为k节点三序节点电压,为三相负荷序导纳矩阵。
优选的,所述步骤(2)中,所述生成正序、负序和零序节点导纳矩阵包括:
正序导纳阵Y1中k行m列的元素为,
式中,A为与k节点连接的输电线的总数,B为与k节点连接的变压器总数,C为与k节点连接的恒阻抗负荷总数,A’为km节点间输电线路总数,B’为km节点间变压器总数,分别为输电线路的串联、并联正序自导纳,为变压器的正序导纳,为负荷的正序导纳;
负序和零序导纳阵Y2、Y0中k行m列的元素为,
式中,D为连接到k节点三相同步电机总数,分别为输电线路的串联、并联负序和零序自导纳,为变压器的负序和零序导纳,为负荷的负序和零序导纳,为三相同步电机的负序和零序导纳。
优选的,所述步骤(3)中,所述序分量潮流算法包括如下步骤:
步骤6-1、设节点k的正序、负序、零序电压初值为
步骤6-2、更新负序、零序网络中节点k的注入电流;
式中,为k节点第ct次迭代的负序、零序注入电流,为km线路上第ct次迭代的负序、零序补偿注入电流,为k节点负荷第ct次迭代的负序、零序补偿注入电流,ct为迭代次数;
步骤6-3、更新正序网络中节点k的注入功率;
式中,为分布式电源正序注入功率,为km配网线路上的正序补偿注入功率,为节点k的正序负荷注入功率;
步骤6-4、由于负序、零序网络为无源网,节点导纳矩阵和节点注入电流已由上面步骤求出,所以采用节点电压方程求得ct+1次各节点负序,零序电压
步骤6-5、正序网络中采用单相潮流算法,求得ct+1次各节点正序电压
步骤6-6、生成ct+1次各节点相电压
式中,α=ej120°,为第ct+1次各节点三相电压,为第ct+1次各节点正序、负序、零序的电压;
步骤6-7、判断收敛,令当ΔU≤ε时,结束迭代,当ΔU>ε时,更新ct=ct+1,自步骤6-2开始进行下一次迭代;
其中,ΔU为三相各节点电压变化量绝对值中的最大值,ε为收敛精度。
与现有技术相比,本发明的有益效果在于:
本发明中在建模中增加了虚拟线路,简化了非全相线路的处理,消除了变压器不同接线方式的影响,建立了实际网络中不同种类、不同接线方式负荷的序分量模型,精确度很高,正序网络可以根据网络的特点灵活选择算法,并且三序网络间可并行计算。
附图说明
图1是本发明提供的并行计算主动配电网三相不平衡潮流的方法的流程图
具体实施方式
下面结合附图对本发明作进一步详细说明。
如图1所示,本发明提供的一种并行计算主动配电网三相不平衡潮流的方法,该方法具体步骤如下:
步骤1、基于补偿法建立主动配电网中各元件序分量模型;
步骤2、由所述各元件序分量模型建立网络的正序、负序和零序拓扑图,生成正序、负序和零序节点导纳矩阵;
步骤3、根据序分量潮流算法请求全网络的不平衡潮流。
在步骤1中,非全相线路建模由以下方式确定。
对于只有两相、单相的非全相线路,它们的串联阻抗阵为3×3阶奇异矩阵,不能直接求逆,为了使两相、单相线路也能使用序分量法解耦,可增加虚拟线路,只需满足增加虚拟线路后与原线路等值即可。假设距离已知相足够远处存在虚拟线路,虚拟线路的串联电抗足够大,并联电容为0,由于虚拟的线路距离已知线路足够远,可认为虚拟线路与已知线路间互阻抗为0,虚拟线路后,串联阻抗阵变为3×3阶非奇异矩阵。
经虚拟线路后,配网中线路都为三相线路,一段三相不平衡配网线路km,可采用π型等值电路表示。序分量模型满足:
式中,1,2,0分别表示正序、负序、零序,表示线路并联序导纳矩阵,表示并联自导纳,表示并联互导纳(i,j=0,1,2;i≠j);表示线路串联序导纳矩阵,表示串联自导纳,表示串联互导纳(i,j=0,1,2;i≠j)。为k节点序电压,为m节点序电压,表示线路km上注入序电流。
由上式可以看出,序分量之间也是耦合的,当把矩阵非对角元分离出来,只保留矩阵对角元时,可以将后面分离出的两部分看为补偿量,此时,序分量之间已无耦合关系。
式中,分别为线路km上的正序、负序、零序注入电流, 分别为k(m)点正序、负序、零序电压,表示线路补偿注入序电流。
在步骤1中,负荷建模由以下方式确定。
电力***负荷可以分为恒功率、恒电流、恒阻抗三种基本类型,配电***负荷也用这三种类型表示,配电负荷的连接方式有多种,包括接地星形或不接地三角形接法的三相平衡或不平衡负荷,甚至有单相或两相接地负荷,后者可以看成是前者的特殊情况。
(1)恒功率负荷/恒电流负荷:
假设k点接有恒功率/恒电流负载。
当负荷为星型接线时,k节点三相注入电流为,
式中,Sx为k节点x相负荷,为k节点x相电压。
当负荷为三角型接线时,k节点三相注入电流为,
Ia,b,c=CIab,bc,ca
其中,Iab,bc,ca为k节点ab,bc,ca相间电流。
(2)恒阻抗负荷:
当负荷采用三角形接线方式时,可通过星三角变换将三角形连接的阻抗变为星型连接的阻抗,所以,恒阻抗负荷的相分量模型可看为3×3阶的对角阵。
三相负荷导纳为
式中:yxx为x相的负荷导纳,为x相的负荷功率因数角。
序分量电流和电压存在以下关系:
式中,为k节点三序注入电流,为k节点三序节点电压,为三相负荷序导纳矩阵。
在步骤2中,网络的正序、负序、零序节点导纳矩阵由以下方式确定。
正序导纳阵Y1中k行m列的元素为,
其中,A为与k节点连接的输电线的总数,B为与k节点连接的变压器总数,C为与k节点连接的恒阻抗负荷总数,A’为km节点间输电线路总数,B’为km节点间变压器总数。分别为输电线路的串联、并联正序自导纳,为变压器的正序导纳,为负荷的正序导纳。
负序(零序)导纳阵Y2、Y0中k行m列的元素为,
其中,D为连接到k节点三相同步电机总数。分别为输电线路的串联、并联负序(零序)自导纳,为变压器的负序(零序)导纳,为负荷的负序(零序)导纳,为三相同步电机的负序(零序)导纳。
在步骤3中,序分量潮流算法如下:
1)设节点k的正序、负序、零序电压初值为
2)更新负序、零序网络中节点k的注入电流;
其中,为k节点第ct次迭代的负序、零序注入电流,为km线路上第ct次迭代的负序、零序补偿注入电流,为k节点负荷第ct次迭代的负序、零序补偿注入电流,ct为迭代次数。
3)更新正序网络中节点k的注入功率;
其中,为分布式电源正序注入功率,为km配网线路上的正序补偿注入功率,为节点k的正序负荷注入功率。
4)由于负序、零序网络为无源网,节点导纳矩阵和节点注入电流已由上面步骤求出,所以可采用节点电压方程求得ct+1次各节点负序,零序电压
5)正序网络中采用单相潮流算法,可求得ct+1次各节点正序电压 为第ct+1次各节点正序、负序、零序的电压;
6)生成ct+1次各节点相电压
其中,α=ej120°,为第ct+1次各节点三相电压。
7)判断收敛,令当ΔU≤ε时,结束迭代,当ΔU>ε时,更新ct=ct+1,自第2步开始进行下一次迭代。
其中,ΔU为三相各节点电压变化量绝对值中的最大值,ε为收敛精度。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求范围当中。
Claims (4)
1.一种并行计算主动配电网三相不平衡潮流的方法,其特征在于,所述方法包括如下步骤:
(1)基于补偿法建立主动配电网中各元件序分量模型;
(2)由所述各元件序分量模型建立网络的正序、负序和零序拓扑图,生成正序、负序和零序节点导纳矩阵;
(3)根据序分量潮流算法请求全网络的不平衡潮流;
所述步骤(1)中,所述各元件序分量模型包括非全相线路序分量模型和负荷序分量模型;
建立所述非全相线路序分量模型包括如下步骤:
步骤3-1、增加虚拟线路,满足增加虚拟线路后与原线路等值;
步骤3-2、增加虚拟线路后,配电网中线路都为三相线路,一段三相不平衡配电网线路km,采用π型等值电路表示,序分量模型满足:
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
</mrow>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>Y</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mi>S</mi>
</msubsup>
<mn>2</mn>
</mfrac>
<msubsup>
<mi>V</mi>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Y</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mi>Z</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>V</mi>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mi>m</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
式中,1,2,0分别表示正序、负序、零序,表示线路并联序导纳矩阵,表示并联自导纳,表示并联互导纳,其中i,j=0,1,2;i≠j;表示线路串联序导纳矩阵,表示串联自导纳,表示串联互导纳,其中i,j=0,1,2;i≠j,为k节点序电压,为m节点序电压,表示线路km上注入序电流;
步骤3-3、把矩阵非对角元分离出来,只保留矩阵对角元时,将后面分离出的两部分看为补偿量,此时,序分量之间已无耦合关系,
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>I</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
</mrow>
<mn>1</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>I</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>I</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
</mrow>
<mn>0</mn>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>y</mi>
<mn>11</mn>
<mi>Z</mi>
</msubsup>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msubsup>
<mi>y</mi>
<mn>22</mn>
<mi>Z</mi>
</msubsup>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msubsup>
<mi>y</mi>
<mn>00</mn>
<mi>Z</mi>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>(</mo>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>V</mi>
<mi>k</mi>
<mn>1</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>V</mi>
<mi>k</mi>
<mn>2</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>V</mi>
<mi>k</mi>
<mn>0</mn>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>V</mi>
<mi>m</mi>
<mn>1</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>V</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>V</mi>
<mi>m</mi>
<mn>0</mn>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>y</mi>
<mn>11</mn>
<mi>S</mi>
</msubsup>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msubsup>
<mi>y</mi>
<mn>22</mn>
<mi>S</mi>
</msubsup>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msubsup>
<mi>y</mi>
<mn>00</mn>
<mi>S</mi>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>V</mi>
<mi>k</mi>
<mn>1</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>V</mi>
<mi>k</mi>
<mn>2</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>V</mi>
<mi>k</mi>
<mn>0</mn>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<msubsup>
<mi>&Delta;I</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
</mrow>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msubsup>
</mrow>
式中,分别为线路km上的正序、负序、零序注入电流, 分别为k点正序、负序、零序电压,分别为m点正序、负序、零序电压,表示线路补偿注入序电流。
2.根据权利要求1所述方法,其特征在于,所述负荷序分量模型分为恒功率、恒电流和恒阻抗三种基本类型,配电负荷的连接方式包括接地星形或不接地三角形接法的三相平衡或不平衡负荷,建立所述负荷序分量模型包括如下步骤:
步骤4-1、恒功率负荷/恒电流负荷:
假设k点接有恒功率/恒电流负载,
当负荷为星型接线时,k节点三相注入电流为,
<mrow>
<msup>
<mi>I</mi>
<mi>x</mi>
</msup>
<mo>=</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>S</mi>
<mi>x</mi>
</msup>
<mo>/</mo>
<msubsup>
<mi>V</mi>
<mi>k</mi>
<mi>x</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>*</mo>
</msup>
<mo>,</mo>
<mi>x</mi>
<mo>=</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
式中,Sx为k节点x相负荷,为k节点x相电压;
当负荷为三角型接线时,k节点三相注入电流为,
Ia,b,c=CIab,bc,ca
式中,Iab,bc,ca为k节点ab,bc,ca相间电流;
步骤4-2、恒阻抗负荷:
当负荷采用三角形接线方式时,通过星三角变换将三角形连接的阻抗变为星型连接的阻抗,恒阻抗负荷的相分量模型看为3×3阶的对角阵;
三相负荷导纳为
式中:yxx为x相的负荷导纳,为x相的负荷功率因数角;
序分量电流和电压存在以下关系:
<mrow>
<msubsup>
<mi>I</mi>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msubsup>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>Y</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mi>l</mi>
</msubsup>
<msubsup>
<mi>V</mi>
<mi>k</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msubsup>
</mrow>
式中,为k节点三序注入电流,为k节点三序节点电压,为三相负荷序导纳矩阵。
3.根据权利要求1所述方法,其特征在于,所述步骤(2)中,所述生成正序、负序和零序节点导纳矩阵包括:
正序导纳阵Y1中k行m列的元素为,
<mrow>
<msubsup>
<mi>y</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
</mrow>
<mn>1</mn>
</msubsup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>A</mi>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mn>11</mn>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mi>Z</mi>
</msubsup>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>A</mi>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mn>11</mn>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mi>S</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>B</mi>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>r</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>C</mi>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mn>11</mn>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mi>l</mi>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mi>m</mi>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msup>
<mi>A</mi>
<mo>&prime;</mo>
</msup>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mn>11</mn>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mi>Z</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msup>
<mi>B</mi>
<mo>&prime;</mo>
</msup>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>r</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>k</mi>
<mo>&NotEqual;</mo>
<mi>m</mi>
<mo>.</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中,A为与k节点连接的输电线的总数,B为与k节点连接的变压器总数,C为与k节点连接的恒阻抗负荷总数,A’为km节点间输电线路总数,B’为km节点间变压器总数,分别为输电线路的串联、并联正序自导纳,为变压器的正序导纳,为负荷的正序导纳;
负序和零序导纳阵Y2、Y0中k行m列的元素为,
<mrow>
<msubsup>
<mi>y</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>A</mi>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mn>22</mn>
<mrow>
<mo>(</mo>
<mn>00</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mi>Z</mi>
</msubsup>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>A</mi>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mn>22</mn>
<mrow>
<mo>(</mo>
<mn>00</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mi>S</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>B</mi>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>r</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>C</mi>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mn>22</mn>
<mrow>
<mo>(</mo>
<mn>00</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mi>l</mi>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>D</mi>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mn>22</mn>
<mrow>
<mo>(</mo>
<mn>00</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mi>D</mi>
</msubsup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mi>m</mi>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msup>
<mi>A</mi>
<mo>&prime;</mo>
</msup>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mn>22</mn>
<mrow>
<mo>(</mo>
<mn>00</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mi>Z</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msup>
<mi>B</mi>
<mo>&prime;</mo>
</msup>
</msubsup>
<msubsup>
<mi>y</mi>
<mrow>
<mn>22</mn>
<mrow>
<mo>(</mo>
<mn>00</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>r</mi>
</mrow>
</msubsup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>k</mi>
<mo>&NotEqual;</mo>
<mi>m</mi>
<mo>.</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中,D为连接到k节点三相同步电机总数,分别为输电线路的串联、并联负序和零序自导纳,为变压器的负序和零序导纳,为负荷的负序和零序导纳,为三相同步电机的负序和零序导纳。
4.根据权利要求1所述方法,其特征在于,所述步骤(3)中,所述序分量潮流算法包括如下步骤:
步骤6-1、设节点k的正序、负序、零序电压初值为
步骤6-2、更新负序、零序网络中节点k的注入电流;
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>c</mi>
<mi>t</mi>
</mrow>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&Sigma;&Delta;I</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
<mo>,</mo>
<mi>c</mi>
<mi>t</mi>
</mrow>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&Sigma;&Delta;I</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>c</mi>
<mi>t</mi>
</mrow>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msubsup>
</mrow>
式中,为k节点第ct次迭代的负序、零序注入电流,为km线路上第ct次迭代的负序、零序补偿注入电流,为k节点负荷第ct次迭代的负序、零序补偿注入电流,ct为迭代次数;
步骤6-3、更新正序网络中节点k的注入功率;
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>c</mi>
<mi>t</mi>
</mrow>
<mn>1</mn>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&Sigma;S</mi>
<mrow>
<mi>D</mi>
<mi>E</mi>
<mi>R</mi>
</mrow>
<mn>1</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&Sigma;&Delta;S</mi>
<mrow>
<mi>k</mi>
<mi>m</mi>
<mo>,</mo>
<mi>c</mi>
<mi>t</mi>
</mrow>
<mn>1</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&Sigma;&Delta;S</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>c</mi>
<mi>t</mi>
</mrow>
<mn>1</mn>
</msubsup>
</mrow>
式中,为分布式电源正序注入功率,为km配网线路上的正序补偿注入功率,为节点k的正序负荷注入功率;
步骤6-4、由于负序、零序网络为无源网,节点导纳矩阵和节点注入电流已由上面步骤求出,所以采用节点电压方程求得ct+1次各节点负序,零序电压
步骤6-5、正序网络中采用单相潮流算法,求得ct+1次各节点正序电压
步骤6-6、生成ct+1次各节点相电压
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<mi>c</mi>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msup>
<mi>T</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>V</mi>
<mrow>
<mi>c</mi>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msubsup>
</mrow>
式中,α=ej120°,为第ct+1次各节点三相电压,为第ct+1次各节点正序、负序、零序的电压;
步骤6-7、判断收敛,令当ΔU≤ε时,结束迭代,当ΔU>ε时,更新ct=ct+1,自步骤6-2开始进行下一次迭代;
其中,ΔU为三相各节点电压变化量绝对值中的最大值,ε为收敛精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510379552.2A CN104917197B (zh) | 2015-07-01 | 2015-07-01 | 一种并行计算主动配电网三相不平衡潮流的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510379552.2A CN104917197B (zh) | 2015-07-01 | 2015-07-01 | 一种并行计算主动配电网三相不平衡潮流的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104917197A CN104917197A (zh) | 2015-09-16 |
CN104917197B true CN104917197B (zh) | 2018-01-19 |
Family
ID=54086018
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510379552.2A Active CN104917197B (zh) | 2015-07-01 | 2015-07-01 | 一种并行计算主动配电网三相不平衡潮流的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104917197B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105720594B (zh) * | 2016-04-08 | 2017-12-05 | 黑龙江省电力科学研究院 | 纯有功三相不平衡负荷补偿容量的补偿方法 |
CN108134392B (zh) * | 2017-12-29 | 2021-11-26 | 国网天津市电力公司电力科学研究院 | 一种基于节点电流注入的单三相混合潮流计算方法 |
CN110880764B (zh) * | 2019-10-30 | 2023-05-23 | 华南理工大学 | 一种含逆变型分布式电源不平衡配电网的故障处理方法 |
CN111027161B (zh) * | 2019-11-20 | 2023-08-01 | 积成电子股份有限公司 | 一种大规模输电网潮流图自动生成方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102842908A (zh) * | 2012-09-11 | 2012-12-26 | 河海大学 | 含多变压器支路的配电网三相解耦潮流计算方法 |
CN102842907A (zh) * | 2012-09-11 | 2012-12-26 | 河海大学 | 基于道路矩阵的配电网三相解耦潮流计算方法 |
CN103066593A (zh) * | 2012-12-20 | 2013-04-24 | 河海大学 | 含多类型分布式电源的弱环配电网三相潮流计算方法 |
CN103094900A (zh) * | 2012-12-20 | 2013-05-08 | 河海大学 | 考虑相序混合法的分布式发电配电网三相潮流计算方法 |
-
2015
- 2015-07-01 CN CN201510379552.2A patent/CN104917197B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102842908A (zh) * | 2012-09-11 | 2012-12-26 | 河海大学 | 含多变压器支路的配电网三相解耦潮流计算方法 |
CN102842907A (zh) * | 2012-09-11 | 2012-12-26 | 河海大学 | 基于道路矩阵的配电网三相解耦潮流计算方法 |
CN103066593A (zh) * | 2012-12-20 | 2013-04-24 | 河海大学 | 含多类型分布式电源的弱环配电网三相潮流计算方法 |
CN103094900A (zh) * | 2012-12-20 | 2013-05-08 | 河海大学 | 考虑相序混合法的分布式发电配电网三相潮流计算方法 |
Non-Patent Citations (3)
Title |
---|
一种基于序分量解耦的谐波潮流算法;喻荟铭;《福建电脑》;20090101;第58-59页 * |
不对称三相潮流的对称分量分析法;***;《中国电机工程学报》;19931130;第13卷(第6期);第1-12页 * |
基于序分量的弱环配电网三相解耦潮流算法;杨雄;《中国电机工程学报》;20140305;第34卷(第7期);第1-3节 * |
Also Published As
Publication number | Publication date |
---|---|
CN104917197A (zh) | 2015-09-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ju et al. | An extension of FBS three-phase power flow for handling PV nodes in active distribution networks | |
CN102609598B (zh) | 一种对大型电力***进行电磁暂态仿真的方法 | |
CN107563779B (zh) | 一种节点边际电价求解方法 | |
CN104917197B (zh) | 一种并行计算主动配电网三相不平衡潮流的方法 | |
CN103094900B (zh) | 考虑相序混合法的分布式发电配电网三相潮流计算方法 | |
CN104578159B (zh) | 一种含分布式电源的配电网三相潮流修正方法 | |
CN105260949B (zh) | 一种基于矩阵变换的配电网络短路计算方法 | |
JP2012139090A (ja) | 平衡配電系統のための電力潮流解析 | |
CN109802392B (zh) | 大规模配电网潮流计算方法及装置 | |
CN106130025B (zh) | 一种基于回路分析法计算三相配电网潮流的线性方法 | |
CN103066593A (zh) | 含多类型分布式电源的弱环配电网三相潮流计算方法 | |
CN107121604B (zh) | 一种不对称短路电流直流分量衰减时间常数获取方法 | |
CN102842907B (zh) | 基于道路矩阵的配电网三相解耦潮流计算方法 | |
CN104794531A (zh) | 基于网损等值负荷模型的改进直流最优潮流方法 | |
CN105184418B (zh) | 基于灵敏度和约束等值的互联电网最优潮流计算方法 | |
CN110880764B (zh) | 一种含逆变型分布式电源不平衡配电网的故障处理方法 | |
CN104252567A (zh) | 一种电力变压器三相漏磁导纳矩阵通用模型的建模方法 | |
CN104881566B (zh) | 一种用于电网不对称潮流计算的变压器三相建模方法 | |
CN106786540A (zh) | 一种短路电流计算方法及*** | |
Alsaadi et al. | An effective approach for distribution system power flow solution | |
Yang et al. | Three-phase power flow calculations using initial voltage estimation method for unbalanced distribution networks | |
CN110208634B (zh) | 一种复杂电力***不对称短路电流直流分量获取方法 | |
CN108595819B (zh) | 基于转换矩阵的变压器三相模型建立方法 | |
CN106410811A (zh) | 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法 | |
CN105098776A (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 |