一种基于电流波形相似度的发电机故障诊断系统及方法与流程

1.本发明是一种应用于风电机组电气/机械故障的诊断方法，具体涉及一种基于电流波形相似度的发电机故障诊断系统及方法。


背景技术：

2.风电机组运行工况较为复杂，常常面临高温、高湿、盐雾腐蚀、雷电等更加复杂的环境条件，风电功率波动剧烈。其故障表现形式与运行状态密切相关，实现故障辨识具有一定的难度，因此，必须掌握有效的风电机组部件故障特征提取方法，才可以提高机组的故障诊断和状态预测精度，预先判断机组健康衰退趋势，合理调整运行并安排维修，对提高风电机组运行安全性和可靠性，降低运行与维护费用。传统电流波形分析采用快速傅里叶变换、小波变换、小波基函数等手段频域获取电流信息。然而频域分析仍然存在频谱泄漏、分辨率受限、频谱模糊、小波基函数选择对分析结果影响大等问题。


技术实现要素：

3.本发明的目的在于提供一种基于电流波形相似度的发电机故障诊断系统及方法，解决了现有技术中存在的上述不足。
4.为了达到上述目的，本发明采用的技术方案是：
5.本发明提供的一种基于电流波形相似度的发电机故障诊断方法，包括以下步骤：
6.步骤1，获取发电机定子三相电流的波形数据，并对三相电流的波形数据进行归一化处理，得到归一化处理后的波形数据；
7.步骤2，从步骤1中得到的归一化处理后的波形数据中分别获取设定时间内的a相电流的波形数据和c相电流的波形数据，分别记为时间序列a和时间序列c；
8.步骤3，计算时间序列c和时间序列a之间的相似度；
9.步骤4，根据步骤3中得到的相似度判断发电机状态。
10.优选地，步骤3中，计算时间序列c和时间序列a之间的相似度，具体方法是：
11.s31，分别计算得到的矩阵d规整路径的最小规整路径距离、时间序列a的最小规整路径距离，以及时间序列a与噪声波形之间的最小规整路径距离；
12.s32，将得到的矩阵d规整路径的最小规整路径距离、时间序列a与自身之间的最小规整路径距离，以及时间序列a与噪声波形之间的最小规整路径距离进行标准化处理，得到标准化后的距离；
13.s33，根据标准化后的距离计算时间序列c和时间序列a之间的相似度。
14.优选地，s31中，计算得到的矩阵d规整路径的最小规整路径距离，具体方法是：
15.s311，利用时间序列a和时间序列c构建矩阵d；
16.s312，对得到的矩阵d进行规整路径处理，得到多条规整路径；
17.s313，根据得到的多条规整路径计算得到的矩阵d规整路径的最小规整路径距离。
18.优选地，s311中，利用时间序列a和时间序列c构建矩阵d，具体方法是：
19.利用动态时间归整方法对齐时间序列a与时间序列c，得到矩阵d。
20.优选地，s31中，计算时间序列a与自身之间的最小规整路径距离，具体方法是：
21.利用动态时间归整方法对齐时间序列a与时间序列a，得到矩阵d1；
22.对得到的矩阵d1进行规整路径处理，得到多条规整路径；
23.根据得到的多条规整路径计算得到时间序列a与自身之间的最小规整路径距离。
24.优选地，s31中，计算时间序列a与噪声波形之间的最小规整路径距离，具体方法是：
25.利用动态时间归整方法对齐时间序列a与噪声波形，得到矩阵d2；
26.对得到的矩阵d2进行规整路径处理，得到多条规整路径；
27.根据得到的多条规整路径计算得到时间序列a与噪声波形之间的最小规整路径距离。
28.优选地，步骤4中，根据步骤3中得到的相似度判断发电机状态，具体方法是：
29.当0.980＜相似度≤1.000时，则判断发电机状态为正常；
30.当0.900＜相似度≤0.980时，则判断发电机状态为正常；
31.当0.500＜相似度≤0.900时，则判断发电机状态为异常；
32.当相似度≥0.5时，则判断发电机状态为严重异常。
33.一种基于电流波形相似度的发电机故障诊断系统该系统，能够运行所述的一种基于电流波形相似度的发电机故障诊断方法，包括采集单元、计算单元和判断单元，具体地：
34.采集单元，用于获取发电机定子三相电流的波形数据，并对三相电流的波形数据进行归一化处理，得到归一化处理后的波形数据；
35.从得到的归一化处理后的波形数据中分别获取设定时间内的a相电流的波形数据和c 相电流的波形数据，分别记为时间序列a和时间序列c；
36.计算单元，用于计算时间序列c和时间序列a之间的相似度；
37.判断单元，用于根据得到的相似度判断发电机状态。
38.与现有技术相比，本发明的有益效果是：
39.本发明提供的一种基于电流波形相似度的发电机故障诊断方法，基于原始时域波形，通过对比实际采集的电流波形与典型波形的相似度分析设备状态，该方法能够利用相似度预测故障的发生；本发明针对风电机组电气监测参量变化复杂快速、特征量提取不易等特点；利用采集电流波形相似度实时反映风电机组设备状态，分析设备故障概率。
40.进一步的，利用时间序列a与自身之间的最小规整路径距离，时间序列a与噪声波形之间的最小规整路径距离计算出的相似度作为相似度的上限与下限，如此则规定了任一时间序列与序列a的相似度可能区间，使得计算获得的相似度更加准确，更能反映实际状态。
附图说明
41.图1为本发明的流程示意图。
42.图2为正常状态下1相似度展示图。
43.图3为正常状态2下相似度展示图。
44.图4为异常状态下相似度展示图。
具体实施方式
45.下面结合附图，对本发明进一步详细说明。
46.本发明提供的一种基于电流波形相似度的发电机故障诊断方法，针对风电机组电气监测参量变化复杂快速、特征量提取不易等特点；利用采集电流波形相似度实时反映风电机组设备状态，分析设备故障概率；具体包括以下步骤：
47.步骤1，选择电流传感器收集的发电机定子三相电流波形数据，风机运行过程中电定子电流幅值变化较为明显，因此必须对采集数据进行归一化处理，避免幅值的影响，通过下式对采集三相电流每一相的幅值进行归一化处理：
[0048][0049]
其中，i
norm
为归一化后的幅值，i为电流数据中实际采集值，i
min
为一相电流采集数据中幅值最小的值，i
max
为数据中最大值。
[0050]
步骤2，选择一段时间内经过步骤1归一化后采集的a相电流(至少包含10个完整周期) 做为对比基准，记为时间序列a，则该序列可以表示为：
[0051]
a＝{a1,a2,a3,a4,a5,
…
,a
n
}
[0052]
同理，b相电流时间序列与c相电流时间序列(至少包含10个完整周期)可表示为：
[0053]
b＝{b1,b2,b3,b4,b5,
…
,b
m
}
[0054]
c＝{c1,c2,c3,c4,c5,
…
,c
j
}
[0055]
n、m、j分别为时间序列a，b，c元素个数，即序列长度。
[0056]
步骤3，利用动态时间归整方法对齐时间序列a与时间序列c(序列a与序列b对齐方法相同)，具体方法为：
[0057]
采用欧式距离计算序列a中每个元素与序列c中每个元素的距离，根据计算结果构造一个n
×
m的矩阵d，其中，矩阵d中元素d(p，q)表示a
p
和c
q
两个点的距离d(a
p
，c
q
)心
[0058]
d(a
p
，c
q
)＝(a
p
‑
c
q
)2[0059]
其中，a
p
为时间序列a中的元素；c
q
为时间序列c中的元素。
[0060]
步骤4，对步骤3中获得的矩阵d进行warping path规整路径，用w来表示，表达式为：
[0061]
w＝{w1，w2，w3，w4，...，w
k
}，max(m，n)≤k＜m+n
‑1[0062][0063]
其中，w的第k个元素定义为w
k
＝(p，q)
k
，定义了序列a和c的映射；dtw(a，c)表示由矩阵d的(0，0)点开始至矩阵d的(n，m)点终止的路径；k表示路径长度。
[0064]
步骤5，根据步骤3与步骤4的计算结果，计算最小规整路径距离，首先定义一个累加距离cd(cumulative distances)，从矩阵d的(0，0)点开始沿着步骤4中得到的多条路径匹配这两个序列a和c，每到一个点，之前所有的点计算的距离都会累加。到达终点(n，m) 后，这个累积距离就是最后总的距离cd，提取根据多条路径计算而得距离cd中的最小值cd
ag
即为序列a与序列c的最小规整路径距离。累加距离cd计算公式如下：
[0065]
cd(p，q)＝d(a
p
，c
q
)+min{cd(p
‑
1，q)，cd((p
‑
1，q
‑
1))，cd(p，q
‑
1)}
[0066]
步骤6，重复步骤3至步骤5，计算a相时间序列与自身的最小规整路径距离cd1，这
个距离定义为cd
min
；特别的，选取噪声波形计算与a相标准波形最小规整路径距离cd0，这个距离定义为cd
max
。则cd
ac
一定位于[cd
min
，cd
max
]区间内。
[0067]
步骤7，将步骤5，6计算所得cd
min
、cd
max
与cd
ac
标准化在[0，1]区间内，cd
acnor
为标准化后的距离，具体计算方式如下：
[0068][0069]
则有cd
maxnor
＝1，cd
minnor
＝0；计算c序列与标准a相序列相似度，计算方式如下：
[0070][0071]
其中，fi
ac
表示c序列与标准a相序列相似度。
[0072]
步骤8，根据步骤7中计算获得相似度fi判断电机状态，按照以下标准：
[0073][0074]
若相似度(fi)数值在0.9
‑
0.5范围内，则提出预警，提示工作人员关注电机状态；若相似度(fi)数值在0.5以下，则立马提出报警，通知工作人员及时处理。
[0075]
附图2
‑
4分别展示了利用波形相似度开展分析的不同结果，其中，图2中，c相与a 相相似度长期保持在0.999以上，表明状态正常，无需关注；图3中c相与a相相似度长期保持在0.993以上，表明状态正常，但略有下降趋势；图4中c相与a相相似度在0.9
‑
0.5 之间波动，表明两相之间相似度低，已经存在异常，需要关注。