[发明专利]一种电力系统低频振荡模态辨识方法

说明书

技术领域

本发明涉及电力系统低频振荡分析技术领域，特别是一种基于改进EMD和数位演算法的电力系统低频振荡模态辨识方法。

背景技术

随着电力技术的不断进步，大电网互联已逐步实现，加之快速高放大倍数励磁装置的广泛使用，由此带来的低频振荡问题也越来越成为影响大规模电力系统稳定运行的重要因素，其不仅对电网运行安全上产生影响，同时也成为制约互联电网输送能力的瓶颈之一。如何通过广域测量数据分析电力系统低频振荡的局部动态行为，提取各个振荡模态的瞬时参数，进行广域阻尼控制是近年来的研究热点。

目前针对低频振荡问题所采用的分析方法主要有傅立叶变换、小波变换、Kalman滤波法、矩阵束辨识法、Prony算法和希尔伯特-黄(Hilbert-Huangtransform，HHT)算法等。傅立叶变换无法分析出阻尼特性和局部特性，所以不适于非线性、非平稳信号。小波变换存在频率交叠和自适应基选取问题，只适合于瞬态和非平稳信号。Kalman滤波法能消除噪声的影响，对不同输入信号的适应性较好，但计算精度和收敛速度受初始参数设置的影响很大。矩阵束辨识法能够准确估计系统的振荡模态，并具有较强的抗噪声能力，但若信号存在时变特性，该算法的计算误差较大，无法揭示振荡的动态特性。Prony方法虽然可以提取出振荡信号模式和阻尼等信息，但存在受噪声影响大、计算速度慢和定阶问题不确定等问题。HHT算法是近年发展起来的一种新型的适于非平稳、非线性信号的分析方法，传统的HHT算法受端点效应的影响，虽然能得到振荡模态的瞬时频率、瞬时幅值和衰减因子，难以达到较高的计算精度。但通过对该算法存在的端点效应问题进行改进，利用柯尔摩科洛夫-斯米洛夫检验(K-S)法剔除IMF分量中与低频振荡信号相似概率低的伪分量，可以精确提取出电力系统低频振荡的振荡模式和阻尼特性，有效的实现电力系统低频振荡信号的模态辨识。

发明内容

本发明的目的在于提供一种基于改进EMD和数位演算法的电力系统低频振荡模态辨识方法，该方法能够有效的改善EMD分解的端点效应，准确的剔除伪分量，精确的提取低频振荡信号模态瞬时参数，提高电力系统低频振荡模态辨识能力。

为实现上述目的，本发明的技术方案是：采用端点优化对称延拓法对EMD端点效应进行改进，利用改进之后的EMD对电力系统低频振荡信号进行分解，并利用K-S法对经EMD分解得到的IMF分量进行伪分量的剔除，利用数位演算法对剔除伪IMF分量后的有效IMF分量进行低频振荡模态瞬时参数的提取，得到电力系统低频振荡信号x(t)的瞬时频率、瞬时幅值和衰减因子等信息，其具体步骤如下：

步骤1：对滤波后的电力系统低频振荡信号x(t)进行离散化处理，采样步长取为Δt，得到离散化处理信号：

x_i＝[x₁,x₂,...,x_n]，i＝1,2,...,n

式中，x_i为离散化之后的信号值的集合，x₁,x₂,...,x_n为离散化之后t₁,t₂,...,t_n时刻的信号值；

步骤2：信号x_i两端点值用β、γ替代，构造新的数据序列x'_i：

x'_i＝[x'₀,x'_k,x'_n-1]

式中，x'₀＝β，x'_k＝x_k，x'_n-1＝γ，k＝1,2,...,n-2，i＝0,1,...,n-1；

步骤3：以端点x'₀、x'_n-1为中心对x'_i分别向两端进行对称延拓一个周期得到新的数据序列：

h＝[h_l,h_i,h_r]＝[2x'₀-x'_i,x'_i,2x'_n-1-x'_n-1-i]

式中，i＝0,1,...,n-1；