1.本发明涉及一种燃气轮机燃料控制方法，属于燃气轮机技术领域。
背景技术：2.当燃气轮机中燃料燃烧为扩散燃烧时，燃烧温度可以达到很高的理论燃烧温度，高于热nox的起始生成温度1650℃，在燃烧过程中极易产生数量较多的热nox污染物。因此为了实现燃气轮机的低nox燃烧，则需要改进燃料的燃烧方式。干式低nox燃烧技术(dln)采用均相贫预混的湍流火焰传播燃烧方法，通过将燃烧器改为两级串联式，从而将燃料燃烧由最初的扩散燃烧通过模式切换和改变燃料分配的方式逐步过渡至预混燃烧，通过调节燃/空比，可以保证燃烧温度始终低于1 650℃，故能有效抑制热nox的生成。
3.如何计算不同工况下燃烧所取燃料量以及合理分配一区、二区燃料量，而不对燃气轮机正常安全稳定运行产生影响，如导致回火、脱火、振荡燃烧、负荷波动、亦或是co、nox排量高等问题，燃料控制算法的设计至关重要。
技术实现要素：4.本发明的目的在于克服现有技术中的不足，提供一种燃气轮机燃料控制方法，能够准确计算燃气轮机在不同工况下燃烧所取燃料量以及合理分配一区、二区燃料量。
5.为达到上述目的，本发明是采用下述技术方案实现的：
6.第一方面，本发明提供了一种燃气轮机燃料控制方法，包括：
7.当燃料阀初次开启时，在给燃料阀分配流量的同时额外增加预充流量；
8.获取燃料行程基准、燃料分配比例、燃料热值并计算燃料阀分配流量；
9.获取压力机排气绝对压力、燃料温度并计算燃料阀预充流量；
10.将燃料阀分配流量和预充流量相加获取燃料阀的燃料质量流量；
11.获取燃料的摩尔质量、绝热系数比、温度、压缩因子，燃料阀的几何系数、阀门阻塞流压差系数、燃料阀间的绝对压力并结合燃料质量流量计算燃料阀流量系数；
12.根据燃料阀流量系数采用分段线性插值计算获取燃料阀开度；
13.其中，所述燃料阀包括二级燃料阀和切换燃料阀。
14.可选的，所述二级燃料阀和切换燃料阀的分配流量wf为：
[0015][0016]
式中，f为燃料行程基准，为燃料分配比例，h
l
为燃料热值。
[0017]
可选的，所述二级燃料阀的预充流量w0为：
[0018][0019]
所述切换燃料阀的预充流量w1为：
[0020]
[0021]
式中，p2为压力机排气绝对压力，t2为燃料温度。
[0022]
可选的，所述燃料阀：
[0023]
当预充开始时，将预充流量按照预设第一速率增加到目标值；
[0024]
当预充完成后，按照预设时间延迟后，将预充流量按照预设第二速率减小到0。
[0025]
可选的，所述燃料阀阀门流量系数c为：
[0026][0027]
式中，w为燃料质量流量，f
p
为燃料阀的几何系数，p2为燃料阀间的绝对压力，m为燃料摩尔质量，f
γ
为燃料绝热系数比，x
t
为燃料阀的阀门阻塞流压差系数，t2为燃料温度，z为燃料压缩因子。
[0028]
可选的，所述绝热系数比f
γ
为：
[0029][0030]
式中，γ为燃料的等熵指数；
[0031]
所述燃料压缩因子z为：
[0032]
z＝1-(β1×
t
22-β2×
t2+β3)
×
p2[0033]
式中，β1、β2、β3均为拟合系数；
[0034]
所述燃料阀的几何系数f
p
为：
[0035]fp
＝1+ε1×c′2+ε2×c′
[0036]
式中，c
′
为上一周期燃料阀阀门流量系数，ε1、ε2为燃料阀阀门特性系数。
[0037]
可选的，所述燃料阀还包括一级燃料阀，所述一级燃料阀的分配流量与二级燃料阀和切换燃料阀的分配流量计算方法相同，且所述一级燃料阀无预充流量；
[0038]
将一级燃料阀分配流量作为一级燃料阀的燃料质量流量；
[0039]
获取燃料的摩尔质量、绝热系数比、温度、压缩因子，燃料阀的几何系数、阀门阻塞流压差系数、燃料阀间的绝对压力并结合一级燃料阀的燃料质量流量计算一级燃料阀流量系数；
[0040]
根据一级燃料阀流量系数采用分段线性插值计算获取一级燃料阀开度。
[0041]
第二方面，本发明提供了一种燃气轮机燃料控制装置，所述装置执行时实现上述方法的步骤。
[0042]
与现有技术相比，本发明所达到的有益效果：
[0043]
本发明提供了一种燃气轮机燃料控制方法及装置，能够准确计算燃气轮机在不同工况下燃烧所取燃料量以及合理分配一区、二区燃料量，从而使燃料燃烧由最初的扩散燃烧顺利过渡至预混燃烧，有效抑制热nox的生成，并保证燃气轮机正常安全经济稳定运行。
附图说明
[0044]
图1是本发明实施例提供的一种燃气轮机燃料控制方法的流程图；
[0045]
图2是本发明实施例提供的一级燃料阀阀门流量系数对比结果一示意图；
[0046]
图3是本发明实施例提供的一级燃料阀阀门流量系数对比结果二示意图；
[0047]
图4是本发明实施例提供的二级燃料阀阀门流量系数对比结果一示意图；
[0048]
图5是本发明实施例提供的二级燃料阀阀门流量系数对比结果二示意图；
[0049]
图6是本发明实施例提供的切换燃料阀阀门流量系数对比结果示意图。
具体实施方式
[0050]
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
[0051]
实施例一：
[0052]
如图1所示，本发明提供了一种燃气轮机燃料控制方法，包括以下步骤：
[0053]
1、当燃料阀初次开启时，在给燃料阀分配流量的同时额外增加预充流量；
[0054]
1.1、获取燃料行程基准、燃料分配比例、燃料热值并计算燃料阀分配流量；
[0055]
1.2、获取压力机排气绝对压力、燃料温度并计算燃料阀预充流量；二级燃料阀和切换燃料阀的分配流量wf为：
[0056][0057]
式中，f为燃料行程基准，为燃料分配比例，h
l
为燃料热值。数值2.326为热值换算系数；数值1690为燃烧器单位热负荷；277.7778为单元换算系数；数值0.454为流量换算系数。
[0058]
二级燃料阀的预充流量w0为：
[0059][0060]
切换燃料阀的预充流量w1为：
[0061][0062]
式中，p2为压力机排气绝对压力，t2为燃料温度；数值6为预充时间常数，即以w0的流量预充6秒可将管路充满燃料；数值17.0644为燃料相对分子质量；数值6.863、7.35为二级燃料阀和切换燃料阀的燃料管路的体积常数；数值10.72917为换算系数；6.896为压力单位换算系数。
[0063]
燃料阀当预充开始时，将预充流量按照预设第一速率(0.05％/s)增加到目标值；
[0064]
燃料阀当预充完成后，按照预设时间(3s)延迟后，将预充流量按照预设第二速率(0.02％/s)减小到0。
[0065]
2、将燃料阀分配流量和预充流量相加获取燃料阀的燃料质量流量。
[0066]
3、获取燃料的摩尔质量、绝热系数比、温度、压缩因子，燃料阀的几何系数、阀门阻塞流压差系数、燃料阀间的绝对压力并结合燃料质量流量计算燃料阀流量系数；
[0067]
燃料阀阀门流量系数c为：
[0068][0069]
式中，w为燃料质量流量，f
p
为燃料阀的几何系数，p2为燃料阀间的绝对压力，m为燃料摩尔质量，f
γ
为燃料绝热系数比，x
t
为燃料阀的阀门阻塞流压差系数，t2为燃料温度，z为
燃料压缩因子。
[0070]
绝热系数比f
γ
为：
[0071][0072]
式中，γ为燃料的等熵指数；数值1.40为空气的等熵指数；
[0073]
燃料压缩因子z为：
[0074]
z＝1-(β1×
t
22-β2×
t2+β3)
×
p2[0075]
式中，β1、β2、β3均为拟合系数，分别取1.0798
×
10-11
、9.409
×
10-9
、0.0002；
[0076]
燃料阀的几何系数f
p
为：
[0077]fp
＝1+ε1×c′2+ε2×c′
[0078]
式中，c
′
为上一周期燃料阀阀门流量系数，ε1、ε2为燃料阀阀门特性系数，阀门特性系数如表1所示；
[0079]
表1：
[0080]
表2各燃料阀的阀门特性系数ε1和ε2[0081][0082]
燃料阀的膨胀系数y为：
[0083][0084]
式中，x为阻塞压差系数，y一般取2/3；阀门开度与燃料阀的阀门阻塞流压差系数x
t
的关系如表2所示；
[0085]
表2：
[0086][0087]
4、根据燃料阀流量系数采用分段线性插值计算获取燃料阀开度；
[0088]
其中，燃料阀包括二级燃料阀和切换燃料阀。
[0089]
此外，燃料阀还包括一级燃料阀，一级燃料阀的分配流量与二级燃料阀和切换燃料阀的分配流量计算方法相同，且一级燃料阀无预充流量；一级燃料阀开度的获取与上述燃料阀开度获取方法相同：
[0090]
将一级燃料阀分配流量作为一级燃料阀的燃料质量流量；获取燃料的摩尔质量、绝热系数比、温度、压缩因子，燃料阀的几何系数、阀门阻塞流压差系数、燃料阀间的绝对压力并结合一级燃料阀的燃料质量流量计算一级燃料阀流量系数；根据一级燃料阀流量系数采用分段线性插值计算获取一级燃料阀开度。
[0091]
燃料阀开度与燃料阀流量系数的关系如表3所示；
[0092]
表3：
[0093][0094][0095]
如图2-6所示，为本发明实施例提供的一种燃气轮机燃料控制方法与原系统燃料控制算法的对比结果；
[0096]
图2和图3分别为燃气轮机从燃气轮机点火到满负荷之间期间以及燃气轮机正常停机惰走期间一级燃料阀阀门流量系数对比结果，由于一级燃料阀在整个燃气轮机启停过程中均发生动作，因而时间跨度较长，故数据采样时间取2s；
[0097]
图4和图5分别为燃气轮机从贫贫模式(二级燃料阀为初次开启)至满负荷期间以及燃气轮机正常停机惰走期间二级燃料阀阀门流量系数对比结果，由于时间跨度相对较短，故数据采样时间取1s；
[0098]
图6为切换燃料阀阀门流量系数对比结果，由于切换燃料阀只在从贫贫模式切换至预混模式(切换燃料阀为初次开启)期间动作，其动作时间大概为1min，故采样时间取100ms。通过对比可以看出基于本发明实施例提供的一种燃气轮机燃料控制方法与原系统燃料控制算法计算结果完全吻合。
[0099]
实施例二：
[0100]
本发明实施例提供了一种燃气轮机燃料控制装置，所述装置执行时实现实施例一中方法的步骤。所述方法步骤为：
[0101]
1、当燃料阀初次开启时，在给燃料阀分配流量的同时额外增加预充流量；
[0102]
1.1、获取燃料行程基准、燃料分配比例、燃料热值并计算燃料阀分配流量；
[0103]
1.2、获取压力机排气绝对压力、燃料温度并计算燃料阀预充流量；
[0104]
燃料阀当预充开始时，将预充流量按照预设第一速率(0.05％/s)增加到目标值；
[0105]
燃料阀当预充完成后，按照预设时间(3s)延迟后，将预充流量按照预设第二速率(0.02％/s)减小到0。
[0106]
2、将燃料阀分配流量和预充流量相加获取燃料阀的燃料质量流量；
[0107]
3、获取燃料的摩尔质量、绝热系数比、温度、压缩因子，燃料阀的几何系数、阀门阻塞流压差系数、燃料阀间的绝对压力并结合燃料质量流量计算燃料阀流量系数；
[0108]
4、根据燃料阀流量系数采用分段线性插值计算获取燃料阀开度；
[0109]
其中，燃料阀包括二级燃料阀和切换燃料阀。
[0110]
此外，燃料阀还包括一级燃料阀，一级燃料阀的分配流量与二级燃料阀和切换燃料阀的分配流量相同，且一级燃料阀无预充流量；一级燃料阀开度的获取与上述燃料阀开度获取相同：
[0111]
将一级燃料阀分配流量作为一级燃料阀的燃料质量流量；获取燃料的摩尔质量、绝热系数比、温度、压缩因子，燃料阀的几何系数、阀门阻塞流压差系数、燃料阀间的绝对压力并结合一级燃料阀的燃料质量流量计算一级燃料阀流量系数；根据一级燃料阀流量系数采用分段线性插值计算获取一级燃料阀开度。
[0112]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0113]
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。
技术特征：式中，β1、β2、β3均为拟合系数；所述燃料阀的几何系数f
p
为：f
p
＝1+ε1×
c
′2+ε2×
c
′
式中，c
′
为上一周期燃料阀阀门流量系数，ε1、ε2为燃料阀阀门特性系数。7.根据权利要求1所述的一种燃气轮机燃料控制方法，其特征在于，所述燃料阀还包括一级燃料阀，所述一级燃料阀的分配流量与二级燃料阀和切换燃料阀的分配流量计算方法相同，且所述一级燃料阀无预充流量；将一级燃料阀分配流量作为一级燃料阀的燃料质量流量；获取燃料的摩尔质量、绝热系数比、温度、压缩因子，燃料阀的几何系数、阀门阻塞流压差系数、燃料阀间的绝对压力并结合一级燃料阀的燃料质量流量计算一级燃料阀流量系数；根据一级燃料阀流量系数采用分段线性插值计算获取一级燃料阀开度。8.一种燃气轮机燃料控制装置，其特征在于，所述装置执行时实现权利要求1-7任一项所述方法的步骤。
技术总结本发明公开了一种燃气轮机燃料控制方法，包括：当燃料阀初次开启时，在给燃料阀分配流量的同时额外增加预充流量；获取燃料行程基准、燃料分配比例、燃料热值并计算燃料阀分配流量；获取压力机排气绝对压力、燃料温度并计算燃料阀预充流量；将燃料阀分配流量和预充流量相加获取燃料阀的燃料质量流量；获取燃料的摩尔质量、绝热系数比、温度、压缩因子，燃料阀的几何系数、阀门阻塞流压差系数、燃料阀间的绝对压力并结合燃料质量流量计算燃料阀流量系数；根据燃料阀流量系数采用分段线性插值计算获取燃料阀开度；其中，燃料阀包括二级燃料阀和切换燃料阀；本发明能够准确计算燃气轮机在不同工况下燃烧所取燃料量以及合理分配一区、二区燃料量。二区燃料量。二区燃料量。
技术研发人员：陈子静 王勇 曹伟平 胡程斌 黄月丽 王新保 熊超 肖艳藏 刘云
受保护的技术使用者：南京国电南自维美德自动化有限公司
技术研发日：2022.05.10
技术公布日：2022/7/5
