Next Article in Journal
Wind Power Generation Scheduling Accuracy in Europe: An Overview of ENTSO-E Countries
Previous Article in Journal
Submodule Fault-Tolerant Strategy for Modular Multilevel Converter with Scalable Control Structure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Grey Multivariable Prediction Model of Energy Consumption with Different Fractional Orders

College of Economics and Management, Nanjing University of Aeronautics and Astronautics, Nanjing 211100, China
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(24), 16426; https://doi.org/10.3390/su142416426
Submission received: 6 September 2022 / Revised: 1 December 2022 / Accepted: 6 December 2022 / Published: 8 December 2022
(This article belongs to the Section Energy Sustainability)

Abstract

:
The scientific prediction of energy consumption plays an essential role in grasping trends in energy consumption and optimizing energy structures. Energy consumption will be affected by many factors. In this paper, in order to improve the accuracy of the prediction model, the grey correlation analysis method is used to analyze the relevant factors. First, the factor with the largest correlation degree is selected, and then a new grey multivariable convolution prediction model with dual orders is established. Different fractional orders are used to accumulate the target data sequence and the influencing-factor data sequence, and the model is optimized by particle swarm optimization algorithm. The model is used to fit and test the energy consumption of Shanghai, Guizhou and Shandong provinces in China from 2011 to 2020 compared with other multivariable grey prediction models. Experimental results with the MAPE and RMSPE measurements show that our improved model is reasonable and effective in energy consumption prediction. At the same time, the model is applied to forecast the energy consumption of the three regions from 2021 to 2025, providing reliable information for future energy distribution.

1. Introduction

With the fast development of the social economy and industry, the energy consumed by humankind is also growing. Currently, energy is the indispensable material basis for human survival and development, and it is also important for supporting the normal operation of society and economy. At present, China is the world’s second largest economy. Because of its rapidly developing economy, China’s energy consumption is also increasing year by year. However, the world’s resources are not unlimited. Excessive energy consumption will cause a certain degree of damage to the global ecological environment and is not conducive to the sustainable development of all countries. In order to make a more rational use of energy, energy consumption prediction is particularly important. Through the scientific prediction of energy consumption, we can predict future trends in energy use, providing the basis for relevant departments to formulate reasonable energy consumption plans, optimize energy structures, and promote the realization of social, economic and ecological sustainable development.
Recently, many domestic and foreign scholars have adopted different models and methods to study different kinds of energy consumption. Chaturvedi et al. [1] compared the performances of Indian’s Central Energy Authority’s (CEA) existing trend-based model, SARIMA, LSTM RNN and Facebook prophet in predicting monthly total energy demand and peak energy demand of India. Peng et al. [2] proposed an optimized long short-term memory-based model and applied it to the prediction of energy consumption. Rick et al. [3] proposed a deep learning approach based on LSTM, CNN and an auto-encoder to predict data with different time series lengths. Jin et al. [4] combined singular spectrum analysis (SSA) and parallel long- and short-term memory (PLSTM) neural networks to form a new hybrid artificial intelligence prediction model and predicted the power consumption of five families in the UK. Somu et al. [5] presented an energy consumption forecasting model based on long short-term memory networks and improved a sine cosine optimization algorithm to accurately predict building energy consumption, including the energy consumption of the teaching building in Indian Institute of Technology. Kazemzadeh et al. [6] proposed a hybrid forecasting method that combined the auto-regressive integrated moving average (ARIMA), an artificial neural network (ANN) and the proposed support vector regression technique, and used this method to forecast the annual peak load and total energy demand of Iran’s national electric energy system. Adedeji et al. [7] used a non-linear autoregressive neural network (NARNET) to predict the energy consumption of a university in South Africa. Zeng et al. [8] first identified 17 factors that affect the structure of Chinese energy consumption. On this basis, they used a copula function to develop a multi-factor dynamic support vector machine model to predict the advanced index of China’s energy consumption structure in the future. Hu et al. [9] developed a stacked hierarchy of reservoirs (DeepESN) for forecasting energy consumption and wind power generation by introducing the deep learning framework into the basic echo state network. Bedi et al. [10] proposed a memory network prediction model based on decomposition and automatic encoder integration, which can achieve an accurate prediction of power demand. These prediction models mentioned above are deep learning models and time series methods that need a large quantity of data for training to obtain more accurate prediction results.
However, in practice, energy consumption statistics are often small data samples. When the amount of sample data is not enough, the deep learning model loses its application value. The grey prediction model can play an important role because it requires fewer data samples. The grey prediction model is a key component of the grey system theory proposed by Professor Deng, which has wide applicability to the prediction problems of small samples, a lack of information and uncertainty [11]. Therefore, many scholars usually use a grey prediction model to solve the problem of energy prediction with less data.
Xiong et al. [12] proposed a new incomplete gamma grey forecasting model based on fractional-order accumulation and used this model to simulate and predict natural gas consumption in the Asia–Pacific region from 2008 to 2018. Tong et al. [13] used an extrapolation method to extend the background value of the GM (1,1) model, optimized the model with a simulated annealing algorithm, and applied the model to forecast the coal consumption of China, India and the United States in the next five years. Ding et al. [14] established an optimized structure-adaptative grey model to predict the nuclear energy consumption of China and America. Ma et al. [15] proposed a new fractional time-delayed grey model, which is optimized by the Grey Wolf Optimizer, and applied the model to forecast the coal and natural gas consumption in Chongqing. Wu et al. [16] combined fractional-order accumulation and a grey multivariable convolution model and established an FGMC (1,N) model that was used to forecast the power consumption of Shandong Province. Liu et al. [17] proposed an optimized nonlinear grey Bernoulli prediction model with weighted fractional accumulation generation operation and used five optimization algorithms to optimize the model’s parameters to predict the natural gas production of Germany, Italy, Canada and China. Wang et al. [18] introduced adjacent accumulation into the grey multivariable convolution model and established the AGMC (1,N) model to predict the energy consumption of 13 provinces (cities) in seven regions of China. Li et al. [19] selected a genetic algorithm to optimize the background value and established a new nonlinear multivariable Verhulst model to predict China’s oil consumption. Duan et al. [20] used a novel multivariable Verhulst grey model (MVGM (1,N)) to predict the coal consumption of the Inner Mongolia and Gansu provinces in China. Zhao et al. [21] studied the prediction of non-renewable energy consumption in the Asia–Pacific Economic Cooperation countries by using the adjacent accumulation discrete grey model. Chen et al. [22] proposed a modeling method based on a grey prediction model and cloud model to predict the energy consumption of Shanghai, Zhejiang, Jiangsu and Anhui provinces. Zhou et al. [23] established a discrete grey seasonal model by cycle accumulation generation and predicted monthly industrial electricity consumption and quarterly natural gas consumption. Zhou et al. [24] proposed a new optimized fractional grey Holt–Winters model (NOFGHW) and applied it to the prediction of monthly crude oil production and quarterly industrial electricity consumption. Cheng et al. [25] used an improved GM (1,N) model to simulate and predict China’s clean energy consumption. Zhang et al. [26] established a new incomplete gamma grey model (IGGM) to predict solar energy consumption in Japan. Wu et al. [27] proposed a novel nonlinear grey Bernoulli model with fractional-order accumulation (FANGBM (1,1) model) to predict China’s renewable energy consumption from 2016 to 2020. Wang et al. [28] proposed a new fractional time-delayed grey Bernoulli model for data with scarcity, complexity and nonlinear characteristics, and predicted renewable energy, crude oil and fossil fuels usage. Liu et al. [29] established a discrete fractional grey model with a time power term and made a reliable prediction of China’s natural gas consumption. Wang et al. [30] proposed a new fractional-order grey model FENBGM (1,1) to forecast China’s oil consumption.
From the above research, we can see that scholars have improved the prediction accuracy of the grey prediction model from different aspects, such as optimizing the background value of the model, proposing new accumulation operators, etc. Although scholars have made great progress in energy consumption prediction at present, there are still problems, such as single factors affecting energy consumption, model improvements mainly focusing on background value optimization, a lack of sequence smoothness improvements, etc. These problems restrict improvements in model prediction accuracy. Therefore, in order to make the prediction of energy consumption more scientific and reasonable, considering that energy consumption is affected by many factors, such as economic development, national policies, population and the living standards of residents, a new grey multivariable prediction model is proposed in this paper, which can provide an effective reference for the government to formulate reasonable energy use policies. Tien (2005) [31] improved the basic GM (1,N) model and proposed a grey multivariable convolution model (GMC (1,N)), which effectively improved the prediction accuracy. The basic GMC (1,N) model uses first-order accumulation to generate sequences, which cannot make full use of old and new information and has certain limitations. In order to predict energy consumption more accurately, this paper introduces the different-order accumulation generation sequence proposed by Yin [32] into the GMC (1,N) model and then constructs a new dual-order multivariable grey convolution prediction model (FGMC (1,N,2r)). At the same time, a particle swarm optimization algorithm is selected to optimize the two fractional orders, which further increases the prediction accuracy of the model.
The main contributions of this paper are as follows:
(1)
In order to obtain the objective factors affecting energy consumption and establish a reasonable grey multivariable model, six factors affecting energy consumption are selected for grey relation analysis, from which the impact index with the largest correlation degree is selected;
(2)
On the basis of combining the different-order accumulation method and the grey multivariable convolution model, the grey multivariable convolution model based on different fractional-order accumulation is established. Different fractional orders are used to accumulate the target sequence and the influencing factor sequence, which can make full use of the information of different data sequences and reflect the influence of influencing factors on the change trend of the target sequence. Additionally, a particle swarm optimization algorithm is used to optimize fractional orders r 1 and r 2 . The model is used to fit, test and forecast the energy consumption of some provinces (cities) in China, which verifies the effectiveness of the model.
The other main parts of this paper are as follows: In Section 2, the relevant methods involved in the paper are introduced. In Section 3, the grey multivariable convolution model based on different fractional-order accumulation is proposed and the optimization method of the model parameters and model inspection index are given. In Section 4, we compare and discuss the energy consumption prediction results of the three provinces and use the model to simulate and test the energy consumption of other provinces (cities). In Section 5, the conclusions and suggestions of this paper are drawn.

2. Related Methods

2.1. Grey Relation Analysis

The main idea of grey relation analysis is to transform the system’s data series of influencing factors into sequence polylines and then construct a model to measure the degree of relevance on the basis of the geometric characteristics of the polylines. The more similar the geometric shape of the broken line is, the greater the relational degree between the corresponding sequences, and vice versa. Different from the traditional mathematical statistics method, it has no strict requirements on the amount and law of data, and the amount of calculation is small. Therefore, it is scientific and reasonable to choose the grey relation analysis method to analyze the impact of economic conditions, population, environment and other factors on energy consumption. In this paper, Deng’s grey correlation degree [33] is used to calculate the correlation degree between relevant factors and energy consumption. The specific calculation steps are as described below.
Step 1: Assuming that X 0 = ( x 0 ( 1 ) , x 0 ( 2 ) , , x 0 ( n ) ) is the sequence of dependent variable (target variable) and X i = ( x i ( 1 ) , x i ( 2 ) , , x i ( n ) ) is the sequence of independent variables (influencing factors).
Step 2: The sequence X 0 and sequence X i are processed dimensionless, and the initial value images of each sequence are obtained by Equations (1) and (2).
X 0 = { x 0 ( k ) x 0 ( 1 ) , k = 1 , 2 , , n }
X i = { x i ( k ) x i ( 1 ) , k = 1 , 2 , , n }
Step 3: Use the Equation (3) to calculate the absolute difference sequence of X 0 and X i .
Δ i ( k ) = { x 0 ( k ) - x i ( k ) , k = 1 , 2 , , n }
Step 4: Finding that the maximum value of Δ i ( k ) is M = m a x i m a x k Δ i ( k ) and the minimum value of Δ i ( k ) is m = m i n i m i n k Δ i ( k ) .
Step 5: Use the Equation (4) to calculate the correlation coefficient and use Equation (5) to calculate Deng’s grey correlation degree.
γ 0 i ( k ) = m + ξ M Δ i ( k ) + ξ M , k = 1 , 2 , , n
γ 0 i = 1 n k = 1 n γ 0 i ( k ) , k = 1 , 2 , , n

2.2. Particle Swarm Optimization

Particle swarm optimization (PSO) is a classical swarm intelligence algorithm, which originates from research on the predatory behavior of birds. The core of the algorithm is to find the optimal solution of the problem by using cooperation and information sharing among individuals in the group [34]. The specific description of the particle swarm optimization algorithm is as follows: First, a group of particles is initialized, assuming that there are m particles in the D - d i m e n s i o n a l search space, where the position of the i t h particle in the D - d i m e n s i o n a l space is expressed as X i = ( x i 1 , x i 2 , , x i D ) and the velocity is expressed as V i = ( v i 1 , v i 2 , , v i D ) . We can obtain its fitness value by substituting the position of each particle into the objective function and determine the quality of the position according to the fitness value. The historical optimal position of the i t h particle is expressed as P i = ( p i 1 , p i 2 , , p i D ) and the historical optimal position of the global particle is P g = ( p g 1 , p g 2 , , p g D ) . The velocity and position of each generation of particles can be iteratively calculated according to Equations (6) and (7).
v i d ( t + 1 ) = w v i d ( t ) + c 1 r 1 ( p i d - x i d ( t ) ) + c 2 r 2 ( p g d - x i d ( t ) )
x i d ( t + 1 ) = x i d ( t ) + v i d ( t + 1 )
Among them, w is the inertia factor, and its value is non-negative. The bigger the value of w is, the more the particles focus on the global search, and the smaller the value of w is, the more the particles focus on the local search. c 1 and c 2 are the learning factors and constants, which determine the maximum step size of particles moving to the individual optimal position and the global optimal position. r 1 and r 2 are two random numbers in the range of [0, 1], which are used to enhance the randomness of the search.

3. FGMC (1,N,2r) Model

3.1. Establishment of the Model

The fractional order is one of the major factors that influence the results of the fractional grey multivariable prediction model. In many existing grey multivariable prediction models, the target sequence and influencing factor sequence are accumulated by the same order, and different data sequences often have different characteristics. Using the same accumulation operator ignores the differences between data sequences. Therefore, in this paper, the different fractional orders that are used to accumulate the target sequence and influencing factor sequence are introduced into the GMC (1,N) model to construct a FGMC (1,N,2r) model. The modeling process is as shown below.
Step 1: Assume the target sequence is X 1 ( 0 ) = ( x 1 ( 0 ) ( 1 ) , x 1 ( 0 ) ( 2 ) , , x 1 ( 0 ) ( n ) ) and the influencing factor sequence is X i ( 0 ) = ( x i ( 0 ) ( 1 ) , x i ( 0 ) ( 2 ) , , x i ( 0 ) ( n ) ) .
Step 2: X 1 ( r 1 ) = ( x 1 ( r 1 ) ( 1 ) , x 1 ( r 1 ) ( 2 ) , , x 1 ( r 1 ) ( n ) ) is the fractional-order-accumulated generation of the target sequence of X 1 ( 0 ) by using Equation (8). X i ( r 2 ) = ( x i ( r 2 ) ( 1 ) , x i ( r 2 ) ( 2 ) , , x i ( r 2 ) ( n ) ) is the fractional-order-accumulated generation of the influencing factor sequence of X i ( 0 ) by using Equation (9).
x 1 ( r 1 ) ( k ) = i = 1 k Γ ( r 1 + k - i ) Γ ( k - i + 1 ) Γ ( r 1 ) x 1 ( 0 ) ( i ) , k = 1 , 2 , , n
x i ( r 2 ) ( k ) = i = 1 k Γ ( r 2 + k - i ) Γ ( k - i + 1 ) Γ ( r 2 ) x i ( 0 ) ( i ) , k = 1 , 2 , , n
Step 3: The convolution grey model of double fractional orders can be established as d x 1 ( r 1 ) ( t ) d t + b 1 x 1 ( r 1 ) ( t ) = i = 2 N b i x i ( r 2 ) ( t ) + μ , where b 1 , b 2 , , b N , μ are parameters, which can be solved by the least square method according to Equation (10).
[ b 1 , b 2 , , b n , μ ] Τ = ( B Τ B ) - 1 B Τ Y
where, Y = x 1 ( r 1 ) ( 2 ) - x 1 ( r 1 ) ( 1 ) x 1 ( r 1 ) ( 3 ) - x 1 ( r 1 ) ( 2 ) x 1 ( r 1 ) ( n ) - x 1 ( r 1 ) ( n - 1 )
B = ( - x 1 ( r 1 ) ( 2 ) + x 1 ( r 1 ) ( 1 ) 2 x 2 ( r 2 ) ( 2 ) + x 2 ( r 2 ) ( 1 ) 2 x N ( r 2 ) ( 2 ) + x N ( r 2 ) ( 1 ) 2 1 - x 1 ( r 1 ) ( 3 ) + x 1 ( r 1 ) ( 2 ) 2 x 2 ( r 2 ) ( 3 ) + x 2 ( r 2 ) ( 2 ) 2 x N ( r 2 ) ( 3 ) + x N ( r 2 ) ( 2 ) 2 1 - x 1 ( r 1 ) ( n ) + x 1 ( r 1 ) ( n - 1 ) 2 x 2 ( r 2 ) ( n ) + x 2 ( r 2 ) ( n - 1 ) 2 x N ( r 2 ) ( n ) + x N ( r 2 ) ( n - 1 ) 2 1 )
Step 4: The time response function is obtained as Equation (11).
x 1 ( r 1 ) ( k ) = x 1 ( 0 ) ( 1 ) e - b 1 ( k - 1 ) + j = 2 k { e - b 1 ( k - j + 0.5 ) f ( j ) + f ( j - 1 ) 2 }
where f ( u ) = b 2 x 2 ( r 2 ) ( u ) + b 3 x 3 ( r 2 ) ( u ) + + b N x N ( r 2 ) ( u ) + μ .
Step 5: The final prediction sequence is obtained by Equation (12).
x 1 ( 0 ) ( k ) = i = 0 k - 1 ( - 1 ) i Γ ( r 1 + 1 ) Γ ( i + 1 ) Γ ( r 1 - i + 1 ) x 1 ( r 1 ) ( k )
In order to more intuitively represent the FGMC (1,N,2r) model, the process of using the model for prediction is shown in a flowchart, as shown in Figure 1.

3.2. Parameter Optimization

An accumulation generation operator is an approach that makes a grey system turn from grey to white by using data mining and preprocessing technology. For the grey convolution multivariable prediction model based on different fractional-order accumulations, the orders r 1 and r 2 are important parameters that affect the original data preprocessing effect and model performance and act on the target sequence and the influencing factor sequence, respectively. The values of these two parameters have a decisive effect on the accuracy of the model, and the prediction accuracy of the model will change with their changes. Therefore, we select a particle swarm optimization algorithm to find the optimal value of the parameters, so as to improve the prediction accuracy and reduce the prediction error of the model. This optimization problem can be expressed by Equation (13).
a r g m i n r 1 , r 2 1 n k = 1 n | x 1 ( 0 ) ( k ) - x 1 ( 0 ) ( k ) x 1 ( 0 ) ( k ) | × 100 % s . t . { r 1 , r 2 [ 0 , 2 ] [ b 1 , b 2 , , b n , μ ] Τ = ( B Τ B ) - 1 B Τ Y x 1 ( r 1 ) ( k ) = x 1 ( 0 ) ( 1 ) e - b 1 ( k - 1 ) + j = 2 k { e - b 1 ( k - j + 0.5 ) f ( j ) + f ( j - 1 ) 2 } f ( j ) = b 2 x 2 ( r 2 ) ( j ) + b 3 x 3 ( r 2 ) ( j ) + + b N x N ( r 2 ) ( j ) + μ x 1 ( 0 ) ( k ) = x 1 ( r 1 ) ( 1 - r 1 ) ( k ) - x 1 ( r 1 ) ( 1 - r 1 ) ( k - 1 )
The objective of the optimization problem is to find a set of parameters to minimize the MAPE value of the model. The constraints include the value range of the two parameter values and the relevant formulas involved in the modeling process of the model.

3.3. Model Evaluation Index

To verify different prediction models’ precision more reasonably, we choose three evaluation indexes, namely the absolute percentage error (APE), average absolute percentage error (MAPE) and root mean square percentage error (RMSPE), to evaluate the prediction accuracy of the model. The specific equations are shown in Equations (14)–(16). The accuracy tested by MAPE is shown in Table 1.
A P E = x 1 ( 0 ) ( k ) - x 1 ( 0 ) ( k ) x 1 ( 0 ) ( k ) × 100 %
M A P E = 1 n k = 1 n x 1 ( 0 ) ( k ) - x 1 ( 0 ) ( k ) x 1 ( 0 ) ( k ) × 100 %
R M S P E = 1 n k = 1 n ( x 1 ( 0 ) ( k ) - x 1 ( 0 ) ( k ) x 1 ( 0 ) ( k ) ) 2 × 100 %

4. Model Application

In order to improve energy efficiency and push forward the construction of a modern energy system, relevant policies focusing on green and low-carbon energy development have been formulated in all regions and provinces of China. In the process of promoting green energy consumption, the analysis and prediction of energy consumption in recent years is an essential step. Therefore, we selected Shanghai, Guizhou and Shandong because the development levels of these three regions are different and representative. The development level of Shanghai has basically reached the level of developed countries, the development level of Shandong is basically at the level of developing countries, and the development level of Guizhou is slightly lower than that of Shandong. The geographical locations of these three regions are shown in Figure 2. Then, we used the model proposed in this paper to fit and test the energy consumption of these three regions to forecast their energy consumption in the next few years.

4.1. Experimental Data Collections

Table 2 shows the total energy consumption data of Shanghai, Guizhou and Shandong from 2011 to 2020 that are used in this study. The energy consumption in Shanghai accounts for 2.51% of China’s energy consumption, that in Guizhou Province accounts for 2.04%, and that in Shandong Province accounts for 8.51% [35]. The data of the first eight years are used for model construction, and the data of the second two years are used for model testing.
According to the relevant literature, six factors that affect energy consumption are preliminarily selected in this paper, namely the added value of the primary industry, the added value of the secondary industry, the added value of the tertiary industry, the total population, the urban per capita disposable income, and the rural per capita disposable income. The added value of the three major industries constitutes GDP, which reflects the economic condition of the province. Economic development is an important factor affecting energy consumption, and rapid economic development promotes the growth of energy consumption. Population also affects energy consumption. The increase in population means the increase in demand for various resources and, accordingly, energy consumption will also increase. In addition, urban per capita disposable income and rural per capita disposable income are the main indicators to reflect the living standards of residents. The growth in the living standards of residents will increase the consumption levels of high-quality energy, such as gasoline, electricity and natural gas, thus affecting energy consumption. On this basis, we use the grey correlation analysis method to select the factors most related to energy consumption for modeling.
All data are from the statistical yearbooks of all provinces (cities) and the official website of the National Bureau of Statistics (http://www.stats.gov.cn/ accessed on 6 July 2022). The statistical official website of each province (city) and their corresponding statistical yearbook can be accessed:

4.2. Model Comparison and Analysis

4.2.1. Energy Consumption in Shanghai

The data series of energy consumption is the target sequence, which is recorded as X 1 ( 0 ) . The data series of the added value of the primary industry, the added value of the secondary industry, the added value of the tertiary industry, the total population, the urban per capita disposable income and the rural per capita disposable income are the influencing factor sequences, which are recorded as X 2 ( 0 ) , X 3 ( 0 ) , X 4 ( 0 ) , X 5 ( 0 ) , X 6 ( 0 ) , X 7 ( 0 ) ,respectively. According to the calculation formula in Section 2.1, the grey correlation degrees between the energy consumption of Shanghai and the six influencing factors are calculated. The calculated results are shown in Table 3.
It can be seen from Table 3 that the total population is the most significant factor affecting Shanghai’s energy consumption, with a correlation value of 0.9705. Therefore, on this basis, the FGMC (1,2,2r) model is built, and the PSO algorithm is used to optimize the values of the two fractional orders. Meanwhile, in order to better reflect the effectiveness and accuracy of the new model, the model is compared with four models, including the GM (1,N) model, GMC (1,N) model [24] optimized AGMC (1,N) model [18] and optimized FGMC (1,N) model [16]. The data sizes used in the five models involved in the comparison are identical, and a PSO algorithm is used to optimize the parameters of the models that need to be optimized.
Table 4 shows the modeling and testing results of energy consumption in Shanghai from 2011 to 2020 and gives the parameter values of each model after optimization. It can be seen from Table 4 that in the modeling stage, the MAPE of the five grey prediction models are 6.32%, 0.85%, 0.85%, 0.85%, and 0.40%, and the RMSPE are 8.16%, 1.09%, 1.27%, 1.27% and 0.49%. In the testing stage, the MAPE of the five models are 2.71%, 3.77%, 3.80%, 3.27% and 2.40%, and the RMSPE are 3.19%, 5.27%, 5.17%, 4.15% and 2.45%. From the two evaluation indexes, the grey multivariable convolution model based on different fractional orders presented in this paper has good results in modeling and testing and meets the model accuracy requirements, indicating that the model is excellent. In order to more clearly show the fitting and prediction effects of each model, the result comparison diagram is shown in Figure 3.
From the above analysis, we can see that the new model has a high prediction accuracy. Additionally, we can use this model to predict the energy consumption of Shanghai from 2021 to 2025. Meanwhile, the optimized FGM (1,1) model is used to predict the sequence of influencing factors. The final results are shown in Table 5. Table 5 shows that in the next five years, the energy consumption of Shanghai will follow a downward trend, which may be related to the slower growth rate of the population.

4.2.2. Energy Consumption in Guizhou

In the modeling and testing of energy consumption in Guizhou Province, the grey relation analysis method is also used to calculate the grey correlation degrees between energy consumption and the six influencing factors. The results are shown in Table 6. According to Table 6, the grey correlation degree between energy consumption and population in Guizhou Province is the largest, which is 0.8641. Therefore, the FGMC (1,2,2r) model is established with population as the influencing factor sequence, which is compared with the results of the other four grey models. The results are shown in Table 7 and the prediction results are shown in Figure 4.
Table 7 shows that in the modeling stage, the MAPE of the five grey prediction models are 6.94%, 6.83%, 1.68%, 2.16% and 1.40%, and the RMSPE are 8.99%, 7.61%, 2.82%, 2.61% and 2.07%. In the testing stage, the MAPE of the five models are 5.76%, 10.09%, 8.96%, 5.35% and 5.07%, and the RMSPE are 5.76%, 10.09%, 8.98%, 5.35% and 5.07%. In the modeling stage and the testing stage, the MAPE and RMSPE of the proposed model are the smallest of the five models, which shows that this model is superior to other models.
In addition, the FGMC (1,2,2r) model is used to continuously forecast the energy consumption of Guizhou Province from 2021 to 2025. The results are shown in Table 8. The results reflect that the energy consumption of Guizhou Province generally presents an upward trend. As population is the most influential factor on the energy consumption of Guizhou Province, the energy consumption of Guizhou Province will also increase as the population increasing. However, the increasing trend is relatively gentle and there will be no excessive consumption of energy, which is of great significance for achieving sustainable development.

4.2.3. Energy Consumption in Shandong

Similarly, when studying the energy consumption of Shandong Province, the grey correlation degrees between the energy consumption and the six influencing factors are calculated first, and the results are shown in Table 9. Table 9 shows that the added value of the primary industry has the greatest correlation with the energy consumption of Shandong Province, with a value of 0.9420. Therefore, taking the added value of the primary industry as the influencing factor sequence, the five grey models are established, and Table 10 and Figure 5 are obtained. Compared with other models, the new model has better performance in fitting and predicting the trends of data, and some data points can almost coincide with the actual value.
In order to analyze the development trend of energy consumption in Shandong Province in the next few years, we apply the FGMC (1,2,2r) model to predict the energy consumption of Shandong Province from 2021 to 2025. The results are shown in Table 11. According to the prediction results, the energy consumption of Shandong Province generally presents a downward trend. Shandong is a large agricultural province, and the added value of the primary industry has the greatest impact on energy consumption. In the future, the energy consumption of Shandong Province will be inversely proportional to the added value of the primary industry, probably because Shandong Province has carried out a green and low-carbon transformation of its primary industry, which has improved its energy utilization efficiency.
After analyzing and predicting the energy consumption of the three regions, we can draw the following conclusions:
(1)
Compared with other grey models, the FGMC (1,N,2r) model developed in this paper is more accurate in its predictions;
(2)
The new model has better prediction performance after parameter optimization;
(3)
The new model can make up for the shortcomings of other models in predicting energy consumption.

4.3. Further Discussion

The FGMC (1,N,2r) model is applied to fit and test the energy consumption of other provinces (cities) in China so as to achieve the purpose of testing the role of the new model in the prediction of energy consumption. The model parameter values of provinces (cities) are shown in Table 12. The results are shown in Table 13 and Table 14.
From these data, this model has produced good results in the prediction of energy consumption in these six regions. In the modeling process from 2011 to 2018, all MAPE values are less than 1%. In the testing process from 2019 to 2020, the MAPE values are also controlled below 4%. When the model effect is good, the MAPE value can reach 0.17%.
In addition, we found that the influencing factors of each region are different, and the correlation between the influence factors and energy consumption is also different, which is one of the reasons for the model’s different prediction accuracies. However, in general, according to the accuracy of the model tested by MAPE given in Section 3.3 (see Table 1), the newly proposed model is excellent in the prediction of energy consumption, and its prediction performance indicators exceed the existing latest methods.
The development trends in energy consumption in China’s provinces (cities) are different; therefore, the governments and energy departments in different regions will need to take different measures. An accurate prediction of energy consumption can provide useful information for decision makers, help them formulate appropriate policies according to local conditions, and promote the optimization of energy structure. It has important practical significance.

5. Conclusions and the Future Work

To accurately predict energy consumption in different regions, a novel grey multivariable convolution prediction model based on different fractional-order accumulation is proposed. The new grey model uses different fractional orders to accumulate the target sequence and the influencing factor sequence, fully mining the information of different data sequences, and effectively combines the fractional-order accumulation and the grey multivariable convolution model. The main conclusions of this study are as follows:
(1)
The novel model is used to forecast the energy consumption of Shanghai, Guizhou and Shandong provinces. Compared with the simulation results of different grey models, the proposed grey multivariable convolution prediction model based on different fractional-order accumulation has higher prediction accuracy than the other grey multivariable prediction models, including the GM (1,N) model, GMC (1,N) model, AGMC (1,N) model and FGMC (1,N) model. From the two measurements of MAPE and RMSPE, the new model has the lowest MAPE and RMSPE values, making it superior to the other models;
(2)
Different fractional orders are used to accumulate the target sequence and the influencing factor sequence, which makes full use of the information of the data sequence with different characteristics. Furthermore, the optimized fractional-order r 1 and r 2 also greatly improve the prediction accuracy of the model;
(3)
The new model is a successful optimization of the grey multivariable convolution prediction model. It has a good performance in predicting energy consumption and it also provides a certain reference for other energy prediction problems.
Moreover, the limits of the study are as follows:
(1)
Grey relation analysis is applied to find out the most influential factors on energy consumption, and only one influencing factor is used to establish a grey multivariable model. In the future, more influencing factors can be considered to establish the model, which may lead to more accurate prediction results;
(2)
In this paper, a particle swarm optimization algorithm is used to optimize the parameters of the model. In the future, other optimization algorithms can be considered for optimization, and other parameters of the model can also be considered for optimization, so as to improve the prediction accuracy of the model;
(3)
The energy consumption data of 2020 is used and the impact of COVID-19 on energy consumption is temporarily ignored. This may be an important research direction in the future.

Author Contributions

Conceptualization, Y.S. and F.Z.; investigation, Y.S.; methodology, Y.S.; software, Y.S.; supervision, Y.S. and F.Z.; validation, Y.S.; writing—original draft, Y.S.; writing—review and editing, Y.S. and F.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Chaturvedi, S.; Rajasekar, E.; Natarajan, S.; McCullen, N. A comparative assessment of SARIMA, LSTM RNN and Fb Prophet models to forecast total and peak monthly energy demand for India. Energy Policy 2022, 168, 113097. [Google Scholar] [CrossRef]
  2. Peng, L.; Wang, L.; Xia, D.; Gao, Q. Effective energy consumption forecasting using empirical wavelet transform and long short-term memory. Energy 2022, 238, 121756. [Google Scholar] [CrossRef]
  3. Rick, R.; Berton, R. Energy forecasting model based on CNN-LSTM-AE for many time series with unequal lengths. Eng. Appl. Artif. Intell. 2022, 113, 104998. [Google Scholar] [CrossRef]
  4. Jin, N.; Yang, F.; Mo, Y.; Zeng, Y.; Zhou, X.; Yan, K.; Ma, X. Highly accurate energy consumption forecasting model based on parallel LSTM neural networks. Adv. Eng. Inform. 2022, 51, 101442. [Google Scholar] [CrossRef]
  5. Somu, N.; Raman, G.M.R.; Ramamritham, K. A hybrid model for building energy consumption forecasting using long short term memory networks. Appl. Energy 2020, 261, 114131. [Google Scholar] [CrossRef]
  6. Kazemzadeh, M.R.; Amjadian, A.; Amraee, T. A hybrid data mining driven algorithm for long term electric peak load and energy demand forecasting. Energy 2020, 204, 117948. [Google Scholar] [CrossRef]
  7. Adedeji, P.A.; Akinlabi, S.; Ajayi, O.; Madushele, N. Non-linear autoregressive neural network (NARNET) with SSA filtering for a university energy consumption forecast. Procedia Manuf. 2019, 33, 176–183. [Google Scholar] [CrossRef]
  8. Zeng, S.; Su, B.; Zhang, M.; Gao, Y.; Liu, J.; Luo, S.; Tao, Q. Analysis and forecast of China’s energy consumption structure. Energy Policy 2021, 159, 112630. [Google Scholar] [CrossRef]
  9. Hu, H.; Wang, L.; Lv, S. Forecasting energy consumption and wind power generation using deep echo state network. Renew. Energy 2020, 154, 598–613. [Google Scholar] [CrossRef]
  10. Bedi, J.; Toshniwal, D. Energy load time-series forecast using decomposition and autoencoder integrated memory network. Appl. Soft Comput. 2020, 93, 106390. [Google Scholar] [CrossRef]
  11. Deng, J.L. Control problem of grey system. Syst. Control. Lett. 1982, 1, 288–294. [Google Scholar]
  12. Xiong, P.; Li, K.; Shu, H.; Wang, J. Forecast of natural gas consumption in the Asia-Pacific region using a fractional-order incomplete gamma grey model. Energy 2021, 237, 121533. [Google Scholar] [CrossRef]
  13. Tong, M.; Dong, J.; Luo, X.; Yin, D.; Duan, H. Coal consumption forecasting using an optimized grey model: The case of the world’s top three coal consumers. Energy 2022, 242, 122786. [Google Scholar] [CrossRef]
  14. Ding, S.; Tao, Z.; Zhang, H.; Li, Y. Forecasting nuclear energy consumption in China and America: An optimized structure-adaptative grey model. Energy 2022, 239, 121928. [Google Scholar] [CrossRef]
  15. Ma, X.; Mei, X.; Wu, W.; Ma, X.; Wu, W.; Nie, R.; Chi, P.; Zhang, Y. A novel fractional time delayed grey model with Grey Wolf Optimizer and its applications in forecasting the natural gas and coal consumption in Chongqing China. Energy 2019, 178, 487–507. [Google Scholar] [CrossRef]
  16. Wu, L.; Gao, X.; Xiao, Y.; Yang, Y.; Chen, X. Using a novel multi-variable grey model to forecast the electricity consumption of Shandong Province in China. Energy 2018, 157, 327–335. [Google Scholar] [CrossRef]
  17. Liu, C.; Lao, T.; Wu, W.; Xie, W.; Zhu, H. An optimized nonlinear grey Bernoulli prediction model and its application in natural gas production. Expert Syst. Appl. 2022, 194, 116448. [Google Scholar] [CrossRef]
  18. Wang, M.; Wang, W.; Wu, L. Application of a new grey multivariate forecasting model in the forecasting of energy consumption in 7 regions of China. Energy 2022, 243, 123024. [Google Scholar] [CrossRef]
  19. Li, H.; Liu, Y.; Luo, X.; Duan, H. A novel nonlinear multivariable Verhulst grey prediction model: A case study of oil consumption forecasting in China. Energy Rep. 2022, 8, 3424–3436. [Google Scholar] [CrossRef]
  20. Duan, H.; Luo, X. A novel multivariable grey prediction model and its application in forecasting coal consumption. ISA Trans. 2022, 120, 110–127. [Google Scholar] [CrossRef]
  21. Zhao, H.; Wu, L. Forecasting the non-renewable energy consumption by an adjacent accumulation grey model. J. Clean. Prod. 2020, 275, 124113. [Google Scholar] [CrossRef]
  22. Chen, H.; Yang, Z.; Peng, C.; Qi, K. Regional energy forecasting and risk assessment for energy security: New evidence from the Yangtze River Delta region in China. J. Clean. Prod. 2022, 361, 132235. [Google Scholar] [CrossRef]
  23. Zhou, W.; Pan, J.; Tao, H.; Ding, S.; Chen, L.; Zhao, X. A novel grey seasonal model based on cycle accumulation generation for forecasting energy consumption in China. Comput. Ind. Eng. 2022, 163, 107725. [Google Scholar] [CrossRef]
  24. Zhou, W.; Tao, H.; Jiang, H. Application of a Novel Optimized Fractional Grey Holt-Winters Model in Energy Forecasting. Sustainability 2022, 14, 3118. [Google Scholar] [CrossRef]
  25. Cheng, M.; Li, J.; Liu, Y.; Liu, Y. Forecasting Clean Energy Consumption in China by 2025: Using Improved Grey Model GM (1, N). Sustainability 2020, 12, 698. [Google Scholar] [CrossRef] [Green Version]
  26. Zhang, P.; Ma, X.; She, K. Forecasting Japan’s Solar Energy Consumption Using a Novel Incomplete Gamma Grey Model. Sustainability 2019, 11, 5921. [Google Scholar] [CrossRef] [Green Version]
  27. Wu, W.; Ma, X.; Zeng, B.; Wang, Y.; Cai, W. Forecasting short-term renewable energy consumption of China using a novel fractional nonlinear grey Bernoulli model. Renew. Energy 2019, 140, 70–87. [Google Scholar] [CrossRef]
  28. Wang, Y.; He, X.; Zhang, L.; Ma, X.; Wu, W.; Nie, R.; Chi, P.; Zhang, Y. A novel fractional time-delayed grey Bernoulli forecasting model and its application for the energy production and consumption prediction. Eng. Appl. Artif. Intell. 2022, 110, 104683. [Google Scholar] [CrossRef]
  29. Liu, C.; Wu, W.; Xie, W.; Chen, X.; Xie, H.; Wang, J. Forecasting natural gas consumption of China by using a novel fractional grey model with time power term. Energy Rep. 2021, 7, 788–797. [Google Scholar] [CrossRef]
  30. Wang, Y.; Zhang, Y.; Nie, R.; Chi, P.; He, X.; Zhang, L. A novel fractional grey forecasting model with variable weighted buffer operator and its application in forecasting China’s crude oil consumption. Petroleum 2022, 8, 139–157. [Google Scholar] [CrossRef]
  31. Tien, T. The indirect measurement of tensile strength of material by the grey prediction model GMC (1, n). Measurement Science & Technology 2005, 16, 1322. [Google Scholar]
  32. Yin, F.; Zeng, B. A novel multivariable grey prediction model with different accumulation orders and performance comparison. Appl. Math. Model. 2022, 109, 117–133. [Google Scholar] [CrossRef]
  33. Zeng, B.; Luo, C.; Liu, S.; Li, C. A novel multi-variable grey forecasting model and its application in forecasting the amount of motor vehicles in Beijing. Comput. Ind. Eng. 2016, 101, 479–489. [Google Scholar] [CrossRef]
  34. Du, B.; Huang, S.; Guo, J.; Tang, H.; Wang, L.; Zhou, S. Interval forecasting for urban water demand using PSO optimized KDE distribution and LSTM neural networks. Appl. Soft Comput. 2022, 122, 108875. [Google Scholar] [CrossRef]
  35. Statistics Bureau of the People’s Republic of China. China Statistical Yearbook; China Statistics Press: Beijing, China, 2021. [Google Scholar]
Figure 1. The prediction process of FGMC (1,N,2r) model.
Figure 1. The prediction process of FGMC (1,N,2r) model.
Sustainability 14 16426 g001
Figure 2. Geographical distribution map of the study provinces.
Figure 2. Geographical distribution map of the study provinces.
Sustainability 14 16426 g002
Figure 3. Comparison of prediction results of energy consumption in Shanghai.
Figure 3. Comparison of prediction results of energy consumption in Shanghai.
Sustainability 14 16426 g003
Figure 4. Comparison of prediction results of energy consumption in Guizhou.
Figure 4. Comparison of prediction results of energy consumption in Guizhou.
Sustainability 14 16426 g004
Figure 5. Comparison of prediction results of energy consumption in Shandong.
Figure 5. Comparison of prediction results of energy consumption in Shandong.
Sustainability 14 16426 g005
Table 1. The accuracy of the model tested by MAPE.
Table 1. The accuracy of the model tested by MAPE.
MAPEModel Accuracy
<10%Superior
10–20%Good
20–50%General
>50%Poor
Table 2. Energy consumption of provinces (cities) during 2011–2020 (unit:10000 tce).
Table 2. Energy consumption of provinces (cities) during 2011–2020 (unit:10000 tce).
Region2011201220132014201520162017201820192020
Shanghai10,489.0910,57310,890.3910,639.8610,930.5311,241.7311,381.8511,453.7311,696.4611,099.59
Guizhou7295.117873.578715.489015.189319.68831.279379.59656.7310,074.1410,111.12
Shandong31,211.832,686.734,234.935,362.639,331.640,137.940,097.740,580.541,390.041,826.8
Table 3. Grey correlation degree results of Shanghai.
Table 3. Grey correlation degree results of Shanghai.
γ ( X 1 ( 0 ) , X 2 ( 0 ) ) γ ( X 1 ( 0 ) , X 3 ( 0 ) ) γ ( X 1 ( 0 ) , X 4 ( 0 ) ) γ ( X 1 ( 0 ) , X 5 ( 0 ) ) γ ( X 1 ( 0 ) , X 6 ( 0 ) ) γ ( X 1 ( 0 ) , X 7 ( 0 ) )
Value0.86100.91410.57960.97050.63400.6178
Table 4. Modeling and testing results of energy consumption in Shanghai.
Table 4. Modeling and testing results of energy consumption in Shanghai.
YearActual
Values
GM (1,N)GMC (1,N) ξ = 0.0497
AGMC (1,N)
r = 1.0826
FGMC (1,N)
r 1 = 0.6319 , r 2 = 0.000045
FGMC (1,N,2r)
Predicted
Values
APE
(%)
Predicted
Values
APE
(%)
Predicted
Values
APE
(%)
Predicted
Values
APE
(%)
Predicted
Values
APE
(%)
Modeling
201110,489.0910,489.090.0010,489.090.0010,489.090.0010,489.090.0010,489.090.00
201210,573.009054.3114.3610,574.230.0110,585.160.1210,597.230.2310,619.530.44
201310,890.3912,264.7412.6210,697.671.7710,830.790.5510,564.992.9910,795.910.87
201410,639.8611,599.359.0210,819.351.6910,950.172.9210,727.330.8210,715.440.71
201510,930.5311,204.852.5110,954.860.2211,073.861.3110,926.540.0410,927.180.03
201611,241.7311,149.550.8211,108.431.1911,222.530.1711,120.721.0811,197.820.39
201711,381.8511,120.922.2911,280.230.8911,381.850.0011,301.990.7011,363.760.16
201811,453.7311,156.752.5911,472.870.1711,555.60.8911,465.270.1011,474.780.18
MAPE (%) 6.32 0.85 0.85 0.85 0.40
RMSPE (%) 8.16 1.09 1.27 1.27 0.49
Testing
201911,696.4611,182.474.3911,687.280.0811,730.340.2911,613.860.7111,468.481.95
202011,099.5911,213.881.0311,926.917.4511,911.167.3111,746.815.8311,416.892.86
MAPE (%) 2.71 3.77 3.80 3.27 2.40
RMSPE (%) 3.19 5.27 5.17 4.15 2.45
Table 5. Forecasting results of energy consumption from 2021 to 2025 in Shanghai.
Table 5. Forecasting results of energy consumption from 2021 to 2025 in Shanghai.
Year20212022202320242025
Energy consumption11,371.2911,346.4511,297.9911,238.5411,170.68
Table 6. Grey correlation degree results of Guizhou.
Table 6. Grey correlation degree results of Guizhou.
γ ( X 1 ( 0 ) , X 2 ( 0 ) ) γ ( X 1 ( 0 ) , X 3 ( 0 ) ) γ ( X 1 ( 0 ) , X 4 ( 0 ) ) γ ( X 1 ( 0 ) , X 5 ( 0 ) ) γ ( X 1 ( 0 ) , X 6 ( 0 ) ) γ ( X 1 ( 0 ) , X 7 ( 0 ) )
Value0.58150.67640.63560.86410.79270.6905
Table 7. Modeling and testing results of energy consumption in Guizhou.
Table 7. Modeling and testing results of energy consumption in Guizhou.
YearActual
Values
GM (1,N)GMC (1,N) ξ = 0.4576
AGMC (1,N)
r = 0.0673
FGMC (1,N)
r 1 = 1.3243 , r 2 = 0.2525
FGMC (1,N,2r)
Predicted
Values
APE
(%)
Predicted
Values
APE
(%)
Predicted
Values
APE
(%)
Predicted
Values
APE
(%)
Predicted
Values
APE
(%)
Modeling
20117295.117295.110.007295.110.007295.110.007295.110.007295.110.00
20127873.577317.127.076866.3012.797744.021.657917.720.567901.350.35
20138715.4810,323.8418.457983.648.408714.070.028445.983.098715.480.00
20149015.189988.0910.798429.426.509014.490.018819.332.178954.760.67
20159319.69564.682.638638.827.309319.580.009074.152.639096.532.39
20168831.279435.936.858772.810.669399.856.449250.544.759229.194.51
20179379.59443.280.688892.525.199402.420.249379.560.009379.500.00
20189656.739451.282.138983.056.989329.593.399470.201.939478.101.85
MAPE (%) 6.94 6.83 1.68 2.16 1.40
RMSPE (%) 8.99 7.61 2.82 2.61 2.07
Testing
201910,074.149501.605.689048.8310.189240.048.289531.695.389549.955.20
202010,111.129521.535.839099.2310.019136.709.649572.685.339612.064.94
MAPE (%) 5.76 10.09 8.96 5.35 5.07
RMSPE (%) 5.76 10.09 8.98 5.35 5.07
Table 8. Forecasting results of energy consumption from 2021 to 2025 in Guizhou.
Table 8. Forecasting results of energy consumption from 2021 to 2025 in Guizhou.
Year20212022202320242025
Energy consumption9666.289732.079790.989846.459899.17
Table 9. Grey correlation degree results of Shandong.
Table 9. Grey correlation degree results of Shandong.
γ ( X 1 ( 0 ) , X 2 ( 0 ) ) γ ( X 1 ( 0 ) , X 3 ( 0 ) ) γ ( X 1 ( 0 ) , X 4 ( 0 ) ) γ ( X 1 ( 0 ) , X 5 ( 0 ) ) γ ( X 1 ( 0 ) , X 6 ( 0 ) ) γ ( X 1 ( 0 ) , X 7 ( 0 ) )
Value0.94200.93200.58820.78320.69530.6491
Table 10. Modeling and testing results of energy consumption in Shandong.
Table 10. Modeling and testing results of energy consumption in Shandong.
YearActual
Values
GM (1,N)GMC (1,N) ξ = 0.0495
AGMC(1,N)
r = 0.2628
FGMC (1,N)
r 1 = 1.5668 , r 2 = 0.9899
FGMC (1,N,2r)
Predicted
Values
APE
(%)
Predicted
Values
APE
(%)
Predicted
Values
APE
(%)
Predicted
Values
APE
(%)
Predicted
Values
APE
(%)
Modeling
201131,211.831,211.80.0031,211.80.0031,211.80.0031,211.80.0031,211.80.00
201232,686.726,709.018.2930,446.76.8531,945.82.2731,987.32.1432,686.70.00
201334,234.939,516.415.4332,808.34.1734,234.70.0034,469.00.6833,725.31.49
201435,362.639,158.210.7335,260.60.2936,226.62.4436,615.43.5436,357.72.81
201539,331.640,074.21.8937,403.14.9037,919.53.5938,266.62.7138,621.01.81
201640,137.939,195.82.3538,799.33.3439,130.72.5139,410.21.8139,955.70.45
201740,097.739,145.62.3739,329.41.9239,797.90.7540,082.30.0440,283.70.46
201840,580.540,084.21.2239,841.91.8240,332.50.6140,483.70.2440,580.50.00
MAPE (%) 7.47 3.33 1.74 1.59 1.00
RMSPE (%) 10.03 3.90 2.11 2.01 1.41
Testing
201941,390.041,428.80.0940,693.21.6840,958.31.0440,758.11.5341,103.20.69
202041,826.843,426.13.8241,982.40.3741,775.50.1241,000.51.9841,927.30.24
MAPE (%) 1.96 1.03 0.58 1.75 0.47
RMSPE (%) 2.70 1.22 0.74 1.77 0.52
Table 11. Forecasting results of energy consumption from 2021 to 2025 in Shandong.
Table 11. Forecasting results of energy consumption from 2021 to 2025 in Shandong.
Year20212022202320242025
Energy consumption42,118.8241,567.8141,059.8340,449.8539,744.12
Table 12. Model parameter values of provinces (cities).
Table 12. Model parameter values of provinces (cities).
RegionZhejiangHebeiChongqingJiangxiShanxiJiangsu
r 1 1.028430.001200.785770.098141.026811.35040
r 2 0.000050.408850.003810.230740.613420.84192
Table 13. Modeling and testing results of energy consumption in Zhejiang, Hebei and Chongqing.
Table 13. Modeling and testing results of energy consumption in Zhejiang, Hebei and Chongqing.
YearZhejiang Hebei Chongqing
Actual ValuesPredicted
Values
APE (%)Actual ValuesPredicted
Values
APE (%)Actual ValuesPredicted
Values
APE (%)
201117,82717,827.000.0028,075.0328,075.030.005516.155516.150.00
201218,07618,080.340.0228,762.4728,459.261.055834.845824.190.18
201318,82018,582.541.2629,664.3829,433.660.786225.926264.250.62
201418,82619,103.411.4729,320.2130,192.212.976603.616589.010.22
201519,61019,609.220.0031,036.7330,844.900.626924.776883.650.59
201620,27620,284.580.0431,458.0531,397.550.197099.717103.780.06
201721,03020,992.060.1832,082.5631,841.150.757251.597251.580.00
201821,67521,642.270.1532,185.2432,232.100.157452.727452.720.00
MAPE (%) 0.45 0.93 0.24
201922,39322,399.590.0332,545.4332,573.610.097687.257688.520.02
202024,66023,227.845.8132,782.7632,863.650.257621.877921.133.93
MAPE (%) 2.92 0.17 1.97
Table 14. Modeling and testing results of energy consumption in Jiangxi, Shanxi and Jiangsu.
Table 14. Modeling and testing results of energy consumption in Jiangxi, Shanxi and Jiangsu.
YearJiangxi Shanxi Jiangsu
Actual ValuesPredicted ValuesAPE (%)Actual ValuesPredicted ValuesAPE (%)Actual ValuesPredicted ValuesAPE (%)
20116847.16847.100.009107.489107.480.0027,588.9727,588.970.00
20127148.37148.310.009914.539907.320.0728,849.8428,849.840.00
20137582.97580.860.0310,610.4810,610.480.0029,205.3829,116.500.30
20148055.48004.840.6311,222.4611,230.590.0729,863.0329,863.030.00
20158423.48392.660.3611,745.9311,724.410.1830,374.1430,536.530.53
20168730.18728.310.0212,146.4712,149.150.0231,209.7131,068.760.45
20178971.99021.680.5512,548.5212,540.250.0731,602.0931,476.990.40
20189285.79285.790.0012,900.3812,900.380.0031,635.2031,742.770.34
MAPE (%) 0.23 0.06 0.29
20199665.29428.892.4413,748.0613,334.611.0632,525.9731,886.421.97
20209808.69313.775.0413,512.2613,958.433.3032,672.4931,924.922.29
MAPE (%) 3.74 2.18 2.13
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, Y.; Zhang, F. Grey Multivariable Prediction Model of Energy Consumption with Different Fractional Orders. Sustainability 2022, 14, 16426. https://doi.org/10.3390/su142416426

AMA Style

Sun Y, Zhang F. Grey Multivariable Prediction Model of Energy Consumption with Different Fractional Orders. Sustainability. 2022; 14(24):16426. https://doi.org/10.3390/su142416426

Chicago/Turabian Style

Sun, Yijue, and Fenglin Zhang. 2022. "Grey Multivariable Prediction Model of Energy Consumption with Different Fractional Orders" Sustainability 14, no. 24: 16426. https://doi.org/10.3390/su142416426

APA Style

Sun, Y., & Zhang, F. (2022). Grey Multivariable Prediction Model of Energy Consumption with Different Fractional Orders. Sustainability, 14(24), 16426. https://doi.org/10.3390/su142416426

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