Multiobjective Optimization of a Fractional-Order PID Controller for Pumped Turbine Governing System Using an Improved NSGA-III Algorithm under Multiworking Conditions

In order to make the pump turbine governing system (PTGS) adaptable to the change of working conditions and suppress the frequency oscillation caused by the “S” characteristic area running at middle or low working water heads, the traditional singleobjective optimization for fractional-order PID (FOPID) controller under single working conditions is extended to amultiobjective framework in this study. To establish themultiobjective FOPID controller optimization (MO-FOPID) problemundermultiworking conditions, the integral of the time multiplied absolute error (ITAE) index of PTGS running at low and high working water heads is adopted as objective functions. An improved nondominated sorting genetic algorithm III based on Latin hypercube sampling and chaos theory (LCNSGA-III) is proposed to solve the optimization problem. The Latin hypercube sampling is adopted to generate well-distributed initial population and take full of the feasible domain while the chaos theory is introduced to enhance the global search and local exploration ability of the NSGA-III algorithm. The experimental results on eight test functions and a real-world PTGS have shown that the proposed multiobjective framework can improve the Pumped storage units’ adaptability to changeable working conditions and the proposed LCNSGA-III algorithm is able to solve the MO-FOPID problem effectively.


Introduction
With the continuous expansion of the power system scale, the power grid's requirements for power quality, safety, and intelligence have been constantly improving.Pumped storage units (PSUs) have played an important role in maintaining the balance of power supply and demand because of their fast start-up and shutdown speed, flexible working condition conversion, excellent peak-load regulation, and frequency regulation ability [1][2][3].Pump turbine governing system is the core control system of the pumped storage power station which is responsible for stabilizing the unit frequency and regulating the unit power [4,5].Due to the huge flow inertia of the long-distance water pipeline and the existence of the unstable "S" characteristic area, the optimal control of PTGS is highly complex [6].Therefore, it is of great theoretical value and practical significance to explore optimization methods for PTGS and research new control laws.The control quality of PSU and the dynamic response performance of PTGS can then be improved.
The classical Proportional-Integral-Derivative (PID) controller is widely used in the optimal control of PTGS because of its simple and reliable structure and easy adjustment of control parameters [7][8][9][10].However, due to the strong nonlinear characteristics of different parts of PSU and the changeable working conditions, the traditional PID controller often 2 Complexity fails to realize the global optimization of PTGS.The control problems of phase modulation instability and no-load frequency fluctuation are becoming increasingly prominent.As an extension of the classical PID controller, the fractionalorder PID controller (FOPID) has attracted the attention of many scholars for its better adaptability and flexibility and greater potential to obtain better control performance [11][12][13][14][15]. Li et al. [16] proposed an improved gravitational search algorithm using the Cauchy and Gaussian mutations to adjust the parameters of the FOPID controller automatically.A number of tests have shown that the FOPID controller can improve the dynamic characteristics and stability of regulation frequency of PSUs.In Xu et al. [17], a robust nonfragile FOPID controller was proposed for PTGS.The parameters of the FOPID controller were selected using the bacterial foraging algorithm and multiscenario analysis functions.The FOPID controller turned out to have obtained higher robustness and stability compared with the traditional PID controller.And it can fully track the nonlinear characteristics of PTGS in the "S" area.Xu et al. [18] proposed an adaptively fast fuzzy fractional-order PID (AFFFOPID) control method for PSUs by combining a fuzzy fractional-order PD controller with a fuzzy fractional-order PI controller.Experiments of PTGS at various water heads under unload running condition have shown that the controller can effectively improve the performance and control quality of PSUs during the transient process.
Apart from the advantages of the FOPID compared with the PID controller, one of the most important and challenge issues of the employment of FOPID controller in PTGS is the optimal optimization of its parameters.With the continuous development of optimization algorithms and control theory in recent years, scholars have combined the FOPID controller with intelligent algorithms to achieve the optimal tuning of control parameters and improvement of control laws for PTGS [19].The related works on optimal optimization for PID controller or hydroturbine governing system (HTGS) can also provide meaningful reference for the research of FOPID controller for PTGS.Fang et al. [20] developed an improved particle swarm optimization (PSO) algorithm for optimal tuning of PID control parameters for water turbine governor.Simulation results have demonstrated the stable convergence characteristic and good computational ability of the developed optimization stagey.Kou et al. [21] proposed a novel BFO-PSO algorithm by introducing PSO into the framework of the bacterial foraging optimization (BFO) algorithm to improve the control performance of the PID controller for HTGS.The advantages of the proposed algorithm to the BFO and PSO algorithms have been demonstrated through experiments at real working conditions.Wang et al. [22] proposed a three-stage start-up strategy for PSUs by opening guide vanes to a large opening degree firstly and then reducing the opening degree and finally switching to the PID controller.The switching time and PID control parameters are optimized synchronously using an integrated optimization scheme based on artificial sheep algorithm (ASA).Simulation results under various water heads have shown that the control strategy can shorten the start-up time and reduce the speed oscillation.Some researchers have paid attention to multiobjective designing to consider multiobjectives that reflect the specific characteristics of the control system [23,24].Zamani et al. [25] developed a multiobjective cuckoo search approach to optimize the parameters of a FOPID controller.Sánchez et al. [26] proposed a multiobjective optimization strategy for identifying the optimal solution for a robust FOPID controller.Zhao et al. [27] proposed a parameter tuning scheme based on two lbests multiobjective PSO (2LB-MOPSO) to minimize the integral squared error and balanced robust performance criteria of a robust PID controller simultaneously.Chen et al. [28,29] put forward a parameter tuning scheme to optimize the integral of the squared error (ISE) and the integral of the time multiplied squared error (ITSE) performance indices of PTGS simultaneously.The adaptive grid PSO (AGPSO) and the chaotic nondominated sorting genetic algorithm II (NSGA II) were adopted to realize the optimal control of the PID and the FOPID controllers, respectively.The proposed multiobjective PID and FOPID controller turned out to have achieved better control effects than the compared methods.
It is noticed that most of the research works on optimal control of the PID or FOPID controllers are designed under a single working condition.The traditional multiobjective optimization control usually adopts a set of objective functions to obtain the Pareto optimal solutions for a certain condition and then select the compromise optimal solution.The traditional multiobjective optimization control is essentially an optimization scheme for single working conditions.However, the characteristics of the controlled object of PTGS are not only related to the nonlinear characteristics of the pump turbine, the penstock system, the generator, and the other parts of PTGS but also vary with the changes of the working conditions.The optimal control of a single working condition is often at the expense of the deterioration of some other working conditions.The optimal control of PTGS under multiworking conditions should be a process of overall tradeoff and cannot be limited to the pursuit of the optimal control under a certain working condition.In order to enhance the adaptability of PTGS to the change of working conditions, a multiobjective optimization framework which takes into account the integral of the time multiplied absolute error (ITAE) index of PTGS running at multiworking conditions is constructed.The nondominated sorting genetic algorithm III based on the Latin hypercube sampling and chaos theory (LCNSGA-III) is proposed to optimize the control parameters of the FOPID controller for PTGS under multiworking conditions.
The rest of this paper is arranged as follows: Section 2 gives a brief introduction to the multiobjective optimization problem; Section 3 builds the mathematical model of MO-FOPID under multiworking conditions; Section 4 proposes the LCNSGA-III algorithm based on Latin hypercube sampling and the chaos theory to solve the MO-FOPID problem; Section 5 employs eight benchmark functions and a real-world PTGS to verify the effectiveness of the LCNSGA-III algorithm and the developed multiobjective optimization framework; Section 6 gives the conclusions.

Multiobjective Optimization Problem (MOOP)
The purpose of the single-objective (SO) optimization problem is to obtain the optimal solution by searching for the minimum or maximum value of one single objective function.However, optimization problems in scientific research or engineering application usually contain not only one objective.These objectives are sometimes in concordance with each other, but sometimes there are conflicts between them.The objective of the multiobjective optimization problem (MOOP) is to search for the optimal solutions of all objectives, which is also known as the Pareto optimal solution [30].A typical MOOP with D decision variables, N objective functions, and m+k constraints can be described as follows: min  () = (min where Pareto Optimal Solution.Pareto optimal solution is a solution that cannot be dominated by any solution in the feasible region  which means that if and only if ¬∃ ∈  :  ≺  * ,  * is the Pareto optimal solution.
Pareto Optimal Set.For a given MOOP, the Pareto optimal set  * can be defined as Pareto Optimal Front.For a given MOOP, the Pareto optimal front  * can be described as

Problem Formulation of the MO-FOPID Problem under Multiworking Conditions
To formulate the multiobjective FOPID controller optimization problem (MO-FOPID) under multiworking conditions, the fractional calculus theory and fractional-order PID (FOPID) controller are first introduced.

Introduction to Fractional Calculus Theory and FOPID
Controller.Fractional calculus theory is an extension of the traditional calculus theory to the fractional systems.The fractional calculus theory can describe complex systems in an easy way with a clear physical meaning.Among the various fractional calculus operators, the Laplacian transformation defined by Caputo is the most usual mathematical expression for calculating the fractional-order time derivatives and has been widely used in fractional-order PID controller [28,31].Given a continuous derivable function (), Caputo's fractional derivative of order  can be defined as where 0    denotes the fractional calculus operator and Γ(⋅) denotes the Euler Gamma function.
Fractional calculus equations usually need to be transformed into algebraic equations.The Laplacian transformation of (5) under zero initial condition can be expressed as In recent years, with the further research and development of fractional calculus theory, the combination of fractional calculus theory and modern control theory is becoming more and more popular.Controllers based on fractional calculus have been implemented and applied in many research fields [12,25].The Oustaloup recursive filter and its improved version have been widely adopted to realize the discretization approximation of the fractional calculus operator.The expression of the Oustaloup filter is as follows: where   =   ( ℎ /  ) (++(1+)/2)/(2+1) ;    =   ( ℎ /   ) (++(1−)/2)/(2+1) ;  =   ℎ ;  denotes the order of the fractional calculus; (2N+1) is the order of filter; (  ,  ℎ ) denotes the expected fitting range; ,   , and  ℎ are determined according to the accuracy requirement of the numerical approximation.
The FOPID controller proposed by Professor Podlubny [32] is an extension of the classical PID controller.Compared with the classical PID controller, the range of control rate of the FOPID controller is much wider.The transfer function of the FOPID controller is as follows:  where  denotes the control deviation;  denotes the controller output;   ,   , and   represent the gain parameter, the integral order, and the differential order, respectively.The traditional PID controller is a special case of the FOPID controller when  = 1 and  = 1.Because the integral and differential orders are variational, the FOPID controller has better adaptability and flexibility and bigger potential to obtain better control performance [33].The structure of the FOPID controller for PTGS is shown in Figure 1.In Figure 1,   denotes the given unit speed;  denotes the unit speed;   denotes the permanent slip coefficient;   denotes the differential time constant;   ,   , and   represent the proportional, integral, and differential gain coefficients of the controller, respectively;  and  denote the integral and differential orders of the controller, respectively.

Description of the PTGS
System.PTGS is a complex nonlinear time-varying system of hydraulic, mechanical, and electrical connections.PTGS mainly contains five parts, namely, a controller, an electrohydraulic servomechanism system, a pump turbine, a generator, and a penstock system, where the controller and the servomechanism system consist of the speed governor of PTGS [4].A PID controller has always been employed as part of the speed governor of PTGS.In this study, the FOPID controller is designed for PTGS.The FOPID controller has been described in Section 3.1.In what follows, transfer functions of the other four connectors are illuminated.
(1) Servomechanism System.Servomechanism is the actuator of the governor of PSU.It is made up of an auxiliary servomotor and a main servomotor of which the transfer functions are as follows: where  1 and  V are the time constant and scale factor of the auxiliary servomotor, respectively, and   represents the main servomotor open-loop time constant.
(2) Pump Turbine.To reflect the complex nonlinear characteristics among water, machine, and electricity during the operation process accurately, the pump turbine model based on characteristic curves is constructed.The model of pump turbine is as follows: where  11 ,  11 , and  11 represent the unit torque, unit flow, and unit speed, respectively;  is the guide vane opening;   and   denote the functions of the moment and flow characteristic curves, respectively.Because of the strong nonlinear characteristics of the flow and moment characteristic curves, the improved Suter transform is introduced to transfer the flow and moment characteristic curves into WH and WM characteristic curves, respectively [4].The WH and WM characteristic curves of a pump turbine in a pumped storage power station in China using the improved Suter transform have been illustrated in Zhou et al. [4].In this study, the WH and WM characteristic curves in two dimensions in [4] have been changed to three dimensions.The three-dimensional WH and WM characteristic curves are illustrated as in Figure 2.
(3) Generator.The common first-order model [6,16,18] is adopted in this study to balance the pump turbine torque and the generator torque.The transfer function of the first-order model is as follows: where   and   are the inertia time constant and selfadjusting factor of the generator, respectively.
(4) Penstock System.Because of the fluid and tube wall elastic effects on the penstock, the second-order elastic water hammer model is exploited in this study by applying the second-order Taylor expansion on the nonlinear hyperbolic tangent function.The transfer function of the second-order elastic water hammer model is as follows:  where   denotes water flow inertia time constant,  represents the water head loss coefficient, and   is the reflection time of water hammer wave.

Multiobjective Optimization of FOPID Controller for PTGS.
For single-optimization of the FOPID controller under a certain working condition, the ITAE [34,35] is usually employed as the objective function to obtain satisfactory transient dynamic performance of the system.The ITAE index is defined as follows: where () denotes the relative deviation of the rotational speed of PTGS.The ITAE index considers the stable time and overshoot of the dynamic response of PTGS simultaneously.It can evaluate the speed and stability of the system at the same time.The smaller the ITAE, the better the speed and stability of the system.When operating at the working condition of low water head, the PSU is easy to fall into the "S" characteristic region, resulting in oscillation of the unit speed near the rated frequency.The working condition of medium or high water head, in the other way, is the most common working condition in the operation process of PSU.The optimal control of PTGS under different working conditions should be considered and researched to make the PTGS system better adapt to the changeable working environment.However, the ITAE for PTGS at working condition of low head (referred as ITAE 1 ) and that of high head (referred as ITAE 2 ) usually influence and restrict each other.In this study, the single-objective optimization of the FOPID controller for PTGS is expanded to multiobjective theoretical framework to find a compromise solution and ensure that the PTGS system can achieve relatively better control performance under changeable working conditions.To formulate the MO-FOPID problem, the ITAE 1 and ITAE 2 are adopted as objective functions.The five parameters of the FOPID controller for PTGS including the proportional coefficient   , the integral coefficient   , the differential coefficient   , the integral order , and the differential order  are taken as decision variables.The MO-FOPID problem can then be formulated as subject to where  1 (⋅) and  2 (⋅) are functions of   ,   ,   , ,  at working condition of low water head and high water head, respectively. min and  max are the lower and upper bounds of   ,   ,   , ,  of the FOPID controller, respectively.

Nondominated Sorting Genetic Algorithm-III Based on Latin-Hypercube Sampling and Chaos Theory (LCNSGA-III)
An improved version of NSGA-III based on Latin-hypercube sampling and chaos theory is developed to solve the proposed MO-FOPID problem under multiworking conditions.

Brief Introduction to NSGA-III.
The NSGA-III algorithm [36], first introduced by Deb and Jain in 2014, is a novel reference-point-based nondominated sorting genetic algorithm following the NSGA-II framework.Unlike the crowding distance operator exploited in NSGA-II [37], NSGA-III employed a reference point-based mechanism to make the Pareto optimal front well-distributed.The step-to-step Complexity procedures of the NSGA-III algorithm can be expressed as follows.
Step 1. Calculate the number of reference points () to be placed on the hyper-plane.
where  represents the number of objective functions and  denotes the number of divisions (For  = 3 and  = 4,  is calculated as 15).
Step 2. Generate NP individuals in the feasible region randomly to form the initial population of the NSGA-III algorithm, and record it as   .Set the iteration number as  = 1.
Step 3. Generate the offspring population   using the simulation binary crossover (SBX) operator and polynomial mutation operator.
Step 4. Let   =   ∪  and calculate the fitness value of each individual in   .
Step 6. Normalize the objectives and associate the solutions in   with the reference points.Delete the useless reference points and preserve solutions with higher rankings according to the niche preservation strategy to construct the next generation of population  +1 .
Step 7. If  <  max , skip to Step 3; else stop iteration and output the Pareto optimal set.
The crossover and mutation operators in Step 3 and the nondominated sorting operator in Step 5 of the NSGA-III algorithm are the same as those of the NSGA-II algorithm.In Step 6, the NSGA-III algorithm generates the next population using the reference point-based selection mechanism other than the crowding distance operator of NSGA-II.Readers can refer to [36] for more details about NSGA-III.The improvements of the NSGA-III algorithms including the Latin hypercube sampling-based initialization technique and the chaotic crossover and mutation operators are introduced in the following subsections.

Latin Hypercube Sampling Based Initialization Technique.
NSGA-III generates the initial values of the decision variables in the feasible region randomly to form the initial population of the algorithm.However, a large number of individuals may assemble into a local area of the feasible region because of the random initiation, leading to premature convergence in the iterative process.In order to make the individuals of the initial population well-distributed in the feasible region, the Latin hypercube sampling based initialization technique [38] is introduced to NSGA-III to improve its search performance.
The brief steps to generate the initial population of NSGA-III using Latin hypercube sampling are as follows.
Step 1. Suppose the size of the population is  and each individual in the population contains  discrete elements.The range of the discrete element   ,  = 1, 2, . . .,  of individual  can be divided into  equal mini zones: where (   <  <  +1  ) = 1/ and the value space of  can be divided into   small hypercubes finally.
Step 2. Generate a matrix  of which the dimension is  × , and every column of  is the full permutation of {1, 2, . . ., }.
Step 3. Generate an individual in each row of  randomly; an initial population with  individuals is then generated.

Chaotic Crossover and Mutation
where  represents the crossover coefficient. can be calculated according to the following: where  = rand (⋅) is a random number uniformly generated in [0, 1].
The NSGA-III algorithm generates a random number   to determine whether to apply the mutation operator or not.The polynomial mutation operator is implemented when   <   , where   represents the mutation distribution index.For the feasible solution   , the mutation individual is generated using the polynomial mutation operator: where  *  represents the mutation individual;    and    represent the upper and lower bounds of   , respectively;   is the mutation coefficient;   can be calculated as follows: where   = rand(⋅) is a random number uniformly generated in [0, 1].
The standard NSGA-III algorithm has excellent computational efficiency and stability.However, the NSGA-III algorithm sometimes may fall into the local optimal solution because of the insufficient exploration of the feasible region.Because of the ergodicity and stochasticity of chaotic sequences [39], the chaotic map which can generate chaotic sequences is introduced into the crossover and mutation operators of NSGA-III to enhance its global search and local exploration ability [40].In this study, the tent map [41] is employed to improve the crossover and mutation operators of NSGA-III.
The tent chaotic map with uniform distribution probability density can be expressed as where  () ∈ (0, 1) denotes the chaotic variable generated in the th iteration.
According to ( 19) and ( 21), the crossover and mutation operators of the standard NSGA-III algorithm need to generate two random numbers  and   in [0, 1], respectively.The two numbers  and   are generated using the tent chaotic map in the LCNSGA-III algorithm.The two random numbers  and   are generated according to the following: Since the frequency of the rotational speed is 50 Hz and the frequency disturbance is 2 Hz, the frequency perturbation is set as 4% of the rated frequency.and chaotic map, a total number of eight test functions including the ZDT1-4, 6 [42] and the DTLZ1-2, 5 [43] (shown in Table 1) are employed to test its performance.Among the eight test functions for MOOPs, the ZDT benchmark functions are two-objective MOOPs.The dimension of the decision variables and the number of the test functions of the DTLZ benchmark functions can be adjusted.In this study, the number of the objective functions of the DTLZ test functions is selected as three to display the Pareto front of threeobjective MOOPs.Four other typical multiobjective algorithms including NSGA-II, NSGA-III, MOEA/D, and PESA-II [30,44] are adopted as control group.The performances of the multiobjective evolutionary algorithms (MOEAs) are evaluated and compared using the widely used evaluation metrics of generational distance (GD) and Spread.GD is adopted to measure the mean value of the distance between the Pareto solution set and the real Pareto front.Spread is employed to describe the distribution uniformity of the Pareto optimal set [40].GD can be expressed as follows: where   denotes the Euclidean distance between the ith nondominated solution and the nearest nondominated solution on the real Pareto front and  denotes the size of the Pareto optimal set.The smaller the GD, the closer the Pareto solution set to the real Pareto front and the better the convergence accuracy.

Numerical Experiments and Analysis
Spread can be expressed as follows: where   denotes the distance between ith solution with its neighboring solution;  is the mean of   ;    denotes the Euclidean distance between the extreme solution of the real Pareto front and the boundary solution of the obtained Pareto optimal set; M denotes the number of objective functions.The smaller the Spread, the better the distribution of the solutions.

Results
Analysis.The parameters of the five algorithms for MOOPs with two objectives are set as follows: the population size is set as 100 while the number of iterations is set as 300.The parameters of the five algorithms for MOOPs with three objectives are set as follows: the population size is set as 150 while the number of iterations is 500.The crossover and mutation probability of the five algorithms are set as 0.7 and 0.3, respectively.The neighborhood size of MOEA/D is set as 20 [45,46].All the experiments are implemented in Matlab environment.Because of the random initialization of the MOOPs, all the experiments have been repeated ten times independently to eliminate the effectiveness of randomness.The average GD and Spread of the five algorithms for the eight benchmark functions are given in Tables 2 and 3, respectively.Numbers in bold represent the optimal GD and Spread among the five algorithms for the eight benchmark functions.
As can be seen from Table 2, the GD of the Pareto optimal solution obtained by the newly developed LCNSGA-III is the best out of the total eight benchmark functions.Although the GD of LCNSGA-III algorithm is slightly worse than that of NSGA-II and NSGA-III to optimize ZDT4 and DTLZ1, the Spread of LCNSGA-III is prominent.The Spread of LCNSGA-III performs best to optimize the ZDT1, ZDT4, and DTLZ5 benchmark problems.And the Spread for LCNSGA-III does not differ much from those for the other five MOOPs except DTLZ1, which indicates that the LCNSGA-III algorithm can obtain good convergence performance for the eight benchmark functions.The Pareto optimal front obtained by LCNSGA-III to optimize the eight benchmark functions is shown in Figure 5 (ZDT1-4) and Figure 6 (ZDT6, DTLZ1-2, 5) (Grey points denote the true Pareto front while blue ones denote the obtained Pareto front), respectively.As depicted in Figures 5 and 6, the Pareto front obtained by LCNSGA-III to optimize the eight benchmark MOOPs can approximate the real Pareto front perfectly and the distributions of the Pareto optimal sets are uniform.As a result, the LCNSGA-III algorithm has prominent convergence performance and optimization ability compared with the other algorithms.

Experiments Design and Results
. In order to fully verify the performance and effectiveness of the LCNSGA-III algorithm in solving the MO-FOPID problem under multiworking conditions, the nonlinear model of PTGS with different controllers (PID and FOPID) and water heads (198m, 205m, and 210m) under no-load conditions is simulated on MATLAB environment.The frequency perturbation is set as 4% of the rated frequency and the simulation time is set as 50s.A total number of ten schemes have been designed to obtain the optimal parameters.The ten schemes can be divided into four categories as follows: (1) The backtracking search algorithm (BSA) [47] is exploited to optimize the parameters of PID controller under a single working condition of different water heads (198m, 205m, and 210m).The ITAE of PTGS under a certain working condition is selected as the objective function for (  ) = ∑ * ( − 0.5)     (2) The BSA algorithm is exploited to optimize the parameters of FOPID controller under a single working condition of different water heads (198m, 205m, and 210m).The experiments are simplified as S-198-FOPID, S-205-FOPID, and S-210-FOPID, respectively.The optimized control parameters are applied on the other two working conditions to test the adaptability of different schemes in tracking the dynamic responses of PTGS; (3) The NSGA-III algorithm is adopted to optimize the parameters of PID or FOPID controller under two working conditions (198m and 210m).The ITAE of PTGS at working conditions of 198m and 210m water head is selected as the objective functions for multiobjective optimization, and the compromise optimal solution among the Pareto optimal set is selected.The experiments are simplified as NSGA-III-PID and NSGA-III-FOPID for the PID and FOPID controller, respectively.The compromise optimal control parameters are applied on the working condition of 205m water head to test the adaptability of different schemes in tracking the dynamic responses of PTGS; (4) The LCNSGA-III algorithm is adopted to optimize the parameters of PID or FOPID controller under two working conditions (198m and 210m).The experiments are simplified as LCNSGA-III-PID and LCNSGA-III-FOPID for the PID and FOPID controller, respectively.The compromise optimal control parameters are applied on the working condition of 205m water head to test the adaptability of different schemes in tracking the dynamic responses of PTGS.
The parameters of PTGS are set as follows: the ranges of   ,   , and   of the PID and FOPID controllers are all set as [0, 15] and the ranges of  and  are all set as [0, 2].The parameters of the BSA algorithm are set as follows: the population size is set as 50, the number of iterations is 200, and the control parameter  is set as the default value.The parameters of the NSGA-III and LCNSGA-III algorithms are as follows: the population size is set as 50; the iteration number is 300; the crossover probability is 0.7; the mutation probability is 0.3.
Apart from ITAE, the integral of the ITSE [29], the stable time (ST) and the overshoot (OSO) is adopted to evaluate the performance of different schemes.The ITSE index is defined as follows: In this study, the best compromise solution is selected according to the subjective weighting method based on experts' preferences of weights.The optimal control parameters of the ten schemes for PTGS are shown in Table 4.The performance indices including ITAE, ITSE, ST, and OSO of the ten schemes under different working water heads are shown in Table 5.In Table 5, for scheme S-198-PID, the results for 198m water head are training results while those for 205m and 210m heads are testing results.For scheme S-205-PID, the Complexity  at complex operating conditions simultaneously using multiobjective optimization algorithms should be studied and investigated.LCNSGA-III-PID.Similar phenomenon also exists between the solutions obtained by S-198-FOPID and S-210-FOPID and the Pareto optimal solutions obtained by LCNSGA-III-FOPID, which demonstrates that the pursing for a single objective (high or low water head) is at the cost of the other objective (corresponding low or high water head).What is more, compared with the schemes for single working conditions, the adaptability and robustness of the controllers can be greatly improved using multiobjective optimization schemes.A set of Pareto optimal solutions are obtained using multiobjective optimization schemes; it is convenient for the operator to select the most appropriate control parameters when the working condition changes or much more attention should be paid to the extremely low water head or high water head working conditions.

Analysis of Controllers
The effectiveness of the developed LCNSGA-III algorithm has been verified using the eight test functions which have been described in Section 5. 1.In what follows the superiority of the developed LCNSGA-III algorithm is further demonstrated by applying it to PTGS.The Pareto fronts of PID and FOPID controllers optimized by the two algorithms have been illustrated in Figures 11(a 5 that the performances indices for the compromise Pareto optimal solution of LCNSGA-III are all smaller than NSGA-III, which further demonstrates the superiority of LCNSGA-III in optimizing the MO-FOPID problem.For example, the improved percentages of scheme LCNSGA-III-PID to scheme NSGA-III-PID are 62.9%, 50.0%, and 84.9% in terms of ITAE, ITSE, and OSO, respectively in the training stage (198m).The improved percentages in the testing stage (205m) are 46.1%,33.3%, and 91.3% in terms of ITAE, ITSE, and OSO, respectively.

Conclusions
In order to make PSUs adaptable to the changes of working environment and improve their control quality and stability, this study constructs a multiobjective optimization framework to optimize the FOPID controller for PTGS under multiworking conditions.An LCNSGA-III algorithm based on Latin-hypercube sampling and chaos theory is proposed to solve the MO-FOPID problem under multiworking  (1) The classical NSGA-III algorithm is improved using the Latin hypercube sampling-based initialization technique and the chaotic crossover and mutation operators.Experiments of the eight test functions show that the improvement strategies can effectively improve the convergence and diversity of the Pareto optimal front of the NSGA-III algorithm.
(2) Compared with the traditional PID controller, the FOPID can effectively suppress the frequency oscillation of PSUs in the "S" characteristic area running at middle or low working water heads and can enhance the dynamic response performance of PTGS.
(3) This study extends the single-objective optimization under single working conditions to the multiobjective optimization framework under multiworking conditions.The multiobjective framework has provided better dynamic performances than the single-objective optimization methods.The multiobjective implementation of FOPID makes it much more convenient for a operator to select the most appropriate control parameters for a certain working condition.
This work sets a basis for research on multiobjective optimization of a FOPID controller of PTGS under multiworking conditions.The extension of the single-objective optimization under a single working-condition to the multiobjective framework under complex working conditions can provide new control law and optimization algorithms for the optimal control of PTGS.What is more, the proposed LCNSGA-III algorithm can be easily extended to optimization problems in other fields of scientific research and industrial application.

Figure 1 :
Figure 1: The structure of the FOPID controller for pumped turbine governing system.

Figure 2 :
Figure 2: Three-dimensional surface of the characteristic curves of pump turbine.

) 4 . 4 . 3 . 4 . 5 .
The Flowchart of LCNSGA-III for MOOPs.Based on the above introduction of LCNSGA-III, the Latin hypercube sampling, and the chaotic crossover and mutation operators, the flowchart of LCNSGA-III for MOOPs is shown in Figure Implementation of LCNSGA-III for Solving the MO-FOPID Problem under Multiworking Conditions.The above LCNSGA-III algorithm is used to optimize the FOPID controller for PTGS under multiworking conditions.The schematic diagram of the MO-FOPID problem optimized by LCNSGA-III under multiworking conditions is shown in Figure 4.In Figure 4,   denotes the rotational speed (or frequency) and   denotes the given guide vane opening.

Figure 4 :
Figure 4: Schematic diagram of the MO-FOPID problem optimized by LCNSGA-III under multiworking conditions.

Figure 7 :
Figure 7: Unit frequency obtained at 198m working water head.

Figure 8 :
Figure 8: Unit frequency obtained at 205m working water head.

Figure 9 :
Figure 9: Unit frequency obtained at 210m working water head.

Figure 11 :
Figure 11: Comparison of Pareto fronts obtained by LCNSGA-III and NSGA-III for PID and FOPID controllers.

Table 1 :
Eight test functions for MOOPs.

Table 2 :
Comparison of GD for eight benchmark functions.

Table 3 :
Comparison of Spread for eight benchmark functions.

Table 4 :
Optimal control parameters for nonlinear PTGS using different schemes.

Table 5 :
Performance indices for nonlinear PTGS using different schemes at different water heads.