Abstract

In this study, a nonisothermal kinetics analysis of petcoke was performed at heating rates of 10, 15, and 20°C/min using thermal gravimetric analysis (TGA). The behaviour of petcoke at different gasification stages (dewatering, volatilization, char burning, and burnout) was studied. The effect of heating rate on the activation energy of petcoke gasification was also investigated. The activation energy of petcoke was estimated using different kinetic models that include volume reaction model (VRM), shrinking core model (SCM), random pore model (RPM), Coats and Redfern model (CRM), and normal distribution function (NDF). The NDF model was modified in this study. It was found that the experimental data were best fitted with the modified normal distribution function (MNDF) and SCM. The results also showed that activation energy decreases as heating rate increases, leading to reduction in gasification completion time.

1. Introduction

Petcoke is a heavy crude oil refining coproduct. It is identified as a black-colored, carbon-rich solid. Despite the few human health or environmental risks posed by the exploitation of petcoke, it has a high economic value and many industrial applications. It is mostly used as a boiling and combusting fuel in industrial, power generation, and cement plants. Moreover, it can be used as a raw material in manufacturing processes. In addition, petcoke is a promising substitute for steam coal in power plants because of its higher heating value, carbon content, and low ash, compared to bituminous coals [1, 2]. However, petcoke gasification is a difficult process because of its high content of fixed carbon [3, 4] and low volatile matter [5]. Moreover, the mass transfer of petcoke is influenced by the porosity, pore size, and volume as well as diffusivity and tortuosity of the carbon substrate [6, 7]. To tackle this issue, different solutions have been studied such as modifying particle size, inclusion of catalysts, and cogasification [2, 8]. In cogasification, as an example, petcoke is combined with additional fuels such as coal or biomass to improve its low reactivity [9, 10]. Appropriate catalysts can also be incorporated into the gasification of petcoke to improve their low reactivity [4, 11, 12]. The other alternative approach to achieve the high carbon conversion is to increase the residence time of the fuel particles [7].

The above discussion demonstrates the importance of petcoke gasification development. Therefore, understanding the gasification kinetics of fuels is key to finding solutions to potential problems of gasification as well as improving operating conditions to develop efficient gasification processes. In addition, a detailed understanding of reaction kinetics is vital for the feasibility, design, and scaling of gasification applications. It will also provide valuable information for proper design and operation of gasifiers. Therefore, the combustion characteristics, chemical elements, and technical analysis of the fuels fed into the gasification reactor must be understood. For this purpose, thermal gravimetric analysis (TGA) and derivative thermogravimetry (DTG) curves provide information on fuel reactivity properties such as ignition, peak, and burnout temperatures [13, 14]. The change in the mass of the sample caused by devolatilization during thermal decomposition is monitored by TGA as a function of temperature or time. The maximum reaction rate of the sample can then be obtained through the first derivative of the TGA curve (dx/dt). This curve is known as DTG.

In 1975, Tyler and Smith [15] studied the reactivity of petcoke with CO2 over a temperature range of 1018–1178 K. A rate order of approximately 0.6 and activation energies in the range of 203–237 kJ/mol were obtained for the petcoke with respect to carbon dioxide concentration. A constant rate order and activation energy of about 45% with a burn-off range of 21%–45% were reported. Revankar et al. [16] investigated the effect of particle size, porosity, and thickness on the steam gasification of petcoke, where they reported increase in the rate constant with decrease in both grain size and pellet size, with or without a catalyst. In contrast, activation energy was found to be independent of particle size for the noncatalytic gasification, but the frequency factor decreased with particle size. However, this trend is in reverse for catalytic gasification.

Zou et al. [17] investigated the kinetic characteristics of petcoke gasification with CO2 at 1248–1323 K and 0.1 MPa. They proposed a normal distribution function model (NDF) to fit the kinetics data. The gasification rate was observed to increase with the conversion rate (X) up to the peak of X = 0.3. The activation energy and reaction order range of 198 kJ/mol and 0.54–0.88 were reported, respectively, for the CO2 gasification of petcoke. Yoon et al. [18] conducted a TGA study of coal and petroleum coke for cogasification at 15°C/min heating rate. The petcoke activation energies of 53.73 and 46.03 kJ/mol were reported using kinetic models of the shrinking core model (SCM) and the integrated model or the modified volume reaction model (MVRM), respectively. In a similar study, Nemanova et al. [19] explored the cogasification of petcoke and biomass at a rate of 10°C/min using TGA. They analyzed TGA data using the volume reaction model (VRM) and determined an activation energy value of 121.5 kJ/mol for the petcoke steam gasification. Recently, Jayaraman and Gokalp [8] studied the effect of particle size on steam gasification of petcoke using TGA and mass spectrometry analysis. Relatively high reaction rates and 90% conversion efficiency were observed for smaller-sized particles (30 μm). The study proposed an operation temperature of ≥950°C for petcoke gasification under steam or blended steam ambiences with efficient fuel conversion.

A summary of the studies discussed above is given in Table 1. Without doubt, the review shows there is no discrete activated state as confirmed by the wide variations of calculated activation energies [20]. However, the existing models have some limitations. For example, the assumptions of the SCM model may not accurately match the real conditions, even though the model is considered the most appropriate simple representation for the majority of reacting gas-solid systems [21]. According to SCM, the reaction occurs first at the ash layer and the outer skin of the particle and then moves into the nonreacted core. Alternatively, the reaction may occur along a diffuse front, which is a kind of intermediate behaviour between the SCM and CRM models. In addition, considering the fast reactions, the heat release rate is sufficiently high to cause significant temperature gradients within the core of the particle. However, there is a lack of kinetic models that accurately analyze the reactions of more complex heat distributions in gas-solid systems. Therefore, this study investigates the gasification characteristics and kinetic analysis of petcoke using TGA and proposed modified normal distribution function (MNDF) at differential heating rates. Based on the kinetic data collected by TGA, some of the common kinetic models were utilized to estimate the activation energy of petcoke. The results were then comparatively analyzed and discussed.

2. Materials and Methods

2.1. Sample Preparation

For this study, petcoke obtained from PETRONAS Melaka Refinery with a Hardgrove Grindability Index (HGI) of 111 and a particle size range of 5–20 mm was used. The petcoke was crushed and milled to a size of 0.7 mm. It was then dried for 24 hours at 110°C using an industrial oven. Afterward, the combustion characteristics of petcoke were investigated using PerkinElmer STA 6000. The initial weights of the petcoke samples were varied since the TGA results were unaffected by the sample weight. To investigate the effect of heating rate on the gasification phases, TGA experiments were carried out at heating rates of 10, 15, and 20°C/min. In addition, an untreated petcoke sample (as received) was used to investigate the influence of sample preparation conditions on TGA curves. The descriptions of samples used in this study are outlined in Table 2.

2.2. Ultimate and Proximate Analysis

Ultimate analysis was carried out to determine the percentage of carbon, hydrogen, nitrogen, and sulphur chloride of the petcoke using a CHNOS elemental analyzer under consideration of ASTM D5373. Proximate analysis was carried out using ASTM D5142-90 to acquire preliminary results on the moisture content, volatile matter, fixed carbon, and ash content of petcoke samples. In addition, a bomb calorimeter with standard reference of the ASTM D5865 test method was used to determine the calorific value of the samples. Dried petcoke sample was used for both ultimate and proximate analysis. However, the calorific value test was carried out on both dried and untreated petcoke samples.

Table 3 shows the results for ultimate, proximate, and calorific value analyses. Different heating rates of 10, 15, and 20°C/min were used, while the TGA experiment temperature was kept at 800°C. The nitrogen and oxygen flow rates were kept at 100 ml/min for all experiments. The change in sample weight was recorded every second. A continuous supply of nitrogen was maintained for 40 minutes to eliminate moisture until the reaction temperature specified for the generation of char was attained. Afterward, the effect of heating rate on the ignition temperature and ignition time of petcoke was analyzed by burning each petcoke sample under different heating rates in the presence of oxygen. Finally, the combustion characteristics of the char were measured for each sample.

2.3. Kinetic Analysis

Similar to other organic or chemical matters like biomass or coal, the kinetic of petcoke decomposition can also be predicted by the following equation:where dX/dt, K, f(X), and X are the isothermal reaction rate, the reaction rate constant, the reaction model equation, and the extent of reaction (or the conversion rate), respectively.

The reaction rate K is presented by the Arrhenius expression as follows:

By substituting equation (2) into equation (1), the rate law can be expressed as

The extent of reaction, X, can be obtained using the following expression:where m0, mt, and mf are the initial sample mass, sample mass at time t, and final sample mass, respectively.

Some of the important rate models, f(X), used to describe the kinetic behaviour of solid-state reactions were listed in ref. [22]. Among them, the VRM [23], SCM [21], random pore model (RPM) [24], MVRM [25, 26], NDF [17], and Coats and Redfern (CR) model [27, 28] have been used by researchers to predict the kinetic parameters of petcoke. Based on the VRM (or homogeneous model) assumption, the particle density changes uniformly due to the distribution of gas inside the solid particles. According to the VRM, the reaction kinetic expression can be described aswhere (1 − X) is the remaining fraction of volatile material in the sample.

The SCM was first established by Yagi and Kunii [29]. SCM assumes that the reaction first occurs at the outer skin of the spherical particle, and it then moves into the nonreacted core which shrinks in size gradually during reaction process. According to SCM, the reaction kinetic model is as follows:

The RPM takes into account the pore structure and its evolution in the course of reaction. Based on the RPM, the gasification rate can be written aswhere ψ is the particle structure parameter. The RP model with the value of ψ equal to zero gives the VR model.

The volume reaction model was improved to MVRM (or integrated model) by Kasaoka et al. [25] and Yang et al. [26]. It was extended for all reaction orders by adding a new parameter (n) as follows:

Zou et al. [17] proposed the normal distribution function (NDF) as a model to predict kinetic parameters. The conventional distribution function is approximated by a Gaussian distribution that yields a mean value and standard deviation. Based on NDF, the reaction kinetic model can be written aswhere Xm and ω are the maximal gasification rate and the width of the curve, respectively.

On the contrary, the reaction rate can be presented as a function of temperature for nonisothermal solid-state reactions:where β is the heating rate and dX/dT is the nonisothermal reaction rate.

Considering equations (3) and (10), the reaction rate can be represented as

By taking natural logarithms of each side,

Integration of equation (11) yields

Equation (13) is the integral expression of the rate law. If we consider (Ea/RT) as u, equation (13) becomes

The p(u) in equation (14) known as the temperature integral cannot be integrated by one of the methods of calculus. However, this is not a serious limitation because it can be approximated via empirical interpolation formulas like an integral method based on the Coats and Redfern (CR) equation [27, 28]. Therefore, p(u) in equation (14) can be estimated using a Taylor series expansion to yield the following expression [22]:

To begin, we studied the weight reduction of petcoke for different heating rates. Afterward, some of the abovementioned models were used to fit the experimental data obtained by TGA. A modified form of NDF was then proposed. Finally, using an Arrhenius plot, the accuracy of the kinetic models in petcoke activation energy approximation was comparatively evaluated.

3. Results and Discussion

3.1. TGA and DTG Analysis

Petcoke was evaluated using TGA analysis to understand its ignition temperature and combustion characteristics. The weight reduction for different heating rates of 10, 15, and 20°C/min is demonstrated in Figure 1. Since there was no evidence of weight loss within the first 30 minutes of the TGA, Figures 1(a) and 1(b) only depict the values that exceed 240°C and 30 min, respectively. It is discernible from the TGA that the dewatering stage starts around 243°C (at 32.13 min), 463°C (at 39.40 min), 419°C (at 30.00 min), and 465°C (at 39.51 min) and continues till the temperature reached approximately 518°C (at 57.71 min), 523°C (at 42.08 min), 525°C (at 34.25 min), and 519°C (at 41.98 min) for samples 1 to 4, respectively. Taking into cognizance the dewatering stage, the weight of all samples reduced by nearly 6%. However, the differential weight values can be attributed to the different applied heating rates.

It is also important to note that the time required for complete petcoke dewatering for treated and untreated samples (samples 2 and 4) is almost the same. Therefore, since the effect of sample drying is negligible, more energy can be saved by using untreated samples. Volatilization is initiated immediately after the dewatering stage for all samples and is completed at about 542°C (at 58.38 min), 548°C (at 42.65 min), 554°C (at 34.85 min), and 542°C (at 42.45 min) for samples 1 to 4, respectively. Subsequent to the volatilization stage is the char combustion phase. The combustion process was finally completed at different temperatures of 572°C (at 61.16 min), 595°C (at 45.78 min), 620°C (at 38.21 min), and 583°C (at 44.83 min) for samples 1 to 4, respectively. The gray area in Figure 1(b) shows our implementations for petcoke combustion, which can be used as a tool to predict the petcoke combustion behaviour under certain heating rates. This specifically implies the rate of petcoke weight loss during the combustion reaction at heating rates between 10 and 20°C/min can be predicted from the gray area.

To study the effect of reaction rate on different reaction stages, the start and end points of time and temperature for each phase were calculated and the results are presented in Table 4. Irrespective of the heating rates employed, the time required for total volatilization was the same. A similar trend can be observed for the combustion stage. However, there is a slight disparity between the temperature values of the volatilization and combustion stages, which can be attributed to the different heating rates applied. This indicates that the effect of heating rate on the volatilization and combustion stages is negligible. Significant deviations are only observed for the time required for the dewatering stage to be completed. Therefore, the most pertinent deduction is that the time required for moisture removal from petcoke can be significantly reduced by increasing the heating rate. Moreover, comparison of the heating rates shows 15°C/min is the most effective for all stages of petcoke gasification.

The TGA and DTG curves for the thermal decomposition of petcoke at a heating rate of 15°C/min are shown in Figure 2(a). As earlier explained, the TGA contrasts the weight loss of the petcoke with temperature change, whereas the DTG represents the time derivative of the weight loss for the same temperature change as TGA. It can be inferred from Figure 2(a) that the rate of weight loss is maximum at 558°C (peak temperature). Cursory evaluation of the TGA curve indicates the peak temperature occurred in the char combustion phase. Thus, the rapid weight reduction corresponds to the combustion phase.

It is also important to mention that the change in the values of derivative weight, as observed from the DTG curve, indicates that the progressive thermal breakdown of the organic matter present in the petcoke occurs within the range of 480–620°C. In addition, the burnout temperature of 620°C represents the temperature where sample oxidation is complete. The derivative weight curves for the thermal decomposition of petcoke at different heating rates are shown in Figure 2(b). It is observed that the maximum weight losses for heating rates between 10 and 15°C/min occur at almost the same temperature of 558°C. Moreover, the maximum weight loss for the treated sample (sample 1) and the untreated sample (sample 2) occurred at 558 and 556°C, respectively. However, at the higher heating rate of 20°C/min, the maximum weight loss was recorded at 568°C.

Conversely, the DTG profiles clearly show that the rate of weight loss slightly decreases with the increase in heating rate. The maximum conversion rates of −32.15, −29.87, and −28.52% min−1 were observed for the heating rates of 10, 15, and 20°C/min, respectively. Remarkably, the DTG profile of the untreated sample (sample 4) shows a higher rate of weight loss (−37.47% min−1) than the others. It is evident that the reaction sensitivity of petcoke to heating rate and sample treatment is very low.

3.2. Arrhenius Plot and Activation Energy

The reaction rates, K, were calculated using selected kinetic models that include, VRM, SCM, RPM, NDF, and CRM. The Arrhenius plots for the gasification reaction of petcoke (Figure 3) were derived using the logarithmic reaction rate ln (K) min−1 versus the reciprocal gasification temperature 1000/T K−1. It is discernible from the Arrhenius plots that the logarithmic gasification reaction rate exhibits a linear relationship with temperature for almost all kinetic models and reaction rates, with the exception of 10°C/min. As observed in Figure 3, the slope of ln (K) versus 1000/T becomes steeper as the heating rate declines. Given the proportionality of the activation energy (Ea) with the slope of the Arrhenius plot, it can be inferred that the activation energy increases as heating rate decreases.

The activation energies Ea for different kinetic models were determined using equation (16), and results are presented in Table 5 for comparison.where slope is the slope of ln (K) versus 1/T from the Arrhenius plot.

Based on the data outlined in Table 5, there is a significant difference between the activation energy values obtained by CRM and NDF and those acquired by the other kinetic models. This difference can be attributed to systematic errors associated with the calculation of Ea using the integral isoconversional kinetic models. Considering each kinetic model, it can be noted that activation energies obtained from different heating rates were not the same and they all increase with decreased heating rates. This is because more energy is required to initiate release of the refractory molecules enclosed in petcoke in the form of volatiles when using lower heating rates. The comparative analysis of Ea values obtained in this study with those reported in earlier studies (Table 1) shows that more reliable Ea values can be estimated using SCM, RPM, and VRM, while those obtained using NDF and CRM are overestimated. The results indicated that increasing the heating rate leads to reduction in activation energy. As a result, less energy is needed to activate atoms or molecules to a state where they can undergo a chemical reaction. For example, using SCM, the activation energy of petcoke is reduced by 71.53% with the increase in heating rate from 10 to 20°C/min.

3.3. Modified Normal Distribution Function (MNDF)

As stated above, the NDF overestimated the values of activation energy. Nonetheless, NDF has some features which make it a valuable approximation function for complex distributions of variables. The NDF model was then modified by calculating the right values of maximal gasification rate (Xm) and the curve width (ω) of the normal distributed function. Considering the normal distribution function (equation (9)), its exponential function increases to zero while the value of ω decreases. Meanwhile, the value of conversion (X) changes from 0 to 1. Consequently, the graph of NDF for Xm equal to mean is symmetric with respect to the line X equal to mean (X = mean). Considering these features of NDF, choosing the right values of ω and Xm is necessary to best fit the graph of experimental data with the straight line. To this end, the values of 1 and 0.736 are proposed in this study for Xm and ω, respectively. Using these values, the following modified NDF (MNDF) is proposed:

Plots of −ln (k) versus 1000/T, derived from different heating rates and MNDF and NDF kinetic models, are illustrated in Figure 4.

It is worthwhile to mention that using NDF, the logarithmic reaction rate results in a straight line with positive slope. Therefore, the negative logarithmic function (−ln) was used in this study to make the slope of the Arrhenius graph negative. This change in sign does not affect the results of activation energy. However, this finding was not reported in the literature [17]. It can be clearly seen in Figure 4 that the MNDF yields much tighter curves that establish a better fit with experimental data than the conventional NDF model. Furthermore, compared with the NDF model, the MNDF provided a better minimization of the discrepancy in calculated activation energies. The activation energies are consistent with those extracted from the literature.

4. Conclusions

In this study, TGA was performed on petcoke to study its combustion characteristics under different heating rates. The different stages of combustion, namely, dewatering, volatilization, char burning, and burnout for petcoke were examined. Some common kinetic models (VRM, SCM, RPM, CRM, and NDF) were utilized to estimate the activation energy of the combustion of petcoke. The NDF model was modified to best fit the experimental data. The kinetic analysis showed that the combustion results of petcoke under the heating rate of 10°C/min show a longer transition stage between dewatering and volatilization. It was also found that the activation energy of petcoke is significantly affected by increase in heating rate. The TGA results revealed that the total time needed for gasification of petcoke was shortened with the use of 20°C/min heating rate, which indicated that increasing the heating rate leads to reduction in activation energy. Thus, shorter gasification finishing time is expected while using higher heating rate. The proposed modified kinetic model (MNDF) was found to accurately predict the gasification kinetics the most and achieved 50% reduction in activation energy. In addition, the results show that the duration for complete petcoke combustion can vary from 38 to 61 minutes. The proposed kinetic model provides valuable information for proper design and operation of the gasifier reactors. It also improves the understanding and development of the gasification process.

Data Availability

The TGA data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

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

Acknowledgments

This research was funded by University Technology PETRONAS (UTP), YUTP grant number 015LC0-078.