首页 / 专利库 / 空气动力学 / 旋转失速 / 基于单分类超限学习机的轴流压气机气动失稳检测方法

基于单分类超限学习机的轴流压气机气动失稳检测方法

阅读:794发布:2020-05-11

专利汇可以提供基于单分类超限学习机的轴流压气机气动失稳检测方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种基于单分类超限学习机的轴流 压气机 气动 失稳检测方法。首先对采集到的动态压 力 信号 ,进行去趋势项与平滑处理和滑动窗口法做数据截取;其次选取多个时域统计特征来描述所截取的每段窗口数据,通过自 编码器 采用ELM 算法 作为 优化算法 ,对动态压力数据特征进行编码和解码处理,得到描述该工况下压力数据的二次重构特征;将只包含正常工况下的数据特征作为单分类超限学习机的输入,进行模型训练;最后将包含正常和气动失稳工况的测试数据按照时间顺序输入已训练的单分类模型,通过对数据特征的分类,实现对气动失稳的检测。本发明通过计算简单且高效的OC-ELM进行快速识别分类,准确率高,为轴流压气机的气动失稳检测提供了有效且便利的工具。,下面是基于单分类超限学习机的轴流压气机气动失稳检测方法专利的具体信息内容。

1.基于单分类超限学习机的轴流压气机气动失稳检测方法,其特征在于对采集的动态压信号预处理之后,利用滑动窗口方法实现数据截取,从每个窗口数据提取多个时域统计特征;将正常工况下压力数据的时域统计特征经过多层ELM-AE自编码器编码和解码处理,得到二次重构的特征;将重构的特征作为ELM最后一层的输入,获得ELM的实际输出;通过实际输出与正常类标签的误差,设置最优阈值参数来区分正常和异常数据;最后将包含正常和气动失稳工况下的测试样本数据作为该多层单分类超限学习机的输入,实现对轴流压气机气动失稳检测。
2.根据权利要求1所述的基于单分类超限学习机的轴流压气机气动失稳检测方法,其特征在于具体包含如下步骤:
步骤1、数据预处理:
步骤2、选取时域统计特征:对每个滑动窗口中的数据选取12个时域统计特征;
步骤3、输入训练样本归一化处理;
步骤4、样本数据二次特征提取;
步骤5、计算分类学习的输出权重β;
步骤6、计算单分类器阈值θ
步骤7、输入测试数据进行测试;
步骤8、旋转失速时间检测。
3.根据权利要求1或2所述的基于单分类超限学习机的轴流压气机气动失稳检测方法,其特征在于步骤1具体实现如下:
1-1、去趋势项与平滑处理:采集压气机机匣壁面周向动态压力数据,采用常规的滑动平均法和五点三次平滑法做去趋势项和平滑处理;动态压力数据包含正常和气动失稳两部分工况数据;
1-2、数据截取:采用滑动窗口方法对数据进行截取,选取窗口大小和滑动距离为转子周期T的倍数,1T为64个采样点,滑动距离为32个采样点。
4.根据权利要求3所述的基于单分类超限学习机的轴流压气机气动失稳检测方法,其特征在于步骤2所选取的时域统计特征为:平均值、均方值、均方根值、方均根值、方差、标准差、峰值、波形指标、峰值指标、脉冲指标、裕度指标、峭度指标。
5.根据权利要求4所述的基于单分类超限学习机的轴流压气机气动失稳检测方法,其特征在于步骤步骤3所述的输入训练样本归一化处理,具体实现如下:
3-1.将步骤2提取的时域统计特征作为输入训练样本:
其中 表示第i个样本的第α个时域统计特征,α=12表示时域统计特征个数, 作为样本标签,表示其为正常转速样本,训练集只包含正常转速样本,N为总的训练样本数量。
6.根据权利要求5所述的基于单分类超限学习机的轴流压气机气动失稳检测方法,其特征在于步骤4所述的样本数据二次特征提取,具体实现步骤如下:
4-1、归一化后的训练样本 即ELM-AE的输入和输出矩阵;
4-2、随机生成隐层输入权值矩阵 和正交化偏置向量矩阵
将输入数据映射到相同或者不同的数据维度空间:hk=g(akxα+bk),(ak)Tak=I,(bk)Tbk=1,其中:g()表示激活函数,k=1,2,...,K为ELM-AE个数;
4-3、求解ELM-AE的输出权值矩阵
设ELM-AE个数为K,输入输出层神经元数量为d,隐含层神经元数量为 以及每个隐含层的 正则化参数 若 或 即对于稀疏或压缩的特征表达,
若 即对于等维度的特征映射,βk=H-1XA,(βk)Tβk=I,k=
1,2,...,K;其中: 表示隐含层输出矩阵,即经过K个ELM-AE提取的特征矩阵。
7.根据权利要求6所述的基于单分类超限学习机的轴流压气机气动失稳检测方法,其特征在于步骤5所述的计算分类学习的输出权重β,具体实现如下:
5-1、经过多层ELM-AE系统得到输出X'=[x′1,x′2,...,x′N],输入层神经元数量为d,ELM分类层隐含神经元数量为L以及ELM分类层的正则化参数C;


隐藏层输出节点矩阵为 且
5-2、得到实际输出
8.根据权利要求7所述的基于单分类超限学习机的轴流压气机气动失稳检测方法,其特征在于步骤6所述的计算单分类器阈值θ,具体实现如下:
6-1、计算实际输出Y与样本标签T的误差,
6-2、将得到的误差 从大到小排序,得到 其中 和 分别
表示最大和最小的误差距离;
6-3、设置一个阈值参数μ,得到阈值为θ=εfloor(μ·N)。
9.根据权利要求8所述的基于单分类超限学习机的轴流压气机气动失稳检测方法,其特征在于步骤7所述的输入测试数据进行测试,具体实现如下:
7-1、给定一系列测试样本 其中 表示第i个样本,
表示其为正常转速样本, 表示其旋转失速样本,P为总的测试样本数量;
7-2、输入到多层ELM-AE得到每层输出 令最后一层ELM-AE输出为Y'=
[y′1,y′2,...,y′N];
7-3、输入到ELM分类层得到Yβ=β·Y',计算实际输出与样本标签T的误差
7-4、将得到的误差 与单分类器阈值θ比较:
10.根据权利要求9所述的基于单分类超限学习机的轴流压气机气动失稳检测方法,其特征在于步骤8所述的旋转失速时间检测具体实现如下:
将包含正常和气动失稳工况下的测试数据集输入已经训练完成的多层单分类超限学习机,设输出的第一个为-1的样本标签为 P为总测试样本数量;由步骤1的滑动窗口方法得到气动失稳的窗口节点为 即为轴流压气机气动失稳的检测时间。

说明书全文

基于单分类超限学习机的轴流压气机气动失稳检测方法

技术领域

[0001] 本发明涉及轴流压气机气动失稳检测方法,具体涉及基于多层单分类超限学习机的轴流压气机气动失稳检测方法。

背景技术

[0002] 轴流压气机作为航空发动机的核心部件之一,广泛应用于军用和大型民用飞机的动装置中。当前航空发动机的发展趋势是追求更高的单级压比和更少的级数。单级压比的提高和级数的减少,必然要求压气机叶片负荷越来越高,导致发动机稳定裕度的减少和抗干扰能力的下降,更容易引起流动分离,即气动失稳现象的发生。气动失稳的发生不但使发动机性能(推力、经济性)大为恶化,限制发动机的工作范围,而且更严重的是它们会引起发动机突然熄火,或引起压气机叶片剧烈振动以致叶片断裂而造成整台发动机的损坏。因此,轴流压气机气动失稳检测的研究对提高系统的安全性和可靠性具有重要的研究意义和工程应用价值。
[0003] 目前,关于轴流压气机气动失稳检测面临诸多挑战。一、轴流压气机从正常转速进入气动失稳状态的数据并没有明确的界限;二、轴流压气机气动失稳状态的异常数据难以获取,且数据量很少,导致正常与异常数据分布不平衡;三、现有的异常检测单分类器针对高维大数据样本时,往往存在表征能力不够,识别率低等问题;四、数据本身存在噪声,噪声对正常数据产生影响,而且和异常数据难以区分。
[0004] 已知现有的异常检测方法中,基于传统的决策树支持向量机、最近邻等方法进行异常检测不仅准确率不高,而且面对高维不均衡数据集时,对数据的表征能力不够,导致训练效果差且无法满足实时性的需求。因此,如何在更一般的数据集中通过一种优异的算法获得表征能力强泛化性能好识别率高训练速度快的异常检测模型是一个核心问题。本文采用了一种基于多层自编码器和超限学习机(Extreme Learning Machine,ELM)算法的多层单分类超限学习机来进行轴流压气机气动失稳检测。
[0005] ELM是一种简单易用、有效的单隐层前馈神经网络学习算法。它在异常检测分类上有以下几点优势:(1)ELM在输入层隐藏层之间采用随机权重。我们可以多次训练相同的数据集,这给不同的分类精度不同的输出空间。(2)ELM是一个更简单的前馈神经网络学习算法。传统的神经网络学习算法(如BP算法)需要人为设置大量的网络训练参数,这样就会非常容易产生局部最优解。而ELM在确定网络参数的过程中,只需要设置网络的隐层节点个数,在算法执行过程中不需要调整网络的输入权值以及隐元的偏置,并且产生唯一的最优解。因此,ELM学习速度比传统人工神经网络更快且泛化性能更好,可以最快地实现模型的训练和测试。
[0006] ELM的输出为 其中:βi是隐层节点和输出节点之间的权重,G(ai,bi,x)是隐层输出函数。h(x)=[G(a1,b1,x),...,G(aL,bL,x)]T是隐层相对于输入x的输出向量。ELM的关键在于最小化训练误差和输出权重范数。
[0007] 即最小化 以及||β||。
[0008] ELM算法总结如下:给定训练集{(xi,ti)|xi∈Rn,ti∈Rm,i=1,2,...N},隐层节点输出函数g(w,b,x)和隐层节点数L。
[0009] (1)随机分配隐层节点的参数(wi,bi),i=1,2...,L。
[0010] (2)计算隐层输出矩阵Η。
[0011] (3)计算隐层节点和输出节点之间的权重β=H+T。
[0012] H+是隐层输出矩阵H的Moore-Penrose广义逆矩阵,可以使用正交投影法、正交化法和奇异值分解等方法计算得到。

发明内容

[0013] 本发明的目的是针对现有检测算法中存在的问题,提供一种基于单分类超限学习机的轴流压气机气动失稳检测方法,是一种更快更高效的异常检测方法。
[0014] 本发明的具体技术方案通过如下步骤实现:
[0015] 步骤1、数据预处理
[0016] 1-1、去趋势项与平滑处理:采集压气机机匣壁面周向动态压力数据(包含正常和气动失稳两部分工况数据),采用常规的滑动平均法和五点三次平滑法做去趋势项和平滑处理。
[0017] 1-2、数据截取:为了便于快速处理和实时检测,采用滑动窗口方法对数据进行截取(选取窗口大小和滑动距离为转子周期T的倍数,1T为64个采样点,滑动距离为32个采样点)。
[0018] 步骤2、特征选取
[0019] 2-1、选取时域统计特征:对每个滑动窗口中的数据选取12个时域统计特征。时域统计特征分别为:平均值、均方值、均方根值、方均根值、方差、标准差、峰值、波形指标、峰值指标、脉冲指标、裕度指标、峭度指标。
[0020] 步骤3、输入训练样本归一化处理:
[0021] 3 - 1 、将 步 骤 2 提 取 的 时 域 统 计 特 征 作 为 输 入 训 练 样 本其中 表示第i个样本的第α个时域统计特征,α=12表示时域统计特征个数, 作为样本标签,表示其为正常转速样本,训练集只包含正常转速样本,N为总的训练样本数量。
[0022] 步骤4、样本数据二次特征提取
[0023] 4-1、归一化后的训练样本 即ELM-AE的输入和输出矩阵;
[0024] 4-2、随机生成隐层输入权值矩阵 和正交化偏置向量矩阵将输入数据映射到相同或者不同的数据维度空间:hk=g(akxα+bk),(ak)Tak=I,(bk)Tbk=1,其中:g()表示激活函数,k=1,2,...,K为ELM-AE个数;
[0025] 4-3、求解ELM-AE的输出权值矩阵
[0026] 设ELM-AE个数为K,输入输出层神经元数量为d,隐含层神经元数量为 以及每个隐含层的正则化参数 若 或 即对于稀疏及压缩的特征表达:
[0027]
[0028] 若 即对于等维度的特征映射,βk=H-1XA,(βk)Tβk=I,(k=1,2,...,K);则表示隐含层输出矩阵,即经过K个ELM-AE提取的特征矩阵。
[0029] 步骤5、分类学习的输出权重β计算
[0030] 5-1、经过多层ELM-AE网络得到输出X'=[x'1,x'2,...,x'N],输入层神经元数量为d,ELM分类层隐含神经元数量为L以及ELM分类层的正则化参数C;
[0031] 若
[0032] 若
[0033] 隐藏层输出节点矩阵为 且
[0034]
[0035] 5-2、得到实际输出
[0036] 步骤6、计算单分类器阈值θ
[0037] 6-1、计算实际输出Y与样本标签T的误差,
[0038]
[0039] 6-2、将得到的误差 进行从大到小排序,得到 其中和 分别表示最大和最小的误差距离;
[0040] 6-3、设置一个阈值参数μ,得到阈值θ=εfloor(μ·N);
[0041] 步骤7、输入测试数据进行测试
[0042] 7-1、给定一系列测试样本 其中 表示第i个样本, 表示其为正常转速样本, 表示其旋转失速样本,P为总共的测试样本数量;
[0043] 7-2、输入到多层ELM-AE得到每层输出 令最后一层ELM-AE输出为Y'=[y'1,y'2,...,y'N];
[0044] 7-3、输入到ELM分类层得到Yβ=β·Y',计算实际输出与样本标签T的误差距离[0045] 7-4、将得到的误差 与单分类器阈值θ比较;
[0046]
[0047] 步骤8、旋转失速时间检测
[0048] 将包含正常和气动失稳工况下的测试数据集作为训练完成的多层单分类超限学习机模型的输入,设测试输出的第一个为-1的样本标签为 P为总测试样本数量。由步骤1的滑动窗口方法得到旋转失速的窗口节点为 即为轴流压气机进入旋转失速的时间。
[0049] 本发明与现有技术相比,具有如下优点和有益效果:
[0050] 本发明采用的是超限学习机自编码器(ELM-AE)算法进行异常检测数据的特征提取,该算法是一种比普通自编码(AE)算法更为高效的自编码算法,它利用ELM替换原始自编码器进行特征提取的BP梯度下降法,可以快速处理较高维度的输入数据,提取其主要信息,实现对原始数据的特征提取。通过多个ELM-AE神经网络叠加增强了特征提取能力,特别是针对异常检测中的高维不均衡大数据集。最后将提取得到的特征通过ELM单分类器进行检测识别,分类速度快泛化性能更好,满足实时性和处理更一般性数据要求。
[0051] 本发明相对于基于传统决策树、支持向量机和最近邻的异常检测算法,可以更为快速高效地提取主要信息并降维,再进行识别分类。准确率更高,泛化性能更好。不仅适用于小数据集也同样适用于高维大数据集,具有普适性。对于今后的实际应用具有重要的意义。
[0052] 本发明测试了轴流压气机实测信号的异常检测效果,具有复杂情况下较强的抗干扰能力和实时性,基于ELM-AE和ELM单分类器相结合的多层单分类超限学习机算法能准确识别轴流压气机从正常工况进入气动失稳的样本数据,确定轴流压气机即将进入气动失稳的时间。附图说明
[0053] 图1为本发明实施例中轴流压气机气动失稳检测方法的流程图
[0054] 图2为本发明实施例中某一轴流压气机从正常到气动失稳工况下的动态压力数据图;
[0055] 图3为单隐层前馈神经网络示意图;
[0056] 图4为ELM-AE网络结构图;
[0057] 图5为多层单分类超限学习机网络结构图;
[0058] 图6为多层单分类超限学习机流程示意图;
[0059] 图7为检测到轴流压气机进入气动失稳的时间点

具体实施方式

[0060] 为了使本发明所要解决的技术问题、技术方案及有益效果更加清楚、明白,以下结合附图和实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
[0061] 本实施例是一种基于单分类超限学习机的轴流压气机气动失稳检测方法。通过采集轴流压气机从正常转速到气动失稳过程的多组动态压力数据,经过预处理之后提取每个窗口数据的时域特征指标,将只包含正常转速的样本数据作为训练集进行模型训练,然后将包含正常转速和旋转失速的样本数据作为测试集,可以快速、准确地检测气动失稳发生的时间。
[0062] 以下结合图1对本发明的实施方法作进一步说明。
[0063] 步骤1、数据预处理
[0064] 1-1、首先在压气机机匣内壁安装动态压力传感器,在2500r/min、2700r/min、2900r/min的转速下获取失稳前和失稳过程的动态压力信号,采样频率为1000HZ,图2显示了采集的原始数据信号。在数据采集过程中,一方面由于传感器频率范围外低频性能的不稳定性以及传感器周围的环境干扰,往往会偏离基线,偏离基线随时间变化的整个过程被称为信号的趋势项,直接影响信号的正确性,常用滑动平均法消除趋势项。另一方面由于不规则的随机干扰信号的频带较宽,有时高频成分所占比例还很大,为了削弱干扰信号的影响,常用五点三次滑动平均法对采样数据做平滑处理。
[0065] 1-2、经数据处理后的动态压力信号利用滑动窗口实现数据截取,每个窗口长度为一个转子周期包含64个采样点,窗口滑动距离为32个采样点。
[0066] 步骤2、特征提取
[0067] 2-1、对步骤2截取的每个窗长数据提取12个时域特征指标,每隔一个旋转周期提取时域统计特征作为输入训练样本。
[0068] 步骤3、输入训练样本归一化处理:
[0069] 3 - 1 、对 于 步 骤 2 提 取 的 时 域 统 计 特 征 作 为 输 入 训 练 样 本做标准归一化处理,其中 表示第i个样本,α=12表示时域统计特征个数, 作为样本标签表示其为正常转速样本,训练集只包含正常转速样本,N为总的训练样本数量。
[0070] 步骤4、样本数据二次特征提取
[0071] 4-1、归一化后的训练样本 即ELM-AE的输入和输出矩阵;图3显示了单隐层前馈神经网络的基本结构,也就是ELM算法的框架。图4显示了单个ELM-AE的网络结构,也就是自编码器采用了ELM算法作为优化算法,加快了自编码器的训练速度,增强了泛化能力。
[0072] 4-2、随机生成隐层输入权值矩阵 和正交化偏置向量矩阵将输入数据映射到相同或者不同的数据维度空间:hk=g(akxα+bk),(ak)Tak=I,(bk)Tbk=1,其中:g()表示激活函数,(k=1,2,...,K)为ELM-AE个数;
[0073] 4-3、求解ELM-AE的输出权值矩阵
[0074] 设ELM-AE个数为K=3,输入输出层神经元数量为d,隐含层神经元数量为 以及每个隐含层的正则化参数 取 即对于压缩的特征表达,βk=H-1XA,(βk)Tβk=I,(k=1,2,...,K),其中: 表示隐含层输出矩阵,即经过K个ELM-AE提取的特征矩阵。
[0075] 步骤5、分类学习的输出权重β计算
[0076] 5-1、图5即为本发明采用方法的网络结构,三层ELM-AE堆叠后再通过一个无隐含层的ELM分类层。
[0077] 经过多层ELM-AE网络得到输出X'=[x'1,x'2,...,x'N],输入层神经元数量为d,ELM 分类 层隐 含 神经元 数量为 L以及 ELM 分类 层的正 则化参数 C ,取隐藏层输出节点矩阵为且
[0078]
[0079] 5-2、得到实际输出
[0080] 步骤6、计算单分类器阈值θ
[0081] 6-1、计算实际输出Y与样本标签T的误差,
[0082]
[0083] 6-2、将得到的误差 进行从大到小排序,得到 其中和 分别表示最大和最小的误差距离;
[0084] 6-3、设置一个阈值参数μ,得到阈值为θ=εfloor(μ·N);
[0085] 步骤7、输入测试数据进行测试
[0086] 7-1、给定一系列测试样本 其中 表示第i个样本, 表示其为正常转速样本, 表示其气动失稳样本,P为总共的测试样本数量;
[0087] 7-2、输入到多层ELM-AE得到每层输出 令最后一层ELM-AE输出为Y'=[y'1,y'2,...,y'N];
[0088] 7-3、输入到ELM分类层得到Yβ=β·Y',计算实际输出与样本标签T的误差距离[0089] 7-4、将得到的误差距离 与单分类器阈值θ比较;
[0090]
[0091] 步骤8、检测时间
[0092] 如图6所示,首先取训练数据(仅正常时域特征数据集)输入到多层ELM-AE进行二次特征提取,再通过ELM分类层(无隐含层)进行分类输出实际结果,根据得到的实际输出与样本标签的误差排序,通过阈值参数确定阈值。之后将待测试的样本输入训练好的异常检测模型,得到测试数据实际输出和样本标签的误差,大于阈值的归为异常类,小于等于阈值的为正常类。图7显示了由多层随机神经网络检测到的第一个为-1的异常类样本标签为P为总测试样本数量,由滑动窗口方法得到气动失稳的窗口节点即为轴流压气机即将进入气动失稳的时间。
[0093] 通过实施例的实验结果可以看出,本发明能直观、有效地判断出即将进入气动失稳的时间。本发明的识别结果可以作为轴流压气机发生气动失稳检测,为提前采取主动控制,避免事故的恶化提供重要参考。
[0094] 上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈