Geochemical characteristics and genetic mechanism of the high-N2 shale gas reservoir in the Longmaxi Formation, Dianqianbei Area, China

As an important pilot target for shale gas exploration and development in China, the Longmaxi Formation shale in the Dianqianbei Area is characterized by high content of nitrogen, which severely increases exploration risk. Accordingly, this study explores the genesis of shale gas reservoir and the mechanism of nitrogen enrichment through investigating shale gas compositions, isotope features, and geochemical characteristics of associated gases. The high-nitrogen shale gas reservoir in the Longmaxi Formation is demonstrated to be a typical dry gas reservoir. Specifically, the alkane carbon isotope reversal is ascribed to the secondary cracking of crude oil and the Rayleigh fractionation induced by the basalt mantle plume. Such a thermogenic oil-type gas reservoir is composed of both oil-cracking gas and kerogen-cracking gas. The normally high nitrogen content (18.05%–40.92%) is attributed to organic matter cracking and thermal ammoniation in the high-maturity stage. Specifically, the high heat flow effect of the Emeishan mantle plume exacerbates the thermal cracking of organic matter in the Longmaxi Formation shale, accompanied by nitrogen generation. In comparison, the abnormally high nitrogen content (86.79%–98.54%) is ascribed to the communication between the atmosphere and deep underground fluids by deep faults, which results in hydrocarbon loss and nitrogen intrusion, acting as the key factor for deconstruction of the primary shale gas reservoir. Results of this study not only enrich research on genetic mechanism of high-maturity N2 shale gas reservoirs, but also provide theoretical guidance for subsequent gas reservoir resource evaluation and well-drilling deployment in this area.


Introduction
Compared with conventional hydrocarbons, shale gas is essentially different in the occurrences, accumulation characteristics, and enrichment mechanisms under the low porosity and low permeability conditions (Wang and Li 2017;Wu et al. 2018;Chen et al. 2018a;Daniel 2018). Specifically, the Longmaxi Formation shale in the Dianqianbei Area is featured by over-high thermal maturity and high nitrogen content, which increases the difficulty and risk of shale gas exploration. Previous studies focus on the evaluation of shale gas richness from perspectives of sedimentary environment (Zhao et al. 2016;Deng et al. 2018), microscopic reservoir space (Hu et al. 2019;Wang et al. 2019) as well as reservoir-forming and diagenetic characteristics Wang et al. 2015). In comparison, this study investigates the shale gas accumulation process in the Dianqianbei Area using geochemical Edited by Jie Hao * Ji-Lin Li ljlgeology@163.com * Ting-Shan Zhang zts_3@126.com methods. Specifically, there are many key issues to be solved, including the genetic type of shale gas, the occurrence and enrichment law of shale gas, and the source of nitrogen in some high-nitrogen, low-carbon gas reservoirs. Elaboration of these issues will greatly promote the exploration and development process of shale gas in this area. It is generally easy to retain shale gas isotope distribution characteristics in the inherited simple basins such as Appalachian and Ordos Basins (Biryol et al. 2016;Feng et al. 2017), while it is much more difficult to reconstruct the shale gas isotope features in the foreland basins or multi-stage superimposed basins, such as the case in the Dianqianbei Area (Chen et al. 2018b), which makes it very challenging to distinguish shale gas genesis. Moreover, whether it is reliable to study high to over-high maturity shale gas using the traditional alkane isotope identification is still controversial (Feng et al. 2016;Liu et al. 2018).
Recently, much attention has been paid to the high-nitrogen marine shale reservoirs in southern China (Guo and Liu 2013;Liu et al. 2016). As an important indicator to reflect preservation conditions of oil and gas reservoirs (Chen 2005;Jiao et al. 2017), nitrogen can also help identify shale gas source and components ). However, too high nitrogen content will increase the exploration risk and affect the accuracy of subsequent gas reservoir evaluation and the reliability of the development plan. Existing research acknowledges many nitrogen sources, including atmosphere , biodegradation (Baxby et al. 1994), thermal decomposition (Wellman et al. 1968), cracking of nitrogen-containing minerals (Jenden and Kaplan 1989), and crust or mantle . So far, rare in-depth research has been conducted to explore the cause of the high nitrogen content in the shale gas reservoirs of the Dianqianbei Area. Accordingly, this study investigates the rule of the multistage alkane isotope reversal in the high to over-high maturity shale gas through investigating the geochemical characteristics of the Longmaxi Formation shale gas, accompanied by establishment of a relatively stable non-hydrocarbon gas characteristic chart that can be used to clarify gas genesis. Moreover, through implementing thermal simulation experiments, this study quantitatively identifies proportions of various types of gas components with different sources. Furthermore, regional tectonics and thermal history are combined with the above results to systematically explore the cause of the "high nitrogen content and low hydrocarbon content" in some shale gas reservoirs. Finally, the genetic mechanism of the high maturity shale gas reservoir is elaborated, accompanied by construction of a model illustrating the shale gas accumulation process in the Dianqianbei Area. Results of this study are expected to provide an important geological basis for subsequent exploration and development of shale gas reservoirs in this area.

Geological setting
As a multi-stage superimposed basin developed on the Upper Yangtze platform, the Dianqianbei Depression is adjacent to the southern margin of the Sichuan Basin in the north, the Dianqian uplift in the south (Fig. 1), and the Diandongbei Depression in the west (Chen et al. 2018c). The depression is featured by strong tectonics in the south and west, and it is tectonically a typical slotted east-west extended fold belt composed of wide anticlines and closed synclines (Liang et al. 2011;He et al. 2019). As the main target layer for shale gas exploration, the Silurian Longmaxi Formation was deposited in the epicontinental sea environment bounded by paleouplifts and local highlands, which resulted from compression of the Middle and Upper Yangtze Block. Generally, the Silurian Longmaxi Formation is a set of deep-water organicrich shale deposits formed in the constrained environment with depleted oxygen (Chen et al. 2018a;Li et al. 2019). This set of shale is about 200 m thick, with a TOC range of 0.19%-8.46% and a vitrinite reflectance (Ro) range of 2.2%-3.96%, which indicates high to over-high maturity that corresponds to the gas generation stage. The content of brittle minerals such as quartz and feldspar in shale is high, reaching 46%-84%. Large quantities of micropores and microfractures provide adsorption space and seepage channels during the shale gas accumulation (Wang 2017). Compared with geological features of large-scale commercial shale gas reservoirs such as those in fields of Barnett, Ohio, and Antrim in the USA (Fu et al. 2015;Zou et al. 2016;Jan et al. 2015), the Longmaxi Formation shale gas in the study area is of good resource basis and high exploration potential.

Samples and experimental methods
Shale gas samples were taken from four Longmaxi Formation shale gas industrial wells in the study area using high-pressure steel cylinders. Both the steel cylinders and the steel pipes can bear a pressure of 15 MPa. For highpressure wells, samples were taken at the wellhead separator, and gas samples of about 5 MPa were filled in the cylinders after multiple replacements.
The shale gas composition analysis was completed at the Research Institute of Petroleum Exploration and Development. A MAT271 gas mass spectrometer was used to analyze the non-hydrocarbon components in the shale gas and the absolute contents of CH 4 and C 2 H 6 . The hydrocarbon gas was analyzed by a GC-9160 gas chromatograph, and the final data were normalized by mass spectrometry and chromatographic data of CH 4 , which has an analysis error of smaller than 0.5%. Determination of carbonate isotope was completed by first separating the gas components using an Agilent 6890 gas chromatograph, then oxidizing corresponding components to CO 2 in the high-temperature furnace, and finally obtaining the value by a Delta Plus XP mass spectrometer. Analysis errors of CH 4 and C 2 H 6 were smaller than 0.5‰, and that of CO 2 was better than 1‰, with V-PDB as the calibration standard. Similarly, determination of hydrogen isotope involved separation of gas components by an Agilent 6890 gas chromatograph, reduction of components to H 2 in a hightemperature furnace, and acquisition of the hydrogen isotope by a MAT253 mass spectrometer. The analysis error was smaller than 0.5‰, and V-SMOW was used as the calibration standard. The helium isotope was analyzed by a Noblesse rare gas mass spectrometer, during which the two-stage separation system was employed for purification and enrichment. The 3 He/ 4 He analysis error was less than ± 1.5%. Nitrogen isotopes are measured on a MAT252 isotope mass spectrometer using the GC-C-MS method, which can achieve an accuracy of 0.5‰ with reference to V-AIR as the calibration standard.

Shale gas compositions
According to experimental tests of the Longmaxi Formation shale samples from Wells Y2 and B1 (Table 1   0.13%-98.75%. Except for Well Y2 that has high nitrogen content, all the other wells are seen with a small amount of ethane, with a content range from 0% to 2.16%, and only a few samples contain trace propane. The humidity coefficient (C 2+ /C 1 ) is between 0.23% and 2.33%, with an average of 0.98%, which indicates typical dry gas reservoirs. The non-hydrocarbon gas is mainly CO 2 (0%-8.95%) and N 2 (18.05%-98.54%).

Isotope compositions and reversal characteristics of alkanes in the shale gas
The methane carbon isotope (δ 13 C 1 ) varies from − 27.9‰ to − 20.8‰, with an average of − 23.6‰; the hydrogen isotope (δD) ranges from − 141‰ to − 104‰, with an average of − 129‰; the ethane carbon isotope (δ 13 C 2 ) varies between − 31.8‰ and − 29.3‰, averaging − 30.8‰; and the propane carbon isotope (δ 13 C 3 ) ranges from − 32.9‰ to − 27.8‰, with an average of − 29.3‰. Overall, δ 13 C 1 is much higher than δ 13 C 2 , while δ 13 C 2 shares a similar range with δ 13 C 3 (Table 2). Conventional natural gas is generally characterized by the monotonic positive sequence of carbon isotope. In comparison, the alkane isotope anomaly in shale gas is mainly reflected in the change of ethane carbon isotope (δ 13 C 2 ) with the increase in kerogen thermal maturity, presenting an anti-"S" shape curve. The δ 13 C 2 reversal that is defined by Tilley and Muehlenbachs (2013) refers to the phenomenon that δ 13 C 2 reaches the maximum at a reversal point, after which δ 13 C 2 starts to decrease. This reversal point is ascribed to the mixture of secondary cracking gas from liquid hydrocarbons with the pre-existing kerogen-cracking gas (Wu et al. 2015). A secondary reversal is defined at the moment that δ 13 C 2 becomes lighter than δ 13 C 1 , which is speculated to be caused by excessively high thermal maturity   (Cao et al. 2015). Samples from representative Wells Y8 and Y9 are all in the reversal parts at the end of the anti-"S" shape curve. This is attributed to the cracking of heavy hydrocarbons and early-formed crude oil in the high to overhigh maturity stage, which result in ethane with much larger volume and lighter isotopes than the contemporary ethane (content less than 2%), leading to deconstruction of the positive carbon isotope sequence (Fig. 3). On the other hand, it is also ascribed to the Rayleigh fractionation effect, which occurs in the reduction reactions between mature ethane and formation water or iron-containing metals that are catalyzed by transition metals under high-temperature conditions (Michelsen et al. 2015) (Fig. 4). In the Middle and Late Permian, the Dianqianbei Area was affected by the thermal effect of the Emeishan basalt super mantle plume Wu et al. 2014), which led to the presence of many high-density methane inclusions in the Longmaxi Formation shale, with the homogenization temperature reaching 205 °C. Moreover, the overlain strata on the Longmaxi Formation shale are Silurian water-bearing clastic rocks, which also favor the occurrence of Rayleigh fractionation. In this process, a small amount of in situ 13 C 2 H 6 was involved in the reaction and converted into 13 CH 4 , which resulted in 13 C 2 deficiency and 13 C 1 enrichment, forming the carbon isotope reversal (δ 13 C 1 > δ 13 C 2 ). Moreover, such reversal was further enhanced by the thermochemical reaction of residual kerogen and formation water.

Genesis of gas hydrocarbons in the shale gas
The carbon isotope distribution in alkane gas is commonly used to identify sources of conventional gas and coalbed  Fig. 4 Evolution process of carbon isotope in gas system with different maturity methane, which has achieved good results in gas reservoir exploration practice. However, it is not applicable in the high to over-high maturity shale gas reservoirs where alkane carbon isotopic reversal is common. Many identification charts for conventional isotopes are not accurate any more. For instance, when classified using the detailed Bernard chart calibrated by conventional natural gas reservoir samples (Chen et al. 1995), the gas in the Dianqianbei Area falls into the coal-forming gas category and its mixed source area, which is contradictory to the fact that the gas originates from marine sapropel shales (Fig. 5a). Therefore, the Bernard series charts are almost ineffective for genesis identification of high to over-high maturity natural gas. Compared with shale gas in the North America, the Longmaxi Formation shale gas in the Dianqianbei Area has abnormal heavy δ 13 C values and high C 1 /(C 2 + C 3 ) ratios in the high to over-high maturity stage. As shown in Fig. 5b, samples from Wells Y9 and Y8 fall in the vicinity of the lower boundary of the thermogenic zone in the δ 13 C(CO 2 ) and δ 13 C(CH 4 ) cross-plot, which is related to the high thermal maturity of the alkane gas. Moreover, the average content of CO 2 is 2.28%, and the δ 13 C(CO 2 ) is between − 10‰ and − 8‰, with an average of − 8.14‰ (Table 2). Therefore, it can be inferred that CO 2 should mainly come from cracking of organic matter, with some mixed inorganic gas resulting from dissolution and cracking of carbonate minerals. According to migration of alkane isotopes in the kerogen cracking process, δ 13 C(C 2 H 6 ) is mainly controlled by kerogen type. The δ 13 C(C 2 H 6 ) value of the Longmaxi Formation shale

Fig. 5
Diagram for identification of genetic types of Longmaxi Formation shale gas components in the Dianqianbei Area. a Modified Bernard diagram about relationship between δ 13 C and C 1 /(C 2 + C 3 ), base map from reference (Chen et al. 1995); b relationship between δ 13 C-CO 2 and δ 13 C-CH 4 , base map from reference (Kotarba 2001;Liu et al. 2016); c relationship between δ 13 C 1 , δ 13 C 2 , and δ 13 C 3 , base map from reference (Dai et al. 2014;Liu et al. 2016) gas in the study area is lower than − 28‰, which indicates sapropel-related oil-type gas. As shown in Fig. 5c, all the shale gas samples fall in the oil-type gas zone in the δ 13 C 1 -δ 13 C 2 -δ 13 C 3 gas source chart and are close to the isotope reversal zone. This is in line with above-stated results about the negative carbon isotope sequence of the secondary reversal, namely δ 13 C 1 > δ 13 C 2 > δ 13 C 3 . The unconventional carbon isotope partitioning of the Longmaxi Formation shale gas is partly caused by secondary cracking of pre-existing liquid hydrocarbons in the shale during the high maturity stage. According to thermal simulation results of the Longmaxi Formation kerogen and normal C 22 alkane (crude oil equivalent) by the previous study , considering the fractionation mechanism of thermogenic gas isotope (Wei et al. 2015), linear trend models are, respectively, established for δ 13 C 1 and δD features of cracking products of crude oil and kerogen (Qu 2015) (Fig. 6).
The two nearly parallel trend lines in the isotope crossplot are located just in the vicinity of the upper and lower boundaries of the thermogenic area. A vertical line is drawn through the sample point, resulting in two intersection points with the two trend lines. These two points are used to, respectively, represent methane carbon isotope end values of cracking gas of crude oil and kerogen, denoted as δ 13 C oil and δ 13 C kero , respectively. Assuming the (1) Crude oil-cracking gas: δ 13 C = 0.1316δD − 22.405 R 2 = 0.9087 (2) Kerogen-cracking gas: δ 13 = 0.1048δD − 6.269 R 2 = 0.9190 proportion of oil-cracking gas is X, the methane carbon isotope of the mixture of two end members (δ 13 C mix ) can be expressed as follows: As shown in Table 3, the proportion of oil-cracking gas ranges between 9.32% and 28.90% for the Longmaxi Formation shale gas in the study area, with an average of 16.7%. As the thermal maturity increases to 2.0%, early-formed heavy hydrocarbons are fully cracked, contributing to the largest proportion of oil-cracking gas. With further increased maturity, wet components in the oil-cracking gas are continuously cracked, and no massive oil-cracking gas is formed, which result in the decreasing proportion of oil-cracking gas.

Genesis of high-content nitrogen
It is critical to investigate genesis and enrichment of N 2 in shale gas reservoirs so as to better understand hydrocarbon generation processes and hydrocarbon resource potentials of mature shale. Moreover, it has practical significance for predicting gas components in gas reservoirs and reducing exploration risks. In terms of gas components in the shale gas reservoirs of the study area, N 2 and CH 4 are dominant in most samples, with a few alkanes, rare gases, and CO 2 . As demonstrated in isotopic experiments, δ 15 N in the shale gas reservoir of the Dianqianbei Area greatly varies between − 2.09‰ and 7.65‰ (Table 4), and such large variation range suggests highly diversified types of N 2 genesis.
The conventional high nitrogen gas reservoir in Well B1 is representative of a certain type of gas reservoir, with the N 2 content ranging between 18.05% and 40.92% (averaging 27.54%). As shown in the δ 15 N-N 2 cross-plot (Fig. 7a), most samples from Well B1 have a δ 15 N range of 4‰-8‰, with an average of 4.77‰. According to the experimental thermal simulation results, N 2 resulting from organic matter cracking generally has a δ 15 N range between 4‰ and 20‰, which becomes heavier as the maturity increases (Kotarba and Lewan 2013). Therefore, it can be preliminarily inferred that the high N 2 content in the gas reservoir of Well B1 is mainly ascribed to organic matter cracking. According to the δ 15 N-δ 13 C 1 cross-plot (Fig. 7b), the δ 13 C 1 value is between − 30.5‰ and − 25.8‰. As the δ 13 C 1 value grows, δ 15 N becomes heavier, presenting features of organic matter cracking during the high to over-high maturity stage.
During the Permian (about 259 Ma ago), the Southwestern Sichuan Region experienced the Emeishan basalt eruption event, and the Dianqianbei Area was in the middle to outer zone of the Emeishan mantle plume (Fig. 8), which provided thermodynamics for maturation of the Longmaxi Formation source rocks. In order to evaluate the influence of the mantle plume on the regionally differential thermal (3) δ 13 C mix = Xδ 13 C oil + (1 − X)δ 13 C kero  Fig. 6 Relationship between methane carbon and hydrogen isotopes and distinguishing model of mixed gas components (Qu 2015).
(Thermal simulation results of kerogen and C 22 are from Qu et al. 2015.) evolution of source rocks, two typical wells located in different zones are selected, namely the Well H1 near the middle zone of the mantle plume and the Well CY84 outside the mantle plume zone. Subsequently, thermal evolution of source rocks is correlated with the burial depth (Fig. 9). Compared with source rocks in Well CY84 that experienced   Fig. 8 Geology of the Emeishan large igneous province and the boreholes in Dianqianbei Area. Base map modified from reference (Xu et al. 2004) normal thermal evolution, those in Well H1 underwent ultra-high temperature caused by the mantle plume thermal activity in a short period of time, reflected in the rightward jumping of the actual maturity measurement from the normal trend line. Maturities of shales in Wells Y9, Y8, and B1 are between 1.96% and 3.42%, which are higher than those in the normal burial conditions. During the overmature stage (R o > 2%), organic matter is affected by the heat flow of the mantle plume, which results in full cracking and denitrification, accompanied by generation of massive N 2 . Moreover, the thermal evolution history of the Longmaxi Formation source rocks is established (Fig. 10) (Rao et al. 2013), according to the T-t thermal history path and the chemical kinetic model (Easy %R o ). Specifically, the Well CF82 in the Northern Sichuan Area is far from the mantle plume-affected zone, and thus Silurian source rocks there are characterized by gradual and continuous thermal evolution, becoming mature at the end of the Permian and finalized in the Cretaceous, with the highest paleotemperature of 180 °C. In comparison, the Well LS1 in the Southern Sichuan Area is located in the middle of the mantle plume-affected zone, and Silurian source rocks there were immature during the Caledonian period, but quickly entered the overmature stage during the Permian due to the anomalous thermal effect of the mantle plume, which resulted in a sharp temperature increase of 250 °C and a short oil window of only 15 Ma. The maturity of Silurian source rocks in the Well LS1 finally reached 3.1% in the Late Permian. Well B1 shares similar thermal history evolution with Well LS1, and source rocks there were quickly heated to be more than 300 °C during the Late Permian, meeting the chemical kinetic conditions of organic matter cracking and high N 2 production . A few samples from Well B1 fall out of the thermal cracking genesis in the identification chart. However, they are less likely to originate from volcanic activity and the deep part of the upper mantle or from thermal cracking of nitrogencontaining clay minerals. The 3 He/ 4 He ratio of the Longmaxi Formation shale is between 1.86 × 10 −8 and 3.96 × 10 −8 (equivalent to R/Ra ratio of 0.01-0.03), which is much lower than the critical condition of mantle-originated helium (R/ Ra > 0.1). Moreover, NH 4 + ions in nitrogen-rich clay minerals are of strong thermal stability, which can only be decomposed when the temperature exceeds 1000 °C (Whelan et al. 1988). In that case, generated N 2 needs to survive the differentiation process during migration, before it can be accumulated as the high-N 2 reservoir. Therefore, it can be concluded that the high-nitrogen Longmaxi Formation shale reservoir in the study area is not a product of inorganic mineral thermal cracking. In this study, it is considered to be related to thermal ammoniation in the high maturity stage of organic matter. When the activation energy required for the amino group (-NH 2 ) is met, NH 3 is produced, which subsequently reacts with CO 2 and Fe 2 O 3 , generating N 2 with a δ 15 N range of − 10‰ to − 1‰. Meanwhile, the associated methane has a δ 13 C range of − 55‰ to − 30‰. These features are consistent with the trend line in the δ 15 N-δ 13 C cross-plot  (Rao et al. 2013) ( Fig. 7b). To conclude, the high content of nitrogen in the shale gas reservoir of Well B1 is mainly ascribed to thermal cracking of organic matter and partly to thermal amination. The abnormally high nitrogen gas reservoir in Well Y2 is representative of another type of gas reservoir, with the N 2 content ranging between 86.79% and 98.54% (averaging 93.27%) and the δ 15 N value between − 2.09‰ and 0.79‰ (Fig. 7a). These features are basically consistent with geochemical characteristics of the atmospheric N 2 (δ 15 N ≈ 0). Through investigation into tectonic deformation in the Dianqianbei Area, it is found that the Indosinian movement in the Middle and Late Triassic is accompanied by strong tectonic compression that results in a series of massively deformed tightly closed fault belts. Specifically, located at the west flank of the Niujie Anticline in the Dianqianbei Area, the Well Y2 is at the flank of the gentle anticline, where inherited active deep faults are present. These faults enhance the circulation between atmospheric leaching water on the ground surface and fluids deep in the subsurface, adversely affecting preservation of early formed shale gas reservoirs. In addition, those main faults that cut through the basement are associated with derivative secondary faults, forming densely distributed fault belts that subsequently act as primary channels for hydrocarbon gas loss and external fluid intrusion. The alkane gas in the Lower Triassic shale gas reservoirs gets washed multiple times and finally releases along those deep faults to the surface. In this context, methane with heavier δ 13 C is more soluble in water and thus preferentially washed away, leading to light δ 13 C signatures (− 36.7‰ to − 29.5‰) of methane in the shale gas reservoir (Fig. 7b). Meanwhile, the downward infiltrating atmospheric water enters the gas reservoir with carried N 2 , resulting in a relative large proportion of atmospheric N 2 and a relative small proportion of thermogenic N 2 as observed in Well Y2.

Shale gas formation mechanism
There are mainly three episodes of tectonic movements after the deposition of the Lower Silurian Longmaxi Formation shale in the study area, namely the Caledonian, the Hercynian-Indosinian-Yanshanian, and the Sichuan-Himalayan. During the Caledonian, marine organic-rich shale rapidly subsided, accompanied by gradually growing maturity. During the immature stage (R o < 0.55%), organic matter was decomposed by microbes to form a small amount of biogas. Then, the source rocks entered the low maturity stage (0.55% < R o < 0.7%) and thus the hydrocarbon generation threshold during the Late Caledonian (420 Ma), associated with oil generation. The hydrocarbon generation process was intermittently interrupted by the Hercynian uplifting. Then, source rocks entered the oil window in the Late Hercynian (280 Ma ago), which was accompanied by generation of massive CO 2 in addition to liquid hydrocarbons. At the end of the Hercynian during the Middle and Late Permian, the Emeishan basalt erupted. The continuous, rapid burial and the high anomalous heat flow of basalt magma exacerbated the maturation of hydrocarbon-generating organic matter, contributing to significantly earlier (90-110 Ma) occurrence of hydrocarbon generation peaks. Moreover, the greatly shortened oil window was followed by an extended gas formation process. In the early stage of the Indosinian Movement, organic matter entered into a high maturity stage (R o > 1.3%), associated with thermal cracking of residual kerogen and liquid hydrocarbons to wet gas, gradually decreasing CO 2 yield, and generation of a small amount of N 2 by thermal ammoniation. During this period, the δ 13 C 1 value was relatively stable. With formation of more and more oil-cracking gas, more light δ 13 C ethane was generated, and carbon isotope reversal occurred when R o reached 1.8%. In the middle of the Indosinian, all heavy hydrocarbon components such as crude oil were cracked into gaseous hydrocarbons, and wet gas started to crack when R o reached 2.0%. During the Yanshanian, the study area was in a stable subsidence stage, when the carbon isotope reversal reached the maximum. At that time, wet gas was cracked in a large amount, and the yield of organic matter cracking-related nitrogen gradually increased, approaching 50% of the hydrocarbon generation rate. Influenced by both the long-term subsidence and thus deep burial as well as the regional basalt mantle plume, the Longmaxi Formation shale reached the temperature conditions required for Rayleigh fractionation, accompanied by partial conversion of 13 C 2 H 6 to 13 CH 4 and maintained carbon isotope reversal. Subsequently, during the Sichuan-Himalayan tectonic cycles, stably distributed Longmaxi Formation shale becomes elevated due to rapid uplifting. Moreover, junctions of anticlines and synclines in the strongly compressed fold bet are subjected to development of deep faults. Furthermore, the flank of an anticline is vulnerable to tectonic uplift and denudation. All these factors result in poor sealing and storage conditions in the peripheries of deep faults and the tops of shallow buried anticlines, leading to loss of early-formed hydrocarbon gas and intrusion of atmospheric gas (Fig. 11).
In this study, it is suggested to first locate relatively lowmaturity shale gas targets in the Longmaxi Formation of the study area. Moreover, the central and northern parts of the study area, where tectonics are weak and deep faults are absent, are favorable for shale gas enrichment, and thus drilling risk of meeting low-abundance, high-nitrogen shale gas reservoir is effectively reduced there.

Conclusions
The high-nitrogen shale gas reservoir of the Longmaxi Formation in the Dianqianbei Area is a typical dry gas reservoir. The complete carbon isotope reversal (δ 13 C 1 > δ 13 C 2 ) is speculated to result from the secondary cracking of earlyformed liquid hydrocarbons in the high to over-high maturity stage as well as the high-temperature Rayleigh fractionation induced by the mantle plume heat flow. According to isotope migration features of alkanes and CO 2 , the Longmaxi Formation shale gas is identified as a type of sapropel-related thermogenic oil-type gas. Through application of the thermal simulation model of crude oil and kerogen, proportions of various types of gases in the shale gas are determined. As the thermal maturity increases, the proportion of oil-cracking gas first increases and then decreases. As for the normally high nitrogen content (18.05%-40.92%), it is believed to result from organic matter thermal cracking and thermal ammoniation in the high-maturity stage. Specifically, the high heat flow effect of the Emeishan mantle plume exacerbates the thermal cracking of organic matter in the Longmaxi Formation shale, accompanied by nitrogen generation. In comparison, the abnormally high nitrogen  Fig. 11 Thermal history and evolution model of Longmaxi Formation shale gas in Dianqianbei Area content (86.79%-98.54%) is ascribed to the communication between the atmosphere and deep underground fluids by deep faults, which results in hydrocarbon loss and nitrogen intrusion. Influenced by the basalt mantle plume, the shale gas accumulation process is featured by accelerated organic matter maturation, much earlier (90-110 Ma) arrival of the hydrocarbon generation peak, entrance into the wet gas generation stage and subsequent cracking stage during the Early Indosinian, and occurrence of carbon isotope reversal during the Yanshanian continuous subsidence.
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 holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.