Next Article in Journal
Impact Analysis of Slope on the Head Loss of Gas-Liquid Two-Phase Flow in Siphon Pipe
Previous Article in Journal
Optimization of Horseradish Peroxidase Catalytic Degradation for 2-Methyl-6-Ethylaniline Removal Using Response Surface Methodology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Optimization-Evaluation Agricultural Water Planning Approach Based on Interval Linear Fractional Bi-Level Programming and IAHP-TOPSIS

1
Centre for Agricultural Water Research in China, China Agricultural University, Beijing 10083, China
2
Wuwei Experimental Station for Efficient Water Use in Agriculture, Ministry of Agriculture and Rural Affairs, Wuwei 733000, China
3
Faculty of Graduate Studies and Research, University of Regina, Regina, SK S4S 0A2, Canada
*
Author to whom correspondence should be addressed.
Water 2019, 11(5), 1094; https://doi.org/10.3390/w11051094
Submission received: 5 April 2019 / Revised: 3 May 2019 / Accepted: 20 May 2019 / Published: 24 May 2019
(This article belongs to the Section Water Resources Management, Policy and Governance)

Abstract

:
In this study, an interval linear fractional bi-level programming (ILFBP) model was developed for managing irrigation-water resources under uncertainty. The ILFBP can fully address system fairness, uncertainties, and the leader–follower relationship of decision makers in the optimization process, which can better reflect the complexity of real decision-making process and help formulate reasonable water policies. An interactive fuzzy coordination algorithm based on satisfaction degree was introduced to solve the ILFBP model. In order to evaluate the applicability of optimization schemes, the interval analytic hierarchy process (IAHP) and the interval technique for order preference by similarity to an ideal solution (TOPSIS) method were integrated as IAHP-TOPSIS. To verify its validity, the developed optimization-evaluation framework was applied to an irrigation water management case study in the middle reaches of the Shiyang River Basin, located in the northwest China. The ILFBP model results show that the total water allocation is [6.73, 7.37] × 108 m3, saving nearly 0.9 × 108 m3 more than the current situation. The benefit per unit of water is [2.38, 2.95] yuan/m3, nearly 0.4 yuan/m3 more than the status quo, and the Gini coefficient is within a reasonable range of [0.35, 0.38]. The ILFBP model can well balance economic benefits and system fairness. Through the evaluation bases on IAHP-TOPSIS, the results of ILFBP show better water allocation effects and applicability than the other two models in this study area. Furthermore, due to various characteristics such as geographical location, population and area, there are three irrigation districts, Xiying, Donghe, and Qinghe, showing higher importance than others when considering regional water allocation. These findings can provide useful information for limited water resource managers and help decision makers determine effective alternatives of water resource planning under uncertainty.

1. Introduction

Water is an indispensable and essential resource for human survival and development. With the development of social economy, the conflicts between supplies and demands of water resources have become increasingly prominent [1]. As the large freshwater consumer, agriculture accounts for more than 70% of the total freshwater consumption in the world [2,3]. Especially in the semi-arid and arid regions, water deficits limit sustainable development of the local agriculture, and thus more efforts should be made to improve the usage efficiency of limited irrigation water [4].
The agricultural water management systems are complicated with multiple uncertainties, including inexactness, fuzziness, and randomness existed in precipitation and river runoff, crop market prices variations, and irrigation target adjustments. The necessity of dealing with these uncertainties brings a lot of difficulties in the planning processes of irrigation water management [5,6]. To address these uncertainties in practical problems, many methods have been proposed for helping decision makers formulate better water policies [7,8,9,10,11]. Generally, these methods can be categorized into fuzzy programming (FP), stochastic programming (SP), and interval programming (IP). Among these categories, IP is the most widely used planning method under uncertainty due to its low data requirements [12,13]. Although IP can deal with the uncertainties expressed as interval numbers, it cannot fully reflect multiple decision-making hierarchy in the practical problems. To obtain the satisfied irrigation water allocation schemes, the various requirements from decision makers of different levels should be taken into account within one framework when formulating allocation models [14]. To reflect the tradeoff between two decision-making levels with different concerns, the bi-level programing (BP) approach with a two-level structure was developed for supporting regional water policies formulation [15]. A leader–follower decision-making strategy was incorporated into the optimization process of BP to generate satisfied decision-making plan for both the upper-level and the lower-level decision makers. However, more challenges exist when attempting to tackle interval uncertainty and coordinate the preferences of different-level stakeholders simultaneously.
In recent years, some researches were conducted to integrate interval linear programming (ILP) and BP into one framework to address such problems in irrigation water management [16]. An interval linear bi-level programming (ILBP) approach was thus developed and applied to real-world problems of water resources planning to address different-level decision-making preferences and uncertainties expressed as intervals [17]. These researches proved that ILBP would be an effective tool for decision makers to formulate water policies. However, the concerns of decision makers in these researches mainly focused on the economic and environmental factors such as benefits and water volumes water resource allocation systems. It was neglected that water resources are not only accounted as economic resources but also represent a kind of social resource. Unfair allocations among water users sometimes may lead to serious social conflicts. Therefore, the fairness issue needs to be considered in the water resource allocation process [18,19]. Many indicators have been proposed for measuring system fairness, such as the Gini coefficient, weighted proportional fairness, and so on [20,21,22]. Among them, the Gini coefficient is the most widely used measurement of inequality [23,24]. The Gini coefficient was proposed by the Italian economist Corrado Gini in 1912 to measure the balance of income in a country or region. It is usually calculated via the Lorenz curve [25], and has been applied to the fairness studies in many fields of education, economy, and resources management [22,26,27,28]. Therefore, the Gini coefficient is considered as a proper configuration goal of water resources planning to guarantee the fairness of water allocation. However, because the Gini coefficient is always presented as fractional, it is necessary to introduce the interval linear fractional programming (ILFP) model into the ILBP framework to solve such a problem. This results in the development of an interval linear fractional bi-level programming (ILFBP) model, which can simultaneously tackle interval uncertainty, coordinate the preferences of different-level decision makers and reflect system fairness in one framework to obtain satisfied water allocation schemes.
Moreover, although the mathematical optimization models can generate water allocation schemes under uncertainty, few of them are incorporated with the results evaluation methods to measure their applicability and effectiveness in specific cases. This limitation may reduce the guiding significance of these optimization planning results in real-world practice [29,30,31]. Therefore, the evaluation methods based on actual regional conditions are necessary for measuring impacts of implementing the optimization results in various agricultural regions [32]. However, compared to general evaluation methods proposed for assessing the deterministic optimization results, the evaluation of optimization results under interval uncertainty would be more difficult. In order to address such a challenge, an interval analytic hierarchy process (IAHP) method has been proposed to address system uncertainties associated with systematic errors, incomplete information and fuzzy subjective judgment, and shown good performance in evaluating multiple schemes under interval uncertainty [33,34]. Generally, the calculation of IAHP often involves the ordering of interval numbers and the typical method is to directly use the mean of the interval numbers (i.e., the average of the upper and lower bounds) [35]. The limitation of this approach is that it ignores the effects of interval lengths and might lead to unreasonable comparison results. Besides, the possibility matrix has been used in many researches as a possible way to sort interval numbers [36,37], but it becomes complicated and inefficient when dealing with more than two indicators. More effective methods are thus needed to sort interval numbers in the results evaluation process. As the interval-parameter TOPSIS (ITOPSIS) could sort interval numbers via calculating the distances of alternative schemes to the ideal scheme, the integrated method of IAHP and ITOPSIS (IAHP-TOPSIS) is considered as an effective way to deal with interval numbers’ sorting in the evaluation of optimization results [38,39,40]. Through the IAHP-TOPSIS method, the effects evaluation of the optimization schemes can be carried out and would further provide decision-making support for local water policies development. There have been few studies that can address both optimization and evaluation processes within one framework and simultaneously consider interval characteristics associated with agriculture water management systems.
The main objective of this study is to develop an optimization-evaluation approach for agricultural water planning based on the ILFBP and IAHP-TOPSIS methods to help decision makers formulate reasonable regional water policies. The optimization model (ILFBP) proposed for optimal irrigation water allocation and the evaluation model (IAHP-TOPSIS) used for evaluating the effectiveness of optimization schemes will be integrated within one framework under interval uncertainty to address the complexities of practical problems. Such an approach can: (1) Tackle interval numbers existed in both optimization and evaluation processes; (2) consider the fairness issue in water allocation; (3) make trade-offs between different concerns of two-level decision makers; (4) generate multiple water resource allocation schemes under different runoff levels; and (5) provide the corresponding evaluation results of the obtained allocation schemes. The proposed approach will also be applied to a real-world case study to demonstrate its applicability and effectiveness, wherein the regional water managers need to identify more effective irrigation water allocation policies to support the sustainable development of local agriculture.

2. Methodology

This section presents the development of an optimization-evaluation framework for an agricultural water allocation system, which entails the following three major components: (1) The parameter preprocessing and development of interval linear fractional bi-level programming (ILFBP); (2) the ILFBP model for irrigation water resources optimal allocation; and (3) the IAHP-TOPSIS method for evaluating the effectiveness of the optimization schemes. These components and their connections are shown in Figure 1 in detail.

2.1. Interval Linear Fractional Bi-Level Programming (ILFBP)

The general ILFBP problem can be described follows [41,42]:
(Upper-level objective)
min f 1 ± ( x 1 ± ,    x 2 ± )    =    c 11 ± x 1 ± + c 12 ± x 2 ± + α 1 ± d 11 ± x 1 ± + d 12 ± x 2 ± + β 1 ±
(Lower-level objective)
min f 2 ± ( x 1 ± ,    x 2 ± )    =    c 21 ± x 1 ± + c 22 ± x 2 ± + α 2 ± d 21 ± x 1 ± + d 22 ± x 2 ± + β 2 ±
Subject to:
G ±    =    { ( x 1 ± ,    x 2 ± ) | g i ± ( x 1 ± ,    x 2 ± )       0 ,      i    =    1 ,    2 ,    ,    m ,      x 1 ± ,    x 2 ±       0 }
where the upper-level variables are x 1 R n 1 and the lower-level variables are x 2 R n 2 ; c 11 ,    c 21 ,    c 12 ,    c 22 ,    d 11 ,    d 21 ,    d 12 ,    d 22 ,    α 1 ,    β 1 ,    α 2 ,    β 2 are interval parameters; G is the set of bi-level constraints; f 1 and f 2 are the upper-level and lower-level objective functions, respectively; the upper-level decision maker (ULDM) and the lower-level decision maker (LLDM) accept the leader–follower Stackelberg game [43].
To solve model (1), an interactive fuzzy coordination algorithm is proposed through analyzing the interrelationship between the objective functions and the constraints [16,41]. The upper- and lower-bound submodels of the ILFBP can be decomposed [44]. In particular, when the number of interval parameters in the ILFBP model is relatively small, the upper- or lower-bound parameters corresponding to two ILFBP submodels can be determined respectively by comparing all the combinations of the intervals. If the objective functions of the ULDM and LLDM problems are formulated as the maximization goals, they need to be converted as the equivalent minimization goals, i.e., max f    =    min ( f ) , and then the general ILFBP model can be solved according to following steps.
(1)
The ILFBP submodel corresponding to the lower-bound objectives can be formulated as:
min f 1 ( x 1 , x 2 )
min f 2 ( x 1 , x 2 )
Subject to:
G =    { ( x 1 ,    x 2 ) | g i ( x 1 ,    x 2 )       0 ,        i    =    1 ,    2 ,    ,    m ,        x 1 ,    x 2       0 }
(2)
The ULDM problem and the LLDM problem can be solved independently with the solutions of ( x 1 U ,    x 2 U ,    f 1 U ) and ( x 1 L ,    x 2 L ,    f 1 L ) . The ranges of the decision variables ( x 1 ) should be around within the maximum tolerances ( r 1 ) of their most preferred values. The membership functions that specify the satisfaction degrees of x 1 can be presented as follows:
μ x 1 ( x 1 )    =    { x 1 ( x 1 U r 1 ) r 1 , x 1 U r 1 < x 1 < x 1 U ( x 1 U + r 1 ) x 1 r 1 , x 1 U < x 1 < x 1 U + r 1 0 , o t h e r w i s e
where x 1 U are the most preferred values of upper-level decision variables while ( x 1 U + r 1 ) and ( x 1 U r 1 ) are the least preferred values of upper-level decision variables; μ x 1 ( x 1 ) are the satisfaction degrees of the upper-level decision variables x 1 , which increase within the intervals of [ x 1 U r 1 , x 1 U ] and decrease within the intervals of [ x 1 U , x 1 U +    r 1 ] , even though the other decision variable values are not acceptable.
(3)
For the ULDM problem, the objective function can be considered acceptable when f 1 f 1 U and unacceptable when f 1 f 1 = f 1 ( x 1 L , x 2 L ) . The LLDM problem can be solved with the optimal solutions of ( x 1 L , x 2 L ) , which in turn provides the ULDM problem with the objective value of f 1 by considering f 1 f 1 unacceptable in practice. The membership function that specifies the satisfaction degree of the ULDM objective value can be stated as:
μ f 1 [ f 1 ( x ) ] = { 1 , f 1 ( x ) < f 1 U f 1 f 1 ( x ) f 1 f 1 U , f 1 U < f 1 ( x ) < f 1 0 , f 1 ( x ) > f 1
In contrast, the LLDM objective value has the following membership function:
μ f 2 [ f 2 ( x ) ] = { 1 , f 2 ( x ) < f 2 L f 2 f 2 ( x ) f 2 f 2 L , f 2 L < f 2 ( x ) < f 2 0 , f 2 ( x ) > f 2
where f 2    =    f 2 ( x 1 U , x 1 U ) and the preference decreases within the interval of [ f 2 L ,    f 2 ] .
(4)
Consequently, the solution of submodel (2) with maximum overall satisfaction for both the ULDM and LLDM can be obtained through the following model:
max λ
Subject to:
x G
μ x 1 ( x 1 ) λ I
μ f 1 [ f 1 ( x ) ] λ
μ f 2 [ f 2 ( x ) ] λ
x 1 , x 2 0
λ [ 0 , 1 ]
where λ is the overall satisfactory degree of the lower-bound ILFBP submodel; I is a column vector with the same dimension as μ x 1 or x 1 and all elements of I are equal to 1. Thus, the solutions of the lower-bound ILFBP submodel and the corresponding objective values (i.e., x 1 o p t , x 2 o p t , f 1 o p t ,   and   f 2 o p t ) can be obtained through model (6).
(5)
Accordingly, the upper-bound ILFBP submodel can be formulated as follows:
(Upper-level objective)
min f 1 + ( x 1 + , x 2 + )
(Lower-level objective)
min f 2 + ( x 1 + , x 2 + )
Subject to:
G +   =   { ( x 1 + , x 2 + ) | g i + ( x 1 + ,    x 2 + )       0 ,          i    =    1 ,    2 ,    ,    m ,          x 1 + , x 2 + 0 }
x 1 + x 1 opt
x 2 + x 2 opt
(6)
Similar to the solution process of submodel (2), the solution of submodel (7) with maximum overall satisfaction for both the ULDM and LLDM can be obtained through the following model:
max λ +
Subject to:
x + G +
μ x 1 + ( x 1 + ) λ + I
μ f 1 + [ f 1 + ( x ) ] λ +
μ f 2 + [ f 2 + ( x ) ] λ +
x 1 + x 1 opt
x 2 + x 2 opt
λ + [ 0 , 1 ]
where:
μ x 1 + ( x 1 + ) = { x 1 + ( x 1 U + r + ) r 1 + , x 1 U + r 1 + < x 1 + < x 1 U + ( x 1 U + + r 1 + ) x 1 + r 1 , x 1 U + < x 1 + < x 1 U + + r 1 + 0 , o t h e r w i s e
μ f 1 + [ f 1 + ( x ) ] = { 1 , f 1 + ( x ) < f 1 U + f 1 + f 1 + ( x ) f 1 + f 1 U + , f 1 U + < f 1 + ( x ) < f 1 + 0 , f 1 + ( x ) > f 1 +
μ f 2 + [ f 2 + ( x ) ] = { 1 , f 2 + ( x ) < f 2 L + f 2 + f 2 + ( x ) f 2 + f 2 L + , f 2 L + < f 2 + ( x ) < f 2 + 0 , f 2 + ( x ) > f 2 +
The solutions of the upper-bound ILFBP submodel and the objective values (i.e., x 1 o p t + , x 2 o p t + , f 1 o p t + ,   and f 2 o p t + ) can be obtained from model (8).
(7)
Finally, the solutions for the ILFBP model can be generated as follows:
x 1 o p t ± = [ x 1 o p t , x 1 o p t + ]
x 2 o p t ± = [ x 2 o p t , x 2 o p t + ]
f 1 o p t ± = [ f 1 o p t , f 1 o p t + ]
f 2 o p t ± = [ f 2 o p t , f 2 o p t + ]

2.2. The ILFBP Model for Agricultural Water Resources Management

To achieve effective management of limited irrigation water resources, the decision makers often desire optimal allocation schemes to avoid high costs and large water consumptions, especially for the arid and semi-arid areas. Moreover, the decision makers need to consider the system fairness, balance the different-level preferences, and address interval uncertainties in practical problems [45,46,47]. To deal with these challenges of agricultural water resources management, an ILFBP model can be formulated as follows:
(Upper-level objective)
max f 1 ± = i = 1 I j = 1 J h = 1 H p h ( B i j ± Y i j ± ( W i j h ± k i + P i j h ± ) J L i W i j h ± ) i = 1 I J B i A i
where f 1 is the upper-level objective of maximizing the system’s net economic benefits (yuan); i is the index for irrigation districts (i = 1, 2, …, 10); j denotes the type of crops (j = 1, 2, 3); h denotes the levels of water availability (h = 1, 2, 3, 4, 5); W i j h are the decision variables, representing the amount of irrigation water allocated to the crop j in the irrigation district i under the flow level h (m3); p h are the possibilities of different water levels; B i j are the prices of crop j in the irrigation district i (yuan/kg); Y i j are water productivities of crop j in the irrigation district i (kg/m3); A i are irrigated areas of subarea i (ha); k i are the irrigation water use efficiencies; P i j h are effective precipitations (m3); J B i are basic water prices (yuan/ha); and J L i are the metering water prices (yuan/m3).
(Lower-level objective)
min f 2 ± = l = 1 I k = 1 I | j = 1 J h = 1 H p h ( W i j h ± + P i j h ± ) R l j = 1 J h = 1 H p h ( W k j h ± + P k j h ± ) R k | 2 I i = 1 I j = 1 J h = 1 H p h ( W i j h ± + P i j h ± ) R i
where f 2 is the lower-level objective function to reflect the fairness of water allocation among irrigation districts via the Gini coefficient; and R i are populations of the irrigation district i.
The Gini coefficient varies within the range of [0, 1]. The smaller Gini coefficient represents the higher fairness of water allocation among irrigation districts. Conversely, the larger Gini coefficient represents the lower fairness of water allocation. In regard to income distribution levels, a Gini coefficient less than 0.2 stands for absolute evenness, that between 0.2 and 0.3 represents relative evenness, that between 0.3 and 0.4 indicates basically reasonable distribution, and that between 0.4 and 0.5 indicates unreasonable distribution. When a Gini coefficient reaches 0.5 or more, it means a significant gap in income levels. Generally, a Gini coefficient higher than 0.4 is used as a sign of unfairness.
Subject to:
(1)
Water availability constraint
i = 1 I j = 1 J W i j h ±       T W h ±        h
j = 1 J W i j h ± T W i h ±     i ,    h
where T W h are the total water supplies (m3) under the flow level h; and T W i h denote the water supplies (m3) for the irrigation district i under flow level h.
(2)
Irrigation water demand constraint
I R i j h , min ±    W i j h ± k i + P i j h ± I R i j h , max ±      i ,    j ,    h
where I R i j h , min ± and I R i j h , max ± are respectively the minimum and maximum crop water requirement (m3) of crop j in the irrigation district i under flow level h.
(3)
Non-negative constraint
W i j h ± 0      i ,    j ,    h

2.3. The IAHP-TOPSIS Evaluation Approach Coupling the IAHP and ITOPSIS Methods

The IAHP method is widely used for multi-criteria decision making (MCDM) under uncertainty [48]. It can build an evaluation framework based on the analysis of interval parameters, and obtain the interval weights of the indicators. However, more effective methods are needed to sort interval numbers when using the IAHP. The interval-parameter TOPSIS method (ITOPSIS) can analyze the uncertain information of the scheme results and sort interval numbers through measuring the distances between solution alternatives and the ideal solution. Therefore, for improving the assessment accuracy under interval uncertainty, the IAHP and ITOPSIS methods are integrated to evaluate the optimal allocation results of water resources. Particularly, the IAHP method is used for calculating the indicators’ interval weights and the ITOPSIS method is used for ranking the uncertain water resources allocation schemes. The specific steps can be shown as follows [33,40,49]:
(1)
Suppose a MCDM problem with m alternatives X m ( m = 1 ,    2 ,    ,    M ) , and n criteria Y n ( n = 1 ,    2 ,    ,    N ) . All indicators in MCDM are interval numbers, and thus an interval-number decision matrix A ± = ( a m n ± ) M × N can be established.
(2)
Normalize the decision matrix to obtain a standardized interval-number decision matrix R ± = ( r m n ± ) M × N as follows:
[ r m n , r m n + ] = [ a m n , a m n + ] A n , m = 1 ,    2 ,    ,    M ,        n = 1 ,    2 ,       , N
where A n is the norm of the nth column vector and we have: A n = max ( max ( | a 1 n | , | a 1 n + | ) ,    max ( | a 2 n | , | a 2 n + | ) ,    ,    max ( | a M n | , | a M n + | ) ) ,         n = 1 ,    2 ,    ,    N . Please note the reverse indicators need to be converted to positive indicators as: [ c m n ,    c m n + ] = [ a m n + ,    a m n ] . Therefore, all the metrics are positive metrics.
(3)
Invite experts to compare the importance of indicators. According to the experts’ opinions, a judgment matrix can be built as follows with the interval numbers [31]:
B = [         [ 1 , 1 ]             [ b 12 , b 12 + ]                           [ b 1 N ,    b 1 N + ] [ 1 b 12 + , 1 b 12 ]           [ 1 , 1 ]                                   [ b 2 N , b 2 N + ]                                                                                                       [ 1 b 1 N + , 1 b 1 N ]     [ 1 b 2 N + , 1 b 2 N ]                               [ 1 , 1 ] ]
where b m n = [ b m n , b m n + ] , 1 b m n b m n + 9 , B ± = ( b m n ± ) N × N .
(4)
Calculate the eigenvectors x n , x n + , which corresponding to the largest eigenvalues of the judgment matrices B and B + , respectively. Then calculate the interval weights using the interval number eigenvalue method as follows [50].
w n ± = [ k x n , l x n + ]
where k = n = 1 N 1 m = 1 N b m n + , l = n = 1 N 1 m = 1 N b m n .
(5)
Judgment matrix consistency test. If 0 k 1 l , the consistency of the interval judgment matrix B is reasonable [36]. If k > 1 or l < 1, the consistency of the judgment matrix is poor. The results should be sent back to the experts, and the judgment matrix should be corrected or re-judged until the consistency requirements are met.
(6)
Multiply the columns of normalized decision matrix by the associated weights:
S ± = ( s m n ± ) M × N , [ s m n , s m n + ] = [ w n , w n + ] × [ r m n , r m n + ] ,     m = 1 ,    2 , ,    M ,      n = 1 , 2 , , N
It can be calculated based on the calculation principle of interval number multiplication [49].
(7)
Respectively determine the positive ideal and negative ideal solutions:
t n = max 1 m M s m n ,    t n + = max 1 m M s m n +
p n = min 1 m M s m n ,    p n + = min 1 m M s m n +
The positive ideal solution is X + = ( [ t 1 , t 1 + ] , [ t 2 , t 2 + ] , , [ t N , t N + ] ) ; the negative ideal solution is X = ( [ p 1 , p 1 + ] ,    [ p 2 , p 2 + ] , , [ p N , p N + ] ) .
(8)
Calculate the distances between each solution alternative and the ideal solution.
(9)
The distance from the scheme alternative X m to the positive ideal solution is:
d m + = d ( X m , X + ) = n = 1 N ( d m n + ) 2
where d m n + = max ( | s m n t m n | ,    | s m n + t m n + | ) .
The distance from the scheme alternative X m to the negative ideal solution is:
d m = d ( X m , X ) = n = 1 N ( d m n ) 2
where d m n = max ( | s m n p m n | , | s m n + p m n + | ) .
(10)
Calculate the proximity of each scenario to a positive ideal.
C m = d m d m + d m +
(11)
Rank the alternatives according to the relative closeness to the ideal alternatives, the larger value of C m means better effects of the corresponding scheme.

3. Case Study

3.1. Study Area

Shiyang River Basin (101°41’–104°16’ E, 36°29’–39°27’ N), as shown in Figure 2, is one of three inland rivers in the Hexi Corridor of Gansu Province in northwest China. The Shiyang River Basin has a typical arid continental climate, which is featured by low and irregular rainfalls, high evaporation rates, and long drought periods [51,52,53]. Since the 1960s, the population of the Shiyang River Basin has increased significantly, the scale of agricultural production has expanded, and the underground water resources in the basin have been overexploited. These changes caused the deterioration of the ecological environment and water shortage problems. In order to restore the ecological environment of the Shiyang River Basin, the local government formulated a series of policies such as the “Shiyang River Basin Key Management Plan” [54] and a series of supportive policies. However, in recent years, due to the decreased river runoff, the water shortage problems in the middle reaches of the Shiyang River Basin are becoming more and more serious. Irrigation is the largest water consumer in the middle reaches of the Shiyang River Basin, accounting for 88.28% of the total water consumptions [29]. The low efficiency of irrigation water usage further deepens the gap between water demands and supplies. Therefore, the practical allocation methods of irrigation water resources are highly desired to improve water usage efficiency in the middle reaches of the Shiyang River Basin.
The study area was a mountain plain located in the middle reaches of the Shiyang River Basin, which is called the Wuwei Basin. It includes two districts: The Liangzhou District of Wuwei City and the Yongchang County of Jinchang City. The dominated industry of the Wuwei Basin is agriculture and the local economy is relatively undeveloped. There are six rivers in the Wuwei Basin, namely the Dongda River, the Xiying River, the Jinta River, the Zamu River, the Huangyang River, and the Gulang River. These rivers can provide agricultural irrigation water for six irrigation districts located in the south of the study area, including Gulang (GL), Huangyang (HY), Zamu (ZM), Jinta (JT), Xiying (XY), and Donghe (DH). In addition, there are four groundwater irrigation districts in the north of the study area, including Qingyuan (QY), Jingyang (JY), Yongchang (YC), and Qinghe (QH) [55]. The surface water flows and groundwater withdrawals can be represented as stochastic distributions. Aiming to allocate the limited irrigation water resources in the Wuwei Basin, the decision makers from different hierarchical levels are responsible for making decision preferences from different viewpoints. The upper-level decision makers always rank the economy development as the primary goal, while the lower-level decision makers always focus on the fairness of water allocation among the irrigation districts to avoid social conflicts. Therefore, a comprehensive planning approach based on the ILFBP model and the IAHP-TOPSIS method would be applied to the study area to help water managers formulate well-rounded water allocation schemes under uncertainty.

3.2. Data Collection and Processing

The inputs of the developed model include the estimated water demands for irrigation districts, the available supplies of surface water and groundwater, the relevant crop production data, the regional meteorological data, and the economic data. In detail, the available supplies of groundwater are based on the collected data from the Shiyang River Basin Key Governance Projects [56]. The available supplies of surface water are based on historical data of the six rivers’ annual runoffs from 1955 to 2009. By fitting P-III curve based on the original runoff data (Figure 3), the runoffs of six rivers can be classified as different flow levels as: Specifically, P < 12.5 % corresponds to the very high flow level, 12.5 % < P < 37.5 % corresponds to the high flow level, 37.5 % < P < 62.5 % corresponds to the middle flow level, 62.5 % < P < 87.5 % corresponds to the low flow level, and P > 87.5 % corresponds to the very low flow level (where P = m n + 1 × 100 % , m represents the number of years with annual runoffs not less than a certain flow value, and n represents the total number of years with the flow data). Therefore, the probability levels (i.e., 0.125, 0.25, 0.25, 0.25, and 0.125) and the corresponding water amounts under the five flow levels can be obtained. In this study, a confidence interval with a significance level of 0.05 is used to obtain the interval ranges of water supplies. Table 1 shows the available surface water supplies of six rivers under different flow levels.
The irrigation water demands can be obtained according to E T c = K c × E T 0 , where E T c is the actual crop evapotranspiration (mm/day), K c is the crop coefficient, and E T 0 is the reference crop evapotranspiration (mm/day). The E T 0 value can be calculated based on the FAO 56 Penman-Monteith method [57]. The meteorological data of the Wuwei Sub-basin from 1956 to 2017 are obtained from the China Meteorological Data Sharing Service System. Annual rainfall data from 1956 to 2017 are collected from the observation stations in Wuwei and the effective rainfall coefficient is 0.8 [58]. The effective rainfall is estimated by the product of actual rainfall and the effective precipitation coefficient. The main crops in the Wuwei Basin include wheat, corn, and economic crops, and their market prices come from the 2013–2016 China Agricultural Products Price Survey Yearbook. The price of wheat is [2.44, 2.57] yuan/kg, the price of corn is [2.09, 2.47] yuan/kg, and the price of economic crops is [4.94, 5.34] yuan/kg. Effective irrigation areas of crops are from the previous case study [55].
Table 2 shows the other parameters of the irrigation districts needed to solve the developed optimization model, such as the population data, irrigation water usage efficiencies, water productivities of different crops, areas of irrigation districts, and water prices. According to the water price adjustment policy of Wuwei City, the water price consists of two parts: The metering water price and the basic water price. For surface water supplies, the basic and metering water prices are respectively 30 yuan/ha and 0.157 yuan/m3; for groundwater supplies, the basic and metering water price are respectively 120 yuan/ha and 0.01 yuan/m3.

4. Results Analysis and Discussion

4.1. The ILFBP Model Results of the Water Allocation

Figure 4 shows the optimization allocation results of each irrigation district under different flow levels. Obviously, with the decreased incoming water, the configuration target was reduced, and the water allocation in each irrigation district was also gradually reduced.
Except for the irrigation districts of Xiying, Donghe, and Qinghe, the water allocation results for other irrigation districts could reach the allocation target. The Xiying irrigation district has the largest planting area, thus faces the most serious water shortage. The main water source for irrigation of the Xiying irrigation district is the Xiying River. The Xiying River is the largest tributary of the Shiyang River system, and the Xiying Reservoir is located at the exit of the mountain. The Xiying Reservoir not only provides irrigation water for the Xiying irrigation district, but also has the functions of preventing drought and flood, supplying safe drinking water for humans and animals, and transferring water to downstream. The current irrigation schemes cannot meet the water demands, resulting in serious waste of water resources and low irrigation water usage efficiency. These reasons basically answer the key question why water shortage often occurs in the Xiying irrigation district.
Figure 5 shows the water allocation results for different crops in different irrigation districts under the medium flow level. The regional water allocation for wheat, maize, and the economic crops were respectively [5.25, 5.69] × 108 m3, [1.23, 1.34] × 108 m3, and [0.23, 0.25] × 108 m3. Aggressive decision makers may choose the upper bounds of water allocations to obtain higher food production and economic benefits, but they also need to bear relatively higher risks; conservative decision makers may choose the lower bounds of water allocations with relatively lower levels of benefits and risk. In general, the water allocation of food crops production is far greater than that economic crops, which has the similar trends with other previous studies in this study area [47,55]. For modern agricultural development, the water managers of irrigation districts should be able to adjust the regional planting structure and the irrigation modes in time. The water management department should guide local farmers to plant more economic crops to improve water usage efficiency. Additionally, the water supplies of economic crops should be increased to improve the agricultural economic benefits in irrigation districts.
Figure 6 shows the comparison of the total regional water allocation among the current situation, the results of the ILFBP model, and the results of the economic model without considering the system fairness. The total water allocation of the ILFBP model in the study area was [6.73, 7.37] × 108 m3, while that of the economic model is [7.04, 7.97] × 108 m3 and that of the status quo was [7.57, 8.37] × 108 m3. The total water allocation of the ILFBP model was the least one among the three schemes. Except for Xiying and Donghe irrigation districts, the water allocation of the economic model and ILFBP model was almost the same. The crop planting area and irrigation water allocation in the Xiying and Donghe irrigation districts were the largest among the study areas. However, the population of these two irrigation districts was not the largest. Therefore, in order to improve the water allocation fairness between irrigation districts, the ILFBP model appropriately reduced the water allocation of the Xiying and Donghe irrigation district.
Figure 7 shows the economic benefits of the three water allocation schemes. The economic benefits of the economic model, ILFBP model, and the status quo were respectively [16.49, 23.06] × 108 yuan, [16.01, 21.73] × 108 yuan, and [15.66, 22.04] × 108 yuan. The economic model only maximized the regional overall economic benefits brought by the water allocation, so more water was allocated to the irrigation areas with more economic crops to improve the overall economic benefits. Therefore, the economic model could bring higher economic benefits. The benefits of the actual situation and the ILFBP model were similar, but the interval ranges of the system benefits were different. The upper-bound benefits of the status quo were larger than those of the ILFBP model, indicating that the current situation could obtain greater benefits. However, the larger interval range of the benefits under the status quo was larger indicates greater level of uncertainties existed in the benefits of the status quo. It seems that more irrigation water could bring more economic benefits, but the marginal benefits existed in the relationship between economic benefits and water allocation amounts, which may lead to the decreased water usage efficiency when using more irrigation water. Therefore, the water allocation schemes could not be evaluated only through comparing the amount of irrigation water and economic benefits. An appropriate evaluation method is still needed to compare the effects of water allocation schemes and provide more useful decision-making information.

4.2. Evaluation Results

In order to reflect the economic, social, and resource effects of water-allocation schemes, the evaluation model indicators such as the Gini coefficient, the irrigation water productivity, the benefits per unit of water usage, the irrigation loss per unit of area, and the water shortage rate were calculated based on the results obtained from three optimization models. The values of the indicators of the three schemes are shown in Table 3.
Y 1 ,    Y 5 represent the above-mentioned five indicators, respectively, and the weights of the five indicators are represented by the weight vector w. The IAHP was used for comparing the estimation indicators to build an interval judgment matrix and calculate the results of the index weights as shown in Table 4. The calculation results of k = 0.92 < 1 and l = 1.06 > 1 could pass the consistency test. The weighted normalization decision matrix is shown in Table 5.
The positive ideal solution determined the ITOPSIS method is X + = { [ 0.21 , 0.17 ] , [ 0.13 , 0.15 ] , [ 0.35 , 0.50 ] , [ 0.07 , 0.05 ] , [ 0 , 0 ] } , and the negative ideal solution is X = { [ 0.28 , 0.23 ] , [ 0.11 , 0.13 ] , [ 0.28 , 0.45 ] , [ 0.08 , 0.06 ] , [ 0.04 , 0.02 ] } . The distance between each solution to the positive and negative ideal solutions is: d 1 + = 7.06 × 10 2 , d 1 = 8.15 × 10 2 , d 2 + = 6.54 × 10 2 , d 2 = 7.97 × 10 2 , d 3 + = 7.69 × 10 2 , d 3 = 8.06 × 10 2 . The relative closeness of each scheme is C 1 = 5.36 × 10 1 , C 2 = 5.49 × 10 1 , C 3 = 5.12 × 10 1 . From the above results, the ranking result of the three schemes can be obtained as C 2 > C 1 > C 3 .
The evaluation results show that the applicability of the ILFBP model is better than other two models in this study area. The regional economic, resource, and social indicators can be improved by adjusting the water allocation schemes for different irrigation districts. Among the three schemes, the state quo showed the highest level of fairness among the three schemes with the Gini coefficient of only [0.29, 0.31], but the regional economic growth and local irrigation water productivity were also limited by strong fairness concerns. In the economic model, in order to obtain greater economic benefits, more water was distributed to the irrigation districts with more economic crops, such as Xiying, Donghe, and Qinghe. However, this leads to unfair water distribution between irrigation districts. The Gini coefficient of the economic model was [0.39, 0.41], which was close to triggering the warning sign. The ILFBP model coordinates the decision-making preferences of two-level stakeholders with the considerations of system economic benefits and system fairness. These efforts contribute to the social contradiction caused by the extremely unfair water allocation among irrigation districts. The Gini coefficient of the ILFBP scheme was lower than 0.4. The superiority of the ILFBP model could also be reflected through the indicator comparison, such as irrigation loss, irrigation water productivity, and benefit per unit water. These results indicate that the water allocation scheme obtained by ILFBP could effectively reduce unnecessary water waste and make most of the indicators closer to the positive ideal solution.

4.3. Discussion

This study used the ILFBP model as a planning management tool to improve irrigation-water use efficiency and further increase available ecological water in the Shiyang River Basin. The results obtained from ILFBP model provided specific adjustment and allocation recommendations under different flow levels, which contributed to the sustainable development of local agriculture.
The effects evaluation of the water allocation schemes could test the applicability of the ILFBP model to the specific regions, and further improve the implementation significance of the optimization schemes for practical problems. The IAHP-TOPSIS evaluation method could deal with uncertainties and provide more practical strategies than existing deterministic evaluation models. The uncertainties associated with the water allocation variables and the subjective preferences could be fully reflected in the expert evaluation process.
In order to evaluate the effects of the ILFBP model comprehensively, two models, respectively the status quo situation and the single economic goal without the fairness concern, were built and solved. The obtained results showed that the ILFBP model could coordinate the economic benefits and fairness of system under uncertainty. Moreover, it could improve the water allocation effects from multiple perspectives, and obtained an irrigation-water allocation scheme with a high satisfaction degree for different stakeholders. The water allocation schemes obtained from the ILFBP model were obviously more reasonable than the status quo, and could meet the policy making requirements of the Shiyang River Basin. These efforts could help local water managers of the Shiyang River Basin formulate more suitable water policies.
Compared with other previous studies on water resources allocation in the Shiyang River Basin [29,47,51], the ILFBP model could obtain a satisfied irrigation-water allocation scheme by considering the system fairness and making tradeoff between the upper- and lower-level decision makers under uncertainty. Moreover, the optimization-evaluation framework developed in this study could not only provide specific water allocation schemes, but also evaluate the effectiveness and the applicability of the optimization schemes, making the results more convincing and more realistic.
Furthermore, the optimization-evaluation framework developed in this study could be applied to other regions for improving water resource management at a government, regional, local, and basin level. The specific planning benefits and the evaluation results are associated with irrigation water allocation, ecological water allocation, different flow levels, and interval range of social and economic parameters. These data can be obtained through field research or model programming. Therefore, the proposed optimization-evaluation framework under uncertainty is universal and has significant implications for supporting regional water allocation practices, especially in arid and semi-arid regions.

5. Conclusions

In this study, an interval linear fractional bi-level programming (ILFBP) model was developed for irrigation-water allocation planning under uncertainty. The bi-level planning (BP), interval planning (IP) and linear fractional programming (LFP) techniques are integrated to deal with the problems in reality. The interactive fuzzy coordination algorithm was used for solving the developed model, and the interval solutions under different flow levels could be obtained. In order to reasonably evaluate the effects of the optimization schemes under uncertainty, the interval analytic hierarchy process (IAHP) method and the interval TOPSIS method were simultaneously applied for evaluating the obtained allocation schemes. The developed optimization-evaluation approach was applied to a real-world case study in the middle reaches of the Shiyang River Basin. Compared with the previous studies, the developed framework had the following advantages. (1) It could consider the decision-making preferences of multi-level stakeholders and reflect different priorities existed in the decision-making process. The goals of the upper-level decision makers could be prioritized, and then be weighed among the two levels of decision makers, and finally the decision-making plan with higher overall satisfaction could be obtained; (2) the fairness goal could be added to the water resources optimization model to ensure system fairness of allocation schemes under uncertainty; (3) it could quantify the subjective intention of the experts and reflect the uncertainty in this process; and (4) it could provide the corresponding effects evaluation for optimization schemes under interval uncertainty. The proposed approach could help decision makers develop and identify more appropriate water-irrigation schemes. This framework could also be applied to other practical problems of resources and environmental management for developing well-rounded strategies under various system complexities.

Author Contributions

Conceptualization, F.Z.; Data curation, Y.Z.; Formal analysis, Y.Z., F.Z. and H.Z.; Methodology, Y.Z.; Project administration, P.G.; Resources, P.G.; Software, Y.Z.; Validation, F.Z. and P.G.; Writing – review & editing, Y.Z., F.Z., H.Z. and P.G.

Funding

This research was funded by the National Natural Science Foundation of China grant number [41871199, 51621061].

Acknowledgments

The authors are grateful to the editor and the anonymous reviewers for their insightful comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AbbreviationMeaning and Description
ILFBPInterval linear fractional bi-level programming
IAHPInterval analytic hierarchy process
ITOPSISInterval-parameter TOPSIS method
MCDMMulti-criteria decision making
ULDMUpper-level decision maker
LLDMLower-level decision maker
FPFuzzy programming
SPStochastic programming
IPInterval programming
BPBi-level programming
ILPInterval linear programming
ILFPInterval linear fractional programming
ILBPInterval linear bi-level programming
GLGulang irrigation district
HYHuangyang irrigation district
ZMZamu irrigation district
JTJinta irrigation district
XYXiying irrigation district
DHDonghe irrigation district
QYQingyuan irrigation district
JYJingyang irrigation district
YCYongchang irrigation district
QHQinghe irrigation district

Parameters

The ILFBP Methods
Parameters and VariablesMeaning and Description
UUpper-level
LLower-level
optThe optimal values
x 1 ,    x 2 The upper- and lower-level decision variables
f 1 ,   f 2 The upper- and lower-level objective functions
c11, c21, c12, c22, d11, d21, d12, d22Interval parameters
GThe set of bi-level constraints
r 1 The maximum tolerances of the upper-level decision variables
f 1 , f 2 The worst values of the upper- and lower-level objective functions
μ x 1 ( x 1 ) The satisfaction degrees of the upper-level decision variables
μ f 1 [ f 1 ( x ) ] , μ f 2 [ f 2 ( x ) ] ,The satisfaction degrees of the upper- and lower-level objective functions
λ The overall satisfactory degree of the ILFBP model
IA column vector with the same dimension as x and all elements of I equal to 1
The application of the ILFBP model
Parameters and VariablesMeaning and Description
iIrrigation districts
jCrop types
hFlow levels
W i j h The amount of irrigation water allocated to the crop j in the irrigation district i under the flow level h (m3), the decision variables
p h The possibilities of different flow levels
B i j The prices of crop j in the irrigation district i (yuan/kg)
Y i j The water productivities of crop j in the irrigation district i (kg/m3)
A i Irrigated areas of subarea i (ha)
k i The irrigation water use efficiency
P i j h Effective precipitation (m3)
J B i Basic water price (yuan/ha)
J L i The metering water price (yuan/m3)
R i Populations of the irrigation district i
T W h The total water supplies under the flow level h (m3)
T W i h Water supplies for the irrigation district i under flow level h (m3)
I R i j h , min ± , I R i j h , max ± The minimum and maximum crop water requirement of crop j in the irrigation district i under flow level h (m3)
The IAHP-TOPSIS evaluation methods
Parameters and VariablesMeaning and Description
X m Alternatives ( m = 1 , 2 , , M )
Y n Indicators ( n = 1 ,    2 , , N )
A ± = ( a m n ± ) M × N Interval-number decision matrix
R ± = ( r m n ± ) M × N Standardized interval-number decision matrix
A n The norm of the jth column vector
B ± = ( b m n ± ) N × N Interval judgment matrix
w n ± The interval weights
x n ± Eigenvectors
k, lParameters required to calculate interval weights
S ± = ( s m n ± ) M × N Weighted normalized decision matrix
d m + , d m The distance from the scheme alternative X m to the positive and the Negative ideal solution, respectively
X + , X The positive ideal and negative ideal solutions
C m The relative closeness

References

  1. Xu, Q.; Chen, Q.; Zhao, S.; Liu, K.; Ma, J. Saving Water and Associated Energy from Distribution Networks by Considering Landscape Factors in Pressure Management and Use of District Metered Areas. J. Environ. Inform. 2017, 31, 65–73. [Google Scholar] [CrossRef]
  2. Tan, Q.; Zhang, S.; Li, R. Optimal use of agricultural water and land resources through reconfiguring crop planting structure under socioeconomic and ecological objectives. Water 2017, 9, 488. [Google Scholar] [CrossRef]
  3. Kang, S.; Hao, X.; Du, T.; Tong, L.; Su, X.; Lu, H.; Li, X.; Huo, Z.; Li, S.; Ding, R. Improving agricultural water productivity to ensure food security in China under changing environment: From research to practice. Agric. Water Manag. 2017, 179, 5–17. [Google Scholar] [CrossRef]
  4. Hrozencik, R.A.; Manning, D.T.; Suter, J.F.; Goemans, C.; Bailey, R.T. The heterogeneous impacts of groundwater management policies in the Republican River Basin of Colorado. Water Resour. Res. 2017, 53, 10757–10778. [Google Scholar] [CrossRef]
  5. Al-Isawi, R.H.; Scholz, M.; Al-Faraj, F.A. Assessment of diesel-contaminated domestic wastewater treated by constructed wetlands for irrigation of chillies grown in a greenhouse. Environ. Sci. Pollut. Res. 2016, 23, 25003–25023. [Google Scholar] [CrossRef] [Green Version]
  6. Al-Isawi, R.H.K.; Almuktar, S.A.A.-A.N.; Scholz, M. Recycling of river, rain, gully pot and grey waters for irrigating chillies. Environ. Monit. Assess. 2016, 188, 1–12. [Google Scholar] [CrossRef] [PubMed]
  7. Sethi, L.N.; Kumar, D.N.; Panda, S.N. Optimal crop planning and conjunctive use of water resources in a coastal river basin. Water Resour. Manag. 2002, 16, 145–169. [Google Scholar] [CrossRef]
  8. Guo, P.; Huang, G.H.; He, L.; Sun, B.W. ITSSIP: Interval-parameter two-stage stochastic semi-infinite programming for environmental management under uncertainty. Environ. Model. Softw. 2008, 23, 1422–1437. [Google Scholar] [CrossRef]
  9. Li, M.; Guo, P.; Fang, S.Q.; Zhang, L.D. An inexact fuzzy parameter two-stage stochastic programming model for irrigation water allocation under uncertainty. Stoch. Environ. Res. Risk Assess. 2013, 27, 1441–1452. [Google Scholar] [CrossRef]
  10. Li, Y.P.; Huang, G.H.; Guo, P.; Yang, Z.F.; Nie, S.L. A dual-interval vertex analysis method and its application to environmental decision making under uncertainty. Eur. J. Oper. Res. 2010, 200, 536–550. [Google Scholar] [CrossRef]
  11. Huang, G.H. A hybrid inexact-stochastic water management model. Eur. J. Oper. Res. 1998, 107, 137–158. [Google Scholar] [CrossRef]
  12. Zhang, C.; Li, M.; Guo, P. An interval multistage joint-probabilistic chance-constrained programming model with left-hand-side randomness for crop area planning under uncertainty. J. Clean. Prod. 2017, 167, 1276–1289. [Google Scholar] [CrossRef]
  13. Zhang, D.; Guo, P. Integrated agriculture water management optimization model for water saving potential analysis. Agric. Water Manag. 2016, 170, 5–19. [Google Scholar] [CrossRef]
  14. Madani, K. Game theory and water resources. J. Hydrol. 2010, 381, 225–238. [Google Scholar] [CrossRef]
  15. Camacho-Vallejo, J.F.; Gonzalez-Rodriguez, E.; JavierAlmaguer, F.; Gonzaez-Ramirez, R.G. A bi-level optimization model for aid distribution after the occurrence of a disaster. J. Clean. Prod. 2015, 105, 134–145. [Google Scholar] [CrossRef]
  16. Lv, Y.; Huang, G.; Li, Y.; Yang, Z. Planning regional water resources system using an interval fuzzy bi-level programming method. J. Environ. Inform. 2010, 16, 43–56. [Google Scholar] [CrossRef]
  17. Chen, Y.; Hongwei, L.U.; Jing, L.I.; Lixia, R.E.N.; Li, H.E. Empirical research on bi-level allocation model of water resources under uncertainty (in Chinese). Acta Sci. Circumstantiae 2016, 36, 2252–2261. [Google Scholar]
  18. Hu, Z.; Wei, C.; Yao, L.; Li, L.; Li, C. A multi-objective optimization model with conditional value-at-risk constraints for water allocation equality. J. Hydrol. 2016, 542, 330–342. [Google Scholar] [CrossRef]
  19. Zhang, C.; Engel, B.A.; Guo, P.; Zhang, F.; Guo, S.; Liu, X.; Wang, Y. An inexact robust two-stage mixed-integer linear programming approach for crop area planning under uncertainty. J. Clean. Prod. 2018, 204, 489–500. [Google Scholar] [CrossRef]
  20. Feng, T.; Yang, Z. A New Algorithm for Weighted Proportional Fairness Based Spectrum Allocation of Cognitive Radios. In Proceedings of the 8th Acis International Conference on Software Engineering, Melbourne, Australia, 11–13 July 2007. [Google Scholar]
  21. Uchida, M.; Kurose, J. An Information-Theoretic Characterization of Weighted alpha-Proportional Fairness. In Proceedings of the IEEE INFOCOM, Rio de Janeiro, Brazil, 19–25 April 2009. [Google Scholar]
  22. Alvaredo, F. A note on the relationship between top income shares and the Gini coefficient. Econ. Lett. 2011, 110, 274–277. [Google Scholar] [CrossRef]
  23. Karagiannis, E.; Kovacevic, M. A Method to Calculate the Jackknife Variance Estimator for the Gini Coefficient. Oxf. Bull. Econ. Stat. 2010, 62, 119–122. [Google Scholar] [CrossRef]
  24. Zhang, Y.B.; Mai, Z.Q.; Chen, X.G. Analysis of city resource-environment Gini coefficient in Guangdong Province. Acta Ecol. Sin. 2008, 49, 29–38. [Google Scholar]
  25. Chotikapanich, D. A comparison of alternative functional forms for the Lorenz curve. Econ. Lett. 1993, 41, 129–138. [Google Scholar] [CrossRef]
  26. Qi, L. Gini-coefficient of Chinese manufacturing industry. Stat. Res. 2003, 9, 21–25. [Google Scholar]
  27. Gregorio, J.D.; Lee, J.W. Education and income inequality: New evidence from cross-country data. Rev. Income Wealth 2002, 48, 395–416. [Google Scholar] [CrossRef]
  28. Zhang, C.; Huan, Z.; Li, H. An Empirical Study on Education Equality Based on Education Gini Coefficient in China: 1978–2004. Tsinghua J. Educ. 2006, 27, 10–14. (In Chinese) [Google Scholar]
  29. Ren, C.; Guo, P.; Tian, Q.; Zhang, L. A multi-objective fuzzy programming model for optimal use of irrigation water and land resources under uncertainty in Gansu Province, China. J. Clean. Prod. 2017, 164, 85–94. [Google Scholar] [CrossRef]
  30. Zhang, C.; Li, M.; Guo, P. Two-stage stochastic chance-constrained fractional programming model for optimal agricultural cultivation scale in an arid area. J. Irrig. Drain. Eng. 2017, 143, 05017006. [Google Scholar] [CrossRef]
  31. Zhang, F.; Guo, S.; Ren, C.; Guo, P. Integrated IMO-TSP and AHP Method for Regional Water Allocation under Uncertainty. J. Water Resour. Plan. Manag. 2018, 144, 04018025. [Google Scholar] [CrossRef]
  32. Zhang, Y.; Guo, P.; Zhang, F.; Tang, Y. Applicability comparison of water resources optimization model based on effect evaluation. J. Drain. Irrig. Mach. Eng. 2019, 37, 154–159. (In Chinese) [Google Scholar]
  33. Zhang, S.; Sun, B.; Yan, L.; Wang, C. Risk identification on hydropower project using the IAHP and extension of TOPSIS methods under interval-valued fuzzy environment. Nat. Hazards 2013, 65, 359–373. [Google Scholar] [CrossRef]
  34. Sugihara, K.; Tanaka, H. Interval evaluations in the analytic hierarchy process by possibility analysis. Comput. Intell. 2001, 17, 567–579. [Google Scholar] [CrossRef]
  35. Zhang, C.; Su, X.; Cai, H. Harmony Evaluation of Water Resources Allocation System in Yuyang Area Based on IAHP. J. Nat. Resour. 2017, 32, 105–1063. (In Chinese) [Google Scholar]
  36. Xu, Z.; Da, Q. A Possibility-based Method for Priorities of Interval Judgement Matrices. Chin. J. Manag. Sci. 2003, 11, 63–65. [Google Scholar]
  37. Zhang, Q.; Fan, Z.; Pan, D. A Ranking Approach for Interval Numbers in Uncertain Multiple Attribute Decision Making Problems. Syst. Eng. Theory Pract. 1999, 19, 129–133. (In Chinese) [Google Scholar]
  38. Jahanshahloo, G.R.; Lotfi, F.H.; Izadikhah, M. An algorithmic method to extend TOPSIS for decision-making problems with interval data. Appl. Math. Comput. 2006, 175, 1375–1384. [Google Scholar] [CrossRef]
  39. Dymova, L.; Sevastjanov, P.; Tikhonenko, A. A direct interval extension of TOPSIS method. Expert Syst. Appl. 2013, 40, 4841–4847. [Google Scholar] [CrossRef]
  40. Tsaur, R. Decision risk analysis for an interval TOPSIS method. Appl. Math. Comput. 2011, 218, 4295–4304. [Google Scholar] [CrossRef]
  41. Jin, S.W.; Li, Y.P.; Nie, S. An integrated bi-level optimization model for air quality management of Beijing’s energy system under uncertainty. J. Hazard. Mater. 2018, 350, 27–37. [Google Scholar] [CrossRef]
  42. Ahlatcioglu, M.; Tiryaki, F. Interactive fuzzy programming for decentralized two-level linear fractional programming (DTLLFP) problems. Omega 2007, 35, 432–450. [Google Scholar] [CrossRef]
  43. Emam, O.E. A fuzzy approach for bi-level integer non-linear programming problem. Appl. Math. Comput. 2006, 172, 62–71. [Google Scholar] [CrossRef]
  44. Zhu, H.; Huang, W.W.; Huang, G.H. Planning of regional energy systems: An inexact mixed-integer fractional programming model. Appl. Energy 2014, 113, 500–514. [Google Scholar] [CrossRef]
  45. Guo, P.; Huang, G.H.; He, L.; Li, H.L. Interval-parameter Fuzzy-stochastic Semi-infinite Mixed-integer Linear Programming for Waste Management under Uncertainty. Environ. Model. Assess. 2009, 14, 521–537. [Google Scholar] [CrossRef]
  46. Zhang, F.; Tan, Q.; Zhang, C.; Guo, S.; Guo, P. A Regional Water Optimal Allocation Model Based on the Cobb-Douglas Production Function under Multiple Uncertainties. Water 2017, 9, 923. [Google Scholar] [CrossRef]
  47. Wang, H.; Zhang, C.; Guo, P. An Interval Quadratic Fuzzy Dependent-Chance Programming Model for Optimal Irrigation Water Allocation under Uncertainty. Water 2018, 10, 684. [Google Scholar] [CrossRef]
  48. Xiao, J.; Wang, C.; Zhou, M. An IAHP-based MADM method in urban power system planning. Proc. CSEE 2004, 4, 50–57. (In Chinese) [Google Scholar]
  49. Zhang, J.; Liu, J. The Study of Decision Making Method to Interval Number Multiple Attribute Decision Making Problem. Forecasting 2002, 21, 73–75. (In Chinese) [Google Scholar]
  50. Wei, Y.; Liu, J.; Wang, X. Concept of consistence and weights of the judgement matrix in uncertain type of AHP. Syst. Eng. Theory Pract. 1994, 14, 16–22. [Google Scholar]
  51. Zhang, F.; Li, M.; Guo, S.; Zhang, C.; Guo, P. Integrated uncertain models for runoff forecasting and crop planting structure optimization of the Shiyang River Basin, north-west China. Front. Agric. Sci. Eng. 2018, 5, 177–187. [Google Scholar] [CrossRef]
  52. Wang, Q.; Shi, J.; Chen, G.; Xue, L. Environmental effects induced by human activities in arid Shiyang River basin, Gansu province, northwest China. Environ. Geol. 2002, 43, 219–227. [Google Scholar]
  53. Shi, Q.; Chen, F.H.; Zhu, Y.; Madsen, D. Lake evolution of the terminal area of Shiyang River drainage in arid China since the last glaciation. Quat. Int. 2002, 93, 31–43. [Google Scholar] [CrossRef]
  54. Chen, M. Thoughts concerning the Shiyanghe River Basin improvement. China Water Resour. 2010, 1, 36–38. (In Chinese) [Google Scholar]
  55. Zhang, L.; Guo, P.; Fang, S.; Li, M. Monthly Optimal Reservoirs Operation for Multicrop Deficit Irrigation under Fuzzy Stochastic Uncertainties. J. Appl. Math. 2014, 1–11. [Google Scholar] [CrossRef]
  56. Gansu Provincial Department of Water Resources. Shiyang River Basin Key Governance Projects; Gansu Provincial Department of Water Resources: Gansu, China, 2007. (In Chinese)
  57. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements-FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998; Volume 156, p. 178. [Google Scholar]
  58. Xu, X.; Zhou, H.; Wang, Z.; Jiapaer, R. Study on effective rainfall use efficiency in arid irrigation district. Water Sav. Irrig. 2010, 12, 44–46. [Google Scholar]
Figure 1. Model framework.
Figure 1. Model framework.
Water 11 01094 g001
Figure 2. The study area.
Figure 2. The study area.
Water 11 01094 g002
Figure 3. P-III frequency distribution curves of runoff from the six rivers.
Figure 3. P-III frequency distribution curves of runoff from the six rivers.
Water 11 01094 g003aWater 11 01094 g003b
Figure 4. Water allocation at different flow levels in different irrigation districts.
Figure 4. Water allocation at different flow levels in different irrigation districts.
Water 11 01094 g004
Figure 5. Water allocation for all the crops in different irrigation districts under the medium flow level.
Figure 5. Water allocation for all the crops in different irrigation districts under the medium flow level.
Water 11 01094 g005
Figure 6. The total water allocation of the three water allocation schemes in each irrigation district. Note: m = 1, m = 2, and m = 3 represent the economic model, ILFBP model, and status quo, respectively.
Figure 6. The total water allocation of the three water allocation schemes in each irrigation district. Note: m = 1, m = 2, and m = 3 represent the economic model, ILFBP model, and status quo, respectively.
Water 11 01094 g006
Figure 7. Economic benefits of the three water allocation schemes.
Figure 7. Economic benefits of the three water allocation schemes.
Water 11 01094 g007
Table 1. The runoff of the six rivers in different flow levels (108 m3).
Table 1. The runoff of the six rivers in different flow levels (108 m3).
Water LevelGulang RiverHuangyang RiverZamu RiverJingta RiverXiying RiverDongda River
h = 1[1.02, 1.39][1.73, 1.87][3.00, 3.97][1.54, 1.98][4.53, 4.80][3.86, 4.03]
h = 2[0.73, 0.80][1.39, 1.54][2.60, 2.75][1.49, 1.57][4.07, 4.31][3.29, 3.43]
h = 3[0.59, 0.62][1.23, 1.27][2.26, 2.33][1.27, 1.31][3.57, 3.68][2.87, 2.96]
h = 4[0.47, 0.52][0.99, 1.07][2.00, 2.07][1.09, 1.15][3.10, 3.24][2.65, 2.73]
h = 5[0.34, 0.39][0.76, 1.01][1.53, 1.71][0.87, 0.97][2.68, 2.89][2.07, 2.55]
Table 2. The main parameters of the developed model.
Table 2. The main parameters of the developed model.
Irrigation DistrictPopulation (104)Irrigation Water Usage EfficiencyWater Productivity/(kg/m3)Irrigation Areas (104 ha)Water Price
Food CropEconomic CropsBasic Water Price (yuan/ha)Metering Water Price (yuan/m3)
Gulang8.770.53[1.01, 1.14][1.61, 1.81]0.58300.157
Huangyang10.900.59[1.47, 1.65][1.47, 1.65]1.11300.157
Zamu12.880.58[1.39, 1.56][1.14, 1.27]0.94300.157
Jinta9.260.71[1.13, 1.26][1.11, 1.24]0.65300.157
Xiying10.720.61[1.73, 1.94][1.44, 1.61]1.86300.157
Donghe4.600.52[1.20, 1.36][1.81, 2.05]1.69300.157
Qingyuan11.650.75[1.29, 1.44][0.94, 1.04]0.741200.01
Jinyang10.520.72[0.97, 1.08][0.82, 0.91]0.511200.01
Yongchang11.470.60[1.25, 1.40][1.10, 1.23]0.691200.01
Qinghe5.140.69[0.77, 0.86][0.98, 1.09]0.831200.01
Table 3. Decision Matrix.
Table 3. Decision Matrix.
ModelsGini CoefficientIrrigation Water ProductivityBenefits Per Unit of Water UsageIrrigation Loss Per Unit of AreaWater Shortage Rate
/Kg/m3yuan/m3(103 m3)/ha/
Economic model[0.39, 0.41][1.27, 1.44][2.34, 2.89][2.90, 3.27]0
ILFBP model[0.35, 0.38][1.38, 1.51][2.38, 2.95][2.75, 3.01][0.04, 0.07]
Status quo[0.29, 0.31][1.21, 1.34][1.87, 2.63][2.97, 3.28]0
Table 4. Indicator judgment matrix and weight.
Table 4. Indicator judgment matrix and weight.
YY1Y2Y3Y4Y5w
Y1[1, 1][2, 3][1/4, 1/2][4, 6][6, 7][0.24, 0.28]
Y2[1/3, 1/2][1, 1][1/5, 1/4][2, 3][5, 6][0.14, 0.15]
Y3[2, 4][4, 5][1, 1][5, 7][6, 8][0.43, 0.50]
Y4[1/6, 1/4][1/3, 1/2][1/7, 1/5][1, 1][2, 3][0.07, 0.08]
Y5[1/7, 1/6][1/6, 1/5][1/8, 1/6][1/3, 1/2][1, 1][0.04, 0.04]
Table 5. Weighted normalized decision matrix.
Table 5. Weighted normalized decision matrix.
Distance to the Ideal SolutionY1Y2Y3Y4Y5
Economic model[−0.28, −0.23][0.12, 0.15][0.35, 0.50][−0.08, −0.06]0
ILFBP model[−0.26, −0.21][0.13, 0.15][0.35, 0.50][−0.07, −0.05][−0.04, −0.02]
Status quo[−0.21, −0.17][0.11, 0.13][0.28, 0.45][−0.08, −0.06]0

Share and Cite

MDPI and ACS Style

Zhang, Y.; Zhang, F.; Zhu, H.; Guo, P. An Optimization-Evaluation Agricultural Water Planning Approach Based on Interval Linear Fractional Bi-Level Programming and IAHP-TOPSIS. Water 2019, 11, 1094. https://doi.org/10.3390/w11051094

AMA Style

Zhang Y, Zhang F, Zhu H, Guo P. An Optimization-Evaluation Agricultural Water Planning Approach Based on Interval Linear Fractional Bi-Level Programming and IAHP-TOPSIS. Water. 2019; 11(5):1094. https://doi.org/10.3390/w11051094

Chicago/Turabian Style

Zhang, Yan, Fan Zhang, Hua Zhu, and Ping Guo. 2019. "An Optimization-Evaluation Agricultural Water Planning Approach Based on Interval Linear Fractional Bi-Level Programming and IAHP-TOPSIS" Water 11, no. 5: 1094. https://doi.org/10.3390/w11051094

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop