Microscopic pore-fracture combination types and their impact on exploration and development: in case of deep shale gas reservoirs in the Luzhou area of the Sichuan Basin

Abstract

Micro-pores and micro-fractures constitute the most vital reservoir spaces and seepage pathways for shale gas. Conducting an in-depth study of pore-fracture combination types within shale gas reservoirs is beneficial for understanding shale gas enrichment and exploitation. In this study, we conducted a joint experiment involving MICP and WM impregnation on the main gas-producing layer of Longmaxi formation shale reservoir from five shale gas wells in the Luzhou area. This study achieved the visualization of MICP experiments, enabling a quantitative classification of pore volume and micro-fracture volume within the shale reservoir. A Q_P-F parameter system for classifying pore-fracture combination types has been established to evaluate pore-fracture connectivity, which are categorized into three types. The strongly connected pore-fracture type (SCPFT) shale gas reservoir in the study area has the characteristics of high porosity, high gas content and strong pore-fracture connectivity, which is the preferred target for shale gas exploration and development. Moderately connected pore-fracture type (MCPFT) and weakly connected pore-fracture type (WCPFT) shale gas reservoirs exhibit significant variations in porosity and gas content, coupled with weak pore-fracture connectivity. Despite some of these reservoirs possessing high porosity and gas content, they fail to achieve high production rates. Therefore, during the development of these two types of shale gas reservoirs, it is imperative to enhance and optimize fracturing techniques to create a more intricate fracture network structure, thereby enabling the release of shale gas trapped in poorly connected pores. This research holds significant theoretical and practical implications for identifying optimal shale gas development zones and improving gas recovery efficiency.

1 Introduction

Shale gas refers to natural gas that is formed through the hydrocarbon generation and evolution in shale, and is predominantly stored in the micro-fractures and pores of shale reservoirs in a free or adsorbed state, without significant migration (Jia et al., 2012; Zhang et al., 2008; Yang et al., 2025). After more than a decade of exploration and development, China has become a major global shale gas producer, with its share in the energy mix continuously increasing. In the future, it is poised to become one of the most reliable alternative energy sources (Zou et al., 2020; Zhao et al., 2020). The factors influencing production levels include geological and engineering factors, but geological factors such as the level of gas content or the degree of enrichment of shale gas are fundamental. Some scholars have recognized that shale gas can be redistributed to higher locations through relay transmission between pores and micro-fractures, resulting in high production from shale gas wells in higher locations, such as the Zhaotong, Ningxi, and Yuxi regions (Liang et al., 2021; Chen et al., 2021; Fu et al., 2021; Gao et al., 2024). Therefore, the coupled pore-fracture system has attracted the attention of academia and industry.

Shale gas reservoirs have extremely low porosity and permeability, with a large number of micro- and nano-pores developed (Curtis et al., 2012; Gou et al., 2023). This has attracted numerous scholars to analyze the microscopic pore development characteristics, pore size distribution, and pore connectivity of shale gas through experimental methods such as field emission scanning electron microscope (FE-SEM), focused ion beam-scanning electron microscope (FIB-SEM), low-pressure gas adsorption, nuclear magnetic resonance (NMR), asequence and sedimentary characteristics of Lower Silurind mercury injection capillary pressure (MICP) (Klaver et al., 2015; Clarkson et al., 2013; Guo Q. Y. et al., 2019; Hu et al., 2012; Fu et al., 2019; Gao and Hu, 2016). Therefore, numerous scholars have analyzed the impact of organic matter type, total organic carbon (TOC) content, mineral composition, diagenesis, and hydrocarbon generation evolution on the pore development of shale reservoirs using the above methods (Ma et al., 2020; Yang et al., 2020; Cui et al., 2024; Huang et al., 2023; Wang et al., 2025). It can be seen that the research on pore characteristics and development control factors of shale reservoir has become quite mature.

Foreign scholars have conducted early research on fractures. Hill and Nelson (2000) believed that natural fracture systems are the main storage spaces for free gas in shale reservoirs, and that organic matter abundance is an important factor affecting fracture development. Curtis (2002) believed that for shale reservoirs with low matrix permeability, the degree of fracture development is an important controlling factor for gas production capacity. Domestic scholars have concluded through fracture research that different types of micro fractures can be formed in the process of hydrocarbon generation and diagenetic evolution, which provides a large amount of reservoir space for shale gas (Cui et al., 2024) and is conducive to shale reservoir fracturing and shale gas production (Zhou et al., 2020; Huang et al., 2023). However, in complex structural areas, due to excessive tectonic activity, excessive fracture development can lead to the dissipation of shale gas, resulting in smaller or disappearing pores, which has a detrimental effect on shale gas preservation (Liu et al., 2016; Shi and Chen, 2023). Overall, current fracture research mainly focuses on fracture genesis, classification, activity periods, development control factors, fracture distribution prediction (Wang et al., 2016; Wang et al., 2015; Yang et al., 2022; He et al., 2017).

In order to deeply analyze the influence of pore structure development on shale gas enrichment and exploitation, scholars pay more and more attention to the micro connectivity of shale reservoir (Sun et al., 2023; Yang et al., 2022). For example, the connectivity of shale reservoir pores is characterized by imbibition media with different wettability (Hu et al., 2012; Fu et al., 2019; Guo Q. Y. et al., 2019), and the characteristics and permeability of shale reservoir pores and fractures are jointly characterized using multiple methods to reveal the connectivity of shale reservoirs (Gao et al., 2024; Sun et al., 2023; Yang et al., 2022). Although these methods all focus on the connectivity of shale reservoirs, current research on shale pores and fractures remains fragmented, with studies typically examining these components in isolation rather than as an integrated system (Chen et al., 2021; Fu et al., 2019; Gao and Hu, 2016). Therefore, it is necessary to establish new parameters for pore-fracture connectivity based on the degree of micro-pore and micro-fracture development, classify pore-fracture combination types, reveal the dynamic process of shale gas migration and enrichment, and guide the selection of shale gas areas, layers, and scientific exploitation.

This study selects downhole cores from the Longmaxi Formation in the Sichuan Basin and innovatively conducts a joint experimental test involving MICP, wood’s metal (WM) impregnation, and FE-SEM. This achieves visualization of the high-pressure mercury injection process and quantitatively identifies the micro-pore and fracture volumes in shale reservoirs. Based on the characteristics of the MICP curve, the combination types of micro-pores and fractures in shale reservoirs are classified, and the impact of different pore and fracture combination types on exploration and development is explored.

2 Geological background

According to the characteristics of regional tectonic deformation and oil and gas distribution in the basin, the Sichuan Basin is divided into six secondary structural units: the North Sichuan low-flat belt, the Middle Sichuan low-flat belt, the Eastern Sichuan high-steep belt, the Western Sichuan low-steep belt, the Southwestern Sichuan low-steep belt, and the South Sichuan low-steep belt. The study area is located in Luzhou area in the South Sichuan low-steep belt (Figure 1a). The current burial depth of Longmaxi formation in the study area is more than 3,500 m (Figure 1b), which is a typical deep shale gas reservoir. In this study, the downhole cores of five typical shale gas wells in the study area were selected for research. After nearly two decades of exploration and development, it has been verified that the organic-rich siliceous shale within the high GR (natural gamma ray log value) interval at the base of the Longmaxi Formation possesses commercial exploitation value, serving as the target interval of this study (Figure 1c).

FIGURE 1

3 Sampling and methodologies

3.1 Sample information

In this study, the cores of deep shale gas reservoirs in 5 shale gas wells in Luzhou area of Sichuan Basin were selected for experimental test and analysis. These five wells are respectively L211, L217, Y101H3-8, Y101H91-4, and Y101H65-5, whose locations are displayed in Figure 1b. In order to enhance the comparability between samples, the samples in this study were taken from the high GR section at the bottom of Longmaxi formation in five wells. Among them, two core samples were selected from each well for experimental testing. The mineral composition, TOC content, and porosity of these samples are shown in Table 1. The samples have a TOC content of over 4%, with quartz content exceeding 50% and clay content less than 20%, indicating good gas generation potential and fracturing conditions. However, there are inter-well differences in porosity values, ranging from 3.52% to 5.81%, with an average of 4.88%. Among them, the porosity of wells L211 and Y101H91-4 is relatively low, falling below the average level.

3.2 Experimental testing analysis

3.2.1 Mercury injection capillary pressure (MICP) measurement

The Micromeritics AutoPore IV 9520 porosimeter was used for MICP testing, with a maximum pressure of 413.6 MPa. According to the Washburn equation (P = 2σcosθ/r, where P is the pressure of mercury intrusion; θ is the contact angle; σ is the surface tension; r is the radius of pore throat), the minimum detectable pore throat radius is approximately 3 nm under maximum pressure. Before conducting the MICP test, the shale samples were cut into cubes of approximately 1 cm³, polished on each surface to minimize surface roughness effects, and then dried at 60 °C. The bulk density was obtained by mercury encapsulation under initial pressure (0.034 MPa). To distinguish between residual mercury and actual pore intrusion, repeated MICP measurements were conducted on the same sample.

3.2.2 Wood’s metal (WM) impregnation

3.2.3 FE-SEM

4 Results

4.1 MICP test results

As the mercury injection pressure increases from 0.035 MPa to 413.6 MPa, the cumulative mercury injection volume of 10 shale samples gradually increases (Figure 2a). For some samples (such as L217-1/Y101H91-4-1/Y101H65-5-1/Y101H65-5-2, etc.), when the pressure reaches its maximum value (413.6 MPa), it begins to decrease, but the cumulative mercury injection volume continues to increase before gradually decreasing. This “over topping” phenomenon mainly occurs in low-porosity and low-permeability shale samples, which is due to the recovery of mercury volume expansion during the pressure reduction process (Sun et al., 2020). The MICP porosity values of the samples, calculated based on the cumulative mercury intrusion volume under maximum pressure and the bulk density measured under initial pressure, are presented in Table 2. The samples show significant variation in MICP porosity values, ranging from 0.42% to 3.16%, with an average of 1.56%. Well L217 exhibits the highest average MICP porosity (2.68%), while Well L211 shows the lowest (0.56%). According to the Washburn equation, the connected pore volume related to the corresponding pore throat size was calculated. The pore throat size distribution curve shows that more than half of the total pore volume is connected by pore throats smaller than 20 nm, and the pores and micro-fractures larger than 150 nm to several micrometers are also developed (Figure 2). In this study, the porosity of overmature marine shale samples with pore throats ranging from 3 to 20 nm accounts for 50.34%∼82.56% of the MICP porosity, providing substantial storage space for shale gas.

FIGURE 2

4.2 FE-SEM photos of WM immersion

To gain a clear understanding of the phenomenon of mercury molecules injecting into pores and micro-fractures during MICP process, we designed WM immersion experiments under different pressure conditions, and observed the locations of WM filling into pores and micro-fractures under FE-SEM. Sample L217-1 was subjected to WM impregnation at pressures of 500, 1700, 8,000, and 30,000 psi, respectively.

At a low pressure of 500 Psi, FE-SEM observations show that WM mainly occupied micro-fractures with large width and long length parallel to the lamination (Figure 3a). However, micro-fractures with both narrow width and short length remain unfilled by WM, particularly those diagenetic contraction fractures that exhibit almost no WM infiltration (Figure 3b). As the pressure increases to 1700 Psi, corresponding to a pore size of 105 nm, WM filled micro-fractures and a small number of pores, but organic pores and smaller pore size inorganic pores are basically not filled (Figure 3c). As the pressure increases to 8,000 Psi, corresponding to a pore size of 25 nm, WM successfully filled small organic pores, small-diameter inorganic pores, and diagenesis-related micro-fractures in sample L217-1 (Figures 3d,e). As the pressure increases to 30,000 psi, the pore size accessible to WM decreases to 7 nm. Under a low magnification, it is clearly visible that the pores and micro-fractures are filled with a large amount of WM (Figure 3f). Upon continuous magnification and observation of micro-fractures and pores, it was found that some interconnected micro-fractures with both narrow width and short length exhibited a partial filling with WM, which may be related to the volume shrinkage of the WM after cooling (Figure 3g). The pores (organic and inorganic pores) connected to the micro-fractures were fully filled with WM (Figures 3g,h). However, due to poor connectivity, a large number of intragranular pores did not exhibit WM filling (Figure 3i). It can be seen that with the increase of pressure, WM is regularly immersed in the pores and micro-fractures of shale. WM typically penetrates into larger-scale micro-fractures, then fills smaller-scale micro-fractures and inorganic pores, subsequently entering inorganic and organic pores with smaller diameters until most of the connected pores and micro-fractures are occupied by WM.

FIGURE 3

5 Discussion

5.1 Stage division of liquid mercury into pores and micro-fractures

The pore-fracture network in shale consists of a fracture system connected to the matrix pores (Davudov et al., 2020). Regarding the study of pore connectivity, scholars have employed methods such as spontaneous imbibition, CM-SANS, and FIB-SEM 3D reconstruction to analyze the pore connectivity within shale reservoirs, and these studies indicate that the organic pore connectivity in high-overmature marine shale gas reservoirs is relatively good (Sun et al., 2020; Sun et al., 2023; Hu et al., 2015; Hu et al., 2017), but the connectivity between micro-fractures and pores remains poorly understood. To address this knowledge gap, this study quantitatively analyses the types of pore-fracture combinations through FE-SEM observations after WM immersion experiments and MICP.

A large number of nano pores and micro nano scale micro-fractures are developed in shale gas reservoir. Numerous studies have shown that in the overmature Longmaxi Formation shale reservoirs, the pore diameters of inorganic pores are larger than those of organic pores, but both are smaller than those of micro-fractures (Yu et al., 2019; Yang et al., 2019). Among the inorganic pores, there are numerous isolated intragranular pores, and the pore diameters and quantities of organic pores in kerogen are significantly smaller than those in pyrobitumen (Teng et al., 2022; Liu et al., 2022; Guo X. W. et al., 2019). The MICP and WM immersion experiments involve injecting liquid metal fluid into the shale pore-fracture system under external force. This process follows the Washburn equation, which indicates that the mercury first enters fractures and pores with larger pore sizes, and then moves into pores with relatively smaller pore sizes. Therefore, by taking the partial derivative of the cumulative mercury injection volume with respect to the injection pressure, we can effectively distinguish the rate of liquid mercury filling into the pore-fracture system, which reflects the connectivity of pore-fracture system in shale samples. Figure 4 shows that the mercury injection rate (dV/dP) of the shale reservoir in the study area gradually decreases with pressure increases. The pressure corresponding to the first minimum value after the highest peak of dV/dP is used to distinguish micro-fractures, and the pressure corresponding to the minimum value of dV/dP tending to 0 for the first time is used to distinguish pores. Taking L217 sample as an example, according to the value of dV/dP, mercury injection process can be roughly divided into three intervals: high mercury injection rate (<1.2 MPa), moderate mercury injection rate (1.2∼16 MPa), and low mercury injection rate (>16 MPa) (Figure 4). It is worth noting that due to differences in the pore structure among shale samples, the pressure corresponding to the mercury injection rate range at different stages also varies. For example, the two samples from well Y101H65-5 have the lowest mercury injection rate (Figure 4) and also can be divided into three intervals: high mercury injection rate (<3 MPa), moderate mercury injection rate (3∼16 MPa), and low mercury injection rate (>16 MPa), indicating poor connectivity of pore-fracture system in the shale reservoir of this well. On the contrary, the dV/dP of shale samples from Well L217 and Well Y101H3-8 is relatively high, indicating good connectivity of pore-fracture system in the shale reservoir.

FIGURE 4

To better distinguish the stages of mercury molecules entering pores and micro-fractures during the MICP experimental process, the L217-1 sample that had completed the WM immersion experiment was selected for detailed analysis. Based on the change in dV/dP value, the mercury injection stage of the MICP is divided into three stages (Figure 5). The dV/dP value corresponding to point A is the first minimum value after the highest peak, indicating a turning point where the mercury injection rate changes from low to high and then back to low, and the pressure corresponding to point A is 1.2 MPa. Point B represents the minimum value of dV/dP tending to 0 for the first time, revealing the process where the mercury injection rate changes from high to low and then slight increase. The pressure value corresponding to point B is 16 MPa. So, when the mercury injection pressure is less than 1.2 MPa, the main space where mercury molecules penetrate is micro-fractures, with little penetration into pores and smaller-scale diagenetic contraction fractures (Figures 3a,b). When the injection pressure ranges from 1.2 to 16 MPa, micro-fractures are completely filled with WM, and then some inorganic pores begin to be filled with WM (Figure 3c). As pressure gradually increases, WM fills a large number of intergranular pores and organic pores, but intragranular pores and organic pores with poor connectivity are not filled (Figures 3g–i). Therefore, the mercury injection curve is categorized into three stages: the micro-fracture filling stage with mercury (<1.2 MPa), the stage of mercury filling pores and micro-fractures (1.2∼16 MPa), and the pore filling stage with mercury (>16 MPa). This sequential penetration pattern, where mercury infiltrates micro-fractures first followed by pores, is consistent with previous studies (Sun et al., 2023; Yu et al., 2019). The same analytical approach can be applied to other samples for pore-scale classification during mercury injection processes.

FIGURE 5

5.2 Quantitative characterization of pores-fractures combination type

Some scholars have classified the types of micro pore-fracture combinations in shale gas reservoirs, mainly through qualitative understanding of the contact relationships between micro-fractures and pores based on electron microscopy observations (Fu et al., 2021), but quantitative characterization has not been achieved. In the previous section, the mercury injection curve was divided into three stages, revealing the process of mercury entering pores and fractures, and achieving quantitative division of shale reservoir pores and fractures using the mercury injection curve, which lays the foundation for quantitative analysis of pore-fracture combination type.

The pore-fracture combination types in shale reservoirs aim to understand the connectivity between pores and micro-fractures, and further enhance the understanding of the difficulty of shale gas enrichment and extraction. MICP can characterize pores ranging from 3 nm to several micrometers, providing significant advantages in understanding pore-fracture structure (Yu et al., 2019). MICP involves injecting fluid into the pores and micro-fractures of shale reservoirs under external force. When characterizing pore connectivity using MICP, the main focus is on the mercury injection volume and rate. Based on this, we have defined the micro pore-fracture combination parameter (Q_P-F), which can be calculated using the following formula:

Where, Q_P is pore connectivity parameter; Q_F is fracture connectivity parameter; C is total connectivity parameter. C, Q_P and Q_F can express as Equations 2–4, respectively.

Where, Φ_T is total porosity of shale core, %; Φ_C is connective porosity calculated by MICP, %; V_P is mercury injection volume after pressure exceeds point A, mL/g; V_T is accumulated mercury injection volume, mL/g; V_F is mercury injection volume at pressure less than point A, mL/g; K₁ is average mercury injection rate after pressure exceeds point A, (mL/g)/MPa; K₂ is average mercury injection rate at pressure less than point A, (mL/g)/MPa.

From Equation 1, it can be understood that Q_P-F represents the pore volume connected by the unit micro-fracture volume. The larger the value of Q_P-F, the better the connectivity of the pore-fracture system. Relevant calculation parameters were listed in Table 2. Our calculations reveal that when the Q_P value exceeds the Q_F value by an order of magnitude, the Q_P-F value is greater than 1; conversely, the Q_P-F value is less than 1. From the geological perspective of shale gas, it is believed that shale reservoirs are characterized by a large number of organic pores and a relatively small number of micro-fractures (Li et al., 2025; Chen et al., 2019). However, these micro-fractures serve as connecting pores, enhancing the migration of shale gas (Li et al., 2024; Liu et al., 2024). Studies also indicate that in areas with fault development, shale reservoirs with a high degree of natural fractures and micro-fractures exhibit very poor pore development, often resulting in poor reservoir performance and gas-bearing properties (Hu et al., 2024; Fu et al., 2021; Fan et al., 2020). It is evident that a higher degree of micro-fractures development hinders pore development. This serves as indirect evidence that our calculation results align with the geological features of shale gas reservoirs. Therefore, greater attention should be directed towards the connectivity between micro-fractures and pores, particularly in scenarios of limited micro-fractures development.

The Q_P-F values of all shale samples are distributed between 0.01 and 9.67, with an average value of 1.82, indicating that the degree of pore development in the shale reservoir is higher than that of micro-fractures. Among them, the Q_P-F values of the four samples from Well L217 and Well Y101H3-8 are greater than 1, with an average value of 4.02, indicating that the shale reservoirs of these two wells exhibit good pore-fracture connectivity. The Q_P-F value of Well L211 is the smallest, with two samples being 0.20 and 0.10 respectively, indicating that the shale reservoirs of this well exhibit poor pore-fracture connectivity. Based on the magnitude of the Q_P-F value, the pore-fracture combination types in shale reservoirs are categorized into three types: strongly connected pore-fracture type (SCPFT) (Q_P-F>1), moderately connected pore-fracture type (MCPFT) (0.5<Q_P-F<1), and weakly connected pore-fracture type (WCPFT) (Q_P-F<0.5). The Q_P value of SCPFT shale sample is much greater than Q_F value (Table 2), indicating that the degree of pore development is much better than that of micro-fractures. The Q_P value of WCPFT shale sample is equal to or slightly greater than the Q_F value (Table 2), reflecting that the pores of shale reservoir are developed, but the connectivity between pores and micro-fractures is poor. MCPFT shale reservoir is between the first two types. It is worth noting that the Q_P-F value of the shale reservoir of well L217 is much higher than that of other shale gas wells, especially the L217-2 sample, because the well is located far away from other wells and is relatively flat (Figure 1b), with weak structural deformation, high reservoir overpressure and high degree of pore development. Therefore, Q_P-F can indirectly distinguish the quality of shale reservoir and guide the optimization of favorable areas for shale gas exploration.

5.3 Influence of pores-fractures combination type on exploration and development

Micro pores and micro-fractures in shale reservoirs are important reservoir spaces and seepage channels of shale gas, which affect the enrichment and single well recovery of shale gas. Therefore, the pore-fracture combination type has an important impact on the exploration and development of shale gas. In the previous section, according to the Q_P-F value, the pore-fracture combination types of shale reservoir are divided into three categories. The following mainly analyzes the reservoir capacity and gas bearing property of shale reservoirs with different pore-fracture combination types.

Figure 6a shows that the porosity of all SCPFT shale samples is high, with an average of 5.28%; but the porosity of MCPFT and WCPFT shale samples varies greatly. It shows that SCPFT shale reservoir has strong reservoir capacity, while MCPFT and WCPFT shale reservoirs have strong heterogeneity of reservoir capacity. The corresponding shale gas reservoir has a similar variation law that the gas content of all SCPFT shale samples is greater than 7.5 m³/t, with an average of 8.28 m³/t, and the gas content of MCPFT and WCPFT shale samples is great difference (Figure 6b). Shale gas reservoirs exhibit typical characteristics of being both source and reservoir, and shale gas does not need to be migrated or migrated for a short distance in situ (Curtis et al., 2012; Liu et al., 2016; Hu et al., 2024). This will cause some disconnected pores to still contain gas in MCPFT and WCPFT shale reservoirs. However, it is difficult to exploit shale gas in these pores during the exploitation of shale gas. SCPFT shale gas reservoirs have high porosity (5.28%), high gas content (8.28 m³/t), and good pore-fracture connectivity which is conducive to shale gas exploitation. At the same time, this kind of reservoir has high quartz content, which is conducive to fracturing operation and forms complex fracture network structure. Therefore, SCPFT shale gas reservoir is the most important target in the process of shale gas exploration and development.

FIGURE 6

In order to further analyze the differences of shale gas enrichment and percolation in different pore-fracture combination types, the correlation between pore-fracture combination type parameters and gas bearing property and permeability was drawn in this study. The results show that there is a good logarithmic relationship between Q_P and gas content (Figure 7a), there is no correlation between Q_F and gas content (Figure 7b), and there is a positive correlation between Q_P-F and gas content to some extent (Figure 7c), suggesting that shale gas in shale reservoirs is predominantly stored within pores, with fractures exhibiting a relatively weaker storage capacity. A shale reservoir characterized by more developed pores exhibits superior productivity (Fu et al., 2021; Liu et al., 2024), aligning with our established understanding. The correlation between Q_P-F and gas content is not as good as Q_P (Figures 7a,c). Q_P-F merely reflects the degree of pore development to a certain extent, but more crucially, it embodies the connectivity of pore-fracture system. So there is a good positive correlation between Q_P-F and shale permeability (Figure 7d). Generally, the higher the permeability of a shale reservoir, the greater its flow capacity becomes, resulting in better productivity for a single well (Davudov et al., 2020). Figure 8 demonstrates a positive correlation between the Q_P-F value and the estimated ultimate recovery (EUR) of a single well, indicating that shale gas wells with superior pore-fracture connectivity tend to exhibit higher gas production rates.

FIGURE 7

FIGURE 8

Among the five shale gas wells selected for this study, Wells L217 and Y101H3-8 exhibit Q_P-F values for their shale reservoirs exceeding 1, classifying them as SCPFT. This kind of shale reservoir has good development of pores and micro-fractures, and the contact degree between pores and micro-fractures is high (Figure 9a), which is conducive to the enrichment and exploitation of shale gas. Well L211, on the other hand, has a Q_P-F value for its shale reservoir below 0.5, categorizing it as WCPFT. The development degree of pores and micro-fractures in this kind of shale reservoir is poor, and the connectivity of pores-fractures is poor (Figure 9c), which is not conducive to the enrichment and exploitation of shale gas. For Wells Y101H91-4 and Y101H65-5, their shale reservoir Q_P-F values are within the range of 0.16–0.97 but above 1 (note: this should be a clarification or correction as per context, possibly meaning “above 0.16 but some values are less than 1″), with average values of 0.4 and 0.64, respectively. Consequently, Well Y101H91-4 falls into the WCPFT category, while Well Y101H65-5 is classified as MCPFT, which shows good development characteristics of pores and micro-fractures, but the connectivity is relatively weak (Figure 9b). This elucidates the reason why Wells Y101H91-4 and Y101H65-5 exhibit high porosity and gas content; however, despite undergoing the same fracturing process, they failed to yield satisfactory gas production results. Therefore, in shale reservoirs where Q_P-F is less than 1, the pores-fracture connectivity is inadequate. To address this, the fracturing operation method must be optimized, including measures such as increasing the number of perforation clusters, enhancing the intensity of fracturing fluid, and augmenting the sand concentration. These adjustments aim to achieve a more intricate fracture network structure, thereby unlocking more pores with limited connectivity and ultimately enhancing shale gas recovery.

FIGURE 9

6 Conclusion

References

• 1

  ChenG. S.WuJ. F.LiuY.HuangH. Y.ZhaoS. X.ChangC.et al (2021). Geology-engineering integration key technologies for ten billion cubic meters of shale gas productivity construction in the Southern Sichuan Basin. Nat. Gas. Ind.41 (1), 72–82. 10.3787/j.issn.1000-0976.2021.01.006

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  ChenZ. Y.SongY.JiangZ. X.LiuS. B.LiZ.ShiD. S.et al (2019). Identification of organic matter components and organic pore characteristics of marine shale: a case study of wufeng-longmaxi shale in Southern Sichuan basin, China. Mar. Petroleum Geol.109, 56–69. 10.1016/j.marpetgeo.2019.06.002

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  ClarksonC. R.SolanoN.BustinR.BustinA.ChalmersG.HeL.et al (2013). Pore structure characterization of North American shale gas reservoirs using USANS/SANS, gas adsorption, and mercury intrusion. Fuel103, 606–616. 10.1016/j.fuel.2012.06.119

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  CuiZ.SunM. D.MohammadianE.HuQ. H.LiuB.OstadhassanM.et al (2024). Characterizing microstructural evolutions in low-mature lacustrine shale: a comparative experimental study of conventional heat, microwave, and water-saturated microwave stimulations. Energy294, 130797. 10.1016/j.energy.2024.130797

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  CurtisJ. B. (2002). Fractured shale-gas systems. AAPG Bull.86 (11), 1921–1938. 10.1306/61EEDDBE-173E-11D7-8645000102C1865D

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  CurtisM.CardottB.SondergeldC.RaiC. S. (2012). Development of organic porosity in the woodford shale with increasing thermal maturity. Int. Jounal Coal Geol.103 (23), 26–31. 10.1016/j.coal.2012.08.004

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  DavudovD.MoghanlooR. G.ZhangY. X. (2020). Interplay between pore connectivity and permeability in shale sample. Int. J. Coal Geol.220, 103427. 10.1016/j.coal.2020.103427

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  FanC. H.LiH.QinQ. R.HeS.ZhongC. (2020). Geological conditions and exploration potential of shale gas reservoir in wufeng and longmaxi formation of southeastern sichuan basin, China. J. Petroleum Sci. Eng.191, 107138–4105. 10.1016/j.petrol.2020.107138

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  FuY. H.JiangY. Q.DongD. Z.HuQ. H.LeiZ. A.PengH.et al (2021). Microscopic pore-fracture configuration and gas-filled mechanism of shale reservoirs in the Western Chongqing area, Sichuan Basin, China. Petroleum Explor. Dev.48 (05), 1063–1076. 10.1016/s1876-3804(21)60091-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  FuY. H.JiangY. Q.WangZ. L.HuQ. H.XieJ.NiG. S.et al (2019). Non-connected pores of the longmaxi shale in southern Sichuan Basin of China. Mar. Petroleum Geol.110, 420–433. 10.1016/j.marpetgeo.2019.07.014

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  GaoP.XiaoX. M.HuD. F.LashG. G.LiuR. B.ZhangB. Y.et al (2024). Comparison of silica diagenesis between the lower Cambrian and lower Silurian shale reservoirs in the middle-upper yangtze platform (southern China). AAPG Bull.108 (6), 971–1003. 10.1306/01242422096

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  GaoZ. Y.HuQ. H. (2016). Initial water saturation and imbibition fluid affect spontaneous imbibition into barnett shale samples. J. Nat. Gas Sci. Eng.34, 541–551. 10.1016/j.jngse.2016.07.038

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  GouQ.XuS.HaoF.ZhangL.ChenZ.BaiN. (2023). Petrography and mineralogy control the nm-μm-scale pore structure of saline lacustrine carbonate-rich shales from the Jianghan Basin, China. Mar. Petroleum Geol.155, 106399. 10.1016/j.marpetgeo.2023.106399

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  GuoQ. Y.XuS.HaoF.YangF.ZhangB. Q.ShuZ. G.et al (2019). Full-scale pores and micro-fractures characterization using FE-SEM, gas adsorption, nano-CT and micro-CT: a case study of the Silurian longmaxi formation shale in the fuling area, Sichuan Basin, China. Fuel253, 167–179. 10.1016/j.fuel.2019.04.116

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  GuoX. W.QinZ. J.YangR.DongT.HeS.HaoF.et al (2019). Comparison of pore systems of clay-rich and silica-rich gas shales in the lower Silurian longmaxi formation from the jiaoshiba area in the eastern Sichuan Basin, China. Mar. Petroleum Geol.101, 265–280. 10.1016/j.marpetgeo.2018.11.038

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  HeZ. L.NieH. K.ZhaoJ. H.LiuW. X.FangB.ZhangW. T. (2017). Types and origin of nanoscale pores and fractures in wufeng and longmaxi shale in Sichuan Basin and its periphery. J. Nanosci. Nanotechnol.17 (9), 6626–6633. 10.1166/jnn.2017.14425

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  HillD. G.NelsonC. R. (2000). Gas productive fractured shales: an overview and update. EurekaMag6 (2), 4–13. 10.2523/29580-ms

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  HuQ. H.EwingR. P.DultzS. (2012). Low pore connectivity in natural rock. J. Contam. Hydrology133, 76–83. 10.1016/j.jconhyd.2012.03.006

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  HuQ. H.EwingR. P.RoweH. D. (2015). Low nanopore connectivity limits gas production in barnett formation. J. Geophys. Res. Solid Earth120, 8073–8087. 10.1002/2015JB012103

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  HuQ. H.LiuH.YangR.ZhangY. X.KibriaG.SahiS.et al (2017). Applying molecular and nanoparticle tracers to study wettability and connectivity of longmaxi formation in southern China. J. Nanosci. Nanotechnol.17 (9), 6284–6295. 10.1166/jnn.2017.14428

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  HuZ. Q.LiuZ. Q.DuW.BaiZ.WangR. Y.SunC. X.et al (2024). Key geological factors governing sweet spots in the wufeng-longmaxi shales of the Sichuan basin, China. Energy Geosci.5 (3), 100268. 10.1016/j.engeos.2023.100268

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  HuangX.KangZ.ZhaoJ.WangG.ZhangH.YangD. (2023). Experimental investigation on micro-fracture evolution and fracture permeability of oil shale heated by water vapor. Energy277, 127677. 10.1016/j.energy.2023.127677

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  JiaC. Z.ZhengM.ZhangY. F. (2012). Unconventional hydrocarbon resources in China and the prospect of exploration and development. Petroleum Explor. Dev.39 (2), 139–146. 10.1016/S1876-3804(12)60026-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  KlaverJ.DesboisG.LittkeR.UraiJ. (2015). BIB-SEM characterization of pore space morphology and distribution in postmature to overmature samples from the haynesville and bossier shales. Mar. Petroleum Geol.59, 451–466. 10.1016/j.marpetgeo.2014.09.020

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  LiH.HeS.RadwanA. E.XieJ. T.QinQ. R. (2024). Quantitative analysis of pore complexity in lacustrineorganic-rich shale and comparison to marine shale: insights from experimental tests and fractal theory. Energy and Fuel38 (17), 16171–16188. 10.1021/acs.energyfuels.4c03095

  • CrossRef
  • Google Scholar

• 26

  LiJ.ZhangQ.JiangW.LiH.LiL.XueT. F.et al (2025). Lithological controls on pore structure and their implications for deep shale gas reservoir quality in the longmaxi formation, Luzhou area, Southern Sichuan Basin, China. Energy and Fuel39, 1541–1558. 10.1021/acs.energyfuels.4c05247

  • CrossRef
  • Google Scholar

• 27

  LiangX.ShanC. A.ZhangZ.XuZ. Y.XuJ. B.WangW. X.et al (2021). Three-dimensional closed system accumulation model of taiyang anticline mountain shallow shale gas in the zhaotong demonstration area. Acta Geol. Sin.95 (11), 3380–3399. 10.3969/j.issn.0001-5717.2021.11.015

  • CrossRef
  • Google Scholar

• 28

  LiuB.MastalerzM.SchieberJ. (2022). SEM petrography of dispersed organic matter in black shales: a review. Earth-Science Rev.224, 103874. 10.1016/j.earscirev.2021.103874

  • CrossRef
  • Google Scholar

• 29

  LiuJ. F.HeX.XueF. J.DaiJ. J.YangY. X.HuangH. Y.et al (2024). The influence of natural fractures of multi-feature combination on seepage behavior in shale reservoirs. J. Min. Strata Control Eng.6 (1), 013437. 10.13532/j.jmsce.cn10-1638/td.20240018.001

  • CrossRef
  • Google Scholar

• 30

  LiuS. G.DengB.ZhongY.RanB.YongZ. Q.SunW.et al (2016). Unique geological features of burial and superimposition of the lower Paleozoic shale gas across the Sichuan Basin and its periphery. Earth Sci. Front.23 (1), 11–28. 10.13745/j.esf.2016.01.002

  • CrossRef
  • Google Scholar

• 31

  MaX. H.XieJ.YongY.ZhuY. Q. (2020). Geological characteristics and high production control factors of shale gas reservoirs in Silurian longmaxi formation, southern Sichuan Basin, SW China. Petroleum Explor. Dev.47 (5), 901–915. 10.1016/s1876-3804(20)60105-7

  • CrossRef
  • Google Scholar

• 32

  ShiQ.ChenP. (2023). “Rigid-elastic chimera” pore model and overpressure porosity measurement method for shale: a case study of the deep overpressure siliceous shale of Silurian longmaxi formation in southern Sichuan Basin, SW China. Petroleum Explor. Dev.50 (1), 113–124. 10.1016/S1876-3804(22)60374-4

  • CrossRef
  • Google Scholar

• 33

  SunM. D.DuanX. G.LiuQ.BlachT.OstadhassanM.LiuB.et al (2023). The importance of pore-fracture connectivity in overmature marine shale for methane occurrence and transportation. Mar. Petroleum Geol.157, 106495. 10.1016/j.marpetgeo.2023.106495

  • CrossRef
  • Google Scholar

• 34

  SunM. D.ZhangL. H.HuQ. H.PanZ. J.YuB. S.SunL. W.et al (2020). Multiscale connectivity characterization of marine shales in southern China by fluid intrusion, small-angle neutron scattering (SANS), and FIB-SEM. Mar. Petroleum Geol.112, 104101. 10.1016/j.marpetgeo.2019.104101

  • CrossRef
  • Google Scholar

• 35

  TengJ.LiuB.MastalerzM.SchieberJ. (2022). Origin of organic matter and organic pores in the overmature ordovician-silurian wufeng-longmaxi shale of the Sichuan Basin, China. Int. J. Coal Geol.253, 103970. 10.1016/j.coal.2022.103970

  • CrossRef
  • Google Scholar

• 36

  WangL. C.ZhangW.CaoZ. Z.XueY.XiongF. (2025). Coupled effects of the anisotropic permeability and adsorption-induced deformation on the hydrogen and carbon reservoir extraction dynamics. Phys. Fluids37, 066608. 10.1063/5.0270765

  • CrossRef
  • Google Scholar

• 37

  WangR. Y.DingW. L.GongD. J.ZengW. T.WangX. H.ZhouX. H.et al (2016). Development characteristics and major controlling factors of shale fractures in the Lower Cambrian niutitang formation, southeastern Chongqing-northern Guizhou area. Acta Pet. Sin.37 (07), 832–845. 10.7623/syxb201607002

  • CrossRef
  • Google Scholar

• 38

  WangY. M.HuangJ. L.LiX. J.DongD. Z.WangS. F.GuanQ. Z. (2015). Quantitative characterization of fractures and pores in shale beds of the Lower Silurian, longmaxi formation, Sichuan Basin. Nat. Gas. Ind. B35 (09), 481–488. 10.1016/j.ngib.2015.12.002

  • CrossRef
  • Google Scholar

• 39

  YangG. Q.ZengJ. H.QiaoJ. C.LiuY. Z.LiuS. N.JiaW. T.et al (2022). Role of organic matter, mineral, and rock fabric in the full-scale pore and fracture network development in the mixed lacustrine source rock system. Energy and Fuels36 (15), 8161–8179. 10.1021/acs.energyfuels.2c01484

  • CrossRef
  • Google Scholar

• 40

  YangR.HuQ. H.HeS.HaoF.GuoX. S.YiJ. Z.et al (2019). Wettability and connectivity of overmature shales in the fuling gas field, Sichuan Basin (China). AAPG Bull.103 (3), 653–689. 10.1306/09051817138

  • CrossRef
  • Google Scholar

• 41

  YangW.CaiJ. F.WangQ. Y.CuiZ. J.CuiZ.XuL.et al (2020). The controlling effect of organic matter coupling with organic matter porosity on shale gas enrichment of the wufeng-longmaxi marine shale. Petroleum Sci. Bull.5 (2), 148–160. 10.3969/j.issn.2096-1693.2020.02.015

  • CrossRef
  • Google Scholar

• 42

  YangW.YangR.WangM.BaoS. J.YangL. L.HouH. D.et al (2025). Impact of lithofacies categories on inflection-point behaviors in micro-nano pore-structure evolution: implications for differential reservoir-forming mechanisms and “sweet-spot intervals” identification for marine gas shales. Mar. Petroleum Geol.176, 107336. 10.1016/j.marpetgeo.2025.107336

  • CrossRef
  • Google Scholar

• 43

  YuY. X.LuoX. R.WangZ. X.ChengM.LeiY. H.ZhangL. K.et al (2019). A new correction method for mercury injection capillary pressure (MICP) to characterize the pore structure of shale. J. Nat. Gas Sci. Eng.68, 102896. 10.1016/j.jngse.2019.05.009

  • CrossRef
  • Google Scholar

• 44

  ZhangJ. C.XuB.NieH. K.WangZ. Y.LinT.JiangS. L.et al (2008). Exploration potential of shale gas resources in China. Nat. Gas. Ind.28 (6), 136–140. 10.3787/j.issn.1000-0976.2008.06.040

  • CrossRef
  • Google Scholar

• 45

  ZhaoW. Z.JiaA. L.WeiY. S.WangJ. L.ZhuH. Q. (2020). Progress in shale gas exploration in China and prospects for future development. China Pet. Explor.25 (1), 31–44. 10.3969/j.issn.1672-7703.2020.01.004

  • CrossRef
  • Google Scholar

• 46

  ZhouX. J.YongR.FanY.ZengB.SongY.GuoX. W.et al (2020). Influence of natural fractures on fracturing of horizontal shale gas wells and process adjustment. China Pet. Explor.25 (06), 94–104. 10.3969/j.issn.1672-7703.2020.06.010

  • CrossRef
  • Google Scholar

• 47

  ZouC. N.PanS. Q.JingZ. H.GaoJ. L.YangZ.WuS. T.et al (2020). Shale oil and gas revolution and its impact. Acta Pet. Sin.41 (1), 1–12. 10.7623/syxb202001001

  • CrossRef
  • Google Scholar