## Abstract

The demand for water resources in Shaanxi Province increases greatly due to the continuous growth of its population and the rapid development of the social economy. Water demand forecasting is a significant issue in the designing, maintaining and operating of a reliable and economical water supply system. An explicit mathematical method was presented in this study, based on the indicators of industrial output value, such as the gross output value of agriculture, forestry, animal husbandry and fishery. The impact of total retail sales and year trends in the domestic or industrial water demands, can accurately forecast the water demand fluctuations for a municipality. Adopt RAGA-PP optimal selection model through a grouping method of data handling for water demand management to test for the case study Shaanxi, China. Results showed that the prediction effect of multivariate logarithmic model accuracy can reach 99.50%, and it is estimated that the demand for water resources in Shaanxi would exceed 10 billion m^{3} by 2030. The average relative error of the water consumption from 2015 to 2017 is 3.05% for the model of multiple linear and 0.50% for the model of multivariate logarithm model. Our framework can assist in developing sustainable solutions.

## HIGHLIGHTS

Using the RAGA-PP model to realize the optimal selection of water resources demand indicators can improve model accuracy compared with model solving.

Adopt a curve fitting method and regression analysis method to establish the water demand model and solve it.

Due to the huge consumption of agricultural water, this paper puts forward corresponding agricultural water-saving measures.

### Graphical Abstract

## INTRODUCTION

There is a great demand for freshwater because of the rapid development of the social economy and its growing urban regions in Shaanxi Province. According to statistics, the total water reserve of Shaanxi from 2009 to 2018 was lower than that of China, only 40.498 billion m^{3}, and the average water resources per capita was only half the national average level. The potable water supply for Shaanxi depends on its rivers and groundwater sources. In recent years, due to the over-exploitation of groundwater in Shaanxi, its level in some areas has dropped sharply, and some areas are even threatened with groundwater depletion. The water resource shortage in Shaanxi is becoming serious. Reasonable control of the water demand among the residents, businesses and other customers in the community through forecasting is a key measurement in creating a reliable and economic system.

Researchers have successfully applied time series modeling based on stochastic modeling (linear) and/or artificial intelligence methods (non-linear) for water demand forecasting. Zubaidi *et al.* (2020) adopted a new methodology combining discrete wavelet transform (DWT) with an adaptive neuro-fuzzy inference system (ANFIS) which was proposed to predict monthly urban water demand based on several intervals of historical water consumption. Lu *et al.* (2020) adopted a hybrid model to improve the accuracy of predicting water demand. The model combined a simulated annealing (SA) algorithm, genetic cross factor (GCF), fruit fly optimization algorithm (FFOA), and support vector machine (SVM) to predict monthly water demand. Yuan *et al.* (2019) employed a grey power model [GPM (1, 1) model] which was a good prediction method for predicting urban water consumption. The traditional GPM (1, 1) model generated its grey sequence by a first-order accumulating generation operator (1-AGO) and produced the predicted results by a first-order inverse accumulating generation operator (1-IAGO). Smolak *et al.* (2020) compared classical and adapted machine learning algorithms used for water usage predictions including autoregressive integrated moving average (ARIMA), support vector regression, random forests and extremely randomized trees. Chen *et al.* (2017) adopted a multiple random forests model, integrated wavelet transform and random forests regression (W-RFR) which were proposed for the prediction of daily urban water consumption in southwest of China. Jain & Ormsbee (2002) adopted conventional regression and time series analysis methods which all were artificial neural network models to model and forecasted short-term water demand. They considered the rainfall, the maximum air temperature and the past water demands as effective parameters for water demand forecast model. Herrera *et al.* (2010) employed different learning machine methods to forecast hourly urban water demand. Wu & Yan (2010) utilized two genetic programming (GP) approaches to find a model for water demand forecast for a real water system, which showed that the average temperature had a significant impact on water demand. Gharabaghi *et al.* (2019) introduced a linear transfer function forecasting model based on the ARIMA to forecast monthly water consumption and considered several climatic, economic and flow rates as inputs for the model development in El Paso city in Texas, USA. As a result, changes in climatic and economic parameters were compared to changes in water rates parameters, which had a greater impact on consumption per customer**.**Lopez *et al.* (2018) applied a multi-model predictor architecture for water demand forecasting in water distribution networks. Yin *et al.* (2018) developed a feed-forward artificial neural network model to forecast urban water resources and energy demand. Ghalehkhondabi *et al.* (2017) reviewed soft computing techniques for water demand forecast. More recently, Bonakdari *et al.* (2018) introduced a new insight in time series modeling by detailing three linear/nonlinear/hybrid methods. Zeynoddin *et al.* (2018) proposed a new methodology for modeling of a time series. This methodology was based on a combination of stationarization techniques, and normality transforms into three different scenarios. Zeynoddin *et al.* (2018) and Ebtehaj *et al.* (2019) highlighted the impact of appropriate pre-processing techniques for elimination of deterministic term in time series, leading to a reliable linear approach with higher accuracy compared to non-linear approaches. Sanchez *et al.* (2020) used a Geographically Weighted Regression model informed by socio-economic, environmental and landscape pattern metrics to evaluate how future scenarios of population densities and climate warming would jointly affect water demand across two rapidly growing US States (North Carolina and South Carolina). Results showed that future water demand would increase under rising temperatures, but could be ameliorated by policies that promoted higher density development and urban infill. These water-efficient land-adopted policies led to a 5% regional reduction in water demand and up to 25% reduction locally for counties with the highest expected population growth by 2065. For rural counties experiencing depopulation, we considered that the land policies carried out were insufficient to significantly reduce water demand.

Balha *et al.* (2020) analyzed the impact of urban land-use change on Delhi's water resources during time period 2005–2016 and for 2031. Using Monte-Carlo simulations in Land Change Modeler (LCM), they predicted the future LULC for year 2031 which exhibited a similar pattern of LULC change as observed during 2005–2016. An increase of 36.49% in urban area was observed during time period 2005–2016 and an increase of 14.05% was predicted for 2016–2031. Adopting in-situ groundwater measurements, the groundwater abstract (−18.75%) was found to be greater than the amount of total groundwater recharge (14.67%).

In this paper, we aimed to propose a model with an explicit equation with high accuracy for forecasting the water demand. The projection pursuit (PP) model optimized by the accelerated genetic algorithm (RAGA) was used to study the contribution degree of each evaluation index of water consumption.

## DATA DESCRIPTION

Shaanxi Province is located in the northwest of China. It spans the middle of the Yellow River and the Yangtze River, connecting the eastern and central regions of China with important hubs in the northwest and southwest. The values of total water resources in Shaanxi from 2009 to 2018 are shown on Figure 1. The total amount of water resources in Shaanxi increased from 2009 to 2011 and has increased since 2017, but it is slightly lower than that of a decade ago. According to the *Statistical Yearbook of Shaanxi Province*, the total amount of water resources in Shaanxi Province was closely related to the amount of precipitation. In 2011, the amount of precipitation was the greatest, which made the total amount of water resources in that year the largest.

The water demand of Shaanxi mainly includes four parts: industrial water demand, agricultural water demand, domestic water demand, and ecological water demand. According to the average water consumption percentage of Shaanxi Province from 2004 to 2014, the highest agricultural water consumption was 67.95%, while the proportion of industrial water consumption, domestic water consumption and the ecological water consumption were 15.39%, 14.57%, and 2.08% respectively. The results of annual average water consumption ratio of four parts in Shaanxi during 2004–2014 are shown on Figure 2. On the premise of the model accuracy less affected, the ecological water consumption with the least share could be ignored in order to simplify the influencing factors of water resources demand in Shaanxi Province. It was reasonable to ignore ecological water consumption in terms of the accuracy of prediction results. The following analysis mainly considered industrial water demand, agricultural water demand and domestic water demand.

## METHODS

### Projection pursuit model based on accelerated genetic algorithm

Projection pursuit was proposed by Friedman & Tukey (Liu *et al.* 2019). It is mainly used to project high-dimensional data into low-dimensional subspace. The best projection direction and projection value can be obtained after conversion.

In this paper, the contribution rate of each evaluation index and the change trend of water resources demand in Shaanxi Province can be judged respectively.

A PP model optimized by accelerated genetic algorithm was used to study the contribution of water consumption evaluation indexes in Shaanxi Province, and to determine the evaluation indexes of water resource demand model.

### Normalization of evaluation indicators

Let the sample set to be evaluated as { | *i* = 1,2, …,n; *j* = 1,2, …,m}, where is the *j*th index of the *i*th sample. In order to unify the variation range of each index value and eliminate the influence of each index dimension, it is necessary to normalize each index.

### Constructive projection function

### Based on RAGA optimized projection function

The following steps are taken to optimize the projection function using accelerated genetic algorithm:

*m*of the

*n*group {}, which is converted from equation (8) to decimal, and then from equation (7) to binary. Equation (8) is shown as follows:

Step 3: Evaluation of individual adaptability of parent generation: substitute *i* group parameters into equation (5) and equation (6) to get f_{i} projection function. The adaptability of the th group of parameters is inversely proportional to f_{i}.

Step 4: The probability selection of the parent generation: the generated *n* set of random numbers {}, if, select the th individual of the parent generation.

Step 5: Crossing the parent generation: the th pair of parents is transformed into binary numbers to produce a new th pair of offspring.

Step 6: Subgenerational individual variation: when it is less than the probability of offspring variation, the offspring individual will perform a flip operation.

Step 7: Iterative evolution: take th offspring as the parent and substitute Step 3 into the next round until the optimization criterion of the optimal individual is less than the specific threshold.

Step 8: Accelerating the cycle: use the excellent individuals generated by each iteration to adjust the variable interval to make it closer to the optimal point.

### Selection of assessment indicators

Many indicators can reflect the demand for water resources in Shaanxi Province, such as the quota of 10,000 Yuan (100 Yuan = 15 US dollars) output value, industrial added value, water consumption for industrial increase, industrial output value, irrigation area, total output value of agriculture, forestry, animal husbandry and fishery, total sowing area, population quantity, per capita disposable income of urban residents, total retail sales of consumer goods, water consumption quota, etc. Regional water demand assessment indicators are shown in Table 1.

The PP model parameters of accelerated genetic algorithm were set as follows: population size N was 400, crossover probability Pc was 0.8, mutation probability Pm was 0.2, random number required for variation direction M was 10, acceleration times Ci was 7, number of variables n was 9, and limit number set for acceleration once after two evolution was D2. Optimal projection direction and projection values of the nine indexes C1–C9 the obtained index layer are shown in Tables 2 and 3. The evaluation index of water resource demand selects the maximum of the best projection direction from each of the three system layers. That was to choose C3, C5 and C9 as the evaluation index of water resources demand.

Target layer . | System layer . | Indicator layer . |
---|---|---|

Water resources demand in Shaanxi A | Industrial water demand B1 | 10,000 Yuan output value quota C1 (×10^{8} Yuan) |

Industrial value added C2 (×10^{8} Yuan) | ||

Industrial output value C3 (×10^{8} Yuan) | ||

Agricultural water demand B2 | Irrigation area C4 (1,000 Ha) | |

Gross output value of agriculture, forestry, animal husbandry and fishery C5 (×10^{8} Yuan) | ||

Gross sowing area C6 (1,000 Ha) | ||

Living water demand B3 | Population C7 (10,000) | |

C8 (Yuan per capita disposable income of urban residents) | ||

Gross retail sales of consumer goods C9 (×10^{8} Yuan) |

Target layer . | System layer . | Indicator layer . |
---|---|---|

Water resources demand in Shaanxi A | Industrial water demand B1 | 10,000 Yuan output value quota C1 (×10^{8} Yuan) |

Industrial value added C2 (×10^{8} Yuan) | ||

Industrial output value C3 (×10^{8} Yuan) | ||

Agricultural water demand B2 | Irrigation area C4 (1,000 Ha) | |

Gross output value of agriculture, forestry, animal husbandry and fishery C5 (×10^{8} Yuan) | ||

Gross sowing area C6 (1,000 Ha) | ||

Living water demand B3 | Population C7 (10,000) | |

C8 (Yuan per capita disposable income of urban residents) | ||

Gross retail sales of consumer goods C9 (×10^{8} Yuan) |

Index . | B1 . | B2 . | B3 . | ||||||
---|---|---|---|---|---|---|---|---|---|

C1 . | C2 . | C3 . | C4 . | C5 . | C6 . | C7 . | C8 . | C9 . | |

Optimal projection direction | 0.33 | 0.36 | 0.37 | 0.34 | 0.38 | 0.13 | 0.31 | 0.35 | 0.36 |

Index . | B1 . | B2 . | B3 . | ||||||
---|---|---|---|---|---|---|---|---|---|

C1 . | C2 . | C3 . | C4 . | C5 . | C6 . | C7 . | C8 . | C9 . | |

Optimal projection direction | 0.33 | 0.36 | 0.37 | 0.34 | 0.38 | 0.13 | 0.31 | 0.35 | 0.36 |

Year . | 2004 . | 2005 . | 2006 . | 2007 . | 2008 . | 2009 . | 2010 . | 2011 . | 2012 . | 2013 . | 2014 . |
---|---|---|---|---|---|---|---|---|---|---|---|

Projection value | 0.36 | 0.49 | 0.53 | 0.73 | 1.04 | 1.15 | 1.46 | 1.84 | 2.15 | 2.43 | 2.67 |

Water consumption (×10^{8}m^{3}) | 75.53 | 78.76 | 84.07 | 81.54 | 85.46 | 84.34 | 83.40 | 87.76 | 88.04 | 89.21 | 89.81 |

Year . | 2004 . | 2005 . | 2006 . | 2007 . | 2008 . | 2009 . | 2010 . | 2011 . | 2012 . | 2013 . | 2014 . |
---|---|---|---|---|---|---|---|---|---|---|---|

Projection value | 0.36 | 0.49 | 0.53 | 0.73 | 1.04 | 1.15 | 1.46 | 1.84 | 2.15 | 2.43 | 2.67 |

Water consumption (×10^{8}m^{3}) | 75.53 | 78.76 | 84.07 | 81.54 | 85.46 | 84.34 | 83.40 | 87.76 | 88.04 | 89.21 | 89.81 |

Figure 3 shows the demand for water resources in Shaanxi. From the overall trend of change, the trend of projection value was basically in line with that of water resources demand. This can indicate that the selection of C1-C9 index was reasonable, and was appropriate to use the RAGA-PP model to select the optimal index of water resources demand prediction in Shaanxi.

### Establishment and solution of water resource demand forecasting model

#### Model establishment

The total output value of agriculture, forestry, animal husbandry and fishery (), industrial output value (), total retail sales of consumer goods () and water consumption (*Y*) of the year were fitted to the curves in Table 4, showing that the best fit between the three indexes and water consumption was logarithmic function. Many fit results showed that the Log3P1 and Bradley functions fitted best. All correlation coefficient Rs were above 0.9, and the fit results are shown in Figure 4.

Year . | Water consumption . | Industrial output value . | Gross output value of agriculture, forestry, animal husbandry and fishery . | Total retail sales of consumer goods . |
---|---|---|---|---|

2004 | 75.53 | 3,389.88 | 651.21 | 1,162.80 |

2005 | 78.76 | 4,109.32 | 730.72 | 1,331.35 |

2006 | 84.07 | 5,248.39 | 821.54 | 1,542.37 |

2007 | 81.54 | 6,587.41 | 1,002.85 | 1,837.25 |

2008 | 85.46 | 8,358.86 | 1,277.86 | 2,317.11 |

2009 | 84.34 | 9,553.70 | 1,337.22 | 2,725.67 |

2010 | 83.40 | 12,421.80 | 1,666.06 | 3,257.54 |

2011 | 87.76 | 15,811.48 | 2,058.60 | 3,900.58 |

2012 | 88.04 | 18,591.89 | 2,303.20 | 4,581.62 |

2013 | 89.21 | 20,820.07 | 2,562.51 | 5,245.04 |

2014 | 89.81 | 21,944.58 | 2,741.82 | 5,918.71 |

Year . | Water consumption . | Industrial output value . | Gross output value of agriculture, forestry, animal husbandry and fishery . | Total retail sales of consumer goods . |
---|---|---|---|---|

2004 | 75.53 | 3,389.88 | 651.21 | 1,162.80 |

2005 | 78.76 | 4,109.32 | 730.72 | 1,331.35 |

2006 | 84.07 | 5,248.39 | 821.54 | 1,542.37 |

2007 | 81.54 | 6,587.41 | 1,002.85 | 1,837.25 |

2008 | 85.46 | 8,358.86 | 1,277.86 | 2,317.11 |

2009 | 84.34 | 9,553.70 | 1,337.22 | 2,725.67 |

2010 | 83.40 | 12,421.80 | 1,666.06 | 3,257.54 |

2011 | 87.76 | 15,811.48 | 2,058.60 | 3,900.58 |

2012 | 88.04 | 18,591.89 | 2,303.20 | 4,581.62 |

2013 | 89.21 | 20,820.07 | 2,562.51 | 5,245.04 |

2014 | 89.81 | 21,944.58 | 2,741.82 | 5,918.71 |

Annual water consumption unit: ×10^{8} m^{3}; Indicator unit: ×10^{8} Yuan.

According to the fitting of three indexes and water consumption, this thesis builds a multivariate logarithmic model to predict water resources demand in Shaanxi, and establishes a multivariate linear model to make a comparison.

#### Model solution

Regression analysis of multivariate linear and multivariate logarithmic models using R software, and relevant output results were obtained by stepwise analysis, as shown in Tables 5 and 6.

Models . | R^{2}
. | Errors in standard estimates . | Errors in Durbin–Watson statistic . |
---|---|---|---|

Multivariate linear models | 0.891 | 2.411 | 1.593 |

Multivariate logarithmic models | 0.951 | 1.650 | 2.559 |

Models . | R^{2}
. | Errors in standard estimates . | Errors in Durbin–Watson statistic . |
---|---|---|---|

Multivariate linear models | 0.891 | 2.411 | 1.593 |

Multivariate logarithmic models | 0.951 | 1.650 | 2.559 |

Models . | Unstandardized coefficient . | Standardized coefficient . | t . | Significance . | |
---|---|---|---|---|---|

B . | Beta . | ||||

Multivariate linear models | Constant total retail sales of consumer goods | 74.279 | – | 10.965 | 0.000 |

Industrial output value | − 0.000448 | − 0.685 | − 0.168 | 0.871 | |

Gross output value of agriculture, forestry, animal husbandry and fishery | 0.015123 | 2.570 | 0.613 | 0.559 | |

Total retail sales of consumer goods | − 0.002712 | − 1.003 | − 0.505 | 0.629 | |

Multivariate logarithmic models | Constant total retail sales of consumer goods | 1.519 | – | 0.143 | 0.890 |

Industrial output value | 3.332908 | 3.149 | 1.020 | 0.342 | |

Gross output value of agriculture, forestry, animal husbandry and fishery | − 2.354693 | − 2.210 | − 0.821 | 0.439 | |

Total retail sales of consumer goods | 0.003782 | 0.004 | 0.001 | 0.999 |

Models . | Unstandardized coefficient . | Standardized coefficient . | t . | Significance . | |
---|---|---|---|---|---|

B . | Beta . | ||||

Multivariate linear models | Constant total retail sales of consumer goods | 74.279 | – | 10.965 | 0.000 |

Industrial output value | − 0.000448 | − 0.685 | − 0.168 | 0.871 | |

Gross output value of agriculture, forestry, animal husbandry and fishery | 0.015123 | 2.570 | 0.613 | 0.559 | |

Total retail sales of consumer goods | − 0.002712 | − 1.003 | − 0.505 | 0.629 | |

Multivariate logarithmic models | Constant total retail sales of consumer goods | 1.519 | – | 0.143 | 0.890 |

Industrial output value | 3.332908 | 3.149 | 1.020 | 0.342 | |

Gross output value of agriculture, forestry, animal husbandry and fishery | − 2.354693 | − 2.210 | − 0.821 | 0.439 | |

Total retail sales of consumer goods | 0.003782 | 0.004 | 0.001 | 0.999 |

The results of solving the multivariate linear and multivariate logarithmic models are shown in the equation (13) and equation(14), respectively.

## VERFICATION INDICES

The total output value of agriculture, forestry, animal husbandry and fishery (), industrial output value () and total retail sales of consumer goods () in Table 7 were substituted into multivariate linear and multivariate logarithmic models to obtain the and of water consumption in 2015 and 2017, respectively. Table 7 shows the actual water consumption and three index values during 2015–2017, while Table 8 shows the prediction results and relative errors of two models.

Year . | Water consumption . | Industrial output . | Total value of agriculture, forestry, animal husbandry and fishery . | Total retail sales of consumer goods . |
---|---|---|---|---|

2015 | 91.16 | 21,060.03 | 2,813.50 | 6,578.14 |

2016 | 90.83 | 22,549.34 | 2,985.76 | 7,367.57 |

2017 | 92.99 | 24,726.17 | 3,070.45 | 8,113.84 |

Year . | Water consumption . | Industrial output . | Total value of agriculture, forestry, animal husbandry and fishery . | Total retail sales of consumer goods . |
---|---|---|---|---|

2015 | 91.16 | 21,060.03 | 2,813.50 | 6,578.14 |

2016 | 90.83 | 22,549.34 | 2,985.76 | 7,367.57 |

2017 | 92.99 | 24,726.17 | 3,070.45 | 8,113.84 |

Annual water consumption unit: ×10^{8} m^{3}; Indicator unit: ×10^{8} Yuan.

Year . | Actual value . | Multiple linear . | Multiple logarithm . | ||||
---|---|---|---|---|---|---|---|

Forecast . | Error value . | Relative error (%) . | Multiple logarithm . | Error value . | Relative error (%) . | ||

2015 | 91.16 | 89.55 | 1.61 | 1.76 | 91.00 | 0.16 | 0.17 |

2016 | 90.83 | 89.35 | 1.48 | 1.63 | 91.34 | − 0.51 | 0.56 |

2017 | 92.99 | 87.63 | 5.36 | 5.76 | 92.27 | 0.72 | 0.78 |

Average error | 3.05 | 0.50 |

Year . | Actual value . | Multiple linear . | Multiple logarithm . | ||||
---|---|---|---|---|---|---|---|

Forecast . | Error value . | Relative error (%) . | Multiple logarithm . | Error value . | Relative error (%) . | ||

2015 | 91.16 | 89.55 | 1.61 | 1.76 | 91.00 | 0.16 | 0.17 |

2016 | 90.83 | 89.35 | 1.48 | 1.63 | 91.34 | − 0.51 | 0.56 |

2017 | 92.99 | 87.63 | 5.36 | 5.76 | 92.27 | 0.72 | 0.78 |

Average error | 3.05 | 0.50 |

In Table 8, the average relative errors of water consumption during 2015–2017 were 3.05% and 0.50% respectively. The multivariate logarithmic model was used to predict the water resources demand in Shaanxi Province with high accuracy and better prediction results.

## RESULTS AND DISCUSSION

### The results of water demand forecast in Shaanxi

In order to predict the water resources demand in Shaanxi Province in the future, it was necessary to establish a functional relationship between the year and each index, and the index relation with the year obtained on Figure 5. With the exponential function increasing rapidly, the precision of the coefficient goes higher. So it was necessary to logarithmic the index.

^{2}were close to 1. The functions are listed as follows:where, the , , stands the total output value of agriculture, forestry, animal husbandry and fishery, the industrial output value, and the total retail sales of consumer goods, the

*X*indicates the year, respectively

By substituting equations (15–17) into equation (14), the functional relationship between year and water resources demand in Shaanxi can be obtained. The water resources demand of Shaanxi in 2020, 2025 and 2030 was shown in Table 9.

Year . | Water Resource Demand (×10^{8}m^{3})
. |
---|---|

2020 | 95.78 |

2025 | 99.79 |

2030 | 103.81 |

Year . | Water Resource Demand (×10^{8}m^{3})
. |
---|---|

2020 | 95.78 |

2025 | 99.79 |

2030 | 103.81 |

According to the data above, Shaanxi Province will need 9.578 billion m^{3} of water in 2020 and 9.979 billion m^{3} in 2025 which is 0.401 billion m^{3} more than that of 2020. Compared with 2025, the demand for water resources will be 10.381 billion m^{3} in 2030, which is 0.402 billion m^{3} higher.

Thereafter, this paper concludes that the demand for water resources in Shaanxi is increasing year by year. However, on Figure 1, the total amount of water resources has decreased in the past ten years which would cause many water resources shortage problems.

### Various water-saving measures

Figure 2 shows that the agricultural water demand accounts for the largest proportion, followed by industrial water demand, domestic water demand, and ecological water demand. So, this paper proposes the following suggestions:

- a.
Adopt water-saving irrigation technology, such as sprinkler irrigation, drip irrigation, seepage irrigation. As one of the most advanced water utilization technologies, this kind of irrigation technology has no surface evaporation and deep leakage loss.

- b.
Use soil water storage and moisture conservation technology. Increase soil water storage capacity by fertilizer.

- c.
Utilize water surface seepage prevention treatment. In order to, increase runoff rate and reduce surface storage and infiltration capacity, water surface seepage prevention treatment is an effective technical means.

- d.
Adopt solid water planting technology.

## CONCLUSIONS

- (1)
The RAGA-PP model to optimize the selection of water resources demand indicators is reasonable.

- (2)
The average relative error of water resources demand in Shaanxi during 2015–2017 predicted by multivariate logarithmic model was only 0.50%, and the accuracy of prediction results was 6 times than that of the multivariate linear model.

- (3)
Demand for water resources in Shaanxi will continue to increase in the future, and it is expected to reach close to 10 billion m

^{3}in 2025 and 12 billion m^{3}in 2030. With the decrease of total amount of water resources year by year, it would face a series of problems of water resources shortage, which means that various water-saving measures should be considered.

## FUNDING

This research is supported by the fund of Industrial Support Plan Project of Gansu Province colleges and universities, No. 2020C-38.

## CONFLICTS OF INTEREST

The authors declare no conflict of interest.

## DATA AVAILABILITY STATEMENT

All relevant data are included in the paper or its Supplementary Information.