Multi-Objective Optimization for Reservoir Operation Considering Water Diversion and Power Generation Objectives

: Due to the uneven distribution of water resources in time and space, the problem of water shortage has become increasingly serious in some areas. To optimize use of water resources, it is urgent to establish multi-objective models and apply effective optimization algorithms to guide reservoir management. This study proposed a model of multi-objective optimization for reservoir operation (MORO) with the objectives of maximizing water diversion and power generation. The multi-objective evolutionary algorithm based on decomposition with adaptive weight vector adjustment (MOEA/D-AWA) was applied to solve the MORO problem. In addition, the performance of the MOEA/D-AWA was compared with two other algorithms based on the hyper-volume index. Huangjinxia reservoir, which is located in Shaanxi, China, was selected as the case study. The results show that: (1) the proposed model is effective and reasonable in theory; (2) the optimization results obtained by MOEA/D-AWA demonstrate this algorithm can be applied to the MORO problem, providing a set of evenly distributed non-dominated solutions; and (3) water diversion and power generation are indeed contradictory objectives. The MORO strategy can be used to efﬁciently utilize water resources, improve the comprehensive beneﬁts of reservoirs, and provide decision support for actual reservoir operation.


Introduction
With the progress of society, economic development, and the growth of the population, the problems of the water shortage and the imbalance of water supply and demand have become increasingly prominent. Human activity and land use can affect water resources. The effects can be explained by connectivity thinking: that is, human activities directly affect connectivity by land cover changes and different land use practices altering sediment delivery from hillslopes and floodplains to the channel system [1][2][3][4][5]. Water resource systems are very complicated, and the benefits of water resources should be maximized in terms of management and utilization of water resources. Thus, it is necessary to engage in management actions to fulfill consumer needs to achieve sustainable development and maximize benefits, especially in areas where the imbalance problem is serious [6]. Reservoirs are one of the most efficient key infrastructure components in water resources development and management [6,7]. Since water conservation projects often perform many tasks, such as flood control, power generation, meeting ecological needs, and providing water supply, etc., a single-objective operation based solely on water supply or power generation can no longer meet actual needs. Therefore, it has become increasingly obvious that the needs of multiple objectives and MOEA/D-AWA is an improved algorithm for MOEA/D, the results obtained by MOEA/D-AWA were compared with the results of NSGA-II and MOEA/D, and algorithm performance was measured based on the Hyper-volume index [35][36][37]. The Huangjinxia reservoir, which is located in Shaanxi, China, was the focus of the case study. It's meaningful for improving the comprehensive benefits of water conservancy projects to verify whether the proposed model is successful and the introduced algorithm has potential to solve the MORO problem. Better models and appropriate algorithms can help ensure that water resources are used efficiently and reasonably, and can effectively guide decisions on actual reservoir operation.
The rest of this paper is organized as follows. Section 2 proposes and establishes the reservoir operation optimization model, and considers some constraints. Section 3 introduces the MOEA/D-AWA algorithm, and presents the basic idea and the procedure of application of MOEA/D-AWA for the MORO problem. Section 4 describes the Huangjinxia reservoir as the focus of the case study, presenting the parameters and data of the model and the algorithm. Section 5 compares the performance of three algorithms, and analyzes and discusses the rationality of the operation results. Section 6 concludes the paper.

Reservoir Operation Optimization Model
Generally, reservoirs can be used to irrigation, municipal and industrial water supply, hydropower generation, flood control, and other uses. Therefore, the management of a multi-objective reservoir is complex due to conflicts between these objectives [38]. Reservoir optimization operation can be used for formulating, analyzing, and solving operation optimization problems in water resources planning and management [39]. The MORO problem aims to maximize water resources benefits as much as possible by determining the optimal release for each reservoir over the entire operation period, while satisfying appropriate physical and operational constraints [40]. The objective functions and associated constraints of the MORO problem studied in the paper can be expressed as follows.

Objective Functions
In order to rationally and effectively use the water resources of a reservoir, a multi-objective optimization model is required that considers both water diversion and power generation objectives for reservoirs with water diversion and power generation tasks [9]. The monthly water levels are taken as the decision variables. Due to the fact that the water storage is reflected by the water level, the fluctuation of the water level represents the variation of the outflow [7]. The amount of water for water diversion and power generation is dynamically determined by the water level for each period, and also affects the water level of the reservoir in turn. Once the upstream water level is determined, the amount of available water can be determined through the water balance equation shown in Equation (3). Simultaneously, the downstream water level is connected with the outflow of the reservoir. In terms of the hydroelectric power station, the increase of the outflow will result in a higher downstream water level, thus leading to a lower hydraulic head. In contrast, as for the pump station, the increase of the outflow will result in a lower upstream water level, thus leading to a higher hydraulic head for water diversion. Thus, the two objective functions have obvious conflict. These two objective functions are expressed as follows: 1. Maximum water diversion: where Q wd (t) is the water diversion flow; ∆T(t) is the time step; t is the index of operational periods; and T is the number of operational periods in one year.
where K is the power coefficient; Q f d (t) is the release flow passing through the turbines; H(t) is the hydraulic head of a hydropower plant.

Constraints
The above objective functions of the model on reservoir optimization operation are subject to the following constraints: (1) Water balance constraint: where V(t + 1) and V(t) are the final and initial storage volumes during period; respectively; I(t) is the inflow into the reservoir during each period ∆T(t); and S(t) is the amount of the loss of the reservoir during each period ∆T(t).
(2) Storage of the reservoir constraint: where V min (t) and V max (t) are the minimum and maximum storage allowed in each period, respectively. The minimum storage V min (t) is considered equal to the dead storage. V max (t) is constrained by flood protection rules during the flood periods when the maximum storage V max (t) is equivalent to the corresponding storage capacity of the flood-limited level of the reservoir. In the remaining periods, the maximum storage V max (t) is equivalent to the utilizable capacity of the reservoir.
(3) Power output constraint: where N min and N max are the minimum and maximum power output of the hydropower station, respectively.
(4) Release constraint: where Q min and Q max are the minimum and maximum turbine capacity, respectively.
(5) Pumping station design flow constraint: where Q design is the pumping station design flow.

Methodology
With the above model and constraints, the appropriate optimization algorithm is applied to optimize the MORO problem. In this study, MOEA/D-AWA was adopted to the MORO problem with the proposed operation model, and the monthly water level values of the reservoir were selected as the decision variables.

MOEA/D-AWA
MOEA/D-AWA was adopted to address the MORO problem to obtain a uniformly distributed optimal solution set on the Pareto front of the objectives. MOEA/D-AWA has the same framework as MOEA/D, but is enhanced with inclusion of a new weight vector initialization method and an adaptive weight vector adjustment (AWA) strategy. Here, the new weight vector initialization method is based on the geometric analysis of the Chebyshev decomposition approach, and the adaptive weight vector adjustment (AWA) strategy is to deal with the multi-objective optimization problems with complex Pareto fronts [33]. A more detailed description of MOEA/D-AWA is written in Qi et al. (2014) [33].
The simplified flow chart of the MOEA/D-AWA is shown in Figure 1. In the initialization stage, the data, parameters and studied problem are input. After the utilizing the simulated binary crossover (SBX) operator, polynomial mutation operator, adaptive weight adjustment, and other steps, the stopping condition is finally satisfied and the optimal solutions are obtained. Here, the generation number is termed as gen, the external population is termed as EP, the evolving population is termed as evol_pop, and the adaptive weight vector adjustment strategy is termed as wag.
Water 2018, 10, x FOR PEER REVIEW 5 of 15

MOEA/D-AWA
MOEA/D-AWA was adopted to address the MORO problem to obtain a uniformly distributed optimal solution set on the Pareto front of the objectives. MOEA/D-AWA has the same framework as MOEA/D, but is enhanced with inclusion of a new weight vector initialization method and an adaptive weight vector adjustment (AWA) strategy. Here, the new weight vector initialization method is based on the geometric analysis of the Chebyshev decomposition approach, and the adaptive weight vector adjustment (AWA) strategy is to deal with the multi-objective optimization problems with complex Pareto fronts [33]. A more detailed description of MOEA/D-AWA is written in Qi et al. (2014) [33].
The simplified flow chart of the MOEA/D-AWA is shown in Figure 1. In the initialization stage, the data, parameters and studied problem are input. After the utilizing the simulated binary crossover (SBX) operator, polynomial mutation operator, adaptive weight adjustment, and other steps, the stopping condition is finally satisfied and the optimal solutions are obtained. Here, the generation number is termed as gen , the external population is termed as EP , the evolving population is termed as _ evol pop , and the adaptive weight vector adjustment strategy is termed as wag .

Application of MOEA/D-AWA for the MORO Problem
Based on the basic ideas and basic steps of MOEA/D-AWA, the MOEA/D-AWA algorithm is applied to the MORO problem. The established functions and constraints of the MORO problem are added into the algorithm, and the established MORO model is used as the evaluation function. The population size is N, the maximum iteration time Maxgen is 800, and the operation periods are 12 months. Water levels are taken as decision variables, and N individuals are randomly generated within the feasible ranges of water levels (the upper and lower limits of water levels). Under certain constraints, the quality of the solutions is evaluated according to the proposed objective functions. Through a series of steps, mainly including allocation of computing resources, evolution and adaptive weight adjustment, constantly update the solutions, so repeatedly, until the stopping criterion is satisfied. The steps of MOEA/D-AWA are shown in Algorithm 1.
Choose other N 5 − m indices using 10-tournament selection according to π i , and add them to I; Step3: Evolution For each i ∈ I, do 3.1 {indices r 1 , indices r 2 , indices r 3 }←Selection (B(i)); 3.2 y←Genetic Operator (r 1 , r 2 , r 3 ); 3.3 Repair y; 3.4 Update the reference z←Compare f i (y) with z j * ; 3.5 Update the solutions ←Compare g y λ j , z * ≤ g x j λ j , z * ; Step4: Adaptive Weight Adjustment If gen ≥ rate_evol * Maxgen and genmodwag = 0, adaptively adjust the weight vectors as follows: 4.

Study Area
Shaanxi province in China has the two main drainage systems of the Yellow River and Yangtze River, bounded by the Qinling Mountains. As the largest tributary of the Yangtze River, the Hanjiang River originates from the Qinling Mountains in Shaanxi, China, with a total length of 1577 km, covering 159,000 km 2 . The average annual precipitation of the river basin is 700-1100 mm, of which the upstream is about 800-1200 mm. The runoff of the Hanjiang River is mainly supplied by precipitation, that mainly occurs in May-October. The Hanjiang river has abundant water resources, with 58.2 billion m 3 . The Weihe River, originating from the Gansu province in China and passing through the Shaanxi province in China, is the largest tributary of the Yellow River in China [41,42]. The total length of the Weihe River is 818 km and its drainage area is 135,000 km 2 , accounting for 17.9% of the total amount of the Yellow River basin. The main stream of the Weihe River is in Shaanxi, with a total length of 502.4 km and a drainage area of 67,108 km 2 . The average annual runoff of the Weihe River is 7.57 billion m 3 , and that of Shaanxi is 5.38 billion m 3 . The Weihe River plays an important role in Shaanxi, but its runoff has decreased significantly in recent decades [42].
The Hanjiang to Weihe River Water Diversion Project is located on the upper reach of the Hanjiang River in Shaanxi, China. The water diversion project is designed to divert water from the Hanjiang River with relatively abundant water resources into the Weihe River with severely deficient water resources. The main infrastructures of the water diversion project include two reservoirs (Huangjinxia reservoir and Sanhekou reservoir), two hydro plants, and two pumping stations.
Because the Huangjinxia reservoir with abundant inflow is the main source of water for Hanjiang to Weihe River Water Diversion Project, the case study focus is the Huangjinxia reservoir. The reservoir is located on the mainstream of the Hanjiang River, as shown in Figure 2. The Huangjinxia reservoir is a multi-objective reservoir, mainly taking on water supply, and also taking on power generation, flood control, and shipping. The reservoir supplies water to the Weihe River using the pumping station with a pumping capacity of 70 m 3 /s. Moreover, the activity of water diversion is carried out by direct intake from the reservoir area by the pumping station. The reservoir will be dispatched out water from July to June next year, and the water level will be stored at a limited water level during the flood periods, otherwise, it will be stored at the normal water level. The variational water level will drop to the dead water level at the end of June next year. The basic parameters of the Huangjinxia reservoir, hydro plant, and pump station are presented in Table 1.

Study Area
Shaanxi province in China has the two main drainage systems of the Yellow River and Yangtze River, bounded by the Qinling Mountains. As the largest tributary of the Yangtze River, the Hanjiang River originates from the Qinling Mountains in Shaanxi, China, with a total length of 1577 km, covering 159,000 km 2 . The average annual precipitation of the river basin is 700-1100 mm, of which the upstream is about 800-1200 mm. The runoff of the Hanjiang River is mainly supplied by precipitation, that mainly occurs in May-October. The Hanjiang river has abundant water resources, with 58.2 billion m 3 . The Weihe River, originating from the Gansu province in China and passing through the Shaanxi province in China, is the largest tributary of the Yellow River in China [41,42]. The total length of the Weihe River is 818 km and its drainage area is 135,000 km 2 , accounting for 17.9% of the total amount of the Yellow River basin. The main stream of the Weihe River is in Shaanxi, with a total length of 502.4 km and a drainage area of 67,108 km 2 . The average annual runoff of the Weihe River is 7.57 billion m 3 , and that of Shaanxi is 5.38 billion m 3 . The Weihe River plays an important role in Shaanxi, but its runoff has decreased significantly in recent decades [42].
The Hanjiang to Weihe River Water Diversion Project is located on the upper reach of the Hanjiang River in Shaanxi, China. The water diversion project is designed to divert water from the Hanjiang River with relatively abundant water resources into the Weihe River with severely deficient water resources. The main infrastructures of the water diversion project include two reservoirs (Huangjinxia reservoir and Sanhekou reservoir), two hydro plants, and two pumping stations.
Because the Huangjinxia reservoir with abundant inflow is the main source of water for Hanjiang to Weihe River Water Diversion Project, the case study focus is the Huangjinxia reservoir. The reservoir is located on the mainstream of the Hanjiang River, as shown in Figure 2. The Huangjinxia reservoir is a multi-objective reservoir, mainly taking on water supply, and also taking on power generation, flood control, and shipping. The reservoir supplies water to the Weihe River using the pumping station with a pumping capacity of 70 m 3 /s. Moreover, the activity of water diversion is carried out by direct intake from the reservoir area by the pumping station. The reservoir will be dispatched out water from July to June next year, and the water level will be stored at a limited water level during the flood periods, otherwise, it will be stored at the normal water level. The variational water level will drop to the dead water level at the end of June next year. The basic parameters of the Huangjinxia reservoir, hydro plant, and pump station are presented in Table 1.

Model and Algorithm Application
A monthly time scale was selected with a total time period of 12 months. The operation period was from July to June next year. The inflow data was selected with a probability of 75%.
In the study, the MOEAD-AWA algorithm was introduced to solve the MORO problem. The results obtained by MOEAD-AWA algorithm were compared with two other multi-objective evolutionary algorithms, MOEA/D and NSGA-II. The optimization problem included 13 decision variables, which were the 13 water level values of the reservoir for the 12 time periods, including the starting water level and the ending water level.
MOEA/D-AWA, MOEA/D, and NSGA-II all used a simulated binary crossover operator and polynomial mutation. After many experiments and comparisons, the following parameters were adopted: the three compared algorithms had the same population size of 100, the neighborhood size was set as 20 for MOEA/D-AWA and MOEA/D; the maximum iteration time was 800. The other parameter settings of the three algorithms are listed in Table 2, where n is the number of variables. The single-objective problem is to find an optimal and unique point in a space of decision variables, while a multi-objective optimization problem is not likely to have a unique solution that optimizes all the objectives, but instead has a set of optimal solutions. The basic parameters of the reservoir, hydropower station and pumping station and monthly inflow data are input, and this study proposes a model with two objectives, including maximizing water diversion and maximizing power generation, and various constraints are used as boundary conditions. MOEA/D-AWA is adopted to obtain optimal solutions. The algorithm parameters are set, and the optimal results are finally obtained until the stop condition is reached. The results of MOEA/D-AWA are compared with the results of NSGA-II and MOEA/D. The results include the Pareto fronts and the Hyper-volume curves of the three algorithms, the variable values, the performance evaluation results and the comparison results. Figure 3 shows three basic steps for the application of the mathematical optimization model and the introduced MOEA/D-AWA algorithm.

Results
The effectiveness of the proposed model and the feasibility of the introduced MOEA/D-AWA algorithm could be verified, by comparing the performance of different algorithms and analyzing the rationality of the operation results.

Results
The effectiveness of the proposed model and the feasibility of the introduced MOEA/D-AWA algorithm could be verified, by comparing the performance of different algorithms and analyzing the rationality of the operation results. Figure 4 showed the Pareto fronts of MOEA/D, NSGA-II, and MOEA/D-AWA algorithms for the MORO problem after 800 iterations. Overall, the trends of the three algorithms were roughly similar. However, as can be seen in Figure 4a, the Pareto front obtained by MOEA/D was discontinuously distributed. The non-dominated solutions obtained by MOEA/D were mainly concentrated at the bottom, but were discontinuously distributed at the top. As shown in Figures 4b and 4c, NSGA-II and MOEA/D-AWA performed significantly better than MOEA/D in terms of uniformity, and even the Pareto front distribution of MOEA/D-AWA was slightly more uniform than NSGA-II, due to slightly less uniformity at the top and bottom of NSGA-II's Pareto front compared to that of MOEA/D-AWA.

Analysis of the Pareto distribution of different algorithms
Water 2018, 10, x FOR PEER REVIEW 10 of 15 Figure 4 showed the Pareto fronts of MOEA/D, NSGA-II, and MOEA/D-AWA algorithms for the MORO problem after 800 iterations. Overall, the trends of the three algorithms were roughly similar. However, as can be seen in Figure 4a, the Pareto front obtained by MOEA/D was discontinuously distributed. The non-dominated solutions obtained by MOEA/D were mainly concentrated at the bottom, but were discontinuously distributed at the top. As shown in Figure 4b and Figure 4c, NSGA-II and MOEA/D-AWA performed significantly better than MOEA/D in terms of uniformity, and even the Pareto front distribution of MOEA/D-AWA was slightly more uniform than NSGA-II, due to slightly less uniformity at the top and bottom of NSGA-II's Pareto front compared to that of MOEA/D-AWA.

Analysis of convergence speed of different algorithms
The hyper-volume (HV) is a comprehensive index to measure the quality of multi-objective optimization algorithms. Especially when the ideal Pareto front is unknown, the HV can objectively evaluate the convergence, breadth, and uniformity of the Pareto front [35,36]. For the minimization problem of two-objective optimization, the larger the value of the HV, the better the convergence, breadth, and uniformity of the Pareto front. Additionally, a high HV value means that the nondominated solution set is a good approximation of the unknown Pareto front [37].
As shown in Figure 5, with increased number of iterations, the convergence curves of HV of three algorithms became more stable. However, there were differences in the convergence speed and HV values of the three algorithms. The curve of MOEA/D started to stabilize after about 150 iterations, that of MOEA/D-AWA gradually stabilized after about 250 iterations, and the convergence rate for NSGA-II was slower than that of the other two. Before 600 iterations, the stability of the distribution of the non-dominated solution set of NSGA-II was relatively poor, and it gradually stabilized after 600 generations. After 600 or so iterations, the corresponding HV values of the three algorithms no longer increased with the increased number of iterations, and the distribution of the non-dominated solution sets was stable. The HV values are a measure of the quality of the nondominated solution sets of algorithms. As shown in Figure 5, after 600 iterations, the HV value of MOEA/D-AWA was slightly higher than the HV values of MOEA/D and NSGA-II, which indicated that the distribution of the non-dominated solutions of MOEA/D-AWA was slightly more similar to the ideal Pareto front.
In general, although the convergence speed of MOEA/D was faster than that of the two other algorithms, the non-dominated solution distribution of MOEA/D was not uniform. By the above comparison and analysis, the convergence speed of MOEA/D-AWA was relatively fast, and the distribution of the non-dominated solution set of MOEA/D-AWA algorithm was uniform and closer

Analysis of convergence speed of different algorithms
The hyper-volume (HV) is a comprehensive index to measure the quality of multi-objective optimization algorithms. Especially when the ideal Pareto front is unknown, the HV can objectively evaluate the convergence, breadth, and uniformity of the Pareto front [35,36]. For the minimization problem of two-objective optimization, the larger the value of the HV, the better the convergence, breadth, and uniformity of the Pareto front. Additionally, a high HV value means that the non-dominated solution set is a good approximation of the unknown Pareto front [37].
As shown in Figure 5, with increased number of iterations, the convergence curves of HV of three algorithms became more stable. However, there were differences in the convergence speed and HV values of the three algorithms. The curve of MOEA/D started to stabilize after about 150 iterations, that of MOEA/D-AWA gradually stabilized after about 250 iterations, and the convergence rate for NSGA-II was slower than that of the other two. Before 600 iterations, the stability of the distribution of the non-dominated solution set of NSGA-II was relatively poor, and it gradually stabilized after 600 generations. After 600 or so iterations, the corresponding HV values of the three algorithms no longer increased with the increased number of iterations, and the distribution of the non-dominated solution sets was stable. The HV values are a measure of the quality of the non-dominated solution sets of algorithms. As shown in Figure 5, after 600 iterations, the HV value of MOEA/D-AWA was slightly higher than the HV values of MOEA/D and NSGA-II, which indicated that the distribution of the non-dominated solutions of MOEA/D-AWA was slightly more similar to the ideal Pareto front.
In general, although the convergence speed of MOEA/D was faster than that of the two other algorithms, the non-dominated solution distribution of MOEA/D was not uniform. By the above comparison and analysis, the convergence speed of MOEA/D-AWA was relatively fast, and the distribution of the non-dominated solution set of MOEA/D-AWA algorithm was uniform and closer to the theoretical Pareto front. Thus, the MOEA/D-AWA algorithm can be applied to the MORO problem and can achieve better results.
Water 2018, 10, x FOR PEER REVIEW 11 of 15 to the theoretical Pareto front. Thus, the MOEA/D-AWA algorithm can be applied to the MORO problem and can achieve better results.

Rationality of the Operation Results
The MOEA/D-AWA algorithm was used to optimize the MORO problem, providing a nondominated solution set with uniform distribution. The solution set can provide decision support to coordinate the contradictory needs of the two objectives. Figure 6 shows the opposing objectives of water diversion and power generation: as water diversion increased, power generation decreased, and vice versa.
As shown in Figure 6, all the points on the Pareto front were optimal solutions. Selection among these solutions depends on the actual needs of reservoir operation and the choice of the decision makers. It is even possible to add another step in the process to weigh the optimal solutions. As evident in Figure 6, the AB and CD segments changed slowly, and the solutions corresponding to the four endpoints were obtained, namely solution A, solution B, solution C, and solution D. Solution A and solution D correspond to the maximum power generation and the maximum water diversion. As can be seen from Table 3, the water diversion from solution A to solution B underwent a large change, increasing by 3.275%, while the power generation decreased only by 0.414%. The water diversion from solution C to solution D underwent a very large change, with an increase of 1.866%, and power generation decreased only by 0.387%. To assess the comprehensive benefits of the reservoir, the economic benefits of the four solutions were calculated based on the water price of the water source of the Hanjiang to Weihe River Water Diversion Project, at 2 yuan/m 3 [9], and the water price of the hydropower station, at 0.3 yuan/kWh [43]. From solution A to solution B, the average annual economic benefits of the reservoir increased by 3.198%; from solution C to solution D, the average annual economic benefits of the reservoir increased by 1.825%. Although all solutions on the Pareto front were optimal solutions, in terms of the AB and CD gentle segments, solution B was better than solution A and other solutions in the AB segments, and solution C was better than solution D and other solutions in the CD segments in the actual reservoir operation, with greater comprehensive benefits of the reservoir.

Rationality of the Operation Results
The MOEA/D-AWA algorithm was used to optimize the MORO problem, providing a non-dominated solution set with uniform distribution. The solution set can provide decision support to coordinate the contradictory needs of the two objectives. Figure 6 shows the opposing objectives of water diversion and power generation: as water diversion increased, power generation decreased, and vice versa.
As shown in Figure 6, all the points on the Pareto front were optimal solutions. Selection among these solutions depends on the actual needs of reservoir operation and the choice of the decision makers. It is even possible to add another step in the process to weigh the optimal solutions. As evident in Figure 6, the AB and CD segments changed slowly, and the solutions corresponding to the four endpoints were obtained, namely solution A, solution B, solution C, and solution D. Solution A and solution D correspond to the maximum power generation and the maximum water diversion. As can be seen from Table 3, the water diversion from solution A to solution B underwent a large change, increasing by 3.275%, while the power generation decreased only by 0.414%. The water diversion from solution C to solution D underwent a very large change, with an increase of 1.866%, and power generation decreased only by 0.387%. To assess the comprehensive benefits of the reservoir, the economic benefits of the four solutions were calculated based on the water price of the water source of the Hanjiang to Weihe River Water Diversion Project, at 2 yuan/m 3 [9], and the water price of the hydropower station, at 0.3 yuan/kWh [43]. From solution A to solution B, the average annual economic benefits of the reservoir increased by 3.198%; from solution C to solution D, the average annual economic benefits of the reservoir increased by 1.825%. Although all solutions on the Pareto front were optimal solutions, in terms of the AB and CD gentle segments, solution B was better than solution A and other solutions in the AB segments, and solution C was better than solution D and other solutions in the CD segments in the actual reservoir operation, with greater comprehensive benefits of the reservoir.

Conclusions and Discussion
Efficient reservoir operation can make full use of storage capacity of the reservoir and improve the comprehensive benefits of the reservoir. In this study, a model of the MORO problem was proposed to maximize water diversion and power generation objectives, with consideration of various constraints. The MOEA/D-AWA algorithm was applied to solve and optimize the MORO problem and compared with two other algorithms. This model was applied to the Huangjinxia reservoir in Shaanxi, China. The following conclusions can be drawn: (1) The proposed reservoir operation model is effective and reasonable in theory, and can be used to improve the comprehensive benefits of the reservoir. The effectiveness of the proposed model and the introduced algorithm is verified successfully. MOEA/D-AWA has the ability to solve the MORO problem. The results of this work provide a reference for MORO problems, can improve the comprehensive benefits of reservoirs and the utilization efficiency of water resources, thus alleviate the imbalance of water resources in time and space. This is a study for reservoirs with water diversion and power generation tasks. In future work,

Conclusions and Discussion
Efficient reservoir operation can make full use of storage capacity of the reservoir and improve the comprehensive benefits of the reservoir. In this study, a model of the MORO problem was proposed to maximize water diversion and power generation objectives, with consideration of various constraints. The MOEA/D-AWA algorithm was applied to solve and optimize the MORO problem and compared with two other algorithms. This model was applied to the Huangjinxia reservoir in Shaanxi, China. The following conclusions can be drawn: (1) The proposed reservoir operation model is effective and reasonable in theory, and can be used to improve the comprehensive benefits of the reservoir. The effectiveness of the proposed model and the introduced algorithm is verified successfully. MOEA/D-AWA has the ability to solve the MORO problem. The results of this work provide a reference for MORO problems, can improve the comprehensive benefits of reservoirs and the utilization efficiency of water resources, thus alleviate the imbalance of water resources in time and space. This is a study for reservoirs with water diversion and power generation tasks. In future work, it will also be possible to combine with other objectives of reservoirs to build models and involve the multi-objective operation of cascade reservoirs. According to the Sustainable Development Goals (OSS/SDGs) of the UN 2030 agenda, the operation model can be continuously supplemented and improved. The decision support system for reservoir operation can be constantly integrated, and can even include spatial planning tools.