Effects of the Pore Morphology on Multiphase Fluid Displacement in Porous Media-A High-Resolution Modeling Investigation.

Advances in computational and image acquisition capabilities have made direct simulation of multiphase fluid flow through porous media possible. An example is application of volume of fluid modeling on images produced using the X-ray computed micro-tomography technique. Analysis of such high-resolution (both temporal and spatial) data sets provides new insights into pore-scale dynamics of previously less well-known processes. We present the outcomes of a high-resolution direct-simulation two-phase fluid displacement study performed on a series of five two-dimensional images of sandy porous media produced using erosion and dilation algorithms. This has enabled us to study the pore-scale dynamics systematically in models that are similar in connectivity but different in the morphology (pore sizes and aspect ratio). Our results show that the drainage and imbibition processes result in very distinct fluid displacement patterns in these models at the pore scale. As a result of drainage, the more open (eroded grains) models accommodate large oil clusters, while the tighter (dilated grains) models trap smaller oil clusters. The imbibition process is dominated by oil trapping in two ways: (i) bypassing larger oil clusters in the eroded models and (ii) local trapping of smaller clusters in the dilated models. This behavior is shown to arise from the relatively larger average aspect ratios of the dilated models compared to those of the eroded models. This promotes snap-off at pore throats (in competition with piston-like displacement), resulting in local trapping of the non-wetting phase. Both of these pore-scale trapping regimes seen here allow trapping of oil in as much as 50% of the pore space. Through the use of erosion and dilation operations, we show that a power-law relationship exists between the average pore size and the average grain size for these sandy media. This relationship is useful in designing engineered porous materials where pore sizes need to be estimated based on grain sizes.

Inertia Controlled Capillary Pressure at the Juncture between Converging and Uniform Channels.

In this research, we reveal the transient behavior of capillary pressure as the fluid-fluid interface travels across the juncture between a converging and uniform capillary, via high-resolution CFD (Computational Fluid Dynamics) simulations. Simulations were performed at different wetting conditions (strong-wet and intermediate-wet) and capillary wall convergence angles. Our results demonstrate that as the angle of convergence increases, capillary pressure at the junction decreases commensurately. Moreover, in contrast to strong-wet conditions, the profile of capillary pressure at the converging-uniform capillary juncture under intermediate-wet conditions is highly non-monotonic, being characterized by a parabola-like form. This non-monotonic behavior is a manifestation of strong inertial forces governing dynamic fluid-fluid interface morphology. This yields conditions that promote the advancement of the fluid-fluid interface, as inertial forces partially nullify the capillary pressure required for the immiscible interface to enter the uniform capillary. In addition to numerical analysis detailed above, a novel theoretical stability criteria that is capable of distinguishing between stable (capillary dominated) and unstable (inertia dominated) interfacial regimes at the converging-uniform capillary juncture is also proposed. In summary, this fundamental study offers new insights into the interface invasion protocol, and paves the way for the re-evaluation of capillary junction controlled interfacial dynamics.

Pore geometry control of apparent wetting in porous media.

Wettability, or preferential affinity of a fluid to a solid substrate in the presence of another fluid, plays a critical role in the statics and dynamics of fluid-fluid displacement in porous media. The complex confined geometry of porous media, however, makes upscaling of microscopic wettability to the macroscale a nontrivial task. Here, we elucidate the contribution of pore geometry in controlling the apparent wettability characteristics of a porous medium. Using direct numerical simulations of fluid-fluid displacement, we study the reversal of interface curvature in a single converging-diverging capillary, and demonstrate the co-existence of concave and convex interfaces in a porous medium-a phenomenon that we also observe in laboratory micromodel experiments. We show that under intermediate contact angles the sign of interface curvature is strongly influenced by the pore geometry. We capture the interplay between surface chemical properties and pore geometry in the form of a dimensionless quantity, the apparent wettability number, which predicts the conditions under which concave and convex interfaces co-exist. Our findings advance the fundamental understanding of wettability in confined geometries, with implications to macroscopic multiphase-flow processes in porous media, from fuel cells to enhanced oil recovery.

Suppressing viscous fingering in structured porous media.

Finger-like protrusions that form along fluid-fluid displacement fronts in porous media are often excited by hydrodynamic instability when low-viscosity fluids displace high-viscosity resident fluids. Such interfacial instabilities are undesirable in many natural and engineered displacement processes. We report a phenomenon whereby gradual and monotonic variation of pore sizes along the front path suppresses viscous fingering during immiscible displacement, that seemingly contradicts conventional expectation of enhanced instability with pore size variability. Experiments and pore-scale numerical simulations were combined with an analytical model for the characteristics of displacement front morphology as a function of the pore size gradient. Our results suggest that the gradual reduction of pore sizes act to restrain viscous fingering for a predictable range of flow conditions (as anticipated by gradient percolation theory). The study provides insights into ways for suppressing unwanted interfacial instabilities in porous media, and provides design principles for new engineered porous media such as exchange columns, fabric, paper, and membranes with respect to their desired immiscible displacement behavior.

New insights on the complex dynamics of two-phase flow in porous media under intermediate-wet conditions.

Multiphase flow in porous media is important in a number of environmental and industrial applications such as soil remediation, CO2 sequestration, and enhanced oil recovery. Wetting properties control flow of immiscible fluids in porous media and fluids distribution in the pore space. In contrast to the strong and weak wet conditions, pore-scale physics of immiscible displacement under intermediate-wet conditions is less understood. This study reports the results of a series of two-dimensional high-resolution direct numerical simulations with the aim of understanding the pore-scale dynamics of two-phase immiscible fluid flow under intermediate-wet conditions. Our results show that for intermediate-wet porous media, pore geometry has a strong influence on interface dynamics, leading to co-existence of concave and convex interfaces. Intermediate wettability leads to various interfacial movements which are not identified under imbibition or drainage conditions. These pore-scale events significantly influence macro-scale flow behaviour causing the counter-intuitive decline in recovery of the defending fluid from weak imbibition to intermediate-wet conditions.

Effects of intermediate wettability on entry capillary pressure in angular pores.

Entry capillary pressure is one of the most important factors controlling drainage and remobilization of the capillary-trapped phases as it is the limiting factor against the two-phase displacement. It is known that the entry capillary pressure is rate dependent such that the inertia forces would enhance entry of the non-wetting phase into the pores. More importantly the entry capillary pressure is wettability dependent. However, while the movement of a meniscus into a strongly water-wet pore is well-defined, the invasion of a meniscus into a weak or intermediate water-wet pore especially in the case of angular pores is ambiguous. In this study using OpenFOAM software, high-resolution direct two-phase flow simulations of movement of a meniscus in a single capillary channel are performed. Interface dynamics in angular pores under drainage conditions have been simulated under constant flow rate boundary condition at different wettability conditions. Our results shows that the relation between the half corner angle of pores and contact angle controls the temporal evolution of capillary pressure during the invasion of a pore. By deviating from pure water-wet conditions, a dip in the temporal evolution of capillary pressure can be observed which will be pronounced in irregular angular cross sections. That enhances the pore invasion with a smaller differential pressure. The interplay between the contact angle and pore geometry can have significant implications for enhanced remobilization of ganglia in intermediate contact angles in real porous media morphologies, where pores are very heterogeneous with small shape factors.