Abstract
In recent years, the nonhomogeneous grey model has received much attention owing to its flexibility and applicability of forecasting small samples. To improve further the prediction accuracy of the nonhomogeneous grey model, this paper is to introduce a new whitening equation with variable coefficient into the original nonhomogeneous grey model, which is abbreviated as . First of all, the detailed computational steps of the time response function of the novel model and the restored values of the raw data sequence are deduced through grey modelling techniques. Secondly, two empirical examples from the previous literature are conducted to prove the validity of the novel model. Finally, the novel model is applied to forecast natural gas demand of China, and the results show that the novel model has a better prediction performance compared with other commonly used grey models, including , , , and .
1. Introduction
Grey prediction model is one of the key branches in grey system theory [1], which has been widely used in various fields, including industry [2, 3], engineering [4, 5], and especially energy [6–8] on account of its abilities of sufficiently dealing with small sample-sized problems. Besides, many extensive models of the benchmark grey model have been developed in practical applications. For example, Xie and Liu [9] investigated the discrete grey model which is abbreviated as ; in their work, they also revealed the connection of to the traditional grey model. Kumar and Jain [10] utilized grey Markov, grey model with the rolling mechanism and singular spectrum analysis to forecast energy consumption in India. Hamzacebi et al. [11] presented an improved grey model through combining iterative and direct grey forecasting approaches in order to predict electricity demand of Turkey. Ding et al. [12] combined a new initial condition and rolling mechanism to improve the prediction performance of the traditional grey model and then employed the model to forecast electricity consumption of China. Chen et al. [13] put forward a grey approach to achieve time-series interval forecasting. Inspired by existing approaches to forecasting time-series intervals [13], Wu and Zhang [14] further optimized the grey method for forecasting time-series intervals. Ma and Liu [15] raised a novel time-delayed polynomial grey model and applied it to forecast natural gas consumption of China. Intharathirat et al. [16] used multivariate grey models to forecast municipal solid waste quantity in developing countries.
However, it should be noticed that these above models are all linear, which are unapplicable for nonlinear problems. As a result, Chen et al. [17] firstly propounded a novel nonlinear grey Bernoulli model by introducing the Bernoulli equation into the grey action quality. Afterward, Wang et al. [18] proposed an improved nonlinear grey Bernoulli model which is based on optimization of background value. Under the inspiration of the literature [18], Wu et al. [19] investigated a new approach to improve the prediction performance by simultaneously considering the new initial condition and dynamic background value. Nguyen et al. [20] proposed a nonlinear grey Bernoulli model based on Fourier transformation. Ma et al. [21] proposed a novel nonlinear multivariable grey Bernoulli model to forecast tourism income of China, among them. It is worth noting the accumulation appeared in these models is all integer-order accumulation that impairs the prediction accuracy of forecasting models. For this, Wu et al. [22] originally proposed a grey model with fractional-order accumulation. Later, they also proposed a novel fractional multivariate grey model to forecast electricity consumption of Shandong Province [23]. Considering that computation of the existing fractional accumulation is complicated, therefore, Ma et al. [24] introduced a novel definition of conformable fractional accumulation, which is helpful in developing grey forecasting models.
Besides, the nonhomogeneous grey model, proposed by Cui et al. [25], also has received extensive attention owing to its ease of construction and flexibility. Chen and Yu [26] proposed a novel nonhomogeneous grey model in which the grey action quality was replaced by . Meanwhile, Ma et al. [27] investigated a novel kernel regularized nonhomogeneous grey model. Zeng et al. [28] developed a self-adaptive intelligence grey model. Based on the aforementioned knowledge, this paper introduces a new whitening equation with variable coefficient into the nonhomogeneous grey model to further enhance the prediction ability of the traditional nonhomogeneous grey model, which is abbreviated . The novelties of this paper are drawn as follows: (i) a new whitening equation with variable coefficient is considered into the nonhomogeneous grey model; (2) the time response function and the restored values of the original series are deduced in detail; (3) the two empirical examples are used to validate the effectiveness of the novel model; and (4) the novel model is applied in natural gas demand of China.
The rest of this paper is organized as follows: Section 2 briefly depicts the model. In Section 3, the modelling procedure of the novel model is discussed in detail. Section 4 provides two empirical examples to validate the novel model. Section 5 utilizes the novel model to forecast natural gas demand of China, and the main conclusions are listed in Section 6.
2. Description of the Traditional Nonhomogeneous Grey Model
To further improve the prediction performance of the model, Chen and Yu presented the nonhomogeneous discrete grey model, namely, . The computational steps of are described as follows.
Step 1. Assume a nonnegative series to beand then the first-order accumulative generating operation () series iswhere .
Step 2. The whitening equation of can be given bywhere and represent the development parameter and grey action quality, respectively.
Step 3. The discrete form of equation (3) is easily written aswhere is called the background value and .
Step 4. The model parameters , , and can be estimated through the least square method; there iswhere
Step 5. The time response function of equation (3) with is given by
Step 6. The restored value of , can be acquired by using the first-order inverse accumulative generating operation ; there is
3. Presentation of the Model
Definition 1. Let ; the whitening equation of a novel grey forecasting model isIt is worth noting that equation (9) will be degenerated to equation (3) when . In other words, the novel model is a more general form of .
Theorem 1. The discrete formula of the whitening equation of the novel model iswhere .
Proof. According to equation (9), it is easily yieldedIntegrating both sides of equation (11) over the interval ,Using the two-point trapezoidal formula, equation (12) becomesThis completes the proof.
Similar to Section 2, the model parameters of the novel model can be given bywhere
Theorem 2. Given the model parameters , , and of the novel model, the corresponding time response function can be given by
Proof. Multiply both sides of equation (9) by :That is,Integrating both sides of equation (18) over the interval ,Furthermore, by simplifying equation (19),where .
In particular, can be approximately calculated by the following equation:where is taken as 0.001 for any .
It is proved.
4. Validation of the Model
This section provides two empirical examples from the previous literature to demonstrate the effectiveness of the novel model; the competitive models including the traditional grey model , the discrete grey model , the nonlinear grey Bernoulli model , and the model are established thereof. To assess the prediction accuracies of these models, two statistical indices are used, which are the root mean square error (RMSE) and the mean absolute percentage error (MAPE):where is the simulative error and . The MAPE criteria of forecasting accuracies of these models are also given in Table 1.
Case 1. The data from paper [29] are taken as an example to verify the novel model, which describes wind turbine capacity of Europe from 2007 to 2017, the data from 2007 to 2015 are used to build prediction models, and the left two samples are used to examine the accuracies of these models. The prediction results by these models are given in Table 2, and curves of raw data and simulated and predicted values are graphed in Figure 1.

Case 2. The data from paper [30] are taken as the second example which depicts the carbon dioxide emissions of India from 2000 to 2018. Similar to Case 1, the first 17 data are used to build the prediction models, and the other data are used to examine the prediction accuracies of these models. Table 3 lists the prediction results, and Figure 2 exhibits the curves of raw data and predictive values.
From the above computational results, the desired conclusions can be drawn as follows: In Table 2, for simulated period, the RMSE values of competing models including , , , , and the proposed model are 1877.41, 1880.05, 1909.43, 4325.85, and 586.52; the MAPE values of these models are 1.58%, 1.54%, 1.79%, 4.22%, and 0.44%, respectively. For predicted period, the RMSE values of competing models are 7037.15, 7378.88, 5467.26, 6920.28, and 859.14; the MAPE values of these models are 3.99%, 4.19%, 3.07%, 4.00%, and 0.50%, respectively. It is seen that the proposed model outperforms among these models owing to its lowest RMSE value and MAPE value either in simulated period or in predicted period. In Figure 1, it is known that simulated and predicted values are more close to actual data. In Table 3, for simulated period, the RMSE values of these models are 30.83, 30.86, 31.66, 25.37, and 19.83; the MAPE values of these models are 1.59%, 1.60%, 2.03%, 1.50%, and 0.99%, respectively. For predicted period, the RMSE values of these models are 105.93, 107.18, 83.74, 50.7323, and 35.25; the MAPE values of them are 4.42%, 4.47%, 3.48%, 1.93%, and 1.09%, respectively. It is obviously found that the proposed model has a better prediction performance with the lowest RMSE values and MAPE values both in simulated period and predicted period. From Figure 2, it is easily seen that the predicted values of the carbon dioxide emissions of India are more close to actual data.In summary, the proposed model outperforms among these competing models listed in this paper comparing the prediction accuracies of these models. Due to the feasibility of the proposed model, it will be used to predict China’s natural gas demand in the next section.

5. Application
Natural gas plays an important role in the energy market on account of its high-quality, efficacy, and low-carbon feature. Therefore, development of natural gas has been written in the 13th Five Plan of China which is the most important official policy in China.
This paper collects the data of natural gas demand of China from 2001 to 2013 (cf. [28]), as shown in Table 4. The first ten data are used for building the novel model and other competing models, and the left 2 samples are used to examine the prediction accuracies of these models.
Furthermore, the simulated and predicted values of natural gas demand of China and statistical indices by the five grey models are listed in Table 5.
Ignoring the first item of predicted value, the simulated and predicted values are shown in Table 5, and these predicted values are graphed in Figure 3. It should be notable that the RMSE value of competitive models including , , , , and the novel model is 30.83, 30.86, 31.66, 25.37, and 19.83; the MAPE values of these models are 1.59%, 1.60%, 2.03%, 1.50%, and 0.99%, respectively. What is more, the forecasting ability in predicted period is fact more important; for the predicted stage, the RMSE values of these models are 15.07, 15.30, 13.91, 3.40, and 2.80; the MAPE values of them are 13.19%, 13.39%, 12.18%, 2.24%, and 1.97%, respectively. It is obviously seen that the proposed model has best prediction performance than that of other grey system models. From the perspective of plotting curves of raw data and simulated and predicted values, it is shown that the fitting values of the proposed model are more close to actual data. Therefore, the proposed model can improve the prediction performance significantly and should be regarded as an appreciate method to predict natural gas demand of China.

6. Conclusion
Aimed to improve further the prediction performance of the nonhomogeneous grey model, a new whitening equation with variable coefficient is introduced into the nonhomogeneous grey model; as a result of this paper, is proposed. The main contributions of this paper can be summarized as follows:(1)A new whitening equation with variable coefficient is introduced into the traditional nonhomogeneous grey model to further generalize the mathematical form of the nonhomogeneous grey model.(2)By using grey modelling techniques, the time response function of the novel model and restored values of the original series are deduced in detail.(3)Two empirical examples are used to prove the validity of the novel, and the novel model is applied to predict natural gas demand of China. The results indicate that the novel model has a better prediction performance compared with other competitive models.
Up to this point, the advantages of the proposed model have been discussed. However, there still exist some issues, which we should concentrate and solve. For instance, the model is regarded as a kind of univariate model, that is, the other influential factors could be ignored in applications, and the multivariable model with shape parameter will be discussed in our future work; the fractional-order accumulation should also be introduced in this study to further enhance the prediction performance.
Data Availability
Raw data used for this work can be found in the existing literature.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grant no. 11661001) and the Natural Science Foundation of Guangxi (Grant no. 2016GXNSFBA380069).