一种基于工况识别和相似性匹配的变工况下航空发动机整机剩余寿命预测方法
【专利摘要】本发明提出了一种基于工况识别和相似性匹配的变工况下航空发动机整机剩余寿命预测方法，该方法着力解决复杂的工况变化掩盖航空发动机真实性能退化趋势致使整机剩余寿命难于预测这个问题，包括如下步骤：首先，识别航空发动机的运行工况；其次，对不同工况下的历史退化数据进行数据标准化；再次，对标准化后的数据进行传感器选择和参数维度约减；接着，使用相似性匹配方法匹配每一台参考发动机和待预测发动机的退化轨迹，获得待预测发动机的剩余寿命估计值和两者之间的相似性距离；最后，通过相似性距离生成加权权重，加权这些寿命估计以获得待预测发动机的剩余寿命。经验证，本发明对不同的发动机测试样本均具有较高的寿命预测精度。
【专利说明】
-种基于工况识别和相似性匹配的变工况下航空发动机整机 剩余寿命预测方法
技术领域
[0001] 本发明设及航空发动机剩余寿命预测的技术领域，具体设及一种基于工况识别和 相似性匹配的变工况下航空发动机整机剩余寿命预测方法。
【背景技术】
[0002] 目前我国航空业正处于蓬勃发展阶段，但依然存在维护和维修成本高等问题。作 为飞机核屯、动力源的航空发动机，是一种安全和可靠性要求极高，设计生产和维修保养花 费极大的大型复杂设备。资料显示，自1963年至1975年期间，美国空军的战斗机共发生了飞 行事故3824起，其中由于发动机原因导致的事故有1664起，占到总数的43.5%;自1994年至 2000年世界航空运输机事故共279,其中因发动机的故障导致飞行事故约20%，因此作为飞 机屯、脏的发动机的维修是航空安全保障的重点。世界范围内航空公司每年大致花费31亿美 元用于飞机维修，其中用于航空发动机的维修费用约占总费用的30%，因此需要不断发展 与改进现有的航空发动机健康管理技术和维修策略，W减少飞行事故的发生并削减维修花 费。随着传感器及监测技术的迅速发展，在航空发动机的使用过程中，积累了大量与寿命信 息相关的性能监测数据，采用运类性能退化数据对高可靠和长寿命设备进行寿命预测成为 发展趋势。
[0003] 当前，国内对航空发动机寿命预测的研究还停留在单个零部件或滑油系统上，针 对航空发动机整机剩余寿命预测的研究基本没有，并且现有研究都是基于单工况假设下进 行的。然而，航空发动机常工作在复杂多变的工况环境下，变化的工况掩盖了系统真实的性 能退化趋势，在此情况下现存的技术方法不能进行准确的剩余寿命预计。
【发明内容】
[0004] 本发明的目的在于:针对变工况下航空发动机整机剩余寿命难于预测的问题，提 出一种基于工况识别和相似性匹配的变工况下航空发动机整机剩余寿命预测方法。
[0005] 本发明采用的技术方案为:一种基于工况识别和相似性匹配的变工况下航空发动 机整机剩余寿命预测方法，其特征在于：
[0006] (1)基于K均值聚类的航空发动机工况识别方法，识别系统在不同时刻所处的工况 后将不同工况下的数据标准化到统一的参数空间；
[0007] (2)为最大限度的挖掘和借鉴参考发动机退化轨迹的衰退规律，通过时间定位和 时间尺度伸缩最优地匹配待预测发动机与参考发动机的性能退化轨迹；
[000引（3)将参考发动机伸缩后轨迹的剩余寿命作为待预测发动机的剩余寿命估计，并 使用相似性距离加权运些估计获得待预测发动机的剩余寿命。
[0009] 本发明与现有技术相比的优点在于：
[0010] (1)、本发明通过工况识别方法解决了航空发动机寿命预测领域目前尚未很好解 决、而实际中常常存在的变工况问题；
[0011] (2)、本发明提出时间定位和时间尺度伸缩来最大化地匹配参考发动机轨迹与待 预测发动机轨迹，提高了历史全寿命数据在预测中的参考价值，算法具有较高的预测精度；
[0012] (3)、本发明借鉴集成学习中的某些思想，对待预测发动机寿命进行多次预测后加 权集成获得待预测的寿命，算法稳健性好，抗干扰能力强。
【附图说明】
[0013] 图1为变工况下航空发动机整机剩余寿命预测流程图；
[0014] 图2为K均值聚类流程图；
[0015] 图3为工况识别流程图；
[0016] 图4为变工况工作示意图；
[0017] 图5为多工况下的退化特性融合示意图，其中，图5(a)为坐标平移示意图，图5(b) 为尺度伸缩示意图；
[0018] 图6为多工况数据标准化流程图；
[0019] 图7为时间定位示意图；
[0020] 图8为伸缩匹配原理图，其中，图8(a)为时间尺度伸缩原理图，图8(b)为最优匹配 示意图；
[0021 ]图9为满轮风扇发动机气路部件简图；
[0022] 图10为工况聚类图；
[0023] 图11为传感器1，2和3标准化前后的退化数据对比（训练集1#发动机），其中，图11 (a)为标准化前的退化数据，图11(b)为标准化之后的退化数据。
[0024] 图12为维度约减后的性能退化轨迹(训练集1#发动机）；
[0025] 图13为测试集1#发动机的3个寿命估计值，其中，图13(a)为使用训练集1#发动机 的估计结果，图13(b)为使用训练集2#发动机的估计结果，图13(c)为使用训练集3#发动机 的估计结果。
【具体实施方式】
[0026] 下面结合附图W及【具体实施方式】进一步说明本发明。
[0027] 如图1所示，基于工况识别和相似性匹配的变工况下航空发动机整机剩余寿命预 测方法，具体步骤如下：
[002引 1工况识别
[0029] 复杂的工况变化常常掩盖航空发动机的性能退化，为还原真实的退化趋势，采用 先辨识各状态的工况，然后进行数据标准化的策略。本发明采用K均值聚类算法对工况参数 进行聚类和分类，从而辨识出系统所处各个状态的运行工况信息。在数据挖掘中，K均值聚 类算法应用广泛，是经典聚类算法之一。
[0030] 设X=陆，X2, . . .，Xn}为已知数据集，X中的Xl，X2, . . .，Xn是n个数据对象并且每个 数据对象都是N维的，即Xi= (Xii，Xi2, ...，XinKK均值聚类算法就是要找到含有K个聚类中 屯、的集合C= (Cl ,〇2 , . . . ,Ck} = { (C11 , C12 , . . . , C1N) , (C21 , C22 , . . . , C2N) , . . . , (CK1 , CK2 ,..., CKN)}使得目标函数：
[0031]
(I)
[0032] 其中，m被归为类Cl的数据对象点数，d(Ci，、)表示聚类中屯、与数据对象的欧几里 德距离，巧巧女化下：
[0033]
(2)
[0034] K均值聚类算法的核屯、思想是把数据集划分成使目标函数达到最小值的K个类。具 体步骤如图2所示。
[0035] 针对训练集数据通过聚类得到各个工况的聚类中屯、和半径。对于测试数据首先计 算实时工况参数与各个聚类中屯、的距离，距离较近者即为该工况。工况识别的流程如图3所 /J、- O
[0036] 2数据标准化
[0037] 复杂系统运行在变工况下时，由于测量尺度的不一致性，状态特征监测数据并不 能真实反映出系统的退化趋势。如图4所示，假设随着系统的退化，系统的某个特征参数是 呈下降趋势的，3条实线(红、黄、蓝)展示了系统在巧中运行工况下不同退化的曲线。情况1曲 线(一种包含多个工况的运行模式)是上下波动的，完全不能反映出系统随着时间的真实退 化趋势;情况2曲线在一定程度上反映了系统的退化趋势，但是是虚假的、不真实的；情况3 曲线表象地完全呈现了与真实退化完全相反的趋势。
[0038] 为实现寿命预测的目的，需要采用数据处理方法，将工况造成的系统监测特征数 据的不同尺度特征去掉，使不同工况特征之间具有可比性、呈现一致的退化特征。直观的， 就是将图4中所示的3种不同工况下的退化曲线合并（融合映射巧Ij 一起，如图5所示。至此， 不同运行工况下的特征数据被转化到了同一个尺度之下，具有切实的可比性，并能反映出 系统的退化趋势。从图形上看，将=种工况合并的两个核屯、步骤是：首先，将=条曲线的重 屯、移到一起(与X轴重合）；其次，对坐标平移后的数据进行尺度缩放，归一化到一个可比的 特征空间。
[0039] 多工况数据标准化方法是使用特征层预测方法预测剩余寿命时的预处理环节。使 用该方法，不同工况模式下的数据的N个特征将被分别标准化到一个标定的范围，然后形成 一个与原始数据特征维数相同的新的腺隹时间序列。图6展示了多工况数据标准化的流程。
[0040] 本发明采用正态标准化(Z-Score)方法将多工况特征数据标准化。令在时刻ti的 工况参数为111，性能状态参数化，合在一起构成系统的状态参数集(1：1，111，义〇4 = 1，2，..，]1。 每个系统状态都将被聚类到P个工况模式Op之一中，因而形成性能状态样本子集lx} W = Ixi UiG〇p}，i = l，2，..，化，p = l，2，. .，P。每个工况模式的样本均值和样本标准差为：
[0041]
(3)
[0042] 对属于工况P的化个状态参数向量进行标准化：
[0043]
(4)
[0044] 式中，np为属于工况P的数据状态数，N为性能参数维度。
[0045] 该标准化方法假定在任一使用循环中，每一个工况发生的概率相等，W确保在同 一运行工况下计算得到的样本均值和标准差是给定工况模式下特征均值和标准差的无偏 估计。最后，合并不同工况模式下的标准化状态参数得到系统退化N维状态时间序列Y = {yi,Y2,. ..,yn}〇
[0046] 3状态参数维度约减
[0047] 高维状态监测参数通常很难处理或处理时计算量较大，实际中通常需要对其进行 数据维度的约减。主成分分析(PCAiPrinciple Component Analysis)是应用最广泛的一种 降维方法之一，是一种统计学方法，在信号处理、模式识别、数字图像处理等领域已经得到 了广泛的应用。主成分分析方法基本思想是提取出空间原始数据中的主要特征（主元），减 少数据冗余，使得数据在一个低维的特征空间被处理，同时保持原始数据的绝大部分的信 息，从而解决数据空间维数过高的瓶颈问题。
[004引设航空发动机退化状态矩阵为X，其列向量Xk=(Xlk，X2k，. . .，xnk)%n维状态参数， 发动机的某一性能状态可由Xk描述，Xk的协方差矩阵为：
[0049]
(5)
[0化0]式中，N为退化状态的采样数，玄为各状态参数的均值：
[0化1]
(6)
[0052]求解Rx的全部特征值、Q = I，2, ...，n)和特征向量Vi,将特征值、按照从大到小 的顺序排列:MH2〉. . .Hn,则对应的特征向量为ViQ = I，2, ..，n)。样本Xi投影到特征向量 Vi得到该方向对应的主分量为：
[005；3]
(7)
[0054]所有的特征向量张成一个n维正交空间，X投影到该正交空间得到相应的n维主分 量。特征向量所对应的特征值越大，它在重构时的贡献也越大，特征值越小的特征向量在重 构时的贡献就越小。设正交空间中前m个主分量为yi，y2,...，ym，其累计方差贡献率为：
[0 化 5]
(8)
[0056]当前面少数几个主分量的累计方差贡献率足够大，如h(m)〉95%，即95% W上的原 始数据信息保留在前面几个主分量中，可只取前m(m<n)个主分量来表征原始信息，在保证 信息完整的情况下，达到数据降维的目的。
[0化7] 4相似性匹配
[005引 （1)时间定位
[0059] 如图7所示，由于不知道待预测发动机轨迹起点的位置，首先需要进行待预测轨迹 的时间定位，找到参考轨迹中与待预测轨迹初始状态参数最接近的时间点。然后将两个轨 迹从该时间点截断，取截断点后的两个轨迹进行后续的进一步匹配计算。设L台参考发动机 的全寿命循环数分别为{1C}，i = l，2, ...，L，待预测发动机的采样点数为no,将待预测发动 轨迹与各参考发动机进行时间定位得到一个集合{Ho}，i = l，2,...，L。
[0060] (2)时间尺度伸缩匹配
[0061] 如图8(a)所示，令离散函数x(t)，t = 0，l, . . .，m和y(t)，t = 0，l, . . .，n分别表示参 考发动机和待预测发动机的退化轨迹。两个轨迹需要满足W下条件：（A)相同的起点；（B)自 变量为正整数。
[0062] 首先，用一个因子A去伸展或压缩x(t)的时间尺度，得到一个临时轨迹x'（t)=x(A t)，t = 0，lA，. . .，mA;AGr。由于临时轨迹的采样点不符合整数采样要求，使用PCHIP插 值算法插值序列x'（t)从而获得一个新轨迹x"(t)，t = 0，l,. . .，[mA]。
[0063] PCHIP算法具体为:令a = xo<:?a<. . .<xn = b为区间I = [a,b]的一个划分，fk，k = 0, 1，. .，n为划分点的函数值，令hk=xk+i-xk，pk = f'（xk)JCHIP函数Ih(X)满足如下限制：
[0064] A.定义域为I = [a，b];
[0065] B.在每个子区间Ik=[祉，祉+1]化=0，1，...，11-1)，14^)可^用一个；次样条多项 式等式Ik(X)表示；
[0066] C.插值条件：Ih(xk)=fk，Ih'（xk)=pk;(k = 0,l,2,...,n)。
[0067] 在每个子区间[Xk，祉+1 ]，PCHIP函数可表示为：
[006引 Ik(X) =fk+Ck,l(X-Xk)+Ck,2(X-Xk)2+Ck,3(X-Xk)3 (9)
[0069] 式中；
[0070]
(10)
[0071] 接着，定义一个相似性距离去量化x"(t)和y(t)之间的相似性：
[0072]
'IP
[0073] 式中，a是一个0-1之间的平滑参数。在本发明中，a = 〇.98。因为近期时间点对未来 退化具有更为重要的影响，因而通过
^距当期时刻较近的时间点分配更大的 权重。[?]表示取整算子。
[0074] 如图8(b)所示，随着伸缩因子A的变化，相似性距离随之变化。显然运是一个凸优 化问题，也就是说存在一个最优的^DPt是相似性距离化(x(t)，y(t))最小。最优函数值被定 义为最优匹配距离Dom。
[00 巧]
U2)
[0076] 式中，DA(x(t)，y(t))为式（11)定义的相似性距离。
[0077] 5集成加权
[007引L台参考发动机在时间定位点之后的寿命分别记为{iM}，1 = 1，2, . . .，L，其为：
[0079] Im=Ic-Ho,1 = 1,2, . . . ,L (13)
[0080] 式中，Ic为第1个参考发动机的全寿命，Ho为待预测发动机与第1个参考发动机的 时间定位位置。
[0081] 通过相似性匹配，得到L个最优伸缩因子和最优匹配距离{iADpt，iD〇M}，l = l， 2,...，L，通过最优伸缩因子计算L个待预测发动机的剩余寿命估计值：
[0082] Ir=I入。Pt ? iM-n〇，l = l，2, . . .，L (14)
[008引式中，no为待预测发动机的数据点数，iAopt为最优匹配因子。
[0084] 至此，得至化个相似性距离和对应的L个抓L估计值，用集合{1d〇m，1R}，1 = 1，2,...， L来表示。为获取待预测发动机的剩余寿命，通过相似性距离生成相似性权重。
r U5)
[0085]
[0086] 式中，Is是第1个参考发动机的相似性的分，其夫
[0087] 加权集成的目的是集成多个寿命估计最终获得待预测发动机的寿命。最简单的集 成方法是加权和，待预测发动机RUL的点估计为：
[008引
（16)
[0089] 实验验证
[0090] 本发明使用美国航空航天局NASA提供的满扇发动机数据集的化taset 2(包括测 试集、训练集）对提出的预测方法进行有效性验证。该数据集由C - M A P S S (C O m m e r C i a 1 Modular Aero-propulsion System Simulation)大型满轮风扇发动机仿真模型得到。模型 通过14个参数的输入来模拟包括风扇、低压压气机、高压压气机、高压满轮、低压满轮(Fan, LPC，HPC，HPT，LPT)在内的五大旋转部件失效影响和退化过程，发动机仿真模型主要部件简 图如图9所示。
[0091] 数据由发动机不同运行循环时间点的3个工况参数(飞行高度、飞行马赫数和油口 杆角度)21个传感器监测的性能参数组成。传感器监测值是针对研究模块，通过发动机热动 力仿真模型获取，且包含噪声。训练数据包括多组独立同一单元的多元时间序列样本，反映 了每个单元从初始到失效的全寿命周期的变化过程。每个部件的初始时刻的退化情况是随 机且未知的，是否发生退化W及退化程度不同，随着发动机的运行，当性能降级到一个阔值 的后整个系统发生失效。用于衡量发动机性能状态的21个监测参数值如表1所示。
[0092] 表1模型输入参数 「nnno1
[0097] I)工况识别
[0098] 使用训练集数据通过K均值聚类得到6个工况的聚类中屯、和半径分别如表3所示， 工况聚类图如图10所示。对于测试数据首先计算实时工况参数与6个聚类中屯、的距离，距离 最近者即判为该工况。
[0099] 表3聚类中屯、和半径
[0100]
[0101] 2)变工况数据标准化
[0102] 即使是同一种状态监测参数，不同工况下的数据不具有可比性，使用本发明提出 的方法对各工况数据进行标准化处理，去除工况的干扰信息，训练集1#发动机的传感器1，2 和3标准化前后的数据如图11所示。
[0103] 从图11可W看出，传感器2和3的监测参数在标准化后呈现出了一定的趋势性，而 传感器1仍然没有明显的退化趋势，说明传感器1完全不能表征任何系统的退化信息。因此， 接下来要进行参数的的选择或降维，去除没有趋势性的传感器参数，并对有趋势的参数进 行维度的约减。
[0104] 3)多元参数维度约减
[0105] 训练集1#发动机状态参数维度约减后的轨迹如图12所示。
[0106] 4)相似性匹配计算
[0107] 用训练集1#、2#、3#发动机对测试集1#发动机进行剩余寿命估计的估计结果如图 13所示，其中测试机1#发动机的真实剩余寿命为18个循环。
[0108] 5)预测结果分析
[0109] 为验证所提出方法更多发动机样本的适用性，使用本发明所提出的方法对测试样 本集中随机选择的10台变工况下运行的航空发动机进行算法的进一步验证，预测结果如表 4所示。
[0110] 表4预测结果 「01111
[0112]从表4可知，本发明提出的航空发动机整机剩余寿命预测方法的平均预测百分比 误差为17%左右，对不同的发动机样本均具有较高的预测精度，提出的预测方法可为航空 发动机的整机剩余寿命预测提供一个新的研究思路。
【主权项】
1. 一种基于工况识别和相似性匹配的变工况下航空发动机整机剩余寿命预测方法，其 特征在于： (1) 基于K均值聚类的航空发动机工况识别方法，识别系统在不同时刻所处的工况后将 不同工况下的数据标准化到统一的参数空间； (2) 为最大限度的挖掘和借鉴参考发动机退化轨迹的衰退规律，通过时间定位和时间 尺度伸缩最优地匹配待预测发动机与参考发动机的性能退化轨迹； (3) 将参考发动机伸缩后轨迹的剩余寿命作为待预测发动机的剩余寿命估计，并使用 相似性距离加权这些估计获得待预测发动机的剩余寿命。
【文档编号】G06F17/50GK106021826SQ201610543483
【公开日】2016年10月12日
【申请日】2016年7月11日
【发明人】刘红梅, 李连峰, 吕琛, 马剑
【申请人】北京航空航天大学