Flux-dependent percolation transition in immiscible two-phase flows in porous media

Development of stochastically reconstructed 3D porous media micromodels using additive manufacturing: numerical and experimental validation

We propose an integrated methodology for the design and fabrication of 3D micromodels that are suitable for the pore-scale study of transport processes in macroporous materials. The micromodels, that bear the pore-scale characteristics of sandstone, such as porosity, mean pore size, etc, are designed following a stochastic reconstruction algorithm that allows for fine-tuning the porosity and the correlation length of the spatial distribution of the solid material. We then construct a series of 3D micromodels at very fine resolution (i.e. 16 μ m) using a state-of-the-art 3D printing infrastructure, specifically a ProJet MJP3600 3D printer, that utilizes the Material Jetting technology. Within the technical constraints of the 3D printer resolution, the fabricated micromodels represent scaled-up replicas of natural sandstones, that are suitable for the study of the scaling between the permeability, the porosity and the mean pore size. The REV- and pore-scale characteristics of the resulting physical micromodels are recovered using a combination of X-ray micro-CT and microfluidic studies. The experimental results are then compared with single-phase flow simulations at pore-scale and geostatistic models in order to determine the effects of the design parameters on the intrinsic permeability and the spatial correlation of the velocity profile. Our numerical and experimental measurements reveal an excellent match between the properties of the designed and fabricated 3D domains, thus demonstrating the robustness of the proposed methodology for the construction of 3D micromodels with fine-tuned and well-controlled pore-scale characteristics. Furthermore, a pore-scale numerical study over a wider range of 3D digital domain realizations reveals a very good match of the measured permeabilities with the predictions of the Kozeny-Carman formulation based on a single control parameter, k 0 , that is found to have a practically constant value for porosities ϕ ≥ 0.2 . This, in turn, enables us to customize the sample size to meet REV constraints, including enlarging pore morphology while considering the Reynolds number. It is also found that at lower porosities there is a significant increase in the fraction of the non-percolating pores, thus leading to different k 0 , as the porosity approaches a numerically determined critical porosity value, ϕ c , where the domain is no longer percolating.

Effective Rheology of Two-Phase Flow in Three-Dimensional Porous Media: Experiment and Simulation

We present an experimental and numerical study of immiscible two-phase flow of Newtonian fluids in three-dimensional (3D) porous media to find the relationship between the volumetric flow rate (Q) and the total pressure difference ([Formula: see text]) in the steady state. We show that in the regime where capillary forces compete with the viscous forces, the distribution of capillary barriers at the interfaces effectively creates a yield threshold ([Formula: see text]), making the fluids reminiscent of a Bingham viscoplastic fluid in the porous medium. In this regime, Q depends quadratically on an excess pressure drop ([Formula: see text]). While increasing the flow rate, there is a transition, beyond which the overall flow is Newtonian and the relationship is linear. In our experiments, we build a model porous medium using a column of glass beads transporting two fluids, deionized water and air. For the numerical study, reconstructed 3D pore networks from real core samples are considered and the transport of wetting and non-wetting fluids through the network is modeled by tracking the fluid interfaces with time. We find agreement between our numerical and experimental results. Our results match with the mean-field results reported earlier.

Ensemble distribution for immiscible two-phase flow in porous media

We construct an ensemble distribution to describe steady immiscible two-phase flow of two incompressible fluids in a porous medium. The system is found to be ergodic. The distribution is used to compute macroscopic flow parameters. In particular, we find an expression for the overall mobility of the system from the ensemble distribution. The entropy production at the scale of the porous medium is shown to give the expected product of the average flow and its driving force, obtained from a black-box description. We test numerically some of the central theoretical results.

Oil displacement through a porous medium with a temperature gradient

We investigate the effect of a temperature gradient on oil recovery in a two-dimensional pore-network model. The oil viscosity depends on temperature as μ(o) [Please see text] e(B/T), where B is a physicochemical parameter, depending on the type of oil, and T is the temperature. A temperature gradient is applied across the medium in the flow direction. Initially, the porous medium is saturated with oil, and then another fluid is injected. We have considered two cases representing different injection strategies. In the first case, the invading fluid viscosity is constant (finite viscosity ratio), while in the second one, the invading fluid is inviscid (infinite viscosity ratio). Our results show that for the case of finite viscosity ratio, recovery increases with ΔT independent of strength or sign of the gradient. For an infinite viscosity ratio, a positive temperature gradient is necessary to enhance recovery. Moreover, we show that for ΔT>0, the percentage of oil recovery generally decreases (increases) with B for a finite (infinite) viscosity ratio. Finally, we also extend our results for infinite viscosity ratio to a three-dimensional porous media geometry.