Abstract

Organic matter is the material basis of shale hydrocarbon generation. The current organic matter content in shale is controlled by the original sedimentary organic matter abundance. Therefore, the study of the enrichment mechanism of sedimentary organic matter in shale has become an important issue to be solved. The Upper Yangtze area is the important exploration and exploitation area of marine shale gas in China. The shale of the Upper Ordovician Wufeng Formation-Lower Silurian Longmaxi Formation in the Yangtze area is the research object. Choosing redox indicator and biological productivity indicator, the study explores the enrichment mechanism of sedimentary organic matter from two aspects, sealing of water and volcanic activity. The results show that excess siliceous mineral in the shale of the Wufeng Formation-Longmaxi Formation in the Upper Yangtze area is bioorigin. Excess siliceous mineral can be used as one of the indicators of biological productivity. On the one hand, layer phenomenon occurred since the strong water sealing during the sedimentary period of Wufeng and the lower section of the Longmaxi Formation, which results in the high content of oxygen in surface water. On the other hand, the active volcanic activity brought volcanic ash which was beneficial to biological reproduction. Both of these factors led to higher biological productivity during this period. At the same time, the strong sealing of water made the lower layer of the water more reductive, and the active volcanic activity caused climate change, enhancing the reduction of the lower layer of the water, which made the rich organic matter deposited from the surface water well preserved. In the sedimentary period of the upper section of the Longmaxi Formation 1st member in the Upper Yangtze area, on the one hand, due to the weakened sealing of water, the oxygen content of the upper water decreased. On the other hand, the volcanic activity weakened until it stopped, and the source of volcanic ash rich in nutrient elements decreased. These two aspects led to lower biological productivity during this period. At the same time, the weaker water sealing could lead to a decrease in the reduction of the lower layer of the water, and the gradual cessation of volcanic activity no longer affected the climate, causing the destruction of sedimentary organic matter by oxidation.

1. Introduction

With the change of exploration thinking and the advancement of hydraulic fracturing technology, shale gas exploration has achieved great success in North America such as Permian Basin, Fort Worth Basin, Appalachian Basin, and Williston Basin [16]. China also has huge potential for shale gas, especially the Upper Ordovician Wufeng Formation-Lower Silurian Longmaxi Formation 1st member in the Sichuan Basin of the Upper Yangtze area. After the oil companies, such as PetroChina and Sinopec, explored the marine shale gas in this stratum, shale gas fields such as Weiyuan, Changning, Zhaotong, Fushun-Yongchuan, Dingshan, and Jiaoshiba have been built since 2010, and good industrial production capacity has been obtained [710]. Predecessors did researches on shale reservoir characteristics [1117], shale gas accumulation mechanism [1822], and shale gas preservation conditions [2328] and have yielded a series of scientific research results. Organic matter is the material basis for hydrocarbon generation in shale and is also the main reservoir space and seepage channel for shale gas. The sedimentary environment affects the enrichment of organic matter. So, the analysis of the geological factors affecting the abundance of sedimentary organic matter and the study of the enrichment mechanism of sedimentary organic matter in shale have become the important issues to be solved.

The enrichment mechanism of deposited organic matter in shale has been studied by some researchers. Xia et al. [29] studied the Lower Cambrian Niutitang Formation in the southeastern Chongqing in periphery of Sichuan Basin in the Upper Yangtze area and concluded that the different strata of the Niutitang Formation were affected by the deep hydrothermal fluid brought by the rising ocean current. When the environmental condition reached the reducing condition which was favorable for the preservation of organic matter, the degree of reduction was no longer the main factor determining the degree of organic matter enrichment. Besides, the rising ocean current promoted the enrichment of organic matter and the change of paleoaltitude had little effect on the enrichment of organic matter. Zhang et al. [30] studied the sedimentary environment of the Lower Cambrian shale rich in organic matter in the Upper Yangtze area and concluded that the organic shale in the shelf was mainly controlled by seawater redox conditions, and the organic matter was more easily preserved in the reducing environment. In addition, the organic matter content in the organic-rich shale in the slope-basin area was mainly controlled by the biological primary productivity and sedimentary environment. Under the same sedimentary conditions, the organic matter abundance was higher than in the shelf. Qiu et al. [31], after studying the sedimentary environment of the Wufeng Formation-Longmaxi Formation shale in the Wuxi Block of the Upper Yangtze area, believed that the eustacy controlled the degree of hypoxia in the sedimentary water. The sea level decline made the bottom water rich in oxygen, while the sea level rise made the bottom water lack oxygen, thus controlling the degree of organic matter enrichment in shale.

The large-scale exploration of shale gas in recent years has provided new data for the study of the mechanism of deposited organic matter. Predecessors did researches on the source of siliceous mineral. Holdaway and Clayton [32] defined the concept of excess siliceous mineral, which refers to siliceous minerals beyond the source of normal terrigenous debris, and proposed a quantitative calculation of excess siliceous mineral. Wedepohl [33], Adachi et al. [34], and Yamamoto [35] proposed a method for determining whether siliceous mineral is of hydrothermal origin or bioorigin by using the Al-Fe-Mn ternary diagram. Li et al. [36], Shu et al. [37], and Wu et al. [38], after giving detailed description of the shale core of the Upper Ordovician Wufeng Formation-Lower Silurian Longmaxi Formation, found bentonite in some layers, which was derived from volcanic activity. This paper firstly determines the redox indicator and biological productivity indicator, and then discusses the degree of water sealing and the influence of volcanic activity on sedimentary organic matter enrichment, and establishes the Late Ordovician-Early Silurian sedimentary organic matter enrichment model in the Upper Yangtze area.

2. Geological Settings

2.1. Sedimentary and Stratum Characteristics

According to previous studies [3941], during the period of Cambrian to Silurian, the Yangtze Plate and the Cathaysian Plate gradually collided. As shown in Figure 1, the Upper Yangtze area has a wide area. In the Late Ordovician-Early Silurian, after being squeezed by the Cathaysian Plate, it formed the intracratonic depression basin. This paper takes a set of strata which is widely deposited in the Yangtze Plate during the Late Ordovician-Early Silurian as research layer. Because the Yangtze Plate is large, it has different names in different areas. In the Upper Yangtze area, the Late Ordovician sedimentary strata are called as the Wufeng Formation, and the Early Silurian sedimentary strata are called as the Longmaxi Formation. The target layer of the study is the Wufeng Formation-Longmaxi Formation 1st member, and its lithology is bipartite. The Wufeng Formation-the lower section of the Longmaxi Formation 1st member is mainly black siliceous organic-rich shale and the upper section of the Longmaxi Formation 1st member is a combination of dark grey shale, silty shale, and siltstone.

2.2. Tectonic Characteristics

According to previous studies [4648], the primitive continental crust of southern China was separated into two ancient plates, Yangtze Plate and Cathaysian Plate, in the Early Mesoproterozoic. In the Early Cambrian, the two plates were in tension and large-scale transgression occurred, depositing a set of organic-rich shale almost throughout the plate. Afterwards, the water gradually became shallower, and the lithology gradually changed from fine shale and silty shale to coarse-grained clastic rock such as siltstone and sandstone. In the Ordovician, being affected by the compression collision of the Cathaysian Plate, the water continued to become shallower and lithology was transformed from a clastic rock sedimentary system to a carbonate rock sedimentary system. A large-scale transgression occurred in the Late Ordovician-Early Silurian, lithology was changed into a clastic sedimentary system, and a set of organic-rich shale was deposited in a deep shelf surrounded by ancient land. During the period of Cambrian to Silurian, the Cathaysian Plate gradually subducted and collided with the Yangtze plate. At the end of the Silurian, the Yangtze Plate and the Cathaysian Plate were integrated into one and became the unified South China Plate [49].

3. Samples, Experiments, and Source of Data

Because of the large sedimentary range of the deep shelf, the small variation of the lithofacies, and the high homogeneity of marine shale, in this study, Jiaoye-1 well, the typical well in the Jiaoshiba Block of the Sichuan Basin, was selected to analyze the enrichment mechanism of the Late Ordovician-Early Silurian sedimentary organic matter in the Upper Yangtze area. Fifteen debris samples were taken at intervals of 5~6 m from the Upper Ordovician Wufeng Formation-Lower Silurian Longmaxi Formation shale in the Jiaoye-1 well and Al, Fe, and Mn were analyzed by Axios-MAX X-ray Fluorescence. Forty-five core samples were taken from the Wufeng Formation-Longmaxi Formation in the Jiaoye-1 well, and TOC content and whole rock mineral composition were experimented by OG-2000V total organic carbon analyzer and Ultima IV X-ray diffraction full rock analyzer.

In this paper, 20 samples were taken from the bentonite in the shale of the Jiaoye-1 well. The main element and trace element were analyzed by PE Elan 6000 inductively coupled plasma mass spectrometer (ICP-MS). Some of the data were referenced in the data of Wang et al. [50], Luo et al. [51], Hu et al. [52], and Tang et al. [53]. Besides, 2 samples from the Wufeng Formation and 10 samples from the Longmaxi Formation 1st member were selected to measure content of molybdenum by PE Elan 6000 inductively coupled plasma mass spectrometer (ICP-MS) and excess molybdenum (MoXS) was calculated. This paper collected the calculation data of excess molybdenum (MoXS) of the biological productivity of the Wufeng Formation-Longmaxi Formation in the Jiaoye-1 well in Guo et al. [54] and collected the logging data of Si, Al, U, and Th of the Jiaoye-1 well provided by Schlumberger.

4. Results and Discussion

4.1. Indicator of Redox Environment and Biological Productivity
4.1.1. Indicator of Redox Environment

Determining sedimentary environments by using elemental geochemical indicators is a common method. By analysis of the whole rock, Jones and Manning [55] had a proposal that U/Th ratio can reflect the sedimentary redox conditions. Generally speaking, reflects an anoxic environment, U/Th ratio between 0.75 and 1.25 reflects a dysoxic environment, and reflects an oxidizing environment. Chen et al. [56], Zhao et al. [57], and Wang et al. [58] also used this ratio to determine the redox environment of south of the Sichuan Basin. This paper used the same method. Using this ratio and classification, results are as follows. As is shown in Table 1, U/Th ratios of the Wufeng Formation (2411 m~2415 m) are mostly 0.75~1.25, meaning an oxygen deficiency oxidation environment. The U/Th ratios of the bottom of the lower section of the Longmaxi Formation 1st member (2406 m~2411 m) are principally more than 1.25, with a maximum of 3.73, showing an anaerobic environment. Most of the U/Th ratios of the lower section of the Longmaxi Formation 1st member (2386 m~2406 m) are 0.75~1.25, suggesting an oxygen deficiency oxidation environment. Most of the U/Th ratios of the upper section of the Longmaxi Formation 1st member (2326 m~2386 m) are less than 0.75, with a minimum of 0.17, reflecting an oxidation environment.

4.1.2. Indicator of Biological Productivity

(1) Excess Molybdenum. The content of molybdenum is a widely used indicator reflecting the productivity of paleoocean organism. Molybdenum has two origins, terrestrial detrital and biological origin. Mo, which is only derived from biological action, is called excess molybdenum (MoXS) and content can be calculated by the following formula.

Mos is the content of molybdenum in the sample, and Als is the content of aluminum in the sample. The (Mo/Al)PAAS value is 0.00001368, which is the average content in Post-Archaean Australian shale.

Besides, 7 data points are from Guo et al. [54]. And all MoXS values are as follows (Table 2).

(2) Siliceous Mineral from Bioorigin. The sources of siliceous mineral can be divided into 3 types, terrestrial detrital origin formed in normal condition, hydrothermal origin, and bioorigin formed in special condition [5963]. Excess siliceous mineral (abbreviation: Siex) refers to siliceous mineral excepting terrestrial detrital origin, and content can be calculated by the following formula.

Sis is the content of silicon in the sample, and Als is the content of aluminum in the sample. The (Si/Al)bg value is 3.11, which is the average content in shale [32].

The excess siliceous mineral content of the Upper Ordovician Wufeng Formation-Lower Silurian Longmaxi Formation 1st member in the Jiaoye-1 well was calculated, and the result is shown in Figure 2. The siliceous mineral from the upper section of the Longmaxi Formation 1st member is basically of terrigenous detrital origin, while the siliceous mineral from the Wufeng Formation-the lower section of the Longmaxi Formation 1st member contains excess siliceous mineral. By analysis of the mineral composition of shale, the siliceous mineral content from the Wufeng Formation-the lower section of the Longmaxi Formation 1st member is higher, 40%~70%. And the siliceous mineral content from the upper section of the Longmaxi Formation 1st member gradually decreases to 20%~35%. In the intervals with excess siliceous mineral of the Wufeng Formation-the lower section of the Longmaxi Formation, the excess siliceous mineral content in more than half of the intervals is 0%~5%, and the excess siliceous mineral content in a part of the intervals is 5%~15%, and the maximal content is 15%~20%.

Wedepohl [33], Adachi et al. [34], and Yamamoto [35] put forward a method for determining whether siliceous minerals are of hydrothermal origin or bioorigin by using the Al-Fe-Mn ternary diagram. In this paper, the measured values of Al, Fe, and Mn element of the Wufeng Formation-Longmaxi Formation 1st member with excess siliceous mineral of the Jiaoye-1 well were placed on the ternary diagram. As shown in Figure 3, it can be found that the values basically fell in the bioorigin area, indicating that the excess siliceous mineral is of bioorigin. According to this, the sources of excess siliceous mineral in the vertical dimension are accurately described in Figure 2. The siliceous mineral derived from bioorigin is from plankton, such as diatoms [6367]. As a result, similar as excess molybdenum (MoXS), the siliceous mineral derived from bioorigin can also reflect biological productivity.

4.2. Analysis of Water Sealing and Its Effect on Sedimentary Organic Matter Enrichment

The degree of sealing of water can be reflected by the relationship between Mo and TOC content [30, 40, 68, 69]. Based on the base map from Algeo and Lyons [68] and Zhang et al. [30], 11 data points of the upper section of the Longmaxi Formation 1st member and 9 data points of the Wufeng Formation-the lower section of the Longmaxi Formation 1st member of the Jiaoye-1 well are placed on the diagram, respectively. Figure 4 illustrates that Mo and TOC content are mainly in strong restriction. The distribution range of the 11 data points from the upper section of the Longmaxi Formation 1st member indicates that its sealing degree is similar to today’s Framvaren (Figure 4(a)). The distribution range of 9 data points from the Wufeng Formation-the lower section of the Longmaxi Formation 1st member indicates that the sealing degree is similar to that of the Black Sea whose sealing is stronger nowadays (Figure 4(b)). This shows that the sealing of water during the Late Ordovician-Early Silurian is stronger, and the sealing during the deposition of the Wufeng Formation-the lower section of the Longmaxi Formation 1st member is stronger than that during the deposition of the upper section of the Longmaxi Formation 1st member.

The analysis results of sealing are consistent with the tectonic background of the Late Ordovician-Early Silurian. Figure 1 demonstrates that during the deposition of the Wufeng Formation-the lower section of the Longmaxi Formation 1st member, there was a strong collision between the Yangtze Plate and the Cathaysian Plate, and shale rich in organic matter was deposited in the deep shelf limited by ancient land, poor connection with the open sea, leading to strong sealing; while during the deposition of the upper Longmaxi Formation 1st member, the collision and compression between the Yangtze Plate and the Cathaysian Plate weakened, and the water became shallower and sealing became weaker.

The sealing of water is one of the major reasons for the difference in redox environment and biological productivity between the Wufeng Formation-the lower section of the Longmaxi Formation and the upper section of the Longmaxi Formation. Due to the strong sealing of water during the deposition of the Wufeng Formation-the lower section of the Longmaxi Formation, layer phenomenon occurred, making the upper water rich in oxygen and the lower water lack in oxygen. The abundant oxygen in the upper water supplied the growth, living, and reproduction of plankton, producing rich source of organic matter, while the anoxic lower water kept sedimentary organic matter from decomposition, which was beneficial for the preservation of organic matter [7072]. During the deposition of the Longmaxi Formation 1st member, the sealing of water weakened, the anoxic lower water gradually mixed with warm upper water rich in oxygen, bringing that the lower water became warmer, and reducing environment was damaged and sedimentary organic matter was diluted and oxidative decomposed, while the oxygen content of the upper water decreased, the upper water became colder, and biological productivity declined. These two aspects made organic matter abundance descend.

4.3. Volcanic Activity Analysis and Its Effect on the Enrichment of Sedimentary Organic Matter

Bentonite occurs at the Wufeng Formation-Longmaxi Formation 1st member. As the product of sedimentation and alteration of tuffaceous materials, bentonite reveals information of volcanic activity. The thickness of bentonite can be used as direct evidence of volcanic eruption interval and reflect the activity frequency of volcanic eruption [37, 38, 73, 74]. As shown in Figure 2, according to the density of bentonite development, the Wufeng Formation-Longmaxi Formation is divided to 3 sections. The Wufeng Formation-the lower section of the Longmaxi Formation is the dense bentonite layer, representing frequent volcanic activity. The middle section of the Longmaxi Formation is the rare bentonite layer, reflecting weakened volcanic activity. The upper section of the Longmaxi Formation is no bentonite layer, showing the stop of volcanic activity. The major elements and trace elements of 20 bentonite samples (Table 3) are put on the discrimination diagram, and the base map is from Winchester and Floyd [75]. As can be seen from Figure 5, the volcanic materials forming bentonite are mainly trachyandesite, rhyodacite dacite, and andesite, all of which are intermediate-acidic magmatic rocks originating from deep in the crust.

Volcanic activity is another reason for the difference in redox environment and biological productivity between the Wufeng Formation-the lower section of the Longmaxi Formation 1st member and the upper section of the Longmaxi Formation 1st member. Influenced by the collision between the Yangtze Plate and the Cathaysian Plate, during the deposition of the Wufeng Formation-the lower section of the Longmaxi Formation 1st member, volcanic activity was frequent and intense, leading to dramatic climate changes, improving the reducibility of the lower water [37, 38, 7680]. After the ash materials formed by volcanic eruption subsiding into surface water, minerals rich in Fe, P, and other nutrient elements for the propagation of plankton could be provided by hydrolysis, which was conducive to the improvement of biological productivity [37, 38, 8184].

Data from the Jiaoye-1 well proves the above text. As described in Figure 2, data from the Jiaoye-1 well shows that due to the strong collision between the Yangtze Plate and the Cathaysian Plate, the sealing of water is strong and the volcanic activity is active during the deposition of the Wufeng Formation-the lower section of the Longmaxi Formation 1st member, causing the U/Th ratio, the indicator of redox environment, to be the highest around 3 at the bottom of Longmaxi Formation. With the collision between two plates weakening, the sealing of water reduces and volcanic activity weakens and even stops during the deposition of the upper section of the Longmaxi Formation 1st member, causing the U/Th ratio to reduce to 0.5 from the bottom upward. The environment is gradually changed from anoxic environment to dysoxic environment, and then to oxidizing environment. Two indicators of biological productivity, excess molybdenum (MoXS) and silicon derived from bioorigin, show the same trend. During the sedimentation period of the Wufeng Formation-the lower section of the Longmaxi Formation 1st member, due to the strong water sealing and active volcanic activity, the two biological productivity indicators are both high; while during the sedimentation period of the upper section of the Longmaxi Formation 1st member, due to the weak water sealing and weakened and even stopped volcanic activity, these two biological productivity indicators are significantly reduced.

The strong sealing of water and active volcanic activity result in high biological productivity and reducing environment during the deposition of the Wufeng Formation-the lower section of the Longmaxi Formation 1st member, making the TOC content range in 2%~5% in this interval. With the water sealing weakening and volcanic activity gradually stopping, the TOC content of the upper section of the Longmaxi Formation 1st member descends to 0%~2% (Figure 2).

4.4. Models of Organic Matter Enrichment during the Deposition of the Upper Ordovician-Lower Silurian in the Upper Yangtze Area

As reflected in Figure 6(a), during the deposition of the Upper Ordovician Wufeng Formation-the lower section of the Lower Silurian Longmaxi Formation 1st member, due to the strong collision between the Yangtze Plate and the Cathaysian Plate, the deep shelf surrounded by ancient land is formed in the Upper Yangtze area, leading to poor connection with the sea and the strong sealing of water, which cause the stratification phenomenon of water. Meanwhile, due to the collision, the volcanic activity occurs in the ancient land surrounding the deep shelf, producing much ash. Volcanic materials are neutral-acidic magmatic rocks derived from deep in the crust. On the one hand, strong sealing led to high content of oxygen in upper water. On the other hand, lots of ash bring nutrient elements that are beneficial for plankton to grow and breed. Both factors cause the high biological productivity. Strong sealing and active volcanic activity can also make the lower water more reducible, which is conducive to the preservation of rich organic matter deposited from the upper layer, leading to the enrichment of sedimentary organic matter.

As depicted in Figure 6(b), during the deposition of the upper section of the Lower Silurian Longmaxi Formation 1st member, on the one hand, because of the weak collision between the Yangtze Plate and the Cathaysian Plate, the water becomes shallower, weakening the sealing, causing the hypoxic water in the lower layer to be mixed with the warm oxygen-rich water in the upper layer, and then, the reduction environment in the lower layer is damaged, and the oxygen content in the upper layer decreases and the water temperature becomes colder. On the other hand, with the collision weakening, volcanic activity becomes weak and even stops. So the source of ash declines; in other words, the nutrient elements beneficial for plankton to breed decrease. These two factors make the biological productivity of surface water decrease. Weakened sealing and gradually stopped volcanic activity reduce the reducing property of the lower water, causing the sedimentary organic matter from the upper layer to be damaged, which is not conducive to the accumulation of sedimentary organic matter.

5. Conclusions

The paper takes the marine shale of the Upper Ordovician Wufeng Formation-Lower Silurian Longmaxi Formation 1st member in the Upper Yangtze area, southern China, as the research objective interval. As the typical well, Jiaoye-1 well located in Sichuan Basin, Upper Yangtze area is the source of samples in the following experiments. The TOC content and mineral composition, major elements, and trace elements are analyzed by experiments on shale cores, shale debris, and bentonite cores, collecting the logging data to study the mechanism of sedimentary organic matter enrichment, and conclusions are as follows. (1)The siliceous minerals of the Upper Ordovician Wufeng Formation-Lower Silurian Longmaxi Formation 1st member shale in the Jiaoye-1 well, the typical well in the Upper Yangtze area, are derived from bioorigin except for the terrigenous clastic deposits, and the content of bioorigin siliceous minerals can be used as an indicator of biological productivity(2)During the deposition period of the Wufeng Formation and the lower section of the Longmaxi Formation 1st member, on the one hand, the oxygen content of surface water is high due to the strong sealing ability of water. On the other hand, the active volcanic activity brings ash which is beneficial for biology to breed. Both factors lead to high biological productivity of the water surface. At the same time, because the sealing ability of water is strong, reducibility of the lower water is strong. Besides, active volcanic activity causes climate change, increasing the sealing of water, which makes the organic matter subsiding well preserved, leading to the high TOC content of the Wufeng Formation-the lower section of the Longmaxi Formation 1st member today(3)During the deposition period of the upper section of the Longmaxi Formation 1st member in the Upper Yangtze area, on the one hand, the water becomes shallower, weakening the sealing of water, and the reduced water from the lower layer of water is brought to the surface, which caused the decrease of oxygen content in surface water. On the other hand, the volcanic activity weakens and even stops, making the source of ash with rich nutrient elements decrease. Both factors lead to low biological productivity of the water surface. At the same time, because the sealing of water weakens, the reducibility of the lower layer of water declines, and the inactive volcanic activity no longer affects the climate, which makes the little subsiding organic matter damaged, leading to the low TOC content of the upper section of the Longmaxi Formation 1st member today

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.

Additional Points

Highlights. (1) The excess siliceous mineral content of the Upper Ordovician Wufeng Formation-Lower Silurian Longmaxi Formation 1st member in the Jiaoye-1 well, a typical well in the Upper Yangtze area, is accurately calculated and the excess siliceous mineral as a biological productivity indicator is derived from bioorigin by elemental analysis. (2) The strong sealing and active volcanic activity during the deposition of the Wufeng Formation-the lower section of the Longmaxi Formation 1st member are conducive to the improvement of biological productivity and the accumulation of sedimentary organic matter. (3) The weak sealing and gradual cessation of volcanic activity during the deposition of the upper section of the Longmaxi Formation 1st member lead to the decrease of biological productivity and the oxidative destruction of sedimentary organic matter.

Conflicts of Interest

There are no conflicts of interest with respect to the results of this paper.

Authors’ Contributions

Xiaoxue Liu and Kun Zhang contributed equally to this work.

Acknowledgments

First and foremost, thanks for the support from the National Science and Technology Major Project (No. 2017ZX05035-002), the Science Foundation of the Ministry of Land and Resources of China (No. 12120114046701), the Natural Science Foundation of China (No. 41472112 and No. 41728004), and the open fund from Sinopec Key Laboratory of Shale Oil/Gas Exploration and Production Technology. Furthermore, I sincerely appreciate all coworkers who work hard with me.