CN108021876B - 一种基于参数辨识的压气机故障诊断方法 - Google Patents

一种基于参数辨识的压气机故障诊断方法 Download PDF

Info

Publication number
CN108021876B
CN108021876B CN201711210934.8A CN201711210934A CN108021876B CN 108021876 B CN108021876 B CN 108021876B CN 201711210934 A CN201711210934 A CN 201711210934A CN 108021876 B CN108021876 B CN 108021876B
Authority
CN
China
Prior art keywords
components
compressor fault
component
correlation
signal
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
Application number
CN201711210934.8A
Other languages
English (en)
Other versions
CN108021876A (zh
Inventor
徐搏超
阮圣奇
吴仲
王松浩
许昊煜
李强
胡中强
任磊
蒋怀锋
陈开峰
邵飞
徐钟宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Datang Corp Science and Technology Research Institute Co Ltd
Original Assignee
China Datang Corp Science and Technology Research Institute Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Datang Corp Science and Technology Research Institute Co Ltd filed Critical China Datang Corp Science and Technology Research Institute Co Ltd
Priority to CN201711210934.8A priority Critical patent/CN108021876B/zh
Publication of CN108021876A publication Critical patent/CN108021876A/zh
Application granted granted Critical
Publication of CN108021876B publication Critical patent/CN108021876B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M99/00—Subject matter not provided for in other groups of this subclass
    • G01M99/005—Testing of complete machines, e.g. washing-machines or mobile phones
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/23—Clustering techniques
    • G06F18/231—Hierarchical techniques, i.e. dividing or merging pattern sets so as to obtain a dendrogram

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Computation (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Signal Processing (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

本发明公开了一种基于参数辨识的压气机故障诊断方法，包括以下步骤：首先对原始压气机故障信号进行信号分解，得到一系列的分量；然后分别求出每一个分量与原始压气机故障信号间的K‑L散度、互信息和相关系数值作为该分量综合指标的三个维度值；再将综合指标作为各分量的特征向量，将所有分量的特征向量组成一个集合，对集合内的元素进行分层聚类，使真实分量与虚假分量聚为两类；最后消除虚假分量。本发明的优点在于：基于本发明，可以降低信号分解方法的分析误差，构建故障诊断特征阈值库。

Description

一种基于参数辨识的压气机故障诊断方法
技术领域
本发明涉及旋转机械振动响应的时频域分析技术领域，尤其涉及一种基于参数辨识的压气机故障诊断方法。
背景技术
燃机压气机气动特性恶化机理较为复杂，可能引起故障的原因众多。在生产实际中为了精确定位故障激振源，需要对故障参数进行深入的处理分析。在提取故障参数特征，构建实时诊断平台的过程中，信号处理方法的误差往往对诊断结果有着极为重大的影响。
集成经验模态分解(EEMD)，希尔伯特振动分解(HVD)等信号分解方法在各领域得到了广泛应用，但这些方法在分解过程中存在着虚假分量的问题。为了解决虚假分量问题，使HVD等常用信号分解方法具有更高的工程实用价值，需要引入合适的指标对虚假分量进行鉴别。针对信号分解过程中的虚假分量问题，许多学者进行了相关研究并提出了一些鉴别指标，包括K-L散度、互信息和相关系数等。这些鉴别指标需要通过预先设定阈值区分真实和虚假分量，不能够实现真假分量的自动识别；同时，设定的阈值只是对试验样本规律的总结，缺乏理论支持，容易受到试验样本集片面性和单一指标误诊率的影响，难以推广。上述指标在使用过程中往往并不能全面反映分量与原始信号的内在关联。
发明内容
本发明所要解决的技术问题是提供一种可以实现真假分量的自动识别的基于参数辨识的压气机故障诊断方法。
本发明是通过以下技术方案解决上述技术问题的：一种基于参数辨识的压气机故障诊断方法，包括以下步骤：首先对原始压气机故障信号进行信号分解，得到一系列的分量；然后分别求出每一个分量与原始压气机故障信号间的K-L散度、互信息和相关系数值作为该分量综合指标的三个维度值；再将综合指标作为各分量的特征向量，将所有分量的特征向量组成一个集合，对集合内的元素进行分层聚类，使真实分量与虚假分量聚为两类；最后消除虚假分量。
作为优化的技术方案，
求出分量与原始压气机故障信号间的K-L散度的方法包括以下步骤：
设X＝{x1，x2，…，xn}，Y＝{y1，y2，…，yn}，X和Y分别为其中一个分量和原始压气机故障信号；
首先通过非参数估计法求出X概率分布的核密度估计p(x)和Y概率分布的核密度估计q(x)，
Figure BDA0001484579260000023
Figure BDA0001484579260000024
其中，平滑参数h是给定的正数，K(·)是核函数；
然后求出X和Y的K-L距离δ(p，q)和δ(q，p)，
Figure BDA0001484579260000021
Figure BDA0001484579260000022
最后求出X和Y的K-L散度值D(p,q)，
D(p,q)＝δ(p,q)+δ(q,p)。
本发明的优点在于：通过构造一种综合指标向量，结合多种误差方法的优势，能够更准确地识别虚假分量，更真实的反映故障特征；基于本发明，可以降低信号分解方法的分析误差，构建故障诊断特征阈值库，为电厂后续的处理提供了参考依据。
附图说明
图1是本发明基于参数辨识的压气机故障诊断方法的流程图。
具体实施方式
虚假分量是由于分解算法的计算误差产生的，虚假分量与在原始信号中真实存在的分量相比，二者与原始信号的相关度差异很大。通过计算各分量与原始信号之间的相似程度，可以将真实分量与虚假分量自动划分为两类。
一种基于参数辨识的压气机故障诊断方法，包括以下步骤：首先对原始压气机故障信号进行信号分解，得到一系列的分量；然后分别求出每一个分量与原始压气机故障信号间的K-L散度、互信息和相关系数值作为该分量综合指标的三个维度值；再将综合指标作为各分量的特征向量，将所有分量的特征向量组成一个集合，对集合内的元素进行分层聚类，使真实分量与虚假分量聚为两类；最后消除虚假分量。
综合指标具有三个维度，记作F，F＝[K-L散度，互信息，相关系数]。
1、K-L散度
K-L散度又称为相对熵，是在信息论中用于衡量两种概率分布差异的方法，K-L散度值越大，表示两种分布的差别越大，反之则表示差别越小。
设p(x)、q(x)表示两种概率分布，则K-L距离δ(p,q)为：
Figure BDA0001484579260000031
由于K-L距离不具有对称性，不满足真正意义上距离的概念，所以不适合作为衡量p(x)、q(x)差异的定量指标。
故选用K-L散度值D(p，q)作为评价指标：
D(p，q)＝δ(p，q)+δ(q，p)。
对于两个信号X＝{x1，x2，…，xn}和Y＝{y1，y2，…，yn}，K-L散度值的计算方法如下：
首先通过非参数估计法求出X概率分布的核密度估计p(x)和Y概率分布的核密度估计q(x)：
Figure BDA0001484579260000041
Figure BDA0001484579260000042
其中，平滑参数h是给定的正数，K(·)是核函数，最常用的核函数是高斯核函数，即：
Figure BDA0001484579260000043
然后求出X和Y的K-L距离δ(p，q)和δ(q，p)，
Figure BDA0001484579260000044
Figure BDA0001484579260000045
最后求出X和Y的K-L散度值D(p,q)，
D(p,q)＝δ(p,q)+δ(q,p)。
2、互信息
在信息论中，假定变量中蕴含了一类样本集中的某些信息，这个变量就被称为这个样本集的有信息变量。互信息(mutual information,MI)定量表示任意两个随机变量之间含有相同信息量的多少，可以从非线性的角度对随机变量间的相关程度进行度量。
对于两个随机变量X和Y，用下式计算X和Y之间的互信息：
Figure BDA0001484579260000046
当X和Y为离散变量时，它们之间的互信息为：
Figure BDA0001484579260000051
当变量之间相互独立时，它们之间的互信息为零。
互信息与信息熵之间的关系如下：
I(X；Y)＝H(X)-H(X|Y)
＝H(Y)-H(Y|X)，
＝H(X)+H(Y)-H(XY)
其中H(X)和H(Y)分别是X和Y的信息熵，H(XY)是X和Y的联合信息熵，
Figure BDA0001484579260000052
对于离散变量，计算如下：
Figure BDA0001484579260000053
3、相关系数
互相关函数可以用来描述两个信号之间的相关性，如果一个信号x(t)不变，另一个信号y(t)延迟一段时间τ，互相关分析描述了它们之间的相关程度。这种相关程度是随时间τ的变化而变化，所以也称为互相关函数。
对于x(t)和y(t)两个平稳信号，它们的互相关函数为：
Figure BDA0001484579260000054
由于工业生产中传感器采集到的信号都是离散的，如果两个离散信号x(t)和y(t)的均值为零，则它们相关系数的求解公式为：
Figure BDA0001484579260000055
式中n为信号的长度。
相关系数的绝对值在0和1之间。ρxy越接近0表示两者之间的相关性越小，越接近1表示两者之间的相关性越大。
4、聚类分析
聚类的严格定义为：对于确定数据集合X＝{X1,X2,...Xn}，遵循数据对象间的相似性将样本集进行无监督划分，分成m个簇{Cl,C2,...,Cm}，这个过程称作聚类，其中
Figure BDA0001484579260000061
聚类分析是在数据集中自动搜寻不同簇的隐藏的数据模型。聚类不同于分类，不需要预先设定类别，也不需要对样本个体进行类别标记。聚类分析适用于数据集合内部蕴含规律但数据类别未知的集合元素划分过程。基于综合指标的特征向量可以较为准确地反映各分量和原始信号之间的相似程度，故对特征向量进行聚类能够实现真实与虚假分量的自动分类。
本发明基于参数辨识的压气机故障诊断方法通过聚类分析对HVD等常用信号分解方法信号分解过程中的虚假分量进行识别。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (2)

1.一种基于参数辨识的压气机故障诊断方法，其特征在于：包括以下步骤：首先对原始压气机故障信号进行信号分解，得到一系列的分量；然后分别求出每一个分量与原始压气机故障信号间的K-L散度、互信息和相关系数值作为该分量综合指标的三个维度值；再将综合指标作为各分量的特征向量，将所有分量的特征向量组成一个集合，对集合内的元素进行分层聚类，使真实分量与虚假分量聚为两类；最后消除虚假分量；
求出分量与原始压气机故障信号间的K-L散度的方法包括以下步骤：
设X＝{x1，x2，…，xn}，Y＝{y1，y2，…，yn}，X和Y分别为其中一个分量和原始压气机故障信号；
首先通过非参数估计法求出X概率分布的核密度估计p(x)和Y概率分布的核密度估计q(x)，
Figure FDA0002564712990000011
Figure FDA0002564712990000012
其中，平滑参数h是给定的正数，K(·)是核函数；
然后求出X和Y的K-L距离δ(p，q)和δ(q，p)，
Figure FDA0002564712990000013
Figure FDA0002564712990000014
最后求出X和Y的K-L散度值D(p,q)，
D(p,q)＝δ(p,q)+δ(q,p)；
对于两个随机变量X和Y，用下式计算X和Y之间的互信息：
Figure FDA0002564712990000021
当X和Y为离散变量时，它们之间的互信息为：
Figure FDA0002564712990000022
当变量之间相互独立时，它们之间的互信息为零；
互信息与信息熵之间的关系如下：
Figure FDA0002564712990000023
其中H(X)和H(Y)分别是X和Y的信息熵，H(XY)是X和Y的联合信息熵，
Figure FDA0002564712990000024
对于离散变量，计算如下：
Figure FDA0002564712990000025
求出分量与原始压气机故障信号间的相关系数的方法包括以下步骤：
对于x(t)和y(t)两个平稳信号，它们的互相关函数为：
Figure FDA0002564712990000026
由于工业生产中传感器采集到的信号都是离散的，如果两个离散信号x(t)和y(t)的均值为零，则它们相关系数的求解公式为：
Figure FDA0002564712990000027
式中n为信号的长度，
相关系数的绝对值在0和1之间，ρxy越接近0表示两者之间的相关性越小，越接近1表示两者之间的相关性越大。
2.如权利要求1所述的基于参数辨识的压气机故障诊断方法，其特征在于：对特征向量集合内的元素进行分层聚类的方法包括以下步骤：
对于确定数据集合X＝{X1,X2,...Xn}，遵循数据对象间的相似性将样本集进行无监督划分，分成m个簇{Cl,C2,...,Cm}，其中
Figure FDA0002564712990000031
CN201711210934.8A 2017-11-28 2017-11-28 一种基于参数辨识的压气机故障诊断方法 Active CN108021876B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711210934.8A CN108021876B (zh) 2017-11-28 2017-11-28 一种基于参数辨识的压气机故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711210934.8A CN108021876B (zh) 2017-11-28 2017-11-28 一种基于参数辨识的压气机故障诊断方法

Publications (2)

Publication Number Publication Date
CN108021876A CN108021876A (zh) 2018-05-11
CN108021876B true CN108021876B (zh) 2020-09-01

Family

ID=62077510

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711210934.8A Active CN108021876B (zh) 2017-11-28 2017-11-28 一种基于参数辨识的压气机故障诊断方法

Country Status (1)

Country Link
CN (1) CN108021876B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109587652B (zh) * 2018-12-31 2022-03-08 宁波工程学院 一种农业物联网故障诊断方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7552005B2 (en) * 2004-03-16 2009-06-23 Honeywell International Inc. Method for fault diagnosis of a turbine engine
CN103868692B (zh) * 2014-03-18 2016-02-10 电子科技大学 基于核密度估计和k-l散度的旋转机械故障诊断方法

Also Published As

Publication number Publication date
CN108021876A (zh) 2018-05-11

Similar Documents

Publication Publication Date Title
CN107451102B (zh) 一种基于改进自训练算法的半监督高斯过程回归软测量建模对脱丁烷塔底丁烷浓度进行预测的方法
Wang et al. An improved K-Means clustering algorithm
CN110632572A (zh) 基于无意调相特性的雷达辐射源个体识别方法及装置
CN110352349A (zh) 异常音检测装置、异常度计算装置、异常音生成装置、异常音检测学习装置、异常信号检测装置、异常信号检测学习装置、及它们的方法以及程序
CN109543720B (zh) 一种基于对抗生成网络的晶圆图缺陷模式识别方法
CN106257488B (zh) 一种基于邻域特征空间鉴别分析的雷达目标识别方法
CN110632546B (zh) 基于全网域证据集的电子式互感器可信度评估方法和装置
Yang Monitoring and diagnosing of mean shifts in multivariate manufacturing processes using two-level selective ensemble of learning vector quantization neural networks
Arif et al. Cascade quality prediction method using multiple PCA+ ID3 for multi-stage manufacturing system
CN113642666B (zh) 一种基于样本扩充和筛选的主动增强软测量方法
Zhao et al. A novel deep fuzzy clustering neural network model and its application in rolling bearing fault recognition
Wang et al. Production quality prediction of multistage manufacturing systems using multi-task joint deep learning
CN108021876B (zh) 一种基于参数辨识的压气机故障诊断方法
Buabeng et al. Hybrid intelligent predictive maintenance model for multiclass fault classification
Liu et al. A roughness measurement method based on genetic algorithm and neural network for microheterogeneous surface in deep-hole parts
Hassanvand et al. A robust control chart for simple linear profiles in two‐stage processes
CN111160464B (zh) 基于多隐层加权动态模型的工业高阶动态过程软测量方法
Wang et al. FCM algorithm and index CS for the signal sorting of radiant points
Lennox et al. Deep metric learning for proteomics
CN114021275A (zh) 一种基于深度卷积模糊系统的滚动轴承故障诊断方法
CN113220960A (zh) 一种基于自主学习的不平衡时间序列数据分类方法
Azmer et al. Comparative analysis of classification techniques for leaves and land cover texture.
Corchado et al. Optimizing the operating conditions in a high precision industrial process using soft computing techniques
Masti et al. Learning virtual sensors for estimating the scheduling signal of parameter-varying systems
Laosiritaworn et al. Classification techniques for control chart pattern recognition: A case of metal frame for actuator production

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