Abstract

Sulfur dioxide is an important source of atmospheric pollution. Many countries are developing policies to reduce sulfur dioxide emissions. In this paper, a novel prediction model is proposed, which could be used to forecast sulfur dioxide emissions. To improve the modeling procedure, fractional order accumulating generation operator and fractional order reducing generation operator are introduced. Based on fractional order operators, a discrete grey model with fractional operators is developed, which also makes use of genetic algorithms to optimize the modeling parameter . The improved performance of the model is demonstrated via comparison studies with other grey models. The model is then used to predict China’s sulfur dioxide emissions. The forecast result shows that the amount of sulfur dioxide emissions is steadily decreasing and the policies of sulfur dioxide reduction in China are effective. According to the current trend, by 2020, the value of China’s sulfur dioxide emissions will be only 86.843% of emissions in 2015. Fractional order generation operators can be used to develop other fractional order system models.

1. Introduction

Sulfur dioxide (chemical formula SO2) is the most common form of sulfur oxides. The colorless gas has a strong irritating odor and is a major air pollutant. Inhaling it will lead to respiratory inflammation, bronchitis, emphysema, conjunctivitis, and other health problems. Sulfur dioxide will also weaken the immunity of young people, making them less resistant to infections. When sulfur dioxide reacts with water, sulfurous acid is formed [1]. If sulfurous acid is further oxidized by particulate matter in the atmosphere, sulfuric acid (the main component of acid rain) is rapidly and efficiently produced. With the presence of oxidants, light can catalyze the production of sulfate aerosols, which have an adverse effect on human health and could increase the mortality of patients [2]. Sulfur dioxide is the main cause of global “acid rain.” It is estimated that the direct economic losses caused by metal corrosion in industrialized countries account for between 2 and 4 per cent of the gross domestic product (GDP). Due to the serious negative impact of sulfur dioxide on flora, fauna, and buildings, sulfur dioxide emissions and the emission concentrations are under strict regulation. The concentration of sulfur dioxide is an important indicator of whether the air is polluted. Using coal as the main energy source, China’s sulfur dioxide emissions are ranked the first in the world [3]. Hence, it is of great practical significance to forecast the emissions of sulfur dioxide in order to formulate emission reduction policies.

The world has exerted great efforts to reduce sulfur dioxide emissions. There has also been a lot of research focusing on sulfur dioxide emissions. Lu and colleagues have studied sulfur dioxide emissions in China and sulfur trends in East Asia since 2000 [4] and also analyzed sulfur dioxide and primary carbonaceous aerosol emissions in China and India from 1996 to 2010 [5]. Smith et al. investigated global and regional anthropogenic sulfur dioxide emissions [6]. Ma et al. developed scenario analysis to study the reduction potential of sulfur dioxide emissions in China’s iron and steel industry [7]. Yang and Hu conducted research on the inventory system of sulfur dioxide emissions in China [8]. Smith et al. discussed the methods and results of historical sulfur dioxide emissions from 1850 to 2000 [9]. Li et al. observed recent large reduction in sulfur dioxide emissions from China’s power plants by using ozone monitoring instruments [10]. Hao et al. designated sulfur dioxide and acid rain pollution control zones and studied the impact on energy industry in China [11]. Mohajan studied China’s sulfur dioxide emissions and local environmental pollution [12]. Zhou and Zhang established an improved metabolism grey model for predicting small samples with a singular datum and applied it to study sulfur dioxide emissions in China [13]. An adaptive neurofuzzy logic method has been developed by Yildirim et al. to estimate the impact of meteorological factors on sulfur dioxide pollution levels [14]. Hassanzadeh and coworkers developed statistical models and time series to forecast sulfur dioxide [15]. Li et al. looked into the contribution of China’s emissions to global climate forcing [16]. Liu et al. believed that air pollutant emissions from Chinese households are a major and underappreciated source of ambient pollution [17].

To improve China’s policies on sulfur dioxide emission reduction, it is necessary to predict China’s future sulfur dioxide emissions. Grey prediction model is an important part of grey system theory [18, 19]. It has been studied by many researchers and improved over the years [2022]. Recently, Liu et al. [23], Wu et al. [24, 25], and Xiao et al. [26, 27] published results on fractional order accumulation and grey system model with the fractional order accumulation. However, all the existing research [28, 29] does not deduce the analytical expression of fractional order accumulating generation operator and fractional order reducing generation operator, nor does it prove that these fractional order generation operators satisfy mutual invertibility. Accumulating generation operator is similar to summation operator and reducing generation operator is similar to difference operator. Cheng and Wu [28, 29] proposed the analytic expression of fractional summation operator but did not propose fractional difference operator and believed that the fractional difference operator generally does not satisfy the exponential rule. In this paper, both the fractional order accumulating generation operator and fractional order reducing generation operator are studied. Furthermore, a discrete grey model with fractional operators is used to predict China’s future sulfur dioxide emissions.

The rest of the paper is organized as follows: in Section 2, the analytical expressions of fractional order accumulating generation operator and fractional order reducing generation operator are derived. In Section 3, the discrete grey model with fractional operators and order optimization algorithm is studied. This is followed by comparisons of the proposed model with the other discrete grey models, demonstrated via three case studies, in Section 4. In Section 5, the proposed discrete grey model with fractional operators is used to predict China’s sulfur dioxide emissions. Then, conclusions are drawn in Section 6.

2. Fractional Order Generation Operators

Definition 1. Assume that is a sequence of raw data and a sequence operator satisfying thatwhereand then the sequence operator is called a (first-order) accumulating generator of , denoted as 1-AGO.

Definition 2. For , we define as the th-order accumulating generation sequence of , where

Definition 3. Assume that the sequence of raw data is and is the first-order accumulating generation sequence of , whereand it follows that From Definition 2, it follows thatSo

Theorem 4. is the second-order accumulating generation sequence of , where

Proof. By induction on , consider the following.
First, for , , so it is true.
Next, assume that, for some and , the statement is true. That is, Then, for , it follows thatSo, Theorem 4 is true.

Theorem 5. Assume that the sequence of raw data iswhere , and is the th integer order accumulating generation sequence of , where

Proof. By induction on , consider the following.
For , . It is true.
For each , if is true, it follows thatSo Theorem 5 is true for all .

Definition 6. For and , is Gamma function for .
Using integration by parts, the Gamma function satisfies the functionIn particular, if is a positive integer,For example, , , ,
The Gamma function is an extension of the factorial function.
By Gamma function, it follows that

Definition 7. Assume that the sequence of raw data iswhere , and is the th integer order accumulating generation sequence of , where

Definition 8. Assume that is a sequence of raw data and a sequence operator satisfying thatwhereand the sequence operator is called the first-order forward reducing generation operator.
Ifthe sequence operator is called the first-order backward reducing generation operator.
By default, the reducing generation operator is a backward reducing generation operator. The first-order reducing generation operator is the first-order inverse accumulating generation operator (IAGO), written as 1-RGO.

Definition 9. For , we write as the th-order reducing generation sequence of , where

Definition 10. Assume that the sequence of raw data isand is the first-order reducing generation sequence of , whereand then it follows that From Definition 9, it obviously follows thatThen, it follows that

Theorem 11. is the second-order reducing generation sequence of , where

Proof. By induction on , consider the following.
For , For , For , So, Theorem 11 is true.

Theorem 12. Assume thatwhere , and is the th integer order reducing generation sequence of , where

Proof. By induction on , consider the following.
For , if ,if , For , if , if , if , For each , if is true, then it follows thatThen, it follows thatThe coefficient of for the expansion of is So,
So, Theorem 12 is true.

Proposition 13. For , the th integer order reducing generation operatorcan be simplified as

Proof. According to the property of Gamma function,
When , When ,for , then ;for , then .Thus, it follows thatSo, Proposition 13 is true.

Definition 14. Assume thatwhere , and is the th integer order reducing generation sequence of , where

Theorem 15. Assume that is a sequence of raw data, where and . is the sequence of the th-order accumulating generation of . is the sequence of the th-order reducing generation of . is the sequence of the th-order reducing generation of . is the sequence of the th-order accumulating generation of . The following holds true [19]:(1)If , is the order accumulating generation of .(2)If , is the order accumulating generation of .(3)The fractional accumulating generation operator and the fractional reducing generation operator satisfy the exchange law and the exponential rate.

Theorem 16. Assume that is a sequence of raw data and . is the sequence of the th-order accumulating generation of . is the sequence of the th-order reducing generation of . The th order accumulating generation operator and the th-order reducing generation operator are inverse operations [19].

3. Discrete Grey Model with Fractional Operators

Definition 17. Assume that and are defined as in Definition 7, andis the discrete grey model with fractional operators. In particular, consider the following.If , is .If , is direct modeling .

Theorem 18. Assume that is a sequence of raw data and is the sequence of the th-order accumulating generation of , where .is the th-order accumulating generation of , whereThe parameter vector of discrete grey model with fractional operators can be calculated by the least squares method.It follows thatwhereSo, it follows that

Theorem 19. Assume that , , and are the same as in Theorem 18. Ifthen the following is true.
(1) The time response sequence of the discrete grey model with fractional operators is given byor(2) The restored values of can be given by

Theorem 20. Discrete grey model DGM is the special case of the discrete grey model with fractional operators, where .

The conclusion is clearly true.

To determine the optimal order of the discrete grey model with fractional operators, we can solve the following optimization problem:

Genetic algorithms have obvious advantages at solving nonlinear optimization problems [30]. So, we can apply genetic algorithms to solve the optimized value of order .

Step 1 (start). Initialize the parameters of population : initial value of order , range of value for , population size , crossover probability , mutation probability , value accuracy represented as , fitness function results accuracy , and generations of terminating evolution .

Step 2 (fitness). Determine the fitness of population .
Step 2.1. Compute the order accumulating generation sequence .
Step 2.2. Compute the parameter vector .
Step 2.3. Compute the time response sequence .
Step 2.4. Compute the forecast value .
Step 2.5. Compute the fitness .

Step 3. New population.
Step 3.1 (selection). Select parents from population .
Step 3.2 (crossover). Perform crossover on parents creating population .
Step 3.3 (mutation). Perform mutation of population .
Step 3.4 (accepting). Place new offspring in a new population.

Step 4 (replace). Use new generated population for a further run of algorithm.

Step 5 (test). If the end condition is satisfied, stop, and return the best solution in current population; otherwise go to Step 2.

Step 6 (result). Display the result data of the model, the optimized order , the parameter vector , and the forecast value .

The overall modeling steps of the discrete grey model with optimized fractional operators are shown in Figure 1.

4. Verification of Discrete Grey Model with Fractional Operators

In this section, the advantage of the discrete grey model with fractional operators over the other discrete grey models is demonstrated by three case studies. Mean absolute percentage error (MAPE) is used to compare the actual values with forecasted values to evaluate the precision. MAPE is defined as

Case 1. In [31], a modified was proposed to predict China’s foreign exchange reserves. We consider the example from this paper [31] which provides the sample data. The actual values and the forecast values of the three compared models are presented in Table 1. As can be seen from Table 1, the discrete grey model with fractional operators where and achieved the lowest MAPE compared to and the modified .

Case 2. In [32], the initial condition of was optimized and the results of direct modeling were tested against the original data, in respect of a geometric shape with up, down, concave, and convex sequences. Sample data from [32] were used in this case study. The actual values and the forecast values of three compared models are presented in Table 2. Again, the discrete grey model with fractional operators achieved the lowest MAPE compared to (where ) and direct modeling .

Case 3. The authors of [25] proposed a discrete grey model based on fractional order accumulation and studied the turnover of goods in Jiangsu Province, China. We consider the sample data from this paper as an example. The actual values and the forecast values of the different DGM models being compared are presented in Tables 3 and 4. The results show that the discrete grey model with fractional operators has a lower MAPE than the discrete grey model based on fractional order accumulate.

5. Prediction of China’s Sulfur Dioxide Emissions

According to China’s environmental status report published in 2016, from 2007 to 2015, the sulfur dioxide emission data (unit, million tons) are as follows:

The order of discrete grey model with fractional order operators can be resolved by the algorithm of particle swarm optimization. The minimum mean relative error is thus obtained. The results of discrete grey model with fractional order operators where , model (discrete grey model with fractional order operators where ), direct modeling (discrete grey model with fractional order operators where ), and are shown in Table 5. From Tables 5 and 6, we can see that the discrete grey model with fractional order operators where achieves the best fitting accuracy and test accuracy among the four models. The average relative error of discrete grey model with fractional order operators where is 1.118%. By using discrete grey model with fractional order operators where , the predicted values of sulfur dioxide emissions in China in the next five years are shown in Table 7.

According to the predicted results of China’s sulfur dioxide emissions shown in Table 7, the emissions show a steady decline in the five years from 2016 to 2020. By 2020, the amount of sulfur dioxide emissions is 16.145 million tons, only 86.843% of emissions in 2015. The result indicates that China’s current sulfur dioxide emission reduction policies have been well implemented and achieved a positive outcome.

6. Conclusion

In this paper, the discrete grey model with fractional operators is firstly studied and then used to predict China’s sulfur dioxide emissions. This paper introduces a new prediction model by changing and optimizing the values of fractional order, giving rise to a discrete grey model with fractional operators. The steps of modeling and genetic algorithms of order optimization for the proposed model are explained. Case studies show that the discrete grey model with fractional operators can achieve the best precision with optimized order. Then, China’s future sulfur dioxide emissions are predicted by the discrete grey model with fractional operators. Genetic algorithm can be used to determine the optimal order under the condition of minimum mean error. The results show that the proposed discrete grey model with fractional operators can obtain a higher fitting accuracy than the grey model , discrete grey model , and direct modeling . The forecast results show that China’s sulfur dioxide emissions will decrease steadily, indicating that China’s sulfur dioxide reduction policies are effective. It is noted that the fractional order accumulating generation operator and fractional order reducing generation operator may have the same expression. The unified expression of fractional order accumulating/reducing generation operators is an attractive topic in the grey system theory.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (71771033), Humanity and Social Science Foundation of the Ministry of Education of China (14YJAZH033), Project Foundation of Chongqing Municipal Education Committee (KJ1706166), and Research Fund of Chongqing Technology and Business University (670101485, KFJJ2017059).