Dispersion Measurements in Highly Heterogeneous Laboratory Scale Porous Media

Laboratory experiments and numerical simulations investigate conservative contaminant transport in a heterogeneous porous medium. The laboratory experiments were performed in cylindrical columns 1m long and 3.5cm inside diameter filled with spherical glass beads. Concentration breakthrough curves are measured at a scale much finer than the size of the heterogeneity. Numerical simulations are based on a random walk in a known constant velocity field. The heterogeneity is a distinct, discontinuous change in the local permeability field. Fluid flow is miscible, flowing in a saturated porous medium. Previous work has shown this to be a very poorly understood phenomenon. The measurements reported here help to better understand how dispersion evolves through and past a heterogeneity.     

A Pore Scale Study of Non-Newtonian Effect on Foam Propagation in Porous Media

CO₂ injection is one of the most promising techniques to enhance oil recovery. However, an unfavorable mobility ratio, reservoir heterogeneity and gravity segregation can reduce the macroscopic sweep efficiency. In situ foaming of injected CO₂ is the method that has the most potential for improving sweep efficiency based on controlling CO₂ mobility. This study investigates the foaming behavior of N,N,N′-trimethyl-N′-tallow-1,3-diaminopropane (DTTM) surfactant with CO₂ in a transparent porous microflow model with natural rock pore structures. It focuses on the effect of the salinity induced non-Newtonian behavior of DTTM solution on foam propagation. The performance of foams stabilized by 0.5 wt% DTTM solution over the viscosity range from 0.71 (at 5 wt% NaCl) to 41 cp (at 20 wt% NaCl) was compared with conventional polymer-enhanced foams whose liquid phase contained a commonly used foaming surfactant, C15–18 Internal Olefin Sulfonate (C15–18 IOS) and a hydrolyzed polyacrylamide. Such comparisons have also provided insight into the respective impacts of liquid phase viscosification by worm-like surfactant micelles and polymer on foam texture associated with its rheological characteristics. It was found that at low aqueous phase viscosity (injection liquid viscosity of 0.71 cp) the maximum achievable viscosity of DDTM foam was around 1000 cp, which was 80 times IOS stabilized foam. The interfacial tension of DTTM was higher than that of IOS, resulting coarser foam texture and higher individual lamella resistance. An increase in DTTM solution viscosity by a factor of 33 decreased foam generation and viscosity for gas injection. This was not observed for the simultaneous injection of gas and DTTM solution. Overall, the effect of liquid phase viscosity on transient foam behavior during gas injection is similar for both DTTM and IOS regardless of the difference in the nature of viscosifying agents (WLM vs 3330 s polymer). An increase in gas injection pressure without liquid injection delayed foam propagation and reduced the magnitude of foam viscosity. The results from this study indicated that DTTM surfactant is an important alternative to commercially available polymers that have been used to enhance foam performance in porous media. This particular surfactant type also overcomes several disadvantages of polymers such as limited temperature and salinity tolerance, shear degradation, and filtering in low permeability formations.

     

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.     

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.     

Dispersion and Chemical Reactions in Porous Media

The effect of chemical reactions on the process of admixture transport by a flow through a porous medium is considered. On the basis of a number of examples it is shown that the dispersion coefficient depends on the chemical reaction rate constant.     

Effect of Vibrations on Pore Fluid Distribution in Porous Media

Understanding the role of shuttle vibrations in pore fluid distribution is an essential task in the exploration of plant growth in root modules aboard space flights. Results from experimental investigations are reported in this paper on the distribution of immiscible fluid phases in glass beads under vibrations. Hexadecane, a petroleum compound immiscible with and lighter than water, was used in the experiments. The higher freezing point of Hexadecane (18 °C) allowed the solidification of the entrapped blobs in the presence of water in porous media, so that their size distribution can be obtained. van Genuchten function, commonly used to express moisture retention curves, is found to be an adequate fit for blob size distribution at residual saturation. The effect of vibrations on the fate (mobilization, stranding, or breakup) of a solitary ganglion in porous media was studied using a network model. A mobility criterion considering viscous, gravity, and capillary forces was developed to determine the fate of a solitary ganglion in a porous medium. It is concluded that the effect of vibrations is to increase the likelihood of breakup and mobilization of blobs entrapped in porous media at residual saturation. The pore fluid distributions after vibrations are less uniform than those before vibrations.     

Pore Network Modeling of the Effects of Viscosity Ratio and Pressure Gradient on Steady-State Incompressible Two-Phase Flow in Porous Media

We perform steady-state simulations with a dynamic pore network model, corresponding to a large span in viscosity ratios and capillary numbers. From these simulations, dimensionless steady-state time-averaged quantities such as relative permeabilities, residual saturations, mobility ratios and fractional flows are computed. These quantities are found to depend on three dimensionless variables, the wetting fluid saturation, the viscosity ratio and a dimensionless pressure gradient. Relative permeabilities and residual saturations show many of the same qualitative features observed in other experimental and modeling studies. The relative permeabilities do not approach straight lines at high capillary numbers for viscosity ratios different from 1. Our conclusion is that this is because the fluids are not in the highly miscible near-critical region. Instead they have a viscosity disparity and intermix rather than forming decoupled, similar flow channels. Ratios of average mobility to their high capillary number limit values are also considered. Roughly, these vary between 0 and 1, although values larger than 1 are also observed. For a given saturation, the mobilities are not always monotonically increasing with the pressure gradient. While increasing the pressure gradient mobilizes more fluid and activates more flow paths, when the mobilized fluid is more viscous, a reduction in average mobility may occur.

     

Pore Scale Modeling of Rate Effects in Imbibition

We use pore scale network modeling to study the effects of flow rate and contact angle on imbibition relative permeabilities. The model accounts for flow in wetting layers that occupy roughness or crevices in the pore space. Viscous forces are accounted for by solving for the wetting phase pressure and assuming a fixed conductance in wetting layers. Three-dimensional simulations model granular media, whereas two-dimensional runs represent fracture flow.We identify five generic types of displacement pattern as we vary capillary number, contact angle, and initial wetting phase saturation: flat frontal advance, dendritic frontal advance, bond percolation, compact cluster growth, and ramified cluster growth. Using phase diagrams we quantify the range of physical properties under which each regime is observed. The work explains apparently inconsistent experimental measurements of relative permeability in granular media and fractures.     

Dispersion in Heterogeneous Porous Media: One-Equation Non-equilibrium Model

In this paper, the method of large-scale averaging is used to develop two different one-equation models describing dispersion in heterogeneous porous media. The first model represents the case of large-scale mass equilibrium, while the second represents the asymptotic behavior of a two-equation model obtained by large-scale averaging. It is shown that a one-equation, non-equilibrium model can be developed even when the intrinsic large-scale averaged concentrations for each region are not equal. The solution of this non-equilibrium model is equivalent to the asymptotic behavior of the two-equation model.     

Review on Scale Dependent Characterization of the Microstructure of Porous Media

The paper discusses local porosity theory and its relation with other geometric characterization methods for porous media such as correlation functions and contact distributions. Special emphasis is placed on the charcterization of geometric observables through Hadwigers theorem in stochastic geometry. The four basic Minkowski functionals are introduced into local porosity theory, and for the first time a relationship is established between the Euler characteristic and the local percolation probabilities. Local porosity distributions and local percolation probabilities provide a scale dependent characterization of the microstructure of porous media that can be used in an effective medium approach to predict transport.     

Dispersion Stability and Transport of Nanohybrids through Porous Media

Single-walled carbon nanotube-silica nanohybrid particles are a very promising material that could be used for enhanced oil recovery because of their interfacial activity. To demonstrate the basic principle, aqueous nanohybrid particle dispersions were evaluated by looking at the effect of pH, surfactant, and polymer. The results showed that pH did not have significant effect on the dispersion stability of nanohybrid particles. Although surfactant could improve the dispersion stability, it reduced the interfacial activity of the nanohybrid particles, causing them to stay in the aqueous phase. The best nanohybrid particle dispersion stability was found upon polymer addition, where the dispersions were stable for more than a week even at low polymer concentration (50?ppm). One-dimensional sand-pack studies were performed to evaluate the flow of the nanohybrid particles through porous media. The results showed that most of the nanohybrid particles (>99%) could pass through a column packed with glass beads while a measurable fraction of the particles was retained in the column packed with crushed Berea. When the columns contained a residual saturation of decane, additional nanohybrid particles were retained at the oil/water interface in both glass beads and crushed Berea sand media. The sand pack studies showed that not only can the nanohybrid particles flow through porous media but also about half of the particles injected will go the O/W interface when the porous medium contains a residual saturation of hydrocarbon, where they could be used to support a catalytic conversion of components of the oil in reservoirs.     

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.     

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.     

Numerical Algorithms for Network Modeling of Yield Stress and other Non-Newtonian Fluids in Porous Media

Many applications involve the flow of non-Newtonian fluids in porous, subsurface media including polymer flooding in enhanced oil recovery, proppant suspension in hydraulic fracturing, and the recovery of heavy oils. Network modeling of these flows has become the popular pore-scale approach for understanding first-principles flow behavior, but strong nonlinearities have prevented larger-scale modeling and more time-dependent simulations. We investigate numerical approaches to solving these nonlinear problems and show that the method of fixed-point iteration may diverge for shear-thinning fluids unless sufficient relaxation is used. It is also found that the optimal relaxation factor is exactly equal to the shear-thinning index for power-law fluids. When the optimal relaxation factor is employed it slightly outperforms Newton??s method for power-law fluids. Newton-Raphson is a more efficient choice (than the commonly used fixed-point iteration) for solving the systems of equations associated with a yield stress. It is shown that iterative improvement of the guess values can improve convergence and speed of the solution. We also develop a new Newton algorithm (Variable Jacobian Method) for yield-stress flow which is orders of magnitude faster than either fixed-point iteration or the traditional Newton??s method. Recent publications have suggested that minimum-path search algorithms for determining the threshold pressure gradient (e.g., invasion percolation with memory) greatly underestimate the true threshold gradient when compared to numerical solution of the flow equations. We compare the two approaches and reach the conclusion that this is incorrect; the threshold gradient obtained numerically is exactly the same as that found through a search of the minimum path of throat mobilization pressure drops. This fact can be proven mathematically; mass conservation is only preserved if the true threshold gradient is equal to that found by search algorithms.