CN104135017A - 一种基于振荡能量的失步解列判别方法 - Google Patents
一种基于振荡能量的失步解列判别方法 Download PDFInfo
- Publication number
- CN104135017A CN104135017A CN201410353928.8A CN201410353928A CN104135017A CN 104135017 A CN104135017 A CN 104135017A CN 201410353928 A CN201410353928 A CN 201410353928A CN 104135017 A CN104135017 A CN 104135017A
- Authority
- CN
- China
- Prior art keywords
- interconnection
- oscillation
- disconnection device
- installation place
- device installation
- 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
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的电压幅值积分初始值。
子系统的储能是子系统振荡剧烈程度的标尺,振荡越剧烈,则其储能数值越大;子系统的储能受到端口交换能量和内部供给能量以及本身耗散的影响。若将内部供给能量也视为能量供给元件的端口注入,则储能即为系统所有端口供给能量之和减去其本身耗散。若所有端口交换能量都为负,则对象系统内所有子系统的储能一定会减少到平衡点时的储能数值,意即振荡的平息。
将交换能量推至联络线,定义联络线上的振荡能量(Oscillation Energy Flow)为:
公式(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,PhasorMeasurement Unit)采样的连续设定次数据,计算每次采样后联络线上解列装置安装处的振荡能量。
设定次数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 (3)
1.一种基于振荡能量的失步解列判别方法,其特征是所述方法包括:
步骤1:判断系统是否发生故障,如果系统发生故障,则执行步骤2;否则,继续执行步骤1;
步骤2:根据联络线上解列装置安装处的连续设定次采样数据,计算每次采样后联络线上解列装置安装处的振荡能量;
步骤3:根据每次采样后联络线上解列装置安装处的振荡能量判断系统联络线是否发生失步振荡,如果系统联络线发生失步振荡,则执行步骤4;否则,返回步骤1;
步骤4:控制联络线解列装置进行解列。
2.根据权利要求1所述的方法,其特征是所述计算每次采样后联络线上解列装置安装处的振荡能量采用公式:
其中,OEF(t)为时刻t联络线上解列装置安装处振荡能量;
P为采样时联络线上解列装置安装处的有功功率;
Q为采样时联络线上解列装置安装处的无功功率;
U为采样时联络线上解列装置安装处的电压幅值;
为采样时联络线上解列装置安装处的电压幅值的微分;
θ为采样时联络线上解列装置安装处的电压相角;
为采样时联络线上解列装置安装处的电压相角的微分。
3.根据权利要求1或2所述的方法,其特征是所述根据每次采样后联络线上解列装置安装处的振荡能量判断系统联络线是否发生失步振荡具体为:
如果每次采样后联络线上解列装置安装处的振荡能量均为正数,则系统联络线未发生失步振荡;
如果每次采样后联络线上解列装置安装处的振荡能量均为负数,则系统联络线未发生失步振荡;
如果每次采样后联络线上解列装置安装处的振荡能量由正值变为负值,之后又变为正值并且稳定为正值,则系统联络线未发生失步振荡;
如果每次采样后联络线上解列装置安装处的振荡能量由负值变为正值,之后又变为负值并且稳定为负值,则系统联络线未发生失步振荡;
如果每次采样后联络线上解列装置安装处的振荡能量呈现两周以上的正值和负值交替变化,则系统联络线发生失步振荡。
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 true CN104135017A (zh) | 2014-11-05 |
CN104135017B 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) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104569820A (zh) * | 2015-01-15 | 2015-04-29 | 华南理工大学 | 一种基于fft变换的多机群失步解列判断方法 |
CN106556776A (zh) * | 2016-11-10 | 2017-04-05 | 华北电力大学 | 一种基于振荡能量的失步振荡中心捕获方法 |
CN107102192A (zh) * | 2017-04-28 | 2017-08-29 | 许继集团有限公司 | 一种全相异步振荡判别方法及装置 |
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 | 国家电网公司 | 一种具有自适应功能的失步解列方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104569820A (zh) * | 2015-01-15 | 2015-04-29 | 华南理工大学 | 一种基于fft变换的多机群失步解列判断方法 |
CN104569820B (zh) * | 2015-01-15 | 2017-07-28 | 华南理工大学 | 一种基于fft变换的多机群失步解列判断方法 |
CN106556776A (zh) * | 2016-11-10 | 2017-04-05 | 华北电力大学 | 一种基于振荡能量的失步振荡中心捕获方法 |
CN107102192A (zh) * | 2017-04-28 | 2017-08-29 | 许继集团有限公司 | 一种全相异步振荡判别方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104135017B (zh) | 2016-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103559347B (zh) | 一种大规模交直流电力系统电磁暂态仿真模型的构建方法 | |
CN102636728B (zh) | 基于力矩分解法识别负阻尼低频振荡的方法 | |
CN101794998A (zh) | 基于多机电力系统单一发电机电磁功率简洁表达方式的暂态稳定在线分析方法 | |
CN103606922B (zh) | 一种基于典型故障集的电力系统功角稳定性近似判定方法 | |
CN103094901B (zh) | 一种基于实测响应信息的暂态功角失稳实时判别的方法 | |
CN102611128B (zh) | 一种高压直流输电系统直流功率调制方法 | |
CN103632043A (zh) | 一种基于实测响应信息的电力系统主导失稳模式识别方法 | |
CN108023362A (zh) | 多回直流换相失败后系统的暂态稳定评价指标的计算方法 | |
CN104135017A (zh) | 一种基于振荡能量的失步解列判别方法 | |
CN108075478A (zh) | 一种暂态稳定紧急控制方法 | |
Heo et al. | An adaptive autoreclosure scheme with reference to transient stability for transmission lines | |
CN104505827B (zh) | 基于响应信息的复杂电力系统闭环控制方法 | |
CN104316827B (zh) | 一种电力系统振荡中心定位方法 | |
CN104124673B (zh) | 一种确定振荡中心位置的方法 | |
Padhi et al. | Solution of swing equation for transient stability analysis in dual-machine system | |
Zamani et al. | A PMU-augmented stability power limit assessment for reliable arming of special protection systems | |
Ilamparithi et al. | On the use of real lime simulators tor the lest and validation of protection and control systems of micro grids and smart grids | |
CN104331571A (zh) | 基于二阶非均匀Kuramoto模型频率同步理论的电力系统扰动极限切除时间方法 | |
Deng et al. | Real time transient instability detection based on trajectory characteristics and transient energy | |
Li et al. | New delay-dependent steady state stability analysis for WAMS assisted power system | |
CN106100125A (zh) | 基于广域信息的电网系统及电网稳定性分析与控制的方法 | |
Abdelaziz et al. | Phasor measurement units for out of step detection | |
Kyriakidis et al. | A transient stability assessment method using post-fault trajectories | |
Li et al. | Decentralized coordinated robust adaptive control for AC/DC interconnected power systems based on WAMS | |
Abiola et al. | Critical clearing time evaluation of Nigerian 330kV transmission system |
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 |