Abstract

How the geochemical characteristics of organic matter shale for the carbonate platform facies remain uncertain, which restricts an integrated reconstruction of the model of organic matter accumulation and preservation. Here, we present new results from element geochemical fingerprinting and integrated analyses of paleoclimate, paleoproductivity, paleoredox environment, and terrigenous input of the targeted Middle Devonian Tanshanping and Dacaozi Formations in the Ninglang-Yanyuan Basin, Southwest China. It is worth noting that although the carbonate platform connects with the open sea partially, the redox environment will not be completely controlled by relative sea level variations. Specially, the paleoclimate, paleoproductivity, and paleoredox conditions are the main controlling factors of the accumulation and preservation of organic matter. In view of the paleoclimate indexes, we suggest that both a relatively warm-humid climate characterized by intensified chemical weathering conditions and a higher terrigenous input are identified as two major drivers forcing the reductive environment in the sedimentary waterbody. Finally, a comprehensive model is established for providing new insights into the mechanism of organic matter accumulation and preservation for the carbonate platform facies. The paleoredox environment, paleoproductivity, paleoclimate, and terrigenous input are believed to have exerted a very considerable force on reconstructing the model of organic matter accumulation and preservation for the carbonate platform facies. Specially, the coupling interactions between the paleoproductivity and redox condition are thus also stressed. We found that the preservation condition is much more important than the paleoproductivity, resulting in the degree of organic matter enrichment. Even if the paleoproductivity of a sedimentary waterbody of a depositional period of the Dacaozi Formation was higher, the TOC concentrations were relatively low due to the poor preservation condition by fall of the sea level and increase of the terrigenous input. In another aspect, the better preservation condition of the Tanshanping Formation makes the TOC concentrations higher in the case of lower paleoproductivity in the sedimentary waterbody.

1. Introduction

The mud shale of the Tanshanping Formation in the Ninglang-Yanyuan area is mainly composed of clay minerals and carbonate minerals, and the mud shale of the Dacaozi Formation is mainly composed of clay minerals and siliceous minerals. As organic geochemical characteristic analysis shows, the organic matter types in the middle Devonian of the Ninglang-Yanyuan region are type II1 kerogen, reflecting that the maturity is high which is in the overmature dry gas phase, and the maturity increases with the increase of burial depth [1, 2]. Previous studies have proved that the Tanshanping Formation and the Dacaozi Formation are the shale-concentrated sections in the Ninglang-Yanyuan Basin, and they have been important favorable exploration layers for shale gas, which is characterized by large thickness, high maturity, and high organic carbon content. As the development of the shale exploration in the Tanshanping Formation and the Dacaozi Formation, combined with the deployment of two-dimensional seismic data, we chose a structure favorable site (Lizihe anticline) to deploy a geological survey well called Well Yunningdi-1. Through the analysis and testing, the well drilling by the stable natural gas in the Devonian strata shows that the average total content of analytical gas is 3.25 m3/t. Specially, the total organic carbon content in the Tanshanping Formation was between 0.09% and 1.83% (with an average of 1.14%), and the organic carbon content in the Dacaozi Formation was between 0.32% and 3.48% (with an average of 1.05%).

A series of studies have reported that the Ninglang area in Yunnan is an important area for shale gas exploration and development in the complex tectonic area of Southern China. The Tanshanping Formation and the Dacaozi Formation in the Middle Devonian are components of the main source rocks. Previous studies have shown that the hydrocarbon generation capacity of shale reservoirs depends on the content, type, and thermal evolution of organic carbon [3, 4]. The development of the shale reservoir space is also affected by organic carbon content to some extent [3]. At the same time, the shale which is rich in organic matter tends to have high gas content, which indicates that the organic carbon content has a good correlation with the adsorption performance of shale [59]. Therefore, these important scientific issues put forward the necessity to look for more effective evidences to correctly clarify the main control factors of organic matter accumulation and preservation and the resource potential evaluation of Paleozoic shale gas in the Ninglang-Yanyuan area, Yunnan. Many scholars all over the world believe that the favorable conditions for organic matter accumulation in black shale mainly include three modes: high productivity, anoxic preservation, and the combination of the two [1016]. The model of high productivity is a hypoxic environment which is conducive to the preservation of organic matter when high organic matter is formed under the influence of paleoclimatology, terrestrial organic matter supply, biological productivity of waterbody, and other factors by consuming oxygen in the sedimentary waterbody [1719]. To sum up, the three models are in the final analysis to discuss the effects of the paleosedimentary environment formed by different factors on the accumulation and preservation of organic matter. Therefore, the reconstruction of the paleosedimentary environment is very important. Many scholars have used the geochemical index of shale to study the sedimentary environment of the organic rich shale [2025].

At present, most scholars have studied the shale outcrops of the Tanshanping Formation and the Dacaozi Formation of the Devonian in the Ninglang-Yanyuan area in terms of geochemical characteristics, generation and storage characteristics, sedimentary characteristics, and shale gas geological conditions. However, the mechanism of organic matter accumulation and preservation in black shale has not been studied in the Tanshanping Formation and the Dacaozi Formation. In this paper, we took the WellYunningdi-1 as an example, and systematic elemental geochemistry analysis was conducted, including the analysis of main elements by X-ray fluorescence spectrometry, the measurement of trace elements and rare earth elements by inductively coupled plasma-mass spectrometry (ICP-MS), the combustion oxidation-nondispersive infrared absorption method, and the measurement of TOC concentrations by the total organic carbon analyzer and summary of its variation characteristics. Moreover, we have studied the elemental chemical characteristics and paleosedimentary environment of the Tanshanping and Dacaozi Formations. On this basis, the main controlling factors of accumulation and preservation of organic matter in shale rocks were discussed in order to provide a theoretical basis for the mechanism of accumulation and preservation of organic matter and the reconstruction of sedimentary model.

2. Geological Background

2.1. Tectonic Position

The Yanyuan Basin is located in the west of Sichuan Basin, at the junction of the Sichuan-Yunnan plateau, and is structurally located in the south-north tectonic belt between western Sichuan and eastern Tibet. Specially, the basin was sandwiched between the Songpan-Ganzi Indosinian fold belt and the Kang-Yunnan ancient land, which is part of the Zhenzhon gnappe thrust belt. The northern part of the Yanyuan Basin is controlled by the Xiaojinhe fault, which is adjacent to the Sanjiang island arc orogeny. The west side is controlled by the Guobaoshan-Yanfeng fault, which divides the Yanyuan Basin, the Yongning fault block on the west side and the Ninglang Basin on the southwest side. The east side is controlled by the Jinhe-Jinghe fault and adjacent to the Kangdian ancient land. The Lijiang Basin lies to the south of the Kang-Dian ancient land [26].

2.2. Tectonic Characteristics

The tectonic evolution of the Upper Yangtze plate is mainly controlled by the tectonic activities of the Tethys tectonic domain, which has undergone three major development stages: the proto-Tethys Oceanic basin (Z-O-S), the Guthys Oceanic basin (D-T3), and the New Tethys Oceanic basin (T3-N), reflecting the formation and evolution of the Tethys Ocean basin showed the change of becoming more and more earlier from south to north [27, 28]. Evidence of the disappearance of the New Tethys basin can be seen in the Yaluzangbu river and the Bangongcuo-Nu river Ophiolitic jumble zone in Tibet. At this time, the marine deposition in the western part of the Upper Yangtze river was terminated by the Indochina movement, which was no longer in the category of the tectonic evolution of the New Tethys Basin but was mainly located in the process of sea-land transition facies-continental sedimentary evolution. Paleozoic and Mesozoic marine sedimentary events in the western part of the Upper Yangtze Plate are bounded by the end of Silurian, and the geological period can be divided into two main periods:(1) Proto-Tethyan period of the Sinian-Early Paleozoic and (2) Late Paleozoic-Early Mesozoic Paleo-Tethyan period. Based on the observation and analysis of a large number of drilling data and outcrops, the sedimentary sequence and structural characteristics of the Yanyuan Basin have been discussed. It is believed that the Yanyuan Basin was formed during the Miocene to Pliocene by the collision and extrusion of the Indian plate and the Eurasian plate. The emergence of this basin indirectly supported and recorded an important collision event of the Indo-Eurasian Plate in the Late Cenozoic [29]. It is important to focus on the tectonic position of the Yanyuan Basin which is in the arc structural belt in front of the Sichuan-Yunnan block, and it is the most obvious area of Cenozoic extrusion in the Qinghai-Tibet Plateau, which is sandwiched between the Xianshuihe fault and Honghe fault. The Sichuan-Yunnan block and its interior were cut by the Xiaojinhe fault, which can be divided into two secondary part blocks called the northwestern Sichuan Block and central Yunnan Block. The sliding velocities of the block boundaries are different between the blocks in northwestern Sichuan and central Yunnan, which reflects the composite motion model of the secondary blocks sliding horizontally and rotating clockwise about the vertical axis [30].

2.3. Sedimentary and Stratum Characteristics
2.3.1. The Ninglang Area

According to the “Yunnan Provincial Regional Geology Records of Ministry of Geology and Mineral Resources,” the Silurian (S) strata are mainly distributed in the Dacaozi area in the northwest and east of the Ninglang Basin, with thickness ranging from 182 m to 659 m. The Lower Silurian is mainly composed of shallow marine facies, with argillaceous limestone and gray limestone shale with a small amount of black siliceous layers. The Upper Silurian is mainly composed of lagoon facies, and a large number of carbonate rocks were developed. As for the Devonian system, Devonian strata are well developed in this area, which mainly outcrops in the Dacaozi, Xilaping, Wenxiang, and other areas, and there are deposits in the Upper, Middle, and Lower Series (Figure 1). The measurement area is divided from bottom to top as follows: The Lower Devonian Daguaping Formation (D1d), with a thickness of 92 m~687 m, is exposed in the eastern Dacaozi and Maoniuping areas, which belongs to coastal and shallow marine facies deposition and be dominated by dark gray shale and siltstone at the top and conglomeratic sandstone at the bottom. It is worth noting that the lithology is in a regular cycle and in parallel unconformable contact with the Silurian system, reflecting an obvious depositional discontinuity. The thickness of the outcrops in the eastern part of the basin is about 91 m, and it is missing in the Xilaping area in the southeast. The Middle Devonian Dacaozi Formation (D2d) is exposed at Dacaozi and Maoniuping in the east, with a thickness of 266 m~378 m, in which thickness is 266 m at the Dacaozi and missing at the Xilaping in the southeast. The strata of the Dacaozi Formation belong to shallow marine sedimentary facies, which is dominated by dark gray and gray biolimestone, which is interbedded with shale and integrated over the Daguaping Formation of Lower Devonian. The Middle Devonian Tanshanping Formation (D2t) is exposed in the Dacaozi area in the east with a thickness range of 618 m~640 m, and the thickness of the Dacaozi area is 141.3 m, which belongs to shallow marine sedimentary facies. The lithology is quite different from bottom to top. In particular, the bottom is mainly interbedded with sandstone, conglomeratic sandstone, and black shale, the upper part has experienced a transition that it has less sandy and more calcareous, and the top part is mainly grey limestone, bioclastic limestone, marl, and shale. The thickness of the Tanshanping Formation is about 617 m, and it is in a conformable contact with the Dacaozi Formation. On the whole, the Devonian strata overlay the Lower Ordovician strata. The Upper Devonian Lanniqing Formation (D3l), with a thickness ranging from 476 m to 569 m, is exposed in the Dacaozi, Xilaping, and Luniqing areas in the east, composed of limestone at the top and dolomitic limestone at the bottom, which is integrated above the Tanshanping Formation. The Upper Devonian Gangou Formation (D3g), with a thickness of 298 m~527 m, is exposed in the Dacaozi, Maoniuping, and other areas in the east, which belongs to shallow marine sedimentary facies, and its lithology is dominated by limestone and oolitic limestone, which is integrated over the Lanniqing Formation. Moreover, the Carboniferous system (C), with a thickness range of 266 m~964 m, is exposed in Yaoshan, Zhushan, and other areas in the northwest of China, in which lithology is dominated by dark gray bioclastic limestone and oolitic limestone, with pebbly sandstone and conglomerate at the bottom and chert nodule oolitic limestone at the top [31].

2.3.2. The Yanyuan Area

The strata in the Yanyuan Basin are all exposed from Sinianto Quaternary, except for the absence of Jurassic and Cretaceous. The strata below Triassic are all marine sediments, and the upper Permian has a set of Emei mountain basalt intrusion. The Sinian strata in the basin are exposed in the southeast margin of the basin, which only has a part of the upper Sinian strata; the lithology is mainly dolomite, and the bottom of which is marl, sandy mudstone interlayer. The Lower Paleozoic strata are exposed in the southeast and southwest sides of the basin, and the upper part of the Lower Cambrian to upper Cambrian and upper Ordovician strata are missing, so the Cambrian and Ordovician, Ordovician and Silurian are in parallel unconformable contact. The stratigraphic lithology of the Cambrian is dominated by fine silty clastic rocks, while the dolomite and limestone of the Middle and Upper Ordovician belong to a relative high level of enrichment, and the top is in parallel unconformable contact with the Silurian. As for the Lower Silurian Longmaxi Formation, it belongs to the marine deposit, which consists of graptolite-bearing siliceous argillaceous shale, intercalated argillaceous limestone, argillaceous sandstone, and sandy mudstone. The Middle Silurian is mainly argillaceous limestone, while the Upper Silurian is sand-mudstone and siliceous rock, and the Middle and Upper Silurian is locally dolomite. In terms of the Devonian system, Devonian and Silurian, Carboniferous and Devonian, and Permian and Carboniferous are in parallel unconformity contacts. The lithology of the Lower Devonian is mainly quartz sandstone, fine siltstone, limestone, siliceous rock intercalated with shale, argillaceous siltstone, and black shale. In the upward direction, the gray and sandy matter are reduced, while the argillaceous matter is increased. The Middle and Upper Devonian is dominated by carbonate deposits interbedded with calcareous mudstone and argillaceous limestone, while it contains a small amount of sandstone conglomerate. Then, the Carboniferous system is mainly composed of limestone deposits with light gray limestone lithology and a large number of coral fossils. In the end, the lower Permian is dominated by sandy mudstone with local coal intercalation lines, and the middle Permian is dominated by limestone with weathering crust developed at the top, and the upper Permian is invaded by basalt with a thickness of about 2800 m. The top is dominated by fine clastic rocks with black carbonaceous shale mainly.

The Ninglang-Yanyuan Basin is located at the junction of the Songpan-Ganzi fold system and the western margin of the Yangtze platform. It is a component of the Yanyuan-Lijiangtai fold belt on the western margin of the Upper Yangtze block, adjacent to the Sanjiang tectonic belt in the west and the Kangdian ancient land in the east. Accordingly, this study focuses on Well Yunningdi-1,which is selected as the favorable site for structural preservation in this study (Figure 2).

3. Samples and Methods

3.1. Samples

A total of 52 shale samples in the Tanshanping and Dacaozi Formations were obtained from Well Yunningdi-1 in the Ninglang-Yanyuan Basin. The depth of these shale samples was in a range of 750.28 m~1486.5 m. Specially, the data of shale gas content was provided by the Oil & Gas Resource Survey, China Geological Survey. In this study, 52 shale samples were systematically collected from the black shales of the Devonian Tanshanping and Dacaozi Formations in Well Yunningdi-1 in the Ninglang-Yanyuan Basin. In order to take the next step, pretreatment was performed firstly. After drying the sample, which had been cleaned by ultrasonic waves with deionized water several times, the sample was ground to 200 mesh with agate mortar.

3.2. Analytical Methods
3.2.1. TOC Analysis

In this study, a total of 52 shale samples from the black shales of the Devonian Tanshanping and Dacaozi Formations were selected for TOC analysis. About 100 mg of powdered samples was prepared and then acidified with 5% HCI in a crucible when the temperature reached approximately 80°C to remove carbonates. The TOC concentration was determined using a Leco CS-230 apparatus at the Oil & Gas Resource Survey, China Geological Survey.

3.2.2. Major Element

All 52 samples were air-dried and crushed into powders before geochemical analysis. The concentration of major elements were measured at the Oil & Gas Resource Survey, China Geological Survey, using a Leeman Prodigy (USA) inductively coupled plasma-optical emission spectrometer (ICP-OES) with high-dispersion Echelle optics. Samples were prepared by dilution with tetraborate/lithium metaborate fusion and nitric dissolution processes. Chinese National Rock Standards (GSR-1 and GSR-5) and US Geological Survey (USGS) standard materials (BCR-1, AVG-2, and SCO-1) were used to correct the analytical uncertainties. Analytical precision for most of the major element oxides is more than 98%.

3.2.3. Trace and Rare Earth Elements

Trace element and rare earth element concentrations were determined with the Thermo Scientific Element XR Sector Field Inductively Coupled Plasma Source Mass Spectrometer (ICP-MS). Firstly, 52 test samples were added to lithium metaborate and lithium tetraborate solvent. Then, the powder samples were mixed well and dried at 105°C. Finally, ICP⁃MS was used to complete the test, which accuracy of for most of the trace and rare earth elements was better than 5%, with higher precision.

4. Results

4.1. Distribution of the TOC

The total organic carbon concentrations of the Tanshanping Formation and the Dacaozi Formation in Well Yunningdi-1 range from 0.10% to 3.48%, with an average of 1.08%. In particular, the total organic carbon concentrations of the Tanshanping Formation are higher than that of the Dacaozi Formation means 1.13% and 1.06%, respectively). Among them, it is worth noticing that the total organic carbon content of the organic-rich shale (1070.6~1088.6 m) in the Tanshanping Formation ranges from 0.92% to 1.43%, with an average of 1.21%. And the total organic carbon content of organic-rich shale (1320.8~1486.5 m) in the Dacaozi Formation ranges from 0.37% to approximately 3.48%, with an average of 1.14%. They were thought to be a good source rock because of the relatively high organic carbon content of these Formation shales, which has reached the standard of good source rock.

4.2. Distribution of Major Elements

The contents of major-element oxides in the shale samples from the Tanshanping Formation and the Dacaozi Formation are presented in Table 1. Table 1 shows that the main elements of the Tanshanping shale are SiO2(with an average of 42.55%), Al2O3 (with an average of 12.17%), and CaO (with an average of 11.94%), and the total content of these 3 components ranges from 48.96% to 78.34%. In addition, the concentrations of other major element oxides in the Tanshanping Formation are obviously less than 4% ,which include Fe2O3 (with an average of 3.41%), MgO (with an average of 1.79%), K2O (with an average of 2.87%), Na2O (with an average of 1.26%), P2O5, SO3, TiO2, and MnO (less than 1%). It is worth noting that the main components of the organic-rich shale in the Tanshanping Formation are SiO2 (with an average of 44.74%), Al2O3 (with an average of 12.95%), and CaO (with an average of 11.97%), and the total content of these three components ranges from 63.98% to 73.61%.

The main elements of the Dacaozi shale are SiO2 (with an average of 56.07%), Al2O3 (with an average of 15.95%), Fe2O3 (with an average of 5.89%), and CaO (with an average of 4.57%). And the total content of these 4 components ranges from 80.74% to 96.92%. The concentrations of other major element oxides in the Dacaozi Formation are obviously less than 6%, which include K2O (with an average of 5.53%), MgO (with an average of 1.54%), Na2O (with an average of 1.44%), P2O5, SO3, TiO2, and MnO (less than 1%).

The variation trend of the mass fraction of major elements and TOC is similar to the study section of the Tanshanping Formation and the Dacaozi Formation (Figure 3). In the Tanshanping organic-rich shale, the contents of SiO2, Al2O3, Fe2O3, Na2O, TiO2, and MnO decrease with the increase of TOC, while the contents of CaO and MgO increase with the increase of TOC. In the Tanshanping organic-poor shale, the contents of SiO2, Al2O3, Fe2O3, Na2O, TiO2, and MnO increase with the increase of TOC, while the contents of CaO decrease with the increase of TOC. In particular, TiO2 content is relatively stable, and K2O content increases with the increase of TOC in both organic-rich and organic-poor shale segments of the Tanshanping Formation. The variation trend of the mass fraction of major elements and TOC in the Dacaozi Formation is similar to that of the Tanshanping Formation.

4.3. Distribution of Trace Elements

Because the total number of samples from the Dacaozi Formation is obviously larger than the Tanshanping Formation, the Dacaozi Formation is divided into Dacaozi1 and Dacaozi2 in order to indicate the transition of burial depth. Trace and rare earth element (REE) contents were normalized to North America shale composite (NASC) (for trace elements) and chondrite (for REES) in Figures 4(a) and 4(b) [32, 33]. As shown in Figure 4 and Table 2, the shale samples of the study area were standardized by North American shale (NASC). In the samples of the Tanshanping Formation and Dacaozi Formation of Well Yunningdi-1 (52 samples studied), trace elements Zn (ranges from 0.18 to 2.26) and Sr (ranges from 0.22 to 4.45) are relatively rich, while Rb (ranges from 0.02 to 1.30) is relatively poor (Figure 4(a)). In the two shale members, the trace elements have a similar trend of change, but the relative content of the trace elements in the Tanshanping Formation fluctuates slightly more greatly than the Dacaozi Formation (Figure 4(a)). Compared with the Dacaozi1 Formation and the Dacaozi2 Formation, the contents of trace elements in some samples of the Tanshanping Formation fluctuate more greatly, and the peak values of Zn and Sr are lower, while the relative contents of Rb are lower.

4.4. Distribution of Rare Earth Elements

To begin with, the rare earth elements refer to the lanthanides from La to Lu, which exhibit similar chemical properties. Rare earth elements are good tracers and indicators, which are often used in the analysis of paleosedimentary environment of organic-rich shale [34, 35]. REE content of black shales of the Tanshanping Formation and the Dacaozi Formation in Well Yunningdi-1 in the Ninglang-Yanyuan Basin is shown in Table 3. The REE (∑REE) concentrations of the Tanshanping Formation ranges from 6.00 to 321.24 μg/g, with an average of 198.00 μg/g, the average is significantly higher than the abundance of North American shale REE (173.21 μg/g) and upper crustal abundance (146.40 μg/g) [36, 37]. And the REE (∑REE) concentrations of the Dacaozi Formation ranges from 103.95 to 898.57 μg/g, with an average of 337.47 μg/g; the value is also significantly higher than the North American shale ∑REE, which is nearly twice as much. From the data above, the REE content of the Dacaozi Formation samples is higher than that of the Tanshanping Formation samples, indicating that the shale of the Dacaozi Formation is farther from the provenance area.

The ratio of ∑LREE (light rare earth elements) to ∑HREE (heavy rare earth elements) (), reflecting the REE differentiation degree to some extent. The content of the Tanshanping Formation ranges from 1.12 to 2.29, with an average of 1.85. And the content of the Dacaozi Formation ranges from 1.17 to 2.59, with an average of 2.01, and light rare earth elements are more enriched. However, the of both shale formations is smaller than that of North American shale (7.44), indicating that there was not much distance between these 2 groups and the provenance area.

The degree of differentiation of rare earth elements can be characterized by the slope of the partition curve of rare earth elements to chondrites, reflecting the higher the slope is, the slower the deposition rate is. LaN/YbN represents the slope of rare earth elements in the normalized chondrite diagram (Figure 4(b)), which can be used to evaluate the sedimental rate of shale. Meanwhile, the higher the deposition rate is, the closer LaN/YbN is to 1. After the test, we found that the LaN/YbN of the Tanshanping Formation ranges from 0.85 to 1.37, with an average of 1.13. And the LaN/YbN of the Dacaozi Formation ranges from 0.75 to1.58, with an average of 1.16. It is not difficult to see that the LaN/YbN of the samples of the Tanshanping Formation is closer to 1, leading the deposition rate of shale in the Tanshanping Formation is higher. What is more, the transport distance of the Tanshanping shale is shorter than that of the Dacaozi shale, which can reduce the loss of organic matter effectively.

The geochemical composition of sediments provides important information about the parent rocks from which they derived. The rare earth elements (REEs) are considered to be best indicators of provenance, as they are unlikely to be lost during postdepositional or sedimentary transport processes [38, 39]. In order to eliminate the serrated change of REE abundance with increasing atomic number caused by the even-odd rule of element abundance, REE concentrations were normalized to chondrite in Figure 4(b). Because it is generally accepted that there is no fractionation between LREE and HREE in chondrites, the chondrite standardization enables any separation between REE in the sample to be clearly shown.

The REE distribution pattern of shale samples from the Tanshanping Formation and the Dacaozi Formation in the study area is relatively stable and similar (Figure 4(b)). Figure 4(b) shows that LREE is relatively enriched and HREE is relatively depleted, and the distribution curve has a large slope in the light REE member, while it is flat in the heavy REE member. In addition, Ce and Eu are slightly depleted. On the whole, the REE distribution pattern is right tilted, which means LREE is enriched and HREE is flat. In another aspects, it also reflects the consistency of the provenance of the Tanshanping Formation and Dacaozi Formation in the study area.

5. Discussion

5.1. Paleoredox Condition Control on Deposition of the Tanshanping and Dacaozi Formations

Redox condition was the indispensable factor for organic matter preservation. Considering that particular self-stable and redox-sensitive elements (such as V, Cr, Co, Ni, Sc, and Zr) have been a bit affected or unaffected by diagenetic alternation due to their significantly lower correlation coefficients relative to other element contents [40, 41], the associated redox-sensitive trace elements, especially the trace-element ratios (e.g., V/Ni, Ni/Co, V/Cr, V/Sc, and V/ (V + Ni)) are considered empirical indicators of paleoredox conditions. For the above reasons, the ratios of these elements (e.g., V/Ni, V/Cr, V/Sc, Ni/Co, and V/(V+Ni)) are usually considered effective proxies for distinguishing redox conditions of seawater to eliminate the underlying limitation utilizing a single trace element. Previous studies have suggested that high ratios of these elements reflect reducing conditions.

Previous studies have demonstrated that the redox-sensitive elements were commonly enriched in the anoxic sediments, which could be used as powerful proxies for redox conditions [42]. Redox conditions have an important effect on the preservation of organic matter, the formation of sediments, and the formation and evolution of organisms. Specially, the redox condition indicators including V/Sc, V/Cr, Ni/Co, and V/(V+Ni) were widely used [4345]. Constraints on the paleoredox conditions and paleoenvironment are provided by the above indicators, which show the discrimination of trace elements in different redox environments of shale in Well Yunningdi-1 (Figure 5). As the indexes to quantitatively evaluate the redox environment of the sedimentary waterbody, it is obvious that the sedimentary environment of the Tanshanping and Dacaozi Formations is different, suggesting that the Tanshanping Formation mainly deposited in a relatively anoxic environment.

On the whole, the redox environments changed greatly in the sedimentary history from the Dacaozi Formation to the Tanshanping Formation. The analysis of this study shows that, influenced by the paleoclimate and terrigenous input, the preservation conditions are triggered to change, reflecting on the TOC concentrations in the end, showing reasonable agreement with the decisive effect of the preservation condition of the mechanism of accumulation and preservation of the organic matter.

5.2. Paleoclimate Control on Deposition of the Tanshanping and Dacaozi Formations

Our investigations on CIA and REE further strengthened overall dysoxic to anoxic conditions of the waterbody during accumulation of the Dacaozi Formation to Tanshanping Formation, which provides the foundation for the following interpretations of the paleoclimate effect on redox environments.

The climatic conditions control the temperature, Eh, and pH of the sedimentary environment [46]. Chemical weathering is a crucial process to control the evolution of the earth surface system, by shaping landscapes, supplying nutrients and trace elements from the lithosphere to the biosphere, and regulating global chemical cycles. In the process of sediment formation, several geochemical indexes have been proposed to assess the degree of chemical weathering [4750], among which the chemical index of alteration (CIA) is one of the most frequently used indexes. The chemical index of alteration (CIA) is calculated as the ratio in molecular proportions, where CaO represents Ca content in silicate minerals only. It is worth noting that the CIA has been widely used to quantify the degree of chemical weathering for (paleo-)weathering profiles and to reconstruct source region chemical weathering conditions for ancient and modern sediments [51]. In addition, the climate indexes (CIA) display the slowly warming environment, which is often used to reflect the paleoclimate conditions. It means the higher CIA values experienced relatively intense chemical weathering, which indicates a warm and humid climate [52]. It is also consistent with previous studies on the relationship between the CIA and paleoclimate which provides a better standard. As discussed in previous studies, when the value of CIA ranges from 50 to 65, it indicates arid climate. When the value of CIA ranges from 65 to 85, it indicates semiarid climate. When the value of CIA ranges from 85 to 100, it indicates warm and humid climate [53, 54]. Specially, the shale-dominated the Dacaozi Formation and Tanshanping Formation experienced a moderate to weak degree of chemical weathering as evidenced by the significant decrease of CIA (mean 58.5 and 43.2, respectively) (Figure 6 and Table 1). Besides, the can also indicate climate conditions [55]. It is noticed that the climate became more dry-cold from the Dacaozi Formation to the Tanshanping Formation as evidenced by the significant decrease of ∑REE (mean 337.5 and 198.0, respectively) (Figure 6 and Table 1).

Therefore, these considerable variations likely reflect the impact of alternating climatic conditions of warm-humid and cool-arid, which confirm a relatively intensive weathering and its associated warm-humid setting that occur during the deposition of shales of the Dacaozi Formation.

In summary, there are several lines of evidence that it is more conducive to form a higher paleoproductivity in the context of a relatively warm-humid climate which promote biological growth in the sedimentary waterbody. From the above discussion of the values (CIA, ∑REE) from the Dacaozi Formation to the Tanshanping Formation, we have found that although the climate is relatively warm-humid during the deposition of the Dacaozi Formation, the weak preservation condition leads to destroy the accumulation of organic matter, which shows the lower TOC concentrations in the end.

5.3. Paleoproductivity Control on Deposition of the Tanshanping and Dacaozi Formations

In particular, the initial productivity of waterbody can be characterized by the nutrient richness and organic carbon flux of surface water, and the change of initial productivity plays a key role in the accumulation of organic matter.

The Ba content is an indicator widely used to reflect the paleoocean bioproductivity [5660]. There is a high concentration of SO42- ions reoxidized by H2S on the surface of decayed organic matter. These ions will react with the Ba2+ in seawater to form BaSO4 and result in sedimentation. Therefore, the area with a higher bioproductivity also has a higher content of BaSO4 [61, 62]. The trace elements in rocks are composed of terrigenous input constituents and authigenic ones. Only the latter can reflect the characteristics of the paleosedimentary environment [6365]. The Ba produced under the biological action is known as excess Ba (Baxs). The Late Archean Australian Shale standard is commonly used to calculate excess barium (Baxs) to characterize paleoproductivity. It is worth noticing that the higher values of Baxs are indicators of higher paleoproductivity in the sedimentary waterbody of the Dacaozi Formation (Figure 6).

Furthermore, phosphorus (P) is considered to be one of the ultimately limited elements in marine environment, and also, it belongs to the most important nutrient elements for plankton [66] and a major component of skeletal material. Therefore, P is widely used as an indicator of paleoproductivity. Organic matter and authigenic minerals may have a dilutive effect on the absolute P content in the terrigenous detrital matter. To eliminate the effect, the P/Ti ratio is used to evaluate the paleoproductivity rather than the absolute P content [67, 68] because Ti is usually derived from the terrigenous detrital matter [69, 70]. The higher ratios of P/Ti of the Dacaozi Formation (which is higher than the P/Ti ratio of the black flint in the Ubara profile with medium productivity (with an average of 0.34)) reported in Figure 6 can also imply relatively high paleoproductivity throughout the depositional sequence of shales, which coincides with the better provenance supply as argued above in a relatively warm-humid sedimentary setting, which provided favorable conditions for the accumulation of organic matter.

Since the depth of water controls the redox conditions and Ce depletion degree of water, Ce anomaly can also be used to indicate the relative change of the sea level. With the rise of sea level, the oxygen content of bottom water and the δCe of sediments are reduced [71].

From the above discussion, we can get a stronger idea that detailed analysis of the paleoredox conditions and further investigation of the paleoclimatic and paleoproductivity forcing mechanisms operating during sedimentation indicate that the Dacaozi organic-rich shales was deposited in the more warm-humid environments. Climatic fluctuations resulted in the difference degree of the terrigenous input of the Dacaozi and Tanshanping Formations, which coincide well with the water paleoproductivity changes and redox environmental evolution [72]. Specifically, during the deposition of the Dacaozi Formation, significant increases in the temperature and humidity led to the relatively high terrigenous inputs with more nutrient inputs. The terrigenous inputs promoted abundant planktons and favoring organic matter accumulation [73]. In addition, the shallow waterbody made it easier for sunlight to penetrate the waterbody, which coincides with the relatively high paleoproductivity. However, the preservation condition of the organic matter in the Dacaozi Formation was generally poor, as often the relative sea level was so low that it trend to accumulate more oxygen [74]. Specially, the higher terrigenous input diluted more organic matter, and the longer distance of transportation and lower deposition rate have accelerated the destruction of preservation conditions. Due to the worse preservation conditions, the TOC concentrations of the Dacaozi Formation were lower. On the contrary, we have found that the paleoproductivity during the sedimentary period of the Tanshanping Formation was lower than that of the Dacaozi Formation, reflecting the relative cold-dry climate and lower terrigenous input. It is obvious that the paleoproductivity in the waterbody was worse, and that is mainly because it did not have enough terrigenous input, reflecting the limited provenance supply. In another aspects, we found that the range of the penetration from sunlight was relatively narrow, for the sea level was relatively high. Nevertheless, the higher sea level supports a reductive and stable bottom waterbody (Figure 7). And the lower terrigenous input indicates the decrease of fresh water input, leading to be more reductive in the waterbody. In addition, it is worth to emphasis on the shorter distance of transportation during the sedimentary period of the Tanshanping Formation, which enhance the preservation ability in the sedimentary center.

6. Mechanism of Accumulation and Preservation of Organic Matter

The study area was in the sedimentary environment of the carbonate platform. The high biological productivity, anoxic, and reducing environment were conducive to the accumulation and preservation of organic matter in shale [75]. In order to reconstruct the model of accumulation and preservation of organic matter in the carbonate platform facies, it is essential to clear up the relationship between TOC and relative sea level, paleoproductivity, and redox conditions of the unique sedimentary facies which has many unconventional features. Therefore, this study preliminarily discusses the accumulation and preservation mechanism of shale organic matter in Well Yunningdi-1 in the Ninglang-Yanyuan area by using the evolution of water redox conditions combined with the research results of sea level change, paleoproductivity, paleoclimate, and terrigenous input.

6.1. Relationship between TOC and Relative Sea Level Variations

In considering sea level variations during the sedimentary period, the relative sea level and regional paleoclimate are suggested here as two distinct factors for the influence of accumulation and preservation of organic matter at different scales.

On the one hand, the difference of the redox environment of the Tanshanping and Dacaozi Formations can be analyzed from the relative sea level. Previous studies have shown that water depth is one of the elements which controls the redox condition and Ce depletion degree. With the rise of sea level, both of the oxygen content of bottom water and the δCe of sediments decrease [71]. The value of δCe in the Tanshanping Formation is lower than that in the Dacaozi Formation (mean 0.763 and 0.824, respectively), reflecting the rise of the relative sea level. The research of the redox environment emphasized the importance of the relative sea level response in the process of deposition of the carbonate platform and throw light on whether the relative sea level variations control the redox environment of the sedimentary waterbody in the carbonate platform.

On the other hand, the anoxic environment in the sedimentary waterbody of the Tanshanping Formation is ascribed to the regional paleoclimate effect during the sedimentary period. At the same time, the terrigenous inputs promote planktons’ propagation and growth.

This implies that the combined effects of relative sea level variations and regional paleoclimate driving mechanisms favored the accumulation and preservation of organic matter in the study area. To summarize, the combined impact of these factors is ultimately reflected in the TOC.

6.2. Relationship between TOC and Paleoproductivity

Based on the above researches, Baxs indicates that the paleoproductivity of the sedimentary waterbody in the Dacaozi Formation is higher than that of the Tanshanping Formation (mean 856.29 and 254.49, respectively), which are shown in Figure 6 and Table 2. Considering that the accumulation and preservation of the organic matter in the Dacaozi Formation is largely controlled by paleoproductivity due to their significantly higher correlation coefficients between TOC and paleoproductivity indexes (e.g., Pb and Ni/Al) [76], which are shown in Figures 7(d) and 8(e). In addition, the associated indexes of climate and terrigenous input (e.g., REE and CIA values) are considered important indicators of paleoproductivity. However, due to the damage of the preservation condition, there is more oxygen enrichment of the sedimentary waterbody of the Dacaozi Formation, which led to the lower TOC concentrations of the Dacaozi Formation.

6.3. Relationship between TOC and Redox Conditions and Terrigenous Inputs

The redox conditions are crucial in the mechanism of accumulation and preservation of the organic matter, which shown in Figure 7. As most scholars believe, this is in close agreement that the redox environment of sedimentary waterbody is not only affected by redox indexes but also affected by terrigenous input.

On the one hand, we have supported the idea that the accumulation and preservation of the organic matter in the Tanshanping Formation are largely controlled by paleoredox due to their significantly higher correlation coefficients between TOC and paleoredox indexes (e.g., Cr, Cu/Zn, and δU), which have been considered reliable redox indexes[77, 78], reflecting in Figures 8(a)8(c).

On the other hand, terrigenous input is also considered an important element, coinciding reasonably well with the paleoproductivity and paleoredox environment.

From the above discussion, we can get a stronger idea that CIA and ∑REE can reflect paleoclimate and terrigenous input together. The CIA and ∑REE values of the Dacaozi Formation are higher, reflecting the relatively warm-humid paleoclimate and higher terrigenous input. At the same time, the higher terrigenous input diluted more organic matter. According to the previous studies, both of the paleoclimate and terrigenous input led to oxidation of the sedimentary waterbody of the Dacaozi Formation. On the contrary, it was relatively dry-cold during the sedimentary period of the Tanshanping Formation, which reflects lower terrigenous input. The lower terrigenous input has effectively decreased the dilution of organic matter, leading a reductive environment of sedimentary waterbody in the Tanshanping Formation.

6.4. Model of Accumulation and Preservation of Organic Matter
6.4.1. The Interactions between the Relative Sea Level, Paleoproductivity, Redox Condition, and Carbonate Platform Sedimentation

As one of the key research points of the Ninglang-Yanyuan Basin, the Tanshanping and Dacaozi Formations of Well Yunningdi-1 have the following features.

According to previous studies, we have done a lot of works about the marine facies and pure continental facies. Specially, the sedimentary position of the study area is different from the previous researches, which belongs to the carbonate platform facies (Figure 7).

Unlike marine deposits, sedimentary water bodies of the carbonate platform are partially connected with the open sea, showing reasonable agreement that the accumulation and preservation of organic matter will be partially affected by the open sea (Figure 7). However, the connection area with the open sea is really limited in the sedimentary waterbody here, it will be greatly affected by other multiple factors (such as paleoclimate, paleoproductivity, terrigenous input, and redox conditions) within the sedimentary waterbody.

Aiming at the mechanism of accumulation and preservation of the carbonate platform, this paper takes the Tanshanping Formation and Dacaozi Formation as an excellent example and focuses on the analysis of the dialectical relation between the paleoproductivity and paleoredox environment and the accumulation and preservation model of shale in the study area.

In terms of paleoproductivity and redox condition, it is obvious to know that both of them affect the accumulation and preservation of organic matter. In the study on the accumulation and preservation mechanism of organic matter in the Tanshanping and Dacaozi Formations, the comprehensive results of the combined effects of paleoproductivity and redox conditions form a marked contrast.

From the above discussion, we have summarized a view. Based on the example of the Dacaozi Formation, the TOC of the study shale member is not necessarily higher in the context of higher productivity waterbody. The reason why it shows lower TOC concentrations is that the sedimentary waterbody of the Dacaozi Formation had more oxygen because of the fall of the sea level and the dilution from the terrigenous input. In another aspect, based on the example of the Tanshanping Formation, the paleoproductivity was lower, for the relatively dry-cold climate has blocked biological growth. As a complement, the range of the penetration from sunlight was relatively narrow because of the rise of the sea level. However, the waterbody was more reductive during this depositional period, because the lower terrigenous input which was reflected by the dry-cold climate indicates the decrease of fresh water input. All of the conditions guarantee a better preservation condition of organic matter. On the whole, the TOC concentrations of the study shale are not necessarily lower although the paleoproductivity is relatively poor.

6.4.2. Depositional Model for Organic-Rich Shale of Carbonate Platform

In terms of sedimentary facies, the Dacaozi Formation developed saline lagoon facies, which belongs to shallow water deposits farther from the sedimentary center. During the deposition period of the Dacaozi Formation, the water energy was low, which was mainly due to weak tidal action. It is noticed that only a limited number of organisms developed in the waterbody because of the abnormal salinity of the waterbody. Based on the above factors, it can be seen that the organic matter enrichment capacity of the Dacaozi Formation was relatively poor. As for the sedimentary period during the Tanshanping Formation, it developed the reef deposit with rich reservoir performance. It developed closer to the sedimentary center. Due to the rich nutrients and clean water, there were a lot of organisms growing in the waterbody during the sedimentary period of the Tanshanping Formation.

According to the detailed analysis of the paleoredox conditions and further investigation of the multifaceted studies about relative sea level, paleoclimate, terrigenous input, and paleoproductivity during sedimentation, we have purposed the conceptual model illustrating the influences of the relative sea level variations, paleoredox condition, paleoproductivity condition, and paleoclimate condition on organic matter accumulation and preservation for the Tanshanping Formation and the Dacaozi Formation.

Detailed analysis of the paleoredox conditions and further investigation of the paleoproductivity forcing mechanisms operating during sedimentation indicate that the shales of the Dacaozi and Tanshanping Formations were deposited under the controlling of paleoredox conditions. During the deposition period of the Dacaozi Formation, which was deposited in the more warm-humid shallow water environments, we found that significant increases in the temperature and humidity led to more fresh water inputs and a high degree of chemical weathering [7983]. The higher terrigenous inputs, promoting abundant planktons and better condition of organic matter accumulation, coincides with the higher level of penetration from sunlight. From the above factors, it is obvious to exhibit a higher paleoproductivity in the sedimentary waterbody of the Dacaozi Formation. However, the preservation condition of the organic matter in the Dacaozi Formation was poor, as often the relative sea level was so low that it tend to accumulate more oxygen. Specially, the higher terrigenous input diluted more organic matter, and the longer distance of transportation and lower deposition rate have accelerated the destruction of preservation conditions. As we can see in Figure 7(b), there were a few dead planktons left in the bottom water, in which their bodies were damaged to varying degrees because of the poor preservation condition. Due to the worse preservation conditions, the TOC concentrations of the Dacaozi Formation were lower. During our investigation of the mechanism of accumulation and preservation of organic matter in the Tanshanping Formation, the depositional model for organic-rich shale of the Tanshanping Formation has been discovered. Our previous work has summarized that the paleoproductivity during the sedimentary period of the Tanshanping Formation was lower than that of the Dacaozi Formation, which reflects the relative cold-dry climate and lower terrigenous input. As the relative sea level rise, increased water depth would decrease the mixing chance of surface oxygenated water into sediment-water interface, and increased sea level also trapped clastic flux inshore, which reflects lower terrigenous input [84]. In addition, the range of the penetration from sunlight was relatively narrow, for the sea level was so high that it blocks sunlight from entering into the water. Nevertheless, it supports a reductive and stable bottom water because of the rise of the sea level (Figure 7(a)), showing reasonable agreement with the influence of the decrease of fresh water input which is attributed to the fall of terrigenous input. Moreover, the shorter distance of transportation of terrigenous detritus during the sedimentary period of the Tanshanping Formation enhances the preservation ability in the sedimentary center. As is shown in Figure 7(a), the number of plankton carcasses greatly increases due to the lack of oxygen in the waterbody or catastrophic events. Decomposition of these carcasses consumes large amounts of oxygen and form H2S gas, which further ensures the anoxic condition of the bottom water [85]. It is worth noting that there are a large number of plankton carcasses in the bottom water, which were hardly destroyed by external forces because of the stable anoxic waterbody and the shorter distance of transportation of terrigenous detritus.

To sum up, the significance of establishing the accumulation and preservation model of organic matter lies in discussing the coupling interactions between the paleoproductive forces and paleoredox environment. Because of the unique features of the carbonate platform facies, the examples of the Dacaozi Formation and the Tanshanping Formation are used to illustrate: the TOC concentrations of the study shale member are not necessarily higher although the paleoproductivity is relatively high, and the TOC concentrations of the study shale is not necessarily lower although the paleoproductivity is relatively poor. The considerable fluctuations in the TOC concentrations are reasonably attributed to the redox environment of the Tanshanping and Dacaozi Formations, which determines the preservation conditions.

These findings highlight the importance of coupling between the paleoproductivity and redox environment in the more specific carbonate platforms in further research.

7. Conclusions

On the basis of our results from TOC concentrations, major, trace, and rare earth elemental fingerprints and comprehensive analyses on the special features of the carbonate platform facies in the Tanshanping and Dacaozi Formations, this current study provides a good opportunity for documenting the major controlling factors of the accumulation and preservation of organic matter and for revealing new insights into the response of the coupling interactions between the paleoproductivity and redox environment. The achievements are thus conducted to generate a comprehensive and specific depositional model for the carbonate platform facies as evident in the case of the Tanshanping and Dacaozi Formations, Ninglang-Yanyuan Basin, Southwest China, and the following conclusions have been drawn: (1)The variation trend of the content of major elements and TOC is similar to the study shales of the Tanshanping and Dacaozi Formations. Compared to the Dacaozi Formation, the value of the trace elements (Ni/Co, V/Cr, and V/Sc) in our study confirms that the shale of the Tanshanping Formation was deposited in an reductive environment. As for the rare earth element, the shale members of the Middle Devonian Tanshanping and Dacaozi Formations are both enriched in LREE and depleted in HREE, with negative anomalies of Ce and Eu(2)Specially, the sedimentary position of the study area is different from the previous researches, which belongs to the carbonate platform facies. We have concluded that the main controlling factors are the paleoclimate, paleoproductivity, terrigenous input, and paleoredox conditions in the sedimentary waterbody. Aiming at the mechanism of accumulation and preservation of the carbonate platform, this paper takes the Tanshanping Formation and Dacaozi Formation as an excellent example and focuses on the analysis of the coupling interactions between the paleoproductivity and paleoredox environment. On the whole, these above features show its uniqueness, which deserved to be analyzed(3)Dacaozi organic-rich shales were deposited in the more warm-humid environments. In addition, the shallow waterbody made it easier for sunlight to penetrate the waterbody, which coincides with the relatively high paleoproductivity. However, the preservation condition of the organic matter was generally poor, as often the relative sea level was so shallow that it tend to accumulate more oxygen. Specially, the higher terrigenous input diluted more organic matter, leading to the destruction of preservation conditions. Due to the worse preservation conditions, the TOC concentrations of the Dacaozi Formation were lower. On the contrary, the paleoproductivity during the sedimentary period of the Tanshanping Formation was lower. Nevertheless, the higher sea level and the lower terrigenous input reduced input of fresh water and supported a reductive bottom waterbody. Through a combination of the above factors, the TOC concentrations of the Tanshanping Formation shale were relatively high(4)The significance of constructing the accumulation and preservation model of organic matter lies in discussing the coupling interactions between paleoproductive forces and redox environment. The examples of the Dacaozi and Tanshanping Formations are used to illustrate the following: the TOC concentrations of the study shale member are not necessarily higher although the paleoproductivity is relatively high. In the same way, the TOC concentrations of the study shale are not necessarily lower although the paleoproductivity is relatively poor. As a result, the preservation conditions have been the decisive factor in the accumulation and preservation model of organic matter of the carbonate platform

Data Availability

Some of the data are contained in a published source cited in the references. All the data in this article is accessible to the readers.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

A critical review and constructive comments from Dr. Rusi Zuo undoubtedly enhanced the quality of interpretations of chemostratigraphic data and are highly appreciated. Many thanks are given to Dr. Liang Xu for his sincere help during the analysis of the major and trace elements. We are especially grateful to Dr. Rong Chen of Oil & Gas Resource Survey, China Geological Survey, for the help during TOC content analysis. This work was supported by the Science Foundation for Distinguished Young Scholars of China University of Petroleum, Beijing (No. 2462020QNXZ004) and the National Natural Science and Technology Major Project (No. 2016ZX05034-001 and 2017ZX05035-002).