一种基于多变量状态估计的状态监测方法及系统与流程

本发明涉及旋转设备故障诊断领域，特别是涉及一种基于多变量状态估计的状态监测方法及系统。
背景技术：
：近年来，国内航空航天、轨道交通、电力等工业生产向系统化、数字化和自动化发展，机械系统逐渐复杂，对工艺参数要求指标和设备可靠性要求日益提高。旋转设备承载着长时间的运转工作，极易出现异常，进一步地影响生产效率和设备性能。因此对旋转机械进行监测，及时监测出旋转设备异常，实现旋转设备早期异常判断，避免工厂自动化生产线由于故障和维护导致的意外停机所带来的直接经济损失，实现近零计划外停机和无忧生产，提高生产效率。旋转设备发生异常情况时，温度升高，振动指标增大，频域产生故障频带。目前常用的旋转设备异常判断方法是单个特征指标进行阈值判断方法，根据iso2372标准阈值，对比计算振动信号的有效值，若超过阈值则判断设备出现异常情况。iso2372是一种通用的阈值判断标准，在不用应用场景下振动情况不同，采用通用的阈值判断引起的误报率和漏报率偏高。一般情况下，通用阈值偏大，振动超标超过通用阈值时，旋转设备已经出现不可逆损伤；并且现有方法仅能对当前时刻的设备状态判断，无法实现预测，现有状态监测系统中，设备早期异常报警机制的误报率与漏报率比较高。因此需要选择一种能够针对设备详情，构造可反映设备运行情况的特征，进一步地，基于多特征参数实现设备状态评估，能够实现设备异常预警，并降低旋转设备异常判断的误报率。技术实现要素：本发明是为了解决现有旋转设备异常判断方法中阈值固定，无法针对设备工况自适应调整阈值大小，并且特征参数根据通用经验选择，没有针对设备选择合适的特征参数；没有将多变量特征参数融合，不能进行综合评估设备状态的问题，提供一种基于多变量状态估计的状态监测方法，本发明能有效选择出可反映设备状态的关键参数，降低特征维度，减小状态估计时间，避免无效特征的干扰，能够提高状态估计的准确率，并且采用马氏距离方法将多维特征参数融合为一个估计值，并通过指数加权平均方法进行状态估计，可以降低因随机误差引发的波动，并能吸收瞬时突发的能力，时效性强。为了实现上述目的，本发明采用以下技术方案：一种基于多变量状态估计的状态监测方法，包括步骤：a)采集设备工况参数信号和振动信号，获得m组设备信号样本，对每组设备信号样本分别进行特征提取，获得每组设备振动信号样本的n1个特征和工况参数信号的n2个特征，获得m×n维的特征参数矩阵，n＝n1+n2；b)采用gbdt方法计算每个特征的重要度，根据重要度确定特征选择个数q，得到特征选择后的特征参数矩阵；c)采用马氏距离方法计算设备的状态估计值；d)采用指数加权平均方法建立设备异常判断机制，判断设备状态情况。本发明通过对设备振动信号进行时域、频域和时频域上的特征提取，能最大限度地提取反映设备运行状态的多维特征。通过gbdt方法，可提取出能真实反映设备运行状态的特征参数。通过马氏距离方法实现多维特征参数估计设备运行状态，实现设备异常预警。根据状态估计值建立异常报警机制，并通过指数加权平均方式，可降低异常预警误报率，判断设备状态情况准确率高。进一步的，步骤a)中，对每组设备信号样本分别进行特征提取，获得每组设备振动信号样本的n1个特征和工况参数信号的n2个特征，获得m×n维的特征参数矩阵，包括步骤：a1)提取各组设备振动信号样本的f1个时域特征；a2)提取各组设备振动信号样本的f2个频域特征；a3)对各组振动信号进行小波包变换，得到p个分量信号，分别计算p个分量信号的频带能量和能量熵；频带能量计算公式为能量熵的计算公式为其中i表示第i个分量信号，n表示分量信号数据点数，j表示分量信号数据点数序数，表示第k组设备振动信号样本的第i个分量信号的频带能量；a4)获得第k组设备振动信号的能量总和a5)获得每组设备振动信号样本的n1个特征，所述n1个特征包括时域特征、频域特征、频带能量、能量熵和能量总和，所述n2个工况参数信号特征包括电流、电压、转速等信息，获得特征参数矩阵，获得特征参数矩阵本发明通过对设备振动信号进行小波包变换，不仅能实现等频段的频带划分，而且能够提取各频带内能量，能最大限度地提取反映设备运行状态的多维特征。进一步的，步骤b)中采用gbdt(gradientboostingdescisiontree)方法计算每个特征的重要度，根据重要度确定特征选择个数q，得到特征选择后的特征参数矩阵，包括：b1)对n个特征分别计算在每颗树上的重要度，第j个特征在单颗树中的重要度的计算公式为其中，l为树的叶子节点数量，l-1为树的非叶子节点数量，vt是和节点t相关联的特征，是节点t分裂之后平方损失的减少值；b2)计算n个特征在每个树上的重要度的平均值，第j个特征在每个树上的重要度的平均值的计算公式为其中m为树的个数；b3)计算n个特征的重要度的平均值除以最大的重要度值，获得各个特征的重要度值，第j个特征的重要度值的计算公式为表示第j个特征的最大的重要度值；b4)将重要度值按从大到小顺序依次排列，获得排列后的重要度排列集合{i1,i2,…,in}，根据重要度确定特征选择个数q；b5)根据特征个数q，得到选择后的特征参数矩阵本发明采用gbdt方法计算每个特征的重要度，根据重要度确定特征选择个数q，对特征参数矩阵进行降维，获取前q个重要度最大的特征，再根据前q个重要度最大的特征重新构建特征参数矩阵，即得到选择后的特征参数矩阵v。通过采用gbdt重要度计算方法能有效选择出可反映设备状态的关键参数，通过参数选择方法一方面可以降低特征维度，减小状态估计时间；另一方面避免无效特征的干扰，可以提高状态估计的准确率。进一步的，步骤c)中采用马氏距离方法计算设备的状态估计值，包括步骤：c1)根据所述选择后的特征参数矩阵v计算各维矩阵的期望值μ＝[μ1,μ2,…,μq]t；c2)根据所述期望值μ计算特征参数矩阵的协方差矩阵∑；c3)根据所述特征参数矩阵的协方差矩阵∑计算马氏距离d，将所述马氏距离d作为设备的状态估计值。马氏距离(mahalanobisdistance)表示数据的协方差距离，是一种有效的计算两个未知样本集的相似度的方法，与欧氏距离不同的是它考虑到各种特性之间的联系，马氏距离不受量纲的影响，两点之间的马氏距离与原始数据的测量单位无关，由标准化数据和中心化数据(即原始数据与均值之差)计算出的二点之间的马氏距离相同。马氏距离还可以排除变量之间的相关性的干扰。对于一个均值为μ，协方差矩阵为∑-1的多变量向量x，其马氏距离为∑-1表示协方差矩阵对应的逆矩阵。进一步的，在步骤d)采用指数加权移动平均方法建立设备异常判断机制，判断设备状态情况，包括步骤：d1)获取一段时间内正常样本数据的马氏距离集合，计算所述马氏距离集合的均值μ1和方差σ1，根据3σ法则确定正常样本空间范围为[μ1-3σ，μ1+3σ]；d2)获取新进来的观测样本xt，计算观测样本xt的马氏距离dt；d3)计算加权移动平均后的观测样本估计值dt＝w0dt+w1dt-1+w3dt-2+…+wxdt-n，其中dt-n表示观测样本xt-n的马氏距离，wn代表第t-n时刻的马氏距离的权重，w0+w1+…wn＝1；d4)判断设备估计值dt是否在正常样本空间[μ1-3σ，μ1+3σ]内，若是，则表明设备正常；若否，则表明设备异常。本发明采用的指数加权平均估计方法能跟踪实际数据发生突变变化的能力，具备减少短期波动、保留长期发展趋势的能力。一种基于多变量状态估计的状态监测系统，包括数据获取模块、特征提取模块、特征选择模块、状态估计模块和异常判断模块；数据获取模块，用于采集设备工况参数信号和振动信号；特征提取模块，用于对所述数据获取模块采集到的设备工况参数信号和振动信号进行特征提取，所述特征包括时域特征、频域特征、频带能量、能量熵和能量总和；特征选择模块，用于采用gbdt方法计算每个特征的重要度对所述特征提取模块中提取到的设备信号特征计算重要度，并跟据重要度进行选取，获得特征选择后的特征参数矩阵；状态估计模块，用于采用马氏距离方法计算设备的状态估计值，获得设备的状态估计值；异常判断模块，用于采用指数加权移动平均方法建立设备异常判断机制，判断设备状态情况。本发明具有如下有益效果：本发明采用gbdt重要度计算方法，能有效选择出可反映设备状态的关键参数，通过参数选择方法一方面可以降低特征维度，减小状态估计时间；另一方面避免无效特征的干扰，可以提高状态估计的准确率。通过采用马氏距离方法，将多维特征参数融合为一个估计值，并通过指数加权平均方法进行状态估计，可以降低因随机误差引发的波动，并能吸收瞬时突发的能力，时效性较强。附图说明图1是实施例一流程示意图。图2是实施例一前100组样本马氏距离计算结果图。图3是实施例一后300组样本评估结果图。图4是实施例一采用37维特征的评估结果图。具体实施方式下面结合附图与具体实施方式对本发明做进一步的描述。实施例一，如图1所示，一种基于多变量状态估计的状态监测方法，以某型号电机为例，电机的转动速度为2500rpm，包括步骤：a)采集设备工况参数信号和振动信号，采样频率为25600hz，共m组振动信号样本，m＝400，对每组设备信号样本分别进行特征提取，获得每组设备振动信号样本的n1个特征和工况参数信号的n2个特征，获得m×n维的特征参数矩阵，n＝n1+n2，包括步骤：a1)提取各组设备振动信号样本的15个时域特征，15个时域特包含均值、绝对均值、方差、标准差、最大值、最小值、峰峰值、有效值、偏斜度、峰值指标、脉冲指标、裕度指标、波形指标、偏斜度指标、峭度和峭度指标；a2)提取各组设备振动信号样本的4个频域特征，4个频域特征包含谱有效值、谱方差、谱均值和谱中心；a3)对各组振动信号进行3层小波包变换，得到8个分量信号，分别计算8个分量信号的频带能量和能量熵；频带能量计算公式为能量熵的计算公式为其中i表示第i个分量信号，n表示分量信号数据点数，j表示分量信号数据点数序数，表示第k组设备振动信号样本的第i个分量信号的频带能量；a4)获得第k组设备振动信号的能量总和a5)获得每组设备振动信号样本的n1个特征和工况参数信号的n2个特征，n＝37，特征包括15个时域特征、4个频域特征、8个频带能量、8个能量熵和1个能量总和，获得400×37维的特征参数矩阵b)采用gbdt方法计算每个特征的重要度，根据重要度确定特征选择个数q，包括步骤：b1)对37个特征分别计算在每颗树上的重要度，第j个特征在单颗树中的重要度的计算公式为其中，l为树的叶子节点数量，l-1为树的非叶子节点数量，vt是和节点t相关联的特征，是节点t分裂之后平方损失的减少值；b2)计算37个特征在每个树上的重要度的平均值，第j个特征在每个树上的重要度的平均值的计算公式为其中m为树的个数；b3)计算37个特征的重要度的平均值除以最大的重要度值，获得各个特征的重要度值，第j个特征的重要度值的计算公式为表示第j个特征的最大的重要度值；b4)将重要度值按从大到小顺序依次排列，获得排列后的重要度排列集合{i1,i2,…,in}，根据重要度确定特征选择个数q。按照重要度排序，选择重要度超过50％的特征，即选择波形指标、谱重心、各频带能量熵值和全频带能量熵，共记11个特征，即q＝11。b5)根据特征个数得到选择后的特征参数矩阵c)采用马氏距离方法计算设备的状态估计值，包括步骤：c1)从特征参数矩阵v中选择前100组样本，计算各维矩阵的期望值μ＝[μ1，μ2，…，μq]t，期望值结果如表1所示。表1前100组样本的特征参数对应均值表特征参数均值特征参数均值特征参数均值谱重心6154.1波形指标1.3328h0.8953h10.1070h20.0985h30.1221h40.1172h50.0982h60.1023h70.1260h80.1241表1c2)根据所述期望值μ计算特征参数矩阵的协方差矩阵∑；c3)根据所述特征参数矩阵的协方差矩阵∑计算马氏距离d，将所述马氏距离d作为设备的状态估计值。图2是前100组信号的马氏距离结果。d)采用指数加权平均方法建立设备异常判断机制，判断设备状态情况，包括步骤：d1)获取前100正常样本数据的马氏距离集合，计算前100组数据的马氏距离的均值μ1和方差σ1，根据3σ法则确定正常样本空间范围为[μ1-3σ，μ1+3σ]；根据计算得到μ1＝7.315，σ1＝5.096，则正常样本空间范围为[0,22.6037]。d2)计算第101组至第400组观测样本的马氏距离dt，获取新进来的观测样本xt，计算观测样本xt的马氏距离dt；d3)计算加权移动平均后的观测样本估计值dt＝w0dt+w1dt-1+w2dt-2+…+wxdt-n，其中dt-n表示观测样本xt-n的马氏距离，wn代表第t-n时刻的马氏距离的权重，w0+w1+…wn＝1；d4)判断设备估计值dt是否在正常样本空间[0,22.6037]内，若是，则表明设备正常；若否，则表明设备异常。图3是计算300组样本的马氏距离估计值，在第160组样本时，估计值超出报警线，表明设备状态发生异常。图4是选择所有的37维特征进行状态估计的特征评估结果图，在第178组样本时设备估计值超出报警线，晚于本发明在第160组样本时就发现设备异常的报警时间。因此，本发明可提前发现设备异常情况，并实现提前报警。本发明通过对设备振动信号进行小波包变换，不仅能实现自适应的频带划分，而且能够提取各频带内能量，通过对设备振动信号进行时域、频域和时频域上的特征提取，能最大限度地提取反映设备运行状态的多维特征。通过gbdt方法能有效选择出可反映设备状态的关键参数，通过参数选择方法一方面可以降低特征维度，减小状态估计时间；另一方面避免无效特征的干扰，可以提高状态估计的准确率。通过马氏距离方法实现多维特征参数估计设备运行状态，实现设备异常预警。根据状态估计值建立异常报警机制，并通过指数加权平均方式，可降低异常预警误报率，对设备状态情况进行判断到的准确率高。本发明采用马氏距离方法，将多维特征参数融合为一个估计值，并通过指数加权平均方法进行状态估计，可以降低因随机误差引发的波动，并能吸收瞬时突发的能力，时效性较强。上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明保护范围以内。当前第1页12