基于动态模式分解的电力系统振荡的识别方法技术方案

一种基于动态模式分解的电力系统振荡识别方法，该方法针对严重影响电力系统安全稳定、电能质量的振荡问题，通过动态模式分解(Dynamic Mode Decomposition，DMD)方法提取系统振荡的模态信息，并建立相应的压缩降维动力学模型，分析计算特征值和特征向量，实现数据驱动下的系统振荡模式准确识别，并在时域上对动态特性进行重构。这种新型振荡识别方法适用于电力系统振荡模式识别研究，能够实现振荡的频率、幅值的识别，且具有良好的精确性。因此，本发明专利技术意义重大。本发明专利技术意义重大。本发明专利技术意义重大。

【技术实现步骤摘要】
基于动态模式分解的电力系统振荡的识别方法


[0001]本专利技术涉及电力系统，特别是一种基于动态模式分解的电力系统的振荡的识别方法。

技术介绍

[0002]自电力系统诞生以来，振荡就是其动态或稳定性研究的重要问题之一。现代电网在本质上是一个“强制”工作在50/60Hz(交流)和0Hz(直流)的电能系统。电力系统振荡问题通常是指在工作频率之外的或机械、或电磁或其耦合的循环往复的能量交换，这种能量交换严重时会危及电力系统的正常运行，将造成稳定性或电能的质量问题。
[0003]随着国家深化推进能源改革，电力系统的复杂度不断提高。新能源发电大规模并入电网，在电网中渗透率日益增加，同时，高压直流输电技术迅速发展，电网中还存在大量高压直流输电设备、动态无功补偿设备以及串联补偿设备。电力系统的输电方式和电网结构都发生重大变化，发电侧、负荷侧及输电网络之间的耦合形式也更加多样化。从而使电力系统振荡更加复杂，存在很多薄弱环节，将对电网安全稳定运行具有较大威胁性。因此，对电力系统振荡动态模式的准确识别变得愈加重要。
[0004]近年来，人们利用多种方法对电力系统振荡进行辨识分析。方法主要有两类，第一类是基于模型的分析方法，第二类是基于量测的分析方法。
[0005]特征值分析法是典型的基于模型的方法，它通过建立系统的状态空间方程，在平衡点线性化后得到系统的小信号模型，求解系统状态矩阵的特征值以判断系统的稳定性。该方法理论严谨、准确性高，但对于高阶系统，建模难度高、计算量大，无法表征振荡的非线性特性。
[0006]基于量测的分析方法不依赖模型，方法有快速傅里叶变换(FFT)、Prony算法、希尔伯特
‑
黄变换(HHT)等。FFT变换能直接获取信号的频率，但其余振荡信息，如衰减系数等无法得到，只能进行大致的定性分析。Prony算法通过指数函数的线性组合，对扰动信号数据进行等间隔采样的拟合，并由此获取信号中包含的振荡幅值、振荡相位、系统阻尼、扰动频率等信息，但是Prony算法在存在噪声的情况下，计算结果偏差较大。HHT算法包括两个步骤，首先利用经验模特分解获得有限数目的固有模态函数(IMF)，然后对IMF进行Hilbert变换，得到每一个IMF随时间变化的瞬时频率和瞬时幅值，拟合得到各个振荡参数。HHT算法能够分析非线性、非平稳信号，抗扰动性强，然而无法分辨频率靠近而能量差异大的两个分量，同时这种方法的计算负担相对较重。
[0007]综上所述，如何克服空间维度高难以计算、抗噪声性能差和无法分辨频率靠近而能量差异大的问题，实现电力系统振荡的准确识别，具有很重要的意义。

技术实现思路

[0008]针对上述问题，本专利技术提出一种基于动态模式分解的电力系统振荡识别方法。动态模式分解(Dynamic Mode Decomposition，DMD)算法不需要知道非线性系统的数学模型，
首先通过奇异值分解(Singular Value Decomposition，SVD)对海量数据的低维近似处理，进而准确获取其特征值和特征向量；然后，利用获得的特征值及特征向量对系统振荡模式进行准确识别，并对原非线性系统进行线性时域重构。可以有效解决上述问题。
[0009]本专利技术的技术解决方案如下：
[0010]一种基于动态模式分解的电力系统振荡识别方法，其特点在于，对电力系统状态变量数据进行采样，形成数据矩阵；
[0011]基于采样所得的数据矩阵，进行动态模式分解，并建立相应的压缩降维动力学模型，提取系统振荡的模态信息；
[0012]基于得到的动态模式分解结果，分析计算特征值和特征向量，实现数据驱动下的系统振荡模式准确识别，并在时域上对动态特性进行重构。
[0013]具有的步骤如下：
[0014]建立电力系统的离散状态方程。连续系统通常表示为：
[0015]dx/dt＝f(x,t,y)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0016]其中x表示n维状态变量，t表示时间，y表示系统的运行参数。一般来说，非线性系统很难求解。DMD算法不需要知道系统的方程，即f(x,t,y)是未知方程。将非线性系统构造为每个小区间内的近似线性系统，假设系统状态矩阵为A，式(1)可以等效为：
[0017]dx/dt＝Ax
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0018]计算过程中采用离散系统，对x以Δt为时间间隔进行m次采样得到数据矩阵X。
[0019][0020]3、对数据矩阵X进行动态模式分解，取数据矩阵X中的前m
‑
1个采样数据和后m
‑
1个采样数据分别构成两个新矩阵X1和X2，
[0021][0022][0023]系统状态方程的离散形式如式(6)所示，
[0024][0025]其中，A
d
是系统状态矩阵A的离散形式。是广义逆。如果系统规模较小，数据矩阵X维数较低，则离散状态矩阵A
d
可由(6)直接求解。但在日益复杂的现代电力系统中，A
d
通常是一个高维矩阵，直接求解广义逆是非常困难的。因此，可以通过考虑系统低阶近似矩阵的特征值分析来研究原系统的振荡特性。
[0026]通过奇异值分解(SVD)得到系统的低阶近似矩阵。首先，对数据矩阵进行奇异值分
解如式(7)所示：
[0027]X1≈UΣV
*
ꢀꢀꢀꢀꢀꢀ
(7)
[0028]其中*表示共轭转置，U和V是酉矩阵,有UU
*
＝I，VV
*
＝I，∑是对角矩阵，U∈C
n
×
r
，∑∈C
r
×
r
，V∈C
m
×
r
。r是X1奇异值分解并降阶后的秩，所以r＜n。联立式(6)和式(7)，得到：
[0029]A
d
＝X2VΣ
‑1U
*
ꢀꢀꢀꢀꢀꢀ
(8)
[0030]将式(8)等式的两侧同时左乘U
*
并右乘U得到：
[0031]U
*
A
d
U＝U
*
X2VΣ
‑1ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)
[0032]设A
′
d
＝U
*
A
d
U，矩阵A
′
d
是r维的，它是得到的低阶近似矩阵。进行特征值分解，结果如下：
[0033]A
′
d
Ω
′
＝Λ
′
Ω
′ꢀꢀꢀꢀꢀ
(10)
[0034本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于动态模式分解的电力系统振荡识别方法，其特征在于，包括步骤；对电力系统状态变量数据进行采样，形成数据矩阵；基于采样所得的数据矩阵，进行动态模式分解，并建立相应的压缩降维动力学模型，提取系统振荡的模态信息；基于得到的动态模式分解结果，分析计算特征值和特征向量，实现数据驱动下的系统振荡模式准确识别，并在时域上对动态特性进行重构。2.根据权利要求1所述的基于动态模式分解的电力系统振荡的识别方法，其特征在于，该方法包括下列步骤：步骤1)对电力系统状态变量数据进行采样：采用离散系统，对电力系统x以Δt为时间间隔进行m次采样：第i次采样值为x
1i
、x
2i
、、、、x
ni
；构建数据矩阵X：步骤2)基于数据矩阵X，进行动态模式分解：取数据矩阵X中的前m
‑
1个采样数据和后m
‑
1个采样数据分别构成两个新矩阵X1和X2，，系统状态方程的离散形式如式(6)所示，其中，A
d
是系统状态矩阵A的离散形式，是广义逆，如果系统规模较小，数据矩阵X维数较低，则离散状态矩阵A
d
可由(6)直接求解；通过奇异值分解得到系统的低阶近似矩阵：首先，对数据矩阵进行奇异值分解如式(7)所示：X1≈UΣV
*
ꢀꢀꢀꢀꢀꢀꢀ
(7)其中，*表示共轭转置，U和V是酉矩阵,有UU
*
＝I，VV
*
＝I，∑是对角矩阵，U∈C
n
×
r
，∑∈C
r
×
r
，V∈C
m
×
r
，r是X1奇异值分解并降阶后的秩，所以r＜n；求离散状态矩阵A
′
d
：联立式(6)和式(7)，得到：A
d
＝X2VΣ
‑1U
*
ꢀꢀꢀꢀꢀ
(8)
将式(8)等式的两侧同时左乘U
*
并右乘U得到：U
*
A
d
U＝U
*
X2VΣ
‑1ꢀꢀꢀꢀꢀ
(9)设A
′
d
＝U
*
A
d
U，矩阵A
′
d
是r维的，它是得到的低阶近似矩阵；特征值分解：A
d
′
Ω
′
＝Λ
′
Ω
′ꢀꢀꢀꢀ
(10)其中，矩阵Ω
′
的每列向量ξ
′
i
是特征向量，Λ
′
是包含相应特征值λ
i
的对角矩...

【专利技术属性】
技术研发人员：王茂海，张传云，陈志同，訾鹏，赵峰，袁峥，
申请(专利权)人：国家电网有限公司华北分部，
类型：发明
国别省市：