一种基于相量的发电机强迫扰动源定位方法
【专利摘要】本发明公开了一种基于相量的发电机强迫扰动源定位方法，包括以下步骤：S1：根据电力系统振荡预警信号，选定功率振荡幅值超过阈值区域内的发电机组；S2：获取选定发电机组实测的输出电气功率、频率或转速、机端电流和励磁电流数据；S3：对实测数据进行预处理并获得主导强迫振荡模式，并计算主导强迫振荡模式下的相应电气波动量相量；S4：根据发电机相应电气波动量相量之间的关系进行发电机强迫扰动源的定位，并进一步确定扰动源所在的控制系统。本发明通过分析两大控制通道上的相关电气波动量的相量关系，准确定位扰动源所在的发电机组并且判断提供的强迫扰动是来自于调速系统还是励磁系统，这对于快速切除扰动源、有效抑制强迫振荡有重大意义。
【专利说明】
一种基于相量的发电机强迫扰动源定位方法
技术领域
[0001] 本发明涉及电力系统运行与控制技术领域，特别是涉及一种基于相量的发电机强 迫扰动源定位方法。
【背景技术】
[0002] 在大区电网互联的局面下，电网可靠性得到提高、资源配置优化的同时，低频振荡 问题也成为影响电力系统安全稳定运行的关键问题。当系统受到持续的周期性的扰动时， 扰动频率接近或与系统自然频率相等时，会引起大幅度的功率振荡，这就是强迫功率振荡 理论。该理论指出，强迫扰动所引起的响应不仅与电力系统本身的特性有关，也与扰动的变 化规律有关。强迫振荡具有明确的扰动源，起振快，表现为持续的等幅振荡，当扰动源切除 后，振幅迅速衰减。因此，对于强迫功率振荡最有效的抑制和处理措施就是快速准确地找到 扰动源并切除。
[0003] 发电机作为扰动源引起系统发生强迫振荡的可能性更高，引起的电网功率振荡幅 值也越大，对电网安全稳定运行的威胁更大。现有的扰动源定位方法主要是基于能量函数 来计算网络中振荡能量，根据能量流向能够定位到扰动源所在的区域以及相关机组。然而， 现有技术中大多数只能确定到扰动源存在的发电机组，但对于发电机中具体哪一控制系统 产生的扰动引发功率振荡却无法判断。为了确定发电机内哪种控制设备产生的扰动，需要 进行大量的推理、排查和仿真工作，增加了实现快速切除扰动源，抑制强迫功率振荡的复杂 程度。
【发明内容】
[0004] 发明目的:本发明的目的是提供一种能够解决现有技术中存在的缺陷的基于相量 的发电机强迫扰动源定位方法。
[0005] 技术方案:为达到此目的，本发明采用以下技术方案：
[0006] 本发明所述的基于相量的发电机强迫扰动源定位方法，包括以下步骤：
[0007] SI:根据电力系统振荡预警信号，选定功率振荡幅值超过阈值区域内的发电机组；
[0008] S2:获取选定发电机组实测的输出电气功率、频率或转速、机端电流和励磁电流数 据；
[0009] S3:对实测数据进行预处理并获得主导强迫振荡模式，并计算主导强迫振荡模式 下的相应电气波动量相量；
[0010] S4:根据发电机相应电气波动量相量之间的关系进行发电机强迫扰动源的定位， 并进一步确定扰动源所在的控制系统:如果机组输出电气功率波动与频率波动相位相差小 于90度，则扰动源位于该机组且由调速系统提供机械功率扰动；如果机端电流波动量超前 于励磁电流波动，则扰动源位于该机组且由励磁系统提供励磁扰动;如果以上两个判据都 不满足，则该机组不是扰动源。
[0011] 进一步，所述步骤Sl包括以下的步骤：
[0012] SI. 1:确定功率振荡预警阙值；
[0013] SI.2:电力系统发生强迫振荡时，系统根据阙值发出预警信号，并选取振荡明显区 域内的相应发电机组。
[0014] 进一步，所述步骤S3中，对实测数据进行的预处理包括剔除异常数据、填补数据、 滤波、去直流以及对称化处理;采用TLS-ESPRIT算法获取主导强迫振荡模式下的相应电气 波动量相量。
[0015] 进一步，所述步骤S4中的判据通过以下分析得到：
[0016] S4.1:当单机无穷大系统发生强迫功率振荡时，发电机的微分方程增量化后可表 示为：
[0017] (1)
[0018] 式（1)中，Td'o为发电机d轴短路暂态时间常数，M为发电机惯性常量，D为发电机阻 尼转矩系数，AEq'为电机q轴瞬变电动势波动，AE f为定子励磁电动势波动，AEqS空载电 动势波动，A S为发电机内节点电压相角波动，A ω为转速波动，APmS电机机械功率波动， A Pe3为输出电气功率波动，ρ = d/dt为对时间的导数算子；
[0019 ] S4.2:记调速器的传递函数为Gcqv (s )、励磁系统的传递函数为Ge (s)及电力系统稳 定PSS的传递函数为Gpss(S),三种判据分别由如下的方法得到：
[0020] 1)当发电机存在机械功率扰动时，此时励磁系统正常工作，如式(2)所示：
[0021]
[0022]其中，Ke+sDe为发电机存在机械功率扰动时的发电机复转矩，Ke和De都为正值，
:为机端电压，Te为一阶励磁 系统的时间常数，Ke为一阶励磁系统的放大倍数；
[0023] 因此，Δ Pdg位介于Δ δ和Δ ω之间，g卩Δ ?0与Δ ω之间相位相差小于90度；
[0024]
(3)
[0025] 由式(3)可知，Δ E'q超前于Δ Pe，Δ Eq超前于Δ ?0且滞后于Δ E'q，因此Δ E'q相位超 前于AEq且小于90度；
[0026] 2)当发电机存在励磁扰动时，此时可得：
[0027]
[0028] 其中，K'e+D'e为发电机存在励磁扰动时的发电机复转矩，此时调速系统正常工作， 可得：
[0029] APe= APm-DA ω-sM Δ ω =Gpe(S) Δ ω (5)
[0030] 其中，Gpe(S) =Ggqv(s)-D，由于Ggqv(S)具有负实部，因此Gpe(S)实部为负，Δ ?0与Δ ω之间相位相差大于90度；
[0031] 由于APdg位介于Δδ和-Δ ω之间，根据式（3)可得，AE'q滞后于APe，ΔEq滞后 于A Pe且超前于Δ E ' q，因此△ Eq相位超前于Δ E ' q且小于90度；
[0032] 3)当发电机不是扰动源时，此时可得：
[0033] APe= APeO+APei (6)
[0034] 其中，APe3Q即式⑵中的APe3，是发电机控制系统产生的电气功率波动；APe3l是系 统发生强迫振荡，施加给发电机的一个电气功率波动强迫项；因此，根据式(5)可知，^匕与 A ω之间相位相差大于90度;根据式⑶可知，△ E 'q相位超前于△ Eq且小于90度；
[0035] S4.3:将步骤S4.2分析得出的三种相量关系转换为易直接测量的电气波动量之间 的相量关系，并应用于发电机强迫扰动源定位:如果机组输出电气功率波动与频率波动相 位相差小于90度，则扰动源位于该机组且由调速系统提供机械功率扰动;如果机端电流波 动量超前于励磁电流波动，则扰动源位于该机组且由励磁系统提供励磁扰动;如果以上两 个判据都不满足，则该机组不是扰动源。
[0036] 有益效果:与现有技术相比，本发明具有如下的有益效果：
[0037] (1)分析两大控制通道上的相关电气波动量的相量关系，准确定位扰动源所在的 发电机组并且判断提供的强迫扰动是来自于调速系统还是励磁系统，这对于快速切除扰动 源、有效抑制强迫振荡有重大意义；
[0038] (2)采用TLS-ESPRIT算法进行识别分析，提取出可用于扰动定位的主导频率下的 相量信息，能够减少异常数据以及与强迫扰动无关分量对判别结果的影响；
[0039] (3)由于具体到电网中某台机组进行分析，所需数据皆可以在发电机组本地测量， 不需要远方的实测数据，避免引入全网模型导致分析误差变大，且只需实测数据信息无需 进行额外计算，提高了定位的准确性；
[0040] (4)为电力系统安全稳定运行提供保障，能够快速定位系统扰动源，减少了低频振 荡带来的经济损失及社会问题，具有经济和社会效益。
【附图说明】
[0041] 图1为本发明的流程图；
[0042]图2为本发明的单机无穷大系统振荡传递函数图；
[0043]图3是本发明的WAMS实测数据预处理并获取相量信息流程图；
[0044]图4是本发明的发电机内电气波动相量关系转换图；
[0045]图5是本发明【具体实施方式】的四机两区系统接线图。
【具体实施方式】
[0046] 下面结合【具体实施方式】对本发明的技术方案作进一步的介绍。
[0047] 下面以四机两区系统发生强迫振荡的扰动源定位为例，说明本发明的具体实施方 式，该系统接线图如图5所示。其中，发电机内控制系统的传递函数框图如图2所示，调速器 的传递函数为Gcov(S)、励磁系统的传递函数为Ge(S)及电力系统稳定PSS的传递函数为G pss (s )，系数Ki~K6是和系统结构、参数、运行工况有关的常数，通常为正值。其中，
发电机经过小干 扰稳定分析可知，该系统包含1个区间振荡模式，频率为〇. 64Hz。首先对该模式的强相关机 组Gl额外施加持续性的机械功率扰动，扰动频率为0.64Hz，扰动幅值为0.0 Ip. u .，发生强迫 振荡后进行扰动源定位;接着对该模式的强相关机组Gl额外施加持续性的励磁电压扰动， 扰动频率为0.64Hz，扰动幅值为0.8p. u.，发生强迫振荡后进行扰动源定位。本发明公开了 一种基于相量的发电机强迫扰动源定位方法，如图1所示，包括以下步骤：
[0048] SI:根据电力系统振荡预警信号，选定功率振荡幅值超过阈值区域内的发电机组； [0049 ] S2:获取选定发电机组实测的输出电气功率、频率或转速、机端电流和励磁电流数 据；
[0050] S3:对实测数据进行预处理并获得主导强迫振荡模式，并计算主导强迫振荡模式 下的相应电气波动量相量;与扰动相关的振荡频率下的相量信息通过以下步骤得到：
[0051] S3.1:对实测数据进行剔除、填补处理，并且经过滤波、去直流以及对称化，数据处 理流程图如图3所不；
[0052] 33.2:通过11^45?1?11'算法计算出预处理后数据中各分量的频率大小、衰减系数 及阻尼比；
[0053] S3.3:根据强迫振荡过程中波动稳态量不衰减特点，确定强迫振荡的主导频率，即 共振频率；
[0054] S3.4:通过最小二乘法获取共振频率下，电气波动相量的幅值和相位；
[0055] 以上方法能够滤去异常数据以及与振荡模式无关的高次分量的影响，使所获得数 据更加理想，不会对其相量信息造成影响，采用TLS-ESPRIT获取主导振荡频率下的相量幅 值和相位信息，使结果更加准确。具体实施例的数据分析如下：
[0056] 1)当发电机调速系统提供机械功率扰动发生强迫振荡时，发电机相关电气波动量 相位信息如表1所示：
[0057] 表1主导振荡模式相位参数
[0059] 2)当发电机励磁系统提供励磁电压扰动发生强迫振荡时，发电机相关电气波动量 相位信息如表2所示：
[0060]表2主导振荡模式相位参数
[0062] S4:根据发电机相应电气波动量相量之间的关系进行发电机强迫扰动源的定位， 并进一步确定扰动源所在的控制系统。
[0063] 根据八？0与八ω以及之间的相量关系得出的判据，将之转换为易直接 测量的电气波动量之间的相量关系，转换过程包括：
[0064] a.发电机角频率的波动量相量△ ω与机端母线电压频率波动量相量△ f存在Δ ω =2π △ f这一关系，因此发电机角频率与输出电气有功功率之间的相量关系可转换为机端 电压频率与输出电气有功功率之间的相量关系。
[0065] b.发电机内部各电气量存在关系式：
[0066]
(7)
[0067] 其中xd，x'd分别为发电机的d轴同步电抗和d轴瞬态同步电抗且xd-x'd大于0, A id 为d轴电流波动，xadS发电机d轴电枢反应电抗，△ if为励磁电流波动。根据相量关系，可将 AE'q与AEq之间的相量关系转换为Δ id与Δ if之间的相量关系，相量关系图如图4所不。
[0068] 另外，发电机d轴电流波动量△ id与发电机机端电流波动量A I同相位，因此△ id与 Δ if之间的相量关系转化为△ I与Δ if之间的相量关系。
[0069] 综上，如果机组输出电气功率波动相位辟％与频率波动相位你相差小于90度，即 Δ炉=I卿e-得I <90°，则扰动源位于该机组且由调速系统提供机械功率扰动；如果机端电流 波动量相位件超前于励磁电流波动量相位物'，即夠>0°，则扰动源位于该机组且由励 磁系统提供励磁扰动；以上两个判据都不满足，则该机组不是扰动源。
[0070] 1)当发电机调速系统提供机械功率扰动发生强迫振荡时，根据表1可知，机组Gl满 足式(5)而不满足式(7)，且其他三台机组均不满足于，因此可以得出，机组Gl是扰动源且调 速系统提供机械功率扰动。这与实际情况相符合。
[0071] 2)当发电机励磁系统提供励磁电压扰动发生强迫振荡时，根据表2可知，机组Gl满 足式(7)而不满足式(5)，且其他三台机组均不满足于，因此可以得出，机组Gl是扰动源且励 磁系统提供励磁电压扰动。这与实际情况相符合。
[0072] 以上分析结果与实际情况皆相符合，验证了本发明的有效性及实用性。
[0073] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人 员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形 也应视为本发明的保护范围。
【主权项】
1. 一种基于相量的发电机强迫扰动源定位方法，其特征在于:包括以下步骤： S1:根据电力系统振荡预警信号，选定功率振荡幅值超过阈值区域内的发电机组； S2:获取选定发电机组实测的输出电气功率、频率或转速、机端电流和励磁电流数据； S3:对实测数据进行预处理并获得主导强迫振荡模式，并计算主导强迫振荡模式下的 相应电气波动量相量； S4:根据发电机相应电气波动量相量之间的关系进行发电机强迫扰动源的定位，并进 一步确定扰动源所在的控制系统:如果机组输出电气功率波动与频率波动相位相差小于90 度，则扰动源位于该机组且由调速系统提供机械功率扰动;如果机端电流波动量超前于励 磁电流波动，则扰动源位于该机组且由励磁系统提供励磁扰动;如果以上两个判据都不满 足，则该机组不是扰动源。2. 根据权利要求1所述的一种基于相量的发电机强迫扰动源定位方法，其特征在于:所 述步骤S1包括以下的步骤： S1.1:确定功率振荡预警阙值； S1.2:电力系统发生强迫振荡时，系统根据阙值发出预警信号，并选取振荡明显区域内 的相应发电机组。3. 根据权利要求1所述的一种基于相量的发电机强迫扰动源定位方法，其特征在于:所 述步骤S3中，对实测数据进行的预处理包括剔除异常数据、填补数据、滤波、去直流以及对 称化处理;采用TLS-ESPRIT算法获取主导强迫振荡模式下的相应电气波动量相量。4. 根据权利要求1所述的一种基于相量的发电机强迫扰动源定位方法，其特征在于:所 述步骤S4中的判据通过以下分析得到： S4.1:当单机无穷大系统发生强迫功率振荡时，发电机的微分方程增量化后可表示为：式（1)中，f do为发电机d轴短路暂态时间常数，Μ为发电机惯性常量，D为发电机阻尼转 矩系数，AEq'为电机q轴瞬变电动势波动，AEf为定子励磁电动势波动，AEq为空载电动势 波动，△ δ为发电机内节点电压相角波动，△ ω为转速波动，△ Pm*电机机械功率波动，Δ Pe 为输出电气功率波动，P = d/dt为对时间的导数算子； S4.2:记调速器的传递函数为GCQV(s)、励磁系统的传递函数为Ge(s)及电力系统稳定PSS 的传递函数为GPSS(s)，三种判据分别由如下的方法得到： 1)当发电机存在机械功率扰动时，此时励磁系统正常工作，如式(2)所示：其中，Ke+sDe为发电机存在机械功率扰动时的发电机复转矩，Ke和D e都为正值， ，ut为机端电压，TE为一阶励磁系统的时间常数，KE为一阶励磁系统的放大倍数；因此，APjg位介于Δδ和Δ ω之间，即八？(3与八ω之间相位相差小于90度；由式(3)可知，Δ E'q超前于Δ Pe，Δ Eq超前于Δ ?0且滞后于Δ E'q，因此Δ E'q相位超前于 A Eq且小于90度； 2) 当发电机存在励磁扰动时，此时可得：其中，K'e+D'e为发电机存在励磁扰动时的发电机复转矩，此时调速系统正常工作，可 得： Δ Pe= Δ Pm-D Δ ω -sM Δ ω =Gpe(s) Δ ω (5) 其中，GPe3(s)=GGQV(s)-D，由于Ggqv(s)具有负实部，因此G Pe(s)实部为负，ΔΡΑΔ ω之 间相位相差大于90度； 由于Δ Pjg位介于Δ δ和-Δ ω之间，根据式(3)可得，Δ E'q滞后于Δ Pe，Δ Eq滞后于Δ Pe 且超前于Δ E ' q，因此Δ Eq相位超前于Δ E ' q且小于90度； 3) 当发电机不是扰动源时，此时可得： Δ Pe= Δ Pe〇+ Δ Pel (6) 其中，Δ 即式(2)中的Δ pe，是发电机控制系统产生的电气功率波动；Δ pel是系统发 生强迫振荡，施加给发电机的一个电气功率波动强迫项；因此，根据式(5)可知，△?(3与八ω 之间相位相差大于90度;根据式⑶可知，△ Ε ' q相位超前于△ Eq且小于90度； S4.3:将步骤S4.2分析得出的三种相量关系转换为易直接测量的电气波动量之间的相 量关系，并应用于发电机强迫扰动源定位:如果机组输出电气功率波动与频率波动相位相 差小于90度，则扰动源位于该机组且由调速系统提供机械功率扰动；如果机端电流波动量 超前于励磁电流波动，则扰动源位于该机组且由励磁系统提供励磁扰动;如果以上两个判 据都不满足，则该机组不是扰动源。
【文档编号】H02J3/24GK106058897SQ201610601468
【公开日】2016年10月26日
【申请日】2016年7月28日 公开号201610601468.5, CN 106058897 A, CN 106058897A, CN 201610601468, CN-A-106058897, CN106058897 A, CN106058897A, CN201610601468, CN201610601468.5
【发明人】蒋平, 郑斌青, 冯双
【申请人】东南大学