CN102522759B - 一种基于扰动后电压轨迹的电力系统暂态失稳实时判别方法 - Google Patents
一种基于扰动后电压轨迹的电力系统暂态失稳实时判别方法 Download PDFInfo
- Publication number
- CN102522759B CN102522759B CN 201110417168 CN201110417168A CN102522759B CN 102522759 B CN102522759 B CN 102522759B CN 201110417168 CN201110417168 CN 201110417168 CN 201110417168 A CN201110417168 A CN 201110417168A CN 102522759 B CN102522759 B CN 102522759B
- Authority
- CN
- China
- Prior art keywords
- voltage
- tenesmus
- power system
- value
- system transient
- 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
- 230000001052 transient effect Effects 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 title claims abstract description 84
- 238000005259 measurement Methods 0.000 claims abstract description 21
- 150000001875 compounds Chemical class 0.000 claims description 63
- 230000010354 integration Effects 0.000 claims description 59
- 206010057071 Rectal tenesmus Diseases 0.000 claims description 57
- 208000012271 tenesmus Diseases 0.000 claims description 57
- 230000008569 process Effects 0.000 claims description 36
- 238000005070 sampling Methods 0.000 claims description 17
- 230000005540 biological transmission Effects 0.000 claims description 15
- 230000004069 differentiation Effects 0.000 claims description 14
- 230000004044 response Effects 0.000 claims description 11
- 238000004088 simulation Methods 0.000 claims description 7
- 238000012544 monitoring process Methods 0.000 claims description 3
- 230000035772 mutation Effects 0.000 claims description 3
- 230000006641 stabilisation Effects 0.000 claims description 2
- 238000011105 stabilization Methods 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 8
- 230000007123 defense Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 230000009471 action Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 230000033228 biological regulation Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及电力系统故障判别技术领域,具体涉及一种基于扰动后电压轨迹的电力系统暂态失稳实时判别方法,包括遴选电压观测节点、故障清除后电压下坠识别、电压轨迹复合积分、系统暂态稳定性判别、局部电压过低识别五个步骤。本发明基于电网中现有的相量测量单元(PMU)量测配置以及广域量测系统(WAMS),着力解决二道防线控制范围以外的大扰动造成的系统暂态失稳。本发明算法简单,计算量小,对数据精度依赖性低,可靠性高。
Description
技术领域
本发明涉及电力系统故障判别技术领域,具体涉及一种基于扰动后电压轨迹的电力系统暂态失稳实时判别方法。
背景技术
随着电网的互联,电力系统规模日益增大,系统的动态行为日益复杂,局部的故障影响范围增大,可能带来连锁性的故障,进而导致大面积的停电。电力系统现有的基于事件发生的常规控制策略在适应性、预测速度、控制效率、可靠性方面还存在很多不足。GPS技术和同步相量测量单元(PMU)在电力系统中的应用为实现基于实测轨迹的电力系统暂态失稳预警创造了条件。
电力系统暂态失稳实时判别是指通过分析电力系统各状态量在扰动后的受扰轨迹,快速判断系统在当前的动态行为下是否能够恢复到稳定。如果判定系统将失稳,则可以及时采取切机、切负荷、解列等紧急控制措施。当前暂态失稳的判别,大多基于发电机功角受扰轨迹的分析。文献一《Fastgeneration shedding equipment based on the observation of swings of generators》(IEEE Transactions on Power Systems,1988年第3卷第2期第439页)提出以系统中任意两台发电机间的功角差超过一个门槛值作为失稳判据。算法简单直接,容易操作。但是比较粗糙,对于非自治系统门槛值会随系统运行状态而变化,难以适应复杂的系统情况。文献二《基于轨迹几何特征的暂态不稳定识别研究》(中国电机工程学报2008年第28卷第4期第16页)提出了基于相平面轨迹几何特征的暂态不稳定判别方法,此类算法属于微分类算法,虽能快速判别系统暂态稳定性,但对数据质量要求高,数据误差有可能会导致误判。文献三《实测摇摆曲线的暂态稳定量化分析》(电力系统自动化2004年第28卷第20期第17页)提出基于扩展等面积法则(EEAC)的暂态稳定实时判别方法,但是它依赖于发电机分群结果以及对数据量要求较大,对于基于实测轨迹的大电网的暂态失稳判别而言较为困难。
因此,传统的基于发电机功角受扰轨迹的暂态稳定性分析存在一定的局限性,且工程实用性较差。文献四《WACS-Wide-Area Stability and VoltageControl System:R&D and Online Demonstration》(Proceedings of The IEEE,2005年第93卷第5期第892页)提出了以电压作为观测信息,利用直接积分运算分析电力系统暂态稳定性的新方法。为基于实测轨迹的电力系统暂态失稳预警提供了新的视野。但是本方法采用直接积分算法,若电压快速下坠时会造成漏判,危机电网的安全稳定,且原文中针对美国实际电网设置布点及判别方案,仅提供了开放性思路,并未给出完整的运行流程,尚存在许多问题需要解决。
发明内容
本发明解决的技术问题是克服现有技术的不足,提供一种算法简单、计算量小、对数据精度依赖性低且可靠性高的基于扰动后电压轨迹的暂态失稳实时判别方法。
常规的电力系统第二道防线采用的是基于事件的稳定控制技术,它的优点是针对性强、可靠性高,但是若系统发生的故障超出了预先设定的事件范围,将不会采取控制措施,故障严重时则会发生大面积停电事故。本发明是基于电网中现有的相量测量单元(PMU)量测配置以及广域量测系统(WAMS),着力解决二道防线控制范围以外的大扰动造成的系统暂态失稳,当电网发生大的扰动后,根据PMU上传的变电站电压量测信息,自动识别电压在故障清除后的下坠轨迹,快速挖掘轨迹所蕴含的系统稳定信息,计算得到暂态稳定指标,超前地判别系统的暂态稳定性。
为解决上述技术问题,本发明采用的技术方案如下:
一种基于扰动后电压轨迹的电力系统暂态失稳实时判别方法,包括以下步骤:
步骤一,遴选电压观测节点:按照电网的实际运行方式,选取真实反映系统暂态稳定特性的重要变电站,将其电压作为观测信息;
步骤二,故障清除后电压下坠识别:以步骤一中得到的节点作为实际电压观测点,在电网运行过程中,通过相量测量单元采集实时电压数值,差分得到电压微变化量,电压下坠过程中通过分析相邻时刻电压微变化量是否发生突变,判断是否为故障清除后电压正常下坠;若判断出是故障清除后电压正常下坠,转至下一步骤,若不成立则进入下一采样时刻,重复步骤二;
步骤三,电压轨迹复合积分:在步骤二中判断故障清除后电压为正常下坠之后,以故障清除后电压下坠的实时数值为积分基数,加权考虑电压下坠速率进行复合积分得到复合积分基础值;若故障清除后,电压轨迹下坠的初值低于非正常电压基准,则对复合积分面积进行线性补偿得到复合积分最终值,否则,复合积分基础值即作为复合积分最终值;
步骤四,系统暂态稳定性判别:若电压下坠过程结束前,复合积分最终值大于暂态失稳门槛值,则判定系统暂态失稳,启动相应稳控装置;若电压下坠过程结束时,积分数值仍小于暂态失稳门槛值,则判定系统暂态稳定,转至步骤五;
步骤五,局部电压过低识别:步骤四中判定系统暂态稳定后,监测电压微变化量的数值,若在一定时间内小于设定值,则认为节点电压将逐步稳定在此时刻数值;若该时刻电压值低于非正常电压基准,则判定为局部电压过低,启动安全预警操作。
所述步骤一中遴选电压观测节点方法为:
根据电网实际运行方式及潮流图,分析系统主要的交流与直流传输线路;针对上述线路设置不同类型、不同地域的故障进行N-1连锁故障时域仿真,观测各变电站电压扰动响应曲线,在电网主要交流与直流传输断面附近选取真实反映系统暂态稳定特性的重要变电站,将该变电站的电压作为观测信息。
所述步骤二中故障清除后电压下坠识别方法为:
通过相量测量单元采集得到步骤一中电压观测节点的实时电压序列:{V(t),V(t+T),V(t+2T)。。。},t为采样起始时刻,T为采样周期,对连续时刻电压数值进行差分运算,得到电压微变化量:
ΔV(t+mT)=V(t+mT)-V(t+(m-1)T)
式中:m为大于1的整数;
并判断下式是否成立,如成立,则判断此时刻为故障清除后电压正常下坠,转至步骤三进行电压轨迹复合积分,若不成立则进入下一采样时刻,重复步骤二:
式中:k为线性约束因子,VU为非正常电压基准。
所述线性约束因子k取值范围优选为5-10,非正常电压基准VU优选为0.8p.u.。
所述步骤三中电压轨迹复合积分方法为:
步骤二中判别故障清除后电压正常下坠之后,对实时电压数值与电压下降速率进行复合积分运算:
式中:t0为积分启始时间,即是步骤二中判别正常下坠的时刻,tend为积分终止时间,A1为复合积分基础值;
若故障清除后的电压下坠初值低于非正常电压基准VU,则对基础复合积分值进行线性补偿:
式中:VD为电压下坠下限,V0为故障清除后电压下坠初值,A2为复合积分最终值;
若故障清除后的电压下坠初值V0高于非正常电压基准VU,则A2=A1。
所述电压下坠下限VD优选为0.2p.u.。
所述步骤四中系统暂态稳定性判别方法为:
在步骤三电压轨迹复合积分过程中,若满足连续三个采样时刻的电压微变化值:ΔV(t+(n+1)T),ΔV(t+(n+2)T),ΔV(t+(n+3)T)均大于零,则判定t+nT时刻的电压到达最低点,完成下坠过程,电压轨迹复合积分结束,其中n为大于1的整数;
在复合积分结束前判断下式是否成立,若成立,则判定系统暂态失稳:
A2>Aset
反之,则系统暂态稳定;式中,Aset为给定的暂态失稳门槛值。所述暂态失稳门槛值是电压稳定轨迹中的最大积分值。
所述步骤五中局部电压过低识别方法为:
步骤四系统暂态稳定性判别结束后,若从时刻t+aT开始,在一定时间域度t′内,t′∈[t+aT,t+(a+1)T,。。。t+(a+b)T],其中a,b为大于1的整数,若满足下式,则判定该观测节点附近存在局部电压过低问题,进行安全预警操作,
式中,ε为电压回稳约束项。
所述电压回稳约束项ε优选为0.05p.u.,连续采样个数b优选为100。
上述的p.u.的具体含义为:p.u.=实际电压值/基准电压值,无量纲,其中基准值可以选电网标称电压,也可以选电网平均电压。例如:当以标称电压作基准值时,如果p.u.=1,就表示实际电压值等于电网标称电压,即实际值=基准值×p.u.值。
与现有技术相比,本发明技术方案的有益效果是:本发明基于电网中现有的相量测量单元(PMU)量测配置以及广域量测系统(WAMS),着力解决二道防线控制范围以外的大扰动造成的系统暂态失稳。采用电压复合积分技术识别系统暂态失稳,尽早触发切机、切负荷等稳定控制措施,避免系统解列或大面积停电。具体具有如下优点:
1.在量测的选取上,本发明相对于功角类的判别方法具有一定的优势,可根据实际故障情况设置有针对性的布点方案。
2.本发明避免对受扰轨迹进行拟合,无需计算系统惯量中心,算法简单,对数据具有较好的容错性,在一定程度上减少了计算量和对数据精度的依赖性。
3.本发明用于暂态失稳预警,所需要量测少,计算量小,对于暂态失稳的识别率高,提高了判别的可靠性,并且相对于现有技术缩短了判别时间。
4.本发明还可用于局部电压较低的识别。
附图说明
图1是随故障时长变化的电压受扰轨迹示意图;
图2是本发明方法的总体流程框图;
图3是直流传输故障测点电压受扰轨迹示意图;
图4是交流线路故障测点电压受扰轨迹示意图;
图5是本发明方法的积分运算示意图;
图6是中国南方电网直流传输故障测点电压响应曲线(稳定算例);
图7是中国南方电网直流传输故障测点电压响应曲线(失稳算例);
图8是中国南方电网交流线路故障测点电压响应曲线(稳定算例);
图9是中国南方电网交流线路故障测点电压响应曲线(失稳算例)。
具体实施方式
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
图1给出了随故障时长变化的电压受扰轨迹示意图,其中11为稳定轨迹,轨迹12为临界轨迹,轨迹13为失稳轨迹。对于同一故障,电压轨迹的下降幅度以及下降的速度随着扰动强度的增大逐渐增大,系统暂态特性从稳定逐渐变为失稳。研究表明,通过观测一定时间窗内的电压数据,建立特定的积分办法,分析轨迹由稳定情况到失稳情况下的跃变特性,即可得到区分稳定轨迹和失稳轨迹的积分值,以此为门槛实现电力系统暂态失稳预警,称为暂态失稳门槛值。
下面结合图2的流程框图具体介绍本发明的基于扰动后电压轨迹的电力系统暂态失稳实时判别方法的步骤:
步骤一,遴选电压观测节点
根据电网实际运行方式及潮流图,分析系统主要的交流与直流传输线路;针对上述线路设置不同类型、不同地域的故障进行N-1连锁故障时域仿真,观测各变电站电压扰动响应曲线,在电网主要交流与直流传输断面附近选取真实反映系统暂态稳定特性的重要变电站,将该变电站的电压作为观测信息。
其中,主要的交流与直流传输线路的选取方式为:建立网络拓扑结构及等值模型,将电网作一定的简化处理并进行系统仿真,调节电网的运行方式为极限方式。在极限运行方式下得到潮流图,若线路的传输功率大于一定数值,则定义其为电网主要的交流与直流传输线路。
步骤二,故障清除后电压下坠识别
图3和图4分别给出了系统在交流线路故障及直流传输故障下,相关测点电压轨迹示意图。当交流线路发生三相短路接地故障,电压轨迹将瞬间跌落至低点。进行故障切除后,电压瞬间又会发生瞬间的回升现象,而直流传输故障不存在这一现象。由于短路接地造成的电压跌落只与故障点与观测点的距离相关,若计及故障过程中的电压积分将造成一定的误判。积分算法首先需要正确识别故障清除后的电压下坠。
以步骤一中得到的节点作为实际电压观测点,在电网运行过程中,通过相量测量单元采集得到实时电压序列:{V(t),V(t+T),V(t+2T)。。。},t为采样起始时刻,T为采样周期,对连续时刻电压数值进行差分运算,得到电压微变化量:
ΔV(t+mT)=V(t+mT)-V(t+(m-1)T)
式中:m为大于1的整数;
通过分析电压下坠过程中相邻时刻电压微变化量是否发生突变,判断是否为故障清除后电压正常下坠,即判断下式是否成立,如果成立,则判断为故障清除后电压正常下坠,转至步骤三。否则,即为非正常下坠,不进行积分运算,进入下一采样时刻,重复步骤二:
式中:k为线性约束因子,VU为非正常电压基准。
所述线性约束因子k取值范围为5-10,非正常电压基准VU为0.8p.u.。
步骤三,电压轨迹复合积分
通过大量离线仿真发现,同一故障下的扰动后电压轨迹的下降速率会随故障持续时间的增加而增大,且失稳轨迹的下降深度大致相同。因此,本发明采用以电压下降速率作为积分系数的复合积分方法。
在步骤二中判断故障清除后电压为正常下坠之后,以故障清除后电压下坠的实时数值为积分基数,加权考虑电压下坠速率进行复合积分运算:
式中:t0为积分启始时间,即是步骤二中判别正常下坠的时刻,tend为积分终止时间,A1为复合积分基础值;
若故障清除后的电压下坠初值低于非正常电压基准VU,则对基础复合积分值进行线性补偿:
式中:VD为电压下坠下限,V0为故障清除后电压下坠初值,A2为复合积分最终值;
若故障清除后的电压下坠初值V0高于非正常电压基准VU,则A2=A1。
所述电压下坠下限VD为0.2p.u.。
以稳定算例为例,电压积分过程如图5所示。
步骤四,系统暂态稳定性判别
为达到快速识别的目的,本发明采用故障清除后电压第一次摇摆至最低点作为复合积分结束时间,力求在电压第一次下坠过程中识别系统的暂态稳定性。
在步骤三电压轨迹复合积分过程中,若满足连续三个采样时刻的电压微变化值:ΔV(t+(n+1)T),ΔV(t+(n+2)T),ΔV(t+(n+3)T)均大于零,则判定t+nT时刻的电压到达最低点,完成下坠过程,电压轨迹复合积分结束,其中n为大于1的整数;
在复合积分结束前,判断下式是否成立,若成立,则判定系统暂态失稳,启动相应稳控装置:
A2>Aset
式中,Aset为给定的暂态失稳门槛值。
若电压下坠过程结束时,积分数值仍小于暂态失稳门槛值,则判定系统暂态稳定,转至步骤五;
步骤五,局部电压过低识别
除电力系统暂态失稳外,电网还存在着局部电压过低等安全隐患。本发明利用观测节点的电压逐步结束动态过程时的电压数值,分析是否存在局部电压过低的问题,达到地区电压安全预警的目的。
步骤四中判定系统暂态稳定后,监测电压微变化量的数值,若在一定时间内小于设定值,则认为节点电压将逐步稳定在此时刻数值;若该时刻电压值低于非正常电压基准,则判定为局部电压过低,启动安全预警操作。具体为:
步骤四系统暂态稳定性判别结束后,若从时刻t+aT开始,在一定时间域度t′内,t′∈[t+aT,t+(a+1)T,。。。t+(a+b)T],其中a,b为大于1的整数,若满足下式,则判定该观测节点附近存在局部电压过低问题,进行安全预警操作,
式中,ε为电压回稳约束项。
所述电压回稳约束项ε为0.05p.u.,连续采样个数b为100。
结果验证:为了测试本发明所提方法的有效性,应用本发明方法对中国南方电网系统不同故障情况进行了仿真验证。
建立中国南方电网的系统等值模型,利用BPA软件进行仿真分析,设置工作方式为2010年丰大极限运行方式。根据中国南方电网实际存在的“西电东送”运行情况,针对广西-广东主要的八条交流线路以及四条直流线路,设置传输断面附近重要变电站及换流站为电压观测点,如表1所示。
表1 中国南方电网电压观测点
名称 | 所属地区 | 电压等级 |
桂林 | 广西 | 500 |
贤令山 | 广西 | 500 |
柳东 | 广西 | 500 |
贺州 | 广西 | 500 |
罗洞 | 广东 | 500 |
沙塘 | 广西 | 500 |
来宾 | 广西 | 500 |
梧州 | 广西 | 500 |
百色 | 广西 | 500 |
南宁 | 广西 | 500 |
玉林 | 广西 | 500 |
茂名 | 广东 | 500 |
蝶岭 | 广东 | 500 |
在上述测点的基础上,对系统中一组不同类型、不同地域的故障进行大量离线仿真,得到扰动后电压轨迹,特别是临界稳定和临界失稳情况下的电压轨迹。采用本发明方法,通过大量仿真计算电压轨迹复合积分数值,积分时域为:故障清除后电压正常下坠到第一次到达最低点的过程。以稳定轨迹中的最大积分值作为暂态失稳判别的门槛值,得到可以区分稳定轨迹和失稳轨迹的复合积分值为2,即为暂态失稳门槛值Aset。
通过系统潮流图锁定主要的交直流传输线路,设置下列N-1连锁故障。
算例1:云广直流单极闭锁+天广直流单极闭锁。(稳定算例)
故障描述:0周波,云广直流发生单极闭锁,5周波,滤波器动作;0周波,天广直流发生单极闭锁,5周波,滤波器动作。其中,贤令山为受扰最严重测点。受扰最严重测点电压响应曲线如图6所示。
具体判别过程:
将实时电压数据带入下式计算,
得到直至1.04s这个时刻,上式才成立,即判断此时刻为故障清除后电压正常下坠过程,然后按照下式进行电压轨迹积分运算:
其中t0为1.04s,在之后的开始进行复合积分运算。因为此次下坠过程的电压下坠初值为0.98p.u.>0.8p.u.,即大于非正常电压基准VU的值,所以不必进行线性补偿,此复合积分初始值A1即为复合积分最终值,即A2=A1,将A2的值与暂态失稳门槛值Aset=2进行实时比较。该复合积分值直到此次下坠过程的最低点时,即1.75s时,其复合积分最终值(即此次下坠过程的最大复合积分值)A2=A1=0.0975仍小于暂态失稳门槛值Aset=2,即A2<Aset,可判断出此时系统暂态稳定。此判断结果与系统真实情况相符合,判断准确。
算例2:云广直流单极闭锁+新安直流单极闭锁。(失稳算例)
故障描述:0周波,云广直流发生单极闭锁,5周波,滤波器动作;0周波,新安直流发生单极闭锁,5周波,滤波器动作。其中,贺州为受扰最严重测点。受扰最严重测点电压响应曲线如图7所示。
具体判别过程同算例1类似:
利用本发明的方法,识别电压为正常下坠后,即以此时刻为积分起始时刻,对各测点电压进行复合积分操作。贺州测点电压积分开始时间t0为0.78s,当积分运算至1.44s时(还未达到此次下坠过程的最低点),复合积分数值就已超过了暂态失稳门槛值Aset=2,在电压第一次下坠结束前已正确识别系统暂态失稳。
算例3:柳东-贺州三永+来宾-梧州跳双回。(稳定算例)
故障描述:0周波,柳东-贺州发生三相永久短路,5周波柳东-贺州跳双回以及玉林-茂名双回线。其中,梧州为受扰最严重测点。受扰最严重测点电压响应曲线如图8所示。
具体判别过程同算例1类似:
根据本发明的方法,识别出电压在0s-0.1s时间段的短路过程中非正常下坠,避免了0s-0.1s的复合积分过程,有效防止误判的发生。而判断出在0.61s-1.98s的电压下坠过程为故障清除后电压正常下坠,按照本发明的方法进行复合积分运算。计算得到,直到电压下坠最低点时,即1.98s时,其复合积分最终值(即此次下坠过程的最大复合积分值)为0.688,小于暂态失稳门槛值Aset=2,判定系统稳定,与实际情况相符合。本算例中的电压下坠初值高于非正常电压基准VU的值,所以不必进行线性补偿,积分电压基准值即作为积分电压最终值。
算例4:来宾-梧州三永+桂林-贤令山跳双回。(失稳算例)
故障描述:0周波,来宾到梧州发生三相永久短路,5周波来宾-梧州跳双回线以及玉林-茂名双回线。其中,茂名为受扰最严重测点。受扰最严重测点电压响应曲线如图9所示。
具体判别过程同算例1类似:
根据本发明的方法,识别出电压在0s-0.1s时间段为非正常下坠,避免了0s-0.1s的复合积分过程,有效防止误判的发生。而判断出从0.63s开始的电压下坠过程为故障清除后电压正常下坠,按照本发明的方法,进行复合积分运算。电压积分运算开始时间t0为0.63s,当积分运算至1.19s时(还未达到此次下坠过程的最低点),复合积分数值就已超过了暂态失稳门槛值Aset=2,在电压第一次下坠结束前已正确识别系统暂态失稳。
利用本发明方法分析上述算例系统的暂态稳定性,并与直接积分法及功角门槛值法作对比,判定结果如表2所示。表2中的时间表示判断出系统暂稳失态所需的时间。
表2
故障编号 | 本发明 | 直接积分法 | 功角门槛值法 |
算例1 | 稳定 | 稳定 | 稳定 |
算例2 | 1.44秒 | 1.54秒 | 1.49秒 |
算例3 | 稳定 | 稳定 | 稳定 |
算例4 | 1.19秒 | 1.32秒 | 1.65秒 |
从上述判别结果中可以发现,本发明方法识别准确且较传统方法更为快速。
本发明可用于各级电网的调度控制中心,基于广域量测系统,实现对电力系统暂态稳定的实时预测和预警。
上述所述是本发明最佳实施例,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于扰动后电压轨迹的电力系统暂态失稳实时判别方法,其特征在于,包括以下步骤:
步骤一,遴选电压观测节点:按照电网的实际运行方式,选取真实反映系统暂态稳定特性的重要变电站,将其电压作为观测信息;
步骤二,故障清除后电压下坠识别:以步骤一中得到的节点作为实际电压观测点,在电网运行过程中,通过相量测量单元采集实时电压数值,差分得到电压微变化量,电压下坠过程中通过分析相邻时刻电压微变化量是否发生突变,判断是否为故障清除后电压正常下坠;若判断出是故障清除后电压正常下坠,则转至下一步骤,若不成立则进入下一采样时刻,重复步骤二;
步骤三,电压轨迹复合积分:在步骤二中判断故障清除后电压为正常下坠之后,以故障清除后电压下坠的实时数值为积分基数,加权考虑电压下坠速率进行复合积分得到复合积分基础值;若故障清除后,电压轨迹下坠的初值低于非正常电压基准,则对复合积分面积进行线性补偿得到复合积分最终值,否则,复合积分基础值即作为复合积分最终值;
步骤四,系统暂态稳定性判别:若电压下坠过程结束前,复合积分最终值大于暂态失稳门槛值,则判定系统暂态失稳,启动相应稳控装置;若电压下坠过程结束时,复合积分最终值仍小于暂态失稳门槛值,则判定系统暂态稳定,转至步骤五;
步骤五,局部电压过低识别:步骤四中判定系统暂态稳定后,监测电压微变化量的数值,若在一定时间内小于设定值,则认为节点电压将逐步稳定在此时刻数值;若该时刻电压值低于非正常电压基准,则判定为局部电压过低,启动安全预警操作。
2.根据权利要求1所述的基于扰动后电压轨迹的电力系统暂态失稳实时判别方法,其特征在于,所述步骤一中遴选电压观测节点方法为:
根据电网实际运行方式及潮流图,分析系统主要的交流与直流传输线路;针对上述线路设置不同类型、不同地域的故障进行N-1连锁故障时域仿真,观测各变电站电压扰动响应曲线,在电网主要交流与直流传输断面附近选取真实反映系统暂态稳定特性的重要变电站,将该变电站的电压作为观测信息。
3.根据权利要求1所述的基于扰动后电压轨迹的电力系统暂态失稳实时判别方法,其特征在于,所述步骤二中故障清除后电压下坠识别方法为:
通过相量测量单元采集得到步骤一中电压观测节点的实时电压序列:{V(t),V(t+T),V(t+2T)…},t为采样起始时刻,T为采样周期,对连续时刻电压数值进行差分运算,得到电压微变化量:
ΔV(t+mT)=V(t+mT)-V(t+(m-1)T)
式中:m为大于1的整数;
并判断下式是否成立,如成立,则判断此时刻为故障清除后电压正常下坠,转至步骤三进行电压轨迹复合积分,若不成立则进入下一采样时刻,重复步骤二:
式中:k为线性约束因子,VU为非正常电压基准。
4.根据权利要求3所述的基于扰动后电压轨迹的电力系统暂态失稳实时判别方法,其特征在于,所述线性约束因子k取值范围为5-10,非正常电压基准VU为0.8p.u.。
5.根据权利要求3所述的基于扰动后电压轨迹的电力系统暂态失稳实时判别方法,其特征在于,所述步骤三中电压轨迹复合积分方法为:
步骤二中判别故障清除后电压正常下坠之后,对实时电压数值与电压下坠速率进行复合积分运算:
式中:t0为积分启始时间,即是步骤二中判别正常下坠的时刻,tend为积分终止时间,A1为复合积分基础值;
若故障清除后的电压下坠初值低于非正常电压基准VU,则对复合积分基础值进行线性补偿:
式中:VD为电压下坠下限,V0为故障清除后电压下坠初值,A2为复合积分最终值;
若故障清除后的电压下坠初值V0高于非正常电压基准VU,则A2=A1。
6.根据权利要求5所述的基于扰动后电压轨迹的电力系统暂态失稳实时判别方法,其特征在于,所述电压下坠下限VD为0.2p.u.。
7.根据权利要求5所述的基于扰动后电压轨迹的电力系统暂态失稳实时判别方法,其特征在于,所述步骤四中系统暂态稳定性判别方法为:
在步骤三电压轨迹复合积分过程中,若满足连续三个采样时刻的电压微变化值:ΔV(t+(n+1)T),ΔV(t+(n+2)T),ΔV(t+(n+3)T)均大于零,则判定t+nT时刻的电压到达最低点,完成下坠过程,电压轨迹复合积分结束,其中n为大于1的整数;
在复合积分结束前判断下式是否成立,若成立,则判定系统暂态失稳:
A2>Aset
反之,则系统暂态稳定;式中,Aset为暂态失稳门槛值,为电压稳定轨迹中的最大积分值。
8.根据权利要求7所述的基于扰动后电压轨迹的电力系统暂态失稳实时判别方法,其特征在于,所述步骤五中局部电压过低识别方法为:
步骤四系统暂态稳定性判别结束后,若从时刻t+aT开始,在一定时间域度t′内,t'∈[t+aT,t+(a+1)T,…t+(a+b)T],其中a,b为大于1的整数,若满足下式,则判定该观测节点附近存在局部电压过低问题,进行安全预警操作,
式中,ε为电压回稳约束项。
9.根据权利要求8所述的基于扰动后电压轨迹的电力系统暂态失稳实时判别方法,其特征在于,所述电压回稳约束项ε为0.05p.u.,连续采样个数b为100。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110417168 CN102522759B (zh) | 2011-12-13 | 2011-12-13 | 一种基于扰动后电压轨迹的电力系统暂态失稳实时判别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110417168 CN102522759B (zh) | 2011-12-13 | 2011-12-13 | 一种基于扰动后电压轨迹的电力系统暂态失稳实时判别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102522759A CN102522759A (zh) | 2012-06-27 |
CN102522759B true CN102522759B (zh) | 2013-11-06 |
Family
ID=46293582
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110417168 Active CN102522759B (zh) | 2011-12-13 | 2011-12-13 | 一种基于扰动后电压轨迹的电力系统暂态失稳实时判别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102522759B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104750994B (zh) * | 2015-04-03 | 2017-05-17 | 国家电网公司 | 基于外区电力比例的暂态失稳主导性判别方法 |
CN105811437A (zh) * | 2016-04-19 | 2016-07-27 | 贵州电网有限责任公司 | 一种基于广域响应的暂态稳定在线预判系统实现方法 |
CN106769646B (zh) * | 2016-11-11 | 2020-03-03 | 青岛科技大学 | 一种基于反演思想的精馏过程双变量扰动原因诊断方法 |
CN107895936B (zh) * | 2017-09-28 | 2019-05-03 | 国网辽宁省电力有限公司 | 一种数字化虚拟失步继电器 |
CN107919666B (zh) * | 2017-12-07 | 2020-11-17 | 贵州电网有限责任公司电力科学研究院 | 一种基于广域响应的地区电网暂态稳定在线综合预判方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1110838A (zh) * | 1995-02-25 | 1995-10-25 | 电力部电力自动化研究院 | 电力系统暂态稳定的在线监视和预防控制系统 |
US5566085A (en) * | 1992-09-24 | 1996-10-15 | Hydro-Quebec | Stability transfer limit calculation in a power transmission network |
CN101359026A (zh) * | 2008-09-17 | 2009-02-04 | 中国电力科学研究院 | 一种全时域电压稳定性综合评估方法 |
CN101625389A (zh) * | 2009-08-06 | 2010-01-13 | 中国电力科学研究院 | 一种考虑负荷变化特性的电力系统稳定性分析方法 |
-
2011
- 2011-12-13 CN CN 201110417168 patent/CN102522759B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5566085A (en) * | 1992-09-24 | 1996-10-15 | Hydro-Quebec | Stability transfer limit calculation in a power transmission network |
CN1110838A (zh) * | 1995-02-25 | 1995-10-25 | 电力部电力自动化研究院 | 电力系统暂态稳定的在线监视和预防控制系统 |
CN101359026A (zh) * | 2008-09-17 | 2009-02-04 | 中国电力科学研究院 | 一种全时域电压稳定性综合评估方法 |
CN101625389A (zh) * | 2009-08-06 | 2010-01-13 | 中国电力科学研究院 | 一种考虑负荷变化特性的电力系统稳定性分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102522759A (zh) | 2012-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103279639B (zh) | 基于响应的受端电网电压稳定全过程态势评估及防控方法 | |
Sodhi et al. | A simple scheme for wide area detection of impending voltage instability | |
CN102522759B (zh) | 一种基于扰动后电压轨迹的电力系统暂态失稳实时判别方法 | |
CN103065016A (zh) | 确定引起直流换相失败的交流系统故障范围的方法及系统 | |
CN103632043B (zh) | 一种基于实测响应信息的电力系统主导失稳模式识别方法 | |
Wu et al. | Review on transient stability prediction methods based on real time wide-area phasor measurements | |
CN103592575A (zh) | 一种基于多传感器系统的自适应加权数据融合故障测距方法 | |
CN103995948B (zh) | 一种基于多项式模型的振荡中心电压预测方法 | |
CN103050933B (zh) | 基于单端电流突变量的大型电池储能电站接口保护方法 | |
CN104252686A (zh) | 一种电网安全综合指标确定方法 | |
CN104505855B (zh) | 一种风电场并网特性在线评价系统 | |
CN107919666B (zh) | 一种基于广域响应的地区电网暂态稳定在线综合预判方法 | |
CN109428327B (zh) | 基于响应的电网关键支路与主导稳定模式识别方法及系统 | |
CN109217263A (zh) | 一种实现带并联电抗器输电线路的自适应重合闸方法 | |
CN105447774A (zh) | 一种电网安全事故风险等级在线评估与预警方法及装置 | |
CN103592528A (zh) | 一种基于动态轨迹灵敏度的光伏逆变器模型参数辨识方法 | |
CN108845225B (zh) | 电力电容器、电抗器二次电流回路接线正确性分析方法 | |
CN102510072B (zh) | 一种电网系统暂态失稳判别方法 | |
CN102723710B (zh) | 一种用于电力系统暂态稳定分析的临界机群识别方法 | |
CN108614189A (zh) | 双馈风电场送出线三相短路故障方向判断方法及方向元件 | |
CN101825677B (zh) | 基于环流量电流相角差的同塔四回输电线路故障选相方法 | |
CN106599425A (zh) | 电力系统暂态电压稳定检测方法和装置 | |
CN103647284B (zh) | 一种解决单时间断面问题的电压稳定预测方法 | |
CN104502800A (zh) | 一种电力系统故障因子特征提取方法 | |
CN105896525B (zh) | 基于wams的电网大扰动识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |