Experimental Study on the Development Characteristics and Controlling Factors of Microscopic Organic Matter Pore and Fracture System in Shale

Research on microscopic pore and fracture system of shale is a hot spot in the field of unconventional petroleum geology. Micro- and nano-scale organic matter pores in shale play a vital role in the accumulation of hydrocarbons. Research on the types, evolution rules and controlling factors of organic matter pore-fracture system in shale reservoirs can provide scientific guidance for the prediction of shale “sweet spots”. In this paper, taking the shales from the Napo Formation, Oriente Basin, the Shahejie Formation, Zhanhua Sag, Bohai Bay Basin, and the Longmaxi and Wufeng Formations, Sichuan Basin as an example, the developmental characteristics of organic matter pore-fracture system were systematically studied using thin section, argon ion profiling scanning electron microscopy, X-ray diffraction, N2 adsorption and desorption, geochemistry experiments, and image processing technology. The types of shale organic matter pores were divided into kerogen-hosted pores, organic matter microfractures (intra-organic matter and organic matter edge microfractures), and asphalt pores (or intra-asphalt pores). The circumferences of organic pore were generally within 100 nm, and the areas of most pores were smaller than 1,000 nm2. Face rates of the organic pores were generally less than 1.5%, and the proportion of shale samples with a shape factor of 1 reached more than 50%. In addition, the deviation angles of organic matter pores at (0°, 45°) reached 90%, which showed that most of the organic matter pores tended to be oriented pores. The increase in the degree of thermal evolution provided driving force for the formation of circular pores in the organic matter. Internal factors (abundance of organic matter, kerogen types, and maturity) and external factors (diagenesis and mineral composition) controled the development of shale organic matter pores. Maturity, TOC content and inorganic minerals such as clay and pyrite content were positively correlated with the development of organic matter pores. However, brittle minerals caused a decrease in the face ratio of organic matter pores. Diagenetic compaction caused the organic matter pores being deformed or eventually disappeared. This research can provide scientific guidance for the high-efficiency exploration of hydrocarbons in shale.


INTRODUCTION
With the successful development of the North American shale gas industry, research on the microscopic pore structures in shale has become an emerging hot topic (Curtis, 2002;Zou et al., 2010;Dong et al., 2012;Santosh and Feng., 2020). According to the previous studies, unconventional gas shale has special reservoir structures and hydrocarbon accumulation mechanism (Zhang et al., 2004;Dong et al., 2009). The microscopic pore structures of shale affected its hydrocarbon storage capacities significantly Ko et al., 2016;Sun et al., 2016;Nie et al., 2018). Nano-scale pores and microfractures were well developed in organic-rich shales. Since microfractures were mainly developed in organic matter, therefore, the evolution of microfractures was affected by the evolution of organic matter significantly (Ko et al., 2016;Sun et al., 2016). The developmental characteristics of organic pores in shales with different sedimentary systems, maturities and depths were quite different (Chalmers and Bustin, 2007;. At present, the quantitative characterization methods of pores in porous media mainly include image analysis, fluid injection, nuclear magnetic resonance (NMR) and three-dimensional CT scanning methods . Nano-scale pores in shales can be observed through image analysis (SEM, TEM, etc.) (Reed et al., 2007;. However, the image analysis method can only be used to observe the features of pores on the surface of the sample. On the contrary, fluid injection (MIP, nitrogen and carbon dioxide adsorption etc.), NMR, and nano-CT methods have effectively overcame this deficiency. Different fluid injection measurements can be applied to satisfy requirements of different pore size (Chen et al., 2013;Yang et al., 2013;Gai et al., 2016;Zhao et al., 2020). NMR and nano-CT tests can be used to characterize the morphology, connectivity, pore size, and distribution of pores in shale reservoirs (Chen et al., 2012;Li et al., 2012;Chen, 2014;Tian et al., 2016).
The previous have proposed some classification standards for shale micro-pores (Loucks et al., 2009;Slatt and O'Brien, 2011;Yu, 2013). For example, according to the scheme proposed by Loucks et al. (2009), shale pores were classified into three categories: Intergranular, intragranular, and organic matter pores. However, there was no more detailed classification scheme for the types of organic pores. For organic-rich shales, organic pores play a vital role in the total pore space of shale (Jarvie et al., 2007;Chalmers et al., 2009;Wang et al., 2018). Organic matter coexisted with inorganic minerals, and they can influence each other (Zhang et al., 1996). For example, montmorillonite has strong adsorption and ion exchange ability, it can be transformed into a clay-organic matter complex with a strong catalytic ability in the hydrocarbon generation process (Gao et al., 1990;Liu, 1995).
In addition, high organic matter contents, high maturity, and high-quality kerogen types all can promote the development of organic pores .
In this paper, taking the typical organic-rich shales in the Bohai Bay Basin, the Oriente Basin and the Sichuan Basin as an example, the developmental characteristics of organic matter pore-fracture system were systematically studied using thin section, argon ion profiling scanning electron microscopy, X-ray diffraction, N 2 adsorption and desorption, geochemistry experiments and image processing technology. This research can provide scientific guidance for the extraction of reservoir parameters, gas-bearing indicators and the rational deployment of exploration and development plans of shale reservoirs.

Materials
In this study, 21 groups of shale samples were collected from three regions and three layers, namely the Napo Formation, Oriente Basin, Ecuador, the Shahejie Formation, Zhanhua Sag, Bohai Bay Basin, and the Longmaxi and Wufeng Formations, Sichuan Basin ( Figure 1). These three types of shale samples were collected from different regions and different depositional environments. However, all these shale samples have good hydrocarbon generation conditions. Among them, the Napo Formation shales are stable shallow sea shelf deposits with high organic matter content and dominated by Type II kerogen (R o is greater than 0.5%); the Shahejie Formation shales are shallow lake-semideep lake deposits with large shale thickness and high organic matter abundance; the Longmaxi and Wufeng Formations are deeper shallow shelf facies deposits, with high organic matter abundance, large shale rock thickness and high maturity. The differences in sedimentary environment will cause differences in the composition, type and abundance of organic matter within the rock. Moreover, the thermal evolution degrees of shales in different regions are quite different, which may be the controlling factors that cause the differences in the development of organic matter pore systems.
A Quanta 200F field emission environmental scanning electron microscope was used to observe the microscopic pore structures of the shale samples. The samples were cut into blocks of 1.5 cm × 1 cm × 3 mm and subjected to argon ion polishing treatment. The image information of the organic pores on the shale surface was obtained by scanning electron microscope observations. Moreover, the composition of minerals were determined through energy spectrum analysis.
The mineral composition detections of shale samples were realized via a D8 ADVANCE X-ray diffractometer; and the XRD experiments were used to obtain the content of main minerals and clay minerals in the samples. In addition, a CS230 analyzer was used to obtain the organic carbon content of shale after removing inorganic carbon. Geochemical parameters such as pyrolysis peak temperature T max and hydrogen index HI were obtained using a Rock-Eval 6 Standard pyrolysis apparatus. The maturity of organic matter was obtained by the extraction of kerogen from shale combined with a vitrinite reflectance measuring instrument.
The isothermal N 2 adsorption experiments used a Quantachrome ® Autosorb-iQ2-MP instrument from China University of Petroleum (Beijing). This instrument was used to measure the volume, specific surface area, and diameter of shale pores. BET (Brunauer Emmett Teller) equation was used to calculated the pore surface area, while the pore volume and pore size distribution were obtained using the Barrett, Johner and Halenda (BJH) methods.
In this paper, ImageJ Processing Software was used to perform pore segmentation on SEM images after argon ion polishing. The software can perform threshold segmentation of organic matter pores, and the pore structure parameters such as the area shape coefficient and face ratio of the organic pores in the field of view can be obtained. Among them, the Feret diameter is the distance between the two furthest points of the organic pores. The two-FIGURE 1 | Location and geological background of the study area. Notes: The study area is located in the Sichuan Basin in China, the Orient Basin in Ecuador and the Zhanhua Sag in the Bohai Bay Basin in China (A). The shale samples in the Sichuan Basin (B) were collected from the Silurian Longmaxi and Wufeng Formations; the shale samples in the Orient Basin (C) were collected from the Cretaceous Napo Formation; and the shale samples in the Zhanhua Sag of the Bohai Bay Basin (D) were collected from the Paleogene Shahejie Formation.
Frontiers in Earth Science | www.frontiersin.org October 2021 | Volume 9 | Article 773960 dimensional morphological features of pores were used to describe the roughness of organic pores, which was defined as: Where ff is the shape factor, S is the pore area, and C is the pore circumference (Jiao, 2015). The cross-sectional shapes of different regular aggregates represent pores (Figure 2) (Yao and Zhao, 2010). At present, it is difficult to accurately extract the proportion of pores on organic matter. Thus, the percentage of organic matter pore area to the total area under the same view field was used to characterize the pore content.
In addition, the measured 135°circular organic matter pore angle was used as a statistical baseline. The difference between the measuring angle of Feret diameter and 135°was used to describe the directional properties of organic pores, which was defined as: Where G is the deviation angle and F is the measurement angle of Feret's diameter.

Petrology and Geochemical Parameters
The results of whole rock and clay content XRD experiments were shown in Table 1. Experiments showed that the black shale of the Shahejie Formation was dominated by clay minerals and calcareous minerals; among them, the content of calcareous minerals ranged from 28 to 82%; while the content of clay minerals was only 7-35%; the content of the Imon mixed layers in clay minerals was between 51 and 77%. The Napo Formation samples were mainly composed of clay minerals and quartz; the content of clay minerals was high, FIGURE 2 | Classification of pore shape factors (After Yao and Zhao, 2010, modified).  Figure 3A) and clay minerals ( Figure 3B) of the shale samples in the three regions has a certain difference. Shale samples from 15 wells in three study areas were selected to determine the abundance, type and maturity of organic matter. The results are shown in Table 2. The TOC content of the Napo Formation samples ranged from 1.69 to 2.15%; Type II kerogen (mainly Type II 2 ) was developed; R o ranged from 0.646 to 1.087% (mature stage). The TOC content of the Shahejie Formation samples ranged from 0.65 to 4.8%; Type I and II kerogen (II 1 and II 2 ) were developed; R o ranged from 0.524 to 3.314%. The Longmaxi and Wufeng Formation samples have a high TOC content ranging from 3.41 to 3.99%; Type I kerogen was developed; R o ranged from 3.08 to 3.26%. On the whole, the Longmaxi and Wufeng Formation shales have the strongest hydrocarbon-generating ability, followed by the Shahejie Formation, while the hydrocarbon-generating ability of the Napo Formation is relatively poor.

Pore Structure Parameters
Studies have shown that the higher the thermal evolution degree of organic matter, the more developed the organic matter pores (Milliken et al., 2013). Therefore, shale samples with different thermal evolution degrees were selected for comparative research on the development characteristics of organic pores, which can help us better understand the evolution mode of organic pores (Milliken et al., 2013). Data statistics and analysis of 29 representative images in 6 wells showed that the organic pore diameter ( Figure 4A) and Frontiers in Earth Science | www.frontiersin.org October 2021 | Volume 9 | Article 773960 6 circumference ( Figure 4B) of different shale samples generally within 100 nm. In this interval, the proportion of organic matter pores are all greater than 50%. Meanwhile, the statistical results of the pore area ( Figure 4C) showed that the area of most pores was smaller than 1,000 nm 2 , but there were also many pores with an area greater than 5,000 nm 2 . The difference in pore diameter, circumference and area content of organic matter pores is obviously affected by the degree of thermal evolution (Cui et al., 2012).
The face rate of the organic matter pores in the argon ion polished scanning electron microscope images was calculated ( Figures 4D-F). The results showed that the face rate of the Longmaxi and Wufeng Formation shales was relatively low, which was generally less than 1.5%. Except for the Lb-1 well, the shale samples of the other wells have a face rate of less than 1%. The shale samples from Well BS5 of Shahejie Formation have a significantly higher face rate than samples from the Longmaxi and Wufeng Formations. This may be due to the fact that the field of view is less, and there is a larger error (Loucks et al., 2009;Li et al., 2020).
The range of pore diameters of the shale samples measured by the nitrogen adsorption experiments is shown in Table 2. The pore diameters of the Napo shale samples ranged from 1.688 to 5.199 nm; the pore diameters of Longmaxi and Wufeng Formation (L, W) shales ranged from 1.432 to 3.794 nm, and their pore volume and specific surface area were larger than those of the other two regions; the pore diameters of the Shahejie Formation (S) shales ranged from 1.432 to 29.396 nm, and the pore diameter varies greatly.
The pore volume of all shale samples ranged from 0.001 to 0.019 cc/g. Among them, the samples of Longmaxi and Wufeng Formations have a larger average pore volume, ranging from 0.011 to 0.019 cc/g, while those of the Napo Formation have the smallest average pore volume, ranging from 0.001 to 0.011 cc/g. The test results of the specific surface area are similar to the pore volume. The samples of the Longmaxi and Wufeng Formations have the largest specific surface area, ranging from 15.24 to 18.085 m 2 /g, while the samples of the Napo Formation have the smallest specific surface area, ranging from 0.167 to 4.66 m 2 /g.

Developmental Characteristics of Organic Matter Pores
Observations under the microscope showed that organic pores have various shapes, including round, elliptical, triangular, slitshaped, honeycomb, beaded, and irregular ( Figure 5). The structural parameters of organic pores in shale can well characterize the development characteristics of pores. Moreover, the shape, connectivity, concentration and directional properties of pores all have a significant influence on the large-scale accumulation of shale gas and the formulation of development plans (Hao et al., 2006;Houben et al., 2013;Kuila et al., 2014;Wu et al., 2014;Wu et al., 2015;Jiao et al., 2017).
Micron and nano-scale organic pores are mostly round (shape factor 0.89-1), indicating that their initial shapes are usually round ( Figures 5A-E). The larger the pores, the more complex their morphology. Observations showed that there were few circular pores for the larger pores, and the morphology of these pores tended to be elliptical or irregular (Figures 5F-H).
As the maturity of shale increased, organic pores of various shapes were formed, which was also the main reason for the formation of irregular pores. In addition, the organic pores may merge together, and then became flower-like ( Figure 5E) or irregular pores with good connectivity ( Figure 5B). Some pores have even undergone severe deformation or collapse due to strong compaction ( Figure 5D).
Organic matter pores have strong heterogeneity. The more pores develop inside the rock, the better its connectivity ( Figure 5D). For low maturity shale, there were few organic pores, so the connectivity of pores and microfractures in the rock was poor. The spatial distribution of organic pores were classified into two categories: directional type and chaotic type. When the external pressure was high, the long axis direction of the organic pores was consistent with the extension direction affected by plastic deformation ( Figure 5E).
The distribution of pores in the same piece of organic particles was usually similar ( Figure 5A). However, there were also cases where even within the same organic matter, the size and development of pores were quite different ( Figure 5H). This was because, as the thermal evolution continued, the small pores became larger and merged into large pores ( Figure 5H). Then, the shape and development of pores became complex, and the heterogeneity of shale pore structures gradually became significant. In some microscopic views, large numbers of organic pores merged to form pore clusters, while in other microscopic views, isolated organic pores may be still the main ones. The symbiosis of organic matter and pyrite was easy to be observed in shale samples, and the pyrite intercrystalline pores and oval organic pores were mainly developed ( Figure 5I). In the shale sample of Well Lc-1 (2,483 m), it was also found that some extremely small organic pores appeared on the boundary between organic matter and titanium dioxide ( Figure 5J).
In addition to organic matter pores, organic matter microfractures were also very developed in shale samples. In this study, organic matter microfractures were further divided into intra-organic matter microfractures ( Figure 5K) and organic matter edge microfractures ( Figures 5L,M). Organic matter microfractures generally had a long strip or slit-shaped structure, and their formation was closely related to the diagenetic shrinkage of kerogen ( Figures 5N,O).
The statistical results showed that the proportion of shale samples with a shape factor of 1 reached more than 50%. Moreover, the proportion of pores with a shape factor of 0.89-1.00 was close to 10%, that is, the pores in the samples were mainly round or close to round. This showed that during the thermal evolution process, the pore morphology tended to develop into circular pores. The proportion of pores that tended to be square (0.60-0.89) was between 8.32 and 37.57%. While the proportion of organic pores with a shape factor of 0-0.60 was between 9.86 and 22.75%, and its content was similar to that of pores that tended to be square ( Figure 4E).
Frontiers in Earth Science | www.frontiersin.org October 2021 | Volume 9 | Article 773960 Frontiers in Earth Science | www.frontiersin.org October 2021 | Volume 9 | Article 773960 8 According to Figure 6B, for shale samples in the same area or in different areas, the deviation angle of the organic matter pores was mainly less than 45°. The proportion of shale samples with deviation angles less than 45°was absolutely dominant, with an average content of more than 80%. Oriented pores refer to pores that are round and had a deviation angle of 0°, and their proportion was generally greater than that of pores whose deviation angle was between 0°and 45°, but their proportions were relatively close ( Figure 4F).
Statistics found that the deviation angle of organic matter pores in shale reservoirs at (0°, 45°) can reach 90% ( Figures 4E,F). This showed that most of the organic matter pores tended to be oriented pores as a whole, and only part of the pores are randomly distributed. When the maturity of the organic matter was low, the organic matter pores were mainly small, round pores, and the degree of orientation was very good. However, as R o continued to increase, the disordered distribution of organic pores gradually became significant. When the shale reached a large depth, the organic pores will undergo strong compaction and deformation, and the strong directional environmental stress may improve the directional properties of the pores. The increase of brittle minerals can increase the anti-compaction ability of organic pores, thereby weakening the directional properties of pores (Chalmers et al., 2009).

Classification of Organic Matter Pores
The previous have conducted preliminary studies on the microscopic pore classification in shale. For example, Zhang et al. (2016) divided shale pore system into pores and microfractures, while shale pores were further divided into matrix pores (intercrystalline pores and intragranular pores) and organic pores. However, there is no further fine division for the types of shale organic matter pores.
In this paper, organic pores in shale were divided into three major categories and four sub-categories combining research in this paper and the classification schemes given by the previous (Jarvie et al., 2007;Loucks et al., 2009;Milliken et al., 2013;Yu., 2013;Lu et al., 2015;Zhang et al., 2016). The classification scheme in this study comprehensively considered the shapes, structures, location, origin and size of shale pores ( Table 3).
According to Table 3, organic pores were divided into organic matter pores, organic matter microfractures and asphalt pores. Organic matter pores were kerogen-hosted pores. They have round, elliptical and irregular structures, and were mainly formed in the main hydrocarbon generation period. In addition, organic matter microfractures were divided into intra-organic matter microfractures and organic matter edge microfractures, both of which were related to organic matter shrinkage and are mainly formed in shales with low maturity. Moreover, the second cracking of asphalt for hydrocarbon generation can also form a large number of organic pores, and the shape, size and occurrence of the pores were similar to those of kerogen-hosted pores ( Table 3).

Influence of Abundance and Maturity of Organic Matter on Development Degree of Organic Matter Pores
For the coupling relationship between organic pores and the maturity, it was found that the higher the shale maturity, the more organic pores develop. For example, in Figure 6A, the maturity of the shale sample from Well TP1-1 (1,586.2 m) is low, and its R o is 0.65%. Observations under the microscope showed that the organic pores of the shale sample from Well TP1-1were not developed. In contrast, for the shale sample from Well Lc-1 (2,477 m) in Figure 6B, the R o is 3.02%. Observations under the microscope showed that the organic pores were very developed.
The study also found that when the R o of the shale samples was equivalent to immature and low-maturity conditions, there were   Frontiers in Earth Science | www.frontiersin.org October 2021 | Volume 9 | Article 773960 10 few organic pores in the shale, but large numbers of organic matter edge microfractures existed in the shale, such as the Well Y120 sample shown in Figure 7A (3,570.50 m, R o 0.59%) and the Well Y117 sample shown in Figure 7B (3,401.93 m, R o 0.56%).
Under mature and high-maturity conditions, strong hydrocarbon generation will promote the formation of a large number of organic matter shrinkage microfractures and intraorganic matter microfractures, such as the Well BS5 sample shown in Figure 7C (4,764.45 m, R o 1.73%).
When the shale reached the over-mature stage, almost all the pores are intra-organic matter pores, and almost all organic matter microfractures disappear. Meanwhile, organic matter edge microfractures can be transformed into organic pores, such as the Well Lc sample shown in Figure 7D (2,477.00 m, R o 3.02%).
The enrichment degree of organic pores are affected by many factors (Pommer and Milliken, 2015). These factors can be divided into two categories: internal factors (kerogen type, abundance, and maturity) and external factors (diagenesis and mineral composition). The study showed that as the shale maturity increased, the pore volume and specific surface area of the shale samples in the study area had the same changing trend (Figure 8). However, the pore structure parameters were in a fluctuating state as the shale maturity increased, that is, they didn't show a single upward or downward changing trend. In addition, the morphology of the correlation curves between pore diameter and cumulative pore volume can also be divided into three regions with significant differences (Figure 9). Therefore, the use of TOC or R o alone cannot determine the changing trend of pore structure parameters well. This feature also reflected that the development characteristics of organic pores were not only controlled by a single factor (Li., 2021).
From Figure 10, it can be found that the relationship between TOC content and pore structure parameters of shale samples in different regions has a certain difference. The correlation of L-W sample data was very poor (Figures 10A,B), but when all Longmaxi and Wufeng Formation samples were integrated, the data showed a good linear relationship ( Figures 10C,D). For the samples from the Napo and Shahejie Formations dominated by inorganic pores, the correlation between specific surface area and pore volume and TOC was not obvious, and there were many abnormal points ( Figures 10A-D). The abundance of organic matter provides a material basis for the evolution of organic matter pores. Generally, the higher the abundance of organic matter, the more abundant the organic matter pores. However, organic matter pores can develop in large quantities only when the organic matter reached mature stage. Therefore, the poor correlation showed that TOC content was not the only factor that affects the pore structures.
The types of kerogen mainly affect the development degree and morphology of organic matter pores. Many scholars have studied the differential evolution of shale organic matter pores (Bernard et al., 2012;Dong et al., 2015;Ma et al., 2017). The pore evolution process of different types of kerogen was different. For example, Dong et al. (2012) proposed that Type I kerogen was a "shrinkage type" and Type III kerogen was a "porous type", and partial sapropel kerogen was much more conducive to the development of organic pores. Jarvie et al. (2007) found that Type II kerogen was more prone to produce organic pores than Type I and III kerogens. Some scholars have found that when the chloroform bitumen "A" in shale was extracted, the pores increased significantly. Therefore, the retained oil will occupy the pore space generated by the hydrocarbon generation of organic matter. The kerogen in the Zhanhua Sag of the study area was dominated by Type I kerogen. During the rapid oil generation stage, as the degree of thermal evolution increased, the specific surface area and pore volume didn't show an increasing trend but decreased gradually (Figure 7). This may be due to the fact that the generated oil directly filled the pores. However, when R o was high enough, different types of organic matter can generate a lot of pores, and then the kerogen type was no longer the main controlling factor of organic matter pores. The correlation between R o and pore structure parameters showed that the pore volume and specific surface area of Longmaxi and Wufeng Formation shales increased with the increase of R o ; however, this relationship was poor in Napo and Shahejie Formation shale samples ( Figures 10E,F). Moreover, the pore volume and specific surface area of the Longmaxi and Wufeng Formations were much larger than those of the other two regions. The organic matter abundance and R o of the Longmaxi and Wufeng Formation shale samples were high, and the R o of the samples were all above 2%, which belong to the high-mature-over-mature shale. High-matureover-mature shale mainly generated dry gas, and the escape of hydrocarbons will cause a significant increase in pore volume and specific surface area. The Napo and Shahejie Formation shale samples had a low degree of thermal evolution. In the immaturelow maturity stage, the organic matter pores were relatively underdeveloped, and the organic matter shrinkage microfracture was an important pore type.
Some samples have higher R o values. With the increase of R o , the organic pores in shale develop from underdeveloped to highly developed, from small size to large size, from organic microfractures to organic pores, and from good orientation to disorderly development. Some correlations in Figure 10 were poor, indicating that the influencing factors of the development of organic matter pores were complex, which were affected by the geological factors such as organic matter abundance, kerogen type, and different mineral composition.

Influence of Mineral Components and Diagenesis on Development Degree of Organic Matter Pores
Diagenesis and mineral composition also have a significant impact on the development of organic pores. Diagenesis has a great impact on the preservation of organic matter pores. Strong compaction can cause the organic pores to be compressed, deformed, elongated, or even disappeared (Hao et al., 2006;Jiao et al., 2017;Wang et al., 2020). According to Figure 11, the pore volume and specific surface area values of the Longmaxi and Wufeng Formation samples were inversely proportional to the content of brittle minerals ( Figures 11A,B), and directly proportional to the content of clay minerals ( Figures  11C,D); however, these relationships were not significant for the Napo Formation shale samples.
A weaker hydrodynamic environment was conducive to the growth and reproduction of organisms, and thus the organic matter components will be relatively enriched. However, the content of clay minerals in this kind of sedimentary environment will be higher, and the content of brittle minerals dominated by terrigenous clastics will be less. Therefore, the increase in the content of brittle minerals reflected a gradually stronger hydrodynamic condition, which was not conducive to the evolution and preservation of kerogen. Meanwhile, the clay minerals coexisting with organic matter had a catalytic effect on the maturation of organic matter, which in turn facilitated the evolution of kerogen-hosted pores and will cause the specific surface area and pore volume being greatly improved. Therefore, the development of organic pores was positive to the content of clay minerals, and it was negative to the content of brittle minerals ( Figure 11). It can also be seen from Figure 11 that the correlation between the pore structure parameters and mineral components of the high maturity shales (Longmaxi and Wufeng Formation shales) was very good; however, the correlations between these parameters was poor in low-maturity shales (Napo Formation shales). This was because the Napo Formation shales have low maturity and Frontiers in Earth Science | www.frontiersin.org October 2021 | Volume 9 | Article 773960 underdeveloped organic pores, so the controlling factors of pore structure parameters were much more complicated.

Development Mode of Organic Matter Pores
In summary, the development of pores in shale reservoirs are affected by internal and external factors. Internal factors include organic matter abundance, kerogen type and maturity; and external factors include mineral components (brittle minerals, clay minerals and FeS) and diagenesis. Organic carbon is the material basis of organic pores (Morozov et al., 2021). The higher the degree of thermal evolution of organic matter, the more developed the organic pores in shales (Wu et al., 2021). The kerogen type of organic matter is an important factor controlling the differential evolution of organic matter pores. Under the same or similar TOC and R o conditions, the specific surface area and pore volume of the microscopic pores of shales dominated by Type II kerogen were better than those of Type I kerogen, however, the organic pores in shale dominated by Type III kerogen were relatively underdeveloped ( Figure 12). Organic matter maturity provides the driving force for the evolution of organic matter pores. Under the same TOC content, the higher the maturity, the more organic matter pores will develop.
However, under immature and low-maturity conditions, organic matter microfractures were developed in all types of kerogen. These microfractures were mainly organic-edge shrinkage microfractures, and organic pores were relatively underdeveloped ( Figures 7A,B). With the increase of R o , during the early-to-mid mid-diagenesis period, the microfractures developed at the edge of organic matter gradually transitioned into the interior of organic matter. However, in the mid-to-late mid-diagenesis period, the organic FIGURE 11 | Relationship between pore structure parameters and content of different mineral components of the shale samples. Notes: (A) Correlation between brittle mineral content and specific surface area; (B) Correlation between brittle mineral content and pore volume; (C) Correlation between clay mineral content and specific surface area; (D) Correlation between clay mineral content and pore volume; "Napo"-Napo Formation shales; "L-W"-Longmaxi and Wufeng Formation shales; "L"-Longmaxi Formation shale data from Ji et al., 2016. matter was in the mature-over-mature stage, a large number of organic matter pores were formed and the organic matter microfractures no longer developed (Figure 12).
Diagenesis and different mineral compositions are external factors that affect the development of organic pores. Strong compaction can cause the organic pores deformed or even disappeared. Clay minerals and pyrite can promote the development of organic pores, while brittle minerals (such as quartz, feldspar and carbonate minerals) can inhibit the development of organic pores.
Overall, the development mode of organic pores in shale reservoirs is the coupling results of multiple geological factors ( Figure 12). TOC and kerogen type are the most direct controlling factors. Under the premise of high degree of thermal degree, the higher the abundance of organic matter, the more organic pores developed; and under the same thermal degree, the development of organic pores satisfies: Type II kerogen > Type I kerogen > Type III kerogen. Thermal evolution or R o is the source of organic matter pores. Only when the organic matter reaches the mature to high maturity range, the organic matter pores will develop in large numbers in the shales and have good connectivity.
CONCLUSION 1) In this paper, taking the typical shales from the Bohai Bay Basin, the Oriente Basin, and the Sichuan Basin as an example, argon ion profiling scanning electron microscopy, N 2 adsorption and desorption, geochemistry experiments and ImageJ software image processing technology were used to study the types and evolution of micro-and nano-scale organic matter pores in shale reservoirs.
2) The circumferences of organic pore were generally within 100 nm, and the areas of most pores were smaller than 1,000 nm 2 . Face rates of the organic pores were generally less than 1.5%, and the proportion of shale samples with a shape factor of 1 reached more than 50%. In addition, the deviation angles of organic matter pores at (0°, 45°) can reach 90%, which showed that most of the organic matter pores tended to be oriented pores.
3) The organic pores of shale were divided into kerogen-hosted pores, organic matter microfractures (intra-organic matter and organic matter edge microfractures) and asphalt pores (intra-asphalt pores). As the degree of thermal evolution increased, the organic pores gradually tend to be round, and the directional properties of the pores was strengthened. 4) The development of shale organic pores was controlled by internal factors (abundance of organic matter, kerogen type and maturity) and external factors (diagenesis and mineral composition). Maturity promoted the development of organic pores; TOC content was positively correlated with the development of organic pores; and there were significant differences in the degree of development of organic pores in shale organic matter of different sedimentary conditions. 5) Clay minerals and pyrite can promote the development of organic pores, while brittle minerals can cause a decrease in the face ratio of organic pores. The diagenetic compaction will cause the organic pores to deform or eventually disappear. Finally, the development mode of shale organic matter pores was clarified combined with the development characteristics and influencing factors of organic matter pores.
FIGURE 12 | Influence of thermal maturity-related parameters on the development degree of shale organic matter pores. With the increase of R o , during the earlyto-mid mid-diagenesis period, the microfractures were developed. In the mid-to-late mid-diagenesis period, the organic matter was in the mature-over-mature stage, a large number of organic matter pores were formed and the organic matter microfractures disappeared gradually.
Frontiers in Earth Science | www.frontiersin.org October 2021 | Volume 9 | Article 773960 6) Through this study, the development characteristics, evolution process and controlling factors of organic pores in shales formed in different sedimentary environments were compared. This research can provide scientific guidance for the extraction of evaluation parameters and gas-bearing indicators and the rational deployment of exploration and development plans of shale reservoirs.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
QZ is responsible for the idea and writing of this article, and ZL, CL, XZ, RS, HT, KW, and ZS are responsible for the experimental part.