Foam Drainage in Porous Media

In this paper we present a simple analysis of liquid drainage in foams confined in porous media. First we derive the equation for the evolution of the liquid saturation using general mass and momentum conservation arguments and phenomenological relations between the transport parameters and liquid saturation. We find an unusual foam drainage equation in which the determinant terms express the competition between the external force field, represented here by the gravity field, and capillary pressure gradient. We present analytical solutions of the drainage equation in three cases: (a) gravity forces are dominant over capillary forces, (b) capillary forces are dominant over gravity forces, and (c) capillary and gravity forces are comparable in order of magnitude.     

The Spreading of Immiscible Fluids in Porous Media under the Influence of Gravity

We use an approach based on invasion percolation in a gradient (IPG) to describe the displacement patterns that develop when a fluid spreads on an impermeable boundary in a porous medium under the influence of gravity (buoyancy) forces in a drainage process. The approach is intended to simulate applications, such as the spreading of a DNAPL in the saturated zone and of a NAPL in the vadose zone on top of an impermeable layer, or the classical problems of gravity underruning and gravity override in reservoir engineering. As gravity acts in a direction transverse to the main displacement direction, a novel form of IPG develops. We study numerically the resulting patterns for a combination of transverse and parallel Bond numbers and interpret the results using the concepts of gradient percolation. A physical interpretation in terms of the capillary number, the viscosity ratio and the gravity Bond number is also provided. In particular, we consider the scaling of the thickness of the spreading gravity tongue, for the cases of gravitydominated and viscousunstable displacements, and of the propagating front in the case of stabilized displacement at relatively high rates. It is found that the patterns have percolation (namely fractallike) characteristics, which cannot be captured by conventional continuum equations. These characteristics will affect, for example, mass transfer and must be considered in the design of remediation processes.     

Coupling of Foam Drainage and Viscous Fingering in Porous Media Revealed by X-ray Computed Tomography

We investigate the development foam in granular porous media and the subsequent flow of the surfactant solution, where the fluid fraction variations are visualized and quantified using X-ray computed tomography. It is found that foam flows in a front like manner leading to a residual liquid fraction of 0.18±0.01, far from the inlet surface of the porous sample. A desaturation backward wave is also observed during foam development. We provided direct evidence that the flow of surfactant solution in porous media containing foam gives rise to superposition of a drainage wave and a characteristic viscous fingering pattern. In the wave the liquid fraction ranges from the above residual value to nearly 0.25±0.01. The liquid fraction associated with the viscous fingering decays as a function of distance but the inlet value increases up to 0.06±0.01. Certain ideas about the physics of foam flow in porous media are revised in the light of our findings.     

1-D Modeling of Hydrate Depressurization in Porous Media

A thermal, three-phase, one-dimensional numerical model is developed to simulate two regimes of gas production from sediments containing methane hydrates by depressurization: the dissociation-controlled regime and the flow-controlled regime. A parameter namely dissociation-flow time-scale ratio, R, is defined and employed to identify the two regimes. The numerical model uses a finite-difference scheme; it is implicit in water and gas saturations, pressure and temperature, and explicit in hydrate saturation. The model shows that laboratory-scale experiments are often dissociation-controlled, but the field-scale processes are typically flow-controlled. Gas production from a linear reservoir is more sensitive to the heat transfer coefficient with the surrounding than the longitudinal heat conduction coefficient, in 1-D simulations. Gas production is not very sensitive to the well temperature boundary condition. This model can be used to fit laboratory-scale experimental data, but the dissociation rate constant, the multiphase flow parameters and the heat transfer parameters are uncertain and should be measured experimentally.     

Percolation Thresholds in 2-Dimensional Prefractal Models of Porous Media*

Considerable effort has been directed towards the application of percolation theory and fractal modeling to porous media. We combine these areas of research to investigate percolation in prefractal porous media. We estimated percolation thresholds in the pore space of homogeneous random 2-dimensional prefractals as a function of the fractal scale invariance ratio b and iteration level i. The percolation thresholds for these simulations were found to increase beyond the 0.5927l... porosity expected in Bernoulli (uncorrelated) percolation networks. Percolation in prefractals occurs through large pores connected by small pores. The thresholds increase with both b (a finite size effect) and i. The results allow the prediction of the onset of percolation in models of prefractal porous media and can be used to bound modeling efforts. More fundamental applications are also possible. Only a limited range of parameters has been explored empirically but extrapolations allow the critical fractal dimension to be estimated for a large combination of b and i values. Extrapolation to infinite iterations suggests that there may be a critical fractal dimension of the solid at which the pore space percolates. The extrapolated value is close to 1.89 – the well-known fractal dimension of percolation clusters in 2-dimensional Bernoulli networks.     

Pore Network Analysis of Resistivity Index for Water-Wet Porous Media

Pore network analysis is used to investigate the effects of microscopic parameters of the pore structure such as pore geometry, pore-size distribution, pore space topology and fractal roughness porosity on resistivity index curves of strongly water-wet porous media. The pore structure is represented by a three-dimensional network of lamellar capillary tubes with fractal roughness features along their pore-walls. Oil-water drainage (conventional porous plate method) is simulated with a bond percolation-and-fractal roughness model without trapping of wetting fluid. The resistivity index, saturation exponent and capillary pressure are expressed as approximate functions of the pore network parameters by adopting some simplifying assumptions and using effective medium approximation, universal scaling laws of percolation theory and fractal geometry. Some new phenomenological models of resistivity index curves of porous media are derived. Finally, the eventual changes of resistivity index caused by the permanent entrapment of wetting fluid in the pore network are also studied.Resistivity index and saturation exponent are decreasing functions of the degree of correlation between pore volume and pore size as well as the width of the pore size distribution, whereas they are independent on the mean pore size. At low water saturations, the saturation exponent decreases or increases for pore systems of low or high fractal roughness porosity respectively, and obtains finite values only when the wetting fluid is not trapped in the pore network. The dependence of saturation exponent on water saturation weakens for strong correlation between pore volume and pore size, high network connectivity, medium pore-wall roughness porosity and medium width of the pore size distribution. The resistivity index can be described succesfully by generalized 3-parameter power functions of water saturation where the parameter values are related closely with the geometrical, topological and fractal properties of the pore structure.     

Some Considerations on Modeling Heat and Mass Transfer in Porous Media

In this paper some considerations are presented about the equations needed to set up a model of the process of heat and mass transfer in porous media. A clear classification is made of the various types of equations used and of their physical meaning. Special attention is paid to the thermodynamic equilibrium equations and to their derivation since they are too often taken for granted. The importance of the various transport mechanisms (of mass and energy) is analyzed and the consequences that can arise when some term is neglected are indicated.     

The Physical Origin of Interfacial Coupling in Two-Phase Flow through Porous Media

Recently developed transport equations for two-phase flow through porous media usually have a second term that has been included to account properly for interfacial coupling between the two flowing phases. The source and magnitude of such coupling is not well understood. In this study, a partition concept has been introduced into Kalaydjian's transport equations to construct modified transport equations that enable a better understanding of the role of interfacial coupling in two-phase flow through natural porous media. Using these equations, it is demonstrated that, in natural porous media, the physical origin of interfacial coupling is the capillarity of the porous medium, and not interfacial momentum transfer, as is usually assumed. The new equations are also used to show that, under conditions of steady-state flow, the magnitude of mobilities measured in a countercurrent flow experiment is the same as that measured in a cocurrent flow experiment, contrary to what has been reported previously. Moreover, the new equations are used to explicate the mechanism by which a saturation front steepens in an unstabilized displacement, and to show that the rate at which a wetting fluid is imbibed into a porous medium is controlled by the capillary coupling parameter, . Finally, it is argued that the capillary coupling parameter, , is dependent, at least in part, on porosity. Because a clear understanding of the role played by interfacial coupling is important to an improved understanding of two-phase flow through porous media, the new transport equations should prove to be effective tools for the study of such flow.     

Effect of Wetting on Gravity Drainage in Porous Media

A series of benchmark experiments on the effect of the wetting state on the flow properties in porous media were performed, allowing us to relate the wetting properties at the pore scale to the macroscale hydrodynamics. Drainage of n-alkanes (oils) displaced by air in a model porous medium consisting of water-wet sand was studied using gamma-ray densitometry and weight measurements. The enormous advantage of our system is that we know and control the wetting properties perfectly: we can tune the wetting properties by changing the salinity of the water. This allows us to perform porous medium flow experiments for the different wetting states without changing the transport properties (viscosity and density) of the oil. Drainage is found to be more efficient, and consequently oil recovery more important for partial wetting.     

Simulation of capillary-dominated displacements in microtomographic images of reservoir rocks

Displacement simulation in realistic pore networks, such as those derived from X-ray microtomography, is presented for the regime where capillarity controls fluid motions and spatial distributions. Complex displacement sequences involving both imbibition and drainage are constructed to extract wettability indices. The percolation properties of predicted equilibrium phase distributions are analyzed. Equilibrium fluid distributions are used to model transport properties for each phase.     

The Role of Capillarity in Two-Phase Flow through Porous Media

When determining experimentally relative permeability and capillary pressure as a function of saturation, a self-consistent system of macroscopic equations, that includes Leverett's equation for capillary pressure, is required. In this technical note, such a system of equations, together with the conditions under which the equations apply, is formulated. With the aid of this system of equations, it is shown that, at the inlet boundary of a vertically oriented porous medium, static conditions pertain, and that potentials, because of the definition of potential, are equal in magnitude to pressures. Consequently, Leverett's equation is valid at the inlet boundary of the porous medium, provided cocurrent flow, or gravity-driven, countercurrent flow is taking place, and provided the porous medium is homogeneous. Moreover, it is demonstrated that Leverett's equation is valid for flow along the length of a vertically oriented porous medium, provided cocurrent flow, or gravity-driven, countercurrent flow is taking place, and provided the porous medium is homogeneous and there are no hydrodynamic effects. However, Leverett's equation is invalid for horizontal, steady-state, forced, countercurrent flow. When such flow is taking place, it is the sum of the pressures, and not the difference in pressures, which is related to capillary pressure.     

Electrical Conductivity and Percolation Aspects of Statistically Homogeneous Porous Media

A method of 3-D stochastic reconstruction of porous media based on statistical information extracted from 2-D sections is evaluated with reference to the steady transport of electric current. Model microstructures conforming to measured and simulated pore space autocorrelation functions are generated and the formation factor is systematically determined by random walk simulation as a function of porosity and correlation length. Computed formation factors are found to depend on correlation length only for small values of this parameter. This finding is explained by considering the general percolation behavior of a statistically homogeneous system. For porosities lower than about 0.2, the dependence of formation factor on porosity shows marked deviations from Archie's law. This behavior results from the relatively high pore space percolation threshold (0.09) of the simulated media and suggests a limitation to the applicability of the method to low porosity media. It is additionally demonstrated that the distribution of secondary porosity at a larger scale can be simulated using stochastic methods. Computations of the formation factor are performed for model media with a matrix-vuggy structure as a function of the amount and spatial distribution of vuggy porosity and matrix conductivity. These results are shown to be consistent with limited available experimental data for carbonate rocks.     

Pore-Network Modeling of Isothermal Drying in Porous Media

In this paper we present numerical results obtained with a pore-network model for the drying of porous media that accounts for various processes at the pore scale. These include mass transfer by advection and diffusion in the gas phase, viscous flow in the liquid and gas phases and capillary effects at the liquid--gas interface. We extend our work by studying the effect of capillarity-induced flow in macroscopic liquid films that form at the pore walls as the liquid--gas interface recedes. A mathematical model that accounts for the effect of films on the drying rates and phase distribution patterns is presented. It is shown that film flow is a major transport mechanism in the drying of porous materials, its effect being dominant when capillarity controls the process, which is the case in typical applications.     

A Two-Dimensional Network Simulator for Two-Phase Flow in Porous Media

We investigate a two-dimensional network simulator that model the dynamics of drainage dominated flow where film flow can be neglected. We present a new method for simulating the temporal evolution of the pressure due to capillary and viscous forces in the displacement process. To model the dynamics, we let the local capillary pressure change as if the menisci move in and out of hour-glass shaped tubes. Furthermore, a method has been developed to allow simultaneous flow of two liquids into one tube. The model is suitable to simulate different time dependencies in two-phase drainage displacements. In this paper, we simulate the temporal evolution of the fluid pressures and analyze the time dependence of the front between the two liquids. The front width was found to be consistent with a scaling relation w t h(t/t_s). The dynamical exponent, , describing the front width evolution as function of time, was estimated to = 1.0. The results are compared to experimental data of Frette and co-workers.     

Effect of Wetting on the Dynamics of Drainage in Porous Media

We examine the consequences of the wettability properties on the dynamics of gravity drainage in porous media. The relation between the wetting properties at the pore scale and the macroscale hydrodynamics is studied. Model porous media consisting of hydrophilic and hydrophobic glass beads or sand with well defined wetting properties, are prepared for this study. Gravity drainage experiments with air displacing water (two-phase flow), are performed for different Bond numbers, and using different techniques such as gamma-ray densitometry, magnetic resonance imaging (MRI) and weight measurements. The dynamics of drainage is found to be different for hydrophilic and hydrophobic porous media in the transition zone (funicular regime). Moreover, for hydrophilic (water-wet) porous media, MRI experiments reveal the importance of drainage through the continuous water film, which leads to an increase of the residual quantity of water in the transition zone with time.     

Effective Diffusivities of Point-Like Molecules in Isotropic Porous Media by Monte Carlo Simulation

Monte Carlo simulations of point-like molecules in random and structured media are used to determine and characterize the effective diffusion coefficients of the molecules in the media. Simulations were carried out in 2D and 3D media. Monte Carlo simulation results in 2D and 3D media are compared with those obtained by analytical techniques. Simulation results indicate that for the structured, isotropic media the effective diffusivities can be characterized according to percolation thresholds in addition to porosity. The effective diffusivities in two isotropic media with the same porosity but different percolation thresholds can differ significantly. The effects of dimensionality on the effective diffusivities can also be significant. It is shown that in general the effective diffusion coefficients obtained from 2D simulation are not a good approximation to those of 3D, especially when the percolation thresholds of the 2D media and the 3D media are very different.     

Foam Mobility in Heterogeneous Porous Media

It is shown experimentally that in situ generation of foam is an effective method for achieving gas mobility control and diverting injected fluid to low permeability strata within heterogeneous porous media. The experimental system is composed of a 0.395 porosity, 5.35 µm² synthetic sandstone and a 0.244 porosity, 0.686 µm² natural sandstone. The cores are arranged in parallel and communicate through common injection and production conditions. Nitrogen is the gas phase and alpha-olefin sulfonate (AOS 1416) in brine is the foamer. Three types of experiments were conducted. First, gas alone was injected into the system after presaturation with the foamer solution. Second, gas and foamer solution were coinjected at an overall gas fraction of 90% into cores presaturated with surfactant. Each core accepted a portion of the injected gas and liquid according to the mobility within the core. Lastly, gas and foamer solution were coinjected into the individual, isolated porous media in order to establish baseline behavior. The results are striking. It is possible to achieve total diversion of gas injection to the low permeability medium in some cases. The results also confirm previous predictions that foamed gas can be more mobile in lower permeability porous media.     

Network Modeling of Non-Darcy Flow Through Porous Media

Darcy's law is inadequate for describing high-velocity gas flow in porous media, which occurs in the near well-bore region of high capacity gas and condensate reservoirs. This study is directed at understanding the non-Darcy flow behavior. A pore-level network model has been developed to describe high velocity flow. The inputs to the model are pore size distributions and network coordination numbers. The outputs are permeability, non-Darcy coefficient, tortuousity and porosity. The additional pressure gradient term is found to be proportional to the square of the velocity in accordance with the Forchheimer's equation. The correlation between the non-Darcy coefficient and other flow properties (the permeability, the porosity and the tortuousity) is found to depend on the morphological parameters being changed. General correlations are derived between these flow properties.     

Modeling Transverse Dispersion and Variable Density Flow in Porous Media

A two-dimensional numerical model is used to study the nonlinear behavior of density gradients on transverse dispersion. Numerical simulations are conducted using d ³ f, a computer code for simulation of density-dependent flow in porous media. Considering a density-stratified horizontal flow in a heterogeneous porous media, a series of simulations is carried out to examine the effect of the density gradient on macro-scale transverse dispersivity. Changing salt concentration significantly affects fluid properties. This physical behavior of the fluid involves a non-linearity in modeling the interaction between salt and fresh water. It is concluded that the large-scale transport properties for high density flow deviate significantly from the tracer case due to the spatial variation of permeability, described by statistical parameters, at the local-scale. Indeed, the presence of vertical flow velocities induced by permeability variations is responsible for the reduction of the mixing zone width in the steady state in the case of a high density gradient. Uncertainties in the model simulations are studied in terms of discretization errors, boundary conditions, and convergence of ensemble averaging. With respect to the results, the gravity number appears to be the controlling parameter for dispersive flux. In addition, the applicability and limitations of the nonlinear model of Hassanizadeh (1990) and Hassanizadeh and Leijnse (1995) (Adv Water Resour 18(4):203–215, 1995) in heterogeneous porous media are investigated. We found that the main cause of the nonlinear behavior of dispersion, which is the interaction between density contrast and vertical velocity, needs to be explicitly accounted for in a macro-scale model.