Heterogeneity Characteristics of Lacustrine Shale Oil Reservoir Under the Control of Lithofacies: A Case Study of the Dongyuemiao Member of Jurassic Ziliujing Formation, Sichuan Basin

The lacustrine shale in the Dongyuemiao Member of the Fuling area, Sichuan Basin, is widely distributed and has huge shale oil resource potential. It is one of the important replacement areas for shale oil exploration in China. To investigate the key shale oil evaluation well, Well FY10, in the Fuling area, X-ray diffraction (XRD) mineral analysis, Rock-Eval, argon ion polishing-scanning electron microscope (SEM), Mercury injection capillary pressure (MICP), and low pressure nitrogen adsorption were launched to determine the heterogeneity of the pore system in the lacustrine shale of the Dongyuemiao Member. The mineral composition exhibits a high degree of heterogeneity, and the shale can be divided into two main lithofacies: argillaceous shale and mixed shale. The porosity ranges from 2.95 to 8.43%, and the permeability ranges from 0.05 to 1.07 × 10−3 μm2. The physical properties of mixed shale are obviously better than those of argillaceous shale. Inorganic mineral pores, such as linear pores between clay minerals and calcite dissolution pores, are mainly developed, while a small amount of organic pores can be observed. The average total pore volume (Vp) is 0.038 ml/g with an average specific surface area of 5.38 m2/g. Mesopores provide the main Vp (average 61.72%), and micropores provide mostly specific surface area. TOC imposes a strong controlling effect on the development of micropores. Clay minerals are the main contributors to mesopores and macropores. The organic-inorganic interaction during the process of diagenesis and hydrocarbon generation controls the formation of shale pore systems.


INTRODUCTION
With the commercial exploitation of shale oil in the Bakken Formation in the Williston Basin, the Eagle Ford Formation in Southern Texas, and the Barnett Formation in the Fort Worth Basin, shale oil has altered the American energy structure and imposes a profound impact on the world energy pattern (Hill et al., 2007;Zhang et al., 2017;Saputra et al., 2019). In recent years, with the rapid development of China's economy, the demand for oil and gas has increased substantially, and the degree of dependence on foreign oil and natural gas has continued to increase. There is an urgent requirement to find new types of oil and gas resources, achieve theoretical and technological breakthroughs, and form a contribution to domestic oil and gas production (Jia et al., 2012;Zou et al., 2020). Among the known new resource types, shale oil is one of the most realistic exploration and development targets (Zou et al., 2013;Jin et al., 2019). Shale oil refers to the oil contained in the organic-rich shale strata (including the thin interbedded tight carbonate and clastic rock layers in the shale). There is no production capacity under natural conditions, and unconventional technology is required to realize economies of scale exploitation (Jin et al., 2019).
Organic rich shale is widely distributed in China, which not only provides abundant conventional oil and gas resources, but also contains plentiful shale oil and gas resources (Jin et al., 2019;Zou et al., 2019). Many breakthroughs have been achieved in the Permian Lucaogou Formation in the Junggar Basin, Triassic Yanchang Formation in the Ordos Basin, Cretaceous Qingshankou Formation in the Songliao Basin, middle-lower Jurassic formations in the Sichuan Basin, and Paleogene Shahejie Formation in the Bohai Bay Basin, and industrial production has been realized in some regions Liu et al., 2019;Fu et al., 2021;Zhang et al., 2021).
Lacustrine shale is characterized by very developed strata, wide distribution, and massive exploration potential (Jia et al., 2012;Jin et al., 2019;Zou et al., 2019;. However, there are still some unresolved problems such as strong heterogeneity, an unclear occurrence mechanism, and the main controlling factors of enrichment (Guo et al., 2016;Li et al., 2020;.
Based on the comprehensive study of the geological conditions of Dongyuemiao shale of Jurassic Ziliujing Formation in Sichuan Basin, the heterogeneity of lacustrine shale under lithofacies constraints in Fuling area is studied, in order to provide support for promoting the exploration and development process of lacustrine shale oil and gas in the Sichuan Basin and enrich the geological theory of lacustrine shale oil and gas.

GEOLOGICAL SETTING
The Jurassic formation in the Sichuan Basin is a set of important petroleum-bearing strata containing both oil and gas (Liu et al., 2019). During the early and middle Jurassic, most of the Sichuan Basin was composed of shore shallow lacustrine and semi-deep lacustrine sediments (Li and He, 2014). Vertically, three sets of organic rich shales are developed in the Dongyuemiao Member, Da'anzhai Member of Ziliujing formation, and Lianggaoshan Formation (Zheng et al., 2013). During the deposition period of Dongyuemiao Member, the sedimentary center of Sichuan Basin is mainly located around Liangping-Guang'an-Fuling area in the east of Sichuan Province. The black shale exhibits a relatively high abundance of organic matter and a wide distribution area (Shu et al., 2021).
The Fuling area is located in the eastern of the Sichuan Basin ( Figure 1). During the Dongyuemiao sedimentary period, the Fuling area experienced a large-scale lake transgression, depositing a set of shore shallow lacustrine to semi-deep lacustrine black shale interbedded with thin shell limestone, with regional thickness of 60-70 m. The Dongyuemiao Member in the area can be further divided into three sub members: the first sub member is black shale interbedded with thin layers of shell limestone, with a thickness of 28 m. The second sub member is composed of black argillaceous limestone with a thickness of 6.5 m. The third sub member is dark mudstone with a thickness of 27 m. In general, the lacustrine water body of the Dongyuemiao Member in the study area has experienced a change process from shallow to deep to shallow.
In March 2020, the FY10 well deployed by Sinopec Jianghan Oilfield Company was successfully fractured and became the first production well of shale oil and gas in the area (Figure 1). It obtained a daily gas production of 5.58 × 10 4 m 3 and oil production of 17.6 m 3 . The total organic carbon (TOC) content of the favorable interval of Well FY10 is 2.11%, the porosity is 6.11%, and the permeability is 0.127 × 10 −3 μm 2 , showing good reservoir performance in general. The Ro of the shale in Dongyuemiao Member is between 1.23 and 2.09%. The organic matter has evolved to a high mature stage, that is, the condensate to wet gas stage (Shu et al., 2021).

SAMPLES AND METHODS
A total of 13 shale samples were carefully selected from the shale strata in the Dongyuemiao Member of Well FY10. Each sample is divided into several parts for different analyses.

Experimental Method
The analyses in this study were all launched at the Wuxi Research Institute of Sinopec Petroleum Exploration and Development Research Institute. First, 13 samples were analyzed for TOC content, rock pyrolysis, and XRD mineral composition. The analyses process comply with the operating steps outlined in the National Standard (SY/T) 14506.28-2010 and Chinese Oil and Gas Industry Standard (GB/T) 19145-2003. Nine of the samples were selected for thin slice identification analysis. Meanwhile, eight shale samples were launched for argon ion polishing-scanning electron microscopy analysis to observe the development characteristics of pores.
To quantitatively characterize the full-scale pore size distribution characteristics of shale, eight samples were selected for Mercury injection capillary pressure (MICP) and low pressure nitrogen adsorption, under the experimental procedure proposed in the National Standard (GB/T) 21650.1-2008 and (GB/T) 19587-2017 (Li et al., , 2017. Theoretically speaking, the measurement range of the nitrogen adsorption method is from 0.35 to 500 nm. Due to the limitations of various formulas and models, the nitrogen adsorption method only has high measurement accuracy for pores below 50 nm, while the precision of the test results for large pores is reduced. The measurement range of the MICP method is from 6.5 nm to 63 μm. Because high pressure can damage the shale structure and then affect the results, the MICP has high measurement accuracy for pores above 50 nm (Li et al., 2017). In general, only by combining the two methods can it get the continuous pore distribution of shale. The method used in this paper is to combine the nitrogen adsorption data below 50 nm with the MICP data above 50 nm. In the meantime, the Brunauer Emmett Teller (BET) method was applied to calculate the specific surface area.

Lithofacy Classification Method
According to the quantitative analysis of the total mineral composition of lacustrine shale, the two-level lithofacies' division is carried out on the basis of the graphic method of clay minerals, QMF (quartz, mica, and feldspar), and carbonate minerals as three terminal elements. In the first-level division process, the traditional petrological classification method is respected, and the name of vary rock is determined with 50% as the boundary. At the three ends of the triangle diagram, the shale area (I), sandstone area (III), and limestone area (IV) are determined, respectively. In the meantime, the middle area where the content of the three components does not exceed 50% is named the mixed rock area (II) (Figure 2A). On this basis for the I, III, and IV areas, taking the mineral content of 25% as the boundary, the secondary division is carried out. For example, the shale area (I) is further divided into three secondary areas: argillaceous shale, silty shale, and carbonate shale. For the mixed rock area, it is named mixed shale ( Figure 2B).
On the basis of the above-mentioned three-terminal division, it is also necessary to identify the structure (including grain structure and morphology), genesis, and related sedimentary characteristics of the two non-clay minerals, quartz and calcite in shale, to improve the division of lithofacies.

Mineral Contents
The detailed results of XRD mineral analysis are shown in Table 1. The mineral composition of the shale in the Dongyuemiao Member is highly heterogeneous, mainly including clay minerals, quartz, calcite, and a small amount of feldspar, dolomite, and pyrite. The content of clay minerals ranges from 25.9 to 69.0%, with an average value of 55.0%. The quartz content ranges from 17.8 to 29.4%, with an average value of 23.5%. The calcite content is between 0 and 51.2%, with an average value of 14.0% ( Table 1). The three primary minerals accounted for 92.4% of the contents of the total minerals. Except for the relatively stable quartz content, the contents of clay minerals and calcite are both highly heterogeneous and fluctuate with each other.

Types of Lithofacies
On account of the ternary diagram of the lithofacies classification, two lithofacies of the Dongyuemiao shale can be identified, namely, mixed shale and silty shale ( Figure 2). But for quartz, thin-section observation shows that the grain size structure of quartz in the shale is different from that of normal sandstone. Taking the two samples of the same silty shale (FY-74 and FY-77) as an example, the clay mineral contents are greater than 50%, and the quartz contents are 21 and 25.5%, respectively, which are all argillaceous structures, and it is difficult to identify obvious quartz particles (particle size <0.039 mm) ( Figures 3A,B).
Obviously, it is inappropriate to define it as silty shale, and it is more suitable to identify it as argillaceous shale. For the calcite mineral in mixed shale, core observation and thin-section identification analysis exhibit that it is mainly the product of calcification and cementation of the original biological shell, which is caused by bio-chemical sedimentation. Under the microscope, two-stage calcite cementation can be seen in individual shells ( Figure 3C), and there exist different degrees of recrystallization ( Figure 3D), but in most cases the original biological structure is still better.
By determining the structure and genesis of minerals, the definition of certain rocks can be effectively corrected, and the lithofacies of the same interval in the triangle diagram of mineral quantitative analysis can be distinguished. According to the results of XRD mineral analysis and thin slice observation, the Dongyuemiao shale of Well FY10 mainly develops argillaceous shale and mixed shale. The vertical lithofacies change and the water depth fluctuation have a certain synergism.

Organic Geochemical and Physical Characteristics
Total organic carbon (TOC) is usually the main indicator for evaluating the abundance of source rocks Liu et al., 2019;Feng et al., 2020). The TOC content of the shale in the Dongyuemiao shale of Well FY10 ranges from 0.77 to 3.13%, with an average value of 1.67% ( Table 2). Rock pyrolysis hydrocarbons are composed of adsorbed hydrocarbons (S 1 ) and pyrolysis hydrocarbons (S 2 ). Their contents vary with the abundance and type of organic matter, and are direct indicators for  evaluating the quality of source rocks (Huang et al., 2020). S 1 ranges from 0.48 to 1.10 mg/g, with an average value of 0.70 mg/g. Hydrocarbon generation potential (PG S 1 +S 2 ) is an important indicator to characterize the abundance of organic matter (Huang et al., 2020). With the increase of maturity, the PG will obviously decrease. The PG of the Dongyuemiao shale of Well FY10 ranges from 1.16 to 3.08 mg/g, with an average value of 1.91 mg/g, which makes it a generally medium source rock. PG and TOC exhibit an obviously positive relationship, but there exist certain differences between argillaceous shale and mixed shale ( Figure 4). The range  To express the oiliness of shale, Jarvie (2012); Jarvie (2014) proposed the Oil Saturation Index (OSI), that is, the ratio of the TOC of the shale to its pyrolysis S 1 (S 1 /TOC). The OSI of shale in Dongyuemiao Member of Well FY10 ranges from 30.85 to 62.34 mg/gTOC, with an average value of 43.92 mg/ gTOC, and there exist no significant differences between different types of shale. The analysis suggests that S 1 reflects the remaining oil in the source rock, but is affected by sampling conditions, evolution stages, and oil and gas migration. The shale in the study area has evolved to a high maturity stage of the condensate to wet gas stage, so the light hydrocarbon components are easily lost, which has a certain impact on the results. The δ 13 C of kerogen ranges from −25.0‰ to −30.7‰, which indicates that type II kerogen is dominant in organic matter.
Through physical experiments, the porosity and permeability of the Dongyuemiao shale can be obtained. The porosity ranges from 2.95 to 8.43%, with an average value of 6.84%; the permeability ranges from 0.05 to 1.07 × 10 −3 μm 2 , with an average value of 0.26 × 10 −3 μm 2 ( Table 2). The physical properties of mixed shale are obviously better than those of argillaceous shale, but the relation between porosity and permeability of mixed shale is unclear, and the porosity and permeability of argillaceous shale exhibit an obvious positive correlation ( Figure 5).

Pore Types and Full-Scale Pore Distribution
A combination of argon ion polishing-SEM, low-temperature nitrogen adsorption, and MICP was launched to carry out qualitative observation and quantitative analysis of shale micro-pores. To fully consider the contact and genetic relationship between inorganic mineral particles, organic matter, and pores, the pores of the Dongyuemiao shale in the study area are finally divided into organic pores and inorganic pores. Inorganic pores can be further divided into clay mineral pores and calcite-related pores.

Pore Morphology From Scanning Electron Microscope Inorganic Pores
The development of inorganic pores, in essence, is similar to conventional reservoirs. There are two main formation mechanisms for the development of inorganic pores. One is the residual pores developed after the early intergranular pores undergo diagenetic compaction and cementation, and the other is the dissolution pores formed by the later dissolution modification (which can occur between and within grains); in addition, some pores also develop during mineral diagenetic alteration. The inorganic mineral pores in the shale of Dongyuemiao Member are mainly calcite dissolved pores and linear pores between clay minerals.
Compared with the marine Wufeng-Longmaxi shale (Dai et al., 2014;Feng et al., 2016;Nie et al., 2018), the Jurassic lacustrine shale has a higher content of clay minerals (>50%), and the average content of clay minerals in the Dongyuemiao shale of Well FY 10 can reach 55.0%. SEM observation reveals that the shale of Dongyuemiao Member widely develop linear micropores related to clay minerals ( Figures 6A,B), that mainly exist between illite flakes or illite and mineral particles. Linear pores vary in length and width, and some of the pores are filled with solid asphalt. Therefore, the organic matter content in different samples has a significant effect on the development of linear pores. This type of pore is mainly composed of primary pores between early clay minerals. With the continuous increase of burial depth, due to the weak anti-compaction ability between clay minerals, the pores rapidly shrink to form residual linear pores. This type of pore basically develops along the cleavage direction of clay minerals, forming slit-shaped or linear pores, which can extend several microns.
Calcite dissolved pores are developed in the Dongyuemiao shale of Well FY10. Calcite is mainly composed of a crystalline granular structure, and the development of dissolution pores is various, mainly in three forms ( Figures 6C-E). The first is the dissolution pores at the edge of calcite grains with irregular shapes. The second is the dissolution of macropores in the calcite grains, with various shapes, having clay minerals or dissolution residual calcite in the pores. The third is the dissolution needle-like pores in the calcite grains, and the individual is significantly smaller than the formers. The  formation of calcite-related dissolved pores may be mainly related to the dissolution of organic acids during hydrocarbon generation.

Organic Pores
Observations of most shale samples from the Dongyuemiao Member reveal that the organic pores are well developed in the solid asphalt filled within inorganic minerals. The pore shapes vary, and include oval, irregular, slit, and honeycomb shapes ( Figures 6F-H). The pore diameter is mainly in nanometers, but the size is not uniform, from tens of nanometers to hundreds of nanometers (the long axis of the pore size) are developed, and some samples show the development of micron-scale organic pores. The development of organic pores is mainly controlled by organic matter components. The organic pores in the solid asphalt are generally developed. For vitrinite, most of them did not develop organic pores except for a few which were modified by bacteria ( Figure 6I).

Pore-Throat Structure Parameters and Pore Size Distribution
The entry pressure of the lacustrine shale is 0.46 MPa on average. The corresponding maximum pore-throat radius is 580.61 nm. The average median pressure reaches 20.88 MPa. The corresponding median radius is 8.09 nm. There exists no obvious difference in pore-throat structure parameters of different lithofacies shale.
The total pore volume (V p ) distribution range of Dongyuemiao shale is 0.0097-0.0163 cm 3 /g, the average value is 0.038 cm 3 /g, and the specific surface area ranges from 2.93 to 9.65 m 2 /g, with the average value of 5.38 m 2 /g ( Table 3). According to the IUPAC pore classification standard, pores are classified into micropores (<2.0 nm), mesopores (2-50 nm), and macropores (>50 nm) according to their size (Loucks et al., 2012;Thommes et al., 2015). In general, the average proportions of micropores (<2 nm), mesopores (2-50 nm), and macropores (>50 nm) are 15.35, 61.72, and 22.93%, respectively. It can be seen that the proportion of mesopores is dominant and the proportion of macropores is also relatively high. The shale of the marine Wufeng-Longmaxi Formation is dominated by micropores and mesopores, the contents of which are similar, and together they can account for more than 95% of the total V p Feng et al., 2018Feng et al., , 2019. In general, compared with the marine shale, the proportion of meso-macropores in the Dongyuemiao Member of the lacustrine shale is higher, which is very beneficial for improving the reservoir performance of shale. From the full pore size distribution of shale, it can be detected that the distribution of different samples still has certain differences. The pore size distribution of the argillaceous shale is similar to that of the mixed shale. An obvious main peak of the pore size distribution exits and the main body is in the range of mesopores (2-20 nm) ( Figures 7A,B). The pores of the two types of shale are mainly mesopores, which accounts for more than 60% on average. The proportion of micropores in mixed shale is slightly lower than that of macropores, with average values of 18.28 and 19.17%, respectively ( Figure 7C; Table 3). The proportion of macropores in argillaceous shale is significantly higher than that of micropores, with average values of 25.18 and 13.60%, respectively ( Figure 7D; Table 3). This difference in pore size distribution may be due to the different organic and inorganic mineral compositions of the two types of shale.
It can be detected from the intersection of pore volume and specific surface area that the total V p and specific surface area have a certain positive correlation, but the correlation is not strong ( Figure 8A). The relationship between the corresponding V p of pores in different scales and specific surface areas is very different. The V p of micropores and specific surface areas show a significant positive correlation, R 2 can reach 0.9787 ( Figure 8B), while V p of mesopores and macropores have no obvious correlation with specific surface area ( Figures 8C,D), indicating that micropores are the main contributor to the specific surface area, which is consistent with marine shale. Different types of shale also exhibit obvious differences. The V p of micropores of mixed shale is only positively correlated with the specific surface area, while the V p of micropores and mesopores of argillaceous shale exhibit an obvious positive correlation with specific surface area, indicating that the development of pores in argillaceous shale imposes a stronger controlling effect on surface area than mixed shale, which achieves more obvious regularity characteristics.

DISCUSSION
To probe the differences between shale porosity and pore size distribution of different lithofacies, based on the material composition of shale, the reason for these differences is analyzed comprehensively from the composition of organic matter and inorganic mineral. The Relationship Between Total Organic Carbon and Pore Structure Parameters As an important part of the shale pore system, organic pores are formed during the process of hydrocarbon generation. The development of organic pores is of great significance to evaluate the gas storage capacity of shale intervals. Therefore, the formation mechanism of organic pores has always been a hot issue in the field of shale reservoirs. The development of organic pores in the Dongyuemiao shale is unique, which exhibits that some organic matter in the same sample is compact and uniform and develops no internal pores, while some organic matter has relatively developed internal pores. The organic matter with undeveloped pores is mostly in the shape of a strip or block. It is a microscopic component formed by the humification and gelation of the lignocellulosic tissue of higher plants (mainly vitrinite). The organic matter with organic pores is mostly unfixed, controlled by the original pore carrier (surrounding inorganic minerals), and is solid asphalt. As shown in Figure 9, the V p of micropores and specific surface areas are positively correlated with TOC content, indicating that micropores are mainly developed in organic matter. TOC content seems to be the main parameter controlling micropores' development. However, there exists no obvious correlation between the V p of mesopores and macropores and the TOC content, which indicates that the pores in organic matter are not the main contributors of mesopores and macropores.
In general, TOC content was positively correlated with micropore volume, indicating that organic matter could be used as the main index of micropore development. The higher the TOC content is, the more micropores and the larger the specific surface area in shale, which is consistent with the previous research results of marine shale (Loucks et al., 2012;Feng et al., 2019).

The Relationship Between Mineral Composition and Pore Structure Parameters
There are mainly two types of inorganic mineral pores in the Dongyuemiao shale of Well FY10, namely linear pores between clay minerals and calcite dissolved pores (Figure 6), which are closely related to the mineral composition of the shale. The average content of clay minerals in the Dongyuemiao shale is about 55%. The high clay content leads to weak compaction resistance, and the clay minerals are mostly oriented and distributed in layers, with residual clay minerals interlayers. The interlayer linear micropores of residual clay minerals are the main pores of lacustrine shale. Although the intergranular pores between clay minerals and mineral particles are developed,  Frontiers in Earth Science | www.frontiersin.org September 2021 | Volume 9 | Article 736544 10 the number of pores is relatively small, and most of the large pores have been filled by solid asphalt. According to the statistics of the relationship between clay minerals and pore volume, it can be detected that with the increase of clay minerals, the V p of mesopores and macropores tends to become larger, and the correlation coefficient with macropores is higher, showing a better relevance (Figure 10). It is shown that its contribution to micropores and mesopores is relatively apparent, and the linear pores between clay minerals provide the main reservoir space for the shale in the study area.
The calcite of Dongyuemiao Member is dominated by a crystalline granular structure with various dissolved pores, but there is no obvious regularity in its development. It is speculated that the formation of calcite-related dissolved pores mainly related to the organic acid dissolution during the hydrocarbon generation process. According to the statistics of the relationship between calcite and pore volume, it can be detected that, with the increase of calcite content, the V p of mesopores and macropores tends to decrease, but there is no obvious correlation with the V p of micropores. Although calcite may develop various types of pores, its contribution to pores in different sizes is relatively limited, so calcite dissolved pore is not the main contributor of reservoir space of lacustrine shale (Figure 11).

Diagenesis-Hydrocarbon Generation Controls the Formation and Evolution of Pores
The pore development and evolution of shale is complex, and is closely related to the compaction during the diagenesis process, the hydrocarbon generation and expulsion of organic matter, and the precipitation and dissolution of carbonate minerals. In addition to various inorganic minerals (clay minerals, quartz,  Frontiers in Earth Science | www.frontiersin.org September 2021 | Volume 9 | Article 736544 feldspar, calcite, etc.), shale is also rich in organic matter, and organic pores that contribute to the accumulation of oil and gas are found inside. With the increase of the buried depth of the shale, the organic matter also evolves with the increase of the temperature in the process of the inorganic mineral diagenesis. Inorganic pores formed by the diagenetic evolution of inorganic minerals coexist with the pores formed by the hydrocarbon generation of organic matter, which reveals the formation and evolution of diagenesis, hydrocarbon generation, and pores as an integrated whole. The shale in the Dongyuemiao Member of the Fuling area has entered the middle diagenetic stage B, with strong compaction; a large amount of montmorillonite was transformed into lamellar illite, the intergranular pores and intragranular pores have become smaller and the number decreased, but a large number of mesopores and macropores still formed, which are the main contributors to the shale reservoir space. The maturity (Ro) of shale is 1.5%, and it has entered a high evolutionary stage. Organic pores begin to form during the secondary pyrolysis of petroleum. Observation under the SEM exhibits that organic pores, mainly micropores, are widely developed in the asphalt. To some extent, it increases the porosity and contributes to the total porosity. Therefore, the difference in the reservoir space of different lithofacies is only a characteristic phenomenon. The organic-inorganic interaction controls the formation of shale pore systems in different lithofacies during diagenesishydrocarbon generation.

CONCLUSION
According to the mineral contents and thin slice analysis, the Dongyuemiao shale is divided into two main shale lithofacies: argillaceous shale and mixed shale. The porosity ranges from 2.95 to 8.43%, and the permeability ranges from 0.05 to 1.07 × 10 −3 μm 2 . The physical properties of the mixed shale are obviously better than those of the argillaceous shale.
The shale in the Dongyuemiao Member is mainly composed of inorganic mineral pores, such as linear pores, between clay minerals and calcite dissolved pores, while a few organic pores are also developed. It is characterized by a high proportion of macropores and a predominance of mesopores. The mesopores provide the main pore volume (61.72% on average), while the micropores provide the main specific surface area.
The TOC content of the shale in the Dongyuemiao Member imposes the main parameter controlling the development of micropores. Clay minerals influence the development of mesopores and macropores, and the combined interaction of diagenesis and hydrocarbon generation controls the formation of pore systems in shale of different lithofacies.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.