专利名称：基于mm和arma算法的次同步振荡模态辨识方法
技术领域：
本发明涉及电力系统的安全稳定运行，特别涉及一种基于MM和ARMA算法的次同步振荡模态辨识方法。
背景技术：
随着电网规模的扩大化和结构的复杂化，电力系统中的次同步振荡问题亟待解决。快速准确地辨识次同步振荡模态参数，对于次同步振荡的监测、预警以及控制措施的制定都起着至关重要的作用。随着广域测量系统(WAMS)相关技术的快速发展，基于实测数据提取振荡模态的方法由于其无需详细的系统模型和大规模特征值计算而逐渐得到广泛应用。目前，基于实测数据的次同步振荡模态辨识的典型方法有:普罗尼(Prony)分析法、小波分析法、快速傅里叶(FFT)法、希尔伯特-黄变化法(HHT)等。传统的Prony分析法抗噪能力差；小波分析法难区分信号中的相近频率，不利于扭振模态参数的提取；FFT法无法提取振荡的瞬时频率和衰减因子；HHT法对信号的采样率要求较高。上述方法都难以满足复杂系统强噪声背景下的次同步振荡模态辨识的要求。数学形态学(mathematical morphology, MM)是基于积分几何和随机集论发展而来，具有计算快速、简便、去噪和重构信号能力强的特点，被广泛应用于图像处理、形状分析、模式识别等领域。董超等[1]采用数学形态学设计出一种多结构元素的并行复合形态滤波器，可有效滤除多种噪声，保留更多的有用信息。对消噪后的信号采用基于总体最小二乘法-旋转不变技术的信号参数估计(TLS — ESPRIT)算法进行辨识，从而获取低频振荡各个模式参数。自回归滑动平均模型(ARMA)是建立在线性模型基础上，是以参数化模型处理动态随机数据的一种实用方法，也是系统辨识和预测的经典方法。ARMA模型以白噪声为输入，解决了系统识别中输入未知时不可识别的问题，拓宽了系统识别的应用领域。吴超等[2]基于类噪声信号，采用自回归滑动平均(ARMA)法进行低频振荡模式辨识，从而实现电网正常运行状态下的动态稳定性预警。文中涉及如下参考文献:[I]董超等.基于数学形态学滤波技术和TLS-ESPRIT算法的低频振荡模式辨识研究[J].电力系统保护与控制，2012，40 (3):114-118, 123.
[2]吴超等.计及模型定阶的低频振荡模式类噪声信号辨识[J].电力系统自动化，2009，33(21):1-6.
发明内容
针对现有技术存在的问题，本发明将数学形态滤波技术与自回归滑动平均模型结合，提出了一种基于数学形态学(MM)和自回归滑动平均模型(ARMA)的次同步振荡模态辨识方法。为解决上述技术问题，本发明采用如下技术方案:
一种基于丽和ARMA算法的次同步振荡模态辨识方法，包括步骤:步骤1、采用数学形态滤波器对次同步振荡信号进行消噪；步骤2、对消噪后的次同步振荡信号构建自回归滑动平均模型，并对信号进行预处理；步骤3，对构建的自回归滑动平均模型进行定阶，并对定阶后的自回归滑动平均模型进行模型参数估计；步骤4，基于步骤3所获取的自回归滑动平均模型参数估计值进行次同步振荡信号模态辨识。步骤I中的数学形态滤波器为形态开-闭滤波器和形态闭-开滤波器的组合，具
体为
权利要求
1.一种基于MM和ARMA算法的次同步振荡模态辨识方法，其特征在于，包括步骤: 步骤1、采用数学形态滤波器对次同步振荡信号进行消噪； 步骤2、对消噪后的次同步振荡信号构建自回归滑动平均模型，并对信号进行预处理；步骤3，对构建的自回归滑动平均模型进行定阶，并对定阶后的自回归滑动平均模型进行模型参数估计； 步骤4，基于步骤3所获取的ARMA模型参数估计值进行次同步振荡信号模态辨识。
2.如权利要求1所述的基于MM和ARMA算法的次同步振荡模态辨识方法，其特征在于: 步骤I中所述的数学形态滤波器为形态开-闭滤波器和形态闭-开滤波器的组合。
3.如权利要求2所述的基于MM和ARMA算法的次同步振荡模态辨识方法，其特征在于: 所述的数学形态滤波器为
4.如权利要求1所述的基于MM和ARMA算法的次同步振荡模态辨识方法，其特征在于: 步骤I中所述的数学形态滤波器的结构元素为半圆结构元素和三角结构元素。
5.如权利要求1所述的基于MM和ARMA算法的次同步振荡模态辨识方法，其特征在于: 采用相似度指标S评估步骤I的滤波效果，所述的相似度指标
6.如权利要求1所述的基于MM和ARMA算法的次同步振荡模态辨识方法，其特征在于: 步骤3中采用BIC准则对构建的自回归滑动平均模型进行定阶。
7.如权利要求1所述的基于MM和ARMA算法的次同步振荡模态辨识方法，其特征在于: 步骤3中采用加权递推最小二乘法对定阶后的自回归滑动平均模型进行模型参数估计。
8.如权利要求7所述的基于MM和ARMA算法的次同步振荡模态辨识方法，其特征在于: 所述的采用加权递推最小二乘法对定阶后的自回归滑动平均模型进行模型参数估计，具体为: 采用加权递推最小二乘法对定阶后的自回归滑动平均模型进行模型参数估计，根据所得模型参数估计值计算加权准则函数，判断加权准则函数是否满足控制精度要求，若满足，则执行步骤4 ;否则，基于当前模型参数估计值对自回归滑动平均模型继续进行模型参数估计。
9.如权利要求1所述的基于MM和ARMA算法的次同步振荡模态辨识方法，其特征在于: 步骤4进一步包括子步骤: 4-1基于自回归滑动平均模型参数估计值构建次同步振荡离散信号的特征方程，并获取特征方程的共轭特征根； 4-2基于共轭特征根分析自回归滑动平均模型参数，并对分析所得自回归滑动平均模型参数进行z域到s域的变换，从而获得次同步振荡模型参数。
全文摘要
本发明公开了一种基于MM和ARMA算法的次同步振荡模态辨识方法，该方法实现了噪声干扰下的次同步振荡模态准确辨识。本发明方法包括步骤采用数学形态滤波器对次同步振荡信号进行消噪处理，保留信号的主要特征信息；对消噪后的信号建立ARMA模型，并经数据预处理、模型定阶、模型参数估计后获得次同步振荡模态参数。本发明方法能快速、准确地辨识出次同步振荡的模态参数，具有抗噪能力强、辨识精度高等优点，在基于实测数据的次同步振荡分析、监测、预警及阻尼控制器设计等方面具有较好的应用前景。
文档编号G06F17/50GK103207931SQ201310076469
公开日2013年7月17日 申请日期2013年3月11日 优先权日2013年3月11日
发明者廖清芬, 董飞飞, 唐飞, 查晓明, 陈恩泽, 魏大千, 涂炼, 朱振山 申请人:武汉大学