一种基于振动相似度的故障检测方法及装置制造方法及图纸

本发明专利技术涉及振动测量技术领域，特别涉及一种基于振动相似度的故障检测方法及装置，装置包括数据获取模块、数据清洗模块、聚类建模模块、异常检测模块和相似度计算模块，检测方法为：S1，通过加速度传感器采集波形数据；S2，数据清洗；S3，对通过S2清洗后的数据进行聚类建模；S4，对聚类建模后的数据进行异常检测；S5，若S4中检测为非异常，进行相似度计算判定待测设备的故障劣化趋势。本发明专利技术通过自适应确定最优工况数，并根据波形特征自动确定类别归属，从而进行异常检测和相似度的计算。避免了正常波形参照基准选择困难的问题，同时异常检测利用概率统计方法自动确定检测阈值，提高了阈值设置的可靠性。设置的可靠性。设置的可靠性。

A fault detection method and device based on vibration similarity

【技术实现步骤摘要】
一种基于振动相似度的故障检测方法及装置


[0001]本专利技术涉及振动测量
，特别涉及一种基于振动相似度的故障检测方法及装置。

技术介绍

[0002]机械设备运行工况复杂,对设备进行故障诊断需要专业的诊断分析人员对采集的信号波形、频谱进行分析，通过对设备故障特征的了解，通过波形特征比对，确定设备可能存在的故障。
[0003]但在对设备结构信息不了解、故障特征库不完备情况下，如何利用正常数据及时发现异常情况，并触发预警检修成为避免设备进一步劣化的关键。并且机械设备采集的数据可能受到传感器、设备运行工况变化导致的数据异常情况，可能会造成较多的误报警。
[0004]现有技术大多基于时频域特征计算相似度值，并未考虑数据异常和设备运行工况变化可能对相似度计算造成的影响，如专利CN 112836574、CN 107024361，均为与预设特征量阈值进行比较，但实际上预设阈值的设定比较困难，且实际旋转机械设备的振动受到环境噪声和工况变化的影响，正常基准振动数据的选取比较困难。
[0005]为此，本申请设计了一种基于振动相似度的故障检测方法及装置，以解决上述问题。

技术实现思路

[0006]本专利技术为了弥补现有技术中机械设备运行环境复杂、技术人员对设备故障特征不了解的不足，提供了一种基于振动相似度的故障检测方法及装置，通过数据清洗自动过滤无效数据对检测可能造成的影响，利用聚类算法，自动确定所属工况，并基于特征值实现振动的异常检测和相似度的计算。
[0007]一种基于振动相似度的故障检测装置，其特征在于：包括数据获取模块、数据清洗模块、聚类建模模块、异常检测模块和相似度计算模块，所述数据获取模块为模拟压电加速度传感器或数字式加速度传感器，数据获取模块设置在待测设备的传递轴或齿轮箱位置。
[0008]基于上述的基于振动相似度的故障检测装置，其检测方法为：S1，通过加速度传感器采集波形数据，采集波形的采样频率表示为fs，采样点数表示为N；S2，数据清洗，对S1中采集到的波形数据进行计算，判定数据是否存在异常；S3，对通过S2清洗后的数据进行聚类建模；S4，对聚类建模后的数据进行异常检测；S5，若S4中检测为非异常，进行相似度计算判定待测设备的故障劣化趋势。
[0009]进一步地，为了更好的实现本专利技术，所述S2具体为，S21，计算波形的过零率指标，过零率为波形序列中任意两相邻点异号，并且其绝对值差值大于一设定值thr，所述thr大于待测设备静止状态时信号噪声的幅度；
S22，计算S21中信号通过零点的次数，每通过一次，计数一次，总计count次，则过零率pass_zeros_per=count/N；S23，记录S21过程中过零点的所有位置到一记录数组location_i[],并取出数组中第一个元素的位置pass_zeros_first_index；S24，对记录数组location[]求取其前向差分数组的最大值，记录为pass_zeros_maxdiff；S25，给定过零率阈值thr2、起始过零率阈值thr3、过零的最大差分值thr4,Ifpass_zeros_per<thr2||pass_zeros_first_index>th3||pass_zeros_maxdiff>th4,则判定数据为异常无效。
[0010]进一步地，为了更好的实现本专利技术，所述S3具体为，S31，获取待测设备在不同运行工况下，在相同采样频率、采样点数采集的振动波形数据若干条；S32，对采集的任一波形数据进行傅里叶变换，可得到信号的频域幅值表示Y；S33，根据采集数据的频率分辨率，对频域幅值序列Y进行平均分段累加处理，得到信号频域在不同频段的幅值和数组,并对幅值和数组进行归一化处理，得到特征数组z[]；根据信号的最大分析频率、分辨率确定最适宜分段数=10~20；S34，预设聚类个数2~10,利用kmeans算法对特征数组z[]进行聚类操作,并记录每一聚类个数对应的轮廓系数；所述轮廓系数计算为：计算任一个波形特征值数组z[]与同类的其他样本特征值数组的平均距离a;同时计算其与不同类的所有样本特征值数组的平均距离的最小值b;则此样本对应的轮廓系数:，计算所有波形特征值数组对应的轮廓系数的平均值，得到全局轮廓系数S；S35，找到轮廓系数最大值对应的聚类个数，为最优聚类类别个数m，并保存最优聚类类别数为m的每一类所对应的聚类中心和聚类方差。
[0011]进一步地，为了更好的实现本专利技术，所述S4具体为，S41，对实时新采集到的波形，计算其对应的特征数组z1[]，并利用S35建模得到的聚类中心，计算z1[]到m个聚类中心的距离，并将其归属到距离最小的簇类中心所对应的类别m1中；所述距离的定义为欧式距离：为与对应聚类中心特征值差的平方开根号；S42，取出m1类别建模对应的聚类中心mean_cluster_center、聚类方差std_cluster_center，确定异常检测的阈值为：thres_category_max=mean_cluster_center+3*std_cluster_center;S43，求取z1[]特征数组超出thres_category_max对应特征阈值的个数，为防止误报警，设置超出个数如大于3，则判定为异常，机械设备存在故障隐患。
[0012]进一步地，为了更好的实现本专利技术，所述S5具体为，若S4检测为非异常，则对特征数组z1[]对应的波形A和m1类别中一组波形B,求取波形A和波形B对应的fft幅值序列：A_fft,B_fft和fft对应频率序列F；相似度定义为：Corr{A,B}=abs((F*A_fft)
‑
F*B_fft)/(F*B_fft)；
为了避免由于频谱噪声导致累加干扰主频的统计，对相似计算的A_fft和B_fft 进行如下预处理，以A_fft为例：首先，求取A_fft序列的均值和标准差，得到过滤阈值threshold1 = 均值+3倍标准差；然后，对A_fft序列中小于threshold1的幅值执行置0操作；最后，对处理后的A_fft幅值序列进行上述相似度计算；定义相似度检测阈值，对波形A和波形B计算得到的结果进行判定，若相似度计算值大于检测阈值，则判定设备振动正常，波形频谱结构相似度较高，未发生明显的特征改变；若相似度计算值小于检测阈值，证明波形频谱结构相似度较差，机械设备振动的特征发生改变，存在故障劣化的趋势。
[0013]进一步地，为了更好的实现本专利技术，所述S34中的kmeans算法可替换为密度聚类DBSCAN或层次聚类。
[0014]进一步地，为了更好的实现本专利技术，所述S3中的聚类建模可采用相似度先行计算得到相应类别的参照，具体如下：Sa，首先利用振动加速度传感器，采集获取机械振动设备月或季度的历史数据若干；Sb，然后随机从历史数据中选择一条数据, 对剩余所有数据依次利用相似度计算模块，计算得到对应的相似度数组corr_list1[];Sc，再通过预定义的相似度阈值，找出相似度数组c本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于振动相似度的故障检测装置，其特征在于：包括数据获取模块、数据清洗模块、聚类建模模块、异常检测模块和相似度计算模块，所述数据获取模块为模拟压电加速度传感器或数字式加速度传感器，数据获取模块设置在待测设备的传递轴或齿轮箱位置。2.一种基于权利要求1所述的基于振动相似度的故障检测装置的检测方法，其特征在于，包括以下步骤：S1，通过加速度传感器采集波形数据，采集波形的采样频率表示为fs，采样点数表示为N；S2，数据清洗，对S1中采集到的波形数据进行计算，判定数据是否存在异常；S3，对通过S2清洗后的数据进行聚类建模；S4，对聚类建模后的数据进行异常检测；S5，若S4中检测为非异常，进行相似度计算判定待测设备的故障劣化趋势。3.根据权利要求2所述的基于振动相似度的故障检测方法，其特征在于：所述S2具体为，S21，计算波形的过零率指标，过零率为波形序列中任意两相邻点异号，并且其绝对值差值大于一设定值thr，所述thr大于待测设备静止状态时信号噪声的幅度；S22，计算S21中信号通过零点的次数，每通过一次，计数一次，总计count次，则过零率pass_zeros_per=count/N；S23，记录S21过程中过零点的所有位置到一记录数组location_i[], 并取出数组中第一个元素的位置pass_zeros_first_index；S24，对记录数组location[]求取其前向差分数组的最大值，记录为pass_zeros_maxdiff；S25，给定过零率阈值thr2、起始过零率阈值thr3、过零的最大差分值thr4,If pass_zeros_per < thr2 || pass_zeros_first_index >th3 || pass_zeros_maxdiff > th4, 则判定数据为异常无效。4.根据权利要求2所述的基于振动相似度的故障检测方法，其特征在于：所述S3具体为，S31，获取待测设备在不同运行工况下，在相同采样频率、采样点数采集的振动波形数据若干条；S32，对采集的任一波形数据进行傅里叶变换，可得到信号的频域幅值表示Y；S33，根据采集数据的频率分辨率，对频域幅值序列Y进行平均分段累加处理，得到信号频域在不同频段的幅值和数组, 并对幅值和数组进行归一化处理，得到特征数组z[]；根据信号的最大分析频率、分辨率确定最适宜分段数=10~20；S34，预设聚类个数2~10, 利用kmeans算法对特征数组z[]进行聚类操作,并记录每一聚类个数对应的轮廓系数；所述轮廓系数计算为：计算任一个波形特征值数组z[]与同类的其他样本特征值数组的平均距离a; 同时计算其与不同类的所有样本特征值数组的平均距离的最小值b;则此样本对应的轮廓系数:，计算所有波形特征值数组对应的轮廓系数的平均值，得到全局轮廓系数S；
S35，找到轮廓系数最大值对应的聚类个数，为最优聚类类别个数m，并保存最优聚类类别数为m的每一类所对应的聚类中心和聚类方差。5.根据权利要求2所述的基于振动相似度的故障检测方法，其特征在于：所...

【专利技术属性】
技术研发人员：谭晓林，张盛英，李小良，李玉奎，
申请(专利权)人：济南嘉宏科技有限责任公司，
类型：发明
国别省市：