1 Introduction
In recent years, the development concept “green mountains are gold mountains” has been widely accepted in China, and the importance of natural gas in the energy consumption structure has been increasingly improved. To date, significant achievements have been made in terms of the generation mechanism, origin discriminant, distribution of gas reservoirs, and migration and accumulation of natural gas (
Song et al., 2014;
Li et al., 2018,
2020;
Wang and Li, 2020;
Li et al., 2021a), which greatly facilitate the exploration of natural gas in China, such as in the Tarim Basin (
Liu et al., 2018), Sichuan Basin (
Wu et al., 2016), Junggar Basin (
Hu et al., 2018;
Chen et al., 2019;
Li et al., 2019a).
The Jiyang Depression has been explored and developed for more than 50 years and still shows huge exploration potential (
Vincent et al., 2011;
Zhang et al., 2011;
Cao et al., 2016). Although the Jiyang Depression possesses the second largest oilfield in China (Shengli Oilfield), its associated gas has been identified as small-scale and scattered. In the Jiyang Depression, the natural gas reservoirs are distributed in multiple sets of gas-bearing layers and occupy different tectonic positions (
Zhang, 1991;
Gan et al., 2005;
Hu et al., 2021a). Previous studies have mainly focused on the genetic types, accumulation characteristics, enrichment model, and migration pathways of hydrocarbons in Minfeng Subsag (
Li et al., 2005;
Yang et al., 2014;
Ping et al., 2017;
Ci et al., 2020), Gubei-Bonan Subsags (
Lin et al., 2007;
Wang et al., 2007a;
Guo et al., 2009;
Jiang et al., 2009,
2010;
Zhang et al., 2009), and Huagou-Pingfangwang-Pingnan Area (
Hu et al., 2009;
Meng et al., 2015), that is, most studies have focused on limited areas with certain gas types. Based on the large amount of analysis data from 30 gas fields of different scales, this study attempted to improve the understanding of genetic type, distribution, and main controlling factors of natural gas discovered in the Jiyang Depression in order to facilitate further natural gas exploration.
2 Geological setting
Jiyang Depression, located in the southwest of the Bohai Bay Basin, is the biggest sub-unit of the Bohai Bay Basin (Fig.1). It is the main petroleum-producing base for Shengli Oilfield. From 1960s to the present, the Jiyang Depression has proven crude oil reserves of about 50 × 10
8 t, and proven natural gas reserves of 2500 × 10
8 m
3, total production of oil is 10.7 × 10
8 t, and the natural gas is 460 × 10
8 m
3 (
Li et al., 2003;
Zhang, 2012;
Song and Li, 2020). The Jiyang Depression is separated from Ludong Uplift in the east, and it is bounded by the Luxi Uplift in the south, and the Qingcheng-Linpanjia-Binxian highs in the west, the Chenjiazhuang high in the north. Many sub-depressions and highs were developed in the depression, covering from the southwest to the northeast (
Chen et al., 2009). The sub-depression, i.e., Zhanhua, Dongying, Huimin, and Chezhen, respectively, can be further divided into different ternary structural units, which all showed good oil and gas generation and accumulation conditions (Fig.1).
The tectonic evolution history of the Jiyang Depression comprised of three main stages. 1) Crystalline basement forming stage: the basement, Archean Taishan group, is a set of metamorphic rocks formed after Taishan Movement. 2) Platform sediments development stage (Palaeozoic-Mesozoic): in this stage, the depression experienced intensive uplifting, and marine and continental sediments alternately deposited. 3) Faulted depression-depression development stage: during this stage, the Jiyang Depression experienced Yanshanian movement and Himalayan movement successively, and developed a number of faults, accompanied by multiple episodic magmatic activity (
Zhang, 2012;
Zhang et al., 2020). The faults in the Jiyang Depression can be divided into three groups according to their occurrence, i.e., NE-, NW-, and near-E trending. Among them, NW-trending faults are mainly basin-controlling faults in the Late Jurassic–Early Cretaceous and Early Cenozoic. However, NE and near-E Faults are mainly active in the Cenozoic and Paleogene, especially in the Middle Eocene (
Cai, 2008).
In terms of oil and gas forming conditions, five sets of high-quality gas source rocks were found in the Jiyang Depression, such as Carboniferous–Permian, Paleogene Kong 2 Member (E
k2), Sha 4 Member (E
s4), Sha 3 Member (E
s3) and Sha 1 Member (E
s1). These source rocks possessed different organic matter type, burial depth, and thermal evolution, which make the types of gas generated complex and diverse. The natural gas reservoirs consisted of both hydrocarbon gas and non-hydrocarbon gas (mainly carbon dioxide). At present, 13 industrial gas-bearing layers have been discovered (
Su et al., 2009;
Li et al., 2015) covering Archaean, Paleozoic, Mesozoic, and Cenozoic intervals (Fig.2). From the perspective of the genesis of natural gas, Jiyang Depression possessed oil-associated gas, coal-derived gas, and biogas, as well as inorganic mantle-derived carbon dioxide gas (
Yang et al., 2019). The existence of various types of natural gas makes the Jiyang Depression become a promising exploration area.
3 Samples collection and methods
In this study, 139 natural gas samples from 29 gas fields in the Jiyang Depression were geochemically investigated. All samples were adopted from the database of the Shengli Oilfield Company and published literature (
Song and Zhang, 2004;
Zhu et al., 2005;
Lin et al., 2007;
Shen et al., 2007;
Wang et al., 2007a;
Wang et al., 2007b;
Wang, 2008;
Guo et al., 2009;
Hu et al., 2009;
Wang et al., 2009;
Zhang et al., 2009;
Li et al., 2010;
Yang et al., 2014;
Meng et al., 2015;
Li W et al., 2015;
Li Y et al., 2017). The natural gas samples involved in this study cover almost all the main gas production fields and formations in the Jiyang Depression. The characteristics of source rocks (including E
s and Carboniferous–Permian) associated with natural gas samples are summarized in Tab.1 (
Zhou et al., 2006;
Zhu et al., 2006;
Wang et al., 2010;
Zhang et al., 2011;
Gao et al., 2020), mainly including total organic carbon (TOC), organic matter (OM), and vitrinite reflectance (
Ro).
The natural gas compositions were measured using an Agilent 6890N gas chromatograph (GC). The stable carbon isotope compositions of the gases were measured using a Finnigan MAT-252 instrument. The TOC was measured using a LECOCS-230 analyzer after removing carbonate and washing with distilled water to remove residual hydrochloric acid (HCl). The Ro data were obtained using an oil immersion lens and a Leitz MVP-ST reflected-light microscope.
4 Results
4.1 Gas components characteristics
The results of 109 hydrocarbon gas samples from 29 gas fields are summarized in Tab.2 and Tab.3, and the results of 30 carbon dioxide (CO
2) gas samples are summarized in Tab.4. As shown in Tab.1 and Tab.2, for the hydrocarbon gases, the methane (CH
4) content ranges from 59.60% to 99.26%, with an average of 87.41%, and the content of heavy hydrocarbons (C
2+) ranged from 0.11% to 38.79%, averaging 8.33%. The other gases include CO
2, accounting for 0.00 to 19.27%, with an average of 2.87%, followed by nitrogen gas (N
2), accounting for 0.00 to 24%, averaging 1.51%. The drying coefficients of those samples vary from 0.67 to 0.99 (average 0.91), suggesting a variable maturity of their source rocks. The samples’ CH
4 content, C
2+ content and drying coefficient with burial depth exhibited a good correlation (Fig.3): whereas in the shallow and middle layers (< 3500 m), the CH
4 content and drying coefficient of the samples decreased with burial depth, and varied from typical dry gas to wet gas. It could be seen that as the depth increases, the CH
4 content and drying coefficient showed an increasing tendency in the deeper layers (3500–5500 m). Although the gas was still dominated by wet gas, it gradually changes to dry gas. As shown in Tab.4, in CO
2 gas reservoirs, the content ranged from 57.41% to 98.59%, averaging 82.8%, showing typical inorganic genesis (
Dai, 1993), while the other non-hydrocarbon gases were N
2, with contents ranging from 0.06% to 5.43%. The CH
4 content accounted for 0.44% to 42.59%, with an average of 12.21%.
4.2 Stable isotopes characteristics
The carbon isotopic values of all hydrocarbon gas sample exhibited a wide range: the carbon isotope of methane (
δ13C
1) ranged from −63.2‰ to −32.7‰, while the
δ13C
2 values fell between −38.3‰ and −16.8‰, the
δ13C
3 ranged from −34.9‰ to −7.64‰, and the
δ13C
4 varied between −32.1‰ and −14.9‰. The overall carbon isotopic composition of the alkanes follows the trend of
δ13C
1 <
δ13C
2 <
δ13C
3 <
δ13C
4, indicating a typical organic origin (
Liu et al., 2018), and in some cases (i.e., the Gudao, Shanjiasi, Chenjiazhuang, Caoqiao, Kenxi, Bonan-Gubei, Zhuangxi and Gaoqing-Huagou Gas Fields) varied with the
δ13C
3 >
δ13C
4 (Fig.4). By correlating the
δ13C
1 value with the sample burial depth (Fig.5), in the shallow strata (< 1500 m), the
δ13C
1 value of most sample was between −40‰ and −55‰, and only a few samples were less than −55‰. Reaching the middle strata (1500−3500 m), the value is between −45‰−60‰, which is significantly reduced. Reaching the deep depth (3500−5500 m), the value increased obviously, ranging from −30‰ to −55‰, implying that the genetic types of gases vary with their burial depth.
For the CO
2 gas reservoirs, the
δ13C
1 value was distributed from −54.39‰ to −42.51‰, and the carbon isotopes of the C
1 to C
4 alkane presented a positive-sequence distribution. The value of
δ13 ranged from −9.8‰ to −3.35‰. It has been established that the isotope characteristics of the associated inert gases helium (He) and argon (Ar) are important indicators for the study of CO
2 gas reservoirs (
Li et al., 2008;
Hu et al., 2009;
Ni et al., 2014;
Liu et al., 2016a). For the samples of this study, the
R/
Ra value (where
R is the
3He/
4He value of the sample and
Ra is the
3He/
4He value of the atmosphere) of CO
2 was generally greater than 2, and the
40Ar/
36Ar value was relatively high, ranging from 317 to 3000, and mostly higher than 1000.
5 Discussion
To facilitate the study, we divided gas samples into four groups: 1) Group I: comprising the samples from Dongfenggang and Dawangzhuang Gas Fields, located in the northwest part of the depression; 2) Group II: comprising the samples from Chengdao, Laohekou, Zhuangxi, Chengdong, Bonan-Gubei, Gudao, Gudong, Kenxi, Chenjiazhuang, Yanjia, and Yonganzhen Gas Fields, located in northeast part of the depression; 3) Group III: comprising the samples from Minfeng, Guangli, Lijin, Xianhe, Shengtuo, Shanjiasi, Liangjialou, Pingfangwang, Pingnan, Gaoqing, and Huagou Gas Fields, located in southeast part of the depression; 4) Group IV: comprising the samples from Qudi, Yuhuangmiao, Balipo, and Yangxin gas fields, located in the southwest part of the depression.
5.1 Genetic types of natural gas
Gas components and stable carbon isotope characteristics are often widely used as measurements for the classification of natural gas genetic groups, gas–gas and gas–source rock correlation, evaluation of maturity, and determination of gas migration direction and distance (
Dai, 1990,
1992;
Clayton, 1991;
Dai et al., 2005,
2012,
2014,
2016;
Ni et al., 2013;
Shi et al., 2019,
2020a,
2020b;
Hu et al., 2021b).
5.1.1 Genetic types of samples in north depression
For Group I, the CH
4 content was greater than 80%, and the corresponding drying coefficients were also higher than 0.85, while the
δ13C
1 values were lower than −40‰ (Table 2 ). These characteristics indicate that Group I gas is oil-associated gas. For Group II, as listed in Tab.1, the characteristics vary greatly. The minimum CH
4 content reached 56.9%, the minimum value of
δ13C
1 reached −63.2‰, and the maximum value of
δ13C
2 was −16.8‰. Some samples from Group II showed carbon isotope reversion, that is,
δ13C
3 >
δ13C
4 (Fig.4). It is thought that the reversed alkane carbon isotope trend can result from a variety of processes related to thermal effects, or mixing of the oil-associated gas and coal-derived gas, or the effects of secondary alteration (
Cai et al., 2001;
Dai et al., 2004;
Hao et al., 2008;
Liu et al., 2018). The
δ13C
3 value of Group II, whose gas samples originate from the Gudao, Shanjiasi, Chenjiazhuang, Caoqiao, Kenxi and Gaoqing-Huagou Gas Fields, was nearly 5‰ higher than the value of
δ13C
4; meanwhile, the propane content of Group II was lowest, and the drying coefficient was higher than 0.90, indicating that these gases are influenced by biological processes. However, for the Group II samples from the Bonan-Gubei and Zhuangxi Gas Fields, although the reversal phenomenon was similar to that in the above areas, the burial depth of these gases was sufficiently deep so as to exclude the effects of biodegradation, hence the real reason for partial reversal is related to the several sets of source rocks being comprised of matured E
s (Types I and II), Mesozoic coal-measure source rocks (Type III) with different OM types, which produced mixed gas.
Typically,
δ13C
2 values of −28‰ and
δ13C
3 values are generally considered as the criteria for identifying oil-associated gas from coal-derived gas (
Liu et al., 2019). Fig.6 displays the correlation between
δ13C
2 and
δ13C
3, Family II gas from Zhuangxi, Chengdao, Yanjia, Yonganzhen, Chenjiazhuang, and part of the Bonan Gas Fields (Wells BoS4, Yi11, XinYi12, Yi37, Yi52, Yi99, Yi170, Luo36-1, and BoS6-5) have values of
δ13C
2 less than −28‰ and a
δ13C
3 value less than −25‰, which demonstrates a typical pattern of thermogenic oil-associated gas. Furthermore, gases from the Gubei Gas Field, part of Bonan (Well Yi115, Yi121, BoG4, and BoG403) Gas Field, and well Gux1-11, reveal characteristics contrary to the above samples, in that they show typical patterns of thermogenic coal-derived gas.
These conclusions are further supported by the cross-plots of
δ13C
1,
δ13C
2, and
δ13C
3 (Fig.7) of natural gas and the C
1/C
2+3 and
δ13C
1 identification diagram of natural gas (Fig.8), which were proposed by Dai (
1990). The above samples confirm the carbon isotope threshold of oil-associated gas fall in zones II and IV in Fig.7, and those confirming the carbon isotope threshold of coal-derived gas fall in zones I and V in Fig.7. Group II gases from Wells Yi11-381, Yi11-G25, Yi17, and Yi171 had
δ13C
1 values lower than −55‰, and the values of
δ13C
2 were between −40‰ and −35‰, which falls in the VI zone, implying that these samples are biogas. However, these samples were wet gas, as their drying coefficient was lower than 0.90, indicating that the gases seem to be a mixture of thermogenic gas and biogas. Fig.8 provides further insight into the genetic types of the samples: some Group II gases from the Gudong Gas Filed had high C
1/C
2+3 ratios that plot outside the boundary of the D zone in Fig.8, while the rest fall well within the D and G zones, showing a typical pattern of oil-associated gas and coal-derived gas. For the samples that have high C
1/C
2+3 ratios, most over 100, the identification diagram established by
Bernard et al. (1978) can give the interpretation, as shown in Fig.9, in which the samples’ C
1/C
2+3 ratios are higher than thermogenic gases and fall within the secondary transformation zone. This is compared with typical oil-associated gas, in which the
δ13C
1 and
δ13C
2 values are 3–8‰ heavier than those of oil-associated gas, and show the distribution of
δ13C
1 <
δ13C
2 <
δ13C
3 >
δ13C
4. Thus, it can be stated that the high C
1/C
2+3 value is caused by secondary effects such as long-distance migration or biological action.
Oil-related gases include mature oil gas (oil-associated gas) and high-maturity oil gas. High-mature oil-related gas has high
δ13C
2 and
δ13C
3 values, and is commonly located near the boundary of the oil-associated gas zone in Fig.6. Using the differentiation chart (Fig.10) proposed by Wang et al. (
2009), it can be distinguished as to the genetic type of these samples, as shown in Fig.10, the
δ13C
3–
δ13C
2 value of Group II gas from the Laohekou, Bonan (Wells BoS6, Bos601, BoS3) and Zhuangxi (Well ZhuangG23) Gas Fields varies from 10 to 15, while the
δ13C
1 values are −45‰ to −35‰, plotting in zone III, which demonstrates a typical pattern of high-mature oil-related gas. In summary, the genetic types of Group II are very complex, except for non-inorganic gas, and other genetic types of natural gas are all involved.
5.1.2 Genetic types of samples in south depression
For Group III, the geochemical characteristics of the samples are also complicated, and the
δ13C
2 and
δ13C
3 values were generally lower than those of Group II, but some samples had higher
δ13C
3 values, up to −8.52‰. Biogas, whose
δ13C
1 value is typically less than −55‰, was only distributed in the Huagou Gas Field. The samples from Guangli, Xianhe, and Liangjialou Gas Fields had high contents of C
2+ components, and their drying coefficients were less than 0.95, which is typical of wet gas. Moreover, the values of
δ13C
1 were between −55‰ and −50‰, which shows the characteristics of typical oil-associated gases. Samples from the Panhe area all conformed to the oil-associated gas boundary in Fig.6 and fell into the zone II in Fig.7 and zone D in Fig.8. The samples from Caoqiao, Huagou, Gaoqing and Shanjiasi Gas Fields had shallow burial depth, and were significantly affected by secondary changes, as reflected in Fig.9. The figure shows that their C
1/C
2+3 and
δ13C
3 values were very high and, all fell in the secondary transformation zone; moreover, the distribution of
δ13C
3 >
δ13C
4 further confirm this conclusion. The samples from Shengtuo and Lijin, whose burial depth was more than 4000 m, had a slightly higher dryness coefficient than the oil-associated gas in Xianhe and other areas. Their
δ13C
2 and
δ13C
3 values were located near the boundary of the oil-associated gas zone, and the
δ13C
3–
δ13C
2 values and
δ13C
1 values fell within the range of high-mature oil-related gas zone in Fig.10, which shows that they are typical high-maturity oil-associated gases. Although the natural gas from Minfeng was ever distinguished as kerogen cracking gas by a few samples (
Luo et al., 2008;
Hu et al., 2009), according to the judging chart (Fig.11) proposed by
Prinzhofer and Huc (1995), the variation trend of ln(C
2/C
3) and
δ13C
3–
δ13C
2 of the gas samples is more consistent with that of oil-cracked gas. In addition, the
δ13C
1 value was lower than that of typical kerogen pyrolysis gas, so it can be determined that the gas samples in the Minfeng area are mainly oil-cracked gas, which is consistent with the conclusions of Li et al. (
2010) and Yang et al. (
2014).
In addition to the above hydrocarbon reservoirs, Group III also contain CO
2 gas samples from the Huagou, Pingfangwang, Pingnan Gas Fields. Further, it is known that the genetic type of CO
2 can be classified as inorganic or organic. Considering this, these samples had high
δ13 values, most over −8‰, ranging from −9.8‰ to −3.35‰, indicating inorganic characteristics. The cross-plot of CO
2 content vs.
δ13 also confirms that these samples are inorganic gas, as shown in Fig.12 (
Dai et al., 1992), in which all samples fall into zone II. Next, we consider that inorganic gas can be further divided into mantle-derived gas and crust-derived gas, which can provide meaningful information about the accumulation model of CO
2 reservoirs (
Dai et al., 1992;
Li et al., 2004;
He et al., 2005a,
2005b). According to the relationship between
δ13 and
R/
Ra,
He et al. (
2005a,
2005b) proposed an identification chart for inorganic CO
2 (Fig.13). For this study, as shown in Fig.13, the
R/
Ra values of most samples were higher than 2 and fell in the mantle-derived zone. In addition, their
40Ar/
36Ar values were also relatively high, ranging from 317 to 3000, most of which were greater than 1000 (Tab.3), revealing the characteristics of mantle-derived Ar.
For Group IV, as shown in Tab.2 and Tab.3, the samples from the Yangxin Gas Field had low δ13C1 value, ranging from −60.9‰ to −56.5‰ (av. −59.3‰), and a high drying coefficient (all over 0.99), showing a typical pattern of biogas. As shown in Fig.7, the cross-plots of δ13C1 versus δ13C2 and δ13C3 demonstrated that the samples from the Qudi and Yuhuangmiao Gas Fields fell into various regions. The former fells into zone I, showing coal-derived gas characteristics, while the latter fell into zone II, showing oil-associated gas characteristics. In addition, the CO2 gas samples from the Balipo Gas Field had high R/Ra values, which fell into the volcanic mantle-derived inorganic CO2 zone in Fig.13.
It has been established that, based on the characteristics of natural gas components and carbon isotopes, natural gas samples can typically be divided into five types: oil-associated gas, high-mature oil-related gas, coal-derived gas, biogas, and secondarily generated gas, among which, the oil-associated gas is predominant. In this study, we find specifically that, Group I is the only oil-associated gas, while Group II has all genetic types, which are dominated by oil-associated gases. The gas samples from Group IV are mainly biogas, while Group III is largely oil-associated gas and high-mature oil-associated gas.
5.1.3 Distribution of gas genetic types
As shown for the distribution of natural gas genetic types in the Jiyang Depression (Fig.14), the law of type changes from north to south and has two changing trends as follows. First, in the east, when reaching the Chenjiazhuang uplift, the complex types of gas changes to only oil-associated gas, and then the genetic types become complicated again. In the west, the changing trend of the genetic type of natural gas becomes gradually, increasingly complicated from north to south.
Furthermore, combining the genetic type of the gas reservoir with its burial depth (Fig.15) reveals that the vertical distribution of natural gas in the Jiyang Depression has no obvious progressive change. Oil-associated gases are distributed in all layers, while high-mature oil-related gas is mainly distributed in deep layers, coal-derived gas is mainly distributed in the middle and deep layers, and the inorganic CO2 gas reservoirs are mainly distributed in the middle and shallow layers.
According to the
Ro calculation formulas in general for oil-associated gas:
δ13C
1 = 15.8 lg
Ro−42.2 and
δ13C
1 = 21.72 lg
Ro−43.31 and for coal-derived gas:
δ13C
1 = 14.12·lg
Ro−34.39 (
Dai, 1993;
Wang, 2008;
Li et al., 2015), it is implied that the values of
Ro do not increase with the increase in burial depth, which is confirmed by the findings of this study (Fig.16). However, the samples also show that, in a certain depth range, incremental changes still occur. This indicates that in this study area, the vertical distribution of genetic types of natural gas is not only controlled by the maturity of the source rock, but also by other conditions such as tectonic evolution (
Wang et al., 2017;
Wu et al., 2017)
5.2 Main controlling factors for the occurrence of natural gas
5.2.1 Distribution of gas genetic types controlled by source rock conditions
Previous studies have demonstrated that the source rocks exert an important control on hydrocarbon accumulation and distribution (
Guo et al., 2014;
Liu et al., 2016b;
Tian et al., 2018;
Li et al., 2019b;
Zhang et al., 2020). The spatiotemporal changes of gas reservoirs shown in Fig.14 and Fig.15 demonstrate that gas reservoirs in the Jiyang Depression are closely related to the qualities and thermal evolution of the source rocks. The source rocks in Jiyang Depression mainly include Paleogene E
s, E
k, and Upper Paleozoic Carboniferous–Permian strata.
The northwest area, where Group I gas is located, extending to the Yihezhuang Uplift in the east and the Chenjiazhuang Uplift in the south, belongs to lacustrine deposits. The main source rocks are Es1, Es3 and Es4, and the main hydrocarbon-bearing strata are Paleogene strata. Although the Lower Paleozoic carbonate strata can also generate hydrocarbons, they are not the main source rocks. The results of this study show that the TOC value of Es4 is between 1.0% and 8.0%, Es3 is between 2.5% and 12.1%, and Es1 is approximately 5.5%, implying that they are all excellent source rocks, and lay a material foundation for the generation of natural gas. The burial depths of Es4 and Es3 are greater than 3000 m, and the value of Ro is over 0.5%, indicating that they have entered the peak of hydrocarbon generation, with both Type I and Type II1 kerogen, mainly producing oil and oil-associated gas. The Es1 Formation is distributed throughout the entire area, but its burial depth is 1600–2800 m, the value of Ro < 0.5‰, in the immature stage, so it is not the main gas production unit. Carboniferous to Permian coal-measure strata is also distributed in this area, but the average thickness of the coal seam is very small, and its Ro value is low (Fig.17), indicating that these strata have just entered the early mature stage without a large amount of hydrocarbon generation. Although the Mesozoic strata in the Jiyang Depression have a secondary hydrocarbon generation process during the Himalayan period, they were not found in this area. Therefore, only oil-associated gas reservoirs are found in this area.
The northeast area, where Group II gas is distributed and which is bordered by the Yihezhuang Uplift in the west and the Chenjiazhuang Uplift in the south, includes many oil-rich subsags, for example, Bonan subsag. Compared with the northwest area, except for E
s1, E
s3, and E
s4, the Carboniferous–Permian coal measure strata are also the main source rock. The E
s41 in this area has a high TOC value (ranging 0.77% to 15.18%), and a burial depth of up to 5200 m, corresponding to a high value of Ro (between 0.8% and 1.56 ‰), implying that it is in the high-mature stage and in the oil generation-condensate stage. It produces high-mature oil-related gases. The characteristics of E
s3 in this area are similar to those of E
s3 in the northwest, the corresponding TOC values are between 0.55% and 12.44%, and the main OM is Type I and Type II
1. Its
Ro value ranges from 0.64% to 1.12%, inferring that it is a mature source rock that mainly generates oil-associated gas. The coal-derived gas is mainly produced by the Carboniferous–Permian coal measure strata in the Bonan-Gubei area, which is mainly distributed in the buried hill area (Fig.17). The thickness of the coal seam is approximately 30 m, and the
Ro value is between 0.94% and 1.77%, indicating that it is a mature source rock, which is conducive to a large amount of gas generation. Moreover, it is thought to have two hydrocarbon generation processes, the first hydrocarbon generation period being at the end of the Middle Triassic, and the second being in the Yanshanian–Hishanian Period (
Wang et al., 2010;
Zhu et al., 2010). Typically, since the first generation of hydrocarbons is difficult to preserve through tectonic evolution, the gas generated by the second generation of hydrocarbons can form gas reservoirs, and has higher research value. The biogas in this area was mainly generated from E
s1. The corresponding TOC value was between 0.9% and 7.0%, and the OM was Type I, the buried depth was less than 3000 m, and the
Ro was between 0.43% and 0.60%, implying that it is in the immature stage, generating biogas reservoirs. In addition, the shallow biological action in this area is relatively strong, which is another main reason for the formation of the biogas reservoirs.
In the southeast area, where Group III gas is distributed, the main source rocks in this area are E
s4 and E
s3. In recent years, studies have determined that the second member of the E
k Formation has a certain hydrocarbon expulsion potential, but it is not the main source rock (
Ren et al., 2006;
Chen et al., 2011). The burial depths of E
s3 and E
s41 are relatively shallow (between 2000 and 3000 m) corresponding to the
Ro value ranging from 0.58% to 1.03% and their OM is mainly Type I and Type II
1 mixed with a small amount of kerogen II
2, indicating that they mainly produce oil-associated gas reservoirs. The oil-cracked gas mainly comes from the lower subsection of E
s4, which belongs to salt-lake sediment, and its burial depth is more than 4200 m in the Minfeng area, with the ground temperature exceeding 160°C, showing that it has the conditions for cracking crude oil into gas. The Carboniferous–Permian strata are also distributed in this area, and the thickness of the coal measure is approximately 20 m, which has the potential for hydrocarbon generation. However, the hydrocarbon generation process was mainly in the Indosinian movement, and the natural gas reservoirs that were generated at this time were difficult to preserve after structural changes. After the Yanshanian and Himalayan Periods, the ground temperature increased with increasing the burial depth, but the secondary hydrocarbon generation process was not found in this area (
Fan et al., 2008;
Wang et al., 2010). Owing to the influence of strong volcanic movement, it is not easy to generate large-scale gas reservoirs, so no coal-derived gas reservoirs are distributed in this area.
The southernmost part of this area is the Huagou and Gaoqing Gas Fields, which are connected to the Qingchengshan Uplift in the north. This area contains both CO2 gas reservoirs and hydrocarbon gas reservoirs. The CO2 gas reservoir is of volcanic mantle origin, which is related to the development of multiple magma eruptions. The main source rocks of the hydrocarbon gas reservoirs in this area are Es1 and Es3. The burial depth of Es3 is between 1800 and 3500 m and the corresponding Ro value is between 0.36% and 0.8%, indicating that it mainly produces oil-associated gas. Meanwhile, the Ro value of Es1 is < 0.5%, which is in the immature stage and generates biogas reservoirs.
The southwest area, where Group IV gas is distributed and whose source rocks are also Es4 and Es3, reaches to the Wudi Uplift in the north and Linfanjia Uplift in the east. The TOC value of them is about 3%, the burial depth is approximately 3000 m, and the Ro value is greater than 0.65%, indicating the early mature stage, while the OM is mainly Type I and Type II1, generating oil-associated gas reservoirs. The main characteristic of this area is the coal-derived gas reservoirs in the Qudi area. In the buried hill zone of the Qudi area, the Upper Paleozoic Carboniferous–Permian coal-measure strata are rich in coal resources (Fig.17), the thickness of the coal seam is large, and the OM type is Type III, which is conducive to gas generation. According to the history of hydrocarbon generation, the Carboniferous–Permian strata experienced two peaks of hydrocarbon generation after the deposition of the Early Triassic and the Jurassic to the Cretaceous. In addition, the Qudi buried hill area is close to the secondary gas generation center with sufficient gas injection. Reservoir types include sandstone, carbonate rock and igneous rock, all of which perform good storage conditions, so coal-derived gas reservoirs can be generated. The Yangxin subsag is located in the northeast part of this area. Biogas reservoirs are also distributed in this area, primarily due to the relatively well-developed bio-limestone in the Es1 formation, and its shallow burial depth (approximately 1750 m), low thermal evolution (the Ro value is less than 0.4%).
In general, by summarizing the source rock characteristics of each part (Tab.2), it is demonstrated that the lateral distribution of various genetic types of natural gas in the Jiyang Depression is controlled by the characteristics and distribution of the Es1, Es3, Es4 and Carboniferous–Permian coal-measure strata. To further explore and develop certain types of natural gas in the Jiyang Depression, future research should look for potential zones based on the distribution of these source rocks.
5.2.2 Distribution of gas genetic types controlled by tectonic conditions
Tectonic conditions in general refer to various transport systems and tectonic movements in various periods that affect the formation of natural gas reservoirs (
Zhang, 2012;
Li et al., 2021b). They also play a decisive role in the distribution of various types of gas reservoirs, especially with regard to the vertical distribution.
In the Jiyang Depression, it has been established that tectonic conditions control the thickness and distribution of source rocks, thereby affecting the distribution of gas reservoirs (
Wu, 2018;
Chen, 2019). According to the occurrence of faults, fault regions can be divided into three groups: NE-, NW-, and near-E trending. The number and location of faults are well correlated with the thickness of the source rocks (Fig.18). As shown in Fig.18, source rocks are mainly distributed near active concave-controlling faults, and faults clearly control the extent of source rocks, that is, E
s3 extends along the NE direction and E
s4 is distributed in two directions along the NE and NW directions.
It is known that faults’ sealing ability controls the process and the model of gas accumulation, hence for the study area it can be said as follow. 1) The biogas that was distributed in the mid-deep layer in the Bonan area was formed when the shallow biogas was transported to the middle layer under the favorable transportation condition of the Guxi fault. Some samples in this area are mixed sources that are inseparable from the development of multistage unconformities and good lateral migration channels. 2) The distribution of shallow gas reservoirs is directly related to tectonic conditions. Most shallow gas reservoirs are developed in the Neogene Minghuazhen Formation (N
m) and Guantao Formation (N
g). One of the factors for this is that those areas have good migration conditions in which many faults can communicate with the lower source rocks E
s. The convex shallow gas accumulation model in the Gudao-Gudong area is a coincident example (
Cheng et al., 1995). Another factor is that following the late Tertiary, the activity rate of faults had been very low, so the damage to the caprock was slight, which is conducive to the migration and accumulation of natural gas. 3) The distribution of coal-derived gas reservoirs is related to the secondary hydrocarbon generation of the Carboniferous–Permian strata, in which the evolution of tectonic activity determines its process. At the end of the Indosinian period, the first hydrocarbon generation stopped owing to the uplift of the depression. Reaching the rapid subsidence stage of the Yanshanian–Himalayan Period, the burial depth of the Carboniferous–Permian strata continued to increase, and the strata was affected by thermal anomalies, so the second hydrocarbon generation was initiated (
Wang et al., 2010;
Zhu et al., 2010). Because it is limited by the distribution of subsidence centers, secondary hydrocarbon generation was generated only in the northeast and southwest areas, and coal-derived gas reservoirs were only distributed in these two areas.
The finding shows that the distribution of inorganic CO
2 gas reservoirs in the study area is mainly affected by tectonic conditions. For example, in the southern area, which has undergone been multiple periods of magmatic activity since the Mesozoic and Cenozoic, the Balibo, Gaoqing, and Huagou Gas Fields are comprised of tertiary alkaline basalts located in this area. These alkaline basalts are facies conducive to the release of CO
2, and they provide a good material basis for the accumulation of CO
2 gas reservoirs (Shen et al.,
2007;
Hu et al., 2009;
Wang et al., 2018). In general, it is understood that mantle-derived CO
2 gas generation requires deep connections to transport CO
2 in the mantle. The Tanlu Fault is the only fault that connects to the upper mantle lithosphere in the Jiyang Depression, thus mantle-derived CO
2 gas can be transported through it. In this area, only the Yangxin and Shicun Faults in the south connect the Tanlu Fault, and the Linshang, Xiakou, and Gaoqing-Pingnan Faults in the southern area have large cutting depths, which are also favorable channels for the transport of deep CO
2. Therefore, in the southern region, CO
2 gas reservoirs are mostly distributed at the intersection of the NE-, E-, and NW-trending faults (Fig.19). In summary, the CO
2 gas reservoirs which are further favorable exploration targets are those representing structural and lithological traps near the fault zone in the southern area.
6 Conclusions
The hydrocarbon gas in the Jiyang Depression can be divided into four groups. Group I gas distributed in the northwest is only oil-associated gas. Group II gas distributed in the northeast area, is dominated by oil-associated gas, and involves biogas, coal-derived gas, and high-mature oil-relate gas. Group III, distribute in southeast area, has all genetic types gas and that also dominated by oil-associated gas, samples from the shallow layer have high δ13C3 values, showing the distribution of δ13C1 < δ13C2 < δ13C3 > δ13C4. And Group III also contains volcanic mantle-derived carbon dioxide gas. Group IV gas distributed in the southwest area, is dominated by biogas and involves coal-derived gas and oil-associated gas. In short, natural gas samples in the eastern and southern depressions have more genetic types. However, in the vertical direction, the distribution of natural gas does not have an obvious evolutionary law.
There are two main factors controlling factors of the spatiotemporal changes of gas are in two aspects. First, the characteristics of the Shahejie Formation and Carboniferous Permian source rocks are various in different area. The occurrence of oil-associated gas reservoirs is mainly controlled by the maturity of Es1, Es3, and Es4 and its burial depth. Second, the tectonic conditions, mainly the sealing properties of various faults, affected the lateral and vertical migration of hydrocarbon gases and the upwelling of mantle-derived CO2. Moreover, the secondary hydrocarbon generation process of Carboniferous Permian is controlled by the uplift of strata, which has affected affecting the lateral distribution of coal-derived gas reservoirs. In general, the main controlling factor analysis of the spatiotemporal changes of gas reservoirs revealed that the synergy of characteristics, thermal evolution of the Shahejie Formation and Carboniferous–Permian source rocks, and sealing properties of various faults are jointly responsible for determining the spatiotemporal changes of gas reservoirs.