Next Article in Journal
Visible Light Photocatalytic Degradation of Environmental Pollutants Using Zn-Doped NiO Nanoparticles
Previous Article in Journal
Effects of Aquaculture and Thalassia testudinum on Sediment Organic Carbon in Xincun Bay, Hainan Island
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Flood Simulation of the Modified Muskingum Model with a Variable Exponent Based on the Artificial Rabbit Optimization Algorithm

1
School of Water Conservancy, North China University of Water Resources and Electric Power, Zhengzhou 450046, China
2
School of Materials Science and Engineering, North China University of Water Resources and Electric Power, Zhengzhou 450046, China
3
College of Water Resources, North China University of Water Resources and Electric Power, Zhengzhou 450046, China
*
Author to whom correspondence should be addressed.
Water 2024, 16(2), 339; https://doi.org/10.3390/w16020339
Submission received: 14 December 2023 / Revised: 16 January 2024 / Accepted: 18 January 2024 / Published: 19 January 2024

Abstract

:
In order to further improve the accuracy of flood routing, this article uses the Variable Exponential Nonlinear Muskingum Model (VEP-NMM), combined with the Artificial Rabbit Optimization (ARO) algorithm for parameter calibration, to construct the ARO-VEP-NMM flood routing model. Taking Wilson’s (1974) flood as an example, the model calculation results were compared and analyzed with the Muskingum model constructed with seven optimization algorithms. At the same time, six measured floods in the Zishui Basin were selected for model applicability testing. The results show that the ARO algorithm exhibits stronger robustness and search ability compared with other optimization algorithms and can better solve the parameter optimization problem of the Muskingum model. The use of the ARO-VEP-NMM model for flood routing accurately reflects the movement patterns of floods. The Nash coefficient of the Wilson section reached 0.9983, and the average Nash coefficient during the flood validation period in the Zishui Basin was 0.9, further verifying the adaptability and feasibility of the ARO-VEP-NMM model in flood routing. The research results can provide certain references and a theoretical basis for improving the accuracy of flood forecasting.

1. Introduction

Flood disasters are currently the most significant natural calamities, causing substantial losses to humanity worldwide. With the impact of global climate change, frequent occurrences of intense rainfall events, tsunamis, and other disasters have led to an increasing number of extreme flooding incidents, highlighting the growing prominence of flood control and disaster reduction [1,2,3,4]. In order to mitigate the impact of flood disasters, using real-time and accurate flood forecasting is an effective and crucial measure. Among these, river flood simulation aims to predict the propagation and evolution of floods within river channels based on upstream inflow forecasts, thereby providing services for flood prediction and flood control scheduling [5,6]. At present, flood simulation is mainly divided into hydraulic methods and hydrological methods. The former, based on the Saint-Venant equations, requires a substantial amount of data and involves lengthy computation times. The latter, relying on principles of mass and momentum conservation, offers relatively simpler models and quicker computations.
The Muskingum model proposed by McCarthy, G.T., is a hydrological method widely applied in flood simulation problems [7,8,9]. It is based on the assumption of a one-dimensional model, mainly considering the flow in the downstream direction to predict the propagation and evolution process of floods in a river channel. Although the Muskingum model performs well in many applications, there are also some assumptions and limitations that need to be considered. Firstly, the Muskingum model is a one-dimensional model, which means that it only focuses on the longitudinal flow of the river during simulation. This results in the model being unable to consider lateral flow in areas such as floodplains and unable to capture the complexity of spatial changes. In some complex river systems, especially those with complex geometric shapes and flow characteristics, the Muskingum model may not perform as expected. However, in reality, rivers are numerous, and with human influences, the shapes of river cross-sections and hydraulic characteristics undergo constant changes, posing challenges to the application of the Muskingum model in flood simulation.
Over the past few decades, researchers have improved the model by introducing additional parameters [10,11,12,13]. These efforts aim to overcome some limitations of the Muskingum model, especially in its application under complex hydraulic conditions. The original model may perform well in some river systems, but as the complexity of the river increases, a model with more degrees of freedom and fitting ability is needed.
Among them, the Variable Exponent Nonlinear Muskingum Model (VEP-NMM) [12], with increased degrees of freedom and better fitting capabilities compared with the original model, has shown enhanced performance under more complex hydraulic conditions. However, in practical applications, the accuracy of the Muskingum model prediction relies on the definition of the model structure and the determination of model parameters. As the number of parameters increases, the estimation of these parameters becomes more complex. Therefore, determining the optimal parameter estimates for the model has become a crucial and challenging aspect of flood simulation.
With the deepening of research, various optimization algorithms are used to precisely calculate the parameters of the Muskingum model [14,15,16,17,18]. However, conventional optimization algorithms often yield similar non-optimal estimates, all falling into local optima. Therefore, the need for more robust algorithms arises to overcome the local optimum dilemma and obtain optimal parameter estimates [19,20,21].
The Artificial Rabbit Optimization algorithm (ARO), proposed by Liying Wang, is a metaheuristic algorithm inspired by the foraging behavior of rabbits [22]. In comparison with other classical algorithms, ARO has demonstrated outstanding performance in optimizing model parameters and handling nonlinear problems, exhibiting stronger global search capabilities and adaptability to complex engineering issues [23,24,25].
Therefore, this study uses the VEP-NMM model, utilizing the ARO algorithm for parameter calibration. The model is solved using various case studies, and comparative analyses with other methods are conducted to validate the model’s advanced and effective performance. The aim is to better simulate the evolution of riverine floods.

2. Materials and Methods

2.1. Variable Exponent Nonlinear Muskingum Model (VEP-NMM)

The Linear Muskingum Model (LMM) proposed by McCarthy, G.T. imposes relatively low requirements on channel roughness, slope, and morphology while still accurately simulating changes in river channel form [8]. The equation for the LMM is as follows:
S t = K [ X I t + ( 1 X ) Q t ]
d S d t Δ S Δ t = I Q
where I t and Q t represent the upstream and downstream flows at time t, respectively; S t is the channel storage; K is the slope of the relationship line between the storage capacity and the indicated storage flow, which has a time dimension and actually represents the propagation time of floods in the river section; and X is the proportion factor of flow, representing the relative weight of upper and lower section flow in the tank storage capacity, reflecting the role of wedge storage in flow calculation.
Research [26,27,28] has demonstrated that when there is a strong nonlinear relationship between channel storage and discharge, the accuracy of the Linear Muskingum Model in flood simulation is not satisfactory, necessitating the treatment of this nonlinear relationship. In this context, Gill [10] proposed the widely applied Nonlinear Muskingum Method (NMM), which introduces a nonlinear parameter m. The equation for NMM is as follows:
S t = K [ X I t + ( 1 X ) Q t ] m
As research advanced, Easa et al. [12] suggested that the parameter m should be a continuous variable exponent parameter whose variation is related to the current inflow and the inflow process. Therefore, by introducing the dimensionless inflow variable u t , the Variable Exponent Nonlinear Muskingum Model (VEP-NMM) was proposed, further enhancing the model’s performance, making the simulation results more consistent with natural river channel characteristics, and exhibiting time-varying properties. The specific computational equation is as follows:
m t = β ( u t ) = a + b e e c u t u t = I t / I m a x I m a x = m a x 1 t T I t
where a, b, and c are the model parameters to be optimized and I m a x is the maximum inflow to the river channel. Then, Equation (3) transforms into:
S t = K [ X I t + ( 1 X ) Q t ] β u t
After initializing the parameters, the value of m t for each time period is obtained. Assuming the initial outflow value is equal to the inflow value, Q0 = I0, the initial channel storage is then calculated using the following equation:
S 0 = K X I 0 + ( 1 X ) I 0 β u 0
The equation for the cumulative channel storage in the next time period is expressed as follows:
S t + 1 = S t + Δ t q i
q i = ( 1 1 X ) S i K 1 β u i 1 + 1 1 X I i
The equation for the calculated outflow is then expressed as follows:
Q t + 1 = 1 1 X S t + 1 K 1 / β ( u t + 1 ) X 1 X I t + 1

2.2. Artificial Rabbit Optimization Algorithm (ARO)

Originating from the survival strategies of rabbits in nature, the Artificial Rabbit Optimization Algorithm (ARO) is inspired by two distinct behaviors: the detour foraging strategy and the random hiding strategy, each associated with a transition between these strategies driven by their energy levels [22]. The detour foraging strategy is used by rabbits to avoid predators, involving foraging away from the burrow with the advantage of their broad field of vision, which facilitates the detection of food over a large area. This foraging strategy is considered the exploration phase of the algorithm. On the other hand, the random hiding strategy involves rabbits digging multiple burrows and randomly selecting one as a refuge. Rabbits may also execute sudden turns or abrupt stops while running, confusing predators and evading pursuit. This hiding strategy is seen as the exploitation phase of the algorithm.
Given the necessity for rabbits to run rapidly to escape danger, resulting in energy changes, they adaptively switch between foraging and evading danger by leveraging their energy variations. This adaptive switching enhances their chances of survival by balancing food acquisition and avoiding threats.
The ARO algorithm incorporates these survival traits of rabbits by simulating the transitions between detour foraging and random hiding to solve problems. This simulation corresponds to three behavioral phases: detour foraging, random hiding, and energy transformation. The mathematical model of the algorithm is represented as follows.

2.2.1. Detour Foraging Strategy (Exploration Phase)

In ARO, each rabbit in the population has a burrow area, and rabbits randomly visit the locations of other individuals for foraging. In reality, while foraging, rabbits often disrupt food sources to obtain an adequate amount of food. Therefore, the detour foraging behavior in ARO implies that each searching individual will randomly select another individual in the population to search and update its position, introducing a disturbance. Thus, the mathematical model for detour foraging is expressed as:
v i t + 1 = x j t + R x i t x j t + r o u n d 0.5 0.05 + r 1 n 1
R = L × c L = e e t 1 T 2 sin ( 2 π r 2 ) c k = 1 ,   i f   k = = g l g = randperm ( D )   0 ,   e l s e
where v i t + 1 is the position of the i-th rabbit at time t + 1; i and j are both numbers from 1 to N, and N is the population size of rabbits; r1, r2, and r3 are random numbers from 0 to 1; n1 is a random number that follows a standard normal distribution; D is the dimension; T is the maximum number of iterations; t is the current number of iterations; and k = 1...D and l = 1... [r3 *D].

2.2.2. Random Hiding Strategy (Exploitation Phase)

To evade predators, rabbits typically dig various burrows around their burrows. In each iteration of ARO, each rabbit generates holes along every dimension of the search space and randomly chooses one of these holes to hide in, reducing the likelihood of being preyed upon. The mathematical model for generating burrows is as follows:
b i , j t = x i t + H g x i t ,   i = 1 , , n   and   j = 1 , , D
H = T t + 1 T r 4
g k = 1 , i f   k = j 0 , e l s e k = 1 , , D
where b i , j t is the j-th cave generated by the i-th rabbit and H is a parameter that gradually decreases linearly with the number of algorithm iterations. Influenced by this parameter, the rabbit’s initial burrow occurs in a larger domain, which decreases as the iteration progresses.
As mentioned above, predators often chase and attack rabbits, so rabbits must find a safe hiding place to survive. The modeling of this random hiding strategy is as follows:
v i t + 1 = x i t + R r 4 b i , r t x i t , i = 1 , , N
b i , r t = x i t + H g r x i t
g r k = 1 , i f   k = = r 5 D 0 , e l s e k = 1 , , D
where r4 and r5 are both random numbers from 0 to 1. The random selection of hidden caves is determined using Equation (17), and the meanings of other parameters are the same as before.
After selecting either the detour foraging strategy or the random hiding strategy, the update of the position for the i-th rabbit is as follows:
x i t + 1 = x i t ,                 i f   f x i t f v i t + 1 v i t + 1 ,   i f   f x i t > f v i t + 1

2.2.3. Energy Factor A

In ARO, rabbits tend to use the detour foraging strategy in the early iterations and switch to the random hiding strategy in the later iterations. This mechanism is related to the energy level within the rabbits. Therefore, factor A is designed to implement the transition from the exploration phase to the exploitation phase. The energy coefficient is defined as follows:
A t = 4 1 t T ln 1 r
where the energy conversion factor A indicates the rabbit’s endurance and energy for detour foraging when higher and the need for random hiding when lower. In the equation, r is a random number between 0 and 1, and A decreases to zero with the growth in the oscillation amplitude during the iteration process.
In ARO, when A > 1, rabbits tend to randomly explore different areas for detour foraging; when A < 1, rabbits prefer the random hiding strategy. Thus, based on the energy conversion factor A, ARO transitions between detour foraging and random hiding strategies. To demonstrate the balance between exploration and exploitation in the algorithm, the probability of executing the detour foraging strategy is calculated using the following equation:
P { A t > 1 } = 0 2 0 e 1 2 k d r d θ 1 2 = 1 4 1 4 e t t d t + e 1 4 0.5177
Therefore, during the iteration process, the probability of ARO executing the detour foraging strategy is approximately 0.5. This means that the probabilities of executing detour foraging and random hiding are nearly equal, contributing to balanced exploration and exploitation in the algorithm.
Hence, compared with other optimization algorithms, the ARO algorithm is characterized by its simplicity in principles, fewer parameters, and effective balancing of exploration and exploitation capabilities with the energy transformation factor A. It uses corresponding detour foraging and random hiding strategies, achieving better global exploration and local exploitation to seek optimal global solutions.

2.3. ARO-VEP-NMM Flood Simulation Model

Based on the analysis of the VEP-NMM model in Section 2.1, given the inflow I from upstream and outflow Q from downstream of the river, the ARO algorithm is used to calibrate the model parameters K, X, a, b, and c. After calibration, these parameters are substituted into the model to obtain the simulated flow Qs from upstream to downstream. The simulated flow is then compared with the observed flow, and if the criteria are not met, parameter calibration continues until satisfactory results are achieved. The workflow of the ARO-VEP-NMM flood simulation model is illustrated in Figure 1.
To achieve a more accurate simulated outflow, the optimization process aims to minimize the sum of squared differences (SSQ) between the observed and simulated outflow. The mathematical model for this optimization objective function is as follows:
S S Q = t + 1 T Q ( t ) Q s ( t ) 2
where Q t is the measured flow rate at time t; and the flow rate is calculated at time t for Q s t .

2.4. Model Evaluation Metrics

To assess the reliability of flood simulation results, multiple metrics are used. Tailoring the evaluation to the characteristics of different flood instances enhances the intuitiveness of the results. For case 1, the sum of squared differences (SSQ), the sum of absolute differences (SAD) between the observed and simulated flow, and the Nash coefficient (NSE) are used as evaluation metrics. For case 2, which involves multiple flood events, metrics such as peak flow, the time lag of peak flow, the peak flow relative error (QE), and flood forecasting accuracy based on peak flow and peak time are used. The Nash–Sutcliffe Efficiency (NSE) coefficient is used as an indicator to judge the fitting degree between forecasted and observed data. The mathematical models for each evaluation metric are as follows:
S A D = t + 1 T Q ( t ) Q s ( t )
Q E = Q c ( i ) Q 0 ( i ) Q 0 ( i ) × 100 %
N S E = 1 i = 1 N Q c i Q 0 i 2 i = 1 N Q 0 i Q 0 ¯ 2
where Q 0 i and Q c i are the measured and predicted flow data, respectively, and the definitions of other parameters are the same as before.

3. Results and Discussion

The establishment of the ARO-VEP-NMM model ultimately aims to effectively improve accuracy in river flood simulation. Building upon the analyses of the VEP-NMM model and the ARO algorithm in Section 2, an ARO-VEP-NMM flood simulation model is constructed. Two instances are considered for validation: Case 1 uses flood data from the Wilson River (1974) and Case 2 uses recurrent floods in the Zishui River basin. These instances are used to assess the flood simulation capabilities and practical effectiveness of the ARO-VEP-NMM model.

3.1. Case 1—Wilson River (1974)

The Wilson River (1974) data represent a smooth single-peak hydrograph with a pronounced non-linear relationship, making it widely used to validate the effectiveness of various parameter optimization estimation techniques in the Muskingum model [29,30].
To verify the feasibility and effectiveness of the ARO algorithm in parameter calibration for the VEP-NMM model, flood data from the Wilson River segment is chosen as an example for testing. A comprehensive comparison is conducted with a series of advanced optimization algorithms, including the Arithmetic Optimization Algorithm [31] (AOA), the Harris Hawks Optimization Algorithm [32] (HHO), and Aquila Optimize (AO) [33]. In order to thoroughly examine the performance of the ARO-VEP-NMM model, the simulation results of the ARO-VEP-NMM model are compared with a range of optimization algorithms used by past scholars, including the Genetic Algorithm (GA) [34], the Harmony Search Algorithm (HS) [35], Particle Swarm Optimization (PSO) [36], and Differential Evolution (DE) [37]. This aims to validate the adaptability and superiority of the ARO-VEP-NMM model in flood simulation.
The parameter calibration for the ARO-VEP-NMM model and other algorithms is conducted based on the mentioned theory, and the parameter settings for each algorithm are consistent with the references. Given the relatively fewer control parameters for algorithms like ARO, HHO, and AO, and their strong stochastic nature, each model is independently run 10 times to avoid falling into the local optimum dilemma and to seek the optimal simulation results.
Based on the introduction of the ARO-VEP-NMM model in the previous text, the ARO algorithm is used to calibrate the model parameters K, X, a, b, and c, given the inflow I upstream and outflow Q downstream of the Wilson River. After calibration, the parameters are input into the model to obtain the calculated flow rate Qs from upstream to downstream. The calculated flow rate is then evaluated and compared with the measured flow rate. If it does not meet the requirements, parameter calibration is continued until the results meet the evaluation criteria. Using parameter calibration, the optimal parameters for the model are determined as follows: K = 0.7091, X = 0.2750, a = 1.8057, b = 12.1784, and c = 8.9845.
To evaluate the performance of the model, the best value, average value, standard deviation, and solution time of the statistical model after 10 independent runs of the model are compiled as evaluation indicators. The comparative results are presented in Table 1. For a more intuitive understanding of the simulation effects of each model, Figure 2 provides box plots and normal distribution plots comparing the SSQ results after 10 runs for each model.
From Table 1, it can be observed that the optimal SSQ values for ARO-VEP-NMM, AOA-VEP-NMM, HHO-VEP-NMM, and AO-VEP-NMM are 20.47, 20.55, 27.69, and 36.63, respectively. Among these, ARO-VEP-NMM yields the optimal computation results, with the average value closely approaching the optimal value. Compared with the other algorithms, it demonstrates superior optimization accuracy and convergence capability. Additionally, the standard deviation of ARO-VEP-NMM is 1.63, indicating a relatively small degree of dispersion and higher stability compared with the other three algorithms. However, it is worth noting that compared with the other algorithms, the ARO algorithm requires a longer solving time.
As depicted in Figure 2, the SSQ values for HHO-VEP-NMM and AO-VEP-NMM fluctuate significantly across multiple runs, exhibiting substantial variations between the highest and lowest values, reflecting lower model stability. AOA-VEP-NMM exhibits relatively greater stability. In contrast, the ARO-VEP-NMM model maintains consistent SSQ values across multiple runs, indicating superior stability. Furthermore, its SSQ values are lower compared to the other three algorithms, highlighting its superior applicability.
Given that the calculations for the Genetic Algorithm (GA), Harmony Search Algorithm (HS), Particle Swarm Algorithm (PSO), and Differential Evolution algorithm (DE) were already conducted in reference [26], this paper directly references the numerical results from that study for a convenient comparison of the performance of the ARO-VEP-NMM model with other methods. Table 2 and Table 3 present the outflow process and calibration results for various methods in the Wilson case.
Based on the data in Table 2 and Table 3, it is evident that the ARO-VEP-NMM model outperforms other methods significantly in both SSQ and SAD aspects. Furthermore, the simulation results of ARO-VEP-NMM exhibit the highest accuracy, with the smallest simulation error and the best fit to the observed inflow curve. This outcome further validates the outstanding performance and efficient computational capabilities of the ARO algorithm in parameter calibration.
Additionally, it is noteworthy that the VEP-NMM model’s performance in flood simulation surpasses that of the traditional NMM model. This indicates that the chosen variable exponent parameter model is feasible and effective in enhancing flood simulation accuracy. Therefore, the ARO-VEP-NMM model demonstrates significant advantages in flood simulation and proves to be a viable approach for flood forecasting.

3.2. Case 2—Site Floods in the Zishui River Basin

The Zishui River Basin is located within the territory of Shaoyang City, Hunan Province, China. Due to the complex topography within the river channel and the interference of tributaries, the basin exhibits intricate channel morphology and irregular flood propagation characteristics. To further validate the practical simulation performance of the ARO-VEP-NMM model, the Luo Jiamiao hydrological station is chosen as the outlet section of the basin. Six sets of observed flood data from Tangdukou and Longhui hydrological stations from 2014 to 2015 were selected for flood simulation.
After analyzing the data from the six flood events, they are divided into two different periods. The first four flood events are used for the calibration period, while the last two are reserved for the validation period. The ARO algorithm is used to calibrate the parameters, resulting in the optimal parameter values: K = 1.1039, X = 0.5000, a = 1.2618, b = 0.3487, and c = 0.8914. The simulation results for each flood event are presented in Table 4.
From Table 4, it is evident that for the ARO-VEP-NMM flood simulations at the Luo Jiamiao station, the relative errors in the flood peaks during the calibration period are within an acceptable range, with an average of 7.53%. Additionally, the Nash coefficient averages 0.73 for the calibration period. Despite one flood event (“27 June 2014”) during the calibration period having a peak time difference exceeding two time intervals, failing that particular flood simulation, the overall pass rate for the floods during the calibration period reaches 75%. Furthermore, both flood events during the validation period passed the criteria, with an average relative error of 8.93% and a Nash coefficient average of 0.90. Overall, the flood simulation results are reasonably accurate.
Although only six flood events were selected in this study, these events can reflect the overall characteristics of the basin’s floods, and the simulation results can serve as important references for downstream flood control. The simulation hydrographs for the validation period are shown in Figure 3 and Figure 4.
From Figure 3 and Figure 4, it is evident that both floods during the calibration period are of the single-peak type, and their simulation results are excellent, closely fitting the actual flow curves. The overall trend in flood simulation closely resembles the original Luo Jiamiao flow curve, indicating good simulation performance of the ARO-VEP-NMM model for these two floods. Additionally, the small differences between the simulated and actual peak values further demonstrate the accuracy and adaptability of the model’s flood simulation.

3.3. Parameter Sensitivity Analysis

In order to comprehensively evaluate the performance of the ARO algorithm in the VEP-NMM model, this section describes a parameter sensitivity analysis of the ARO algorithm for each case. Due to the fact that the ARO algorithm itself does not involve additional parameters, the population size N of the algorithm was considered the main parameter for sensitivity analysis when solving the case.
Multiple sets of experiments were conducted on N, and each set of experiments was run independently 10 times. Table 5 shows the SSQ of Case 1 (Wilson River section), recording the average, best, worst, standard deviation, and solution time of the output results. Meanwhile, Table 6 shows the QE, NSE, and solution time of flood events in the Zishui Basin during the regular and validation periods.
From the data in the tables, it can be observed that as N gradually increases, all evaluation indicators improve, indicating that the increase in population size increases the range of the ARO algorithm in finding the optimal parameters. Especially at N = 100, all results reach the optimal level. Further increasing the number of groups at this point does not improve the indicators, but rather leads to a continued increase in the solving time and the consumption of additional computing resources. Therefore, in this article, it is found that selecting the ARO algorithm with a population size of N of 100 can achieve beneficial results. It is worth noting that since the ARO algorithm does not have additional parameters, it does not need to consider the rationality of the proportional parameter values during its own evolution process. This is also one of the important reasons why this article chooses ARO as the tool for finding the optimal parameters.

4. Conclusions

This study aimed to improve the accuracy of river flood simulations by introducing the Variable Exponent Nonlinear Muskingum Model (VEP-NMM). To address the parameter optimization challenge of the VEP-NMM model, this study used the Artificial Rabbit Optimization algorithm (ARO) to construct a variable exponent Nonlinear Muskingum Model based on ARO (ARO-VEP-NMM). The feasibility and adaptability of the ARO-VEP-NMM model were tested using two specific case studies, and the conclusions are as follows:
  • Using the VEP-NMM model can enhance flood simulation accuracy, with parameter optimization being a crucial aspect of the model’s flood simulation. By utilizing the ARO algorithm for parameter optimization and testing with the Wilson (1974) River segment data for flood simulation, a comparison of seven different optimization methods revealed that the ARO algorithm provides higher optimization accuracy and better robustness. This offers a more effective method for Muskingum model parameter optimization.
  • The ARO-VEP-NMM model, applied to simulate measured flood events in the Zishui River basin, accurately reproduces the movement and characteristics of floods. This further validates the excellent applicability of the ARO-VEP-NMM model in flood simulation.
Therefore, the proposed ARO-VEP-NMM flood simulation method combines the advantages of the VEP-NMM model and the ARO algorithm. The research results indicate that their integration can effectively conduct flood simulation, offering a new approach to flood prediction.
However, the model studied in this article may not have fully taken into account the backwater effect caused by hydraulic structures such as bridges and culverts. The impact of these structures on flood routing may be an important factor, and future research will further explore this aspect. In addition, it was noted in this study that the ARO algorithm may have higher solving time and may consume more computing resources compared with other algorithms. Therefore, efforts will be made to improve the ARO algorithm or consider combining optimization with other techniques (such as machine learning) to optimize its convergence speed while improving the accuracy of the results. This is also a new direction for future research.

Author Contributions

M.L.: methodology, program implementation, and writing—original draft. Z.C.: conceptualization, data curation, and writing—original draft preparation. T.F.: writing and editing—original draft. All authors have read and agreed to the published version of this manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to privacy.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Yin, J.B.; Guo, S.L.; Gentine, P.; Sullivan, S.C.; Gu, L.; He, S.K.; Chen, J.; Liu, P. Does the Hook Structure Constrain Future Flood Intensification Under Anthropogenic Climate Warming? Water Resour. Res. 2021, 57, e2020WR028491. [Google Scholar] [CrossRef]
  2. Wang, W.; Tian, W.; Chau, K.; Zang, H.; Ma, M.; Feng, Z.; Xu, D. Multi-Reservoir Flood Control Operation Using Improved Bald Eagle Search Algorithm with ε Constraint Method. Water 2023, 15, 692. [Google Scholar] [CrossRef]
  3. Kundzewicz, Z.W.; Su, B.; Wang, Y.; Xia, J.; Huang, J.; Jiang, T. Flood risk and its reduction in China. Adv. Water Resour. 2019, 130, 37–45. [Google Scholar] [CrossRef]
  4. Wang, Y.; Imai, K.; Miyashita, T.; Ariyoshi, K.; Takahashi, N.; Satake, K. Coastal tsunami prediction in Tohoku region, Japan, based on S-net observations using artificial neural network. Earth Planets Space 2023, 75, 154. [Google Scholar] [CrossRef]
  5. Sun, K.; Hu, L.; Guo, J.; Yang, Z.; Zhai, Y.; Zhang, S. Enhancing the understanding of hydrological responses induced by ecological water replenishment using improved machine learning models: A case study in Yongding River. Sci. Total Environ. 2021, 768, 145489. [Google Scholar] [CrossRef] [PubMed]
  6. Wang, W.; Tian, W.; Xu, L.; Liu, C.; Xu, D. Mε-OIDE algorithm for solving constrained optimization problems and its application in flood control operation of reservoir group. Shuili Xuebao/J. Hydraul. Eng. 2023, 54, 148–158. [Google Scholar] [CrossRef]
  7. Norouzi, H.; Bazargan, J. Calculation of Water Depth during Flood in Rivers using Linear Muskingum Method and Particle Swarm Optimization (PSO) Algorithm. Water Resour. Manag. 2022, 36, 4343–4361. [Google Scholar] [CrossRef]
  8. McCarthy, G.T. The unit hydrograph and flood routing. In Proceedings of the Conference of North Atlantic Division, Wahsington, DC, USA, 24 June 1938. [Google Scholar]
  9. Lee, E.H. Development of a New 8-Parameter Muskingum Flood Routing Model with Modified Inflows. Water 2021, 13, 3170. [Google Scholar] [CrossRef]
  10. Gill, M.A. Flood routing by the Muskingum method. J. Hydrol. 1978, 36, 353–363. [Google Scholar] [CrossRef]
  11. Perumal, M.; Price, R.K. A fully mass conservative variable parameter McCarthy–Muskingum method: Theory and verification. J. Hydrol. 2013, 502, 89–102. [Google Scholar] [CrossRef]
  12. Easa, S.M. Evaluation of nonlinear Muskingum model with continuous and discontinuous exponent parameters. Ksce J. Civ. Eng. 2015, 19, 2281–2290. [Google Scholar] [CrossRef]
  13. Easa, S.M. Improved Nonlinear Muskingum Model with Variable Exponent Parameter. J. Hydrol. Eng. 2013, 18, 1790–1794. [Google Scholar] [CrossRef]
  14. Moradi, E.; Yaghoubi, B.; Shabanlou, S. A new technique for flood routing by nonlinear Muskingum model and artificial gorilla troops algorithm. Appl. Water Sci. 2023, 13, 49. [Google Scholar] [CrossRef]
  15. Okkan, U.; Kirdemir, U. Locally tuned hybridized particle swarm optimization for the calibration of the nonlinear Muskingum flood routing model. J. Water Clim. Chang. 2020, 11, 343–358. [Google Scholar] [CrossRef]
  16. Yuan, G.; Lu, J.; Wang, Z. The modified PRP conjugate gradient algorithm under a non-descent line search and its application in the Muskingum model and image restoration problems. Soft Comput. 2021, 25, 5867–5879. [Google Scholar] [CrossRef]
  17. Zhang, S.; Kang, L.; Zhou, L.; Guo, X. A new modified nonlinear Muskingum model and its parameter estimation using the adaptive genetic algorithm. Hydrol. Res. 2017, 48, 17–27. [Google Scholar] [CrossRef]
  18. Wang, W.; Xu, Z.; Qiu, L.; Xu, D. Hybrid Chaotic Genetic Algorithms for Optimal Parameter Estimation of Muskingum Flood Routing Model. In Proceedings of the 2009 International Joint Conference on Computational Sciences and Optimization, Sanya, China, 24–26 April 2009; pp. 215–218. [Google Scholar]
  19. Wang, W.; Tian, W.; Chau, K.-W.; Xue, Y.; Xu, L.; Zang, H. An Improved Bald Eagle Search Algorithm with Cauchy Mutation and Adaptive Weight Factor for Engineering Optimization. CMES—Comput. Model. Eng. Sci. 2022, 136, 1603–1642. [Google Scholar] [CrossRef]
  20. Ehteram, M.; Othman, F.B.; Yaseen, Z.M.; Afan, H.A.; Allawi, M.F.; Malek, M.B.A.; Ahmed, A.N.; Shahid, S.; Singh, V.P.; El-Shafie, A. Improving the Muskingum Flood Routing Method Using a Hybrid of Particle Swarm Optimization and Bat Algorithm. Water 2018, 10, 807. [Google Scholar] [CrossRef]
  21. Wang, W.-C.; Xu, L.; Chau, K.-w.; Zhao, Y.; Xu, D.-M. An orthogonal opposition-based-learning Yin–Yang-pair optimization algorithm for engineering optimization. Eng. Comput. 2022, 38, 1149–1183. [Google Scholar] [CrossRef]
  22. Wang, L.Y.; Cao, Q.J.; Zhang, Z.X.; Mirjalili, S.; Zhao, W.G. Artificial rabbits optimization: A new bio-inspired meta-heuristic algorithm for solving engineering optimization problems. Eng. Appl. Artif. Intell. 2022, 114, 105082. [Google Scholar] [CrossRef]
  23. Abd Elaziz, M.; Dahou, A.; Mabrouk, A.; El-Sappagh, S.; Aseeri, A.O. An Efficient Artificial Rabbits Optimization Based on Mutation Strategy For Skin Cancer Prediction. Comput. Biol. Med. 2023, 163, 107154. [Google Scholar] [CrossRef]
  24. Dangi, D.; Chandel, S.T.; Dixit, D.K.; Sharma, S.; Bhagat, A. An efficient model for sentiment analysis using artificial rabbits optimized vector functional link network. Expert Syst. Appl. 2023, 225, 119849. [Google Scholar] [CrossRef]
  25. Ozkaya, B.; Duman, S.; Kahraman, H.T.; Guvenc, U. Optimal solution of the combined heat and power economic dispatch problem by adaptive fitness-distance balance based artificial rabbits optimization algorithm. Expert Syst. Appl. 2024, 238, 122272. [Google Scholar] [CrossRef]
  26. Niazkar, M.; Afzali, S.H. Parameter estimation of an improved nonlinear Muskingum model using a new hybrid method. Hydrol. Res. 2017, 48, 1253–1267. [Google Scholar] [CrossRef]
  27. Swain, R.; Sahoo, B. Variable parameter McCarthy–Muskingum flow transport model for compound channels accounting for distributed non-uniform lateral flow. J. Hydrol. 2015, 530, 698–715. [Google Scholar] [CrossRef]
  28. Easa, S.M.; Barati, R.; Shahheydari, H.; Nodoshan, E.J.; Barati, T. Discussion: New and improved four-parameter non-linear Muskingum model. Proc. Inst. Civ. Eng. Water Manag. 2014, 167, 612–615. [Google Scholar] [CrossRef]
  29. Wang, W.C.; Tian, W.C.; Xu, D.M.; Chau, K.W.; Ma, Q.; Liu, C.J. Muskingum Models’ Development and their Parameter Estimation: A State-of-the-art Review. Water Resour. Manag. 2023, 37, 3129–3150. [Google Scholar] [CrossRef]
  30. Lu, C.; Ji, K.; Wang, W.; Zhang, Y.; Ealotswe, T.K.; Qin, W.; Lu, J.; Liu, B.; Shu, L. Estimation of the Interaction Between Groundwater and Surface Water Based on Flow Routing Using an Improved Nonlinear Muskingum-Cunge Method. Water Resour. Manag. 2021, 35, 2649–2666. [Google Scholar] [CrossRef]
  31. Abualigah, L.; Diabat, A.; Mirjalili, S.; Abd Elaziz, M.; Gandomi, A.H. The Arithmetic Optimization Algorithm. Comput. Methods Appl. Mech. Eng. 2021, 376, 113609. [Google Scholar] [CrossRef]
  32. Heidari, A.A.; Mirjalili, S.; Faris, H.; Aljarah, I.; Mafarja, M.; Chen, H.L. Harris hawks optimization: Algorithm and applications. Future Gener. Comput. Syst. Int. J. Escience 2019, 97, 849–872. [Google Scholar] [CrossRef]
  33. Abualigah, L.; Yousri, D.; Abd Elaziz, M.; Ewees, A.A.; Al-qaness, M.A.A.; Gandomi, A.H. Aquila Optimizer: A novel meta-heuristic optimization algorithm. Comput. Ind. Eng. 2021, 157, 107250. [Google Scholar] [CrossRef]
  34. Guojing, Z. Exploration of parameter estimation methods for the Maskingen model. J. China Hydrol. 1997, 3, 41–43+37. [Google Scholar] [CrossRef]
  35. Rui, W.; Jun, X.; Wenhua, Z. Application of Harmory Search Method in Parameter Estimation Based on Nonlinear Muskingum Model. Water Resour. Power 2008, 4, 36–39. [Google Scholar]
  36. Xixia, M.; Dandan, S.; Yugui, H. Parameter Estimation Method of Nonlinear Muskingum Model Based on PSO. J. Zhengzhou Univ. (Eng. Sci.) 2007, 4, 122–125. [Google Scholar]
  37. Xu, D.M.; Qiu, L.; Chen, S.Y. Estimation of Nonlinear Muskingum Model Parameter Using Differential Evolution. J. Hydrol. Eng. 2012, 17, 348–353. [Google Scholar] [CrossRef]
Figure 1. Flood routing flowchart of the ARO-VEP-NMM model.
Figure 1. Flood routing flowchart of the ARO-VEP-NMM model.
Water 16 00339 g001
Figure 2. Box plot comparison of SSQ for four models after 10 operations.
Figure 2. Box plot comparison of SSQ for four models after 10 operations.
Water 16 00339 g002
Figure 3. Flood routing results for the “10 June 2015” site.
Figure 3. Flood routing results for the “10 June 2015” site.
Water 16 00339 g003
Figure 4. Flood routing results for the “2 July 2015” site.
Figure 4. Flood routing results for the “2 July 2015” site.
Water 16 00339 g004
Table 1. Comparison of SSQ results after 10 independent runs of each model.
Table 1. Comparison of SSQ results after 10 independent runs of each model.
IndexAROAOAHHOAO
20.4721.4647.3477.13
25.2135.6327.6965.28
20.4722.8143.9283.04
20.4738.5068.00122.33
20.4721.06104.86106.29
20.5342.2952.6380.51
23.1420.5572.2147.51
20.4724.3837.4889.49
20.4720.74153.53124.75
20.4721.0229.9836.63
Best20.4720.5527.6936.63
Average21.2126.8463.7683.30
Standard deviation1.638.4839.1029.13
Time (s)6.352.542.973.98
Table 2. Flow results calculated with different methods in the Wilson River (1974) case.
Table 2. Flow results calculated with different methods in the Wilson River (1974) case.
Time Period
/(h)
Actual Inflow
/(m3/s)
Actual Measured Outflow
/(m3/s)
Flood Routing Calculates Flow Rate
/(m3/s)
VEP-NMMNMM
AROAOAHHOAODEGAHSPSO
022222222222222222222
6232122.722.722.922.8222222.722
12352123.723.724.22422.422.423.522.6
18712626.526.526.125.926.626.326.128.1
241033433.633.632.432.434.534.235.932.2
301114443.443.442.542.944.244.245.245
361095556.156.155.756.556.956.95757
421006667.267.267.268.368.168.267.567.5
48867576.276.276.277.477.177.176.275.9
54718282.382.382.283.483.383.281.181.2
60598584.984.884.585.485.985.74.585.6
66478483.683.58383.584.584.283.584.2
72398079.879.87979.180.680.280.179.6
78327373.373.272.572.273.773.373.373.3
84286465.465.364.66465.46565.465
90245455.155.254.253.45655.855.756.2
96224444.844.943.943.546.746.746.146.5
102213636.636.636.23637.83836.637.3
108203029.729.729.829.730.530.929.529.7
114192524.524.424.624.925.225.724.324.3
120192222.922.823.423.421.722.120.820.6
126181919.319.119.519.92020.22019.6
NSE0.99830.99830.99780.99700.99690.99700.46690.9970
Table 3. Comparison of parameter calibration results for various methods in the Wilson River (1974) case.
Table 3. Comparison of parameter calibration results for various methods in the Wilson River (1974) case.
ModelARO-VEP-NMMAOA-VEP-NMMHHO-VEP-NMMAO-VEP-NMMDE-NMMGA-NMMHS-NMMPSO-NMM
SSQ20.4720.5527.6936.6336.7738.2336.7836.89
SAD16.4416.6418.2321.5523.462323.424.1
Table 4. Flood simulation results for Tangdu, Longhui, and Luo Jiamiao River sections.
Table 4. Flood simulation results for Tangdu, Longhui, and Luo Jiamiao River sections.
Flood SiteMeasured Peak Value
/(m3/s)
Simulated Peak/(m3/s)QE
/(%)
Peak Time Difference
/(h)
NSEQualified or Not
Regular rate20 June 201451804715.058.9800.96Yes
27 June 20141378.71290.266.42−40.39No
5 July 201418601761.215.31−20.74Yes
8 June 201531703468.679.42−30.82Yes
Average 7.53 0.73
Validation period10 June 201517011726.871.5200.93Yes
2 July 201533383883.2416.33−20.87Yes
Average 8.93 0.9
Table 5. Case 1 parameter sensitivity analysis.
Table 5. Case 1 parameter sensitivity analysis.
NBestAverageWorstStandard DeviationTime(s)
3020.465721.213925.20261.63323.17
5020.465720.551321.32150.27064.28
10020.465720.465720.46576.86 × 10−146.35
20020.465720.465720.46572.23 × 10−1413.35
30020.465720.465720.46575.32 × 10−1415.67
Table 6. Case 2 parameter sensitivity analysis.
Table 6. Case 2 parameter sensitivity analysis.
Regular RateValidation Period
NQENSEQENSETime(s)
307.730.719.120.898.27
507.670.719.010.899.98
1007.530.738.930.9011.13
2007.530.738.930.9014.54
3007.530.738.930.9016.63
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

Li, M.; Cui, Z.; Fan, T. The Flood Simulation of the Modified Muskingum Model with a Variable Exponent Based on the Artificial Rabbit Optimization Algorithm. Water 2024, 16, 339. https://doi.org/10.3390/w16020339

AMA Style

Li M, Cui Z, Fan T. The Flood Simulation of the Modified Muskingum Model with a Variable Exponent Based on the Artificial Rabbit Optimization Algorithm. Water. 2024; 16(2):339. https://doi.org/10.3390/w16020339

Chicago/Turabian Style

Li, Min, Zhirui Cui, and Tianyu Fan. 2024. "The Flood Simulation of the Modified Muskingum Model with a Variable Exponent Based on the Artificial Rabbit Optimization Algorithm" Water 16, no. 2: 339. https://doi.org/10.3390/w16020339

APA Style

Li, M., Cui, Z., & Fan, T. (2024). The Flood Simulation of the Modified Muskingum Model with a Variable Exponent Based on the Artificial Rabbit Optimization Algorithm. Water, 16(2), 339. https://doi.org/10.3390/w16020339

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