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.     

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.     

Effects of Precursor Wetting Films in Immiscible Displacement Through Porous Media

A computer-aided simulator of immiscible displacement in strongly water-wet consolidated porous media that takes into account the effects of the wetting films is developed. The porous medium is modeled as a three-dimensional network of randomly sized unit cells of the constricted-tube type. Precursor wetting films are assumed to advance through the microroughness of the pore walls. Two types of pore wall microroughness are considered. In the first type of microroughness, the film advances quickly, driven by capillary pressure. In the second type, the meniscus moves relatively slowly, driven by local bulk pressure differences. In the latter case, the wetting film often forms a collar that squeezes the thread of oil causing oil disconnection. Each pore is assumed to have either one of the aforementioned microroughness types, or both. The type of microroughness in each pore is assigned randomly. The simulator is used to predict the residual oil saturation as a function of the pertinent parameters (capillary number, viscosity ratio, fraction of pores with each type of wall microroughness). These results are compared with those obtained in the absence of wetting films. It is found that wetting films cause substantial increase of the residual oil saturation. Furthermore, the action of the wetting films causes an increase of the mean volume of the residual oil ganglia.     

An Analysis of the Movement of Wetting and Nonwetting Fluids in Homogeneous Porous Media

The movement of wetting and nonwetting fluid flow in columns packed with glass beads is used to understand the more complicated flows in homogeneous porous media. The motion of two immiscible liquids (oil and water) is observed with different surfactants. Through dimensional analyses, fluid velocity is well correlated with interfacial tension and less dependent on particle size. In water–oil (W/O) experiments, finger pattern flows are observed if water is the displacing fluid that flows in an oil-filled porous media, whereas oil ganglia tend to form if oil is the displacing fluid in the water-wetted porous media. The results are well described by a simple model based on an earlier theory of flow in a tube.     

Stokesian Dynamics Simulation of Suspension Flow in Porous Media

Transport in Porous Media - Suspension flow through porous medium was studied using the Stokesian dynamics simulation method. Stokesian dynamics is an efficient tool to carry out numerical...     

Pore-Level Observation of Free Gravity Drainage of Oil in Fractured Porous Media

This work presents results from two sets of experiments conducted to study, in pore level, the role of fracture aperture and tilt angle on the stability of liquid bridges and the shape of a front during free gravity drainage process. Glass micromodels of two different aperture sizes were used to monitor the mechanism of gravity drainage of air?Ccrude oil system, rotating around a bottom corner to create different tilting angles. Oil content within the matrix blocks was determined as a function of time using a series of images obtained during the experiments, from which net drainage rate from the upper and lower matrix blocks is calculated. Liquid bridges are more frequent but less stable at early time of drainage. The liquid bridges, which have widths as thin as 50 ??m, can resist instability to maintain continuity. Liquid bridges formed in stacks with higher tilt angles are more stable, enhancing oil drainage from the upper matrix block and causing higher recoveries. Quantitative analysis of the results shows that a wider fracture aperture increases the oil production rate, but reduces the ultimate recovery. Furthermore, stacks with higher tilt angles present larger ultimate recoveries and smaller production rates. The front geometry in the lower block deviates from linearity due to formation of liquid bridges in the middle fracture. The results of this work can be helpful to better understand the interaction between fractures and matrix blocks.     

Wave Dynamics of Saturated Porous Media and Evolutionary Equations

Nonlinear wave dynamics of an elastically deformed saturated porous media is investigated following the Biot approach. Mathematical models under research are the Biot model and its generalization by consideration of viscous stresses inside liquids. Using two-scales and linear WKB methods, the classical Biot system is transformed to a first-order wave equation. To construct the solution of the other system, an asymptotic modified two-scales method is developed. Initial system of equations is transformed to a nonlinear generalized Korteweg–de Vries–Burgers equation for quick elastic wave. Distinctions of wave propagation in the context of the Biot model and its generalization are shown.     

Accurate and Efficient Implementation of Pore-Morphology-based Drainage Modeling in Two-dimensional Porous Media

We have developed an efficient and accurate numerical implementation for pore-morphological modeling of drainage in two-dimensional, totally wetting porous media. The new numerical method uses level sets to describe the fluid distribution and polygons that can be defined with subgrid scale accuracy for the pore boundaries, while a previously developed approach represents the phases by pixels arranged on a square lattice. We analyze and compare the previous and new method. For both approaches, the simulated fluid saturations are first-order accurate. For the level-set approach, the simulated interfacial lengths converge to the real values, while the pixel approach yields biased results. The level-set method is orders of magnitudes faster than the pixel method.     

A New Mathematical Model for Force Gravity Drainage in Fractured Porous Media

In force gas/oil gravity drainage process in fractured porous media, gas is flowing in both matrix and fractures leading to produce a finite gas pressure gradient. Consequently, viscous force plays an important role for displacing matrix oil toward fractures in addition to gravity force that is required to be modeled appropriately. A new analytical model for estimation of steady state oil saturation distribution with assumption of fixed gas pressure gradient throughout the matrix is presented. Moreover, based on some results of this analytical model a different numerical formulation is developed to predict the performance of oil production process. Comparison of the results obtained from this numerical model with the results of a conventional simulator demonstrates that the newly developed model can be applied with satisfactory accuracy. Numerical simulations show that the viscous displacement in fractured porous media can reduce the capillary threshold height, and thus it suggests the force gravity drainage as a favorable production mechanism when the matrix length is close to the threshold height.     

Mechanistic Model of Steady-State Two-Phase Flow in Porous Media Based on Ganglion Dynamics

Recent experimental work has shown that the pore-scale flow mechanism during steady-state two-phase flow in porous media is ganglion dynamics (GD) over a broad and practically significant range of the system parameters. This observation suggests that our conception and theoretical treatment of fractional flow in porous media need careful reconsideration. Here is proposed a mechanistic model of steady-state two-phase flow in those cases where the dominant flow regime is ganglion dynamics. The approach is based on the ganglion population balance equations in combination with a microflow network simulator. The fundamental information on the cooperative flow behavior of the two fluids at the scale of a few hundred pores is expressed through the system factors, which are functions of the system parameters and are calculated using the simulator. These system factors are utilized by the population balance equations to predict the macroscopic behavior of the process. The dependence of the conventional relative permeability coefficients not only on the wetting fluid saturation S_wbut also on the capillary number, Ca, the viscosity ratio the wettability (⁰ a, ⁰ r), the coalescence factor, Co, as well as the porous medium geometry and topology is explained and predicted on a mechanistic basis. Sample calculations have been performed for steady-state fully developed (SSFD) and steady-state nonfully developed (SSnonFD) flow conditions. The number distributions of the moving and the stranded ganglia, the mean ganglion size, the fraction of the nonwetting fluid in the form of mobile ganglia, the ratio of the conventional relative permeability coefficients and the fractional flows are studied as functions of the system parameters and are correlated with the flow phenomena at pore level and the system factors.     

Dynamics of Viscous Entrapped Saturated Zones in Partially Wetted Porous Media

As a typical multiphase fluid flow process, drainage in porous media is of fundamental interest both in nature and in industrial applications. During drainage processes in unsaturated soils and porous media in general, saturated regions, or clusters, in which a liquid phase fully occupies the pore space between solid grains, affect the relative permeability and effective stress of the system. Here, we experimentally study drainage processes in unsaturated granular media as a model porous system. The distribution of saturated clusters is analysed by optical imaging under different drainage conditions, with pore-scale information from Voronoi and Delaunay tessellation used to characterise the topology of saturated cluster distributions. By employing statistical analyses, we describe the observed spatial and temporal evolution of multiphase flow and fluid entrapment in granular media. Results indicate that the distributions of both the crystallised cell size and pore size are positively correlated to the spatial and temporal distribution of saturated cluster sizes. The saturated cluster size is found to follow a lognormal distribution, in which the generalised Bond number (\( Bo^{*} \)) correlates negatively to the scale parameter (μ) and positively to the shape parameter (σ). With further consideration of the total surface energy obtained based on liquid–air interfaces, we were able to include additional grain-scale information in the constitutive modelling of unsaturated soils using both the degree of saturation and generalised Bond number. These findings can be used to connect pore-scale behaviour with overall hydro-mechanical characteristics in granular systems.     

An Experimental Study of the Factors Influencing Apparent Wetting Phase Residual Saturation in Dynamically Drained Porous Media

When a porous medium drains, the wetting phase saturation in the medium [the fractional volume of wetting phase (often water) in the pores] is typically observed to approach a minimum value, referred to as the wetting phase residual saturation ( $S_{\text {wr}}$ ). While many simulators of unsaturated and multiphase flow assume $S_{\text {wr}}$ to be a single value for a given fluid and porous medium system, there is considerable evidence that $S_{\text {wr}}$ may be a dynamic property, with its value a function of the rate of saturation change. The objective of this work was to study this relationship, with emphasis on understanding the effect of drainage rate on the apparent residual ( $S_{\text {wr-apparent}}$ ) at moderate capillary pressures and over short time frames. Dynamic unsaturated drainage experiments were conducted to explore this relationship. A total of 61 experiments were conducted, involving 363 individual secondary drainages. Experiments covered a total of four different porous media, and three wetting phases. Results indicate that systems that exhibit greater resistance to flow (lower mean grain size and higher kinematic viscosity) also exhibit a greater sensitivity in $S_{\text {wr-apparent}}$ to drainage rate. An imaging experiment was conducted to examine whether the source of the relationship could be observed through continuum-scale imaging. Results of the imaging were able to rule out membrane desaturation artifacts as a cause of the relationship, but were unable to identify significant differences in fluid configurations between fast and slow drainages that might lead to the observed effects.     

Radiative Effect on Natural Convection Flows in Porous Media

A regular two-parameter perturbation analysis based upon the boundary layer approximation is presented here to study the radiative effects of both first- and second-order resistances due to a solid matrix on the natural convection flows in porous media. Four different flows have been studied, those adjacent to an isothermal surface, a uniform heat flux surface, a plane plume and the flow generated from a horizontal line energy source on a vertical adiabatic surface. The first-order perturbation quantities are presented for all these flows. Numerical results for the four conditions with various radiation parameters are tabulated.     

Dynamics of Flow Through Porous Media with Unsteady Phase Permeabilities

The influence of nonequilibrium effects developing due to the formation of an emulsion of each phase (gas bubbles in the water and water dust in the gas) on the flow dynamics is investigated with reference to the displacement of water by a gas. The nonequilibrium effects manifest themselves in a change in the shape of the phase permeability curves (they become convex) and the threshold phase saturations in the course of flow through the porous medium. A kinetic equation in which the relaxation time is proportional to the seepage rate is used to describe such effects. The case in which the liquid displaced by the gas is itself gassed and the volume concentration of the gas bubbles is constant is considered.     

Unified Model of Drainage and Imbibition in 3D Fractionally Wet Porous Media

We develop a grain-based model for capillarity controlled displacement within 3D fractionally wet porous media. The model is based on a novel local calculation of the position of stable fluid–fluid interfaces in contact with multiple spherical grains of arbitrary contact angles. The interface is assumed to be locally spherical between bulk phases; the interface is assumed to be toroidal between pairs of grains (surfaces of pendular rings). Because the calculation of interface position is entirely local and grain-based, it provides a single, generalized, geometric basis for computing pore-filling events during drainage as well as imbibition with both Melrose events (merging of two interfaces) and Haines events (geometric instability). The model is validated against a series of drainage/imbibition experiments (oil/water) on fractionally wet porous media prepared by mixing oil-wet grains with water-wet grains.     

Efficiency of Capillary Imbibition Dominated Displacement of Nonwetting Phase by Wetting Phase in Fractured Porous Media

The critical and optimum injection rates as well as the critical fracture capillary number for an efficient displacement process are determined based on the experimental and numerical modeling of the displacement of nonwetting phase (oil) by wetting phase (water) in fractured porous media. The efficiency of the process is defined in terms of the nonwetting phase displaced from the system per amount of wetting phase injected and per time. Also, the effects of injection rate on capillary imbibition transfer dominated two-phase flow in fractured porous media are clarified by visualizing the experiments. The results reveal that as the injection rate is increased, fracture pattern begins to become an effective parameter on the matrix saturation distribution. As the rate is lowered, however, the system begins to behave like a homogeneous system showing a frontal displacement regardless the fracture configuration.     

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.     

Effect of Organosilicon-Acrylic Emulsion Treatment on Wettability of Porous Media

To investigate the influence of the organosilicon-acrylic on wetting properties of porous media, contact angle tests were performed on two different sandstones. In addition, the effectiveness of the emulsion on wettability alteration of porous media was validated by capillary rise and spontaneous imbibition tests. The results of wettability tests showed that the wettability of two sandstones was altered from water-wet to gas-wet after treatment with the emulsion. The principle that the critical radius of pore throats and wettability of porous media affect liquids flow was derived analytically and verified experimentally. Coreflood results demonstrated that the latex resulted in increasing the water permeability through altering the rock wettability to gas-wetting, then decreasing the friction drag between liquids and rocks surface. Thereby, the emulsion treatment could increase the flowback rate of trapped liquids. Experimental results were in good agreement with the theoretical analysis. In conclusion, all results indicated that the emulsion could alter the wettability from water-wet to intermediate gas-wet and enhance water permeability in porous media. It was extrapolated that the emulsion had the tremendous potential to be applied in field conditions, enhancing gas productivity through the cleanup of trapped water in the vicinity of the wellbore.     

Effect of Efflorescence Formation on Drying Kinetics of Porous Media

We study experimentally the evaporation of an aqueous NaCl solution with efflorescence formation in hydrophilic or hydrophobic two-dimensional model porous media. The efflorescence formation only marginally affects the invasion patterns but greatly modifies the drying kinetics compared to pure water. Two mechanisms are identified for explaining the impact of efflorescence formation on drying kinetics. It is shown that the efflorescence contributes to increase the evaporation rate compared to drying with pure water, a surprising result since the water activity is reduced in the presence of dissolved salt. This effect is explained by the efflorescence liquid capillary pumping effect associated with the porous nature of the efflorescence. Then, we identify a second phase where the efflorescence dries out and acts as a vapour diffusion barrier, leading to a dramatic reduction of evaporation rate.     

The Effect of Isopropyl Alcohol and Non-Ionic Surfactant Mixtures on the Wetting of Porous Coated Paper

The influence of isopropyl alcohol and non-ionic surfactant solutions on aqueous droplet wetting behaviour on porous coated paper was determined. Paper coatings provide a micro- and nano-porous surface structure, which strictly speaking cannot be described in simple roughness terms as sub-surface lateral absorption directly impacts on the apparent contact angle. It is this very deviation from an idealised system that leads to novel wetting phenomena. Isopropyl alcohol and surfactant-based systems, both of which are commonly used in the printing industry, show differences in wetting behaviour, on both short and long timescales, with changes in the relative composition of the mixtures. Small variations of 0.1?wt% in surfactant concentration have a dramatic influence on the dynamic surface tension, and thus the wetting. It was observed that the wetting kinetics for isopropyl alcohol and surfactant solutions were different in terms of both wetting area and the penetration rate, even in cases where the dynamic surface tension of the solutions was kept the same. Different stages in the wetting and following drying processes could be observed with near infrared spectral imaging. In addition, the surfactant chemistries such as their degrees of hydrophilicity and molecular weights generated comparative differences in the wetting kinetics. The dominating factor affecting the wetting was, as expected, the solid?Cliquid interfacial energy defined on the practical porous substrate, which differed from the direct comparison with dynamic surface tension, thus exemplifying the deviation from idealised surface roughness behaviour when considering porous materials. An apparent ??equivalent?? surface roughness value for the porous material was determined, and it was seen that an increase in this equivalent parameter enhanced the rate of wetting behaviour with decreasing solution surface tension, and so also affected the wetting evolution. The wetting was enhanced by cavities in the coating layer, which were enlarged by the penetrating liquids.