Next Article in Journal
A Values-Based Approach to Exploring Synergies between Livestock Farming and Landscape Conservation in Galicia (Spain)
Next Article in Special Issue
Going Public and Industrial Upgrading of Traditional Clusters in Developing Countries: Rethinking the Dynamics of the ‘Jinjiang Model’ in China
Previous Article in Journal
The Behavioral Response to Location Based Services: An Examination of the Influence of Social and Environmental Benefits, and Privacy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Forecasting the Energy Consumption of China’s Manufacturing Using a Homologous Grey Prediction Model

1
College of Business Planning, Chongqing Technology and Business University, Chongqing 400067, China
2
Chongqing Key Laboratory of Electronic Commerce & Supply Chain System, Chongqing Technology and Business University, Chongqing 400067, China
*
Author to whom correspondence should be addressed.
Sustainability 2017, 9(11), 1975; https://0-doi-org.brum.beds.ac.uk/10.3390/su9111975
Submission received: 13 September 2017 / Revised: 24 October 2017 / Accepted: 27 October 2017 / Published: 31 October 2017
(This article belongs to the Special Issue Transition from China-Made to China-Innovation )

Abstract

:
With the rapid development of China’s manufacturing, energy consumption has increased rapidly, and this has become a major bottleneck affecting the sustainable development of China’s economy. This paper deduces and constructs a homologous grey prediction model with one variable and one first order equation (HGEM(1,1)) for forecasting the total energy consumption of China’s manufacturing based on the Grey system theory. Both parameter estimation (PE) and the deduction of the final restored expression (FRE) of the HGEM(1,1) model are all from the time response expression of the whitenization differential equation, which solves the ‘non-homologous’ defects of PE and FRE with traditional grey prediction models. HGEM(1,1) has good performance and can unbiasedly simulate a homogeneous/non-homogeneous exponential function sequence and a linear function sequence. Then, the HGEM(1,1)model is used to simulate and forecast the total energy consumption of China’s energy manufacturing, and the results show that the comprehensive performance of this model is much better than that of the classic Grey Model with one variable and single order equation, GM(1,1) for short and the frequently-used Discrete Grey Model with one variable and single order equation, DGM(1,1) for short. Finally, we forecast the total energy consumption of China’s manufacturing industry during the years 2018–2024. The results show that the total energy consumption in China’s manufacturing is slowing down but is still too large. For this, some measures, such as optimizing the manufacturing structure and speeding up the development and promotion of energy saving and emission reduction technologies, to ensure the effective supply of energy in China’s manufacturing industry are suggested.

1. Introduction

Since the 1990s, China’s manufacturing industry has continued to develop at a high speed and has become the main driving force of the continued rapid development of China’s economy [1,2,3]. The average annual growth rate of China’s manufacturing industry was 15.8% and the total value of the output of China’s manufacturing industry increased by 407.4% during 1995–2006 [4,5,6,7]. Meanwhile, the average annual energy consumption of the manufacturing industry accounts for 56.7% of the total annual energy consumption in China. Clearly, the manufacturing sector is the foremost energy consumer in China [8,9]. Hence, energy is the basis for the survival and development of China’s manufacturing industry. The high energy consumption characteristics of the manufacturing industry has greatly increased the overall level of energy consumption in China and intensified the dependence of China’s economic development on energy [10,11].
A scientific forecast of the energy consumption of China’s manufacturing industry and targets for preventive control measures will be provided to ensure the effective supply of energy to China’s manufacturing industry, which is has a positive significance in promoting stable and healthy development and plays an important role in guaranteeing the sustainable development of China’s economy [12].
The energy consumption of the manufacturing industry is affected by many uncertain factors, such as industry structure, technology level, energy price, economic scale and national policy [13] and has the typical characteristic of uncertainty, that is ‘grey cause’ [14]. An econometric regression model (ERM) is an important and frequently-used prediction model, which operates under the premise of a large sample of data (not less than 30), and mainly by studying data statistical laws to find the functional relation among variables. When the size of sample data is small, or sometimes even with a large quantity of data, there might not be any statistical laws to be found, in these cases, an ERM cannot be used to forecast [15]. The Markov prediction model is effective in the prediction of the state of a process, but it is not suitable for medium and long-term predictions for a system [16]. The neural network model implements the mapping function from input to output, but ‘over fitting’ often results in poor prediction performance [17]. Since the grey prediction model has the advantage of single variable modelling, we will use it to solve the issue of the modelling and prediction of the energy consumption of China’s manufacturing industry [18,19,20].
The GM (1,1) model [21] is the first grey prediction model with a single variable and one order derivative, and its final restored expression shows a homogeneous exponential function. Hence, when a modelling sequence has the characteristic of approximately homogeneous exponential growth, the model shows better performance for simulation and prediction [22,23,24]. However, the real world is full of complexity and uncertainty, and a sequence with approximately exponential growth is only a special case. More systematic behaviour sequences exhibit the characteristic of approximately inhomogeneous exponential growth [25,26]. In this case, if the GM(1,1) model is used to simulate or forecast the approximate inhomogeneous exponential growth sequence, the inherent modelling mechanism and model structure will lead to unsatisfactory simulation and prediction accuracy. However, the GM(1,1) model estimates the model parameters via differential equations, and the model time response is derived by differential equations. Therefore, the model has the properties of partial differential (smooth) and partial differential (jump). The ‘inconsistency’ between parameter estimation and model expression leads to poor simulation and prediction performance even in the face of strictly homogeneous exponential sequences [27,28].
In this paper, a new grey prediction model is proposed to solve the prediction issue of the energy consumption of China’s manufacturing industry. The new model has better simulation and prediction performance than those of the other grey prediction models because it solves the ‘inconsistency’ defect between parameter estimation and model expression of the classical GM(1,1) model and can unbiasedly simulate a homogeneous exponential sequence, non-homogeneous exponential sequence, and linear function sequence.
The major contributions of the paper include two aspects, as follows:
(I) A new grey prediction model named homologous grey energy prediction model, HGEM(1,1) is proposed, which solves the ‘misplaced replacement’ issue of the classical GM(1,1) model; (II) The HGEM(1,1) model is used to simulate and forecast the total energy consumption of China’s energy manufacturing, and some measures are suggested to ensure the effective supply of energy in China’s manufacturing industry.
The remainder of the paper is organized as follows. In Section 2, we build a new grey prediction model, HGEM(1,1). In Section 3, we study the error-checking method for the HGEM(1,1) model. In Section 4, we use the HGEM(1,1) model to simulate and forecast the energy consumption of China’s manufacturing industry. In Section 5, we provide some countermeasures and suggestions in the field of energy-saving and emission-reduction around the prediction results in Section 4. Our conclusions are presented in Section 6.

2. Homologous Grey Energy Prediction Model

Definition 2.1.
Assume that T E ( 0 ) = ( t e ( 0 ) ( 1 ) , t e ( 0 ) ( 2 ) , , t e ( 0 ) ( n ) ) is the time series data of the total energy consumption of China’s manufacturing industry, then T E ( 1 ) = ( t e ( 1 ) ( 1 ) , t e ( 1 ) ( 2 ) , , t e ( 1 ) ( n ) ) is called the accumulating generation sequence [19] with one order of T E ( 0 ) , where
t e ( 1 ) ( k ) = i = 1 k t e ( 0 ) ( k ) , k = 1 , 2 , , n
Definition 2.2.
Assume T E ( 0 ) and T E ( 1 ) are stated as Definition 2.1, then
d t e ( 1 ) d t + a t e ( 1 ) = b t + c
is named the whitenization differential equation of the homologous grey energy prediction model, where a ^ = ( a , b , c ) T denotes the undetermined parameters.
Now, we deduce the solution of the differential equation d t e ( 1 ) ( t ) / d t + a t e ( 1 ) ( t ) = b t + c . The homogeneous equation of Equation (1) is as follows,
d t e ( 1 ) ( t ) d t + a t e ( 1 ) ( t ) = 0 d t e ( 1 ) ( t ) d t = a t e ( 1 ) ( t )
Then
ln | t e ( 1 ) | = a t + ln | C 1 |
The general solution of the homogeneous Equation (2)
t e ( 1 ) ( t ) = C 1 e a t
In Equation (3), we replace C 1 with u ( t ) based on the constant variation method, and let
t e ( 1 ) ( t ) = u ( t ) e a t
By derivation of t on both sides of Equation (4), we obtain
d t e ( 1 ) ( t ) d t = u ( t ) e a t a u ( t ) e a t
Substituting Equation (5) into Equation (1)
u ( t ) e a t a u ( t ) e a t = b t + c a t e ( 1 )
Since u ( t ) = ( b t + d ) e a t , then
u ( t ) = ( b t + c ) e a t d t = b t e a t d t + c e a t d t
that is
u ( t ) = b a ( t e a t 1 a e a t ) + h a e a t + C = b a t e a t b a 2 e a t + h a e a t + C
Substituting Equation (6) into Equation (4)
t e ( 1 ) ( t ) = ( b a t e a t b a 2 e a t + h a e a t + C ) e a t
By re-arranging Equation (7), we get
t e ( 1 ) ( t ) = b a t b a 2 + c a + C e a t
when t = 1 , we get
t e ( 1 ) ( 1 ) = C e a + b a b a 2 + c a
and the value of C is as follows,
C = t e ( 1 ) ( 1 ) b a + b a 2 c a e a
Substituting Equation (10) into Equation (8), we get
t e ( 1 ) ( t ) = ( t e ( 1 ) ( 1 ) b a + b a 2 c a ) e a ( t 1 ) + b a t b a 2 + c a
The final restored expression of Equation (11) is as follows,
t e ( 0 ) ( t ) = t e ( 1 ) ( t ) t e ( 1 ) ( t 1 ) = ( 1 e a ) ( t e ( 0 ) ( 1 ) b a + b a 2 c a ) e a ( t 1 ) + b a
where t = 2 , 3 , , n , . In Equation (12), when t = 2 , 3 , , n , t ^ e ( 0 ) ( t ) is called the simulation data; when t = n + 1 , n + 2 , , t ^ e ( 0 ) ( t ) is called the prediction data.
From Equation (12),
t e ( 1 ) ( t + 1 ) = ( t e ( 0 ) ( 1 ) b a + b a 2 c a ) e a t + b a ( t + 1 ) b a 2 + c a = e a ( t e ( 0 ) ( 1 ) b a + b a 2 c a ) e a ( t 1 ) + b a ( t + 1 ) b a 2 + c a = e a [ ( t e ( 0 ) ( 1 ) b a + b a 2 c a ) e a ( t 1 ) + b a t b a 2 + c a ]
that is
t e ( 1 ) ( t + 1 ) = e a ( b a t b a 2 + c a ) + b a ( t + 1 ) b a 2 + c a = e a t e ( 1 ) ( t ) + b a ( 1 e a ) t + ( 1 e a ) ( c a b a 2 ) + b a
Equation (13) is the solution of the differential equation d t e ( 1 ) ( t ) / d t + a t e ( 1 ) ( t ) = b t + c , and it is also called the time response function of the whitenization differential equation.
To solve the ‘misplaced replacement’ issue of the classical GM(1,1) model, we use the ordinary least-squares (OLS) method and Cramer’s rule to estimate parameters a, b and c according to the time response function of the whitenization differential equation, that is Equation (13). Let
α = e a , β = b a ( 1 e a )   and   γ = ( 1 e a ) ( c a b a 2 ) + b a
Then, Equation (13) can be simplified as follows,
t e ( 1 ) ( t + 1 ) = α t e ( 1 ) ( t ) + β t + γ
We employ the OLS method and Cramer’s rule to estimate parameters α ^ , β ^ , γ ^ in Equation (14). After this, parameters a, b and c in Equation (13) can be calculated.
Assume that α ^ , β ^ , γ ^ are the parametric estimated values of Equation (14), replacing t e ( 1 ) ( t + 1 ) , t = 1 , 2 , , n 1 with simulation values t ^ e ( 1 ) ( t + 1 ) = α ^ t e ( 1 ) ( t ) + β ^ t + γ ^ , to minimise the simulation error, the following condition needs to be satisfied:
S = t = 1 n 1 [ t e ( 1 ) ( t + 1 ) α ^ t e ( 1 ) ( t ) β ^ t γ ^ ] 2
According to OLS, we minimise S with respect to parameters α ^ , β ^ and γ ^ to obtain
{ S α ^ = 2 t = 1 n 1 [ t e ( 1 ) ( t + 1 ) α ^ t e ( 1 ) ( t ) β ^ t γ ^ ] t e ( 1 ) ( t ) = 0 S β ^ = 2 t = 1 n 1 [ t e ( 1 ) ( t + 1 ) α ^ t e ( 1 ) ( t ) β ^ t γ ^ ] t = 0 S γ ^ = 2 t = 1 n 1 [ t e ( 1 ) ( t + 1 ) α ^ t e ( 1 ) ( t ) β ^ t γ ^ ] = 0
From the above formulas, an equation set can be obtained:
{ α ^ t = 1 n 1 t e ( 1 ) ( t ) 2 + β ^ t = 1 n 1 t t e ( 1 ) ( t ) + γ ^ t = 1 n 1 t e ( 1 ) ( t ) = t = 1 n 1 t e ( 1 ) ( t + 1 ) t e ( 1 ) ( t ) α ^ t = 1 n 1 t t e ( 1 ) ( t ) + β ^ t = 1 n 1 t 2 + γ ^ t = 1 n 1 t = t = 1 n 1 t t e ( 1 ) ( t + 1 ) α ^ t = 1 n 1 t e ( 1 ) ( t ) + β ^ t = 1 n 1 t + γ ^ ( n 1 ) = t = 1 n 1 t e ( 1 ) ( t + 1 )
Next, the calculation of unknown parameters α ^ , β ^ and γ ^ in Equation Set (2.15) is presented. According to Cramer’s rule, we can obtain the following results:
B = | t = 1 n 1 t e ( 1 ) ( t ) 2 t = 1 n 1 t t e ( 1 ) ( t ) t = 1 n 1 t e ( 1 ) ( t ) t = 1 n 1 t t e ( 1 ) ( t ) t = 1 n 1 t 2 t = 1 n 1 t t = 1 n 1 t e ( 1 ) ( t ) t = 1 n 1 t n 1 |
B 1 = | t = 1 n 1 t e ( 1 ) ( t + 1 ) t e ( 1 ) ( t ) t = 1 n 1 t t e ( 1 ) ( t ) t = 1 n 1 t e ( 1 ) ( t ) t = 1 n 1 t t e ( 1 ) ( t + 1 ) t = 1 n 1 t 2 t = 1 n 1 t t = 1 n 1 t e ( 1 ) ( t + 1 ) t = 1 n 1 t n 1 |
B 2 = | t = 1 n 1 t e ( 1 ) ( t ) 2 t = 1 n 1 t e ( 1 ) ( t + 1 ) t e ( 1 ) ( t ) t = 1 n 1 t e ( 1 ) ( t ) t = 1 n 1 t t e ( 1 ) ( t ) t = 1 n 1 t t e ( 1 ) ( t + 1 ) t = 1 n 1 t t = 1 n 1 t e ( 1 ) ( t ) t = 1 n 1 t e ( 1 ) ( t + 1 ) n 1 |
B 3 = | t = 1 n 1 t e ( 1 ) ( t ) 2 t = 1 n 1 t t e ( 1 ) ( t ) t = 1 n 1 t e ( 1 ) ( t + 1 ) t e ( 1 ) ( t ) t = 1 n 1 t t e ( 1 ) ( t ) t = 1 n 1 t 2 t = 1 n 1 t t e ( 1 ) ( t + 1 ) t = 1 n 1 t e ( 1 ) ( t ) t = 1 n 1 t t = 1 n 1 t e ( 1 ) ( t + 1 ) |
In light of Cramer’s rule, parameters α ^ , β ^ and γ ^ can be computed, as shown below:
α ^ = B 1 B , β ^ = B 2 B , γ ^ = B 3 B .
Since
α = e a , β = b a ( 1 e a ) , γ = ( 1 e a ) ( c a b a 2 ) + b a ,
parameters a, b and c can be obtained, as follows:
a ^ = ln α ^ ,   b ^ = a ^ β ^ 1 α ^   and   c ^ = a ^ γ ^ b ^ 1 α ^ + b ^ a ^
Substituting a, b and c into Equation (12), the new model is established.
Equation (12) is called the homologous grey energy prediction model with a single variable and one order derivate, HGEM(1,1) for short. Compared with the classical GM(1,1) model, the parameter estimation and model time response expression of HGEM(1,1) are derived from Equation (14), which ensures the consistency of their sources (homologous). The proposed HGEM(1,1) model can unbiasedly simulate a homogeneous exponential sequence, non-homogeneous exponential sequence, and linear function sequence (the detailed proof is omitted here), which shows it has good performance.

3. Error Checking Method for the HGEM(1,1) Model

A model’s performance includes two aspects: simulation performance and prediction performance. Normally, only the models that pass various tests can be meaningfully employed to make predictions.
Definition 3.1.
Assume a raw sequence
T E ( 0 ) = ( t e ( 0 ) ( 1 ) , t e ( 0 ) ( 2 ) , , t e ( 0 ) ( n ) , t e ( 0 ) ( n + 1 ) , , t e ( 0 ) ( n + p ) )
A subsequence composed of the first n elements of sequence T E ( 0 ) is used to build the HGEM(1,1) model, and simulation sequence S ^ ( 0 ) is as follows,
S ^ ( 0 ) = ( t ^ e ( 0 ) ( 1 ) , t ^ e ( 0 ) ( 2 ) , , t ^ e ( 0 ) ( n ) ) .
We use the HGEM(1,1) model to forecast the latter p -step data and the prediction sequence F ^ ( 0 ) is as follows:
F ^ ( 0 ) = ( t ^ e ( 0 ) ( n + 1 ) , t ^ e ( 0 ) ( n + 2 ) , , t ^ e ( 0 ) ( n + p ) ) .
The error sequences of S ^ ( 0 ) and F ^ ( 0 ) are ε s and ε F respectively, which as follows
ε S = ( ε S ( 1 ) , ε S ( 2 ) , , ε S ( n ) ) ,   ε F = ( ε F ( n + 1 ) , ε F ( n + 2 ) , , ε F ( n + p ) )
where
ε S ( u ) = | t e ( 0 ) ( u ) t ^ e ( 0 ) ( u ) | , u = 1 , 2 , , n
and
ε F ( v ) = | t e ( 0 ) ( v ) t ^ e ( 0 ) ( v ) | , v = n + 1 , n + 2 , , n + p
The relative simulation percentage error (RSPE) of the simulation sequence is Δ S
Δ S = ( Δ S ( 1 ) , Δ S ( 2 ) , , Δ S ( n ) )
where
Δ S ( u ) = | ε S ( u ) t e ( 0 ) ( u ) × 100 % | , u = 1 , 2 , , n
The mean relative simulation percentage error (MRSPE) of simulation sequence Δ ¯ S is as follows:
Δ ¯ S = 1 n u = 1 n Δ S ( u ) .
The relative prediction percentage error (RPPE) of prediction sequence Δ F is as follows:
Δ F = ( Δ F ( n + 1 ) , Δ F ( n + 2 ) , , Δ F ( n + p ) )
where
Δ F ( v ) = | ε F ( v ) t ^ e ( 0 ) ( v ) × 100 % | , v = n + 1 , n + 2 , , n + p .
The mean relative prediction percentage error (MRPPE) of prediction sequence Δ ¯ F is as follows:
Δ ¯ F = 1 p v = n + 1 n + p Δ F ( v ) .
The comprehensive mean relative percentage error (CMRPE) of the HGEM(1,1) model is Δ
Δ = Δ ¯ S + Δ ¯ F 2
For giving α and β when Δ ¯ S < α and Δ ¯ F < β hold true, the HGEM(1,1) model is said to be error-satisfactory. However, when the size of modelling data is small, the original sequence cannot be divided into “simulation subsequence” and “prediction subsequence”. At this time, we only test the simulated error of the model, and the test of prediction error is omitted.
The modelling flowchart of the HGEM(1,1) model can be seen in Figure 1, as follows.

4. Forecasting the Energy Consumption of China’s Manufacturing with HGEM(1,1)

In recent years, China’s manufacturing industry has risen rapidly and has become a major manufacturing country in the world. However, the development of the manufacturing industry is mainly based on the consumption of large amounts of coal, petroleum and other non-renewable energy sources, which results in an energy shortage in China. China’s manufacturing industry consumes 54% of China’s total energy consumption. Hence, it is of great theoretical and practical significance to forecast the energy consumption situation of China’s manufacturing industry in the future. Such a forecast would be instrumental to ensuring the effective supply of energy in China’s manufacturing industry and promote the sustainable development of China’s economy. The energy consumptions of China’s manufacturing (CMEC) during 2006–2012 are as shown in Table 1, as follows.
From Table 1
T E ( 0 ) = ( t e ( 0 ) ( 1 ) , t e ( 0 ) ( 2 ) , t e ( 0 ) ( 3 ) , t e ( 0 ) ( 4 ) , t e ( 0 ) ( 5 ) , t e ( 0 ) ( 6 ) , t e ( 0 ) ( 7 ) ) = ( 143051.46 , 156218.80 , 172106.52 , 180595.96 , 88497.85 , 200403.37 , 205667.69 )
Data t e ( 0 ) ( 1 ) ~ t e ( 0 ) ( 6 ) are applied to build the HGEM(1,1) model of the energy consumption of China’s manufacturing; and the last data t e ( 0 ) ( 7 ) is used to check the prediction performance of the HGEM(1,1) model.
The detailed modelling process of the HGEM(1,1) model contains four steps, which are parameter estimation, model construction, model performance test and data prediction, as follows.
Step 1 Parameter estimation
From Definition 2.1,
T E ( 1 ) = ( t e ( 1 ) ( 1 ) , t e ( 1 ) ( 2 ) , t e ( 1 ) ( 3 ) , t e ( 1 ) ( 4 ) , t e ( 1 ) ( 5 ) , t e ( 1 ) ( 6 ) ) = ( 143051.460 , 299270.260 , 471376.780 , 651972.740 , 840470.590 , 1040873.960 )
Constructing determinants B , B 1 , B 2 , B 3 from T E ( 1 ) , the parameters of the HGEM(1,1) can be estimated by Cramer’s rule as follows:
α ^ = B 1 B = 0.861037 β ^ = B 2 B = 34760.480867 γ ^ = B 3 B = 142156.224711
That is
a ^ = ln α ^ = 0.14962 b ^ = a ^ β ^ 1 α ^ = 37425.70438 c ^ = a ^ γ ^ b ^ 1 α ^ + b ^ a ^ = 133876.59532
Step 2 Model construction
Substituting a ^ = ( a , b , c ) T = ( 0.14962 , 37425.70438 , 133876.59532 ) T into Equation (2), we can obtain
t e ( 0 ) ( t ) = ( 1 e a ) ( t e ( 0 ) ( 1 ) b a + b a 2 c a ) e a ( t 1 ) + b a
Then
t e ( 0 ) ( t ) = 108129.384713 e 0.14962 × ( t 1 ) + 250141.126576
Equation (17) is just the HGEM(1,1) model for forecasting the energy consumption of China’s manufacturing.
Step 3 Model performance comparisons and tests
When t = 1 , 2 , , 6 , we can simulate the energy consumption of China’s manufacturing; when t = 7 , the energy consumption of China’s manufacturing can be forecast. In order to compare the simulation and prediction performance of the HGEM(1,1), we also apply the classic GM(1,1) model and the frequently-used DGM(1,1) model to simulate and forecast the energy consumption of China’s manufacturing industry. All the simulation and prediction results of the above three models are shown in Table 2, as follows:
Table 2 shows the HGEM(1,1) model to be better than the DGM(1,1) and GM(1,1) models in both simulation and prediction accuracy. In order to clearly compare the performance of the three models, we drew the simulation and prediction curves and errors of the three models based on the data in Table 2 in MATLAB, as shown in Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6, as follows.
According to Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6, the performances of the simulation and prediction of the HGEM(1,1) model are the best among the above three models, and those of the DGM(1,1) model and the GM(1,1) model are very close to each other. By checking the grey model error level reference table, we can see that the comprehensive grade of the HGEM(1,1) model is I, which can be used for prediction.
Step 4 Prediction
According to Equation (17), the total energy consumption of China’s manufacturing industry in the next 2018–2024 years can be predicted, as shown in Table 3.
From Table 3, we can see that the total energy consumption of China’s manufacturing industry will slow down during 2018–2024, but the total amount is still too large. In order to ensure the effective supply of China’s manufacturing energy, and promote the sustainable development of China’s economy, some countermeasures and suggestions are put forward in Section 5.
From Table 2, the out-of-sample data is only one, and it is difficult to test the prediction performance of the HGEM(1,1) model usefully. To this end, we adjusted the proportion of in-sample and out-of-sample and built a new HGEM(1,1) model, HGEM(1,1)II to simulate and forecast the energy consumption of China’s manufacturing. The parameters of HGEM(1,1)II are shown in Table 4, as follows.
We built the HGEM(1,1)II model with the parameters in Table 4, and the simulation and prediction results of the HGEM(1,1)II model are shown in Table 5, as follows:
It is obvious from Table 2 and Table 5 that the comprehensive mean relative percentage error (CMRPE) of the HGEM(1,1) model is only 0.494%, which is better than that of the other seven models. Therefore, it is reliable and reasonable that we apply the HGEM(1,1) model to forecast the energy consumption of China’s manufacturing (CMEC) during 2013–2024.

5. Suggestions

By 2024, the energy consumption of China’s manufacturing will be as high as 242824.144 ten thousand tons of standard coal, which accounts for approximately 60% of China’s total energy consumption in 2016 according to Table 3. Therefore, it is a major concern of China’s government that some measures need to be taken to control the energy consumption of China’s manufacturing industry. In order to control the high energy consumption of China’s manufacturing industry, the main measures are to optimize the industrial structure, speed up the elimination of backward production capacity, and promote the transformation and upgrading of traditional industries. After this, the green development of the manufacturing industry can be achieved.
First, it should adjust and optimize China’s manufacturing structure. The excessive growth trend of high energy consumption and pollution industries need be controlled. Some policies and measures to promote industrial restructuring should be improved and should actively promote the adjustment of the energy structure in China’s manufacturing industry. Meanwhile, we should formulate policies and measures to promote the development of high-tech industries. Secondly, we will eliminate backward production capacity and control the development of ‘two high’ industries which are highly polluting and high energy consuming enterprises. China’s government should eliminate backward technology, technology and equipment in small and medium-sized enterprises, especially high energy consumption and heavy pollution industries, such as iron and steel, nonferrous metals, chemicals, building materials, power and other industries. China’s government departments should eliminate “two high” and “five small” enterprises and the backward production capacity of enterprises be in accordance with the relevant requirements, step by step. Thirdly, they should speed up energy-saving emission reduction technology development and generalization, and promote the development of a circular economy. In view of the characteristics of small and medium-sized enterprises, we will speed up the development of common, key and cutting-edge energy-saving emission reduction technologies, and foster a technological innovation system of energy saving and emission reduction, which combines enterprises as the main body and combines production, production and research, and speeds up the transformation of scientific and technological achievements. A variety of methods are used to speed up the promotion and application of the efficient energy-saving new technologies, new processes, new products and new equipment.

6. Conclusions

China is the largest manufacturing energy consuming country in the world. Energy supply and price are the two key factors that influence the sustainable development of China’s manufacturing industry. Hence, a scientific prediction of the total energy consumption of China’s manufacturing industry has a positive significance on the smooth and healthy development of China’s manufacturing industry. To this end, the HGEM(1,1) model specially used for energy prediction was constructed, and we studied the parameter estimation method, time response formula and performance test method of the new model. Finally, the HGEM(1,1) model was applied to simulate and forecast the energy of China’s manufacturing industry, and corresponding countermeasures and suggestions were put forward according to the prediction results.
Optimizing the initial value, background value and accumulating order of the HGEM(1,1) model, and then constructing a better grey prediction model of total energy consumption are the next research target of our team.

Acknowledgments

Our work was supported by National Natural Science Foundation of China (71771033, 71271226), Chongqing Municipal Social Science Planning Commission Project of China (2016WT37), Science and technology research project of Chongqing Education Commission (KJ1706166) and the open research fund of Chongqing Key Laboratory of Electronic Commerce & Supply Chain System (1456026). We would like to thank the anonymous referees for their constructive comments that helped to improve the clarity and completeness of this paper.

Author Contributions

Bo Zeng conceived and designed the study. Meng Zhou performed the modeling process and error checking methods for the HGEM(1,1) model. Jun Zhang provided the MATLAB code with HGEM(1,1) and performed the calculations of Section 4. All authors read and approved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wu, J.; Mohamed, R.; Wang, Z.; Wu, J.; Mohamed, R.; Wang, Z. An agent-based model to project China’s energy consumption and carbon emission peaks at multiple levels. Sustainability 2017, 9, 893. [Google Scholar] [CrossRef]
  2. Liu, P.; Zhou, Y.; Zhou, D.K.; Xue, L. Energy performance contract models for the diffusion of green-manufacturing technologies in China: A stakeholder analysis from SMEs’ perspective. Energy Policy 2017, 106, 59–67. [Google Scholar] [CrossRef]
  3. Wan, Z. Analysis of economic growth and energy consumption of China’s manufacturing industry. East China Sci. Technol. Acad. Ed. 2015, 3, 421. (In Chinese) [Google Scholar]
  4. Sun, W.; Ji, S.S. Research on the Energy Consumption Index of Manufacturing Based on Maximum Deviation. Sci. Technol. Manag. Res. 2016, 354, 233–239. (In Chinese) [Google Scholar]
  5. Hou, J.; Chen, H.; Xu, J. External knowledge sourcing and green innovation growth with environmental and energy regulations: Evidence from manufacturing in China. Sustainability 2017, 9, 342. [Google Scholar] [CrossRef]
  6. Pan, H.; Zhang, H.; Zhang, X. China’s provincial industrial energy efficiency and its determinants. Math. Comput. Model. 2013, 58, 1032–1039. [Google Scholar]
  7. Li, K.; Lin, B. Impact of energy conservation policies on the green productivity in China’s manufacturing sector: Evidence from a three-stage dea model. Appl. Energy 2016, 168, 351–363. [Google Scholar] [CrossRef]
  8. Yang, M.; Yang, F. Energy-efficiency policies and energy productivity improvements: Evidence from China’s manufacturing industry. Emerg. Mark. Financ. Trade 2016, 52, 1395–1404. [Google Scholar] [CrossRef]
  9. Zhu, J.; Ruth, M. Relocation or reallocation: Impacts of differentiated energy saving regulation on manufacturing industries in China. Ecol. Econ. 2015, 110, 119–133. [Google Scholar] [CrossRef]
  10. Wang, Z.X.; Zheng, H.H.; Pei, L.L.; Jin, T. Decomposition of the factors influencing export fluctuation in China’s new energy industry based on a constant market share model. Energy Policy 2017, 109, 22–35. [Google Scholar] [CrossRef]
  11. Qu, Y.; Yu, Y.; Appolloni, A.; Li, M.R.; Liu, Y. Measuring Green Growth Efficiency for Chinese Manufacturing Industries. Sustainability 2017, 9, 637. [Google Scholar] [CrossRef]
  12. Chen, X.; Gong, Z.W. DEA Efficiency of Energy Consumption in China’s Manufacturing Sectors with Environmental Regulation Policy Constraints. Sustainability 2017, 9, 210. [Google Scholar] [CrossRef]
  13. Lin, B.; Liu, W. Scenario Prediction of Energy Consumption and CO2 Emissions in China’s Machinery Industry. Sustainability 2017, 9, 87. [Google Scholar] [CrossRef]
  14. Liu, S.F.; Forrest, J.; Yang, Y.J. A brief introduction to grey systems theory. Grey Syst. Theory Appl. 2012, 2, 89–104. [Google Scholar] [CrossRef]
  15. Helaleh, A.H.; Alizadeh, M. Performance prediction model of miscible surfactant-CO2, displacement in porous media using support vector machine regression with parameters selected by ant colony optimization. J. Nat. Gas Sci. Eng. 2016, 30, 388–404. [Google Scholar] [CrossRef]
  16. Lakshmanan, G.T.; Shamsi, D.; Doganata, Y.N.; Unuvar, M.; Khalaf, R. A markov prediction model for data-driven semi-structured business processes. Knowl. Inf. Syst. 2015, 42, 97–126. [Google Scholar] [CrossRef]
  17. Chae, Y.T.; Horesh, R.; Hwang, Y.; Lee, Y.M. Artificial neural network model for forecasting sub-hourly electricity usage in commercial buildings. Energy Build. 2016, 111, 184–194. [Google Scholar] [CrossRef]
  18. Deng, J.L. The Control problem of grey systems. Syst. Control Lett. 1982, 1, 288–294. [Google Scholar]
  19. Liu, S.F.; Lin, Y. Grey Systems Theory and Applications; Springer-Verlag: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  20. Zeng, B.; Chen, G.; Liu, S.F. A novel interval grey prediction model considering uncertain information. J. Frankl. Inst. 2013, 350, 3400–3416. [Google Scholar] [CrossRef]
  21. Li, C.; Cabrera, D.; de Oliveira, J.V.; Sanchez, R.V.; Cerrada, M.; Zurita, G. Extracting repetitive transients for rotating machinery diagnosis using multiscale clustered grey infogram. Mech. Syst. Signal Process. 2016, 76–77, 157–173. [Google Scholar] [CrossRef]
  22. Wang, Z.X.; Li, Q.; Pei, L.L. Grey forecasting method of quarterly hydropower production in China based on a data grouping approach. Appl. Math. Model. 2017, 51, 302–316. [Google Scholar] [CrossRef]
  23. Wang, Y.H.; Dang, Y.G.; Li, Y.Q.; Liu, S.F. An approach to increase prediction precise of GM(1,1) model based on optimization of the initial condition. Expert Syst. Appl. 2010, 37, 5640–5644. [Google Scholar] [CrossRef]
  24. Zeng, B.; Li, C.; Long, X.J. Equivalency and unbiasedness of grey prediction models. J. Syst. Eng. Electron. 2015, 26, 110–118. [Google Scholar] [CrossRef]
  25. Zeng, B.; Meng, W.; Tong, M.Y. A self-adaptive intelligence grey predictive model with alterable structure and its application. Eng. Appl. Artif. Intell. 2016, 50, 236–244. [Google Scholar] [CrossRef]
  26. Xie, N.M.; Liu, S.F.; Yang, Y.J.; Yuan, C.Q. On novel grey forecasting model based on non-homogeneous index sequence. Appl. Math. Model. 2013, 37, 5059–5068. [Google Scholar] [CrossRef]
  27. Zeng, B.; Meng, W.; Liu, S.F. Research on prediction model of oscillatory sequence based on GM(1,1) and its application in electricity demand prediction. J. Grey Syst. 2013, 25, 31–40. [Google Scholar]
  28. Zeng, B. Forecasting the relation of supply and demand of natural gas in China during 2015–2020 using a novel grey model. J. Intell. Fuzzy Syst. 2017, 32, 141–155. [Google Scholar] [CrossRef]
Figure 1. The modelling flowchart of the HGEM(1,1) model.
Figure 1. The modelling flowchart of the HGEM(1,1) model.
Sustainability 09 01975 g001
Figure 2. The simulated/forecast curve of the HGEM(1,1) model.
Figure 2. The simulated/forecast curve of the HGEM(1,1) model.
Sustainability 09 01975 g002
Figure 3. The simulated/forecast curve of the DGM(1,1) model.
Figure 3. The simulated/forecast curve of the DGM(1,1) model.
Sustainability 09 01975 g003
Figure 4. The simulated/forecast curve of the GM(1,1) model.
Figure 4. The simulated/forecast curve of the GM(1,1) model.
Sustainability 09 01975 g004
Figure 5. The errors of simulated/forecast curve of HGEM(1,1), DGM(1,1) and GM(1,1).
Figure 5. The errors of simulated/forecast curve of HGEM(1,1), DGM(1,1) and GM(1,1).
Sustainability 09 01975 g005
Figure 6. The mean relative simulation/prediction percentage error of HGEM(1,1), DGM(1,1) and GM(1,1).
Figure 6. The mean relative simulation/prediction percentage error of HGEM(1,1), DGM(1,1) and GM(1,1).
Sustainability 09 01975 g006
Table 1. The energy consumption of China’s manufacturing (CMEC) during 2006–2012.
Table 1. The energy consumption of China’s manufacturing (CMEC) during 2006–2012.
Year2006200720082009201020112012
k 1234567
CMEC143,051.46156,218.80172,106.52180,595.96188,497.85200,403.37205,667.69
Unit: Ten thousand tons of standard coal.
Table 2. Simulated/forecast values and errors of HEGM(1,1), DGM(1,1) and GM(1,1).
Table 2. Simulated/forecast values and errors of HEGM(1,1), DGM(1,1) and GM(1,1).
Serial Number k Raw Data t e ( 0 ) ( k ) Model HGEM(1,1)Model DGM(1,1)Model GM(1,1)
t ^ e ( 0 ) ( k ) ε S ( k ) Δ S ( k ) t ^ e ( 0 ) ( k ) ε S ( k ) Δ S ( k ) t ^ e ( 0 ) ( k ) ε S ( k ) Δ S ( k )
1143,051.46143,051.4600.0000.000%143,051.460.0000.000%143,051.460.0000.000%
2156,218.80157,037.777818.9770.524%159,377.0783158.2782.022%159,317.0443098.2441.983%
3172,106.52169,975.7422130.7781.238%168,891.0703215.4501.868%168,831.3173275.2031.903%
4180,595.96181,115.803519.8430.288%178,972.9981622.9620.899%178,913.7741682.1860.931%
5188,497.85190,707.8022209.9521.172%189,656.7661158.9160.615%189,598.3461100.4960.584%
6200,403.37198,966.8631436.5070.717%200,978.299574.92900.287%200,920.991517.6210.258%
MRSPE ( Δ ¯ S )0.788%1.138%1.132%
In-Sample (Simulated Data)
k t e ( 0 ) ( k ) t ^ e ( 0 ) ( k ) ε F ( k ) Δ F ( k ) t ^ e ( 0 ) ( k ) ε F ( k ) Δ F ( k ) t ^ e ( 0 ) ( k ) ε F ( k ) Δ F ( k )
7205,667.69206,078.217410.5270.199%212,975.6697307.9793.553%212,919.8157252.1253.526%
MRFPE ( Δ ¯ F )0.199%4.572%4.570%
Out-of-Sample (Forecast Data)
CMRPE( Δ ¯ )0.494%2.855%2.851%
Table 3. Prediction data of the energy consumption of China’s manufacturing (CMEC) during 2013–2024.
Table 3. Prediction data of the energy consumption of China’s manufacturing (CMEC) during 2013–2024.
YearCMECYearCMECYearCMECYearCMEC
2013212,201.3522016225,921.9452019234,680.6042022240,271.769
2014217,473.5952017229,287.5272020236,829.0522023241,643.249
2015222,013.1892018232,185.4152021238,678.9442024242,824.144
Unit: Ten thousand tons of standard coal.
Table 4. Value of the parameters of the HGEM(1,1)II model.
Table 4. Value of the parameters of the HGEM(1,1)II model.
Parameter α ^ β ^ γ ^ a ^ b ^ c ^
Value0.66167968,008.99593068,748.1271360.66167968,008.99593068,748.127136
Table 5. Simulated/forecast values and errors of HEGM(1,1)II, DGM(1,1) and GM(1,1).
Table 5. Simulated/forecast values and errors of HEGM(1,1)II, DGM(1,1) and GM(1,1).
Serial Number k Raw Data t e ( 0 ) ( k ) Model HGEM(1,1) IIModel GM(1,1) IIModel DGM(1,1) II
t ^ e ( 0 ) ( k ) ε S ( k ) Δ S ( k ) t ^ e ( 0 ) ( k ) ε S ( k ) Δ S ( k ) t ^ e ( 0 ) ( k ) ε S ( k ) Δ S ( k )
1143,051.46143,051.4600.0000.000%143,051.460.0000.000%143,051.460.0000.000%
2156,218.80156,368.819150.0190.096%158,990.838−2772.0381.774%159,060.6952841.8951.819
3172,106.52171,474.973−631.5470.367%168,791.7173314.8031.926%168,855.452−3251.0681.889
4180,595.96181,470.400874.440.484%179,196.7641399.1960.775%179,253.359−1342.6010.743
5188,497.85188,497.850150.0190.096%190,243.223−1745.3730.926%190,291.5561793.7060.952
MRSPE ( Δ ¯ S )0.237%1.350%1.351%
In-Sample (Simulated Data)
k t e ( 0 ) ( k ) t ^ e ( 0 ) ( k ) ε F ( k ) Δ F ( k ) t ^ e ( 0 ) ( k ) ε F ( k ) Δ F ( k ) t ^ e ( 0 ) ( k ) ε F ( k ) Δ F ( k )
6200,403.37192,460.354−7943.0163.964%201,970.6331567.2630.782%202,009.4711606.1010.801%
7205,667.69195,355.988−10,311.7025.014%214,420.9718753.2814.256%214,448.9618781.2714.270%
MRFPE ( Δ ¯ F )4.489%2.519%2.536%
Out-of-Sample (Forecast Data)
CMRPE( Δ ¯ )2.564%1.935%1.944%

Share and Cite

MDPI and ACS Style

Zeng, B.; Zhou, M.; Zhang, J. Forecasting the Energy Consumption of China’s Manufacturing Using a Homologous Grey Prediction Model. Sustainability 2017, 9, 1975. https://0-doi-org.brum.beds.ac.uk/10.3390/su9111975

AMA Style

Zeng B, Zhou M, Zhang J. Forecasting the Energy Consumption of China’s Manufacturing Using a Homologous Grey Prediction Model. Sustainability. 2017; 9(11):1975. https://0-doi-org.brum.beds.ac.uk/10.3390/su9111975

Chicago/Turabian Style

Zeng, Bo, Meng Zhou, and Jun Zhang. 2017. "Forecasting the Energy Consumption of China’s Manufacturing Using a Homologous Grey Prediction Model" Sustainability 9, no. 11: 1975. https://0-doi-org.brum.beds.ac.uk/10.3390/su9111975

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