An Ameliorative Whale Optimization Algorithm for Multi-Objective Optimal Allocation of Water Resources in Handan, China

: With the deepening discrepancy between water supply and demand caused by water shortages, alleviating water shortages by optimizing water resource allocation has received extensive attention. How to allocate water resources optimally, rapidly, and effectively has become a challenging problem. Thus, this study employs a meta-heuristic swarm-based algorithm, the whale optimization algorithm (WOA). To overcome drawbacks like relatively low convergence precision and convergence rates, when applying the WOA algorithm to complex optimization problems, logistic mapping is used to initialize swarm location, and inertia weighting is employed to improve the algorithm. The resulting ameliorative whale optimization algorithm (AWOA) shows substantially enhanced convergence rates and precision than the WOA and particle swarm optimization algorithms, demonstrating relatively high reliability and applicability. A water resource allocation optimization model with optimal economic efﬁciency and least total water shortage volume is established for Handan, China, and solved by the AWOA. The allocation results better reﬂect actual water usage in Handan. In 2030, the p = 50% total water shortage is forecast as 404.34 × 10 6 m 3 or 14.8%. The shortage is mainly in the primary agricultural sector. The allocation results provide a reference for regional water resources management.


Introduction
Water is one of the most indispensable natural resources for human survival and development [1].Rapid social and economic development and population growth, along with climate change, have aggravated water shortages [2][3][4][5][6][7][8].The recent situation is unlike early water supply management.There are fewer exploitable water resources and exploitation costs are increasingly expensive, often leading to more serious environmental problems.By optimizing water resource allocation, water resources are transferred from low-to high-value uses, and the efficiency of water usage is enhanced.These are important ways to solve water shortages.
Optimal allocation of water resources has been studied for years all over the world.The first research focused on single-objective water allocation management in an irrigation area.However, it has grown to cover areas like mega-scale multi-objective regional water resource allocation optimization.The coverage and scope of water resource allocation have expanded greatly.The earliest research on water resource allocation optimization dates to Masse et al. '

s study on reservoir optimization
Water 2018, 10, 87 3 of 29 low convergence precision and convergence rates.In order to overcome these issues, this study proposes the ameliorative whale optimization algorithm (AWOA), and compares its performance to WOA and PSO's performance in solving multi-objective water resource allocation optimization models.The results indicate that the AWOA yields higher convergence precision and rates.

Whale Optimization Algorithm
The WOA is a novel swarm intelligence optimization algorithm proposed by Australian researchers Mirjalili and Lewis in 2016 [19].This algorithm is inspired by the hunting mechanism of humpback whales in nature, and simulates the shrinking encircling, spiral updating position, and random hunting mechanisms of humpback whale pods.This model consists of the following three stages: Encircling prey, bubbling-net attacking, and search for prey [19].

Encircling Prey
For encircling prey, the WOA optimization algorithm assumes that the current best solution is the target prey.Other individual whales then try to update their positions toward the optimal position.This behavior is represented by the following equations: a = 2 − 2 t t max (5) where t represents the current iteration, → X(t) is the position vector at iteration t, → X p (t) is the position vector of the best solution obtained so far, D represents the distance between the position of → X p (t) and → X(t),

→
A and → B are coefficient vectors, ⊗ is an element-by-element multiplication operator (e.g., (x 1 , 1], a is a control parameter that is linearly decreased from 2 to 0 over the course of iterations, and t max is the maximum number of iterations.

Bubble-Net Attacking Method (Exploitation Phase)
The bubble-net behavior of humpback whales includes the shrinking, encircling, and spiral updating position.

• Shrinking encircling mechanism
This mechanism is mainly achieved by decreasing the value of control parameter a.

→
A is a random value in the interval [−a, a] (a is linearly decreased from 2 to 0).When random values for → A are in the interval [−1, 1], the new position of the individual whales can be defined anywhere between the original position and the current best position.The mathematical modeling is expressed by Equations ( 1) and (2).

• Spiral updating position
Water 2018, 10, 87 4 of 29 For the spiral updating position, a logarithmic spiral equation is created between the position of the whales and the current best solution to simulate the spiral movement of humpback whales.The mathematical model is shown as → X(t + 1) = D e bl cos(2πl) + → X p (t) where D represents the distance between the i-th whale and its prey (current best solution), b is a constant for defining the shape of the logarithmic spiral, and l is a random number in [−1, 1].
When humpback whales attack their prey, they move simultaneously within a shrinking encircling circle and along a spiral-shaped path.The WOA assumes there is a 50% probability of choosing between the shrinking encircling mechanism and the spiral model to update the position of whales during optimization.The mathematical model is as follows: → X(t + 1) = D e bl cos(2πl) + → X p (t) p ≥ 0.5 (9) where p is a random probability in [0, 1].

Search for Prey (Exploration Phase)
In the search for prey or exploration phase, the positions of other individual whales are updated according to a randomly chosen individual whale.By setting || ≥ 1, individual whales are forced to move far away from the reference whale in search of other better prey, thereby allowing the WOA algorithm to perform a global search.The model can be expressed as where → X rand is a position vector randomly selected from the current whale population.

Ameliorative Whale Optimization Algorithm
Compared to other optimization algorithms, the advantages of the WOA are its simple operation and limited number of adjustment parameters, as there are only two adjustment parameters ( Through the adjustment of → A, the algorithm can achieve a good balance between the capabilities of exploration and exploitation, thereby increasing the possibility of getting away from the local optimum.However, the exploitation phase in the WOA is completely dependent on randomness, thereby resulting in low convergence accuracy and slow convergence speed.
In response to the aforementioned drawbacks of the WOA, the AWOA is proposed in this study.The improved method is described in detail below.
First, in order to effectively improve the quality of the initial population position, the logistic map [37] is used to obtain the initial solution.
where X k ∈ [0, 1] and µ is a control parameter, µ = 4. Second, after analyzing the weaknesses of the PSO algorithms, Shi and Eberhart [38] introduced the inertia weight ω to facilitate rapid convergence to a global optimal solution.Their analysis indicated that a larger inertia weight facilitates a global search, while a smaller one facilitates a local search.In order to improve local search capability, increase convergence accuracy, and speed up convergence, inertia weight is introduced in the WOA to improve the algorithm.The inertia weight introduced is expressed as follows: ω is the maximum value of inertia weight, ω is the minimum value of inertia weight, t is the current iteration number, and T max is the maximum number of iterations.The improved position vector update formula is shown as ω is progressively decreased with the number of iterations, thereby making it useful for global search and local optimization at the start and end of the iteration process, respectively.As the rate of decrement in the introduced ω is large, it facilitates local optimization in the algorithm, increases convergence accuracy, and speeds up convergence.

Best Individual Selection Based on Pareto
In a single-objective problem, the fitness function value is the objective function value.In a multi-objective problem, based on the solution strategy of the Pareto solution set, all individuals in the population can be sorted into different non-inferior frontiers using a fast non-dominated sorting algorithm.With a crowded-comparison operator, the best individual can be selected by setting the crowding distance as the winning standard among the same level groups after sorting.

Fast Non-Dominated Sorting
With fast non-dominated sorting, relatively good individuals in the population can be selected, thereby reflecting the superiority or inferiority of the individuals.The population size Q is assumed to be N; i is an individual within Q; and two parameters of the individual i, namely, n i (the number of dominating individuals in the population) and S i (the set of individuals dominated by individual i in the population), have to be calculated.First, all individuals in the population with n i = 0 are sorted into F 1 , and F 1 is taken as the current set.For each individual k in set F 1 , those dominated by individual k are sorted into S k .For each individual in the set S k , 1 is subtracted from each n k value, if n k = n k − 1 = 0, and then, the corresponding individual is saved into set H. Finally, all individuals in set F 1 are recorded as the first level of non-dominating frontiers, while set H is taken as the current set.Then, the steps mentioned above are repeated until the entire population has been sorted [39][40][41][42].

• Crowding distance
To improve the diversity and distribution of a population, the concept of crowding distance was proposed by Srinivas and Deb [43].Crowding distance indicates the density surrounding an individual n in the population.The available space around individual n is represented by the length of the largest rectangle-containing individual n, which is expressed as n d .As shown in Figure 1, the formula for calculating the crowding distance of individual i [39] is Water 2018, 10, 87 where f m (i) represents the objective function value of individual i. M is the number of objective functions.The objective function corresponding to individual i needs to be normalized because of its different units.According to [44], Equation ( 16) can be rewritten as where f mmax represents the maximum value obtained by all individuals in the population in the i-th objective function, and f mmin represents the minimum value obtained by all individuals in the population in the i-th objective function.For each objective function, after all individuals in the population have undergone fast non-dominated sorting, the crowding distance of two individuals at the boundary of each level is infinite.[44], Equation ( 16) can be rewritten as where ) and (18) When the non-dominated ranks between two individuals are different, the individual with the lower rank is considered superior; when the non-dominated ranks between two individuals are the same, the individual with the greater crowding distance is considered superior [39].

AWOA Solving Steps
1. Initialization of parameters: Whale population size: N, with random generation of the initial population; maximum number of iterations: Tmax; maximum value of inertia weight:  ; minimum value of inertia weight:  .
2. The location of whale population is initialized using logistic mapping, and initial variables are transformed into a range of values corresponding to initial variables of the objective function.3. Fast non-dominated sorting and calculation of individual crowding degree are carried out for all individuals in the population, and the best individual is then selected as the optimal position based on Equation (18).4. When  • Crowded-comparison operator After carrying out the fast non-dominated sorting operations and calculation of crowding degree for each individual i in the population, each individual acquires two attributes, i rank (non-dominated rank) and i d (crowding degree).i rank and i d can be used to differentiate between the dominance and non-dominance relationship of any two individuals in the population.The crowded-comparison operator is expressed as ≥ n : i ≥ n j(i f (i rank < j rank ) or (i rank = j rank and i d > j d ))

p
When the non-dominated ranks between two individuals are different, the individual with the lower rank is considered superior; when the non-dominated ranks between two individuals are the same, the individual with the greater crowding distance is considered superior [39].

1.
Initialization of parameters: Whale population size: N, with random generation of the initial population; maximum number of iterations: T max ; maximum value of inertia weight: ω ; minimum value of inertia weight: ω .

2.
The location of whale population is initialized using logistic mapping, and initial variables are transformed into a range of values corresponding to initial variables of the objective function.

3.
Fast non-dominated sorting and calculation of individual crowding degree are carried out for all individuals in the population, and the best individual is then selected as the optimal position based on Equation (18).4.
When p < 0.5 and |A| < 1, the current position of individual whales are updated with Equation ( 14); when |A| ≥ 1, a random individual whale → X rand is selected, and the current positions of individual whales are updated with Equations ( 10) and (11).
When p ≥ 0.5, the current positions of individual whales are updated with Equation (15).6.
All individual whales are checked to determine if their updated locations are beyond the search space, and whale locations beyond the search space are corrected.7.
The number of iterations is checked.If the maximum number of iterations has been reached, the algorithm is stopped; otherwise, proceed to step (3) to continue the algorithm.

Background of Study Area
Handan is located in the south of Hebei Province (latitude 36 • 04 -37 • 01 N, longitude 113 • 28 -115 • 28 E), and covers an area of 12,047 km 2 .Taking 100 m contours as the boundary, the western part is covered by middle to low mountains, hills, and basins over an area of 4460 km 2 , while the eastern part consists of plains, over an area of 7587 km 2 .Handan has a temperate semi-humid and -arid continental monsoon climate with an average annual precipitation of 548.9 mm concentrated in June to September.Handan is divided into 16 administrative divisions (Figure 2)-Handan city, Shexian, Wu'an, Fengfeng, Cixian, Yongnian, Quzhou, Jize, Linzhang, Cheng'an, Weixian, Guangping, Feixiang, Daming, Guantao, and Qiuxian.
Water 2018, 10, 87 7 of 27 6.All individual whales are checked to determine if their updated locations are beyond the search space, and whale locations beyond the search space are corrected.7. The number of iterations is checked.If the maximum number of iterations has been reached, the algorithm is stopped; otherwise, proceed to step (3) to continue the algorithm.

Background of Study Area
Handan is located in the south of Hebei Province (latitude 36°04′-37°01′N, longitude 113°28′-115°28′E), and covers an area of 12,047 km 2 .Taking 100 m contours as the boundary, the western part is covered by middle to low mountains, hills, and basins over an area of 4460 km 2 , while the eastern part consists of plains, over an area of 7587 km 2 .Handan has a temperate semi-humid and -arid continental monsoon climate with an average annual precipitation of 548.9 mm concentrated in June to September.Handan is divided into 16 administrative divisions (Figure 2)-Handan city, Shexian, Wu'an, Fengfeng, Cixian, Yongnian, Quzhou, Jize, Linzhang, Cheng'an, Weixian, Guangping, Feixiang, Daming, Guantao, and Qiuxian.Handan is a region with water shortages because of insufficient water resources.Average water resources per capita are less than 180 m 3 , equivalent to only one-twelfth of the national average and far lower than the global average.With the rapid socioeconomic developments in the region, water demand has surged.Reduced precipitation and uneven spatial and temporal distribution of precipitation due to climate change have intensified regional water shortages, which can be alleviated only by over-exploitation of groundwater.Groundwater over-exploitation causes associated levels to drop continuously, and cones of depression to expand.Other severe consequences are scrapped wells, ground subsidence, and earth fissures.To safeguard economic development, apart from the South-North Water Transfer Project, Handan has imported large amounts of water from the Yellow River to alleviate regional water shortages problem.In order to ensure water supply for the public, industry, and ecology, as well as attain maximum socioeconomic benefits, it is necessary to examine how imported water can be rationally utilized and coordinated with the local water supply.Handan is a region with water shortages because of insufficient water resources.Average water resources per capita are less than 180 m 3 , equivalent to only one-twelfth of the national average and far lower than the global average.With the rapid socioeconomic developments in the region, water demand has surged.Reduced precipitation and uneven spatial and temporal distribution of precipitation due to climate change have intensified regional water shortages, which can be alleviated only by over-exploitation of groundwater.Groundwater over-exploitation causes associated levels to drop continuously, and cones of depression to expand.Other severe consequences are scrapped wells, ground subsidence, and earth fissures.To safeguard economic development, apart from the South-North Water Transfer Project, Handan has imported large amounts of water from the Yellow River to alleviate regional water shortages problem.In order to ensure water supply for the public, industry, and ecology, as well as attain maximum socioeconomic benefits, it is necessary to examine how imported water can be rationally utilized and coordinated with the local water supply.
Sources for the water supply in Handan include local surface water, groundwater, recycled water, the Yuecheng reservoir, the Dongwushi reservoir, water transferred through the South-North Water Transfer Project, and water imported from the Yellow River.Medium-and small-scale reservoirs, impoundment, and ponds in the administrative divisions support the local surface water supply.Groundwater refers to the exploitable groundwater in the divisions.Recycled water is the water supply that sewage treatment plants in all divisions can provide.The Yuecheng reservoir and the Dongwushi reservoir are two essential water supply sources in the region.They are responsible mainly for the agricultural, industrial, and ecological water supply in the eastern plain area.In particular, the Yuecheng reservoir supplies water to eight counties, including Cixian, Linzhang, Cheng'an, Weixian, Guangping, Feixiang, Daming, and Guantao, while the Dongwushi reservoir supplies water to Handan city and five counties, namely, Yongnian, Quzhou, Jize, Feixiang, and Qiuxian.Water transferred through the South-North Water Transfer Project is relatively expensive because it is transported over a long distance.According to the project requirements, this water can be used only as domestic and industrial water.It is supplied to Handan city and 12 counties-Cixian, Yongnian, Quzhou, Jize, Linzhang, Cheng'an, Weixian, Guangping, Feixiang, Daming, Guantao, and Qiuxian.Water imported from the Yellow River is sourced through the Dongfeng drainage in the eastern plains, with an average volume of 100 × 10 6 m 3 .Due to the water quality, this water source is mainly used as agricultural water and supplements industrial and ecological water.It is supplied to seven counties: Quzhou, Weixian, Guangping, Feixiang, Daming, Guantao, and Qiuxian.

Optimal Water Resource Allocation Model
According to the characteristics of the Handan water resources system, the region is divided into 16 sub-regions based on administrative divisions.Socioeconomic benefits are selected as objective functions, and are combined with corresponding constraint conditions to establish a regional multi-objective water resource allocation optimization model.Water sources in the region can be classified into public water and independent sources according to the scope of their supplies.Public water sources can supply water to two or more sub-regions simultaneously, like the Yuecheng reservoir, the Dongwushi reservoir, water transferred via the South-North Water Transfer Project, and water imported from the Yellow River.Independent water sources are only capable of supplying water in their own sub-regions (e.g., local surface water, groundwater, and recycled water).Water uses in Handan include those of urban and rural domestic use; primary, secondary, and tertiary sector use; and ecological uses.The primary sector covers agriculture, forestry, livestock farming, and fisheries; the secondary sector comprises industries that produce finished and usable products and construction industries; the tertiary sector includes catering and services.

Objective Functions
• Social benefits target Since it is difficult to measure social benefits, the least regional water shortage volume is used to represent social benefits indirectly in this study.
where D kj is the water demand of user j in sub-region k (×10 6 m 3 ).

• Economic benefits target
Water 2018, 10, 87 Economic benefits are expressed by the maximum sum of net efficiency yielded through the supplies by regional independent and public water sources.
where x kij and x kcj denote the water supplies by independent source i and public source c to user j in sub-region k (×10 6 m 3 ), respectively; b kij and b kcj are the coefficients of efficiency per unit water supply by independent source i and public source c to user j in sub-region k, respectively, (CNY/m 3 ); c kij and c kcj are the coefficients of cost per unit water supply by independent source i and public source c to user j in sub-region k, respectively (CNY/m 3 ); α ki and α kc denote coefficients of water supply sequence by independent source i and public source c in sub-region k, respectively; and β kj is the coefficient of water supply fairness of user j in sub-region k.

Constraint Conditions
• Water supply capacity constraints by the water supply system Independent water source Public water source where W ki and W c are the available water supplies by independent source i and public source c in sub-region k, respectively.
• Water delivery capacity constraints by the water delivery system Independent water source Public water source W kc ≤ Q kc (24) where Q ki and Q kc denote the maximum water delivery capacities of independent source i and public source c in sub-region k, respectively.
• Users' water demand constraints where D kjmin and D kjmax are the minimum and maximum water demands of user j in sub-region k, respectively.

Water Demand and Supply Forecasting
Taking the situations in 2015 as the status quo, according to the long-term development planning of Handan, the water demand of six user types in the region in 2030 under a guaranteed rate of p = 50% is predicted.The results are presented in Table 1.The prediction was carried out using the quota method.The quotas were determined with reference to Hebei Province Water Use Quota DB13/T 1161-2016 [44].The upper limit of the water demand of each user in different sub-regions in the target year is taken as the predicated water demand.The lower limit is 80% of the predicted volume.The available water supply in Handan in 2030 is estimated using the prediction on three types of independent water sources and four types of public water sources.For independent sources, the available water supply forecast on local surface water and groundwater is based on the results of the Second Water Resources Evaluation in Handan [45].The recycled water is calculated by multiplying the urban domestic and industrial water use volumes by their corresponding sewage conversion factor, followed by the sewage collection factor and then the recycled water conversion factor.For public sources, the available supply by the South-North Water Transfer Project directly uses the results in Supporting Works Plan in Handan for the South-North Water Transfer Project (Middle Line) in Hebei Province [46].The available supply from water imported from the Yellow River uses the available allocated volume listed in the Proposal for the Yellow River-to-Baiyangdian Water Transfer Project [47].The available supplies by the Yuecheng reservoir and Dongwushi reservoir are based on research by Wang and Li [48].The amount of available water supply in Handan is shown in Table 2.The efficiency coefficients for domestic and environmental water uses are, in general, difficult to determine.Since domestic water use is given the highest priority, the efficiency coefficients for urban and rural domestic water use are taken as 600 and 500 CNY/m 3 , respectively.Those for primary sector water use are defined by multiplying the yield increase efficiency by the water allocation coefficient and are assumed to be 10 CNY/m 3 in this study.Those for secondary and tertiary sector water use are based on the reciprocals of the water quotas and are 450 and 400 CNY/m 3 , respectively.To maintain sustainable development, a relatively high value is assigned for the efficiency coefficient for ecological water use; in this study, it is 350 CNY/m 3 [49,50].

Coefficient of Water Supply Cost c kij
According to the latest water cost in Handan, the water supply coefficient is defined.The coefficients for urban and rural domestic water use are 3.75 and 1.2 CNY/m 3 , respectively.The coefficients for primary, secondary, and tertiary sector water use are 0.38, 9.2, and 5.8 CNY/m 3 , respectively, and 1.8 CNY/m 3 for ecological water use.

Coefficient of Water Supply Sequence α ki
The water supply sequence coefficient α ki reflects the water supply priority of water source i in sub-region k, relative to other sources.The priority of each source is converted into a coefficient in the [0, 1] interval, that is, the coefficient of water supply sequence.n ki denotes the water supply order of source i in sub-region k and n kmax represents the maximum order in sub-region k.The determination of α ki is based on the following formula [51]: There are seven types of water supply sources in the Handan water resource allocation optimization model.Different users in each sub-region have different needs for the water supply sequence of these seven sources.First, the supply sequences of these seven sources for different users are determined.Second, the sequence coefficients are calculated using Equation (28).The coefficients for different water users are shown in Table 3.The water supply fairness coefficient β kj measures the priority of user j receiving water supply in sub-region k relative to other users.β kj and α ik are similar; both are associated with the water supply sequence of users.First, according to the nature and importance of a user, the corresponding water supply sequence is obtained.Then, the coefficient is determined through the formula for α ik .

Results Analysis
The AWOA is employed to solve the Handan water resource allocation optimization model.The initialization swarm size is 150 and there are 180 iterations.There are 20 runs, and each run starts with a different random initialization swarm.Through numerous simulation experiments, the maximum and minimum inertia weights of the AWOA are finally determined to be ω = 0.9 and ω = 0.2, respectively.
The aforementioned parameters are entered into the MATLAB (MATLAB 9.0, R2016a, MathWorks, Natick, MA, USA) for solving in order to obtain the Pareto front.Nine solutions are included on the Pareto front, as shown in Table 4.The Pareto optimal solution set cannot explain further how to choose a scenario as the final decision on the Pareto front.For any non-inferior solution of multi-objective models, there is no other possibility that the two objective function values corresponding to the non-inferior solutions can be improved from the two objective functions corresponding to other non-inferior solutions.Policymakers can sacrifice only one goal in a certain range in return for another goal improvement.They show their preferences in the alternative tradeoff between such goals.After one operation, a number of Pareto optimal solutions are obtained under the AWOA, which results in higher solution efficiency, and provides greater choices for policymakers.Table 4 shows that if there is a special preference for the minimum amount of water, we can choose scenario 1; if there is a special preference for the best economic benefits, we can choose scenario 9; and if we consider the minimum amount of water shortage and the best economic benefits, we can choose one of the scenarios from 2 to 8. As there is a serious shortage of water resources in Handan, meeting the water demand of different water sectors is an urgent challenge, which must be solved to improve current water resources management.Thus, scenario 1 better reflects the fairness of water consumption among different households.Therefore, it is chosen as an effective scenario, and this study analyzes it in detail.The optimal water resource allocation results in each sub-region in Handan are shown in Figure 3.The total water demand for all users in the region is 2726.08 × 10 6 m 3 , and the total allocated volume is 2321.74× 10 6 m 3 , resulting in a total water shortage of 404.34 × 10 6 m 3 and a shortage rate of 14.8%.The four sub-regions with relatively high shortage rates are Weixian, Quzhou, Cheng'an, and Handan city, with rates of 43.3%, 38.5%, 32.4%, and 22.3%, respectively.In Shexian, Fengfeng, Cixian, Jize, and Feixiang, the allocated volume can meet the water demand, and there is no water shortage.The allocated volumes by local surface water, groundwater, water transferred through the South-North Water Transfer Project, the Dongwushi reservoir, the Yuecheng reservoir, imported water from the Yellow River, and recycled water in Handan are 153.06× 10 6 m 3 , 1007.32 × 10 6 m 3 , 259.11 × 10 6 m 3 , 221.79 × 10 6 m 3 , 197.56 × 10 6 m 3 , 100.00 × 10 6 m 3 , and 382.90 × 10 6 m 3 , respectively.They account for 6.6%, 43.4%, 11.2%, 9.6%, 8.5%, 4.3%, and 16.5%, respectively, of the total water allocation in the region.It is noteworthy that groundwater resources are the main water source of Handan.Figure 4 illustrates the urban domestic water allocation in each sub-region in Handan.The total urban domestic water demand is 387.16 × 10 6 m 3 while the total allocated volume is 387.16 × 10 6 m 3 .The water supply can fully satisfy the urban domestic water demand in all sub-regions.The water demand is mainly supported by groundwater and water transferred through the South-North Water Transfer Project, the allocated volumes of which are 263.07× 10 6 m 3 and 124.09 × 10 6 m 3 , respectively, accounting for 67.9% and 32.1%, respectively, of the total supply.It is noteworthy that groundwater is the dominant water source for urban domestic water demand.Figure 4 illustrates the urban domestic water allocation in each sub-region in Handan.The total urban domestic water demand is 387.16 × 10 6 m 3 while the total allocated volume is 387.16 × 10 6 m 3 .The water supply can fully satisfy the urban domestic water demand in all sub-regions.The water demand is mainly supported by groundwater and water transferred through the South-North Water Transfer Project, the allocated volumes of which are 263.07× 10 6 m 3 and 124.09 × 10 6 m 3 , respectively, accounting for 67.9% and 32.1%, respectively, of the total supply.It is noteworthy that groundwater is the dominant water source for urban domestic water demand.Figure 4 illustrates the urban domestic water allocation in each sub-region in Handan.The total urban domestic water demand is 387.16 × 10 6 m 3 while the total allocated volume is 387.16 × 10 6 m 3 .The water supply can fully satisfy the urban domestic water demand in all sub-regions.The water demand is mainly supported by groundwater and water transferred through the South-North Water Transfer Project, the allocated volumes of which are 263.07× 10 6 m 3 and 124.09 × 10 6 m 3 , respectively, accounting for 67.9% and 32.1%, respectively, of the total supply.It is noteworthy that groundwater is the dominant water source for urban domestic water demand.Figure 5 gives the rural domestic water allocation in each sub-region in Handan.The rural domestic water supply is able to meet the demand in all sub-regions.The total rural domestic water demand is 35.91 × 10 6 m 3 while the total allocated volume is 35.91 × 10 6 m 3 , supported wholly by groundwater.Figure 6 shows the primary sector water allocation in each sub-region in Handan.The total water demand of the primary sector in the region is 1591.82× 10 6 m 3 , whereas the total allocated volume is 1210.66× 10 6 m 3 , with a shortage of 381.15 × 10 6 m 3 , and an average shortage rate of 23.9%.The four sub-regions with relatively great shortages are Handan city, Weixian, Quzhou, and Cheng'an and their shortage rates are 62.4%, 60.3%, 50.4%, and 47.5%, respectively.In Shexian, Wu'an, Fengfeng, Cixian, Jize, and Feixiang, water supply can meet the primary sector water demand.The allocated volume supported by local surface water, groundwater, the Dongwushi reservoir, the Yuecheng reservoir, imported water from the Yellow River, and recycled water are 127.18× 10 6 m 3 , 477.69 × 10 6 m 3 , 189.08 × 10 6 m 3 , 165.70 × 10 6 m 3 , 85.83 × 10 6 m 3 , and 165.19 × 10 6 m 3 , respectively, accounting for 10.5%, 39.5%, 15.6%, 13.7%, 7.1%, and 3.6% of the total supply, respectively.Again, groundwater is the main water source for the primary sector.The secondary sector water allocation in each sub-region in Handan is shown in Figure 7.The total water demand of the secondary sector is 307.06 × 10 6 m 3 , and the total allocated volume is 292.60 × 10 6 m 3 , leading to a shortage of 14.46 × 10 6 m 3 , and an average shortage rate of 4.7%.In all Figure 6 shows the primary sector water allocation in each sub-region in Handan.The total water demand of the primary sector in the region is 1591.82× 10 6 m 3 , whereas the total allocated volume is 1210.66× 10 6 m 3 , with a shortage of 381.15 × 10 6 m 3 , and an average shortage rate of 23.9%.The four sub-regions with relatively great shortages are Handan city, Weixian, Quzhou, and Cheng'an and their shortage rates are 62.4%, 60.3%, 50.4%, and 47.5%, respectively.In Shexian, Wu'an, Fengfeng, Cixian, Jize, and Feixiang, water supply can meet the primary sector water demand.The allocated volume supported by local surface water, groundwater, the Dongwushi reservoir, the Yuecheng reservoir, imported water from the Yellow River, and recycled water are 127.18× 10 6 m 3 , 477.69 × 10 6 m 3 , 189.08 × 10 6 m 3 , 165.70 × 10 6 m 3 , 85.83 × 10 6 m 3 , and 165.19 × 10 6 m 3 , respectively, accounting for 10.5%, 39.5%, 15.6%, 13.7%, 7.1%, and 3.6% of the total supply, respectively.Again, groundwater is the main water source for the primary sector.  Figure 6 shows the primary sector water allocation in each sub-region in Handan.The total water demand of the primary sector in the region is 1591.82× 10 6 m 3 , whereas the total allocated volume is 1210.66× 10 6 m 3 , with a shortage of 381.15 × 10 6 m 3 , and an average shortage rate of 23.9%.The four sub-regions with relatively great shortages are Handan city, Weixian, Quzhou, and Cheng'an and their shortage rates are 62.4%, 60.3%, 50.4%, and 47.5%, respectively.In Shexian, Wu'an, Fengfeng, Cixian, Jize, and Feixiang, water supply can meet the primary sector water demand.The allocated volume supported by local surface water, groundwater, the Dongwushi reservoir, the Yuecheng reservoir, imported water from the Yellow River, and recycled water are 127.18× 10 6 m 3 , 477.69 × 10 6 m 3 , 189.08 × 10 6 m 3 , 165.70 × 10 6 m 3 , 85.83 × 10 6 m 3 , and 165.19 × 10 6 m 3 , respectively, accounting for 10.5%, 39.5%, 15.6%, 13.7%, 7.1%, and 3.6% of the total supply, respectively.Again, groundwater is the main water source for the primary sector.The secondary sector water allocation in each sub-region in Handan is shown in Figure 7.The total water demand of the secondary sector is 307.06 × 10 6 m 3 , and the total allocated volume is 292.60 × 10 6 m 3 , leading to a shortage of 14.46 × 10 6 m 3 , and an average shortage rate of 4.7%.In all The secondary sector water allocation in each sub-region in Handan is shown in Figure 7.The total water demand of the secondary sector is 307.06 × 10 6 m 3 , and the total allocated volume is 292.60 × 10 6 m 3 , leading to a shortage of 14.46 × 10 6 m 3 , and an average shortage rate of 4.7%.In all sub-regions except Wu'an, the supply meets the demand.In Wu'an, the demand is 71.00 × 10 6 m 3 , and the allocated volume is 56.54 × 10 6 m 3 .There is a shortage of 14.46 × 10 6 m 3 .Water sources for the secondary sector water supply include groundwater, water transferred through the South-North Water Transfer Project, the Dongwushi reservoir, the Yuecheng reservoir, imported water from the Yellow River, and recycled water.Their allocated volumes are 81.02× 10 6 m 3 , 71.26 × 10 6 m 3 , 14.03 × 10 6 m 3 , 19.35 × 10 6 m 3 , 8.40 × 10 6 m 3 , and 98.55 × 10 6 m 3 , and contribute to total supply by 26.4%, 23.2%, 4.6%, 6.3%, 2.7%, and 32.1%, respectively.The main water source for the secondary sector is recycled water.
Water 2018, 10, 87 15 of 27 sub-regions except Wu'an, the supply meets the demand.In Wu'an, the demand is 71.00 × 10 6 m 3 , and the allocated volume is 56.54 × 10 6 m 3 .There is a shortage of 14.46 × 10 6 m 3 .Water sources for the secondary sector water supply include groundwater, water transferred through the South-North Water Transfer Project, the Dongwushi reservoir, the Yuecheng reservoir, imported water from the Yellow River, and recycled water.Their allocated volumes are 81.02× 10 6 m 3 , 71.26 × 10 6 m 3 , 14.03 × 10 6 m 3 , 19.35 × 10 6 m 3 , 8.40 × 10 6 m 3 , and 98.55 × 10 6 m 3 , and contribute to total supply by 26.4%, 23.2%, 4.6%, 6.3%, 2.7%, and 32.1%, respectively.The main water source for the secondary sector is recycled water.Figure 8 shows the tertiary sector water allocation in each sub-region in Handan.As shown in Figure 8, the total water demand and allocated volume are 213.41× 10 6 m 3 and 213.41 × 10 6 m 3 , respectively.In all sub-regions, the supply can satisfy the demand.The tertiary sector water demand is wholly supported by groundwater and water transferred through the South-North Water Transfer Project, the supplies of which are 8.40 × 10 6 m 3 and 98.55 × 10 6 m 3 , respectively, and account for 70.1% and 29.9% of the total, respectively.Groundwater is the dominant water source for the tertiary sector.  Figure 8 shows the tertiary sector water allocation in each sub-region in Handan.As shown in Figure 8, the total water demand and allocated volume are 213.41× 10 6 m 3 and 213.41 × 10 6 m 3 , respectively.In all sub-regions, the supply can satisfy the demand.The tertiary sector water demand is wholly supported by groundwater and water transferred through the South-North Water Transfer Project, the supplies of which are 8.40 × 10 6 m 3 and 98.55 × 10 6 m 3 , respectively, and account for 70.1% and 29.9% of the total, respectively.Groundwater is the dominant water source for the tertiary sector.
Water 2018, 10, 87 15 of 27 sub-regions except Wu'an, the supply meets the demand.In Wu'an, the demand is 71.00 × 10 6 m 3 , and the allocated volume is 56.54 × 10 6 m 3 .There is a shortage of 14.46 × 10 6 m 3 .Water sources for the secondary sector water supply include groundwater, water transferred through the South-North Water Transfer Project, the Dongwushi reservoir, the Yuecheng reservoir, imported water from the Yellow River, and recycled water.Their allocated volumes are 81.02× 10 6 m 3 , 71.26 × 10 6 m 3 , 14.03 × 10 6 m 3 , 19.35 × 10 6 m 3 , 8.40 × 10 6 m 3 , and 98.55 × 10 6 m 3 , and contribute to total supply by 26.4%, 23.2%, 4.6%, 6.3%, 2.7%, and 32.1%, respectively.The main water source for the secondary sector is recycled water.Figure 8 shows the tertiary sector water allocation in each sub-region in Handan.As shown in Figure 8, the total water demand and allocated volume are 213.41× 10 6 m 3 and 213.41 × 10 6 m 3 , respectively.In all sub-regions, the supply can satisfy the demand.The tertiary sector water demand is wholly supported by groundwater and water transferred through the South-North Water Transfer Project, the supplies of which are 8.40 × 10 6 m 3 and 98.55 × 10 6 m 3 , respectively, and account for 70.1% and 29.9% of the total, respectively.Groundwater is the dominant water source for the tertiary sector.Water 2018, 10, 87 16 of 29 Figure 9 shows the ecological water allocation in each sub-region in Handan.The total water demand and total allocated volume are 190.72 × 10 6 m 3 and 182.00 × 10 6 m 3 , respectively.The shortage is 8.72 × 10 6 m 3 , and the shortage rate is 4.6%.In all sub-regions, except Handan city, the supply can satisfy the demand.In Handan city, the demand is 114.84 × 10 6 m 3 , while the allocated volume is 106.12 × 10 6 m 3 , resulting in a deficiency of 8.72 × 10 6 m 3 .The supplies by local surface water, the Dongwushi reservoir, the Yuecheng reservoir, imported water from the Yellow River, and recycled water are 25.88 × 10 6 m 3 , 18.68 × 10 6 m 3 , 12.50 × 10 6 m 3 , 5.77 × 10 6 m 3 , and 119.17 × 10 6 m 3 , respectively, accounting for 14.2%, 10.3%, 6.9%, 3.2%, and 65.5%, respectively, of the total.Recycled water is the main source for the ecological water demand.
Water 2018, 10, 87 16 of 27 is 8.72 × 10 6 m 3 , and the shortage rate is 4.6%.In all sub-regions, except Handan city, the supply can satisfy the demand.In Handan city, the demand is 114.84 × 10 6 m 3 , while the allocated volume is 106.12 × 10 6 m 3 , resulting in a deficiency of 8.72 × 10 6 m 3 .The supplies by local surface water, the Dongwushi reservoir, the Yuecheng reservoir, imported water from the Yellow River, and recycled water are 25.88 × 10 6 m 3 , 18.68 × 10 6 m 3 , 12.50 × 10 6 m 3 , 5.77 × 10 6 m 3 , and 119.17 × 10 6 m 3 , respectively, accounting for 14.2%, 10.3%, 6.9%, 3.2%, and 65.5%, respectively, of the total.Recycled water is the main source for the ecological water demand.Table 5 lists the volumes of local surface water allocated to each user in each sub-region in Handan.Due to the water quality, local surface water is allocated only to the primary sector and ecology.The corresponding allocated volumes are 127.18× 10 6 m 3 and 25.88 × 10 6 m 3 , respectively, and the total is 153.06 × 10 6 m 3 .Since the total available water supply by local surface water in the region is 236.58 × 10 6 m 3 , 83.52 × 10 6 m 3 of water not utilized.From Figures 6 and 9, local surface water accounts for 10.5% and 14.2% of the total volumes allocated to the primary sector and ecology, respectively.In Wu'an, the highest volumes of local surface water are allocated to the primary sector and ecology.The corresponding volumes are 80.68 × 10 6 m 3 and 13.93 × 10 6 m 3 , which account for 63.4% and 53.8%, respectively, of the total volumes of local surface water allocated.This is because Wu'an has abundant local surface water, or 48.6% of the total in Handan.Table 5 lists the volumes of local surface water allocated to each user in each sub-region in Handan.Due to the water quality, local surface water is allocated only to the primary sector and ecology.The corresponding allocated volumes are 127.18× 10 6 m 3 and 25.88 × 10 6 m 3 , respectively, and the total is 153.06 × 10 6 m 3 .Since the total available water supply by local surface water in the region is 236.58 × 10 6 m 3 , 83.52 × 10 6 m 3 of water not utilized.From Figures 6 and 9, local surface water accounts for 10.5% and 14.2% of the total volumes allocated to the primary sector and ecology, respectively.In Wu'an, the highest volumes of local surface water are allocated to the primary sector and ecology.The corresponding volumes are 80.68 × 10 6 m 3 and 13.93 × 10 6 m 3 , which account for 63.4% and 53.8%, respectively, of the total volumes of local surface water allocated.This is because Wu'an has abundant local surface water, or 48.6% of the total in Handan.Table 6 shows the volumes of groundwater allocated to each user in each sub-region in Handan.Groundwater is allocated to urban domestic, rural domestic, primary sector, secondary sector, and tertiary sector uses in volumes of 263.07 × 10 6 m 3 , 35.91 × 10 6 m 3 , 477.69 × 10 6 m 3 , 81.02 × 10 6 m 3 , and 149.63 × 10 6 m 3 , respectively.These account for 26.1%, 3.6%, 47.4%, 8.0%, and 14.9% of the total volume of groundwater allocated, respectively.Groundwater is mostly supplied to the primary sector.This is in line with actual local groundwater usage.The total available water supply by groundwater and total allocated volume are 1130.42× 10 6 m 3 and 1007.32 × 10 6 m 3 , respectively.There is 123.10 × 10 6 m 3 of unallocated water remaining.The allocation results of water transferred through the South-North Water Transfer Project to each user in each sub-region in Handan are listed in Table 7.The transferred water is mostly supplied for urban domestic, secondary, and tertiary sector uses.The total available water supply and total allocated volumes are 388.62 × 10 6 m 3 and 259.11 × 10 6 m 3 , respectively, resulting in a surplus of

Discussion
The WOA, proposed by Mirjalili and Lewis in 2016, has less application to water resources management The WOA was used to solve long-term optimal operation of a single reservoir and cascade reservoirs by Cui Donglai [52] in 2017.The results were compared with those obtained from six algorithms, including the cuckoo search algorithm, differential evolution algorithm, shuffled frogleaping algorithm, PSO algorithm, firefly algorithm, and SCE-UA algorithm.The results show that the WOA for the optimal operation of a single reservoir and cascade reservoirs is superior to the other six algorithms.The WOA is feasible and effective for reservoir optimal operation [52].This study applied the WOA to the multi-objective optimization of water resources allocation model for the first time in the literature.When the number of iterations is set to 180 times, the algorithm cannot converge, and attempts to run 300 times, 500 times, 1000 times, 1500 times, and 2000 times while the WOA still cannot converge.In order to solve this problem, an AWOA is proposed in this study, where by logistic mapping is used to initialize swarm location, and inertia weighting is employed to strengthen local search capability.It can be concluded from Figure 11 and Table 13 that after the AWOA is applied to the multi-objective optimization of water resources allocation model in Handan, it speeds up the algorithm convergence and improves the accuracy of the algorithm.
The AWOA is compared with the PSO to verify which is superior when applied to the multi-objective water resources allocation model.As a mature swarm-based algorithm, PSO has obtained a large number of successful applications to the optimal allocation of water resources [29,30,[53][54][55][56][57][58][59][60].According to previous research results, efficiency of the water resources allocation problem can be effectively improved under PSO.The AWOA out performs the PSO in terms of convergence rate and solution accuracy through Figures 11 and 13, and Table 13.The AWOA provides a new way to solve the complex multi-objective optimization of water resources allocation model.

Conclusions
Due to population growth and rapid socio-economic development, regions with water shortages are expanding, making shortages increasingly severe in China.This has developed into a critical issue that significantly hinders social and economic development.Alleviating water shortages by optimizing water resource allocation has become increasingly important.Water resource allocation optimization is a multi-objective problem, which cannot be solved easily with traditional programming methods.Based on previous research, this study introduces the WOA into multi-objective water resource allocation optimization models.To overcome drawbacks like relatively low convergence precision and convergence when applying the WOA algorithm to complex optimization problems, logistic mapping is used to initialize swarm location, and enhance the quality of initialization swarm location.In addition, inertia weighting is employed to strengthen the local search capability.Compared to the WOA and PSO algorithm, the AWOA shows substantially enhanced convergence rates and precision.
The AWOA is utilized to solve the Handan water resource allocation optimization model.The results indicate that the total water demand for all users in the city is 2726.08 × 10 6 m 3 , and the total volume is 2321.74× 10 6 m 3 , resulting in a total water shortage of 404.34 × 10 6 m 3 , and a shortage rate of 14.8%.Water shortages occur mainly in the primary sector (agriculture), with a shortage volume of 381.15 × 10 6 m 3 .This accounts for 94.3% of the total shortage.The allocation optimization results show good agreement with the reality of water resource exploitation and utilization.This verifies the relatively high reliability and applicability of the AWOA in solving multi-objective multi-source multi-user water resource allocation optimization models.
The AWOA proposed in this study can effectively solve multi-objective water resource allocation optimization models.It provides a Pareto non-inferior solution set and hence, several reliable multi-source joint control plans.According to the preferences for the objective functions, the authority can select the appropriate plan to control and manage multiple water sources jointly.
Since the sewage disposal data are not complete, this study establishes only a water resource allocation optimization model that aims for the greatest economic benefit and smallest shortage volume.In future, with complete sewage disposal data, a water resource allocation optimization model targeting relatively optimal economic, environmental, and social benefits should be constructed and solved using the AWOA.This could verify the applicability of the AWOA to more multi-objective water resource allocation optimization problems.
value obtained by all individuals in the population in the i-th objective function, and min m f represents the minimum value obtained by all individuals in the population in the i-th objective function.For each objective function, after all individuals in the population have undergone fast non-dominated sorting, the crowding distance of two individuals at the boundary of each level is infinite.

Figure 1 .
Figure 1.Crowding degree of individual i.  Crowded-comparison operator After carrying out the fast non-dominated sorting operations and calculation of crowding degree for each individual i in the population, each individual acquires two attributes, rank i (non-dominated

and 1 
 , the current position of individual whales are updated with Equation (14); when 1   , a random individual whale rand X  is selected, and the current positions of individual whales are updated with Equations (10) and (11). 5.When 0.5 p  , the current positions of individual whales are updated with Equation (15).

Figure 1 .
Figure 1.Crowding degree of individual i.

Figure 2 .
Figure 2. The map of administrative division and water supply system in Handan.

Figure 2 .
Figure 2. The map of administrative division and water supply system in Handan.

Figure 5 Figure 3 .
Figure5gives the rural domestic water allocation in each sub-region in Handan.The rural domestic water supply is able to meet the demand in all sub-regions.The total rural domestic water

Figure 5 Figure 4 .
Figure5gives the rural domestic water allocation in each sub-region in Handan.The rural domestic water supply is able to meet the demand in all sub-regions.The total rural domestic water

Figure 9 Figure 7 .
Figure9shows the ecological water allocation in each sub-region in Handan.The total water demand and total allocated volume are 190.72 × 10 6 m 3 and 182.00 × 10 6 m 3 , respectively.The shortage

Figure 9 Figure 8 .
Figure9shows the ecological water allocation in each sub-region in Handan.The total water demand and total allocated volume are 190.72 × 10 6 m 3 and 182.00 × 10 6 m 3 , respectively.The shortage
indicates that the administrative division can be allocated a certain amount of water.

Table 3 .
Coefficients of water supply sequence for different water users.

Table 4 .
Results of Pareto front under AWOA.

Table 5 .
Results for local surface water allocated to every user (10 6 m 3 ).

Table 5 .
Results for local surface water allocated to every user (10 6 m 3 ).

Table 6 .
Results for groundwater allocated to every user (10 6 m 3 ).

Table 13 .
Results for objective function values of AWOA, WOA, and PSO.Water allocation based on three algorithms for each sub-region in Handan.