CN111122191A - 一种基于ewma控制的设备安康报警阈值设定方法 - Google Patents

Abstract

本发明涉及智能控制领域，旨在提供一种基于EWMA控制的设备安康报警阈值设定方法。该方法是利用数据采集模块实时采集运转设备的振动加速度信号，将其存储于数据存储模块中的历史数据库；数据处理模块提取历史数据库中的振动加速度值，去除异常值后，根据指数加权移动平均值控制图划分四个报警区，实现设备报警阈值设定。本发明所述方法可以根据具体设备的实际运行情况确定报警阈值，具有较高的时效性和较强的适应性；对故障的辨识度高，能够及时发现设备故障的早期征兆，有效延长设备运行寿命，降低维修费用，提高设备可靠性，保证设备全生命周期的安全、高效；能够存储有效数据，使得监测数据可信度高，有利于科学决策和进一步分析。

Description

一种基于EWMA控制的设备安康报警阈值设定方法

技术领域

本发明涉及智能控制领域，具体地说，涉及到一种基于EWMA(指数加权移动平均，Exponentially Weighted Moving Average)控制的设备安康报警阈值设定方法。

背景技术

近年来，随着科学技术的进步和生产效率的提高，运转设备进入了机、电结合的新阶段，不断向高速和高精度方向发展，结构日趋复杂。在工作过程中运转设备会经历正常-退化-最终失效三个阶段，其可用性、安全性和可靠性在不同程度上影响着企业的生产效率。因此，对运转设备的运行状态进行实时监测和故障报警，对于企业降低维护成本和生产损失、提高生产效益具有重要意义。

设备状态报警技术是设备状态监测和故障诊断的核心技术之一，设备运行状态等级的划分会影响到报警方法的确定及报警精度。现有的设备安康报警技术均是依据振动标准进行状态等级划分，存在划分标准的不确定、报警阈值与设备实际运行状态无关等问题。

如能根据具体设备的实际运行情况确定报警阈值，使得报警系统有较高的时效性和较强的适应性，能够及时发现设备故障的早期征兆，有效延长设备运行寿命，降低维修费用，同时提高设备可靠性，保证设备全生命周期的安全、高效。因此，提供一种新的设备安康报警阈值设定方法是十分有必要的。

发明内容

本发明要解决的技术问题是，克服现有技术中的不足，提供一种基于EWMA控制的设备安康报警阈值设定方法，以解决运转设备的故障智能报警问题。

为解决技术问题，本发明的解决方案为：

提供一种基于EWMA控制的设备安康报警阈值设定方法，是利用数据采集模块实时采集运转设备的振动加速度信号，将其存储于数据存储模块中的历史数据库；数据处理模块提取历史数据库中的振动加速度值，去除异常值后，根据指数加权移动平均值控制图划分四个报警区，实现设备报警阈值设定；

其中，数据处理模块的处理过程具体包括：

(1)数据处理服务器从历史数据库中提取数据X₁，X₂，X₃……X₁₀₂₄，进行快速傅里叶变换以获取振动加速度频谱图，然后根据频谱得到指定频率分量的幅值U_i(i为采样大小，1≤i≤1024)；

(2)根据运转设备的运行机理设定关机值a和异常值b，去除小于关机值a或者大于异常值b的幅值U_i；

(3)利用最近的n个振动加速度测量值的算术平均值构建评估过程水平的控制图；根据处于受控状态的频率分量的幅值U_i确定分组大小n和分组间隔m，根据设备运行时的受控数据和失控数据，通过指定可接受的最小稳态平均运行链长(ARL)确定控制限参数L值和λ值，以实现所需的稳态平均运行链长；最后根据公式(10)-(12)得到控制上限UCL、中心线CL、控制下限LCL的值，确定指数加权移动平均控制图；

CL＝μ (11)

各符号的定义：假设随机变量为X～N(μ，σ²)，μ为过程均值；λ为平滑系数，常数，其取值范围为0＜λ≤1；

(4)设定A、B、C、D四个报警区分别代表“安全运行”、“可长期运行”“不可长期运行”和“停机检修”；四个报警区分别对应：

“(-∞，LCL]U[UCL，+∞)”这四个范围。

本发明中，所述快速傅里叶变换是指离散傅里叶变换快速算法，按下述计算过程将振动加速度由时域转换为频域：

DFT公式：

其中，X(k)为经过傅里叶变换的频域上加速度序列值；X(n)为采集的时域上振动加速度序列值，

为旋转因子，由公式(14)计算得到；k为该采样时刻点的序数，N为采样总点数，N＝1024；

X_(n)＝X_1(n)+X_2(n) (15)

X_1(n)、X_2(n)分别为偶数序列和奇数序列，其长度均为N/2；

所以

由于

所以

k＝0，……，N-1 (19)

其中，X_1(k)、X_2(k)分别为X_1(n)、X_2(n)的N/2点DFT；

所以，

依此类推，得到每个采样点的频域幅值，进而得到振动加速度频谱图。

本发明中，该方法基于下述系统实现：系统数据采集模块、数据存储模块、数据处理模块和数据输出模块；其中，数据采集模块用于实时采集设备的振动数据，将其存储在数据存储模块；数据存储模块用于存储记载设备各种运行状态的历史数据库；数据处理模块用于提取和处理数据存储模块中的数据，根据指数加权移动平均值控制确定报警阈值；数据输出模块用于显示数据处理结果，实现报警信号发布；

数据处理模块将振动加速度信号经过处理后与预设的报警阈值进行比较后，若结果显示振动烈度值在预设的报警阈值以内，则将数据发送至数据存储模块进行存储，否则进行故障报警并根据振动烈度所在区域更新历史数据库。

本发明中，所述运转设备是水泵、风机、空压机、电机、电梯曳引机或电梯轿厢中的任意一种。

本发明中，所述数据存储模块是基于本地或云端的数据存储服务器，所述数据处理模块是数据处理服务器；所述数据输出模块是手机、电脑或智能显示屏中的任意一种。

本发明中，所述数据采集模块由振动检测传感器和无线网关组成，其中振动检测传感器用于采集运转设备的振动信号并通过无线网关上传至数据处理服务器。

本发明中，所述振动检测传感器是有线振动传感器、LoRa振动传感器、NB-IoT振动传感器或WiFi振动传感器中的任意一种；所述无线网关是LoRa无线网关或Wifi无线网关中的任意一种。

本发明还进一步提供了基于所述方法的运转设备安康报警方法，包括：由数据处理模块对实时采集到的振动加速度信号进行分析，如其频率分量的幅值U_i落在A或B区域内，认为设备正常运行；若频率分量的幅值U_i落在C或D区域内，则认为设备处于故障状态，此时由数据输出模块进行报警，并由数据处理模块并将落在C区域的频率分量的幅值U_i对应的原始数据通过数据处理模块更新至历史数据库；对于落在D区域的原始数据则直接丢弃。

本发明还进一步提供了基于所述方法的运转设备安康报警方法，包括：由数据处理模块根据预置规则对经过处理的频率分量的幅值U_i进行比较和判断，确认并输出故障类型诊断结果。

与现有技术相比，本发明的有益效果是：

1、本发明所述方法提供了一种全新的设备安康报警阈值设定方法，可以根据具体设备的实际运行情况确定报警阈值，具有较高的时效性和较强的适应性；

2、本发明对故障的辨识度高，能够及时发现设备故障的早期征兆，有效延长设备运行寿命，降低维修费用，提高设备可靠性，保证设备全生命周期的安全、高效；

3、本发明能够存储有效数据，使得监测数据可信度高，有利于科学决策和进一步分析。

附图说明

图1是本发明整体设计图。

图2是本发明安康报警阈值设定方法流程图。

图3是本发明指数加权移动平均值控制图模型构建示意图。

图4是本发明故障类型诊断示意图。

具体实施方式

下面结合附图与具体实施例对本发明做进一步的说明。

首先需要说明的是，本发明的实现会涉及信号检测技术，系计算机技术在工业生产领域的应用。在本发明的具体应用过程中，可能涉及到软件功能模块的应用。申请人认为，如在仔细阅读申请文件、准确理解本发明的实现原理和发明目的以后，在结合现有公知技术的情况下，本领域技术人员完全可以运用其掌握的软件编程技能实现本发明。凡本发明申请文件提及的均属此范畴，申请人不再一一列举。

移动平均控制图是利用最近n个观测值(其中现行观测值已取代n+1个观测值中的最初观测值)的算术平均值，是评估过程水平的控制图。移动平均控制图可以改善在生产过程中发现小的变动的能力。指数加权移动平均(Exponentially Weighted MovingAverage，即EWMA)控制图是移动平均控制图的一种变形图，考虑了当前样本的信息比先前样本提供的信息更重要的移动平均控制图就称为指数加权移动平均控制图，它不仅能够及时检测出生产过程均值所发生的较小波动(小于2σ)，而且对当前过程的突变性变异具有一定的检出效果。

如图1所示，一种基于EWMA控制的设备安康报警阈值设定方法，基于下述系统实现：包括数据采集模块、数据存储模块、数据处理模块、数据输出模块；数据采集模块用于实时采集设备的振动数据，将其存储在数据存储模块；数据存储模块用于存储记载设备各种运行状态的历史数据库；数据处理模块用于提取和处理数据存储模块中的数据，根据指数加权移动平均值控制确定报警阈值；数据输出模块用于显示数据处理结果，实现报警信号发布；数据处理模块将振动加速度信号经过处理后与预设的报警阈值进行比较后，若结果显示振动烈度值在预设的报警阈值以内，则将数据发送至数据存储模块进行存储，否则进行故障报警并根据振动烈度所在区域更新历史数据库。

基于该系统的设备安康报警阈值设定方法，是利用数据采集模块实时采集运转设备的振动加速度信号，将其存储于数据存储模块中的历史数据库；数据处理模块提取历史数据库中的振动加速度值，去除异常值后，根据指数加权移动平均值控制图划分四个报警区，实现设备报警阈值设定；

其中，数据处理模块的处理过程具体包括：

(1)数据处理服务器从历史数据库中提取数据X₁，X₂，X₃……X₁₀₂₄，进行快速傅里叶变换以获取振动加速度频谱图，然后根据频谱得到指定频率分量的幅值U_i(i为采样大小，1≤i≤1024)；

(2)根据运转设备的运行机理设定关机值a和异常值b，去除小于关机值a或者大于异常值b的幅值U_i；

(3)利用最近的n个振动加速度测量值的算术平均值构建评估过程水平的控制图；根据处于受控状态的频率分量的幅值U_i确定分组大小n和分组间隔m，根据设备运行时的受控数据和失控数据，通过指定可接受的最小稳态平均运行链长(ARL)确定控制限参数L值和λ值，以实现所需的稳态平均运行链长；最后根据公式(10)-(12)得到控制上限UCL、中心线CL、控制下限LCL的值，确定指数加权移动平均控制图；

CL＝μ (11)

各符号的定义：假设随机变量为X～N(μ，σ²)，μ为过程均值；λ为平滑系数，常数，其取值范围为0＜λ≤1；

(4)设定A、B、C、D四个报警区分别代表“安全运行”、“可长期运行”“不可长期运行”和“停机检修”；四个报警区分别对应：

“(-∞，LCL]∪[UCL，+∞)”这四个范围。

基于上述方法的运转设备安康报警方法包括：由数据处理模块对实时采集到的振动加速度信号进行分析，如其频率分量的幅值U_i落在A或B区域内，认为设备正常运行；若频率分量的幅值U_i落在C或D区域内，则认为设备处于故障状态，此时由数据输出模块进行报警，并根据振动烈度所在区域将当前结果更新至历史数据库。根据当前结果更新至历史数据库，是指若频率分量的幅值U_i在A、B、C任一区域，则将该频率分量的幅值U_i所对应的振动加速度值保存至数据存储服务器。

进一步地，由数据处理模块根据预置规则对经过处理的频率分量的幅值U_i进行比较和判断，确认并输出故障类型诊断结果。故障类型可设定为轴承故障、基座松动故障、轴承松动故障、不平衡动故障、不对中故障、转子弯曲故障、动静件摩擦故障和早期潜在故障。

下面结合具体实施例子，对上述内容进一步阐释：

为得到频域数据，涉及的公式如下：

快速傅里叶变换(Fast Fourier Transform，FFT)，是离散傅里叶变换(DiscreteFourier Transform，DFT)的快速算法，它是根据离散傅里叶变换的奇、偶、虚、实等特性，对离散傅里叶变换的算法进行改进获得的，可将振动加速度由时域转换为频域。

快速傅里叶变换公式推算过程如下：

DFT公式：

其中，X(k)为经过傅里叶变换的频域上加速度序列值；X(n)为采集的时域上振动加速度序列值，

为旋转因子，可由公式(14)计算得到；k为该采样时刻点的序数，N为采样总点数，本发明中N＝1024。

X_(n)＝X_1(n)+X_2(n) (15)

X_1(n)、X_2(n)分别为偶数序列和奇数序列，其长度均为N/2。

所以

由于

所以

(k＝0，……，N-1) (19)

其中，X_1(k)、X_2(k)分别为X_1(n)、X_2(n)的N/2点DFT。

所以，

依此类推，得到每个采样点的频域幅值，进而得到振动加速度频谱图。

指数加权移动平均控制图模型建立中涉及到的公式如下：

假设随机变量频率分量幅值U_i为X～N(μ，σ²)，μ为过程均值，μ₀为过程均值目标值，μ的取值根据受控状态与否而变化，当过程处于受控状态，μ＝μ₀；当过程处于失控状态，μ＝μ₀+δσ，δ为过程均值漂移。

过程均值：

标准差：

指数加权移动平均控制图的定义公式为：

Z_i＝λx_i+(1-λ)Z_i-1 (24)

其中，λ为平滑系数，常数，其取值范围为0＜λ≤1，Z_i是EWMA统计量，即所有之前样本均值的加权平均值。初始值Z₀＝μ₀。

因为假设观测值xi是独立的随机变量，而其方差为σ²，所以z_i的方差为

因此，指数加权移动平均控制图的控制限为：

CL＝μ (27)

当i逐渐增大时，(1-λ)²ⁱ将很快收敛到0，。因此当i逐渐增大时，UCL和LCL将稳定到下面两个值：

CL＝μ (30)

其中，UCL为控制上限，CL为中心线，LCL为控制下限；n为采样分组数量，L为控制限系数，λ为平滑系数。通常，λ取值在0.05≤λ≤0.25范围在实践中有较好的表现，实践中通常取λ＝0.05，λ＝0.10或者λ＝0.20；而对于L，L＝3的使用效果不错，特别是当λ较大的时候。另外，当λ≤0.1时，2.6≤L≤2.8能有效地使控制限变窄。最优的控制限参数L值和λ值通过指定可接受的较小的稳态平均运行链长(ARL)确定。

如图2所示，基于EWMA的设备安康报警阈值设定方法，包括：数据处理服务器从数据库存储服务器中提取数据X₁，X₂，X₃......X₁₀₂₄，进行快速傅里叶变换，得到振动加速度信号的频谱，根据频谱得到指定频率分量的幅值U_i大小，并根据相应设备的运行机理，设定设备的关机值a和异常数据b，去除小于关机值a或者异常大的数据b值。根据指数加权移动平均值控制图，分别设定

“(-∞，LCL]∪[UCL，+∞)”为A、B、C、D四个区；并判断当前实时数据U_i处于哪个区；若指定频率分量的幅值U_i在A、B区域，认为设备正常运行，若指定频率分量的幅值U_i在C、D区域则进行故障类型诊断和报警，并根据当前结果更新历史数据库。(对于落在C区域的频率分量的幅值U_i对应的原始数据，更新至历史数据库；对于落在D区域的原始数据则直接丢弃。)

如图3所示，指数加权移动平均值控制图模型构建包括：

根据处于受控状态的频率分量的幅值U_i，确定分组大小n和分组间隔m，然后根据设备运行时的受控数据和失控数据，通过指定可接受的较小的稳态平均运行链长(ARL)确定控制限参数L值和λ值，来实现所需的ARL；最后根据公式得到UL、UCL和LCL的值，确定指数加权移动平均控制图。

如图4所示，故障类型诊断过程如下：

分别对0.5倍频、1倍频、1.5倍频、2倍频、2.5倍频、3倍频、4倍频和高倍频的频率分量的幅值U_i，根据设定的8个规则进行故障诊断，输出的诊断结果主要有轴承故障、基座松动故障、轴承松动故障、不平衡动故障、不对中故障、转子弯曲故障、动静件摩擦故障和早期潜在故障。

对设备状态报警阈值进行设定是一个技术难点，本发明提出基于指数加权移动平均值控制图的设备安康报警阈值设定方法，具有较高的时效性和较强的适应性，为防止关键设备的异常提供预警，有效延长设备运行寿命，降低维修费用，提高设备的可靠性，保证设备全生命周期的安全、高效，在机械设备报警处理系统中具有良好的应用前景。

上述内容虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，任何本领域技术人员，在不脱离本发明的精神和范围内所做出的应用或改动，均应属于本发明的保护范围。

Claims (9)

1.一种基于EWMA控制的设备安康报警阈值设定方法，其特征在于，是利用数据采集模块实时采集运转设备的振动加速度信号，将其存储于数据存储模块中的历史数据库；数据处理模块提取历史数据库中的振动加速度值，去除异常值后，根据指数加权移动平均值控制图划分四个报警区，实现设备报警阈值设定；

其中，数据处理模块的处理过程具体包括：

(1)数据处理服务器从历史数据库中提取数据X₁，X₂，X₃……X₁₀₂₄，进行快速傅里叶变换以获取振动加速度频谱图，然后根据频谱得到指定频率分量的幅值U_i，i为采样大小，1≤i≤1024；

(2)根据运转设备的运行机理设定关机值a和异常值b，去除小于关机值a或者大于异常值b的幅值U_i；

(3)利用最近的n个振动加速度测量值的算术平均值构建评估过程水平的控制图；根据处于受控状态的频率分量的幅值U_i确定分组大小n和分组间隔m，根据设备运行时的受控数据和失控数据，通过指定可接受的最小稳态平均运行链长(ARL)确定控制限参数L值和λ值，以实现所需的稳态平均运行链长；最后根据公式(10)-(12)得到控制上限UCL、中心线CL、控制下限LCL的值，确定指数加权移动平均控制图；

CL＝μ (11)

各符号的定义：假设随机变量为X～N(μ，σ²)，μ为过程均值；λ为平滑系数，常数，其取值范围为0＜λ≤1：

(4)设定A、B、C、D四个报警区分别代表“安全运行”、“可长期运行”“不可长期运行”和“停机检修”；四个报警区分别对应：

“(-∞，LCL]U[UCL，+∞)”这四个范围。

2.根据权利要求1所述的方法，其特征在于，所述快速傅里叶变换是指离散傅里叶变换的快速算法，按下述计算过程将振动加速度由时域转换为频域：

离散傅里叶变换公式：

其中，X(k)为经过傅里叶变换的频域上加速度序列值；X(n)为采集的时域上振动加速度序列值，

为旋转因子，由公式(14)计算得到；k为该采样时刻点的序数，N为采样总点数，N＝1024；

X_(n)＝X_1(n)+X_2(n) (15)

X_1(n)、X_2(n)分别为偶数序列和奇数序列，其长度均为N/2；

所以

由于

所以

其中，X_1(k)、X_2(k)分别为X_1(n)、X_2(n)的N/2点DFT；

所以，

依此类推，得到每个采样点的频域幅值，进而得到振动加速度频谱图。

3.根据权利要求1所述的方法，其特征在于，该方法基于下述系统实现：包括数据采集模块、数据存储模块、数据处理模块和数据输出模块；其中，

数据采集模块用于实时采集设备的振动数据，将其存储在数据存储模块；

数据存储模块用于存储记载设备各种运行状态的历史数据库；

数据处理模块用于提取和处理数据存储模块中的数据，根据指数加权移动平均值控制确定报警阈值；

数据输出模块用于显示数据处理结果，实现报警信号发布；

数据处理模块将振动加速度信号经过处理后与预设的报警阈值进行比较后，若结果显示振动烈度值在预设的报警阈值以内，则将数据发送至数据存储模块进行存储，否则进行故障报警并根据振动烈度所在区域更新历史数据库。

4.根据权利要求1所述的方法，其特征在于，所述运转设备是水泵、风机、空压机、电机、电梯曳引机或电梯轿厢中的任意一种。

5.根据权利要求1所述的方法，其特征在于，所述数据存储模块是基于本地或云端的数据存储服务器，所述数据处理模块是数据处理服务器；所述数据输出模块是手机、电脑或智能显示屏中的任意一种。

6.根据权利要求1所述的方法，其特征在于，所述数据采集模块由振动检测传感器和无线网关组成，其中振动检测传感器用于采集运转设备的振动信号并通过无线网关上传至数据处理服务器。

7.根据权利要求1所述的方法，其特征在于，所述振动检测传感器是有线振动传感器、LoRa振动传感器、NB-IoT振动传感器或WiFi振动传感器中的任意一种；所述无线网关是LoRa无线网关或Wifi无线网关中的任意一种。

8.基于权利要求1所述方法的运转设备安康报警方法，其特征在于，进一步包括：由数据处理模块对实时采集到的振动加速度信号进行分析，如其频率分量的幅值U_i落在A或B区域内，认为设备正常运行；若频率分量的幅值U_i落在C或D区域内，则认为设备处于故障状态，此时由数据输出模块进行报警，并将落在C区域的频率分量的幅值U_i对应的原始数据通过数据处理模块更新至历史数据库；对于落在D区域的原始数据则直接丢弃。

9.基于权利要求1所述方法的运转设备安康报警方法，其特征在于，进一步包括：由数据处理模块根据预置规则对经过处理的频率分量的幅值U_i进行比较和判断，确认并输出故障类型诊断结果。

Images