Prony算法在电力系统小干扰稳定分析中应用.docVIP

Prony算法在电力系统小干扰稳定分析中应用 　　摘要：随着大规模电力系统中区域电网之间的互联，传统用于电力系统小干扰稳定分析的特征值分析方法受到系统规模及完整参数获取困难等因素的限制。本文介绍了一种不依赖于数学模型，用于分析和研究电力系统低频振荡的有力工具――Prony算法。在阐述Prony算法基本原理的基础上，结合电力系统实际，探讨Prony算法在实际应用时应如何确定模型阶次、采样步长、数据长度及数据预处理，并给出了将Prony辨识的结果应用于小干扰稳定分析的方法。 　　关键词：Prony算法；小干扰稳定分析；时域拟合精度；能量级 　　作者简介：蔡燕春（1984-），男，湖北浠水人，广东电网公司广州供电局调度通信中心，助理工程师；姚煜（1985-），男，江西抚州人，广东电网公司广州供电局调度通信中心，助理工程师。（广东#8194;广州#8194;510620） 　　中图分类号：TM711#8195;#8195;#8195;#8195;#8195;文献标识码：A#8195;#8195;#8195;#8195;#8195;文章编号：1007-0079（2011）06-0164-03 　　 　　随着大规模电力系统中区域电网之间的互联，电力系统低频振荡问题日益突出。传统的电力系统小干扰稳定分析的方法，通常是在稳态运行点将描述电力系统动态特性的微分-代数方程线性化，消去代数方程，得到常系数线性微分方程组： 　　（1） 　　其中为描述电力系统特性的状态向量，由各发电机的转角、转速等状态变量组成。因此，分析电力系统小干扰稳定特性转化为对A矩阵进行特征值分析，以判断整个系统的小干扰稳定特性。随着电网规模的不断扩大，用来描述电力系统特性的状态向量也不断增多，这直接导致用于稳定分析的常微分方程组的阶数很大。实际应用中，对一个中等规模的省网而言，在外部网络等值的情况下，其常微分方程组的阶数已达到数千阶。大型电网A矩阵的特征值分析会存在“维数灾”问题，研究者对其进行了一定的研究，但仍没有很好的解决方法。 　　Prony分析是近年来得到关注和应用的一种信号处理方法，它可以通过给定输入信号下的响应直接估计系统的振荡频率、衰减、幅值和相对相位。由于它不依赖于数学模型，也就不存在“维数灾”问题。Prony算法已被引入电力系统低频振荡辨识，并实际应用于电力系统如南方电网进行低频振荡分析。[1]本文结合电力系统的实际情况，全面介绍了Prony算法数学原理及实现过程，探讨算法在实际应用时如何确定模型阶次、采样步长、数据长度及数据预处理以及分析评价Prony辨识的结果。 　　一、Prony算法介绍 　　1795年，法国数学家Prony提出了用复指数函数的一种线性组合来描述等间距采样数据的数学模型，从而形成了Prony算法。[2]Prony算法是一种能够根据采样值直接估算出信号频率、衰减、幅值和初相位的分析方法。其数学描述如下：令作为测量数据的模型，且#8195;n=0，1……N-1 （2） 　　它假设等间距采样点的估计值是一系列具有任意振幅、相位、频率和衰减因子的指数函数的线性组合。式中p为Prony模型的阶数，N为采样数据点的个数，N≥2p，bi=Aiexp(jθi)，其中Ai为信号幅值；θi为初相；zi=exp[(αi+j2πfi)△t]，αi为衰减因子；fi表示振荡频率；△t代表采样间隔。 　　Prony方法的关键是认识到式（2）的拟合是一常系数线性差分方程的齐次解，通过一系列的数学变换，可推出差分方程式： 　　 n=0,1,……N-1（3） 　　定义，i, j=0,1,……p，根据最小二乘法，式（3）的法方程如式（4）所示： 　　（4） 　　解此法方程，即可得到系数a1、a2……ap的估计值，将其代入式（3）的特征方程：1+a1z-1+…+apz-p=0（5） 　　解方程（5）求解出式（2）中的Z值。于是，数学模型（2）可简化为： 　　（6） 　　求解方程（6）并利用式（7）即可算出幅值Ai、相位θi、频率fi和衰减因子αi。 　　（7） 　　需要指出的是，根据常系数线性微分方程组的知识，可得式（1）中状态向量的时域解： 　　#8195;j=1，2……n（8） 　　它与测量数据的模型在形式上完全相同，表明式（2）建立的模型是合理的。 　　二、Prony算法在电力系统小干扰稳定分析中的应用 　　Prony算法本身意义明确，但直接按照其数学过程编制的程序对电力系统实际响应曲线或者仿真曲线进行辨识时，效果并不理想。本文结合电力系统的实际情况，对Prony算法的辨识过程进行研究，探讨算法在应用时需要注意的问题并提出解决方法。 　　1.Prony算法模型阶次的确定 　　Prony模型的阶数p的选择是对实际系