# A Grey Seasonal Index Model for Forecasting Groundwater Depth of Ningxia Plain

Authors: Ziqi Yin; Kai Zhang
Journal: Discrete Dynamics in Nature and Society
 (2021)

---

## Abstract

Forecasting the depth of groundwater in arid and semiarid areas is a great challenge because these areas are complex hydrogeological environments and the observational data are limited. To deal with this problem, the grey seasonal index model is proposed. The seasonal characteristics of time series were represented by indicators, and the grey model with fractional-order accumulation was employed to fit and forecast different periodic indicators and long-term trends, respectively. Then, the prediction results of the two were combined together to obtain the prediction results. To verify the model performance, the proposed model is applied to groundwater prediction in Yinchuan Plain. The results show that the fitting error of the proposed model is 2.08%, while for comparison, the fitting error of the grey model of data grouping and Holt–Winters model is 3.94% and 5%, respectively. In the same way, it is concluded that the fitting error of groundwater in Weining Plain by the proposed model is 2.26%. On the whole, the groundwater depth in Ningxia Plain including Yinchuan Plain and Weining Plain will increase further.


---

## Body

## 1. Introduction

### 1.1. Motivation

Groundwater is closely related to the development of human civilization. More than 2 billion people worldwide drink groundwater directly. This number is rising as rivers are polluted and depleted. The overuse of groundwater destroys the natural water cycle and causes problems such as ground collapse. China’s groundwater resources have been drastically reduced due to excessive exploitation of groundwater. At present, the accumulated deficit of groundwater deficit in North China is 180 billion cubic meters. In order to reduce the impact of groundwater depth decline, it is meaningful to achieve accurate groundwater prediction. Especially in dry areas, the reliance on groundwater is greater.Ningxia Hui Autonomous Region (Ningxia) is located in Northwest China and is a typical dry area. Ningxia has the fewest water resources in China. The average annual runoff in Ningxia is 18.3 mm deep, only 1/15 of the national average, with annual rainfall from 150 mm to 600 mm. To solve the problem of farmland irrigation, the ancient people of Ningxia built water conservancy projects to divert the Yellow River. The water conservancy project is still in use today. In 2017, the ancient irrigation area in Ningxia was listed on the world irrigation engineering heritage list. Although irrigation projects have eased the drought on farmland, groundwater is still the main source of water resources in Ningxia. Moreover, with the continuous development of China’s western development strategy, Ningxia’s economy has developed rapidly, which has intensified the demand for water resources. According to Ningxia water resource bulletin, the groundwater consumption in Ningxia has been on the rise from 2014 to 2018. The content comes from Ningxia Hydrological Information Network.In order to provide reference for groundwater management in Ningxia, the grey seasonal index model (GSIM (1,1)) was established to predict the groundwater depth of Ningxia Plain and compared with the results of Holt–Winters model (Holt–Winters) and grey model of data grouping (DGGM (1,1)).This paper is organized as follows. The model is introduced in Section2. Experimental results are given in Section 3. In Section 4, some conclusions and directions for future research are provided.


### 1.2. Literature Review

Under the background of water resource depletion, groundwater management has gradually developed into a discipline to be studied by scholars [1]. Accurate prediction of groundwater is important for groundwater management. There are many studies on prediction method [2]. The most direct way is through the surface vegetation to predict groundwater depth [3]. However, changes in groundwater are most affected by humans. Industrial water, agricultural water, and domestic water are considered as the main factors affecting groundwater [4]. Based on different influencing factors, stochastic time series and artificial neural network models are considered to have better effect in groundwater depth assessment [5]. Rivers are the main natural way to replenish groundwater and play a regulating role in the water cycle. A groundwater prediction model based river stage was built [6]. Based on groundwater dynamic data and related factors, the quantile regression method was used to predict groundwater depth [7]. Other factors affecting groundwater depth have also been investigated [8]. The development of science and technology provides convenience for the study of groundwater. The conceptual model, the Bayesian network, and artificial network models were applied to groundwater depth prediction [9]. Support vector machines have good performance in groundwater depth prediction [10]. The prediction effects of grey self-memory model, radial basis function network, and adaptive neuro fuzzy inference system models on groundwater were compared [11]. The mathematical models have their own characteristics, and different models are combined and applied to groundwater depth prediction [12]. Hybrid support vector machine regression and artificial neural network models were applied to groundwater depth prediction [13]. A hybrid model using entropy spectral analysis with the optimal input is proved to be effective in predicting groundwater depth [14]. The empirical mode decomposition (EMD) was combined with Elman neural network, and the coupling forecasting model was constructed to groundwater depth forecasting [15]. The EMD is also used in combination with phase space reconstruction, particle swarm optimization, and extreme learning machine [16].Groundwater is a huge system. Due to the influence of geological conditions, hydrological environment, and other factors [17–20], the information about groundwater is ambiguous. Therefore, groundwater is a complex grey system. Moreover, the above model needs a lot of data to guarantee the accuracy of prediction [21]. Without a large number of training samples, the prediction accuracy of the model cannot reach the best [22]. The grey prediction model has high precision in small sample prediction [23, 24]. Grey models are widely used in air pollution, energy, electricity, biology, computer, water resources, and other fields [25–33]. The fractional-order grey prediction model improves the stability of the model from the perspective of fractional order [34, 35]. On this basis, the seasonal index is introduced into the fractional-order grey prediction model, and the processing ability of the model to nonlinear data is improved. Thus, PSO-GSIM (1,1) is used to predict the groundwater depth of Ningxia Plain.


## 1.1. Motivation

Groundwater is closely related to the development of human civilization. More than 2 billion people worldwide drink groundwater directly. This number is rising as rivers are polluted and depleted. The overuse of groundwater destroys the natural water cycle and causes problems such as ground collapse. China’s groundwater resources have been drastically reduced due to excessive exploitation of groundwater. At present, the accumulated deficit of groundwater deficit in North China is 180 billion cubic meters. In order to reduce the impact of groundwater depth decline, it is meaningful to achieve accurate groundwater prediction. Especially in dry areas, the reliance on groundwater is greater.Ningxia Hui Autonomous Region (Ningxia) is located in Northwest China and is a typical dry area. Ningxia has the fewest water resources in China. The average annual runoff in Ningxia is 18.3 mm deep, only 1/15 of the national average, with annual rainfall from 150 mm to 600 mm. To solve the problem of farmland irrigation, the ancient people of Ningxia built water conservancy projects to divert the Yellow River. The water conservancy project is still in use today. In 2017, the ancient irrigation area in Ningxia was listed on the world irrigation engineering heritage list. Although irrigation projects have eased the drought on farmland, groundwater is still the main source of water resources in Ningxia. Moreover, with the continuous development of China’s western development strategy, Ningxia’s economy has developed rapidly, which has intensified the demand for water resources. According to Ningxia water resource bulletin, the groundwater consumption in Ningxia has been on the rise from 2014 to 2018. The content comes from Ningxia Hydrological Information Network.In order to provide reference for groundwater management in Ningxia, the grey seasonal index model (GSIM (1,1)) was established to predict the groundwater depth of Ningxia Plain and compared with the results of Holt–Winters model (Holt–Winters) and grey model of data grouping (DGGM (1,1)).This paper is organized as follows. The model is introduced in Section2. Experimental results are given in Section 3. In Section 4, some conclusions and directions for future research are provided.


## 1.2. Literature Review

Under the background of water resource depletion, groundwater management has gradually developed into a discipline to be studied by scholars [1]. Accurate prediction of groundwater is important for groundwater management. There are many studies on prediction method [2]. The most direct way is through the surface vegetation to predict groundwater depth [3]. However, changes in groundwater are most affected by humans. Industrial water, agricultural water, and domestic water are considered as the main factors affecting groundwater [4]. Based on different influencing factors, stochastic time series and artificial neural network models are considered to have better effect in groundwater depth assessment [5]. Rivers are the main natural way to replenish groundwater and play a regulating role in the water cycle. A groundwater prediction model based river stage was built [6]. Based on groundwater dynamic data and related factors, the quantile regression method was used to predict groundwater depth [7]. Other factors affecting groundwater depth have also been investigated [8]. The development of science and technology provides convenience for the study of groundwater. The conceptual model, the Bayesian network, and artificial network models were applied to groundwater depth prediction [9]. Support vector machines have good performance in groundwater depth prediction [10]. The prediction effects of grey self-memory model, radial basis function network, and adaptive neuro fuzzy inference system models on groundwater were compared [11]. The mathematical models have their own characteristics, and different models are combined and applied to groundwater depth prediction [12]. Hybrid support vector machine regression and artificial neural network models were applied to groundwater depth prediction [13]. A hybrid model using entropy spectral analysis with the optimal input is proved to be effective in predicting groundwater depth [14]. The empirical mode decomposition (EMD) was combined with Elman neural network, and the coupling forecasting model was constructed to groundwater depth forecasting [15]. The EMD is also used in combination with phase space reconstruction, particle swarm optimization, and extreme learning machine [16].Groundwater is a huge system. Due to the influence of geological conditions, hydrological environment, and other factors [17–20], the information about groundwater is ambiguous. Therefore, groundwater is a complex grey system. Moreover, the above model needs a lot of data to guarantee the accuracy of prediction [21]. Without a large number of training samples, the prediction accuracy of the model cannot reach the best [22]. The grey prediction model has high precision in small sample prediction [23, 24]. Grey models are widely used in air pollution, energy, electricity, biology, computer, water resources, and other fields [25–33]. The fractional-order grey prediction model improves the stability of the model from the perspective of fractional order [34, 35]. On this basis, the seasonal index is introduced into the fractional-order grey prediction model, and the processing ability of the model to nonlinear data is improved. Thus, PSO-GSIM (1,1) is used to predict the groundwater depth of Ningxia Plain.


## 2. Methodology

### 2.1. The Holt–Winters Model

The Holt–Winters model is a traditional statistical model and has been widely used since it was proposed by Winters in 1960. Scholars have conducted in-depth research on the model [36]. The Holt–Winters model is used to process data with seasonal characteristics. Time series data are divided into three parts, namely, residual data, trend data, and seasonal data. A prediction model is established based on the Holt–Winters model. There are two kinds of calculation methods for the Holt–Winters model, namely, addition and multiplication. The multiplication Holt–Winters model works better with solving the seasonal time series prediction problem. Therefore, the multiplication Holt–Winters model is used for comparison in this paper. The establishment process of the model is as follows.For the time seriesX0=x01,x02,…,x0n, the multiplication Holt–Winters model is(1)St=αx0tIt−L+1−αSt−1+bt−1,0<α<1,bt=γSt−St−1+1−γbt−1,0<γ<1,It=βx0tSt1−βIt−L,0<β<1,where L is the seasonal length, such as four quarters, 12 months, and so on, It is the seasonal index, and bt represents the tendency.(2)It=x0L¯x0t¯,SL+1=x0L+1,bL+1=x0L+1−x01+x0L+2−x02+x0L+3−x033L,where x0L¯ represents the average of same month in different years and x0t¯ is general average.The Holt–Winters model has three parameters, and the application of the grey wolf optimization algorithm (GWO) improves the efficiency of multiparameter optimization. GWO is an optimization algorithm inspired by the hunting behavior of wolves. It is widely used for multiparameter model optimization. The behavior of a pack of wolves is represented digitally as follows:(3)D=C⋅xpi−xi,xi+1=xpi−A⋅D,where i represents the literation at this point, A and D are coefficient vectors, and xp represents the position vector of the prey, i.e. objective vector. xi+1 is the position of one of the wolves in the pack. Coefficient vectors A and D are obtained by the following formula:(4)A=2a⋅r1−a,C=2⋅r2,where a decreases linearly from 2 to 0 along with iteration and r1 and r2 are randomly generated from 0,1.In GWO, there are three optimization parameters,α, β, and δ, corresponding to the three variables in the Holt–Winters model, respectively, where alpha is the dominant player in GWO. Suppose that the three wolves knew where their prey was; therefore, the first three best hunting sites were preserved. Other searcher is forced to update their locations on this basis. The other search location update formula is as follows:(5)Dα=C1⋅xαi−xi,Dβ=C2⋅xβi−xi,Dδ=C3⋅xδi−xi,x1i+1=xαi−A1⋅Dα,x2i+1=xβi−A2⋅Dβ,x3i+1=xδi−A3⋅Dδ,xi+1=x1i+1+x2i+1+x3i+13,where xi is the position of other searcher, xα, xβ, and xδ are the search locations of alpha, beta, and delta, and xi+1 is where the other searchers are in the next iteration. To speed up search efficiency, the maximum number of iterations is set to 30. The value range of parameters is 0,1. Detailed description and source codes of GWO can be found in the literature [37].


### 2.2. Grey Model of Data Grouping Method

The traditional grey prediction model is not applicable to the data of seasonal fluctuation. For the grey model to process the seasonal data, DGGM (1,1) model was proposed [38]. Two more steps are added to the traditional GM (1,1) model. The calculation process of DGGM (1,1) model is described below.Step 1. The time series data are divided into twelve monthly groups and are given by(6)x0τ=x0τ,1,x0τ,2,…,x0τ,n,τ=1,2,…,12.Step 2. GM (1,1) models are calculated based on grouped data separately. An accumulative series x1τ is obtained by accumulative calculation:(7)x1τ=x1τ,1,x1τ,2,…,x1τ,n=x1τ,1,x1τ,1+x0τ,2,…,x1τ,n−1+x0τ,n−1,τ=1,2,…,12,wherex1τ,k=∑i=1kx0τ,i,k=1,2,…,n;τ=1,2,…,12, and the mean series is z1τ,k=0.5x1τ,k+0.5x1τ,k−1k=2,3,…,n;τ=1,2,…,12.The first-order differential equation of a single variable is used as the prediction model, namely, GM (1,1), and the grey differential equation is as follows:(8)x0τ,k+az1τ,k=b,k=2,3,…,n;τ=1,2,…,12.The corresponding whitening differential equation isdx1τ,t/dt+ax1τ,t=b, where t,a, and b represent the time variable, grey development coefficient, and grey action quantity, respectively. Data matrix B and data vector Y are given by(9)B=−12x1τ,1+x1τ,21−12x1τ,2+x1τ,31⋮⋮−12x1τ,n−1+x1τ,n1,Y=x0τ,2x0τ,3⋮x0τ,n,τ=1,2,…,12.The undetermined coefficient vector isa∧b∧=BTB−1BTY, and the prediction equation is obtained:(10)x∧1τ,t+1=x0τ,1−bae−at+ba,τ=1,2,…,12.Step 3. Reductive treatment is conducted. Through a subtraction calculation, the predicted values of original series x0 are obtained:(11)x∧0τ,t+1=x∧1τ,t+1−x∧1τ,t,t=1,2,…,n−1;τ=1,2,…,12.Furthermore, the twelve groups of monthly predictions are combined into a new time series:(12)x∧0τ,t+1=x∧01,1,x∧02,1,x∧03,1,x∧04,1,…,x∧01,n,x∧02,n,x∧03,n,x∧04,n.


### 2.3. The GSIM (1,1) Model

The data varied substantially among different months, but the traditional grey prediction model is only suitable for time series with an exponential trend, and it is not capable of effectively predicting data with large fluctuations. Thus, in order to enhance the forecasting accuracy, the seasonal data are summarized into annual data in this study. The annual data are predicted by GM (1,1) with fractional-order accumulation (FGM (1,1)) [24]. The prediction results were calculated by seasonal index and annual data.The time series data on the depth of groundwater areX0=x01,x02,…,x0n with a cycle of L, for example, L=12, and the modeling process for the GSIM (1,1) model is as follows.Step 1. The yearly data are Y0=y01,y02,…,y0ϕ,ϕ=n/L, where y01=∑i=1Lx0i,y02=∑i=L+12Lx0i,…,y0ϕ=∑i=ϕ−1L+1nx0i. φ is the number of periods.Step 2. By using yrk=∑j=1kCk−j+r−1k−jy0j, the r-order accumulation sequence is(13)Yr=yr1,yr2,…,yrφ,whereCr−10=1,Ckk+1=0,Ck−j+r−1k−j=k−j+r−1k−j+r−2…r+1r/k−j!. r is the fractional order, and the optimal value can be obtained through particle swarm optimization (PSO) algorithm.Step 3. The differential equation with one variable of the r-order accumulation sequence Yr (i.e., the FGM (1,1) model) can be defined as(14)dyrdt+ayr=b,wherea and b represent grey development coefficient and grey action quantity, respectively. The solution to equation (14) is(15)yrt+1=y01−bae−at+ba.The parameters of the model are calculated by the least square method. The calculation process ofa∧, b∧ is shown in the following equation:(16)a∧b∧=BTB−1BTY,where(17)B=−0.5yr1+yr21−0.5yr2+yr31⋮⋮−0.5yrφ−1+yrφ1,Y=yr2yr3⋮yrφ.Step 4. a∧ and b∧ are entered in the time response function:(18)y∧rk+1=y01−b∧a∧e−a∧k+b∧a∧,andy∧rk+1 is the fitting value at time k+1.Step 5. For Y∧r=y∧r1,y∧r2,…,y∧rφ, the restored sequence is(19)αrY∧r=α1y∧r1−r1,α1y∧r1−r2,…,α1y∧r1−rφ,whereα1y∧r1−rk=y∧r1−rk−y∧r1−rk−1. Then, the restored value is(20)y∧01,y∧02,…,y∧0φ,y∧0φ+1,y∧0φ+2,….Step 6. The seasonal index di,ti=1,2,…,L,t=1,2,…,ϕ is calculated by using the average method:(21)d0i,t=d01,1,d02,1,…,d0L,1,d01,2,…,d0L,φ,where(22)d01,1=x01,1∑i=1Lx0i,1,d02,1=x02,1∑i=1Lx0i,1,…,d0L,1=x0L,1∑i=1Lx0i,1,d01,2=x01,2∑i=1Lx0i,2,⋯,d0L,φ=x0L,φ∑i=1Lx0i,φ.Step 7. According to the calculation process in Step 2, the accumulation series dri,t is obtained.(23)dri,t=dr1,1,dr1,2,…dr1,φ;dr2,1,dr2,2,…dr2,φ;⋮drL,1,drL,2,…drL,φ.The accumulation of the seasonal index is no longer carried out in accordance with the time series, but data of the same month are combined into a group for the accumulation. The purpose of this is to reduce the dimension of the data so as to make full use of the grey model to process the small data. At the same time, it is also preparation for the prediction of the next cycle.The next calculation refers to Steps 3–5, and the forecasting value is obtained:(24)d∧01,φ+1,d∧02,φ+1,…,d∧0L,φ+1,d∧01,φ+2,d∧02,φ+2,….Step 8. The fitting value of the groundwater depth is obtained.ForX∧0=x∧01,x∧02,…,x∧0n,x∧0n+1,…, the solution formula is(25)x∧0k=y∧0ω+1d∧0υ,ω+1L,modk,L≠0,x∧0k=y∧0ωd∧0L,ωL,modk,L=0,whereω and υ are the quotient and remainder of k/L, respectively.Step 9. The mean absolute percentage error (MAPE) is used to test model accuracy, and the calculation process is as follows:(26)MAPE=1n∑k=1nx∧0k−x0kx0k×100%.In the paper, PSO is applied to obtain the optimal value of the parameterr of the GSIM (1,1) model. PSO has been widely used since it was proposed [39]. Compared with other optimization algorithms, PSO is simpler and easier to implement without too much parameter adjustment. To increase the search effect, the position and velocity of the particle need to be constantly updated. The updating process depends on the values Pbest and gbest, where Pbest represents the current optimal solution and gbest represents the current optimal solution for any particle. The updating process is shown in the following equations:(27)vi+1=ωvi+c1∗rand∗Pbesti−xi+c2∗rand∗gbesti−xi,(28)xi+1=xi+vi+1,where vi and xi represent the velocity and position of the particle at i, respectively, ω represents the inertia weight value which is generally set as 0.8 [40], and rand is a random function; the optimal value of r of most experiments is generally within the range of [0, 1]; to ensure the accuracy of the experiment, the value range of GSIM (1,1) parameter r is expanded to [0, 2.5]; therefore, rand represents the random number generated from [0, 2.5]. c1 and c2 are acceleration constants with equal values and set to 2 [41]. Suppose the number of individuals in the initial population is N=50, the maximum number of iterations is 200. More detailed description of PSO can be found in literature [42].


## 2.1. The Holt–Winters Model

The Holt–Winters model is a traditional statistical model and has been widely used since it was proposed by Winters in 1960. Scholars have conducted in-depth research on the model [36]. The Holt–Winters model is used to process data with seasonal characteristics. Time series data are divided into three parts, namely, residual data, trend data, and seasonal data. A prediction model is established based on the Holt–Winters model. There are two kinds of calculation methods for the Holt–Winters model, namely, addition and multiplication. The multiplication Holt–Winters model works better with solving the seasonal time series prediction problem. Therefore, the multiplication Holt–Winters model is used for comparison in this paper. The establishment process of the model is as follows.For the time seriesX0=x01,x02,…,x0n, the multiplication Holt–Winters model is(1)St=αx0tIt−L+1−αSt−1+bt−1,0<α<1,bt=γSt−St−1+1−γbt−1,0<γ<1,It=βx0tSt1−βIt−L,0<β<1,where L is the seasonal length, such as four quarters, 12 months, and so on, It is the seasonal index, and bt represents the tendency.(2)It=x0L¯x0t¯,SL+1=x0L+1,bL+1=x0L+1−x01+x0L+2−x02+x0L+3−x033L,where x0L¯ represents the average of same month in different years and x0t¯ is general average.The Holt–Winters model has three parameters, and the application of the grey wolf optimization algorithm (GWO) improves the efficiency of multiparameter optimization. GWO is an optimization algorithm inspired by the hunting behavior of wolves. It is widely used for multiparameter model optimization. The behavior of a pack of wolves is represented digitally as follows:(3)D=C⋅xpi−xi,xi+1=xpi−A⋅D,where i represents the literation at this point, A and D are coefficient vectors, and xp represents the position vector of the prey, i.e. objective vector. xi+1 is the position of one of the wolves in the pack. Coefficient vectors A and D are obtained by the following formula:(4)A=2a⋅r1−a,C=2⋅r2,where a decreases linearly from 2 to 0 along with iteration and r1 and r2 are randomly generated from 0,1.In GWO, there are three optimization parameters,α, β, and δ, corresponding to the three variables in the Holt–Winters model, respectively, where alpha is the dominant player in GWO. Suppose that the three wolves knew where their prey was; therefore, the first three best hunting sites were preserved. Other searcher is forced to update their locations on this basis. The other search location update formula is as follows:(5)Dα=C1⋅xαi−xi,Dβ=C2⋅xβi−xi,Dδ=C3⋅xδi−xi,x1i+1=xαi−A1⋅Dα,x2i+1=xβi−A2⋅Dβ,x3i+1=xδi−A3⋅Dδ,xi+1=x1i+1+x2i+1+x3i+13,where xi is the position of other searcher, xα, xβ, and xδ are the search locations of alpha, beta, and delta, and xi+1 is where the other searchers are in the next iteration. To speed up search efficiency, the maximum number of iterations is set to 30. The value range of parameters is 0,1. Detailed description and source codes of GWO can be found in the literature [37].


## 2.2. Grey Model of Data Grouping Method

The traditional grey prediction model is not applicable to the data of seasonal fluctuation. For the grey model to process the seasonal data, DGGM (1,1) model was proposed [38]. Two more steps are added to the traditional GM (1,1) model. The calculation process of DGGM (1,1) model is described below.Step 1. The time series data are divided into twelve monthly groups and are given by(6)x0τ=x0τ,1,x0τ,2,…,x0τ,n,τ=1,2,…,12.Step 2. GM (1,1) models are calculated based on grouped data separately. An accumulative series x1τ is obtained by accumulative calculation:(7)x1τ=x1τ,1,x1τ,2,…,x1τ,n=x1τ,1,x1τ,1+x0τ,2,…,x1τ,n−1+x0τ,n−1,τ=1,2,…,12,wherex1τ,k=∑i=1kx0τ,i,k=1,2,…,n;τ=1,2,…,12, and the mean series is z1τ,k=0.5x1τ,k+0.5x1τ,k−1k=2,3,…,n;τ=1,2,…,12.The first-order differential equation of a single variable is used as the prediction model, namely, GM (1,1), and the grey differential equation is as follows:(8)x0τ,k+az1τ,k=b,k=2,3,…,n;τ=1,2,…,12.The corresponding whitening differential equation isdx1τ,t/dt+ax1τ,t=b, where t,a, and b represent the time variable, grey development coefficient, and grey action quantity, respectively. Data matrix B and data vector Y are given by(9)B=−12x1τ,1+x1τ,21−12x1τ,2+x1τ,31⋮⋮−12x1τ,n−1+x1τ,n1,Y=x0τ,2x0τ,3⋮x0τ,n,τ=1,2,…,12.The undetermined coefficient vector isa∧b∧=BTB−1BTY, and the prediction equation is obtained:(10)x∧1τ,t+1=x0τ,1−bae−at+ba,τ=1,2,…,12.Step 3. Reductive treatment is conducted. Through a subtraction calculation, the predicted values of original series x0 are obtained:(11)x∧0τ,t+1=x∧1τ,t+1−x∧1τ,t,t=1,2,…,n−1;τ=1,2,…,12.Furthermore, the twelve groups of monthly predictions are combined into a new time series:(12)x∧0τ,t+1=x∧01,1,x∧02,1,x∧03,1,x∧04,1,…,x∧01,n,x∧02,n,x∧03,n,x∧04,n.


## 2.3. The GSIM (1,1) Model

The data varied substantially among different months, but the traditional grey prediction model is only suitable for time series with an exponential trend, and it is not capable of effectively predicting data with large fluctuations. Thus, in order to enhance the forecasting accuracy, the seasonal data are summarized into annual data in this study. The annual data are predicted by GM (1,1) with fractional-order accumulation (FGM (1,1)) [24]. The prediction results were calculated by seasonal index and annual data.The time series data on the depth of groundwater areX0=x01,x02,…,x0n with a cycle of L, for example, L=12, and the modeling process for the GSIM (1,1) model is as follows.Step 1. The yearly data are Y0=y01,y02,…,y0ϕ,ϕ=n/L, where y01=∑i=1Lx0i,y02=∑i=L+12Lx0i,…,y0ϕ=∑i=ϕ−1L+1nx0i. φ is the number of periods.Step 2. By using yrk=∑j=1kCk−j+r−1k−jy0j, the r-order accumulation sequence is(13)Yr=yr1,yr2,…,yrφ,whereCr−10=1,Ckk+1=0,Ck−j+r−1k−j=k−j+r−1k−j+r−2…r+1r/k−j!. r is the fractional order, and the optimal value can be obtained through particle swarm optimization (PSO) algorithm.Step 3. The differential equation with one variable of the r-order accumulation sequence Yr (i.e., the FGM (1,1) model) can be defined as(14)dyrdt+ayr=b,wherea and b represent grey development coefficient and grey action quantity, respectively. The solution to equation (14) is(15)yrt+1=y01−bae−at+ba.The parameters of the model are calculated by the least square method. The calculation process ofa∧, b∧ is shown in the following equation:(16)a∧b∧=BTB−1BTY,where(17)B=−0.5yr1+yr21−0.5yr2+yr31⋮⋮−0.5yrφ−1+yrφ1,Y=yr2yr3⋮yrφ.Step 4. a∧ and b∧ are entered in the time response function:(18)y∧rk+1=y01−b∧a∧e−a∧k+b∧a∧,andy∧rk+1 is the fitting value at time k+1.Step 5. For Y∧r=y∧r1,y∧r2,…,y∧rφ, the restored sequence is(19)αrY∧r=α1y∧r1−r1,α1y∧r1−r2,…,α1y∧r1−rφ,whereα1y∧r1−rk=y∧r1−rk−y∧r1−rk−1. Then, the restored value is(20)y∧01,y∧02,…,y∧0φ,y∧0φ+1,y∧0φ+2,….Step 6. The seasonal index di,ti=1,2,…,L,t=1,2,…,ϕ is calculated by using the average method:(21)d0i,t=d01,1,d02,1,…,d0L,1,d01,2,…,d0L,φ,where(22)d01,1=x01,1∑i=1Lx0i,1,d02,1=x02,1∑i=1Lx0i,1,…,d0L,1=x0L,1∑i=1Lx0i,1,d01,2=x01,2∑i=1Lx0i,2,⋯,d0L,φ=x0L,φ∑i=1Lx0i,φ.Step 7. According to the calculation process in Step 2, the accumulation series dri,t is obtained.(23)dri,t=dr1,1,dr1,2,…dr1,φ;dr2,1,dr2,2,…dr2,φ;⋮drL,1,drL,2,…drL,φ.The accumulation of the seasonal index is no longer carried out in accordance with the time series, but data of the same month are combined into a group for the accumulation. The purpose of this is to reduce the dimension of the data so as to make full use of the grey model to process the small data. At the same time, it is also preparation for the prediction of the next cycle.The next calculation refers to Steps 3–5, and the forecasting value is obtained:(24)d∧01,φ+1,d∧02,φ+1,…,d∧0L,φ+1,d∧01,φ+2,d∧02,φ+2,….Step 8. The fitting value of the groundwater depth is obtained.ForX∧0=x∧01,x∧02,…,x∧0n,x∧0n+1,…, the solution formula is(25)x∧0k=y∧0ω+1d∧0υ,ω+1L,modk,L≠0,x∧0k=y∧0ωd∧0L,ωL,modk,L=0,whereω and υ are the quotient and remainder of k/L, respectively.Step 9. The mean absolute percentage error (MAPE) is used to test model accuracy, and the calculation process is as follows:(26)MAPE=1n∑k=1nx∧0k−x0kx0k×100%.In the paper, PSO is applied to obtain the optimal value of the parameterr of the GSIM (1,1) model. PSO has been widely used since it was proposed [39]. Compared with other optimization algorithms, PSO is simpler and easier to implement without too much parameter adjustment. To increase the search effect, the position and velocity of the particle need to be constantly updated. The updating process depends on the values Pbest and gbest, where Pbest represents the current optimal solution and gbest represents the current optimal solution for any particle. The updating process is shown in the following equations:(27)vi+1=ωvi+c1∗rand∗Pbesti−xi+c2∗rand∗gbesti−xi,(28)xi+1=xi+vi+1,where vi and xi represent the velocity and position of the particle at i, respectively, ω represents the inertia weight value which is generally set as 0.8 [40], and rand is a random function; the optimal value of r of most experiments is generally within the range of [0, 1]; to ensure the accuracy of the experiment, the value range of GSIM (1,1) parameter r is expanded to [0, 2.5]; therefore, rand represents the random number generated from [0, 2.5]. c1 and c2 are acceleration constants with equal values and set to 2 [41]. Suppose the number of individuals in the initial population is N=50, the maximum number of iterations is 200. More detailed description of PSO can be found in literature [42].


## 3. Experimental Results

### 3.1. Study Area and Data

Ningxia spans 35°14′–39°14′N and 104°17′–109°39′E and includes five cities. The climatic types of Ningxia are temperate, continental, arid, and semiarid. Ningxia Plain is considered as the largest plain in Ningxia. The study region has an area of 66.400 km2 and is bordered by the Helan Mountains and Liupan Mountains to the northwestward and south, respectively. Ningxia is 1100–1200 meters above sea level and the terrain slopes from southwest to northeast. The Yellow River enters from the southwest of Ningxia and leaves from the northeast. Due to the different physiognomy types in northcentral Ningxia and south Ningxia, Ningxia Plain can be divided into Yinchuan Plain and Weining Plain (Figure 1). The recoverable groundwater accounts for 88% of the total groundwater in Ningxia. Ningxia Plain is surrounded by mountains, water flows to the plain, and the water table changes slowly and steadily, creating a predictable gradient of water [43].Figure 1
The general area of Ningxia Plain.The data are from the official website of the Ministry of Water Resources of the People’s Republic of China (https://www.mwr.gov.cn/sj/tjgb/dxsdtyb/). There are a total of 47 groundwater monitoring wells in Ningxia Plain. Yinchuan Plain contains three cities and 28 groundwater monitoring wells. Weining Plain contains one city and 11 groundwater monitoring wells (by the end of the 2019). The groundwater depth data for Yinchuan Plain and Weining Plain are the average result of multiple groundwater monitoring wells.


### 3.2. Forecasting the Groundwater Depth of Yinchuan Plain

The groundwater depth of Yinchuan Plain is predicted in this section. The depth of groundwater is one of the contents of hydrology. Groundwater is subject to changes in rainfall or in rivers. As we all know, changes in rainfall and rivers are cyclical and seasonal. In recent years, the seasonal variation of groundwater depth is more obvious under the influence of human production activities. For example, the original groundwater depth data of Yinchuan Plain are summarized in Figure2. Figure 2 shows that the groundwater depth changes with time, and groundwater depth fluctuations are similar at the same time in different years. Therefore, the seasonal model is used to study the variation of groundwater depth.Figure 2
The groundwater depth of Yinchuan Plain and Weining Plain.


#### 3.2.1. The Calculation Results of the Holt–Winters Model

The groundwater depth data of Yinchuan Plain from 2014 to 2018 were fitted to three models. The models do not fit the data of the first period, so the data of the first period are omitted. The fitting results of the Holt–Winters model are listed in Table1. The calculation process is completed by MATLAB 2016b. The GWO was used to find optimal parameters α,β, and γ, and the optimization result of parameters is [0.51639, 0.05160, 0.00752]. The fitting error (MAPE) of the Holt–Winters model is 5.00%. It can be seen that the fitting accuracy of the traditional model is lower than 10%, and the fitting results are accepted. As the first cycle, the data of 2014 were omitted.Table 1
The fitting value of Yinchuan Plain groundwater depth by Holt–Winters model (unit: m).
MonthYear201520162017201812.061.692.002.0922.272.192.482.6532.632.532.712.8942.692.622.752.9252.752.682.813.0761.941.902.052.2771.711.701.812.1081.551.591.692.0191.911.812.072.23102.462.182.582.94112.612.573.023.16121.531.561.791.85MAPE5.00


#### 3.2.2. The Calculation Results of the DGGM (1,1) Model

The DGGM (1,1) model was applied to the data calculation on a monthly basis. The data grouping method is adopted to improve the fitting accuracy of the traditional GM (1,1) model, and the model can be applied to process seasonal fluctuation data. The fitting results of DGGM (1,1) model are listed in Table2. MAPE is 3.94%. Compared with the Holt–Winters model, the new model has higher fitting accuracy. Moreover, there are no fitting results for the first period of the two models, increasing contrast between models.Table 2
The fitting value of Yinchuan Plain groundwater depth by DGGM (1,1) model (unit: m).
MonthYear201520162017201812.031.972.112.1822.382.512.452.6232.562.742.792.9142.642.672.843.0952.562.783.073.2461.832.082.142.2871.641.671.892.0781.611.741.671.7891.821.932.092.26102.332.412.622.77112.452.672.913.05121.311.451.631.8MAPE3.94The parameters generated by the DGGM (1,1) model in the fitting process are shown in Table3. The time response function of each set of parameters can be obtained from equation (16).Table 3
The parameters of the DGGM (1,1) model fitting process.
Modeling the serial numberGrey development coefficienta^Grey action quantityb^1−0.0391.8382−0.0192.3813−0.0262.5394−0.0472.3955−0.0722.2796−0.0851.5747−0.1011.3378−0.0431.4959−0.0861.5310−0.0752.00311−0.0692.26312−0.1160.992


#### 3.2.3. The Calculation Results of the GSIM (1,1) Model

The seasonal index in this paper refers to the index calculated on a monthly basis. The seasonal index reflects the fluctuation of the data in the time series. Most models do not deal well with fluctuating data. For the purpose of stabilizing the sequence, scholars used seasonal factors to eliminate seasonal changes in time series. This method improves the ability of the model to process the fluctuation data. However, the time series that eliminate seasonal effects may lose a lot of information. Its accuracy is also difficult to improve. To reduce information loss in the sequence, the seasonal index was extracted and studied separately. The results of yearly data are listed in Table4. The MAPE is 0.38%, and the fitting accuracy is high.Table 4
The annual fitting value of Yinchuan Plain groundwater depth (unit: m).
YearActual valueFitting value201426.6826.68201525.8325.81201625.7125.70201728.0227.70201830.6630.74MAPE0.38Then, the seasonal index of monthly data is fitted. During the fitting process, the GSIM (1,1) model is used sometimes, so the multiple sets of parameters are generated. Each set of parameters in Table5 can be substituted into equation (20) to obtain the corresponding time response function.Table 5
The parameters of the GSIM (1,1) model fitting process in Yinchuan Plain.
Modeling the serial numberThe optimal order numberra^b^11.21−0.051.0820.990.051.2430.150.891.4041.13−0.031.3551.10−0.041.2861.10−0.050.9071.11−0.070.7781.25−0.070.8791.58−0.191.11100.540.180.83110.570.170.95120.98−0.050.55The fitting values of the seasonal indices are listed in Table6. The average MAPE from 2015 to 2018 is 2.01%. Groundwater is consumed more in February, March, April, May, October, and November. Affected by dry and wet seasons, the groundwater depth varies seasonally throughout the year. The GSIM (1,1) model can not only predict accurately but also quantify the periodic variation of data within the time span. With the allowable error, it is beneficial to understand the variation law of groundwater depth.Table 6
Seasonal index fitting value of groundwater depth in Yinchuan Plain.
Month20142015201620172018Actual valueActual valueFitting valueActual valueFitting valueActual valueFitting valueActual valueFitting value10.920.960.960.890.890.890.880.880.8921.121.171.171.131.121.071.081.041.0331.231.251.251.241.221.171.181.141.1441.231.251.251.221.211.181.201.201.1951.221.241.241.231.231.251.251.271.2760.850.880.880.900.880.870.900.950.9370.730.770.770.800.780.780.810.870.8580.710.780.780.700.730.780.720.710.7490.810.920.920.770.770.870.840.940.96101.041.041.071.121.111.191.121.071.10111.201.181.201.251.241.281.241.201.22120.930.570.630.760.660.670.700.720.73MAPE2.01The groundwater depth fitting results from 2014 to 2018 are shown in Table7. The MAPE of fitting results is 2.08%. By calculating the seasonal index and annual fitting value of groundwater depth, the monthly fitting value is obtained. As the trend parameter, the annual fitting value controls the change direction of groundwater. The seasonal index breaks down the units of time from years into months.Table 7
The fitting value of Yinchuan Plain groundwater depth by GSIM (1,1) model (unit: m).
MonthYear2014201520162017201812.052.061.912.032.2722.492.522.412.492.6532.732.682.612.722.9242.732.692.602.763.0552.712.672.642.883.2661.891.891.892.082.3971.621.661.681.882.1881.581.681.551.671.9091.801.981.651.932.45102.312.292.382.582.83112.672.582.652.863.12122.071.351.421.611.88MAPE2.08In the calculation process of the above three models, the GSIM (1,1) model extracts the most information, and the seasonal variation of the data is shown in the prediction process. The prediction accuracy of the GSIM (1,1) model is also higher than that of the GWO-Holt–Winters model and the DGGM (1,1) model. The MAPE of the GSIM (1,1) model is 2.08%, the MAPE of the Holt–Winters and the DGGM (1,1) models is 5% and 3.94%, respectively (Table8). This shows that the GSIM (1,1) model is more practical for seasonal fluctuation data.Table 8
The fitting error comparison of three models.
ModelMAPEHolt–Winters5DGGM (1,1)3.94GSIM (1,1)2.08The fitting process of the three models is visualized in Figure3. The consumption of groundwater in Ningxia keeps increasing, and seasonal fluctuation trend is obvious. The fitting line of the GSIM (1,1) model is close to the original data line. The three models did not fit the data of the first period. It can be seen that the original calculation of the three models has the same characteristics, and the three models are comparable.Figure 3
The fitting curve of groundwater depth in Yinchuan Plain.


#### 3.2.4. Prediction of Groundwater Depth in Yinchuan Plain

In the fitting results of groundwater depth in Yinchuan Plain, the GSIM (1,1) model has the smallest error. Therefore, the GSIM (1,1) model is selected to predict the groundwater depth of Yinchuan Plain from 2019 to 2020 (Table9). From the prediction results, groundwater consumption in Yinchuan Plain will increase further. The seasonal law of groundwater depth remains unchanged. Ningxia is inland in Northwest China. The inland areas are short of water. Yinchuan Plain is the most economically dynamic area in Ningxia. Continued economic growth has increased water consumption. In addition, the tradition farmland irrigation model consumes most of the water resources. Strengthening water resources management is crucial to the economic development and agricultural production of Yinchuan Plain.Table 9
Prediction of groundwater depth in Yinchuan Plain (unit: m).
MonthYearMonthYear201920202019202012.603.0472.603.1222.863.1282.232.6533.193.5393.244.3743.443.92103.103.4053.764.40113.423.7462.793.30122.222.65The GSIM (1,1) model is the best among three models in Yinchuan Plain. This model can also predict the groundwater depth of Weining Plain. These results in Weining Plain and Yinchuan Plain help us to know the groundwater situation in the whole Ningxia Plain. The prediction of groundwater depth in Weining Plain is described in the next section.


### 3.3. Forecasting the Groundwater Depth of Weining Plain

Yinchuan, the capital of Ningxia, is located on the Yinchuan Plain. Zhongwei city is located on the Weining Plain. The Yellow River enters Ningxia from Zhongwei city. From a geographical point of view, Weining Plain has more water resources than Yinchuan Plain. The Yellow River flows from Weining Plain to Yinchuan Plain. Weining Plain needs to coordinate the local water resources to support the water shortage of Yinchuan Plain. Therefore, the accurate prediction of groundwater in Weining Plain is conductive not only to the regulation of local water resources but also to the diversion of water to Yinchuan Plain.The annual groundwater depth in Weining Plain changes slowly. The annual data of groundwater depth were fitted by the FGM (1,1) model. The fitting results are listed in Table10, and the MAPE of fitting value is 0.10% and has high precision. Since 2016, the groundwater depth of Weining Plain has gradually increased. It shows that water consumption in Weining Plain has increased in recent years.Table 10
Annual fitting value of groundwater depth in Weining Plain (unit: m).
YearActual valueFitting value201421.2421.24201521.5421.52201620.6920.66201721.1921.18201822.3122.28MAPE0.10The variation law of groundwater consumption in Weining Plain is known by seasonal index. The parameters generated by the GSIM (1,1) model during the fitting process are listed in Table11. The fitting results of the seasonal index are listed in Table 12, and the MAPE of fitting value is 2.23%. The groundwater consumption in Weining Plain is large from January to May and November to December. The seasonal index in December showed a downward trend. It indicates that the groundwater consumption in December was reduced. But the seasonal index in October was above 1 several times. Changes in groundwater consumption are also dynamic, and it is important to predict the next year based on the index of same season.Table 11
The parameters of the GSIM (1,1) model fitting process in Weining Plain.
Modeling the serial numberThe optimal order numberra^b^11.39−0.081.6220.520.210.9430.550.170.9140.520.170.8650.140.500.8161.59−0.181.2271.56−0.170.9781.77−0.240.9291.44−0.180.81100.420.441.08110.470.371.13120.880.010.60Table 12
Seasonal index fitting value of Weining Plain groundwater depth.
MonthYear2014201520162017201811.201.321.101.051.0521.171.181.211.201.1731.211.201.241.261.2441.231.191.221.231.2251.171.181.191.181.1760.880.990.820.870.9970.770.780.670.720.8180.720.690.560.670.8290.790.700.690.760.87100.811.061.091.050.99111.041.171.181.131.06121.010.710.720.740.75MAPE2.23The MAPE of the GSIM (1,1) model for the groundwater depth of Weining Plain is marked in Figure4. The error is obtained by averaging the monthly data fitting error of a year. The original value of groundwater depth in Weining Plain is represented by the line in Figure 4. The seasonal variation of groundwater depth is obvious. The largest error occurred in 2016, which is 5.6%. The overall level of error is low.Figure 4
Annual fitting error value by the GSIM (1,1) model.From the fitting value of the year and the seasonal index, the fitting value of the month from 2014 to 2018 is calculated (Table13). The prediction section in Table 13 is visualized in Figure 5. It can be seen that the predicted value has a small increase compared with the original data. The groundwater depth in 2019 to 2020 will still fluctuate seasonally. The groundwater consumption in the first half is more than the second half. Government departments should focus on water management in the first half of the year.Table 13
The fitting value and prediction value of groundwater depth in Weining Plain (unit: m).
MonthFitting valuePrediction value201420152016201720182019202012.122.371.891.851.952.142.3922.072.112.082.132.182.242.3032.142.152.142.222.312.412.5142.182.132.112.182.272.372.4752.072.112.042.092.172.272.4061.561.781.401.541.842.282.8971.361.401.151.261.501.842.3181.271.240.971.171.532.042.7891.401.251.191.341.612.002.52101.431.891.881.851.831.821.83111.842.102.032.001.971.961.95121.791.271.251.311.401.511.65MAPE2.26Figure 5
The results of groundwater depth fitting and prediction in Weining Plain.The groundwater depth of Yinchuan Plain and Weining Plain is calculated and analyzed. It shows that the groundwater depth in Ningxia Plain is increasing year by year. Maintaining groundwater stability is important in dry areas. Water resources are coordinated between the two plains by the Yellow River. Weining Plain in the upper reaches of the Yellow River can allocate more water resources to Yinchuan Plain. The concept of coordinated development is also applied to water resources management.


## 3.1. Study Area and Data

Ningxia spans 35°14′–39°14′N and 104°17′–109°39′E and includes five cities. The climatic types of Ningxia are temperate, continental, arid, and semiarid. Ningxia Plain is considered as the largest plain in Ningxia. The study region has an area of 66.400 km2 and is bordered by the Helan Mountains and Liupan Mountains to the northwestward and south, respectively. Ningxia is 1100–1200 meters above sea level and the terrain slopes from southwest to northeast. The Yellow River enters from the southwest of Ningxia and leaves from the northeast. Due to the different physiognomy types in northcentral Ningxia and south Ningxia, Ningxia Plain can be divided into Yinchuan Plain and Weining Plain (Figure 1). The recoverable groundwater accounts for 88% of the total groundwater in Ningxia. Ningxia Plain is surrounded by mountains, water flows to the plain, and the water table changes slowly and steadily, creating a predictable gradient of water [43].Figure 1
The general area of Ningxia Plain.The data are from the official website of the Ministry of Water Resources of the People’s Republic of China (https://www.mwr.gov.cn/sj/tjgb/dxsdtyb/). There are a total of 47 groundwater monitoring wells in Ningxia Plain. Yinchuan Plain contains three cities and 28 groundwater monitoring wells. Weining Plain contains one city and 11 groundwater monitoring wells (by the end of the 2019). The groundwater depth data for Yinchuan Plain and Weining Plain are the average result of multiple groundwater monitoring wells.


## 3.2. Forecasting the Groundwater Depth of Yinchuan Plain

The groundwater depth of Yinchuan Plain is predicted in this section. The depth of groundwater is one of the contents of hydrology. Groundwater is subject to changes in rainfall or in rivers. As we all know, changes in rainfall and rivers are cyclical and seasonal. In recent years, the seasonal variation of groundwater depth is more obvious under the influence of human production activities. For example, the original groundwater depth data of Yinchuan Plain are summarized in Figure2. Figure 2 shows that the groundwater depth changes with time, and groundwater depth fluctuations are similar at the same time in different years. Therefore, the seasonal model is used to study the variation of groundwater depth.Figure 2
The groundwater depth of Yinchuan Plain and Weining Plain.


### 3.2.1. The Calculation Results of the Holt–Winters Model

The groundwater depth data of Yinchuan Plain from 2014 to 2018 were fitted to three models. The models do not fit the data of the first period, so the data of the first period are omitted. The fitting results of the Holt–Winters model are listed in Table1. The calculation process is completed by MATLAB 2016b. The GWO was used to find optimal parameters α,β, and γ, and the optimization result of parameters is [0.51639, 0.05160, 0.00752]. The fitting error (MAPE) of the Holt–Winters model is 5.00%. It can be seen that the fitting accuracy of the traditional model is lower than 10%, and the fitting results are accepted. As the first cycle, the data of 2014 were omitted.Table 1
The fitting value of Yinchuan Plain groundwater depth by Holt–Winters model (unit: m).
MonthYear201520162017201812.061.692.002.0922.272.192.482.6532.632.532.712.8942.692.622.752.9252.752.682.813.0761.941.902.052.2771.711.701.812.1081.551.591.692.0191.911.812.072.23102.462.182.582.94112.612.573.023.16121.531.561.791.85MAPE5.00


### 3.2.2. The Calculation Results of the DGGM (1,1) Model

The DGGM (1,1) model was applied to the data calculation on a monthly basis. The data grouping method is adopted to improve the fitting accuracy of the traditional GM (1,1) model, and the model can be applied to process seasonal fluctuation data. The fitting results of DGGM (1,1) model are listed in Table2. MAPE is 3.94%. Compared with the Holt–Winters model, the new model has higher fitting accuracy. Moreover, there are no fitting results for the first period of the two models, increasing contrast between models.Table 2
The fitting value of Yinchuan Plain groundwater depth by DGGM (1,1) model (unit: m).
MonthYear201520162017201812.031.972.112.1822.382.512.452.6232.562.742.792.9142.642.672.843.0952.562.783.073.2461.832.082.142.2871.641.671.892.0781.611.741.671.7891.821.932.092.26102.332.412.622.77112.452.672.913.05121.311.451.631.8MAPE3.94The parameters generated by the DGGM (1,1) model in the fitting process are shown in Table3. The time response function of each set of parameters can be obtained from equation (16).Table 3
The parameters of the DGGM (1,1) model fitting process.
Modeling the serial numberGrey development coefficienta^Grey action quantityb^1−0.0391.8382−0.0192.3813−0.0262.5394−0.0472.3955−0.0722.2796−0.0851.5747−0.1011.3378−0.0431.4959−0.0861.5310−0.0752.00311−0.0692.26312−0.1160.992


### 3.2.3. The Calculation Results of the GSIM (1,1) Model

The seasonal index in this paper refers to the index calculated on a monthly basis. The seasonal index reflects the fluctuation of the data in the time series. Most models do not deal well with fluctuating data. For the purpose of stabilizing the sequence, scholars used seasonal factors to eliminate seasonal changes in time series. This method improves the ability of the model to process the fluctuation data. However, the time series that eliminate seasonal effects may lose a lot of information. Its accuracy is also difficult to improve. To reduce information loss in the sequence, the seasonal index was extracted and studied separately. The results of yearly data are listed in Table4. The MAPE is 0.38%, and the fitting accuracy is high.Table 4
The annual fitting value of Yinchuan Plain groundwater depth (unit: m).
YearActual valueFitting value201426.6826.68201525.8325.81201625.7125.70201728.0227.70201830.6630.74MAPE0.38Then, the seasonal index of monthly data is fitted. During the fitting process, the GSIM (1,1) model is used sometimes, so the multiple sets of parameters are generated. Each set of parameters in Table5 can be substituted into equation (20) to obtain the corresponding time response function.Table 5
The parameters of the GSIM (1,1) model fitting process in Yinchuan Plain.
Modeling the serial numberThe optimal order numberra^b^11.21−0.051.0820.990.051.2430.150.891.4041.13−0.031.3551.10−0.041.2861.10−0.050.9071.11−0.070.7781.25−0.070.8791.58−0.191.11100.540.180.83110.570.170.95120.98−0.050.55The fitting values of the seasonal indices are listed in Table6. The average MAPE from 2015 to 2018 is 2.01%. Groundwater is consumed more in February, March, April, May, October, and November. Affected by dry and wet seasons, the groundwater depth varies seasonally throughout the year. The GSIM (1,1) model can not only predict accurately but also quantify the periodic variation of data within the time span. With the allowable error, it is beneficial to understand the variation law of groundwater depth.Table 6
Seasonal index fitting value of groundwater depth in Yinchuan Plain.
Month20142015201620172018Actual valueActual valueFitting valueActual valueFitting valueActual valueFitting valueActual valueFitting value10.920.960.960.890.890.890.880.880.8921.121.171.171.131.121.071.081.041.0331.231.251.251.241.221.171.181.141.1441.231.251.251.221.211.181.201.201.1951.221.241.241.231.231.251.251.271.2760.850.880.880.900.880.870.900.950.9370.730.770.770.800.780.780.810.870.8580.710.780.780.700.730.780.720.710.7490.810.920.920.770.770.870.840.940.96101.041.041.071.121.111.191.121.071.10111.201.181.201.251.241.281.241.201.22120.930.570.630.760.660.670.700.720.73MAPE2.01The groundwater depth fitting results from 2014 to 2018 are shown in Table7. The MAPE of fitting results is 2.08%. By calculating the seasonal index and annual fitting value of groundwater depth, the monthly fitting value is obtained. As the trend parameter, the annual fitting value controls the change direction of groundwater. The seasonal index breaks down the units of time from years into months.Table 7
The fitting value of Yinchuan Plain groundwater depth by GSIM (1,1) model (unit: m).
MonthYear2014201520162017201812.052.061.912.032.2722.492.522.412.492.6532.732.682.612.722.9242.732.692.602.763.0552.712.672.642.883.2661.891.891.892.082.3971.621.661.681.882.1881.581.681.551.671.9091.801.981.651.932.45102.312.292.382.582.83112.672.582.652.863.12122.071.351.421.611.88MAPE2.08In the calculation process of the above three models, the GSIM (1,1) model extracts the most information, and the seasonal variation of the data is shown in the prediction process. The prediction accuracy of the GSIM (1,1) model is also higher than that of the GWO-Holt–Winters model and the DGGM (1,1) model. The MAPE of the GSIM (1,1) model is 2.08%, the MAPE of the Holt–Winters and the DGGM (1,1) models is 5% and 3.94%, respectively (Table8). This shows that the GSIM (1,1) model is more practical for seasonal fluctuation data.Table 8
The fitting error comparison of three models.
ModelMAPEHolt–Winters5DGGM (1,1)3.94GSIM (1,1)2.08The fitting process of the three models is visualized in Figure3. The consumption of groundwater in Ningxia keeps increasing, and seasonal fluctuation trend is obvious. The fitting line of the GSIM (1,1) model is close to the original data line. The three models did not fit the data of the first period. It can be seen that the original calculation of the three models has the same characteristics, and the three models are comparable.Figure 3
The fitting curve of groundwater depth in Yinchuan Plain.


### 3.2.4. Prediction of Groundwater Depth in Yinchuan Plain

In the fitting results of groundwater depth in Yinchuan Plain, the GSIM (1,1) model has the smallest error. Therefore, the GSIM (1,1) model is selected to predict the groundwater depth of Yinchuan Plain from 2019 to 2020 (Table9). From the prediction results, groundwater consumption in Yinchuan Plain will increase further. The seasonal law of groundwater depth remains unchanged. Ningxia is inland in Northwest China. The inland areas are short of water. Yinchuan Plain is the most economically dynamic area in Ningxia. Continued economic growth has increased water consumption. In addition, the tradition farmland irrigation model consumes most of the water resources. Strengthening water resources management is crucial to the economic development and agricultural production of Yinchuan Plain.Table 9
Prediction of groundwater depth in Yinchuan Plain (unit: m).
MonthYearMonthYear201920202019202012.603.0472.603.1222.863.1282.232.6533.193.5393.244.3743.443.92103.103.4053.764.40113.423.7462.793.30122.222.65The GSIM (1,1) model is the best among three models in Yinchuan Plain. This model can also predict the groundwater depth of Weining Plain. These results in Weining Plain and Yinchuan Plain help us to know the groundwater situation in the whole Ningxia Plain. The prediction of groundwater depth in Weining Plain is described in the next section.


## 3.2.1. The Calculation Results of the Holt–Winters Model

The groundwater depth data of Yinchuan Plain from 2014 to 2018 were fitted to three models. The models do not fit the data of the first period, so the data of the first period are omitted. The fitting results of the Holt–Winters model are listed in Table1. The calculation process is completed by MATLAB 2016b. The GWO was used to find optimal parameters α,β, and γ, and the optimization result of parameters is [0.51639, 0.05160, 0.00752]. The fitting error (MAPE) of the Holt–Winters model is 5.00%. It can be seen that the fitting accuracy of the traditional model is lower than 10%, and the fitting results are accepted. As the first cycle, the data of 2014 were omitted.Table 1
The fitting value of Yinchuan Plain groundwater depth by Holt–Winters model (unit: m).
MonthYear201520162017201812.061.692.002.0922.272.192.482.6532.632.532.712.8942.692.622.752.9252.752.682.813.0761.941.902.052.2771.711.701.812.1081.551.591.692.0191.911.812.072.23102.462.182.582.94112.612.573.023.16121.531.561.791.85MAPE5.00


## 3.2.2. The Calculation Results of the DGGM (1,1) Model

The DGGM (1,1) model was applied to the data calculation on a monthly basis. The data grouping method is adopted to improve the fitting accuracy of the traditional GM (1,1) model, and the model can be applied to process seasonal fluctuation data. The fitting results of DGGM (1,1) model are listed in Table2. MAPE is 3.94%. Compared with the Holt–Winters model, the new model has higher fitting accuracy. Moreover, there are no fitting results for the first period of the two models, increasing contrast between models.Table 2
The fitting value of Yinchuan Plain groundwater depth by DGGM (1,1) model (unit: m).
MonthYear201520162017201812.031.972.112.1822.382.512.452.6232.562.742.792.9142.642.672.843.0952.562.783.073.2461.832.082.142.2871.641.671.892.0781.611.741.671.7891.821.932.092.26102.332.412.622.77112.452.672.913.05121.311.451.631.8MAPE3.94The parameters generated by the DGGM (1,1) model in the fitting process are shown in Table3. The time response function of each set of parameters can be obtained from equation (16).Table 3
The parameters of the DGGM (1,1) model fitting process.
Modeling the serial numberGrey development coefficienta^Grey action quantityb^1−0.0391.8382−0.0192.3813−0.0262.5394−0.0472.3955−0.0722.2796−0.0851.5747−0.1011.3378−0.0431.4959−0.0861.5310−0.0752.00311−0.0692.26312−0.1160.992


## 3.2.3. The Calculation Results of the GSIM (1,1) Model

The seasonal index in this paper refers to the index calculated on a monthly basis. The seasonal index reflects the fluctuation of the data in the time series. Most models do not deal well with fluctuating data. For the purpose of stabilizing the sequence, scholars used seasonal factors to eliminate seasonal changes in time series. This method improves the ability of the model to process the fluctuation data. However, the time series that eliminate seasonal effects may lose a lot of information. Its accuracy is also difficult to improve. To reduce information loss in the sequence, the seasonal index was extracted and studied separately. The results of yearly data are listed in Table4. The MAPE is 0.38%, and the fitting accuracy is high.Table 4
The annual fitting value of Yinchuan Plain groundwater depth (unit: m).
YearActual valueFitting value201426.6826.68201525.8325.81201625.7125.70201728.0227.70201830.6630.74MAPE0.38Then, the seasonal index of monthly data is fitted. During the fitting process, the GSIM (1,1) model is used sometimes, so the multiple sets of parameters are generated. Each set of parameters in Table5 can be substituted into equation (20) to obtain the corresponding time response function.Table 5
The parameters of the GSIM (1,1) model fitting process in Yinchuan Plain.
Modeling the serial numberThe optimal order numberra^b^11.21−0.051.0820.990.051.2430.150.891.4041.13−0.031.3551.10−0.041.2861.10−0.050.9071.11−0.070.7781.25−0.070.8791.58−0.191.11100.540.180.83110.570.170.95120.98−0.050.55The fitting values of the seasonal indices are listed in Table6. The average MAPE from 2015 to 2018 is 2.01%. Groundwater is consumed more in February, March, April, May, October, and November. Affected by dry and wet seasons, the groundwater depth varies seasonally throughout the year. The GSIM (1,1) model can not only predict accurately but also quantify the periodic variation of data within the time span. With the allowable error, it is beneficial to understand the variation law of groundwater depth.Table 6
Seasonal index fitting value of groundwater depth in Yinchuan Plain.
Month20142015201620172018Actual valueActual valueFitting valueActual valueFitting valueActual valueFitting valueActual valueFitting value10.920.960.960.890.890.890.880.880.8921.121.171.171.131.121.071.081.041.0331.231.251.251.241.221.171.181.141.1441.231.251.251.221.211.181.201.201.1951.221.241.241.231.231.251.251.271.2760.850.880.880.900.880.870.900.950.9370.730.770.770.800.780.780.810.870.8580.710.780.780.700.730.780.720.710.7490.810.920.920.770.770.870.840.940.96101.041.041.071.121.111.191.121.071.10111.201.181.201.251.241.281.241.201.22120.930.570.630.760.660.670.700.720.73MAPE2.01The groundwater depth fitting results from 2014 to 2018 are shown in Table7. The MAPE of fitting results is 2.08%. By calculating the seasonal index and annual fitting value of groundwater depth, the monthly fitting value is obtained. As the trend parameter, the annual fitting value controls the change direction of groundwater. The seasonal index breaks down the units of time from years into months.Table 7
The fitting value of Yinchuan Plain groundwater depth by GSIM (1,1) model (unit: m).
MonthYear2014201520162017201812.052.061.912.032.2722.492.522.412.492.6532.732.682.612.722.9242.732.692.602.763.0552.712.672.642.883.2661.891.891.892.082.3971.621.661.681.882.1881.581.681.551.671.9091.801.981.651.932.45102.312.292.382.582.83112.672.582.652.863.12122.071.351.421.611.88MAPE2.08In the calculation process of the above three models, the GSIM (1,1) model extracts the most information, and the seasonal variation of the data is shown in the prediction process. The prediction accuracy of the GSIM (1,1) model is also higher than that of the GWO-Holt–Winters model and the DGGM (1,1) model. The MAPE of the GSIM (1,1) model is 2.08%, the MAPE of the Holt–Winters and the DGGM (1,1) models is 5% and 3.94%, respectively (Table8). This shows that the GSIM (1,1) model is more practical for seasonal fluctuation data.Table 8
The fitting error comparison of three models.
ModelMAPEHolt–Winters5DGGM (1,1)3.94GSIM (1,1)2.08The fitting process of the three models is visualized in Figure3. The consumption of groundwater in Ningxia keeps increasing, and seasonal fluctuation trend is obvious. The fitting line of the GSIM (1,1) model is close to the original data line. The three models did not fit the data of the first period. It can be seen that the original calculation of the three models has the same characteristics, and the three models are comparable.Figure 3
The fitting curve of groundwater depth in Yinchuan Plain.


## 3.2.4. Prediction of Groundwater Depth in Yinchuan Plain

In the fitting results of groundwater depth in Yinchuan Plain, the GSIM (1,1) model has the smallest error. Therefore, the GSIM (1,1) model is selected to predict the groundwater depth of Yinchuan Plain from 2019 to 2020 (Table9). From the prediction results, groundwater consumption in Yinchuan Plain will increase further. The seasonal law of groundwater depth remains unchanged. Ningxia is inland in Northwest China. The inland areas are short of water. Yinchuan Plain is the most economically dynamic area in Ningxia. Continued economic growth has increased water consumption. In addition, the tradition farmland irrigation model consumes most of the water resources. Strengthening water resources management is crucial to the economic development and agricultural production of Yinchuan Plain.Table 9
Prediction of groundwater depth in Yinchuan Plain (unit: m).
MonthYearMonthYear201920202019202012.603.0472.603.1222.863.1282.232.6533.193.5393.244.3743.443.92103.103.4053.764.40113.423.7462.793.30122.222.65The GSIM (1,1) model is the best among three models in Yinchuan Plain. This model can also predict the groundwater depth of Weining Plain. These results in Weining Plain and Yinchuan Plain help us to know the groundwater situation in the whole Ningxia Plain. The prediction of groundwater depth in Weining Plain is described in the next section.


## 3.3. Forecasting the Groundwater Depth of Weining Plain

Yinchuan, the capital of Ningxia, is located on the Yinchuan Plain. Zhongwei city is located on the Weining Plain. The Yellow River enters Ningxia from Zhongwei city. From a geographical point of view, Weining Plain has more water resources than Yinchuan Plain. The Yellow River flows from Weining Plain to Yinchuan Plain. Weining Plain needs to coordinate the local water resources to support the water shortage of Yinchuan Plain. Therefore, the accurate prediction of groundwater in Weining Plain is conductive not only to the regulation of local water resources but also to the diversion of water to Yinchuan Plain.The annual groundwater depth in Weining Plain changes slowly. The annual data of groundwater depth were fitted by the FGM (1,1) model. The fitting results are listed in Table10, and the MAPE of fitting value is 0.10% and has high precision. Since 2016, the groundwater depth of Weining Plain has gradually increased. It shows that water consumption in Weining Plain has increased in recent years.Table 10
Annual fitting value of groundwater depth in Weining Plain (unit: m).
YearActual valueFitting value201421.2421.24201521.5421.52201620.6920.66201721.1921.18201822.3122.28MAPE0.10The variation law of groundwater consumption in Weining Plain is known by seasonal index. The parameters generated by the GSIM (1,1) model during the fitting process are listed in Table11. The fitting results of the seasonal index are listed in Table 12, and the MAPE of fitting value is 2.23%. The groundwater consumption in Weining Plain is large from January to May and November to December. The seasonal index in December showed a downward trend. It indicates that the groundwater consumption in December was reduced. But the seasonal index in October was above 1 several times. Changes in groundwater consumption are also dynamic, and it is important to predict the next year based on the index of same season.Table 11
The parameters of the GSIM (1,1) model fitting process in Weining Plain.
Modeling the serial numberThe optimal order numberra^b^11.39−0.081.6220.520.210.9430.550.170.9140.520.170.8650.140.500.8161.59−0.181.2271.56−0.170.9781.77−0.240.9291.44−0.180.81100.420.441.08110.470.371.13120.880.010.60Table 12
Seasonal index fitting value of Weining Plain groundwater depth.
MonthYear2014201520162017201811.201.321.101.051.0521.171.181.211.201.1731.211.201.241.261.2441.231.191.221.231.2251.171.181.191.181.1760.880.990.820.870.9970.770.780.670.720.8180.720.690.560.670.8290.790.700.690.760.87100.811.061.091.050.99111.041.171.181.131.06121.010.710.720.740.75MAPE2.23The MAPE of the GSIM (1,1) model for the groundwater depth of Weining Plain is marked in Figure4. The error is obtained by averaging the monthly data fitting error of a year. The original value of groundwater depth in Weining Plain is represented by the line in Figure 4. The seasonal variation of groundwater depth is obvious. The largest error occurred in 2016, which is 5.6%. The overall level of error is low.Figure 4
Annual fitting error value by the GSIM (1,1) model.From the fitting value of the year and the seasonal index, the fitting value of the month from 2014 to 2018 is calculated (Table13). The prediction section in Table 13 is visualized in Figure 5. It can be seen that the predicted value has a small increase compared with the original data. The groundwater depth in 2019 to 2020 will still fluctuate seasonally. The groundwater consumption in the first half is more than the second half. Government departments should focus on water management in the first half of the year.Table 13
The fitting value and prediction value of groundwater depth in Weining Plain (unit: m).
MonthFitting valuePrediction value201420152016201720182019202012.122.371.891.851.952.142.3922.072.112.082.132.182.242.3032.142.152.142.222.312.412.5142.182.132.112.182.272.372.4752.072.112.042.092.172.272.4061.561.781.401.541.842.282.8971.361.401.151.261.501.842.3181.271.240.971.171.532.042.7891.401.251.191.341.612.002.52101.431.891.881.851.831.821.83111.842.102.032.001.971.961.95121.791.271.251.311.401.511.65MAPE2.26Figure 5
The results of groundwater depth fitting and prediction in Weining Plain.The groundwater depth of Yinchuan Plain and Weining Plain is calculated and analyzed. It shows that the groundwater depth in Ningxia Plain is increasing year by year. Maintaining groundwater stability is important in dry areas. Water resources are coordinated between the two plains by the Yellow River. Weining Plain in the upper reaches of the Yellow River can allocate more water resources to Yinchuan Plain. The concept of coordinated development is also applied to water resources management.


## 4. Conclusions

In this paper, the GSIM (1,1) model, the Holt–Winters model, and the DGGM (1,1) model were established and applied to groundwater depth prediction in Ningxia Plain. By comparison, the GSIM (1,1) model is proved to be an effective one, and its prediction accuracy is higher than that of other models. The following conclusions were obtained.The GSIM (1,1) model is more suitable for dealing with seasonal fluctuation data. In the process of calculation, the seasonal variation of data is presented. The GSIM (1,1) model provides a new choice for studying the data with the characteristics of seasonal variation.GSIM (1,1) retains the characteristics of the traditional grey model for processing small sample data and has high prediction accuracy. Seasonal variations in the data were quantified and used for the first time to predict research. Data for the same month for different years are used as the object of the GSIM (1,1) model, and the processing method of traditional series data is changed.In the paper, the GSIM (1,1) model parameterr is reserved for two decimal places. During the experiment, we found that the more the bits reserved for parameter r, the higher the accuracy of the model. Scholars can adjust parameter r according to the experimental needs to meet higher accuracy requirements.



---

*Source: 6872538-2021-08-20.xml*