Development of an Optimal Model for the Xiluodu-Xiangjiaba Cascade Reservoir System Considering the Downstream Environmental Flow

: To explore the inﬂuence of the Xiluodu-Xiangjiaba cascade reservoir system on the appropriate environmental ﬂow (AEF) of the Jinsha River, a multiobjective optimal cascade reservoir model was established with the aim of maximizing power generation while minimizing the downstream degree of AEF alteration. The AEF was determined using the range of variability approach (RVA). The optimal model was solved using an improved version of NSGA-II called INSGA2-DS. Inﬂows in typical normal and dry years were selected for optimization. The results show that in a normal year, power generation can be increased by 1.28% compared with that under the current regular operation conditions by prioritizing the maximization of power generation, in which case the degree of AEF alteration will increase by 13.86%. In contrast, the degree of AEF alteration will decrease by 22.53% if ecological protection is prioritized, but power generation will decrease by 0.62%. Similarly, in a dry year, power generation can be increased by 1.76% compared with that under the current regular operation conditions to maximize economic beneﬁt, in which case, the degree of AEF alteration will increase by 4.95%. By contrast, the degree of AEF alteration can be decreased by 13.70% if the objective is AEF minimization, but power generation will decrease by 0.48%. These research results provide useful information for the formulation of ecological operation schemes involving cascade reservoirs on the Jinsha River.


Introduction
Among all renewable energy sources, hydropower is the most commonly used for electricity production, with the corresponding global capacity having increased by more than 30% between 2007 and 2015 [1,2]. In the past three decades, reservoirs and hydropower stations have been rapidly developed to provide clean energy and reduce greenhouse gas emissions in China [3,4]. One of the most prevalent ways to take full advantage of available hydropower resources is to build cascade reservoirs along a river [5,6]. Relative to the total revenue generated by the operation of independent reservoirs, synergistic benefits can be achieved through the coordinated operation of cascade reservoirs [7]. Generally, the main objective of cascade reservoirs is to maximize energy and revenue yields [8,9]. However, this strategy has an adverse impact on downstream river ecosystems. In addition to

Study Area
The lower reaches of the Jinsha River (from Panzhihua to Yibin) represent some of the most abundant river resources in China. Four cascade reservoirs, namely, the Wudongde (WDD), Baihetan (BHT), Xiluodu (XLD) and Xiangjiaba (XJB) Reservoirs, either are under construction or have been completed. (Figure 1). Xiluodu and Xiangjiaba were first-phase projects, for which impoundment began in 2013 and 2012, respectively. The installed capacities of these two hydropower stations are 12,600 MW and 6000 MW, respectively [34,35]. The hydraulic connection between the two reservoirs is tight, forming a cascade reservoir system with effective joint operation conditions, the main parameters of these two reservoirs are listed in Table 1. Baihetan (BHT), Xiluodu (XLD) and Xiangjiaba (XJB) Reservoirs, either are under construction or have been completed. (Figure 1). Xiluodu and Xiangjiaba were first-phase projects, for which impoundment began in 2013 and 2012, respectively. The installed capacities of these two hydropower stations are 12,600 MW and 6000 MW, respectively [34,35]. The hydraulic connection between the two reservoirs is tight, forming a cascade reservoir system with effective joint operation conditions, the main parameters of these two reservoirs are listed in Table 1. .

Panzhihua Yibin
Lower reaches of the Jinsha River Middle reaches of the Jinsha River . .  The objectives of the model were to maximize total hydropower generation and to minimize the degree of AEF alteration. These two objectives are mutually conflicting because maximizing hydropower production requires increasing the water level in the reservoir, while minimizing the degree of AEF alteration requires more water to be released from the reservoir for downstream needs [36] These two competing objectives of the system are expressed as follows.

XLD-XJB
(1) Maximizing the hydropower generation of the overall system 1 1 max ( ) (2) Minimizing the degree of AEF alteration where T is the total optimization period; t is the time step index; XLD t P and XJB t P are the outputs of the Xiluodu and Xiangjiaba hydropower stations, respectively, in the th t time step; t  is the  The objectives of the model were to maximize total hydropower generation and to minimize the degree of AEF alteration. These two objectives are mutually conflicting because maximizing hydropower production requires increasing the water level in the reservoir, while minimizing the degree of AEF alteration requires more water to be released from the reservoir for downstream needs [36] These two competing objectives of the system are expressed as follows.
(1) Maximizing the hydropower generation of the overall system (2) Minimizing the degree of AEF alteration where T is the total optimization period; t is the time step index; P XLD t and P XJB t are the outputs of the Xiluodu and Xiangjiaba hydropower stations, respectively, in the tth time step; ∆t is the length of a time step; Q XJB,out t is the release from the Xiangjiaba Reservoir in the tth time step; and Q AEF,t is the downstream AEF. The detailed solution process is described in Section 3.2.

Constraint Conditions
The operation constraints of a reservoir include constraints on the water level, water balance, output and outflow. The specific constraints considered in this study are as follows.
(1) Restriction on the reservoir water level: where Z i t , Z i,min t , and Z i,max t are the actual, minimum, and maximum water levels, respectively, of reservoir i in time step t.
(2) Water balance constraint: where V i t+1 and V i t are the capacities of reservoir i at the beginning of time steps t and t + 1, respectively; I i t is the inflow to reservoir i; QT i t is the outflow from the turbines; and QS i t is the spillway or outlet flow that cannot be utilized for power generation. Together, QT i t and QS i t reflect the total outflow from the reservoir.
(3) Spatial transition flow constraint: where IP i+1 t denotes the inflow between power stations i and i + 1 generated by precipitation. This equation indicates that the total inflow to reservoir i + 1 consists of the inflow generated between reservoirs i and i + 1, and the total outflow from reservoir i.
(4) Output constraint: P i,min necessary calculations and then to remove infeasible solutions by means of a penalty function [38,39]. Such a penalty function penalizes individuals based on violations of constraints. In this study, the penalty function approach was used to reduce the applicability of infeasible solutions in order to encourage the selection of feasible solutions. Initially, the extent of constraint violation is calculated as follows: where c 1 and c 2 represent the constraint violations, C is a constraint, and C min and C max represent the allowable constraint scope of each solution. Thus, c can be expressed as follows: where c represents the average extent of violation of constraint C. Then, the corresponding penalty value is calculated as follows: where P is the penalty value, β is the penalty parameter, N X is the number of genes, N C is the number of constraints, and c ji is the penalty function for the jth constraint on the ith variable. When the objective is to determine the minimum value, the objective function value is multiplied by the penalty value: When the objective is to determine the maximum value, the objective function value is divide by the penalty value.
F m = f m /P (14) In Equations (13) and (14), F m is the penalized objective function, and f m is the objective function. Constraints can be divided into three categories: boundary constraints, soft constraints, and hard constraints. Boundary constraints (Equation (3)) limit the search space for the design task. Soft constraints (Equations (6) and (7)) are enforced by means of penalties, which will slightly reduce the applicability of solutions with low levels of violation and significantly reduce the applicability of solutions with high levels of violation; consequently, such constraints are not strictly related to certain levels of violation. The degree of violation of soft constraints will guide the optimization algorithm towards the feasible solution space during the search process. The constant coefficient of a penalty term cannot be set to too large a value because doing so would make an optimal fitness value difficult to obtain.
During the calculation process, if various genes have reached their limiting conditions, then corresponding penalties must be imposed in accordance with the degrees of violation of these limits. Consequently, the fitness value will change rapidly. The penalty functions for the various limiting conditions are defined as described below.
When Q t,i < 0, the corresponding solution is impossible, and the penalty function associated with the violation of this physical condition involves multiplying the obtained negative release value by the appropriate penalty value: where c 1t represents the average extent of violation of the release constraint. ).c 2t (16) where c 2t represents the average extent of violation of the release constraint.
For power generation, we establish penalty functions based on the requirements that Q g t < Q g,min t and P t < P min t . These penalty functions are formulated as follows: where c XLD ), respectively, for the Xiangjiaba Reservoir. To maintain the ecological base flow demand in the downstream area, if the reservoir release is less than the ecological base flow demand, we apply the following penalty function: where c 5t represents the average extent of violation of the release constraint associated with the minimum downstream environmental flow. The penalty parameter β is the sum of all levels of constraint violation: The coefficients defined in the above functions are mainly subjective, and several trial-and-error processes must be conducted to effectively obtain suitable solutions.

Model Solution
The goal of a single-objective optimization problem is to obtain an optimal solution by searching for the minimum or maximum value of the objective function. However, in real life, most water resource optimization problems involve conflicting objectives, for which there is no efficient method of finding optimal solutions with multiple trade-offs. The goal of such a multiobjective problem is to find an optimal solution for all targets, also known as a Pareto-optimal solution.
The nondominated sorting genetic algorithm II (NSGA-II) is a fast elite multiobjective algorithm that follows the NSGA framework [40][41][42]. This method has been effectively implemented to address various water management issues. However, this algorithm has many deficiencies when solving for the Pareto boundaries of high-dimensional multivariable and complex nonlinear multiobjective problems, such as an inability to identify pseudo-nondominated solutions, an unreasonable crowding distance formula, and high computational complexity.
This paper uses an improved algorithm based on NSGA-II, namely, INSGA2-DS. The pseudodominated solutions in the nondominated solution set are selected by introducing the concept of a dominance intensity order, and a new crowding distance formula that considers the variance is adopted. For calculating the crowding distance, an adaptive elite retention strategy is introduced to control the scale of elite retention, and other means are also adopted to speed up the calculation efficiency. The details of the improved algorithm are introduced as follows.

Fast Nondominated Sorting
NSGA-II uses the fast nondominated sorting method to determine the nondominated level and sort by comparing the dominance relations between each pair of individuals. However, this method results in a large number of pseudo-nondominated solutions in the nondominated set, which reduces the convergence efficiency and solution quality. To remove these pseudo-nondominated solutions, a measure called the dominance intensity is introduced as follows: where ξ i is the dominance intensity of the ith member of the population, m is the number of subobjectives, and f max k and f min k are the maximum and minimum values, respectively, of the kth subobjective in the nondominated set.

Distribution Retention Method
The concept of distribution here refers to the degrees of similarity and difference along the real Pareto boundary, reflecting the uniformity and diversity of the solution set. NSGA-II considers only the distances between neighboring individuals. However, individuals with different degrees of crowding distance with respect to different subobjectives cannot yield more genetic opportunities, and this situation is not conducive to maintaining the solution distribution. Therefore, a new crowding distance formula that considers the variance is introduced as follows: where i D is the crowding distance obtained with the new method, i d is the crowding distance in NSGA-II, and f i+1 j and f i−1 j are the values of the jth subobjective for individuals i+1 and i-1, respectively. After the above calculations, each individual member of the population is associated with three parameters: a nondominated set rank (rank), a dominance intensity (ξ) and a crowding distance (i D ). The priorities of these parameters are ordered as follows: rank (smaller is better) > i D (larger is better) > ξ (smaller is better).

Adaptive Elite Retention Strategy
The elite retention scale of NSGA-II is a fixed value, which is not conducive to solution set convergence. In the improved algorithm, an adaptive elite retention strategy is adopted. In the early stage of evolution, the elite retention scale is set to increase the diversity of the solutions. In the later stage of evolution, the elite retention scale is set to improve the convergence of the algorithm. The number of elites in the population of the gth generation, E g , is calculated as follows: where α g is the impact factor for the number of gth-generation elites and P is the population size. The value of α g is adaptively and iteratively calculated as follows: where α g+1 is the impact factor for the number of (g+1)th-generation elites and ρ is the ratio of the number of solutions in the nondominated set in the gth-generation population to the total population size P.
The improved algorithm, INSGA2-DS, consists of five main operations: initialization, crossover and mutation, fast nondominated sorting, implementation of elite crowding comparison operators, and termination and output of results.
The loop structure of INSGA2-DS is shown in Figure 2. P t is the parent population, and P c and P m are the offspring populations. F 1 , F 2 , · · · , F k are solution classes. F 1 is the best solution set among the parents and offspring, F 2 is the second-best solution set, etc. The approach for using INSGA2-DS to solve the problem considered in this paper is illustrated in Figure 3.

Parameter Setting
For this study, the INSGA2-DS optimization model was programmed based on MATLAB. The population size was set to 30, the maximum number of iterations was 1500, the crossover probability was 0.75, the mutation probability was 0.01, and 30 Pareto-optimal schemes were obtained through the corresponding calculations. The approach for using INSGA2-DS to solve the problem considered in this paper is illustrated in Figure 3. The approach for using INSGA2-DS to solve the problem considered in this paper is illustrated in Figure 3.

Parameter Setting
For this study, the INSGA2-DS optimization model was programmed based on MATLAB. The population size was set to 30, the maximum number of iterations was 1500, the crossover probability was 0.75, the mutation probability was 0.01, and 30 Pareto-optimal schemes were obtained through the corresponding calculations.

Parameter Setting
For this study, the INSGA2-DS optimization model was programmed based on MATLAB. The population size was set to 30, the maximum number of iterations was 1500, the crossover probability was 0.75, the mutation probability was 0.01, and 30 Pareto-optimal schemes were obtained through the corresponding calculations.

Optimization Period and Typical Years
In accordance with the typical operation and regulation of cascade reservoirs, the Xiluodu Reservoir operates at a limiting water level of 560.0 m during the flood season (from June 1 to September 10); then, it begins to store water on September 10, and the water level reaches 600.0 m at the end of September. The Xiangjiaba Reservoir operates at 370.0 m during the flood season and then also starts to store water on September 10 and the water level reaches 380.0 m at the end of September. However, recent research has shown that the flood control capacity of cascade reservoirs may somewhat decrease from the end of August to early September. Therefore, to improve the reservoir storage fill rate and reduce the impact of water storage downstream, the impoundment time should be advanced to the end of August.
Considering both the downstream environmental flow requirements and the calculation time, we divided the data to be analyzed into subperiods. Reservoir managers often use cycles of 5 days when making operation plans during the impoundment period. In terms of calculation efficiency, either 5 or 6 days is a suitable time interval; therefore, the entire period of interest was divided into 8 subperiods, with the first subperiod being Aug. 21-25, the second subperiod being Aug. 26-31, the third subperiod being Sep. 1-5, and so on.
The optimal results strongly depend on the choice of the typical hydrological year, as inflows often vary widely among different regulatory years. Typical years include wet, normal and dry years. In this study, the inflow to the Xiluodu Reservoir in wet years is very high, in which case impoundment has little effect on the downstream environmental flow; thus, we selected a normal year and dry year as the typical years of interest. The Pingshan hydrological station is a controlled station in the lower reaches of the Jinsha River (Figure 1). Based on the daily average flow data of Pingshan station from 1956 to 2007 (52 years), the total flows during the impoundment period for every year were accumulated and ranked in descending order (Figure 4). Considering both the downstream environmental flow requirements and the calculation time, we divided the data to be analyzed into subperiods. Reservoir managers often use cycles of 5 days when making operation plans during the impoundment period. In terms of calculation efficiency, either 5 or 6 days is a suitable time interval; therefore, the entire period of interest was divided into 8 subperiods, with the first subperiod being Aug. 21-25, the second subperiod being Aug. 26-31, the third subperiod being Sep. 1-5, and so on.
The optimal results strongly depend on the choice of the typical hydrological year, as inflows often vary widely among different regulatory years. Typical years include wet, normal and dry years. In this study, the inflow to the Xiluodu Reservoir in wet years is very high, in which case impoundment has little effect on the downstream environmental flow; thus, we selected a normal year and dry year as the typical years of interest. The Pingshan hydrological station is a controlled station in the lower reaches of the Jinsha River (Figure 1). Based on the daily average flow data of Pingshan station from 1956 to 2007 (52 years), the total flows during the impoundment period for every year were accumulated and ranked in descending order (Figure 4).  According to the definition of typical hydrological years in hydrology, the top 25% of the ranked total flows correspond to wet years, the total flows ranked between the top 25% and 75% correspond to normal years, and the total flows ranked below the top 75% correspond to dry years. Based on this principle, the typical dry year (ranked at 75%) and the typical normal year (ranked at 50%) were selected from among the dry years and the normal years, respectively. The inflows to the Xiluodu Reservoir in the corresponding typical years are shown in Figure 5.  According to the definition of typical hydrological years in hydrology, the top 25% of the ranked total flows correspond to wet years, the total flows ranked between the top 25% and 75% correspond to normal years, and the total flows ranked below the top 75% correspond to dry years. Based on this principle, the typical dry year (ranked at 75%) and the typical normal year (ranked at 50%) were selected from among the dry years and the normal years, respectively. The inflows to the Xiluodu Reservoir in the corresponding typical years are shown in Figure 5.
According to the definition of typical hydrological years in hydrology, the top 25% of the ranked total flows correspond to wet years, the total flows ranked between the top 25% and 75% correspond to normal years, and the total flows ranked below the top 75% correspond to dry years. Based on this principle, the typical dry year (ranked at 75%) and the typical normal year (ranked at 50%) were selected from among the dry years and the normal years, respectively. The inflows to the Xiluodu Reservoir in the corresponding typical years are shown in Figure 5.

Determination of the Appropriate Environmental Flow
The AEF concept was historically developed as a response to the degradation of aquatic ecosystems caused by human intervention [43,44]. The AEF is a crucial phenomenon in terms of electricity production and sustaining river integrity and is a key driver of ecological processes. Any flow disturbance may significantly affect and alter fluvial ecosystem dynamics [45,46].
More than 200 methods of determining environmental flows have been developed, and they can be classified into four principal groups: (1) hydrological methods [47,48], including the Tennant, 7Q10, flow duration curve (FDC), aquatic biological base flow, and RVA methods [49]; (2) hydraulic methods [50,51], including the wetted perimeter and R2CROSS methods; (3) habitat simulation methods [52], including the instream flow incremental methodology (IFIM), the physical habitat method, the weighted effective width method, and the minimum demand method for a biological space; and (4) holistic methods, including the building block method and benchmarking.
Currently, most AEFs obtained via either hydrological or hydraulic methods are calculated as fixed values which do not reflect the temporal variations in the river flow demand. For example, in the FDC method, only the flow value at a cumulative frequency of either 90% or 95% on the daily flow process curve is taken as the AEF, which is an approach that lacks practical rationality. To obtain an AEF that considers the temporal variations in demand, in this paper, the RVA method was used to calculate the AEF for the lower reaches of the Jinsha River. We calculated the AEF based on the runoff data collected at the Pingshan station from 1956 to 2007. The RVA method is a comprehensive method for assessing changes in river hydrological regimes by considering the size, frequency and variation characteristics of river flows. The calculation formula is as follows: where Q AEF,t is the AEF in time step t, Q mean,t is the mean annual flow in time step t, and Q max,t and Q min,t are the maximum and minimum threshold values, respectively, in time step t as calculated using the RVA method. The formulas for calculating the maximum and minimum threshold values are as follows: where η is a percentage that is generally equal to 17%; Q median,t is the median flow in time step t; µ is an integer greater than 0, usually taken to be 1.0; and Q sd,t is the standard variance corresponding to the average flow in time step t. The AEF results calculated are plotted in Figure 6.
t m e d i a n tm e a n t s d t t m e d i a n t m e a n t s d t where  is a percentage that is generally equal to 17%;

Optimization Results for a Normal Year
For the selected normal year, under the current regular operation policy, the total power generated by the cascade reservoir system during the impoundment period is 16.011 billion kW·h, and the degree of AEF alteration is 11.76%. The Pareto-optimal solution distribution for this normal year is shown in Figure 7. Three typical optimal schemes were selected for analysis of the reservoir water levels and outflows: scheme 1 (minimum degree of AEF alteration), scheme 2 (an unbiased optimal solution) and scheme 3 (maximum power generation). The red circle represents the reference scheme (regular operation). To achieve the maximum power generation, we can choose scheme 3, in which power generation is increased by 1.28% at the cost of increasing the degree of AEF alteration by 13.86%. To achieve the maximum ecological benefit, scheme 1 can be selected, in which the degree of AEF alteration is reduced by 22.53% with a 0.62% reduction in power generation. If there is no demand preference, we can choose scheme 2; while reducing the degree of AEF alteration by 6.38%, this scheme can also increase power generation by 0.89%. The details of these optimal results are shown in Table 2.

Optimization Results for a Normal Year
For the selected normal year, under the current regular operation policy, the total power generated by the cascade reservoir system during the impoundment period is 16.011 billion kW•h, and the degree of AEF alteration is 11.76%. The Pareto-optimal solution distribution for this normal year is shown in Figure 7. Three typical optimal schemes were selected for analysis of the reservoir water levels and outflows: scheme 1 (minimum degree of AEF alteration), scheme 2 (an unbiased optimal solution) and scheme 3 (maximum power generation). The red circle represents the reference scheme (regular operation). To achieve the maximum power generation, we can choose scheme 3, in which power generation is increased by 1.28% at the cost of increasing the degree of AEF alteration by 13.86%. To achieve the maximum ecological benefit, scheme 1 can be selected, in which the degree of AEF alteration is reduced by 22.53% with a 0.62% reduction in power generation. If there is no demand preference, we can choose scheme 2; while reducing the degree of AEF alteration by 6.38%, this scheme can also increase power generation by 0.89%. The details of these optimal results are shown in Table 2.      Figure 8 shows the cascade reservoir water levels and outflows for these three optimal schemes in a normal year. Under all three optimal schemes, the water storage capacity in the first two periods (intervals 1 and 2) is increased compared with that under regular operation, and the water level is increased earlier to enhance power generation. The water level of the Xiangjiaba Reservoir in the fourth period (interval 4) reaches 379.86 m and 379.85 m under schemes 3 and 2, respectively, and then remains at 379.80-379.90 m until the end of the water storage period under scheme 3. Considering the downstream environmental flow, under scheme 2, the outflow of the Xiangjiaba Reservoir is increased at the beginning of the seventh period (interval 7). To ensure that the downstream environmental water demand is met, in scheme 1, the outflow of the Xiangjiaba Reservoir is maintained as close as possible to the AEF; thus, the reservoir water level is slowly increased in the fourth and fifth periods. The Xiluodu Reservoir water level reaches 595.37 m in the 5th period; then, the outflow is increased in the seventh and eighth periods to meet the downstream environmental flow demand. Reservoir is maintained as close as possible to the AEF; thus, the reservoir water level is slowly increased in the fourth and fifth periods. The Xiluodu Reservoir water level reaches 595.37 m in the 5th period; then, the outflow is increased in the seventh and eighth periods to meet the downstream environmental flow demand.

Optimization Results for a Dry Year
For the selected typical dry year, under the regular operation scheme, the total power generated by the cascade reservoir system during the impoundment period is 12.156 billion kW•h, and the degree of AEF alteration is 27.30%. The Pareto-optimal distribution for this dry year is shown in Figure 9. Again, three typical optimal schemes are selected, namely, scheme 1 (minimum degree of AEF alteration), scheme 2 (unbiased optimal solution), and scheme 3 (maximum power generation), to analyze the reservoir water levels and outflows. To maximize power generation, we can choose scheme 3, which increases power generation by 1.76% while increasing the degree of AEF alteration by 4.95%. To minimize the degree of AEF alteration, we can choose scheme 1, which reduces the degree of AEF alteration by 13.70% while decreasing power generation by 0.48%. If there is no demand preference, scheme 2 can increase power generation by 1.08% while decreasing the degree

Optimization Results for a Dry Year
For the selected typical dry year, under the regular operation scheme, the total power generated by the cascade reservoir system during the impoundment period is 12.156 billion kW·h, and the degree of AEF alteration is 27.30%. The Pareto-optimal distribution for this dry year is shown in Figure 9. Again, three typical optimal schemes are selected, namely, scheme 1 (minimum degree of AEF alteration), scheme 2 (unbiased optimal solution), and scheme 3 (maximum power generation), to analyze the reservoir water levels and outflows. To maximize power generation, we can choose scheme 3, which increases power generation by 1.76% while increasing the degree of AEF alteration by 4.95%. To minimize the degree of AEF alteration, we can choose scheme 1, which reduces the degree of AEF alteration by 13.70% while decreasing power generation by 0.48%. If there is no demand preference, scheme 2 can increase power generation by 1.08% while decreasing the degree of AEF alteration by 7.51%. The details of the optimal results are shown in Table 3.   Figure 10 shows the cascade reservoir water levels and outflows for these three optimal schemes in a dry year. Due to the small inflow in dry years, the differences between the optimal schemes are smaller than those for a normal year. Under scheme 3, which prioritizes raising the water level, the Xiangjiaba Reservoir water level reaches 379.83 m in the fourth period and then remains above 379.00 m. The Xiluodu Reservoir water level reaches 599.36 m in the sixth period, and the water level is then decreased to 594.52 m in the seventh period to maintain a hydraulic head difference for power generation. In scheme 1, the outflow of the Xiluodu Reservoir begins to be increased in the fifth period to meet the demand for the downstream environmental flow. Compared to those in scheme 3, the average outflows of the Xiluodu Reservoir in the fifth and sixth periods in scheme 1 are increased by 775.46 m 3 /s and 817.27 m 3 /s, respectively. In both scheme 1 and scheme 2, the water level of the Xiluodu Reservoir is reduced in the seventh period so that the outflow is as close as possible to the AEF of the downstream river. In scheme 1, the average outflow to Xiangjiaba in the seventh period is 8481.55 m 3 /s, whereas in scheme 2, it is 8616.98 m 3 /s to meet the downstream AEF demand in the dry season.   Figure 10 shows the cascade reservoir water levels and outflows for these three optimal schemes in a dry year. Due to the small inflow in dry years, the differences between the optimal schemes are smaller than those for a normal year. Under scheme 3, which prioritizes raising the water level, the Xiangjiaba Reservoir water level reaches 379.83 m in the fourth period and then remains above 379.00 m. The Xiluodu Reservoir water level reaches 599.36 m in the sixth period, and the water level is then decreased to 594.52 m in the seventh period to maintain a hydraulic head difference for power generation. In scheme 1, the outflow of the Xiluodu Reservoir begins to be increased in the fifth period to meet the demand for the downstream environmental flow. Compared to those in scheme 3, the average outflows of the Xiluodu Reservoir in the fifth and sixth periods in scheme 1 are increased by 775.46 m 3 /s and 817.27 m 3 /s, respectively. In both scheme 1 and scheme 2, the water level of the Xiluodu Reservoir is reduced in the seventh period so that the outflow is as close as possible to the AEF of the downstream river. In scheme 1, the average outflow to Xiangjiaba in the seventh period is 8481.55 m 3 /s, whereas in scheme 2, it is 8616.98 m 3 /s to meet the downstream AEF demand in the dry season.
to those in scheme 3, the average outflows of the Xiluodu Reservoir in the fifth and sixth periods in scheme 1 are increased by 775.46 m 3 /s and 817.27 m 3 /s, respectively. In both scheme 1 and scheme 2, the water level of the Xiluodu Reservoir is reduced in the seventh period so that the outflow is as close as possible to the AEF of the downstream river. In scheme 1, the average outflow to Xiangjiaba in the seventh period is 8481.55 m 3 /s, whereas in scheme 2, it is 8616.98 m 3 /s to meet the downstream AEF demand in the dry season.

Discussion
In the above analysis, operation optimization was performed for a typical normal hydrological year and a typical dry year, and three different optimal schemes and their levels of optimization were investigated for each typical year. For the typical normal year, the differences among the different optimal schemes are more obvious due to the larger inflow to the reservoirs. When the maximum power generation is sought, the water level of each reservoir increases quickly and then, once the highest water level is reached, it is maintained at this level for power generation. Therefore, the range of outflow fluctuation is large, leading to an increase in the environmental flow. By contrast, when the environmental flow is prioritized to minimize the degree of AEF alteration, the stability of the outflow must be maintained, resulting in a relatively slow rise in the reservoir water levels and, consequently, a reduction in the total power generation of the cascade reservoir system.
For the typical dry year, because the incoming flow is generally lower than the environmental flow, the differences among the various optimal schemes are smaller than those in the typical normal year. In the regular operation scheme, the reservoir water levels rise slowly. When the power generation is sought, the corresponding optimal scheme still gives priority to increasing the reservoir water levels as soon as possible, whereas when the goal is to minimize the degree of AEF alteration, the corresponding optimal scheme is designed to improve the ecological benefit of the cascade reservoir system by reducing the water levels during a certain period. In general, the optimal schemes obtained here are reasonably consistent with actual operations. During real operation, the manager can choose a suitable scheme in accordance with specific needs.

Conclusions
This paper establishes an optimal model for the impoundment period of the Xiangjiaba-Xiluodu cascade reservoir system considering both economic and ecological benefits, where the AEF of the downstream river is determined using the RVA method. In the proposed method, the improved INSGA2-DS algorithm is used as a viable means of determining the optimal reservoir scheduling by balancing objectives related to human needs, water management and river ecological protection. When the optimization is focused on maximizing economic benefit in a normal hydrological year and a dry year, the level of power generation can be increased by 1.28% and 1.76%, respectively, but the degree of change in the AEF is simultaneously increased by 13.86% and 4.95%, respectively (compared to the current regular operation scheme). Conversely, when the optimization is focused on ecological needs, the degree of AEF alteration can be reduced by 22.53%

Discussion
In the above analysis, operation optimization was performed for a typical normal hydrological year and a typical dry year, and three different optimal schemes and their levels of optimization were investigated for each typical year. For the typical normal year, the differences among the different optimal schemes are more obvious due to the larger inflow to the reservoirs. When the maximum power generation is sought, the water level of each reservoir increases quickly and then, once the highest water level is reached, it is maintained at this level for power generation. Therefore, the range of outflow fluctuation is large, leading to an increase in the environmental flow. By contrast, when the environmental flow is prioritized to minimize the degree of AEF alteration, the stability of the outflow must be maintained, resulting in a relatively slow rise in the reservoir water levels and, consequently, a reduction in the total power generation of the cascade reservoir system.
For the typical dry year, because the incoming flow is generally lower than the environmental flow, the differences among the various optimal schemes are smaller than those in the typical normal year. In the regular operation scheme, the reservoir water levels rise slowly. When the power generation is sought, the corresponding optimal scheme still gives priority to increasing the reservoir water levels as soon as possible, whereas when the goal is to minimize the degree of AEF alteration, the corresponding optimal scheme is designed to improve the ecological benefit of the cascade reservoir system by reducing the water levels during a certain period. In general, the optimal schemes obtained here are reasonably consistent with actual operations. During real operation, the manager can choose a suitable scheme in accordance with specific needs.

Conclusions
This paper establishes an optimal model for the impoundment period of the Xiangjiaba-Xiluodu cascade reservoir system considering both economic and ecological benefits, where the AEF of the downstream river is determined using the RVA method. In the proposed method, the improved INSGA2-DS algorithm is used as a viable means of determining the optimal reservoir scheduling by balancing objectives related to human needs, water management and river ecological protection.
When the optimization is focused on maximizing economic benefit in a normal hydrological year and a dry year, the level of power generation can be increased by 1.28% and 1.76%, respectively, but the degree of change in the AEF is simultaneously increased by 13.86% and 4.95%, respectively (compared to the current regular operation scheme). Conversely, when the optimization is focused on ecological needs, the degree of AEF alteration can be reduced by 22.53% and 13.70%, respectively, while power generation is reduced by 0.62% and 0.48%, respectively. In addition, the wide distribution of the Pareto fronts produced using the proposed optimization techniques provides decision makers with considerable flexibility in establishing appropriate rules for reservoir operations by balancing human and ecological needs. The results not only verify the effectiveness of the proposed method but also provide a quantitative relationship between power generation and the degree of AEF alteration as well as operational references.