Next Article in Journal
On the Optimal Input Rate in Queues with Batch Service
Next Article in Special Issue
A New Reverse Extended Hardy–Hilbert’s Inequality with Two Partial Sums and Parameters
Previous Article in Journal
Analysis of Water Infiltration under Impermeable Dams by Analytical and Boundary Element Methods in Complex
Previous Article in Special Issue
Robust Consensus in a Class of Fractional-Order Multi-Agent Systems with Interval Uncertainties Using the Existence Condition of Hermitian Matrices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Forecasting Water Consumption in the Yangtze River Delta Based on Deformable Cumulative Multivariable Grey Model

1
School of Economics and Management, Shanxi University, Taiyuan 030006, China
2
Institute of Management and Decision, Shanxi University, Taiyuan 030006, China
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(7), 655; https://doi.org/10.3390/axioms12070655
Submission received: 22 May 2023 / Revised: 27 June 2023 / Accepted: 28 June 2023 / Published: 1 July 2023
(This article belongs to the Special Issue Advances in Analysis and Control of Systems with Uncertainties II)

Abstract

:
The intensified contradiction between water resources and social development has restricted the development of the Yangtze River Delta. Due to the importance of water consumption in relieving this contradiction, this paper proposes a novel cumulative multivariable grey model with a high performance to predict the water consumption. Firstly, the grey correlation analysis is applied to study the influencing factors, and then the DGM(1,N) with deformable accumulation (DDGM(1,N) model) is constructed and used to predict the water consumption. The results show that the resident population has a significant impact on the water consumption, and the performance of the DDGM(1,N) model is better than the other two grey models. Secondly, the proposed novel grey model is applied to predict the water consumption in 17 cities in the Yangtze River Delta, and the predicted water consumption in Zhejiang and Shanghai indicates a downward trend, while the predicated water consumption in some cities of the Anhui Province presents an upward trend, such as Chizhou, Chuzhou, Wuhu and Tongling. Finally, some policy implications are provided that correspond to the population growth and three major industries in different situations. This paper enriches the research method and prediction analysis used for the water consumption, and the findings can provide some decision-making references for water resources management.

1. Introduction

Water occupies an important strategic position in the development of society, so it is undoubtedly an extremely precious resource [1]. However, rapid economic development leads to an increased demand for resources, and a new problem arises regarding how to seek a balance between resource protection and economic development. Due to the impact of external complex factors, many cities in China have experienced varying degrees of water shortages [2]. For example, the functionality of water resources is reduced in cities such as Hangzhou and Ningbo due to the inadequate management of water resources, which has caused adverse consequences that can lead to the local area bearing heavier pressure in future developments. Therefore, accurately predicting the water consumption can help to mitigate possible future conflicts during sustainable development and environmental protection.
In recent years, many scholars have been devoted to the forecasting research of water resources and have achieved many innovative results, which provided some new ideas for the prediction of water consumption, such as neural networks [3] and support vector machines [4]. Meanwhile, many well-performing models were used to study the issues related to water resources. A multi-objective model was proposed to study the sustainable management of water resources in the Zarrineh Basin in Iran [5]. A support vector regression model was used for short-term predictions for eight water supply areas in the Netherlands and Belgium [6]. Aiming at the prediction problem of urban water demand, the bat algorithm is introduced to improve the prediction accuracy [7]. The interaction between Xinjiang’s population, water resources and economic development was studied based on the system dynamics model [8]. A long- and short-term memory model was used to predict the short-term demand in Hefei, China [9]. A water scarcity pricing model was proposed to promote more effective water resources management in areas with scarce water resources [10].
Among the studies on water consumption requirements, the influential factors and model specification will have important impacts on the analysis results. Considering the climate and groundwater conditions, the Colne watershed in eastern England was analyzed under different drought conditions, and the population factors of the basin were integrated to predict the water demand [11]. By combining the experimental design and Monte Carlo simulation, the prediction uncertainty of the water quality was evaluated [12]. A quantitative method considering water resource distribution was proposed to study the relationship between climate change and water demand [13]. A new demand forecasting method combining machine learning and statistical methods was proposed to give short-term forecasts for the water demand of British households [14]. The evaluation and planning model was used to predict the water demand in the central and western regions of Côte d’Ivoire, which assisted with formulating water management optimization measures [15].
In addition, the problem of predicting water resources using grey models has continued to attract the attention of scholars, which produced a fruitful contribution to this field [16]. Considering the effect of various factors, some policy suggestions for the safety of agricultural water use were put forward based on the predictive analysis of the change in agricultural water use [17]. The grey models, by optimizing objectives and parameters, are used to predict the water level quantity and quality [18,19]. An improved grey model, which utilizes the exponentially weighted moving average and particle swarm algorithm, was established to predict water emissions from faults [20]. By combining multiple factors and grey relational analysis, a safety evaluation model was established to analyze the safety problems on water resources in Guizhou Province [21]. A new grey model was put forward to evaluate the rural water environment quality by comprehensively using the network search information and principal component analysis [22]. The grey water footprint index was used to study the impact of changes in residents’ consumption levels and patterns on water pollution [23]. By integrating the grey decision method and entropy weight method, a multivariable model was established to offer a novel way to study mine water sources [24]. The fractional grey model was used to make water consumption recommendations for 31 provinces [25].
Since Professor Deng founded the grey system theory in 1982, the grey system theory has been widely used in various fields [26]. The grey prediction model is an important part of the grey system theory. It uses the grey generating operator or sequence operator to weaken randomness, excavate the potential law, establish the continuous dynamic differential equation based on the discrete data series, and make the quantitative prediction of the research object according to the GM model [27]. The traditional grey model (GM(1,1)) is represented with one order difference and one variable. In order to improve the prediction performance of the model in forecasting problems, many scholars have improved the grey model in the following aspects: (1) the optimization of the model background value [28]; (2) the extension of the modeling equation [29]; (3) the improvement of the cumulative generation method [30]. In this study, the grey generating operator of the model is optimized by introducing the deformable accumulation generating operator. The detailed modeling process of the prediction model is introduced in Section 3.2.
Because the traditional statistical model has relatively high requirements regarding sample data, it requires a large sample size of the research object and certain sampling probability distribution [31]. When insufficient information is used to meet the conditions, the prediction effect of the statistical model is often unsatisfactory. Additionally, traditional econometrics has some shortcomings in the research of forecasting problems. Although traditional econometrics, which is based on economic theory, has excellent interpretability, its forecasting ability is relatively inferior [32]. The grey prediction model has low requirements for the amount of information, as low as four data points, and has the advantages of easy operation and a high modeling accuracy [33]. With the development of the grey system theory, the grey prediction model has been widely used in many forecasting fields and has become an effective tool to deal with forecasting problems [34]. The grey prediction model provides a new way to solve the prediction problem, and the research on water resources prediction has been a hot topic. Finally, the grey prediction model is used as a research tool in this study. We optimize and improve the traditional grey model and use it to forecast water consumption.
Therefore, the accurate prediction of water consumption is particularly important, as traditional models are complicated and require too much relevant information. The forecasting model utilized in this paper integrates a novel type of accumulation method [35], which has a better forecasting effect through the processing of the original data. According to the principle of using small samples and new information first in the grey system [36], a new cumulative grey prediction model was established for the data and information from recent years. Because the finiteness of the new information and the validity of the initial value are fully considered, the new model has a better prediction accuracy. At the same time, we learn from the research methods of other scholars and use grey correlation analysis in the selection of factors affecting water consumption [37]. Finally, through the combination of the new cumulative multivariable grey model and grey correlation analysis, this paper predicts the water consumption in the Yangtze River Delta region.
The rest of this paper is organized as follows. Section 2 states the research area and data sources. The methods are provided in Section 3. The forecasting results and discussion are presented in Section 4. Section 5 gives the conclusion and policy implications.

2. Research Area and Data Source

The research area covered in this paper is shown in Figure 1, which involves 26 cities in three provinces (Zhejiang, Anhui and Jiangsu) and Shanghai, with a land area of 211.7 thousand square kilometers, accounting for about 2.2% of China’s area. This area has many rivers, seas and coastal ports along the river, and owns the best basis for urbanization in China. However, regional resources and environmental pressures have restricted the economic development. In particular, the contradiction between freshwater resources and urban consumption is gradually aggravating. Therefore, it is imperative to forecast the water consumption in this area to ensure the sustainable development of regional economy.
The data sources of this paper are the Statistical Yearbook of local statistics websites. The data of Shanghai, Zhejiang and Anhui are from the official website of the Municipal Bureau of Statistics (http://tjj.sh.gov.cn/, http://tjj.zj.gov.cn/, http://tjj.ah.gov.cn/) (accessed on 15 July 2022). Since the water consumption in each city in Jiangsu Province cannot be found in the Statistical Yearbook of Jiangsu Provincial Bureau of Statistics, Jiangsu Province is not included in the study area of this paper. Therefore, 17 cities are selected as the research areas of this paper.

3. Methods

The grey correlation in this paper is applied to analyze the relationship between water consumption and influence factors. Then, the DGM(1,N) with deformable accumulation (DDGM(1,N) model) is established to study the change trend of water consumption. Finally, the validity of the model is proven via model comparative analysis.

3.1. Grey Correlation Analysis

The changing trend of water consumption is affected by many factors. Grey correlation analysis is used to explore the correlation between research objects and influencing factors, and the calculation process is shown as follows.
(1) Set   T 0 = t 0 1 , t 0 2 , , t 0 n as the system characteristic behavior sequence, and T i = t i 1 , t i 2 , , t i n i = 1 , 2 , , m as the sequence of relevant factors. Calculate the mean terms of each sequence using s i k = t i p t ¯ i 1 i = 1 , 2 , , m . p = 1 , 2 , , n .
(2) Calculate the grey relational coefficients using Equation (1). For φ 0 , 1 ,
θ s 0 p , s i p = min i   min k s 0 p s i p + φ max i   max k s 0 p s i p s 0 p s i p + φ max i   max k s 0 p s i p .
(3) Calculate the grey relational degree using Equation (2).
θ T 0 , T i = 1 n k = 1 n θ s 0 p , s i p i = 1 , 2 , , m .

3.2. The DDGM(1,N) Model

In this paper, the DGM(1,N) with deformable accumulation (DDGM(1,N) model) is used to predict the water consumption (billion cubic meters) of 17 cities selected in the Yangtze River Delta Economic Zone. The process is introduced as follows.
(1) Assume that the system principal variable data sequence is Equation (3).
T 1 0 = t 1 0 1 , t 1 0 2 , , t 1 0 m .
The sequences of relevant factors are Equation (4).
T i 0 = t i 0 1 , t i 0 2 , , t i 0 m i = 2 , , N .
T i r is a sequence of cumulative generation operations of T i r . Consider that the definition of the deformable derivative (Wu and Zhao, 2019) r -order accumulation sequence is Equation (5).
T i r = t i r 1 , t i r 2 , , t i r m i = 1 , 2 , , N ,
The cumulative generation of the model is calculated using Equation (6).
t r p = t 0 p + r t 1 p 1 .
The particle swarm optimization algorithm is used to find the optimal order r ( r 0 , 1 ) in MATLAB (R2016b). r is the deformable accumulation parameter; it can adjust the weight of old and new data during sequence generation.
(2) The grey model DDGM(1,N) of the r -order is Equation (7).
t 1 r p = b 0 + b 1 t 1 r p 1 + i = 2 N b i t i r p .
b 0 is the grey function system parameter of the model; b 1 is the principal variable parameter. b 2 , , b N are grey action system parameters of different independent variables in the model.
Then, the least squares method is used to solve the parameters in the model (Equation (8)). Therefore, the parameters B and Y can be obtained by using Equations (9) and (10).
b ^ = b 0 , b 1 , b 2 , , b N T = B T B 1 B T Y ,
where
B = 1 t 1 r 1 t 2 r 2 t N r 2 1 t 1 r 2 t 2 r 3 t N r 3 1 t 1 r m 1 t 2 r m t N r m ,
Y = t 1 r 2 t 1 r 3 t 1 r m .
t ^ 1 r p = b 0 + b 1 t 1 r p 1 + i = 2 N b i t i r p
is recursive function of the DDGM(1,N) model. Then, the simulated value of t 1 r can be obtained by using Equation (11).
t ^ 0 p = t ^ 1 r p r t ^ 1 1 p 1 .
According to Equation (12), we can see that the simulation value of T 1 0 is T ^ 1 0 (see Equation (13)).
T ^ 1 0 = t ^ 1 0 1 , t ^ 1 0 2 , , t 1 0 m .
(4) The fitting error and prediction error are solved by calculating the mean absolute percentage error value (MAPE) according to Equation (14).
M A P E = 1 n p = 1 n t 1 0 p t ^ 1 0 p t 1 0 p 100 %

3.3. Properties of the Cumulative Generating Operator

Property 1. 
When  t 0 p + 1 > 1 r t 0 p , p = 1 , 2 , , n , the accumulation sequence is a monotonically increasing sequence. When  t 0 p + 1 < 1 r t 0 p , p = 1 , 2 , , n , the accumulation sequence is a monotonically decreasing sequence. Otherwise, the accumulation sequence is not monotonic.
Proof. 
Knowing that t r p = t 0 p + r i = 1 p 1 t 0 i ,
t r p + 1 = t 0 p + 1 + r i = 1 p t 0 i .
Then,
t r p + 1 t r p = t 0 p + 1 t 0 p + r t 0 p     = t 0 p + 1 + r 1 t 0 p     = t 0 p + 1 1 r t 0 p .
Therefore, when t 0 p + 1 > 1 r t 0 p , p = 1 , 2 , , n , the accumulation sequence is a monotonically increasing sequence; when t 0 p + 1 < 1 r t 0 p , p = 1 , 2 , , n , the accumulation sequence is a monotonically decreasing sequence. Otherwise, the accumulation sequence is not monotonic.
The proof of Property 1 is completed. □
Property 2. 
The larger the parameter  r , the larger the value obtained by accumulating the sequence.
Proof. 
By assuming that 0 < r 1 < r 2 < 1 , we can obtain the following formulas.
t r 1 p = t 0 p + r 1 i = 1 p 1 t 0 i ,
t r 2 p = t 0 p + r 2 i = 1 p 1 t 0 i .
t r 2 p t r 1 p = r 2 r 1 i = 1 p 1 t 0 i > 0 .
Therefore, t r 2 p > t r 1 p is true. This shows that the larger the parameter r , the larger the value obtained by accumulating the sequence.
The proof of Property 2 is completed. □

3.4. Analysis and Discussion of Model Stability

Lemma 1. 
Let  B C m × n , Δ B C n , Y C m × n , Δ Y C n ; the vector norm  ·  and matrix norm  ·  are compatible. If a matrix norm  ·  has  B 1 Δ B < 1 , then the solutions of the inhomogeneous linear equations  A x = Y  and  ( B + Δ B ) ( x + Δ x ) = Y + Δ Y  satisfy 
Δ x x κ γ Δ B B + Δ Y Y
where
κ = B B 1 , γ = 1 B B 1 Δ B B > 0
It can be known from Lemma 1 that matrix B and matrix Y are original matrices, and matrix Δ B and matrix Δ Y are perturbation matrices, which are perturbations of the solution. According to the matrix perturbation analysis, the perturbation bound of the solution can be obtained.
Theorem 1. 
Assume that the solution of the DDGM(1,N) model is  x . If only  x i 0 1  is perturbed by  ε i , then the perturbation bound of the solution in this case is
L x i 0 1 = ε κ γ m 2 r + 1 ε 1 + m 1 r i = 2 N ε i B + m 1 r ε 1 Y
Proof. 
If only x i 0 1 is perturbed by ε i , it is equivalent to when x ^ i ( 0 ) ( 1 ) = x i ( 0 ) ( 1 ) + ε i .
where
B ^ = B + Δ B = 1 x 1 r 1 x 2 r 2 x N r 2 1 x 1 r 2 x 2 r 3 x N r 3 1 x 1 r m 1 x 2 r m x N r m + 0 ε 1 r ε 2 r ε N 0 r ε 1 r ε 2 r ε N 0 r ε 1 r ε 2 r ε N
B ^ = B + Δ B = B + 0 ε 1 r ε 2 r ε N 0 r ε 1 r ε 2 r ε N 0 r ε 1 r ε 2 r ε N ,   Δ B = 0 ε 1 r ε 2 r ε N 0 r ε 1 r ε 2 r ε N 0 r ε 1 r ε 2 r ε N
It is obvious that we can obtain the value of Δ B .
Δ B = m 2 r + 1 ε 1 + m 1 r i = 2 N ε i
In the same way,
Y ^ = Y + Δ Y = x 1 r 2 x 1 r 3 x 1 r m + r ε 1 r ε 1 r ε 1
Y ^ = Y + Δ Y = Y + r ε 1 r ε 1 r ε 1 ,   Δ Y = r ε 1 r ε 1 r ε 1
We can also obtain the value of Δ Y .
Δ Y = m 1 r ε 1
Therefore, Δ B and Δ Y can be obtained.
According to Lemma 1, they can be obtained as follows:
Δ x x κ γ m 2 r + 1 ε 1 + m 1 r i = 2 N ε i B + m 1 r ε 1 Y
Then, the perturbation bound of the solution can be obtained.
L x i 0 1 = ε κ γ m 2 r + 1 ε 1 + m 1 r i = 2 N ε i B + m 1 r ε 1 Y
Theorem 2. 
If only  x ^ i 0 k = x i 0 k + ε i k = 2 , , m 1  is perturbed by  ε i , then the perturbation bound of the solution in this case is
L x i 0 k = κ γ m 3 r + 1 ε 1 + m 2 r + 1 i = 2 N ε i B + m 2 r + 1 ε 2 Y
Proof. 
If only x i 0 t is perturbed by ε i , it is equivalent to when x ^ i ( 0 ) ( t ) = x i ( 0 ) ( t ) + ε i , where
B ^ = B + Δ B = 1 x 1 r 1 x 2 r 2 x N r 2 1 x 1 r 2 x 2 r 3 x N r 3 1 x 1 r m 1 x 2 r m x N r m + 0 0 ε 2 ε N 0 ε 1 r ε 2 r ε N 0 r ε 1 r ε 2 r ε N
B ^ = B + Δ B = B + 0 0 ε 2 ε N 0 ε 1 r ε 2 r ε N 0 r ε 1 r ε 2 r ε N ,   Δ B = 0 0 ε 2 ε N 0 ε 1 r ε 2 r ε N 0 r ε 1 r ε 2 r ε N
The value of Δ B can be obtained.
Δ B = m 3 r + 1 ε 1 + m 2 r + 1 i = 2 N ε i
In the same way,
Y ^ = Y + Δ Y = x 1 r 2 x 1 r 3 x 1 r m + ε 2 r ε 2 r ε 2
Y ^ = Y + Δ Y = Y + ε 2 r ε 2 r ε 2 ,   Δ Y = ε 2 r ε 2 r ε 2
The value of Δ Y can be obtained.
Δ Y = m 2 r + 1 ε 2
Therefore, Δ B and Δ Y can be obtained.
According to Lemma 1, they can be obtained as follows:
Δ x x κ γ m 3 r + 1 ε 1 + m 2 r + 1 i = 2 N ε i B + m 2 r + 1 ε 2 Y
Then, the perturbation bound of the solution can be obtained.
L x i 0 k = κ γ m 3 r + 1 ε 1 + m 2 r + 1 i = 2 N ε i B + m 2 r + 1 ε 2 Y
Theorem 3. 
If only  x i 0 m  is perturbed by  ε i , then the perturbation bound of the solution in this case is
L x i 0 m = κ γ ε 2 + ε 3 + + ε N B + ε N Y
Proof. 
If only x i 0 m is perturbed by ε i , it is equivalent to when x ^ i ( 0 ) ( m ) = x i ( 0 ) ( m ) + ε i , where
B ^ = B + Δ B = 1 x 1 r 1 x 2 r 2 x N r 2 1 x 1 r 2 x 2 r 3 x N r 3 1 x 1 r m 1 x 2 r m x N r m + 0 0 0 0 0 0 0 0 0 0 ε 2 ε N
B ^ = B + Δ B = B + 0 0 0 0 0 0 0 0 0 0 ε 2 ε N ,   Δ B = 0 0 0 0 0 0 0 0 0 0 ε 2 ε N
The value of Δ B can be obtained.
Δ B = ε 2 + ε 3 + + ε N = i = 2 N ε i
In the same way,
Y ^ = Y + Δ Y = x 1 r 2 x 1 r 3 x 1 r m + 0 0 ε N
Y ^ = Y + Δ Y = Y + 0 0 ε N ,   Δ Y = 0 0 ε N
The value of Δ Y can be obtained.
Δ Y = ε N
Therefore, Δ B and Δ Y can be obtained.
According to Lemma 1, they can be obtained as follows:
Δ x x κ γ ε 2 + ε 3 + + ε N B + ε N Y
Then, the perturbation bound of the solution can be obtained.
L x i 0 m = κ γ ε 2 + ε 3 + + ε N B + ε N Y
According to the analysis and proof presented above:
The disturbance bound of the model L x i ( 0 ) ( k ) is correlated with the size of the sample m . In the case of equivalent disturbance, the larger the number of samples, the larger the disturbance bound of the model. This indicates that as the number of samples increases, the disturbance bound of the model increases and stability becomes worse, which means that the grey model is suitable for the modeling of the “less data” problem.
When the sample size n is unchanged, the particle swarm optimization algorithm can be used to adjust the parameter r , so that the disturbance bound of the solution satisfies L x i ( 0 ) ( 2 ) < L x i ( 0 ) ( 3 ) < < L x i ( 0 ) ( k 1 ) < L x i ( 0 ) ( k ) . That means x i ( 0 ) ( k ) is more sensitive to the solution than x i ( 0 ) ( k 1 ) . It also indicates that new information is given more weight than old information. Therefore, the DDGM(1,N) model conforms to the new information priority principle.

3.5. The Case Study

This paper takes the relationship between water consumption and different factors in Jinhua City as an example. The original data of Jinhua presented in Table 1 are from 2010 to 2017. The relationship between population and water consumption is analyzed according to the results, as well as the industrial added value and water consumption. Following the method described in Section 3.1 above, we can obtain the results that θ T 1 , T 2 = 0.8247 , θ T 1 , T 3 = 0.7595 ,   θ T 1 , T 4 = 0.6615 , θ T 1 , T 5 = 0.5347 . The results show that the population has the largest impact on the water consumption, followed by the added value of the primary industry and the secondary industry, and finally, the added value of the tertiary industry has the least impact on the water consumption.
DDGM(1,N) model is constructed according to the original data in Table 1, and the calculation process is shown as follows.
(1) The water consumption sequence of Jinhua City from 2010 to 2017 is
T 1 0 = 19.44 , 18.91 , 18.96 , 18.78 , 17.88 , 17.51 , 17.19 , 17.10
The population of Jinhua City from 2010 to 2017 is
T 2 0 = 536.2 , 538.6 , 539.9 , 542.8 , 543.7 , 545.4 , 552.0 , 556.4
Then, the optimal order r = 0.02 is obtained through the particle swarm optimization (PSO) algorithm in MATLAB (R2016b). The r -order r = 0.02 cumulative sequences are as follows:
T 1 0.02 = 19.44 , 19.30 , 19.73 , 19.93 , 19.41 , 19.40 , 19.44
T 2 0.02 = 536.2 , 549.40 , 561.55 , 575.33 , 587.16 , 599.82 , 617.40 , 632.92
(2)
b ^ = b 0 , b 1 , b 2 T = B T B 1 B T Y = 0.0744 0.0003 17.9512 T
According to Equations (9) and (10), the parameters B and Y can be obtained.
B = 19.4400 549.4017 1 19.3016 561.5518 1 19.7326 575.3281 1 19.9345 587.1628 1 19.4128 599.8156 1 19.4030 617.4026 1 19.4358 632.9227 1
Y = 19.3016 19.7326 19.9345 19.4128 19.4030 19.4358 19.6921
(3) The simulated value of T 1 0 is obtained using Equation (11) as follows:
T ^ 1 0 = 19.44 , 19.1478 , 18.7619 , 18.3694 , 17.976 , 17.5826 , 17.1906 , 16.7981 , 16.4015 , 16.0077
Finally, the fitting and prediction results can be obtained, and the validity can be tested using MAPE with fitting MAPE = 0.9022% and prediction MAPE = 1.36%, respectively.
The fitting and prediction results of DDGM(1,2), GM(1,1) and GM(1,2) are shown in Table 2. Compared with the two classical grey models, the fitting MAPE and prediction MAPE of DDGM(1,2) model are both the smallest. There is no significant difference between the fitting MAPE of DDGM(1,2) and GM(1,1), but the prediction MAPE of DDGM(1,2) is smaller than that of GM(1,1). As for GM(1,2), both the fitting MAPE and the prediction MAPE of DDGM(1,2) are much smaller than GM(1,2). Therefore, DDGM(1,2) model has better fitting and prediction performance, indicating that this model is more effective in predicting water consumption.
Similarly, the DDGM(1,2) model is constructed according to the water consumption and the added value of the tertiary industry shown in Table 1. The fitting and prediction results of the three models are shown in Table 3. It can be found that the fitting MAPE and prediction MAPE of the DDGM(1,2) model are smaller than those of GM(1,1) and GM(1,2). This further illustrates the effectiveness of the DDGM(1,2) model to predict water consumption.

4. Forecasting Results and Discussion

In order to accurately predict the water consumption, we should consider the impact of three major industries and the population. However, the future changes of these variables are unknown. Therefore, this paper makes a hypothesis on the growth rate based on the government work report and urban development planning from recent years. Combined with different influencing factors, the DDGM(1,N) model was established to predict the urban water consumption. The calculation steps are the same as those in Section 3.2 in this paper. In this section, we predict the water consumption in the years of 2020–2025 using data from 2010 to 2019 for cities within Zhejiang Province and Shanghai. Due to the incompleteness of the data obtained in Anhui Province, the original data we had are from 2013 to 2019, so the forecasted range of water consumption in Anhui Province covers the years of 2020–2023.

4.1. Discussion on the Forecasting Results Combined with Population

The forecasting results for water consumption in the cities of Zhejiang Province are shown in Table 4 and Figure 2. As shown in Figure 2, the water consumption in Zhoushan will have an upward trend. However, the change in water consumption in Zhoushan is not significant compared with other cities. In the other seven cities, the water consumption shows a descending trend with the rise in population. In terms of numerical analysis, Hangzhou and Ningbo have a comparatively large population, so water consumption will also be at a high level. Water consumption in Huzhou, Jiaxing, Shaoxing, Jinhua and Taizhou will be at a moderate level. Due to the continuous increase in population and the increase in per capita income, personal and family consumption continue to increase, which has resulted in an increasing discharge of urban domestic sewage. Therefore, the government should adopt a series of control measures such as a tiered pricing mechanism for the management of water prices.
At present, many cities in the Yangtze River Delta have become water-scarce cities. The conservation and protection of water resources has become an urgent matter. The government should aim to improve water quality and protect the water system by establishing a water pollution prevention mechanism. It is strictly forbidden to develop high-risk and high-polluting industries at the source of rivers and upstream. At the same time, the population of these cities should be taken seriously. The government should also actively promote water conservation knowledge and carry out thematic activities related to the protection of resources so they can finally achieve the goal of improving people’s water conservation awareness.
Table 5 and Figure 3 show the prediction results for Shanghai at different growth rates of population. With the rise in the population growth rate, the water consumption in Shanghai will decline. During the process of rapid industrialization and urbanization, the influence of environmental factors on regional economic development is increasingly apparent. With the continuous emission and accumulation of pollutants from production activities and people’s daily lives, Shanghai is facing the pressure of regional environmental quality. At the same time, Shanghai as a densely populated city, and the problem of population control also brings challenges to the government’s management of water resources. Taihu Lake, the lower reaches of the Yangtze River, Qiantang River and other water resources are polluted in different degrees, while Taihu Lake is the most serious. At present, only 1% of the surface water in Shanghai meets the national standard of drinking water sources. Water pollution makes this area, which has abundant water resources, gradually become a water shortage region.
The forecasting results for Anhui Province under different population growth rates are shown in Table 6 and Figure 4. The water consumption in Anqing, Chuzhou, Hefei and Maanshan in Anhui Province decrease with the increase in population. This shows that under the correct guidance of these city governments and with the reasonable and effective control of water resources, the water conservation awareness of residents and the protection of resources improve. With the change in the population growth rate, the tendency of water consumption in Xuancheng inclines to present a consistently downward trend. This shows that the effect of demographic factors on the changes in water consumption in Xuancheng is not obvious, and it can be analyzed from other aspects. When the population shows a negative growth, the water consumption in Tongling declines. When the population shows a positive growth trend, the water consumption in Tongling begins to show an upward trend. This indicates that the population factor has an impact on the change in water consumption. So, we should pay attention to the population management of the city to achieve better results in the management of water resources.

4.2. Discussion on the Forecasting Results Combined with the Primary Industry

The forecasting results for the cities in Zhejiang Province under different growth rates of the added value of the primary industry are shown in Table 7 and Figure 5. The figure shows that the prediction results for six cities will descend with the rise in the added value of the primary industry. On the one hand, people’s water conservation awareness in these areas was improved, and the water-saving technology industry was improved. On the other hand, due to the growth of the primary industry, the government has strengthened the control and protection of water resources. The trend of water consumption in Jiaxing and Zhoushan is different from that in other cities. The faster the primary industry develops, the greater the annual water consumption. However, compared with other cities, the population base of Zhoushan is relatively small, and the expansion of the primary industry is relatively backward. Therefore, the change in the water consumption is not obvious.
Table 8 and Figure 6 show the prediction results for Shanghai. With the accelerated growth in the primary industry, the water consumption in the city has shown a downward trend. This shows that the government’s management is effective. The agriculture and animal husbandry industries that emit large amounts of pollutants are carefully governed to reduce the pressure on the environment. Shanghai is a leader in the economic integration of the Yangtze River Delta and, accompanied with the economic transformation and construction of an international metropolis, Shanghai has a huge spillover effect on the development of the surrounding cities. In terms of water management and control, Shanghai should continue to maintain its current momentum. Agriculture that pollutes the environment should continue to be paid attention to and given management measures. Furthermore, Shanghai should promote the construction of a water-saving society in an all-around and in-depth manner.
The prediction results for water consumption in the cities of Anhui Province are shown in Table 9 and Figure 7. When the primary industry increases, the annual water consumption in Chuzhou, Hefei and Ma’anshan will show a downward trend. The greater the growth rate of the primary industry is, the faster the reduction in water consumption will be. The change trend of water consumption in Tongling and Wuhu is very similar. A faster growth rate of water consumption will emerge when the growth rate of the primary industry is higher. Water consumption in Chizhou shows an upward trend when the assumed growth rate increases. This shows that there is an important relationship between Chizhou’s primary industry and water consumption. When the growth rate changes, Xuancheng will always show a downward trend. Anqing’s annual water consumption has shown a downward trend. Therefore, the cities should continue to strengthen the evolution and upgrading of technology on the preservation of the environment, which ultimately helps to achieve the aim of improving efficiency.

4.3. Discussion on the Forecasting Results Combined with the Secondary Industry

The forecasting results for the cities in Zhejiang Province under different change rates of the added value of the secondary industry are shown in Table 10 and Figure 8. With the rise in the added value of the secondary industry, the annual water consumption in cities except for Zhoushan all show a downward trend. This shows that under the guidance of the government’s intervention in the secondary industry, good water-saving results were achieved in water resources management. For example, Hangzhou issued a “water limit order” to 212 companies that use more than 10,000 tons of water per month. Shaoxing launched a water-saving campaign for the eight major water-consuming industries, including thermal power, chemical industry, and papermaking; other cities adopted various water-saving measures. It is recommended that the water-saving management departments in various cities should strengthen their communication and exchanges of experience in building a water-saving society. Therefore, we should jointly seek the key points of conservation in cities with “water-quality water shortage”.
The water consumption in Zhoushan has shown an upward trend, but the base and change in water consumption are very small. This is because compared with other cities, Zhoushan’s industry and construction are not well developed, and the speed of economic growth is comparatively slow, which is limited by the geographical conditions of Zhoushan. Zhoushan is composed of many islands, so it faces many problems in terms of transportation.
Table 11 and Figure 9 show the forecasting results for Shanghai under different growth rates of the secondary industry added value. The faster the secondary industry grows, the faster the water consumption will decrease. This shows that the efforts made in Shanghai to promote the construction of a water-saving society have achieved excellent results. The government implements and improves a strict water resources management system to fundamentally facilitate water conservation. The water resources management department conducts regular inspections of water-using enterprises. Through the formulation and implementation of water-saving measures, the water-saving potential of enterprises and industrial water-saving technologies can be tapped. Meanwhile, through the innovation and application of technology, the recycling utilization rate of industrial water was improved.
The forecasting results for the cities in Anhui Province under different growth rates of the secondary industry added value are shown in Table 12 and Figure 10. The results show that the water consumption in Anqing, Ma’anshan and Xuancheng decrease with the rise in the added value of the secondary industry. This is because these cities have high industrial water utilization and advanced industrial water-saving technologies. The change trends of water consumption in Tongling and Wuhu are very similar. Water consumption will decline with the development of the secondary industry in these cities. Water consumption in Chizhou will rise with the increase in the growth rate. This shows that with the expansion of the city’s size, the scale and number of industries have increased, which further led to an increase in industrial water consumption, which, in turn, triggered the rise in Chizhou’s water consumption. With the development of the secondary industry added value, the water consumption in Hefei will show a downward trend because Hefei increased the utilization rate of industrial water by improving water-saving technologies. Therefore, the local government should strengthen the management and control of industrial water and increase the utilization rate of industrial water.

4.4. Discussion on the Forecasting Results Combined with the Tertiary Industry

The forecasting results for the cities in Zhejiang Province are shown in Table 13 and Figure 11. The results show that except for Zhoushan, the water consumption in the other seven cities present a downward trend because these cities pay more attention to water management, e.g., in the technical transformation and upgrading of water-using equipment for industries with high water consumption such as accommodation, catering and real estate. At the same time, the adjustment of water prices can also help to achieve the goal of improving water efficiency. With the development of the tertiary industry, Zhoushan’s annual water consumption has shown a gentle upward trend. Compared with other cities, Zhoushan’s economic development is lagging, and the growth of the tertiary industry is also at a low level. Therefore, the water consumption in Zhoushan is kept at a low level.
Table 14 and Figure 12 show the forecasting results for Shanghai under different growth rates of the tertiary industry. It is shown that a downtrend will present in the next few years. This is because Shanghai has strengthened their management of water resources by launching a number of water-saving campaigns to raise people’s awareness of the importance of saving water. Regarding the remediation of industries with a high water consumption, the water resources management department has carried out measures by adjusting water prices, updating water-saving equipment and upgrading water-saving technologies.
The forecasting results for water consumption in the cities of Anhui Province are shown in Table 15 and Figure 13. We can find that the trends of water consumption in Chizhou, Chuzhou, Hefei and Tongling are very similar. The rise in water consumption accelerates as the growth rate of the tertiary industry increases. It is suggested that these cities should control the growth speed of the tertiary industry so as to avoid excessively fast economic growth from putting a huge pressure on the resources and environment. Traditional high-polluting industries such as the printing and dyeing, papermaking and chemical industries still account for a relatively high share, while the proportion of high-tech industries and service industries is still relatively small in these cities. The annual water consumption in Anqing, Ma’anshan, Wuhu and Xuancheng decrease. The government departments of these cities are actively carrying out water-saving-related activities. The water efficiency has been improved by raising water prices and improving reclaimed water utilization technology. Through these methods, the coordinated development between water use and economic development can be achieved.

5. Conclusions

This article firstly uses the grey relational analysis to study the related factors affecting water consumption. Then, through the model comparisons, we found that the DDGM(1,N) model has a better fitting and prediction effect. Finally, the DDGM (1,N) model was used to predict the water consumption in 17 cities selected in the Yangtze River Delta region by considering different influencing factors. The results show that as the growth rate of the population increases, the water consumption in 15 cities show a downward trend, and 2 cities show an upward trend. With the increase in the growth rate of the primary industry, the water consumption in 12 cities show a downward trend, and 5 cities show an upward trend. With the rise in the growth rate of the secondary industry, the water consumption in 14 cities show a downward trend, and 3 cities show an upward trend. With the rise in the growth rate of the tertiary industry, the water consumption in 12 cities show a downward trend, and 5 cities show an upward trend. In general, the water consumption in most cities will continue to show a downward trend in the future.
According to the prediction results obtained in this paper, we have some policy recommendations as follows. First of all, urban populations should be controlled. We should avoid a situation where the population density is concentrated, which will cause great pressure on the water consumption. Secondly, for cities where water consumption has increased with the development of the primary industry, the water resource utilization efficiency of agriculture and animal husbandry should be improved. Moreover, for cities where the water consumption increases with the development of the secondary industry, it is suggested that they adjust the internal structure of the secondary industry to reduce the proportion of high water-consuming industries, and meanwhile, strengthen the supervision of industrial water efficiency. It is possible to raise the utilization rate by adjusting the water prices and upgrading the water-saving technologies of water-using equipment. Finally, for cities whose water consumption has increased with the development of the tertiary industry, the service industries of these cities should shift to the direction of eco-environmental protection. We should focus on the accommodation and catering industry, carry out technical transformation of its water equipment and eliminate outdated high water-consuming equipment simultaneously. Through the government’s management and policy guidance, the goal of improving the utilization effect of water resources is realized.

Author Contributions

Conceptualization, Z.Q. and W.Y.; methodology, Z.Q.; software, Z.Q.; validation, Z.Q., W.Y.; formal analysis, Z.Q.; investigation, Z.Q.; resources, Z.Q.; data curation, Z.Q.; writing—original draft preparation, Z.Q.; writing—review and editing, W.Y.; visualization, Z.Q.; supervision, W.Y.; project administration, W.Y.; funding acquisition, W.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (NSFC) under grant No. 71971131; the MOE (Ministry of Education in China) Youth Project of Humanities and Social Sciences Fund under grant number 17YJC630121; the Research Project Supported by Shanxi Scholarship Council of China under grant No. 2022–027. And The APC was funded by Wei Yang.

Data Availability Statement

The data sources are the Statistical Yearbook of local statistics web-sites. The data of Shanghai, Zhejiang and Anhui are from the official website of the Municipal Bureau of Statistics (http://tjj.sh.gov.cn/, http://tjj.zj.gov.cn/, http://tjj.ah.gov.cn/) (accessed on 15 July 2022).

Acknowledgments

The relevant research is supported by the National Natural Science Foundation of China (NSFC) under grant No. 71971131; the MOE (Ministry of Education in China) Youth Project of Humanities and Social Sciences Fund under grant number 17YJC630121; the Research Project Supported by Shanxi Scholarship Council of China under grant No. 2022−027.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fan, J.-L.; Wang, J.-D.; Zhang, X.; Kong, L.-S.; Song, Q.-Y. Exploring the changes and driving forces of water footprints in China from 2002 to 2012: A perspective of final demand. Sci. Total Environ. 2019, 650, 1101–1111. [Google Scholar] [CrossRef] [PubMed]
  2. Li, S.; Han, W. Performance evaluation for urban water supply services in China. Water Supply 2020, 20, 3511–3516. [Google Scholar] [CrossRef]
  3. Bata, M.H.; Carriveau, R.; Ting, D.S.-K. Short-Term Water Demand Forecasting Using Nonlinear Autoregressive Artificial Neural Networks. J. Water Resour. Plan. Manag. 2020, 146, 04020008. [Google Scholar] [CrossRef]
  4. Mensik, P.; Marton, D. Hybrid Optimization Method for Strategic Control of Water Withdrawal from Water Reservoir with Using Support Vector Machines. Procedia Eng. 2017, 186, 491–498. [Google Scholar] [CrossRef]
  5. Tarebari, H.; Javid, A.H.; Mirbagheri, S.A.; Fahmi, H. Multi-Objective Surface Water Resource Management Considering Conflict Resolution and Utility Function Optimization. Water Resour. Manag. 2018, 32, 4487–4509. [Google Scholar] [CrossRef]
  6. Vonk, E.; Cirkel, D.G.; Blokker, M. Estimating Peak Daily Water Demand under Different Climate Change and Vacation Scenarios. Water 2019, 11, 1874. [Google Scholar] [CrossRef] [Green Version]
  7. Pei, X.; Sun, Y.; Ren, Y. Demand estimation of water resources via bat algorithm. Int. J. Wirel. Mob. Comput. 2020, 18, 16–21. [Google Scholar] [CrossRef]
  8. Shao, Z.; Wu, F.; Li, F.; Zhao, Y.; Xu, X. System Dynamics Model for Evaluating Socio-Economic Impacts of Different Water Diversion Quantity from Transboundary River Basins—A Case Study of Xinjiang. Int. J. Environ. Res. Public Health 2020, 17, 9091. [Google Scholar] [CrossRef]
  9. Mu, L.; Zheng, F.; Tao, R.; Zhang, Q.; Zoran Kapelan, Z. Hourly and Daily Urban Water Demand Predictions Using a Long Short-Term Memory Based Model. J. Water Resour. Plan. Manag. 2020, 146, 05020017. [Google Scholar] [CrossRef]
  10. Araujo, W.C.; Esquerre, K.P.O.; Sahin, O.; Silva, B.B.S. Scarcity Pricing for Water Demand Management: A Case Study in Campina Grande-Brazil. J. Environ. Account. Manag. 2021, 9, 75–91. [Google Scholar] [CrossRef]
  11. Anderson, B.; Manouseli, D.; Nagarajan, M. Estimating scenarios for domestic water demand under drought conditions in England and Wales. Water Sci. Technol. Water Supply 2018, 18, 2100–2107. [Google Scholar] [CrossRef] [Green Version]
  12. Gibson, J.; Karney, B.; Guo, Y. Effects of Demand, Mixing Fraction, and Rate Coefficient Uncertainty on Water Quality Models. J. Water Resour. Plan. Manag. 2020, 146, 06020004. [Google Scholar] [CrossRef]
  13. Massoud, B.; Reza, K. Evaluating the resilience of water resources management scenarios using the evidential reasoning approach: The Zarrinehrud river basin experience. J. Environ. Manag. 2021, 284, 112025. [Google Scholar] [CrossRef]
  14. Maria, X.; Chris, H.; Jan, H.; Zoran, K. Short-Term Forecasting of Household Water Demand in the UK Using an Interpretable Machine-Learning Approach. J. Water Resour. Plan. Manag. 2021, 147, 04021004. [Google Scholar] [CrossRef]
  15. Berthe, Y.A.; Jules, M.O.M.; Sampah, G.E.; Alioune, K.; Albert, G.B.T. Using “Water Evaluation and Planning” (WEAP) Model to Simulate Water Demand in Lobo Watershed (Central-Western Cote d’Ivoire). J. Water Resour. Prot. 2021, 13, 216–235. [Google Scholar] [CrossRef]
  16. Tu, L.; Chen, Y. An unequal adjacent grey forecasting air pollution urban model. Appl. Math. Model. 2021, 99, 260–275. [Google Scholar] [CrossRef]
  17. Rees, P.L.S. Advancing Agricultural Water Security and Resilience Under Nonstationarity and Uncertainty: Evolving Roles of Blue, Green, and Grey Water. J. Contemp. Water Res. Educ. 2018, 165, 1–3. [Google Scholar] [CrossRef]
  18. Kose, E.; Tasci, L. Geodetic deformation forecasting based on multi-variable grey prediction model and regression model. Grey Syst. Theory Appl. 2019, 9, 464–471. [Google Scholar] [CrossRef]
  19. Jamshidi, S. An approach to develop grey water footprint accounting. Ecol. Indic. 2019, 106, 105477. [Google Scholar] [CrossRef]
  20. Ma, D.; Duan, H.; Li, W.; Zhang, J.; Liu, W.; Zhou, Z. Prediction of water inflow from fault by particle swarm optimization-based modified grey models. Environ. Sci. Pollut. Res. Int. 2020, 27, 42051–42063. [Google Scholar] [CrossRef]
  21. Li, B.; Wu, Q.; Zhang, W.; Liu, Z. Water resources security evaluation model based on grey relational analysis and analytic network process: A case study of Guizhou Province. J. Water Process Eng. 2020, 37, 101429. [Google Scholar] [CrossRef]
  22. Zhang, K.; Zhong, Q.; Qu, P.; Yin, Y.; Zuo, Y. Rural water environment quality grey prediction model based on network search information. Chin. Manag. Sci. 2020, 28, 222–230. [Google Scholar]
  23. Liao, X.; Chai, L.; Liang, Y. Income impacts on household consumption’s grey water footprint in China. Sci. Total Environ. 2021, 755, 142584. [Google Scholar] [CrossRef] [PubMed]
  24. Ju, Q.; Hu, Y. Source identification of mine water inrush based on principal component analysis and grey situation decision. Environ. Earth Sci. 2021, 80, 157. [Google Scholar] [CrossRef]
  25. Meng, X.; Wu, L. Prediction of per capita water consumption for 31 regions in China. Environ. Sci. Pollut. Res. 2021, 28, 29253–29264. [Google Scholar] [CrossRef]
  26. Deng, J. Control problems of grey systems. Syst. Control Lett. 1982, 1, 288–294. [Google Scholar] [CrossRef]
  27. Guo, J.; Tu, L.; Qiao, Z.; Wu, L. Forecasting the air quality in 18 cities of Henan Province by the compound accumulative grey model. J. Clean. Prod. 2021, 310, 127582. [Google Scholar] [CrossRef]
  28. Wang, Z.; Dang, Y.; Liu, S. Optimization of Background Value in GM(1,1) Model. Syst. Eng. Theory Pract. Online 2008, 28, 61–67. [Google Scholar] [CrossRef]
  29. Ma, X.; Liu, Z. Application of a novel time-delayed polynomial grey model to predict the natural gas consumption in China. J. Comput. Appl. Math. 2017, 324, 17–24. [Google Scholar] [CrossRef]
  30. Zeng, B.; Duan, H.; Bai, Y.; Meng, W. Forecasting the output of shale gas in China using an unbiased grey model and weakening buffer operator. Energy 2018, 151, 238–249. [Google Scholar] [CrossRef]
  31. Chen, Y.; Yang, Y. A Multivariate Statistical Model of Water Leakage in Urban Water Supply Networks Based on Random Matrix Theory. Math. Probl. Eng. 2022, 2022, 2314972. [Google Scholar] [CrossRef]
  32. Li, M.; Cheng, Z.; Lin, W.; Wei, Y.; Wang, S. What can be learned from the historical trend of crude oil prices? An ensemble approach for crude oil price forecasting. Energy Econ. 2023, 123, 106736. [Google Scholar] [CrossRef]
  33. Liu, L.; Chen, Y.; Wu, L. The damping accumulated grey model and its application. Commun. Nonlinear Sci. Numer. Simul. 2021, 95, 105665. [Google Scholar] [CrossRef]
  34. Duan, H.; Pang, X. A novel grey prediction model with system structure based on energy background: A case study of Chinese electricity. J. Clean. Prod. 2023, 390, 136099. [Google Scholar] [CrossRef]
  35. Wu, L.; Zhao, H. Discrete grey model with the weighted accumulation. Soft Comput. 2019, 23, 12873–12881. [Google Scholar] [CrossRef]
  36. Li, Y.; Li, J. Study on unbiased interval grey number prediction model with new information priority. Grey Syst. Theory Appl. 2019, 10, 1–11. [Google Scholar] [CrossRef]
  37. Huynh, N.T.; Nguyen, T.V.T.; Tam, N.T.; Nguyen, Q.M. Optimizing Magnification Ratio for the Flexible Hinge Displacement Amplifier Mechanism Design. In Proceedings of the 2nd Annual International Conference on Material, Machines and Methods for Sustainable Development (MMMS2020), MMMS 2020, Lecture Notes in Mechanical Engineering, Nha Trang, Vietnam, 12–15 November 2020; Long, B.T., Kim, Y.H., Ishizaki, K., Toan, N.D., Parinov, I.A., Vu, N.P., Eds.; Springer: Cham, Switzerland, 2021; pp. 769–778. [Google Scholar] [CrossRef]
Figure 1. The geographical location of the study area.
Figure 1. The geographical location of the study area.
Axioms 12 00655 g001
Figure 2. Forecasting results for Zhejiang Province.
Figure 2. Forecasting results for Zhejiang Province.
Axioms 12 00655 g002aAxioms 12 00655 g002b
Figure 3. Forecasting results for Shanghai.
Figure 3. Forecasting results for Shanghai.
Axioms 12 00655 g003
Figure 4. Forecasting results for Anhui Province.
Figure 4. Forecasting results for Anhui Province.
Axioms 12 00655 g004aAxioms 12 00655 g004b
Figure 5. Forecasting results for Zhejiang Province.
Figure 5. Forecasting results for Zhejiang Province.
Axioms 12 00655 g005aAxioms 12 00655 g005b
Figure 6. Forecasting results for Shanghai.
Figure 6. Forecasting results for Shanghai.
Axioms 12 00655 g006
Figure 7. Forecasting results for Anhui Province.
Figure 7. Forecasting results for Anhui Province.
Axioms 12 00655 g007aAxioms 12 00655 g007b
Figure 8. Forecasting results for Zhejiang Province.
Figure 8. Forecasting results for Zhejiang Province.
Axioms 12 00655 g008aAxioms 12 00655 g008b
Figure 9. Forecasting results for water consumption in Shanghai.
Figure 9. Forecasting results for water consumption in Shanghai.
Axioms 12 00655 g009
Figure 10. Forecasting results for water consumption in Anhui Province.
Figure 10. Forecasting results for water consumption in Anhui Province.
Axioms 12 00655 g010
Figure 11. Forecasting results for Zhejiang Province.
Figure 11. Forecasting results for Zhejiang Province.
Axioms 12 00655 g011
Figure 12. Forecasting results for Shanghai.
Figure 12. Forecasting results for Shanghai.
Axioms 12 00655 g012
Figure 13. Forecasting results for Anhui Province.
Figure 13. Forecasting results for Anhui Province.
Axioms 12 00655 g013aAxioms 12 00655 g013b
Table 1. Data on water consumption and related factors in Jinhua City from 2010 to 2017.
Table 1. Data on water consumption and related factors in Jinhua City from 2010 to 2017.
YearWater Consumption
(Billion Cubic Meters)
Population
(10,000 People)
Value Added by the Primary Industry
(CNY 100 Million)
Value Added by the Secondary Industry
(CNY 100 Million)
Value Added by the Tertiary Industry
(CNY 100 Million)
201019.44536.2104.811054.64941.68
201118.91538.6120.661180.601131.85
201218.96539.9128.361285.131276.36
201318.78542.8128.991402.261449.69
201417.88543.7130.391482.061602.72
201517.51545.4131.901550.901798.28
201617.19552.0137.731561.292006.21
201717.10556.4134.971622.742204.45
Table 2. Comparison results between three models based on the population factor (units: billion cubic meters).
Table 2. Comparison results between three models based on the population factor (units: billion cubic meters).
YearActual ValueDDGM(1,2)GM(1,1)GM(1,2)
201019.4419.4419.4419.44
201118.9119.1519.1614.70
201218.9618.7618.7818.97
201318.7818.3718.418.15
201417.8817.9818.0317.99
201517.5117.5817.6718.02
201617.1917.1917.3218.24
201717.1016.8016.9718.39
Simulation MAPE (%) 0.9021.085.36
201816.2816.4117.4118.52
201915.7116.0116.8218.58
Forecast MAPE (%) 1.367.0016.01
Table 3. Comparison results between three models based on the tertiary industry factor (units: billion cubic meters).
Table 3. Comparison results between three models based on the tertiary industry factor (units: billion cubic meters).
YearActual ValueDDGM(1,2)GM(1,1)GM(1,2)
201019.4419.4419.4419.44
201118.9119.1619.1614.70
201218.9618.7718.7818.97
201318.7818.3718.418.15
201417.8817.9818.0317.99
201517.5117.5817.6718.02
201617.1917.1917.3218.24
201717.1016.8016.9718.39
Simulation MAPE (%) 0.9031.085.36
201816.2816.4017.4118.52
201915.7116.0016.8218.58
Forecast MAPE (%) 1.317.0016.01
Table 4. Forecasting results at different growth rates of population in Zhejiang (units: billion cubic meters).
Table 4. Forecasting results at different growth rates of population in Zhejiang (units: billion cubic meters).
CityGrowth Rate202020212022202320242025
Hangzhou1.0%23.5620.2917.5015.0512.8210.75
2.0%23.1019.0515.2211.567.984.46
3.0%22.6617.8012.928.003.01−2.06
Ningbo1.0%20.2620.1219.9919.8719.7619.65
3.0%20.2219.9819.7319.4719.2118.93
5.0%20.1419.7919.4018.9918.5518.08
Huzhou0.5%12.0811.2510.449.658.888.11
1.0%11.9710.949.858.707.516.28
1.5%11.8610.629.247.746.134.42
Jiaxing1.0%17.8217.5617.3117.0516.8016.54
2.0%17.6517.1716.6616.1515.6215.08
3.0%17.4916.7816.0315.2414.4213.58
Shaoxing0.5%17.1216.6116.0715.5114.9314.34
1.0%17.2516.6015.8715.0614.1813.24
1.5%16.9716.1615.2114.1312.9211.59
Jinhua0.5%15.2314.7214.2013.6613.1212.57
1.0%15.1914.5813.9313.2212.4811.68
1.5%15.1414.4513.6612.7911.8310.79
Zhoushan0.5%1.621.641.661.681.691.71
1.0%1.631.661.691.721.751.78
1.5%1.641.671.711.761.811.86
Taizhou0.5%14.7514.2313.6713.0612.4211.74
1.0%14.6613.9813.1912.2811.2610.15
1.5%14.5813.7412.7011.4810.098.53
Table 5. Forecasting results at different growth rates of population in Shanghai (units: billion cubic meters).
Table 5. Forecasting results at different growth rates of population in Shanghai (units: billion cubic meters).
CityGrowth Rate202020212022202320242025
Shanghai1.0%101.89100.8699.7198.4597.1395.74
3.0%101.1598.8796.0392.7489.0985.12
5.0%100.4896.9292.3286.8480.6273.74
Table 6. Forecasting results at different growth rates of population in Anhui (units: billion cubic meters).
Table 6. Forecasting results at different growth rates of population in Anhui (units: billion cubic meters).
CityGrowth Rate2020202120222023
Hefei0.5%30.4930.4430.3930.34
1.0%30.4830.4230.3630.30
1.5%30.4730.4030.3230.25
Tongling−0.5%12.6012.5512.4412.30
0.5%12.7212.8713.0013.13
1.0%12.7913.0413.3013.57
Chizhou−0.5%10.1210.1210.1210.12
0.5%10.1010.0910.0810.07
1.0%10.1010.0810.0610.04
Wuhu−0.5%29.9729.9729.9629.96
0.5%29.9629.9529.9429.94
1.0%29.9529.9429.93329.92
Anqing−0.5%23.6323.1122.5021.97
0.5%23.6823.1722.6022.14
1.0%23.7023.2022.6522.14
Xuancheng−0.5%14.1514.2214.1014.13
0.5%14.1314.1914.0714.08
1.0%14.1214.1814.0514.06
Maanshan−0.5%31.3731.1130.8930.71
0.5%31.2530.7830.3229.87
1.0%31.1830.6230.0329.44
Chuzhou−0.5%22.4222.4122.4022.39
0.5%22.4022.3722.3522.32
1.0%22.3822.3522.3222.28
Table 7. Forecasting results at different growth rates of the primary industry in Zhejiang (units: billion cubic meters).
Table 7. Forecasting results at different growth rates of the primary industry in Zhejiang (units: billion cubic meters).
CityGrowth Rate202020212022202320242025
Hangzhou5%26.6923.4820.0616.4612.678.69
10%25.6420.5514.487.50−0.40−9.21
15%24.6517.588.56−2.42−15.47−30.74
Ningbo5%20.3120.1720.0319.8719.7119.53
10%20.2820.0619.8019.5119.1818.82
15%20.2219.9119.5219.0618.5317.91
Huzhou5%10.539.097.415.483.260.72
10%10.468.866.914.581.80−1.50
15%10.398.626.383.590.14−4.09
Jiaxing−5%18.7718.7318.7018.6618.6318.59
5%18.7918.8018.8318.8618.9018.94
10%18.8018.8418.9018.9819.0719.18
Shaoxing5%17.2416.6916.0515.3414.5613.70
10%17.0716.2015.0813.7312.1310.29
15%16.9115.6914.0411.939.346.21
Jinhua5%15.2014.5813.8913.1412.3211.43
10%15.1114.3213.3712.2310.909.35
15%15.0314.0612.8011.229.286.91
Zhoushan5%1.641.671.691.721.741.77
10%1.661.711.761.821.881.95
15%1.681.751.831.932.052.18
Taizhou5%13.9613.0912.1411.119.998.80
10%13.8312.6811.309.677.775.60
15%13.6912.2510.408.085.241.82
Table 8. Forecasting results at different growth rates of the primary industry in Shanghai (units: billion cubic meters).
Table 8. Forecasting results at different growth rates of the primary industry in Shanghai (units: billion cubic meters).
CityGrowth Rate202020212022202320242025
Shanghai−5%106.75108.42110.42112.70115.20117.89
5%105.68105.31104.39102.93100.9698.48
10%105.14103.67101.0597.2492.2185.93
Table 9. Forecasting results at different growth rates of the primary industry in Anhui (units: billion cubic meters).
Table 9. Forecasting results at different growth rates of the primary industry in Anhui (units: billion cubic meters).
CityGrowth Rate2020202120222023
Hefei−5%30.4430.4630.4830.50
5%30.2630.1129.9729.81
10%30.1729.9229.6629.38
Tongling−5%12.9512.6412.3111.98
5%13.2213.2813.3413.42
10%13.3613.6213.9314.29
Chizhou5%10.1210.1210.1110.11
10%10.1510.1910.2310.27
15%10.1810.2610.3510.46
Wuhu−5%30.0529.9229.7929.66
5%30.3430.4530.5630.68
10%30.4930.7431.0231.33
Anqing−5%23.7723.4822.5922.26
5%24.2323.9623.4923.22
10%24.4624.2524.0223.85
Xuancheng−5%14.1814.2114.1514.16
5%14.1114.1314.0214.00
10%14.0914.1013.9513.91
Maanshan−5%30.1729.7829.5229.34
5%29.6528.4027.0825.70
10%29.3927.6825.7223.54
Chuzhou−5%22.4322.4722.5722.65
5%22.1522.0021.8921.76
10%22.0021.7521.4921.19
Table 10. Forecasting results at different growth rates of the secondary industry in Zhejiang (units: billion cubic meters).
Table 10. Forecasting results at different growth rates of the secondary industry in Zhejiang (units: billion cubic meters).
CityGrowth Rate202020212022202320242025
Hangzhou5%27.2824.4721.5218.4215.1811.78
10%26.4622.1016.9511.044.38−3.06
15%25.6219.6412.042.81−8.15−20.97
Ningbo5%20.2920.1620.0219.8819.7419.59
10%20.2420.0319.7919.5219.2218.90
15%20.2019.9019.5419.1218.6318.06
Huzhou5%12.1911.4410.719.979.228.45
10%12.0310.999.838.537.105.51
15%11.8810.538.896.944.651.97
Jiaxing5%17.5617.3117.0716.8316.5916.33
10%17.4517.0316.5616.0615.5114.90
15%17.3416.7316.0115.1914.2313.14
Shaoxing5%16.9016.2815.6314.9614.2513.52
10%16.7215.7614.6413.3511.9010.27
15%16.5415.2313.5811.569.176.36
Jinhua5%15.3314.8214.2913.7313.1512.54
10%15.2214.5113.6712.7011.6010.37
15%15.1214.1913.0111.569.837.78
Zhoushan5%1.681.721.751.791.821.86
10%1.691.751.821.891.972.06
15%1.711.791.892.012.152.31
Taizhou5%14.1913.3612.5111.6310.719.76
10%14.0712.9911.7410.338.746.96
15%13.9512.5910.918.896.473.62
Table 11. Forecasting results at different growth rates of the secondary industry in Shanghai (units: billion cubic meters).
Table 11. Forecasting results at different growth rates of the secondary industry in Shanghai (units: billion cubic meters).
CityGrowth Rate202020212022202320242025
Shanghai−5%99.67100.19101.46103.15105.06107.08
5%97.7794.7991.6888.4385.0381.47
10%96.6791.7786.0379.4972.1564.00
Table 12. Forecasting results at different growth rates of the secondary industry in Anhui (units: billion cubic meters).
Table 12. Forecasting results at different growth rates of the secondary industry in Anhui (units: billion cubic meters).
CityGrowth Rate2020202120222023
Hefei−5%30.6230.7530.8730.99
5%30.3230.1830.0329.88
10%30.1729.8829.5529.18
Tongling−5%12.0511.5010.9810.48
5%12.7412.9713.2513.54
10%13.0813.7614.5515.43
Chizhou−5%10.029.909.779.66
5%10.1110.0810.0610.03
10%10.1610.1810.2210.26
Wuhu−5%29.6829.4529.2128.98
5%30.0430.0930.1530.21
10%30.2330.4430.6930.98
Anqing−5%24.8324.5524.3624.10
5%24.5124.1223.6823.27
10%24.3523.8923.2922.75
Xuancheng−5%14.2214.1414.0714.01
5%14.1814.0813.9913.89
10%14.1614.0413.9313.82
Maanshan−5%26.2223.5020.6017.50
5%25.8622.4818.6814.51
10%25.6721.9417.6112.75
Chuzhou−5%22.3722.4022.4322.46
5%22.2922.2522.2122.17
10%22.2422.1722.0821.99
Table 13. Forecasting results at different growth rates of the tertiary industry in Zhejiang (units: billion cubic meters).
Table 13. Forecasting results at different growth rates of the tertiary industry in Zhejiang (units: billion cubic meters).
CityGrowth Rate202020212022202320242025
Hangzhou5%26.1323.6821.3719.0816.7614.36
10%25.2421.2716.9712.287.141.51
15%24.3418.7712.234.64−4.12−14.20
Ningbo10%20.1819.9319.6719.3919.0818.75
15%20.1319.7919.4018.9518.4417.85
20%20.0819.6419.1118.4717.7016.77
Huzhou10%11.9810.959.848.627.295.83
15%11.8410.538.997.195.082.65
20%11.7010.108.095.602.55−1.15
Jiaxing10%17.6717.3316.9716.5716.1315.65
15%17.5917.1216.5515.9015.1514.28
20%17.5116.8816.0915.1413.9912.60
Shaoxing5%16.8616.3515.8615.3614.8514.33
10%16.7015.9215.0514.0713.0011.81
15%16.5515.4814.1212.6410.848.76
Jinhua5%15.4115.0314.6714.3013.9313.55
10%15.2914.7114.0713.3512.5711.70
15%15.1914.3913.4412.3111.009.48
Zhoushan5%1.651.671.701.721.751.77
10%1.671.711.761.821.881.94
15%1.681.751.831.932.032.16
Taizhou5%13.8213.0012.1411.2610.369.43
10%13.7212.6611.4610.108.576.88
15%13.6012.2910.698.776.483.79
Table 14. Forecasting results at different growth rates of the tertiary industry in Shanghai (units: billion cubic meters).
Table 14. Forecasting results at different growth rates of the tertiary industry in Shanghai (units: billion cubic meters).
CityGrowth Rate202020212022202320242025
Shanghai5%94.7792.4790.3588.2286.0483.76
10%93.7689.7685.5380.9575.9570.46
15%92.5286.7080.0672.4663.7453.73
Table 15. Forecasting results at different growth rates of the tertiary industry in Anhui (units: billion cubic meters).
Table 15. Forecasting results at different growth rates of the tertiary industry in Anhui (units: billion cubic meters).
CityGrowth Rate2020202120222023
Hefei10%30.7330.7030.6830.67
15%30.8530.9531.0931.29
20%30.9731.2131.5431.99
Tongling10%12.6012.7312.9313.18
15%12.8213.2813.8914.65
20%13.0313.8514.9316.31
Chizhou10%10.1110.1110.1210.13
15%10.1510.2010.2610.35
20%10.1910.2910.4310.61
Wuhu10%29.9729.9729.9629.95
15%29.9729.9629.9629.95
20%29.9729.9629.9529.93
Anqing10%24.3523.5023.1922.49
15%24.1423.1422.5821.63
20%23.9222.7721.9220.64
Xuancheng10%14.1313.9713.8513.70
15%14.0513.8313.6213.36
20%13.9713.6813.3612.97
Maanshan10%26.1823.3520.2616.94
15%26.0022.7819.1014.96
20%25.8122.1917.8612.77
Chuzhou10%22.3822.3622.3622.36
15%22.4422.4822.5522.65
20%22.5022.6122.7722.99
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Qiao, Z.; Yang, W. Forecasting Water Consumption in the Yangtze River Delta Based on Deformable Cumulative Multivariable Grey Model. Axioms 2023, 12, 655. https://doi.org/10.3390/axioms12070655

AMA Style

Qiao Z, Yang W. Forecasting Water Consumption in the Yangtze River Delta Based on Deformable Cumulative Multivariable Grey Model. Axioms. 2023; 12(7):655. https://doi.org/10.3390/axioms12070655

Chicago/Turabian Style

Qiao, Zhengran, and Wei Yang. 2023. "Forecasting Water Consumption in the Yangtze River Delta Based on Deformable Cumulative Multivariable Grey Model" Axioms 12, no. 7: 655. https://doi.org/10.3390/axioms12070655

APA Style

Qiao, Z., & Yang, W. (2023). Forecasting Water Consumption in the Yangtze River Delta Based on Deformable Cumulative Multivariable Grey Model. Axioms, 12(7), 655. https://doi.org/10.3390/axioms12070655

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