1 Introduction
Shale gas is a crucial hydrocarbon resource, which is distributed worldwide and has enormous potential for development (
Rimmer et al., 2004;
Ding et al., 2013;
Tang et al., 2014;
Khaled et al., 2022). Chinese shale is characterized by its wide distribution, high organic matter content, and large resources. According to their depositional environment, they can be classified into three kinds of types: continental, marine-continental transitional, and marine shales, among which marine-continental transitional shale accounts for a higher proportion and has better development prospects (
Wei et al., 2020;
Zhang et al., 2021). The Carboniferous to Permian was a vital period in the transformation of the depositional environment of China into a marine-continental transitional facies. Transitional shale is widely deposited in the Ordos, Tarim, and Junggar Basins and other regions. The marine-continental transitional shale in the Ordos Basin is developed in the Benxi Formation, Taiyuan Formation, and Shanxi Formation (
Song et al., 2016;
Yang et al., 2019). The Benxi Formation, as a typical organic-rich shale of the marine-continental transitional facies, has enormous potential for development.
Over the last few decades, the geochemical characteristics, depositional environments, organic matter enrichment modes, deposition modes, and gas-in-place content of organic matter-rich shale have been extensively studied (
Li et al., 2019a;
Chen et al., 2020;
Li et al., 2021;
Yu et al., 2022). In general, organic matter enrichment modes can be divided into two types. The “production mode” emphasizes the influence of factors such as palaeoclimate, palaeoproductivity, and terrestrial clastic input on organic matter enrichment, while the “preservation mode” emphasizes the control of organic matter by factors such as redox conditions, palaeosalinity, and sedimentation rate (
Ma et al., 2019;
Awan et al., 2020;
Khaled et al., 2022;
Wu et al., 2022a;
Zhou et al., 2022). In terms of the enrichment of shale organic matter by palaeosalinity, it has been widely concluded that with the increase of palaeosalinity, organic matter increases first and then decreases. Shale developed in semi-saline-alkaline environment is the most ideal target for exploration and development. (
He et al., 2018;
Hu et al., 2018). Commonly, geochemical parameters are used to reconstruct the palaeoenvironment, while CIA, Sr/Cu, and C
value can be widely used in evaluating palaeoclimatic conditions quantificationally. P/Ti and Mo are used as geochemical indicators of palaeoproductivity, and Fe/Ti and (Fe + Mn)/Ti are used to restore volcanic hydrothermal fluid (
Nesbitt and Young, 1982;
Yamamoto, 1987;
Bernárdez et al., 2008;
Stock et al., 2017;
Ma et al., 2019;
Wang et al., 2020). As a result of previous studies, Aluminum (Al), Titanium (Ti), and Thorium (Th), which are transported into the water column and are less susceptible to diagenesis and weathering, can often indicate the input of terrestrial clastic, while (La/Yb)
N can be used to evaluate the rate of sedimentation (
Tyson, 2001;
Ding et al., 2015;
Doner et al., 2019). Rare earth and trace elements vary in enrichment intensity under different redox conditions, and trends in the ratios of Ni/Co, U/Th, and V/(V + Ni) are commonly used to characterize changes in the redox environment of the deposition site. Palaeosalinity is usually used to identify the extent and duration of water intrusion, while the ratios of Sr/Ba and Th/U, which are composed of trace elements sensitive to water column salinity, are common indicators of palaeosalinity (
Jones and Manning, 1994;
Wignall and Twitchett, 1996;
Stüben et al., 2002;
Rimmer, 2004;
Zhang et al., 2008;
Jenkyns et al., 2017). The organic matter enrichment modes of marine and continental shales have been extensively studied. The continental shale is mainly developed in lake basins, where the paleoclimate was cold and dry or warm and humid during the deposition period, and the water reducibility and palaeosalinity were low. Organic matter enrichment is mainly influenced by redox conditions and palaeoproductivity, with palaeoproductivity being the main controlling factor (
Yu et al., 2022). Marine shale is mainly deposited on continental margins, where the climate is warm and wet and the water column is dominated by reducing conditions and high palaeosalinity. There is also possible submarine hydrothermal influences, where organic matter enrichment is influenced by multi-factor and is usually dominated by a preservation model (
Zhang et al., 2018a;
Li et al., 2019a;
Li et al., 2020). The transitional sedimentary environment is rich in shale, but its organic matter enrichment pattern is lacking, mainly concentrated in the Shanxi–Taiyuan Formation in the Ordos Basin, the Longtan Formation in the Sichuan Basin, and the Cuizhuang Formation in the Yuncheng Basin (
Liu et al., 2018;
Li et al., 2019c;
Chen et al., 2020;
Wei et al., 2020;
Khaled et al., 2022;
Yang et al., 2022). In terms of gas-in-place content in shales, previous work has developed a gas-in-place content calculation model based on a carbon isotope fractionation model and evaluated gas-adsorbed ratio, providing a powerful tool for evaluating key parameters of unconventional gas resources (
Li et al., 2022). The above studies have systematically investigated the enrichment mode of organic matter shales, but the geochemical characteristics of the marine-continental transitional shale of the Benxi Formation remain to be clarified. Their palaeoenvironmental conditions and the mechanisms influencing organic matter enrichment are yet to be revealed, while the comparison of the differences in organic matter enrichment mechanisms between different sedimentary facies is open to debate.
This research systematically sampled shales of the Upper Carboniferous Benxi Formation in the east-central part of the Ordos Basin. Organic geochemical tests and elemental geochemical analyses were conducted to recover the palaeoenvironment of the stratigraphic deposition period in this area, explore the influences on the palaeo-depositional environment that affected organic matter enrichment (for instance, palaeoclimate, water redox conditions, palaeosalinity, palaeoproductivity, volcanic hydrothermal fluid, and sedimentation rate), and the enrichment mechanism of organic matter was revealed. The differences in organic matter enrichment mechanisms of shales under different sedimentary facies are also compared and analyzed. These results provide a reference for shale gas exploration and reservoir selection.
2 Geological setting
The Early Palaeozoic Middle Ordovician Ordos Basin was in the stage of marginal rifting and intra-land depression, and the main body of the basin is a marine sedimentary environment with carbonate rocks deposited on metamorphic rocks at the base of the basin during the Taurus and Metamorphic periods. In the Early-Late Palaeozoic, the basin entered a stage of circumferential collisional orogeny, and the Ordos basin was uplifted and subjected to long-term exfoliation. In the end of Late Paleozoic, circumferential rifting occurred and the basin receded, and the sedimentary background was a marine and land interactive sedimentary system, when the North China Plate on which the Ordos Basin is located was in the northern hemisphere, closer to the equator (Fig.1(a)). The Mesozoic basin underwent intra-land depression. The basin margins were uplifted and lifted off as a whole and it had a terrestrial sedimentary environment. The Cenozoic basin perimeter was gradually fractured to form the present-day tectonic pattern (
Yang et al., 2005;
Ding et al., 2013;
Li et al., 2019b). The basin consists of six primary tectonic units: the Yimeng uplift, the Weibei uplift, the Tianhuan sag, the Yishaan slope, the Jinxi folding belt, and the western edge thrust belt (
Tang et al., 2014). The Ordos Basin underwent transformation through marine to marine-continental transitional to continental-phase sedimentary environments during the Palaeozoic. The Carboniferous-Permian stratigraphy of the marine-continental transitional facies has developed, from bottom to top, in five Formations: Benxi, Taiyuan, Shanxi, Shihezi, and Shiqianfeng, among which the Benxi Formation is the marine-continental transitional facies with deltaic and barrier coastal sedimentation, and the research area is situated at the middle-eastern part of the Yishaan slope, mainly developing gray and gray-brown shales and light gray sandstones (Fig.1(b)−Fig.1(d)) (
Peters and Cassa, 1994;
Li et al., 2019b;
Guo, 2020).
3 Sample preparation and testing
3.1 Sample preparation
The 16 shale core samples in this study were obtained from 10 wells including M35, M109, and Q20 of the Benxi Formation in the central-eastern part of the Ordos Basin (Fig.1(d)). The rock samples were dried after several washes with deionised water in an ultrasonic environment, then ground to 200 mesh in an agate mortar and subjected to geochemical testing and analysis, including total organic carbon (TOC) and major and trace element tests.
3.2 Test procedure
TOC was tested using a Vario TOC Total Organic Carbon analyzer from Element, Germany, using the infrared absorption method. A 200 mesh powdered sample of 20mg was weighed and subjected to a stream of oxygen at 950°C to oxidise the solid carbon to carbon dioxide. The resulting gas is tested through a carbon dioxide infrared detection cell and converted to the relative content of organic carbon. The test method is performed according to GB/T 19145-2003.
The experimental procedure of the main element tests was referenced from GB/T 14506.28-2010 (
Yang et al., 2020). The 200 mesh powdered shale samples were placed in an open muffle furnace at 950°C and dissolved in a mixture of HClO
4 and HF. X-ray fluorescence (XRF) measurements were carried out using a combination of AA-6800 atomic absorption spectroscopy and UV-2600 UV-Vis spectrophotometer.
For trace and rare earth elements, 200 mesh powdered samples were dissolved in a mixture of HNO3, HF, and HClO4 in a microwave oven and extracted with aqua regia after 5 h of digestion at 230°C to make a final 50 mL solution. The prepared liquid sample was measured by a Perkin Elmer SciexElan 6000 inductively coupled plasma mass spectrometer (ICP-MS). The test accuracy was strictly controlled and the analytical uncertainty was kept within 1% (the test procedure was following GB/T 14506.30-2010 standard).
4 Results
4.1 Organic matter abundance
The TOC content of shale samples in the study area ranges from 0.79 wt% to 9.83 wt%, the mean value was 3.37 wt% (Tab.1). The mean TOC content of the shale samples was higher than 2 wt%, and a TOC between 2 wt% and 4 wt% is considered a good hydrocarbon source rock abundance criterion based on the previously proposed criteria for evaluating organic matter abundance in hydrocarbon source rocks (
Peters and Cassa, 1994). Therefore, the Benxi Formation shale has excellent hydrocarbon source rock potential.
4.2 Main element characteristics
The main elements of the shale of the Benxi Formation are SiO2 as the main component, with a content of 51.36−92.30 wt% (average 61.75 wt%), then Al2O3, Fe2O3, and K2O with average contents of 27.64 wt%, 6.30 wt%, and 1.90 wt%, respectively. The remaining main element concentrations are less than 1 wt% (Tab.2). Relative to the Upper Continental Crust (UCC) mean counterparts, the Na2O, MgO, P2O5, K2O, CaO, and MnO contents of the Benxi Formation shale show significant deficits, while Al2O3, TiO2, and Fe2O3 are relatively enriched (Fig.2(a)).
4.3 Trace element characteristics
Trace elements in sedimentary rocks are often used to reconstruct their palaeoenvironment during the time of deposition. This study mainly considered 14 trace elements, For instance, V, Cr, and Co (Tab.3). The elemental enrichment factor (EF) is used to quantitatively describe the enrichment or deficit of trace elements in shales. It is calculated as EF = (
Xsample/Al
sample)/(
XUCC/Al
UCC) (
Wedepohl and Earth, 1971;
McLennan, 2001). Compared to UCC standards, the Benxi Formation shale is relatively enriched in V, Cr, Co, Ni, Mo, U, and Y (EF > 1), while Zn, Sr, Zr, Ba, Hf, Ta, Th, and U are deficient (EF < 1; Tab.4).
4.4 Rare earth element characteristics
The rare earth elements (REE) results for Benxi Formation Shale and corresponding standard values, such as Chondrite, are shown in Tab.5 (
Taylor and McLennan, 1995). Total rare earth elements (ΣREE) variability in the Benxi Formation shale ranges from 50.48 to 621.49 ppm (mean 274.23 ppm), significantly greater than that of the North American Shale Complex (NASC; 173.21 ppm;
Haskin et al., 1968) and the Upper Continental Crust (UCC; 146.37 ppm;
McLennan, 2001), a little more than the Post-Archean Australian Shale (PAAS; 204.12 ppm;
Taylor and McLennan, 1985)). The light rare earth elements (ΣLREE) and heavy rare earth elements (ΣHREE) averaged 247.66 ppm and 26.56 ppm, respectively, and the ratio of light rare earth elements / heavy rare earth elements (ΣLREE/ΣHREE) ranges from 5.63 to 19.52 (average 10.04; Tab.5), which is higher than those of North American Shale (7.50), indicating that the LREEs in the shale of Benxi Formation is more enriched than the HREEs. The values of (La/Yb)
n, (La/Sm)
n, and (Gd/Yb)
n can be used to indicate the enrichment or deficiency of light and heavy REEs. High values of (La/Yb)
n > 1 reflect the enrichment of heavy rare earths, while low values of (La/Yb)
n < 1 are a deficit, while the two values of (La/Sm)
n and (Gd/Yb)
n reflect the fractionation of light and heavy rare earths, with higher values indicating a higher degree of fractionation (
Hoyle et al., 1984). The (La/Yb)
n, (La/Sm)
n, and (Gd/Yb)
n of the Benxi Formation shale ranges from 6.41 to 36.54 (average 13.70), 2.45−8.87 (average 5.06), and 0.73−2.98 (average 1.83), respectively, all indicating relatively high fractionation of light and heavy rare earths (Tab.5). After normalized by chondrite, the distribution model of REEs in the shale samples of the study area showed a right-sloping trend for LREEs and a flattening trend for HREEs (Fig.2(c)).
5 Discussion
5.1 Palaeo-depositional environment of the Benxi Formation
5.1.1 Depositional site, palaeodepth, and material source
The aqueous environment of the deposition site is important for the preservation and production of organic matter in shale. The Al
2O
3/(Al
2O
3 + TFe
2O
3) values of oceanic ridges, ocean basins, and land margins are < 0.4, 0.4−0.7, and 0.7−0.9, respectively (
Roser and Korsch, 1986). The Al
2O
3/(Al
2O
3 + TFe
2O
3) of the Benxi Formation ranges from 0.41 to 0.97 (mean 0.80), indicating that the Benxi Formation was deposited mainly on the terrestrial margin (Tab.6) (
Roser and Korsch, 1986;
Khaled et al., 2022). The La-Th-Sc diagram and ratio of K
2O/Na
2O to SiO
2 can be used to create a discriminant map of depositional sites. The La-Th-Sc diagram shows that most of the samples are located in or near the active continental arc and passive-active continental arc margin regions (Fig.3(b)). The K
2O/Na
2O-to-SiO
2 ratio discrimination plot shows that the vast majority of the samples are located in passive continental environments (Fig.3(c)). Both results are consistent with the Al
2O
3/(Al
2O
3 + TFe
2O
3) ratio results, indicating that the Benxi Formation Shale was deposited on a passive continental margin. In this paper, the shale data from the Yichuan area and the Benxi Formation of the Gaoqiao area in the southern Ordos Basin were selected for comparison with the present study. Fig.3(b) shows that the depositional position of the Yichuan area is mainly located in the continental arc, while the main depositional position of the Gaoqiao area is consistent with the present test results, with most samples falling into or adjacent to the active continental arc and passive active continental arc margin areas (
Hu et al., 2022).
The stability of manganese (Mn), Ferrum (Fe), and titanium (Ti) elements in sediment transport varies. In general, smaller Mn/Fe and Mn/Ti ratios reflect shallow palaeowater depths (
Wang et al., 2020). Mn/Fe and Mn/Ti ratios for the shale samples ranged from 0.002 to 0.189 (mean 0.017) and 0.007−8.853 (mean 0.631), respectively, which are low and indicate relatively shallow palaeowater depths during this period at the continental margin and a marine-continental transitional depositional environment (Tab.6).
The clastic composition reflects the area of origin, stripping, and transport and controls the hydrocarbon potential of the shale, and the cross-plot of ∑REE and La/Yb can be used as a common indicator of sediment origin (
Floyd and Leveridge, 1987). The majority of shale samples from the study area fall into granite and sedimentary rock areas, and most of the Benxi Formation samples in the Yichuan and Gaoqiao areas were also obtained from granite and sedimentary rock areas (Fig.3(a)). This suggests that the source rocks of the shale of the Benxi Formation are granites and sedimentary rocks of the upper crust, consistent with the source rocks of the Yichuan and Gaoqiao areas in the southern part of the Benxi Formation studied by previous authors. The Y/Ho ratio is often used to reduce the origin of detrital material in the sediment. These values are about 28 for terrestrial clastic rocks and 44−74 for marine chemical sediments. The Y/Ho ratios of the present study range from 23.55 to 32.25 (mean 27.20; Tab.6), indicating that the Benxi Formation sediments are dominated by terrestrial clastic rocks (
Xu et al., 2014). Metallic elements such as Titanium (Ti) and Aluminum (Al) are often used as indicators to assess the extent of surface clastic flux (
Hatch and Leventhal, 1992;
Zhao et al., 2021). The TiO
2 and Al
2O
3 contents of the Benxi Formation shale range from 0.07% to 1.22% and 2.97%−39.76%, averages of 0.71% and 27.64% (Tab.2), respectively, which are greater than the UCC standard values (0.64% and 15.40%, respectively) and close to the TiO
2 (1.00) and Al
2O
3 (18.90) contents of the PAAS, indicating that the Benxi Formation shale has a high clastic input. TiO
2 has a high positive correlation with Al
2O
3 (
R2 = 0.75; Fig.3(d)), indicating that the titanium (Ti) in Benxi Formation shale mainly comes from terrestrial clastic material, further implying that the deposition of these shales was accompanied by strong terrestrial clastic input.
5.1.2 Palaeoclimate
Previous studies have shown that palaeoclimate influences to some extent the intensity of chemical weathering and debris influx, which in turn influences the mineralogy and chemical composition of the sediments.
Nesbitt and Young (1982) put forward an erratics index (CIA = Al
2O
3/(A1
2O
3 + CaO* + Na
2O + K
2O) × 100) for the evaluation of palaeoclimatic changes. CaO in silicates as CaO* and is corrected by the expression CaO* = CaO − (10/3 × P
2O
5). The result is taken as the difference between CaO* and Na
2O content the smaller value (
Nesbitt and Young, 1982). 50 < CIA < 70, 70 < CIA < 85, and 85 < CIA < 100 indicate cold and dry, warm and humid, and hot and humid climates, respectively, and the degree of chemical weathering gradually increases with the CIA value (
Nesbitt and Young, 1982;
Li et al., 2019c). The Index of Constituent Variability (ICV) is commonly used to represent sediment recycling and is a prerequisite for the CIA. ICV = (Fe
2O
3 + K
2O + Na
2O + CaO* + MgO + MnO + TiO
2)/Al
2O
3, with ICV > 1.0 indicating that deposition occurred during a tectonically active period and that the sediments contain small amounts of clay minerals. Meanwhile ICV < 1.0 indicate that debris was redeposited or deposited under strong weathering conditions, at which point the CIA is indicative of the degree of chemical weathering (
Fedo et al., 1995;
Cullers and Podkovyrov, 2002). The ICV of the Benxi Formation ranges from 0.12 to 3.04 (mean 0.61) and the CIA ranges from 78.56 to 95.28 (mean 90.34), this indicates that the climate was hot and humid during the deposition of the Benxi Formation shale in the study area, with strong chemical alteration (Tab.6). The relative contents of some climate-sensitive trace elements are often used to restore paleoclimate, of which C
value is a common indicator, where C
value = ∑(Fe + Mn + Cr + Ni + V + Co)/∑(Ca + Mg + Sr + Ba + K + Na). C
value > 0.8, 0.2 < C
value < 0.8, C
value < 0.2 for humid climate, semi-arid to semi-humid climate, and arid climate, respectively (
Awan et al., 2020;
Wu et al., 2022a). The C
value of the Benxi Formation shale ranges from 0.29 to 4.97 (mean 1.91), indicating that the Benxi Formation shale was deposited in a humid environment, which is similar to the results measured by the CIA parameters (Tab.6).
5.1.3 Redox conditions
Sedimentary rocks vary in their enrichment of trace elements under different redox conditions, and previous authors have commonly used U/Th and Ni/Co to reduce the redox conditions of the sedimentary environment (
Tribovillard et al., 2006;
Jenkyns et al., 2017;
Shi et al., 2022). Common redox condition proxies and discrimination criteria are shown in Tab.7. The mean values of five parameters of the Benxi Formation, Ni/Co, U/Th, δU, V/(V + Ni), and Ce/La are 3.08, 0.44, 0.98, 0.71, and 1.68, respectively, indicating that the Benxi Formation shale depositional environment as a whole exhibits a depleted oxygen reduction state (Tab.7). The cerium anomaly index (Ce
anom) can be used as a marker to judge the redox conditions of the palaeoaqueous medium, where Ce
anom = Log
10 [3Ce
n/(2La
n + Nd
n)] (n represents the standardization of Post-Archean Australian Shale). When Ce
anom > −0.1 indicates a Cerium (Ce) enriched, anoxic reducing environment, while Ce
anom < −0.1 indicates a Cerium (Ce) deficit, reflecting an oxygen-rich oxidizing environment in the water column (
Elderfield and Greaves, 1982). The Ce
anom of the Benxi Formation shale is −0.12−0.01 (mean −0.06), indicating a depositional environment dominated by reducing conditions (Tab.7).
5.1.4 Palaeosalinity
Palaeosalinity determination can be used to identify the extent and duration of seawater intrusion, analyze sedimentary facies zones and characterize the depositional environment (
Abanda and Hannigan, 2006). The ratio of Strontium (Sr) and Barium (Ba), as elements that are more sensitive to the salinity of the water column, is a more common indicator of paleosalinity (
Meng et al., 2012;
Chen et al., 2020). In fine-grained sedimentary rocks of terrestrial origin, potassium feldspar and black mica act as carriers to elemental Barium (Ba), which is relatively enriched, while the marine sedimentary environment is more favorable to Strontium (Sr) enrichment, which is due to the similar radii of Strontium ion (Sr
2+ ) and Calcium ion (Ca
2+ ). In Ca-rich carbonate rocks, Strontium ion (Sr
2+ ) similarly replaces Calcium ion (Ca
2+ ), resulting in Strontium (Sr) enrichment, and the higher the Sr/Ba ratio, the greater the palaeosalinity. Typically, the deposition environment is identified as saline water environment when Sr/Ba > 1, semi-saline water environment when 0.5 < Sr/Ba < 1, and freshwater environment when Sr/Ba < 0.5 (
Stüben et al., 2002;
Zhang et al., 2018b). Shale Sr/Ba of the Benxi Formation between 0.29 and 3.71 (mean 0.98), indicating overall high palaeosalinity and a semi-saline environment during deposition of the Benxi Formation shale (Tab.6). Th/U values can often be used to identify shale depositional environments. In general, Th/U > 7, 2 < Th/U < 7, and Th/U < 2 indicate terrestrial freshwater environments, brackish-saline sedimentary environments, and marine brackish environments respectively (
Zhang et al., 2008). Th/U values of Benxi shale are between 0.54 and 6.18 (mean 3.26), further demonstrating that the Benxi Formation has high palaeosalinity and was in a semi-saline depositional environment (Tab.6). The ratio of the content of MgO in sedimentary rocks can be used to discern the palaeosalinity of water bodies, which is usually calculated by
m = 100 × (MgO)/(Al
2O
3). Freshwater depositional environments are indicated by
m < 1, marine-continental transitional depositional environments by
m = 1−10, and seawater depositional environments by
m > 10 (
Guo, 2020). The
m values of the Benxi Formation shale range from 0.76 to 22.56 (mean 3.53; Tab.6), indicating that the salinity of the water column was high during the deposition of the shale in the Benxi Formation, consistent with the results measured for both Sr/Ba and Th/U.
5.1.5 Palaeoproductivity and volcanic hydrothermal fluid
Palaeoproductivity is often defined as the rate at which ancient marine organisms fix energy during the energy cycle. A sedimentary environment with high palaeoproductivity not only provides a large amount of organic matter, but this organic matter can consume oxygen in the water column, thus further creating an anoxic reduction environment that is conducive to the preservation of organic matter (
Bernárdez et al., 2008;
Stock et al.,2017). The P/Ti can be used as a geochemical indicator of palaeoproductivity. Phosphorus (P) is an important planktonic nutrient and is the element limiting palaeoproductivity, so phosphorus in sediments is often used as a proxy for palaeoproductivity (
Latimer and Filippelli, 2002). Benxi Formation shale with P/Ti of 0.01−1.3 (mean 0.20; Tab.6), which is higher than those of the PAAS (0.12) and UCC (0.17), indicating high palaeoproductivity of the Benxi Formation during the depositional period (
Taylor and McLennan, 1985). In recent years, numerous studies have confirmed that TOC is associated with Molybdenum (Mo) content in organic-rich sediments, and Molybdenum (Mo) content can be used to reflect the magnitude of palaeoproductivity (
Algeo and Tribovillard, 2009). The Benxi Formation shale has a Molybdenum (Mo) content of 0.75 × 10
−6−13.11 × 10
−6 (mean 3.35 × 10
−6; Tab.3), which is greater than the PAAS standard value of 1.0 × 10
−6. Both indicators suggest a high level of palaeoproductivity during its depositional period. It is known from previous experience that the main body of the shale of the Benxi Formation in the Ordos Basin is Type III organic matter derived from higher plants, and the high P/Ti values of the shale of the Benxi Formation suggest that large amounts of plankton were also present in their deposition (
Li et al., 2019c;
Chen et al., 2020;
Yang et al., 2022).
Volcanic hydrothermal fluid is rich in nutrients and can increase water temperatures, providing good conditions for the flourishing of biological communities such as bacteria and algae. Thus, hydrothermal fluid enhances biological productivity (
Qi et al., 2004;
Chen et al., 2021). The hydrothermal environment favors the deposition of Ferrum (Fe) and Manganese (Mn), while high values of Titanium (Ti) may reflect terrestrial input. Fe/Ti and (Fe + Mn)/Ti values are widely used indicators of hydrothermal sources, when Fe/Ti > 20 or (Fe + Mn)/Ti > 25 ± 5, the deposit is usually considered to be influenced by hydrothermal fluid (
Yamamoto, 1987;
Ma et al., 2019;
He et al., 2022 ). The Fe/Ti and (Fe + Mn)/Ti values of the Benxi Formation shale range from 1.21 to 77.45 (mean 15.38) and 0.23 to 78.08 (mean 16.02), respectively (Tab.6), indicating that the Benxi Formation sediments were deposited without an obvious influence from volcanic hydrothermal fluids. However, it is noteworthy that the Fe/Ti and (Fe + Mn)/Ti values of samples BX08 and BX13 are higher than the standard values, indicating local volcanic hydrothermal fluid during the deposition of the Benxi Formation shale.
5.1.6 Sedimentation rate and stagnation environment
The shale organic enrichment process is more significantly influenced by the sedimentation rate. A higher sedimentation rate can reduce the degree of rare earth differentiation. At lower sedimentation rates, the differentiation of light and heavy rare earths is stronger, and rare earth and organic matter are adsorbed by clays and other materials, while too high sedimentation rates have a dilutive effect on organic matter (
Tyson, 2001;
Tenger et al., 2006). (La/Yb)
N is a parameter used for evaluating the sedimentation rate of shales. When (La/Yb)
N is similar to 1, it indicates a high sedimentation rate and low rare earth differentiation, and When the difference between (La/Yb)
N values compared to 1 is large, it indicates a low shale deposition rate and high rare earth differentiation, the subscript N represents the standardization of NASC (
Doner et al., 2019). The (La/Yb)
N distribution of the shale samples ranges from 0.73 to 2.98 (mean 1.83; Tab.5), indicating a high deposition rate and low rare earth differentiation in the Benxi Formation shale.
In stagnant deepwater environments, bacterial sulfate-reducing environments form hydrogen sulfide, which enriches Mo in the sediments while depleting it from seawater. Three types of stagnation environment can be derived from previous experience, Mo/TOC < 4.5, 4.5 < Mo/TOC < 45 and Mo/TOC > 45 representing strong, semi-stagnation and weak stagnation environment, respectively (
Algeo and Rowe, 2012). The Mo/TOC of Benxi Formation shale samples ranges from 0.36 to 9.93 (mean 1.31; Tab.6), indicating deposition in a strong stagnant environment. This is consistent with the report of a stagnant environment for the Benxi Formation shale in the Zhaoxian section and well M115 of Linxian, Shanxi, eastern Ordos Basin (Fig.4) (
Cui et al. 2022).
5.2 Organic matter enrichment mode of the Benxi Formation shale
It is clear from the above that the shale of the Benxi Formation was deposited in a hot and humid climate with high palaeoproductivity and local volcanic hydrothermal fluid, accompanied by strong terrestrial clastic input and a high sedimentation rate and stagnation degree. The bottom water as a whole was in oxygen-poor reducing conditions and a semi-saline deposition environment. The organic matter enrichment of shale is influenced by multiple factors. The palaeoenvironment (palaeoclimate, redox conditions, palaeosalinity, etc.) is a fundamental factor influencing the diversity, growth, and decomposition of organisms, controlling the production and preservation of organic matter (
Talbot, 1988;
Chen et al., 2019). To determine which or which factors control their organic enrichment, this study used TOC to represent the degree of organic matter enrichment and analyzed its relationships with terrestrial clastic input indicators (Aluminum (Al), Titanium (Ti), Thorium (Th)), palaeoproductivity indicators (P/Ti, Molybdenum (Mo)), volcanic hydrothermal fluid indicators (Fe/Ti, (Fe + Mn)/Ti), palaeoclimatic indicators (CIA, C
value), water body redox condition indicators (U/Th, δU), palaeosalinity indicators (Sr/Ba,
m (
m = 100 × (MgO)/(Al
2O
3)), sedimentation rate indicators ((La/Yb)
N) and stagnant environment indicators (Mo/TOC).
Terrestrial clastic and sedimentation rate are important factors affecting organic matter enrichment. Terrestrial clastic input can transport terrestrial organic matter into water and enrich organic matter accumulation. However, excessive input of terrestrial-sourced debris material can have a diluting effect on productivity. Aluminum (Al), Titanium (Ti), and Thorium (Th) elements are often used to indicate the degree of influence of clastic input from terrestrial sources (
Tyson, 2001;
Ding et al., 2015). Aluminum (Al), Titanium (Ti), and Thorium (Th) elements of the Benxi Formation shale samples all show negative correlations with TOC (
R2 = 0.36, 0.51, and 0.13, respectively; Fig.5(a)). The influx of detritus provides abundant organic matter (plant decomposition) on the one hand, and acts as a diluent of organic matter, on the other hand, indicating that the excessive terrestrial detrital material infusion has a dilution effect while providing sufficient organic matter, and the dilution effect exceeds its organic matter supply. Thus, the terrestrial clastic is not an influencing factor for its organic matter enrichment. The sedimentation rate indicators ((La/Yb)
N) are not significantly correlated with the terrestrial clastic input indicators (Aluminum (Al), Titanium (Ti), Thorium (Th)) and is weakly correlated with TOC, while the stagnant environment indicators (Mo/TOC) are weakly negatively correlated with TOC (
R2 = 0.25; Fig.5(b)−(c)). This suggests that the input of terrestrial clastic from the Benxi Formation shale is independent of the sedimentation rate and that both sedimentation rate and sedimentation rate have a weak influence on shale organic matter enrichment.
The palaeoclimatic indicators (CIA, C
value) of the Benxi Formation shale correlate extremely weakly with TOC (Fig.5(d)), suggesting that hot and humid climatic conditions do not play a significant role in promoting biological growth and reproduction, therefore palaeoclimate is not the main controlling factor for organic matter enrichment in the Benxi Formation shale. This trend differs from those reported for the Shanxi Formation and Longtan Formation of the marine-continental transitional facies (
Liu et al., 2018;
Wei et al., 2020), which indicates that under equally humid conditions, excessive temperatures could probably inhibit organic enrichment in the shale of the Benxi Formation.
The TOC of the Benxi Formation shale is positively correlated with the P/Ti ratio and Molybdenum (Mo) content (
R2 = 0.39 and 0.14, respectively; Fig.6(a)), indicating that higher palaeoproductivity provided a certain material basis for the organic matter enrichment during the deposition of the shales in the Benxi Formation and that palaeoproductivity is one of the influences on the organic-rich shale. This trend is consistent with previous reports on Carboniferous-Permian marine-continental transitional facies rocks of Ordos Basin (
Wei et al., 2020;
Zhang et al., 2021). In this study, the TOC content was weakly positively correlated with (Fe + Mn)/Ti and Fe/Ti (
R2 = 0.18 and 0.27, respectively; Fig.6(b)), indicating that local volcanic hydrothermal fluid provided certain nutrients to Benxi Formation shale and elevated water temperatures, which could provide better conditions for biological communities such as mycorrhizal algae and, in turn, promote organic matter enrichment.
The palaeowater depth indicators (Mn/Fe, Mn/Ti) and redox condition indicators U/Th of the Benxi Formation shale are significantly positively correlated (
R2 = 0.79 and 0.84, respectively; Fig.7(a)) and the water body redox condition indicators (U/Th, δU) show a strong correlation with TOC (
R2 = 0.63 and 0.47, respectively; Fig.7(b)). This indicates that with the deepening of palaeowater, the aqueous reducing conditionality of the Benxi Formation shale was enhanced. The moderate palaeowater depths of the marine-continental transitional shale depositional environment, result in an oxygen-poor reducing aqueous environment. To some extent, the reducing environment occurred between the continental and marine shale, and the oxygen-poor reducing environment favors the preservation of organic matter and promotes organic matter enrichment. It is consistent with the trends of the Niutitang and Longmaxi Formations measured by Li et al. (2019a) and the Shahejie Formation measured by
Hu et al. (2018). Currently, redox-sensitive elements (RSEs) enrichments, for instance, Mo and U, can be used to analyze the depositional environment of shales. The Mo
EF of the Benxi Formation shale ranges from 0.20 to 39.92 (mean 3.69) and the U
EF ranges from 0.44 to 7.57 (mean 1.72), with a high positive correlation between U
EF and Mo
EF (
R2 = 0.79; Fig.7(c)). This reveals that Mo and U elements are highly enriched and co-occur, forming in an anoxic-sulfidic deposition environment
(Tang et al., 2015). The TOC shows a good positive correlation with both Mo
EF and U
EF (
R2 = 0.55 and 0.65, respectively; Fig.7(d)), and multiple indicators suggest that bottom-water anoxic reduction conditions favor organic enrichment in the shales of the Benxi Formation.
The palaeosalinity indicators (Sr/Ba,
m (
m = 100 × (MgO)/(Al
2O
3))) of the Benxi Formation shale are both proportional to TOC (
R2 = 0.47, 0.49, respectively; Fig.8), manifesting that higher palaeosalinity has a significant impact on the growth and preservation of palaeontological. The higher salinity water environment favors the preservation of organic matter, which is close to previous work on the organic matter enrichment of transitional shales in the north-eastern Ordos Basin (
Chen et al., 2020).
From the above analysis, it is clear that the enrichment of organic matter in the marine-continental transitional shale of the Benxi Formation is the result of the interaction of several factors, such as palaeoproductivity, volcanic hydrothermal fluid, water body redox conditions, and palaeosalinity. All of the above factors directly or indirectly influence the production and preservation of organic matter in shale, further contributing to organic enrichment. Among these factors, redox conditions have the greatest influence, while palaeosalinity and palaeoproductivity promotion are second, and volcanic hydrothermal fluids have a weaker influence (Tab.8).
Combined with the above analyses and the regional sedimentary tectonic setting, a organic matter enrichment mode of the Benxi Formation shale can be restored (Fig.9(a)). These shales are of a marine-continental transitional facies and were deposited on the continental margin. The hot and humid environment during the deposition period promoted biological growth and accelerated rock weathering. Input of terrestrial clastic and local volcanic hydrothermal fluid also promoted biological reproduction, so that the Benxi Formation shale had high palaeoproductivity and formed abundant organic matter. An oxygen-poor and high-salinity bottom-water environment is conducive to organic matter preservation, while biological respiration can increases the consumption of oxygen in the water column. The high sedimentation rate and strong stagnation environment also reduce the time for oxidative decomposition of organic matter in the water column, allowing for rapid deposition and burial and creating a unique mode of organic enrichment in this formation and sedimentary facies.
5.3 Differential enrichment mechanism of organic-rich shale under sedimentary facies
In this study, the geochemical characteristics of organic matter-rich shale under three different sedimentary facies were cited, compared, and analyzed (where the Longtan and Shanxi formations are marine-continental transitional facies, the Chang 7, Chang 9, and Shahejie formations are derived from continental facies, and the Dalong, Sowa, and Niutitang formations are marine facies) (
Zeng et al., 2015;
Zhang et al., 2017;
Chen et al., 2021;
Wu et al., 2021;
Zhang et al., 2021;
Wu et al., 2022b;
Xu et al., 2022). The specific parameters are shown in Tab.9. Combining their production factors (terrestrial clastic input, palaeoclimate, palaeoproductivity, and volcanic hydrothermal fluid) with preservation factors (redox conditions, palaeosalinity, sedimentation rate, and stagnation environment), so the differences in depositional environments and organic matter enrichment mechanisms of shale in three different sedimentary facies are proposed (Fig.9 and Fig.10).
The organic-rich shale of the continental, marine-continental transitional, and marine facies are often deposited in the semi-deep lacustrine – deep lacustrine, lagoonal, and Shelf sedimentary facies. The marine-continental transitional tidepool-lagoon deposits are located closer to land and are more influenced by the land, with the highest Al
2O
3 and TiO
2 contents among the three. The Al
2O
3 and TiO
2 contents of organic-rich shale under the three sedimentary facies are as follows: marine-continental transitional shale (average of 22.45 wt% and 0.78 wt%) > continental shale (average of 14.23 wt% and 0.49 wt%) > marine shale (average of 8.36 wt% and 0.34 wt%; Tab.9, Fig.10(a)). In terms of paleoclimate (CIA, C
value), the values under each depositional facies are as follows: marine-continental transitional shale (average of 86.44 and 0.89) > marine shale (average of 71.64 and 0.42) ≈ continental shale (average of 68.47 and 0.38). This indicates that the main body of the marine-continental transitional shale was deposited in a warm, humid and hot, humid climate with strong chemical weathering, whereas the terrestrial and marine shale deposition environments were in a cold, dry to warm, humid state, but both were dominated by the warm, humid environment (Tab.9, Fig.10(b)). It can be seen that the stronger chemical weathering and warm and humid climatic conditions are more conducive to shale organic matter enrichment. The global climate change curve shows that the palaeoclimate during the deposition of continental and marine shales was in the vicinity of the warm and cold climate divide and the warm region, in agreement with the CIA parameter results(Fig.10(j)) (
Liang et al., 2020). It is important to note that the Carboniferous-Permian period is predominantly located in a cold region and the period as a whole was a major ice age, but global mean temperatures rebounded near the junction of the late Carboniferous and early Permian, and the marine-continental transitional shale in the study area developed mainly during this facies. Previous studies have shown that the southern hemisphere had a cold glacier-covered climate during this period, while the northern hemisphere had a warm and humid tropical and subtropical climate, and northern and southern China were in the northern hemisphere and close to the equator at this time (Fig.1(a)) (
Liang et al., 2020;
Scotese et al., 2021;
Zhang et al., 2022). So the climate was warm and humid to hot and humid, which is consistent with the results shown by geochemical parameters in this study. This suggests that although the palaeoclimate of organic matter-rich shale in the study area in the marine-continental transitional facies is not consistent with the global mean temperature, it is still an important factor in the enrichment of organic matter in shale (Tab.9, Fig.10(j)). In terms of palaeoproductivity (P/Ti, Molybdenum (Mo)) and volcanic hydrothermal fluid (Fe/Ti, (Fe + Mn)/Ti), the differences between the three sedimentary facies are relatively small, with P/Ti and Molybdenum (Mo) values greater than 0.1 and 1 ppm, respectively, and Fe/Ti and (Fe + Mn)/Ti values both close to 10 on average, with some values greater than 25, indicating that the palaeoproductivity corresponding to their organic-rich shale is at a high level, and the volcanic hydrothermal fluid is often associated with relative plate movement and developed at plate boundaries, which is not significantly affected by differences in sedimentary facies (Tab.9, Fig.10(c) and Fig.10(d)). The marine shale is in the shelf facies, where the water column is the most reductive, and the bottom water is in a dysoxic-anoxic state during its deposition. In contrast, during the deposition of marine-continental transitional and continental shale, the water column is in weakly oxic-dysoxic conditions, but both are predominantly dysoxic. It is clear that the strongly reductive water column is more conducive to shale organic matter enrichment. The redox conditions (Ni/Co, δU) of the three are as follows: marine shale (average of 5.23 and 1.53) > marine-continental transitional shale (average of 3.00 and 0.94) ≈ continental shale (average of 2.46 and 1.31). Meanwhile, the paleosalinity (Sr/Ba) of the shale under the three sedimentary facies is freshwater, semi-saline, and saline respectively, showing: marine shale (average 4.85) > marine-continental transitional shale (average 0.74) > continental shale (average 0.49; Fig.10(e) and 10(f))). From the analysis of both deposition rate and stagnation degree, it can be seen that the organic-rich shale (La/Yb)
N in all three sedimentary facies lies around 1 (Fig.10(g)), and the stagnation environment indicators (Mo/TOC) is less than 45, which indicating high deposition rate and strong stagnation and semi-stagnation characteristics (Tab.9, Fig.10(h)). These mean that the high deposition rate and strong stagnation environment can promote shale organic matter enrichment to some extent. As can be seen from Tab.9 and Fig.10(i), there is no significant correlation between palaeowater depth and the various depositional facies, which may be related to the location of deposition of the organic-rich shale.
From the above analysis, it can be seen that the differences in depositional environments of organic-rich shale under the three sedimentary facies are mainly reflected in four aspects of production factors (terrestrial clastic input and palaeoclimate) and preservation factors (palaeosalinity and redox conditions). Based on the differences between shale depositional environments, it is then possible to construct modes of organic matter enrichment during periods of shale deposition under different sedimentary phases. The continental shale was mainly developed in lake basins, where the climate was warm and humid during its deposition. The main body of bottom water was dysoxic, and the palaeosalinity was low. Previous experience showed that continental shale was strongly influenced by the continent at the lake margins and therefore has a high input of terrestrial clastic, which decreases as it extends toward the center of the lake and is moderate overall. The period of deposition was characterized by high levels of productivity and local volcanic hydrothermal fluid. By comparing the organic matter enrichment mode of the continental shale at different formations, it is clear that they are mainly influenced by palaeoproductivity, terrestrial clastic, and volcanic hydrothermal fluid, from which it can be grouped into “production mode” (Fig.9(b)) (
Zhang et al., 2017;
Chen et al., 2021;
Wu et al., 2021). The marine shale was mainly developed in the land shelf, and its deposition period had a warm and humid climate, while the water body was highly reductive and in a saline state. From previous analyses, it is clear that the marine shale was weakly influenced by continental influences, resulting in a low input of terrestrial clastic, while palaeoproductivity was high and possibly influenced by volcanic hydrothermal fluid in the plate junction region. The organic matter enrichment is influenced by several factors, including redox conditions and palaeosalinity, which are more controlled, while the influence of palaeoproductivity and volcanic hydrothermal fluid is relatively weak, showing a predominantly “preservation mode” (Fig.9(c)) (
Zeng et al., 2015;
Chen et al., 2021; Lan and Shen, 2022;
Wu et al., 2022b;
Xu et al., 2022). The main depositional facies of the marine-continental transitional shale was lagoons and tidal flats, and therefore close to land, resulting in high input of terrestrial clastic and an overall hot and humid palaeo-climate for the marine-continental transitional shale. Likewise, the palaeosalinity of the marine-continental transitional shale was high and intermediate between the continental and marine shale, being semi-saline, and the water bodies were in a reduced state of dysoxic during deposition. Volcanic hydrothermal fluid was present at the plate junction at the same time as the palaeoproductivity was generally high. Organic matter enrichment is influenced by palaeoproductivity, palaeosalinity, and redox conditions, showing a “hybrid mode of preservation and production” (Fig.9(a)) (
Liu et al., 2018;
Zhang et al., 2021;
Wu et al., 2022b).
6 Conclusions
This study focused on the geochemical characteristics of the Upper Carboniferous Benxi Formation shale in the east-central Ordos Basin. The palaeoenvironment during their depositional period was reconstructed and the organic matter enrichment mechanism was revealed. The conclusions reached are shown below.
1) The Benxi Formation shale main elements in the study area are dominated by SiO2, with Al2O3, TiO2, and Fe2O3 being relatively enriched. Trace elements such as V, Cr, Co, Ni, Mo, U, and Y are enriched relative to UCC, while Zn, Zr, Sr, Ba, Ta, Hf, and Th are relatively deficient. Total rare earth elements (ΣREE) vary from 50.48 to 621.49 ppm (average 274.23 ppm), which is higher than in the NASC and UCC, and slightly higher than in the PAAS. The REE distribution pattern shows a right-sloping trend for light REEs and a flattening trend for heavy REEs.
2) The tectonic setting of the Benxi Formation shale source area in the study area is mainly continental margin with a shallow palaeowater depth. It belongs to marine-continental transitional facies, where the deposition period was accompanied by strong terrestrial clastic input. The source rocks are upper crustal granite and sedimentary rocks.
3) The Benxi Formation shale was deposited in a hot and humid climate, with high palaeoproductivity and local volcanic hydrothermal fluid. At the same time, the sedimentation rate and stagnation environment of the water column were high, and the bottom water was generally in dysoxic conditions and a semi-saline sedimentary environment.
4) Multiple factors, such as palaeoproductivity, volcanic hydrothermal fluid, redox conditions, and palaeosalinity are considered to interact in influencing the enrichment and preservation of organic matter of the marine-continental transitional facies. Among these, the redox conditions have the greatest influence on organic matter preservation, while palaeosalinity and palaeoproductivity facilitate the preservation and enrichment of organic matter. The organic matter enrichment mode is a composite model combining the “preservation mode” of oxygen-poor reductive and semi-saline environment, and the “production mode” of palaeoproductivity and local volcanic hydrothermal fluid.
5) There are significant differences in the depositional environments and organic matter enrichment mechanisms of terrestrial, marine transitional, and marine organic-rich shale. The differences in depositional environments are mainly reflected in the four aspects of terrestrial detritus, palaeoclimate, redox conditions, and palaeosalinity, and the organic matter enrichment modes can be categorised as “production mode”, “hybrid mode of preservation and production”, and “preservation mode”, respectively.