Next Article in Journal
Evaluating the Hydrologic Risk of n-Year Floods According to RCP Scenarios
Previous Article in Journal
Effects of Artificial Reefs on Phytoplankton Community Structure in Baiyangdian Lake, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reservoir Scheduling Using a Multi-Objective Cuckoo Search Algorithm under Climate Change in Jinsha River, China

1
Changjiang Water Resources Commission, Changjiang River Scientific Research Institute, Wuhan 430010, China
2
Hubei Provincial Key Laboratory of Basin Water Resources and Ecological Environment, Changjiang River Scientific Research Institute, Wuhan 430010, China
3
School of Earth and Space Sciences, Peking University, Beijing 100871, China
4
China Institute of Water Resources and Hydropower Research, Beijing 100038, China
*
Author to whom correspondence should be addressed.
Water 2021, 13(13), 1803; https://doi.org/10.3390/w13131803
Submission received: 30 May 2021 / Revised: 27 June 2021 / Accepted: 28 June 2021 / Published: 29 June 2021

Abstract

:
Changes in rainfall and streamflow due to climate change have an adverse impact on hydropower generation reliability and scheduling of cascade hydropower stations. To estimate the impact of climate change on hydropower, a combination of climate, hydrological, and hydropower scheduling models is needed. Here, we take the Jinsha River as an example to estimate the impact of climate change on total power generation of the cascade hydropower stations and residual load variance of the power grid. These two goals are solved by applying an improved multi-objective cuckoo search algorithm, and a variety of strategies for the optimal dispatch of hydropower stations are adopted to improve the efficiency of the algorithm. Using streamflow prediction results of CMIP5 climate data, in conjunction with the Xinanjiang model, the estimated results for the next 30 years were obtained. The results indicated that the negative correlation between total power generation and residual load variance under the RCP 2.6 scenario was weaker than that under the RCP 8.5. Moreover, the average power generation and the average residual load variance in RCP 2.6 was significantly larger than that in RCP 8.5. Thus, reducing carbon emissions is not only beneficial to ecological sustainability, but also has a positive impact on hydropower generation. Our approaches are also applicable for cascade reservoirs in other river catchments worldwide to estimate impact of climate change on hydropower development.

1. Introduction

Hydropower is one of the most important renewable, and environmentally sustainable energy resources. It can store energy at low costs, maintains lower operating and maintenance costs, and can operate with great stability [1,2]; thus, when integrated with intermittent energy sources, such as solar and wind, can make a significant contribution to the consistency of the energy grid [3,4]. As shown in the 2019 BP Statistical Review of World Energy, global primary energy consumption has grown rapidly in recent years, with hydroelectric generation increasing by 3.1%; however, carbon emissions grew by 2.0% in 2018, the fastest over the previous seven years. The emissions of greenhouse gases from anthropogenic activity are the primary mechanism driving the changing global climate, imposing severe stress on the global ecosystem with potentially immeasurable consequences. Climate change is having a direct effect on meteorological factors such as temperature, rainfall, and evaporation, and indirectly affects other factors such as soil moisture content and runoff, resulting in the spatiotemporal redistribution of water resources [5]. It is projected that annual mean streamflow will increase in the high-latitudes and wet tropical regions, decrease in most of the dry tropics, and global hydrological extremes, such as rainstorms, floods, and droughts, will continue to increase significantly, preventing investors from accurately forecasting future energy production [6,7]. Variability in streamflow resulting from climate change is likely to have a significant impact on hydropower generation [8]. The impact factors of climate change on river systems, such as temperature fluctuations, rainfall patterns, floods, and droughts [9,10], have a direct or indirect influence on hydropower generation [11]. Accordingly, hydropower planning must be careful and take into account the potential changes in river streamflow as a result of a changing climate [12], and such research is critical for the development of mitigation plans.
Increasing anthropogenic activities over the past century have resulted in increased atmospheric concentrations of greenhouse gases CO2, CH4, and N2O, the primary gaseous drivers of climate change [13]. Increases in global temperature are self-accelerating, as the melting snow in polar regions and high-altitude mountains decreases surface albedo as low-elevation snow cover, glaciers, and permafrost are in decline [14]. The melting of glaciers and snow in the polar regions, and the resulting rise in sea level, will increase oceanic area, which along with higher temperatures, will increase evaporation and correlated large-scale changes in the global water cycle. Water vapor from the ocean is converted into precipitation, and when rainfall increases or decreases in wet or dry areas, respectively, the probability of natural disasters, such as floods and droughts, increases significantly. Changes in glacial and snow-covered areas of mountains have also altered river runoff [15]; however, streamflow changes by region have neither consistently increased nor decreased.
With respect to the indirect impacts of climate change on hydropower, higher surface air temperatures will change the temperature and stratification of a reservoir [16], then increase evaporation over the reservoir, reducing the amount of available water [17]. Second, stronger rainfall will bring increased pollutants into rivers, potentially exacerbating reservoir contamination, and the resulting wear and corrosion of turbines by pollutants will reduce power generation efficiency [18]. Third, as reservoirs require constant water supply regulation, rising temperatures will increase agricultural and domestic water demand, reducing reservoir storage levels, resulting in a decrease in the water head and generation flow.
The research here investigated the impacts of climate change on hydropower production and proposed a method to quantitatively analyze the economic effects on power generation of the cascade hydropower stations along the Jinsha River, China. The remainder of this paper is organized as follows: Section 2.1 describes the available methods for quantitative analysis; Section 2.2 presents the improved multi-objective cuckoo search algorithm (MoCS), a tool for finding optimal solutions; Section 2.3 and Section 2.4 describes the multi-objective hydropower generation model (MLTHG) and the application of MoCS to deal with its limitations; Section 3 describes the impact of future climate changes on reservoir scheduling in the Jinsha River; and Section 4 provides a summary conclusion.

2. Methods

2.1. Estimating Climate Change Impacts on Hydropower

New tools, such as the latest generation of general circulation models (GCMs) and downscaling methods, have advanced the ability to study the potential impacts of climate change on hydropower generation [19]. It is possible to identify the impact of climate change on water resource by simulating climate change scenarios [20]. The most common method is to run a calibrated baseline hydrological model under changing hydro-climatic conditions to estimate streamflow [21] and using this value as the input of the power system model [22,23,24].
Climate change is a complex process characterized by massive uncertainty; accordingly, future research scenarios based on a series of scientific assumptions are generally adopted, using GCMs to simulate the future climate conditions. In combination with downscaling technology, GCM climate data can be converted to that required by hydrological models, such as empirical statistical, conceptual hydrological, and distributed hydrological models, and ultimately be used to assess the impact of climate change on water resources. The conceptual hydrological model is a physical approach to describing the relationship between hydrological factors by generalizing complex phenomena [25,26]. The distributed hydrological model, such as TOP model, SWAT, Mike SHE, and VIC, divides the watershed into multiple sub-simulation units, and performs separate hydrological simulations to better reflect the characteristics of medium-and large-scale watersheds [27,28].
After completing a rainfall runoff simulation, a power system model can be applied to estimate the impacts on energy production. For example, Queiroz et al. introduced spatiotemporal information of future climate change into the operation plan of Brazil’s hydropower system to evaluate its impact on the revenue of hydropower stations [29]. By applying GCM projections, the trend of hydropower demand over the next 30 years can be obtained [30]. Other objectives can also be estimated through this method. For example, Zhai et al. predicted the trend of inflow of the Jinsha river through two climate models and carried out a qualitative analysis of the possible strategies of hydropower stations for flood control purposes in the future [31]. Zhang et al. studied the derivation and the adaption of operating rules for the irrigation reservoir under climate change [32].
To best estimate the impact of climate change on hydropower, a combination of climate, hydrological, and power system models is needed. The accuracy of the hydrological model, and the efficiency of the power system model, will determine the accuracy of the climate change impact prediction. Uncertainty of the results mainly comes from three aspects: Observational Uncertainty, Model Structural Uncertainty, and Parametric Uncertainty. In order to evaluate the uncertainty, probabilistic analysis methods, such as Expectancy Method, Utility Function Method, and Model Analysis, are mainly used. Given the uncertain factors a random change within a certain range, analyze and determine the probability distribution of this change, so as to calculate its expected value and standard deviation. The flowchart of the method used in this paper is shown in Figure 1. Details will be explained in the following sections.

2.2. Improved Multi-Objective Cuckoo Search Algorithm

2.2.1. NGSA-II

In general, real-world engineering problems must consider multiple goals, and tradeoffs are required for goal optimization; therefore, one must consider the relationship between conflicting targets to obtain the most reasonable multi-objective solution that maximizes total benefits. The concept of “dominance” is commonly used in multi-objective optimization algorithms. Pareto optimal solutions are feasible, and not dominated by other arbitrary solutions in the set. The multidimensional space formed by the target values of all Pareto optimal solutions is the Pareto optimal front; at present, the main methods for obtaining the frontal space are mathematical programming and intelligent algorithms.
Elitist nondominated sorting genetic algorithm v. II (NSGA-II) is a popular genetic algorithm-based multi-objective optimization method. It uses an elite strategy and maintains population diversity through the operator of crowding comparison. The basic process of NSGA-II is:
(1)
Randomly generate an initial population of size N. Then, sort the initial population using the non-dominated method.
(2)
Generate new subpopulations through selection, crossover, and mutation operations of the evolutionary algorithm.
(3)
Merge the parent and child populations and select the most fit individuals to form the next-generation population based on a dominance relationship and crowding degree.
(4)
Repeat the above processes until the termination condition is met.

2.2.2. Improved Cuckoo Search

The cuckoo search (CS) algorithm is a meta-heuristic method developed by Yang and Deb (2009). It is inspired by the obligate brood parasitism of some cuckoo species [33].
The most important feature of the CS algorithm is the use of Lévy flight characteristics, where a number of studies have shown that the behavior of many insects and larger animals exhibited the typical characteristics of a Lévy flight, shown by a power-law behavior. Inspired by this action, a new solution was generated as follows:
x i t + 1 = x i t + α L e v y ( λ )
Research has shown that the strength of the local convergence ability of basic CS insufficiently and several improvements have been proposed to improve its performance:
(1)
Dynamic parameter adjustment strategy
In the basic CS algorithm, the abandon probability (pa) is a fixed value that affects the convergence speed of the algorithm. The improved solution turns pa into a dynamic adjustment using Equation (2):
p a = p a s + ( p a e p a s ) × C E N E
where p a s and p a e are the start and end of p a , respectively, and CE and NE are the current and maximum evaluation times, respectively.
(2)
Differential strategy for Lévy flight
Originally, Lévy flights were described according to Equation (3), although the details of its implementation are not provided. Different interpretations of this formula have been suggested by researchers, and in the present study, we implemented the equation using the differential strategy as follows:
x i d = x i d + ( s l × L e v y ( u , c ) × ( x j d x i d ) )
where xi and xj are randomly chosen nests, sl is the step size, and Levy(u,c) is the sample value generated by the Lévy distribution. The probability density function of Levy(u,c) is calculated according to Equation (4):
f ( x ) = c 2 π e c 2 ( x u ) ( x u ) 3 / 2
(3)
Revised solution
The basic CS generates a new solution and compares it with randomly chosen ones from the nests. If the randomly chosen solution is inferior, it is replaced; however, the two solutions are irrelevant, and the conversion is inadequate, and the source solution of the new solution should be replaced accordingly. The pseudo-code of the improved cuckoo search (ICS) Algorithm 1 is as follows:
Algorithm 1. Improved cuckoo search
Objective function f ( x ) , x = ( x 1 , , x d ) T
Initialize default parameters
Generate initial population of n host nests x i ( i = 1 , 2 , , n )
 While (t < MaxEvaluation) or (stop criterion)
 Select two solution x i , x j from host nests randomly
 For d=1,…, D do
   x i d = x i d + ( s l × L e v y ( u , c ) × ( x j d x i d ) )
 End for
 If f ( x i ) < f ( x i )
  Replace xi by the new solution x i
 End if
 If r a n d ( 0 , 1 ) < p a
  Init the worst nest x w o r s t
 End if
End while

2.2.3. Multi-Objective Cuckoo Search

The multi-objective cuckoo search (MoCS) in the present study was based on the NSGA-II algorithm and the improved cuckoo search algorithm described above. The MoCS procedure is explained as follows (Figure 2):
(1)
Generate a random initial population, and classify the individuals using the non-dominated sorting method. A new population can then be generated using the ICS algorithm.
(2)
Merge the parent and child populations, employ a fast, non-dominated sorting and crowding degree calculation in the mixed population, and select the most fit individuals to form the next generation.
(3)
Repeat the processes above until the termination conditions are met.

2.3. Gradient Multi-Objective Cuckoo Search for Reservoir Scheduling

A reasonable scheduling plan can achieve efficient trade-offs for various conflicting objectives [34]. The objectives and constraints of the multi-objective long-term hydropower generation (MLTHG) are described in this section, as is the application of the MoCS for MLTHG.

2.3.1. Power Generation Objective

One of the optimization goals of the MLTHG is to maximize total power generation of cascade hydropower stations, using the following formula Equation (5):
f 1 = m a x t = 1 T i = 1 N k i Q i , t H i , t Δ t
where f1 (kW∙h) is the total energy production of the cascade hydropower system, T is the period count, N is the number of reservoirs, ki is the output coefficient of reservoir i, Qi,t (m³/s) is the generation of outflow through hydropower units of reservoir i at time t; Hi,t (m) is the net water head of hydropower reservoir i at time t, and ∆t is the length of the time interval.

2.3.2. Residual Load Variance Objective

An important function of a hydropower station is to optimally respond to the power grid, thus minimizing the residual load; however, this objective is not easily resolved, and researchers have carried out a variety of targeted improvements towards this end. In the present study, the objective was to minimize the residual load variance of the power grid, according to the following formula Equation (6):
{ f 2 = m i n ( t T ( P t P ¯ ) 2 ) P t = i = 1 N k i Q i , t H i , t
where f2 (kW) is the residual load variance of the power grid, Pt (kW) is the total power output of the cascade hydropower stations over period t, and P ¯ (kW) is the average power output of all periods. Notably, it is necessary to maintain a low water level during the flood season due to flood control requirements; thus, this season was not considered in the formula.

2.3.3. Constraints

(1)
Hydraulic connection (Equation (7)):
{ I i , t = O i 1 , t + R i , t O i 1 , t = Q i 1 , t + S i 1 , t
where I i , t (m³/s) is the inflow of reservoir i at time t;  O i 1 , t (m³/s) is the outflow of reservoir i − 1 at time t, R i , t (m³/s) is the interval inflow between reservoir i − 1 and i;  S i 1 , t (m³/s) is the spillage of the upstream reservoir i − 1, and Q i 1 , t (m³/s) is the generation flow of the upstream reservoir i − 1.
(2)
Water-balance constraint (Equation (8)):
V i , t = V i , t 1 + ( I i , t O i , t ) Δ t
where Vi,t (m³) is the storage of reservoir i in time t.
(3)
Water-level constraints (Equations (9) and (10)):
Z i , t m i n Z i , t Z i , t m a x
| Z i , t Z i , t 1 | Z i , t s t e p
where Z i , t m i n (m) and Z i , t m a x (m) are the lower and upper bounds of the water level, respectively; and Z i s t e p (m) is the limitation of the water-level variation.
(4)
Outflow constraint (Equation (11)):
O i , t m i n O i , t O i , t m a x
where O i , t m i n (m³/s) and O i , t m a x (m³/s) are the minimum and maximum outflows of reservoir i, respectively.
(5)
Output constraint (Equation (12)):
N i , t m i n N i , t N i , t m a x
where N i , t m i n (kW∙h) and N i , t m a x (kW∙h) are the minimum and maximum power output of reservoir i, respectively.
(6)
Boundary condition (Equation (13)):
Z i , 0 = Z i , s t a r t , Z i , T = Z i , e n d
where Z i , s t a r t (m) and Z i , e n d (m) are the initial and terminal water levels of reservoir i, respectively.

2.3.4. Solution Encoding and Initialization

When applying the MoCS to address multi-objective power generation scheduling, the water levels or outflows of the hydropower stations are often chosen as the decision variable. In the present study, the water level was employed, and encoded as (Equation (14)):
Z = [ Z 1 Z 2 Z N ] = [ z 1 1 , z 1 d , z 1 T z 2 1 , z 2 d , z 2 T z N 1 , z N d , z N T ]
where z i d (m) is the water level at hydropower station i and time d. The initial value of z i d is randomly generated as:
z i d = z i m i n + R A N D ( 0 1 ) · ( z i m a x z i m i n )
where R A N D ( 0 1 ) is a randomly generated value between 0 and 1; and z i m a x (m) and z i m i n (m) are the upper and lower boundaries of reservoir i, respectively. When a new individual is randomly generated, there is a high probability that it will not meet the constraints, and it will be replaced with new individuals until the new individuals meet the constraints.
In the evolution stage, the individual does not meet the constraints and is corrected by adjusting the water level to the feasible range as follows:
z i = { z i m i n ,   i f   z i < z i m i n z i m a x ,   i f   z i > z i m a x z i ,   o t h e r w i s e
where z i m i n and z i m a x are determined by the water level, outflow, and power output constraints, respectively according to Equations (17)–(19):
z i m a x = min ( z i 1 + z m a x , Z ( z i 1 , Q i m i n ) ,   Z ( z i 1 , Q i P   m i n ) , z m a x )
z i + 1 m i n = max ( z i z , Z ( z i , Q i m a x ) , z m i n )
Z ( z i , Q i ) = Z ( V ( z i ) + ( I i Q i ) Δ t )
where Q t m i n (m³/s) and Q t m a x (m³/s) are the minimum and maximum outflows of the reservoir; Q t P m i n (m³/s) is the minimum outflow for a guaranteed output; z (m) is the limitation of the water-level change, Z ( V ) and V ( Z ) are the relationship between the water level and the storage.
A combination of forward and reverse corrections was used when applying the correction to an unfeasible range. First, the water level was adjusted from time period 1 to T. If the correction failed midway, attempts were made to correct the water level from period T to the breakpoint. If the correction was unsuccessful, the fitness of the individual was marked as 0.

2.3.5. Gradient Search Strategy

Following the gradient search strategy of the MoCS algorithm, the solution was adjusted by a small gradient, as shown in Figure 3. If the water level of a reservoir changed, both its power generation capacity and its downstream reservoirs changed.
{ Z i , t a d j u s t = Z i , t + Δ l Δ Q i , t = ( V i ( Z i , t + Δ l ) V i ( Z i , t ) ) / Δ T t Δ H i , t = Δ l / 2 Z d i ( Q i , t + Δ Q i , t ) + Z d i ( Q i , t ) Δ P i , t = k i ( Q i , t + Δ Q i , t ) ( H i , t + Δ H i , t ) k i Q i , t H i , t
Accordingly, the partial derivative of power generation E (kW∙h), with respect to gradient Δl is (Equation (21)):
E l = Δ T t P i , t l + Δ T t + 1 P i , t + 1 l + j = i + 1 N ( Δ T t P j , t l + Δ T t + 1 P j , t + 1 l )
where if E l > 0 , the water level is adjusted by Δl; otherwise, the adjustment is rejected. Detail derivation process can be found in [35]. A flowchart of the gradient-based search strategy for the cascade system is shown in Figure 4.

2.3.6. Single Entry External Archive

Elite individuals were inserted into the concept of an external archive that was introduced in the MoCS algorithm. When the number of individuals in the archive exceeded the limit, the non-optimal solution was removed based on the non-dominated sorting and crowding degree calculation.
Notably, if all of the elite individuals were inserted into the external archive at the same time, some may be concentrated. The crowding degrees of these individuals will be small, and there is a high probability that they will be removed during selection, which is not conducive to the generation of the Pareto front. Therefore, it was necessary to insert only one elite individual at a time, with the subsequent individual inserted only after the selection was completed. The Pareto optimal front will be more uniform following this strategy shown in Figure 5.

2.3.7. Self-Tuning Divergent Operator Strategy

A divergent operator strategy was introduced to increase population diversity. During new individual generation, a certain number from the population were randomly selected and underwent single-objective evolution separately. Simultaneously, the remaining individuals executed the original evolution selection strategy.
The number of individuals performing single-target evolution (n) was determined by the number of targets (Figure 6). If the objective maximum changed, the amount of the corresponding single-target individuals increased, and if no change was observed, this number decreased. The purpose of this strategy was to expand the algorithm’s ability to search for edges.

2.3.8. Procedures of Solving MLTHG with Gradient Based MoCS

Accordingly, the procedures of solving MLTHG with Gradient based MoCS (GMoCS) were as follows (Figure 7):
(1)
Set the initial conditions of the MLTHG model, including incoming water, water level of each hydropower station in the initial time period and the end time period.
(2)
Set up the constraints of the MLTHG model, including the maximum outflow, minimum outflow, maximum water level, minimum water level, and water level variation in each period of each hydropower station. The setting of water level constraints needs to consider the flood control requirements in the flood season.
(3)
Set the parameters of the GMoCS algorithm.
(4)
Generate a random individual according to Equations (14) and (15), calculate its objective values and check whether it meets all constraints. If the constraints are met, added it to the initial population. Generate other new individuals using the same method until the population size is reached.
(5)
Generate new individual by Lévy flight.
(6)
Calculate objective values of the new individual and check whether it meets all constraints. If the new individual does not meet the constraints, try to adjust the water level to make it meet the constraints.
(7)
Mix the old and new populations. Perform non-dominated sorting and crowding calculations on the mixed population. Select better individuals to form the next generation population.
(8)
Select the individual with the largest power generation and perform gradient search for the power generation objective. Then, select the individual with the smallest residual load variance, and perform gradient search for the residual load variance objective.
(9)
Select non-dominated individuals from the next generation population. Insert the non-dominated individuals into the external archives one by one, while replacing the worst individual with a better random individual with a certain probability.
(10)
Repeat steps 3–9 until the termination condition is met.
(11)
Select non-dominated individuals from the external archive set as the Pareto optimal frontier. Export time series data of these individuals, including water level, outflow, power output, etc.

2.4. MLTHG in Jinsha River

2.4.1. The Projection of Streamflow in the Context of Climate Change

The historical observed daily precipitation and temperature (maximum and minimum) data are collected from China’s Ground Precipitation 0.5° × 0.5° Gridded Dataset (V2.0) and China’s Ground Temperature 0.5° × 0.5° Gridded Dataset (V2.0) (http://data.cma.cn/ accessed on 25 May 2021), respectively. The dataset covers 2474 national meteorological stations over the Chinese mainland. It consists of a series of daily observed data from 1961 to the present. We use data from 208 grid boxes associate with the Jinsha River, shown in Figure 8.
The future climate change scenarios are obtained from the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP) [36]. Continuous daily precipitation and temperature data on a horizontal grid with 0.5° × 0.5° resolution for the period 1960–2099 were offered by five GCMs, including GFDL-ESM2M, HADGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, and NORESM1-M. Output of the climate model was spatially interpolated using a bilinear interpolation method, and the data has already been bias-corrected using a trend-preserving bias correction method [37]. The bias correction method modifies the monthly mean and daily variability of the simulated data to match the observations. Change of average daily precipitation and temperature (maximum and minimum) between each GCM and observations are shown in Figure 9, Figure 10 and Figure 11.
The main hydrological model used for streamflow forecasting is the Xinanjiang (XAJ) model. The daily precipitation of each sub-catchment is evaluated following:
P = i = 1 N P i S i i = 1 N S i
where P is the average precipitation of the sub-catchment, Pi is the precipitation of the grid box i, and Si is the area of grid box i. N is the number of grid boxes in the sub-catchment.
The daily potential evapotranspiration of each grid is calculated using Hargreaves method [38]:
E T 0 = H A · R e ( T + 17.8 ) · T H E
where HA and HE are the parameters with standard values of 0.0023 and 0.5, respectively. Re is the extraterrestrial radiation. T is the mean temperature (T = (Tmax + Tmin)/2) and ΔT is the air temperature range (ΔT = Tmax – Tmin). The daily evapotranspiration of each sub-catchment is calculated in the same way with precipitation.
The parameters of XAJ model are calibrated using the observed hydrological data and the calibrate method is based on ICS. The outflow of each sub-catchment is routed using the Muskingum method. The precipitation and temperature data from 1961 to 2050 given by five GCMs in three different scenarios (RCP2.6, RCP4.5, RCP8.5) was input to the parameterized XAJ model, and the daily streamflow in the context of climate change can be simulated. Detailed prediction methods can be found in [39].

2.4.2. Modeling of MLTHG in Jinsha River

The Jinsha River is in the upper reaches of the Yangtze River. It is born at the junction of Qinghai and Sichuan Province, and ends at the mainstream of the Yangtze River in Yibin. There is strong hydrological potential along the river, particularly in the lower reaches where elevation decreases dramatically. Four large cascade hydropower station—Wudongde, Baihetan, Xiluodu, and Xiangjiaba—are located along the lower reaches of the river (Figure 12). The main primary station parameters are listed in Table 1, and the topology structure is shown in Figure 13.
There are two hydrological stations in the basin—Huatan and Pinshan, and the streamflow of the four hydropower stations was projected using the Xinanjiang model. The inflow of Wudongde was generalized using the streamflow of the Huatan station, and the interval inflow between Baihetan and Xiluodu was generalized using the streamflow of the Pingshan station according to Equations (22) and (23):
I 1 = Q Huatan
Δ I 3 = Q Pingshan Q Huatan
The power generation at time t can be obtained by Equation (24):
E i , t = k i Q i , t H i , t Δ t
where Ei,t (kW∙h) is the hydropower generation of reservoir i at time t; Qi,t (m³/s) is the generation flow through the hydropower units of reservoir i at time t; ki is the comprehensive benefit coefficient of reservoir i, which reflects the hydro-generating unit efficiency; Δt is the length of time interval, and Hi,t (m) is the net water head of reservoir i at time t, defined by Equation (25):
H i , t = Z i , t 1 + Z i , t 2 Z d i , t Δ H i
where Zi,t−1 (m) and Zi,t (m) are the beginning and end water levels of reservoir i at time t, respectively; Zdi,t (m) is the water level under the dam; and ΔHi (m) is the water head loss of reservoir i.
Overall, the objective of the MLTHG scheduling model was to find an optimal set of water releases or storage volumes to maximize the MLTHG of cascade hydropower stations, which can be described mathematically according to Equation (26):
E = m a x t = 1 T i = 1 N E i , t
where E (kW∙h) is the total hydropower generation of the cascade reservoirs, T is the number of time periods, and N is the number of reservoirs.

3. Results and Discussion

3.1. Performance of GMoCS

To evaluate the performance of the GMoCS, the optimal results of the cascade hydropower stations were compared with those obtained using MoCS. The parameters of GMoCS were: n s = 40 ,   p a s = 0.3 ,   p a e = 0.1 ,   s l = 0.01 ,   u = 0 ,   c = 1.5 , and the maximum evaluation time was set to 12000. The initial and end water levels of the reservoirs were set to their normal levels (Table 1), and each month was divided into three periods of ~10 days each, resulting in 36 time intervals per year. Each method was run 100 times, independently.
Three typical streamflows of Pingshan station (Figure 14) chosen by the annual runoff frequency curve —dry (frequency 75%), normal (50%), and wet (25%) years—are used as the input of MLTHG. Water levels of the four reservoirs at different time intervals (excluding the 36th) were chosen as the optimization parameters, making the problem dimension = 140. Equations (5) and (6) are the fitness functions of MLTHG, and the greater the power generation and the smaller the residual load variance, the better the solution. The simulation results are shown in Figure 15.
The Pareto optimal frontier obtained by GMoCS was superior to the results of MoCS (Figure 15). MoCS fell into a local optimum when solving the MLTHG; whereas the results obtained by GMoCS indicated that the strategies proposed in the present study effectively improved the algorithm’s search ability. It can be seen that total power generation of the cascade hydropower stations and residual load variance of the power grid were inversely correlated. These two objectives were conflicting and mutual. The scheduling processes of the 2nd and 49th solutions of wet year are shown in Figure 16 as an example. The difference between the two solutions is mainly in the period before the flood season (June, July). Due to the large inflow during the flood season, the power generation will increase significantly. The 2nd solution increases the power generation capacity for a period of time before the flood season by lowering the water level in advance. This makes the growth of power generation before the flood season relatively smooth, which is conducive to reducing the residual load variance of the power grid. The 49th solution maintains high water level for the period before flood season, resulting in an increase in the total power generation. While the residual load variance of the power grid is larger than that in 2nd solution.
In this section, we focus on the trade-off between the two objectives of power generation and residual load variance. Research done by Geressu et al. [34] shows that when the financial benefit of the reservoirs is increased, the downstream release will be adversely affected. This is similar to our research. Only a compromise plan can be chosen if we consider two conflicting goals. Moreover, as explained by Zhai et al. [31], during the flood season the cascade hydropower stations are operated below flood-limiting water level in order to reserve adequate reservoir storage for flood prevention. We also considered the same restriction, the change of the water level is mainly in the non-flood season as shown in Figure 16.

3.2. Climate Change on Multi-Objective Scheduling of Cascade Hydropower Stations

The impact of climate change on the total power generation of the cascade hydropower stations and residual load variance of the power grid over the next 30 years was estimated. First, streamflow prediction data of Pingshan station for the years 2021–2050 were prepared (Figure 17), and prediction results were obtained by applying the CMIP5 climate data and the Xinanjiang model (detailed prediction methods can be found in [39]). Then, the multi-objective scheduling results under three climate change scenarios—RCP2.6, RCP4.5, and RCP8.5—using five climate models—GFDL-ESM2M, HADGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, and NORESM1-M—were obtained with 30 years of daily streamflow data as the input. Accordingly, 450 solution sets were obtained. Owing to the inaccuracy of the predictions of various climate models, the Pareto frontiers also showed large levels of uncertainty.

3.2.1. Power Generation Objective

As explained in Section 2.3.7, the boundary of the Pareto frontier is close to that obtained through single-objective optimization. First, the optimal solution for power generation was obtained with the operator. Figure 18 displays average total power generation results per 10 years for the next 30 years under the three climate change scenarios, and Figure 18a shows the average generation results for scenario RCP2.6. The average generation results of the five GCMs were 2266, 2238, and 2188 (108 kWh), respectively, while the uncertainty ranges were similar to each other: (−4.4%, 4.0%), (−5.2%, 3.9%), and (−4.3%, 5.4%). Figure 18b,c show the results under scenarios RCP4.5 and RCP8.5, which both displayed much larger uncertainty ratios than RCP2.6. Moreover, the average power generation in RCP2.6 was significantly larger than that in RCP8.5, and when compared with the average power generation in RCP4.5, RCP2.6 was larger for the first 20 years, but smaller following decade.

3.2.2. Residual Load Variance Objective

The optimal solution for the residual load variance was obtained by using the respective operator. Figure 19 displays the predicted results for the next 30 years under each scenario. Figure 19a shows the results under scenario RCP2.6. The average residual load variance results for the five GCMs were 910, 916, and 867 (108 kWh), for 2021–2030, 2031–2040, 2041–2050, respectively; while the uncertainty ranges were (13.6%, 6.0%), (6.1, 6.5%) and (−9.7%, 17.5%). Figure 19b,c show the results under scenarios RCP4.5 and RCP8.5, respectively, with relatively larger uncertainty ranges. Similar to the characteristics of power generation in the preceding section, the average residual load variance in RCP2.6 was larger than RCP8.5 at each interval, and larger than RCP4.5 for the first 20 years, but smaller following decade.

3.2.3. Combination of the Two Objectives

The impact of climate change on the combination of the two objectives was evaluated using the contradiction between the two goals (Equation (29)):
C = f 2 i , m a x f 2 i , m i n f 2 m a x f 2 m i n · ( f 1 i , m a x f 1 i , m i n f 1 m a x f 1 m i n ) 1
where C is the variation ratio; f 1 m a x and f 1 m i n are the maximum and minimum values of the power generation objective in the solution set, respectively; f 2 m a x and f 2 m i n are the maximum and minimum values of the residual load variance objective in the solution set, respectively; f 1 i , m a x ( m i n ) is the maximum (minimum) value of the solution i. When C was larger, the reduced power generation due to the increase in the residual load variance was larger, implying a clear contradiction between the two goals; however, when C was smaller, the reduced power generation due to the increase in the residual load variance was decreased, implying a more subtle contradiction between the two objectives.
Figure 20 shows the predicted variation ratios of the total power generation and residual load variance over the next 30 years, at decadal intervals. The results indicated that: under scenario RCP2.6 scenario indicated the average ratio of change between 2021–2030 and 2041–2050 was similar; under scenario RCP4.5, the average variation ratio in 2021–2030 was significantly higher than that in 2031–2040 or 2041–2050; however, the variation ratios in 2031–2040 and 2041–2050 were similar to those under RCP2.6; and the average variation ratio under the RCP8.5 scenario was even higher than that of RCP2.6.
In this section, we focus on the impact of climate change on two objectives: total power generation and residual load variance. The key factor associated with climate change and hydropower generation is streamflow, which is the same with the research done by Zhang [32] and Zhai [31]. However, due to the different climate scenarios selected, there is a difference in the predicted results. Zhang used a climate scenario with reduced rainfall and increased evaporation and concluded that agricultural profits are declining. Zhai adopted a climate scenario with increased discharge and concluded that extreme hydrological events are increasing. Trend of the objectives that are greatly affected by streamflow is generally depended on the trend of streamflow. This is also reflected in our article.
As can be seen in this paper, the results of RCP2.6 were more positive than those of RCP8.5. This is the same with the research done by Nam [20]. Nam carried out the assessment of vulnerable seasons for paddy irrigation under two climate change scenarios: RCP4.5 and RCP8.5. Results show that the total duration of the vulnerable irrigation seasons is longer in RCP8.5 than that in RCP4.5. Therefore, it is necessary to reduce carbon emissions to mitigate the adverse effects of climate change.

4. Conclusions

In the present research, a method to estimate the impact of climate change on hydropower was discussed, and a comprehensive impact of climate change assessment was made on the power generation of cascade hydropower stations in the lower reaches of the Jinsha River and residual load variance of the power grid. Climate change is affecting the amount of available water for power generation in a variety of ways, thereby altering the cascade of hydropower stations. The most common method to estimate the impacts of climate change on hydropower generation is to run a hydrological model under climate change conditions, obtain the changes in streamflow, and then assess the correlated impact on hydropower. An improved multi-objective cuckoo search algorithm has been proposed, and a variety of strategies for the optimal dispatch of hydropower stations, such as the gradient search strategy, single-entry external archive, and self-turning divergent operator strategy, were adopted to improve the efficiency of the algorithm.
In the case studies, the improved and unimproved algorithms were compared. The optimal results of the three typical years obtained by MoCS and GMoCS showed that the Pareto optimal frontier obtained by GMoCS was vastly superior to the results obtained by MoCS. The results also showed that the trends of total power generation of cascade hydropower stations and residual load variance of power grid were inversely correlated. After simulating the average variation ratios of five climate models under the three climate change scenarios, it was revealed that the strength of this relationship under the RCP 2.6 scenario was weaker than that under the RCP 8.5.
Overall, the impact of climate change on hydropower is reflected in the total power generation, as well as other indicators, including residual load variance. Although the results obtained through different climate models did not produce a unified conclusion, the results of RCP2.6 were more positive than those of RCP8.5, indicating that reducing carbon emissions is not only beneficial to ecological sustainability, but also has a positive impact on hydropower generation. Our approaches are also applicable for cascade reservoirs in other river catchments worldwide to estimate impact of climate change on hydropower development.

Author Contributions

Y.F. and C.W. conceived and designed the study; Y.F. and C.W. proposed the MoCS algorithm and the multi-objective hydropower generation model; Z.Y. and Y.W. wrote the predict model and performed the simulations. J.X. and Y.H. advised on the work and performed the final checks; All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No. U2040212, No. 52009005) and the National Public Research Institutes for Basic R and D Operating Expenses Special Project (no. CKSF2019212/SZ).

Data Availability Statement

Not applicable.

Acknowledgments

We acknowledge the entire development team, whose help throughout this research was indispensable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mendes, C.A.B.; Beluco, A.; Canales, F.A. Some important uncertainties related to climate change in projections for the Brazilian hydropower expansion in the Amazon. Energy 2017, 141, 123–138. [Google Scholar] [CrossRef]
  2. Fan, J.-L.; Hu, J.-W.; Zhang, X.; Kong, L.-S.; Li, F.; Mi, Z. Impacts of climate change on hydropower generation in China. Math. Comput. Simul. 2020, 167, 4–18. [Google Scholar] [CrossRef]
  3. Gaudard, L.; Romerio, F.; Valle, F.; Gorret, R.; Maran, S.; Ravazzani, G.; Stoffel, M.; Volonterio, M. Climate change impacts on hydropower in the Swiss and Italian Alps. Sci. Total Environ. 2013, 493. [Google Scholar] [CrossRef]
  4. Chang, J.; Wang, X.; Li, Y.; Wang, Y.; Zhang, H. Hydropower plant operation rules optimization response to climate change. Energy 2018, 160, 886–897. [Google Scholar] [CrossRef]
  5. Zhong, W.; Guo, J.; Chen, L.; Zhou, J.; Zhang, J.; Wang, D. Future hydropower generation prediction of large-scale reservoirs in the upper Yangtze River basin under climate change. J. Hydrol. 2020, 588, 125013. [Google Scholar] [CrossRef]
  6. Qin, P.; Xu, H.; Liu, M.; Du, L.; Xiao, C.; Liu, L.; Tarroja, B. Climate change impacts on Three Gorges Reservoir impoundment and hydropower generation. J. Hydrol. 2020, 580, 123922. [Google Scholar] [CrossRef]
  7. Kim, K.; Jeong, H.; Ha, S.; Bang, S.; Bae, D.-H.; Kim, H. Investment timing decisions in hydropower adaptation projects using climate scenarios: A case study of South Korea. J. Clean. Prod. 2017, 142, 1827–1836. [Google Scholar] [CrossRef]
  8. Boehlert, B.; Strzepek, K.M.; Gebretsadik, Y.; Swanson, R.; McCluskey, A.; Neumann, J.E.; McFarland, J.; Martinich, J. Climate change impacts and greenhouse gas mitigation effects on U.S. hydropower generation. Appl. Energy 2016, 183, 1511–1519. [Google Scholar] [CrossRef] [Green Version]
  9. Huang, S.; Krysanova, V.; Hattermann, F. Projections of climate change impacts on floods and droughts in Germany using an ensemble of climate change scenarios. Reg. Environ. Chang. 2014, 15. [Google Scholar] [CrossRef]
  10. Mezosi, G.; Bata, T.; Meyer, B.; Blanka, V.; Ladányi, Z. Climate Change Impacts on Environmental Hazards on the Great Hungarian Plain, Carpathian Basin. Int. J. Disaster Risk Sci. 2014, 5, 136–146. [Google Scholar] [CrossRef] [Green Version]
  11. Hasan, M.M.; Wyseure, G. Impact of climate change on hydropower generation in Rio Jubones Basin, Ecuador. Water Sci. Eng. 2018, 11, 157–166. [Google Scholar] [CrossRef]
  12. Falchetta, G.; Gernaat, D.E.H.J.; Hunt, J.; Sterl, S. Hydropower dependency and climate change in sub-Saharan Africa: A nexus framework and evidence-based review. J. Clean. Prod. 2019, 231, 1399–1417. [Google Scholar] [CrossRef] [Green Version]
  13. Lu, S.; Dai, W.; Tang, Y.; Guo, M. A review of the impact of hydropower reservoirs on global climate change. Sci. Total Environ. 2020, 711, 134996. [Google Scholar] [CrossRef]
  14. IPCC; McDowell, G.; Barr, I. ‘High Mountain Areas’ Chapter—IPCC Special Report on the Oceans and Cryosphere in a Changing Climate (SROCC); IPCC: Geneva, Switzerland, 2019. [Google Scholar]
  15. Mishra, S.K.; Hayse, J.; Veselka, T.; Yan, E.; Kayastha, R.B.; LaGory, K.; McDonald, K.; Steiner, N. An integrated assessment approach for estimating the economic impacts of climate change on River systems: An application to hydropower and fisheries in a Himalayan River, Trishuli. Environ. Sci. Policy 2018, 87, 102–111. [Google Scholar] [CrossRef]
  16. Feldbauer, J.; Kneis, D.; Hegewald, T.; Berendonk, T.U.; Petzoldt, T. Managing climate change in drinking water reservoirs: Potentials and limitations of dynamic withdrawal strategies. Environ. Sci. Eur. 2020, 32, 48. [Google Scholar] [CrossRef] [Green Version]
  17. Helfer, F.; Lemckert, C.; Zhang, H. Impacts of climate change on temperature and evaporation from a large reservoir in Australia. J. Hydrol. 2012, 475, 365–378. [Google Scholar] [CrossRef] [Green Version]
  18. Shu, J.; Qu, J.; Motha, R.; Xu, J.; Dong, D. Impacts of climate change on hydropower development and sustainability: A review. IOP Conf. Ser. Earth Environ. Sci. 2018, 163, 012126. [Google Scholar] [CrossRef]
  19. Turner, S.W.D.; Hejazi, M.; Kim, S.H.; Clarke, L.; Edmonds, J. Climate impacts on hydropower and consequences for global electricity supply investment needs. Energy 2017, 141, 2081–2090. [Google Scholar] [CrossRef]
  20. Nam, W.-H.; Kim, T.; Hong, E.-M.; Choi, J.-Y. Regional Climate Change Impacts on Irrigation Vulnerable Season Shifts in Agricultural Water Availability for South Korea. Water 2017, 9, 735. [Google Scholar] [CrossRef] [Green Version]
  21. Carvajal, P.E.; Li, F.G.N.; Soria, R.; Cronin, J.; Anandarajah, G.; Mulugetta, Y. Large hydropower, decarbonisation and climate change uncertainty: Modelling power sector pathways for Ecuador. Energy Strategy Rev. 2019, 23, 86–99. [Google Scholar] [CrossRef]
  22. Chilkoti, V.; Bolisetti, T.; Balachandar, R. Climate change impact assessment on hydropower generation using multi-model climate ensemble. Renew. Energy 2017, 109, 510–517. [Google Scholar] [CrossRef]
  23. Sample, J.E.; Duncan, N.; Ferguson, M.; Cooksley, S. Scotland׳s hydropower: Current capacity, future potential and the possible impacts of climate change. Renew. Sustain. Energy Rev. 2015, 52, 111–122. [Google Scholar] [CrossRef]
  24. Pereira-Cardenal, S.; Madsen, H.; Arnbjerg-Nielsen, K.; Riegels, N.; Jensen, R.; Mo, B.; Wangensteen, I.; Bauer-Gottwein, P. Assessing climate change impacts on the Iberian power system using a coupled water-power model. Clim. Chang. 2014, 126, 351–364. [Google Scholar] [CrossRef]
  25. González-Villela, R.; Martínez, M.J.M.; Sepúlveda, J.S.S. Effects of climate change on the environmental flows in the Conchos River (Chihuahua, Mexico). Ecohydrol. Hydrobiol. 2018, 18, 431–440. [Google Scholar] [CrossRef]
  26. Givati, A.; Thirel, G.; Rosenfeld, D.; Paz, D. Climate change impacts on streamflow at the upper Jordan River based on an ensemble of regional climate models. J. Hydrol. Reg. Stud. 2019, 21, 92–109. [Google Scholar] [CrossRef]
  27. Clifton, C.F.; Day, K.T.; Luce, C.H.; Grant, G.E.; Safeeq, M.; Halofsky, J.E.; Staab, B.P. Effects of climate change on hydrology and water resources in the Blue Mountains, Oregon, USA. Clim. Serv. 2018, 10, 9–19. [Google Scholar] [CrossRef]
  28. Bhatta, B.; Shrestha, S.; Shrestha, P.K.; Talchabhadel, R. Evaluation and application of a SWAT model to assess the climate change impact on the hydrology of the Himalayan River Basin. CATENA 2019, 181, 104082. [Google Scholar] [CrossRef]
  29. de Queiroz, A.R.; Faria, V.A.D.; Lima, L.M.M.; Lima, J.W.M. Hydropower revenues under the threat of climate change in Brazil. Renew. Energy 2019, 133, 873–882. [Google Scholar] [CrossRef]
  30. Arango-Aramburo, S.; Turner, S.W.D.; Daenzer, K.; Ríos-Ocampo, J.P.; Hejazi, M.I.; Kober, T.; Álvarez-Espinosa, A.C.; Romero-Otalora, G.D.; van der Zwaan, B. Climate impacts on hydropower in Colombia: A multi-model assessment of power sector adaptation pathways. Energy Policy 2019, 128, 179–188. [Google Scholar] [CrossRef]
  31. Zhai, M.Y.; Lin, Q.G.; Huang, G.H.; Zhu, L.; An, K.; Li, G.C.; Huang, Y.F. Adaptation of Cascade Hydropower Station Scheduling on A Headwater Stream of the Yangtze River under Changing Climate Conditions. Water 2017, 9, 293. [Google Scholar] [CrossRef] [Green Version]
  32. Zhang, W.; Liu, P.; Wang, H.; Lei, X.; Feng, M. Operating rules of irrigation reservoir under climate change and its application for the Dongwushi Reservoir in China. J. Hydro-Environ. Res. 2017, 16, 34–44. [Google Scholar] [CrossRef]
  33. Yang, X.S.; Deb, S. Cuckoo Search via Levy Flights. Mathematics 2010, 210–214. [Google Scholar]
  34. Geressu, R.; Harou, J. Reservoir system expansion scheduling under conflicting interests. Environ. Model. Softw. 2019, 118, 201–210. [Google Scholar] [CrossRef]
  35. Feng, Y.; Zhou, J.; Mo, L.; Wang, C.; Yuan, Z.; Wu, J. A Gradient-Based Cuckoo Search Algorithm for a Reservoir-Generation Scheduling Problem. Algorithms 2018, 11, 36. [Google Scholar] [CrossRef] [Green Version]
  36. Warszawski, L.; Frieler, K.; Huber, V.; Piontek, F.; Serdeczny, O.; Schewe, J. The Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP): Project framework. Proc. Natl. Acad. Sci. USA 2014, 111, 3228–3232. [Google Scholar] [CrossRef] [Green Version]
  37. Hempel, S.; Frieler, K.; Warszawski, L.; Schewe, J.; Piontek, F. A trend-preserving bias correction—The ISI-MIP approach. Earth Syst. Dyn. 2013, 4, 219–236. [Google Scholar] [CrossRef] [Green Version]
  38. Berti, A.; Tardivo, G.; Chiaudani, A.; Rech, F.; Borin, M. Assessing reference evapotranspiration by the Hargreaves method in north-eastern Italy. Agric. Water Manag. 2014, 140, 20–25. [Google Scholar] [CrossRef]
  39. Feng, Y.; Zhou, J.; Mo, L.; Yuan, Z.; Zhang, P.; Jiang, W.; Wang, C.; Wang, Y. Long-Term Hydropower Generation of Cascade Reservoirs under Future Climate Changes in Jinsha River in Southwest China. Water 2018, 10, 235. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flowchart for assessing the impact of climate change on hydropower generation.
Figure 1. Flowchart for assessing the impact of climate change on hydropower generation.
Water 13 01803 g001
Figure 2. Flowchart of MoCS.
Figure 2. Flowchart of MoCS.
Water 13 01803 g002
Figure 3. Gradient search strategy.
Figure 3. Gradient search strategy.
Water 13 01803 g003
Figure 4. Flowchart of gradient-based search strategy.
Figure 4. Flowchart of gradient-based search strategy.
Water 13 01803 g004
Figure 5. Single entry external archive.
Figure 5. Single entry external archive.
Water 13 01803 g005
Figure 6. Self-tuning divergent operator strategy.
Figure 6. Self-tuning divergent operator strategy.
Water 13 01803 g006
Figure 7. Flowchart of GMoCS for solving MLTHG.
Figure 7. Flowchart of GMoCS for solving MLTHG.
Water 13 01803 g007
Figure 8. Locations of major cities and hydrological stations along the Jinsha River.
Figure 8. Locations of major cities and hydrological stations along the Jinsha River.
Water 13 01803 g008
Figure 9. Change of precipitation.
Figure 9. Change of precipitation.
Water 13 01803 g009
Figure 10. Change of maximum temperature.
Figure 10. Change of maximum temperature.
Water 13 01803 g010
Figure 11. Change of minimum temperature.
Figure 11. Change of minimum temperature.
Water 13 01803 g011
Figure 12. Locations of major cities and hydrological stations along the Jinsha River.
Figure 12. Locations of major cities and hydrological stations along the Jinsha River.
Water 13 01803 g012
Figure 13. Topology of hydrological stations in the Jinsha River.
Figure 13. Topology of hydrological stations in the Jinsha River.
Water 13 01803 g013
Figure 14. Three typical streamflows of Pingshan station.
Figure 14. Three typical streamflows of Pingshan station.
Water 13 01803 g014
Figure 15. Optimal results obtained by MoCS and GMoCS: (a) Wet year; (b) Normal year; (c) Dry year.
Figure 15. Optimal results obtained by MoCS and GMoCS: (a) Wet year; (b) Normal year; (c) Dry year.
Water 13 01803 g015aWater 13 01803 g015b
Figure 16. Comparison between the 2nd and 49th solutions in wet year. (a) Water level of Wudongde (b) Power generation of Wudongde (c) Water level of Baihetan (d) Power generation of Baihetan (e) Water level of Xiluodu (f) Power generation of Xiluodu (g) Water level of Xiangjiaba (h) Power generation of Xiangjiaba.
Figure 16. Comparison between the 2nd and 49th solutions in wet year. (a) Water level of Wudongde (b) Power generation of Wudongde (c) Water level of Baihetan (d) Power generation of Baihetan (e) Water level of Xiluodu (f) Power generation of Xiluodu (g) Water level of Xiangjiaba (h) Power generation of Xiangjiaba.
Water 13 01803 g016
Figure 17. Annual streamflow in Pingshan Station under IPCC climate scenarios: (a) RCP2.6, (b) RCP4.5, and (c) RCP8.5.
Figure 17. Annual streamflow in Pingshan Station under IPCC climate scenarios: (a) RCP2.6, (b) RCP4.5, and (c) RCP8.5.
Water 13 01803 g017
Figure 18. Power generation projections for the next three decades under scenarios: (a) RCP2.6, (b) RCP4.5, and (c) RCP8.5.
Figure 18. Power generation projections for the next three decades under scenarios: (a) RCP2.6, (b) RCP4.5, and (c) RCP8.5.
Water 13 01803 g018aWater 13 01803 g018b
Figure 19. Residual load variance for the next three decades under scenarios: (a) RCP2.6, (b) RCP4.5, and (c) RCP8.5.
Figure 19. Residual load variance for the next three decades under scenarios: (a) RCP2.6, (b) RCP4.5, and (c) RCP8.5.
Water 13 01803 g019aWater 13 01803 g019b
Figure 20. Predicted, decadal variation ratios under different climate scenarios: (a) RCP2.6, (b) RCP4.5, and (c) RCP8.5.
Figure 20. Predicted, decadal variation ratios under different climate scenarios: (a) RCP2.6, (b) RCP4.5, and (c) RCP8.5.
Water 13 01803 g020aWater 13 01803 g020b
Table 1. Main parameters of the four cascade hydropower stations along the lower reaches of the Jinsha River.
Table 1. Main parameters of the four cascade hydropower stations along the lower reaches of the Jinsha River.
ParametersWudongdeBaihetanXiluoduXiangjiaba
Dead water level (m)945765540370
Normal water level (m)977825600380
Flood limit water level (m)952785560370
Installed capacity (104 kw)102016001260600
Total capacity (108 m3)74.08206.27126.751.63
Minimum outflow (m3·s−1)90690515001500
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Feng, Y.; Xu, J.; Hong, Y.; Wang, Y.; Yuan, Z.; Wang, C. Reservoir Scheduling Using a Multi-Objective Cuckoo Search Algorithm under Climate Change in Jinsha River, China. Water 2021, 13, 1803. https://doi.org/10.3390/w13131803

AMA Style

Feng Y, Xu J, Hong Y, Wang Y, Yuan Z, Wang C. Reservoir Scheduling Using a Multi-Objective Cuckoo Search Algorithm under Climate Change in Jinsha River, China. Water. 2021; 13(13):1803. https://doi.org/10.3390/w13131803

Chicago/Turabian Style

Feng, Yu, Jijun Xu, Yang Hong, Yongqiang Wang, Zhe Yuan, and Chao Wang. 2021. "Reservoir Scheduling Using a Multi-Objective Cuckoo Search Algorithm under Climate Change in Jinsha River, China" Water 13, no. 13: 1803. https://doi.org/10.3390/w13131803

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop