CN109829657B - 一种考虑关键节点的电网临近连锁跳闸的危险度计算方法 - Google Patents
一种考虑关键节点的电网临近连锁跳闸的危险度计算方法 Download PDFInfo
- Publication number
- CN109829657B CN109829657B CN201910150822.0A CN201910150822A CN109829657B CN 109829657 B CN109829657 B CN 109829657B CN 201910150822 A CN201910150822 A CN 201910150822A CN 109829657 B CN109829657 B CN 109829657B
- Authority
- CN
- China
- Prior art keywords
- power grid
- power
- branch
- grid
- formula
- 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
Images
Classifications
-
- 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
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Abstract
本发明公开了一种考虑关键节点的电网临近连锁跳闸的危险度计算方法,步骤一,针对电网的连锁跳闸场景,利用直流潮流法,给出电网在初始故障发生后的支路电流与节点电气控制量与扰动量之间的映射关系;步骤二,形成一个新的向量;步骤三,根据式(5)给出的数学模型,定义惩罚函数,并选择粒子群优化算法进行求解。本发明的主要优点是在严格分析的基础上,只考虑电网中影响连锁跳闸的关键节点,以此计算电网当前状态距离连锁跳闸事件发生的危险程度,可大大降低计算的时间,提高分析效率。
Description
技术领域
本发明涉及电气领域,尤其涉及一种考虑关键节点的电网临近连锁跳闸的危险度计算方法。
背景技术
现如今的电力系统,除了不断地采用新技术增加电力的可靠性外,通过互联将其规模扩大是一个必然的趋势,可以预计,在未来,电力系统规模还将扩大,而其运行则变得越来越复杂,在这样的背景下,一方面电网可以使能源得到更加高效合理的利用,但另一方面,电力系统的安全运行问题变得日益复杂,因为在大电网中发生的故障往往相互具有连锁效应,元件之间相互波及,一旦连锁效应放大,导致大范围的停电,其损失往往是很严重的。近年来世界范围内的一些大停电事故已经很深刻地说明了这一点,所以,近年来,电网的连锁跳闸问题或与之有关的连锁故障问题受到了许多电力工作者的关注。
目前,对电网连锁跳闸的分析和处理技术主要可以分为两大类,其中第一大类是将连锁跳闸置于长过程的连锁故障过程中,分析和处理整个连锁故障过程的相关问题;另一类是针对初始故障发生后的局部时段分析和研究连锁跳闸以及相关的预防措施。
(1)针对连锁故障长过程的分析方法和技术
长过程的连锁故障模型,包括连锁故障的发生机理,电网结构对连锁故障的影响以及长过程连锁故障的模拟。其中研究和分析研究故障的发生机理,主要是运用自组织临界理论来解释[1-6],目前,一些研究工作者正在开展将此理论用于指导实际的连锁故障分析的研究中,如文献[7]给出了通过连锁故障的搜索判别系统进入自组织临界态的算法,文献[8]]给出了运用自组织临界理论判别系统可能的停电规模的方法。这些研究和技术对于连锁故障场景中的任何一级故障一般都按照连锁跳闸的模式进行分析,一般不涉及连锁故障中的复杂动态现象。
研究和分析电网结构对连锁故障的影响,目前主要是基于复杂网络理论分析电网的拓扑结构故障传播的影响[9-12],如文献[13]给出了小世界电网的故障传播规律,文献[14-15]给出辨识影响电网故障传播的关键结构参数的方法。这些研究和技术对于连锁故障场景中的任何一级故障一般也都按照连锁跳闸的模式进行分析,一般也不涉及连锁故障中的复杂动态现象。
对于现有长过程连锁故障的模拟技术,目前,既有基于多级连锁跳闸模拟的分析技术[16],也有在涉及复杂动态的情况进行模拟的技术[17-19],不过无论是哪种技术,目前还很难完全真实地再现完整的连锁故障动态过程,主要原因是这些动态过程过于复杂,涉及到的自动装置也没有完全统一的标准。
(2)针对连锁故障某一局部时段的分析方法和技术
针对连锁故障的某一局部时段的分析方法和技术,一般主要针对连锁故障的早期阶段的连锁跳闸事件进行分析,之所以这样做,是因为连锁故障早期的发展动因是故障线路停运造成的电网潮流的重新调整,其发展速度较为缓慢,有较为充裕的时间采取措施,而且对连锁故障的预防较为有利。近年来的主要方法和技术有:1)基于广域信息的潮流转移识别技术[20-21],以及基于潮流转移识别的广域后备保护技术[22];2)基于本地信息的潮流转移识别和系统保护技术[22,23];3)对初始故障产生的潮流转移所关联的输电断面进行识别的技术[24]以及输电断面保护技术[25];4)根据初始故障发生后电网是否因潮流转移诱发下一轮连锁过载跳闸,辨识电网的脆弱支路的技术[26-28];5)根据初始故障,根据电网的节点注入功率计算电网针对连锁跳闸的安全裕度[29]。而传统分析和应对连锁跳闸的技术主要是从九十年代以后开始在电力系统广泛使用的电力系统静态安全分析和控制技术[20-31]。
与本发明最相近的技术方法主要是文献[29]提出的计算电网对于连锁跳闸事件的安全裕度计算方法,该文献中所考虑的安全裕度计算是针对电网中的所有节点而进行的,计算效率受到一定的限制,与本发明中提到的方法是不同的。
术语解释:
连锁跳闸:当电网中某条线路因故障被跳闸切除后,引起电力系统其他线路的电气量发生变化,当这些线路中的某一些因电气量进入继电保护的动作区而被保护切除,则这样的现象称为连锁跳闸。
初始故障:电网连锁跳闸现象中出现的第一个故障,或在分析电网连锁跳闸现象时所考虑的第一个故障。
电网连锁跳闸临界运行状态:当电网中发生初始故障后,该初始故障是否可以引发连锁故障,这与电网的运行状态是密切相关的,同样的一个初始故障,在有的运行状态下可以引发连锁故障,而在有的运行状态下则不能引发连锁故障,刚好使得该初始故障能够触发连锁故障的运行状态即是电网临界运行状态。
发明内容
本发明的目的在于提供一种计算时间短、分析高效的考虑关键节点的电网临近连锁跳闸的危险度计算方法。
为实现上述目的,本发明采用以下技术方案:
一种考虑关键节点的电网临近连锁跳闸的危险度计算方法,包括以下步骤:
步骤一,针对电网的连锁跳闸场景,利用直流潮流法,给出电网在初始故障发生后的支路电流与节点电气量如下式所示的数学关系:
I=f(u,p)(1)
式(1)中,I为初始故障发生后电网中剩余支路的电流向量,u为电网节点上的控制量,p为电网节点上的扰动量;
进一步,利用式(1)及各剩余支路的受扰情况,筛选出对严重受扰支路贡献较大的节点及其相应的电气量u1和p1,并将对严重受扰支路贡献不大的节点电气量记为u2和p2;
步骤二,以u1和p1形成一个新的向量,表示为:
O=[u1,p1](2)
并分别用OS和OC作为O在电网当前运行状态及连锁跳闸临界运行状态时的表示形式,由此,进一步给出电网当前运行状态与连锁跳闸临界运行状态之间的距离表示形式:
D=||OC-OS||(3)
式(3)中,D为OC与OS之差的范数;
再进一步,将如下的约束条件考虑进来:
式(4)中,x为电网的状态变量;h1代表初始故障发生前电网的潮流约束;h2代表初始故障切除后电网的潮流约束;g1为初始故障发生前电网的不等约束关系,包括由电网的节点电压和发电机所发的有功、无功功率以及支路功率参量的上下限所确定的不等式约束关系;h3和g2代表电网处于连锁跳闸临界运行状态时其支路电流的等式约束和不等式约束关系;
综合式(3)~(4),本发明给出如下用以计算电网当前运行状态与连锁跳闸临界运行状态之间的最短距离的数学模型:
此数学模型也即计算电网当前状态距离连锁跳闸事件发生的危险程度的模型;
步骤三,根据式(5)给出的数学模型,定义惩罚函数,并选择粒子群优化算法进行求解。
所述步骤一中,利用式(1)中,I,u,p以及I与u,p之间的映射关系中的元素,并利用K均值聚类算法筛选出对严重受扰支路贡献较大的节点。
本发明的主要优点是在严格分析的基础上,只考虑电网中影响连锁跳闸的关键节点,以此计算电网当前状态距离连锁跳闸事件发生的危险程度,可大大降低计算的时间,提高分析效率。
附图说明
以下结合附图和具体实施方式对本发明做进一步详细说明:
图1为本发明的流程图。
具体实施方式
如图1所示,本发明一种考虑关键节点的电网临近连锁跳闸的危险度计算方法,包括以下步骤:
步骤一,针对电网的连锁跳闸场景,利用直流潮流法,给出电网在初始故障发生后的支路电流与节点电气量如下式所示的数学关系:
I=f(u,p)(1)
式(1)中,I为初始故障发生后电网中剩余支路的电流向量,u为电网节点上的控制量,p为电网节点上的扰动量;
进一步,根据各剩余支路的受扰情况,利用式(1)中的I,U,P以及I与U,P之间的映射关系中的元素,利用K均值聚类算法筛选出对严重受扰支路贡献较大的节点及其相应的电气量u1和p1,并将对严重受扰支路贡献不大的节点电气量记为u2和p2;
步骤二,以u1和p1形成一个新的向量,表示为:
O=[u1,p1](2)
并分别用OS和OC作为O在电网当前运行状态及连锁跳闸临界运行状态时的表示形式,由此,进一步给出电网当前运行状态与连锁跳闸临界运行状态之间的距离表示形式:
D=||OC-OS||(3)
式(3)中,D为OC与OS之差的范数;
再进一步,将如下的约束条件考虑进来:
式(4)中,x为电网的状态变量;h1代表初始故障发生前电网的潮流约束;h2代表初始故障切除后电网的潮流约束;g1为初始故障发生前电网的不等约束关系,包括由电网的节点电压和发电机所发的有功、无功功率以及支路功率参量的上下限所确定的不等式约束关系;h3和g2代表电网处于连锁跳闸临界运行状态时其支路电流的等式约束和不等式约束关系;
综合式(3)~(4),本发明给出如下用以计算电网当前运行状态与连锁跳闸临界运行状态之间的最短距离的数学模型:
此数学模型也即计算电网当前状态距离连锁跳闸事件发生的危险程度的模型;
步骤三,根据式(5)给出的数学模型,定义惩罚函数,并选择粒子群优化算法进行求解。
本发明的主要优点是在严格分析的基础上,只考虑电网中影响连锁跳闸的关键节点,以此计算电网当前状态距离连锁跳闸事件发生的危险程度,可大大降低计算的时间,提高分析效率。
参考文献:
[1]Dobson I,Carreras B A,Lynch V E et al.An initial model for complexdynamics in electric power system blackouts[C].Hawaii InternationalConference on System Science,Hawaii,2001.
[2]Dobson I,Carreras B A,Newman D E.A probabilistic loading dependentmodel of cascading failure and possible implications for blackouts[C].HawaiiInternational Conference on System Science,Hawaii,2003.
[3]Dobson I,Carreras B A,Newman D E.A branching process approximationto cascading load-dependent system failure[C].35th Hawaii InternationalConference on System Science,Hawaii,2004.
[4]曹一家,江全元,丁理杰.电力系统大停电的自组织临界现象[J].电网技术,2005,29(15):1-5.
[5]梅生伟,翁晓峰,薛安成,等.基于最优潮流的停电模型及自组织临界性分析[J].电力系统自动化,2006,30(13):1-5,32.
[6]Carreras B A,Newman D E,Dobson I,et al.Initial evidence for self-organized criticality in electric power system blackouts[C].HawaiiInternational Conference on System Science,Hawaii,2000.
[9]石立宝,史中英,姚良忠,等.现代电力系统连锁性大停电事故机理研究综述[J].电网技术,2010,34(3):48-54.
[10]薛禹胜,谢云云,文福拴,等.关于电力系统相继故障研究的评述[J].电力系统自动化,2013,37(19):1-9,40.
[11]刘友波,胡斌,刘俊勇,等.电力系统连锁故障分析理论与应用(一)——相关理论方法与应用[J].电力系统保护与控制,2013,41(9):148-155.
[12]曹一家,王光增,包哲静,等.一种复杂电力网络的时空演化模型[J].电力自动化设备,2009,29(1):1-5.
[13]Ding Jian,Bai Xiaomin,Zhao Wei,et al.The improvement of thesmall-world network model and its application research in bulk power system[C].Proceedings of 2006International Conference on Power System Technology,2006,Chongqing,China.
[14]张国华,张建华,杨京燕,等.基于有向权重图和复杂网络理论的大型电力系统脆弱性评估[J].电力自动化设备,2009,29(4):21-26.
[15]徐林,王秀丽,王锡凡.电气介数及其在电力系统关键线路识别中的应用[J].中国电机工程学报,2010,30(1):33-39.
[16]韩学军,石磊,朱岩,等.考虑多重故障的连锁过载分析[J].电网技术,2008,32(16):86-89.
[17]王英英,罗毅,涂光瑜,等.采用图论的电网连锁故障模式搜索方法[J].高电压技术,2010,36(2):401-405.
[18]王安斯,罗毅,涂光瑜,等.用于预防控制的电力系统连锁故障事故链在线生成方法[J].高电压技术,2009,35(10):2446-2551.
[19]刁塑,刘友波,刘俊勇,许立雄,Gareth Taylor,Maozhen Li.电力系统连锁故障的多层时序运行演化模型与应用[J].中国电机工程学报,2015,(S1):82-92.
[20]徐慧明,毕天姝,黄少锋,等.基于WAMS的潮流转移识别算法[J].电力系统自动化,2006,30(14):14-19.
[21]徐慧明,毕天姝,黄少锋,等.基于广域同步测量系统的预防连锁跳闸控制策略[J].中国电机工程学报,2007,27(19):32-38.
[22]董新洲,丁磊,刘琨,等.基于本地信息的系统保护[J].中国电机工程学报,2010,30(22):7-13.
[23]董新洲,丁磊.基于变电站本地信息的电网系统保护方法:中国,200910089563.1[P].2009-07.
[24]姚峰,张保会,周徳才,等.输电断面有功安全性保护及其快速算法[J].中国电机工程学报,2006,26(13):31-36.
[25]张保会,姚峰,周德才,等.输电断面安全性保护及其关键技术研究[J].中国电机工程学报,2006,26(21):1-7.
[26]邓慧琼,王帅,孙克军,等.一种考虑二级连锁过载的电网故障关联区域划分方法[J].河北科技大学学报,2012,33(3):237-243.
[27]李勇,刘俊勇,刘晓宇,等.基于潮流熵测度的连锁故障脆弱线路评估及其在四川主干电网中的应用[J].电力自动化设备,2013,33(10):40-46.
[28]李杨,苏慧玲.N-K故障下影响电力系统脆弱性的关键线路研究[J].2015,35(3):60-67.
[29]邓慧琼.一种考虑连锁故障的电网运行安全裕度计算方法[P].2015,ZL201510295192.8
[30]诸俊伟.电力系统分析[M].北京:中国电力出版社,1995.
[31]吴际舜.电力系统静态安全分析[M].上海:上海交通大学出版社,1985.
Claims (2)
1.一种考虑关键节点的电网临近连锁跳闸的危险度计算方法,其特征在于:包括以下步骤:
步骤一,针对电网的连锁跳闸场景,利用直流潮流法,给出电网在初始故障发生后的支路电流与节点电气量如下式所示的数学关系:
I=f(u,p) (1)
式(1)中,I为初始故障发生后电网中剩余支路的电流向量,u为电网节点上的控制量,p为电网节点上的扰动量;
进一步,利用式(1)及各剩余支路的受扰情况,筛选出对严重受扰支路贡献大的节点及其相应的电气量u1和p1,并将对严重受扰支路贡献不大的节点电气量记为u2和p2;
步骤二,以u1和p1形成一个新的向量,表示为:
O=[u1,p1] (2)
并分别用OS和OC作为O在电网当前运行状态及连锁跳闸临界运行状态时的表示形式,由此,进一步给出电网当前运行状态与连锁跳闸临界运行状态之间的距离表示形式:
D=||OC-OS|| (3)
式(3)中,D为OC与OS之差的范数;
再进一步,将如下的约束条件考虑进来:
式(4)中,x为电网的状态变量;h1代表初始故障发生前电网的潮流约束;h2代表初始故障切除后电网的潮流约束;g1为初始故障发生前电网的不等约束关系,包括由电网的节点电压和发电机所发的有功、无功功率以及支路功率参量的上下限所确定的不等式约束关系;h3和g2代表电网处于连锁跳闸临界运行状态时其支路电流的等式约束和不等式约束关系;
综合式(3)~(4),给出如下用以计算电网当前运行状态与连锁跳闸临界运行状态之间的最短距离的数学模型:
此数学模型也即计算电网当前状态距离连锁跳闸事件发生的危险程度的模型;
步骤三,根据式(5)给出的数学模型,定义惩罚函数,并选择粒子群优化算法进行求解。
2.根据权利要求1所述的一种考虑关键节点的电网临近连锁跳闸的危险度计算方法,其特征在于:所述步骤一中,利用式(1)中,I,u,p以及I与u,p之间的映射关系中的元素,并利用K均值聚类算法筛选出对严重受扰支路贡献大的节点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910150822.0A CN109829657B (zh) | 2019-02-28 | 2019-02-28 | 一种考虑关键节点的电网临近连锁跳闸的危险度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910150822.0A CN109829657B (zh) | 2019-02-28 | 2019-02-28 | 一种考虑关键节点的电网临近连锁跳闸的危险度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109829657A CN109829657A (zh) | 2019-05-31 |
CN109829657B true CN109829657B (zh) | 2022-05-31 |
Family
ID=66864854
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910150822.0A Active CN109829657B (zh) | 2019-02-28 | 2019-02-28 | 一种考虑关键节点的电网临近连锁跳闸的危险度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109829657B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110601189B (zh) * | 2019-09-20 | 2022-09-16 | 福建工程学院 | 一种考虑连锁跳闸的电网安全性评估方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104242314B (zh) * | 2014-08-26 | 2017-01-11 | 清华大学 | 一种风电场汇集区域电网电压安全域值的确定方法 |
CN104901306B (zh) * | 2015-06-02 | 2017-03-22 | 福建工程学院 | 一种考虑连锁故障的电网运行安全裕度计算方法 |
CN108539719B (zh) * | 2018-03-15 | 2019-07-09 | 福建工程学院 | 一种针对配置电流型后备保护的电网的连锁跳闸预防方法 |
CN109116184B (zh) * | 2018-07-05 | 2020-11-03 | 福建工程学院 | 一种计算电网对于连锁跳闸事件危险度的方法 |
-
2019
- 2019-02-28 CN CN201910150822.0A patent/CN109829657B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109829657A (zh) | 2019-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rios et al. | Value of security: modeling time-dependent phenomena and weather conditions | |
CN103279808B (zh) | 恶劣天气下的电力系统线路状态转移概率预测分析方法 | |
CN104268410B (zh) | 一种连锁故障多阶段动态博弈防御方法 | |
WO2011101189A2 (en) | Energy generating system and control thereof | |
CN103440400A (zh) | 考虑灾害因素的电力系统短期风险判定方法 | |
CN108389002A (zh) | 一种基于n-1严重故障集的多重故障生成方法 | |
CN101923685A (zh) | 一种基于断线故障概率预测的电力切负荷决策系统及方法 | |
Wehenkel | Emergency control and its strategies | |
Panwar et al. | Integration of flow battery for resilience enhancement of advanced distribution grids | |
CN106548265A (zh) | 一种基于连锁故障事故链搜索的输电网可靠性评估方法 | |
Raj et al. | Controlled islanding of transmission system using synchrophasor measurements | |
CN103344881A (zh) | 基于组合型交叉熵算法的电网故障诊断方法 | |
CN109829657B (zh) | 一种考虑关键节点的电网临近连锁跳闸的危险度计算方法 | |
Ma et al. | A measurement-simulation hybrid method for transient stability assessment and control based on the deviation energy | |
CN111209535A (zh) | 一种电力设备相继故障风险识别方法及系统 | |
Li et al. | A comprehensive method for fault location of active distribution network based on improved matrix algorithm and optimization algorithm | |
Biswal et al. | Integrated wide-area backup protection algorithm during stressed power system condition in presence of wind farm | |
CN112613676A (zh) | 一种高分辨度电网弹性评估方法 | |
CN110190596B (zh) | 一种考虑连锁跳闸的电网临界状态搜索方法 | |
He et al. | Power system frequency situation prediction method based on transfer learning | |
CN111404163A (zh) | 一种电磁环网开环方法 | |
Yunus et al. | Reliability and availability study on substation automation system based on IEC 61850 | |
CN106159940A (zh) | 基于电网负荷特性分析的pmu最优布点方法 | |
Zima et al. | Wide area monitoring and control as a tool for mitigation of cascading failures | |
Ashok et al. | A comprehensive review on wide-area protection, control and monitoring systems |
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 |