[发明专利]基于递推随机子空间的电力系统低频振荡在线辨识方法

说明书

技术领域

本发明涉及电力系统低频振荡分析技术领域，特别是一种基于递推随机子空间的电力系统低频振荡在线辨识方法。

背景技术

随着大区电力系统互联以及大规模风、光发电系统并网运行且在电网装机容量中比重不断提升，电网弱阻尼低频振荡问题日益凸显，严重影响电网的稳定运行和电能质量控制。实时、准确辨识振荡各模态信息，定位振荡区域和机组，进而施加有效低频振荡抑制措施就显得极为关键。当前，广域测量系统(WAMS)的广泛应用，特别相量测量单元(PMU)采样频率和精度的不断提高，为低频振荡在线辨识提供了支持。

近年来，多种基于WAMS量测数据的系统模态辨识方法在电力系统低频振荡在线辨识中得到广泛应用和发展。Prony及其改进算法可准确提取主导振荡模态，获得低频振荡详细参数信息，但算法对噪声敏感，这一缺点使得算法定阶困难，计算过程极易产生虚假特征根或漏根，不利于工程运用。ARMA算法是一种基于线性模型的信号处理方法，具有一定的抗噪性能，但是在低频振荡模态辨识中，无法直接辨识振型，与Prony算法一样存在着模型定阶困难问题。希尔伯特-黄变换(HHT)适用于非线性、非平稳信号处理和辨识，但其理论框架尚需进一步完善，算法本身存在着模态混叠、边界处理等问题，在噪声干扰下经验模态(EMD)筛选结果也不尽人意。随机子空间辨识算法在振型辨识方面具有优势，是系统辨识领域最成功的算法之一。该算法直接以WAMS量测数据构成数据矩阵，利用信号子空间和噪声子空间的正交性，通过奇异值分解将信号空间分解为信号子空间和噪声子空间，具有很高的抗噪性能，是一种高精度的系统辨识方法。但该算法在每次辨识中都需要进行高阶数据矩阵奇异值分解(SVD)，算法费时间、冗余大，在电力系统低频振荡模态辨识应用中，不仅振荡模态动态辨识效果差，而且在线辨识工程实现困难。因此需要寻求能够显著减少算法复杂度、提高辨识实时性的方法，能够实现电力系统低频振荡模态信息的准确在线辨识，也可以为多时空尺度下电力系统的稳定性分析提供了有效支持。

发明内容

本发明的目的在于提供一种基于递推随机子空间的电力系统低频振荡在线辨识方法，具有良好的抗干扰能力和实时性，从而实现对电力系统低频振荡模态信息的在线辨识，为电力系统实时监测和稳定性分析提供有效支持。

为实现上述目的，本发明的技术方案是：一种基于递推随机子空间的电力系统低频振荡在线辨识方法，利用随机子空间方法辨识电力系统低频振荡模态信息，并引入遗忘因子对Hankel协方差矩阵更新，基于投影近似子空间跟踪方法对子空间进行递推，避免SVD运算，降低了算法复杂度、提升辨识实时性，包括如下步骤，

步骤1：由有限长度的振荡电气量实测数据构造Hankel协方差矩阵：

式中，为振荡电气量实测时间序列数据，，，为数据协方差，为数据窗的长度；

步骤2：对矩阵H做一次SVD分解计算：

将奇异值从大到小排列：

令，，和分别为左、右奇异向量，、和、分别为相应信号子空间和噪声子空间的左右奇异向量，是包含全部奇异值的对角矩阵，包含个主奇异值，为张成的信号子空间，为张成的噪声子空间；

步骤3：进行以下计算，其结果作为步骤4随机子空间迭代递归运算的初始值：

其中，符号表示广义逆；为权值矩阵，为协方差矩阵估计；

步骤4：利用时间序列数据，广义可观测矩阵的计算更新可由以下基于投影近似子空间跟踪方法的递归计算实现：

其中，遗忘因子可取值0.95-0.99；

步骤5:依据系统的最小实现理论，系统的广义可观测矩阵可表示为：

其中，C为线性状态空间输出矩阵，为系统矩阵，矩阵由下式求出：

上式中，和分别表示矩阵去除最后一行和第一行后的矩阵；

步骤6:求取矩阵的特征值，设采样时间间隔为，可辨识出各振荡模态的频率衰减因子和阻尼比为：

步骤7:随机子空间模态辨识方法由下式可以直接确定振荡模态振型：

上式中，和分别表示算法模型振型和系统振型，也是矩阵特征值对应的右特征向量；

步骤8:通过最小二乘法求取各模态的幅值和相角；对于个采样数据构造向量，模态幅值和相角辨识通过求解以下线性系统方程：

其中，，最小二乘解为：