CN110070286A - 电网多灾种耦合连锁故障分析方法及*** - Google Patents
电网多灾种耦合连锁故障分析方法及*** Download PDFInfo
- Publication number
- CN110070286A CN110070286A CN201910318624.0A CN201910318624A CN110070286A CN 110070286 A CN110070286 A CN 110070286A CN 201910318624 A CN201910318624 A CN 201910318624A CN 110070286 A CN110070286 A CN 110070286A
- Authority
- CN
- China
- Prior art keywords
- power grid
- grid
- power
- generator
- disaster
- 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.)
- Granted
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 42
- 230000008878 coupling Effects 0.000 title claims abstract description 32
- 238000010168 coupling process Methods 0.000 title claims abstract description 32
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 32
- 238000004088 simulation Methods 0.000 claims abstract description 35
- 239000011159 matrix material Substances 0.000 claims abstract description 8
- 230000005611 electricity Effects 0.000 claims description 24
- 230000005540 biological transmission Effects 0.000 claims description 20
- 238000000034 method Methods 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 15
- 230000001052 transient effect Effects 0.000 claims description 12
- 230000001788 irregular Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 230000007257 malfunction Effects 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 239000003990 capacitor Substances 0.000 claims description 5
- 230000007423 decrease Effects 0.000 claims description 5
- 230000005284 excitation Effects 0.000 claims description 5
- 230000006870 function Effects 0.000 claims description 5
- 238000002271 resection Methods 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 4
- 230000004888 barrier function Effects 0.000 claims description 4
- 230000001010 compromised effect Effects 0.000 claims description 4
- 238000010248 power generation Methods 0.000 claims description 4
- 230000001681 protective effect Effects 0.000 claims description 4
- 230000001012 protector Effects 0.000 claims description 4
- 238000011835 investigation Methods 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 230000001351 cycling effect Effects 0.000 claims description 2
- 230000033228 biological regulation Effects 0.000 claims 1
- 235000013399 edible fruits Nutrition 0.000 claims 1
- 230000008859 change Effects 0.000 description 4
- 208000025274 Lightning injury Diseases 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000003795 chemical substances by application Substances 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000000547 structure data Methods 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Computer Hardware Design (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Water Supply & Treatment (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Educational Administration (AREA)
Abstract
本发明公开了一种电网多灾种耦合连锁故障分析方法及***,包括:获取电网一次参数和电网二次参数,建立动态电网***模型;获取外部环境数据,输入动态电网***模型,判断如果存在外部电网故障事件,得到初始电网运行方式下由于外部环境引发电网元件故障的序列;建立电网混合微分代数方程;根据电网连锁故障仿真分析的初始条件和建立的电网混合微分代数方程对动态电网***模型进行仿真;仿真时,根据导纳矩阵Ybus判断电网是否解列,如果电网不解列,则继续运算;判断是否结束仿真;形成电网相继故障事故链;输出仿真数据,综合评价电网运行风险。本发明能考虑继电保护动作特性的多灾种耦合连锁故障对电网元件的影响。
Description
技术领域
本发明涉及电力***分析控制领域,尤其涉及电网多灾种耦合连锁故障分析方法及***。
背景技术
电网各种事故灾难呈现高频次发生的态势。大量案例表明,事故灾难的各种动力学过程之间、各动力学过程与环境之间不断进行物质和能量的交换,往往会形成灾害链,对灾难演变效应的不断放大起着关键作用。因此,面向重大事故灾难综合风险防范的总体目标,亟需由单一灾害的研究向灾害链复杂演变及其防治技术研究转变。
重大事故灾难往往包含着多灾种耦合与相互次生衍生动力学过程。重大事故灾难次生衍生及并发多灾种间的关系复杂多样,灾种间相互作用和时空组多灾种之间可能存在发生时间先后、影响范围叠加、致灾效果消长、相关或触发等复杂次生衍生及耦合关系。忽视这些关系会导致无法捕捉到多灾种并发作用下的事故灾难致灾机理。当前研究多采用单灾种风险直接相加或者赋权相加的方法来评估多灾种风险,然而对次生衍生及耦合作用的多灾种而言,综合风险并非简单的叠加,厘清灾害之间的相互作用关系和灾害链模式对多灾种风险评估极为重要。如何建立多灾种耦合关系的数学描述并在此基础上发展多灾种风险评估方法,是需要解决的重大科学问题。重大事故灾难次生衍生及并发多灾种间的关系复杂多样,灾种间相互作用和时空组合方式具有极强的非线性和不确定性,使用静态评估模型预测、评估重大事故灾难次生衍生和多灾害耦合的复杂***,会产生很大误差。
具体到电网而言,电网连锁故障分析方面,现有研究采用直流潮流方法计算电网的潮流分布,无法获得电压崩溃的节点;其次,未建立考虑一次拓扑和二次继电保护的电网连锁故障模型;现有方法均是基于历史电网元件故障或者预想事故集条件下的连锁故障分析,未考虑实际电网运行外部环境(如大范围山火、大范围冰灾、大范围雷击、频繁外力破坏等)的多种灾害耦合对电网元件的影响。最后,电网潮流与稳定分析计算量巨大,难以同时保证计算精度和计算准确率。
因此,迫切需要提出一种考虑继电保护动作特性的多灾种耦合连锁故障分析方法。
发明内容
本发明提供了一种电网多灾种耦合连锁故障分析方法及***,用以解决现有电网连锁故障分析未考虑多种灾害耦合对电网元件的影响的技术问题。
为解决上述技术问题,本发明提出的技术方案为:
一种电网多灾种耦合连锁故障分析方法,包括以下步骤:
获取电网一次参数和电网二次参数,根据电网一次参数和电网二次参数建立动态电网***模型;
获取外部环境数据,输入动态电网***模型,判断是否可能存在外部电网故障事件;如果存在外部电网故障事件,得到初始电网运行方式下由于外部环境引发电网元件故障的序列;
建立电网混合微分代数方程;根据电网连锁故障仿真分析的初始条件和建立的电网混合微分代数方程对动态电网***模型进行仿真;仿真时,计算获得输电线路及变压器元件是否由于电流、温度、节点电压以及频率参数和各元件设定的运行阈值而产生过流、超温、低压及低频的不正常运行状态,假如电网的元件存在不正常运行状态,则继电保护动作切除该元件;每发生一次切除元件的动作事件,根据导纳矩阵Ybus判断电网是否解列,如果电网不解列,则继续运算;如果电网解列,切除对应的支路,返回重新获取电网一次参数和电网二次参数;
根据仿真时间或电网混合微分代数方程的解算结果判断是否结束仿真;
根据由于外部环境引发电网元件故障的序列,形成电网相继故障事故链;
输出仿真数据,综合评价电网运行风险。
作为本发明的方法的进一步改进:
优选地,电网一次参数包括拓扑结构数据、线路参数、变压器参数、负荷参数、发电机参数、电容器参数和电抗器参数;电网二次参数包括:输电线路、发电机、母线继电保护装置配置及过流保护、距离保护、温度保护、低压保护以及低频保护的整定定值。
优选地,电网数据包括电网一次参数和电网二次参数;电网一次参数包括拓扑结构数据、线路参数、变压器参数、负荷参数、发电机参数、电容器参数和电抗器参数;电网二次参数包括:输电线路、发电机、母线继电保护装置配置及过流保护、距离保护、温度保护、低压保护以及低频保护的整定定值。
优选地,外部环境数据包括受威胁的输电线路的名称、受威胁的灾害因素及跳闸概率;
通过以下步骤得到初始电网运行方式下由于外部环境引发电网元件故障的序列:判断是否可能存在外部电网故障事件;如果存在外部电网故障事件,计算电网元件故障概率,当电网元件故障概率大于故障概率阈值时,判定元件故障,切除故障元件,并记录。
优选地,电网连锁故障仿真分析的初始条件,包括:电网仿真分析总时长Tend;电网潮流计算方式及收敛条件;电网仿真时间步长△t;过流保护动作时间Toc及动作方式Aoc;距离保护动作时间Tz及动作方式Az;温度保护动作时间Ttemp及动作方式Atemp;低压保护动作时间Tu及动作方式Au;以及低频保护动作时间Tf及动作方式Af。
优选地,动态电网***模型包括以下四种的任意组合:
(1)电网发电机模型:
有功功率和无功功率的非线性方程分别如下:
其中,i为发电机节点编号,|E′a,i|为暂态开环电压幅值,|Vi|为发电机节点的电压幅值,δm,i为发电机i的相角,X′d,i为直轴暂态电抗,Xq,i为交轴电抗,Xd,i为发电机直轴同步电抗;|E′a,i|为暂态开环电压幅值,其由以下微分方程决定:
T′do,i直轴暂态电抗时间常数,Efd,i是发电机励磁机的输出电压;
(2)电网转子设备模型:
其中,i为所有发电机节点集合中的任意节点,M为发电机的惯性常数;ωi为归一化的转子转速,D为抑制系数;Pm,i为发电机机械输入功率;Pg,i发电机的功率输出;
(3)电网励磁设备模型:
其中,Efd是发电机励磁机的输出电压,Vref为参考电压,Vt为节点实际电压,TA为自励绕组时间常数,TB为它励绕组时间常数,KE为励磁设备时间常数,sigm(·)为取值为[Emin,Emax]的S型函数;
(4)电网调速机模型:
调速机控制机械力从而使发电机产生一个转速偏差,本发明调速模型包括等比调速和二级代理调速控制,从而使***能保持频率偏差接近于0;
Pm为发电机输出功率,Pref为发电机参考输出功率,ωi为归一化的发电机转速,R为发电机允许在负荷突然增大或***电压下降情况下发电机的机端电压短暂的下降而不起强励或停机的时间常数,Ti为PI控制时间常数,Tt为伺服电机控制时间常数,P3为PI控制过程中的瞬时差分变量。
优选地,根据仿真时间和电网混合微分代数方程的解算结果判断是否结束仿真,包括:
判断仿真时间T是否达到设定总时长Tend,如果是,则终止仿真;如果不是则继续;
判断电网混合微分代数方程解算结果是否为h<0,如果h<0,则找出导致h<0的元件编号及对应的继电保护装置编号,切除该元件,并更新动态电网***模型。
优选地,仿真数据包括:时间t、发电机Ui的频率fui、发电机频率标幺值、发电机功率、发电机电势、母线电压、母线相角、励磁电压、励磁机终端电压、调压机功率和输电线路温度中的任意几种的组合。
优选地,以导纳矩阵Ybus判断电网是否解列为循环条件,包括如果电网不解列,则继续运算;如果电网解列,则返回读取电网数据步骤。
优选地,综合评价电网运行风险包括:根据电网相继故障事故链,结合事故链产生后的潮流分布计算获得电网支路切除率、电网负荷损失率和电网发电损失率三个指标综合判断电网运行风险,潮流分布包括:元件输送的功率和节点的功率。
本发明还提供一种电网多灾种耦合连锁故障分析***,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现上述任一方法的步骤。
本发明具有以下有益效果:
本发明的一种电网多灾种耦合连锁故障分析方法及***,综合考虑电网运行方式、继电保护动作特性和外部环境对电网影响,可准确仿真外部大范围山火、冰灾、雷击等灾害引发电网故障后对电网连锁故障链发展的分析。在优选方案中,结合了电网一次电网运行方式和二次继电保护动作特性,实现了电网连锁故障的高精度和高准确率仿真。根据电网连锁故障事故链结论,可及时做好相应的应急处置预案,实现多灾种耦合条件下的提前应对,减少大范围灾害耦合对电网所造成的损失。准确率高,实用性强。
除了上面所描述的目的、特征和优点之外,本发明还有其它的目的、特征和优点。下面将参照附图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本发明优选实施例的电网多灾种耦合连锁故障分析方法的流程示意图。
具体实施方式
以下结合附图对本发明的实施例进行详细说明,但是本发明可以由权利要求限定和覆盖的多种不同方式实施。
参见图1,本发明的一种电网多灾种耦合连锁故障分析方法,包括以下步骤:
获取电网一次参数和电网二次参数,根据电网一次参数和电网二次参数建立动态电网***模型;
获取外部环境数据,输入动态电网***模型,判断是否可能存在外部电网故障事件;如果存在外部电网故障事件,得到初始电网运行方式下由于外部环境引发电网元件故障的序列;
建立电网混合微分代数方程;根据电网连锁故障仿真分析的初始条件和建立的电网混合微分代数方程对动态电网***模型进行仿真;仿真时,计算获得输电线路及变压器元件是否由于电流、温度、节点电压以及频率参数和各元件设定的运行阈值而产生过流、超温、低压及低频的不正常运行状态,假如电网的元件存在不正常运行状态,则继电保护动作切除该元件;每发生一次切除元件的动作事件,根据导纳矩阵Ybus判断电网是否解列,如果电网不解列,则继续运算;如果电网解列,切除对应的支路,返回重新获取电网一次参数和电网二次参数;
根据仿真时间或电网混合微分代数方程的解算结果判断是否结束仿真;
根据由于外部环境引发电网元件故障的序列,形成电网相继故障事故链;
输出仿真数据,综合评价电网运行风险。
上述步骤,综合考虑电网运行方式、继电保护动作特性和外部环境对电网影响,可准确仿真外部大范围山火、冰灾、雷击等灾害引发电网故障后对电网连锁故障链发展的分析。
实际实施时,以上的方法还能进行以下的扩充或应用,以下实施例中的技术特征都能相互组合,实施例仅作为示例,不作为对技术特征的正常组合限制。
实施例1:
参见图1,本实施例的一种电网多灾种耦合连锁故障分析方法,包括以下步骤:
S1:设置电网连锁故障仿真分析初始条件;
电网连锁故障仿真分析初始条件包括:电网仿真分析总时长Tend;电网潮流计算方式及收敛条件;电网仿真时间步长△t;过流保护动作时间Toc及动作方式Aoc;距离保护动作时间Tz及动作方式Az;温度保护动作时间Ttemp及动作方式Atemp;低压保护动作时间Tu及动作方式Au;以及低频保护动作时间Tf及动作方式Af。
S2:读取电网数据。
2.1读取电网一次参数包括:拓扑结构数据、线路参数、变压器参数、负荷参数、发电机参数、电容器参数和电抗器参数。
2.2读取电网二次参数包括:输电线路、发电机、母线继电保护装置配置及过流保护、距离保护、温度保护、低压保护以及低频保护的整定定值。
S3:建立电网混合微分代数方程。
由一系列的微分代数方程表示动态电网***,本发明考虑到继电保护装置产生的离散事件,增加了离散方程,从而产生混合微分代数方程。
假如电网***在t时刻的状态可以由x(t)、y(t)、z(t)三个矢量描述:
其中:x(t)为连续状态矢量,按照如下微分方程随着时间变化。
y(t)为连续状态变量,它与其他变量为纯粹的数学关系。
g(t,x(t),y(t),z(t))=0 (2)
z(t)为状态矢量,它仅仅在整数状态。
h(t,x(t),y(t),z(t))<0 (3)
当约束hi(...)<0不满足时,相关的计数器di激活。其中,hi(...)<0为判断电网***是否解算成功的依据,当hi(…)>0时,则计数器di激活,当多次hi(…)>0时,di计数会达到其极限值,则认为该次解算潮流不收敛;改变z(t)的值,di清零,重新开始仿真计算。直至z(t)改变的次数达到初始条件设置的收敛仿真次数,假如最后一次还不收敛,则认为潮流无法收敛,停止计算。对于复杂电网***而言,需要计算的参量特别多,当功率平衡或初始值选择不正确都可能造成这种情况。
此时则判断解算失败,停止解算,分析无结果。对于复杂电网***而言,需要计算的参量特别多,当功率平衡或初始值选择不正确都可能造成解算失败的情况。
本实施例实施时,动态电网***模型包括下面的电网元件模型,从而得到上述的电网混合微分代数方程。以下模型可根据实际的电网仿真的精度和时间要求选取,本实施例为提高仿真精度,选择同时使用下述的四种模型。
所述电网发电机模型:
采用发电机凸子和极子模型,其有功功率和无功功率的非线性方程:
i为发电机节点编号,|E′a,i|为暂态开环电压幅值,|Vi|为发电机节点的电压幅值,δm,i为发电机i的相角,X′d,i为直轴暂态电抗,Xq,i为交轴电抗,Xd,i为发电机直轴同步电抗。|E′a,i|为暂态开环电压幅值,其由以下微分方程决定:
T′do,i直轴暂态电抗时间常数,Efd,i是发电机励磁机的输出电压。
所述电网转子设备模型:
转子转速与母线i的关联关系用典型的二阶转子方程表示。用于描述同步发电机的转动特性。标准的转子转速w在正常运行状态下是固定的,当产生有振荡的时候,w会加速或者减速。
i为所有发电机节点集合中的任意节点,M为发电机的惯性常数;ωi为归一化的转子转速,D为抑制系数;Pm,i为发电机机械输入功率;Pg,i发电机的功率输出;
所述电网励磁设备模型:
发电机励磁***的利用一般的二阶微分方程模型。
Efd是发电机励磁机的输出电压,Vref为参考电压,Vt为节点实际电压,TA为自励绕组时间常数,TB为它励绕组时间常数,KE为励磁设备时间常数,sigm(·)为取值为[Emin,Emax]的S型函数。
所述电网调速机模型:
调速机控制机械力从而使发电机产生一个转速偏差,本发明调速模型包括等比调速和二级代理调速控制,从而使***能保持频率偏差接近于0。
Pm为发电机输出功率,Pref为发电机参考输出功率,ωi为归一化的发电机转速,R为发电机允许在负荷突然增大或***电压下降等情况下发电机的机端电压短暂的下降而不起强励或停机的时间常数,Ti为PI控制时间常数,Tt为伺服电机控制时间常数,P3为PI控制过程中的瞬时差分变量。
S4:读取外部环境数据,判断是否可能存在外部电网故障事件;
4.1如果存在外部电网故障事件,计算电网元件故障概率,如跳闸概率超过0.8,则认为元件故障。
S5:根据电网连锁故障仿真分析的初始条件和建立的电网混合微分代数方程对动态电网***模型进行仿真;仿真时,计算获得输电线路、变压器等元件是否由于电流、温度、节点电压以及频率参数,再基于各元件设定的运行阈值而产生过流、超温、低压、低频等的不正常运行状态,假如电网的元件存在所述不正常运行状态,则继电保护动作切除该元件;每发生一次切除元件的动作事件。
S6:根据导纳矩阵Ybus判断电网是否解列;
6.1如果电网不解列,则继续运算;
6.2如果电网解列,则切除对应的支路,返回重新获取电网一次参数和电网二次参数;进入步骤S2。
S7:根据仿真时间或电网混合微分代数方程的解算结果判断是否结束仿真。
解算电网混合微分代数方程根据梯形模型。当在时间td发生跳闸事件时,解算以下方程式:
g(td,x+,y+,z)=0 (5)
h(td,x+,y+)=0 (6)
d(td,x+,y+)=0 (7)
式中,t为前一步长的时间点,即td=t+△td,d是计数函数,x、y、z是公式(1)、(2)、(3)中x(t)、y(t)、z(t)的简写。x+、y+为x、y在t时刻下一个瞬间的初值。
7.1判断仿真时间T是否达到设定总时长Tend,如果是,则终止仿真;如果不是则继续。需要每计算完成一个步长就会判断时间。另外,计算完一个步长就会对应有一个结果,当仿真结果对应有元件发生连锁故障时,需要先修改电网拓扑,然后再继续仿真,即使没有发生元件连锁故障,还需要判断是否外部环境中是否有故障输入。
7.2判断电网混合微分代数方程解算结果是否为h<0,如果h<0,则找出导致h<0的元件编号及对应的继电保护装置编号,切除该元件,即引发的连锁故障元件,并进入步骤S3;如果h≥0,则进入步骤S4。
S8:根据所述由于外部环境引发电网元件故障的序列,形成电网相继故障事故链。记录在初始电网运行方式下由于外部环境引发电网元件故障的序列和外部故障引发的电网由于过流、超温、低压、低频等原因导致继电保护动作的切断元件的序列,形成电网相继故障事故链。
S9:输出仿真数据到excel,综合评价电网运行风险。
仿真数据包括:时间t、发电机Ui的频率fui、发电机频率标幺值、发电机功率、发电机电势、母线电压、母线相角、励磁电压、励磁机终端电压、调压机功率和输电线路温度几种的组合。根据电网支路切除率、电网负荷损失率、电网发电损失率三个指标判断电网运行风险。
步骤S9中,电网支路切除率表征在连锁故障结束时保持有源网络的支路与初始电网支路数的比值。
其中,ls为在连锁故障的最后一层时独立子电网***中有源网络的支路数;n为连锁故障的最后一层时独立子电网***个数。li为初始电网支路数。ηSl值总为非负值,值越大表示电网***过载程度越小。
电网发电损失率:
其中Cl,i为在连锁故障的结束时独立子电网***中发电机的有功输出功率,k为最后一层时连锁故障独立子电网***的数量;n为初始独立子电网***个数,Cl为初始电网发电机输出功率总数。ηSc值总为非负值,值越大表示电网***中输电容量保持程度越高。
电网负荷损失率:
对于连锁故障的结束时每一个电网子***,获得每一个电网子***中保持负荷量Ls。
其中m为最后一层时连锁故障独立子电网***集合Send的数量;Pjk,gen为集合Send中第j个独立电网子***Gend,j中第k个发电机的输出功率有功,Djl,bus为集合Send中第j个独立电网子***Gend,j中第l个节点的负荷有功。ηSd值总为非负值,值越大表示电网***中负荷需求保持程度越高。
实施例2:
本实施例为实施例1的应用例,包括以下步骤:
S1:读取电网数据,包括电网一次参数中1711条线路的参数、1323条母线参数、121台电机(含励磁和控制器)、233个电容和电感、667个负荷。
S2:读取电网二次参数包括:输电线路、发电机、母线继电保护装置配置及过流保护、距离保护、温度保护、低压保护以及低频保护的整定定值;
S3:建立电网混合微分代数方程。
S4:设置电网连锁故障仿真分析初始条件;
电网仿真分析总时长Tend为2小时;
电网潮流计算方式为牛顿拉弗逊潮流,收敛条件为误差小于等于10-9;
电网仿真时间步长△t=0.5s;
过流保护动作时间Toc=0.5s及动作方式Aoc为切除线路;
距离保护动作时间Tz=0.5s及动作方式Az为切除线路;
温度保护动作时间Ttemp=0s及动作方式Atemp为切除线路;
低压保护动作时间Tu=0.5s及动作方式Au为切除所在母线25%的负荷;
低频保护动作时间Tf=0.5s及动作方式Af为切除所在母线25%的负荷。
S5:读取外部环境数据,在t=2800s时电网发现2个输电线路山火Fire1和Fire2,其中Fire1威胁500kV船苏I线,其跳闸概率超过0.92,认为在t=2800s时500kV船苏I线跳闸;
S6:根据导纳矩阵Ybus判断电网未发生解列,步骤7、判断电网混合微分代数方程解算结果为h<0,220kV输电线路耒真I线#1继电保护的过流保护动作,经查,该线路此时的输送电流为546A,超过了其输电容量380A。
S7:输出仿真数据到excel。
仿真数据包括:时间t、发电机Ui的频率fui、发电机频率标幺值、发电机功率、发电机电势、母线电压、母线相角、励磁电压、励磁机终端电压、调压机功率、输电线路温度。
S8:记录在初始电网运行方式下由于外部环境引发电网元件故障的序列和外部故障引发的电网由于过流、超温、低压、低频等原因导致继电保护动作的切断元件的序列,形成电网相继故障事故链。
S9:根据电网支路切除率、电网负荷损失率、电网发电损失率三个指标判断电网运行风险。
在上述的实施例中,本发明通过以下两个要点实现,电网的多灾种耦合分析,即多种灾害相继作用于电网的元件,其对电网运行产生的影响通过叠加而不断扩大,比独立的电网灾害(单一灾种)对电网影响大得多。
1)通过电网故障事件的输入,即在一次性输入多个外部故障或者在短时间内输入多个外部故障,即体现了多灾种耦合;
2)在运算过程中,通过计算电网各节点、各元件的运行参数时,通过微分代数方程考虑了电网多种灾害的故障叠加的影响。
实施例3:
本实施例提供一种电网多灾种耦合连锁故障分析***,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现上述任一实施例的步骤。
综上可知,本发明通过综合电网运行方式、继电保护动作特性及外部环境的跳闸概率等多因素的多灾种耦合连锁故障分析方法。针对电网实际跳闸的继电保护类型,建立继电保护模型,实现电网元件的故障模拟;基于电网外部环境的发展过程,实时引入外部多种灾害对电网的故障概率;构建电网元件物理模型的电网拓扑及特性,提高电网仿真的准确率,采用变步长的电网仿真分析方法,提高了仿真效率。该方法思路新颖、流程清晰、准确率高和实用性强。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种电网多灾种耦合连锁故障分析方法,其特征在于,包括以下步骤:
获取电网一次参数和电网二次参数,根据所述电网一次参数和电网二次参数建立动态电网***模型;
建立电网混合微分代数方程;
获取外部环境数据,输入所述动态电网***模型,判断是否可能存在外部电网故障事件;如果存在外部电网故障事件,得到初始电网运行方式下由于外部环境引发电网元件故障的序列;
根据电网连锁故障仿真分析的初始条件和建立的电网混合微分代数方程对动态电网***模型进行仿真;仿真时,计算获得输电线路及变压器元件是否由于电流、温度、节点电压以及频率参数和各元件设定的运行阈值而产生过流、超温、低压及低频的不正常运行状态,假如电网的元件存在所述不正常运行状态,则继电保护动作切除该元件;每发生一次切除元件的动作事件,根据导纳矩阵Ybus判断电网是否解列,如果电网不解列,则继续运算;如果电网解列,切除对应的支路,返回重新获取电网一次参数和电网二次参数;
根据仿真时间或电网混合微分代数方程的解算结果判断是否结束仿真;
根据所述由于外部环境引发电网元件故障的序列,形成电网相继故障事故链;
输出仿真数据,综合评价电网运行风险。
2.根据权利要求1所述的电网多灾种耦合连锁故障分析方法,其特征在于,所述电网一次参数包括拓扑结构数据、线路参数、变压器参数、负荷参数、发电机参数、电容器参数和电抗器参数;所述电网二次参数包括:输电线路、发电机、母线继电保护装置配置及过流保护、距离保护、温度保护、低压保护以及低频保护的整定定值。
3.根据权利要求1所述的电网多灾种耦合连锁故障分析方法,其特征在于,所述外部环境数据包括受威胁的输电线路的名称、受威胁的灾害因素及跳闸概率;
通过以下步骤得到初始电网运行方式下由于外部环境引发电网元件故障的序列:判断是否可能存在外部电网故障事件;如果存在外部电网故障事件,计算电网元件故障概率,当所述电网元件故障概率大于故障概率阈值时,判定元件故障,切除故障元件,并记录。
4.根据权利要求1所述的电网多灾种耦合连锁故障分析方法,其特征在于,所述电网连锁故障仿真分析的初始条件,包括:电网仿真分析总时长Tend;电网潮流计算方式及收敛条件;电网仿真时间步长△t;过流保护动作时间Toc及动作方式Aoc;距离保护动作时间Tz及动作方式Az;温度保护动作时间Ttemp及动作方式Atemp;低压保护动作时间Tu及动作方式Au;以及低频保护动作时间Tf及动作方式Af。
5.根据权利要求1至4中任一项所述的电网多灾种耦合连锁故障分析方法,其特征在于,所述动态电网***模型包括以下四种的任意组合:
(1)电网发电机模型:
有功功率和无功功率的非线性方程分别如下:
其中,i为发电机节点编号,|E′a,i|为暂态开环电压幅值,|Vi|为发电机节点的电压幅值,δm,i为发电机i的相角,X′d,i为直轴暂态电抗,Xq,i为交轴电抗,Xd,i为发电机直轴同步电抗;|E′a,i|为暂态开环电压幅值,其由以下微分方程决定:
T′do,i直轴暂态电抗时间常数,Efd,i是发电机励磁机的输出电压;
(2)电网转子设备模型:
其中,i为所有发电机节点集合中的任意节点,M为发电机的惯性常数;ωi为归一化的转子转速,D为抑制系数;Pm,i为发电机机械输入功率;Pg,i发电机的功率输出;
(3)电网励磁设备模型:
其中,Efd是发电机励磁机的输出电压,Vref为参考电压,Vt为节点实际电压,TA为自励绕组时间常数,TB为它励绕组时间常数,KE为励磁设备时间常数,sigm(·)为取值为[Emin,Emax]的S型函数;
(4)电网调速机模型:
调速机控制机械力从而使发电机产生一个转速偏差,本发明调速模型包括等比调速和二级代理调速控制,从而使***能保持频率偏差接近于0;
Pm为发电机输出功率,Pref为发电机参考输出功率,ωi为归一化的发电机转速,R为发电机允许在负荷突然增大或***电压下降情况下发电机的机端电压短暂的下降而不起强励或停机的时间常数,Ti为PI控制时间常数,Tt为伺服电机控制时间常数,P3为PI控制过程中的瞬时差分变量。
6.根据权利要求5所述的电网多灾种耦合连锁故障分析方法,其特征在于,根据仿真时间和电网混合微分代数方程的解算结果判断是否结束仿真,包括:
判断仿真时间T是否达到设定总时长Tend,如果是,则终止仿真;如果不是则继续;
判断电网混合微分代数方程解算结果是否为h<0,如果h<0,则找出导致h<0的元件编号及对应的继电保护装置编号,切除该元件,并更新动态电网***模型。
7.根据权利要求6所述的电网多灾种耦合连锁故障分析方法,其特征在于,所述仿真数据包括:时间t、发电机Ui的频率fui、发电机频率标幺值、发电机功率、发电机电势、母线电压、母线相角、励磁电压、励磁机终端电压、调压机功率和输电线路温度中的任意几种的组合。
8.根据权利要求6所述的电网多灾种耦合连锁故障分析方法,其特征在于,以导纳矩阵Ybus判断电网是否解列为循环条件,包括如果电网不解列,则继续运算;如果电网解列,则返回读取电网数据步骤。
9.根据权利要求8所述的电网多灾种耦合连锁故障分析方法,其特征在于,综合评价电网运行风险包括:根据所述电网相继故障事故链,结合事故链产生后的潮流分布计算获得电网支路切除率、电网负荷损失率和电网发电损失率三个指标综合判断电网运行风险,所述潮流分布包括:元件输送的功率和节点的功率。
10.一种电网多灾种耦合连锁故障分析***,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现上述权利要求1至9中任一所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910318624.0A CN110070286B (zh) | 2019-04-19 | 2019-04-19 | 电网多灾种耦合连锁故障分析方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910318624.0A CN110070286B (zh) | 2019-04-19 | 2019-04-19 | 电网多灾种耦合连锁故障分析方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110070286A true CN110070286A (zh) | 2019-07-30 |
CN110070286B CN110070286B (zh) | 2022-05-27 |
Family
ID=67368122
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910318624.0A Active CN110070286B (zh) | 2019-04-19 | 2019-04-19 | 电网多灾种耦合连锁故障分析方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110070286B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110489729A (zh) * | 2019-08-08 | 2019-11-22 | 国网湖南省电力有限公司 | D5000-matpower的电网多灾种耦合连锁故障模型的自动转换方法及*** |
CN110675043A (zh) * | 2019-09-17 | 2020-01-10 | 深圳供电局有限公司 | 基于连锁故障模型确定电网停电关键线路的方法及*** |
CN111488675A (zh) * | 2020-03-18 | 2020-08-04 | 四川大学 | 电力***连锁故障潜在触发模式的挖掘方法 |
CN111641360A (zh) * | 2020-06-11 | 2020-09-08 | 上海外高桥造船有限公司 | 一种瞬态压降获取、设备起动的方法和装置、设备和介质 |
CN111798112A (zh) * | 2020-06-28 | 2020-10-20 | 中国矿业大学(北京) | 一种城市社区社会公共安全综合风险评估方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105701288A (zh) * | 2016-01-12 | 2016-06-22 | 国网山东省电力公司电力科学研究院 | 一种极端灾害气象条件下电网复杂相继故障模拟与仿真方法 |
CN106356846A (zh) * | 2016-10-18 | 2017-01-25 | 国网山东省电力公司烟台供电公司 | 一种基于时间的初期电网连锁故障模拟仿真方法 |
CN106548265A (zh) * | 2016-09-07 | 2017-03-29 | 清华大学深圳研究生院 | 一种基于连锁故障事故链搜索的输电网可靠性评估方法 |
CN106655167A (zh) * | 2016-12-06 | 2017-05-10 | 国网北京市电力公司 | 电网连锁故障的仿真方法及装置 |
EP3299828A1 (en) * | 2016-09-21 | 2018-03-28 | Rolls-Royce plc | Electrical fault detection |
-
2019
- 2019-04-19 CN CN201910318624.0A patent/CN110070286B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105701288A (zh) * | 2016-01-12 | 2016-06-22 | 国网山东省电力公司电力科学研究院 | 一种极端灾害气象条件下电网复杂相继故障模拟与仿真方法 |
CN106548265A (zh) * | 2016-09-07 | 2017-03-29 | 清华大学深圳研究生院 | 一种基于连锁故障事故链搜索的输电网可靠性评估方法 |
EP3299828A1 (en) * | 2016-09-21 | 2018-03-28 | Rolls-Royce plc | Electrical fault detection |
CN106356846A (zh) * | 2016-10-18 | 2017-01-25 | 国网山东省电力公司烟台供电公司 | 一种基于时间的初期电网连锁故障模拟仿真方法 |
CN106655167A (zh) * | 2016-12-06 | 2017-05-10 | 国网北京市电力公司 | 电网连锁故障的仿真方法及装置 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110489729A (zh) * | 2019-08-08 | 2019-11-22 | 国网湖南省电力有限公司 | D5000-matpower的电网多灾种耦合连锁故障模型的自动转换方法及*** |
CN110489729B (zh) * | 2019-08-08 | 2022-12-06 | 国网湖南省电力有限公司 | D5000-matpower的电网多灾种耦合连锁故障模型的自动转换方法及*** |
CN110675043A (zh) * | 2019-09-17 | 2020-01-10 | 深圳供电局有限公司 | 基于连锁故障模型确定电网停电关键线路的方法及*** |
CN111488675A (zh) * | 2020-03-18 | 2020-08-04 | 四川大学 | 电力***连锁故障潜在触发模式的挖掘方法 |
CN111488675B (zh) * | 2020-03-18 | 2023-04-25 | 四川大学 | 电力***连锁故障潜在触发模式的挖掘方法 |
CN111641360A (zh) * | 2020-06-11 | 2020-09-08 | 上海外高桥造船有限公司 | 一种瞬态压降获取、设备起动的方法和装置、设备和介质 |
CN111798112A (zh) * | 2020-06-28 | 2020-10-20 | 中国矿业大学(北京) | 一种城市社区社会公共安全综合风险评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110070286B (zh) | 2022-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110070286A (zh) | 电网多灾种耦合连锁故障分析方法及*** | |
Kundur | Power system stability | |
Singh et al. | Decision tree based anomaly detection for remedial action scheme in smart grid using PMU data | |
US20140032138A1 (en) | Apparatus and method for out-of-step protection using the analysis of trajectories of electrical measurements in state plane | |
CN106159945B (zh) | 一种基于多场景模式下的风电连锁脱网路径链搜寻方法 | |
CN104978482B (zh) | 基于psd‑bpa批处理变压器故障累积效应评估方法 | |
Adhikari et al. | Development of power system test bed for data mining of synchrophasors data, cyber-attack and relay testing in RTDS | |
Anagnostou et al. | Impact of overexcitation limiters on the power system stability margin under stressed conditions | |
Mu et al. | The mechanism of DFIGs grouping tripped off from power grid | |
Wang et al. | Adaptive corrective control strategies for preventing power system blackouts | |
Ramasubramanian et al. | Viability of synchronous reference frame phase locked loop inverter control in an all inverter grid | |
CN105701288B (zh) | 一种极端灾害气象条件下电网复杂相继故障模拟与仿真方法 | |
Liu et al. | An adaptive inverse time overcurrent relay model implementation for real time simulation and hardware-in-the-loop testing | |
Lu et al. | Validation of real-time system model in western interconnection | |
CN109546650A (zh) | 输电通道暂稳极限的估算方法及装置 | |
CN106451441B (zh) | 一种电力***N-k故障集的快速筛选方法 | |
Ruhle et al. | Ranking of system contingencies in DSA systems-first experiences | |
Zhu et al. | Identifying critical load locations for power system voltage, angular and frequency stability | |
Rovnyak et al. | Real-time transient stability prediction-possibilities for on-line automatic database generation and classifier training | |
Clark et al. | Development of a Base Model in RMS and EMT Environment to Study Low Inertia System | |
Kim et al. | Vulnerability indices for power systems | |
Mola-Jimenez et al. | PowerFactory-Python based assessment of frequency and transient stability in power systems dominated by power electronic interfaced generation | |
Rovnyak et al. | Performance index and classifier approaches to real-time, discrete-event control | |
Izuegbunam et al. | Dynamic security assessment of 330kV Nigeria power system | |
Yuan et al. | Simulation model and characteristics of variable frequency transformers used for grid interconnection |
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 |