Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Forecasting China's energy intensity by using an improved DVCGM (1, N) model considering the hysteresis effect

Forecasting China's energy intensity by using an improved DVCGM (1, N) model considering the... Purpose – The purpose of this paper is to examine the effectiveness of an improved dummy variables control grey model (DVCGM) considering the hysteresis effect of government policies in China’s energy intensity (EI) forecasting. Design/methodology/approach – Energy consumption is considered as an important driver of economic development. China has introduced policies those aim at the optimization of energy structure and EI. In this study, EI is forecasted by an improved DVCGM, considering the hysteresis effect of energy-saving policies of the government. A nonlinear optimization method based on particle swarm optimization (PSO) algorithm is constructed to calculate the hysteresis parameter. A one-step rolling mechanism is applied to provide input data of the prediction model. Grey model (GM) (1, N), DVCGM (1, N) and ARIMA model are applied to test the accuracy of the improved DVCGM (1, N) model prediction. Findings – The results show that the improved DVCGM provides reliable results and works well in simulation and predictions using multivariable data in small sample size and time-lag virtual variable. Accordingly, the improved DVCGM notes the hysteresis effect of government policies and significantly improves the prediction accuracy of China’s EI than the other three models. Originality/value – This study estimates the EI considering the hysteresis effect of energy-saving policies in China by using an improved DVCGM. The main contribution of this paper is to propose a model to estimate EI, considering the hysteresis effect of energy-saving policies and improve forecasting accuracy. Keywords Grey system theory, Energy intensity forecasting, Grey prediction model, Hysteresis effect Paper type Research paper © Zhaosu Meng, Xiaotong Liu, Kedong Yin, Xuemei Li and Xinchang Guo. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) licence. Anyone may reproduce, distribute, translate and create derivative works of this article (for both Grey Systems: Theory and Application commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this licence may be seen at http://creativecommons.org/licences/by/4.0/ pp. 372-393 Emerald Publishing Limited legalcode. 2043-9377 DOI 10.1108/GS-02-2020-0022 This work was supported by the National Social Science Foundation of China (Grant No. 19BJY151). 1. Introduction Hysteresis Energy consumption is considered as an important driver of economic growth. Energy effect of intensity (EI), which is the ratio of total energy consumption to gross domestic product (GDP) China’sEI of a country or region in a certain period of time, measures the performance of energy forecasting utilization. According to the energy section of the National Bureau of Statistics (NBS), China’s largest energy consumption is 4,640 million tons of standard coal, and China has become the world largest energy consumer. Entering a new normal period, the government focuses on a crucial rebalancing and diversifying economy, with higher requirements for sustainable development. In the 13th Five-Year Plan for economic and social development, a total reduction of 15% has been set as the energy performance target. It is imperative to estimate and predict the EI to evaluate energy conservation and emission reduction. However, due to the complexity and dynamics of social economy, the implementation of energy conservation and emission reduction policies will not immediately reduce EI, and there is often a certain time lag. This inevitable time-lag should be taken into account to estimate and forecast the EI in China accurately. Numerous research studies have studied the influencing factors and driving forces of EI through econometrics methods such as cointegration analysis, metrology and decomposition analysis and scenario analysis. Zhu et al. (2015) applied cointegration analysis method on the large and state-owned enterprises and found that energy-saving regulations in China are one of the most important factors in reducing aggregate EI. Karimu et al. (2017) studied the EI and convergence of Swedish industry by combining metrology and decomposition analytical methods. Ma and Yu (2017) used panel data model to discuss the driving factors which lead to EI decline. It is revealed that industrial structure, energy conservation regulations and EI are closely related. Tan et al. (2018) used index decomposition analysis and production decomposition analysis methods to analyze the factors which are related to the decline in the EI and pointed out that technology improvement effect is the most significant factor. In EI forecasting, Pao et al. (2012) used improved grey models (GMs) to predict China’sCO emissions, energy consumption and economic growth. Dong et al. (2018) estimated the driving force of regional EI in China and forecasted the potential of regional energy conservation with scenario analysis. Wu et al. (2018) used a new multivariable GM to predict energy consumption in Shandong Province. The grey system theory is an interdisciplinary theory proposed by Deng (1982) and the grey prediction method applies well to small size data forecasting. As an important part of grey prediction theory, the GM (1, N) model is the basic model of multivariable grey system modeling approach. In recent years, numerous scholars have thoroughly discussed the model parameters optimization (Tien, 2005, 2010), the model accuracy improvement (Tien, 2011; Wang et al., 2016) and the expansion of the GM (Guo et al., 2013; Kose and Tasci, 2015, 2019; Ding et al., 2017). Based on grey multivariable model with time-lagged system, Zhai et al. (1996) introduced the lag term into the GM (1, 2) model and determined the delay parameters with the goal of minimizing the modeling error. Hao (2011) used grey correlation analysis to determine time-lag period between variables and then on this basis to establish forecast model of GM (1, N). Zhang et al. (2015) constructed a time-delay multivariable discrete GM,DDGM (1, N) model, by introducing a time-delay control factor and solved the time-delay parameters by using the grey dimensionally expanding identification method, which obtained a good application effect. Ma and Yu (2017) used a novel time-delay multivariable GM to predict the natural gas consumption in China. Dang et al. (2017) constructed the discrete delay grey multivariable DDGMD (1, N) model by introducing the driving information control adjustment coefficient T and the action coefficient λ and solved the coefficient respectively by grey dimension expansion method and particle swarm optimization (PSO) algorithm. Xiong (2019) built a multivariable time-delay discrete MGM (1, m, t) model and studied the mechanism of modeling and the process of modeling, and the GS calculation method of time delay is given. The hysteresis effect is discussed by an example 11,3 to verify the validity of the model. GM (1, N) model, in spite of successfully applying in various fields, sometimes ignores the influence of virtual variables such as policy on the main system in practical applications. Zhang (2016) considered the influence of dummy variables on system behavior variables and built a discrete multivariable prediction model based on dummy variables, which further expanded the application scope of the model. Ding et al. (2018) introduced the dummy variables into the GM (1, N) model, gave the concrete model construction method in mathematics and verified the effectiveness of the new model with cases. Generally speaking, we have abundant literature discussing the influencing factors of EI and its forecasting. Furthermore, there are some practices targeting the time-lag phenomenon using grey theory. The existing literature has explored and studied the GM from multiple angles, but there is still room for the GM to expand in the combination of dummy variables and time-delay systems. As the common GM (1, N) model does not take into account the time delay between variables, and there are dummy variables in the system which are difficult to be measured by quantity, our optimized DVCGM (1, N) model, which is short for dummy variables control grey model of N variables, incorporates the above two features and improves prediction accuracy to a satisfying level. It is practical significant because the time delay of dummy variables, i.e. government policy, is often seen in real world but is difficult to measure. The optimized DVCGM (1, N) model takes these variables into consideration, solves the practical problem and expands the grey system theory and the grey prediction method system, which improves the accuracy of grey prediction model. This paper’s main interest is to estimate and forecast EI by considering the influence of government policies and to test the accuracy of the improved GM through a comparison study. The main contribution of the study to the literature is to consider hysteresis effect and increase the forecast accuracy. It is imperative for optimizing energy structure, improving energy utilization efficiency and ensuring energy security. The results show that the improved GM produces better results than the other three conventional models. The rest of the paper is organized as follows. Section 2 briefly introduces grey theory and multivariable grey prediction models. A nonlinear optimization method based on PSO algorithm is constructed to calculate the hysteresis parameter in the improved model. In Section 3,EIis estimated by considering the hysteresis effect of energy-saving policies by an improved DVCGM (1, N) model. We compare the results with the other two GMs and one econometric model. The improved DVCGM (1, N) model has the best performance in the estimation comparison and is applied to forecast future EI. Section 4 is the conclusion of the study with the limitation and future path. 2. Methodology The three grey prediction models used in this paper are interrelated and in a progressive order. The GM (1, N) model is a traditional multivariable grey prediction model. The DVCGM (1, N) model introduces virtual variables, taking into account the influence of policy and other factors on the basis of the GM (1, N) model. Time-delay parameter is introduced in the improved DVCGM (1, N) model, considering the hysteresis effects of historical variables, which further enriches the existing grey prediction theory. 2.1 GM (1, N) model Grey prediction model can be regarded as two levels of work. At a lower level, the original sequence produces the sequence of generation by one accumulative generation (1-AGO), and then it forms the sequence of mean generation of consecutive neighbors; similarly, the Hysteresis sequence of influencing factors generates the sequence of generation by 1-AGO; Constructing effect of B and Y matrix and calculating system parameters through ordinary least squares (OLS) China’sEI regression. Once the system parameters a and b are determined, we can obtain the time forecasting response function (TRF) by solving the whitenization equation. At a higher level, the continuous differential equation with initial values is used as the reflection equation, and the discrete data are mapped to a manageable function, which is further restored to the TRF as the simulation basis. Typical procedures can be described briefly by the program in Figure 1. Definition 1. The original sequence is ð0Þ ð0Þ ð0Þ ð0Þ X ¼ x ð1Þ; x ð2Þ; :::; x ðnÞ 1 1 1 1 ð1Þ ð0Þ ð1Þ X is the 1 - AGO sequence from X , and X is the sequences of the relevant factors, where 1 1 i ð1Þ ð0Þ X ðkÞ¼ x ðmÞ i i m¼1 ð1Þ ð1Þ Z is the sequence of mean generation of consecutive neighbors from X , where 1 1 ð1Þ ð1Þ ð1Þ Z ðkÞ¼ 0:5X ðk  1Þþ 0:5X ðkÞ; ðk ¼ 2; 3; ... ; nÞ 1 1 1 Definition 2. Denote Eq (1) as the definition formula of GM (1, N) model: ð0Þ ð1Þ ð1Þ X ðkÞþ aZ ðkÞ¼ b X ðkÞ (1) 1 1 i i¼2 where k 5 2, 3, n, a is the development coefficient and b ; b ;  ; b are the grey input 1 2 n coefficients obtained by the least squares method. To determine these coefficients, the matrix B and Y are defined as follows: 2 3 2 3 ð0Þ ð1Þ ð1Þ ð1Þ z ð2Þ X ð2Þ X ð2Þ x ð2Þ 1 1 2 n 6 7 6 7 ð0Þ 6 ð1Þ ð1Þ 7 ð1Þ 6 7 x ð3Þ z ð3Þ X ð3Þ X ð3Þ 6 7 1 1 2 n 6 7 6 7 Y ¼ 6 7; B ¼ 6 7 . . . . 6 7 . . . . 6 7 4 . 5 . . . 4 5 ð0Þ ð1Þ ð1Þ ð1Þ x ðnÞ z ðnÞ X ðnÞ X ðmÞ 1 1 2 n The values of the coefficients a and b ; b ;  ; b can be determined by the following equation: 1 2 n 2 3 6 b 7 T −1 T b 6 7 P ¼ ¼ðB BÞ B Y 4 5 Definition 3. Let Eq (2) be defined as the differential equation (or) called grey reflection equation: ð0Þ n dx ð0Þ ð1Þ þ ax ðtÞ¼ b x ðtÞ (2) dt i¼1 GS 11,3 Figure 1. Principal diagram of GM (1, N) model Initialization Immediate original Generating mean sequence sequence1 sequence Grey equation Influencing Generating factor Time response Whitenization sequence2 function sequence Whitenization equation system BMatrix OLS parameters a and Y Matrix and b Prediction After determining the coefficients of a and b ; b ;  ; b , the differential equation of the GM 1 2 n Hysteresis can be determined by Eq (2). The solution of the above differential equation is as follows: effect of ( ) n n X X ð1Þ China’sEI ð0Þ ð1Þ b ð1Þ −ak b b b X ðk þ 1Þ¼ X ð1=aÞ b X ðkÞ 3e þð1=aÞ b X ðkÞ (3) i−1 i−1 1 1 i i forecasting i¼2 i¼2 ð1Þ where X ðkÞ is the prediction of the AGO of the original sequence. By considering that the estimation of the first element of the first AGO of a sequence is equal to the first element of the ð1Þ ð0Þ sequence, the following relation is determined: X ð1Þ¼ X ð1Þ 1 1 Finally, in order to predict the elements of the original sequence, the inverse accumulated generating operation should be performed. Therefore, the predicted values can be determined as follows: ð0Þ ð1Þ ð1Þ b x ðkÞ¼ b x ðkÞ b x ðk  1Þ; k≥ 2; 1 1 1 ð0Þ where x ðnÞ is an estimation of the original sequence, which is simulation values, ð0Þ ð0Þ b x ðn þ 1Þ;b x ðn þ 2Þ;  are predictive values. 1 1 2.2 DVCGM (1, N) model Traditional GM (1, N) model ignores the influence of virtual variables, which will inevitably lead to significant errors in practical applications. Therefore, it is necessary to construct a new multivariable predictive model with virtual variable control, based on the traditional GM (1, N) model, i.e. the DVCGM (1, N) model. The modeling steps for DVCGM (1, N) can be illustrated in Figure 2. Definition 4. The original sequence is ð0Þ ð0Þ ð0Þ ð0Þ X ¼ x ð1Þ; x ð2Þ;  ; x ðnÞ : i i i i Virtual variable sequence is ð0Þ ð0Þ ð0Þ ð0Þ ð0Þ D ¼ d ð1Þ; d ð2Þ; d ðnÞ ;  ; d ðnÞ¼ 0or 1: j j j j j ð1Þ ð1Þ ð0Þ X and D are the 1-AGO sequence, X is the behavior sequence of the system, i j 1 ð1Þ ð1Þ X ði ¼ 2;  ; MÞ and D ðj ¼ M þ 1;  ; NÞ is the driving factor sequence. Then GM i j (1, N) model with dummy variable can be expressed as: M N X X ð0Þ ð1Þ ð1Þ ð1Þ x ðkÞþ az ðkÞ¼ b x ðkÞþ b d ðkÞ (4) i j 1 1 i j i¼2 j¼Mþ1 named DVCGM (1, N) model (dummy variables control grey model of N variables). ð1Þ ð1Þ Z is the sequence of mean generation of consecutive neighbors from X where, 1 1 ð1Þ ð1Þ ð1Þ Z ðkÞ¼ 0:5X ðk− 1Þþ 0:5X ðkÞ; ðk ¼ 2; 3; ... ; nÞ 1 1 1 M N P P ð1Þ ð1Þ where b x ðkÞ is independent quantization variable driver, b d ðkÞ is the dummy i j i j i¼2 j¼Mþ1 variable driver. If virtual variables are not considered or b 50, the model can be transformed to a traditional GM (1, N) model. GS 11,3 Figure 2. Principal diagram of DVCGM (1, N) mode Initialization (0) (1) (1) 1 1 Z (0) (1) Grey equation X (0) (1) Time response D D Whitenization j j function TT 1 OLS B Y Matrix PB =()B BY Whitenization n equation Prediction Theorem 1. Assuming Hysteresis ð0Þ ð1Þ ð0Þ ð1Þ effect of X ; X ; D ðkÞ; D ðkÞ i i j j China’sEI as mentioned in Definition 4. The parameter column of the model is forecasting b ¼½a; b ;  b  : 2 N The matrix B and Y are defined as follows: 2 3 2 3 ð1Þ ð1Þ ð1Þ ð1Þ ð1Þ ð0Þ z ð2Þ X ð2Þ X ð2Þ d ð2Þ d ð2Þ x ð2Þ 1 2 1 M Mþ1 N 6 7 6 7 6 ð1Þ ð1Þ ð1Þ ð1Þ ð1Þ 7 ð0Þ 6 7 x ð3Þ z ð3Þ X ð3Þ X ð3Þ d ð3Þ d ð3Þ 6 7 1 1 2 M Mþ1 N 6 7 6 7 Y ¼ 6 7; B ¼ 6 7 . . . . . . 6 7 6 . . . . . 7 4 . 5 . . . . . 4 5 ð0Þ ð1Þ ð1Þ ð1Þ ð1Þ ð1Þ x ðnÞ z ðnÞ X ðnÞ X ðnÞ d ðnÞ d ðnÞ 1 1 2 M Mþ1 N Definition 5. Assuming ð0Þ ð1Þ ð0Þ ð1Þ X ; X ; D ðkÞ; D ðkÞ i i j j As mentioned in Definition 4, the parameter column of the model is b ¼½a; b ;  ; b  : 2 N Let Eq (5) be defined as the differential equation of DVCGM (1, N) model: ð1Þ M N X X dx ð1Þ ð1Þ ð1Þ þ ax ðtÞ¼ b x ðtÞþ b d ðtÞ (5) i j 1 i j dt i¼2 j¼Mþ1 After determining the coefficients of a and b ; b ;  ; b , the differential equation of the GM 1 2 n can be determined by Eq (5). The solution of the above differential equation is as follows: ( " # M N X X ð1Þ ð0Þ ð1Þ ð1Þ −at x ðtÞ¼ e X  t b x ð0Þþ b d ð0Þ i j 1 1 i j i¼2 j¼Mþ1 " # ) N M N X X X ð1Þ ð1Þ at þ b x ðtÞþ b d ðtÞ e dt (6) i j i j i¼2 i¼2 j¼Mþ1 Finally, in order to predict the elements of the original sequence, the inverse accumulated generating operation should be performed. Therefore, the predicted values can be determined as follows: ð0Þ ð1Þ ð1Þ b b b x ðkÞ¼ x ðkÞ x ðk  1Þ; k≥ 2; 1 1 1 ð0Þ where b x ðnÞ is an estimation of the original sequence, which is simulation values, ð0Þ ð0Þ b b x ðn þ 1Þ; x ðn þ 2Þ;  are predictive values. 1 1 2.3 Improved DVCGM (1, N) model The classic multivariable grey prediction models, such as traditional GM (1, N) and DVCGM (1, N) models, can reflect the influences of current driving-variables on the present system behavior and innately ignore the hysteresis effect of historical variables. Therefore, an GS improved DVCGM (1, N) model is proposed, integrating these above prediction model. 11,3 2.3.1 Construction of the improved DVCGM (1, N) model. In this section, the hysteresis parameter λ is innovatively introduced into the DVCGM (1, N) model to improve the prediction accuracy. Supported by PSO algorithm, detailed process and algorithm can be described as following: Theorem 2. Assuming ð0Þ ð1Þ ð0Þ ð1Þ X ; X ; D ðkÞ; D ðkÞ i i j j As mentioned in Definition 4. The parameter column of the model is b ¼½a; b ;  ; b  : 2 N The matrix B and Y are defined as follows: 2 3 ð0Þ x ð2Þ 6 7 ð0Þ 6 7 x ð3Þ 6 7 Y ¼ 6 7; 6 7 4 5 ð0Þ x ðnÞ 2 3 2 2 2 2 X X X X ð1Þ 2−j ð1Þ 2−j ð1Þ 2−j ð1Þ 2−j ð1Þ Z ð2Þ; λ x ðjÞ;  λ x ðjÞ; λ d ðjÞ λ d ðjÞ 1 2 2 M M Mþ1 Mþ1 N N 6 7 j¼1 j¼1 j¼1 j¼1 6 7 6 7 3 3 3 3 6 7 X X X X 6 ð1Þ ð1Þ ð1Þ ð1Þ ð1Þ 7 3−j 3−j 3−j 3−j Z ð3Þ; λ x ðjÞ;  λ x ðjÞ; λ d ðjÞ λ d ðjÞ 6 7 1 2 2 M M Mþ1 Mþ1 N N 6 7 B ¼ j¼1 j¼1 j¼1 j¼1 6 7 6 7 6 . . . . . 7 . . . . . 6 7 . . . . . 6 7 n n n n 4 X X X X 5 ð1Þ ð1Þ ð1Þ ð1Þ ð1Þ n−j n−j n−j n−j Z ðnÞ; λ x ðjÞ; ... λ x ðjÞ; λ d ðjÞ ... λ d ðjÞ 1 2 2 M M Mþ1 Mþ1 N N j¼1 j¼1 j¼1 j¼1 The least square estimation of the parameter column satisfies the following requirements: −1 (1) When n 5 Nþ1, b ¼ B Y ; jBj≠ 0; −1 T T T (2) When n > Nþ1, b ¼ðB BÞ B Y ;B B≠ 0; −1 T T T (3) When n < Nþ1, b ¼ B ðB BÞ Y ;B B≠ 0; Proof: Substitute k 5 2, 3,..., n into the model, you can get M k N k X X X X ð1Þ ð1Þ ð1Þ k−j k−j ð1Þ x ð2Þ¼ −ax ð2Þþ b λ x ð2Þþ b λ d ð2Þ i q 1 1 i i q q i¼2 j¼1 q¼Mþ1 j¼1 M k N k X X X X ð1Þ ð1Þ k−j ð1Þ k−j ð1Þ x ð3Þ¼ −ax ð3Þþ b λ x ð3Þþ b λ d ð3Þ i q 1 1 i i q q i¼2 j¼1 q¼Mþ1 j¼1 M k N k X X X X ð1Þ ð1Þ ð1Þ k−j k−j ð1Þ x ðnÞ¼ −ax ðnÞþ b λ x ðnÞþ b λ d ðnÞ i q 1 1 i i q q i¼2 j¼1 q¼Mþ1 j¼1 That is, by the least square method, Y ¼ Bb (1) When n5 Nþ1, B has an inverse matrix, the equations have a unique solution, we can Hysteresis −1 get b ¼ B Y ; jBj≠ 0; effect of China’sEI (2) When n > Nþ1, B is column full rank, the full rank decomposition of B is B 5 DC. forecasting Then the generalized inverse matrix of B can be obtained: −1  −1 −1 −1 T T T T T T T T B ¼ C CC D D D ; β ¼ C CC D D D Y ; Because B is a full rank matrix, C can be taken as a unit matrix, B 5 D,so −1 −1 −1 T T T T T T b ¼ C CC D D D Y ¼ B B B Y ; (3) When n < Nþ1, B is a row full rank matrix, D can be taken as a unit matrix, B5 C,so −1 −1 −1 T T T T T T b ¼ C CC D D D Y ¼ B B B Y ; Definition 6. Let Eq (7) be defined as the differential equation of improved DVCGM (1, N) model: Z Z ð1Þ M T N Q X X dx ðtÞ ð1Þ ð1Þ 1 t−s t−s ð1Þ þ ax ðtÞ¼ b λ x ðsÞds þ b λ d ðsÞds (7) i Q 1 i q q dt 0 0 i¼1 q¼Mþ1 Where k 5 2, 3, n, a is the development coefficient and b ; b ;  ; b are the grey input 1 2 n coefficients, λ is the hysteresis parameter. PSO algorithm is used to determine the coefficients of a, b ; b ;  ; b , and λ . Then the differential equation of the GM can be determined by 1 2 n i Eq (7). The solution of Eq (7) is as follows: ( " # M N X X ð1Þ ð0Þ ð1Þ −at ð1Þ x ðtÞ¼ e X  t b x ð0Þþ b d ð0Þ i q 1 1 i q i¼2 q¼Mþ1 " # ) N M k N k X X X X X k−j ð1Þ k−j ð1Þ at þ b λ x ðtÞþ b λ d ðtÞ e dt (8) i q i i q q i¼2 i¼2 j¼1 q¼Mþ1 j¼1 When the range of the driving factor sequence is small, the driver term can be viewed as a grey constant, and then the approximate TRF sequence of the grey differential equation of the model is " # M k N k X X X X ð1Þ ð1Þ k−j k−j ð1Þ x ðkÞ¼ b λ x ðtÞþ b λ d ðtÞ i q 1 i i q q i¼2 j¼1 q¼Mþ1 j¼1 " # M k N k −aðk−1Þ X X X X k−j ð1Þ ð1Þ k−j ð1Þ −aðk−1Þ b λ x ðtÞþ b λ d ðtÞ þ X e (9) i q i i q q 1 i¼2 j¼1 q¼Mþ1 j¼1 Finally, in order to predict the elements of the original sequence, the inverse accumulated GS generating operation should be performed. The predicted values can be determined as 11,3 follows: ð0Þ ð1Þ ð1Þ b x ðkÞ¼ b x ðkÞ b x ðk  1Þ; k≥ 2; 1 1 1 ð0Þ where b x ðnÞ is an estimation of the original sequence, which is simulation values, ð0Þ ð0Þ b x ðn þ 1Þ;b x ðn þ 2Þ;  are predictive values. 1 1 2.3.2 Estimation and optimization of hysteresis parameter in the improved DVCGM (1, N) model. The most important part of the improved DVCGM (1, N) model is estimating the time-lag parameter, which directly affects the accuracy of the model. However, the time-lag parameters must be determined in advance, followed by B and Y matrix construction and system parameters calculation through OLS. Once the system parameters a and b are determined, we can obtain the TRF and the simulation andpredictionvalue of the model. In this paper, a nonlinear optimization model is established by using the Least One Multiplication. Then, the time-lag parameter is determined. When the range of the driving factor sequence is small, Eq (9) is used as the TRF, λ can be solved by the following nonlinear programming model: " # M k N k X X X X ð1Þ ð1Þ k−j k−j ð1Þ b x ðkÞ¼ b λ x ðtÞþ b λ d ðtÞ i q 1 i i q q i¼2 j¼1 q¼Mþ1 j¼1 min (10) " # M k N k −aðk−1Þ X X X X k¼2 k−j ð1Þ k−j ð1Þ ð1Þ −aðk−1Þ b λ x ðtÞþ b λ d ðtÞ þ X e i q i i q q 1 i¼2 j¼1 q¼Mþ1 j¼1 8 ð0Þ ð1Þ ð1Þ b x ðkÞ¼ b x ðkÞ b x ðk  1Þ; k ≥ 2 > 1 1 1 " # M k N k > X X X X > 1 ð1Þ ð1Þ > k−j k−j ð1Þ b x ðkÞ¼ b λ x ðtÞþ b λ d ðtÞ i q 1 i i q q > a > i¼2 j¼1 q¼Mþ1 j¼1 " # M k N k −aðk−1Þ s:t: X X X X (11) k−j ð1Þ k−j ð1Þ ð1Þ −aðk−1Þ >  b λ x ðtÞþ b λ d ðtÞ þ x ð1Þe i q > i i q q 1 > i¼2 j¼1 q¼Mþ1 j¼1 b ¼½a; b ;  ; b > 2 N 0 < λ < 1; i ¼ 2; 3; ... ; M; M þ 1; ... ; N The model takes the relationship between structural parameters as the constraint condition and minimizes the average simulation relative error of the system characteristic variables, which can improve accuracy to the greatest extent. The above optimization problem can be solved by PSO (Kiran and Mustafa, 2017; Mason et al., 2018; Chen et al., 2018). According to Eq (10), a nonlinear optimization method based on PSO algorithm is constructed to obtain the hysteresis parameter. PSO sets a certain number of particles in feasible region to find the best location and can be used to seek optimal values of λ .Denote λ in Eq (10) and construct the fitness function of each particle, i i according to Eq (12). " # M k N k X X X X  Hysteresis ð1Þ ð1Þ k−j k−j ð1Þ b x ðkÞ¼ b λ x ðtÞþ b λ d ðtÞ i q 1 i i q q  effect of i¼2 j¼1 q¼Mþ1 j¼1 China’sEI " # forecasting M k N k −aðk−1Þ X X X X Fitnes λ ¼ (12) i e k−j ð1Þ k−j ð1Þ b λ x ðtÞþ b λ d ðtÞ i q k¼2 i i q q i¼2 j¼1 q¼Mþ1 j¼1 ð1Þ −aðk−1Þ þX e Obviously, the average simulation relative error of the system characteristic variable sequence varies depending on the number of lag periods. The value of lag period should be selected to make the average simulation relative error as small as possible. Therefore, the improved DVCGM (1, N) model can well describe the hysteresis effect between the system characteristic variables. Once the hysteresis parameter is determined, the structural parameters of the model are set accordingly, and the simulation and prediction results can be obtained according to Eq (9). 2.3.3 Modeling procedure. Detailed procedure of the improved DVCGM (1, N) model is illustrated as follows. ð0Þ ð1Þ Step 1. Collect raw data and establish original sequence X , X is the 1-AGO sequence 1 1 ð0Þ ð1Þ ð1Þ from X , X is the sequence of the relevant factors. Z is the sequence of mean 1 i 1 ð1Þ generation of consecutive neighbors from X . Then, determining virtual variable sequence ð0Þ ð0Þ ð0Þ ð0Þ ð0Þ D ¼ d ð1Þ; d ð2Þ;  d ðnÞ ; d ðnÞ¼ 0or 1: j j j j j ð1Þ D is also the sequence of the relevant factors. Step 2. Solving delay parameters λ by PSO according (10), then constructing vector Y and matrix B. Using the least square method, the values of the coefficients a and b ; b ;  ; b 1 2 n can be determined. Step 3. After considering the lag effect, the differential equation of the grey model can be determined by Eq (7). Step 4. The time response function of the new model is established to generate prediction data according to Eq (9). 3. Application Forecasting EI can be considered as a grey problem, because EI is greatly affected by technological progress, population factor, industrial structure and so on. These factors influence EI through a dynamic and complicated mechanism. The uncertain impacts and limited number of data provides a good basis for grey theory application. There are four procedures in this part, including data collection, parameter estimation, result comparisons and future forecasts. Three competing models, namely GM (1, N), DVCGM (1, N) and ARIMA model, are employed to test the accuracy of the improved DVCGM (1, N) on EI forecasting. The three GMs used in this paper are interrelated and are from basic to the advanced. The GM (1, N) model is a traditional grey multivariate prediction model. On its basis, DVCGM (1, N) model introduces dummy variables, taking into account the influence of policies and other factors. The improved DVCGM (1, N) model introduces time-lag parameters, which further enriches the existing grey prediction theory by considering the hysteresis effect of GS policies. The above GMs are used to estimate China’s EI, and the optimal model is selected by 11,3 comparing their performance, to predict China’s EI in the next five years. In addition, we use ARIMA model as a comparison of grey methods to show that the improved DVCGM (1, N) model is not only better than the traditional GM, but also better than the non-grey econometric model. 3.1 Variables selection and data collection The indicators selected in this paper are as follows. EI is the ratio of total energy consumption to GDP. Population factor is the employed population. The ratio of the added value of industrial production to the added value of energy consumed by the industry is used as the substitution variable of technological progress (Yan, 2011). Industrial structure is measured by the ratio of output value of each industry to GDP. For the consistency of the statistical scope, we choose a time scale of 2001–2017, and all data are collected from China Statistical Yearbook. With small sample size (17 periods’ real measurement values) and insufficient information, this case fits well with the grey system. First, we calculate the grey correlation between EI and population, technological progress and industrial structure. As shown in Table 1, technological progress has the highest correlation with EI. Therefore, we select technological progress as the driving variable, and EI is the system behavior variable. Then, we can build GM (1, N) model for EI forecast. As shown in Figure 3, China’s EI increased rapidly from 2001, rising sharply to a peak of 1.4542 tons standard coal per 10,000 Yuan in 2005 and declined steadily afterward. In 2017, the EI has decreased by 62.67%, leaving only 0.5428 tons of standard coal per 10,000 Yuan. The 11th Five-Year Plan in 2006 is a kind of watershed for Chinese EI. The government introduced a series of strict policies of energy conservation and emission reduction, and Relative correlation Energy intensity Population factors 0.6126 Table 1. Technological progress 0.9550 Relative correlation Industrial structure 0.7652 analysis tons standard coal per 10,000 Yuan 1.6 0.39 0.37 1.4 0.35 0.33 1.2 0.31 0.29 0.27 0.25 0.8 Figure 3. 0.23 Trends of energy intensity (dash line) 0.21 0.6 and technological 0.19 progress (solid line) 0.4 0.17 of China 2001 2003 2005 2007 2009 2011 2013 2015 2017 Energy intensity Technological progress China’s EI declined obviously. As a result, we need to put policy and the hysteresis effect of Hysteresis policy into consideration when estimate EI of China. effect of China’sEI 3.2 Simulation of energy intensity in China forecasting Four models are applied in the simulation of EI of China. First, we use the GM (1, N) model mentioned in Section 2.1 to predict EI. We select technological progress as the driving variable, EI as the system behavior variable and then establish the GM (1, N) model. According to 2.1, the TRF is obtained as follows: ð1Þ ð1Þ −1:8419ðk−1Þ −1:8419ðk−1Þ X ðkÞ¼ 4:2429X ðkÞ3 1  e þ e 1 2 The simulation results of GM (1, N) model are illustrated in Table 3, and the relative errors and average relative errors can be calculated, as shown in Table 3. Secondly, we build DVCGM (1, N) model for EI estimate. We select technological progress as the driving variable, EI as the system behavior variable. Energy conservation policy (P)is introduced as a dummy variable. Before 2006, as the strict energy-saving policies had not been implemented, P value is 0; after 2006, P value is 1. According to Section 2.2, the TRF is obtained as follows: ð1Þ ð1Þ ð1Þ −1:7372ðk−1Þ −1:7372ðk−1Þ X ðkÞ¼ 4:2830X ðkÞ 0:0232d ðkÞ 3 1  e d þ e 1 2 3 The simulation results of DVCGM (1, N) model are illustrated in Table 3, and the relative errors and average relative errors can be calculated, as shown in Table 3. Thirdly, we use the improved DVCGM (1, N) model to estimate EI. As discussed in Section 2.3, DVCGM (1, N) model takes consideration of the hysteresis effect of dummy variables. We select technological progress as the driving variable, EI as the system behavior variable and energy-saving policy as virtual variable. Different from DVCGM (1, N) model, we need to determine the delay parameters, along with the structure parameters and build the TRF to get the simulation results. (1) Determination of time-delay parameter. As stated by Eq (10-12), a nonlinear optimization method based on PSO algorithm is constructed to obtain the hysteresis parameter. According to the optimization model shown in Eq (9), the average relative percentage errors (APE) of the model under different lag periods is calculated in Table 2. When the hysteresis parameter is 2, the average relative percentage errors (APE) of the model is the smallest (0.9754%). Therefore, the value of time lag is two years. In the view of the complexity of socio-economic ecology, adaptive adjustments are made according to the changes of policies. These energy-saving policies influence economic operation through a certain transmission mechanism and gradually affect the EI. (2) Determination of structural parameters. The hysteresis parameter is substituted into the B matrix and Y matrix. According to Theorem 2, the matrix operation of least square regression is used to obtain the structural parameter values b , b . B 5 (1.251669, 0.938994). 1 2 Table 2. Average percentage Lag period 0 1 2 3 4 5 errors (APEs) of the model under different APE(%) 1.1107 1.0023 0.9754 1.0667 1.0989 0.9836 lag periods GS Improved 11,3 DVCGM (1, N) DVCGM (1, N) Years GM (1, N) model model model ARIMA model In- Original Simulation PE Simulation PE Simulation PE Simulation PE sample data values % values % values % values % 2001 1.3731 1.3731 0 1.3731 0 1.3731 0 1.3731 0 2002 1.3932 1.4108 1.26 1.40783 1.05 1.4195 1.89 1.4381 3.22 2003 1.4341 1.5118 5.42 1.50494 4.94 1.4866 3.66 1.5134 5.53 2004 1.4229 1.4697 3.29 1.45264 2.09 1.4502 1.92 1.4536 2.16 2005 1.4542 1.4955 2.84 1.48736 2.28 1.4847 2.1 1.5111 3.91 2006 1.3054 1.3319 2.03 1.34535 3.06 1.3271 1.66 1.361 4.26 2007 1.1525 1.1965 3.82 1.19042 3.29 1.1681 1.35 1.235 7.16 2008 1.0034 1.0611 5.75 1.05487 5.13 1.044 4.05 1.0605 5.69 2009 0.9629 1.0025 4.11 0.99564 3.4 0.9729 1.04 1.0488 8.92 2010 0.8732 0.9179 5.12 0.91127 4.36 0.8875 1.64 0.9616 10.1 2011 0.791 0.82 3.66 0.81457 2.98 0.804 1.64 0.8661 9.49 2012 0.7442 0.8106 8.92 0.77739 4.46 0.7518 1.02 0.7913 6.33 2013 0.7004 0.7648 9.19 0.76 8.51 0.7112 1.54 0.7575 8.15 2014 0.6612 0.678 2.54 0.6852 3.63 0.6683 1.08 0.6966 5.36 APE% 4.14 3.51 1.76 5.74 Table 3. Simulation of energy Out- Original Prediction PE Prediction PE Prediction PE Prediction PE intensity in China by sample data value % value % value % value % GM (1, N), DVCGM (1, N), the Improved 2015 0.6239 0.6577 5.41 0.64118 2.77 0.6334 1.52 0.6718 7.68 DVCGM (1, N) and 2016 0.5861 0.6076 3.67 0.61136 4.31 0.5943 1.4 0.6157 5.05 ARIMA model (Unit: 2017 0.5428 0.5676 4.57 0.56033 3.23 0.5582 2.83 0.5677 4.59 tons of standard coal per 10,000 Yuan) APE% 4.55 3.44 1.92 5.78 (3) Calculate the simulated and predicted values. Putting the values of the estimated structural parameter and hysteresis parameter in Eq (7-9), we can obtain the optimal TRF as follows: ð1Þ ð1Þ ð1Þ −1:7372ðk−1Þ −1:7372ðk−1Þ X ðkÞ¼ 4:2830X ðkÞ 0:0232d ðkÞ 3 1  e d þ e 1 2 3 Then, the simulation results of improved DVCGM (1, N) model are illustrated in Table 3, and the relative percentage errors (PE) and average relative percentage errors (APE) can be calculated, as shown in Table 3. Finally, we employ ARIMA (autoregressive composite moving average) model to test our results from the non-grey perspective. The estimation results of ARIMA model are shown in Table A1-4 in the Appendix. All the coefficients are statistically significant and the model is well fitted (R-squared is 0.850452). The ARIMA (2,1,1) model of time series is determined as follows: ΔEI ¼ 0:0311 þ 1:52ΔEI  0:94ΔEI þ ε ε t1 t2 t0:99 t1 The performance of the ARIMA model and the calculated relative errors and average relative errors are presented in Table 3. As illustrated in Table 3, there is great consistency between simulated values and real values for the improved DVCGM (1, N) model. For APE, which is the performance prediction index, its values of the improved DVCGM (1, N) is the smallest (1.76% in the in-sample periods and 1.92% in the out-sample periods) among all four models. 3.3 Comparison and discussion of the results Hysteresis As shown in Figure 4a, the overall trends of the GM (1, N), DVCGM (1, N) and improved effect of DVCGM (1, N) model and ARIMA model comply with the true curve to some extent. However, China’sEI the performance evaluation is distinct (Figure 4b). The APE value of the improved DVCGM forecasting (1, N) model is the smallest with the minimum fluctuation, demonstrating the efficacy and reliability of the model. It also suggests that the improved DVCGM (1, N) model is better than the traditional grey multivariable models. The ARIMA model has the largest values of APE among these four competing models. GM (1, N) model obtains the second largest values of APE, suggesting DVCGM (1, N) is the suboptimal choice. In view of the hysteresis effect of energy-saving policies, the improved DVCGM (1, N) model has a much lower error than the other three models. Therefore, we choose the improved DVCGM (1, N) model as the best model for predicting EI. 3.4 Forecasting the future energy intensity from 2018 to 2022 Based on the improved DVCGM (1, N) model, we can forecast the output value of China’sEI from 2018 to 2022 using data of technological progress. However, the data of technological progress from 2018 to 2022 is unknown and needs to be predicted in advance. The one-step rolling GM (1, 1) model is considered for the data prediction. We used the GM (1, 1) model for in-sample simulation and out-of-sample prediction of technological progress. The results (Table 4) show that the in-sample simulation error is 2.97%, and the out-of-sample hind cast error is 2.11%, which indicates that the GM (1, 1) model gives satisfying results and can be used to predict technological progress from 2018 to 2022. It is also noted that energy conservation policy (P) is introduced as a dummy variable. Before 2020, P value is 0; after 2020, P value is 1. It is because that the 13th Five-Year Plan is from 2016 to 2020, new performance target and new policies for the next stage are Original data ARIMA model GM(1,N)model DVCGM(1,N)model Improved DVCGM(1,N)model (APE) (Energy intensity values) 1.5 1.3 1.1 0.9 Figure 4. Performance 0.7 evaluation of the four models. (a) Estimated values of energy 0.5 intensity (Unit: tons of 2001 2004 2007 2010 2013 2016 2001 2004 2007 2010 2013 2016 standard coal per Years Years 10,000 Yuan). (b) APEs of the four models (a) (b) GS Technical progress 11,3 Years Original value Simulated value Error 2001 0.3464 0.3464 0 2002 0.3465 0.385443 0.112389 2003 0.3616 0.376489 0.041174 2004 0.3821 0.367743 0.037575 2005 0.384 0.3592 0.064584 2006 0.3724 0.350855 0.057854 2007 0.3563 0.342705 0.038157 2008 0.3367 0.334743 0.005812 2009 0.3297 0.326967 0.00829 2010 0.3191 0.319371 0.00085 2011 0.3086 0.311952 0.010862 2012 0.3075 0.304705 0.009089 2013 0.2958 0.297627 0.006175 2014 0.2843 0.290713 0.022555 APE 0.029669 2015 0.2750 0.283959 0.032578 2016 0.2769 0.277362 0.00167 Table 4. 2017 0.2633 0.270919 0.028937 Simulation results of technical progress APE 0.021062 Years Technical progress Energy conservation policy 2018 0.26462549 0 2019 0.25847804 0 2020 0.25247341 1 Table 5. 2021 0.24660827 1 Prediction results of technological progress 2022 0.24087938 1 Table 6. Predicted value of China’s energy intensity by using the improved DVCGM Years 2018 2019 2020 2021 2022 (1, N) with one-step Energy intensity 0.5245 0.4893 0.4674 0.4358 0.4258 rolling mechanism formulating and will be effective after 2020. Therefore, the data we need to substitute into the improved DVCGM (1, N) model is shown in Table 5. The predicted results of China’s EI are shown in Table 6. In addition, we also draw a line graph to make the results more iconic in Figure 5. As illustrated in Figure 5, there is a downward trend of the EI in the next five years. By 2020, the EI is expected to decrease by 20% or more than it was in 2016. That is to say, during the 13th Five-Year Plan period (2016–2020), EI will drop by more than 15%, meeting the country’s energy performance target. Therefore, government policies have a profound influence on EI. When formulating energy conservation and emission reduction policies, we should consider the hysteresis effect of the policies and make adjustments accordingly to achieve the goal. (tons standard coal/10,000 yuan) Hysteresis 1.6 effect of 1.4 China’sEI Forecasted data 1.2 forecasting Original data 0.8 0.6 0.4 Figure 5. 0.2 The growth trend of energy intensity 2001 2003 2005 2007 2009 2011 2013 2015 2017 2019 2021 of China Year 4. Conclusions Over the past 20 years, China is gradually shifting from a resource-intensive and energy- driven economy to a more sustained economy. EI in China has fallen almost continuously while China focuses on the industrial upgrading and promoting transformation of the economic structure. Energy-saving policies and regulations are introduced to help China reach its energy performance target. However, few studies have been carried out to consider the hysteresis effects of policies on the estimation of EI. Therefore, to address such a challenge problem, an improved grey multivariable model is designed to forecast China’sEI considering the hysteresis effect of government policies. To further improve its forecasting capability, a nonlinear optimization method based on PSO algorithm is constructed to calculate the hysteresis parameter. In addition, three conventional models, namely GM (1, N), DVCGM (1, N) and ARIMA models, are applied to test the accuracy of this improved DVCGM (1, N) model. The empirical results demonstrate that the proposed model considering the hysteresis effects of energy conservation policies performs best and matches well with the actual observations. Accordingly, this proposed model is used to forecast EI value from 2018 to 2022. The main conclusions are as follows: (1) The improved DVCGM (1, N) model can solve the modeling problem of small sample systems with time-delay causality. A nonlinear optimization method based on PSO algorithm is constructed to calculate the hysteresis parameter. It overcomes the defects of traditional GMs and econometric models. (2) GM (1, N), DVCGM (1, N) and ARIMA model are taken as comparative models. The accuracy of improved DVCGM (1, N) model was tested by the average relative percentage errors. The results show that the Improved DVCGM (1, N) model notes the hysteresis effect of government policies and significantly improves the prediction accuracy of China’s EI than the other three models. As suggested by APEs, the overall fitting in descending order is improved DVCGM (1, N) model, DVCGM (1, N), GM (1, N) and ARIMA model. (3) China’s EI is greatly influenced by technological progress and is much of policy- driven. When formulating energy conservation and emission reduction policies, we should fully consider the hysteresis effect of the policies, so as to make adjustment of the relative policies and better achieve the national energy performance target. A few caveats are appropriate. It is an interesting further path to work out the hysteresis parameter directly from the nonlinear programming model. Besides, the sustainability of the Energy intensity hysteresis effect of policy is worth considering. Furthermore, population factors and GS industrial structure also have good correlation with EI. These will be investigated in our 11,3 further studies. References Chen, K., Zhou, F. and Liu, A. (2018), “Chaotic dynamic weight particle swarm optimization for numerical function optimization”, Knowledge-Based Systems, Vol. 8 No. 10, pp. 23-40. Dang, Y.G., Wei, L. and Ding, S. (2017), “Grey multivariable discrete delay model based on driver information control and its application”, (in Chinese with English abstract), Control and Decision making, Vol. 32 No. 9, pp. 1672-1680. Deng, J.L. (1982), “Control problems of grey system”, Systems and Control Letters, Vol. 1 No. 5, pp. 288-294. Ding, S., Dang, Y.G., Li, X.M., Wang, J.J. and Zhao, K. (2017), “Forecasting Chinese CO emissions from fuel combustion using a novel grey multivariable model”, Journal of Cleaner Production, Vol. 162 No. 3, pp. 1527-1538. Ding, S., Dang, Y.G. and Xu, N. (2018), “Construction and application of GM (1, N) model based on virtual variable control”, (in Chinese with English abstract), Control and Decision making, Vol. 33 No. 2, pp. 309-315. Dong, K., Sun, R., Hochman, G. and Li, H. (2018), “Energy intensity and energy conservation potential in China: a regional comparison perspective”, Energy, Vol. 155 No. 2, pp. 782-795. Guo, H., Xiao, X. and Forrest, J. (2013), “A research on a comprehensive adaptive grey prediction model CAGM (1, N)”, Applied Mathematics and Computation, Vol. 225 No. 3, pp. 216-227. Hao, Y., Wang, Y., Zhao, J. and Li, H. (2011), “Grey system model with time lag and application to simulation of karst spring discharge”, Grey Systems: Theory and Application, Vol. 1 No. 1, pp. 47-56. Karimu, A., Br€annlund, R. and Lundgren, T. (2017), “Energy intensity and convergence in Swedish industry: a combined econometric and decomposition analysis”, Energy Economics, Vol. 62 No. 5, pp. 347-356. Kiran and Mustafa, S.. (2017), “Particle swarm optimization with a new update mechanism”, Applied Soft Computing, Vol. 604, pp. 670-678. Kose, E. and Tasci, L. (2015), “Prediction of the vertical displacement on the crest of Keban Dam”, Journal of Grey System, Vol. 27 No. 1, pp. 12-20. Kose, E. and Tasci, L. (2019), “Geodetic deformation forecasting based on multi-variable grey prediction model and regression model”, Grey Systems: Theory and Application, Vol. 9 No. 4, pp. 464-471. Ma, B. and Yu, Y. (2017), “Industrial structure, energy-saving regulations and energy intensity: evidence from Chinese cities”, Journal of Cleaner Production, Vol. 141 No. 3, pp. 1539-1547. Mason, K., Duggan, J. and Howley, E. (2018), “A meta optimisation analysis of particle swarm optimisation velocity update equations for watershed management learning”, Applied Soft Computing, Vol. 62 No. 10, pp. 148-161. Pao, H.T., Fu, H.C. and Tseng, C.L. (2012), “Forecasting of CO emissions, energy consumption and economic growth in China using an improved grey model”, Energy, Vol. 40 No. 1, pp. 400-409. Tan, R. and Lin, B. (2018), “What factors lead to the decline of energy intensity in China’s energy intensive industries?”, Energy Economics, Vol. 71 No. 10, pp. 213-221. Tien, T.L. (2005), “The indirect measurement of tensile strength of material by the grey prediction model GMC(1,N)”, Measurement Science and Technology, Vol. 16 No. 6, pp. 1322-1328. Tien, T.L. (2010), “Forecasting CO output from gas furnace by grey prediction model IGMC(1,N)”, Journal of the Chinese Society of Mechanical Engineers, Vol. 31 No. 1, pp. 55-65. Tien, T.L. (2011), “The indirect measurement of tensile strength by the new model FGMC(1,N)”, Hysteresis Measurement, Vol. 44 No. 10, pp. 1884-1897. effect of Wang, Z.X. and Hao, P. (2016), “An improved grey multivariable model for predicting industrial China’sEI energy consumption in China”, Applied Mathematical Modeling, Vol. 40 No. 11, pp. 5745-5758. forecasting Wu, L.F., Gao, X., Xiao, Y., Yang, Y. and Chen, X. (2018), “Using a novel multi-variable grey model to forecast the electricity consumption of Shandong Province in China”, Energy, Vol. 157 No. 10, pp. 327-335. Xiong, P.P., Zhang, Y., Xing, Z. and Chen, F. (2019), “Multivariable time-delay discrete MGM (1, m, tau) model and its application”, (in Chinese with English abstract), Statistics and Decision, Vol. 35 No. 8, pp. 18-22. Yan, Q.Y. and Liu, F. (2011), “Analysis and prediction of industrial energy consumption intensity in China from 2010 to 2030”, (in Chinese with English abstract), East China electric power, Vol. 39 No. 12, pp. 1858-1862. Zhai, X. (1996), “GM (1, 2) model with time delay and its application”, Systems Engineering, Vol. 14 No. 6, pp. 56-62. Zhang, K. (2015), “Time-delayed multivariable discrete grey model and its application”, Systems Engineering Theory and Practice, Vol. 35 No. 8, pp. 2092-2103. Zhang, K. (2016), “Multivariate discrete grey model base on dummy drivers”, Grey Systems: Theory and Application, Vol. 6 No. 2, pp. 246-258. Zhu, J.M. and Ruth, M. (2015), “Relocation or reallocation: impacts of differentiated energy saving regulation on manufacturing industries in China”, Ecological Economics, Vol. 110 No. 1, pp. 119-133. Appendix GS 11,3 Null hypothesis: EI has a unit root Exogenous: constant Lag length: 1 (automatic - based on SIC, maxlag 5 1) t-statistic Prob.* AuGMented Dickey–Fuller test statistic 2.587646 0.1293 Test critical values 1% level 4.420595 Table A1. 5% level 3.259808 Unit root test results of EI of ARIMA model 10% level 2.771129 Null hypothesis: DEI has a unit root Exogenous: constant, linear trend Lag length: 2 (automatic - based on SIC, maxlag 5 2) t-statistic Prob.* AuGMented Dickey–Fuller test statistic 5.277473 0.0102 Test critical values 1% level 5.295384 Table A2. 5% level 4.008157 Unit root test results of DEI of ARIMA model 10% level 3.460791 Autocorrelation Partial correlation AC PAC Q-stat Prob . j***** j . j***** j 1 0.714 0.714 8.2763 0.004 . j*. j *****j . j 2 0.195 0.641 8.9492 0.011 .**j . j . j . j 3 0.219 0.035 9.8859 0.020 .***j . j . j . j 4 0.366 0.010 12.796 0.012 .**j . j .*j . j 5 0.297 0.066 14.940 0.011 .*j . j . j . j 6 0.124 0.030 15.367 0.018 . j . j .**j . j 7 0.039 0.267 15.417 0.031 . j . j . j*. j 8 0.024 0.139 15.439 0.051 .*j . j .**j . j 9 0.087 0.284 15.805 0.071 Table A3. .*j . j . j*. j 10 0.136 0.079 17.004 0.074 The autocorrelation .*j . j . j . j 11 0.100 0.033 17.982 0.082 and partial correlation of ARIMA model . j . j .*j . j 12 0.017 0.188 18.037 0.115 Hysteresis Variable Coefficient Std. Error t-statistic Prob effect of C 0.031125 0.007696 4.044042 0.0037 China’sEI AR(1) 1.519400 0.137896 11.01843 0.0000 forecasting AR(2) 0.936168 0.059890 15.63139 0.0000 MA(1) 1.000000 39728.31 2.52E-05 1.0000 SIGMASQ 0.000442 0.403940 0.001093 0.9992 R-squared 0.850452 Mean dependent var 0.018404 Adjusted R-squared 0.775678 S.D. dependent var 0.056562 S.E. of regression 0.026789 Akaike info criterion 3.664240 Sum squared resid 0.005741 Schwarz criterion 3.446952 Log likelihood 28.81756 Hannan-Quinn criterion 3.708902 Table A4. F-statistic 11.37363 Durbin–Watson stat 1.064574 Estimation of Prob(F-statistic) 0.002202 ARIMA model Corresponding author Xinchang Guo can be contacted at: guo201064@163.com For instructions on how to order reprints of this article, please visit our website: www.emeraldgrouppublishing.com/licensing/reprints.htm Or contact us for further details: permissions@emeraldinsight.com http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Grey Systems: Theory and Application Emerald Publishing

Forecasting China's energy intensity by using an improved DVCGM (1, N) model considering the hysteresis effect

Loading next page...
 
/lp/emerald-publishing/forecasting-china-s-energy-intensity-by-using-an-improved-dvcgm-1-n-Z0OGwsNbic
Publisher
Emerald Publishing
Copyright
© Zhaosu Meng, Xiaotong Liu, Kedong Yin, Xuemei Li and Xinchang Guo
ISSN
2043-9377
DOI
10.1108/gs-02-2020-0022
Publisher site
See Article on Publisher Site

Abstract

Purpose – The purpose of this paper is to examine the effectiveness of an improved dummy variables control grey model (DVCGM) considering the hysteresis effect of government policies in China’s energy intensity (EI) forecasting. Design/methodology/approach – Energy consumption is considered as an important driver of economic development. China has introduced policies those aim at the optimization of energy structure and EI. In this study, EI is forecasted by an improved DVCGM, considering the hysteresis effect of energy-saving policies of the government. A nonlinear optimization method based on particle swarm optimization (PSO) algorithm is constructed to calculate the hysteresis parameter. A one-step rolling mechanism is applied to provide input data of the prediction model. Grey model (GM) (1, N), DVCGM (1, N) and ARIMA model are applied to test the accuracy of the improved DVCGM (1, N) model prediction. Findings – The results show that the improved DVCGM provides reliable results and works well in simulation and predictions using multivariable data in small sample size and time-lag virtual variable. Accordingly, the improved DVCGM notes the hysteresis effect of government policies and significantly improves the prediction accuracy of China’s EI than the other three models. Originality/value – This study estimates the EI considering the hysteresis effect of energy-saving policies in China by using an improved DVCGM. The main contribution of this paper is to propose a model to estimate EI, considering the hysteresis effect of energy-saving policies and improve forecasting accuracy. Keywords Grey system theory, Energy intensity forecasting, Grey prediction model, Hysteresis effect Paper type Research paper © Zhaosu Meng, Xiaotong Liu, Kedong Yin, Xuemei Li and Xinchang Guo. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) licence. Anyone may reproduce, distribute, translate and create derivative works of this article (for both Grey Systems: Theory and Application commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this licence may be seen at http://creativecommons.org/licences/by/4.0/ pp. 372-393 Emerald Publishing Limited legalcode. 2043-9377 DOI 10.1108/GS-02-2020-0022 This work was supported by the National Social Science Foundation of China (Grant No. 19BJY151). 1. Introduction Hysteresis Energy consumption is considered as an important driver of economic growth. Energy effect of intensity (EI), which is the ratio of total energy consumption to gross domestic product (GDP) China’sEI of a country or region in a certain period of time, measures the performance of energy forecasting utilization. According to the energy section of the National Bureau of Statistics (NBS), China’s largest energy consumption is 4,640 million tons of standard coal, and China has become the world largest energy consumer. Entering a new normal period, the government focuses on a crucial rebalancing and diversifying economy, with higher requirements for sustainable development. In the 13th Five-Year Plan for economic and social development, a total reduction of 15% has been set as the energy performance target. It is imperative to estimate and predict the EI to evaluate energy conservation and emission reduction. However, due to the complexity and dynamics of social economy, the implementation of energy conservation and emission reduction policies will not immediately reduce EI, and there is often a certain time lag. This inevitable time-lag should be taken into account to estimate and forecast the EI in China accurately. Numerous research studies have studied the influencing factors and driving forces of EI through econometrics methods such as cointegration analysis, metrology and decomposition analysis and scenario analysis. Zhu et al. (2015) applied cointegration analysis method on the large and state-owned enterprises and found that energy-saving regulations in China are one of the most important factors in reducing aggregate EI. Karimu et al. (2017) studied the EI and convergence of Swedish industry by combining metrology and decomposition analytical methods. Ma and Yu (2017) used panel data model to discuss the driving factors which lead to EI decline. It is revealed that industrial structure, energy conservation regulations and EI are closely related. Tan et al. (2018) used index decomposition analysis and production decomposition analysis methods to analyze the factors which are related to the decline in the EI and pointed out that technology improvement effect is the most significant factor. In EI forecasting, Pao et al. (2012) used improved grey models (GMs) to predict China’sCO emissions, energy consumption and economic growth. Dong et al. (2018) estimated the driving force of regional EI in China and forecasted the potential of regional energy conservation with scenario analysis. Wu et al. (2018) used a new multivariable GM to predict energy consumption in Shandong Province. The grey system theory is an interdisciplinary theory proposed by Deng (1982) and the grey prediction method applies well to small size data forecasting. As an important part of grey prediction theory, the GM (1, N) model is the basic model of multivariable grey system modeling approach. In recent years, numerous scholars have thoroughly discussed the model parameters optimization (Tien, 2005, 2010), the model accuracy improvement (Tien, 2011; Wang et al., 2016) and the expansion of the GM (Guo et al., 2013; Kose and Tasci, 2015, 2019; Ding et al., 2017). Based on grey multivariable model with time-lagged system, Zhai et al. (1996) introduced the lag term into the GM (1, 2) model and determined the delay parameters with the goal of minimizing the modeling error. Hao (2011) used grey correlation analysis to determine time-lag period between variables and then on this basis to establish forecast model of GM (1, N). Zhang et al. (2015) constructed a time-delay multivariable discrete GM,DDGM (1, N) model, by introducing a time-delay control factor and solved the time-delay parameters by using the grey dimensionally expanding identification method, which obtained a good application effect. Ma and Yu (2017) used a novel time-delay multivariable GM to predict the natural gas consumption in China. Dang et al. (2017) constructed the discrete delay grey multivariable DDGMD (1, N) model by introducing the driving information control adjustment coefficient T and the action coefficient λ and solved the coefficient respectively by grey dimension expansion method and particle swarm optimization (PSO) algorithm. Xiong (2019) built a multivariable time-delay discrete MGM (1, m, t) model and studied the mechanism of modeling and the process of modeling, and the GS calculation method of time delay is given. The hysteresis effect is discussed by an example 11,3 to verify the validity of the model. GM (1, N) model, in spite of successfully applying in various fields, sometimes ignores the influence of virtual variables such as policy on the main system in practical applications. Zhang (2016) considered the influence of dummy variables on system behavior variables and built a discrete multivariable prediction model based on dummy variables, which further expanded the application scope of the model. Ding et al. (2018) introduced the dummy variables into the GM (1, N) model, gave the concrete model construction method in mathematics and verified the effectiveness of the new model with cases. Generally speaking, we have abundant literature discussing the influencing factors of EI and its forecasting. Furthermore, there are some practices targeting the time-lag phenomenon using grey theory. The existing literature has explored and studied the GM from multiple angles, but there is still room for the GM to expand in the combination of dummy variables and time-delay systems. As the common GM (1, N) model does not take into account the time delay between variables, and there are dummy variables in the system which are difficult to be measured by quantity, our optimized DVCGM (1, N) model, which is short for dummy variables control grey model of N variables, incorporates the above two features and improves prediction accuracy to a satisfying level. It is practical significant because the time delay of dummy variables, i.e. government policy, is often seen in real world but is difficult to measure. The optimized DVCGM (1, N) model takes these variables into consideration, solves the practical problem and expands the grey system theory and the grey prediction method system, which improves the accuracy of grey prediction model. This paper’s main interest is to estimate and forecast EI by considering the influence of government policies and to test the accuracy of the improved GM through a comparison study. The main contribution of the study to the literature is to consider hysteresis effect and increase the forecast accuracy. It is imperative for optimizing energy structure, improving energy utilization efficiency and ensuring energy security. The results show that the improved GM produces better results than the other three conventional models. The rest of the paper is organized as follows. Section 2 briefly introduces grey theory and multivariable grey prediction models. A nonlinear optimization method based on PSO algorithm is constructed to calculate the hysteresis parameter in the improved model. In Section 3,EIis estimated by considering the hysteresis effect of energy-saving policies by an improved DVCGM (1, N) model. We compare the results with the other two GMs and one econometric model. The improved DVCGM (1, N) model has the best performance in the estimation comparison and is applied to forecast future EI. Section 4 is the conclusion of the study with the limitation and future path. 2. Methodology The three grey prediction models used in this paper are interrelated and in a progressive order. The GM (1, N) model is a traditional multivariable grey prediction model. The DVCGM (1, N) model introduces virtual variables, taking into account the influence of policy and other factors on the basis of the GM (1, N) model. Time-delay parameter is introduced in the improved DVCGM (1, N) model, considering the hysteresis effects of historical variables, which further enriches the existing grey prediction theory. 2.1 GM (1, N) model Grey prediction model can be regarded as two levels of work. At a lower level, the original sequence produces the sequence of generation by one accumulative generation (1-AGO), and then it forms the sequence of mean generation of consecutive neighbors; similarly, the Hysteresis sequence of influencing factors generates the sequence of generation by 1-AGO; Constructing effect of B and Y matrix and calculating system parameters through ordinary least squares (OLS) China’sEI regression. Once the system parameters a and b are determined, we can obtain the time forecasting response function (TRF) by solving the whitenization equation. At a higher level, the continuous differential equation with initial values is used as the reflection equation, and the discrete data are mapped to a manageable function, which is further restored to the TRF as the simulation basis. Typical procedures can be described briefly by the program in Figure 1. Definition 1. The original sequence is ð0Þ ð0Þ ð0Þ ð0Þ X ¼ x ð1Þ; x ð2Þ; :::; x ðnÞ 1 1 1 1 ð1Þ ð0Þ ð1Þ X is the 1 - AGO sequence from X , and X is the sequences of the relevant factors, where 1 1 i ð1Þ ð0Þ X ðkÞ¼ x ðmÞ i i m¼1 ð1Þ ð1Þ Z is the sequence of mean generation of consecutive neighbors from X , where 1 1 ð1Þ ð1Þ ð1Þ Z ðkÞ¼ 0:5X ðk  1Þþ 0:5X ðkÞ; ðk ¼ 2; 3; ... ; nÞ 1 1 1 Definition 2. Denote Eq (1) as the definition formula of GM (1, N) model: ð0Þ ð1Þ ð1Þ X ðkÞþ aZ ðkÞ¼ b X ðkÞ (1) 1 1 i i¼2 where k 5 2, 3, n, a is the development coefficient and b ; b ;  ; b are the grey input 1 2 n coefficients obtained by the least squares method. To determine these coefficients, the matrix B and Y are defined as follows: 2 3 2 3 ð0Þ ð1Þ ð1Þ ð1Þ z ð2Þ X ð2Þ X ð2Þ x ð2Þ 1 1 2 n 6 7 6 7 ð0Þ 6 ð1Þ ð1Þ 7 ð1Þ 6 7 x ð3Þ z ð3Þ X ð3Þ X ð3Þ 6 7 1 1 2 n 6 7 6 7 Y ¼ 6 7; B ¼ 6 7 . . . . 6 7 . . . . 6 7 4 . 5 . . . 4 5 ð0Þ ð1Þ ð1Þ ð1Þ x ðnÞ z ðnÞ X ðnÞ X ðmÞ 1 1 2 n The values of the coefficients a and b ; b ;  ; b can be determined by the following equation: 1 2 n 2 3 6 b 7 T −1 T b 6 7 P ¼ ¼ðB BÞ B Y 4 5 Definition 3. Let Eq (2) be defined as the differential equation (or) called grey reflection equation: ð0Þ n dx ð0Þ ð1Þ þ ax ðtÞ¼ b x ðtÞ (2) dt i¼1 GS 11,3 Figure 1. Principal diagram of GM (1, N) model Initialization Immediate original Generating mean sequence sequence1 sequence Grey equation Influencing Generating factor Time response Whitenization sequence2 function sequence Whitenization equation system BMatrix OLS parameters a and Y Matrix and b Prediction After determining the coefficients of a and b ; b ;  ; b , the differential equation of the GM 1 2 n Hysteresis can be determined by Eq (2). The solution of the above differential equation is as follows: effect of ( ) n n X X ð1Þ China’sEI ð0Þ ð1Þ b ð1Þ −ak b b b X ðk þ 1Þ¼ X ð1=aÞ b X ðkÞ 3e þð1=aÞ b X ðkÞ (3) i−1 i−1 1 1 i i forecasting i¼2 i¼2 ð1Þ where X ðkÞ is the prediction of the AGO of the original sequence. By considering that the estimation of the first element of the first AGO of a sequence is equal to the first element of the ð1Þ ð0Þ sequence, the following relation is determined: X ð1Þ¼ X ð1Þ 1 1 Finally, in order to predict the elements of the original sequence, the inverse accumulated generating operation should be performed. Therefore, the predicted values can be determined as follows: ð0Þ ð1Þ ð1Þ b x ðkÞ¼ b x ðkÞ b x ðk  1Þ; k≥ 2; 1 1 1 ð0Þ where x ðnÞ is an estimation of the original sequence, which is simulation values, ð0Þ ð0Þ b x ðn þ 1Þ;b x ðn þ 2Þ;  are predictive values. 1 1 2.2 DVCGM (1, N) model Traditional GM (1, N) model ignores the influence of virtual variables, which will inevitably lead to significant errors in practical applications. Therefore, it is necessary to construct a new multivariable predictive model with virtual variable control, based on the traditional GM (1, N) model, i.e. the DVCGM (1, N) model. The modeling steps for DVCGM (1, N) can be illustrated in Figure 2. Definition 4. The original sequence is ð0Þ ð0Þ ð0Þ ð0Þ X ¼ x ð1Þ; x ð2Þ;  ; x ðnÞ : i i i i Virtual variable sequence is ð0Þ ð0Þ ð0Þ ð0Þ ð0Þ D ¼ d ð1Þ; d ð2Þ; d ðnÞ ;  ; d ðnÞ¼ 0or 1: j j j j j ð1Þ ð1Þ ð0Þ X and D are the 1-AGO sequence, X is the behavior sequence of the system, i j 1 ð1Þ ð1Þ X ði ¼ 2;  ; MÞ and D ðj ¼ M þ 1;  ; NÞ is the driving factor sequence. Then GM i j (1, N) model with dummy variable can be expressed as: M N X X ð0Þ ð1Þ ð1Þ ð1Þ x ðkÞþ az ðkÞ¼ b x ðkÞþ b d ðkÞ (4) i j 1 1 i j i¼2 j¼Mþ1 named DVCGM (1, N) model (dummy variables control grey model of N variables). ð1Þ ð1Þ Z is the sequence of mean generation of consecutive neighbors from X where, 1 1 ð1Þ ð1Þ ð1Þ Z ðkÞ¼ 0:5X ðk− 1Þþ 0:5X ðkÞ; ðk ¼ 2; 3; ... ; nÞ 1 1 1 M N P P ð1Þ ð1Þ where b x ðkÞ is independent quantization variable driver, b d ðkÞ is the dummy i j i j i¼2 j¼Mþ1 variable driver. If virtual variables are not considered or b 50, the model can be transformed to a traditional GM (1, N) model. GS 11,3 Figure 2. Principal diagram of DVCGM (1, N) mode Initialization (0) (1) (1) 1 1 Z (0) (1) Grey equation X (0) (1) Time response D D Whitenization j j function TT 1 OLS B Y Matrix PB =()B BY Whitenization n equation Prediction Theorem 1. Assuming Hysteresis ð0Þ ð1Þ ð0Þ ð1Þ effect of X ; X ; D ðkÞ; D ðkÞ i i j j China’sEI as mentioned in Definition 4. The parameter column of the model is forecasting b ¼½a; b ;  b  : 2 N The matrix B and Y are defined as follows: 2 3 2 3 ð1Þ ð1Þ ð1Þ ð1Þ ð1Þ ð0Þ z ð2Þ X ð2Þ X ð2Þ d ð2Þ d ð2Þ x ð2Þ 1 2 1 M Mþ1 N 6 7 6 7 6 ð1Þ ð1Þ ð1Þ ð1Þ ð1Þ 7 ð0Þ 6 7 x ð3Þ z ð3Þ X ð3Þ X ð3Þ d ð3Þ d ð3Þ 6 7 1 1 2 M Mþ1 N 6 7 6 7 Y ¼ 6 7; B ¼ 6 7 . . . . . . 6 7 6 . . . . . 7 4 . 5 . . . . . 4 5 ð0Þ ð1Þ ð1Þ ð1Þ ð1Þ ð1Þ x ðnÞ z ðnÞ X ðnÞ X ðnÞ d ðnÞ d ðnÞ 1 1 2 M Mþ1 N Definition 5. Assuming ð0Þ ð1Þ ð0Þ ð1Þ X ; X ; D ðkÞ; D ðkÞ i i j j As mentioned in Definition 4, the parameter column of the model is b ¼½a; b ;  ; b  : 2 N Let Eq (5) be defined as the differential equation of DVCGM (1, N) model: ð1Þ M N X X dx ð1Þ ð1Þ ð1Þ þ ax ðtÞ¼ b x ðtÞþ b d ðtÞ (5) i j 1 i j dt i¼2 j¼Mþ1 After determining the coefficients of a and b ; b ;  ; b , the differential equation of the GM 1 2 n can be determined by Eq (5). The solution of the above differential equation is as follows: ( " # M N X X ð1Þ ð0Þ ð1Þ ð1Þ −at x ðtÞ¼ e X  t b x ð0Þþ b d ð0Þ i j 1 1 i j i¼2 j¼Mþ1 " # ) N M N X X X ð1Þ ð1Þ at þ b x ðtÞþ b d ðtÞ e dt (6) i j i j i¼2 i¼2 j¼Mþ1 Finally, in order to predict the elements of the original sequence, the inverse accumulated generating operation should be performed. Therefore, the predicted values can be determined as follows: ð0Þ ð1Þ ð1Þ b b b x ðkÞ¼ x ðkÞ x ðk  1Þ; k≥ 2; 1 1 1 ð0Þ where b x ðnÞ is an estimation of the original sequence, which is simulation values, ð0Þ ð0Þ b b x ðn þ 1Þ; x ðn þ 2Þ;  are predictive values. 1 1 2.3 Improved DVCGM (1, N) model The classic multivariable grey prediction models, such as traditional GM (1, N) and DVCGM (1, N) models, can reflect the influences of current driving-variables on the present system behavior and innately ignore the hysteresis effect of historical variables. Therefore, an GS improved DVCGM (1, N) model is proposed, integrating these above prediction model. 11,3 2.3.1 Construction of the improved DVCGM (1, N) model. In this section, the hysteresis parameter λ is innovatively introduced into the DVCGM (1, N) model to improve the prediction accuracy. Supported by PSO algorithm, detailed process and algorithm can be described as following: Theorem 2. Assuming ð0Þ ð1Þ ð0Þ ð1Þ X ; X ; D ðkÞ; D ðkÞ i i j j As mentioned in Definition 4. The parameter column of the model is b ¼½a; b ;  ; b  : 2 N The matrix B and Y are defined as follows: 2 3 ð0Þ x ð2Þ 6 7 ð0Þ 6 7 x ð3Þ 6 7 Y ¼ 6 7; 6 7 4 5 ð0Þ x ðnÞ 2 3 2 2 2 2 X X X X ð1Þ 2−j ð1Þ 2−j ð1Þ 2−j ð1Þ 2−j ð1Þ Z ð2Þ; λ x ðjÞ;  λ x ðjÞ; λ d ðjÞ λ d ðjÞ 1 2 2 M M Mþ1 Mþ1 N N 6 7 j¼1 j¼1 j¼1 j¼1 6 7 6 7 3 3 3 3 6 7 X X X X 6 ð1Þ ð1Þ ð1Þ ð1Þ ð1Þ 7 3−j 3−j 3−j 3−j Z ð3Þ; λ x ðjÞ;  λ x ðjÞ; λ d ðjÞ λ d ðjÞ 6 7 1 2 2 M M Mþ1 Mþ1 N N 6 7 B ¼ j¼1 j¼1 j¼1 j¼1 6 7 6 7 6 . . . . . 7 . . . . . 6 7 . . . . . 6 7 n n n n 4 X X X X 5 ð1Þ ð1Þ ð1Þ ð1Þ ð1Þ n−j n−j n−j n−j Z ðnÞ; λ x ðjÞ; ... λ x ðjÞ; λ d ðjÞ ... λ d ðjÞ 1 2 2 M M Mþ1 Mþ1 N N j¼1 j¼1 j¼1 j¼1 The least square estimation of the parameter column satisfies the following requirements: −1 (1) When n 5 Nþ1, b ¼ B Y ; jBj≠ 0; −1 T T T (2) When n > Nþ1, b ¼ðB BÞ B Y ;B B≠ 0; −1 T T T (3) When n < Nþ1, b ¼ B ðB BÞ Y ;B B≠ 0; Proof: Substitute k 5 2, 3,..., n into the model, you can get M k N k X X X X ð1Þ ð1Þ ð1Þ k−j k−j ð1Þ x ð2Þ¼ −ax ð2Þþ b λ x ð2Þþ b λ d ð2Þ i q 1 1 i i q q i¼2 j¼1 q¼Mþ1 j¼1 M k N k X X X X ð1Þ ð1Þ k−j ð1Þ k−j ð1Þ x ð3Þ¼ −ax ð3Þþ b λ x ð3Þþ b λ d ð3Þ i q 1 1 i i q q i¼2 j¼1 q¼Mþ1 j¼1 M k N k X X X X ð1Þ ð1Þ ð1Þ k−j k−j ð1Þ x ðnÞ¼ −ax ðnÞþ b λ x ðnÞþ b λ d ðnÞ i q 1 1 i i q q i¼2 j¼1 q¼Mþ1 j¼1 That is, by the least square method, Y ¼ Bb (1) When n5 Nþ1, B has an inverse matrix, the equations have a unique solution, we can Hysteresis −1 get b ¼ B Y ; jBj≠ 0; effect of China’sEI (2) When n > Nþ1, B is column full rank, the full rank decomposition of B is B 5 DC. forecasting Then the generalized inverse matrix of B can be obtained: −1  −1 −1 −1 T T T T T T T T B ¼ C CC D D D ; β ¼ C CC D D D Y ; Because B is a full rank matrix, C can be taken as a unit matrix, B 5 D,so −1 −1 −1 T T T T T T b ¼ C CC D D D Y ¼ B B B Y ; (3) When n < Nþ1, B is a row full rank matrix, D can be taken as a unit matrix, B5 C,so −1 −1 −1 T T T T T T b ¼ C CC D D D Y ¼ B B B Y ; Definition 6. Let Eq (7) be defined as the differential equation of improved DVCGM (1, N) model: Z Z ð1Þ M T N Q X X dx ðtÞ ð1Þ ð1Þ 1 t−s t−s ð1Þ þ ax ðtÞ¼ b λ x ðsÞds þ b λ d ðsÞds (7) i Q 1 i q q dt 0 0 i¼1 q¼Mþ1 Where k 5 2, 3, n, a is the development coefficient and b ; b ;  ; b are the grey input 1 2 n coefficients, λ is the hysteresis parameter. PSO algorithm is used to determine the coefficients of a, b ; b ;  ; b , and λ . Then the differential equation of the GM can be determined by 1 2 n i Eq (7). The solution of Eq (7) is as follows: ( " # M N X X ð1Þ ð0Þ ð1Þ −at ð1Þ x ðtÞ¼ e X  t b x ð0Þþ b d ð0Þ i q 1 1 i q i¼2 q¼Mþ1 " # ) N M k N k X X X X X k−j ð1Þ k−j ð1Þ at þ b λ x ðtÞþ b λ d ðtÞ e dt (8) i q i i q q i¼2 i¼2 j¼1 q¼Mþ1 j¼1 When the range of the driving factor sequence is small, the driver term can be viewed as a grey constant, and then the approximate TRF sequence of the grey differential equation of the model is " # M k N k X X X X ð1Þ ð1Þ k−j k−j ð1Þ x ðkÞ¼ b λ x ðtÞþ b λ d ðtÞ i q 1 i i q q i¼2 j¼1 q¼Mþ1 j¼1 " # M k N k −aðk−1Þ X X X X k−j ð1Þ ð1Þ k−j ð1Þ −aðk−1Þ b λ x ðtÞþ b λ d ðtÞ þ X e (9) i q i i q q 1 i¼2 j¼1 q¼Mþ1 j¼1 Finally, in order to predict the elements of the original sequence, the inverse accumulated GS generating operation should be performed. The predicted values can be determined as 11,3 follows: ð0Þ ð1Þ ð1Þ b x ðkÞ¼ b x ðkÞ b x ðk  1Þ; k≥ 2; 1 1 1 ð0Þ where b x ðnÞ is an estimation of the original sequence, which is simulation values, ð0Þ ð0Þ b x ðn þ 1Þ;b x ðn þ 2Þ;  are predictive values. 1 1 2.3.2 Estimation and optimization of hysteresis parameter in the improved DVCGM (1, N) model. The most important part of the improved DVCGM (1, N) model is estimating the time-lag parameter, which directly affects the accuracy of the model. However, the time-lag parameters must be determined in advance, followed by B and Y matrix construction and system parameters calculation through OLS. Once the system parameters a and b are determined, we can obtain the TRF and the simulation andpredictionvalue of the model. In this paper, a nonlinear optimization model is established by using the Least One Multiplication. Then, the time-lag parameter is determined. When the range of the driving factor sequence is small, Eq (9) is used as the TRF, λ can be solved by the following nonlinear programming model: " # M k N k X X X X ð1Þ ð1Þ k−j k−j ð1Þ b x ðkÞ¼ b λ x ðtÞþ b λ d ðtÞ i q 1 i i q q i¼2 j¼1 q¼Mþ1 j¼1 min (10) " # M k N k −aðk−1Þ X X X X k¼2 k−j ð1Þ k−j ð1Þ ð1Þ −aðk−1Þ b λ x ðtÞþ b λ d ðtÞ þ X e i q i i q q 1 i¼2 j¼1 q¼Mþ1 j¼1 8 ð0Þ ð1Þ ð1Þ b x ðkÞ¼ b x ðkÞ b x ðk  1Þ; k ≥ 2 > 1 1 1 " # M k N k > X X X X > 1 ð1Þ ð1Þ > k−j k−j ð1Þ b x ðkÞ¼ b λ x ðtÞþ b λ d ðtÞ i q 1 i i q q > a > i¼2 j¼1 q¼Mþ1 j¼1 " # M k N k −aðk−1Þ s:t: X X X X (11) k−j ð1Þ k−j ð1Þ ð1Þ −aðk−1Þ >  b λ x ðtÞþ b λ d ðtÞ þ x ð1Þe i q > i i q q 1 > i¼2 j¼1 q¼Mþ1 j¼1 b ¼½a; b ;  ; b > 2 N 0 < λ < 1; i ¼ 2; 3; ... ; M; M þ 1; ... ; N The model takes the relationship between structural parameters as the constraint condition and minimizes the average simulation relative error of the system characteristic variables, which can improve accuracy to the greatest extent. The above optimization problem can be solved by PSO (Kiran and Mustafa, 2017; Mason et al., 2018; Chen et al., 2018). According to Eq (10), a nonlinear optimization method based on PSO algorithm is constructed to obtain the hysteresis parameter. PSO sets a certain number of particles in feasible region to find the best location and can be used to seek optimal values of λ .Denote λ in Eq (10) and construct the fitness function of each particle, i i according to Eq (12). " # M k N k X X X X  Hysteresis ð1Þ ð1Þ k−j k−j ð1Þ b x ðkÞ¼ b λ x ðtÞþ b λ d ðtÞ i q 1 i i q q  effect of i¼2 j¼1 q¼Mþ1 j¼1 China’sEI " # forecasting M k N k −aðk−1Þ X X X X Fitnes λ ¼ (12) i e k−j ð1Þ k−j ð1Þ b λ x ðtÞþ b λ d ðtÞ i q k¼2 i i q q i¼2 j¼1 q¼Mþ1 j¼1 ð1Þ −aðk−1Þ þX e Obviously, the average simulation relative error of the system characteristic variable sequence varies depending on the number of lag periods. The value of lag period should be selected to make the average simulation relative error as small as possible. Therefore, the improved DVCGM (1, N) model can well describe the hysteresis effect between the system characteristic variables. Once the hysteresis parameter is determined, the structural parameters of the model are set accordingly, and the simulation and prediction results can be obtained according to Eq (9). 2.3.3 Modeling procedure. Detailed procedure of the improved DVCGM (1, N) model is illustrated as follows. ð0Þ ð1Þ Step 1. Collect raw data and establish original sequence X , X is the 1-AGO sequence 1 1 ð0Þ ð1Þ ð1Þ from X , X is the sequence of the relevant factors. Z is the sequence of mean 1 i 1 ð1Þ generation of consecutive neighbors from X . Then, determining virtual variable sequence ð0Þ ð0Þ ð0Þ ð0Þ ð0Þ D ¼ d ð1Þ; d ð2Þ;  d ðnÞ ; d ðnÞ¼ 0or 1: j j j j j ð1Þ D is also the sequence of the relevant factors. Step 2. Solving delay parameters λ by PSO according (10), then constructing vector Y and matrix B. Using the least square method, the values of the coefficients a and b ; b ;  ; b 1 2 n can be determined. Step 3. After considering the lag effect, the differential equation of the grey model can be determined by Eq (7). Step 4. The time response function of the new model is established to generate prediction data according to Eq (9). 3. Application Forecasting EI can be considered as a grey problem, because EI is greatly affected by technological progress, population factor, industrial structure and so on. These factors influence EI through a dynamic and complicated mechanism. The uncertain impacts and limited number of data provides a good basis for grey theory application. There are four procedures in this part, including data collection, parameter estimation, result comparisons and future forecasts. Three competing models, namely GM (1, N), DVCGM (1, N) and ARIMA model, are employed to test the accuracy of the improved DVCGM (1, N) on EI forecasting. The three GMs used in this paper are interrelated and are from basic to the advanced. The GM (1, N) model is a traditional grey multivariate prediction model. On its basis, DVCGM (1, N) model introduces dummy variables, taking into account the influence of policies and other factors. The improved DVCGM (1, N) model introduces time-lag parameters, which further enriches the existing grey prediction theory by considering the hysteresis effect of GS policies. The above GMs are used to estimate China’s EI, and the optimal model is selected by 11,3 comparing their performance, to predict China’s EI in the next five years. In addition, we use ARIMA model as a comparison of grey methods to show that the improved DVCGM (1, N) model is not only better than the traditional GM, but also better than the non-grey econometric model. 3.1 Variables selection and data collection The indicators selected in this paper are as follows. EI is the ratio of total energy consumption to GDP. Population factor is the employed population. The ratio of the added value of industrial production to the added value of energy consumed by the industry is used as the substitution variable of technological progress (Yan, 2011). Industrial structure is measured by the ratio of output value of each industry to GDP. For the consistency of the statistical scope, we choose a time scale of 2001–2017, and all data are collected from China Statistical Yearbook. With small sample size (17 periods’ real measurement values) and insufficient information, this case fits well with the grey system. First, we calculate the grey correlation between EI and population, technological progress and industrial structure. As shown in Table 1, technological progress has the highest correlation with EI. Therefore, we select technological progress as the driving variable, and EI is the system behavior variable. Then, we can build GM (1, N) model for EI forecast. As shown in Figure 3, China’s EI increased rapidly from 2001, rising sharply to a peak of 1.4542 tons standard coal per 10,000 Yuan in 2005 and declined steadily afterward. In 2017, the EI has decreased by 62.67%, leaving only 0.5428 tons of standard coal per 10,000 Yuan. The 11th Five-Year Plan in 2006 is a kind of watershed for Chinese EI. The government introduced a series of strict policies of energy conservation and emission reduction, and Relative correlation Energy intensity Population factors 0.6126 Table 1. Technological progress 0.9550 Relative correlation Industrial structure 0.7652 analysis tons standard coal per 10,000 Yuan 1.6 0.39 0.37 1.4 0.35 0.33 1.2 0.31 0.29 0.27 0.25 0.8 Figure 3. 0.23 Trends of energy intensity (dash line) 0.21 0.6 and technological 0.19 progress (solid line) 0.4 0.17 of China 2001 2003 2005 2007 2009 2011 2013 2015 2017 Energy intensity Technological progress China’s EI declined obviously. As a result, we need to put policy and the hysteresis effect of Hysteresis policy into consideration when estimate EI of China. effect of China’sEI 3.2 Simulation of energy intensity in China forecasting Four models are applied in the simulation of EI of China. First, we use the GM (1, N) model mentioned in Section 2.1 to predict EI. We select technological progress as the driving variable, EI as the system behavior variable and then establish the GM (1, N) model. According to 2.1, the TRF is obtained as follows: ð1Þ ð1Þ −1:8419ðk−1Þ −1:8419ðk−1Þ X ðkÞ¼ 4:2429X ðkÞ3 1  e þ e 1 2 The simulation results of GM (1, N) model are illustrated in Table 3, and the relative errors and average relative errors can be calculated, as shown in Table 3. Secondly, we build DVCGM (1, N) model for EI estimate. We select technological progress as the driving variable, EI as the system behavior variable. Energy conservation policy (P)is introduced as a dummy variable. Before 2006, as the strict energy-saving policies had not been implemented, P value is 0; after 2006, P value is 1. According to Section 2.2, the TRF is obtained as follows: ð1Þ ð1Þ ð1Þ −1:7372ðk−1Þ −1:7372ðk−1Þ X ðkÞ¼ 4:2830X ðkÞ 0:0232d ðkÞ 3 1  e d þ e 1 2 3 The simulation results of DVCGM (1, N) model are illustrated in Table 3, and the relative errors and average relative errors can be calculated, as shown in Table 3. Thirdly, we use the improved DVCGM (1, N) model to estimate EI. As discussed in Section 2.3, DVCGM (1, N) model takes consideration of the hysteresis effect of dummy variables. We select technological progress as the driving variable, EI as the system behavior variable and energy-saving policy as virtual variable. Different from DVCGM (1, N) model, we need to determine the delay parameters, along with the structure parameters and build the TRF to get the simulation results. (1) Determination of time-delay parameter. As stated by Eq (10-12), a nonlinear optimization method based on PSO algorithm is constructed to obtain the hysteresis parameter. According to the optimization model shown in Eq (9), the average relative percentage errors (APE) of the model under different lag periods is calculated in Table 2. When the hysteresis parameter is 2, the average relative percentage errors (APE) of the model is the smallest (0.9754%). Therefore, the value of time lag is two years. In the view of the complexity of socio-economic ecology, adaptive adjustments are made according to the changes of policies. These energy-saving policies influence economic operation through a certain transmission mechanism and gradually affect the EI. (2) Determination of structural parameters. The hysteresis parameter is substituted into the B matrix and Y matrix. According to Theorem 2, the matrix operation of least square regression is used to obtain the structural parameter values b , b . B 5 (1.251669, 0.938994). 1 2 Table 2. Average percentage Lag period 0 1 2 3 4 5 errors (APEs) of the model under different APE(%) 1.1107 1.0023 0.9754 1.0667 1.0989 0.9836 lag periods GS Improved 11,3 DVCGM (1, N) DVCGM (1, N) Years GM (1, N) model model model ARIMA model In- Original Simulation PE Simulation PE Simulation PE Simulation PE sample data values % values % values % values % 2001 1.3731 1.3731 0 1.3731 0 1.3731 0 1.3731 0 2002 1.3932 1.4108 1.26 1.40783 1.05 1.4195 1.89 1.4381 3.22 2003 1.4341 1.5118 5.42 1.50494 4.94 1.4866 3.66 1.5134 5.53 2004 1.4229 1.4697 3.29 1.45264 2.09 1.4502 1.92 1.4536 2.16 2005 1.4542 1.4955 2.84 1.48736 2.28 1.4847 2.1 1.5111 3.91 2006 1.3054 1.3319 2.03 1.34535 3.06 1.3271 1.66 1.361 4.26 2007 1.1525 1.1965 3.82 1.19042 3.29 1.1681 1.35 1.235 7.16 2008 1.0034 1.0611 5.75 1.05487 5.13 1.044 4.05 1.0605 5.69 2009 0.9629 1.0025 4.11 0.99564 3.4 0.9729 1.04 1.0488 8.92 2010 0.8732 0.9179 5.12 0.91127 4.36 0.8875 1.64 0.9616 10.1 2011 0.791 0.82 3.66 0.81457 2.98 0.804 1.64 0.8661 9.49 2012 0.7442 0.8106 8.92 0.77739 4.46 0.7518 1.02 0.7913 6.33 2013 0.7004 0.7648 9.19 0.76 8.51 0.7112 1.54 0.7575 8.15 2014 0.6612 0.678 2.54 0.6852 3.63 0.6683 1.08 0.6966 5.36 APE% 4.14 3.51 1.76 5.74 Table 3. Simulation of energy Out- Original Prediction PE Prediction PE Prediction PE Prediction PE intensity in China by sample data value % value % value % value % GM (1, N), DVCGM (1, N), the Improved 2015 0.6239 0.6577 5.41 0.64118 2.77 0.6334 1.52 0.6718 7.68 DVCGM (1, N) and 2016 0.5861 0.6076 3.67 0.61136 4.31 0.5943 1.4 0.6157 5.05 ARIMA model (Unit: 2017 0.5428 0.5676 4.57 0.56033 3.23 0.5582 2.83 0.5677 4.59 tons of standard coal per 10,000 Yuan) APE% 4.55 3.44 1.92 5.78 (3) Calculate the simulated and predicted values. Putting the values of the estimated structural parameter and hysteresis parameter in Eq (7-9), we can obtain the optimal TRF as follows: ð1Þ ð1Þ ð1Þ −1:7372ðk−1Þ −1:7372ðk−1Þ X ðkÞ¼ 4:2830X ðkÞ 0:0232d ðkÞ 3 1  e d þ e 1 2 3 Then, the simulation results of improved DVCGM (1, N) model are illustrated in Table 3, and the relative percentage errors (PE) and average relative percentage errors (APE) can be calculated, as shown in Table 3. Finally, we employ ARIMA (autoregressive composite moving average) model to test our results from the non-grey perspective. The estimation results of ARIMA model are shown in Table A1-4 in the Appendix. All the coefficients are statistically significant and the model is well fitted (R-squared is 0.850452). The ARIMA (2,1,1) model of time series is determined as follows: ΔEI ¼ 0:0311 þ 1:52ΔEI  0:94ΔEI þ ε ε t1 t2 t0:99 t1 The performance of the ARIMA model and the calculated relative errors and average relative errors are presented in Table 3. As illustrated in Table 3, there is great consistency between simulated values and real values for the improved DVCGM (1, N) model. For APE, which is the performance prediction index, its values of the improved DVCGM (1, N) is the smallest (1.76% in the in-sample periods and 1.92% in the out-sample periods) among all four models. 3.3 Comparison and discussion of the results Hysteresis As shown in Figure 4a, the overall trends of the GM (1, N), DVCGM (1, N) and improved effect of DVCGM (1, N) model and ARIMA model comply with the true curve to some extent. However, China’sEI the performance evaluation is distinct (Figure 4b). The APE value of the improved DVCGM forecasting (1, N) model is the smallest with the minimum fluctuation, demonstrating the efficacy and reliability of the model. It also suggests that the improved DVCGM (1, N) model is better than the traditional grey multivariable models. The ARIMA model has the largest values of APE among these four competing models. GM (1, N) model obtains the second largest values of APE, suggesting DVCGM (1, N) is the suboptimal choice. In view of the hysteresis effect of energy-saving policies, the improved DVCGM (1, N) model has a much lower error than the other three models. Therefore, we choose the improved DVCGM (1, N) model as the best model for predicting EI. 3.4 Forecasting the future energy intensity from 2018 to 2022 Based on the improved DVCGM (1, N) model, we can forecast the output value of China’sEI from 2018 to 2022 using data of technological progress. However, the data of technological progress from 2018 to 2022 is unknown and needs to be predicted in advance. The one-step rolling GM (1, 1) model is considered for the data prediction. We used the GM (1, 1) model for in-sample simulation and out-of-sample prediction of technological progress. The results (Table 4) show that the in-sample simulation error is 2.97%, and the out-of-sample hind cast error is 2.11%, which indicates that the GM (1, 1) model gives satisfying results and can be used to predict technological progress from 2018 to 2022. It is also noted that energy conservation policy (P) is introduced as a dummy variable. Before 2020, P value is 0; after 2020, P value is 1. It is because that the 13th Five-Year Plan is from 2016 to 2020, new performance target and new policies for the next stage are Original data ARIMA model GM(1,N)model DVCGM(1,N)model Improved DVCGM(1,N)model (APE) (Energy intensity values) 1.5 1.3 1.1 0.9 Figure 4. Performance 0.7 evaluation of the four models. (a) Estimated values of energy 0.5 intensity (Unit: tons of 2001 2004 2007 2010 2013 2016 2001 2004 2007 2010 2013 2016 standard coal per Years Years 10,000 Yuan). (b) APEs of the four models (a) (b) GS Technical progress 11,3 Years Original value Simulated value Error 2001 0.3464 0.3464 0 2002 0.3465 0.385443 0.112389 2003 0.3616 0.376489 0.041174 2004 0.3821 0.367743 0.037575 2005 0.384 0.3592 0.064584 2006 0.3724 0.350855 0.057854 2007 0.3563 0.342705 0.038157 2008 0.3367 0.334743 0.005812 2009 0.3297 0.326967 0.00829 2010 0.3191 0.319371 0.00085 2011 0.3086 0.311952 0.010862 2012 0.3075 0.304705 0.009089 2013 0.2958 0.297627 0.006175 2014 0.2843 0.290713 0.022555 APE 0.029669 2015 0.2750 0.283959 0.032578 2016 0.2769 0.277362 0.00167 Table 4. 2017 0.2633 0.270919 0.028937 Simulation results of technical progress APE 0.021062 Years Technical progress Energy conservation policy 2018 0.26462549 0 2019 0.25847804 0 2020 0.25247341 1 Table 5. 2021 0.24660827 1 Prediction results of technological progress 2022 0.24087938 1 Table 6. Predicted value of China’s energy intensity by using the improved DVCGM Years 2018 2019 2020 2021 2022 (1, N) with one-step Energy intensity 0.5245 0.4893 0.4674 0.4358 0.4258 rolling mechanism formulating and will be effective after 2020. Therefore, the data we need to substitute into the improved DVCGM (1, N) model is shown in Table 5. The predicted results of China’s EI are shown in Table 6. In addition, we also draw a line graph to make the results more iconic in Figure 5. As illustrated in Figure 5, there is a downward trend of the EI in the next five years. By 2020, the EI is expected to decrease by 20% or more than it was in 2016. That is to say, during the 13th Five-Year Plan period (2016–2020), EI will drop by more than 15%, meeting the country’s energy performance target. Therefore, government policies have a profound influence on EI. When formulating energy conservation and emission reduction policies, we should consider the hysteresis effect of the policies and make adjustments accordingly to achieve the goal. (tons standard coal/10,000 yuan) Hysteresis 1.6 effect of 1.4 China’sEI Forecasted data 1.2 forecasting Original data 0.8 0.6 0.4 Figure 5. 0.2 The growth trend of energy intensity 2001 2003 2005 2007 2009 2011 2013 2015 2017 2019 2021 of China Year 4. Conclusions Over the past 20 years, China is gradually shifting from a resource-intensive and energy- driven economy to a more sustained economy. EI in China has fallen almost continuously while China focuses on the industrial upgrading and promoting transformation of the economic structure. Energy-saving policies and regulations are introduced to help China reach its energy performance target. However, few studies have been carried out to consider the hysteresis effects of policies on the estimation of EI. Therefore, to address such a challenge problem, an improved grey multivariable model is designed to forecast China’sEI considering the hysteresis effect of government policies. To further improve its forecasting capability, a nonlinear optimization method based on PSO algorithm is constructed to calculate the hysteresis parameter. In addition, three conventional models, namely GM (1, N), DVCGM (1, N) and ARIMA models, are applied to test the accuracy of this improved DVCGM (1, N) model. The empirical results demonstrate that the proposed model considering the hysteresis effects of energy conservation policies performs best and matches well with the actual observations. Accordingly, this proposed model is used to forecast EI value from 2018 to 2022. The main conclusions are as follows: (1) The improved DVCGM (1, N) model can solve the modeling problem of small sample systems with time-delay causality. A nonlinear optimization method based on PSO algorithm is constructed to calculate the hysteresis parameter. It overcomes the defects of traditional GMs and econometric models. (2) GM (1, N), DVCGM (1, N) and ARIMA model are taken as comparative models. The accuracy of improved DVCGM (1, N) model was tested by the average relative percentage errors. The results show that the Improved DVCGM (1, N) model notes the hysteresis effect of government policies and significantly improves the prediction accuracy of China’s EI than the other three models. As suggested by APEs, the overall fitting in descending order is improved DVCGM (1, N) model, DVCGM (1, N), GM (1, N) and ARIMA model. (3) China’s EI is greatly influenced by technological progress and is much of policy- driven. When formulating energy conservation and emission reduction policies, we should fully consider the hysteresis effect of the policies, so as to make adjustment of the relative policies and better achieve the national energy performance target. A few caveats are appropriate. It is an interesting further path to work out the hysteresis parameter directly from the nonlinear programming model. Besides, the sustainability of the Energy intensity hysteresis effect of policy is worth considering. Furthermore, population factors and GS industrial structure also have good correlation with EI. These will be investigated in our 11,3 further studies. References Chen, K., Zhou, F. and Liu, A. (2018), “Chaotic dynamic weight particle swarm optimization for numerical function optimization”, Knowledge-Based Systems, Vol. 8 No. 10, pp. 23-40. Dang, Y.G., Wei, L. and Ding, S. (2017), “Grey multivariable discrete delay model based on driver information control and its application”, (in Chinese with English abstract), Control and Decision making, Vol. 32 No. 9, pp. 1672-1680. Deng, J.L. (1982), “Control problems of grey system”, Systems and Control Letters, Vol. 1 No. 5, pp. 288-294. Ding, S., Dang, Y.G., Li, X.M., Wang, J.J. and Zhao, K. (2017), “Forecasting Chinese CO emissions from fuel combustion using a novel grey multivariable model”, Journal of Cleaner Production, Vol. 162 No. 3, pp. 1527-1538. Ding, S., Dang, Y.G. and Xu, N. (2018), “Construction and application of GM (1, N) model based on virtual variable control”, (in Chinese with English abstract), Control and Decision making, Vol. 33 No. 2, pp. 309-315. Dong, K., Sun, R., Hochman, G. and Li, H. (2018), “Energy intensity and energy conservation potential in China: a regional comparison perspective”, Energy, Vol. 155 No. 2, pp. 782-795. Guo, H., Xiao, X. and Forrest, J. (2013), “A research on a comprehensive adaptive grey prediction model CAGM (1, N)”, Applied Mathematics and Computation, Vol. 225 No. 3, pp. 216-227. Hao, Y., Wang, Y., Zhao, J. and Li, H. (2011), “Grey system model with time lag and application to simulation of karst spring discharge”, Grey Systems: Theory and Application, Vol. 1 No. 1, pp. 47-56. Karimu, A., Br€annlund, R. and Lundgren, T. (2017), “Energy intensity and convergence in Swedish industry: a combined econometric and decomposition analysis”, Energy Economics, Vol. 62 No. 5, pp. 347-356. Kiran and Mustafa, S.. (2017), “Particle swarm optimization with a new update mechanism”, Applied Soft Computing, Vol. 604, pp. 670-678. Kose, E. and Tasci, L. (2015), “Prediction of the vertical displacement on the crest of Keban Dam”, Journal of Grey System, Vol. 27 No. 1, pp. 12-20. Kose, E. and Tasci, L. (2019), “Geodetic deformation forecasting based on multi-variable grey prediction model and regression model”, Grey Systems: Theory and Application, Vol. 9 No. 4, pp. 464-471. Ma, B. and Yu, Y. (2017), “Industrial structure, energy-saving regulations and energy intensity: evidence from Chinese cities”, Journal of Cleaner Production, Vol. 141 No. 3, pp. 1539-1547. Mason, K., Duggan, J. and Howley, E. (2018), “A meta optimisation analysis of particle swarm optimisation velocity update equations for watershed management learning”, Applied Soft Computing, Vol. 62 No. 10, pp. 148-161. Pao, H.T., Fu, H.C. and Tseng, C.L. (2012), “Forecasting of CO emissions, energy consumption and economic growth in China using an improved grey model”, Energy, Vol. 40 No. 1, pp. 400-409. Tan, R. and Lin, B. (2018), “What factors lead to the decline of energy intensity in China’s energy intensive industries?”, Energy Economics, Vol. 71 No. 10, pp. 213-221. Tien, T.L. (2005), “The indirect measurement of tensile strength of material by the grey prediction model GMC(1,N)”, Measurement Science and Technology, Vol. 16 No. 6, pp. 1322-1328. Tien, T.L. (2010), “Forecasting CO output from gas furnace by grey prediction model IGMC(1,N)”, Journal of the Chinese Society of Mechanical Engineers, Vol. 31 No. 1, pp. 55-65. Tien, T.L. (2011), “The indirect measurement of tensile strength by the new model FGMC(1,N)”, Hysteresis Measurement, Vol. 44 No. 10, pp. 1884-1897. effect of Wang, Z.X. and Hao, P. (2016), “An improved grey multivariable model for predicting industrial China’sEI energy consumption in China”, Applied Mathematical Modeling, Vol. 40 No. 11, pp. 5745-5758. forecasting Wu, L.F., Gao, X., Xiao, Y., Yang, Y. and Chen, X. (2018), “Using a novel multi-variable grey model to forecast the electricity consumption of Shandong Province in China”, Energy, Vol. 157 No. 10, pp. 327-335. Xiong, P.P., Zhang, Y., Xing, Z. and Chen, F. (2019), “Multivariable time-delay discrete MGM (1, m, tau) model and its application”, (in Chinese with English abstract), Statistics and Decision, Vol. 35 No. 8, pp. 18-22. Yan, Q.Y. and Liu, F. (2011), “Analysis and prediction of industrial energy consumption intensity in China from 2010 to 2030”, (in Chinese with English abstract), East China electric power, Vol. 39 No. 12, pp. 1858-1862. Zhai, X. (1996), “GM (1, 2) model with time delay and its application”, Systems Engineering, Vol. 14 No. 6, pp. 56-62. Zhang, K. (2015), “Time-delayed multivariable discrete grey model and its application”, Systems Engineering Theory and Practice, Vol. 35 No. 8, pp. 2092-2103. Zhang, K. (2016), “Multivariate discrete grey model base on dummy drivers”, Grey Systems: Theory and Application, Vol. 6 No. 2, pp. 246-258. Zhu, J.M. and Ruth, M. (2015), “Relocation or reallocation: impacts of differentiated energy saving regulation on manufacturing industries in China”, Ecological Economics, Vol. 110 No. 1, pp. 119-133. Appendix GS 11,3 Null hypothesis: EI has a unit root Exogenous: constant Lag length: 1 (automatic - based on SIC, maxlag 5 1) t-statistic Prob.* AuGMented Dickey–Fuller test statistic 2.587646 0.1293 Test critical values 1% level 4.420595 Table A1. 5% level 3.259808 Unit root test results of EI of ARIMA model 10% level 2.771129 Null hypothesis: DEI has a unit root Exogenous: constant, linear trend Lag length: 2 (automatic - based on SIC, maxlag 5 2) t-statistic Prob.* AuGMented Dickey–Fuller test statistic 5.277473 0.0102 Test critical values 1% level 5.295384 Table A2. 5% level 4.008157 Unit root test results of DEI of ARIMA model 10% level 3.460791 Autocorrelation Partial correlation AC PAC Q-stat Prob . j***** j . j***** j 1 0.714 0.714 8.2763 0.004 . j*. j *****j . j 2 0.195 0.641 8.9492 0.011 .**j . j . j . j 3 0.219 0.035 9.8859 0.020 .***j . j . j . j 4 0.366 0.010 12.796 0.012 .**j . j .*j . j 5 0.297 0.066 14.940 0.011 .*j . j . j . j 6 0.124 0.030 15.367 0.018 . j . j .**j . j 7 0.039 0.267 15.417 0.031 . j . j . j*. j 8 0.024 0.139 15.439 0.051 .*j . j .**j . j 9 0.087 0.284 15.805 0.071 Table A3. .*j . j . j*. j 10 0.136 0.079 17.004 0.074 The autocorrelation .*j . j . j . j 11 0.100 0.033 17.982 0.082 and partial correlation of ARIMA model . j . j .*j . j 12 0.017 0.188 18.037 0.115 Hysteresis Variable Coefficient Std. Error t-statistic Prob effect of C 0.031125 0.007696 4.044042 0.0037 China’sEI AR(1) 1.519400 0.137896 11.01843 0.0000 forecasting AR(2) 0.936168 0.059890 15.63139 0.0000 MA(1) 1.000000 39728.31 2.52E-05 1.0000 SIGMASQ 0.000442 0.403940 0.001093 0.9992 R-squared 0.850452 Mean dependent var 0.018404 Adjusted R-squared 0.775678 S.D. dependent var 0.056562 S.E. of regression 0.026789 Akaike info criterion 3.664240 Sum squared resid 0.005741 Schwarz criterion 3.446952 Log likelihood 28.81756 Hannan-Quinn criterion 3.708902 Table A4. F-statistic 11.37363 Durbin–Watson stat 1.064574 Estimation of Prob(F-statistic) 0.002202 ARIMA model Corresponding author Xinchang Guo can be contacted at: guo201064@163.com For instructions on how to order reprints of this article, please visit our website: www.emeraldgrouppublishing.com/licensing/reprints.htm Or contact us for further details: permissions@emeraldinsight.com

Journal

Grey Systems: Theory and ApplicationEmerald Publishing

Published: Jun 18, 2021

Keywords: Grey system theory; Energy intensity forecasting; Grey prediction model; Hysteresis effect

References