Influence of phase connectivity on the relationship among capillary pressure, fluid saturation, and interfacial area in two-fluid-phase porous medium systems

A Dataset of 3D Structural and Simulated Transport Properties of Complex Porous Media

Physical processes that occur within porous materials have wide-ranging applications including - but not limited to - carbon sequestration, battery technology, membranes, oil and gas, geothermal energy, nuclear waste disposal, water resource management. The equations that describe these physical processes have been studied extensively; however, approximating them numerically requires immense computational resources due to the complex behavior that arises from the geometrically-intricate solid boundary conditions in porous materials. Here, we introduce a new dataset of unprecedented scale and breadth, DRP-372: a catalog of 3D geometries, simulation results, and structural properties of samples hosted on the Digital Rocks Portal. The dataset includes 1736 flow and electrical simulation results on 217 samples, which required more than 500 core years of computation. This data can be used for many purposes, such as constructing empirical models, validating new simulation codes, and developing machine learning algorithms that closely match the extensive purely-physical simulation. This article offers a detailed description of the contents of the dataset including the data collection, simulation schemes, and data validation.

Capillary equilibrium of bubbles in porous media

In geologic, biologic, and engineering porous media, bubbles (or droplets, ganglia) emerge in the aftermath of flow, phase change, or chemical reactions, where capillary equilibrium of bubbles significantly impacts the hydraulic, transport, and reactive processes. There has previously been great progress in general understanding of capillarity in porous media, but specific investigation into bubbles is lacking. Here, we propose a conceptual model of a bubble's capillary equilibrium associated with free energy inside a porous medium. We quantify the multistability and hysteretic behaviors of a bubble induced by multiple state variables and study the impacts of pore geometry and wettability. Surprisingly, our model provides a compact explanation of counterintuitive observations that bubble populations within porous media can be thermodynamically stable despite their large specific area by analyzing the relationship between free energy and bubble volume. This work provides a perspective for understanding dispersed fluids in porous media that is relevant to CO₂ sequestration, petroleum recovery, and fuel cells, among other applications.

Pore-Scale Modeling of Fluid-Fluid Interfacial Area in Variably Saturated Porous Media Containing Microscale Surface Roughness

A pore-scale model is developed to simulate fluid-fluid interfacial area in variably saturated porous media, with a specific focus on incorporating the effects of solid-surface roughness. The model is designed to quantify total (film and meniscus) fluid-fluid interfacial area (A_nw ) over the full range of wetting-phase fluid saturation (S_w ) based on the inherent properties of the porous medium. The model employs a triangular pore space bundle of cylindrical capillaries (BCC) framework, modified with three surface roughness-related parameters. The first parameter (surface roughness factor) represents the overall magnitude of surface roughness, whereas the other two parameters (interface growth factor and critical adsorptive film thickness) reflect the micro-scale structure of surface roughness. A series of sensitivity analyses was conducted for the controlling variables, and the efficacy of the model was tested using air-water interfacial area data measured for three natural porous media. The model produced good simulations of the measured A_nw data over the full range of saturation. The results demonstrate that total interfacial areas for natural media are typically much larger than those for ideal media comprising smooth surfaces due to the substantial contribution of surface roughness to wetting-film interfacial area. The degree to which fluid-fluid interfacial area is influenced by roughness is a function of fluid-retention characteristics and the nature of the rough surfaces. The full impact of roughness may be masked to some degree due to the formation of thick wetting films, which is explicitly quantified by the model. Application of the model provides insight into the importance of the interplay between pore-scale distribution and configuration of wetting fluid and the surface properties of solids.

Toward a New Generation of Two-Fluid Flow Models Based on the Thermodynamically-Constrained Averaging Theory

Traditional models of two-fluid flow through porous media at the macroscale have existed for nearly a century. These phenomenological models are not firmly connected to the microscale; thermodynamic constraints are not enforced; empirical closure relations are well known to be hysteretic; fluid pressures are typically assumed to be in a local equilibrium state with fluid saturations; and important quantities such as interfacial and curvilinear geometric extents, tensions, and curvatures, known to be important from microscale studies, do not explicitly appear in traditional macroscale models. Despite these shortcomings, the traditional model for two-fluid flow in porous media has been extensively studied to develop efficient numerical approximation methods, experimental and surrogate measure parameterization approaches, and convenient pre- and post-processing environments; and they have been applied in a large number of applications from a variety of fields. The thermodynamically constrained averaging theory (TCAT) was developed to overcome the limitations associated with traditional approaches, and we consider here issues associated with the closure of this new generation of models. It has been shown that a hysteretic-free state equation exists based upon integral geometry that relates changes in volume fractions, capillary pressure, interfacial areas, and the Euler characteristic. We show an analysis of how this state equation can be parameterized with a relatively small amount of data. We also formulate a state equation for resistance coefficients that we show to be hysteretic free, unlike traditional relative permeability models. Lastly, we comment on the open issues remaining for this new generation of models.