A Novel Approach to Supply the Water Reservoir Demand Based on a Hybrid Whale Optimization Algorithm

Managing water resources requires the optimum operation of dam reservoirs. To satisfy the downstream water demand in the operational optimization of Boostan dam reservoir, the improved whale optimization algorithm (IWOA) performance was compared in the present study with that of its constituents (i.e., the whale optimization and diﬀerential evolution) based on GAMS nonlinear programming results. The model evaluative indicators and an objective function were used to select the optimal algorithm. The ﬁndings suggested that IWOA resulted in the lowest computational duration and fastest convergence rate compared to the other algorithms. Additionally, the average water demand and discharge volume of IWOA were 3.21 × 10 6 m 3 and 3.03 × 10 6 m 3 , respectively. In contrast, the other algorithms yielded lower water release volumes. IWOA enhanced the WOA performance by 21.7% through reducing the variation coeﬃcient by 78% in optimizing the objective function. The water demand was therefore more eﬀectively satisﬁed by the IWOA compared to the other algorithms. Furthermore, the IWOA resulted in a lower amount of errors. The hybrid algorithm performance increased in terms of all the evaluative indicators. Developing multicriteria decision-making models using TOPSIS and the Shannon entropy also suggested the IWOA excels the other algorithms in optimizing the reservoir operational problem.


Introduction
Managing water resources constitutes an essential domain in hydrological sciences for the sustainable utilization and allocation of valuable resources such as water [1]. Managers of water resources have identified water shortage as a globally major problem in various areas.
is practical problem can be solved using state-of-the-art technology, innovative ideas, and strategic decisions made on water utilization [2,3]. Water scarcity can also be solved through the operational optimization of reservoirs. Water resource policymakers therefore require appropriate optimization algorithms for lowering the difference between the demand for water and water discharge and generating optimal operation rules for water usage [4]. A reservoir operation is controlled mainly to reduce water shortage by developing the rules for the controlled discharge of water for downstream users [5]. e minimum permissible volume of the water stored in reservoirs in critical conditions and dry seasons must also be ensured [6]. Appropriate optimization algorithms should be therefore selected to take the necessary measures for the controlled discharge of water, provide various water users, including agricultural and industrial sectors, with a sustainable water supply, and therefore solve the operational problems of reservoirs [7,8].
Metaheuristic and evolutionary optimization algorithms have been developed during the past 2 decades to solve complicated engineering problems. e main advantages of these methods include their simplicity of adjusting system constraints and objective functions within their mathematical domain. Obtaining a globally optimal solution can be, however, time-intensive, which is unacceptable in real-time operation. e computational time required by metaheuristic and evolutionary techniques for solving complicated problems has been demonstrated to be the least, making them an optimal choice for main decision-makers who need real-time solutions for their problems [4].

Review of the Literature
Genetic programming was applied by FallahMehdipour et al. [9] to a reservoir operation to decrease shortages of irrigation by considering the downstream discharge of water, the decision variable, and testing various mathematical and arithmetic operators for genetic programming. ey found the demand for irrigation to be satisfied with the maximum reliability indicator compared to the reliability indicators of genetic algorithm and particle swarm optimization.
A case study was conducted in Iran by Ostadrahimi et al. [10] by applying particle swarm optimization to a multireservoir system to obtain various rule curves and decrease shortages of irrigation. e convergence rate of this new version of PSO was increased by making it function on the basis of a modified inertia weight. e data obtained were compared with those of nonlinear programming and genetic algorithm.
e released water met the irrigation demand using particle swarm optimization, whose convergence was found to be faster than that of GA. e reservoir operation was optimized and the benefits of power generation were improved by Afsahr [5] taking advantage of a modified particle swarm optimization adapted to hydraulic and hydrological boundary conditions of the problem in a way that the algorithm modeled the reservoir operation through a drought. e solution obtained using this version of particle swarm optimization was on average close to the globally optimal solution, and the power generated per year was higher than the power obtained using harmony search and genetic algorithm. Employing ant colony optimization to operationally optimize multireservoir systems helped satisfy the demand for irrigation with the minimum risk and minimize the average irrigation deficit [11]. Afshar and Moeini (2013) increased the diversity of the initial population by developing ant colony optimization on the basis of a mutation operator.
Genetic algorithm and novel dynamic programming proposed by [12] was applied to multireservoir operation to increase the power generated in China. Genetic algorithm outperformed dynamic programming despite the limited power generated in dry years, which required a rationing process. e position of drops was considered the initial population and a decision variable in the water cycle algorithm to reduce deficit in power based on the drop movements in the environment [13]. e running time of the water cycle algorithm was lower than that of genetic algorithm and particle swarm optimization, and this algorithm minimized the average deficit. A case study in Iran by Bozorg-Hadad et al. [14] was aimed at reducing power deficits using the bat algorithm based on surrounding frequencies received by bats. As a decision variable, the position of a bat was updated according to frequency and velocity. e rate of convergence of this algorithm based on random walk and local search was found to be higher than that of particle swarm optimization and genetic algorithm. Moreover, its power generation received a higher index of reliability.
Deficits in supplying demand for irrigation were minimized in a reservoir operation, and the released water obtained from a third-order nonlinear rule curve on the basis of genetic algorithm was found to satisfy irrigation requirements given the high reliability index obtained [15]. Biogeography-based optimization used in a case study in Iran to reduce deficits in power [16] on the basis of migration processes and habitats in the environment resulted in solutions that approached the globally optimal solution through fewer iterations. e simple algorithm process used an initial species population in the environment. Applying the weed optimization algorithm to increase power generation in a multireservoir operation by Asghari et al. [17] showed the superiority of this algorithm by convergence rate and generated benefits compared to the water cycle algorithm and genetic algorithm as well as particle swarm optimization.
Employing fixed length genetic programming to decrease the demand for irrigation in multireservoir operation [18] found this method to contribute to meeting the demand with a larger convergence rate than that of other algorithms and a higher reliability index than that of genetic programming. Utilizing the gravitational search algorithm to decrease power generation in a multireservoir system [19] suggested the need of this algorithm for accurate analyses of sensitivity to its random parameters. In comparison with particle swarm optimization and genetic algorithm, this algorithm resulted in a higher power generated annually and a larger convergence rate. Applying the shark optimization algorithm to reservoir operation demonstrated the effectiveness of this algorithm in satisfying the demand for irrigation with a larger convergence rate and index of reliability than those obtained from particle swarm optimization or genetic algorithm [20]. Applying monarch butterfly optimization by [5] Ehteram et al. [8] to multireservoir operation in China based on a migration operator for the butterflies showed increases of 14% and 12% in the average annual power generation than that of particle swarm optimization and genetic algorithm, respectively.
A case study by Mousavi et al. [21] using the crow algorithm for extracting rule curves to manage irrigation on the basis of the potential of crows for finding and hiding foods showed significant decreases in the vulnerability index. e average annual deficit in irrigation obtained by applying the krill algorithm to irrigation management and reservoir operation by Karami et al. [22] using the initial krill's position as the decision variable, respectively, decreased by of 20% and 23% compared to the deficit of the particle swarm and genetic algorithms. Applying the kidney algorithm based on the kidney function in the body to a reservoir operation to decrease power generation [23] showed that the kidney algorithm provided the best performance against the other evolutionary algorithms. It was stated that the KA achieved the first rank for the optimization problem and, thus, it showed a high potential to be applied for different problems in the field of water resources management. e authors of [24] developed the Multiobjective Grey Wolf Optimization (MOGWO) algorithm to obtain the optimum rules on the operation of the Golestan Dam in Golestan province, Iran, under the climate change conditions. eir results showed that the river flow would decline by 0.17 percent of the baseline period under climate change conditions in addition to increasing the temperature by 20% as well as decreasing the rainfall by 21.1%. In another research in the same year, they studied the Multiobjective Farmland Fertility Optimization (MOFFA) algorithm to derive optimum rules on the operation of the Golestan Dam under climate change conditions. e study results indicated an increase in release rates for climate change conditions in comparison with the baseline ones and stronger dam efficiency in climate change conditions [25].

Study Objectives, Novelty, and Problem Statement.
Great potential for solving water resources management problems has been reported for evolutionary algorithms.
Recently developed hybrid models based on 2 rather than a single metaheuristic have been applied to reservoir and dam water systems [24]. e hybrid particle swarm optimization and artificial fish optimization was used for a single-objective single-reservoir system problem. Establishing proper communication between the algorithms constitutes a major challenge in developing the hybrid models [24]. An easily implemented hybrid bat-particle swarm algorithm developed using a certain communication process was inefficient in terms of convergence rate [24]. Moreover, the shortcomings of evolutionary algorithms include being trapped in local optimal points, early convergence, producing immature solutions, and a low convergence rate [26]. Inspired from the social behaviors of humpback whales, Lewis and Mirjalili developed the whale optimization in 2016 as a swarm intelligence algorithm [27]. Inefficiency in examining the entire search space constitutes a major disadvantage of this algorithm. Given the acceptable potential of the differential evolution for exploration [28], the present study developed and implemented the improved whale optimization algorithm as a hybrid whale optimization-differential evolution algorithm in MATLAB and analyzed its application to managing water resources.

Differential Evolution.
Despite basic similarities between the differential evolution introduced by [29] and the structure of evolutionary algorithms, their difference lies in the way mutation and crossover are applied in the former to update solutions based on a standard evolutionary algorithm. In fact, the solution vectors randomly selected from the population undergo component swapping, addition, and subtraction before updating the generation. e differential evolution performs mutation as a constant and a weighting factor for differential vectors before crossover and mutation. Differential vectors are calculated as the average of several differential vectors obtained from different randomly selected pairs or the difference between 2 of randomly selected solution vectors from the population. An updated solution vector known as the noisy vector is obtained by adding the weighted difference vector to a solution vector that is selected randomly. e trial vector then obtained through a crossover between the target vector in the population and the noisy vector is ultimately compared to the target vector to pass the better one to the next generation. is process continues to fill up the next generation by the number of newly generated candidate solutions. e stages of the differential evolution are shown in Figure 1.
e randomly selected vectors need to be discriminated from one another. In contrast to genetic algorithm, no selection mechanisms are used in the differential evolution and minimum-cost vectors in minimization problems associated with trial vectors or parent target vectors are advanced to the next generation. e differential evolution is based on P (G) as the population of generation G including n pop individuals as candidate solutions. Position matrix P (G) is expressed as follows: where G max represents the maximum number of generations defined by users, n var the number of variables or diameters to be selected, X G i,j the jth component, i.e., that is, the jth link diameter, of the ith candidate solution, and X G i the ith candidate solution or individual in the Gth generation.
An arbitrarily generated initial population represented as P (0) in the differential evolution is as follows: in which X (L) j and X (U) j , respectively, represent the lower and upper bounds of variable jand randij a random variable uniformly distributed on [0, 1]. Subsequent generations P (G+1) or new vectors are generated successively through mutation as a combination of vectors that are selected randomly from the current population. Afterwards, a predetermined target vector is combined with the noisy vector. Moreover, P (G+1) as the trial vectors population is generated through mutation and recombination as follows: ]. F represents a mutation constant or weighting factor defined by users on [0, 1] and C r the minimization problem defined by users.

Shock and Vibration 3
Otherwise, the target vector is selected for the following generation, whose population P (G+1) is selected using the following equation: in which f(X (G) i ) denotes the cost of the i th individual of the G th generation [30]. e computational modules applied to optimally designing a water distribution network using the differential evolution are shown in Figure 2. Given the proper explorative potential of mutation in the differential evolution, this algorithm can effectively identify the zone of a global optimum despite its inappropriate exploitation [28]. No

Whale Optimization
Otherwise Figure 1: Flowchart of the differential evolution.
Mirjalil proposed the whale optimization as metaheuristic describing unique hunting behaviors of humpback whales. According to this algorithm, whales follow typical bubbles and create a circular or 9-shaped route while encircling prey during their hunting. e bubble-net feeding/hunting behavior can be explained by the fact that humpback whale go down in water by approximately 10-15 meters followed by producing spiral bubbles encircling the prey while following the bubbles and moving toward the surface. e following expression is  Computational modules applied to optimally design a water distribution network using the differential evolution.
Shock and Vibration 5 used to mathematically present the whale optimization algorithm.

Equation for
Encircling a Prey. After encircling small fish as the prey by humpback whales, they update their position towards the optimal solution as the number of iterations is maximized: in which X(t) represents the vector of position,X → * (t) the vector of position of the most recent optimal solution, t the current iteration, and A → and D coefficient vectors which can be calculated as follows: where r represents a random value on [0, 1] and a → a linearly decreasing variable from 2 to 0 throughout the iterations.

Bubble-Net
Attacks. Two methods proposed as the mathematical representation of the bubble-net behavior of humpback whales, as illustrated in Figure 3, are as follows.

Shrinking Encirclement.
is method is based on a linear decrease from two to zero in a → . e vector random value lies on [−1, 1].

Spiral Position Updates.
e spiral position between humpback whales and the prey following a helix shape can be mathematically expressed as follows: in which b is a constant defining the logarithmic shape, l a random value over [−1, 1], and D → � |X * � �→ (t) − X(t)| the distance between the prey and the ith whale used to update the solution. e probability of the whale following a logarithmic path or perform shrinking encircling during optimization is assumed to be 50%. e related mathematic equation is expressed in the following: with p being a random value over [0, 1].

Searching for the Prey.
A → is used to find the prey on the basis of the following 2 conditions while taking values above 1 or below −1: e following conditions are ultimately met: (i) | A → | > 1 ensures that traps in local optimal solutions are avoided throughout the search by the whale optimization algorithm. (ii) | A → | < 1 is selected to update the current search agent/ best solution in terms of position [27]. e whale optimization algorithm flowchart is illustrated in Figure 4.

Improved Whale Optimization.
is section presents improvement in the whale optimization algorithm. e present study integrated the mutation of the differential evolution with the whale optimization algorithm to improve and balance the exploration and exploitation phases and increase its convergence rate. e algorithm proposed is explained as follows.
e trade-off existing between the exploration and exploitation of the whale optimization algorithm was solved in the improved whale optimization algorithm using a hybrid model comprising the operators of this algorithm and the mutation operator of the differential evolution. e main hybrid operator of the improved algorithm (see lines 10-28 of proposed algorithm) is a combination of mutation in the differential evolution and the components of the whale optimization algorithm, namely, encircling prey, search for prey, and spiral updating position. e improved whale optimization algorithm consists of two main parts, that is, exploration (see lines 11-18 of proposed Algorithm) and exploitation (see lines 19-27 of proposed Algorithm). If rand < k, the exploration part changes the individuals. Equation (8) was used to adjust k.

Bubbles
Fish colony Humpback whale 6 Shock and Vibration in which t max and t, respectively, show the maximum number of generations and the current generation and k is utilized to control the exploration and exploitation of the improved whale optimization algorithm. According to (8), decreasing k over time from 1 to 0 allows the individuals to explore the initial generation, while doing exploitation as time passes. e exploration part contains the mutation of the differential evolution and search for the prey of the whale optimization algorithm. e improved algorithm is integrated with mutation to enhance exploration throughout the search space [28]. e exploitation component of the improved whale optimization algorithm resembles that of the whale optimization, although the improved algorithm is an elitist method. In fact, the updated position of the ith individual in the next generation is selected as the fitter one of offspring U i and parent X i . e solution is obtained by considering the boundary conditions, which should not be violated; otherwise, the repairing rule is applied according to (9).
with rndreal (0, 1) being a random value between 0 and 1, X i (j) the jth dimension of the ith solution, while δ j and μ j , respectively, show the lower and upper limits of the jth dimension. Improved whale optimization algorithm is as follows.
(1) Generating an initial population X i , i � 1, 2, . . ., N P (2) Measuring the individuals fitness in X i (3) Calculating X * as the best individual (4) While the stopping criterion is not satisfied (5) Mapping the fitness to the number of species for every individual (6) For i � 1 to N P (7) Selecting a uniformly-distributed random value r1-r2-r3-i (8) Updating a � 2 − t(2/t Max ), A � 2ara, C � 2r, l � rndreal(−1, 1), p � rndreal(0, 1) and k � 1 − (t/t Max ) (9) j rand � rndint(1, n) (10) For j � 1 : n (11) If p_k then (12) If rndrealj (0, 1)_CR or j � � j rand then (13) (38) End while e structure of the proposed improved whale optimization algorithm is simple. e only difference between the improved algorithm and the whale optimization lies in integrating the mutation of the differential evolution with  Shock and Vibration 7 the latter algorithm. e selection function of the improved algorithm resembles that of the differential evolution. e improved algorithm therefore computationally resembles the whale optimization based on O(MaxGen × Np × O(fitness)). e maximum number of generations is denoted as maxgen and O(fitness) is calculated by the application [31]. Figure 5 shows the different stages of the improved whale optimization algorithm.

Case Study.
Boostan is a clay-core Earth dam in northeastern Gonbad-eKavus, Golestan Province, Iran, for satisfying an irrigation demand of 4200 ha of land and environmental and industrial water demand. e dam's length is 665 meters, its top width 10 meters, the minimum storage volume of its reservoir 0.1 × 10 6 m 3 , and the maximum storage volume 43.7 × 10 6 m 3 . According to Figure 6 showing the monthly inputs, the present study was conducted during 151 months from October 2005 to April 2018.
Boostan Reservoir Dam lies in the Gorganrood-Gharasoo Basin as part of the 13061 km 2 Caspian Basin in the southeast of the Caspian Sea. Climate change affects the average annual precipitation in the watershed from the east to west and from the north to south. e average annual temperature in this basin ranges from approximately 17°C at low altitudes to 7.5°C in the southern highland area. e geographical situation of the present research area is illustrated in Figure 7.
Findings obtained from the present case study are recommended to be used by policymakers to satisfy the irrigation demand in agricultural regions. e released water was considered a decision variable, and the water volume should be adequate for supplying the demand of the downstream farmers. e shortage of water was minimized using the following objective function: in which De max represents the maximum demand for water during the operation, Re t is the discharged water, De t is the required volume, and F Release is the objective function. e following continuity equation was also considered: in which Sp (t) represents overflow at time t, Re (t) the discharged water, Ev (t) evaporation from the reservoir surface, S (t) the reservoir storage volume at time t as well as an input to the reservoir, and S t+1 the reservoir storage volume at time (t + 1). e following constraints were applied to overflow from the reservoir: e other constraints were as follows: Re min ≤ Re t ≤ Re max , (16) in which Re min represents the minimum discharge volume of the reservoir, Re max the maximum discharge volume, S min the minimum volume of the reservoir, and S max its maximum volume. e following penalty functions were included in the objective function to take account of potentially unsatisfied constraints [32].
e indicators used to measure the performance of various algorithms in solving water resource management problems were as follows. e volumetric reliability indicator is defined as the volume of water discharged in the whole period divided by the overall demand for water as follows [33]: (1) e time-based reliability indicator showing the percentage of time the demand is completely satisfied by the system is expressed as follows [32].
in which NDef represents the total number of failures that occurs. Higher values of this index suggest a higher system performance. (2) e vulnerability index refers to the highest failure rate throughout a reservoir operation, with lower values suggesting a higher performance [34]: (3) e resiliency indicator reflects the recovery rate of a system after failures; for instance, the occurrence of failures in 4 months of a 12-month operation 8 Shock and Vibration suggests the significant effect of the sequence of failure periods on the system. Moreover, larger values of this indicator are more desirable.
in which F i represents the total frequency of failure periods and f si the frequency of occurrence of failure series [26]. (4) e sustainability index defined by performance criteria as follows is defined as the sum of system performance measures in one general index to facilitate comparison and decision-making between different management alternatives and water resources plans [35,36]: 3.7. Calculating the Weights of the Criteria. Subjective weights were mainly calculated using expert comments according to experiences and subjective judgments such as Delphi and AHP, whereas objective weights were obtained reservoir from actual alternative information associated with techniques, including the entropy weight, which decreases the subjective effect of decision makers and improves objectivity. As a concept of a thermodynamic origin, entropy was used to calculate system disorders and degree of confusion. Entropy was also used by Shannon to solve information theory problems and resolve uncertainty [37]. e entropy theory suggests lower values of the entropy result in more available data. A higher weight can be therefore assigned to the criterion [38]. Numerous studies have used entropy weight in different fields. e Shannon entropy weight can be calculated through the following steps [39] by assuming a decision matrix with n criteria and m alternatives. e first stage is normalizing the decision matrix: Crossover and Selection operation using equation (3 and 4) Obtain the individual of next generation Search for the prey Update the position of the search agent using equation (11) and (12) Update the position of the search agent using equation (5) and (6) Encircling the prey Update the best solution X * Evaluate the parameters r, a, A, C, and P

Bubble-net attacking
Update the position of the search agent using equation (9)    Shock and Vibration 9 e second stage is calculating entropy: in which K � 1/In m. e third stage is calculating the weight of individual criteria: 3.8. TOPSIS. Yoon and Hwang proposed TOPSIS as a multicriteria decision-making approach [40] for selecting an optimal alternative. According to the main rule of TOPSIS, the best alternative is, respectively, nearest and farthest from positive ideal and negative ideal solutions [41]. Multicriteria decision-making problems have been commonly solved in numerous fields using this technique. According to Yoon and Huang, the following algorithm is presented for TOPSIS [42]. e first stage is building R as the normalized decision matrix: e second stage is building the weighted normalized decision matrix: in which w j represents the j th criterion weight. e third stage is obtaining the positive ideal and negative ideal solutions, respectively, shown as A + and A − as follows: in which J and J′, respectively, represents the sets of benefit and cost criteria. e fourth stage is calculating the distance between individual alternatives and positive ideal and negative ideal solutions: e fifth stage is calculating the closeness coefficient and prioritizing the alternatives: in which cl * i ∈ [0, 1], i � 1, 2, . . . , m. An optimal alternative is selected based on the priority of cl * i whose higher values denote a better alternative. e closer the cl * i to unity, the closer the alternative A i to the positive ideal solution.

Discussion
Tables 1-3 present the results of a sensitivity analysis on the random variables of various algorithms. e optimal population size was obtained as 60 for the improved whale optimization algorithm with an objective function value of 25.74 being the lowest compared to the other population sizes. e optimal mutation factor f was obtained as 0.4 and the optimal cross over rate c as 0.3 for the differential evolution and its objective function value of 25.84 was found for these two values.
According to Table 4, after randomly running the individual algorithms for ten times, their performance, convergence rate, and statistical parameters were compared with one another and with the absolute optimal values obtained of GAMS as a nonlinear optimization technique. Ten times of running the improved whale algorithm, the differential evolution and the whale optimization algorithm, respectively, resulted in average responses of 25.74, 25.84, and 32.91. e objective function based on the improved whale optimization algorithm therefore minimized water supply shortages. e calculation duration was obtained as 521 seconds for the improved whale algorithm, 527 seconds for the whale optimization algorithm, and 1040 seconds for the differential evolution, suggesting that the improved algorithm yields the minimum computational duration compared to those of the other two algorithms. e lower variation coefficient obtained after running the improved whale algorithm for 10 times compared to that of the other algorithms suggested the adequately reliable results of even one execution of the improved whale algorithm. Although the globally optimal solution obtained with GAMS was 12.553, the hybrid whale optimization outperformed the whale optimization algorithm given that the former improved the performance of the latter by up to 21.7% in deriving the optimal values of the water supply function through decreasing the variation coefficient by up to 78% and the standard deviation by up to 82%.
Comparing the three algorithms in terms of convergence rate in Figure 8 showed that the improved whale algorithm convergence is the fastest and therefore outperforms in terms of computation duration and response quality. Replacing a poor solution obtained from an algorithm with a good solution obtained from another algorithm improved the rate of convergence of the improved whale algorithm. In fact, despite the poor individual performance of the differential evolution and the whale optimization algorithm, their hybrid performance was higher.
Given the integration of the mutation operator of the differential evolution with the whale optimization, Figure 8 confirmed the performance of the improved whale algorithm by demonstrating the lower minimum values obtained using this algorithm compared to those obtained from the whale optimization algorithm. e three algorithms were compared in Table 5 in terms of satisfying the demand for water in the 151-month study period. e correlation coefficient between the required and actual water discharge obtained as 0.94 for the improved whale algorithm was higher than that of the other algorithms.
e root-mean-square errors obtained as 0.72 × 10 6 m 3 for the improved algorithm, 2.64 × 10 6 m 3 for    Shock and Vibration the differential evolution, and 1.99 × 10 6 m 3 for the whale optimization algorithm suggested the higher performance of the improved whale optimization algorithm in satisfying demand for water. e lower mean absolute error presented in Table 5 for the improved algorithm compared to that of the other evolutionary algorithms showed that this algorithm outperforms in terms of meeting the downstream demand for water.
With an average demand for water of 3.21 × 10 6 m 3 , Figure 9 shows the improved whale optimization algorithm resulted in a water discharge volume of 3.03 × 10 6 m 3 from Boostan dam for supplying the demand for the reservoir water over the study period, whereas the differential evolution yielded 1.18 × 10 6 m 3 and the whale optimization 1.99 × 10 6 m 3 . e improved algorithm was therefore the best algorithm in terms of meeting the demand for water. Table 6 compares the performance of the study algorithms by the water resources management indicators, suggesting that the improved whale optimization algorithm, with a reliability indicator of 94%, outperformed the other evolutionary algorithms in terms of satisfying the downstream demand for water and with the largest resiliency indicator (89%) and the minimum objective function (25.74) has the best status. e optimal decision can therefore be made through considering several indicators and using multicriteria decision-making models. In addition to objective function, the other indicators considered in the present case study to select the best algorithm for supplying the downstream water demand included reliability, vulnerability, and resiliency. Comparing the whale optimization algorithm and the improved whale algorithm confirms performance improvements in terms of all the evaluative indicators. Table 7 presents the positive and negative criteria and the criteria weights obtained using the entropy method. e results showed that two criteria, that is, resiliency and time reliability, have, respectively, the highest weights  1  29  57  85  113  141  169  197  225  253  281  309  337  365  393  421  449  477  505  533  561  589  617  701  729  757  785  813  841  869  925  953  981  645 673 897 DE WOA Figure 8: Comparing the three algorithms in terms of convergence rate.  1  5  9  13  17  21  25  29  33  37  41  45  49  53  57  61  65  69  73  77  81  85  89  93  97  101  105  109  113  117  121  125  129  133  137  141  145  149 Month Demand RE-WOA Re-IWOA Figure 9: Volume of water supplied using the different algorithms. followed by the other criteria prioritized as the third to the sixth. e results obtained from TOPSIS on the basis of the decision matrix are shown in Table 8, suggesting that the improved whale optimization algorithm outperformed the other algorithms with the maximum similarity index, whereas the whale optimization and the differential evolution were, respectively, ranked second to third.

Conclusion
e present research proposed the improved whale optimization as a novel hybrid whale optimization differential evolution algorithm to optimize the operation of Boostan dam reservoir through satisfying the downstream demand for water by Boostan dam reservoir in Golestan province and reducing the shortage of water using a single-reservoir system. Exploration was performed in the improved algorithm by combining the mutation operator of the differential evolution and the search for prey of the whale optimization algorithm. e convergence rate, performance, and statistical parameters of the three algorithms were compared with one another and with the globally optimal solution obtained from nonlinear programming in GAMS by running the individual algorithms for ten times. e obtained results confirmed the performance superiority of the improved whale optimization algorithm over the differential evolution with the objective function value of 25.74, minimum computation duration, and fastest convergence. In addition, the improved whale optimization improved the objective function, standard deviation, and variation coefficient by 21.7%, 82%, and 78% compared to those of the whale optimization, respectively. Lower values of mean absolute error and root-mean-square error obtained for the improved algorithm also suggest this algorithm can more effectively help supply the demand for water than the other evolutionary algorithms. Given an average demand for water of 3.21 × 10 6 m 3 over the 151-month study period, the improved whale optimization algorithm outperformed the other evolutionary algorithm with a water discharge volume of 3.03 × 10 6 m 3 , which was higher than that of the other algorithms. In other words, the improvements in performance achieved can be explained by the integration of the mutation operator of the differential evolution with the whale optimization in the hybrid algorithm, which caused the mutation operator to search the exploration space more effectively.
A multicriteria decision-making model was utilized based on TOPSIS and the Shannon entropy to identify the optimal algorithm. Calculating the evaluative indicators for the three algorithms suggested that the improved whale optimization algorithm, with a reliability index of 94%, outperforms the other evolutionary algorithms in terms of meeting the downstream demand for water. is algorithm was also found to be the best in terms of objective function value (25.74) and resiliency index (89%). e multicriteria decision-making model employed in the present research confirmed the functional superiority of the improved algorithm over the other two algorithms in optimizing the operation of Boostan dam reservoir. Moreover, an increase was observed in the volume of Boostan dam reservoir during the exploitation using the improved whale optimization compared to the whale optimization. In line with the results obtained by Kai et al. [4,32], the present findings showed the superiority of hybrid algorithms over their constituent algorithms. e effect of climate change on this novel hybrid algorithm is recommended to be examined for multireservoir operation.