Geometric and Topological Analysis of Three-Dimensional Porous Media: Pore Space Partitioning Based on Morphological Skeletonization

Abstract

This article presents a versatile, rigorous, and efficient methodology for extracting various geometric and topological parameters of 3D discrete porous media. The new approach takes advantage of the morphological skeleton of the pore structure—a lower dimensional representation of the pore space akin to the topological “deformation retract”. The skeleton is derived by a fully parallel thinning algorithm that fulfils two essential requirements: it generates a medial axis and preserves the connectivity of the pore space. Topological analysis is accomplished by classifying all skeleton points as node or link (branch) points according to the concept of λ-adjacency in 3D discrete space. In this manner, node coordination number and link length distributions are directly obtained from the skeleton. Pore necks (throats) are identified through a search for minima in the hydraulic radius of individual pore space channels outlined by skeleton links. In addition to the determination of the size distribution of the constrictions (pore necks) that control nonwetting phase invasion, improved estimates of the distributions of effective hydraulic and electric conductivity of individual pore space channels are obtained. Furthermore, erection of planes at the location of pore necks results in partitioning of the pore space into its constituent pores. This enables the characterization of the pore space in terms of a pore volume distribution. The new methodology is illustrated by application to a regular cubic pore network and irregularly shaped 2D and 3D pore networks generated by stochastic simulation. In the latter case, important new results are obtained concerning the sensitivity of geometric and topological properties of the microstructure to the parameters of stochastic simulation, namely, the porosity and correlation function. It is found that model porous media reconstructed from the same porosity and correlation function can exhibit marked differences in geometry and connectivity, which correlate with differences in specific surface area.