# Structural Analysis and Total Coal Demand Forecast in China

Authors: Qing Zhu; Zhongyu Zhang; Rongyao Li; Kin Keung Lai; Shouyang Wang; Jian Chai
Journal: Discrete Dynamics in Nature and Society
 (2014)

---

## Abstract

Considering the speedy growth of industrialization and urbanization in China and the continued rise of coal consumption, this paper identifies factors that have impacted coal consumption in 1985–2011. After extracting the core factors, the Bayesian vector autoregressive forecast model is constructed, with variables that include coal consumption, the gross value of industrial output, and the downstream industry output (cement, crude steel, and thermal power). The impulse response function and variance decomposition are applied to portray the dynamic correlations between coal consumption and economic variables. Then for analyzing structural changes of coal consumption, the exponential smoothing model is also established, based on division of seven sectors. The results show that the structure of coal consumption underwent significant changes during the past 30 years. Consumption of both household sector and transport, storage, and post sectors continues to decline; consumption of wholesale and retail trade and hotels and catering services sectors presents a fluctuating and improving trend; and consumption of industry sector is still high. The gross value of industrial output and the downstream industry output have been promoting coal consumption growth for a long time. In 2015 and 2020, total coal demand is expected to reach 2746.27 and 4041.68 million tons of standard coal in China.


---

## Body

## 1. Introduction

As the basic source of energy, coal has facilitated the rapid development of the Chinese economy and has had a positive effect on the stability of the market economy. Therefore, in order to avoid indiscriminate production effectively, improving the efficiency of the coal industry and ensuring national energy security and accurate prediction of coal demand are necessary.In the past, scholars have established a variety of energy demand forecasting models for different countries and regions as well as different kinds of energy. Traditional methods such as time series, regression, econometric, decomposition, unit root test and cointegration, ARIMA (autoregressive integrated moving average), and input-output, as well as soft computing techniques such as fuzzy logic, GA (genetic algorithm), and ANN (artificial neural networks) are being extensively used for demand side management. Support vector regression, ACO (ant colony), and PSO (particle swarm optimization) are new techniques being adopted for energy demand forecasting. Bottom up models such as MARKAL (acronym for MARKet ALlocation) and LEAP (long-range energy alternatives planning model) are also being used at the national and regional level for energy demand management [1].For foreign energy consumption, Ediger and Akar [2] considered that estimated economic and demographic parameters usually deviate from the realizations. Time-series forecasting appears to give better results, so they used the ARIMA and SARIMA (seasonal ARIMA) methods to estimate the future primary energy demand of Turkey from 2005 to 2020. Mohr and Evans [3] considered the model of worldwide coal production developed for three scenarios (Hubbert Linearisation method scenario, reserves plus cumulative production scenario, best guess scenario). The ultimately recoverable resources (URR) estimates used in the scenarios ranged from 700 Gt to 1243 Gt. Unler [4] proposed a model to forecast the energy demand of Turkey until 2025 by using PSO-based energy demand forecasting (PSOEDF) and made a comparison with the ant colony optimization (ACO) energy demand estimation model. The result showed that the former algorithm had better accuracy. Lee and Chiu [5] applied a newly developed panel smooth transition regression model with the error-correction term (PSECM) to estimate the nonlinear relationships among energy consumption, real income, and real energy prices for 24 OECD countries.For domestic energy consumption, various approaches to forecasting have been used in the literature, including cointegration and error correction models [6] and demand equations [7]. In 2002, Yu and Zhu [8] proposed an improved hybrid algorithm called PSO-GA (particle swarm optimization-genetic algorithm) for energy demand forecasting in China with higher precision compared with single optimization methods, such as GA, PSO or ant colony optimization, and multiple linear regressions. Results of this study show that China’s energy demand will be 4.70 billion tons coal equivalent in 2015. In the same year, they also proposed mix-encoding particle swarm optimization and radial basis function (MPSO-RBF) network-based model to forecast China’s energy consumption until 2020, based on GDP, population, proportion of industry in GDP, urbanization rate, and share of coal energy for the period from 1980 to 2009 [9]. Zhang et al. [10] forecast transport energy demand for 2010, 2015, and 2020 based on partial least square regression (PLSR) method under two scenarios by analyzing gross domestic product (GDP), urbanization rate, passenger-turnover, and freight-turnover for the period of 1990–2006. Crompton and Wu [11] applied the Bayesian vector autoregressive methodology to forecast China’s energy consumption and to discuss potential implications. The results of this paper suggested that total energy consumption should increase to 2173 MtCE in 2010, an annual growth rate of 3.8%, which is slightly slower than the average rate in the past decade due to structural changes in the Chinese economy.Among them, some of the latest techniques such as Bayesian vector autoregression (BVAR), support vector regression, ant colony, and particle swarm optimization models are being used in energy demand analysis. Chai et al. [12] established a VARX (vector autoregressive model with exogenous variables) model of crude oil market structure to study the effect of every variable on oil price by screening variables with price, supply, demand of oil, dollar index and China oil net import as the endogenous variables and reserve, speculative factors as the exogenous variables. Then, based on this VARX model and the Bayes theory, a MSBVAR (Markov switching bayes vector autoregression) model is established to identify and analyze the structural changing of oil price system structure within the study period. BVAR model and Granger-causality are applied to study growth in energy demand and the relationship between energy consumption and real gross domestic product per capita in selected few Caribbean countries [13]. Bayesian neural network approach is used for short term electric load forecasting [14, 15]. The BVAR model can avoid the rigid inclusion/exclusion restrictions of VAR models by allowing inclusion of many coefficients while simultaneously controlling the extent of mathematical expectation and standard deviation to which they can be influenced by the data. This reduces the spurious correlations captured by the model, thereby reducing forecast error and improving forecasting performance.On the basis of the background and demand mentioned above, this paper attempts to analyze structural changes of coal consumption during the past 30 years and forecast future coal demand in China up to 2020. This paper includes seven parts: introduction, describing the realistic background and the academic background for research topic; variables selection and data processing, screening the core effect factors for total coal demand; methodology, introducing VAR, BVAR, and ETS models; forecast of total coal consumption, based on VAR models, applying the impulse response function and variance decomposition to portray the dynamic correlations between coal consumption and economic variables, and forecasting total coal consumption by BVAR models; structural analysis of total coal consumption, establishing ETS models based on coal consumption of seven sectors to analyze structural changes of coal consumption and forecast total coal consumption; results and discussion, comparing forecast results of total coal consumption by BVAR models and ETS models; conclusions.


## 2. Variables Selection and Data Processing

Many factors that affect future demand for coal, such as domestic and international price of coal, the national policy of saving energy and reducing consumption, efficiency of resource utilization, coal consumption of downstream industries, railway and highway capacity for coal transportation, technology development, urbanization, industrial growth, industrial structure, the supply and price of other alternative energy sources (such as oil, natural gas, hydropower, nuclear power, etc.), coal production cost (pollution control costs, coal mine safety production costs, exit cost, etc.), and consumer habits.Based upon the China Statistical Yearbook, China Energy Statistical Yearbook, China Coal Industry Yearbook, and WIND Database, this paper has identified the following influence factors of coal demand from 1985 to 2011:(1) unit of coal consumption is ten thousand tons of standard coal; (2) since most coal is consumed in industry and contribution of the secondary industry to gross domestic product (GDP) fluctuates sharply, gross industrial output (100 million yuan) and industrial structure (%) are considered important factors that affect coal consumption; (3) electric power, steel, building materials, and chemical industry are four major downstream industries that consume coal, so output of thermal power (100 million kilowatt hour), crude steel (10000 tons), and cement (10000 tons) are also factors used to forecast coal consumption; (4) effect of international prices and supply position of different substitute sources of energy on domestic coal usage and consumption, such as coal price (Newcastle in Australia/Ken blah FOB), oil price (Europe Brent Spot Price FOB), natural gas price (Louisiana spot prices), and producer price indices for domestic industrial products by sector (power, coal, and petroleum); (5) per capita GDP (yuan) and urbanization rate (%) data are collected since the improvement of consumer capacity and change of urban structure affect the energy consumption per capita in China.First of all, in order to eliminate the heteroscedasticity of the economic time series data and to linearize its trend, the above variables are put into the natural logarithms. Then, to prevent false regression leading to an invalid conclusion, ADF (Augment Dikey-Fuller) unit root test and cointegration test are used to examine stationarity and long-run equilibrium relationship of logarithmic time series, respectively. Finally, only the following 5 groups’ data meet the conditions of the VAR model: coal consumption (COAL), gross value of industrial output (IND), output of cement industrial products (CEM), output of crude steel industrial products (STE), and output of thermal power industrial products (POW).


## 3. Methodology

### 3.1. VAR

A vector autoregressive (VAR) model based on statistical properties of the data is established, without exerting a theory of a priori constraints for the data mechanism, so it is an unstructured time series model, not a structural econometric model with a theoretical foundation. Specifically, all economic variables are considered as endogenous variables in VAR model, through multistage lag regression, to estimate their relationships and for establishing forecast models.In 1980, Sims [16] introduced the VAR model and promoted its application for dynamic analysis of the economic system extensively. VAR model is used to forecast interconnected time series system and analyze dynamic shocks from stochastic disturbances to variables, so as to explain the influence of various economic shocks on economic variables.An unrestricted vector autoregression (UVAR) model containingn time series variables and a lag length of p has the following general form:
(1)yt=c+A1yt-1+⋯+Apyt-p+μt,μt~i.i.d.N(0,σ2),t=1,2,…,n,
where yt=(y1,y2,…,yn)′ is an (n×1) vector of endogenous variables; c is an (n×1) vector of intercept terms; Ai,  i>1, are (n×n) coefficient matrices; and μt is an (n×1) vector of the independent normal random error vector.Based on the VAR model, the dynamic relationship between the variables can be analyzed by impulse response function and variance decomposition.Impulse response function describes the effect on the current and future variables from information shocks of a random disturbance which equals a standard deviation, and vividly depicts the path changes of the dynamic interaction between the variables and tests the intensity and duration of the impact of economic variables on coal consumption. This function implies that decomposition does not have to completely depend on the order of the variables in the VAR system, so as to improve the stability and reliability of the estimation results.Variance decomposition method decomposes mean square error (MSE) of the forecast system into contributions of each variable, facilitating examination of MSE decomposition of any endogenous variables and measurement of the relative importance of random disturbances for the variables by comparing the variance contribution rate. In this paper, through the variance decomposition method, the role of the economic variables fluctuation in coal consumption growth can be determined.


### 3.2. BVAR

In (1), the vector c contains n intercept terms and each matrix Ai contains n2 coefficients; hence, with the increase of variables and lag, n+pn2 coefficients must be estimated, which increase exponentially with the number of variables in the system. A major problem in estimation of VAR models when p is large is over-parameterization, where too many coefficients must be estimated, relative to sample size. This can lead to an overfitting problem: the large number of coefficients in unrestricted VAR models tends to fit the data unrealistically well, while performing poorly in out-of-sample forecasting due to the effects of spurious correlations in the data set. Over-fitting of the data can distort the long-run relationships between variables in the model and inflate coefficient values on distant lags due to low degrees of freedom [17]. In addition, VAR model ignores the a priori information and presets the same importance for all estimated parameters, which can lead to the wrong model. There are many ways to overcome the over-fitting problem. By applying certain constraints on the parameters, such as reducing the lag length, or removing some of the variables in an individual equation, can solve the problem of lower degrees of freedom effectively. But from the view of Bayesian, this means that forecasters believe that the probability of removing lags whose coefficient is zero is 100%. But it is impossible to know whether this constraint is established.An alternate approach is the BVAR method proposed by Litterman [18, 19], Doan et al. [20], Sims [21], and Sims and Zha [22]. The BVAR approach modifies the OLS estimates of (1) by treating all coefficients as random variables around their Bayesian prior mean, such that the model has the flexibility to impose these priors, to varying degrees, on the coefficient estimates.


### 3.3. ETS

ETS technology appeared and was used in the 1950s. After decades of development, it became mature. However, selection of “optimum” model from different models for use began few years ago. Hyndman and Khandakar [23] summed up exponential smoothing models and divided them into fifteen types. For more details, see Table 1.Table 1
Classification and summary of ETS models.
Trend component
Seasonal component
N (none)
A (additive)
M (multiplicative)
N (none)
N, N
N, A
N, M
A (additive)
A, N
A, A
A, M
Ad (additive damped)
Ad, N
Ad, A
Ad, M
M (multiplicative)
M, N
M, A
M, M
Md (multiplicative damped)
Md, N
Md, A
Md, MError, Trend, and Seasonal of ETS model, respectively, represent error term, trend term, and season term, in which, (Trend, Seasonal) group includes fifteen models listed in Table1. Residual error term is classified as overlapping forms (addition forms) and multiplication forms. On the assumption of yt=μt+εt, it is an additive errors model. On the assumption of yt=μt(1+εt), it is a multiplication errors model. Under the circumstances of consideration of different error forms, the above fifteen models can be extended to thirty types. From the forecast results alone, addition form or multiplication form hardly influence residual error term. However, for different sampled data, various residual error forms seem superior or inferior. Because dividing by 0 is involved, attention should be paid to select different addition or multiplication forms for residual error, trend term, or season term. When samples are all positive values, using residual error multiplication form presents an advantage. However, when samples are zero value or negative value, multiplication form model is not applicable.


## 3.1. VAR

A vector autoregressive (VAR) model based on statistical properties of the data is established, without exerting a theory of a priori constraints for the data mechanism, so it is an unstructured time series model, not a structural econometric model with a theoretical foundation. Specifically, all economic variables are considered as endogenous variables in VAR model, through multistage lag regression, to estimate their relationships and for establishing forecast models.In 1980, Sims [16] introduced the VAR model and promoted its application for dynamic analysis of the economic system extensively. VAR model is used to forecast interconnected time series system and analyze dynamic shocks from stochastic disturbances to variables, so as to explain the influence of various economic shocks on economic variables.An unrestricted vector autoregression (UVAR) model containingn time series variables and a lag length of p has the following general form:
(1)yt=c+A1yt-1+⋯+Apyt-p+μt,μt~i.i.d.N(0,σ2),t=1,2,…,n,
where yt=(y1,y2,…,yn)′ is an (n×1) vector of endogenous variables; c is an (n×1) vector of intercept terms; Ai,  i>1, are (n×n) coefficient matrices; and μt is an (n×1) vector of the independent normal random error vector.Based on the VAR model, the dynamic relationship between the variables can be analyzed by impulse response function and variance decomposition.Impulse response function describes the effect on the current and future variables from information shocks of a random disturbance which equals a standard deviation, and vividly depicts the path changes of the dynamic interaction between the variables and tests the intensity and duration of the impact of economic variables on coal consumption. This function implies that decomposition does not have to completely depend on the order of the variables in the VAR system, so as to improve the stability and reliability of the estimation results.Variance decomposition method decomposes mean square error (MSE) of the forecast system into contributions of each variable, facilitating examination of MSE decomposition of any endogenous variables and measurement of the relative importance of random disturbances for the variables by comparing the variance contribution rate. In this paper, through the variance decomposition method, the role of the economic variables fluctuation in coal consumption growth can be determined.


## 3.2. BVAR

In (1), the vector c contains n intercept terms and each matrix Ai contains n2 coefficients; hence, with the increase of variables and lag, n+pn2 coefficients must be estimated, which increase exponentially with the number of variables in the system. A major problem in estimation of VAR models when p is large is over-parameterization, where too many coefficients must be estimated, relative to sample size. This can lead to an overfitting problem: the large number of coefficients in unrestricted VAR models tends to fit the data unrealistically well, while performing poorly in out-of-sample forecasting due to the effects of spurious correlations in the data set. Over-fitting of the data can distort the long-run relationships between variables in the model and inflate coefficient values on distant lags due to low degrees of freedom [17]. In addition, VAR model ignores the a priori information and presets the same importance for all estimated parameters, which can lead to the wrong model. There are many ways to overcome the over-fitting problem. By applying certain constraints on the parameters, such as reducing the lag length, or removing some of the variables in an individual equation, can solve the problem of lower degrees of freedom effectively. But from the view of Bayesian, this means that forecasters believe that the probability of removing lags whose coefficient is zero is 100%. But it is impossible to know whether this constraint is established.An alternate approach is the BVAR method proposed by Litterman [18, 19], Doan et al. [20], Sims [21], and Sims and Zha [22]. The BVAR approach modifies the OLS estimates of (1) by treating all coefficients as random variables around their Bayesian prior mean, such that the model has the flexibility to impose these priors, to varying degrees, on the coefficient estimates.


## 3.3. ETS

ETS technology appeared and was used in the 1950s. After decades of development, it became mature. However, selection of “optimum” model from different models for use began few years ago. Hyndman and Khandakar [23] summed up exponential smoothing models and divided them into fifteen types. For more details, see Table 1.Table 1
Classification and summary of ETS models.
Trend component
Seasonal component
N (none)
A (additive)
M (multiplicative)
N (none)
N, N
N, A
N, M
A (additive)
A, N
A, A
A, M
Ad (additive damped)
Ad, N
Ad, A
Ad, M
M (multiplicative)
M, N
M, A
M, M
Md (multiplicative damped)
Md, N
Md, A
Md, MError, Trend, and Seasonal of ETS model, respectively, represent error term, trend term, and season term, in which, (Trend, Seasonal) group includes fifteen models listed in Table1. Residual error term is classified as overlapping forms (addition forms) and multiplication forms. On the assumption of yt=μt+εt, it is an additive errors model. On the assumption of yt=μt(1+εt), it is a multiplication errors model. Under the circumstances of consideration of different error forms, the above fifteen models can be extended to thirty types. From the forecast results alone, addition form or multiplication form hardly influence residual error term. However, for different sampled data, various residual error forms seem superior or inferior. Because dividing by 0 is involved, attention should be paid to select different addition or multiplication forms for residual error, trend term, or season term. When samples are all positive values, using residual error multiplication form presents an advantage. However, when samples are zero value or negative value, multiplication form model is not applicable.


## 4. Forecast of Total Coal Consumption by BVAR Models

In order to analyze the dynamic relationship between coal consumption and economic variables, 5 d VAR model is built based on (1). According to the lag length criteria (likelihood ratio, final prediction error, AIC, SC, and Hannan-Quinn information criterion), VAR (2) is the most appropriate model. VAR models are set up as follows:
(2)[LNCOALLNINDLNCEMLNSTELNPOW]=[4.4768-20.03133.1506-2.45372.3387]+[0.67850.1126-0.11740.05340.49130.13050.55861.4616-0.48480.6649-0.68340.04701.07560.01980.62380.4434-0.23510.06060.89570.0600-0.25710.01390.15780.34820.9504]×[LNCOAL(-1)LNIND(-1)LNCEM(-1)LNSTE(-1)LNPOW(-1)]+[-0.3453-0.1248-0.00670.1699-0.20822.0973-0.4018-0.3793-0.5128-0.48530.3309-0.0134-0.7911-0.04540.3111-0.29740.0022-0.0554-0.26830.7020-0.0600-0.0030-0.1852-0.19370.0594]×[LNCOAL(-2)LNIND(-2)LNCEM(-2)LNSTE(-2)LNPOW(-2)]+[ε1ε2ε3ε4ε5].TheR-squared of the models are 0.9981, 0.9951, 0.9966, and 0.9986, which indicate good simulation results. At the same time, the reciprocals of all roots locate inside the unit circle when calculating the AR characteristic polynomial (Figure 1), which shows that the established VAR (2) model is stable. That is to say, when a variable changes in the model (i.e., to generate a shock), it leads to changes in other variables too. But as time goes on, the effect gradually disappears. In Figure 1, the horizontal axis represents the real number axis; the vertical axis represents the imaginary numbers axis; the unit circle is a circle of radius 1 and center at the origin; the points in the unit circle represent roots of the characteristic equation on VAR (2) model.Figure 1
Inverse roots of AR characteristic polynomial.Based on the VAR model, the dynamic relationship between the variables can be analyzed by impulse response function and variance decomposition.In Figure2, the horizontal axis represents the lag periods of shock from innovation (unit: years); the vertical axis represents the response degree from the dependent variables to the explaining variables; the solid line is the calculated value of the impulse response function. The lag period is set to 15 years.Figure 2
Responses of LNCOAL to LNCOAL, LNIND, LNCEM, LNSTE, and LNPOW.According to Figure2, the line graph shows that when the explaining variables (coal consumption, gross value of industrial output, cement output, crude steel industrial output, and thermal power output) are shocked from innovation, respectively, how does the dependent variable (coal consumption) respond during the future 15 periods. That is to say, when these five factors change, respectively, we can forecast how much and how long will coal consumption respond.Firstly, coal consumption volatility responds strongly to its own shock which is shown by the blue line, and in the first 15 periods it shows a fluctuating trend. Early, as a strong positive response, in the second period, it reaches the peak at 0.0310% and then plunges to the lowest point at −0.0145% in the eighth period. From this point onwards, the period between 9 and 15 experiences a gradual growth until the numbers become positive response again. This shows that coal consumption and its lagging values have a strong correlation. Secondly, after suffering a positive shock, gross value of industrial output only brings positive response to coal consumption in the second period, and then it quickly drops to a negative response, to the lowest point of −0.0379% in the sixth period. Afterwards, it rises to zero slightly. This indicates that the increase of the gross value of industrial output immediately leads to an increase in coal consumption, but in the long term it reduces consumption. The reasons may be that more output provides more capital for internal restructuring and technological progress, thus reducing coal consumption costs and increasing profits. Thirdly, coal consumption volatility responds modestly to cement output shocks in the early stage. From the fifth period to the eleventh period, the response shows a large fluctuation from negative to positive. Fourthly, coal consumption volatility has increasing positive response to crude steel industrial output shocks gradually over time, and in the seventh period that reaches the peak of 0.0219%. Fifthly, a positive shock on thermal power output brings positive response to coal consumption, which has weak intensity and stable tendency. These results indicate that the development of downstream industries drives an increase in demand for coal, whose long-term relation is in the same direction.In Figure3, the horizontal axis represents the lag periods of impact from innovation (unit: years); the vertical axis represents the contribution rate to coal consumption from the economic variables.Figure 3
Factors variance decomposition on the effect of LNCOAL.As shown in Figure3, in the future 15 periods, coal consumption has the greatest impact during the first five periods with its own contribution rate above 40%. Since then, rates in the rest of the period level off at about 22%. The contribution rate of the gross value of industrial output to coal consumption increases steeply from the third period and then exceeds all other variables at the sixth period and remains stable and high between 51 and 53%. Thermal power industry had the highest contribution in the second period, but after that its contribution rates decline to 3-4%. Contribution rates of crude steel industry have obvious rise from the third period and make the largest contribution in the downstream industries between 12–14%. The contribution rate of cement industry has an obvious rise in the sixth period, after that, the contribution rate remains in the 7-8% level. Therefore, in the long run, industrial output has the biggest contribution to the changes of coal consumption; in a short time, coal consumption have the biggest contribution to own changes.The system of priors commonly used in the specification of BVAR models includes conjugate prior distribution, maximum entropy prior distribution, ML-II prior distribution, and multilayer prior distribution. Different BVAR prior distributions have different impacts on the prediction results. This paper uses theR software to forecast coal consumption based on three kinds of prior distribution in the MSBVAR package. Command is run as follows: R > fcast < -szbvar (coal, p=2, lambda 0=0.6, lambda 1=0.1, lambda 3=2, lambda 4=0.25, lambda 5=0, mu 5=0, mu 6=0, prior = 0). There are three different prior distributions: 0 = Normal-Wishart prior, 1 = Normal-flat prior, and 2 = flat-flat prior (i.e., akin to MLE). The forecasting results from 2012 to 2020 are shown in Table 4.


## 5. Structural Analysis of Total Coal Consumption

A univariate ETS (exponential smoothing) forecast model is used to predict coal consumption of each sector in 2020 in the forecast package forR. Error, Trend, and Seasonal of ETS model denote error term (A, M, Z), trend term (N, A, M, Z), and season term (N, A, M, Z), respectively, in which the group includes thirty models.From the coal balance sheet in China Statistical Yearbook and WIND Database, coal consumption data of seven sectors are selected for the period 1980–2010. The balance equation for coal consumption: Total coal consumption (TO) = Agriculture, Forestry, Animal Husbandry, Fishery and Water Conservancy (AFAFW) + Industry (IND) + Construction (CON) + Transport, Storage and Post (TSP) + Wholesale and Retail Trades, Hotels and Catering Services (WRHC) + Other Sectors (OS) + Household Consumption (HC). Unit is 10000 tons.The final models (Table2) are selected from thirty models by comparing information criterion (AIC, BIC, AICc) and forecast precision (mean error, root mean square error, mean absolute error, maximum permissible errors, mean absolute percentage error, mean absolute error square) (Table 3). Command is run as follows: R>fcast<-forecast(est(x)).Table 2
Formulae for recursive calculations, point forecasts, and parametric estimated value on sector by ETS method.
Variables
Forms
Models and parameters
AFAFW
(M, N, N)
l
t
=
α
y
t
+
(
1
-
α
)
l
t
-
1, y^t+h∣t=lt, yt+1=y^t+1∣t(1+εt), εt~NID(0,σ2), 0<α<1,0<β<α, α=0.0001,l=1793.7859,σ=0.1809
IND
(A, A, N)
l
t
=
α
y
t
+
(
1
-
α
)
(
l
t
-
1
+
b
t
-
1
), bt=β(lt-lt-1)+(1-β)bt-1, y^t+h∣t=lt+hbt yt+1=y^t+1∣t+εt, εt~NID(0,σ2), 0<α<1, 0<β<α, α=0.9999,β=0.2697, l=54782.1109,b=9338.5498,σ=12910.3600
CON
(A, A, N)
l
t
=
α
y
t
+
(
1
-
α
)
(
l
t
-
1
+
b
t
-
1
), bt=β(lt-lt-1)+(1-β)bt-1, y^t+h∣t=lt+hbt yt+1=y^t+1∣t+εt, εt~NID(0,σ2), 0<α<1,0<β<α, α=0.0009,β=0.0009,l=448.3868,b=8.7382,σ=53.6530
TSP
(M, M, N)
l
t
=
a
y
t
+
(
1
-
a
)
l
t
-
1
b
t
-
1, bt=β(lt/lt-1)+(1-β)bt-1, y^t+h|t=ltbth yt+1=y^t+1∣t(1+εt), εt~NID(0,σ2), 0<α<1, 0<β<α, α=0.0001,β=0.0001,l=2305.0531,b=0.9371,σ=0.1049
WRHC
(M, M, N)
l
t
=
a
y
t
+
(
1
-
a
)
l
t
-
1
b
t
-
1, bt=β(lt/lt-1)+(1-β)bt-1, y^t+h∣t=ltbth yt+1=y^t+1∣t(1+εt), εt~NID(0,σ2), 0<α<1, 0<β<α, α=0.0001,β=0.0001,l=765.1436,b=1.0328,σ=0.2562
OS
(A, N, N)
l
t
=
a
y
t
+
(
1
-
a
)
l
t
-
1, y^t+h∣t=lt yt+1=y^t+1∣t+εt, εt~NID(0,σ2), 0<α<1, 0<β<α, α=0.8917,l=1148.7978,σ=497.8256
HC
(M, N, N)
l
t
=
a
y
t
+
(
1
-
a
)
l
t
-
1, y^t+h∣t=lt, yt+1=y^t+1∣t(1+εt), εt~NID(0,σ2), 0<α<1,0<β<α, α=0.9999,l=11373.6307,σ=0.1327
TO
(A, A, N)
l
t
=
α
y
t
+
(
1
-
α
)
(
l
t
-
1
+
b
t
-
1
), bt=β(lt-lt-1)+(1-β)bt-1, y^t+h∣t=lt+hbt yt+1=y^t+1∣t+εt, εt~NID(0,σ2), 0<α<1, 0<β<α, α=0.9999,β=0.1665,l=78421.9911,b=9365.1674,σ=14646.2100
In the above formulae for recursive calculations and point forecasts:lt denotes the series level at time t; bt denotes the slope at time t; α and β are smoothness index; h denotes the lag period.Table 3
Accuracy of forecast models on sector by ETS method.
AFAFW
IND
CON
TSP
WRHC
OS
HC
TO
AIC
310.7869
469.4980
239.2015
272.4609
308.5274
328.7655
371.5486
474.7963
AICc
311.4536
471.9980
241.7015
274.9609
311.0274
329.4321
372.2153
477.2963
BIC
312.8760
473.6761
243.3796
276.6390
312.7055
330.8545
373.6377
478.9744
ME
−0.1636
1297.1479
0.2783
5.3547
−21.8557
45.5007
−105.4612
1903.6901
RMSE
324.4691
12910.3640
53.6530
154.4569
309.8943
497.8256
1540.4774
14646.2092
MAE
255.6755
10635.7804
39.6210
115.8266
272.7594
365.0402
1081.5817
11811.2234
MPE
−4.0490
−0.2898
−1.0556
−0.8693
−10.5543
−6.4377
−1.9142
−0.1609
MAPE
15.7232
9.8982
7.7865
8.7576
27.1766
30.0246
9.6963
9.3149
MASE
0.7907
0.7024
0.8702
0.7315
1.1268
0.9915
0.9609
0.7435Table 4
Forecasting result of coal consumption on sector by ETS method.
ETS
TO
AFAFW
IND
CON
TSP
WRHC
OS
HC
2011
328820.36
1793.79
313117.60
640.68
552.60
1555.68
2000.86
9159.17
2012
345931.25
1793.79
330203.50
649.42
517.87
1606.66
2000.86
9159.17
2013
363045.90
1793.79
347289.30
658.16
485.32
1659.31
2000.86
9159.17
2014
380164.41
1793.79
364375.20
666.91
454.81
1713.68
2000.86
9159.17
2015
397286.63
1793.79
381461.10
675.65
426.23
1769.84
2000.86
9159.17
2016
414412.38
1793.79
398546.90
684.39
399.44
1827.84
2000.86
9159.17
2017
431541.82
1793.79
415632.80
693.14
374.33
1887.74
2000.86
9159.17
2018
448674.80
1793.79
432718.70
701.88
350.80
1949.61
2000.86
9159.17
2019
465811.18
1793.79
449804.50
710.62
328.75
2013.50
2000.86
9159.17
2020
482951.14
1793.79
466890.40
719.37
308.09
2079.48
2000.86
9159.17
Unit: 10000 tons.As can be seen from the Table4, the total coal consumption is the sum of coal consumption in seven sectors, which is 3.9728663 billion tons in 2015. This value is similar to the goal in the development of coal industry in the twelfth five-year plan which is 3.9 billion tons in 2015. By 2020, China’s total coal consumption will reach 4.83 billion tons with an annual growth rate of 4.36%, of which industry and household consumption account for the largest proportion. In Figure 4, the industry coal consumption was 438 million tons, 811 million tons, 1.278 billion tons, and 2.96 billion tons in 1980, 1990, 2000, and 2010, respectively, and it is predicted to reach 4.669 billion tons in 2020 accounting for 96.67% of total coal consumption. Therefore, to reduce coal consumption in China, industrial energy conservation and emissions reduction are the fastest and most effective breakthrough points, including the adjustment of industrial structure and optimization of technology.Figure 4
Changes of industry and household coal consumption.However, as the second largest main channel, the proportion of household consumption is decreasing year by year. Figure4 shows that household consumption was up to 167 million tons in 1990, which accounts for 15.83% of the total coal consumption. By 2020, household consumption will be 0.92 million tons accounting for 1.89% of the total coal consumption. In the field of household energy consumption, as the diversity of the type of energy and development of technology, electricity and heating oil are more available and convenient than coal for residents, so less and less coal is used in the household rapidly. This is beneficial to environmental protection and nonrenewable resources saving.As shown in Figure5, transport, storage, and post sectors were the largest coal consumption sectors in 1980, but they became the lowest coal consumption sectors after thirty years in 2010, and their coal consumption will decline constantly during the next decade by predictions. The main reason is that coal-fired steam locomotive engines have been replaced by internal combustion and electric locomotives by the railways. So, most of the raw coal has been replaced by fuel and electricity at the same time. Nowadays, road transportation mainly consumes gasoline and diesel oil; waterways transportation mainly consumes fuel oil and diesel oil; air transport mainly consumes aviation kerosene. Therefore, although the development of transport, storage, and post sectors still consumes large amounts of energy, the usage of coal is insignificant.Figure 5
Changes of other sectors’ coal consumption.As the pillar industry of national economy in our country, construction sector is also an energy-intensive industry. Coal consumption increased from 5.56 to 7.1892 million tons between 1980 and 2010, representing a relatively small growth. By forecasting, the amount of coal consumption is still increasing in the future (Figure5). At the same time, the construction sector does not have a high proportion of coal consumption and accounted for only 0.23% in 2010. That is because there are various types of energies consumed by the construction industry, including raw coal, gasoline, diesel oil, fuel oil, heat, electricity, and other petroleum products. With the development of technology, the construction sector will become dependent on power and all kinds of oil products increasingly instead of coal, which can improve the energy efficiency and control environmental pollution gradually.Wholesale and retail trade and hotels and catering services consume coal, electricity, heat, gasoline, diesel oil, liquefied petroleum gas, natural gas, and other varieties of energy. Influenced by price changes and economic development, this sector has gone through booms and busts repeatedly. From 1980 to 1990, the sector had a rising trend due to adapting to the market demand quickly; from 1990 to 1997 with high price, the sector was in the integration stage; from 1997 to 1999, the sector was in a low price and low growth stage; from 2000 to 2004, the sector was in a low price and high growth stage; domestic prices began to soar in 2005 and since the American subprime crisis and the domestic policy adjustment in 2008, the sector appears to have had violent fluctuations. However, Figure5 shows that the wholesale and retail trade, hotels, and catering services will have a growth trend of coal consumption in the future. Industry scale will continue to expand, which is consistent with our country’s policy to develop the tertiary industry vigorously.Agriculture, forestry, animal husbandry, fishery, and water conservancy are the foundation of the national economy. Coal consumption was in the second place before 1990, but in 1993, 2000, 2005, and 2008 there was a great fall (Figure5); in contrast, coal consumption of wholesale and retail trade, hotels, and catering services, and other sectors had an obvious increase in these three years. Considerable contact with the change of the prices is guessed.


## 6. Results and Discussion

Table5 shows prediction results of BVAR models under the three different prior distributions and the ETS models. The unit is 10000 tons in ETS models, which is 10000 tons of standard coal in BVAR models. So the results of ETS model are multiplied by the converted coefficient of 0.7143 as tons of standard coal.Table 5
Total coal consumption forecast.
BVAR Prior = 0 or 1
BVAR Prior = 2
ETS aggregation
ETS forecast
2011
238033.37
238033.37
234876
234475
2012
252686.75
254034.07
247099
245921
2013
268152.68
257667.41
259324
257366
2014
284494.14
257113.39
271551
268811
2015
301762.9
274627.29
283782
280256
2016
320015.13
314720.65
296014
291701
2017
339311.96
360700.12
308250
303146
2018
359720.03
401190.03
320489
314591
2019
381312.12
444913.39
332729
326037
2020
404167.84
504823.69
344972
337481
Unit: 10000 tons of standard coal.The result of BVAR models under the normal-Wishart prior distribution (prior = 0) is the same as under the normal-flat prior distribution (prior = 1). Using the ETS model to predict coal consumption of every sector and then adding the seven values to get 3.44972 billion tons of standard coal in 2020; only based on the historical total coal consumption, 3.37481 billion tons of standard coal are obtained in 2020 directly. From Table5, in 2015 forecast results of four methods have small differences. But in 2020, the results of BVAR model are larger than ETS model (Figure 6).Figure 6
Total coal consumption forecast.The reason may be that in BVAR model, the gross value of industrial output and the downstream industrial production are influencing factors, whose values are increasing over the years. So they drive the rise in total coal consumption in the long run. However, in ETS model, total coal consumption is the sum of coal consumption in seven sectors, and the above has introduced that the coal consumption of some sectors has showed a decreasing trend, such as agriculture, forestry, animal husbandry, fishery and water conservancy, transport, storage, and post. So the ETS model prediction results are smaller.


## 7. Conclusions

China is a big developing country in the process of industrialization and urbanization, and therefore coal consumption will remain high in the future, for a long time. Based on guaranteeing the stability of the VAR system, coal consumption, gross industrial output value, and the downstream industrial production (crude steel, cement, thermal power) are selected as variables to establish the BVAR (2) model. At the same time, this paper analyses the changes of coal consumption structure based on sector division. The forecast models are selected through comparing information criteria and forecast accuracy comparison from thirty ETS models. By forecasting the total coal consumption and analyzing structural changes, this paper gets the following main conclusions.(1)
Results of the impulse response function indicate that the increased output of the downstream industry will become a driving force of coal consumption. In the first two periods, the same response has been there, with the increase of the gross value of industrial output. Later there will be a negative response.(2)
The variance decomposition results show that the contribution of gross value of industrial output accounted for 50% of growth in coal consumption and has long-term stable influence; in the downstream industry, the contribution rate of crude steel output is greater than that of the cement output and thermal power output.(3)
In the next decade, from a sector division point of view, industry and household consumption account for the largest proportion of total coal consumption; the proportion of industrial coal consumption will continue to rise, but the proportion of household coal consumption will continue to decline. Transport, storage, and post sectors were the largest consumers of coal consumption in 1980, which will become the lowest in the future. Construction sector, despite its lower reliance on coal, has experienced a rise in coal consumption in recent years. The scales of consumption of wholesale and retail trade, hotels and catering services, agriculture, forestry, animal husbandry, fishery, and water conservancy sectors were influenced by price fluctuations in the past few decades, causing the fluctuation of coal consumption. But in the next decade, coal consumption of both sectors will have a rising trend.(4)
The BVAR (2) forecast models predict that in 2015 coal demand will reach 2.75–3.02 billion tons of standard coal, and in 2020 it will reach 4.04–5.05 billion tons of standard coal in China. That is to say, as the momentum of rapid development of China’s economy in the future, the coal consumption in 2020 is almost double that of in 2012 which is 2.46–2.53 billion tons of standard coal. However, because the unrenewable characteristics of the coal and the pollution caused by burning are bound to draw the attention of the country and the residents, how to improve the coal utilization efficiency and promote the use of clean energy will become a main research direction in the future.



---

*Source: 612064-2014-06-05.xml*