CN105067275A - 基于排温偏离指数的燃气轮机燃烧系统在线监测诊断方法 - Google Patents
基于排温偏离指数的燃气轮机燃烧系统在线监测诊断方法 Download PDFInfo
- Publication number
- CN105067275A CN105067275A CN201510443528.0A CN201510443528A CN105067275A CN 105067275 A CN105067275 A CN 105067275A CN 201510443528 A CN201510443528 A CN 201510443528A CN 105067275 A CN105067275 A CN 105067275A
- Authority
- CN
- China
- Prior art keywords
- thermopair
- gas turbine
- row
- temperature
- data
- 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.)
- Pending
Links
Landscapes
- Control Of Turbines (AREA)
Abstract
基于排温偏离指数的燃气轮机燃烧系统在线监测诊断方法,涉及燃气轮机燃烧系统在线监测诊断方法。为了解决现有的燃气轮机燃烧系统在线监测方法存在严重“事后”诊断现象的问题。本发明首先获取燃气轮机正常运行tm时间内的排温数据T,计算tj时刻n个热电偶的平均温度;并将tj时刻不同热电偶的排温数据与分别组成一一对应的关系;求得第i个热电偶排温数据和平均温度的关系;针对第i个热电偶,将βi和排温数据和平均温度数据带入到中,求得第i个热电偶tj时刻的,统计出,设定的阈值为[μi-3σi,μi+3σi];当燃气轮机运行时,求解某个时间段每个热电偶对应排温数据的;若超出了对应的阈值[μi-3σi,μi+3σi],判断该热电偶对应的燃烧室内发生故障;反之,正常。本发明适用于燃气轮机燃烧系统的在线监测诊断。
Description
技术领域
本发明涉及燃气轮机燃烧系统在线监测诊断方法。
背景技术
随着技术的进步,燃气轮机以其结构紧凑、启动快、运行平稳、热效率高等优势已成为21世纪核心动力装备之一。异常检测与故障诊断技术是保障燃气轮机安全运行的重要手段。燃烧室是燃气轮机的重要组成部分,其工作条件恶劣,故障频发,且一旦损坏后果十分严重,因此加强燃烧室的异常检测和故障诊断有着十分重要的意义。
然而,由于燃烧室工作时,温度极高,现行的测温手段无法在这种条件下长期工作。而在实际运行时,当燃烧筒出现异常的时候,排温的结果也会出现异常,所以用户通常是通过排温异常来判断燃烧筒的工作情况是否出现异常。可以通过均匀布置在机组透平排气通道中的热电偶测量燃气轮机的排温,来间接反映燃烧室内的燃烧情况。
《Real-timeon-lineperformancediagnosticsofheavy-dutyindustrialgasturbines》中是利用形状因子来实现排温的在线监测,其形状因子的定义为t时刻下,热电偶测得的最高温与平均温度的比值,由于温度的突变会影响形状因子,所以便可以用一段时间内的形状因子的变化来推测燃机燃烧室是否发生了异常,如当形状因子突然高于正常范围,则说明机组的排温出现了异常。
该方法仅仅利用了各个温度测点之间的最高值来判断排温是否发生异常,不能全面的描述燃机燃烧室工作的情况及异常演变的过程,当某些异常发生时,会导致一些热电偶的测得的温度出现问题,若这些排温不是最高温度,则该方法不能检测出系统异常。
GE公司开发的MARKVI燃烧监测系统定义S为排气温度的允许排温分散度,认为S是燃气轮机出口的平均排气温度T4 *、压气机出口温度T2 *的函数,具体函数是个经验公式:
在该公式里,温度均是以℉为计量单位的。公式右端的(100)带有括号,表示变工况条件下才加入该项。
此外,MARKVI燃烧监测系统还定义:S1为排气温度热电偶的最高读数与最低读数之间的差;S2为排气温度热电偶的最高读数与第2个低读数之间的差;S3为排气温度热电偶的最高读数与第3个低读数之间的差。
基于上述的公式和定义,MARKVI燃烧监测保护系统的判别原理见图1。图1中,K1,K2,K3是三个依据经验定义的参数。典型情况下:K1=1.0;K2=5.0;K3=0.8;
但是,在实际应用中发现,当检测系统发出报警时燃烧系统已经损坏较严重,存在严重的“事后”诊断现象。
发明内容
本发明为了解决现有的燃气轮机燃烧系统在线监测方法存在严重“事后”诊断现象的问题,进而提出了基于排温偏离指数的燃气轮机燃烧系统在线监测诊断方法。
基于排温偏离指数的燃气轮机燃烧系统在线监测诊断方法,包括下述步骤:
步骤1、设燃气轮机在排气端均匀分布着n个热电偶(测点),获取燃气轮机正常运行时tm时间内的排温数据每行数据为第i个热电偶tm时间内各个时刻的排温数据,每列为某个时刻n个热电偶的排温数据;其中,表示第i个热电偶tj时刻测得的温度,i=1,2,...,n,tj=1,2,...,tm;
计算tj时刻n个热电偶的平均温度并将tj时刻不同热电偶的排温数据与分别组成一一对应的关系:
步骤2、提取第一个热电偶不同时刻对应关系数据:并根据这些数据进行拟合,得到如公式(1)的排温与平均温度关系;
用相同方法求得第i个热电偶(每个热电偶)的排温与平均温度关系,如公式(2)所示;
将αi定义为第i个热电偶排温偏离指数;
根据每个热电偶的拟合公式,得到第i个热电偶(每个热电偶)对应的βi;
步骤3、针对第i个热电偶,在不同时刻有
为第i个热电偶tj时刻(不同时刻)排温偏离指数;
将步骤2得到的每个热电偶对应的βi的值带入到公式(3)中,并将对应热电偶的tj时刻排温数据和平均温度数据带入到公式(3)中,得到第i个热电偶tj时刻(每个热电偶不同时刻)的 构成一个时间序列;
当在燃气轮机正常运转情况下,的值保持恒定,但由于噪声的存在,值会在一定范围内波动,由该波动范围可确定正常运行时的阈值;若燃气轮机发生异常,用该段时间燃气轮机的运行数据求出的会超出正常运行数据所确定的阈值范围,若燃气轮机运行正常,则求出的在该阈值范围内;
在燃气轮机正常运行时,针对第i个热电偶,统计的正态分布情况:其中μi为均数、σi为标准差;则的阈值为[μi-3σi,μi+3σi],如图3所示;
步骤4、当燃气轮机运行时,采集某一时间段内热电偶对应的排温数据,并求解该时间段每个热电偶对应排温数据的若超出了对应的阈值[μi-3σi,μi+3σi],则判断该热电偶对应的燃烧室内发生故障,若未超出,则燃气轮机运行正常,为如图4所示,。
本方法可以充分利用排温各个测点的数据信息,准确检测出异常演变过程,在故障发生的早期就能及时发现,减小经济损失,提高安全性。
附图说明
图1为MARKVI燃烧监测保护系统的判别原理图;
图2为利用最小二乘法进行线性回归拟合的拟合图;
图3为统计出的正态分布图;
图4为基于排温偏离指数的燃气轮机燃烧系统在线监测诊断方法检测示意图。
具体实施方式
具体实施方式一:
1、基于排温偏离指数的燃气轮机燃烧系统在线监测诊断方法,其特征在于它包括下述步骤:
步骤1、设燃气轮机在排气端均匀分布着n个热电偶,获取燃气轮机正常运行时tm时间内的排温数据每行数据为第i个热电偶tm时间内各个时刻的排温数据,每列为某个时刻n个热电偶的排温数据;其中,表示第i个热电偶tj时刻测得的温度,i=1,2,...,n,tj=1,2,...,tm;
计算tj时刻n个热电偶的平均温度并将tj时刻不同热电偶的排温数据与分别组成一一对应的关系:
步骤2、提取第一个热电偶不同时刻对应关系数据:并根据这些数据进行拟合,得到如公式(1)的排温与平均温度关系;
用相同方法求得第i个热电偶(每个热电偶)的排温与平均温度关系,如公式(2)所示;
将αi定义为第i个热电偶排温偏离指数;
根据每个热电偶的拟合公式,得到第i个热电偶(每个热电偶)对应的βi;
步骤3、针对第i个热电偶,在不同时刻有
为第i个热电偶tj时刻(不同时刻)排温偏离指数;
将步骤2得到的每个热电偶对应的βi的值带入到公式(3)中,并将对应热电偶的tj时刻排温数据和平均温度数据带入到公式(3)中,得到第i个热电偶tj时刻(每个热电偶不同时刻)的 构成一个时间序列;
当在燃气轮机正常运转情况下,的值保持恒定,但由于噪声的存在,值会在一定范围内波动,由该波动范围可确定正常运行时的阈值;若燃气轮机发生异常,用该段时间燃气轮机的运行数据求出的会超出正常运行数据所确定的阈值范围,若燃气轮机运行正常,则求出的在该阈值范围内;
在燃气轮机正常运行时,针对第i个热电偶,统计出的正态分布情况:其中μi为均数、σi为标准差;则的阈值为[μi-3σi,μi+3σi];如图3所示。
步骤4、当燃气轮机运行时,采集某一时间段内热电偶对应的排温数据,并求解该时间段每个热电偶对应排温数据的若超出了对应的阈值[μi-3σi,μi+3σi],则判断该热电偶对应的燃烧室内发生故障,若未超出,则燃气轮机运行正常,如图4所示。
具体实施方式二:本实施方式所述的步骤2所述进行拟合的过程采用最小二乘法进行拟合。利用最小二乘法进行线性回归拟合的拟合图如图2所示。
其它步骤和参数与具体实施方式一相同。
Claims (2)
1.基于排温偏离指数的燃气轮机燃烧系统在线监测诊断方法,其特征在于它包括下述步骤:
步骤1、设燃气轮机在排气端均匀分布着n个热电偶,获取燃气轮机正常运行时tm时间内的排温数据其中,表示第i个热电偶tj时刻测得的温度,i=1,2,…,n,tj=1,2,…,tm;
计算tj时刻n个热电偶的平均温度并将tj时刻不同热电偶的排温数据与分别组成一一对应的关系:
步骤2、提取第一个热电偶不同时刻对应关系数据:并根据这些数据进行拟合,得到如公式(1)的排温与平均温度关系;
用相同方法求得第i个热电偶的排温与平均温度关系,如公式(2)所示;
将αi定义为第i个热电偶排温偏离指数;
根据每个热电偶的拟合公式,得到第i个热电偶对应的βi;
步骤3、针对第i个热电偶,在不同时刻有
为第i个热电偶tj时刻排温偏离指数;
将步骤2得到的每个热电偶对应的βi的值带入到公式(3)中,并将对应热电偶的tj时刻排温数据和平均温度数据带入到公式(3)中,得到第i个热电偶tj时刻的构成一个时间序列;
在燃气轮机正常运行时,针对第i个热电偶,统计的正态分布情况:其中μi为均数、σi为标准差;则的阈值为[μi-3σi,μi+3σi];
步骤4、当燃气轮机运行时,采集某一时间段内热电偶对应的排温数据,并求解该时间段每个热电偶对应排温数据的若超出了对应的阈值[μi-3σi,μi+3σi],则判断该热电偶对应的燃烧室内发生故障,若未超出,则燃气轮机运行正常。
2.根据权利要求1所述的基于排温偏离指数的燃气轮机燃烧系统在线监测诊断方法,其特征在于步骤2所述进行拟合的过程采用最小二乘法进行拟合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510443528.0A CN105067275A (zh) | 2015-07-24 | 2015-07-24 | 基于排温偏离指数的燃气轮机燃烧系统在线监测诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510443528.0A CN105067275A (zh) | 2015-07-24 | 2015-07-24 | 基于排温偏离指数的燃气轮机燃烧系统在线监测诊断方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105067275A true CN105067275A (zh) | 2015-11-18 |
Family
ID=54496699
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510443528.0A Pending CN105067275A (zh) | 2015-07-24 | 2015-07-24 | 基于排温偏离指数的燃气轮机燃烧系统在线监测诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105067275A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105510045A (zh) * | 2015-11-27 | 2016-04-20 | 哈尔滨工业大学 | 基于系数矩阵的燃气轮机燃烧系统在线监测方法 |
CN108827643A (zh) * | 2018-06-21 | 2018-11-16 | 哈尔滨工业大学 | 一种考虑排温温度场旋转的燃气轮机高温部件故障预警方法 |
CN109556876A (zh) * | 2018-11-07 | 2019-04-02 | 国网浙江省电力有限公司电力科学研究院 | 一种区分燃气轮机燃烧故障和热通道设备故障的诊断方法 |
CN110672332A (zh) * | 2019-09-10 | 2020-01-10 | 上海电力大学 | 一种基于sarima模型的燃气轮机故障预警系统 |
CN111855221A (zh) * | 2020-07-28 | 2020-10-30 | 青岛中科国晟动力科技有限公司 | 一种燃气轮机燃烧状态监测方法及系统 |
CN111879522A (zh) * | 2020-07-24 | 2020-11-03 | 山东大学 | 基于时序概率的汽轮机运行监测和故障判别方法及系统 |
CN112460634A (zh) * | 2020-11-23 | 2021-03-09 | 西安热工研究院有限公司 | 一种确定燃气轮机故障燃烧室的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1267823A (zh) * | 1999-03-19 | 2000-09-27 | 中国船舶工业总公司第七研究院第七一一研究所 | 气体发动机气缸不正常燃烧现象的检测方法及报警装置 |
CN2625876Y (zh) * | 2002-12-30 | 2004-07-14 | 沈阳黎明航空发动机(集团)有限责任公司 | 发动机燃烧室温度场旋转测试装置 |
JP2005249420A (ja) * | 2004-03-01 | 2005-09-15 | Mazda Motor Corp | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム |
CN101403654A (zh) * | 2008-11-06 | 2009-04-08 | 西安交通大学 | 一种用于燃气轮机涡轮叶片的双工质冷却实验系统 |
CN103195583A (zh) * | 2013-04-07 | 2013-07-10 | 北京华清燃气轮机与煤气化联合循环工程技术有限公司 | 一种采用排气温度分散度监测保护燃气轮机燃烧的方法 |
CN103674555A (zh) * | 2012-09-14 | 2014-03-26 | 通用电气公司 | 用于燃气涡轮排气中的物质曲线测量的方法和系统 |
-
2015
- 2015-07-24 CN CN201510443528.0A patent/CN105067275A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1267823A (zh) * | 1999-03-19 | 2000-09-27 | 中国船舶工业总公司第七研究院第七一一研究所 | 气体发动机气缸不正常燃烧现象的检测方法及报警装置 |
CN2625876Y (zh) * | 2002-12-30 | 2004-07-14 | 沈阳黎明航空发动机(集团)有限责任公司 | 发动机燃烧室温度场旋转测试装置 |
JP2005249420A (ja) * | 2004-03-01 | 2005-09-15 | Mazda Motor Corp | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム |
CN101403654A (zh) * | 2008-11-06 | 2009-04-08 | 西安交通大学 | 一种用于燃气轮机涡轮叶片的双工质冷却实验系统 |
CN103674555A (zh) * | 2012-09-14 | 2014-03-26 | 通用电气公司 | 用于燃气涡轮排气中的物质曲线测量的方法和系统 |
CN103195583A (zh) * | 2013-04-07 | 2013-07-10 | 北京华清燃气轮机与煤气化联合循环工程技术有限公司 | 一种采用排气温度分散度监测保护燃气轮机燃烧的方法 |
Non-Patent Citations (2)
Title |
---|
俞立凡等: "9FA燃机排气温度场偏转规律研究", 《电力建设》 * |
刘娇: "燃气轮机排温异常检测研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105510045A (zh) * | 2015-11-27 | 2016-04-20 | 哈尔滨工业大学 | 基于系数矩阵的燃气轮机燃烧系统在线监测方法 |
CN105510045B (zh) * | 2015-11-27 | 2017-12-01 | 哈尔滨工业大学 | 基于系数矩阵的燃气轮机燃烧系统在线监测方法 |
CN108827643A (zh) * | 2018-06-21 | 2018-11-16 | 哈尔滨工业大学 | 一种考虑排温温度场旋转的燃气轮机高温部件故障预警方法 |
CN108827643B (zh) * | 2018-06-21 | 2020-04-07 | 哈尔滨工业大学 | 考虑排温温度场旋转的燃气轮机高温部件故障预警方法 |
CN109556876A (zh) * | 2018-11-07 | 2019-04-02 | 国网浙江省电力有限公司电力科学研究院 | 一种区分燃气轮机燃烧故障和热通道设备故障的诊断方法 |
CN109556876B (zh) * | 2018-11-07 | 2020-09-04 | 国网浙江省电力有限公司电力科学研究院 | 一种区分燃气轮机燃烧故障和热通道设备故障的诊断方法 |
CN110672332A (zh) * | 2019-09-10 | 2020-01-10 | 上海电力大学 | 一种基于sarima模型的燃气轮机故障预警系统 |
CN111879522A (zh) * | 2020-07-24 | 2020-11-03 | 山东大学 | 基于时序概率的汽轮机运行监测和故障判别方法及系统 |
CN111855221A (zh) * | 2020-07-28 | 2020-10-30 | 青岛中科国晟动力科技有限公司 | 一种燃气轮机燃烧状态监测方法及系统 |
CN112460634A (zh) * | 2020-11-23 | 2021-03-09 | 西安热工研究院有限公司 | 一种确定燃气轮机故障燃烧室的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105067275A (zh) | 基于排温偏离指数的燃气轮机燃烧系统在线监测诊断方法 | |
CN103195583B (zh) | 一种采用排气温度分散度监测保护燃气轮机燃烧的方法 | |
CN105241669A (zh) | 基于比较编码的燃气轮机燃烧系统在线监测方法 | |
CN111365083B (zh) | 一种基于热力参数的汽轮机通流部分故障诊断方法 | |
CA2865200A1 (en) | Method and system for real time dry low nitrogen oxide (dln) and diffusion combustion monitoring | |
CN112257943A (zh) | 一种基于mset和偏离度的电厂引风机故障预警方法 | |
CN102682862B (zh) | 压水堆核电站堆芯热电偶有效性状态在线诊断方法 | |
US9651457B2 (en) | Method for detecting deterioration in a turbomachine by monitoring the performance of said turbomachine | |
JP2011106467A (ja) | ガスタービンの性能診断システムにおける表示画面の表示方法 | |
CN103675355A (zh) | 风速仪监测方法和系统 | |
CN106641045A (zh) | 联合循环电厂同步离合器的安全监测系统和方法 | |
CN109298700B (zh) | 一种实时判断火电机组运行参数异常变化的方法及系统 | |
CN111929579A (zh) | 发电机在线故障诊断方法、设备及计算机设备 | |
CN109556876B (zh) | 一种区分燃气轮机燃烧故障和热通道设备故障的诊断方法 | |
CN105114977A (zh) | 一种基于排温测点相关性的燃机燃烧系统在线监测方法 | |
CN105134386A (zh) | 基于测点加权值的燃气轮机燃烧系统在线监测方法 | |
JP2016009352A (ja) | プラントの異常予兆診断装置及び方法 | |
CN111312414A (zh) | 核电厂疲劳监测和寿命评估系统的疲劳评价方法 | |
CN212618234U (zh) | 一种锅炉过热器寿命的实时监测系统 | |
CN106768000A (zh) | 一种风力发电机组变流器水冷系统压力异常检测方法 | |
CN116484267B (zh) | 变压器故障特征提取、判定方法、计算机设备和存储介质 | |
CN105510045A (zh) | 基于系数矩阵的燃气轮机燃烧系统在线监测方法 | |
CN110333265B (zh) | 一种预测机车发动机排气波纹管剩余寿命的方法及系统 | |
CN112983570A (zh) | 基于相关性的汽轮机轴承温度高跳机征兆捕捉方法及装置 | |
CN114199505B (zh) | 一种基于相关性分析的发电机定子线棒流通性评价方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20151118 |
|
WD01 | Invention patent application deemed withdrawn after publication |