Dual models of pore spaces

We present a model for pore spaces that consists of two parts related by duality: (1) a decomposition of an open polyhedral pore space into open contractible pore bodies separated by relatively open interfaces and (2) a pore network that is homotopy equivalent to the pore space. The dual model is unique and free of parameters, but it relies on regularity conditions for the pore space. We show how to approximate any pore space by the interior of a polyhedral complex such that the regularity conditions are fulfilled. Thus, we are able to calculate the dual model from synthetic porous media and images of real porous media. The pore bodies are unions of relatively open Delaunay cells with respect to the corners of the pore boundary, and the pore network consists of certain at most two-dimensional (2D) Voronoi cells with respect to the corners of the pore boundary. The pore network describes the neighborhood relations between the pore bodies. In particular, any relatively open 2D Delaunay face f separating two pore bodies has a unique (relatively open) dual network edge. In our model, f is a pore throat only if it is hit by its dual network edge. Thus, as opposed to widespread intuition, any pore throat is convex, and adjacent pore bodies are not necessarily separated by pore throats. Due to the duality between the pore network and the decomposition of the pore space into pore bodies it is straightforward to store the geometrical properties of the pore bodies [pore throats] as attributes of the dual network vertices [edges]. Such an attributed network is used to perform 2D drainage simulations. The results agree very well with those from a pore-morphology based modeling approach performed directly on the digital image of a porous medium. Contractibility of the pore bodies and homotopy equivalence of the pore space and the pore network is proven using discrete Morse theory and the nerve theorem from combinatorial topology.