Comparison of Transboundary Water Resources Allocation Models Based on Game Theory and Multi-Objective Optimization

: Transboundary water resources allocation is an effective measure to resolve water-related conﬂicts. Aiming at the problem of water conﬂicts, we constructed water resources allocation models based on game theory and multi-objective optimization, and revealed the differences between the two models. We compare the Pareto front solved by the AR-MOEA method and the NSGA-II method, and analyzed the difference between the Nash–Harsanyi Leader–Follower game model and the multiobjective optimization model. The Huaihe River basin was selected as a case study. The results show that: (1) The AR-MOEA method is better than the NSGA-II method in terms of the diversity metric ( ∆ ); (2) the solution of the asymmetric Nash–Harsanyi Leader–Follower game model is a non-dominated solution, and the asymmetric game model can obtain the same water resources allocation scheme of the multi-objective optimal allocation model under a speciﬁc preference structure; (3) after the multi-objective optimization model obtains the Pareto front, it still needs to construct the preference information of the Pareto front for a second time to make the optimal solution of a multi-objective decision, while the game model can directly obtain the water resources allocation scheme at one time by participating in the negotiation. The results expand the solution method of water resources allocation models and provide support for rational water resources allocation.


Introduction
Transboundary water resources allocation usually involves multi-stakeholders. When water resources are insufficient to meet the water demand of the whole basin, water conflicts may arise among multi-stakeholders. In the case of water resources shortage, transboundary water resources allocation become an important means to solve water conflicts. Transboundary water resources allocation is a kind of complex decision-making problem with nonlinear, multi-objective and multi-stage characteristics.
In transboundary water resources allocation based on game theory, the way various stakeholders directly participate in the negotiation makes the water resources allocation results more acceptable. Roger [1] originally applied game theory to solve water resources conflicts. Since then, water resources allocation based on game theory has been studied. Bogardi and Szidarovsky [2] applied the oligopoly game to different water resources management problems. Parrachino et al. [3] applied cooperative game theory to solve water conflicts. Carraro et al. [4][5][6] systematically described the application of non-cooperative negotiation theory in water resources conflicts. Madani et al. [7] demonstrated that the application of water resources allocation based on game theory can be divided into five The rest of this paper is organized as follows: Section 2 describes the material and methods. Section 3 describes the results of the asymmetric Nash-Harsanyi Leader-Follower game model and multi-objective optimization model based on the AR-MOEA method in the Huaihe River basin. Section 4 discusses the comparison of water resources allocation results based on the NSGA-II method and the AR-MOEA method, and reveals the differences between water resources allocation based on the asymmetric Nash-Harsanyi Leader-Follower game model and the multi-objective optimization model based on the AR-MOEA method. Section 5 concludes the paper.

Study Area
The Huaihe River basin located in the east (longitude from 111.9° to 121.4°) and in the north (latitude from 30.9° to 36.6°) of China, with a total area of 19,000 km 2 . The Huaihe River basin contains Henan province, Anhui province, and Jiangsu province (three provinces). The problem of water resources shortage is prominent, and the inter-provincial water resources conflict has become a serious problem. Though the study area is not international, an inter-provincial basin is effectively equivalent to an international basin as long as its boundaries do not match political boundaries. The stakeholders in the process of the Huaihe River basin water resources allocation are the watershed management agency, Henan province, Anhui province, and Jiangsu province. The spatial distribution and location of the associated province of Huaihe River basinis described in Figure 1.

Transboundary Water Resources Allocation Based on Asymmetric Nash-Harsanyi Leader-Follower Game Model and Its Solution
Transboundary water resources allocation based on the asymmetric Nash-Harsanyi Leader-Follower game model is proposed by Fu et al. [21]. The game model considering fairness, efficiency, and sustainable principles can not only guarantee the basic water demand, but also better balance the economic development levels among followers. Therefore, the game model is used to compare with the multi-objective optimization model.
Transboundary water resources allocation based on the asymmetric Nash-Harsanyi Leader-Follower game model is a combination of the leader-follower concepts and the Nash-Harsanyi theory. This model involves a single leader and multiple followers. In this

Transboundary Water Resources Allocation Based on Asymmetric Nash-Harsanyi Leader-Follower Game Model and Its Solution
Transboundary water resources allocation based on the asymmetric Nash-Harsanyi Leader-Follower game model is proposed by Fu et al. [21]. The game model considering fairness, efficiency, and sustainable principles can not only guarantee the basic water demand, but also better balance the economic development levels among followers. Therefore, the game model is used to compare with the multi-objective optimization model.
Transboundary water resources allocation based on the asymmetric Nash-Harsanyi Leader-Follower game model is a combination of the leader-follower concepts and the Nash-Harsanyi theory. This model involves a single leader and multiple followers. In this study, a watershed management agency is considered to be the leader, who is responsible for the eco-environmental benefit of the whole basin. The associated areas serve as the followers. The eco-environmental water used by the whole basin in the allocation process is the public water. If the public water allocation is unreasonable, the eco-environmental situation will deteriorate, and sustainable development is challenged. Therefore, it is necessary to rationally distribute public water. The transboundary river water resources allocation process can be divided into two stages. The first stage involves allocating the public water. In the second stage, followers need to allocate water resources according to the asymmetric Nash-Harsanyi game model after the decisions are given by the watershed management agency.
The model of the watershed management agency (leader) can be described as the following: Subject to: Water balance constraints: Public water constraints: where V is the objective function of the watershed management agency, s is the public water and the solution of the leader, s * is the ideal eco-environmental water demand, −(s − s * ) 2 is the eco-environmental benefit, w = {w 1 , w 2 , · · · , w n } is the decision vector of water allocated to followers, w i is the water in the ith follower(i = 1, 2, · · · , n), Q is the total allocable water, and λ i is the minimum survival water demand in the ith follower. The model of the regions (followers) can be described as the following: Subject to: Optimal solution existence constraints: Water balance constraints: Water 2021, 13, 1421

of 18
The bargaining weight constraints: Water resources allocation constraints in each region: where u i (w i ) is the objective function for follower i, d i is the disagreement point for follower i, α i is the bargaining weight for follower i, u i (I i ) is the benefit of minimum water allocation for follower i, r i is the water demand for follower i, I i is the minimum water allocation for follower i, η is the weight of the equity principle, δ 1i is the bargaining weight for follower i when only considering the equity principle that represents the equity of water for follower i, (1 − η) is the weight of efficiency principle, δ 2i is the bargaining weight for follower i when only the efficiency principle is considered, β i is the water use correction coefficient for follower i, D i is the integrated water utility for follower i, r i,k is the water demand of the kth water user in the ith follower, and D i,k is the integrated water utility of the kth water user in the ith follower.
The solution of the game model can refer to Fu's article [21]. After the watershed management agency obtains the public water, the followers can obtain the water resources allocation results by using the successive linear programming method.
The successive linear programming method was first proposed in 1961 [23]. It is an iterative optimization method using linear programming technology to solve nonlinear problems. At present, this method has been widely used in practical engineering and achieved good results. The steps for solving successive linear programming are as follows:

1.
Let x 0 be the initial feasible solution of the optimization problem, and the corresponding objective function value is f (x 0 ); 2.
The nonlinear objective function is linearized by using the first order descriptive form point (x 0 , f (x 0 )) of Taylor series, and the optimal solution (x 1 , f (x 1 )) of linearized objective function in the neighborhood (x 0 − δ, x 0 − δ) is obtained by using the linear programming method; 3.
Taking the optimal solution (x 1 , f (x 1 )) as a new feasible solution, the nonlinear objective function is linearized again by using the first order description point (x 1 , f (x 1 )) of Taylor series, and the optimal solution (x 2 , f (x 2 )) of linearized objective function in the neighborhood (x 1 − δ, x 1 − δ) is obtained by using the linear programming method; 4.
Repeat the above step 3 for the optimal solution and continue to iterate until the obtained optimal solution and the optimal objective function meet the condition of iteration termination.

Transboundary Water Resources Allocation Based on Multi-Objective Optimization Model
Transboundary water resources allocation is a multi-stakeholders utilization conflict issue, involving watershed management agency and multiple areas.
The watershed management agency is responsible for the eco-environmental water utilization of the entire basin in the process of water resources allocation. The optimization goal of a watershed management agency is to maximize the eco-environmental benefit. The specific formula is as follows: Water 2021, 13, 1421 6 of 18 The optimization goal of area i is to maximize the comprehensive benefit of the area. The specific calculation formula is as follows: Constraints conditions for n + 1 targets: Water balance constraints: Water resources allocation constraints in each region: In the multi-objective optimization model of transboundary water resources allocation, the ecological environmental water used in the entire basin is public water. If the water resources allocation result is not reasonable, the eco-environmental situation is severe, and sustainable development is faced with challenges.
In the process of actual water resources allocation, due to the lack of consideration of eco-environmental water requirements, the eco-environmental water consumption is largely consumed by the production of water. This could not meet the sustainable development of water resources, which brings about water resources and eco-environmental problems, so the eco-environmental benefit of the watershed management agency needs to be prioritized. Therefore, the objective function of the watershed management agency in transboundary water resources allocation can be converted into constraint conditions, and the n + 1 objective function converted into n objective function.

Multi-Objective Optimization Method
(1) AR-MOEA method. The AR-MOEA method is an evolutionary algorithm by Tian, which can be used to solve multi-objective optimization problems [22]. This method uses the enhanced inverted generational distance indicator as the main environmental selection strategy. In this method, a uniformly distributed group of reference points is used as the reference points set for calculating the enhanced inverted generational distance indicator. Therefore, the AR-MOEA method can be applied to a different front, and the distribution of reference points can be adjusted automatically according to the current population distribution in the calculation process.
The algorithm steps of AR-MOEA are as follows: 1.
An initial population P l of size M is randomly generated, save the population P l to an external document E, and multi-objective optimization method based on objective decomposition is used to generated a set of reference points R with a scale of M R ; 2.
A mating pool selection strategy based on the enhanced inverted generational distance with noncontributing solution detection is used to select the initial population P l , mating pool p l is selected and the offspring population P l is generated; 3.
The archive E and the adapted reference point set R are updated through the offspring population P l , then the next generation population P l+1 is selected by the enhanced inverted generational distance with a noncontributing solution detection selection strategy; 4.
The above steps 1, 2, and 3 are repeated for the new initial population P l+1 until a termination criterion is reached.
The AR-MOEA method contains four solutions sets: initial population P l , initial reference point set R, archive set E, and adapted reference point set R . The population P l contains the candidate solutions as final output, the initial reference point set R is used to guarantee uniform distribution of the candidate solution in P l , and the archive set E provides the shape of the Pareto front of the adapted reference point set R . An enhanced termination criterion is reached.
The AR-MOEA method contains four solutions sets: initial population l P , initial reference point set R , archive set E , and adapted reference point set ' R . The population l P contains the candidate solutions as final output, the initial reference point set R is used to guarantee uniform distribution of the candidate solution in l P , and the archive set E provides the shape of the Pareto front of the adapted reference point set ' R . An enhanced inverted generational distance indicator is used to select populations l P to change the adapted reference point set ' R . The relationships between the four solution sets are shown in Figure 2.  (2) NSGA-II method. The NSGA-II method is a multi-objective genetic algorithm proposed by Deb, which can be used to solve multi-objective optimization problems [24]. The NSGA-II method has the advantages of reducing the complexity of the algorithm, fast operation speeds, and good distribution of the solution set, so the method has become a widely used multi-objective genetic algorithm.
The algorithm steps of AR-MOEA are as follows: 1.
An initial population p l of size M is randomly generated, and the initial population is stratified by non-inferiority, and the genetic operator (selection operator, crossover operator, mutation operator) is used to obtain an offspring population q l . Then, the initial population p l and the offspring population q l are mixed together to form a new population N l of size 2M; 2.
The new population N l of size 2M is rapidly sorted to form a non-dominated set E l , then the crowding degree of the individuals in the non-dominated set E l is calculated. According to the non-dominant relationship and the crowding degree of individuals, a new initial population p l+1 with the appropriate individual composition size of M is selected; 3.
The above steps 1 and 2 are repeated for the new initial population p l+1 until a termination criterion is reached. Figure 3 shows the specific calculation flow chart of NSGA-II.
(3) Linear weighted sum method. The linear weighted sum method transforms the multi-objective optimization problem into a single objective problem to solve. The specific solving steps are as follows: According to the importance of each goal, a set of non-negative weights τ i (i = 1, 2, · · · , n) corresponding to the objective function F i is given, The non-negative weight of the objective function is closely related to the preference structure information of the decision maker. Each objective function is multiplied by its corresponding non-negative weight, then added to obtain a single objective function. The formula for the linear weighted sum method can be written as the following: Subject to: Inequality constraints: Equality constraints: The linear weighted sum method can be solved by the successive linear programming method. According to the characteristics of the solution set measured by the performance metrics, the existing algorithm performance metrics can be divided into three categories: (1) Convergence metrics-evaluate the approximation degree between the true Pareto front and non-dominated solution.

Performance Metrics for Multi-Objective Optimization Evolutionary Algorithms
According to the characteristics of the solution set measured by the performance metrics, the existing algorithm performance metrics can be divided into three categories: (1) Convergence metrics-evaluate the approximation degree between the true Pareto front and non-dominated solution.
Generational distance (GD) can be used to measure the convergence of non-dominated solutions obtained by multi-objective evolutionary algorithms [25]. The smaller the GD is, the better the convergence of non-dominated solutions obtained by multi-objective evolutionary algorithms is. It is assumed that P * is a set of solutions uniformly distributed on the real Pareto front and A is the non-dominated solution set obtained by multi-objective evolutionary algorithms. The expression of GD can be defined as: where d(x, y) is the Euclidean distance between point x in P * and point y in A, and |A| is the number of non-dominated solutions.
(2) Diversity metrics-evaluate scatter of non-dominated solutions. Diversity metrics can be subdivided into the evenness and spread of the non-dominated solutions in the whole feasible region. Evenness in diversity metrics reflects the evenness of distribution in non-dominated solutions. Spread in diversity metrics reflects the spread degree of solution set in target space.
A range metric (∆) is employed to measure the diversity of non-dominated solutions obtained by the multi-objective evolutionary algorithm [24]. The mathematical definition of diversity metrics is described as: where d f and d l are the Euclidean distance between the extreme non-dominated solution and the boundary solutions of the obtained non-dominated solution set, N is the number of non-dominated solutions obtained by multi-objective evolutionary algorithms, d i is the distance between continuous solutions in obtained non-dominated solution set, and d is the average of all d i (i = 1, 2, . . . , N − 1). When ∆ is zero, it indicates that the non-dominated solution set calculated by the multi-objective optimization algorithm has good diversity.
(3) Composite metrics considering convergence and diversity. Inverted generational distance (IGD) and hypervolume (HV) are employed to measure the convergence and diversity of a non-dominated solution obtained by the multi-objective evolutionary algorithm.
IGD can be used to measure the mean of the minimum distance between point x in P * and point y in A. The smaller IGD is, the better the convergence and diversity of non-dominated solutions obtained by multi-objective evolutionary algorithms is. The expression of IGD can be described as the following: where d(x, y) is the Euclidean distance between point x in P * and point y in A, |P * | is the number of true Pareto solutions. HV is the volume in the target space region enclosed by the non-dominated solutions A and the reference point set R [25]. The larger the HV, the better the convergence and diversity of non-dominated solutions obtained by multi-objective evolutionary algorithms. The calculation formula of HV is: where volume(·) is Lebesgue measure, N is the number of non-dominated solutions obtained by multi-objective evolutionary algorithms, and v i is the volume enclosed between the ith non-dominated solution and the reference point. We selected the above representative performance metrics to evaluate the AR-MOEA method and the NSGA-II method.

Data of Nash-Harsanyi Leader-Follower Game Model and Multi-Objective Optimization Model
The stakeholders in the Huaihe River basin are the watershed management agency and the three previously mentioned provinces. According to the two water resources allocation models mentioned in Sections 2.2 and 2.3, the data required for the two models are total allocable water Q, ideal eco-environmental water demand s * , weight of the equity principle η, weight of efficiency principle (1 − η), water demand in each province r i , minimum survival water demand in each province λ i , minimum water allocation in each province I i , disagreement point in each province d i , bargaining weight in each province α i , objective function in each province u i (w i ), and integrated water utility in each province D i . The data used in the two models are all taken from Fu et al. (2018) [21]. Table 1 shows the data required for two transboundary water resources allocation models.

Water Resources Allocation Based on Asymmetric Nash-Harsanyi Leader-Follower Game Model and Its Solution in Huaihe River Basin
For the four stakeholders in the Huaihe River basin, the watershed management agency serves as leader and the three provinces serve as followers. The model of the watershed management agency is as follows: Subject to: Water balance constraints: Public water constraints: The public water s obtained by the above model of watershed management agency is 100.1 hundred million m 3 .
After the watershed management agency determines that the public water to be 100.1 hundred million m 3 , the model of the three provinces is as follows: Water balance constraints: The water allocation constraints in Henan province: The water allocation constraints in Anhui province: The water allocation constraints in Jiangsu province: We used the successive linear programming method to solve the above model for the three provinces. The amount of water allocated to Henan province, Anhui province, and Jiangsu province were 95.5 hundred million m 3 , 97.5 hundred million m 3 , and 107.4 hundred million m 3 , respectively, and the satisfaction of rates in the corresponding provinces were 68.6%, 61.6%, and 65.6%, respectively. In addition, the benefit of the three provinces were 734.6 hundred billion yuan, 600.8 hundred million yuan, and 683.3 hundred million yuan; the total benefit was 2018.7 hundred million yuan.

Water Resources Allocation Based on Multi-Objective Optimization Model and Its Solution in Huaihe River Basin
In the water resources allocation based on multi-objective optimization model, the eco-environmental benefit F 0 of the watershed management agency is given priority, so the public water allocated to watershed management agency is 100.1 hundred million m 3 . The objective function of the watershed management agency in transboundary water resources allocation can be converted into constraint conditions, then four objective functions converted into three objective functions. Water resources allocation based on the multi-objective optimization model is established as follows: The objective function of Henan province is: The objective function of Anhui province is: The objective function of Jiangsu province is: Subject to: Water balance constraints: The water allocation constraints in Henan province: The water allocation constraints in Anhui province: The water allocation constraints in Jiangsu province: where w 1 is the amount of water allocated to Henan province, w 2 is the amount of water allocated to Anhui province, and w 3 is the amount of water allocated to Jiangsu province. The NSGA-II method and the AR-MOEA method are used to solve the multi-objective optimization model established above. The population size of the two methods (NSGA-II method and AR-MOEA method) are both 100, and the maximum iterations of the two methods are both 1000. The multi-objective optimization model based on the NSGA-II method and the AR-MOEA method are both solved by using PlatEMO platform [26].
The Pareto front of multi-objective optimization based on NSGA-II method in the Huaihe River basin is shown in Figure 4. As can be seen from Figure 4, the points on the Pareto front solved by NSGA-II method are not uniformly distributed.    Figure 5 shows the Pareto front of multi-objective optimization based on the AR-MOEA method in the Huaihe River basin. In Figure 5, the points on the Pareto front solved by the AR-MOEA method are uniformly distributed.

Comparison of Water Resource Allocation Results Based on NSGA-II Method and AR-MOEA Method
In order to measure the effect of the NSGA-II method and the AR-MOEA method on multi-objective optimization in the three provinces of the Huaihe River basin, we introduced three categories of performance metrics (convergence metrics, diversity metrics, and the metrics considering convergence and diversity), as mentioned in Section 2.3.3. The average value and significance test of algorithm performance metrics of the two methods running for 30 times are shown in Table 2. there is no significant difference in the values of algorithm performance metrics of two methods; 2 indicates that significance test rejects the null hypothesis at the 5% significance level, this means that there is significant difference in the values of algorithm performance metrics of two methods. Table 2 shows: 1.
The convergence metric (GD) and composite metrics (IGD, HV) pass a 5% significance level; this means that there is no significant difference in convergence metric (GD) and composite metrics (IGD, HV). Therefore, we think that the two methods have almost the same effect on convergence metric (GD) and composite metrics (IGD, HV); 2.
The diversity metric (∆) cannot pass a 5% significance level; this means that there is a significant difference in the diversity metric (∆). As for the diversity metric (∆), when ∆ tends to zero, it indicates that the non-dominated solution set calculated by the multi-objective optimization algorithm has good diversity. The average value of the diversity metric (∆) in the AR-MOEA method is lower than that in the NSGA-II method. Therefore, the AR-MOEA method is better than the NSGA-II method in terms of the diversity metric (∆).
In conclusion, the AR-MOEA method is better than the NSGA-II method in terms of the diversity metric (∆).

Comparison of Water Resource Allocation Results Based on Multi-Objective Optimization Model and Game Model
Because the AR-MOEA method is superior in the diversity metric (∆) when solving the multi-objective optimization model of the Huaihe River basin water resources allocation, we choose the AR-MOEA method to solve the multi-objective optimization model of water resource allocation in the Huaihe River basin. Then, we compared the water resources allocation results based on the the asymmetric Nash-Harsanyi Leader-Follower game model and the multi-objective optimization model based on the AR-MOEA method. The benefits for the three provinces of Huaihe River obtained by the two model are shown in Figure 6.  Figure 6 shows that the water resources allocation result of the asymmetric Nash-Harsanyi Leader-Follower game model is on the Pareto front solved by the AR-MOEA method, that is, the solution of the asymmetric Nash-Harsanyi Leader-Follower game model is a non-dominated solution.
The multi-objective optimization solution is the closest solution to the ideal optimal Figure 6. The benefits of three provinces of the Huaihe River obtained by the two models. Figure 6 shows that the water resources allocation result of the asymmetric Nash-Harsanyi Leader-Follower game model is on the Pareto front solved by the AR-MOEA method, that is, the solution of the asymmetric Nash-Harsanyi Leader-Follower game model is a non-dominated solution.
The multi-objective optimization solution is the closest solution to the ideal optimal point, and the ideal optimal point is different under different preference information. The asymmetric Nash-Harsanyi Leader-Follower game model considers the principles of fairness and efficiency, so the solution nearest to the ideal optimal point under the principle of efficiency and the principle of fairness is chosen for the multi-objective optimization solution. The multi-objective optimization solution under the principle of efficiency can be converted into solving non-dominated solution of maximum comprehensive benefit. The non-dominated solution with the maximum comprehensive benefit can be obtained by the linear weighted sum method with equal non-negative weights. The multi-objective optimization solution under the principle of fairness can be converted into solving nondominated solution of equal water demand satisfaction rate. The non-dominated solution with an equal water demand satisfaction rate can be solved by the water balance equation and an equal water demand equation in the three provinces of the Huaihe River basin.
The water resources allocation results in the three provinces of the Huaihe River basin solved by the multi-objective optimization solution under different preference information (focusing on overall comprehensive benefit or overall fairness) and the asymmetric Nash-Harsanyi Leader-Follower game model (considering the principle of efficiency and fairness) were compared and analyzed, as shown in Table 3.  Table 3 shows: 1.
The water demand satisfaction rate of the three provinces of the Huaihe River basin is not balanced in the multi-objective optimization solution under the principle of efficiency, so the fairness of water resources allocation result in the multi-objective optimization solution under the principle of efficiency is ignored. The water demand satisfaction rate of Henan province is higher than that of Anhui province and Jiangsu province, which makes Anhui province and Jiangsu province unable to accept the water resources allocation scheme; 2.
The multi-objective optimization solution under the principle of fairness satisfies the fairness of water resources allocation result in the three provinces, but the total benefit of the multi-objective optimization solution under the principle of fairness is lower than the water resources allocation result solved by the asymmetric Nash-Harsanyi Leader-Follower game model. The above analysis makes the three provinces unable to accept the water resources allocation scheme; 3.
The water resources allocation result of the asymmetric Nash-Harsanyi Leader-Follower game model is a compromise solution between two kinds of multi-objective optimization solutions under different preference information (focusing on overall comprehensive benefit or overall fairness). The efficiency and fairness are considered comprehensively in the game model, which makes it easier for the three provinces to accept the water resources allocation scheme.
In order to explore the relationship between the water resources allocation results obtained by the linear weighted sum method with different objective function non-negative weights and the water resource allocation results based on the asymmetric Nash-Harsanyi Leader-Follower game model, the non-negative weights of objective function in linear weighted sum method are fitted by using a trial and error method.
When the non-negative weights of Henan province, Anhui province, and Jiangsu province are 0.281, 0.364, and 0.355, respectively, the water resources allocation obtained by the linear weighted sum method are equal to the water resources allocation results obtained by the asymmetric Nash-Harsanyi Leader-Follower game model. This result shows that the multi-objective optimization model can obtain the same water resources scheme as the Nash-Harsanyi Leader-Follower game model under a specific preference structure.
After the multi-objective optimization model obtains the Pareto front, it still needs to construct the preference information of Pareto front for the second time to make the optimal coordinated solution of the multi-objective decision. In this kind of integrated water resources allocation, each stakeholder cannot directly participate in water resources allocation, and the behavior influence of various stakeholders cannot be taken into account, which may make the obtained water resources unacceptable to various stakeholders. The asymmetric Nash-Harsanyi Leader-Follower game model can not only guarantee basic water demand, but can also better balance the economic development level among the three provinces. Moreover, in the asymmetric Nash-Harsanyi Leader-Follower game model, the three provinces of the Huaihe River basin can directly participate in the negotiation to make the water resources allocation results acceptable to all stakeholders, and the obtained result has a practical guiding role in water resources allocation.

Conclusions
In this study, we analyzed the difference between the asymmetric Nash-Harsanyi Leader-Follower game model and the multi-objective optimization model in water resources allocation. First, we established a transboundary water resources allocation model based on multi-objective optimization. Second, the AR-MOEA method was used to solve the Pareto front, and the Pareto front solved by AR-MOEA method was compared with the Pareto front solved by NSGA-II. Finally, we analyzedthe relationship between the water resources allocation result of the Nash-Harsanyi Leader-Follower game model and the Pareto front of the multi-objective optimization solution, and illustrated the difference between the asymmetric Nash-Harsanyi Leader-Follower game model and the multi-objective optimization model in water resources allocation. The asymmetric Nash-Harsanyi Leader-Follower game model and the multi-objective optimization model based on AR-MOEA method were applied to the Huaihe River basin. The main conclusions are as follows: 1.
The AR-MOEA method and the NSGA-II method were applied to solve the multiobjective optimization model of the Huaihe River basin water resources allocation. The results show that the AR-MOEA method is better than the NSGA-II method in terms of the diversity metric (∆); 2.
The solution of the asymmetric Nash-Harsanyi Leader-Follower game model is a non-dominated solution, and the asymmetric Nash-Harsanyi Leader-Follower game model can obtain the same water resources allocation scheme of the multi-objective optimal allocation model under a specific preference structure; 3.
After the multi-objective optimization model obtains the Pareto front, it still needs to construct the preference information of the Pareto front for the second time to make the optimal solution of multi-objective decision, while the game model can directly obtain the water allocation scheme in one time by participating in the negotiation.