CN105260949A - 一种基于矩阵变换的配电网络短路计算方法 - Google Patents
一种基于矩阵变换的配电网络短路计算方法 Download PDFInfo
- Publication number
- CN105260949A CN105260949A CN201510690724.8A CN201510690724A CN105260949A CN 105260949 A CN105260949 A CN 105260949A CN 201510690724 A CN201510690724 A CN 201510690724A CN 105260949 A CN105260949 A CN 105260949A
- Authority
- CN
- China
- Prior art keywords
- distribution network
- node
- phase
- admittance matrix
- load
- 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
本发明涉及一种基于矩阵变换的配电网短路计算方法,属于电力***配电网短路计算领域。该方法执行如下步骤:1)根据配电网三相不对称的特性,建立配电网三相***模型;2)利用配电网三相***模型,形成配电网络节点导纳矩阵;3)利用三相牛顿拉夫逊法对配电网络进行潮流计算;4)将集中负荷等效为恒阻抗并入配电网络模型,修改节点导纳矩阵;5)建立配电网络的节点电压方程;6)选择短路节点以及短路类型,以及平衡节点中性点是否接地来确定所对应的边界情况;7)对网络节点电压方程进行矩阵变换并求解以判断配电网络中的继电保护装置是否动作。该方法将负荷视作阻抗恒定的恒阻抗负荷,将负荷的等效阻抗并入节点导纳矩阵,极大地减少了方程中变量个数。
Description
技术领域
本发明涉及一种配电网短路计算方法,属于电力***配电网短路计算领域。
背景技术
短路故障是电力***中危害最严重的一种故障,在电力***设计和运行的许多工作中都须有短路计算的结果作为依据。以往人们比较注重发输电的管理和研究,现在己有一套成熟的输电网短路计算方法。
自20世纪80年代中期以来,随着国内外电力企业对配电网管理重视程度的提高,对配电网短路电流计算日益重视起来,但由于配电网与输电网在网络结构及***运行方式上有许多不同之处,如配电网通常以环型结构设计、开环方式运行,正常运行时又是辐射型树状结构:支路参数R/X的比值较大;分支线较多;单相负荷、主干馈线上的单相和两相“辐射状支线”以及不相等的三相负荷的存在,而使***通常是不对称的。
传统的短路电流计算一般采用对称分量法,其做法是:将***分解为三序网络并分别列出方程,结合故障点处的边界条件,计算出故障点处的电压、电流的各序分量,若要分析计算网络中任意处的电压、电流,必须先在各序网中求得该处电压、电流的各序分量,然后再合成三相电压、电流。对称分量法的优点是可以分别处理3个序网矩阵和方程,所以在过去的几十年里对称分量法得到了广泛的应用。然而,对称分量法在应用于配电网的短路电流计算时却遇到了困难,这是由配电网的上述特征导致网络参数不对称,对称分量法解耦失效,所以,相关计算必须采用abc全耦合模型。
发明内容
本发明涉及一种基于矩阵变换的配电网络短路计算方法,执行如下步骤:
1)根据所述配电网三相不对称的特性,建立配电网三相***模型;
在所述配电网模型中,所述配电网的上级变电站出口母线视为无穷大电源,并作为所述潮流计算中作为平衡节点,等值为电压幅值和相角为恒定已知量并假定三相电压对称;将所述配电网络的支路上所有负荷等效为集中负荷,所述集中负荷为恒功率PQ负荷;所述支路采用集中参数模型,三相之间的耦合电抗为Zli;
2)利用所述配电网三相***模型,形成配电网络节点导纳矩阵,所述配电网络节点导纳矩阵节点个数为3n×3n;
3)利用三相牛顿拉夫逊法对所述配电网络进行潮流计算;
利用步骤2)中的节点导纳矩阵,列出网络节点电压方程I=YU,进一步写出功率不平衡方程,并求解雅克比矩阵,利用三相牛顿拉夫逊法进行潮流计算,以求得所述配电网三相系通过模型中各节点电压情况以及负荷功率;
4)利用步骤3)中的计算结果,将所述集中负荷等效为恒阻抗并入所述配电网络模型,修改所述节点导纳矩阵;
5)利用步骤4)中修改后的节点导纳矩阵建立所述配电网络的节点电压方程;
6)选择短路节点以及短路类型,以及所述平衡节点中性点是否接地来确定所对应的边界情况;
所述短路类型包括单相接地短路、两相短路、两相接地短路和三相短路;
7)对所述网络节点电压方程进行矩阵变换并求解,求得所述配电网络中各节点电压,将算出的结果与所述配电网络的整定值比较,以此判断所述配电网络中的继电保护装置是否动作。
上述技术方案的改进是,步骤2)中通过如下方法形成所述配电网络导纳节点矩阵,将导纳矩阵扩展到三相的形式;将单相网络节点导纳矩阵中的各个元素用3×3的矩阵替换,获得新的三相节点导纳矩阵;
采用公式和Yik=-yik(k≠i)计算节点导纳矩阵;yi0为节点i对地导纳矩阵yij为节点i与节点j之间的支路导纳矩阵,即yik为节点i与节点k之间的支路导纳矩阵,即若两节点之间没有支路直接关联,则支路导纳矩阵为0。
上述技术方案的改进是,所述步骤7)中对节点电压方程进行矩阵变换并求解,求得网络各节点电压的方法包括:对于节点电压方程I=YU,由步骤6)所确定的边界条件,可获取短路计算结果。
本发明提供的一种基于矩阵变换的配电网短路计算方法的有益效果包括:
1、本发明提供的一种基于矩阵变换的配电网短路计算方法,考虑了配电网三相不对称的特性,分别对配电网平衡节点接地和不接地两种情况作了分析,总结不同情况的边界条件。
2、在短路情况下,将负荷等效为恒阻抗并入网络,修改节点导纳矩阵。结合边界条件,通过矩阵变换将未知量移动到方程一边,进而通过简单的矩阵计算求解网络方程,简化了求解步骤,缩短了计算时间。
3、将负荷视作阻抗恒定的恒阻抗负荷,将负荷的等效阻抗并入节点导纳矩阵,从而整个网络除了平衡节点外,注入电流均为0,极大地减少了方程中变量个数;再利用短路时的边界条件,将节点导纳矩阵处理成可以直接求解的方程,从而减小求解难度,缩短计算时间。
附图说明
图1为本发明一种基于矩阵变换的配电网短路计算方法实施例的流程图。
图2为本发明一种基于矩阵变换的配电网短路计算方法实施例的配电网三相***测试案例模型图。
图3为本发明一种基于矩阵变换的配电网短路计算方法实施例的单相接地短路情况。
图4为本发明一种基于矩阵变换的配电网短路计算方法实施例的两相短路情况。
图5为本发明一种基于矩阵变换的配电网短路计算方法实施例的两相接地短路情况。
图6为本发明一种基于矩阵变换的配电网短路计算方法实施例的三相短路情况。
具体实施方式
实施例
本实施例的一种基于矩阵变换的配电网络短路计算方法,如图1所示,执行如下步骤:
1)根据配电网三相不对称的特性,建立配电网三相***模型;
在配电网模型中,配电网的上级变电站出口母线视为无穷大电源,并作为潮流计算中作为平衡节点,等值为电压幅值和相角为恒定已知量并假定三相电压对称;将配电网络的支路上所有负荷等效为集中负荷,集中负荷为恒功率PQ负荷;支路采用集中参数模型,三相之间的耦合电抗为Zli;
2)利用配电网三相***模型,形成配电网络节点导纳矩阵,配电网络节点导纳矩阵节点个数为3n×3n;
3)利用三相牛顿拉夫逊法对配电网络进行潮流计算;
利用步骤2)中的节点导纳矩阵,列出网络节点电压方程I=YU,进一步写出功率不平衡方程,并求解雅克比矩阵,利用三相牛顿拉夫逊法进行潮流计算,以求得配电网三相系通过模型中各节点电压情况以及负荷功率;
4)利用步骤3)中的计算结果,将集中负荷等效为恒阻抗并入配电网络模型,修改节点导纳矩阵;
5)利用步骤4)中修改后的节点导纳矩阵建立配电网络的节点电压方程;
6)选择短路节点以及短路类型,以及平衡节点中性点是否接地来确定所对应的边界情况;
短路类型包括单相接地短路、两相短路、两相接地短路和三相短路;
7)对网络节点电压方程进行矩阵变换并求解,求得配电网络中各节点电压,将算出的结果与配电网络的整定值比较,以此判断配电网络中的继电保护装置是否动作。
本实施例的步骤2)中通过如下方法形成配电网络导纳节点矩阵,将导纳矩阵扩展到三相的形式;将单相网络节点导纳矩阵中的各个元素用3×3的矩阵替换,获得新的三相节点导纳矩阵;
采用公式和Yik=-yik(k≠i)计算节点导纳矩阵;yi0为节点i对地导纳矩阵yij为节点i与节点j之间的支路导纳矩阵,即yik为节点i与节点k之间的支路导纳矩阵,即若两节点之间没有支路直接关联,则支路导纳矩阵为0。
本实施例的步骤7)中对节点电压方程进行矩阵变换并求解,求得网络各节点电压的方法包括:对于节点电压方程I=YU,由步骤6)所确定的边界条件,可获取短路计算结果。
本发明提供的一种基于矩阵变换的配电网短路计算方法的实施例一中,优选的,步骤S1中建立配电网三相***模型包括:
将所述配电网的上级变电站出口母线视为无穷大电源,作为所述潮流计算中作为平衡节点,等值为电压幅值和相角为恒定已知量并假定三相电压对称,考虑中性点是否接地;将所述支路上所有负荷等效为支路末节点集中负荷,模型为可考虑不对称的恒功率PQ负荷;所述支路采用集中参数模型,考虑三相之间的耦合电抗,即
步骤S2中形成配电网络节点导纳矩阵包括:
参考单相网络节点导纳形成方法,将导纳矩阵扩展到三相的形式,即考虑abc三相模型;将单相网络节点导纳矩阵中的各个元素用3×3的矩阵替换,获得新的三相节点导纳矩阵。
步骤S3中利用三相牛顿拉夫逊法进行配电网潮流计算的方法包括:
利用步骤S2中的节点导纳矩阵,列出网络节点电压方程I=YU,进一步写出功率不平衡方程,求解雅克比矩阵,利用三相牛顿拉夫逊法进行潮流计算,求得各节点电压情况以及注入功率。
步骤S4中利用潮流计算结果,将网络负荷等效为恒阻抗并入网络,修改节点导纳矩阵包括:
由步骤S3得到网络各节点电压及注入功率,利用电压与功率间的关系求得网络负荷的等效阻抗,将等效阻抗并入网络节点导纳矩阵,并将负荷节点的注入电流视为0。
步骤S5中利用节点导纳矩阵建立网络节点电压方程的方法,即写出修改后的I=YU方程。
步骤S6中选择短路节点以及短路类型,列写对应的边界条件的方法包括:
短路类型有单相接地短路、两相短路、两相接地短路和三相短路,在考虑边界条件的情况下还需要考虑平衡节点中性点接地情况,中性点接地情况会影响接地短路情况下的对地电流,例如不接地***在单相接地情况下不存在对地电流,所以需要结合中性点接地情况和短路类型列写边界条件,即推出短路点各相电流、电压的值或关系式。
步骤S7中对节点电压方程进行矩阵变换并求解,求得网络各节点电压的方法包括:
对于节点电压方程I=YU,其中Y是已知的,而I、U中都含有已知量和未知量,所以不能直接求解;步骤S4将负荷等效为恒阻抗后。除平衡节点,各节点注入电流为0,即为已知量,电压中大部分为未知量;结合边界条件,方程可转换为左边为已知量,右边为Y乘以未知量,再对Y求逆直接求得方程中的未知量,从而整理获取短路计算结果。
如图2所示为本实施例的测试案例模型图。由图2可知,该配电网三相简化***网络即为纯辐射网,看作一棵树,源节点编号为0,其他节点按树的广度进行编号,各支路编号取支路末节点的节点编号,包括L1~L32。
按节点编号顺序存储配电网各节点三相负荷参数,如节点i应当存储有功功率Pia、Pib、Pic和无功功率Qia、Qib、Qic;按支路编号顺序存储配电网各支路三相阻抗参数,如支路Li应当存储阻抗参数为各相阻抗含实部电阻和虚部电抗。
对各个节点的三相负荷参数和各个支路的三相阻抗参数进行标幺化,本实施例中取UB=12.67kV,SB=10MW,对上一步中存储的各个节点的三相负荷参数和各个支路的三相阻抗参数进行标幺化,考虑节点的负荷功率是单相的,其基准值应当取阻抗基准值取
利用标幺化后的阻抗参数,建立节点导纳矩阵Y。建立网络功率方程,并进一步推导出功率不平衡方程,利用牛顿拉夫逊法计算网络潮流,得到各节点三相电压的幅值和相角。
短路情况下电压不能维持在正常水平,所以负荷并不能按潮流计算时的方法处理,大量文献中将负荷处理成恒定阻抗,便于短路计算。利用潮流计算结果以及负荷本身的参数,采用计算负荷的等值阻抗,其中U为潮流计算得到的电压,S为负荷参数P+jQ。
将得到的阻抗值并入节点导纳矩阵,即此时除平衡节点外各节点将不存在负荷,即各节点的注入电流均为0。新的网络方程为:
其中,节点导纳矩阵为恒阻抗并入网络后的新的节点导纳矩阵;左侧只有平衡节点注入电流I1a、I1b、I1c是未知量,其余均为0;右侧电压矩阵中除了平衡节点电压U1a、U1b、U1c是已知量,其余均为未知量。
发生短路故障时,节点导纳矩阵不变,只改变网络边界条件,只要方程个数等于未知量个数,便可通过矩阵变换进行求解。
1)单相接地短路(假设为第k-1号节点a相短路),如图3所示,
①在平衡节点中性点直接接地的情况下,Uka=0,Ika未知,利用矩阵变换得到易求解的形式如式(2)。
该方程左侧为已知量,右侧为一个已知方阵乘以未知量,两边同乘已知方阵的逆,便可直接求得未知量,从而得到各节点电压值以及注入电流值。进一步利用节点电压和线路参数,可以求得各线路流过的电流值。
②在平衡节点中性点不接地的情况下,平衡节点的电压变为对称电压叠加一个偏移量Uref,此外,短路对地电流因为不存在对地回路而为0。矩阵变换后的形式为式(3)。
2)两相短路(假设为第k-1号节点bc相短路),如图4所示,
因为没有对地支路,平衡节点中性点是否接地不影响短路计算。此时,边界条件为在节点电压方程中用-Ikb代替Ikc,用Ukb代替Ukc,再进行矩阵变换,得到方程组(4)。
3)两相接地短路(假设为第k-1号节点bc相短路),如图5所示,
①在平衡节点中性点直接接地的情况下,Ukb=Ukc=0,Ikb和Ikc未知,经过矩阵变换将Ikb和Ikc移到方程右边得方程(5)。
②在平衡节点中性点不接地的情况下,平衡节点的电压变为对称电压叠加一个偏移量Uref,此外,对地电流满足Ikb+Ikc=0。用-Ikb代替Ikc,利用矩阵变换将Uref和Ikb移到方程右侧得方程(6)。
4)三相短路(假设为第k-1号节点三相短路),如图6所示,
边界条件为用Ukc代替Uka和Ukb,用-(Ika+Ikb)代替Ikc,利用矩阵变化将Ika、Ikb和Ukc移动到方程右侧得方程(7)。
对于以上方程,都可以直接求解,最终解得各节点电压值,利用线路电流和两端电压的关系,结合线路阻抗,利用式(8)计算求得线路电流。
其中I12流过线路的电流,为3×1的向量;U1、U2分别为首末节点电压,为3×1的向量;Z12为线路阻抗,为3×1的矩阵。
案例采用IEEE33节点标准配电网测试***,每一条母线代表abc三相三个节点。结合案例仿真,利用Matlab编程处理,对表一中的六种情况进行仿真。
表1短路情况测试案例
1 | 中性点不接地,9号节点a相接地短路 |
2 | 中性点直接接地,9号节点a相接地短路 |
3 | 中性点不接地,9号节点bc两相短路 |
4 | 中性点不接地,9号节点bc两相接地短路 |
5 | 中性点直接接地,9号节点bc两相接地短路 |
6 | 中性点不接地,9号节点三相短路 |
经过仿真计算,该短路计算方法完成了各种情况的计算,并与输电网同种类型下的短路结果基本一致,表明了该算法的可行性。此处对单相接地短路两种情况的结果加以展示。
表2为平衡节点中性点直接接地***利用该短路算法求得的各节点电压幅值和相角,其中电压为标幺值,相角为角度值。
表2中性点直接接地
节点号 | a相幅值 | b相幅值 | c相幅值 | a相相角 | b相相角 | c相相角 |
0 | 1.0000 | 1.0000 | 1.0000 | 0.0000 | -120.0000 | 120.0000 |
1 | 0.9652 | 0.9992 | 0.9994 | 0.4511 | -120.0125 | 120.0257 |
2 | 0.8602 | 0.9961 | 0.9967 | 1.5290 | -120.0460 | 120.0989 |
3 | 0.8078 | 0.9951 | 0.9959 | 2.0709 | -120.0674 | 120.1302 |
4 | 0.7534 | 0.9941 | 0.9952 | 2.7138 | -120.0901 | 120.1627 |
5 | 0.6146 | 0.9921 | 0.9933 | 1.7321 | -120.1777 | 120.2205 |
6 | 0.5474 | 0.9925 | 0.9928 | -3.9615 | -120.2241 | 120.2570 |
7 | 0.4474 | 0.9923 | 0.9935 | -1.3730 | -120.2789 | 120.3037 |
8 | 0.2287 | 0.9921 | 0.9938 | -2.2162 | -120.3754 | 120.4291 |
9 | 0.0000 | 0.9929 | 0.9954 | 0.0000 | -120.4846 | 120.5665 |
10 | 0.0334 | 0.9929 | 0.9950 | -23.3281 | -120.4649 | 120.5472 |
11 | 0.0970 | 0.9928 | 0.9942 | -23.3158 | -120.4273 | 120.5089 |
12 | 0.1767 | 0.9913 | 0.9926 | -11.1038 | -120.3954 | 120.4548 |
13 | 0.2131 | 0.9906 | 0.9922 | -5.8793 | -120.3882 | 120.4220 |
14 | 0.2482 | 0.9903 | 0.9919 | -4.1164 | -120.3743 | 120.4012 |
15 | 0.2957 | 0.9895 | 0.9909 | -4.0490 | -120.3434 | 120.3753 |
16 | 0.3994 | 0.9876 | 0.9893 | 1.0964 | -120.2966 | 120.2882 |
17 | 0.4475 | 0.9870 | 0.9885 | 0.7729 | -120.2619 | 120.2650 |
18 | 0.9304 | 0.9988 | 0.9991 | 0.4956 | -120.0349 | 120.0487 |
19 | 0.6167 | 0.9957 | 0.9970 | 2.3129 | -120.2260 | 120.2702 |
20 | 0.5493 | 0.9949 | 0.9963 | -1.2579 | -120.2885 | 120.3367 |
21 | 0.4151 | 0.9948 | 0.9960 | -4.5767 | -120.3812 | 120.4541 |
22 | 0.8259 | 0.9944 | 0.9952 | 1.7659 | -120.0626 | 120.1188 |
23 | 0.7538 | 0.9912 | 0.9922 | 1.9569 | -120.1090 | 120.1496 |
24 | 0.6830 | 0.9891 | 0.9904 | 2.2287 | -120.1420 | 120.1971 |
25 | 0.6165 | 0.9917 | 0.9931 | 1.7133 | -120.1645 | 120.2208 |
26 | 0.6192 | 0.9913 | 0.9926 | 1.6875 | -120.1573 | 120.2233 |
27 | 0.6311 | 0.9894 | 0.9908 | 1.8512 | -120.1503 | 120.2078 |
28 | 0.6402 | 0.9882 | 0.9895 | 1.9691 | -120.1412 | 120.1997 |
29 | 0.6111 | 0.9872 | 0.9888 | 2.7102 | -120.1333 | 120.2402 |
30 | 0.5395 | 0.9866 | 0.9882 | 2.4287 | -120.2088 | 120.2506 |
31 | 0.5150 | 0.9866 | 0.9882 | 2.1043 | -120.2319 | 120.2573 |
32 | 0.4835 | 0.9869 | 0.9884 | 1.0896 | -120.2529 | 120.2751 |
表3为平衡节点中性点不接地***求得的各节点电压幅值和相角,其中,其中电压为标幺值,相角为角度值,求得的平衡节点中性点电压为-1.004+j0.000。
表3中性点不接地
节点号 | a相幅值 | b相幅值 | c相幅值 | a相相角 | b相相角 | c相相角 |
0 | 0.0004 | 1.7324 | 1.7324 | 179.7637 | -150.0062 | 150.0061 |
1 | 0.0003 | 1.7307 | 1.7307 | -179.0163 | -150.0003 | 150.0104 |
2 | 0.0002 | 1.7247 | 1.7245 | -172.6376 | -149.9835 | 150.0215 |
3 | 0.0002 | 1.7228 | 1.7224 | -169.3344 | -149.9814 | 150.0233 |
4 | 0.0001 | 1.7208 | 1.7205 | -165.0128 | -149.9792 | 150.0243 |
5 | 0.0001 | 1.7160 | 1.7159 | -153.9864 | -150.0125 | 149.9893 |
6 | 0.0000 | 1.7154 | 1.7154 | -166.9129 | -150.0441 | 149.9595 |
7 | 0.0000 | 1.7148 | 1.7148 | -171.0603 | -150.0494 | 149.9527 |
8 | 0.0000 | 1.7129 | 1.7130 | -37.4470 | -150.0563 | 149.9433 |
9 | 0.0000 | 1.7126 | 1.7129 | 0.0000 | -150.0641 | 149.9371 |
10 | 0.0000 | 1.7126 | 1.7129 | 143.5852 | -150.0649 | 149.9362 |
11 | 0.0000 | 1.7127 | 1.7130 | 147.4356 | -150.0673 | 149.9340 |
12 | 0.0000 | 1.7107 | 1.7110 | -5.4188 | -150.0651 | 149.9323 |
13 | 0.0000 | 1.7100 | 1.7102 | -4.1269 | -150.0691 | 149.9263 |
14 | 0.0000 | 1.7098 | 1.7100 | -7.6161 | -150.0672 | 149.9281 |
15 | 0.0001 | 1.7088 | 1.7089 | -11.8073 | -150.0595 | 149.9339 |
16 | 0.0001 | 1.7066 | 1.7068 | -8.7839 | -150.0571 | 149.9327 |
17 | 0.0001 | 1.7060 | 1.7062 | -11.5713 | -150.0455 | 149.9426 |
18 | 0.0003 | 1.7297 | 1.7297 | -178.9951 | -150.0044 | 150.0056 |
19 | 0.0001 | 1.7214 | 1.7214 | -174.1691 | -150.0259 | 149.9779 |
20 | 0.0001 | 1.7191 | 1.7191 | -175.2966 | -150.0407 | 149.9619 |
21 | 0.0000 | 1.7169 | 1.7171 | 166.5125 | -150.0594 | 149.9429 |
22 | 0.0002 | 1.7215 | 1.7213 | -167.0823 | -149.9832 | 150.0170 |
23 | 0.0000 | 1.7153 | 1.7151 | -114.8659 | -149.9956 | 149.9941 |
24 | 0.0001 | 1.7110 | 1.7110 | -33.1200 | -149.9947 | 149.9888 |
25 | 0.0001 | 1.7152 | 1.7154 | -145.0894 | -150.0073 | 149.9848 |
26 | 0.0000 | 1.7144 | 1.7147 | -129.5686 | -150.0015 | 149.9899 |
27 | 0.0001 | 1.7115 | 1.7116 | -52.4682 | -149.9997 | 149.9883 |
28 | 0.0001 | 1.7095 | 1.7096 | -33.4667 | -149.9946 | 149.9907 |
29 | 0.0001 | 1.7076 | 1.7077 | -31.7012 | -149.9709 | 150.0139 |
30 | 0.0001 | 1.7058 | 1.7060 | -16.4554 | -150.0144 | 149.9714 |
31 | 0.0001 | 1.7055 | 1.7057 | -13.1560 | -150.0272 | 149.9592 |
32 | 0.0001 | 1.7057 | 1.7059 | -11.5015 | -150.0360 | 149.9512 |
从以上结果可以看出,中性点是否接地对网络短路计算结果的影响较大。中性点直接接地情况下发生单相短路,短路相电压从平衡节点到短路点逐渐降低,其余相变化较小;中性点不接地情况下发生单相短路,短路相电压全部降低到0,非短路相电压升高为原来的倍,中性点电压由原来的0变为-1.004+j0.000,尽管网络相电压值发生了变化,但是各相与平衡节点之间的电压差并没有多大变化,线电压基本不变,这就是为什么不接地***发生单相短路时还能继续运行一段时间。
以上虽然根据附图对本发明的实施例进行了详细说明,但不仅限于此具体实施方式,本领域的技术人员根据此具体技术方案进行的各种等同、变形处理,也在本发明的保护范围之内。
Claims (3)
1.一种基于矩阵变换的配电网络短路计算方法,其特征在于执行如下步骤:
1)根据所述配电网三相不对称的特性,建立配电网三相***模型;
在所述配电网模型中,所述配电网的上级变电站出口母线视为无穷大电源,并作为所述潮流计算中作为平衡节点,等值为电压幅值和相角为恒定已知量并假定三相电压对称;将所述配电网络的支路上所有负荷等效为集中负荷,所述集中负荷为恒功率PQ负荷;所述支路采用集中参数模型,三相之间的耦合电抗为Zli;
2)利用所述配电网三相***模型,形成配电网络节点导纳矩阵,所述配电网络节点导纳矩阵节点个数为3n×3n;
3)利用三相牛顿拉夫逊法对所述配电网络进行潮流计算;
利用步骤2)中的节点导纳矩阵,列出网络节点电压方程I=YU,进一步写出功率不平衡方程,并求解雅克比矩阵,利用三相牛顿拉夫逊法进行潮流计算,以求得所述配电网三相系通过模型中各节点电压情况以及负荷功率;
4)利用步骤3)中的计算结果,将所述集中负荷等效为恒阻抗并入所述配电网络模型,修改所述节点导纳矩阵;
5)利用步骤4)中修改后的节点导纳矩阵建立所述配电网络的节点电压方程;
6)选择短路节点以及短路类型,以及所述平衡节点中性点是否接地来确定所对应的边界情况;
所述短路类型包括单相接地短路、两相短路、两相接地短路和三相短路;
7)对所述网络节点电压方程进行矩阵变换并求解,求得所述配电网络中各节点电压,将算出的结果与所述配电网络的整定值比较,以此判断所述配电网络中的继电保护装置是否动作。
2.如权利要求1所述的基于矩阵变换的配电网络短路计算方法,其特征在于步骤2)中通过如下方法形成所述配电网络导纳节点矩阵,将导纳矩阵扩展到三相的形式;将单相网络节点导纳矩阵中的各个元素用3×3的矩阵替换,获得新的三相节点导纳矩阵;
采用公式和Yik=-yik(k≠i)计算节点导纳矩阵;yi0为节点i对地导纳矩阵 yij为节点i与节点j之间的支路导纳矩阵,即yik为节点i与节点k之间的支路导纳矩阵,即若两节点之间没有支路直接关联,则支路导纳矩阵为0。
3.如权利要求1所述的方法,其特征在于,所述步骤7)中对节点电压方程进行矩阵变换并求解,求得网络各节点电压的方法包括:对于节点电压方程I=YU,由步骤6)所确定的边界条件,可获取短路计算结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510690724.8A CN105260949B (zh) | 2015-10-22 | 2015-10-22 | 一种基于矩阵变换的配电网络短路计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510690724.8A CN105260949B (zh) | 2015-10-22 | 2015-10-22 | 一种基于矩阵变换的配电网络短路计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105260949A true CN105260949A (zh) | 2016-01-20 |
CN105260949B CN105260949B (zh) | 2019-07-09 |
Family
ID=55100624
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510690724.8A Active CN105260949B (zh) | 2015-10-22 | 2015-10-22 | 一种基于矩阵变换的配电网络短路计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105260949B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105977962A (zh) * | 2016-01-29 | 2016-09-28 | 中国电力科学研究院 | 一种基于改进节点法的输电网-配电网联合故障分析方法 |
CN107506553A (zh) * | 2017-08-28 | 2017-12-22 | 国家电网公司 | 适用于新能源电源接入不平衡配电网的短路电流计算方法 |
CN109254187A (zh) * | 2017-07-14 | 2019-01-22 | 中国电力工程顾问集团华东电力设计院有限公司 | 一种确定变电站***短路电流的***及方法 |
CN110490340A (zh) * | 2019-08-30 | 2019-11-22 | 深圳供电局有限公司 | 配网接地故障计算方法 |
CN111900738A (zh) * | 2020-08-14 | 2020-11-06 | 广东电网有限责任公司广州供电局 | 基于补偿算法的三相不平衡潮流计算方法 |
CN112865110A (zh) * | 2021-01-27 | 2021-05-28 | 国网安徽省电力有限公司电力科学研究院 | 一种配电线路与负荷间的阻抗式电气距离计算方法 |
CN113447803A (zh) * | 2021-06-08 | 2021-09-28 | 国网上海市电力公司 | 断路器开断能力校核的短路电流计算电压系数取值方法 |
CN113884811A (zh) * | 2021-10-08 | 2022-01-04 | 邓朝尹 | 一种基于直算法的配网线路短路故障定位方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103353918A (zh) * | 2013-04-24 | 2013-10-16 | 国家电网公司 | 一种配电网三相牛顿拉夫逊法潮流计算初值选取方法 |
CN103529275A (zh) * | 2013-10-25 | 2014-01-22 | 国家电网公司 | 地区电网短路电流分析方法和装置 |
-
2015
- 2015-10-22 CN CN201510690724.8A patent/CN105260949B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103353918A (zh) * | 2013-04-24 | 2013-10-16 | 国家电网公司 | 一种配电网三相牛顿拉夫逊法潮流计算初值选取方法 |
CN103529275A (zh) * | 2013-10-25 | 2014-01-22 | 国家电网公司 | 地区电网短路电流分析方法和装置 |
Non-Patent Citations (5)
Title |
---|
CÉSAR AUGUSTO PEÑUELA MENESES 等: ""Probabilistic Algorithms for Power Load Flow and Short-Circuit Analysis in Distribution Networks with Dispersed Generation"", 《JOURNAL OF CONTROL, AUTOMATION AND ELECTRICAL SYSTEMS》 * |
CHANGQING KANG: ""A Power Flow Solution by Newton-Raphson and Time Domain Simulation"", 《2010 INTERNATIONAL CONFERENCE ON BROADCAST TECHNOLOGY AND MULTIMEDIA COMMUNICATION》 * |
张笑: ""配电网短路故障计算方法研究"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
盛爽: ""基于相分量法的配电***短路故障分析研究"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
车仁飞 等: ""基于叠加原理的配电网短路电流计算"", 《电力***自动化》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105977962B (zh) * | 2016-01-29 | 2021-08-20 | 中国电力科学研究院 | 一种基于改进节点法的输电网-配电网联合故障分析方法 |
CN105977962A (zh) * | 2016-01-29 | 2016-09-28 | 中国电力科学研究院 | 一种基于改进节点法的输电网-配电网联合故障分析方法 |
CN109254187A (zh) * | 2017-07-14 | 2019-01-22 | 中国电力工程顾问集团华东电力设计院有限公司 | 一种确定变电站***短路电流的***及方法 |
CN107506553A (zh) * | 2017-08-28 | 2017-12-22 | 国家电网公司 | 适用于新能源电源接入不平衡配电网的短路电流计算方法 |
CN110490340A (zh) * | 2019-08-30 | 2019-11-22 | 深圳供电局有限公司 | 配网接地故障计算方法 |
CN110490340B (zh) * | 2019-08-30 | 2022-04-26 | 深圳供电局有限公司 | 配网接地故障计算方法 |
CN111900738A (zh) * | 2020-08-14 | 2020-11-06 | 广东电网有限责任公司广州供电局 | 基于补偿算法的三相不平衡潮流计算方法 |
CN111900738B (zh) * | 2020-08-14 | 2022-02-11 | 广东电网有限责任公司广州供电局 | 基于补偿算法的三相不平衡潮流计算方法 |
CN112865110A (zh) * | 2021-01-27 | 2021-05-28 | 国网安徽省电力有限公司电力科学研究院 | 一种配电线路与负荷间的阻抗式电气距离计算方法 |
CN113447803A (zh) * | 2021-06-08 | 2021-09-28 | 国网上海市电力公司 | 断路器开断能力校核的短路电流计算电压系数取值方法 |
CN113447803B (zh) * | 2021-06-08 | 2024-05-28 | 国网上海市电力公司 | 断路器开断能力校核的短路电流计算电压系数取值方法 |
CN113884811A (zh) * | 2021-10-08 | 2022-01-04 | 邓朝尹 | 一种基于直算法的配网线路短路故障定位方法 |
CN113884811B (zh) * | 2021-10-08 | 2024-04-19 | 邓朝尹 | 一种基于直算法的配网线路短路故障定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105260949B (zh) | 2019-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105260949A (zh) | 一种基于矩阵变换的配电网络短路计算方法 | |
Kamh et al. | Unbalanced model and power-flow analysis of microgrids and active distribution systems | |
CN104578159B (zh) | 一种含分布式电源的配电网三相潮流修正方法 | |
Li et al. | A supply restoration method of distribution system based on soft open point | |
CN102104252B (zh) | 一种适用于电磁暂态分析的电力***动态等值方法 | |
CN103353918A (zh) | 一种配电网三相牛顿拉夫逊法潮流计算初值选取方法 | |
CN109802392B (zh) | 大规模配电网潮流计算方法及装置 | |
CN103605829A (zh) | 对交直流混联电网进行电磁暂态仿真的等值建模方法 | |
CN103795057A (zh) | 配电网潮流计算中基于搜索的配电网拓扑编号生成方法 | |
CN105162099A (zh) | 一种确定分布式发电接入电网不对称短路电流的运算曲面法 | |
CN106532710A (zh) | 计及电压稳定约束的微电网潮流优化方法 | |
CN104113061A (zh) | 一种含分布式电源的配电网三相潮流计算方法 | |
CN105162144A (zh) | 一种配电网电压无功优化的线性逼近求解方法 | |
CN105958429A (zh) | 一种降低交直流***中工频故障分量距离保护误动率的单相重合时序整定方法 | |
CN107196307B (zh) | 一种变压器故障后电网有功潮流快速估算的方法 | |
CN103560494B (zh) | 配电网保护整定的短路电流获取方法 | |
CN104917197A (zh) | 一种并行计算主动配电网三相不平衡潮流的方法 | |
CN102938010A (zh) | 平行线路不平衡度分析方法 | |
CN107104463A (zh) | 换流站在黑启动过程中考虑无功的有功投入量优化方法 | |
Ustun et al. | Modeling and simulation of a microgrid protection system with central protection unit | |
Zhang et al. | A three-phase power flow algorithm for ungrounded network based on constraints of zero-sequence components | |
CN103390890B (zh) | 基于电流分布系数的配电网潮流分析方法 | |
CN110048427A (zh) | 一种多端柔性直流配电网潮流保护定值计算方法 | |
Sun et al. | Application of GUPFC in China's Sichuan power grid-modeling, control strategy and case study | |
CN107436995A (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 |