Characterization of Pore Heterogeneity in Lacustrine Shale Based on MIP, LTNA, NMR, and Multifractal Characteristics: A Case Study of the Jurassic Dongyuemiao Member, China

Pore structure plays a critical role in evaluating shale “sweet spots”. Compared to marine shale, lacustrine shale has more diverse lithofacies types and greater heterogeneity in pore structure due to frequently changing environmental conditions. Using methods such as mercury intrusion porosimetry (MIP), field emission scanning electron microscopy (FE-SEM), nuclear magnetic resonance (NMR), and X-ray diffraction (XRD), this study investigates the micropore structures and heterogeneity of different lithofacies in the Jurassic Dongyuemiao Member lacustrine shale. Image processing and multifractal theory were employed to identify the controlling factors of pore structure heterogeneity. The key findings are as follows. (1) Based on mineral content and laminae types, the lithofacies types of Dongyuemiao lacustrine shale are classified into four types: shell–laminae mixed shale (SLMS), silty–laminae clay shale (SLCS), clast–laminae clay shale (CLCS), and clay shale (CS). (2) Based on genesis, shale reservoirs’ pore and permeability space are categorized into inorganic pores, organic pores, and micro-fractures. Inorganic pores consist of inter-particle pores and intra-particle pores. Pore size distribution curves for all four lithofacies exhibit two main peaks, with pore sizes concentrated in the ranges of 2–10 nm and 50–80 nm. Mesopores and macropores dominate, accounting for over 80% of the total pore volume. Mesopores are most developed in CLCS, representing 56.3%. (3) Quartz content is positively correlated with the multifractal dimension, while clay content shows a negative correlation. Higher quartz content, coupled with lower clay content, weakens pore structure heterogeneity. A negative correlation exists between total organic carbon (TOC) and the multifractal dimension, indicating that higher organic matter content enhances organic pore development and increases microscopic heterogeneity. (4) Porosity heterogeneity in SLMS is effectively characterized by D0-Dmax, while in the other three lithofacies, it is characterized by Dmin-D0. Permeability across all lithofacies correlates with D0-Dmax. In CS, SLMS, and SLCS, permeability is positively correlated with D0-Dmax, with higher values indicating greater permeability heterogeneity. In CLCS, permeability is negatively correlated with D0-Dmax, such that lower values reflect stronger heterogeneity.