1.本发明属于航空发动机性能退化技术领域，具体涉及一种基于多源传感器状态距离的航空发动机性能退化状态评价方法。
背景技术：2.航空发动机退化分析是航空工业领域研究的重点之一，如何有效分析发动机的传感器退化数据并得到描述发动机性能退化的综合指标是退化状态评价的关键。随着传感器技术和大数据的发展，基于数据驱动的航空发动机退化评价方法成为重要趋势。基于数据驱动的方法重点在于如何将退化过程中的传感器监测数据有效利用起来以准确描述发动机的退化过程。
3.目前，现有的方法采用与退化过程最相关的传感器数据近似描述发动机的退化过程。这种简化的方法虽然可以描述退化的整体趋势，但是单一传感器所包含的信息较少，这些信息不足以准确描述退化过程。随着数据融合技术的发展，许多研究开始通过融合退化过程中的多传感器数据得到一个综合指标描述发动机退化状态。融合多传感器数据的方法与单一传感器数据方法相比较有很大的提升，但是目前大多数研究都没有考虑退化过程中传感器的重要性和传感器之间的相关性。在退化过程的多个传感器中，一部分传感器对性能退化影响很大，另一部分传感器对性能退化过程影响很小。此外，部分传感器之间还存在某些关系，例如温度和压力等。因此，航空发动机的传感器退化监测数据具有不同的特点和相关性。为了解决上述问题，提出了一种新的描述发动机性能退化状态的方法，即基于多源传感器状态距离的性能退化状态评价方法。
技术实现要素：4.本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种基于多源传感器状态距离的性能退化状态评价方法。通过多源传感器状态距离的方法，将退化过程中的多个传感器数据融合为一个综合变量表征退化过程，提出的方法充分考虑了发动机退化过程中不同传感器对退化过程影响不同、传感器之间的关联性以及初始运行条件不同等问题，所以可以更加准确描述退化过程。
5.为解决上述技术问题，本发明采用的技术方案是：一种基于多源传感器状态距离的性能退化状态评价方法，其特征在于，包括以下步骤：
6.步骤一、数据预处理：
7.步骤101：整个退化过程的时刻为n，输入完整退化过程中m个传感器退化序列数据x
o1
＝[x1,x2,
…
,xm]，第j个传感器的退化序列表示为剔除性能退化过程中的无关传感器数据，最终保留k个传感器退化数据x
o2
＝[x1,x2,
…
,xk]，k≤m；
[0008]
步骤102：对输入的k个传感器退化数据进行归一化处理，根据公式
计算各个传感器归一化后的数据，其中x表示归一化前的数据，x' 表示归一化后的数据；第j个传感器归一化后的数据为j个传感器归一化后的数据为
[0009]
步骤二、计算传感器的权重：
[0010]
步骤201：根据公式计算整个退化过程中每个传感器的单调性指标；
[0011]
步骤202：根据公式计算整个退化过程中每个传感器的预测性指标，其中初始时刻的均值；为失效时刻的均值；σ(xf)为失效时刻的标准差；
[0012]
步骤203：根据公式计算整个退化过程中每个传感器的相关性指标，其中ti表示退化过程的时间序列， i＝1,2,
…
,n；
[0013]
步骤204：将单调性、预测性和相关性三个指标组成向量(mon,pre,corr) 与最优向量(1,1,1)计算欧式距离d，d描述了多评价指标的坐标向量与最优点之间的绝对差异；
[0014]
步骤205：根据公式将绝对差异转化为相似度s＝[s1,s2,
…
,sk]，再根据公式归一化相似度得到传感器的权重w＝[w1,w2,
…
,wk]；
[0015]
步骤三、计算多源传感器状态距离表征退化状态：
[0016]
步骤301：生成权重矩阵，将权重w＝[w1,w2,
…
,wk]转化为对角线权重矩阵
[0017]
步骤302：根据公式计算传感器数据之间的相关系数，其中x和y表示不同的传感器数据，生成相关系数矩阵
[0018]
步骤303：根据公式计算性能退化状态，其中dii表示第i个时刻的性能退化状态指标；表示第i个时刻的传感器数据向量；表示健康状态的传感器数据向量；表示失效状态的传感器数据向量；t 表示向量转置；w和r表示权重矩阵和相关系数矩阵；
[0019]
步骤四、计算退化过程的近似数学模型：
[0020]
步骤401：选择非线性维纳过程描述退化过程，设定的方程为 x(t)＝x(0)+ae
bt
+σb(t)其中x表示退化状态指标，b(t)表示标准布朗运动，a、b 和σ表示求解的参数，令θ＝{a,b,σ}；
[0021]
步骤402：设发动机从健康状态运行到完全失效状态共得到n个时刻表示退化性能状态的指标di＝{x1,x2,
…
,xn},退化过程中的时间序列为t＝{t1,t2,
…
,tn}，退化指标的增量数据集为δdi＝{δx1,δx2,
…
,δxn}，则根据非线性维纳过程的性质可以计算得到关于θ＝{a,b,σ}的似然函数为
[0022]
步骤403：采用多维搜索优化的方法求解参数的极大似然估计值。
[0023]
本发明与现有技术相比具有以下优点：
[0024]
1、本发明提出的基于单调性、预测性和相关性的传感器特征加权方法，可以更合理地计算退化过程中各个传感器的权重。
[0025]
2、本发明提出的基于多源传感器状态距离的航空发动机性能退化状态评价方法，充分考虑了发动机退化过程中不同传感器对退化过程影响不同、传感器之间的关联性以及初始运行条件不同等问题，可以更合理地融合退化过程的传感器数据以表征退化过程。
[0026]
综上所述，本发明通过基于多源传感器状态距离的方法，将退化过程中的多个传感器数据融合为一个综合变量表征退化过程，提出的方法充分考虑了发动机退化过程中不同传感器对退化过程影响不同、传感器之间的关联性以及初始运行条件不同等问题，更加准确地描述退化过程。
[0027]
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
[0028]
图1为本发明的流程框图
[0029]
图2为100台发动机退化状态指标图
[0030]
图3为5种方法生成的退化状态指标
[0031]
图4为计算的近似数学模型
[0032]
图5为15个传感器的单调性、预测性和相关性指标计算结果
[0033]
图6为15个传感器计算的权重结果
[0034]
图7为100台涡扇发动机5种不同方法的评价标准平均值
[0035]
图8为100台涡扇发动机5种不同评价状态退化方法的在翼时间误差平均值。
具体实施方式
[0036]
下面结合实施例对本发明的方法作进一步详细的说明。
[0037]
如图1所示，本发明包括以下步骤：
[0038]
步骤一、数据预处理：
[0039]
为了验证提出了方法的有效性，采用nasa公开的c-mapss数据集进行验证。c-mapss数据集包含4个子数据集，每个子集包含不同的工况数和故障数，采用单工况和单故障模式条件下的fd001子集展示实施过程。
[0040]
步骤101：fd001数据集包含100台发动机的完整退化过程传感器数据，其中包含3个工况操作数据，21维传感器监测数据。剔除无关传感器数据和 3个工况操作参数数据，最终保留15维传感器监测数据；
[0041]
步骤102：对最终保留的15传感器数据进行归一化处理，根据公式计算各个传感器归一化后的数据x
nor
＝[x1,x2,
…
,xk]，其中第j 个传感器归一化后的数据为
[0042]
步骤二、计算传感器的权重：
[0043]
步骤201：根据公式计算整个退化过程中每个传感器的单调性指标。实际使用过程中，分别计算15维传感器数据的单调性指标值。
[0044]
步骤202：根据公式计算整个退化过程中每个传感器的预测性指标，其中为初始时刻的均值；为失效时刻的均值；σ(xf)为失效时刻的标准差。实际使用时，选择前5％的初始时刻的数据计算初始时刻的均值选择后5％的退化末期的数据计算失效时刻的均值分别计算 15维传感器数据的预测性指标值。
[0045]
步骤203：根据公式计算整个退化过程中每个传感器的相关性指标，其中ti表示退化过程的时间序列， i＝1,2,
…
,n。实际使用过程中，分别计算15维传感器数据的相关性指标值。
[0046]
步骤204：将单调性、预测性和相关性三个指标组成向量(mon,pre,corr) 与最优向量(1,1,1)计算欧式距离，所述的计算公式为d描述了多评价指标的坐标向量与最优点之间的绝对差异。实际使用过程中，图5展示了15维传感器的单调性、预测性和相关性三个指标
组成向量(mon,pre,corr)。
[0047]
步骤205：根据公式将绝对差异转化为相似度s＝[s1,s2,
…
,sk]，再根据公式归一化相似度得到传感器的权重w＝[w1,w2,
…
,wk]。实际使用过程中，图6展示了15维传感器计算的权重结果。
[0048]
步骤三、计算多源传感器状态距离表征退化状态：
[0049]
步骤301：生成权重矩阵，将权重w＝[w1,w2,
…
,wk]转化为对角线权重矩阵
[0050]
步骤302：根据公式计算传感器数据之间的相关系数，其中x和y表示不同的传感器数据，生成相关系数矩阵
[0051]
步骤303：根据公式计算性能退化状态，其中dii表示第i个时刻的性能退化状态指标；表示第i个时刻的传感器数据向量；表示健康状态的传感器数据向量；表示失效状态的传感器数据向量；t 表示向量转置；w和r表示权重矩阵和相关系数矩阵；
[0052]
实际使用过程中，计算fd001数据集下100台发动机的退化状态如图2 所示，图中横轴表示时刻，纵轴表示退化状态值。从图2可以看出，基于多源传感器状态距离得出的100台涡扇发动机的退化状态退化趋势基本保持一致，每条退化状态曲线在退化初期变化缓慢，随着运行时间的增加，发动机的退化状态逐渐加快，越靠近退化的末期，退化趋势变化越剧烈。
[0053]
为了证明方法的有效性，选取与退化最相关的单一传感器直接作为性能退化评价方法m1、由基于欧式距离构建退化性能评价方法m2、基于马氏距离构建退化性能方法m3和基于粗糙核距离构建性能退化方法m4这四种方法进行对比。图3展示了fd001下一台发动机
通过5种方法得出的退化状态曲线，从图3中可以看出提出的方法提升效果很大，可以更准确的描述发动机的状态退化过程。图7展示了100台涡扇发动机初始健康状态周期个数设定不同下5种不同退化状态评价方法计算得出的预测性，单调性以及相关性评价标准的平均值，同样证明提出方法的有效性。
[0054]
步骤四、计算退化过程的近似数学模型：
[0055]
步骤401：选择非线性维纳过程描述退化过程，设定的方程为 x(t)＝x(0)+ae
bt
+σb(t)其中x表示退化状态指标，b(t)表示标准布朗运动，a、b 和σ表示求解的参数，令θ＝{a,b,σ}；
[0056]
步骤402：设发动机从健康状态运行到完全失效状态共得到n个时刻表示退化性能状态的指标di＝{x1,x2,
…
,xn},退化过程中的时间序列为 t＝{t1,t2,
…
,tn}，退化指标的增量数据集为δdi＝{δx1,δx2,
…
,δxn}，则根据非线性维纳过程的性质可以计算得到关于θ＝{a,b,σ}的似然函数为
[0057]
步骤403：采用多维搜索优化的方法求解参数的极大似然估计值；
[0058]
实际使用过程中，通过非线性维纳过程计算的近似数学模型如图4所示，图中横轴表示时刻，纵轴表示退化状态值，di表示退化状态值。为了进一步验证方法的有效性，通过非线性维纳过程对所提出的方法计算的退化状态指标以及四种对比方法进行方程建模，用以表征退化状态。如图4所示，将每种方法计算的退化状态指标，代入到对应的非线性维纳方程中，计算每个退化状态所对应方程预测的在翼时间与真实的在翼时间的误差。图8展示了100 台涡扇发动机在5种不同的方法建立退化状态指标后在翼时间预测的误差的平均值，同样可以证明提出的方法的有效性。
[0059]
以上所述，仅是本发明的实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技术方案的保护范围内。
技术特征：1.一种基于多源传感器状态距离的航空发动机性能退化状态评价方法，其特征在于，包括以下步骤：步骤一、数据预处理：步骤101：整个退化过程的时刻为n，输入完整退化过程中m个传感器退化序列数据x
o1
＝[x1,x2,
…
,x
m
]，第j个传感器的退化序列表示为剔除性能退化过程中的无关传感器数据，最终保留k个传感器退化数据x
o2
＝[x1,x2,
…
,x
k
]，k≤m；步骤102：对输入的k个传感器退化数据进行归一化处理，根据公式计算各个传感器归一化后的数据，其中x表示归一化前的数据，x'表示归一化后的数据；第j个传感器归一化后的数据为j＝1,2,
…
,k；步骤二、计算传感器的权重：步骤201：根据公式计算整个退化过程中每个传感器的单调性指标；步骤202：根据公式计算整个退化过程中每个传感器的预测性指标，其中为初始时刻的均值；为失效时刻的均值；σ(x
f
)为失效时刻的标准差；步骤203：根据公式计算整个退化过程中每个传感器的相关性指标，其中t
i
表示退化过程的时间序列，i＝1,2,
…
,n；步骤204：将单调性、预测性和相关性三个指标组成向量(mon,pre,corr)与最优向量(1,1,1)计算欧式距离d，d描述了多评价指标的坐标向量与最优点之间的绝对差异；步骤205：根据公式将绝对差异转化为相似度s＝[s1,s2,
…
,s
k
]，再根据公式归一化相似度得到传感器的权重w＝[w1,w2,
…
,w
k
]；步骤三、计算多源传感器状态距离表征退化状态：步骤301：生成权重矩阵，将权重w＝[w1,w2,
…
,w
k
]转化为对角线权重矩阵
步骤302：根据公式计算传感器数据之间的相关系数，其中x和y表示不同的传感器数据，生成相关系数矩阵步骤303：根据公式计算性能退化状态，其中di
i
表示第i个时刻的性能退化状态指标；表示第i个时刻的传感器数据向量；表示健康状态的传感器数据向量；表示失效状态的传感器数据向量；t表示向量转置；w和r表示权重矩阵和相关系数矩阵；步骤四、计算退化过程的近似数学模型：步骤401：选择非线性维纳过程描述退化过程，设定的方程为x(t)＝x(0)+ae
bt
+σb(t)其中x表示退化状态指标，b(t)表示标准布朗运动，a、b和σ表示求解的参数，令θ＝{a,b,σ}；步骤402：设发动机从健康状态运行到完全失效状态共得到n个时刻表示退化性能状态的指标di＝{x1,x2,
…
,x
n
},退化过程中的时间序列为t＝{t1,t2,
…
,t
n
}，退化指标的增量数据集为δdi＝{δx1,δx2,
…
,δx
n
}，则根据非线性维纳过程的性质可以计算得到关于θ＝{a,b,σ}的似然函数为步骤403：采用多维搜索优化的方法求解参数的极大似然估计值。2.所述权利要求1步骤204中计算欧式距离的公式为
技术总结本发明公开了一种基于多源传感器状态距离的航空发动机性能退化状态评价方法，包括以下步骤：步骤一、数据预处理；步骤二、计算传感器的权重；步骤三、计算多源传感器状态距离表征退化状态；步骤四、计算退化过程的近似数学模型；本方法提出了一种基于多评价指标的传感器特征加权的方法，充分利用了多源传感器退化数据，同时提出了一种基于多源传感器状态距离的航空发动机性能退化状态评价方法，充分考虑了传感器之间的关联性以及初始运行条件不同等问题，因此可以更加准确描述航空发动机退化过程。过程。过程。
技术研发人员：邓鑫洋 李新宇 蒋雯 耿杰
受保护的技术使用者：西北工业大学
技术研发日：2022.01.22
技术公布日：2022/7/5
