Optimization of Water and Energy Spatial Patterns in the Cascade Pump Station Irrigation District

: Cascade pump station irrigation districts (CPSIDs) consume large quantities of water and energy. Water- and energy-saving results and income increases are guaranteed under the sustainable development of the CPSID. The CPSID is divided into several sub-districts based on the elevation difference of topography and pump station distributions. The spatial patterns of crops and irrigation technologies can be changed by adjusting crop planting structures and developing drip irrigation in each sub-district. Its optimization will change the spatial patterns of irrigation water and energy consumption to achieve water- and energy-saving results, increase income, and provide an ecological advantage. To obtain the optimal spatial patterns of water and energy in the CPSID, a multiobjective linear programming model of minimum irrigation water consumption, minimum energy consumption, and highest crop output value was established. This model was applied to the Jingdian Phase I Irrigation District in northwest China, and an optimal scheme of water and energy spatial patterns was obtained. Compared with the present situation, the optimal scheme could save water by 26.18%, save energy by 29.38%, and increase income by 29.55%. The increased investment in the drip irrigation project would lead to reduced irrigation water and energy consumption and increased crop output value. The research results provide a scientiﬁc basis for the sustainable development of agriculture and ecological environment protection in the CPSID.


Introduction
Cascade pump stations irrigation district (CPSID) is divided into several sub-districts according to the topography elevation difference and the distribution of pump stations. The CPSID is a complex system consisting of cascade pump stations (CPSs), channels, pipes, and sub-districts, in which each part has a close water-energy nexus. The spatial patterns of water and energy are determined by the spatial patterns of irrigation technologies and crops in each sub-district. Thus, the research objective of this paper is to obtain the optimal spatial pattern of crops and irrigation technologies, in which energy and irrigation water use are minimal, but the crop value is the highest in the CPSID.
Due to population growth and economic development, there is an expected increased demand of 50% more food, 40% more energy, and 30% more water by 2030 [1]. Agriculture is the key link to ensure crop growth, but it consumes large quantities of water and energy. Therefore, much research has been carried out on the water-food-energy nexus in agricultural irrigation to achieve higher crop yields with less water and energy consumption [2][3][4][5][6][7]. solution. A KKM approach is the new algorithm for IO based on capacitated transportation problems in a neutrosophic environment [46]. In application problems, the decision-maker may often encounter a problem to optimize the ratio with both the denominator and numerator being linear. If the ratio is considered as an objective function, then under linear constraints the problem is said to be linear fractional programming (LFP). In recent years, some algorithms have been proposed to solve LFP problems [47][48][49].
This paper is organized as follows. In Section 2, The information of the study area is introduced. The multi-objective linear programming and calculation of water and energy spatial patterns in the CPSID are introduced in Section 3. In Section 4, we obtain the optimal scheme of water and energy spatial patterns of the Jingdian Phase I Irrigation District. Then, we analyze changes in various indexes of the optimal scheme under different drip irrigation project investment, and we further discuss the sustainable development of the CPSID. Section 5 presents the conclusion.

Study Area
Data for the Jingdian Phase I Irrigation District were obtained from the administrative department of the district, namely the Jingtaichuan Electric Lifting Irrigation Water Resources Utilization Center of Gansu Province.
The Jingdian Phase I Irrigation District is located in the central region of Gansu Province, in northwest China (103 • 20 -104 • 04 E and 37 • 26 -38 • 41 N). Figure 1 shows the location of the ID, with an elevation of 1561 to 1757 m. The average annual rainfall is 184 mm, and the average annual evaporation is 2290 mm. Therefore, the Jingdian Phase I Irrigation District is located in an extremely arid area, and irrigation is the primary crop water source.
The terrain elevation difference is very large in the Jingdian Phase I Irrigation District. According to the elevation difference and the distribution of CPS, the ID is divided into 10 sub-districts (Table 1). Water from the Yellow River is transported to the farmland through the CPS, water conveyance pipe and channel (WCP&C), as well as the water distribution pipe and channel system (WDP&CS) of each sub-district. The main crop irrigation techniques for the farmland in the ID are surface and drip irrigation. The Jingdian Phase I Irrigation District comprises a complex system ( Figure 2) of a water source (the Yellow River), 13 CPSs, 13 WCP&C segments, and 10 sub-districts. The WCP&C is composed of pressure pipes and open channels. Each sub-district contains farmland and WDP&CS; the WDP&CS consists of open channel systems and pipe networks. The water from the WCP&C is diverted to the farmland by the open channel system for surface irrigation. The water from the WCP&C is sent to the drip irrigation emitter near the crops through pipe networks.  -1580  72  94  357  2  1550  1580-1596  69  93  2306  3  4439  1596-1620  67  93  6533  4  1707  1620-1638  70  95  2375  5  600  1638-1658  69  94  806  6  3140  1658-1676  68  92  4054  7  2319  1676-1694  71  93  3454  8  1545  1694-1711  72  94  1943  9  31  1711-1738  74  95  35  10  73  1738-1757  73  The terrain elevation difference is very large in the Jingdian Phase I Irrigation District. According to the elevation difference and the distribution of CPS, the ID is divided into 10 sub-districts (Table 1). Water from the Yellow River is transported to the farmland through the CPS, water conveyance pipe and channel (WCP&C), as well as the water distribution pipe and channel system (WDP&CS) of each sub-district. The main crop irrigation techniques for the farmland in the ID are surface and drip irrigation. The Jingdian Phase I Irrigation District comprises a complex system ( Figure 2) of a water source (the Yellow River), 13 CPSs, 13 WCP&C segments, and 10 sub-districts. The WCP&C is composed of pressure pipes and open channels. Each sub-district contains farmland and WDP&CS; the WDP&CS consists of open channel systems and pipe networks. The water from the WCP&C is diverted to the farmland by the open channel system for surface irrigation. The water from the WCP&C is sent to the drip irrigation emitter near the crops through pipe networks.  3  4439  1596-1620  67  93  6533  4  1707  1620-1638  70  95  2375  5  600  1638-1658  69  94  806  6  3140  1658-1676  68  92  4054  7  2319  1676-1694  71  93  3454  8  1545  1694-1711  72  94  1943  9  31  1711-1738  74  95  35  10  73  1738-1757  73 95 89 The parameters of the pump stations are shown in Table 2. Table 1 shows the irrigation water use efficiency of different irrigation technologies in each sub-district. The water use efficiency of the WCP&C segment is shown in Table 3. The parameters of the irrigation technologies are shown in Table 4.  The parameters of the pump stations are shown in Table 2. Table 1 shows the irrigation water use efficiency of different irrigation technologies in each sub-district. The water use efficiency of the WCP&C segment is shown in Table 3. The parameters of the irrigation technologies are shown in Table 4. The spatial patterns of crops and irrigation technologies for the Jingdian Phase I Irrigation District in 2020 are listed in Table 5. The annual irrigation water consumption was 14,219 × 10 4 m 3 , the annual energy consumption was 18,770 × 10 4 kWh. The crop output value was 65,244 × 10 4 yuan, and the irrigated area was 21,952 ha. The drip irrigation area comprised 7.56% of the total irrigated area.    Table 4. To reduce the investment, two principles were used to optimize the spatial patterns of crops and irrigation technologies in the CPSID. The first principle was to keep the present situation of each pump station and WCP&C unchanged; the second principle was to keep the area of various crops in the present drip irrigation of each sub-district unchanged.

The Mathematical Model and Calculation
According to the type of irrigation technology, each sub-district is divided into several parts ( Figure 3). The present area of the sub-district consists of two parts-one is the present area of surface irrigation and the other is the present area of drip irrigation where d ik is the area of crop k in sub-district i, and d ik is a constant. In the optimization scheme, x ijk is taken as the decision variable, which is the area of crops k of types of irrigation technologies j (j = 1 for surface irrigation, j = 2 for drip irrigation) in sub-district i. Thus, in the optimization scheme, the sub-district comprises three parts, the optimized surface irrigation area, the optimized drip irrigation area, and the unchanged present drip irrigation area.
In the zation scheme, xijk is taken as the decision variable, which is the area of crops k of irrigation technologies j (j = 1 for surface irrigation, j = 2 for drip irrigation) in sub i. Thus, in the optimization scheme, the sub-district comprises three parts, the op surface irrigation area, the optimized drip irrigation area, and the unchanged pres irrigation area.

The Water-Energy-Crops-Irrigation Nexus in the CPSID
As shown in Figure 2, the first pump station pumps water from the Yellow and the design discharge of each pump station decreases successively (Table 2) du reasons. One reason is the loss of water that occurs in each WCP&C segment con the two pump stations. Another reason is that the WCP&C supplies water to vario districts for crop irrigation. The water use efficiency of the WCP&C segment rep the water loss degree of the segment between pump stations (Table 3), which can resented as the ratio of its outflow and inflow. From the fourth pump station, wa the WCP&C is fed into the WDP&CS of the sub-districts to irrigate the crops in th land. The water supply of each sub-district should cover two parts. The first pa irrigation water for crop growth in the sub-district, which is related to both the

The Water-Energy-Crops-Irrigation Nexus in the CPSID
As shown in Figure 2, the first pump station pumps water from the Yellow River, and the design discharge of each pump station decreases successively ( Table 2) due to two reasons. One reason is the loss of water that occurs in each WCP&C segment connecting the two pump stations. Another reason is that the WCP&C supplies water to various sub-districts for crop irrigation. The water use efficiency of the WCP&C segment represents the water loss degree of the segment between pump stations (Table 3), which can be represented as the ratio of its outflow and inflow. From the fourth pump station, water from the WCP&C is fed into the WDP&CS of the sub-districts to irrigate the crops in the farmland. The water supply of each sub-district should cover two parts. The first part is the irrigation water for crop growth in the sub-district, which is related to both the type of irrigation technologies and the type of crops. The other part is the water loss through the WDP&CS in the sub-district and the water loss in the farmland, which is related to the irrigation technology.
The annual irrigation water consumption in the CPSID is related to crop, irrigation technology, the irrigated area, and irrigation water use efficiency. The annual irrigation water consumption W 1 of the CPSID is calculated by Equation (1): where −F 1 is the annual irrigation water consumption of the CPSID (m 3 ); i is the sub-district number; I is the total number of sub-districts; j is the number of irrigation technology types, j = 1 is surface irrigation, j = 2 is drip irrigation; J is the total number of irrigation technology types; k is the number of crop types; K is the total number of crop types; q jk is the net irrigation quota (m 3 /ha); n is the number of WCP&C segments; η n is the water use efficiency of WCP&C segments; and γ ij is the water use efficiency of the sub-district (which is related to the irrigation technology). The annual energy consumption of the CPSID is related to crop, irrigation technology, irrigated area, irrigation water efficiency, pump station efficiency, and pump station head. The annual energy consumption W 2 of the CPSID is calculated according to Equation (2): where F 2 is the annual irrigation water consumption of the CPSID (kWh); p is the pump station number; P is the total number of pump stations; Hp is the pump station head (m); α p is the pump station efficiency (%); h d is the pump station head in the drip irrigation system in the sub-district (m; h d = 50 m), which comprises a head drip irrigation emitter of 10 m and a head loss in the drip irrigation system of 40 m; and ε d is the efficiency of the pump station of the drip irrigation system (%; ε d = 82%). The total output value of various crops F 3 in the CPSID is calculated according to Equation (3): where u jk is the output value of the crops per unit area (yuan/ha), which is related to the type of crop and irrigation technology.

Objective Function
The irrigation water amount, energy consumption, and crop output value in an ID depend on the spatial patterns of crops and irrigation technologies; in other words, the area distribution of crops and irrigation technologies in each sub-district. In the CPSID, the irrigation water consumption, pump station energy consumption, and crop output value vary with different crops and irrigation technologies. Because of differences in planting density of different crops, the investment per unit area of drip irrigation also varies. Thus, improving the investment efficiency of the drip irrigation project could be a solution. The research objective of this paper is to seek an optimal scheme of the spatial patterns of crops and irrigation technologies by making full use of existing engineering facilities in the ID. The optimal scheme should improve the efficiency of drip irrigation project investment: minimize the annual irrigation water and energy consumption and maximize the crop output value in IDs. These conditions indicate a multi-objective optimization problem.
The first objective function W 1 (X) is the minimum annual irrigation water consumption in the CPSID. The second term in Equation (1) is a constant. Therefore, W 1 (X) is calculated according to Equation (4): where X = {x ijk }, X is the set of decision variables. The second objective function W 2 (X) is the minimum annual energy consumption in the CPSID, and the third term in Equation (2) is a constant, so W 2 (X) is expressed as Equation (5): The third objective function Z 1 (X) is the maximum output value of crops in the CPSID, since the second term in Equation (3) is a constant, Z 1 (X) is expressed as Equation (6):

Constraints
(1) Area constraint of each sub-district To prevent the ecological environment degradation caused by blindly expanding the farmland area, the area of each sub-district in the CPSID remains unchanged, as demonstrated in Equation (7): (2) Constraints of the annual water amount for a sub-district The annual water amount for the sub-district is limited by the amount of water from the source, as well as by the capacity of the CPS, the WCP&C, and the WDP&CS of the sub-district. The constraint on the upper limit of annual water amount for the sub-district is shown in Equation (8).
where Q i is the upper limit of annual water amount for the sub-district (m 3 ), as shown in Table 1. (

3) Area constraints of various crops in the CPSID
To both ensure food security and consider the supply and demand relationship of the agricultural products market, the planting area of each crop in the ID should be controlled between its upper and lower limits, as shown in Equations (9) and (10).
where b max k and b min k are the upper limit and lower limit of crop k planting area, respectively ( Table 6). (4) Constraint conditions of drip irrigation project investment The investment of new drip irrigation projects in the ID is constrained by the total investment.
where c k is the investment per unit area of the new drip irrigation project (yuan/ha), which depends on the crop type, because different crops are planted in different densities. T is the total investment quota of the drip irrigation projects in the ID (yuan).

The Algorithm of the Model Solution
This is a multi-objective linear programming model. The objective functions are Equations (4)- (6). The equality constraint is Equation (7), and the inequality constraints are Equations (8)- (12). The decision variable is x ijk . The multi-objective linear programming model is as follows: subject to AX ≤ B, where AX ≤ B represents the inequality constraints (7)- (11), and X ≥ 0 represents the nonnegative constraints (12). Zimmerman solved the above multiple objective programming problem by using the concept of fuzzy set [32]. The membership functions for the objective are defined as where Z * 1 , W * 1 , W * 2 are the ideal solutions of problems (13)- (15) and (16), (17); are the anti-ideal solutions of problems (13)- (15) and (16), (17). The two-stage algorithm of multi-objective linear programming is used to solve Equations (13) and (14). The experimental environment is Windows 10 and MATLAB R2014a. The algorithm steps are as follows: Step 1: Obtain the ideal and the anti-ideal solutions.
Step 2: First stage of the algorithm.
Step 3: Second stage of the algorithm.
where the problems (26)-(30) are a linear programming problem. The objective function is problem (26). The constraints are problems (27)- (30). The data in Tables 1-6 were substituted into MATLAB linear programming toolbox, and obtained the values of λ 1 , λ 2, λ 3, and xijk. According to the above algorithm steps, the flowchart of the algorithm is shown in Figure 4.

Results and Discussion
The drip irrigation investment in the Jingdian Phase Irrigation District was 40,000 × 10 4 yuan, and the value and data presented in Tables 1-6 were input to the objective functions (4)-(6) and the constraint conditions (7)- (12). The ideal solutions (Z1 * , W1 * , W2 * ) and the anti-ideal solutions ( In the first step, the problems (18) The calculation result of the second step is better than that of the first step. The calculation results of the example in literature [34,35] are similar to those in our paper, so we choose the calculation results of the second step.
According to the xijk obtained in the second step and the present area of various crops in drip irrigation ( Table 5). The optimal scheme of the spatial patterns of crops and irrigation technologies in the Jingdian Phase I Irrigation District were obtained (Table 7).

Results and Discussion
The drip irrigation investment in the Jingdian Phase Irrigation District was 40,000 × 10 4 yuan, and the value and data presented in Tables 1-6 were input to the objective functions (4)- (6) and the constraint conditions (7)- (12 In the first step, the problems (18) The calculation result of the second step is better than that of the first step. The calculation results of the example in literature [34,35] are similar to those in our paper, so we choose the calculation results of the second step.
According to the x ijk obtained in the second step and the present area of various crops in drip irrigation ( Table 5). The optimal scheme of the spatial patterns of crops and irrigation technologies in the Jingdian Phase I Irrigation District were obtained ( Table 7).
The data for the optimal scheme (Table 7) can be substituted into Equations (1)-(3) for calculation of the annual irrigation water consumption, annual energy consumption, and crop output value of the ID (Table 8). Compared with the 2020 conditions (Table 5), the optimal scheme could reduce the annual irrigation water and energy consumption by 3722 × 10 4 m 3 and 5515 × 10 4 kWh, respectively. Furthermore, it could increase the crop output value by 19,277 × 10 4 yuan; that is, the scheme could reduce water and energy consumption by 26.18% and 29.38%, and increase income by 29.55%. Table 7. The optimal scheme of spatial patterns of crops and irrigation technologies in the Jingdian Phase Irrigation District (ha).

Sub-District
Surface Irrigation Drip Irrigation Crop abbreviation as shown in Table 4. A comparison of the spatial patterns in 2020 ( Table 5) and that of the optimal scheme (Table 7) reveals that: (1) The spatial distribution of crops in each sub-district of the present situation is less concentrated than that of the optimal scheme. For example, there are currently eight crops in surface irrigation in sub-district 7, but only one crop in surface irrigation in the optimal sub-district 7 ( Table 7). The centralized distribution of crops in the optimal scheme is beneficial to modern agriculture's large-scale and intensive management.
(2) In the optimal scheme, the area of drip irrigation increased from 1660 ha to 9810 ha (Table 8), which is 44.7% of the ID's total area. After the optimization, the areas of drip irrigation in the low elevation sub-districts 1-4 remained almost unchanged, while those in the high elevation sub-districts 5-10 increased obviously ( Table 8). The increase of the drip irrigation area in the high elevation sub-districts was much bigger than in the low elevation sub-districts. The reason is that the higher the elevation of the sub-district, the longer the distance of water transportation, which leads to a greater water loss. For the same crop, the net irrigation quota of drip irrigation is lower than that of surface irrigation (Table 4); therefore, the area of drip irrigation is greatly increased in the high elevation sub-districts, resulting in significant water saving and energy saving effects. (3) In the optimal scheme, the planting structure changed considerably ( Table 7). The areas of grain crops with high water consumption and low output values decreased significantly. For example, the areas of wheat and corn reached the lower limits of the constraint conditions (Tables 6 and 7). By contrast, the areas of crops with low water consumption and high output values increased significantly. For example, the areas of fruit trees, vegetables, and potatoes all reached their upper limits (Tables 6 and 7). In the optimal scheme, the spatial patterns of crops and irrigation technologies were adjusted synchronically. For instance, fruit trees provide a high output value with a low net irrigation quota per unit area and low investment in drip irrigation per unit area (Table 3). Therefore, drip irrigation was applied to all fruit trees while the area reached the upper limits of constraint conditions ( Table 6). The above adjustment is beneficial to water saving, energy saving, and output increasing.
A comparison of the indicators (Table 8) from the present situation ( Table 5) and those of the optimal scheme (Table 7) reveals that: (1) For sub-districts 1 and 2, the two indexes of each sub-district (irrigation water consumption and annual energy consumption) increased after the optimization. However, for sub-districts 3-10, these two indexes were smaller than the present situation (Table 8). If the whole ID is taken into consideration, the two indexes were lower than the present situation. The optimal scheme can save 26.18% water and 29.38% energy compared with the present situation. (2) For sub-districts 2-4, the output value of crops in each sub-district for the optimal scheme was lower than the present situation, while for sub-districts 5-10, the output value of crops increased after the optimization (Table 8). From the perspective of the whole ID, the output value of crops in the optimal scheme was higher than the present situation ( Table 8), indicating that the optimal scheme adjusts crops with a high output value to the high elevation sub-districts. (3) Although the three indexes of the optimal scheme (annual irrigation water consumption, annual energy consumption, and crop output value) were not as good as the present situation in some sub-districts, they were better than the present situation in the rest of the sub-districts (Table 8). From the perspective of the whole ID, the three indexes of the optimal scheme were better than the present situation ( Table 8), indicating that the optimization reached a global solution instead of a local one, which also shows the rationality and correctness of the optimization model and solution method. To increase the drip irrigation area, it is necessary to increase drip irrigation project investment. Different optimal schemes can be obtained when the drip irrigation project investment is 10,000 × 10 4 yuan, 20,000 × 10 4 yuan, 30,000 × 10 4 yuan, 40,000 × 10 4 yuan, and 50,000 × 10 4 yuan, respectively. The annual irrigation water consumption, annual energy consumption, and crop output value of different investment optimization schemes can be obtained by using Equations (1)-(3). The annual irrigation water and energy consumption in the ID demonstrated a decreasing trend with an increase in drip irrigation investment (Figure 5a,b). Energy intensity is an important indicator to measure energy conservation in the CPSID [7]. The energy intensity in the ID demonstrated a downward trend with an increase in drip irrigation investment (Figure 5c), indicating that the CPSID can achieve significant energy conservation by saving water. The crop output value demonstrated a rising trend (Figure 5d).
The whole Jingdian Irrigation District consists of two-phase projects: Jingdian Phase I Irrigation District and Jingdian Phase II Irrigation District. The two projects are independent irrigation systems adjacent to each other. Non-point source pollution and soil salinization are environmental problems closely related to irrigation [50], especially in semi-arid areas where water resources are scarce due to excessive irrigation. The Jingdian Phase I Irrigation District drainage is into the Yellow River. The Jingdian Phase II Irrigation District drainage is into the desert. The irrigation water consumption of the Jingdian Phase I Irrigation District is relatively large, mainly due to the large proportion of high water-consuming crops planted by surface irrigation. In addition, the ID uses many chemical fertilizers; part of the chemical fertilizers enters the Yellow River with the drainage of the ID and pollutes the water source. By optimizing the crop planting structure in Jingdian Irrigation District [27,28], irrigation water consumption can be reduced, and farmland drainage can be reduced by 5.21~8.93%, effectively reducing agricultural non-point source pollution and soil salinization. In the Violada Irrigation District, the nitrate load decreased by 70.4% after replacing the surface irrigation with sprinkler irrigation [51]. If the model proposed here is adopted to optimize crop planting structures and the spatial patterns of irrigation technologies are considered simultaneously, the agricultural irrigation water consumption and farmland drainage in the Jingdian Irrigation District will be reduced, and the pollution degree of farmland drainage to the Yellow River will also be reduced. The whole Jingdian Irrigation District consists of two-phase projects: Jingdian Phase I Irrigation District and Jingdian Phase II Irrigation District. The two projects are independent irrigation systems adjacent to each other. Non-point source pollution and soil salinization are environmental problems closely related to irrigation [50], especially in semi-arid areas where water resources are scarce due to excessive irrigation. The Jingdian Phase I Irrigation District drainage is into the Yellow River. The Jingdian Phase II Irrigation District drainage is into the desert. The irrigation water consumption of the Jingdian Phase I Irrigation District is relatively large, mainly due to the large proportion of high water-consuming crops planted by surface irrigation. In addition, the ID uses many chemical fertilizers; part of the chemical fertilizers enters the Yellow River with the drainage of the ID and pollutes the water source. By optimizing the crop planting structure in Jingdian Irrigation District [27,28], irrigation water consumption can be reduced, and farmland drainage can be reduced by 5.21~8.93%, effectively reducing agricultural non-point source pollution and soil salinization. In the Violada Irrigation District, the nitrate load decreased by 70.4% after replacing the surface irrigation with sprinkler irrigation [51]. If the model proposed here is adopted to optimize crop planting structures and the spatial patterns of irrigation technologies are considered simultaneously, the agricultural irrigation water consumption and farmland drainage in the Jingdian Irrigation District will be reduced, and the pollution degree of farmland drainage to the Yellow River will also be reduced. Drip irrigation can reduce harmful gas emissions [52]. The experimental results show that drip irrigation reduced nitrous oxide and nitric oxide by 13.0% and 16.7%, respectively, compared with surface irrigation [26,53].
Developing a water-saving irrigation method and controlling water consumption are effective measures to protect the environment in arid regions [8][9][10]. Moreover, the area constraint of each sub-district in the optimization model in this paper avoids the trend of ecological deterioration caused by blindly expanding irrigation areas in arid regions [10,11]. The research shows that under the condition of keeping the irrigation area unchanged, the optimization scheme can save water, save energy, and increase crop income. This is conducive to the simultaneous improvement of economic, social, and ecological benefits and to guarantee the sustainable development of the CPSID.
The quantitative analysis of the ecological and environmental benefits of the optimized spatial patterns of crops and irrigation technologies in the Jingdian Phase I Irrigation District should be further explored. The contents to be further studied include soil salinization, non-point source pollution, and harmful gas emissions reduction by drip irrigation in the CPSIDs. These problems are also the limitations of this study. Further study could supplement the ecological environment benefit maximization of the target function. Drip irrigation has the effect of reducing the salt alkalinity, non-point source pollution, and harmful gas discharge. Based on this effect, the constraint is added in the model. Of course, these further studies require a lot of experiments, using measured data, obtaining the parameters required in the model, and improving the model. We will also carry out further research regarding the uncertainty of the model parameters [35][36][37][44][45][46].

Conclusions
The multiple objective programming problem models is applied to the Jingdian Phase I Irrigation District. Under the condition of 40,000 × 10 4 yuan of drip irrigation project investment, the optimal scheme of the spatial pattern of crops and irrigation technologies in the ID was obtained. Compared with the present situation, the annual irrigation water consumption was reduced by 3722 × 10 4 m 3 , and the energy consumption was reduced by 5515 × 10 4 kWh, whereas the crop output value was increased by 19,277 × 10 4 yuan. This optimization saves water by 26.18%, saves energy by 29.38%, and increases income by 29.55%. The optimal spatial pattern of the ID suggested that crops with low water consumption and high output values should be distributed in the sub-districts with high elevations. Meanwhile, crops with high water consumption and low output values should be distributed in the sub-districts with low elevations. The results showed that the annual irrigation water consumption, energy consumption, and energy intensity decreased with the increasing investment in drip irrigation projects, while the crop output value increased.
Optimizing the spatial patterns of crops and irrigation technologies in the CPSID is vital to promoting the sustainable development of agriculture and ecological environment protection in arid regions.