Research Article  Open Access
Forecasting Primary Energy Requirements of Territories by Autoregressive Integrated Moving Average and Backpropagation Neural Network Models
Abstract
Forecasting energy data, especially the primary energy requirement, is the key part of policymaking. For those territories of different developing types, seeking a knowledgebased and dependable forecasting model is an essential prerequisite for the prosperous development of policymaking. In this paper, both autoregressive integrated moving average and backpropagation neural network models which have been proved to be very efficient in forecasting are applied to the forecasts of the primary energy consumption of three different developing types of territories. It is shown that the average relative errors between the actual data and simulated value are from 4.5% to 5.9% by the autoregressive integrated moving average and from 0.04% to 0.47% by the backpropagation neural network. Specially, this research shows that the backpropagation neural network model presents a better prediction of primary energy requirement when considering gross domestic product, population, and the particular values as predictors. Furthermore, we indicate that the singleinput backpropagation neural network model can still work when the particular values have contributed most to the energy consumption.
1. Introduction
Since the emergence of the issue about applying mathematical models to predict energy consumption, relevant researches have been carried out [1] and several related experiments have been implemented [2]. The recent study [3] has combined the nonlinear metabolic grey model (NMGM) and autoregressive integrated moving average (ARIMA) model and has used the linear ARIMA to correct NMGM forecasting residuals, which has improved forecasting accuracy steadily and given useful policy recommendation. The timeseries forecasting techniques based on the metabolic grey model, autoregressive integrated moving average modelgrey model, and induced ordered weighted geometric averaging operator have been investigated in [4], which has found a way to provide reliable information and has indicated that the results from the timeseries and econometric forecasting technique are consistent. The above two studies have shown that the prediction of relevant models is effective and also indicate the research prospects and directions in this field. In addition, the singlelinear, hybridlinear, and nonlinear forecasting techniques based on grey theory have been presented in [5] to forecast energy demand in both China and India more accurately. In a novel study [6], researchers have found that the NMGMARIMA technique can significantly improve forecasting effectiveness and outperform other related forecasting models, which have been a guide to practical applications. In the field of primary energy consumption (PEC), many mathematical models have been successfully applied in prediction. For example, a new hybrid method (HAP) has been proposed for estimating energy demand in Turkey using particle swarm optimization (PSO) and ant colony optimization (ACO), and the research has found that different models have different advantages [7]. Two new models based on artificial bee colony (ABC) and particle swarm optimization (PSO) techniques have been proposed in [8] to estimate electrical energy demand in Turkey. In addition, the research methods in [9, 10] also have made good predictions. Among these models, the ARIMA model [11] and backpropagation (BP) neural networks [12] are the two most commonly used predictive models.
The forecast of energy consumption, especially the PEC, has always been a significant reference point for a territory [7]. Specifically, on the one hand, it is always important for a territory to make the future energy policies, to achieve a better development and to reduce greenhouse gas emissions [13]. In the recent study [14], particle swarm optimization (PSO) and artificial bee colony (ABC) techniques have been applied in estimating CO_{2} emission in Turkey based on socioeconomic indicators. On the other hand, the PEC refers to the overall energy consumption within a geographic territory and represents the total supply of energy available to the territory which supports all the requirements for energy transformation and final consumption in that territory. In general, the PEC of a territory includes both its indigenous energy sources and imported energy commodities consumed within the territory [15].
As for models, both economists and applied mathematicians have built a large number of models and accumulated experiences in forecasting energy [16]. For instance, the macroeconomic analysis of energy [17] and the analytic network process in energy policy planning [18] are two methods which are based on econometrics transforming the forecasting problem into an economic one and are seemingly in accordance with the thought of “big data.” However, applied mathematicians have developed serial linear or nonlinear mathematical models to complete the prediction [19].
Motivated by previous studies which show that the consumption of primary energy could be approximated by using linear or nonlinear econometrics with economic and/or noneconomic indicators, mathematical models [20] and especially the ARIMA model [21] and the BP neural network [22] are common models used for prediction. Therefore, in this paper, we will adopt the ARIMA model and the BP neural network one to make predictions and draw conclusions through making comparisons of models’ correctness.
This paper also focuses on solving some gaps in previous research and filling these gaps. Compared with the previous research concerning the single model [1], this paper adopts two common and effective mathematical models and then makes a comparison. In contrast to the research taking into account multiple models [2], this paper not only explores the accuracy of comparison but also opens up a new way to use neural networks to reveal the possible time of economic transformation.
The novelty and usefulness of this paper can be summed up as follows: (a) Compared with the available studies focusing on only one territory or one mathematical model which leads to the lack of comparison between different models and different types of territories, this study, however, includes three cities of different developing types from Greater Bay Area to make a comparison and then obtain more reasonable conclusions. It is shown in this paper that both the ARIMA model and the BP neural network model can be applied in different developing types of territories. (b) As for the novelty, this paper concludes that the singleinput BP neural network model can be used for judging whether the territory has transformed completely. (c) In terms of practicality, this paper explores that ARIMA can be applied to predict the PEC in addition to manufacturing territories, and the BP neural network model can be applied to predict the PEC in different types of territories.
2. Methodology
2.1. Models and Territories
It is shown in [23] that both the ARIMA model and the BP neural network model are suitable for timeseries data and can better catch on the related data and predict future points in the series. In this paper, we consider three cities including Hong Kong, Shunde, and Zhaoqing from Greater Bay Area in China as involved representative territories. The reason why we choose these territories is that they have varying development models; that is, they are three different industrydominated cities. Hong Kong has been a servicedominated city since the 2000s or earlier. Shunde has turned to be a manufacturingdominated city since the beginning of the Chinese economic reform. But Zhaoqing is famous for its agriculture, so it is regarded as an agriculturedominated one.
The ARIMA is a commonly used timeseries model which achieves the object’s characteristics of selfsimilarity, periodicity, suddenness, and trends [23] and has a better achievement in the shortterm subject’s forecast. Therefore, it has been applied in the prediction of the stock price index, the blood glucose concentrations [24], the current blockchain technology [25], the wind generation [26], and so on. The BP neural network is a commonly used timeseries and nonlinear prediction model applied in prediction of shortterm wind power, indoor temperature, wind speed [27], and hydraulic press machine.
In short, we deliberately consider the agricultural territory, the manufacturing territory, and the servicedominated territory combined with effective ARIMA and BP neural network models.
2.2. Datasets
The monthly data are adopted to be the time step of the data in this paper. Four monthly data series adopted in this paper are addressed in detail below.
The first one is the PEC which means the total energy consumed by a geographic territory. Generally, it represents the energy including the transformation parts, the final consumption ones [28], the energy produced locally, and the imported energy sources consumed locally. Therefore, to calculate the PEC, we take Hong Kong as an example; first, we need to sum up the renewable energy produced locally, net import of coal, oil products, and electricity and then minus the net usage of energy storage needed, and finally, we need to adjust the results by the supply from stock. This type of data comes from 2000 to 2017 for the agricultural territory, comes from 2005 to 2017 for the manufacturing one, and comes from 1979 to 2017 for the servicedominated one. Data are obtained from Statistical Yearbook of Shunde published by Development Planning and Statistics Department of Shunde, Statistical Yearbook of Guangdong Province published by Statistics Department of Guangdong, and Hong Kong Energy Statistics published by Census and Statistics Department of Hong Kong [28].
The second one and the third one are population and GDP. Population is used as a proxy for the demands of humans, and GDP can be regarded as a symbol of the growth situation of the economy. These two types of data are selected from 2000 to 2017 for the agricultural territory, from 2005 to 2017 for the manufacturing one, and from 1979 to 2017 for the servicedominated one. It should be noticed that the way we obtain these data is the same as the one we obtain PEC [29].
Last but not least, additional three different forecasting parameters are employed in this paper, which are the value added from primary products, the manufacturing industryadded value, and the values of total services. These values are of great importance for cities of different developing types. The value added from primary products indicates the growth of products directly from the natural sector (including plantation, forestry, animal husbandry, and fisheries) in this liquidation cycle compared to the previous liquidation cycle. It is shown in [29] how the territory gains from the agricultural products. In addition, the manufacturing industryadded value represents the final result of the production activities of the secondary industry in the form of money during the reporting period, which is the total result of all production activities of the production unit minus the value of the physical products and services consumed or transferred in the production process. Finally, the values of total services include the exports of services and imports of services, where the exports of services refer to the services that Hong Kong have sold to other entities and the imports of services mean the services purchased from the other entities. The data are obtained from Hong Kong Trade in Services Statistics published by Census and Statistics Department of Hong Kong [30].
2.3. ARIMA Modeling
Developed by Box and Jenkins, the ARIMA model is a linear regression model and is widely applied in forecasting when it is with timeseries data [31]. ARIMA(, d, q), an autoregressive integrated moving average model where parameters , d, and q are nonnegative integers, consists of three parts including the autoregressive (AR) model, differencing (I) model, and moving average (MA) model. Specially, represents the order of the AR model, d means the degree of I, and q is the order of the MA, respectively. To achieve the prediction through ARIMA(, d, q), these three steps must be completed, that is, to smooth the sequence through dorder difference calculation, to obtain a more suitable and d by calculating and comparing the autocorrelation coefficient with the partial autocorrelation one, and to carry out the prediction by the selected model. Mathematically, the expression of ARIMA(, d, q) is as follows [31]:where and are the actual value and random error at every year and and are two parameters of the ARIMA model. To be clear, and q are integers and the orders of the ARIMA model. By extracting the common factor, a new expression can be obtained as follows:
It can be expressed asby settingwhere and are the order and qorder characteristic polynomials of the independent variable “L.”
To conclude, in ARIMA(, d, q), d is the integral sum order, is the autoregressive coefficient, and q is the moving average coefficient [32].
2.4. BP Neural Network Modeling
Suppose that the network has R nodes [33] and the transfer function of each layer is of sigmoid type [34]. The following notations are used throughout this section: denotes the output of the ith hiddenlayer node; denotes the output of the kth outputlayer node; and represent the weight between node i and node j and the weight between node k and node i, respectively; denotes the threshold for hiddenlayer node i; and denotes the threshold for outputlayer node k [35].
The BP neural network can be understood in detail below [36].
Firstly, the sample forward propagation is input, which is a way of calculating the sample from the input layer, passing through the hidden layer, and finally to the output layer. The input samples should be known, and then the corresponding output is obtained by the above relevant steps.
In the hidden layer, the output of the ith neuron iswhere is the input value, is an integer denoting the number of neurons in the hidden layer, and is the activation function defined as follows:
In the output layer, the output of the kth neuron iswhere is an integer denoting the number of neurons in the output layer and is defined as follows:
Secondly, the error backpropagation is output.
Through the above forward propagation calculation, the actual output can be obtained. However, in general, the actual expectations are different. When two values are different or the errors of two values exceed some specified values, corresponding learning corrections should be made for the network.
In particular, the error function E is defined as follows:where is the predicted value.
In the output layer, the weights from the ith input to the kth output arewhere is the input value and can be described aswhere is the range fault.
After implicating the layer weight correction, the weight from the jth input to the ith output can be described as follows:where η is the learning coefficient.
The error of the output is reversed from the output layer, through the intermediate layer and finally to the input layer, where each layer is corrected once. The threshold can be seen as one of the weights and also contributes to the adjustment of the weights.
Thirdly, the memory training is made to be cyclic. In order to improve the accuracy of the network and reduce the appearance of the output error, it is necessary to carry on the loop memory training for all the samples input to the network. In addition, the number of loops should not be too few so that this sample mode can be effectively remembered by the network.
Last but not least, the end of the study is checked. The output error is checked to see if it meets the required standard whenever the loop memory training is completed. If it meets the requirements, the process is ended. Otherwise, the loop training is made again until it meets the requirements [37].
3. Results and Discussion
3.1. ARIMA Modeling
Following the above methodology, the PEC data of three different developing territories are used as the training dataset in the ARIMA model. In addition, unit root tests are stochastic trends in time series which are executed for the stationarity in time series. The stationarity in time series is that a shift in time does not cause a change in the shape of the distribution [38]. To obtain a smooth sequence and “” parameter of each model, we conduct three different unit root tests from zero orders to three orders for each model and make a comparison of three types of results. The unit root tests and results are shown in Tables 1–6.






From Tables 1–6, the augmented Dickey–Fuller (ADF) statistic is a statistical method to check whether a time series is stationary or not [37]. And the test critical value means threshold values, comparing the ADF statistic with the following three values at three significant levels. It is observed from the data from Tables 1–6 that the secondorder difference is suitable for both the agricultural territory’s model and the model of the manufacturing territory and the zeroorder difference is suitable for the one of the servicedominated territory. Therefore, we obtain the correlogram which shows the autocorrelation and the partial correlation with the help of EViews, where the autocorrelation is the correlation of a signal with its delayed copy and the partial correlation means the degree of association between two random variables [39]. Both the autocorrelation and the partial correlation are shown in Figures 1–3. In addition, the AC and PAC in Figures 1–3 represent the autocorrelation coefficient of the sequence and the partial correlation coefficient of the sequence, respectively. QStat represents the output of the statistical test about whether any of a group of autocorrelations of a time series are different from zero which obeys the chisquare distribution, and Prob represents its possibilities [40]. Representatively, it is presented in Figure 3 that the smoothness of the correlation data is good without difference. Furthermore, the autocorrelation function image is decreased, and the partial autocorrelation function image is truncated after the first image. Therefore, we can determine the specific values of the coefficients based on these facts. The specific determination method can be shown as follows.
The determination of the coefficients requires experience, but the theoretical support is still very meaningful. According to some available literatures like [41], the determination of the values of can be summed up as follows: (a) Since the sequence requires a firstorder difference to achieve smoothness, we determine that the value of d is 1. (b) Next, based on the fact that the autocorrelation images are truncated and the partial correlation images are smeared, we conclude that the sequence is suitable for the AR model. Among them, since none of the partial correlation images has been outside the confidence interval, we determine that the value of is 0. (c) We also determine that the value of q is 1 because the autocorrelation images are truncated. (d) Finally, we find out that the best model is obtained by taking a combination of several values of and q and simulating the prediction. To conclude, as it is shown in the correlogram, ARIMA(0,1,1) can be used for the prediction of the primary energy consumption of the servicedominated territory when both AC and PAC are considered. Similarly, it is seen that the applicability of ARIMA(2, 0, 2) and ARIMA(2, 0, 5) for the agricultural territory and the manufacturing territory can be guaranteed, respectively. Then, we construct three models, and the summaries of these three models are shown in Table 7. The stationary in Table 7 represents the measurement which compares the stationary part of the model with a simple mean model. A positive stationary means the model under consideration is better than the baseline one. In addition, the in Table 7 means the goodness of fit. The closer the value of to 1, the better the fit of the regression line to the observations [42].

It is noticed that the of ARIMA(2, 0, 2) is 0.941 and the corresponding one of ARIMA(0, 1, 1) is 0.955, which means that both of them are greater than 0.80. So the models are considered as the fitting ones for the agricultural and servicedominated territories [37]. However, ARIMA is not suitable for the manufacturing territory because the of ARIMA(2, 0, 5) is 0.721 which is less than 0.80. Finally, the results of predicting numbers of primary energy consumption of three territories are shown in Figure 4, where UCL and LCL denote the upper control line and the lower control line, respectively. Their mathematical definitions are as follows:where μ is the average value of the output and σ is the standard deviation of the output.
(a)
(b)
(c)
3.2. BP Neural Network Modeling
A threelayer backpropagation neural network with three hidden layers is developed for predicting three territories’ primary energy consumption. Based on Kolmogorov’s theorem [43], for any continuous function, it can be implemented with a threelayer network, where the input layer has m units, the hidden layer has nodes, and the output layer has n units. So the number of hidden nodes is set to be three. And the following results show that our models are feasible under this condition. Two activation functions, i.e., the pure linear function and the logsigmoid transfer function, are considered in the hidden layers. With the help of MATLAB software, we establish two types of models including three inputs and single input. The singleinput model means only one among four predictors (population/primary GDPs/manufacturing industryadded value/values of total services) is adopted. The singleinput model can be used to observe whether one variable has a significant impact on the PEC and can also be used to determine whether the prediction can be completed with only one variable. In contrast, the threeinput model includes three predictors in the model. The schematic diagram is presented in Figure 5.
After debugging repetitions, we set the number of training iterations to be 50000. It should be noticed that the number of hiddenlayer nodes can influence the network structure to some extent. On the one hand, if the number of hiddenlayer nodes is too large, an overmatching phenomenon will occur. On the other hand, if the number of hiddenlayer nodes is too small, the useful information obtained by the network from the original input may be scanty, which is not enough to discover the characteristics of the data, and the generalized nonlinear learning ability of the model may be weaker. Therefore, according to Kolmogorov’s theorem, for any continuous function, it can be implemented with a threelayer network, where the input layer has m units, the hidden layer has nodes, and the output layer has n units [44]. So the number of hidden nodes is set to be 3. As mentioned above, we selected the relevant data from 1999 to 2014 as the training sample for the servicedominated territory and for other territories similarly. Different models’ predictions are presented in Figures 6–9.
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
We can see that the predicting data of the threeinput BP neural network model are more close to the actual data compared to the three other singleinput models. Therefore, we consider the threeinput model as the best BP neural network for prediction. Furthermore, we think some values have a significant impact on the PEC of both the agricultural and manufacturing territories, which means that the prediction of the PEC can be conducted with only one value. However, we see from Figure 8 that the singleinput model is not suitable for the servicedominated territory. The regression of the threeinput models is presented in Figure 10. As we can see in Figure 10, the regression line fits the observations very well, and the values of are very close to 1, which indicates that the three models for three different developing types of cities perform well.
(a)
(b)
(c)
3.3. Comparison of Models
We examine two forecasting models in three different developing types of cities in this paper. From Figures 4, 9, and 10, we draw an elementary conclusion that the forecasting capability of BP neural network models is better than that of ARIMA models on the whole.
In order to verify accurately the correctness, we make a comparison of these two types of models through calculating the mean absolute percentage error (MAPE) and the root mean square error (RMSE). The mathematical definition of the MAPE and RMSE is as follows:where is the accurate value of the data, is the predictive value of the data, and n is the number of data we employed. The result of the correctness is shown in Table 8.

4. Conclusion
In this paper, the autoregressive integrated moving average (ARIMA) model and the backpropagation (BP) neural network model are applied to forecast the primary energy requirement of three different developing types of territories, and a comparison of the accuracy is made.
Three different models including ARIMA(2, 0, 2), ARIMA(2, 0, 5), and ARIMA(0, 1, 1) are conducted to make predictions. Their accuracies are within the range from 0.721 to 0.941. We believe that the ARIMA model is not suitable for the manufacturing territories since the of ARIMA(2, 0, 5) is 0.721 which is less than 0.80. Furthermore, we conclude that ARIMA can predict accurately without employing too much data.
As for the BP neural network model, the data concerning primary products, manufacturing industryadded value, and values of total services are employed as the predictor in the model, and the most accurate prediction is made. In addition, we also use only one predictor to build three singeinput models, and the corresponding result shows that not only the accuracy of the singleinput models is acceptable but also the singleinput models can demonstrate whether the territory has been totally supported by certain industry. This conclusion is supported by economics. In the economic literature [45], experts have pointed out that Hong Kong has completely transformed into a servicedominated territory since 2009. And in the literature [46], Shunde was regarded as a manufacturingdominated territory in 2005. The above results support the conclusions of this paper. As for limitations, we believe that if researchers aim to obtain more accurate predictions, they need to collect multiple types of and a large amount of data, preferably monthly data, which will make the research work cumbersome. As for policymakers, this paper presents a judgment to determine whether the territory is fully transformed. When the territory has completely transformed into a specific development model, the corresponding energy policy should be more inclined to the industry.
In the process of this research, we also find that some new technologies can be well applied in predicting energy.In particular, the following papers deserve further research in the future: the literature [8],in which swarm intelligence approaches have achieved excellent predictions and comparable advantages over artificial neural networks (ANNs), and the literature[14], in which swarm intelligence approaches have also been applied in predicting greenhouse gas emissions, which is very important for sustainable development.
Data Availability
Previously reported data (four monthly data series) were used to support this study, as addressed in Section 2.2. The prior studies (and datasets) are cited at relevant places as references [18, 28, 30].
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Natural Science Foundation of China (Grant no. 11801077) and the Natural Science Foundation of Guangdong Province (Grant nos. 2018A1660005 and 2018A030310357).
References
 C. R. Rivero, V. Sauchelli, H. D. Patio, J. A. Pucheta, and S. Laboret, “Longterm power consumption demand prediction: a comparison of energy associated and bayesian modeling approach,” in Proceedings of the 2015 Latin America Congress on Computational Intelligence (LACCI), pp. 1–6, Curitiba, Brazil, October 2015. View at: Publisher Site  Google Scholar
 Y. He, J. Y. Wu, Y. Ge, D. Z. Li, and H. G. Yan, “Research on prediction of energy demand in china based on grey theory and system dynamics,” in Proceedings of the 2017 IEEE Conference on Energy Internet and Energy System Integration (EI2), pp. 1–6, Beijing, China, November 2017. View at: Publisher Site  Google Scholar
 Q. Wang, S. Y. Li, and R. R. Li, “China’s dependency on foreign oil will exceed 80% by 2030: developing a novel NMGMARIMA to forecast China’s foreign oil dependence from two dimensions,” Energy, vol. 163, pp. 151–167, 2018. View at: Publisher Site  Google Scholar
 Q. Wang, S. Li, and R. Li, “Will Trump’s coal revival plan work?—comparison of results based on the optimal combined forecasting technique and an extended IPAT forecasting technique,” Energy, vol. 169, pp. 762–775, 2019. View at: Publisher Site  Google Scholar
 Q. Wang, S. Li, and R. Li, “Forecasting energy demand in China and India: using singlelinear, hybridlinear, and nonlinear time series forecast techniques,” Energy, vol. 161, pp. 821–831, 2018. View at: Publisher Site  Google Scholar
 Q. Wang, X. Song, and R. Li, “A novel hybridization of nonlinear grey model and linear ARIMA residual correction for forecasting U.S. shale oil production,” Energy, vol. 165, pp. 1320–1331, 2018. View at: Publisher Site  Google Scholar
 M. S. Kiran, E. Özceylan, M. Gündüz, and T. Paksoy, “A novel hybrid approach based on particle swarm optimization and ant colony algorithm to forecast energy demand of Turkey,” Energy Conversion and Management, vol. 53, no. 1, pp. 75–83, 2012. View at: Publisher Site  Google Scholar
 M. S. Kiran, E. Özceylan, M. Gündüz, and T. Paksoy, “Swarm intelligence approaches to estimate electricity energy demand in Turkey,” KnowledgeBased Systems, vol. 36, pp. 93–103, 2012. View at: Publisher Site  Google Scholar
 X. Wang, D. Luo, X. Zhao, and Z. Sun, “Estimates of energy consumption in China using a selfadaptive multiverse optimizerbased support vector machine with rolling crossvalidation,” Energy, vol. 152, pp. 539–548, 2018. View at: Publisher Site  Google Scholar
 D. V. D. Meer, G. R. C. Mouli, G. M. E. Mouli, L. R. Elizondo, and P. Bauer, “Energy management system with PV power forecast to optimally charge EVs at the workplace,” IEEE Transactions on Industrial Informatics, vol. 14, no. 1, pp. 311–320, 2018. View at: Publisher Site  Google Scholar
 S. L Lai, M. Liu, K. C. Kuo, and R. Chang, “Energy consumption forecasting in Hong Kong using ARIMA and artificial neural networks models,” Applied Mechanics and Materials, vol. 672674, pp. 2085–2097, 2014. View at: Publisher Site  Google Scholar
 C. Wan, J. Zhao, Y. Song, Z. Xu, J. Lin, and Z. Hu, “Photovoltaic and solar power forecasting for smart grid energy management,” CSEE Journal of Power and Energy Systems, vol. 1, no. 4, pp. 38–46, 2015. View at: Publisher Site  Google Scholar
 T. Ma, P. A. Østergaard, H. Lund, H. Yang, and L. Lu, “An energy system model for Hong Kong in 2020,” Energy, vol. 68, pp. 301–310, 2014. View at: Publisher Site  Google Scholar
 E. Özceylan, “Forecasting CO_{2} emission of Turkey: swarm intelligence approaches,” International Journal of Global Warming, vol. 9, no. 3, pp. 337–361, 2016. View at: Publisher Site  Google Scholar
 Hong Kong Energy Statistics, Census and Statistics Department Hong Kong Special Administrative Region, Hong Kong Energy Statistics, Hong Kong, China, 2017.
 F. Jiang, X. Yang, and S. Li, “Comparison of forecasting India’s energy demand using an MGM, ARIMA model, MGMARIMA model, and BP neural network model,” Sustainability, vol. 10, no. 2, p. 2225, 2018. View at: Publisher Site  Google Scholar
 J. Xia and W. Wu, “Application of the gray neural network combination model in forecasting energy demand of China,” Computer Science and Application, vol. 8, no. 2, pp. 181–189, 2018. View at: Publisher Site  Google Scholar
 R. P. Hämäläinen and T. O. Seppäläinen, “The analytic network process in energy policy planning,” SocioEconomic Planning Sciences, vol. 20, no. 6, pp. 399–405, 1986. View at: Publisher Site  Google Scholar
 X. F. Song, G. Liang, C. Z. Li, and W. W. Chen, “Electricity consumption prediction for Xinjiang electric energy replacement,” Mathematical Problems in Engineering, vol. 2019, Article ID 3262591, 11 pages, 2019. View at: Publisher Site  Google Scholar
 H. R. Cui, R. R. Wu, and T. Zhao, “Dynamic decomposition analysis and forecasting of energy consumption in shanxi Province based on VAR and GM (1, 1) models,” Mathematical Problems in Engineering, vol. 2018, Article ID 4234206, 11 pages, 2018. View at: Publisher Site  Google Scholar
 S. A. Javed and S. Liu, “Predicting the research output/growth of selected countries: application of even GM (1, 1) and NDGM models,” Scientometrics, vol. 115, no. 1, pp. 395–413, 2018. View at: Publisher Site  Google Scholar
 J. Y. Liu, H. X. Wang, Y. Y. Sun, C. Q. Fu, and J. Guo, “Realcoded quantuminspired genetic algorithmbased BP neural network algorithm,” Mathematical Problems in Engineering, vol. 2015, Article ID 571295, 10 pages, 2015. View at: Publisher Site  Google Scholar
 E. Egrioglu, U. Yolcu, C. H. Aladag, and C. Kocak, “An arma type fuzzy time series forecasting method based on particle swarm optimization,” Mathematical Problems in Engineering, vol. 2013, Article ID 935815, 12 pages, 2013. View at: Publisher Site  Google Scholar
 J. Yang, L. Li, Y. Shi, and X. Xie, “An ARIMA model with adaptive orders for predicting blood glucose concentrations and hypoglycemia,” IEEE Journal of Biomedical and Health Informatics, vol. 23, no. 3, pp. 1251–1260, 2019. View at: Publisher Site  Google Scholar
 B. Wang, X. Zhu, Q. S. He, and G. D. Gu, “The forecast on the customers of the member point platform built on the blockchain technology by arima and lstm,” in Proceedings of the 2018 IEEE 3rd International Conference on Cloud Computing and Big Data Analysis (ICCCBDA), Chengdu, China, April 2018. View at: Publisher Site  Google Scholar
 R. S. F. Ferraz, R. S. F. Ferraz, F. C. Cruz, E. F. S. Filho, and A. F. S. Correia, “Multi step forecasting of the wind generation by ARIMA and KNN,” in Proceedings of the 2018 Simposio Brasileiro de Sistemas Eletricos (SBSE), Niteroi, Brazil, May 2018. View at: Publisher Site  Google Scholar
 J. Li, S. W. Yin, G. S. Shi, and L. Wang, “Optimization of indoor thermal comfort parameters with the adaptive networkbased fuzzy inference system and particle swarm optimization algorithm,” Mathematical Problems in Engineering, vol. 2017, Article ID 3075432, 13 pages, 2017. View at: Publisher Site  Google Scholar
 Hong Kong Energy Statistics, Census and Statistics Department Hong Kong Special Administrative Region, Hong Kong Energy Statistics, Hong Kong, China, 1979.
 B. S. Liang, B. Gao, and F. Liu, “Research on the dynamic relationship between fiscal expenditure on supporting agriculture and development of the primary industry in Hunan Province,” Journal of Hunan Agricultural University, vol. 43, pp. 450–458, 2009, in Chinese. View at: Google Scholar
 Hong Kong Trade in Services Statistics, Census and Statistics Department Hong Kong Special Administrative Region, Hong Kong Trade in Services Statistics, Hong Kong, China, 1999.
 A. Jeyasekar, S. V. Kasmir, and R. A. Uthra, “Congestion avoidance algorithm using ARIMA(2,1,1) modelbased RTT estimation and RSS in heterogeneous wiredwireless network,” Journal of Network and Computer Applications, vol. 93, no. 1, pp. 91–109, 2017. View at: Publisher Site  Google Scholar
 A. Rath, S. Samantaray, and P. P. C. Swain, “Flow forecasting of hirakud reservoir with ARIMA model,” in Proceedings of the 2017 International Conference on Energy, Communication, Data Analytics and Soft Computing (ICECDS), Chennai, India, August 2017. View at: Publisher Site  Google Scholar
 C. Ren, N. An, J. Wang, L. Li, B. Hu, and D. Shang, “Optimal Parameters Selection for BP neural network based on particle swarm optimization: a case Study of wind speed forecasting,” KnowledgeBased Systems, vol. 56, pp. 226–239, 2014. View at: Publisher Site  Google Scholar
 B. Y. Xu, H. J. Zhang, Z. T. Wang, H. X. Wang, and Y. L. Zhang, “Model and algorithm of BP neural network based on expanded multichain quantum optimization,” Mathematical Problems in Engineering, vol. 2015, Article ID 362150, 11 pages, 2015. View at: Publisher Site  Google Scholar
 H. Z. Li and S. Guo, “External economies evaluation of wind power engineering project based on analytic hierarchy process and matterelement extension model,” Mathematical Problems in Engineering, vol. 2013, Article ID 848901, 11 pages, 2013. View at: Publisher Site  Google Scholar
 P. Zhou, M. Yuan, H. Wang, and T. Y. Chai, “Datadriven dynamic modeling for prediction of molten iron silicon content using ELM with selffeedback,” Mathematical Problems in Engineering, vol. 2015, Article ID 326160, 11 pages, 2015. View at: Publisher Site  Google Scholar
 P. G. Li, J. Zhang, C. B. Li et al., “Dynamic similar subseries selection method for time series forecasting,” IEEE Access, vol. 6, pp. 32532–32542, 2018. View at: Publisher Site  Google Scholar
 A. Bhargava, “On the theory of testing for unit roots in observed time series,” The Review of Economic Studies, vol. 53, no. 3, pp. 369–384, 1986. View at: Publisher Site  Google Scholar
 S. Vergura and M. Carpentieri, “Phase coherence index, HHT and wavelet analysis to extract features from active and passive distribution networks,” Applied Sciences, vol. 8, no. 1, p. 71, 2018. View at: Publisher Site  Google Scholar
 C. Guarnaccia, J. Quartieri, E. R. Rodrigues, and C. Tepedino, “Acoustical noise analysis and prediction by means of multiple seasonality time series model,” International Journal of Mathematical Models and Method in Applied Sciences, vol. 8, pp. 384–393, 2014. View at: Google Scholar
 J. Liu, G. Y. Shi, and K. G. Zhu, “Highprecision combined tidal forecasting model,” Algorithms, vol. 12, no. 3, p. 65, 2019. View at: Publisher Site  Google Scholar
 Z. S. Wu and R. Xue, “A cyclical nonlinear inertiaweighted teachinglearningbased optimization algorithm,” Algorithms, vol. 12, no. 5, p. 94, 2019. View at: Publisher Site  Google Scholar
 A. N. Kolmogorov, “On tables of random numbers,” Theoretical Computer Science, vol. 207, no. 2, pp. 387–395, 1998. View at: Publisher Site  Google Scholar
 M. Paluch and L. JackowskaStrumillo, “Hybrid models combining technical and fractal analysis with ANN for shortterm prediction of close values on the warsaw stock exchange,” Applied Sciences, vol. 8, no. 12, p. 2473, 2018. View at: Publisher Site  Google Scholar
 Y. Li, Returning to 20 Years: Review and Prospect of Hong Kong’s Economic Development, Chinese Finance, Beijing, China, 2017.
 R. N. Su, “Study on the Path of Division of Labor, cooperation and competition in the international value chain in the development sstrategy of “Guangzhoufoshanonecity,” Special Economic Zone, vol. 3, pp. 152–161, 2017. View at: Google Scholar
Copyright
Copyright © 2019 NingKang Pan and Chunwan Lv. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.