CN107274110A - 在信息层网络影响下的电力网络脆弱性评估方法 - Google Patents

在信息层网络影响下的电力网络脆弱性评估方法 Download PDF

Info

Publication number
CN107274110A
CN107274110A CN201710546839.9A CN201710546839A CN107274110A CN 107274110 A CN107274110 A CN 107274110A CN 201710546839 A CN201710546839 A CN 201710546839A CN 107274110 A CN107274110 A CN 107274110A
Authority
CN
China
Prior art keywords
mrow
msub
node
munder
mfrac
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.)
Pending
Application number
CN201710546839.9A
Other languages
English (en)
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.)
Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd
Original Assignee
Electric Power Dispatch Control Center 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 Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd filed Critical Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd
Priority to CN201710546839.9A priority Critical patent/CN107274110A/zh
Publication of CN107274110A publication Critical patent/CN107274110A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种在信息层网络影响下的电力网络脆弱性评估方法,该方法是在传统的单独对电力网络脆弱性分析和传统的对信息层网络的脆弱性基础上,通过对信息层网络和电力网络的交互影响分析,来综合信息层网络影响下的电力网络脆弱度分析。本方法能够通过求出负荷系统脆弱性指标来评判某一时间断面电力‑信息层复合系统的脆弱性程度。能够有效地分析复合系统的薄弱部分,提高复合系统的安全性和稳定性。

Description

在信息层网络影响下的电力网络脆弱性评估方法
技术领域
本发明属于电力系统自动化技术领域,具体涉及一种考虑信息层网络影响下的电力系统的脆弱性方法。
背景技术
随着电力系统的不断发展,我国电网正往跨区域、大规模的电网发展。这同时也提高了对电网和信息层网络的安全要求,电力技术的不断发展,同时信息层和电力系统的相互影响变得更加多,则电力系统在保持自身的稳定性的同时,对信息层网络安全的要求越来越高,网络脆弱性在一定程度上能够测量网络可靠性,可反映不同网络单元失效时网络性能的下降程度。以便于找出网络中的薄弱环节。关于对电力网络脆弱性和信息层网络脆弱性的发明和研究已经有不少,但对于信息层系统影响下的电力网络脆弱性发明几乎没有,为了进一步研究复合系统的脆弱性,需要对传统的的研究方法进行筛选和改进,提出新的脆弱性评估方法。本发明就是基于电气介数求出电力网络脆弱性和基于复杂网络求信息层网络脆弱性及负荷节点间的信息通道的脆弱性基础上,求出电力-信息层负荷系统的脆弱性。
本发明以脆弱度为指标,在综合电力网络和信息层网络的基础上,提出复合系统的脆弱性指标。
发明内容
本发明的目的在于针对现有技术的不足,提供一种在信息层网络影响下的电力网络脆弱性评估方法,本发明考虑了信息层网络对电力系统的影响因素,提出一个良好有效的考虑信息层网络-电力系统复合网络评估方法,具体的是先通过考虑单独的电力系统和信息层网络脆弱性,再通过研究信息层网络-电力系统复合系统的脆弱性。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤S1:先对原始数据输入后的电网进行拓扑化建模,再对所有需要评估的线路,节点和相应的保护进行分组,以求出电力系统线路和节点的保护脆弱度。基于基尔霍夫定律得到线路和节点的电气介数。利用平均电压传输能力E、电压相对裕度H、无功变化率W这三个指标来构成互补性指标集,E指标从全局性角度评估了事故链发生后电网的输电能力,H指标和W指标衡量了事故链发生对系统局部无功平衡的影响,并对E、H指标求倒后,将全部指标归一化处理定义保护k的第i条事故链的后果函数Ik,j。由事故链的后果函数可定义保护k的风险保护度R(k)。并由此可定义线路保护脆弱度λ(i,j)和节点保护脆弱度λ(i)
步骤S2:结合电力系统连锁故障的扩大和继电保护装置的不正确动作,提出结合电气介数和单元保护脆弱度的评估模型。
步骤S3:基于复杂网络计算信息层网络的效能函数。定义任两个节点间的距离为连接这两个节点的最短路径dij计算系统的总效能E(G),考虑连通性改进效能函数。由业务矩阵Sij计算网络的总业务量S,再将最短路径矩阵Aij,其中aij是节点i和j之间的最短路径。通过统计每次仿真去除1个点时失去连接的业务数量,来计算改进后的效能函数E′,即考虑了业务联通性及业务分布对网络传输特性的影响。
步骤S4:考虑电力系统层和信息层网络之间的业务信息交互通道层的通道脆弱性,融合信息传输的路由选择概率后计算复合节点间信息通道的脆弱性。
步骤S5:综合考虑电力脆弱性指标、信息层脆弱性指标和信息交互层的通道脆弱性,可得到考虑信息层网络影响的电力系统脆弱性评估指标。在计算过程中电力脆弱性指标和信息层网络脆弱性评估指标和复合节点间信息通道的脆弱值均归一化。
所述的步骤S1中对线路lij的电气介数定义为
其中:G和L分别为发电机节点集合和负荷节点集合;(m,n)为发电机-负荷节点对;Imn(i,j)为在发电机-负荷节点对(m,n)间加上单位注入电流元后,在线路lij上引起的电流。
节点i的电气介数定义为
其中:(m,n)为发电-负荷节点对;Bmn(i)为(m,n)间加入单位注入电流源后在节点i上产生的电气介数。
其中j为所有与节点i有支路直接相连的节点。表征了线路、节点的结构脆弱性。
平均传输能力是衡量事故链发生后的最大连通子网的输电能力:
式中:M为负荷节点的集合;Pv为负荷节点吸收的有功功率;l为任意一条线路;Pl为线路l上的有功功率;Xl为线路l的电抗值。
相对裕度是衡量事故链发生后的负荷电压裕度水平。
式中:P为所有点的集合;Vcr为临界状态电压标幺值;Vs为事故链发生后的电压标幺值。
无功变化率为衡量事故链发生后的发电机无功出力的变化率。
Qg和Qg0分别表示事故链发生后和事故链发生前发电机g的无功出力,将E.H求倒后,所有指标归一化处理后定义保护k的第i条事故链的后果函数Ik,i为:
式中指标E求导后归一化的数值,指标H求导后归一化的数值,N(W)指标W归一化的数值,式中Ei为保护k的第i条事故链发生后的平均传输能力;Hi为保护k的第i条事故链发生后的电压相对裕度;Wi为保护k的第i条事故链发生后的无功变化率。
定义保护k的风险重要度
式中:Lk为保护k的事故链集合;PL(i)为保护k的第i条事故链的概率;Ik,i为保护k的第i条事故链的后果函数.若风险重要度越大,则保护k的不正确动作对电网安全影响越大。
线路保护脆弱度是线路lij上的保护对其脆弱性的平均贡献程度。
节点保护脆弱度是节点i的邻接线路上的保护对其脆弱性的平均贡献程度
式中:R(k)和R(k+1)为线路lij两端保护的风险重要度;j为所有与节点i有线路直接相连节点;Ki为节点i的度。
所述的步骤S2中,结合了电气介数与单元保护脆弱度的电网脆弱性评估模型:
Q(i,j)=B(i,j)λ(i,j) (11)
Q(i)=B(i)λ(i) (12)
式中:Q(i,j)为线路lij的脆弱度;Q(i)为节点i的脆弱度。
所述的步骤S3中,平均距离L,即反映了网络中节点之间信息传播的平均长度,即2个节点之间的最短距离:
n为节点总个数,定义整个网络G的功效性指标为:
假设εij表示2节点之间的通信效率,与这2节点之间的最短距离成反比,εij=(1/dij)。
改进后的效能函数为:
其中E是原来的功效性指标效能函数,为正在传输的业务所占的比例。可看出修改后的信息层网络的脆弱性评估指标更全面。
所述的步骤S5中,定义复合节点间信息通道的脆弱性时考虑了路由选择概率:
式中:Cf[Ef(i,j)]是通信节点i和电力节点j间的通脆弱性,psk是第k类业务信息在该通道的路由选择概率。εt为通道传输实际延时和固有时延之比;n为通道中的业务量,vk为第k类业务的传输速率;Bsk为第k类通信业务的重要度指标。
电力节点脆弱性和信息层节点脆弱值和信息交互通道的脆弱值三者叠加可以得到考虑信息层影响下的电力节点脆弱性指标:
Cp'(Npi)=Q(i)+[1-Q(i)]·maxj{max{Cf[Ef(i,j)]},E'} (17)
在式中计算时,各指标均需要进行归一化处理,其中maxj就是与电力i节点相连的信息层通道的最大脆弱值。
考虑到信息层的融入会对电力支路的脆弱性产生间接的影响,因此考虑了信息层网络的影响下的电力支路脆弱性指标:
式中wik表示节点i的度数。
本发明的有益效果是:本发明在传统的单独对电力网络脆弱性分析和传统的对信息层网络的脆弱性基础上,通过对信息层网络和电力网络的交互影响分析,来综合信息层网络影响下的电力网络脆弱度分析。本方法能够通过求出负荷系统脆弱性指标来评判某一时间断面电力-信息层复合系统的脆弱性程度。能够有效地分析复合系统的薄弱部分,提高复合系统的安全性和稳定性。
附图说明
图1是本发明流程图。
具体实施方式
下面将结合附图对本发明作进一步说明。
步骤S1:先对原始数据输入后的电网进行拓扑化建模,再对所有需要评估的线路,节点和相应的保护进行分组,以求出电力系统线路和节点的保护脆弱度。基于基尔霍夫定律得到线路和节点的电气介数。利用平均电压传输能力E、电压相对裕度H、无功变化率W这三个指标来构成互补性指标集,E指标从全局性角度评估了事故链发生后电网的输电能力,H指标和W指标衡量了事故链发生对系统局部无功平衡的影响,并对E、H指标求倒后,将全部指标归一化处理定义保护k的第i条事故链的后果函数Ik,j。由事故链的后果函数可定义保护k的风险保护度R(k)。并由此可定义线路保护脆弱度λ(i,j)和节点保护脆弱度λ(i);
步骤S2:结合电力系统连锁故障的扩大和继电保护装置的不正确动作,提出结合电气介数和单元保护脆弱度的评估模型。
步骤S3:基于复杂网络计算信息层网络的效能函数。定义任两个节点间的距离为连接这两个节点的最短路径dij计算系统的总效能E(G),考虑连通性改进效能函数。由业务矩阵Sij计算网络的总业务量S,再将最短路径矩阵Aij,其中aij是节点i和j之间的最短路径。通过统计每次仿真去除1个点时失去连接的业务数量,来计算改进后的效能函数E′,即考虑了业务联通性及业务分布对网络传输特性的影响。
步骤S4:考虑电力系统层和信息层网络之间的业务信息交互通道层的通道脆弱性,融合信息传输的路由选择概率后计算复合节点间信息通道的脆弱性。
步骤S5:综合考虑电力脆弱性指标、信息层脆弱性指标和信息交互层的通道脆弱性,可得到考虑信息层网络影响的电力系统脆弱性评估指标。在计算过程中电力脆弱性指标和信息层网络脆弱性评估指标和复合节点间信息通道的脆弱值均归一化。
所述的步骤S1中对线路lij的电气介数定义为
其中:G和L分别为发电机节点集合和负荷节点集合;(m,n)为发电机-负荷节点对;Imn(i,j)为在发电机-负荷节点对(m,n)间加上单位注入电流元后,在线路lij上引起的电流。
节点i的电气介数定义为
其中:(m,n)为发电-负荷节点对;Bmn(i)为(m,n)间加入单位注入电流源后在节点i上产生的电气介数。
其中j为所有与节点i有支路直接相连的节点。表征了线路、节点的结构脆弱性。
平均传输能力是衡量事故链发生后的最大连通子网的输电能力:
式中:M为负荷节点的集合;Pv为负荷节点吸收的有功功率;l为任意一条线路;Pl为线路l上的有功功率;Xl为线路l的电抗值。
相对裕度是衡量事故链发生后的负荷电压裕度水平。
式中:P为所有点的集合;Vcr为临界状态电压标幺值;Vs为事故链发生后的电压标幺值。
无功变化率为衡量事故链发生后的发电机无功出力的变化率。
Qg和Qg0分别表示事故链发生后和事故链发生前发电机g的无功出力,将E.H求倒后,所有指标归一化处理后定义保护k的第i条事故链的后果函数Ik,i为:
式中指标E求导后归一化的数值,指标H求导后归一化的数值,N(W)指标W归一化的数值,式中Ei为保护k的第i条事故链发生后的平均传输能力;Hi为保护k的第i条事故链发生后的电压相对裕度;Wi为保护k的第i条事故链发生后的无功变化率。
定义保护k的风险重要度
式中:Lk为保护k的事故链集合;PL(i)为保护k的第i条事故链的概率;Ik,i为保护k的第i条事故链的后果函数.若风险重要度越大,则保护k的不正确动作对电网安全影响越大。
线路保护脆弱度是线路lij上的保护对其脆弱性的平均贡献程度。
节点保护脆弱度是节点i的邻接线路上的保护对齐脆弱性的平均贡献程度
式中:R(k)和R(k+1)为线路lij两端保护的风险重要度;j为所有与节点i有线路直接相连节点;Ki为节点i的度。
所述的步骤S2中,结合了电气介数与单元保护脆弱度的电网脆弱性评估模型:
Q(i,j)=B(i,j)λ(i,j) (11)
Q(i)=B(i)λ(i) (12)
式中:Q(i,j)为线路lij的脆弱度;Q(i)为节点i的脆弱度。
所述的步骤S3中,平均距离L,即反映了网络中节点之间信息传播的平均长度,即2个节点之间的最短距离:
n为节点总个数,定义整个网络G的功效性指标为:
假设εij表示2节点之间的通信效率,与这2节点之间的最短距离成反比,εij=(1/dij)。
改进后的效能函数为:
其中E是原来的功效性指标效能函数,为正在传输的业务所占的比例。可看出修改后的信息层网络的脆弱性评估指标更全面。
所述的步骤S5中,定义复合节点间信息通道的脆弱性时考虑了路由选择概率:
式中:Cf[Ef(i,j)]是通信节点i和电力节点j间的通脆弱性,psk是第k类业务信息在该通道的路由选择概率。εt为通道传输实际延时和固有时延之比;n为通道中的业务量,vk为第k类业务的传输速率;Bsk为第k类通信业务的重要度指标。
电力节点脆弱性和信息层节点脆弱值和信息交互通道的脆弱值三者叠加可以得到考虑信息层影响下的电力节点脆弱性指标:
Cp'(Npi)=Q(i)+[1-Q(i)]·maxj{max{Cf[Ef(i,j)]},E'} (17)
在式中计算时,各指标均需要进行归一化处理,其中maxj就是与电力i节点相连的信息层通道的最大脆弱值。
考虑到信息层的融入会对电力支路的脆弱性产生间接的影响,因此考虑了信息层网络的影响下的电力支路脆弱性指标:
式中wik表示节点i的度数。
本发明的有益效果是:本发明在传统的单独对电力网络脆弱性分析和传统的对信息层网络的脆弱性基础上,通过对信息层网络和电力网络的交互影响分析,来综合信息层网络影响下的电力网络脆弱度分析。本方法能够通过求出负荷系统脆弱性指标来评判某一时间断面电力-信息层复合系统的脆弱性程度。能够有效地分析复合系统的薄弱部分,提高复合系统的安全性和稳定性。

Claims (5)

1.在信息层网络影响下的电力网络脆弱性评估方法,其特征在于,包括如下步骤:
步骤S1:对原始数据输入后的电网进行拓扑化建模,对所有需要评估的线路、节点和相应的保护进行分组,求出电力系统线路和节点的保护脆弱度:
基于基尔霍夫定律得到线路的电气介数B(i,j)和节点的电气介数B(i),利用平均电压传输能力E、电压相对裕度H、无功变化率W三个指标构成互补性指标集,对E、H指标求倒后,将全部指标归一化后,获得保护k的第i条事故链的后果函数Ik,j,再由事故链的后果函数获得保护k的风险保护度R(k),由此获得线路保护脆弱度λ(i,j)和节点保护脆弱度λ(i);
步骤S2:结合电力系统连锁故障的扩大和继电保护装置的不正确动作,建立结合电气介数和单元保护脆弱度的评估模型;
步骤S3:基于复杂网络计算信息层网络的效能函数,定义任两个节点间的距离为连接这两个节点的最短路径dij,计算系统的总效能E(G),并考虑连通性改进效能函数;
步骤S4:考虑电力系统层和信息层网络之间的业务信息交互通道层的通道脆弱性,融合信息传输的路由选择概率后计算复合节点间信息通道的脆弱性;
步骤S5:综合考虑电力脆弱性指标、信息层脆弱性指标和信息交互层的通道脆弱性,得到考虑信息层网络影响的电力系统脆弱性评估指标。
2.根据权利要求1所述的在信息层网络影响下的电力网络脆弱性评估方法,其特征在于,步骤S1中,对线路lij的电气介数定义为
<mrow> <mi>B</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>&amp;Element;</mo> <mi>G</mi> <mo>,</mo> <mi>n</mi> <mo>&amp;Element;</mo> <mi>L</mi> </mrow> </munder> <mo>|</mo> <msub> <mi>I</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>|</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
其中:G和L分别为发电机节点集合和负荷节点集合;(m,n)为发电机-负荷节点对;Imn(i,j)为在发电机-负荷节点对(m,n)间加上单位注入电流元后,在线路lij上引起的电流;
对节点i的电气介数定义为
<mrow> <mi>B</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>&amp;Element;</mo> <mi>G</mi> <mo>,</mo> <mi>n</mi> <mo>&amp;Element;</mo> <mi>L</mi> </mrow> </munder> <msub> <mi>B</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
其中:(m,n)为发电-负荷节点对;Bmn(i)为(m,n)间加入单位注入电流源后在节点i上产生的电气介数,
<mrow> <msub> <mi>B</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <munder> <mo>&amp;Sigma;</mo> <mi>j</mi> </munder> <mo>|</mo> <msub> <mi>I</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> </msub> <mo>|</mo> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>,</mo> <mi>i</mi> <mo>&amp;NotEqual;</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>1</mn> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
其中j为所有与节点i有支路直接相连的节点;
平均电压传输能力E为:
<mrow> <mi>E</mi> <mo>=</mo> <mo>(</mo> <mrow> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>v</mi> <mo>&amp;Element;</mo> <mi>M</mi> </mrow> </munder> <msub> <mi>P</mi> <mi>v</mi> </msub> </mrow> <mo>)</mo> <mo>/</mo> <mrow> <mo>(</mo> <munder> <mo>&amp;Sigma;</mo> <mi>l</mi> </munder> <msub> <mi>P</mi> <mi>l</mi> </msub> <msub> <mi>X</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
式中:M为负荷节点的集合;Pv为负荷节点吸收的有功功率;l为任意一条线路;Pl为线路l上的有功功率;Xl为线路l的电抗值;
电压相对裕度H为:
<mrow> <mi>H</mi> <mo>=</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>s</mi> <mo>&amp;Element;</mo> <mi>P</mi> </mrow> </munder> <mfrac> <mrow> <mo>|</mo> <msub> <mi>V</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>V</mi> <mrow> <mi>c</mi> <mi>r</mi> </mrow> </msub> <mo>|</mo> </mrow> <msub> <mi>V</mi> <mrow> <mi>c</mi> <mi>r</mi> </mrow> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow> 1
式中:P为所有点的集合;Vcr为临界状态电压标幺值;Vs为事故链发生后的电压标幺值;
无功变化率W为:
<mrow> <mi>W</mi> <mo>=</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>g</mi> <mo>&amp;Element;</mo> <mi>G</mi> </mrow> </munder> <mfrac> <mrow> <mo>|</mo> <msub> <mi>Q</mi> <mi>g</mi> </msub> <mo>-</mo> <msub> <mi>Q</mi> <mrow> <mi>g</mi> <mn>0</mn> </mrow> </msub> <mo>|</mo> </mrow> <msub> <mi>Q</mi> <mrow> <mi>g</mi> <mn>0</mn> </mrow> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
Qg和Qg0分别表示事故链发生后和事故链发生前发电机g的无功出力,保护k的第i条事故链的后果函数Ik,i为:
<mrow> <msub> <mi>I</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>3</mn> </mfrac> <mo>&amp;lsqb;</mo> <mi>N</mi> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <msub> <mi>E</mi> <mi>i</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <mi>N</mi> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <msub> <mi>H</mi> <mi>i</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <mi>N</mi> <mrow> <mo>(</mo> <msub> <mi>W</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
式中指标E求导后归一化的数值,指标H求导后归一化的数值,N(W)指标W归一化的数值,式中Ei为保护k的第i条事故链发生后的平均传输能力;Hi为保护k的第i条事故链发生后的电压相对裕度;Wi为保护k的第i条事故链发生后的无功变化率。
保护k的风险保护度
<mrow> <mi>R</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>L</mi> <mi>k</mi> </msub> </mrow> </munder> <msub> <mi>P</mi> <mi>L</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <msub> <mi>I</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
式中:Lk为保护k的事故链集合;PL(i)为保护k的第i条事故链的概率;Ik,i为保护k的第i条事故链的后果函数;
线路保护脆弱度为:
<mrow> <mi>&amp;lambda;</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mo>&amp;lsqb;</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
节点保护脆弱度为:
<mrow> <mi>&amp;lambda;</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>j</mi> </munder> <mi>&amp;lambda;</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> <msub> <mi>K</mi> <mi>i</mi> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
式中:R(k)和R(k+1)为线路lij两端保护的风险重要度;j为所有与节点i有线路直接相连节点;Ki为节点i的度。
3.根据权利要求1所述的在信息层网络影响下的电力网络脆弱性评估方法,其特征在于,步骤S2中,建立结合电气介数与单元保护脆弱度的电网脆弱性评估模型:
Q(i,j)=B(i,j)λ(i,j) (11)
Q(i)=B(i)λ(i) (12)
式中:Q(i,j)为线路lij的脆弱度;Q(i)为节点i的脆弱度。
4.根据权利要求1所述的在信息层网络影响下的电力网络脆弱性评估方法,其特征在于,步骤S3中,获得网络中节点之间信息传播的平均长度,即平均距离L:
<mrow> <mi>L</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>n</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mfrac> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>&amp;NotEqual;</mo> <mi>j</mi> </mrow> </munder> <msub> <mi>d</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
n为节点总个数,定义整个网络G的功效性指标为:
<mrow> <mi>E</mi> <mrow> <mo>(</mo> <mi>G</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>n</mi> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </mrow> </mfrac> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>&amp;Element;</mo> <mi>G</mi> <mo>,</mo> <mi>i</mi> <mo>&amp;NotEqual;</mo> <mi>j</mi> </mrow> </munder> <msub> <mi>&amp;epsiv;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>n</mi> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </mrow> </mfrac> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>&amp;Element;</mo> <mi>G</mi> <mo>,</mo> <mi>i</mi> <mo>&amp;NotEqual;</mo> <mi>j</mi> </mrow> </munder> <mfrac> <mn>1</mn> <msub> <mi>d</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow> 2
假设εij表示2节点之间的通信效率,与这2节点之间的最短距离成反比,εij=(1/dij)。
改进后的效能函数为:
<mrow> <msup> <mi>E</mi> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <mfrac> <mrow> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>&amp;NotEqual;</mo> <mi>j</mi> <mo>,</mo> <msub> <mi>a</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&amp;NotEqual;</mo> <mn>0</mn> </mrow> </munder> <msub> <mi>S</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mrow> <mrow> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>&amp;NotEqual;</mo> <mi>j</mi> </mrow> </munder> <msub> <mi>S</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mrow> </mfrac> <mi>E</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
其中E是原功效性指标效能函数,为正在传输的业务所占的比例。
5.根据权利要求1所述的在信息层网络影响下的电力网络脆弱性评估方法,其特征在于,步骤S5中,定义复合节点间信息通道的脆弱性时考虑了路由选择概率:
<mrow> <msub> <mi>C</mi> <mi>f</mi> </msub> <mo>&amp;lsqb;</mo> <msub> <mi>E</mi> <mi>f</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>=</mo> <msub> <mi>&amp;epsiv;</mi> <mi>t</mi> </msub> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mfrac> <msub> <mi>v</mi> <mi>k</mi> </msub> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>v</mi> <mi>j</mi> </msub> </mrow> </mfrac> <msub> <mi>p</mi> <mrow> <mi>s</mi> <mi>k</mi> </mrow> </msub> <msub> <mi>B</mi> <mrow> <mi>s</mi> <mi>k</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
式中:Cf[Ef(i,j)]是通信节点i和电力节点j间的通脆弱性;psk是第k类业务信息在该通道的路由选择概率;εt为通道传输实际延时和固有时延之比;n为通道中的业务量,vk为第k类业务的传输速率;Bsk为第k类通信业务的重要度指标。
电力节点脆弱性和信息层节点脆弱值和信息交互通道的脆弱值三者叠加得到考虑信息层影响下的电力节点脆弱性指标:
Cp'(Npi)=Q(i)+[1-Q(i)]·maxj{max{Cf[Ef(i,j)]},E'} (17)
在式中计算时,各指标均需要进行归一化处理,其中maxj就是与电力i节点相连的信息层通道的最大脆弱值;
考虑到信息层的融入会对电力支路的脆弱性产生间接的影响,因此考虑了信息层网络的影响下的电力支路脆弱性指标:
<mrow> <msubsup> <mi>C</mi> <mi>p</mi> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msub> <mi>E</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>=</mo> <mi>Q</mi> <mrow> <mo>(</mo> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>{</mo> <mrow> <mn>1</mn> <mo>-</mo> <mi>Q</mi> <mrow> <mo>(</mo> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mo>&amp;CenterDot;</mo> <mi>max</mi> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>w</mi> <mrow> <mi>i</mi> <mi>k</mi> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>&amp;CenterDot;</mo> <msubsup> <mi>C</mi> <mi>p</mi> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mrow> <mi>p</mi> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>w</mi> <mrow> <mi>j</mi> <mi>k</mi> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>&amp;CenterDot;</mo> <msubsup> <mi>C</mi> <mi>p</mi> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mrow> <mi>p</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
式中wik表示节点i的度数。
CN201710546839.9A 2017-07-06 2017-07-06 在信息层网络影响下的电力网络脆弱性评估方法 Pending CN107274110A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710546839.9A CN107274110A (zh) 2017-07-06 2017-07-06 在信息层网络影响下的电力网络脆弱性评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710546839.9A CN107274110A (zh) 2017-07-06 2017-07-06 在信息层网络影响下的电力网络脆弱性评估方法

Publications (1)

Publication Number Publication Date
CN107274110A true CN107274110A (zh) 2017-10-20

Family

ID=60072248

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710546839.9A Pending CN107274110A (zh) 2017-07-06 2017-07-06 在信息层网络影响下的电力网络脆弱性评估方法

Country Status (1)

Country Link
CN (1) CN107274110A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108711862A (zh) * 2018-06-07 2018-10-26 国网江苏省电力有限公司淮安供电分公司 一种间歇式能源功率波动对多层级配电网电压影响的评估方法
CN109066650A (zh) * 2018-07-16 2018-12-21 国网河北省电力有限公司经济技术研究院 电力系统脆弱性评估方法及终端设备
CN109145428A (zh) * 2018-08-14 2019-01-04 国网四川省电力公司信息通信公司 一种连锁故障模式下信息物理融合系统脆弱性评估方法
CN109936133A (zh) * 2019-02-28 2019-06-25 国网山东省电力公司青岛供电公司 考虑信息物理联合攻击的电力系统脆弱性分析方法
CN110991839A (zh) * 2019-11-21 2020-04-10 中国电建集团国际工程有限公司 一种电-气综合能源系统整体脆弱性的评估方法
CN114124728A (zh) * 2021-11-11 2022-03-01 国网四川省电力公司经济技术研究院 电网控制系统与电力一次设备间脆弱性耦合链路识别方法
CN114598612A (zh) * 2022-03-18 2022-06-07 西华大学 电力通信融合网络级联失效模型及节点脆弱性评估方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104156769A (zh) * 2013-05-31 2014-11-19 贵州电网公司电力调度控制中心 电力系统脆弱性评估方法
CN105389629A (zh) * 2015-11-10 2016-03-09 国网四川省电力公司经济技术研究院 一种结合电网结构脆弱性的电网规划方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104156769A (zh) * 2013-05-31 2014-11-19 贵州电网公司电力调度控制中心 电力系统脆弱性评估方法
CN105389629A (zh) * 2015-11-10 2016-03-09 国网四川省电力公司经济技术研究院 一种结合电网结构脆弱性的电网规划方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张海翔等: "计及保护故障因素的电网脆弱性评估模型", 《电网技术》 *
汤奕等: "考虑通信系统影响的电力系统综合脆弱性评估", 《中国电机工程学报》 *
郭静等: "基于复杂网络理论的电力通信网脆弱性分析", 《电力系统通信》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108711862A (zh) * 2018-06-07 2018-10-26 国网江苏省电力有限公司淮安供电分公司 一种间歇式能源功率波动对多层级配电网电压影响的评估方法
CN108711862B (zh) * 2018-06-07 2021-06-11 国网江苏省电力有限公司淮安供电分公司 间歇式能源功率波动对多层级配电网电压影响的评估方法
CN109066650A (zh) * 2018-07-16 2018-12-21 国网河北省电力有限公司经济技术研究院 电力系统脆弱性评估方法及终端设备
CN109145428A (zh) * 2018-08-14 2019-01-04 国网四川省电力公司信息通信公司 一种连锁故障模式下信息物理融合系统脆弱性评估方法
CN109936133A (zh) * 2019-02-28 2019-06-25 国网山东省电力公司青岛供电公司 考虑信息物理联合攻击的电力系统脆弱性分析方法
CN110991839A (zh) * 2019-11-21 2020-04-10 中国电建集团国际工程有限公司 一种电-气综合能源系统整体脆弱性的评估方法
CN110991839B (zh) * 2019-11-21 2022-11-29 中国电建集团国际工程有限公司 一种电-气综合能源系统整体脆弱性的评估方法
CN114124728A (zh) * 2021-11-11 2022-03-01 国网四川省电力公司经济技术研究院 电网控制系统与电力一次设备间脆弱性耦合链路识别方法
CN114598612A (zh) * 2022-03-18 2022-06-07 西华大学 电力通信融合网络级联失效模型及节点脆弱性评估方法

Similar Documents

Publication Publication Date Title
CN107274110A (zh) 在信息层网络影响下的电力网络脆弱性评估方法
Zhu et al. Comprehensive weight method based on game theory for identify critical transmission lines in power system
Yu et al. An UHV grid security and stability defense system: Considering the risk of power system communication
CN100438251C (zh) 一种大型电力系统脆弱线路辨识的方法
CN103279808B (zh) 恶劣天气下的电力系统线路状态转移概率预测分析方法
CN114665498A (zh) 一种计及新能源影响的主动配电网脆弱节点辨识方法
CN105912835B (zh) 一种计及风电场接入的电力线路脆弱性辨识模型建立方法
CN105893740A (zh) 一种用于挖掘电力网络中脆弱节点和线路的方法
CN104573939A (zh) 一种电网安全稳定分析方法
CN101895124B (zh) 一种判断交直流并/混联系统相互影响的耦合度评估方法
CN112072647B (zh) 考虑通信故障影响的配电网cps安全性评估方法及装置
CN110096764B (zh) 一种电-气耦合系统脆弱线路辨识和优化方法
Nasiruzzaman et al. Transient stability assessment of smart power system using complex networks framework
Guohua et al. Vulnerability assessment of bulk power grid based on complex network theory
CN108053116A (zh) 一种城市弹性配电网恢复力评估方法
CN106600458A (zh) 基于改进功率介数的电网风险综合评估方法
CN107622360A (zh) 一种综合考虑主客观因素的关键线路识别方法
CN106779269A (zh) 基于变异系数法的智能变电站二次系统可靠性评估方法
CN107272645A (zh) 中性点经电阻接地的光伏电站并网故障模型与分析方法
CN107292481A (zh) 一种基于节点重要度的电网关键节点评估方法
CN102570451A (zh) 一种输电网静态可靠性评估方法
CN104462767B (zh) 一种基于元件传递特性的配电网可靠性快速评估方法
CN103559287A (zh) 基于scd文件的智能变电站保护系统可靠性自动分析方法
CN107196307A (zh) 一种变压器故障后电网有功潮流快速估算的方法
CN105830230B (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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20171020