CN107623319B - 一种基于多评价指标的电网关键线路辨识方法 - Google Patents

一种基于多评价指标的电网关键线路辨识方法 Download PDF

Info

Publication number
CN107623319B
CN107623319B CN201710710271.XA CN201710710271A CN107623319B CN 107623319 B CN107623319 B CN 107623319B CN 201710710271 A CN201710710271 A CN 201710710271A CN 107623319 B CN107623319 B CN 107623319B
Authority
CN
China
Prior art keywords
line
index
power
power grid
load
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
CN201710710271.XA
Other languages
English (en)
Other versions
CN107623319A (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 CN201710710271.XA priority Critical patent/CN107623319B/zh
Publication of CN107623319A publication Critical patent/CN107623319A/zh
Application granted granted Critical
Publication of CN107623319B publication Critical patent/CN107623319B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种基于多评价指标的区域电网关键线路辨识方法,包括步骤:输入电力系统参数;计算线路在各项评价指标值;标准化处理后计算各项指标在主客观赋权方法下的权重;建立基于矩估计理论的最优组合权重模型来确定指标综合权重;根据雷达图计算线路的综合评价结果,识别出电网中的关键线路。本发明方法采用多指标来分析线路重要性,比单一评价指标更加全面合理。采用最优组合权重模型确定指标权重,能较好地兼顾主客观赋权方法的优点。采用雷达图法对线路进行综合评价,能形象直观地反映各方面因素对关键线路辨识结果的影响,帮助决策者制定正确的策略。本发明方法能够较好地辨识出电网的关键线路,评估结果与电网实际运行状况吻合。

Description

一种基于多评价指标的电网关键线路辨识方法
技术领域
本发明涉及电力系统领域,特别是涉及一种基于多评价指标的电网关键线路辨识方法。
背景技术
近年来国内外多次发生大停电事故。研究表明,大停电事故一般由单个元件首先发生故障,造成潮流重新分配、设备过载,继而引发连锁故障,最终导致电网崩溃。其中少数起关键作用的线路,在故障的传播过程中起着推波助澜的作用。因此,有效辨识电网中的关键线路,加强网架结构建设,对于提高电力系统可靠性具有重要的研究意义和应用价值。
现有的方法一般从网络拓扑结构和电网运行状态两个角度来辨识关键线路,例如潮流熵指标、潮流介数指标、最大流指标、计及线路电压等级和潮流熵的关键线路辨识方法、基于有功输电介数和改进潮流转移熵指标的线路重要性评估方法和基于事故链的电网连锁故障模型的关键线路辨识方法等。但目前关于电网关键线路辨识的方法大多仅考虑了线路对拓扑结构和运行状态的影响,而较少虑及线路故障对电网经济性的影响,以及线路参数及外部环境对线路故障概率的影响。此外,大部分方法采用单个评价指标来辨识关键线路,未能较为全面地考虑各方面因素对线路辨识结果的影响。
如何提取复杂电力系统的关键特征参数,选取科学有效的评价方法来对线路进行综合评价,是目前电网关键线路辨识需要深入研究的问题。可见,现有关于电网关键线路辨识的方法还有待改进。
发明内容
本发明主要解决的技术问题是采用多评价指标来定量分析电网中线路的重要性,提供一种电网关键线路的综合辨识方法。
本发明采用以下的技术方案:
综合考虑网络拓扑结构、电网实时运行参数、线路故障所导致的经济损失、线路自身特性及外部环境的影响,提出了五项评价指标来综合辨识电网关键线路:加权线路介数、改进潮流熵、局部变化量、经济损失度、线路故障率;
计算各输电线路在各项评价指标下的指标值,采用效用理论对加权线路介数和改进潮流熵的指标值进行标准化处理,根据五项指标生成电网线路重要度的评估决策矩阵;
综合考虑四种主客观赋权方法:层次分析法、改进序关系分析法、CRITIC 法和熵权法,以集成权重与主客观权重间的离差平方和最小为目标,建立多指标的主客观赋权的最优组合权重模型,用以确定指标的综合权重;
采用改进雷达图法对线路进行综合评价,根据评价结果来确定电网中的关键线路。
具体的,包括步骤:
输入电力系统拓扑结构参数、电气特性参数、线路参数、外部环境参数和经济性指标;
综合考虑网络拓扑结构、电网实时运行参数、线路故障所导致的经济损失、线路自身特性及外部环境的影响,提出了五项评价指标来综合辨识电网关键线路:加权线路介数、改进潮流熵、局部变化量、经济损失度、线路故障率;
计算各输电线路在各项评价指标下的指标值,采用效用理论对加权线路介数和改进潮流熵的指标值进行标准化处理,根据五项指标生成电网线路重要度的评估决策矩阵;
综合考虑四种主客观赋权方法:层次分析法、改进序关系分析法、CRITIC 和熵权法,以集成权重与主客观权重间的离差平方和最小为目标,建立多指标的主客观赋权的最优组合权重模型,用以确定指标的综合权重;
采用改进雷达图法对线路进行综合评价,根据评价结果来确定电网中的关键线路。
采用加权线路介数指标来辨识电网关键线路,包括:
复杂网络理论在研究电力系统的结构脆弱性和连锁故障传播机制等方面具有重要的参考价值;由于高介数线路对电网脆弱性影响较大,故采用第k条线路的加权线路介数Ak来辨识在拓扑结构中处于核心地位的线路;实际电力线路间存在着电气联系,系统有功的分布受到线路电抗的影响;故将线路电抗值Xij作为线路Li-j的权重wij,则最短电气路径为线路电抗和最小的路径;
Figure GDA0001481728530000021
式中:Zij(k)表示连接节点i和j的最短电气路径经过线路k的次数;加权线路介数属于电网的固有属性,可用于辨识电力系统网架结构的关键环节;加权线路介数越大,表示该线路发生故障对电网拓扑结构的影响越显著。
采用改进潮流熵指标来辨识电网关键线路,包括:
大部分停电事故是由于单一元件故障导致支路潮流转移,进而引起其他支路相继过载切除;由于不同支路的断开对电网功率传输的影响不同,故定义潮流熵指标来衡量线路切除后电网潮流分布的不均衡程度,进而辨识在功率传输中处于重要地位的线路;
设第k条线路正常运行时传输的有功功率为Pk0,所能承载的最大有功功率为Pkmax,其负载率为μk=|Pk0/Pkmax|;给定常数序列W=[W1,W2,…,Wr],共r-1个常数序列区间;对第t个常数序列区间(Wt,Wt+1],lt表示线路负载率μk∈(Wt,Wt+1]的线路条数,过载线路置于(1,∞)区间;则切除第k条线路后系统的潮流熵Bk′为:
Figure GDA0001481728530000031
Figure GDA0001481728530000032
式中:P(t)表示负载率μk∈(Wt,Wt+1]的线路条数占总线路数量的比例;可以看出,潮流分布越均衡,则潮流熵越小;因此,切除线路后电网的潮流熵越大,表明该线路在功率传输中地位越重要;
实际运行中,负载率越高的线路发生越界的可能性越大;传统的潮流熵指标仅以潮流转移量来衡量冲击大小,未考虑高负载率线路的脆弱性,故将负载率加入潮流熵指标中进行分析;假定负载率超过一定数值a(0<a<1)后,负载率区间的地位随着负载率的增大而上升,则第t个负载率区间(Wt,Wt+1]的重要性αt为:
Figure GDA0001481728530000033
式中:β为负载率区间的加权系数;Wt为第t个负载率区间的区间下限值。
计及负载率的影响,则第k条线路故障后电网的改进潮流熵可表示为:
Figure GDA0001481728530000041
采用局部变化量指标来辨识电网关键线路,包括:
为分析线路故障对系统局部电压稳定和无功平衡的影响,定义线路k的局部变化量指标Ck为:
Ck=max(CUk,CQk)
Figure GDA0001481728530000042
Figure GDA0001481728530000043
式中:Vd表示负荷节点集;CUk用于判断线路k断开后负荷节点是否出现电压越限;Uik表示线路k故障后节点i的电压;
Figure GDA0001481728530000044
Uic 分别为节点i的临界电压上限和下限;VG表示发电机节点集;CQk用于分析线路k故障对发电机节点无功出力的影响;Qi0和Qik分别为线路k故障前和故障后发电机节点i的无功出力;Qic为发电机节点i的无功容量;若线路k故障,导致发电机节点无功出力的增加量超过其无功容量的ε倍,说明发电机节点受到严重影响,因而也说明该线路越重要。
采用经济损失度指标来辨识电网关键线路,包括:
某些输电线路断开后系统会形成孤网,各孤网内功率可能无法平衡,此时需要采取切除负荷或调整发电机出力的措施;由于电网中不同负荷节点的单位经济损失不同,故引入ζi-L来衡量节点i处切除单位负荷的经济损失;电力系统正常运行时是一个互联的系统,当出现孤网后整个电力系统的可靠性水平将下降,从而导致系统的机组备用容量上升;假设互联系统的机组备用率为最大发电负荷的 6%,以机组的容量为比例,将系统备用容量分摊到各机组;若线路断开后形成孤网,则孤网内系统的机组备用率为最大发电负荷的10%,引入ζi-G来衡量节点i处的机组提高单位备用容量的经济损失;采用经济损失度指标Dk来分析线路k 故障所造成的经济损失:
Figure GDA0001481728530000051
式中:ψk-L和ψk-G分别为线路k退出运行后系统切负荷和机组备用容量上升导致的经济损失;ψ0为设置的损失基准,用于衡量经济损失的相对大小;ΩL和ΩG分别为需要切负荷或提高机组备用容量的节点集;Li-L和Li-G分别为节点i处的负荷损失量和机组备用容量增加量。
采用线路故障率指标来辨识电网关键线路,包括:
故障概率高的线路对电力系统的可靠运行影响较大,故辨识电网关键线路还应考虑线路自身的特性和外部环境的影响;经验表明,当气象状况的指标超过阈值时,气象条件成为影响线路故障概率最重要的因素;以台风天气为例,主要考虑线路k的断线故障率和倒塔故障率;采用指数函数来拟合断线故障率Ek1和导线最大应力σm的关系:
Figure GDA0001481728530000052
式中:σ1为导线设计应力;μ1为导线安全系数;K1和T1为与线路参数相关的常数;导线最大应力σm可根据导线风荷载、线路档距、悬点间高差角等参数求出;
考虑杆塔所受的不平衡力及杆塔风荷载超过杆塔的承受能力,得到倒塔故障率Ek2与风总荷载Qt的关系:
Figure GDA0001481728530000053
式中:Ht为杆塔设计荷载;μt为杆塔安全系数;K2和T2为与线路参数相关的常数;风总荷载Qt可根据导线风荷载、杆塔沿线路两侧的不平衡力等参数求出;
综合考虑断线故障率和倒塔故障率,得到线路k的故障率Ek为:
Ek=1-(1-Ek1)(1-Ek2)。
采用效用理论对加权线路介数和改进潮流熵的指标值进行标准化处理,包括:
从指标的定义可以看出,加权线路介数和改进潮流熵的量纲与其他三项指标不同,因此指标赋权前需先对这两项评价指标值进行标准化处理;传统的线性归一化方法会产生遮蔽现象,即各项指标值比较适中的线路的归一化值,会超过某一项指标突出而其他指标值较小的线路;这在关键线路辨识中将不利于筛选在某方面特别脆弱的线路;因此,这里采用效用理论来对指标进行归一化处理,其处理过程为:先对指标进行极差变换,再进行指数处理,这样可以避免上述的遮蔽现象,其归一化公式为:
Figure GDA0001481728530000061
Figure GDA0001481728530000062
式中:rki表示第k条线路的第i个评价指标的指标值。
建立多指标的主客观赋权的最优组合权重模型,用以确定指标的综合权重,包括:
假设m种赋权方法中有v种主观赋权方法、m-v种客观赋权方法,第i个指标在第s种赋权方法下的权重为γsi。针对不同的评价指标,主客观权重的相对重要程度不同;因此,根据矩估计理论可得第i个指标的主、客观权重的期望分别为:
Figure GDA0001481728530000063
Figure GDA0001481728530000064
进而得到第i个指标的主、客观加权系数分别为τi
Figure GDA0001481728530000074
Figure GDA0001481728530000071
Figure GDA0001481728530000072
以集成权重与主客观权重间的离差平方和最小为目标,则确定第i个指标的综合权重Mi的最优组合权重模型为:
Figure GDA0001481728530000073
采用雷达图法来对线路进行综合评价,包括:
雷达图法是一种基于形似雷达显示屏图形而构建的多变量对比分析技术;该方法通过对图像特征的数学处理,能够直观形象地反映各评价指标的独立权重、指标间的相互影响及均衡性;雷达图法通过计算所绘制的雷达图的特征参数,来给出被评估对象的综合评价结果;传统雷达图法各指标轴的夹角是等分关系,弱化了指标自身权重的影响;且扇形区域为相邻两个指标共同拥有,难以清晰地划分指标在综合评价中的重要性的不同;基于此,对传统雷达图法进行改进,利用扇形的角平分线来绘制线路的雷达图。
本发明提供的技术方案的有益效果:
本发明提供的区域电网关键线路的辨识方法,综合考虑了网络拓扑结构、实时运行状态、线路故障导致的经济损失、线路自身特性和外部环境的影响。采用多指标来分析线路重要性,比单一评价指标更加全面合理。基于最优组合权重模型来求取指标的综合权重,能较好地兼顾主客观赋权方法的优点。采用雷达图法对线路进行综合评价,不仅比简单加权法更加科学,还能形象直观地反映各方面因素对关键线路辨识结果的影响,帮助决策者制定正确的策略。所提出的方法能够较好地辨识出电网的关键线路,评估结果与电网实际运行状况吻合。
附图说明
图1为基于多评价指标的电网关键线路辨识方法流程图;
图2为多指标综合评价的雷达图;
图3为IEEE 39节点系统接线图;
图4为部分线路的雷达图。
具体实施方式
为了更好地理解本发明的目的、技术方案以及技术效果,以下结合附图对本发明进行进一步的讲解说明。
本发明提出了一种基于多评价指标的电网关键线路辨识方法,其实施流程包括如下详细步骤:
步骤1、输入电力系统拓扑结构参数、电气特性参数、线路参数、外部环境参数和经济性指标;
步骤2、计算各线路在各评价指标上的指标值,具体指标包括:
a.加权线路介数指标Ak
Figure GDA0001481728530000081
式中:Zij(k)表示连接节点i和j的最短电气路径经过线路k的次数。
b.改进潮流熵指标Bk
Figure GDA0001481728530000082
Figure GDA0001481728530000083
Figure GDA0001481728530000084
式中:W=[W1,W2,…,Wr]表示常数序列;lt表示线路负载率μk∈(Wt,Wt+1]的线路条数,过载线路置于(1,∞)区间;P(t)表示负载率μk∈(Wt,Wt+1]的线路条数占总线路数量的比例;αt为则第t个负载率区间(Wt,Wt+1]的重要性;β为负载率区间的加权系数;Wt为第t个负载率区间的区间下限值。
c.局部变化量指标Ck
Ck=max(CUk,CQk)
Figure GDA0001481728530000091
Figure GDA0001481728530000092
式中:Vd表示负荷节点集;CUk用于判断线路k断开后负荷节点是否出现电压越限;Uik表示线路k故障后节点i的电压;
Figure GDA0001481728530000093
Uic 分别为节点i的临界电压上限和下限;VG表示发电机节点集;CQk用于分析线路k故障对发电机节点无功出力的影响;Qi0和Qik分别为线路k故障前和故障后发电机节点i的无功出力;Qic为发电机节点i的无功容量。
d.经济损失度指标Dk
Figure GDA0001481728530000094
式中:ψk-L和ψk-G分别为线路k退出运行后系统切负荷和机组备用容量上升导致的经济损失;ψ0为设置的损失基准,用于衡量经济损失的相对大小;ΩL和ΩG分别为需要切负荷或提高机组备用容量的节点集;ζi-L为节点i处切除单位负荷的经济损失;ζi-G为节点i处的机组提高单位备用容量的经济损失;Li-L和 Li-G分别为节点i处的负荷损失量和机组备用容量增加量。假设互联系统的机组备用率为最大发电负荷的6%,以机组的容量为比例,将系统备用容量分摊到各机组。若线路断开后形成孤网,则孤网内系统的机组备用率为最大发电负荷的 10%。
e.线路故障率指标Ek
Ek=1-(1-Ek1)(1-Ek2)
Figure GDA0001481728530000101
Figure GDA0001481728530000102
式中:σ1为导线设计应力;μ1为导线安全系数;K1和T1为与线路参数相关的常数;导线最大应力σm可根据导线风荷载、线路档距、悬点间高差角等参数求出;Ht为杆塔设计荷载;μt为杆塔安全系数;K2和T2为与线路参数相关的常数;风总荷载Qt可根据导线风荷载、杆塔沿线路两侧的不平衡力等参数求出。
步骤3、采用效用理论对加权线路介数和改进潮流熵的指标值进行标准化处理,先对指标进行极差变换,再进行指数处理:
Figure GDA0001481728530000103
Figure GDA0001481728530000104
式中:rki表示第k条线路的第i个评价指标的指标值。
步骤4、根据各线路的五项指标值生成电网线路重要度的评估决策矩阵;
步骤5、计算各项指标在四种主客观赋权方法下的权重,四种主客观赋权方法分别是:层次分析法、改进序关系分析法、CRITIC和熵权法;
步骤6、根据矩估计理论求取各项指标的主、客观权重的期望。其中,第i个指标的主、客观权重的期望分别为U1i)和U2i):
Figure GDA0001481728530000105
Figure GDA0001481728530000106
式中:m为赋权方法的总数;v为主观赋权方法的数目;γsi为第i个指标在第s种赋权方法下的权重。
步骤7、根据各项指标的主、客观权重的期望,求取各项指标的主、客观加权系数。其中,第i个指标的主、客观加权系数分别为τi
Figure GDA0001481728530000119
Figure GDA0001481728530000111
Figure GDA0001481728530000112
步骤8、以集成权重与主客观权重间的离差平方和最小为目标,建立多指标的主客观赋权的最优组合权重模型,用以确定指标的综合权重Mi
Figure GDA0001481728530000113
步骤9、用于评价线路重要度的各项指标按综合权重从大到小进行排序,得到排序结果为uk=(uk1,uk2,…,uk5),排序后的第i项指标在雷达图中对应的扇形圆心角为θki=2πuki
步骤10、作单位圆,通过圆心Ok作射线OkPk1,与圆交于点Pk1;作射线OkPk2,与圆交于点Pk2,使∠Pk1OkPk2=θk1;同理,根据圆心角θki依次绘制射线OkPk3、OkPk4、 OkPk5;依次作∠Pk1OkPk2、∠Pk2OkPk3、…、∠Pk5OkPk1的角平分线。
步骤11、排序后的第i项指标值按比例换算为距离圆点的长度
Figure GDA0001481728530000114
所有指标值换算为长度
Figure GDA0001481728530000115
根据长度
Figure GDA0001481728530000116
在角平分线上标出相应的点Ak、Bk、…、Ek
步骤12、依次连接Ak、Pk1、Bk、Pk2、…、Ek,得到如图2所示的第k条线路综合评价的雷达图。
步骤13、根据绘制的雷达图,计算线路的综合评价结果。其中,第k条线路的综合评价结果Zk为:
Figure GDA0001481728530000117
Figure GDA0001481728530000118
Figure GDA0001481728530000121
式中:Sk为封闭多边形面积;Ck为多边形周长。
为了进一步理解本发明,以下以IEEE(Institute of Electrical andElectronics Engineers电气和电子工程师协会)39节点系统为例,来解释本发明的实际应用。
IEEE39节点系统有39个节点,其中发电机节点10个,34条线路,12条变压器支路,系统如图3所示,其中,数字表示节点编号,带圆圈的字母G表示发电机节点。参数选取a=0.8,β=1.5,K1=10-6,T1=μ1σ1/ln104,K2=10-6,T2=μtHt/ln104
首先,归一化处理线路的加权线路介数和改进潮流熵的评价指标值。通过层次分析法、G1法、CRITIC法和熵权法求出五项指标的权重后,再根据最优组合权重模型确定指标综合权重,其结果如表1所示。
表1 各项评价指标的权重
Figure GDA0001481728530000122
从表1可以看出,指标Ak和Ek的综合权重比较小。因为Ak反映的是网络的拓扑结构,忽略了电气特性;而Ek受线路参数和天气影响,不确定性高,故这两项指标权重较小是比较合理的。
接着采用改进雷达图法对线路进行综合评价,图3中红色的输电线路为辨识出的关键线路,其评价指标的归一化结果如表2所示。
表2 线路评价指标的标准化结果
Figure GDA0001481728530000131
从表2可以看出:
●评价指标间具有相似性:如L16-19、L15-16和L26-27同时具有较大的Ak、Bk和Ck指标,这说明某些关键线路的失效会同时影响电网的拓扑结构和运行特性。
●评价指标间具有互补性:如L21-22同时具有较小的Bk值和较大的Ck值,说明该线路对系统全局的有功传输分布有较大影响,而对局部电压稳定和无功平衡的影响较小。
结合表1和表2可知:指标Ck和Dk的数据间差异大,这类数据特征会导致采用客观赋权法求出的指标权重比较大。例如,熵权法中Ck和Dk的权重分别为 0.4657和0.3895,淹没了其他指标,而采用最优组合权重模型求出Ck和Dk的综合权重分别为0.3091和0.2294,避免了因数据特征而造成权重设置不合理的现象。此外,采用最优组合权重模型求出指标Ak和Ek的综合权重也比较小,与层次分析法和G1法的相应结果相类似,这说明最优组合权重模型能在一定程度上考虑专家经验的影响。
综上,评价指标间既相互联系又互相补充,故指标选取比较全面合理。采用最优组合权重模型来求取指标的综合权重,能较好地兼顾主客观赋权方法的优缺点,使权重设置更加合理。
下面采用效用风险熵指标(模型1)、加权介数指标(模型2)、电气介数指标(模型3)与本发明模型进行对比,比较结果如表3所示。同时绘制相应线路的雷达图如图4所示。
表3 关键线路辨识结果比较
Figure GDA0001481728530000141
从表3中可以看出:
●按本发明模型识别出的关键线路排序与其他模型基本相似,且有4 到5条关键线路相同。本发明模型和模型1均认为L16-19是最重要的线路;本发明模型和其他模型均认为L15-16是比较重要的线路。
●L16-19是本发明模型识别出的最重要的线路。从图4可以看出,L16-19的雷达图的多边形面积大且分布比较均衡,Ak、Bk、Ck、Dk都比较大。从图3 也可以看出,L16-19处于网络拓扑结构的核心位置,若该线路断开,会形成两个孤岛,导致机组备用容量上升而产生经济损失。从电网的运行状态看,L16-19承担了502.67MW的输送功率;此外,若该线路故障还将导致节点31的无功出力由198.89MVar上升为356.95MVar,增加了79.4%的无功出力。
●L14-15不在本发明模型识别出的前八位重要线路中,但却在其他模型辨识的前八位重要线路中。从图4可以看出,L14-15的雷达图的多边形面积小,除Ak外的其他指标值都比较小。L14-15虽在拓扑结构中地位比较重要,但该线路故障不会形成孤岛。从电网运行特性来看,L14-15承担较小的传输功率 (18.97MW)。此外,L14-15故障不会引起节点电压越限或较大的发电机无功出力变化量,而且其发生故障的概率也比较低。因此,该线路排序靠后相对合理些。
●本发明模型中排序靠前的L26-27并不在模型1和模型2识别的排序前八的线路中。由图2和图4可知,L26-27的Bk比较大、Ck为1,承担较重的功率传输任务(263.3MW),且位于发电机节点37、38向其他区域传输功率的主要通道上,这说明该线路故障将会引起全网潮流分布发生较大的改变量,增加潮流分布的不均衡性。此外,该线路故障会引起容量较小的发电机节点30 的无功出力由145.89MVar上升为190.44MVar。相比较而言,模型1和模型 2中没有考虑线路故障对全网潮流分布均衡程度和局部无功平衡的影响,因此L26-27并不在模型1和模型2识别的排序前八的线路中。
综上,本发明方法辨识电网关键线路有效,且考虑比较全面。采用雷达图法对线路进行综合评价,还能形象直观地反映线路各项指标值的特点。

Claims (9)

1.一种基于多评价指标的电网关键线路辨识方法,其特征在于,包括步骤:
输入电力系统拓扑结构参数、电气特性参数、线路参数、外部环境参数和经济性指标;
综合考虑网络拓扑结构、电网实时运行参数、线路故障所导致的经济损失、线路自身特性及外部环境的影响,提出了五项评价指标来综合辨识电网关键线路:加权线路介数、改进潮流熵、局部变化量、经济损失度、线路故障率;
计算各输电线路在各项评价指标下的指标值,采用效用理论对加权线路介数和改进潮流熵的指标值进行标准化处理,根据五项指标生成电网线路重要度的评估决策矩阵;
综合考虑四种主客观赋权方法:层次分析法、改进序关系分析法、CRITIC和熵权法,以集成权重与主客观权重间的离差平方和最小为目标,建立多指标的主客观赋权的最优组合权重模型,用以确定指标的综合权重;
采用改进雷达图法对线路进行综合评价,根据评价结果来确定电网中的关键线路。
2.根据权利要求1所述的基于多评价指标的电网关键线路辨识方法,其特征在于,采用加权线路介数指标来辨识电网关键线路,包括:
复杂网络理论在研究电力系统的结构脆弱性和连锁故障传播机制方面具有重要的参考价值;由于高介数线路对电网脆弱性影响较大,故采用第k条线路的加权线路介数Ak来辨识在拓扑结构中处于核心地位的线路;实际电力线路间存在着电气联系,系统有功的分布受到线路电抗的影响;故将线路电抗值Xij作为线路Li-j的权重wij,则最短电气路径为线路电抗和最小的路径;
Figure FDA0002616329150000011
式中:Zij(k)表示连接节点i和j的最短电气路径经过线路k的次数;加权线路介数属于电网的固有属性,可用于辨识电力系统网架结构的关键环节;加权线路介数越大,表示该线路发生故障对电网拓扑结构的影响越显著。
3.根据权利要求1所述的基于多评价指标的电网关键线路辨识方法,其特征在于,采用改进潮流熵指标来辨识电网关键线路,包括:
大部分停电事故是由于单一元件故障导致支路潮流转移,进而引起其他支路相继过载切除;由于不同支路的断开对电网功率传输的影响不同,故定义潮流熵指标来衡量线路切除后电网潮流分布的不均衡程度,进而辨识在功率传输中处于重要地位的线路;
设第k条线路正常运行时传输的有功功率为Pk0,所能承载的最大有功功率为Pkmax,其负载率为μk=|Pk0/Pkmax|;给定常数序列W=[W1,W2,…,Wr],共r-1个常数序列区间;对第t个常数序列区间(Wt,Wt+1],lt表示线路负载率μk∈(Wt,Wt+1]的线路条数,过载线路置于(1,∞)区间;则切除第k条线路后系统的潮流熵Bk′为:
Figure FDA0002616329150000021
Figure FDA0002616329150000022
式中:P(t)表示负载率μk∈(Wt,Wt+1]的线路条数占总线路数量的比例;可以看出,潮流分布越均衡,则潮流熵越小;因此,切除线路后电网的潮流熵越大,表明该线路在功率传输中地位越重要;
实际运行中,负载率越高的线路发生越界的可能性越大;传统的潮流熵指标仅以潮流转移量来衡量冲击大小,未考虑高负载率线路的脆弱性,故将负载率加入潮流熵指标中进行分析;假定负载率超过一定数值a(0<a<1)后,负载率区间的地位随着负载率的增大而上升,则第t个负载率区间(Wt,Wt+1]的重要性αt为:
Figure FDA0002616329150000023
式中:β为负载率区间的加权系数;Wt为第t个负载率区间的区间下限值;
计及负载率的影响,则第k条线路故障后电网的改进潮流熵可表示为:
Figure FDA0002616329150000024
4.根据权利要求1所述的基于多评价指标的电网关键线路辨识方法,其特征在于,采用局部变化量指标来辨识电网关键线路,包括:
为分析线路故障对系统局部电压稳定和无功平衡的影响,定义线路k的局部变化量指标Ck为:
Ck=max(CUk,CQk)
Figure FDA0002616329150000031
Figure FDA0002616329150000032
式中:Vd表示负荷节点集;CUk用于判断线路k断开后负荷节点是否出现电压越限;Uik表示线路k故障后节点i的电压;
Figure FDA0002616329150000033
Uic 分别为节点i的临界电压上限和下限;VG表示发电机节点集;CQk用于分析线路k故障对发电机节点无功出力的影响;Qi0和Qik分别为线路k故障前和故障后发电机节点i的无功出力;Qic为发电机节点i的无功容量;若线路k故障,导致发电机节点无功出力的增加量超过其无功容量的ε倍,说明发电机节点受到严重影响,因而也说明该线路越重要。
5.根据权利要求1所述的基于多评价指标的电网关键线路辨识方法,其特征在于,采用经济损失度指标来辨识电网关键线路,包括:
某些输电线路断开后系统会形成孤网,各孤网内功率可能无法平衡,此时需要采取切除负荷或调整发电机出力的措施;由于电网中不同负荷节点的单位经济损失不同,故引入ζi-L来衡量节点i处切除单位负荷的经济损失;电力系统正常运行时是一个互联的系统,当出现孤网后整个电力系统的可靠性水平将下降,从而导致系统的机组备用容量上升;假设互联系统的机组备用率为最大发电负荷的6%,以机组的容量为比例,将系统备用容量分摊到各机组;若线路断开后形成孤网,则孤网内系统的机组备用率为最大发电负荷的10%,引入ζi-G来衡量节点i处的机组提高单位备用容量的经济损失;采用经济损失度指标Dk来分析线路k故障所造成的经济损失:
Figure FDA0002616329150000034
式中:ψk-L和ψk-G分别为线路k退出运行后系统切负荷和机组备用容量上升导致的经济损失;ψ0为设置的损失基准,用于衡量经济损失的相对大小;ΩL和ΩG分别为需要切负荷或提高机组备用容量的节点集;Li-L和Li-G分别为节点i处的负荷损失量和机组备用容量增加量。
6.根据权利要求1所述的基于多评价指标的电网关键线路辨识方法,其特征在于,采用线路故障率指标来辨识电网关键线路,包括:
故障概率高的线路对电力系统的可靠运行影响较大,故辨识电网关键线路还应考虑线路自身的特性和外部环境的影响;经验表明,当气象状况的指标超过阈值时,气象条件成为影响线路故障概率最重要的因素;以台风天气为例,主要考虑线路k的断线故障率和倒塔故障率;采用指数函数来拟合断线故障率Ek1和导线最大应力σm的关系:
Figure FDA0002616329150000041
式中:σ1为导线设计应力;μ1为导线安全系数;K1和T1为与线路参数相关的常数;导线最大应力σm可根据导线风荷载、线路档距、悬点间高差角参数求出;
考虑杆塔所受的不平衡力及杆塔风荷载超过杆塔的承受能力,得到倒塔故障率Ek2与风总荷载Qt的关系:
Figure FDA0002616329150000042
式中:Ht为杆塔设计荷载;μt为杆塔安全系数;K2和T2为与线路参数相关的常数;风总荷载Qt可根据导线风荷载、杆塔沿线路两侧的不平衡力参数求出;
综合考虑断线故障率和倒塔故障率,得到线路k的故障率Ek为:
Ek=1-(1-Ek1)(1-Ek2)。
7.根据权利要求1所述的基于多评价指标的电网关键线路辨识方法,其特征在于,采用效用理论对加权线路介数和改进潮流熵的指标值进行标准化处理,包括:
从指标的定义可以看出,加权线路介数和改进潮流熵的量纲与其他三项指标不同,因此指标赋权前需先对这两项评价指标值进行标准化处理;传统的线性归一化方法会产生遮蔽现象,即各项指标值比较适中的线路的归一化值,会超过某一项指标突出而其他指标值较小的线路;这在关键线路辨识中将不利于筛选在某方面特别脆弱的线路;因此,这里采用效用理论来对指标进行归一化处理,其处理过程为:先对指标进行极差变换,再进行指数处理,这样可以避免上述的遮蔽现象,其归一化公式为:
Figure FDA0002616329150000051
Figure FDA0002616329150000052
式中:rki表示第k条线路的第i个评价指标的指标值。
8.根据权利要求1所述的基于多评价指标的电网关键线路辨识方法,其特征在于,建立多指标的主客观赋权的最优组合权重模型,用以确定指标的综合权重,包括:
假设m种赋权方法中有v种主观赋权方法、m-v种客观赋权方法,第i个指标在第s种赋权方法下的权重为γsi;针对不同的评价指标,主客观权重的相对重要程度不同;因此,根据矩估计理论可得第i个指标的主、客观权重的期望分别为
Figure FDA0002616329150000053
Figure FDA0002616329150000054
进而得到第i个指标的主、客观加权系数分别为τi
Figure FDA0002616329150000057
Figure FDA0002616329150000055
Figure FDA0002616329150000056
以集成权重与主客观权重间的离差平方和最小为目标,则确定第i个指标的综合权重Mi的最优组合权重模型为:
Figure FDA0002616329150000061
9.根据权利要求1所述的基于多评价指标的电网关键线路辨识方法,其特征在于,采用雷达图法来对线路进行综合评价,包括:
雷达图法是一种基于形似雷达显示屏图形而构建的多变量对比分析技术;该方法通过对图像特征的数学处理,能够直观形象地反映各评价指标的独立权重、指标间的相互影响及均衡性;雷达图法通过计算所绘制的雷达图的特征参数,来给出被评估对象的综合评价结果;传统雷达图法各指标轴的夹角是等分关系,弱化了指标自身权重的影响;且扇形区域为相邻两个指标共同拥有,难以清晰地划分指标在综合评价中的重要性的不同;基于此,对传统雷达图法进行改进,利用扇形的角平分线来绘制线路的雷达图。
CN201710710271.XA 2017-08-17 2017-08-17 一种基于多评价指标的电网关键线路辨识方法 Active CN107623319B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710710271.XA CN107623319B (zh) 2017-08-17 2017-08-17 一种基于多评价指标的电网关键线路辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710710271.XA CN107623319B (zh) 2017-08-17 2017-08-17 一种基于多评价指标的电网关键线路辨识方法

Publications (2)

Publication Number Publication Date
CN107623319A CN107623319A (zh) 2018-01-23
CN107623319B true CN107623319B (zh) 2021-03-02

Family

ID=61088892

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710710271.XA Active CN107623319B (zh) 2017-08-17 2017-08-17 一种基于多评价指标的电网关键线路辨识方法

Country Status (1)

Country Link
CN (1) CN107623319B (zh)

Families Citing this family (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108493979B (zh) * 2018-03-28 2021-01-26 国电南瑞科技股份有限公司 一种特高压直流接入省级电网分区优化运行控制方法
CN108763161B (zh) * 2018-05-31 2019-07-09 哈尔滨工业大学 一种基于多层次指标体系的电梯安全等级评价方法
CN109375050B (zh) * 2018-08-22 2021-01-01 山东科技大学 一种基于初始故障线路割集的电力系统关键断面辨识方法
CN108879708B (zh) * 2018-08-28 2021-06-22 东北大学 一种主动配电网的无功电压分区方法及系统
CN108899905B (zh) * 2018-09-03 2021-08-03 云南电网有限责任公司电力科学研究院 一种复杂电网中关键节点的辨识方法及装置
CN109685374A (zh) * 2018-12-26 2019-04-26 国网江苏省电力有限公司 一种基于复杂网络的电网关键单元辨识评估方法
CN109687410B (zh) * 2018-12-27 2020-11-13 国网福建省电力有限公司 一种基于最大熵增的连锁故障集构建方法
CN109830963B (zh) * 2019-01-08 2020-07-07 国网山东省电力公司滨州供电公司 一种计及电压等级影响的电网关键节点识别方法
CN110288489B (zh) * 2019-07-01 2022-04-12 华北电力大学(保定) 一种电网脆弱线路分类方法及装置
CN110277835B (zh) * 2019-07-03 2020-12-15 浙江大学 一种基于用电信息采集系统的低压用户过电压风险监测方法
CN110390472B (zh) * 2019-07-03 2022-03-25 浙江大学 一种基于数据驱动方法的低压用户用电异常评估方法
CN110768260B (zh) * 2019-09-12 2022-09-06 南京邮电大学 一种基于电气介数的电网级联失效模型建立方法
CN111092453B (zh) * 2019-12-30 2023-04-07 南京工程学院 一种多点接入光伏电源的电网关键线路辨识方法
CN111191946A (zh) * 2020-01-02 2020-05-22 国网北京市电力公司 故障线路识别方法及装置、电子设备
CN111275312B (zh) * 2020-01-17 2024-07-12 中国电力科学研究院有限公司 一种用于对电网多重故障严重性进行综合评价排序的方法及系统
CN111444593B (zh) * 2020-03-02 2022-05-03 浙江大学 一种电-气综合能源系统元件脆弱性改善方法
CN111313408B (zh) * 2020-03-05 2022-04-29 西南交通大学 一种考虑暂态能量关联的电网脆弱线路辨识方法
CN111553565B (zh) * 2020-04-10 2021-11-16 中国南方电网有限责任公司超高压输电公司检修试验中心 一种变电站在线监测装置性能评估方法及系统
CN112003734B (zh) * 2020-07-28 2022-06-07 长沙理工大学 一种基于改进结构熵的信息物理系统关键节点辨识方法
CN112421638B (zh) * 2020-11-19 2022-07-22 国家电网公司西南分部 一种基于潮流贡献度的直流近区资产划分方法
CN112701711B (zh) * 2020-11-30 2021-09-28 中国南方电网有限责任公司超高压输电公司广州局 一种交直流合建换流站重要线路运维安全性检测评价方法
CN112488336A (zh) * 2020-12-09 2021-03-12 深圳供电局有限公司 一种用于极端天气下的电网故障处理方法及系统
CN112783939A (zh) * 2020-12-31 2021-05-11 广东电网有限责任公司 一种基于数据挖掘的低压配电网运行状态评估方法
CN112836950A (zh) * 2021-01-26 2021-05-25 国网浙江省电力有限公司经济技术研究院 一种台风灾害条件下电力线路的故障率评估方法和装置
CN113705986B (zh) * 2021-08-12 2023-11-07 武汉大学 一种基于风险理论的核心骨干网架构建方法
CN113762751A (zh) * 2021-08-30 2021-12-07 国网冀北电力有限公司电力科学研究院 机组功率调节参数权重确定方法及装置
CN114066015A (zh) * 2021-10-22 2022-02-18 山东旗帜信息有限公司 一种部门人员组成优化方法、设备及介质
CN114998051B (zh) * 2022-05-30 2024-04-16 西南石油大学 一种基于最大流理论的电-气耦合系统关键线路辨识方法
CN115545358A (zh) * 2022-12-01 2022-12-30 中国民航工程咨询有限公司 一种冰雪状态下机场跑道适航性评估方法、系统及应用
CN115865740B (zh) * 2023-03-02 2023-05-12 中国人民解放军国防科技大学 一种基于网络结构的关键链路识别方法及装置
CN117913827B (zh) * 2024-03-18 2024-05-28 广东电网有限责任公司广州供电局 一种计及触发功能的复杂配电网的优化方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105574614A (zh) * 2015-12-15 2016-05-11 东北大学 基于多属性分析的电力系统线路脆弱度评估系统及方法
CN106529791A (zh) * 2016-10-27 2017-03-22 云南电网有限责任公司 一种电力系统中支路重要度的评估方法
CN106779309A (zh) * 2016-11-24 2017-05-31 华北电力大学(保定) 关键线路多角度多层次的辨识方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105574614A (zh) * 2015-12-15 2016-05-11 东北大学 基于多属性分析的电力系统线路脆弱度评估系统及方法
CN106529791A (zh) * 2016-10-27 2017-03-22 云南电网有限责任公司 一种电力系统中支路重要度的评估方法
CN106779309A (zh) * 2016-11-24 2017-05-31 华北电力大学(保定) 关键线路多角度多层次的辨识方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于信息熵权和层次分析法的电网关键节点识别;林鸿基等;《广东电力》;20161231;第29卷(第12期);50-56 *

Also Published As

Publication number Publication date
CN107623319A (zh) 2018-01-23

Similar Documents

Publication Publication Date Title
CN107623319B (zh) 一种基于多评价指标的电网关键线路辨识方法
CN107563680B (zh) 一种基于ahp和熵权法的配电网可靠性评估方法
CN106850254B (zh) 一种电力通信网中关键节点识别方法
CN107992962B (zh) 一种基于熵权法的输电线路防雷措施优化选择方法
CN110034581B (zh) 风电并网下电力系统的区间电气介数脆弱性评估方法
CN107895947A (zh) 一种配电网脆弱线路辨识方法
CN104063750B (zh) 基于改进ahp‑反熵权的电力系统受灾影响的预测方法
CN107633320B (zh) 一种基于气象预测和风险评估的电网线路重要度评估方法
CN104392391A (zh) 一种电网运行安全风险量化方法
CN111210102B (zh) 一种极端灾害条件下电网故障建模方法
CN103489133B (zh) 一种自动推荐输电线路防雷措施的方法
CN102801732A (zh) 一种基于网络拓扑结构的电力通信骨干网安全风险评估方法
CN110705879A (zh) 一种高比例可再生能源接入下的电网脆弱性评估方法
CN108681786B (zh) 基于供电社团结构的分布式发电选址规划方法
CN110535144A (zh) 风沙天气下含多类型负荷的智能配电网韧性定量分析方法
CN110378549B (zh) 一种基于fahp-熵权法的输电杆塔鸟害等级评估方法
CN110428191B (zh) 配电网脆弱节点辨识的方法
CN104156890A (zh) 一种风电并网方案决策方法
CN103093097A (zh) 基于规范切的电力系统脆弱断面识别方法
CN110826842A (zh) 一种台风场景下的配电网风险控制方法
CN112668821A (zh) 一种基于风沙区域绝缘子故障概率的配电线路风险分析方法
CN107239896A (zh) 一种继电保护安全性评估方法
CN112418461A (zh) 一种智慧输电的巡检系统及方法
CN115545479A (zh) 一种配电网重要节点或重要线路的确定方法及装置
WO2021077596A1 (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