航空发动机性能图像纹理片段划分方法.pdf

收稿日期：2021-12-17作者简介：蔡舒妤（1985），女，硕士，副教授。通信作者：师利中（1984），男，硕士，副教授。引用格式：蔡舒妤，史涛，师利中.航空发动机性能图像纹理片段划分方法J.航空发动机，2023，49（4）：86-94.CAI Shuyu，SHI Tao，SHI Lizhong.Texture segmentation method for aeroengine performance imageJ.Aeroengine，2023，49（4）：86-94.第 49 卷 第 4 期2023 年 8 月Vol.49 No.4Aug.2023航空发动机Aeroengine航空发动机性能图像纹理片段划分方法蔡舒妤，史涛，师利中（中国民航大学 航空工程学院，天津 300300）摘要：为了解决深度学习在航空发动机性能异常检测中出现的数据不平衡的问题，在高频次纹理片段对应性能稳定状态、低频次纹理片段对应异常波动状态的合理假设下，提出了一种航空发动机性能图像纹理片段划分方法。通过研究性能数据空间与RGB彩色空间的映射，提出高维性能数据图像化。引入分形盒维理论，定义纹理特征维数表征性能图像纹理特征，提出基于高频次生长树的图像纹理片段划分方法，实现性能图像按纹理出现频次的精细划分，建立正常样本和异常样本数量接近的数据集。选取多个航段发动机性能数据，与时频图和采样法生成的数据集进行对比验证。验证结果表明：该方法可以生成数量相当且真实客观的正常样本和异常样本，在ResNet50模型上的准确率达91.72%，减小了数据不平衡对准确率的影响。关键词：性能图像；纹理特征维数；高频次生长树；图像纹理片段划分；深度学习；航空发动机中图分类号：V239文献标识码：Adoi：10.13477/ki.aeroengine.2023.04.011Texture Segmentation Method For Aeroengine Performance ImageCAI Shu-yu，SHI Tao，SHI Li-zhong（College of Aeronautical Engineering，Civil Aviation University of China，Tianjin 300300，China）Abstract：In order to solve the problem of data imbalance in aeroengine performance anomaly detection by deep learning，a texturesegmentation method for aeroengine performance images was proposed under the reasonable assumption that high-frequency texture segments correspond to stable performance states and low-frequency texture segments correspond to abnormal fluctuation states.By studyingthe mapping between performance data space and RGB color space，the high dimensional performance data pictorialization was proposed.The fractal box dimension theory was introduced to define the texture feature dimension to characterize the performance image texture features，and the image texture segmentation method based on a high-frequency growth tree was proposed to realize the fine segmentation ofthe performance image by the frequency of texture occurrence，achieving the balance of normal and abnormal samples in the data set.Theengine performance data of several flight segments were selected and compared with the data sets generated by the time-frequency plotsand sampling methods for verification.The verification results show that the method can generate a considerable number of realistic and objective normal samples and abnormal samples，with an accuracy of 91.72%on the ResNet50 model，reducing the impact of data imbalance on accuracy.Key words：performance image；texture feature dimension；high-frequency growth tree；image texture segmentation；deep learning；aeroengine0引言航空发动机作为飞机的核心部件，其健康状态和性能极大地影响着飞机的安全，发动机性能数据作为反映其运行状态的重要依据，为异常诊断与预警1-3、飞行品质监控与改善4-6、排放与环保评测7-9等方面的研究提供了重要的数据支持。近年来，深度学习智能算法为航空发动机的故障诊断提供了新思路，基于性能数据的智能诊断方法成为研究热点。车畅畅等10提出了基于深度学习的发动机故障融合诊断模型，利用深度学习模型得出故障分类置信度，后对多次故障分类结果进行决策融合，得出更准确的诊断结果；钟诗胜等11分析民航发动机飞行历史数据，结合深度置信网络建立基于不均衡样蔡舒妤等：航空发动机性能图像纹理片段划分方法第 4 期图1发动机性能数据图像化映射本驱动的民航发动机故障诊断模型；Zhang等12针对基于卷积神经网络忽略数据的时间序列信息问题，提出了一种基于递归神经网络的旋转机械故障类型识别方法，并在2个公共数据集上取得了很好的效果；Park等13提出了一种基于深度学习的液体火箭发动机起动瞬态故障检测与诊断方法，建立发动机起动过程中潜在故障类型的数据集，用于训练深层神经网络，并利用实际发射地面试验数据验证了该方法的有效性；Guo等14提出基于混合特征模型和深度学习的无人机传感器故障诊断方法，利用短时傅立叶变换将故障信号变换为相应的时频图，然后利用卷积神经网络提取图像特征，实现了无人机传感器的故障诊断；Chen等15针对电气传动系统，提出了一种基于数据驱动和深度学习相结合的方法来处理早期故障；Manjurul I等16提出了一种自适应深度卷积神经网络，利用试验台采集的数据验证了所提出的故障诊断方法的有效性。然而深度学习在故障诊断领域仍然存在不适应性。基于深度学习的故障诊断方法需要大量高质量的数据，才能取得良好的诊断性能。然而实际的发动性能数据故障样本远远小于正常样本，会影响深度学习模型的检测结果，使诊断模型结果向大样本倾斜17。目前解决不平衡故障诊断问题时通常有2种手段：第1种是从数据增强方面入手，如利用生成对抗网络（Generative Adversarial Network，GAN）18-19生成和原始数据相似的数据样本；第2种是对数据进行预处理，如采用过采样和欠采样技术来平衡故障数据和正常数据之间的数量差异。然而，通过GAN进行训练对样本数量要求严苛，训练样本太小会导致模型无法学习到真实样本的分布，训练样本太多会脱离发动机实际工况；对于过采样预处理，会复制少数样本以扩充样本量，但容易出现过度拟合问题，而欠采样预处理减少了大多数类的样本数，难免造成信息丢失。针对以上问题，本文提出航空发动机性能图像纹理片段划分方法。1发动机性能图像的纹理特征表示1.1发动机性能数据图像化发动机性能数据图像化通过建立各发动机性能数据值域与RGB彩色空间的映射关系，将高数据维度的性能参数数据转换为2维性能图像，保持性能数据的细节，为性能参数之间和性能数据之间的关联分析提供高效、直观的支持。对于任意发动机性能参数，其性能数据值域空间imin,imax与 RGB彩色空间0,255的映射关系，如图 1所示。假设vij表示第i个性能参数在第j个时刻采集的数据值，则vijimin,imax其中imin为第i个性能参数取值的最小值；imax为第i个性能参数取值的最大值，则发动机性能数据图像化公式为RGB=(vij-imin)/(imax-imin)/256)（1）式中：RGB为计算所得颜色分量值。由此转化的性能图像横轴为时间，纵轴为各性能参数取值。1.2基于分形盒维的性能图像纹理特征表示分形维数理论是对非光滑、非规则、非完整等的极其复杂形体定量分析的重要参数，是形体的复杂度、粗糙度的度量20。可将分形维数理论引入发动机性能图像的纹理特征表示21-22。假设大小为M N的发动机性能图像中，任意像素点(x,y)定义其灰度值函数为I(x,y)，均满足I(x,y)0,255，且1 x M，1 y N。构建图像的灰度空间，如图 2 所示。性能图像灰度空间中，X轴、Y轴分别表示图像像素点的坐标值，Z轴表示图像像素点的灰度值。取边长为的分形盒，像素点(x,y)的分形盒维数为N(x,y)=I(x,y)/（2）定 义 像 素 点(x,y)处 的 纹 理 特 征 维 数 为N(x,y)，即以像素点(x,y)为中心，边长为的邻域图2灰度空间87航空发动机第 49 卷内像素分形盒维数最大值和最小值的差值N(x,y)=NMAX(x+i,y+i)i (-/2,/2),j (-/2,/2)-NMIN(x+i,y+i)i (-/2,/2),j (-/2,/2)（3）式中：NMAX(x+i,y+i)i (-/2,/2),j (-/2,/2)表示以像素点(x,y)为中心()区域内像素分形盒维数的最大值，记为NMAX (x,y)；NMIN(x+i,y+i)i (-/2,/2),j (-/2,/2表示以像素点(x,y)为中心()区域内像素分形盒维数的最小值，记为NMIN (x,y)。纹理特征维数如图 3所示。纹理特征维数表示了像素点-邻域内的纹理特征，通过图像化方法将独立且离散的性能数据相融合。这种融合既包含了不同参数间的关联，同时包含了不同时刻参数的关联。2基于高频次生长树的图像纹理片段划分发动机性能图像中，纹理片段出现的不同频次表征了性能状态的不同。高频次纹理片段体现性能的稳定，低频次纹理片段体现性能的波动。可通过对性能状态频次的提取，实现高低频次纹理图像纹理片段的划分。2.1高频次生长树理论依据发动机性能数据特点，性能图像中的每一列对应于飞机某时刻各参数的性能状态，每一行对应于某个参数各飞行时刻的性能状态。可通过性能图像每列像素的纹理特征维数表示图像性能状态。Sk=N(0,k),N(1,k),N(M-1,k),N(M,k)（4）式中：k为图像性能状态的列数，1 k N。高频次生长树理论如图4所示。将发动机性能图像中每个性能状态定义为1颗种子，其初始状态为S()0k()1 k N，计算该性能状态在图像中出现的频次Pk。当频次Pk大于等于一定阈值T时，种子开始第1轮生长S()1k=S()0k-1+S()0k+S()0k+1；第n轮生长可表示为S()nk=S()n-1k-1+S()n-1k+S()n-1k+1。当频次Pk小于一定阈值T时，种子则停止生长。该过程持续进行，直至所有种子停止生长。2.2基于高频次生长树的图像纹理片段划分算法本文提出一种基于高频次生长树的图像纹理片段划分算法建立数据集，算法流程如图5所示。算法的具体步骤：（1）性能数据图像化。输入发动机性能参数，确图3纹理特征维数图4高频次生长树Sk=N(0,k),N(1,k),N(M-1,k),N(M,k)RGB=（vij-imin）/（imax-imin）/256）I（x,y）0,255N（x,y）=I（x,y）/N(x,y)=NMAX(x+i,y+i)i (-/2,/2),j (-/2,/2)-NMIN(x+i,y+i)i (-/2,/2),j (-/2,/2)Pk,Pk TPk,PkSk=Sk-1+Sk+Sk+1Pk图5算法流程88蔡舒妤等：航空发动机性能图像纹理片段划分方法第 4 期航段ABCDEFGHIJKL数据量11085677166819521104135936541154217868734559476840运行状态正常正常正常正常正常正常正常正常正常正常振动异常振动异常表1各航段性能数据定其性能数据值域空间imin,imax。（2）依据式（1）建立与彩色空间的映射关系，实现发动机性能数据图像化。（3）计算分形盒维数。取边长为的分形盒，根据式（2）计算性能图像中每个像素点的盒维数N(x,y)。（4）计算纹理特征维数。计算像素点邻域内分形盒维数的最大值和最小值，根据式（3）计算纹理特征维数N。（5）计算性能状态的频次Pk。记图像中每列特征维数为一种性能状态，计算各性能状态在性能图像中出现的频次Pk。（6）高频次生长。当性能状态的频次Pk小于阈值T则淘汰，大于阈值T时输出Pk，然后依据式（4）开始生长，计算新性能状态的频次Pk再与阈值T比较，重复此过程直至所有性能状态生长完成；（7）依据性能状态频次，对图像进行划分。3实例验证与分析3.1性能数据图像化方法验证为验证基于高频次生长树纹理分析的发动机性能图像纹理片段划分方法，选取12个航段的发动机性能参数数据，包括左右发低压转子振动值N1V（lowpressure compressor rotor vibrate）和高压转子振动值N2V（high pressure compressor rotor vibrate），本文采用的数据均来自于某航空公司译码后的QAR数据，其数据信息见表1。3.1.1本文图像化方法验证依据性能数据与颜色空间映射关系，将性能数据转化成2维性能图像，各航段性能数据对应性能图像如图6所示。从图中可见，对于正常航段AJ所对应的性能图像，其存在明显的主色调区域，局部区域呈现相对统一的颜色，如图中红框所示，说明该性能状态平稳。对于异常航段K和L，其性能图像纹理细密，说明性能状态波动较大，平稳性差，使得图像颜色变化较大，难以维持局部稳定色调。3.1.2时频图像化方法验证常用的图像化方法有短时傅里叶变换法，利用短时傅里叶变换法对发动机性能数据进行图像化，生成的时频图像如图7所示。（c）航段C（d）航段D（e）航段E（f）航段F（g）航段G（h）航段H（i）航段I（j）航段J（k）航段K（l）航段L图6各航段性能图像（a）航段A（b）航段B89航空发动机第 49 卷航段ABCDEF本文方法2.4271.0540.9331.9692.2800.807时频法0.9320.6180.6080.7690.8290.564航段GHIJKL本文方法0.6780.7021.4281.2210.8561.124时频法0.5290.5360.6750.6410.5660.635表2图像化时间结果对比航段ABCDEF本文方法7.52997.77277.71597.48597.62127.7492时频法6.94376.87856.91086.77176.94736.8045航段GHIJKL本文方法7.52057.48837.55117.37417.90507.9016时频法6.50576.64646.76806.78946.94856.8221表3信息熵计算结果对比航段ABCDEFG低频次图像纹理片段数量1931451421671249971高频次图像纹理片段数量1931441411661239870航段HIJKL合计低频次图像纹理片段数量609195111189高频次图像纹理片段数量609195001181表4性能图像划分2种方法图像化所需要的时间见表2。从表中可见，本文提出的图像化方法的运算效率优于时频法的。图像的细节纹理可以用图像信息熵的大小来衡量，反映了图像中平均信息量的多少，信息熵越大，说明图像中所包含的细节纹理越多。H(i,j)=i=0255pijlogpij（5）式中：i表示某像素点的灰度值；j表示某像素点领域灰度值；pij表示某像素点的2维特征。信息熵计算结果对比见表3。从表中可见，本文方法生成的图像信息量远远大于时频法所得的图像，在融合不同时刻不同性能参数的情况下，保留性能数据更多的细节，减小信息损失。3.2数据集建立验证3.2.1本文方法建立的数据集在性能图像的基础上，依据基于高频次生长树的图像纹理片段划分算法，对性能图像纹理片段进行划分，结果见表4，从表中可见，本文提出的算法可以对完整的性能图像进行精细划分，将其划分为高频次图像纹理片段和低频次图像片段。其中低频次图像纹理片段1189张，高频次图像纹理片段1181张，消除了数据不平衡。对于正常航段，划分后的低频次图像纹理片段数量与高频次图像纹理片段数量接近，增强了低频次图像纹理片段的数量，并且真实客观，不损失信息细节；对于异常航段，均为低频次纹理片段，印证了前文的假设，即高频次纹理片段正常，低频次纹理片段异常。高、低频次图像纹理片段部分样本分别如图8、9所示。从图中可见，高频次图像纹理片段与低频次图像纹理片段差异明显，高频次图像纹理片段多数较宽，表明该性能状态平稳，低频次图像纹理片段窄小，表明该性能状态波动大。（a）航段A（b）航段B（c）航段C（d）航段D（e）航段E（f）航段F（g）航段G（h）航段H（i）航段I（j）航段J（k）航段K（l）航段L图7各航段时频图像90蔡舒妤等：航空发动机性能图像纹理片段划分方法第 4 期参数名称batch_sizeepoch学习率对应设置165000.001参数名称优化方法激活函数对应设置AdamRelu表8网络模型参数设置根据式（5）计算高频次和低频次图像纹理片段的信息熵，结果见表5。从表中可见，高频次图像纹理片段的熵值均小于低频次图像纹理片段熵值，这是因为高频图像纹理片段在纹理变化上较为稳定，能够维持稳定色调，低频图像纹理片段颜色变化波动较大，难以维持稳定色调，进而导致信息熵增加。选取正常航段 AJ与异常航段 LL中的高频次和低频次图像纹理片段，按照不同比例构建数据集，见表6。3.2.2时频方法建立的数据集时频图像中横轴每一点对应一列性能数据，以100列数据为一组对图像进行划分，建立数据集d，结合过采样方法建立数据集e，结合欠采样方法建立数据集f，见表7。3.3应用验证与分析基于VGG16、ResNet18和ResNet50深度学习网络模型对本文方法建立的数据集 a、b和 c，和传统方法建立的数据集d、e和f分别进行对比验证，网络模型参数设置见表8。3.3.1本文方法建立的数据集上的训练准确率VGG16、ResNet18 和ResNet50模型在数据集a、b 和 c 上训练过程的准确率变化如图1012所示。从图中可见，VGG16模型在数据集a上，随着训练次数增加，准确率略微上升，在数据集b和数据集c上，准确率先突增后保持不变。造成这种结果是因为数据集a中正常样本和异常样本数量几乎一致，浅层卷积神经网络无法提取更深层次的图像特征，最终导致准确率很低。数据集b中正常样本和异常样本数量的比例为1 2，数据集c中为2 1，VGG16模型将所有输入样本识别为数量占多数的样本类型，反而导致准确率高于数据集a，可以看出浅层卷积神经网络效果不好。ResNet18和ResNet50这2种模型的准确率变化趋势相似，首先随着训练次数增加，准确率先不断上升，上升到一定程度后开始稳定，最终达到收敛。其中数据集a的准确率最高，数据集b和数据集c的准图8高频次图像纹理片段部分样本图9低频次图像纹理片段部分样本高频次图像纹理片段低频次图像纹理片段序号（I）（II）（III）（IV）平均值（I）（II）（III）（IV）平均值信息熵7.07696.63527.21136.89206.95397.50037.55057.45207.64137.5360表5信息熵计算结果数据集abc比例1 1.011 22 1高频次图像纹理片段数11815951181低频次图像纹理片段数11891189591表6本文方法构建的数据集（）（）（）（）（）（）（）（）数据集def正常纹理片段数770770129异常纹理片段数129770129表7时频法构建的数据集图10VGG16在数据集a、b和c的准确率曲线图11ResNet18在数据集a、b和c的准确率曲线图12ResNet50在数据集a、b和c的准确率曲线91航空发动机第 49 卷确率相近且都低于数据集a，这是因为数据集a中的高频和低频图像纹理片段数量接近，深度学习模型能够准确的学习高频次和低频次图像纹理片段的特征，消除了数据不平衡对识别精度的影响。3.3.2时频法建立的数据集上的训练准确率VGG16、ResNet18 和ResNet50模型在数据集e、b 和 c 上训练过程的准确率变化如图1315所示。从图中可见，VGG16模型在数据集d和f上，准确率先迅速上升后保持不变，在数据集e上准确率逐渐上升，最终达到收敛。数据集d中正常样本和异常样本数量的比例接近6 1，正常样本数量远远大于异常样本数量，VGG16模型将所有输入样本识别为正常样本，反而达到较高的准确率。ResNet18和ResNet50这2种模型的准确率变化趋势相似，首先随着训练次数增加，准确率快速上升，到一定程度后开始稳定，最终达到收敛。经过采样方法建立的数量平衡数据集e和f相较于数量不平衡的数据集 d 在准确率上提升不明显，且 ResNet18 和ResNet50模型在数据集d、e和f上的准确率均低于数据集a、b和c。选取各数据集 80%作为训练集，20%作为测试集，训练准确率和测试准确率作为评价指标，分别为训练集准确率=训练集中检测正确样本数训练集样本总数（6）测试集准确率=测试集中检测正确样本数测试集样本总数（7）以数据集a为例，在VGG16模型上，随着训练次数不断增加，训练集中检测正确样本数量不断增多，训练集样本总数不变，当训练次数分别达到50、100、200和500次时，根据式（6）计算训练训练准确率，分别为52.38%、53.23%、53.60%和54.03%，保存训练次数为500次的模型对测试集进行测试，测试集样本总数不变，根据式（7）计算得到测试准确率为48.73%。在ResNet18模型上，随着训练次数不断增加，训练集中检测正确样本数量不断增多，训练集样本总数不变，训练准确率逐渐增加，当训练次数分别达到50、100、200 和 500 次时，训练准确率分别为 76.60%、81.96%、83.65%和84.08%，保存训练次数为500次的模型对测试集进行测试，测试准确率为80.34%。在ResNet50 模 型 上，其 训 练 准 确 率 的 变 化 趋 势 与ResNet18模型相似，但最终训练准确率测试准确率较高，分别为 91.72%和 88.37%，这是因为 ResNet50 网络层数多于ResNet18网络的，提取信息量更多。上述3种模型在6种数据集上的训练准确率、测试准确率和运算时间分别见表9、并如图16所示。在准确率方面，VGG16在数据集b、c和d上的训练准确率和测试准确率相同，这是因为数据集中正常样本和异常样本数量不平衡，VGG16模型无法进行有效学习，导致将输入样本都识别为数量占多的样本类型，在数据集a、e和f上的准确率均低于ResNet18图13VGG16在数据集d、e和f的准确率曲线图14ResNet18在数据集d、e和f的准确率曲线图15ResNet50在数据集d、e和f的准确率曲线模型名称VGG16ResNet 18ResNet 50数据集abcdefabcdefabcdef训练集上的准确率/%54.0366.6766.6785.6353.2150.0084.0878.1277.0171.8673.9472.5191.7287.9386.9374.9676.8975.02测试集上的准确率/%48.7366.7666.7685.6344.9045.5180.3474.7972.6663.3358.8765.3888.3778.2677.1265.3754.4661.54运算时间/min141.77108.48108.4850.1791.7316.7341.7033.3533.3516.7733.378.35116.7791.7291.7241.7875.1216.73表93种模型在各数据集上的准确率和运算时间92蔡舒妤等：航空发动机性能图像纹理片段划分方法第 4 期和ResNet50；这2种模型在数据a的训练集和测试集上的准确率都达到最高，其中ResNet50模型的准确率最高，分别为91.72%和88.37%，这是因为数据集a中的正常样本和异常样本数量接近，减少了数据不平衡对模型准确率带来的影响；数据集e训练集和测试集的准确率相差最大，这是因为过采样的方法只是对少数样本进行随机复制以增加样本数量，造成过拟合现象，从而导致这训练集上的准确率远高于测试集；数据集f是通过欠采样方法得到，欠采样方法随机丢弃多数样本以达到数据平衡，这样会丢失有用的信息，导致准确率难以提升。在运算时间上，由于本文所采用的 ResNet中的残差结构减少了运算量，50层的ResNet模型运算时间仍略低于 16 层的 VGG 模型的，且准确率远高于VGG16模型的，ResNet18模型运算时间最短，相较于ResNet50模型的运算时间最多提升65%，但准确率略低于ResNet50模型的。传统的过采样或欠采样方法得到的平衡数据集不仅准确率提升有限而且容易发生过拟合问题，本文提出的方法不但避免了过拟合问题，而且准确率提升明显。4结论（1）在运算速度方面，本文方法能够快速完成性能数据向性能图像的转化，时间复杂度较低。（2）在划分效果方面，本文方法在保留原始数据信息的同时，对性能数据进行图像化降维并精细划分，为建立包含数量相近的正常样本和异常样本的数据集提供了保障。（3）在应用效果方面，利用本文方法划分生成的数据集相较于采样方法得到的数据集在深度学习模型的性能异常检测中效果更好。参考文献：1 周春华，叶子龙，芦吉云，等.基于瞬时能量密度谱的航空发动机双发差异信号监测J.科学技术与工程，2020，20（23）：9655-9659.ZHOU Chunhua，YE Zilong，LU Jiyun，et al.Dual differential signalmonitoring of aeroengine based on instantaneous energy density spectrumJ.Science Technology and Engineering，2020，20（23）：9655-9659.（in Chinese）2 姜广义，武鹏，张洁.航空发动机腔内积油引发振动故障分析J.航空发动机，2013，39（2）：75-78.JIANG Guangyi，WU Peng，ZHANG Jie.Vibration failure caused byaeroengine cavity filled with oilJ.Aeroengine，2013，39（2）：75-78.（in Chinese）3 梁坤，左洪福，孙见忠，等.民机引气系统快速存取记录器数据健康监测方法J.机械工程学报，2015，51（8）：191-197.LIANG Kun，ZUO Hongfu，SUN Jianzhong，et al.Remaining usefullife estimation method for the turbine blade of a civil aircraft enginebased on the QAR and field failure dataJ.Journal of Mechanical Engineering，2015，51（8）：191-197.（in Chinese）4 汪磊，王帅.基于QAR数据和小波变换的飞行操作分析方法J.飞行力学，2020，38（5）：77-83.WANG Lei，WANG Shuai.Flight operation analysis method based onQAR data and wavelet transformJ.Flight Dynamics，2020，38（5）：77-83.（in Chinese）5 陈思，孙有朝，郑敏.基于支持向量机的飞机重着陆风险预警模型J.兵器装备工程学报，2019，40（9）：154-158.CHEN Si，SUN Youzhao，ZHENG Min.Aircraft heavy landing riskearly warning model based on support vector machineJ.Journal ofOrdnance Equipment Engineering，2019，40（9）：154-158.（in Chinese）6 Wang L，Zhang J Y，Dong C T，et al.A method of applying flight data to evaluate landing operation performanceJ.Aeroengine，2019，62（2）：171-180.7 李文朝，林宇震，秦皓.典型工况下低排放燃烧室的压力振荡特性J.航空发动机，2016，42（2）：61-66.LI Wenzhao，LIN Yuzhen，QIN Hao.Pressure oscillation characteristics of low emissions combustor in typical conditionsJ.Aeroengine，2016，42（2）：61-66.（in Chinese）8 曹惠玲，汤鑫豪.基于QAR数据的碳当量值适航符合性验证方法J.北京航空航天大学学报，2020，46（3）：474-480.CAO Huiling，TANG Xinhao.Airworthiness compliance verificationmethod of carbon equivalent based on QAR dataJ.Journal of BeijingUniversity of Aeronautics and Astronautics，2020，46（3）：474-480.（in Chinese）9 魏志强，韩孝兰，胡杨，等.基于航迹数据的空中交通绿色绩效计算分析J.中国环境科学，2019，39（3）：988-993.图163种模型在各数据集上的运算时间对比93航空发动机第 49 卷WEI Zhiqiang，HAN Xiaolan，HU Yang，et al.Calculation and analysis of green performance of air traffic based on track dataJ.China Environmental Science，2019，39（3）：988-993.（in Chinese）10 车畅畅，王华伟，倪晓梅，等.基于深度学习的航空发动机故障融合诊断J.北京航空航天大学学报，2018，44（3）：621-628.CHE Changchang，WANG Weihua，NI Xiaomei，et al.Aeroenginefault fusion diagnosis based on deep learningJ.Journal of BeijingUniversity of Aeronautics and Astronautics，2018，44（3）：621-628.（in Chinese）11 钟诗胜，李旭，张永健.基于DBN的不均衡样本驱动民航发动机故障诊断J.航空动力学报，2019，34（3）：708-716.ZHONG Shisheng，LI Xu，ZHANG Yongjian.Fault diagnosis of civilaviation engine driven by unbalanced ample based on DBNJ.Journal of Aerospace Power，2019，34（3）：708-716.（in Chinese）12 Zhang Y H，Zhou T T，Huang X F，et al.Fault diagnosis of rotatingmachinery based on recurrent neural networksJ.Measurement，2021（1）：171.13 Park S Y，Ahn J.Deep neural network approach for fault detectionand diagnosis during startup transient of liquid-propellant rocket engineJ.Acta Astronautica，2020，177：714-730.14 Guo D F，Zhong M Y，Ji H Q，et al.A hybrid feature model anddeep learning based fault diagnosis for unmanned aerial vehicel sensorsJ.Neurocomputing，2018，319：155-163.15 Chen H T，Jiang B，Zhang T Y，et al.Data-driven and deep learning-based detection and diagnosis of incipient faults with applicationto electrical traction systemsJ.Neurocomputing，2020，396：429-437.16 Manjurul I，Kim J M.Automated bearing fault diagnosis scheme using 2D representation of wavelet packet transform and deep convolutional neural networkJ.Computers in Industry，2019，106：142-153.17 Mateusz B，Atsuto M，Maciej A，et al.A systematic study of the classimbalance problem in convolutional neural networksJ.Neural Networks，2018，106：249-259.18 Guo Q，Li Y，Song Y，et al.Intelligent fault diagnosis method basedon full 1-D convolutional generative adversarial networkJ.IEEETransactions on Industrial Informatics，2020，16（3）：2044-2053.19 Huang N，Chen Q，Cai G