A new approach for the characterization of the pore structure of dual porosity rocks

For the realistic representation of the pore space of dual porosity rocks, a new method of pore structure characterization is developed by combining experimental Hg intrusion/retraction curves with back-scattered scanning electron microscope (BSEM) images and inverse modeling algorithms. The pore space autocorrelation function measured by processing the digitized BSEM images is combined with the surface fractal dimension estimated from the high pressure Hg intrusion (MIP) data to derive a synthetic small-angle neutron scattering (SANS) intensity function, the inversion of which provides a volume-based pore body radius distribution (PBRD). The volume-based PBRD is fitted with a multimodal number-based PBRD consisting of two component distributions: one representing the macroporosity and another one representing the microporosity. Based on arguments of percolation theory, analytical mathematical models are developed to describe the Hg intrusion in and retraction from dual pore networks in terms of the complete PBRD, pore throat radius distribution (PTRD) of macroporosity, drainage accessibility functions (DAFs) of both porosities, and imbibition accessibility functions (IAFs) of both porosities. Inverse modeling of the Hg intrusion data set enables us to estimate the PTRD and DAFs. Inverse modeling of the Hg retraction datasets enables us to estimate a set of primary and secondary IAFs. The method is demonstrated by the pore structure characterization of four outcrop samples of carbonate and sandstone rocks. Analytic approximate equations developed from the critical path analysis (CPA) of percolation theory enable us to calculate explicitly the absolute permeability and the formation factor of the porous rocks using the estimated parameters (PBRD, PTRD, DAF) of the macroporosity. The measured permeability of cores is predicted satisfactorily and observed discrepancies may be attributed to large length-scale macro-heterogeneities which are not evident in BSEM images and Hg porosimetry data.