Darcian Dynamics: A New Approach to the Mobilization of a Trapped Phase in Porous Media

We develop a new approach, which we term Darcian Dynamics, to simulate two-phase (liquid-gas) flow in porous media, when the gas phase is disconnected in the form of ganglia. The method is based on the assumption of homogeneous fluid flow for the liquid, although it does allow for heterogeneous capillary thresholds due to the pore microstructure. Using techniques from potential theory, the hydrodynamic interaction between liquid and gas is expressed through an integral representation over the ganglia interfaces. We use a numerical method to solve the resulting integral equation, and explore conditions for the onset of ganglia mobilization as well as for subsequent events, such as break-up, coalescence and stranding. The interaction between the ganglia and the flowing phase is influenced by the capillary and gravity (Bond) numbers, and by geometric factors, such as size, orientation, and ganglia density. The latter effect depends on the hydrodynamic interaction in addition to the intuitively expected crowding effect.     

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.     

Surfactant Concentration and End Effects on Foam Flow in Porous Media

Foaming injected gas is a useful and promising technique for achieving mobility control in porous media. Typically, such foams are aqueous. In the presence of foam, gas and liquid flow behavior is determined by bubble size or foam texture. The thin-liquid films that separate foam into bubbles must be relatively stable for a foam to be finely textured and thereby be effective as a displacing or blocking agent. Film stability is a strong function of surfactant concentration and type. This work studies foam flow behavior at a variety of surfactant concentrations using experiments and a numerical model. Thus, the foam behavior examined spans from strong to weak.Specifically, a suite of foam displacements over a range of surfactant concentrations in a roughly 7m², one-dimensional sandpack are monitored using X-ray computed tomography (CT). Sequential pressure taps are employed to measure flow resistance. Nitrogen is the gas and an alpha olefin sulfonate (AOS 1416) in brine is the foamer. Surfactant concentrations studied vary from 0.005 to 1wt%. Because foam mobility depends strongly upon its texture, a bubble population balance model is both useful and necessary to describe the experimental results thoroughly and self consistently. Excellent agreement is found between experiment and theory.     

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.     

Monte Carlo Simulations of Observation Time-Dependent Self-Diffusion in Porous Media Models

Observation time-dependent self-diffusion coefficients can be used to obtain microstructural information of porous media. This paper presents two different kinds of Monte Carlo simulations of the self diffusion process of fluids like water in porous systems, a lattice-free method and a lattice-based method. The results for simple porous media model geometries agree well with each other and with published analytical as well as semi-analytical equations. The use of these equations, which are important for the interpretation of Pulsed Field Gradient-Nuclear Magnetic Resonance (PFG-NMR) time-dependent diffusion data with respect to properties of porous media, is discussed.     

Thermodynamic Automaton Simulations of Fluid Flow and Diffusion in Porous Media

A thermodynamic automaton model of fluid flow in porous media is presented. The model is a nonrelativistic version of a Lorentz invariant lattice gas model constructed by Udey et al. (1998). In the previous model it was shown that the energy momentum tensor and the relativistic Boltzman equation can be rigorously derived from the collision and propagation rules. In the present paper we demonstrate that this nonrelativistic model can be used to accurately simulate well known results involving single phase flow and diffusion in porous media. The simulation results show that (1) one-phase flow simulations in porous media are consistent with Darcy's law; (2) the apparent diffusion coefficient decreases with a decrease in permeability; (3) small scale heterogeneity does not affect diffusion significantly in the cases considered.     

Network Model of Flow,Transport and Biofilm Effects in Porous Media

In this paper, we develop a network model to determine porosity and permeability changes in a porous medium as a result of changes in the amount of biomass. The biomass is in the form of biofilms. Biofilms form when certain types of bacteria reproduce, bond to surfaces, and produce extracellular polymer (EPS) filaments that link together the bacteria. The pore spaces are modeled as a system of interconnected pipes in two and three dimensions. The radii of the pipes are given by a lognormal probability distribution. Volumetric flow rates through each of the pipes, and through the medium, are determined by solving a linear system of equations, with a symmetric and positive definite matrix. Transport through the medium is modeled by upwind, explicit finite difference approximations in the individual pipes. Methods for handling the boundary conditions between pipes and for visualizing the results of numerical simulations are developed. Increases in biomass, as a result of transport and reaction, decrease the pipe radii, which decreases the permeability of the medium. Relationships between biomass accumulation and permeability and porosity reduction are presented.     

Role of Thermal Diffusion on Heat and Mass Transfer in Porous Media

In this paper, thermal diffusion phenomena in a porous cavity are investigated. The Brinkman model, coupled with the energy and the mass balance equations was solved numerically using a finite element techniques. A two-component system was included in the model. Different models were investigated to demonstrate the importance of the Soret effect with the presence of gravity vector. We do not take into consideration the pressure effect in the thermal diffusion. Even with such simplification to the problem, results reveal that the thermal diffusion is important and drives a strong convection. A series of convection cells are observed and steady-state solutions are obtained. Asymmetric solutions are obtained for various cases of dual-porosity porous media. Variations in the gravity vector indicated that the convection patterns, as well as the role of Soret coefficient, are profoundly impacted. Finally, the importance of including thermal diffusion in petroleum reservoir simulation is discussed.     

Pore-Scale Analysis of NAPL Blob Dissolution and Mobilization in Porous Media

A pore-scale analysis of nonaqueous phase liquid (NAPL) blob dissolution and mobilization in porous media was presented. Dissolution kinetics of residual NAPLs in an otherwise water-saturated porous medium was investigated by conducting micromodel experiments. Changes in residual NAPL volume were measured from recorded video images to calculate the mass transfer coefficient, K and the lumped mass transfer rate coefficient, k. The morphological characteristics of the blobs such as specific and intrinsic area were found to be independent of water flow rate except at NAPL saturations below 2%. Dissolution process was also investigated by separating the mass transfer into zones of mobile and immobile water. The fractions of total residual NAPL perimeters in contact with mobile water and immobile water were measured and their relationship to the mass transfer coefficient was discussed. In general, residual NAPLs are removed by dissolution and mobilization. Although these two mechanisms were studied individually by others, their simultaneous occurrence was not considered. Therefore, in this study, mobilization of dissolving NAPL blobs was investigated by an analysis of the forces acting on a trapped NAPL blob. A dimensional analysis was performed to quantify the residual blob mobilization in terms of dimensionless Capillary number (Ca _I). If Ca _I is equal to or greater than the trapping number defined as , then blob mobilization is expected.     

Automaton Simulations of Dispersion in Porous Media

A thermodynamic lattice gas (automaton) model is used to simulate dispersion in porous media. Simulations are constructed at two distinctly different scales, the pore scale at which capillary models are constructed and large scale or Darcy scale at which probabilistic collision rules are introduced. Both models allow for macroscopic (pore scale) phase separation. The pore scale models clearly show the effect of pore structure on dispersion. The large scale (mega scale) simulations indicate that when the pressure difference between the displacing phase and displaced phase is properly chosen (representing the average pressure gradient between the phases). The simulation results are consistent with both theoretical predictions and experimental observations.     

Onset of Buoyancy-driven Instability in Porous Media Melted from Below

When a nonhomogeneous solid is melting from below, convection may be induced in a thermally–unstable melt layer. In this study, the onset of buoyancy-driven convection during time-dependent melting is investigated by using similarly transformed disturbance equations. The critical Darcy–Rayleigh numbers based on the melt-layer thickness, Ra _H,c, are found numerically for various conditions. For small superheats, the present predictions show that Ra _H,c is located between 27.1 and 4π ² and it approaches the well-known results of the original Horton–Rogers–Lapwood problem. However, for high superheats, it is dependent on the phase change rate λ and the relation of Ra _H,c λ = 25.89 is shown.     

A New Model for Immiscible Displacement in Porous Media

Immiscible displacement is regarded as the superposition of forward flows of both water and oil, due to injection of water into the medium, and of additional forward flow of water coupled with reverse flow of oil, caused by the existence of capillary pressure gradients. The model has been evaluated numerically for the prediction of the evolution of saturation profiles in waterfloods covering a wide range of water injection rates. In agreement with experimentation, saturation profiles ranging from a completely flat shape to piston-shape, depending on the injection rate, have been obtained. Also in agreement with experimentation, numerical evaluation of the model for the case of a closed system with an initial step-function saturation profile has predicted a gradual spreading of the piston front into S-shaped profiles with an increasing variance. The final profile corresponds to uniform saturation everywhere in the medium.     

Modeling Nanoparticle Transport in Porous Media in the Presence of a Foam

Transport in Porous Media - Nano-remediation is a promising in situ remediation technology. It consists in injecting reactive nanoparticles (NPs) into the subsurface for the displacement or the...     

A Multiscale Network Model for Simulating Moisture Transfer Properties of Porous Media

A multiscale network model is presented to model unsaturated moisture transfer in hygroscopic capillary-porous materials showing a broad pore-size distribution. Both capillary effects and water sorption phenomena, water vapour and liquid water transfer are considered. The multiscale approach is based on the concept of examining the porous space at different levels of magnification. The conservation of the water vapour permeability of dry material is used as scaling criterion to link the different pore scales. A macroscopic permeability is deduced from the permeabilities calculated at the different levels of magnification. Each level of magnification is modelled using an isotropic nonplanar 2D cross-squared network. The multiscale network simulates the enhancement of water vapour permeability due to capillary condensation, the hysteresis phenomenon between wetting and drying, and the steep increase of moisture permeability at the critical moisture saturation level. The calculated network permeabilities are compared with experimental data for calcium silicate and ceramic brick and a good agreement is observed.     

Modeling and Experimental Validation of Rheological Transition During Foam Flow in Porous Media

Transport in Porous Media - Flow of nitrogen foam stabilized by alpha olefin sulfonate (C14-16 AOS) was studied in a natural sandstone porous media using X-ray Computed Tomography. Foam was...     

Diffusion of Randon in Porous Media Saturated with Gels and Emulsions

The effective diffusion coefficient of radon was determined in polymer/silicate gels and clay suspension used as sealing materials in environmental protection. On the basis of the experimental findings, it was concluded that both materials drastically decrease the convective mass transport in porous media. Simultaneously, the effective diffusion coefficient was reduced significantly. Thus, the radon flux might be decreased by 5 to 6 orders of magnitude in porous systems originally having gas or low water saturation by injection of gel-forming materials or placement of clay suspensions. At high water saturation, however, the diffusion transport of radon can be slightly restricted in consolidated and unconsolidated porous media. The laboratory studies may firmly allow us to conclude that hydrogels and clay suspensions are prospective candidates in an integrated environmental technology to be used for restriction of radon migration in subsurface regions.     

The Use of Fick's Law for Modeling Trace Gas Diffusion in Porous Media

Two models for combined gas-phase diffusion and advection in porous media, the advective-diffusive model (ADM) and the dusty-gas model (DGM), are commonly used. The ADM is based on a simple linear addition of advection calculated by Darcy's law and ordinary diffusion using Fick's law with a porosity–tortuosity–gas saturation multiplier to account for the porous medium. The DGM applies the kinetic theory of gases to the gaseous components and the porous media (or dust) to develop an approach for combined transport due to diffusion and advection that includes porous medium effect. The ADM and Fick's law are considered to be generally inferior for gas diffusion in porous media, and the more mechanistic DGM is preferred. Under trace gas diffusion conditions, Fick's law overpredicts the gas diffusion flux compared to the DGM. The difference between the two models increases as the permeability decreases. In addition, the difference decreases as the pressure increases. At atmospheric pressure, the differences are minor (<10%) for permeabilities down to about 10^–13 m². However, for lower permeabilities, the differences are significant and can approach two orders of magnitude at a permeability of 10^–18 m². In contrast, at a pressure of 100 atm, the maximum difference for a permeability of 10^–18 m² is only about a factor of 2. A molecule–wall tortuosity coefficient based on the DGM is proposed for trace gas diffusion using Fick's law. Comparison of the Knudsen diffusion fluxes has also been conducted. For trace gases heavier than the bulk gas, the ADM mass flux is higher than the DGM. Conversely, for trace gases lighter than the bulk gas, the ADM mass flux is lower than the DGM. Similar to the ordinary diffusion variation, the differences increase as the permeability decreases, and get smaller as the pressure increases. At atmospheric pressure, the differences are small for higher permeabilities (>10^–13 m²) but may increase to about 2.7 for He at lower permeabilities of about 10^–18 m². A modified Klinkenberg factor is suggested to account for differences in the models.     

Nonlinear Pressure Diffusion in Flow of Compressible Liquids Through Porous Media

In the flow of liquids through porous media, nonlinear effects arise from the dependence of the fluid density, porosity, and permeability on pore pressure, which are commonly approximated by simple exponential functions. The resulting flow equation contains a squared gradient term and an exponential dependence of the hydraulic diffusivity on pressure. In the limiting case where the porosity and permeability moduli are comparable, the diffusivity is constant, and the squared gradient term can be removed by introducing a new variable y, depending exponentially on pressure. The published transformations that have been used for this purpose are shown to be special cases of the Cole–Hopf transformation, differing in the choice of integration constants. Application of Laplace transformation to the linear diffusion equation satisfied by y is considered, with particular reference to the effects of the transformation on the boundary conditions. The minimum fluid compressibilities at which nonlinear effects become significant are determined for steady flow between parallel planes and cylinders at constant pressure. Calculations show that the liquid densities obtained from the simple compressibility equation of state agree to within 1% with those obtained from the highly accurate Wagner-Pru?  equation of state at pressures to 20 MPa and temperatures approaching 600 K, suggesting possible applications to some geothermal systems.     

Onset of Buoyancy-Driven Convection in Porous Media Saturated with Cold Water Cooled from Above

When porous media saturated with initially stagnant cold water around the density maximum temperature are cooled from above, convection may be induced in an unstable lower layer. In this study, the onset of buoyancy-driven convection during time-dependent cooling is investigated using the propagation theory, which transforms disturbance equations similarly, and also considering the density inversion effect. The critical Darcy–Rayleigh number Ra _D,c is found as a function of the dimensionless density maximum temperature θ _max. For Ra _D > Ra _D,c the dimensionless critical time τ _c to mark the onset of instability is presented as a function of Ra _D and θ _max. These critical conditions are compared with previous theoretical results.