CN104135017B - 一种基于振荡能量的失步解列判别方法 - Google Patents
一种基于振荡能量的失步解列判别方法 Download PDFInfo
- Publication number
- CN104135017B CN104135017B CN201410353928.8A CN201410353928A CN104135017B CN 104135017 B CN104135017 B CN 104135017B CN 201410353928 A CN201410353928 A CN 201410353928A CN 104135017 B CN104135017 B CN 104135017B
- Authority
- CN
- China
- Prior art keywords
- interconnection
- disconnection device
- installation place
- device installation
- oscillation
- 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.)
- Expired - Fee Related
Links
- 230000010355 oscillation Effects 0.000 title claims abstract description 114
- 238000000034 method Methods 0.000 title claims abstract description 14
- 238000009434 installation Methods 0.000 claims abstract description 68
- 238000005070 sampling Methods 0.000 claims abstract description 36
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000004146 energy storage Methods 0.000 description 8
- 230000010354 integration Effects 0.000 description 6
- 230000001052 transient effect Effects 0.000 description 5
- 230000003111 delayed effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000001360 synchronised effect Effects 0.000 description 4
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了电力***规划与运行控制技术领域中的一种基于振荡能量的失步解列判别方法。包括:判断***是否发生故障,如果***发生故障,则根据联络线上解列装置安装处的连续设定次采样数据,计算每次采样后联络线上解列装置安装处的振荡能量;再根据每次采样后联络线上解列装置安装处的振荡能量判断***联络线是否发生失步振荡,如果***联络线发生失步振荡,则控制联络线解列装置进行解列。本发明能够根据振荡能量的变化对联络线上的失步振荡做出正确判断,进而为失步解列提供准确判据,避免了由于失步振荡识别错误导致的大规模停电事故的发生。
Description
技术领域
本发明属于电力***规划与运行控制技术领域,尤其涉及一种基于振荡能量的失步解列判别方法。
背景技术
随着互联电力***规模的不断扩大,失步振荡发生的频率不断增加。作为电力***第三道防线的重要组成部分(前两道防线为继电保护和安全自动装置),失步解列判据对于***失步振荡的识别及防止事故蔓延导致大规模停电至关重要。在电力***遭受严重干扰,电网的整体性无法保持时,解列装置对大电网进行解列能够避免电网全面崩溃,保证重要负荷的持续供电。
目前,大多数复杂***建模的方法都是基于网络互联化思想,把复杂***看成是由许多简单的带端口的子***按照能量守恒的原则逐级递归互联形成。电力***是复杂的***,而振荡问题是子***间的相互动态影响问题。如果能如上所述,把电力***这个复杂的大***分成许多带端口的子***,那么振荡问题的研究就能转化成为端口间动态影响的研究。
在这个思路上,面对的两个关键问题是:1)能否严谨直观地描述***的端口互联结构;2)能否恰当的提出能量的概念,满足电力***的应用要求。
由此,本发明引入了Dirac结构和端口Hamilton***。Dirac结构是数学领域上的一种本质上表征能量守恒性质的几何结构;端口Hamilton***应用于机械领域,给出的Hamilton函数可以表征***的储能。在Dirac结构和端口Hamilton***联立的基础上,提出能量结构的概念,可以证明能量结构上的能量满足李雅普诺夫能量函数,可以应用于电力***的动态分析,端口上流动的能量可以用来表征与对象子***相连的子***对对象子***的动态影响。
发明内容
本发明的目的在于,提供一种基于振荡能量的失步解列判别方法,利用联络线上的振荡能量变化判断***是否发生失步振荡,并在***发生失步振荡时启动解列装置进行解列操作,从而确保电网的安全运行。
为了实现上述目的,本发明提出的技术方案是,一种基于振荡能量的失步解列判别方法,其特征是所述方法包括:
步骤1:判断***是否发生故障,如果***发生故障,则执行步骤2;否则,继续执行步骤1;
步骤2:根据联络线上解列装置安装处的连续设定次采样数据,计算每次采样后联络线上解列装置安装处的振荡能量;
步骤3:根据每次采样后联络线上解列装置安装处的振荡能量判断***联络线是否发生失步振荡,如果***联络线发生失步振荡,则执行步骤4;否则,返回步骤1;
步骤4:控制联络线解列装置进行解列。
所述计算每次采样后联络线上解列装置安装处的振荡能量采用公式:
其中,OEF(t)为时刻t联络线上解列装置安装处振荡能量;
P为采样时联络线上解列装置安装处的有功功率;
Q为采样时联络线上解列装置安装处的无功功率;
U为采样时联络线上解列装置安装处的电压幅值;
为采样时联络线上解列装置安装处的电压幅值的微分;
θ为采样时联络线上解列装置安装处的电压相角;
为采样时联络线上解列装置安装处的电压相角的微分。
所述根据每次采样后联络线上解列装置安装处的振荡能量判断***联络线是否发生失步振荡具体为:
如果每次采样后联络线上解列装置安装处的振荡能量均为正数,则***联络线未发生失步振荡;
如果每次采样后联络线上解列装置安装处的振荡能量均为负数,则***联络线未发生失步振荡;
如果每次采样后联络线上解列装置安装处的振荡能量由正值变为负值,之后又变为正值并且稳定为正值,则***联络线未发生失步振荡;
如果每次采样后联络线上解列装置安装处的振荡能量由负值变为正值,之后又变为负值并且稳定为负值,则***联络线未发生失步振荡;
如果每次采样后联络线上解列装置安装处的振荡能量呈现两周以上的正值和负值交替变化,则***联络线发生失步振荡。
本发明能够根据振荡能量的变化对联络线上的失步振荡做出正确判断,进而为失步解列提供准确判据,避免了由于失步振荡识别错误导致的大规模停电事故的发生。
附图说明
图1是等值两机***结构图;
图2是等值两机***间发生失步振荡的一个振荡周期内整条联络线上的振荡能量分布图;
图3是基于振荡能量的失步解列判别方法流程图。
具体实施方式
下面结合附图,对优选实施例作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
本发明提供的方法以基于Dirac结构和端口Hamilton***的电力***能量结构为基础,该能量结构上的能量满足李雅普诺夫能量函数的条件,子***间的交互能量可以用来表征包括振荡问题在内的***间的动态影响。
基于Dirac结构和端口Hamilton***的电力***能量结构如下:
公式(1)中,V为基于Dirac结构和端口Hamilton***的电力***能量结构上的能量,该能量满足李雅普诺夫能量函数。n为电力***中母线数量,nG为电力***中发电机母线数量。公式(1)前三项为发电机节点的能量,且有:
公式(2)中,Mi为第i台发电机的惯性时间常数,ω0为发电机的同步转速,ωi为第i台发电机的转速,Eqi为第i台发电机的同步电动势q轴分量,E′qi为第i台发电机的暂态电动势q轴分量,Edi为第i台发电机的同步电动势d轴分量,E′di为第i台发电机的暂态电动势d轴分量,xdi为第i台发电机的暂态电抗d轴分量,x′di为第i台发电机的次暂态电抗d轴分量,xqi为第i台发电机的暂态电抗q轴分量,x′qi为第i台发电机的次暂态电抗q轴分量,Idi为第i个台发电机的注入电流d轴分量,Iqi为第i个台发电机的注入电流q轴分量,Pmi为第i台发电机的机械功率,Pei为第i台发电机的电磁功率,δi为第i台发电机的功角,δi0为第i台发电机的功角积分初始值,Uqi为第i台发电机的电压幅值q轴分量,Uqi0为第i台发电机的电压幅值q轴分量的积分初始值,Udi为第i台发电机的电压幅值d轴分量,Udi0为第i台发电机的电压幅值d轴分量的积分初始值,Eqi0为第i台发电机的同步电动势q轴分量积分初始值,E′qi0为第i台发电机的暂态电动势q轴分量积分初始值,Efi为第i台发电机的励磁电压。
公式(1)中,第四项为内部网络中的储能能量,且有:
公式(3)中,Bii为第i个节点的自电纳,Bij为支路ij的互电纳,θij为母线ij的电压相角差,Ui为第i个节点的电压幅值,Uj为第j个节点的电压幅值。
公式(1)中,第五项为与外部***交换的能量,且有:
公式(4)中,Pi为第i个节点与外部网络交换的有功功率,Qi为第i个节点与外部网络交换的无功功率,θi为第i个节点的电压相角,Ui为第i个节点的电压幅值,θi0为母线i的电压相角积分初始值,Ui0为母线i的电压幅值积分初始值。
子***的储能是子***振荡剧烈程度的标尺,振荡越剧烈,则其储能数值越大;子***的储能受到端***换能量和内部供给能量以及本身耗散的影响。若将内部供给能量也视为能量供给元件的端口注入,则储能即为***所有端口供给能量之和减去其本身耗散。若所有端***换能量都为负,则对象***内所有子***的储能一定会减少到平衡点时的储能数值,意即振荡的平息。
将交换能量推至联络线,定义联络线上的振荡能量(OscillationEnergyFlow)为:
公式(5)中,P为联络线上解列装置安装处的有功功率,Q为联络线上解列装置安装处的无功功率,θ为联络线上解列装置安装处的电压相角,为联络线上解列装置安装处的电压相角的微分,U为联络线上解列装置安装处的电压幅值,为联络线上解列装置安装处的电压幅值的微分。
失步解列判别的工作原理为:以图1所示的等值两机***为例。联络线两侧等值电动势分别为:领前侧滞后侧为参考侧。等值功角差为δ=Δωt+δ0,令初始功角差δ0=0,则δ=Δωt,x为解列装置安装处与领前侧之间的电抗,x∑为两侧电势源之间的总电抗,设代入解列装置安装处测得的有功功率、无功功率、电压相角与电压幅值,可以得到此处的振荡能量如下:
公式(6)中,Δδ为功角变化量,k为等值两机***两侧电动势幅值的比值且E1为等值两机***的一侧(领前侧)电动势幅值,E2为等值两机***的另一侧(滞后侧)电动势幅值。zΣ为两电势源之间的总电抗,包括两个发电机的内阻抗与线路的阻抗。
由公式(6)可知,在一个失步振荡周期内,联络线两侧的等值功角差δ由0°变化到180°后又变化到360°。在δ∈(0°,180°)时,线路上的振荡能量方向由领前侧向滞后侧传输;δ∈(180°,360°)时,线路上的振荡能量方向由滞后侧向领前侧传输。图2进一步说明了这一点,图2中,x轴表示等值功角差的变化,y轴表示考察点在联络线上的变化(用ρ来表示),z轴表示振荡能量的幅值和方向。等值两机***的电压,线路的阻抗等参数均使用的标幺值,计算出的有功无功与振荡能量也符合标幺制。
基于上述原理分析,本发明提供的一种基于振荡能量的失步解列判别方法包括:
步骤1:判断***是否发生故障,如果***发生故障,则执行步骤2;否则,继续执行步骤1。
通常,***是否发生故障通过联络线上的电流变化率进行判别。当联络线上的电流变化率ΔI>ΔIset时,判定***发生故障。其中,ΔI为联络线上的电流变化率,ΔIset为设定阈值。
步骤2:利用联络线上解列装置安装处的相量测量单元(PMU,PhasorMeasurementUnit)采样的连续设定次数据,计算每次采样后联络线上解列装置安装处的振荡能量。
设定次数N的取值为其中T为经验值。由于失步振荡的频率远低于电力***的工频(50hz),因此失步振荡的周期比0.02秒长得多。即使按工频每个周波采两个点,也可以保证有足够的点。如果对这些点拟合,可以得到类似正弦的曲线。所以,经验值T的取值只要确保设定次数N能够涵盖两个以上的失步振荡周期即可。
每次通过PMU采样后,都可以得到解列装置安装处的有功功率、无功功率、电压幅值和电压相角,根据公式(5)即可计算出该采样时刻的解列装置安装处的振荡能量。
步骤3:根据每次采样后联络线上解列装置安装处的振荡能量判断***联络线是否发生失步振荡,判断方式如下:
(1)如果每次采样后联络线上解列装置安装处的振荡能量均为正数,则***联络线未发生失步振荡。
根据N次采样数据得到的结果,经过计算得到N个解列装置安装处的振荡能量值。如果N个解列装置安装处的振荡能量值都大于零,则***联络线未发生失步振荡。
(2)如果每次采样后联络线上解列装置安装处的振荡能量均为负数,则***联络线未发生失步振荡。
根据N次采样数据得到的结果,经过计算得到N个解列装置安装处的振荡能量值。如果N个解列装置安装处的振荡能量值都小于零,则***联络线未发生失步振荡。
(3)如果每次采样后联络线上解列装置安装处的振荡能量由正值变为负值,之后又变为正值并且稳定为正值,则***联络线未发生失步振荡。
根据N次采样数据得到的结果,经过计算得到N个解列装置安装处的振荡能量值。如果N个解列装置安装处的振荡能量值前若干个值大于零,之后若干个值小于零,再之后的值都大于零,则***联络线未发生失步振荡。
(4)如果每次采样后联络线上解列装置安装处的振荡能量由负值变为正值,之后又变为负值并且稳定为负值,则***联络线未发生失步振荡。
根据N次采样数据得到的结果,经过计算得到N个解列装置安装处的振荡能量值。如果N个解列装置安装处的振荡能量值前若干个值小于零,之后若干个值大于零,再之后的值都小于零,则***联络线未发生失步振荡。
(5)如果每次采样后联络线上解列装置安装处的振荡能量呈现两周以上的正值和负值交替变化,则***联络线发生失步振荡。
根据N次采样数据得到的结果,经过计算得到N个解列装置安装处的振荡能量值。
如果N个解列装置安装处的振荡能量值有若干个值大于零,之后的若干个值小于零,接下来若干个值又大于零,再接下来的若干个值又小于零,则***联络线发生失步振荡。
或者,如果N个解列装置安装处的振荡能量值有若干个值小于零,之后的若干个值大于零,接下来若干个值又小于零,再接下来的若干个值又大于零,则***联络线发生失步振荡。
依据上述(1)-(5)的判断方式,如果判断***联络线发生失步振荡,则执行步骤4。如果判断***联络线未发生失步振荡,则返回步骤1。
步骤4:控制联络线解列装置进行解列。
本发明利用联络线上的振荡能量变化判断***是否发生失步振荡,并在***发生失步振荡时启动解列装置进行解列操作,确保了电网的安全运行。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (2)
1.一种基于振荡能量的失步解列判别方法,其特征是所述方法包括:
步骤1:判断***是否发生故障,如果***发生故障,则执行步骤2;否则,继续执行步骤1;
步骤2:在联络线上解列装置安装处连续采集N次数据,计算每次采样后联络线上解列装置安装处的振荡能量,设定次数N的取值为其中T为经验值;
步骤3:根据每次采样后联络线上解列装置安装处的振荡能量判断***联络线是否发生失步振荡,如果***联络线发生失步振荡,则执行步骤4;否则,返回步骤1;
步骤4:控制联络线解列装置进行解列;
所述计算每次采样后联络线上解列装置安装处的振荡能量采用公式:
其中,OEF(t)为时刻t联络线上解列装置安装处振荡能量;
P为采样时联络线上解列装置安装处的有功功率;
Q为采样时联络线上解列装置安装处的无功功率;
U为采样时联络线上解列装置安装处的电压幅值;
为采样时联络线上解列装置安装处的电压幅值的微分;
θ为采样时联络线上解列装置安装处的电压相角;
为采样时联络线上解列装置安装处的电压相角的微分。
2.根据权利要求1所述的方法,其特征是所述根据每次采样后联络线上解列装置安装处的振荡能量判断***联络线是否发生失步振荡具体为:
如果每次采样后联络线上解列装置安装处的振荡能量均为正数,则***联络线未发生失步振荡;
如果每次采样后联络线上解列装置安装处的振荡能量均为负数,则***联络线未发生失步振荡;
如果每次采样后联络线上解列装置安装处的振荡能量由正值变为负值,之后又变为正值并且稳定为正值,则***联络线未发生失步振荡;
如果每次采样后联络线上解列装置安装处的振荡能量由负值变为正值,之后又变为负值并且稳定为负值,则***联络线未发生失步振荡;
如果每次采样后联络线上解列装置安装处的振荡能量呈现两周以上的正值和负值交替变化,则***联络线发生失步振荡;
所述呈现两周以上的正负值交替变化是指根据N次采样数据计算得到的N个解列装置安装处的振荡能量值中有若干个值大于零,之后的若干个值小于零,接下来若干个值又大于零,再接下来的若干个值又小于零;或者根据N次采样数据计算得到的N个解列装置安装处的振荡能量值中有若干个值小于零,之后的若干个值大于零,接下来若干个值又小于零,再接下来的若干个值又大于零。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410353928.8A CN104135017B (zh) | 2014-07-23 | 2014-07-23 | 一种基于振荡能量的失步解列判别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410353928.8A CN104135017B (zh) | 2014-07-23 | 2014-07-23 | 一种基于振荡能量的失步解列判别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104135017A CN104135017A (zh) | 2014-11-05 |
CN104135017B true CN104135017B (zh) | 2016-04-20 |
Family
ID=51807593
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410353928.8A Expired - Fee Related CN104135017B (zh) | 2014-07-23 | 2014-07-23 | 一种基于振荡能量的失步解列判别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104135017B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104569820B (zh) * | 2015-01-15 | 2017-07-28 | 华南理工大学 | 一种基于fft变换的多机群失步解列判断方法 |
CN106556776A (zh) * | 2016-11-10 | 2017-04-05 | 华北电力大学 | 一种基于振荡能量的失步振荡中心捕获方法 |
CN107102192B (zh) * | 2017-04-28 | 2019-07-23 | 许继集团有限公司 | 一种全相异步振荡判别方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1419133A (zh) * | 2002-09-30 | 2003-05-21 | 南京南瑞继保电气有限公司 | 基于ucosφ的电力***失步检测判别方法 |
CN101345409A (zh) * | 2008-08-19 | 2009-01-14 | 国网电力科学研究院 | 一种基于就地有功无功变化量判断失步解列的设计方法 |
CN102064524A (zh) * | 2010-11-10 | 2011-05-18 | 中电普瑞科技有限公司 | 一种电力***失步快速解列判断方法 |
CN103337831A (zh) * | 2013-06-17 | 2013-10-02 | 国家电网公司 | 一种具有自适应功能的失步解列方法 |
-
2014
- 2014-07-23 CN CN201410353928.8A patent/CN104135017B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1419133A (zh) * | 2002-09-30 | 2003-05-21 | 南京南瑞继保电气有限公司 | 基于ucosφ的电力***失步检测判别方法 |
CN101345409A (zh) * | 2008-08-19 | 2009-01-14 | 国网电力科学研究院 | 一种基于就地有功无功变化量判断失步解列的设计方法 |
CN102064524A (zh) * | 2010-11-10 | 2011-05-18 | 中电普瑞科技有限公司 | 一种电力***失步快速解列判断方法 |
CN103337831A (zh) * | 2013-06-17 | 2013-10-02 | 国家电网公司 | 一种具有自适应功能的失步解列方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104135017A (zh) | 2014-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sexauer et al. | Phasor measurement units for the distribution grid: Necessity and benefits | |
CN102636728B (zh) | 基于力矩分解法识别负阻尼低频振荡的方法 | |
CN104135017B (zh) | 一种基于振荡能量的失步解列判别方法 | |
WO2012163979A2 (en) | Assessment of power systems | |
CN103606922B (zh) | 一种基于典型故障集的电力***功角稳定性近似判定方法 | |
CN101794998A (zh) | 基于多机电力***单一发电机电磁功率简洁表达方式的暂态稳定在线分析方法 | |
Hu et al. | An adaptive three-bus power system equivalent for estimating voltage stability margin from synchronized phasor measurements | |
CN102611128B (zh) | 一种高压直流输电***直流功率调制方法 | |
CN104950196A (zh) | 一种识别电力***振荡中心的复合判据方法 | |
CN103217896A (zh) | 基于自由权矩阵方法的多facts抗时滞协调控制方法 | |
CN103217623A (zh) | 与过渡电阻和负荷电流无关的线路双端故障测距方法 | |
You et al. | Us eastern interconnection (ei) electromechanical wave propagation and the impact of high pv penetration on its speed | |
CN108075478A (zh) | 一种暂态稳定紧急控制方法 | |
Xu et al. | Real-time damping estimation on nonlinear electromechanical oscillation | |
Byrne et al. | Optimal locations for energy storage damping systems in the Western North American interconnect | |
Saunders et al. | Transient stability analysis using potential energy indices for determining critical generator sets | |
Wu et al. | Decentralized control for satellite formation using local relative measurements only | |
CN104316827A (zh) | 一种电力***振荡中心定位方法 | |
Azim et al. | Design of a general droop controller for islanded microgrids | |
CN104809512B (zh) | 一种电力***不同紧急控制措施性能指标的快速评估方法 | |
KR20090064796A (ko) | 지역측정 페이저 정보를 이용하여 과도 안정도를 감시하기위한 시스템 및 그 방법 | |
Abdelaziz et al. | Phasor measurement units for out of step detection | |
Zhou et al. | The transient stability analysis based on WAMS and online admittance parameter identification | |
Farantatos et al. | Real-time power system dynamic equivalencing to preserve system center of oscillations via PMU-based dynamic state estimator | |
Zhang et al. | Robust H∞ Load Frequency Control of future power grid with energy storage considering parametric uncertainty and time delay |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160420 Termination date: 20160723 |