Journal of Fuels

Volume 2014 (2014), Article ID 392698, 11 pages

http://dx.doi.org/10.1155/2014/392698

## Multivariable Regression and Adaptive Neurofuzzy Inference System Predictions of Ash Fusion Temperatures Using Ash Chemical Composition of US Coals

^{1}Department of Mining Engineering, Tehran Science and Research Branch, Islamic Azad University, Tehran 1477893855, Iran^{2}Young Researchers and Elites Club, Tehran Science and Research Branch, Islamic Azad University, Tehran 1477893855, Iran

Received 12 September 2013; Accepted 16 April 2014; Published 22 May 2014

Academic Editor: Kaustubha Mohanty

Copyright © 2014 Shahab Karimi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

In this study, the effects of ratios of dolomite, base/acid, silica, SiO_{2}/Al_{2}O_{3}, and Fe_{2}O_{3}/CaO, base and acid oxides, and 11 oxides (SiO_{2}, Al_{2}O_{3}, CaO, MgO, MnO, Na_{2}O, K_{2}O, Fe_{2}O_{3}, TiO_{2}, P_{2}O_{5}, and SO_{3}) on ash fusion temperatures for 1040 US coal samples from 12 states were evaluated using regression and adaptive neurofuzzy inference system (ANFIS) methods. Different combinations of independent variables were examined to predict ash fusion temperatures in the multivariable procedure. The combination of the “11 oxides + (Base/Acid) + Silica ratio” was the best predictor. Correlation coefficients () of 0.891, 0.917, and 0.94 were achieved using nonlinear equations for the prediction of initial deformation temperature (IDT), softening temperature (ST), and fluid temperature (FT), respectively. The mentioned “best predictor” was used as input to the ANFIS system as well, and the correlation coefficients () of the prediction were enhanced to 0.97, 0.98, and 0.99 for IDT, ST, and FT, respectively. The prediction precision that was achieved in this work exceeded that reported in previously published works.

#### 1. Introduction

Usually, coal quality for coal used for the generation of electricity refers to differences in heating value, grindability, sulfur content, and ash fusion characteristics. The potential for slagging on furnace walls is related to the ash composition and temperature. Ash fusibility is a factor that can be used to determine the performance of coals related to slagging [1].

ASTM Standard D1857, which specifies the experimental method for determining ash fusion temperatures (AFTs), is based on the gradual thermal deformation of a pyramid-shaped ash sample in either an oxidizing or reducing atmosphere. The test results can be reported for four temperatures: (1) the initial deformation temperature (IDT), which is the temperature at which the first rounding of the apex of the cone occurs and the pyramid begins to demonstrate evidence of deformation; (2) the softening temperature (ST) or fusion temperature (FT), which is the temperature at which the cone has fused and the height equals the width; (3) the hemispherical temperature (HT), which is the temperature at which the cone has fused into a hemispherical lump and the height is equal to half of the width at the base; and (4) the fluid temperature (FT), which is the temperature at which the fused mass has spread out in a nearly flat layer [2].

The temperature difference between the initial deformation temperature and the fluid temperature gives information on the type of deposit to be expected on the surfaces of the furnace tubes [3]. A small temperature difference indicates a thin and difficult-to-remove slag, whereas greater temperature differences indicate less adhesive deposits. The softening temperature is a criterion used to indicate the ease with which deposits can be removed from heat-transfer surfaces. If ash particles arrive at heat-absorbing surfaces at temperatures below their softening temperature, ash removal is relatively easy. If the ash particles arrive at these surfaces at temperatures above their softening temperature, the resulting deposit is more difficult to remove [3].

Although AFTs can be used to predict the slagging and fouling properties of coals, the reproduction of AFTs for a given coal in different laboratories may differ by °C [4–6].

Many researchers have tried to predict ash fusion temperatures based on the chemical composition of the coal [4, 7–11]. If reliable, correlative results could be established, and they could serve as a guide for coal combustion operation and eliminate the need for duplicate laboratory tests. Winegartner and Rhodes (1975) used regression analysis in an attempt to correlate ash fusibility with its chemical composition for a large number of ash samples from coal in United States (eastern and western coal), but their correlations were poor when they were applied to other coal ash samples from different sources [4]. Gray (1987) predicted ash fusion temperatures from the ash compositions of some New Zealand coals using three techniques of ternary equilibrium phase diagrams, stepwise regression analyses, and multiple regression analyses; these techniques were capable of predicting ash fusion temperatures to within 40°K [7]. Seggiani (1999) used a database that included about 300 ash samples from coals of different sources and biomasses. The proposed correlations predicted ash fusion temperatures and the temperatures of critical viscosity with a correlation () of 0.8. It was also mentioned that a small number of samples with similar compositions can give good results by the statistical approach, but the application to a large number of samples from different sources and with a wide range of compositions can decrease the effectiveness of the method [8].

Lolja et al. (2002) analyzed the ash fusion temperatures for 17 Albanian coals using oxides analysis from various perspectives, such as discrete species, acids and bases, crystal components, fluxing agents, and cement constituents, and in accordance with the periodic table. It was shown that deformation temperature and flow temperature correlated well with basic and acidic oxides, respectively. Also, AFTs decreased as basic oxide content increased [9].

G. Özbayoğlu and M. E. Özbayoğlu (2006) used the chemical composition of ash (eight oxides) and coal parameters, such as ash content, specific gravity, Hardgrove index, and mineral matter content, to predict AFTs and showed that nonlinear correlations are better than linear correlations for estimating ash fusion temperatures [11].

Liu et al. (2007) used the Ant Colony Optimization-Back Propagation (ACO-BP) neural network approach to model coal ash fusion temperatures based on the chemical compositions of the ashes. Data on 80 typical Chinese coal ash samples were used for training and testing. The results showed that the ACO-BP neural network can obtain better performance than empirical formulas and BP neural networks [12].

The aim of this research is to provide new correlations for the prediction of AFTs by using an extensive coal ash composition database. The specific questions we sought to answer in this research are listed below.(i)Is it possible to predict AFTs that have acceptable correlations using an extensive database that consists of data from the chemical analyses of 1040 coal ash samples from 12 states in the USA?(ii)Which variables have significant positive and negative effects?(iii)Which combination of predictive indices is the best for predicting AFTs for a large number of samples from different sources?(v)Is it possible to achieve better correlations and decrease errors using adaptive neurofuzzy inference system (ANFIS) modeling?

#### 2. Experimental

The database that was used to test the proposed approaches was obtained from the U.S. Geological Survey Coal Quality (COALQUAL) project of the USGS. A total of 1040 sets of coal ash analyses were used. The database, including the major and minor constituents from the oxide analyses, that is, SiO_{2}, Al_{2}O_{3}, CaO, MgO, MnO, Na_{2}O, K_{2}O, Fe_{2}O_{3}, TiO_{2}, P_{2}O_{5}, and SO_{3}, as well as ash fusion temperatures, is also available in the electronic Appendix to this paper.

The sampling procedures and the chemical analytical methods used can be found at the following web address: http://energy.er.usgs.gov/products/databases/CoalQual/index.htm. The number of samples and the range of AFTs for different states are shown in Table 1.

#### 3. Methods

##### 3.1. Input Variables

Table 2 shows the range of input variables, which consists of 11 oxides, base, and acid, as well as base/acid, dolomite, silica, SiO_{2}/Al_{2}O_{3}, and Fe_{2}O_{3}/CaO ratios, that was used to predict AFTs; these are the variables that were commonly used in previously published works [4, 7–12]. Each variable was examined separately and in combination with the other variables. The input variables in separate forms are as follows:

##### 3.2. Regression

Linear regression estimates the coefficients of the linear equation, involving one or more independent variables, to have a reliable prediction of the value of the dependent variable. All variables must pass the tolerance criterion to be entered in the equation, regardless of the entry method specified. The default tolerance level is 0.0001. Also, a variable is not entered if it would cause the tolerance of another variable already in the model to drop below the tolerance criterion. All independent variables selected are added to a single regression model. However, different entry methods can be specified for different subsets of variables. Method selection allows specifying how independent variables are entered into the analysis. Using different methods, a variety of regression models can be selected from the same set of variables [13].

Nonlinear regression is a method of finding a nonlinear model of the relationship between the dependent variable and a set of independent variables. Unlike traditional linear regression, which is restricted to estimating linear models, nonlinear regression can estimate models with arbitrary relationships between independent and dependent variables. This is accomplished using iterative estimation algorithms [13].

In this study, both linear and nonlinear regression analyses were used to develop correlations between the chemical compositions of various ashes and their fusion temperatures by using the stepwise regression in the SPSS software package. A stepwise procedure for selecting variables was used, and the variables were entered sequentially into the model. The first variable considered for use in the equation was the one with the largest positive or negative correlation with the dependent variable. This variable was entered into the equation only if it satisfied the criterion for entry. The next variable, with the largest partial correlation, was considered as the second input to the equation. The procedure stops when there are no variables that meet the entry criterion [13].

##### 3.3. Adaptive Neurofuzzy Inference System (ANFIS)

In the artificial intelligence field, the term “neurofuzzy” refers to combinations of artificial neural networks and fuzzy logic. Fuzzy modeling and neural networks have been recognized as powerful tools that can facilitate the effective development of models and integrate information from different sources, such as empirical models, physical laws, or measurements and heuristics [14]; these two tools were combined in order to achieve readability and learning ability at the same time [15].

The neurofuzzy in the fuzzy modeling research field is divided into two areas: linguistic fuzzy modeling that is focused on interpretability, mainly the Mamdani model, and precise fuzzy modeling that is focused on accuracy, mainly the Takagi-Sugeno-Kang (TSK) model [16].

ANFIS is an architecture that is functionally equivalent to a Takagi-Sugeno-Kang-type fuzzy rule base [17]; it is a class of adaptive, multilayer, feed-forward networks that is functionally equivalent to a fuzzy inference system.

A fuzzy rule in a Sugeno fuzzy model has the form of: if is and is , then ,where and are input fuzzy set in antecedent and, usually, is a zero- or first-order polynomial function in the consequent. The fuzzy reasoning procedure for the first-order Sugeno fuzzy model and equivalent ANFIS structure is shown in Figure 1.

Here, the defuzzification procedure in the Mamdani fuzzy model is replaced by the operation of the weighted average in order to avoid the time-consuming procedure of defuzzification. Defuzzification refers to the way a crisp value is extracted from a fuzzy set as a representative value [17].

More details about the ANFIS architecture, learning algorithm, and training methods were described by [15, 17].

#### 4. Results and Discussion

##### 4.1. Relationships between AFTs and Individual Input Variables

By the least squares mathematical method, the correlation coefficients () of SiO_{2}, Al_{2}O_{3}, CaO, MgO, MnO, Na_{2}O, K_{2}O, Fe_{2}O_{3}, TiO_{2}, P_{2}O_{5}, SO_{3}, and base and acid oxides, as well as base/acid, dolomite, silica, SiO_{2}/Al_{2}O_{3}, and Fe_{2}O_{3}/CaO ratios for AFTs were determined. Table 3 shows the relationship between individual input variables and AFTs in both linear and nonlinear (quadratic and cubic) procedures. Exponential, logarithmic, logistic, inverse, compound, power, and growth nonlinear procedures were examined as well; the results were not better than the results obtained using quadratic and cubic procedures.

It can be seen that the correlations of SiO_{2}, Al_{2}O_{3}, Fe_{2}O_{3}, base, and acid, as well as base/acid and silica ratios, with AFTs are significant. The other parameters are not significant. Relationships of the above-mentioned significant variables with IDT, which is representative of AFTs, are shown in Figure 2. For the ST and FT temperatures, the effects of the variables were similar. According to Table 3 and Figure 2, it can be seen that the correlations between the AFTs and significant variables should be nonlinear in the prediction of AFTs. As a general result in IDT, ST, and FT predictions, the variables of acid oxide, silica ratio, Al_{2}O_{3}, and SiO_{2} are the most significant variables, in that order, with positive effects on AFTs, and the variables of base oxide, base/acid, and Fe_{2}O_{3} are the most significant variables, with negative effects on AFTs.

##### 4.2. Multivariable Regression

The multivariable procedure was used to evaluate the possibility of achieving better prediction accuracy than the single-variable procedure. Different combinations of the most significant variables (Table 3 and Figure 2) were examined. Finally, the following eight combinations were determined as the suitable predictors:

Table 4 shows the correlation coefficient, standard deviation, mean, and residual (% difference between actual and predicted AFTs) less than 120°F for each of the inputs (see (9) to (16)) in both linear and nonlinear (quadratic) situations. The cubic equations did not improve the accuracy significantly.

The results show that input 8 (see (16)) was the best predictor based on the correlation coefficient and a residual of less than 120°F in the nonlinear condition. Equations (17) to (19) show the predictions of AFTs in nonlinear forms using input 8 in a multivariable, stepwise, regression procedure:

The distributions of the differences between values predicted from (17) to (19) and the actual determined amounts of IDT, ST, and FT are shown in Figures 3, 4, and 5, respectively.

##### 4.3. ANFIS Prediction

Input 8, the best predictor in the multivariable regression procedure, was used to determine whether ANFIS is able to predict AFTs more accurately than regression. The determination was conducted by using the MATLAB software package, ANFIS menu, to identify the relationships between the “11 oxide + (Base/Acid) + Sr” input variable and IDT, ST, and FT.

The first step in a neurofuzzy inference system is the determination of system input and output. For this study, input 8 was used as the input to the system. The outputs of the system include three temperatures, IDT, ST, and FT, and each of these outputs has been studied separately.

The fuzzy inference system used in this research was the Sugeno system. The output functions in the Sugeno system are linear or constant. A rule in the fuzzy Sugeno model is:

In the Sugeno system, for a zero-order model, the plane is constant (). The plane of output of any rule is weighted by . The system’s final output is the weighted average of all outputs, which is calculated as follows:

The subtractive clustering scheme was used to cluster data. The subtractive clustering algorithm is a fast, one-pass algorithm for estimating the number of clusters and the cluster centers in a set of data. Subtractive clustering is based on densitometry in space. In this work, the best designed neurofuzzy system for each of the three temperatures was a system with three clusters that could predict the AFTs with the highest correlation coefficients. For IDT, we selected range of influence, squash factor, accept ratio, and reject ratio as 0.5, 1.6, 0.5, and 0.15, respectively. For ST, we selected the values for these same four factors as 0.5, 1.5, 0.5, and 0.15, respectively. For FT, we selected the values of 0.4, 1.5, 0.5, and 0.15, respectively. The structure of the model based on this neurofuzzy system designed for prediction of FT, as an example, is shown in Figure 6. Gaussian membership functions were used in this work. Each of the input variables was assigned to one of three linguistic classes, that is, low, medium, and high. For training of the ANFIS, the hybrid method was used. The data were divided into two subsets. We used 780 sets of data for training, and the remaining 260 sets were used for testing. For the training stage, we selected 200 epochs. After the training was completed, specific rules for predicting each temperature were obtained by the software.

The prediction results achieved by the nurofuzzy system are very promising. The correlation coefficient (), standard deviation (S.D.) of the residual, and residual (% difference between actual and predicted AFTs) less than 120°F (°F) for neurofuzzy and nonlinear regression (quadratic) are shown in Table 5. It can be seen that the predictions of AFTs for ANFIS are significantly better than those from multivariable regression.

The IDT estimated by ANFIS in the testing stage compared to actual determined values are shown in Figure 7. The distribution of the differences between actual and estimated IDT is shown in Figure 8.

Figure 9 shows ST estimated by ANFIS compared to actual determined values. The distribution of the differences between actual and estimated ST is shown in Figure 10.

FT estimated by ANFIS in the testing stage compared to actual determined values is shown in Figure 11. The distribution of differences between actual and estimated FT is shown in Figure 12.

##### 4.4. Technical Considerations

Table 4 shows that nonlinear regression analysis produced more accurate results than linear analysis, with regard to correlation coefficient and percent of residual Re less than 120°F (49°C). An value of 0.7 is generally acceptable, 0.8 is good, and 0.9 or higher is excellent [8]. In the present study, it can be seen that, for nonlinear correlations (see (17) to (19)), values reached 0.943 ( of 0.891), 0.957 ( of 0.917), and 0.97 ( of 0.94) for the prediction of IDT, ST, and FT, respectively; whereas, the ANFIS prediction improved the values to 0.985, 0.99, and 0.995 for the prediction of IDT, ST, and FT, respectively (Table 5).

As mentioned before, test results from different laboratories for the AFT of a given coal may differ by 20–100°C [4–6]; therefore, the residual can be used as a criterion for the comparison of prediction error and common experimental errors in the laboratory for the regression and ANFIS approaches. In the current work, the nonlinear correlations (see (17) to (19)) can estimate IDT, ST, and FT with a difference less than 120°F (49°C) for about 83, 90, and 98%, respectively. The standard deviations for these nonlinear equations are between 56–93°F (13–34°C) (Input 8, Table 4). The differences less than 120°F for ANFIS prediction of IDT, ST, and FT are 99, 100, and 100%, respectively (Table 5); the standard deviation is between 23 and 45°F (−5 to 7°C) as well (Table 5). These fall within the experimental errors of the measurements. The above-mentioned results indicate that the final proposed nonlinear correlations and ANFIS model were satisfactorily successful; however, the ANFIS prediction is much better than nonlinear correlations.

In other works related to this study, G. Özbayoğlu and M. E. Özbayoğlu (2006) applied linear and nonlinear regression analyses for 30 coal samples collected from different parts of Turkey to develop correlations between the chemical composition of ash, ash content, mineral matter content, specific gravity, and Hardgrove index and ash fusion temperatures; the accuracies of of 0.955, 0.963, and 0.934 were archived for ST, MT (melting temperature), and FT, respectively [11]. Lolja et al. (2002) predicted AFTs for 17 Albanian coals using the chemical compositions of various coal ashes [9]. Seggiani (1999) used a database that included about 300 ash samples from coals of different sources and biomasses. The proposed equations could predict the ash fusion temperatures with correlation of 0.8 to 0.92 and standard deviation of 45 to 80°C [8].

In the current work, as mentioned before, we used a database that was spread out over a large area, and we used 1040 data sets. Good correlations of of 0.943, 0.957, and 0.97 in regression and 0.985, 0.99, and 0.995 in ANFIS modeling were achieved for the predictions of IDT, ST, and FT, respectively; this is not in agreement with the Seggiani (1999) suggestion [8] that the statistical approach can give good results in a small number of samples with similar composition. This work shows that a suitable combination of variables can improve the accuracy of the predictions of AFTs, although it is used in widespread databases.

#### 5. Conclusions

(i)Single-variable regression studies show that the SiO_{2}, Al_{2}O_{3}, and Fe_{2}O_{3}, as well as base, acid, base/acid, and silica ratios, are significant variables for predicting AFTs in the examined database.(ii)The basic oxide is the most significant variable with negative effect on AFTs; consequently, the coal ash is favorable with the least basic oxide.(iii)The acidic oxide is the second significant variable with positive effect on AFTs. Thus, coal ash with a high percentage of acidic oxide is desirable. (iv)The “11 oxide + (Base/Acid) + Silicate ratio” combination is the best predictor of AFTs in nonlinear equations. This prediction is more precise for FT than for ST or IDT.(v)The nonlinear equation with “11 oxide + (Base/Acid) + Silicate ratio” input variable can predict IDT, ST, and FT with acceptable correlation coefficients () of 0.891, 0.917, and 0.94 and standard deviation of about 56–93°F (13–34°C).(vi)The ANFIS prediction improved the correlation coefficients () to 0.972, 0.98, and 0.99 for the prediction of IDT, ST, and FT, respectively. Correlation coefficients of these magnitudes have not been reported in previously published works.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

The authors gratefully acknowledge the financial support provided by the Tehran Science and Research Branch of Islamic Azad University.

#### References

- B. G. Miller,
*Coal Energy Systems*, Elsevier, 2005. - American Society for Testing and Materials Standards,
*Standard Test Method for Fusibility of Coal and Coke Ash*, D1857-03, ASTM, 2003. - B. Caylor, “Coal’s role in the national energy plan,”
*Energia*, vol. 13, no. 3, article 5, 2001. View at Google Scholar - E. C. Winegartner and B. T. Rhodes, “An empirical study of the relation of chemical properties to ash fusion temperatures,”
*Journal of Engineering for Gas Turbines and Power*, vol. 97, no. 3, pp. 395–406, 1975. View at Google Scholar · View at Scopus - F. E. Huggins, D. A. Kosmack, and G. P. Huffman, “Correlation between ash-fusion temperatures and ternary equilibrium phase diagrams,”
*Fuel*, vol. 60, no. 7, pp. 577–584, 1981. View at Google Scholar · View at Scopus - R. W. Bryers and T. E. Taylor, “An examination of the relationship between ash chemistry and ash fusion temperatures in various coal size and gravity fractions using polynomial regression analysis,”
*ASME Paper*75-WA/CD-3, 1975. View at Google Scholar - V. R. Gray, “Prediction of ash fusion temperature from ash composition for some New Zealand coals,”
*Fuel*, vol. 66, no. 9, pp. 1230–1239, 1987. View at Google Scholar · View at Scopus - M. Seggiani, “Empirical correlations of the ash fusion temperatures and temperature of critical viscosity for coal and biomass ashes,”
*Fuel*, vol. 78, no. 9, pp. 1121–1125, 1999. View at Publisher · View at Google Scholar · View at Scopus - S. A. Lolja, H. Haxhi, R. Dhimitri, S. Drushku, and A. Malja, “Correlation between ash fusion temperatures and chemical composition in Albanian coal ashes,”
*Fuel*, vol. 81, no. 17, pp. 2257–2261, 2002. View at Publisher · View at Google Scholar · View at Scopus - E. Jak, “Prediction of coal ash fusion temperatures with the F*A*C*T thermodynamic computer package,”
*Fuel*, vol. 81, no. 13, pp. 1655–1668, 2002. View at Publisher · View at Google Scholar · View at Scopus - G. Özbayoğlu and M. E. Özbayoğlu, “A new approach for the prediction of ash fusion temperatures: a case study using Turkish lignites,”
*Fuel*, vol. 85, no. 4, pp. 545–552, 2006. View at Publisher · View at Google Scholar · View at Scopus - Y. P. Liu, M. G. Wu, and J. X. Qian, “Predicting coal ash fusion temperature based on its chemical composition using ACO-BP neural network,”
*Thermochimica Acta*, vol. 454, no. 1, pp. 64–68, 2007. View at Publisher · View at Google Scholar · View at Scopus - SPSS, “Version 13,” SPSS Inc., Help Files, 2004.
- R. Babuska,
*Fuzzy Modeling for Control*, Kluwer Academic, Boston, Mass, USA, 1998. - J. Jantzen, “Neurofuzzy modelling,” Tech. Rep. 98-H-874, Department of Automation, Technical University of Denmark, 1998. View at Google Scholar
- 2009, http://en.wikipedia.org/wiki/Neuro-fuzzy.
- J.-S. R. Jang and C.-T. Sun, “Neuro-fuzzy modeling and control,”
*Proceedings of the IEEE*, vol. 83, no. 3, pp. 378–406, 1995. View at Publisher · View at Google Scholar · View at Scopus