CN113435258A - 转子***异常智能检测方法、***、计算机设备、终端 - Google Patents
转子***异常智能检测方法、***、计算机设备、终端 Download PDFInfo
- Publication number
- CN113435258A CN113435258A CN202110628157.9A CN202110628157A CN113435258A CN 113435258 A CN113435258 A CN 113435258A CN 202110628157 A CN202110628157 A CN 202110628157A CN 113435258 A CN113435258 A CN 113435258A
- Authority
- CN
- China
- Prior art keywords
- rotor system
- data
- discriminator
- frequency
- matrix
- 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.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 46
- 239000011159 matrix material Substances 0.000 claims abstract description 57
- 230000002159 abnormal effect Effects 0.000 claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000012549 training Methods 0.000 claims abstract description 27
- 238000006243 chemical reaction Methods 0.000 claims abstract description 24
- 230000005856 abnormality Effects 0.000 claims abstract description 22
- 238000006073 displacement reaction Methods 0.000 claims abstract description 9
- 238000010606 normalization Methods 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 9
- 238000000354 decomposition reaction Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000011478 gradient descent method Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 230000017105 transposition Effects 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 abstract description 3
- 230000007547 defect Effects 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 230000002194 synthesizing effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000003042 antagnostic effect Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Computational Linguistics (AREA)
- Software Systems (AREA)
- Biophysics (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Signal Processing (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明属于机械技术领域,公开了一种转子***异常智能检测方法、***、计算机设备、终端,所述转子***异常智能检测方法包括:采集转子***支承界面相互垂直方向的电涡流位移信号;利用谐波窗分解提取转子***振动信号的特定频率成分,构建转子***的合成轴心轨迹,将合成轴心轨迹数据转换为一个数值方阵;将所述数值方阵作为生成对抗网络模型的输入,通过正常运行状态下的数据对生成对抗网络进行训练,利用待检测数据输入网络后所得的损失作为异常指标实现转子***运行异常检测。本发明利用全息谱技术,能够更全面地反映转子***转频、分倍频以及高倍频故障相关信息,使得转子***异常检测范围更全面,能够有效避免现有技术中漏检率高的问题。
Description
技术领域
本发明属于机械技术领域,尤其涉及一种转子***异常智能检测方法、***、计算机设备、终端。
背景技术
目前,转子***作为旋转机械中的核心设备,影响着旋转机械的可靠运行。因此,准确对转子***的运行异常情况进行检测对旋转机械的运维管理方面十分重要。随着机器学习和深度学习的快速发展,基于数据驱动的异常检测方法快速发展。神经网络、支持向量回归等都在异常检测领域得到有效的应用。然而,传统的数据驱动的方法需要大量的带标签数据进作为支撑,由于大多数的设备都不允许“带病运行”,这就导致了处于异常运行状态下的数据的缺乏,因此传统的数据驱动方法难以用于转子***异常检测中,无监督方法是解决这个问题最好的方法之一。
戴俊等人在其发表的论文“基于生成对抗网络和自动编码器的机械***异常检测”(仪器仪表学报,2019,40(09):16-26.)中提出了一种基于生成对抗网络的转子***异常检测方法。该方法的步骤是:首先,采集转子***在正常运行状态下的数据;其次,提取这些数据的频谱特征并归一化;然后,利用这些数据训练生成对抗网络与自编码器;最后,将待检测数据输入至网络中,计算异常值实现异常检测。该方法虽然在转子***异常检测方面进行了研究,但是,该方法存在的不足之处是,利用频谱特征作为网络的输入来进行异常检测,只能检测出频域方面的异常,这就导致了容易出现漏检的情况。
通过上述分析,现有技术存在的问题及缺陷为:
(1)传统的数据驱动的方法需要大量的带标签数据进作为支撑,由于大多数的设备都不允许“带病运行”,导致处于异常运行状态下的数据的缺乏。
(2)现有基于生成对抗网络的转子***异常检测方法,利用频谱特征作为网络的输入来进行异常检测,只能检测出频域方面的异常,导致容易出现漏检的情况。
解决以上问题及缺陷的难度为:较难。
解决以上问题及缺陷的意义为:可以充分检测出各种表征在幅值、相位和频率方面的异常运行情况,检测准确率高。
发明内容
针对现有技术存在的问题,本发明提供了一种转子***异常智能检测方法、***、计算机设备、终端,尤其涉及一种基于全息谱与生成对抗网络的转子***异常智能检测方法、***、计算机设备、终端,旨在解决现有转子***运行异常检测容易漏检、准确率不高以及智能化水平不足的问题。
本发明是这样实现的,一种转子***异常智能检测方法,所述转子***异常智能检测方法包括:
首先采集转子***支承界面相互垂直方向的电涡流位移信号;利用谐波窗分解提取转子***振动信号的特定频率成分,构建转子***的合成轴心轨迹,将合成轴心轨迹数据转换为一个数值方阵;将所述数值方阵作为生成对抗网络模型的输入,通过正常运行状态下的数据对生成对抗网络进行训练,利用待检测数据输入网络后所得的损失作为异常指标实现转子***运行异常检测。
进一步,所述转子***异常智能检测方法包括以下步骤:
步骤一,数据采集;
步骤二,构建转子***的合成轴心轨迹;充分考虑转子***运行过程种的幅值、频率和相位信息,可以较为全面地覆盖转子***运行过程种地各种异常情况,可有效避免现有技术中只能检测频域方面异常的缺陷;
步骤三,构建生成对抗网络模型;这是一个无监督模型,该模型地训练不需要依赖带标签的数据,可有效避免现有技术中要依赖大量带标签数据的缺陷;
步骤四,训练生成对抗网络;
步骤五,计算异常值。
进一步,步骤一中,所述数据采集,包括:
采集转子***运行过程中支承截面相互垂直X方向和Y方向的电涡流振动数据,把采样频率记为f,转子的转速记为n,所采集的X方向振动数据记为Dx,所采集的Y方向数据记为Dy,X方向和Y方向数据开始采集与结束采集的时间需要相同,且数据开始采集的时间点为转子***的初相点。
进一步,步骤二中,所述构建转子***的合成轴心轨迹,包括:
(1)采用谐波窗分解分别提取X方向和Y方向转频和高倍频频率成分,记转频为fr,高倍频为fh,fh=kfr,其中k=1,2,3,...,N为正整数,N可根据分析需求自行设定,谐波窗窗口的上限分别为fr+10和fh+10,下限分别为fr-10和fh-10;
(2)采用二进谐波小波变换分别提取X方向和Y方向振动信号分倍频区域内的频率成分,将大于fr的分量滤掉,保留(0,fr)区间内的所有频率成分;
(3)将所得X方向与Y方向振动信号的分倍频、转频和高倍频故障特征频率成分数据分别整合成一个2列m行的矩阵形式,其中,矩阵的第一列为所提取的X方向故障特征频率成分,第二列为Y方向的故障特征频率成分;
(4)将上述矩阵中的元素进行重新排列,得到整合后的转子***合成轴心轨迹数值矩阵,排列方式为按列从原始的矩阵中进行取数,按行填入数值方阵中;
(5)将分倍频、转频和高倍频数值方阵合成为一个三维数值矩阵,在实际运用过程中,也可结合所关注的故障模式,根据实际分析需求选择单一或两个通道的数值矩阵作为后续模型输入。
进一步,步骤三中,所述构建生成对抗网络模型,包括:
(1)构建生成器;
(2)初始化生成器;
(3)构建判别器;
(4)初始化判别器;
其中,所构建的生成器包含9层,依次为转置卷积层、批归一化层、转置卷积层、批归一化层、转置卷积层、批归一化层、转置卷积层、批归一化层、转置卷积层和批归一化层;
所构建的判别器包含9层,依次为卷积层、卷积层、批归一化层、卷积层、批归一化层、卷积层、批归一化层、卷积层和批归一化层。
进一步,步骤四中,所述训练生成对抗网络,包括对所述生成器和判别器进行交替训练:
(1)固定生成器参数,将生成器生成的数据与所述轴心轨迹数据组成一个数据集输入至判别器,采用梯度下降法对判别器进行训练来更新判别器的参数;
(2)固定判别器的参数,将当前判别器的损失的梯度反向传播至生成器,对生成器的参数进行更新;
(3)重复执行步骤(1)-步骤(2),直到网络收敛。
进一步,步骤五中,所述计算异常值,包括:
(1)随机生成噪声向量z;
(2)将噪声向量z输入至生成器,得到生成样本G(z);
(3)根据以下公式计算残差损失:
LR(z)=∑|x-G(z)|;
其中,x为待检测数据;
(4)将生成样本与待检测数据分别输入至判别器,提取中间特征层;
(5)根据以下公式计算判别损失:
LD=∑|f(x)-f(G(z))|;
其中,f(x)表示待检测数据输入至判别器后的中间层特征,f(G(z))表示生成样本输入至判别器后的中间层特征;
(6)根据以下公式计算总损失:
L=(1-λ)LR+λLD;
其中,λ在区间(0,1)中取值;
(7)将总损失L反向传播至噪声z,对噪声z进行更新;
(8)重复上述步骤100次,将第100次所得的总损失作为异常值。
本发明的另一目的在于提供一种应用所述的转子***异常智能检测方法的转子***异常智能检测***,所述转子***异常智能检测***包括:
数据采集模块,用于采集转子***支承界面相互垂直方向电涡流位移信号;
合成轴心轨迹构建模块,用于利用谐波窗分解提取转子***振动信号的特定频率成分,构建转子***的合成轴心轨迹;
对抗网络模型构建模块,用于将合成轴心轨迹数据转换为一个数值方阵,并将所述数值方阵作为生成对抗网络模型的输入;
生成对抗网络训练模块,用于通过正常运行状态下的数据对生成对抗网络进行训练;
异常值计算模块,用于利用待检测数据输入网络后所得的损失作为异常指标实现转子***运行异常检测。
本发明的另一目的在于提供一种计算机设备,所述计算机设备包括存储器和处理器,所述存储器存储有计算机程序,所述计算机程序被所述处理器执行时,使得所述处理器执行如下步骤:
首先采集转子***支承界面相互垂直方向的电涡流位移信号;利用谐波窗分解提取转子***振动信号的特定频率成分,构建转子***的合成轴心轨迹,将合成轴心轨迹数据转换为一个数值方阵;将所述数值方阵作为生成对抗网络模型的输入,通过正常运行状态下的数据对生成对抗网络进行训练,利用待检测数据输入网络后所得的损失作为异常指标实现转子***运行异常检测。
本发明的另一目的在于提供一种信息数据处理终端,所述信息数据处理终端用于实现所述的转子***异常智能检测***。
结合上述的所有技术方案,本发明所具备的优点及积极效果为:本发明提供的转子***异常智能检测方法,具体涉及机械异常检测领域中的一种基于全息谱与生成对抗网络的转子***异常智能检测方法,可用于对转子***的异常运行情况进行智能检测,具有不依赖带标签数据,漏检率低的优点。
本发明利用全息谱技术所衍生的合成轴心轨迹来表征转子***的运行特性,全息谱技术综合利用了振动信号的幅值、频率、相位信息,从而更准确地识别转子运行特性的变化。
本发明采用谐波窗分解法提取若干故障特征频率成分,分别构建了多个合成轴心轨迹,能够更全面地反映转子***转频故障、分倍频故障、高倍频故障相关信息,使得转子***异常检测范围更全面,可有效避免现有技术中漏检率高的问题。
本发明在实现异常检测的过程中,采用生成对抗网络作为异常检测模型,生成对抗网络是一种无监督学习方法,不依赖带标签数据,因此可避免现有技术中需要依赖带标签数据所带来的局限性。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对本发明实施例中所需要使用的附图做简单的介绍,显而易见地,下面所描述的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的转子***异常智能检测方法流程图。
图2是本发明实施例提供的转子***异常智能检测方法原理图。
图3是本发明实施例提供的转子***异常智能检测***结构框图;
图中:1、数据采集模块;2、合成轴心轨迹构建模块;3、对抗网络模型构建模块;4、生成对抗网络训练模块;5、异常值计算模块。
图4是本发明实施例提供的待检测数据的异常值折线图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
针对现有技术存在的问题,本发明提供了一种转子***异常智能检测方法、***、计算机设备、终端,下面结合附图对本发明作详细的描述。
如图1所示,本发明实施例提供的转子***异常智能检测方法包括以下步骤:
S101,采集转子***支承界面相互垂直方向的电涡流位移信号;
S102,利用谐波窗分解提取转子***振动信号的特定频率成分,构建转子***的合成轴心轨迹;
S103,将合成轴心轨迹数据转换为一个数值方阵,并将所述数值方阵作为生成对抗网络模型的输入;
S104,通过正常运行状态下的数据对生成对抗网络进行训练;
S105,利用待检测数据输入网络后所得的损失作为异常指标实现转子***运行异常检测。
本发明实施例提供的转子***异常智能检测方法原理图如图2所示。
如图3所示,本发明实施例提供的转子***异常智能检测***包括:
数据采集模块1,用于采集转子***支承界面相互垂直方向的电涡流位移信号;
合成轴心轨迹构建模块2,用于利用谐波窗分解提取转子***振动信号的特定频率成分,构建转子***的合成轴心轨迹;
对抗网络模型构建模块3,用于将合成轴心轨迹数据转换为一个数值方阵,并将所述数值方阵作为生成对抗网络模型的输入;
生成对抗网络训练模块4,用于通过正常运行状态下的数据对生成对抗网络进行训练;
异常值计算模块5,用于利用待检测数据输入网络后所得的损失作为异常指标实现转子***运行异常检测。
下面结合实施例对本发明的技术方案作进一步描述。
实施例1
参照图2,对本发明实现的步骤作进一步的详细说明。
步骤1,数据采集。
采集转子***运行过程中支承截面相互垂直X方向和Y方向的电涡流振动数据,把采样频率记为f,转子的转速记为n,所采集的X方向振动数据记为Dx,所采集的Y方向数据记为Dy,X方向和Y方向数据开始采集与结束采集的时间需要相同,且数据开始采集的时间点为转子***的初相点。
步骤2,构建转子***的合成轴心轨迹;
第1步,采用谐波窗分解分别提取X方向和Y方向转频和高倍频频率成分,记转频为fr,高倍频为fh,fh=kfr,其中k=1,2,3,...,N为正整数,N可根据分析需求自行设定,谐波窗窗口的上限分别为fr+10和fh+10,下限分别为fr-10和fh-10;
第2步,将上述矩阵中的元素进行重新排列,得到整合后的转子***合成轴心轨迹数值矩阵,排列方式为按列从原始的矩阵中进行取数,按行填入数值方阵中;
第3步,将所得X方向与Y方向的分倍频、转频和高倍频故障特征频率成分数据都分别整合成一个2列m行的矩阵的形式,其中,矩阵的第一列为所提取的X方向故障特征频率成分,第二列为Y方向的故障特征频率成分,得到整合后的转子***合成轴心轨迹数据矩阵;
第4步,上述矩阵中的元素进行重新排列,得到整合后的转子***合成轴心轨迹数值矩阵,排列方式为按列从原始的矩阵中进行取数,按行填入数值方阵中;
第5步,将分倍频、转频和高倍频数值方阵合成为一个三维数值矩阵,在实际运用过程中,也可结合所关注的故障模式,根据实际分析需求选择单一或两个通道的数值矩阵作为后续模型输入;
步骤3,构建生成对抗网络模型:
第1步,构建生成器,该生成器包含9层,依次为转置卷积层、批归一化层、转置卷积层、批归一化层、转置卷积层、批归一化层、转置卷积层、批归一化层、转置卷积层、批归一化层;
第2步,初始化生成器;
第3步,构建判别器,该判别器包含9层,依次为卷积层、卷积层、批归一化层、卷积层、批归一化层、卷积层、批归一化层、卷积层、批归一化层;
第4步,初始化判别器;
步骤4,训练生成对抗网络。
第1步,固定生成器参数,将生成器生成的数据与步骤(2)所得轴心轨迹数据组成一个数据集输入到判别器中,采用梯度下降法对判别器进行训练来更新判别器的参数。
第2步,固定判别器的参数,将当前判别器的损失的梯度反向传播至生成器,对生成器的参数进行更新。
第3步,重复执行上述步骤,直到网络收敛。
步骤5,计算异常值。
第1步,随机生成噪声向量z;
第2步,将噪声向量z输入至生成器,得到生成样本G(z);
第3步,根据以下公式计算残差损失:
LR(z)=∑|x-G(z)|
其中,x为待检测数据;
第4步,将生成样本与待检测数据分别输入至判别器,提取中间特征层;
第5步,根据以下公式计算判别损失:
LD=∑|f(x)-f(G(z))|
其中,f(x)表示待检测数据输入至判别器后的中间层特征,f(G(z))表示生成样本输入至判别器后的中间层特征;
第6步,根据以下公式计算总损失:
L=(1-λ)LR+λLD
第7步,将总损失L反向传播至噪声z,对噪声z进行更新;
第8步,重复上述步骤100次,将第100次所得的总损失作为异常值。
实施例2
本发明实施例中的转子***为鼓风机转子***。
1)数据采集
采集鼓风机转子***运行过程中支承截面相互垂直X方向和Y方向的电涡流振动数据,把采样频率记为f,转子的转速记为n,所采集的X方向振动数据记为Dx,所采集的Y方向数据记为Dy。在鼓风机正常运行阶段共分别采集得到X方向和Y方向数据256个样本,每一个样本包含1024个数据点;在鼓风机故障发展阶段,共采集到X方向和Y方向数据287个样本,每一个样本包含1024个数据点;在鼓风机故障发生阶段共采集到X方向和Y方向数据15个样本,每一个样本包含1024个数据点
2)构建转子***的合成轴心轨迹;
2.1)采用谐波窗分解分别提取X方向和Y方向转频和高倍频频率成分,记转频为fr,高倍频为fh,fh=kfr,其中k=1,2,3,...,N为正整数,N可根据分析需求自行设定,在本实施例中N=2,3,4,谐波窗窗口的上限分别为fr+10和fh+10,下限分别为fr-10和fh-10;
2.2)采用二进谐波小波变换分别提取X方向和Y方向振动信号分倍频区域内的频率成分,将大于fr的分量滤掉,保留(0,fr)区间内的所有频率成分;
2.3)将所得X方向与Y方向振动信号的分倍频、转频和高倍频故障特征频率成分数据分别整合成一个2列m行的矩阵形式,其中,矩阵的第一列为所提取的X方向故障特征频率成分,第二列为Y方向的故障特征频率成分;
2.4)将上述矩阵中的元素进行重新排列,得到整合后的转子***合成轴心轨迹数值矩阵,排列方式为按列从原始的矩阵中进行取数,按行填入数值方阵中;
2.5)将分倍频、转频和高倍频数值方阵合成为一个三维数值矩阵,在实际运用过程中,也可结合所关注的故障模式,根据实际分析需求选择单一或两个通道的数值矩阵作为后续模型输入;
步骤1)所采集的鼓风机正常运行阶段,故障发展阶段和故障发生阶段的数据都采用上述步骤进行处理,本例中关注的是转频以及高倍频故障,因此不知用分倍频数据,高倍频频率成分包含了二倍转频、三倍转频、四倍转频,转频和高倍频频率成分所得每一个合成轴心轨迹数值方阵的大小为64*64,将其合成后所得每一个数值矩阵的大小为2*64*64。
3)构建生成对抗网络模型:
3.1)构建生成器,该生成器包含9层,依次为转置卷积层、批归一化层、转置卷积层、批归一化层、转置卷积层、批归一化层、转置卷积层、批归一化层、转置卷积层、批归一化层;
3.2)初始化生成器;
3.3)构建判别器,该判别器包含9层,依次为卷积层、卷积层、批归一化层、卷积层、批归一化层、卷积层、批归一化层、卷积层、批归一化层;
3.4)初始化判别器;
4)训练生成对抗网络
4.1)固定生成器参数,将生成器生成的数据与步骤(2)所得轴心轨迹数据组成一个数据集输入到判别器中,采用梯度下降法对判别器进行训练来更新判别器的参数。
第2步,固定判别器的参数,将当前判别器的损失的梯度反向传播至生成器,对生成器的参数进行更新。
第3步,重复执行上述步骤,直到网络收敛。
5)计算异常值
5.1)随机生成噪声向量z;
5.2)将噪声向量z输入至生成器,得到生成样本G(z);
5.3)根据以下公式计算残差损失:
LR(z)=∑|x-G(z)|
其中,x为待检测数据,本实施例中使用的待检测数据为步骤2)处理后所得的所有的数值矩阵;
5.4)将生成样本与待检测数据分别输入至判别器,提取中间特征层;
5.5)根据以下公式计算判别损失:
LD=∑|f(x)-f(G(z))|
其中,f(x)表示待检测数据输入至判别器后的中间层特征,f(G(z))表示生成样本输入至判别器后的中间层特征;
5.6)根据以下公式计算总损失:
L=(1-λ)LR+λLD
5.7)将总损失L反向传播至噪声z,对噪声z进行更新;
5.8)重复上述步骤100次,将第100次所得的总损失作为异常值。
本发明实施例中将所有待检测数据的异常值都计算完毕后,将正常运行阶段、故障发展阶段和故障发生阶段这三个阶段共595个异常值绘制成图4所示的3条折线。图4中的横坐标表示数据序号,纵坐标表示异常值。
在上述实施例中,可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用全部或部分地以计算机程序产品的形式实现,所述计算机程序产品包括一个或多个计算机指令。在计算机上加载或执行所述计算机程序指令时,全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中,或者从一个计算机可读存储介质向另一个计算机可读存储介质传输,例如,所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL)或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输)。所述计算机可读取存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质,(例如,软盘、硬盘、磁带)、光介质(例如,DVD)、或者半导体介质(例如固态硬盘SolidState Disk(SSD))等。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种转子***异常智能检测方法,其特征在于,所述转子***异常智能检测方法包括:采集转子***支承界面相互垂直方向的电涡流位移信号;利用谐波窗分解提取转子***振动信号的特定频率成分,构建转子***的合成轴心轨迹,将合成轴心轨迹数据转换为一个数值方阵;将所述数值方阵作为生成对抗网络模型的输入,通过正常运行状态下的数据对生成对抗网络进行训练,利用待检测数据输入网络后所得的损失作为异常指标实现转子***运行异常检测。
2.如权利要求1所述转子***异常智能检测方法,其特征在于,所述转子***异常智能检测方法包括以下步骤:
步骤一,数据采集;
步骤二,构建转子***的合成轴心轨迹;
步骤三,构建生成对抗网络模型;
步骤四,训练生成对抗网络;
步骤五,计算异常值。
3.如权利要求2所述转子***异常智能检测方法,其特征在于,步骤一中,所述数据采集,包括:
采集转子***运行过程中支承截面相互垂直X方向和Y方向的电涡流振动数据,把采样频率记为f,转子的转速记为n,所采集的X方向振动数据记为Dx,所采集的Y方向数据记为Dy,X方向和Y方向数据开始采集与结束采集的时间需要相同,且数据开始采集的时间点为转子***的初相点。
4.如权利要求2所述转子***异常智能检测方法,其特征在于,步骤二中,所述构建转子***的合成轴心轨迹,包括:
(1)采用谐波窗分解分别提取X方向和Y方向转频和高倍频频率成分,记转频为fr,高倍频为fh,fh=kfr,其中k=1,2,3,...,N为正整数,N可根据分析需求自行设定,谐波窗窗口的上限分别为fr+10和fh+10,下限分别为fr-10和fh-10;
(2)采用二进谐波小波变换分别提取X方向和Y方向振动信号分倍频区域内的频率成分,将大于fr的分量滤掉,保留(0,fr)区间内的所有频率成分;
(3)将所得X方向与Y方向振动信号的分倍频、转频和高倍频故障特征频率成分数据分别整合成一个2列m行的矩阵形式,其中,矩阵的第一列为所提取的X方向故障特征频率成分,第二列为Y方向的故障特征频率成分;
(4)将上述矩阵中的元素进行重新排列,得到整合后的转子***合成轴心轨迹数值矩阵,排列方式为按列从原始的矩阵中进行取数,按行填入数值方阵中;
(5)将分倍频、转频和高倍频数值方阵合成为一个三维数值矩阵,在实际运用过程中,也可结合所关注的故障模式,根据实际分析需求选择单一或两个通道的数值矩阵作为后续模型输入。
5.如权利要求2所述转子***异常智能检测方法,其特征在于,步骤三中,所述构建生成对抗网络模型,包括:
(1)构建生成器;
(2)初始化生成器;
(3)构建判别器;
(4)初始化判别器;
其中,所构建的生成器包含9层,依次为转置卷积层、批归一化层、转置卷积层、批归一化层、转置卷积层、批归一化层、转置卷积层、批归一化层、转置卷积层和批归一化层;
所构建的判别器包含9层,依次为卷积层、卷积层、批归一化层、卷积层、批归一化层、卷积层、批归一化层、卷积层和批归一化层。
6.如权利要求2所述转子***异常智能检测方法,其特征在于,步骤四中,所述训练生成对抗网络,包括对所述生成器和判别器进行交替训练:
(1)固定生成器参数,将生成器生成的数据与所述轴心轨迹数据组成一个数据集输入至判别器,采用梯度下降法对判别器进行训练来更新判别器的参数;
(2)固定判别器的参数,将当前判别器的损失的梯度反向传播至生成器,对生成器的参数进行更新;
(3)重复执行步骤(1)-步骤(2),直到网络收敛。
7.如权利要求2所述转子***异常智能检测方法,其特征在于,步骤五中,所述计算异常值,包括:
(1)随机生成噪声向量z;
(2)将噪声向量z输入至生成器,得到生成样本G(z);
(3)根据以下公式计算残差损失:
LR(z)=∑|x-G(z)|;
其中,x为待检测数据;
(4)将生成样本与待检测数据分别输入至判别器,提取中间特征层;
(5)根据以下公式计算判别损失:
LD=∑|f(x)-f(G(z))|;
其中,f(x)表示待检测数据输入至判别器后的中间层特征,f(G(z))表示生成样本输入至判别器后的中间层特征;
(6)根据以下公式计算总损失:
L=(1-λ)LR+λLD;
其中,λ在区间(0,1)中取值;
(7)将总损失L反向传播至噪声z,对噪声z进行更新;
(8)重复上述步骤100次,将第100次所得的总损失作为异常值。
8.一种实施权利要求1~7任意一项所述转子***异常智能检测方法的转子***异常智能检测***,其特征在于,所述转子***异常智能检测***包括:
数据采集模块,用于采集转子***支承界面相互垂直方向电涡流位移信号;
合成轴心轨迹构建模块,用于利用谐波窗分解提取转子***振动信号的特定频率成分,构建转子***的合成轴心轨迹;
对抗网络模型构建模块,用于将合成轴心轨迹数据转换为一个数值方阵,并将所述数值方阵作为生成对抗网络模型的输入;
生成对抗网络训练模块,用于通过正常运行状态下的数据对生成对抗网络进行训练;
异常值计算模块,用于利用待检测数据输入网络后所得的损失作为异常指标实现转子***运行异常检测。
9.一种计算机设备,其特征在于,所述计算机设备包括存储器和处理器,所述存储器存储有计算机程序,所述计算机程序被所述处理器执行时,使得所述处理器执行如下步骤:
首先采集转子***支承界面相互垂直方向的电涡流位移信号;利用谐波窗分解提取转子***振动信号的特定频率成分,构建转子***的合成轴心轨迹,将合成轴心轨迹数据转换为一个数值方阵;将所述数值方阵作为生成对抗网络模型的输入,通过正常运行状态下的数据对生成对抗网络进行训练,利用待检测数据输入网络后所得的损失作为异常指标实现转子***运行异常检测。
10.一种信息数据处理终端,其特征在于,所述信息数据处理终端用于实现如权利要求8所述转子***异常智能检测***。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110628157.9A CN113435258B (zh) | 2021-06-06 | 2021-06-06 | 转子***异常智能检测方法、***、计算机设备、终端 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110628157.9A CN113435258B (zh) | 2021-06-06 | 2021-06-06 | 转子***异常智能检测方法、***、计算机设备、终端 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113435258A true CN113435258A (zh) | 2021-09-24 |
CN113435258B CN113435258B (zh) | 2024-04-23 |
Family
ID=77803776
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110628157.9A Active CN113435258B (zh) | 2021-06-06 | 2021-06-06 | 转子***异常智能检测方法、***、计算机设备、终端 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113435258B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110823576A (zh) * | 2019-11-18 | 2020-02-21 | 苏州大学 | 基于生成对抗网络的机械异常检测方法 |
WO2020044495A1 (en) * | 2018-08-30 | 2020-03-05 | Nec Corporation | Anomaly detection apparatus, method and computer-readable medium |
US20200202221A1 (en) * | 2018-12-20 | 2020-06-25 | Shandong University Of Science And Technology | Fault detection method and system based on generative adversarial network and computer program |
CN111337243A (zh) * | 2020-02-27 | 2020-06-26 | 上海电力大学 | 基于acgan的风电机组行星轮齿轮箱故障诊断方法 |
CN112200244A (zh) * | 2020-10-09 | 2021-01-08 | 西安交通大学 | 一种基于递阶对抗训练的航天发动机异常智能检测方法 |
US20210049452A1 (en) * | 2019-08-15 | 2021-02-18 | Intuit Inc. | Convolutional recurrent generative adversarial network for anomaly detection |
KR102229638B1 (ko) * | 2019-10-31 | 2021-03-18 | 한국철도기술연구원 | 고장 예측 장치 및 방법 |
CN112561383A (zh) * | 2020-12-24 | 2021-03-26 | 航天科工网络信息发展有限公司 | 一种基于生成对抗网络的实时异常检测方法 |
WO2021064833A1 (ja) * | 2019-09-30 | 2021-04-08 | 日本電気株式会社 | 異常検知装置、異常検知方法、及び、記録媒体 |
CN112766342A (zh) * | 2021-01-12 | 2021-05-07 | 安徽容知日新科技股份有限公司 | 一种电气设备的异常检测方法 |
-
2021
- 2021-06-06 CN CN202110628157.9A patent/CN113435258B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020044495A1 (en) * | 2018-08-30 | 2020-03-05 | Nec Corporation | Anomaly detection apparatus, method and computer-readable medium |
US20200202221A1 (en) * | 2018-12-20 | 2020-06-25 | Shandong University Of Science And Technology | Fault detection method and system based on generative adversarial network and computer program |
US20210049452A1 (en) * | 2019-08-15 | 2021-02-18 | Intuit Inc. | Convolutional recurrent generative adversarial network for anomaly detection |
WO2021064833A1 (ja) * | 2019-09-30 | 2021-04-08 | 日本電気株式会社 | 異常検知装置、異常検知方法、及び、記録媒体 |
KR102229638B1 (ko) * | 2019-10-31 | 2021-03-18 | 한국철도기술연구원 | 고장 예측 장치 및 방법 |
CN110823576A (zh) * | 2019-11-18 | 2020-02-21 | 苏州大学 | 基于生成对抗网络的机械异常检测方法 |
CN111337243A (zh) * | 2020-02-27 | 2020-06-26 | 上海电力大学 | 基于acgan的风电机组行星轮齿轮箱故障诊断方法 |
CN112200244A (zh) * | 2020-10-09 | 2021-01-08 | 西安交通大学 | 一种基于递阶对抗训练的航天发动机异常智能检测方法 |
CN112561383A (zh) * | 2020-12-24 | 2021-03-26 | 航天科工网络信息发展有限公司 | 一种基于生成对抗网络的实时异常检测方法 |
CN112766342A (zh) * | 2021-01-12 | 2021-05-07 | 安徽容知日新科技股份有限公司 | 一种电气设备的异常检测方法 |
Non-Patent Citations (2)
Title |
---|
包萍 等: "不均衡数据集下基于生成对抗网络的改进深度模型故障识别研究", 电子测量与仪器学报, no. 03, 15 March 2019 (2019-03-15) * |
尹海: "离心式压缩机轴心轨迹自动识别及其故障诊断***", 中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑, 15 December 2013 (2013-12-15) * |
Also Published As
Publication number | Publication date |
---|---|
CN113435258B (zh) | 2024-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112785016B (zh) | 基于机器学习的新能源汽车保养维护与故障监测诊断方法 | |
US20170356936A1 (en) | Enhanced system and method for conducting pca analysis on data signals | |
Isham et al. | Variational mode decomposition: mode determination method for rotating machinery diagnosis | |
CN111504676A (zh) | 基于多源监控数据融合的设备故障诊断方法、装置及*** | |
Kaya et al. | Classification of bearing vibration speeds under 1D-LBP based on eight local directional filters | |
Ma et al. | Early fault diagnosis of rotating machinery based on composite zoom permutation entropy | |
CN112084375A (zh) | 一种车辆的故障诊断方法、装置、终端设备及存储介质 | |
Ince et al. | Early bearing fault diagnosis of rotating machinery by 1D self-organized operational neural networks | |
CN111881594B (zh) | 一种核动力设备的非平稳信号状态监测方法及*** | |
CN111291918B (zh) | 平稳子空间外源矢量自回归的旋转机械退化趋势预测方法 | |
CN114819315A (zh) | 基于多参数融合健康因子和时间卷积神经网络的轴承退化趋势预测方法 | |
Xin et al. | Intelligent fault diagnosis method for rotating machinery based on vibration signal analysis and hybrid multi‐object deep CNN | |
Lee et al. | An explainable deep learning-based prognostic model for rotating machinery | |
CN116108346A (zh) | 一种基于生成特征重放的轴承增量故障诊断终身学习方法 | |
Wang et al. | Fault diagnosis of industrial robots based on multi-sensor information fusion and 1D convolutional neural network | |
CN115600136A (zh) | 一种基于多传感器的高压套管故障诊断方法、***及介质 | |
CN115481657A (zh) | 一种基于电信号的风力发电机组通信滑环故障诊断方法 | |
CN114462446A (zh) | 基于振动信号的滚动轴承故障诊断方法及电子设备 | |
Wang et al. | Fast Spectral Correlation Based on Sparse Representation Self‐Learning Dictionary and Its Application in Fault Diagnosis of Rotating Machinery | |
CN113435258A (zh) | 转子***异常智能检测方法、***、计算机设备、终端 | |
CN116977708A (zh) | 一种基于自适应聚合可视图的轴承智能诊断方法及*** | |
CN115356599A (zh) | 一种多模态城市电网故障诊断方法及*** | |
CN114692694B (zh) | 一种基于特征融合和集成聚类的设备故障诊断方法 | |
Ghosh | Comparative DNN model analysis for detection of various types of optical noise | |
CN114186593A (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 |