CN113239493A - 燃气轮机启机过程NOx排放性能建模和软件开发设计方法 - Google Patents
燃气轮机启机过程NOx排放性能建模和软件开发设计方法 Download PDFInfo
- Publication number
- CN113239493A CN113239493A CN202110570553.0A CN202110570553A CN113239493A CN 113239493 A CN113239493 A CN 113239493A CN 202110570553 A CN202110570553 A CN 202110570553A CN 113239493 A CN113239493 A CN 113239493A
- Authority
- CN
- China
- Prior art keywords
- gas turbine
- heavy
- nox emission
- state change
- cantera
- 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
- 238000000034 method Methods 0.000 title claims abstract description 82
- 230000008569 process Effects 0.000 title claims abstract description 51
- 238000013461 design Methods 0.000 title claims abstract description 26
- 230000008859 change Effects 0.000 claims abstract description 86
- 238000002485 combustion reaction Methods 0.000 claims abstract description 76
- 239000000126 substance Substances 0.000 claims abstract description 24
- 239000000446 fuel Substances 0.000 claims description 61
- 238000009826 distribution Methods 0.000 claims description 16
- 238000006243 chemical reaction Methods 0.000 claims description 14
- 230000007246 mechanism Effects 0.000 claims description 6
- 238000010792 warming Methods 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 5
- 238000002347 injection Methods 0.000 claims description 4
- 239000007924 injection Substances 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 4
- 239000003344 environmental pollutant Substances 0.000 abstract description 15
- 231100000719 pollutant Toxicity 0.000 abstract description 15
- 238000002474 experimental method Methods 0.000 abstract description 4
- 238000001514 detection method Methods 0.000 abstract description 2
- 239000007789 gas Substances 0.000 description 73
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 7
- 239000000243 solution Substances 0.000 description 6
- 238000011160 research Methods 0.000 description 5
- 230000003068 static effect Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 241000894007 species Species 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 239000000571 coke Substances 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
Abstract
本发明提供了一种燃气轮机启机过程NOx排放性能建模和软件开发设计方法，包括：选取目标为重型燃气轮机燃烧室并采集所述重型燃气轮机燃烧室选型设计参数的初始数据信息：将所述重型燃气轮机燃烧室划分为若干区域；根据若干所述区域的特点利用开源程序包Cantera建立燃烧室的化学反应器网络模型；根据重型燃气轮机启机过程中各个参数的变化，划分出若干个状态变化点；建立NOx预估模型。设计人员能根据检测过程中的准确数据，自动化的算出启机各个阶段污染物的排放值，较原有的实验方法省时、省力、省资金，有较强的准确性和可实施性，这个方法可以合理、准确给出该重型燃气轮机在启机过程中污染物排放的预估值。
Description
技术领域
本发明涉及燃气轮机污染物减排技术领域，具体而言，涉及一种燃气轮机启机过程NOx排放性能建模和软件开发设计方法。
背景技术
随着环境问题不断增长以及对能源需求的增加，人类社会正朝着保护地球生态系统的方向发展。近年来我国一直在大力发展可再生能源，但风能、太阳能和氢能等前期成本高，研究周期长，这些技术进入能源市场相对缓慢。与新能源相比，传统燃料的燃烧依然是目前最为常见的能源使用方式。近些年来随着我国西气东输工程的全面竣工，使得更多的地区将天然气这种清洁能源代替传统的燃煤，这大大改善了能源消耗的整体结构，在全国范围内使用清洁能源的比例显著提升。
燃烧室是燃气轮机重要组成部分之一，同时也是产生污染物的源头，燃烧室性能指标，如热力性能、污染物排放、使用寿命和运行可靠性等，都会对燃汽轮机产生重要的影响。由于近年来对燃气轮机的性能要求不断增加的需要，燃气轮机的功率和流量也持续增加，使燃烧室内的工作压力与温度也不断提高，最终温度和压力的增加会加剧NOx的产生。伴随人类环保意识的增强和排放标准的日趋严格，降低污染物排放特别是NOx的排放已成为燃气轮机技术发展的一个突出问题，因此对燃烧室内污染物排放预估和软件开发变得尤为重要。低NOx技术的发展可为我国燃气轮机的设计提供研究方向，对于研究重型燃机低污染物排放的性能建模也具有重要的指导作用。预测软件在设计过程中的意义十分重大，因此研究重型燃机的污染物排放的模型预估十分重要。
发明内容
本发明旨在至少解决现有技术中存在的技术问题之一。
为此，本发明提供了一种重型燃气轮机启机过程NOx排放性能建模和软件开发设计方法。
本发明提供了一种燃气轮机启机过程NOx排放性能建模和软件开发设计方法，所述方法包括：
步骤1：选取目标为重型燃气轮机燃烧室并采集所述重型燃气轮机燃烧室选型设计参数的初始数据信息：
步骤2：根据与所述重型燃气轮机相近的温度分布云图和速度分布云图，将所述重型燃气轮机燃烧室划分为若干区域；根据若干所述区域的特点利用开源程序包Cantera建立燃烧室的化学反应器网络模型；
步骤3：根据重型燃气轮机启机过程中各个参数的变化，划分出若干个状态变化点；
步骤4：根据若干所述状态变化点不同的工况，并结合所述化学反应器网络模型，计算出各个所述状态变化点的NOx排放值，将各个所述状态变化点的NOx排放值连成曲线，建立NOx预估模型。
根据本发明上述技术方案的燃气轮机启机过程NOx排放性能建模和软件开发设计方法，还可以具有以下附加技术特征：
在上述技术方案中，所述步骤1中的燃烧室选型设计参数为重型燃气轮机燃烧室总体设计人员在燃烧室选型过程中会用到的选型设计参数，所述参数包括：
进口主燃料量、进口值班燃料量、进口空气流量、二次空气流量、压气机压比、燃烧室体积、燃烧室温度和燃烧筒个数。
在上述技术方案中，所述步骤2中的温度分布云图和速度分布云图通过 CFD仿真计算流体力学软件建立。其中，所述CFD仿真计算流体力学软件可以是ANSYS Fluent、cfx、Phoenics、OpenFOAM以及FloEFD中的一种或几种。用于建立所述温度分布云图和速度分布云图。
在上述技术方案中，所述步骤3中的启机过程中各个参数包括：
机组转速、机组负荷、主燃料量、值班燃料量、二次空气流量、旁路阀开度和NOx排放值。
在上述技术方案中，所述步骤3中的状态变化点通过所述机组转速、所述机组负荷、所述主燃料量以及所述旁路阀开度的变化建立，即通过上述4个参数的变化为基础建立状态变化点，且状态变化点数量为12个；
其中，所述状态变化点至少包括：点火点、机组定速点、并网点、暖机负荷点以及满负荷点；
在上述技术方案中，所示步骤4具体包括如下步骤：
步骤4-1：根据所述化学反应器模型建立对应的Cantera模型，调用各个区域所对应反应器的Cantera程序包；
步骤4-2：在每个反应器的Cantera程序包中，输入反应条件；
步骤4-3：在所述Cantera程序包中选择适配于所述重型燃气轮机燃烧室内燃烧的反应机理；
步骤4-4：根据所述化学反应器网络模型，将各个反应器的Cantera程序段建立连接，得到整个化学反应器模型的Cantera程序包；
步骤4-5：根据步骤3的各个工况，利用步骤4-4得到的Cantera程序包计算出若干所述状态变化点的NOx排放值，将若干所述状态变化点NOx排放值数据进行整合，得出重型燃气轮机启机过程中NOx排放预估模型。
在上述技术方案中，所述步骤4-2中的反应条件包括：
温度、压力、燃料流量、空气流量、反应器内物种组分、反应器体积和停留时间。
在上述技术方案中，所述步骤4-4中的各个反应器的Cantera程序段通过 Python软件建立连接。
本发明的附加方面和优点将在下面的描述部分中变得明显，或通过本发明的实践了解到。
本发明提出的一种燃气轮机启机过程NOx排放性能建模和软件开发设计方法，设计人员能根据检测过程中的准确数据，自动化的算出启机各个阶段污染物的排放值，较原有的实验方法省时、省力、省资金，有较强的准确性和可实施性，这个方法可以合理、准确给出该重型燃气轮机在启机过程中污染物排放的预估值。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1为本发明一种燃气轮机启机过程NOx排放性能建模和软件开发设计方法流程图；
图2为本发明实施步骤4的子流程图；
图3为本发明重型燃气轮机的温度分布云图；
图4为本发明重型燃气轮机的速度分布云图；
图5为本发明重型燃气轮机燃烧室的区域划分图；
图6为本发明重型燃气轮机燃烧室化学反应器网络模型图；
图7为本发明重型燃气轮机启机过程中机组转速的变化；
图8为本发明重型燃气轮机启机过程中机组负荷的变化；
图9为本发明重型燃气轮机启机过程中燃烧室主燃料量的变化；
图10为本发明重型燃气轮机启机过程中旁路阀开度的变化；
图11为本发明重型燃气轮机启机过程NOx排放量变化与现有实验值对比图。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是，在不冲突的情况下，本申请的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述的方式来实施，因此，本发明的保护范围并不受下面公开的具体实施例的限制。
下面参照图1至图11来描述根据本发明一种重型燃气轮机启机过程 NOx排放性能建模和软件开发设计方法。
本发明第一个实施例提供了一种燃气轮机启机过程NOx排放性能建模和软件开发设计方法，所述方法包括：
步骤1：选取目标为重型燃气轮机燃烧室并采集所述重型燃气轮机燃烧室选型设计参数的初始数据信息：
步骤2：根据与所述重型燃气轮机相近的温度分布云图和速度分布云图，将所述重型燃气轮机燃烧室划分为若干区域；根据若干所述区域的特点利用开源程序包Cantera建立燃烧室的化学反应器网络模型；
步骤3：根据重型燃气轮机启机过程中各个参数的变化，划分出若干个状态变化点；
步骤4：根据若干所述状态变化点不同的工况，并结合所述化学反应器网络模型，计算出各个所述状态变化点的NOx排放值，将各个所述状态变化点的NOx排放值连成曲线，建立NOx预估模型。
在本实施例中，所述步骤1中的燃烧室选型设计参数为重型燃气轮机燃烧室总体设计人员在燃烧室选型过程中会用到的选型设计参数，所述参数包括：
进口主燃料量、进口值班燃料量、进口空气流量、二次空气流量、压气机压比、燃烧室体积、燃烧室温度和燃烧筒个数。
在本实施例中，所述步骤2中的温度分布云图和速度分布云图通过CFD 仿真计算流体力学软件建立。其中，所述CFD仿真计算流体力学软件可以是 ANSYS Fluent、cfx、Phoenics、OpenFOAM以及FloEFD中的一种或几种。用于建立所述温度分布云图和速度分布云图。
在本实施例中，所述步骤3中的启机过程中各个参数包括：
机组转速、机组负荷、主燃料量、值班燃料量、二次空气流量、旁路阀开度和NOx排放值。
在本实施例中，所述步骤3中的状态变化点通过所述机组转速、所述机组负荷、所述主燃料量以及所述旁路阀开度的变化建立，即通过上述4个参数的变化为基础建立状态变化点，且状态变化点数量为12个；
其中，所述状态变化点至少包括：点火点、机组定速点、并网点、暖机负荷点以及满负荷点；
在本实施例中，所示步骤4具体包括如下步骤：
步骤4-1：根据所述化学反应器模型建立对应的Cantera模型，调用各个区域所对应反应器的Cantera程序包；
步骤4-2：在每个反应器的Cantera程序包中，输入反应条件；
步骤4-3：在所述Cantera程序包中选择适配于所述重型燃气轮机燃烧室内燃烧的反应机理；
步骤4-4：根据所述化学反应器网络模型，将各个反应器的Cantera程序段建立连接，得到整个化学反应器模型的Cantera程序包；
步骤4-5：根据步骤3的各个工况，利用步骤4-4得到的Cantera程序包计算出若干所述状态变化点的NOx排放值，将若干所述状态变化点NOx排放值数据进行整合，得出重型燃气轮机启机过程中NOx排放预估模型。
在本实施例中，所述步骤4-2中的反应条件包括：
温度、压力、燃料流量、空气流量、反应器内物种组分、反应器体积和停留时间。
在本实施例中，所述步骤4-4中的各个反应器的Cantera程序段通过Python 软件建立连接。
本发明第二个实施例提供了一种燃气轮机启机过程NOx排放性能建模和软件开发设计方法的具体实施例，包括：
步骤1：选取某一型号的重型燃气轮机为研究对象，采集该重型燃气轮机燃烧室选型设计参数的初始数据信息；
所述燃烧室选型设计参数为重型燃气轮机燃烧室总体设计人员在燃烧室选型过程中会用到的选型设计参数，所述参数包括：
进口主燃料量、进口值班燃料量、进口空气流量、二次空气流量、压气机压比、燃烧室体积、燃烧室温度和燃烧筒个数；
其中，本具体实施例以该重型燃气轮机满负荷工作状态为例，其选型设计参数包括：
进口主燃料量0.8969kg/s，进口值班流量0.076kg/s，进口空气流量 3.2895kg/s，二次空气流量5.3659kg/s，压气机压比17，燃烧筒个数20。
步骤2：以该重型燃气轮机满负荷工作下的温度分布云图和速度分布云图为参考对象，将燃烧室划分为若干区域，根据各个区域的特点利用开源程序包 Cantera建立燃烧室的化学反应器网络模型。
步骤3：分析启机过程中各个参数的变化曲线，划分出点火点、自持转速点、机组定速点、并网点、暖机负荷点等各种状态变化点共12个；
所述启机过程中参数包括：机组转速、机组负荷、主燃料量、值班燃料量、二次空气流量、旁路阀开度；
其中，各个状态变化点的建立以及工作状态描述如下：
1.状态变化点1，即点火点：如图7所示，为该重型燃气轮机启机过程转速变化，启动开始后一小段时间转速成线性增加，增加至700r/min，随后一段时间，转速开始减小，减小至560r/min，这一状态变化点为点火点，记为状态变化点1；
工作状态描述：清扫结束后，燃气轮机转速降至560r/min，主燃料以及值班燃料压力调节阀、流量调节阀开启，IGV(重型燃气轮机的进气导叶)由全开0％到半开39.4％，开始点火，燃料由于刚刚喷入，燃烧室内燃料很少，为了燃烧稳定，旁路阀开度保持100％开。
2.状态变化点2：点火成功后转速在一段时间内增长缓慢，但由于点火成功后，温度瞬间增加，燃烧室燃烧开始变剧烈，点火成功后这一点记为状态变化点2；
工作状态描述：点火成功后，当机组转速达到600r/min时，SFC(静态变频器)和透平共同做功，机组开始进入了升速阶段，这是主燃流量由 FLCSO(燃气轮机在升速过程中的燃料量)控制，IGV(重型燃气轮机的进气导叶)由全开0％到半开39.4％。
3.状态变化点3：如图10所示，为该重型燃气轮机旁路阀开度变化，在机组点火前，旁路阀开度为100％，随后旁路阀开度很快减小到30％，这时转速继续缓慢增加，将点火后旁路阀降至30％这一点记为状态变化点3；
工作状态描述：燃气轮机内喷入燃料升高，燃烧室内燃烧开始稳定，旁路阀开度由100％降至30％。
4.状态变化点4：当转速再次升至700r/min时，转速的增速开始变快达到135r/min，这时机组转速变化主要是SFC(静态变频器)的作用，将转速第二次达到700r/min这一点记为状态变化点4；
工作状态描述：转速再次达到700r/min时，透平做功开始增大，SFC (静态变频器)做功开始减小，转速的增速开始变快达到135r/min。
5.状态变化点5：点火成功后主燃流量开始增加，在机组升速期间主燃料量由FLCSO(燃气轮机在升速过程中的燃料量)控制，变化量由预定程度来决定，当主燃料量增加至0.127kg/s时，主燃料增长速度开始变快，这时转速仍以135r/min的速率增加，旁路阀开度正在逐渐增加，将主燃料流量0.127kg/s记为状态变化点5；
工作状态描述：由于FLCSO(燃气轮机在升速过程中的燃料量)控制的燃料是按程序设定变化的，当转速在0-1100r/min时，主燃流量流量增长较慢，转速在1100r/min时，喷入燃料的速度开始变快，由于燃烧室内喷入燃料逐渐增多，旁路阀开度已升高至60％。
6.状态变化点6：随后旁路阀开度再次增至100％，这一状态变化点记为状态变化点6；
工作状态描述：由于主燃料流量增长变快，燃烧室内燃烧剧烈，旁路阀开度瞬间增至100％，机组开始带负荷，此时转速达到1500r/min。
7.状态变化点7，即并网点：当机组转速增长到2170r/min时，SFC(静态变频器)不起作用，这一点为自持转速点，记为状态变化点7。
工作状态描述：启机过程达到30min时，机组负荷达到200MW，这一点成为并网点，转速达到自持转速2170r/min，燃气轮机轴仅靠透平做功即可加速，SFC(静态变频器)推出工作；
8.状态变化点8，即机组定速点：燃机启机后，机组温度未达到燃机带负荷条件，当温度达到条件后，机组负荷瞬间从0增加至20MW，机组为并网状态，此时转速达到2170r/min，与自持转速点相同，随后机组负荷开始呈线性增加，当机组转速增加至3000r/min时，机组负荷此时增加到 60MW，这一状态变化点为机组定速点，记为状态变化点8；
工作状态描述：转速达到2750r/min时，IGV(重型燃气轮机的进气导叶)开度由半开39.4％变为0％，进口空气流量发生变化。
9.状态变化点9：当主燃料量增至0.89kg/s后，燃气轮机内燃料开始稳定一段时间，这时负荷继续增大，转速为3000r/min，旁路阀开度依旧为 100％，将主燃料点0.89kg/s这一点记为状态变化点9；
工作状态描述：转速达到3000r/min时，机组定速，燃气轮机负荷达到 60MW，这是主燃料控制由LDCSO(负荷控制下决定的燃料量)进行，这一阶段主燃流量的变化和机组负荷的变化相同，且燃烧依旧剧烈，旁路阀开度保持在100％。
10.状态变化点10：启机达到状态变化点9后，主燃料量保持稳定一段时间后，燃料量突升至1.017kg/s，负荷继续增加，旁路阀开度由100％减小至30％，这一点记为状态变化点10；
工作状态描述：机组定速后，由于旁路空气长时间的作用，燃烧室内燃烧开始稳定，燃烧室内转为预混燃烧为主导，旁路阀开度开始降低。
11.状态变化点11，即暖机负荷点：当负荷增长至120MW时，旁路阀开度依旧保持30％，这时负荷暂时稳定在120MW，整个机组处于暖机状态，将负荷刚好增长至120MW这一点成为暖机负荷点，记为状态变化点11；
工作状态描述：当负荷增加到120MW时，机组进入了暖机负荷状态，旁路阀开度保持在30％，IGV(重型燃气轮机的进气导叶)开度变为86％。
12.状态变化点12，即满负荷点：随后负荷继续增长至200MW并稳定，主燃流量增长至0.9869kg/s，旁路阀开度降至20％，机组稳定运行，这一状态变化点为满负荷点，记为状态变化点12；
工作状态描述：暖机结束后，机组负荷增长至200MW，旁路阀开度降至20％，机组完全进入了全负荷稳定运行状态。
将上述状态变化点进行整理，整理出启机过程中共12个状态变化点的各个参数，启机过程12个状态变化点的参数如表1所示。
表1启机过程中各个状态变化点的运行参数
步骤4：根据各个状态变化点不同的运行参数，结合本具体实施例中步骤2建立的化学反应器网络模型，利用开源程序包Cantera和Python编程软件计算出各个状态变化点的NOx排放值，将各个状态变化点的NOx 排放值连成曲线，与先前的实验结果进行比较，验证这种方法的可行性并建立NOx预估模型。
本领域普通人员可以理解上述NOx预估模型是基于开源程序包 Cantera和Python的数值计算；
具体步骤如下：
步骤4-1：基于步骤2所建立的化学反应器网络模型，利用开源程序包 Cantera和Python编程软件编写化学反应器网络模型中每个反应器的 Cantera程序段；
步骤4-2：在每个反应器的Cantera程序中，输入反应条件，包括：温度、压力、燃料流量、空气流量、反应器内物种组分、反应器体积和停留时间；
具体实施方式中利用步骤4-1的Cantera程序段，在程序段对应位置中输入表1中12个状态变化点所对应的参数。燃烧采用的燃料为焦炉煤气，在气体头部程序段燃料组分处输入焦炉煤气的组分；气体头部的计算结果为下一段燃烧的输入结果，根据图5的区域划分图，可以计算出每个燃烧区域的体积；
其中，图5中各编号为：1-混合区、2-外火焰区、3-中心火焰区、4-内火焰区、5-二次燃烧区、6-二次掺混区、7-火焰锋面区、8-后火焰区、9-火焰过渡区、10-二次火焰锋面区；
通过燃烧区的长度和图4的速度分布图计算出该区域的停留时间；燃烧室内属于恒压燃烧，同一状态变化点，不同燃烧区域的压力都是恒定的。
步骤4-3：在Cantera程序包中选择适合该重型燃气轮机燃烧室内燃烧的反应机理；
焦炉煤气的主要燃料成分是CH4、C2H6、H2和CO等气体，其中碳氢化物燃烧起主要作用，因此在程序中选择Gri Mech3.0(甲烷反应机理)，根据不同的燃烧室类型，在程序段中输入甲烷反应机理不同的形式，为后续计算提供理论基础。
步骤4-4：根据步骤2的化学反应器网络模型，将各个反应器的Cantera 程序段利用Python软件将各个反应器的Cantera程序段建立连接，得到整个化学反应器模型的Cantera程序包；
其中，包括：
(1)在每个状态变化点的程序段中，利用Python编程软件，建立外部程序，将整个化学反应器网络模型程序中所用的参数均定义在一个 Python脚本文件中，将这个脚本文件运行后，导入每一段程序中；
(2)根据化学反应网络模型，所建立的每段程序中，利用import语言，将脚本文件和上一程序段连接在一起，实现每个状态变化点内程序段连接，实现程序自动化。
步骤4-5：根据步骤3的各个工况，利用步骤4-4连接得到的Cantera 程序包计算出步骤3各个状态变化点的污染物排放值，将各个状态变化点污染物排放值数据进行整合，得出某重型燃气轮机启机过程中NOx排放预估模型；
其中，包括：
(1)检验步骤4-4连接的程序，将程序段中靠近出口的反应器程序段运行，最后生成scv文件中的各个物种的组分为该状态变化点燃烧后的产物；
(2)将1-12状态变化点的NOx排放量做成曲线，得到三种污染物排放的数学模型与西安热工研究院实验得到的数值相比较。
总结：如图11所示，为重型燃气轮机启机过程NOx排放量变化与现有实验值对比图，二者数值差距不大且各种污染物启机过程中变化趋势相同，本次采用的数学计算方法可行。
在本说明书的描述中，术语“一个实施例”、“一些实施例”、“具体实施例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或实例。而且，描述的具体特征、结构、材料或特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种燃气轮机启机过程NOx排放性能建模和软件开发设计方法，其特征在于，所述方法包括：
步骤1：选取目标为重型燃气轮机燃烧室并采集所述重型燃气轮机燃烧室选型设计参数的初始数据信息：
步骤2：根据与所述重型燃气轮机相近的温度分布云图和速度分布云图，将所述重型燃气轮机燃烧室划分为若干区域；根据若干所述区域的特点利用开源程序包Cantera建立燃烧室的化学反应器网络模型；
步骤3：根据重型燃气轮机启机过程中各个参数的变化，划分出若干个状态变化点；
步骤4：根据若干所述状态变化点不同的工况，并结合所述化学反应器网络模型，计算出各个所述状态变化点的NOx排放值，将各个所述状态变化点的NOx排放值连成曲线，建立NOx预估模型。
2.根据权利要求1所述的重型燃气轮机启机过程NOx排放性能建模和软件开发设计方法，其特征在于，所述步骤1中的燃烧室选型设计参数为重型燃气轮机燃烧室总体设计人员在燃烧室选型过程中会用到的选型设计参数，所述参数包括：
进口主燃料量、进口值班燃料量、进口空气流量、二次空气流量、压气机压比、燃烧室体积、燃烧室温度和燃烧筒个数。
3.根据权利要求2所述的重型燃气轮机启机过程NOx排放性能建模和软件开发设计方法，其特征在于，所述步骤2中的温度分布云图和速度分布云图通过CFD仿真计算流体力学软件建立。
4.根据权利要求3所述的重型燃气轮机启机过程NOx排放性能建模和软件开发设计方法，其特征在于，所述步骤3中的启机过程中各个参数包括：
机组转速、机组负荷、主燃料量、值班燃料量、二次空气流量、旁路阀开度和NOx排放值。
5.根据权利要求4所述的重型燃气轮机启机过程NOx排放性能建模和软件开发设计方法，其特征在于，所述步骤3中的状态变化点通过所述机组转速、所述机组负荷、所述主燃料量以及所述旁路阀开度的变化建立，且数量为12个；
其中，所述状态变化点至少包括：点火点、机组定速点、并网点、暖机负荷点以及满负荷点。
6.根据权利要求5所述的重型燃气轮机启机过程NOx排放性能建模和软件开发设计方法，其特征在于，所示步骤4具体包括如下步骤：
步骤4-1：根据所述化学反应器模型建立对应的Cantera模型，调用各个区域所对应反应器的Cantera程序包；
步骤4-2：在每个反应器的Cantera程序包中，输入反应条件；
步骤4-3：在所述Cantera程序包中选择适配于所述重型燃气轮机燃烧室内燃烧的反应机理；
步骤4-4：根据所述化学反应器网络模型，将各个反应器的Cantera程序段建立连接，得到整个化学反应器模型的Cantera程序包；
步骤4-5：根据步骤3的各个工况，利用步骤4-4得到的Cantera程序包计算出若干所述状态变化点的NOx排放值，将若干所述状态变化点NOx排放值数据进行整合，得出重型燃气轮机启机过程中NOx排放预估模型。
7.根据权利要求6所述的重型燃气轮机启机过程NOx排放性能建模和软件开发设计方法，其特征在于，所述步骤4-2中的反应条件包括：
温度、压力、燃料流量、空气流量、反应器内物种组分、反应器体积和停留时间。
8.根据权利要求6所述的重型燃气轮机启机过程NOx排放性能建模和软件开发设计方法，其特征在于，所述步骤4-4中的各个反应器的Cantera程序段通过Python软件建立连接。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110570553.0A CN113239493B (zh) | 2021-05-25 | 2021-05-25 | 重型燃气轮机启机过程NOx排放性能建模和软件开发设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110570553.0A CN113239493B (zh) | 2021-05-25 | 2021-05-25 | 重型燃气轮机启机过程NOx排放性能建模和软件开发设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113239493A true CN113239493A (zh) | 2021-08-10 |
| CN113239493B CN113239493B (zh) | 2023-12-26 |
Family
ID=77138693
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110570553.0A Active CN113239493B (zh) | 2021-05-25 | 2021-05-25 | 重型燃气轮机启机过程NOx排放性能建模和软件开发设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113239493B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115292882A (zh) * | 2022-06-26 | 2022-11-04 | 哈尔滨工程大学 | 一种基于化学反应器网络法的燃烧室污染物排放预测方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1930569A1 (en) * | 2006-11-01 | 2008-06-11 | ALSTOM Technology Ltd | System for controlling a combustion process for a gas turbine |
| US20190063332A1 (en) * | 2017-08-22 | 2019-02-28 | General Electric Company | Systems and methods for nox prediction in a power plant |
| CN110765698A (zh) * | 2019-10-15 | 2020-02-07 | 哈尔滨工程大学 | 一种燃气轮机燃烧室变工况排放性能预测方法 |
| CN111046577A (zh) * | 2019-12-24 | 2020-04-21 | 王国峰 | 一种重型燃气轮机污染物排放预估模型设计方法 |
-
2021
- 2021-05-25 CN CN202110570553.0A patent/CN113239493B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1930569A1 (en) * | 2006-11-01 | 2008-06-11 | ALSTOM Technology Ltd | System for controlling a combustion process for a gas turbine |
| US20190063332A1 (en) * | 2017-08-22 | 2019-02-28 | General Electric Company | Systems and methods for nox prediction in a power plant |
| CN110765698A (zh) * | 2019-10-15 | 2020-02-07 | 哈尔滨工程大学 | 一种燃气轮机燃烧室变工况排放性能预测方法 |
| CN111046577A (zh) * | 2019-12-24 | 2020-04-21 | 王国峰 | 一种重型燃气轮机污染物排放预估模型设计方法 |
Non-Patent Citations (2)
| Title |
|---|
| 乔征: "燃气轮机燃烧室NOx排放的数值研究", 中国优秀硕士学位论文全文数据库, pages 1 - 63 * |
| 赵巧男等: "轻型燃气轮机燃烧室NOx排放估算", 航空发动机, vol. 46, no. 1, pages 49 - 53 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115292882A (zh) * | 2022-06-26 | 2022-11-04 | 哈尔滨工程大学 | 一种基于化学反应器网络法的燃烧室污染物排放预测方法及系统 |
| CN115292882B (zh) * | 2022-06-26 | 2023-04-07 | 哈尔滨工程大学 | 一种基于化学反应器网络法的燃烧室污染物排放预测方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113239493B (zh) | 2023-12-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2608533C2 (ru) | Газотурбинная система, способ изменения выходной мощности газотурбинной системы, способ расширения диапазона регулирования газотурбинной системы, способ и система для повышения эффективности газовой турбины | |
| Chacartegui et al. | Real time simulation of medium size gas turbines | |
| Nielsen et al. | Modeling and validation of the thermal effects on gas turbine transients | |
| CN113239493B (zh) | 重型燃气轮机启机过程NOx排放性能建模和软件开发设计方法 | |
| Ntonas et al. | Integrated simulation framework for assessing turbocharger fault effects on diesel-engine performance and operability | |
| Hao et al. | Off-design performance of 9F gas turbine based on gPROMs and BP neural network model | |
| Pourhedayat et al. | An improved semi-analytical model for evaluating performance of gas turbine power plants | |
| Cameretti et al. | Cycle optimization and combustion analysis in a low-NO x micro-gas turbine | |
| Omar et al. | Modelling of Dry-Low Emission Gas Turbine Fuel System using First Principle Data-Driven Method. | |
| Long et al. | Research on active modulation of gas turbine cooling air flow | |
| Al-Mamoori et al. | Evaluation of gas fuel and biofuel usage in turbine | |
| Singh et al. | Transient performance analysis of an industrial gas turbine operating on low-calorific fuels | |
| Pan et al. | Off-design performance of gas turbine power units with alternative load-control strategies | |
| Petkovic et al. | Modelling the transient behaviour of gas turbines | |
| Bozza et al. | Adapting the micro-gas turbine operation to variable thermal and electrical requirements | |
| Iliescu et al. | Gas turbine modeling for load-frequency control | |
| Messina et al. | Aergia number: A new non‐dimensional group for gas turbine power estimation | |
| Wang et al. | Research on the Numerical Method of Pollutant Emission during Gas Turbine Start-up | |
| Stuen | Influence of Hydrogen Use As a Fuel on Aeroderivative Gas Turbine Performance | |
| Bozza et al. | Adapting the micro-gas turbine operation to variable thermal and electrical requirements | |
| Boksteen et al. | Performance modeling of a modern gas turbine for dispatch optimization | |
| Crosa et al. | Heavy-duty gas turbine plant aerothermodynamic simulation using simulink | |
| Cameretti et al. | A cfd based off-design study of micro-gas turbines combustors | |
| Mantelli et al. | Surge Prevention Techniques for a Turbocharged Solid Oxide Fuel Cell Hybrid System | |
| Nicolosi | Bio-fueled micro gas turbine: fuel flexibility and management optimization |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |