Gas generation from coal: taking Jurassic coal in the Minhe Basin as an example

The gas generation features of coals at different maturities were studied by the anhydrous pyrolysis of Jurassic coal from the Minhe Basin in sealed gold tubes at 50 MPa. The gas component yields (C1, C2, C3, i-C4, n-C4, i-C5, n-C5, and CO2); the δ13C of C1, C2, C3, and CO2; and the mass of the liquid hydrocarbons (C6+) were measured. On the basis of these data, the stage changes of δ13C1, δ13C2, δ13C3, and δ13CO2 were calculated. The diagrams of δ13C1–δ13C2 vs ln (C1/C2) and δ13C2–δ13C1 vs δ13C3–δ13C2 were used to evaluate the gas generation features of the coal maturity stages. At the high maturity evolution stage (T > 527.6 °C at 2 °C/h), the stage change of δ13C1 and the CH4 yield are much higher than that of CO2, suggesting that high maturity coal could still generate methane. When T < 455 °C, CO2 is generated by breaking bonds between carbons and heteroatoms. The reaction between different sources of coke and water may be the reason for the complicated stage change in δ13CCO2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta^{{{13}}} {\text{C}}_{{{\text{CO}}_{{2}} }}$$\end{document} when the temperature was higher than 455 °C. With increasing pyrolysis temperature, δ13C1–δ13C2 vs ln (C1/C2) has four evolution stages corresponding to the early stage of breaking bonds between carbon and hetero atoms, the later stage of breaking bonds between carbon and hetero atoms, the cracking of C6+ and coal demethylation, and the cracking of C2–5. The δ13C2–δ13C1 vs δ13C3–δ13C2 has three evolution stages corresponding to the breaking bonds between carbon and hetero atoms, demethylation and cracking of C6+, and cracking of C2–5.


Introduction
The natural gas generated from different origins provides its formation characteristics during its evolution (Behar et al. 1992(Behar et al. , 2008(Behar et al. , 2010Wang et al. 2013). Coal-formed gas has been an important field of natural gas and plays an important role in China's natural gas resources (Dai et al. 2014). The heterogeneity of coal is strong, and the evolution features of coal-formed gas are more complex than those of other source rocks (Sun et al. 2013). The kinetics of hydrocarbon generation extend the hydrocarbon generation of coal under laboratory conditions to geologic history (Butala et al. 2000;Ping'an et al. 2009;Shuai et al. 2006), which simplifies the gas evolution of the coal. d 13 C is one of the most important properties of natural gas. The d 13 C of coal-formed methane has obvious evolution stage features (Cramer 2004;Cramer et al. 1998;Liu and Xu 1999), and the d 13 C of natural gas could be used to trace the gas origin (Schoell 1980;Song and Xu 2005). The process of gas formation and evolution could be reflected by the d 13 C relationship of different carbon numbers (Chung et al. 1988;Peng et al. 2009). The diagram of d 13 C 2 -d 13 C 1 versus d 13 C 3 -d 13 C 2 could reflect the maturity of natural gas (Jenden et al. 1993), and the diagrams of ln (C 2 /C 3 ) vs ln (C 1 /C 2 ) and d 13 C i -d 13 C j vs ln (C i /C j ) can be used to distinguish the gas origin (Prinzhofer and Huc 1995); these indicators are largely associated with thermal evolution (Tian 2006).
In recent years, there has been increasing attention on the CO 2 generated during the thermal evolution of coal (Fu et al. 2019a;Killops et al. 1996;Kotarba and Lewan 2004;Lewan and Kotarba 2014;Shuai et al. 2013a, b). CO 2 could be generated from organic matter at any maturity stage (Shuai et al. 2013b); at low maturity, CO 2 could be generated from the cracking of kerogen itself, and the disproportionation between organic matter and water could generate CO 2 at the high maturity stage (Seewald 2003;Seewald et al. 1998). The d 13 C of CO 2 in natural gas could be used for gas source contrast (Song and Xu 2005), and some post-reformation of natural gas can be reflected by d 13 CO 2 (Rice 1993;Zhang et al. 2004).
The evolution of coal-formed gas abundance and the d 13 C of gas can be measured by pyrolysis experiments (Behar et al. 1995;Hill et al. 2007;Lewan et al. 1985). Coal has a favorable gas generation ability and low oil generation potential (Hill et al. 2003;Shuai et al. 2006), and because of its molecular sieve nature and chemical adsorption property, liquid hydrocarbons can be generated when R o % = 0.7, but oil expulsion cannot occur (Cook 1988;Johnston et al. 1991). Therefore, a closed system simulation experiment could reflect the gas generation process from coal, and anhydrous pyrolysis in sealed gold tubes under pressure could simulate the geological condition very well (Behar et al. 2008;Hill et al. 2007;Tang et al. 1996).
Previous hydrocarbon-generation kinetics research on Jurassic coal from the Minhe Basin shows that the mean values of the activation energies of CH 4 and C 1-5 are 64.55 kcal/mol and 63.93 kcal/mol, respectively, with a frequency factor of 1.0 9 10 14 s -1 (Fu et al. 2019b). However, the stage evolution characteristics of gas generation have not been studied. In this study, the anhydrous pyrolysis of Jurassic coal from the Minhe Basin in sealed gold tubes was simulated. The gas component yields (C 1 , C 2 , C 3 , i-C 4 , n-C 4 , i-C 5 , n-C 5 , and CO 2 ); the d 13 C of C 1 , C 2 , C 3 , and CO 2 ; and the mass of liquid hydrocarbons (C 6? ) were measured. Then, the stage d 13 C values of C 1 , C 2 , C 3 , and CO 2 were calculated, and the diagrams of d 13 C 1 -d 13 C 2 versus ln (C 1 /C 2 ) and d 13 C 2 -d 13 C 1 versus d 13 C 3 -d 13 C 2 were obtained to study the gas evolution of coal.

Sample and closed system pyrolysis
The immature Jurassic coal used in this study was from the Minhe Basin, Gansu Province, China (Fig. 1), and the basic geochemical characterization of the sample is listed in Table 1.
The coal was pulverized via a 100 mesh, and no other chemical treatment was conducted. The pyrolysis experiments were conducted in sealed gold tubes (40 mm 9 5 mm). The sample (15-60 mg) was loaded in a gold tube with one end welded closed, and the air was replaced by argon. The open end was then welded closed by TIG (Pulse arc welding) with the closed end protected in room temperature water. The samples in sealed gold tubes were heated in one furnace in 12 separated stainless steel autoclaves. To avoid the influence of phase differentiation, the experiment was conducted at a constant pressure of 50 MPa, in which the fluid phase is basically in a single phase (Fu et al. 2019a, b). Two series of experiments were conducted, and the gold tubes were heated for 10 h from room temperature to 250°C and then heated to 600°C with heating rates of 20°C/h and 2°C/h. Additional descriptions of the experimental procedure can be found in (Pan et al. 2007;Shuai et al. 2013b;Wang et al. 2013).

Gas analysis
After the pyrolysis experiment, the gold tube was punctured in a closed system and connected to an Agilent 7890A gas chromatograph (GC) that has a Wasson ECE module for analyzing the gas composition (C 1-5 , CO 2 ) using the external standard method to calculate the quantities of each gas component. The GC employed an HP-AL/ S capillary column (25 m 9 0.32 mm 9 0.8 lm) and used helium as the carrier gas. The column temperature was programmed from 60°C (held for 3 min) to 190°C (held for 3 min) at 25°C/min.

Stable carbon isotope analysis
The d 13 C levels of C 1 , C 2 , C 3 , and CO 2 analysis were conducted on an Isoprime 100 mass spectrometer that interfaced an Agilent 6890 GC. The GC was equipped with a CP-Poraplot Q column (27.5 m 9 0.32 mm 9 10 lm), and helium was used as the carrier gas. The column temperature was programmed from 50°C (held for 3 min) to 190°C (held for 5 min) at 25°C/min. The d 13 C of each sample was measured twice, and the average deviation was less than 0.3%.

C 61 analysis
The C 6? composition was separated to C 6-14 and C 14? for analysis. After the gas composition was analyzed, the C 6-14 liquids were collected by a liquid nitrogen cold trap (4 mL, quartz bottle) (Behar et al. 1995) and injected with 3 mL of n-hexane. Then, the gold tube was cut into pieces and put into the bottle to ensure that the C 6? liquids could be dissolved completely. Deuterated-24-alkanes were used as the internal standard to determine the quantities of the liquids. The analysis was conducted on an Agilent 7890A GC that employed a DM-5 capillary column (30 mm 9 0.32 mm 9 0.25 lm) and used helium as the carrier gas. The initial oven temperature was 40°C and held for 5 min, after which the oven was heated to 290°C with a heating rate of 4°C/min and held at 290°C for 15 min. The C 14?   3 Results

Liquid (C 61 ) yield
The yields of C 6? should be the sum of C 6-14 and C 14? . Table 2 and Fig. 2a show the cumulative yields of the liquid hydrocarbons (C 6? ) in the experiments. As Fig. 2a shows, the yields of C 6? increased rapidly at low temperature (T \ 410°C), and the maximum yield was 78.17 mg/ g TOC (20°C/h, 406.9°C). The yield began to decrease rapidly when the temperature [ 410°C because of cracking.

Gas yields
The yields of C 1 , C 2 , C 3 , C 4 , C 5 , and CO 2 generated from coal pyrolysis are shown in Fig. 2 and Table 2. As a whole, C 2-5 ( Fig. 2b) peaked at 431°C (2°C/h) and 478°C (20°C/h), which is close to the end of the C 6? liquids cracking. The volumetric productivities of C 2 , C 3 , C 4 and C 5 were in descending order, and the relationship with temperature was similar. However, the cracking temperature was different and decreased with increasing carbon number; in addition, the cracking temperature was different at different heating rates, which manifested as the temperature became lower during the slow heating rate than during the rapid heating rate. The yields of methane were low before 410°C and increased rapidly after 410°C, which was close to the cracking temperature of the C 6? liquids. The maximum yield of methane was 257.35 mL/g TOC (599.3°C, 2°C/ h), which was much higher than that of all the other gases. At the same temperature, the gas yield during the slow heating rate was higher than during the rapid heating rate, which was the result of the complementary relationship between time and temperature on the chemical kinetics (Connan 1974).
CO 2 is the most important nonhydrocarbon gas of the coal thermal maturity process (Shuai et al. 2013b;Tang et al. 1996). Figure 2h shows that the yield of CO 2 was significant at 336°C and much higher than that of the hydrocarbon gases. The yield of methane exceeded that of CO 2 when the temperature was higher than 455°C (20°C/ h) and 431°C (2°C/h), and the maximum CO 2 was Table 2 Yields of gases and their d 13 C value and liquid hydrocarbon generated by pyrolysis experiment  Fig. 2 The C 6? yields of (a), C 2-5 (b), C 5 (c), C 4 (d), C 3 (e), C 2 (f), C 1 (g), and CO 2 (h) during the pyrolysis experiments Gas generation from coal: taking Jurassic coal in the Minhe Basin as an example 615 112.13 mL/g TOC (599.3°C, 2°C/h). The XRD results (Table 3) indicated that the sample contains very small amounts of carbonate minerals, and the detection of liquid HCl showed the same results. The carbonate minerals contained in the sample could not generate such a large amount of CO 2 , which means that the CO 2 was mainly from an organic origin.

The d 13 C of gases
The d 13 C values of C 1 , C 2 , C 3 , and CO 2 are shown in Fig. 3 and Table 2. The d 13 C of C 1-3 decreased with increasing pyrolysis temperature and then began to enrich 13 C, manifesting as d 13 C 1 \ d 13 C 2 \ d 13 C 3 at the same temperature and heating rate. Taking 2°C/h as an example, the d 13 C 1 decreased over the range of 334.8-383.5°C and began to enrich the 13 C after 383.5°C and the lightest d 13 C 1 was -37.71%. Many other pyrolysis experiments of kerogen and crude oil samples have similar changes (Tang et al. 2000;Tian et al. 2007;Wang et al. 2013;Xiong et al. 2004). The d 13 C of CO 2 ranged from -24.03% to 20.62%, and this is much less than that of hydrocarbon gases and is different from the report of Shuai et al. (2013b), showing that CO 2 became enriched in 13 C with increasing temperature but increased to -20.74% at 358.7°C and then fluctuated between -20.63% and -21.03% to 479.0°C; in addition, with increasing temperature, it became lighter after 479.0°C. The instantaneous changes in d 13 C could not be obtained because of the closed pyrolysis system. To solve this problem, the following formula was used to calculate the stage value of d 13 C according to the law of isotopic conservation and based on the gas yields and d 13 C (Shuai et al. 2013b): where, d 13 C Ti ' is the stage value of d 13 C from temperature T i-1 to T i ; d 13 C Ti and d 13 C Ti-1 are the cumulative d 13 C values of gases at temperatures T i and T i-1 , respectively; and V Ti and V Ti-1 are the cumulative volumetric yields of gases at temperature T i and T i-1 , respectively.
The stage d 13 C values of C 1 , C 2 , C 3 , and CO 2 are shown in Fig. 4. Ethane and propane are generated and cracked during pyrolysis, so d 13 C Ti ' , calculated by Eq. (1), can be divided into two processes. There is a certain temperature Fig. 4 The different stage d 13 C values of C 1 , C 2 , C 3 , and CO 2 during the pyrolysis experiment Gas generation from coal: taking Jurassic coal in the Minhe Basin as an example 617 overlap between ethane and propane at the end of the formation and at the beginning of the cracking (Shuai et al. 2006;Tian 2006), so the calculated values that were close to the temperature when cracking began were discarded. As shown in Fig. 4, with increasing pyrolysis temperature, the stage d 13 C values of C 2 and C 3 first decreased and then increased during the generation stage. During the cracking processes, the light carbon ethane and propane were cracked prior because the bond energy of 12 C-13 C is higher than that of 12 C-12 C (Arneth and Matzigkeit 1986;Stevenson et al. 1948;Tang et al. 2000). Therefore, the cumulative values of d 13 C 2 and d 13 C 3 in Fig. 3 began to increase when the cracking temperature was reached. The stage d 13 C value of C 1 is significantly different from its cumulative value. Cramer (2004) explained the evolution of the d 13 C value of coal-generated methane according to kinetics and divided methane generation into the following four reaction stages: (1) cleavage of heteroatoms, (2) demethylation reactions, (3) second cracking of long-chain alkanes and cyclic compounds generated, and (4) polycondensation reactions. Taking 2°C/h as an example, with the increase of pyrolysis temperature, the stage d 13 C 1 value decreased to -38.85% at 383.5°C, which means it was close to the end of the reaction stage (1) (Cramer 2004;Liao et al. 2007). After that, the stage value of d 13 C 1 began to increase, and the peak value of -17.59% was reached at 527.6°C. This stage corresponded to the second cracking of liquid hydrocarbon (C 6? ) and wet gas (C 2-5 ) one after another (Fig. 2a, b), which reflects that reactions stage (2) and (3) occurred (Cramer 2004). Then, the stage value of d 13 C 1 started to decrease with increasing temperature, and the yield continued to increase which corresponding to the reaction stage (4). However, the changes in d 13 C 2 and d 13 C 3 (Figs. 3, 4) and their remaining amounts (Fig. 2) showed that the cracked gases were almost complete and that d 13 C became heavier, which means that the methane generated at this stage came from not only the cracking of C 2-5 but also other reactions that may be the main source of methane. The reaction between water and coke in the higher evolution stage could generate some of the CH 4 , but as Eq. (2) below shows The gas generation yields of CO 2 and CH 4 of this evolution stage should be equal; however, as Table 2 shows, the increase in CH 4 is much higher than that in CO 2 . Therefore, there must be another origin of CH 4 . Cramer (2004) points out that the polycondensation reaction is the reason that d 13 C 1 decreases during this stage and that coal has gas generation potential at this stage.
The stage d 13 C value of CO 2 was complex and changed between -37.04% and -10.67%, and it is generally believed that the d 13 CO 2 of an organic origin is less than -10% (Dai et al. 1996;Wycherley et al. 1999); in addition, the calculated stage d 13 C value of CO 2 shows that the CO 2 in the experiment was of an organic origin and that the XRD result (Table 3) also supported this view. Shuai et al. (2013b) suggest that the d 13 C of organic origin CO 2 could be as high as 18%. However, in this paper, the data do not support this conclusion, which may be due to the complex composition of coal. The stage d 13 C value of CO 2 changed smoothly before 455°C (remaining at approximately -20%), which was significantly higher than the stage d 13 C value of methane at the same stage. This is because during this stage, the CO 2 is mainly sourced from the cleavage of the heteroatom reaction of coal (Seewald et al. 1998), and the carbon that is attached to the heteroatom is more enriched in 13 C than the aliphatic side chains (Cheng et al. 2009). At the high maturity stage, CO 2 is generated from the reaction between H 2 O and organic matter (Lewan 1997;Seewald et al. 1998). Because of anhydrous pyrolysis, H 2 O may originate from the pyrolysis of coal or from the bond water in clay minerals (Wang et al. 2013). When the temperature was higher than 455°C, the stage d 13 C value of CO 2 changed substantially, which may be related to the different causes of coke (cracking of aromatization carbon or liner carbon) reacting with H 2 O. In fact, as the experimental temperature increases, the number of free radicals in the reaction system increases accordingly. The generated CO 2 should result of the reaction of oxygen free radicals and carbon free radicals. The oxygen free radicals and carbon free radicals here are from the cracking of H 2 O and coal, respectively. At the same time, we believe that the production of C radicals is mainly affected by temperature and pressure, and the isotope fractionation process is more complicated, which directly leads to the complex characteristics of the d 13 C value of CO 2 .

d 13 C i -d 13 C j vs ln C i /C j
The d 13 C of hydrocarbon gases has a good correlation with thermal maturity and source rock (Dai and Qi 1989;Schoell 1983), and the gas composition is very sensitive to thermal maturity (Wang et al. 2013). Therefore, the d 13 C id 13 C j vs ln C i /C j diagram can be used to indicate the thermal maturity and source rock of natural gas (Hill et al. 2003;Prinzhofer and Huc 1995;Tian et al. 2010). The d 13 C 2 -d 13 C 3 vs ln C 2 /C 3 diagram based on the pyrolysis experimental data (Fig. 5) shows a significant feature of kerogen pyrolysis gas (Prinzhofer and Huc 1995), which is characteristic of coal.
The d 13 C 1 -d 13 C 2 versus ln C 1 /C 2 diagram could be used to reflect the maturity of natural gas, the leakage of natural gas and the mixing of thermogenic gas and biogenic gas (Jenden et al. 1993;Prinzhofer and Huc 1995). The closed system pyrolysis experiment could exclude the influences of leakage and biogenic gas, which should only be correlated with the thermal maturity and the sample. The d 13 C 1d 13 C 2 versus ln C 1 /C 2 diagram based on the pyrolysis experimental data (Fig. 6) shows obvious stage changes. Taking 2°C/h as an example, the first stage was 334.8-358.7°C, and C 1 /C 2 decreased with increasing temperature, which means that the increasing ratio of ethane was higher than that of methane (Wang et al. 2013).
The d 13 C 1 -d 13 C 2 decreased slightly, which was possibly because both methane and ethane were generated from cleavage of the heteroatom side chain, and the isotope fractionation effect of methane was higher than that of ethane. The second stage was from 358.7 to 407.2°C, corresponding to the end of the cleavage heteroatom reaction. The C 1 /C 2 almost remained unchanged, and the d 13 C 1 -d 13 C 2 decreased rapidly because of high hereditability from the source rock of ethane (James 1983;Xie et al. 1999). The next stage was from 407.2 to 455.3°C, both the methane and ethane yields were increased because of the cracking of C 6? (Fig. 2) during this interval, and the d 13 C 1 and d 13 C 2 were also increased. C 1 /C 2 increased, which means that the increase ratio of methane was higher than that of ethane, and the d 13 C 1 -d 13 C 2 increased slightly, showing that the degree of increasing d 13 C 1 heaviness was larger than that of d 13 C 2 . When the temperature was higher than 455.3°C, ethane began to crack (Fig. 2), and C 1 /C 2 increased rapidly with decreasing d 13 C 1 -d 13 C 2 .
4.3 d 13 C 3 -d 13 C 2 vs d 13 C 2 -d 13 C 1 The d 13 C 3 -d 13 C 2 and d 13 C 2 -d 13 C 1 values decrease with increasing thermal maturity (Jenden et al. 1993;Prinzhofer and Huc 1995), so the d 13 C 3 -d 13 C 2 vs d 13 C 2 -d 13 C 1 diagram can be used to reflect the maturity of natural gas (Jenden et al. 1993). However, many simulation experiments show the opposite change (Guo et al. 2009(Guo et al. , 2011Tian 2006). The indexes above largely depend on the maturity of natural gas (Tian 2006). As Fig. 7 shows, in this experiment, the d 13 C 3 -d 13 C 2 vs d 13 C 2 -d 13 C 1 diagram shows obvious stage changes. Taking 2°C/h as an example, when the temperature was lower than 383.5°C, the d 13 C 3 -d 13 C 2 increased slightly from 2.16% to 2.94%, and Fig. 6 d 13 C 1 -d 13 C 2 vs ln C 1 /ln C 2 ((1) the increase ratio of ethane was higher than that of methane, and the isotope fractionation effect of methane was higher than that of ethane; (2) the end of the cleavage heteroatom reaction, where C 1 /C 2 remained almost unchanged and d 13 C 1d 13 C 2 decreased rapidly; (3) both the methane and ethane yields increased because of C 6? undergoing cracking; (4) the ethane begins to crack) the d 13 C 2 -d 13 C 1 increased rapidly from 5.36% to 9.82%. This is because during this stage, d 13 C 1 decreased, but d 13 C 2 first decreased and then increased, and d 13 C 3 remained almost unchanged (Figs. 2, 3). Then, when the temperature was between 383.5 and 431.2°C, the d 13 C 3d 13 C 2 first decreased and then increased, changing over the range of 2% to 3%, and the d 13 C 2 -d 13 C 1 remained almost unchanged. This result indicates that the variations in d 13 C 1 , d 13 C 2 and d 13 C 3 were basically the same, and the main sources of these gases were the same during this stage. In contrast to Fig. 2, we conclude that the cracking of C 6? was the main source of these gases. The last stage was from 431.2 to 455.3°C, where d 13 C 2 -d 13 C 1 remained at approximately 10% and d 13 C 3 -d 13 C 2 increased rapidly from 1.95% to 5.08% because of propane undergoing cracking before ethane.

Conclusions
Anhydrous pyrolysis experiments were conducted in sealed gold tubes with Jurassic coal from the Minhe Basin. The hydrocarbon composition and CO 2 yields and the d 13 C of gas compositions were obtained. The following conclusions were drawn: The stage d 13 C value of methane shows a change trend of the ''S'' type that first decreased, then increased and finally decreased. The stage d 13 C values of ethane and propane first decreased and then increased during the generating stage and became enriched in 13 C during the cracking stage because the bond energy of 12 C-13 C was higher than that of 12 C-12 C. The stage d 13 C value of CH 4 decreased when T [ 520°C, but the d 13 C values of cracked ethane and propane increased, which proved that coal still had the potential to generate CH 4 at high maturity. The stage d 13 C value of CO 2 changes when T [ 455°C since the difference caused coke to be reacted with H 2 O.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright Fig. 7 d 13 C 3 -d 13 C 2 vs d 13 C 2 -d 13 C 1 ((1) d 13 C 1 decreased, but d 13 C 2 first decreased and then increased, and d 13 C 3 remained almost unchanged; (2) the variations in d 13 C 1 , d 13 C 2 and d 13 C 3 were basically the same, and the main sources of these gases were generated from C 6? cracking; (3) propane cracked before ethane) holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.