3.1 The influence of the O2/CO2/N2 atmosphere on the combustion characteristic parameters
Figure 2 gives the calculation method of combustion characteristic parameters, including ignition temperature (Ti), maximum heat release rate (vp) and the temperature corresponding to maximum heat release rate (Tp). Ti corresponds to the temperature which coal samples begin to burn, and its value reflects the propensity of coal to spontaneous combustion. (Tp, vp) is the point that the heat release rate of coal reaches to the maximum. It shows as the valley point on the DSC curve. The corresponding results are shown in Table 3. When the heating rate was constant, the values of Ti and Tp in the O2/CO2/N2 atmosphere visibly increased compared with that in the O2/N2 atmosphere. This indicated that a delay of ignition and heat release existed during coal oxygen-lean combustion in the O2/CO2/N2 atmosphere. This result was consistent with the literature [34–35, 39–40].
In the O2/CO2/N2 atmosphere, the values of vp for BLT coal samples decreased by 2.55 mW/mg, 4.18 mW/mg, and 16.51 mW/mg at 10°C/min, 15°C/min, and 20°C/min, respectively, compared with that in the O2/N2 atmosphere, because the decreasing O2 concentration leaded to a reduction of local reactivity [28]. The values of vp for YW coal samples increased by -3.52 mW/mg, 1.43 mW/mg, and 4.25 mW/mg at 10°C/min, 15°C/min, and 20°C/min, respectively, because the lower heat conduction coefficient of CO2 than N2 caused a thermal lag effect [41]. It can be seen that the influence of the O2/CO2/N2 atmosphere on the maximum heat release rate was restricted by the coal rank. The low-rank coal burned faster due to its low carbon content, and O2 had a significant impact on the maximum heat release rate. The high-rank coal contained more carbon and burned slowly, and CO2 had a significant impact on the maximum heat release rate.
Table 3
Characteristic temperatures of BLT and YW coal
Experiment number
|
Ti (°C)
|
Tp (°C)
|
Vp (mW/mg)
|
BLT 1
|
401.16
|
491.94
|
32.75
|
BLT 2
|
394.45
|
522.18
|
41.01
|
BLT 3
|
370.52
|
539.18
|
45.36
|
BLT 4
|
431.16
|
521.05
|
30.20
|
BLT 5
|
425.75
|
549.20
|
36.83
|
BLT 6
|
385.53
|
573.96
|
28.85
|
YW 1
|
513.40
|
568.39
|
24.84
|
YW 2
|
519.18
|
608.91
|
26.68
|
YW 3
|
529.82
|
641.79
|
29.74
|
YW 4
|
530.42
|
598.98
|
21.32
|
YW 5
|
546.19
|
636.54
|
28.11
|
YW 6
|
544.22
|
654.27
|
35.99
|
3.2 The influence of the O2/CO2/N2 atmosphere on the kinetic parameters by KAS method
Figure 3 shows the changes in the values of apparent activation energy and correlation coefficients (R2) by KAS method, in the two atmospheres. For BLT coal, as the conversion rate increased, the values of apparent activation energy all first decreased, then increased, and finally decreased. In the range of 0.5 ~ 0.15 conversion rate, it was in the low-temperature oxidation process. Since the R2 was lower than 0.80, the values of apparent activation energy were not accurate and cannot be compared. When the conversion rate was higher than 0.15, the coal sample was ignited, and the values of apparent activation energy values in the two atmospheres appeared a sudden increase. During the combustion process, the values of apparent activation energy in the O2/CO2/N2 atmosphere were approximately 33%~58% lower than that in the O2/N2 atmosphere. This was because the heat released by coal combustion accumulated more easily in the O2/CO2/N2 atmosphere than in the O2/N2 atmosphere, as a result of the reduction of 6 vol.% O2 and the addition of 5 vol.% CO2 (low heat conduction coefficient).
For YW coal, as the conversion rate increased, the values of apparent activation energy in the two atmospheres kept decreasing. When the conversion rate was 0.05, it was in the low-temperature oxidation process, and the values of apparent activation energy in the O2/CO2/N2 atmosphere were approximately 40% higher than that the O2/N2 atmosphere. This has been confirmed in the research of others [25, 35]. When the conversion rate was higher than 0.05, the coal was ignited, and the values of apparent activation energy in the two atmospheres were close. The influence of the atmosphere was no longer obvious.
In addition, the R2 in the two atmospheres was greater than 0.99. However, for BLT coal, a decrease behavior was showed in the conversion rates ranges of 0.10 ~ 0.25 in O2/N2 and 0.15 ~ 0.45 in O2/CO2/N2, respectively. The reason was that the precipitation of the remaining volatiles was promoted by the heat release of separated volatiles combustion [28, 42–43], and the precipitation and combustion of volatile was significantly deferred in the O2/CO2/N2 atmosphere compared with that in the O2/N2 atmosphere, as a result of the slightly lower diffusivity of volatiles in CO2 than in N2 [28, 44–45] and the lower mass flux of oxygen to the volatiles flame [28, 46].
3.3 The influence of the O2/CO2/N2 atmosphere on the heat release
The reaction rate between oxygen and coal is the key factor influencing the heat release rate [47]. Studying the relationship between the heat release rate and reaction rate is beneficial to understand in the heat release process during coal oxygen-lean combustion in the O2/CO2/N2 atmosphere, which can provide a theoretical foundation for revealing the law of coalfield fire spreading. Since the value of the reaction rate constant can directly reflect the reaction rate, the reaction rate constant was used instead of the reaction rate in this study. According to our previous research [35], the kinetic mechanism functions of BLT and YW coal were Jander (Diffusional (3-D)) and three-level chemical reaction, respectively. The values of pre-exponential factor were calculated though Eq. (5), and then the values of reaction rate constant were obtained though Eq. (3).
Figure 4 and 5 show the DSC-k(T) curves of BLT and YW coal, respectively. The conversion rate corresponding to the maximum heat release rate was taken as a segment point, and the DSC-k(T) curves were divided into two stages: the increasing stage and the decreasing stage of the heat release rate. The conversion rate corresponding to the maximum heat release rate of BLT and YW coal was always about 0.80 and 0.50, respectively., indicating that the conversion rate corresponding to the maximum heat release rate was only related to the coal rank, and not corrected to the atmosphere.
At the increasing stage of the heat release rate, the heat release rate of two coals increased exponentially with the increasing reaction rate constant. At the decreasing stage of the heat release rate, the heat release rate of YW coal decreased exponentially with the increasing reaction rate constant, whereas the heat release rate of BLT coal decreased linearly with the increasing reaction rate constant. This was because the ash on the coal surface obviously hindered the continued diffusion of O2 into the coal pores, resulting in a decrease in the subsequent heat release rate. BLT coal contained less carbon and less ash, so its heat release rate decreased more slowly than that of YW coal and decreased linearly with the increasing reaction rate constant. ExpGro1 exponential model (see Eq. (6)) was selected to fit the DSC-k(T) curves at the increasing stage of heat release rate for the two coal samples. The model showed a high degree of fit, with the R2 for both BLT and YW coal sample above 0.94. Therefore, the relationship between the heat release rate and reaction rate constant for both BLT and YW coal sample can be effectively expressed by the model. The relationship between the heat release rate and the reaction rate constant is approximately as Eq. (7).
$$y={y}_{0}+{A}_{1}\text{e}\text{x}\text{p}\frac{x}{{t}_{1}}$$
6
$$\left|\text{D}\text{S}\text{C}\right|={y}_{0}+{A}_{1}\text{e}\text{x}\text{p}\frac{\text{d}\alpha /\text{d}t}{{t}_{1}f\left(\alpha \right)}$$
7
Where, y0 is the offset. A1 is the amplitude, t1 is the width.
In order to quantitatively analyze the relationship between y0, A1 and t1 and heating rate, Figs. 6 and 7 show the changes in y0, A1 and t1 with the heating rate, respectively. There was a linear relationship between y0, A1, t1 and heating rate for YW coal. For BLT coal, y0, A1 and t1 were basically linear with the heating rate in the O2/N2 atmosphere, whereas there was a non-linear relationship between y0, A1, t1 and heating rate in the O2/CO2/N2 atmosphere. Furthermore, y0-β, y0-A1 and y0-t1 curves were fitted respectively. The reaction rate constant was calculated using Eq. (1). On the whole, there was a following relationship between the heat release rate and reaction rate of coal, as follows
$$\left|\text{D}\text{S}\text{C}\right|={a}_{1}\beta +{b}_{1}+({a}_{2}\beta +{b}_{2})\times \text{e}\text{x}\text{p}\frac{\text{d}\alpha /\text{d}t}{\left({a}_{3}\beta +{b}_{3}\right)f\left(\alpha \right)}$$
8
Where, a1 and b1 are constants related to y0, a2 and b2 are constants related to A1, a3 and b3 are constants related to t1, as seen in Table 4.
Table 4
The value of constants in equation. (8)
Coal samples | Atmosphere | Conversion rate range | a1 | b1 | a2 | b2 | a3 | b3 |
---|
BLT | O2/N2 | 0.05–0.80 | 1.59 | 20.78 | -1.67 | -13.35 | -0.09 | 0.19 |
0.80–0.95 | 1.28 | 21.04 | -0.04 | 0.83 | / | / |
YW | O2/N2 | 0.05–0.50 | 0.33 | 22.12 | 0.08 | -25.54 | -0.00 | -0.01 |
0.50–0.95 | 0.75 | 0.84 | -0.35 | 21.32 | -0.42 | 2.89 |
O2/CO2/N2 | 0.05–0.50 | 1.47 | 6.42 | -1.23 | -12.41 | -0.00 | 0.00 |
0.50–0.95 | 1.31 | -4.48 | 0.05 | 13.32 | -0.42 | 3.07 |