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.     

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.     

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.     

Numerical Modeling of Momentum Dispersion in Porous Media Based on the Pore Scale Prevalence Hypothesis

Transport in Porous Media - A macroscopic model that accounts for the effect of momentum dispersion on flows in porous media is proposed. The model is based on the pore scale prevalence hypothesis...     

A Resistance Model for Newtonian and Power-Law Non-Newtonian Fluid Transport in Porous Media

A theoretical model based on the force balance between pressure, viscous force, and inertia force is proposed to predict the flow resistance of Newtonian and power-law non-Newtonian fluids through porous packed beds. The present model takes inertia effect into consideration, and the flow regime can be extended from Darcy flow to non-Darcy flow. It is demonstrated that the present model can predict most available experimental data well. The present results are also compared to the Ergun equation and other drag correlations.     

Sound Propagation in Rigid Porous Media: Non-local Macroscopic Effects Versus Pores Scale Regime

This paper investigates how the regime (quasi-static, transient, out of equilibrium) of the phenomena occurring at pores scale determine the nature of the (non)-local effect—in time and/or space—involved in the macroscopic behavior of a porous medium. The study focuses on sound propagation examining—through the homogenization method of periodic media—situations of single porosity, Rayleigh scattering and double porosity. Non-locality effects reveals the loss of a perfect quasi-static equilibrium free of volume loading at the local scale. The non-locality in time is due to phenomena in transient regime at the ERV scale, while non-locality in space is due to the non-homogeneity in space of the macrofields. The generality of the arguments lead to infer that the conclusions about non-locality versus pores scale regime, could be extended to other physical phenomena in heterogeneous media.     

膨胀性非牛顿流体多孔介质流动的模拟分析

在表征体元尺度采用格子Boltzmann方法分析膨胀性非牛顿流体在多孔介质中的流动,基于二阶矩模型在演化方程中引入表征介质阻力的作用力项,求解描述渗流模型的广义Navier-Stokes方程．采用局部法计算形变速率张量,通过循环迭代得到非牛顿粘度和松弛时间．对多孔介质的Poiseuille流动进行分析,通过比较发现结果与孔隙尺度的解析解十分吻合,并且收敛较快,表明方法合理有效．分析了渗透率和幂律指数对速度和压力降的影响,研究结果表明,膨胀性流体的多孔介质流动不符合达西规律,压力降的增加幅度小于渗透率的减小幅度．当无量纲渗透率Da小于10-5时,流道中的速度呈现均匀分布,并且速度分布随着幂律指数的减小趋于平滑．压力降随着幂律指数的增加而增加,Da越大幂律指数对压力降的影响越明显．     

Foam Generation,Propagation and Stability in Porous Medium

Transport in Porous Media - There have been several foam field applications in recent years. Foam treatments targeting gas mobility control in injectors as well as gas blocking in production wells...     

A New Stochastic Bubble Population Model for Foam Flow in Porous Media

Foam has been widely used as a mobility control agent for Improved and Enhanced Oil Recovery IOR/EOR, gas blocking, and acid diversion during matrix stimulation. The prediction of foam performance relies on macroscopic modeling. Traditionally, foam modeling approaches include fractional flow theories and population balance models. However, fractional foam models assume implicitly that foam is incompressible and do not account directly for the evolution of bubble population. The population balance models, instead, rely on the idea that foam mobility depends on bubble density and are more comprehensive. Yet, population balance models did not gain full acceptance thus far, because of their perceived complexity, with parameters that are hard to obtain experimentally. This article presents an improved foam model based on a simpler but realistic foam rheology and stochastic bubble generation ideas. Physical ideas in agreement with pictures emerging from recent foam studies using X-ray computed tomography form the basis for the new model.     

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.     

A Numerical Method for Wave Propagation in Viscoelastic Stratified Porous Media

This paper presents a numerical method, a transmission matrix method, for the wave propagation in viscoelastic stratified saturated porous media. The wave propagation in saturated media, based on Biot theory, is a coupled problem. In this stratified three-dimensional model we do the Laplace transform for the time variable and the Fourier transform for the horizontal space coordinate. The original problem is transformed into ordinary differential equations with six independent unknown variables, which are only the function of the coordinate of depth. Thus, we get a transmission matrix of the wave problem for each layer. In the process of solution we use numerical method to calculate the eigenvalues and the eigenvectors of the transmission matrices. In the first step of the solution process we can obtain the wave field in the transformed space. The fast Fourier transform (FFT) method is used to do the inverse Laplace and the inverse Fourier transforms to get the solution in the time space. The detailed formulae are derived and some numerical examples are given.     

Nonsimilar Solutions for Mixed Convection in Non-Newtonian Fluids Along Horizontal Surfaces in Porous Media

A nonsimilar boundary layer analysis is presented for the problem of mixed convection in power-law type non-Newtonian fluids along horizontal surfaces with variable heat flux distribution. The mixed convection regime is divided into two regions, namely, the forced convection dominated regime and the free convection dominated regime. The two solutions are matched. Numerical results are presented for the details of the velocity and temperature fields. A discussion is provided for the effect of viscosity index on the surface heat transfer rate.     

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.     

Propagation and Evolution of Wave Fronts in Two-Phase Porous Media

An investigation is made into the propagation and evolution of wave fronts in a porous medium which is intended to contain two phases: the porous solid, referred to as the skeleton, and the fluid within the interconnected pores formed by the skeleton. In particular, the microscopic density of each real material is assumed to be unchangeable, while the macroscopic density of each phase may change, associated with the volume fractions. A two-phase porous medium model is concisely introduced based on the work by de Boer. Propagation conditions and amplitude evolution of the discontinuity waves are presented by use of the idea of surfaces of discontinuity, where the wave front is treated as a surface of discontinuity. It is demonstrated that the saturation condition entails certain restrictions between the amplitudes of the longitudinal waves in the solid and fluid phases. Two propagation velocities are attained upon examining the existence of the discontinuity waves. It is found that a completely coupled longitudinal wave and a pure transverse wave are realizable in the two-phase porous medium. The discontinuity strength of the pore-pressure may be determined by the amplitude of the coupled longitudinal wave. In the case of homogeneous weak discontinuities, explicit evolution equations of the amplitudes for two types of discontinuity waves are derived.     

Squirt-Flow in Fluid-Saturated Porous Media: Propagation of Rayleigh Waves

Propagation of attenuated waves is studied in a squirt-flow model of porous solid permeated by two different pore regimes saturated with same viscous fluid. Presence of soft compliant microcracks embedded in the grains of stiff porous rock defines the double-porosity formation. Microcracks and pores respond differently to the compressional effect of a propagating wave, which induces the squirt-flow from microcracks to pores. Elastodynamics of constituent particles in porous aggregate is represented through a single-porosity formulation, which involves the frequency-dependent complex moduli. This formulation is deduced as a special case of double-porosity formation allowing the wave-induced flow of pore-fluid. This squirt-flow model of porous solid supports the attenuated propagation of two compressional waves and one shear wave. Superposition of these body waves, subject to stress-free surface, defines the propagation of Rayleigh wave. This wave is governed by a complex irrational dispersion equation, which is solved numerically after rationalising into an algebraic equation. For existence of Rayleigh wave, a complex solution of the dispersion equation should represent a leaky wave, which decays for propagation along any direction in the semi-infinite medium. A numerical example is solved to analyse the effects of squirt-flow on phase velocity, attenuation and polarisation of the Rayleigh waves, for different combinations of parameters. Numerical results suggest the existence of an additional (second) Rayleigh wave in the squirt-flow model of dissipative porous solids.     

The Peculiarities of Linear Wave Propagation in Double Porous Media

The features of propagation of longitudinal and transverse waves (LW and TW) in fractured porous medium (FPM) saturated with liquid are investigated by methods of multiphase mechanics. The mathematical model of FPM accounting for inequality of velocities and pressures of liquid in pores and fractures, liquid mass exchange and nonstationary interaction forces is developed. Processes of monochromatic wave propagation are studied. The dispersion relation is obtained and the effect of model parameters on wave propagation is analysed. It is established that one transverse and three longitudinal waves propagate in FPM saturated with liquid. The fastest LW is a deformational wave and the two others are filtrational. Filtrational waves attenuate much stronger than deformational and transverse waves. Distinction of velocities and pressures in liquid in various pore systems provides an explanation for the existence of the two filtrational waves in porous medium with two different characteristic sizes of pores.     

Foam Generation Hysteresis in Porous Media: Experiments and New Insights

Foam application in subsurface processes including environmental remediation, geological carbon-sequestration, and gas-injection enhanced oil recovery (EOR) has the potential to enhance contamination remediation, secure \(\hbox {CO}_{2}\) storage, and improve oil recovery, respectively. Nanoparticles are a promising alternative to surfactants in creating foam in harsh environments. We conducted \(\hbox {CO}_{2}\)-in-brine foam generation experiments in Boise sandstones with surface-treated silica nanoparticle in high-salinity conditions. All the experiments were conducted at the fixed \(\hbox {CO}_{2}\) volume fraction and fixed flow rate which changed in steps. The steady-state foam apparent viscosity was measured as a function of injection velocity. The foam flowing through the cores showed higher apparent viscosity as the flow rate increased from low to medium and high velocities. At very high velocities, once foam bubbles were finely textured, the foam apparent viscosity was governed by foam rheology rather than foam creation. A noticeable hysteresis occurred when the flow velocity was initially increased and then decreased, implying multiple (coarse and strong) foam states at the same superficial velocity. A normalized generation function was combined with CMG-STARS foam model to cover full spectrum of foam behavior in the experiments. The new model successfully captures foam generation and hysteresis trends in presented experiments in this study and data from the literature. The results indicate once foam is generated in porous media, it is possible to maintain strong foam at low injection rates. This makes foam more feasible in field applications where foam generation is limited by high injection rates that may only exist near the injection well.     

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...     

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.