Calculus with combinatorial differential forms for fluid flow analysis in porous and fractured media

TL;DR

This paper introduces a novel mathematical framework using combinatorial differential forms to model and analyze fluid flow in complex porous and fractured media with variable geometries, based on experimental fabric data.

Contribution

The work develops a flexible, structure-preserving method for representing multi-dimensional voids in porous media and formulates conservation laws directly in matrix form for computational analysis.

Findings

Validated with XCT images of four rocks

Accurately models diverse void geometries

Enables precise computation of transport properties

Abstract

The fabric of porous and fractured media contains solid regions (grains) and voids. The space conducting fluids is a system of connected voids with variable geometries. Relative to the grain sizes, the voids can be voluminous with three comparable large extensions, narrow expansive with two comparable large extensions and one smaller extension, and thin long with one comparable large extension and two smaller extensions. The widely used representation of void spaces by systems of spheres connected by cylinders (pore network models) is an acceptable approximation for some special cases, but not for most porous and fractured media. We propose a flexible method for modelling such media by mapping their measured fabric's characteristics - void and grain volume distributions and shapes - onto polyhedral tessellations of space. The map assigns voluminous voids and grains to polyhedrons (3D), narrow expansive voids to some polyhedral faces (2D), and thin long voids to some polyhedral edges (1D), as dictated by experimental data. The analysis of transport through such discrete structures with components of different dimensions is performed by a novel mathematical method, which uses combinatorial differential forms to represent physical properties and their fluxes, as well as structure-preserving operators on such forms to formulate the conservation laws exactly and directly in matrix form, ready for computation. The method allows for individual material properties, such as conductivity, to be assigned to voids of all dimensions, so that the three types of voids are suitably represented. Publicly available XCT images of four different rocks are used to test the method.