CN104993480A - 基于递推随机子空间的电力系统低频振荡在线辨识方法 - Google Patents
基于递推随机子空间的电力系统低频振荡在线辨识方法 Download PDFInfo
- Publication number
- CN104993480A CN104993480A CN201510432943.6A CN201510432943A CN104993480A CN 104993480 A CN104993480 A CN 104993480A CN 201510432943 A CN201510432943 A CN 201510432943A CN 104993480 A CN104993480 A CN 104993480A
- Authority
- CN
- China
- Prior art keywords
- matrix
- subspace
- identification
- frequency oscillation
- power system
- 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
- Complex Calculations (AREA)
Abstract
本发明涉及一种基于递推随机子空间的电力系统低频振荡在线辨识方法。本发明方法针对随机子空间辨识算法中要进行算法复杂度高的SVD分解,从而造成低频振荡模态辨识的实时性和动态辨识效果差这一问题,提出引入遗忘因子对Hankel协方差矩阵进行更新,利用投影近似子空间跟踪方法对子空间进行递推,避免SVD运算,显著减少计算复杂度,在每一次递归计算中,本发明算法的复杂度为 远小于SVD计算的复杂度,能够有效提高辨识的实时性,适合低频振荡模态的在线辨识,也可以为多时空尺度下电力系统的在线监测和稳定性分析提供有效支持。
Description
技术领域
本发明涉及电力系统低频振荡分析技术领域,特别是一种基于递推随机子空间的电力系统低频振荡在线辨识方法。
背景技术
随着大区电力系统互联以及大规模风、光发电系统并网运行且在电网装机容量中比重不断提升,电网弱阻尼低频振荡问题日益凸显,严重影响电网的稳定运行和电能质量控制。实时、准确辨识振荡各模态信息,定位振荡区域和机组,进而施加有效低频振荡抑制措施就显得极为关键。当前,广域测量系统(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:通过最小二乘法求取各模态的幅值和相角;对于个采样数据构造向量,模态幅值和相角辨识通过求解以下线性系统方程:
其中,,最小二乘解为:
由此,各振荡模态幅值和相位为:
步骤9:通过时间序列数据的更新,进行步骤4-8递推计算实现振荡模态信息的在线辨识。
相较于现有技术,本发明具有以下有益效果:
1、利用信号子空间和噪声子空间的正交性,递推随机子空间辨识方法能有效辨识出噪声环境下电力系统低频振荡主导模态:幅值、频率、相位、阻尼比、衰减因子、振型,具有很好的抗噪性能和辨识精度高;
2、引入遗忘因子对Hankel协方差矩阵进行更新,利用投影近似子空间跟踪方法对子空间进行递推,避免SVD运算,显著减少计算复杂度,在每一次递归计算中,本发明算法的复杂度为远小于SVD计算的复杂度,能够有效提高辨识的实时性,适合低频振荡模态的在线辨识,也可以为多时空尺度下电力系统的稳定性分析提供了有效支持。
附图说明
图1为本发明实施例的工作流程图。
图2是本发明方法对振荡模态辨识的效果图。
具体实施方式
下面结合附图,对本发明的技术方案进行具体说明。
本发明的一种基于递推随机子空间的电力系统低频振荡在线辨识方法,具体即电力系统低频振荡模态信息辨识方法方法结合图1进行说明并对电力系统低频振荡理想信号上叠加信噪比为10dB的白噪声,数据采样时间0-20秒,采样频率20Hz。采用本发明方法对此振荡模态的主要信息进行辨识。其中选择5s采样初始化数据构造Hankel矩阵,窗数据为0.1s。具体步骤如下:
步骤1:由采样离散时间序列数据构造协方差Hankel矩阵:
式中,为振荡电气量实测时间序列数据,,,为数据协方差,为数据窗的长度。
步骤2:对矩阵H做一次SVD分解计算:
将奇异值从大到小排列:
令,,和分别为左、右奇异向量,、和、分别为相应信号子空间和噪声子空间的左右奇异向量,是包含全部奇异值的对角矩阵,包含个主奇异值,为张成的信号子空间,为张成的噪声子空间。
步骤3:进行以下计算,其结果作为步骤4随机子空间迭代递归运算的初始值,其中符号“”表示广义逆。
其中, 为权值矩阵,为协方差矩阵估计;
步骤4:利用时间序列数据,广义可观测矩阵的计算更新可由以下基于投影近似子空间跟踪方法的递归计算实现:
其中遗忘因子可取值0.95-0.99。
步骤5:依据系统的最小实现理论,系统的广义可观测矩阵可表示为:
系统矩阵由下式求出:
上式中,和分别表示矩阵去除最后一行和第一行后的矩阵;符号“”表示广义逆。
步骤6: 求取矩阵的特征值,设采样时间间隔为,可辨识出各振荡模态的频率衰减因子和阻尼比为:
步骤7:随机子空间模态辨识方法由式(9-10)可以直接确定振荡模态振型:
上式中,和分别表示算法模型振型和系统振型,也是矩阵特征值对应的右特征向量。
步骤8:通过最小二乘法求取各模态的幅值和相角。对于个采样数据构造向量,模态幅值和相角辨识通过求解以下线性系统方程。
其中,最小二乘解为:
由此,各振荡模态幅值和相位为:
步骤9:通过时间序列数据的更新,通过步骤4-8递推计算实现了振荡模态信息的在线辨识。
为了让本领域人员更容易理解本发明,利用本发明提出的方法对下列振荡理想信号进行辨识,并叠加信噪比为10dB的白噪声,数据采样时间0-20秒,采样频率20Hz。图1给出了本发明所提出方法的流程图,图2给出了振荡理想信号模态信息辨识的仿真图。
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。
Claims (1)
1.一种基于递推随机子空间的电力系统低频振荡在线辨识方法,利用随机子空间方法辨识电力系统低频振荡模态信息,并引入遗忘因子对Hankel协方差矩阵更新,基于投影近似子空间跟踪方法对子空间进行递推,避免SVD运算,降低了算法复杂度、提升辨识实时性,其特征在于:包括如下步骤,
步骤1:由有限长度的振荡电气量实测数据构造Hankel协方差矩阵:
式中,为振荡电气量实测时间序列数据,,,为数据协方差,为数据窗的长度;
步骤2:对矩阵H做一次SVD分解计算:
将奇异值从大到小排列:
令,,和分别为左、右奇异向量,、和、分别为相应信号子空间和噪声子空间的左右奇异向量,是包含全部奇异值的对角矩阵,包含个主奇异值,为张成的信号子空间,为张成的噪声子空间;
步骤3:进行以下计算,其结果作为步骤4随机子空间迭代递归运算的初始值:
其中,符号表示广义逆;为权值矩阵,为协方差矩阵估计;
步骤4:利用时间序列数据、,广义可观测矩阵的计算更新可由以下基于投影近似子空间跟踪方法的递归计算实现:
其中,遗忘因子可取值0.95-0.99;
步骤5:依据系统的最小实现理论,系统的广义可观测矩阵可表示为:
其中,C为线性状态空间输出矩阵,为系统矩阵,矩阵由下式求出:
上式中,和分别表示矩阵去除最后一行和第一行后的矩阵;
步骤6:求取矩阵的特征值,设采样时间间隔为,可辨识出各振荡模态的频率衰减因子和阻尼比为:
步骤7:随机子空间模态辨识方法由下式可以直接确定振荡模态振型:
上式中,和分别表示算法模型振型和系统振型,也是矩阵特征值对应的右特征向量;
步骤8:通过最小二乘法求取各模态的幅值和相角;对于个采样数据构造向量,模态幅值和相角辨识通过求解以下线性系统方程:
其中,,最小二乘解为:
由此,各振荡模态幅值和相位为:
步骤9:通过时间序列数据的更新,进行步骤4-8递推计算实现振荡模态信息的在线辨识。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510432943.6A CN104993480B (zh) | 2015-07-22 | 2015-07-22 | 基于递推随机子空间的电力系统低频振荡在线辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510432943.6A CN104993480B (zh) | 2015-07-22 | 2015-07-22 | 基于递推随机子空间的电力系统低频振荡在线辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104993480A true CN104993480A (zh) | 2015-10-21 |
CN104993480B CN104993480B (zh) | 2017-03-08 |
Family
ID=54305255
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510432943.6A Active CN104993480B (zh) | 2015-07-22 | 2015-07-22 | 基于递推随机子空间的电力系统低频振荡在线辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104993480B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105678000A (zh) * | 2016-01-14 | 2016-06-15 | 上海交通大学 | 汽车横向动力学的子空间辨识方法 |
CN106786514A (zh) * | 2016-12-01 | 2017-05-31 | 西南交通大学 | 一种电力系统低频振荡模式在线辨识方法 |
CN107368125A (zh) * | 2017-09-08 | 2017-11-21 | 东北大学 | 一种基于cbr与rbr并行混合推理的高炉炉温控制系统及方法 |
CN107632522A (zh) * | 2017-08-31 | 2018-01-26 | 南京理工大学 | 一种质子交换膜燃料电池非线性状态空间模型辨识方法 |
WO2019080393A1 (zh) * | 2017-10-25 | 2019-05-02 | 东北电力大学 | 一种电力系统动态稳定协同辨识方法 |
CN111797500A (zh) * | 2020-06-02 | 2020-10-20 | 上海卫星工程研究所 | 基于规范变量分析与改进ssi的太阳电池阵模态辨识方法 |
CN113010844A (zh) * | 2021-03-09 | 2021-06-22 | 东北电力大学 | 一种基于子空间动态模式分解的参与因子计算方法 |
CN113036780A (zh) * | 2021-03-09 | 2021-06-25 | 东北电力大学 | 一种基于子空间动态模式分解的机电振荡参数辨识方法 |
CN113158785A (zh) * | 2021-03-11 | 2021-07-23 | 复旦大学 | 一种振荡信号模态参数的识别方法 |
CN113625555A (zh) * | 2021-06-30 | 2021-11-09 | 佛山科学技术学院 | 一种基于递推子空间辨识的自适应逆控制agv转速控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090204355A1 (en) * | 2006-06-27 | 2009-08-13 | Ata Engineering, Inc. | Methods and apparatus for modal parameter estimation |
CN103944174A (zh) * | 2014-04-28 | 2014-07-23 | 四川大学 | 基于互相关函数滤噪算法的低频振荡在线辨识方法 |
CN104242325A (zh) * | 2014-09-18 | 2014-12-24 | 国家电网公司 | 一种电力系统低频振荡模式参数辨识方法 |
-
2015
- 2015-07-22 CN CN201510432943.6A patent/CN104993480B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090204355A1 (en) * | 2006-06-27 | 2009-08-13 | Ata Engineering, Inc. | Methods and apparatus for modal parameter estimation |
CN103944174A (zh) * | 2014-04-28 | 2014-07-23 | 四川大学 | 基于互相关函数滤噪算法的低频振荡在线辨识方法 |
CN104242325A (zh) * | 2014-09-18 | 2014-12-24 | 国家电网公司 | 一种电力系统低频振荡模式参数辨识方法 |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105678000A (zh) * | 2016-01-14 | 2016-06-15 | 上海交通大学 | 汽车横向动力学的子空间辨识方法 |
CN106786514A (zh) * | 2016-12-01 | 2017-05-31 | 西南交通大学 | 一种电力系统低频振荡模式在线辨识方法 |
CN107632522B (zh) * | 2017-08-31 | 2020-06-19 | 南京理工大学 | 一种质子交换膜燃料电池非线性状态空间模型辨识方法 |
CN107632522A (zh) * | 2017-08-31 | 2018-01-26 | 南京理工大学 | 一种质子交换膜燃料电池非线性状态空间模型辨识方法 |
CN107368125B (zh) * | 2017-09-08 | 2019-06-25 | 东北大学 | 一种基于cbr与rbr并行混合推理的高炉炉温控制系统及方法 |
CN107368125A (zh) * | 2017-09-08 | 2017-11-21 | 东北大学 | 一种基于cbr与rbr并行混合推理的高炉炉温控制系统及方法 |
WO2019080393A1 (zh) * | 2017-10-25 | 2019-05-02 | 东北电力大学 | 一种电力系统动态稳定协同辨识方法 |
US11271401B2 (en) | 2017-10-25 | 2022-03-08 | Northeast Electric Power University | Synergic identification method for dynamic stability of power system |
CN111797500A (zh) * | 2020-06-02 | 2020-10-20 | 上海卫星工程研究所 | 基于规范变量分析与改进ssi的太阳电池阵模态辨识方法 |
CN111797500B (zh) * | 2020-06-02 | 2023-07-14 | 上海卫星工程研究所 | 基于规范变量分析与改进ssi的太阳电池阵模态辨识方法 |
CN113010844A (zh) * | 2021-03-09 | 2021-06-22 | 东北电力大学 | 一种基于子空间动态模式分解的参与因子计算方法 |
CN113036780A (zh) * | 2021-03-09 | 2021-06-25 | 东北电力大学 | 一种基于子空间动态模式分解的机电振荡参数辨识方法 |
CN113010844B (zh) * | 2021-03-09 | 2022-11-11 | 东北电力大学 | 一种基于子空间动态模式分解的参与因子计算方法 |
CN113158785A (zh) * | 2021-03-11 | 2021-07-23 | 复旦大学 | 一种振荡信号模态参数的识别方法 |
CN113625555A (zh) * | 2021-06-30 | 2021-11-09 | 佛山科学技术学院 | 一种基于递推子空间辨识的自适应逆控制agv转速控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104993480B (zh) | 2017-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104993480A (zh) | 基于递推随机子空间的电力系统低频振荡在线辨识方法 | |
Djurovic et al. | A hybrid CPF-HAF estimation of polynomial-phase signals: Detailed statistical analysis | |
Zhou et al. | Estimation of the dynamic states of synchronous machines using an extended particle filter | |
Xiao et al. | Dynamic tracking of low-frequency oscillations with improved Prony method in wide-area measurement system | |
Caro et al. | Multiple bad data identification considering measurement dependencies | |
Chen et al. | An adaptive TLS-ESPRIT algorithm based on an SG filter for analysis of low frequency oscillation in wide area measurement systems | |
Ning et al. | Two-level ambient oscillation modal estimation from synchrophasor measurements | |
CN108462191A (zh) | 一种基于o3kid算法的电力系统低频振荡辨识方法 | |
Mérigaud et al. | Incorporating ocean wave spectrum information in short-term free-surface elevation forecasting | |
CN106529700B (zh) | 一种基于时空相关的风电场预测模型构建方法和系统 | |
Gao et al. | A three-layer hybrid model for wind power prediction | |
Peng et al. | Adaptive sampling scheme for monitoring oscillations using Prony analysis | |
Phan et al. | Short-term Solar Power Forecasting Using XGBoost with Numerical Weather Prediction | |
Rashed et al. | False data detection in a clustered smart grid using unscented Kalman filter | |
Zhang et al. | Lodim: A novel power system state estimation method with dynamic measurement selection | |
CN111342478B (zh) | 一种基于最优变量投影的电力系统动态稳定评估方法 | |
Khandelwal et al. | Power system state estimation comparison of Kalman filters with a new approach | |
Yan et al. | Adaptive normalized lmp estimation for graph signal processing | |
Lu et al. | The iterative completion method of the spectrum map based on the difference of measurement values | |
Hassak et al. | Implementing a pmu based linear state estimator on a single board computer | |
CN110879311A (zh) | 一种计算电网功率的方法 | |
Lugnani et al. | Real-time Coherency Identification using a Window-Size-Based Recursive Typicality Data Analysis | |
Manoj et al. | Reduced-rank sigma-point Kalman filter and its application in ENSO model | |
Andersson et al. | On an iterative method for direction of arrival estimation using multiple frequencies | |
Jatoth et al. | A hybrid GA-adaptive particle swarm optimization based tuning of unscented Kalman filter for harmonic estimation |
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 |