Next Article in Journal
Research on Cascade Reservoirs’ Short-Term Optimal Operation under the Effect of Reverse Regulation
Next Article in Special Issue
Flood Hydrograph Prediction Using Machine Learning Methods
Previous Article in Journal
Water Quality Prediction Model of a Water Diversion Project Based on the Improved Artificial Bee Colony–Backpropagation Neural Network
Previous Article in Special Issue
Physical Hybrid Neural Network Model to Forecast Typhoon Floods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improving the Muskingum Flood Routing Method Using a Hybrid of Particle Swarm Optimization and Bat Algorithm

1
Department of Water Engineering and Hydraulic Structures, Faculty of Civil Engineering, Semnan University, Semnan 35131-19111, Iran
2
Department of Civil Engineering, Faculty of Engineering, University of Malaya, Kuala Lumpur 50603, Malaysia
3
Sustainable Developments in Civil Engineering Research Group, Faculty of Civil Engineering, Ton Duc Thang University, Ho Chi Minh City, Vietnam
4
Civil and Structural Engineering Department, Faculty of Engineering and Built Environment, Universiti Kebangsaan Malaysia, Bangi 43600, Malaysia
5
Department of Civil Engineering, College of Engineering, University Tenaga National, Kajang 43000, Malaysia
6
Institute of Policy Energy Research (IPERe), Universiti Tenaga National, Kajang 43000, Malaysia
7
Institute for Energy Infrastructures, University Tenaga National, Kajang 43000, Malaysia
8
Faculty of Civil Engineering, Universiti Teknologi Malaysia, Johor Bahru 81310, Malaysia
9
Department of Biological and Agricultural Engineering and Zachry Department of Civil Engineering, Texas A&M University, 2117 TAMU, College Station, TX 77843-2117, USA
*
Author to whom correspondence should be addressed.
Water 2018, 10(6), 807; https://doi.org/10.3390/w10060807
Submission received: 21 May 2018 / Revised: 12 June 2018 / Accepted: 15 June 2018 / Published: 19 June 2018
(This article belongs to the Special Issue Flood Forecasting Using Machine Learning Methods)

Abstract

:
Flood prediction and control are among the major tools for decision makers and water resources planners to avoid flood disasters. The Muskingum model is one of the most widely used methods for flood routing prediction. The Muskingum model contains four parameters that must be determined for accurate flood routing. In this context, an optimization process that self-searches for the optimal values of these four parameters might improve the traditional Muskingum model. In this study, a hybrid of the bat algorithm (BA) and the particle swarm optimization (PSO) algorithm, i.e., the hybrid bat-swarm algorithm (HBSA), was developed for the optimal determination of these four parameters. Data for the three different case studies from the USA and the UK were utilized to examine the suitability of the proposed HBSA for flood routing. Comparative analyses based on the sum of squared deviations (SSD), sum of absolute deviations (SAD), error of peak discharge, and error of time to peak showed that the proposed HBSA based on the Muskingum model achieved excellent flood routing accuracy compared to that of other methods while requiring less computational time.

1. Introduction

Floods cause huge economic and social effects on the surrounding environment [1,2], such as breaking levees [3], inundating houses, disrupting transportation systems [4], damaging crops and eroding fertile lands [5]. Thus, flood prediction and flood control are important issues for policy makers and designers [6,7]. Hydrological and hydraulic models are used for flood prediction. The estimation of flood discharge hydrograph at a downstream location given the discharge hydrograph upstream is known as flood routing [8,9]. Hydraulic models based on numerical methods can be used for flood routing, but they involve complex unsteady flow equations [10]. Hydrological models use the spatially lumped continuity equation and a storage equation for flood routing. These models need a small amount of data to predict floods [11]. The Muskingum model is an important hydrological model for flood routing. This model has multiple parameters that should be obtained to accurately predict floods [12], and different versions of the model have been applied for flood routing. One strategy is the optimization method by which the parameters are computed as decision variables [13]. The evolutionary algorithms exhibit a high degree of ability for solving complex optimization problems. These algorithms are efficient, accurate, and flexible [14]. Thus, parameters of the Muskingum model can be computed using evolutionary algorithms.

1.1. Background

Luo and Xie [15] applied the clonal selection algorithm (CLA) to flood routing in China, and CLA decreased the sum of absolute deviations (SAD) by approximately 20% compared to that of the genetic algorithm (GA). The difference between the simulated peak discharge and the observed discharge was small, approximately 0.5 m3/s. Using the initial values of the Muskingum parameters as decision variables, CLA optimized the objective function to obtain the best parameter values.
The Nelder-Mead simplex algorithm was used for flood routing in the USA [16]. This method obtained the best values for the 3 Muskingum model parameters. Thus, the peak discharges were predicted to be better than those of other nonlinear programming methods.
Barati et al. [17] applied a dynamic wave method for flood routing, and the results indicated that some of the parameters, such as Manning roughness or bed slope, were effective. Thus, accurate determination of these parameters is important for flood routing equations.
Karahan et al. [18] applied hybrid harmony search (HS) and particle swarm optimization (PSO) to estimate the Muskingum parameters and indicated that the computational time for the new hybrid model was less than those of PSO and HS. In addition, peak discharge was estimated with the least difference with the observed discharge.
Fallah-Mehdipour [19] applied PSO, GA and nonlinear programming methods for flood routing and indicated that PSO simulated the flood discharges with lower SAD and sum of squared deviation (SSD) values than those of other algorithms. In PSO, the Muskingum parameters were considered as decision variables.
Easa [20] introduced a four-parameter nonlinear Muskingum model for flood routing and showed that the four-parameter model estimated flood discharge better than the three-parameter model did. The GA was used to obtain the optimal values of the Muskingum parameters.
Nelder-simplex and hybrid PSO were evaluated as new methods to estimate Muskingum parameters [21]. The results showed that the new method decreased the computational time compared to that of the simple PSO and GA. Additionally, there was a small difference between the estimated discharge and the observed discharge.
Honey bee mating (HBM) optimization was used for flood routing by a three-parameter nonlinear Muskingum model for one flood in the USA [22]. Results indicated that the SSQ (sum of squared deviations) and SAD values of the model significantly decreased compared to those of the GA and PSO methods. Additionally, the time of peak discharge was predicted well. The convergence of HBM was faster than those of GA and PSO. A new charged search system (CSS) and PSO hybrid method was used for flood routing based on two- and three-parameter nonlinear Muskingum models [23]. Results indicated that the new evolutionary hybrid algorithm needed a sensitivity analysis for the accurate determination of parameters, and the three-parameter model predicted peak discharge more accurately than GA and HS. Additionally, basing PSO and CSS on more population diversity increased the convergence speed of the algorithm.
A nonlinear Muskingum model that considers lateral flow was used for flood routing [24]. This model used the cuckoo algorithm and predicted peak discharge more accurately than other models. The shuffled frog leaping algorithm (SFLA) based on a three-parameter Muskingum model was used for flood routing [25]. Results showed that the SSD and SAD values for the SFLA were lower than those of the GA, PSO and nonlinear programming methods. Additionally, the correlation coefficient based on the SFA between the observed and estimated discharges was more than those of the GA, PSO and nonlinear programming methods.
A new version of the Muskingum model with nine parameters was used for flood routing in another study [26]. The Karun River in Iran was considered as an important case study for the research. The study showed that the results were better than those of previous Muskingum models, but the new model needed more computational time.
A real-coded adaptive GA was used for flood routing, considering lateral flow along the river reach [14]. The study considered a river in China with important floods. Results indicated that the three-parameter nonlinear Muskingum model based on adaptive GA simulated flood discharges, yielded the least values of SSQ and SAD.
Gene-expression programming (GEP) and the weed algorithm (WA) were used for flood routing [1]. Some important floods in the USA were considered as case studies for the study. The results indicated that GEP had a convergence speed that was approximately 100 times higher than that of WA; in addition, the computed SSD for the estimated discharges was lower than that of WA.
The literature review shows that the evolutionary algorithms have high ability for predicting flood routing, but there are some limitations for these algorithms; they all use the same procedure within the Muskingum model for flood routing. The Muskingum models have parameters with unknown values, and these parameters are applied to the algorithms as decision variables. Computation of the objective function can show the best value for each parameter. In fact, the parameters are input as the initial population to the algorithms, and then, the different operators are applied based on the process of each algorithm [1,10,22]. Basically, there are two major limitations: First, there is no assurance that the computed value for Muskingum’s parameters are optimally achieved, second, the convergence rate was relatively high as the required computed values are four parameters.
For example, some algorithms such as the GA can become trapped in the local optimums and cannot compute the best values for the parameters [10]. Some literature reviews have reported that some algorithms have immature solutions because the convergence process happens abruptly as with PSO [21]. Some algorithms (e.g., GA [10]) require more computational time for the convergence process, and some algorithms such as Anti Bee Colony and Shark Algorithm need to accurately determine many random parameters, which leads to complex processes for optimization. Thus, the presentation of a better method that could hybridize the advantages of two different methods and hence has the potential to overcome those drawbacks and achieve better results with high convergence rate is necessary. The present study attempts to develop one of the known algorithms based on a hybrid process for the flood routing. The Bat Algorithm (BA) is developed based on the motivation that the algorithm has problems such as trapping in the local optimums and slow convergence, and the algorithm is based on the high ability of bats to perceive sounds [21]. Thus, the PSO is suggested for the hybrid process to improve the BA, and the hybrid process prevents the abrupt convergence of the PSO. The next section presents the innovation and objectives of this study. Then, the optimization methods are explained, and the results and conclusions for three case studies are explained in the following sections.

1.2. Innovation and Objectives

The bat algorithm (BA), as an optimization method, is based on living bats and the powerful ability of bats to receive sounds from their surroundings. It is used widely in different fields of image processing [27], data sensing systems [28], the determination of the seismic safety of structures [29], the design of wireless sensors [30], and water resource management [31]. However, the algorithm has some weaknesses, such as the probability of being trapped in the local optimums and slow convergence in some complex engineering problems, so it is necessary to modify the BA. This paper reports on a hybrid algorithm (HA) based on PSO and BA. The new HA substitutes the weaker BA solution with the best PSO solution, which can prevent trapping in the local optimums and increase the convergence speed. Additionally, with the hybrid model, the BA population diversity and exploration ability increase as the optimal solution is obtained. These improvements are the motivations underlying the introduction of the new HA in this paper. The new HA was used for flood routing by a four-parameter Muskingum model in two case studies. These two case studies have been investigated widely as benchmarks in the literature, so there is comprehensive information for the comparison of the new HA with other evolutionary algorithms. This study used the HA to find the optimal values of the four parameters of the Muskingum model, extracted the output hydrographs for two case studies, and then compared the results with those generated by other types of Muskingum models and other evolutionary algorithms used in previous studies.

2. Methods

2.1. Muskingum Model

The linear Muskingum model is based on the continuity equation and a storage equation [32,33,34]:
d s t d t = I t O t
S t = K [ x I t + ( 1 x ) O t ]
where S t is the storage (L3); I t is the inflow (L3·T−1); O t is the outflow (L3·T−1); K is the storage time constant, which varies from 0 to 30 (T); and x is the weighting factor, which varies from 0 to 0.5. Previous studies have shown that the linear Muskingum model does not perform well for some rivers; thus, nonlinear Muskingum models have been suggested [33]:
S t = K t [ x I t + ( 1 x ) O t ] m
S t = K [ x I t m + ( 1 x ) O t m ]
The models of Equations (3) and (4) have an additional parameter (m), and the dimension of K is (L3(1−m)·Tm). The current study uses a four-parameter nonlinear Muskingum model based on the model introduced by Easa [20], with a reported high flood routing ability:
S t = K [ x I t α + ( 1 x ) O t α ] m
Equation (5) is similar to Equation (3). However, in Equation (5), parameter m is related to the nonlinear form of the storage equation, while in Equation (3), m is related to the linear form of the storage equation. Equation (5) can cover previous equations; that is, if α = m = 1 , Equation (5) will be the same as Equation (2). If m = 1 ( o r ) α = 1 , Equation (5) covers Equations (3) and (4). The outflow for Equation (3) is computed as follows:
O t = [ ( 1 1 x ) ( S t K ) 1 m ( x 1 x ) ] 1 α
Then, Equation (6) is inserted into Equation (1), and the change in storage is computed based on the following equation:
Δ S t Δ t = I t [ ( 1 1 x ) ( S t K ) 1 m ( x 1 x ) I t α ] 1 α
Finally, the storage for the next step time is computed based on the following equation:
S t + 1 = S t + Δ t

2.2. Bat Algorithm

The BA is based on the powerful echolocation ability of bats, which can generate loud sounds and receive the echoes of the sounds as they return from the surroundings. The BA is based on the following assumptions [31]:
(1)
All bats use echolocation to identify prey and obstacles based on received sound frequencies.
(2)
All bats fly randomly with the velocity (vl) at position (yl), and the frequency, loudness and wavelength values are f l , A 0 and λ , respectively.
(3)
The loudness changes from a large positive (A0) to a small positive value (Amin).
The sounds generated by bats have a pulsation rate (rl) that varies from 0 to 1. The value 1 means that the pulsation rate has reached a maximum value, and 0 means that the pulsation rate has reached a minimum value. The velocity, frequency and position are updated based on the following equations [31]:
f l = f min + ( f max f min ) × β
v l ( t ) = [ y l ( t 1 ) Y ] × f l
y l ( t ) = y l ( t 1 ) + v l ( t ) × t
where f l is the frequency, f min is the minimum frequency, f max is the maximum frequency, Y is the best position for the bats, t is the time step, y l ( t 1 ) is the position of the bats at time t − 1, v l ( t ) is the velocity, and β is the random vector.
A random walk is considered as the local search for the BA:
y ( t ) = y ( t 1 ) + ε A ( t )
where ε is the random value between −1 and 1; and A ( t ) is the loudness.
The pulsation rate and loudness are updated for each level. When the bats find prey, the pulsation rate increases and the loudness decreases for each level. The pulsation rate is updated, based on the following equation:
r l t + 1 = r l 0 [ 1 exp ( γ t ) ] A l t + 1 = α A l t
where α and γ are the constant parameters. The BA is shown in Figure 1.

2.3. Particle Swarm Optimization

For a d-dimensional search space, the population is formed by X i = ( x i 1 , x i 2 , . . , x i D ) T , and the velocity is given by V i = ( v i 1 , v i 2 , . . , v i D ) T . The best previous position is given by P i = ( p i 1 , p i 2 , . . , p i D ) T , and the g index is used for the global solution. The velocity and the position are updated, based on the following equations:
v i d n + 1 = χ [ w v i d n + c 1 r 1 ( p i d n x i d n ) Δ t + c 2 r 2 ( p g d n x i d n ) Δ t ]
x i d n + 1 = x i d n + Δ t v i d n + 1
where d is the number of dimensions, χ is the constriction coefficient, N is the size of the swarm, w is the inertial weight, c1 and c2 are the acceleration coefficients, r1 and r2 are the random numbers, Δ t is the time step, x i d n + 1 is the new position, v i d n + 1 is the new velocity vector, and n is time index.
First, the random population of the swarm is initialized, and then the objective function is computed to determine the local and global leaders. Finally, the velocity and the position are updated, and the cycle continues until the convergence criterion is satisfied.

2.4. Hybrid PSO and BA

The new HA acts based on a communication strategy between two algorithms. The idea of this strategy is that the parallel performance of the two algorithms allows the weaker solutions of one algorithm to be substituted with those of the other algorithm. The parallel structure has subgroups that are based on the division of the population. The subgroups act independently for each iteration. Thus, the objective function is computed for each group, and the weaker solutions of each group are selected and substituted with the better solutions of the subgroups of the other algorithm. The total number of iterations is R, and the total populations is N. Moreover, N1 and N2 are equal, i.e., N/2. Figure 2 shows the performance of the HA.
The algorithm acts based on the following levels:
(1)
The random parameters for both algorithms (PSO + BA) are initialized, and the initial populations for the two algorithms are considered.
(2)
The first initial values for the hydrological parameters (K, x, m and α ) are considered at the start of the algorithm.
(3)
The variation in storage is computed based on Equation (7). The initial outflow is the same as inflow.
(4)
The accumulated storage is computed based on Equation (8).
(5)
The outflow is computed based on Equation (6).
(6)
The time step is compared with the total flood time. If it is less than the total time, the algorithm goes to step 3; otherwise, the algorithm goes to the next level.
(7)
The objective function is computed for the two algorithms and all members that can be seen in the algorithms.
(8)
The velocity and position for the PSO algorithm are updated based on Equations (14) and (15), and the velocity, frequency and position are updated based on Equations (9)–(11).
(9)
The best particles migrate from the PSO algorithm to the BA, and there is a condition for BA similarity. In fact, the specific number of best members for each algorithm is known and is substituted for the worst solutions of the other algorithm.
(10)
The convergence criteria are considered. If the criteria are satisfied, the algorithm finishes; otherwise, the algorithm returns to the second step.
The following indexes are used for comparison of the different methods [35,36,37,38]:
(1)
The sum of the squared deviations between observed and estimated discharges is considered the objective function and is computed based on the following equation:
M i n i m i z e ( S S Q ) = t = 1 n ( O s t O b t ) 2
(2)
The SAD between estimated and observed discharges is computed based on the following equation:
M i n i m i z e ( S A D ) = t = 1 n | O b t O s t |
(3)
The mean absolute error (MARE) between estimated and observed discharges is computed based on the following equation:
M A R E = 1 N i = 1 N | O b t O s t | O b t
(4)
The error for peak discharge (EO) is computed based on the following equation:
E O = | O p e a k , b t O p e a k , s t | O p e a k , b t
(5)
The error for peak time is computed based on the following equation:
E T = | T p e a k , b t T p e a k , s t |
where O s t is the simulated discharge, O b t is the observed discharge, T p e a k , b t is the peak time of observed discharge, T p e a k , s t is the peak time of simulated discharge, and N is the number of data.

3. Case Studies

In this section, brief introduction on the selected case studies have been reported. The first one is the Wilson flood event. The Wilson flood was selected for comparative analysis because it has been widely studied in the literature, resulting in a comprehensive body of information for comparison between the new hybrid method and other methods. This is a benchmark experimental problem that was considered by Wilson [33,34,35,36,37]. The time flood is equal to 120 h, and the peak occurs at a time step of 60 h with a value of 85 cm.
On the other hand, the second case study is the Karahan flood event. The Karahan flood, known as a benchmark problem in flood routing, is based on the 1960 flood of the Wye River in the United Kingdom [14,33,34,35,36]. The Wye River from Everwood to Belmont does not have any tributaries, and it has a small lateral flow, which is considered an important problem for flood routing. The flood time is 198 h, and the peak time is at the step time of 102 h.
Finally, the third case study is the Viessman and Lewis flood event. The Viessman and Lewis [39] multi-peak flood hydrograph was also selected as a benchmark study in this research.

4. Results and Discussion

4.1. Wilson Flood

4.1.1. Sensitivity Analysis for Different Algorithms for the Wilson Flood

Table 1, Table 2 and Table 3 show the inflow and outflow values calculated by different methods. The evolutionary algorithm parameters do not have specific values at the start of the algorithm. Thus, sensitivity analysis was used to obtain the parameters. The objective function selected for this study was SSQ, and the variation in the objective function value was computed for various values of parameters. The best value for each parameter was selected when the objective function value was the minimum. For example, the population size for the HA varied from 20 to 80. When the size was 60, the SSQ value was 4.234 and led to the lowest objective function value. Further, the maximum frequency varied from 3 to 9. When the maximum frequency was 7 Hz, the objective function had the best value. The minimum frequency varied from 0.1 to 0.4, and the best value for the minimum frequency was 0.2. In addition, the maximum loudness varied from 0.2 to 0.8, and the best value based on the objective function value was 0.6 dB. The other parameters are listed in Table 1. For example, the acceleration coefficient varied from 1.6 to 2.2. The best value for this parameter was 2, because the SSQ, which was 4.233, was the lowest. Further, the inertia coefficient varied from 0.3 to 0.9. The best value was 0.7 because the SSQ was 4.234.

4.1.2. Ten Random Results for Different Algorithms for Wilson Flood

Table 4 shows 10 random results for different algorithms. The average solution for the HA was 4.234, while the average solutions were 5.554 and 5.342 for the PSO algorithm and BA, respectively. Thus, the objective function or error value decreased by 23.71% and 20.7% compared to those of the PSO algorithm and BA, respectively. The computational time for the HA was 20 s, while it was 27 s and 25 s for the PSO algorithm and BA; thus, the computational time decreased by 25% and 20% in the HA compared to those in the PSO algorithm and BA, respectively. Additionally, the coefficient of variation for the HA was smaller than that for the PSO algorithm and BA and was therefore nonsignificant, which implied that the one-run program for the HA was reliable. Figure 3 shows the convergence method for the three algorithms, and it is apparent that the HA converged faster than the PSO algorithm and BA did. In fact, the researchers present the 10 or 15 random solutions to show the variation of solutions. It is apparent that the all runs for the HA are close to the average solution and that all runs are reliable as the final solution. Although the BA and PSO have the solutions that are closest to the average solution, the other factors are important for the selection of the best methods. Table 4 shows that the hybrid method based on soft computing is the best method because the minimum, maximum and average solutions for the 10 results are close to each other. Additionally, the computation time for the HA is less than other methods (based on a PC with an i5 2.4 GHz CPU and 500 GB RAM), and the average of the error objective function based on HA is less than those of the other two methods.

4.1.3. Discussion of the Wilson Flood Results

Table 5 shows the performance of the HA, PSO algorithm and BA based on the 4PMM (four-parameter Muskingum model). The computed SSQ based on the HA was 4.234, and the HA method decreased the SSQ by 23.7% and 20.74% compared to those of the PSO and BA methods, respectively. The SAD value based on the HA decreased by 25% and 24% compared to those of the PSO algorithm and BA, respectively. The error between the simulated discharge and the observed discharge was shown by the EO, and the minimum value, 0.0111, was related to the HA. Additionally, the three methods predicted the peak time well, with an ET equal to zero. Table 6 compares the performances of different models based on the current research and previous studies. The performance of the HA based on the 4PMM was better than that based on the 3PMM because the error index for the 4PMM HA was lower than that for the 3PMM HA. For example, SSQ and SAD for the 4PMM were 4.234 and 3.125, respectively, while they were 12.25 and 10.95 for the 3PMM. Table 6 shows the performance of GA, HS, and the imperial competitive algorithm (ICA) based on the literature. The highest value of SSQ was related to GA based on the 3PMM; thus, this model had the worst performance. The SAD for the ICA was 23.46, which indicated that this model had the worst performance based on SAD and 3PMM (three-parameter Muskingum model). When the performances of the BA and PSO algorithm based on the 3PMM were compared with that based on the 4PMM, the models based on the 4PMM had the best performance because of their smallest index errors. All the models in Table 6 predicted the peak time well, with an ET value equal to zero. Table 7 shows the inflow and outflow and the peak value for the observed discharge, which was 85 cm when the time was 60 h. The results showed that peak discharge based on the HA was closer to the observed value. Table 8 shows the values of different coefficients for the Muskingum models. Generally, the results showed that the HA based on the 4PMM performed better than the other models did, with an SSQ value that decreased by 65%, 72% and 47% compared to those of the BA (3PMM), PSO algorithm (3PMM) and HA (3PMM), respectively. The values for other indexes supported this trend. Figure 3 shows that the HA based on the 4PMM performed better than the PSO and BA based on the 4PMM during the flood. Additionally, the comparison of results with other research shows the superiority of HA compared to the other algorithms. The honey bee mating algorithm was considered for flood routing [34]. SSQ based on 3PMM was 37.451, while the 3PMM and HA have smaller values for SSQ. Thus, HA has a better performance than HBMO.
Thus, all sections for the HA show the superiority of the method compared to the other methods. Although, the Muskingum method based on four parameters has more parameters than the Muskingum model based on three parameters, it is important to have a simulated hydrograph that matches with the observed flood. This section shows that prediction of peak discharge is the first priority in flood hydrographs. The observed peak in the hydrograph is 85 cm, while the simulated peak discharge is 85.011 cm. Clearly, the new method has good performance for this issue. The method is the best method for flood routing among other evolutionary algorithms because the peak time can be computed accurately. Table 5 shows the peak time predicted based on HA so that there is no delay between the computed time based on HA and the real time. The flood-emergency management authorities can make the best decisions based on the importance of the different projects.

4.2. Karahan Flood

4.2.1. Discussion of the Karahan Results

The parameters used for the HA were maximum frequency, minimum frequency, maximum loudness, initial population, acceleration coefficient, and inertia coefficient; which were computed based on sensitivity analysis described in the previous section (maximum frequency: 7; minimum frequency: 3; maximum loudness: 0.95; initial population: 50; acceleration coefficient: 2; and inertial weight: 0.7).
Table 9 compares the performances of different algorithms based on the 4PMM for the Karahan flood. The SSQ of the HA and the 4PMM was 30,235, and the SSQ of the HA decreased by 5.85% and 2.8% compared to those of the PSO algorithm and BA, respectively. The SAD for the HA based on the 4PMM was 625, which were 0.32% and 7.54% lower than those of the PSO algorithm and BA, respectively. The difference between the simulated peak discharge and the observed discharge was nonsignificant; the EO of the HA equaled 0.101, a value that was 7.3% and 6.4% lower than those of the PSO algorithm and BA. The best performance based on MARE was related to the HA, which produced the lowest MARE value. Additionally, the three methods predicted the peak discharge well, with very low EO values. The GA had the worst performance of the methods, as shown in Table 10. The SSQ for the HS method was 39,944, which was the highest SSQ value produced by the models based on the 4PMM and 3PMM. The HA based on the 3PMM performed better than the other methods based on the 3PMM, as the SSQ of the HA decreased by 11.38%, 18%, 7.1%, 3.4% and 6.3% compared to those of the GA, HS, ICA, BA and PSO algorithm, respectively. The HA based on the 3PMM produced the lowest SAD value of all the 3PMM-based models. The performances of the HA, PSO algorithm and BA based on the 4PMM were comparable to the corresponding performances based on the 3PMM, but the 4PMM algorithm performed better than did the 3PMM algorithm. For example, the SSQ values of the HA, PSO algorithm and BA based on the 4PMM were 2.8%, 3.4% and 0.27% lower than those for the HA, PSO algorithm and BA based on the 3PMM, respectively. Table 11 shows the inflow and outflow for different methods and the peak discharge, which was 830 cm when the step time was 96 h. These results show that the discharges estimated by the HA were near to the observed discharges. Table 12 shows the extracted coefficients for different algorithms. Figure 4 shows the HA based on the 4PMM had the best performance during the flood, with a good match between the observed discharges and the discharges estimated by the HA. The general results for this section show that the HA based on 3PMM and 4PMM can predict the peak time well such that the ET parameter for HA based on 3PMM and 4PMM equals zero. Another point is related to the value of predicted discharge such that the EO index shows the ability of the different methods. 4PMM models based on HA, PSO and BA have a better performance: If the EO based on the table equals to 0.701, the EO based on Table 9 and FPMM models is less than 0.701. Although the 3PMM needs less parameters, the time and value peak can be predicted accurately based on 4PMM. Clearly, the importance of the model accuracy to decision makers can be considered in the selection of 3PMM and 4PMM.

4.2.2. Ten Random Results for Karahan Flood

Table 13 shows the averages of 10 random results for different methods. The average solution of the HA was lower than that of the PSO algorithm and BA. Additionally, the computational time for the HA was 17% and 29% lower than those for the BA and PSO algorithm, respectively. Additionally, the coefficient of variation for the HA was a small value, which proved that the HA based on one computer program run can be reliable, producing high-quality solutions. Table 13 helps determine the best decision based on time, value of objective function and variation coefficient. The least value for the error objective function is provided by the HA based on the least value of the error objective function. Additionally, the least probable time is related to the HA, and all the runs for the HA have the smallest probable variations. Thus, the method encourages the decision makers to select HA based on the mentioned indexes. Additionally, Table 14 shows the sensitivity analysis for the HA and 4PMM model. The best population size for the HA is 60, and the best value for the maximum frequency and minimum frequency is 7 and 3. Other parameters can be seen in Table 14.

4.3. Discussion of the Viessman and Lewis Flood Results

The Viessman and Lewis Fi [39] multi-peak flood hydrograph was selected for this section. Table 15 shows the performances of different models in analysing the Viessman and Lewis flood. The results were compared with those of the WA in the literature, and this is a synthetic problem [33,34,35,36,37]. The highest SSQ for this flood, 73,312, was produced by the WA based on the 3PMM. The HA based on the 4PMM performed better than did the other models. The SSQ, SAD and MARE values of the HA based on the 4PMM had the lowest values. The BA based on the 4PMM performed better than did the WA and PSO algorithm based on the 3PMM. For example, the SSQ of the BA based on the 4PMM was 47,224, which was 14%, 3.8%, 16.72% and 35% lower than those of the PSO algorithm (4PMM), BA (3PMM), PSO algorithm (3PMM) and WA (3PMM), respectively. Additionally, the HA, PSO algorithm and BA based on the 4PMM outperformed the PSO algorithm, BA and HA based on the 3PMM. Figure 5 shows that the estimated discharge of the HA based on the 4PMM matched the observed discharge during the flood well. Furthermore, Table 16 shows the sensitivity analysis for the 4PMM and HA. The best size population for the HA is 60 because it has the least value for the objective function. The maximum loudness for the HA is 0.60, and the other parameters can be seen in Table 16. The comparison of results with other research studies shows the superiority of 4PMM and HA. For example, one study considered flood routing based on improved PSO based on correction of Wight inertia and a Muskingum flood with two parameters. The MARE index is 0.911, while the MARE for the HA and 4PMM is 0.794. HA acts better than the improved PSO and the Muskingum model with two parameters [34].
Moreover, Table 17 shows the 10 random results for different methods, and it can be seen that the value of the error objective function for the 4PMM and HA is less than PSO and BA. Additionally, the computational time based on HA and 4PMM is less than those of other methods. The solution for the HA is high quality because the variation coefficient is small compared to those of the other methods.

5. Conclusions

The present study considers flood routing based on the new hybrid BA and PSO algorithm. The new algorithm is based on the substitution of the weaker BA solutions with PSO. The structure of the new HA, with the elimination of the weaker BA solutions, enables the escape from local optimums. Three case studies that are benchmark studies in the flood routing field were used in this study, and the 4PMM was selected for analysis. For the Wilson flood case study, the computational time of the HA method was 25% and 20% lower than those of the PSO algorithm and BA methods, respectively. The SSQ computed by the HA was 4.234, which was 23.7% and 20.74% lower than those computed by the PSO algorithm and BA, respectively. Generally, the Wilson flood results showed that the HA based on the 4PMM performed better than the other models, with an SSQ value 65%, 72% and 47% lower than those of the BA (3PMM), PSO algorithm (3PMM) and HA (3PMM), respectively; the values of the other indexes agreed with this result. The results for the Karahan flood case study showed that the performances of different models based on the BA, PSO algorithm and HA with the 4PMM were better than those of the same models with the 3PMM. In the Karahan flood case study, the HA based on the 3PMM performed better than the other methods based on the 3PMM did, with an SSQ value that was 11.38%, 18%, 7.1%, 3.4%, and 6.3% lower than those of the GA, HS, ICA, BA and PSO methods, respectively. For the Viessman and Lewis flood case study, the SSQ, SAD and MARE of the HA based on the 4PMM were the lowest. The HA based on the 4PMM performed better than did the BA and PSO based on the 3PMM. The results indicated that the HA based on substitution of weaker solutions of each algorithm with the strong solutions can decrease the computational time, and the chance of obtaining the best solutions increased significantly. With the new strategy, there is no problem with trapping in the local optimums because the elimination of weaker solutions can lead to exits from the local optimums.
However, the HA based on the 4PMM had the best performance of all the models, suggesting that it should be used in future studies with other advanced Muskingum models and more parameters, such as 7- or 9-parameter Muskingum models, to evaluate the skill of HA for flood routing prediction modelling.

Author Contributions

A.E.-S. and V.P.S. initiated the research point of the hydrological problem and supervised the study. M.E. collected the data and performed the modeling. M.B.A.M., A.N.A. and M.F.A. handled the writing up of the introduction and methodology. The analysis of the modeling outcomes have been handled by F.B.O., S.S., Z.M.Y., H.A.A. and M.E.; Z.M.Y. and H.A.A. organized the whole manuscript and managed the paper submission and revision.

Funding

This research was funded by University of Malaya research grants RP017C-15SUS and RP025A-18SUS.

Acknowledgments

The authors appreciate the technical and financial support received from the University of Malaya.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bagatur, T.; Onen, F. Development of predictive model for flood routing using genetic expression programming. J. Flood Risk Manag. 2018, 11, S444–S454. [Google Scholar] [CrossRef]
  2. El-shafie, A.; Noureldin, A.; Taha, M.; Hussain, A.; Mukhlisin, M. Dynamic versus static neural network model for rainfall forecasting at Klang River Basin, Malaysia. Hydrol. Earth Syst. Sci. 2012, 16, 1151–1169. [Google Scholar] [CrossRef] [Green Version]
  3. Mazzoleni, M.; Noh, S.J.; Lee, H.; Liu, Y.; Seo, D.J.; Amaranto, A.; Alfonso, L.; Solomatine, D.P. Real-time assimilation of streamflow observations into a hydrological routing model: Effects of model structures and updating methods. Hydrol. Sci. J. 2018, 63, 386–407. [Google Scholar] [CrossRef]
  4. Zhou, Y.; Guo, S.; Chang, F.J.; 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]
  5. El-shafie, A.; Mukhlisin, M.; Najah, A.A.; Taha, M.R. Performance of artificial neural network and regression techniques for rainfall-runoff prediction. Int. J. Phys. Sci. 2011, 6, 1997–2003. [Google Scholar] [CrossRef]
  6. Vatankhah, A.R. Discussion of “Assessment of Modified Honey Bee Mating Optimization for Parameter Estimation of Nonlinear Muskingum Models” by Majid Niazkar and Seied Hosein Afzali. J. Hydrol. Eng. 2018, 23, 07018002. [Google Scholar] [CrossRef]
  7. Yaseen, Z.M.; El-shafie, A.; Jaafar, O.; Afan, H.A.; Sayl, K.N. Artificial intelligence based models for stream-flow forecasting: 2000–2015. J. Hydrol. 2015, 530, 829–844. [Google Scholar] [CrossRef]
  8. Pati, A.; Sen, S.; Perumal, M. Modified Channel-Routing Scheme for SWAT Model. J. Hydrol. Eng. 2018, 23, 1–8. [Google Scholar] [CrossRef]
  9. Yaseen, Z.M.; Jaafar, O.; Deo, R.C.; Kisi, O.; Adamowski, J.; Quilty, J.; El-shafie, A. Stream-flow forecasting using extreme learning machines: A case study in a semi-arid region in Iraq. J. Hydrol. 2016. [Google Scholar] [CrossRef]
  10. Barati, R. Comment on “Evaluation of a physically based quasi-linear and a conceptually based nonlinear Muskingum methods” by Muthiah Perumal, Gokmen Tayfur, C. Madhusudana Rao and Gurhan Gurarslan. J. Hydrol. 2017, 550, 218–219. [Google Scholar] [CrossRef]
  11. Barbetta, S.; Moramarco, T.; Perumal, M. A Muskingum-based methodology for river discharge estimation and rating curve development under significant lateral inflow conditions. J. Hydrol. 2017, 554, 216–232. [Google Scholar] [CrossRef]
  12. Barati, R.; Badfar, M.; Azizyan, G.; Akbari, G.H. Discussion of “Parameter Estimation of Extended Nonlinear Muskingum Models with the Weed Optimization Algorithm” by Farzan Hamedi, Omid Bozorg-Haddad, Maryam Pazoki, Hamid-Reza Asgari, Mehran Parsa and Hugo A. Loáiciga. J. Irrig. Drain. Eng. 2018, 144, 07017021. [Google Scholar] [CrossRef]
  13. Yoo, C.; Lee, J.; Lee, M. Parameter Estimation of the Muskingum Channel Flood-Routing Model in Ungauged Channel Reaches. J. Hydrol. Eng. 2017, 22, 05017005. [Google Scholar] [CrossRef]
  14. Zhang, S.; Kang, L.; Zhou, B. Parameter estimation of nonlinear Muskingum model with variable exponent using adaptive genetic algorithm. In Environmental Conservation, Clean Water, Air & Soil (CleanWAS); IWA Publishing: London, UK, 2017; pp. 231–237. [Google Scholar]
  15. Luo, J.; Xie, J. Parameter Estimation for Nonlinear Muskingum Model Based on Immune Clonal Selection Algorithm. J. Hydrol. Eng. 2010, 15, 844–851. [Google Scholar] [CrossRef]
  16. Barati, R. Parameter Estimation of Nonlinear Muskingum Models Using Nelder-Mead Simplex Algorithm. J. Hydrol. Eng. 2011, 16, 946–954. [Google Scholar] [CrossRef]
  17. Barati, R.; Rahimi, S.; Akbari, G.H. Analysis of dynamic wave model for flood routing in natural rivers. Water Sci. Eng. 2012, 5, 243–258. [Google Scholar] [CrossRef]
  18. Karahan, H.; Gurarslan, G.; Geem, Z.W. Parameter Estimation of the Nonlinear Muskingum Flood-Routing Model Using a Hybrid Harmony Search Algorithm. J. Hydrol. Eng. 2013, 18, 352–360. [Google Scholar] [CrossRef]
  19. Orouji, H.; Haddad, O.B.; Fallah-Mehdipour, E.; Mariño, M.A. Estimation of Muskingum parameter by meta-heuristic algorithms. In Proceedings of the Institution of Civil Engineers: Water Management; Institution of Civil Engineers (ICE): London, UK, 2013; Volume 166, pp. 315–324. [Google Scholar]
  20. 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] [Green Version]
  21. Ouyang, A.; Li, K.; Truong, T.K.; Sallam, A.; Sha, E.H.M. Hybrid particle swarm optimization for parameter estimation of Muskingum model. Neural Comput. Appl. 2014, 25, 1785–1799. [Google Scholar] [CrossRef]
  22. Niazkar, M.; Afzali, S.H. Assessment of modified honey bee mating optimization for parameter estimation of nonlinear Muskingum models. J. Hydrol. Eng. 2015, 20, 04014055. [Google Scholar] [CrossRef]
  23. Talatahari, S.; Sheikholeslami, R.; Farahmand Azar, B.; Daneshpajouh, H. Optimal Parameter Estimation for Muskingum Model Using a CSS-PSO Method. Adv. Mech. Eng. 2013, 5, 480954. [Google Scholar] [CrossRef]
  24. Karahan, H.; Gurarslan, G.; Geem, Z.W. A new nonlinear Muskingum flood routing model incorporating lateral flow. Eng. Optim. 2015, 47, 737–749. [Google Scholar] [CrossRef]
  25. Bozorg Haddad, O.; Hamedi, F.; Orouji, H.; Pazoki, M.; Loáiciga, H.A. A Re-Parameterized and Improved Nonlinear Muskingum Model for Flood Routing. Water Resour. Manag. 2015, 29, 3419–3440. [Google Scholar] [CrossRef]
  26. Niazkar, M.; Afzali, S.H. Parameter estimation of new improved Muskingum model using new hybrid model. Hydrol. Res. 2016, 49, nh2016089. [Google Scholar] [CrossRef]
  27. Dhar, S.; Kundu, M.K. A novel method for image thresholding using interval type-2 fuzzy set and Bat algorithm. Appl. Soft Comput. 2018, 63, 1–37. [Google Scholar] [CrossRef]
  28. Cui, Z.; Cao, Y.; Cai, X.; Cai, J.; Chen, J. Optimal LEACH protocol with modified bat algorithm for big data sensing systems in Internet of Things. J. Parallel Distrib. Comput. 2018. [Google Scholar] [CrossRef]
  29. Bekdaş, G.; Nigdeli, S.M.; Yang, X.S. A novel bat algorithm based optimum tuning of mass dampers for improving the seismic safety of structures. Eng. Struct. 2018, 159, 89–98. [Google Scholar] [CrossRef]
  30. Mihoubi, M.; Rahmoun, A.; Lorenz, P.; Lasla, N. An effective Bat algorithm for node localization in distributed wireless sensor network. Secur. Priv. 2017, 1, e7. [Google Scholar] [CrossRef]
  31. Bozorg-Haddad, O.; Karimirad, I.; Seifollahi-Aghmiuni, S.; Loáiciga, H.A. Development and Application of the Bat Algorithm for Optimizing the Operation of Reservoir Systems. J. Water Resour. Plan. Manag. 2015, 141, 04014097. [Google Scholar] [CrossRef]
  32. Gill, M.A. Flood routing by the Muskingum method. J. Hydrol. 1978, 36, 353–363. [Google Scholar] [CrossRef]
  33. Tung, Y. River Flood Routing by Nonlinear Muskingum Method. J. Hydraul. Eng. 1985, 111, 1447–1460. [Google Scholar] [CrossRef]
  34. Geem, Z.; Kim, J.; Loganathan, G. A new heuristic optimization algorithm: Harmony search. Simulation 2001, 76, 60–68. [Google Scholar] [CrossRef]
  35. Elzwayie, A.; El-shafie, A.; Yaseen, Z.M.; Afan, H.A.; Allawi, M.F. RBFNN-based model for heavy metal prediction for different climatic and pollution conditions. Neural Comput. Appl. 2017, 28, 1991–2003. [Google Scholar] [CrossRef]
  36. Yaseen, Z.M.; Fu, M.; Wang, C.; Hanna, W.; Wan, M.; Deo, R.C.; El-shafie, A. Application of the Hybrid Artificial Neural Network Coupled with Rolling Mechanism and Grey Model Algorithms for Streamflow Forecasting Over Multiple Time Horizons. Water Res. Manag. 2018, 32, 1883–1899. [Google Scholar] [CrossRef]
  37. Diop, L.; Bodian, A.; Djaman, K.; Yaseen, Z.M.; Deo, R.C.; El-shafie, A.; Brown, L.C. The influence of climatic inputs on stream-flow pattern forecasting: Case study of Upper Senegal River. Environ. Earth Sci. 2018, 77, 182. [Google Scholar] [CrossRef]
  38. Ghorbani, M.A.; Deo, R.C.; Yaseen, Z.M.; Kashani, M.H. Pan evaporation prediction using a hybrid multilayer perceptron-firefly algorithm (MLP-FFA) model : Case study in North Iran. Theor. Appl. Climatol. 2017, 1–13. [Google Scholar] [CrossRef]
  39. Warrren, V.; Gary, L.L. Introduction to Hydrology; Prentice Hall India (P) Ltd.: Delhi, India, 2003. [Google Scholar]
  40. Moghaddam, A.; Behmanesh, J.; Farsijani, A. Parameters Estimation for the New Four-Parameter Nonlinear Muskingum Model Using the Particle Swarm Optimization. Water Resour. Manag. 2016, 30, 2143–2160. [Google Scholar] [CrossRef]
Figure 1. Bat Algorithm procedure.
Figure 1. Bat Algorithm procedure.
Water 10 00807 g001
Figure 2. Hybrid algorithm procedure.
Figure 2. Hybrid algorithm procedure.
Water 10 00807 g002
Figure 3. Simulated discharges based on different methods for the Wilson flood.
Figure 3. Simulated discharges based on different methods for the Wilson flood.
Water 10 00807 g003
Figure 4. Simulated discharges by different algorithms for Karahan flood.
Figure 4. Simulated discharges by different algorithms for Karahan flood.
Water 10 00807 g004
Figure 5. Extracted hydrograph for Viessman and Lewis flood based on 4PMM.
Figure 5. Extracted hydrograph for Viessman and Lewis flood based on 4PMM.
Water 10 00807 g005
Table 1. Sensitivity analysis for hybrid algorithm.
Table 1. Sensitivity analysis for hybrid algorithm.
Population SizeObjective FunctionMaximum Frequency (Hz)Objective FunctionMinimum Frequency (Hz)Objective FunctionMaximum Loudness (dB)Objective Functionc1Objective Functionc2Objective FunctionwObjective Function
205.12335.25415.5650.204.9991.65.1451.65.0110.35.133
404.78954.88424.9870.404.8451.84.9331.84.9870.54.654
604.23474.23334.2340.604.2342.04.2342.04.2340.704.235
804.31294.67944.7890.804.5652.24.5552.24.4450.904.512
Table 2. Sensitivity analysis for PSO.
Table 2. Sensitivity analysis for PSO.
Population SizeObjective Functionc1Objective Functionc2Objective FunctionwObjective Function
205.9811.605.8911.605.9540.305.845
405.7851.805.6541.805.8780.505.764
605.5552.05.5542.05.5540.705.555
705.8942.25.8922.25.8910.905.789
Table 3. Sensitivity analysis for BA.
Table 3. Sensitivity analysis for BA.
Population SizeObjective FunctionMaximum Frequency (Hz)Objective FunctionMinimum Frequency (Hz)Objective FunctionLoudness (dB)Objective Function
205.76535.81215.9110.35.912
405.45555.69125.7830.55.678
605.34275.34235.3430.705.343
705.69495.61145.4550.905.678
Table 4. Ten random results for average solutions for Wilson flood.
Table 4. Ten random results for average solutions for Wilson flood.
Run NumberHABAPSO
14.2345.3425.555
24.2335.3485.555
34.2345.3425.555
44.2345.3425.555
54.2345.3425.559
64.2335.3425.560
74.2345.3425.555
84.2345.3425.555
94.2345.3425.555
104.2345.3425.555
Average4.2345.3425.555
Computational time20 s27 s25 s
Variation coefficient0.000070.00030.0004
Table 5. Comparison of results based on four-parameter Muskingum model for Wilson flood.
Table 5. Comparison of results based on four-parameter Muskingum model for Wilson flood.
MethodSSQSADMAREEOET
HA4.2343.1250.0120.001110
PSO5.5554.1280.0170.002510
BA5.3424.1170.0150.001670
Table 6. Comparison of results for Wilson flood.
Table 6. Comparison of results for Wilson flood.
MethodSSQSADMAREEOET
GA [40]
(Three-parameter Muskingum)
38.23023.000.09120.00830
HS [40]
(Three-parameter Muskingum)
36.78023.400.08780.01070
ICA [40]
(Three-parameter Muskingum)
36.80123.460.07450.01050
BA (current research)
(Three-parameter Muskingum)
12.2510.950.02150.00790
PSO (current research)
(Three-parameter Muskingum)
14.7812.720.03250.00810
HA (current research)
(Three-parameter Muskingum)
8.2156.5150.02050.00430
Table 7. Inflow and outflow for flood routing.
Table 7. Inflow and outflow for flood routing.
TimeInflow (cm)Outflow (Observed-cm)Hybrid Method (cm)BA (cm)PSO
0222222.022.022.0
6232122.023.023.0
12352121.022.523.5
18712625.025.026.0
241033434.035.035.5
301114443.544.044.0
361095554.055.055.5
421006666.067.068.0
48867574.074.075.0
54718281.582.083.0
60598585.001185.0025185.00267
66478484.084.084.0
72398081.080.581.0
78327374.073.074.0
84286464.065.066.0
90245454.055.056.0
96224444.044.045.0
102213636.037.038.0
108203030.531.031.0
114192525.526.226.9
120192223.024.025.0
126181920.021.022.0
Table 8. Extracted coefficients for the four-parameter Muskingum model.
Table 8. Extracted coefficients for the four-parameter Muskingum model.
MethodKxm α
HA0.1640.28793.7810.4678
BA0.1520.27683.5670.4567
PSO0.1440.26453.1230.3789
Table 9. Comparison of results based on the four-parameter Muskingum model for Karahan flood.
Table 9. Comparison of results based on the four-parameter Muskingum model for Karahan flood.
MethodSSQSADMAREEOET
HA30,2356250.220.1010
PSO32,1196970.250.1090
BA31,1126760.240.1080
Table 10. Comparison of results for Karahan flood.
Table 10. Comparison of results for Karahan flood.
MethodSSQSADMAREEOET
GA [40]
(Three-parameter Muskingum)
35,12319800.9100.7010
HS [40]
(Three-parameter Muskingum)
37,94421610.9240.7980
ICA [40]
(Three-parameter Muskingum)
37,82520540.9140.7870
BA (current research)
(Three-parameter Muskingum)
32,2287120.4200.1150
PSO (current research)
(Three-parameter Muskingum)
33,2297350.4540.1250
HA (current research)
(Three-parameter Muskingum)
31,1256970.2540.1050
Table 11. The inflow and outflow for flood routing.
Table 11. The inflow and outflow for flood routing.
TimeInflow (cm)Outflow (Observed-cm)Hybrid Method (cm)BA (cm)PSO (cm)
0154102102.0102.0102.0
6150140139.23138.23154.2
12219169170.21171.24152.1
18182190185.12183.24179.4
24182209202.34200.11190.9
30192218212.23198.23185.4
36165210207.11192.32186.9
42150194192.12189.23180.2
48128172170.21169.24164.1
54168149147.21146.74143.7
60260136137.21139.23152.8
66471228219.21212.23196.3
72717303300.11298.21267.3
781092366358.11354.23351.4
841145456436.32426.73431.8
90600615612.21623.24617.4
96365830830.101830.108830.109
102277969894.12879.12836.70
108227665665.101665.108665.109
114187519519.21523.12549.10
120161444435.68424.32416.90
126143321315.23312.11305.0
132126208210.21212.21221.40
138115176169.21166.24163.38
144102148142.12139.23131.20
15093125119.21115.67110.0
15688114109.21100.2196.40
16282106110.21112.1189.20
168769792.2189.2382.70
174738982.1279.4376.30
180708180.2378.1273.00
186677679.1475.1269.80
192637170.1470.1166.7
198596670.2369.1262.40
Table 12. The extracted coefficients for the four-parameter Muskingum model.
Table 12. The extracted coefficients for the four-parameter Muskingum model.
MethodKxm α
HA0.6100.4043.7811.125
BA0.5780.3112.8961.112
PSO0.5780.3092.7891.105
Table 13. Investigation of different methods for Karahan flood.
Table 13. Investigation of different methods for Karahan flood.
Run NumberHABAPSO
130,23531,11232,119
230,23731,11732,119
330,23531,11232,119
430,23531,11232,119
530,23531,11232,119
630,23531,11232,119
730,23731,11232,122
830,23531,11732,112
930,23531,11232,119
1030,23531,11232,119
Average30,235.431,11332,119
Computational time19 s23 s27 s
Variation coefficient0.000020.000050.00007
Table 14. Sensitivity analysis for Karahan flood.
Table 14. Sensitivity analysis for Karahan flood.
Population SizeObjective FunctionMaximum Frequency (Hz)Objective FunctionMinimum Frequency (Hz)Objective FunctionMaximum Loudness (dB)Objective Functionc1Objective Functionc2Objective FunctionwObjective Function
2034,231333,278132,2780.2031,1241.632,1121.632,1140.331,127
4032,278532,211231,1120.4030,2981.831,2141.831,2890.531,119
6030,235730,235330,2350.6030,2352.030,2352.030,2350.7030,235
8031,112931,265431,1120.8030,2362.231,1122.231,1120.9030,254
Table 15. Evaluation of different methods for Viessman and Lewis flood.
Table 15. Evaluation of different methods for Viessman and Lewis flood.
MethodSSQSADMAREEOET
HA (4PMM)45,225998.240.7940.1110
PSO (4PMM)55,1241012.220.8120.2090
BA (4PMM)47,2241001.140.7980.1180
HA (3PMM)48,2251002.230.8120.1150
PSO (3PMM)56,7121014.450.8670.2880
BA (3PMM)49,1121009.230.7240.2020
WA (3PMM)73,3121037.250.9940.4880
Table 16. Sensitivity analysis for Viessman and Lewis flood.
Table 16. Sensitivity analysis for Viessman and Lewis flood.
Population SizeObjective FunctionMaximum Frequency (Hz)Objective FunctionMinimum Frequency (Hz)Objective FunctionMaximum Loudness (dB)Objective Functionc1Objective Functionc2Objective FunctionwObjective Function
2047,229347,312149,2780.2046,1241.647,1191.648,1240.348,119
4046,214547,001247,1120.4045,2981.846,2241.847,2110.547,015
6045,225745,225345,2250.6045,2252.045,2252.045,2250.7045,225
8049,112945,287448,1120.8047,1192.246,1792.246,1170.9046,119
Table 17. Investigation of different methods for Karahan flood.
Table 17. Investigation of different methods for Karahan flood.
Run NumberHAPSOBA
145,22555,12447,224
245,22655,12447,226
345,22555,12747,224
445,22555,12447,224
545,22555,12447,224
645,22555,12447,224
745,22555,12447,224
845,22555,12447,224
945,22555,12447,224
1045,22555,12447,224
Average45,2253111347,224
Computational time15 s17 s19 s
Variation coefficient0.0000040.0000060.00005

Share and Cite

MDPI and ACS Style

Ehteram, M.; Binti Othman, F.; Mundher Yaseen, Z.; Abdulmohsin Afan, H.; Falah Allawi, M.; Bt. Abdul Malek, M.; Najah Ahmed, A.; Shahid, S.; P. Singh, V.; El-Shafie, A. Improving the Muskingum Flood Routing Method Using a Hybrid of Particle Swarm Optimization and Bat Algorithm. Water 2018, 10, 807. https://doi.org/10.3390/w10060807

AMA Style

Ehteram M, Binti Othman F, Mundher Yaseen Z, Abdulmohsin Afan H, Falah Allawi M, Bt. Abdul Malek M, Najah Ahmed A, Shahid S, P. Singh V, El-Shafie A. Improving the Muskingum Flood Routing Method Using a Hybrid of Particle Swarm Optimization and Bat Algorithm. Water. 2018; 10(6):807. https://doi.org/10.3390/w10060807

Chicago/Turabian Style

Ehteram, Mohammad, Faridah Binti Othman, Zaher Mundher Yaseen, Haitham Abdulmohsin Afan, Mohammed Falah Allawi, Marlinda Bt. Abdul Malek, Ali Najah Ahmed, Shamsuddin Shahid, Vijay P. Singh, and Ahmed El-Shafie. 2018. "Improving the Muskingum Flood Routing Method Using a Hybrid of Particle Swarm Optimization and Bat Algorithm" Water 10, no. 6: 807. https://doi.org/10.3390/w10060807

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