An Integrated Start-Up Method for Pumped Storage Units Based on a Novel Artiﬁcial Sheep Algorithm

: Pumped storage units (PSUs) are an important storage tool for power systems containing large-scale renewable energy, and the merit of rapid start-up enable PSUs to modulate and stabilize the power system. In this paper, PSU start-up strategies have been studied and a new integrated start-up method has been proposed for the purpose of achieving swift and smooth start-up. A two-phase closed-loop startup strategy, composed of switching Proportion Integration (PI) and Proportion Integration Differentiation (PID) controller is designed, and an integrated optimization scheme is proposed for a synchronous optimization of the parameters in the strategy. To enhance the optimization performance, a novel meta-heuristic called Artiﬁcial Sheep Algorithm (ASA) is proposed and applied to solve the optimization task after a sufﬁcient veriﬁcation with seven popular meta-heuristic algorithms and 13 typical benchmark functions. Simulation model has been built for a China’s PSU and comparative experiments are conducted to evaluate the proposed integrated method. Results show that the start-up performance could be signiﬁcantly improved on both indices on overshoot and start-up, and up to 34%-time consumption has been reduced under different working condition. The signiﬁcant improvements on start-up of PSU is interesting and meaning for further application on real unit.


Introduction
Nowadays, the increase of energy demand and the depletion of fossil fuel have necessitated the large-scale utilization of renewable energy (RE). More and more REs, like wind and solar power have access to power grid [1]. Due to the intermittent of the wind and solar power, the stability of the power grid that contains large capacity of RE has become a severe problem. Energy storage devices like batteries [2,3], flywheels [4] and pumped storage units (PSUs) [5] are indispensable in RE-connected power systems to mitigate this problem. Among these, PSUs might be the most economical and mature tool for power storage [6]. In addition, PSUs have good performance in peak load shaving, frequency regulation and emergency response in power systems [7]. In recent years, with the increasing scale of RE-integrated power systems, a considerable amount of work [8][9][10][11] has been carried out to study large-scale pumped storage facilities integrated with the grid-connected wind or solar power systems. In [8], joint operation of REs and energy storage devices like PSUs is studied in the energy and ancillary service market. In [9], the scheduling problem of a hybrid energy system containing intermittent solar power and pumped hydro storage (PHES) system has been investigated. In [10], a flexibility-based reserve scheduling method for energy storage system with PSUs has been developed to improve the flexibility of the power system. In [11], a model for transient stability analysis of a combined wind-pumped storage generation system is studied. Furthermore, this new generation model will be proposed and applied to solve different optimization problems, there is no specific algorithm which can solve all optimization problems. For optimization of start-up strategy, a powerful algorithm is essential to achieve the desired goal.
Flocks of sheep, which is a common swarm in nature, also show intelligence in social activities. In foraging ground, the strolling of sheep implies individual exploration, and the bellwether is the leading elite that affects the swarm movement. Many swarm intelligence (SI) techniques are inspired by foraging and search behaviors. However, there is no SI technique in the literature mimicking the Herding Effect, which is a well-known phenomenon in the sheep flock. Motivated by this interesting behavior, a new meta-heuristic optimization algorithm is proposed in this paper. The mathematical model of the social behavior of a sheep flock, as well as investigation of its abilities in solving benchmark problems, will be fully discussed before applying it in start-up strategy optimization of PSU.
Motivated by the above discussion, a new integrated start-up method for PSU is proposed in this paper, while a two phase's closed-loop startup strategy is designed, and optimization scheme is built for parameter optimization of the strategy. Moreover, in order to promote the optimization performance, new meta-heuristic algorithm is studied. This paper presents a novel meta-heuristic called Artificial Sheep Algorithm (ASA) based on social behaviors of sheep flock. The mathematically model of the social behavior of sheep flock, as well as investigation of its abilities in solving complicated problems, will be fully discussed.
The remaining part of this paper is organized as follows: Section 2 establishes a mathematical of Pump-Turbine Governing System which composed of PID controller, governor servo-mechanism, water diversion system, pump-turbine and generator. Section 3 introduces the theoretical knowledge of ASA. Section 4 delineates specific operations of traditional start-up and the proposed integrated start-up strategy. In Section 5, the performances of optimization schemes and integrated start-up strategy are analyzed by comparison experiment. The conclusions are summarized in the Section 6.

Modelling of Control System of PSU
The pump-turbine governing system (PTGS) is a main control system that undertakes the modulation of frequency (rotational speed) and power output. The simulation of PSU start-up is also related to PTGS. PTGS is composed of governor controller, governor servo-mechanism, water diversion system, pump-turbine and generator, as presented in Figure 1. In this paper, a method of characteristic (MOC) and an improved Suter transformation method have been studied for the modelling and simulation of PTGS [30,31].

Modelling of Water Diversion System
In this paper, the fundamental equations are applied to describe the water diversion pipelines, as shown below [32]:

Modelling of Water Diversion System
In this paper, the fundamental equations are applied to describe the water diversion pipelines, as shown below [32]: The details of all the symbols in these equations are given in the nomenclature. The method of characteristics (MOC) is applied to the partial differential equations (PDEs) as shown in Equations (1) and (2), and the partial differential terms associated with the flow velocity and the pressure are reduced to ordinary differential ones compatible with two characteristics lines C + and C − as shown in Figure 2.
Energies 2018, 11, 151 4 of 29 Continuity equation: The details of all the symbols in these equations are given in the nomenclature. The method of characteristics (MOC) is applied to the partial differential equations (PDEs) as shown in Equations (1) and (2), and the partial differential terms associated with the flow velocity and the pressure are reduced to ordinary differential ones compatible with two characteristics lines C + and C − as shown in Figure 2. The fixed-grid MOC requires that a common time step ( t Δ ) being used for the solution of the governing equations in all pipelines [33]. Consider L is one of the pipelines total length subdivided into equal sections N, each of which / r L N Δ = . If we start with known steady state conditions at t = 0, then we know Q and H at the N + 1 sections of the pipeline. If we specify the time interval t Δ defined as the / t r c Δ = Δ , the characteristic lines from the sections A and B intersect at P [34]. In these conditions, the final equations can be written in the following forms [32]:

ΔL ΔL
where coefficients a C gA c = In addition, the elastic water hammer effect is considered and different forms of pipelines, channels, and surge tanks are included.

Modelling of Pump-Turbine
Pump-turbines are the key component of PTGS, and the modeling of a pump turbine is based on its complete characteristic curves. The original performance curves of the pump-turbine unit are shown in Figure 3. The "S" area of characteristic curves presents an uneven distribution, crossing and aggregation phenomenon, which brings difficulties for modeling of pump-turbines. In order to overcome these problems, an improved Suter transformation method [35,36] have been proposed, and experiments have proved that the improved Suter transformation method could eliminate the uneven distribution, crossing and aggregation of "S" zone. Different processing methods for the curves can establish the pump-turbine linear model under different conditions as well as a nonlinear model by interpolation or fitting. The fixed-grid MOC requires that a common time step (∆t) being used for the solution of the governing equations in all pipelines [33]. Consider L is one of the pipelines total length subdivided into equal sections N, each of which ∆r = L/N. If we start with known steady state conditions at t = 0, then we know Q and H at the N + 1 sections of the pipeline. If we specify the time interval ∆t defined as the ∆t = ∆r/c, the characteristic lines from the sections A and B intersect at P [34]. In these conditions, the final equations can be written in the following forms [32]: where coefficients C a = gA/c, C f = f ∆t/2DA. In addition, the elastic water hammer effect is considered and different forms of pipelines, channels, and surge tanks are included.

Modelling of Pump-Turbine
Pump-turbines are the key component of PTGS, and the modeling of a pump turbine is based on its complete characteristic curves. The original performance curves of the pump-turbine unit are shown in Figure 3. The "S" area of characteristic curves presents an uneven distribution, crossing and aggregation phenomenon, which brings difficulties for modeling of pump-turbines. In order to overcome these problems, an improved Suter transformation method [35,36] have been proposed, and experiments have proved that the improved Suter transformation method could eliminate the uneven distribution, crossing and aggregation of "S" zone. Different processing methods for the curves can establish the pump-turbine linear model under different conditions as well as a nonlinear model by interpolation or fitting. By setting the unit parameters N11~Q11 and N11~M11 as a cross-ordinate, a series of twodimensional curves could be built to describe the complete characteristics of a pump turbine. While setting the gate opening y as reference variable, the flow characteristic curves and the torque characteristic curves of the complete characteristics could be established for a pump-turbine, as shown in Figure 3. In this figure, the range of guide vane opening is from 0° to 26°. From Figure 3, it can be noticed that the "S" area exists in characteristic curves. A multiple-valued problem occurs in interpolation calculating for pump-turbine modeling. In order to overcome this obstacle, an improved Suter transformation is used, as follows [36]: , k2 = 0.5~1.2, Cy = 0.1~0.3, Ch = 0.4~0.6. A transformed WH and WM curves based on the improved Suter transformation method have been shown in Figure 4. It is obvious that the improved Suter transformation method eliminate the "S" characteristics, the uneven distribution, cross and aggregation of the curves.  By setting the unit parameters N 11~Q11 and N 11~M11 as a cross-ordinate, a series of two-dimensional curves could be built to describe the complete characteristics of a pump turbine. While setting the gate opening y as reference variable, the flow characteristic curves and the torque characteristic curves of the complete characteristics could be established for a pump-turbine, as shown in Figure 3. In this figure, the range of guide vane opening is from 0 • to 26 • . From Figure 3, it can be noticed that the "S" area exists in characteristic curves. A multiple-valued problem occurs in interpolation calculating for pump-turbine modeling. In order to overcome this obstacle, an improved Suter transformation is used, as follows [36]: , k 2 = 0.5~1.2, C y = 0.1~0.3, C h = 0.4~0.6. A transformed WH and WM curves based on the improved Suter transformation method have been shown in Figure 4. It is obvious that the improved Suter transformation method eliminate the "S" characteristics, the uneven distribution, cross and aggregation of the curves. By setting the unit parameters N11~Q11 and N11~M11 as a cross-ordinate, a series of twodimensional curves could be built to describe the complete characteristics of a pump turbine. While setting the gate opening y as reference variable, the flow characteristic curves and the torque characteristic curves of the complete characteristics could be established for a pump-turbine, as shown in Figure 3. In this figure, the range of guide vane opening is from 0° to 26°. From Figure 3, it can be noticed that the "S" area exists in characteristic curves. A multiple-valued problem occurs in interpolation calculating for pump-turbine modeling. In order to overcome this obstacle, an improved Suter transformation is used, as follows [36]: where,  Figure 4. It is obvious that the improved Suter transformation method eliminate the "S" characteristics, the uneven distribution, cross and aggregation of the curves.   According to the above content, the calculation model of pump-turbine could be summarized as Equation (6): where equations h n+1 and m n+1 represent the interpolation of the characteristic curves of the pump turbine unit.

Modelling of Power System
This model calculates the unit running speed from the unbalance between the load moment M g and the shaft mechanical moment M t as defined in follows [37]: Under start-up operation, the value of M g is 0 and the corresponding equation can be considered as a special case of Equation (7), as follows:

Modelling of Pump-Turbine Governor
The pump-turbine governor contains a controller and servomechanism. A conventional PID controller is often used to eliminate the speed deviations from a reference speed in the later phases of a start-up process and under no-load conditions [38]. The transfer function of a PID controller is described as: G PID (s) = K p + K i s + K d s (9) and the output of the controller is shown as: where K p , K i and K d are the proportional gain, integral gain, differential gain of PID controller. u(s), a ref (s) and a(s) are the Laplace transform of controller output u, reference speed a ref and synchronous generator speed a. The servomechanism is the actuator of the speed governor which consists of main and assistant servomotors. The servomechanism is used to amplify the control signal and provide power to operate the guide vane of a pump-turbine. The structure of the servomechanism is shown as Figure 5. According to the above content, the calculation model of pump-turbine could be summarized as Equation (6) where equations hn+1 and mn+1 represent the interpolation of the characteristic curves of the pump turbine unit.

Modelling of Power System
This model calculates the unit running speed from the unbalance between the load moment Mg and the shaft mechanical moment Mt as defined in follows [37]: Under start-up operation, the value of Mg is 0 and the corresponding equation can be considered as a special case of Equation (7), as follows:

Modelling of Pump-Turbine Governor
The pump-turbine governor contains a controller and servomechanism. A conventional PID controller is often used to eliminate the speed deviations from a reference speed in the later phases of a start-up process and under no-load conditions [38]. The transfer function of a PID controller is described as: and the output of the controller is shown as: where Kp, Ki and Kd are the proportional gain, integral gain, differential gain of PID controller. u(s), aref (s) and a(s) are the Laplace transform of controller output u, reference speed aref and synchronous generator speed a. The servomechanism is the actuator of the speed governor which consists of main and assistant servomotors. The servomechanism is used to amplify the control signal and provide power to operate the guide vane of a pump-turbine. The structure of the servomechanism is shown as Figure 5.

Simulation of PTGS
While simulating the PTGS, the coupling relationships of water diversion system, pump-turbine and generator should be considered. Thus the iterative computation to calculate the water flow Qn+1

Simulation of PTGS
While simulating the PTGS, the coupling relationships of water diversion system, pump-turbine and generator should be considered. Thus the iterative computation to calculate the water flow Q n+1 and turbine speed N n+1 of the next time n + 1 for the simulation of PTGS before switching PID regulation is adopted, as the flow chart shown in Figure 6. and turbine speed Nn+1 of the next time n + 1 for the simulation of PTGS before switching PID regulation is adopted, as the flow chart shown in Figure 6.
Assuming the speed Nn+1,k=Nn and the flow Qn+1,k=Qn，intialize the variables ɛ1, ɛ2, λ and k max , k=1 Solve the equations of the water diversion system by characteristic method, calculate the flow Q n+1,k+1 .  Figure 6. The water flow and turbine speed calculation of the PTGS.

Artificial Sheep Algorithm
Optimization of parameters in the start-up strategy is an important task in this paper. A novel meta-heuristic called Artificial Sheep Algorithm (ASA) based on the social behaviors of sheep flocks is proposed and used to optimize the parameters. A series of test experiments are conducted to evaluate the performance of the proposed ASA.
Sheep flocks are always regarded as a scattered organization. Although they are loose, and can collide with each other in group activities, they have the habitual nature of following and moving towards the bellwether who is the strongest in the group. As sheep appear to blindly follow, people

Artificial Sheep Algorithm
Optimization of parameters in the start-up strategy is an important task in this paper. A novel meta-heuristic called Artificial Sheep Algorithm (ASA) based on the social behaviors of sheep flocks is proposed and used to optimize the parameters. A series of test experiments are conducted to evaluate the performance of the proposed ASA. Sheep flocks are always regarded as a scattered organization. Although they are loose, and can collide with each other in group activities, they have the habitual nature of following and moving towards the bellwether who is the strongest in the group. As sheep appear to blindly follow, people always associate them with conformity. Actually, the behavior of this social group has undergone ten million years of natural testing. From the perspective of SI, the loose individuals indicate random diversity, the blind herd behavior indicates high-efficiency positive feedback, and the bellwether leading implements the sharing of information in the whole group. This paper is inspired by the herd behavior of sheep flock and tries to simulate a novel meta-heuristic optimization algorithm with some artificial measures.
The social behaviors of sheep flock are attractive, and two main behaviors, namely the free strolling of individuals and strong leading of bellwether. In a sheep flock, the strongest sheep is called the bellwether, acting as the leader of the swarm. When fleeing from a predator or foraging, the individual will follow the bellwether by moving as close as possible. When strolling or playing, individuals are loose and always move randomly in their own local region.

Theoretical Knowledge of Algorithm
Consider a sheep flock with N sheep, the position of the ith sheep at specific time "t" is defined as: where x d i represents the position of ith sheep in the dth dimension, D is the dimension of the position. For solving a minimization problem: where B L = (b 1 l , . . . , b d l , . . . , b D l ) is the lower boundary of the searing space, and B U = (b 1 u , . . . , b d u , . . . , b D u ) is the upper boundary of the searing space. The objective function value of the ith agent at time "t" is expressed ).

Leading of bellwether
The influence of the bellwether is decisive. When the bellwether moves with a big stride, individuals will adjust their motion trajectory to follow the bellwether closely. The position of the bellwether should be recorded and inherited, and this position is denoted as X B (t) = X B d (t) 1×D . The influence of the bellwether acting on the ith sheep is expressed as bellwether vector, denoted by The bellwether vector that affects the movement of the ith agent, i = 1, . . . , N, is defined as: where δ i,d is the influence scope of the bellwether playing on the ith sheep on the dth dimension, , w is a dynamic weight that linearly decreased from 1 to 0 over the course of iterations. The coefficient c 1 is a random value whose center is 1, and its radius is determined by the parameter α, which is selected from [0, 1]. The α is a control parameter, which influences the consensus effect of the bellwether in defining the distance vector. When α tends to 1, it emphasizes the consensus influence of the bellwether; when α tends to 0, it enhances the stochastic components. The coefficient c 2 is a random dynamic number that is automatically generated, and its random range is linearly decreased over the course of iterations.

Individual strolling
Every individual of flock forages autonomously in a local area, and this behavior is called "self-awareness". The mathematical model of "individual strolling" to represent shelf-awareness of sheep individual in the process of foraging is proposed as follow. We define a location vector to denote self-awareness in foraging.
The shelf-awareness vector that affects the movement of the ith agent, i = 1, . . . , N, is defined as: where rand 2 is a random number generated from [0, 1], the term cos(2π · rand 1 ) is used to generate periodic random.
The vector X sel f i,d represents the self-driven behavior and local random search of a sheep, which consists of a series of nonlinear operations. The β parameter is a positive number that modulates the amplitude of the strolling step. As β increases, the sheep's jump steps increase exponentially and vice versa. As a result, this parameter controls the resolution of individual exploration. The value of β should be chosen according to the search scope of the optimization problem.
In a sheep flock, the direction of the sheep flock is determined by the lead of the bellwether and the autonomous foraging. Based on the discussion above, the movement of a sheep is affected by its self-awareness and the summoning of the bellwether. The individuals in artificial sheep flock will automatically update their positions as follows: where w is linearly decreased from 1 to 0 over the course of iterations and r 3 is random number generated in [0, 1].

Competition strategy
In ASA, competition mechanism is designed to keep the diversity of the flock. For a minimum problem, at a specific time "t", calculate the average value of N objective function values of the flock F t ave , and the minimal objective function value F t min . For the ith sheep, if the elimination criteria as following is satisfied: Then, the ith sheep X i is eliminated and reinitialized between [B L , B U ].

The Optimized Procedures of Algorithm
The whole idea and the specific architecture of ASA have been clearly presented in above paragraphs. The flow chart of ASA optimization is illustrated in Figure 7. Accordingly, the main steps of the ASA is summarized as follows.

•
Step 1: Initialization. Initialize locations X i (0) of sheep flock with N sheep in the solution space with boundaries [B L , B U ], and calculate the objective function value F 0 i = f (X i (0)), i = 1, . . . ,N; set the first sheep as the bellwether X B = X 1 (0), F B = F 0 1 ; set other control parameters: the initial scope coefficient of leading α and the modulation coefficient of strolling β; set the total number of iteration T, and the current number of iteration t = 0.

•
Step 5: Updating position of flock according to the Equation (15).

•
Step 6: Judge whether sheep need to be eliminated according to the Equation (16) and reinitialize the same number of new sheep between [B L , B U ].

•
Step 7: t = t + 1; if t > T max , end and output bellwether's position as the final solution; else, go to Step 2.
The ASA is inspired by the behavior of the sheep flock and the complete optimization mechanism of ASA algorithm has been established now. , i = 1,…,N, according to the Equation (14).

•
Step 5: Updating position of flock according to the Equation (15).

•
Step 6: Judge whether sheep need to be eliminated according to the Equation (16)  L U B B .

•
Step 7: t = t + 1; if t > Tmax, end and output bellwether's position as the final solution; else, go to Step 2.
The ASA is inspired by the behavior of the sheep flock and the complete optimization mechanism of ASA algorithm has been established now.  Figure 7. Flowchart of ASA optimization.

Experimental Study and Comparison Results on Benchmark Functions
To evaluate the performance of the proposed ASA, as well as other seven meta-heuristic optimization algorithms, i.e., Particle Swarm Optimization (PSO) [25], Differential Evolution (DE) [39], Ant Colony Optimization for continuous domain (ACOR) [40], Artificial Bee Colony (ABC) [41], Cuckoo Search (CS) [42], Gravitational Search Algorithm (GSA) [27], and Grey Wolf Optimizer (GWO) [43], were applied to solve 13 standard benchmark functions. The selected algorithms are representatives of evolutionary, physics-based and SI based meta-heuristics. In these algorithms, the newly developed meta-heuristics are included. The following experiments are conducted on MATLAB (R2016a, MathWorks, Natick, MA, USA).

Experimental Study and Comparison Results on Benchmark Functions
To evaluate the performance of the proposed ASA, as well as other seven meta-heuristic optimization algorithms, i.e., Particle Swarm Optimization (PSO) [25], Differential Evolution (DE) [39], Ant Colony Optimization for continuous domain (ACO R ) [40], Artificial Bee Colony (ABC) [41], Cuckoo Search (CS) [42], Gravitational Search Algorithm (GSA) [27], and Grey Wolf Optimizer (GWO) [43], were applied to solve 13 standard benchmark functions. The selected algorithms are representatives of evolutionary, physics-based and SI based meta-heuristics. In these algorithms, the newly developed meta-heuristics are included. The following experiments are conducted on MATLAB (R2016a, MathWorks, Natick, MA, USA).

Benchmark Functions
Tables 1 and 2 present the benchmark functions used in our experimental study. The employed benchmark functions can be divided into two groups: unimodal, multimodal test functions with high dimension. The unimodal test functions have only one global optimal value, so they are often used to test the ability of exploitation of algorithms, as shown in Table 1. The multimodal test functions with high dimension have many local optimal values expect for the global optimum, so they are most difficult to optimize, as shown in Table 2.
In these tables, D is the dimension of function, f opt is the minimum value of the function, and searching space is a subset of R D . The minimum value (f opt ) of the functions of Tables 1 and 2

Parameters of Algorithms
A set of fair parameters obtained by a simple trail-and-error procedure was adopted for ASA. For a fair comparison, the recommended parameters of PSO, DE, ACO R , ABC, CS, GSA and GWO were used to tackle these problems. In order to conduct a fair comparison, the maximum number of iteration was chosen to 1000 and the population size of all algorithms was set to 30 (i.e., N = 30) in all experiences in this section. Other parameters of the comparative algorithms were set as follows: • PSO: the inertia weight w was decreased linearly from 1 to 0.2. c 1 and c 2 are learning genes and chosen to be 2 in this section in this paper [44]. • ACO R : the number of ants used in an iteration m = 3, an archive size k = 50, the Locality of the search process, q = 0.05 and the speed of convergence, E = 0.85 [44]. • ABC: the number of onlooker bees, employed bees and food sources was 30 [44]. • CS: the discovery rate p = 0.25 [42]. • GSA: G 0 = 100, α = 20 [27].
• GWO: → a are linearly decreased from 2 to 0 over the course of iterations and r 1 , r 2 are random vectors in [0, 1] [43].

Comparison Results on Benchmark Functions
All experiences were repeated 30 times independently. The mean value, the best value, as well as the standard deviation, of the optimal objective function values of 30 runs, were reported for unimodal functions in Table 3 and for multimodal test functions with high dimension in Table 4.   The unimodal test functions are relatively simple to optimize because they only have one global optimal value. From Table 3, it is clear that the ASA has a very good performance for the unimodal test functions. The ASA algorithm shows obvious advantage over the competitors in the remainder of high dimension unimodal functions, except for F 4 and F 6 . The multimodal test functions with high dimension are relatively difficult to optimize because they have many local optimal values and one optimal global value. The ASA also has significant performance, as shown in Table 4. Generally, the ASA algorithm has shown an obvious advantage over the competitors on high dimensional benchmark functions F 1 -F 13 . Furthermore, optimization processes of the comparative algorithms were given in Figures 8 and 9. The values shown in these figures were the average of "best so far" in iterations achieved from 30 runs. The "best so far" is the best objective function value searched by the artificial sheep flock at certain iteration. The unimodal test functions are relatively simple to optimize because they only have one global optimal value. From Table 3, it is clear that the ASA has a very good performance for the unimodal test functions. The ASA algorithm shows obvious advantage over the competitors in the remainder of high dimension unimodal functions, except for F4 and F6. The multimodal test functions with high dimension are relatively difficult to optimize because they have many local optimal values and one optimal global value. The ASA also has significant performance, as shown in Table 4. Generally, the ASA algorithm has shown an obvious advantage over the competitors on high dimensional benchmark functions F1-F13. Furthermore, optimization processes of the comparative algorithms were given in Figures 8 and 9. The values shown in these figures were the average of "best so far" in iterations achieved from 30 runs. The "best so far" is the best objective function value searched by the artificial sheep flock at certain iteration.       Figure 9 shows the average iteration processes obtained by the eight methods on the multimodal functions F8-F13. From Figure 9, it is demonstrated that ASA performs the best compared with another seven methods on F8-F12. ASA has advantages over other algorithms, except for DE, in the optimization process of F13. From these figures, it was seen that the ASA had a faster converge speed compared with other approaches. In particular, ASA had a very good performance on high dimensional test functions.

Traditional Start-Up Strategies of a PSU
The traditional start-up method of a PSU contains two phases. In the first phase, the control system of the PSU is open-loop, in which the feedback signal of rotational speed of pump-turbine is not adopted and used for control. In this phase, the control system undertakes a kind of direct guide vane control (DGVC). After receiving the start-up command, governor will control the guide vane to tract a DGVC trajectory. As the rational speed of pump-turbine increase to a threshold value, the second phase starts. At that moment, the control system is closed-loop, while the feedback signal of pump-turbine speed is used for control. The system is then under a speed control mode. In the second  Figure 9 shows the average iteration processes obtained by the eight methods on the multimodal functions F 8 -F 13 . From Figure 9, it is demonstrated that ASA performs the best compared with another seven methods on F 8 -F 12 . ASA has advantages over other algorithms, except for DE, in the optimization process of F 13 . From these figures, it was seen that the ASA had a faster converge speed compared with other approaches. In particular, ASA had a very good performance on high dimensional test functions.

Traditional Start-Up Strategies of a PSU
The traditional start-up method of a PSU contains two phases. In the first phase, the control system of the PSU is open-loop, in which the feedback signal of rotational speed of pump-turbine is not adopted and used for control. In this phase, the control system undertakes a kind of direct guide vane control (DGVC). After receiving the start-up command, governor will control the guide vane to tract a DGVC trajectory. As the rational speed of pump-turbine increase to a threshold value, the second phase starts. At that moment, the control system is closed-loop, while the feedback signal of pump-turbine speed is used for control. The system is then under a speed control mode. In the second phase, PID control is adopted for the purpose of adjusting and maintaining the stability of pump-turbine speed referring to the rated speed.

One-stage DGVC
The guide value opens as the highest speed to no-load opening, while the rotational speed of unit will increase. As the speed reaches 90% of the rated value, a closed-loop PID control will be switched on to adjust and maintain the stability of speed for grid connection. The one-stage DGVC is illustrated in Figure 10.
phase, PID control is adopted for the purpose of adjusting and maintaining the stability of pumpturbine speed referring to the rated speed.

One-stage DGVC
The guide value opens as the highest speed to no-load opening, while the rotational speed of unit will increase. As the speed reaches 90% of the rated value, a closed-loop PID control will be switched on to adjust and maintain the stability of speed for grid connection. The one-stage DGVC is illustrated in Figure 10.

Two-stage DGVC
The guide vane opens rapidly to start opening (opening is about twice as much as the no-load opening) with the fastest speed, and keep the opening for a period, until the speed rises to a set value nc (nc is usually 60% of the rated speed). The guide vane opening is then adjusted to the no-load opening immediately. After the speed reaches 90% of the rated value, the PID control is switched on. The illustration of the two-stage DGVC is shown in Figure 11.

The Second Phase: PID Control
As mentioned above, following the open-loop direct guide vane control, the closed-loop PID control will be switched on when the pump-turbine speed reaches a threshold value, which is usually set as 90% of the rated speed. Though the PID controller is a very traditional design, it is still one of the favorite and most widely used controllers for many industrial process control applications. This

Two-stage DGVC
The guide vane opens rapidly to start opening (opening is about twice as much as the no-load opening) with the fastest speed, and keep the opening for a period, until the speed rises to a set value n c (n c is usually 60% of the rated speed). The guide vane opening is then adjusted to the no-load opening immediately. After the speed reaches 90% of the rated value, the PID control is switched on. The illustration of the two-stage DGVC is shown in Figure 11. phase, PID control is adopted for the purpose of adjusting and maintaining the stability of pumpturbine speed referring to the rated speed.

One-stage DGVC
The guide value opens as the highest speed to no-load opening, while the rotational speed of unit will increase. As the speed reaches 90% of the rated value, a closed-loop PID control will be switched on to adjust and maintain the stability of speed for grid connection. The one-stage DGVC is illustrated in Figure 10.

Two-stage DGVC
The guide vane opens rapidly to start opening (opening is about twice as much as the no-load opening) with the fastest speed, and keep the opening for a period, until the speed rises to a set value nc (nc is usually 60% of the rated speed). The guide vane opening is then adjusted to the no-load opening immediately. After the speed reaches 90% of the rated value, the PID control is switched on. The illustration of the two-stage DGVC is shown in Figure 11.

The Second Phase: PID Control
As mentioned above, following the open-loop direct guide vane control, the closed-loop PID control will be switched on when the pump-turbine speed reaches a threshold value, which is usually set as 90% of the rated speed. Though the PID controller is a very traditional design, it is still one of the favorite and most widely used controllers for many industrial process control applications. This

The Second Phase: PID Control
As mentioned above, following the open-loop direct guide vane control, the closed-loop PID control will be switched on when the pump-turbine speed reaches a threshold value, which is usually set as 90% of the rated speed. Though the PID controller is a very traditional design, it is still one of the favorite and most widely used controllers for many industrial process control applications. This is due to its simple structure, satisfactory control effect and acceptable robustness [45,46]. The PID controller is easier to understand due to intuitive simplicity of the algorithm and simple meaning of its tuning parameters proportional (K p ), integral (K i ) and derivative (K d ).
In the following, the start-up strategy composed of open-loop one-stage DGVC and closed-loop PID control is named the one-stage start-up strategy. The start-up strategy composed of open-loop two-stage DGVC and closed-loop PID control is named the two-stage start-up strategy.

Optimization Variables for Start-Up Strategies
In order to improve the start-up performance, the parameters in the start-up strategies could be optimized. Traditionally, PID control parameters are often optimized to enhance the control performance for a quick and smooth PSU start-up. The tuning parameter vector for a PID controller is The parameters of the DGVC could also be optimized, though this has seldom been researched. For a one-stage DGVC, the maximum value y c and the rising slope of guide vane k c can be used as parameters to be optimized, as θ 2 = [y c , k c ]. The parameters to be optimized of a two-stage DGVC are the maximum value of the guide vane opening y c1 , the minimum value after guide vane closure y c2 and the set frequency f c when the guide vane began to decrease, as θ 3 = [y c1 , y c2 , f c ].
In this paper, optimization experiments of the traditional two-phase start-up strategies are grouped into two categories: Scheme A: optimization of PID controller. In practice, parameters of the DGVC are set as a rule of thumb. The PID parameter vector θ is optimized. In these categories, one-stage DGVC and two-stage DGVC can be selected corresponding.
Scheme B: optimization of DGVC. DGVC parameters are optimized, the parameters of PID controller are kept the same as the optimized parameters of Scheme A.

An Integrated Start-Up Method
In order to improve the performance of the start-up process of PSUs, a new integrated start-up strategy and the relevant optimization method are proposed. This strategy is composed of two phases, in which a closed-loop PI control, instead of the open-loop GVO control, is applied in the first phase to open GVO, and a closed-loop PID control is switched on once the speed reaches the threshold. In the first phase, a closed-loop PI control is conducted to open the GVO to accelerate the rotational speed, while the control target is to keep the ratio of differential speed and speed deviation to be a constant: where d∆a/dt is the speed differential speed; ∆a is deviation of speed. When the speed is close to the rated speed (about 98% of the rated speed), the governor will switch to the PID controller to track the rated speed. Parameters of the integrated start-up strategy can be optimized synchronously. Parameters of the first phase include parameters K p1 , K i1 , of PI controller, and the constant C in Equation (17). Parameters of the second phase are that of PID controller, K p , K i, and K d . Therefore, the optimization vector is θ 4 = K p1 , K i1 , C, K p , K i , K d .

Objective Function
The objective function for start-up strategy optimization is an integral performance index to evaluate the control performances, including overshoot and start-up time. The Integral Time Absolute Error (ITAE) [47] index is used as the objective function as defined by: where θ is parameter vector to be optimized, k is the sample number, N s is the number of samples, T(k) is the sample time, a(k) is the relative value of unit speed. In order to assess the performance of the proposed integrated start-up methods quantitatively, three indices about overshoot, start-up time and steady-state error are presented: a max is the relative value of the peak speed, t s called the start-up time is the corresponding time when the relative value of the speed reaches rated relative speed 1 and the steady-state error means the deviation of the actual rotational unit speed and rated relative speed 1 that can't be eliminate. Here, considering the constraints of the experiment, we choose the steady-state error at 100 s. The specific description is shown in Figure 12. where θ is parameter vector to be optimized, k is the sample number, Ns is the number of samples, T(k) is the sample time, a(k) is the relative value of unit speed. In order to assess the performance of the proposed integrated start-up methods quantitatively, three indices about overshoot, start-up time and steady-state error are presented: (19) Steady-state error = (20) amax is the relative value of the peak speed, ts called the start-up time is the corresponding time when the relative value of the speed reaches rated relative speed 1 and the steady-state error means the deviation of the actual rotational unit speed and rated relative speed 1 that can't be eliminate. Here, considering the constraints of the experiment, we choose the steady-state error at 100 s. The specific description is shown in Figure 12.

Procedures
In the proposed method, the ASA is used to optimize the objective function, while optimization parameters are elements in θ = [ , , , , , ]. As for optimization of traditional start-up strategies, i.e., Scheme A and Scheme B, the objective function is also selected as the ITAE function, while the optimization vectors are θ , θ , θ , respectively. Take optimization of the integrated startup strategy as an example, the detail flow chart of start-up strategy optimization based on the ASA is shown in Figure 13. To illuminate the procedures briefly, the main steps are summarized as follows.

•
Step 1: Initialization. Initialize locations Step 2.1: Decode Xi, and get the control parameters of the first phase Kp1, Ki1, and C, and parameters of second phase Kp, Ki, Kd.
Step 2.2: Set these parameters for controllers and start the established PTGS simulation plant to simulate the start-up process of PSU; sample and record system outputs including the unit speed a(k), opening of guide vanes y(k).
Step 2.3: Calculate the objective function value by Equation (18), and let

•
Step 3: Implement other iterative procedures of ASA as described in Section 3.2;

Procedures
In the proposed method, the ASA is used to optimize the objective function, while optimization parameters are elements in θ 4 = K p1 , K i1 , C, K p , K i , K d . As for optimization of traditional start-up strategies, i.e., Scheme A and Scheme B, the objective function is also selected as the ITAE function, while the optimization vectors are θ 1 , θ 2 , θ 3 , respectively. Take optimization of the integrated start-up strategy as an example, the detail flow chart of start-up strategy optimization based on the ASA is shown in Figure 13. To illuminate the procedures briefly, the main steps are summarized as follows.

•
Step 1: Initialization. Initialize locations X i (0) of a sheep flock with N sheep in the solution space with boundaries [B L , B U ], set the first sheep as the bellwether X B = X 1 (0), F B = F 0 1 ; set other control parameters: the initial scope coefficient of leading α and the modulation coefficient of strolling β; set the total number of iteration T max , and the current number of iteration t = 0.

•
Step 2: Objective function calculation. For the ith (i = 1, .., N) agent, the objective function value F t i is calculated as three steps: Step 2.1: Decode X i , and get the control parameters of the first phase K p1 , K i1 , and C, and parameters of second phase K p , K i , K d .
Step 2.2: Set these parameters for controllers and start the established PTGS simulation plant to simulate the start-up process of PSU; sample and record system outputs including the unit speed a(k), opening of guide vanes y(k).
Step 2.3: Calculate the objective function value by Equation (18), and let F t i = F ITAE .

•
Step 3: Implement other iterative procedures of ASA as described in Section 3.2;

•
Step 4: Set, t = t + 1; if t > T max , stop simulation and output bellwether's position as the final solution; or else, go to Step 2. Objective function calculation Figure 13. The detail flow chart of parameters optimization.

Experiments
To verify the effectiveness of the integrated intelligent start-up strategy, a PSHP in Jiangxi Province of China has been investigated and used as the study target. The structure of the water diversion system of the PSHP is shown in Figure 14. In this study, the following condition is considered: unit 1 is working and unit 2 stops. The control system of the working PSU is simulated in MATLAB (R2016a, MathWorks, Natick, MA, USA).

Experiments
To verify the effectiveness of the integrated intelligent start-up strategy, a PSHP in Jiangxi Province of China has been investigated and used as the study target. The structure of the water diversion system of the PSHP is shown in Figure 14. In this study, the following condition is considered: unit 1 is working and unit 2 stops. The control system of the working PSU is simulated in MATLAB (R2016a, MathWorks, Natick, MA, USA). Objective function calculation Figure 13. The detail flow chart of parameters optimization.

Experiments
To verify the effectiveness of the integrated intelligent start-up strategy, a PSHP in Jiangxi Province of China has been investigated and used as the study target. The structure of the water diversion system of the PSHP is shown in Figure 14. In this study, the following condition is considered: unit 1 is working and unit 2 stops. The control system of the working PSU is simulated in MATLAB (R2016a, MathWorks, Natick, MA, USA).

Model Parameters
The values of initial and specified parameters of the simulated PTGS are listed in Table 5, while t s is sampling time. Experiments are conducted under same water head condition which the water level of upper reservoir is 735.45 m and downstream reservoir is 181 m. Table 5. Parameters of the simulated model of PTGS.

Components Value
Suter transformation k 1 = 10 k 2 = 0.9 C y = 0.2 C h = 0.5 Generator J = 96.84 Servo-mechanism T yB = 0.05 T y = 0.3 k 0 = 1 Simulation setting k max = 20 t s = 0.02 s As discussed in Section 4, different PSU start-up strategies, including traditional start-up strategies and the proposed one, are optimized by applying ASA. In the following experiments, the parameters of ASA is set as: the initial scope coefficient of leading α = 0 and the modulation coefficient of strolling β = 1, the population size is 30, the total number of iteration T max = 200. The boundaries [B L , B U ] of each optimized parameter are listed in Table 6. Initial DGVC parameters are selected by experience as some popular rules, which are presented in Table 7. Table 6. The boundaries of each optimized parameter.

Parameters Boundary
Value  Table 7. DGVC parameters selected by experience.

Comparative Analysis of Traditional Start-Up Strategies
In this part, experiments on optimization of the traditional start-up strategies are conducted. At first, optimization experiments on Scheme A are executed, while one-stage and two-stage DGVC is applied, respectively. In this optimization scheme, the parameters of the first phase DGVC are given by experience, as shown in Table 7, and the second phase PID controller is optimized by ASA.
The results of optimized PID controller parameters are listed in Table 8. The dynamic processes of one-stage and two-stage start-up strategies are shown in Figures 15 and 16, where the rotational speed and guide vane opening curves are presented respectively. The performance indices calculated by Equations (19) and (20), including start-up time, overshoot and steady-state error, are listed in Table 9.
From the results in Figure 5 and Table 9, it is found that one-stage strategy achieves a better performance on overshoot at the price of a longer start-up time. Compared to the one-stage start-up strategy, the two-stage strategy possesses a quicker rotational speed rise ratio. From Figure 6, the change process of guide vane opening of one-stage strategy has the same slope with the two-stage strategy to make the rotational speed rises quickly within the first 4.5 s, but the two-stage start-up strategy has a larger guide vane opening as time elapses, therefore, less time is consumed for the rotational speed to reach 90% of its rated value. strategy has a larger guide vane opening as time elapses, therefore, less time is consumed for the rotational speed to reach 90% of its rated value.    Subsequently, one-stage DGVC and two-stage DGVC are compared in experiments on optimization Scheme B, while the parameters of the PID controller are kept the same as in the previous results on Scheme A. In this scheme, the parameters of DGVC are optimized by the ASA.
The optimized DGVC parameters are listed in Table 10. The dynamic processes of one-stage and two-stage start-up strategies are shown in Figures 17 and 18, where the rotational speed and guide vane opening curves are presented. The performance indices are listed in Table 11. From these results, a similar conclusion is drawn that a smaller overshoot is obtained with the one-stage DGVC strategy and a shorter start-up time is achieved with the two-stage DGVC strategy. strategy has a larger guide vane opening as time elapses, therefore, less time is consumed for the rotational speed to reach 90% of its rated value.    Subsequently, one-stage DGVC and two-stage DGVC are compared in experiments on optimization Scheme B, while the parameters of the PID controller are kept the same as in the previous results on Scheme A. In this scheme, the parameters of DGVC are optimized by the ASA.
The optimized DGVC parameters are listed in Table 10. The dynamic processes of one-stage and two-stage start-up strategies are shown in Figures 17 and 18, where the rotational speed and guide vane opening curves are presented. The performance indices are listed in Table 11. From these results, a similar conclusion is drawn that a smaller overshoot is obtained with the one-stage DGVC strategy and a shorter start-up time is achieved with the two-stage DGVC strategy.   Subsequently, one-stage DGVC and two-stage DGVC are compared in experiments on optimization Scheme B, while the parameters of the PID controller are kept the same as in the previous results on Scheme A. In this scheme, the parameters of DGVC are optimized by the ASA.
The optimized DGVC parameters are listed in Table 10. The dynamic processes of one-stage and two-stage start-up strategies are shown in Figures 17 and 18, where the rotational speed and guide vane opening curves are presented. The performance indices are listed in Table 11. From these results, a similar conclusion is drawn that a smaller overshoot is obtained with the one-stage DGVC strategy and a shorter start-up time is achieved with the two-stage DGVC strategy.     Some interesting facets are worth discussing. The guide vane responses are similar to those in scheme A. However, in zone A of Figure 18, the guide vane opening doesn't remain invariant as a horizontal line like in Figure 16. The reason is that the rotational speed increases fast as the guide vane opening continues increasing in this case. The relative value of the speed has risen to 0.64 before guide vane opening reaches the inflection point yc1 = 0.413 to become horizontal. Moreover, point b is the moment the PID controller is switched on and there exists a tiny fluctuation.     Some interesting facets are worth discussing. The guide vane responses are similar to those in scheme A. However, in zone A of Figure 18, the guide vane opening doesn't remain invariant as a horizontal line like in Figure 16. The reason is that the rotational speed increases fast as the guide vane opening continues increasing in this case. The relative value of the speed has risen to 0.64 before guide vane opening reaches the inflection point yc1 = 0.413 to become horizontal. Moreover, point b is the moment the PID controller is switched on and there exists a tiny fluctuation.  Some interesting facets are worth discussing. The guide vane responses are similar to those in scheme A. However, in zone A of Figure 18, the guide vane opening doesn't remain invariant as a horizontal line like in Figure 16. The reason is that the rotational speed increases fast as the guide vane opening continues increasing in this case. The relative value of the speed has risen to 0.64 before guide vane opening reaches the inflection point y c1 = 0.413 to become horizontal. Moreover, point b is the moment the PID controller is switched on and there exists a tiny fluctuation.

Results of the Integrated Start-Up Method
In order to verify the superiority of the proposed integrated start-up method, a comparison of different start-up strategies is made. At first, the results of the integrated start-up method are presented, while the optimized parameters of the integrated method are listed in Table 12, and the corresponding dynamic responses of start-up process are illuminated in Figure 19. From these results, it is obvious that the start-up process achieved by the integrated method shows significant performance improvement, with a short start-up time and negligible overshoot.

Results of the Integrated Start-Up Method
In order to verify the superiority of the proposed integrated start-up method, a comparison of different start-up strategies is made. At first, the results of the integrated start-up method are presented, while the optimized parameters of the integrated method are listed in Table 12, and the corresponding dynamic responses of start-up process are illuminated in Figure 19. From these results, it is obvious that the start-up process achieved by the integrated method shows significant performance improvement, with a short start-up time and negligible overshoot.  The results of different strategies are presented in Table 13 and Figure 20, where the indices of overshoot, start-up time and steady-state error are shown, and the rotational speed curves of startup process are compared. From these results, it is easy to find out that the start-up time obtained by integrated start-up strategy is shorter than with other methods. This is especially important from the practical implementation point of view because it is a key index to evaluate the ability of quick connection to the grid. Because the guide vane opening is under closed-loop control in the whole start-up process, the control system has a better dynamic quality. Indices of overshoot and steadystate error obtained by the integrated method are still satisfactory, although the overshoot obtained by the integrated model is a little bigger than one-stage start-up strategy.  The results of different strategies are presented in Table 13 and Figure 20, where the indices of overshoot, start-up time and steady-state error are shown, and the rotational speed curves of start-up process are compared. From these results, it is easy to find out that the start-up time obtained by integrated start-up strategy is shorter than with other methods. This is especially important from the practical implementation point of view because it is a key index to evaluate the ability of quick connection to the grid. Because the guide vane opening is under closed-loop control in the whole start-up process, the control system has a better dynamic quality. Indices of overshoot and steady-state error obtained by the integrated method are still satisfactory, although the overshoot obtained by the integrated model is a little bigger than one-stage start-up strategy.

Performance under Different Initial Water Head Condition
In practical operation, the water head of a PSU might vary frequently. Hence other two initial conditions (T2 and T3) with different water levels in the upper and downstream reservoir have been adopted to test the performance of the discussed start-up methods. Three initial working conditions are listed in Table 14. Traditional start-up strategies, including the one-stage strategy and two-stage strategy are optimized in Scheme B. The performance indices of different shut-up methods under various head conditions are presented in Table 15. Compared with the traditional start-up strategies, the proposed method has achieved the best performance on all indices overall. Especially, the start-up time has been significantly shortened by the proposed method under all working conditions.

Performance under Different Initial Water Head Condition
In practical operation, the water head of a PSU might vary frequently. Hence other two initial conditions (T2 and T3) with different water levels in the upper and downstream reservoir have been adopted to test the performance of the discussed start-up methods. Three initial working conditions are listed in Table 14. Traditional start-up strategies, including the one-stage strategy and two-stage strategy are optimized in Scheme B. The performance indices of different shut-up methods under various head conditions are presented in Table 15. Compared with the traditional start-up strategies, the proposed method has achieved the best performance on all indices overall. Especially, the start-up time has been significantly shortened by the proposed method under all working conditions. In order to show the results more explicitly, bar graphs of the indices of overshoot and start-up time are exhibited in Figures 21 and 22. From these figures, it's clearly shown that a smooth and swift start-up of PSU could be realized by applying the proposed integrated start-up strategy, with an extremely low overshoot and a short start-up time. Compared with one-stage strategy, the proposed method has shortened the start-up time 22.8%, 34.4%, 26.9% under the three working conditions. Energies 2018, 11, 151 25 of 29 In order to show the results more explicitly, bar graphs of the indices of overshoot and start-up time are exhibited in Figures 21 and 22. From these figures, it's clearly shown that a smooth and swift start-up of PSU could be realized by applying the proposed integrated start-up strategy, with an extremely low overshoot and a short start-up time. Compared with one-stage strategy, the proposed method has shortened the start-up time 22.8%, 34.4%, 26.9% under the three working conditions.

Conclusions
In order to improve the start-up process indices, a new integrated start-up method is proposed. In this method, a two-phase closed-loop start-up strategy is designed, while PI control is applied to quickly open the guide vane with a specially designed control target in the first phase, and PID control is adopted to adjust the rotational speed of a PSU to track the rated value in the second phase. What's more, an integrated optimization method is used to synchronously tune the parameters of the PI and PID controllers.
To solve this complicated optimization problem, the ASA, a novel meta-heuristic that mimics the social behaviors of sheep, has been proposed and fully verified by comparing it with seven popular meta-heuristics on 13 typical benchmark functions. The results show that the ASA significantly outperforms all competitors for high dimension numerical functions.
Based on the mathematical model of PSUs, the control system, PTGS, has been modeled and simulated in MATLAB, and simulation experiments on different start-up strategies have been conducted. The proposed method has been compared with the traditional one-stage DGVC, and twostage DGVC strategies with different optimization schemes. The experimental results have revealed that the proposed integrated start-up method shows great advantages in terms of performance indices, i.e., overshoot, start-up time, compared with the traditional start-up strategies. The start-up time could be improved by as much as 34%, while maintaining the overshoot under a low level. The significant improvements on these key indices is interesting and merit further application in real PSUs. In order to show the results more explicitly, bar graphs of the indices of overshoot and start-up time are exhibited in Figures 21 and 22. From these figures, it's clearly shown that a smooth and swift start-up of PSU could be realized by applying the proposed integrated start-up strategy, with an extremely low overshoot and a short start-up time. Compared with one-stage strategy, the proposed method has shortened the start-up time 22.8%, 34.4%, 26.9% under the three working conditions.

Conclusions
In order to improve the start-up process indices, a new integrated start-up method is proposed. In this method, a two-phase closed-loop start-up strategy is designed, while PI control is applied to quickly open the guide vane with a specially designed control target in the first phase, and PID control is adopted to adjust the rotational speed of a PSU to track the rated value in the second phase. What's more, an integrated optimization method is used to synchronously tune the parameters of the PI and PID controllers.
To solve this complicated optimization problem, the ASA, a novel meta-heuristic that mimics the social behaviors of sheep, has been proposed and fully verified by comparing it with seven popular meta-heuristics on 13 typical benchmark functions. The results show that the ASA significantly outperforms all competitors for high dimension numerical functions.
Based on the mathematical model of PSUs, the control system, PTGS, has been modeled and simulated in MATLAB, and simulation experiments on different start-up strategies have been conducted. The proposed method has been compared with the traditional one-stage DGVC, and twostage DGVC strategies with different optimization schemes. The experimental results have revealed that the proposed integrated start-up method shows great advantages in terms of performance indices, i.e., overshoot, start-up time, compared with the traditional start-up strategies. The start-up time could be improved by as much as 34%, while maintaining the overshoot under a low level. The significant improvements on these key indices is interesting and merit further application in real PSUs.

Conclusions
In order to improve the start-up process indices, a new integrated start-up method is proposed. In this method, a two-phase closed-loop start-up strategy is designed, while PI control is applied to quickly open the guide vane with a specially designed control target in the first phase, and PID control is adopted to adjust the rotational speed of a PSU to track the rated value in the second phase. What's more, an integrated optimization method is used to synchronously tune the parameters of the PI and PID controllers.
To solve this complicated optimization problem, the ASA, a novel meta-heuristic that mimics the social behaviors of sheep, has been proposed and fully verified by comparing it with seven popular meta-heuristics on 13 typical benchmark functions. The results show that the ASA significantly outperforms all competitors for high dimension numerical functions.
Based on the mathematical model of PSUs, the control system, PTGS, has been modeled and simulated in MATLAB, and simulation experiments on different start-up strategies have been conducted. The proposed method has been compared with the traditional one-stage DGVC, and two-stage DGVC strategies with different optimization schemes. The experimental results have revealed that the proposed integrated start-up method shows great advantages in terms of performance indices, i.e., overshoot, start-up time, compared with the traditional start-up strategies. The start-up time could be improved by as much as 34%, while maintaining the overshoot under a low level. The significant improvements on these key indices is interesting and merit further application in real PSUs.