基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部件自适应性能诊断方法
【专利摘要】本发明的目的在于提供基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部件自适应性能诊断方法，利用部件相对折合参数，建立燃气轮机部件级热力模型，用相对折合参数重新定义压气机和透平的气路健康指数，通过设置热力模型中的各个部件的气路健康指数，来模拟单部件和多部件性能衰退时的气路测量参数，通过灰色关联理论算法实现性能衰退模式识别，从而隔离实际真正发生性能衰退的部件，从而确定了待诊断的实际气路部件健康指数的子集。本发明能更简单方便地用于气路诊断时设置各个气路部件健康指数，更准确地表征由于部件性能衰退而导致的气路健康指数的变化，能更好解决由于测量噪音和部件数目较多导致诊断结果出现模糊效应的问题。
【专利说明】
基于灰色关联理论与热力模型相结合的混合型燃气轮机气路 部件自适应性能诊断方法
技术领域
[0001] 本发明设及的是一种诊断方法，具体地说是燃气轮机的诊断方法。
【背景技术】
[0002] 在20世纪下半叶，随着燃气轮机在航空工业中广泛应用，越来越多的受到工业电 站领域、石油和天然气管道运输W及舰船工业领域的关注。在运行中，由于高溫、高压、高转 速及高应力的恶劣工作条件及环境污染等影响，各种类型的燃气轮机都会逐渐性能衰退。 燃气轮机的主要气路部件包含压气机、燃烧室和透平。运些主要部件会随着时间遭受不同 的退化现象，如污垢、泄漏、腐蚀、热崎变、外来物损坏等，将会引起性能恶化并易导致各种 严重的故障发生运行安全问题。因此，对于燃气轮机用户来说，当前燃气轮机的性能健康状 况是非常重要的信息。
[0003] 目前大多数燃气轮机的维修策略是预防性维修保养，通常按照燃机制造商指示的 当量运行小时数化OH)来考虑是否需要小修、中修、大修。燃机停运，无论是计划内或计划外 的，总是意味着昂贵的成本代价。为了节省维修费用，用户需要根据燃机实际的性能健康状 况采取维修策略，即预测性维修保养。
[0004] 基于热力模型决策的气路部件性能诊断方法已经广泛应用于燃气轮机性能健康 状态监测中，并且已经成为支持维修策略改革的关键技术之一。通常的基于热力模型决策 的气路部件性能诊断方法使用气路部件的性能参数(绝对参数)来定义部件健康指数，因此 在诊断前气路测量参数需要进行数据预处理来消除环境条件或操作条件及工质组分变化 而导致燃气轮机运行性能变化的影响。并且由于W部件性能参数为自适应变量，所W诊断 时通常需要分两步:第一步为根据气路测量参数通过线性或非线性牛顿-莱普生迭代算法 计算得到部件性能参数，如空气流量、压气机压比、压气机等赌效率、透平前溫、透平等赌效 率等绝对性能参数;第二步为在同一部件特性图上比较实际性能衰退情况下的部件运行点 与健康基准情况下的部件运行点，从而观测此时部件特性图上的特性线发生偏移的程度 (即得到气路部件健康指数），从而来评估当前部件的性能健康状况。当燃气轮机中参与诊 断的部件数目增多时，故障系数矩阵的维数会随之增大，加之受到测量噪音的干扰，模糊效 应(即，尽管某些部件实际上并无发生性能退化，但诊断出的性能衰退情况几乎分布在所有 的部件气路健康指数上)可能会很强，导致实际性能衰退的部件未被识别出来。
【发明内容】
[0005] 本发明的目的在于提供可W更好解决由于测量噪音和部件数目较多导致诊断结 果出现模糊效应问题的基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部件 自适应性能诊断方法。
[0006] 本发明的目的是运样实现的：
[0007] 本发明基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部件自适应 性能诊断方法，其特征是：
[0008] (1)基于对象燃气轮机健康时的气路测量参数，建立燃气轮机非线性热力模型，压 气机和透平都用相似折合参数形式表示；
[0009] (2)用相似折合参数重新定义压气机和透平的气路健康指数，消除由于环境条件 变化而导致燃气轮机运行性能变化的影响；
[0010] (3)通过设置热力模型中各个部件的气路健康指数，来模拟单部件和多部件性能 衰退时的气路测量参数，用于积累性能衰退模式和衰退征兆的知识数据库，作为灰色关联 理论实现对实测气路参数进行模式识别的依据，所述热力模型的部件包括压气机、透平和 燃烧室；
[0011] (4)采集当前燃气轮机稳定运行时的一个时段的气路测量参数，进行降噪处理后 作为待离线诊断的气路测量参数；
[0012] (5)设置已建立的燃气轮机热力模型的环境输入条件和操作输入条件与采样时的 对象燃气轮机运行工况一致，消除由于环境条件和操作条件变化而导致燃气轮机运行性能 变化的影响；
[0013] (6)将待离线诊断的气路测量参数与热力模型的气路参数数据之间的相对偏差输 入灰色关联理论算法中实现性能衰退模式识别，隔离已发生性能衰退的部件，从而实现故 障系数矩阵的降维；
[0014] (7) W待离线诊断的气路测量数据与热力模型的气路参数数据之间的偏差作为残 差，通过迭代算法计算得到当前压气机、透平和燃烧室的气路健康指数，用W评估对象燃气 轮机实际的性能健康状况。
[001引本发明还可W包括；
[0016] 1、步骤(1)中建立燃气轮机非线性热力模型的具体步骤如下：
[0017] (a)压气机特性线整理成通用的相似折合参数形式如下：
[001 引 Gc,cor,rel = f (ncor'rel,化,巧1)
[0019] nc,rel = f (ncor,rel,3IC,rel)
[0020] 其中
为相对折合转速，n为实际转速，攻为压气机进口 滞止溫度，Rg为流经压气机工质的气体常数，下角标0表示设计点；
[0021]
为相对折合流量，Ge为实际压气机进口流量， 巧:为压气机进口滞止压力，
％相对压比，町为实际压气机压比，nc,rel =化Alco为 相对等赌效率，化为实际压气机等赌效率；
[0022] 透平特性线整理成通用的相似折合参数形式如下：
[0023 ] Gt, cor ' rel = f ( ncor ' rel,化 rel)
[0024] nx, rel = f (ricor, rel, TlJ, rel )
[00巧]式中
为相对折合转速，n为实际转速，瑞为透平进口滞 止溫度，Rg为流经透平工质的气体常数，
[0026]
为相对折合流量，Gt为实际透平进口流量，
村目对压比，化rei = TItAit日为相对等赌效率，巾为实际透平等赌效率，下角标0表 示设计点；
[0027] (b)根据采集的燃气轮机健康时的气路测量参数，修正各个部件的特性线数据，使 所建的燃气轮机热力模型的计算值与气路实测参数相匹配，从而消除热力模型计算误差给 诊断结果带来的负影响。
[0028] 2、步骤(2)用相似折合参数重新定义压气机和透平的气路健康指数的具体步骤如 下：
[0029] 压气机气路健康指数定义如下：
[0030] SFc, FC = Gc, cor, rel, deg/Gc, cor, rel
[0031 ] A SFc.FC= (Gc,cor,rel,deg-Gc,cor,rel)/Gc,cor,rel
[0032] SFc,Eff = nC,degAlC
[0033] A SFc'Eff = (nc,deg-nc)Alc
[0034] 其中SFC,FC为压气机流量特性指数;Ge, cor,rel,deg为压气机性能衰退时的相对折合流 量;Gg, Gor, rel为压气机健康时的相对折合流量;SFg, Eff为压气机效率特性指数;化,deg为压气 机性能衰退时的等赌效率;化为压气机健康时的等赌效率；
[0035] 透平气路健康指数定义如下：
[0036] SFt,FC = Gt,cor,deg/GT,cor
[0037] A SFt.FC= (GT,cor,deg-GT,cor)/Gx,cor
[003引 SFt, Eff = WidegAlT
[0039] A SFx.Eff = (nX,deg-riT)/nT
[0040] 其中SFt,FC为透平流量性能指数;Gt,cor,deg为透平性能衰退时的折合流量;Gt,cor为 透平健康时的折合流量;SFt,Eff为透平效率性能指数;riT,deg为透平性能衰退时的等赌效率； %为透平健康时的等赌效率。
[0041] 3、步骤(3)积累性能衰退模式和衰退征兆的知识数据库作为灰色关联理论实现对 实测气路参数进行模式识别的依据的具体步骤如下：
[0042] (a)通过设置热力模型中各个部件的气路健康指数，来模拟单部件和多部件性能 衰退时的气路测量参数，并整理为性能衰退时的气路测量参数相对于健康时的相对偏差形 式：
[0043]
[0044] 其中，三为通过设置热力模型中的各个部件的气路健康指数来模拟单部件和多部 件性能衰退时的气路参数计算值，％为在相同环境条件和操作条件下部件性能健康时的热 力模型气路参数计算值，为整理为性能衰退时的气路测量参数相对于健康时的相对偏 差形式；
[0045] (b)积累性能衰退模式和衰退征兆的知识数据库作为灰色关联理论实现对实测气 路参数进行模式识别的依据，其中灰色关联理论中的知识数据库形式为：
[0046]
[0047] 其中(：^^' = 1，2^-111)为已知的部件性能衰退模式;(3^^' = 1，2，-，）为每一个气路可 测参数;m为总的已知的部件性能衰退模式数目；n为总的对象燃气轮机气路可测参数数目。 [004引4、步骤(6)将待离线诊断的气路测量参数与热力模型的气路参数数据之间的相对 偏差输入灰色关联理论算法中实现性能衰退模式识别，隔离已发生性能衰退的部件的具体 步骤如下：
[0049] (a)将待离线诊断的气路测量参数与热力模型的气路参数数据之间的相对偏差输 入灰色关联理论算法中实现性能衰退模式识别：
[0化0]对于pG(〇，l):
[0化1 ]
[0化2]
[0053] 其中P为分辨系数，Xi为一组待离线诊断的实测气路数据，CUi化）向化)）为Xi和Cj 的第k个特征参数间的灰色关联系数，C(Xi，Cj)为Xi和C北勺灰色关联度。
[0054] Xi被分类至当C(Xi，Cj)(j = l，2,…，m)取得最大值时所属的部件性能衰退模式类；
[0055] (b)隔离已发生性能衰退的部件，从而实现故障系数矩阵的降维，即确定待诊断的 实际气路部件健康指数A京的子集。
[0056] 5、步骤(7) W待离线诊断的气路测量参数与热力模型的气路参数数据之间的偏差 作为残差，通过迭代算法计算得到当前的各个部件的气路健康指数，用W评估对象燃气轮 机实际的性能健康状况的具体步骤如下：
[0化7] (a) W待离线诊断的气路测量参数J与热力模型计算的气路参数数据J之间的偏 差作为残差累= 5-;，通过牛顿-拉普森迭代算法进行迭代，当预测的气路测量参数与实际 的测量参数之差小于阔值，即p|<S，e是设定的阔值，则得到最终的部件健康指数(A涼U ;
[005引（b)根据预测的已隔离部件的性能衰退程度，评估当前对象燃气轮机总体性能健 康状况。
[0059] 本发明的优势在于：
[0060] (1)本发明根据燃气轮机气动热力学特性，利用部件相对折合参数，建立燃气轮机 部件级热力模型，其中压气机和透平特性线数据整理成通用的相对折合参数形式，相比于 现有技术，能更简单方便地用于气路诊断时设置各个气路部件健康指数。
[0061] (2)本发明用相对折合参数重新定义压气机和透平的气路健康指数(代表了部件 特性线的偏移），相比于现有技术，能更准确地表征由于部件性能衰退而导致的气路健康指 数的变化，消除由于环境条件(大气压力、溫度和相对湿度)变化而给诊断结果带来的负影 响。
[0062] (3)本发明通过设置热力模型中的各个部件(压气机、透平和燃烧室）的气路健康 指数，来模拟单部件和多部件性能衰退时的气路测量参数，用于积累性能衰退模式和衰退 征兆(整理为性能衰退时的气路测量参数相对于健康时的相对偏差形式)的知识数据库，作 为灰色关联理论实现对实测气路参数进行模式识别的依据，相比于现有技术，既能节省高 昂的试验成本，还能准确详尽地积累各种性能衰退模式和衰退征兆的知识数据库。
[0063] (4)本发明通过灰色关联理论算法实现性能衰退模式识别，从而隔离实际真正发 生性能衰退的部件，从而确定了待诊断的实际气路部件健康指数的子集，从故障系数矩阵 降维的角度来改善诊断结果的准确性，相比于现有技术，能更好解决由于测量噪音和部件 数目较多导致诊断结果出现模糊效应的问题。
【附图说明】
[0064] 图1为本发明的诊断过程流程图；
[0065] 图2为典型非线性气路诊断算法预测隔离部件性能衰退程度的迭代计算过程；
[0066] 图3为某型=轴船用燃气轮机气路工作截面标识图；
[0067] 图4为通过灰色关联算法(GRA)识别、隔离发生性能衰退部件的结构图；
[0068] 图5为通过典型非线性气路诊断算法预测隔离部件的性能衰退程度的结构图；
[0069] 图6某型=轴船用燃气轮机诊断案例的诊断结果。
【具体实施方式】
[0070] 下面结合附图举例对本发明做更详细地描述：
[0071] 结合图1-6,本发明基于灰色关联理论与热力模型相结合的混合型燃气轮机气路 部件自适应性能诊断方法，包括W下步骤：
[0072] 步骤1)，基于对象燃气轮机新投运(或健康)时的气路测量参数建立能完全反映各 个部件特性的燃气轮机非线性热力模型，其中压气机和透平都用相似折合参数形式表示；
[0073] 步骤2)，用相似折合参数重新定义压气机和透平的气路健康指数，消除由于环境 条件(大气压力、溫度和相对湿度)变化而导致燃气轮机运行性能变化的影响；
[0074] 步骤3)，通过设置热力模型中的各个部件(压气机、透平和燃烧室）的气路健康指 数，来模拟单部件和多部件性能衰退时的气路测量参数，用于积累性能衰退模式和衰退征 兆(整理为性能衰退时的气路测量参数相对于健康时的相对偏差形式)的知识数据库，作为 灰色关联理论实现对实测气路参数进行模式识别的依据；
[0075] 步骤4)，采集当前对象燃气轮机稳定运行时的某一时段的气路测量参数，进行降 噪处理后作为待离线诊断的气路测量参数；
[0076] 步骤5)，设置已建立的燃气轮机热力模型的环境输入条件(大气压力、溫度和相对 湿度)和操作输入条件与采样时的对象燃气轮机运行工况一致，消除由于环境条件(大气压 力、溫度和相对湿度)和操作条件变化而导致燃气轮机运行性能变化的影响。
[0077] 步骤6)，将待离线诊断的实测气路数据与热力模型计算的气路参数数据之间的相 对偏差输入灰色关联理论算法中实现性能衰退模式识别，隔离已发生性能衰退的部件，从 而实现故障系数矩阵的降维（即确定了待诊断的实际气路部件健康指数的子集）；
[0078] 步骤7)，W待离线诊断的气路测量数据与热力模型计算的气路参数数据之间的偏 差作为残差，通过牛顿-拉普森迭代算法计算得到当前的各个部件(压气机、透平和燃烧室） 的气路健康指数，用W评估对象燃气轮机实际的性能健康状况。
[0079] 作为本发明基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部件自 适应性能诊断方法进一步的优化方案，步骤1)中所述基于对象燃气轮机新投运(或健康)时 的气路测量参数建立能完全反映各个部件特性的燃气轮机非线性热力模型的具体步骤如 下：
[0080] 步骤1.1)，根据燃气轮机气动热力学特性，利用部件相对折合参数，建立燃气轮机 部件级热力模型，其中压气机和透平特性线数据整理成通用的相对折合参数形式。
[0081 ]其中压气机特性线整理成通用的相似折合参数形式如下：
[0082] Gc,cor,rel = f (ncor,rel,3IC,rel)
[008；3] ric,rel = f (ncor,rel,化,rel)
[0084] 其4
为相对折合转速，n为实际转速，瑞为压气机进口 滞止溫度，Rg为流经压气机工质的气体常数，下角标0表示设计点；
[0085]
为相对折合流量，Ge为实际压气机进口流量， 每为压气机进口滞止压力，
为相对压比，化为实际压气机压比，化rel = %AlC0为 相对等赌效率，化为实际压气机等赌效率。
[0086] 透平特性线整理成通用的相似折合参数形式如下：
[0087] Gt, cor, rel 二 f (Hcor'rel,巧,rel)
[008引 TlT, rel 二 f (Ilcor, rel,叫,rel )
[0089] 式中：
为相对折合转速，n为实际转速，每为透平进口滞 止溫度，Rg为流经透平工质的气体常数
对相对折合流 量，Gt为实际透平进口流量，
为相对压比，riT,rei = TItAt日为相对等赌效率，TlT为实 际透平等赌效率，下角标0表示设计点。
[0090] 步骤1.2)，根据采集的对象燃气轮机新投运(或健康)时的气路测量参数（降噪处 理后），如总溫、总压、转速等，逐步修正各个部件的特性线数据(包括设计工况和变工况）， 使所建的燃气轮机热力模型的计算值与气路实测参数相匹配，从而消除热力模型计算误差 给诊断结果带来的负影响。
[0091] 作为本发明基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部件自 适应性能诊断方法进一步的优化方案，步骤2)中所述用相似折合参数重新定义压气机和透 平的气路健康指数的具体步骤如下：
[0092] 步骤2.1 )，燃气轮机总体性能健康状况可W由各主要部件的气路健康指数，如压 气机和透平的流量特性指数和效率特性指数、燃烧室的效率特性指数来表示。运里用相对 折合参数重新定义压气机和透平的气路健康指数(代表了部件特性线的偏移），W消除由于 环境条件(大气压力、溫度和相对湿度)变化而给诊断结果带来的负影响。
[0093] 步骤2.2)，压气机气路健康指数定义如下：
[0094] SFc, FC = Gc, cor,rel, deg/Gc, cor, rel
[0095] A SFc.FC= (Gc,cor,rel,deg-Gc,cor,rel)/Gc,cor,rel
[0096] SFc.Eff = nc,deg/nc
[0097] A SFc.Eff = (nc,deg-nc)/nc
[0098] 其中SFC, FC为压气机流量特性指数;Gc, cor, rel, deg为压气机(性能衰退时)相对折合流 量;Ge, cor, rel为压气机(健康时）相对折合流量；SFc, Ef f为压气机效率特性指数;nc, deg为压气 机(性能衰退时)等赌效率;化为压气机(健康时)等赌效率。
[0099] 燃烧室气路健康指数定义如下：
[0100] SFb,Eff = TiB,deg/riB
[0101] A SFb,Eff = (nB,deg-nB)/nB
[0102] 其中SFb,Eff为燃烧室燃烧效率性能指数;riB,deg为燃烧室(性能衰退时)燃烧效率;％ 为燃烧室(健康时)燃烧效率。
[0103] 透平气路健康指数定义如下：
[01 04] SFt, FC二Gt, cor, deg/Gl", cor
[010引 A SFt'FC= (GT,cor,deg-GT,cor)/GT,cor
[0106] SFt, Eff = TlMegAlT
[0107] A SFt,Eff = (TlMeg-TlT)AlT
[010引其中SFt,FC为透平流量性能指数;GT,cnr,deg为透平(性能衰退时）折合流量;Gt,cdt为 透平(健康时）折合流量;SFt, Ef f为透平效率性能指数;化deg为透平(性能衰退时）等赌效率； %为透平(健康时)等赌效率。
[0109] 作为本发明基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部件自 适应性能诊断方法进一步的优化方案，步骤3)中所述积累性能衰退模式和衰退征兆(整理 为性能衰退时的气路测量参数相对于健康时的相对偏差形式)的知识数据库作为灰色关联 理论实现对实测气路参数进行模式识别的依据的具体步骤如下：
[0110] 步骤3.1)，通过设置热力模型中的各个部件(压气机、透平和燃烧室）的气路健康 指数，来模拟单部件和多部件性能衰退时的气路测量参数，并整理为性能衰退时的气路测 量参数相对于健康时的相对偏差形式：
[0111]
[0112]其中，胃吉为通过设置热力模型中的各个部件(压气机、透平和燃烧室）的气路健康指 数来模拟单部件和多部件性能衰退时的气路参数计算值，云为在相同环境条件和操作条件 下部件性能健康时的热力模型气路参数计算值，aJ为整理为性能衰退时的气路测量参数 相对于健康时的相对偏差形式。
[0113] 步骤3.2)，积累性能衰退模式和衰退征兆(整理为性能衰退时的气路测量参数相 对于健康时的相对偏差形式)的知识数据库作为灰色关联理论实现对实测气路参数进行模 式识别的依据，其中灰色关联理论中的知识数据库形式为：
[0114]
[0115] 其中(：^^' = 1，2^-111)为已知的部件性能衰退模式;(3^^' = 1，2，-，）为每一个气路可 测参数;m为总的已知的部件性能衰退模式数目；n为总的对象燃气轮机气路可测参数数目。
[0116] 作为本发明基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部件自 适应性能诊断方法进一步的优化方案，步骤6)所述中将待离线诊断的实测气路数据与热力 模型计算的气路参数数据之间的相对偏差输入灰色关联理论算法中实现性能衰退模式识 另IJ，隔离已发生性能衰退的部件，从而实现故障系数矩阵的降维（即确定了待诊断的实际气 路部件健康指数的子集)的具体步骤如下：
[0117] 步骤6.1)，将待离线诊断的实测气路数据与热力模型(部件性能健康时)计算的气 路参数数据之间的相对偏差输入灰色关联理论算法中实现性能衰退模式识别：
[011 引对于 pG(0，l):
[0119]
[0120]
[0121] 其中P为分辨系数（其值通常设置为0.5)，Xi为某一组待离线诊断的实测气路数 据，Uxi化），C^k))为Xi和C撕第k个特征参数间的灰色关联系数，UXi尼)为Xi和C撕灰色 关联度。
[0122] Xi被分类至当|化尼Kj = I，2,…，m)取得最大值时所属的部件性能衰退模式类。
[0123] 步骤6.2)，隔离已发生性能衰退的部件，从而实现故障系数矩阵的降维（即确定了 待诊断的实际气路部件健康指数A恋的子集京)、,,,，）。
[0124] 作为本发明基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部件自 适应性能诊断方法进一步的优化方案，步骤7)所述中W待离线诊断的气路测量数据与热力 模型计算的气路参数数据之间的偏差作为残差，通过牛顿-拉普森迭代算法计算得到当前 的各个部件(压气机、透平和燃烧室）的气路健康指数，用W评估对象燃气轮机实际的性能 健康状况。具体步骤如下：
[01巧]步骤7.1)，W待离线诊断的气路测量数据热力模型计算的气路参数数据；之间 的偏差作为残差忘= 5-;，通过牛顿-拉普森迭代算法(New-Ra地son)计算得到当前的已隔 离部件(压气机、透平和燃烧室）的气路健康指数(A亦X,，。运里由于燃气轮机热力性能的强 非线性特性，一个迭代计算过程(牛顿-拉普森算法)用于计算部件健康指数(A承U。随着迭 代，当预测的气路测量参数非常接近实际的测量参数时（即，满足p||<f，e是设定的一个相 对较小的阔值），则可得到最终的部件健康指数，迭代过程如图2所示。
[0126] 步骤7.2)，根据预测的已隔离部件的性能衰退程度，评估当前对象燃气轮机总体 性能健康状况。
[0127] 本发明【具体实施方式】W某型=轴船用燃气轮机的气路部件健康诊断为例，其该型 =轴船用燃气轮机气路工作截面标识图如图3所示。该型=轴燃气轮机包括两个压气机(即 一个低压压气机化C)和一个高压压气机化C))、一个燃烧室(B)和S个透平（即一个高压透 平化T)、一个低压透平化T)和一个动力透平(PT))，其中发电机通过一个减速齿轮箱与动力 透平(PT)相连接。低压透平化T)的输出功通过低压轴驱动低压压气机化C)来压缩从进气道 出来的空气，高压透平化T)的输出功通过高压轴驱动高压压气机化C)来继续压缩从低压压 气机化C)出来的空气。从高压压气机化C)出来的高压空气进入燃烧室(B)与燃料发生燃烧 化学反应生成高溫、高压的燃气，燃气依次进入高压透平化T)、低压透平化T)和动力透平 (PT)来驱动透平输出功。最终，动力透平(PT)通过减速齿轮箱驱动发电机来产生电功率。同 时，从压气机中抽取的冷却空气流入热端气流通道去冷却各个透平前几级的静叶、动叶和 轮盘。当燃气轮机稳定运行时，发电机的电功率和动力透平(PT)的转速通常作为主要控制 参数而维持定常。该机组的气路测量参数如表1所示，各个气路部件的健康指数如表2所示。
[0128] 表1该型燃气轮机机组的气路测量参数
[0129]
[0130] 表2主要气路部件的健康指数
[0131]
[0132] 基于该型燃气轮机新投运(或健康）时的气路测量参数建立能完全反映各个部件 特性的燃气轮机非线性热力模型，其中压气机和透平都用相对折合参数形式表示。
[0133] 其中压气机特性线整理成通用的相似折合参数形式如下：
[0134] Gc,cor,rel = f (ncor,rel,3IC,rel)
[01 35] 化，rel 二 f (Ilcor, rel,化，rel )
[0136] 其4
'为相对折合转速，n为实际转速，C为压气机进口 滞止溫度，Rg为流经压气机工质的气体常数，下角标0表示设计点；
[0137]
.为相对折合流量，Ge为实际压气机进口流量， 感为压气机进口滞止压力
％相对压比，化为实际压气机压比，％,rel = %AlC0为 相对等赌效率，化为实际压气机等赌效率。
[0138] 透平特性线整理成通用的相似折合参数形式如下：
[0139] Gt, cor, rel = f (Hcor,rel, JiT, rel)
[0140] nx, rel = f (ricor, rel, TlJ, rel )
[0141] 式中
为相对折合转速，n为实际转速，C为透平进口滞 止溫度，Rg为流经透平工质的气体常数，
^相对折合流 量，Gt为实际透平进口流量
为相对压比，化Tel = IVrAlT日为相对等赌效率，TlT为实 际透平等赌效率，下角标0表示设计点。
[0142] 对象燃气轮机的热力模型在Matlab仿真平台上建立。热力模型的输入条件为环境 条件(大气溫度to、压力Po、相对湿度RH)、发电机输出功率化(作为操作条件）、燃料组分、燃 料低位热值、气路部件健康指数A恋(对于新投运机组，A：豪二0)。热力模型的计算输出为燃 料流量Gf、各个部件进出口气路截面处的热力参数(如总压、总溫)及转速等。
[0143] 用相对折合参数重新定义压气机和透平的气路健康指数(代表了部件特性线的偏 移），消除由于环境条件(大气压力、溫度和相对湿度)变化而导致燃气轮机运行性能变化的 影响。
[0144]压气机气路健康指数定义如下：
[01 45] SFc, FC二Ge, cor,rel, deg/Gc, cor, rel
[0146] A SFc.FC= (Gc,cor,rel,deg-Gc,cor,rel)/Gc,cor,rel
[0147] SFc.Eff = nc,deg/nc
[014引 A SFc,Eff = (nc,deg-nc)Alc
[0149] 其中SFC,FC为压气机流量特性指数;Ge, cor,rel, deg为压气机(性能衰退时)相对折合流 量;Ge, cor, rel为压气机(健康时）相对折合流量；SFc, Ef f为压气机效率特性指数;nc, deg为压气 机(性能衰退时)等赌效率;化为压气机(健康时)等赌效率。
[0150] 燃烧室气路健康指数定义如下：
[0151] SF B, Ef f = n,B, deg/n,B
[0152] A SFe'Eff 二（HB'deg-riB)/nB
[0153] 其中SFb,Eff为燃烧室燃烧效率性能指数;riB, deg为燃烧室(性能衰退时)燃烧效率;％ 为燃烧室(健康时)燃烧效率。
[0154] 透平气路健康指数定义如下：
[0155] SFt , FC = Gt , cor, deg/Gx, cor
[0156] A SFt.FC= (GT,cor,deg-GT,cor)/Gx,cor
[0157] SFt, Eff = HTidegAlT
[015引 A SF^Eff= (TlMeg-TlT)AlT
[0159] 其中SFt,FG为透平流量性能指数;GT,eDr,deg为透平(性能衰退时）折合流量;GT,cDr为 透平(健康时）折合流量;SFt, Ef f为透平效率性能指数;化deg为透平(性能衰退时）等赌效率； %为透平(健康时)等赌效率。
[0160] 通过设置热力模型中的各个部件(压气机、透平和燃烧室）的气路健康指数，来模 拟单部件和多部件性能衰退时的气路测量参数，用于积累性能衰退模式和衰退征兆(整理 为性能衰退时的气路测量参数相对于健康时的相对偏差形式A三）的知识数据库，作为灰色 关联理论实现对实测气路参数进行模式识别的依据。
[0161]
[0162] 其中，三为通过设置热力模型中的各个部件(压气机、透平和燃烧室）的气路健康指 数来模拟单部件和多部件性能衰退时的气路参数计算值，；为在相同环境条件和操作条件 下部件性能健康时的热力模型气路参数计算值，Ai为整理为性能衰退时的气路测量参数相 对于健康时的相对偏差形式。
[0163] 其中灰色关联理论中的知识数据库形式为：
[0164]
[0165] 其中C^j = I，2,…m)为已知的部件性能衰退模式;m为总的已知的部件性能衰退 模式数目。
[0166] 采集当前对象燃气轮机稳定运行时的某一时段的气路测量参数，进行降噪处理后 作为待离线诊断的气路测量参数，设置已建立的燃气轮机热力模型的环境输入条件(大气 压力、溫度和相对湿度)和操作输入条件与采样时的对象燃气轮机运行工况一致，消除由于 环境条件(大气压力、溫度和相对湿度)和操作条件变化而导致燃气轮机运行性能变化的影 响。将待离线诊断的实测气路数据与热力模型计算的气路参数数据之间的相对偏差输入灰 色关联理论算法(GRA)中实现性能衰退模式识别，隔离已发生性能衰退的部件，从而实现故 障系数矩阵的降维（即确定了待诊断的实际气路部件健康指数的子集），如图3所示。
[0167] 图4中三为待离线诊断的实测气路数据，1^为对象燃气轮机环境输入条件（大气压 力、溫度和相对湿度)和操作输入条件，么三为待离线诊断的实测气路数据;^与热力模型（气 路部件健康时）计算的气路参数数据^之间的相对偏差。识别、隔离发生性能衰退的部件 后，即可得到部件健康指数A添的子集(A涼)、,,6,从而大大降低故障系数矩阵的维数。
[0168] W待离线诊断的气路测量数据^与热力模型计算的气路参数数据三之间的偏差作 为残差，通过牛顿-拉普森迭代算法计算得到当前的已隔离部件(压气机、透平和燃烧室）的 气路健康指数(A赤U，用W评估对象燃气轮机实际的性能健康状况，如图5所示。
[0169] 运里I=；-三。由于燃气轮机热力性能的强非线性特性，一个迭代计算过程(牛顿- 拉普森算法)用于计算部件健康指数掉恋kb。随着迭代，当预测的气路测量参数非常接近实 际的测量参数时（即，满足||^||<€，e是设定的一个相对较小的阔值），则可得到最终的部件 健康指数(A承U，如图6所示。
[0170] 从图6可知，基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部件自 适应性能诊断方法(GRA-GPA)能够有效地识别、隔离实际性能衰退部件，并能准确地预测性 能衰退程度。然而，传统诊断方法(GPA)诊断的结果可能无法接受。基于W上对本发明方法 的诊断案例分析可知，当部件发生性能退化时，由于存在气路测量噪音，且燃气轮机中的部 件数目较多，通过传统诊断方法(GPA)诊断会出现一定程度的模糊效应和误导性性能衰退 预测。而通过基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部件自适应性能 诊断方法(GRA-GPA)诊断则能有效识别、隔离性能衰退的部件，并准确地预测性能衰退程 度。
[0171] 本发明的核屯、技术内容在于为更好解决由于测量噪音和部件数目较多导致诊断 结果的模糊效应问题，借鉴模式识别技术的优点，从故障系数矩阵降维的角度来改善诊断 结果的准确性，提出了一种基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部 件自适应性能诊断方法。根据燃气轮机气动热力学特性，利用部件相对折合参数，建立燃气 轮机部件级热力模型，其中压气机和透平特性线数据整理成通用的相对折合参数形式，相 比于现有技术，能更简单方便地用于气路诊断时设置各个气路部件健康指数；用相对折合 参数重新定义压气机和透平的气路健康指数(代表了部件特性线的偏移），相比于现有技 术，能更准确地表征由于部件性能衰退而导致的气路健康指数的变化，消除由于环境条件 (大气压力、溫度和相对湿度)变化而给诊断结果带来的负影响;通过设置热力模型中的各 个部件(压气机、透平和燃烧室）的气路健康指数，来模拟单部件和多部件性能衰退时的气 路测量参数，用于积累性能衰退模式和衰退征兆(整理为性能衰退时的气路测量参数相对 于健康时的相对偏差形式）的知识数据库，作为灰色关联理论实现对实测气路参数进行模 式识别的依据，相比于现有技术，既能节省高昂的试验成本，还能准确详尽地积累各种性能 衰退模式和衰退征兆的知识数据库;通过灰色关联理论算法实现性能衰退模式识别，从而 隔离实际真正发生性能衰退的部件，从而确定了待诊断的实际气路部件健康指数的子集， 从故障系数矩阵降维的角度来改善诊断结果的准确性，相比于现有技术，能更好解决由于 测量噪音和部件数目较多导致诊断结果出现模糊效应的问题。对于保证燃气轮机安全可靠 运行、降低维护费用有着积极促进作用。
【主权项】
1. 基于灰色关联理论与热力模型相结合的混合型燃气轮机气路部件自适应性能诊断 方法，其特征是： (1) 基于对象燃气轮机健康时的气路测量参数，建立燃气轮机非线性热力模型，压气机 和透平都用相似折合参数形式表示； (2) 用相似折合参数重新定义压气机和透平的气路健康指数，消除由于环境条件变化 而导致燃气轮机运行性能变化的影响； (3) 通过设置热力模型中各个部件的气路健康指数，来模拟单部件和多部件性能衰退 时的气路测量参数，用于积累性能衰退模式和衰退征兆的知识数据库，作为灰色关联理论 实现对实测气路参数进行模式识别的依据，所述热力模型的部件包括压气机、透平和燃烧 室； (4) 采集当前燃气轮机稳定运行时的一个时段的气路测量参数，进行降噪处理后作为 待离线诊断的气路测量参数； (5) 设置已建立的燃气轮机热力模型的环境输入条件和操作输入条件与采样时的对象 燃气轮机运行工况一致，消除由于环境条件和操作条件变化而导致燃气轮机运行性能变化 的影响； (6) 将待离线诊断的气路测量参数与热力模型的气路参数数据之间的相对偏差输入灰 色关联理论算法中实现性能衰退模式识别，隔离已发生性能衰退的部件，从而实现故障系 数矩阵的降维； (7) 以待离线诊断的气路测量数据与热力模型的气路参数数据之间的偏差作为残差， 通过迭代算法计算得到当前压气机、透平和燃烧室的气路健康指数，用以评估对象燃气轮 机实际的性能健康状况。2. 根据权利要求1所述的基于灰色关联理论与热力模型相结合的混合型燃气轮机气路 部件自适应性能诊断方法，其特征是:步骤(1)中建立燃气轮机非线性热力模型的具体步骤 如下： (a)压气机特性线整理成通用的相似折合参数形式如下：%相对折合转速，η为实际转速，$为压气机进口滞止温 度，心为流经压气机工质的气体常数，下角标O表示设计点；为相对折合流量，Ge为实际压气机进口流量，IC为压 气机进口滞止压力/3相对压比，叱为实际压气机压比，= 为相对等熵 效率，ne为实际压气机等熵效率； 透平特性线整理成通用的相似折合参数形式如下： Gt, cor,rel - f (Ilcor, rel，冗T,rel) Π .Τ, rel - f (Ilcor，rel，冗T, rel) 式中％相对折合转速，η为实际转速，< 为透平进口滞止温 度，Rg为流经透平工质的气体常数，I相对折合流量，Gt 为实际透平进口流量，7r^rel - "为相对压比，％,ι·θι=ητ/%()为相对等烟效率，％为实际透平 v^TO 等熵效率，下角标〇表示设计点； (b)根据采集的燃气轮机健康时的气路测量参数，修正各个部件的特性线数据，使所建 的燃气轮机热力模型的计算值与气路实测参数相匹配，从而消除热力模型计算误差给诊断 结果带来的负影响。3. 根据权利要求2所述的基于灰色关联理论与热力模型相结合的混合型燃气轮机气路 部件自适应性能诊断方法，其特征是:步骤(2)用相似折合参数重新定义压气机和透平的气 路健康指数的具体步骤如下： 压气机气路健康指数定义如下： SFc，FC - Ge, cor, rel, deg/Gc, cor,rel Δ SFc,FC - (Ge, cor, rel, deg-Ge, cor, rel ) /Ge, cor,rel SFc, Eff - He, deg/^C Δ SFcjEff= (ric,deg-ric)/ric 其中SFC,FC为压气机流量特性指数;Ge, 为压气机性能衰退时的相对折合流量； Gn, rel为压气机健康时的相对折合流量;SFe, Eff为压气机效率特性指数;ne,deg为压气机性 能衰退时的等熵效率;取为压气机健康时的等熵效率； 透平气路健康指数定义如下： SFt, FC - Gt, cor, deg/GT, cor Δ SFt,fc - (Gt, c〇r, deg-Gt, c〇r)/Gt ,cor SFTjEff = Ht, deg/ητ Δ SFTjEff= (rIT, deg-ητ) At 其中SFt^s透平流量性能指数;Gt, CC^deg为透平性能衰退时的折合流量;Gt, 为透平 健康时的折合流量;SFt,Eff为透平效率性能指数;nT,deg为透平性能衰退时的等熵效率;ητ为 透平健康时的等熵效率。4. 根据权利要求3所述的基于灰色关联理论与热力模型相结合的混合型燃气轮机气路 部件自适应性能诊断方法，其特征是:步骤(3)积累性能衰退模式和衰退征兆的知识数据库 作为灰色关联理论实现对实测气路参数进行模式识别的依据的具体步骤如下： (a)通过设置热力模型中各个部件的气路健康指数，来模拟单部件和多部件性能衰退 时的气路测暈参数，并整理为性能衰退时的气路测量参数相对于健康时的相对偏差形式：其中，?为通过设置热力模型中的各个部件的气路健康指数来模拟单部件和多部件性 能衰退时的气路参数计算值，为在相同环境条件和操作条件下部件性能健康时的热力模 型气路参数计算值，Ai为整理为性能衰退时的气路测量参数相对于健康时的相对偏差形 式； (b)积累性能衰退模式和衰退征兆的知识数据库作为灰色关联理论实现对实测气路参 数进行模式识别的依据，其中灰色关联理论中的知识数据库形式为：其中(^(」=1，2，一!11)为已知的部件性能衰退模式;(^(」=1，2，一）为每一个气路可测参 数;m为总的已知的部件性能衰退模式数目；η为总的对象燃气轮机气路可测参数数目。5. 根据权利要求4所述的基于灰色关联理论与热力模型相结合的混合型燃气轮机气路 部件自适应性能诊断方法，其特征是:步骤(6)将待离线诊断的气路测量参数与热力模型的 气路参数数据之间的相对偏差输入灰色关联理论算法中实现性能衰退模式识别，隔离已发 生性能衰退的部件的具体步骤如下： (a) 将待离线诊断的气路测量参数与热力模型的气路参数数据之间的相对偏差输入灰 色关联理论算法中实现性能衰退模式识别： 对于 pe(〇，l):其中P为分辨系数，X1S-组待离线诊断的实测气路数据，I(X1U)，Cj(k))为XjPQ的第 k个特征参数间的灰色关联系数，IU1, Cj)为XdPCj的灰色关联度。 Xi被分类至当ξ(Xi，Cj) (j = 1，2，…，m)取得最大值时所属的部件性能衰退模式类； (b) 隔离已发生性能衰退的部件，从而实现故障系数矩阵的降维，即确定待诊断的实际 气路部件健康指数人尿的子集。6. 根据权利要求5所述的基于灰色关联理论与热力模型相结合的混合型燃气轮机气路 部件自适应性能诊断方法，其特征是:步骤(7)以待离线诊断的气路测量参数与热力模型的 气路参数数据之间的偏差作为残差，通过迭代算法计算得到当前的各个部件的气路健康指 数，用以评估对象燃气轮机实际的性能健康状况的具体步骤如下： (a) 以待离线诊断的气路测量参数i与热力模型计算的气路参数数据纟之间的偏差作为 残差通过牛顿-拉普森迭代算法进行迭代，当预测的气路测量参数与实际的测量 参数之差小于阈值，即Ρ|<^，ε是设定的阈值，则得到最终的部件健康指数(Δ远 (b) 根据预测的已隔离部件的性能衰退程度，评估当前对象燃气轮机总体性能健康状 况。
【文档编号】G06F17/50GK106021757SQ201610362580
【公开日】2016年10月12日
【申请日】2016年5月26日
【发明人】李淑英, 应雨龙, 曹云鹏
【申请人】哈尔滨工程大学