Abstract

In this work, a novel time-delayed polynomial grey prediction model with the fractional order accumulation is put forward, which is abbreviated as TDPFOGM(1,1), based on the new grey system theory to predict the small sample in comparison with the existing forecasting models. The new model takes into account the nonhomogeneous term and the priority of new information can be better reflected in the in-sample model. The data in this paper all come from the existing literatures. The results demonstrate that the TDPFOGM(1,1) model outperforms the TDPGM(1,1) and FOGM(1,1) model.

1. Introduction

Shale gas is one kind of fossil fuel which is clean, efficient, and cheap. With the shale gas revolution in the United States [1, 2], shale gas has been developed in many countries. And China has achieved some success [3]. The development of shale gas industry has been regarded as an important mission of the government [4, 5]. In such circumstances, it is very important to predict the natural gas production and consumption in China for the decision-makers, such as the government and energy companies.

Grey model proposed by Deng [6] is used for small sample forecasting. The grey prediction models play an important role in the grey system theory. These models are often named as grey models. GM(1,1), that is, first-order Grey model with one variable, has been successfully applied in many disciplines [713]. However, the existing GM(1,1) model cannot be accurate prediction for many actual systems. In recent literatures [1423], the improved grey prediction models are popularly used for energy consumption. Ma and Liu [14] have proposed a novel time-delayed polynomial grey model to predict the natural gas consumption in China. Akay and Atak [15] have proposed a novel method based on the basic GM(1,1) model with rolling mechanism to forecast the electricity demand of Turkey. Wu and Shen [16] have proposed a grey-related least squares support vector machine optimization model and applied it, predicting the natural gas consumption demand. The authors in [17] have proposed a novel kernel regularized nonhomogeneous grey model and applied it, forecasting the petroleum production. Kumar and Jain [18] have used the Grey-Markov and GM(1,1) model with rolling mechanism to predict the energy consumption in India. In [19], a prediction method using grey model for cumulative plastic deformation under cyclic loads was demonstrated. Pao and Tsai [22] have used the GM(1,1) model to predict the energy consumption in Brazil compared with the ARIMA model.

Discrete grey model also has good effect on prediction. Ayvaz and Kusakci [24] have used the nonhomogeneous discrete grey model to predict the electricity consumption forecasting for Turkey. Xie and Pearman [25] have used the discrete GM(1,1) model (DGM(1,1)) to predict the energy consumption of China. Simultaneously, the fractional gray model is proposed. In [26], the authors have proposed a grey system model with the fractional order accumulation. Yuan et al. [27] have pointed out the advantages and disadvantages of GM(1,1) and the Autoregressive Integrated Moving Average (ARIMA) model and then proposed a novel hybrid model based on the GM(1,1) model and the ARIMA model to predict primary energy consumption in China. Wu et al. [28] have proposed a novel GM(1,1) model with the Principle of New Information Priority (NIGM(1,1)) to predict the natural gas consumption in China, and the results show that NIGM(1,1) outperforms several existing models. Mao et al. [29] present a new fractional grey model, in which first-order differential equations are transformed into fractional differential equations. And it has high modeling precision and can overcome the GM(1,1) model class ratio test restrictions. Yang et al. [30] modified optimized fractional grey model using the error feedback, and the performance is evaluated and greatly improved in modeling. Meng et al. [31] develop a discrete grey model with fractional operators, which also makes use of genetic algorithms to optimize the modeling parameter. Yang et al. [32] proposed that the geometric coordinate features are used by the coordinates of area and middle point lines and established the grey prediction model for interval grey number by the fractional order accumulation calculus. Wu et al. [33] developed a novel multivariable grey forecasting model that considered the total population to forecast the electricity consumption. Li et al. [34] proposed a novel grey forecasting model with full-order time power terms (FOTP-GM(1,1)). Wang et al. [35] propose a data grouping approach based grey modeling method DGGM (1,1) to predict quarterly hydropower production in China. All these researches indicate that the grey models are efficient to predict the consumption of many kinds of energy for many countries.

Motivated by the above, in this paper, we propose a novel time-delayed polynomial grey model with the fractional order accumulation. Simultaneously, the model considered the character time-delayed polynomial and fractional order. The main contributions of this paper are summarized as follows: (1) this paper is based on predecessor’s research about grey model; a new model is put forward; (2) using the same data, the TDPFOGM(1,1) model makes predictions just as shown in figures. And the results were compared with the TDPGM(1,1) model and the FOGM(1,1) model; (3) more intuitively, the mean absolute percentage error (MAPE) is shown in the tables.

The rest of this paper is organized as follows: the details of modeling procedures of the TDPFOGM(1,1) model are given in Section 2; the application of TDPFOGM(1,1) to small sample prediction is given in Section 3, including the comparison with the other commonly used prediction models presented based on the same data sample, and the conclusions are drawn in Section 4.

2. Description of the Problem

In this section, we will present the time-delayed polynomial grey system model with the fractional order accumulation, abbreviated as the TDPFOGM(1,1) or model, including the principles and the computational steps.

2.1. Grey Model with the Fractional Order Accumulation

Let the th order accumulated generating operator of the original nonnegative sequence be ,

Set ; , ; thenwhere

Fractional derivatives accumulate the whole history of the system in weighted form. in grey system theory denotes the weight of as 1. The larger of is the larger the weight of old data is; the smaller of is the smaller the weight of old data is. Reducing can reduce the weights of old data, which can put more emphasis on the newer data (more details and the properties of the fractional order accumulation generation operation can be seen in [25]).

The original form of the model is as follows:

If , then order inverse accumulated generating operator of ; we write

2.2. The Representation of the TDPFOGM(1,1) Model

The new form of the model is as follows:where is called the background value. For model, the first number “1” stands for “first order” and the second number “1” means “one sequence,” as only one sequence is considered in this model. It is the traditional model [14], when . The ordinary least squares estimate sequence of the model is satisfied.

2.3. The Solution of the TDPFOGM(1,1) Model

is the given original sequence, and the unconstrained optimization equation (6) can describe the least squares criteria for the TDPFOGM(1,1) model:

The linear system in (7) can solve the above optimization problem:where

The equation

is called a whitenization differential equation of TDPFOGM(1,1), and the solution of the whitenization equation is given by the following equation:where

Equation (10) is called the discrete response function, and within the initial condition , the discrete response function can be obtained as follows:

The values of the series can be computed using the response function (12), and the predicted values of the original series can be obtained using the fractional order inverse accumulative generation operation as follows:The derivation of the parameter estimation of TDPFOGM(1,1) is similar to the other first order grey prediction models, such as TDPGM(1,1) (see [14]).

Remark 1. The TDPFOGM(1,1) model has the properties of fractional order accumulation and the polynomial function, and it is not applicable for the periodic series and the random series and so forth according to its mathematical formulation. Exponential function and a discrete integral with an exponential and the polynomial function combination explain the discrete response function of the TDPFOGM(1,1) model. The priority of new information can be better reflected when the accumulation order number becomes smaller in the model.

2.4. The Computational Steps

According to the principles of time-delayed polynomial grey system model with the fractional order accumulation, the computational steps can be summarized as follows.

Step 1. Within a given original sequence, compute the series of the given time series using (1) and compute the background values using (5).

Step 2. Substitute the original sequence along with its series and the background values into (7) and compute the parameters by solving the linear system in (7).

Step 3. Substitute the parameters into the discrete response function in (12) and then compute the series .

Step 4. Compute the predicted values of using the (13).

3. Model Application and Comparison

3.1. Raw Data Collection

The raw data of the natural gas consumption () of China are collected from [14] (page 20) as shown in Table 1. We build the prediction models using the data from 1995 to 2004 and validate the modeling accuracy using the data from 2005 to 2013.

We consider an example from paper [36]. The same sample is applied here to compare the precision. Actual values are presented in Table 2.

3.2. Evaluation Indices of the Modeling Accuracy

The mean absolute percentage error (MAPE) is used to evaluate the overall forecast performance of the prediction models, which is defined as follows:

3.3. Results and Analysis

Case 1. We consider the data from Table 1. The same sample is applied here to compare the precision. Actual values and the fitting values of four compared models are presented in Table 3 and the predicted values are also plotted in Figure 1. From Table 3, yielded lower MAPE compared with the traditional TDPGM(1,1), GM(1,1), and .

Case 2. We consider the data from Table 2. The same sample is applied here to compare the precision. Actual values and the fitting values of four compared models are presented in Table 4 and the predicted values are also plotted in Figure 2. From Table 4, yielded lower MAPE compared with the traditional GM(1,1) and .

Remark 2. In this paper, we established the novel grey prediction model that considers the fractional order accumulation while considering the polynomial based on the TDPGM(1,1) model which is different from the existing literature [14]. There are three main reasons; firstly, in [14], the authors proposed a time-delay polynomial grey model, but fractional order accumulation is not mentioned. Secondly, The TDPFOGM(1,1) model is a further study on the basis of the TDPGM(1,1) model and TDPGM(1,1) model is a special case when fractional order . Finally, The prediction performance of TDPFOGM(1,1) model is better than that of the TDPGM(1,1) model according to the numerical results and analysis by comparing the TDPFOGM(1,1) and TDPGM(1,1) models based on the same data sample.

4. Conclusion

Small sample forecasting is a difficult and important problem. We have improved the grey prediction model after studying the grey system model with the fractional order accumulation and the time-delayed polynomial grey model. Based on a series of analysis and derivation for grey prediction system, the following conclusions could be drawn:

(1) Time-delayed polynomial grey system model with the fractional order accumulation (TDPFOGM(1,1)) is established and it can be used for forecasting.

(2) The anticipated effect using the TDPFOGM(1,1) model is better than the TDPGM model and the FOGM model under the same data condition.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors acknowledge the support of the applied fundamental research projects (major frontier projects) of Sichuan Province (no. 2016JC0314).