A Novel Fuzzy Model Predictive Control of a Gas Turbine in the Combined Cycle Unit

The complex characteristics of the gas turbine in a combined cycle unit have brought great difficulties in its control process. Meanwhile, the increasing emphasis on the efficiency, safety, and cleanliness of the power generation process also makes it significantly important to put forward advanced control strategies to satisfy the desired control demands of the gas turbine system. Therefore, aiming at higher control performance of the gas turbine in the gas-steam combined cycle process, a novel fuzzy model predictive control (FMPC) strategy based on the fuzzy selection mechanism and simultaneous heat transfer search (SHTS) algorithm is presented in this paper. The objective function of rolling optimization in this novel FMPC consists of two parts which represent the state optimization and output optimization. In the weight coefficient selection of those two parts, the fuzzy selection mechanism is introduced to overcome the uncertainties existing in the system. Furthermore, on account of the rapidity of the control process, the SHTS algorithm is used to solve the optimization problem rather than the traditional quadratic programming method. The validity of the proposed method is confirmed through simulation experiments of the gas turbine in a combined power plant. The simulation results demonstrate the remarkable superiorities of the adopted algorithm with higher control precision and stronger disturbance rejection ability as well as less optimization time.


Introduction
At present, thermal power generation still occupies the leading position in the power generation industry in China [1,2].Meanwhile, the increasing attention paid to energy conservation and emission reduction [3] has promoted the emergence and development of clean energy-based power generation.As a readily available clean energy, natural gas shows outstanding superiorities in environmental sustainability, energy efficiency, and security [4].Therefore, there are certain development prospects using natural gas as fuel in power generation.
The combined cycle unit which consists of the gas turbine, heat recovery boiler, and steam turbine [5] is one of the maturely developed power generation devices using natural gas.In the generation process, the generator is driven by the gas turbine; meanwhile, turbine exhaust enters the heat recovery boiler and the residual heat of the gas is made full use of through the steam circulation system [6].In this way, the characteristics of the gas turbine cycle with high average endothermic temperature and the steam turbine cycle with low average exothermic temperature are combined based on the principle of thermal cascade utilization.Considering the aforementioned process, the gas turbine is regarded as one of the most critical components in the combined cycle unit.Accordingly, it needs to be precisely controlled.However, because of the complex characteristics of the gas turbine with nonlinearity, variable coupling, and uncertainty, conventional control methods have difficulty in achieving the desired performance [7,8].Thus, it is remarkably significant to select one of the advanced control strategies for the control of the gas turbine.
With the flourishing development of control methods in recent years, some advanced control algorithms have emerged continuously and obtained great achievements, such as fuzzy control [9][10][11], predictive control [12][13][14][15], and neural network control [16][17][18][19].In all of those methods, predictive control is developed in the industrial process directly [20].The predictive control consists of these three main components: predictive model, rolling optimization, and feedback correction [21].Compared with the other aforementioned control strategies, predictive control has stronger ability to handle the time-varying processes with variable constraints and disturbances.Furthermore, this high-performance strategy also can combine the practical requirements of the control process with theoretical derivation maximally.Consequently, the predictive control will be the selected control method in this paper.Since its appearance, predictive control has attracted the great attention of experts and scholars of the world.Among all of the existing predictive control algorithms, the generalized predictive control (GPC) method [22,23] and the model predictive control (MPC) method [24,25] are commonly used.The most difference between these two predictive control strategies is that the GPC algorithm is based on the controlled autoregressive integrated moving average (CARIMA) model while the MPC method is based on the state space model [26].In view of the state space model, it can not only reflect the internal state of the system but also reveal the relationship between the internal state and external input/output variables [27].Therefore, the MPC algorithm will be used as the starting point of the novel algorithm adopted in this paper.In El-Ferik et al. [28], a distributed nonlinear model predictive control (NMPC) method was put forward for formation control of constrained autonomous vehicles in consideration of the communication bandwidth limitation and transmission delays.The collision was avoided effectively through the spatially filtered potential field, and the satisfactory performance of the presented scheme was also demonstrated through simulation experiments in either strongly connected or weakly connected networks.However, the model uncertainty and disturbances were not taken into consideration which would cause the reduction of control performance and universality of the proposed algorithm.In addition to that, the model predictive control (MPC) algorithm was applied to the coalconveying process of a coal-fired power plant [29].In that control scheme, the energy models of the belt conveyor and crusher were constructed and employed to formulate an open-loop energy efficiency optimization problem.Then, based on the open-loop optimization problem, a close-loop MPC strategy was adopted and applied to handle the nonlinearity of the controlled system.This control strategy showed strong robustness, adaptability, and energy-saving ability.Unfortunately, the performance reduction caused by the uncertainties of the system was inevitable.Considering the control of the complex nonlinear process with strong coupling and model mismatch, an extended nonminimal state space predictive control (ENMSSPC) strategy combined with a modified linear quadratic regulator (LQR) was put forward for the improvement of the control performance [30].Although this method showed better properties with ideal tracking performance and strong disturbance rejection ability, the accompanying increased calculation burden was ignored at the same time.Furthermore, the distributed stochastic MPC [31,32] was adopted for these control systems with uncertainty and probabilistic constraints as well as disturbances.According to the distribution manner, only one subsystem was selected at each time interval.Meanwhile, the probabilistic constraints were transformed into a set of deterministic constraints through the probabilistic distribution.The outstanding advantages of this strategy with assurances of recursive feasibility and close-loop stability for any updated sequence had been verified through numerical experiments.Nevertheless, the rapidity of the control process should also be taken into account in the design of the controller.
Considering that there are many values that cannot be expressed accurately in most control processes, the concept of "fuzzy" was introduced and applied in the construction of the control strategy to handle uncertainties existing in the controlled system [33,34].Aiming at a better dynamic response and reduction in the low-frequency ripple of integrated LED drivers, a fast fuzzy logic control algorithm based on a Takagi-Sugeno-Kang fuzzy controller was put forward [35].The high-speed mathematical engine of a digital signal controller was used to increase the calculation speed which contributed to a less settling time of the control process.Unfortunately, this method had a less positive effect on the improvement of control precision.The accurate mobile positioning in the cellular network is a matter of great challenge when it comes to the uncertainty in the signal intensity measurement process.Therefore, the fuzzy reasoning was added to a Kalman filter for improving the accuracy of localization [36], and the formed fuzzy extended filter showed remarkable feasibility and superiority compared with single Kalman filter method or fuzzy control strategy [37].What is more, the fuzzy logic also can be satisfactorily combined with the MPC algorithm.In Boulkaibet et al. [38], the T-S fuzzy model of each typical operating condition was identified by particle swarm optimization (PSO) clustering method and Kernel ridge regression algorithm initially.Then, the predictive control was introduced based on the fuzzy model, and the resulting adaptive fuzzy predictive control strategy could deal with the disturbances and parameter variations effectively.Nevertheless, the output variations and operating time of the control process were not taken into account.In view of the nonlinearities and various disturbances inevitably existing in a complex multivariable system, the fuzzy model predictive control (FMPC) method has become a desired candidate for the promotion of the control performance [39].The application of extended state observer-based FMPC algorithms on the boiler-turbine system in the ultrasupercritical unit of the power plant had handled these problems and obtained better tracking ability [40], whereas the rapidity of this method needed to be verified.Thus, considering the complex characteristics of the controlled system and the significant superiorities of the FMPC mentioned above, a new FMPC strategy will be constructed for the control of the gas turbine system in a combined cycle unit.
As one of the emerging intelligent algorithms in recent years, the bionic algorithm has developed rapidly with many branches, in which the swarm-based optimization technique [41,42] stands out and has been extensively applied in modeling and control processes.Model identification of an 2 Complexity unknown plant is important for the subsequent controller design to meet desired control specifications.Therefore, in the modeling process of an unknown digital system, a multiple subpopulation-based modified particle swarm optimization (PSO) algorithm was introduced to solve the coefficients of the infinite impulse response (IIR) digital filter [43].Although this method was able to overcome the multimodal error surface problem effectively, its precise identification would decrease when the order of the controlled system has increased.In addition, an artificial bee colony (ABC) algorithm and a grey wolf optimization (GWO) algorithm were used for solving portfolio optimization problems with cardinality constrains [44] and smart wind speed forecasting [45], respectively.The solution approach combined with the ABC algorithm had reduced the computational difficulties from large-scale optimization problems to a certain extent.Meanwhile, the proposed new hybrid strategy with the addition of GWO algorithm had ensured high accuracy and stability in the prediction process of wind speed.However, similar to most other swarm intelligence algorithms, it is difficult to achieve a win-win situation considering both the accuracy and the rapidity of the control process.Heat transfer search (HTS) algorithm is a novel method for solving the optimization problem which was presented in 2015 [46].On the basis of the law of thermodynamics in the heat transfer process, the molecule which consisted of the decision variables of the optimization problem made every effort to maintain the thermal equilibrium between itself and the surrounding environment through the three manners of heat transfer: conduction, convection, and radiation.This algorithm had won the first place in the Friedman rank test compared with some other optimization methods, such as PSO algorithm, ABC algorithm, and Yin-Yang pair optimization.Furthermore, in view of the expensive calculating burden in the iterative optimization process, a simultaneous heat transfer search (SHTS) was adopted and its lower computational complexity had been verified through simulation experiments [47].Consequently, a SHTS algorithm considering the variable constraints will be utilized in this paper to decrease the computational complexity and optimization time of the control process.On account of the content mentioned above, a novel fuzzy model predictive control (FMPC) which combined fuzzy selection mechanism with the SHTS algorithm is proposed for the promotion of control performance.Firstly, the objective function of the presented algorithm is made up of two parts which represent the state optimization and output optimization.These two parts occupy different weights according to the variations of state variables and output error of the system.Then, considering the uncertainties of the control system, the fuzzy selection mechanism is introduced for the adaptive selections of the above two weighting coefficients in each time interval.What is more, aiming at the reductions of computational burden and optimization time of the control process, the quadratic programming solving process in rolling optimization of conventional MPC is replaced by a SHTS algorithm while taking the variable constraints into account.This SHTS method is a less subjective strategy with two artificially defined parameters.Unlike the original HTS method, the whole population of SHTS is divided into three groups randomly, and each group is assigned to one of three modes of heat transfer.In this way, a great decline of the optimization time of the control process is obtained easily.Finally, the proposed algorithm will be utilized in the control of the gas turbine system in the combined cycle unit and compared with some other similar algorithms.The main contributions of this paper are summarized as follows briefly: (i) A novel FMPC algorithm is put forward based on the fuzzy selection of parameters while taking the state optimization problem and output optimization problem into consideration simultaneously.The addition of the fuzzy selection mechanism can overcome the uncertainties of the complex control system effectively.Furthermore, the precision of the control process also has great promotion (ii) Based on the FMPC algorithm described in (i), a HSTS method which takes variable constraints into account is adopted to replace the quadratic programming solving process in rolling optimization of the traditional MPC algorithm.This replacement can decrease the optimization time of the control system to a certain extent (iii) The presented algorithm is firstly applied in the control of the gas turbine in a combined cycle unit.The simulation results demonstrate that this control strategy can overcome the complex characteristics of a multivariable system effectively.Meanwhile, it is proved to be valuable in the promotions of control precision, disturbance rejection ability, and optimization speed.
The rest of paper is organized as follows.The gas turbine system in the combined cycle unit is described briefly in Section 2. The construction process of the proposed novel FMPC algorithm is introduced in Section 3. Next, the simulation results are shown to verify the superiorities of this algorithm in Section 4. Finally, some general conclusions are presented in Section 5.

Problem Description
2.1.The Combined Cycle Unit in a Power Plant.In order to increase the overall efficiency of electric power plants, multiple thermodynamic cycles can be combined to fully utilize energy.The term "combined cycle" refers to the combination of different thermodynamic processes.Compared with the simple cycle, the combined cycle can achieve the utmost utilization of energy and produce more electricity from the same fuel.
The combined cycle is always assembled by a topping cycle and a bottoming cycle, and its most common type is combined with the Brayton cycle operated by a gas turbine and the Rankine cycle operated by a steam turbine separately, which is usually called the combined cycle gas turbine (CCGT).The principle of the combined cycle is that the 3 Complexity temperature of the working fluid is still high enough to supply power for the subsequent engine after completing the cycle in the first engine.The operating process of the combined cycle can be described briefly as follows.The gas is burned in a gas turbine and becomes exhaust gases with temperature range from 550 °C to 610 °C.Meanwhile, the heat obtained in the burning process is converted to electric power through a coupled generator.Then a heat recovery steam generator is employed to produce steam and hot water by capturing exhaust heat.Next, a steam turbine is driven by the steam to deliver additional electricity by driving the corresponding generator.Finally, these two electric powers are transported to the grid synchronously.The diagram of the combined cycle unit described above is structured in Figure 1.
The basic combined cycle unit consists of two cycles as shown in Figure 2: the cycle 1-2-3-4-1 depicts the heat transfer process existing in the high temperature region while the cycle a-b-c-d-e-f-a represents the Rankine steam cycle existing in the lower temperature region.In addition, the heat recovery steam generator captures the exhaust energy to heat water into steam and deliver it to the steam turbine.Furthermore, the steam turbine transfers its energy to the generator, where the energy is converted into additional electricity.

Description of the Gas Turbine
System.The efficiency of the combined cycle gas turbine depends on the performance of the gas turbine which is regarded as the core equipment of the combined cycle power plant.However, the existing conventional control strategies cannot deal with the complex characteristics of the gas turbine effectively.As a result, it is necessary to develop advanced control schemes based on the dynamic characteristics of the gas turbine to meet grid dispatch requirements quickly.
There are three main components of the gas turbine system diagram which are depicted briefly in Figure 3: compressor, combustion chamber, and turbine.First, the air is compressed by the gas turbine and mixed with fuel.Then, the mixture is heated to a high temperature and high pressure gas to move through the turbine blades and make the turbine spin.Finally, a generator is driven by the mechanical energy from high-speed rotation of the turbine to produce electricity.
High generating efficiency and stable rotor speed, as well as the ideal exhaust temperature and exhaust flow are main control objectives of a gas turbine system.On the one hand, the turbine inlet temperature cannot be too high when it is limited by blade thermal stress.On the other hand, if this temperature decreases excessively, the efficiency of the combined cycle power plant is also reduced.Therefore, turbine inlet temperature should be kept as close as possible to the ideal value.What is more, to ensure the efficiency of the power plant, not only the fuel flow but also the inlet guide vane of the gas turbine should be adjusted to maintain a certain fuel-air ratio when the load of the combined cycle unit has changed.
Considering the most significant factors during the operating process, the gas turbine system can be simplified as a system with two input variables and four output variables as illustrated in Figure 4.
In Figure 4, G f and IGV are two input variables which represent the fuel flow and inlet guide vane, respectively; in addition, the output variables N op , n rs , T ex , and G ex stand for output power, rotor speed, exhaust temperature, and exhaust flow, respectively.Although this simplified multivariable model can reveal the main features of the gas turbine, obviously, its complex characteristics with nonlinearity, uncertainties, and coupling are also difficulties in the controller design of the gas turbine system.Therefore, it is remarkably significant to make further improvement of the original control strategy to satisfy the aforementioned control objectives of the gas turbine system in a combined cycle unit.

Simultaneous Heat Transfer Search-Based Fuzzy Model Predictive Control
The local linear models of the gas turbine system can be obtained through the modeling method proposed in Hou et al. [48] which overcomes the nonlinearity of the complex  4 Complexity system effectively.Therefore, on the basis of the identified model, a novel simultaneous heat transfer search-(SHTS-) based fuzzy model predictive control (FMPC) algorithm is adopted in this part.Firstly, a novel FMPC strategy is constructed to improve the control precision while considering the uncertainties in the control system.Then, the SHTS optimization method is introduced for reducing optimization time.
3.1.The Novel FMPC Algorithm.The proposed FMPC algorithm is an improvement of the original MPC combined with fuzzy selection mechanism; it mainly has the structure shown in Figure 5.
In the diagram shown above, r k is the set point of the k-th time interval, u k is the input variable vector of the controlled process while y k is the output variable vector, and x k is the predictive state variable vector.Furthermore, ŷ k + j stands for the predictive output of the j-th steps later after the k-th time interval, e k is the error vector between predictive outputs and actual outputs, c 1 and c 2 are the weight coefficients of output optimization and state optimization, respectively, obtained by fuzzy reasoning from fuzzy selection mechanism.
A multivariable system with m input variables, n state variables, and q output variables can be expressed as follows: where N is the sampling time, x n k T is the state vector and u k = u 1 k , u 2 k , … , u m k T and y k = y 1 k , y 2 k , … , y q k T are the input variable vector and output variable vector, respectively; meanwhile, matrices A, B, and C represent the state matrix, input matrix, and output matrix, respectively, of the system.
Considering that Δx k = x k − x k − 1 and Δu k = u k − u k − 1 , the increment vector of the output variables is exhibited as Constructing a new state variable vector as x k = Δx k y k T , the corresponding state space model shown as (3) can be obtained.5 Complexity Next, the X, ΔU, and Ŷ can be calculated through where Then, the objective function should be structured and minimized to obtain the control increment ΔU.
In view of the output optimization problem, the objective function has the following expression.
in which, W k is the reference trajectory of the system and Q y is the error weight coefficient while R y is the control weight coefficient of output optimization.
Furthermore, the state optimization-based objective function can be defined as where Q x is the state weight coefficient and R x is the control weight coefficient of state optimization.Substituting ( 6) into ( 8) and ( 9), we omit items that are not related to the control effect.The simplified objective functions are formulated as Next, on account of constraints of control variables, a new objective function J will be constituted based on the subfunctions J 1 k and J 2 k . Complexity where 12 c 1 and c 2 represent the weight coefficients of J 1 and J 2 , respectively.
Consequently, on the basis of objective function (11), the following control law of presented FMPC can be achieved by calculating the partial derivation of ΔU.
The first step of the control law will be used in the control system in each time interval.Considering the parameter uncertainties existing in the controlled system and inherent fuzzy properties of variables, the fuzzy selection mechanism described in the next part is introduced in the selection process of weight coefficients c 1 and c 2 .

Fuzzy Selection of the Weight Coefficients.
In the fuzzy selection mechanism of weight coefficients adopted in this part, the sum of absolute value of system output error e i k i = 1, … , q, k = 1, … , N and the sum of absolute value of state variable variation Δx i k i = 1, … , q, k = 1, … , N shown in (14) are taken as inputs of fuzzy rules while c 1 and c 2 are regarded as corresponding outputs. 14 Set the fuzzy domains of E k and ΔX k as [0, 7] and [0, 100], respectively.Meanwhile, the domains of fuzzy outputs c 1 and c 2 are both assigned as [0, 1].In addition, all of the membership functions of inputs and outputs are defined as Gaussian functions for its ideal properties and the simplicity of the calculation process, and the shown seven fuzzy degrees are used in the description of each variable.

S = NB, NM, NS, ZO, PS, PM, PB 15
Then, membership function curves of input variables in this fuzzy selection system can be drawn as Figure 6.The Gaussian function curves of the output variables are equally available.
In the weight coefficient selection process, the value of c 1 will be enlarged when the sum of output errors of the controlled system is increased; at the same time, the value of c 2 will be promoted when the sum of state variable variations is increased.Combining this basic regulation principle with expert experiences, the resulting fuzzy rules are illustrated as In order to ensure the integrity of the fuzzy rules, all of the rules are listed in following Table 1.
Consequently, the weight factors of output optimization and state optimization can be chosen reasonably when given the inputs of fuzzy rules.However, this parameter selection process will increase the settle time of the whole control system.Under this circumstance, a SHTS algorithm will be presented for the replacement of the conventional quadratic programming solving method to reduce the optimization time.

Simultaneous Heat Transfer Search Algorithm in
Optimization Process.In order to solve a certain optimization problem through the quadratic programming, an objective function with standard quadratic programming form needs to be structured.Nevertheless, this procedure has increased the solving complexity to a great extent.On the contrary, the standardization step is not necessary in the optimization process of the simultaneous heat transfer search (SHTS) 7 Complexity algorithm mentioned.Therefore, the SHTS algorithm will be applied in every iterative optimization process to realize the promotions of optimization speed and control performance.The main steps of the SHTS method are summarized as follows.
Step 1. Initialize the population size NP, search dimension D and maximum iteration number K, as well as the upper bound ΔU max and lower bound ΔU min of the variable vector to be optimized, in other words, the control variable vector ΔU k k = 1, … , N mentioned above.
Step 2. Generate the initial population ΔU k with size of N P × D through where r is a random NP × D matrix with values between 0 and 1.
Step 3. Calculate the value of fitness function J k , in other words, the aforementioned objective function, based on the initial population.Then, select the member with minimum fitness as the best one and denote as ΔU best k , J best k k = 1, … , N .
Step 4. Divide the population into three equal parts randomly and distribute each part to one of the three heat transfer modes.Furthermore, denote the subpopulations with conduction phase, convection phase, and radiation phase as ΔU D k , ΔU V k , and ΔU R k , respectively.
Step 5. Implement following three heat transfer modes of each subpopulation.
Conduction phase: obtain new population through the equation described as where r D is a random NP/3 × D matrix with values between 0 and 1.
Convection phase: define c k = r V ΔU best k R, in which r V is a random matrix of 1 × D dimension with values between 0 and 1 and R is a matrix of 1 × D dimension with values between 1 and 2.Then, generate a new matrix Cc k which consists of NP/3 copies of c k ; therefore, the size of matrix Cc k is NP/3 × D.Then, obtain the new population through.
Radiation phase: randomly divide the subpopulation ΔU R k into two equal parts ΔU R1 k and ΔU R2 k .Next, assign every member of these two parts to a partner which is selected from ΔU R k and denote the partner populations as ΔU R1 ′ k and ΔU R2 ′ k , respectively.Then, generate the new population determined as where r R1 and r R2 are both the random matrices of the NP/6 × D dimension with values between 0 and 1.
Step 7. Modify the new population through (21) considering its feasibility in the search region.
Step 8. Calculate the fitness J * k of the new population and make a greedy choice with the initial optimal fitness J best k .If the new fitness J * k is smaller than the original J best k , the corresponding new population ΔU * k will be selected as the updated optimum solution ΔU best k and J best k = min J * k , J best k , else, the new population will be abandoned and search for another population will begin.
Step 9.If the control precision has satisfied the ideal demands or the number of iterations has reached the limitation, end this optimization process and apply the first step of ΔU best k to the controller.Then, start a new optimization process.Else, go back to step 2.
The above discussed optimization process combined with SHTS can be simplified as Figure 7.

Illustrative Example
The remarkable advantages of the proposed control strategy should be demonstrated through illustrative examples to prove its feasibility and practicability.Therefore, the model of the gas turbine system of a combined cycle unit in Taiyanggong power plant which was identified through the T-S fuzzy modeling [48] is introduced as ( 22) to verify the performance of the adopted algorithm.The verification process is made up of the comparison experiments of tracking performance and disturbance rejection ability with conventional MPC (MPC), fuzzy MPC (FMPC), FMPC with PSO (FMPC + PSO), and FMPC with GWO (FMPC + GWO).

Parameter Selections.
Parameter selection is a significantly important part in the controller design.The empirical method is adopted in this paper because of its easy operation.In addition, some parameters in PSO, GWO, and SHTS can be obtained or fixed automatically through the stochastic technique.Thus, there are only two parameters that need to be determined, namely, the population size NP and maximum iteration number K. Considering the fairness of comparison experiment, the corresponding parameters in each algorithm are selected as equal values in every simulation process and listed in Table 2.
What is more, the error weight coefficient vector and control weight coefficient vector of output optimization are set as Q y = diag 1,10,1,10 and R y = diag 0 1,0 1 , respectively; meanwhile, the state weight coefficient vector and control weight coefficient vector of state optimization are set as

Simulation
Results.This part is made up of the tracking performance test and disturbance rejection verification of the proposed controller.At first, based on the selected parameters, the following simulation experiment is executed to test the tracking performance of the controller.The initial stable operating point of the gas turbine system, in other words, the output variable vector of the controlled system, is set as 180 79, 3000, 642 27, 1106 8 T .Then, the set point of the system is changed to 164 98, 3000, 652 7, 1040 8 T at 20 s; the responses of the   In Figures 8 and 9, the black dotted lines present the reference trajectories of the system and the blue solid lines stand for responses of variables controlled by MPC while the magenta dotted lines are the responses of variables controlled by the FMPC.Besides, the cyan solid lines, green solid lines, and red solid lines present the responses of variables controlled by the FMPC which are combined with PSO, GWO, and SHTS, respectively.It is clear from these two figures that when the set point of the system has changed, the variables of the system controlled by the presented algorithm show high tracking performance with lower overshoot and shorter settling time compared with other   11 Complexity methods.On account of the state stability of the proposed algorithm, the state variable responses of the system controlled by MPC-, FMPC-, and SHTS-based FMPC are selected for comparison and illustrated in Figure 10.x 1 , x 2 , … , x 8 stand for eight state variables of the system.Although not all of the used control methods are included in Figure 10, it is enough for the state stability verification.Because the number of state variables in this gas turbine system is set as 8; therefore, Figure 10 consists of four subfigures and each subfigure obtains two state variables of the system controlled by the aforementioned three strategies.Moreover, the dot-dash lines, dotted lines, and solid lines stand for the state variables of the system controlled by MPC, FMPC, and FMPC which is combined the SHTS (FMPC + SHTS), respectively.From Figure 10, it is obvious that the responses of state variables of the system controlled by the presented algorithm illustrate smoother trends and reach a new steady state in a shorter period of time.Thus, the high state stability of the proposed control method is demonstrated through this comparative experiment.
Apart from tracking performance, the disturbance rejection ability is also an indispensable index to judge the performance of a controller.As a result, the following disturbance rejection experiments are finished for the further performance verification of the adopted algorithm.
(i) Apply a negative step signal with amplitude of 5 to the power output of the system while keeping the other three output variables unchanged.Then, the responses of the output variables and input variables are shown in Figure 11.
(ii) Apply a negative step signal with amplitude of 1 to the rotor speed of the system while keeping the other three output variables unchanged.The responses of the output variables and input variables are shown in Figure 12.
(iii) Apply a negative step signal with amplitude of 10 to the temperature of the exhaust gas while keeping the other three output variables unchanged.Then, the responses of the output variables and input variables are depicted in Figure 13.
(iv) Apply a negative step signal with amplitude of 50 to the exhaust gas flow of the system while keeping the other three output variables unchanged.The responses of the output variables and input variables are shown in Figure 14.
It is very clear from  that the output variables of the system controlled by the presented control method have the smallest overshoot, and the three other 12 Complexity output variables except for the output variable with disturbance can return to the original stable states over the shortest period of time.What is more, the responses of input variables also illustrate the stronger disturbance rejection ability of the controller.It is worth mentioning that the FMPC which is combined with GWO also has the highest settling speed among these five algorithms; however, there are high stable errors existing in the input and output variables which would not be ignored in the control process.Therefore, the presented algorithm is proved to be the best one with remarkable high control performance compared with the other four control methods.
In order to verify the rapidity of the rolling optimization process of the proposed algorithm, the optimization times of above five control strategies in tracking performance tests and disturbance rejection experiments are recorded and listed in following table for comparison.
As is shown in Table 3, the optimization times of MPC are smallest among all of the used controllers because of the absence of the fuzzy selection mechanism.Although the fuzzy selections of the weight coefficients in the objective function have decreased the optimization speed, the higher control performance of FMPC with lower overshoot and shorter settling time as well as higher accuracy has exhibited 13 Complexity its remarkable superiority compared with traditional MPC.Furthermore, aiming at promoting the optimization speed of FMPC while maintaining its priorities, the quadratic programming solving process of the initial FMPC is replaced by a swarm-based optimization method such as PSO, GWO, and SHTS listed in this table.It is obvious that both the FMPC + GWO and the FMPC + SHTS have shorter periods of optimization time compared with FMPC and FMPC + PSO, whereas the control process of FMPC + GWO is a poor regulation process with high steady state errors.Thus, the proposed improved FMPC which is combined with SHTS has demonstrated the shortest optimization time, the highest tracking performance, and stronger disturbance rejection ability simultaneously.

Conclusions
A novel FMPC algorithm which combined fuzzy selection mechanism and the SHTS while considering the variable constraints is put forward in this paper.Initially, a new objective function is constructed based on the comprehensive consideration of the state optimization problem and output optimization problem.Then, in view of the weight coefficient selections of these two parts in the objective function, a fuzzy selection mechanism is established to deal with the uncertainties existing in the control process and improve the control precision of the controller.Next, in order to promote the optimization speed for the compensation of time increase caused by the fuzzy selection process, the SHTS is adopted in the rolling optimization of FMPC as a replacement of the initial quadratic programming method.Furthermore, the proposed novel FMPC strategy is applied to the gas turbine system of a combined cycle unit in a power plant.The
If E = NB and ΔX = NB, then c 1 = PM and c 2 = NM, ⋮ If E = PB and ΔX = PB, then c 1 = PB and c 2 = ZO 16

Figure 6 :
Figure 6: Gaussian membership function curves of fuzzy input variables.
Initialize the population size NP, search dimension D, maximum iteration number K as well as the upper and lower bound of optimal variables Calculate the fitness J(k) of population ΔU(k) and select the best member (ΔU best (k), J best (k))Divide the population into three equal parts randomly and assign one of the three heat transfer modes to each partConduction phaseConvection phaseRadiation phaseGenerate the new population ΔU ⁎ (k)Modify the new population ΔU ⁎ (k) into search region Calculate the fitness J ⁎ (k) of the new population and make greedy choice with initial minimum fitness J best (k) population ΔU(k)

Figure 7 :
Figure 7: Flow chart of the SHTS algorithm.

Figure 8 :
Figure 8: Responses of output variables when the set point has changed.

Figure 9 :
Figure 9: Responses of input variables when the set point has changed.

Figure 10 :
Figure 10: Responses of state variables when the set point has changed.

Figure 11 :
Figure 11: Responses of (a) output variables and (b) input variables when a negative step signal with amplitude of 5 is applied to the power output of the system.

Figure 12 :
Figure 12: Responses of (a) output variables and (b) input variables when a negative step signal with amplitude of 1 is applied to the rotor speed of the system.

Figure 13 :
Figure 13: Responses of (a) output variables and (b) input variables when a negative step signal with amplitude of 10 is applied to the temperature of the exhaust gas.

Figure 14 :
Figure 14: Responses of (a) output variables and (b) input variables when a negative step signal with amplitude of 50 is applied to the exhaust gas flow of the system.

Table 1 :
Fuzzy rules of weight coefficients selection process.

Table 2 :
Parameter selection results of the control process.
system output variables, input variables, and corresponding state variables are shown as respectively.

Table 3 :
Optimization times of tracking performance tests and disturbance rejection experiments.