当前位置：文档库 › 基于改进遗传算法的燃气轮机自适应建模仿真

基于改进遗传算法的燃气轮机自适应建模仿真

Dynamical Systems and Control 动力系统与控制, 2019, 8(2), 140-147

Published Online April 2019 in Hans. https://www.wendangku.net/doc/de14576993.html,/journal/dsc

https://https://www.wendangku.net/doc/de14576993.html,/10.12677/dsc.2019.82016

Adaptive Modeling and Simulation of Gas

Turbine Based on Improved Genetic

Algorithms

Jianfeng Zhao1, Qianchao Liang1, Huazhi Li2*

1College of Power Engineering, Naval University of Engineering, Wuhan Hubei

2College of Information and Communication, National University of Defense Technology, Wuhan Hubei

Received: Mar. 12th, 2019; accepted: Mar. 23rd, 2019; published: Apr. 4th, 2019

Abstract

In order to guarantee the accurate description and diagnosis of gas turbine performance and fault, the establishment of accurate gas turbine model plays a key role. Under variable working condi-tions, the accuracy of gas turbine model calculation depends on the accuracy of the characteristics of gas turbine components (compressor and turbine). In order to ensure the accuracy of gas tur-bine model characteristics, an improved genetic algorithm and adaptive technology are estab-lished in this paper to obtain gas turbine characteristics. Adaptive genetic algorithm (AGA) and gray scale prediction (GM) were combined to optimize the component characteristic correction factor to ensure the optimal gas turbine characteristics and improve the convergence speed of the genetic algorithm.

Keywords

Gas Turbine, Adaptive Genetic Algorithm, Gray Scale Prediction, Modeling and Simulation

基于改进遗传算法的燃气轮机自适应建模仿真

赵建锋1，梁前超1，李华志2*

1海军工程大学动力工程学院，湖北武汉

2国防科技大学信息通信学院，湖北武汉

收稿日期：2019年3月12日；录用日期：2019年3月23日；发布日期：2019年4月4日

*通讯作者。

赵建峰等

摘要

为了保证对燃气轮机的性能与故障进行精准的描述及诊断，建立精确的燃气轮机模型便起到关键作用。燃气轮机处于变工况条件下，其模型计算准确性取决于燃气轮机各部件(压气机与透平)特性的精度。为保证获得燃气轮机模型特性的精度，本文建立一种改进的遗传算法与自适应技术对燃气轮机的特性进行获取。采用自适应遗传算法(AGA)和灰度预测(GM)相结合的方法，对部件特向修正因子进行优化，保证在获得最优的燃气轮机特性的同时，能提高遗传算法的收敛速度。

关键词

燃气轮机，自适应遗传算法，灰度预测，建模仿真

This work is licensed under the Creative Commons Attribution International License (CC BY). https://www.wendangku.net/doc/de14576993.html,/licenses/by/4.0/

1. 引言

燃气轮机模型是燃气动力系统中的关键环节，因此燃气轮机模型的准确度直接影响动力系统仿真、故障诊断、检测控制等工作的质量与效率。同时燃气轮机是一种非线性强、时变性高、构成复杂的气动热力系统，因此如何建立高精度的燃机轮机数学模型是燃气动力系统仿真研究的热点问题[1]。

燃气轮机由压气机、燃烧室、涡轮、附属设备及机匣等组成，其中压气机、涡轮又可根据压力、气动结构等细分为高、低压压气机、高、低压涡轮和动力涡轮，而燃气轮机中的压气机、涡轮以及附属系统中的燃油系统、引气系统等均为强非线性、动态部件，因此自适应控制建模的关键是对以上部件特性参数进行合理的预测与校正。段守付等[2]通过改进的加权函数自适应建模方法对部件特性参数进行校正。谢光华等[3]运用改进的共轭梯度算法对部件特性参数进行校正。陈玉春等[4]根据实验测量参数合理选择部件特性修正因子并通过部件特性删除法求解部件特性修正因子。王永华等[5]通过在数学模型中引入压力比函数Ｚ和辅助参数β，对部件特性进行自适应修正。Li YG 等[6]采用遗传算法对燃气轮机进行了自适应建模，并研究自适应参数的敏感度和自适应的鲁棒性。刘永葆等[7]用测量基数、部件匹配、测量参数相关性、以及传感器测量精度等四个因素对部件特性修正，从而进一步改进了遗传算法，提高了算法的收敛精度与速度。

获得最优解的方法从传统的优化算法(如线性规划的单纯形法，非线性规划的基于梯度的各类迭代算法)改进到现代智能算法(如粒子群寻优算法、遗传算法)，不但在寻优方法上不断获得突破，同时在目标函数构造上也不断改进，通过设定合理的加权，使得目标函数更加合理[8]。

智能算法的引入使得能够实现参数更多的多目标优化，对初值不敏感，限制条件更加宽松，但是也带来了计算耗费资源较重，计算时间较长的缺点。因此如果能够合理地预测最优解的取值空间，对智能算法中的限制条件进一步进行限制，则能够缩短最优解的解算时间，同时也能一定程度提高最优解的精度。本文建立一种改进的遗传算法与自适应技术对燃气轮机的特性进行获取。采用自适应遗传算法(AGA)和灰度预测(GM)相结合的方法(本文设定为GM-AGA 算法)，对部件特向修正因子进行优化，保证在获得最优的燃气轮机特性的同时，能提高遗传算法的收敛速度。

Open Access

赵建峰等

2. 燃气轮机自适应控制

由于通过实验来获取燃气轮机部件特性耗费昂贵、周期长，因此不可能对每台燃气轮机进行特性实验，但是可以根据同一型号燃气轮机，其特性相近的情况，用已知改型燃气轮机装置的通用标准特性建立数学模型，在根据所研究的特定燃气轮机装置实际运行的有关测量参数，对模型中的部件特性进行修正，以达到获取精确的燃气轮机数学模型，这种方法就是燃气轮机部件的自适应建模[9]。

在数学模型中发动机目标参数Y 可通过函数用部件特性参数数据X 、输入参数Uin 及环境参数Uam 表示：

(),,Y f X Uin Uam = (1-1)

其中部件特性参数X 可由解析式获图标法获得。但由于部件特性参数是由通用特性数据获得，因此与实测数据存在偏差e 。用Xact 表示部件的真实特性，用Xcal 表示部件通过通用特性获得的特性，定义修正因子M ，则M 计i 算如式(2)所示：

Xact

M Xcal

(1-2) 将式(2)代入式(1)中可得修正后的数学模型；

(),,Y f X Uin Uam = (1-3)

根据数学模型的需求将需要修正的一系列部件特性参数的修正因子设为()1,2i M i n = ，为了获得最优的M i 解，需要通过构造特定的数学模，因此设发动机试车平台中能够实测的参数值为

()1,21i Ca i n = ，对应的模型计算出的参数值为()1,21i Cc i n = ，求解发动机模型平衡方程产生的

残量为(),1,22i err i i n == ，则建立函数F 如下：

()()2

n n i i i i

i F Ca Cc err =

==?+∑∑ (1-4)

能够使F 值最小的M i 解则为最优解。

获得最优解的方法

获得最优解的方法从传统的优化算法(如线性规划的单纯形法，非线性规划的基于梯度的各类迭代算法)改进到现代智能算法(如粒子群寻优算法、遗传算法)，不但在寻优方法上不断获得突破，同时在目标函数构造上也不断改进，通过设定合理的加权，使得目标函数更加合理[10]。

3. 最优算法的构造

3.1. 灰度预测

造成单个个体燃气轮机部件的特性与通用特性的差异的原因有多方面，例如生产中几何尺寸的偏差，使用环境与测试环境的偏差以及在使用过程中的零部件的损耗等。由于根本无法获得完整的信息能够保证将造成偏差的所有原因都进行量化并形成合理的修正系数，因此为了能够有效合理的预测部件特性的修正因子，本文引入简单()1,1GM 灰度预测模型[11]。

根据文献[12]，()1,1GM 模型定义为：

(

)

()()()()()()

011

1112

x k a x k x k u +=?

+++ (2-1) 其中，,a u 为待定系数，()0X 为原始数列：

赵建峰等

()(

)

()()()()(){

0001,2,,X x x x n = (2-2)

取2-2式中()

()0i x

i M =，对()0X 进行一次累加，得到新数列：

()()()()()()(){}

1111

1,2,,X x x x n = (2-3)

其中()

()()()101k i x

k x k ==∑。

构造矩阵：

()()()()()()()()

()()()()()()()

()()()

()()

110110(0)1111212212313,21112x x x x x x Y B x n x n x n ?+

?+ == ??+

(2-4) [],T

a u ?=为待辨识参数向量，则(2-1)可写成：

*Y B =? (2-5)

参数向量?可用最小二乘法求取，即

(

)

T T a u B B

B Y ? ?==

(2-6)

把求取的参数代入(2-1)式，并求出其近似时间响应式为：

()()()()111ak u u x k ax k e a a

?+=?+ (2-7) 累减还原为：

()()()()()()

()

()()01111111a

x k x k x k e x k x e ?+=+?=?? (2-8)

()()01i x k + 即为预估的部件修正因子 ()1,2,,i M

i n = 。取() 10%1i M ±+为遗传算法中的种群空间进行进一步优化求解。

3.2. 自适应遗传算法(AGA)

遗传算法(GA)是一类借鉴生物界的进化规律(适者生存，优胜劣汰遗传机制)演化而来的随机化搜索方法。它是由美国的J. Holland 教授1975年首先提出，其主要特点是直接对结构对象进行操作，不存在求导和函数连续性的限定；具有内在的隐并行性和更好的全局寻优能力；采用概率化的寻优方法，能自动获取和指导优化的搜索空间，自适应地调整搜索方向，不需要确定的规则。遗传算法的这些性质，已被人们广泛地应用于组合优化、机器学习、信号处理、自适应控制和人工生命等领域。

简单遗传算法局部搜索能力不强，容易陷入局部最优解，且进化运算的后期速度比较慢。因此可引入自适应机制[10]，形成AGA 算法，AGA 旨在通过以不同的方式实现搜索和随机性之间的权衡，根据适合度值自适应地改变Pc 和Pm 的值，当群体倾向于停留在局部最优时(也就是群体适应度集中，多样性比较差时) Pc 和Pm 的值增加，并且当群体在解空间中散布时(也就是群体适应度分散，多样性比较高时)减小。所以需要一个针对Pc 和Pm 的计算公式如下所示，让其符合动态变化。

赵建峰等

max 1avg

max avg

3avg 1,,f f K f f f f Pc K f f ′? ′≥ ?=

′< (2-9) max 2avg

max avg

avg 2,,f f K f f f f Pm K f f

′? ′≥ ?= ′< (2-10) f max ，群体中的最大适应值f max ，群体中的最大适应值，f avg ，群体平均适应值f avg ，群体平均适应值，

1f ′要交叉的两个个体中较大的适应值，2f ′要变异个体的适应值，K 1，K 2，K 3和K 4是常数。

至此本文中燃气轮机自适应建模流程如图1所示。

Figure 1. Flow chart of gas turbine adaptive modelling

图1. 燃气轮机自适应建模流程

4. 燃气轮机数学模型精确解

本文的研究对象为某型单轴带回热燃气轮机，该燃气轮机装置结构示意图如图2所示。图1中，各个截面用0~6这7个阿拉伯数字表示：0表示大气环境，1表示压气机进口截面，2表示压气机出口截面，3表示回热器冷端进口与热端出口截面，4表示回热器冷端出口热端进口截面，5表示燃烧室出口(涡轮进口)截面，6表示涡轮出口截面，对该型燃气轮机，其部件特性采用6个修正因子，见表1所示。

根据划分的7个截面中描述每个截面气动特性的参数需求，确定共需要9个测量参数，分别为：转子转速n 、压气机进口流量G 1、压气机出口温度T 1、压气机出口压力P 1、回热器冷端出口温度T 2、燃烧室出口温度T 3、涡轮出口温度T 4、涡轮出口压力P 2、回热器热端出口温度T 5 [13]。

3个输入参数分别为：环境温度T 0、环境压力P 0、燃油质量流量G 2；1个平衡检验误差参数为：轴功率平衡误差e 。

适应度F 的定义为[7]：

赵建峰等

F F =

+ (3-1)

Table 1. Modified factors 表1. 修正因子

部件名称修正因子

压气机 1cact

ccal M ηη= 2cact

ccal

G M G =

换热器 3eact

ecal

M ηη=

燃烧室 4bact

bcal

M ηη= 涡轮

5tact

tcal

M ηη=

MM 6=

G tt aaaatt

G tt aaaacc

Figure 2. Structure diagram of single-shaft regenerator gas turbine 图2. 单轴回热燃气轮机结构图

(3-1)式中c F 由式(1-4)获得，由于式(1-4)中的定义可知函数c F 为正数，随着M i 逐步优化，c F 逐步减少，适应度F 逐步增大，适应度极值等于1。

首先获得每个测量参数的多组测量数据，由测量数据计算每次测量的修正因子集()0X 作为原始数列，

根据本文2.1章节中所述的()1,1GM 模型获得预测的修正因子 ()1,2,i M i n = 。取() 10%1i M ±+为自适应遗传算法中的种群空间。

根据自适应建模原理，运用自适应遗传算法(AGA)编程．其中GM-AGA 算法种群空间设为

() 10%1i M

±+，AGA 遗传算法种群空间设为通用形式，两种算法的其他参数设置如下：最大进化代数为100，种群规模为20，染色体的选择方法为轮盘赌法，编码方式为浮点法，自适应交叉和变异机制中K 1 = 0.8，K 2 = 0.6，K 3 = 0.005，K 4 = 0.002。

燃气轮机设计工况点的自适应建模结果见表2所示。在收敛精度方面，优化前的相对误差最大值约为?1.357%，AGA 算法的最大误差约为?0.879%，GM-AGA 的最大误差约为?0.604%；而目标函数c F 由0.144%下降到0.0224%，改进后的自适应遗传算法效果较理想。在收敛的速度方面，收敛时间由0.537下降到0.235 s ，耗时减小了约56%，改进后的算法收敛速度大大提高。由表2可知，相对于AGA 算法进行自适应建模而言，GM-AGA 算法进行的自适应建模效果更好。

赵建峰等

Table 2. Relative error between calculated value and measured value at design operating condition 表2. 设计工况点计算值与测量值相对误差

参数优化前误差/%

优化后误差/%

AGA算法GM?AGA算法

n 0 0.699 ?0.299

G1 ?0.449 ?0.426 0.314

T1 ?0.550 ?0.133 0.194

P1 ?1.010 ?0.0665 ?0.217

T2 ?0.942 ?0.879 ?0.101

T3 ?1.357 ?0.857 ?0.604

T4 0.604 0.541 ?0.142

P2 ?0.628 0.553 0.221

T5 ?1.036 ?0.208 0.369

e 0.307 0.345 ?0.191

F0.144 0.0460 0.0224 收敛时间------ 0.537 0.235

5. 结束语

本文结合实际工程应用提出了灰度预测机制，综合考虑单个部件因制造、使用及环境等因素导致的特性偏差，进而将遗传优化算法的种群空间进一步合理收缩．改进的自适应遗传算法能很快接近最优解。但是该算法也存在以下几个问题需要进一步进行研究：一是可能存在不合理的缩小种群空间导致只能获得局部最优；二是本算法未与其他改进的自适应遗传算法进行对比，其优势有待进一步进行评估；三是试验数据可能会对算法的精度产生较大的影响，需要进一步研究试验数据质量对算法影响程度。

基金项目

国防科技创新特区基金资助项目(17-H-05-ZT-002-041-01)。

参考文献

[1]刘前智. 多级轴流压气机非设计性能的数值预估[J]. 航空动力学报, 2004, 19(1): 108-112.

[2]段守付, 樊思齐, 卢燕. 航空发动机自适应建模技术研究[J]. 航空动力学报, 1999, 14(4): 440-443.

[3]谢光华, 曾庆福, 张燕东. 航空发动机仿真模型参数自适应校正[J]. 航空动力学报, 1998, 13(1): 37-40.

[4]陈玉春, 黄兴, 徐思远, 等. 涡轮发动机部件特性自适应模型的确定方法[J]. 推进技术, 2008, 29(2): 214-218.

[5]王永华, 李本威, 孙涛, 等. 基于部件特性自适应的涡扇发动机仿真[J]. 燃气涡轮试验与研究, 2009, 35(2):

11-13.

[6]Li, Y.G. and Pilidis, P. (2010) GA-Based Design-Point Performance Adaptation and its Comparison with ICM-Based

Approach. Applied Energy, 87, 340-348. https://https://www.wendangku.net/doc/de14576993.html,/10.1016/j.apenergy.2009.05.034

[7]刘永葆, 贺星, 黄树红. 基于改进遗传算法的燃气轮机自适应建模[J]. 航空动力学报, 2012, 27(3): 695-700.

[8]应雨龙, 李淑英. 一种基于粒子群优化算法的燃气轮机自适应热力计算方法[J]. 燃气轮机技术, 2015, 28(4):

48-54.

[9]Stamatis, A., Mathioudakis, K. and Papailiou, K.D. (1990) Adaptive Simulation of Gas Turbine Performance. Journal

of Gas Turbine and Power, 112, 168-175. https://https://www.wendangku.net/doc/de14576993.html,/10.1115/1.2906157

[10]凌玲, 胡于进, 王青青, 等. 基于改进遗传算法的孔群加工路径优化[J]. 华中科技大学学报(自然版), 2009, 37(8):

88-91.

[11]吕东洋. 等时序沉降观测GM(1, 1)预测精度分析[J]. 测绘与空间地理信息, 2018, 41(11): 194-197.

赵建峰等

[12]孔雪, 王丽, 冯益华. 灰色预测GM(1, 1)模型应用现状与展望[J]. 齐鲁工业大学学报, 2018, 32(6): 49-53.

[13]张乔斌, 王文华, 贺星. 基于RBF神经网络的燃气轮机监测参数非线性时间序列预测[J]. 海军工程大学学报,

2018, 30(2): 96-99.

1. 打开知网页面https://www.wendangku.net/doc/de14576993.html,/kns/brief/result.aspx?dbPrefix=WWJD

下拉列表框选择：[ISSN]，输入期刊ISSN：2325-677X，即可查询

2. 打开知网首页https://www.wendangku.net/doc/de14576993.html,/

左侧“国际文献总库”进入，输入文章标题，即可查询

投稿请点击：https://www.wendangku.net/doc/de14576993.html,/Submission.aspx

期刊邮箱：dsc@https://www.wendangku.net/doc/de14576993.html,