CN107482626B - 一种区域电网关键节点识别方法 - Google Patents

一种区域电网关键节点识别方法 Download PDF

Info

Publication number
CN107482626B
CN107482626B CN201710708797.4A CN201710708797A CN107482626B CN 107482626 B CN107482626 B CN 107482626B CN 201710708797 A CN201710708797 A CN 201710708797A CN 107482626 B CN107482626 B CN 107482626B
Authority
CN
China
Prior art keywords
node
nodes
power
power network
index
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
Application number
CN201710708797.4A
Other languages
English (en)
Other versions
CN107482626A (zh
Inventor
吕云锋
韩畅
林冠强
叶晓君
莫天文
李惠松
王志奎
赵昱宣
林振智
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Huizhou Power Supply Bureau of Guangdong Power Grid Co Ltd
Original Assignee
Zhejiang University ZJU
Huizhou Power Supply Bureau of Guangdong Power Grid Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU, Huizhou Power Supply Bureau of Guangdong Power Grid Co Ltd filed Critical Zhejiang University ZJU
Priority to CN201710708797.4A priority Critical patent/CN107482626B/zh
Publication of CN107482626A publication Critical patent/CN107482626A/zh
Application granted granted Critical
Publication of CN107482626B publication Critical patent/CN107482626B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Secondary Cells (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种区域电网关键节点识别方法,基本实施流程包括步骤:输入电力系统拓扑结构和电气特性参数;计算节点在各评价指标上的指标值;根据指标值得到节点重要度评估决策矩阵;生成标准化区域电网的节点重要度评估决策矩阵;考虑基尼系数和肯德尔相关系数,采用CRITIC法求取指标客观权重;采用TOPSIS法求取节点重要度的评估结果。本发明提供的技术方案充分考虑了网络拓扑结构的物理特征和电网运行状态的参数。所提出的方法能够较好地识别出电网的关键节点,避免了人为确定权重的主观性,且评估结果更加贴近电网运行的实际情况。

Description

一种区域电网关键节点识别方法
技术领域
本发明涉及电力系统的技术领域,更具体地,涉及一种区域电网关键节点识别方法。
背景技术
随着电力系统装机容量的不断增大、区域电网间的互联以及具有随机性和不确定性的高比例可再生能源的并网,使得电力系统运行和维护的复杂性不断增大、电力系统经常运行于更接近极限的状态。此外,电力系统也不可避免地面对极端天气和自然灾害的挑战。这些因素导致电力系统发生大面积停电事故的可能性有所增大。因此,准确快速地识别出电网关键节点,对于预防电网发生大面积停电事故具有非常重要的意义。
目前,在电网关键节点识别方面,目前已提出了一些指标和方法,例如有功网络流介数指标、无功支撑度和无功贡献度指标、节点电气耦合连接度等指标、抗灾变性指标以及重要度评价矩阵、潮流追踪、正规形方法、PageRank、熵权法等方法,但这些指标或算法大多是从单个或两个方面进行节点重要度评价,未能从多个角度或维度对节点进行综合评估。
可见,现有关于区域电网关键节点识别的方法还有待改进。
发明内容
本发明为克服上述现有技术所述的至少一种缺陷,提供一种区域电网关键节点识别方法,能够较好地识别出电网的关键节点,避免了人为确定权重的主观性,且评估结果更加贴近电网运行的实际情况。
本发明的技术方案是:基于复杂区域电网的复杂网络拓扑特性和电气特性,提出了评估区域电网关键节点的指标,即电力网络节点的功率集中度、电力网络节点的电气介数、电力网络的传输效能变化率、电力网络的凝聚度变化率、电力网络的生成树变化率;
采用TOPSIS法来解决区域电网的节点重要度评估问题,按照各个节点与理想解的相对逼近度大小,对区域电网的节点进行重要性程度排序,从而得到节点重要度的评估结果;
采用CRITIC法确定TOPSIS模型中指标的客观权重,利用基尼系数和肯德尔相关系数来分别衡量不同评价对象的取值差异性和评价指标之间的冲突性。
具体的,本发明中一种区域电网关键节点识别方法,其中,包括以下步骤:
输入电力系统拓扑结构参数、电气特性参数;
基于区域电网的复杂网络拓扑特性和电气特性,提出了评估区域电网关键节点的指标,即电力网络节点的功率集中度、电力网络节点的电气介数、电力网络的传输效能变化率、电力网络的凝聚度变化率、电力网络的生成树变化率;
计算各节点在各项评价指标下的指标值,生成标准化区域电网的节点重要度评估决策矩阵;
采用理想点法来解决区域电网的节点重要度评估问题,按照各个节点与理想解的相对逼近度大小,对区域电网的节点进行重要性程度排序,从而得到节点重要度的评估结果;
采用CRITIC法确定理想点法模型中指标的客观权重,利用基尼系数和肯德尔相关系数来分别衡量不同评价对象的取值差异性和评价指标之间的冲突性。
进一步的,采用电力网络节点的功率集中度来识别电网关键节点,包括:
在电力网络中,节点的度在拓扑结构方面定义为与该节点连接的输电线路数目;然而对于电力网络而言,节点的重要性程度不仅只与电网的拓扑结构有关,而且还与该节点所连接的重要负荷量或者发电机组的装机容量密切相关;如果某一个节点上所连接的发电机组装机容量越大或重要负荷量越大,则该节点的重要性程度就越大;故综合电力网络的复杂网络拓扑特性及其电气特性,定义了电力网络节点的功率集中度指标:
Figure BDA0001382173300000031
式中:Aij为网络邻接矩阵A的第i行第j列元素;Ω(i)为与节点i直接相连的节点的集合;N为电力网络的节点总数;PG,max为发电机组的最大有功发电出力;PG,i和PL,i分别为节点i上所连接的发电机有功出力和重要有功负荷量;综上,如果某一个节点所连接的有功负荷和发电机组有功出力越大,以及与该节点连接的输电线路数目越多,则该节点的重要性程度就越高。
进一步的,采用电力网络节点的电气介数来识别电网关键节点,包括:
节点的介数定义为网络中所有节点对的最短路径经过该节点的条数;节点的介数越大,表明该节点在整个网络中起到枢纽作用,因此该节点的重要性程度也越高;然而,在电力网络中,电能会通过电源节点和负荷节点之间所有可能的路径进行传输,而不是复杂网络理论中所定义的最短路径;因此,复杂网络理论中所定义的节点介数不能直接用于分析电力网络;在电力网络中,功率传输分布系数用于表征电能在电力系统中的传输特性,该系数不仅考虑了电力网络的拓扑特性,而且还计及了输电线路的电气参数;因此,这里基于功率传输分布系数定义了节点电气介数:
Figure BDA0001382173300000032
Figure BDA0001382173300000033
Figure BDA0001382173300000034
式中:ΩSE和ΩSK分别表示电源节点和负荷节点的集合;
Figure BDA0001382173300000035
表示从电源节点g通过线路i-j传输到负荷节点l的有功功率;xij表示线路i-j的电抗;
Figure BDA0001382173300000036
表示节点i和g之间的等效阻抗;Zig表示从节点i到节点g的转移阻抗;Zii和Zgg分别表示节点i和g之间的输入阻抗。电力网络中节点的电气介数越大,表明该节点的重要性程度越高。
进一步的,采用电力网络传输效能的变化率来识别电网关键节点,包括:
在电力网络中,电力网络传输效能定义为电力网络中节点间最短电气距离倒数的平均值,其表达式为:
Figure BDA0001382173300000041
式中:ΩN为电力网络的节点集合;
Figure BDA0001382173300000042
表示电力网络中节点j和k之间的最短电气距离,即节点j和k之间电能传输路径上输电线路的电抗值之和最小;
电力网络的连通性受到每一个节点的影响,电力网络如果失去部分关键节点将使其连通性以及电能传输效能大大降低;因此,电力网络中节点的重要性程度可以通过移除该节点及其所连接的输电线路前后电力网络传输效能的变化率来衡量;因此,节点i的电力网络传输效能可以定义为该节点移除前后电力网络传输效能的变化率,即
Figure BDA0001382173300000043
Figure BDA0001382173300000044
式中:EB,i表示电力网络在节点i失效后的电力网络传输效能;因此,一个节点的电力网络传输效能变化率越大,表明该节点失效后对电网的电能传输影响越大,则该节点在电网中的重要性程度也就越高。
进一步的,采用电力网络的凝聚度变化率来识别电网关键节点,包括:
在电力网络中,节点的凝聚度定义为电力网络中该节点到达其他所有节点的最短电气距离之和的倒数,该指标反映了该节点在电力网络中位置;因而,电力网络的凝聚度可以定义为所有节点之间最短电气距离之和的倒数,其表达式为:
Figure BDA0001382173300000045
在电力网络中,一个节点的失效对电力网络的凝聚度将产生影响;节点收缩法常被用于衡量节点失效对与其直接相连的邻居节点的影响;假设原始电力网络G中节点i收缩后得到了新的电力网络Gi,则电力网络Gi的凝聚度为
Figure BDA0001382173300000051
式中:NC为电力网络Gi的节点数;节点i的重要性程度可以通过节点i及其邻居节点收缩成一个节点后,新电力网络Gi的凝聚度变化率来衡量,故电力网络的凝聚度变化率可表示为:
Figure BDA0001382173300000052
该指标值越大,表明新电力网络Gi的凝聚度增加得越多,则该节点的重要性程度就越大。
进一步的,采用电力网络的生成树变化率来识别电网关键节点,包括:
在复杂网络中,网络的生成树为包括了所有节点的子图;定义网络G的拉普拉斯矩阵为
Figure BDA0001382173300000053
式中:D为网络的度矩阵,即对角元素为
Figure BDA0001382173300000054
的对角矩阵;根据Kirchhoff矩阵树理论可得,网络的生成树数目等于拉普拉斯矩阵的任何一个N-1阶主子式,即
Figure BDA0001382173300000055
式中:λ12,...,λN-1为拉普拉斯矩阵的非零特征值;
在电力网络中,一个节点的失效使得电力网络的生成树数目减小,甚至可能使得电力网络成为一个不连通的网络;因此,该节点的重要性程度可以通过移除该节点及其直接相连的输电线路前后的生成树数目的变化率来衡量;因此,节点i的生成树数目的变化率可表示为:
Figure BDA0001382173300000056
式中:τ(G-Vi)为图G删除第i节点及与其相连的输电线路后电力网络的生成树数目;τ(G-Vi)的数值越小,IRST(i)的数值越大,则节点i的失效对整个电力网络的影响程度越大,该节点越重要。
进一步的,包括建立多属性决策问题的决策矩阵;在区域电网的节点重要度评估中,假设该区域电网共有N个节点并选取了M个用于衡量节点重要度的评价指标;因此,区域电网的节点重要度评估问题的决策矩阵可以表示为
Figure BDA0001382173300000061
式中:Iij表示节点i对应的指标j的指标值,其中i∈{1,2,...,N},j∈{1,2,...,M};M等于5。
采用TOPSIS法来解决区域电网的节点重要度评估问题,决策流程如下:
步骤1:生成标准化区域电网的节点重要度评估决策矩阵
Y=(yij)N×M
式中:
Figure BDA0001382173300000062
i∈{1,2,...,N},j∈{1,2,...,M};
步骤2:生成加权后的区域电网节点重要度评估矩阵
Z=(zij)N×M
式中:zij=ωjyij;ωj为指标j的权重,即指标权重向量ω=(ω12,…,ωM)的第j个元素,其中ωj∈[0,1],且
Figure BDA0001382173300000063
i∈{1,2,...,N},j∈{1,2,...,M};指标的权重采用CRITIC法确定;
步骤3:确定节点重要度评估问题的理想解
Figure BDA0001382173300000064
和反理想解
Figure BDA0001382173300000065
其中理想解是期望的最重要的电网节点,该节点的各个指标值都比电网中其它各个节点的相应指标的值都好,即
Figure BDA0001382173300000066
反理想解是假想的最不重要的电网节点,该节点的各个指标值都比电网中其它各个节点的相应指标的值都差,即
Figure BDA0001382173300000067
j∈{1,2,...,M};
步骤4:计算区域电网各节点分别与理想解和反理想解的欧氏距离
Figure BDA0001382173300000068
Figure BDA0001382173300000069
Figure BDA00013821733000000610
Figure BDA0001382173300000071
式中:
Figure BDA0001382173300000072
为加权后的节点重要度评估矩阵Z的第i行;
步骤5:计算区域电网中各个节点与理想解的相对逼近度
Figure BDA0001382173300000073
节点与理想解的相对逼近度的取值范围为[0,1],该值愈逼近1,则表明该节点越重要;反之,该值愈逼近0,则表明该节点越不重要;
步骤6:按照各个节点与理想解的相对逼近度大小,对区域电网的节点进行重要性程度排序,从而得到节点重要度的评估结果。
进一步的,采用CRITIC法来确定指标的权重;采用基尼系数和肯德尔系数分别来衡量指标在不同评价对象的取值差异性和评价指标之间的冲突性;包括:
第j个指标的基尼指数可以定义为:
Figure BDA0001382173300000074
该值的取值范围为[0,1],1表示该指标的取值完全不相似,0表示该指标的取值完全相似;在电网的节点重要度评估问题中,用于评估节点重要性的5个指标可以看作为基尼指数计算中的社会成员,而各个指标的属性值可以看作为基尼指数计算中社会成员的财富值;
具有N个元素的两列变量
Figure BDA0001382173300000075
Figure BDA0001382173300000076
其第i个变量值分别为zij和zik(1≤i≤N),其中j∈{1,2,...,M},k∈{1,2,...,M};假设XjR和XkR分别为zij和zik
Figure BDA0001382173300000077
Figure BDA0001382173300000078
的排序值,则XjR和XkR对应排序值变量可以搭配组成变量对的集合XR,该变量对集合的第i个变量对为(XjR,XkR);因此,第j个和k个指标之间的肯德尔相关系数可以定义为:
Figure BDA0001382173300000081
式中:A为常数,其值为N(N-1)/2;Ncc和Ndc分别表示变量对的集合XR中变量排序值相等的变量对数目和变量值排序值不相等的变量对数目;
Figure BDA0001382173300000082
Figure BDA0001382173300000083
分别表示变量
Figure BDA0001382173300000084
Figure BDA0001382173300000085
中具有相同变量值的个数;因此,第j个指标与其他指标的整体肯德尔相关系数可以定义为:
Figure BDA0001382173300000086
当指标j肯德尔系数为1时,则表明该指标与其他指标具有一致的等级相关性;而肯德尔系数为0时,则表明该指标与其他指标是相互独立的;
综合考虑基尼系数和肯德尔系数,即指标j的客观权重可以表示为:
Figure BDA0001382173300000087
与现有技术相比,有益效果是:本发明提供的区域电网关键节点的识别方法,充分考虑了网络拓扑结构的物理特征和电网运行状态的参数。所提出的方法能够较好地识别出电网的关键节点,避免了人为确定权重的主观性,且评估结果更加贴近电网运行的实际情况。
附图说明
图1是本发明整体流程示意图。
图2是本发明广东某区域电网拓扑图。
图3是本发明基于本发明的广东某区域电网节点的重要度评估结果图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。附图中描述位置关系仅用于示例性说明,不能理解为对本专利的限制。
如图1所示,一种区域电网关键节点识别方法,其基本实施流程包括如下详细步骤:
步骤1、输入电力系统拓扑结构参数(线路连接关系)、电气特性参数(节点参数、线路参数);
步骤2、计算各节点在各评价指标上的指标值,具体指标包括:
a.电力网络节点的功率集中度指标IPCD(i):
Figure BDA0001382173300000091
式中:Aij为网络邻接矩阵A的第i行第j列元素;Ω(i)为与节点i直接相连的节点的集合;N为电力网络的节点总数;PG,max为发电机组的最大有功发电出力;PG,i和PL,i分别为节点i上所连接的发电机有功出力和重要有功负荷量。
b.电力网络节点的电气介数指标
Figure BDA0001382173300000092
Figure BDA0001382173300000093
Figure BDA0001382173300000094
Figure BDA0001382173300000095
式中:ΩSE和ΩSK分别表示电源节点和负荷节点的集合;
Figure BDA0001382173300000096
表示从电源节点g通过线路i-j传输到负荷节点l的有功功率;xij表示线路i-j的电抗;
Figure BDA0001382173300000097
表示节点i和g之间的等效阻抗;Zig表示从节点i到节点g的转移阻抗;Zii和Zgg分别表示节点i和g之间的输入阻抗。
c.电力网络的传输效能变化率指标ITER(i):
Figure BDA0001382173300000098
Figure BDA0001382173300000099
Figure BDA00013821733000000910
式中:ΩN为电力网络的节点集合;
Figure BDA00013821733000000911
表示电力网络中节点j和k之间的最短电气距离,即节点j和k之间电能传输路径上输电线路的电抗值之和最小;EB为电力网络的传输效能;EB,i表示电力网络在节点i失效后的电力网络传输效能。
d.电力网络的凝聚度变化率指标IRNC(i):
Figure BDA0001382173300000101
Figure BDA0001382173300000102
Figure BDA0001382173300000103
式中:NC为电力网络Gi的节点数。ΨNC(G)为原电力网络的凝聚度;ΨNC(Gi)为原电力网络G中节点i收缩后得到的新电力网络Gi的凝聚度。
e.电力网络的生成树变化率指标IRST(i):
Figure BDA0001382173300000104
Figure BDA0001382173300000105
Figure BDA0001382173300000106
式中:D为网络的度矩阵,即对角元素为
Figure BDA0001382173300000107
的对角矩阵;λ12,...,λN-1为拉普拉斯矩阵的非零特征值;τ(G-Vi)为图G删除第i节点(即Vi)及与其相连的输电线路后电力网络的生成树数目。
步骤3、根据各节点在各项评价指标下的指标值,得到区域电网的节点重要度评估决策矩阵IDM
Figure BDA0001382173300000108
式中:Iij表示节点i对应的指标j的指标值,其中i∈{1,2,...,N},j∈{1,2,...,M}。在本发明中M等于5。
步骤4、对各项评价指标下的指标值进行标准化处理,生成标准化区域电网的节点重要度评估决策矩阵Y:
Y=(yij)N×M
式中:
Figure BDA0001382173300000111
i∈{1,2,...,N},j∈{1,2,...,M}。
步骤5、综合考虑基尼系数和肯德尔相关系数,采用CRITIC法求取指标的客观权重ωj
Figure BDA0001382173300000112
Figure BDA0001382173300000113
Figure BDA0001382173300000114
Figure BDA0001382173300000115
式中:Gj为第j个指标的基尼指数;τjk为第j个和k个指标之间的肯德尔相关系数;A为常数,其值为N(N-1)/2;Ncc和Ndc分别表示变量对的集合XR中变量排序值相等的变量对数目和变量值排序值不相等的变量对数目;
Figure BDA0001382173300000116
Figure BDA0001382173300000117
分别表示变量
Figure BDA0001382173300000118
Figure BDA0001382173300000119
中具有相同变量值的个数。τj为第j个指标与其他指标的整体肯德尔相关系数。
步骤6、生成加权后的区域电网节点重要度评估矩阵Z
Z=(zij)N×M
式中:zij=ωjyij;ωj为指标j的权重,即指标权重向量ω=(ω12,…,ωM)的第j个元素,其中ωj∈[0,1],且
Figure BDA00013821733000001110
i∈{1,2,...,N},j∈{1,2,...,M}。
步骤7、确定节点重要度评估问题的理想解
Figure BDA00013821733000001111
和反理想解
Figure BDA0001382173300000121
其中理想解
Figure BDA0001382173300000122
反理想解
Figure BDA0001382173300000123
j∈{1,2,...,M}。
步骤8、计算区域电网各节点分别与理想解和反理想解的欧氏距离
Figure BDA0001382173300000124
Figure BDA0001382173300000125
Figure BDA0001382173300000126
Figure BDA0001382173300000127
式中:
Figure BDA0001382173300000128
为加权后的节点重要度评估矩阵Z的第i行。
步骤9、计算区域电网中各个节点与理想解的相对逼近度
Figure BDA0001382173300000129
步骤10、按照各个节点与理想解的相对逼近度大小,对区域电网的节点进行重要性程度排序,从而得到节点重要度的评估结果。
为了进一步理解本发明,以下以简化后的广东某区域电网为例,来解释本发明的实际应用。
简化后的广东某区域电网共有20台发电机组、66个节点、95条输电线路支路,其拓扑图如图2所示,其中,数字表示节点编号。
首先,根据各个指标的定义计算节点各个指标的数值,从而形成区域电网的节点重要度评估问题的决策矩阵IDM,然后对其归一化处理;接着,根据基尼系数和肯德尔系数的定义,分别计算各个指标的基尼系数、肯德尔系数以及客观综合权重,其结果如表1所示。
表1指标的基尼系数、肯德尔系数和客观综合权重
Figure BDA00013821733000001210
从表1可以看出:指标ITER具有最大的基尼系数,其值为0.7017,这表明各个节点在该指标上的取值差别最大,该指标给节点的综合重要度评估提供了非常有用的信息,故该指标在节点的综合重要度评估中的比重应较大;指标IRST具有最小的基尼系数,其值为0.2677,这表明各个节点在该指标上的取值差别最小,该指标提供了较少的有用信息,故该指标所占的比重应较小。
此外,从表1还可以看出:指标IEBC具有最小的肯德尔系数,其值为0.3588,这表明该指标和其他指标的相关性最小,即该指标提供的有用信息和其他指标重合度不大,因此该指标在节点的综合重要度评估中的比重应较大;指标IRST具有最大的肯德尔系数,其值为0.6034,这表明该指标和其他指标的相关性较大,提供的有用信息重合度较大,因此该指标在节点的综合重要度评估中的比重应较小。从5个指标的定义可以看出,指标IEBC主要是从电网的运行状况来衡量节点的重要性程度,而其他指标主要是从拓扑结构和电气参数的角度来衡量节点的重要性程度,因此指标IEBC与其他四个指标的相关性不大,这与计算得到的该指标的肯德尔相关系数是相吻合的。综合基尼系数和肯德尔系数后,可以得到5个指标的综合客观权重分别为:0.1555,0.3316,0.2641,0.1608和0.0880。从中可以看出:经过综合后指标IEBC具有最大的权重;而指标IRST具有最小的权重。
接着,计算区域电网节点重要度评估问题的理想解和反理想解,进而得到各节点分别与理想解和反理想解的欧氏距离,在此基础上分别计算各个节点与理想解的相对逼近度,其结果如图2和表2所示。从图3和表2可以看出:对于发电厂节点,节点1在所有发电厂节点中排序最高且也位于节点重要度排名前十的序列中,因为节点1连接有抽水蓄能电厂,其有功发电容量为2400MW,且与2个变电站相连,因此其在整个电网中处于非常重要的地位;而对于变电站节点,节点29排在最重要的位置,因为其带有310MW的有功负荷且比较靠近发电厂节点。
此外,为了验证本发明方法的有效性和精确性,表2给出了采用抗灾变性指标(即针对输电网的特性,这里采用了考虑输电网架结构限制和电气极限容量限制后的系统失负荷电量指标)、Gini法、熵权法评估得到的节点重要度结果。
表2基于CRITIC和Gini法的节点重要度评估结果比较
Figure BDA0001382173300000141
从表2可以看出:采用Gini法、熵权法和本发明方法得到的节点重要度排名前十的节点编号集合是一致的,只是排序稍微有所区别;而采用输电网的抗灾变性指标进行评估得到的结果中有50%节点不在其他3种方法得到的评估结果中,故其得到结果差异较大,这是因为其只采用了单一指标,无法较为全面地考虑影响节点重要度评估的各方面因素。基于CRITIC法得到的最重要的变电站节点为节点29,而基于Gini法和熵权法得到的最重要的变电站节点为节点32。从该电网的实际情况看,由于节点24和29分别连接有125MW和310MW的有功负荷,且这2个节点附近的发电厂节点9和11分别接有400MW和60MW的有功发电容量。因此,从电网拓扑和电气特性来看,节点29的重要性应该较节点24重要些,即节点29应该排在节点24前面。
从表2可以看出:基于CRITIC法的结果更加符合分析的结果。此外,基于Gini法和熵权法得到的节点重要性结果一样,这是因为这2种方法都只考虑了评价指标的对比强度,没有考虑评价指标之间的冲突性;而本发明方法同时考虑评价指标的对比强度和评价指标之间的冲突性,因此从理论上也较能精确地评估电网节点的重要性。综上所述,相对而言基于CRITIC法的节点重要度评估结果更加符合电网实际运行的状况,即能够更好、更精确地识别出电网节点的重要性程度。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。

Claims (9)

1.一种区域电网关键节点识别方法,其特征在于,包括以下步骤:
输入电力系统拓扑结构参数、电气特性参数;
基于区域电网的复杂网络拓扑特性和电气特性,提出了评估区域电网关键节点的指标,即电力网络节点的功率集中度、电力网络节点的电气介数、电力网络的传输效能变化率、电力网络的凝聚度变化率、电力网络的生成树变化率;
计算各节点在各项评价指标下的指标值,生成标准化区域电网的节点重要度评估决策矩阵;
采用理想点法来解决区域电网的节点重要度评估问题,按照各个节点与理想解的相对逼近度大小,对区域电网的节点进行重要性程度排序,从而得到节点重要度的评估结果;
采用CRITIC法确定理想点法模型中指标的客观权重,利用基尼系数和肯德尔相关系数来分别衡量不同评价对象的取值差异性和评价指标之间的冲突性。
2.根据权利要求1所述的一种区域电网关键节点识别方法,其特征在于:采用电力网络节点的功率集中度来识别电网关键节点,包括:
在电力网络中,节点的度在拓扑结构方面定义为与该节点连接的输电线路数目;然而对于电力网络而言,节点的重要性程度不仅只与电网的拓扑结构有关,而且还与该节点所连接的重要负荷量或者发电机组的装机容量密切相关;如果某一个节点上所连接的发电机组装机容量越大或重要负荷量越大,则该节点的重要性程度就越大;故综合电力网络的复杂网络拓扑特性及其电气特性,定义了电力网络节点的功率集中度指标:
Figure FDA0002575607700000011
式中:Aij为网络邻接矩阵A的第i行第j列元素;Ω(i)为与节点i直接相连 的节点的集合;N为电力网络的节点总数;PG,max为发电机组的最大有功发电出力;PG,i和PL,i分别为节点i上所连接的发电机有功出力和重要有功负荷量;综上,如果某一个节点所连接的有功负荷和发电机组有功出力越大,以及与该节点连接的输电线路数目越多,则该节点的重要性程度就越高。
3.根据权利要求2所述的一种区域电网关键节点识别方法,其特征在于:采用电力网络节点的电气介数来识别电网关键节点,包括:
节点的介数定义为网络中所有节点对的最短路径经过该节点的条数;节点的介数越大,表明该节点在整个网络中起到枢纽作用,因此该节点的重要性程度也越高;然而,在电力网络中,电能会通过电源节点和负荷节点之间所有可能的路径进行传输,而不是复杂网络理论中所定义的最短路径;因此,复杂网络理论中所定义的节点介数不能直接用于分析电力网络;在电力网络中,功率传输分布系数用于表征电能在电力系统中的传输特性,该系数不仅考虑了电力网络的拓扑特性,而且还计及了输电线路的电气参数;因此,这里基于功率传输分布系数定义了节点电气介数:
Figure FDA0002575607700000021
Figure FDA0002575607700000022
Figure FDA0002575607700000023
式中:ΩSE和ΩSK分别表示电源节点和负荷节点的集合;
Figure FDA0002575607700000024
表示从电源节点g通过线路i-j传输到负荷节点l的有功功率;xij表示线路i-j的电抗;
Figure FDA0002575607700000025
表示节点i和g之间的等效阻抗;Zig表示从节点i到节点g的转移阻抗;Zii和Zgg分别表示节点i和g之间的输入阻抗,电力网络中节点的电气介数越大,表明该节点的重要性程度越高。
4.根据权利要求3所述的一种区域电网关键节点识别方法,其特征在于:采用电力网络传输效能的变化率来识别电网关键节点,包括:
在电力网络中,电力网络传输效能定义为电力网络中节点间最短电气距离 倒数的平均值,其表达式为:
Figure FDA0002575607700000031
式中:ΩN为电力网络的节点集合;
Figure FDA0002575607700000032
表示电力网络中节点j和k之间的最短电气距离,即节点j和k之间电能传输路径上输电线路的电抗值之和最小;
电力网络的连通性受到每一个节点的影响,电力网络如果失去部分关键节点将使其连通性以及电能传输效能大大降低;因此,电力网络中节点的重要性程度可以通过移除该节点及其所连接的输电线路前后电力网络传输效能的变化率来衡量;因此,节点i的电力网络传输效能可以定义为该节点移除前后电力网络传输效能的变化率,即
Figure FDA0002575607700000033
Figure FDA0002575607700000034
式中:EB,i表示电力网络在节点i失效后的电力网络传输效能;因此,一个节点的电力网络传输效能变化率越大,表明该节点失效后对电网的电能传输影响越大,则该节点在电网中的重要性程度也就越高。
5.根据权利要求4所述的一种区域电网关键节点识别方法,其特征在于:采用电力网络的凝聚度变化率来识别电网关键节点,包括:
在电力网络中,节点的凝聚度定义为电力网络中该节点到达其他所有节点的最短电气距离之和的倒数,该指标反映了该节点在电力网络中位置;因而,电力网络的凝聚度可以定义为所有节点之间最短电气距离之和的倒数,其表达式为:
Figure FDA0002575607700000035
在电力网络中,一个节点的失效对电力网络的凝聚度将产生影响;节点收缩法常被用于衡量节点失效对与其直接相连的邻居节点的影响;假设原始电力网络G中节点i收缩后得到了新的电力网络Gi,则电力网络Gi的凝聚度为
Figure FDA0002575607700000041
式中:NC为电力网络Gi的节点数;节点i的重要性程度可以通过节点i及其邻居节点收缩成一个节点后,新电力网络Gi的凝聚度变化率来衡量,故电力网络的凝聚度变化率可表示为:
Figure FDA0002575607700000042
该指标值越大,表明新电力网络Gi的凝聚度增加得越多,则该节点的重要性程度就越大。
6.根据权利要求5所述的一种区域电网关键节点识别方法,其特征在于:采用电力网络的生成树变化率来识别电网关键节点,在复杂网络中,电力网络的生成树为包括了所有节点的子图;在电力网络中,一个节点的失效使得电力网络的生成树数目减小,甚至可能使得电力网络成为一个不连通的网络;因此,该节点的重要性程度可以通过移除该节点及其直接相连的输电线路前后的生成树数目的变化率来衡量,其计算步骤为:
步骤1:定义电力网络的拉普拉斯矩阵为
Figure FDA0002575607700000043
式中:D为电力网络的度矩阵,即对角元素为
Figure FDA0002575607700000044
的对角矩阵;
步骤2:根据Kirchhoff矩阵树理论可得,电力网络的生成树数目等于拉普拉斯矩阵的任何一个N-1阶主子式,即
Figure FDA0002575607700000045
式中:λ12,...,λN-1为拉普拉斯矩阵的非零特征值;N为电力网络节点总数;
步骤3:节点i的生成树数目的变化率可表示为:
Figure FDA0002575607700000051
式中:τ(G-Vi)为电力网络删除第i节点及与其相连的输电线路后的生成树数目;τ(G-Vi)的数值越小,IRST(i)的数值越大,则节点i的失效对整个电力网络的影响程度越大,该节点越重要。
7.根据权利要求6所述的一种区域电网关键节点识别方法,其特征在于:包括建立多属性决策问题的决策矩阵;在区域电网的节点重要度评估中,假设该区域电网共有N个节点并选取了M个用于衡量节点重要度的评价指标;因此,
区域电网的节点重要度评估问题的决策矩阵可以表示为
Figure FDA0002575607700000052
式中:Iij表示节点i对应的指标j的指标值,其中i∈{1,2,...,N},j∈{1,2,...,M};M等于5。
8.根据权利要求7所述的一种区域电网关键节点识别方法,其特征在于:采用理想点法来解决区域电网的节点重要度评估问题,决策流程如下:
步骤1:生成标准化区域电网的节点重要度评估决策矩阵
Y=(yij)N×M
式中:
Figure FDA0002575607700000053
i∈{1,2,...,N},j∈{1,2,...,M};
步骤2:生成加权后的区域电网节点重要度评估矩阵
Z=(zij)N×M
式中:zij=ωjyij;ωj为指标j的权重,即指标权重向量ω=(ω12,L,ωM)的第j个元素,其中ωj∈[0,1],且
Figure FDA0002575607700000054
i∈{1,2,...,N},j∈{1,2,...,M};指标的权重采用CRITIC法确定;
步骤3:确定节点重要度评估问题的理想解
Figure FDA0002575607700000061
和反理想解
Figure FDA0002575607700000062
其中理想解是期望的最重要的电网节点,该节点的各个指标值都比电网中其它各个节点的相应指标的值都好,即
Figure FDA0002575607700000063
反理想解是假想的最不重要的电网节点,该节点的各个指标值都比电网中其它各个节点的相应指标的值都差,即
Figure FDA0002575607700000064
j∈{1,2,...,M};
步骤4:计算区域电网各节点分别与理想解和反理想解的欧氏距离
Figure FDA0002575607700000065
Figure FDA0002575607700000066
Figure FDA0002575607700000067
Figure FDA0002575607700000068
式中:
Figure FDA0002575607700000069
为加权后的节点重要度评估矩阵Z的第i行;
步骤5:计算区域电网中各个节点与理想解的相对逼近度
Figure FDA00025756077000000610
节点与理想解的相对逼近度的取值范围为[0,1],该值愈逼近1,则表明该节点越重要;反之,该值愈逼近0,则表明该节点越不重要;
步骤6:按照各个节点与理想解的相对逼近度大小,对区域电网的节点进行重要性程度排序,从而得到节点重要度的评估结果。
9.根据权利要求8所述的一种区域电网关键节点识别方法,其特征在于:采用CRITIC法来确定指标的权重;采用基尼系数和肯德尔系数分别来衡量指标在不同评价对象的取值差异性和评价指标之间的冲突性;包括:
第j个指标的基尼指数可以定义为:
Figure FDA00025756077000000611
式中:该值的取值范围为[0,1],1表示该指标的取值完全不相似,0表示该指 标的取值完全相似;在电网的节点重要度评估问题中,用于评估节点重要性的5个指标可以看作为基尼指数计算中的社会成员,而各个指标的属性值可以看作为基尼指数计算中社会成员的财富值;
具有N个元素的两列变量
Figure FDA0002575607700000071
Figure FDA0002575607700000072
其第i个变量值分别为zij和zik(1≤i≤N),其中j∈{1,2,...,M},k∈{1,2,...,M},;假设XjR和XkR分别为zij和zik
Figure FDA0002575607700000073
Figure FDA0002575607700000074
的排序值,则XjR和XkR对应排序值变量可以搭配组成变量对的集合XR,该变量对集合的第i个变量对为(XjR,XkR);因此,第j个和k个指标之间的肯德尔相关系数可以定义为:
Figure FDA0002575607700000075
式中:A为常数,其值为N(N-1)/2;Ncc和Ndc分别表示变量对的集合XR中变量排序值相等的变量对数目和变量值排序值不相等的变量对数目;
Figure FDA0002575607700000076
Figure FDA0002575607700000077
分别表示变量
Figure FDA0002575607700000078
Figure FDA0002575607700000079
中具有相同变量值的个数;因此,第j个指标与其他指标的整体肯德尔相关系数可以定义为:
Figure FDA00025756077000000710
当指标j肯德尔系数为1时,则表明该指标与其他指标具有一致的等级相关性;而肯德尔系数为0时,则表明该指标与其他指标是相互独立的;
综合考虑基尼系数和肯德尔系数,即指标j的客观权重可以表示为:
Figure FDA00025756077000000711
CN201710708797.4A 2017-08-17 2017-08-17 一种区域电网关键节点识别方法 Active CN107482626B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710708797.4A CN107482626B (zh) 2017-08-17 2017-08-17 一种区域电网关键节点识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710708797.4A CN107482626B (zh) 2017-08-17 2017-08-17 一种区域电网关键节点识别方法

Publications (2)

Publication Number Publication Date
CN107482626A CN107482626A (zh) 2017-12-15
CN107482626B true CN107482626B (zh) 2020-09-25

Family

ID=60601775

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710708797.4A Active CN107482626B (zh) 2017-08-17 2017-08-17 一种区域电网关键节点识别方法

Country Status (1)

Country Link
CN (1) CN107482626B (zh)

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3074916B1 (fr) * 2017-12-13 2019-12-27 Safran Electrical & Power Procede et systeme de surveillance d'un reseau electrique maille de retour de courant d'un aeronef en materiau composite
CN108281959B (zh) * 2018-01-03 2019-12-31 浙江大学 一种高比例可再生能源电力系统的骨干网架优化方法
CN108090696A (zh) * 2018-01-04 2018-05-29 浙江大学 考虑运行拓扑特性的电力信息物理融合系统作用关系排序方法
CN109389280B (zh) * 2018-08-17 2023-07-04 国网浙江省电力有限公司宁波供电公司 基于多属性决策模型的电能表生产厂商评价方法
CN109389145B (zh) * 2018-08-17 2023-10-10 国网浙江省电力有限公司宁波供电公司 基于计量大数据聚类模型的电能表生产厂商评价方法
CN109409629B (zh) * 2018-08-17 2023-07-25 国网浙江省电力有限公司宁波供电公司 基于多属性决策模型的采集终端生产厂商评价方法
CN109409628B (zh) * 2018-08-17 2023-07-04 国网浙江省电力有限公司 基于计量大数据聚类模型的采集终端生产厂商评价方法
CN110046797A (zh) * 2019-01-04 2019-07-23 国网浙江省电力有限公司 基于critic和理想点法的计量设备运行质量评估方法
CN110070256B (zh) * 2019-02-28 2023-12-08 国网浙江省电力有限公司 基于critic方法的零电量用户排查优先度权重计算方法
CN110188473A (zh) * 2019-05-30 2019-08-30 国网甘肃省电力公司经济技术研究院 一种电力骨干通信网节点重要度识别方法
CN110263791B (zh) * 2019-05-31 2021-11-09 北京京东智能城市大数据研究院 一种识别功能区的方法和装置
CN110336273A (zh) * 2019-06-10 2019-10-15 国网天津市电力公司 一种220千伏电网分区方法及装置
CN110350522B (zh) * 2019-07-10 2022-04-29 西南交通大学 一种基于加权h指数的电力系统脆弱线路辨识方法
CN110719191A (zh) * 2019-08-30 2020-01-21 中国电力科学研究院有限公司 一种面向次生失效的网络可靠性评估方法
CN110768260B (zh) * 2019-09-12 2022-09-06 南京邮电大学 一种基于电气介数的电网级联失效模型建立方法
CN110661341A (zh) * 2019-10-16 2020-01-07 佳源科技有限公司 一种分布式配网终端
CN111105145B (zh) * 2019-11-26 2022-07-26 广西电网有限责任公司 一种考虑间歇性能源的电网生存性评估方法
CN111339491A (zh) * 2020-02-28 2020-06-26 上海电力大学 一种城市配电网改造方案的评估方法
CN111725801B (zh) * 2020-05-06 2022-05-24 国家计算机网络与信息安全管理中心 基于注意力机制的配电系统脆弱节点辨识方法及系统
CN112348351B (zh) * 2020-11-05 2022-06-10 招商局重庆交通科研设计院有限公司 一种路网中桥梁抗震加固优先级评估方法
CN112508251A (zh) * 2020-11-27 2021-03-16 国网浙江省电力有限公司杭州供电公司 一种基于电力驾驶舱的负荷预测方法及装置
CN112487658B (zh) * 2020-12-14 2022-09-16 重庆邮电大学 一种电网关键节点的识别方法、装置及系统
CN113285828B (zh) * 2021-05-19 2022-04-29 湖南经研电力设计有限公司 复杂网络关键节点辨识方法及其电网关键节点辨识方法
CN113239250B (zh) * 2021-06-09 2022-09-02 河南科技大学 一种数据采集关键节点的筛选方法
CN113507389B (zh) * 2021-07-08 2022-05-03 华北电力大学 一种电网关键节点识别方法及系统
CN113569961B (zh) * 2021-07-30 2024-01-26 上海电机学院 一种电网节点分类方法及计算机可读介质
CN114143204B (zh) * 2021-11-19 2023-11-17 南方电网科学研究院有限责任公司 电网继电保护整定关键节点筛选方法、装置、终端及介质
CN114567562B (zh) * 2022-03-01 2024-02-06 重庆邮电大学 一种电网与通信网耦合网络关键节点识别的方法
CN116050701B (zh) * 2022-11-16 2023-12-15 华北电力大学 一种计及碳流的电力系统节点重要度评估方法
CN117729058B (zh) * 2024-02-18 2024-04-26 四川大学 一种应对网络攻击的信息物理系统关键节点辨识方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105721207A (zh) * 2016-01-29 2016-06-29 国家电网公司 电力通信网中通信节点重要度的确定方法和装置
CN106655174A (zh) * 2017-01-03 2017-05-10 昆明理工大学 一种配电网综合重构优化方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105721207A (zh) * 2016-01-29 2016-06-29 国家电网公司 电力通信网中通信节点重要度的确定方法和装置
CN106655174A (zh) * 2017-01-03 2017-05-10 昆明理工大学 一种配电网综合重构优化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Division algorithm and interconnection strategy of restoration subsystems based on complex network theory;Lin, Z.Z.等;《 IET Generation, Transmission and Distribution》;20111231;第5卷(第6期);674-683 *
基于信息熵权和层次分析法的电网关键节点识别;林鸿基等;《广东电力》;20161231;第29卷(第12期);50-56 *

Also Published As

Publication number Publication date
CN107482626A (zh) 2017-12-15

Similar Documents

Publication Publication Date Title
CN107482626B (zh) 一种区域电网关键节点识别方法
CN106850254B (zh) 一种电力通信网中关键节点识别方法
CN108038300B (zh) 基于改进的隶属度函数结合神经网络的光纤状态评估方法
CN110969347B (zh) 一种输电网结构形态评估方法
CN107169628B (zh) 一种基于大数据互信息属性约简的配电网可靠性评估方法
CN107508297B (zh) 一种配网拓扑结构校验与维护方法
CN112699913A (zh) 一种台区户变关系异常诊断方法及装置
CN112561251B (zh) 配电网异常点检测方法、装置、计算机设备及存储介质
CN112149873B (zh) 一种基于深度学习的低压台区线损合理区间预测方法
CN106329516A (zh) 一种基于典型场景识别的配电网动态重构方法
CN110417011A (zh) 一种基于互信息与迭代随机森林的在线动态安全评估方法
CN107437135B (zh) 一种新型储能选型方法
CN108336739A (zh) 一种基于rbf神经网络的概率潮流在线计算方法
CN106384302B (zh) 配电网可靠性评估方法和系统
CN111460374A (zh) 一种考虑节点差异性的配电网d-pmu优化配置方法
CN110601173B (zh) 基于边缘计算的配网拓扑识别方法及装置
CN107295537A (zh) 一种测评无线传感器网络可靠性的方法及系统
CN116148753A (zh) 一种智能电能表运行误差监测系统
CN109686402A (zh) 基于动态加权相互作用网络中关键蛋白质识别方法
Ren et al. A multiple randomized learning based ensemble model for power system dynamic security assessment
CN111027841A (zh) 一种基于梯度提升决策树的低压台区线损计算方法
CN112990776A (zh) 一种配网设备健康度评价方法
CN112989526B (zh) 一种基于核极限学习机的航空网络关键节点识别方法
CN108959652A (zh) 基于密度峰值与社区归属度的重叠社区发现方法
CN115081933B (zh) 一种基于改进谱聚类的低压用户拓扑构建方法和系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant