Pore Structure and Multifractal Characteristics of the Upper Lianggaoshan Formation in the Northeastern Sichuan Basin, China

The Upper Lianggaoshan (LGS) Formation in the northeastern Sichuan Basin, composed of shale with interbedded siltstone, is a promising target layer for shale oil. Accurate evaluation of pore structures is essential for effective exploration of shale oil. This study investigated pore structures of siltstone and shale samples from the Upper LGS Formation using low-pressure CO2 adsorption (LTCA), low-temperature N2 adsorption (LTNA), high-pressure mercury intrusion (HPMI), and nuclear magnetic resonance (NMR) methods. The single-exponent and multifractal dimensions of samples were determined, and the relationships between fractal dimensions and pore structures were explored. Results show that the pore size distribution (PSD) of siltstone and shale samples exhibits multi-peak characteristics, with mesopores (2–50 nm) being dominant in the total pore volumes. The multi-scaled pores in shale and siltstone samples exhibit fractal characteristics. The average values of single-fractal dimensions (D1, D2) obtained by LTNA data are 2.39 and 2.62 for shale samples, and 2.24 and 2.59 for siltstone samples, respectively. Compared to siltstones, the pore structures of shale samples exhibit greater complexity, indicated by larger fractal dimensions. The samples from subsections Liang 2 and Liang 3 exhibit greater heterogeneity compared to subsection Liang 1. The single-fractal dimensions of micropores and mesopores show positive correlations with specific surface area (SSA) and pore volume (PV), while the fractal dimension of macropores shows a negative correlation with average pore diameter and median radius. The average values of single-fractal dimension D3 obtained from HPMI data are 2.9644 and 2.9471 for shale and siltstone samples, respectively, indicating more complex structures of macropores in shale samples compared to siltstone samples. The average value of ΔDNMR and singularity strength range Δα obtained by a multifractal model for core samples from subsection Liang 1 are 1.868 and 2.155, respectively, which are the smallest among all of the three subsections, indicating that the heterogeneity of pore structures of subsection Liang 1 is the weakest. This research provides valuable guidance for shale oil development in the northeastern Sichuan Basin, China.