Application of the novel fractional grey model FAGMO(1,1,k) to predict China's nuclear energy consumption

by   Wenqing Wu, et al.

At present, the energy structure of China is shifting towards cleaner and lower amounts of carbon fuel, driven by environmental needs and technological advances. Nuclear energy, which is one of the major low-carbon resources, plays a key role in China's clean energy development. To formulate appropriate energy policies, it is necessary to conduct reliable forecasts. This paper discusses the nuclear energy consumption of China by means of a novel fractional grey model FAGMO(1,1,k). The fractional accumulated generating matrix is introduced to analyse the fractional grey model properties. Thereafter, the modelling procedures of the FAGMO(1,1,k) are presented in detail, along with the transforms of its optimal parameters. A stochastic testing scheme is provided to validate the accuracy and properties of the optimal parameters of the FAGMO(1,1,k). Finally, this model is used to forecast China's nuclear energy consumption and the results demonstrate that the FAGMO(1,1,k) model provides accurate prediction, outperforming other grey models.



There are no comments yet.


page 7


Future Energy Consumption Prediction Based on Grey Forecast Model

We use grey forecast model to predict the future energy consumption of f...

On the Energy Consumption Forecasting of Data Centers Based on Weather Conditions: Remote Sensing and Machine Learning Approach

The energy consumption of Data Centers (DCs) is a very important figure ...

A Complete Model for Modular Simulation of Data Centre Power Load

Data centres are very fast growing structures with significant contribut...

Evaluation of Deep Learning-based prediction models in Microgrids

It is crucial today that economies harness renewable energies and integr...

Optimal Stochastic Evasive Maneuvers Using the Schrodinger's Equation

In this paper, preys with stochastic evasion policies are considered. Th...

Prediction of Energy Consumption for Variable Customer Portfolios Including Aleatoric Uncertainty Estimation

Using hourly energy consumption data recorded by smart meters, retailers...

IrEne: Interpretable Energy Prediction for Transformers

Existing software-based energy measurements of NLP models are not accura...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

Energy is the most important strategic resource and provides a key material basis for economic development and social progress. Energy consumption prediction constitutes an important aspect of energy policies for countries globally, particularly developing countries such as China, where the energy consumption structure is changing at a rapid speed. Numerous models have been introduced for forecasting energy consumption, such as dynamic causality analysis Mirza2017EnergyRSER, nonlinear and asymmetric analysis Shahbaz2017EnergyEE, time-series analysis Bekhet2017CO2RSEE; Amri2017TheRSEE

, machine learning models

fan2018comparison, the coupling mathematical model Wang2017TransientJPM; Wang2018FlowIJNSNS; Hu2018streamlineSPEPO, autoregressive distributed lag model Brini2017RenewableRSEE, hybrid forecasting system Du2018MultistepRE; Cai2018developmentE, machining system Cai2018EnergyE, fuzzy systems Wu2017TopologicalIJBC, LEAP model schnaars1987houLRRP; Dong2017APS, TIMES model Shi2016modellingAE; Zhang2016timesAE, NEMS model Gabriel2001theOR; Soroush2017aIJSSD and grey model zeng2018forecasting; zeng2017self; Feng2012ForecastingESPBEPP; Cui2013AAPM; duan2018forecasting; Zeng2018ForecastingE; Zeng2018improvedCIE; Wu2018usingE; Wang2018AnASC

. Among these prevalent methods, simple linear regression, multivariate linear regression, and time-series analysis are often significant in accurately demonstrating the phenomena of long-term trends. However, these exhibit the limitations of requiring a large amount of observed data, at least 50 or more sets, to construct models. The computational intelligence method requires a substantial amount of training data to derive the optimised parameters. However, in many practical situations, it is very difficult and sometimes even impossible to obtain complete information. Therefore, it is important to identify a favourable method for forecasting the trend of an analysed system using scarce information with less errors.

The grey forecasting theory, proposed by Professor Deng Deng1982ControlSCL, offers a feasible and efficient method for dealing with uncertain problems containing poor information. The main advantage of this theory is that only four or more samples are required to describe the behaviour and evolution of the analysed system. In Deng’s pioneering work, the first-order one variable grey model GM(1,1) was discussed in detail. Over three decades of development, the classical continuous GM(1,1) model has been studied extensively; for example, by Xie et al. Xie2009DiscreteAMM; Xie2013OnAPMAMM; Xie2017AJGS, Wang et al. Wang2017ForecastingJCP; Wang2017GreyAMM; Wang2017DecompositionEP, Ma et al. Ma2017ApplicationJCAAM; Ma2017ACNSNS; Ma2018ThekernelAMM; Ma2018PredictingNCA and others. However, we note that these generalised grey models all include integer-order accumulation, which results in less flexibility in time-series forecasting. Thus, the fractional-order accumulation grey model is considered in this paper.

By extending the integer accumulated generating operation into the fractional accumulated generating operation, Wu et al. Wu2013GreyCNSNS first proposed the fractional accumulation GM(1,1) model known as the FAGM(1,1) model. The computational results demonstrated that the novel model outperformed the conventional GM(1,1) model. Later, Wu and his peers successfully applied fractional accumulation to the fuel production of China Wu2014NonNCA, tourism demand Wu2015UsingSC and electricity consumption Wu2015PropertiesAMAC. Subsequently, Xiao et al. Xiao2014TheAMM studied the GM(1,1) model, in which they regarded the fractional accumulated generator matrix as a type of generalised accumulated generating operation. Gao et al. Gao2015EstimationJGS presented a new discrete fractional accumulation GM(1,1) model known as FAGM(1,1,D) and applied it to China’s CO emissions. Mao et al. Mao2016AAMM investigated a novel fractional grey model FGM(,1). Interested readers can refer to Shen2014OptimizationJGS; Yang2016ContinuousE; Yang2017ModifiedJAS; Liu2017NonJIM for further details on fractional accumulation grey models.

A further significant issue in grey system theory is that the solution applied for prediction does not match the grey difference equation. In 2009, Kong and Wei Kong2009OptimizationJGS proposed a parameter optimisation technique to study the DGM(2,1) model. Later, Chen et al. used a similar technique to improve the GM(1,1) Yu2017ApplyingJGS and ONGM(1,1) models Chen2014FoundationMPE, in which the basic structure of the original models remain in the optimised ones. Recently, Ma and Liu Ma2016PredictingJCTN studied the exact non-homogeneous grey prediction model (ENGM) with an exact basic equation and background value. Thereafter, Ma and Liu Ma2017TheJGS considered the GMC(1,) model with optimised parameters and applied it to forecasting the urban consumption per capita and industrial power consumption of China. Following the concept of fractional accumulation and the parameter optimisation method, we propose a novel FAGMO(1,1,) model.

In this paper, we study the nuclear energy consumption of China by means of the FAGMO(1,1,) model. The computational results indicate that the proposed grey model outperforms the existing ENGM model, optimised non-homogeneous grey model abbreviated as the ONGM(1,1,) model, FAGM(1,1) model and FAGM(1,1,) model. The main contributions of our paper are listed below. 1) A fractional accumulation grey model with optimised parameters is developed. 2) Detailed properties of optimised parameters are studied according to two theorems. These indicate that the first parameter is the most important factor affecting the accuracy of the FAGM(1,1,) model. 3) Simulation results and two practical cases are considered to assess the effectiveness of the FAGMO(1,1,) model compared to other models. 4) The FAGMO(1,1,) grey forecasting model is implemented to forecast the nuclear energy consumption of China. It is demonstrated in the results that the newly proposed model offers higher precision than other grey models.

The remainder of this paper is organised as follows. Section 2 provides a compendium of China’s energy consumption. Section 3 discusses several preliminaries. A detailed discussion of the FAGM(1,1,) model is provided in section 4. Section 5 discusses the optimised parameters. Modelling evaluation criteria and detailed steps are provided in section 6. Section 7 discusses the validation of the FAGMO(1,1,) model. Applications are explained in section 8 and conclusions are drawn in the final section.

2 Brief overview of China’s energy consumption

This section presents a systematic and comprehensive investigation of China’s energy consumption using five fuels, namely coal, oil, natural gas, nuclear energy and renewables. In China, renewables include hydroelectricity, wind, solar, geothermal, biomass and others. According to the statistical data of British Petroleum (BP) Statistical Review of World Energy 2018 ( statisticalreview), the International Energy Agency (IEA) World Energy Outlook 2017 (, Asia-Pacific Economic Cooperation (APEC) Energy Overview 2017 (, and National Bureau of Statistics of China (NBS) China Statistical Yearbook 2017 (, China’s primary energy consumption increased from 142.9 million tonnes oil equivalent (Mtoe) in the first year of the third Five-Year Plan (1966 to 1970) to 3132.2 Mtoe in the second year of the Five-Year Plan (2016 to 2020), and increased dramatically since the turn of the millennium owing to continuous economic growth. According to the statistical data of BP, China’s primary energy consumption from 1966 to 2017 is plotted in Fig. 1.

Figure 1: Total primary energy consumption of China from 1966 to 2017

It is well known that China is the world’s largest energy consumer, accounting for 23% and 23.2% of the global energy consumption in 2016 and 2017, respectively. While coal remains the dominant fuel, its share of total energy consumption was 62% in 2016 and 60.4% in 2017. China’s Five-Year Plan set an ambitious target for adjusting the primary energy consumption structure. The energy plan set by China for the Five-Year Plan can met the adjustment target of the primary energy consumption structure. A brief overview of China’s primary energy consumption from the perspective of five fuel types is provided below.

2.1 Coal

Since the foundation of the People’s Republic of China, coal has always been the primary energy fuel, owing to abundant domestic reserves and its low cost Dong2017APS. From Figs. 2 and 3, it can be observed that coal soared from 122.4 Mtoe in 1966 to 1892.6 Mtoe in 2017, although the percentage of coal in the total primary energy consumption decreased from 85.7% in 1966 to 60.4% in 2017. Specifically, despite a continuous increase in coal consumption during the third and fourth Five-Year Plan periods, the proportion of coal in the total primary energy consumption has decreased from 85.7% to 72.5%. Following China’s Reform and Opening-Up Policy in 1978, the coal consumption has expanded rapidly from 282.8 Mtoe in 1978 to 1685.8 Mtoe in 2010, while the share of coal in the total primary energy consumption is around at 73.8%. At the beginning of the Five-Year Plan (2011 to 2015), the Chinese government has been stepping up its efforts to reduce coal consumption to deal with air pollution and climate change. The “supply-side reform” removes unnecessary and out-dated production capacity to avoid supply overcapacity in the coal mining industry. During this period, the consumption of coal decreased from 1903.9 Mtoe in 2011 to 1892.6 Mtoe in 2017, while the share of coal decreased from 70.8% to 60.4%.

2.2 Oil

Oil is a major component of primary energy resources globally and plays a strategic role in economic growth. It can be observed in Figs. 2 and 3 that the oil consumption in China increased from 14.3 Mtoe in 1966 to 608.4 Mtoe in 2017, with an average annual growth rate of 7.5%. Owing to China’s oil reserves accounting for only 2% of the global amount, China is highly dependent on overseas oil imports of more than 60%. China became a net importer of crude oil in 1993 and the world’s second largest oil consumer in 2002 Dong2017APS. In April 2015, China surpassed the US as the world’s largest oil importer, with imports of 7.4 million barrels per day (Mbbl/D), thereby exceeding the US imports of 7.2 Mbbl/D. Following the Five-Year Plan, the oil consumption increased rapidly owing to economic growth and the improved quality of life.

Figure 2: China’s primary energy consumption under fuel types
Figure 3: Percentage of China’s primary energy consumption under fuel types

2.3 Natural gas

Natural gas is a fossil fuel for electricity generation, chemical feedstock, heating and cooking, among others. Chinese organisations have estimated that the technically and ultimately recoverable resources of natural gas are 6.1 trillion cubic meters (tcm) and 37 tcm

Hou2015unconventionalEES, respectively. However, natural gas has not become a major energy resource in China because the domestic natural gas industry has developed slowly. In recent years, the Chinese government has set the stable natural gas supply as one of the country’s energy strategies and encourages gas transportation from areas with significant resources to East China. The National Development and Reform Commission constructed three west-east gas pipelines in 2004, 2007 and 2015, respectively. Furthermore, the “shifting from coal to gas” policy has a significant impact on the natural gas market. The natural gas consumption has increased from 116.2 Mtoe in 2011 to 206.7 Mtoe in 2017, with an average annual growth rate of 8.6%.

2.4 Nuclear energy

Nuclear energy is almost always used to generate electricity. To reduce the air pollution from coal-fired power plants, nuclear energy is an inevitable strategic option for China. In fact, China began to develop nuclear energy in the 1980s and the Qinshan Nuclear Power Plant began operating in 1991. In 2012, the State Council set a goal of 58 GW nuclear capacity by 2020. At the beginning of the Five-Year Plan, 38 nuclear power reactors were in operation with a production of 213.3 TWh, while 19 nuclear power reactors were under construction. At present, the Chinese government focuses on fourth-generation reactors with increased safety. From the Five-Year Plan (2006 to2010), nuclear energy consumption has soared rapidly from 12.4 Mtoe in 2006 to 56.2 Mtoe in 2017, with an average annual growth rate of 13.4%.

2.5 Renewables

China’s renewable energy has been expanding rapidly in recent decades, owing to the development of the modern renewable energy industry. In 2017, China’s renewables consumption accounted for 21.9% of the total global amount, increasing by 31% and accounting for 36% of the global renewables consumption growth. Meanwhile, the renewables consumption increased from 101.1 Mtoe in 2006 to 368.3 Mtoe in 2017; the share has increased from 5.1% to 11.8% with an average annual growth rate of 11.4%. The Five-Year Plan set targets for an installed wind power generation capacity of 250 GW, solar power generation capacity of 110 GW, and hydropower generation capacity of 350 GW by 2020.

In summary, China’s primary energy consumption using five fuels for the period of 1966 to 2017 can be provided below. The coal consumption has gradually declined, the oil consumption has gradually increased, and the natural gas, nuclear energy and renewables have rapidly increased. China’s primary energy consumption structure exhibits a diversified trend, and the clean energy has increased yearly.

3 Definitions and properties of fractional accumulation

This section provides the fractional accumulated generating operation (AGO), which can reduce the randomness of raw data in grey theory. Correspondingly, the inverse operation of accumulated generation is known as the inverse accumulated generating operation (IAGO). The AGO and IAGO are provided below, which can be found in paper Wu2013GreyCNSNS; Mao2016AAMM.

Definition 1

Let be an original sequence and be the accumulated generating operation (-AGO) sequence of , where . Denote by the -AGO matrix that satisfies , and


Obviously, the 1-AGO sequence , namely with .

Definition 2

The inverse accumulated generation is defined as . Denote by the inverse accumulated generating operation (-IAGO) matrix, which satisfies , and


Similarly, the 1-IAGO sequence ; that is, , with .

Theorem 1

The expression is a function of and ; for any value ,

is a monotonically decreasing function of ;

; and

is a monotonically increasing function of .

Proof 1

We consider the difference

From the difference results, we complete the proof.

To gain an improved understanding of Theorem 1, two figures are displayed in the following Fig. 4

Figure 4: Function versus values and : left , right

It follows from that


which means that is the weight of .

From Theorem 1, when , the weight of the old data is smaller than that of the new data. When , the weights of the old and new data are all 1. When , the weight of the old data is larger than that of the new data.

Theorem 2

The values of -AGO and -IAGO satisfy .

Proof 2

From the definition of , it is easy to calculate the determinant , which means that is reversible.

Employing mathematical induction, when , we obtain

Assuming that the properties hold true when , this means that

Then, when , we obtain

so the result is proven.

4 Fractional grey FAGM(1,1,) model

Definition 3

The first-order differential equation


is known as the whitening differential equation of the FAGM(1,1,) model. The parameter is a development coefficient, while is the grey action quantity.

The discrete differential equation


is referred to as the basic equation of the FAGM(1,1,). , .

The least-squares estimation for of the FAGM(1,1,) model satisfies



in which is the number of samples used to construct the model.

Theorem 3

The time response function of the FAGM(1,1,) model is


and the restored value of can be expressed by

Proof 3

From Eq. (2), we have


Let ; then, Eq. (7) is transformed into


To perform the indefinite integral on Eq. (8) and reduce it, we obtain


where is a constant to be determined.

Substituting and into Eq. (9), we obtain


It follows from Eqs. (9) and (10) that


Thus, the time response function of the FAGM(1,1,) model is

and the restored value of can be expressed by

Setting in Eq. (2), the fractional FAGM(1,1,) model is reduced to the fractional FAGM(1,1) model Wu2013GreyCNSNS with the form


Setting in Eq. (2), the fractional FAGM(1,1,) model is reduced to the GM(1,1,,) model Chen2014FoundationMPE with the form


Setting , in Eq. (2), the fractional FAGM(1,1,) model is reduced to the GM(1,1,) model Cui2013AAPM with the form


Setting , , in Eq. (2), the fractional FAGM(1,1,) model is reduced to the GM(1,1) model Deng1982ControlSCL with the form


Thereafter, the flaw of the FAGM(1,1,) model is provided. Integrating both sides of Eq. (2) in the interval , we obtain


With the knowledge of , and , the exact discrete differential equation is expressed by


A comparison between Eq. (17) and the basic Eq. (3) indicates that differences exist in the background value and . It is highly inaccurate to compute the integration utilising the trapezoid formula if is not a linear function. Thus, the basic form and whitenisation differential equation of the FAGM(1,1,) model do not strictly match.

5 Parameter optimisation of FAGM(1,1,) model

It can easily be verified that the parameters , derived by the least-squares estimation in Eq. (3) and the parameters of the time response function derived by Eq. (2), have different meanings. When the response function dose not satisfy the basic equation, large errors may arise. To match the basic Eq. (3) and response function (5), the system parameters are optimised in this system.

Setting the optimised parameters of the grey system as and replacing the parameters in Eq. (2), the whitening differential equation is rewritten as


Similarly, the general solution of Eq. (18) is given by


Furthermore, we have


Substituting Eq. (20) into the left side of Eq. (3), we obtain


Owing to the left side and right side equivalence, namely , it is implied that


It follows from Eqs. (22) to (24) that


Thus, the optimised parameters are obtained by Eqs. (25) to (27), and they also indicate that the parameters derived by the least-squares estimation satisfy the relationship in Eqs. (25) to (27). In this paper, the FAGM(1,1,) model with optimised parameters is referred to as the FAGMO(1,1,) model.

Theorem 4

Assuming that the original data satisfy Eq. (20) with the given parameters , the parameters of FAGMO(1,1,) obtained by Eqs. (4) and (25) to (27) satisfy the relationship , and the predicted values are equal to the given data .

Proof 4

Substituting the original data into Eq. (4), the parameters can be derived. The parameters of FAGMO(1,1,) can be obtained from Eqs. (25) to (27). Obviously, . Therefore, the predicted values of the FAGMO(1,1,) are equal to the given data .

Theorem 4 demonstrates that the FAGMO(1,1,) model is accurate for predicting arbitrary sequences that can be modelled by Eq. (20), while the FAGM(1,1,) model cannot describe the sequences accurately owing to there always being a non-zero difference between the real parameters and the parameters .

Theorem 5

The optimised parameters are approximately equivalent to the parameters when the value of is very small; that is,

Proof 5

We first consider the difference between parameter and , which is


It is known that and the first-order derivative is


When , the derivative of is positive, which indicates that the function is a monotonically increasing function in the interval [-2,2]. Thus, the value approaches zero as decreases. Therefore, when is very small.

Secondly, the difference between and is expressed as


Owing to , we know that when .

The first-order derivative of is


which is also positive when . Thus, decreases when the value of decreases and when is very small.

Thirdly, the difference between and is expressed as


It follows from and that when is very small.

From Theorem 5, we know that the differences between the parameters and are decrease along with smaller . Table 1 provides the values of and under different values of .

0.1 0.2 0.3 0.5 0.7 1.0 1.3 1.6 1.9 0.0001 0.0007 0.0023 0.0108 0.0309 0.0986 0.2506 0.5972 1.7636 0.0008 0.0034 0.0076 0.0217 0.0441 0.0986 0.1928 0.3733 0.9282

Table 1: Values of and under different values of

6 Modelling evaluation criteria and detailed modelling steps

To evaluate forecasting accuracy of the FAGMO(1,1,) model, the root mean squared percentage error (RMSPE) is applied to the prior-sample period (RMSPEPR) and post-sample period (RMSPEPO). In general, the RMSPEPR, RMSPEPO and RMSPE are defined as


where is the number of samples used to construct the model and is the total number of samples.

The index of agreement of the forecasting results is defined as


which is also a useful performance measure for sensitivity to differences in the observed and predicted data, where is the average sample value.

The average forecasting error (AE) and the mean absolute forecasting error (MAE) are


where AE reflects the positive and negative errors between the predicted and observed values, while MAE is applied for estimating the change in the forecasting model.

The detailed modelling steps of the fractional FAGMO(1,1,) are provided below.

Step 1: Determine the original data series , and -AGO series .

Step 2: Calculate the matrices and to determine using Eq. (4).

Step 3: Compute the parameters by employing Eqs. (25) to (27).

Step 4: Substitute the values of and into Eq. (20) to compute the predicted values .

Step 5: Apply the -IAGO matrix to obtain the restored values .

7 Validation of FAGMO(1,1,) model

This section provides numerical examples to validate the accuracy of the FAGMO(1,1,) model compared to the FAGM(1,1,) model and others.

7.1 Validation of FAGMO(1,1,) and FAGM(1,1,) models

This subsection presents a numerical example to validate the accuracy of the FAGMO(1,1,) and FAGM(1,1,) models. The values and are provided in the interval [0.01, 2] and [-1.99, 1.99], respectively. The initial point

is randomly generated in the interval [1, 2] by the uniform distribution, while the parameters

and are randomly generated in the intervals [0, 5] and [0, 100], respectively, by the uniform distribution. The other are generated with the aid of Eq. (20). All data used for the example are explained in Fig. 5.

Figure 5: Diagram of data for validation

We define the notation in the following analysis


where are the provided parameters of Eq. (20) and are the estimated parameters of the FAGM(1,1,) or FAGMO(1,1,) model.

When applying the above parameters, the graphs are displayed in Figs. 6 and 7. We observe from Fig. 6 that the maximum of FAGMO(1,1,) and FAGM(1,1,) are and 489.9434, respectively, where the magnitude is approximately 9034932. Furthermore, the of the FAGM(1,1,) model is very small when is near zero, which is coincident with Theorem 5. From Fig. 7, the maximum RMSPEs of FAGMO(1,1,) and FAGM(1,1,) are 0.0103% and 814.3864%, respectively, where the magnitude is approximately 79100.

It is known that the parameters , , and initial points are all randomly generated, which implies that the values of parameters , and have no influence on the output series. Here, the values and are the most important factors affecting the accuracy of the grey models.

Figure 6: Values of of FAGMO(1,1,) (left) and FAGM(1,1,) (right) models
Figure 7: Values of RMSPE of FAGMO(1,1,) (left) and FAGM(1,1,) (right) models

7.2 Validation of FAGMO(1,1,) model and other grey models

This subsection further demonstrates the advantage of the FAGMO(1,1,) model using two real cases.

Case 1: (Predicting cumulative oil field production). We consider an example from the paper Ma2016PredictingJCTN that provides sample data. The data from 1999 to 2009 are applied to construct the grey model, while the data from 2010 to 2012 are used for prediction. The values are listed in Table 2, indicating that the FAGMO(1,1,) model outperforms the other models in this case.

Year Data ENGM FAGM(1,1) FAGM(1,1,) FAGMO(1,1,)
1999 73.8217 73.8217 73.8217 73.8217 73.8217
2000 136.8817 138.4900 138.1621 137.1758 136.4573
2001 195.0590 195.4541 195.5377 196.1598 195.7633
2002 247.8547 247.9776 247.7638 249.3183 249.1781
2003 297.0902 296.4067 295.7629 297.2895 297.2750
2004 342.6394 341.0604 340.1238 341.0008 341.0322
2005 382.4312 382.2332 381.2700 381.2882 381.3320
2006 420.0399 420.1964 419.5291 418.8204 418.8699
2007 454.0430 455.2001 455.1670 454.1099 454.1712
2008 485.1171 487.4752 488.4068 487.5452 487.6290
2009 519.8508 517.2342 519.4402 519.4217 519.5393
2010 552.6569 544.6734 548.4350 549.9665 550.1281
2011 581.6092 569.9736 575.5400 579.3572 579.5714
2012 608.1863 593.3015 600.8887 607.7346 608.0086
RMSPEPR 0.4521% 0.4582% 0.3539% 0.3259%
RMSPEPO 2.0066% 1.0185% 0.3617% 0.3332%
Table 2: Results of ENGM, FAGM(1,1), FAGM(1,1,) and FAGMO(1,1,) models

Case 2: (Predicting foundation settlement close neighbouring Yangtze River). We consider an example from the paper Chen2014FoundationMPE, which provides sample data to construct the grey model. The values are presented in Table 3, indicating that the FAGMO(1,1,) model outperforms the other models in this case.

Day Data ONGM(1,1,,) FAGM(1,1) FAGM(1,1,) FAGMO(1,1,)
10 23.36 23.3600 23.3600 23.3600 23.3600
20 43.19 42.1779 43.3517 43.0586 43.0644
30 58.73 59.2549 59.4403 58.8205 58.8124
40 70.87 72.8374 72.4009 71.9763 71.9545
50 83.71 83.6405 82.8451 83.0247 82.9932
60 92.91 92.2330 91.2620 92.2158 92.1789
70 99.73 99.0672 98.0442 99.6885 99.6491
80 105.08 104.5030 103.5079 105.5215 105.4805
90 109.73 108.8264 107.9079 109.7568 109.7127
100 112.19 112.2652 111.4497 112.4117 112.3598
110 113.45 115.0002 114.2991 113.4857 113.4181
RMSPE 1.2730% 1.3257% 0.6030% 0.6011%
Table 3: Results of ONGM(1,1,,), FAGM(1,1), FAGM(1,1,) and FAGMO(1,1,) models

8 Applications

In this section, the FAGMO(1,1,) model is applied to forecast the nuclear energy consumption of China. The computational results of the FAGMO(1,1,) model are compared to the ENGM Ma2016PredictingJCTN, ONGM(1,1,) Chen2014FoundationMPE, FAGM(1,1) Wu2013GreyCNSNS and FAGM(1,1,) models.

8.1 Raw data

Raw data of the nuclear energy consumption of China were collected from the report of the BP Statistical Review of World Energy 2018. The first 10 samples belonging to the and the Five-Year Plans are applied to construct the prediction model, while the remaining samples of the Five-Year Plan are used to validate and compare the forecasting results (see Table 4).

Year Data Year Data Year Data
2006 12.4 2011 19.5 2016 48.2
2007 14.1 2012 22.0 2017 56.2
2008 15.5 2013 25.3
2009 15.9 2014 30.0
2010 16.7 2015 38.6
Table 4: Raw data of nuclear energy consumption of China, Mtoe

8.2 Simulation and prediction results

The simulation and prediction results are listed in Table 5 and Fig. 8, while the errors are listed in Table 6 and Fig. 9.

The nuclear energy consumption of China from 2016 to 2017 is predicted according to the established grey models. It can be observed in Table 5 and Fig. 8 that five grey models, namely ENGM, ONGM(1,1,), FAGM(1,1), FAGM(1,1,) and FAGMO(1,1,), successfully identify the trend of China’s nuclear energy consumption. However, these grey models differ from one another in terms of the prediction values from 2016 to 2020. From Fig. 8, China’s nuclear energy consumption is overestimated by the ENGM, ONGM(1,1,) and FAGM(1,1,) models, and underestimated by the FAGM(1,1) model. The values predicted by FAGMO(1,1,) are substantially closer to the raw data than those predicted by the other models.

We can observe from Table 6 and Fig. 9 that the RMSPEPR, RMSPEPO and RMSPE of FAGMO(1,1,) are 3.1409%, 4.1502% and 3.3304%, respectively. The RMSPEPR, RMSPEPO and RMSPE of ENGM are as high as 8.3788%, 30.3663% and 14.5667%, those of ONGM(1,1,) are 2.0494%, 12.0510% and 5.2635%, those of FAGM(1,1) are 4.8680%, 11.7968% and 6.5529%, and those of FAGM(1,1,) are 2.3299%, 6.3828% and 3.3636%, respectively. The IA, AE and MAE of FAGMO(1,1,) are 0.9985, 0.2526 and 0.7513, those of ENGM are 0.9538, 4.0536 and 4.0536, those of ONGM(1,1,) are 0.9911, 1.0225 and 1.1896, those of FAGM(1,1) are 0.9887, -1.1818 and 1.7105, and those of FAGM(1,1,) are 0.9971, 0.2736 and 0.8043, respectively. The computational results indicate that the FAGMO(1,1,) model outperforms ENGM, ONGM(1,1,), FAGM(1,1) and FAGM(1,1,), while ENGM exhibits the most inferior performance.

Year Data ENGM ONGM(1,1,) FAGM(1,1) FAGM(1,1,) FAGMO(1,1,)
2006 12.4 12.4000 12.4000 12.4000 12.4000 12.4000
2007 14.1 14.9788 14.4788 15.0242 14.7054 15.0891
2008 15.5 15.6744 15.1057 13.9808 15.0121 14.8608
2009 15.9 16.6846 16.0032 15.0566 15.8012 15.5886
2010 16.7 18.1520 17.2884 16.9219 17.0700 16.9760
2011 19.5 20.2831 19.1286 19.3953 18.9344 19.0534
2012 22.0 23.3785 21.7635 22.4687 21.5861 21.9432
2013 25.3 27.8741 25.5363 26.1951 25.3029 25.8633
2014 30.0 34.4036 30.9383 30.6625 30.4740 31.1013
2015 38.6 43.8871 38.6732 35.9872 37.6390 38.0473
2016 48.2 57.6610 49.7483 42.3129 47.5433 47.2178
2017 56.2 77.6662 65.6063 49.8133 61.2149 59.2933
2018 106.7219 88.3125 58.6959 80.0704 75.1679
2019 148.9226 120.8244 69.2074 106.0614 96.0147
2020 210.2150 167.3766 81.6403 141.8758 123.3723
Table 5: Simulation and prediction results of nuclear energy consumption by grey models

Figure 8: Comparison among five grey models for nuclear energy consumption
Year ENGM ONGM(1,1,) FAGM(1,1) FAGM(1,1,) FAGMO(1,1,)
2006 0 0 0 0 0
2007 0.0623 0.0269 0.0655 0.0429 0.0701
2008 0.0112 0.0254 0.0980 0.0315 0.0412
2009 0.0493 0.0065 0.0530 0.0062 0.0196
2010 0.0869 0.0352 0.0133 0.0222 0.0165
2011 0.0402 0.0190 0.0054 0.0290 0.0231
2012 0.0627 0.0108 0.0213 0.0188 0.0026
2013 0.1017 0.0093 0.0354 0.0001 0.0222
2014 0.1468 0.0313 0.0221 0.0158 0.0367
2015 0.1370 0.0019 0.0677 0.0249 0.0143
2016 0.1963 0.0321 0.1221 0.0136 0.0204
2017 0.3820 0.1674 0.1136 0.0892 0.0550
RMSPEPR 8.3788% 2.0494% 4.8680% 2.3299% 3.1409%
RMSPEPO 30.3663% 12.0510% 11.7968% 6.3828% 4.1502%
RMSPE 14.5667% 5.2635% 6.5529% 3.3636% 3.3304%
IA 0.9538 0.9911 0.9887 0.9971 0.9985
AE 4.0536 1.0225 -1.1818 0.2736 0.2526
MAE 4.0536 1.1896 1.7105 0.8043 0.7513
Table 6: Relative error values of nuclear energy consumption by five grey models
Figure 9: Errors among five grey models for nuclear energy consumption

8.3 Further discussions

As demonstrated by the case study, the novel FAGMO(1,1,) model outperforms other grey models. Moreover, it should be noted that in this paper we only conduct short-term forecasting, while it is well known that several existing energy models can perform long-term forecasting, such as LEAP, TIMES and NEMS. We will discuss the difference between our model and these models further, following a very brief introduction to such models.

LEAP (long-range energy alternatives planning system) schnaars1987houLRRP; Dong2017APS is a scenario-based energy environment modelling tool for climate change mitigation and energy policy analysis. It can be applied to examine energy production and consumption, as well as resource extraction in all sectors. The model studies the effects of various factors on energy consumption under different scenarios given an objective. LEAP is generally used for forecasting studies of between 20 and 50 years.

TIMES (The Integrated MARKAL-EFOM System) Shi2016modellingAE; Zhang2016timesAE

is an evolution of MARKAL, which was developed by the Energy Technology Systems Analysis Programme of the IEA. It combines technical engineering and economic approaches, and uses linear programming to produce a least-cost energy system under numerous user-specified constraints. The software is used to analyse energy, economic and environmental issues at different levels over several decades.

NEMS (National Energy Modeling System) Gabriel2001theOR; Soroush2017aIJSSD is a long-standing US government policy model, which computes equilibrium fuel prices and quantities for the US energy sector. NEMS is used to model the demand side explicitly; in particular, to determine consumer technology choices in the residential and commercial building sectors.

These models can perform long-term energy consumption projections. However, they may require a large amount of data, such as population growth, GDP, urbanisation, energy policies and energy strategies. In numerous practical situations, it is very difficult to obtain complete information because of time and cost limitations. The grey prediction model is an efficient method for conducting accurate forecasting with at least four samples. Compared to the major energy models, the grey model is an effective choice for predicting China’s nuclear energy consumption.

This paper collected 12 samples of China’s nuclear energy consumption from the BP Statistical Review of World Energy 2018. Thus, the LEAP, TIMES and NEMS models are all inapplicable owing to poor information. By employing the grey system theory and actual data from the and the Five-Year Plans, the FAGMO(1,1,) model was constructed. It can be observed in Table 5 that the prediction value of FAGMO(1,1,) is 123.3723 Mtoe in 2020, which is larger than the 84.6318 Mtoe provided in the BP energy outlook 2018. The main reasons for this are as follows.

i) It is infeasible to consider factors such as energy policies and China’s energy strategies, which affect the current situation of China’s nuclear energy consumption, in our proposed model because the FAGMO(1,1,) model is univariate. However, the forecasting models of institutions including BP, the IEA and APEC are based on widely collected data. Furthermore, grey models are mainly used for short-term forecasting in the calculation process, such as Feng2012ForecastingESPBEPP; Cui2013AAPM; duan2018forecasting; Wu2018usingE; Chen2008ForecastingCNSNS. Therefore, the forecasting results are relatively acceptable, reflecting the growth trend of future nuclear energy consumption in China.

ii) In China’s nuclear energy market, 38 nuclear power reactors are in operation, 19 nuclear power reactors are under construction and more are to be constructed by the end of 2016. This is the reason for the increase in nuclear energy consumption in recent years. However, no new nuclear projects have been approved for construction in 2016. Moreover, the State Council approved new safety rules and a nuclear power development plan following Japan’s Fukushima Daiichi crisis in 2011. These factors have also resulted in a slight slowdown in China’s nuclear energy consumption.

In the future, nuclear energy could provide an important alternative to fossil fuels such as coal and oil, and its proportion of the total primary energy consumption will increase yearly. Based on our forecasting results using the FAGMO(1,1,) model, the future nuclear energy consumption of China will increase rapidly if no certain restrictions are placed thereon. This implies that higher management and technical levels are necessary to meet the safety and quality requirements. Therefore, China’s government and policy makers should pay additional attention to the safety and quality issues of nuclear energy to achieve long-term, environmentally friendly and low-carbon energy goals and lay the foundation for the sustainable development of China’s energy and economy.

9 Conclusions

By applying the grey modelling technique and parameter optimisation method, the fractional FAGMO(1,1,) model was proposed to predict China’s nuclear energy consumption of the Five-Year Plan, based on the updated data from 2006 to 2015. The forecasting results provide the growth trend of the future nuclear energy consumption of China, and also offer a guideline for policymaking and project planning.

It can be observed that FAGMO is quite easy to use, with satisfactory accuracy in short-term nuclear consumption forecasting. For long-term prediction, its error will be larger because only 10 samples are used for modelling. This study is expected to be able to forecast the energy consumption of other countries that share similar patterns of economic development and energy consumption structures, among others. Furthermore, the optimised method applied to improve the FAGM(1,1,) model can be used to improve other first-order grey models, such as NGBM(1,1), GMC(1,) and RDGM(1,). These are possible extensions and suggested directions for our future research.


This research was supported by the National Natural Science Foundation of China (No. 71771033), Longshan academic talent research supporting program of SWUST (No. 17LZXY20), Doctoral Research Foundation of Southwest University of Science and Technology (No. 15zx7141, 16zx7140) and the Open Fund (PLN201710) of the State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Southwest Petroleum University).