CN111368885B - 一种航空发动机气路故障诊断方法 - Google Patents

一种航空发动机气路故障诊断方法 Download PDF

Info

Publication number
CN111368885B
CN111368885B CN202010112320.1A CN202010112320A CN111368885B CN 111368885 B CN111368885 B CN 111368885B CN 202010112320 A CN202010112320 A CN 202010112320A CN 111368885 B CN111368885 B CN 111368885B
Authority
CN
China
Prior art keywords
model
sample
evidence
data
training
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
CN202010112320.1A
Other languages
English (en)
Other versions
CN111368885A (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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN202010112320.1A priority Critical patent/CN111368885B/zh
Publication of CN111368885A publication Critical patent/CN111368885A/zh
Application granted granted Critical
Publication of CN111368885B publication Critical patent/CN111368885B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/25Fusion techniques
    • G06F18/254Fusion techniques of classification results, e.g. of results related to same input data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Software Systems (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Mathematical Physics (AREA)
  • Computational Linguistics (AREA)
  • Health & Medical Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种基于深度学习和信息融合的航空发动机气路故障诊断方法,属于动力机械故障诊断领域。首先,对数据进行预处理;其次,利用处理后的样本数据将卷积神经网络模型和前馈神经网络分别训练至最佳,继而得到两个模型对样本数据所属类别的判定结果;再次,在证据体的基本概率分布的构造阶段,两种模型的判别结果被视为两个证据体,该阶段将证据体中样本所属各类别评分转换为样本的基本概率分布,同时计算样本不确定度;最后,在证据体的合成与决策阶段,根据D‑S证据理论的合成规则计算两个证据体合成后的概率分布,根据决策规则给出每条样本所属类别。本发明融合了两种深度神经网络从不同角度学习样本数据后给出的决策结果,获得一种精度较高的航发故障诊断模型,同时可以克服非深度学习模型对不同来源的航发大数据判别能力不稳定、鲁棒性低的问题。

Description

一种航空发动机气路故障诊断方法
技术领域
本发明属于动力机械故障诊断领域,涉及一种基于深度学习和信息融合的航空发动机气路故障诊断方法,用于航空发动机或燃气轮机实时故障检测与故障分类。
背景技术
航空发动机(以下简称航发)被誉为“工业皇冠上的明珠”,其内部系统结构具有极高的复杂性。正因如此,研究人员难以利用航发领域的机械、流体及热力学方面的专业知识对其进行故障诊断。现有的航发故障诊断技术大多是通过大数据与人工智能的方法来实现的。故障检测与分类属于人工智能领域的模式识别问题。
在航发故障诊断这类模式识别的实际应用领域中还存在一些共性问题:航发大数据的数据格式不规范。从不同航空企业、研究单位获取的数据在维度、连续性、数量级及质量上有很大差异。所以,目前很多故障诊断模型的应用场景难以泛化,可能只针对某些特定采样方式下的航发参数数据集才有很高的精度。若将其推广应用,效果可能会下降很多。
当前应用信息技术方法进行航空发动机故障诊断的方法可归为以下几类:
基于信号处理的诊断技术,如应用小波变换、傅里叶变换、卡尔曼滤波等方法,对现场采集的连续的信号波形进行特征提取,分析得出结果。此类方法需要用一定时间间隔内的一段波形进行分析,才能得出数据特征。故只适用于连续采样的数据。
基于位置和距离进行分类的算法,如支持向量机。这类方法的思想是将样本的各维数据视为多维向量的坐标。不适用于较多分类的应用和数值区分度不高的样本数据。
基于专家经验与数学模型相融合的诊断方法。专家经验通常表现为依据样本属性或特征的阈值进行分支的决策树,每一叶节点对应一个故障类别的概率值。此类方法适合低维样本数据,样本维数过高时或多分类时,专家经验概率组合会十分复杂。此外,专家经验概率主观性较强,其正确与否将直接影响判别结果。
基于神经网络的诊断模型,如前馈神经网络、径向基神经网络、自组织竞争网络等基本的深度学习模型。这些模型通过拟合样本数据与其所属分类间的复杂函数关系,来做出正确的分类判断。相比于以上方法,神经网络由于有很强的泛化能力和一定的容错能力,因此神经网络在应用时,其对数据质量和规格要求相对较低,适用范围更广。此外还有些方法是利用卷积神经网络或循环神经网络提取样本数据特征,再用其他分类方法进行分类,但都是依据单一类别特征进行的。
基于信息融合思想的复合模型。该思想通过综合考量多信息源数据、样本的多特征或多模型的输出,利用信息间的互补来提高分类器的判别精度。信息融合学习方式可分为数据层融合、特征层融合、决策层融合。目前三种融合方式在航发故障判别技术上均有应用。其中,决策层融合是将不同基本分类模型的决策结果进行融合,也有一些方法对同一模型多次训练的决策结果进行融合。Dempster-Shafer证据理论(D-S证据理论)、模糊积分和贝叶斯网络等是常用的决策层信息融合方法。
集成学习思想和多视图学习方法,也属于信息融合理论,但对于航发故障诊断这种复杂的非线性问题来说,用于集成的基本分类器也应是比较复杂的,如神经网络、支持向量机等。因此其集成的轮次不能过高,否则训练难度较大。此外,一些集成学习方法选择的基础弱分类器间同质性过高,导致其互补性不强,集成的提升效果不明显。
发明内容
现有的航发故障诊断机制大多基于单一模型的决策信息进行判别,未将优质模型的进行整合。信息融合与集成学习方法在此领域也有应用,但选择和构造的待融合基础模型存在一定的缺陷。
针对现有技术存在的问题,本发明提出一种基于深度学习和信息融合的航空发动机气路故障诊断方法。利用泛化能力较强的深度卷积神经网络(CNN)和深度前馈神经网络(BPNN) 分别对样本数据进行判别,获得模型对样本所属各个类别概率的评分,将其作为Dempster-Shafer证据理论的证据体。之后利用证据理论对两种证据体的分类结果进行融合,得到最终的分类结果。
为了达到上述目的,本发明采用的技术方案为:
一种基于深度学习和信息融合的航空发动机气路故障诊断方法,包括四个阶段:数据预处理;神经网络训练;构造证据体的基本概率分布;证据体的合成与决策。首先,在数据预处理阶段,对时序样本数据进行清洗、归一化处理、添加标签等操作;其次,在神经网络训练阶段,利用处理后的样本数据将卷积神经网络模型和前馈神经网络分别训练至最佳,继而得到两个模型对样本数据所属类别的判定结果;再次,在证据体的基本概率分布的构造阶段,两种模型的判别结果被视为两个证据体,这一阶段要将证据体中样本所属各类别评分转换为样本的基本概率分布,同时计算出样本的不确定度;最后,在证据体的合成与决策阶段,根据D-S证据理论的合成规则计算两个证据体合成后的概率分布,并根据决策规则,给出每条样本最后所属的类别。具体步骤如下:
(1)对数据进行预处理:本阶段将初始的大量的航发故障数据样本处理为能够直接输入神经网络模型运行、并便于神经网络对其进行特征感知和数据拟合的形式。
1)清除原始样本的序号、列名等无关数据,舍去序列样本数据首部和尾部的零值较多的样本值。
2)对样本按列进行归一化处理,消除指标量纲的影响,增强不同属性指标间的可比性。
采用的归一化公式如下:
Figure BDA0002390451140000031
式中:Xij表示样本中第i个序列数据的第j个属性;N表示样本中序列数据的个数;n表示每条样本的属性维度。max{·}表示取最大值;min{·}表示取最小值。
3)对样本逐条添加标签,按其所属的故障类别F1、F2...Fm分类,用1,2,3,...,m作为对应标签值,使样本的第一列数据为标签。
4)对不同故障类型的序列数据进行合并,合为一个文件后,进行随机打乱。
5)将上述样本按比例分割为训练集D、验证集V和测试集T,分别用于神经网络的训练、验证和测试。
(2)对神经网络模型进行充分训练,神经网络模型包括CNN模型和BPNN模型。两个模型计算单元中的权重矩阵W和偏置系数b要在模型初始化过程中给出,并利用样本数据进行反复的训练得以调整。模型层数、每层单元数等超参数也需要人工进行调整,以使得模型的判别精度达到最佳。CNN模型的训练按照以下步骤进行:
1)设定模型中的激活函数和损失函数。所述的常用的激活函数包括线性修正单元(RELU)、泄露型线性修正单元(Leaky-RELU)、挤压函数(Sigmoid)等,损失函数包括交叉熵损失和平方损失等。
2)设定初始的超参数。对于CNN,超参数为:卷积层数、全连接层数、每层卷积核尺寸、卷积核数量、卷积步长、各层的单元数和迭代次数。对航发故障诊断问题,神经网络的全连接层数采用二层或三层结构。对于训练的迭代次数,根据训练过程中损失函数的下降情况确定。对于航发数据这种时间序列数据,进行一维卷积,卷积核尺寸应小于数据的维度。
3)在给定的超参数下训练CNN模型。将训练集D和验证集V的样本数据输入模型,运行模型,等待模型训练完毕,记录本次训练后模型对测试集T的判别精度。每组超参数下,要进行3~5次上述操作,取判别精度的平均值作为模型对测试集的平均判别精度。
4)遵循控制变量法对模型中的超参数进行调整。将本组超参数下与上一组超参数下的模型平均判别精度作对比,之后按使平均判别精度增大的方向对模型某一超参数进行小幅度调整,保持其他参数不变。
5)重复步骤3)~4),保留使模型平均判别精度达到最高的一组超参数设置。
6)固定已调整完毕的模型的初始化参数及超参数设定,保存已训练好的模型。用训练好的模型对测试集T中样本进行判别,保存判别结果和从输入样本数据中剥离的测试集T的标签。
对BPNN模型进行训练,其需要设置和调整的超参数为:层数、各层单元数和迭代次数。训练的步骤和方法与上述的CNN模型训练过程完全一致。训练完毕后,用BPNN模型对测试集T中样本进行判别,并保存其判别结果。
(3)采用神经网络的判别结果构造证据体的基本概率分布。本发明构造证据体基本概率分布的方法与步骤如下:
1)将神经网络模型判别结果中对每个样本属于的m个基本类别的概率评分作为初始的概率分布。计算公式如下:
Figure BDA0002390451140000041
式中,R代表神经网络输出的判别结果。EM0代表样本的初始概率分布。角标i表示输入样本中的第i条序列。上标M代表某一证据体,M={A:CNN,B:BPNN}。l=1,2,3,...,m 指故障类别之一。
2)设定样本的不确定性度量指标及其计算方式。利用初始概率分布计算样本不确定性度量值。选取每个证据体初始概率分布到两证据体的平均初始概率分布值的距离平方和作为不确定性度量。计算公式如下:
Figure BDA0002390451140000042
Figure BDA0002390451140000043
式中:AVG_Ei表示两证据体对第i条样本的平均判别结果;EA0、EB0是EM0的具体取值,表示两个证据体的初始概率分布;
Figure BDA0002390451140000044
表示证据体M的判定结果到平均判定结果的距离平方和。
3)将不确定性度量转换为样本的不确定度U,依据U和步骤1)中的初始概率分布,按照公式(10)和(11),对其余m项命题的概率进行重新分配,使得每条样本的m+1个可能项的概率之和为1,能够得到证据体对每条样本的每个命题的信度函数(概率值)。
Figure BDA0002390451140000045
Figure BDA0002390451140000046
Figure BDA0002390451140000047
式中:
Figure BDA0002390451140000051
表示在证据体M中,样本i的基本概率分布中的第l项的概率值;
Figure BDA0002390451140000052
表示在证据体M中,样本i的不确定度,即样本归属于全集项的概率值;
Figure BDA0002390451140000053
表示证据体M对第i 条样本的判别结果的不确定度。
(4)确定证据体的融合方式和决策准则。D-S理论在融合证据体时,以识别框架Θ中的m个基础命题为出发视角,对不同证据体中各分布项的信度函数间的组合进行归类。依据合成后的信度函数,按决策规则可以获得样本最终的诊断决策。
1)将两个证据体的基本概率分布中各分布项两两组合。利用存在交集的证据组合,计算出合成后的各命题的信度函数。对于交集为空的证据组合,予以舍弃。如公式(12)和(13):
Figure BDA0002390451140000054
Figure BDA0002390451140000055
式中:k,h∈{F1、F2...Fm}∪{Θ},表示m+1个分布项中的两项;Ai(k)、Bi(h)表示针对第i条序列在A、B两证据体中分别取到的分布项;Ui表示样本第i条序列数据的合成判别结果的不确定度;Ei(l)表示按融合后的信度函数,第i条序列属于第l类的概率。
2)以信度函数为依据制定决策规则,取其概率最大的命题为该样本的决策结果。
3)利用步骤(2)中保存的测试集标签和经D-S理论融合模型的给出的决策结果,得到本发明的诊断精度。
本发明的有益效果:本发明针对航空发动机气路故障诊断问题,设计了一种基于卷积神经网络的特征感知能力、前馈神经网络的映射拟合能力和D-S证据理论的信息融合能力的判别方法。模型综合了两种深度神经网络从不同角度学习样本数据后给出的决策结果,二者互为补充,融合后可以提高判别精度。同时,可以克服非深度学习模型对不同来源的航发大数据判别能力不稳定、鲁棒性低的问题。
附图说明
图1是本发明基本框架。
图2是本发明的工作流程图。
具体实施方式
以下结合附图和技术方案,进一步说明本发明的具体实施方式。
一种基于深度学习和信息融合的航空发动机气路故障诊断方法,包括四个阶段:数据预处理阶段;神经网络训练阶段;构造证据体的基本概率分布阶段;证据体的合成与决策阶段。
第一步,在数据预处理阶段,将航发数据处理为能够直接输入神经网络模型运行的形式。同时,通过预处理,要使得样本更佳适用于神经网络对其进行特征感知和数据拟合。
1)清除原始样本的序列号、表头列名等无关数据,舍去序列样本数据首部和尾部第一列参数值为0的样本值。
2)初始样本数据的数值可能数量级过高、过低或不同属性间的单位量纲不统一。对样本按属性进行归一化处理,可消除指标之间量纲的影响,增强不同属性指标间的可比性,有助于模型更好的学习和判别数据。采用归一化公式如下:
Figure BDA0002390451140000061
式中:Xij表示样本中第i个序列数据的第j个属性;N表示样本中序列数据的个数;n表示每条样本的属性维度。max{·}表示取最大值;min{·}表示取最小值。
3)对样本逐条添加标签,将数据集改造为UCI标准时间序列数据集形式。按样本所属的故障类别F1、F2...Fm分类,用l=1,2,3,...,m作为对应标签值添加到样本的第一列。
4)对多种故障序列数据进行合并,合为一个文件后将样本序列进行随机打乱。利用CNN 和BPNN模型做序列分类问题时,不依赖样本间的顺序关系,打乱样本可以增加训练和测试样本的随机性,训练出更好的模型。
5)将样本按比例0.7、0.1和0.2的比例分割为训练集、验证集和测试集,保存为三个文件。用于神经网络的训练、验证和测试。
第二步,神经网络训练阶段,通过神经网络的超参数调整等工作,使每个作为证据体的神经网络的性能达到最佳。
1)设定模型中的激活函数和损失函数。常用的激活函数有:线性修正单元(RELU)、泄露型线性修正单元(Leaky-RELU)、挤压函数(Sigmoid)等。损失函数有交叉熵损失和平方损失等。激活函数一般默认采用RELU,分类问题一般采用交叉熵损失函数,分别如公式(2)和(3)所示:
Figure BDA0002390451140000062
Figure BDA0002390451140000063
式中:x表示某条序列数据;y表示样本真实标签值;
Figure BDA0002390451140000064
表示判别模型的输出值;W,b是各神经元中的参数;σ(·表示激活函数;L(·表示损失函数。
2)设定初始的超参数。CNN超参数为:卷积层数、全连接层数、每层卷积核尺寸、卷积核数量、卷积步长和各层的单元数。对于BPNN,要设置的超参数为:层数、各层单元数。此外学习率α,迭代次数等,在两种深度网络中都应进行设置。一般神经网络的全连接层数不应超过三层,对航发故障诊断问题,由于输入输出数据间的映射较为复杂,应采用二层或三层结构。对于训练的迭代次数,要根据训练过程中损失函数的下降情况确定,训练充分时,损失函数与训练精度会在一个值附近震荡,达到接近稳定的状态。对于航发数据这种时间序列数据,应进行一维卷积,即卷积核的第二维度尺寸为1。一维卷积核提取样本各属性间的内在联系。卷积核尺寸应小于数据的维度。
3)训练模型与调节超参数,超参数的调整遵循控制变量法。每次,将超参数设定后将样本数据输入模型,用公式(4)-(5)给出的误差反向传播的梯度下降法和链式法则对模型的各计算单元中参数W和b进行不断修正,训练完毕后记录该次模型对测试集的判别精度。每组超参数下进行3次实验,并把平均判别精度与上一组超参数下的模型平均判别精度作对比,按使模型平均判别精度增大的方向对某一超参数进行小幅度调整,保持其他参数不变。保留使模型平均判别精度达到最高的一组超参数设置。
Z(l)=W(l).x(l)+b(l) (4)
Figure BDA0002390451140000071
式中:Z表示净输入;α表示学习率;上标l表示第l层的数据;下标n表示第n轮训练后的数据。
4)固定已调整完毕的模型的初始化参数及超参数设定,可以保存已训练好的模型。用训练好的模型对测试集T中样本进行判别,保存模型对测试集的判别结果和从输入样本数据中剥离的测试集标签。
第三步,构造证据体的基本概率分布阶段。本发明构造证据体基本概率分布的方法与步骤如下:
1)将神经网络模型输出的对每个样本序列数据的所属m个基本类别命题的评分转化为初始的概率分布。计算公式如下:
Figure BDA0002390451140000072
式中,R代表神经网络输出的判别结果。EM0代表样本的初始概率分布。角标i表示输入样本中的第i条序列。上标M代表某一证据体,M={A:CNN,B:BPNN}。l=1,2,3,...,m指故障类别之一。
2)设定样本的不确定性度量指标以及计算方式。利用初始概率分布计算样本不确定性度量值。选取每个证据体初始概率分布到两证据体的平均初始概率分布值的距离平方和作为不确定性度量。计算公式如下:
Figure BDA0002390451140000081
Figure BDA0002390451140000082
式中:AVG_Ei表示两证据体对第i条样本的平均判别结果;EA0、EB0是EM0的具体取值,表示两个证据体的初始概率分布;
Figure BDA0002390451140000083
表示证据体M的判定结果到平均判定结果的距离平方和。
3)将不确定性度量转换为样本的不确定度U,依据U和1)中的初始概率分布,按公式(10) 和(11),对其余m项命题的概率进行重新分配,使得每条样本的m+1个可能项的概率之和为1。可以得到证据体对每条样本的每个命题的信度函数。
Figure BDA0002390451140000084
Figure BDA0002390451140000085
Figure BDA0002390451140000086
式中:
Figure BDA0002390451140000087
表示在证据体M中,样本i的基本概率分布中的第l项的概率值;
Figure BDA0002390451140000088
表示在证据体M中,样本i的不确定度,即样本归属于全集项的概率值;
Figure BDA0002390451140000089
表示证据体M对第i 条样本的判别结果的不确定度。
第四步,证据体的合成与决策阶段。确定证据体的融合方式和决策准则。D-S理论在融合证据体时,以识别框架Θ中的m个基础命题为出发视角,对不同证据体中各分布项的信度函数间的组合进行归类。为给出每个样本的确切判别结果,最终决策结果不能再存在表示不确定类别的全集项。依据合成后信度函数和决策规则可以得到样本最终的诊断决策。合理的证据体融合方式和决策规则可以对最终决策精度的提高起到积极作用。
1)将两个证据体中的各证据进行两两组合。利用存在交集的证据组合,计算出合成后的各证据的信度函数。对于交集为空的证据组合,对其予以舍弃。如公式(12)和(13):
Figure BDA00023904511400000810
Figure BDA0002390451140000091
式中:k,h∈{F1、F2...Fm}∪{Θ},表示m+1个分布项中的两项;Ai(k)、Bi(h)表示针对第i条序列在A、B两证据体中分别取到的分布项;Ui表示样本第i条序列数据的合成判别结果的不确定度;Ei(l)表示按融合后的信度函数,第i条序列属于第l类的概率。
2)决策规则:为了使每一条样本都有一个确切的决策结果,将不对各命题间的信度函数差值设阈值限制,直接取每条样本融合后的各证据的最大的项作为该样本的决策结果,如公式(14):
Yi=argMax(Ei(l)) (14)
式中:Yi表示模型对第i个样本序列的判别结果。
3)利用第二步中保存的测试集标签和经D-S理论融合后的决策结果,计算本发明的诊断精度。
结合本发明的方案,进行实验分析如下:
由于航发故障诊断领域目前还没有公开的标准数据集,且基于人工智能技术判别航发故障的方法大都未公开源代码,故难以与其他研究人员所提出的航发故障方法的判别精度进行对比。
但航发故障数据本质上是时间序列数据,对时间序列数据进行分类的方法都可以用于航发故障数据。可以借助在时间序列标准数据集上的实验及相关对比来证明本模型的有效性。
因此,对一种基于深度学习和信息融合的航空发动机气路故障诊断方法的实验分析过程可以分两部分进行,一是将方法应用在标准时间序列数据集上,与已有的其他时间序列分类模型判别结果进行横向对比;二是使用真实的航发故障数据,验证此方法应用在目标领域时的有效性。
(1)在UCI标准时间序列数据集上的实验
UCI标准系列数据集是由研究人员汇总的,在各实际应用领域采集的时间序列数据,如心电图数据、传感器数据及模拟数据等。这组数据集主要用来研究时间序列分类问题。
本方法在UCI标准数据集组中的Chlor.Conc、Cinc_ECG等数据集上进行了实验验证。上述数据集的基本情况如表1所示。
表1 实验使用的UCI标准数据集基本情况
Figure BDA0002390451140000101
本方法实验结果与时间序列数据分类领域的其他优秀判别模型的对比情况如表2所示。其中,CTN-T是一种应用了迁移学习思想的基于CNN提取出的单特征深度学习方法,BOSS 模型是一种基于距离的判别方法,这两种模型是研究人员于2019年提出的用于时间序列分类的新方法,具有较好的性能。
表2 UCI数据集上的实验结果比较
Figure BDA0002390451140000102
通过对实验结果的横向对比,可以看出深度学习模型总体上性能优于距离进行分类的方法,ResNet、CTN-s、CTN-T模型对上述数据集的总体正确率在90%左右,高于SP-HI、BOTSW-BD、BOSS方法。在上述数据集中,基础的深度神经网络模型也能达到超出其他复杂方法的判别精度。经过模型融合后可以得到更好的结果,即使在基于单特征的判别模型正确率达到99%的超高精度情况下,融合模型仍能够提升判别精度。
(2)在航发故障数据集上的实验
1)航发故障数据集情况:
本实验使用的航发故障数据集来源于中国航发沈阳发动机研究所(606所),数据集中有 A、B两种发动机机型的故障数据,每种机型数据有23维属性、5种故障模式。经过实验对比,共取每种机型60000条各种故障类别样本时,训练出的模型效果最好。
2)数据处理:
对上述样本进行如下处理:
数据清洗:去除各类故障数据中首尾的0元较多的样本序列;
数据归一化:将样本各维属性按公式进行归一化,增加属性间的可比性;
数据集分割:取70%的样本数据作为训练和验证集,30%留作测试集;
加入噪声模拟真实环境:本数据集中数据的采集环境为温度和气压相对稳定的实验室,为还原高空真实的恶劣工作环境的对数据采集的影响,随机选取了四分之一的样本数据加入了高斯白噪声。
实验步骤按具体实施方式中描述的步骤进行,得到的两个单特征模型的最优状态下超参数及其对故障数据集的判别精度如表3和表4所示。表4中也给出了利用D-S证据理论融合后的模型的正确率。
表3 航发故障数据集实验两种神经网络最佳状态下超参数
Figure BDA0002390451140000111
表4 航发故障数据集实验结果
Figure BDA0002390451140000112
在两种机型的故障数据集上,与单特征的最大值相比,融合后模型的精度均提高了0.25%。实验结果表明利用D-S证据理论将两种单特征模型进行决策融合可以有效提升判别精度,此模型对于航发故障分类问题具有很好的效果。
以上所述实施例仅表达本发明的实施方式,但并不能因此而理解为对本发明专利的范围的限制,应当指出,对于本领域的技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些均属于本发明的保护范围。

Claims (4)

1.一种基于深度学习和信息融合的航空发动机气路故障诊断方法,其特征在于,所述航空发动机气路故障诊断方法,首先,对数据进行预处理;其次,在神经网络训练阶段,利用处理后的样本数据将卷积神经网络模型和前馈神经网络分别训练至最佳,继而得到两个模型对样本数据所属类别的判定结果;再次,在证据体的基本概率分布的构造阶段,两种模型的判别结果被视为两个证据体,该阶段将证据体中样本所属各类别评分转换为样本的基本概率分布,同时计算出样本的不确定度;最后,在证据体的合成与决策阶段,根据D-S证据理论的合成规则计算两个证据体合成后的概率分布,并根据决策规则,给出每条样本最后所属的类别;步骤如下:
(1)对数据进行预处理:将初始的大量的航发故障数据样本处理为能够直接输入神经网络模型运行、并便于神经网络对其进行特征感知和数据拟合的形式,将处理后的样本按比例分割为训练集D、验证集V和测试集T,分别用于神经网络的训练、验证和测试;
(2)对神经网络模型进行充分训练,神经网络模型包括CNN模型和BPNN模型;两个模型计算单元中的权重矩阵W和偏置系数b要在模型初始化过程中给出,并利用样本数据进行反复的训练得以调整;模型层数、每层单元数也需要人工进行调整,以使得模型的判别精度达到最佳;CNN模型的训练按照以下步骤进行:
1)设定模型中的激活函数和损失函数;
2)设定初始的超参数;对于CNN,超参数为:卷积层数、全连接层数、每层卷积核尺寸、卷积核数量、卷积步长、各层的单元数和迭代次数;
3)在给定的超参数下训练CNN模型;将训练集D和验证集V的样本数据输入模型,运行模型,等待模型训练完毕,记录本次训练后模型对测试集T的判别精度;每组超参数下,进行多次上述操作,取判别精度的平均值作为模型对测试集的平均判别精度;
4)遵循控制变量法对模型中的超参数进行调整;将本组超参数下与上一组超参数下的模型平均判别精度作对比,之后按使平均判别精度增大的方向对模型某一超参数进行小幅度调整,保持其他参数不变;
5)重复步骤3)~4),保留使模型平均判别精度达到最高的一组超参数设置;
6)固定已调整完毕的模型的初始化参数及超参数设定,保存已训练好的模型;采用训练好的模型对测试集T中样本进行判别,保存判别结果和从输入样本数据中剥离的测试集T的标签;
对BPNN模型进行训练,其需要设置和调整的超参数为:层数、各层单元数和迭代次数;训练方法与CNN模型训练过程完全一致;训练完毕后,采用BPNN模型对测试集T中样本进行判别,并保存其判别结果;
(3)采用神经网络的判别结果构造证据体的基本概率分布;
1)将神经网络模型判别结果中对每个样本属于的m个基本类别的概率评分作为初始的概率分布;计算公式如下:
Figure FDA0003308971180000021
式中,R代表神经网络输出的判别结果;EM0代表样本的初始概率分布;角标i表示输入样本中的第i条序列;上标M代表某一证据体,M={A:CNN,B:BPNN};l=1,2,3,...,m指故障类别之一;
2)设定样本的不确定性度量指标及其计算方式;利用初始概率分布计算样本不确定性度量值;选取每个证据体初始概率分布到两证据体的平均初始概率分布值的距离平方和作为不确定性度量;计算公式如下:
Figure FDA0003308971180000022
Figure FDA0003308971180000023
式中:AVG_Ei表示两证据体对第i条样本的平均判别结果;EA0、EB0是EM0的具体取值,表示两个证据体的初始概率分布;
Figure FDA0003308971180000024
表示证据体M的判定结果到平均判定结果的距离平方和;
3)将不确定性度量转换为样本的不确定度U,依据U和步骤1)中的初始概率分布,按照公式(10)和(11),对其余m项命题的概率进行重新分配,使得每条样本的m+1个可能项的概率之和为1,能够得到证据体对每条样本的每个命题的信度函数;
Figure FDA0003308971180000025
Figure FDA0003308971180000026
Figure FDA0003308971180000027
式中:
Figure FDA0003308971180000028
表示在证据体M中,样本i的基本概率分布中的第l项的概率值;
Figure FDA0003308971180000029
表示在证据体M中,样本i的不确定度,即样本归属于全集项的概率值;
Figure FDA00033089711800000210
表示证据体M对第i条样本的判别结果的不确定度;
(4)确定证据体的融合方式和决策准则;D-S理论在融合证据体时,以识别框架Θ中的m个基础命题为出发视角,对不同证据体中各分布项的信度函数间的组合进行归类;依据合成后的信度函数,按决策规则可以获得样本最终的诊断决策;
1)将两个证据体的基本概率分布中各分布项两两组合;利用存在交集的证据组合,计算出合成后的各命题的信度函数;对于交集为空的证据组合,予以舍弃;如公式(12)和(13):
Figure FDA0003308971180000031
Figure FDA0003308971180000032
式中:k,h∈{F1、F2...Fm}∪{Θ},表示m+1个分布项中的两项;Ai(k)、Bi(h)表示针对第i条序列在A、B两证据体中分别取到的分布项;Ui表示样本第i条序列数据的合成判别结果的不确定度;Ei(l)表示按融合后的信度函数,第i条序列属于第l类的概率;
2)以信度函数为依据制定决策规则,取其概率最大的命题为该样本的决策结果;
3)利用步骤(2)中保存的测试集标签和经D-S理论融合模型的给出的决策结果,得到最终诊断精度。
2.根据权利要求1所述的一种基于深度学习和信息融合的航空发动机气路故障诊断方法,其特征在于,步骤(1)中所述的对数据进行预处理具体为:
1)清除原始样本的无关数据,舍去序列样本数据首部和尾部的零值较多的样本值;
2)对样本按列进行归一化处理;
3)对不同故障类型的序列数据进行合并,合为一个文件后,进行随机打乱;
4)将上述样本按比例分割为训练集D、验证集V和测试集T,分别用于神经网络的训练、验证和测试。
3.根据权利要求1所述的一种基于深度学习和信息融合的航空发动机气路故障诊断方法,其特征在于,步骤(2)的步骤1)中所述的激活函数包括线性修正单元、泄露型线性修正单元、挤压函数,所述的损失函数包括交叉熵损失、平方损失。
4.根据权利要求1所述的一种基于深度学习和信息融合的航空发动机气路故障诊断方法,其特征在于,步骤(2)的步骤2)中对航发故障诊断问题,所述CNN模型超参数的全连接层数采用二层或三层结构;所述迭代次数根据训练过程中损失函数的下降情况确定;对于航发数据这种时间序列数据,进行一维卷积,卷积核尺寸小于数据的维度。
CN202010112320.1A 2020-02-24 2020-02-24 一种航空发动机气路故障诊断方法 Active CN111368885B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010112320.1A CN111368885B (zh) 2020-02-24 2020-02-24 一种航空发动机气路故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010112320.1A CN111368885B (zh) 2020-02-24 2020-02-24 一种航空发动机气路故障诊断方法

Publications (2)

Publication Number Publication Date
CN111368885A CN111368885A (zh) 2020-07-03
CN111368885B true CN111368885B (zh) 2021-12-03

Family

ID=71211530

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010112320.1A Active CN111368885B (zh) 2020-02-24 2020-02-24 一种航空发动机气路故障诊断方法

Country Status (1)

Country Link
CN (1) CN111368885B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111814870B (zh) * 2020-07-06 2021-05-11 北京航空航天大学 一种基于卷积神经网络的cps模糊测试方法
CN111858328B (zh) * 2020-07-15 2021-11-12 南通大学 一种基于有序神经网络的软件缺陷模块严重程度预测方法
CN112036479A (zh) * 2020-08-29 2020-12-04 大连海事大学 一种船舶空调系统故障识别方法、装置及存储介质
CN112052628B (zh) * 2020-09-06 2022-09-16 上海电力大学 一种基于知识+数据的燃气轮机气路故障预测诊断方法
CN114663333A (zh) * 2020-12-08 2022-06-24 重庆云石高科技有限公司 一种机车小尺度零部件丢失的智能检测系统
CN112699597A (zh) * 2020-12-09 2021-04-23 西安交通大学 一种核电启动水泵滚动轴承故障检测方法及系统
CN112613584A (zh) * 2021-01-07 2021-04-06 国网上海市电力公司 一种故障诊断方法、装置、设备及存储介质
CN112749764B (zh) * 2021-01-28 2022-09-13 厦门大学 一种基于qar数据的航空发动机运行状态分类方法
CN113865859B (zh) * 2021-08-25 2024-05-14 西北工业大学 多尺度多源异构信息融合的齿轮箱状态故障诊断方法
CN113780292B (zh) * 2021-08-31 2022-05-06 北京交通大学 一种基于证据推理的语义分割网络模型不确定性量化方法
CN114330517A (zh) * 2021-12-15 2022-04-12 南京航空航天大学 一种基于神经网络的航空发动机传感器故障自诊断方法
CN116807435A (zh) * 2023-08-28 2023-09-29 嘉兴学院G60科创走廊产业与创新研究院 一种便携式心音传感装置及心率监测系统
CN117807521A (zh) * 2023-12-12 2024-04-02 四川航空股份有限公司 基于航空发动机疲劳可靠性分析的故障预测方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105911975A (zh) * 2016-04-14 2016-08-31 南京航空航天大学 一种航空发动机气路故障融合诊断方法
CN107655693A (zh) * 2017-11-02 2018-02-02 哈尔滨理工大学 一种船舶发动机故障诊断系统及方法
CN109669087A (zh) * 2019-01-31 2019-04-23 国网河南省电力公司 一种基于多源信息融合的电力变压器故障诊断方法
JP2019122008A (ja) * 2018-01-11 2019-07-22 株式会社Nttドコモ 電波伝搬推定装置、電波伝搬推定方法、及び、電波伝搬推定プログラム
CN110109005A (zh) * 2019-05-24 2019-08-09 电子科技大学 一种基于序贯测试的模拟电路故障测试方法
CN110321603A (zh) * 2019-06-18 2019-10-11 大连理工大学 一种用于航空发动机气路故障诊断的深度计算模型

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105911975A (zh) * 2016-04-14 2016-08-31 南京航空航天大学 一种航空发动机气路故障融合诊断方法
CN107655693A (zh) * 2017-11-02 2018-02-02 哈尔滨理工大学 一种船舶发动机故障诊断系统及方法
JP2019122008A (ja) * 2018-01-11 2019-07-22 株式会社Nttドコモ 電波伝搬推定装置、電波伝搬推定方法、及び、電波伝搬推定プログラム
CN109669087A (zh) * 2019-01-31 2019-04-23 国网河南省电力公司 一种基于多源信息融合的电力变压器故障诊断方法
CN110109005A (zh) * 2019-05-24 2019-08-09 电子科技大学 一种基于序贯测试的模拟电路故障测试方法
CN110321603A (zh) * 2019-06-18 2019-10-11 大连理工大学 一种用于航空发动机气路故障诊断的深度计算模型

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Aero engine gas path fault prediction based on multi-sensor information fusion;Yang Xiaohong 等;《2016 IEEE Chinese Guidance, Navigation and Control Conference (CGNCC)》;20170123;全文 *
Bearing Fault Diagnosis with a Feature Fusion Method Based on an Ensemble Convolutional Neural Network and Deep Neural Network;Hongmei Li 等;《Sensors》;20190430;全文 *
基于多源异构数据融合的深度学习故障诊断;何一帆;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》;20200115;第2020年卷(第01期);全文 *
基于深度学习的轨道转辙机故障检测系统;孙迪钢;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》;20190115;第2018年卷(第12期);全文 *

Also Published As

Publication number Publication date
CN111368885A (zh) 2020-07-03

Similar Documents

Publication Publication Date Title
CN111368885B (zh) 一种航空发动机气路故障诊断方法
AU2020104006A4 (en) Radar target recognition method based on feature pyramid lightweight convolutional neural network
US7362892B2 (en) Self-optimizing classifier
CN109190665A (zh) 一种基于半监督生成对抗网络的通用图像分类方法和装置
CN112766303B (zh) 一种基于cnn的航空发动机故障诊断方法
CN106897738A (zh) 一种基于半监督学习的行人检测方法
Özbılge et al. Tomato disease recognition using a compact convolutional neural network
CN111680875B (zh) 基于概率基线模型的无人机状态风险模糊综合评价方法
CN110210380B (zh) 基于表情识别与心理学测试生成性格的分析方法
CN113919220B (zh) 一种振动图像驱动的滚动轴承智能故障诊断方法
CN112052933B (zh) 基于粒子群优化的深度学习模型的安全性测试方法和修复方法
Arzar et al. Butterfly species identification using convolutional neural network (CNN)
US7233692B2 (en) Method and computer program product for identifying output classes with multi-modal dispersion in feature space and incorporating multi-modal structure into a pattern recognition system
CN116881841A (zh) 一种基于F1-score多级决策分析的混合模型故障诊断方法
CN115051864B (zh) 基于pca-mf-wnn的网络安全态势要素提取方法及系统
CN115994302A (zh) 基于柯西核贝叶斯自适应共振网络的航空发动机故障诊断方法
CN116977834B (zh) 一种开放条件下分布内外图像识别方法
Wang et al. Interpret neural networks by extracting critical subnetworks
Liu et al. MRD-NETS: multi-scale residual networks with dilated convolutions for classification and clustering analysis of spacecraft electrical signal
CN108898157B (zh) 基于卷积神经网络的数值型数据的雷达图表示的分类方法
CN115526227A (zh) 一种直升机飞行数据瞬时异常检测方法
CN112990279B (zh) 基于自动编码器的雷达高分辨距离像库外目标拒判方法
Khotanzad et al. Hand written digit recognition using BKS combination of neural network classifiers
Kohsasih et al. Comparison of CNN Architecture for White Blood Cells Image Classification
Watanabe et al. Understanding community structure in layered neural networks

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