基于抽象对偶网络和连锁故障超网络的脆弱线路辨识方法

1.本发明涉及电力系统安全稳定运行技术领域，特别是涉及到一种基于抽象对偶网络和连锁故障超网络的脆弱线路辨识方法。


背景技术：

2.电力系统作为现代社会覆盖面最广的公共基础性网络之一，为所有系统稳定运行提供电力保障。随着全球经济的不断增长，大电网向着远距离、特高压以及智能互联的方向发展，网络规模结构不断复杂化，电力网络运行经济性和高效性得到提高。但在电力网络迅速发展的同时，网络结构复杂多变使得网络存在固有的脆弱环节，一旦发生故障并引发一系列的连锁反应，将会造成严重的大停电事故。因此加强电网安全稳定分析对预防大停电事故尤为重要，评估电网脆弱性，并找出其薄弱环节，是建立“坚强智能电网”的前提与基础。随着复杂网络理论研究的不断发展及其研究结果的普遍应用，该方法突破了传统方法分析连锁故障的局限性，为研究大电网脆弱性提供了新的视角和方向。
3.目前，针对电力系统脆弱性研究可分为两个方面，结构脆弱性和状态脆弱性。对于电力系统结构脆弱性的研究主要是基于复杂网络理论，对电力网络拓扑结构存在的固有脆弱性进行研究，建立结构脆弱性评估指标；电力系统状态脆弱性评估，主要是分析系统运行状态，重点考虑电力网络运行时各状态参数的变化情况，以及系统节点线路发生连锁故障对全局和局部网络所造成的影响，建立状态脆弱性评估指标。
4.但是，纯粹的电网结构脆弱性评估指标往往忽略了系统的电气特性和物理本质，仅依靠复杂网络指标对电力网络线路进行脆弱辨识，不能真实地反映输电线路在电网中的关键作用。虽然分析系统运行状态建立状态脆弱性评估指标能真实有效地反应元件在电网中所占据的作用，但是分析高纬度、非线性的大规模复杂电网具有较高的复杂度，很难实现在线评估，且忽略了复杂电网拓扑结构固有的脆弱性。因此，电力系统脆弱线路辨识需要进一步综合考虑电网拓扑结构和运行状态，以提高输电线路脆弱性评估精确度，搭建一种新的电力网络输电线路脆弱性评估体系具有重要的理论学术和工程应用意义。


技术实现要素：

5.为克服现有技术中存在的不足，本发明提出了一种基于抽象对偶网络和连锁故障超网络的脆弱线路辨识方法，该方法通过综合考虑抽象对偶网络和连锁故障超网络的结构特性和状态特性，建立互网络节点脆弱度对脆弱线路进行有效辨识。
6.本发明的目的可以通过以下技术方案来实现：
7.基于抽象对偶网络和连锁故障超网络的脆弱线路辨识方法，所述辨识方法具体包括以下步骤：
8.s1：利用复杂网络理论将电力网络简化为网络拓扑结构模型，搭建抽象对偶网络，获取抽象对偶网络邻接矩阵；
9.s2：选取发电机与负荷之间的供需平衡规则以及输电线路中断规则作为电力网络
连锁故障的运行机理，遍历电力网络线路作为初始故障，获取引发的连锁故障过程中的过载失效线路和系统甩负荷；
10.s3：根据s1中搭建的抽象对偶网络以及s2中遍历所有线路作为初始故障获取的连锁故障过程中的过载失效线路，搭建连锁故障超网络模型；
11.s4：根据s3中搭建的连锁故障超网络模型，获取连锁故障超网络关联矩阵，进而得到连锁故障超网络全连接矩阵以及连锁故障超网络邻接矩阵；
12.s5：根据s4中连锁故障超网络关联矩阵、连锁故障超网络全连接矩阵以及连锁故障超网络邻接矩阵，分析所述连锁故障超网络模型的网络拓扑特性，分别计算连锁故障超网络中节点的节点超度，连锁故障超网络中节点的节点度以及连锁故障超网络中节点的节点重合度；
13.s6：采用复杂网络理论粗粒化分解k-shell算法，基于抽象对偶网络邻接矩阵以及连锁故障超网络全连接矩阵，分别对抽象对偶网络以及连锁故障超网络进行k-shell分解，获取所述抽象对偶网络和连锁故障超网络中对应节点的k-shell分层位置；
14.s7：考虑所述抽象对偶网络以及连锁故障超网络中节点的全局特性，获取s6中各k-shell分层位置所分配的全局特性节点k核中心值；
15.s8：考虑所述抽象对偶网络中节点的局部特性，分析抽象对偶网络节点邻域节点集对该节点的依靠性，定义局部特性节点邻域贡献度；
16.s9：综合考虑全局特性节点k核中心值以及局部特性节点邻域贡献度，定义抽象对偶网络节点重要因子；
17.s10：根据s5中连锁故障超网络中节点的节点超度，连锁故障超网络中节点的节点度以及连锁故障超网络中节点的重合度，定义节点故障效率因子；
18.s11：根据s7中全局特性节点k核中心值以及s10中节点故障效率因子，定义连锁故障超网络节点重要因子；
19.s12：根据s9中抽象对偶网络节点重要因子以及s11中连锁故障超网络节点重要因子，考虑电力网络传输能力(网络效率)以及故障后系统甩负荷情况，定义互网络节点脆弱度；
20.s13：根据所述互网络节点脆弱度大小，对电力网络中的脆弱输电线路进行辨识。
21.进一步的，将电力网络中的输电线路和变压器支路抽象为网络拓扑结构中的节点，电力网络中支路邻接关系抽象为网络拓扑结构中的边，搭建的抽象对偶网络表示为：
22.g＝(v,e)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
23.式(1)中，v为抽象对偶网络中的节点集合，表示电力网络中的输电线路和变压器支路集合；
24.v＝{vi|vi＝bi,i＝1,2...,nb}
ꢀꢀꢀꢀꢀ
(2)
25.式(2)中，vi为抽象对偶网络中的节点，bi表示电力网络中的输电线路或变压器支路，nb表示电力网络中的输电支路总数；
[0026][0027]
式(3)中，e
ij
表示当电力网络中支路bi与支路bj存在邻接关系，则抽象对偶网络中节点vi和节点vj存在连边e
ij
，b为电力网络中的支路集合；
[0028]
进一步的，所述发电机与负荷之间的供需平衡规则包括：电力网络中发电机和负荷的供需量减少由一个共同因子决定，在发电机和负荷组件中，如果发电机发出的有功大于负荷的有功需求，则会限制发电机的有功输出，反之，如果发电机发出的有功不足以满足负荷的有功需求，则会进行甩负荷处理；
[0029]
所述输电线路中断规则包括：
[0030]
根据公式(4)，判断输电支路运行状态os
bi
：
[0031][0032]
式(4)中，表示输电支路的潮流大小，表示输电支路所允许的最大传输容量，当时，表示该输电支路处于正常运行状态，时，表示该输电支路应遵循输电线路中断规则进行切除处理，则该线路中断。
[0033]
进一步的，将电力网络中的输电线路和变压器支路抽象为网络拓扑结构中的节点，将连锁故障后提取的失效线路构成连锁故障超网络超边，进而搭建连锁故障超网络模型，所述连锁故障超网络模型可表示为：
[0034]
h＝(vs,es)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0035]
式(5)中，vs为连锁故障超网络的节点集合，可表示为：
[0036][0037]
式(6)中，为连锁故障超网络的节点，与网络对偶抽象图的节点物理本质相同，都表示电力网络中的输电支路；
[0038]
式(5)中，es为连锁故障超网络中的超边集合，可表示为：
[0039][0040]
式(7)中，为连锁故障超网络超边，i＝1～n，n为连锁故障超网络的超边总数，可表示为：
[0041][0042]
的物理本质是电力网络某条线路故障后，引发一系列的连锁故障，可能会导致其他线路存在过载失效的情况，此时提取连锁故障后直至系统稳定时所有过载失效的线路作为连锁故障超网络中的一条超边，该超边中所有节点将具有很大的连锁故障传播相关性；每条超边是：遍历线路作为初始故障，引发的连锁故障链，其中失效的线路即节点就作为超边里的节点，每条超边相当于是失效节点集。
[0043]
进一步的，所述连锁故障超网络关联矩阵ms可表示为：
[0044][0045]
式(9)中，m
ij
是ms中第i行第j列的元素，nb为连锁故障超网络的节点总数，若连锁故障超网络中节点包含在连锁故障超网络超边中，则有m
ij
＝1，否则m
ij
＝0；
[0046]
所述连锁故障超网络全连接矩阵cs可表示为：
[0047][0048]
式(10)中，c
ij
是cs中第i行第j列的元素，所述连锁故障超网络中同一条连锁故障超网络超边内的所有节点全连接后，形成全连接图，当连锁故障超网络中节点与节点存在于同一条连锁故障超网络超边时，则有c
ij
＝1，否则c
ij
＝0；
[0049]
所述连锁故障超网络邻接矩阵as可表示为：
[0050][0051]
式(11)中，a
ij
为as中第i行第j列的元素，d为连锁故障超网络中节点的节点度对角矩阵，a
ij
为连锁故障超网络中包含节点与节点的连锁故障超网络超边的数量。
[0052]
进一步的，所述连锁故障超网络中节点的节点超度的计算公式为：
[0053][0054]
所述连锁故障超网络中节点的节点度的计算公式为：
[0055][0056]
所述连锁故障超网络中节点的重合度的计算公式为：
[0057][0058]
进一步的，所述粗粒化分解k-shell算法流程包括以下步骤：
[0059]
s1：确定初始网络的初始最小度值；
[0060]
s2：删除度值等于初始网络的初始最小度值的节点及其连边；
[0061]
s3：重复s2，直至网络中不再出现度值等于初始网络的初始最小度值的节点，此时将所有被删节点归为1-shell层，并为其分配k核中心值，1-shell层节点的k核中心值等于初始网络的初始最小度值；
[0062]
s4：确定1-shell层后，更新网络情况形成一次刷新网络，获取一次刷新网络的初始最小度值；
[0063]
s5：删除度值等于一次刷新网络的初始最小度值的节点及其连边；
[0064]
s6：重复s5，将所删节点归为2-shell层，并为其分配k核中心值，2-shell层节点的k核中心值等于一次刷新网络的初始最小度值；
[0065]
s7：重复s1～s6，直至k次刷新网络后所有节点均被k-shell分层并分配k核中心值。
[0066]
进一步的，所述节点领域贡献度表现了节点在其邻域内的关键性，如果移除该节点会使得邻域内其他节点失去信息交流的枢纽通道，所述局部特性节点邻域贡献度d
vi,con
的计算公式为：
[0067][0068]
式(15)中，nh表示节点vi的邻域节点集合，包括节点vi，节点vi的邻居节点集合以及节点vi的次邻居节点集合，表示邻域节点集合内，节点之间的信息传输流总和，由于本发明建立的邻域贡献度是从网络拓扑的角度出发，所以在步长小于等于2的两个节点vn和vm之间就存在信息传输流表示将邻域节点集合nh中的vi节点移除后网络刷新，得到的邻域节点集合内的信息传输流总和。
[0069]
所述抽象对偶网络节点重要因子表现了从网络结构拓扑出发，综合考虑网络拓扑的全局特性k核中心值和局部特性节点邻域贡献度，所述抽象对偶网络节点重要因子的计算公式为：
[0070][0071]
为抽象对偶网络中节点vi的k核中心值，ks
max
为抽象对偶网络所有节点中的最大k核中心值，d
max,con
和d
min,con
分别为抽象对偶网络所有节点中的最大邻域贡献度和最小邻域贡献度。
[0072]
进一步的，所述节点故障效率因子综合考虑了连锁故障超网络中节点在超网络超边中出现的次数，节点与其他节点在连锁故障中的交互关系以及节点在连锁故障中的活跃影响程度，所述节点故障效率因子的计算公式为：
[0073][0074]
所述连锁故障超网络节点重要因子的计算公式为：
[0075][0076]
式(18)中，为连锁故障超网络中节点的k核中心值，分别表示连锁故障超网络所有节点中的最大k核中心值。
[0077]
进一步的，所述互网络节点脆弱度vdi的计算公式为：
[0078][0079]
式(19)中，ag为抽象对偶网络移除传输支路前的网络效率，为抽象对偶网络移除节点vi后的网络效率，为连锁故障超网络发生故障之前的系统总负荷，为连锁故障超网络中节点故障直至系统稳定后的总负荷；根据所述互网络节点脆弱度vdi的大小进行降序排序，获取排序靠前的脆弱线路，添加至重点保护线路集合中。
[0080]
与现有的技术相比，本发明有益的效果为：
[0081]
(1)所建立的连锁故障超网络摒弃了传统的线路邻接关系，选取连锁故障运行机理提取过载失效线路构成超边集合，更能直接捕获连锁故障中的节点传播相关性；
[0082]
(2)采用复杂网络理论k-shell算法，结合全局和局部拓扑特性建立抽象对偶网络节点重要因子和连锁故障超网络节点重要因子，多层次、多网络、多尺度地构建捕获脆弱性的重要因子，使得脆弱性辨识更为准确；
[0083]
(3)结合电力网络的运行特征与拓扑特征构建更为接近电力系统真实特性的脆弱线路辨识。
附图说明
[0084]
图1为本发明流程图；
[0085]
图2为ieee14节点抽象对偶网络图；
[0086]
图3为ieee14节点测试系统中20条输电线路的互网络节点脆弱度图。
具体实施方式
[0087]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
[0088]
如图1所示，基于抽象对偶网络和连锁故障超网络的脆弱线路辨识方法，具体包括以下步骤：
[0089]
基于抽象对偶网络和连锁故障超网络的脆弱线路辨识方法，所述辨识方法具体包括以下步骤：
[0090]
s1：利用复杂网络理论将电力网络简化为网络拓扑结构模型，搭建抽象对偶网络，获取抽象对偶网络邻接矩阵；
[0091]
s2：选取发电机与负荷之间的供需平衡规则以及输电线路中断规则作为电力网络连锁故障的运行机理，遍历电力网络线路作为初始故障，获取引发的连锁故障过程中的过载失效线路和系统甩负荷；
[0092]
s3：根据s1中搭建的抽象对偶网络以及s2中遍历所有线路作为初始故障获取的连锁故障过程中的过载失效线路，搭建连锁故障超网络模型；
[0093]
s4：根据s3中搭建的连锁故障超网络模型，获取连锁故障超网络关联矩阵，进而得到连锁故障超网络全连接矩阵以及连锁故障超网络邻接矩阵；
[0094]
s5：根据s4中连锁故障超网络关联矩阵、连锁故障超网络全连接矩阵以及连锁故障超网络邻接矩阵，分析所述连锁故障超网络模型的网络拓扑特性，分别计算连锁故障超网络中节点的节点超度，连锁故障超网络中节点的节点度以及连锁故障超网络中节点的节点重合度；
[0095]
s6：采用复杂网络理论粗粒化分解k-shell算法，基于抽象对偶网络邻接矩阵以及连锁故障超网络全连接矩阵，分别对抽象对偶网络以及连锁故障超网络进行k-shell分解，获取所述抽象对偶网络和连锁故障超网络中对应节点的k-shell分层位置；
[0096]
s7：考虑所述抽象对偶网络以及连锁故障超网络中节点的全局特性，获取s6中各k-shell分层位置所分配的全局特性节点k核中心值；
[0097]
s8：考虑所述抽象对偶网络中节点的局部特性，分析抽象对偶网络节点邻域节点集对该节点的依靠性，定义局部特性节点邻域贡献度；
[0098]
s9：综合考虑全局特性节点k核中心值以及局部特性节点邻域贡献度，定义抽象对偶网络节点重要因子；
[0099]
s10：根据s5中连锁故障超网络中节点的节点超度，连锁故障超网络中节点的节点度以及连锁故障超网络中节点的重合度，定义节点故障效率因子；
[0100]
s11：根据s7中全局特性节点k核中心值以及s10中节点故障效率因子，定义连锁故障超网络节点重要因子；
[0101]
s12：根据s9中抽象对偶网络节点重要因子以及s11中连锁故障超网络节点重要因子，考虑电力网络传输能力(网络效率)以及故障后系统甩负荷情况，定义互网络节点脆弱度；
[0102]
s13：根据所述互网络节点脆弱度大小，对电力网络中的脆弱输电线路进行辨识。
[0103]
将电力网络中的输电线路和变压器支路抽象为网络拓扑结构中的节点，电力网络中支路邻接关系抽象为网络拓扑结构中的边，搭建的抽象对偶网络表示为：
[0104]
g＝(v,e)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0105]
式(1)中，v为抽象对偶网络中的节点集合，表示电力网络中的输电线路和变压器支路集合；
[0106][0107]
式(2)中，vi为抽象对偶网络中的节点，bi表示电力网络中的输电线路或变压器支路，nb表示电力网络中的输电支路总数；
[0108][0109]
式(3)中，e
ij
表示当电力网络中支路bi与支路bj存在邻接关系，则抽象对偶网络中节点vi和节点vj存在连边e
ij
，b为电力网络中的支路集合。
[0110]
如图2所示为ieee14节点抽象对偶网络图，本发明选用ieee14节点测试系统进行抽象对偶网络的搭建，ieee14节点测试系统的输电支路数量优选为20，所搭建的抽象对偶网络节点数量为20；输电支路的邻接关系抽象为节点之间的连边，所搭建的抽象对偶网络连边数量为46。
[0111]
所述发电机与负荷之间的供需平衡规则包括：电力网络中发电机和负荷的供需量减少由一个共同因子决定，在发电机和负荷组件中，如果发电机发出的有功大于负荷的有功需求，则会限制发电机的有功输出，反之，如果发电机发出的有功不足以满足负荷的有功需求，则会进行甩负荷处理；
[0112]
所述输电线路中断规则包括：
[0113]
根据公式(4)，判断输电支路运行状态
[0114][0115]
式(4)中，表示输电支路的潮流大小，表示输电支路所允许的最大传输容量，当时，表示该输电支路处于正常运行状态，时，表示该输电支路应遵循输电线路中断规则进行切除处理，则该线路中断。
[0116]
将电力网络中的输电线路和变压器支路抽象为网络拓扑结构中的节点，将连锁故障后提取的失效线路构成连锁故障超网络超边，进而搭建连锁故障超网络模型，所述连锁故障超网络模型可表示为：
[0117]
h＝(vs,es)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0118]
式(5)中，vs为连锁故障超网络的节点集合，可表示为：
[0119][0120]
式(6)中，为连锁故障超网络的节点，与网络对偶抽象图的节点物理本质相同，都表示电力网络中的输电支路；
[0121]
式(5)中，es为连锁故障超网络中的超边集合，可表示为：
[0122][0123]
式(7)中，为连锁故障超网络超边，i＝1～n，n为连锁故障超网络的超边总数，可表示为：
[0124][0125]
的物理本质是电力网络某条线路故障后，引发一系列的连锁故障，可能会导致其他线路存在过载失效的情况，此时提取连锁故障后直至系统稳定时所有过载失效的线路作为连锁故障超网络中的一条超边，该超边中所有节点将具有很大的连锁故障传播相关性；每条超边是：遍历线路作为初始故障，引发的连锁故障链，其中失效的线路即节点就作为超边里的节点，每条超边相当于是失效节点集。
[0126]
本发明针对ieee14节点测试系统的所有输电支路进行遍历连锁故障仿真，每条支路作为初始故障时，引发一系列的连锁反应，提取系统稳定后的所有过载失效线路作为连锁故障超网络的超边。
[0127]
表1为ieee14节点测试系统连锁故障超网络中的超边集合表，由于复杂网络理论中超边的节点数量大于等于2，因此总共存在9条连锁故障超边。
[0128]
表1.ieee14节点连锁故障超网络超边集合表
[0129][0130]
所述连锁故障超网络关联矩阵ms可表示为：
[0131][0132]
式(9)中，m
ij
是ms中第i行第j列的元素，nb为连锁故障超网络的节点总数，若连锁故障超网络中节点包含在连锁故障超网络超边中，则有m
ij
＝1，否则m
ij
＝0；
[0133]
所述连锁故障超网络全连接矩阵cs可表示为：
[0134][0135]
式(10)中，c
ij
是cs中第i行第j列的元素，所述连锁故障超网络中同一条连锁故障超网络超边内的所有节点全连接后，形成全连接图，当连锁故障超网络中节点与节点存在于同一条连锁故障超网络超边时，则有c
ij
＝1，否则c
ij
＝0；
[0136]
所述连锁故障超网络邻接矩阵as可表示为：
[0137][0138]
式(11)中，a
ij
为as中第i行第j列的元素，d为连锁故障超网络中节点的节点度对角矩阵，a
ij
为连锁故障超网络中包含节点与节点的连锁故障超网络超边的数量。
[0139]
所述连锁故障超网络中节点的节点超度的计算公式为：
[0140][0141]
所述连锁故障超网络中节点的节点度的计算公式为：
[0142][0143]
所述连锁故障超网络中节点的重合度的计算公式为：
[0144]
[0145]
所述粗粒化分解k-shell算法流程包括以下步骤：
[0146]
s1：确定初始网络的初始最小度值；
[0147]
s2：删除度值等于初始网络的初始最小度值的节点及其连边；
[0148]
s3：重复s2，直至网络中不再出现度值等于初始网络的初始最小度值的节点，此时将所有被删节点归为1-shell层，并为其分配k核中心值，1-shell层节点的k核中心值等于初始网络的初始最小度值；
[0149]
s4：确定1-shell层后，更新网络情况形成一次刷新网络，获取一次刷新网络的初始最小度值；
[0150]
s5：删除度值等于一次刷新网络的初始最小度值的节点及其连边；
[0151]
s6：重复s5，将所删节点归为2-shell层，并为其分配k核中心值，2-shell层节点的k核中心值等于一次刷新网络的初始最小度值；
[0152]
s7：重复s1～s6，直至k次刷新网络后所有节点均被k-shell分层并分配k核中心值。
[0153]
所述节点领域贡献度表现了节点在其邻域内的关键性，如果移除该节点会使得邻域内其他节点失去信息交流的枢纽通道，所述局部特性节点邻域贡献度的计算公式为：
[0154][0155]
式(15)中，nh表示节点vi的邻域节点集合，包括节点vi，节点vi的邻居节点集合以及节点vi的次邻居节点集合，表示邻域节点集合内，节点之间的信息传输流总和，由于本发明建立的邻域贡献度是从网络拓扑的角度出发，所以在步长小于等于2的两个节点vn和vm之间就存在信息传输流表示将邻域节点集合nh中的vi节点移除后网络刷新，得到的邻域节点集合内的信息传输流总和。
[0156]
所述抽象对偶网络节点重要因子表现了从网络结构拓扑出发，综合考虑网络拓扑的全局特性k核中心值和局部特性节点邻域贡献度，所述抽象对偶网络节点重要因子的计算公式为：
[0157][0158]
为抽象对偶网络中节点vi的k核中心值，ks
max
为抽象对偶网络所有节点中的最大k核中心值，d
max,con
和d
min,con
分别为抽象对偶网络所有节点中的最大邻域贡献度和最小邻域贡献度。
[0159]
所述节点故障效率因子综合考虑了连锁故障超网络中节点在超网络超边中出现的次数，节点与其他节点在连锁故障中的交互关系以及节点在连锁故障中的活跃影响程
度，所述节点故障效率因子的计算公式为：
[0160][0161]
所述连锁故障超网络节点重要因子的计算公式为：
[0162][0163]
式(18)中，为连锁故障超网络中节点的k核中心值，分别表示连锁故障超网络所有节点中的最大k核中心值。
[0164]
如表2所示为ieee14节点测试系统中输电线路所对应的抽象对偶网络节点重要因子和连锁故障超网络节点重要因子。
[0165]
表2.ieee14节点抽象对偶网络和连锁故障超网络节点重要因子
[0166][0167][0168]
所述互网络节点脆弱度vdi的计算公式为：
[0169][0170]
式(19)中，ag为抽象对偶网络移除传输支路前的网络效率，为抽象对偶网络移除节点vi后的网络效率，为连锁故障超网络发生故障之前的系统总负荷，为连锁故障超网络中节点故障直至系统稳定后的总负荷；根据所述互网络节点脆弱度vdi的大小进行降序排序，获取排序靠前的脆弱线路，添加至重点保护线路集合中。
[0171]
本发明根据所搭建的互网络节点脆弱度评估指标对ieee14节点测试系统的输电线路进行脆弱性评估，仿真计算电力网络中每条输电线路的互网络节点脆弱度，得到如图3
所示的ieee14节点测试系统中20条输电线路的互网络节点脆弱度，并进行降序排序，将排名前10的线路添加至重点保护线路集合中，如表3所示为ieee14节点测试系统中的重点保护线路集合的互网络节点脆弱度排序结果。
[0172]
表3 ieee14节点测试系统中重点保护线路集合的互网络节点脆弱度排序
[0173][0174]
与参考文献[1]线路综合脆弱度表中排序前10的线路对比可知，本发明得到的重点保护线路集合存在5条相同线路；与参考文献[2]线路综合脆弱度表中排序前10的线路对比，本发明得到的重点保护线路集合存在7条相同线路；并且本发明所获得的重点保护线路集合中的线路都属于产生9条超边中的初始故障线路。因此，对本发明所获得的重点保护线路集合加以重点保护，能够很好的预防系统连锁故障传播影响范围放大。仿真结果表明，采用本发明提出的互网络节点脆弱度评估模型能准确评估电力网络输电线路的脆弱性，可有效运用于电力系统安全保护中。
[0175]
参考文献包括：
[0176]
[1]陈召阳,吕鹏飞.考虑运行方式的电网关键线路辨识.电力系统及其自动化学报,2014.
[0177]
[2]parviz khaledian,brian k. johnson,saied hemati.power grid security improvement by remedial action schemes using vulnerability assessment based on fault chains and power flow,2018ieee international conference on probabilistic methods applied to power systems(pmaps),2018.