一种区域互联电网强迫功率振荡扰动源位置判断方法

申请号 CN201110390520.4 申请日 2011-12-01 公开(公告)号 CN102411118B 公开(公告)日 2013-06-26
申请人 武汉华中电力电网技术有限公司; 国家电网公司; 发明人 杨东俊; 李继升; 丁坚勇; 许汉平; 罗纯坚;
摘要 本 发明 公开了一种区域互联 电网 强迫功率振荡扰动源 位置 判断方法,涉及电 力 系统调度自动化及电网的安全运行。在区域电网中发生强迫功率振荡时,该方法通过计算区域互联电网中各区域能流方向因子的代数和,实时在线地判断出扰动源所在的位置。与传统的基于 能量 函数的扰动源 定位 方法相比,本发明提出基于计算能流方向因子的扰动源定位方法,能够减少周期性扰动分量和初始常数对支路 势能 非周期性分量判断的影响,具有更好的准确性。且减少了积分环节,简化了计算过程,能够更好的适应电网实时计算的要求。
权利要求

1.一种区域互联电网强迫功率振荡扰动源位置判断方法,其特征在于,该方法包含下列步骤:
a、在设有广域测量系统,由 个区域组成的交流区域互联电网中,区域间联络线支路有 条,各联络线支路的起始节点和终止节点装设同步相量测量单元PMU;
b、当广域测量系统监测到电网中振荡的主导振荡频率 、有功功率振荡幅度 、持续的周波数 均满足低频振荡告警阈值时,立即发出低频振荡告警;
c、广域测量系统发出低频振荡告警后,以 个振荡周波为一组,持续传递各条联络线支路的有功功率 ,起始节点频率 及主导振荡频率 的数据,其中:m=1,2,3,…,L;
d、将接收到的数据以 个振荡周波为一组,采用Prony分析方法,对第m条联络线支路的有功功率 、起始节点频率 数据进行参数辨识,得到主导振荡频率 模式下,第m条联络线支路的功率幅值 、起始节点频率振荡幅值 、功率初相位 、频率振荡初相位 ,其中:m=1,2,3,…,L;
e、将Prony分析方法辨识的 、 、 、 结果代入下式:

其中: 为基准频率, , 为电网基准频率;得到第m条联络线支路起始
节点的能流方向因子 ,其中:m=1,2,3,…,L;以功率流出节点方向为正,当能流方向因子时,扰动源位于第m条联络线支路起始节点所在的区域;当 时,扰动源位于第m条联络线支路终止节点区域;
f、求出 个区域中每个区域各条联络线支路的能流方向因子,如果第m条联络线支路起始节点位于该区域,能流方向因子取 ;如果第m条联络线支路终止节点位于该区域,能流方向因子取 ;其中:m=1,2,3,…,L;
g、求出每个区域的能流方向因子代数和,能流方向因子代数和最大的区域即为扰动源所在区域。

说明书全文

一种区域互联电网强迫功率振荡扰动源位置判断方法

技术领域

[0001] 本发明涉及电系统调度自动化及电网的安全运行,更具体涉及区域互联电网强迫功率振荡扰动源位置的判断。

背景技术

[0002] 电力系统在扰动下会发生电机转子间的相对摇摆并在缺乏阻尼的情况下引起持续振荡,振荡频率范围在0.1~2.5Hz,故称为低频振荡。低频振荡问题属于小扰动稳定范畴,随着互联电力系统的规模扩大、远距离重负荷输电系统的投入运行、快速自动励磁调节器和快速励磁系统的应用,国内外不少电力系统出现了低频振荡问题,低频振荡是影响电力系统的安全稳定运行的重要因素之一。
[0003] 强迫功率振荡是低频振荡的一种形式,即当系统受到外界持续周期性功率扰动的频率接近系统功率振荡的固有频率时,会引起大幅度的功率振荡。汤涌在电网技术2006,30(10) :29-33发表的《电力系统强迫功率振荡的基础理论》;杨东俊等在电力系统自动化
2009,33(23):24-28发表的《基于WAMS量测数据的低频振荡机理分析》;以及杨东俊等在电力系统自动化2011,35(10):99-103发表的《同步发电机非同期并网引起电力系统强迫功率振荡分析》等论文中,分别通过理论及电网实际案例分析,论证了对于强迫功率振荡最有效的处理措施是迅速找到并切除扰动源,但未提出找到扰动源的具体方法。
[0004] 广域测量系统由基于全球定位系统(GPS)的同步相量测量单元(PMU)及其通信系统组成,能够在广域电力系统中同步、高速采集机组和运行设备的有功功率、无功功率电压电流、相以及重要的开关信号,是一种能对电力系统动态过程进行监测和分析的工具。广域测量系统为电网的频振荡的监测、振荡事件分析和振荡预防及抑制等几个方面提供了新的技术手段。
[0005] 对于区域互联电网强迫功率振荡的扰动源定位问题,目前研究成果尚不多见,其中余一平等在电力系统自动化2010,34(5):1-6发表的《基于能量函数的强迫功率振荡扰动源定位》采用基于线性化的系统运动方程建立能量函数,借助系统中的能量转换特性识别强迫功率振荡扰动源所在位置。文献指出在多机系统中由外施扰动引发的强迫功率振荡所产生的能量,只有通过扰动源所在的机组才能注入系统,而被系统中各网络元件的阻尼所消耗,因此通过扰动源所在机组的势能变化与其他机组不同的特点,就可以找到扰动源。由于外施扰动做功注入系统的能量通过势能在网络中传播,并且系统中任一节点的势能守恒,即流入节点的势能之和等于流出节点的势能,发电机的势能可以用与其相连的支路势能表示。系统中i节点和j节点之间的支路为Li-j,其i端势能函数为:
[0006] (1)
[0007] 式中: 为支路Li-j电功率改变量, 为i节点相角偏移, 为支路i端角频率改变量,基准角频率 , 为电网基准频率,t为扰动时间。根据支路势能函数计算结果,若大于0则支路Li-j扰动能量从i节点流向j节点,扰动源在i节点方向;若小于0则支路Li-j扰动能量从j节点流向i节点,扰动源在j节点方向,据此判断扰动源所在的位置,经论证该方法可以适用于多机系统。但是,该方法需要通过积分求出支路势能函数,计算过程比较复杂不利于实时计算,并且在势能函数中主要通过非周期性分量来判断扰动源的位置,而势能函数中的周期性分量和积分初值的选取则会对判断的准确性产生一定的干扰。
[0008] 由于在强迫振荡的稳态阶段,各状态量都以扰动频率周期性的变化,对于电网中节点 和节点 之间的支路 ,设在 端观测到的振荡功率为 ,振荡频率 ,其中 、 分别为支路 功率和 端频率变化幅度,
、 为支路 功率和 端频率变化初相位, 为扰动频率;代入(1)式可得势能函数为:
[0009] (3)
[0010] 令:
[0011] (4)
[0012] (5)
[0013] (6)
[0014] 即有:
[0015] (7)
[0016] 可以看到势能函数是由一条幅值为 ,频率为 ,初相角为 的正弦曲线和一条斜率为 的直线叠加而成。由外施扰动注入系统并在网络中传播消耗的能量主要由非周期分量体现。定义势能流出节点 为正,流入节点 为负,由于 为不随时间变化的初始常数,通过斜率 即可判断支路势能的方向,当 时势能趋势由节点 流向节点 ,当时势能趋势由节点 流向节点 。由此, 表征了支路势能中非周期性分量的大小和流动方向,本发明定义 为能流方向因子。
[0017] 王强等在中国电力2001,11(34):38-41上发表的《Prony算法分析低频振荡的有效性研究》论证了Prony分析方法在电网低频振荡分析中的有效性。基于轨迹线的Prony分析方法直接对现场记录的数据波形进行信号分析,将实测轨线视为某些频率固定,幅值按指数规律变化的正弦信号(振荡模式)的线性组合,其模型表示为:
[0018] (2)
[0019] 其中:为振荡模式的个数;对于第 个振荡模式, 为振荡幅值; 为阻尼因子; 为振荡角频率; 为初始相位。从而将问题归纳为对各频率、幅值与阻尼系数的识别。Prony算法能根据采样得到的以 个振荡周波为一组的有功功率、节点频率及主导振荡频率数据进行参数辨识,估算得到主导振荡频率模式下的功率幅值、频率振荡幅值、功率初相位和频率振荡初相位。在确定振荡特征方面是一种较好的分析方法,在工程实践中被较为广泛应用。

发明内容

[0020] 本发明的目的是提供一种区域互联电网强迫功率振荡扰动源位置判断方法。在区域电网中发生强迫功率振荡时,该方法通过计算区域互联电网中各区域能流方向因子的代数和,实时在线地判断出扰动源所在的位置。本发明简单、易行,经实际案例检验,能够为电网调度人员快速查找和解决故障提供参考。
[0021] 为了达到上述目的,本发明采用如下技术方案:通过对广域测量系统实时测量的低频振荡有功和频率变化量,采用Prony分析方法进行参数辨识,根据公式(5)求解支路的能流方向因子,据此判断扰动源所在的位置。
[0022] 一种区域互联电网强迫功率振荡扰动源位置判断方法,该方法包含下列步骤:
[0023] a、在设有广域测量系统,由 个区域组成的交流区域互联电网中,区域间联络线支路有 条,各联络线支路的起始节点和终止节点装设同步相量测量单元PMU; [0024] b、当广域测量系统监测到电网中振荡的主导振荡频率 、有功功率振荡幅度、持续的周波数 均满足低频振荡告警阈值时,立即发出低频振荡告警;
[0025] c、广域测量系统发出低频振荡告警后,以 个振荡周波为一组,持续传递各条联络线支路的有功功率 ,起始节点频率 及主导振荡频率 的数据,其中:m=1,2,3,…,L;
[0026] d、将接收到的数据以 个振荡周波为一组,采用Prony分析方法,对第m条联络线支路的有功功率 、起始节点频率 数据进行参数辨识,得到主导振荡频率 模式下,第m条联络线支路的功率幅值 、起始节点频率振荡幅值 、功率初相位 、频率振荡初相位 ,其中:m=1,2,3,…,L;
[0027] e、将Prony分析方法辨识的 、 、 、 结果代入下式:
[0028] ,
[0029] 其中: 为基准角频率, , 为电网基准频率;得到第m条联络线支路起始节点的能流方向因子 ,其中:m=1,2,3,…,L;以功率流出节点方向为正,当能流方向因子 时,扰动源位于第m条联络线支路起始节点所在的区域;当 时,扰动源位于第m条联络线支路终止节点区域;
[0030] f、求出 个区域中每个区域各条联络线支路的能流方向因子,如果第m条联络线支路起始节点位于该区域,能流方向因子取 ;如果第m条联络线支路终止节点位于该区域,能流方向因子取 ;其中:m=1,2,3,…,L;
[0031] g、求出每个区域的能流方向因子代数和,能流方向因子代数和最大的区域即为扰动源所在区域。
[0032] 本发明具有以下优点:
[0033] 1、与传统的基于能量函数的扰动源定位方法相比,本发明提出基于计算能流方向因子的扰动源定位方法,能够减少周期性扰动分量和初始常数对支路势能非周期性分量判断的影响,具有更好的准确性。
[0034] 2、本方法减少了积分环节,简化了计算过程,能够更好的适应电网实时计算的要求。附图说明
[0035] 图1为一种区域互联电网强迫功率振荡扰动源位置判断方法的流程示意图。
[0036] 图2为区域互联电网地理位置示意图。

具体实施方式

[0037] 下面结合附图对本发明作进一步的说明。
[0038] 实施例
[0039] 本方法通过对广域测量系统实时测量的低频振荡有功和频率变化量,采用Prony分析方法进行参数辨识,根据公式(5)求解支路的能流方向因子,据此判断扰动源所在的位置。该方法包含下列步骤:
[0040] a、在有 个区域组成的交流区域互联电网中,区域间联络线支路有 条,设在其中一个区域存在一个周期性功率振荡的扰动源,引发整个区域互联电网产生强迫功率振荡,各区域间联络线至少一端装设有PMU,数据同步传输到电网调度中心的广域测量系统中。
[0041] b、当广域测量系统监测到电网中振荡的主导振荡频率 、有功功率振荡幅度、持续的周波数 均满足低频振荡告警阈值,即发出低频振荡告警。
[0042] c、在广域测量系统发出振荡告警后,以 个振荡周波为一组将广域测量系统采集到的有功功率,频率及主导振荡频率数据持续传递给分析程序。设:节点 和节点 分别为两区域间联络线支路 的起始节点和终止节点,则需要将联络线任一端(如 端)所采集联络线支路 个振荡周波的有功功率 ,端的频率 及主导振荡频率 数据持续传递给扰动源定位分析程序。
[0043] d、扰动源定位分析程序以 个振荡周波为一组数据,采用Prony分析方法对联络线支路 的有功功率 、端的频率 数据进行参数辨识,得到主导振荡频率 模式下联络线支路功率幅值 、端频率振荡幅值 、支路功率初相位 、频率振荡初相位,并依次对每条联络线一端数据进行分析。
[0044] e、将Prony分析方法辨识的 、 、 、 计算结果代入公式(5)求解各联络线的能流方向因子。对于节点 与节点 之间的联络线支路 ,以 端测量数据判断,当能流方向因子 时扰动源位于节点 所在区域侧,当 时扰动源位于节点 所在区域侧,分别计算各条联络线的能流方向因子 。
[0045] f、分区域求出一条联络线两端区域的能流方向因子。若测量节点 位于区域 ,对端节点 位于区域 ,则 , 。
[0046] g、区域能流方向因子为每一区域联络线能流方向因子的代数和,求出各区域的能流方向因子。如:与第 个区域相连的联络线有 条( ; ),则第 个区域能流方向因子为 ;能流方向因子代数和最大的区域即为扰动源所在区域,如:区域E的能流方向因子 ,则区域E即为扰动源所在区
域。
[0047] 实施例二
[0048] 一种区域互联电网强迫功率振荡扰动源位置判断方法,以附图2所示的5区域互联电网为例,说明本发明的实际应用过程。具体步骤如下:
[0049] a、在附图2所示的交流互联电网包括:S1、S2、S3、S4、S5共5个供电区域,分别通过4条区域间联络线:m1、m2、m3、m4将 S1和S3、S3和S5、S3和S4、S3和S2区域相连。联络线m1两端节点为i1和j1,m2两段节点为i2和j2,m3两段节点为i3和j3,m4两段节点为i4和j4;其中j1为区域S1中的节点,j2为区域S5中的节点,j3为区域S4中的节点,i4为区域S2中的节点,i1、i2、i3、j4为区域S3中的节点。设在区域S1中存在一个振荡频率为0.62Hz的周期性功率振荡的扰动源,引发整个区域互联电网产生强迫功率振荡,在节点i1、i2、i3、i4装设有PMU设备,振荡发生后数据数据同步传输到电网调度中心的广域测量系统。
[0050] b、设广域测量系统低频振荡的阈值需同时满足:主导振荡频率在0.1-2.5Hz,有功功率振荡幅度10MW,持续5个周波,当监测到电网中发生主导频率0.62Hz的功率振荡满足所设定的阈值,即发出低频振荡告警。
[0051] c、在广域测量系统发出振荡告警后,将节点i1、i2、i3、j4采集的联络线支路m1、m2、m3、m4的有功功率,频率及主导振荡频率 数据的以5个周波为一组持续传递给分析程序。
[0052] d、分析程序以5个振荡周波为一组数据,采用Prony分析方法对联络线m1、m2、m3、m4有功功率和频率数据进行参数辨识,得到主导振荡频率0.62Hz模式下联络线功率振荡幅值 、频率振荡幅值 、支路功率初相位 、频率振荡初相位 。
[0053] 表1 联络线振荡参数辨识结果
[0054]
[0055] e、将Prony分析方法辨识的 、 、 、 计算结果代入公式(5)求解能流方向因子 ,取电网基准频率 Hz,计算结果如下:
[0056] 表2 联络线能流方向因子计算结果
[0057]
[0058] 由于联络线m1、m2、m3、m4的广域测量数据取自节点i1、i2、i3、j4,节点i1、i2、i3属于区域S3,节点i4属于区域S2,因此,联络线m1能流方向因子 说明振荡能量从区域S1流向区域S3; 说明振荡能量从区域S3流向区域S5; 说明振荡能量从区域S3流向区域S4;联络线m4能流方向因子 说明振荡能量从区域S3流向区域S2。
[0059] f、分区域求出联络线m1、m2、m3、m4两端区域的能流方向因子,计算结果如表3所示:
[0060] 表3 联络线能流方向因子计算结果
[0061]
[0062] g、求出S1、S2、S3、S4、S5各区域的能流方向因子,计算结果如表4所示:
[0063] 表4区域能流方向因子计算结果
[0064]
[0065] 根据计算结果,区域S1为能流方向因子代数和最大的区域,即扰动源位于区域S1内。
[0066] 经案例证明,本发明能够实现区域互联电网强迫功率振荡扰动源位置的快速判断。
QQ群二维码
意见反馈