CN108363844A - 一种航空发动机启动过程排气温度预测方法 - Google Patents
一种航空发动机启动过程排气温度预测方法 Download PDFInfo
- Publication number
- CN108363844A CN108363844A CN201810075177.6A CN201810075177A CN108363844A CN 108363844 A CN108363844 A CN 108363844A CN 201810075177 A CN201810075177 A CN 201810075177A CN 108363844 A CN108363844 A CN 108363844A
- Authority
- CN
- China
- Prior art keywords
- data
- aero
- engine
- prediction
- para
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Abstract
本发明属于航空发动机预测技术领域,提供了一种航空发动机启动过程排气温度预测方法。本发明所要解决的技术问题是填补航空发动机在启动过程中排气温度预测的空缺,提供一种航空发动机启动过程排气温度预测方法,利用航空发动机地面试车数据,采用机器学习的方法得到发动机启动过程排气温度预测模型,该模型预测精度高、泛化性能好,预测结果可进一步用于发动机控制方面等,降低了发动机出现超温的可能性。相对于传统的单参数预测,本发明由于采用融合预测,包含了更多信息,使得预测误差降低;相对于单一的预测算法,本发明由于采用了AdaBoost.RT集成算法,对弱学习机进行集成,预测误差更小。
Description
技术领域
本发明属于航空发动机预测技术领域,具体涉及一种航空发动机启动过程排气温度预测方法。
背景技术
在飞机启动时,航空发动机处于高温、高负荷和高转速的状态,发动机出现超温的可能性较大,增加了飞行风险,因此,需要对排气温度进行预测以及时对航空发动机进行控制防止出现超温。对航空发动机排气温度预测的方法主要有三种:基于模型的方法、基于回归的方法和基于机器学习的方法。其中,基于模型的方法计算复杂、实时计算可能出现迭代不收敛等问题;基于回归的方法有时变量之间不一定有明显的线性或者其他函数关系,模型很难选择;基于机器学习的方法具有非常强的非线性映射能力,训练时间短。在文献《基于支持过程向量机的航空发动机排气温度预测》中,于广斌等人提出一种支持过程向量机模型,并运用到航空发动机排气温度预测中以预测航空发动机气路性能衰退规律,预测精度高。在文献《Application of Neural Networks in Forecasting Engine SystemsReliability》中,Xu K等人利用神经网络对航空发动机排气温度进行预测以预测发动机系统故障和可靠性。以上方法都是以航空发动机多个飞行循环的排气温度为数据进行预测,以显示航空发动机的性能状态及退化情况,但没有涉及到整个启动过程的排气温度预测,因此无法在有超温现象前对发动机进行控制。
发明内容
本发明所要解决的技术问题是填补航空发动机在启动过程中排气温度预测的空缺,提供一种航空发动机启动过程排气温度预测方法,利用航空发动机地面试车数据,采用机器学习的方法得到发动机启动过程排气温度预测模型,该模型预测精度高、泛化性能好,预测结果可进一步用于发动机控制方面等。
按照本发明,提供了一种航空发动机启动过程排气温度预测方法,本发明的技术方案如下:
首先,对传感器采集到的航空发动机地面试车数据,如高压压气机转速、低压压气机转速、主燃油流量、低压涡轮后温度等进行预处理,主要包括对异常数据的识别和处理、数据的平滑处理及数据的归一化处理。然后,基于信息融合的思想,通过合适的相关性方法选择与排气温度相关性大的参数作为输入参数对其进行预测。另外,对选择的参数进行相空间重构构建输入输出样本。最后,采用机器学习的算法对排气温度进行预测,获得预测精度高、泛化能力强、鲁棒性好的航空发动机启动过程排气温度预测模型。
优选地,在本发明中,对异常数据的识别采用基于密度的方法然后进行剔除。对数据的平滑处理采用专用函数平滑法。相关性分析方法采用互信息法。相空间重构的参数分别采用互信息法和Cao法。机器学习算法采用AdaBoost.RT集成算法,通过对弱学习机极限学习机(ELM)进行集成,得到效果优越的强学习机。
本发明的有益效果:本发明的预测模型预测精度高、泛化能力强、鲁棒性强,可对航空发动机启动过程的排气温度进行实时预测,预测结果可用于航空发动机控制等方面,降低了发动机出现超温的可能性。相对于传统的单参数预测,本发明由于采用融合预测,包含了更多信息,使得预测误差降低;相对于单一的预测算法,本发明由于采用了AdaBoost.RT集成算法,对弱学习机进行集成,预测误差更小。
附图说明
图1是本发明的流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及技术方案,对本发明进行进一步详细说明。
一、航空发动机启动过程地面试车数据预处理
设传感器采集到的航空发动机启动过程地面试车数据Data为
Data=[Para1,Para2,…,Paral,…,ParaN] (1)
其中Para为高压压气机转速、低压压气机转速、主燃油流量、低压涡轮后温度等航空发动机性能参数数据,为相对应的时间序列,N为参数个数,n为样本个数。航空发动机启动过程地面试车数据预处理包括异常点的识别与处理、数据平滑处理及数据归一化处理。
1、基于密度的异常点识别
在时间序列中,点xli临近点的数量越少,即其周围点密度较小,则点xli就更加可能为一异常点。对于时间序列中任意两点形成的有效点对(xli,xlj),其欧几里得距离可表示为:
对有效点对(xli,xlj)中的点xli,定义其k近邻点距离(k>0,k∈N)为k-dist(xli),且k-dist(xli)满足:
(1)在中,满足dist(xli,xlj)≤k-dist(xli)的数据点的数量至少为k个;(2)在中,满足dist(xli,xlj)<k-dist(xli)的数据点的数量至多为k-1个。
对有效点对(xli,xlj)中的点xli,称
r-distk(xli,xlj)=max(dist(xli,xlj),k-dist(xli)) (4)
为有效点xlj与xli的k近邻点极限距离。
为衡量xli点周围点数量的多少,定义局部极限密度这一概念,称
为点xlj的k局部极限密度,其中K(xlj)为点xlj的k近邻点集合。
定义为点xlj的异常系数,若λ<ρ,则点xlj为正常数据点,若λ≥ρ,则点xlj为时间序列的一个异常点,其中ρ为异常因子上限。
对识别出的异常点进行剔除,并用左右两端相邻数据值的均值填充。
2、基于二次指数平滑法的数据平滑处理
采用二次指数平滑法对异常点识别及处理后的数据Data1进行处理,以去除在信号采集过程中可能出现的噪声或数据污染,其中
Data1=[Para1 1,Para2 1,…,Paral 1,…,ParaN 1] (6)
Para1为Para异常点识别及处理后的航空发动机性能参数数据,为异常点识别及处理后的时间序列。
二次指数平滑算法如下:
其中,α为平滑系数;Sli (1),Sli (2)分别为一次和二次平滑值,定义平滑初始值Sl0为
3、数据归一化处理
对平滑处理后的数据Data2进行归一化处理,转化为[0,1]范围内的数据,其中
Data2=[Para1 2,Para2 2,…,Paral 2,…,ParaN 2] (10)
Para2为Para1平滑处理后的航空发动机性能参数数据,为平滑处理后的时间序列。
二、航空发动机启动过程地面试车数据相关性分析
对预处理后的数据Data3采用互信息法进行相关性分析,其中
Data3=[Para1 3,Para2 3,…,Paral 3,…,ParaN 3] (12)
Para3为Para2归一化处理后的航空发动机性能参数数据,为归一化处理后的时间序列。
设有两组航空发动机性能参数数据和其中p,q∈N且p≠q。xpi 3和xqj 3的概率密度分别为Pxp[xpi 3]和Pxq[xqj 3],联合概率密度为Pxpq[xpi 3,xqj 3],则互信息函数MI(xp,xq)为
MI(xp,xq)=H(xp)+H(xq)-H(xp,xq) (14)
计算各个参数与排气温度的互信息函数值,考虑到各个参数与排气温度相关性大小的差异及训练预测模型的时间要求,取互信息函数值最大的三个参数及排气温度(EGT)作为预测模型的输入参数。
三、航空发动机启动过程地面试车数据相空间重构
因为航空发动机启动过程地面试车数据是一组时间序列数据,为把其隐含的信息充分显露出来,对一维时间序列数据进行相空间重构。对相关性分析后的数据Data4进行相空间重构,其中
Data4=[Para1 3,Para2 3,Para3 3,Para4 3] (18)
特别地,Para4 3=EGT3。对于时间序列对其重构后的相空间为
Xl=[Xl1,Xl2,…,XlI,…,XlM]T (20)
其中,
XlI=[xlI,xl(I+τ),…,xl(I+(m-1)τ)],I=1,2,…,M;M=n-(m-1)τ (21)
m为嵌入维数,τ是延迟时间,分别用互信息法和Cao法进行求解。根据相空间重构构造输入输出样本,如表1所示,其中h为预测步长。
表1基于相空间重构的输入输出数据
样本个数 | 输入数据X | 输出数据Y |
1 | [X11,X21,X31,X41] | EGT1+(m-1)τ+h 3 |
… | … | … |
I | [X1I,X2I,X3I,X4I] | EGTI+(m-1)τ+h 3 |
… | … | … |
M | [X1M,X2M,X3M,X4M] | EGTM+(m-1)τ+h 3 |
四、航空发动机启动过程排气温度预测模型
本发明采用AdaBoost.RT_ELM算法对航空发动机启动过程的排气温度进行预测。AdaBoost.RT_ELM具体算法如下:
(1)输入
相空间重构后的输入输出数据
选择弱学习算法
指定迭代次数T(也表示最终生成弱学习机的数目);
指定相对误差绝对值的阈值φ,在训练时根据φ将训练样本分为预测正确的样本和预测错误的样本。
(2)初始化
令初始迭代次数t=1;
令第一次训练时,训练样本权值分布为Dt(I)=1/M,I=1,…,M;
令初始误差率εt=0。
(3)迭代过程
for t=1,…,T:
Step 1:在权值为Dt的训练样本上训练第t个弱学习机;
Step 2:记第t个学习机ft对第I个样本XI的预测结果为ft(XI),而实际真值为YI。计算ft的误差率:
Step 3:设置a可以为1,2或3;
Step 4:更新样本权值Dt:
其中Zt是标准化因子。通过调整各样本的权值,即增加预测误差大的样本的权值,减少预测误差小的样本的权值,使得误差大的样本在下一次迭代中更受关注。
end
(4)输出
强学习机:
本发明以学习速度快、泛化性能好的极限学习机(ELM)为弱学习机,通过设定合适的迭代次数以及阈值,得到预测精度高的强学习机,即航空发动机启动过程排气温度预测模型。
Claims (1)
1.一种航空发动机启动过程排气温度预测方法,其特征在于,对传感器采集到的航空发动机启动过程地面试车数据,采用基于密度的方法对数据进行异常点的识别和处理,采用二次指数平滑法对数据中的噪声或数据污染进行平滑或滤波处理,并对数据进行归一化处理使其转化为[0,1]范围内的数据;
基于信息融合的思想,利用互信息法进行相关性分析,计算各个参数与排气温度的互信息函数值,考虑到各个参数与排气温度相关性大小的差异及训练预测模型的时间要求,取互信息函数值最大的三个参数及排气温度作为预测模型的输入参数;设预处理后的数据为Data,其中
Data=[Para1,Para2,…,Paral,…,ParaN] (1)
其中,Para航空发动机性能参数数据,为相对应的时间序列,N为参数个数,n为样本个数;
设有两组航空发动机性能参数数据和其中p,q∈N且p≠q;xpi和xqj的概率密度分别为Pxp[xpi]和Pxq[xqj],联合概率密度为Pxpq[xpi,xqj],则互信息函数MI(xp,xq)为
MI(xp,xq)=H(xp)+H(xq)-H(xp,xq) (3)
对选择的参数进行相空间重构构建输入输出样本以充分显露时间序列数据中隐含的信息;设相关性分析后的数据为Data1,其中,
Data1=[Para1,Para2,Para3,Para4] (7)
特别地,Para4=EGT;对于时间序列对其重构后的相空间为
Xl=[Xl1,Xl2,…,XlI,…,XlM]T (9)
其中,
XlI=[xlI,xl(I+τ),…,xl(I+(m-1)τ)],I=1,2,…,M;M=n-(m-1)τ (10)
其中,m为嵌入维数,τ是延迟时间,分别用互信息法和Cao法进行求解,根据相空间重构构造输入输出样本,如表1所示,其中h为预测步长;
表1 基于相空间重构的输入输出数据
采用AdaBoost.RT_ELM算法对航空发动机启动过程的排气温度进行预测,AdaBoost.RT_ELM具体算法如下:
(1)输入
相空间重构后的输入输出数据
选择弱学习算法
指定迭代次数T;
指定相对误差绝对值的阈值φ,在训练时根据φ将训练样本分为预测正确的样本和预测错误的样本;
(2)初始化
令初始迭代次数t=1;
令第一次训练时,训练样本权值分布为Dt(I)=1/M,I=1,…,M;
令初始误差率εt=0;
(3)迭代过程
for t=1,…,T;
Step 1:在权值为Dt的训练样本上训练第t个弱学习机;
Step 2:记第t个学习机ft对第I个样本XI的预测结果为ft(XI),而实际真值为YI;计算ft的误差率:
Step 3:设置a可以为1,2或3;
Step 4:更新样本权值Dt:
其中,Zt是标准化因子;通过调整各样本的权值,即增加预测误差大的样本的权值,减少预测误差小的样本的权值,使得误差大的样本在下一次迭代中更受关注;
end
(4)输出
强学习机:
以学习速度快、泛化性能好的极限学习机为弱学习机,通过设定合适的迭代次数以及阈值,得到预测精度高的强学习机,即航空发动机启动过程排气温度预测模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810075177.6A CN108363844B (zh) | 2018-01-26 | 2018-01-26 | 一种航空发动机启动过程排气温度预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810075177.6A CN108363844B (zh) | 2018-01-26 | 2018-01-26 | 一种航空发动机启动过程排气温度预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108363844A true CN108363844A (zh) | 2018-08-03 |
CN108363844B CN108363844B (zh) | 2020-04-24 |
Family
ID=63007146
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810075177.6A Active CN108363844B (zh) | 2018-01-26 | 2018-01-26 | 一种航空发动机启动过程排气温度预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108363844B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109272174A (zh) * | 2018-11-14 | 2019-01-25 | 上海交通大学 | 基于循环神经网络的燃气轮机排气系统工况预测方法 |
CN109598371A (zh) * | 2018-11-21 | 2019-04-09 | 北京航空航天大学 | 一种飞行器载电子设备的温度预测方法及系统 |
CN109858714A (zh) * | 2019-04-03 | 2019-06-07 | 杭州安脉盛智能技术有限公司 | 基于改进神经网络的烟丝质检指标预测方法、装置及系统 |
CN109886460A (zh) * | 2019-01-16 | 2019-06-14 | 中国计量大学 | 基于adaboost的隧道沉降时间序列的预测方法 |
CN111649951A (zh) * | 2020-04-15 | 2020-09-11 | 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) | 航空发动机故障检测方法、装置、计算机设备及存储介质 |
CN114165338A (zh) * | 2020-09-11 | 2022-03-11 | 中国航发商用航空发动机有限责任公司 | 发动机冷却气流温度预测方法和装置、控制方法和设备 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102207015A (zh) * | 2011-05-20 | 2011-10-05 | 潍柴动力股份有限公司 | 一种scr催化器的温度预测装置和方法 |
CN103195583A (zh) * | 2013-04-07 | 2013-07-10 | 北京华清燃气轮机与煤气化联合循环工程技术有限公司 | 一种采用排气温度分散度监测保护燃气轮机燃烧的方法 |
CN103955742A (zh) * | 2014-04-28 | 2014-07-30 | 淮阴工学院 | 一种基于集成学习的pm2.5预报方法 |
CN104951836A (zh) * | 2014-03-25 | 2015-09-30 | 上海市玻森数据科技有限公司 | 基于神经网络技术的发帖预测系统 |
CN105005583A (zh) * | 2015-06-17 | 2015-10-28 | 清华大学 | 一种社交网络中信息转发增量的预测方法及系统 |
CN105354363A (zh) * | 2015-10-09 | 2016-02-24 | 上海大学 | 基于核极限学习机的脉动风速预测方法 |
US20170061305A1 (en) * | 2015-08-28 | 2017-03-02 | Jiangnan University | Fuzzy curve analysis based soft sensor modeling method using time difference Gaussian process regression |
-
2018
- 2018-01-26 CN CN201810075177.6A patent/CN108363844B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102207015A (zh) * | 2011-05-20 | 2011-10-05 | 潍柴动力股份有限公司 | 一种scr催化器的温度预测装置和方法 |
CN103195583A (zh) * | 2013-04-07 | 2013-07-10 | 北京华清燃气轮机与煤气化联合循环工程技术有限公司 | 一种采用排气温度分散度监测保护燃气轮机燃烧的方法 |
CN104951836A (zh) * | 2014-03-25 | 2015-09-30 | 上海市玻森数据科技有限公司 | 基于神经网络技术的发帖预测系统 |
CN103955742A (zh) * | 2014-04-28 | 2014-07-30 | 淮阴工学院 | 一种基于集成学习的pm2.5预报方法 |
CN105005583A (zh) * | 2015-06-17 | 2015-10-28 | 清华大学 | 一种社交网络中信息转发增量的预测方法及系统 |
US20170061305A1 (en) * | 2015-08-28 | 2017-03-02 | Jiangnan University | Fuzzy curve analysis based soft sensor modeling method using time difference Gaussian process regression |
CN105354363A (zh) * | 2015-10-09 | 2016-02-24 | 上海大学 | 基于核极限学习机的脉动风速预测方法 |
Non-Patent Citations (1)
Title |
---|
丁刚等: "基于过程神经网络的航空发动机排气温度预测", 《航空动力学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109272174A (zh) * | 2018-11-14 | 2019-01-25 | 上海交通大学 | 基于循环神经网络的燃气轮机排气系统工况预测方法 |
CN109598371A (zh) * | 2018-11-21 | 2019-04-09 | 北京航空航天大学 | 一种飞行器载电子设备的温度预测方法及系统 |
CN109598371B (zh) * | 2018-11-21 | 2021-02-02 | 北京航空航天大学 | 一种飞行器载电子设备的温度预测方法及系统 |
CN109886460A (zh) * | 2019-01-16 | 2019-06-14 | 中国计量大学 | 基于adaboost的隧道沉降时间序列的预测方法 |
CN109858714A (zh) * | 2019-04-03 | 2019-06-07 | 杭州安脉盛智能技术有限公司 | 基于改进神经网络的烟丝质检指标预测方法、装置及系统 |
CN111649951A (zh) * | 2020-04-15 | 2020-09-11 | 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) | 航空发动机故障检测方法、装置、计算机设备及存储介质 |
CN114165338A (zh) * | 2020-09-11 | 2022-03-11 | 中国航发商用航空发动机有限责任公司 | 发动机冷却气流温度预测方法和装置、控制方法和设备 |
CN114165338B (zh) * | 2020-09-11 | 2023-06-23 | 中国航发商用航空发动机有限责任公司 | 发动机冷却气流温度预测方法和装置、控制方法和设备 |
Also Published As
Publication number | Publication date |
---|---|
CN108363844B (zh) | 2020-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108363844A (zh) | 一种航空发动机启动过程排气温度预测方法 | |
CN115358432B (zh) | 用于智慧燃气的燃气表计量故障确定方法和物联网系统 | |
US11863920B2 (en) | Methods and internet of things systems for gas usage safety warning based on smart gases | |
EP3550447B1 (en) | Performing cache update adaption | |
CN110969290A (zh) | 一种基于深度学习的径流概率预测方法及系统 | |
CN105550368A (zh) | 一种高维数据的近似最近邻检索方法及检索系统 | |
CN113516291A (zh) | 电力负荷预测方法、装置及设备 | |
Rusin et al. | Practical algorithms for online thermal stress calculations and heating process control | |
CN110147820A (zh) | 航班额外油量的推荐方法、装置、设备及存储介质 | |
CN110598336B (zh) | 热水器用水量预测方法及装置、热水器及电子设备 | |
CN115681821A (zh) | 用于智慧燃气设备管理的加臭自动控制方法和物联网系统 | |
CN110866682A (zh) | 基于历史数据的地下电缆预警方法和装置 | |
CN115329823A (zh) | 延迟焦化加热炉炉管温度的动态深度迁移预测方法及装置 | |
Yang et al. | RS-SVM forecasting model and power supply-demand forecast | |
Peng et al. | Short‐term wind power prediction based on stacked denoised auto‐encoder deep learning and multi‐level transfer learning | |
CN113837467B (zh) | 点对点能源项目评估方法、装置、计算机设备和存储介质 | |
Wang et al. | Trajectory similarity-based prediction with information fusion for remaining useful life | |
Sun et al. | Research on Satellite Fault Detection Method Based on MSET and SRPRT | |
CN111160614B (zh) | 资源转移预测模型的训练方法、装置及计算设备 | |
Hu et al. | Remaining useful life prediction for aircraft engines based on phase space reconstruction and hybrid VNS-SVR model | |
CN113297254A (zh) | 一种概念化查询方法和装置 | |
CN112614006A (zh) | 负荷预测方法、装置、计算机可读存储介质以及处理器 | |
CN108304449B (zh) | 基于自适应数据集划分方式的大数据Top-k查询方法 | |
Shah et al. | Comparative study on estimation of remaining useful life of turbofan engines using machine learning algorithms | |
Xu et al. | LS-SVM Based Two Step Commercial Short-Term Load Prediction |
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 |