一种调速系统引发低频振荡的机理识别方法
【专利摘要】本发明公开了属于电力系统安全稳定领域的一种调速系统引发低频振荡的机理识别方法。该方法通过对电网有关参数和调速系统相关信号的测量，对系统进行初步模态分析，结合对电网功率振荡信号频率和波形的分析，形成低频振荡产生机理的判据。结合低频振荡产生机理的判据和电网功率信号与调速系统信号的振荡频率的一致性计算，形成判断调速系统低频振荡机理的方法。
【专利说明】
-种调速系统引发低频振荡的机理识别方法
技术领域
[0001] 本发明属于电力系统安全稳定领域，特别设及电力系统低频振荡机理分析。尤其 设及汽轮发电机调速系统引发低频振荡的机理识别方法，适用于生产实践。
【背景技术】
[0002] 互联电网可W提高电力系统运行的经济性，然而也导致电力系统动态稳定问题非 常突出，国内外多次发生过低频振荡现象。现代快速、高放大倍数励磁系统是引起负阻尼的 主要原因，加装电力系统稳定器PSS可W对运类低频振荡现象进行比较有效的抑制。随着调 速系统的快速发展，调速系统动态对电力系统动态稳定性的影响日益突出，系统中多次出 现由调速系统引发的低频振荡现象，严重影响的系统的安全稳定运行。
[0003] 现代大型汽轮机调速系统通常采用功频电液控制，其响应速度比W往机械液压式 调速系统大大提高，已经能够对电力系统的动态产生较大影响。调速系统包括转速测量机 构、控制器、电液转换器、油动机W及阀口控制等，某个环节参数设置不合理或者出现故障 均可能导致调节阀口摆动，进而造成电网低频功率振荡。
[0004] 目前，电力系统低频振荡的机理主要包括负阻尼机理和强迫振荡机理。负阻尼机 理是解释自由振荡稳定性的成熟理论，而强迫振荡是指系统在外施扰动源激励作用下产生 的振荡。对于汽轮发电机调速系统引起的振荡，究竟是负阻尼机理还是强迫振荡机理并没 有完善的判断方法。本发明专利通过对调速系统和电网有关参数的测量，形成判断调速系 统低频振荡机理的方法。
【发明内容】
[0005] 本发明针对调速系统引发的低频振荡，提出了一种调速系统引发低频振荡的机理 识别方法，其特征在于，所述方法包括
[0006] 步骤1、计算互联电网的各振荡模式，根据特征值求各振荡模式的振荡频率与阻尼 比;计算参与因子，参与因子Pki表示第i个模式中第k个变量的参与程度；
[0007] 步骤2、实测信号，信号包括电网功率、发电机转速、调节阀口指令、等效阀位、阀位 开度W及调节级压力；
[000引步骤3、提取实测信号的振荡频率及其衰减因子；
[0009] 步骤4、基于上述步骤和判据判断调速系统低频振荡的机理。
[0010] 所述步骤4的具体过程为
[0011] 步骤401、对系统进行初步模态分析；
[0012] 根据广域测量系统同步相量测量装置的监测数据，若是监测到至少两个区域的机 组群发生低频振荡现象，则可判断系统振荡模式为区域振荡模式，振荡产生机理为负阻尼 振荡;若是仅监测到单个发电厂内的机组发生低频振荡现象，则可初步判断系统振荡模式 为本地振荡模式，振荡产生机理须进一步判断；
[OOU]步骤402、判据Fi:电网功率信号振荡频率；
[0014] 由步骤I求得互联电网本地振荡模式的振荡频率，其中频率最小值为fmin，最大值 为fmax ;当电网功率信号的振荡频率f Pe在[fmin，fmax]范围W外时，判据Fl = 0，振荡产生机理 为强迫振荡；当时e在[fmin，fmax]范围W内时，判据Fl = I，振荡产生机理须进一步判断；
[0015] 步骤403、判据F2:电网功率信号振荡波形；
[0016] 在Fl = I的基础上，起振阶段为增幅振荡而稳态阶段为等幅振荡时，判据F2 = 0,振 荡产生机理为强迫振荡;起振阶段与稳态阶段均为增幅振荡时，判据F2=l，振荡产生机理 为负阻尼振荡；
[0017] 步骤404、判据F3:电网功率信号与调速系统信号的振荡频率的一致性；
[0018] 计算电网实测信号与调速系统实测信号的振荡频率方差：
[0019]
[0020] 其中，fPe为电网功率信号的振荡频率，为调节阀口指令的振荡频率，f。为发电 机转速信号的振荡频率，fsv为等效阀位信号的振荡频率，fgv为阀位开度信号的振荡频率， 时"为调节级压力信号的振荡频率;y为所有实测信号的平均振荡频率；
[0021] 当D(f)《0.005时，判据F3=l，判断该强迫振荡是由调速系统引起；当D(f)> 0.005时，判据F3 = 0,判断该强迫振荡不是由调速系统引起。
[0022] 有益效果
[0023] 本发明方法结合系统结构、运行参数W及广域测量装置对系统各信号的检测，形 成识别系统低频振荡产生机理的判据，能够判断系统低频振荡产生的原因，从而能够更有 针对性地采取控制措施抑制低频振荡。
【附图说明】
[0024] 图1为调速系统引发低频振荡的机理识别方法原理图；
[0025] 图2为电网功率信号的在线检测图；
[0026] 图3为调节阀口指令信号的在线检测图；
[0027] 图4为等效阀位信号的在线检测图；
[0028] 图5为阀位开度信号的在线检测图；
[0029] 图6为调节级压力信号的在线检测图；
[0030] 图7为发电机转速信号的在线检测图。
【具体实施方式】
[0031] 本发明提出了一种判断调速系统低频振荡机理的方法，图1为调速系统引发低频 振荡的机理识别方法原理图，具体步骤为
[0032] (1)计算互联电网的各振荡模式、参与因子；
[0033] (a)计算互联电网的振荡模式；
[0034] 线性系统的状态空间方程表示如下：
[0035]
(1)
[003引其中A X是系统的状态向量，A U是系统的输入向量，A y为输出向量。矩阵A与电网 系统的发电机参数和系统网络参数有关，通过求取矩阵A的特征值和特征向量，可W获得系 统模态。
[0037]系统的特征方程与特征向量求解如下：
[003引 det(A-人 1)=0 (2)
[0039] A4)i = Ai(l)i (3)
[0040] 由特征方程可W求得状态矩阵的特征值，其中每一对共辆复数特征值〇± jo对应 系统的一个振荡模式，并可根据特征值求得振荡的频率f与阻尼比C如下：
[0041] (4)
[0042] 皎）
[0043] (b)计算参与因子；
[0044] 由于式(2)中矩阵A为非对角阵，变量之间存在禪合，为解禪作如下变换，取A X = 巫A Z，其中O为矩阵A对应的右模态矩阵。取A Z = ^A X，其中4为矩阵A对应的左模态矩阵。 则参与矩阵P=[Pi P2，，，Pn]结合左模态矩阵与右模态矩阵，作为状态变量与模式之间的一 种度量。
[0045] 其中
[0046]
(6)
[0047] 参与因子Pki表示第i个模式中第k个变量的参与程度。
[0048] (2)电网功率、发电机转速、调节阀口指令、等效阀位、阀位开度、调节级压力等信 号的检测；
[0049] Prony算法是针对等间距采样点，假设模型是一些具有任意幅值、相位、频率和衰 减因子的指数函数的线性组合，即认为测量输入X(O),…，x(n-l)的估计值可W表示为：
[0050]
(7)
[ooy] 其中Am为幅值，fm为频率，am为衰减因子，0m为初化A t为时间间隔。
[0化2]算法的主要步骤为：
[0053] (a)构造样本函数矩阵；
[0056] (b)确定R阵的有效秩P;[0057] (C)根据式(9)求解a;
[0054； (8)
[0化5:
[005引
（9)
[0059] (d)求多项式(10)的根z;
[0060] l+aiz-i+...+SpZ-P = O (10)
[0061] 并根据式an递推出i(")
[0065] (f)利用式(13)计算振幅、相位、频率和衰减因子；
[006^ (U)
[0063]
[0064] (12)
[0066；
(13)
[0067] (3)实测信号的振荡频率和衰减因子提取；
[0068] W调速系统的电网功率信号为例，通过Prony算法对在线监测的电网功率信号进 行分析后，得到一组具有不同幅值、相位、频率和衰减因子的指数函数信号的线性组合。指 数信号的幅值Aj最大时对应的信号为主导振荡信号。取主导振荡信号的频率。取为振荡频 率，取主导振荡信号的衰减因子ai为振荡衰减因子，即
[0069] (14)
[0070] (巧）
[0071] 尼比由式(16)求得
[0072] 触
[0073] 同理可得调节阀口指令信号的振荡频率fuT，发电机转速信号的振荡频率f。、等效 阀位信号的振荡频率fsv、阀位开度信号的振荡频率fgv、调节级压力信号的振荡频率fpm。
[0074] (4)判断调速系统低频振荡的机理；
[0075] (a)对系统进行初步模态分析；
[0076] 根据广域测量系统(WAMS)同步相量测量装置(PMU)的监测数据，若是监测到两个 或多个区域的机组群发生低频振荡现象，则可判断系统振荡模式为区域振荡模式，振荡产 生机理为负阻尼振荡；若是仅监测到单个发电厂内的机组发生低频振荡现象，则可初步判 断系统振荡模式为本地振荡模式，振荡产生机理须进一步判断。
[0077] (b)判据Fi:电网功率信号振荡频率；
[0078] 在初步分析系统为本地振荡模式的基础上，进一步判断电网功率信号的振荡频率 fPe是否在本地振荡模式的频率范围内。由步骤（1)求得互联电网的各个振荡模式的振荡频 率，对本地模式振荡频率进行排序，其中本地振荡模式频率最大值为fmax，最小值为fmin。若 电网功率信号的振荡频率时e在[fmin，fmax]范围W外，判断振荡产生机理为强迫振荡;若电网 功率信号的振荡频率时e在[fmin，fmax]范围W内，振荡产生机理须进一步判断。
[0079] f Pe在[fmin，fmax ]范围W外时，Fl = 0，振荡产生机理为强迫振荡；
[0080] f Pe在[fmin，fmax ]范围W内时，Fl = 1，振荡产生机理须进一步判断。
[0081] (C)判据F2:电网功率信号振荡波形；
[0082] 当电网功率信号的振荡频率时e在本地模式范围内时，通过振荡信号的波形特征判 断振荡机理。若振荡波形的起振阶段为增幅振荡，阻尼呈现为负值，而稳态阶段为等幅振 荡，阻尼呈现为弱阻尼，振荡产生机理为强迫振荡;若振荡波形的起振阶段和稳态阶段均为 增幅振荡，阻尼一直呈现为负值，振荡产生机理为负阻尼振荡。
[0083] 起振阶段为增幅振荡而稳态阶段为等幅振荡时，F2 = 0,振荡产生机理为强迫振 荡；
[0084] 起振阶段与稳态阶段均为增幅振荡时，F2=1，振荡产生机理为负阻尼振荡。
[0085] (d)判据F3:电网功率信号与调速系统相关信号的振荡频率的一致性；
[0086] 在判断系统功率振荡产生机理为强迫振荡的基础上，利用电网功率信号与调速系 统信号的振荡频率的一致性来判断该强迫振荡是否由调速系统引起。
[0087] 方差是在概率论和统计方差衡量随机变量或一组数据时离散程度的度量，概率论 中方差用来度量随机变量和其数学期望之间的偏离程度，方差的计算如式(17)所示
[008引 (17)
[0089] 其中，W = E(X)
[0090] 令电网功率信号与各调速系统信号的振荡频率为变量，则频率的方差如式（18)所 示，用来度量电网功率的振荡频率与调速系统信号的振荡频率的一致性。
[0091]
侧
[0092] 其中，y = E(f);
[0093] 方差的大小与信号振荡频率的一致程度呈负相关性，即方差越小，则数据间的波 动越小，信号间的振荡频率一致性越高;方差越大，则数据间的波动越大，信号间的振荡频 率一致性越低。
[0094] D(f)《0.005时，判据F3= 1，判断该强迫振荡是由调速系统引起；
[00M] D( f) > 0.005时，判据F3 = 0，判断该强迫振荡不是由调速系统引起。
[0096]南方电网某电厂共有2台额定功率330MW的汽轮发电机组。2013年5月8日，机组检 修后并网运行，1号机组带220MW负荷运行，2号机组带230MW负荷运行，13:50分1号机组开始 进行单顺阀切换操作。阀切换过程中，1号机组发生低频功率振荡，有功功率在186MW~ 279MW范围内摆动，振荡持续77秒，振荡频率为0.171监，PMU/WAMS系统记录到了机组有功振 荡，其实测信号如图2所示。
[0097] (a)实测信号的振荡频率提取；
[0098] 图3-图7分别为调速系统的调节阀口指令、等效阀位、阀位开度、调节级压力、发 电机转速等实测信号图，通过Prony算法，可依次得出各实测信号的振荡频率。其中，fuT = 0.163监，f CV = 0.16細Z，f gv=0.163监，fpm=0.180HZ，f。= 0.164HZ。
[0099] (b)电网功率信号振荡频率；
[0100] 功率振荡的频率fPe为0.171监，电网本地振荡模式的最小振荡频率fmin为0 .細Z，最 小振荡频率fmax为1.2监。由于时e在[fmin，fmax]范围W外时，判据Fl = 0，振荡产生机理为强迫 振荡。
[0101] (C)电网功率信号与调速系统信号的振荡频率的一致性；
[0102] 根据公式（18)，对信号的振荡频率进行方差计算，求得实测信号振荡频率的期望 值
[0103] E(f)=0.168 监
[0104] 求得实测信号振荡频率的方差
[0105] D(f)=3.714X10-已
[0106] 由于0(門《0.005，判据。3 = 1，因此，可^判断本次低频振荡是由调速系统引起的 强迫振荡。
【主权项】
1. 一种调速系统引发低频振荡的机理识别方法，其特征在于，所述方法包括 步骤1、计算互联电网的各振荡模式，根据特征值求各振荡模式的振荡频率与阻尼比； 计算参与因子，参与因子Pki表示第i个模式中第k个变量的参与程度； 步骤2、实测信号，信号包括电网功率、发电机转速、调节阀口指令、等效阀位、阀位开度 W及调节级压力； 步骤3、提取实测信号的振荡频率及其衰减因子； 步骤4、基于上述步骤和判据判断调速系统低频振荡的机理。2. 根据权利要求1所述的一种调速系统引发低频振荡的机理识别方法，所述步骤4的具 体过程为 步骤401、对系统进行初步模态分析； 根据广域测量系统同步相量测量装置的监测数据，若是监测到至少两个区域的机组群 发生低频振荡现象，则可判断系统振荡模式为区域振荡模式，振荡产生机理为负阻尼振荡； 若是仅监测到单个发电厂内的机组发生低频振荡现象，则可初步判断系统振荡模式为本地 振荡模式，振荡产生机理须进一步判断； 步骤402、判据Fi:电网功率信号振荡频率； 由步骤1求得互联电网本地振荡模式的振荡频率，其中频率最小值为fmin，最大值为 fmax ;当电网功率信号的振荡频率f P庙[f min，fmax]范围W外时，判据Fl = 0，振荡产生机理为 强迫振荡；当时e在[fmin，fmax]范围W内时，判据Fl = l，振荡产生机理须进一步判断； 步骤403、判据F2 :电网功率信号振荡波形； 在Fl = l的基础上，起振阶段为增幅振荡而稳态阶段为等幅振荡时，判据F2 = 0,振荡产 生机理为强迫振荡;起振阶段与稳态阶段均为增幅振荡时，判据F2=l，振荡产生机理为负 阻尼振荡； 步骤404、判据F3:电网功率信号与调速系统信号的振荡频率的一致性； 计算电网实测信号与调速系统实测信号的振荡频率方差：其中，fPe为电网功率信号的振荡频率，为调节阀口指令的振荡频率，f。为发电机转速 信号的振荡频率，fsv为等效阀位信号的振荡频率，fgv为阀位开度信号的振荡频率，fpm为调 节级压力信号的振荡频率;μ为所有实测信号的平均振荡频率； 当D(f)《0.005时，判据F3 = l，判断该强迫振荡是由调速系统引起；当D(f)>0.005时， 判据F3 = 0,判断该强迫振荡不是由调速系统引起。
【文档编号】G01R31/08GK106099952SQ201610515402
【公开日】2016年11月9日
【申请日】2016年7月1日
【发明人】徐衍会, 伍双喜, 张莎, 吴国炳, 杨银国, 钱峰
【申请人】华北电力大学, 广东电网有限责任公司电力调度控制中心