#### Abstract

Accurately forecasting China’s total electricity consumption is of great significance for the government in formulating sustainable economic development policies, especially, China as the largest total electricity consumption country in the world. The calculation method of the background value of the GM(1, 1) model is an important factor of unstable model performance. In this paper, an extrapolation method with variable weights was used for calculating the background value to eliminate the influence of the extreme values on the performance of the GM(1, 1) model, and the novel extrapolation-based grey prediction model called NEGM(1, 1) was proposed and optimized. The NEGM(1, 1) model was then used to simulate the total electricity consumption in China and found to outperform other grey models. Finally, the total electricity consumption of China from 2018 to 2025 was forecasted. The results show that China’s total electricity consumption will be expected to increase slightly, but the total is still very large. For this, some corresponding recommendations to ensure the effective supply of electricity in China are suggested.

#### 1. Introduction

With the rapid growth of the Chinese economy, the total electricity consumption is increasing year by year, and China is the largest energy consumer in the world currently. From China Statistical Yearbook 2019 [1], China’s total electricity consumption reached 6282.1 billion kilowatt-hours (billion kw/h) in 2017; compared with 2016 and 2007, the electricity consumption increased by 4.67% and 98.16%, respectively. Due to the increasing demand for electricity consumption in China, vigorously developing and efficiently using electricity to ensure the effective supply of energy is the inevitable choice in China. Therefore, the establishment of a corresponding mathematical model for accurately forecasting future electricity consumption is an important contribution to the formulation of China’s energy security strategy and sustainable development program.

Since the end of the last century, numerous models have been introduced for forecasting electricity consumption, such as the time series analysis model [2, 3], autoregressive integrated moving average (ARIMA) model [4, 5], support vector machine (SVM) [6, 7], Bayesian statistics [8], random forest [9], composite system method [10], artificial neural networks (ANNs) [11, 12], and deep learning models [13, 14]. However, these exhibit the limitations of requiring a great many high-quality data samples.

The grey system theory, proposed by Professor Deng [15] in 1982, plays a key role in processing uncertain problems with poor information and small samples. After more than 30 years of development, many scholars have applied it to industry, agriculture, society, economy, energy, and other fields [16–20]. The prediction problem of total electricity consumption is a complex uncertainty problem, which can be thought of as a grey system problem. On the one hand, electricity consumption is closely related to a series of uncertain factors, such as the degree of economic development, industrial structure, population, and distribution losses; on the other hand, as China’s statistical departments did not release annual data relating to the electricity consumption before 2000 years ago, the sample size is limited. However, the grey forecasting model has a good effect on such a “small sample, poor information” system.

The grey forecasting model is an important component of grey system theory, which has attracted attention to improving the stability of simulation performance and the reliability of prediction results by many scholars. The traditional GM(1, 1) model [16] is the most significant constituent of grey prediction. However, it can only simulate homogeneous exponential growth sequences. The main reason is that it is difficult for a sequence to completely satisfy the law of homogeneous exponential growth. Many studies were conducted on the improvement of the grey model mainly on the following four aspects: (1) converting the raw data for improvement of smoothness by some buffer operators and effective methods [21, 22]; (2) enhancing the computational methods of parameters by some intelligence algorithms to optimize the initial value and background value [23, 24]; (3) extending the structure of traditional grey models [25, 26] for improving traditional model compatibility with different data, such as non-homogeneous sequence or fluctuation sequence; and (4) combining the grey prediction models with other models, such as combining neural network and Markov model with the grey prediction model [27, 28].

The background value is an important consideration of the grey model [29]. Currently, the background value is usually taken as the mean value of and , which is defined as . Geometrically, the trapezoidal area with a straight edge is used instead of the trapezoidal area with a curved edge. This is the main cause of the low prediction accuracy of the model. Many scholars have put forward methods to improve the traditional method of calculation of grey model background value to increase modeling accuracy. These calculation methods can be roughly divided into two types: one is to define the background value as [30]. The parameter is an unknown interpolation coefficient, , which is designed to adjust the background value to enhance the performance of grey forecasting model by some intelligence algorithms. The other one is to reconstruct the background value by function integral [31] or numerical calculation method [32]. The above research studies on the background value greatly enriched the theoretical system of the grey prediction models. However, there exists a defect that the interpolation method was used to calculate the background value of these models. It can be determined from the above structures of background value that the extreme values and are important factors that affect the smoothness of the background value and even affect the performance of the model. Therefore, it is necessary to expand the calculating range of the background value to weaken the influence of extreme values on .

Inspired by the literature [33], a novel extrapolation method for calculating the background value of the grey forecasting model (NEGM(1, 1)) is proposed to predict China’s total electricity consumption in this study. The calculation process of background value for the new model is composed of three values with variable weights , , and . Compared with the traditional grey models, NEGM(1, 1) not only can significantly improve the smoothness effect of the background value but also can effectively weaken the influence of extreme values on the model’s performance. In addition, some classic grey models are included in NEGM(1, 1), and then, it can simulate the homogeneous exponential sequence, non-homogeneous exponential sequence, and linear function sequence unbiased.

As compared with the existing research, this paper is unique in the following two aspects:(i)A novel grey forecasting model, NEGM(1, 1), is proposed, which solves the problem that many grey prediction models cannot overcome the influence of extreme values and the data type is single.(ii)The total electricity consumption in China is simulated and predicted by using the NEGM(1, 1) model, and the measures to ensure the effective supply of electricity in China are put forward.

The rest of the paper is arranged as follows. We propose a novel grey prediction model NEGM(1, 1) in Section 2. Section 3 presents modeling the evaluation criteria. In Section 4, the NEGM(1, 1) model is used to compare simulation and errors with the other two grey models and predict the total electricity consumption of China. In Section 5, countermeasures and suggestions are put forward based on the prediction results of Section 4. Conclusions are drawn in Section 6.

#### 2. The Novel Extrapolation-Based Grey Prediction Model

In this section, the novel extrapolation-based grey prediction model called NEGM(1, 1) is proposed and optimized. The modeling mechanism and prediction function of the new model are deduced.

*Definition 1. *(see [16]). Assume that a non-negative raw sequence is , where , for . is called the 1-AGO (Accumulating Generation Operator) sequence of , i.e.,where

*Definition 2. *Assume that and are the same as in Definition 1 and are constants; when , the equationis the novel extrapolation-based grey model with a single variable and one first-order equation, NEGM(1, 1) for short.

According to the inverse process of Accumulating Generation Operator in Definition 1,There are some special cases of the NEGM(1, 1) model as follows:(i)When and in equation (3), the NEGM(1, 1) model reduces to the SAIGM model [34] with the following form:(ii)Setting in equation (5), the SAIGM model reduces to the NGM(1, 1, k) model [35] with the following form:(iii)Setting in equation (5), the SAIGM model reduces to the classic GM(1, 1) [16] with the following form:It is easy to verify that the NEGM(1,1) model can unbiasedly simulate the homogeneous exponential sequence, non-homogeneous exponential sequence, and linear function sequence according to the above structural features.

Let the set of data sequence and be the same as in Definition 1; the background value is written as , . Then, the least-squares estimation for of the NEGM(1,1) model satisfieswhereAccording to Definition 2 and equation (3),Then, we haveThat is,Arranging equation (12), we getLetThen, we getand equation (15) is called the unbiased estimation form of the NEGM(1, 1) model. According to Definition 1, the final restored expression is as follows:Because equation (15) is a second-order difference equation, it is complicated to deduce its time response function by the direct substitution method. In addition, since we only make a short-term prediction of China's total electricity demand, we can directly calculate and by equations (15) and (16). Notice that the main purpose of constructing the NEGM(1, 1) model is to simulate or forecast and , and the time response function of the NEGM (1, 1) model is not important for us.

#### 3. Optimization of System Parameters

The performance of one prediction model is evaluated by its simulation and prediction errors. To test the performance of NEGM(1, 1), the mathematical meanings of symbols in Table 1 are defined as follows:(i) is the simulation or prediction data of the real data .(ii) is the **r**elative simulation percentage error of (RSPE): .(iii) is the relative prediction percentage error of (RPPE): .(iv) is the mean relative simulation percentage error (MRSPE): .(v) is the mean relative prediction percentage error (MRPPE): .(vi) is the comprehensive mean relative percentage error (CMRPE): .

When using the NEGM(1,1) model to forecast the original data, we first need to determine the parameters of the new model and then use equations (8) and (14) to obtain the parameters and , respectively, and the prediction values can be calculated by equation (15). In summary, the following minimum problem of the square of MRSPE can be defined to obtain the optimal parameters .

Since problem (17) is the convex quadratic programming problem, the optimal values of can be computed by the nonlinear solver of MATLAB. Figure 1 gives the calculation flowchart from which we can clearly understand the modeling process.

#### 4. Forecasting China’s Total Electricity Consumption

To verify the effectiveness of the NEGM(1, 1) model in Section 3, we establish grey models for the total electricity consumption in China from 2000 to 2017 in Table 2. First, the original data to use in the empirical analysis are divided into two parts. The first part (years 2000–2015) is the training data used for modeling, and the other part (years 2016-2017) is used as a test set. In addition, two competing models, including SAIGM [34] and traditional GM(1, 1) [16], are employed for simulating and forecasting the total electricity consumption in China. Finally, the model which has the best performance is used to predict future values from 2018 to 2025.

From Table 2,

The detailed modeling process of the NEGM(1, 1) model contains four steps: parameter estimation, model construction, model performance test, and data prediction.

*Step 1. *Parameter estimation.

From Definition 1,The parameter sequence can be calculated by equation (8), and the unknown parameters of NEGM(1, 1) are determined by the constraint optimization problem (17). The parameters of the NEGM(1, 1) model are shown in Table 3.

*Step 2. *Model construction.

According to Table 3 and equation (14), the modeling parameters , and are 1.0088, 0.0006, 305.0882, and 633.0891, respectively. We obtainand . Equation (20) is just the unbiased estimation form of the NEGM(1, 1) model for forecasting the total electricity consumption in China.

*Step 3. *Model performance comparisons and tests.

The NEGM(1, 1) model can be constructed and the simulation and prediction data of the NEGM(1, 1) can be computed. Then, the relative simulation/prediction percentage errors (RSPE/RPPE), the mean RSPE/RPPE (MRSPE/MRPPE), and the comprehensive mean relative percentage errors (CMRPE) can be obtained from the simulation and prediction data. The calculation results are shown in Table 1.

As shown in Table 1, the proposed NEGM(1, 1) model, both in-sample and out-of-sample, has the smallest mean relative percentage error (MRPE). Specifically, the values of MRSPE and MRPPE are 1.9780% and 2.8027%, respectively.

To show a direct comparison of the performances of NEGM(1, 1), SAIGM, and GM(1, 1), the scatter diagrams with the simulation and prediction results and percentage errors of the three models in Table 1 are shown in Figures 2 and 3, respectively.

According to 2 and 3, it is easy to observe that the simulation-prediction curve of NEGM(1, 1) is the closest to the real data, which indicates that the new model has the best comprehensive performance among the above three models. Furthermore, the comprehensive mean relative percentage error (CMRPE) is only 2.32%, which is close to grade I of the grey model error level reference table [16]. Therefore, NEGM(1, 1) is more suitable for forecasting the total electricity consumption in China.

*Step 4. *Prediction.

According to equation (20), the total electricity consumption in China from 2018 to 2025 can be predicted, as shown in Table 4.

From Table 4, we can see that in 2025, China’s total electricity consumption is expected to reach 9736.9 (billion kw/h). Overall, China’s total electricity consumption shows an increasing trend year by year in the next eight years.

#### 5. Suggestions

According to the prediction results in Table 4, though the annual average growth rate of total electricity consumption in China is about 4.79%, the total amount is still large. It is a major concern of China that some measures need to be taken to ensure the effective supply of electricity and rapid development of the following aspects.

Firstly, in terms of technology, China has continuously optimized the structure of the industry and controlled the excessive growth trend of industries with high energy consumption and high pollution. Simultaneously, it has promoted the development of high-tech industries and increased the research of electricity storage core technologies.

Secondly, in terms of policies, the power company has continuously accelerated the research and implementation of the renewable energy generation grid-connecting policy to promote the construction of microgrid facilities such as photovoltaic solar power. In addition, the electricity-saving awareness of residents should be cultivated, and it can employ the proper use of price leverage to strive for electricity conservation.

Thirdly, in terms of management, the Chinese government has continuously deepened the reform of power supply measurement according to market demands and scientifically promoted the construction of renewable energy. Moreover, it gradually established a management mechanism on energy storage for the power grid to increase the proportion of renewable energy in peaking alternative trading.

#### 6. Conclusion

In order to weaken the influence of extreme values in the background value on the accuracy of forecasting models, a novel univariate grey model called NEGM(1, 1) was proposed based on the SAIGM model, which adopted an extrapolation method with variable weights to calculate the background value. The new model was applied in forecasting China’s total electricity consumption and compared with two other models, SAIGM and GM(1, 1). Results showed that the NEGM(1, 1) model outperformed the other models, which further verified that the new model enriched the theoretical framework of grey prediction models. Lastly, NEGM(1, 1) was used to forecast the total electricity consumption in China from 2018 to 2025, and corresponding policy suggestions were given based on the prediction results.

Although the performance of the NEGM(1, 1) model was better than that of SAIGM and GM(1, 1), there are still some limitations. Because NEGM(1, 1) is a univariate model, the influence of external factors on China’s total electricity consumption is not modeled. Therefore, by fully mining the ways of these factors affecting total electricity consumption and expanding NEGM(1, 1) to a multivariate model, the performance of prediction may be further improved. The establishment of a multivariable NEGM(1, N) model is our next research direction.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The author declares that there are no conflicts of interest.

#### Acknowledgments

This study was supported by the National Social Science Foundation of China (no. 17ZDA065), the Natural Science Foundation of Chongqing Municipal Education Commission, China (no. KJQN201901601), and the Science and Technology Research Project of Chongqing University of Education (no. KY201910B).