一种次同步振荡传播关键影响因素识别方法及系统与流程

1.本发明涉及电力系统及其自动化领域，特别是涉及一种次同步振荡传播关键影响因素识别方法及系统。


背景技术：

2.风电场在电力系统电源侧广泛投建，其与传统设备完全不同的动态响应特性会重塑系统的动态行为，引发次同步振荡问题。同时由于电力资源分布与需求的逆向分布，负荷中心与可再生能源丰富地区相距甚远，为解决供电用电需求，建设了大规模的交直流混联电网。随着大规模风电外送经由交直流线路输送功率，各区域电网之间联系加强，导致振荡能量会经由联络线在各区域电网之间传播，导致振荡传播的影响范围进一步扩大。而直流输电在新型电力系统要求下将逐步扩大应用范围，其电力电子装置控制系统会与新能源场站相互作用从而引发新的振荡问题，新能源场站设备与交直流电网之间的相互作用较为复杂，其传播影响因素对于电网安全可靠运行具有重大影响。
3.目前在次同步振荡的研究中，关于次同步振荡在电网中传输规律的问题尚有不足，研究的场景较为狭隘，多数是进行局部传播规律的研究，针对次同步振荡的全局性问题的研究尚缺少相关理论支撑，故考虑多因素次同步振荡传播的关键影响因素的识别具有重要的意义，有利于及时切除振荡源、维持新型电力系统的安全稳定。


技术实现要素：

4.本发明的目的是提供一种次同步振荡传播关键影响因素识别方法及系统，以识别次同步振荡传播关键影响因素，在次同步振荡发生时更有效的识别、抑制次同步振荡源。
5.为实现上述目的，本发明提供了如下方案：
6.一种次同步振荡传播关键影响因素识别方法，包括：
7.以交直流线路参数中的任意一种参数为识别对象，多次调整所述识别对象的值，并在每次调整后获得所述识别对象对应的一组交直流线路有功功率信号；一组交直流线路有功功率信号包括交直流线路上不同位置处的有功功率信号；
8.对所述识别对象对应的每个有功功率信号进行ceemdan分解，获得每个有功功率信号的多个imf分量；
9.根据每个有功功率信号的多个imf分量，采用prony算法分析每个有功功率信号对应的系统阻尼比；
10.根据每个位置处的多个系统阻尼比和所述识别对象多次调整的值，利用相关性评价组合模型，确定用于度量每个位置与所述识别对象相关性的评价指标；
11.从用于度量每个位置与所有识别对象相关性的评价指标中选取处于评价指标阈值范围内的评价指标，并将选取的评价指标对应的识别对象作为每个位置发生次同步振荡时的关键影响因素。
12.可选的，所述对所述识别对象对应的每个有功功率信号进行ceemdan 分解，获得
每个有功功率信号的多个imf分量，具体包括：
13.向有功功率信号中多次加入高斯白噪声信号；
14.在每次加入高斯白噪声信号后进行一次emd分解，获得多个模态分量；
15.将多个模态分量的平均值作为本征模态分量；
16.计算有功功率信号去除第一个本征模态分量后的残差，并判断所述残差是否为单调函数，获得判断结果；
17.若所述判断结果表示否，则在残差中多次加入经emd分解后的高斯白噪声信号，并将所述残差替换所述有功功率信号，同时返回步骤“在每次加入高斯白噪声信号后进行一次emd分解，获得多个模态分量”；
18.若所述判断结果表示是，则输出所有的本征模态分量，作为每个有功功率信号的多个imf分量。
19.可选的，所述对所述识别对象对应的每个有功功率信号进行ceemdan 分解，获得每个有功功率信号的多个imf分量，之后还包括：
20.计算每个imf分量与每个imf分量对应的有功功率信号的相关系数；
21.将相关系数小于相关系数阈值的imf分量确定为目标去噪imf分量；
22.采用小波阈值去噪方法对每一个目标去噪imf分量进行去噪，获得去噪后的imf分量；
23.将每个有功功率信号的所有去噪后的imf分量和每个有功功率信号的未去噪imf分量合并重构，获得每个去噪后的有功功率信号。
24.可选的，所述根据每个有功功率信号的多个imf分量，采用prony算法分析每个有功功率信号对应的系统阻尼比，具体包括：
25.根据每个有功功率信号的多个imf分量，采用prony算法确定次同步振荡频率下的特征值；
26.根据次同步振荡频率下的特征值，利用公式确定每个有功功率信号对应的系统阻尼比；其中，λi表示次同步振荡频率下的特征值，ξ表示系统阻尼比，ωi表示振荡角频率。
27.可选的，所述根据每个位置处的多个系统阻尼比和所述识别对象多次调整的值，利用相关性评价组合模型，确定用于度量每个位置与所述识别对象相关性的评价指标，具体包括：
28.根据每个位置处的多个系统阻尼比和所述识别对象多次调整的值，利用相关性评价组合模型，分别计算皮尔逊相关系数、斯皮尔曼相关系数和互信息量；
29.对每个位置对应的皮尔逊相关系数、斯皮尔曼相关系数和互信息量进行加权求和，获得用于度量每个位置与所述识别对象相关性的评价指标。
30.一种次同步振荡传播关键影响因素识别系统，包括：
31.信号获取模块，用于以交直流线路参数中的任意一种参数为识别对象，多次调整所述识别对象的值，并在每次调整后获得所述识别对象对应的一组交直流线路有功功率信号；一组交直流线路有功功率信号包括交直流线路上不同位置处的有功功率信号；
32.分解模块，用于对所述识别对象对应的每个有功功率信号进行ceemdan 分解，获
得每个有功功率信号的多个imf分量；
33.分析模块，用于根据每个有功功率信号的多个imf分量，采用prony算法分析每个有功功率信号对应的系统阻尼比；
34.评价指标确定模块，用于根据每个位置处的多个系统阻尼比和所述识别对象多次调整的值，利用相关性评价组合模型，确定用于度量每个位置与所述识别对象相关性的评价指标；
35.识别模块，用于从用于度量每个位置与所有识别对象相关性的评价指标中选取处于评价指标阈值范围内的评价指标，并将选取的评价指标对应的识别对象作为每个位置发生次同步振荡时的关键影响因素。
36.可选的，所述分解模块，具体包括：
37.噪声加入子模块，用于向有功功率信号中多次加入高斯白噪声信号；
38.emd分解子模块，用于在每次加入高斯白噪声信号后进行一次emd分解，获得多个模态分量；
39.平均子模块，用于将多个模态分量的平均值作为本征模态分量；
40.判断子模块，用于计算有功功率信号去除第一个本征模态分量后的残差，并判断所述残差是否为单调函数，获得判断结果；
41.循环子模块，用于若所述判断结果表示否，则在残差中多次加入经emd 分解后的高斯白噪声信号，并将所述残差替换所述有功功率信号，同时返回步骤“在每次加入高斯白噪声信号后进行一次emd分解，获得多个模态分量”；
42.输出子模块，用于若所述判断结果表示是，则输出所有的本征模态分量，作为每个有功功率信号的多个imf分量。
43.可选的，还包括：
44.相关系数计算模块，用于计算每个imf分量与每个imf分量对应的有功功率信号的相关系数；
45.分量筛选模块，用于将相关系数小于相关系数阈值的imf分量确定为目标去噪imf分量；
46.去噪模块，用于采用小波阈值去噪方法对每一个目标去噪imf分量进行去噪，获得去噪后的imf分量；
47.重构模块，用于将每个有功功率信号的所有去噪后的imf分量和每个有功功率信号的未去噪imf分量合并重构，获得每个去噪后的有功功率信号。
48.可选的，所述分析模块，具体包括：
49.特征值确定子模块，用于根据每个有功功率信号的多个imf分量，采用 prony算法确定次同步振荡频率下的特征值；
50.系统阻尼比确定子模块，用于根据次同步振荡频率下的特征值，利用公式确定每个有功功率信号对应的系统阻尼比；其中，λi表示次同步振荡频率下的特征值，ξ表示系统阻尼比，ωi表示振荡角频率。
51.可选的，所述评价指标确定模块，具体包括：
52.评价子模块，用于根据每个位置处的多个系统阻尼比和所述识别对象多次调整的
值，利用相关性评价组合模型，分别计算皮尔逊相关系数、斯皮尔曼相关系数和互信息量；
53.加权子模块，用于对每个位置对应的皮尔逊相关系数、斯皮尔曼相关系数和互信息量进行加权求和，获得用于度量每个位置与所述识别对象相关性的评价指标。
54.根据本发明提供的具体实施例，本发明公开了以下技术效果：
55.本发明公开一种次同步振荡传播关键影响因素识别方法及系统，基于 ceemdan-小波阈值去噪技术和prony分析计算不同位置处参数对应的系统阻尼比，根据相关性评价组合模型判断交直流线路参数与不同位置系统阻尼比的线性相关性，形成基于相关性评价组合模型的次同步振荡传播关键影响因素识别，从而能在次同步振荡发生时更有效的识别、抑制次同步振荡源。
附图说明
56.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
57.图1为本发明实施例提供的次同步振荡传播关键影响因素识别方法的原理图；
58.图2为本发明实施例提供的两区域四机系统等效电路图；
59.图3为本发明实施例提供的vsc-hvdc定直流电压控制策略示意图；
60.图4为本发明实施例提供的vsc-hvdc定有功功率控制策略示意图；
61.图5为本发明实施例提供的ceemdan分解后的imf分量示意图；
62.图6为本发明实施例提供的ceemdan分解后的各imf分量与原信号的相关系数示意图；
63.图7为本发明实施例提供的原信号去噪前后对比示意图；图7中的(a) 为含噪信号示意图，图7中的(b)为去噪后信号示意图；
64.图8为本发明实施例提供的系统阻尼比的相关性组合模型的评价指标总图；
65.图9为本发明实施例提供的dfig出口系统阻尼比的相关性组合模型的评价指标图；
66.图10为本发明实施例提供的交流线路系统阻尼比的相关性组合模型的评价指标图；
67.图11为本发明实施例提供的直流线路系统阻尼比的相关性组合模型的评价指标图。
具体实施方式
68.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
69.本发明的目的是提供一种次同步振荡传播关键影响因素识别方法及系统，以识别次同步振荡传播关键影响因素，在次同步振荡发生时更有效的识别、抑制次同步振荡源。
70.为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
71.本发明针对目前大规模风电场并网次同步振荡传播规律研究的不足，提出了一种次同步振荡传播关键影响因素识别方法，如图1所示，包括以下步骤：
72.步骤s1，以交直流线路参数中的任意一种参数为识别对象，多次调整所述识别对象的值，并在每次调整后获得所述识别对象对应的一组交直流线路有功功率信号。一组交直流线路有功功率信号包括交直流线路上不同位置处的有功功率信号。
73.交直流线路参数包括交流线路潮流、交流线路阻抗特性以及直流线路的电缆参数与直流线路中vsc控制参数。其中，交流线路阻抗特性具体指次同步频率交流线路阻抗。直流线路的电缆参数具体指直流线路电缆长度。支流线路中vsc控制参数包括定电压电流内环比例、定电压电流内环积分、定电压电流外环比例、定电压电流外环积分、定功率电流内环比例、定功率电流内环积分、定功率功率外环比例和定功率功率外环积分。
74.举例来说，调整识别对象的值的方式：控制其他参数不变，在仿真模型数值稳定的前提下，识别对象取100个，仿真得到100个样本所对应的有功功率。
75.步骤s2，对识别对象对应的每个有功功率信号进行ceemdan分解，获得每个有功功率信号的多个imf分量。
76.ceemdan分解方法如下：
77.(1)vj(t)为满足标准正态分布的高斯白噪声信号，j＝1,2,
…
,n为加入白噪声次数，将高斯白噪声vj(t)加入到待分解信号y(t)得到新信号εj表示第j次加入高斯白噪声的强度，表示第j次加入噪声后经过emd分解的第一个分量，rj表示第j次加入噪声后经过emd 分解的残差值。其中，ε：高斯白噪声的强度，第j次加入噪声后经过 emd分解的第一个分量，rj：第j次加入噪声后经过emd分解的残差值。
78.(2)对产生的n个模态分量进行总体平均就得到ceemdan分解的第1 个本征模态分量：
[0079][0080]
(3)计算去除第一个模态分量后的残差：
[0081][0082]
(4)在r1(t)中加入经emd处理后的n次噪声，定义ej(
·
)表示信号经emd 分解后的第j个imf分量，得到n个新信号r1(t) e1(ε1vj(t))，以新信号 r1(t) e1(ε1vj(t))为载体进行emd分解，得到第二阶模态分量c2(t)，由此可以得到ceemdan分解的第2个本征模态分量：
[0083][0084]
(5)计算去除第二个模态分量后的残差：
[0085][0086]
(6)重复(3)至(5)步骤，直到获得的残差为单调函数，第k阶的imf 表达式如下所示：
[0087][0088]
算法结束后，此时得到的本征模态imf分量数量为k，则原始信号y(t)被分解为：
[0089][0090]
可见，原信号分解为k个imf分量和一个残差，后续主要是针对k个imf 分量进行小波阈值去噪。
[0091]
下面利用相关系数选取目标去噪分量，利用小波阈值去噪方法对目标去噪分量进行去噪方法如下：
[0092]
(1)利用相关系数筛从经过ceemdan分解后的k个imf分量中，选出h个高频目标去噪imf分量优选方式为：将相关系数小于0.6的imf 分量选取为高频目标去噪imf分量。
[0093]
(2)对每一个目标去噪分量进行小波变换，得到一组小波分解系数 w
j,k
；
[0094]
(3)选择合适的阈值处理小波分解系数，得到小波估计系数使得尽可能的小；
[0095]
(4)利用估计的小波系数进行小波重构，得到估计信号ck'(t)，该信号即为去噪后的信号；
[0096]
(5)将h个去噪后的分量ck'(t)、k-h个未去噪分量以及残差重新合并，即可得到去噪后的有效信号，如下式所示：
[0097][0098]
步骤s3，根据每个有功功率信号的多个imf分量，采用prony算法分析每个有功功率信号对应的系统阻尼比。
[0099]
prony分析不同位置系统阻尼比方法如下：
[0100]
prony算法旨在利用复指数函数的线性组合拟合等间距采样数据的数学模型，一般把采用p个指数项的线性组合模型成为p阶模型。
[0101]
设x(0)，x(1)，
……
，x(n-1)为n个采样数据，则：
[0102][0103]
式中：为拟合后的prony模型，bi和zi为复数，定义：
[0104][0105]
由欧拉公式可得：
[0106]
[0107]
即只需求解bi和zi，可用p个具有任意幅值ai、频率fi、相位θi和衰减因子αi的余弦分量来拟合采样数据，prony模型与采样数据的误差平方为：
[0108][0109]
当该误差平方数值最小时，即可得到幅值ai、频率fi、相位θi和衰减因子αi的解，注意到拟合的prony模型是常系数差分方程的齐次解，可得常系数线性差分方程为：
[0110][0111]
对参数ai进行最小二乘估计可得方程为：
[0112][0113]
令可得prony方法的法方程式：
[0114][0115]
其中解该法方程，即可得到系数a1,
···ap
的估计值，之后可得到特征多项式
[0116]
1 a1z-1

···
a
p
z-p
＝0
[0117]
的根zi,i＝1,2,
···
,p，于是指数模型式可以化简为未知参数和bi的线性方程：
[0118][0119]
该方程的最小二乘解为：
[0120][0121]
故可得幅值ai、频率fi、相位θi和衰减因子αi分别为：
[0122]ai
＝|bi|
[0123][0124]
αi＝ln|zi|/δt
1.0之间认为参数极强相关，0.6-0.8强相关，0.4-0.6中等程度相关，0.2-0.4弱相关，0.0-0.2极弱相关或无相关。
[0147]
(2)斯皮尔曼相关系数，其计算方法如下：
[0148]
斯皮尔曼等级相关是一种非参数统计量，其值与两组相关变量的具体值无关，而仅仅与其值之间的大小关系有关。斯皮尔曼等级相关依据两列成对等级的各对等级数之差进行计算，所以又称为“等级差数法”。当变量在至少是有序的尺度上测量时，它是合适的相关分析方法。斯皮尔曼相关系数的数学表达如下：
[0149][0150]
式中，xi为自变量，yi为因变量，为自变量平均值，为因变量平均值。
[0151]
(3)互信息量，其计算方法如下：
[0152]
互信息源于信息论，用于衡量2个变量之间的关联程度。互信息量的计算公式为：
[0153][0154]
式中，i(x,y)为x与y之间的互信息；p(x,y)为x与y共同出现的频率；p(x)、 p(y)分别为x、y单独发生的频率。若i(x,y)0，数值越大表明x与y的关联程度越强；若i(x,y)≈0，则x、y的关联程度比较弱，若i(x,y)＜0，则x与y互补，不存在关联关系。
[0155]
步骤s5，从用于度量每个位置与所有识别对象相关性的评价指标中选取处于评价指标阈值范围内的评价指标，并将选取的评价指标对应的识别对象作为每个位置发生次同步振荡时的关键影响因素。
[0156]
在相关性组合模型的度量中，认为评价指标在0.8-1.0之间认为参数极强相关，0.6-0.8强相关，0.4-0.6中等程度相关，0.2-0.4弱相关，0.0-0.2极弱相关或无相关。
[0157]
本发明公开了属于电力系统安全稳定领域的一种次同步振荡传播关键影响因素识别方法，目前快速识别和切除振荡源仍缺乏有效方法，对风电并网系统次同步振荡关键影响因素进行定量分析研究，有助于规划风电并网系统，同时在系统运行过程中及时调整控制参数，解决工程实际问题。其过程为：对风机出口、交流线路和直流线路三个不同位置的有功功率进行特征提取，对各位置的有功功率进行ceemdan分解，用相关系数选取目标去噪分量；用小波阈值去噪方法对目标去噪分量进行去噪，并与其他未去噪分量重构；建立皮尔逊相关系数、斯皮尔曼相关系数以及次同步振荡特征分量间的互信息量的相关性评价组合模型，综合各类型相关系数计算交直流线路不同特性的相关性评价指标，对比得出次同步振荡传播关键影响因素。该方法对及时切除振荡源、维持新型电力系统的安全稳定有着重要的意义。
[0158]
下面两区域四机系统为例，图2为两区域四机系统等效电路图。直流输电系统控制方法为vsc-hvdc，vsc两侧分别采用定直流电压控制和定有功功率控制，其中定直流电压控制分别电压外环和电流内环，定有功功率控制为功率外环和电流内环，可以通过改变控制器的pi参数调整vsc-hvdc的输电特性，vsc两种控制策略示意图如图3、4所示；交流输电系统为有串联电容补偿的交流线路。两区域四机系统由汽轮发电机和双馈风机组成，在两区域四机的一端接入双馈风机，通过调整双馈风机参数使其发生次同步振荡。
[0159]
1.调整交直流线路参数获取次同步振荡功率数据，经ceemdan分解后得到一系列imf分量，如图5所示；通过相关系数筛选出目标去噪分量，相关系数如图6所示；对目标去噪分量进行小波阈值去噪后重构imf分量，获得去噪后分量，去噪前后对比图如图7所示。
[0160]
2.使用去噪后的次同步振荡功率数据，基于prony分析拟合不同位置的次同步振荡有功功率波形，利用不同指数项的线性组合模型模拟有功功率波形，即用p个具有任意幅值ai、频率fi、相位θi和衰减因子αi的余弦分量来拟合采样数据，通过特征方程得到不同频率对应的相关参数，通过对特征矩阵进行计算求得系统振荡的特征值λi＝σi±
j2πfi。
[0161]
由系统振荡的特征值和阻尼比的关系可以计算不同位置处监测到的系统阻尼比ξ。以交流潮流参数为例其系统阻尼比分析结果如下表所示：
[0162]
表1系统阻尼比分析结果
[0163]
交流潮流参数dfig出口交流线路直流线路184mw-2.3-39-16220mw-4.6-38-10475mw-3-24-15
[0164]
3.计算交直流线路不同参数与对应的系统阻尼比之间的相关性组合模型的评价指标。交直流参数与各位置次同步功率阻尼比的相关性组合模型的评价指标如图8-11所示，其中图8为不同位置处不同参数的相关性组合模型的评价指标总图，交直流参数与各位置系统阻尼比相关性组合模型的评价指标如图9-11所示。图8中
①
为交流线路潮流，
②
为次同步频率交流线路阻抗，
③
为直流线路电缆长度，
④
为定电压电流内环比例，
⑤
为定电压电流内环积分，
⑥
为定电压电流外环比例，
⑦
为定电压电流外环积分，
⑧
为定功率电流内环比例，
⑨
为定功率电流内环积分，
⑩
为定功率功率外环比例，为定功率功率外环积分。
[0165]
在相关性组合模型的度量中，认为评价指标在0.8-1.0之间认为参数极强相关，0.6-0.8强相关，0.4-0.6中等程度相关，0.2-0.4弱相关，0.0-0.2极弱相关或无相关。由评价指标计算结果可以得到：对于dfig出口次同步有功功率阻尼比来说，参数极强相关的为vsc的定直流电压控制中的电流内环积分系数；对于交流线路次同步功率阻尼比来说，参数极强相关的为交流线路潮流；对于直流线路次同步功率阻尼比来说，参数极强相关的为次同步频率下交流线路阻抗特性。
[0166]
所以对于该两区域四机系统来说，不同位置发生次同步振荡时可以利用其对应的关键影响因素调整系统阻尼比，从而控制次同步振荡在交直流电网中的传播。
[0167]
本发明的有益效果是：设置交直流线路参数，基于ceemdan-小波阈值去噪技术和prony分析计算不同位置处参数对应的系统阻尼比，根据相关性评价组合模型判断交直流线路参数与不同位置系统阻尼比的线性相关性，形成基于相关性评价组合模型的次同步振荡传播关键影响因素识别，从而能在次同步振荡发生时更有效的对电力系统进行操作以识别、抑制次同步振荡源。
[0168]
本发明实施例还提供了一种次同步振荡传播关键影响因素识别系统，包括：
[0169]
信号获取模块，用于以交直流线路参数中的任意一种参数为识别对象，多次调整所述识别对象的值，并在每次调整后获得所述识别对象对应的一组交直流线路有功功率信号；一组交直流线路有功功率信号包括交直流线路上不同位置处的有功功率信号；
[0170]
分解模块，用于对所述识别对象对应的每个有功功率信号进行ceemdan 分解，获得每个有功功率信号的多个imf分量；
[0171]
分析模块，用于根据每个有功功率信号的多个imf分量，采用prony算法分析每个有功功率信号对应的系统阻尼比；
[0172]
评价指标确定模块，用于根据每个位置处的多个系统阻尼比和所述识别对象多次调整的值，利用相关性评价组合模型，确定用于度量每个位置与所述识别对象相关性的评价指标；
[0173]
识别模块，用于从用于度量每个位置与所有识别对象相关性的评价指标中选取处于评价指标阈值范围内的评价指标，并将选取的评价指标对应的识别对象作为每个位置发生次同步振荡时的关键影响因素。
[0174]
分解模块，具体包括：
[0175]
噪声加入子模块，用于向有功功率信号中多次加入高斯白噪声信号；
[0176]
emd分解子模块，用于在每次加入高斯白噪声信号后进行一次emd分解，获得多个模态分量；
[0177]
平均子模块，用于将多个模态分量的平均值作为本征模态分量；
[0178]
判断子模块，用于计算有功功率信号去除第一个本征模态分量后的残差，并判断所述残差是否为单调函数，获得判断结果；
[0179]
循环子模块，用于若所述判断结果表示否，则在残差中多次加入经emd 分解后的高斯白噪声信号，并将所述残差替换所述有功功率信号，同时返回步骤“在每次加入高斯白噪声信号后进行一次emd分解，获得多个模态分量”；
[0180]
输出子模块，用于若所述判断结果表示是，则输出所有的本征模态分量，作为每个有功功率信号的多个imf分量。
[0181]
次同步振荡传播关键影响因素识别系统还包括：
[0182]
相关系数计算模块，用于计算每个imf分量与每个imf分量对应的有功功率信号的相关系数；
[0183]
分量筛选模块，用于将相关系数小于相关系数阈值的imf分量确定为目标去噪imf分量；
[0184]
去噪模块，用于采用小波阈值去噪方法对每一个目标去噪imf分量进行去噪，获得去噪后的imf分量；
[0185]
重构模块，用于将每个有功功率信号的所有去噪后的imf分量和每个有功功率信号的未去噪imf分量合并重构，获得每个去噪后的有功功率信号。
[0186]
分析模块，具体包括：
[0187]
特征值确定子模块，用于根据每个有功功率信号的多个imf分量，采用 prony算法确定次同步振荡频率下的特征值；
[0188]
系统阻尼比确定子模块，用于根据次同步振荡频率下的特征值，利用公式确定每个有功功率信号对应的系统阻尼比；其中，λi表示次同步振荡频率下的特征值，ξ表示系统阻尼比，ωi表示振荡角频率。
[0189]
评价指标确定模块，具体包括：
[0190]
评价子模块，用于根据每个位置处的多个系统阻尼比和所述识别对象多次调整的值，利用相关性评价组合模型，分别计算皮尔逊相关系数、斯皮尔曼相关系数和互信息量；
[0191]
加权子模块，用于对每个位置对应的皮尔逊相关系数、斯皮尔曼相关系数和互信息量进行加权求和，获得用于度量每个位置与所述识别对象相关性的评价指标。
[0192]
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
[0193]
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。