Next Article in Journal
Multi-Objective Optimization of a Pitch Point Absorber Wave Energy Converter
Previous Article in Journal
Instantaneous Water Purification by Deep Ultraviolet Light in Water Waveguide: Escherichia Coli Bacteria Disinfection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Optimization Method for the Layout of Pumping Wells in Oases: Application in the Qira Oasis, Northwest China

1
State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, China
2
Cele National Station of Observation and Research for Desert-Grassland Ecosystem, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, China
3
College of Mathematics and System Sciences, Xinjiang University, Urumqi 830046, China
4
School of Water Resources and Environment, China University of Geosciences, Beijing 100083, China
5
Institution of Groundwater and Earth Sciences, Jinan University, Guangzhou 510632, China
6
College of Construction Engineering, Jilin University, Changchun 130400, China
*
Authors to whom correspondence should be addressed.
Water 2019, 11(5), 970; https://doi.org/10.3390/w11050970
Submission received: 11 March 2019 / Revised: 28 April 2019 / Accepted: 6 May 2019 / Published: 9 May 2019
(This article belongs to the Section Hydrology)

Abstract

:
Oases are vital habitat areas for both humans and wild plants and animals in desert areas of arid Northwestern China. The efficient management of oasis water resources, especially groundwater resources, is very important for the environmental sustainability and economic development of the region. Pumping wells play a vital role in the oasis groundwater supply; therefore, optimizing the layout of these wells is essential for water resource management. In this study, we present a novel optimization methodology that implements a genetic algorithm and nonlinear programming model for the layout of pumping wells. The methodology was tested and evaluated in the real oasis case study of Qira Oasis located in southern Xinjiang Province, China. The optimization result shows that only 68 pumping wells are required for irrigation purposes of Qira Oasis, and this layout reduces the number of current pumping wells strongly by 59%. Thus, a large number of pumping wells can be closed to save resources. The optimizing method presented in this research is mathematically general and can be applied to other oasis areas without any obstacles. This method can provide decision-makers and managers with key information to ensure the optimal management and safety of valuable groundwater resources in oases.

1. Introduction

The management of irrigation water is of critical importance for both agricultural manufacturing and water safety, especially in arid regions with limited water resources [1]. China is the largest consumer of the world’s limited water resources, and the proportion of agricultural water in its total water consumption is greater than 70% [2,3]. Groundwater is an essential source for agricultural irrigation, accounting for nearly 10% in China [4,5]. In the vast land of Northwest China, there are numerous large and small oases. The irrigation water of oases mainly derived from surface runoff and groundwater. Compared with other regions in China, agriculture, as the main industry of oases, has a much higher demand for groundwater than that in other regions [6]. Hence, groundwater management, especially oases groundwater management, is an increasingly important issue. To efficiently exploit groundwater, the layout of pumping wells needs to be optimized to achieve the sustainable utilization of groundwater in irrigated areas [7,8,9]. For the overall planning of irrigation wells, determining suitable locations for wells is very important for helping prevent the construction of redundant wells or identifying existing wells that can be closed. Moreover, it is particularly important to evaluate the spatial distribution of wells in areas where water resources are scarce because the groundwater level in such areas is declining significantly and the ecological environment is deteriorating.
The current optimization research in the context of water supply management focuses on the spatial layout of pumping wells, while considering their number and rates [10,11]. In most of the earlier well layout studies, optimization model was built by using the number of pumping wells as the decision variable and using the economic benefit as the objective function; furthermore, the ecological benefits of irrigation water and groundwater depth are set as constraint conditions for the model [12,13]. Linear programming was first used to solve the optimization problems in previous studies, for example, Gorelick [14] considered many factors, such as the water yield, drawdown, and hydraulic gradient, to establish a model solved by linear and quadratic programming. However, the complexity of groundwater problems leads to complex groundwater management models, hence, nonlinear programing method gained popularity for the problems of the well number optimization [15,16,17]. Moreover, the results calculated by the traditional solver of nonlinear programming do not satisfy the global optimal conditions, and this approach requires too much effort in the calculation processes [18]. Therefore, the heuristic optimization algorithms, which demonstrated its feasibility and efficiency in the optimization problems, became the preferred methodology. The genetic algorithm is one example of these algorithms [19,20]. Other examples include the simulated annealing algorithm derived from the principle of solid annealing [21], and the particle swarm optimization based on bird predation behavior [22] (Table 1).
However, above studies only studied the optimized number of pumping wells, few methods have been developed to study their spatial layout. The choice of well location not only affects the input cost of the well construction, but also influences the economic benefits of agricultural production. For the overall planning of the wells in the irrigation area, especially in areas where water resources are scarce and the groundwater level is seriously degraded, it is necessary to determine the appropriate well locations to adjust the layout of existing wells. In general, researchers have considered a single factor, such as groundwater depth, and have divided the irrigation area into centralized and uniform well areas [25]. To obtain the appropriate layout of pumping wells, researchers have often combined groundwater simulation softwares and optimization algorithms to simulate the pumping-well domain and groundwater-level change [26]. Compared with these methods, a more appropriate evaluation method can be used to comprehensively assess the suitability of various spatial layouts of pumping wells from the perspective of the multiple factors that influence those layouts. The main evaluation methods currently available include the fuzzy comprehensive evaluation method [27], analytic hierarchy process method [28], empirical index method [29], and limit condition method [30], with related research concentrated mainly on soil, cultivated land, forests, city construction, and the benefits of resources on rural residential land [31]. For example, Liu et al. [32] selected land use and groundwater depth as the evaluation indices for spatial layout of pumping wells in irrigation areas, and determined the values of these indices by applying information entropy method. However, most of these methods only focused on the basic mechanisms, there is a lack of research from the perspective of water management in irrigation areas regarding the evaluation of the spatial layout of pumping wells, especially for oases in arid areas.
The southern margin of the Tarim basin is a preferable area for oases in Xinjiang and even in China as a whole [33]. The source of water for agricultural irrigation—the main industry of the oases in this region—originates mainly from the surface runoff caused by the melting of ice and snow in the nearby mountains. However, because of the limited surface water resources, groundwater resources have become increasingly important for the water supply in these oases. In recent years, newly cultivated oasis land has been planted with high-economic-value crops, such as walnut, red date, pomegranate, and other fruits. In the meantime, traditional crops, such as wheat, corn and cotton, are being gradually replaced by fruit cultivation. This change in planting structure has increased the demand for water, which has increased the utilization of groundwater and directly led to the construction of more pumping wells. However, the locations of these wells are randomly distributed. Optimizing the layout and number of wells has not been considered comprehensively from the perspective of the economy and the sustainable use of groundwater. Therefore, a large portion of water resources has been wasted and has caused issues related to groundwater safety and other problems.
Long-term monitoring data are essential for modeling. Because of the availability of long-term monitoring data (groundwater depth) provided by the Cele National Field Scientific Observation and Research Station (hereafter simply referred to as “Cele Station”), Qira Oasis—in the central part of the southern Tarim basin—was selected as an appropriate study area. A nonlinear programming model that considers economic and environmental factors was established first, and then, the optimal number of pumping wells was obtained using a genetic algorithm to solve the nonlinear programming model. After obtaining the number of wells, we introduced information entropy weighting to obtain a pumping-well spatial-layout assessment model. Finally, a “comprehensive evaluation value” of each pumping well was calculated, and wells were ranked in descending order. The optimal distribution of pumping wells was then presented based on the GIS data.

2. Methods

This section describes the methodologies used to optimize the number and layout of the pumping wells, and the global flowsheet is shown in Figure 1. First, we introduce the optimization method of the well number in Section 2.1, and then the spatial layout optimization method of the well is described in Section 2.2.

2.1. The Number of Wells Optimization

2.1.1. Nonlinear Programming Methodology

This study aims to optimize the number of pumping wells and the pumping rate in the study area to achieve the stated economic and ecological goals. The minimization of total cost is taken as the objective function. The number of wells (N) in the study area and the flow rate of the wells (Q) are considered as the decision variables. The flow rate is set to be equal to the annual average value. The supply and demand of water, drawdown, groundwater depth, number of wells, and flow rate are regarded as the constraining conditions. Hence, the model of nonlinear planning for the number of wells and the flow rate can be expressed as follows:
min R = R m + R d + R g + R e
where Rm is the annual management and maintenance cost (RMB) ( R m = a C N , in which a is the coefficient of well management and maintenance cost, C is the well construction fee (RMB), and N is the number of wells in the study area); Rd is the depreciation cost designed to recognize quantify the loss of services values that an asset has suffered in any year ( R d = C N / k , in which k is the fixed number of years of depreciation); Rg is the groundwater cost per year (RMB) ( R g = N Q T t r g , in which Q is the well flow rate (m3/h), T is the annual average number of operating days, t is the daily average number of operating hours for each well (h/d), and rg is the price per cubic meter of groundwater (RMB/m3); and Re is the energy cost per year (RMB) ( R e = 6 . 04 y Q T t h r e [34], in which h is the net pumping lift (m), and r e is the price per kWh of electricity (RMB/(kw·h))).
Subject to:
(1) Constraint of water resources:
N Q T t P a r P s w
N Q T t P a s
where P a r is the irrigation requirement per year in the study area (m3/a) ( P a r = i = 1 m I i A i , in which I i is the irrigation water of crop i throughout the growth period (m3/ha·a), A i is the irrigation area of crop i in the study area (ha)), P s w are available surface water resources used for irrigation (m3/a), which mainly includes river water and rainwater, and P a s refers to the water resources for irrigation (m3/a) ( P a s = P s + P r w P i r P d r P e r , in which P s is the extractable amount of groundwater (m3/a), P r w is the return water demand of irrigation and ecology (m3/a), P i r is the industrial water demand (m3/a), P d r is the domestic water demand (m3/a), and P e r is the required amount of ecological water, which is the total amount of water to maintain the ecological environment (m3/a)).
(2) Constraint of drawdown:
s 0 s max
where S 0 is the drawdown (m), and S max is the allowed maximum drawdown (m).
(3) Constraint of groundwater depth:
h + Δ h h min
h + Δ h h max
where Δ h refers to the dynamic change of the groundwater depth from the start moment to the end moment (m); h refers to the groundwater depth at the initial time (m); h min refers to the allowed minimum groundwater depth (m); and h max refers to the allowed maximum groundwater depth (m).
(4) Constraint of the number of wells:
0 < N < N e
where N e is the number of existing wells in the study area.
(5) Constraint of the well flow rate:
0 < Q < Q max
where Q max is the maximum value of the flow rate of individual wells.

2.1.2. Nonlinear Programming Application

In this study, the pumping-well number optimization model of all equations are represented as follows:
Minimum total consumption:
min R = R m + R d + R g + R e
Subject to the following constraints:
(1) Water resources:
N Q T t P a r P s w
N Q T t P a s
(2) Drawdown:
In Northwest China, one of the significant irrigation water sources is the groundwater in the phreatic water zone. The stability of groundwater resources in a short period can be determined based on the balance between groundwater exploitation and recharge.
After a long period of constant pumping rate, a relatively stable con of depression is formed near the well. Based on the Dupuit assumption, which is commonly used in the modeling of groundwater, we assumed that the phreatic water to the well is approximately horizontal and the discharge of groundwater in different cross-sectional area is equal to the flow rate of the well. The equations for the relationship between the flow-rate and drawdown of pumping-wells can be written as follows [5,35,36,37,38], and the detailed derivation process of Equation (12) is demonstrated in the Appendix A:
Q = 1.364 K ( 2 H s 0 ) s 0 / ( lg R / r 0 )
Q = Q ( 1 α )
Based on Equations (12) and (13), the relationship between s 0 and Q can be denoted as:
s 0 = H H 2 Q lg ( R / r 0 ) / ( 1.364 K ( 1 α ) )
where Q is the well flow-rate without interference from other wells (m3/h); K is the permeability coefficient (m/h); H is the distance from the bottom elevation of the unconfined aquifer to the phreatic free surface (m); R is the radius of influence of the pumping well (m); r 0 is the radius of the pumping well (m); and α is the abatement coefficient of total discharge.
Thus, the drawdown in this study can be demonstrated as:
H H 2 Q lg ( R / r 0 ) / ( 1.364 K ( 1 α ) ) s max
(3) Groundwater depth:
Because of the overexploitation of groundwater in the most northwestern areas of China, the constraints are too strict for the groundwater depth. In future planning, the quantity of groundwater exploitation should be lower than the value of groundwater recharge to achieve a gradual rise in the groundwater level. Thus, the maximum value of the difference in groundwater depth between the end moment and start moment must be smaller than zero for several decades to meet the constraint. Thus, Equations (5) and (6) can be changed to the constraint
Δ h 0
The balance of groundwater can be presented as
P d P r = μ Δ h F a
where P d is the groundwater discharge (m3) ( P d = P l d + P e d + N Q T t + P i r + P d r + P e r , in which P l d is the lateral discharge (m3) and P e d is the discharge from evaporation (m3)); and P r is the groundwater recharge (m3) ( P r = P l r + P p r + P r w , in which P l r is the lateral recharge (m3) and P p r is the recharge due to precipitation (m3)); μ is the extent of groundwater recharge; and F a is the area of the irrigated district (m2).
The above equation can be converted to
Δ h = ( P d P r ) / μ F a
Thus, the constraining condition of groundwater depth can be expressed as:
( P d P r ) / μ F a 0
(4) Number of wells:
0 < N < N e
(5) Well flow-rate:
0 < Q < Q max

2.1.3. Genetic Algorithm

This is a nonlinear programming problem with multivariate and multi-constraint conditions, and cannot be solved using the standard mathematical methods. To solve such problems, it can be converted into a linear problem, or dynamic programming, or the application of adaptive search algorithm, such as genetic algorithm, etc. Due to the robustness of genetic algorithm and the speed of computer operations, genetic algorithms are often the preferred algorithms for nonlinear optimization problems. The genetic algorithms take the objective function as the fitness function to carry out the optimization calculation until it satisfies the convergence condition or reaches evolutionary generation.
The conceptual framework of a genetic algorithm was originally developed by Holland [20]. It is a highly parallel, stochastic and adaptive search method that implements natural selection and the natural genetic mechanism of biology for reference. Different from traditional optimization algorithms, which propagate a single solution to the optimum one, the genetic algorithms search for optimum solutions by evaluating and successively improving a set of problem solutions, also denoted as generation. Each individual in the generation of solution is independently evaluated, which allows for the parallelization process.
Figure 2 displays the flow diagram of a genetic algorithm. The objective function can be f ( y ) ; n is the number of decision variables; y k is one of the decision variables, which contains a vector of y ; and the equation of y = [ y 1 , y 2 , , y k , , y n ] T is generated, in which T represents the transpose operator; the number of constraint conditions about inequality and equality are d1 and d2, respectively; g i ( y ) and h j ( y ) are the inequality and equality constraint conditions (i = 1, 2, …, d1; j = 1, 2, …, d2), respectively; and y k , l and y k , u are the lower and upper bounds of y k (k = 1, 2, …, n), respectively. In this study, f ( y ) represents minimum total consumption, y 1 , y 2 respectively represent the number of wells (N) in the study area and the flow rate of the wells (Q); there is no equality constraints, so g i ( y ) respectively represent the constraints water resources, drawdown and groundwater depth for i = 1, 2, 3; corresponding, y k , l and y k , u are the lower and upper bounds of N,Q (k = 1, 2), respectively. After defining the optimization model, the operation steps of the algorithm are as follows:
(1) Encoding the target problem
The internal expression (genotype) of chromosome is a genes combination, which determines the external expression of individuals. Binary coding, as one of the encoding works to map from phenotype to genotype, is implemented. Based on the upper bounds (200 m3/h) and lower bounds (0 m3/h) of the independent variable of the wells flow rate, we determined the code length of the binary code ( L = log 2 ( ( 200 0 ) / 10 2 ) 11 , the accuracy level is 10−2). Then, a set of binary codes were randomly generated, such as b0, b1, …, b10, which were converted to decimal numbers ( ( b 0 b 1 b 10 ) 2 = ( i = 0 10 b i 2 i ) 10 = Q t ), and the numbers were mapped to real numbers in the interval ( Q = Q t 200 0 2 11 1 ).
(2) Population initialization
To ensure the feasibility and diversity of the individuals in the first generation, the initial individuals are randomly produced according to the constraint conditions.
(3) Determination of fitness
The fitness function is transformed by the object function, and the individuals are evaluated in terms of its fitness value. In the search process, on the basis of the fitness function, the genetic algorithm utilizes the fitness value of each individual in the generation; however, it does not use external information.
(4) Genetic manipulation
The evolution of the population depends on the genetic operators acting on the current population and producing a new generation of the population. The common genetic operators encompass selection, crossover and mutation and play a decisive role in the algorithmic performance. Implicit parallelism and effective utilization of global information are two salient features of genetic algorithms.
(5) Optimal saving strategy
The idea behind an optimal preservation strategy is to retain the good traits of individuals in the parent generation in those of the offspring. In other words, the individual with the greatest fitness in the current population is used to replace the individual with the least fitness after crossover, mutation and other genetic manipulations in the population, and it does not carry through the crossover and mutation processes.
Through the above schemes, the optimal solution set can be obtained.

2.2. Well Layout Optimization

The well spatial layout can be optimized using the information entropy method and comprehensive evaluation index method.
Based on the information entropy weighting method, the spatial-layout evaluation model of pumping wells is established [40]. Groundwater resources must be conserved in woodland and grassland, and thus, these areas are unsuitable for groundwater exploitation. Likewise, it is improper to drill wells for irrigation on unused land or in residential areas. On cultivated land, some factors are selected for use as an evaluation index, including the type of land use, the depth of groundwater, the flow rate of each well, the density of a well (reflective of the productivity of the land) and the exploitation condition of groundwater.
Assuming W i j is the assessment matrix W of the assessment index j with different i of the pumping well spatial-layout, where i = 1, 2, …, 166 represents the number of each pumping well, j = 1, 2, 3, 4, respectively represent types of land use, groundwater depth, well density, and flow rate per well. The weights of the evaluation indices can be calculated by information entropy weighting. We can generalize the procedure as follows:
(1)
Normalization of the assessment matrix to resolve the multidimensional indices and multiple data units for different indices:
{ W i j = W i j min 1 i 166 W i j max 1 i 166 W i j min 1 i 166 W i j Benefit   type   index W i j = max 1 i 166 W i j W i j max 1 i 166 W i j min 1 i 166 W i j Cost   type   index
where W i j is the normalized value of the original evaluation matrix, 0 W i j 1 ; and in the evaluation index j, min 1 i 166 W i j is the minimum value and max 1 i 166 W i j is the maximum value. For i = 1, 2, …, 166, y i 1 represents the value of the land type indicator for each well, as measured by unit area demand for groundwater resources shown in Table 6, which is a benefit type index. The maximum and minimum values of y i 1 are 1504.71 and 394.46, respectively; y i 2 represents the value of the groundwater depth indicator for each well, as measured by the mean of the groundwater depth interval in the distribution area of the well shown in Figure 4, which is a cost type index. For example, if the well is distributed in an area where the groundwater depth is in the interval [0, 10], then y i 2 = 5 . The maximum and minimum values of y i 2 are 55 and 10, respectively; y i 3 represents the value of the well density indicator for each well, as measured by the theoretical well distance, which is a benefit type index. The maximum and minimum values of y i 3 are 589.38 and 294.69, respectively; y i 4 represents the value of the well density indicator for each well, which is a cost type index. The maximum and minimum values of y i 4 are 2000 and 30, respectively. A larger “Benefit type index” value indicates more expensive benefits, while the “Cost type index” indicates the opposite.
(2)
Calculation of the entropy value v j of index j:
v j = k i = 1 166 q i j ln q i j
where q i j is the proportion of the evaluation index j within different values of i, q i j > 0 ( q i j = ( W i j + 1 ) / i = 1 m ( W i j + 1 ) ), k = 1 / ln 166 , and 0 v j 1 [41].
(3)
Calculation of the deviation coefficient c j of index j:
c j = 1 v j
(4)
Calculation of the weight ω j of index j:
ω j = c j / j = 1 4 c j
Based on Equations (22)–(25), the index weights (types of land use, groundwater depth, well density, and flow rate per well) of the spatially distributed evaluation model of the pumping wells in the study area can be estimated (Table 2).
Together with the weights of the assessment index, the “comprehensive evaluation index” e i [32,40,42], regarding the rationality of the pumping well spatial distribution, will be estimated via
e i = j = 1 4 ω j q i j
Then, the rationality of the spatial distribution is sorted based on the descending order of e i . Ultimately, the optimal number of wells is combined with the assessment results of the pumping well spatial distribution to close the lower e i of the wells, and according to the GIS data, the optimized distribution of wells can be obtained.

3. Case Study

3.1. Study Area

This research is based on the data collected at the Qira Oasis, Hotan city, southern Xinjiang Province, China (Figure 3). Irrigation plays a crucial role in local agriculture, as the area has high evaporation and low precipitation in the area [33]. In the past few years, the imbalance between the supply and demand of water resources in this area has become increasingly prominent due to rapid population and economic growth. This shortage of the water supply has resulted in a serious overexploitation of groundwater by local residents for their agricultural water usage, which has led to a groundwater recession and a worsening ecological environment. Optimizing the pumping well layout is necessary to protect the groundwater supply and environment in the Qira Oasis.
The climate of Qira Oasis is temperate and arid, with an average annual precipitation of 35.1 mm and an average annual evaporation of 2595.3 mm [43]. The Qira River, which originates from the Kunlun Mountains, is the only exogenous supply of water to the Qira Oasis, and the river has an annual mean runoff of 1.27 × 108 m3. The main crops grown in the Qira Oasis are wheat and maize, but there are also many plantations of fruit crops, such as pomegranate, red date and walnut. The surface runoff of the Qira River is the main irrigation source for the agriculture of the oasis; however, since the water demands have increased due to agricultural expansion, the groundwater has become an important irrigation supplement in the oasis. The Qira Oasis covers an area of 274.63 km2, within which there are 166 wells (Figure 4) with an average flow rate of 135 m3/h.

3.2. Data Description

The data required for this research were obtained from the Xinjiang Statistical Yearbook, Qira National Field Scientific Observation and Research Station, Hotan Hydrology Bureau, and Qira Hydrological Station (initial data not displayed due to their size). Table 3 shows the irrigated area and net irrigation water demand of the main crops.
Table 4 displays the amount of water, including the supply and demand of all water resources and the specific quantities of groundwater resources’ recharge and discharge. The data in Table 4 was estimated according to crop water demand and crop planting ratio, the relevant industrial data on the Internet, and the water demand of humans and animals. The hydrogeological parameters are shown in Table 5, and these parameters were obtained from the reference of Liu et al. (2015). These data are regionally independent.

4. Results

4.1. Well-Number Optimization

In the Qira Oasis, 166 pumping wells operate normally. With economic growth, the maintenance cost of pumping wells is increasing. To solve these problems, we use Equations (9)−(11), (15) and (19)−(21) as the final computational model, select the data in Table 3, Table 4 and Table 5 as input values, and then perform the calculations using the genetic algorithm. Due to the randomness of the genetic algorithm, it is necessary to run the model multiple times (in this study, 100 times) and take the average value. The optimal well flow rate is 146 m3/h. Although this value is higher than the original flow rate (135 m3/h), there is no threat to groundwater safety. The most appropriate number of wells is calculated to be 68, representing a reduction of 59.04%, which is significantly less than the number of existing wells (166). The cost after optimization is 1.77 × 108 RMB, which is a reduction of 54.38% compared with the original cost (3.88 × 108 RMB). The optimization process results can lead to a large decrease in the economic cost.

4.2. Well Spatial-Layout Optimization

The demand for groundwater resources in the unit area changes with different types of land usages, and decision-makers can make different plans for different wells according to variant groundwater demands in different regions. The demand for groundwater resources in town and village areas is equal to the total water requirement for domestic water and industrial water. In woodlands, such as those established to farm red date, walnut, pomegranate and other fruits, as well as the surrounding ecological woodland, the demand for groundwater resources is the water demand of the fruit trees and the surrounding ecological woodland minus the amount of surface water resources. On the dryland plain areas in the study region, land use is mainly a variety of crops and the surrounding protective forests, and the demand for groundwater resources here is equal to the amount of water needed for crops and the surrounding forests minus the amount of surface water resources. After dividing the demand for groundwater resources by the corresponding area for each type of land, the unit area demand for groundwater resources (Table 6) in the town and village areas, woodland and dryland plains is 394.46, 1504.71, and 654.22 mm, respectively.
Groundwater depth is an important factor used to measure the carrying capacity of groundwater resources in an irrigated area. The decline in the groundwater level caused by extraction in the research area has caused serious ecological and environmental problems, such as desertification and shrub/forest loss. Therefore, when selecting locations for pumping wells, we should choose an area where the groundwater depth is very shallow and easily extracted, and the area should be far from areas where groundwater is deep and water resources are exploited excessively.
Due to the kriging interpolation method is more suitable for geoscience, we select the ordinary kriging method to interpolate the groundwater depth in the Qira Oasis. As shown in Figure 5, the groundwater depth of the Qira Oasis shows an increasing trend from the North to the South.
Due to a lack of numerical data on well density, the well density is determined based on a comparison of the length of theoretical and practical well distance. Thus, the calculation of the number of wells by identifying the theoretical distance within alternative wells can be applied to identify which index best assesses the well density. If the number of wells is large, the wells in the study area are densely distributed. According to the distribution shape of the well, it can be divided into plum-shaped layout and square layout. For example, the wells numbered 8, 9, 11, and 12 are a plum-shaped layout, the wells numbered 41, 42, 43, and 48 are a square layout as shown in Figure 4. On the basis of the GB/T 50625-2010 [34], the theoretical well distance can be framed as follows:
L = { 99.9 f a   square   layout 107.7 f a plum-shaped   layout
The theoretical irrigation area of each well is expressed as
f a = Q T t b
where b is the net irrigation quota (m3/ha).
The weights of the evaluation indices, which have an effect on the pumping well spatial layout in the Qira Oasis Irrigation District, can be calculated using the entropy weighting information. The well density had the highest weight (0.5854), while the weights of the type of land-use and flow rate were relatively lower, with values of 0.1086 and 0.0746, respectively. The 12 wells distributed in the shelter forest belt are retained. In the Qira Oasis, the values of e i ranged from 0.000252 to 0.196688. The alternative wells with sequence numbers prior to 57 are retained combined with the optimal number of wells.
Thus, as shown in Figure 6, wells located in areas where groundwater is overexploited will be closed due to a variety of problems, such as low flow rate or inefficient operation. The optimal distributed pattern of wells is relatively concentrated; few of the wells are located in the area of the groundwater depth lower than 14.53 m. Additionally, the spatial distribution of wells is more conducive to the sustainable utilization of groundwater than was the original distribution.

5. Discussion and Conclusions

In this study, nonlinear programming model solved by genetic algorithm and comprehensive index methods, combined with information entropy weighting, are utilized to optimize the wells spatial layout in the Qira oasis, northwest China. The application of this new optimization method in the Qira Oasis demonstrates the significant reduction in the number of pumping wells and the optimization of the spatial layout considering the water demands and pumping rates.
In general, our study provides an approach to optimize the spatial layout of pumping wells. This method is mathematically general and can be further applied to other oasis regions where groundwater is overexploited and the well density is high. Through this method, the number of excess wells can be greatly reduced, and a substantial reduction in the expense of well construction and operational management can be achieved. Furthermore, this study can provide instruction for residents living in oases to utilize groundwater more rationally, which is of great significance to the sustainable utilization of groundwater and helps improve the ecological environments of these areas.
There are some shortcomings of this work. One weakness is that a simplification (Dupuit assumption) and analytical solution were used in this research for the groundwater flow calculation. A more accurate simulation of groundwater flow can be achieved by establishing a numerical groundwater model, such as MODFLOW, GMS. However, the numerical modeling of groundwater movement in this type of large-scale oasis requires excessive computational power and the accurate simulation of groundwater flow is not the focus of this research. More methods, such as particle swarm optimization, simulated annealing algorithm, can also be applied to solve the nonlinear programming model to ensure the optimization results. Additionally, the comprehensive index method is a traditional ranking approach, thus, new and more advanced methods can be adopted to solve our problem, for example, neural network method. However, the comparison of different methods and algorithms is not the focus of this research. In the future work, we will implement multiple optimization algorithms and quantify the uncertainties further to search for the best optimization model and algorithm for the well layout problems.

Author Contributions

Conceptualization, D.G. and H.D.; methodology, J.Z. and X.F.; formal analysis, M.S.; investigation, H.D. and X.Z.; resources, B.H. and D.G.; data curation, Y.L.; writing—original draft preparation, Y.L. and M.S.; writing—review and editing, H.D., J.J. and S.S.; supervision, B.H. and D.G.; project administration, D.G.

Funding

This work was financially supported by the Regional collaborative innovation project for Xinjiang Uygur Autonomous Region (Shanghai cooperation organization science and technology partnership project) (2017E01029), Western Light—western region leading scientists supporting project, Chinese Academy of Sciences (2018-XBYJRC-002), and Chinese National Natural Science Fund (U1603343, 41501107, 41807182, 41601595).

Acknowledgments

We thank Jie Xue and Yunfei Liu for their help.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The Calculation of Well Flow Rate Based on Dupuit Assumption

The wells in this study are located in unconfined aquifer area. After a long period of constant pumping rate, a relatively stable drawdown funnel is formed near the well. For practical purpose, the Dupuit Assumption has been applied into the phreatic water wells [35,36,37,38].
The assumptions are presented as follows:
The groundwater movement can be considered stable in a short period due to the balance of groundwater exploitation and recharge.
(1)
The aquifer is homogeneous and isotropic.
(2)
The groundwater obeys Darcy Law.
(3)
The phreatic water to wells can be considered as horizontal.
(4)
The flow rate of the wells equals to the discharge of groundwater.
According to the assumptions, the steady groundwater governing equation can be denoted by the following partial differential equations:
d d r ( r d h 2 d r ) = 0
With boundary conditions
h | r = R = H h | r = r w = H w
The Equation (A1) can be integrated as:
r d ( h 2 ) d r = C
Given the Darcy Law,
Q = K A J = K · 2 π r · h · d h d r = π r K d ( h 2 ) d r
Based on Equations (A3) and (A4),
r d ( h 2 ) d r = Q π K
Based on the boundary conditions (Equation (A2)), Equation (A5) can be integrated as,
H w H d ( h 2 ) = Q π K r w R d r r H 2 H w 2 = Q π K ln R r w
Because of s w = H H w , Equation (A6) can be represented as,
( 2 H s w ) s w = Q π K ln R r w Q = 1.364 K ( 2 H s w ) s w lg R r w
which is Equation (12) in the text.

References

  1. Storm, H.; Heckelei, T.; Heidecke, C. Estimating irrigation water demand in the Moroccan Drâa Valley using contingent valuation. J. Environ. Manag. 2011, 92, 2803–2809. [Google Scholar] [CrossRef]
  2. García-Garizábal, I.; Causapé, J.; Abrahao, R. Application of the Irrigation Land Environmental Evaluation Tool for flood irrigation management and evaluation of water use. CATENA 2011, 87, 260–267. [Google Scholar] [CrossRef]
  3. Rivers, M.; Coles, N.; Zia, H.; Harris, N.R.; Yates, R.; Harvey, H.W. How could sensor networks help with agricultural water management issues? Optimizing irrigation scheduling through networked soil-moisture sensors. In Proceedings of the 2015 IEEE Sensors Applications Symposium, Zadar, Croatia, 13–15 April 2015. [Google Scholar]
  4. Chu, Y.; Shen, Y.; Yuan, Z. Water footprint of crop production for different crop structures in the Hebei southern plain, North China. Hydrol. Earth Syst. Sci. 2017, 21, 3061–3069. [Google Scholar] [CrossRef] [Green Version]
  5. Xin, L.; Wang, S.; Huo, Z.; Li, F.; Hao, X. Optimizing layout of pumping well in irrigation district for groundwater sustainable use in northwest China. Hydrol. Process. 2015, 29, 4188–4198. [Google Scholar]
  6. Ling, H.; Xu, H.; Shi, W.; Zhang, Q. Regional climate change and its effects on the runoff of Manas River, Xinjiang, China. Environ. Earth Sci. 2011, 64, 2203–2213. [Google Scholar] [CrossRef]
  7. Bekri, E.; Disse, M.; Yannopoulos, P. Optimizing Water Allocation under Uncertain System Conditions for Water and Agriculture Future Scenarios in Alfeios River Basin (Greece)—Part B: Fuzzy-Boundary Intervals Combined with Multi-Stage Stochastic Programming Model. Water 2015, 7, 6427–6466. [Google Scholar] [CrossRef] [Green Version]
  8. Campana, P.E.; Olsson, A.; Li, H.; Yan, J. An economic analysis of photovoltaic water pumping irrigation systems. Int. J. Green Energy 2016, 13, 831–839. [Google Scholar] [CrossRef]
  9. Ketabchi, H.; Ataie-Ashtiani, B. Review: Coastal groundwater optimization—Advances, challenges, and practical solutions. Hydrogeol. J. 2015, 23, 1129–1154. [Google Scholar] [CrossRef]
  10. Gaur, S.; Chahar, B.R.; Graillot, D. Analytic elements method and particle swarm optimization based simulation–optimization model for groundwater management. J. Hydrol. 2011, 402, 217–227. [Google Scholar] [CrossRef]
  11. Singh, R.M.; Datta, B. Identification of groundwater pollution sources using GA-based linked simulation optimization model. J. Hydrol. Eng. 2006, 11, 101–110. [Google Scholar] [CrossRef]
  12. Ayvaz, M.T.; Karahan, H. A simulation/optimization model for the identification of unknown groundwater well locations and pumping rates. J. Hydrol. 2008, 357, 76–92. [Google Scholar] [CrossRef]
  13. Sidiropoulos, E.; Tolikas, P. Well Locations and Constraint Handling in Groundwater Pumping Cost Minimization via Genetic Algorithms. Water Air Soil Pollut. Focus 2004, 4, 227–239. [Google Scholar] [CrossRef]
  14. Gorelick, S. A review of distributed parameter groundwater management modeling method. Water Resour. Res. 1983, 19, 305–319. [Google Scholar] [CrossRef]
  15. Bouwer, H. Integrated water management: Emerging issues and challenges. Agric. Water Manag. 2000, 45, 217–228. [Google Scholar] [CrossRef]
  16. Monaghan, R.M.; Wilcock, R.J.; Smith, L.C.; Tikkisetty, B.; Thorrold, B.S.; Costall, D. Linkages between land management activities and water quality in an intensively farmed catchment in southern New Zealand. Agric. Ecosyst. Environ. 2007, 118, 211–222. [Google Scholar] [CrossRef]
  17. Stanko, Z.P.; Boyce, S.E.; Yeh, W.G. Nonlinear model reduction of unconfined groundwater flow using POD and DEIM. Adv. Water Resour. 2016, 97, 130–143. [Google Scholar] [CrossRef] [Green Version]
  18. Aral, M.M.; Guan, J. Genetic Algorithms in Search of Groundwater Pollution Sources. In Advances in Groundwater Pollution Control and Remediation; Springer: Dordrecht, The Netherlands, 1996. [Google Scholar]
  19. Goldberg, D.E. Genetic Algorithms in Search, Optimization and Machine Learning; Addison-Wesley Longman Publishing Co., Inc.: Boston, MA, USA, 1989. [Google Scholar]
  20. Holland, J.H. Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence; U Michigan Press: Oxford, UK, 1975; Volume viii, p. 183. [Google Scholar]
  21. Kirkpatrick, S.; Vecchi, M.P. Optimization by Simulated Annealing. Science 1983, 220, 671–680. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Kennedy, J. Particle Swarm Optimization. In Proceedings of the 1995 IEEE International Conference on Neural Networks, Perth, Australia, 27 November–1 December 2011; Volume 1944, pp. 1942–1948. [Google Scholar]
  23. Huang, C.; Mayer, A.S. Pump-and-treat optimization using well locations and pumping rates as decision variables. Water Resour. Res. 1997, 33, 1001–1012. [Google Scholar] [CrossRef] [Green Version]
  24. Fragoso, T.; da Conceição Cunha, M.; Lobo-Ferreira, J.O.P. Optimal pumping from Palmela water supply wells (Portugal) using simulated annealing. Hydrogeol. J. 2009, 17, 1935–1948. [Google Scholar] [CrossRef]
  25. Milnes, E. Process-based groundwater salinisation risk assessment methodology: Application to the Akrotiri aquifer (Southern Cyprus). J. Hydrol. 2011, 399, 29–47. [Google Scholar] [CrossRef] [Green Version]
  26. Bayer, P.; Duran, E.; Baumann, R.; Finkel, M. Optimized groundwater drawdown in a subsiding urban mining area. J. Hydrol. 2009, 365, 95–104. [Google Scholar] [CrossRef]
  27. Tang, H.; Huang, B.G.; Dao-Xuan, L.I.; Zeng, F.H. Determing water driving effectiveness of oil reservoir with fuzzy comprehensive evaluation method. Pet. Explor. Dev. 2002, 29, 97–99. [Google Scholar]
  28. Wang, H.F.; Dong, Z.C. Application of Analytic Hierarchy Process(AHP)Method in Regional Water Resources Sustainable Development Evaluation. Water Power 2002, 7, 2–14. [Google Scholar]
  29. Pan, X.A. Study on Village Reclamation Suitability Evaluation Based on the Index Method. Mod. Agric. Sci. Technol. 2013, 2, 340–342. [Google Scholar]
  30. Ziadat, F.M. Land suitability classification using different sources of information: Soil maps and predicted soil attributes in Jordan. Geoderma 2007, 140, 73–80. [Google Scholar] [CrossRef]
  31. Keshavarzi, A.; Sarmadian, F.; Ahmadi, A. Spatially-based model of land suitability analysis using block kriging. Aust. J. Crop Sci. 2011, 5, 1533–1541. [Google Scholar]
  32. Liu, X.; Wang, S.; Hao, X. Suitability assessment on spatial layout of pumping wells in Hongyashan irrigation district. Trans. Chin. Soc. Agric. Eng. 2013, 29, 101–109. [Google Scholar]
  33. Jie, X.; Gui, D.; Lei, J.; Zeng, F.; Mao, D.; Zhang, Z. Model development of a participatory Bayesian network for coupling ecosystem services into integrated water resources management. J. Hydrol. 2017, 554, 50–65. [Google Scholar]
  34. Technical Code for Water Wells. PRC, N.S.o. GB/T 50625-2010. 2010. Available online: http://www.gbstandards.org/mobile/GB_List-m.asp?id=82800 (accessed on 18 August 2010).
  35. Ahlfeld, D.P.; Sawyer, C.S. Well Location in Capture Zone Design Using Simulation and Optimization Techniques. Groundwater 1990, 28, 507–512. [Google Scholar] [CrossRef]
  36. Guan, J.; Aral, M.M. Optimal remediation with well locations and pumping rates selected as continuous decision variables. J. Hydrol. 1999, 221, 20–42. [Google Scholar] [CrossRef]
  37. Lin, H.-T.; Tan, Y.-C.; Chen, C.-H.; Lin, W.-S.; Liu, C.-W. Multivariate approaches optimize locations of groundwater pumping facilities for different hydrogeological scales. Hydrol. Process. 2012, 26, 2985–2996. [Google Scholar] [CrossRef]
  38. Maas, C. Groundwater flow to a well in a layered porous medium: 1. Steady flow. Water Resour. Res. 1987, 23, 1675–1681. [Google Scholar] [CrossRef]
  39. Kanemoto, Y.; Sugawara, R.; Ohmura, M. A genetic algorithm for the rectilinear Steiner tree in 3-D VLSI layout design. In Proceedings of the 2004 47th Midwest Symposium on Circuits and Systems, Hiroshima, Japan, 25–28 July 2004. [Google Scholar]
  40. Li, X.; Wang, K.; Liu, L.; Xin, J.; Yang, H.; Gao, C. Application of the Entropy Weight and TOPSIS Method in Safety Evaluation of Coal Mines. Procedia Eng. 2011, 26, 2085–2091. [Google Scholar] [CrossRef] [Green Version]
  41. Zhang, H.; Gu, C.L.; Gu, L.W.; Zhang, Y. The evaluation of tourism destination competitiveness by TOPSIS & information entropy—A case in the Yangtze River Delta of China. Tour. Manag. 2011, 32, 443–451. [Google Scholar]
  42. Chen, R.J.; Qian, H.L.; Kan, H.D.; Yuan, D. A Review of Comprehensive Index Methods in Water Quality Assessment. J. Environ. Occup. Med. 2009, 26, 581–584. [Google Scholar]
  43. Qian, J.; Zhao, J.; Gui, D.; Feng, X. Effeccts of ecological water use and irrigation on groundwater depth in Cele Oasis. Bull. Soil Water Conserv. 2018, 38, 96–102. [Google Scholar]
Figure 1. Global flowsheet of well layout optimization.
Figure 1. Global flowsheet of well layout optimization.
Water 11 00970 g001
Figure 2. Genetic algorithm flowchart [39].
Figure 2. Genetic algorithm flowchart [39].
Water 11 00970 g002
Figure 3. Location of Qira Oasis.
Figure 3. Location of Qira Oasis.
Water 11 00970 g003
Figure 4. Distribution of existing wells in Qira Oasis.
Figure 4. Distribution of existing wells in Qira Oasis.
Water 11 00970 g004
Figure 5. Groundwater depth in the Qira Oasis.
Figure 5. Groundwater depth in the Qira Oasis.
Water 11 00970 g005
Figure 6. Distribution of wells after optimization.
Figure 6. Distribution of wells after optimization.
Water 11 00970 g006
Table 1. Chart of different heuristic optimization methods in the field of pumping wells optimization.
Table 1. Chart of different heuristic optimization methods in the field of pumping wells optimization.
Name of MethodReferenceExample of Application
Genetic Algorithm (GA)Holland, 1975Huang et al. [23] optimized the total cost and number of pumping wells applying GA method.
Simulated Annealing (SA)Kirkpatrick et al., 1983Fragoso et al. [24] regarded the pipe costs, the well installation, and pumping costs as the objective function applying SA method.
Particle Swarm Optimization (PSO)Kennedy and Eberhart, 1995Gaur et al. [10] made pumping water amount maximum and the total cost minimum using PSO method.
Table 2. Evaluation model of wells spatial layout.
Table 2. Evaluation model of wells spatial layout.
TargetEvaluation Indexes
Evaluation model of wells spatial layout1. Type of land use
2. Depth of groundwater
3. Flow rate of each well
4. Density of a well
Table 3. The amount of irrigated area and the required quantity of net irrigation water in the entire crop growing season in the Qira Oasis Irrigated District.
Table 3. The amount of irrigated area and the required quantity of net irrigation water in the entire crop growing season in the Qira Oasis Irrigated District.
CropIrrigated Area (ha)Net Irrigation Water Demand (m3/ha·a)Water Demand (m3/a)
Wheat133655877.46 × 106
Maize134546806.29 × 106
Pomegranate669.513,5509.07 × 106
Red date198710,6502.12 × 107
Walnut137110,6001.45 × 107
Sum of values6708.545,0675.85 × 107
Table 4. Values of the water resources.
Table 4. Values of the water resources.
ParameterSymbolUnitValue
Extractable amount of groundwaterPsm3/a1500.4 × 104
Return water of irrigation and ecologyPrw213.5 × 104
Industrial water demandPir104.3 × 104
Domestic water demandPdr172.84 × 104
Required amount of ecological waterPer438.6 × 104
Difference of lateral discharge and rechargePlr-Pld1500.4 × 104
Discharge from evaporationPed0
Recharge due to precipitationPpr0
Surface water resources for irrigationPsw 5815.5 × 104
Note: The groundwater depth in the irrigation district is usually below 5 m. Compared with other recharge and discharge processes, infiltration of precipitation and groundwater evaporation will be ignored.
Table 5. Observed hydrogeological parametric values.
Table 5. Observed hydrogeological parametric values.
ParameterSymbolUnitValue
The well management and maintenance cost coefficientaRMB23.72 × 104
The number of wellsn 20
The annual average operating daysTd70
The daily average operating hoursth/d12
Groundwater price per cubic meterrgRMB/m30.20
Net pumping lifthm22.5
Cost of electricity per K W · h re RMB / ( KW · h ) 0.15
Allowed maximum drawdownSmaxm5.00
Permeability coefficientKm/h0.72
Distance from the bottom elevation of the unconfined aquifer to the phreatic free surfaceHm55.00
Radius of influence of the pumping wellRm400
Radius of the pumping wellr0m0.175
abatement coefficient of total dischargeα 0.03
Maximum value of the well flow rateQmaxm3/h200
Existing number of wellsNe 166
Degree of groundwater rechargeμ 0.25
Table 6. The different types of land area and unit area demand for groundwater resources in the Qira Oasis.
Table 6. The different types of land area and unit area demand for groundwater resources in the Qira Oasis.
Type of Land UseArea (ha)Unit Area Demand for Groundwater Resources (mm)
Towns and villages176.5394.46
Woodland4027.51504.71
Dryland plain2681.0654.22

Share and Cite

MDPI and ACS Style

Liu, Y.; Shen, M.; Zhao, J.; Dai, H.; Gui, D.; Feng, X.; Ju, J.; Sang, S.; Zhang, X.; Hu, B. A New Optimization Method for the Layout of Pumping Wells in Oases: Application in the Qira Oasis, Northwest China. Water 2019, 11, 970. https://doi.org/10.3390/w11050970

AMA Style

Liu Y, Shen M, Zhao J, Dai H, Gui D, Feng X, Ju J, Sang S, Zhang X, Hu B. A New Optimization Method for the Layout of Pumping Wells in Oases: Application in the Qira Oasis, Northwest China. Water. 2019; 11(5):970. https://doi.org/10.3390/w11050970

Chicago/Turabian Style

Liu, Yi, Mengyang Shen, Jianping Zhao, Heng Dai, Dongwei Gui, Xinlong Feng, Jiali Ju, Shilei Sang, Xiaoying Zhang, and Bill Hu. 2019. "A New Optimization Method for the Layout of Pumping Wells in Oases: Application in the Qira Oasis, Northwest China" Water 11, no. 5: 970. https://doi.org/10.3390/w11050970

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