Forecasting solar energy consumption using a fractional discrete grey model with time power term

In order to more accurately predict the global solar energy consumption, a new grey prediction model FDGM(1,1,tα\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t^{\alpha }$$\end{document}) is proposed in this paper. The grey wolf optimizer (GWO) is used to optimize the fractional-order r\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r$$\end{document} and the time power α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} in the model. The proposed FDGM(1,1,tα\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t^{\alpha }$$\end{document}) model, two statistical models and other five existing grey models are used to simulate and predict the solar energy consumption in the four economies (France, South Korea, OECD and Asia Pacific region) from 2010 to 2029. The simulation results show that our proposed FDGM(1,1,tα\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t^{\alpha }$$\end{document}) has higher accuracy than the other seven models. The prediction results based on FDGM(1,1,tα\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t^{\alpha }$$\end{document}) show that in 2029, the solar energy consumption in South Korea, the OECD and the Asia Pacific region will reach 33,935.32 Ten-trillion J, 1,222,123.45 Ten-trillion J and 2,297,274.45 Ten-trillion J, respectively, but that in France will slowly increase in the next few years and will gradually decrease after reaching a peak of 14,727.34 Ten-trillion J in 2026. The better forecasting solar energy consumption by using the proposed model can provide useful information for the policy-makers to formulate and improve energy policies and measures.


Introduction
Since the energy crisis and environmental degradation have attracted worldwide attention, the search for new clean energy sources has become a popular topic in energy research. For example, electric vehicles are used to replace traditional fuel vehicles (Ehsani et al. 2021;Singh et al. 2021). Solar energy is the most abundant energy resource on Earth. The energy radiated by the sun to the Earth's atmosphere is only a fraction of the total solar radiation energy, but it accounts for 173,000 TW. In other words, the energy radiated by the sun to the Earth per second is equivalent to 5 million tons of standard coal, which is more than 10,000 times the total energy consumption of the entire world (Pierce 2016). As an important renewable energy source, solar energy is clean, efficient, endless and relatively stable and is therefore considered to be an effective way to solve the energy crisis (Singh et al. 2019). Solar energy has become an emerging industry of common concern and high priority worldwide. Photovoltaic power generation technology has been widely promoted and applied in countries with abundant solar energy resources, such as Japan, France, the UK and the USA. In addition, various countries have implemented incentive measures to promote the use of solar energy, such as electricity price subsidies and bank discount loans, which are conducive to the large-scale promotion of solar energy.
According to data from the International Renewable Energy Agency, the total scale of global photovoltaic installation is increasing year by year. By the end of 2019, the global grid-connected solar energy capacity reached 580.1 GW. According to the BP Statistical Review of World Energy 2020, global solar energy consumption (SEC) increased from 0.32 EJ in 2010 to 6.45 EJ in 2019, a nearly 20-fold increase in ten years. Table 1 shows the SEC of the economies of France, South Korea, the OECD and the Asia Pacific region from 2010 to 2019. From 2010 to 2019, the SEC of France, South Korea, the OECD and the Asia Pacific region increased significantly, with increases of 1697%, 1320%, 1036% and 5419%, respectively. In 2019, the SEC of the OECD and Asia Pacific region accounted for 53.5% and 54.4% of global SEC, respectively. Forecast of energy consumption is an essential basis for policymakers to formulate and improve energy policies and measures. Therefore, accurate prediction of SEC in developed economies, such as France, South Korea, the OECD and the Asia Pacific region, in the next 10 years will help us to understand the global trend of solar energy utilization in the future, provide data support for the formulation of sound solar energy industry development policies and have considerable significance for ensuring the sustainable development of global energy and the economy.
In existing research, various models are used to predict energy consumption, for example, the hybrid forecasting system (Du et al. 2018), computational intelligence technology (Meenal and Selvakumar 2018), Granger causality analysis (Pinzón, 2018), NEMS model (Soroush et al. 2017), LEAP model (Dong et al. 2017), LSTM model (Chen et al. 2019) and grey prediction model (Tsai 2016;Zeng et al. 2018;Wu et al. 2019;Ye et al. 2019;Guo et al. 2020;Wang et al. 2020;Xie et al. 2020;Liu et al. 2021;Şahin 2021). In all the prediction models, the grey prediction model has attracted substantial attention due to its advantages of convenient use, simple modelling process and high accuracy ). An important prediction model, the grey model was proposed by Professor Deng. The grey differential prediction model was established based on a little incomplete information to more accurately describe development laws (Deng 1982). Professor Deng initially proposed a firstorder univariate model GM(1,1) and its whitening equation where a is the development coefficient and b is the grey action quantity. However, GM (1,1) model can only realize the simulation of homogeneous exponential series . In order to overcome the limitations of the GM (1,1) model, scholars have successively proposed grey extended models suitable for various data forms. For example, Cui et al. (2013) built NGM (1,1,k) model with a grey action quantity of bt suitable for nonhomogeneous exponential series. Chen and Yu (2014) proposed NGM (1,1,k,c) model with a grey action quantity of bt + c suitable for approximate nonhomogeneous exponential series. Qian et al. (2012) established GM (1,1,t ) model with a grey action quantity of bt + c suitable for partial exponential series. Luo and Wei (2017) proposed GMP (1,1,N) model with a polynomial grey action quantity of 0 + 1 t + … + N t N suitable for partial exponential and polynomial series. The above studies showed that the model with improved grey action quantity had good prediction accuracy. However, GM (1,1) and the above models still do not overcome the inherent error of grey model from discrete estimation to continuous prediction. The discrete grey prediction model DGM(1,1) proposed by Xie and Liu(2009) avoided the mutation from difference equation to differential equation, and was unbiased for the fitting of pure exponential series. Similarly, Xie et al. (2013) established NDGM(1,1) model on the basis of DGM(1,1) with a grey action quantity of bt + c . Yang and Zhao (2016) proposed a discrete grey power model DGM(1,1,t ) based on DGM(1,1) with a grey action quantity of bt + c . Luo and Wei (2019) proposed DGMP (1,1,N) model based on DGM(1,1) with a polynomial grey action quantity of 0 + 1 t + … + N t N . Luo et al. (2020) introduced a trigonometric function into a discrete grey prediction model and proposed a discrete grey model with a time periodic term that is suitable for composite series data with periodicity and trends.
Although the method of establishing discrete model can improve the prediction accuracy of the model to a certain extent, there is still room for optimization in the model. The new information priority principle holds that new information plays a greater role in system cognition than old information. In the process of modelling, giving new information greater weight can improve the efficiency of grey modelling ). There are two ways to introduce new information: One way is to introduce the fractional accumulation operator into the grey model. Wu et al. (2013) proposed a new fractional order cumulative grey system model, in which the modelling sequence is no longer limited to first-order accumulation, but fractional accumulation, so as to better reflect the priority of new information in the accumulation process. Due to its excellent performance and accuracy, this approach has attracted the attention of many scholars. Mao et al. (2015) proposed a fractional-order cumulative delay model GM(1,N,). Wu et al. (2019) first proposed the fractional-order nonlinear grey Bernoulli model FANGBM(1,1). Liu et al. (2021) established a new fractional order grey polynomial model FPGM (1,1, t ). These models greatly enrich the fractional-order models of continuous equations. Fractional accumulation operator has also been widely used in discrete grey model. For example, Wu et al. (2014) proposed a fractional order cumulative discrete grey model.  proposed a discrete fractional-order accumulation model FAGM(1,1,D) to predict China's carbon emissions. Liu et al. (2016) proposed a fractional-order inverse model RDGM(1,1). Ma et al. (2019) proposed a novel fractional discrete multivariate grey model FDGM(1,1). The other way is to introduce new information in grey models including the rolling mechanism ) and metabolism (Zhao et al. 2012;Ma et al. 2013). The main difference between these approaches is that when the oldest data in the sequence are deleted, the new data added by the rolling grey model are real data, while the new data added by the metabolism model are predicted data (Yuan et al. 2017). Therefore, the rolling mechanism can only be used to predict one period, while a metabolic model can be used to predict multiple periods.
The above research has greatly promoted the development and application of grey models. However, the improved GM(1,1) model group is suitable only for the fitting and prediction of series data with the accumulated rule; it cannot reflect the priority of new information and is rarely used for the prediction of SEC. Therefore, combining the advantages of existing grey models and processing methods, a new grey model FDGM(1,1,t ) is proposed and used to forecast SEC in France, South Korea, the OECD and the Asia Pacific region. The empirical results show that the new model outperforms the Holt, Power, GM(1,1), FGM(1,1) (Wu et al. 2013), DGM(1,1) (Xie and Liu 2009), (Ma and Liu 2017), NGM(1,1,k,c) (Chen and Yu 2014), ARGM(1,1) (Wu et al. 2015) models. The main contributions of this paper are as follows.
(1) In the FDGM(1,1,t ) model, bt + c is used as the grey action quantity of the model to capture the exponential and power term characteristics of the data. The discretization technology is introduced to reduce the error caused by model jump. Fractional order accumulation and subtraction operators are introduced to make better use of new data and improve the prediction accuracy of the model.
(2) The optimal model structure parameters are obtained by GWO algorithm.
(3) The new model is used to fit six groups of numerical sequences with different forms, and compared with the other seven models showing the excellent fitting effect of the new model.
(4) The new model is used to predict the SEC of four countries or regions over the next 10 years.
The rest of this paper is organized as follows. The second section introduces the basic theory of fractional-order accumulation generation and the basic knowledge of the FDGM(1,1,t ) model, including the modelling mechanism, model characteristics and solution method, and describes the process of using the grey wolf optimizer to solve the model and the error measurement. In the third section, the feasibility and effectiveness of the model are verified by numerical examples and cases. The fourth section forecasts the SEC of France, South Korea, the OECD and the Asia Pacific region in the next ten years. Finally, the fifth section concludes the paper.

The basic theory of fractional-order accumulation
The definitions of the fractional-order accumulation operator and subtract operator are as follows (Wu et al. 2013;Meng et al. 2016): Definition 1. Given an original nonnegative sequence When r = 1 , the fractional-order accumulation operator is the first-order accumulation operator (1-FOA), and

Definition 2 Given an original nonnegative sequence
When r = 1 , IFOA is a first-order differential, also known as the first inverse fractional-order accumulation operator (1-IFOA), and , and the r-th order inverse accumulation operator is ; then, they have the following relations: is defined by Definition 1; then, the expression of the FDGM(1,1,t ) model is It can be seen from the expression of the model that when the parameters r, change, the model is transformed into other models.
Assume there are Y, B as mentioned at the end of the current Section, and that The time response function in the DGM(1,1) model Additionally, the reduction value is: When r = 0, = 0 , the FDGM(1,1,t ) model can be converted to the ARGM(1,1) model.
Assume there are Y, B as mentioned at the end of the current Section, and that The time response function of the ARGM(1,1) model t h e o r i g i n a l s e q u e n c e X (0) = {x (0) (1), x (0) (2), … , x (0) (n)}, r ∈ R + , the least square criterion of the FDGM(1,1,t ) model can be described as the following unconstrained optimization problem.

Solution of the model
Inserting Eq. (10) into Eq. (11) yields From the above calculation process, it can be concluded that According to Definition 3, To minimize the error of the model, the optimal values of the parameters r, must be determined.
The above optimization problems are essentially nonlinear programming problems with equality constraints, which can be solved by intelligent optimization algorithms or heuristic algorithms. In this paper, the GWO is used to solve the problem. The following figure shows the solution process of the model (Fig. 1).

Grey Wolf Optimizer (GWO)
The GWO is a population intelligent optimization algorithm proposed by Mirjalili in 2014. Inspired by the preying activity of the grey wolf, this algorithm is an optimized search method with many advantages.
Step 1: Fitness function. This paper evaluates the adaptability of each grey wolf's position (parameter) in the contemporary group and evaluates the grey wolf's position through adaptability. Moreover, MAPE is used as the fitness function to evaluate the position.
Step 2: Social hierarchy. In the design of the GWO, we build the grey wolf social hierarchy model, calculate the fitness of each individual in the population and mark the three grey wolves with the best fitness as and the remaining grey wolves as ω. The optimization process of the GWO is guided mainly by the best three solutions ( ).
Step 3: Surround the prey. Under the guidance of the three currently optimal grey wolves, the other candidate wolves update their positions randomly near their prey. As the grey wolves search for prey (the optimal parameter), they will gradually approach the prey and surround it. The mathematical model is expressed as follows: where t represents the current iteration, ⃗ V p(t) represents the current position of the prey, ⃗ V(t) represents the current position of the grey wolf, and ⃗ D represents the distance between the grey wolf and prey. Furthermore, ⃗ A and ⃗ C are synergy coefficient vectors defined as: In the whole iteration process, ⃗ a is reduced from 2 to 0, and ⃗ r 1 and ⃗ r 2 are random vectors in [0,1].
Step 4: Optimal location update. In each iteration, the best three grey wolves in the current population and their position information are retained to update the positions of the other candidate wolves. The mathematical model is expressed as follows: Step 5: Iteration stop criterion. Steps 3 and 4 are repeated until the specified maximum number of iterations is reached.
The maximum number of iterations is set to T = 200 , and when t > T , the iterative update is stopped.

Metabolic idea
Because the grey model has some error in predicting longterm data, the introduction of metabolism can make the long-term prediction more accurate. The mechanism of metabolism is as follows: Step 1: Build a model with the original sequence and predict the data of two periods.
Step 2: Delete the first two data points of the original sequence, add the predicted two periods of data, establish the model and predict the next two periods of data.
Step 3: Repeat step 2 to obtain all the predicted values (Fig. 2).

Error metric
This section considers several error measures widely used in the prediction model to test the effectiveness and applicability of the grey prediction model, as shown in Table 2.  ...
Step 1 Step 2 Step 3: repeat Step 2 Repeat Step 2 to obtain the forecast value for the next 10 years

Validation of FDGM(1,1,t˛)
This section uses a number of examples to verify the accuracy of the proposed model and the Holt, Power, GM(1,1), FGM(1,1), DGM(1,1), NGM(1,1,k,c) and ARGM(1,1) models. Numerical simulations are conducted in section "The SEC of the OECD", and the raw data and grey model in Table 1

Numerical examples
To study the fitting performance of the FDGM (1,1,t ), Holt, Power, GM(1,1), FGM(1,1), DGM(1,1), NGM(1,1,k,c), ARGM(1,1) models, six different forms of data were fitted: rising convex, rising concave, falling convex, falling concave, strictly nonhomogeneous and nearly nonhomogeneous. The results of these models are shown in Tables 3, 4, 5, 6 and   Fig. 3. It can been seen that the fitting accuracy and prediction accuracy of the FDGM(1,1,t ) model proposed in this paper are higher than those of the other grey prediction models. The FDGM(1,1,t ) model has strong adaptability because it can adapt to different numerical sequences by adjusting parameters and r ( Table 8).

The SEC of France
According to the data of the Statistical Review of World Energy 2020, France's energy consumption has long been at the forefront of the world. In 2019, France's total primary energy consumption was 9.68 EJ, second in Europe, after Germany. However, France's renewable energy consumption is 0.61 EJ, ranking only fifth in Europe. Therefore, in recent years, the French government has focused on the development and application of clean energy. By 2050, wind power will account for 63%, photovoltaic 17%, hydropower 13%, geothermal and other renewable energy 7% of the French power structure. In addition, from 2010 to 2019, the growth rate of SEC in France was as high as 1697%. Therefore, this section takes France as an example to test the accuracy of the proposed model and other grey prediction models. First, we take FDGM(1,1,t ) as an example to describe the modelling process: The GWO algorithm is used to obtain the optimal r = −

3
The parameter ∧ is calculated by the least squares method.
According to the parameter ∧ , the prediction model is as follows.
Then, the prediction sequence can be obtained via the fractional subtraction in Definition 3.
By adopting the same modelling method for GM(1,1), FGM(1,1), DGM(1,1), NGM(1,1,k,c) and ARGM(1,1), the following prediction models are obtained.    (29), the prediction results of the models are shown in Fig. 4 and Table 9, and the error metrics of each model are shown in Fig. 5 and Table 10. The estimated value of the GM(1,1), DGM(1,1) model is significantly higher than the actual value, the estimated value of the Holt, Power, FGM(1,1), NGM(1,1,k,c), ARGM(1,1) model is significantly lower than the actual value, and the predicted value of the model FDGM (1,1,t ) is the closest to the actual value. On the basis of the prediction value and fitting value, the four error metrics of the FDGM(1,1,t ) model are the best among the eight prediction models. Moreover, the FDGM(1,1,t ) model is more precise than the other prediction models in predicting SEC in France.

The SEC of South Korea
According to the Statistical Review of World Energy 2020, in 2019, South Korea was the world's tenth largest energy consumer and eighth largest carbon dioxide emitter. Therefore, the South Korean government is facing the dual challenges of improving energy security and reducing greenhouse gas emissions. One of the most effective solutions is to  strengthen the role of renewable energy in power production. To this end, the South Korean government plans to increase investment in green energy and set a target of 20% renewable energy by 2030 to ensure that renewable energy, such as solar energy and wind energy, plays a decisive role in meeting the energy demand. Therefore, this section takes South Korea as an example to test the accuracy of the proposed model and other grey prediction models. First, we take FDGM(1,1,t ) as an example to describe the modelling process: The GWO algorithm is used to obtain the optimal r = 0.2641, = 2.7608 and generate X (r) = {761. 94, 1137.19, } {1471.16, 2139.09, 3301.13, 5084.92, 6786.94, 9385.02} via fractional accumulation. Then, according to X (r) , matrices Y and B are obtained: The parameter ∧ is calculated by the least squares method.
According to the parameter ∧ , the prediction model can be calculated as follows.
Then, the prediction sequence can be obtained via the fractional subtraction in Definition 3.

The SEC of the OECD
According to the data of the Statistical Review of World Energy 2020, the primary energy consumption of the OECD in 2019 was 233.43 EJ, accounting for 40% of the world's total, and the SEC was 3.45 EJ, accounting for 53.4% of the world's total. According to the International Energy Agency, renewable energy will provide one-third of the OECD's total power generation by 2035. OECD countries have rich experience in the development and utilization of renewable energy. Therefore, this section takes the OECD as an example to test the accuracy of the proposed model and other grey prediction models.
Then, according to X (r) , the matrices Y and B are obtained: The parameter ∧ is calculated by the least squares method.
Additionally, the following is obtained based on parameter ∧ Then, the prediction sequence is obtained via fractional subtraction (Definition 3).

The SEC of the Asia Pacific region
According to the data of the Statistical Review of World Energy 2020, the primary energy consumption in the Asia Pacific region reached 257.56 EJ in 2019, ranking first in the world, and accounting for approximately 44% of the world's total. Furthermore, SEC reached 3.51 EJ, accounting for (1) (t) = 22754.21t + 19308.6 (45) x (0) (t) = 1.0696x (0) (t − 1) + 23397.98 approximately 54.4% of the world's total. In the Asia Pacific region, China, India, Japan and South Korea are among the top solar energy consumers in the world. Therefore, this section takes the Asia Pacific region as an example to test the accuracy of the proposed model and other grey prediction models. First, we take FDGM(1,1,t ) as an example to describe the modelling process: The GWO algorithm is used to obtain the optimal r = 0.6089, = 3.2181 and generate X (r) = {6362. 61, } { 15298.99, 26415.67, 48178.99, 91267.32, 150504.56, } { 1234754.23, 366314.41} via fractional accumulation.
Then, according to X (r) , the matrices Y and B are obtained:  Then, the prediction sequence can be obtained via the fractional subtraction in Definition 3.
Holt model: trend line prediction    (53), the prediction results of the models are shown in Fig. 10 and Table 15, and the error metrics of each model are shown in Fig. 11 and Table 16.  The estimated value of the Power, GM(1,1), FGM(1,1), DGM(1,1) and ARGM(1,1) models is significantly higher than the actual value, the estimated value of the Holt, NGM(1,1,k,c) models is significantly lower than the actual value, and the predicted value of the FDGM(1,1,t ) model is the closest to the actual value. Moreover, the proposed FDGM(1,1,t ) model is more accurate than the other prediction models in predicting the SEC of the Asia Pacific region.

Forecasting SEC in the next ten years
In this section, we build an FDGM(1,1,t ) model that includes metabolism to forecast the SEC of France, South Korea, the OECD and the Asia Pacific region in the next 10 years (2020-2029). The prediction results are shown in Table 17 and Figs. 12, 13, 14 and 15. The SEC in South Korea, the OECD and the Asia Pacific region will gradually increase in the next decade. However, the SEC in France will rise slowly before 2026 and then begin to decline after reaching its peak in 2026. The consumption of the two major economies of the OECD and the Asia Pacific region was basically the same in 2019; however, after 2029, the distance between the two will rapidly increase, and the SEC in the Asia Pacific region will be almost twice that of the OECD. The main reason for this difference may be that the OECD is composed mainly of developed countries, and the utilization and industrial construction of solar energy are relatively mature compared to that in developing countries. The growth rate of SEC will drop from 17.02% in 2020 to 10.17% in 2029. In addition to developed countries, the Asia Pacific region also contains some developing countries, such as China and India, which are in the stage of rapid development of the solar energy industry. Thus, the average growth rate of SEC from 2020 to 2029 in the Asia Pacific region is higher than that of the OECD.

Conclusions
In order to more accurately predict the global SEC, a new discrete grey prediction model FDGM(1,1,t ) is proposed in this paper. Based on the traditional discrete grey model, the new model combines fractional order accumulation operator, subtraction operator and time power term. Specifically, using fractional accumulation operator and subtraction operator  to process the original data sequence can reduce the volatility and randomness of data, enhance the development trend of data and obtain more useful information from limited data. In addition, introducing the time power term into the grey action can better fit the nonlinear trend in the data and improve the prediction accuracy of the model. In order to further improve the accuracy of the model in fitting data, an intelligent optimization algorithm (GWO) is introduced to find the optimal parameters of the model. With the change of parameters, FDGM(1,1,t ) model can be transformed into DGM(1,1) and ARGM(1,1) models, so it has strong adaptability. Subsequently, the proposed FDGM(1,1,t ) model, two statistical models and other five existing grey models are used to simulate and predict the SEC in the four economies (France, South Korea, OECD and Asia Pacific region) from 2010 to 2029. The simulation results show that our proposed  FDGM(1,1,t ) has higher accuracy than the other seven models. Although the grey prediction model has high prediction accuracy, it cannot meet the demand for future data because it is not suitable for long-term prediction. Therefore, this paper establishes a model based on a metabolism mechanism to predict SEC over the next 10 years. The prediction results based on FDGM(1,1,t ) show that the SEC in South Korea, the OECD and the Asia Pacific region will gradually increase from 2020 to 2029, but that in France will slowly increase in the next few years and will gradually decrease after reaching a peak in 2026. Notably, the FDGM(1,1,t ) model proposed in this paper can be used not only for the prediction of SEC but also for the prediction of other energy sources. So far, the FDGM(1,1,t ) model has significant advantages compared with other prediction methods. However, the model can be further optimized in future research. For example, for data with periodic fluctuations, it is necessary to introduce a periodic function into the grey action quantity to capture the periodicity and volatility of the data. In addition, new fractional order accumulation and subtraction methods can be found to make the modelling data smoother and reduce the volatility and randomness of the data more effectively.  Fig. 12 The SEC and growth rate of France in the next ten years Fig. 13 The SEC and growth rate of South Korea in the next ten years Fig. 14 The SEC and growth rate of OECD in the next ten years Fig. 15 The SEC and growth rate of Asia Pacific in the next ten years