CN113507389B - 一种电网关键节点识别方法及系统 - Google Patents
一种电网关键节点识别方法及系统 Download PDFInfo
- Publication number
- CN113507389B CN113507389B CN202110772639.1A CN202110772639A CN113507389B CN 113507389 B CN113507389 B CN 113507389B CN 202110772639 A CN202110772639 A CN 202110772639A CN 113507389 B CN113507389 B CN 113507389B
- Authority
- CN
- China
- Prior art keywords
- node
- power
- relation matrix
- nodes
- influence relation
- 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 44
- 239000011159 matrix material Substances 0.000 claims abstract description 168
- 239000013598 vector Substances 0.000 claims abstract description 69
- 238000004364 calculation method Methods 0.000 claims abstract description 22
- 238000004088 simulation Methods 0.000 claims description 26
- 230000008859 change Effects 0.000 claims description 14
- 238000010276 construction Methods 0.000 claims description 14
- 230000003044 adaptive effect Effects 0.000 claims description 13
- 238000013507 mapping Methods 0.000 claims description 8
- 238000010606 normalization Methods 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 238000007796 conventional method Methods 0.000 abstract 1
- 230000007423 decrease Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
- H04L41/06—Management of faults, events, alarms or notifications
- H04L41/0631—Management of faults, events, alarms or notifications using root cause analysis; using analysis of correlation between notifications, alarms or events based on decision criteria, e.g. hierarchy, tree or time analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/28—Databases characterised by their database models, e.g. relational or object models
- G06F16/284—Relational databases
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L43/00—Arrangements for monitoring or testing data switching networks
- H04L43/08—Monitoring or testing based on specific metrics, e.g. QoS, energy consumption or environmental parameters
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L67/00—Network arrangements or protocols for supporting network services or applications
- H04L67/01—Protocols
- H04L67/12—Protocols specially adapted for proprietary or special-purpose networking environments, e.g. medical networks, sensor networks, networks in vehicles or remote metering networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B70/00—Technologies for an efficient end-user side electric power management and consumption
- Y02B70/30—Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S20/00—Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
- Y04S20/20—End-user application control systems
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Signal Processing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Theoretical Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Environmental & Geological Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种电网关键节点识别方法及系统,通过故障节点对其他节点的功率的改变来衡量节点故障后对其他节点的影响,构造基于潮流分布的影响关系矩阵,并通过节点实际功率值来比较各节点之间的差异性,构造基于电气拓扑结构的影响关系矩阵,进而根据基于潮流分布的影响关系矩阵、基于电气拓扑结构的影响关系矩阵和自适应加权系数,迭代计算电网节点重要度向量,最终根据电网节点重要度向量确定电网关键节点。本发明综合考虑了电网节点在拓扑结构和潮流分布方面的重要性,并通过迭代计算的方式,考虑系统变化后变量的更新,相较于现有的方法,能够更准确地识别出电网关键节点。
Description
技术领域
本发明涉及电网技术领域,特别是涉及一种电网关键节点识别方法及系统。
背景技术
随着能源互联网的不断发展,它对国家经济和安全起着更加重要的影响,电网结构日趋复杂。因此对电力系统的可靠性,安全性提出了更高的要求。电力系统各节点之间的网络连接性和网络功能存在个体差异,其中存在一些关键节点,这些关键节点失效容易导致严重的电力系统事故。因此,需要准确识别出这些关键节点,为电力系统安全稳定控制和风险评估提供有价值的参考。
已有的电网关键节点识别方法中,基于改进PageRank的关键节点识别方法使用静态下的功率进行分析,无法反应节点故障对电网的影响;基于电气介数的关键节点识别方法仅从电网拓扑结构进行分析,忽略了电网实际潮流分布对电网的影响;基于基尼系数的关键节点识别方法考虑电网中各节点的差异性来识别关键节点,但这种方法没有考虑到电网节点发生故障后,由于系统拓扑结构发生改变,节点之间的差异性也发生改变的影响,在实际测试中效果并不理想。
发明内容
本发明的目的是提供一种电网关键节点识别方法及系统,能够更加准确地识别出电网关键节点。
为实现上述目的,本发明提供了如下方案:
一种电网关键节点识别方法,所述方法包括:
根据节点故障前后对电网中其它节点所连接的线路上功率的变化量构造节点相关性矩阵;
根据所述节点相关性矩阵构造基于潮流分布的影响关系矩阵;
根据电网中每个节点的实际功率值,获得基于电气拓扑结构的影响关系矩阵;
根据所述基于潮流分布的影响关系矩阵、所述基于电气拓扑结构的影响关系矩阵和自适应加权系数,迭代计算电网节点重要度向量;
根据所述电网节点重要度向量确定电网关键节点。
可选的,根据所述节点相关性矩阵构造基于潮流分布的影响关系矩阵,具体包括:
根据所述节点相关性矩阵,利用公式构造初始的基于潮流分布的影响关系矩阵;其中,yij为节点i和节点j之间的相关性,cij为节点相关性矩阵中节点i对节点j的影响程度,为节点i失效后,节点j关联线路m上的相对功率变化量,为节点i失效后,节点j关联线路m上的功率变化量,为节点j关联线路m的剩余容量,ckj为节点k对节点j的影响程度,N为电网节点数;
将初始的基于潮流分布的影响关系矩阵按照列进行归一化,获得归一化后的基于潮流分布的影响关系矩阵。
可选的,所述根据电网中每个节点的实际功率值,获得基于电气拓扑结构的影响关系矩阵,具体包括:
根据每个节点的实际功率值,利用高斯映射函数确定每个节点的归一化功率值;其中,fi为节点i的归一化功率值,ui为节点i的实际功率值,umax为所有节点的实际功率值中的最大值,σ2为所有节点的实际功率值的方差;
根据每个节点的归一化功率值,利用公式Z=feT,确定初始的基于电气拓扑结构的影响关系矩阵;其中,Z为初始的基于电气拓扑结构的影响关系矩阵,f为所有节点的归一化功率值组成的向量,eT为全“1”行向量,eT的维数与电网节点数相同;
将初始的基于电气拓扑结构的影响关系矩阵按照列进行归一化,获得归一化后的基于电气拓扑结构的影响关系矩阵。
可选的,根据所述基于潮流分布的影响关系矩阵、所述基于电气拓扑结构的影响关系矩阵和自适应加权系数,迭代计算电网节点重要度向量,具体包括:
根据所述基于潮流分布的影响关系矩阵、所述基于电气拓扑结构的影响关系矩阵和自适应加权系数,利用公式进行迭代计算,直到第k+1次迭代计算得到的电网节点重要度向量R(k+1)中所有的元素与第k次迭代计算得到的电网节点重要度向量R(k)中对应元素的差值的绝对值均小于差值阈值时迭代结束,获得电网节点重要度向量;
可选的,根据所述电网节点重要度向量确定电网关键节点,具体包括:
对所述电网节点重要度向量中节点的重要度降序排序,获得重要度序列;
初始化节点序号q为1;
从电网中对重要度序列的第q个节点进行仿真移除,获得仿真移除第q个节点后电网的网络指标;
判断所述网络指标的变化率是否小于变化率阈值,获得判断结果;
若所述判断结果表示否,则将q的数值增加1,返回步骤“从电网中对重要度序列的第q个节点进行仿真移除,获得仿真移除第q个节点后电网的网络指标”;
若所述判断结果表示是,则将所述重要度序列中第一个节点到第q个节点之间的所有节点确定为电网关键节点。
一种电网关键节点识别系统,所述系统包括:
节点相关性矩阵构造模块,用于根据节点故障前后对电网中其它节点所连接的线路上功率的变化量构造节点相关性矩阵;
基于潮流分布的影响关系矩阵构造模块,用于根据所述节点相关性矩阵构造基于潮流分布的影响关系矩阵;
基于电气拓扑结构的影响关系矩阵构造模块,用于根据电网中每个节点的实际功率值,获得基于电气拓扑结构的影响关系矩阵;
电网节点重要度向量计算模块,用于根据所述基于潮流分布的影响关系矩阵、所述基于电气拓扑结构的影响关系矩阵和自适应加权系数,迭代计算电网节点重要度向量;
电网关键节点确定模块,用于根据所述电网节点重要度向量确定电网关键节点。
可选的,所述基于潮流分布的影响关系矩阵构造模块,具体包括:
初始的基于潮流分布的影响关系矩阵构造子模块,用于根据所述节点相关性矩阵,利用公式构造初始的基于潮流分布的影响关系矩阵;其中,yij为节点i和节点j之间的相关性,cij为节点相关性矩阵中节点i对节点j的影响程度,为节点i失效后,节点j关联线路m上的相对功率变化量,为节点i失效后,节点j关联线路m上的功率变化量,为节点j关联线路m的剩余容量,ckj为节点k对节点j的影响程度,N为电网节点数;
归一化后的基于潮流分布的影响关系矩阵获得子模块,用于将初始的基于潮流分布的影响关系矩阵按照列进行归一化,获得归一化后的基于潮流分布的影响关系矩阵。
可选的,所述基于电气拓扑结构的影响关系矩阵构造模块,具体包括:
归一化功率值确定子模块,用于根据每个节点的实际功率值,利用高斯映射函数确定每个节点的归一化功率值;其中,fi为节点i的归一化功率值,ui为节点i的实际功率值,umax为所有节点的实际功率值中的最大值,σ2为所有节点的实际功率值的方差;
初始的基于电气拓扑结构的影响关系矩阵确定子模块,用于根据每个节点的归一化功率值,利用公式Z=feT,确定初始的基于电气拓扑结构的影响关系矩阵;其中,Z为初始的基于电气拓扑结构的影响关系矩阵,f为所有节点的归一化功率值组成的向量,eT为全“1”行向量,eT的维数与电网节点数相同;
归一化后的基于电气拓扑结构的影响关系矩阵获得子模块,用于将初始的基于电气拓扑结构的影响关系矩阵按照列进行归一化,获得归一化后的基于电气拓扑结构的影响关系矩阵。
可选的,所述电网节点重要度向量计算模块,具体包括:
电网节点重要度向量确定子模块,用于根据所述基于潮流分布的影响关系矩阵、所述基于电气拓扑结构的影响关系矩阵和自适应加权系数,利用公式进行迭代计算,直到第k+1次迭代计算得到的电网节点重要度向量R(k+1)中所有的元素与第k次迭代计算得到的电网节点重要度向量R(k)中对应元素的差值的绝对值均小于差值阈值时迭代结束,获得电网节点重要度向量;
可选的,所述电网关键节点确定模块,具体包括:
重要度序列获得子模块,用于对所述电网节点重要度向量中节点的重要度降序排序,获得重要度序列;
初始化子模块,用于初始化节点序号q为1;
网络指标获得子模块,用于从电网中对重要度序列的第q个节点进行仿真移除,获得仿真移除第q个节点后电网的网络指标;
判断结果获得子模块,用于判断所述网络指标的变化率是否小于变化率阈值,获得判断结果;
步骤返回子模块,用于若所述判断结果表示否,则将q的数值增加1,返回步骤“从电网中对重要度序列的第q个节点进行仿真移除,获得仿真移除第q个节点后电网的网络指标”;
电网关键节点确定子模块,用于若所述判断结果表示是,则将所述重要度序列中第一个节点到第q个节点之间的所有节点确定为电网关键节点。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种电网关键节点识别方法及系统,通过故障节点对其他节点的功率的改变来衡量节点故障后对其他节点的影响,构造基于潮流分布的影响关系矩阵,并通过节点实际功率值来比较各节点之间的差异性,构造基于电气拓扑结构的影响关系矩阵,进而根据基于潮流分布的影响关系矩阵、基于电气拓扑结构的影响关系矩阵和自适应加权系数,迭代计算电网节点重要度向量,最终根据电网节点重要度向量确定电网关键节点。本发明综合考虑了电网节点在拓扑结构和潮流分布方面的重要性,并通过迭代计算的方式,考虑系统变化后变量的更新,相较于现有的方法,能够更准确地识别出电网关键节点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的一种电网关键节点识别方法的流程图;
图2为本发明实施例提供的网络拓扑图;
图3为本发明实施例提供的节点的重要度的示意图;
图4为本发明实施例提供的网络指标的对比结果示意图;图4(a)为剩余容量的对比结果示意图,图4(b)为剩余节点数的对比结果示意图,图4(c)为剩余负荷量的对比结果示意图,图4(d)为综合指标的对比结果示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种电网关键节点识别方法及系统,能够更加准确地识别出电网关键节点。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
一种电网关键节点识别方法,如图1所示,方法包括:
S101,根据节点故障前后对电网中其它节点所连接的线路上功率的变化量构造节点相关性矩阵;
S102,根据节点相关性矩阵构造基于潮流分布的影响关系矩阵;
S103,根据电网中每个节点的实际功率值,获得基于电气拓扑结构的影响关系矩阵;
S104,根据基于潮流分布的影响关系矩阵、基于电气拓扑结构的影响关系矩阵和自适应加权系数,迭代计算电网节点重要度向量;
S105,根据电网节点重要度向量确定电网关键节点。
具体过程如下:
步骤S101,当电网中的电力节点发生故障时,它所连接的所有线路都会失效,其他线路上功率会重新分配。根据电力节点故障后对其它节点所连接的线路上功率的变化量构造节点相关性矩阵C。
式中
其中,cij为节点相关性矩阵中节点i对节点j的影响程度,为节点i失效后,节点j关联线路m上的相对功率变化量,为节点i失效后,节点j关联线路m上的功率变化量,为节点j关联线路m的剩余容量。矩阵C表示节点之间的相互影响。当节点i失效后,节点j关联线路的相对功率变化量越大,则节点i对节点j的影响就越大。
步骤S102,考虑到如果节点对其它节点的影响较大,同时其自身也容易受到其它节点的影响,那么该节点对电网的影响较大,其重要性较大。因此需要同时考虑电网节点对其他节点的影响以及其他节点对该电网节点的影响,所以构造基于潮流分布的影响关系矩阵Y(yij),具体包括:
步骤S103,电力系统正常运行时,与节点i相连线路的功率值的绝对值之和称为节点i的实际功率值。为了减少小功率节点之间的功率值差距,增加大功率节点之间的功率值差距,本发明利用高斯映射函数,将每个节点的实际功率值映射为节点归一化功率值,该值的区间为[0.1,1],所述高斯映射函数为:
节点归一化功率值的大小可以反映此节点对系统中所有节点的潜在影响,节点之间的基于电气拓扑结构的影响关系矩阵Z为:
Z=feT
总结为:
根据每个节点的实际功率值,利用高斯映射函数确定每个节点的归一化功率值;其中,fi为节点i的归一化功率值,ui为节点i的实际功率值,umax为所有节点的实际功率值中的最大值,σ2为所有节点的实际功率值的方差;
根据每个节点的归一化功率值,利用公式Z=feT,确定初始的基于电气拓扑结构的影响关系矩阵;其中,Z为初始的基于电气拓扑结构的影响关系矩阵,f为所有节点的归一化功率值组成的向量,eT为全“1”行向量,eT的维数与电网节点数相同;
将初始的基于电气拓扑结构的影响关系矩阵按照列进行归一化,获得归一化后的基于电气拓扑结构的影响关系矩阵。
步骤S104,为了综合考虑节点在电气拓扑结构与潮流分布两方面的重要性,本发明用自适应加权系数α来为两种影响关系进行加权,得到最终的综合节点重要度值,即考虑电气拓扑结构与潮流分布的电网节点重要度。根据实际网络发电节点出力的差异性,来确定自适应系数α。差异性越大,自适应系数越大,综合节点重要度值主要由基于实际潮流的影响来决定;差异性越小,自适应系数越小,综合节点重要度值主要由基于电气拓扑和基于实际潮流的影响关系的平均值来决定。自适应系数α的计算公式为:
当发电机出力完全均匀时,自适应加权系数α为0.5,α∈[0.5,1]。
令N维列向量R为考虑电气拓扑结构与潮流分布的电网节点重要度向量,N为节点数。R的计算公式为:
初始向量R(0)所有元素设置为相同值1/N,N为节点数,当R(k+1)中所有的元素与上一步R(k)中对应元素的差值的绝对值都小于0.001时迭代结束,得到R。
步骤总结为:
根据基于潮流分布的影响关系矩阵、基于电气拓扑结构的影响关系矩阵和自适应加权系数,利用公式进行迭代计算,直到第k+1次迭代计算得到的电网节点重要度向量R(k+1)中所有的元素与第k次迭代计算得到的电网节点重要度向量R(k)中对应元素的差值的绝对值均小于差值阈值时迭代结束,获得电网节点重要度向量;
步骤S105,依据向量R所计算出的考虑电气拓扑结构与潮流分布的电网节点重要度识别电网中的关键节点。R(i)的值越大,表明节点i重要度越高,节点越关键,具体包括:
对电网节点重要度向量中节点的重要度降序排序,获得重要度序列;
初始化节点序号q为1;
从电网中对重要度序列的第q个节点进行仿真移除,获得仿真移除第q个节点后电网的网络指标;
判断网络指标的变化率是否小于变化率阈值,获得判断结果;
若判断结果表示否,则将q的数值增加1,返回步骤“从电网中对重要度序列的第q个节点进行仿真移除,获得仿真移除第q个节点后电网的网络指标”;
若判断结果表示是,则将重要度序列中第一个节点到第q个节点之间的所有节点确定为电网关键节点。
以下将根据具体实施例对本发明作进一步的描述,需要说明的是,本实施例以本技术方案为前提,给出了详细的实施方式和具体的操作过程,但本发明的保护范围并不限于本实施例。
利用本发明所述方法,计算每个节点的R值得到系统节点的重要度,如图3所示,横坐标为节点编号,纵坐标为节点重要度值。
为了验证本发明的优越性,将本发明所述方法所得到的关键节点排序(为了简便描述,后文用ETPD-CNIA表示本发明所述方法)同基于改进PageRank(MPG),基于电气介数(EB)与基于基尼系数(MABGC)的关键节点识别方法进行对比,排序结果如表1所示:
表1关键节点排序结果
排序 | MPG | EB | MABGC | ETPD-CNIA |
1 | 31 | 26 | 16 | 16 |
2 | 19 | 29 | 17 | 6 |
3 | 6 | 2 | 3 | 29 |
4 | 10 | 1 | 2 | 26 |
5 | 38 | 25 | 4 | 31 |
6 | 4 | 9 | 5 | 19 |
7 | 8 | 28 | 18 | 22 |
8 | 12 | 23 | 6 | 21 |
9 | 15 | 8 | 25 | 20 |
10 | 20 | 3 | 26 | 10 |
然后本发明根据不同方法下关键节点的排序结果,依此移除节点,并观察移除节点对网络性能指标的影响,影响越大,说明排序结果越合理。本发明采用的网络指标包括剩余容量,剩余节点数,剩余负荷量和综合指标共四个指标来验证本发明方法的有效性和正确性。因系统级联失效模型中线路中断的不确定性,对这四种攻击关键节点的方式进行仿真计算50次取平均值。得到的这四种指标的比较结果平均值如图4所示。
其中,图4的横坐标均表示攻击节点数,纵坐标分别表示:图4(a)剩余容量CR;图4(b)剩余节点数NR;图4(c)剩余负荷量LR;图4(d)综合指标CI。
各图说明了系统按关键节点排序依此失去节点后的电网脆弱性。本发明方法通过节点之间的相互影响关系来判别关键节点,节点对其他节点影响大,且节点自身的功率容量大,节点就重要,同时考虑了拓扑和潮流因素,所以排名越前面的节点,对整个系统的稳定性影响越大。图中明显看出,按本发明方法得到的关键节点排序去攻击节点,相比于对比方法更具有破坏性。
观察图(a),看到(a)中剩余容量指标没有单调下降,折线有波动性,因为逐渐攻击节点,系统中存活节点越来越少。存在这种情况:负荷节点损失很多,发电机的发电量随着负荷量的减少而减少,致使线路中流经的功率变小,但是电力网络线路容量不变,剩余容量会增大。所以个别节点会出现剩余容量增大的情况。对比发现,本发明方法能够准确找到电网关键节点,最终使得电网剩余容量下降最多。
攻击某些节点后,系统发生级联故障,系统解列为几部分,有的部分没有发电机或负荷节点,相当于整个区域失效。再次攻击节点时,节点可能在上一轮攻击中已经失效,因此出现仿真结果图中有出现下降极平缓的情况。观察图4(b)(c)(d)可以发现,按本发明方法进行攻击,指标较少出现平滑的走势,下降趋势比较剧烈,说明本方法能够准确地找到关键节点并进行攻击,使得电网性能下降更快,攻击更具有破坏性。
通过以上分析,本发明算法能够对电网关键节点进行有效识别。
本发明通过故障节点对其他节点的功率的改变来衡量节点故障后对其他节点的影响,通过节点实际功率值来比较各节点之间的差异性,综合考虑了电网节点在拓扑结构和潮流分布方面的重要性。并通过迭代计算的方式,考虑系统变化后变量的更新。相较于现有的方法,本方法能够更准确地识别出电网关键节点,并且在测试中有更好的效果。
本发明还提供了一种电网关键节点识别系统,系统包括:
节点相关性矩阵构造模块,用于根据节点故障前后对电网中其它节点所连接的线路上功率的变化量构造节点相关性矩阵;
基于潮流分布的影响关系矩阵构造模块,用于根据节点相关性矩阵构造基于潮流分布的影响关系矩阵;
基于电气拓扑结构的影响关系矩阵构造模块,用于根据电网中每个节点的实际功率值,获得基于电气拓扑结构的影响关系矩阵;
电网节点重要度向量计算模块,用于根据基于潮流分布的影响关系矩阵、基于电气拓扑结构的影响关系矩阵和自适应加权系数,迭代计算电网节点重要度向量;
电网关键节点确定模块,用于根据电网节点重要度向量确定电网关键节点。
基于潮流分布的影响关系矩阵构造模块,具体包括:
初始的基于潮流分布的影响关系矩阵构造子模块,用于根据节点相关性矩阵,利用公式构造初始的基于潮流分布的影响关系矩阵;其中,yij为节点i和节点j之间的相关性,cij为节点相关性矩阵中节点i对节点j的影响程度,为节点i失效后,节点j关联线路m上的相对功率变化量,为节点i失效后,节点j关联线路m上的功率变化量,为节点j关联线路m的剩余容量,ckj为节点k对节点j的影响程度,N为电网节点数;
归一化后的基于潮流分布的影响关系矩阵获得子模块,用于将初始的基于潮流分布的影响关系矩阵按照列进行归一化,获得归一化后的基于潮流分布的影响关系矩阵。
基于电气拓扑结构的影响关系矩阵构造模块,具体包括:
归一化功率值确定子模块,用于根据每个节点的实际功率值,利用高斯映射函数确定每个节点的归一化功率值;其中,fi为节点i的归一化功率值,ui为节点i的实际功率值,umax为所有节点的实际功率值中的最大值,σ2为所有节点的实际功率值的方差;
初始的基于电气拓扑结构的影响关系矩阵确定子模块,用于根据每个节点的归一化功率值,利用公式Z=feT,确定初始的基于电气拓扑结构的影响关系矩阵;其中,Z为初始的基于电气拓扑结构的影响关系矩阵,f为所有节点的归一化功率值组成的向量,eT为全“1”行向量,eT的维数与电网节点数相同;
归一化后的基于电气拓扑结构的影响关系矩阵获得子模块,用于将初始的基于电气拓扑结构的影响关系矩阵按照列进行归一化,获得归一化后的基于电气拓扑结构的影响关系矩阵。
电网节点重要度向量计算模块,具体包括:
电网节点重要度向量确定子模块,用于根据基于潮流分布的影响关系矩阵、基于电气拓扑结构的影响关系矩阵和自适应加权系数,利用公式进行迭代计算,直到第k+1次迭代计算得到的电网节点重要度向量R(k+1)中所有的元素与第k次迭代计算得到的电网节点重要度向量R(k)中对应元素的差值的绝对值均小于差值阈值时迭代结束,获得电网节点重要度向量;
电网关键节点确定模块,具体包括:
重要度序列获得子模块,用于对电网节点重要度向量中节点的重要度降序排序,获得重要度序列;
初始化子模块,用于初始化节点序号q为1;
网络指标获得子模块,用于从电网中对重要度序列的第q个节点进行仿真移除,获得仿真移除第q个节点后电网的网络指标;
判断结果获得子模块,用于判断网络指标的变化率是否小于变化率阈值,获得判断结果;
步骤返回子模块,用于若判断结果表示否,则将q的数值增加1,返回步骤“从电网中对重要度序列的第q个节点进行仿真移除,获得仿真移除第q个节点后电网的网络指标”;
电网关键节点确定子模块,用于若判断结果表示是,则将重要度序列中第一个节点到第q个节点之间的所有节点确定为电网关键节点。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (6)
1.一种电网关键节点识别方法,其特征在于,所述方法包括:
根据节点故障前后对电网中其它节点所连接的线路上功率的变化量构造节点相关性矩阵;
根据所述节点相关性矩阵构造基于潮流分布的影响关系矩阵;
根据电网中每个节点的实际功率值,获得基于电气拓扑结构的影响关系矩阵;
根据所述基于潮流分布的影响关系矩阵、所述基于电气拓扑结构的影响关系矩阵和自适应加权系数,迭代计算电网节点重要度向量;
根据所述电网节点重要度向量确定电网关键节点;
根据所述节点相关性矩阵构造基于潮流分布的影响关系矩阵,具体包括:
根据所述节点相关性矩阵,利用公式构造初始的基于潮流分布的影响关系矩阵;其中,yij为节点i和节点j之间的相关性,cij为节点相关性矩阵中节点i对节点j的影响程度, 为节点i失效后,节点j关联线路m上的相对功率变化量,为节点i失效后,节点j关联线路m上的功率变化量,rjm为节点j关联线路m的剩余容量,ckj为节点k对节点j的影响程度,N为电网节点数;
将初始的基于潮流分布的影响关系矩阵按照列进行归一化,获得归一化后的基于潮流分布的影响关系矩阵;
所述根据电网中每个节点的实际功率值,获得基于电气拓扑结构的影响关系矩阵,具体包括:
根据每个节点的实际功率值,利用高斯映射函数确定每个节点的归一化功率值;其中,fi为节点i的归一化功率值,ui为节点i的实际功率值,umax为所有节点的实际功率值中的最大值,σ2为所有节点的实际功率值的方差;
根据每个节点的归一化功率值,利用公式Z=feT,确定初始的基于电气拓扑结构的影响关系矩阵;其中,Z为初始的基于电气拓扑结构的影响关系矩阵,f为所有节点的归一化功率值组成的向量,eT为全“1”行向量,eT的维数与电网节点数相同;
将初始的基于电气拓扑结构的影响关系矩阵按照列进行归一化,获得归一化后的基于电气拓扑结构的影响关系矩阵。
2.根据权利要求1所述的电网关键节点识别方法,其特征在于,根据所述基于潮流分布的影响关系矩阵、所述基于电气拓扑结构的影响关系矩阵和自适应加权系数,迭代计算电网节点重要度向量,具体包括:
根据所述基于潮流分布的影响关系矩阵、所述基于电气拓扑结构的影响关系矩阵和自适应加权系数,利用公式进行迭代计算,直到第k+1次迭代计算得到的电网节点重要度向量R(k+1)中所有的元素与第k次迭代计算得到的电网节点重要度向量R(k)中对应元素的差值的绝对值均小于差值阈值时迭代结束,获得电网节点重要度向量;
3.根据权利要求1所述的电网关键节点识别方法,其特征在于,根据所述电网节点重要度向量确定电网关键节点,具体包括:
对所述电网节点重要度向量中节点的重要度降序排序,获得重要度序列;
初始化节点序号q为1;
从电网中对重要度序列的第q个节点进行仿真移除,获得仿真移除第q个节点后电网的网络指标;
判断所述网络指标的变化率是否小于变化率阈值,获得判断结果;
若所述判断结果表示否,则将q的数值增加1,返回步骤“从电网中对重要度序列的第q个节点进行仿真移除,获得仿真移除第q个节点后电网的网络指标”;
若所述判断结果表示是,则将所述重要度序列中第一个节点到第q个节点之间的所有节点确定为电网关键节点。
4.一种电网关键节点识别系统,其特征在于,所述系统包括:
节点相关性矩阵构造模块,用于根据节点故障前后对电网中其它节点所连接的线路上功率的变化量构造节点相关性矩阵;
基于潮流分布的影响关系矩阵构造模块,用于根据所述节点相关性矩阵构造基于潮流分布的影响关系矩阵;
基于电气拓扑结构的影响关系矩阵构造模块,用于根据电网中每个节点的实际功率值,获得基于电气拓扑结构的影响关系矩阵;
电网节点重要度向量计算模块,用于根据所述基于潮流分布的影响关系矩阵、所述基于电气拓扑结构的影响关系矩阵和自适应加权系数,迭代计算电网节点重要度向量;
电网关键节点确定模块,用于根据所述电网节点重要度向量确定电网关键节点;
所述基于潮流分布的影响关系矩阵构造模块,具体包括:
初始的基于潮流分布的影响关系矩阵构造子模块,用于根据所述节点相关性矩阵,利用公式构造初始的基于潮流分布的影响关系矩阵;其中,yij为节点i和节点j之间的相关性,cij为节点相关性矩阵中节点i对节点j的影响程度, 为节点i失效后,节点j关联线路m上的相对功率变化量,为节点i失效后,节点j关联线路m上的功率变化量,rjm为节点j关联线路m的剩余容量,ckj为节点k对节点j的影响程度,N为电网节点数;
归一化后的基于潮流分布的影响关系矩阵获得子模块,用于将初始的基于潮流分布的影响关系矩阵按照列进行归一化,获得归一化后的基于潮流分布的影响关系矩阵;
所述基于电气拓扑结构的影响关系矩阵构造模块,具体包括:
归一化功率值确定子模块,用于根据每个节点的实际功率值,利用高斯映射函数确定每个节点的归一化功率值;其中,fi为节点i的归一化功率值,ui为节点i的实际功率值,umax为所有节点的实际功率值中的最大值,σ2为所有节点的实际功率值的方差;
初始的基于电气拓扑结构的影响关系矩阵确定子模块,用于根据每个节点的归一化功率值,利用公式Z=feT,确定初始的基于电气拓扑结构的影响关系矩阵;其中,Z为初始的基于电气拓扑结构的影响关系矩阵,f为所有节点的归一化功率值组成的向量,eT为全“1”行向量,eT的维数与电网节点数相同;
归一化后的基于电气拓扑结构的影响关系矩阵获得子模块,用于将初始的基于电气拓扑结构的影响关系矩阵按照列进行归一化,获得归一化后的基于电气拓扑结构的影响关系矩阵。
5.根据权利要求4所述的电网关键节点识别系统,其特征在于,所述电网节点重要度向量计算模块,具体包括:
电网节点重要度向量确定子模块,用于根据所述基于潮流分布的影响关系矩阵、所述基于电气拓扑结构的影响关系矩阵和自适应加权系数,利用公式进行迭代计算,直到第k+1次迭代计算得到的电网节点重要度向量R(k+1)中所有的元素与第k次迭代计算得到的电网节点重要度向量R(k)中对应元素的差值的绝对值均小于差值阈值时迭代结束,获得电网节点重要度向量;
6.根据权利要求4所述的电网关键节点识别方法,其特征在于,所述电网关键节点确定模块,具体包括:
重要度序列获得子模块,用于对所述电网节点重要度向量中节点的重要度降序排序,获得重要度序列;
初始化子模块,用于初始化节点序号q为1;
网络指标获得子模块,用于从电网中对重要度序列的第q个节点进行仿真移除,获得仿真移除第q个节点后电网的网络指标;
判断结果获得子模块,用于判断所述网络指标的变化率是否小于变化率阈值,获得判断结果;
步骤返回子模块,用于若所述判断结果表示否,则将q的数值增加1,返回步骤“从电网中对重要度序列的第q个节点进行仿真移除,获得仿真移除第q个节点后电网的网络指标”;
电网关键节点确定子模块,用于若所述判断结果表示是,则将所述重要度序列中第一个节点到第q个节点之间的所有节点确定为电网关键节点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110772639.1A CN113507389B (zh) | 2021-07-08 | 2021-07-08 | 一种电网关键节点识别方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110772639.1A CN113507389B (zh) | 2021-07-08 | 2021-07-08 | 一种电网关键节点识别方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113507389A CN113507389A (zh) | 2021-10-15 |
CN113507389B true CN113507389B (zh) | 2022-05-03 |
Family
ID=78012314
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110772639.1A Active CN113507389B (zh) | 2021-07-08 | 2021-07-08 | 一种电网关键节点识别方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113507389B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113987724A (zh) * | 2021-10-22 | 2022-01-28 | 国网安徽省电力有限公司宿州供电公司 | 基于拓扑分析的电网风险辨识方法及系统 |
CN114266475A (zh) * | 2021-12-22 | 2022-04-01 | 深圳供电局有限公司 | 一种基于多属性决策的电力网络关键节点识别方法 |
CN115994136B (zh) * | 2023-01-10 | 2023-09-05 | 广东电网有限责任公司 | 一种基于能源网络拓扑关系的能源数据清洗方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104901308A (zh) * | 2015-06-26 | 2015-09-09 | 太原理工大学 | 电力系统关键线路辨识方法 |
CN107292481A (zh) * | 2017-04-26 | 2017-10-24 | 广东电网有限责任公司惠州供电局 | 一种基于节点重要度的电网关键节点评估方法 |
CN107947151A (zh) * | 2017-11-02 | 2018-04-20 | 广西电网有限责任公司电力科学研究院 | 一种电力系统关键节点的识别方法 |
CN112487658A (zh) * | 2020-12-14 | 2021-03-12 | 重庆邮电大学 | 一种电网关键节点的识别方法、装置及系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
BR112013005681B1 (pt) * | 2010-09-08 | 2020-02-04 | Siemens Ag | método de detecção de mudanças de topologia em redes elétricas e modem de comunicação de rede elétrica |
CA2915674A1 (fr) * | 2015-12-17 | 2017-06-17 | Francois Leonard | Mise a jour d'une topologie d'un reseau de distribution par reattribution successive des compteurs |
CN106529791B (zh) * | 2016-10-27 | 2019-07-19 | 云南电网有限责任公司 | 一种电力系统中支路重要度的评估方法 |
CN107482626B (zh) * | 2017-08-17 | 2020-09-25 | 广东电网有限责任公司惠州供电局 | 一种区域电网关键节点识别方法 |
CN110350524A (zh) * | 2019-07-11 | 2019-10-18 | 南京理工大学 | 一种基于节点重要度的直流潮流优化方法 |
-
2021
- 2021-07-08 CN CN202110772639.1A patent/CN113507389B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104901308A (zh) * | 2015-06-26 | 2015-09-09 | 太原理工大学 | 电力系统关键线路辨识方法 |
CN107292481A (zh) * | 2017-04-26 | 2017-10-24 | 广东电网有限责任公司惠州供电局 | 一种基于节点重要度的电网关键节点评估方法 |
CN107947151A (zh) * | 2017-11-02 | 2018-04-20 | 广西电网有限责任公司电力科学研究院 | 一种电力系统关键节点的识别方法 |
CN112487658A (zh) * | 2020-12-14 | 2021-03-12 | 重庆邮电大学 | 一种电网关键节点的识别方法、装置及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN113507389A (zh) | 2021-10-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113507389B (zh) | 一种电网关键节点识别方法及系统 | |
CN110542819B (zh) | 一种基于半监督dbnc的变压器故障类型诊断方法 | |
CN110705879A (zh) | 一种高比例可再生能源接入下的电网脆弱性评估方法 | |
CN113947016B (zh) | 针对电网紧急控制系统中深度强化学习模型的脆弱性评估方法 | |
CN110889111A (zh) | 一种基于深度置信网络的电网虚拟数据注入攻击的检测方法 | |
CN110034581B (zh) | 风电并网下电力系统的区间电气介数脆弱性评估方法 | |
CN109217296B (zh) | 一种基于加权潮流熵和电压稳定的电网关键节点辨识方法 | |
CN110969262A (zh) | 一种变压器故障诊断方法 | |
CN112217668B (zh) | 基于比较模型的自适应网络故障诊断方法 | |
CN112380762A (zh) | 一种基于vmd-woa-lssvm的输电线路短路故障诊断方法 | |
CN110634082A (zh) | 一种基于深度学习的低频减载系统运行阶段预测方法 | |
CN112821424A (zh) | 一种基于数据-模型融合驱动的电力系统频率响应分析方法 | |
CN113887729A (zh) | 一种低压电力线载波通信系统故障诊断方法 | |
CN117113244A (zh) | 基于自监督和混合专家模型的传感器数据异常检测方法 | |
CN117175548B (zh) | 一种基于随机潮流的连锁故障紧急控制方法 | |
CN114266396A (zh) | 一种基于电网特征智能筛选的暂态稳定判别方法 | |
CN111814394B (zh) | 一种基于相关性和冗余性检测的电力系统安全评估方法 | |
CN111900720B (zh) | 一种基于双层网页排序算法的输电网脆弱线路辨识方法 | |
CN116756939B (zh) | 一种基于anova-iwoa-svm混合模型的变压器故障诊断方法 | |
CN113472738A (zh) | 一种电力内网终端设备的攻击风险耦合建模方法 | |
CN117077342A (zh) | 一种基于复杂网络理论的电力通信网络拓扑优化方法 | |
CN116911653A (zh) | 一种基于博弈论的配电网综合脆弱性分析方法及系统 | |
CN111191867B (zh) | 电力系统复杂网络的可靠性评估方法 | |
CN116628488A (zh) | 风电功率预测模型的训练方法、风电功率预测方法及设备 | |
CN111313408B (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 |