CN102877945A - 基于解析法的燃气轮机及其联合循环的变工况分析方法 - Google Patents
基于解析法的燃气轮机及其联合循环的变工况分析方法 Download PDFInfo
- Publication number
- CN102877945A CN102877945A CN2012101971352A CN201210197135A CN102877945A CN 102877945 A CN102877945 A CN 102877945A CN 2012101971352 A CN2012101971352 A CN 2012101971352A CN 201210197135 A CN201210197135 A CN 201210197135A CN 102877945 A CN102877945 A CN 102877945A
- Authority
- CN
- China
- Prior art keywords
- gas turbine
- igv
- compressor
- condition
- ambient temperature
- 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
Images
Landscapes
- Control Of Turbines (AREA)
Abstract
一种基于解析法的燃气轮机及其联合循环的变工况分析方法,以此建立变工况模型,根据模型运算的结果对燃气轮机的环境温度特性和负荷特性进行研究。获得带有IGV温度控制的燃气轮机恒速带负荷运行时的稳态变工况特性。本发明提供的燃气轮机及其联合循环的变工况分析方法,可建立起变工况模型,获得带有IGV温度控制的燃气轮机恒速带负荷运行时的稳态变工况特性,根据模型运算的结果可对燃气轮机的环境温度特性和负荷特性进行研究。
Description
技术领域
本发明涉及一种基于解析法的典型燃气轮机及其联合循环的变工况分析方法。
背景技术
由于外界需求负荷、环境温度等诸因素变化,燃气轮机出力、效率或热耗率等性能指标亦随之改变,机组常常处于变工况运行状态。考察、比较不同燃气轮机机组的负荷特性及压气机积垢特性,就需要将燃气轮机的相关性能指标折算到某种基准条件。
燃气轮机属于定容式动力机械,环境温度对其功率和效率的影响有相当大的影响。定性分析,这是由于:①随着环境温度升高,空气密度减小,致使吸入压气机的空气质量流量减少,机组的做功能力随之变小;②压气机的耗功量随吸入空气的热力学温度成正比,即环境温度升高时,燃气轮机的净出力减小;③当环境温度升高时,即使机组的转速和燃气初温保持恒定,压气机的压缩比将有所下降,这将导致燃气透平做功量的减少,而燃气透平的排气温度却有所提高。
对于带有可调导叶(IGV)温度控制的燃气轮机,其环境温度特性的研究远比定性分析复杂,压气机设置IGV的燃气轮机,当机组启停或调整负荷时,通过调节IGV叶片角度,限制进入压气机的空气流量,从而达到保护机组安全运行和提高运行效率的目的。由于IGV参与燃气轮机的自动控制,使得燃气轮机各部件的参数及环境温度特性的研究变得复杂。
因为外界负荷、环境温度等因素的影响,燃气轮机的参数变化频繁,很难通过实际采集的现场数据得到一定负荷情况下的环境温度特性。这就需要尽量结合现场数据,通过理论分析、计算的方法,建立燃气轮机的变工况模型,获得燃气轮机的环境温度特性、进排气阻力特性和负荷特性。
发明内容
本发明所要解决的技术问题,就是提供一种燃气轮机及其联合循环的变工况分析方法。以此建立变工况模型,根据模型运算的结果对燃气轮机的环境温度特性和负荷特性进行研究。获得带有IGV温度控制的燃气轮机恒速带负荷运行时的稳态变工况特性。
解决上述技术问题,本发明采用的技术方案为:
一种基于解析法的燃气轮机及其联合循环的变工况分析方法,包括以下步骤:
S1 进行压气机进气歧管差压数据的合理性分析,剔除不合理数据;
S2 根据步骤S1得到的数据样本,结合分散控制系统(DCS)数据进行燃气轮机运行性能测算;
S3 建立燃气轮机的变工况模型,包括以下步骤;
S3.1建立燃气轮机的IGV温控规律模型;
S3.2建立压气机压比π c ~ ( IGV , ta )模型;
S3.3建立压气机流量Gc ~ ( IGV , ta )模型;
S3.4建立压气机进气阻力模型;
S3.5建立IGV温控规律下的压气机效率-流量特性曲线;
S3.6建立燃气透平的膨胀比模型;
S3.7建立燃气透平的流量模型;
S3.8建立燃气透平等熵效率模型;
S4 进行燃气轮机整体变工况性能模拟;
S5 根据步骤S4的结果画出燃气轮机整体变工况曲线图,研究燃气轮机环境温度特性:包括功率-环境温度特性、效率-环境温度特性、排气参数-环境温度特性、压气机的压比-环境温度特性;
S6 根据步骤S4的结果研究进排气阻力对燃气轮机性能的影响;
S7 根据步骤S4的结果研究燃气轮机的负荷特性,包括:
S7.1根据步骤S5的结果计算基准条件下燃气轮机性能的影响因子:基准功率-环境温度折算因子、基准功率-环境压力折算因子、基准功率-相对湿度折算因子、基准效率-环境温度折算因子;
S7.2根据步骤S7.1的结果将燃气轮机功率及效率修正到基准环境条件,绘制燃气轮机折算效率-折算功率曲线,并对曲线数据进行拟合得到燃气轮机折算效率—折算功率的数学模型。
所述的步骤S1中的进行压气机进气歧管差压数据的合理性分析为:
采用压气机进气歧管静压法,计算压气机的空气流量,即按式(2-1)估算:
式中:Gc 为压气机入口空气流量,kg/s;Pin 为压气机进气歧管处的负压;∆Pin 为过滤器阻力;ρ a 为空气密度;A为进气歧管负压测点所在的过流断面面积;Ф为流量系数;Ф *A可成为有效过流断面面积;
由于压气机进气歧管静压测量值可能有较大误差,导致数据波动,故需要剔除明显错误的数据,本发明提出一种用于鉴别进气歧管差压值合理性的判别式形式:
(2-2)
式中:Dp是压气机进气歧管差压,;上标(i)表示第i次数据样本;a 、 b 、 c为常数,与数据样本有关,由曲线的拟合得到;IGV为可调导叶开度;
所述的步骤S1中的鉴别歧管静压数据合理的步骤为:
③重复步骤②,直到数据样本满足要求。
所述的步骤S2为:
根据S1得到的数据样本结合DCS数据进行燃气轮机运行性能测算,根据测算结果进行建模样本的选择。
所述的步骤S3包括以下步骤:
S3.1建立燃气轮机的IGV温控规律模型
综合分析燃气轮机排气温度T4 - IGV规律图和燃气轮机透平计算进气温度T3 - IGV规律图,得到燃气轮机的IGV温控规律模型:在燃气轮机基本负荷向部分负荷过渡的过程中,先关小IGV,减小压气机进气流量(透平排气温度T4 不变,而进气温度T3 略减小),再逐渐降低T3 。
燃气轮机的IGV温控规律模型可采用如下模型:
式中,常数a 、 b 、 c 、 d与机组设定的温控规律有关,常数e为表示IGV开度的百分数,由实际的燃气轮机透平进气温度T3 - IGV关系图拟合得到;IGV表示可调导叶开度;
S3.2建立压气机压比π c ~ ( IGV , ta )模型
分析压气机压比随IGV开度的变化图,采用如下形式建立模型:
式中:常数a 、 b 、 c 、 d 、 e 、 f与机组设计性能及运行状况有关,根据实际的压比π c ~ ( IGV , ta )曲线得到,环境压力认为定值,忽略进气道阻力及叶片积垢的影响;
S3.3建立压气机流量Gc ~ ( IGV , ta )模型
采用歧管静压法分析压气机空气流量测算值随IGV开度的变化图,用如下形式建立模型:
式中:常数a 、 b 、 c 、 d 、 e 、 f与机组设计性能及运行状况有关,根据实际的流量Gc ~ ( IGV , ta )曲线得到,环境压力认为定值,忽略进气道阻力及叶片积垢的影响;
式(2-4)和式(2-5)采用了相同的数学模型形式,其中的常数a 、 b 、 c 、 d 、 e 、 f由两个式子各自对应的曲线进行拟合得到;
S3.4建立压气机进气阻力模型
压气机进气阻力主要考虑两级进气过滤器的流动阻力,采用如下形式建立模型:
式中:ξ为管件的阻力系数;c为流速;ta 为环境温度,℃;Pa为环境压力,Pa;Gc 为压气机空气质量流量,kg/s;常数a与机组设计性能及运行状况有关,由实际运行数据拟合得到;
S3.5建立IGV温控规律下的压气机效率-流量特性曲线
一种轴流式压气机的典型效率特性解析式为:
式中:π表示压比;η表示效率;“ˉ”表示与设计值的比值;C4 为常数,可取0.4。
IGV温控下的燃气轮机,从实际运行数据仅能得到其特性点簇。
将压气机的折合流量和折合转速设计值分别代入式(2-7),可将式(2-7)
整理成(2-8)如下:
式中:常数a1 、 b1 与压气机折合流量和折合转速设计值有关。
基于压气机基元叶栅法,可以估算压气机的变工况特性;根据理论分析的结果,轴流式压气机的等熵效率具有如下规律:
IGV温控规律下、恒速运行时,满足上述条件的一种多元函数形式:
(2-9)
式中:常数a 、 b 、 c 、 d 、 e与机组运行状况有关,由实际运行数据的特性点簇得到,常数a1 、 b1 与压气机折合流量和折合转速设计值有关。
S3.6建立燃气透平的膨胀比模型
分析燃气透平膨胀比π t 与压气机压比π c 关系图,采用如下形式建立模型:
式中:常数a 、 b与机组设计性能及运行状态有关。
S3.7建立燃气透平的流量模型
燃气透平的流量采用修正的Flugel公式,其形式:
其中:ψ为流量修正系数,下标0表示ISO条件下的设计值,3表示透平进气,4表示透平排气。
S3.8建立燃气透平等熵效率模型。
一种燃气透平效率-燃气流量典型特性表达式为:
将燃气透平的折合流量和折合转速设计值分别代入式(1-9),可将式(2-12)
整理成(2-13)如下:
式中:常数a 、 b与燃气透平的折合流量和折合转速设计值有关;η t0 为设计效率。
所述的步骤S4包括以下步骤:
S4.1压气机性能模拟:
S4.1.1给定环境温度ta 、环境压力Pa 、压气机IGV开度、燃料特性、燃气轮机设计参数;
S4.1.2按式(2-5)的模型计算压气机空气流量;
S4.1.3按式(2-4)的模型计算压气机压比;
S4.1.4按式(2-9)的模型计算机组恒速带负荷运行时的压气机等熵效率;
S4.1.8按计算压气机耗功,空气焓i由物性数据库给出;
S4.2燃烧室及燃气透平性能模拟,包括:
S4.2.2按燃烧室热平衡计算天然气流量Gf ;
S4.2.7按式(2-13)的模型计算机组恒速带负荷运行时燃气透平等熵效率;
S4.3燃气轮机性能:
S4.3.1按计算燃气轮机输出功,η m 为机械效率;
所述的步骤S5 研究燃气轮机环境温度特性
根据上述变工况模型的计算结果画成相应的关系曲线图,采用以下模型研究燃气轮机环境温度特性:
(1)功率-环境温度特性
当以ISO条件燃气轮机满载功率(即IGV开度100%)为基准时,燃气轮机相对功率Pgt/Pgt0 与环境温度及IGV开度之间的关系采用如下形式:
采用式(2-14),可将特定环境温度及负荷下的燃气轮机功率Pgt 折算到ISO条件下的满载功率Pgt0 。
当以ISO条件燃气轮机某负荷下的功率(Pgt0,IGV )为基准时,燃气轮机相对功率Pgt/Pgt0,IGV 与环境温度及IGV开度之间的关系采用如下形式:
(2-15)
采用式(2-15),可将各环境温度下的燃气轮机负荷修正到基准温度条件下的负荷。
(2)效率-环境温度特性
当以ISO条件燃气轮机满载效率(即IGV开度100%)为基准时,燃气轮机相对效率η gt / η gt0 与环境温度及IGV开度之间的关系采用如下形式:
(2-16)
采用式(2-16)可将特定环境温度及负荷下的燃气轮机功率η gt 折算到ISO条件下的满载功率η gt0 。
当以ISO条件燃气轮机某负荷下的效率(η gt0,IGV )为基准时,燃气轮机相对效率η gt / η gt0,IGV 与环境温度及IGV开度之间的关系采用如下形式:
采用式(2-17),可将各环境温度下的燃气轮机效率修正到基准温度条件下的对应效率值。
(3)排气参数-环境温度特性
当以ISO条件燃气轮机满载参数(即IGV开度100%)为基准时,燃气轮机相对排气流量Ggt/Ggt0 及相对排气温度t4/t40 与环境温度之间的关系采用如下形式:
(4)压气机的压比-环境温度特性
当以ISO条件燃气轮机满载参数(即IGV开度100%)为基准时,压气机相对压比π c / π c0 与环境温度之间的关系采用如下形式:
式(2-14)~ 式(2-20)中的常数a 、 b 、 c 、 d 、 e 、 f均与相应的关系曲线有关。
S6
研究燃气轮机进、排气阻力特性
根据上述变工况模型的计算结果画成相应的关系曲线图,采用以下模型研究燃气轮机进、排气阻力特性:
(1)进气阻力对燃气轮机功率及效率的影响
当以ISO条件燃气轮机满载参数(即IGV开度100%)为基准时,燃气轮机相对功率Pgt/Pgt0 及相对效率η gt / η gt0 与压气机进气阻力∆Pin 之间的关系采用如下形式:
(2-21)
(2-22)
式中:常数a与相应的关系曲线有关。
(2)排气阻力对燃气轮机功率及效率的影响
当以ISO条件燃气轮机满载参数(即IGV开度100%)为基准时,燃气轮机相对功率Pgt/Pgt0 及相对效率η gt / η gt0 与燃气透平排气阻力∆Pe 之间的关系采用如下形式:
式中:常数a与相应的关系曲线有关。
计算基准条件下燃气轮机性能的影响因子
(1)环境温度折算因子ζ P,ta
燃气轮机功率折算到基准温度条件的折算因子为:
环境压力、环境湿度等对燃气轮机性能的影响,直接采用供货商提供的性能曲线,其影响因子由相应的性能修正曲线得到。
因此,可将燃气轮机运行功率按下式折算到基准气象条件:
(2-26)
(4)基准效率-环境温度折算因子ζη ,ta
燃气轮机效率折算到基准温度条件的折算公式为:
(2-27)
其他因素,如发电机功率因数等,其修正因子忽略不计。
函数F1和F2为由变工况模型得到的关系。
S7
研究燃气轮机的负荷特性
为比较各台机组的负荷特性,需将燃气轮机功率及效率修正到基准环境条件,将变工况模型的计算结果(功率及效率)按式(2-26)和式(2-27)修正到基准环境条件,绘制燃气轮机效率-功率曲线,并对曲线数据进行拟合得到数学模型;
燃气轮机的负荷特性按下式拟合:
式中:常数a 、 b 、 c 、 d与绘制出来的燃气轮机效率-功率曲线有关,下标s表示基准条件。
有益效果:本发明提供的燃气轮机及其联合循环的变工况分析方法,可建立起变工况模型,获得带有IGV温度控制的燃气轮机恒速带负荷运行时的稳态变工况特性,根据模型运算的结果可对燃气轮机的环境温度特性和负荷特性进行研究。
附图说明
图1为本变工况分析方法的总流程图;
图2为建立变工况模型的具体流程图及具体采用的模型;
图3为燃气轮机环境温度特性、进排气阻力特性及负荷特性采用的模型概括图。
具体实施方式
下面以某带IGV温度控制的270MW燃气轮机为例,对本发明的具体实施方式进行说明。
本发明的燃气轮机及其联合循环的变工况分析方法,包括以下步骤:
S1 进行压气机进气歧管差压数据的合理性分析
采用压气机进气歧管静压法,计算压气机的空气流量,即按式(2-1)估算:
式中:Gc 为压气机入口空气流量,kg/s;Pin 为压气机进气歧管处的负压;∆ Pin 为过滤器阻力;ρ a 为空气密度;A为进气歧管负压测点所在的过流断面面积;Ф为流量系数;Ф *A可成为有效过流断面面积。
将压气机进气歧管差压Dp数据按IGV开度自低向高排列,分析其分布情况,得出部分不合理的数据需要剔除的结论。
鉴别歧管静压数据合理的步骤为:
观察Dp( 2 )的数据分布图,可知1次剔除已经足够。
S2
根据上一步得到的数据样本结合DCS数据进行燃气轮机运行性能测算,采用燃气轮机计算功率与DCS系统估算功率相对误差不超过3%的运行数据为建模样本。
S3
建立燃气轮机的变工况模型,包括以下步骤;
S3.1建立燃气轮机的IGV温控规律模型
综合分析燃气轮机排气温度T4 - IGV规律图和燃气轮机透平计算进气温度T3 - IGV规律图,得到燃气轮机的IGV温控规律如下:在燃气轮机基本负荷向部分负荷过渡的过程中,先关小IGV,减小压气机进气流量(透平排气温度T4 不变,而进气温度T3 略减小),再逐渐降低T3 。即:当IGV开度在39%~100%之间时,随着IGV开度的增大,燃气轮机按等排气温度运行,透平进气温度T3 较缓慢升至限定值;当IGV开度在大约39%以下时,随着IGV开度减小,燃气轮机功率降低,透平进气温度相对较快速下降。
机组带负荷稳定运行时,透平进气温度T3 与IGV开度之间的温控规律可描述为:
S3.2建立压气机压比π c ~ ( IGV , ta )模型
分析压气机压比随IGV开度的变化图,机组带负荷稳定运行时,压气机压比π c ~ ( IGV , ta )关系可描述为:
式中:常数a 、 b 、 c 、 d 、 e 、 f的数值由图中曲线拟合得到。
S3.3建立压气机流量Gc ~ ( IGV , ta )模型
分析压气机空气流量测算值(采用歧管静压法)随IGV开度的变化图,机组带负荷稳定运行时,压气机流量Gc ~ ( IGV , ta )关系可描述为:
式中:常数a 、 b 、 c 、 d 、 e 、 f的数值由图中曲线拟合得到。
S3.4建立压气机进气阻力模型
压气机进气阻力主要考虑进气过滤器(两级)的流动阻力,上述阻力可表达为:
式中:ζ为管件的阻力系数;c为流速;ta 为环境温度,℃;Pa为环境压力,Pa;Gc 为压气机空气质量流量,kg/s;常数a与机组设计性能及运行状况有关,由实际运行数据拟合得到。
S3.5建立IGV温控规律下的压气机效率-流量特性曲线
机组在IGV温控规律下、恒速运行时,压气机效率-流量(环境压力定为101.3kPa)关系可描述为:
式中:∆p为进气系统的阻力,kPa;Pa为环境压力,kPa。系数a 、 b 、 c由实际运行数据的特性点簇得到。此式为已把压气机的折合流量和折合转速设计值代入所得到的模型。
S3.6建立燃气透平的膨胀比模型
分析燃气透平膨胀比π t 与压气机压比π c 关系图,它们的关系可以描述为:
式中:常数a 、 b由图中曲线拟合得到。
S3.7建立燃气透平的流量模型
燃气透平的流量采用修正的Flugel公式,其形式:
式中:Gt0 、 T30 、 P30 、 P40 为ISO条件下的设计值。在本实施例中Gt0 =623.6kg/s,T30 =1564K,P30 =1687.7kPa,P40 =101.3kPa。因此,上式转化为:
S3.8建立燃气透平等熵效率模
燃气透平的等熵效率采用下述模型:
上式为已把燃气透平的折合流量和折合转速设计值代入所得到的模型。
S4 进行燃气轮机整体变工况性能模拟,包括以下步骤:
S4.1压气机性能模拟:
S4.1.1给定环境温度ta 、环境压力Pa 、压气机IGV开度、燃料特性、燃气轮机设计参数;
S4.1.2按式(2-5)的模型计算压气机空气流量;
S4.1.3按式(2-4)的模型计算压气机压比;
S4.1.4按式(2-9)的模型计算机组恒速带负荷运行时的压气机等熵效率;
S4.2燃烧室及燃气透平性能模拟:
② 按燃烧室热平衡计算天然气流量Gf ;
⑦ 按式(2-13)的模型计算机组恒速带负荷运行时燃气透平等熵效率;
S4.3燃气轮机性能:
S5研究燃气轮机环境温度特性。
根据上述变工况模型的计算结果画成相应的关系曲线图,采用以下模型研究燃气轮机环境温度特性:
(1)功率-环境温度特性
当以ISO条件燃气轮机满载功率(即IGV开度100%)为基准时,燃气轮机相对功率Pgt/Pgt0 与环境温度及IGV开度之间的关系采用如下形式:
当以ISO条件燃气轮机某负荷下的功率(Pgt0,IGV )为基准时,燃气轮机相对功率Pgt/Pgt0,IGV 与环境温度及IGV开度之间的关系采用如下形式:
(2)效率-环境温度特性
当以ISO条件燃气轮机满载效率(即IGV开度100%)为基准时,燃气轮机相对效率η gt / η gt0 与环境温度及IGV开度之间的关系采用如下形式:
当以ISO条件燃气轮机某负荷下的效率(η gt0,IGV )为基准时,燃气轮机相对效率η gt / η gt0,IGV 与环境温度及IGV开度之间的关系采用如下形式:
(3)排气参数-环境温度特性
当以ISO条件燃气轮机满载参数(即IGV开度100%)为基准时,燃气轮机相对排气流量Ggt/Ggt0 及相对排气温度t4/t40 与环境温度之间的关系采用如下形式:
(4)压气机的压比-环境温度特性
当以ISO条件燃气轮机满载参数(即IGV开度100%)为基准时,压气机相对压比π c / π c0 与环境温度之间的关系采用如下形式:
上述各式中的常数a 、 b 、 c 、 d 、 e 、 f均由相应的关系曲线拟合得到。
S6研究燃气轮机进、排气阻力特性。
根据上述变工况模型的计算结果画成相应的关系曲线图,采用以下模型研究燃气轮机进、排气阻力特性:
(1)进气阻力对燃气轮机功率及效率的影响
当以ISO条件燃气轮机满载参数(即IGV开度100%)为基准时,燃气轮机相对功率Pgt/Pgt0 及相对效率η gt / η gt0 与压气机进气阻力∆Pin 之间的关系采用如下形式:
(2)排气阻力对燃气轮机功率及效率的影响
当以ISO条件燃气轮机满载参数(即IGV开度100%)为基准时,燃气轮机相对功率Pgt/Pgt0 及相对效率η gt / η gt0 与燃气透平排气阻力∆Pe 之间的关系采用如下形式:
上述各式中的常数a均由相应的关系曲线拟合得到。
计算基准条件下燃气轮机性能的影响因子。
对燃气轮机的运行功率依次按照环境温度折算因子ζ P,,ta 、环境压力折算因子ζ p,pa 、环境相对湿度折算因子ζ p,RH 折算到基准条件(环境温度ta,s =27.8℃、环境压力Pa,s =1001.8hPa、相对湿度RHs=82%)。
燃气轮机的测算效率按照环境温度折算因子ζη ,ta 折算到基准条件。
S7研究燃气轮机的负荷特性
将折算后得到的基准条件下的燃气轮机运行功率和测算效率绘制成关系图,它们的关系可以描述为:
式中:常数a 、 b 、 c 、 d由绘制出来的燃气轮机效率-功率曲线拟合得到,下标s表示基准条件。
Claims (6)
1.一种基于解析法的燃气轮机及其联合循环的变工况分析方法,包括以下步骤:
S1 进行压气机进气歧管差压数据的合理性分析,剔除不合理数据;
S2 根据步骤S1得到的数据样本,结合DCS分散控制系统数据进行燃气轮机运行性能测算;
S3 建立燃气轮机的变工况模型,包括以下步骤;
S3.1建立燃气轮机的IGV温控规律模型;
S3.2建立压气机压比π c ~ ( IGV , ta )模型;
S3.3建立压气机流量Gc ~ ( IGV , ta )模型;
S3.4建立压气机进气阻力模型;
S3.5建立IGV温控规律下的压气机效率-流量特性曲线;
S3.6建立燃气透平的膨胀比模型;
S3.7建立燃气透平的流量模型;
S3.8建立燃气透平等熵效率模型;
S4 进行燃气轮机整体变工况性能模拟;
S5 根据步骤S4的结果画出燃气轮机整体变工况曲线图,研究燃气轮机环境温度特性:包括功率-环境温度特性、效率-环境温度特性、排气参数-环境温度特性、压气机的压比-环境温度特性;
S6 根据步骤S4的结果研究进排气阻力对燃气轮机性能的影响;
S7 根据步骤S4的结果研究燃气轮机的负荷特性,包括:
S7.1根据步骤S5的结果计算基准条件下燃气轮机性能的影响因子:基准功率-环境温度折算因子、基准功率-环境压力折算因子、基准功率-相对湿度折算因子、基准效率-环境温度折算因子;
S7.2根据步骤S7.1的结果将燃气轮机功率及效率修正到基准环境条件,绘制燃气轮机折算效率-折算功率曲线,并对曲线数据进行拟合,最后得到燃气轮机折算效率—折算功率的数学模型。
2.根据权利要求1所述的基于解析法的燃气轮机及其联合循环的变工况分析方法,其特征是:所述的步骤S1中的进行压气机进气歧管差压数据的合理性分析为:
采用压气机进气歧管静压法,计算压气机的空气流量,即按式(2-1)估算:
式中:Gc 为压气机入口空气流量,kg/s;Pin 为压气机进气歧管处的负压;∆Pin 为过滤器阻力;ρ a 为空气密度;A为进气歧管负压测点所在的过流断面面积;Ф为流量系数;Ф *A可成为有效过流断面面积;
由于压气机进气歧管静压测量值可能有较大误差,导致数据波动,故需要剔除明显错误的数据,本发明提出一种用于鉴别进气歧管差压值合理性的判别式形式:
(2-2)
所述的步骤S1中的鉴别歧管静压数据合理的步骤为:
③重复步骤②,直到数据样本满足要求。
3.根据权利要求1所述的基于解析法的燃气轮机及其联合循环的变工况分析方法,其特征是:所述的步骤S4包括以下步骤:
S4.1压气机性能模拟:
S4.1.1给定环境温度ta 、环境压力Pa 、压气机IGV开度、燃料特性、燃气轮机设计参数;
S4.1.2按式(2-5)的模型计算压气机空气流量;
S4.1.3按式(2-4)的模型计算压气机压比;
S4.1.4按式(2-9)的模型计算机组恒速带负荷运行时的压气机等熵效率;
S4.1.6按计算压气机等熵压缩终排气温度;
S4.2燃烧室及燃气透平性能模拟,包括:
S4.2.2按燃烧室热平衡计算天然气流量Gf ;
S4.2.7按式(2-13)的模型计算机组恒速带负荷运行时燃气透平等熵效率;
S4.2.10按计算燃气透平输出功,燃气焓i由物性数据库给出;
S4.3燃气轮机性能:
4.根据权利要求1所述的基于解析法的燃气轮机及其联合循环的变工况分析方法,其特征是:所述的步骤S5采用以下模型研究燃气轮机环境温度特性:
(1)功率-环境温度特性
当以ISO条件燃气轮机满载功率(即IGV开度100%)为基准时,燃气轮机相对功率Pgt/Pgt0 与环境温度及IGV开度之间的关系采用如下形式:
采用式(2-14),可将特定环境温度及负荷下的燃气轮机功率Pgt 折算到ISO条件下的满载功率Pgt0 ;
当以ISO条件燃气轮机某负荷下的功率(Pgt0,IGV )为基准时,燃气轮机相对功率Pgt/Pgt0,IGV 与环境温度及IGV开度之间的关系采用如下形式:
(2-15)
采用式(2-15),可将各环境温度下的燃气轮机负荷修正到基准温度条件下的负荷;
(2)效率-环境温度特性
当以ISO条件燃气轮机满载效率(即IGV开度100%)为基准时,燃气轮机相对效率η gt / η gt0 与环境温度及IGV开度之间的关系采用如下形式:
采用式(2-16)可将特定环境温度及负荷下的燃气轮机功率η gt 折算到ISO条件下的满载功率η gt0 ;
当以ISO条件燃气轮机某负荷下的效率(η gt0,IGV )为基准时,燃气轮机相对效率η gt / η gt0,IGV 与环境温度及IGV开度之间的关系采用如下形式:
采用式(2-17),可将各环境温度下的燃气轮机效率修正到基准温度条件下的对应效率值;
(3)排气参数-环境温度特性
当以ISO条件燃气轮机满载参数(即IGV开度100%)为基准时,燃气轮机相对排气流量Ggt/Ggt0 及相对排气温度t4/t40 与环境温度之间的关系采用如下形式:
(4)压气机的压比-环境温度特性
当以ISO条件燃气轮机满载参数(即IGV开度100%)为基准时,压气机相对压比π c / π c0 与环境温度之间的关系采用如下形式:
(2-20)
式(2-14)~式(2-20)中的常数a 、 b 、 c 、 d 、 e 、 f均与相应的关系曲线有关。
5.根据权利要求1所述的基于解析法的燃气轮机及其联合循环的变工况分析方法,其特征是:所述的步骤S6为:
根据上述变工况模型的计算结果画成相应的关系曲线图,采用以下模型研究燃气轮机进、排气阻力特性:
(1)进气阻力对燃气轮机功率及效率的影响
当以ISO条件燃气轮机满载参数-IGV开度100%为基准时,燃气轮机相对功率Pgt/Pgt0 及相对效率η gt / η gt0 与压气机进气阻力∆Pin 之间的关系采用如下形式:
式中:常数a与相应的关系曲线有关;
(2)排气阻力对燃气轮机功率及效率的影响
当以ISO条件燃气轮机满载参数(即IGV开度100%)为基准时,燃气轮机相对功率Pgt/Pgt0 及相对效率η gt / η gt0 与燃气透平排气阻力∆Pe 之间的关系采用如下形式:
式中:常数a与相应的关系曲线有关;
计算基准条件下燃气轮机性能的影响因子
(1)环境温度折算因子ζ P,ta
燃气轮机功率折算到基准温度条件的折算因子为:
环境压力、环境湿度等对燃气轮机性能的影响,直接采用供货商提供的性能曲线,其影响因子由相应的性能修正曲线得到;
因此,可将燃气轮机运行功率按下式折算到基准气象条件:
(2-26)
(4)基准效率-环境温度折算因子ζη ,ta
燃气轮机效率折算到基准温度条件的折算公式为:
其他因素,如发电机功率因数等,其修正因子忽略不计;
函数F1和F2为由变工况模型得到的关系。
6.根据权利要求1所述的基于解析法的燃气轮机及其联合循环的变工况分析方法,其特征是:所述的步骤S7包括:
为比较各台机组的负荷特性,需将燃气轮机功率及效率修正到基准环境条件,将变工况模型的计算结果(功率及效率)按式(2-26)和式(2-27)修正到基准环境条件,绘制燃气轮机效率-功率曲线,并对曲线数据进行拟合得到数学模型;
燃气轮机的负荷特性按下式拟合:
(2-28)
式中:常数a 、 b 、 c 、 d与绘制出来的燃气轮机效率-功率曲线有关,下标s表示基准条件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210197135.2A CN102877945B (zh) | 2012-06-15 | 2012-06-15 | 基于解析法的燃气轮机及其联合循环的变工况分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210197135.2A CN102877945B (zh) | 2012-06-15 | 2012-06-15 | 基于解析法的燃气轮机及其联合循环的变工况分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102877945A true CN102877945A (zh) | 2013-01-16 |
CN102877945B CN102877945B (zh) | 2014-07-16 |
Family
ID=47479423
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210197135.2A Active CN102877945B (zh) | 2012-06-15 | 2012-06-15 | 基于解析法的燃气轮机及其联合循环的变工况分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102877945B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103294892A (zh) * | 2013-04-01 | 2013-09-11 | 深圳市广前电力有限公司 | 基于燃气轮机温控作用下的压气机在线特性分析方法 |
CN103593570A (zh) * | 2013-11-19 | 2014-02-19 | 中国神华能源股份有限公司 | 一种确定燃气轮机最佳离线水洗周期的方法 |
CN105787211A (zh) * | 2016-03-24 | 2016-07-20 | 华南理工大学 | 针对燃气透平劣化的联合循环余热锅炉压力调整方法 |
CN107301268A (zh) * | 2017-05-18 | 2017-10-27 | 哈尔滨工程大学 | 一种船舶燃气轮机可调静叶压气机偏转角度优化方法 |
CN109655275A (zh) * | 2018-11-29 | 2019-04-19 | 沈阳工程学院 | 一种基于敏感系数法的燃气轮机循环热经济性诊断方法 |
CN110905684A (zh) * | 2019-11-25 | 2020-03-24 | 湖南工程学院 | 海上采油平台透平主机回热循环效率优化方法 |
CN110928248A (zh) * | 2018-09-19 | 2020-03-27 | 上海明华电力技术工程有限公司 | 一种确定燃气轮机性能劣化程度的方法 |
CN110989342A (zh) * | 2019-11-19 | 2020-04-10 | 华北电力大学 | 一种联合循环机组重型燃气轮机实时t-s模糊建模方法 |
CN111638003A (zh) * | 2020-05-06 | 2020-09-08 | 中国船舶重工集团公司第七O三研究所无锡分部 | 水力测功器更换后燃气轮机功率测量出现偏差的修正方法 |
CN112685847A (zh) * | 2020-12-15 | 2021-04-20 | 苏州西热节能环保技术有限公司 | 一种用于确定联合循环机组燃气轮机排气压力的方法 |
CN113640115A (zh) * | 2021-08-11 | 2021-11-12 | 中国工程物理研究院流体物理研究所 | 适用于准等熵压缩实验数据逆问题求解的优化方法和系统 |
CN114048423A (zh) * | 2021-09-30 | 2022-02-15 | 华电电力科学研究院有限公司 | 一拖一单轴燃气轮机发电机组压气机瞬时空气流量计算方法 |
CN114893305A (zh) * | 2022-05-20 | 2022-08-12 | 华电电力科学研究院有限公司 | 一种燃气机组的控制方法和系统 |
CN118277937A (zh) * | 2024-05-29 | 2024-07-02 | 华能海南发电股份有限公司南山电厂 | 生产运行数据的分析方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1004760A2 (en) * | 1998-11-27 | 2000-05-31 | Mazda Motor Corporation | Control apparatus of engine having turbo supercharger |
JP2001330121A (ja) * | 2000-05-23 | 2001-11-30 | Toyota Motor Corp | 車両用無段変速機の制御装置 |
CN1704572A (zh) * | 2004-05-31 | 2005-12-07 | 宝山钢铁股份有限公司 | 燃气-蒸汽联合循环发电站运行效率在线解析方法和系统 |
US20110137474A1 (en) * | 2009-12-16 | 2011-06-09 | Einar Vaughn Larsen | Method and systems for operating a wind turbine |
-
2012
- 2012-06-15 CN CN201210197135.2A patent/CN102877945B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1004760A2 (en) * | 1998-11-27 | 2000-05-31 | Mazda Motor Corporation | Control apparatus of engine having turbo supercharger |
JP2001330121A (ja) * | 2000-05-23 | 2001-11-30 | Toyota Motor Corp | 車両用無段変速機の制御装置 |
CN1704572A (zh) * | 2004-05-31 | 2005-12-07 | 宝山钢铁股份有限公司 | 燃气-蒸汽联合循环发电站运行效率在线解析方法和系统 |
US20110137474A1 (en) * | 2009-12-16 | 2011-06-09 | Einar Vaughn Larsen | Method and systems for operating a wind turbine |
Non-Patent Citations (1)
Title |
---|
阚伟民: "《M701F型燃气轮机的在线运行负荷特性》", 《广东电力》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103294892A (zh) * | 2013-04-01 | 2013-09-11 | 深圳市广前电力有限公司 | 基于燃气轮机温控作用下的压气机在线特性分析方法 |
CN103593570A (zh) * | 2013-11-19 | 2014-02-19 | 中国神华能源股份有限公司 | 一种确定燃气轮机最佳离线水洗周期的方法 |
CN103593570B (zh) * | 2013-11-19 | 2016-09-28 | 中国神华能源股份有限公司 | 一种确定燃气轮机最佳离线水洗周期的方法 |
CN105787211A (zh) * | 2016-03-24 | 2016-07-20 | 华南理工大学 | 针对燃气透平劣化的联合循环余热锅炉压力调整方法 |
CN105787211B (zh) * | 2016-03-24 | 2019-01-29 | 华南理工大学 | 针对燃气透平劣化的联合循环余热锅炉压力调整方法 |
CN107301268A (zh) * | 2017-05-18 | 2017-10-27 | 哈尔滨工程大学 | 一种船舶燃气轮机可调静叶压气机偏转角度优化方法 |
CN110928248B (zh) * | 2018-09-19 | 2022-10-25 | 上海明华电力科技有限公司 | 一种确定燃气轮机性能劣化程度的方法 |
CN110928248A (zh) * | 2018-09-19 | 2020-03-27 | 上海明华电力技术工程有限公司 | 一种确定燃气轮机性能劣化程度的方法 |
CN109655275A (zh) * | 2018-11-29 | 2019-04-19 | 沈阳工程学院 | 一种基于敏感系数法的燃气轮机循环热经济性诊断方法 |
CN110989342A (zh) * | 2019-11-19 | 2020-04-10 | 华北电力大学 | 一种联合循环机组重型燃气轮机实时t-s模糊建模方法 |
CN110989342B (zh) * | 2019-11-19 | 2021-04-16 | 华北电力大学 | 一种联合循环机组重型燃气轮机实时t-s模糊建模方法 |
CN110905684A (zh) * | 2019-11-25 | 2020-03-24 | 湖南工程学院 | 海上采油平台透平主机回热循环效率优化方法 |
CN111638003A (zh) * | 2020-05-06 | 2020-09-08 | 中国船舶重工集团公司第七O三研究所无锡分部 | 水力测功器更换后燃气轮机功率测量出现偏差的修正方法 |
CN111638003B (zh) * | 2020-05-06 | 2021-10-15 | 中国船舶重工集团公司第七O三研究所无锡分部 | 水力测功器更换后燃气轮机功率测量出现偏差的修正方法 |
CN112685847A (zh) * | 2020-12-15 | 2021-04-20 | 苏州西热节能环保技术有限公司 | 一种用于确定联合循环机组燃气轮机排气压力的方法 |
CN112685847B (zh) * | 2020-12-15 | 2024-05-24 | 苏州西热节能环保技术有限公司 | 一种用于确定联合循环机组燃气轮机排气压力的方法 |
CN113640115A (zh) * | 2021-08-11 | 2021-11-12 | 中国工程物理研究院流体物理研究所 | 适用于准等熵压缩实验数据逆问题求解的优化方法和系统 |
CN114048423A (zh) * | 2021-09-30 | 2022-02-15 | 华电电力科学研究院有限公司 | 一拖一单轴燃气轮机发电机组压气机瞬时空气流量计算方法 |
CN114893305A (zh) * | 2022-05-20 | 2022-08-12 | 华电电力科学研究院有限公司 | 一种燃气机组的控制方法和系统 |
CN114893305B (zh) * | 2022-05-20 | 2023-05-12 | 华电电力科学研究院有限公司 | 一种燃气机组的控制方法和系统 |
CN118277937A (zh) * | 2024-05-29 | 2024-07-02 | 华能海南发电股份有限公司南山电厂 | 生产运行数据的分析方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN102877945B (zh) | 2014-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102877945B (zh) | 基于解析法的燃气轮机及其联合循环的变工况分析方法 | |
CN103133147B (zh) | 用于改变燃气涡轮系统中的功率输出的方法和设备 | |
CN109145363B (zh) | 燃料电池系统中离心式空气压缩机瞬态建模方法 | |
CN103630364B (zh) | 一种模拟高原环境试验柴油机的方法 | |
CN105370629B (zh) | Pta装置用能量回收控制方法 | |
CN103294892A (zh) | 基于燃气轮机温控作用下的压气机在线特性分析方法 | |
CN104533821A (zh) | 一种压气机控制方法 | |
CN101545495A (zh) | 鼓风机的防喘振控制装置及方法 | |
CN102767529B (zh) | 一种轴流压缩机节能效果的计算方法 | |
CN107849981A (zh) | 燃气轮机的控制装置及方法、燃气轮机的控制程序、燃气轮机 | |
CN109670244B (zh) | 一种涡轴发动机翻修后燃气涡轮导向器面积调整方法 | |
CN114060143B (zh) | 基于旁通阀流通特性的变海拔增压压力稳定性控制方法 | |
CN105736349B (zh) | 罗茨真空泵性能测试系统及其测试方法 | |
CN110617233A (zh) | 天然气长输管道压缩机组负荷分配控制系统 | |
CN114583766B (zh) | 面向风电消纳的海上油气平台能源系统低碳运行调度方法 | |
CN109583012A (zh) | 先进绝热压缩空气储能与风电协同运行调度方法与装置 | |
CN109209979A (zh) | 基于变工况运行的压缩机动态防喘振的控制方法 | |
RU2294460C2 (ru) | Способ оптимирования эксплуатации множества компрессорных агрегатов компрессорной станции природного газа | |
CN103498788A (zh) | 重型压气机空气流量和效率在线测量方法及装置 | |
CN111008491A (zh) | 一种真空系统在真空泵作用下压力性能的优化方法 | |
CN115614304A (zh) | 基于ladrc的高空台抽气机组自动并网控制方法 | |
CN113050450B (zh) | 一种并联变频泵输配系统仿真模块编写方法 | |
CN110472277B (zh) | 一种汽轮机的建模方法 | |
Tajalli et al. | Thermodynamic simulation of two-shaft gas turbine to study invasive weeds optimization and Min-Max controller strategies considering air-cooled blades | |
CN104732445B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |