Next Article in Journal
Does Future Climate Bring Greater Streamflow Simulated by the HSPF Model to South Korea?
Next Article in Special Issue
Efficient Urban Inundation Model for Live Flood Forecasting with Cellular Automata and Motion Cost Fields
Previous Article in Journal
Comparison of the MUSLE Model and Two Years of Solid Transport Measurement, in the Bouregreg Basin, and Impact on the Sedimentation in the Sidi Mohamed Ben Abdellah Reservoir, Morocco
Previous Article in Special Issue
Multipurpose Use of Artificial Channel Networks for Flood Risk Reduction: The Case of the Waterway Padova–Venice (Italy)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Operation of Complex Flood Control System Composed of Cascade Reservoirs, Navigation-Power Junctions, and Flood Storage Areas

1
State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China
2
Jiangxi Provincial Institute of Water Sciences, Nanchang 330029, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(7), 1883; https://doi.org/10.3390/w12071883
Submission received: 19 May 2020 / Revised: 23 June 2020 / Accepted: 30 June 2020 / Published: 1 July 2020
(This article belongs to the Special Issue Research on Mathematical Models of Floods)

Abstract

:
As more and more water projects are built on rivers, the flood control operation becomes more complex. Studies on the optimal flood control operation are very important to safeguard human life and property. This study focused on optimizing the operation of a complex flood control system composed of cascade reservoirs, navigation-power junctions, flood storage areas, and flood control points. An optimal model was established to jointly maximize flood peak reduction rates of downstream flood control points. A hybrid algorithm named the Dynamic Programming-Progressive Optimality Algorithm (DP-POA) was used to solve this model, and the middle and lower reaches of the Ganjiang River were selected as a case study. The results show that flood reduction at three downstream flood control points ranged from 1080 to 5359 m3/s for designed floods with different return periods, which increased by about 333~1498 m3/s in comparison with the conventional operation. Considering that the maximum water level of reservoirs using DP-POA and the conventional operation is the same, this indicated that DP-POA can make full use of the reservoirs’ flood control storage to reduce downstream flood peaks. In addition, the flood diversion volume of the flood storage area using DP-POA ranged from 0.33 × 108 to 1.79 × 108 m3 for designed floods with 200-year, 300-year, and 500-year return periods, which is smaller than that using the conventional operation.

1. Introduction

Flood hazards are the most frequent natural disaster in the world and have become intensified due to global warming [1,2]. Flood control measures including engineering measures (e.g., dam and flood storage area construction) and non-engineering measures (e.g., reservoir operation), have been used to reduce the loss of lives and property from risks related to floods [3]. Optimal flood control operation of reservoirs is one of the most important and efficient non-engineering measures [4], which can enhance the basin flood control capacity. At present, many scholars have carried out a lot of research on optimal flood control operation model construction and optimization algorithms.
Many optimal flood control operation models have been proposed for different objectives. There has been some research on flood control operation for multi-reservoir systems in which the concentration has been on minimization of flood peaks at the downstream flood control points and minimizing the maximal reservoir water level [5,6,7,8]. Some research concentrates on the dynamic flood limiting water level to address the conflict between flood control and reservoir benefit [9,10,11,12]. Some other research addresses flood risk assessment and minimizing the flood control risk for cascade reservoirs [13,14,15]. However, most research has mainly focused on optimizing the flood control operation for multi-reservoir systems such as parallel and cascade reservoirs. Few studies have been carried out on complex flood control systems. Comparing with the multi-reservoir systems, a complex flood control system is composed of other water projects (e.g., navigation-power junctions and flood storage areas), local flows, and flood control points, besides reservoirs. Therefore, optimizing the flood control operation of the complex flood control system needs to consider more factors including operation of water projects, the hydraulic connection between reservoirs, and time-lag issues caused by flood routing during the optimal operation periods.
The optimization algorithms include modern intelligence algorithms and traditional optimization methods, which have been widely used to conduct optimal reservoir operation [16]. Modern intelligence algorithms including Genetic Algorithms (GA) [17,18,19], Particle Swarm Optimization (PSO) [20,21,22], Immune Algorithm [23,24], and other novel optimal algorithms such as Gravity Search Algorithm [25], Water Cycle Algorithm (WCA) [26,27], etc. can be combined with simulation models to search for the optimal reservoir operation. The traditional optimization methods including Linear Programming (LP) [28,29,30] and Dynamic Programming (DP) [31,32,33] are also widely used. However, when the number of reservoirs rises, the computation scale of DP increases greatly, which is called the “curse of dimensionality” and this limits the application of DP. To deal with the dimensionality problems, modified DP algorithms such as Discrete Differential Dynamic Programming (DDDP) [34,35], Dynamic Programming Successive Approximation (DPSA) [36], and Progressive Optimality Algorithm (POA) [37,38,39] have been widely used to find the optimal or near-optimal reservoir release hydrographs.
With the number of reservoirs rising, the number of decision variables and constraints increases and hydraulic connection of multi-reservoir systems becomes more complex, which makes it difficult for a single traditional optimization method or modern intelligence algorithm to optimize reservoir operation [3]. In addition, the results of modern intelligence algorithms are not stable because of randomness in the iteration, which makes them difficult to use in practical application. In this case, some scholars have proposed methods such as allocating flood storage capacity [40] or aggregating multi-reservoirs into a virtual reservoir [41] to simplify the optimal operation model of multi-reservoir systems, but these methods do not reflect the hydraulic connection between reservoirs [3] or consider the impact of other water project operations on flood peaks in the downstream areas. Some scholars have proposed hybrid optimal algorithms to investigate the optimization of multiple reservoir operation. Wan et al. [42], based on PSO and POA, proposed a novel optimization method called Progressive Reservoir Algorithm-Particle Swarm Optimization (PRA-PSO) for multi-reservoir optimal operation in an inter-basin water transfer system. Zhang et al. [16] combined POA with DPSA, proposing an optimization method called modified Progressive Optimality Algorithm-Successive Approximation (POA-SA), and applied this method to tackle optimal flood control operation of a multi-reservoir system considering time-lags caused by Muskingum flood routing.
This study established an optimal flood control operation model of the complex flood control system, which considered hydraulic connection between reservoirs and operation of other water projects (e.g., navigation-power junctions and flood storage areas). The objective of the whole model was to minimize the flood peak at the flood control points. A hybrid algorithm of DP and POA called DP-POA [3] integrated with the simulation model was applied to optimize this flood control operation model. The middle and lower reaches of Ganjiang River, China, were selected as a case study.

2. Optimization Model and Methodology

2.1. Optimization Model of Complex Flood Control System

Take a complex flood control system composed of reservoirs, navigation-power junctions, flood storage areas, local flows, and flood control points for an example. Assume that there are I cascade water projects (including reservoirs and navigation-power junctions), J flood control points, and S flood storage areas, and the skeleton diagram of the sub flood control system composed of the ith water project, the jth flood control points, and the sth flood storage area is shown in Figure 1. The hydraulic connection during time period, t, is calculated as follows:
{ Q i ( t ) = ψ i ( q i ( t ) , V i ( t ) ) Q j ( t ) = φ j ( Q i ( t ) , Q i ( t 1 ) , , Q i ( t τ i ) ) + Δ q j ( t ) R s ( t ) = ρ s ( Q j ( t ) ) ,
where Qi(t), qi(t), Vi(t) are the discharge, inflow, and water storage of the ith water project during time period, t, respectively; τi is the time-lags of discharge from the ith water project flood routing to the jth flood control point; Qj(t) is the streamflow of the jth flood control point during time period, t; Δqj(t) is the local flow of the jth flood control point during time period, t; Rs(t) is the flood diversion volume of the sth flood storage area during time period, t; ψi(·) is the flood control operation rule of the ith water projects; φj(·) is the flood routing function between the ith water project and the jth flood control point; ρs(·) is the flood control operation rule of the sth flood storage area.
Figure 1 and Equation (1) show the complex hydraulic connection and operation of water projects in the sub flood control system. Comparing with a single reservoir flood control system, the complex flood control system needs to consider multi-reservoir joint operation, flood routing, and local flow of the flood control point; comparing with the multi-reservoir system, it needs to consider operation of other water projects besides reservoirs.

2.1.1. Objective Function

In this study, the objective function was to maximize the decrement rate of the flood peak at the flood control points. There were three flood control points needing to be protected by the flood control system, so this was a multi-objective optimization problem. Therefore, assigning weights for each flood control point could turn this into a single-objective optimization problem. The objective function was as follow:
f = max j = 1 N [ α j Q p j Q f j Q p j ]
where N is the number of flood control points; αi is the weight coefficient for the jth flood control point; Qpj is the flood peak at the jth flood control point without flood control operation; and Qfj is the flood peak at the jth flood control point using optimal flood control operation.

2.1.2. Constraints

The constraints were as follows:
(1)
Water balance equation:
[ q i ( t ) + q i ( t + 1 ) 2 Q i ( t ) + Q i ( t + 1 ) 2 ] Δ t = V i ( t + 1 ) V i ( t ) ,
where qi(t) and qi(t + 1) are the inflow of the ith reservoir during time period, t and t + 1, respectively; Qi(t) and Qi(t + 1) are the discharge of the ith reservoir during time period, t and t + 1, respectively; Δt is the time interval; and Vi(t) and Vi(t + 1) are the initial and final storage of the ith reservoir at time period, t.
(2)
Reservoir water level limits:
Z i min Z i ( t ) Z i max ,
where Zimin and Zimax are the minimal and maximal water level of the ith reservoir during the operation period, respectively; Zi(t) is the water level of the ith reservoir during the time period, t.
(3)
Boundary conditions:
Z i ( 1 ) = Z i , s t a r t , Z i ( T + 1 ) = Z i , e n d ,
where Zi,start is the water level of the ith reservoir at the initial time period; Zi,end is the water level of the ith reservoir at the final time period.
(4)
Reservoir discharge limits:
Q i min Q i ( t ) Q i max [ Z i ( t ) ] ,
where Qimax[Zi(t)] is the maximal discharge of the ith reservoir when the water level of the ith reservoir is Zi(t) during time period, t. Qimin is the minimum discharge of the ith reservoir.
(5)
Discharge variation limits:
| Q i ( t + 1 ) Q i ( t ) | Δ Q i ,
where ΔQi is the allowable discharge variation of the ith reservoir between adjacent time periods.
(6)
Flood routing constraints using the Muskingum method:
Q h d m ( t + 1 ) = C 0 m I h d m ( t + 1 ) + C 1 m I h d m ( t ) + C 2 m Q h d m ( t ) + Q q m ( t ) ,
where Qhdm(t) is the streamflow of the down section of the mth sub-river during the time period, t; Ihdm(t) is the streamflow of the upper section of the mth river during the time period, t; Qqm(t) is the local flow of the mth sub-river; C0m, C1m, C2m means the Muskingum method coefficients of the mth sub-river, respectively.
(7)
Other water projects’ constraints:
Other water projects such as navigation-power junctions and flood storage areas, affect flood routing by flood detention. This study simulates their operation according to their flood control rules for these constraints.
(a)
Constraints of navigation-power junctions:
Navigation-power junctions do not function as flood control for the lack of flood control storage capacity. However, to reduce the submergence loss in the reservoir area, the navigation-power junctions are required to lower the water level and open all the discharge gates when encountering large flood events. The flood control operation rules of navigation-power junctions can be generalized into four stages as follows:
(I)
Power generation stage
When the inflow is small, the main purpose of the navigation-power junction is power generation, and discharge is equal to the inflow. This can be calculated as:
q n j , k ( t ) = Q n j , k ( t ) , Z n j , k ( t ) = Z n l , k
where Qnj,k(t) and qnj,k(t) are inflow and discharge of the kth navigation-power junction during time period, t, respectively; Znj,k(t) is the water level of the kth navigation-power junction during time period, t; Znl,k is the normal pool level of the kth navigation-power junction.
(II)
Water level reduction stage
When the inflow continues increasing, the outflow is increased gradually, which can be calculated as:
q n j , k ( t ) = min { Q n j , k ( t ) + Δ Q d , k , q n j , k max } ,
where ΔQd,k is discharge increment of the kth navigation-power junction and qnj,kmax is the maximal discharge volume of the kth navigation-power junction when all the discharge gates are opened, which can be calculated as:
{ V n j , k ( t + 1 ) = V n j , k ( t ) + [ Q n j , k ( t ) + Q n j , k ( t + 1 ) 2 q n j , k ( t ) + q n j , k ( t + 1 ) 2 ] Δ t Z n j , k ( t + 1 ) = φ n j , k [ V n j , k ( t + 1 ) ] Z n j , k = max { Z n j , k ( t + 1 ) , Z n j , k min } q n j , k max = ψ n j , k ( Z n j , k ) ,
where Znj,kmin is the minimum water level of the kth navigation-power junction during the flood control operation period; φnj,k(·) is the relationship curve of water level and reservoir capacity of the kth navigation-power junction;ψnj,k(·) is the relationship curve of water level and discharge of the kth navigation-power junction, which is similar to the relationship of water level and streamflow in the natural river channel.
(III)
All discharge gates opened stage
When the inflow of the navigation-power junction is too large, all the discharge gates are opened, and the water level is close to that of the natural river, which can be calculated as:
q n j , k ( t ) = Q n j , k ( t ) , Z n j , k ( t ) = f n j , k [ q n j , k ( t ) ]
where fnj,k[qnj,k(t)] is the relationship of water level and streamflow in the natural river channel.
(IV)
Water storage stage
In the falling limb, the water level will rise to the normal pool level, which can be calculated as:
q n j , k ( t ) = Q n j , k ( t ) Δ Q r , k ,
where ΔQr,k is the discharge reduction of the kth navigation-power junction.
(b)
Constraints of the flood storage area:
The flood control operation of flood storage areas can be generalized as follows:
Q d , s ( t ) = Q u , s ( t ) ρ s ( Q u , s ( t ) ) ,
ρ s ( Q u , s ( t ) ) = { 0 Q u , s ( t ) L   or   W ( t ) > W max Q u , s ( t ) L L < Q u , s ( t ) Q max   and   W ( t ) W max Q max Q u , s ( t ) > Q max   and   W ( t ) W max ,
W ( t ) = t = 1 T 1 [ Q u , s ( t ) + Q u , s ( t + 1 ) 2 Q d , s ( t ) + Q d , s ( t + 1 ) 2 ] Δ t ,
where Qd,s(t) and Qu,s(t) are streamflow of the sth flood storage area before and after flood diversion during time period, t; ρs(Qu,s(t)) means flood diversion streamflow according to flood control rules of the sth flood storage area during time period, t; L is the indicative flow for the use of the sth flood storage area; W(t) is the water storage volume of the sth flood storage area during time period, t; Qmax and Wmax are the maximum flood diversion discharge and maximum water storage volume, respectively.
(8)
Non-negative constraints:
All the variables are non-negative.

2.2. Optimization Methods

DP is a very efficient optimization method to minimize the discharge peak of a single reservoir system. However, considering difficulties such as the “curse of dimensionality”, time-lag issues, other water projects operation, and hydraulic connection, DP cannot be used to optimize flood control operation of this complex system. POA can alleviate these difficulties, but its final optimal results are influenced by the initial solutions [3]. In this case, DP-POA, a hybrid optimization algorithm of DP and POA [3] combined with the simulation model was applied to find the optimal operation polices for the complex flood control system, in which DP is used first to obtain the initial solutions for POA. The processes of DP-POA are as follows:
Step 1: input the data of reservoirs, flood hydrographs, flood control operation rules of other water projects, and so on. Set the total iteration number = K.
Step 2: decompose the optimal cascade-reservoirs flood control operation into optimal multiple single-reservoir flood operation and optimize the each single-reservoir flood control operation using DP. When using DP, the sub-objective function for each reservoir can be as follows:
F i = min t = 1 T Q i 2 ( t ) ,
where Qi(t) is the release of the ith reservoir during the time period, t; T = the number of the time periods. Constraints are similar to Equations (3)–(7).
Step 3: based on DP, the initial solution u = {Z1(1),…, Zi(t + 1),…, Zi(T + 1),…, ZM(1),…, ZM(T + 1)} is obtained. Here, M is the total number of reservoirs. Set iteration k = 1 and the reservoir number i = 1.
Step 4: using POA to update the initial status of the ith reservoir from time period, t = 2 to T, according to Equation (2). The reservoir number i = i + 1.
Step 5: if i < M, then go to Step 4 until i = M; else, obtain the updated optimal solution u* = {Z1*(1),…, Zi*(t + 1),…, Zi*(T + 1),…, ZM*(1),…, ZM*(T + 1)}. Replace the initial solution u with u*. Then the iteration k = k + 1.
Step 6: if the iteration k < K, then go to Step 4; else output the u*, and the optimization is completed.
For Steps 3 and 4, because the water level of each reservoir during each time period, t, is certain, this means the discharge of each reservoir can be calculated according to Equation (3) and flood hydrographs of flood control points can be also simulated using Equation (8).
The procedure of DP-POA is shown in the Figure 2.

3. Study Area and Data

3.1. Study Area

The Ganjiang River located in the south part of the middle reaches of the Yangtze River is the largest river of the Poyang Lake basin, with a length of 766 km and basin area of 83,500 km2 (Figure 3). This basin is situated in the subtropical humid monsoon climate zone with mean annual discharge of 2130 m3/s and mean annual precipitation of approximately 1400–1800 mm [43]. Runoff in the flood season from April to June, accounts for more than half of the annual runoff [44]. The Ganjiang River has a great influence on the social and economic development of Jiangxi Province, China. However, flood disasters frequently happen in the middle and lower reaches of the Ganjiang River, which harm human lives and property and threaten the safety of society. In addition, more and more water projects have been built in the middle and lower reaches of Ganjiang River, making the hydraulic connection complex and it is worth researching on optimal flood control operation to reduce the loss of flood disasters.
The middle and lower reaches of Ganjiang River flood control system is composed of two-stage cascade reservoirs, four navigation-power junctions, a flood storage area, five local inflows, and three flood control points. The two-stage cascade reservoirs, consisting of Wan’an reservoir and Xiajiang Reservoir, play a major role in the flood control in the middle and lower reaches of the Ganjiang River. Wan’an reservoir, until now, is the largest reservoir in the Jiangxi Province, China, with a total storage of 22.16 × 108 m3 [45]. Xiajiang reservoir is the second largest reservoir in the Jiangxi province with a storage capacity of 11.87 × 108 m3. Quangang flood storage area, with a storage capacity of 8.35 × 108 m3, located in the lower reaches of Ganjiang River, helps protect flood control points Shishang and Waizhou through flood diversion. The two navigation-power junctions, Jinggangshan and Shihutang, are situated between Wan’an reservoir and Ji’an flood control point. The other navigation-power junctions, Xingan and Longtoushan, are located between the Xiajiang reservoir and Shishang flood control point. The skeleton of the middle and lower reaches of Ganjiang River flood control system is shown in Figure 4. The characteristic parameters of these reservoirs and navigation-power junctions are listed in Table 1.

3.2. Input Data

The middle and lower reaches of the Ganjiang River are divided into eight stream segments by reservoirs, navigation-power junctions, and flood control points. The study used the Muskingum method [46] for flood routing between the stream segments. Different stream segments have different Muskingum coefficients, which are calibrated by the historical flood data and provided by the Jiangxi Provincial Institute of Water Sciences. The time interval, Δt, is 6 h. These Muskingum coefficients are listed in Table 2.
The conventional flood control operation rules of the Wan’an and Xiajiang reservoirs are listed in Table 3 and Table 4.
Six flood hydrographs (namely those from 1961, 1968, 1973, 1994, 1998, and 2010) were selected as typical flood hydrographs. Then all the typical flood hydrographs were magnified to design flood hydrographs with 20-year, 50-year, 100-year, 200-year, 300-year, and 500-year return periods using the peak ratio of the Waizhou station.
All data including flood control operation rules of water projects, the reservoir characteristics, flood hydrographs, and Muskingum method coefficients mentioned were provided by the Jiangxi Provincial Institute of Water Sciences and were used as input data for the optimal flood control model.

4. Results

This study applied DP-POA to optimize flood control operation and compared the results with the results from conventional flood control operations. There are three flood control points in the complex flood control system. In this study, the significance of these flood control points was assumed to be the same, so the weighting coefficient for each flood control point was equal, which is 1/3. The number of discrete points between the maximal and minimal water level of Wan’an and Xiajiang reservoirs was 100. The water level of Wan’an and Xiajiang reservoirs at the initial and final time period was equal to their flood limited water level. The maximal and minimal water level during the optimal flood control operation was equal to those of the conventional operation.

4.1. Flood Peak Reduction of Flood Control Points

Table 5 and Table 6 list the values (average of six years) of maximal water level and flood peak reduction of Wan’an–Xiajiang cascade reservoirs during the operation period using conventional flood control operation and DP-POA, respectively. Table 7, Table 8 and Table 9 list the flood peak reduction (average of six years) at the Ji’an, Shishang, and Waizhou flood control points using conventional flood control operation and DP-POA, respectively.
Table 5 and Table 6 show that the maximal water level of the Xiajiang reservoir is equal to its flood limited water level for the typical and the designed 20-year return period flood. This is because the flood peak of the Xiajiang reservoir is less than 20,000 m3/s and its discharge is equal to inflow. For the same designed return period floods, the flood peak reduction of the Wan’an reservoir using DP-POA was larger than that using conventional flood control operation, which means optimization method can make full use of flood control capacity to reduce flood peak. Table 5 shows that the flood peak reduction of the Xiajiang reservoir increases as designed return periods increase from 50 years to 200 years. While for designed floods with 300-year and 500-year return periods, the flood peak reduction of the Xiajiang reservoir decreased from 1819 to 1339 m3/s. This is because flood peak reduction listed in Table 5 and Table 6 is the average value of six years. When flood peak is larger, greater storage of the reservoir is used to detain the inflow without exceeding the acceptable safety level of the reservoir. However, the flood peaks of some designed floods were too large for the Xiajiang reservoir and the discharge was equal to the inflow for the safety of the reservoir, which means that flood peak reduction was zero. For example, the flood peak reduction of the Xiajiang reservoir was zero only for “1961” designed floods with 300-year return periods, while for designed floods with 500-year return periods, the flood peak reduction of “1961”, “1968”, and ”1994” was zero. Therefore, the average flood peak reduction of the Xiajiang reservoir for designed floods with 500-year return periods was less than that for designed floods with 300-year return periods.
Table 7, Table 8 and Table 9 show that the DP-POA can obtain larger flood peak reduction of flood control points comparing with conventional flood control operation for all the designed return period floods. The Ji’an flood control point is located between the Wan’an and Xiajiang reservoir and its flood peak is mainly influenced by the operation of upstream Wan’an reservoir. The flood peak reduction of Ji’an flood control point using the conventional operation ranges from 365 to 1500 m3/s while it increases by 1041~1498 m3/s using DP-POA. The Shishang flood control point is located in the downstream of Xiajiang reservoir and its flood peak is mainly reduced by the operation of the Wan’an-Xiajiang cascade reservoirs and Quangang flood storage area. Comparing with the conventional operation, DP-POA can increase flood peak reduction of Shishang flood control point by 584~1151 m3/s. The Waizhou flood control point is the outlet of the Ganjiang River. DP-POA can decrease its flood peak by 1080~4809 m3/s while the conventional operation can decrease the flood peak by 747~4063 m3/s. Considering that the maximal water level using DP-POA is equal to that using the conventional operation, these results indicate that DP-POA is more efficient in making full use of the flood control capacity of reservoirs to reduce the flood peak of flood control points.
Figure 5, Figure 6 and Figure 7 show the 100-year return period flood peaks of Ji’an, Shishang, and Waizhou flood control points using conventional and optimal flood control operation for different typical years. Because of utilization of the Quangang flood storage area, the “1994” flood peak of the Shishang flood control point was 24,000 m3/s whenever using conventional or optimal operation. Except for that, other flood peaks of three flood control points using DP-POA were smaller than those using the conventional operation.

4.2. Flood Diversion Volume of the Flood Storage Area

Figure 8 shows the flood diversion volume (average of six years) of typical and designed floods using the conventional operation and DP-POA at the Quangang flood storage area.
Figure 8 shows that the flood diversion volume using the conventional operation and DP-POA increased as the designed flood peak (the return period) increased. For typical floods and designed floods with 20-year and 50-year return periods, the flood diversion volume using the conventional operation was zero. For designed floods with 100-year, 200-year, 300-year, and 500-year return periods, the flood diversion volume using the conventional operation ranges from 0.04 × 108 to 2.49 × 108 m3. Comparing with the conventional operation, the flood diversion volume using DP-POA was smaller. For typical floods and designed floods with 20-year, 50-year, and 100-year return periods, the flood diversion volume of the optimal operation was zero. For designed floods with 200-year, 300-year, and 500-year return periods, the flood diversion volume turned from 0.33 × 108 m3 to 1.79 × 108 m3. This means that DP-POA can retain more storage capacity and it is safe for operation in the Quangang flood storage area.

4.3. Flood Control Operation Process of Reservoirs

This study selected the “1961” flood with a 100-year return period for the further analysis of the conventional and optimal flood control operation for the Wan’an–Xiajiang cascade reservoirs. Figure 9 and Figure 10 show conventional and optimal flood control operation of Wan’an–Xiajiang cascade reservoirs.
Figure 9a and Figure 10 show that there are three peaks in the flood hydrograph of the Wan’an reservoir. The maximal flood peak occurs during the 51st time period. As shown in Figure 9a, before the occurrence of the maximal flood peak, the water level of the Wan’an reservoir using the conventional operation rose gradually from 85.0 to 87.5 m, which means the partial flood control capacity was occupied by the first two flood peaks. As shown in Figure 10a, DP-POA can reduce the first two flood peaks greatly compared with the conventional operation and decrease the water level of the Wan’an reservoir from 92.48 to 85.0 m before the occurrence of the maximal flood peak, which indicates that greater flood control capacity is retained to reduce flow from the 48th to 52nd time period.
Figure 9b and Figure 10b show that the inflow of the Xiajiang reservoir fluctuates before the occurrence of the maximal flood peak. The maximal inflow peak of the Xiajiang reservoir using the conventional operation was 22,344 m3/s, which occurs at the 53rd time period. However because DP-POA optimizes the discharge of the Wan’an reservoir, the maximal inflow peak of the Xiajiang reservoir using DP-POA was reduced to 21,060 m3/s, which is even smaller than the discharge of the Xiajiang reservoir using the conventional operation. In addition, DP-POA did not reduce the inflow peak of the Xiajiang reservoir. Comparing the occurrence time of inflow peak with that of discharge peak, it is shown that DP-POA delays the occurrence time of flood peak from the 46th time period to the 50th time period. This indicates that DP-POA can consider the downstream local flow and delay occurrence time of the discharge peak to avoid encountering downstream flood peak.

5. Discussions and Limitations

The optimal flood control operation model and optimization method is deterministic without random factors. Therefore, the optimal flood control operation results of typical and designed floods were stable and could give some guidance on improving conventional flood control operation rules of the Wan’an–Xiajiang cascade reservoirs. When applied in practical flood control operation, if the flood forecast is accurate and the flood forecast period is long enough, the flood hydrograph of the Wan’an and Xiajiang reservoir and all the local flows can be determined. DP-POA can be applied to find optimal flood control operation polices. For example, if the flood forecast period is 10 days, assuming the computing time interval is 6 h, then the flood hydrograph is determined, and the number of the streamflow value is 40. Then DP-POA can be used to optimize discharge of the Wan’an–Xiajiang cascade reservoirs and the water level during each time period is obtained. Considering that DP-POA in this study is programmed with C++ language, it does not take too much time to optimize the flood control operation. Then reservoirs operators can adjust the opening of spillway gates according to the optimal discharge polices in time. However, many uncertainties including flood forecast error and the uncertainties in the discharge capacity curve and water level-storage curve have not been considered, and stochastic modeling, which may better address these uncertainties, will be researched in future study.
This study turned the multi-objective optimization problem into a single-objective optimization problem by assigning the same weighting coefficients for each flood control point turning. It is obvious that the flood peak reduction of the flood control point will be larger if its assigned weight is larger. The sensitive analysis of weighting coefficients and relationship of flood peak reduction of these flood control points will be further researched in future. In addition, it is feasible to redo this study work in future if some external influences such as climate change and human activities can be quantified. The input data such as flood hydrographs, flood routing parameters, characteristics parameters of reservoirs can be recalibrated considering climate change and human activity. Then the similar deterministic flood control model can be established, and DP-POA can be used to find optimal results. If these influences are uncertain or hard to determine, stochastic frameworks and flood risk assessment models are better for the uncertain problems. Moreover, there are two kinds of flood control operation rules in practical flood control operation. One is the conventional flood control operation rules, which are obtained by flood regulating calculation for historical floods and designed flood with different return periods. The other is real-time flood control operation rules, which are obtained according to flood forecast information. Adjusting the given operation rules on an hourly or perhaps daily basis concerns real-time and flood forecast information. This is not considered in this study, and it will be researched in future.

6. Conclusions

Based on the objective of minimizing the maximum streamflow of downstream flood control points, this study established an optimal flood control operation model for a complex flood control system composed of cascade reservoirs, navigation-power junctions, flood storage areas, local flow, and flood control points. An optimal method called DP-POA was used to optimize flood control operation of reservoirs and the Muskingum method was applied for flood routing. A complex flood control system in the middle and lower reaches of the Ganjiang River was selected as a case study, and the conclusions could be drawn as follows:
(1) The optimal flood control operation problem with time-lag issues, multistage and multidimensional decision variables can be effectively solved by DP-POA. DP-POA outperforms the conventional flood control operation in reducing the flood peaks of the downstream flood control points and flood diversion volume of the flood storage area.
(2) For the typical floods and designed floods with 20-year return period, because the flood peak of the Xiajiang reservoir is less than 20,000 m3/s, its discharge is equal to its inflow and the flood diversion volume of the Quangang flood storage area is zero. This means flood peaks of downstream flood control points are mainly reduced by operation of the Wan’an reservoir. For larger floods, such as designed floods with 50-year, 100-year, 200-year, 300-year, and 500-year return periods, the joint flood control operation of the Wan’an-Xiajiang cascade reservoirs and Quangang flood storage area can greatly reduce flood peak of downstream flood control points.
(3) The flood peak reduction of Ji’an, Shishang, and Waizhou flood control points using optimal flood control operation is 1491~2541, 1102~5359, and 1080~4809 m3/s, respectively for typical floods and designed floods with 20-year, 50-year, 100-year, 200-year, 300-year, and 500-year return periods, which is a great progress comparing with the conventional operation. Considering that the maximal water level of Wan’an–Xiajiang cascade reservoirs during the periods of conventional and optimal operation is the same, this means DP-POA was more efficient in making full of flood control storage of reservoirs.
(4) Comparing with the conventional operation, the flood diversion volume of the Quangang flood storage area using DP-POA was smaller for the same designed floods. The maximum flood diversion volume using the conventional operation was 2.49 × 108 m3, while the corresponding value of DP-POA was only 1.79 × 108 m3. The greater the storage retained, the safer it is for the operation of the flood storage area.
(5) For flood hydrographs with more than one flood peak, DP-POA can reduce the water level of reservoirs before the occurrence of the maximal flood peak, retaining more flood control capacity to help reduce flood peak of downstream flood control points.

Author Contributions

Conceptualization: D.Z., Y.M., and X.X.; methodology: D.Z.; formal analysis: D.Z.; investigation: D.Z.; resources: X.X.; writing—original draft preparation, D.Z.; writing—review and editing: Y.M., J.C., and Y.B.; funding acquisition, Y.M. All the authors have approved of the submission of this manuscript.

Funding

The research was financially founded by the National Science Foundation of China (no. 51479140), the Key Research and Development Plan of Jiangxi Province (no. 20181ACG70018), and the Science and Technology Project of Department of Water Resources of Jiangxi Province (no. 202023ZDKT02).

Acknowledgments

The authors thank the Jiangxi Provincial Institute of Water Sciences for providing the data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Valeriano, O.C.S.; Koike, T.; Yang, K.; Yang, D.W. Optimal Dam Operation during Flood Season Using a Distributed Hydrological Model and a Heuristic Algorithm. J. Hydrol. Eng. 2010, 15, 580–586. [Google Scholar] [CrossRef] [Green Version]
  2. Peng, Y.; Chen, K.; Yan, H.G.; Yu, X.L. Improving Flood-Risk Analysis for Confluence Flooding Control Downstream Using Copula Monte Carlo Method. J. Hydrol. Eng. 2017, 22, 04017018. [Google Scholar] [CrossRef]
  3. Zhou, C.; Sun, N.; Chen, L.; Ding, Y.; Zhou, J.Z.; Zha, G.; Luo, G.L.; Dai, L.; Yang, X. Optimal Operation of Cascade Reservoirs for Flood Control of Multiple Areas Downstream: A Case Study in the Upper Yangtze River Basin. Water 2018, 10, 1250. [Google Scholar] [CrossRef] [Green Version]
  4. Simonovic, S. The implicit stochastic model for reservoir yield optimization. Water Resour. Res. 1987, 23, 2159–2165. [Google Scholar] [CrossRef]
  5. Labadie, J.W. Optimal Operation of Multireservoir Systems: State-of-the-Art Review. J. Water Resour. Plan. Manag. 2004, 130, 93–111. [Google Scholar] [CrossRef]
  6. Hui, R.; Lund, J.R. Flood Storage Allocation Rules for Parallel Reservoirs. J. Water Resour. Plan. Manag. 2015, 141, 04014075. [Google Scholar] [CrossRef] [Green Version]
  7. Chen, C.; Yuan, Y.B.; Yuan, X.H. An Improved NSGA-III Algorithm for Reservoir Flood Control Operation. Water Resour. Manag. 2017, 31, 4469–4483. [Google Scholar] [CrossRef]
  8. Hsu, N.-S.; Wei, C.-C. A multipurpose reservoir real-time operation model for flood control during typhoon invasion. J. Hydrol. 2007, 336, 282–293. [Google Scholar] [CrossRef]
  9. Li, X.A.; Guo, S.L.; Liu, P.; Chen, G.Y. Dynamic control of flood limited water level for reservoir operation by considering inflow uncertainty. J. Hydrol. 2010, 391, 124–132. [Google Scholar] [CrossRef]
  10. Zhou, Y.L.; Guo, S.L.; Liu, P.; Xu, C.Y. Joint operation and dynamic control of flood limiting water levels for mixed cascade reservoir systems. J. Hydrol. 2014, 519, 248–257. [Google Scholar] [CrossRef]
  11. Liu, P.; Li, L.P.; Guo, S.L.; Xiong, L.H.; Zhang, W.; Zhang, J.W.; Xu, C.-Y. Optimal design of seasonal flood limited water levels and its application for the Three Gorges Reservoir. J. Hydrol. 2015, 527, 1045–1053. [Google Scholar] [CrossRef]
  12. Gong, Y.; Liu, P.; Cheng, L.; Chen, G.Y.; Zhou, Y.L.; Zhang, X.Q.; Xu, W.F. Determining dynamic water level control boundaries for a multi-reservoir system during flood seasons with considering channel storage. J. Flood Risk Manag. 2019, 13, e12586. [Google Scholar] [CrossRef]
  13. Dai, L.; Zhou, J.Z.; Chen, L.; Huang, K.D.; Wang, Q.S.; Zha, G. Flood-risk analysis based on a stochastic differential equation method. J. Flood Risk Manag. 2018, 12, 10. [Google Scholar] [CrossRef] [Green Version]
  14. Huang, K.D.; Ye, L.; Chen, L.; Wang, Q.S.; Dai, L.; Zhou, J.Z.; Singh, V.P.; Huang, M.T.; Zhang, J.H. Risk analysis of flood control reservoir operation considering multiple uncertainties. J. Hydrol. 2018, 565, 672–684. [Google Scholar] [CrossRef]
  15. Meng, X.J.; ChWang, J.X.-W.; Wang, X.B.; Wang, Y.M.; Wang, Z.Z. Flood control operation coupled with risk assessment for cascade reservoirs. J. Hydrol. 2019, 572, 543–555. [Google Scholar] [CrossRef]
  16. Zhang, W.; Liu, P.; Chen, X.Z.; Wang, L.; Ai, X.S.; Feng, M.Y.; Liu, D.D.; Liu, Y.Y. Optimal Operation of Multi-reservoir Systems Considering Time-lags of Flood Routing. Water Resour. Manag. 2015, 30, 523–540. [Google Scholar] [CrossRef]
  17. Li, X.G.; Wei, X. An Improved Genetic Algorithm-Simulated Annealing Hybrid Algorithm for the Optimization of Multiple Reservoirs. Water Resour. Manag. 2007, 22, 1031–1049. [Google Scholar] [CrossRef]
  18. Chang, F.-J.; Chen, L. Real-Coded Genetic Algorithm for Rule-Based Flood Control Reservoir Management. Water Resour. Manag. 1998, 12, 185–198. [Google Scholar] [CrossRef]
  19. Kalteh, A.M. Wavelet Genetic Algorithm-Support Vector Regression (Wavelet GA-SVR) for Monthly Flow Forecasting. Water Resour. Manag. 2014, 29, 1283–1293. [Google Scholar] [CrossRef]
  20. Afshar, M.H. Extension of the constrained particle swarm optimization algorithm to optimal operation of multi-reservoirs system. Int. J. Electr. Power Energy Syst. 2013, 51, 71–81. [Google Scholar] [CrossRef]
  21. Ostadrahimi, L.; Mariño, M.A.; Afshar, A. Multi-reservoir Operation Rules: Multi-swarm PSO-based Optimization Approach. Water Resour. Manag. 2011, 26, 407–427. [Google Scholar] [CrossRef]
  22. Chang, J.X.; Bai, T.; Huang, Q.; Yang, D.-W. Optimization of Water Resources Utilization by PSO-GA. Water Resour. Manag. 2013, 27, 3525–3540. [Google Scholar] [CrossRef]
  23. Luo, J.G.; Chen, C.; Xie, J.C. Multi-objective Immune Algorithm with Preference-Based Selection for Reservoir Flood Control Operation. Water Resour. Manag. 2014, 29, 1447–1466. [Google Scholar] [CrossRef]
  24. Fu, X.; Li, A.-Q.; Wang, L.P.; Ji, C.M. Short-term scheduling of cascade reservoirs using an immune algorithm-based particle swarm optimization. Comput. Math. Appl. 2011, 62, 2463–2471. [Google Scholar] [CrossRef] [Green Version]
  25. Bozorg-Haddad, O.; Janbaz, M.; Loáiciga, H.A. Application of the gravity search algorithm to multi-reservoir operation optimization. Adv. Water Resour. 2016, 98, 173–185. [Google Scholar] [CrossRef]
  26. Qaderi, K.; Akbarifard, S.; Madadi, M.R.; Bakhtiari, B. Optimal operation of multi-reservoirs by water cycle algorithm. Proc. Inst. Civ. Eng. Water Manag. 2018, 171, 179–190. [Google Scholar] [CrossRef]
  27. Xu, Y.N.; Mei, Y.D. A modified water cycle algorithm for long-term multi-reservoir optimization. Appl. Soft Comput. 2018, 71, 317–332. [Google Scholar] [CrossRef]
  28. Crawley, P.D.; Dandy, G.C. Optimal Operation of Multiple?Reservoir System. J. Water Resour. Plan. Manag. 1993, 119, 1–17. [Google Scholar] [CrossRef]
  29. De Barros, M.T.L.; Tsai, F.T.-C.; Yang, S.-L.; Lopes, J.E.G.; Yeh, W.W.-G. Optimization of Large-Scale Hydropower System Operations. J. Water Resour. Plan. Manag. 2003, 129, 178–188. [Google Scholar] [CrossRef]
  30. Wei, C.-C.; Hsu, N.-S. Multireservoir real-time operations for flood control using balanced water level index method. J. Environ. Manag. 2008, 88, 1624–1639. [Google Scholar] [CrossRef] [PubMed]
  31. Yakowitz, S. Dynamic programming applications in water resources. Water Resour. Res. 1982, 18, 673–696. [Google Scholar] [CrossRef]
  32. Kelman, J.; Stedinger, J.R.; Cooper, L.A.; Hsu, E.; Yuan, S.-Q. Sampling stochastic dynamic programming applied to reservoir operation. Water Resour. Res. 1990, 26, 447–454. [Google Scholar] [CrossRef]
  33. Zhang, J.W.; Liu, P.; Wang, H.; Lei, X.H.; Zhou, Y.L. A Bayesian model averaging method for the derivation of reservoir operating rules. J. Hydrol. 2015, 528, 276–285. [Google Scholar] [CrossRef]
  34. Cheng, C.T.; Wang, S.; CAhmaudi, K.M.W.H.; Wu, X.Y. Parallel discrete differential dynamic programming for multireservoir operation. Environ. Model. Softw. 2014, 57, 152–164. [Google Scholar] [CrossRef] [Green Version]
  35. Heidari, M.; Chow, V.T.; Kokotovic, P.V.; Meredith, D.D. Discrete Differential Dynamic Programing Approach to Water Resources Systems Optimization. Water Resour. Res. 1971, 7, 273–282. [Google Scholar] [CrossRef]
  36. Opan, M. Irrigation-energy management using a DPSA-based optimization model in the Ceyhan Basin of Turkey. J. Hydrol. 2010, 385, 353–360. [Google Scholar] [CrossRef]
  37. Liu, P.; Cai, X.M.; Guo, S.L. Deriving multiple near-optimal solutions to deterministic reservoir operation problems. Water Resour. Res. 2011, 47, 20. [Google Scholar] [CrossRef]
  38. Bai, T.; Chang, J.X.; Chang, F.L.J.-C.; Huang, Q.; Wang, Y.M.; Chen, G.-S. Synergistic gains from the multi-objective optimal operation of cascade reservoirs in the Upper Yellow River basin. J. Hydrol. 2015, 523, 758–767. [Google Scholar] [CrossRef]
  39. Jiang, Z.Q.; Ji, C.M.; Qin, H.; Feng, Z.-K. Multi-stage progressive optimality algorithm and its application in energy storage operation chart optimization of cascade reservoirs. Energy 2018, 148, 309–323. [Google Scholar] [CrossRef]
  40. Fu, X.; Li, A.-Q.; Wang, H. Allocation of Flood Control Capacity for a Multireservoir System Located at the Yangtze River Basin. Water Resour. Manag. 2014, 28, 4823–4834. [Google Scholar] [CrossRef]
  41. Zhou, Y.L.; Guo, S.L.; Chang, F.L.J.-C.; Liu, P.; Chen, A.B. Methodology that improves water utilization and hydropower generation without increasing flood risk in mega cascade reservoirs. Energy 2018, 143, 785–796. [Google Scholar] [CrossRef]
  42. Wan, W.H.; Guo, X.N.; Lei, X.H.; Jiang, Y.Z.; Wang, H. A Novel Optimization Method for Multi-Reservoir Operation Policy Derivation in Complex Inter-Basin Water Transfer System. Water Resour. Manag. 2017, 32, 31–51. [Google Scholar] [CrossRef]
  43. Ye, X.C.; Zhang, Q.; Liu, J.; Li, X.H.; Xu, C.-Y. Distinguishing the relative impacts of climate change and human activities on variation of streamflow in the Poyang Lake catchment, China. J. Hydrol. 2013, 494, 83–95. [Google Scholar] [CrossRef]
  44. Guo, L.-P.; Mu, X.-M.; Hu, J.-M.; Gao, P.; Zhang, Y.-F.; Liao, K.-T.; Bai, H.; Chen, X.-L.; Song, Y.-J.; Jin, N.; et al. Assessing Impacts of Climate Change and Human Activities on Streamflow and Sediment Discharge in the Ganjiang River Basin (1964–2013). Water 2019, 11, 1679. [Google Scholar] [CrossRef] [Green Version]
  45. Wen, T.F.; Xiong, L.H.; Jiang, C.; Hu, J.M.; Liu, Z.J. Effects of Climate Variability and Human Activities on Suspended Sediment Load in the Ganjiang River Basin, China. J. Hydrol. Eng. 2019, 24, 05019029. [Google Scholar] [CrossRef]
  46. Gill, M.A. Flood routing by the Muskingum method. J. Hydrol. 1978, 36, 353–363. [Google Scholar] [CrossRef]
Figure 1. The skeleton diagram of the sub flood control system.
Figure 1. The skeleton diagram of the sub flood control system.
Water 12 01883 g001
Figure 2. The procedure of Dynamic Programming-Progressive Optimality Algorithm (DP-POA).
Figure 2. The procedure of Dynamic Programming-Progressive Optimality Algorithm (DP-POA).
Water 12 01883 g002
Figure 3. The location of reservoirs, navigation-power junctions, and flood control points in the Ganjiang River.
Figure 3. The location of reservoirs, navigation-power junctions, and flood control points in the Ganjiang River.
Water 12 01883 g003
Figure 4. The skeleton of the middle and lower reaches of Ganjiang River basin flood control system.
Figure 4. The skeleton of the middle and lower reaches of Ganjiang River basin flood control system.
Water 12 01883 g004
Figure 5. Flood peaks of the Ji’an flood control point with 100-year return period.
Figure 5. Flood peaks of the Ji’an flood control point with 100-year return period.
Water 12 01883 g005
Figure 6. Flood peaks of the Shishang flood control point with 100-year return period.
Figure 6. Flood peaks of the Shishang flood control point with 100-year return period.
Water 12 01883 g006
Figure 7. Flood peaks of the Waizhou flood control point with 100-year return period.
Figure 7. Flood peaks of the Waizhou flood control point with 100-year return period.
Water 12 01883 g007
Figure 8. The average flood diversion volume of the Quangang flood storage for typical and design floods.
Figure 8. The average flood diversion volume of the Quangang flood storage for typical and design floods.
Water 12 01883 g008
Figure 9. Conventional flood control operation of the Wan’an–Xiajiang cascade reservoirs: (a) discharge of the Wan’an reservoir; (b) discharge of the Xiajiang reservoir.
Figure 9. Conventional flood control operation of the Wan’an–Xiajiang cascade reservoirs: (a) discharge of the Wan’an reservoir; (b) discharge of the Xiajiang reservoir.
Water 12 01883 g009
Figure 10. Optimal flood control operation of the Wan’an–Xiajiang cascade reservoirs: (a) discharge of the Wan’an reservoir; (b) discharge of the Xiajiang reservoir.
Figure 10. Optimal flood control operation of the Wan’an–Xiajiang cascade reservoirs: (a) discharge of the Wan’an reservoir; (b) discharge of the Xiajiang reservoir.
Water 12 01883 g010
Table 1. Characteristic parameters of reservoirs and navigation-power junctions.
Table 1. Characteristic parameters of reservoirs and navigation-power junctions.
ReservoirWan’anXiajiangJinggangshanShishutangXinganLongtoushan
Basin area (km2)36,90062,71040,93743,77064,77672,810
Mean annual runoff (m3/s)94716401050116016901890
Dead water level (m)854467.756.23223.7
Flood limited water level (m)8545----
Flood control highwater level (m)93.649----
Normal pool level (m)96466856.532.524.2
Regulation abilitySeasonalSeasonalDailyDailyDailyDaily
Table 2. The Muskingum method coefficients for flood routing.
Table 2. The Muskingum method coefficients for flood routing.
Stream SegmentMuskingum CoefficientsSub-River NumberΔt/h
C0C1C2
10.23080.53850.230816
20.28570.42860.285716
30.28570.42860.285716
40.28570.42860.285716
50.33330.33330.333316
60.33330.33330.333316
70.33330.33330.333326
80.500000.500046
Table 3. Conventional flood control operation of the Wan’an reservoir.
Table 3. Conventional flood control operation of the Wan’an reservoir.
ConditionInflow (m3/s)Reservoir Water Level (m)Discharge (m3/s)
Rising limbq ≤ 880085 < Z ≤ 93.6Q = q
8800 < q ≤ 9550Q = 8800
9550 < q ≤ 12,000Q = q − 750
q > 12,000Q = q − 4000
-Z > 93.6Open discharge
Falling limbq > 8000Z > 88Open discharge
Z < 88Q = q
4000 < q ≤ 8000Z > 88Q = 8800
Z < 88Q = q
2800 < q ≤ 4000Z > 88Q = 4000
Z < 88Q = q
q ≤ 2800Z > 88Q = 2800
Z < 88Q = q
Table 4. Conventional flood control operation of the Xiajiang reservoir.
Table 4. Conventional flood control operation of the Xiajiang reservoir.
ConditionInflow (m3/s)Reservoir Water Level (m)Outflow (m3/s)
Rising limb20,000 < q ≤ 21,50045 < Z < 48.4Q = q
21,500 < q ≤ 23,500Q = 20,000
q ≥ 23,500Q = 22,000
q ≤ 22,00048.4 < Z < 49Q = q
22,000 < q ≤ 24,000Q = 22,000
q > 24,000Q = 24,000
q < 26,600Z > 49Q = q
q ≥ 26,600Open discharge
Falling limbq > 19,000-Open discharge
Q ≤ 19,000-Q = 19,000
Table 5. Results of Wan’an–Xiajiang cascade reservoirs using conventional flood control operation.
Table 5. Results of Wan’an–Xiajiang cascade reservoirs using conventional flood control operation.
Return Periods (Year)Wan’an ReservoirXiajiang Reservoir
Maximal Water Level (m)Flood Peak Reduction (m3/s)Maximal Water Level (m)Flood Peak Reduction (m3/s)
Typical88.6227145.000
2090.095345.000
5091.3538245.41323
10091.8766346.481044
20092.4478347.862004
30092.6693748.291819
50093.43121548.751339
Table 6. Results of Wan’an–Xiajiang cascade reservoirs using DP-POA.
Table 6. Results of Wan’an–Xiajiang cascade reservoirs using DP-POA.
Return Periods (Year)Wan’an ReservoirXiajiang Reservoir
Maximal Water Level (m)Flood Peak Reduction (m3/s)Maximal Water Level (m)Flood Peak Reduction (m3/s)
Typical88.6248145.000
2090.0982745.000
5091.35112045.41100
10091.87109246.48324
20092.44164947.861058
30092.66164948.291195
50093.43177148.751161
Table 7. Flood peak at the Ji’an flood control point.
Table 7. Flood peak at the Ji’an flood control point.
Return Periods (Year)Designed Flood Peak (m3/s)OperationReduction (m3/s)Rate (%)
Typical16,032Convention3652.28
DP-POA14919.30
2017,748Convention5983.37
DP-POA205911.60
5020,235Convention9934.91
DP-POA249112.31
10021,950Convention12745.80
DP-POA242511.05
20023,665Convention13285.61
DP-POA258310.91
30024,608Convention13515.49
DP-POA262910.68
50025,808Convention15005.81
DP-POA25419.85
Table 8. Flood peak at the Shishang flood control point.
Table 8. Flood peak at the Shishang flood control point.
Return Periods (Year)Designed Flood Peak (m3/s)OperationReduction (m3/s)Rate (%)
Typical18,059Convention4822.67
DP-POA11026.10
2019,856Convention6273.16
DP-POA16298.20
5022,637Convention11945.27
DP-POA234510.36
10024,556Convention19617.98
DP-POA308812.57
20026,474Convention333612.60
DP-POA417115.75
30027,529Convention408114.82
DP-POA466516.95
50028,873Convention457515.85
DP-POA535918.56
Table 9. Flood peak at the Waizhou flood control point.
Table 9. Flood peak at the Waizhou flood control point.
Return Periods (Year)Designed Flood Peak (m3/s)OperationReduction (m3/s)Rate (%)
Typical18,839Convention7473.96
DP-POA10805.73
2020,700Convention9124.40
DP-POA15277.38
5023,600Convention13705.81
DP-POA20908.86
10025,600Convention20768.11
DP-POA292011.41
20027,600Convention308511.18
DP-POA383313.89
30028,700Convention344912.02
DP-POA422814.73
50030,100Convention406313.50
DP-POA480915.98

Share and Cite

MDPI and ACS Style

Zhu, D.; Mei, Y.; Xu, X.; Chen, J.; Ben, Y. Optimal Operation of Complex Flood Control System Composed of Cascade Reservoirs, Navigation-Power Junctions, and Flood Storage Areas. Water 2020, 12, 1883. https://doi.org/10.3390/w12071883

AMA Style

Zhu D, Mei Y, Xu X, Chen J, Ben Y. Optimal Operation of Complex Flood Control System Composed of Cascade Reservoirs, Navigation-Power Junctions, and Flood Storage Areas. Water. 2020; 12(7):1883. https://doi.org/10.3390/w12071883

Chicago/Turabian Style

Zhu, Di, Yadong Mei, Xinfa Xu, Junhong Chen, and Yue Ben. 2020. "Optimal Operation of Complex Flood Control System Composed of Cascade Reservoirs, Navigation-Power Junctions, and Flood Storage Areas" Water 12, no. 7: 1883. https://doi.org/10.3390/w12071883

APA Style

Zhu, D., Mei, Y., Xu, X., Chen, J., & Ben, Y. (2020). Optimal Operation of Complex Flood Control System Composed of Cascade Reservoirs, Navigation-Power Junctions, and Flood Storage Areas. Water, 12(7), 1883. https://doi.org/10.3390/w12071883

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