CN104300490A - 一种继电保护在线断点求取方法 - Google Patents

Abstract

本发明公开了一种继电保护在线断点求取方法。该方法采用连通性分析方法确定了全网分区和最小断点集的维数，在加快断点集的求取速度的同时，保证了最终断点集的正确性；综合考虑了分区内部和分区间的支路变化，根据拓扑变化类型的不同，实现了网络分区的自动更正和拓扑变化时断点集维数的更新；基于追踪技术求取在线最小断点集，无需重新计算拓扑改变后的全网网络结构，减少了大量的重复工作，提高了计算速度。该方法能有效克服传统算法在电网拓扑动态变化时求取最小断点集计算量大、耗时长的缺点，获得满足在线校核实时性要求的在线最小断点集。

Description

一种继电保护在线断点求取方法

技术领域

本发明属于交流输电系统的继电保护技术领域，更具体地，涉及一种继电保护在线断点求取方法。

背景技术

在整定计算过程中，对于处于同一环路、同段配合的保护，会出现相互等待、无法继续整定的情况。为了解决这一问题，必须在整定之前确定一组能断开所有回路的保护作为整定起始点，这些整定起始点就是整定计算中的断点集。在线校核计算是整定计算的逆过程，同样存在断点集求取的问题。而且，与传统整定计算相比，在线校核计算最大的特点之一就是其实时性，因此，在线校核对断点的计算速度提出了更高的要求。

目前，国内外对最小断点集求取的研究主要可以分为三类：基于图论的断点求取方法、基于配合关系的断点求取方法和基于人工智能的断点求取方法。这些方法逐步完善了最小断点集的求取体系，在计算速度上也有了大幅度提高。但是，将这些方法应用于在线校核计算仍存在以下两方面的不足：一是大部分算法在提高断点求取速度时，牺牲了最小断点集的准确度，得到的都不是最优断点集，而是近似最小断点集；二是其考虑静态电网拓扑下断点的求取方法，当电网拓扑发生变化时，需要对全网的断点进行重新计算，计算量大，重复操作多，无法满足在线校核实时性的要求。

发明内容

针对现有技术的以上缺陷或改进需求，本发明提供了一种继电保护在线断点求取方法，其目的在于克服传统算法在电网拓扑动态变化时求取最小断点集计算量大、耗时长的缺点，获得满足在线校核实时性要求的在线最小断点集。

为实现上述目的，本发明提供了一种继电保护在线断点求取方法，其特征在于，包括如下步骤：

(1)初始化网络拓扑信息，形成全网主后备保护配合矩阵P和节点关联矩阵V；

(2)搜寻全网割节点，对全网进行分区；

(3)确定全网最小断点集的维数；

(4)求取各分区的离线断点，合并各分区的断点得到原始全网最小断点集BPS_old；

(5)判断网络拓扑是否发生变化，是则跳至步骤(7)；否则顺序执行步骤(6)；

(6)等待一个时间周期，返回步骤(5)；

(7)在网络拓扑变化发生在分区内部时，不更新全网分区，对全网最小断点集的维数进行修正；在网络拓扑变化发生在两个分区之间时，将这两个分区合并后，对合并后的区域进行割节点搜索和重新分区，计算新分区的最小断点集的维数，再根据分区数目的改变对全网最小断点集的维数进行修正；

(8)更新全网主后备保护配合矩阵P，得到网络拓扑变化后的主后备保护配合矩阵P'，求取最终断点集。

优选地，所述步骤(1)中，主后备保护配合矩阵P的矩阵规模为n×n，n为网络中的保护总数，P中的第i行第j列元素P(i,j)为：

其中，i＝1,2,…,n，j＝1,2,…,n；

节点关联矩阵V的矩阵规模为m×m，m为网络中的节点总数，V中的第i'行第j'列元素V(i',j')为：

其中，i'＝1,2,…,m，j'＝1,2,…,m。

优选地，所述步骤(2)进一步包括如下步骤：

(2-1)搜寻终端割节点，进一步包括如下步骤：

(2-1-1)分别计算节点关联矩阵V的各行之和R_sum(i')；

(2-1-2)获取R_sum(i')＝1时对应的i'值，

(2-1-3)搜索第i'行，得到V(i',j')＝1时对应的j'的值，删除节点关联矩阵V的第i'行和第i'列；

(2-1-4)判断节点关联矩阵V的第j'行之和R_sum(j')是否为1，是则将j'的值赋给i'，返回步骤(2-1-3)；否则判断节点j'是割节点，将节点j'加入割节点集G，顺序执行步骤(2-1-5)；

(2-1-5)判断节点关联矩阵V是否存在和为1的行，是则返回步骤(2-1-2)；否则过程结束；

(2-2)搜寻非终端割节点，进一步包括如下步骤：

(2-2-1)选择节点关联矩阵V中的第一个节点作为当前节点；

(2-2-2)将当前节点所在行中所有不为1的元素存入临时起始链表；

(2-2-3)遍历临时起始链表，判断剔除当前节点后，临时起始链表中元素的连通性，若该链表中的所有元素都连通，则当前节点不是割节点，跳至步骤(2-2-5)；若不连通，则当前节点是割节点，顺序执行步骤(2-2-4)；

(2-2-4)将当前节点存入割节点集G；

(2-2-5)判断是否遍历节点关联矩阵V中的所有节点，是则过程结束；否则选择下一节点作为当前节点，返回步骤(2-2-2)；

(2-3)断开割节点集G中的割节点，得到电网的自动分区。

优选地，所述步骤(3)中，对于不可分连通网络，全网最小断点集的维数N_BPS＝N_Z-N_V+2；对于可分连通网络，全网最小断点集的维数N_BPS＝N_Z-N_V+k+1；其中，N_Z为网络支路数，N_V为网络节点数，k为可分连通网络能分为独立分离的连通子图的个数。

优选地，所述步骤(8)进一步包括如下步骤：

(8-1)判断拓扑变化是否为新增支路，是则顺序执行步骤(8-2)；否则为删除支路，跳至步骤(8-3)；

(8-2)新增支路两端保护构成集合BPS₁，基于网络拓扑变化后的主后备保护配合矩阵P'，追踪被BPS₁覆盖但未被原始全网最小断点集BPS_old覆盖的回路集L_new，并求取回路集L_new的断点集作为新增断点集BPS_new，合并新增断点集BPS_new和原始全网最小断点集BPS_old，得到临时最小断点集BPS_temp，跳至步骤(8-4)；

(8-3)在删除支路为原始断点所在支路时，直接在原始全网最小断点集BPS_old中删除该断点，得到临时最小断点集BPS_temp；在删除支路为非断点支路时，将原始全网最小断点集BPS_old作为临时最小断点集BPS_temp，同时，选取BPS_temp中的任一断点作为新增断点集BPS_new，BPS_temp中其他断点构成集合BPS₂，追踪被BPS_new覆盖但未被BPS₂覆盖的回路集L_new；

(8-4)将原始全网最小断点集BPS_old中受拓扑变化影响的断点存入集合BPS_delete中；

(8-5)判断临时最小断点集BPS_temp的维数是否等于修正后的全网最小断点集的维数，是则将临时最小断点集BPS_temp作为最终断点集BPS_final，过程结束；否则顺序执行步骤(8-6)；

(8-6)剔除BPS_delete中的部分断点，得到新的集合BPS'_delete，剔除的断点构成集合BPS₃，追踪被BPS₃覆盖但未被BPS'_delete中断点覆盖的回路集L_delete，将L_delete与L_new合并，形成新的回路集L′，求取新的回路集L′的最小断点集BPS_new′，合并BPS_new′和BPS'_delete，得到新的临时最小断点集BPS'_temp，令BPS_temp＝BPS'_temp，返回步骤(8-5)。

优选地，所述步骤(8-2)、(8-3)和(8-6)中，按如下步骤追踪被A^*覆盖但未被B^*覆盖的回路集L：

(A)初始化访问标志位数组CF中的所有元素为0；其中，CF中元素的个数等于拓扑变化后的网络的保护数；

(B)令r＝1；

(C)以A^*中的保护start(r)为起始追踪保护，将其对应的保护编号值赋给变量x，其在CF中的标志位CF(x)置1，并将其添加到临时回路链表LOOP_temp；

(D)令y＝1；

(E)判断网络拓扑变化后的主后备保护配合矩阵P'中第x行的第y个元素P'(x,y)是否为1，且CF(y)是否为0，是则顺序执行步骤(F)；否则跳至步骤(H)；

(F)判断保护y是否为B^*中的断点，是则跳至步骤(H)，否则顺序执行步骤(G)；

(G)判断保护y是否为起始追踪保护start(r)，是则将临时回路链表LOOP_temp添加到回路链表集LOOP，顺序执行步骤(H)；否则将保护y添加到临时回路链表LOOP_temp，并将其在CF中的标志位CF(y)置1，将y赋值给x，返回步骤(D)；

(H)判断第x行的所有元素是否遍历完毕，未遍历完毕则y＝y+1，返回步骤(E)；遍历完毕则判断当前临时回路链表LOOP_temp的最后一个元素是否为start(r)，是则剔除该元素；

(I)判断A^*中的保护是否遍历完毕，未遍历完毕则r＝r+1，返回步骤(C)；遍历完毕则过程结束，由回路链表集LOOP得到被A^*覆盖但未被B^*覆盖的回路集L；

其中，对所述步骤(8-2)，A^*为BPS₁，B^*为BPS_old，L为L_new；对所述步骤(8-3)，A^*为BPS_new，B^*为BPS₂，L为L_new；对所述步骤(8-6)，A^*为BPS₃，B^*为BPS'_delete，L为L_delete。

总体而言，通过本发明所构思的以上技术方案与现有技术相比，采用连通性分析方法确定了全网分区和最小断点集的维数，在加快断点集的求取速度的同时，保证了最终断点集的正确性；综合考虑了分区内部和分区间的支路变化，根据拓扑变化类型的不同，实现了网络分区的自动更正和拓扑变化时断点集维数的更新；基于追踪技术求取在线最小断点集，无需重新计算拓扑改变后的全网网络结构，减少了大量的重复工作，提高了计算速度。

附图说明

图1是本发明实施例的继电保护在线断点求取方法流程图；

图2是典型的8节点系统图；

图3是增加支路后的典型的8节点系统图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

如图1所示，本发明实施例的继电保护在线断点求取方法包括如下步骤：

(1)初始化网络拓扑信息，形成全网主后备保护配合矩阵P和节点关联矩阵V。

主后备保护配合矩阵P表征保护之间的配合关系，其矩阵规模为n×n，其中，n为网络中的保护总数，主后备保护配合矩阵P中的第i行第j列元素P(i,j)为：

其中，i＝1,2,…,n，j＝1,2,…,n。由主后备保护配合矩阵P可知：P的每行非0元素之和表示该行保护需要依赖的主保护个数，每列非0元素之和表示依赖于该列的后备保护个数。

节点关联矩阵V表征了节点之间的关联关系，其矩阵规模为m×m，其中，m为网络中的节点总数，节点关联矩阵V中的第i'行第j'列元素V(i',j')为：

其中，i'＝1,2,…,m，j'＝1,2,…,m。

以图2所示系统为例，遍历系统中的所有保护得到主后备保护配合矩阵P为：

遍历系统中的所有节点得到节点关联矩阵V为：

$\left[\begin{array}{cccccccc}0& 1& 1& 1& 0& 0& 0& 0\\ 1& 0& 1& 0& 0& 0& 0& 0\\ 1& 1& 0& 0& 1& 0& 0& 1\\ 1& 0& 0& 0& 0& 1& 0& 1\\ 0& 0& 1& 0& 0& 0& 0& 1\\ 0& 0& 0& 1& 0& 0& 1& 0\\ 0& 0& 0& 0& 0& 1& 0& 0\\ 0& 0& 1& 1& 1& 0& 0& 0\end{array}\right].$

(2)搜寻全网割节点，对全网进行分区。

对任意给定连通网络，若断开网络中某条支路或某个节点(同时包括与该节点直接关联的支路)，该网络变成两个或多个独立分离的连通子图，则被断开的支路或节点称为割支路或割节点。

步骤(2)进一步包括如下步骤：

(2-1)搜寻终端割节点。进一步包括如下步骤：

(2-1-1)分别计算节点关联矩阵V的各行之和R_sum(i')；

(2-1-2)获取R_sum(i')＝1时对应的i'值；

(2-1-3)搜索第i'行，得到V(i',j')＝1时对应的j'的值，删除节点关联矩阵V的第i'行和第i'列；

(2-1-4)判断节点关联矩阵V的第j'行之和R_sum(j')是否为1，是则将j'的值赋给i'，返回步骤(2-1-3)；否则判断节点j'是割节点，将节点j'加入割节点集G，顺序执行步骤(2-1-5)；

(2-1-5)判断节点关联矩阵V是否存在和为1的行，是则返回步骤(2-1-2)；否则过程结束。

(2-2)搜寻非终端割节点。进一步包括如下步骤：

(2-2-1)选择节点关联矩阵V中的第一个节点作为当前节点；

(2-2-2)将当前节点所在行中所有不为1的元素存入临时起始链表；

(2-2-3)遍历临时起始链表，判断剔除当前节点后，临时起始链表中元素的连通性，若该链表中的所有元素都连通，则当前节点不是割节点，跳至步骤(2-2-5)；若不连通，则当前节点是割节点，顺序执行步骤(2-2-4)；

(2-2-4)将当前节点存入割节点集G；

(2-2-5)判断是否遍历节点关联矩阵V中的所有节点，是则过程结束；否则选择下一节点作为当前节点，返回步骤(2-2-2)。

(2-3)断开割节点集G中的割节点，得到电网的自动分区。

以图2所示系统为例，其终端割节点为节点④，无非终端割节点，全网所对应的连通图可以被分为2个独立分离的连通子图。

(3)确定全网最小断点集的维数。

为了确保在线断点求取的正确性，需要对最小断点集的维数进行确定。

对于不可分连通网络，全网最小断点集的维数N_BPS＝N_L+1，其中，N_L为网络连支数；对于可分连通网络，全网最小断点集的维数N_BPS＝N_L+k，其中，k为该连通网络能分为独立分离的连通子图的个数。

根据图论知识，网络连支数N_L与网络支路数N_Z、网络树枝数N_S和网络节点数N_V存在如下关系：N_L＝N_Z-N_S＝N_Z-(N_V-1)，因此，对于不可分连通网络，全网最小断点集的维数N_BPS＝N_Z-N_V+2；对于可分连通网络，全网最小断点集的维数N_BPS＝N_Z-N_V+k+1。对于给定电网，N_Z和N_V均为已知量，因此，求取最小断点集维数的问题转换为求取电网不可分连通子图数的问题。进一步将最小断点集维数的确定问题转化为寻找全网割节点的问题。

以图2所示系统为例，全网最小断点集的维数N_BPS为：

N_BPS＝N_Z-N_V+k+1＝11-8+2+1＝6

(4)求取各分区的离线断点，合并各分区的断点得到原始全网最小断点集BPS_old。

以图2所示系统为例，原始全网最小断点集为{6，7，8，9，17，19}。

(5)判断网络拓扑是否发生变化，是则跳至步骤(7)；否则顺序执行步骤(6)。

(6)等待一个时间周期，返回步骤(5)。

(7)在网络拓扑变化发生在分区内部时，不更新全网分区，根据连枝数N_L的变化对全网最小断点集的维数进行修正；在网络拓扑变化发生在两个分区之间时，将这两个分区合并后，按照步骤(2)所述的方法对合并后的区域进行割节点搜索和重新分区，计算新分区的最小断点集的维数，再根据分区数目的改变对全网最小断点集的维数进行修正。

(8)更新全网主后备保护配合矩阵P，得到网络拓扑变化后的主后备保护配合矩阵P'，求取最终断点集。进一步包括如下步骤：

(8-1)判断拓扑变化是否为新增支路，是则顺序执行步骤(8-2)；否则为删除支路，跳至步骤(8-3)；

(8-2)新增支路两端保护构成集合BPS₁，基于网络拓扑变化后的主后备保护配合矩阵P'，追踪被BPS₁覆盖但未被原始全网最小断点集BPS_old覆盖的回路集L_new，并求取回路集L_new的断点集作为新增断点集BPS_new，合并新增断点集BPS_new和原始全网最小断点集BPS_old，得到临时最小断点集BPS_temp，跳至步骤(8-4)；

(8-3)在删除支路为原始断点所在支路时，直接在原始全网最小断点集BPS_old中删除该断点，得到临时最小断点集BPS_temp；在删除支路为非断点支路时，将原始全网最小断点集BPS_old作为临时最小断点集BPS_temp，同时，选取BPS_temp中的任一断点作为新增断点集BPS_new，BPS_temp中其他断点构成集合BPS₂，追踪被BPS_new覆盖但未被BPS₂覆盖的回路集L_new；

(8-4)将原始全网最小断点集BPS_old中受拓扑变化影响的断点存入集合BPS_delete中；

(8-5)判断临时最小断点集BPS_temp的维数是否等于修正后的全网最小断点集的维数，是则将临时最小断点集BPS_temp作为最终断点集BPS_final，过程结束；否则顺序执行步骤(8-6)；

(8-6)剔除BPS_delete中的部分断点，得到新的集合BPS'_delete，剔除的断点构成集合BPS₃，追踪被BPS₃覆盖但未被BPS'_delete中断点覆盖的回路集L_delete，将L_delete与L_new合并，形成新的回路集L′，求取新的回路集L′的最小断点集BPS_new′，合并BPS_new′和BPS'_delete，得到新的临时最小断点集BPS'_temp，令BPS_temp＝BPS'_temp，返回步骤(8-5)。

其中，按如下步骤追踪被A^*覆盖但未被B^*覆盖的回路集L：

(A)初始化访问标志位数组CF中的所有元素为0，即未被访问过；其中，CF中元素的个数等于拓扑变化后的网络的保护数；

(B)令r＝1；

(C)以A^*中的保护start(r)为起始追踪保护，将其对应的保护编号值赋给变量x，其在CF中的标志位CF(x)置1，并将其添加到临时回路链表LOOP_temp；

(D)令y＝1；

(E)判断网络拓扑变化后的主后备保护配合矩阵P'中第x行的第y个元素P'(x,y)是否为1，且CF(y)是否为0，是则顺序执行步骤(F)；否则跳至步骤(H)；

(F)判断保护y是否为B^*中的断点，是则跳至步骤(H)，否则顺序执行步骤(G)；

(G)判断保护y是否为起始追踪保护start(r)，是则将临时回路链表LOOP_temp添加到回路链表集LOOP，顺序执行步骤(H)；否则将保护y添加到临时回路链表LOOP_temp，并将其在CF中的标志位CF(y)置1，将y赋值给x，返回步骤(D)；

(H)判断第x行的所有元素是否遍历完毕，未遍历完毕则y＝y+1，返回步骤(E)；遍历完毕则判断当前临时回路链表LOOP_temp的最后一个元素是否为start(r)，是则剔除该元素；

(I)判断A^*中的保护是否遍历完毕，未遍历完毕则r＝r+1，返回步骤(C)；遍历完毕则过程结束，由回路链表集LOOP得到被A^*覆盖但未被B^*覆盖的回路集L。

其中，对步骤(8-2)，A^*为BPS₁，B^*为BPS_old，L为L_new；对步骤(8-3)，A^*为BPS_new，B^*为BPS₂，L为L_new；对步骤(8-6)，A^*为BPS₃，B^*为BPS'_delete，L为L_delete。

以图2所示的系统为例，在节点⑤和节点⑦之间增加一条支路和两个保护，得到新的系统如图3所示。

网络拓扑变化后的主后备保护配合矩阵P'和节点关联矩阵V'分别为：

和

$\left[\begin{array}{cccccccc}0& 1& 1& 1& 0& 0& 0& 0\\ 1& 0& 1& 0& 0& 0& 0& 0\\ 1& 1& 0& 0& 1& 0& 0& 1\\ 1& 0& 0& 0& 0& 1& 0& 1\\ 0& 0& 1& 0& 0& 0& 1& 1\\ 0& 0& 0& 1& 0& 0& 1& 0\\ 0& 0& 0& 0& 1& 1& 0& 0\\ 0& 0& 1& 1& 1& 0& 0& 0\end{array}\right].$

由于新增的支路在两分区之间，将这两个分区合并，重新寻找割节点。最终找寻到的割节点数为0，将全网最小断点集的维数更新为：

N'_BPS＝N_Z-N_V+2＝12-8+2＝6

追踪被保护23和保护24覆盖，但未被原始断点集{6，7，8，9，17，19}中的断点覆盖的回路集，最终新增回路为24→22→11→16→18，新增断点为24。临时最小断点集为{6，7，8，9，17，19,24}，维数为7，大于6。剔除断点及寻找最优断点集的顺序如表1所示，最终临时断点集{6，7，8，9，21,22}满足维数要求，为最优断点集，将其作为最终断点集。

表1

依据上述在线断点求取过程，分别对任意节点间添加支路和删除支路的情况进行仿真，仿真结果如表2所示。

表2

分析上述仿真结果可知：(1)本发明的继电保护在线断点求取方法在各种拓扑变化情况下，均能获得满足维数要求的最优断点集；(2)主要计算量集中在追踪被受影响的断点覆盖的回路过程中，剔除断点的顺序不同，将造成计算速度的不同，本发明方法通过自动分区和追踪受影响的断点集等措施，使得电网规模降低，待搜索回路的断点数减少，计算速度明显提高。

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.一种继电保护在线断点求取方法，其特征在于，包括如下步骤：

(1)初始化网络拓扑信息，形成全网主后备保护配合矩阵P和节点关联矩阵V；

(2)搜寻全网割节点，对全网进行分区；

(3)确定全网最小断点集的维数；

(4)求取各分区的离线断点，合并各分区的断点得到原始全网最小断点集BPS_old；

(5)判断网络拓扑是否发生变化，是则跳至步骤(7)；否则顺序执行步骤(6)；

(6)等待一个时间周期，返回步骤(5)；

(7)在网络拓扑变化发生在分区内部时，不更新全网分区，对全网最小断点集的维数进行修正；在网络拓扑变化发生在两个分区之间时，将这两个分区合并后，对合并后的区域进行割节点搜索和重新分区，计算新分区的最小断点集的维数，再根据分区数目的改变对全网最小断点集的维数进行修正；

(8)更新全网主后备保护配合矩阵P，得到网络拓扑变化后的主后备保护配合矩阵P'，求取最终断点集。

2.如权利要求1所述的继电保护在线断点求取方法，其特征在于，所述步骤(1)中，主后备保护配合矩阵P的矩阵规模为n×n，n为网络中的保护总数，P中的第i行第j列元素P(i,j)为：

其中，i＝1,2,…,n，j＝1,2,…,n；

节点关联矩阵V的矩阵规模为m×m，m为网络中的节点总数，V中的第i'行第j'列元素V(i',j')为：

其中，i'＝1,2,…,m，j'＝1,2,…,m。

3.如权利要求2所述的继电保护在线断点求取方法，其特征在于，所述步骤(2)进一步包括如下步骤：

(2-1)搜寻终端割节点，进一步包括如下步骤：

(2-1-1)分别计算节点关联矩阵V的各行之和R_sum(i')；

(2-1-2)获取R_sum(i')＝1时对应的i'值，

(2-1-3)搜索第i'行，得到V(i',j')＝1时对应的j'的值，删除节点关联矩阵V的第i'行和第i'列；

(2-1-4)判断节点关联矩阵V的第j'行之和R_sum(j')是否为1，是则将j'的值赋给i'，返回步骤(2-1-3)；否则判断节点j'是割节点，将节点j'加入割节点集G，顺序执行步骤(2-1-5)；

(2-1-5)判断节点关联矩阵V是否存在和为1的行，是则返回步骤(2-1-2)；否则过程结束；

(2-2)搜寻非终端割节点，进一步包括如下步骤：

(2-2-1)选择节点关联矩阵V中的第一个节点作为当前节点；

(2-2-2)将当前节点所在行中所有不为1的元素存入临时起始链表；

(2-2-3)遍历临时起始链表，判断剔除当前节点后，临时起始链表中元素的连通性，若该链表中的所有元素都连通，则当前节点不是割节点，跳至步骤(2-2-5)；若不连通，则当前节点是割节点，顺序执行步骤(2-2-4)；

(2-2-4)将当前节点存入割节点集G；

(2-2-5)判断是否遍历节点关联矩阵V中的所有节点，是则过程结束；否则选择下一节点作为当前节点，返回步骤(2-2-2)；

(2-3)断开割节点集G中的割节点，得到电网的自动分区。

4.如权利要求1至3中任一项所述的继电保护在线断点求取方法，其特征在于，所述步骤(3)中，对于不可分连通网络，全网最小断点集的维数N_BPS＝N_Z-N_V+2；对于可分连通网络，全网最小断点集的维数N_BPS＝N_Z-N_V+k+1；其中，N_Z为网络支路数，N_V为网络节点数，k为可分连通网络能分为独立分离的连通子图的个数。

5.如权利要求1至4中任一项所述的继电保护在线断点求取方法，其特征在于，所述步骤(8)进一步包括如下步骤：

(8-1)判断拓扑变化是否为新增支路，是则顺序执行步骤(8-2)；否则为删除支路，跳至步骤(8-3)；

(8-2)新增支路两端保护构成集合BPS₁，基于网络拓扑变化后的主后备保护配合矩阵P'，追踪被BPS₁覆盖但未被原始全网最小断点集BPS_old覆盖的回路集L_new，并求取回路集L_new的断点集作为新增断点集BPS_new，合并新增断点集BPS_new和原始全网最小断点集BPS_old，得到临时最小断点集BPS_temp，跳至步骤(8-4)；

(8-3)在删除支路为原始断点所在支路时，直接在原始全网最小断点集BPS_old中删除该断点，得到临时最小断点集BPS_temp；在删除支路为非断点支路时，将原始全网最小断点集BPS_old作为临时最小断点集BPS_temp，同时，选取BPS_temp中的任一断点作为新增断点集BPS_new，BPS_temp中其他断点构成集合BPS₂，追踪被BPS_new覆盖但未被BPS₂覆盖的回路集L_new；

(8-4)将原始全网最小断点集BPS_old中受拓扑变化影响的断点存入集合BPS_delete中；

(8-5)判断临时最小断点集BPS_temp的维数是否等于修正后的全网最小断点集的维数，是则将临时最小断点集BPS_temp作为最终断点集BPS_final，过程结束；否则顺序执行步骤(8-6)；

(8-6)剔除BPS_delete中的部分断点，得到新的集合BPS'_delete，剔除的断点构成集合BPS₃，追踪被BPS₃覆盖但未被BPS'_delete中断点覆盖的回路集L_delete，将L_delete与L_new合并，形成新的回路集L′，求取新的回路集L′的最小断点集BPS_new′，合并BPS_new′和BPS'_delete，得到新的临时最小断点集BPS'_temp，令BPS_temp＝BPS'_temp，返回步骤(8-5)。

6.如权利要求5所述的继电保护在线断点求取方法，其特征在于，所述步骤(8-2)、(8-3)和(8-6)中，按如下步骤追踪被A^*覆盖但未被B^*覆盖的回路集L：

(A)初始化访问标志位数组CF中的所有元素为0；其中，CF中元素的个数等于拓扑变化后的网络的保护数；

(B)令r＝1；

(C)以A^*中的保护start(r)为起始追踪保护，将其对应的保护编号值赋给变量x，其在CF中的标志位CF(x)置1，并将其添加到临时回路链表LOOP_temp；

(D)令y＝1；

(E)判断网络拓扑变化后的主后备保护配合矩阵P'中第x行的第y个元素P'(x,y)是否为1，且CF(y)是否为0，是则顺序执行步骤(F)；否则跳至步骤(H)；

(F)判断保护y是否为B^*中的断点，是则跳至步骤(H)，否则顺序执行步骤(G)；

(G)判断保护y是否为起始追踪保护start(r)，是则将临时回路链表LOOP_temp添加到回路链表集LOOP，顺序执行步骤(H)；否则将保护y添加到临时回路链表LOOP_temp，并将其在CF中的标志位CF(y)置1，将y赋值给x，返回步骤(D)；

(H)判断第x行的所有元素是否遍历完毕，未遍历完毕则y＝y+1，返回步骤(E)；遍历完毕则判断当前临时回路链表LOOP_temp的最后一个元素是否为start(r)，是则剔除该元素；

(I)判断A^*中的保护是否遍历完毕，未遍历完毕则r＝r+1，返回步骤(C)；遍历完毕则过程结束，由回路链表集LOOP得到被A^*覆盖但未被B^*覆盖的回路集L；

其中，对所述步骤(8-2)，A^*为BPS₁，B^*为BPS_old，L为L_new；对所述步骤(8-3)，A^*为BPS_new，B^*为BPS₂，L为L_new；对所述步骤(8-6)，A^*为BPS₃，B^*为BPS'_delete，L为L_delete。