A connected-component-labeling-based approach to virtual porosimetry

Analyzing the pore-size distribution of porous materials, made up of an aggregation of interconnected pores, is a demanding task. Mercury intrusion porosimetry (MIP) is a physical method that intrudes mercury into a sample at increasing pressures to obtain a pore-size histogram. This method has been simulated in-silice with several approaches requiring prior computation of a skeleton. We present a new approach to simulate MIP that does not require skeleton computation. Our method is an iterative process that considers the diameters corresponding to pressures. At each iteration, geometric tests detect throats for the corresponding diameter and a CCL process collects the region invaded by the mercury. Additionally, a new decomposition model called CUDB, is used. This is suitable for computing the throats and performs better with the CCL algorithm than a voxel model. Our approach obtains the pore-size distribution of the porous medium, and the corresponding pore graph.