Capillary pressure at irregularly shaped pore throats: Implications for water retention characteristics

The random shapes of pore throats in geomaterials hinder accurate estimation of capillary pressure, and conventional pore network models that simply use the Young–Laplace equation assuming circular pore throats overestimate the capillary pressure. As a solution to this problem that does not complicate the pore network model or slow its implementation, we propose a new morphological analysis method to correlate the capillary pressure at an irregular pore channel with its cross-sectional geometry using lattice Boltzmann (LB) simulation and Mayer and Stowe–Princen theory. Geometry-based shape factors for pore throats are shown here to correlate strongly with the capillary pressure obtained by LB simulation. Water retention curves obtained by incorporating the morphological calibration into conventional pore network simulation and their correlative scheme agree well with experimental data. The suggested method is relevant to pore-scale processes such as geological CO2 sequestration, methane bubbling from wetlands, and enhanced carbon recovery.