Coevolutionary Dynamics of Organic-Inorganic Interactions, Hydrocarbon Generation, and Shale Gas Reservoir Preservation: A Case Study from the Upper Ordovician Wufeng and Lower Silurian Longmaxi Formations, Fuling Shale Gas Field, Eastern Sichuan Basin

Shale gas deposits are self-sourced, self-accumulating, and self-preserving in the Upper Ordovician Wufeng Formation and Lower Silurian Longmaxi Formation of the Fuling Shale Gas Field in the eastern Sichuan Basin. They were both seemingly mixed by secondary oil cracking and kerogen cracking gases during the high maturation window. The reservoir space primarily consists of mineral pores and organic matter (OM) pores, and the shale gas was mainly trapped by a high-pressure system. In this study, the Fuling O3w-S1l Shale Gas Field in the eastern Sichuan Basin was used as a case study to discuss the coevolutionary process and organic-inorganic interactions of hydrocarbon generation, accumulation, and preservation. The results indicate that the processes and mechanisms of organic-inorganic interactions and coevolution of hydrocarbon generation and reservoir preservation are quite different among the shale graptolite zones (GZ) with respect to hydrocarbon generation, types and characteristics of shale gas reservoirs, seal characteristics, and their spatiotemporal relations. In the WF2-LM4 GZ, the favorable OM, biogenic authigenic quartz and organic-inorganic interactions are highly coupled, leading to the high level of coevolution demonstrated within the field, as well as to the favorable conditions for shale gas accumulation. Conversely, the overlying LM5LM8 GZ seemingly exhibits early densification and late charge and has a reverse mode of reservoir development (i.e., low degree of coevolution). These two coevolutionary processes were conducive to the development of a high degree of spatiotemporal matching between the reservoir (i.e., WF2-LM4 GZ) and the seal (i.e., LM5-LM8 GZ). This is due to underlying differences in their coevolutionary histories. The synthetic work presented here on the coevolutionary processes and mechanisms of formation for organic-inorganic interactions and hydrocarbon generation and reservoir preservation reveals insights into the driving mechanisms of shale gas enrichment, providing a basis for effectively predicting favorable enrichment intervals for shale gas worldwide.


Introduction
During the past 10 years, most of the research on shale gas reservoirs has focused on source rock [1][2][3], reservoir types, porosity, and pore structure [4][5][6][7][8][9][10][11]. Most importantly, the pores within organic matter (OM) are widely recognized as a significant component of pore systems in gas-bearing shales [9,10,[12][13][14][15][16][17][18]. A greater potential for gas storage and flow exists in higher total organic carbon (TOC) content intervals, depending on the abundance and connectivity of the OM pores [8,16,19]. Therefore, shale gas reservoirs are mainly evaluated using TOC content, thermal maturity, rock mineralogy, rock brittleness, and porosity and permeability [1,2,20]. However, the combination of pore types, especially with respect to their components, characteristics, and organicinorganic interactions during the evolution of the shale reservoir, are less well known. This has restricted our understanding of the diagenetic and pore evolution histories, and the accuracy with which we may evaluate reservoir potential. The pore formation and the evolutionary histories and proportion of OM porosity are the key information needed to ascertain the effectiveness of exploration as well as the associated production capacity of shale gas reservoirs.
Porosity-related variances are directly related to differences in the amount and character of the OM and mineralogical composition [21][22][23]. The pore formation history is composed of the inorganic-organic interactions between OM (i.e., hydrocarbon generating organisms and/or bitumen) and inorganic matter (i.e., minerals such as quartz, feldspar, pyrite, and calcite). Organic-inorganic interactions in sandstones have been documented [24][25][26], and many important insights have been gained on the growth of minerals, the preservation of the pores that are affected by hydrocarbon charge [27,28], and the effects of minerals on the hydrocarbon generation of OM [29]. Meanwhile, the driving processes behind diagenetic pathways in shales [10,30,31], the nature of pores in shales, and the controls on porosity evolution remain poorly known [17]. Additionally, this does not take into consideration the effects of hydrocarbon generation and expulsion processes on mineral diagenesis, mineral pore evolution, and the relationship between organic (hydrocarbon generation) and inorganic (reservoir evolution) coevolutionary processes.
The petroleum retained along with kerogen and shale mineralogy may be a critical factor in shale gas generation [32], but their coevolutionary dynamics are little known. More attention is needed to develop our understanding of the coevolutionary processes that drive organic-inorganic interactions, hydrocarbon generation, and reservoir development of shale gas reservoir [20]. OM and inorganic minerals are not considered as unidirectional, simple, and determinate in their relations, but rather as a kind of interactional and coevolutionary relationship. The interaction between OM and inorganic minerals leads to their coevolution. The interaction is the foundation and motivation, and the coevolution is the process and consequence; both of them are achieved through physical and chemical processes. The key issue centers on how OM and inorganic minerals may interact to achieve coevolution, thus forming the reservoir types and characteristics that we observe as well, as the shale gas reservoir. Furthermore, the specific mechanisms and corresponding developmental states of pore systems (organic and inorganic) are deserving of more attention. Due to the particular system of self-contained, self-sourced, and selfaccumulating reservoir, the qualitative and semiquantitative evaluation of the organic-inorganic coevolutionary process may be conducted through the comprehensive analysis of hydrocarbon generation history, reservoir developmental processes, and preservation.
This study systematically analyzed the coevolutionary process and mechanisms of hydrocarbon generation and reservoir preservation based on organic geochemistry, mineralogy, pore characterization, seal characteristics, and the corresponding evolutionary history, combined with shale gas exploration and development in the largest Lower Paleozoic shale gas field in the world at Fuling, in the eastern Sichuan Basin. The results reveal the mechanism of shale gas enrichment in the Fuling Shale Gas Field and provide insights into the organic-inorganic interactive process and pore development, especially of the OM pore network within the studied shales. This may provide the basis for effective predicting favorable enrichment intervals of the Upper Ordovician Wufeng Formation and Lower Silurian Longmaxi Formation shale gas worldwide.

Geological Setting
As a superimposed basin developed on the upper Yangtze Craton, the Sichuan Basin is one of the most tectonically stable sedimentary basins in China [33]. Surrounded by the Hubei-Hunan-Guizhou Fold Belt on the east, the Longmenshan Fold Belt on the west, the Micang Uplift and Dabashan Fold Belt sitting on the north, and the Emei-Liangshan Fold Belt on the south, the basin can be divided into three tectonic regions and six sublevel tectonic units by the Huayingshan and Longquanshan anticlines [34] (Figure 1). Several tectonic movements occurred in the Sichuan Basin which in chronological order are the Caledonian (Cambrian-Silurian), Hercynian (Devonian-Permian), Indosinian (Triassic), Yanshanian (Jurassic-Late Cretaceous), and Himalayan (Paleogene-Quaternary) orogenies [34][35][36]. Sedimentary sequences of the basin are relatively complete, which include the Middle-Upper Proterozoic, Paleozoic, Mesozoic, and Cenozoic strata. Currently, China's largest gas yield comes from the Sichuan Basin which has the most prolific natural gas resources [33].
The black shale in the O 3 w-S 1 l Fuling Shale Gas Field (i.e., Wufeng and Longmaxi Formations) was deposited in a low-energy anoxic sedimentary environment, from the Late Ordovician to Early Silurian [38]. The maximum thicknesses of the formations of the O 3 w-S 1 l Fuling Shale Gas Field are between 30 and 100 m. In the Sichuan Basin and its surrounding areas where main China's shale gas yield comes from, the Wufeng and Longmaxi Formations have been the main targets for shale gas exploration, development, and production in recent years, leading to the establishment of three national shale-gas demonstration areas including the Fuling, Weiyuan-Changning, and Zhaotong.
To get a fine classification of shales with similar content, biostratigraphy can be an effective method when it is difficult to use lithology, mineral composition, geochemical parameters, and petrophysical properties. In a biostratigraphic study, graptolite is one of the leading globally recognized index fossil groups, which is a basis for fine shale division of the shales in the Wufeng Formation and the Longmaxi Formation [39]. Since specific graptolite types were only developed and preserved in the specific sedimentary environment during a specific period, the specific shales can be recognized and clarified 2 Geofluids [39][40][41]. Both the assemblage of hydrocarbon-forming organisms in different shale intervals and the shale quality were of great difference [42]. During the depositional period of the Wufeng Formation and the Longmaxi Formation, shale thickness and organic carbon content are changing; thus, it was more accurate to describe with a timeline the gradual development of shale characteristics based on the graptolite zone [39,43]. According to the GZ divisions in the Wufeng and Longmaxi Formations [39,40], the shales have been divided into WF1-WF4 in the Wufeng Formation and LM1-LM9 in the Longmaxi Formation.

Samples and Methods
More than 100 black shale samples from the Upper Ordovician Wufeng Formation-Lower Silurian Longmaxi Formation are used and sourced mainly from the shale-gas wells in the Sichuan basin and its periphery ( Figure 1) [44]. Measured shale porosity from the JY 1 Well averaged 4.5%, and the gas content in the basal 38 m of black shales was around 2.99 m 3 /t [44]. The main kerogen type in the WF2-LM8 GZ is considered to be dominated by Type II kerogen with some Type I [45], which favors the development of the OM pores [9]. And the development of OM pores between WF2-LM4 GZ and LM5-LM8 GZ is controlled by hydrocarbon-forming organisms [20,37]. Organic-rich sediments that lack terrigenous organic matter indicate a phytoplankton source [46]. Pelagic deposition of phytoplankton is interpreted as a primary mechanism of organic matter deposition, which is converted to hydrocarbons and matrix bitumen [36]. Based on the morphology, inherited texture, and energy spectrum, the OM can be linked to several hydrocarbon-generating organisms. This includes primary OM, such as various algae (multicellular algae, unicellular algae, acritarchs, etc.), graptolites, sponge spicules, and various types of secondary OM, such as bitumen, which occurs as discrete amorphous OM, and small L o n g q u a n s h a n A n t i c l i n e H u a y i n g s h a n A n t i c l i n e  Figure 1: Location map of the study area in the Sichuan Basin, China; the structural subdivisions include the eastern, southern, western, northern, southeastern, and central zones (after [37]). 3 Geofluids organic domains or organoclay aggregates [42]. Different hydrocarbon-generating organisms exist in different sedimentary environments, leading to the diversity of the type and content of OM, hydrocarbon-generating capacity, and OM pore development ability.
The multicellular algae include red, brown, and other planktonic algae. Pores within algal OM, particularly in smaller porous OM domains, are similar in shape, size, and distribution and may also be influenced by the original structure of the algae [16]. The diameter of single multicellular algae (probably red algae) ranges from several μm to tens of μm. Intense compaction may deform the shape of multicellular algae. The multicellular algae bear a mix of OM and silica since some of them have siliceous shells [42]. Other identified algae which present no siliceous shells fill in the mineral pores with their size varying with the size of mineral pores ( Figure 2). Without the protection of the siliceous shell, the shape of the algae can be easily altered. Under strong compaction, the algae shape deformation caused by internalstructure rearrangement of OM exhibits directionalalignment and/or flow-like structures in the narrow mineral pores. No directional alignment of the OM has been observed in circular or approximately equant mineral pores or in the center of narrow mineral pores, due to the insignificant difference in the degree of OM compaction.
The unicellular algae include blue, green, and other benthic algae. Unicellular algae mainly include rhabditiform blue algae (also known as cyanobacteria) and spherical algae and usually gather to form colonies ( Figure 3). The diameter or length of its colony is several μm to tens of μm. The length of single rhabditiform blue algae is a few hundred nm to more than 1 μm, and the width is a few hundred nm, while the spherical algae are circular or oval and have uneven surfaces and rough edges; some developed pores ranging tens of nm in size.
The red, brown, and other planktonic algae which can take advantage of green, yellow, orange, and other short wavelengths dominate the deep water zone, while the green algae, cyanobacteria, and other benthic algae are common in shallow waters [47,48]. The hydrocarbon-generating organisms of the WF2-LM4 GZ are primarily the red, brown, and other planktonic algae living in the aerobic surface layer and preserved in a strongly reducing environment in the bottom water [37]. A small amount of terrigenous OM is found in the hydrocarbon-generating organisms of the WF2-LM4 GZ [49]. In comparison, the hydrocarbon-generating organisms of the LM5-LM8 GZ were predominated by green algae, cyanobacteria, siliceous sponges, and other benthic algae, as well as a small amount of red, brown, and other planktonic algae [37]. TOC content is largely determined by different sources of hydrocarbon-forming organisms. The highest TOC content of 4.97% on average is found in siliceous mudstone of WF2-LM4 GZ; higher TOC content varying from 2.01% to 3.77% with an average of 2.84% appears in the silica-bearing mudstone of WF2-LM4 GZ. TOC contents of LM5-LM8 GZ are from 1.32% to 2.68% with an average of 1.70% in the silty mudstone and 0.60%-2.53% with an average of 1.59% in the clay-rich mudstone ( Figure 4). The WF2-LM4 GZs are where siliceous mudstone and silica-bearing mudstone are prevalent and have the high production of shale gas wells.

Hydrocarbon Generation History.
According to the burial history of Well JY 1, the black shale in the Wufeng and Longmaxi Formations underwent three stages of hydrocarbon generation ( Figure 5): (1) The initial stage was the rapid subsidence and burial in the Early Devonian. Buried at a depth less than 2000 m, thermal maturity (vitrinite reflectance, R o ) of the black shale was about 0.5% generating a small amount of heavy oil. The black shale was uplifted in the middle Devonian to the Late Permian when hydrocarbon generation was suspended. (2) The second stage was in the Early-Middle Triassic when the subsidence occurred. When the burial depth of the black shale went up to 3000 m, its thermal maturity climbed to about 1.0%, reaching the peak oil generation window. In the Middle-Late Triassic, again the black shale was uplifted and hydrocarbon generation was suspended, but because of a rapid, limited-amplitude uplift, oil expulsion efficiency was low, leaving a large amount of crude oil in the shale. (3) Subsidence from the early Jurassic to middle Cretaceous was the third stage. Black shale was buried at a depth greater than 7000 m with thermal maturity greater than 2.5%, turning into the overmature dry gas generation window. Residual kerogen and oil cracked generating plentiful shale gas as well as a large number of porous pyrobitumen which provided a huge number of OM pores for gas storage. After the Late Cretaceous, the black shale was uplifted again during which the generation of shale gas was terminated. The thermal maturity data range from 2.2% to 3.13% was present, with the mean value of 2.65% [20], suggesting the OM is already in the overmature stage.
According to Montgomery et al. [50], the Barnett Shale in the Fort Worth Basin of Texas, USA, was buried in a shallow depth and then uplifted for a long period during the early evolution stage, then was deeply buried in a short time and uplifted quickly in the middle evolution stage, and was quickly buried and uplifted in the late evolution stage. Similar characteristics have been found when comparing the burial and thermal evolution history of the shale from the Wufeng Formation and Longmaxi Formation in the Sichuan Basin ( Figure 5). During the early stages of their burial history from Early Silurian to Late Permian, both experienced long-term shallow burial and uplift, which was favorable for the kerogen preservation allowing a limited amount of hydrocarbon generation and expulsion. During the second period from Early-Middle Triassic to Late Triassic, short-term burial and uplift led to a large amount of crude oil generation and retention. Finally, the late period from the early Jurassic to recent witnessed rapid subsidence and uplift which was unfavorable for gas discharging, but favorable for shale gas generation, accumulation, and preservation.  (Figure 4), which is similar to the distributions observed for the Barnett and Eagle Ford shales [3,52], which favor the accumulation and production of shale gas.

Geofluids
Mineralogy may also be an important factor that controls both pore types and porosity preservation due to the differences in the ability to resist mechanical compaction. Siliceous mineral content is an important geological parameter in evaluating the quality and potential of shale reservoirs [52]. Thus, the content and origin of quartz must be resolved. Two different types of quartz, namely, detrital quartz and porefilling authigenic quartz, have been identified through petrographic observations, cathodoluminescence (CL), and SEM analysis [53]. The quartz content of shales in the JY 1 Well has a positive correlation with the corresponding TOC content. However, this statistic ignores the differences between siliceous minerals and TOC in different GZs. Further analysis reveals that the quartz content of the WF2-LM4 GZ is mainly greater than 40% which is dominated by biogenic silica [44]. It displays a strong positive correlation with TOC, guaranteeing a high matching of siliceous mineral and TOC. An increase in bioclastic components correlates with higher TOC content due to the higher volume of an organic-rich matrix. In contrast, the quartz content of the LM6-LM8 GZ varying from 20% to 40% is dominated by terrigenous silica, displaying no correlation with TOC content. As a continuous depositional transition between LM4 and LM6, the LM5 GZ has variable quartz content as a result of both terrigenous silica as well as biogenic silica. Siliceous shale is the main shale type of the WF2-LM4 gas-producing section, and the siliceous content is generally greater than 40%. Silica is often cryptocrystalline and amorphous, exhibiting spherical structures (Figures 6(a) and 6(b)), and there are usually sponge spicules and skeletal material, including radiolarian fragments (Figures 6(c) and 6(d)).
The detrital versus an authigenic form of quartz can be discriminated firstly by crystal shape and spatial distribution and then by color and intensity in CL. CL can also be used to determine quartz provenance [52,[54][55][56]. Low-luminescent or nonluminescent quartz observed in scanning monocathodoluminescence electron microscope is commonly attributed to a low-temperature, authigenic/diagenetic origin [54,55]. In this study, a resolution SEM with monocathodoluminescence and energy-dispersive Xray spectroscopy (EDS) detectors were used to identify the provenance of quartz within the Wufeng and Longmaxi Formation shales. Two different types of quartz were identified [57]: (1) Biogenic microcrystalline quartz  (1-3 μm) and irregular microcrystalline aggregates among the detrital grains that were weakly luminescent or nonluminescent ( Figure 7). Biogenic microcrystalline quartz produces a broad luminescence peak between 620 nm and 660 nm (Figure 7(c)), which is similar to the luminescence generated by low-temperature authigenic quartz due to oxygen vacancies [55]. (2) Silt-size quartz grains (5-20 μm) are characterized by a uniform luminescence intensity (Figures 7 and 8). The subrounded and subangular silt-sized quartz grains produce two luminescence peaks with wavelengths of between 390 and 420 nm and 620-640 nm (Figure 8(c)), suggesting a detrital origin [56].

Geofluids
The bulk of opaline silica produced by dissolved biogenic silica of radiolarians and sponge spicules in ocean surface waters does not become incorporated into deep-sea sediments because it dissolves while settling down [54]. Sediment accumulation rate is a significant factor in mudrock cementation [58], especially for the formation of biogenic quartz. In the situation in which there was a slow rate of sediment accumulation in the WF2-LM4 GZ, the amount of dissolved biogenic silica from radiolarian tests and sponge spicules is substantial. In contrast, rapidly deposited shales of the LM5-LM8 GZ tend to have only trace amounts of biogenic silica. In the processes of burial and exposure to high temperatures and pressures, the biogenic opal-A is first transformed into opal-CT and eventually into biogenic quartz [52,54], which is mainly composed of aphanitic, microcrystalline, and macrocrystalline quartz. This may be the reason why there is more biogenic silica in the WF2-LM4 GZ than in the LM5-LM8 GZ. Up to 60% of total quartz has a biogenic origin in siliceous shales of the WF2-LM4 GZ [57], resembling the Eagle Ford Formation in Texas, in which the authi-genic quartz comprises 85% percent and is equal to about 12.6% by volume [52]. There is more biogenic silica in the WF2-LM4 GZ than the silt-size detrital quartz in siliceous shales of the LM5-LM8 GZ. Biogenic quartz cement as matrix-dispersed microcrystalline quartz plays an important role in the diagenesis of the Wufeng and Longmaxi Formations.
Microcrystalline quartz in siliceous shales is mainly formed during diagenesis and is similar to the microcrystalline quartz in the Barnett Shale [3] and Eagle Ford Shale [52]. Thin section analysis and other methods reveal that there are different hydrocarbon-generating fossil organisms in the WF2-LM4 GZ. In addition to the graptolites, algae, and acritarchs, there are many siliceous microfossils, such as radiolarians, sponge spicules, and other skeletal fragments ( Figure 6), which indicate the organic origin of silica in the organic-rich shale of the WF2-LM4 GZ. The content of detrital quartz and feldspar in siliceous shales is relatively low, which is the main difference between the organic-rich and the silty shale (the main shale type of LM5-LM8 GZ).

Geofluids
Radiolarians and sponge spicules are the dominant forms of silica in the intrabasinal detrital assemblage. Dissolved silica released by dissolution of the radiolarians and sponge spicules is the dominant source of silica for the precipitation and reprecipitation of authigenic quartz in a variety of forms, including matrix-dispersed microcrystalline quartz cement, and fillings within primary intragranular pores [52]. The matrix-dispersed biogenic quartz is the dominant quartz in WF2-LM4 GZ (Figures 7 and 8). The matrix-dispersed authigenic microcrystalline quartz in the Wufeng and Longmaxi Formations occurs as the grain-binding cement that plays a role in the generation of brittle mechanical rock properties, similar to the Barnett Shale [30]. Additionally, the specific form of matrix-dispersed authigenic quartz has significant implications for the rock mechanical properties and residual oil. Moreover, as a result of precipitating authigenic quartz from the dissolution of biogenic opal [54], preserved quartz content (i.e., excess silica) may also reflect the surface water productivity during deposition [59], which indicates a higher surface productivity with a higher supply of OM during shale deposition.

Pore
Types. The pore types in the Wufeng and Longmaxi Formations include OM pores, inorganic mineral pores (particle-related, pyrite-related, and clay-related pores), and microfractures [10,11,33]; the most important of these are the OM pores. Greater gas storage and flow capacity exist in shale intervals with higher TOC and depend upon the abundance and connectivity of OM pores [7,10,14,60].
Observations demonstrate a complex evolution of the OM pore system in shales [60]. The difference in OM types (hydrocarbon-generating organisms and amorphous OM (bitumen)) has an impact on the type, size, and distribution of the OM pore system [42]. The OM type, controlled by sequence-stratigraphic position [13], is a key factor in determining the pore type, size, and distribution. The OM pores include pores hosted in algae, biodetritus, and bitumen. The algae are a structured OM type, which generally has consistent or predictable porosity styles in the shale section due to the different algal types. Proceeding from the bottom of the Wufeng Formation to the middle and upper parts of the first member of the Longmaxi Formation, the OM pore types and pore size changed regularly.

Geofluids
(1) The shale of the WF2-LM4 GZ was found to be the interval with the most abundant OM present, in which the most abundant type of OM is multicellular algae, as well as small amounts of unicellular algae and bitumen [42]. Under compaction, the OM formed thin layers and/or three-dimensional interconnected OM pore networks ( Figure 2). The OM pores were left behind by various hydrocarbon-generating organisms (mainly multicellular algae) and oil during gas generation and expulsion. The high total amounts of multicellular algae ensure a large number of OM pores, favoring the formation of three-dimensional interconnected pore networks, which is conducive to the enrichment and production of shale gas. Additionally, the three-dimensional reconstruction of kerogen and porosity distributions reveal that OM pores have the potential to form connected pore networks, where 26%-67% of the pore volume may be connected [16]. Moreover, the mineral pores are predominately slitlike in the cross-section, whereas OM pores are typically round [11]. OM pore networks are a combination of pores developed in depositional OM and migrated OM [14]. The depositional OM (mainly multicellular algae) pores are the main reservoir type in the shale of the WF2-LM4 GZ [37].
(2) The OM types of the LM5-LM8 GZ are mainly composed of unicellular algae and amorphous OM (bitumen), and the TOC is lower than the shale of the WF2-LM4 GZ [19]. The OM pores are mainly developed in unicellular algae and/or little in solid bitumen which often occur as poorly developed isolated pores and in unicellular algae which tended to be randomly developed by the arrangement gaps ( Figure 3). Therefore, this kind of OM pore development is limited. The shape is irregular, and the size is tens to more than one hundred nm. The OM pore system changes are less pronounced and cannot form three-dimensional interconnected OM pore networks due to the low TOC content, which is not conducive to the enrichment and production of shale gas.

Direct Seal
A series of tight strata directly overlying the shale gas reservoir is served as the cap rock, such as shale, silty shale, argillaceous siltstone, and limestone. A prerequisite for shale gas source and seal coupling systems is the heterogeneity of shale strata, while key factors determining the capability of the direct seal are lithology and petrophysical properties, as well as the differences of petrophysical parameters between the   10 Geofluids source and seal [44]. The LM5-LM8 GZ is the direct seal of the WF2-LM4 GZ (Figure 4). For shale mineral composition, the quartz content decreases while the clay mineral content increases from the source to seal, the Wufeng Formation to the first member of Longmaxi Formation. With differences in mineralogy, LM5-LM8 GZ is less brittle but more ductile than the WF2-LM4 GZ when compacted by equivalent geological stress. Observably high gas content in the WF2-LM4 GZ and a significantly lower one in the LM5-LM6 GZ sequence in the Well JY 1 are shown in Figure 4, indicating that the upper part of the Longmaxi Formation (the LM5-LM6 GZ) served as a direct seal with its much lower porosity and higher clay content preventing shale gas in the WF2-LM4 GZ from escaping.
Furthermore, the capability of direct seals can be quantitatively estimated from breakthrough pressure and median pressure. The breakthrough pressure and median pressure of the WF2-LM4 GZ is lower than those of the LM5-LM8 GZ according to the microscopic pore structure analysis. There is a positive relationship between the breakthrough pressure and the breakthrough radius as well as the median pressure and the median radius indicating that the LM5-LM8 GZ plays a role as the direct seal for the WF2-LM4 GZ with its strong sealing capacity [44]. The breakthrough pressure and median pressure, together with the breakthrough radius and median pressure between WF2-LM4 GZ and LM5-LM8 GZ of the Guanyinqiao section in Qijiang County, Chongqing City, have similar characteristics to Well JY 2, that is, the LM5-LM8 GZ is the direct seal for the WF2-LM4 GZ (Figure 9).

Coevolutionary Processes and Mechanisms of Shale Hydrocarbon Generation and Reservoir Preservation
7.1. Source and Reservoir Evolution. Different GZs represent different deposition environments, and this is also true for corresponding hydrocarbon-generating organism assemblages and mineral compositions [19,39]. Combined, these lead to the differences observed in the interrelationships and coevolution of hydrocarbon generation and reservoir formation processes.
(1) The shale of the WF2-LM4 GZ is siliciclasticdominated featuring a slow deposition rate, favorable organism types (i.e., multicellular algae), high TOC, and high hydrocarbon generation capacity. All of these are then favorable conditions for the further generation, accumulation, and preservation of shale gas. The timing sequence of hydrocarbon generation and reservoir evolution in the WF2-LM4 GZ exhibits signs of a siliceous grid that was formed early and oil Median radius (nm) LM5 and the above graptolitic zone WF2-LM4 graptolitic zone Figure 9: The breakthrough pressures, breakthrough radius, median pressures, and median radius (mercury injection capillary pressure data) of the Wufeng and Longmaxi Formations, Guanyinqiao section, Qijiang County, Chongqing Municipality. Because the shale of the outcrop sample is weathered to different degrees, the breakthrough radius and the median radius are large, while the breakthrough pressure and the breakthrough radius are small. However, the overall performance is that the breakthrough pressure and median pressure of the WF2-LM4 GZ is lower than those of the LM5-LM8 GZ, showing the stronger sealing capacity of the LM5-LM8 GZ which is a direct seal for the WF2-LM4 GZ.   Table 1): (a) Quartz cement in the WF2-LM4 GZ may originate from various silica sources, such as biogenic authigenic silica, mineral transformations, and redistribution (pressure dissolution/clay-induced dissolution) [57]. Assessing the relative contributions of these potential silica sources would require the determination of the relative timing of the authigenic quartz [52]. Understanding the formation of biogenic quartz cement is of great importance for shale gas reservoir evaluation, especially when the formation time of three-dimensional connected frameworks of biogenic silicon is earlier than that of the maximum oil generation, which can provide storage spaces for holding a large amount of oil.
Most of the quartz in the WF2-LM4 GZ is biogenic, which may originate from the dissolution and recrystallization of biogenic silica [44]. The documented temperatures and porosities in effect during the alteration of former biogenic opal-A through the transformation to opal-CT and eventually to microcrystalline quartz [52,54,61,62] are consistent with the observed burial histories. The transformation from opal-A to crystalline microcrystalline quartz has well-documented impacts on the generation of brittle behavior in mudstones [61]. Unconfined compressive strength (UCS) exhibits an abrupt increase from several MPa above 1000 m to several tens of MPa at a depth of~1000 m, where the silica is opal-CT, and the UCS remains at high values at depths of 1000-1200 m; although displaying a gradual decline over this depth range, the UCS is still 10-15 MPa below the depth of 1200 m [61]. The results show that when the buried depth of shale is at least 1000-1200 m, the silica is opal-CT and the shale has strong brittleness. According to the evolution history of opal-A→opal-CT→quartz, oxygen isotope analysis shows that the temperature of opal-CT that evolved into quartz-layered flint is 17-21°C [63]. The thermal maturity is lower than 0.5% R o when the temperature is 17-21°C and the burial depth is less than 1000 m, and the formation time for authigenic quartz is Late Early Silurian (early diagenetic stage). Thus, the transformation process seemingly occurred at relatively low temperatures before the oil window. Rapid compaction occurs in response to the gradual increase in burial depth and the infilling of pore spaces in preexisting detrital quartz frameworks with biogenic microcrystalline silica [62].
Intergranular pore space is controlled by the preservation of primary pores [8], which are prone to be preserved under the supporting structure of matrix-dispersed authigenic quartz that prevents compaction [57]. This dissolution and reprecipitation process may transform the siliceous particles in place or reprecipitate dissolved silica in the forms of matrix-dispersed microcrystalline quartz and irregular microcrystalline aggregates in pores that were formed by early detrital quartz [57]. The shape of kerogen particles (hydrocarbon-generating organisms) and the amount of retained oil filled in primary framework pores while relying on the presence of the compaction-resistant rigid grains that are determined by detrital and matrix-dispersed biogenic authigenic quartz for pressure maintenance to preserve porosity. These matrix-dispersed microcrystalline quartz aggregates have an early formation time before the infilling of hydrocarbon compaction occurs. Pore space is abundant enough in microcrystalline quartz aggregates for oil and bitumen filling during the oil window. Source and reservoir evolution Strong ability to form a rigid grid and strong anticompaction ability; high degree of coevolution Weak ability to form a rigid grid and strong anticompaction ability; low degree of coevolution Source and seal evolution Usually as the direct seal for shale gas reservoir Promising shale gas interval 13 Geofluids Matrix-dispersed biogenic authigenic quartz occurs as the cement in the primary, irregular, intergranular, and intragranular pores; restrains the compaction; and preserves the internal pore structure together with the detrital quartz as a rigid and three-dimensional framework. This then strengthens the framework and effectively increases the shale brittleness at the same time. The rock matrix formed by high quartz content has a strong ability to resist compaction that provides abundant spaces for maintaining OM, especially for the retention of large amounts of oil. Furthermore, during the hydrocarbon generation window, the organic acids generated can dissolve the unstable minerals [64,65], which also increases the storage spaces. The filling of early oil may have some inhibiting effects on shale diagenesis and slow the compaction process.
(b) The TOC is one of the main controlling factors for hydrocarbon generation. The TOC content is generally high, commonly between 2% and 5%, and the OM mainly exists in flocculent states in the siliceous particles or between the quartz particles. Organic petrology coupled with backscattered and field emission SEM indicates that the OM type (mainly multicellular algae) is the most significant contributor to the TOC content and OM porosity, ensuring the matching of gas generation and OM pore development. The thermal simulation experiment results show that the generated oil yield from modern planktonic algae is almost double or triple that of benthic algae [48]. The WF2-LM4 GZ dominated by the multicellular algae had a much higher oil generation potential than the LM5-LM8 GZ where unicellular algae were common. The organic-rich WF2-LM4 GZ is favorable to oil generation and adsorption, as well as primary oil expulsion to the intergranular pores of biogenic quartz matrices, bringing in high oil retention in TOC-rich and biogenic quartz-rich graptolitic shales. It is generally assumed that multicellular algae have complex cell structures and a stronger ability to generate multiple kinds of OM pores than unicellular algae during thermal maturation [9]. The OM pores in the multicellular algae are well developed within the algae and inherited the original biological structure [42], the sizes of which are usually hundreds of nm, or even μmscale. Additionally, the TOC content also determines the number of OM pores.
(c) The retention of hydrocarbons in shales of the WF2-LM4 GZ is mainly controlled by the sorption capacity of OM and the spaces in threedimensional connected frameworks of biogenic silica [19,44], in which the oil was retained in absorbed and free states, respectively. A similar retention phenomenon related to the presence of porous fossils exists in the Marathon 1 Mesquite Well [66], where siliceous sponge spicules enhance the storage capacity for bitumen [67]. For gas generation, the secondary cracking of the residual oil is the main source [1,32,33,68,69]. Shale gas content is largely determined by the amount of residual oil cracking during high thermal maturity in organic-rich shales [69,70]. The WF2-LM4 GZ has high content of TOC and silica that is associated with multicellular algae [37]. During the oil window, the generated oil was filled in kerogen and/or expelled from the kerogen and migrated into the adjacent interparticle and intraparticle mineral pores. With increasing maturity, the normalized concentration of retained oil (S1/TOC·100) TOC increased in the oil window [67]. During the stage of oil cracking, the generation of gas, the formation of OM pores, and overpressure coevolved to offset the effects of compaction. In particular, the residual oil and overpressure have important effects on reservoir formation, which are conducive to the development and preservation of OM pores.
However, with increasing maturity, the cracking and releasing of bulk labile hydrocarbon moieties from kerogen structures might likely leave behind a more rigid and aromatic kerogen residue [67,71], while the residual oil similarly leaves behind a more rigid and aromatic bitumen. Therefore, the formation of OM pores is most likely a result of the release of massive amounts of hydrocarbons [19,67]. The abundance of OM pores and degree of connectivity in the three-dimensional OM pore network synergistically increased with an increasing authigenic quartz content in the WF2-LM4 GZ. Residual oil principally fills into the aggregates of authigenic microcrystalline pore networks with the deposition OM and would be likely to generate a better connected three-dimensional OM pore network through secondary cracking, which provides favorable storage space for shale gas in shales of the WF2-LM4 GZ.
To summarize the results above, a logical conclusion is that the coevolution of high contents of favorable OM and quartz (multicellular algae and biogenic quartz) provided high hydrocarbon generation potential and connected threedimensional networks that may have resulted in the retention of large amounts of oil and kerogen in the rigid interparticle and intraparticle mineral pores. During the stage of oil cracking, the generation of gas and overpressure coevolved to produce an environment conducive to the development and preservation of OM pores and 14 Geofluids to offset compaction. In particular, the WF2-LM4 GZ offers a promising reservoir with favorable types and characteristics that are the result of multifactor coupling and coevolution. Reservoir potential should thus not be viewed as being determined by a single factor, but rather as the synergy among many factors, such as hydrocarbon generation potential, porosity, TOC, and quartz contents.
(2) The shale of the LM5-LM8 GZ are clay-dominated and characterized by a fast deposition rate, unfavorable organic organism type (mainly unicellular algae), low TOC, and low hydrocarbon generation capacity, all of which are not favorable conditions for the generation, accumulation, and preservation of shale gas. The timing sequence of hydrocarbon generation and reservoir evolution of the LM5-LM8 GZ experienced early densification of the reservoir and oil that was charged late; the main features are as follows: (a) The siliceous minerals of the LM5-LM8 GZ are mainly terrigenous and detrital in origin with trace amounts of biogenic authigenic quartz [19]. The clay mineral content is higher than that of WF2-LM4 GZ. The most unfavorable factor is that the compaction resistance ability of these terrestrial detrital silica and clay minerals is weak and cannot form well-connected three-dimensional frameworks. The radiolarians and sponge spicules slowly gathered and dissolved due to the high deposition rate. Furthermore, the precipitation amount of biogenic authigenic quartz is also low. Without this support, the framework of terrigenous siliciclastic silica with high clay content collapses easily during strong compaction. This densification of the shale reservoir and the reduced porosity and permeability limit the storage space for the retention of oil and migrated OM patches. The transforming period of quartz produced during clay transformation was consistent with the oil generation window that further reduced the pores for the retention of oil and migrated OM patches [57,72].
(b) The hydrocarbon-generating organisms of the LM5-LM8 GZ were mainly unicellular algae, with lower TOC content. The pores in unicellular algae are developed by the arrangement gaps of unicellular algae, in which the size is tens to greater than one hundred nm [42]. These two factors led to a poor hydrocarbon generation ability, fewer OM pores, isolated pore distribution, and an inability to form an interconnected three-dimensional OM pore system, all of which is unfavorable to the enrichment and production of shale gas.
(c) The amount of natural gas generation in the LM5-LM8 GZ largely depends upon the amount of retained oil generated by authigenic OM and the oil that may be migrated from the WF2-LM4 GZ. Due to the small amount of selfgenerated oil and pores for oil retention, the storage space for natural gas was limited during the oil cracking stage. Additionally, both the kerogen and the residual oil have low capacities to generate gas that make it difficult to form OM networks when there is little to no biogenic authigenic quartz present. In this case, only trace amounts of an OM network exist in the LM5-LM8 GZ.

Source and Seal Evolution.
A good match between the shale gas generation from an organic-rich shale source and the development of the direct seal is critical for the assessment of shale gas enrichment in the Sichuan Basin, South China. The pore evolution simulation performed for TOC and mineral compositions of the WF2-LM4 GZ and LM5-LM8 zones in Well JY 1 showed that the porosity of OM pores in the WF2-LM4 GZ increased at a greater burial depth with the thermal maturity of 0.6%, while the porosity of the mineral pores gradually decreased with evolution. The Middle Cretaceous was the peak period for dry gas generation and OM pore formation in the WF2-LM4 GZ. Yet in the LM5-LM8 GZ, the mineral pores are dominant and with little OM pores. The porosity of the LM5-LM8 GZ decreased due to the lack of supplement to increase OM pores during the overmature dry gas generation window, ranging from 3.19% to 5.09% ( Figure 11). The WF2-LM4 GZ is the source and reservoir for shale gas with its high TOC high residual oil content and OM-hosted porosity, while the overlying LM5-LM8 GZ is the direct cap rock. The shape of porosity-decline curves is defined by initial porosity, grain size, and mineral composition [58]. In short, they are governed by the microfabric [73], and the microfabric, in turn, controls the reservoir types and characteristics. To obtain the evolution of porosity of the two GZ, the burial history reconstruction and the thermal history simulation were performed, revealing the differences between the source (WF2-LM4 GZ) and the direct seal (LM5-LM8 GZ). The porosity of mineral pores in the WF2-LM4 GZ exhibited a decreasing trend with burial depth, dropping to a minimum value of 1.0% to 1.5% in the Middle Cretaceous. Then, the porosity of mineral pores presented a slight increase accompanying the uplift, reaching~2.0% to the recent. On the other hand, the porosity of OM pores in the WF2-LM4 GZ showed an obvious increase during the Early Triassic when the thermal maturity reached 0.6% R o , reaching 2.5% in the Early-Middle Cretaceous. Since then, under the decompression during the uplift, it has slightly increased to~3.0% till the present. The average total porosity is about 6% in the WF2-LM4 GZ, much greater than that in the LM5-LM8 GZ ( Figure 11). The porosity of the OM pores is about 1% in the LM5-LM8 GZ. Thus, the LM5-LM8 GZ can be regarded as a direct seal for the shale gas reservoir in the WF2-LM4 GZ. In conclusion, the source rock and direct seal form a strong coevolutionary matching, favorable for shale gas accumulation. (1) In the Early Devonian, an initial stage of rapid subsidence and burial occurred during which only a small amount of heavy oil was generated. With the small amounts of detrital quartz, the matrixdispersed biogenic authigenic quartz in the WF2-LM4 GZ occurs as the cement in primary anomalous intergranular and intragranular pores that restrained the compaction to preserve the pore structure as rigid and connected threedimensional frameworks strengthened the rock framework and increased the shale brittleness effectively at the same time. Without the support of biogenic silica, the framework of terrigenous siliciclastic silicon with high clay content collapses easily during compaction, which results in smaller storage spaces for retaining oil and migrated OM patches.
(2) The second stage of subsidence occurred in the Early-Middle Triassic. A large amount of oil might remain in the shales of the WF2-LM4 GZ because of a rapid, limited-amplitude uplift and a low oil expulsion efficiency. Existing biogenic authigenic quartzsupported three-dimensional framework porosity is filled with migrating hydrocarbons and subsequently results in the development of OM micropores in those places. During this period, the direct seal rock was formed when the shale reached high maturity at~1.2% R o .
(3) The third stage of subsidence extended from the Early Jurassic to Middle Cretaceous. The kerogen and residual oil can generate a great amount of gas and bitumen through cracking, also resulting in large amounts of OM pores for gas storage when the shale reached maximum maturity. Meanwhile, this time was the most productive gas generation period for the WF2-LM4 GZ, that is, peak gas generation and formation of massive OM pores were synchronous. However, the porosity of the LM5-LM8 GZ was at its minimum in this period. During the periods of largest gas generation, OM pore formation, seal formation of direct cap rocks overlapped and formed a strong matching for shale gas accumulation and preservation in the Wufeng and Longmaxi Formations (Figures 7 and 11). Meanwhile, the regional seal of the Lower Triassic gypsum layer was formed during the Late Jurassic-Early Cretaceous [74], before the periods of greatest gas generation and OM pore formation.
(4) The acceleration of the uplift stage from the middle stage of the Late Cretaceous led to the cease of the hydrocarbon generation process and diagenesis process resulting from the thermal maturation evolution. Besides, pores in the minerals resulting from the decompression process turn larger, the organic pores increase to different degrees, and the preexisting fissures reopen [11]. If the tectonic uplift is too strong, the large and deep faults would develop in the shale and exodiagenesis could occur with freshwater pouring into fractures. And all kinds of dissolution processes of inorganic minerals would do great damage to the shale gas reservoir. What is more, the organic   Figure 11: Porosity evolution history of the Wufeng and Longmaxi Formation shale gas reservoir in the Jiaoshiba area (modified from [20]). 16 Geofluids diagenesis and organic-inorganic interaction process turn relatively weak.
It should be noted that the research of coevolutionary dynamics of organic-inorganic interactions, hydrocarbon generation, and shale gas reservoir preservation of the Upper Ordovician Wufeng and Lower Silurian Longmaxi Formations between the WF2-LM4 GZ and the LM5-LM8 GZ is a theoretical study of the static characteristics, coevolution process, evolution mode, and evolution characteristics of two graptolite zones. However, these are controlled by the black shale depositional environment, the tectonic evolution experienced by the shale strata after deposition, and the current structural location, type, and characteristic. The shale gas reservoir characteristics and shale gas enrichment characteristics of different graptolite zones may not be completely consistent with theoretical models, and specific analysis needs to be combined with the sedimentary environment and structural background in the study area. Taking the Upper Ordovician Wufeng and Lower Silurian Longmaxi Formations in the Sichuan Basin and its surrounding areas as an example, it is necessary to determine the thickness and sequence evolution characteristics of the WF2-LM4 GZ based on the shale depositional environment, structural background, and exploration and development costs. To further predict the shale gas enrichment area, mapping the WF2-LM4 GZ will be a major task in future exploration and a lot of work needs to be carried out.

Analysis of Wufeng and Longmaxi Shale
Gas Potential 8.1. The Coevolution Processes between the WF2-LM4 GZ and the LM5-LM8 GZ. Since the hydrocarbon-generating organisms of the WF2-LM4 GZ are mainly multicellular algae that have high TOC content and a high hydrocarbon-generating capacity, the organic acids generated are greater than those of the LM5-LM8 GZ. The generated organic acids preferentially dissolve the unstable minerals in the shale and typically chemically react with feldspar particles to precipitate kaolinite and microcrystalline quartz [75]. The dissolution of unstable minerals further reinforces and preserves the connected three-dimensional framework spaces of biogenic silica formed before the organic acid formation, which is favorable for oil retention in the shale of the WF2-LM4 GZ. Strongly coupled with the large amounts of oil generation that is determined by the multicellular algae, the kerogen and the biogenic silica-connected three-dimensional frameworks can retain much oil, which greatly determines the gas content in shale by generating cracking gas and porous pyrobitumen at high thermal maturities. Previous reports on adsorbed oil in kerogen proposed that 170 mg/g TOC can contribute tõ 22% of free oil, estimating that shale has great potential for oil generation [69]. This constitutes a large portion of the retained oil in shales. All of these favorable factors (i.e., favorable OM type and TOC content, the formation of connected three-dimensional networks of biogenic authigenic quartz before oil charge, and overlapped generation periods of natural gas, OM pores, and overpressure) suggest a high level of source-reservoir coevolution to provide ideal conditions for shale gas accumulation.
Because of the lack of biogenic silica needed to form connected three-dimensional frameworks in the LM5-LM8 GZ, the dissolution of unstable minerals by organic acids and the precipitation of microcrystalline quartz transformed by clay minerals strengthened the compaction and further reduced the porosity. The storage space for oil is smaller than in the WF2-LM4 GZ. By the coupling of a low oil generation potential (driven by unicellular algae) and small storage spaces for oil retention, a small amount of gas may be hosted in various pores and dominantly in the mineral pores, which reveals a low level of source-reservoir coevolution.
The two coevolution processes are conducive to a high spatiotemporal matching of the reservoir (WF2-LM4 GZ) and seal (LM5-LM8 GZ) due to their different coevolving processes of gas generation, OM pore formation, and direct seal formation.
Shale gas displays the feature of being a self-sourced and self-accumulating reservoir and mixes by the cracking gases of oil and kerogen at high maturity [1,2,76]. The WF2-LM4 GZ is the source and reservoir with favorable accumulation conditions for shale gas development, characterizing a slow deposition rate, favorable organism type, high TOC and high hydrocarbon generation capacity, high quartz content, and high degree of coevolution between hydrocarbon generation and the reservoir [19,20]. At present, the WF2-LM4 GZ is the main prospect section of exploration, and three national shale gas exploration areas in the Fuling, Weiyuan, Changning, and Zhaotong, have been constructed where the thickness of the WF2-LM4 GZ is over 20-25 m [20]. The WF2-LM4 GZ is the most profitable GZ for shale gas accumulation, and it displays excellent exploration potentials. Even outside the Sichuan Basin where the tectonics are relatively strong, the WF2-LM4 GZ can also form good normal pressure shale gas reservoirs [77,78].
8.2. The Shale Gas Exploration Potential of the LM5-LM8 GZ. Shale gas exploration also exhibits generally unsatisfactory production from the shale in the LM5-LM8 GZ except where structural highs have enhanced hydrocarbon accumulation [79]. The LM5-LM8 GZ directly overlaps the high-quality source rocks of the WF2-LM4 GZ, and the local poredeveloped areas typically have a good ability to store oil and maybe reservoirs for migrated oil. During the oil cracking stage, the oil will crack to natural gas to form porous bitumen and store the locally generated shale gas. There may also be some natural gas migrated from the WF2-LM4 GZ. The LM5-LM8 GZ may also have shale gas accumulation which has been confirmed by exploration [80]. The JY 11-2-1HF Well and JY 29-S1 Well are located at the apex and the slope of the Jiaoshi box-shaped anticline, respectively, and the depths of the targets are 2200 m and 2400 m, respectively. In June of 2015, the JY 11-2-1HF Well operated by Sinopec reached a test production of 287,000 m 3 /d after hydraulic fracturing (the horizontal interval only has a length of 780 m), which is a strategic milestone in the LM5-LM8 GZ. However, the JY 29-S1 Well reached a test production of 248,000 m 3 /d after the hydraulic fracturing (the horizontal interval has a length of 1502 m). The horizontal interval of the JY 29-S1 Well is almost twice that of the JY 11-2-1HF Well, but the test production is nearly identical. The permeter production of the JY 11-2-1HF Well is thus twice that of the JY 29-S1 Well. It is believed that there may be a primary migration in the LM5-LM8 GZ, resulting in a high degree of shale gas enrichment in the upper part of the Jiaoshiba anticline [79,80]. It is suggested that the development of LM5-LM8 GZ should start from the high part of the structure and gradually expand to the low part.

Conclusions
In our study, the Fuling O 3 w-S 1 l Shale Gas Field in the eastern Sichuan Basin was used as a case study to discuss the coevolutionary processes and organic-inorganic interactions of hydrocarbon generation and reservoir preservation. This comprehensive analysis revealed that the dynamics of organic-inorganic interactions and coevolution in hydrocarbon generation and reservoir preservation are quite different among shale GZs in terms of hydrocarbon generation, types, and characteristics of shale gas reservoirs; seal characteristics; and their spatiotemporal relationships.
The WF2-LM4 GZ contains planktonic algae featuring the majority of multicellular algae, high TOC, strong hydrocarbon generation capacity, and a high-quality shale gas reservoir with connected OM pore networks. Taken together, all of these result in favorable conditions for the development of a high degree of coevolution in hydrocarbon generation and accumulation. The high degree of coupling among favorable OM, biogenic authigenic quartz, and organic-inorganic interactions likely led to coevolution and the favorable conditions for shale gas accumulation.
However, the overlying LM5-LM8 GZ exhibits signs of early densification and late charge and has a reverse mode of reservoir development. This is likely due to the lack of biogenic microcrystalline opal-CT/quartz precipitation to preserve and reinforce the preexisting framework, together with the unfavorable matching between oil generation and pore development, which cannot provide enough storage spaces for oil retention to generate widespread cracking gas and OM pores during the high maturation window.
These two coevolutionary processes are conducive to a high degree of spatiotemporal matching between the reservoir (i.e., WF2-LM4 GZ) and the seal (i.e., LM5-LM8 GZ). This may be due to their different coevolutionary histories of gas generation, OM pore formation, and direct seal formation, which laid the foundation for the coevolutionary dynamics of hydrocarbon generation and reservoir preservation.
We believe that our comprehensive research and the mechanisms formulated for the coevolution of organicinorganic interactions and hydrocarbon generation and reservoir preservation have determined that the WF2-LM4 GZ is the most promising interval for shale gas exploration and development in the Sichuan Basin. Furthermore, the insights provided here may be used as a framework for the geological theory, exploration, and development of shale gas worldwide.

Data Availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest
The authors declare no conflicts of interest.