A Fuzzy-Interval Dynamic Optimization Model for Regional Water Resources Allocation under Uncertainty

: In this study, a fuzzy-interval dynamic programming (FIDP) model is proposed for regional water management under uncertainty by combining fuzzy-interval linear programming (FILP) and dynamic programming (DP). This model can not only tackle uncertainties presented as intervals, but also consider the dynamic characteristics in the allocation process for water resources. Meanwhile, the overall satisfaction from users is considered in the objective function to solve the conﬂict caused by uneven distribution of resources. The FIDP model is then applied to the case study in terms of water resources allocation under uncertainty and dynamics for the City of Handan in Hebei Province, China. The obtained solutions can provide detailed allocation schemes and water shortage rates at different stages. The calculated comprehensive beneﬁts of economy, water users’ satisfaction and pollutant discharge (i.e., COD) are [2264.72, 2989.33] × 10 8 yuan, [87.50, 96.50] % and [1.23, 1.65] × 10 8 kg respectively with a plausibility degree (i.e., λ ± opt ) ranging within [0.985, 0.993]. Moreover, the beneﬁt from FIDP model under consideration of dynamic features is more speciﬁc and accurate than that of FILP model, whilst the water shortage rate from FIDP is [5.10, 9.10] % lower than that of FILP model.


Introduction
Due to population growth, economic development and consumption upgrade, global water consumption has increased by six times, and it has been continuing to grow steadily at an annual rate of about 1% during the past 100 years [1]. All of these would lead to the water shortage problem that is already pessimistically even severer, and seriously hinders the sustainable development of social economy. Managing water resources is an effective way to deal with the above challenges. However, in the process of management, experts and governors have encountered a lot of problems [2][3][4][5][6], such as dynamic variability and uncertainty, which are thorny and inevitable. Besides, in areas with water shortage, when the available water cannot meet the needs, unreasonable water allocation will lead to conflicts among users [7]. Therefore, it is definitely necessary to put forward a comprehensive model to deal with dynamic variabilities and uncertainties in water resources system as well as the contradiction between different users, so as to improve the management efficiency and the users' satisfaction.
The water resources system is of great complexities involving many uncertain factors, such as water use efficiency, water demand, pollutant discharge, water supply capacity and so on, and these uncertain factors could affect the structure for the optimal allocation model of water resources and resulting solutions [8][9][10][11]. Previously, scholars in related fields have got fruitful achievements in dealing with uncertainties in water resources management. For the optimization under uncertainties, mathematical methods that are commonly used Subject to: It is worth noting that: C ± g X ± = Max f ± g g = 1, 2, . . . , m (2a) 1×t , X ± ∈ R ± 4 t×1 , and R ± e means a set of interval numbers (e ∈ [1, 2,3,4]), g and h are core markers for maximizing and minimizing the objective functions individually, and i is the index of the constraints. f − , f + are the lower and upper bounds of f ± , and λ ± is the membership function in fuzzy decision-making. The larger the λ ± is, the more credible the calculation result would be; on the contrary, the smaller λ ± would lead to less credible results.

Dynamic Programming (DP)
The basic idea of dynamic programming is that it not only separates the current stage from the future stages, but also considers the current benefit and the future benefit together. Therefore, the optimal decision selection of each stage is from the overall consideration, which is generally different from the optimal choice of this stage [41]. Concretely, for a multi-stage decision-making problem, dynamic programming can divide it into several stages according to time or other characteristics, and each stage has several states and decision strategies [42]. The system transfers from one stage to the next according to a certain rule, and the purpose is to obtain the optimal strategy combining each stage [43]. The following Equation (3) is the state transition formula of dynamic programming, and it is also the most important part of dynamic programming.
Sustainability 2022, 14, 1096 4 of 20 where S j stands for the state variable at stage j, with l stages in total. x j−1 represents the decision variable at stage j − 1, and T(S j−1 , x j−1 ) is the state transition function [44].

Fuzzy-Interval Dynamic Programming (FIDP)
It is noted that multiple uncertainties and dynamic variability exist in the water resources system, which seriously affect effective planning and management of water resources. Although FILP and DP can efficiently address interval uncertainty, coordinate conflicts among different objective functions and characterize systems' dynamics individually, they are unable to deal with those problems at the same time. Therefore, this paper aims to propose a FIDP model by incorporating FILP and DP into one framework to comprehensively reflect both uncertainties and dynamic features in the water resources system. In addition, the function of the users' satisfaction is considered to solve the contradiction caused by uneven distribution of resources. The developed model is shown as follows.
Subject to: where the symbol G ± means the user's ideal demand for resources, and α is the weight coefficient of different users. And Equation (4c) can reflect the fairness for different users, Equation (4f) realizes the dynamic transition, and the state constraint after phase transition is achieved by Equation (4g). The steps of solving the FIDP model can be summarized as: (i) Establish FIDP model. (ii) Divide the model into two submodels through an improved two-step method [45]. In order to maximize λ ± , the upper bound submodel should be formulated firstly. (iii) Solve the upper bound submodel and obtain x + opt and λ + opt . (iv) Formulate the lower bound submodel for the FIDP model. (v) Solve the lower bound and obtain x − opt and λ − opt . (vi) According to the results of the above two models, the objective function values are calculated by formulate (2). (vii) Combining these two submodels, the optimal solution can be expressed as In general, the presented model can be applicable for the following problems: (i) For those problems with uncertain factors, this method can reflect them in model establishment, solution process and results in the form of interval numbers. (ii) For multi-stage decisionmaking problems, this model can provide specific schemes for every stage and global optimal solutions for the whole process. (iii) For multi-objective and multi-user problems, this model can coordinate the conflicts among different objective functions by maximizing satisfaction of the objective functions, and reduce the contradictions among users by considering the principle of fairness.

Overview of Handan City
The city of Handan is located in the southernmost part of Hebei Province, China, at the eastern foot of Taihang Mountain, bordering Shandong in the east, Henan in the south, Shanxi Province in the west and Xingtai City in the north. Its jurisdiction covers 6 districts, 1 county-level city and 11 counties. Its geographical location ranges 36 • 04 ~37 • 01 N and, 113 • 28 ~115 • 28 E with warm temperate semi-humid and semi-arid continental monsoon climate. The location of the area is shown in Figure 1.

Application of FIDP Model.
In order to primely solve the problems mentioned above, FIDP is applied to opt the allocation of water resources in Handan city. In detail, the established FIDP m would not only considers multiple objectives, such as the maximum economic benef maximum overall satisfaction of water users, and the maximum environmental be but also take the satisfaction of each water users into account. Meanwhile, the constr would refer to the water supply capacity, the minimum guaranteed water demand ideal water demand, the water delivery capacity, and the COD emission limit. In add the uncertain factors involved in this model (e.g., water use benefit coefficient, ideal w demand, minimum guaranteed water demand, weight coefficient, COD discharge c cient, maximum COD discharge, available water supply, water inflow at different st and water delivery capacity) can be expressed as interval parameters. Moreover, th namic factors in the process of water resources optimization, such as the water users' water demand, guaranteed water demand, available water supply and water alloc changing with the stage, would be reflected by dynamic programming. The frame gram of constructed FIDP model can be seen in Figure 2. In order to facilitate man to make decisions, each stage is divided equally by the planning year, in which, Jan March is the first stage, April-June is the second stage, the third stage is from July to tember, and the fourth stage is from October to December. Its formulation would b pressed in the following form: At present, water resource managers in Handan are facing with many water resource problems, such as water resource shortage, uneven distribution of precipitation, and serious water pollution. For example, according to the Water Resources Bulletin [46], the per capita water consumption of the city in 2019 is 2.02 × 10 2 m 3 with a population of 9.55 million. However, the water supply in 2019 is only 1927.84 × 10 6 m 3 , and the water shortage is 1.26 × 10 6 m 3 . In addition, 61.30%~76.50% of the annual precipitation falls between June and September, which is extremely inconsistent with the needs from various water users. Actually, each user's water demand, especially the agricultural water demand, is different with the season changes. The growing period of crops in Handan mainly ranges from March to August, with the largest water demand occurring at the second stage which would account for about 50% of the annual water consumption. It is noted that the development of agricultural cultivation is paid the most attention in Handan City, and its water consumption accounts for about 55% of the total water consumption. So how to provide periodic water allocation for each user is a challenge for managers. Moreover, due to the uncertainties existing in water supply and the temporal variations of the planning horizon, the water resources system also has a number of uncertain factors, such as the water inflows at different stages, water efficiency, water demand, and pollutant discharge, which should be fully considered. Therefore, how to allocate water resources reasonably to ensure the sustainable development of this region is an urgent problem for managers to solve under condition of discordant water supply and demand, as well as various uncertain factors.

Application of FIDP Model
In order to primely solve the problems mentioned above, FIDP is applied to optimize the allocation of water resources in Handan city. In detail, the established FIDP model would not only considers multiple objectives, such as the maximum economic benefit, the maximum overall satisfaction of water users, and the maximum environmental benefit, but also take the satisfaction of each water users into account. Meanwhile, the constraints would refer to the water supply capacity, the minimum guaranteed water demand, the ideal water demand, the water delivery capacity, and the COD emission limit. In addition, the uncertain factors involved in this model (e.g., water use benefit coefficient, ideal water demand, minimum guaranteed water demand, weight coefficient, COD discharge coefficient, maximum COD discharge, available water supply, water inflow at different stages, and water delivery capacity) can be expressed as interval parameters. Moreover, the dynamic factors in the process of water resources optimization, such as the water users' ideal water demand, guaranteed water demand, available water supply and water allocation changing with the stage, would be reflected by dynamic programming. The frame diagram of constructed FIDP model can be seen in Figure 2. In order to facilitate managers to make decisions, each stage is divided equally by the planning year, in which, January-March is the first stage, April-June is the second stage, the third stage is from July to September, and the fourth stage is from October to December. Its formulation would be expressed in the following form: Water transporting capacity constraint: The total amount of water used in each region would be limited by the water tra porting capacity in the region.
denotes the maximum capacity of water source i transporting to the regio in stage t (m 3 ).
The COD emission constraint: Due to serious damages of human activities to the ecological environment in rec years, more and more managers begin to pay attention to the impact of ecological en ronment with the development of economy. Accordingly, each region has formulated discharge capacity of pollutant COD to control environmental pollution. Therefore, optimal allocation of water resources should meet this requirement.

Objective Functions
The goal of this model is to maximize its membership function: Maxλ ± (5a)

Constraints
Economic benefit constraint (maximize the economic benefit): where t denotes stage of the planning year (t = 1, 2, 3, 4), i is water source (i = 1, 2, 3, 4, representing surface water, underground water, diverted water and recycled water), j means region (j = 1, 2, 3, . . . , 16,representing Urban,Wuan,Jize,Qiu,Quzhou,Guantao,She,Guangping,Chengan,Wei,Ci,Linzhang,Daming,Fengfeng,Yongnian and Feixiang), k stands for water user (k = 1, 2, 3, 4, representing agricultural, industrial, domestic and ecological), e ± tijk is net efficiency coefficient of water used by user k in region j in the t stage of planning year (yuan/m 3 ), a tijk denotes water relationship provided by water source i to user k in region j in stage t of planning year (water distribution is 1, unmatched water is 0), x ± tijk means water allocation from water source i to user k in region j in the stage t of planning year (m 3 ).
Social benefit constraint (maximize the overall satisfaction of water users): Considering the principle of fairness and justice, the weight coefficient α is introduced to balance the water satisfaction among water users and reduce the contradictions between water users and water supply departments.
where G ± tjk is the ideal water demand of user k in region j in stage t (m 3 ); α ± tjk means the weight coefficient of user k in the region j of the t stage.
Environmental constraint (minimize the chemical oxygen demand (COD) discharge of major pollutants in the region): While achieving the economic development, the pollution in the water utilization process should be comprehensively considered. The objective function should be established to measure the COD of the main pollutants in the region, so as to realize the balanced development of environment and economy.
where d ± tjk denotes the unit oxygen consumption generated by user k per unit water consumption in region j in stage t (kg/m 3 ).
Water supply constraint: In the tth stage, the sum of water supply from water source i to all water users is less than the maximum water supply of water source i.
where S ± ti stands for the maximum available water supply of water source i in stage t (m 3 ). Water demand constraint: The amount of water supplied to water users should be greater than or equal to the minimum guaranteed water demand of the user and less than or equal to the ideal water storage capacity of the user.
where D ± tjk means the minimum water demand of user k in region j in stage t (m 3 ). State transition equation: The maximum available water supply from different water sources in each stage is taken as the state variable, and the dynamic configuration of the model is realized through the water balance equation.
where C ± ti is the inflow of water source i in stage t (m 3 ). Water transporting capacity constraint: The total amount of water used in each region would be limited by the water transporting capacity in the region.
where Q ± tij denotes the maximum capacity of water source i transporting to the region j in stage t (m 3 ).
The COD emission constraint: Due to serious damages of human activities to the ecological environment in recent years, more and more managers begin to pay attention to the impact of ecological environment with the development of economy. Accordingly, each region has formulated the discharge capacity of pollutant COD to control environmental pollution. Therefore, the optimal allocation of water resources should meet this requirement.
where F ± j is the rated of COD emission in region j (kg). Nonnegative constrains: x ± tijk ≥ 0 (5j)

Data Collection and Analysis
This article takes Handan City as the research region and selects 2030 as the planning year. Due to the administrative adjustment of the city in recent years, this paper merged the Fuxing district, the Congtai district and the Hanshan district into the urban district to facilitate data compilation and calculation. The data needed in this model are related to economy, society, environment and water resources. All of these data are collected from related literature, field surveys, local statistical yearbooks and website information. Specifically, the water distribution relationship between water sources and users is obtained from the water resources bulletin [46]. The weight coefficient α is calculated based on the proportion of the added value of different users in each region in the recent two years' yearbooks [47]. The planned annual water transport capacity is obtained by combining the water conveyance capacity over the years and the pipeline network construction in recent years (https://www.h2o-china.com/news/295843.html, accessed on 14 June 2021). The unit oxygen consumption d and regional COD emissions are derived from related papers [48,49]. The benefit coefficient of agricultural water and industrial water is determined by the method of net output value allocation [48], and the benefit coefficient of domestic and ecological water use is obtained from relevant literatures [50,51]. According to the priority principle of domestic and ecological water use, the benefit coefficient was adjusted appropriately in this pater to rationalize the results, which are shown in Table 1. The available water supply of surface water, groundwater and diverted water in the planning year were predicted by the trend analysis method. The amount of recycled water was obtained according to the predicted regeneration rate of water consumption in the planning year. By comparing the predicted results with the water situations in recent years, it can be seen that there is similar water inflow situation in 2019. Thus, the water inflow situation of each stage in the planning year can be obtained based on the analysis of water supply proportion in 2019. 110% and 90% of the inflow were taken as the upper and lower bounds individually, and the results are shown in Figure 3.
It is necessary to calculate the planned annual water demand for optimal allocation of water resources. In this paper, the quota method was employed to forecast the water demand for agricultural, domestic and ecological use, whilst the equidimensional complementary residuals-residual modified GM (1, 1) model [52] was adopted to forecast the industrial water demand. Then, the water consumption situation in 2019 was analyzed to derive the water demand of every user at different stages of the planning year. Among them, the proportion of water demand at different stages of agriculture in the planning year is 15.80%, 49.80%, 23.40% and 11.00% respectively. The proportion of industrial water demand is 23.50%, 24.00%, 26.60% and 25.90%. The proportion of domestic water demand is 23.00%, 29.00%, 23.00% and 24.00%. The proportion of ecological water demand is 23.40%, 28.20%, 25.00% and 23.40%. In the planning year, 110% and 90% of the predicted water demand of different users in each region are taken as the upper and lower bounds of their water demand, respectively. The predicted results are shown in Table 2. It is necessary to calculate the planned annual water demand for optimal allocation of water resources. In this paper, the quota method was employed to forecast the water demand for agricultural, domestic and ecological use, whilst the equidimensional complementary residuals-residual modified GM (1, 1) model [52] was adopted to forecast the industrial water demand. Then, the water consumption situation in 2019 was analyzed to derive the water demand of every user at different stages of the planning year. Among them, the proportion of water demand at different stages of agriculture in the planning year is 15.80%, 49.80%, 23.40% and 11.00% respectively. The proportion of industrial water demand is 23.50%, 24.00%, 26.60% and 25.90%. The proportion of domestic water demand is 23.00%, 29.00%, 23.00% and 24.00%. The proportion of ecological water demand is 23.40%, 28.20%, 25.00% and 23.40%. In the planning year, 110% and 90% of the predicted water demand of different users in each region are taken as the upper and lower bounds of their water demand, respectively. The predicted results are shown in Table 2.

Results Analysis
In this study, the FIDP model suitable for Handan's water management was established to obtain the objective function values and water resources allocation schemes, which can be seen in Table 3. The λ + , λ − represent the maximum subordinate degree and the minimum subordinate degree respectively. In detail, by solving the model, the value of λ + opt is 0.993, the corresponding economic benefit is 2989.33 × 10 8 yuan, the satisfaction of users is 96.50%, and the social benefit is 1.23 × 10 8 kg. On the contrary, the value of λ − opt Sustainability 2022, 14, 1096 11 of 20 is 0.985, whilst the corresponding economic benefits, satisfaction and social benefits are 2264.72 × 10 8 yuan, 87.50% and 1.65 × 10 8 kg, respectively.  Table 4 shows the total amount of water allocated to different users in different regions of Handan City in the planning year, while Table 5 shows the total amount of water allocated from different water sources to different regions. It can be seen that the total amount of allocated water in Handan City in 2030 will be [175,412.60, 219,210.86] × 10 4 m 3 , and the total water shortage will reach [34,051.91, 36,800.32] × 10 4 m 3 according to the water demand forecasting results. As the minimum water demand in the planning year will reach [167,571.50, 204,809.62] × 10 4 m 3 , which is less than the allocated water, the water allocation in the planning year can meet its minimum guaranteed water demand on the whole.  (10 4 m 3 ).

Districts Agricultural Industrial Domestic Ecological
In 2030, the agricultural water consumption in Handan City will account for [52.90, 53.60] % of the total water distribution with the detailed allocation being [92,754.97, 117,454.65] × 10 4 m 3 . Since the agricultural water demand is affected by the season and climate, the water demand also changes at different stages of the planning year. In detail, the second stage is the main growth period of crops, and the agricultural water demand in this stage also increases correspondingly, accounting for [48.80, 49.80] % of the annual water demand. On the contrary, the amount of water distribution in the fourth stage accounts for the least proportion, which is only [11.10, 11.50] % of the total agricultural water distribution. The difference between these two stages is [35,988.55, 43,789.37] × 10 4 m 3 . The calculation results show that the satisfaction of the agricultural is [80.00, 82.90] % in 2030, and it reaches 80% in each stage, meeting its minimum water demand. Supported by the soil characteristics of each region, the leading agricultural industries in Urban, Weixian, Daming, and Yongnian have been developing rapidly, and the agricultural water consumption in these four regions would account for [54. 40, 56.40] % of the total agricultural water consumption in the city. The agricultural water distribution in the planning year is shown in Figure 4.  By solving the model, the industrial water distribution of Handan City in 2030 will be [36,210.94, 46,203.93] × 10 4 m 3 , accounting for [20.60, 21.10] % of the total water distribution. The satisfaction of industrial water consumption is [80.00, 83.50] % in the whole year, and such satisfaction is higher than 80.00% in each stage, meeting its minimum water demand. According to the solution results as shown in Figure 5, the third stage has the largest industrial water distribution, which is [9631.84, 12,290.27] × 10 4 m 3 , whereas the water distribution in the first stage is least with the allocation amount of [8509.83, 10,859.91] × 10 4 m 3 , and the difference between the two stages is [1122.11, 1432.34] × 10 4 m 3 . Among them, Urban, Wu 'an and Fengfeng are the major industrial water users, making a contribution of [78.90, 82.40] % for the whole city's industrial water consumption. By solving the model, the industrial water distribution of Handan City in 2030 will be [36,210.94, 46,203.93] × 10 4 m 3 , accounting for [20.60, 21.10] % of the total water distribution. The satisfaction of industrial water consumption is [80.00, 83.50] % in the whole year, and such satisfaction is higher than 80.00% in each stage, meeting its minimum water demand. According to the solution results as shown in Figure 5, the third stage has the largest industrial water distribution, which is [9631.84, 12,290.27] × 10 4 m 3 , whereas the water distribution in the first stage is least with the allocation amount of [8509.83, 10,859.91] × 10 4 m 3 , and the difference between the two stages is [1122.11, 1432.34] × 10 4 m 3 . Among them, Urban, Wu'an and Fengfeng are the major industrial water users, making a contribution of [78.90, 82.40] % for the whole city's industrial water consumption.
In 2030, the domestic water distribution in Handan City will be [34,907.31, 41,915.24] × 10 4 m 3 , accounting for [19.10, 19.90] % of the total water distribution. The satisfaction of domestic water use in the whole year is [93. 40,95.10] %, and it is higher than 92.00% in each stage indicating a high degree for guaranteed domestic water. With the change of temperature, the domestic water consumption at different stages also changes slightly to some extent. Specifically, the proportion of domestic water in the four stages of the planning year is [23.01, 23.44] %, [28.80, 29.00] %, [23.40, 24.10] % and [24.00, 24.40] % respectively. Obviously, the second stage consumes the most domestic water, whereas the first stage consumes the least proportion. During the planning year, the population in Urban and Yongnian will reach 3.40 × 10 6 , and the water allocated to these two areas will be [12,953.41, 15,082.32] × 10 4 m 3 correspondingly, accounting for [36.00, 37.10] % of the domestic water distribution to the whole city. The annual domestic water distribution in the planning year is shown in Figure 6. By solving the model, the industrial water distribution of Handan City in 2030 will be [36,210.94, 46,203.93] × 10 4 m 3 , accounting for [20.60, 21.10] % of the total water distribution. The satisfaction of industrial water consumption is [80.00, 83.50] % in the whole year, and such satisfaction is higher than 80.00% in each stage, meeting its minimum water demand. According to the solution results as shown in Figure 5, the third stage has the largest industrial water distribution, which is [9631.84, 12,290.27] × 10 4 m 3 , whereas the water distribution in the first stage is least with the allocation amount of [8509.83, 10,859.91] × 10 4 m 3 , and the difference between the two stages is [1122.11, 1432.34] × 10 4 m 3 . Among them, Urban, Wu 'an and Fengfeng are the major industrial water users, making a contribution of [78.90, 82.40] % for the whole city's industrial water consumption. In 2030, the domestic water distribution in Handan City will be [34,907.31, 41,915.24] × 10 4 m 3 , accounting for [19.10, 19.90] % of the total water distribution. The satisfaction of domestic water use in the whole year is [93.40, 95.10] %, and it is higher than 92.00% in each stage indicating a high degree for guaranteed domestic water. With the change of temperature, the domestic water consumption at different stages also changes slightly to some extent. Specifically, the proportion of domestic water in the four stages of the planning year is [23.01, 23.44] %, [28.80, 29.00] %, [23.40, 24.10] % and [24.00, 24.40] % respectively. Obviously, the second stage consumes the most domestic water, whereas the first stage consumes the least proportion. During the planning year, the population in Urban and Yongnian will reach 3.40 × 10 6 , and the water allocated to these two areas will be [12,953.41, 15,082.32] × 10 4 m 3 correspondingly, accounting for [36.00, 37.10] % of the domestic water distribution to the whole city. The annual domestic water distribution in the planning year is shown in Figure 6. To a certain extent, the development degree of ecological environment reflects the sustainability level of the region. In 2030, the ecological water distribution in Handan will be [11,539.41, 13,638.02] × 10 4 m 3 , with a contribution of [6.20, 6.60] % for the total water To a certain extent, the development degree of ecological environment reflects the sustainability level of the region. In 2030, the ecological water distribution in Handan will be [11,539.41, 13,638.02] × 10 4 m 3 , with a contribution of [6.20, 6.60] % for the total water distribution. The satisfaction of ecological water use is [96.70, 100.00] % in the whole year with the satisfaction degree over 93.00% in each stage, which reflects the priority in ecological development. It can be seen from Figure 7 that the distribution of ecological water reaches the annual maximum amount of [3254.21, 3704.50] × 10 4 m 3 in the second stage, which is [404.34, 554.02] × 10 4 m 3 more than the least water distribution in the first stage. The ecological water consumption in the third and fourth stage is [2884.81, 3333.10] × 10 4 m 3 and [2700.24, 3300.23] × 10 4 m 3 respectively, accounting for [24.40, 25.00] % and [23.40, 24.20] % of the ecological water consumption in the whole year. Among them, the ecological water consumption in the urban area and Daming county is relatively huge, contribution [32.10, 34.30] % and [23.10, 23.90] % to the total ecological water consumption respectively. This indicates that these two regions pay close attention to ecological environment construction. Based on the analysis above, it can be known that FIDP model can provide global optimal solutions for the planned annual water distribution scheme, as well as specific water distribution schemes at different stages of the year under dynamics and uncertainties. In the planning year, the second stage has the largest water distribution of [68,246.70, 84,190.91] × 10 4 m 3 , while the fourth stage has the smallest water distribution of [30,648.40, 39,012.42] × 10 4 m 3 . The water distribution difference between the two stages is [37,598.12, 45,178.58] × 10 4 m 3 , and the difference accounts for [20.60, 21.40] % of the annual water distribution.

Model Comparison
In order to verify the effectiveness of the proposed FIDP model, the application of FILP model to the case study is provided for comparison, which is shown in Appendix A.
The difference between FILP model and FIDP model is that it deletes dynamic programming and parameter t, but their objective functions, constraints, decision variables and solution methods remain the same. Because the FIDP model takes into account the dynamic factors in different stages of water resources system, the solution results cannot only conform to the case study, but also achieve global optimization under the local optimal conditions of each stage. However, the FILP model only aims at optimality over the whole planning year without considering the dynamic variability of regional water resources system, which would imply that the water use efficiency, water consumption, wa- Based on the analysis above, it can be known that FIDP model can provide global optimal solutions for the planned annual water distribution scheme, as well as specific water distribution schemes at different stages of the year under dynamics and uncertainties. In the planning year, the second stage has the largest water distribution of [68,246.70, 84,190.91] × 10 4 m 3 , while the fourth stage has the smallest water distribution of [30,648.40, 39,012.42] × 10 4 m 3 . The water distribution difference between the two stages is [37,598.12, 45,178.58] × 10 4 m 3 , and the difference accounts for [20.60, 21.40] % of the annual water distribution.

Model Comparison
In order to verify the effectiveness of the proposed FIDP model, the application of FILP model to the case study is provided for comparison, which is shown in Appendix A.
The difference between FILP model and FIDP model is that it deletes dynamic programming and parameter t, but their objective functions, constraints, decision variables and solution methods remain the same. Because the FIDP model takes into account the dynamic factors in different stages of water resources system, the solution results cannot only conform to the case study, but also achieve global optimization under the local optimal conditions of each stage. However, the FILP model only aims at optimality over the whole planning year without considering the dynamic variability of regional water resources system, which would imply that the water use efficiency, water consumption, water demand, water supply and other factors in the FILP model cannot be adjusted correspondingly with the seasonal changes. Therefore, the FIDP model has better optimal solutions and stronger applicability than FILP. The detailed analysis is as follows.
Based on the FILP model, the membership function λ ± = [0.952, 0.992], f ± 1 = [2171. 42, 3124.16] × 10 8 yuan, f ± 2 = [75.00, 84.00] %, f ± 3 = [1.17, 1.85] × 10 8 kg, and the comparison result of two models is shown in Figure 8. It can be seen that, compared with the FILP model, the ranges of the solution results of the FIDP model are reduced in different degrees, making the results more specific and accurate. In detail, the λ ± of the FIDP model is not only reduced in scope, but also presents increases in its overall value, with its upper bound and lower bounds increased by 0.10% and 3.30% respectively. For f ± 1 and f ± 3 , not only the ranges of their value are reduced by 23.90% and 38.20%, but also their lower bounds are increased by 93.30 × 10 8 yuan and 0.06 × 10 8 kg respectively, which are more accurate. The value of f ± 2 in FIDP model is improved by 12.50%, which will alleviate the conflicts between local government and users, and among different users more effectively. Consequently, it can be concluded that FIDP model proposed in this paper performs better and is more suitable for the optimization of water resources allocation in this area. and 3 f ± , not only the ranges of their value are reduced by 23.90% and 38.20%, but also their lower bounds are increased by 93.30 × 10 8 yuan and 0.06 × 10 8 kg respectively, which are more accurate. The value of 2 f ± in FIDP model is improved by 12.50%, which will alleviate the conflicts between local government and users, and among different users more effectively. Consequently, it can be concluded that FIDP model proposed in this paper performs better and is more suitable for the optimization of water resources allocation in this area. The water distribution scheme of FILP model is shown in Table 6, and the water shortage rate of the two models is compared in Figure 9. As presented in Table 6, the total water distribution of the FILP model is [167,643.11, 205,259.94] × 10 4 m 3 , which is reduced by [7769.43, 13,951.80] × 10 4 m 3 compared with the FIDP model. As the water distribution decreases, the corresponding water shortage rate will be [19.70, 20.00] %, with an increase of [5.10, 9.10] % compared with [10.60, 14.90] % from the FIDP model. Compared with FILP, the water shortage rate from FIDP has declined in all regions, especially in Qiu, Quzhou, Guangping, Ci and Feixiang, with a decrease of [7.90, 16.70] %, [7.00, 15.00] %, [8.30, 20.00] %, [9.90, 20.00] %, and [9.80, 14.30] % respectively. Therefore, the model has good applicability to water resources allocation in water-scarce areas.  The water distribution scheme of FILP model is shown in Table 6, and the water shortage rate of the two models is compared in Figure 9. As presented in Table 6, the total water distribution of the FILP model is [167,643.11, 205,259.94] × 10 4 m 3 , which is reduced by [7769.43, 13,951.80] × 10 4 m 3 compared with the FIDP model. As the water distribution decreases, the corresponding water shortage rate will be [19.70, 20.00] %, with an increase of [5.10, 9.10] % compared with [10.60, 14.90] % from the FIDP model. Compared with FILP, the water shortage rate from FIDP has declined in all regions, especially in Qiu, Quzhou, Guangping, Ci and Feixiang, with a decrease of [7.90, 16.70] %, [7.00, 15.00] %, [8.30, 20.00] %, [9.90, 20.00] %, and [9.80, 14.30] % respectively. Therefore, the model has good applicability to water resources allocation in water-scarce areas.

Agricultural Industrial Domestic Ecological
In summary, the FIDP model has the following advantages over the FILP model: (i) The obtained target value intervals are more specific and accurate. (ii) This model can improve the overall satisfaction of the water users and alleviate the water contradiction among them. (iii) The water shortage rate of FIDP model is lower than that of FILP, which effectively alleviates the contradiction between water supply and demand. (iv) Last but not least, FIDP model can get the allocation schemes of each stage in the planning year, and provide theoretical basis for water distribution decision-making in more detail way. Therefore, the model has a good performance in dealing with the dynamic changes of water resources system, and has advantages in optimizing the target value and reducing the water shortage rate.

Conclusions
This study combines fuzzy-interval linear programming and dynamic programming to establish a fuzzy-interval dynamic programming (FIDP) method. The system uncertainty is expressed in the form of interval numbers in the model establishment and results presentation. In order to solve the conflicts among users caused by uneven distribution of In summary, the FIDP model has the following advantages over the FILP model: (i) The obtained target value intervals are more specific and accurate. (ii) This model can improve the overall satisfaction of the water users and alleviate the water contradiction among them. (iii) The water shortage rate of FIDP model is lower than that of FILP, which effectively alleviates the contradiction between water supply and demand. (iv) Last but not least, FIDP model can get the allocation schemes of each stage in the planning year, and provide theoretical basis for water distribution decision-making in more detail way. Therefore, the model has a good performance in dealing with the dynamic changes of water resources system, and has advantages in optimizing the target value and reducing the water shortage rate.

Conclusions
This study combines fuzzy-interval linear programming and dynamic programming to establish a fuzzy-interval dynamic programming (FIDP) method. The system uncertainty is expressed in the form of interval numbers in the model establishment and results presentation. In order to solve the conflicts among users caused by uneven distribution of resources, the FIDP model incorporates the overall satisfaction of users into the objective function in order to reflect fairness in the solution results. As for the dynamic variability, the FIDP model introduces dynamic programming to obtain the specific allocation schemes at different stages of the planning year. Therefore, the model is able to handle the problem of uneven resources allocation under uncertainty and dynamics.
In this paper, the water resources allocation under uncertainty and dynamics in Handan City is studied to verify the feasibility of the model. The solutions cannot only get the optimized target value in this area, but also get the specific water allocations and water shortages for each water user at different stages in the form of intervals, so that the government and users can adjust their strategies to deal with the crisis in time. In addition to the economic and environmental goals, the obtained solution also provides the satisfaction of water users by introducing a weight coefficient, which effectively alleviates the contradictions among users. Through model comparison, the FIDP is proved to be superior to FILP model in this area. Consequently, the developed FIDP model would be more rational and applicable for regional water allocation under uncertainty and dynamics, which is of great help to the sustainable development of the region.