The Improved Binary-Real Coded Shuffled Frog Leaping Algorithm for Solving Short-Term Hydropower Generation Scheduling Problem in Large Hydropower Station

The short-termhydro generation scheduling (STHGS) decomposed into unit commitment (UC) and economic load dispatch (ELD) subproblems is complicated problemwith integer optimization,which has characteristics of high dimension, nonlinear and complex hydraulic and electrical constraints. In this study, the improved binary-real coded shuffled frog leaping algorithm (IBR-SFLA) is proposed to effectively solveUC and ELD subproblems, respectively. For IB-SFLA, the new grouping strategy is applied to overcome the grouping shortage of SFLA, and modified search strategies for each type of frog subpopulation based on normal cloud model (NCM) and chaotic theory are introduced to enhance search performance.The initialization strategywith chaos theory and adaptive frog activation mechanism are presented to strengthen performance of IR-SFLA on ELD subproblem. Furthermore, to solve ELD subproblem, the optimal economic operation table is formed using IR-SFLA and invoked fromdatabase.Moreover, reserve capacity supplement and repair, and minimum on and off time repairing strategies are applied to handle complex constraints in STHGS. Finally, the coupled external and internal model corresponding to UC and ELD subproblems is established and applied to solve STHGS problem inThree Gorges hydropower station. Simulation results obtained from IBR-SFLA are better than other compared algorithms with less water consumption. In conclusion, to solve STHGS optimization problem, the proposed IBR-SFLA presents outstanding performance on solution precision and convergence speed compared to traditional SFLA effectively and outperforms the rivals to get higher precision solution with improving the utilization rate of waterpower resources.


Introduction
Short-term hydro generation scheduling (STHGS), i.e., economic operation in hydropower station, including space and time optimization, is a typical dual optimization, multidimensional, noncontinuous, nonlinear problem with numerous complex hydraulic and electrical constraints.STHGS aims to find a profitable way to decrease water consumption and the cost of on and off transition between hydropower units with satisfying the power generation demand.The STHGS problem can be decomposed into two subproblems: the space optimization, namely, economic load dispatch (ELD), and the time optimization, the hydropower unit commitment (UC).The ELD subproblem attempts to reasonably determine the power output of each unit so as to minimize water consumption [1].Efficiency, under different operation conditions, corresponding to each type of units is discrepant especially the large capacity unit.In addition, to keep unit operational within stable operation zone is beneficial to prolong the life of units.Therefore, to effectively distribute output for a startup unit to burden and focus energy on improving integrated units operation efficiency in stable operating zone are of vital importance in STHGS problem.The unit commitment (UC), aiming to pick up a wise startup and shutdown schedule of units during operational period, is 0-1 binary coded subproblem [2].Along with the power generation, a large amount of water will be consumed in water turbine.Moreover, part of the water will be lost due to state change between startup and shutdown and frequent on and off transition is greatly harmful to service life of units [3].Thus, a rational and highly efficient unit commitment can steeply minimize water consumption and simultaneously 2 Mathematical Problems in Engineering improve dependability of safe operation in hydropower system.In summary, to accurately describe and find reasonable unit commitment and economic load distribution can not only save the water resources and operation cost but also improve the turbine efficiency and service life.In order to efficiently utilize hydropower resources in reservoir and save limited water resources, it is meaningful and significant to build the model of STHGS and acquire feasible solution with effective method.
During the past few years, a number of traditionally mathematical methods have been put forward to handle STHGS problem.The traditional methods are focused on branch and bound algorithm [4] (Cohen and Yoshimura 1987), linear programming (LP) [5] (Jabr et al. 2000), mixed integer nonlinear programming [6] (Catalao J 2011), Lagrangian relaxation (LR) [7][8][9][10] (Oliveira et al. 1992; Ohishi T, W. Ongsakul; N. Petcharaks; Sum-Im et al. 2014; Cheng-C-P), quadratic programming (QP) [11] (Granelli and Montagna 2000), and dynamic programming (DP) [12][13][14][15] (Allen and Bridgeman 1986; Siu et al. 2001;Yi et al. 2003; Cheng et al. 2009), which have met with different degrees of success to STHGS.Nevertheless, methods mentioned have defects on mechanism, stability, and efficiency.Lagrangian multiplier is introduced into the objective function to structure dual problem to simplify calculation in the LR method; however, it is still difficult for coupling various constraints and finding a suitable multiplier.DP, though it can obtain the best solution theoretically, tends to face dimension disaster and costs plenty of computation time, which limits the application in dealing with STHGS issue.Therefore, the approximation of dynamic programming (ADP) algorithm is put forward to solve the problem of dimension disaster and overcomes these shortcomings to an extent [16] (Wei et al. 2014).On the whole, traditional methods encountered numerous obstacles when facing large scale power systems with huge installed capacity.
To overcome shortage of traditional methods, along with high-speed computer technology development, intelligent algorithms, i.e., heuristic stochastic methods, have been developed rapidly to exploit application in managing STHGS problem.These methods mainly include evolutionary programming (EP) [17] (Christober et al. 2011;), differential evolution (DE) [18,19] [22] (Bakhta N et al. 2013), artificial neural network [23] (Sasaki H et al. 1992); particle swarm optimization (PSO) [24][25][26] (Yuan et al. 2009; Kumar et al. 2011;Zhao (2006)), chaotic optimization algorithm (COA) [27] (Cai et al. 2007), immune algorithm [28] (Wang and Yang 2010), ant colony optimization (ACO) [29,30] (Vaisakh K and Srinivas L (2011); Chen et al. 2015), gravitational search algorithm (GSA) [31][32][33][34] Li et al. 2015), and fireworks algorithm (FA) [35](K.Srikanth et al. 2016).Owing to flexible search strategies and powerful convergence capability, these algorithms can obtain higher quality solution with accepted optimization time to a certain extent in comparison with traditional methods.Furthermore, the mixed algorithm which combines superiority of several algorithms is introduced to accelerate computational speed and boost search ability for global optimum.However, intelligent heuristic stochastic methods also cause defects and imperfect constituent parts.The improved immune algorithm is applied to optimize the hydro generation scheduling, but the precision and efficiency need to further enhance.Hu et al. (2012) combined ant colony algorithm with equal incremental method; however, the operation efficiency drastically slows down with the increase in the number of units and parameter configuration may bring premature convergence trouble [36].Wu (2015) adopts the measure of discrete particle swarm optimization (DPSO) with relevant improvement; however, global convergence is inefficient which makes it incapable of obtaining convincible solutions [37].The gravitational search algorithm (GSA) proposed by Shen introduces the chaos theory to enhance local search ability, but the overall efficiency and robustness remain to be strengthened [31,32].In summary, the intelligent heuristic stochastic algorithms mentioned present outstanding performance on computational accuracy, stability, and speed compared to that in traditional method.However, few methods may encounter premature convergence with limitations in local and global search strategies, making it difficult to find theoretical optimal solution.Therefore, it is significant and challengeable to develop effective method to realize more efficient and economical hydro generation scheduling.
The shuffled frog leaping algorithm (SFLA), a stochastic global optimization method, is inspired by predation behavior of frogs in wetland and first brought forward by Eusuff and Lansey in 2006 [38].It is found that SFLA can guide the search toward global optimum direction by means of the global information exchange and internal communication mechanism [39].SFLA, which is easy and convenient to code and has less control parameters, is verified to be compatible with handling comprehensive optimization problems, including the nonlinear and high dimensional discrete systems [40,41] (Li et al. 2014;Cao 2014).Li and Gao both adopt SFLA to solve the problem of economic dispatch and unit commitment and put forth targeted improvements to overcome drawback of SFLA, thus enhancing performance in searching and convergence ability.Nonetheless, the local update formula of SFLA still has limitations, leading the search degree of SFLA not enough to find the optimal solution and easy to be premature convergence, especially when all the frogs are concentrated in a local search field.In fact, the research concerning SFLA is still in the initial stage, in both theoretical and practical application especially the STHGS of hydropower station.
The paper focuses on effective improvements and strategies to overcome defects and strengthen performance of SFLA.However, the initial and fundamental SFLA is the real code version with continuous variable value, not appropriate for hydropower systems involving the discrete value (0 and 1) which represents online or offline status of hydropower units, respectively.For that reason, it is essential and necessary to transform real coded SFLA into the binary or discrete version BSFLA.Consequently, the transformational BSFLA is suitable for dealing with the 0/1 programming problem, i.e., the discrete space of unit commitment and scheduling.
Meanwhile, the real coded RSFLA is designed for the ELD problem.The improvements for BSFLA contained population initialization operated by renewed logistic chaotic map, renewed frog subpopulation division (leader, follower, and mutation frog), the cloud model for enhancing evolution and exploitation capability of leader and mutation frog subpopulations, and improved search strategies for follower frog subpopulations.Furthermore, the self-adaptive activation mechanism is introduced to the RSFLA.In summary, the improved discrete-real coded shuffled frog leaping algorithm (IBR-SFLA) is developed to solve the STHGS problem.The improved IB-SFLA is proposed as external search to optimize the UC subproblem; simultaneously the improved real coded SFLA (IR-SFLA) is as the internal search to handle ELD subproblem.Finally, the case study on Three Gorges hydropower station is used to test overall performance of IBR-SFLA.The simulation results show that it can figure out reliable high-quality solutions and outperform contrastive approaches and algorithms in precision and stability, while costing shorter time under the same operation condition.

STHGS Problem Formulation
As mentioned above, the model for STHGS includes two submodules: the UC and ELD.The aim of STGHS is to minimize, to the utmost, water consumption with meeting the power load demand, water balance, spinning reserve constraint, and other system constraints simultaneously in the duration of given operation periods, i.e., an hour for interval in a day.

Economic Load Dispatch (ELD).
As one subproblem of STHGS, the ELD tends to pursue the minimum water consumption of generating power in hydropower station.The objective function can be described as denotes the water consumption during operation periods;   is the th hydro unit;   &   represent load for th unit to burden and water head for th unit, respectively;   (  ,   ) is the water consumption by th hydro unit.
Specifically, the water consumption function is set as quadratic function through curve fitting method.The function is expressed as ,   ,   are coefficients in the fitting quadratic function and other parameters have the same meaning as the former.

Unit Commitment Problem.
Reasonable hydropower units arrangement can bring enormous benefits on power generation efficiency.The basic goal of unit commitment is to effectively schedule the on/off state of hydropower units at the cost of minimum water which is consumed by power generation and switch process of startup and shutdown.The ascertained upstream and downstream water level, load demand, and various constraints are also taken into account [1,30] (Zheng et al. 2013;Chen et al. 2015).Normally, the objective function can be signified as follows.
denotes total water consumption of hydropower units composed of power generation and startup and shutdown costs of units;  up, ,  dn, are startup and shutdown water consumption, respectively.Δ is interval, i.e., one hour;  is the total scheduling periods;    indicates water head for th unit at th interval;    denotes power output for th unit at current th interval;    (  ,    ) is the outflow of the th unit at th interval corresponding to water head   and power output    ;    valuing 1 and 0, respectively, indicates startup and shutdown state of th unit.

Constraints.
As is mentioned above, there are a number of various hydraulic and electrical constraints in STHGS problem.It can be summarized as follows: is the power demand at th interval of hydropower station;    and    have the same meaning as the former.
denotes spinning reserve at th interval.

C Hydropower Unit Output Limits
,min and  ,max signify minimum and maximum load output of th unit, respectively.

D Prohibited Operating Region Constraint
represents total prohibited operation zones, with  to number them.(   )  and (   )  are lower and upper limits corresponding to th prohibited operation zone for the th unit.

E Minimum Up and Down Time Constraint
, and  , are the minimum online (up time) and offline (down time) limits, ensuring stable and safe operation of hydropower units;   i, and   i, denote accumulated up/down time for th unit that had been online/offline till th interval.

G The Water Balance Constraint
is the storage volume of reservoir at th interval;   ,   are inflow and outflow (water release); Δ is interval.

The Stable Operation Region.
As is demonstrated in Figure 1, the operation area of hydropower unit can be classified into three components: prohibition operation, limited operation, and stable operation regions or zones.Constraint D mentioned describes the prohibited operation limits of turbine-unit.Enormous vibration and noise will emerge when hydro unit operates in the prohibited operation zone (POZ) which brings great harm to unit security and stability.Therefore, to ensure long-term security and prolong the service life of units, it is vital to regulate unit operation zone and limit its operation within stable operation zone (SOZ).The dark blue area in Figure 2 shows the unit SOZ of Three Gorges left bank (VGS and ALOSTOM) under different water heads.

The Unit Operation Efficiency.
The unit operation efficiency (UOE), another important turbine-unit characteristic, has a significant impact on generation and total water consumption.Specifically, the relationship between power generation and UOE can be expressed as the following formula: denotes unit operation efficiency (%). is water discharge and  is effective water head. represents the power output (kW).Under a certain water head, different unit power output corresponds to various efficiency which is shown in Figure 3.The curves in Figure 3, displaying UOE in Xiluodu hydropower station, do not present tendency of monotonic increasing but a sharp reduction after the maximum point [42].The variation trend of curves is consistent with Figure 1.By analyzing the UOE curves under different operation zones in Figure 1, we can safely conclude that a single unit will have a relative higher efficiency when working under the stable operation zone (SOZ) than that in other zones.In this paper, the UOE and SOZ are fully considered into the ELD subproblem to fully fulfill purpose of minimizing water consumption.

The Unit Water Consumption
Rate.The definition of hydropower unit water consumption rate (UWCR) is the water consumption corresponding to unit power output.Generally, it can be acquired by the following formula: where   (  ,   ) denotes generation flow. is water consumption rate of th unit when power output is   and net water head is   .From formula (12), we can see the UWCR is concerned with net water head   and  (UOE), of which differences are great among different units.The UWCR is introduced as an influencing factor, aiming to preferentially start up unit with lower UWCR and shut down unit with higher UWCR.Thus, the water consumption can be saved and spinning reserve constraint (constraint B) is satisfied simultaneously.

The Introduction of IBR-SFLA and Model
Establishment for STHGS Problem  incorporates advantage of memetics evolution and particle swarm algorithm which possesses characteristic of less parameters and easiness to comprehension [43] (Zhao and Chen 2016).The frog population is partitioned into a number of parallel subpopulations (communities), of which frogs are arranged to search prey and evolve independently, i.e., memetic evolution.Frogs in each subpopulation can be seen as potential feasible solutions and idea toward a goal and are infected by ideas of other frogs [44] (Zou et al. 2012).Memetic evolution aims to improve quality of individual frog and enhance performance, thus, directing frogs to exploit more excellent ideas and solutions.After a number of memetic evolution steps, the shuffling process, i.e., mingling frogs from different subpopulations, is activated to boost frog quality after exchanging information with each subpopulation so as to pass down the superior ideas [45] (Luo et al. 2015).
Generally, the frog population division strategy can be described as follows: to begin with, the initial population will be ranked in descending order based on fitness of each frog.Then, the frogs are grouped into  subpopulations which satisfy  =  × ℎ; i.e., each subpopulation contains ℎ frogs.Finally, the first frog will be assessed in first subpopulation and the next one is chosen to the next subpopulation, and so forth, until the th frog is grouped in the th subpopulation.The subpopulations are allocated with frogs in the first round, then the (+1)th frog is divided into the first subpopulation, the (+2)th frog is divided into the second subpopulation, and so on [39] (Sun et al. 2016).The rest of subpopulations adopt similar grouping measure to make all frogs classified so as to find the optimal and worst frog   and   within each subpopulation as well as the current global optimal frogs   in the whole population.
The local search in each subpopulation merely renovates the worst frog individual in frog subpopulations.Usually, the position of frog is updated by the following formula: rand( ) denotes random function generating the number in the interval [0, 1];   is leap step of frog,  =1,2,3... ;  min and  max are minimum and maximum frog leap step, respectively;   and   denote the best and worst frog positions of each subpopulation, respectively;   is the new updated worst frog position.
The position of   will be replaced by the new frog position of   if the fitness of renewed frog   is superior to that of previous one   .Otherwise, new frog position is acquired once through (13), in which   is replaced by   .Then, if the process is still incapable of improving performance of worst frog   , a new frog   is obtained randomly to displace the worst frog position   .All frogs are mixed together and rearranged in the frog shuffling process when local search is accomplished in all subpopulations.Then, the local evolution and global shuffling will perform repeatedly until the maximum iterations preset or convergence criteria are satisfied [44].

The Improved Binary Coded SFLA (IB-SFLA).
The UC is a discrete subproblem which involves 0-1 state variables.Therefore, the binary coded SFLA with improvements is adopted to find commitment schedule (unit on/off state).

The Grouping Strategy.
Usually, the frogs are sorted in a descending order according to fitness in the traditional SFLA.Thereafter, the frogs are divided into subpopulations in order.Consequently, relatively worst frogs are sorted into last subpopulation.The population grouping limitation makes the worst frog in last subpopulation have difficulty in acquiring the excellent information and improving its quality.
To avoid this limitation, the frog grouping process is updated.The frog population is divided into three types of memetic subpopulations: leader frog, follower frog, and mutation frog.As the elite individual in whole population, the leader frog subpopulation plays a role in exploring more superior position.Follower frog tends to share information and mutation frog exploits solution scope of unknown area.Frogs in each type of subpopulation experience evolution with different strategies or methods.Therefore, the relative worst frog can better develop new idea with evolution strategy suitable for itself.Meanwhile, the elite information exchange and communication between parallel memetic subpopulations will be more effective.The grouping method is elaborated as follows: Step 1. Select the optimal frogs in all subpopulations and calculate average fitness of these frogs ( V  ).
Step 2. Obtain the optimal frogs whose fitness is superior and inferior to ( V  ) in all subpopulations.Then figure up average fitness of available frogs ( V  ) and ( V  ), respectively.
Step 3. Count number of subpopulations in which fitness of the optimal frog is superior to ( V  ) and inferior to ( V  ), respectively.The number is stored in variables   and   ; specifically,   is the number of leader frog subpopulations and   denotes number of mutation frog subpopulations.
Step 4. Set up total number of subpopulations .The number of leader and mutation frog subpopulations has been acquired in Step 3. Therefore, the number of follower frog subpopulations is  −   −   .

The Modified Search Strategy for Followers Frog Subpopulation.
In the improved binary coded SFLA (IB-SFLA), the frog position  is discretized and coded with values 1 or 0 (on/off), while the frog leap step  need not be discretized.The values of  determine whether  is 1 or 0 in the form of probability.In this paper, the improved sigmoid function is developed to build relationship between frog position  and leap step .The function details are defined in (15) and graph of function is shown in Figure 4.
[, , ] is the leap step of th unit attached to the worst frog in th subpopulation at th period;   , [, , ] is the old location of the worst frog in previous th iteration;   ,+1 [, , ] is location of the renewed worst frog in th subpopulation in the  + 1 th iteration.
Furthermore, in the frog evolution process of traditional SFLA mentioned above, worst frog in subpopulation is updated primarily according to the best frog   in subpopulation.If the new position of the worst frog is still worse than before, frogs are evolved according to the global optimal frog   in (13).Under this evolution strategy, the information is so scarce for worst frog that it is difficult to develop new idea to the fullest extent.Therefore, the old position updating formula is modified to fully pass down the information and novel idea of   and   simultaneously.The new updating strategy is shown in (16) and conducted if the new frog is not to be improved after traditional evolution process.The frog position evolution is still according to (15).
+1 is the renewed frog leap step;   and   are location of the best and worst frogs, respectively;   denotes the best location in whole frog population currently.Both  1 and  2 are rand( ) function which creates random number in the interval of [0, 1].

The Elite Evolution Strategy for Leader
Frog Subpopulation Based on NCM.The cloud theory was put forward by Li et al. in 1995, which possesses the prominent characteristics of randomness and stable tendency, and it has been proved to avoid the local extreme to a certain degree [46].In recent years, the cloud model has been introduced and combined with a number of algorithms such as GA, PSO, and evolution algorithm [47][48][49]  According to the principle of sociology, the more excellent individual often appears around the current optimal individual in subpopulations.In other words, the better frog is found with great probability in the vicinity of the current best frog.Therefore, the normal cloud model is introduced into SFLA to seek the potential best frog within a feasible zone and reduce possibility of trapping in the local optimum.
Normal cloud model (NCM) is brought forth to describe the uncertain conversion relation of qualitative concept or qualitative knowledge with its quantitative expressions.The cloud model is proposed for handling combination of qualitative and quantitative problem.The model reflects the correlation of randomness and fuzziness, establishing the mapping relation with quantitativeness and qualitativeness. is defined as the domain which is one-dimensional or multidimensional. is a qualitative concept corresponding to .For a certain element  ∈ , a random number with stable tendency  = () is defined as certainty of  to .The distribution of  in the domain  is called cloud model and each  is one cloud droplet in the domain space.The expectation   , entropy   , and hyperentropy   are used for characterizing the digital features of cloud [50] (Ma et al. 2013).The details on   ,   , and   are as follows: A  is central of all droplets and the most representative cloud droplet for qualitative concept.The optimal frog in leader frog subpopulations is selected as the representative cloud droplet.
B The entropy   denotes the uncertainty measurement of qualitative concept.It is not merely reflecting emergence randomness of cloud droplet but also revealing the correlation of fuzziness and randomness.The scope of generating cloud droplet and its randomness highly depends on the value of   ; specifically the larger   is, the more obvious the randomness is and the more broader the scope is.  is valued with the standard deviation of the frog fitness in subpopulation during current iteration.
C The hyperentropy   determined by randomness and fuzziness represents measurement of entropy, i.e., the entropy of entropy.To boost search randomness in initial iteration stage and maintain search stability in the later stage,   is relevant to   with value of   /5.
The process of creating cloud droplet according to basic normal cloud generator is detailed as follows [39,51]  Step 1. Generate a normal random number   with expectation   and standard deviation   .
Step 2. Thereafter, generate a normal random number  as a cloud droplet, which is taken   , |  | as expectation and standard deviation, respectively.
Step 3. The calculation   and  are plugged into formula to figure up certainty pertaining to the qualitative concept .
Step 4. Repeat the above-mentioned processes until the cloud droplets satisfy the terminal condition.
The optimal individual in leader frog subpopulation is regarded as a normal droplet  (  ,   ,   ).Then, several cloud droplets are created according to this droplet by normal cloud generator.The number of new created cloud droplets is the same with capacity of subpopulation ℎ.The new droplet replaces original optimal frog if fitness is better than before.Otherwise, the original optimal frog is still to be retained [52] (Zhang and Li 2015).The elite evolution strategy aims to evolve the best individual in subpopulation, while the worst frog adopts the traditional updating mechanism to finish evolution.

The Mutation Strategy for Mutation Frog Subpopulation.
There are variation phenomena in the biosphere, which has characteristic of randomness and burstiness.However, the biomutation is beneficial to track the potential optimal individual in the unknown solution space and avoid falling into local optimum to some extent.
Through analysis of the grouping strategy, it is known that relative worst frogs are sorted into the mutation subpopulations.In comparison with other frog subpopulations, the mutation frog is weaker than others to make contribution to developing novel idea.Therefore, the mutation frogs are designed to find profitable variation and corresponding update strategy for mutation frogs is adaptively changed using chaotic mapping method.
The logistic map is a typical discrete-time dynamical systems, i.e., nonlinear equation.The chaotic sequence is produced by logistic map, which is seemingly chaotic but with regularity and ergodic property.In this paper, the modification is developed for conventional logistic map to strengthen population diversity and distribution uniformity.The chaotic sequence is calculated by the following reformative logistic map equation [53] (Zou et al. 2016): represents decision variable dimension;  is the current iteration;  max is maximum iterations; variable    is defined as the th chaos variable after th iteration, valuing in the interval [-1, 1].
Each dimension of chaotic variables   is reflected to feasible solution space of optimization variable, i.e., frog step   by (18).Then, the conversion relation between mutation frog leap step and position  is established through (15).
If the new position   of mutation frog obtains a better performance than the overall optimal frog, the position of current optimal frog is replaced by this position.Otherwise, calculations in (17) and (18) are repeated until the given terminal condition is satisfied.The new individual created by mutation strategy based on chaotic mapping possesses characteristic of randomness, fuzziness, and determinacy.This mutation strategy contributes to effectively exploit space beyond local convergence domain, thus guiding algorithm to locate the more outstanding frog.

The Improved Real Coded SFLA (IR-SFLA).
The ELD subproblem is a continuous, highly nonlinear, and multidimensional optimization problem, aiming at efficiently dispatching power load of hydropower units.An optimal dispatch schedule is beneficial to minimize water consumption.Because the power load is continuous, the real coded mechanism of SFLA is suitable for handling ELD subproblem.However, the evolution strategy in traditional SFLA is liable to slow down convergence speed and not easily escape from local optimum.Thus, several changes for the original real coded version of SFLA (IR-SFLA) are introduced in this paper.

Initial Population Based on Chaos
Theory.The initialization has a significant effect on convergence speed and final solution quality.Traditional population initialization is accomplished through the random guided search strategy, making it difficult to locate feasible zone.Chaos theory possesses inherent stochastic characteristic and ergodicity, which is highly sensitive to initial value and liable to escape from local optimum.Therefore, the improved chaos logistic map is applied to frog population initialization to strengthen population diversity and distribution uniformity during the initialization process.It is beneficial to promote overall quality of initial population and locate feasible search zone effectively [53,54]  Similarly, the chaotic sequence   creation process is according to (17), and chaotic variables are mapped to feasible solution space of variable   on the basis of (18).In IR-SFLA, variable   is real coded, which denotes the power load in ELD subproblem and the frog position updating formula is in line with (14).Then, calculate fitness (  ) of   which is concerned with water consumption.If (  )  ≤ (  ) +1 , i.e., the new position is superior to the old created in the last generation, then (  )  = (  ) +1 and    =   +1 .The process is repeated in a number of iterations and select the optimal frogs into initial population.

The Improved Search and Position Evolution Strategy.
As mentioned in IB-SFLA, the information is too scarce to guide the update of worst frog with using search strategy of traditional SFLA.Similarly, the information of local optimal frog   and overall optimal frog   is combined to lead the worst frog into novel position.Furthermore, drawing on the PSO, the history frog leap step is taken into consideration to further enrich update information of worst frog.The improved search strategy is shown in the following equation, and the position updating is according to (14).
,   ,   are real coded with continuous variable, i.e., unit power load in ELD.Similarly,   denotes historical frog leap step;  1 and  2 are evenly distributed random number between (0, 1);  is adjustment coefficient, which aims at controlling leap step and convergence speed in order to achieve balance between local deep-searching and global information shuffling process.In this paper, the improved adjustment coefficient equation is shown in is the number of global shuffling iterations;   is total global shuffling iterations;  is the number of local iterations in memetic subpopulations;   represents total local iterations in subpopulations;  max ,  min are the maximum and minimum adjustment coefficient;  is weight coefficient in the interval of (0, 1).
At the initial evolution stage, the adjustment coefficient  is assigned with a larger value to boost search capacity in global solution space and a lower value at the later evolution stage.A lower adjustment coefficient  slows down the evolution speed and promotes depth search in local feasible solution space in the later stage [55].

The Adaptive Frog Activation Mechanism.
The local optimal individual   , overall optimal   , and worst frog   to be updated are very close with great possibility if frog population search and evolve in a narrow space, i.e., local optimum.Assuming that the difference of (  -  ) is close to zero, population diversity is reduced and new frog position update for   will weaken and even stagnate.Therefore, the adaptive frog activation mechanism is introduced into frog update and position evolution process.The activation mechanism replaces inactive frog lack of evolution ability with creating a new individual.Consequently, search ability of   is reactivated and new position tends to be found.The new frog creation method is expounded as follows.
is random number in the range of [0, 1];  denotes adjustment coefficient which uses the arc tangent function to make it changeable dynamically;   and   are initial and final value of  which are evaluated with 0.9 and 0.4, respectively;  represents controlling factor in the range of [0.4,0.7]; ,   ,   , and   are defined in (20).To maintain fast convergence characteristic of IR-SFLA and reduce uncertainty of random number, a relative lower  is used in the later iteration stage.|(  ) − (  )| ≤  is decision condition;  is a positive number defined as similarity threshold; (  ) and (  ) are fitness of   and   , respectively.Equation ( 21) is activated if difference of (  ) and (  ) satisfies the decision condition in a certain number of iterations.

The Numerical Simulation
Experiment.Numerical simulation experiment is designed to demonstrate the superiority of IR-SFLA in solving continuous space optimization problem.The typical Sphere, Rosenbrock, Rastrigin, Griewank, and Schaffer benchmark problems with 30 dimensions are selected as simulation experiment.The IR-SFLA is compared with PPNGA, EACO, SFLA, CGSA, and IBCO, and 20 independent simulations are conducted for all algorithms.Partial results of the 20 simulations are shown in Figure 5.The Mean and Standard Deviation (STD) acquired from 20 independent simulations are listed in     speed and reliability of each algorithm, respectively.The results are also listed in Table 1 where the "NA" denotes "Not Applicable" and the convergence is set as 10 −8 for all the benchmarks.To further compare IR-SFLA with other rivals, the Friedman and Kruskal-Wallis tests are used to depict the algorithm superiority and listed in Table 2.
Through inspection of the iteration curves in Figure 5 and results of Mean and STD in Table 1, it is clearly revealed that IR-SFLA eventually converges to the evidently beat solution for each benchmark problem in comparison with other algorithms.Furthermore, the performance on convergence speed for IR-SFLA is much better than competitors especially in terms of F 2 Rosenbrock, F 3 Rastrigin, and F 4 Griewank.For F 1 Sphere, F 3 Rastrigin, and F 4 Griewank, the performance for IBCO is outstanding and converges to the best solution; however, the convergence speed for IBCO has no advantage compared with IR-SFLA through the iteration curves in Figure 5.Moreover, the IBCO is not skilled in solving F 2 and F 5 benchmarks from the aspects of final precision and stability reflected by STD in comparison with IR-SFLA and PPNGA.The method PPNGA presents instable performances which are reflected in relative good performance on F 2 and F 5 , especially the F 4 benchmarks and poor performance on F 3 .
From the MNIG and SR in the bottom of Table 1, the SR of IR-SFLA are all "100%" which denotes IR-SFLA can obtain the final theoretical optimum of each benchmark in 20 independent simulations which verify the high reliability standard of IR-SFLA.Furthermore, the MNIG of IR-SFLA are the minimum within all the algorithms.However, from the observation of compared algorithms, although the IBCO can get the best results in F 1 , F 3 , and F 4 , its MNIG are much bigger than that in IR-SFLA that indicates the lower convergence speed than the latter.The PPNGA performs well in F 4 but is difficult to converge to optimum of other benchmarks.The rest of compared algorithms can not acquire the theoretical optimal solutions of each benchmark but perform better if the precision demand is relatively low.Therefore, the IR-SFLA displays the remarkable advantages on reliability and convergence speed compared to others in terms of F 1 and F 5 benchmarks.
In terms of test ranks from Table 2, the less test result is, the better performance algorithm has.It can be seen that IR-SFLA is well positioned among all algorithms which denotes better performance and significant advantages on computational precision and reliability in comparison with other algorithms.
From the results summarized above, it can be concluded that the improvements for SFLA present efficient effects on overall search capability and keep a fast convergence speed compared with other rivals.To sum up, the IR-SFLA is more adaptive and effective in solving the optimization problem in continuous space.

Structure and Initialization of Frog Individual for UC.
The frog individual of UC subproblem represents on or off (1 or 0) state of the th unit at th interval.Therefore, the structure of individual for UC can be described as an integermatrix  , with value 0 or 1 and shown as follows: where  and  are the number of units and operation duration, respectively; the upper limit of  is 24, i.e., a day;  , represents the state of th unit at th interval. , is initialized with 0 or 1 at each period through following (23).
where  , is the same meaning with that in (22); rand( ) denotes random number between 0 and 1.

Strategies for Handling the Constraints (1) The Supplement and Repair Strategy for Unit Reserve
Constraint.There are two cases during the process of optimal load distribution between each unit.Case 1 is that the total maximum load between all operating units ∑  =1    ,max is less than the sum of requirement  , and reserve capacity    ; Case 2 is the total minimum load of all operating units; ∑  =1    ,min is more than load requirement  , and reserve capacity    .In terms of the Cases mentioned above, a dynamic processing method based on priority list of unit startup and shutdown is proposed.The priority list is ranked in ascending order according to the average unit water consumption rate (AUWCR).The AUWCR calculation formula is shown in V denotes the average rate of water consumption;   is the distributed power load within stable operation region (SOR);   (  ) represents the corresponding water consumption;  is maximum power load within SOR, setting 1MW as an interval.
If Case 1 is activated, turn up the current offline units satisfying the minimum down time constraint in turn according to ascending order in the priority list until the reserve capacity constraint is satisfied.Otherwise, shut down the online units satisfying the minimum operation time constraint in turn according to descending order in the priority list until the reserve capacity constraint is satisfied.
(2) The Strategy for Repairing Minimum Up and Down Time Constraint.To satisfy the minimum up and down time constraint, a repair strategy is conducted to adjust the units on and off state.The variables (, , ) and (, , ) are introduced to record startup and shutdown duration of each unit attached to th frog at th period.The record process is described in (25).Meanwhile, the startup and shutdown durations stored in (, , ) and (, , ) are verified whether the minimum up and down time constraint is satisfied or not.If not, unit state remains unchanged; otherwise, the frog can realize the state transition in line with (15).

The IR-SFLA to Solve ELD Subproblem Based on Optimal Economic Operation Table.
A wise and feasible unit commitment is acquired with handling the UC subproblem, and proposed IR-SFLA is applied to realize highly efficient load distribution between committed units.

Structure and Initialization of ELD.
The structure of ELD is correlated with that in UC, which consists of power output for all committed hydropower units.The specific structure is illustrated as follows: where  , represents the power output for the th unit at th interval; the meaning of  , is in line with that in UC structure.Different with UC initialization, the ELD is real-coded continuity optimization problem.Therefore, initialization process of ELD is according to method based on chaotic theory.The logistic mapping equation is in line with (17); however, the reflection equation is demonstrated as follows with corresponding changes according to (18).
denotes the chaotic variable generated by improved logistic mapping formula;  ,min and  ,max are lower and upper limit for th unit. .The economic operation for large scale hydropower station, especially in terms of Three Gorges hydropower station, is multidimensional, complex nonlinear optimization problem.In order to simplify the solving process and significantly improve the computation efficiency, the optimal economic operation table is established by the proposed IR-SFLA.The table formulation is on premise that all power units are involved in power generation; i.e., all units are startup.In addition, the SOR of units is taken into consideration to keep unit operated in high efficiency and safe region.The established economic operation table is stored in database, thus making it easier to call the table.The solving formula is demonstrated as follows:

The Optimal Economic Operation Table
1 ,  1 ( 1 , ) = 0 are boundary conditions in the initial stage;   is accumulative load of units;   (  , ) denotes water consumption of power generation when total load and water head are   and , respectively; (  (  , )) is the best frog individual.

The Process of Invoking Optimal Economic Operation Table.
As mentioned above, the optimal economic operation table is stored in database so as to be conveniently invoked.The number of startup units is recorded using the variable  when a unit commitment is obtained by the UC submodule.The variable  is defined as a count variable to locate the th line in optimal economic operation table.The load distribution between units is obtained by locating the specified line in optimal economic operation table according to the load demand  , corresponding to each interval.The linear interpolation method is adopted to acquire the load distribution if it cannot be located directly in the table.The flowchart of invoking optimal economic operation table is demonstrated in Figure 6 and the detailed procedures are shown as follows.
Step 1. Seek the first line (=1) of economic operation table and figure up the total load  1 of  operational units according to (29).Compare  1 with load demand and judge whether it is submitted to  1 =  , or not.If not,  =  + 1 and turn to Step 2. Otherwise, terminate the search and ( 30) is adopted to obtain the load distribution of operational units and corresponding water consumption.
*  () denotes load distribution of th operational unit in th line of optimal economic operation table;  1 is total load distribution of all the operational units in th line.  is intermediate variable to store  *  ().
Step 2. Calculate  1 and  2 in line with (31). 1 ,  2 are the total load distribution of operational units Step 3. Adjust  according to  =  + 1 and turn to Step 2 if the criterion  2 <  , is satisfied.Otherwise, continue to judge whether the criterion  2 =  , is satisfied; if it is, the power load of each operational unit is obtained by   =   * (),  = 1, 2, 3.....Otherwise, the linear interpolation of ( 32) is applied to obtain load distribution of operational units.
Step 4. Calculate water consumed by power generation and startup and shutdown transition between hydropower units based on NHQ curve.

The Coupled External and Internal Model for STHGS
Problem.The STHGS problem consists of the unit commitment (UC) and load distribution subproblems.Therefore, the coupling model consists of the UC (external submodel) and ELD (internal submodel).The UC submodel delivers message to the ELD submodel in order to realize optimal load distribution and compute total water consumption.On the contrary, the information about water consumption due to power generation and unit startup and shutdown pass back to the UC submodel.The specific steps are listed as follows.
Step 1. Initialize the parameters of IBR-SFLA, including number of frog populations , subpopulations , range of frog leap step, inertia weights  st and  en , threshold of heuristic frog activation mechanism, iterations in subpopulations, and whole populations   and   .
Step 2. Generate initial frog population by (23) where  , is initialized with 0 (shutdown) or 1 (startup) at each period using the rand( ) function.
Step 3. Adopt the supplement and repair strategy of unit reserve capacity to satisfy the sinning reserve constraints.Record the duration of unit startup and shutdown according to (25).Regulate the unit state at each period with the repair strategy of minimum startup and shutdown time in order to satisfy constraint E.
Step 4. Invoke and seek the optimal economic operation table to obtain load distribution of operational units and total water consumption including electricity generation and startup and shutdown of units.Calculate fitness of frogs which are ranked in descending order.
Step 5.The frog grouping strategy is applied to divide population into leader frogs, followers frogs, and mutation frogs subpopulations.
Step 6.The leader frogs subpopulation adopts the elite evolution strategy to search and update optimal frog in subpopulation while the worst frog renewal in subpopulations is conducted through traditional frog local search strategy; mutation frogs are designed to exploit potential global optimal frog in feasible region on the basis of mutation strategy; the followers frogs update is based on modified search strategy.The process of handling constraints is still according to strategies in Step 3.Then, finish the load distribution through the optimal economic operation table mentioned in Step 4. All subpopulations belonging to each type of frog are updated with corresponding strategy mentioned and renew the local best frog   and global best   in each iteration.
Judge whether the iteration of local search is up to the number preset or not; if it is, turn to Step 7. Otherwise, go back to Step 6 and continue to search until the local update is finished in all subpopulations.
Step 7. All subpopulations will be mixed to exchange information and idea between different subpopulations.Record the best frog, minimum water consumption, and optimal unit commitment and load distribution in current global shuffling iteration.
Judge whether the global shuffling iteration is achieved or not; if not, turn to Step 5.If it is, the global best frog, unit commitment, and load distribution are found and output the final results.
The flowchart of coupling external and internal model for STHGS is demonstrated in Figure 7. MW and 6 hydro units with 700MW installed capacity are equipped afterwards.Therefore, the total installed capacity of hydropower station is up to 22,500MW.In this paper, the initial 26 hydro units containing 6 left-bank VGS units, 8 leftbank ALSTOM units, 4 right-bank ORIENT power units, 4 right-bank ALSTOM units, and 4 right-bank HEC units are selected as simulation operation objects.In this study, the 75m, 88m, and 107m water head are selected as three simulation cases and scheduling periods are 24 hours.The stable operation regions (SOZ) of generating units under 75m, 88m, and 107m water head are listed in Table 3.Meanwhile, the flow characteristic curves (NHQ curves) of hydro units corresponding to 75m, 88m, and 107m water head are demonstrated in Figure 8.The daily load demands at each interval for hydropower station are drawn in Figure 9.

Parameter Settings.
The parameter settings significantly affect the efficiency of the IB-RSFLA including search ability, stability, and computation speed.In this study, we set a group of parameter combinations and make comparison of simulation results to find the best control parameter combination.The optimal control parameter settings are described as follows: The number of frog populations  is 80; the number of subpopulations  is 10; the number of frogs in each subpopulation ℎ is 8; the frog leap step is valued in the interval [-4, 4]; the value of inertia weights   and   are, respectively, 0.9 and 0.4; the iterations number of local search in subpopulation   is 10 and global search   is 50.In addition, the similarity threshold  used to motivate the heuristic frog activation mechanism is tested to value 10 −3 .

Results and Analysis.
The 75m, 88m, and 107m water heads are selected as the cases and 20 independent simulations are conducted to find the best solution.The unit commitment and economic load dispatch corresponding to the best solution under 75m, 88m, and 107m water head  are demonstrated in Table 4, respectively.In inspection of economic load dispatch in Table 4, all hydropower units are operated under stable operation region at each interval, and the hydropower unit output limits C and prohibited operating region constraint D are rigidly satisfied.Furthermore, Table 4 shows the corresponding unit commitment (UC) strictly satisfies minimum up and down time constraint, which illustrate the state of units (startup and shutdown) is reasonably arranged using "strategies for handling the constraints" in Section 3.4.2.The effective UC contributes to avoiding units startup and shutdown frequently and prolonging the service life of hydropower units.
The iterative curves of best and average water consumption at each generation under 75, 88, and 107m water head are shown in Figure 10.Furthermore, a bit of "fluctuant points" in the curves denote the possible local optimal point.Before "fluctuant points", the search by IBR-SFLA is close to standstill state which indicates algorithm trapped in local optimum.However, the search is reactivated and explored further to gradually converge to best solution after "fluctuant points".The reactivation can be attributed to evolution operator, modified search strategy of followers frog proposed in IBR-SFLA which contribute algorithm to effectively escape from local optimum and motivate overall search capability.The number of periods (h)

Load demand
Figure 11: The daily load demand and real power output curve of hydropower station in the Three Gorges.
Table 5 gives details of total water consumption at each period corresponding to 75m, 88m, and 107m water head in which the water consumed by power generation and startup and shutdown of units is included.
To verify whether the power balance constraint A is satisfied, the comparison between daily load demand and real power output is presented in Figure 11 which reflects that real power output satisfies the load demand at each interval.
In order to profoundly demonstrate the superiority of IBR-SFLA in solving STHGS problem, a comparison experiment between IBR-SFLA and limited adaptive genetic algorithm (LAGA), extended ant colony optimization (EACO), chaotic gravitational search algorithm (CGSA), firefly algorithm (FA), and improved bee colony optimization (IBCO) is carried out under the same independent simulation times, iteration number, and population scale.The final optimal simulation results obtained by IBR-SFLA and other compared algorithms are listed in Figures 12 and 13. Figure 12 presents iterative process of total water consumption for each algorithm.From Figure 12, it can be observed the LAGA is less effective than other algorithms.The uncertainty of control parameters in LAGA dramatically affects the quality of solutions while the method IBR-SFLA proposed has relative less parameters we need to control.Results obtained by IBR-SFLA are obviously more superior than others and water consumption descends drastically at the initial intervals and converges to stable optimal level at the later phase.The IBR-SFLA outperforms the other algorithms through inspection of the convergence precision, which indicates the efficiency and adaption of IBR-SFLA in dealing with STHGS problem.
More details will be revealed if the results obtained by each algorithm in 20 independent simulations are expressed with form of box plot (Figure 13).The box plots present distribution of the simulation results from which IBR-SFLA exhibits significantly lower and narrower box plot with short whiskers than that in other rivals under 75m and 88m water heads.It can be stated that solutions obtained by IBR-SFLA in 20 simulations are higher quality and closer to the optimum than compared algorithms.In terms of 107 water heads, the  boxes obtained using IBR-SFLA are lower than others but are not well at the coverage area of box plot.However, the whisker of box plot for IBR-SFLA is short which indicates that IBR-SFLA outperforms other algorithms in convergence precision and reliability.To sum up, the higher quality results obtained by IBR-SFLA in 20 simulations demonstrate more effective and stable performance in comparison with other methods.
The results analysis simulated by each algorithm is listed in Table 6, in which the water consumption, standard deviation, and mean simulation time are shown simultaneously.The advantages of the IBR-SFLA compared with rivals are presented in Tables 6 and 7 and Figure 14.In terms of convergence precision, it is clear that IBR-SFLA can seek the less water consumption through more reasonable unit commitment and load dispatch.The results comparison demonstrates that the improved strategies introduced to IBR-SFLA clearly improve the search capability and overall efficiency to converge to better quality solution.
Moreover, the simulation time cost by IBR-SFLA is observably less than other compared methods that reflect higher efficiency in dealing with STHGS for large scale hydropower station.The relative higher search speed has close relationship with optimizing mechanism of IBR-SFLA.Only the worst frog in frog subpopulations is updated during each iteration; therefore, only one calculation is needed to finish the iteration in the best case.Even if in the worst case, three calculations are required during each iteration.
In terms of algorithm stability, it can be concluded that IBR-SFLA has significant advantage over standard deviation (STD) compared with other algorithms.The smaller STD manifests that computational stability of IBR-SFLA is obviously improved and reflected by narrower box and shorter whiskers in Figure 13.Consequently, the proposed IBR-SFLA is more likely to converge to final optimal solution in each simulation than other rivals.
In summary, IBR-SFLA shows higher solution quality, convergence speed, and reliability and stability in solving STHGS problem with complex hydraulic constraints.Therefore, the proposed IBR-SFLA is verified as a feasible and effective method to save water resources and further enhance waterpower utilization efficiency.

Conclusions
In this paper, an improved binary-real coded SFLA (IBR-SFLA) is proposed for solving the STHGS problem for large scale hydropower station.To demonstrate the superiority of IBR-SFLA, numerical simulation experiment is conducted and results show IR-SFLA has exceptional performance on search capability and convergence speed.Then, IBR-SFLA is applied to the STHGS problem in Three Gorges hydropower station and heuristic repairing strategies are introduced to handle complex constraints such as minimum startup and shutdown, spinning reserve, and water balance.Simulation results indicate IBR-SFLA can get better water consumption solution with costing less computation time in comparison with compared methods.By further analysis, distribution of solutions obtained by IBR-SFLA is concentrated with less scatter, which verifies that stability and reliability of IBR-SFLA hold superiority compared to other methods.In conclusion, the proposed IBR-SFLA enhances overall Mathematical Problems in Engineering performance of traditional SFLA with more reasonable grouping strategy, effective search strategy based on cloud model, and chaotic theory.It is demonstrated that IBR-SFLA provides more effective and feasible way to solve STHGS problem with less water consumption and higher utilization rate of hydropower resources.In the future, related studies are focused on further improving convergence speed and enhancing comprehensive profitability for hydropower.

Figure 1 :
Figure 1: The operation region of left-bank ALSTOM units in Three Gorges hydropower station.

Figure 2 :
Figure 2: The SOZ of left-bank units in Three Gorges hydropower station.

Figure 3 :
Figure 3: The unit operation efficiency (UOE) of right-bank unit in Xiluodu hydropower station.

Figure 4 :
Figure 4: The graph of improved sigmoid function.
(Dai et al. 2007;Zhang et al. 2008;Zhang et al. 2012).The introduction of cloud model algorithm can effectively improve overall search algorithm performance and contribute to overcome the defects on easily falling into premature trap to some extent.

FFigure 5 :
Figure5: The iterative curve of IR-SFLA and compared algorithms for test benchmark functions.Notes: the algorithms convergent to the optimal result "0" after stagnation point.
are marked in bold.

Figure 6 :
Figure 6: The flowchart of invoking optimal economic operation table.

Figure 7 :
Figure 7: The flowchart of coupled external and internal model based on IBR-SFLA for solving STHGS.

The
NHQ curves for 75 m water head Left-bank VGS-type Left-bank ALSTOM-type Right-bank Orient-type Right-bank ALSTOM-type Right-bank HEC-The NHQ curves for 88m water head Left-bank VGS-type Left-bank ALSTOM-type Right-bank Orient-type Right-bank ALSTOM-type Right-bank HEC-

Figure 8 :
Figure 8: The NHQ curve of left-and right-bank hydro units corresponding to 75m, 88m, and 107m water head.

Figure 10 :
Figure 10: The iterative curves of water consumption at each generation corresponding to 75m, 88m, and 107m water head.

5 :
The details of water consumption at each interval corresponding to 75m, 88m, and 107m water head.Time period (h)

Figure 12 :
Figure12: The iterative process of water consumption obtained by IBR-SFLA and compared algorithms corresponding to 75m, 88m, and 107m water head.

Figure 13 :
Figure 13: The box plots of water consumption in 20 simulations obtained by IBR-SFLA and compared algorithms.

Figure 14 :
Figure 14:  The contrasting result between four algorithms corresponding to three types of water head (75, 88, and 107m).

Table
1, respectively.In addition, the Success Rate and Mean Number of Iteration Generation (MNIG) required to achieve convergence precision preset are introduced to compare convergence

Table 1 :
The optimization solutions obtained from IR-SFLA and compared algorithms.

Table 2 :
The ranks obtained by nonparametric test of Friedman and Kruskal-Wallis for IR-SFLA and compared algorithms.

Table 3 :
Stable operation region of each unit in Three Gorges corresponding to 75m, 88m, and 107m water head (Unit: MW).

Table 4 :
The best computational solutions of economic load distribution (ELD) and unit commitment (UC) corresponding to 75m, 88m, and 107m water head.

Table 6 :
The water consumption and average computational time of each algorithm corresponding to 75, 88, and 107m water head.