一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方法
【专利摘要】本发明的目的在于提供一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方法，建立燃气轮机非线性热力模型，采集当前对象燃气轮机稳定运行时的某一时段的气路测量参数，对其进行数据调和，检测出发生性能衰退的气路传感器，逐一采集当前对象燃气轮机在各个稳定运行工况时的某一时段的气路测量参数，进行降噪处理，消除由于环境条件和操作条件变化而导致燃气轮机运行性能变化的影响，计算得到当前的各个部件的气路健康指数，用以评估对象燃气轮机实际的性能健康状况。本发明解决了传统燃气轮机气路性能诊断方法诊断准确性高度依赖于气路传感器可靠性的问题，能有效适用于存在测量噪音、测量偏差及复杂燃气轮机机组的气路部件性能诊断情况。
【专利说明】
一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断 方法
技术领域
[0001] 本发明涉及的是一种诊断方法，具体地说是燃气轮机的性能诊断方法。
【背景技术】
[0002] 在20世纪下半叶，随着燃气轮机在航空工业中广泛应用，越来越多的受到工业电 站领域、石油和天然气管道运输以及舰船工业领域的关注。在运行中，由于高温、高压、高转 速及高应力的恶劣工作条件及环境污染等影响，各种类型的燃气轮机都会逐渐性能衰退。 燃气轮机的主要气路部件包含压气机、燃烧室和透平。这些主要部件会随着时间遭受不同 的退化现象，如污垢、泄漏、腐蚀、热畸变、外来物损坏等，将会引起性能恶化并易导致各种 严重的故障发生运行安全问题。因此，对于燃气轮机用户来说，当前燃气轮机的性能健康状 况是非常重要的信息。
[0003] 目前大多数燃气轮机的维修策略是预防性维修保养，通常按照燃机制造商指示的 当量运行小时数(Ε0Η)来考虑是否需要小修、中修、大修。燃机停运，无论是计划内或计划外 的，总是意味着昂贵的成本代价。为了节省维修费用，用户需要根据燃机实际的性能健康状 况采取维修策略，即预测性维修保养。
[0004] 基于热力模型决策的气路部件性能诊断方法已经广泛应用于燃气轮机性能健康 状态监测中，并且已经成为支持维修策略改革的关键技术之一。
[0005] 理论上，气路分析方法方法(特别是非线性气路分析方法)能够容易地预测部件性 能衰退程度。通常，准确的测量信息对获取准确的衰退征兆从而得到准确的气路性能诊断 结果是至关重要的。然而，随着燃气轮机运行，气路传感器同部件一样，有可能会发生性能 衰退，甚至故障，此时传感器会产生显著的测量偏差(不符合高斯分布），会导致误导性的诊 断结果。
【发明内容】
[0006] 本发明的目的在于提供能有效适用于存在测量噪音、测量偏差的复杂燃气轮机机 组的一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方法。
[0007] 本发明的目的是这样实现的：
[0008] 本发明一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方法，其特征 是：
[0009] (1)基于对象燃气轮机新投运或健康时的气路测量参数，建立燃气轮机非线性热 力模型，压气机和透平都用相似折合参数形式表示；
[0010] (2)用相似折合参数重新定义压气机和透平的气路健康指数，消除由于环境条件 变化而导致燃气轮机运行性能变化的影响；
[0011] (3)采集当前对象燃气轮机在一运行工况下稳定运行时的一个时段的气路测量参 数，进行降噪处理后作为待离线诊断的气路测量参数；
[0012] (4)对待离线诊断的气路测量参数进行数据调和，检测出发生性能衰退的气路传 感器；
[0013] (5)根据检测出存在问题的传感器数目，按照可信的气路测量参数2 (待诊断的部 件气路健康指数数目+存在问题传感器的数目）原则，确定对象燃气轮机部件性能诊断所需 的稳态运行工况点的数目；
[0014] (6)根据所需的稳态运行工况点数目，通过调整对象燃气轮机操作条件的方式，逐 一采集当前对象燃气轮机在各个稳定运行工况时的一个时段的气路测量参数，进行降噪处 理后作为多个运行工况点的待离线诊断的气路测量数据集；
[0015] (7)分别设置已建立的燃气轮机热力模型的环境输入条件和操作输入条件与采样 时的对象燃气轮机的各个运行工况一致，消除由于环境条件和操作条件变化而导致燃气轮 机运行性能变化的影响；
[0016] (8)以待离线诊断的多个运行工况点的气路测量数据与热力模型的气路参数数据 之间的均方根误差为目标函数，通过粒子群优化算法计算得到当前的各个部件的气路健康 指数，用以评估对象燃气轮机实际的性能健康状况。
[0017] 本发明还可以包括：
[0018] 1、步骤(1)中建立燃气轮机非线性热力模型的具体步骤如下：
[0019] (a)利用部件相对折合参数，建立燃气轮机部件级热力模型，其中压气机和透平特 性线数据整理成通用的相对折合参数形式：
[0020]压气机特性线整理成通用的相似折合参数形式如下：
[0021 ] Gc, cor, rel - f (Ilcor,rel, ^C, rel)
[0022] HC，rel - f (ilcor, rel, ^C, rel )
[0023]
为相对折合转速，n为实际转速，2：为压气机进口滞 止温度，Rg为流经压气机工质的气体常数，下角标〇表示设计点；
[0024]
为相对折合流量，GC为实际压气机进口流量，Λ: 为压气机进口滞止压力
为相对压比，叱为实际压气机压比，η^ι=%/η〇)为相对 等熵效率，ne为实际压气机等熵效率；
[0025]透平特性线整理成通用的相似折合参数形式如下：
[0026] Gt, cor, rel - f (Ilcor,rel，冗T, rel)
[0027] Π.Τ, rel - f (ilcor, rel, rel )
[0028]
为相对折合转速，n为实际转速，ζ为透平进口滞 止温度，Rg为流经透平工质的气体常数：
为相对折合流 量，Gt为实际透平进口流量
为相对压比，Πτ, rel = %/ητο为相对等熵效率，Πτ为实 际透平等熵效率，下角标0表示设计点；
[0029] (b)根据采集的对象燃气轮机新投运或健康时的气路测量参数，逐步修正各个部 件的特性线数据，使所建的燃气轮机热力模型的计算值与气路实测参数相匹配，从而消除 热力模型计算误差给诊断结果带来的负影响。
[0030] 2、步骤(2)用相似折合参数重新定义压气机和透平的气路健康指数的具体步骤如 下：
[0031] (a)燃气轮机总体性能健康状况通过由各部件的气路健康指数来表示，用相对折 合参数重新定义压气机和透平的气路健康指数以消除由于环境条件变化而给诊断结果带 来的负影响。
[0032] (b)压气机气路健康指数定义如下：
[0033] SFc, FC - Gc, cor,rel, deg/Gc, cor, rel
[0034] Δ SFc, FC - (Gc, cor,rel, deg-Gc, cor, rel )/Gc, cor, rel
[0035] SFc.Eff = ric,deg/ric
[0036] Δ SFc.Eff = (nc,deg-nc)/nc
[0037] 其中SFc,fc为压气机流量特性指数;6^。^1,<^为压气机性能衰退时相对折合流 量;为压气机健康时相对折合流量;SFe, Eff为压气机效率特性指数;ne,deg为压气机 性能衰退时等熵效率;ne为压气机健康时等熵效率；
[0038]燃烧室气路健康指数定义如下：
[0039] SFB.Eff = riB,deg/riB
[0040] Δ SFB，Eff = (riB，deg-riB) Aib
[0041] 其中SFB,Eff为燃烧室燃烧效率性能指数;nB,deg为燃烧室性能衰退时燃烧效率;q B为 燃烧室健康时燃烧效率。
[0042] 透平气路健康指数定义如下：
[0043] SFT,FC = GT,cor,deg/GT,cor
[0044] Δ SFt, FC -(Gt, cor, deg-GT, cor )/Gt, cor
[0045] SFT,Eff = riT,degAlT
[0046] Δ SFT.Eff = (Πτ,deg-ητ) At
[0047] 其中SFt,fc为透平流量性能指数;GT,ccir,deg为透平性能衰退时折合流量;G T,ccir为透 平健康时折合流量;SFT,Eff为透平效率性能指数;％,<^为透平性能衰退时等熵效率;ητ为透 平健康时等熵效率。
[0048] 3、步骤(4)基于高斯修正准则的数据调和原理对该工况下的待离线诊断的气路测 量数据进行数据调和，检测出发生性能衰退的气路传感器的具体步骤如下：
[0049] (a)在一个运行工况下，采集当前对象燃气轮机稳定运行时的一个时段的气路测 量参数，进行降噪处理后作为待离线诊断的气路测量参数；
[0050] (b)在对气路测量参数进行调和时，将气路部件的流量特性指数和效率特性指数 引入作为"虚拟"测量参数，一同进行调和，其"虚拟"测量参数都取为1，其不确定度取为 1%；
[0051] (c)基于高斯修正准则的数据调和原理对该工况下的待离线诊断的气路测量数据 及"虚拟"测量参数进行数据调和，检测出发生性能衰退的气路传感器。
[0052] 4、步骤(8)以待离线诊断的多个运行工况点的气路测量数据与热力模型计算的气 路参数数据之间的偏差作为残差，计算得到当前的各个部件的气路健康指数，用以评估对 象燃气轮机实际的性能健康状况的具体步骤如下：
[0053] (a)分别设置已建立的燃气轮机热力模型的环境输入条件和操作输入条件与采样 时的对象燃气轮机的各个运行工况一致，消除由于环境条件和操作条件变化而导致燃气轮 机运行性能变化的影响；
[0054] (b)以待离线诊断的多个运行工况点的气路测量数据i与热力模型计算的气路参 数数据2;之间的均方根误差为目标函数Fitness，通过粒子群优化算法计算得到当前的各个 部件的气路健康指数Δ互 ：，用以评估对象燃气轮机实际的性能健康状况，
[0056] 式中z，为由燃机热力模型在第i个运行工况点计算得到的气路测量参数向量，i为 在第i个运行工况点实际的气路测量参数向量，Μ为某一运行工况点的气路测量参数的总数 目；
[0057] 优化目标函数定义如式：
[0059] 式中Fitness是优化目标函数，当Fitness随着粒子群优化算法(PS0)迭代寻优计算 过程逐渐地趋近于〇时，预测的气路测量参数[為， ζ2，.··，4,]与实测气路测量参数 相匹配，此时输出最优的气路部件健康指数Δ：^和在不同运行工况点的传感器偏差 \t)\ , · > ihm J ο
[0060] 本发明的优势在于：
[0061] (1)本发明根据燃气轮机气动热力学特性，利用部件相对折合参数，建立燃气轮机 部件级热力模型，其中压气机和透平特性线数据整理成通用的相对折合参数形式，相比于 现有技术，能更简单方便地用于气路诊断时设置各个气路部件健康指数。
[0062] (2)本发明用相对折合参数重新定义压气机和透平的气路健康指数(代表了部件 特性线的偏移），相比于现有技术，能更准确地表征由于部件性能衰退而导致的气路健康指 数的变化，消除由于环境条件(大气压力、温度和相对湿度)变化而给诊断结果带来的负影 响。
[0063] (3)本发明利用高斯修正准则对所有气路测量参数进行数据调和，检测出有问题 的传感器，并基于多运行工况点的气路测量参数采用粒子群优化算来识别、隔离性能衰退 的部件，并进一步量化气路传感器和部件的性能衰退程度，得到精确的诊断结果。解决了传 统燃气轮机气路性能诊断方法诊断准确性高度依赖于气路传感器可靠性的问题，能有效适 用于存在测量噪音、测量偏差的复杂燃气轮机机组的性能诊断情况。
【附图说明】
[0064]图1为本发明的流程图；
[0065]图2为本发明的诊断过程示意图；
[0066]图3某型三轴船用燃气轮机气路工作截面标识图；
[0067] 图4某型三轴船用燃气轮机诊断案例的诊断结果；
[0068] 图5该诊断案例的本发明算法迭代计算搜索过程。
【具体实施方式】
[0069] 下面结合附图举例对本发明做更详细地描述：
[0070] 结合图1-5,本发明一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方 法，包括以下步骤：
[0071] 步骤1)，基于对象燃气轮机新投运(或健康)时的气路测量参数建立能完全反映各 个部件特性的燃气轮机非线性热力模型，其中压气机和透平都用相似折合参数形式表示；
[0072] 步骤2)，用相似折合参数重新定义压气机和透平的气路健康指数，消除由于环境 条件(大气压力、温度和相对湿度)变化而导致燃气轮机运行性能变化的影响；
[0073] 步骤3)，在某一运行工况下，采集当前对象燃气轮机稳定运行时的某一时段的气 路测量参数，进行降噪处理后作为待离线诊断的气路测量参数；
[0074]步骤4)，基于高斯修正准则的数据调和原理对该工况下的待离线诊断的气路测量 数据进行数据调和，检测出发生性能衰退的气路传感器；
[0075] 步骤5)，根据检测出存在问题的传感器数目，按照可信的气路测量参数>=(待诊 断的部件气路健康指数数目+存在问题传感器的数目）原则，确定对象燃气轮机部件性能诊 断所需的稳态运行工况点的数目；
[0076] 步骤6)，根据所需的稳态运行工况点数目，通过调整对象燃气轮机操作条件的方 式，逐一采集当前对象燃气轮机在各个稳定运行工况时的某一时段的气路测量参数，进行 降噪处理后作为多个运行工况点的待离线诊断的气路测量数据集；
[0077]步骤7)，分别设置已建立的燃气轮机热力模型的环境输入条件(大气压力、温度和 相对湿度)和操作输入条件与采样时的对象燃气轮机的各个运行工况一致，消除由于环境 条件(大气压力、温度和相对湿度)和操作条件变化而导致燃气轮机运行性能变化的影响；
[0078] 步骤8)，以待离线诊断的多个运行工况点的气路测量数据与热力模型计算的气路 参数数据之间的均方根误差为目标函数，通过粒子群优化算法计算得到当前的各个部件 (压气机、透平和燃烧室）的气路健康指数，用以评估对象燃气轮机实际的性能健康状况。
[0079] 作为本发明一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方法进一 步的优化方案，步骤1)中所述基于对象燃气轮机新投运(或健康)时的气路测量参数建立能 完全反映各个部件特性的燃气轮机非线性热力模型的具体步骤如下：
[0080] 步骤1.1)，根据燃气轮机气动热力学特性，利用部件相对折合参数，建立燃气轮机 部件级热力模型，其中压气机和透平特性线数据整理成通用的相对折合参数形式。
[0081 ]其中压气机特性线整理成通用的相似折合参数形式如下：
[0082] Gc, cor, rel - f (Ilcor,rel, ^C, rel)
[0083] Π-C, rel - f (llcor, rel, ^C, rel )
[0084]
为相对折合转速，n为实际转速，为压气机进口 滞止温度，心为流经压气机工质的气体常数，下角标0表示设计点；
[0085]
I相对折合流量，Ge为实际压气机进口流量，g 为压气机进口滞止压力，
为相对压比，取为实际压气机压比，nc, rei=nc/n〇)为相对 等熵效率，ne为实际压气机等熵效率。
[0086]透平特性线整理成通用的相似折合参数形式如下：
[0087] Gt, cor, rel - f (Ilcor,rel，町,rel)
[0088] Π.Τ, rel - f (llcor, rel, ^·?, rel )
[0089]
为相对折合转速，n为实际转速，< 为透平进口滞 止温度，心为流经透平工质的气体常数，
为相对折合流 量，Gt为实际透平进口流量，
为相对压比，Πτ, rel = Πτ/ητο为相对等熵效率，Πτ为实 际透平等熵效率，下角标0表示设计点。
[0090] 步骤1.2)，根据采集的对象燃气轮机新投运(或健康)时的气路测量参数（降噪处 理后），如总温、总压、转速等，逐步修正各个部件的特性线数据(包括设计工况和变工况）， 使所建的燃气轮机热力模型的计算值与气路实测参数相匹配，从而消除热力模型计算误差 给诊断结果带来的负影响。
[0091] 作为本发明一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方法进一 步的优化方案，步骤2)中所述用相似折合参数重新定义压气机和透平的气路健康指数的具 体步骤如下：
[0092] 步骤2.1 )，燃气轮机总体性能健康状况可以由各主要部件的气路健康指数，如压 气机和透平的流量特性指数和效率特性指数、燃烧室的效率特性指数来表示。这里用相对 折合参数重新定义压气机和透平的气路健康指数(代表了部件特性线的偏移），以消除由于 环境条件(大气压力、温度和相对湿度)变化而给诊断结果带来的负影响。
[0093] 步骤2.2)，压气机气路健康指数定义如下：
[0094] SFc, FC - Gc, cor,rel, deg/Gc, cor, rel
[0095] Δ SFc, FC - (Gc, cor，rel, deg_Gc, cor, rel )/Gc, cor, rel
[0096] SFc.Eff = nc,deg/ric
[0097] Δ SFc.Eff = (ric,deg-ric)/ric
[0098] 其中SFC,FC为压气机流量特性指数;6^_#1,<^为压气机(性能衰退时)相对折合流 量;Ge, , rel为压气机(健康时)相对折合流量;SFe, Ef f为压气机效率特性指数;ne,deg为压气 机(性能衰退时)等熵效率;ne为压气机(健康时)等熵效率。
[0099] 燃烧室气路健康指数定义如下：
[0100] SFB.Eff = riB,deg/riB
[0101] Δ SFB.Eff = (riB.deg-%)/%
[0102] 其中SFb, Eff为燃烧室燃烧效率性能指数;nB, ^为燃烧室(性能衰退时)燃烧效率;nB 为燃烧室(健康时)燃烧效率。
[0103] 透平气路健康指数定义如下：
[0104] SFT,FC = GT,cor,deg/GT,cor
[0105] Δ SFt,FC= (GT,cor,deg-GT,cor)/GT,cor
[0106] SFT,Eff = riT, deg/ητ
[0107] Δ SFT.Eff = (Πτ,deg-ητ) At
[0108] 其中SFt,fc为透平流量性能指数;GT,ccir,deg为透平(性能衰退时）折合流量;G T,ccir为 透平(健康时)折合流量;sft, Eff为透平效率性能指数;ητ, deg为透平(性能衰退时)等熵效率； ητ为透平(健康时)等熵效率。
[0109] 作为本发明一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方法进一 步的优化方案，步骤4)中所述基于高斯修正准则的数据调和原理对该工况下的待离线诊断 的气路测量数据进行数据调和，检测出发生性能衰退的气路传感器的具体步骤如下：
[0110] 步骤4.1)，在某一运行工况下，采集当前对象燃气轮机稳定运行时的某一时段的 气路测量参数，进行降噪处理后作为待离线诊断的气路测量参数；
[0111] 步骤4.2)，在对气路测量参数进行调和时，实际的气路部件特性也可能会与理论 特性发生一定偏差，这里将气路部件的流量特性指数和效率特性指数引入作为"虚拟"测量 参数(部件健康时，其值为1)，一同进行调和，由于未知气路部件性能健康状况，因此其"虚 拟"测量参数都取为1，其不确定度取为1 %。
[0112]步骤4.3)，基于高斯修正准则的数据调和原理对该工况下的待离线诊断的气路测 量数据及"虚拟"测量参数进行数据调和，检测出发生性能衰退的气路传感器（即调节量超 过置信限值的传感器）。
[0113]作为本发明一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方法进一 步的优化方案，步骤8)中所述以待离线诊断的多个运行工况点的气路测量数据与热力模型 计算的气路参数数据之间的偏差作为残差，通过牛顿-拉普森迭代算法计算得到当前的各 个部件(压气机、透平和燃烧室）的气路健康指数，用以评估对象燃气轮机实际的性能健康 状况的具体步骤如下：
[0114]步骤8.1)，分别设置已建立的燃气轮机热力模型的环境输入条件(大气压力、温度 和相对湿度)和操作输入条件与采样时的对象燃气轮机的各个运行工况一致，消除由于环 境条件（大气压力、温度和相对湿度)和操作条件变化而导致燃气轮机运行性能变化的影 响。
[0115]步骤8.2)，以待离线诊断的多个运行工况点的气路测量数据?与热力模型计算的 气路参数数据】，之间的均方根误差为目标函数Fitness，通过粒子群优化算法计算得到当前 的各个部件(压气机、透平和燃烧室）的气路健康指数Λ夺 ：用以评估对象燃气轮机实际的 性能健康状况，如图2所不。
[0117] 式中&--由燃机热力模型在第i个运行工况点计算得到的气路测量参数向量；
[0118] ?--在第i个运行工况点实际的气路测量参数向量；
[0119] Μ一一某一运行工况点的气路测量参数的总数目。
[0120] 优化目标函数定义如式：
[0122]式中Fitness是优化目标函数，当Fitness随着粒子群优化算法(PS0)迭代寻优计算 过程逐渐地趋近于0时，预测的气路测量参数[^2，…,4]与实测气路测量参数ι?,?，...,ζ] 相匹配，此时输出最优的气路部件健康指数Δ；^和在不同运行工况点的传感器偏差
[ii，U",]。
[0123] 本发明的核心技术内容在于解决了传统燃气轮机气路部件性能诊断方法诊断精 度易受环境条件(大气压力、温度和相对湿度)及操作条件变化影响的问题，并解决了传统 燃气轮机气路性能诊断方法诊断准确性高度依赖于气路传感器可靠性的问题，能有效适用 于存在测量噪音和复杂燃气轮机机组的性能诊断情况。根据燃气轮机气动热力学特性，利 用部件相对折合参数，建立燃气轮机部件级热力模型，其中压气机和透平特性线数据整理 成通用的相对折合参数形式，相比于现有技术，能更简单方便地用于气路诊断时设置各个 气路部件健康指数；用相对折合参数重新定义压气机和透平的气路健康指数(代表了部件 特性线的偏移），相比于现有技术，能更准确地表征由于部件性能衰退而导致的气路健康指 数的变化，消除由于环境条件(大气压力、温度和相对湿度)变化而给诊断结果带来的负影 响;利用高斯修正准则对所有气路测量参数进行数据调和，检测出有问题的传感器，并基于 多运行工况点的气路测量参数采用粒子群优化算来识别、隔离性能衰退的部件，并进一步 量化气路传感器和部件的性能衰退程度，得到精确的诊断结果。解决了传统燃气轮机气路 性能诊断方法诊断准确性高度依赖于气路传感器可靠性的问题，能有效适用于存在测量噪 音、测量偏差的复杂燃气轮机机组的性能诊断情况。
[0124] 以某型三轴船用燃气轮机的气路部件健康诊断为例，其该型三轴船用燃气轮机气 路工作截面标识图如图3所示。该型三轴燃气轮机包括两个压气机(即一个低压压气机(LC) 和一个高压压气机(HC))、一个燃烧室(Β)和三个透平（即一个高压透平(ΗΤ)、一个低压透平 (LT)和一个动力透平(ΡΤ))，其中发电机通过一个减速齿轮箱与动力透平(ΡΤ)相连接。低压 透平(LT)的输出功通过低压轴驱动低压压气机(LC)来压缩从进气道出来的空气，高压透平 (ΗΤ)的输出功通过高压轴驱动高压压气机(HC)来继续压缩从低压压气机(LC)出来的空气。 从高压压气机(HC)出来的高压空气进入燃烧室(Β)与燃料发生燃烧化学反应生成高温、高 压的燃气，燃气依次进入高压透平(ΗΤ)、低压透平(LT)和动力透平(ΡΤ)来驱动透平输出功。 最终，动力透平(ΡΤ)通过减速齿轮箱驱动发电机来产生电功率。同时，从压气机中抽取的冷 却空气流入热端气流通道去冷却各个透平前几级的静叶、动叶和轮盘。当燃气轮机稳定运 行时，发电机的电功率和动力透平(PT)的转速通常作为主要控制参数而维持定常。该机组 的气路测量参数如表1所示，各个气路部件的健康指数如表2所示。
[0125]表1该型燃气轮机机组的气路测量参数
[0127] 表2主要气路部件的健康指数
[0128]
[0129] 基于该型燃气轮机新投运(或健康）时的气路测量参数建立能完全反映各个部件 特性的燃气轮机非线性热力模型，其中压气机和透平都用相对折合参数形式表示。
[0130] 其中压气机特性线整理成通用的相似折合参数形式如下：
[01 31 ] Gc, cor, rel - f (Hcor,rel, -Π-C, rel)
[01 32] He, rel - f (llcor, rel, -Π-C, rel )
[0133]
胃为相对折合转速，n为实际转速，ζ为压气机进口 滞止温度，心为流经压气机工质的气体常数，下角标0表示设计点；
[0134]
.为相对折合流量，Ge为实际压气机进口流量，< 为压气机进口滞止压力，
为相对压比，取为实际压气机压比，nc, rei=nc/n〇)为相对 等熵效率，ne为实际压气机等熵效率。
[0135] 透平特性线整理成通用的相似折合参数形式如下：
[0136] GT，cor，rel -f (Ilcor，rel，兀T ，rel )
[01 37] Γ[Τ, rel - f (ricor, rel , ^T, rel )
[0138]
「为相对折合转速，n为实际转速，<为透平进口滞 止温度，Rg为流经透平工质的气体常数
为相对折合流 量，Gt为实际透平进口流量，
为相对压比，ΠΤ, rel = %/ηΤΟ为相对等熵效率，ΠΤ为实 际透平等熵效率，下角标〇表示设计点。
[0139] 对象燃气轮机的热力模型在Matlab仿真平台上建立。热力模型的输入条件为环境 条件(大气温度to、压力Po、相对湿度RH)、发电机输出功率Ne(作为操作条件）、燃料组分、燃 料低位热值、气路部件健康指数Δ夺^对于新投运机组，么冢=0)。热力模型的计算输出为燃 料流量&、各个部件进出口气路截面处的热力参数(如总压、总温)及转速等。
[0140] 用相对折合参数重新定义压气机和透平的气路健康指数(代表了部件特性线的偏 移），消除由于环境条件(大气压力、温度和相对湿度)变化而导致燃气轮机运行性能变化的 影响。
[0141 ]压气机气路健康指数定义如下：
[0142] SFc ，FC - Gc, cor, rel, deg/Gc, cor, rel
[0143] Λ SFc ,FC (Gc, cor,rel, deg-Gc, cor, rel )/Gc, cor, rel
[0144] SFc.Eff = nc,deg/ric
[0145] Δ SFc.Eff = (nc,deg-nc)/nc
[0146] 其中SFC,FC为压气机流量特性指数;6^_#1,<^为压气机(性能衰退时)相对折合流 量;Ge, , rel为压气机(健康时)相对折合流量;SFe, Ef f为压气机效率特性指数;ne,deg为压气 机(性能衰退时)等熵效率;ne为压气机(健康时)等熵效率。
[0147] 燃烧室气路健康指数定义如下：
[0148] SFB,Eff = riB,degAlB
[0149] Δ SFB.Eff = (riB.deg-%)/%
[0150] 其中SFb, Eff为燃烧室燃烧效率性能指数;nB, ^为燃烧室(性能衰退时)燃烧效率;％ 为燃烧室(健康时)燃烧效率。
[0151 ]透平气路健康指数定义如下：
[0152] SFT,FC = GT,cor,deg/GT,cor
[0153] Δ SFt,FC= (GT,cor,deg-GT,cor)/GT,cor
[0154] SFT,Eff = riT,degAlT
[0155] Δ SFT.Eff = (ητ,deg-ητ) At
[0156] 其中SFt,FCS透平流量性能指数;GT^r.deg为透平(性能衰退时)折合流量;Gt.ci为 透平(健康时)折合流量;SF T, Eff为透平效率性能指数;ητ, deg为透平(性能衰退时)等熵效率； ητ为透平(健康时)等熵效率。
[0157] 在某一运行工况下，采集当前对象燃气轮机稳定运行时的某一时段的气路测量参 数，进行降噪处理后作为待离线诊断的气路测量参数。
[0158] 在对气路测量参数进行调和时，实际的气路部件特性也可能会与理论特性发生一 定偏差，这里将气路部件的流量特性指数和效率特性指数引入作为"虚拟"测量参数(部件 健康时，其值为1 )，一同进行调和，由于未知气路部件性能健康状况，因此其"虚拟"测量参 数都取为1，其不确定度取为1%。基于高斯修正准则的数据调和原理对该工况下的待离线 诊断的气路测量数据及"虚拟"测量参数进行数据调和，检测出发生性能衰退的气路传感器 (即调节量超过置信限值的传感器），如表3所示。
[0159] 表3数据调和结果
[0160]
[0162] 根据检测出存在问题的传感器数目，按照可信的气路测量参数> =(待诊断的部件 气路健康指数数目+存在问题传感器的数目）原则，确定对象燃气轮机部件性能诊断所需的 稳态运行工况点的数目。
[0163] 根据所需的稳态运行工况点数目，通过调整对象燃气轮机操作条件的方式，逐一 采集当前对象燃气轮机在各个稳定运行工况时的某一时段的气路测量参数，进行降噪处理 后作为多个运行工况点的待离线诊断的气路测量数据集。
[0164] 分别设置已建立的燃气轮机热力模型的环境输入条件(大气压力、温度和相对湿 度)和操作输入条件与采样时的对象燃气轮机的各个运行工况一致，消除由于环境条件(大 气压力、温度和相对湿度)和操作条件变化而导致燃气轮机运行性能变化的影响。
[0165] 以待离线诊断的多个运行工况点的气路测量数据ζ与热力模型计算的气路参数数 据^,之间的均方根误差为目标函数，通过粒子群优化算法计算得到当前的各个部件(压气 机、透平和燃烧室）的气路健康指数6尸用以评估对象燃气轮机实际的性能健康状况，如 图2所示。
[0167] 式中丨--由燃机热力模型在第i个运行工况点计算得到的气路测量参数向量；
[0168] --在第i个运行工况点实际的气路测量参数向量；
[0169] Μ一一某一运行工况点的气路测量参数的总数目。
[0170] 优化目标函数定义如式：
[0172]式中Fitness是优化目标函数，当Fitness随着粒子群优化算法(PS0)迭代寻优计算 过程逐渐地趋近于〇时，预测的气路测量参数[H,…,九]与实测气路测量参数 相匹配，此时输出最优的气路部件健康指数Δ!?和在不同运行工况点的传感器偏差
[0173] 粒子群优化算法相关参数的选取如表4所示，这里进化代数为200,种群规模为 1〇〇,用于搜索最优的部件健康参数Δ尿（如表2所示)和传感器偏差g。
[0174] 表4粒子群优化算法(PS0)相关参数的选取
[0176] 通过本发明诊断方法上述的主要诊断步骤后，可以得到最终的诊断结果如图4所 不。
[0177] 图3中横坐标的标志符12表示P5，标志符13表示t5，标志符14表示P6，标志符15表示 t6。从图4可知，当气路部件发生性能衰退且存在传感器性能衰退时，通过这种抗传感器故 障的强鲁棒性燃气轮机气路部件性能诊断方法(Robust PS0-GPA)能够成功地识别、隔离性 能衰退的部件及部件，并准确地量化性能衰退程度，避免误导性的诊断结果产生。其中诊断 案例相应的这种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方法迭代计算过程 如图5所示。
[0178] 基于以上诊断案例的分析可知，当气路部件发生性能衰退且存在传感器偏差时， 通过本发明的诊断方法能够成功地识别、隔离性能衰退的部件，并准确地量化性能衰退程 度。
【主权项】
1. 一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方法，其特征是： (1) 基于对象燃气轮机新投运或健康时的气路测量参数，建立燃气轮机非线性热力模 型，压气机和透平都用相似折合参数形式表示； (2) 用相似折合参数重新定义压气机和透平的气路健康指数，消除由于环境条件变化 而导致燃气轮机运行性能变化的影响； (3) 采集当前对象燃气轮机在一运行工况下稳定运行时的一个时段的气路测量参数， 进行降噪处理后作为待离线诊断的气路测量参数； (4) 对待离线诊断的气路测量参数进行数据调和，检测出发生性能衰退的气路传感器； (5) 根据检测出存在问题的传感器数目，按照可信的气路测量参数含（待诊断的部件气 路健康指数数目+存在问题传感器的数目）原则，确定对象燃气轮机部件性能诊断所需的稳 态运行工况点的数目； (6) 根据所需的稳态运行工况点数目，通过调整对象燃气轮机操作条件的方式，逐一采 集当前对象燃气轮机在各个稳定运行工况时的一个时段的气路测量参数，进行降噪处理后 作为多个运行工况点的待离线诊断的气路测量数据集； (7) 分别设置已建立的燃气轮机热力模型的环境输入条件和操作输入条件与采样时的 对象燃气轮机的各个运行工况一致，消除由于环境条件和操作条件变化而导致燃气轮机运 行性能变化的影响； (8) W待离线诊断的多个运行工况点的气路测量数据与热力模型的气路参数数据之间 的均方根误差为目标函数，通过粒子群优化算法计算得到当前的各个部件的气路健康指 数，用W评估对象燃气轮机实际的性能健康状况。2. 根据权利要求1所述的一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方 法，其特征是:步骤(1)中建立燃气轮机非线性热力模型的具体步骤如下： (a)利用部件相对折合参数，建立燃气轮机部件级热力模型，其中压气机和透平特性线 数据整理成通用的相对折合参数形式： 压气机特性线整理成通用的相似折合参数形式如下：其4为相对折合转速，η为实际转速，?；;为压气机进口滞止溫 度，Rg为流经压气机工质的气体常数，下角标0表示设计点；为相对折合流量，Ge为实际压气机进口流量，钱为压 气机进口滞止压力为相对压比，化为实际压气机压比，化,rel =化/化日为相对等赌 效率，化为实际压气机等赌效率； 透平特性线整理成通用的相似折合参数形式如下：式中：为相对折合转速，η为实际转速，7：为透平进口滞止溫 度，Rg为流经透平工质的气体常数％相对折合流量，Gt 为实际透平进口流量妇相对压比，化κ?=ητΑ?τ日为相对等赌效率，ητ为实际透平 等赌效率，下角标0表示设计点； (b)根据采集的对象燃气轮机新投运或健康时的气路测量参数，逐步修正各个部件的 特性线数据，使所建的燃气轮机热力模型的计算值与气路实测参数相匹配，从而消除热力 模型计算误差给诊断结果带来的负影响。3. 根据权利要求2所述的一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方 法，其特征是:步骤(2)用相似折合参数重新定义压气机和透平的气路健康指数的具体步骤 如下： (a) 燃气轮机总体性能健康状况通过由各部件的气路健康指数来表示，用相对折合参 数重新定义压气机和透平的气路健康指数W消除由于环境条件变化而给诊断结果带来的 负影响。 (b) 压气机气路健康指数定义如下：其中SFg,FG为压气机流量特性指数;GG,wr,rel,deg为压气机性能衰退时相对折合流量； Gc, cor, re功压气机健康时相对折合流量;SFG,Eff为压气机效率特性指数;riG, deg为压气机性能 衰退时等赌效率;化为压气机健康时等赌效率； 燃烧室气路健康指数定义如下：其中SFb, Eff为燃烧室燃烧效率性能指数;化deg为燃烧室性能衰退时燃烧效率;％为燃烧 室健康时燃烧效率。 透平气路健康指数定义如下：其中SFt, FC为透平流量性能指数;Gt, cor, deg为透平性能衰退时折合流量;Gt, cor为透平健 康时折合流量;SFT,Eff为透平效率性能指数;ητ,deg为透平性能衰退时等赌效率;ητ为透平健 康时等赌效率。4. 根据权利要求3所述的一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方 法，其特征是:步骤(4)基于高斯修正准则的数据调和原理对该工况下的待离线诊断的气路 测量数据进行数据调和，检测出发生性能衰退的气路传感器的具体步骤如下： (a) 在一个运行工况下，采集当前对象燃气轮机稳定运行时的一个时段的气路测量参 数，进行降噪处理后作为待离线诊断的气路测量参数； (b) 在对气路测量参数进行调和时，将气路部件的流量特性指数和效率特性指数引入 作为"虚拟"测量参数，一同进行调和，其"虚拟"测量参数都取为1，其不确定度取为1 % ; (C)基于高斯修正准则的数据调和原理对该工况下的待离线诊断的气路测量数据及 "虚拟"测量参数进行数据调和，检测出发生性能衰退的气路传感器。5.根据权利要求4所述的一种抗传感器故障的强鲁棒性燃气轮机气路部件性能诊断方 法，其特征是:步骤(8) W待离线诊断的多个运行工况点的气路测量数据与热力模型计算的 气路参数数据之间的偏差作为残差，计算得到当前的各个部件的气路健康指数，用W评估 对象燃气轮机实际的性能健康状况的具体步骤如下： (a) 分别设置已建立的燃气轮机热力模型的环境输入条件和操作输入条件与采样时的 对象燃气轮机的各个运行工况一致，消除由于环境条件和操作条件变化而导致燃气轮机运 行性能变化的影响； (b) W待离线诊断的多个运行工况点的气路测量数据i与热力模型计算的气路参数数 据^之间的均方根误差为目标函数Fitness,通过粒子群优化算法计算得到当前的各个部件 的气路健康指鑽Δ赤，用W评估对象燃气轮机实际的性能健康状况，式中Z,为由燃机热力模型在第i个运行工况点计算得到的气路测量参数向量，i；为在第i 个运行工况点实际的气路测量参数向量，Μ为某一运行工况点的气路测量参数的总数目； 优化目标函数定义如式：式中Fitness是优化目标函数，当Fitness随着粒子群优化算法(PSO)迭代寻优计算过程 逐渐地趋近于0时，预测的气路测量参数口ι，-V..，ミJ与实测气路测量参数[i；，i，.。.，。相匹 配，此时输出最优的气路部件健康指数Δ涼和在不同运行工况点的传感器偏差 -一. 一， [知，化··，扳规]〇
【文档编号】G01M15/14GK106092594SQ201610362577
【公开日】2016年11月9日
【申请日】2016年5月26日
【发明人】李淑英, 应雨龙, 曹云鹏
【申请人】哈尔滨工程大学