Next Article in Journal
Extraction of Rice Bran Oil Using CO2-Expanded Hexane in the Two-Phase Region
Previous Article in Journal
Experimental Study on the Performance of a Dew-Point Evaporative Cooling System with a Nanoporous Membrane
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Method for the Estimation of Higher Heating Value of Municipal Solid Wastes

1
School of Management, China University of Mining and Technology (Beijing), Beijing 100083, China
2
Division of Solid Waste Management, School of Environment, Tsinghua University, Beijing 100084, China
3
Key Laboratory of Industrial Ecology and Environmental Engineering, School of Environmental Science & Technology, Dalian University of Technology, Dalian 116024, China
*
Authors to whom correspondence should be addressed.
Submission received: 24 February 2022 / Revised: 23 March 2022 / Accepted: 24 March 2022 / Published: 2 April 2022
(This article belongs to the Topic Challenge and Research Trends in Pyrolysis for Waste Raw Materials)
(This article belongs to the Section D1: Advanced Energy Materials)

Abstract

:
The measurement of the higher heating value (HHV) of municipal solid wastes (MSWs) plays a key role in the disposal process, especially via thermochemical approaches. An optimized multi-variate grey model (OBGM (1, N)) is introduced to forecast the MSWs’ HHV to high accuracy with sparse data. A total of 15 cities and MSW from the respective city were considered to develop and verify the multi-variant models. Results show that the most accurate model was POBGM (1, 5) of which the least error measured 5.41% MAPE (mean absolute percentage error). Ash, being a major component in MSW, is the most important factor affecting HHV, followed by volatiles, fixed carbon and water contents. Most data can be included by using the prediction interval (PI) method with 95% confidence intervals. In addition, the estimations indicated that the MAPE from estimating the HHV for various MSW samples, collected from various cities, were in the range of 3.06–34.50%, depending on the MSW sample.

1. Introduction

Municipal solid wastes (MSWs) are amalgamations of different solid waste streams produced as a result of domestic, commercial and institutional activities [1]. Owing to rapid economic growth, changing lifestyles and rapid urbanization, China’s MSWs production witnessed an exponential growth in the last decade. The annual generation of MSWs in Chinese urban areas reached ~0.23 billion tons in 2020 [2] and is expected to reach more than ~0.25 billion tons in 2022. The traditional methods for MSWs treatment include landfilling, anaerobic digestion and composting [3]. Landfilling was once commonly adopted due to its technological maturity and feasibility. However, in recent years, due to the unavailability of land and other associated environmental issues, landfilling has not been considered as a sustainable solution for the disposal of solid wastes [4]. While anaerobic digestion and composting, both being considered as technologically advanced and environmentally friendly techniques to recover energy and other valuable resources from MSWs, there are drawbacks, such as applicability only to biodegradable fractions of MSWs, long processing time and requirement of post-treatment for residues, which further hinder their adoptability [5].
Thermochemical treatments, such as incineration, pyrolysis and gasification, offer superiority over the biochemical conversion platforms, which include complete elimination of pathogens, very short processing time and complete conversion of carbon to fuels and other valuable chemicals [6]. Incineration is considered as the most reliable and economic route of waste to energy conversion in developing countries for reasons such as significant reduction in the volume of wastes to be disposed, processing time and better process control. For instance, the volume and mass of waste can be reduced by 80–90% and 70–80%, respectively [7]. The drawbacks, such as high capital cost and release of dioxins can cause detrimental effects on human and animal health, resulting in the development of emerging technologies, such as pyrolysis and gasification, to meet the various energy needs and carbon reduction strategy [8,9]. It is to be noted that the stable operation of biomass thermal conversion relies mostly on the characteristics of biomass, namely, the heating value, ash and water content and elemental composition [10]. ‘Higher’ heating value (HHV) or ‘lower’ heating value (LHV) is the common property of solid waste and fuel and can be experimentally determined by using an adiabatic oxygen bomb calorimeter. The availability and accessibility of the equipment may limit its usage, and, under such circumstances, researchers can use the proximate and ultimate analysis data to predict the HHV of fuel material by employing more practical models [11,12]. The HHV of feedstocks is the key parameter that determines the stability and efficiency of the thermochemical reactions [13].
The correlation between HHV and proximate/ultimate analysis is not new and dates back to late 1800s, with the first correlation proposed by Dulong based on ultimate analysis of coal [14]. Developing models for estimating the HHVs of MSW has attracted much attention owing to the benefits that are offered to handlers during solid waste management handling and disposal. Regression methods and back propagation neural network methods are the most common techniques used to predict the HHV of MSWs with simple composition based on known mechanisms or a mass of data [12,15]. As the thermal decomposition mechanism of the MSWs remains unclear, the HHV cannot be accurately predicted by the aforementioned models.
Grey model (GM), introduced by Deng in 1982, can be used to accurately predict the values of HHV under ambiguous conditions such as unknown mechanisms and scant data [16]. GM (1, 1) is a univariate model with high accuracy, which is most widely used to predict production under complicated situations, such as wind power generation [17] and MSW production forecasting [18]. The GM (1, 1) is used to characterize an unfamiliar system by employing a first-order differential equation. MSWs are characterized with complicated physical and chemical composition, which cannot be described by univariable model [19]. Therefore, to predict the HHV of the MSWs accurately, grey model with multivariable is more suitable [20]. Zhang et al. [21], Intharathirat et al. [22], and Duman et al. [23] studied the prediction of MSW production using the multi-variable GM (1, N) model.
To our knowledge, little or no studies have reported on the prediction of HHV for MSWs. To this end, proximate and elemental analysis could be used [12,15], and the newly developed GM (1, N) (N ≥ 2) model might offer a better method for HHV prediction as more variables are involved.
This study attempts to predict the HHV value of MSW collected from different locations, using correlation models (including 2 traditional regression methods and, GM (1, N) with different N) based on the proximate or ultimate analysis data. In addition, the practicality of multivariate grey prediction model GM (1, N), which employ grey rational analysis to estimate the HHV of MSW is studied in this work.

2. Data and Basic Methods

The flow path of the modelling program is shown in Figure 1 and is carried out as follows: input of the key factors and HHVs, determination of grey factors with Grey Relational Analysis (GRA) method, finding the best prediction model among the alternative 12 methods, and discussing the model influencing factors.

2.1. Data Sources

Shanxi Province, with a total area of 156,579 km2 has a population over 37.18 million as per 2018 census. The generation of MSWs in Shanxi Province has increased rapidly to 5 million tons in 2019 from 3.54 million tons since 2008. The characteristic data related to MSW from 15 cities of Shanxi Province were collected and averaged in time (3–6 points) within three months. The MSW in Shanxi Province is characterized by high moisture and ash content, and a relatively low content of volatile and fixed carbon, as shown in Table 1. The contents of moisture, volatile, ash and fixed carbon ranged from 14.79–39.46%, 8.70–20.99%, 40.08–65.72 and 2.71–10.80%, respectively.
The contents of C and O were higher than other elements (Table 2). The HHVs for all the samples were much lower than coal, and ranged between 764–1267 kcal·kg−1 [24]. Research on the HHV analysis and predictions prior to thermal disposal were critically important to avoid damage to equipment and ensure the operation under stable conditions [25,26]. In this work, the proximate and ultimate analysis results were used to develop the models for HHV prediction.

2.2. Traditional Forcasting Method

The performance of the established models is evaluated by the mean absolute percentage error (MAPE, %) as following [22].
M A P E ( % ) = ( k = 1 m | A k P k | A k × m ) × 100
MAPE can be divided into four categories as follows: Excellent (<10%), Good (10–20%), Reasonable (20–50%) and Incorrect (>50%) [27].
In this work, two traditional, regression prediction models based on proximate and ultimate analysis are introduced and compared to the proposed model, which is discussed in detail in Section 4.2.
Jigisha Parikh et al. [15] introduced a method with a MAPE of 3.74% to calculate the HHV (MJ·kg−1) of solid feedstock, and the correlation is given below:
H H V = 0.3636 × F C + 0.1559 × V M 0.0078 × A S H
where FC, VM, ASH denote fixed carbon, volatile content and ash content, respectively, in wt.% on a dry basis.
Boumanchar et al. [12] established a multiple regression and programming formalism method for HHV calculation based on ultimate analysis. The equation can be expressed as follows:
H H V = 2.775 + H + 0.004027 × C + 0.004027 × C 2 + 0.05706 H 12.97 + 0.02323 H 6.661 + 0.009398 H 5.961 + 12.97 H H 3 5.922 × C

2.3. Multivairant Grey Model, GM (1, N)

In this section, the multivariable grey model GM (1, N) is introduced based on the works of Zeng and Li [20].
GM (1, N) is a first-order grey model with N variables, including N-1 independent variables and one dependent variable. Let X 1 ( 0 ) be the dependent variable that has m initial sequences, as presented in Equation (4).
X 1 ( 0 ) = ( x 1 ( 0 ) ( 1 ) , x 1 ( 0 ) ( 2 ) , , x 1 ( 0 ) ( m ) )
N-1 independent variables denote X i ( 0 ) ( i = 2 , 3 , , N ) with m initial sequences are of strong correlation with X 1 ( 0 ) , as presented in Equation (5).
X i ( 0 ) = ( x i ( 0 ) ( 1 ) , x i ( 0 ) ( 2 ) , , x i ( 0 ) ( m ) )
Then, the first order accumulated generation operation (1-AGO) of the initial series, which is X j ( 1 ) , can be denoted as following Equation (6).
X j ( 1 ) = ( x j ( 0 ) ( 1 ) , x j ( 0 ) ( 2 ) , , x j ( 0 ) ( m ) ) ( j = 1 , 2 , , N )
where X j ( 1 ) can be calculated by X j ( 1 ) ( k ) = g = 1 k x i ( 0 ) ( g ) , k = 1 , 2 , m .  
Z 1 ( 1 ) is the mean sequence that is calculated by consecutive neighbour of the above X 1 ( 1 ) ,
Z 1 ( 1 ) = ( z 1 ( 1 ) ( 2 ) , z 1 ( 1 ) ( 3 ) , , z 1 ( 1 ) ( m ) )
where z 1 ( 1 ) ( k ) = 0.5 × ( x 1 ( 1 ) ( k ) + x 1 ( 1 ) ( k 1 ) ) , k = 2 , 3 , m .
Then,
x 1 ( 0 ) ( k ) + a z 1 ( 1 ) ( k ) = i = 2 N b i x i ( 1 ) ( k )
is the GM (1, N) model.
The parameters column can be presented as a ^ = [ a , b 1 , b 2 , b 3 , , b N ] T , which satisfies the least square estimate equation given below.
P = ( B T B ) 1 B T Y
where,
B = [ z 1 ( 1 ) ( 2 ) x 2 ( 1 ) ( 2 ) x N ( 1 ) ( 2 ) z 1 ( 1 ) ( 3 ) x 2 ( 1 ) ( 3 ) x N ( 1 ) ( 3 ) z 1 ( 1 ) ( m ) x 2 ( 1 ) ( m ) x N ( 1 ) ( m ) ]
Y = [ x 1 ( 0 ) ( 2 ) x 1 ( 0 ) ( 3 ) x 1 ( 0 ) ( m ) ]
Then, the desired grey parameters obtained in the last step can be substituted into Equation (12) to obtain the (k + 1) th 1-AGO result as following:
x ^ 1 ( 1 ) ( k + 1 ) = [ x 1 ( 0 ) 1 a i = 2 N b i x i ( 1 ) ( k + 1 ) ] e a k + 1 a i = 2 N b i x i ( 1 ) ( k + 1 )
Finally, the (k + 1) th predictive value, also x ^ 1 ( 0 ) ( k + 1 ) , can be obtained through the inverse accumulated generating operation using Equation (13)
x ^ 1 ( 0 ) ( k + 1 ) = a x ^ 1 ( 1 ) ( k + 1 ) = x ^ 1 ( 1 ) ( k + 1 ) x ^ 1 ( 1 ) ( k )
The aforementioned grey model has some drawbacks in practice, including mechanism, parameter and structure defects [20]. The optimization of current grey model is discussed in detail in the subsequent sections.

3. Optimized Grey Forecasting Model

3.1. OGM (1, N)

In this section, GM (1, N) optimized to obtain POGM (1, N) based on proximate analysis and EOGM (1, N) based on ultimate analysis was obtained.
The optimized grey model, OGM (1, N) is presented as follows:
x 1 ( 0 ) ( k ) + a z 1 ( 1 ) ( k ) = i = 2 N b i x i ( 1 ) ( k ) + h 1 × ( k 1 ) + h 2
where k = 2, 3, …, m. h 1 × ( k 1 ) and h 2 are linear correlation term and grey action quantity, respectively.
Then, the parameter series p ^ = [ b 1 , b 2 , b 3 , , b N , a , h 1 , h 2 ] T satisfies:
p ^ = ( B T B ) 1 B T Y
where,
B = [ x 2 ( 1 ) ( 2 ) x 3 ( 1 ) ( 2 ) x N ( 1 ) ( 2 ) z 1 ( 1 ) ( 2 ) 1 1 x 2 ( 1 ) ( 3 ) x 3 ( 1 ) ( 3 ) x N ( 1 ) ( 3 ) z 1 ( 1 ) ( 3 ) 2 1 x 2 ( 1 ) ( m ) x 3 ( 1 ) ( m ) x N ( 1 ) ( m ) z 1 ( 1 ) ( m ) m 1 1 ]
Y = [ x 1 ( 0 ) ( 2 ) x 1 ( 0 ) ( 3 ) x 1 ( 0 ) ( m ) ]
Then, the required grey parameters obtained in the last step can be substituted in Equation (18), and the result can be presented as,
x ^ 1 ( 1 ) ( k ) = t = 1 k 1 [ μ 1 i = 2 N μ 2 t 1 b i x i ( 1 ) ( k t + 1 ) ] + μ 2 k 1 x ^ 1 ( 1 ) ( 1 ) + j = 0 k 2 μ 2 j [ ( k j ) μ 3 + μ 4 ]
where, k = 2 , 3 , , m , μ 1 = 1 1 + 0.5 a , μ 2 = 1 0.5 a 1 + 0.5 a , μ 3 = h 1 1 + 0.5 a , μ 4 = h 2 h 1 1 + 0.5 a .
Finally, the inverse accumulated generating operation is introduced to calculate the (k + 1)th predictive value.
x ^ 1 ( 0 ) ( k ) = x ^ 1 ( 1 ) ( k ) x ^ 1 ( 1 ) ( k 1 ) k = 2 , 3 , , m

3.2. OBGM (1, N)

In the grey system, the background-value coefficient, ξ ( 0 < ξ < 1 ) , strangely equates to 0.5. In this study, the background-value coefficient will be further studied after the best grey model has been determined.
With parameter series p ^ = [ b 2 , b 3 , , b N , a , c , d ] T , the OGM (1, N) is presented as Equation (20).
x i ( 0 ) ( k ) + a ξ x 1 ( 1 ) ( k 1 ) = i = 2 N b i x i ( 1 ) ( k ) + k c + d ,   k = 2 , 3 , , m
However, many works showed that dynamic ξ can reduce the mean absolute prediction error (MAPE %). Therefore, the background-value coefficient ( ξ ) of the above OGM (1, N) is optimized using the particle swarm optimization algorithm (PSO) method [28]. With the optimized ξ , OBGM (1, N) is denoted. The introduction of this optimization method including the detailed MATLAB code is given in the work of Zeng and Li [20].
The final time response function and restored expression of the OBGM (1, N) model can be expressed as follows:
x ^ 1 ( 1 ) ( k ) = u = 1 k 1 [ υ 1 i = 2 N υ 2 u 1 b i x i ( 1 ) ( k u + 1 ) ] + υ = 0 k 2 υ 2 v [ ( k v ) υ 3 + υ 4 ] + υ 2 k 1 x ^ 1 ( 1 ) ( 1 ) k = 2 , 3 , , m
and
x ^ 1 ( 0 ) ( k ) = υ 1 i = 2 k [ i = 2 N ( υ 2 1 ) [ d 2 d ] υ 2 [ d 2 d ] [ d 3 d ] ( d 3 ) b i x i ( 1 ) ( k d + 2 ) ] + v = 0 k 3 [ υ 2 k 2 ( 2 υ 3 + υ 4 + 1 ) [ v k 3 ] + υ 2 v υ 3 ]
where k = 2 , 3 , , m , υ 1 = 1 1 + a ξ , υ 2 = 1 a 1 + a ξ , υ 3 = c 1 + a ξ , υ 4 = d 1 + a ξ .
The best multivariable grey model with optimized background-value coefficient will be denoted as POBGM (1, N) or EOBGM (1, N).

3.3. Grey Relational Analysis (GRA)

The method of grey relational analysis (GRA) was used to recognize the relational grade of selected factors in the aforementioned models and carried out on MATLAB software. Thus, the relationship can be determined between the reference series and compared series, which were denoted as x0 = (x0(1), x0(2), …, x0(k)) and xi = (xi(1), xi(2), …, xi(k)) (i = 1, 2,…, n), respectively. GRA was used to rank the selected factors, of which the one with the highest GRA grade (0~1.0) was the most important factor affecting the HHVs of MSWs. The GRA grades can be obtained using the equations listed in the work of Intharathirat et al. [22].

4. Results and Discussions

The relational grades of selected factors were firstly introduced. To identify the best prediction model, POGM (1, N) and EOGM (1, N) were studied with different N values and compared against the two traditional regression prediction methods mentioned in Section 2.2. Moreover, the coefficients ( ξ ) of the aforementioned best model were discussed. The results of using the best model to predict the HHVs of MSW is presented in the last section.

4.1. Identification of the Facotors

Two kinds of multivariable grey models are introduced in this study on the basis of the proximate and ultimate data of MSWs. In the GRA process, firstly, a series of HHVs and selected factors should be normalized to reduce the inaccuracy. In this work, the raw data are normalized by the mean value.
The grey relational grade (Table 3) shows that, in the proximate analysis related model, ash was the most influencing factor (0.7247), followed by volatile content (0.6361), fixed carbon (0.6157) and water content (0.5759). Soponpongpipat et al., [29] reported a prediction equation without the parameters of ash and mentioned that the ash content did not have a significant impact on the HHV of torrefied biomass. However, Qian et al. [11] reported that ash content had a significant impact on the HHV of biochars. Therefore, in view of the high ash content in the selected MSW feedstocks, the HHV of the MSW samples is expected to be greatly influenced by the ash content.
In the elemental analysis related model, C is the most influencing factor (0.8049), followed by N (0.7814), Cl (0.7617), H (0.7502), O (0.7218) and S (0.5292). The factors with values > 0.5 have a remarkable influence on HHVs [23]. Boumanchar et al. [12] reported that the change in HHV of MSW followed two trends: an increasing trend can be witnessed with the increase of carbon and hydrogen while a decreasing trend can be noticed with the increase in nitrogen and sulphur rates.

4.2. Comparative Analysis

The HHV series, including 10 data sets for establishing the model and five datasets for verification, were used to develop the alternative models, which include two multiple regression methods [12,15], and 10 multivariable models (including POGM (1, N) (N = 2, 3, 4, 5) which are based on the proximate analysis of MSWs, and EOGM (1, N) (N = 2,3,4,5,6,7) as based on the elemental analysis). The results of the 12 alternative models are presented in Table 4.
The results show that traditional models still bear good accuracy in forecasting HHVs of MSWs. Generally, POGM (1, N) exhibits better performance than EOGM (1, N) because the proximate analysis was obtained using a greater testing quantity than elemental analysis. It can be observed that multivariable grey models with larger N using proximate analysis data not only displays high modeling accuracy but also shows excellent performance in prediction. Similar analysis was reported in studies reported on multivariable grey models with larger N values for other areas [30,31,32]. Among the alternative models, POGM (1, 4) and POGM (1, 5) had similar performances and displayed excellent accuracy in forecasting HHV with the least error of 2.89% MAPE. Water content can be seen as a vital parameter which greatly affected the HHVs. Therefore, POGM (1, 5) that included water content was chosen to be used for further studies. The model parameters of multivariable grey models are presented in Table 5 (the values in bold font are the parameters of the best model).

4.3. HHV Forecasting

The background-value coefficient is usually noted as 0.5 in grey model for the simplification of modeling [33,34]. In this section, POBGM (1, 5) models with different background-value coefficients were established and their performances was compared in Table 6. The mathematical definitions in Table 6 are as follows: X1(0)—the actual HHV, x1(0)(k)—the predicted HHV, △(k)—the absolute error, re(k)—the relative error. It can be concluded that grey prediction models with background-value coefficient of 0.5 is rather accurate. Similar results were reported in literature [17]. For the following study, ξ = 0.5 and the model parameters are shown in Section 4.2.
The accuracies of both modeling and prediction steps using different sampling scales were determined in this part. As shown in Figure 2 and Table 7, it can be concluded that higher modeling accuracy can be obtained by using more volume of sample data.
In the predicting step, with a MAPE value of 3.0606%, using 10 samples for modeling can be considered as the most accurate scale. It is of great significance to obtain the best grey model using both modeling and prediction accuracy when establishing the grey model. A single number or point is often predicted which may not provide any necessary guidance in the actual application regarding their likely accuracy [22]. Prediction interval (PI) (95% confidence intervals) was used to forecast the upper and lower limit of HHVs of MSWs. The lower limit should meet the needs of stable operation during thermal disposal, while upper limit should fulfil the safety requirements [35]. As shown in Table 8 and Figure 3, the POBGM (1, 5) provides good modeling accuracy of HHVs in this study. The results show that four out of five forecasted data were in the interval of 95% PI.
Data Number 11, of which the HHV is 947.5 kcal kg−1, is outside the prediction curve and the prediction interval. This might be due to the measurement error or model deviation [36]. In the practical application of the established grey model, most of the data can be precisely forecasted with PI method. However, if the data bear a larger deviation, actual tests should be practiced. Normally, a specific mechanism is often proposed in advance. Callejón-Ferre et al. [37] used 20 models based on proximate or elemental analysis to predict the HHVs of crop residue, including univariate and multivariate models. These models belong to data fitting category and are adopted from other’s work. The composition and structure of organic material, especially MSW mentioned in this study, are quite complicated and cannot be simply described by fitting method. Large-scale data are often used to train the model to obtain a high prediction accuracy, e.g., neural network and other multivariate nonlinear regression methods [38]. The proposed grey model in this study can firstly filter the effect factors by using grey rational analysis, e.g., only factors with GRA larger than 0.5 will be considered in the model [17]. The larger the GRA degree, the more important the factor. Secondly, the grey model uses normalization and 1-AGO methods to reduce or even eliminate the effect of assuming mechanism [39]. More importantly, this novel method does not require large-scale data to train the model and also exhibits high accuracy even with scarce data.

5. Conclusions

The exponential growth of MSW generation in China requires highly efficient disposal on a large scale and rely significantly on thermal conversion platforms in particular. The accurate prediction of HHVs of MSWs is of great significance in the planning and management of these thermal and thermochemical operations. In this study, OBGM (1, N), on the basis of proximate or ultimate data, was developed to predict the HHV of MSWs. Thus, 12 models including two traditional regression methods, POGM (1, 2–5) and EOGM (1, 2–7) were studied. For multivariable models, GRA method is used to identify and quantify the influence of proximate and ultimate analysis. The selected best model is further studied by changing the data scale, optimizing the background-value coefficient and using prediction interval techniques.
Among the selected 12 models, POBGM (1, 5) with a background-value coefficient of 0.5 was observed as the best grey model and displayed a high accuracy, both in modeling and prediction step. According to the data scale study, less data led to higher modeling accuracy, but it cannot be concluded that less samples lead to high prediction accuracy. POBGM (1, 5) can ensure high prediction accuracy based on sparse data. With 95% PI, most of the data can be precisely forecasted. If the data bear a larger deviation, an actual test should be implemented. OBGM (1, N) models provided an effective modelling pathway for the prediction of HHV of MSWs under data scarcity conditions.

Author Contributions

The work was planned by M.Z., R.Z. and G.J. Initial models were created by W.D. and simulation was carried out by Z.C. and J.C. Manuscript was written by Z.C and revised by Z.J.T. All authors gave final approval for publication. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key Research and Development Program of China (2018YFC1901300). W.D. and R.Z. are grateful for support by the Major Science and Technology Projects in Shanxi Province [No. 20181102017].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All raw codes (including GRA analysis and grey model) required to run the model and simulation results are freely available online https://datadryad.org/stash/share/XEg68ZObgAwqN3_6lw7zQLfzxN8jFokJppZLSVPm-NA (accessed on 7 December 2021).

Acknowledgments

The authors also acknowledge the support on model optimization provided by Bo Zeng and his student of Chongqing Technology and Business University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, X.B.; Geng, Y. Municipal solid waste management in Dalian: Practices and challenges. Front. Env Sci Eng. 2012, 6, 540–548. [Google Scholar] [CrossRef]
  2. Yearbook National Bureau of Statistics of China. China Statistical Yearbook 2020; China Statistical Press: Beijing, China, 2021. [Google Scholar]
  3. Ali, J.; Rasheed, T.; Afreen, M.; Anwar, M.T.; Nawaz, Z.; Anwar, H.; Rizwan, K. Modalities for conversion of waste to energy—Challenges and perspectives. Sci. Total Environ. 2020, 727, 138610. [Google Scholar] [CrossRef] [PubMed]
  4. Nakhaei, M.; Amiri, V.; Rezaei, K.; Moosaei, F. An investigation of the potential environmental contamination from the leachate of the Rasht waste disposal site in Iran. Bull. Eng. Geol. Environ. 2015, 74, 233–246. [Google Scholar] [CrossRef]
  5. Saint-Joly, C.; Desbois, S.; Lotti, J.P. Determinant impact of waste collection and composition on anaerobic digestion performance: Industrial results. Water Sci. Technol. 2000, 41, 291–297. [Google Scholar] [CrossRef]
  6. Iqbal, A.; Liu, X.; Chen, G.H. Municipal solid waste: Review of best practices in application of life cycle assessment and sustainable management techniques. Sci. Total Environ. 2020, 729, 138622. [Google Scholar] [CrossRef]
  7. Arena, U. Process and technological aspects of municipal solid waste gasification. A review. Waste Manag. 2012, 32, 625–639. [Google Scholar] [CrossRef]
  8. Martínez, I.; Grasa, G.; Callén, M.S.; López, J.M.; Murillo, R. Optimised production of tailored syngas from municipal solid waste (MSW) by sorption-enhanced gasification. Chem. Eng. J. 2020, 401, 126067. [Google Scholar] [CrossRef]
  9. Lu, J.S.; Chang, Y.; Poon, C.S.; Lee, D.J. Slow pyrolysis of municipal solid waste (MSW): A review. Bioresour. Technol. 2020, 312, 123615. [Google Scholar] [CrossRef]
  10. Yin, C.-Y. Prediction of higher heating values of biomass from proximate and ultimate analyses. Fuel 2011, 90, 1128–1132. [Google Scholar] [CrossRef] [Green Version]
  11. Qian, C.; Li, Q.; Zhang, Z.; Wang, X.; Hu, J.; Cao, W. Prediction of higher heating values of biochar from proximate and ultimate analysis. Fuel 2020, 265, 116925. [Google Scholar] [CrossRef]
  12. Boumanchar, I.; Chhiti, Y.; M’hamdi Alaoui, F.E.; Sahibed-dine, A.; Bentiss, F.; Jama, C.; Bensitel, M. Municipal solid waste higher heating value prediction from ultimate analysis using multiple regression and genetic programming techniques. Waste Manag. Res. 2018, 37, 578–589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Vouvoudi, E.C.; Rousi, A.T.; Achilias, D.S. Thermal degradation characteristics and products obtained after pyrolysis of specific polymers found in Waste Electrical and Electronic Equipment. Front. Environ. Sci. Eng. 2017, 11, 1–10. [Google Scholar] [CrossRef]
  14. Maksimuk, Y.; Antonava, Z.; Krouk, V.; Korsakova, A.; Kursevich, V. Prediction of higher heating value based on elemental composition for lignin and other fuels. Fuel 2020, 263, 116727. [Google Scholar] [CrossRef]
  15. Parikh, J.; Channiwala, S.; Ghosal, G. A correlation for calculating HHV from proximate analysis of solid fuels. Fuel 2005, 84, 487–494. [Google Scholar] [CrossRef]
  16. Deng, J. Control problems of grey systems. Syst. Control Leit. 1982, 1, 288–294. [Google Scholar]
  17. Qian, W.; Wang, J. An improved seasonal GM(1,1) model based on the HP filter for forecasting wind power generation in China. Energy 2020, 209, 118499. [Google Scholar] [CrossRef]
  18. Xu, L.; Gao, P.; Cui, S.; Liu, C. A hybrid procedure for MSW generation forecasting at multiple time scales in Xiamen City, China. Waste Manag. 2013, 33, 1324–1331. [Google Scholar] [CrossRef]
  19. Zhou, H.; Meng, A.; Long, Y.; Li, Q.; Zhang, Y. An overview of characteristics of municipal solid waste fuel in China: Physical, chemical composition and heating value. Renew. Sustain. Energy Rev. 2014, 36, 107–122. [Google Scholar] [CrossRef]
  20. Zeng, B.; Li, C. Improved multi-variable grey forecasting model with a dynamic background-value coefficient and its application. Comput. Indust. Eng. 2018, 118, 278–290. [Google Scholar] [CrossRef]
  21. Zhang, Y. The Prediction of the Generation of Municipal Solid Waste Based on Grey Combination Model. Adv. Mater. Res. 2013, 807–809, 1479–1482. [Google Scholar] [CrossRef]
  22. Intharathirat, R.; Abdul Salam, P.; Kumar, S.; Untong, A. Forecasting of municipal solid waste quantity in a developing country using multivariate grey models. Waste Manag. 2015, 39, 3–14. [Google Scholar] [CrossRef] [PubMed]
  23. Duman, G.M.; Kongar, E.; Gupta, S.M. Estimation of electronic waste using optimized multivariate grey models. Waste Manag. 2019, 95, 241–249. [Google Scholar] [CrossRef] [PubMed]
  24. Mansur, F.Z.; Faizal, C.K.M.; Monir, M.U.; Samad, N.A.F.A.; Atnaw, S.M.; Sulaiman, S.A. Co-gasification between coal/sawdust and coal/wood pellet: A parametric study using response surface methodology. Int. J. Hydrog. Energy 2020, 45, 15963–15976. [Google Scholar] [CrossRef]
  25. Bagheri, M.; Esfilar, R.; Sina Golchi, M.; Kennedy, C.A. Towards a circular economy: A comprehensive study of higher heat values and emission potential of various municipal solid wastes. Waste Manag. 2020, 101, 210–221. [Google Scholar] [CrossRef] [PubMed]
  26. Chen, Y.-C.; Jhou, S.-Y. Integrating spent coffee grounds and silver skin as biofuels using torrefaction. Renew. Energy 2020, 148, 275–283. [Google Scholar] [CrossRef]
  27. Lewis, C.D. Industrial and Business Forecasting Methods; Butterworth Scientific: London, UK, 1982. [Google Scholar]
  28. Eberhart, R.; Kennedy, J. A new optimizer using particle swarm theory. In Proceedings of the Sixth International Symposium on Micro Machine and Human Science, Nagoya, Japan, 4–6 October 1995; pp. 39–43. [Google Scholar]
  29. Soponpongpipat, N.; Sittikul, D.; Sae-Ueng, U. Higher heating value prediction of torrefaction char produced from non-woody biomass. Front. Energy 2015, 9, 461–471. [Google Scholar] [CrossRef]
  30. Kayacan, E.; Ulutas, B.; Kaynak, O. Grey system theory-based models in time series prediction. Expert Syst. Appl. 2010, 37, 1784–1789. [Google Scholar] [CrossRef]
  31. Zhang, M.; Guo, H.; Sun, M.; Liu, S.; Forrest, J. A novel flexible grey multivariable model and its application in forecasting energy consumption in China. Energy 2022, 239, 122441. [Google Scholar] [CrossRef]
  32. Ma, X.; Liu, Z.; Wang, Y. Application of a novel nonlinear multivariate grey Bernoulli model to predict the tourist income of China. J. Comput. Appl. Math. 2019, 347, 84–94. [Google Scholar] [CrossRef]
  33. Guo, X.; Liu, S.; Wu, L.; Gao, Y.; Yang, Y. A multi-variable grey model with a self-memory component and its application on engineering prediction. Eng. Appl. Artif. Intell. 2015, 42, 82–93. [Google Scholar] [CrossRef]
  34. Ma, X.; Liu, Z.-B. The kernel-based nonlinear multivariate grey model. Appl. Math. Model. 2018, 56, 217–238. [Google Scholar] [CrossRef]
  35. Qian, X.; Lee, S.; Soto, A.-M.; Chen, G. Regression Model to Predict the Higher Heating Value of Poultry Waste from Proximate Analysis. Resources 2018, 7, 39. [Google Scholar] [CrossRef] [Green Version]
  36. Razaque, A.; Amsaad, F.; Hariri, S.; Almasri, M.; Rizvi, S.S.; Frej, M.B.H. Enhanced Grey Risk Assessment Model for Support of Cloud Service Provider. IEEE Access 2020, 8, 80812–80826. [Google Scholar] [CrossRef]
  37. Callejón-Ferre, A.J.; Velázquez-Martí, B.; López-Martínez, J.A.; Manzano-Agugliaro, F. Greenhouse crop residues: Energy potential and models for the prediction of their higher heating value. Renew. Sustain. Energy Rev. 2011, 15, 948–955. [Google Scholar] [CrossRef]
  38. Dashti, A.; Noushabadi, A.S.; Asadi, J.; Raji, M.; Chofreh, A.G.; Klemeš, J.J.; Mohammadi, A.H. Review of higher heating value of municipal solid waste based on analysis and smart modelling. Renew. Sustain. Energy Rev. 2021, 151, 111591. [Google Scholar] [CrossRef]
  39. Ren, J. GM(1,N) method for the prediction of anaerobic digestion system and sensitivity analysis of influential factors. Bioresour. Technol. 2018, 247, 1258–1261. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the modelling flow.
Figure 1. Schematic diagram of the modelling flow.
Energies 15 02593 g001
Figure 2. Accuracy comparison of different sample scale.
Figure 2. Accuracy comparison of different sample scale.
Energies 15 02593 g002
Figure 3. Forecasting of MSW’s HHVs using POBGM (1, 5) at 95% PI.
Figure 3. Forecasting of MSW’s HHVs using POBGM (1, 5) at 95% PI.
Energies 15 02593 g003
Table 1. The HHVs and proximate analysis of MSWs of 15 cities in Shanxi Province.
Table 1. The HHVs and proximate analysis of MSWs of 15 cities in Shanxi Province.
Num.CityHHVProximate Analysis
kcal·kg−1Mad/%Vad/%Aad/%FCad/%
1Taiyuan76434.4114.8648.042.71
2Datong125823.0420.0850.746.14
3Yangquan980.539.4612.5541.846.17
4Changzhi135315.9112.2661.0410.8
5Jincheng139322.512.6454.310.59
6Shuozhou86835.3820.9940.083.55
7Xinzhou100425.3614.2856.483.89
8Jinzhong121223.5213.5254.948.01
9Linfen104625.1314.9254.915.05
10Yuncheng91918.628.765.726.96
11Lvlinag947.524.0111.4359.594.98
12Yuanping96821.4113.3360.564.72
13Dingxiang94723.6812.858.834.7
14Fenyang125914.7911.4364.259.53
15Xiaoyi126725.314.9552.097.67
HHV—Higher heating value, air dried basis, kcal/kg; Mad—Water content, air dried basis, wt.%; Vad—Volatile content, air dried basis, wt.%; Aad—Ash content, air dried basis, wt.%; FCad—Fixed carbon, air dried basis, wt.%.
Table 2. The ultimate analysis of MSWs of 15 cities.
Table 2. The ultimate analysis of MSWs of 15 cities.
City Num.Ultimate Analysis (wt.%, Air Dried Basis)
CHONSCl
18.611.097.000.290.590.16
215.341.308.890.440.250.20
310.610.916.390.310.510.18
413.281.317.520.450.510.13
515.261.215.910.390.450.18
614.471.168.280.340.290.14
79.801.165.980.370.870.17
812.081.177.030.360.830.19
910.861.126.970.310.710.18
108.281.235.060.400.690.17
118.081.126.600.360.250.18
129.601.256.110.390.700.16
139.151.206.360.410.390.13
1412.011.276.460.350.900.16
1512.951.236.620.391.430.17
Table 3. Grey relational grade of selected factors.
Table 3. Grey relational grade of selected factors.
Proximate Analysis
Water Content Volatile Content Ash Fixed Carbon
0.57590.63610.72470.6157
Elemental Analysis
CHONSCl
0.80490.75020.72180.78140.52920.7617
Table 4. Comparison of the selected models.
Table 4. Comparison of the selected models.
Models VariableModelling
MAPE * (%)
AccuracyPrediction
MAPE (%)
Traditional Models
1Proximate analysis12.634←G G→15.45
2Elemental analysis11.99←G E→7.45
Multivariate Grey Prediction Models
POBGM (1, N)
(1,2)Ash61.04←I I→1555.00
(1,3)Ash, Vad40.40←R I→887.35
(1,4)Ash, Vad, FC5.42←E E→2.89
(1,5)Ash, Vad, FC, MC5.41←E E→3.06
EOBGM (1, N)
(1,2)C16.39←G I→79.16
(1,3)C N11.12←G R→25.16
(1,4)C N Cl9.27←E G→15.98
(1,5)C N Cl H7.72←E G→14.64
(1,6)C N Cl H O4.90←E R→38.33
(1,7)C N Cl H O S0.0026←E I→123,173.23
* MAPE—Mean absolute percentage error, %. E—Excellent. G—Good. R—Reasonable. I—Incorrect. The best fit model is in bold.
Table 5. The model parameters of multivariable grey models.
Table 5. The model parameters of multivariable grey models.
TypeParameters
Nacdb2b3b4b5b6b7
POBGM2−0.54−1532.211951.2116.74
3−0.521348.741928.1015.14−4.6621
42.23152.52−153.498.9156.8333152.37
52.23145,846.98153.491468.461517.39701613.021460.57
2−0.47−331.901659.18−18.24
EOBGM3−0.11−1210.681414.10−27.483816.6772
40.16−3087.301095.488.865543.44566680.63
50.70−2960.29578.0838.628572.74969620.30−1270.18
6−0.05−7208.66267.51−42.3313306.570516865.80−1863.00319.64
7−1.45−2395.253576.3523.905903.7691−6479.40−1631.2715.001907.87
Table 6. POBGM (1,5) model results with different background-value coefficients (0.05–0.95 with 0.1 interval).
Table 6. POBGM (1,5) model results with different background-value coefficients (0.05–0.95 with 0.1 interval).
X1(0)(a): ξ = 0.05(b): ξ = 0.15(c): ξ = 0.25(d): ξ=0.35
x1(0)(k)△(k)re(k)x1(0)(k)△(k)re(k)x1(0)(k)△(k)re(k)x1(0)(k)△(k)re(k)
764.00.764.000.000.00764.000.000.00764.000.000.00764.000.000.00
1258.001377.97119.979.541235.77−22.231.771139.48−118.529.423877.552619.55208.23
980.50983.653.150.32959.05−21.452.19927.54−52.965.401915.93935.4395.40
1353.001499.66146.6610.841414.4561.454.541363.2010.200.752737.341384.34102.32
1393.001333.11−59.894.301276.75−116.258.351246.21−146.7910.542593.841200.8486.21
868.001026.27158.2718.23972.35104.3512.02944.7876.788.852080.161212.16139.65
1004.00988.13−15.871.58933.73−70.277.00874.37−129.6312.911982.58978.5897.47
1212.001298.8486.847.171230.7818.781.551184.37−27.632.282457.401245.40102.76
1046.001018.73−27.272.61968.94−77.067.37926.16−119.8511.462069.771023.7797.88
919.001097.05178.0519.371038.42119.4212.99971.2152.215.682159.221240.22134.95
8.22 6.42 7.48 118.32
(e): ξ = 0.45(f): ξ = 0.5(i): ξ = 0.55(j): ξ = 0.65
x1(0)(k)△(k)re(k)x1(0)(k)△(k)re(k)x1(0)(k)△(k)re(k)x1(0)(k)△(k)re(k)
764.00.0.000.00764.000.000.00764.000.000.00764.000.000.00
1674.55416.5533.111267.459.450.751438.32180.3214.331893.00635.0050.48
1131.13150.6315.36969.79−10.711.091030.4349.935.091523.77543.2755.41
1595.18242.1817.901383.8930.892.281463.41110.418.161810.90457.9033.84
1512.60119.608.591330.65−62.354.481424.1531.152.241855.74462.7433.22
1166.78298.7834.42998.30130.3015.011080.36212.3624.471526.14658.1475.82
1119.62115.6211.52921.37−82.638.231012.708.700.871415.38411.3840.98
1401.58189.5815.641196.81−15.191.251264.8652.864.361585.71373.7130.83
1189.52143.5213.721005.70−40.303.851105.6359.635.701505.64459.6443.94
1234.31315.3134.311026.91107.9111.741128.04209.0422.751542.45623.4567.84
20.51 5.41 9.77 48.04
(k): ξ = 0.75(l): ξ = 0.85(m): ξ = 0.95(n): ξ= 0.500667307066868
x1(0)(k)△(k)re(k)x1(0)(k)△(k)re(k)x1(0)(k)△(k)re(k)x1(0)(k)△(k)re(k)
764.000.000.00764.000.000.00764.000.000.00764.000.000.00
1977.68719.6857.211071.99−186.0114.79944.89−313.1124.891258.210.210.02
3822.852842.35289.891154.71174.2117.771113.45132.9513.56980.950.450.05
9909.008556.00632.371376.2523.251.721281.88−71.125.261371.5618.561.37
28,801.3627,408.361967.581025.45−367.5526.391152.02−240.9817.301326.27−66.734.79
87,582.7786,714.779990.181284.83416.8348.02966.7998.7911.38997.92129.9214.97
272,022.56271,018.5626,993.88431.28−572.7257.04818.97−185.0318.43915.38−88.628.83
849,023.53847,811.5369,951.452457.821245.82102.791376.55164.5513.581191.19−20.811.72
2,654,524.812,653,478.81253,678.66−1512.57−2558.57244.61565.81−480.1945.91999.86−46.144.41
8,304,352.458,303,433.45903,529.215555.044636.04504.471262.34343.3437.361019.57100.5710.94
140,787.83 113.06 20.85 5.23
Table 7. Comparison of POBGM (1, 5) model prediction accuracy when using a different sample scale.
Table 7. Comparison of POBGM (1, 5) model prediction accuracy when using a different sample scale.
Sample ScaleActual DataPredicted ValueAbsolute DeviationRelative Deviation
kcal·kg−1kcal·kg−1kcal·kg−1%
Eight samples (Modeling MAPE % = 12.0145)
9104695.2940−950.706090.8897
109191340.7320421.732045.8903
11947.51262.7320315.232033.2699
129681231.4700263.470027.2180
139471186.2230239.223025.2611
1412591427.3070168.307013.3683
1512671196.3640−70.63605.5751
Prediction MAPE %34.4960
Nine samples (Modeling MAPE % = 2.6756)
109191522.7460603.746065.6960
11947.51157.5170210.017022.1654
129681145.4470177.447018.3313
139471118.0910171.091018.0666
1412591412.0670153.067012.1578
1512671233.406033.59402.6515
Prediction MAPE %23.1781
Ten samples (Modeling MAPE % = 5.4103)
11947.51000.692053.19205.6139
12968973.91405.91400.6110
13947949.96702.96700.3133
1412591283.238024.23801.9252
1512671180.345086.65506.8394
Prediction MAPE %3.0606
Eleven samples (Modeling MAPE % = 6.4045)
129681565.0840597.084061.6822
13947992.207045.20704.7737
1412591333.711074.71105.9342
1512671228.908038.09203.0065
Prediction MAPE %18.8491
Twelve samples (Modeling MAPE % = 5.6208)
13947683.1850263.815038.6155
1412591246.993012.00700.9629
1512671165.4240101.57608.7158
Prediction MAPE %16.0980
Table 8. Prediction interval of MSW’s HHVs at 95% PI (kcal·kg−1).
Table 8. Prediction interval of MSW’s HHVs at 95% PI (kcal·kg−1).
No.Actual DataPOBGM (1, 5) Model
Prediction DataLower LimitUpper Limit
11947.501000.69For Testing
12968.00973.91963.80984.03
13947.00949.97927.71972.22
141259.001283.241261.991304.48
151267.001180.351090.311270.38
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dong, W.; Chen, Z.; Chen, J.; Ting, Z.J.; Zhang, R.; Ji, G.; Zhao, M. A Novel Method for the Estimation of Higher Heating Value of Municipal Solid Wastes. Energies 2022, 15, 2593. https://0-doi-org.brum.beds.ac.uk/10.3390/en15072593

AMA Style

Dong W, Chen Z, Chen J, Ting ZJ, Zhang R, Ji G, Zhao M. A Novel Method for the Estimation of Higher Heating Value of Municipal Solid Wastes. Energies. 2022; 15(7):2593. https://0-doi-org.brum.beds.ac.uk/10.3390/en15072593

Chicago/Turabian Style

Dong, Weiguo, Zhiwen Chen, Jiacong Chen, Zhao Jia Ting, Rui Zhang, Guozhao Ji, and Ming Zhao. 2022. "A Novel Method for the Estimation of Higher Heating Value of Municipal Solid Wastes" Energies 15, no. 7: 2593. https://0-doi-org.brum.beds.ac.uk/10.3390/en15072593

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop