Geochemical Characteristics and Origin of Shale Gases From Sichuan Basin, China

Natural gases from the Taiyang (shallow), Jiaoshiba (middle), and Weirong (deep) shale gas fields in the southern Sichuan Basin were analyzed for molecular and stable carbon isotopic compositions to investigate the geochemical characteristics and gas origins. All the gases belong to shale gas from the Upper Ordovician–Lower Silurian shale and are dominated by methane with gas wetness generally less than 0.83%. The δ13C1 values are −28.5‰, −30.3‰, and −35.2‰ in Taiyang, Jiaoshiba, and Weirong shale gas fields, respectively. The extremely high thermal maturity is the controlling factor for the enrichment of 13C in methane, with a minor contribution from the heavy carbon isotope of the organic matter in the Ordovician Wufeng Formation. Fischer–Tropsch-type synthesis of hydrocarbon gas from CO2 and H2 contributes to the increase of wet gas, which results in the offset from the δ13C1∼wetness linear trend in the Taiyang and Jiaoshiba gas fields. Methane, ethane, and propane in the Taiyang shale gas field have increasing δ13C values with increasing burial depth, which is mainly caused by diffusive migration. All gases are characterized by a complete carbon isotopic reversal trend (δ13C1 > δ13C2 > δ13C3), and it is mainly caused by the reversible free-radical reactions with the conversion from alkane to alkyl groups, with some contribution from the Fischer–Tropsch-type synthesis. The results of this study will improve our understanding of the geochemical characteristics of shale gases from different burial depths and have important implications for future shale gas exploration in the deep and shallow layers.


INTRODUCTION
Carbon isotopic compositions of methane and its homologs are important parameters to investigate the origin and generation history of natural gas (Tang et al., 2000;Ni et al., 2011;Dai et al., 2014a). Several studies have shown that the carbon isotope is dependent on the maturity of source rocks, the type and carbon isotope of source organic matter, and postgenetic processes (Stahl and Carey, 1975;Schoell, 1980;Dai, 1992;Tang et al., 2000;Dai et al., 2009;Ni et al., 2011;Curiale and Curtis, 2016;Ni et al., 2019a). Due to the kinetic isotopic fractionation effects (Lorant et al., 1998;Tang et al., 2000), primary gases are usually characterized by the normal carbon isotopic trend, that is, δ 13 C 1 < δ 13 C 2 < δ 13 C 3 . However, recently, a number of studies have documented a depletion of 13 C in ethane and propane and an enrichment of 13 C in methane at high maturity, resulting in the carbon isotopic reversal between methane and ethane (δ 13 C 1 > δ 13 C 2 ) (Ferworn et al., 2008;Burruss and Laughrey, 2010;Rodriguez and Philp, 2010;Zumberge et al., 2012;Hao and Zou, 2013;Tilley and Muehlenbachs, 2013;Dai et al., 2014aDai et al., , 2016bDai et al., , 2017Zhang et al., 2018b;Liu et al., 2018;Ni et al., 2018;Xia and Gao, 2018;Feng et al., 2020). Zumberge et al. (2012) demonstrated a reversed carbon isotopic maturity trend in ethane and propane in gases from Barnett and Fayetteville shales, with thermal maturity (Ro%) over 1.5%. Dai et al. (2014b) explored this reversed trend to even higher thermal maturity and found the carbon isotopic maturity trend becomes positive again at stages with gas wetness (C 2+ /C 1+ , %) less than 1.4%. Ni et al. (2018) show such isotopic depletion of 13 C in thermogenic ethane and propane. One common feature of these carbon isotopic anomalies is that they occur at relatively high temperatures. To date, many studies have been carried out to investigate the mechanisms for the carbon isotopic characteristics of shale gases (Ferworn et al., 2008;Burruss and Laughrey, 2010;Rodriguez and Philp, 2010;Zumberge et al., 2012;Hao and Zou, 2013;Xia et al., 2013;Gao et al., 2014;Zhang et al., 2018b;Xia and Gao, 2018;Feng et al., 2020;Mi et al., 2022). Moreover, isotopic reversal may correlate with the gas generation potential, and the reversal magnitude may correspond to the amount of gas retained in the shale (Tilley et al., 2011;Feng et al., 2018;Wang et al., 2021;Zhao et al., 2021).
As the largest shale gas producer outside North America, after 10-years effort, China has achieved effective exploitation of marine shale gas resources in the Sichuan Basin and its surrounding areas with a middle burial depth of 2,000-3,500 m. In the Sichuan Basin, ultra-deep shale gas represents shale gas with a burial depth >4,500 m, deep shale gas means shale gas with a burial depth of 3,500-4,500 m, and shallow shale gas means shale gas with a burial depth <2,000 m. With the success of the burial depth of 2,000-3,500 m, recently, great efforts have been made in the shallow and deep layers. In the future, with the improvement of the technologies, shale gas exploration will extend to areas shallower than 2,000 m and deeper than 3,500 m. How about the geochemistry of shale gas from shallow and deep layers? Are there any differences among the shale gases from the shallow, middle, and deep layers? Here, we describe the molecular and stable carbon isotope geochemistry of natural gases from the Wufeng-Longmaxi shale in the Taiyang (vertical burial depth, 500-2,250 m), Jiaoshiba (vertical burial depth, 2,250-3,500 m) and Weirong (vertical burial depth, 3,550-3,880 m) shale gas fields in the Sichuan Basin. Two main objectives of this study will be addressed: 1) Investigating the origin of the shale gases from the three gas fields with different burial depths based on the molecular and stable carbon isotopic compositions and 2) demonstrating the geochemical differences among shale gases from these three gas fields and exploring the potential causes. The results presented in this study will improve our understanding of the geochemical differences among the shale gases with deep, middle, and shallow burial depths and interpretation of the carbon isotopic composition of shale gases in the future.

GEOLOGICAL SETTINGS
The Sichuan Basin, abundant in natural gas resources, is a large superimposed basin and covers an area of about 180 × 10 3 km 2 ( Figure 1). The sedimentary strata are well-developed with thicknesses up to 12,000 m. It developed nine sets of gas source rocks and has found 25 conventional and tight gasproducing zones (18 being marine facies) and one marine facies shale gas-producing zone. The nine sets of gas source rocks include the Sinian Doushantuo and Dengying formations, Cambrian Qiongzhusi Formation, Upper Ordovician Wufeng-Lower Silurian Longmaxi formations, Permian Qixia, Maokou, Dalong, Longtan/Wujiaping formations, and Triassic Xujiahe Formation (Figure 2). The former seven sets belong to marine sapropelic source rocks. The Sinian Doushantuo and Dengying formations were deposited during a sea transgression. It is distributed mainly in the basin margins in the northwest, northeast, and southeast of the Sichuan Basin, with very thin thickness in the basin. The Cambrian Qiongzhusi Formation is characterized by black carbonaceous shale with an average thickness of 180 m in the basin. The organic matter is dominated by Type I and Type II 1 kerogen with total organic content (TOC) up to 6% in the southern and northern Sichuan, and the thermal maturity equivalent Ro% is 2%-5% (Xu et al., 2011). The Permian Qixia and Maokou formations are a set of carbonate deposits and are widely distributed in the basin with thickness between 150-400 m and TOC generally less than 2% . The Permian Dalong Formation is a set of black mudstone and siliceous mudstone. It is mainly distributed in the northern Sichuan Basin with a thickness up to 30 m (Xia et al., 2013). It has a TOC between 0.43%-21.2% (average = 8.3%, n = 38) and δ 13 C values of kerogen from −24‰ to −29‰ . The Longtan Formation/Wujiaping Formation belongs to marine-continental transitional source rocks, where the Longtan Formation is a set of coal measure gas source rocks, similar to the Xujiahe Formation. Wujiaping Formation is contemporaneous with Longtan Formation, but mainly a set of marine sapropelic source rocks. The Triassic Xujiahe Formation is dominated by Type II 2 and Type III organic matter, which is a set of coal measure gas source rocks. The Wufeng-Longmaxi formations are the only industrial exploration and development strata of shale gas. The Wufeng Formation has a thickness of 2-11 m, and the Longmaxi hot shale is concentrated on the bottom 10-40 m of the Lower Longmaxi Formation. The Wufeng-Longmaxi hot shale has a burial depth of 500-5,500 m and is deposited in an environment of deep water shelf facies and stably distributed. By 2020, Jiaoshiba, Changning, Weiyuan, Huangjinba, Taiyang, Weirong, Yongchuan, and other shale gas fields have been discovered ( Figure 1). In China, the cumulative proved geological reserves of shale gas were 2.0 × 10 12 m 3 , the completed production capacity was 300 × 10 8 m 3 /year, and the annual shale gas output was 200.55 × 10 8 m 3 in 2020.

Taiyang Shale Gas Field
The Taiyang shale gas field is located in the south of Xuyong County, Sichuan Province, and is a part of the Zhaotong National Shale Gas Demonstration Zone (Figure 1). It belongs to the Taiyang anticline, which was a narrow and steep anticlinal structure cut and damaged by strike-slip faults. It has been uplifted continuously since the Late Jurassic, and the uplifted range is about 2,000 m, which resulted in the denudation of the strata above Permian and possible loss of free gas. The black organic-rich siliceous shale and carbonaceous shale of the Wufeng-Longmaxi formations are the main targets of the shale gas field development, with a vertical burial depth (excluding the horizontal portion) of 500-2,250 m and formation pressure coefficient of 1.2-1.47, which belongs to shallow, atmospheric-overpressure shale gas reservoirs. The high-quality reservoirs are 13-18 m thick and stable in distribution, with TOC content of 0.42%-9.05% with an average of 2.92%. The organic matter is dominated by Type II 1 , and the equivalent Ro% ranges from 3.0% to 3.9% (average = 3.45%) Liang et al., 2021;Zou et al., 2021). The shale porosity is 2.92%-8.01%, with an average of 4.98%. The average brittleness is 72%, the gas content is 3.3-5.51 m 3 /t, and the average content of clay minerals is 24.7%. By the end of 2020, the accumulative proved reserves had been 1,365 × 10 8 m 3 . The completed production capacity was 8 × 10 8 m 3 /year, and the annual shale gas output was 4.0 × 10 8 m 3 in 2020.

Jiaoshiba Shale Gas Field
The Jiaoshiba shale gas field is located in the east of Chongqing and is a part of the Fuling National Shale Gas Demonstration Zone ( Figure 1). It has an area of about 575.92 km 2 . The main body the of Jiaoshiba gas field is a box anticlinal structure spreading NE-SW, with faults controlling the east and west sides and large-scale faults developing in the southeast. The uplift was late and large, about 3,000 m, and the limestone of the Triassic Jialingjiang Formation is exposed, but the gas reservoir is well-preserved. The Wufeng-Longmaxi organicrich shale is the main target layer of the shale gas field. The optimal reservoir section has a burial depth (excluding the horizontal portion) of 2,250-3,500 m. It is stably distributed with a thickness of about 38-40 m and a pressure coefficient of 1.55. It is dominated by Type I organic matter with an average TOC of 3.76%, porosity of 4.87%, brittleness of 62.4%, gas content of 2.96 m 3 /t, clay mineral content of 31.6%, and equivalent Ro% of 3.1% (2.6%-3.6%) (Dai et al., 2014a;Wang et al., 2019;Zou et al., 2021). By the end of 2020, the accumulative proved reserves had been 6008.14 × 10 8 m 3 . The completed production capacity was 100 × 10 8 m 3 /year, and the annual shale gas production was 67 × 10 8 m 3 in 2020.

Weirong Shale Gas Field
The Weirong shale gas field is located in Neijiang and Zigong counties of Sichuan Province and is a part of the Changning-Weiyuan National Shale Gas Demonstration Zone ( Figure 1). The main body of the gas field belongs to the Baimazhen syncline of the low-fold structural belt in southwestern Sichuan. The structure is gentle, and the fault is not developed. Influenced by the uplift evolution of the Sichuan Basin in the Late Cretaceous, the Jurassic was exposed on the surface, and the burial depth (excluding the horizontal portion) of the gas layer was 3,550-3,880 m. The overall preservation condition of the gas layer is good, and the formation pressure is high, that is, 1.94-2.06. The highquality reservoirs are 25-39 m thick and stable in distribution, with an average TOC of 2.80%, porosity of 6.07%, brittleness of 64%, gas content of 3.15 m 3 /t, clay mineral content of 34%, and equivalent Ro% of 2.25% (1.8%-2.7%) Zou et al., 2021). By the end of 2020, the cumulative proved reserves had been 1,247 × 10 8 m 3 . The completed production capacity was 10 × 10 8 m 3 / year, and the annual shale gas output was 5.4 × 10 8 m 3 in 2020.

SAMPLES AND METHODS
Eleven shale gas samples from the Weirong shale gas field in the Changning-Weiyuan National Shale Gas Demonstration Zone (marked as Weirong), 10 shale gas samples from the Taiyang shale gas field in the Zhaotong National Shale Gas Demonstration Zone (marked as Taiyang), and 10 shale gas samples from the Jiaoshiba shale gas field in the Fuling National Shale Gas Demonstration Zone (marked as Jiaoshiba) were collected ( Figure 1).
Samples were pure gases and collected directly from the wellheads or separators in the shale fields. Double-ended stainless steel bottles (10 cm diameter about 1,000 cm 3 volume) were used to collect the samples, which were equipped with shut-off valves with a maximum pressure of 22.5 MPa. The bottles and the lines were flushed with the sample gases for 15-20 min to remove air contamination. The pressure inside the bottle was kept above atmospheric pressure. Leakage tests by immersing bottles in a water bath were performed for all samples.
The molecular and isotopic compositions of the gas samples were determined at the PetroChina Research Institute of Petroleum Exploration and Development (RIPED). The molecular composition was determined using an Agilent 7890A gas chromatograph, which was equipped with a flame ionization detector and a thermal conductivity detector. Gas components were separated using six columns (2 Foot 12% UCW982 on PAW 80/100 mesh, 15 Foot 25% DC200 on PAW 80/100 mesh, 10 Foot Hayesep 80/100 mesh, 10 Foot Molecular Sieve 13X 45/60 mesh, 8 Foot Molecular Sieve 5A 60/80 mesh, and 3 Foot Hayesep Q 80/100 mesh). GC oven temperature was set at 85°C for 40 min. Immediately prior to the molecular composition analysis, the GC is calibrated by certified gas standards, which is also used to monitor the analytical drift by injecting them before, during, and after the analysis. The analytical precision and accuracy are typically better than 0.01% for methane, ethane, and propane.
Stable carbon isotopes were determined using a Thermo Delta V mass spectrometer, which was interfaced with a Thermo Trace GC Ultra gas chromatograph (GC). Methane, ethane, propane, and CO 2 were separated on a gas chromatograph using a fused silica capillary column (PLOT Q 27.5 m × 0.32 mm × 10 μm), converted into CO 2 in a combustion interface, and then injected into the mass spectrometer. The temperature programming of the GC oven was as follows: heated from 33°C to 80°C at 8°C/min, then to 250°C at 5°C/min, and maintained at 250°C for 10 min. An oilrelated gas from the Tazhong gas field in the Tarim Basin in China is used as our laboratory working standard. More than 800 measurements including both online and offline of the compound-specific carbon isotopic compositions have been carried out by 10 laboratories (Dai et al., 2012b). International measurement standards (NBS19 and L-SVEC CO 2 ) are used to perform the two-point calibrations. The carbon isotopic values were reported in δ-notation per mil (‰) relative to VPDB. All analyses were duplicated, and the analytical precision and accuracy are better than ±0.3‰.

Molecular Composition
Methane is the dominant composition, which accounts for 98.82% (average, n = 10) in the Taiyang gas samples, 98.13% (average, n = 10) in the Jiaoshiba gas samples, and 96.91% (average, n = 11) in the Weirong gas samples ( Table 1). The content of C 2+ heavy hydrocarbon (mainly C 2 H 6 and C 3 H 8 ) is very small, generally less than 1%. It is 0.61% (average, n = 10) in the Taiyang gas samples, 0.58% (average, n = 10) in the Jiaoshiba gas samples, and 0.45% (average, n = 11) in the Weirong gas samples. In general, Taiyang shale gas has the highest content of both methane (98.82%) and C 2+3 heavy hydrocarbon gases (0.59%), while Weirong gas samples have the lowest content of both methane (96.91%) and C 2+3 heavy hydrocarbon gases (0.44%). All the gases have extremely low gas wetness, which has a value of 0.21%-0.83%, indicating they belong to dry gas. The content of nonhydrocarbon gases such as CO 2 and N 2 is also very less ( Table 1). Both the Taiyang and Jiaoshiba gases have a trace amount of CO 2 and N 2 , less than 1%. Weirong gas has similar content of N 2 (<1%) but a relatively higher content of CO 2 , in which CO 2 is generally more than 1% with an average of 2.05% (average, n = 11) and is up to 6.64% in Well Wei 43-1HF.

DISCUSSION
In order to further explore the gas origin, published data about shale gas from the marine facies sapropelic organic matter (Type I and II 1 kerogen) of the Upper Ordovician Wufeng-Lower Silurian Longmaxi shale in the Changning-Weiyuan, Fuling, and Zhaotong national shale gas demonstration zones, Pengshui and Yongchuan area in the Sichuan Basin, shale gas from the marine facies sapropelic organic matter of the Barnett shale in the Fort Worth Basin, Fayetteville shale in the Arkoma Basin, and New Albany shale in the Illinois Basin were also investigated (Cao et al., 2020;Chen et al., 2020;Dai et al., 2016b;Feng et al., 2019;  Note: All the wells are horizontal and comprise vertical and horizontal portions. Well depth represents the depth including the actual vertical portion and the horizontal portion. Burial depth represents the actual vertical portion. FIGURE 3 | Stable carbon isotopic compositions versus 1/n where n is the carbon number of methane, ethane, and propane from Taiyang, Jiaoshiba, and Weirong gas fields in the Sichuan Basin (A). Gases from these three gas fields all have carbon isotopic reversal among the methane, ethane, and propane (δ 13 C 1 > δ 13 C 2 > δ 13 C 3 ). Also shown are published data about shale gases from the Changning-Zhaotong, Fuling, and Weiyuan national shale gas demonstration zones (B) and shale gases from the Barnett, New Albany, and Fayetteville shale (see text). (C) Modified from Chung et al. (1988).

Origin of the Gas
Carbon isotopic reversal occurs widely among the C 1 -C 3 n-alkanes, that is, methane and its homologs get more depleted in 13 C with increasing carbon number (δ 13 C 1 > δ 13 C 2 > δ 13 C 3 ) ( Figure 3A; Table 1). In general, this is consistent with previous studies of the shale gases from the Changning-Zhaotong, Weiyuan, and Fuling national shale gas demonstration zones ( Figure 3B). In contrast, gas samples from the Barnett and New Albany shales mostly have a positive carbon isotopic distribution pattern among the C 1 -C 3 alkane gases (δ 13 C 1 < δ 13 C 2 < δ 13 C 3 ), and those from the Fayetteville shale have similar carbon isotopic reversal ( Figure 3C) (Martini et al., 2008;Rodriguez and Philp, 2010;Stra;poć et al., 2010;Zumberge et al., 2012). Numerical studies have suggested a cut-off δ 13 C 2 value of −28‰ to distinguish the gases sourced from humic and sapropelic organic matter, that is, gases sourced from humic organic matter have δ 13 C 2 > −28‰, while gases sourced from sapropelic organic matter have δ 13 C 2 < −28‰ (Chen et al., 1995;Dai et al., 2009, Dai et al., 2012aDai et al., 2014b;Liang et al., 2002;Ni et al., 2014, Ni et al., 2015, Ni et al., 2019b. Gases from Taiyang, Jiaoshiba, and Weirong gas fields have δ 13 C 2 < −33‰, and should source from sapropelic organic matter according to the criterion of δ 13 C 2 value. However, such a distinguishing criterion is effective for gases with normal carbon isotopic distribution patterns which have not undergone secondary alteration but may not be applicable for gases with carbon isotopic reversal (Dai et al., 2013). Therefore, genetic information should be obtained from indexes and not only the δ 13 C 2 values. "Bernard diagram," based on the carbon isotopic value of methane and the molecular composition of methane and ethane, was often used to determine the gas origin and later modifications (Bernard et al., 1976). This diagram differentiates the thermogenic and bacterial gases and shows the type of kerogen and later secondary alterations such as mixing between thermogenic and bacterial gases, migration, oxidation, and thermal maturation. Except gases from the New Albany in the Illinois Basin that have mixing effects, all other gases fall in the thermogenic gas area and have no impacts from bacterial gas ( Figure 4A). Figure 4B demonstrates that gases from the Weirong gas field are mainly affected by the thermal maturation, but gases from the Jiaoshiba and Taiyang gas fields are likely affected by the migration effects.
In the δ 13 C 1 versus δ 13 C 2 and δ 13 C 2 versus δ 13 C 3 cross-plots, in general, the least mature gases lie in the lower left area and most mature gases lie in the upper right area. Theoretically, gases which have not undergone secondary alteration and are from the same source rocks but with different thermal maturation stages will lie on a straight line. Gases from the Taiyang, Jiaoshiba, and Weirong shale gas fields fall on the bottom right side of the diagram close to the area of sapropelic Type II kerogen and also in the area below the 1:1 line of δ 13 C 1 versus δ 13 C 2 ( Figure 5A). This indicates gases from the Taiyang, Jiaoshiba, and Weirong gas fields were sourced from sapropelic organic matter. Figure 5A also shows a very low carbon isotopic composition of ethane, therefore demonstrating a carbon isotopic reversal between methane and ethane (δ 13 C 1 > δ 13 C 2 ). Similarly, in the plot of δ 13 C 2 versus δ 13 C 3 , gases from the three gas fields fall in the area close to that of the sapropelic Type II kerogen (Berner and Faber, 1996) and demonstrate a carbon isotopic reversal between ethane and propane (δ 13 C 2 > δ 13 C 3 ) ( Figure 5B).

Carbon Isotope of Methane
When plotting the δ 13 C 1 versus gas wetness, gases from Weirong, Weiyuan, Barnett, and Fayetteville have a good linear relationship, while gases from Taiyang and Jiaoshiba are characterized by an enrichment of 13 C in methane and thus shift away from the linear maturity trend ( Figure 7A). The 13 C-enriched methane was considered to be a late-stage methane and result from Rayleigh-type isotopic fractionation of ethane ; however, such a mechanism may not satisfy the mass and isotopic balance requirements (Xia and Gao, 2018). In the southern Sichuan Basin, geological conditions are FIGURE 5 | δ 13 C 1 versus δ 13 C 2 (A) cross-plot and δ 13 C 2 versus δ 13 C 3 cross-plot (B) of gases from Taiyang, Jiaoshiba, and Weirong gas fields in the Sichuan Basin. Lines of Type II kerogen from Berner and Faber (1996) and Huang et al. (2016) are also shown. For the line from Berner and Faber (1996), the carbon isotope of the initial precursor is estimated to be −29‰.
Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 861040 more complex than those in the Appalachian Basin. The Wufeng-Longmaxi shale has experienced very high geotemperature and is at the over-mature stage at present and has suffered extensive uplifts locally. Considering these complex conditions, more influences should be considered in the southern Sichuan Basin. The offset of the δ 13 C 1 values in Taiyang and Jiaoshiba shale gas fields could be caused either by enrichment of 13 C in methane or by an increase in wetness. Compared to the Cambrian and Sinian conventional natural gas in the Weiyuan gas field (−32.7‰ to 32.0‰) (Dai, 2003;Wu et al., 2014), the Silurian Wufeng-Longmaxi shale gas in Taiyang and Jiaoshiba areas is clearly more enriched in 13 C.

Thermal Maturity of the Source Rocks
Various models have demonstrated the linear logarithmic relationship between the carbon isotopic composition of methane and the thermal maturity (Ro%) of source organic matter (Stahl and Carey, 1975;Schoell, 1980;Dai and Qi, 1989;Berner and Faber, 1996). In the oil window, methane is primarily generated from kerogen cracking. When thermal maturity still increases, wet gas cracking begins, which will increase methane but decrease in gas wetness (Prinzhofer and Huc, 1995). Hence, the carbon isotope of methane generally increases with increasing thermal maturity of the source rocks, and methane generated at higher thermal maturity normally has a higher carbon isotope and lower gas wetness. Though the Wufeng-Longmaxi shale in the southern Sichuan Basin is generally all at the over-mature stage, differences still exist in different areas. Shale in the Weiyuan area has a relatively low equivalent Ro% of 2.25% (1.8%-2.7%) Zou et al., 2021); correspondingly, its shale gas has the lowest δ 13 C 1 values (average = −35.2‰) and wetness of 0.47% (0.21%-0.53%), which falls on the linear maturity trend. Shale in the Taiyang area has the highest equivalent Ro% of 3.45% (3.0%-3.9%) Liang et al., 2020a;Liang et al., 2020b;Zou et al., 2021); correspondingly, its shale gas has the highest δ 13 C 1 values (average = −28.5‰) and wetness of 0.61% (0.50%-0.83%), which falls off the linear maturity trend. Compared to the Weirong shale gas, the Taiyang shale gas should have lower gas wetness due to its higher thermal maturity. Similarly, shale in the Jiaoshiba area has an equivalent Ro% of 3.1% (2.6%-3.6%) (Dai et al., 2014a;Wang et al., 2019;Zou et al., 2021), which is higher than that in Weirong but lower than that in Taiyang. Correspondingly, its shale gas also has δ 13 C 1 values (average = −30.3‰) higher those that in Weirong but lower than those in Taiyang, and gas wetness of 0.59% (0.47%-0.67%), which also falls off the linear trend line ( Figure 7A). Compared to the Weirong shale gas, the Jiaoshiba shale gas also should have lower gas wetness due to its higher thermal maturity. As shown in Figure 7A, if thermal maturity is the controlling factor, the Taiyang shale gas should have gas wetness around 0.05%, while the Jiaoshiba shale gas should have gas wetness around 0.07%, which is an order of magnitude lower than the actual determined wetness. According to such low gas wetness, the carbon isotopes of ethane and propane are relatively lower than expected ( Figures 7B, C).
Carbon isotopic differences between methane and ethane [Δ (δ 13 C 2 -δ 13 C 1 )] have been used to reflect the maturity influences on gases. As thermal maturity increases, the isotopic differences between methane and ethane decreases; therefore, the more mature gases normally will have smaller Δ (δ 13 C 2 -δ 13 C 1 ) values (James, 1983;Clayton, 1991;Tang et al., 2000). Shale gases from Taiyang, Jiaoshiba, and Weirong gas fields generally demonstrate a good linear inversed relationship between the δ 13 C 1 value and Δ(δ 13 C 2 -δ 13 C 1 ) value (R 2 = 0.9048). The δ 13 C 1 value increases from −36‰ to −27‰, correspondingly, and the Δ (δ 13 C 2 -δ 13 C 1 ) value changes from −3‰ to −8‰ ( Figure 8A). The carbon isotopic differences between methane and ethane get larger with increasing thermal maturity, which is opposite to the normal trend ( Figure 8A). This is mainly caused by the depletion of 13 C in ethane with increasing thermal maturity. As the thermal maturity increases, methane becomes more enriched at 13 C, while ethane gets more depleted at 13 C; hence, the higher the thermal maturity, the larger the δ 13 C differences between methane and ethane. This linear correlation implies that thermal maturity should be the controlling factor on the carbon isotopes of the shale gas. In addition to the thermal maturity of the source rocks, there should be a process that will increase the amount of wet gas but will not significantly impact the amount of methane ( Figures 8B, C). As shown in Figure 8, methane content generally increases with increasing δ 13 C 1 values in the three gas fields, while shale gas in FIGURE 7 | Plots of wetness versus δ 13 C 1 (A), δ 13 C 2 (B), and δ 13 C 3 (C) of shale gases from Taiyang, Jiaoshiba, and Weirong gas fields. Also shown are shale gases from the published data (see text).
Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 861040 the Taiyang and Jiaoshiba gas fields apparently has higher wet gas.

Fischer-Tropsch-Type Synthesis of Hydrocarbon
Fischer-Tropsch-type synthesis of hydrocarbon from CO 2 and H 2 may have contributed to the increasing amount of wet gas and depletion of 13 C in hydrocarbon gas. The CO 2 content varies apparently in the Taiyang, Jiaoshiba, and Weirong shale gas fields. The CO 2 content is highest in the Weirong shale gas field, 1.19%-6.64% (average = 2.05%); followed by the Fuling shale gas field, 0.21%-0.66% (average = 0.42%); and the lowest in the Taiyang shale gas field, 0.05%-0.15% (average = 0.10%) (Figure 9). Among the 10 gas samples in Taiyang, CO 2 is detected only in five samples. From the Weirong gas field to the Taiyang gas field, the equivalent Ro% increases from 2.25% to 3.45% Liang et al., 2021;Zou et al., 2021), while the CO 2 content decreases from 2.05% to 0.10%. The carbon isotope of CO 2 in the Jiaoshiba area varies from −5.8‰ to 10.4‰ (average = 3.7‰) , which was considered to be the product of thermal metamorphism of carbonate minerals (Dai et al., 2022). Carbonate minerals are abundant in the Wufeng-Longmaxi shale (Dai et al., 2014a;Dai et al., 2016b), and with increasing thermal maturity, thermal metamorphism will generate more CO 2 . There must be a sink for CO 2 with increasing thermal maturity in the Jiaoshiba and Taiyang areas. During the late cantagensis stage, water is likely involved in the generation of methane (Burruss and Laughrey, 2010). A Fischer-Tropsch-type synthesis of hydrocarbon from CO 2 and H 2 is expected to occur through the following two equations: CHx (organic matter) + 2H 2 O → CO 2 + (2 + x/2) H 2 and CO 2 + mH 2 → xCH 4 + yC 2 H 6 + ... + zH 2 O (Tang and Xia, 2011). Among them, CO 2 could be from the thermal metamorphism of carbonate minerals in the shale and also from the water reforming reaction of organic matter. Higher thermal maturity represents longer reaction time and more consumption of CO 2 in the system, so the CO 2 content is the most in the Weirong shale gas field and the least in the Taiyang shale gas field. Such Fischer-Tropsch-type synthesis of hydrocarbon from CO 2 and H 2 has also been reported in the pyrolytic experiments of shale (Gao et al., 2014) and coal samples , and apparent carbon isotopic depletion in gaseous products has been documented.
With continuous cracking of wet gas, the amount of wet gas (mainly C 2 + C 3 ) will decrease, that is, the ethane content is less than 0.5% for Weirong shale gas. Due to the amount of CO 2 and H 2 , the number of hydrocarbon gases generated via Fischer-Tropsch-type synthesis should be minimal. Hence, limited mixing between the late stage shale gas with the Fischer-Tropsch-type synthesized gas will significantly increase the wet gas content but have little impact on the methane content; similarly, the carbon isotopic depletion impact is little on methane, but to some extent significant on ethane and propane. This should be part of the reasons causing the isotopic depletion in ethane and propane.

δ 13 C Value of the Source Organic Matter
In the southern Sichuan Basin, the δ 13 C value of source organic matter may also have influenced on the δ 13 C value of hydrocarbon gases. Generally, hgher δ 13 C values of source rocks correspond to higher δ 13 C values of gases generated FIGURE 8 | Plots of δ 13 C 1 versus Δ(δ 13 C 2 -δ 13 C 1 ), (A), δ 13 C 1 versus CH 4 (B), and δ 13 C 1 versus wetness (C) of gases from the Taiyang, Jiaoshiba, and Weirong gas fields in the Sichuan Basin, China. from them, which is the inheritance effect of the carbon isotope . A prominent positive carbon isotopic excursion of both shale and limestones has been recorded in the Late Ordovician Hirnantian stage (Middle and Upper Wufeng Formation in Sichuan Basin) in North America (Bergström et al., 2006;Orth et al., 1986), Europe (Brenchley et al., 1994;Marshall et al., 1997), and China ( Figure 10) (Chen, 2018;Fan et al., 2009;Wang et al., 1993;Wang et al., 1997), which has been regarded as to be associated with the Gondwana glaciation. The positive carbon isotopic excursion is up to 5‰-7‰ in the Hirnantian limestones (Qing and Veizer, 1994;Marshall et al., 1997;Kump et al., 1999) and is up to −21.1‰ in the Tianjiawan shale in Daozhen county and −26.6‰ in the Nanbazi shale in Tongzi county (Chen, 2018). In addition to the positive carbon isotopic excursion of the organic matter during the Hirnantian stage, several studies have also demonstrated a significant positive carbon isotopic excursion (~8‰) of the organic matter in the early Katian (Lower Wufeng Formation in Sichuan Basin) or the Guttenberg stage in the late Middle Ordovician, which was caused by the decrease in surface-water CO 2 concentrations and the enrichment of microfossil Gloeocapsomorpha Prisca (Jacobson et al., 1988;Pancost et al., 1999;Pancost et al., 2013). Compounds exclusively from Gloeocapsomorpha Prisca may contribute a positive carbon isotopic shift of 3.5‰, indicating the significant contribution from source organic matter (Pancost et al., 2013). The δ 13 C of shale decreased in the Rhuddanian stage (Lower Longmaxi Formation). Actually, the δ 13 C value of kerogen is about −31.0‰ to −27.8‰ in the Changning-Zhaotong areas, −29.9‰ to −28.6‰ in the Jiaoshiba area (Lin et al., 2017), and −36.0‰ to −27.9‰ in the Weiyuan areas (Feng et al., 2018).
Gases sourced from the Wufeng shale are expected to be enriched in 13 C due to the positive carbon isotopic excursion. The continuous sedimentation and thin thickness of Wufeng-Longmaxi hot shale indicate that post-genetic processes such as migration and mixing in the shale layers are possible and contribution from the 13 C-enriched Wufeng shale cannot be ignored. The "sweet spot" of the Wufeng-Longmaxi shale gas is mainly concentrated in the Wufeng Formation and the lower part of the Longmaxi Formation, with a thickness of 15-50 m, and the Rhuddannian stage is the prime stage of the development of the Lower Silurian organic-enriched shale and also the key stage for the formation of "sweet spot" of shale gas (Zou et al., 2016, Zou et al., 2019. According to the graptolite occurrence, sedimentation from the Late Ordovician to the Early Silurian is continuous , Liang et al., 2020aRong et al., 2011, Rong et al., 2019. Particularly, in the Taiyang area, shale with TOC > 1% has a thickness of 50.3-57.0 m, and it is characterized by longitudinal continuous distribution without an interlayer in the middle (Liang et al., 2020b). Diffusive migration has been found in the Jiaoshiba area, where δ 13 C 1 values increased with proximity to the fault zone . In the Taiyang shale gas field, the linear relationship between burial depth and δ 13 C values of methane, ethane, and propane also indicates possible influences from gas migration. However, the contribution from the carbon isotope of source organic matter cannot be the controlling factor on the enrichment of 13 C in methane. Though the Taiyang shale gas field is close to the Qianzhong paleo-uplift and located on the southern slope of the depression, the sea-level fall in the Hirnantian stage was significant, the positive carbon isotopic excursion should also be important, and TOC could even be up to 8%, while the thickness of the Hirnantian deposit is small, generally 2-3 m, so the contribution from the precursor organic matter should not be ignored but cannot be controlling.
Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 861040 linear correlation is not necessarily caused by thermal maturity. When plotting the carbon isotope of shale gas from the Taiyang gas field versus their actual vertical burial depth, it is clear that methane, ethane, and propane become more enriched at 13 C, with increasing burial depth. When the burial depth changes from 2300 to~900 m, the carbon isotopic composition of methane, ethane, and propane all decreases by about~3‰, and a linear correlation of such carbon isotopic changes exists for methane (R 2 = 0.9223), ethane (R 2 = 0.9252), and propane (R 2 = 0.888) ( Figure 11). Due to the geothermal gradient, increasing burial depth usually represents the increase of geo-temperature and thus thermal maturity; hence, δ 13 C values of hydrocarbon gases usually increase with increasing burial depth. However, the Wufeng-Longmaxi shale in the Taiyang area underwent deep burial from Triassic to Jurassic, resulting in the high thermal maturity, and began to uplift in the Cretaceous, so deep burial does not necessarily mean high geo-temperature and shale with different burial depths may have very similar thermal maturity. For example, the Longmaxi shale in Well Y109 has a burial depth of 2162-2202 m and equivalent Ro% of 3.1%, while the Longmaxi shale in Well Y102 has a much shallower burial depth (732-741 m) but still has similar equivalent Ro% (3.0%) (Liang et al., 2021). Therefore, there must be other post-genetic processes that cause the linear correlation between the burial depth and δ 13 C values of hydrocarbon gases. Though the Wufeng-Longmaxi shale in the Taiyang area is typically characterized by shallow burial depth (can be less than 1,000 m), diffusive gas leakage might occur in these shallow layers. In the Taiyang area, the decreasing carbon isotopic compositions of hydrocarbon gas with decreasing burial depth cannot be caused by diffusive gas leakage. In general, diffusive gas leakage is more likely to occur in shallow layers. During the leakage, since 12 C-12 C bonds are weaker than 12 C-13 C bonds and leak faster, the leaking gas will be more enriched in methane and depleted in 13 C, and the residual gas which is still trapped in the reservoir will be more depleted in methane but enriched in 13 C. If gas leakage occurs in the shallow layers, gases in the shallow layers will be more enriched at 13 C compared to those of the deep layers. Apparently, this is opposite to the case in the Taiyang shale gas field.
Carbon isotopic fractionation of methane due to mass transport has been evaluated in many studies, especially in low-permeability shale. During gas transport, 12 C-12 C bonds are weaker than 12 C-13 C bonds and transport faster; thus, it caused a significant isotopically lighter carbon than the source (Pernaton et al., 1996;Zhang and Krooss, 2001). Previous studies documented a large isotopic fractionation during gas migration under a nonsteady state (Prinzhofer and Pernaton, 1997;Zhang and Krooss, 2001). Zhang and Krooss (2001) proposed that the magnitude of depletion in 13 C is the highest during the initial nonsteady state of the diffusion process and decreases to a constant difference when it approached the steady state.
Canister desorption experimental study about the Longmaxi shale samples at reservoir temperatures showed a rapid enrichment of 13 C up to 13.7‰-16.2‰ in methane as desorption proceeds (Ma et al., 2020). Li et al. (2020) proposed to divide the gas transport process during a complete production into four stages and demonstrated that methane became depleted in 13 C during the transition stage (stage II) and became enriched in 13 C during the adsorbed gas desorption stage (stage III), while became depleted in 13 C again during the diffusion stage (stage IV). Carbon isotopic enrichment of methane due to adsorption/desorption during stage IV might contribute to the heavy carbon isotope of methane in Taiyang and Jiaoshiba shale gas fields. However, according to the analytical results and numerical calculations, it is suggested that under geological conditions, post-genetic coupled diffusion and adsorption/desorption due to gas transport normally cause isotopic fractionation of less than 5‰ (Xia and Tang, 2012). In the Taiyang shale gas field, once diffusive migration occurs, the separates will be more enriched in methane and depleted in 13 C, resulting in the relatively low δ 13 C values in the shallow layer, while the residual will be more depleted in methane and enriched in 13 C, resulting in the relatively high δ 13 C values in the deep layer. The C 1 /(C 2 +C 3 )~δ 13 C 1 plot also shows that the Taiyang shale gas has a rapid change in molecular composition with burial depth, indicating potential effects from diffusive migration ( Figure 11).
The geological conditions of the Taiyang gas field also support possible diffusive gas migration in this area. The main part of the Taiyang shale gas field is the Taiyang anticline, with a gas-bearing area of 580 km 2 . It has suffered alteration from the superimposed multistage tectonic activities since the Caledonian period. The Taiyang anticline is characterized by "strong superimposed fold deformation and weak multistage fault reconstruction," and develops nearly east-west reverse faults and nearly south-north strike-slip faults. The nearly east-west faults are developed on the top and two wings of the anticline due to the north-south compression stress in the early stage. With continuous compression stress from north-south, the complex anticline gradually deforms and forms nearly north-south strike-slip faults. The uplift and denudation stage of extrusion strike-slip in the Late Yanshanian-Himalayan period resulted in the shallow burial depth of Wufeng-Longmaxi formations, that is, 500-1,500 m in the main part of the anticlinal structure (Liang et al., 2021). The formation pressure coefficient is 1.22-1.47, natural micro-cracks are developed in the shale, and the porosity is 4%-9%, with an average of 5.67% (Liang et al., 2020a). Therefore, diffusive migration of shale gases is likely to occur during or after the uplift and denudation stages.

Depletion of 13 C in Ethane and Propane
At low mature and mature stages, δ 13 C values of CH 4 are much lower than those of C 2 H 6 and δ 13 C values of C 2 H 6 are much lower than those of C 3 H 8 , that is, δ 13 C 1 > δ 13 C 2 > δ 13 C 3 . When the source rocks reach the oil window with Ro% of 1.5%, the carbon isotopic difference between methane and ethane and between ethane and propane will approach a similar value of −6‰ (Cesar et al., 2020). At this stage, isotopic effects of thermodynamic/equilibrium are predominant controlling factors, and isotopic exchange reactions occur under these conditions (Cesar et al., 2020). When the thermal maturity further increases, isotopic effects of thermodynamic/equilibrium are no longer the predominant controlling factors, and other processes may become more important; thus, carbon isotopic reversal first occurs in ethane and then in propane (Dai et al., 2016a;Cesar et al., 2020). In general, according to the published data, it is clear that with increasing thermal maturity, ethane and propane first become more enriched in 13 C, then become more depleted in 13 C, and finally become more enriched in 13 C again (Ferworn et al., 2008;Burruss and Laughrey, 2010;Rodriguez and Philp, 2010;Zumberge et al., 2012;Hao and Zou, 2013;Tilley and Muehlenbachs, 2013;Dai et al., 2014a;Dai et al., 2016b;Dai et al., 2017;Liu et al., 2018;Ni et al., 2018;Xia and Gao, 2018;Feng et al., 2020) (Figures 7B, C). The Wufeng-Longmaxi shale gases from the Taiyang, Jiaoshiba, and Weirong shale gas fields fall in the post-rollover area, where ethane and propane are still depleted at 13 C, but the carbon isotopic compositions get to be more positive again (Figures 7B, C). Among the mechanisms which have been invoked to account for the carbon isotopic depletion in ethane and propane, most of them are related to high temperatures, for example, 1) in situ cracking of C 2+ gas (Ferworn et al., 2008), in situ cracking of oil or gas (Rodriguez and Philp, 2010), simultaneous cracking of kerogen, oil, and gas in a closed system (Tilley et al., 2011), and cracking of kerogen-entrapped straight-chain aliphatic hydrocarbons ; 2) water reacting with CH 4 to form isotopically light CO 2 and hydrogen and CO 2 and hydrogen, subsequently forming isotopically light ethane and propane (Zumberge et al., 2012); 3) Rayleigh fractionation of ethane and propane involving redox reactions with transition metals and water at temperatures on the order of 250°C-300°C (Burruss and Laughrey, 2010); 4) Diffusive migration during the stage with high maturity (Dai et al., 2016a). As discussed above, Fischer-Tropsch-type synthesis of hydrocarbon from CO 2 and H 2 has contributed to the carbon isotopic depletion in ethane and propane, but due to the limited amount of the reactant, it cannot be the controlling factor. In the southern Sichuan Basin, wet gas cracking cannot explain the depletion of 13 C in ethane and propane. C 2+ wet gas cracking at elevated temperatures will increase C 2 /C 3 ratio and thus a decrease in gas wetness (Prinzhofer and Huc, 1995). The Wufeng-Longmaxi shale gas is characterized by very low gas wetness, and the content of ethane and propane is generally less than 1. Wet gas cracking of pentanes and butanes will cause enrichment in 12 C in ethane and propane (Xia et al., 2013); however, under the over-mature dry gas stage, cracking of ethane and propane possibly occurs and the residual ethane and propane will become more enriched in 13 C due to the kinetic isotopic fractionation effects (Tang et al., 2000). Obviously, this is opposite to the phenomenon of the Wufeng-Longmaxi shale gas.
Actually, Xia and Gao (2018) stated that mechanisms such as mixing of primary and secondary products in source rocks (Rodriguez and Philp, 2010;Ni et al., 2018;Hao and Zou, 2013;Ferworn et al., 2008;Xia et al., 2013) and water-involved reactions (Burruss and Laughrey, 2010;Gao et al., 2014;Zhang et al., 2018b;Zumberge et al., 2012) cannot satisfy the isotopic and mass balance requirements and proposed that the depletion of 13 C in ethane and propane resulted from the reversible conversion from alkane to alkyl groups. The isotope effect is dependent on the rate and kinetic isotope effects of these reversible free-radical reactions catalyzed by minerals, and the isotopic depletion in 13 C is consistent with the extent of the decomposition (Xia and Gao, 2018). However, these free-radical reactions may proceed in very different ways in different locations since the free-radical concentrations may vary greatly in different geological environments. In the sandstone and limestone formations, there is no sufficient hydrogen donors, and the reversal reactions will be suppressed and the fractionation will proceed likely in a Rayleigh fractionation mode; while organic-rich shales have abundant hydrogen donor, so the reversible reaction is possible, and the isotopic reversal is commonly found in either conventional and unconventional natural gases at relatively high thermal mature stages (Xia and Gao, 2018). The Upper Ordovician Wufeng-Lower Silurian Longmaxi shale comprises black shale, siliceous shale, calcareous-siliceous shale, and argillaceous-siliceous shale (Zou et al., 2019). It is deposited in a deep water environment and dominated by Type I and II 1 organic matter and typically has TOC of around 3% in the Taiyang, Jiaoshiba, and Weirong areas (Dai et al., 2014a;Liang et al., 2020a;Liang et al., 2020b;Zou et al., 2021). During the Hirnantian Stage glacial interval (Middle and Upper Wufeng Formation), there was a rapid global fall in sea level, while the P 2 O 5 /TiO 2 ratio reached 0.84 and the southern Sichuan Basin still experienced plankton bloom, high paleo-productivity, and high burial rates of organic matter (Zou et al., 2019). During the Rhuddanian Stage (Lower Longmaxi Formation), sea level increased and plankton bloomed, resulting in the widespread deposition of back shale in deep, oxygen-poor conditions under high-productivity surface waters (Zou et al., 2019). In southern China, Silurian marine source rocks with restored hydrogen index greater than 600 mg/g TOC account for 43%, and samples with restored hydrogen index between 250-600 mg/g TOC account for 57% (Zhang et al., 2013). In addition, the Wufeng-Longmaxi shale reached over 180°C in Late Permian in the Taiyang area and reached over 200°C in Late Jurassic in the Jiaoshiba area (He et al., 2017;Liang et al., 2020b;Pang et al., 2019;Tenger et al., 2017;Xiong et al., 2021). In general, the Wufeng-Longmaxi shale has equivalent vitrinite reflectance Ro% of 2.73%-4.62% Zou et al., 2021) and gas wetness between 0.09% and 0.97%. Therefore, the Wufeng-Longmaxi shale is abundant in hydrogen donors and at a high thermal maturation stage with a high conversion rate. When plotting the Wufeng-Longmaxi shale gases into the δ 13 C 2 value versus gas wetness plot of the Barnett shale gases by Xia and Gao (2018), the carbon isotopic fractionation of ethane of the Wufeng-Longmaxi shale gases agrees well with the calculated results (pink line in Figure 12), which mainly fall in the post-rollover stage and indicates high thermal maturity of the gas. As indicated by the reversible free-radical reactions, the isotopic compositions are positive at the stage of low conversion ratios and then become negative with increasing conversion ratios; When the conversion ratios get even higher and the remaining reactant becomes even less, the isotopic compositions become more positive again (Xia and Gao, 2018). As for the Wufeng-Longmaxi shale gas, the extremely low content of ethane and propane (C 2+3 < 1%) indicates high conversion ratios, so during the late catagenesis stage, decomposition of ethane and propane is accompanied by the depletion of 13 C in the residual ethane and propane, but in general, the isotopic composition becomes more positive again and falls in the post-rollover area. During this over-mature stage, reactions of free-radical decomposition are important controlling factors, resulting in the depletion of 13 C in ethane and propane, which prevail the thermal cracking reactions.

CONCLUSION
Molecular and stable carbon isotopic compositions of gases from the Wufeng-Longmaxi shale in the Taiyang (shallow), Jiaoshiba (middle), and Weirong (deep) shale gas fields in the Sichuan Basin have been investigated. All the gases are dominated by methane and have very low gas wetness (generally <0.83). Methane accounts for 98.82%, 98.13%, and 96.91% in the Taiyang, Jiaoshiba, and Weirong shale gas fields, respectively. The average carbon isotopic value of methane is −35.2‰ in Weirong, −30.3‰ in Jiaoshiba, and −28.5‰ in Taiyang shale gas fields. The enrichment of 13 C in methane is mainly caused by the extremely high thermal maturity. Shale gases from the Jiaoshiba and Taiyang gas fields fall offset from the δ 13 C 1~w etness linear trend line by the Weirong, Barnett, and Fayetteville shale gases, which is due to the contribution from the Fischer-Tropsch-type synthesis of hydrocarbon gas from CO 2 and H 2 . When thermal cracking gas is mixed with the hydrocarbon gas generated by the Fischer-Tropsch-type synthesis, it will significantly increase the amount of wet gas, but has little impact on methane. In the Taiyang shale gas field, when burial depth changes from~2300 tõ 900 m, δ 13 C 1 , δ 13 C 2, and δ 13 C 3 values decrease by~3‰, and a linear correlation exists between carbon isotope and burial depth. This is mainly caused by diffusive migration. Shale gases from the Taiyang, Jiaoshiba, and Weirong gas fields are characterized by a complete carbon isotopic reversal trend among the C 1 -C 3 alkane gases (δ 13 C 1 > δ 13 C 2 > δ 13 C 3 ). This is mainly caused by the reversible free-radical reactions with the conversion from alkane to alkyl groups with a minor contribution from the Fischer-Tropsch-type synthesis.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding authors.