## Abstract

To efficiently develop power generation and solve downstream ecological health protection in Qingjiang basin, multi-objective ecological operation for cascade reservoirs (MOEOCR) model is established in contrast to conventional models that set ecological water requirement as constraint. The basic, suitable and ideal ecological water requirements in Geheyan and Gaobazhou sections are calculated using a requirement level index. Instead of the traditional evolution mode based on population, we introduce a shuffled frog leaping algorithm (SFLA) which evolves independently in sub-populations. Moreover, the SFLA is converted into a modified multi-objective algorithm (MMOSFLA) with strategies including chaotic population initialization, renewed frog grouping method and local search method, and elite frog set evolution based on cloud model. The water level corridor is used to help effectively handle complex constraints. The IGD and GD indexes are used to evaluate quality of solutions acquired by each method. In terms of normal year, the mean IGD and GD of MMOSFLA are 1.2 × 10^{–1} and 2.75 × 10^{–2}, respectively**.** The scheduling results verify efficient search ability and convergence performance in solution diversity and distribution in comparison with other methods. Therefore, MMOSFLA is verified to provide an effective way to fulfill hydropower and ecological benefits facing the MOEOCR problem.

## INTRODUCTION

Reservoir operation and management have played very considerable roles in flood control, electricity generation, water supply and navigation (Wu *et al.* 2012). Conventional reservoir operation predominantly focuses on maximizing social-economic benefits and addresses ecosystem values only as the constraints on reservoir release and elevations (Jager & Smith 2010; Labadie 2014). With the increasing ecological protection awareness and substantial deterioration in the ecological environment, there is now widespread recognition that water resources management must take account of the vulnerability of ecosystems. Accordingly, the multi-objective reservoir scheduling considering ecological factors has attracted much attention (Mao *et al.* 2016). Vogel *et al.* (2007) designed an ecological flow regime to capture the natural flow variability for maintaining the functional integrity of aquatic ecosystems. In China, increasing attention has also been focused on economical benefit and ecosystem health (Cai *et al.* 2013; Hu *et al.* 2014). In traditional reservoir operation, ecological requirements are typically addressed by meeting certain constraints of minimum instream flow (Castelletti *et al.* 2008). Later, a multi-objective ecological operation for reservoirs with the objective of ecological water requirements has been developed to coordinate conflicting objectives. Currently, the multi-objective ecological operation for cascade reservoirs (MOEOCR) is still challenging as the competitive objectives are difficult to reconcile and the problem is often non-convex, non-linear, multi-dimensional and discontinuous, making it harder to optimize.

Historically, multi-objective problems were solved with conventional methods such as linear programming (LP), nonlinear programming (NLP), dynamic programming (DP) and progress optimality algorithm (POA), through converting the multi-objective problem into a single objective (Feng *et al.* 2017). For ecological operation problems, the ecological factors and objective could be converted into constraints of ecological water demand for aquatic species or arranged weight to describe their importance. Li (2016) gives corresponding weights to target of power generation, flood control and ecology, and adopts a genetic algorithm to optimize the multi-objective Longyangxia Reservoir operation. Kang *et al.* (2010) established an ecological scheduling model of Danjiangkou Reservoir and used the POA method to solve it. Yin *et al.* (2013) and Long & Mei (2017) took cascade reservoirs located in downstream Jinsha River as a case study, and analyzed the reservoir scheduling scheme under different ecological flow constraints. Yang *et al.* (2015) built the cascade reservoirs scheduling model that is coupled with modified DP. Optimization is achieved based on analyzing the relationship between power generation and ecological flow. Although massive computation can be avoided in most cases, there are such issues as sensitivity to transformation coefficients, and difficulty to obtain trade-off solutions (Zhang *et al.* 2013).

Due to the presence of the multiple objectives, it is difficult to find a single optimal solution that optimizes all goals with a conventional single-objective method. Instead, multiple solutions exist in the form of trade-offs, also referred to as Pareto optimal solutions. After identifying the trade-offs, the Pareto optimal solutions are ranked to enable decision makers to select the best scheduling solution (Zhu *et al.* 2017). In recent years, intelligent algorithms and multi-objective optimization techniques have been applied to the hydrology area, including reservoir operation (He *et al.* 2014; Luo *et al.* 2015a, 2015b; Jia *et al.* 2016; Moazenzadeh *et al.* 2018). Both Lu *et al.* (2011) and Wang *et al.* (2013) optimized the target with the minimum ecological water shortage and maximum power generation, and the multi-objective differential evolution method was adopted to solve the problem. Zhang *et al.* (2016) proposed the ecological scheduling model based on the ecological flow interval, and optimized multiple objectives including power generation, corresponding guarantee and ecological guarantee rate by NSGAII. Reddy & Kumar (2006, 2007) employed the multi-objective genetic algorithm (MOGA) and particle swarm to generate a pareto optimal set for a multi-objective reservoir system that serves multiple purposes including irrigation, power generation, and downstream water quality requirements in India. Li & Qiu (2016) put forward a multi-objective reservoir optimization model incorporating ecological adaption and applied this model (using NSGAII) to the Three Gorges Reservoir, the operation of which has damaged the downstream ecosystem. However, to our knowledge, few papers are concerned with the multi-objective reservoir ecological operation problem using the meta-heuristic method. Moreover, the studies mainly focus on the single reservoir ecological operation. In view of complex constraints and targets for large cascade reservoirs compared with a single reservoir, it is more difficult to find an optimal trade-off ecological scheduling solution. The algorithms above, including the NSGA, MOGA, MOEA and MODE, accomplish the evolution and update process based on advantages of the entire population (this means that individuals are updated together in the population community). The evolution information obtained by each individual is from the best individual in the population. The Shuffled Frog Leaping Algorithm (SFLA), by contrast, divides the population into several communities. Individuals in each community are permitted to develop independently. The new population will form through the shuffling strategy after a defined number of iterations. This evolution pattern helps improve the solution by sharing information independently obtained by each community rather than the single information from the entire community.

The SFLA was first proposed by Eusuff & Lansey (2003). It is a memetic meta-heuristic method that is developed to seek global optimal solution by performing an informed heuristic search with a heuristic function (Eusuff & Lansey 2003). The SFLA is a combination of deterministic and stochastic methods which was initially applied into the optimization of a water network distribution design. The deterministic method allows algorithms to effectively use response information to guide a heuristic search, while the stochastic element in SFLA ensures flexibility and robustness of the search patterns. It is found that SFLA can perform the search towards a global optimum using global information exchange and internal communication mechanisms (Sun *et al.* 2016). The SFLA, which is easy and convenient to code with less control parameters, is testified compatible with handling comprehensive optimization problems including the non-linear and high dimensional discrete systems (Cao 2014; Li & Yin 2014).

Qingjiang River is the largest tributary of the middle Yangtze River. Following increasing concerns for ecology health in the downstream watershed dramatically influenced by reservoirs, the ecology protection objective should be considered in reservoir operation. For this purpose, the basic, suitable and ideal ecological water requirements downstream of Geheyan and Gaobazhou reservoirs are analyzed based on the requirement level index. Moreover, to coordinate the hydropower and ecological protection benefits, the power generation in the hydropower station and ecological water spill and shortage are selected as competing objectives in the MOEOCR model. Then, the modified multi-objective shuffled frog leaping algorithm (MMOSFLA) is proposed to achieve multiple trade-off solutions (Pareto optimal solutions) for the model. MMOSFLA is a modified version of conventional SFLA which is initially designed to solve the single-objective optimization problem. We accomplish the conversion and performance improvement with strategies such as population initialization update, renewed frog sub-population partition, dynamic maintenance and evolution mechanism for external elite frog set, and local search capability enhancement. Furthermore, the water level corridor combined with a penalty function is used to handle constraints efficiently. Finally, by using the MMOSFLA, we study the optimal trade-off between power generation and ecological water spill and shortage. Scheduling solutions under scenarios of different ecological flow requirements obtained by MMOSFLA is analyzed to demonstrate the practicability and effectiveness of the proposed model and method.

## THE MULTI-OBJECTIVE ECOLOGICAL OPERATION MODEL FOR CASCADE RESERVOIRS (MOEOCR)

The multiple objectives contain the maximum power generation for the hydropower station located in the reservoir and the minimum ecological water deficit and spill. Both ecological water spill and shortage can lead to adverse effects on the health and stability of ecology. The two objectives dominate each other and we may find it difficult to acquire the conventional optimal solution. Instead, we can obtain a relative optimal solution set. The objectives are expressed by Equations (1) and (2), respectively.

### Objective function

- 1.
- 2.Minimum ecological water spill and shortage: where and denote power generation and ecological water spill and shortage; is the scheduling interval; and are the total amount of reservoir and scheduling intervals, respectively;
*C*represents the quantity of control cross-section which is matched with the reservoirs; is the power output of the th hydropower station at the th interval; denotes the efficiency coefficient of unit in the th hydropower station; is the power generation flow for the th reservoir at the th interval; represents the reservoir water release corresponding to the th control cross-section; shows the different levels of ecological water requirements downstream from the reservoir corresponding to the th control cross-section.

### Constraints

## THE MODIFIED MULTI-OBJECTIVE SHUFFLED FROG LEAPING ALGORITHM (MMOSFLA)

### The traditional single-objective SFLA

Inspired by foraging behavior of frogs in the swamp, the SFLA is a metaheuristic optimization which adopts particle swarm intelligence (PSO) as a local search tool and the idea of competitiveness and mixing information from parallel local searches to gradually approach global optima from the shuffled complex evolution (SCE) (Duan *et al.* 1992).

The frog population is partitioned into a number of parallel sub-populations called memeplexes, of which frogs move toward the prey and exchange ideas independently, i.e. memetic evolution. Frogs in each sub-population or memeplex can be seen as potential feasible solutions and are impacted by ideas of other frogs (Eusuff & Lansey 2003; Zou *et al.* 2012). Information is exchanged between memeplexes in a shuffling process after a number of memetic evolutions, and the best ideas are passed down. The details of SFLA are shown below:

**Step 1.** The frog population initialization and objective function fitness of each frog calculation.

**Step 2.** Frogs are ranked in descending order according to the fitness quality.

**Step 3.** The frogs are grouped into *P* sub-populations which contain *n* frogs. In the first round, the frogs are divided into sub-populations one by one until the th frog is grouped in the th sub-population. Then, the ( + 1)th frog is grouped in the first sub-population and so on (Sun *et al.* 2016). The rest of the frogs adopt the same measure to accomplish the grouping.

**Step 4.**The local search is activated to update the worst frogs in all sub-populations in line with frog leap step and position equations: where the generates random number in the interval [0,1] which presents uniform distribution; is frog leap step, = 1,2,3 … ; and are the lower and upper bounds of the frog leap step, respectively; and correspond to the best and worst frog position in each sub-population, respectively; is the updated worst frog position.

**Step 5.** If the new position of the worst frog is better than before, the new position will eliminate the worst one. Otherwise, the in Equation (8) is replaced by the overall best frog position to produce a new frog position. If the new position is still incapable of improving the position quality, a random frog position is generated for the worst frog. After predefined iterations are finished in all sub-populations, the whole population is mixed in the shuffling process to realize the information exchange and transmission in all sub-populations or memeplexes (Luo *et al.* 2015a, 2015b; Yang *et al.* 2018).

### The multi-objective SFLA realization and modification

Conventional SFLA is incapable of handling multiple competing targets. Moreover, conventional multi-objective methods such as MOEAs, MOEAD and MOPSO may have deficiencies in optimization strategies. For example, the performance of most MOEAs is alleviated as the number of conflicting objectives increases. MOEAD guides the evolution based on the uniformly distributed weight vectors. However, a uniform weight vector on the unit hyper plane cannot guarantee the uniform distribution of solutions. For the MOPSO, the best local guide for each particle of the population from Pareto-optimal solutions is selected inefficiently. Thus, the concepts such as non-dominant relation, external elite frog set, and crowding distance are introduced to convert SFLA into MMOSFLA. In addition, improvements are proposed to further strengthen the capability of multi-objective problem optimization.

#### The population initialization based on chaos theory

The initialization process has significant effects on population diversity and quality. The population initialization in evolutionary multi-objective methods including the MOEA, MOPSO, MOGA and SFLA is accomplished through the random guided search strategy. However, random initialization makes algorithms difficult to quickly locate a feasible solution zone and consume a large amount of computer resources to eliminate worst solutions later.

*et al.*2016). It is confirmed by Cheng

*et al.*(2008) that the distribution diversity and homogeneity of the initial population can be improved by chaos initialization. To promptly locate a feasible search zone and improve computation efficiency, the modified chaotic logistic mapping developed in Equation (10) is applied to strengthen the overall population quality during the initialization process: where

*t*is the current th iteration; is the maximum iterations predefined;

*j*represents the chaotic variable dimension, i.e. the dimension of the optimization problem; variable is described as the th chaos variable after the

*t*th iteration, sliding in the interval [–1,1]. It distributes more uniformly in comparison to the conventional method.

The details of chaotic initialization are as follows:

Ensure the frog population scale; set = 0 and generate initial chaos variable .

Produce chaos variable (the number of variables is and represents the optimization target amount).

Calculate the objective function value of and make a judgment on the non-dominant relation to select the best th non-dominant solutions (chaos variable) for the initial frog population.

Terminate the chaos search until the predefined iterations are accomplished and obtain the final initial population.

#### Elite frog set maintenance and update

As the Pareto optimal solution set, the main function of the elite frog set is to save the non-dominated solutions in each iteration. The frogs in the elite frog set are not involved in population evolution. However, these frogs are devoted to guiding and correcting the evolution direction of the entire population. In the process of frog iteration, if the original frog and new elite frog display the non-dominant relationship, both frogs are included in the elite frog set for an alternative optimal solution. The details of update mechanism are as follows:

- 1.
The alternative frog will be added directly to the elite frog set in the case where that frog set is empty.

- 2.
The new alternative frog will be listed into the elite set if it dominates the current frogs in the set. The frogs which are dominated by the new frog will be deleted.

- 3.The elite frog set grows in size based on iterations. It is adverse to the search efficiency if the set scale is out of control. Therefore, we develop a crowding distance strategy to take control of the elite frog quantity. First, the frog set scale is defined as , and we calculate the crowding distance of frogs and rank them in descending order according to the crowding distance. Then, the redundant low-ranking frogs are eliminated from the set. Finally, the frogs with large crowding distance are retained to maintain diversity and distribution uniformity in the elite frog set. The crowding distance is determined by Equation (12): where denotes the crowding distance of the th elite frog. , and are fitness of –1,
*m*, + 1th frog corresponding to the th optimization objective, respectively. Analogously, the meanings of and are the best and worst fitness in the elite frog set.

To eliminate the frogs with lower crowding distance, all the alternative elite frogs are selected into the elite frog set and ranked by the crowding distance in the conventional method. Nevertheless, this method can effectively reduce counts of dominant comparisons, it may lead to the dense distribution of optimal solutions in the narrow space. Eventually, all frogs in this narrow space may be excluded during the elite frog set maintenance. As a result, we propose an improved maintenance strategy whereby alternative elite frogs are added into the set successively and the distribution density of frogs is evaluated. The frogs distributed intensively are eliminated after each addition. The new method proposed contributes to make optimal solutions uniformly distributed on the Pareto front.

#### The renewed grouping strategy

Assuming that represents the dominance relation of frogs in set , frogs are ranked in descending order by and random permutation is adopted if frogs have the same dominance relationship. Finally, all frogs will be assigned to corresponding sub-populations using conventional grouping strategy as mentioned. The method above shows defects on handling the Pareto dominance relationship between frogs which are not dominated by each other. The random permutation may cause extreme cases in the grouping process whereby most frogs assigned into sub-populations are the best or worst frogs. Accordingly, the population diversity significantly decreases and stagnation will be brought to the evolution process. For this purpose, the frog grouping method is modified as follows:

**Step 1.** Evaluate the fitness of each frog to determine the non-dominant frogs for the whole population. Rank these frogs in descending order according to the distribution density.

**Step 2**. For the rest of the dominant frogs, calculate its Euclidean distance from the nearest non-dominant frog and rank these frogs behind the excellent frogs in descending order.

**Step 3**. The method shown in Equation (14) is to further distinguish between frogs whose dominance relationship is difficult to confirm. This method determines the relationship by means of calculating the fitness offset of the frog. The frog with a smaller fitness offset will be selected for the non-dominant one. where is fitness of the th optimization objective corresponding to the th frog. represents the population scale.

*O*is the total number of objectives.

**Step 4**. Frogs are arranged into a respective sub-population (the number of sub-populations is ) successively with the method described above under ‘The traditional single-objective SFLA’. Apparently each sub-population contains frogs.

#### The modified local search strategy

The crowding distance and distance from the optimal Pareto front are introduced concerning the elite frog selection strategy. The elite frog is obtained from the two selection decisions to effectively conduct the evolution and improve the diversity of the frog population. Details are shown in the following steps:

**Step 1.**Define the fitness of as whose value is reciprocal of . The roulette wheel is introduced to choose the lower frog which is the elite frog with high probability. The fitness of is calculated by Equation (17):

**Step 2.** To avoid the excessively dense distribution of frogs selected in the former step, we continue to screen out the individuals with large crowding distance . Similarly, the roulette wheel method is used to locate the targeted frog with large probability.

The new frog obtained from Equations (15) and (16) is evaluated to determine the dominance relationship with the original worst one . The dominant frogs will be eliminated. Here are three cases for the elimination decision:

**Case 1:** If dominates , the former will replace the latter.

**Case 2:**If is dominated by , it indicates that the algorithm may be trapped into search stagnation because of a marked decrease in the population diversity. Therefore, to reactivate the search performance we need to seek more energetic and potential frogs to replace the current one. The new frog creation is achieved by Equation (18): where is the random number in [0,1]. denotes the adjustment coefficient which is changeable dynamically using the arc tangent function shown in Equation (19). The larger makes the algorithm escape from the local area and search within a larger space the initial iteration stage, while in the later iteration stage, a relatively lesser is used for maintaining the fast convergence characteristic of MMOSFLA: where and are the initial and terminal values of (0.9 and 0.4, respectively); represents the controlling factor which is in the range of [0.4,0.7].

*t*and are the current iterations for the local search in the sub-population and global shuffling.

*T*, show the maximum iterations for local search and global shuffling, respectively.

**Case 3:** If and do not dominate each other, it is equal (50% of the time) to make replace or keep unchanged.

After each local search iteration, frogs are reordered according to the dominant relationship. Repeat the above steps until the local iterations in the sub-population are accomplished.

#### Elite frog evolution strategy

The elite frog set plays an important role in guiding the evolution direction of the frog population by using the elite set maintenance strategy based on crowding distance. However, especially in the later iteration, the evolution will gradually slow down and trap into stagnation because of a decrease in the population diversity, even in the elite frog set. To this end, we explore the potential dominated frog to reactivate the search performance using the evolution strategy based on the normal cloud model theory. The normal cloud model first proposed by Li *et al*. (1995) has characteristics of randomness and stability orientation. Recently, the cloud model has been introduced and combined with a number of algorithms such as GA, PSO and an evolution algorithm (Dai *et al.* 2007; Zhang *et al.* 2008, 2012). The cloud model helps the algorithm to profoundly search in the local space and effectively improve overall performance.

The normal cloud model (NCM) is brought forth to describe the uncertain conversion relation of qualitative concept or qualitative knowledge with its quantitative expressions. The expectation , entropy and hyper-entropy are parameters for characterizing the digital features of the cloud (Ma *et al.* 2013). The details on , and are elaborated as follows and are depicted in Figure 1.

The is central to all droplets and the most representative cloud droplet for qualitative concept.

The entropy denotes the uncertainty measurement of the qualitative concept. It reflects emergence randomness of the cloud droplet and correlation of fuzziness and randomness. The scope of the cloud droplet generated and its randomness highly depends on the value of . Specifically, the larger is, the more obvious the randomness and broader the scope is.

The hyper-entropy , determined by randomness and fuzziness, represents the measurement of entropy, i.e. the entropy of entropy. To boost search randomness in the initial iteration stage and maintain search stability later, the value of is defined as /5.

The cloud droplet is created by the normal cloud generator. Usually, the processes can be described as follows (Sun *et al.* 2016):

**Step 1.** Generate a normal random number with expectation and standard deviation .

**Step 2.** Thereafter, generate a normal random number *x* as a cloud droplet, which is taken as , as expectation and standard deviation, respectively.

**Step 3.** The calculation and *x* are plugged into formula to determine certainty pertaining to the qualitative concept *C*.

**Step 4.** Repeat the above-mentioned processes until the total cloud droplets meet the terminal condition.

For the MMOSFLA, we select the excellent dominated frog of the elite set to execute evolution strategy in the later iteration stage. The location of the elite frog is defined as , and the variance of frogs in the elite set is which can dynamically change the search space. Furthermore, to enhance the search randomness in the early iteration stage and guarantee search stability in the later stage, we set the value of as /5. The specified number of cloud droplets is generated by the cloud generator, which represents potential solutions. Through the fitness evaluation, we establish a dominance relationship between new and original elite frogs. The original one will be replaced if the new cloud droplet (the frog) dominates it.

### The flow chart of MMOSFLA

The specific processes of MMOSFLA for solving MOEOCR problem are shown in Figure 2.

### Numerical simulations

Numerical simulations are conducted to demonstrate the performance of the proposed MMOSFLA approach. Simulation test functions are from the ZDT and DTLZ package, of which the optimum is ‘0’ theoretically. The NSGAII (fast elitist non-dominated sorting genetic algorithm), MOEAD (multi-objective evolutionary algorithm based on decomposition), MOEADDE (MOEAD based on differential evolution operators), dMOPSO (multi-objective particle based on decomposition) and MOPSO (multi-objective particle swarm optimization) are selected to compare with MMOSFLA. The frog population and sub-population size are set to 100 and 10, respectively. The number of local iterations in the sub-population and global shuffling iterations are defined as 100 and 10,000. The 30 independent simulations are conducted to offset randomness of the simulation results.

The iteration curves of the optimal Pareto front are partially shown in Figure 3. Furthermore, the IGD and GD indicators are analyzed to evaluate the performance of MMOSFLA quantitatively. The IGD is an indicator for assessing convergence and solution diversity. The GD represents the distance between the true Pareto front and the Pareto front searched by MMOSFLA. The smaller the IGD and GD are, the better the convergence and distribution of the solution. The comparative results of IGD and GD are demonstrated in Tables 1 and 2.

Test function . | Algorithm . | |||||
---|---|---|---|---|---|---|

NSGAII . | MOEAD . | dMOPSO . | MOEADDE . | MOPSO . | MMOSFLA . | |

ZDT1 | 2.3375 × 10^{–1} (1.16 × 10^{–1}) | 1.7506 × 10^{–1} (5.03 × 10^{–2}) | 2.6954 × 10^{–2} (1.11 × 10^{–2}) | 3.8550 × 10^{0} (2.28 × 10^{0}) | 4.2108 × 10^{1} (8.73 × 10^{0}) | 5.1653 × 10 ^{–3} (2.67 × 10^{–4}) |

ZDT2 | 5.4347 × 10^{–1} (1.40 × 10^{–1}) | 3.6358 × 10^{–1} (1.99 × 10^{–1}) | 4.9077 × 10^{–1} (2.38 × 10^{–1}) | 5.2994 × 10^{0} (2.41 × 10^{–1}) | 4.6156 × 10^{1} (1.08 × 10^{1}) | 2.3943 × 10^{–2} (7.18 × 10^{–2}) |

ZDT3 | 1.7938 × 10^{–1} (1.12 × 10^{–1}) | 2.3082 × 10^{–1} (1.05 × 10^{–1}) | 1.2439 × 10^{–1} (2.07 × 10^{–1}) | 4.5374 × 10^{0} (1.87 × 10^{0}) | 4.3630 × 10^{1} (1.00 × 10^{1}) | 2.2902 × 10 ^{–2} (9.42 × 10^{–2}) |

ZDT4 | 2.6914 × 10^{0} (2.24 × 10^{0}) | 4.6781 × 10^{–1} (2.28 × 10^{–1}) | 5.7986 × 10^{–1} (2.37 × 10^{–1}) | 3.5171 × 10^{0} (1.42 × 10^{0}) | 1.8812 × 10^{1} (7.62 × 10^{0}) | 2.5576 × 10 ^{–1} (1.52 × 10^{–1}) |

ZDT6 | 8.1677 × 10^{–2} (3.31 × 10^{–2}) | 9.3728 × 10^{–2} (5.84 × 10^{–2}) | 4.2025 × 10^{–3} (1.49 × 10^{–3}) | 4.7737 × 10^{–2} (9.57 × 10^{–2}) | 1.3558 × 10^{0} (2.58 × 10^{0}) | 3.5911 × 10 ^{–3} (2.03 × 10^{–4}) |

DTLZ1 | 2.1271 × 10 ^{–1} (1.77 × 10^{–1}) | 2.5625 × 10^{–1} (3.21 × 10^{–1}) | 7.2134 × 10^{0} (5.76 × 10^{0}) | 1.1982 × 10^{0} (2.06 × 10^{0}) | 1.0643 × 10^{1} (3.94 × 10^{0}) | 2.1795 × 10^{0} (1.70 × 10^{0}) |

DTLZ2 | 7.0013 × 10^{–2} (2.98 × 10^{–3}) | 5.4862 × 10 ^{–2} (1.84 × 10^{–4}) | 1.4315 × 10^{–1} (1.03 × 10^{–2}) | 7.7725 × 10^{–2} (1.62 × 10^{–3}) | 1.0528 × 10^{–1} (1.25 × 10^{–2}) | 7.2267 × 10^{–2} (2.42 × 10^{–3}) |

DTLZ3 | 8.8711 × 10 ^{0} (5.11 × 10^{0}) | 1.4796 × 10^{1} (9.72 × 10^{0}) | 4.7035 × 10^{1} (5.00 × 10^{1}) | 1.6801 × 10^{1} (1.81 × 10^{1}) | 1.4728 × 10^{2} (6.50 × 10^{1}) | 8.4056 × 10^{1} (3.13 × 10^{1}) |

DTLZ4 | 1.0072 × 10^{–1} (1.60 × 10^{–1}) | 3.8486 × 10^{–1} (3.21 × 10^{–1}) | 3.0915 × 10^{–1} (3.36 × 10^{–2} | 1.7320 × 10^{–1} (6.79 × 10^{–2}) | 3.4990 × 10^{–1} (1.42 × 10^{–1}) | 1.0009 × 10 ^{–1} (1.20 × 10^{–1}) |

DTLZ5 | 6.1237 × 10 ^{–3} (3.56 × 10^{–4}) | 3.2176 × 10^{–2} (1.05 × 10^{–3}) | 4.4669 × 10^{–2} (6.50 × 10^{–3}) | 1.4273 × 10^{–2} (3.35 × 10^{–4}) | 1.4772 × 10^{–2} (4.71 × 10^{–3}) | 6.5428 × 10^{–3} (5.30 × 10^{–4}) |

DTLZ6 | 6.8978 × 10^{–3} (8.34 × 10^{–4}) | 8.5158 × 10^{–2} (1.68 × 10^{–1}) | 3.2884 × 10^{–2} (8.37 × 10^{–4}) | 1.4065 × 10^{–2} (6.04 × 10^{–4}) | 2.8475 × 10^{0} (8.83 × 10^{–1}) | 5.9284 × 10 ^{–3} (3.58 × 10^{–4}) |

Test function . | Algorithm . | |||||
---|---|---|---|---|---|---|

NSGAII . | MOEAD . | dMOPSO . | MOEADDE . | MOPSO . | MMOSFLA . | |

ZDT1 | 2.3375 × 10^{–1} (1.16 × 10^{–1}) | 1.7506 × 10^{–1} (5.03 × 10^{–2}) | 2.6954 × 10^{–2} (1.11 × 10^{–2}) | 3.8550 × 10^{0} (2.28 × 10^{0}) | 4.2108 × 10^{1} (8.73 × 10^{0}) | 5.1653 × 10 ^{–3} (2.67 × 10^{–4}) |

ZDT2 | 5.4347 × 10^{–1} (1.40 × 10^{–1}) | 3.6358 × 10^{–1} (1.99 × 10^{–1}) | 4.9077 × 10^{–1} (2.38 × 10^{–1}) | 5.2994 × 10^{0} (2.41 × 10^{–1}) | 4.6156 × 10^{1} (1.08 × 10^{1}) | 2.3943 × 10^{–2} (7.18 × 10^{–2}) |

ZDT3 | 1.7938 × 10^{–1} (1.12 × 10^{–1}) | 2.3082 × 10^{–1} (1.05 × 10^{–1}) | 1.2439 × 10^{–1} (2.07 × 10^{–1}) | 4.5374 × 10^{0} (1.87 × 10^{0}) | 4.3630 × 10^{1} (1.00 × 10^{1}) | 2.2902 × 10 ^{–2} (9.42 × 10^{–2}) |

ZDT4 | 2.6914 × 10^{0} (2.24 × 10^{0}) | 4.6781 × 10^{–1} (2.28 × 10^{–1}) | 5.7986 × 10^{–1} (2.37 × 10^{–1}) | 3.5171 × 10^{0} (1.42 × 10^{0}) | 1.8812 × 10^{1} (7.62 × 10^{0}) | 2.5576 × 10 ^{–1} (1.52 × 10^{–1}) |

ZDT6 | 8.1677 × 10^{–2} (3.31 × 10^{–2}) | 9.3728 × 10^{–2} (5.84 × 10^{–2}) | 4.2025 × 10^{–3} (1.49 × 10^{–3}) | 4.7737 × 10^{–2} (9.57 × 10^{–2}) | 1.3558 × 10^{0} (2.58 × 10^{0}) | 3.5911 × 10 ^{–3} (2.03 × 10^{–4}) |

DTLZ1 | 2.1271 × 10 ^{–1} (1.77 × 10^{–1}) | 2.5625 × 10^{–1} (3.21 × 10^{–1}) | 7.2134 × 10^{0} (5.76 × 10^{0}) | 1.1982 × 10^{0} (2.06 × 10^{0}) | 1.0643 × 10^{1} (3.94 × 10^{0}) | 2.1795 × 10^{0} (1.70 × 10^{0}) |

DTLZ2 | 7.0013 × 10^{–2} (2.98 × 10^{–3}) | 5.4862 × 10 ^{–2} (1.84 × 10^{–4}) | 1.4315 × 10^{–1} (1.03 × 10^{–2}) | 7.7725 × 10^{–2} (1.62 × 10^{–3}) | 1.0528 × 10^{–1} (1.25 × 10^{–2}) | 7.2267 × 10^{–2} (2.42 × 10^{–3}) |

DTLZ3 | 8.8711 × 10 ^{0} (5.11 × 10^{0}) | 1.4796 × 10^{1} (9.72 × 10^{0}) | 4.7035 × 10^{1} (5.00 × 10^{1}) | 1.6801 × 10^{1} (1.81 × 10^{1}) | 1.4728 × 10^{2} (6.50 × 10^{1}) | 8.4056 × 10^{1} (3.13 × 10^{1}) |

DTLZ4 | 1.0072 × 10^{–1} (1.60 × 10^{–1}) | 3.8486 × 10^{–1} (3.21 × 10^{–1}) | 3.0915 × 10^{–1} (3.36 × 10^{–2} | 1.7320 × 10^{–1} (6.79 × 10^{–2}) | 3.4990 × 10^{–1} (1.42 × 10^{–1}) | 1.0009 × 10 ^{–1} (1.20 × 10^{–1}) |

DTLZ5 | 6.1237 × 10 ^{–3} (3.56 × 10^{–4}) | 3.2176 × 10^{–2} (1.05 × 10^{–3}) | 4.4669 × 10^{–2} (6.50 × 10^{–3}) | 1.4273 × 10^{–2} (3.35 × 10^{–4}) | 1.4772 × 10^{–2} (4.71 × 10^{–3}) | 6.5428 × 10^{–3} (5.30 × 10^{–4}) |

DTLZ6 | 6.8978 × 10^{–3} (8.34 × 10^{–4}) | 8.5158 × 10^{–2} (1.68 × 10^{–1}) | 3.2884 × 10^{–2} (8.37 × 10^{–4}) | 1.4065 × 10^{–2} (6.04 × 10^{–4}) | 2.8475 × 10^{0} (8.83 × 10^{–1}) | 5.9284 × 10 ^{–3} (3.58 × 10^{–4}) |

Data in brackets are standard deviation; optimal results are in bold font; methods compared are from Tian *et al.* (2017).

Test function . | Algorithm . | |||||
---|---|---|---|---|---|---|

NSGAII . | MOEAD . | dMOPSO . | MOEADDE . | MOPSO . | MMOSFLA . | |

ZDT1 | 1.1359 × 10^{–2} (6.66 × 10^{–3}) | 2.1418 × 10^{–2} (1.73 × 10^{–2}) | 2.8028 × 10^{–3} (1.18 × 10^{–3}) | 1.3408 × 10^{0} (6.29 × 10^{–1}) | 2.0180 × 10^{1} (1.94 × 10^{1}) | 2.4647 × 10 ^{–5} (2.15 × 10^{–5}) |

ZDT2 | 4.3380 × 10^{–2} (2.65 × 10^{–2}) | 1.6640 × 10^{–2} (1.86 × 10^{–2}) | 6.0982 × 10^{–4} (1.51 × 10^{–3}) | 1.7839 × 10^{0} (8.63 × 10^{–1}) | 4.2950 × 10^{1} (1.46 × 10^{1}) | 1.4220 × 10 ^{–5} (2.56 × 10^{–5}) |

ZDT3 | 8.2373 × 10^{–3} (4.32 × 10^{–3}) | 2.3681 × 10^{–2} (1.82 × 10^{–2}) | 1.6209 × 10^{–3} (9.81 × 10^{–4}) | 1.6517 × 10^{0} (6.22 × 10^{–1}) | 2.5155 × 10^{1} (2.06 × 10^{1}) | 4.0986 × 10 ^{–5} (1.50 × 10–5) |

ZDT4 | 2.3943 × 10 ^{–2} (2.32 × 10^{–2}) | 7.0046 × 10^{–2} (4.47 × 10^{–2}) | 2.4815 × 10^{–2} (1.30 × 10^{–1}) | 4.3727 × 10^{0} (3.60 × 10^{0}) | 8.2611 × 10^{0} (7.91 × 10^{0}) | 7.3716 × 10^{–1} (8.90 × 10^{–1}) |

ZDT6 | 1.2204 × 10^{–2} (6.18 × 10^{–3}) | 1.1560 × 10^{–2} (4.90 × 10^{–3}) | 9.4839 × 10^{–5} (7.32 × 10^{–4}) | 4.7398 × 10^{–2} (5.41 × 10^{–2}) | 3.4847 × 10^{–1} (4.91 × 10^{–1}) | 5.6006 × 10 ^{–5} (4.52 × 10^{–4}) |

DTLZ1 | 2.9925 × 10 ^{–2} (3.27 × 10^{–2}) | 4.1773 × 10^{–2} (5.70 × 10^{–2}) | 8.9482 × 10^{0} (1.67 × 10^{0}) | 1.0263 × 10^{0} (7.03 × 10^{–1}) | 7.6957 × 10^{0} (1.42 × 10^{0}) | 4.5638 × 10^{–1} (3.59 × 10^{–1}) |

DTLZ2 | 1.3245 × 10^{–3} (1.58 × 10^{–4}) | 6.2695 × 10 ^{–4} (4.89 × 10^{–5}) | 2.0143 × 10^{–2} (2.63 × 10^{–3}) | 1.7718 × 10^{–3} (3.48 × 10^{–4}) | 8.5144 × 10^{–3} (3.41 × 10^{–3}) | 1.7365 × 10^{–3} (5.93 × 10^{–4}) |

DTLZ3 | 1.4477 × 10 ^{0} (6.88 × 10^{–1}) | 2.5282 × 10^{0} (1.88 × 10^{0}) | 5.2792 × 10^{1} (8.67 × 10^{1}) | 5.3161 × 10^{0} (4.37 × 10^{0}) | 5.3308 × 10^{1} (2.19 × 10^{1}) | 1.1846 × 10^{1} (3.45 × 10^{0}) |

DTLZ4 | 1.2615 × 10^{–3} (3.15 × 10^{–4}) | 1.6947 × 10^{–3} (9.89 × 10^{–4}) | 5.5364 × 10^{–2} (1.67 × 10^{–2}) | 1.3925 × 10^{–3} (6.92 × 10^{–4}) | 3.4345 × 10^{–2} (1.76 × 10^{–2}) | 4.3859 × 10 ^{–4} (2.29 × 10^{–4}) |

DTLZ5 | 2.8789 × 10^{–4} (5.77 × 10^{–5}) | 2.8383 × 10^{–4} (1.06 × 10^{–3}) | 3.3932 × 10^{–2} (1.12 × 10^{–2}) | 3.1252 × 10^{–4} (1.12 × 10^{–4}) | 1.3698 × 10^{–3} (1.39 × 10^{–3}) | 2.0039 × 10 ^{–4} (6.57 × 10^{–5}) |

DTLZ6 | 5.0088 × 10^{–6} (2.83 × 10^{–7}) | 1.4951 × 10^{–2} (3.48 × 10^{–2}) | 1.4138 × 10^{–2} (7.79 × 10^{–3}) | 4.7962 × 10^{–6} (2.53 × 10^{–7}) | 4.2219 × 10^{–1} (7.51 × 10^{–2}) | 4.7883 × 10 ^{–6} (2.33 × 10^{–7}) |

Test function . | Algorithm . | |||||
---|---|---|---|---|---|---|

NSGAII . | MOEAD . | dMOPSO . | MOEADDE . | MOPSO . | MMOSFLA . | |

ZDT1 | 1.1359 × 10^{–2} (6.66 × 10^{–3}) | 2.1418 × 10^{–2} (1.73 × 10^{–2}) | 2.8028 × 10^{–3} (1.18 × 10^{–3}) | 1.3408 × 10^{0} (6.29 × 10^{–1}) | 2.0180 × 10^{1} (1.94 × 10^{1}) | 2.4647 × 10 ^{–5} (2.15 × 10^{–5}) |

ZDT2 | 4.3380 × 10^{–2} (2.65 × 10^{–2}) | 1.6640 × 10^{–2} (1.86 × 10^{–2}) | 6.0982 × 10^{–4} (1.51 × 10^{–3}) | 1.7839 × 10^{0} (8.63 × 10^{–1}) | 4.2950 × 10^{1} (1.46 × 10^{1}) | 1.4220 × 10 ^{–5} (2.56 × 10^{–5}) |

ZDT3 | 8.2373 × 10^{–3} (4.32 × 10^{–3}) | 2.3681 × 10^{–2} (1.82 × 10^{–2}) | 1.6209 × 10^{–3} (9.81 × 10^{–4}) | 1.6517 × 10^{0} (6.22 × 10^{–1}) | 2.5155 × 10^{1} (2.06 × 10^{1}) | 4.0986 × 10 ^{–5} (1.50 × 10–5) |

ZDT4 | 2.3943 × 10 ^{–2} (2.32 × 10^{–2}) | 7.0046 × 10^{–2} (4.47 × 10^{–2}) | 2.4815 × 10^{–2} (1.30 × 10^{–1}) | 4.3727 × 10^{0} (3.60 × 10^{0}) | 8.2611 × 10^{0} (7.91 × 10^{0}) | 7.3716 × 10^{–1} (8.90 × 10^{–1}) |

ZDT6 | 1.2204 × 10^{–2} (6.18 × 10^{–3}) | 1.1560 × 10^{–2} (4.90 × 10^{–3}) | 9.4839 × 10^{–5} (7.32 × 10^{–4}) | 4.7398 × 10^{–2} (5.41 × 10^{–2}) | 3.4847 × 10^{–1} (4.91 × 10^{–1}) | 5.6006 × 10 ^{–5} (4.52 × 10^{–4}) |

DTLZ1 | 2.9925 × 10 ^{–2} (3.27 × 10^{–2}) | 4.1773 × 10^{–2} (5.70 × 10^{–2}) | 8.9482 × 10^{0} (1.67 × 10^{0}) | 1.0263 × 10^{0} (7.03 × 10^{–1}) | 7.6957 × 10^{0} (1.42 × 10^{0}) | 4.5638 × 10^{–1} (3.59 × 10^{–1}) |

DTLZ2 | 1.3245 × 10^{–3} (1.58 × 10^{–4}) | 6.2695 × 10 ^{–4} (4.89 × 10^{–5}) | 2.0143 × 10^{–2} (2.63 × 10^{–3}) | 1.7718 × 10^{–3} (3.48 × 10^{–4}) | 8.5144 × 10^{–3} (3.41 × 10^{–3}) | 1.7365 × 10^{–3} (5.93 × 10^{–4}) |

DTLZ3 | 1.4477 × 10 ^{0} (6.88 × 10^{–1}) | 2.5282 × 10^{0} (1.88 × 10^{0}) | 5.2792 × 10^{1} (8.67 × 10^{1}) | 5.3161 × 10^{0} (4.37 × 10^{0}) | 5.3308 × 10^{1} (2.19 × 10^{1}) | 1.1846 × 10^{1} (3.45 × 10^{0}) |

DTLZ4 | 1.2615 × 10^{–3} (3.15 × 10^{–4}) | 1.6947 × 10^{–3} (9.89 × 10^{–4}) | 5.5364 × 10^{–2} (1.67 × 10^{–2}) | 1.3925 × 10^{–3} (6.92 × 10^{–4}) | 3.4345 × 10^{–2} (1.76 × 10^{–2}) | 4.3859 × 10 ^{–4} (2.29 × 10^{–4}) |

DTLZ5 | 2.8789 × 10^{–4} (5.77 × 10^{–5}) | 2.8383 × 10^{–4} (1.06 × 10^{–3}) | 3.3932 × 10^{–2} (1.12 × 10^{–2}) | 3.1252 × 10^{–4} (1.12 × 10^{–4}) | 1.3698 × 10^{–3} (1.39 × 10^{–3}) | 2.0039 × 10 ^{–4} (6.57 × 10^{–5}) |

DTLZ6 | 5.0088 × 10^{–6} (2.83 × 10^{–7}) | 1.4951 × 10^{–2} (3.48 × 10^{–2}) | 1.4138 × 10^{–2} (7.79 × 10^{–3}) | 4.7962 × 10^{–6} (2.53 × 10^{–7}) | 4.2219 × 10^{–1} (7.51 × 10^{–2}) | 4.7883 × 10 ^{–6} (2.33 × 10^{–7}) |

From inspection of Figure 3 and Tables 1 and 2, it is evident that the mean and STD of the IGD indicator for ZDT1-ZDT3 calculated by MMOSFLA are the best. Similarly, MMOSFLA also outperforms in DTLZ4 and DTLZ6 optimization. The results indicate the convergence performance and solution diversity by MMOSFLA is at a higher level than other approaches. For DTLZ2 and DTLZ5, the margin of optimization result is insignificant even if MMOSFLA is not the best. In terms of the GD indicator, MMOSFLA dominates other approaches in optimizing ZDT1-ZDT3, ZDT6 and DTLZ4-DTLZ6. It is clear that solutions gained by MMOSFLA are multiple and of higher precision (closer to true Pareto-front) in comparison to other methods.

In summary, MMOSFLA shows more excellent optimization performance than other algorithms. Even though IGD and GD indicators do not perform well in terms of several test functions, the final solutions on the Pareto front present good distribution uniformity and diversity. Therefore, we verify significant adaptability and effectiveness of MMOSFLA in complex nonlinear multi-objective problems.

## CASE STUDY

### The overview of cascade reservoirs located in Qingjiang River

Qingjiang River is the largest tributary of the middle Yangtze River between Yichang and Jingjiang. It is 423 km long with a basin area of 17,000 km^{2} which flows through Enshi, Badong, Changyang and Yidu located in Hubei province. Along the river, there are three hydropower stations with large hydro-unit capacity and reservoir regulation storage.

The Shuibuya Reservoir, characterized with multi-year regulating storage and the largest installed capacity of three hydropower stations, is the most important stage of the cascade development in Qingjiang River basin. The hydropower station of Geheyan, located in Changyang county, is a huge water conservancy project with the major function of power generation, flood control and shipping benefit. Gaobazhou Reservoir, which is situated in the lower reaches of the Qingjiang River, plays a role on the reverse regulation of the upstream Geheyan Reservoir. The main characteristics of the water level of cascade reservoirs and power station parameters are shown in Table 3,**.** The generalization of cascade reservoirs system in Qingjiang River catchment is illustrated in Figure 4.

Characteristic water level and parameters . | Cascade reservoirs in Qingjiang . | ||
---|---|---|---|

Shuibuya . | Geheyan . | Gaobazhou . | |

Normal storage water level (m) | 400 | 200 | 80 |

Flood control level (m) | 391.8 | 193.6 | 78.5 |

Dead water level (m) | 350 | 160 | 78.09 |

Total storage capacity (10^{8} m^{3}) | 45.89 | 30.18 | 4.03 |

Beneficial reservoir capacity (10^{8} m^{3}) | 23.83 | 19.75 | 0.54 |

Installed capacity (MW) | 1,840 | 1,212 | 270 |

Annual utilization hours of installed capacity (h) | 2,450 | 2,533 | 3,563 |

Guaranteed output (MW) | 310.0 | 241.0 | 77.3 |

Comprehensive efficiency coefficient | 8.50 | 8.50 | 8.40 |

Maximum water head (m) | 203.0 | 121.6 | 40.0 |

Minimum water head (m) | 147.0 | 80.7 | 22.3 |

Mean water head (m) | 186.6 | 112.0 | 35.4 |

Characteristic water level and parameters . | Cascade reservoirs in Qingjiang . | ||
---|---|---|---|

Shuibuya . | Geheyan . | Gaobazhou . | |

Normal storage water level (m) | 400 | 200 | 80 |

Flood control level (m) | 391.8 | 193.6 | 78.5 |

Dead water level (m) | 350 | 160 | 78.09 |

Total storage capacity (10^{8} m^{3}) | 45.89 | 30.18 | 4.03 |

Beneficial reservoir capacity (10^{8} m^{3}) | 23.83 | 19.75 | 0.54 |

Installed capacity (MW) | 1,840 | 1,212 | 270 |

Annual utilization hours of installed capacity (h) | 2,450 | 2,533 | 3,563 |

Guaranteed output (MW) | 310.0 | 241.0 | 77.3 |

Comprehensive efficiency coefficient | 8.50 | 8.50 | 8.40 |

Maximum water head (m) | 203.0 | 121.6 | 40.0 |

Minimum water head (m) | 147.0 | 80.7 | 22.3 |

Mean water head (m) | 186.6 | 112.0 | 35.4 |

### Ecological flow acquisition

The stream ecological base flow is defined as the basic water amount needed to maintain the stability of the biocenosis in the river, and protect the river's ecological environment and functions from damage. The concept of ecological or environmental flows was historically developed as a response to the degradation of aquatic ecosystems caused by human interventions (Tegos *et al.* 2018). In this paper, the annual runoff and ecological health extent is considered to determine the basic, suitable and ideal ecological flow of the control section located in Geheyan and Gaobazhou reservoirs.

*et al.*2018). The suitable and ideal ecological flow are obtained based on the requirement level index. The requirement level index can represent effects exerted by factors including vegetation coverage, biological diversity, and water quality on the ecological flow (Lyu

*et al.*2016). The specific method to calculate the requirement level index of suitable and ideal ecological flow is shown as follows: where , and denote the requirement level index of basic, suitable and ideal ecological flow, respectively. is the correction coefficient of suitable ecological flow; represents the corresponding characteristic index; and signify the membership degree of variation coefficient and the characteristic index of suitable ecological flow to the requirement level index of suitable ecological flow. The membership degree can be acquired through the entropy evaluation method. and are annual variation coefficient and characteristic index of ideal ecological flow, respectively. and signify the membership degree of and to the requirement level index of ideal ecological flow.

**:**where , denotes suitable ideal ecological flow of the control sections; is average monthly runoff.

Scheduling period . | Month . | Control section . | |||||
---|---|---|---|---|---|---|---|

Geheyan . | Gaobazhou . | ||||||

Basic . | Suitable . | Ideal . | Basic . | Suitable . | Ideal . | ||

Fish breeding | May | 93.08 | 229.63 | 363.41 | 77.86 | 214.04 | 336.92 |

June | 105.62 | 263.35 | 419.64 | 94.66 | 254.58 | 405.82 | |

July | 133.65 | 332.32 | 521.29 | 87.80 | 237.92 | 380.78 | |

August | 76.37 | 186.71 | 295.88 | 61.30 | 167.36 | 268.08 | |

September | 65.87 | 166.41 | 263.83 | 44.00 | 120.49 | 192.03 | |

October | 52.37 | 126.51 | 201.01 | 37.21 | 100.18 | 159.99 | |

Common water use (non-flood season) | November | 23.03 | 56.04 | 87.98 | 20.94 | 56.26 | 89.53 |

December | 10.85 | 26.34 | 41.73 | 19.08 | 52.24 | 83.64 | |

January | 8.09 | 19.96 | 31.46 | 21.90 | 58.46 | 93.18 | |

February | 11.72 | 27.62 | 47.94 | 19.31 | 51.41 | 81.86 | |

March | 22.34 | 54.79 | 86.94 | 28.34 | 76.47 | 122.01 | |

April | 41.06 | 99.86 | 156.48 | 37.41 | 100.87 | 161.90 |

Scheduling period . | Month . | Control section . | |||||
---|---|---|---|---|---|---|---|

Geheyan . | Gaobazhou . | ||||||

Basic . | Suitable . | Ideal . | Basic . | Suitable . | Ideal . | ||

Fish breeding | May | 93.08 | 229.63 | 363.41 | 77.86 | 214.04 | 336.92 |

June | 105.62 | 263.35 | 419.64 | 94.66 | 254.58 | 405.82 | |

July | 133.65 | 332.32 | 521.29 | 87.80 | 237.92 | 380.78 | |

August | 76.37 | 186.71 | 295.88 | 61.30 | 167.36 | 268.08 | |

September | 65.87 | 166.41 | 263.83 | 44.00 | 120.49 | 192.03 | |

October | 52.37 | 126.51 | 201.01 | 37.21 | 100.18 | 159.99 | |

Common water use (non-flood season) | November | 23.03 | 56.04 | 87.98 | 20.94 | 56.26 | 89.53 |

December | 10.85 | 26.34 | 41.73 | 19.08 | 52.24 | 83.64 | |

January | 8.09 | 19.96 | 31.46 | 21.90 | 58.46 | 93.18 | |

February | 11.72 | 27.62 | 47.94 | 19.31 | 51.41 | 81.86 | |

March | 22.34 | 54.79 | 86.94 | 28.34 | 76.47 | 122.01 | |

April | 41.06 | 99.86 | 156.48 | 37.41 | 100.87 | 161.90 |

### The scheduling scenarios

In this paper, two scheduling scenarios are proposed based mainly on power generation of the cascade reservoirs and the different ecological flow requirements of Geheyan and Gaobazhou control sections are considered simultaneously. The conflict relationship between power generation and ecological water spill and shortage under each ecological flow requirement are discussed. The specific scenarios are described as follows:

**Scenario 1:** The cascade reservoirs focus on power generation and take suitable ecological water requirements of Geheyan and Gaobazhou into consideration.

**Scenario 2:** The cascade reservoirs focus on power generation and take ideal ecological water requirements of Geheyan and Gaobazhou into consideration.

### The encoding method and constraint handling

#### The encoding method

#### The constraint handling

It is critical to effectively solve the MOEOCR problem with rational constraint handling. For this purpose, we propose the water level corridor combined with a penalty function to handle constraint.

For the reservoir water level constraint, the frog locations (water level) generated are restricted within a feasible zone. With regard to the discharge flow and output constraints which are an implicit function of reservoir water level, we can convert those constraints to a single water level constraint by use of a water level corridor. The specific conversion processes are as follows:

- 1.
If the final scheduling water level is undetermined, start forward recursion from the water level at the initial period until the final scheduling period. Then the upper and lower bound of water level at each period can be obtained. Next, we can build a positive corridor of water level constraint.

- 2.
By contrast, start reverse recursion from the water level at the final scheduling period until the initial period if the initial water level is unknown. Similarly, we can build a corresponding reverse corridor.

- 3.
The forward and reverse water level corridors are intersected to obtain a feasible zone of reservoir water level.

The convergent speed can be boosted when the search zone is narrowed into the water level corridor. However, the frog locations generated in the corridor are still possible to be unfeasible (Zhou *et al.* 2010). Therefore, the penalty function is used to further improve the frog quality. A larger or smaller value is given to the objective function to make feasible frogs generate.

### The main parameter settings for MMOSFLA

Fewer parameters are controlled in MMOSFLA, and the control parameters are selected as follows with contrastive analysis.

The total iterations of chaotic mapping initialization are set to 100, the frog population size *M* is 200 and the corresponding sub-population or memeplex size *P* is predefined as 20. The number of frogs in each sub-population *n* is 10 accordingly. The upper and lower bound of frog leap step is set to [–4,4]. Initial and terminal adjustment coefficient , = 0.90 and 0.40, controlling factor is 0.60. The number of local iterations in the sub-population and global shuffling iteration are set to 20 and 200, respectively.

## RESULTS AND DISCUSSION

In this paper, the inflow processes of the Qingjiang cascade reservoirs corresponding to the typical wet, normal and dry years are selected as the scheduling model input. Then, the MMOSFLA proposed is applied to optimize the multi-objective scheduling model for cascade reservoirs. The initial and terminal scheduling water level of Shuibuya, Geheyan and Gaobazhou reservoirs are 400, 200 and 80 m, respectively. The spatial distribution of multi-objective optimal scheduling solution set (Pareto front) corresponding to scenario 1 and 2 is displayed in Figure 5. Take the scheduling scenario 2 which considers the ideal ecological flow requirement of Geheyan and Gaobazhou for instance, the partial representative scheduling solutions are compared in Table 5.

The solution number . | The wet year . | The normal year . | The dry year . | |||
---|---|---|---|---|---|---|

Power generation (10^{8} kW·h)
. | Ecological water spill and shortage (10^{8} m^{3})
. | Power generation (10^{8} kW·h)
. | Ecological water spill and shortage (10^{8} m^{3})
. | Power generation (10^{8} kW·h)
. | Ecological water spill and shortage (10^{8} m^{3})
. | |

1 | 85.24 | 122.13 | 73.81 | 109.51 | 53.85 | 94.36 |

15 | 85.80 | 126.46 | 74.18 | 113.20 | 53.98 | 96.40 |

30 | 86.21 | 132.62 | 74.48 | 117.57 | 54.08 | 98.48 |

The solution number . | The wet year . | The normal year . | The dry year . | |||
---|---|---|---|---|---|---|

Power generation (10^{8} kW·h)
. | Ecological water spill and shortage (10^{8} m^{3})
. | Power generation (10^{8} kW·h)
. | Ecological water spill and shortage (10^{8} m^{3})
. | Power generation (10^{8} kW·h)
. | Ecological water spill and shortage (10^{8} m^{3})
. | |

1 | 85.24 | 122.13 | 73.81 | 109.51 | 53.85 | 94.36 |

15 | 85.80 | 126.46 | 74.18 | 113.20 | 53.98 | 96.40 |

30 | 86.21 | 132.62 | 74.48 | 117.57 | 54.08 | 98.48 |

In Figure 5 it is demonstrated that power generation in Qingjiang cascade reservoir and ecological water spill and shortage present an interactional and mutually contradictory relationship. It is manifest that maximum power generation and ecological benefits cannot be achieved simultaneously and the competition between the two goals is clear. Therefore, it is difficult to find a single optimal scheduling solution. Instead, the multi-objective solutions exist in the form of the Pareto front (trade-off) and the solutions acquired have better space distribution and diversity. Moreover, the power generation varies at the most by 1% within each Pareto front set of solutions in Figure 5. It seems that the problem can be solved as a single objective ecological problem with constrained fixed power generation. However, multiple calculations are required for this method to obtain a multi-objective non-inferior solution set. Besides, the weight setting between each objective and constraints handling are subjective which results in the large error between the true Pareto front and inferior solutions we obtain. Furthermore, in fact, it is difficult to pursue the fixed power generation for the hydropower station due to the long-term development and the gradual increase in electricity consumption. The decision makers need to evaluate the scheduling solution and select the best one within acceptable limits. For this purpose, it is necessary to provide a set of solutions for them to make a decision.

From Table 5 we can find that scheduling solution 1 is the most representative scheduling solution focusing on ecological benefit, by contrast solution 30 is aimed at the maximum power generation with less attention to ecological objective. As a trade-off, solution 15 contributes to maintaining the balance between the multiple objectives.

Furthermore, we demonstrate the scheduling solution set (shown in Table 6) of Geheyan Reservoir in scenario 2. The results in Table 6 similarly indicate the competitive relation between power generation and ecological benefit. The value ranges corresponding to the normal and dry years are less than that to the wet year. Hence, under circumstances of relative abundant reservoir inflow, the cascade reservoirs have more space to coordinate the optimal scheduling. The specific scheduling processes including water level, unit output and water release in Geheyan Reservoir corresponding to the typical years are presented in Figure 6.

Solution number . | Wet year . | Normal year . | Dry year . | |||
---|---|---|---|---|---|---|

Power generation (10^{8} kW·h)
. | Ecological water spill and shortage (10^{8} m³)
. | Power generation (10^{8} kW·h)
. | Ecological water spill and shortage (10^{8} m³)
. | Power generation (10^{8} kW·h)
. | Ecological water spill and shortage (10^{8} m³)
. | |

1 | 33.9825 | 55.4698 | 29.2178 | 40.7005 | 20.3696 | 23.3429 |

2 | 33.9990 | 55.6153 | 29.2288 | 40.7996 | 20.3773 | 23.3997 |

3 | 34.0164 | 55.7153 | 29.2394 | 40.8793 | 20.3847 | 23.4454 |

4 | 34.0336 | 55.8623 | 29.2504 | 40.9886 | 20.3923 | 23.5081 |

5 | 34.0498 | 55.9928 | 29.2614 | 41.0935 | 20.4000 | 23.5683 |

6 | 34.0624 | 56.1240 | 29.2678 | 41.1477 | 20.4084 | 23.5994 |

7 | 34.0792 | 56.2950 | 29.2780 | 41.2330 | 20.4155 | 23.6482 |

8 | 34.0958 | 56.4671 | 29.2889 | 41.3464 | 20.4231 | 23.7133 |

9 | 34.1114 | 56.6394 | 29.2996 | 41.4491 | 20.4306 | 23.7722 |

10 | 34.1277 | 56.8074 | 29.3103 | 41.5612 | 20.4381 | 23.8365 |

11 | 34.1442 | 56.9704 | 29.3207 | 41.6601 | 20.4453 | 23.8932 |

12 | 34.1598 | 57.0962 | 29.3320 | 41.7636 | 20.4532 | 23.9526 |

13 | 34.1761 | 57.2633 | 29.3416 | 41.8634 | 20.4598 | 24.0098 |

14 | 34.1908 | 57.4306 | 29.3524 | 41.9574 | 20.4674 | 24.0637 |

15 | 34.2050 | 57.6174 | 29.3630 | 42.0711 | 20.4748 | 24.1289 |

16 | 34.2179 | 57.7736 | 29.3712 | 42.1703 | 20.4805 | 24.1858 |

17 | 34.2300 | 57.9214 | 29.3805 | 42.2670 | 20.4870 | 24.2413 |

18 | 34.2352 | 58.0217 | 29.3838 | 42.3551 | 20.4938 | 24.3165 |

19 | 34.2473 | 58.1710 | 29.3928 | 42.4480 | 20.5001 | 24.3698 |

20 | 34.2590 | 58.3400 | 29.4021 | 42.5511 | 20.5066 | 24.4290 |

21 | 34.2707 | 58.4911 | 29.4112 | 42.6445 | 20.5129 | 24.4826 |

22 | 34.2819 | 58.6522 | 29.4194 | 42.7484 | 20.5186 | 24.5423 |

23 | 34.2930 | 58.8284 | 29.4279 | 42.8431 | 20.5245 | 24.5966 |

24 | 34.3060 | 58.9959 | 29.4355 | 42.9283 | 20.5299 | 24.6455 |

25 | 34.3188 | 59.2179 | 29.4435 | 43.0182 | 20.5354 | 24.6972 |

26 | 34.3313 | 59.4428 | 29.4514 | 43.1294 | 20.5410 | 24.7610 |

27 | 34.3397 | 59.6606 | 29.4586 | 43.2760 | 20.5459 | 24.8303 |

28 | 34.3478 | 59.8671 | 29.4668 | 43.4039 | 20.5515 | 24.9009 |

29 | 34.3546 | 60.0053 | 29.4750 | 43.5314 | 20.5572 | 24.9713 |

30 | 34.3691 | 60.2325 | 29.4837 | 43.6959 | 20.5633 | 25.0469 |

Solution number . | Wet year . | Normal year . | Dry year . | |||
---|---|---|---|---|---|---|

Power generation (10^{8} kW·h)
. | Ecological water spill and shortage (10^{8} m³)
. | Power generation (10^{8} kW·h)
. | Ecological water spill and shortage (10^{8} m³)
. | Power generation (10^{8} kW·h)
. | Ecological water spill and shortage (10^{8} m³)
. | |

1 | 33.9825 | 55.4698 | 29.2178 | 40.7005 | 20.3696 | 23.3429 |

2 | 33.9990 | 55.6153 | 29.2288 | 40.7996 | 20.3773 | 23.3997 |

3 | 34.0164 | 55.7153 | 29.2394 | 40.8793 | 20.3847 | 23.4454 |

4 | 34.0336 | 55.8623 | 29.2504 | 40.9886 | 20.3923 | 23.5081 |

5 | 34.0498 | 55.9928 | 29.2614 | 41.0935 | 20.4000 | 23.5683 |

6 | 34.0624 | 56.1240 | 29.2678 | 41.1477 | 20.4084 | 23.5994 |

7 | 34.0792 | 56.2950 | 29.2780 | 41.2330 | 20.4155 | 23.6482 |

8 | 34.0958 | 56.4671 | 29.2889 | 41.3464 | 20.4231 | 23.7133 |

9 | 34.1114 | 56.6394 | 29.2996 | 41.4491 | 20.4306 | 23.7722 |

10 | 34.1277 | 56.8074 | 29.3103 | 41.5612 | 20.4381 | 23.8365 |

11 | 34.1442 | 56.9704 | 29.3207 | 41.6601 | 20.4453 | 23.8932 |

12 | 34.1598 | 57.0962 | 29.3320 | 41.7636 | 20.4532 | 23.9526 |

13 | 34.1761 | 57.2633 | 29.3416 | 41.8634 | 20.4598 | 24.0098 |

14 | 34.1908 | 57.4306 | 29.3524 | 41.9574 | 20.4674 | 24.0637 |

15 | 34.2050 | 57.6174 | 29.3630 | 42.0711 | 20.4748 | 24.1289 |

16 | 34.2179 | 57.7736 | 29.3712 | 42.1703 | 20.4805 | 24.1858 |

17 | 34.2300 | 57.9214 | 29.3805 | 42.2670 | 20.4870 | 24.2413 |

18 | 34.2352 | 58.0217 | 29.3838 | 42.3551 | 20.4938 | 24.3165 |

19 | 34.2473 | 58.1710 | 29.3928 | 42.4480 | 20.5001 | 24.3698 |

20 | 34.2590 | 58.3400 | 29.4021 | 42.5511 | 20.5066 | 24.4290 |

21 | 34.2707 | 58.4911 | 29.4112 | 42.6445 | 20.5129 | 24.4826 |

22 | 34.2819 | 58.6522 | 29.4194 | 42.7484 | 20.5186 | 24.5423 |

23 | 34.2930 | 58.8284 | 29.4279 | 42.8431 | 20.5245 | 24.5966 |

24 | 34.3060 | 58.9959 | 29.4355 | 42.9283 | 20.5299 | 24.6455 |

25 | 34.3188 | 59.2179 | 29.4435 | 43.0182 | 20.5354 | 24.6972 |

26 | 34.3313 | 59.4428 | 29.4514 | 43.1294 | 20.5410 | 24.7610 |

27 | 34.3397 | 59.6606 | 29.4586 | 43.2760 | 20.5459 | 24.8303 |

28 | 34.3478 | 59.8671 | 29.4668 | 43.4039 | 20.5515 | 24.9009 |

29 | 34.3546 | 60.0053 | 29.4750 | 43.5314 | 20.5572 | 24.9713 |

30 | 34.3691 | 60.2325 | 29.4837 | 43.6959 | 20.5633 | 25.0469 |

It can be seen from Figure 6 that the optimization space of the water level in the major flood season (June and July) is smaller than that in other months due to the flood control task conducted by Geheyan Reservoir. The difference between each solution (Solutions 1, 15 and 30) is mainly concentrated in the dry season and water supply period after the flood season. The water release process of Geheyan Reservoir indicates that a relatively abundant natural inflow in typical wet year results in a large water spill in May before the major flood season. However, the ecological water shortage in the dry season and water supply periods at the end of the flood season can be alleviated. For the typical normal year, the ecological impact is mainly reflected in the form of water shortage in the later stage of the flood season and dry season, especially April and May, while in June, July, August and September the form is converted into water spill. As for the typical dry year, the ecological water spill presents a downward trend in the major flood season and water shortage intensifies after September.

More information can be acquired from analyzing scheduling solutions in Figure 6. Due to the emphasis on the ecological benefit, solution 1 increases the quantity of release and lowers the water level quickly to alleviate an ecological water shortage in the earlier stage of the dry season. In the period of water level falling, the water level need to intensely lower is declined and water release is cut down to avoid excessive ecological water spill. During the flood season, the risk of downstream ecology destruction is reduced with a peak flood mitigation effect of the reservoir. The above measures play an important role in protecting the downstream ecological health of Geheyan Reservoir. However, the quick falling of the water level make the reservoir operate at a lower water level. The power generation cannot be centralized with water level falling so the power generation benefit obviously decreases.

Solution 30 focuses on pursuing the power generation benefit. For this purpose, it condenses the reservoir water release to maintain a higher operating water level. Because of centralized power generation and water supply task conducted by the reservoir at the end of the flood season, the ecological water shortage gap inevitably amplifies. Ecological benefit is lost with less consideration for ecology protection in solution 30. Solution 15, by contrast, is a trade-off solution in terms of power generation and ecological benefits. It contributes to achieving balanced development for economic and ecological benefits.

To verify the advantages of MMOSFLA in the MOEOCR problem, the scheduling solution Pareto front optimized by MMOSFLA is compared with MOPSO and MOEAD. The parameter settings, including population scale and maximum iterations, remain the same. The independent simulations are set to 20. Similarly, we choose the scheduling solutions under scenario 2 for comparison. The Pareto front of power generation for the cascade reservoir and ecological water spill and shortage acquired by each method is displayed in Figure 7. In addition, the IGD and GD indexes are used to evaluate the solutions obtained by each algorithm. The results are listed in Tables 7 and 8.

Typical year . | The mean of IGD . | |||
---|---|---|---|---|

NSGAII . | dMOPSO . | MOEADDE . | MMOSFLA . | |

Wet year | 2.95 × 10^{–1} | 4.60 × 10^{–1} | 4.84 × 10^{0} | 2.77 × 10 ^{–1} |

Normal year | 7.04 × 10^{–1} | 9.19 × 10^{–1} | 8.18 × 10^{0} | 1.20 × 10 ^{–1} |

Dry year | 5.86 × 10^{–1} | 1.06 × 10^{0} | 2.33 × 10^{–1} | 2.09 × 10 ^{–1} |

STD of IGD . | ||||

Wet year | 3.05 × 10^{–2} | 1.70 × 10 ^{–3} | 8.23 × 10^{–2} | 2.03 × 10^{–3} |

Normal year | 6.10 × 10^{–3} | 3.40 × 10^{–3} | 2.16 × 10^{–1} | 1.08 × 10 ^{–3} |

Dry year | 8.47 × 10^{–3} | 4.12 × 10^{–3} | 3.53 × 10^{–1} | 1.88 × 10 ^{–3} |

Typical year . | The mean of IGD . | |||
---|---|---|---|---|

NSGAII . | dMOPSO . | MOEADDE . | MMOSFLA . | |

Wet year | 2.95 × 10^{–1} | 4.60 × 10^{–1} | 4.84 × 10^{0} | 2.77 × 10 ^{–1} |

Normal year | 7.04 × 10^{–1} | 9.19 × 10^{–1} | 8.18 × 10^{0} | 1.20 × 10 ^{–1} |

Dry year | 5.86 × 10^{–1} | 1.06 × 10^{0} | 2.33 × 10^{–1} | 2.09 × 10 ^{–1} |

STD of IGD . | ||||

Wet year | 3.05 × 10^{–2} | 1.70 × 10 ^{–3} | 8.23 × 10^{–2} | 2.03 × 10^{–3} |

Normal year | 6.10 × 10^{–3} | 3.40 × 10^{–3} | 2.16 × 10^{–1} | 1.08 × 10 ^{–3} |

Dry year | 8.47 × 10^{–3} | 4.12 × 10^{–3} | 3.53 × 10^{–1} | 1.88 × 10 ^{–3} |

Typical year . | The mean of GD . | |||
---|---|---|---|---|

NSGAII . | dMOPSO . | MOEADDE . | MMOSFLA . | |

Wet year | 1.22 × 10^{–1} | 9.24 × 10^{–2} | 4.42 × 10^{–1} | 5.60 × 10^{–2} |

Normal year | 4.66 × 10^{–2} | 1.80 × 10^{–1} | 3.18 × 10^{–1} | 2.75 × 10^{–2} |

Dry year | 6.00 × 10^{–2} | 2.07 × 10^{–1} | 1.68 × 10^{–1} | 3.53 × 10^{–2} |

The STD of GD . | ||||

Wet year | 6.15 × 10^{–3} | 7.32 × 10^{–2} | 5.51 × 10^{–1} | 4.52 × 10^{–3} |

Normal year | 1.23 × 10^{–2} | 6.46 × 10^{–2} | 6.19 × 10^{–1} | 4.43 × 10^{–3} |

Dry year | 1.74 × 10^{–2} | 9.77 × 10^{–2} | 1.01 × 10^{–1} | 7.86 × 10^{–3} |

Typical year . | The mean of GD . | |||
---|---|---|---|---|

NSGAII . | dMOPSO . | MOEADDE . | MMOSFLA . | |

Wet year | 1.22 × 10^{–1} | 9.24 × 10^{–2} | 4.42 × 10^{–1} | 5.60 × 10^{–2} |

Normal year | 4.66 × 10^{–2} | 1.80 × 10^{–1} | 3.18 × 10^{–1} | 2.75 × 10^{–2} |

Dry year | 6.00 × 10^{–2} | 2.07 × 10^{–1} | 1.68 × 10^{–1} | 3.53 × 10^{–2} |

The STD of GD . | ||||

Wet year | 6.15 × 10^{–3} | 7.32 × 10^{–2} | 5.51 × 10^{–1} | 4.52 × 10^{–3} |

Normal year | 1.23 × 10^{–2} | 6.46 × 10^{–2} | 6.19 × 10^{–1} | 4.43 × 10^{–3} |

Dry year | 1.74 × 10^{–2} | 9.77 × 10^{–2} | 1.01 × 10^{–1} | 7.86 × 10^{–3} |

On inspection of Figure 7 we find that the solutions that MMOSFLA converges to presents better diversity and distribution effect in comparison with other rivals (MOEAD, MOPSO). It is closer to the true Pareto front than others. According to the ecological controlling objective and measures issued by the administrative department, the reservoir decision maker should take ecology protection into consideration in the operation of the cascade reservoirs. For the decision maker or the end-users, it is difficult to sacrifice excessive economical benefits for ecology protection. The decision maker can choose an acceptable way where the reservoir is operated under given ecological water spill and shortage. For instance, the decision maker indicates that the ecological water spill and shortage should not exceed 112 (10^{8}m^{3}) in terms of the normal year. Accordingly, they can determine a feasible and acceptable reservoir scheduling solution acquired by each algorithm and, apparently, the solution found by MMOSFLA has superiority in respect of total power generation than the compared methods.

It is shown in Tables 7 and 8 that for the normal year, the mean of IGD corresponding to each algorithm are 7.04 × 10^{–1}, 9.19 × 10^{–1}, 8.18 × 10^{0} and 1.20 × 10^{–1}. Similarly, mean of GD corresponding to each algorithm are 4.66 × 10^{–2}, 1.80 × 10^{–1}, 3.18 × 10^{–1} and 2.75 × 10^{–2}. It can be concluded that MMOSFLA performs better than other algorithms. The lower IGD and GD demonstrate that MMOSFLA is closer to true Pareto and has advantages in the solution convergence, distribution and diversity compared with other methods. Although the STD of IGD in the wet year is 2.03 × 10^{–3} is inferior to the dMOPSO (1.70 × 10^{–3}), the difference is controlled in a narrow range and MMOSFLA outperforms others at both the IGD and GD evaluation index. Therefore, we have reason to believe that MMOSFLA is more suitable for solving the MOEOCR problem. The decision maker can pursue relatively more power generation with consideration of ecological benefit.

## CONCLUSIONS

In order to pursue higher economical and ecological benefits, multi-objective ecological operation for cascade reservoirs (MOEOCR) is conducted to give comprehensive consideration to water requirements in power generation and ecological health. For the high-dimensional MOEOCR problem with complex hydraulic constraints, the traditional single objective optimization method has difficulty demonstrating a mutually constrained relationship between objectives. For this purpose, we construct the multi-objective ecological operation model for cascade reservoirs and adopt the MMOSFLA proposed which converts SFLA into a multi-objective method to solve the model. The MMOSFLA introduces chaotic population initialization, renewed frog grouping method, local search method, and evolution for an elite frog set based on cloud model.

The test function results verify that MMOSFLA outperforms other methods at search capability, diversity and distribution uniformity for the Pareto solution set. For the ZDT3 and DTLZ 6, the mean and STD of IGD and GD obtained by MMOSFLA are 2.2902 × 10^{–2} (9.42 × 10^{–2}), 4.0986 × 10^{–5} (1.50 × 10^{–5}) and 5.9284 × 10^{–3} (3.58 × 10^{–4}), 4.7883 × 10^{–6} (2.33 × 10^{–7}), respectively. Finally, the proposed MMOSFLA is applied in an MOEOCR problem of hydropower stations in Qingjiang River. The results reflect a competing relationship between the power generation and ecological water spill and shortage. At the same time, with comparison with other methods, the solution optimized by MMOSFLA brings about more economical and ecological benefits conforming to suitable and ideal ecological water requirements. The less mean and STD of IGD and GD indexes demonstrate better solution distribution and diversity. It can provide a feasible scheduling solution set for decision makers to regulate a scheduling mode which coordinates cascade power generation and downstream ecological protection for Geheyan and Gaobazhou reservoirs. Therefore, we provide an effective measure for solving complex and high-dimensional MOEOCR problems.

Future study could emphasize on using stochastic rainfall and distributed run-off models to generate large amounts of inflow scenarios and test how the cascade reservoirs behave under all scenarios (Dimitriadis *et al.* 2016; Dimitriadis & Koutsoyiannis 2018). In addition, the expected Pareto front and its variability due to inflow uncertainty can also be estimated to determine more desirable reservoir operation rules.

## ACKNOWLEDGEMENTS

This research was funded by The National Key Basic Research Program of China (973 Program) (2012CB417006), and the National Science Support Plan Project of China (2009BAC56B03), the Research and Extension Project of Hydraulic Science and Technology in Shanxi Province ‘Study on the Key Technology for Joint Optimal Operation of Complex Multi-Reservoir System and Water Network’, the Science and Technology Project of Yunnan Water Resources Department ‘Comprehensive Water Saving and Unconventional Water Utilization Research’.