CN108731923B - 一种旋转机械设备的故障检测方法及装置 - Google Patents
一种旋转机械设备的故障检测方法及装置 Download PDFInfo
- Publication number
- CN108731923B CN108731923B CN201810267031.1A CN201810267031A CN108731923B CN 108731923 B CN108731923 B CN 108731923B CN 201810267031 A CN201810267031 A CN 201810267031A CN 108731923 B CN108731923 B CN 108731923B
- Authority
- CN
- China
- Prior art keywords
- vibration signal
- fault
- data
- value
- mechanical equipment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本申请实施例公开了一种旋转机械设备的故障检测方法及装置,用于实现旋转机械设备的故障的准确检测,该方法包括:基于历史数据预先构建用于故障判断的专家决策系统以及用于故障分类的故障分类模型,在获得旋转机械设备的原始振动信号后,可以提取到振动信号特征值,将振动信号特征值以及采集到的旋转机械设备的设备数据输入专家决策系统,可以得到旋转机械设备是否存在故障的故障判断结果,进一步的,当存在故障时,将振动信号特征值以及采集到的旋转机械设备的设备数据输入故障分类模型可以得到故障分类结果。本申请实施例利用人工智能算法可以实现旋转机械设备的智能故障诊断及分类,提高了故障检测的效率以及有效性。
Description
技术领域
本申请涉及工业设备故障检测技术领域,具体涉及一种旋转机械设备的故障检测方法及装置。
背景技术
在工业生产背景中,需要对旋转机械设备例如汽轮机的工作状态进行监控。传统的旋转机械设备状态监控采用的是固定阈值监测和故障判断,该方式可以通用性解决一类状态的监测问题。但是随着被监测状态和参数的复杂化、所处行业生产过程的复杂化以及其他非一致性的客观条件,例如设备老化、型号差异、生产过程差异等,基于传统的固定阈值监测以实现故障检测会存在故障误报、故障漏报等情况,无法做到有效地故障检测。
在现有技术中,还可以采用有经验的人员配合现有的故障检测系统实现人工干预和预判从而进行有效的故障检测,但是增加了人力成本且监测效果受主观因素影响较大,具有一定的局限性。
发明内容
有鉴于此,本申请实施例提供一种旋转机械设备的故障检测方法及装置,以解决现有技术中无法有效进行故障检测的技术问题。
为解决上述问题,本申请实施例提供的技术方案如下:
一种旋转机械设备的故障检测方法,包括:
获取旋转机械设备的原始振动信号;
从所述原始振动信号中提取振动信号特征值;
将所述振动信号特征值以及采集到的所述旋转机械设备的设备数据输入专家决策系统,获得所述旋转机械设备是否存在故障的故障判断结果,所述专家决策系统是根据历史数据、所述历史数据对应的故障判断结果以及专家经验知识构建的,所述历史数据包括历史振动信号特征值以及历史设备数据;
当获取到所述旋转机械设备存在故障的故障判断结果,将所述振动信号特征值以及所述设备数据输入故障分类模型,获得故障分类结果,所述故障分类模型是根据所述历史数据以及所述历史数据对应的故障分类结果训练生成的。
在一种可选的实现方式中,所述方法还包括:
将所述振动信号特征值以及所述设备数据作为所述历史数据,更新所述故障分类模型和/或所述专家决策系统。
在一种可选的实现方式中,所述方法还包括:
根据所述原始振动信号以及自回归滑动平均模型对通频振动幅值超过阈值的时间点进行预测,将所述时间点作为所述旋转机械设备的寿命估计值。
在一种可选的实现方式中,所述从所述原始振动信号中提取振动信号特征值,包括:
将所述原始振动信号进行预处理,获得预处理后的振动信号,所述预处理包括单位转换处理以及积分处理;
将所述预处理后的振动信号进行低通滤波处理,获得去除噪声的振动信号,所述去除噪声的振动信号包括振动信号时域数据;
将所述去除噪声的振动信号进行快速傅里叶变换,得到振动信号频域数据,所述振动信号频域数据包括振动信号幅值数据、振动信号相位数据,将所述振动信号幅值数据、所述振动信号相位数据以及所述振动信号时域数据确定为振动信号特征值。
在一种可选的实现方式中,所述专家决策系统采用树状结构,所述树状结构的叶子节点输出所述旋转机械设备是否存在故障的故障判断结果,所述树状结构的各个非叶子节点记录有一个属性特征的程度值,所述树状结构的各个非叶子节点的位置根据各个非叶子节点对应的属性特征的信息增益确定。
对应于上述旋转机械设备的故障检测方法,本申请提出了一种旋转机械设备的故障检测装置,所述装置包括:
获取单元,用于获取旋转机械设备的原始振动信号;
提取单元,用于从所述原始振动信号中提取振动信号特征值;
判断单元,用于将所述振动信号特征值以及采集到的所述旋转机械设备的设备数据输入专家决策系统,获得所述旋转机械设备是否存在故障的故障判断结果,所述专家决策系统是根据历史数据、所述历史数据对应的故障判断结果以及专家经验知识构建的,所述历史数据包括历史振动信号特征值以及历史设备数据;
分类单元,用于当获取到所述旋转机械设备存在故障的故障判断结果,将所述振动信号特征值以及所述设备数据输入故障分类模型,获得故障分类结果,所述故障分类模型是根据所述历史数据以及所述历史数据对应的故障分类结果训练生成的。
在一种可选的实现方式中,所述装置还包括:
更新单元,用于将所述振动信号特征值以及所述设备数据作为所述历史数据,更新所述故障分类模型和/或所述专家决策系统。
在一种可选的实现方式中,所述装置还包括:
预测单元,用于根据所述原始振动信号以及自回归滑动平均模型对通频振动幅值超过阈值的时间点进行预测,将所述时间点作为所述旋转机械设备的寿命估计值。
在一种可选的实现方式中,所述提取单元具体包括:
预处理单元,用于将所述原始振动信号进行预处理,获得预处理后的振动信号,所述预处理包括单位转换处理以及积分处理;
低通滤波处理单元,用于将所述预处理后的振动信号进行低通滤波处理,获得去除噪声的振动信号,所述去除噪声的振动信号包括振动信号时域数据;
确定单元,用于将所述去除噪声的振动信号进行快速傅里叶变换,得到振动信号频域数据,所述振动信号频域数据包括振动信号幅值数据、振动信号相位数据,将所述振动信号幅值数据、所述振动信号相位数据以及所述振动信号时域数据确定为振动信号特征值。
在一种可选的实现方式中,所述专家决策系统采用树状结构,所述树状结构的叶子节点输出所述旋转机械设备是否存在故障的故障判断结果,所述树状结构的各个非叶子节点记录有一个属性特征的程度值,所述树状结构的各个非叶子节点的位置根据各个非叶子节点对应的属性特征的信息增益确定。
由此可见,本申请实施例具有如下有益效果:
本申请实施例基于历史数据预先构建了用于故障判断的专家决策系统以及用于故障分类的故障分类模型,在获得旋转机械设备的原始振动信号后,可以提取到振动信号特征值,将振动信号特征值以及采集到的旋转机械设备的设备数据输入专家决策系统,可以得到旋转机械设备是否存在故障的故障判断结果,进一步的,当存在故障时,将振动信号特征值以及采集到的旋转机械设备的设备数据输入故障分类模型可以得到故障分类结果。本申请实施例利用人工智能算法可以实现旋转机械设备的智能故障诊断及分类,提高了故障检测的效率以及有效性。
附图说明
图1为本申请实施例提供的一种旋转机械设备的故障检测方法实施例的流程图;
图2为本申请实施例提供的从原始振动信号中提取振动信号特征值的流程图;
图3为本申请实施例提供的专家决策系统树状结构示例图;
图4为本申请实施例提供的一种旋转机械设备的故障检测方法的结构示意图;
图5为本申请实施例提供的一种旋转机械设备的故障检测装置实施例的示意图。
具体实施方式
为使本申请的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本申请实施例作进一步详细的说明。
随着现代工业的快速发展,旋转式机械设备作为驱动各种生产设备运转的原动机,一旦其发生故障而停止工作,将导致整个生产过程的停滞,造成不可估量的经济损失,因此,对其作业实时状态监测及故障检测的研究也变得尤其重要。传统的旋转机械设备状态的监测采用的是固定阈值监测和故障判断,可以通用性解决一类状态的监测问题,但是随着设备本身的老化及工业生产过程差异等自身及外界因素的影响,传统的监测方法已无法有效地对旋转式机械设备故障的进行精准检测。
基于此,本申请提出了一种旋转机械设备的故障检测方法和装置,实现对旋转机械设备的智能的故障诊断及分类。
以下将结合附图对本申请实施例提供的旋转机械设备的故障检测方法进行详细说明。参见图1所示,其示出了本申请实施例提供的一种旋转机械设备的故障检测方法实施例的流程图,本实施例可以包括以下步骤:
步骤101:获取旋转机械设备的原始振动信号。
在本申请实施例中,采用人工智能的方式,结合了历史振动信号与专家系统实现了对旋转机械设备的智能故障检测。为了有效检测旋转机械设备的故障,首先获取旋转机械设备的原始振动信号,其中,旋转机械设备的原始振动信号指的是以旋转机械设备的转速为基频的周期信号。根据实际应用情况,采集的原始振动信号可以分为:轴承振动评定,即,被检测点位于轴承基座;轴振动值评定,即被检测点位于基座上,轴的两侧。当故障发生时,原始振动信号是非稳定和非线性的信号,通过对这些信号进行分析和处理,便于后续判断旋转机械设备的故障结果。
在获取到旋转机械设备的原始振动信号后,可继续执行步骤102。
步骤102:从原始振动信号中提取振动信号特征值。
在实际应用中,在获取到旋转机械设备的原始振动信号后,将从这些原始振动信号中提取出振动信号特征值,在进行旋转机械设备的故障检测过程中,提取出的振动信号特征值是否能准确表征旋转机械设备的不同故障对应的不同振动信号对于检测结果是非常重要的。以下结合附图对该振动信号特征值提取过程进行介绍。
参见图2所示,其示出了本申请实施例提供的从原始振动信号中提取振动信号特征值的流程图,具体可以包括以下步骤:
步骤201:将原始振动信号进行预处理,获得预处理后的振动信号,该预处理包括单位转换处理以及积分处理。
在具体实现过程中,在获取到旋转机械设备的原始振动信号后,将结合震动频率及灵敏度系数,对这些原始振动信号进行单位转换和积分等处理,将其变换成为被检测点的相对位移量,速度量,加速度量等物理单位,同时将振动信号的平均值输出,作为轴心位置的一个分量。接着,执行步骤202.
步骤202:将预处理后的振动信号进行低通滤波处理,获得去除噪声的振动信号,去除噪声的振动信号包括振动信号时域数据。
在具体实现过程中,在对原始振动信号进行预处理后,可得到预处理后的振动信号,但由于在振动信号采集过程中,会受到现场生产车间各种因素的影响,采集的信号数据不可避免的混有噪声,有时,噪声甚至可以把有用的信息完全淹没,因此,需要将预处理后的振动信号进行低通滤波处理。本申请通过设计低通滤波器的各个参数,并采用无相位滞后滤波算法,获得去除噪声的振动信号,其中,去除噪声的振动信号包括振动信号时域数据。在获得去除噪声的振动信号后,可执行步骤203。
步骤203:将去除噪声的振动信号进行快速傅里叶变换,得到振动信号频域数据,该振动信号频域数据包括振动信号幅值数据、振动信号相位数据,将该振动信号幅值数据、振动信号相位数据以及振动信号时域数据确定为振动信号特征值。
在具体实现过程中,获得去除噪声的振动信号后,为了提取振动信号的特征值,将对去除噪声的振动信号进行快速傅里叶变换,本申请采用快速傅里叶变换(Fast FourierTransform,简称FFT),实现将振动信号从时域到频域的转换,从而得到振动信号的幅值、相位等频域数据,并将该振动信号离散的幅值数据、振动信号相位数据以及原来的时域信号数据确定为振动信号特征值。
在本申请实施例中,在从原始振动信号中提取到振动信号特征值之后,可执行步骤103。
步骤103:将该振动信号特征值以及采集到的旋转机械设备的设备数据输入专家决策系统,获得旋转机械设备是否存在故障的故障判断结果,其中,该专家决策系统是根据历史数据、历史数据对应的故障判断的结果以及专家经验知识构建的,该历史数据包括历史振动信号特征值以及历史设备数据。
在本申请实施例中,采用人工智能的方式,结合了历史振动信号与专家系统实现了对旋转机械设备的智能故障检测。
在实际应用中,在从原始振动信号中提取到振动信号特征值后,将该振动信号特征值以及采集到的旋转机械设备的设备数据输入专家决策系统,获得旋转机械设备是否存在故障的故障判断结果,其中,旋转机械设备的设备数据指的是设备自身温度以及所处的外界环境温度、设备的压力以及转速等设备数据。而专家决策系统指的是根据历史数据、历史数据对应的故障判断的结果以及专家经验知识构建的,该历史数据包括历史振动信号特征值以及历史设备数据,这些历史振动信号特征值以及历史设备数据对应的是旋转机械设备是否存在故障的故障判断结果。历史振动信号特征值是从历史振动信号中提取的,提取的方式与从原始振动信号中提取振动信号特征值类似,在此不再赘述。
在本申请实施例中,一种可选的实施方式是,本申请中的专家决策系统采用树状结构,实现对故障发生与否的逻辑判断,其中,该树状结构的叶子节点输出旋转机械设备是否存在故障的故障判断结果,该树状结构的各个非叶子节点记录有一个属性特征的程度值,该树状结构的各个非叶子节点的位置根据各个非叶子节点对应的属性特征的信息增益确定。
在具体实现过程中,专家决策系统树状结构的叶子结点值只能为1或0,表示故障发生或运行正常。树上的各个非叶子节点记录有一个属性特征的程度值。并且结点的每个后继分支对应于该属性的一个可能值或范围。属性结点的参数设定依据人工专家的经验,并且可以根据设备实际运行工况或设备型号专门配置。
专家决策系统的决策过程是从这棵树的根结点开始的,测试该结点的指定属性,然后根据该结点分支的逻辑属性向下移动。然后该过程在以新结点为根的子树上重复。最终到达叶子节点是即可得出故障发生与否的判断。
专家决策系统树状结构举例如图3所示,其示出了本申请实施例提供的专家决策系统树状结构示例图,
举例说明:如图3所示,图中非叶子节点记录有一个属性特征为1×频谱,该属性特征的程度值为20,其他各个非叶子节点的位置根据各个非叶子节点对应的1×频谱这个属性特征的信息增益确定,越靠近该结点属性结点的信息增益越大,例如轴承乌金温度则是信息增益最大等,依次类推,直至到达叶子节点,即可得出故障发生与否的判断,例如当推力瓦块金属温度>50度时,可得出故障发生的故障判断结果,当推力瓦块金属温度>50度不成立时,可得出设备是正常运行的故障判断结果。
本申请中专家决策系统决策树的构造过程不依赖领域知识,它使用属性选择度量来选择将元组最好地划分成不同的类的属性。所谓决策树的构造就是进行属性选择度量确定各个特征属性之间的拓扑结构。
在实际应用中,为了加快决策判断过程,本申请的专家决策系统模型是根据专家经验知识判断其他属性特征的信息增益的大小,并根据信息增益的大小来布置结点属性的分布,越靠近根结点属性结点信息增益越大,构建最佳的故障判断决策系统。为了精确地定义信息增益,本算法采用信息论中广泛使用的一个度量标准,称为熵(Entropy),它刻画了任意属性的纯度。
在具体实现过程中,决策树的关键步骤是分裂属性。所谓分裂属性就是在某个节点处按照某一特征属性的不同划分构造不同的分支,其目标是让各个分裂子集尽可能地“纯”。尽可能“纯”就是尽量让一个分裂子集中待分类项属于同一类别。
期望信息越小,信息增益越大,从而纯度越高。算法的核心思想就是以信息增益度量属性选择,选择分裂后信息增益最大的属性进行分裂。下面先定义几个要用到的概念。
设D为用类别对训练元组进行的划分,则D的熵(info)表示为:
其中pi表示第i个类别在整个训练元组中出现的概率,可以用属于此类别元素的数量除以训练元组元素总数量作为估计。熵的实际意义表示是D中元组的类标号所需要的平均信息量。
现在我们假设将训练元组D按属性A进行划分,则A对D划分的期望信息为:
而信息增益即为两者的差值:
gain(A)=inf o(D)-inf oA(D)
每次需要分裂时,计算每个属性的增益率,然后选择增益率最大的属性进行分裂。
先将D中元素按照特征属性排序,则每两个相邻元素的中间点可以看做潜在分裂点,从第一个潜在分裂点开始,分裂D并计算两个集合的期望信息,具有最小期望信息的点称为这个属性的最佳分裂点,其信息期望作为此属性的信息期望。
对于非离散型连续数据,首先要完成“分类”将其完成成为离散数据,形成数据记录元祖,将数据记录元祖泛化形成符合期望的特征值元祖。
实际构建决策树过程中,需要结合实际故障率和故障分布决定整个决策树叶子节点“繁盛”程度,即实际最终状态故障叶子节点的实际对特征占比覆盖应该符合实际故障率的统计值和经验值。
针对记录元祖故障分类后泛化形成的特征值元祖,将会产生不同差异的故障分布,根据该分布决定决策树深度以及分裂过程中粒度,及分裂深度,以拟合实际故障分布。
在本申请实施例中,当通过专家决策诊断系统获取到旋转机械设备存在故障的故障判断结果之后,可执行步骤104。
步骤104:当获取到旋转机械设备存在故障的故障判断结果,将振动信号特征值以及设备数据输入故障分类模型,获得故障分类结果,该故障分类模型是根据历史数据以及历史数据对应的故障分类结果训练生成的。
在实际应用中,当通过专家决策系统获取到旋转机械设备存在故障的故障判断结果后,将振动信号特征值以及设备数据输入故障分类模型,获得故障分类结果,其中,故障分类模型是根据历史数据以及历史数据对应的故障分类结果训练生成的。本申请在故障分类模型中采用的故障分类算法可以是支持向量机(Support Vector Machine,简称SVM)的多分类算法。该算法在学习复杂的非线性模式时提供了一种更为清晰,更加强大的方式。在分类应用中能实现在有足够域度的情况下实现精确且安全的故障模式分类。
在本申请实施例中,采用的故障分类算法的计算公式如下:
hθ(x)=1;ifθTx≥0
hθ(x)=0;ifθTx<0
其中,hθ(x)函数是对于样本x取值的概率函数,表示其相似该特征故障的概率,而hθ(x)=1则表示明确某种故障的分类,参数θ是故障分类模型的待估计参数。
故障分类模型自学习校正算法通过如下的最小化代价函数得到参数θ:
cost1(θTx(i))=-log hθ(x(i))
cost0(θTx(i))=-log(1-hθ(x(i)))
其中,cost函数为代价函数,用来估计样本属于某类的风险程度,其值越小代表越有可能属于这类,x代表每个样本数据点在某一个特征上的值,即特征向量x的某个值,y则表示每个样本数据的所属类别标签。
对于多故障分类问题可以将上诉二分类算法应用多次,实现多次的二分类,最终实现多故障分类。
在具体实现过程中,在获得故障分类结果后,还可以将设备的故障分类报告至工程师,以便于工程师对发生故障的设备进行维修并确认设备的故障类型,同时对故障分类模型进行修正,在故障分类结果确认后,本申请将正确的故障分类结果以及故障振动信号特征值、设备数据作为历史数据加入数据库,实现对数据库的实时更新。
在本申请一些可选的实现方式中,本申请还包括:将所述振动信号特征值以及所述设备数据作为所述历史数据,更新所述故障分类模型和/或所述专家决策系统。
在实际应用中,当通过本申请的旋转机械设备的故障检测方法获得设备故障的分类结果后,可以将各故障分类结果对应的振动信号特征值以及设备数据作为历史数据,结合对应的故障判断结果对专家决策系统进行更新;同时,也可以将各故障分类结果对应的振动信号特征值以及设备数据作为历史数据,结合对应的故障分类结果对故障分类模型进行更新,以便于后续进行更有效的旋转机械设备的故障检测。
在本申请一些可选的实现方式中,本申请还包括:根据所述原始振动信号以及自回归滑动平均模型对通频振动幅值超过阈值的时间点进行预测,将所述时间点作为所述旋转机械设备的寿命估计值。
在实际应用中,当通过本申请的旋转机械设备的故障检测方法获得设备故障的分类结果后,还可以根据原始振动信号以及自回归滑动平均模型对通频振动幅值超过阈值的时间点进行预测,将该时间点作为旋转机械设备的寿命估计值,并且对于关键设备或元件在寿命小于一定阈值时进行提前的预报警。
其中,通频振动幅值表示的是振动原始波形的振动幅值。而自回归滑动平均模型(Auto-Regressive Moving Average Model,简称ARMA)则是一种常用的传统时间序列分析模型。对于依赖于时间t的一组随机变量,此处选择去除周期成分的振动信号,其整个序列的变化具有一定的规律性,用相应的数学模型加以描述,通过对该数学模型的分析研究,能够更本质的认识时间序列的结构与特征,达到最小方差意义下的预测。
自回归滑动平均序列{xt}模型如下公式所示:
接着,采用最小二乘法对上述模型进行参数辨识,然后计算残差序列的自相关函数来对模型进行检验。对于自相关函数值较小的自回归滑动平均模型,即为对原始振动序列随机信号部分的很好近似,因此,基于该模型和振动信号的固定周期成分即可计算对设备故障时间进行预测,从而得到旋转机械设备的寿命估计值,当对于关键设备或元件在寿命小于一定阈值时,进行提前的预报警。
这样,本申请实施例基于历史数据预先构建了用于故障判断的专家决策系统以及用于故障分类的故障分类模型,在获得旋转机械设备的原始振动信号后,可以提取到振动信号特征值,将振动信号特征值以及采集到的旋转机械设备的设备数据输入专家决策系统,可以得到旋转机械设备是否存在故障的故障判断结果,进一步的,当存在故障时,将振动信号特征值以及采集到的旋转机械设备的设备数据输入故障分类模型可以得到故障分类结果。本申请实施例利用人工智能算法可以实现旋转机械设备的智能故障诊断及分类,提高了故障检测的效率以及有效性。
为便于理解,现结合图4所示一种旋转机械设备的故障检测方法的结构示意图,对本申请实施例提供的旋转机械设备的故障检测方法的实现过程进行整体性介绍。
如图4所示,本申请实施例的整体实现过程为:首先,获取旋转机械设备的原始振动信号,接着,采用预处理、低通滤波以及快速傅里叶变换算法从原始振动信号中提取出振动信号的幅频、相频和时频作为振动信号特征值,具体实现过程参见步骤101~步骤102;然后,将振动信号特征值以及采集到的旋转机械设备的设备数据输入专家决策系统,获得旋转机械设备是否存在故障的故障判断结果,具体实现过程参见步骤103。最后,当获取到旋转机械设备存在故障的故障判断结果后,将振动信号特征值以及设备数据输入故障分类模型,获得故障分类结果,具体实现过程参见步骤104。同时,在获得故障分类结果后,可将结果上报给工程师进行确认和维修。在故障类型确认后,本申请还可以将正确故障分类结果与振动信号特征值以及设备数据建立一定的对应关系并加入数据库,进而对故障分类模型进行加强和修正,实现旋转机械设备的故障检测的自学习功能。经过反复的诊断、学习迭代过程,不断完善该检测方法,提高检测精确度。进一步的,在获得故障分类结果后,本申请还可以根据原始振动信号以及自回归滑动平均模型对通频振动幅值超过阈值的时间点进行预测,将该时间点作为旋转机械设备的寿命估计值,对于关键设备或元件在寿命小于一定阈值时进行提前的预报警。最终可以实现精度远高于传统检测方法的旋转机械设备的故障检测方法以及设备(元件)寿命的预测。
参见图5所示,本申请还提供一种旋转机械设备的故障检测装置实施例,该装置包括:
获取单元501,用于获取旋转机械设备的原始振动信号。
提取单元502,用于从所述原始振动信号中提取振动信号特征值。
判断单元503,用于将所述振动信号特征值以及采集到的所述旋转机械设备的设备数据输入专家决策系统,获得所述旋转机械设备是否存在故障的故障判断结果,所述专家决策系统是根据历史数据、所述历史数据对应的故障判断结果以及专家经验知识构建的,所述历史数据包括历史振动信号特征值以及历史设备数据。
分类单元504,用于当获取到所述旋转机械设备存在故障的故障判断结果,将所述振动信号特征值以及所述设备数据输入故障分类模型,获得故障分类结果,所述故障分类模型是根据所述历史数据以及所述历史数据对应的故障分类结果训练生成的。
在本申请一些可能的实现方式中,所述装置还包括:
更新单元,用于将所述振动信号特征值以及所述设备数据作为所述历史数据,更新所述故障分类模型和/或所述专家决策系统。
在本申请一些可能的实现方式中,所述装置还包括:
预测单元,用于根据所述原始振动信号以及自回归滑动平均模型对通频振动幅值超过阈值的时间点进行预测,将所述时间点作为所述旋转机械设备的寿命估计值。
在本申请一些可能的实现方式中,所述提取单元502具体包括:
预处理单元,用于将所述原始振动信号进行预处理,获得预处理后的振动信号,所述预处理包括单位转换处理以及积分处理;
低通滤波处理单元,用于将所述预处理后的振动信号进行低通滤波处理,获得去除噪声的振动信号,所述去除噪声的振动信号包括振动信号时域数据;
确定单元,用于将所述去除噪声的振动信号进行快速傅里叶变换,得到振动信号频域数据,所述振动信号频域数据包括振动信号幅值数据、振动信号相位数据,将所述振动信号幅值数据、所述振动信号相位数据以及所述振动信号时域数据确定为振动信号特征值。
在本申请一些可能的实现方式中,所述专家决策系统采用树状结构,所述树状结构的叶子节点输出所述旋转机械设备是否存在故障的故障判断结果,所述树状结构的各个非叶子节点记录有一个属性特征的程度值,所述树状结构的各个非叶子节点的位置根据各个非叶子节点对应的属性特征的信息增益确定。
由上述实施例可以看出,本申请实施例基于历史数据预先构建了用于故障判断的专家决策系统以及用于故障分类的故障分类模型,在获得旋转机械设备的原始振动信号后,可以提取到振动信号特征值,将振动信号特征值以及采集到的旋转机械设备的设备数据输入专家决策系统,可以得到旋转机械设备是否存在故障的故障判断结果,进一步的,当存在故障时,将振动信号特征值以及采集到的旋转机械设备的设备数据输入故障分类模型可以得到故障分类结果。本申请实施例利用人工智能算法可以实现旋转机械设备的智能故障诊断及分类,提高了故障检测的效率以及有效性。
需要说明的是,本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统或装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
还需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本申请。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下,在其它实施例中实现。因此,本申请将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (8)
1.一种旋转机械设备的故障检测方法,其特征在于,所述方法包括:
获取旋转机械设备的原始振动信号;
从所述原始振动信号中提取振动信号特征值;
将所述振动信号特征值以及采集到的所述旋转机械设备的设备数据输入专家决策系统,获得所述旋转机械设备是否存在故障的故障判断结果,所述专家决策系统是根据历史数据、所述历史数据对应的故障判断结果以及专家经验知识构建的,所述历史数据包括历史振动信号特征值以及历史设备数据;
当获取到所述旋转机械设备存在故障的故障判断结果,将所述振动信号特征值以及所述设备数据输入故障分类模型,获得故障分类结果,所述故障分类模型是根据所述历史数据以及所述历史数据对应的故障分类结果训练生成的;所述故障分类模型包括故障分类算法和故障分类模型自学习校正算法,所述故障分类算法的计算公式如下:
hθ(x)=1;ifθTx≥0
hθ(x)=0;ifθTx<0
其中,hθ(x)函数是对于样本x取值的概率函数,表示其相似特征故障的概率,hθ(x)=1表示明确某种故障的分类,参数θ是故障分类模型的待估计参数;
所述故障分类模型自学习校正算法通过如下的最小化代价函数得到参数θ:
cost1(θTx(i))=-loghθ(x(i))
cost0(θTx(i))=-log(1-hθ(x(i)))
其中,cost函数为代价函数,用来估计样本属于某类的风险程度,其值越小代表越有可能属于这类,x代表每个样本数据点在某一个特征上的值,即特征向量x的某个值,y则表示每个样本数据的所属类别标签;
其中,将所述振动信号特征值以及所述设备数据作为所述历史数据,结合对应的故障分类结果更新所述故障分类模型和/或所述专家决策系统。
2.根据权利要求1所述的方法,其特征在于,所述方法还包括:
根据所述原始振动信号以及自回归滑动平均模型对通频振动幅值超过阈值的时间点进行预测,将所述时间点作为所述旋转机械设备的寿命估计值。
3.根据权利要求1所述的方法,其特征在于,所述从所述原始振动信号中提取振动信号特征值,包括:
将所述原始振动信号进行预处理,获得预处理后的振动信号,所述预处理包括单位转换处理以及积分处理;
将所述预处理后的振动信号进行低通滤波处理,获得去除噪声的振动信号,所述去除噪声的振动信号包括振动信号时域数据;
将所述去除噪声的振动信号进行快速傅里叶变换,得到振动信号频域数据,所述振动信号频域数据包括振动信号幅值数据、振动信号相位数据,将所述振动信号幅值数据、所述振动信号相位数据以及所述振动信号时域数据确定为振动信号特征值。
4.根据权利要求1所述的方法,其特征在于,所述专家决策系统采用树状结构,所述树状结构的叶子节点输出所述旋转机械设备是否存在故障的故障判断结果,所述树状结构的各个非叶子节点记录有一个属性特征的程度值,所述树状结构的各个非叶子节点的位置根据各个非叶子节点对应的属性特征的信息增益确定。
5.一种旋转机械设备的故障检测装置,其特征在于,所述装置包括:
获取单元,用于获取旋转机械设备的原始振动信号;
提取单元,用于从所述原始振动信号中提取振动信号特征值;
判断单元,用于将所述振动信号特征值以及采集到的所述旋转机械设备的设备数据输入专家决策系统,获得所述旋转机械设备是否存在故障的故障判断结果,所述专家决策系统是根据历史数据、所述历史数据对应的故障判断结果以及专家经验知识构建的,所述历史数据包括历史振动信号特征值以及历史设备数据;
分类单元,用于当获取到所述旋转机械设备存在故障的故障判断结果,将所述振动信号特征值以及所述设备数据输入故障分类模型,获得故障分类结果,所述故障分类模型是根据所述历史数据以及所述历史数据对应的故障分类结果训练生成的;所述故障分类模型包括故障分类算法和故障分类模型自学习校正算法,所述故障分类算法的计算公式如下:
hθ(x)=1;ifθTx≥0
hθ(x)=0;ifθTx<0
其中,hθ(x)函数是对于样本x取值的概率函数,表示其相似特征故障的概率,hθ(x)=1表示明确某种故障的分类,参数θ是故障分类模型的待估计参数;
所述故障分类模型自学习校正算法通过如下的最小化代价函数得到参数θ:
cost1(θTx(i))=-loghθ(x(i))
cost0(θTx(i))=-log(1-hθ(x(i)))
其中,cost函数为代价函数,用来估计样本属于某类的风险程度,其值越小代表越有可能属于这类,x代表每个样本数据点在某一个特征上的值,即特征向量x的某个值,y则表示每个样本数据的所属类别标签;
其中,将所述振动信号特征值以及所述设备数据作为所述历史数据,结合对应的故障分类结果更新所述故障分类模型和/或所述专家决策系统。
6.根据权利要求5所述的装置,其特征在于,所述装置还包括:
预测单元,用于根据所述原始振动信号以及自回归滑动平均模型对通频振动幅值超过阈值的时间点进行预测,将所述时间点作为所述旋转机械设备的寿命估计值。
7.根据权利要求5所述的装置,其特征在于,所述提取单元具体包括:
预处理单元,用于将所述原始振动信号进行预处理,获得预处理后的振动信号,所述预处理包括单位转换处理以及积分处理;
低通滤波处理单元,用于将所述预处理后的振动信号进行低通滤波处理,获得去除噪声的振动信号,所述去除噪声的振动信号包括振动信号时域数据;
确定单元,用于将所述去除噪声的振动信号进行快速傅里叶变换,得到振动信号频域数据,所述振动信号频域数据包括振动信号幅值数据、振动信号相位数据,将所述振动信号幅值数据、所述振动信号相位数据以及所述振动信号时域数据确定为振动信号特征值。
8.根据权利要求5所述的装置,其特征在于,所述专家决策系统采用树状结构,所述树状结构的叶子节点输出所述旋转机械设备是否存在故障的故障判断结果,所述树状结构的各个非叶子节点记录有一个属性特征的程度值,所述树状结构的各个非叶子节点的位置根据各个非叶子节点对应的属性特征的信息增益确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810267031.1A CN108731923B (zh) | 2018-03-28 | 2018-03-28 | 一种旋转机械设备的故障检测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810267031.1A CN108731923B (zh) | 2018-03-28 | 2018-03-28 | 一种旋转机械设备的故障检测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108731923A CN108731923A (zh) | 2018-11-02 |
CN108731923B true CN108731923B (zh) | 2020-11-03 |
Family
ID=63940998
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810267031.1A Active CN108731923B (zh) | 2018-03-28 | 2018-03-28 | 一种旋转机械设备的故障检测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108731923B (zh) |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6950670B2 (ja) * | 2018-12-12 | 2021-10-13 | 横河電機株式会社 | 検出装置、検出方法、および、検出プログラム |
DE102019127211A1 (de) * | 2019-03-05 | 2020-09-10 | Computational Systems, Inc. | System zum Trennen von periodischen Amplitudenspitzen von nicht-periodischen Amplitudenspitzen in Maschinenschwingungsdaten |
CN110108431B (zh) * | 2019-05-22 | 2021-07-16 | 西安因联信息科技有限公司 | 一种基于机器学习分类算法的机械设备故障诊断方法 |
CN110119787B (zh) * | 2019-05-23 | 2021-07-20 | 湃方科技(北京)有限责任公司 | 一种旋转型机械设备工况检测方法及设备 |
CN110426220A (zh) * | 2019-05-23 | 2019-11-08 | 中国航空工业集团公司上海航空测控技术研究所 | 基于自适应滤波电路的机械故障监测系统 |
CN110411766A (zh) * | 2019-07-30 | 2019-11-05 | 中国神华能源股份有限公司神朔铁路分公司 | 列车转向架蛇形失稳检测方法、装置、系统和存储介质 |
CN110864887A (zh) * | 2019-11-19 | 2020-03-06 | 北京瑞莱智慧科技有限公司 | 机械设备运转状况确定方法、装置、介质和计算设备 |
CN111189646A (zh) * | 2019-12-16 | 2020-05-22 | 上海蔚来汽车有限公司 | 车辆nvh自诊断方法、装置、车辆以及控制器和介质 |
CN111079705B (zh) * | 2019-12-31 | 2023-07-25 | 北京理工大学 | 一种振动信号分类方法 |
CN111255674B (zh) * | 2020-01-21 | 2022-08-09 | 武汉瑞莱保科技有限公司 | 一种用于检测旋转机械设备状态的系统及检测方法 |
CN113155271B (zh) * | 2020-01-23 | 2023-08-22 | 上海擎动信息科技有限公司 | 声振检测方法、系统、终端及介质 |
CN111256757A (zh) * | 2020-02-25 | 2020-06-09 | 深圳哈维生物医疗科技有限公司 | 基于云计算的医疗设备监控系统及方法 |
CN111337237B (zh) * | 2020-03-17 | 2021-10-15 | 北京必可测科技股份有限公司 | 一种设备故障的诊断方法及系统 |
CN112541563A (zh) * | 2020-09-30 | 2021-03-23 | 国电龙源电力技术工程有限责任公司 | 基于边缘计算技术的旋转设备振动预测管理系统 |
CN112945535B (zh) * | 2021-02-20 | 2023-04-18 | 广东石油化工学院 | 一种基于数值模拟的旋转机械故障检测方法及装置 |
JP2022188345A (ja) * | 2021-06-09 | 2022-12-21 | 富士電機株式会社 | 診断装置、診断方法、診断プログラム |
CN113551764B (zh) * | 2021-07-29 | 2024-05-03 | 西门子工厂自动化工程有限公司 | 振动分析设备及其方法 |
CN114036998A (zh) * | 2021-09-24 | 2022-02-11 | 浪潮集团有限公司 | 一种基于机器学习对工业硬件进行故障检测的方法及系统 |
CN114935357A (zh) * | 2022-03-14 | 2022-08-23 | 浙江倍时信息科技有限公司 | 一种基于熵值计算的设备健康监测系统 |
CN116910680B (zh) * | 2023-09-11 | 2023-12-05 | 江苏优创生物医学科技有限公司 | 一种健身器材的远程故障检测方法及系统 |
CN117113200B (zh) * | 2023-10-24 | 2024-02-02 | 中海石油气电集团有限责任公司 | 转子故障诊断方法、装置、电子设备及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101025856A (zh) * | 2007-02-08 | 2007-08-29 | 浙江大学 | 一种面向低频振动检测的数据采集装置及方法 |
CN101672723A (zh) * | 2009-10-28 | 2010-03-17 | 北京中能联创风电技术有限公司 | 一种风电机组振动分析故障诊断方法和系统 |
CN102156043A (zh) * | 2010-12-31 | 2011-08-17 | 北京四方继保自动化股份有限公司 | 风力发电机组在线状态监测与故障诊断系统 |
CN104503399A (zh) * | 2014-12-09 | 2015-04-08 | 华电电力科学研究院 | 一种集团级风电机组状态监测及故障诊断平台 |
CN104932519A (zh) * | 2015-05-25 | 2015-09-23 | 北京航空航天大学 | 基于专家知识的无人机飞行指挥辅助决策系统及其设计方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102854015B (zh) * | 2012-10-15 | 2014-10-29 | 哈尔滨理工大学 | 一种滚动轴承故障位置及性能退化程度诊断方法 |
CN103033362B (zh) * | 2012-12-31 | 2015-03-25 | 湖南大学 | 一种基于改进多变量预测模型的齿轮故障诊断方法 |
CN103323274B (zh) * | 2013-05-24 | 2015-10-14 | 上海交通大学 | 旋转机械状态监测与故障诊断系统及方法 |
CN103616635B (zh) * | 2013-12-05 | 2017-02-08 | 国家电网公司 | 一种高压断路器机械特性故障诊断装置及方法 |
CN105136454A (zh) * | 2015-10-15 | 2015-12-09 | 上海电机学院 | 一种风电机组齿轮箱故障识别方法 |
-
2018
- 2018-03-28 CN CN201810267031.1A patent/CN108731923B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101025856A (zh) * | 2007-02-08 | 2007-08-29 | 浙江大学 | 一种面向低频振动检测的数据采集装置及方法 |
CN101672723A (zh) * | 2009-10-28 | 2010-03-17 | 北京中能联创风电技术有限公司 | 一种风电机组振动分析故障诊断方法和系统 |
CN102156043A (zh) * | 2010-12-31 | 2011-08-17 | 北京四方继保自动化股份有限公司 | 风力发电机组在线状态监测与故障诊断系统 |
CN104503399A (zh) * | 2014-12-09 | 2015-04-08 | 华电电力科学研究院 | 一种集团级风电机组状态监测及故障诊断平台 |
CN104932519A (zh) * | 2015-05-25 | 2015-09-23 | 北京航空航天大学 | 基于专家知识的无人机飞行指挥辅助决策系统及其设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108731923A (zh) | 2018-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108731923B (zh) | 一种旋转机械设备的故障检测方法及装置 | |
KR101936283B1 (ko) | 기계적 결함에 대한 진단 및 예측 방법 | |
US8630962B2 (en) | Error detection method and its system for early detection of errors in a planar or facilities | |
Wu et al. | Induction machine fault detection using SOM-based RBF neural networks | |
CN111353482A (zh) | 一种基于lstm的疲劳因子隐性异常检测及故障诊断方法 | |
Wang et al. | An evolving fuzzy predictor for industrial applications | |
CN111538311B (zh) | 一种基于数据挖掘的机械设备柔性多状态自适应预警方法及装置 | |
CN109238455A (zh) | 一种基于图论的旋转机械振动信号监测方法及系统 | |
US12007745B2 (en) | Apparatus for predicting equipment damage | |
CN112633098A (zh) | 一种旋转机械故障诊断方法、系统及存储介质 | |
Senanayaka et al. | Towards online bearing fault detection using envelope analysis of vibration signal and decision tree classification algorithm | |
CN112414694A (zh) | 基于多元状态估计技术的设备多级异常状态识别方法及装置 | |
CN112000081A (zh) | 基于多块信息提取和马氏距离的故障监测方法及系统 | |
CN111881594A (zh) | 一种核动力设备的非平稳信号状态监测方法及系统 | |
Heydarzadeh et al. | Gearbox fault diagnosis using power spectral analysis | |
KR20210006832A (ko) | 기계고장 진단 방법 및 장치 | |
CN113505639B (zh) | 一种基于TPE-XGBoost的旋转机械多参数健康状态评估方法 | |
CN118015812A (zh) | 一种设备异常报警方法 | |
McCormick et al. | Application of periodic time-varying autoregressive models to the detection of bearing faults | |
CN113670611A (zh) | 一种轴承早期退化评估方法、系统、介质及电子设备 | |
CN116030955B (zh) | 基于物联网的医疗设备状态监测方法及相关装置 | |
Wang et al. | Explainable machine learning for motor fault diagnosis | |
CN113654651B (zh) | 一种强鲁棒的信号早期退化特征提取及设备运行状态监测方法 | |
Lei et al. | A distance metric learning based health indicator for health prognostics of bearings | |
CN116502043A (zh) | 一种基于孤立森林算法的精轧电机状态分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |