No Arabic abstract
Immiscible fluid displacement in porous media is fundamental for many environmental processes, including infiltration of water in soils, groundwater remediation, enhanced recovery of hydrocarbons and carbon geosequestration. Microstructural heterogeneity, in particular of particle sizes, can significantly impact immiscible displacement. For instance, it may lead to unstable flow and preferential displacement patterns. We present a systematic, quantitative pore-scale study of the impact of spatial correlations in particle sizes on the drainage of a partially-wetting fluid. We perform pore-network simulations with varying flow rates and different degrees of spatial correlation, complemented with microfluidic experiments. Simulated and experimental displacement patterns show that spatial correlation leads to more preferential invasion, with reduced trapping of the defending fluid, especially at low flow rates. Numerically, we find that increasing the correlation length reduces the fluid-fluid interfacial area and the trapping of the defending fluid, and increases the invasion pattern asymmetry and selectivity. Our experiments, conducted for low capillary numbers, support these findings. Our results delineate the significant effect of spatial correlations on fluid displacement in porous media, of relevance to a wide range of natural and engineered processes.
In this paper we present a new model for modeling the diffusion and relative dispersion of particles in homogeneous isotropic turbulence. We use an Heisenberg-like Hamiltonian to incorporate spatial correlations between fluid particles, which are modeled by stochastic processes correlated in time. We are able to reproduce the ballistic regime in the mean squared displacement of single particles and the transition to a normal diffusion regime for long times. For the dispersion of particle pairs we find a $t^{2}$-dependence of the mean squared separation at short times and a $t$-dependence for long ones. For intermediate times indications for a Richardson $t^{3}$ law are observed in certain situations. Finally the influence of inertia of real particles on the dispersion is investigated.
We present a theoretical framework for immiscible incompressible two-phase flow in homogeneous porous media that connects the distribution of local fluid velocities to the average seepage velocities. By dividing the pore area along a cross-section transversal to the average flow direction up into differential areas associated with the local flow velocities, we construct a distribution function that allows us not only to re-establish existing relationships between the seepage velocities of the immiscible fluids, but also to find new relations between their higher moments. We support and demonstrate the formalism through numerical simulations using a dynamic pore-network model for immiscible two-phase flow with two- and three-dimensional pore networks. Our numerical results are in agreement with the theoretical considerations.
We develop a theory for the problem of high pressure air injection into deep reservoirs containing light oil. Under these conditions, the injected fluid (oxygen + inert components) is completely miscible with the oil in the reservoir. Moreover, exothermic reactions between dissolved oxygen and oil are possible. We use Kovals model to account for the miscibility of the components, such that the fractional-flow functions resemble the ones from Buckley-Leverett flow. This allows to decompose the solution of this problem into a series of waves. We then proceed to obtain full analytical solutions in each wave. Of particular interest {is the case where} the combustion wave presents a singularity in its internal wave profile. Evaluation of the variables of the problem at the singular point determines the macroscopic parameters of the wave, i.e., combustion temperature, wave speed and downstream oil fraction. The waves structure was observed previously for reactive immiscible displacement and we describe it here for the first time for reactive miscible displacement of oil. We validate the developed theory using numerical simulations.
Relative permeability theory for immiscible two-phase flow in porous media assumes a linear dependency of the seepage velocity of each fluid on the pressure gradient. This implies that the average fluid velocity also exhibits such a linear dependence. Recent experimental, computational and theoretical work, however, show that the average flow velocity follows a power law in the the pressure gradient with an exponent in the range larger than one up to two over a wide range of parameters. Such a behavior is incompatible with relative permeability theory. A recent theory based on Euler homogeneity of the volumetric flow rates of the fluids generalizes relative permeability theory in such a way that it is capable of handling this non-linear behavior. A central quantity in this theory is the co-moving velocity which is related to, but not equal to the difference between the seepage velocities of the fluids. In order to close the equation set that ensues from the theory, a constitutive equation has to be supplied for the co-moving velocity. We construct this constitutive equation from relative permeability data in the literature. It turns out to exhibit a remarkably simple form when expressed in the right variables. We follow this analysis up by simulating immiscible two-phase flow using a dynamic pore network model finding the same results as those based on the relative permeability data.
Imbibition, the displacement of a nonwetting fluid by a wetting fluid, plays a central role in diverse energy, environmental, and industrial processes. While this process is typically studied in homogeneous porous media with uniform permeabilities, in many cases, the media have multiple parallel strata of different permeabilities. How such stratification impacts the fluid dynamics of imbibition, as well as the fluid saturation after the wetting fluid breaks through to the end of a given medium, is poorly understood. We address this gap in knowledge by developing an analytical model of imbibition in a porous medium with two parallel strata, combined with a pore network model that explicitly describes fluid crossflow between the strata. By numerically solving these models, we examine the fluid dynamics and fluid saturation left after breakthrough. We find that the breakthrough saturation of nonwetting fluid is minimized when the imposed capillary number Ca is tuned to a value Ca$^*$ that depends on both the structure of the medium and the viscosity ratio between the two fluids. Our results thus provide quantitative guidelines for predicting and controlling flow in stratified porous media, with implications for water remediation, oil/gas recovery, and applications requiring moisture management in diverse materials.