Next Article in Journal
Short-Term Effects of Forest Fire on Water Quality along a Headwater Stream in the Immediate Post-Fire Period
Previous Article in Journal
Flood Vulnerability Study of a Roadway Bridge Subjected to Hydrodynamic Actions, Local Scour and Wood Debris Accumulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Flood Control Optimization of Reservoir Group Based on Improved Sparrow Algorithm (ISSA)

School of Water Resources, North China University of Water Resources and Electric Power, Zhengzhou 450011, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(1), 132; https://doi.org/10.3390/w15010132
Submission received: 15 November 2022 / Revised: 9 December 2022 / Accepted: 24 December 2022 / Published: 29 December 2022
(This article belongs to the Section Hydrology)

Abstract

:
The optimal control problem of reservoir group flood control is a complex, nonlinear, high-dimensional, multi-peak extremum problem with many complex constraints and interdependent decision variables. The traditional algorithm is slow and easily falls into the local optimum when solving the problem of the flood control optimization of reservoir groups. The intelligent algorithm has the characteristics of fast computing speed and strong searching ability, which can make up for the shortcomings of the traditional algorithm. In this study, the improved sparrow algorithm (ISSA) combining Cauchy mutation and reverse learning strategy is used to solve the flood control optimization problem of reservoir groups. This study takes Sanmenxia Reservoir and Xiaolangdi Reservoir on the mainstream of the Yellow River as the research object and Huayuankou as the downstream control point to establish a joint flood control optimization operation model of cascade reservoirs. The results of the improved sparrow algorithm (ISSA), particle swarm optimization (POS) and sparrow algorithm (SSA) are compared and analyzed. The results show that when the improved ISSA algorithm is used to solve the problem, the maximum flood peak flow of the garden entrance control point is 11,676.3 m3, and the peak cutting rate is 48%. The optimization effect is obviously better than the other two algorithms. This study provides a new and effective way to solve the problem of flood control optimization of reservoir groups.

1. Introduction

Flooding has always been one of the most serious, frequent and extensive natural disasters faced by human beings, posing a serious threat to human life and property [1]. In recent years, under the influence of global climate change, floods have become more frequent and extreme. The reservoir plays an important role in flood management in flood season, as it helps to reduce flood peak, reduce flood damage, and prevent floods, making it an important flood control engineering measure [2]. Reservoir flood control operation, as a non-engineering measure combined with reservoir, has always been a hot field of research. We can make a reasonable dispatching scheme to make the reservoir play its role in flood control [3]. Therefore, it is necessary to study a suitable method for reservoir group joint operation.
Reservoir flood control optimization is a complex and conflicting optimization problem with many interdependent decision variables [4]. The conflict of optimal operation of reservoir flood control mainly focuses on the safety of the upstream dam and downstream dike. When ensuring the safety of the upstream dam, the reservoir is required to discharge as much flood as possible to ensure that the water level is maintained at the safety line. In order to ensure the safety of the downstream river, the reservoir is required to store as much water as possible to reduce downstream risks. Obviously, in a flood control decision, to take care of both requires the decision maker to use a “balanced” solution [5]. In order to solve the above problems, early scholars used dynamic programming (DP) [6,7] and linear programming (LP) [8,9]. However, with the increase in computing dimensions, problems such as dimension disaster and slow computing speed would appear [10,11].
With the rapid development of computers, some heuristic intelligent optimization algorithms have provided some new ideas for flood control scheduling decision. Intelligent optimization algorithms such as genetic algorithm (GA) [12,13] have been gradually applied to solve the problem of optimal flood control scheduling of reservoirs. Subsequently, some scholars successively applied ant colony algorithm (ACO) [14,15], differential evolution algorithm (DE) [16], particle swarm optimization algorithm (PSO) [17,18], and so on, to solve the flood control optimization and dispatching problem of reservoir. Although the swarm intelligent algorithm makes up for the shortcomings of traditional optimization algorithms, the swarm intelligent algorithm often has the disadvantages of precocious convergence and poor convergence performance when solving problems, and it cannot ensure the convergence to the best advantage [19].
The sparrow algorithm (SSA) is a new swarm intelligence algorithm, proposed by Xue Jiankai [20] in 2020, which is inspired by sparrow’s foraging behavior and anti-predation behavior. This algorithm is novel, has the advantages of strong optimization searching ability and fast convergence speed, and has been widely used in many fields [21]. In view of the problem that the population diversity of the basic sparrow search algorithm (SSA) decreases in the late iteration period and is easy to fall into local extremum, Mao Qinghua proposed an improved sparrow algorithm (ISSA) that combines Cauchy mutation and reverse learning [22]. However, the basic sparrow search algorithm (SSA) also faces the problem of falling into local extremum when the population diversity decreases in the late iteration period. To solve the above problems, Mao Qinghua proposed an improved sparrow algorithm (ISSA) that combines Cauchy mutation and reverse learning.
In this paper, the ISSA algorithm is firstly used to solve the problem of flood control operation of reservoirs. In order to verify the feasibility and effectiveness of the ISSA algorithm in the joint flood control operation of reservoirs, this paper takes Sanmenxia Reservoir and Xiaolangdian Reservoir in the mainstream of the Yellow River as the research objects and establishes the maximum flood peak cutting model of Huayuankou. The ISSA, SSA and POS algorithms are applied to solve this model, and the calculation results of the three algorithms are compared. The results show that the ISSA algorithm is superior to other algorithms selected in this paper, which provides a new method to solve the joint flood control scheduling problem of reservoir groups.
The structure of the paper is as follows. The second chapter introduces the optimal operation model of reservoir joint flood control. The third chapter introduces the sparrow algorithm (SSA) and the improved sparrow algorithm (ISSA). The fourth chapter is the model analysis part, which mainly introduces the solution steps of the flood control operation model, the research area of the flood control operation model, and the interval flood analysis of the flood control operation model. The fifth chapter is the results and discussion. The sixth chapter is the conclusion of the article.

2. Reservoir Joint Flood Control Dispatching Model

2.1. Objective Function

The flood control system consists of embankments, flood diversion projects and reservoirs. During flood control regulation, it is necessary to give full play to the advantages of various projects and control and regulate flood in a planned and unified way. In flood season, reservoir flood control operation rules are usually used to guide the operation of the reservoir in flood season according to the current storage state and inflow flow of the reservoir. In this paper, the maximum peak clipping criterion is used to establish the flood control optimal operation model of cascade reservoirs. Assuming that there are N reservoirs in series, the reservoir inflow and interval flood are known, and the total number of dispatching periods is T. See Figure 1 for the schematic diagram of cascade reservoirs of N level.
The objective function established by the maximum peak clipping criterion of control points:
f ( q ) = min t = 1 T ( q 1 t + Q h 1 ) 2 + ( q 2 t + Q h 2 ) 2 + + ( q N t + Q hN ) 2
where T is the total number of scheduling periods; t is the number of periods; q 1 t is the average discharge of the first upstream reservoir in t period; q 2 t is the average discharge of the second upstream reservoir in t period; q N t is the average discharge of the Nth upstream reservoir in period t; Q h is the average interval inflow during t period.

2.2. Condition of Constraint

(1)
Water balance constraint:
V t = V t 1 + ( ( Q t + Q t + 1 ) / 2 ( q t + q t + 1 ) / 2 ) Δ t
V t is the average storage capacity of time period t; V t 1 is the average storage capacity in the T − 1 period; Q t , Q t + 1 is the average inbound flow in time period t and time period t + 1; q t , q t + 1 is the average discharge volume in time period t and t + 1.
(2)
Constraints on reservoir water level:
Z min , k Z k Z max , k ( k = 1 , 2 )
where T represents the water level of the k reservoir during flood control operations (m); Z min , k stands for the lowest allowable water level during flood control operations, which is the flood limit water level; Z max , k T stands for the maximum allowable water level during flood control operations, which is the check water level;
(3)
Letdown flow constraint:
0 q k q max , k ( k = 1 , 2 )
where q k is the discharge of the kth reservoir; q max , k is the maximum discharge capacity of the kth reservoir.
(4)
Initial water level constraint:
Z 0 = Z min , k ( k = 1 , 2 )
where Z 0 is the initial water level; Z min , k is the lowest allowable water level during flood control operation.
(5)
Terminal water level constraint:
Z e n d = Z min , k ( k = 1 , 2 )
where Z e n d is the water level of the reservoir at the end of flood control regulation; Z min , k is the lowest allowable water level during flood control operation.

3. Sparrow Optimization Algorithm

3.1. Basic Sparrow Search Algorithm (SSA)

The sparrow algorithm (SSA) was first proposed by Xue Jiankai [20] in 2020, which is a new intelligent optimization algorithm based on the feeding and anti-predation behavior of the sparrow population. According to the position and energy level of sparrows in the population, sparrows are mainly divided into three categories: discoverer, joiner and watcher. Among them, the discoverer is responsible for finding food, the joiner is responsible for catching prey after the discoverer finds food, and the watcher is responsible for sensing natural enemies to ensure population safety.
The population matrix form composed of n sparrows can be expressed as follows:
X = x 1 , 1 x 1 , 2 x 1 , d x 2 , 1 x 2 , 2 x 2 , d x n , 1 x n , 2 x n , d
where d represents the dimension of the variable to be optimized, and n represents the number of sparrows. Then, the fitness values of all sparrows can be expressed as follows:
F x = f ( [ x 1 , 1 x 1 , 2 x 1 , n ] ) f ( [ x 2 , 1 x 2 , 2 x 2 , n ] ) f ( [ x n , 1 x n , 2 x n , n ] )
where f represents the fitness value.
In the sparrow algorithm (SSA), the discoverer with better fitness value will give priority to getting food during the search process. In addition, the discoverer is responsible for finding food for the whole sparrow population and providing the direction for all participants. During each iteration, the location of the discoverer is updated as follows:
X i , j t + 1 = X i , j t exp ( i α i t e r max ) i f   R 2 < S T X i , j t + Q L i f   R 2 < S T
where t represents the current iteration algebra, j = 1, 2, 3, …, d. i t e r max is a constant representing the maximum number of iterations. X i , j represents the position information of the ith sparrow in the jth dimension. α ∈ (0, 1] is a random number. R 2 ( R 2 ∈ [0, 1]) and S T ( S T ∈ [0.5, 1]) represent early warning value and safety value, respectively. Q is a random number subject to normal distribution. L represents a 1 × d, where each element in the matrix is 1. When R 2 < S T , it means that there are no predators around the foraging environment at this time, and the finder can perform extensive search operations. If R 2 S T , it means that some sparrows in the population have found predators and sent an alarm to other sparrows in the population. At this time, all sparrows need to quickly fly to other safe places for feeding.
For predators, in the process of foraging, some participants will constantly monitor the finder. Once they realize that the discoverer has found better food, they will immediately leave their current position to compete for food. The location updates of the participants are described as follows:
X i , j t + 1 = Q exp ( X w o r s t X i , j t i 2 ) X p t + 1 + | X i , j t X p t + 1 | A + L i f     i > n / 2 o t h e r w i s e
where X p is the optimal position occupied by the current discoverer; X w o r s t represents the current global worst position. A represents a 1 × d, and each element is randomly assigned to 1 or −1 and A + = AT (AAT)−1. When i > n/2, it indicates that the i-th finder with low fitness value has not got food and is in a very hungry state. At this time, it needs to fly to other places to forage for more energy.
In the simulation, we assume 10 to 20 percent of the total number of watchmen. The initial positions of these sparrows were randomly generated within the population. Its mathematical expression can be shown as follows:
X i , j t + 1 = X b e s t t + β | X i , j t X b e s t t | X i , j t + K ( | X i , j t X w o r s t t | ( f i f w ) + ε ) i f   f i > f g i f   f i = f g
where X b e s t is the current global optimal position. As a step control parameter, β is a random number subject to the normal distribution with mean value of 0 and variance of 1. K ∈ [−1, 1] is a random number, and f i is the fitness value of the current sparrow individual. f g and f w are the current global best and worst fitness values, respectively. ε is the smallest constant to avoid zeros in the denominator.
For simplicity, when f i > f g indicates that the sparrow is at the edge of the population and is extremely vulnerable to predators. X b e s t means that the sparrow in this position is the best and safe position in the population. When f i = f g , it indicates that the sparrows in the middle of the population are aware of the danger and need to get close to other sparrows to minimize their risk of being preyed on. K indicates that the direction of movement of the sparrow is also a step control parameter.

3.2. Improved Sparrow Algorithm (ISSA)

(1)
Sin chaos initialization population
The Sin chaos model is a model with an infinite number of mapping folds. In this paper, Sin chaos is used to initialize the population of the sparrow algorithm. The one-dimensional self-mapping expression of Sin chaos is as follows:
x n + 1 = sin 2 x n , n = 0 , 1 , , N 1 x n 1 , x n 0
In Formula (12), to prevent fixed points and zeros from being generated in [−1, 1] [−1, 1] [−1, 1], the initial value cannot be set to 0.
(2)
Dynamic adaptive weight
With reference to the idea of inertia weight, the dynamic weight factor ω is introduced into the discoverer position updating mode to make it have a larger value in the early iteration, which can be better explored globally, and reduce adaptively in the late iteration, so as to better perform local search and improve the convergence speed. The calculation formula of the weight coefficient ω and the updated method of the discoverer’s position after improvement are as follows:
ω = e 2 ( 1 t / i t e r max ) e 2 ( 1 t / i t e r max ) e 2 ( 1 t / i t e r max ) + e 2 ( 1 t / i t e r max )
X i , j t + 1 = ( X i , j t + ω ( f j , g t X i , j t ) ) r a n d , R 2 < S T X i , j t + Q , R 2 S T
In Formula (14), f j , g t is the global optimal solution of the jth dimension in the previous generation.
(3)
An improved formula for updating the position of the watchman
X i , j t + 1 = X b e s t t + β ( X i , j t X b e s t t ) ,   f i f g X b e s t t + β ( X w o r s t t X b e s t t ) ,   f i = f g
(4)
Fusion of Cauchy Variation and Reverse Learning Strategy
In order to enable individuals to better find the optimal solution, the reverse learning strategy is integrated into the sparrow algorithm. The mathematical representation is as follows:
X b e s t ( t ) = u b + r ( l b X b e s t ( t ) )
X i , j t + 1 = X b e s t ( t ) + b 1 ( X b e s t ( t ) X b e s t ( t ) )
where X b e s t ( t ) is the inverse solution of the optimal solution of generation t. u b , l b is the upper and lower bounds, respectively. r is 1 subject to normal uniform distribution × d (d is the space dimension); b 1 represents the information exchange control parameter, and the formula is as follows:
b 1 = ( i t e r max t i t e r max ) t
To further improve the optimization performance of the algorithm, Cauchy operator perturbation is incorporated into the sparrow algorithm, and the mathematical expression is as follows:
X i , j t + 1 = X b e s t ( t ) + c a u c h y ( 0 , 1 ) X b e s t ( t )
where Cauchy(0, 1) is the standard Cauchy distribution. The generating function of the Cauchy distribution random variable is η = tan[(ξ − 0.5) π].
In order to further improve the optimization performance of the algorithm, a dynamic selection strategy is adopted to update the target position. The reverse learning strategy and Cauchy mutation operator perturbation strategy are alternately executed under a certain probability to dynamically update the target position. In the reverse learning strategy, the reverse solution is obtained through the reverse learning mechanism to expand the search field of the algorithm. In Cauchy mutation strategy, Cauchy mutation operator is used to perturb and mutate at the optimal solution position to obtain a new solution, which improves the defect that the algorithm falls into local area. As for which strategy to adopt for target location update, it is determined by the selection probability P s , and its calculation formula is as follows:
P s = exp ( 1 t i t e r max ) 20 + θ
In the formula, θ is the adjustment parameter, and its value can be 0.05. The specific selection strategy is as follows: if rand < P s , choose (16)~(18) reverse learning strategy to update the position; otherwise, choose (8) Cauchy mutation perturbation strategy to update the target position. Although the above two perturbation strategies can enhance the ability of the algorithm to jump out of local space, it is impossible to determine that the new position obtained after the perturbation mutation is better than the fitness value of the original position. Therefore, after the perturbation mutation update, the greedy rule is introduced to determine whether to update the position by comparing the fitness values of the new and old positions. The greedy rule is shown in Equation (21), and f(x) represents the location fitness value of x.
X b e s t ( t ) = X i , j t + 1 ,   f ( X i , j t + 1 ) < f ( X b e s t ( t ) ) X i , j t + 1 ,   f ( X i , j t + 1 ) f ( X b e s t ( t ) )

4. Model Analysis

4.1. Solution Method of Flood Control Operation Model

This study uses the improved sparrow algorithm (ISSA) in Section 3.2 to solve the optimal operation model of Sanmenxia and Xiaolangdi reservoirs for joint flood control. The solution steps of the model are shown in Figure 2.
(1)
Initialization parameters, such as population number N, maximum iteration number i t e r max , discoverer proportion PD, joiner proportion SD, alert threshold R2 and initialization of sparrow population using the Sin chaotic map of Equation (12) according to the calculation period given by flood flow. The sparrow is constructed with the discharge of each reservoir at the end of each cycle as the control variable. The flood lasts 78 periods and has two reservoirs, so the vector dimension is 156. Then, the matrix form of sparrow population:
Y ( 0 ) = [ Z 1 ( 0 ) , Z 2 ( 0 ) , X n ( 0 ) ] = z 1 , 1 ( 0 ) z 1 , 2 ( 0 ) z 1 , n ( 0 ) z 2 , 1 ( 0 ) z 2 , 2 ( 0 ) z 2 , n ( 0 ) z 156 , 1 ( 0 ) z 156 , 2 ( 0 ) z 156 , n ( 0 )
Z i ( 0 ) is the ith initial sparrow, including 156 elements. 1~78 represents the reservoir water level at the end of the Sanmenxia Reservoir period, and 79~156 represents the end of the Xiaolangdi Reservoir period.
(2)
Calculate the fitness value of each sparrow and find out the current optimal fitness value, the worst fitness value, and the corresponding position.
(3)
From the sparrows with better fitness values, select some sparrows as the discoverers, and update their positions according to Formula (14).
(4)
The remaining sparrows will be the participants, and their positions will be updated in the original way.
(5)
Randomly select some sparrows from the sparrows as watchers and update their positions according to Formula (15).
(6)
According to probability P s , the Cauchy mutation perturbation strategy and reverse learning strategy are selected to perturb the current optimal solution to generate a new solution.
(7)
Determine whether to update the location according to the greedy rule (21).
(8)
Judge whether the end conditions are met. If the end conditions are met, proceed to the next step; otherwise, skip to step (2).
(9)
The program ends and the optimal result is output.

4.2. Study Area of Flood Control Operation Model

This paper selects the section from Sanmenxia Reservoir to Huayuankou in the middle and lower reaches of the Yellow River as the research area. The Yellow River is a big river in northern China, one of the longest rivers in the world and the second longest river in China. The total length of the Yellow River is about 5464 km, and its drainage area is about 752,443 km2. The Yellow River basin has many tributaries and developed water system, which also leads to frequent floods in the Yellow River basin. The water system of the Yellow River basin is shown in Figure 3.
This paper chooses Sanmenxia and Xiaolangdi Reservoir as the research objects. Sanmenxia Reservoir and Xiaolangdi Reservoir are two large reservoirs located in the middle and lower reaches of the mainstream of the Yellow River. They are the backbone projects of the Yellow River hydraulic engineering system. They are of great significance to the development and utilization of water conservancy and hydropower resources in the Yellow River basin and the relief of flood control pressure in the lower reaches of the Yellow River. Under the condition of ensuring the flood control safety of Huayuankou, the optimal operation of Sanmenxia Reservoir and Xiaolangdi Reservoir in the Yellow River basin is studied. See Figure 4 for the locations of Sanmenxia, Xiaolangdi Reservoir and Huayuankou control points.
Sanmenxia Reservoir is a key project in the middle reaches of the Yellow River. It is 120 km away from Tongguan Station and 260 km away from Huayuankou Station. The total storage capacity of the reservoir is 16.2 billion cubic meters. The drainage area controlled by the reservoir is 688,400 km2, accounting for 91.49% of the total drainage area. Sanmenxia Reservoir is a large comprehensive reservoir focusing on flood control and considering irrigation and power generation. Xiaolangdi Reservoir is 130 km away from Sanmenxia Reservoir. The total storage capacity of the reservoir is 12.65 billion cubic meters, and the controlled drainage area is 694,000 km2, accounting for 92.3% of the total area of the Yellow River basin. It is the only large-scale comprehensive water conservancy project with large storage capacity in the middle and lower reaches of the Yellow River except Sanmenxia. See Table 1 for the current situation of Sanmenxia Xiaolangdi Reservoir.

4.3. Interval Water Supply Analysis of Flood Control Operation Model

Taking Huayuankou as the control point, the upstream can be divided into three areas according to the location of Sanmenxia reservoir and Xiaolangdi reservoir. Therefore, there are three regional flood processes in the mainstream of the Yellow River above Huayuankou, including the upstream flood of Sanmenxia Reservoir, the interval flood between Sanmenxia and Xiaolangdi Reservoir, and the interval flood between Xiaolangdi and Huayuankou. Assume that the inflow flow of Sanmenxia is Q1, the interval flow between Sanmenxia and Xiaolangdi is Q2, and the interval flow between Xiaolangdi Reservoir and Huayuankou is Q3. It is assumed that the discharge of Sanmenxia Reservoir is q1, and the discharge will evolve into q2 through the river channel. The sum of interval flow Q2 and flow q2 of Sanmenxia Reservoir and Xiaolangdi Reservoir is the inflow flow of Xiaolangdi Reservoir. The outflow flow q3 of Xiaolangdi reservoir evolves into flow q4 through the river, and the sum of q4 and flow Q3 is used as the flood flow at Huayuan Estuary of the control point. The flood routing process of the flood control system is shown in Figure 5.
For the calculation of evolution flood, the Muskin method is adopted. Its basic principle: ignoring the inertia term, the dynamic equation can be simplified as the tank storage equation, as shown in the following formula:
W = K [ x I + ( 1 x ) Q ] Q
1 2 ( I 1 I 2 ) Δ t 1 2 ( O 1 O 2 ) Δ t = W 2 W 1
where W is the water storage capacity of the river, and Q is the storage flow. K is the slope of the storage flow relationship, and x is the flow specific gravity coefficient. I 1 , I 2 are the inflow of upstream section at the beginning and end of the calculation period. O 1 , O 2 are the downstream section flow at the beginning and end of the calculation period. Δ t is the length of the calculation period. Simultaneous water balance equation and tank storage equation are solved:
O 2 = C 0 I 2 + C 1 I 1 + C 2 O 1
C 0 = 1 2 Δ t K x K K x + 1 2 Δ t
C 1 = 1 2 Δ t + K x K K x + 1 2 Δ t
C 2 = K K x 1 2 Δ t K K x + 1 2 Δ t
C 0 + C 1 + C 2 = 1
Due to the different characteristics of river courses, the values of parameters K and x are different for different research river sections. See Table 2 for the values of relevant parameters.

5. Results and Discussion

The real data of 58 years typical millennium flood in the mainstream of the Yellow River are selected for this study. The authenticity and reliability of the research are guaranteed. The total duration of this flood is 312 h, which is divided into 78 periods of 4 h each. The initial population number of ISSA algorithm is N = 50, the maximum iteration number is itermax = 1000, the proportion of discoverers is PD = 0.7, the proportion of scouts is SD = 0.2, and the alert threshold is R2 = 0.6. Under the same conditions, ISSA, SSA and POS algorithms are used to solve the joint flood control optimal operation model of Sanmenxia and Xiaolangdi reservoirs. Figure 6 and Figure 7 show the operation process of Sanmenxia and Xiaolangdi reservoirs when the ISSA algorithm is used to solve the model. Figure 8 shows the flood process of the Huayuankou control point without reservoir regulation and when three different algorithms are used to solve the model.
It can be seen from Figure 6 and Figure 7 that when the inflow flow of Sanmenxia Reservoir and Xiaolangdi Reservoir is small, the reservoir discharge flow is equal to the inflow flow, and the water level of the reservoir remains basically unchanged. When the inflow of Sanmenxia Reservoir and Xiaolangdi Reservoir is large, considering the discharge capacity of the reservoir and the bearing capacity of the downstream river channel, the reservoir begins to store water and the water level starts to rise, sharing the pressure for the downstream flood control. After the reservoir completes a flood control operation without considering the water storage, the water level should return to the initial regulation level to cope with the arrival of the next flood. From Figure 6 and Figure 7, the water level of Sanmenxia and Xiaolangdi reservoirs returned to the initial regulation level at the end of flood control regulation. It not only ensures the flood control safety of the reservoir itself and the downstream, but also provides guarantee for the next flood control operation of the reservoir.
Figure 8 shows the flood process of Huayuankou, a downstream control point. The black line on the top represents the flood process of control points when there is no reservoir regulation. The red, blue and green lines represent the flood process of control points when the POS algorithm, SSA algorithm and ISSA algorithm are used to solve the model. It can be seen from Figure 8 that when the ISAA, SSA and POS algorithms are used to solve the model, the flood peak at the control point at Huayuankou can be reduced. However, different algorithms have different effects when solving the model. See Table 3 for flood data and algorithm initialization parameters of the Huayuankou control point.
In this paper, the optimal operation model of Sanmenxia and Xiaolangdi reservoirs for joint flood control is established based on the maximum peak cutting criterion of the Huayuankou control point. It can be seen from Table 3 that when ISSA, SSA and POS algorithms are used to solve the problem, the maximum peak flood flow of the garden entrance control point is 11,676.3, 12,673.65 and 12,408.23 m3/s, respectively, meeting the flood control requirements of 22,000 m3/s for this section. When the ISSA algorithm is used to solve the model, the peak clipping rate of the Huayuankou control point is 48%, which is significantly higher than that of the other two algorithms. Additionally, the flood process of the Huayuankou control point is more stable when the ISSA algorithm is used to solve the problem. It shows that the ISSA algorithm is superior to the other two algorithms in solving the joint flood control operation problem of reservoirs.

6. Conclusions

Sanmenxia Reservoir and Xiaolangdi Reservoir are taken as the research object, and Huayuankou is taken as the downstream control point. Based on the maximum peak clipping criterion, the joint flood control operation model of Sanmenxia and Xiaolangdi Reservoirs is established. The main constraints of the model are water level constraints, flow constraints and water balance constraints. The ISSA algorithm is first used to solve the problem of joint flood control operation of reservoirs. The model solution results of the ISSA algorithm, POS algorithm and SSA algorithm under the same conditions are compared and analyzed.
(1)
From the solution results of the model, it can be seen that the construction of this model not only guarantees the flood control safety of the reservoir itself when encountering the flood with the return period of 1000 years, but also guarantees the flood control safety of the downstream, indicating the rationality and applicability of the model.
(2)
In this paper, ISSA, SSA and POS algorithms are used to study the flood control operation of tandem cascade reservoirs including Sanmenxia Reservoir and Xiaolangdi Reservoir on the Yellow River mainstream. The comparison of the solution results of the three optimization algorithms shows that the ISSA algorithm is more efficient than the conventional SSA algorithm and POS algorithm and can fully play the role of reservoir capacity compensation to achieve the best flood control effect. The results show that the ISSA algorithm can effectively solve the flood control optimal operation problem of cascade reservoirs. It provides a new method to solve the problem of joint flood control operation of cascade reservoirs. It also provides a reference for the ISSA algorithm to be applied to other research fields.
(3)
Without the regulation of Sanmenxia Reservoir and Xiaolangdi Reservoir, the maximum flood peak flow at the Huayuankou control point is 24,325 m3/s, far exceeding the flood control requirement of 22,000 m3/s at the Huayuankou control point section. Through the joint flood control operation of Sanmenxia Reservoir and Xiaolangdi Reservoir, the flood peak of the Huayuankou control point is reduced, making the maximum peak flow of the Huayuankou control point 11,676.3 m3/s, ensuring the flood control safety of the Huayuankou control point when encountering the millennium flood. This shows that it is necessary to study the joint flood control operation of Sanmenxia Reservoir and Xiaolangdi Reservoir, and further explains the significance of this study.

Author Contributions

J.H. and H.-T.C., were responsible for the original concept and writing the paper. S.-L.W. and X.-Q.G. were processed the data and conducted the program design. Y.-R.W. revised the manuscript and shared numerous comments and suggestions to improve the study quality. S.-M.L.: Methodology, Writing—Original draft. All authors have read and agreed to the published version of the manuscript.

Funding

The authors are grateful to acknowledge the funding support of Project of key science and technology of the Henan province (No: 202102310259; No: 202102310588), Henan province university scientific and technological innovation team (No: 18IRTSTHN009).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Chen, W.; Guo, X.; Law, T.; Cui, Z. Flood risk analysis in protective areas in Yellow River downstream. China Water Resour. 2017, 68, 56–58. [Google Scholar] [CrossRef]
  2. Su, C.; Wang, P.; Yuan, W.; Cheng, C.; Data, T.Z.; Yan, D.; Wu, Z. An MILP based optimization model for reservoir flood control operation considering spillway gate scheduling. J. Hydrol. 2022, 613, 128483. [Google Scholar] [CrossRef]
  3. Li, J.; Zhong, P.-A.; Wang, Y.; Yang, M.; Fu, J.; Liu, W.; Xu, B. Risk analysis for the multi-reservoir flood control operation considering model structure and hydrological uncertainties. J. Hydrol. 2022, 612, 128263. [Google Scholar] [CrossRef]
  4. Zhang, Z.; He, X.; Geng, S.; Zhuang, S.; Li, H.; Tian, Y. Study on flood control projects joint flood control operation of The Daqing river basin. IOP Conf. Ser.: Earth Environ. Sci. 2021, 826, 012007. [Google Scholar] [CrossRef]
  5. Zhang, J.; Cai, X.; Lei, X.; Liu, P.; Wang, H. Real-time reservoir flood control operation enhanced by data assimilation. J. Hydrol. 2021, 598, 126426. [Google Scholar] [CrossRef]
  6. Li, H.; Liu, J. The Application of Dynamic Programming Algorithm in the Optimal Dispatching of the Xiaolangdi Reservoir. Henan Sci. 2011, 29, 461–465. [Google Scholar] [CrossRef]
  7. Wang, S.; Jiang, Z.; Liu, Y. Dimensionality Reduction Method of Dynamic Programming under Hourly Scale and Its Application in Optimal Scheduling of Reservoir Flood Control. Energies 2022, 15, 676. [Google Scholar] [CrossRef]
  8. Xin, Y.; Zhou, Y.; Gu, Q.; Zhang, T. Optimal Scheduling of Day-ahead Considering the Safety Difference of Energy Storage Device Capacity. Sci. Technol. Eng. 2022, 22, 2291–2297. [Google Scholar]
  9. Lee, S.-Y.; Hamlet, A.F.; Fitzgerald, C.J.; Burges, S.J. Optimized Flood Control in the Columbia River Basin for a Global Warming Scenario. J. Water Resour. Plann. Manag. 2009, 135, 440–450. [Google Scholar] [CrossRef] [Green Version]
  10. Cuevas-Velásquez, V.; Sordo-Ward, A.; García-Palacios, J.H.; Bianucci, P.; Garrote, L. Probabilistic Model for Real-Time Flood Operation of a Dam Based on a Deterministic Optimization Model. Water 2020, 12, 3206. [Google Scholar] [CrossRef]
  11. Sumi, T.; Mitsunari, M.; Hamaguchi, T. Challenges in Flood Control Operation and dissemination of Related Information-Lessons from the record-breaking heavy rain in July 2018, Japan. Water Energy Int. 2021, 64, 14–21. [Google Scholar]
  12. Huang, K.; Ye, L.; Chen, L.; Wang, Q.; Dai, L.; Zhou, J.; Singh, V.P.; Huang, M.; Zhang, J. Risk analysis of flood control reservoir operation considering multiple uncertainties. J. Hydrol. 2018, 565, 672–684. [Google Scholar] [CrossRef]
  13. Ren, M.; Zhang, Q.; Yang, Y.; Wang, G.; Xu, W.; Zhao, L. Research and application of joint optimal operation for flood defense of reservoir group based on improved genetic algorithm. China Flood Drought Manag. 2022, 32, 21–26. [Google Scholar] [CrossRef]
  14. Wu, Z.; Zhou, J.; Yang, J. Flood Dispatching Optimization of Three Gorges Region Based on Ant Colony Algorithms. Water Power 2008, 55, 5–7. [Google Scholar] [CrossRef]
  15. Lin, Z. Research on optimal operation of cascade reservoirs based on improved ant colony algorithm. Telecom World 2016, 23, 162–163. [Google Scholar]
  16. Qin, H.; Zhou, J.; Wang, G.; Zhang, Y. Multi-objective optimization of reservoir flood dispatch based on multi-objective differential evolution algorithm. J. Hydraul. Eng. 2009, 40, 513–519. [Google Scholar] [CrossRef]
  17. Chen, H.-T.; He, J.; Wang, W.-C.; Chen, X.-N. Simulation of maize drought degree in Xi’an City based on cusp catastrophe model. Water Sci. Eng. 2021, 14, 28–35. [Google Scholar] [CrossRef]
  18. Chen, Q.; Zhang, Y. Research on optimal joint operation of cascade reservoirs based on improved PSO parallel algorithm. Techenical Superv. Water Resour. 2018, 26, 147–149+182. [Google Scholar] [CrossRef]
  19. Chen, J.; Zhong, P.-A.; Liu, W.; Wan, X.-Y.; Yeh, W.W.-G. A multi-objective risk management model for real-time flood control optimal operation of a parallel reservoir system. J. Hydrol. 2020, 590, 125264. [Google Scholar] [CrossRef]
  20. Malinowski, K.; Karbowski, A.; Salewicz, K.A. Hierarchical Control Structures for Real Time Scheduling of Releases in a Multireservoir System During Floods. IFAC Proc. Vol. 1987, 20, 405–412. [Google Scholar] [CrossRef]
  21. Nan, L.; Jiankai, X.; Huisheng, S. A sparrow search algorithm with adaptive t distribution mutation-based path planning of unmanned aerial vehicles. J. Donghua Univ. (Nat. Sci.) 2022, 48, 69–74. [Google Scholar] [CrossRef]
  22. Mao, Q.; Qiang, Z. Improved Sparrow Algorithm Combining Cauchy Mutation and Opposition-Based Learning. J. Front. Comput. Sci. Technol. 2021, 15, 1155–1164. [Google Scholar] [CrossRef]
Figure 1. Diagram of N-level cascade reservoirs.
Figure 1. Diagram of N-level cascade reservoirs.
Water 15 00132 g001
Figure 2. Flow chart of model calculation steps.
Figure 2. Flow chart of model calculation steps.
Water 15 00132 g002
Figure 3. Water system of the Yellow River basin.
Figure 3. Water system of the Yellow River basin.
Water 15 00132 g003
Figure 4. Location of Sanmenxia, Xiaolangdi Reservoir and Huayuankou.
Figure 4. Location of Sanmenxia, Xiaolangdi Reservoir and Huayuankou.
Water 15 00132 g004
Figure 5. Flood evolution process.
Figure 5. Flood evolution process.
Water 15 00132 g005
Figure 6. Sanmenxia Reservoir operation process.
Figure 6. Sanmenxia Reservoir operation process.
Water 15 00132 g006
Figure 7. Xiaolangdi Reservoir operation process.
Figure 7. Xiaolangdi Reservoir operation process.
Water 15 00132 g007
Figure 8. Flood process of Huayuankou control point.
Figure 8. Flood process of Huayuankou control point.
Water 15 00132 g008
Table 1. Current situation of Sanmenxia and Xiaolangdi Reservoirs.
Table 1. Current situation of Sanmenxia and Xiaolangdi Reservoirs.
ReservoirSanmenxia ReservoirXiaolangdi Reservoir
Control watershed area688,400 km2694,000 km2
Flood control limited water level305 m220 m
Flood control high water level335 m275 m
Total reservoir capacity16,200,000,000 m312,650,000,000 m3
Table 2. Relevant parameters.
Table 2. Relevant parameters.
Name of River SectionFlood Propagation TimeNumber of Flood Routing SectionsKxΔt
Sanmenxia–Xiaolangdi8 h23.8750.24 h
Xiaolangdi–Huayuankou12 h34.5760.34 h
Table 3. Flood data of Huayuankou control point and initialization parameters of the algorithm.
Table 3. Flood data of Huayuankou control point and initialization parameters of the algorithm.
AlgorithmPopulation Quantity (N)Maximum Iterations (Itemax)Maximum Peak Flood DischargePeak Clipping Rate
ISSA50100011,676.3 m3/s48%
SSA100100012,673.65 m3/s45%
POS80100012,408.23 m3/s44%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

He, J.; Liu, S.-M.; Chen, H.-T.; Wang, S.-L.; Guo, X.-Q.; Wan, Y.-R. Flood Control Optimization of Reservoir Group Based on Improved Sparrow Algorithm (ISSA). Water 2023, 15, 132. https://doi.org/10.3390/w15010132

AMA Style

He J, Liu S-M, Chen H-T, Wang S-L, Guo X-Q, Wan Y-R. Flood Control Optimization of Reservoir Group Based on Improved Sparrow Algorithm (ISSA). Water. 2023; 15(1):132. https://doi.org/10.3390/w15010132

Chicago/Turabian Style

He, Ji, Sheng-Ming Liu, Hai-Tao Chen, Song-Lin Wang, Xiao-Qi Guo, and Yu-Rong Wan. 2023. "Flood Control Optimization of Reservoir Group Based on Improved Sparrow Algorithm (ISSA)" Water 15, no. 1: 132. https://doi.org/10.3390/w15010132

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