一种基于迁移学习的航空发动机寿命预测与健康评估方法

1.本发明属于复杂设备寿命预测领域，尤其涉及一种基于迁移学习的航空发动机寿命预测与健康评估方法。


背景技术：

2.作为近年来重要的民用出行与军事保障手段，航空飞机以及直升机的使用愈加广泛，由于航空发动机关键限寿件故障或衰退引起的通航一般事故也屡见不鲜。而直升机由于其复杂的动力传输结构和极端的服役环境，发生事故的概率远高于固定翼飞机。因此，开发和应用航空发动机智能运维与健康管理系统(prognostics and health management，phm)对于我国的航空安全事业十分重要。
3.在当前航空发动机phm系统中，关键限寿件的剩余使用寿命(rul)预测与健康状态评估是急需解决的问题；关键限寿件的剩余使用寿命(rul)预测的相关方法可粗略划分为基于模型的方法和基于数据驱动的方法。
4.基于模型的方法主要是依托于系统的物理衰退模型，分析机械设备的动力学特性、机械结构特性以及材料特性。然而，建立一个完备的航空发动机物理模型需要考虑复杂的物理、化学和空气热力学过程，需要耗费巨大的人力资源。相反，基于数据驱动的方法则不再依托于机械知识，根据设备的工作状态参数和环境参数即可评估健康状态和预测寿命。基于数据驱动的方法将粗糙的传感器数据转化为可以学习的有用信息，建立数据与标签的对应关系，使关键限寿件的剩余使用寿命预测与健康状态评估更加容易实现。
5.随着人工智能技术的发展和应用，机器学习方法凭借其灵活高效的优点，逐渐被phm系统研发人员采用。传统机器学习方法极依赖复杂的特征工程技术，而深度学习则免去了这一过程。深度学习依托的神经网络能够自动提取原始数据的深层特征。当前国内航空发动机phm系统的研究主要着眼于航空发动机单个零部件的诊断，对发动机整体性能衰退评估的研究较少。对于航空发动机这类复杂的机械系统，几乎无研究提出一套完整的健康状态划分标准。因此建立健全的健康状态影响因子对于模型评估有着极大的促进作用。


技术实现要素：

6.本发明的目的是克服现有技术中的不足，提供一种基于迁移学习的航空发动机寿命预测与健康评估方法。
7.这种基于迁移学习的航空发动机寿命预测与健康评估方法，包括以下步骤：
8.步骤1、数据采集和数据储存：在数据采集模块内布置传感器，由传感器采集航空发动机关键限寿件的工作状态参数与环境参数；数据存储模块以传感器类别为列，以采集数据时间为行，将采集的工作状态参数与环境参数制成表格存入航空发动机工作状态数据库；
9.步骤2、数据预处理：数据预处理模块将数据采集模块采集得到的工作状态参数与环境参数进行缺失值、异常值及标准化处理；以各个飞机的单次旅程飞行时长作为标准来
划分飞行状态；飞行状态为不同飞行长度发动机的损伤情况；
10.步骤3、健康因子构建：首先选择航空发动机关键限寿件的效率和空气质量流量作为性能指标；然后对性能指标进行归一化处理，选取所有性能指标中的最小值作为健康因子，得到健康因子曲线；最终对健康因子曲线进行平滑处理；
11.步骤4、数据集构建模块采用xgboost模型提取关键性能指标，对关键性能指标进行降维处理；
12.步骤5、将健康因子与标准化处理后的减小速度v取交集，将航空发动机的健康状态划分为初始衰退过程、正常衰退过程和异常衰退过程；减小速度v为当前时刻健康因子曲线的斜率值；建立寿命与健康状态标签，最终对数据切片，制作时序数据集；
13.步骤6、搭建并训练堆叠gru神经网络模型，对时间序列数据进行特征提取和预测；
14.步骤7、利用堆叠gru神经网络模型对多维时序数据的自动特征提取能力挖掘关联寿命与健康状态的有用特征，对航空发动机进行剩余寿命预测与健康评估；
15.步骤8、利用迁移学习策略将堆叠gru神经网络模型泛化至不同的飞行状态下：针对其中一个工况训练堆叠gru神经网络模型，然后将训练完成的堆叠gru神经网络模型的底层神经网络层参数冻结，针对不同工况的数据，对高层网络和输出层进行参数微调，并与未微调堆叠gru神经网络模型的识别准确率对比，检验迁移的成功性。
16.作为优选，步骤1中关键限寿件包括风扇、高压压缩机、低压压缩机、高压涡轮和低压涡轮；工作状态参数与环境参数包括风扇入口压力p1、风扇入口温度t1、风扇出口流量w2、风扇出口压力p2、支路导管压力p3、lpc出口压力p4、lpc出口温度t4、hpc出口压力p5、hpc出口静压力p
s5
、hpc出口温度t5、燃料流量wf、hpt冷却液流量w6、lpt冷却液流量w7、燃烧室压力pb、燃烧室温度tb、lpc入口流量w8、hpc入口流量w9、hpt出口流量w
10
、hpt出口温度t
10
、hpt出口压力p
10
、lpt出口流量w
11
、lpt出口压力p
11
、lpt出口温度t
11
、风扇速度nf、物理核心速度nc、风扇装配余量、lpc装配余量、hpc装配余量、hpc入口燃料流量比例、飞机飞行高度、马赫数和油门旋转变压器角度。
17.作为优选，步骤2中数据预处理模块采用均值补全法对采集得到的工作状态参数与环境参数进行缺失值处理，通过缺失值两侧值的平均数补全空缺；对采集得到的工作状态参数与环境参数进行异常值处理时，将异常值直接舍弃，其中异常值指偏离设备正常工作时状态参数范围的值；进行标准化处理时采用z-score标准化方式，公式为：
[0018][0019]
上式中，μ和σ分别是样本x的均值和方差，x表示进行标准化处理前的数据，x`表示进行标准化处理后的数据。
[0020]
作为优选，步骤3中：
[0021]
采用max-min归一化对性能指标进行归一化处理，计算公式为：
[0022][0023]
上式中，x为样本原值，x`为归一化结果，x
max
和x
min
分别是所有性能指标中的最小值和最大值；
[0024]
采用最小二乘法插值对离散的健康因子进行拟合，根据航空发动机性能衰退关系式对健康因子曲线进行平滑处理，得到最终健康因子平滑曲线；其中航空发动机性能衰退关系式为：
[0025]
δa(t)＝1-exp(atb)+δn(ts)+ξ
[0026]
上式中，δa(t)为t时刻的健康因子，tb表示t时刻的b次方，ts表示航空发动机性投入使用的时刻；δn(ts)为航空发动机初始磨损量；其中a＝u(0.001,0.003)，b＝u(1.4,1.6)，ξ＝n(0,0.001)。
[0027]
作为优选，步骤4中，xgboost模型是由n个基模型组成的加法模型，假设第t次迭代要训练的树模型是f
t
(x)，则第t次迭代时的预测结果满足：
[0028][0029]
上式中：为前t-1棵树的预测结果，f
t
(xi)为第t棵树模型。
[0030]
作为优选，步骤5中：
[0031]
初始衰退过程指航空发动机关键零部件出现的轻度初始磨损的过程；正常衰退过程指由于轻度初始磨损引发的航空发动机性能衰退的过程；异常衰退过程指航空发动机由于关键设备故障引起的性能退化速度加快的过程；
[0032]
对数据切片时，取ts为一个时间步长，分别对不同健康状态的数据进行切片操作，最终得到n
×
ts×
c的输入样本集合，其中n为样本数量，c为样本维度。
[0033]
作为优选，步骤6具体包括以下步骤：
[0034]
步骤6.1、采用python语言搭建堆叠gru神经网络模型：
[0035]
gru神经网络模型包括一个输入层，多个中间层和一个输出层，依次连接输入层、中间层和输出层；
[0036]
输入层特征映射组为多维数组，输入样本格式为n
×
ts×
c；ts为一个时间步长，n为样本数量，c为样本维度；
[0037]
隐含层包含3对堆叠gru-dropout层和一个flatten层；gru神经网络模型的结构单元包括更新门z
t
和重置门r
t
，更新门z
t
用来控制当前状态h
t
从历史状态h
t-1
中保留的信息量和从候选状态h
t`
中接受新信息的量；重置门r
t
用来控制候选状态h
t`
的计算是否依赖历史状态h
t-1
；
[0038]
gru神经网络模型的状态更新方式为：
[0039]ht
＝z
t
⊙ht-1
+(1-z
t
)
⊙ht
`
[0040]
其中更新门输出为：
[0041]zt
＝σ(wzx
t
+u
zht-1
+bz)
[0042]
候选状态h
t`
为
[0043]ht
`＝tanh(whx
t
+uh(r
t
⊙ht-1
)+bh)
[0044]
其中重置门输出为：
[0045]rt
＝σ(wrx
t
+u
rht-1
+br)
[0046]
上式中，x
t
为t时刻gru神经网络模型的输入，wz、uz和bz分别为更新门的输入权重、历史状态权重和偏置；wh、uh和bh分别为候选状态的输入权重、历史状态权重和偏置；wr、ur和br
分别为重置门的输入权重、历史状态权重和偏置；σ和tanh为非线性激活函数；
[0047]
堆叠gru-dropout层对上层的参数进行随机舍弃；最后将堆叠gru-dropout层的输出输入flatten层，进行维数的削减后，转化为一维向量；将一维向量经过输出层输出预测序列结果；
[0048]
步骤6.2、训练堆叠gru神经网络模型，对时间序列数据进行特征提取和预测：
[0049]
将建立的时序数据集输入待训练的gru神经网络模型中，按设定比例划分训练集和测试集；设置单次训练送入数据批量；采用adam算法对堆叠gru神经网络模型进行优化，并设置提前中止命令；当测试集准确率不上升时，停止训练并保存堆叠gru神经网络模型；实时监测堆叠gru神经网络模型的损失函数值变化，最后以折线图的形式输出堆叠gru神经网络模型的预测误差与准确率。
[0050]
作为优选，步骤7中对航空发动机进行剩余寿命预测与健康评估时，采用均方根误差e
rmse
、平均绝对误差e
mape
以及相关系数r2来衡量航空发动机寿命预测的预测值与真实值之间的差异，计算式分别如下：
[0051][0052][0053][0054]
上式中，yi为真实值，为预测值，为均值，n是预测值的个数；
[0055]
进行航空发动机健康评估时，使用准确率来表征gru神经网络模型的分类正确的样本个数占整体样本个数的比例。
[0056]
本发明的有益效果是：首先提出一个健康状态划分因子的构建流程；其次采用受极端偏值影响较小、泛化性更高的极端梯度提升(xgboost)回归模型提取关键性能参数，对数据进行降维；然后利用门控循环单元(gru)网络对多维时序数据的自动特征提取能力挖掘关联寿命与健康状态的有用特征，实现剩余使用寿命预测与健康评估；最终利用迁移学习策略将模型泛化至不同的飞行状态下，实现多工况模型迁移；本发明高效利用航空发动机全寿命周期历史运行数据资源，为航空发动机的寿命预测与健康评估提供可靠的依据。
附图说明
[0057]
图1是航空发动机寿命预测与健康评估的流程图；
[0058]
图2是涡轮增压发动机传感器布置点示意图；
[0059]
图3是健康因子构建过程图；
[0060]
图4是基于xgboost关键性能参数重要性分析结果直方图；
[0061]
图5是健康状态划分过程图；
[0062]
图6是堆叠gru网络结构简图；
[0063]
图7-1为短途飞行中航空发动机剩余使用寿命预测结果图；
[0064]
图7-2为中途飞行中航空发动机剩余使用寿命预测结果图；
[0065]
图7-3为长途飞行中航空发动机剩余使用寿命预测结果图；
[0066]
图8是本发明实施例得到的健康状态评估结果图。
[0067]
附图标记说明：低压压缩机1、高压压缩机2、高压涡轮3、低压涡轮4。
具体实施方式
[0068]
下面结合实施例对本发明做进一步描述。下述实施例的说明只是用于帮助理解本发明。应当指出，对于本技术领域的普通人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干修饰，这些改进和修饰也落入本发明权利要求的保护范围内。
[0069]
实施例一
[0070]
本技术实施例一提供了一种如图1所示基于迁移学习的航空发动机寿命预测与健康评估方法：
[0071]
步骤1、数据采集和数据储存：在数据采集模块内布置传感器，由传感器采集航空发动机关键限寿件的工作状态参数与环境参数；数据存储模块利用sql server数据库技术，以传感器类别为列，以采集数据时间为行，将采集的工作状态参数与环境参数制成表格存入航空发动机工作状态数据库，实现数据的交互和有效存储；航空发动机工作状态数据库一方面与用户和云端进行数据交互，实现接受来自用户的数据、提前为用户缓存数据、向云端上传数据，另一方面存储部分历史数据，为深度学习模块提供训练样本；
[0072]
关键限寿件包括风扇、高压压缩机、低压压缩机、高压涡轮和低压涡轮；工作状态参数与环境参数包括风扇入口压力p1、风扇入口温度t1、风扇出口流量w2、风扇出口压力p2、支路导管压力p3、lpc出口压力p4、lpc出口温度t4、hpc出口压力p5、hpc出口静压力p
s5
、hpc出口温度t5、燃料流量wf、hpt冷却液流量w6、lpt冷却液流量w7、燃烧室压力pb、燃烧室温度tb、lpc入口流量w8、hpc入口流量w9、hpt出口流量w
10
、hpt出口温度t
10
、hpt出口压力p
10
、lpt出口流量w
11
、lpt出口压力p
11
、lpt出口温度t
11
、风扇速度nf、物理核心速度nc、风扇装配余量(smfan)、lpc装配余量(smlpc)、hpc装配余量(smhpc)、hpc入口燃料流量比例(phi)、飞机飞行高度(alt)、马赫数(mach)和油门旋转变压器角度(tra)；
[0073]
步骤2、数据预处理：数据预处理模块将数据采集模块采集得到的工作状态参数与环境参数进行缺失值、异常值及标准化处理；飞行状态是由该航空飞机一次飞行周期的总飞行时长和飞行高度来决定的，又由于数据集中飞行高度与飞行时长存在正相关，本发明以各个飞机的单次旅程飞行时长作为标准来划分飞行状态，从而区别不同飞行里程的发动机之间的损伤情况；飞行状态为不同飞行长度发动机的损伤情况；数据预处理模块采用均值补全法对采集得到的工作状态参数与环境参数进行缺失值处理，通过缺失值两侧值的平均数补全空缺；对采集得到的工作状态参数与环境参数进行异常值处理时，将异常值直接舍弃，其中异常值指偏离设备正常工作时状态参数范围的值；进行标准化处理时采用z-score标准化方式，公式为：
[0074]
[0075]
上式中，μ和σ分别是样本x的均值和方差，x表示进行标准化处理前的数据，x`表示进行标准化处理后的数据；不同损伤模式和飞行状态下的数据均需要进行标准化操作；
[0076]
步骤3、健康因子构建：首先选择航空发动机关键限寿件的效率(e)和空气质量流量(f)作为性能指标，伴随航空发动机的服役时长增加，10个性能指标存在不同速度的退化趋势，其中退化量最大的当前性能指标决定了航空发动机的性能；然后对性能指标进行归一化处理，选取所有性能指标(10个性能指标)中的最小值作为健康因子，得到健康因子曲线；最终对健康因子曲线进行平滑处理；
[0077]
采用max-min归一化对性能指标进行归一化处理，计算公式为：
[0078][0079]
上式中，x为样本原值，x`为归一化结果，x
max
和x
min
分别是所有性能指标中的最小值和最大值；
[0080]
采用最小二乘法插值对离散的健康因子进行拟合，根据航空发动机性能衰退关系式对健康因子曲线进行平滑处理，得到最终健康因子平滑曲线；其中航空发动机性能衰退关系式为：
[0081]
δa(t)＝1-exp(atb)+δn(ts)+ξ
[0082]
上式中，δa(t)为t时刻的健康因子，tb表示t时刻的b次方，ts表示航空发动机性投入使用的时刻；δn(ts)为航空发动机初始磨损量，是发动机出厂即存在的不可避免的轻微损伤；其中a＝u(0.001,0.003)，b＝u(1.4,1.6)，ξ＝n(0,0.001)；
[0083]
步骤4、数据集构建模块采用受极端偏值影响较小、泛化性更高的xgboost模型(极端梯度提升回归模型)提取关键性能指标，对关键性能指标进行降维处理；xgboost模型拟合各个性能参数对健康因子的重要性影响，提取最高的12个关键性能参数作为后续深度学习模型的输入；xgboost模型是一种经优化的分布式梯度提升库，内部决策树采用回归树，具有高效、灵活可移植的特点；其根据样本空间内元素的数据特征对元素进行区间分类，每一次分类均形成树状分枝，若干次迭代后的分枝组合便形成了回归树模型；每个树模型均包含若干个内部节点和叶节点，内部节点将当前空间进行二分，叶节点则是划分后对应的空间结果；
[0084]
xgboost模型是由n个基模型组成的加法模型，假设第t次迭代要训练的树模型是f
t
(x)，则第t次迭代时的预测结果满足：
[0085][0086]
上式中：为前t-1棵树的预测结果，f
t
(xi)为第t棵树模型；
[0087]
步骤5、剩余使用寿命是指航空发动机在经过一次重大检修后经飞机正常飞行若干轮次后至下一次重大检修时所经历的时间周期。当关键限寿件的损伤或出现设备故障时，航空发动机的性能衰退加快，其rul也相对变短。本发明以100轮为重大检修的最长间隔点，当发动机在飞机飞行n(n《100)轮出现故障无法保证飞机正常飞行时，取寿命终点为n，当飞机飞行轮次达到100轮但发动机依旧能够满足飞行要求时也需要检修，此时取100轮为
寿命终点。以起始点寿命记作100％，寿命终止记作0，寿命退化曲线与健康因子成正相关。进一步地，健康状态标签构建方面，将健康因子与标准化处理后的减小速度v取交集，将航空发动机的健康状态划分为初始衰退过程、正常衰退过程和异常衰退过程；减小速度v为当前时刻健康因子曲线的斜率值；当斜率越大，表明健康因子减小速度越快，进而说明航空发动机此时性能衰退的速度也更快；建立寿命与健康状态标签，最终对数据切片，制作时序数据集；初始衰退过程指航空发动机关键零部件由于制造和装配公差出现的不可避免的轻度初始磨损的过程；正常衰退过程指由于轻度初始磨损引发的航空发动机性能衰退的过程，该过程发动机还未收到严重的损伤，处于正常工作状态；异常衰退过程指航空发动机由于关键设备故障引起的性能退化速度加快的过程；
[0088]
对数据切片时，取ts为一个时间步长，分别对不同健康状态的数据进行切片操作，最终得到n
×
ts×
c的输入样本集合，其中n为样本数量，c为样本维度，标签对应时序数据最末端时刻rul和健康状态；
[0089]
步骤6、搭建并训练堆叠gru神经网络模型，对时间序列数据进行特征提取和预测；
[0090]
步骤6.1、采用python语言搭建堆叠gru神经网络模型：
[0091]
gru神经网络模型包括一个输入层，多个中间层(隐含层)和一个输出层(具体网络结构可以根据特定的数据规模进行调整)，依次连接输入层、中间层和输出层；
[0092]
输入层特征映射组为多维数组(维度即为关键性能参数个数)，输入样本格式为n
×
ts×
c；ts为一个时间步长，n为样本数量，c为样本维度；
[0093]
隐含层包含3对堆叠gru-dropout层和一个flatten层；gru是一种自动挖掘时序数据特征，被广泛地用于各类设备phm系统中，例如振动、声学与温度信号的识别；gru相对于传统的lstm网络来说结构更加简化，训练时间较少，避免了lstm内部“门”之间的信息冗余，更有利于处理时间长度大的航空数据集；gru神经网络模型的结构单元包括更新门z
t
和重置门r
t
，更新门z
t
用来控制当前状态h
t
从历史状态h
t-1
中保留的信息量和从候选状态h
t`
中接受新信息的量；重置门r
t
用来控制候选状态h
t`
的计算是否依赖历史状态h
t-1
；
[0094]
gru神经网络模型的状态更新方式为：
[0095]ht
＝z
t
⊙ht-1
+(1-z
t
)
⊙ht
`
[0096]
其中更新门输出为：
[0097]zt
＝σ(wzx
t
+u
zht-1
+bz)
[0098]
候选状态h
t`
为
[0099]ht
`＝tanh(whx
t
+uh(r
t
⊙ht-1
)+bh)
[0100]
其中重置门输出为：
[0101]rt
＝σ(wrx
t
+u
rht-1
+br)
[0102]
上式中，x
t
为t时刻gru神经网络模型的输入，wz、uz和bz分别为更新门的输入权重、历史状态权重和偏置；wh、uh和bh分别为候选状态的输入权重、历史状态权重和偏置；wr、ur和br分别为重置门的输入权重、历史状态权重和偏置；σ和tanh为非线性激活函数；
[0103]
堆叠gru-dropout层对上层的参数进行随机舍弃，从而降低神经网络的复杂度，提高训练效率；最后将堆叠gru-dropout层的输出输入flatten层，进行维数的削减后，转化为一维向量；将一维向量经过输出层输出预测序列结果；输出层为一个全连接层(dense)，对于寿命预测模块，dense层设置1个神经元，不设置非线性激活函数；对于健康评估模块，
dense层设置3个神经元(对应3中健康状态)，非线性激活函数为softmax；dense层增添l1、l2正则化项，优化模型收敛过程，防止过拟合；
[0104]
步骤6.2、训练堆叠gru神经网络模型，对时间序列数据进行特征提取和预测：
[0105]
将建立的时序数据集输入待训练的gru神经网络模型中，按设定比例划分训练集和测试集(一般为4:1，数据量较大时可以适当提高测试集占比)；设置单次训练送入数据批量(32个样本，或32的整数倍)；采用adam算法对堆叠gru神经网络模型进行优化，并设置提前中止命令(early stopping)；当测试集准确率不上升时，停止训练并保存堆叠gru神经网络模型；实时监测堆叠gru神经网络模型的损失函数值变化，最后以折线图的形式输出堆叠gru神经网络模型的预测误差与准确率；
[0106]
步骤7、利用堆叠gru神经网络模型对多维时序数据的自动特征提取能力挖掘关联寿命与健康状态的有用特征，对航空发动机进行剩余寿命预测与健康评估；航空发动机寿命预测属于回归问题，而健康评估属于分类问题，采用均方根误差e
rmse
、平均绝对误差e
mape
以及相关系数r2来衡量航空发动机寿命预测(回归问题)的预测值与真实值之间的差异，计算式分别如下：
[0107][0108][0109][0110]
上式中，yi为真实值，为预测值，为均值，n是预测值的个数；
[0111]
进行航空发动机健康评估(分类问题)时，使用准确率(acc)来表征gru神经网络模型的分类正确的样本个数占整体样本个数的比例。
[0112]
步骤8、迁移学习是指将某一领域(源域)学习到的知识迁移至另一相似领域(目标域)中应用的机器学习策略；本发明采用单一飞行状态下的数据来训练模型，利用迁移学习策略将堆叠gru神经网络模型泛化至不同的飞行状态下，实现多工况下的gru网络模型迁移：将不同的航空发动机运行工况视作不同的领域，针对其中一个工况训练堆叠gru神经网络模型，然后基于模型微调的方法利用其他工况的数据对模型高层神经网络层进行参数微调，进而达到节省计算机资源，提高效率的目的；将训练完成的堆叠gru神经网络模型的底层神经网络层(输入层、部分gru层)参数冻结，针对不同工况的数据，对高层网络和输出层进行参数微调，并与未微调堆叠gru神经网络模型的识别准确率对比，检验迁移的成功性，提高计算效率。
[0113]
实施例二
[0114]
在实施例一的基础上，本技术实施例二提供了实施例一中方法的具体应用：
[0115]
step1：数据采集(本文为涡轮增压发动机全寿命周期运行数据，布设32个传感
器)，并将采集的原始数据存储入sql server数据库；
[0116]
step2：数据预处理，对数据进行异常值处理、缺失值补全，标准化处理，并将数据转换为可用于监督学习的数据类型，并针对不同的飞行参数划分飞行状态，具体标准如表1；
[0117]
表1飞行状态划分依据表
[0118]
飞行状态短途飞行中途飞行长途飞行代号fs1fs2fs3划分标准《3h3h～5h》5h
[0119]
step3：依据健康因子建立标准构建健康因子曲线，表征航空发动机的健康状态退化过程，并利用xgboost模型对性能参数进行重要性分析，提取关键性能参数，划分健康状态并建立寿命标签。
[0120]
step4：将源域数据集(单一工况数据)以4：1的比例划分为训练集，测试集，目标域数据集(其他工况数据)设置为验证集；
[0121]
step5：搭建堆叠gru模型并将训练集数据代入模型进行训练，利用测试集检验误差与准确率变化，保存训练完成的模型；
[0122]
step6：将保存的模型低级网络层锁定，利用验证集小批量数据对高级网络层进行参数微调，达到模型迁移的效果，并验证迁移成功性。
[0123]
如图2所示，为本发明传感器布设点，其中p为总压力，ps为静压力，t为总温度，w为流量，下标1～11分别代表风扇入口、风扇出口、支路导管、lpc出口、hpc出口、hpt冷却液、lpt冷却液、lpc入口、hpc入口、hpt出口、lpt出口，b代表燃烧室，wf为燃料流量，nf为风扇速度，nc为物理核心速度。此外还包括风扇装配余量(smfan)、lpc装配余量(smlpc)、hpc装配余量(smhpc)、hpc入口燃料流量比例(phi)以及飞机飞行高度(alt)、马赫数(mach)和油门旋转变压器角度(tra)一共32个参数。
[0124]
如图3所示，为健康因子构建过程图，本发明以5种关键限寿件的效率(e)和空气质量流量(f)作为性能指标，并归一化，然后取每个时刻的最小值为当前健康因子值，最后根据航空发动机退化关系式对所得健康因子离散值进行拟合为一条光滑曲线。
[0125]
如图4所示，为本发明xgboost性能参数重要性分析结果图，xgboost模型需要设定的主要超参数为最大深度、学习率和参数量。本发明在回归模型建立时的超参数优化方法采取网格化搜索，搜索空间为(最大深度：[5,11]，学习率：[0.01,0.1]，参数量：[100,400])，采用均方根误差作为目标函数，并实施三折交叉验证。最终确定最优超参数为最大深度为8，学习率为0.1，参数量为300。
[0126]
通过xgboost回归树拟合后输出各个元素对健康因子的重要性程度如图4，其中虚线(2.5％)以上的部分表示对健康因子影响较大的性能参数，加上健康因子共13个特征，作为后续寿命预测与健康评估的输入变量。
[0127]
如图5所示，为健康状态划分过程图，根据健康因子的衰减速度可以将航空发动机的健康状态划分为初始衰退、正常衰退与异常衰退过程。初始衰退是指发动机关键零部件由于制造和装配公差引起的不可避免的轻度初始磨损量；正常衰退是指由于初始磨损量引发的发动机性能衰退的过程，该过程发动机还未收到严重的损伤，处于正常工作状态；而异常衰退则是指发动机由于关键设备故障引起的性能退化速度加快的阶段。本发明提出将健
fs3，fs2-fs1，fs2-fs3，fs3-fs1，fs3-fs2，前者为源域，后者为目标域。采用迁移策略和未采用迁移策略的结果对比如下表4。结果显示，在跨工况预测和识别方面，未采用迁移学习策略前，单一工况模型对其他工况数据的识别效果很差，甚至存在低于60％，采用了迁移学习后误差明显降低，准确率明显上升，且相邻两个工况由于数据相似度较高，识别效果也更好。这表明，在实际情况下数据稀缺或分布不均衡的情况下采用迁移学习策略对模型进行更新，能够明显地提升预测和识别效果。
[0138]
表4迁移效果比较表
[0139]