CN106874569A - 一种变几何分轴式燃气轮机实时仿真建模方法 - Google Patents

Abstract

本发明涉及一种变几何分轴式燃气轮机实时仿真建模方法，该方法针对的燃气轮机结构特点为压气机导叶可调且是分轴式，由燃气发生器与动力涡轮两部分组成。该方法包括以下步骤：样本点选取、训练样本点获取、利用获取的样本点训练RBF神经网络、利用部件建模方法建立动力涡轮计算模块，将该模块与RBF神经网络相连，RBF神经网络模块负责输出燃气发生器相关参数，动力涡轮计算模块负责输出功率等相关参数，如此可得燃气轮机实时仿真模型。本发明的实时仿真模型的建模方法，该方法考虑了压气机进口导叶角度变化的情况，合理地简化了模型，结合RBF神经网络与部件法二者的优点，在保证模型高实时性和高精度的同时，大大缩小样本规模。

Description

一种变几何分轴式燃气轮机实时仿真建模方法

技术领域

本发明涉及燃气轮机领域，更具体地，涉及一种变几何分轴式燃气轮机实时仿真的建模方法。

背景技术

在进行燃气轮机半物理仿真实验以及燃气轮机控制规律研究过程中，对燃气轮机仿真模型的实时性要求较高。特别是对于变几何的分轴式燃气轮机，该类燃气轮机包括带可调导叶压气机的燃气发生器以及动力涡轮，由于变量较多，基于燃气轮机工作机理建立的非线性模型实时性较差。线性化建模方法虽然有较好的实时性，但对于偏离稳态工况较远的情况仿真精度难以保证。基于神经网络的系统辨识法虽然能达到较高的精度和实时性，但对于变几何的分轴式燃气轮机，由于动态工况繁多，所需样本规模庞大，难以直接建立仿真模型。

发明内容

本发明为解决以上现有技术的缺陷和不足，提供了一种基于RBF神经网络与部件法的复合建模方法。基于燃气轮机工作机理建立起来的部件法具有高精度的优点，但实时性较差，而利用RBF神经网络建立的燃气轮机仿真模型虽然具有优异的实时性，但对于变几何分轴式燃气轮机面临着变量多、样本规模庞大的问题，本发明的变几何分轴式燃气轮机实时仿真建模方法，结合现有RBF神经网络法与部件法二者的优势，利用动力涡轮相对转速NPT对燃气发生器影响很小的特点，用于训练RBF神经网络的样本点中仅选取了两个不同的动力涡轮相对转速NPT值，大大减小了样本规模，在此基础上， 通过部件法单独建立的动力涡轮计算模块模拟动力涡轮工况，保证了样本规模减小的同时模型的计算精度与实时性。

为实现以上发明目的，本发明采用的技术方案是：

一种变几何分轴式燃气轮机实时仿真建模方法，所述变几何分轴式燃气轮机包括可调导叶压气机、燃气发生器和动力涡轮，其特征在于，所述实时仿真建模方法包括以下步骤：

SS1.样本点选取。

对于变几何分轴式燃气轮机，可调导叶压气机的导叶角度VGV、燃油流量Wf，燃气发生器相对转速NH以及动力涡轮相对转速NPT四个变量可以唯一确定燃气轮机某一动态工况。为确定各部件工作状况，可调导叶压气机需要确定转速、流量、压比中的任意两个参数以及进口导叶角度共计三个参数，燃气涡轮与动力涡轮均需确定转速、膨胀比、流量三个参数中任意两个变量，燃烧室部件则需确定燃油流量的大小才可唯一确定各部件工作状态，综上，如果确定压气机流量WaC、导叶角度VGV、燃气发生器相对转速NH、燃油流量Wf、燃气涡轮膨胀比π_PT、动力涡轮相对转速NPT、动力涡轮膨胀比π_PT7个变量，即可确定燃气轮机的工作状态。燃气轮机工作时，各部件并不是相互孤立的，有机械和气动间的联系，有如下三条约束关系：

Wa_C＝Wa_GT-W_f (1)

Wa_C＝Wa_PT-W_f (2)

π_C＝π_GT×π_PT (3)

上面三个约束关系消除了三个变量，最终剩余导叶角度VGV、燃油流量Wf、燃气发生器相对转速NH、动力涡轮相对转速NPT四个变量唯一确定燃机的工作状况。

在燃气轮机工作时，对于导叶角度VGV、燃油流量Wf、燃气发生器相对转速NH三个参数，在其各自变化范围内，分别均匀地选取a个、b个、c个值(a、b、c均为大于等于1的整数)。由于动力涡轮相对转速NPT对动力涡轮流通能力影响较小，故对燃气发生器工况的影响很小。需要注意的是，由于气动关系，动力涡轮相对转速NPT对动力涡轮本身的效率以及输出功率还是有较大影响的。对于影响燃气发生器工况较小的参数，稀疏的取值即可满足神经网络精确地模拟燃气发生器动态工况，而过多的取值只会增加样本规模、加大工作量，因此为减小样本规模，动力涡轮相对转速NPT仅取一大一小两个值，具体取值视具体仿真情况而定。如，仅需模拟动力涡轮转速在NPT₁到NPT₂间的变化时，为防止燃机工况超出样本规模，NPT取NPT₁和NPT₂两个值即可，综上，共选取样本点a×b×c×2组。

SS2.样本点获取。

每一组样本点由输入参数和输出参数组成。导叶角度VGV、燃油流量Wf、燃气发生器相对转速NH以及动力涡轮相对转速NPT四个参数构成输入参数，其余燃气发生器工作参数为输出参数，具体的，输出参数包括：压气机压比、压气机效率、进气流量、燃气涡轮膨胀比、燃气涡轮效率、动力涡轮效率、动力涡轮膨胀比、压气机出口总压、压气机出口总温、燃烧室出口总温、燃烧室出口总压、燃气涡轮出口总压、燃气涡轮出口总温、动力涡轮出口总温、动力涡轮出口总压、动力涡轮输出功率、燃气发生器角加速度，通过将输入参数输入到基于燃气轮机工作机理的燃气轮机总体计算商业软件或已有的通过求解非线性方程组的方式建立起来的燃气轮机总体计算仿真模型中，分别计算出每组样本点的输出参数，从而获取完整的训练样本。

SS3.训练RBF神经网络。

利用现有的RBF神经网络算法，例如MATLAB/SIMULINK内置的RBF 神经网络工具箱，以输出参数作为目标值，用a×b×c×2组样本点训练神经网络，神经网络在这里起到了一种多维拟合的作用，利用它建立起输入参数与输出参数间的对应关系，例如输出功率与输入参数的对应关系可表示为：PWSD＝f(NH,NPT,Wf,VGV)。

SS4.建立实时仿真模型。

通过RBF神经网络算法模拟燃气发生器，部件法建立的动力涡轮仿真模块模拟动力涡轮，燃气轮机的实时仿真模型得以建立。RBF神经网络输入量为VGV、Wf、NH以及NPT，其中VGV、Wf以及NPT为模型的给定输入量，NH则是通过模型计算得到。为减小输入数据的分布范围，加快神经网络的收敛速度，减小训练时间，将输入参数归一化到[-1,1]的区间范围内。将归一化后的参数输入到训练好的RBF神经网络中，计算出燃气发生器各工作参数后经逆归一化过程，将参数映射回归一化前的状态。RBF神经网络计算的燃气发生器工作参数包括当前时刻燃气发生器的转速加速度通过对时间积分即可求出每一时刻的NH值，离散化后的计算公式： 其中，Δt为设置的时间步长，时间步长越小，计算越精确，运算量越大。RBF神经网络算出的燃气发生器的出口参数，包括油气比、燃气发生器出口总压、总温与工质流量，将这些参数加上模型给定输入量NPT作为输入参数输入到基于部件法建立起来的动力涡轮计算模块中去，从而完成动力涡轮部分的仿真计算。

采用两种方法分别计算燃气发生器与动力涡轮，其目的是：上述训练好的RBF神经网络可以预测不同VGV、Wf、NH、NPT所对应的燃气轮机燃气发生器各工作参数。可以看出，由于样本点中仅取了两个NPT的值，对于受NPT影响极小的燃气发生器，神经网络能准确的模拟，但对于受NPT影 响较大的动力涡轮本身来说，较少的NPT取值会直接影响动力涡轮的计算精度，为此，NPT对工况的影响将通过部件法建立起来的动力涡轮计算模块单独实现。该动力涡轮计算模块利用燃气热力性质计算公式以及已有的部件特性完成动力涡轮部分的仿真。具体的，由输入量NPT和动力涡轮膨胀比对部件特性插值(可以是常用的线性插值)得到涡轮效率以及工质流量，计算出若工质可逆绝膨胀，动力涡轮的出口总温和焓降。进而由动力涡轮效率算出动力涡轮实际输出功率与出口总温。

需要注意，NH的初始值应为稳态工况下VGV、Wf以及NPT的初始值所对应的高压轴转速。可以通过传统的燃气轮机仿真模型计算出VGV、Wf以及NPT的初始值所对应的燃气轮机稳态工况下NH的取值。如此，给定VGV、Wf以及NPT即可运行仿真模型。

SS5.模型验证。具体的，通过给定一段VGV、Wf以及NPT的变化规律，对比已有的通过求解非线性方程组的方式建立起来的燃气轮机总体计算仿真模型与本专利建立起来的实时仿真模型各参数计算结果，即可验证模型的准确性。

基于燃气轮机工作机理建立起来的部件法具有高精度的优点，但实时性较差。利用RBF神经网络建立的燃气轮机仿真模型具有优异的实时性，但对于变几何分轴式燃气轮机面临着变量多、样本规模庞大的问题。本发明提供的变几何分轴式燃气轮机实时仿真模型的建模方法，结合了两种建模方法的优点，通过部件法与神经网络分别计算燃气发生器部分与动力涡轮部分得到了如下效果：

1、由于NPT对燃气发生器影响较小，但对动力涡轮本身影响较大。两种方法分别计算，使得燃气发生器部分由原来的VGV、Wf、NH、NPT四变量问题简化为主要与VGV、Wf、NH相关的准三变量问题，如此，用于训练RBF神经网络的样本规模便大大减小了，简化了模型的建立。

2、通过部件法单独计算动力涡轮模块，考虑了NPT对动力涡轮工况的影响，保证了模型的仿真精度。

3、神经网络计算了燃气发生器部分，避免了部件法中非线性方程组的求解以及绝大部分的工质热力计算，大大提高了模型的运算速度与实时性。

与现有技术相比，本发明的有益效果是：

本发明提供了变几何分轴式燃气轮机实时仿真模型的建模方法，该方法考虑了压气机进口导叶角度变化的情况，合理地简化了模型，结合RBF神经网络与部件法二者的优点，在保证模型高实时性和高精度的同时，大大缩小样本规模，减小了工作量。

附图说明

图1为本发明提供的方法的流程图。

图2为本发明建立的仿真模型结构图。

具体实施方式

为使本发明的目的、技术方案及优点更加清楚明白，以下参照附图并举实施例，对本发明进一步详细说明，以使本发明的优点和特征能更易于被本领域技术人员理解。需要说明的是，以下所述仅为本发明的较佳实施例，并不因此而限定本发明的保护范围。实际上，在未背离本发明的范围或精神的情况下，可以在本发明中进行各种修改和变化，这对本领域技术人员来说将是显而易见的。例如，作为一个实施例的一部分示出或描述的特征可以与另一个实施例一起使用来产生又一个实施例。因此，意图是本发明将这样的修改和变化包括在所附的权利要求书和它们的等同物的范围内。

变几何分轴式燃气轮机的结构特点为压气机导叶可调且是分轴式结构，由燃气发生器与动力涡轮两部分组成。如图1所示，本发明提供的变几何分 轴式燃气轮机实时仿真建模方法包括以下步骤：

SS1.样本点选取。本例仅描述该燃气轮机慢车以上转速的建模过程。通过传统的燃气轮机仿真模型，获得不同燃气发生器转速下压气机进口导叶角度VGV的变化范围：f₁(NH)≤VGV≤f₂(NH)，其中VGV的上限f₂(NH)为燃气发生器在该转速下，发动机刚好处于喘振边界处。VGV在其变化范围内取5个值，Wf可按等差数列在额定燃油流量范围内选取8个值。NH在0.75到1转速范围内按等差数列取8个值。NPT取0.4和1.1两个值。综上，样本点规模为5×8×8×2＝640组点。一组样本点包括输入参数和输出参数。输入参数为Wf、VGV、NH以及NPT共计四个参数，输出参数为燃气发生器各截面以及相关工作参数。

SS2.样本点获取。通过实验数据或商业软件获取燃气轮机工况样本点。基于传统的部件法，通过求解非线性方程组，建立传统的燃气轮机仿真模型。将选取好的640组样本点输入到建好的仿真模型中，算出640组燃气轮机动态工况参数，包括燃气发生器角加速度、压气机流量、压气机出口总温、燃烧室出口总温、燃气涡轮出口总温、燃气涡轮出口总压、压气机压比、压气机效率、燃气涡轮膨胀比、燃气涡轮效率共计10个参数。

SS3.利用获取的样本点训练RBF神经网络。本实例中，利用Matlab软件，以Wf、VGV、NH和NPT为输入参数，燃气发生器角加速度压压气机流量W_C、压气机出口总温T₂、燃烧室出口总温T₄、燃气涡轮出口总温T₆、燃气涡轮出口总压p₆、压气机压比π_C、压气机效率η_C、燃气涡轮膨胀比π_HT、燃气涡轮效率η_HT共计10个参数为目标量，通过matlab中的“premnmx”和“newrb”函数分别归一化样本数据和训练RBF神经网络。函数的详细用法参见matlab帮助文件。

SS4.如图2所示，建立实时仿真模型。神经网络输入参数为NPT、Wf、VGV和NH，先利用matlab中的“tramnmx”函数将输入参数归一化到[-1,1]区间内，将归一化后的数据作为输入量输入已训练好的RBF神经网络，从而算出燃气发生器各相关参数。之后利用函数“postmnmx”将神经网络算出的参数映射回归一化前的值。神经网络算出当前时刻燃气发生器的转速加速度 通过公式可算出下一时刻的NH值，如此迭代循环，可算出各时刻的NH，完成燃气发生器的仿真。其中，Δt为设置的时间步长，时间步长越小，计算越精确，运算量越大，本例中Δt取为0.2s。动力涡轮部分的仿真由利用部件建模方法建立的动力涡轮计算模块单独完成。动力涡轮模块的输入参数为NPT以及由神经网络算出的燃气发生器出口总压P₆、总温T₆、比焓h₆以及压气机进气流量W_C，输出参数为动力涡轮效率η_PT、输出功率PWSD与排气温度T₈。利用RBF神经网络算出的工质流量与给定的NPT通过已有的动力涡轮特性数据点插值获得动力涡轮效率η_PT，动力涡轮膨胀比通过公式求得，π_PT、θ₆₇、θ₈₀、P₀分别代表动力涡轮膨胀比、燃气涡轮与动力涡轮间过渡段总压恢复系数、排气道总压恢复系数与大气压力。之后，利用工质热力计算关系式可算出若工质可逆绝热膨胀，动力涡轮的排气总温T'₈与比焓h'₈。通过公式h₈＝(h₆-h'₈)·η_PT+h₆可求出动力涡轮出口实际比焓h₈。根据燃气热力性质计算公式可由h₈求出动力涡轮排气总温T₈。则动力涡轮输出功率可由公式PWSD＝(h₆-h₈)·(W_C+Wf)求得，如此完成动力涡轮部分的建模过程。综上，RBF神经网络模块负责输出燃气发生器相关参数，动力涡轮计算模块负责计算动力涡轮部分，燃气轮机实时仿真模型得以建立。

SS5.将实时仿真模型与之前已建好的传统的仿真模型进行对比验证。向 实时仿真模型和传统仿真模型输入一段相同的输入参数Wf、NPT、VGV，对比两个模型的计算结果，两个模型在各参数上均吻合良好，模型的可靠性得到验证。

对于变几何的分轴式燃气轮机，只需确定燃油流量Wf、燃气发生器相对转速NH、动力涡轮相对转速NPT以及压气机进口导叶角度VGV四个参数即可唯一确定燃气轮机的某一动态工况。通过计算得到的燃气轮机动态工况样本点，利用RBF神经网络建立起这四个变量与燃气发生器其余各参数的对应关系，并据此建立燃气发生器部分的仿真模型。

NPT、VGV以及Wf是模型的输入量，由外部负载和控制系统决定，NH则需要经过计算得到。Δt为离散化后的时间步长，需注意Δt不宜取得过大，以免影响模型的仿真精度。本实例中取为0.2s。通过反复迭代可以模拟燃气发生器部分各时刻的动态工况。

显然，本发明的上述实例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。

Claims (5)

1.一种变几何分轴式燃气轮机实时仿真建模方法，所述变几何分轴式燃气轮机包括可调导叶压气机、燃气发生器和动力涡轮，其特征在于：所述实时仿真建模方法包括以下步骤：

SS1.样本点选取

选取可调导叶压气机的导叶角度VGV、燃油流量Wf、燃气发生器相对转速NH、动力涡轮相对转速NPT四个变量确定燃机的工作状况；

在燃气轮机工作时，对于可调导叶压气机的导叶角度VGV、燃油流量Wf、燃气发生器相对转速NH等三个参数，在其各自变化范围内，分别均匀地选取a个、b个、c个值，其中a、b、c均为大于等于1的整数；对于动力涡轮相对转速NPT，在其变化范围内，选取一大一小两个值；

综上，共选取样本点a×b×c×2组。

SS2.训练样本点获取

每一组训练样本点由输入参数和输出参数组成，其中导叶角度VGV、燃油流量Wf、燃气发生器相对转速NH以及动力涡轮相对转速NPT四个参数构成输入参数，其余燃气发生器工作参数为输出参数，

通过将输入参数输入到现有的通过求解非线性方程组的方式建立起来的燃气轮机总体计算仿真模型中，分别计算出每组样本点的输出参数，从而获取完整的训练样本。

SS3.利用获取的样本点训练RBF神经网络

利用现有的RBF神经网络算法，以输出参数作为目标值，用a×b×c×2组样本点训练RBF神经网络，利用RBF神经网络建立输入参数与输出参数间的对应关系。

SS4.建立实时仿真模型

利用部件建模方法建立动力涡轮计算模块模拟动力涡轮，并利用RBF神经网络模块模拟燃气发生器，将所述动力涡轮计算模块与RBF神经网络模块相连，RBF神经网络模块负责输出燃气发生器相关参数，动力涡轮计算模块负责输出功率等相关参数，如此可得燃气轮机实时仿真模型。

2.根据权利要求1所述的变几何分轴式燃气轮机实时仿真建模方法，其特征在于：由于动力涡轮相对转速NPT只对动力涡轮本身有较大影响，对燃气发生器影响很小，动力涡轮相对转速NPT在其变化范围内，选取一大一小两个值；动力涡轮相对转速NPT对动力涡轮本身的影响通过动力涡轮计算模块单独实现。

3.根据上述权利要求所述的变几何分轴式燃气轮机实时仿真建模方法，其特征在于：一组训练样本点的输入参数为供油量Wf、压气机进口导叶角度VGV、燃气发生器相对转速NH、以及动力涡轮相对转速NPT共计四个参数，输出参数为燃气发生器各截面参数以及各工作参数。

4.根据上述权利要求所述的变几何分轴式燃气轮机实时仿真建模方法，其特征在于：训练样本点的输出参数包括：压气机压比、压气机效率、进气流量、燃气涡轮膨胀比、燃气涡轮效率、动力涡轮效率、动力涡轮膨胀比、压气机出口总压、压气机出口总温、燃烧室出口总温、燃烧室出口总压、燃气涡轮出口总压、燃气涡轮出口总温、动力涡轮出口总温、动力涡轮出口总压、动力涡轮输出功率、燃气发生器角加速度。

5.根据上述权利要求所述的变几何分轴式燃气轮机实时仿真建模方法，其特征在于：所述实时仿真建模方法还包括模型验证步骤，通过给定一段VGV、Wf以及NPT的变化规律，对比已有的通过求解非线性方程组的方式建立起来的燃气轮机总体计算仿真模型与步骤SS4实时仿真模型各参数计算结果，即可验证模型的准确性。