CN109687410B - 一种基于最大熵增的连锁故障集构建方法 - Google Patents

一种基于最大熵增的连锁故障集构建方法 Download PDF

Info

Publication number
CN109687410B
CN109687410B CN201811608477.2A CN201811608477A CN109687410B CN 109687410 B CN109687410 B CN 109687410B CN 201811608477 A CN201811608477 A CN 201811608477A CN 109687410 B CN109687410 B CN 109687410B
Authority
CN
China
Prior art keywords
power
power grid
entropy
line
weighted
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
CN201811608477.2A
Other languages
English (en)
Other versions
CN109687410A (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.)
China Electric Power Research Institute Co Ltd CEPRI
State Grid Fujian Electric Power Co Ltd
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
State Grid Fujian Electric Power 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 China Electric Power Research Institute Co Ltd CEPRI, State Grid Fujian Electric Power Co Ltd filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN201811608477.2A priority Critical patent/CN109687410B/zh
Publication of CN109687410A publication Critical patent/CN109687410A/zh
Application granted granted Critical
Publication of CN109687410B publication Critical patent/CN109687410B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H7/00Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
    • H02H7/26Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured
    • H02H7/261Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured involving signal transmission between at least two stations
    • H02H7/262Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured involving signal transmission between at least two stations involving transmissions of switching or blocking orders
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H3/00Emergency protective circuit arrangements for automatic disconnection directly responsive to an undesired change from normal electric working condition with or without subsequent reconnection ; integrated protection
    • H02H3/02Details
    • H02H3/06Details with automatic reconnection
    • 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

  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种基于最大熵增的连锁故障集构建方法,首先建立所要研究电网的典型运行方式数据;然后计算该运行方式下电网的加权潮流熵;接着逐一开断关键输电通道,比较各输电通道开断后的电网加权潮流熵水平,将引起熵值增加最大的输电通道作为下级故障,列入连锁故障序列,直至电网加权潮流熵大于等于0.5或潮流不收敛。本发明能够实现电网连锁故障集构建过程中关键支路的准确识别。

Description

一种基于最大熵增的连锁故障集构建方法
技术领域
本发明涉及电力系统领域,特别是一种基于最大熵增的连锁故障集构建方法。
背景技术
大停电事故的发生将给社会生产和居民生活造成重大损失,进行电力系统连锁故障传播、发展机理方面的深入研究,对电网规划、电网运行决策以及事故处理等具有重要意义。通过对近年来诸多国内外大停电事故共性特征上分析,一回 (或多回)线路故障导致的系统中大范围的潮流转移,使得并行输电通道上其他线路过载或保护误动而相继开断,一系列的连锁反应激发系统振荡并最终导致系统大面积停电甚至瓦解。例如2003年的美、加大停电就是由一条高压输电线路跳闸逐步蔓延开来的;2005年的莫斯科大停电则是从一个高压电流互感器爆炸开始的。由此可见连锁故障具有十分复杂的演变过程及发展路径,并随着电力系统规模及其复杂性、互联性的增加而加剧。准确预测连锁故障的发展过程,构建能够反映电网实际拓扑、运行特性的连锁故障故障集对于研究大停电事故的预防控制技术具有重大意义。
目前学者对大停电故障集构建技术的研究主要基于两大类方法,一类是基于时域仿真分析法,但其依赖人工经验设定大量故障,难以保证故障集的覆盖范围;另一类是基于复杂理论的分析,其优点是可以从宏观上把握电网发生大停电事故的风险,但仍存在无法具体筛选故障形式等问题。
发明内容
有鉴于此,本发明的目的是提出一种基于最大熵增的连锁故障集构建方法,能够实现电网连锁故障过程中关键支路的准确识别。
本发明采用以下方案实现:一种基于最大熵增的连锁故障集构建方法,具体包括以下步骤:
步骤S1:建立所要研究电网的典型运行方式数据;
步骤S2:计算该运行方式下电网的加权潮流熵;
步骤S3:逐一开断关键输电通道,比较各输电通道开断后的电网加权潮流熵水平,将引起熵值增加最大的输电通道作为下级故障,列入连锁故障序列;
步骤S4:判断此时的电网加权潮流熵是否小于预设的阈值并且潮流收敛,若是则保持当前作为下级故障的输电通道开断,并返回步骤S3;否则,进入步骤S5;
步骤S5:故障集构建过程结束。
进一步地,步骤S4中,所述预设的阈值为0.5。
进一步地,步骤S2中,加权潮流熵的计算包括以下步骤:
步骤S21:设线路i的传输稳定限额为Pmax,电网运行时线路i实时传输功率为 Pi0,则线路i的负载率μi0为:
Figure GDA0002688282370000021
式中,n为电网线路总数;
步骤S22:设线路负载率区间为[0,2.0],按照一定的步长X将负载率区间进行细分,形成负载率细分区间序列:
{[0,X],(X,2X]......((k-1)X,kX]......((m-1)X,2.0]},
其中,m=2/X,k=1,2......m;
步骤S23:设Lk为负载率μk0∈((k-1)X,kX]的线路条数,得到负载率μk0在区间段((k-1)X,kX]的线路占全网线路的概率:
Figure GDA0002688282370000022
步骤S24:建立潮流熵为:
Figure GDA0002688282370000023
式中,C1为常数,m为状态数;
步骤S25:计算加权潮流熵为:
Figure GDA0002688282370000031
式中,Wk为[(k-1)X,kX]区间线路的平均负载率。
进一步地,步骤S25中,平均负载率Wk的计算为:
Figure GDA0002688282370000032
式中,r为[(k-1)X,kX]区间中的线路数,Qk为第k条线路的负载率。
与现有技术相比,本发明有以下有益效果:本发明以加权熵来评价故障后电网所处状态及稳定风险,用最大熵增指标来识别电网关键支路和筛序后续故障。本发明的方法有助于深入研究大停电事故发展机理和故障集构建原则,并在此基础上研究应对大停电风险的电网应急控制方法及恢复技术,进而全面提升电网防御大停电事故能力和事故后恢复能力,防止连锁故障的蔓延,减少故障负荷损失,降低停电事故对社会经济的影响,有效预防大停电事故的发生。本发明能够准确识别电网薄弱环节。
附图说明
图1为本发明实施例的方法流程示意图。
图2为本发明实施例的网架构示意图。
图3为本发明实施例的停电规模与停电概率曲线。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
应该指出,以下详细说明都是示例性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和 /或它们的组合。
如图1所示,本实施例提供了一种基于最大熵增的连锁故障集构建方法,具体包括以下步骤:
步骤S1:建立所要研究电网的典型运行方式数据;
步骤S2:计算该运行方式下电网的加权潮流熵;
步骤S3:逐一开断关键输电通道,比较各输电通道开断后的电网加权潮流熵水平,将引起熵值增加最大的输电通道作为下级故障,列入连锁故障序列;
步骤S4:判断此时的电网加权潮流熵是否小于预设的阈值并且潮流收敛,若是则保持当前作为下级故障的输电通道开断,并返回步骤S3;否则,进入步骤S5;
步骤S5:故障集构建过程结束。
在本实施例中,步骤S4中,所述预设的阈值为0.5。
在本实施例中,步骤S2中,加权潮流熵的计算包括以下步骤:
步骤S21:设线路i的传输稳定限额为Pmax,电网运行时线路i实时传输功率为 Pi0,则线路i的负载率μi0为:
Figure GDA0002688282370000041
式中,n为电网线路总数;
步骤S22:综合考虑电网正常运行情况和可能遇到的故障后极端情况,并考虑一定的裕度,设线路负载率区间为[0,2.0],即线路可能空载,也可能达到负载率2.0 的极端状况;按照一定的步长X将负载率区间进行细分,形成负载率细分区间序列:
{[0,X],(X,2X]......((k-1)X,kX]......((m-1)X,2.0]},
其中,m=2/X,k=1,2......m;
例如,取X=5%,即可形成负载率区间的不同细分段,[0,0.05],(0.05, 0.10],(0.10,0.15]……(1.90,1.95],(1.95,2.00]。
步骤S23:设Lk为负载率μk0∈((k-1)X,kX]的线路条数,得到负载率μk0在区间段((k-1)X,kX]的线路占全网线路的概率:
Figure GDA0002688282370000051
步骤S24:建立潮流熵为:
Figure GDA0002688282370000052
式中,C1为常数,可根据研究需要设定,m为状态数;
步骤S25:考虑到线路负载率因素的影响,进一步提出了加权熵的概念,计算加权潮流熵为:
Figure GDA0002688282370000053
式中,Wk为[(k-1)X,kX]区间线路的平均负载率。
在本实施例中,步骤S25中,平均负载率Wk的计算为:
Figure GDA0002688282370000054
式中,r为[(k-1)X,kX]区间中的线路数,Qk为第k条线路的负载率。
具体的,本实施例以某地实际电网典型方式为例验证本实施例提出的连锁故障集构建及仿真方法,示例区域电网500及1000kV网架结构如图2所示。
首先对示例电网进行初始方式下加权潮流熵指标计算分析,加权潮流熵为0.288,小于设定的门槛值0.5,进行关键输电通道的逐一开断,并计算通道开断后的系统熵值,如下表所示:
开断通道 系统熵值
(26,27) 0.361
(30,31) 0.358
(33,34) 0.307
(34,35) 0.309
(1,9) 0.346
(2,3) 0.341
初始熵值为0.288,选出能够使系统熵值增加最大的通道26-27作为下级连锁故障支路,开断26-27后系统熵值未超过0.5,且潮流收敛,转步骤3,继续进行后续支路的逐一开断熵值分析。
重复上述步骤,得出连锁故障序列为:
1.开断通道26-27
2.开断通道30-31
开断通道26-27,30-31后系统潮流不收敛,故障集构建过程结束。
分析所提方法得出的故障集,通道26-27,30-31开断后,示例电网南北输电通道在东部的500KV层面联络被完全断开,此时南送电力全部通过西部10-15单回500KV线路和东部4回220KV线路联结。潮流大量转移至220KV线路引起相关区域电压下降,系统存在较高电压稳定风险。经过后续仿真分析,发生上述连锁故障后将会导致示例电网的电压失稳。
可以看出,本实施例提出的基于最大熵增的连锁故障集构建方法,通过在某地实际电网应用,验证了所提方法的可行性和有效性。
本实施例步骤S2通过计算系统(电网)加权熵指标,不仅能体现线路在各负载率区间的分布情况,还能综合考虑各区间线路的负载率大小,可以有效地表征潮流变化对电网自组织临界性的影响。采用加权潮流熵指标评估电网运行方式的潮流分布有序性,分析复杂电网自组织临界状态演化规律,对于从宏观上评估电网所处运行状态。
本实施例步骤S3提出的比较各关键支路(输电通道)开断后的电网加权潮流熵水平,将引起熵值增加最大的支路作为关键支路,列入连锁故障序列的故障集构建方法,能够筛选出对系统稳定性影响最大的输电支路,并将其列入连锁故障集。
特别的,本实施例以IEEE39节点算例为例,通过构建数十个运行方式并进行故障扫描计算,统计失稳故障负荷损失情况,能够证明在加权潮流熵较大的情况下,电网运行状态较差,负荷损失与停电累计概率呈现幂律特性,即电网进入自组织临界态,此时发生大停电的概率显著增大,其中,停电规模与停电概率曲线如图3所示。
可以看出,加权潮流熵符合电网实际情况,不仅能体现线路在各负载率区间的分布情况,还能综合考虑各区间线路的负载率大小,可以有效地表征潮流变化对电网自组织临界性的影响。可见,采用加权潮流熵指标评估电网运行方式的潮流分布有序性,分析复杂电网自组织临界状态演化规律,对于从宏观上评估电网所处运行状态以及发生大停电风险具有重要意义。
本实施例通过加权潮流熵指标判断电网所处状态。加权潮流熵指标能够定量评估系统中关键指标的有序稳定程度,通过IEEE39节点系统大量仿真计算表明当加权潮流熵增大到40-50%及以上时,电网易进入自组织临界态,即发生大停电的风险显著增大。对电网中重要输电通道分别做N-2故障,研究各个通道故障开断后的系统熵值,选出使系统熵增最大的支路作为连锁故障序列的后续故障形式,以系统熵值大于等于0.5作为故障及构建结束条件。
本实施例所提出的基于最大熵增的连锁故障集构建方法,以加权熵来评价故障后电网所处状态,以通道开断后引起最大熵增来识别电网关键支路和筛序后续故障的方法经验证能够适用于实际电网大停电的故障集构建,并以某地实际电网为例,验证了所提方法的有效性。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。

Claims (3)

1.一种基于最大熵增的连锁故障集构建方法,其特征在于:包括以下步骤:
步骤S1:建立所要研究电网的典型运行方式数据;
步骤S2:计算该运行方式下电网的加权潮流熵;
步骤S3:逐一开断关键输电通道,比较各输电通道开断后的电网加权潮流熵水平,将引起熵值增加最大的输电通道作为下级故障,列入连锁故障序列;
步骤S4:判断此时的电网加权潮流熵是否小于预设的阈值并且潮流收敛,若是则保持当前作为下级故障的输电通道开断,并返回步骤S3;否则,进入步骤S5;
步骤S5:故障集构建过程结束;
其中,步骤S2中,加权潮流熵的计算包括以下步骤:
步骤S21:设线路i的传输稳定限额为Pmax,电网运行时线路i实时传输功率为Pi0,则线路i的负载率μi0为:
Figure FDA0002688282360000011
式中,n为电网线路总数;
步骤S22:设线路负载率区间为[0,2.0],按照一定的步长X将负载率区间进行细分,形成负载率细分区间序列:
{[0,X],(X,2X]......((k-1)X,kX]......((m-1)X,2.0]},
其中,m=2/X,k=1,2......m;
步骤S23:设Lk为负载率μk0∈((k-1)X,kX]的线路条数,得到负载率μk0在区间段((k-1)X,kX]的线路占全网线路的概率:
Figure FDA0002688282360000012
步骤S24:建立潮流熵为:
Figure FDA0002688282360000013
式中,C1为常数,m为状态数;
步骤S25:计算加权潮流熵为:
Figure FDA0002688282360000021
式中,Wk为[(k-1)X,kX]区间线路的平均负载率。
2.根据权利要求1所述的一种基于最大熵增的连锁故障集构建方法,其特征在于:步骤S4中,所述预设的阈值为0.5。
3.根据权利要求1所述的一种基于最大熵增的连锁故障集构建方法,其特征在于:步骤S25中,平均负载率Wk的计算为:
Figure FDA0002688282360000022
式中,r为[(k-1)X,kX]区间中的线路数,Qk为第k条线路的负载率。
CN201811608477.2A 2018-12-27 2018-12-27 一种基于最大熵增的连锁故障集构建方法 Active CN109687410B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811608477.2A CN109687410B (zh) 2018-12-27 2018-12-27 一种基于最大熵增的连锁故障集构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811608477.2A CN109687410B (zh) 2018-12-27 2018-12-27 一种基于最大熵增的连锁故障集构建方法

Publications (2)

Publication Number Publication Date
CN109687410A CN109687410A (zh) 2019-04-26
CN109687410B true CN109687410B (zh) 2020-11-13

Family

ID=66188581

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811608477.2A Active CN109687410B (zh) 2018-12-27 2018-12-27 一种基于最大熵增的连锁故障集构建方法

Country Status (1)

Country Link
CN (1) CN109687410B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110556832B (zh) * 2019-08-22 2022-03-04 中国电力科学研究院有限公司 确定交直流电网发生连锁故障的影响因素的方法及系统
CN111082402B (zh) * 2019-12-31 2022-01-07 西安理工大学 一种关于输电网连锁故障序列的预测方法
CN111521908B (zh) * 2020-04-30 2021-02-12 华中科技大学 一种应用于四端风电直流电网的交流故障定位方法
CN111799773B (zh) * 2020-05-26 2022-05-20 中国电力科学研究院有限公司 一种用于预测电网连锁故障集的方法及系统
CN111817304B (zh) * 2020-05-29 2022-05-20 中国电力科学研究院有限公司 一种考虑降低连锁故障风险的电网断面调整方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107546734A (zh) * 2017-10-23 2018-01-05 国网天津市电力公司 一种电压型直流潮流控制器变比优化计算方法
CN107623319A (zh) * 2017-08-17 2018-01-23 广东电网有限责任公司惠州供电局 一种基于多评价指标的电网关键线路辨识方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107623319A (zh) * 2017-08-17 2018-01-23 广东电网有限责任公司惠州供电局 一种基于多评价指标的电网关键线路辨识方法
CN107546734A (zh) * 2017-10-23 2018-01-05 国网天津市电力公司 一种电压型直流潮流控制器变比优化计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于加权潮流熵和全动态仿真的故障集构建方法;霍启迪等;《2018 International Conference on Power System Technology》;20181108;正文第225-231页 *
基于联合加权熵的电网自组织临界状态演化;刘文颖等;《中国电机工程学报》;20150320;第35卷(第6期);全文 *

Also Published As

Publication number Publication date
CN109687410A (zh) 2019-04-26

Similar Documents

Publication Publication Date Title
CN109687410B (zh) 一种基于最大熵增的连锁故障集构建方法
Wang et al. The node degree distribution in power grid and its topology robustness under random and selective node removals
Ho et al. Optimal placement of fault indicators using the immune algorithm
CN110971525B (zh) 一种面向电力通信网的业务运行的业务路由寻址方法
CN110232642B (zh) 一种面向停电风险预防的拓扑规划与优化方法及装置
CN109449931B (zh) 结合加权潮流熵和全过程仿真的连锁故障集构建方法
CN108695846B (zh) 一种单元制配电网运行风险评估方法
CN107622360A (zh) 一种综合考虑主客观因素的关键线路识别方法
CN113328437B (zh) 一种智能配电网cps拓扑构建方法及故障恢复方法
CN112329376B (zh) 一种基于蒙特卡洛模拟的变电站系统抗震韧性量化评估方法
CN109494728B (zh) 一种配电终端选址选型方法
CN111882125B (zh) 基于风险元理论的电力信息物理系统连锁故障预测方法
CN114167217A (zh) 一种铁路配电网的多重故障诊断方法
CN106067074A (zh) 一种通过优化链路的开关状态来提升电网系统鲁棒性的方法
CN107808200A (zh) 配电网短时检修计划优化方法和装置
CN109861855B (zh) 一种电力通信网络中节点重要度的确定方法和装置
CN110021933B (zh) 考虑组件故障的电力信息系统控制功能可靠性评估方法
CN111191867B (zh) 电力系统复杂网络的可靠性评估方法
CN111126841A (zh) 一种基于随机潮流的配电网负荷转供能力评价方法
CN107994572B (zh) 一种配电网可靠性提升措施适用性量化分析方法
Fan et al. Risk assessment of power communication network based on node importance
CN112001073B (zh) 一种牵引供电系统的可靠性分析研究方法
CN112329377A (zh) 用于量化多输入多输出变电站系统地震易损性的评估方法
CN110829413B (zh) 基于道路矩阵的辐射状配电网可靠性评估方法
CN114519259A (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