Radon Emanation in Partially Saturated Porous Media

Radon emanation which is the fraction of radon-222 atoms released in the connected pore space of a porous material, increases when the water content increases because of the low recoil range of this atom in water compared with air. This complex phenomenon is studied using reconstructed porous media and random packings where the phase distribution is obtained by a lattice Boltzmann technique incorporating interfacial tension and wetting. The influence of the pore structure, of the recoil range and of saturation, is systematically studied. The numerical results are in good agreement with data from uranium mine tailings.     

NMR Measurements of Tortuosity in Partially Saturated Porous Media

The tortuosity (τ), defined in the present context as the ratio of the free diffusion coefficient to the restricted diffusion coefficient of a contained fluid, is an important but difficult to measure characteristic of a porous medium, particularly when it is partially saturated with water. We develop and apply methodology, based on nuclear magnetic resonance (NMR) pulsed field gradient techniques, to measure τ for various sandstone rock cores as a function of residual water fraction. The NMR methodology requires the use of bipolar pulsed field gradient stimulated echo pulse sequences to avoid systematic errors due to magnetic susceptibility differences and D₂O as a stationary immiscible water phase; this was selected as it provides no ¹H NMR signal. Tortuosity of the free pore space was successfully measured using liquid ethane as a probe fluid for three different sandstones over the full accessible range of residual water saturation. Generally, the tortuosity was observed to increase with residual water (D₂O) content; however, significant variations were observed between the different sandstones.     

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.     

Co-colloidal BTEX and Microbial Transport in a Saturated Porous System: Numerical Modeling and Sensitivity Analysis

Transport in Porous Media - A large number of studies were carried out in the past to analyze the significance of colloids and microbes in remediating groundwater aquifers contaminated with...     

Three-Dimensional Analytical Models for Virus Transport in Saturated Porous Media

Analytical models for virus transport in saturated, homogeneous porous media are developed. The models account for three-dimensional dispersion in a uniform flow field, and first-order inactivation of suspended and deposited viruses with different inactivation rate coefficients. Virus deposition onto solid particles is described by two different processes: nonequilibrium adsorption which is applicable to viruses behaving as solutes; and colloid filtration which is applicable to viruses behaving as colloids. The governing virus transport equations are solved analytically by employing Laplace/Fourier transform techniques. Instantaneous and continuous/periodic virus loadings from a point source are examined.     

A Numerical Study of the Distribution of Water in Partially Saturated Porous Rock

The distribution of water and air phases in small blocks of porous sandstone is examined by using a simulated annealing technique that finds the minimum interfacial energy distributions at different saturations. Simulations are based on existing sandstone microstructures that were determined by X-ray microtomography. At low saturations, some of the water is distributed in films along the walls of larger pore spaces, and connects to pendular structures in the crevices and smaller pores. As the amount of water in the pores increases the water films become thicker and pores fill from the pendular structures. The distribution of water voxels in the pore space is examined by calculating interfacial areas, by classifying water voxels as to whether they lie within films or clusters, and by determining the size and distribution of these film clusters. An exponential relationship is found between the fraction of water voxels in the films and the degree of saturation. In addition, the dependency of small-sample electrical conductivity on saturation is examined by using a random walk method.     

Simulation of Cavitation in Water Saturated Porous Media Considering Effects of Dissolved Air

An extension of a mathematical model for non-isothermal multiphase materials to consider the dissolution of air in liquid water and air mass sources during its desorption at lower water pressure is presented. The solid skeleton is assumed elasto-plastic; heat, water and air flows and water phase changes are taken into account. Physics of air dissolution and desaturation due to the air released from liquid water during cavitation in porous media are discussed. A numerical example where cavitation develops during shear band development in undrained water-saturated dense sands is solved with the developed model as discretized in space and time with the Finite Element Method.     

Stochastic Simulations of NAPL Mass Transport in Variably Saturated Heterogeneous Porous Media

A multiphase flow and transport numerical model is developed to study the effects of porous media heterogeneities on residual NAPL mass partitioning and transport of dissolved and/or volatilized NAPL mass in variably saturated media. The results indicate the significance of porous media heterogeneity in influencing the mass transfer processes and NAPL transport in the subsurface. Among the parameters investigated in this study, the heterogeneity of the permeability field has the most significant influence on the NAPL mass partitioning and transport. In general, the heterogeneity of the porous media properties enhances the NAPL mass plume spreading in both the water phase and the gas phase while the influence on the water phase is much more significant. Overall, the porous media property heterogeneities tend to increase the accumulation of NAPL mass in the water phase. The nonequilibrium mass transfer processes result in the expected trend of decreasing the NAPL mass dissipation rate and causing long-term groundwater contamination.     

Effects of Particle Size Non-Uniformity on Transport and Retention in Saturated Porous Media

Short-pulse injection experiments are investigated to study the effects of particle size non-uniformity on the transport and retention in saturated porous media. Monodisperse particles (3, 10, and 16 \(\upmu \hbox {m}\) latex microspheres) and polydisperse particles (containing 3, 10, and 16 latex microspheres) were explored. The obtained results suggest considering not only the particle sizes but also their polydispersivity (particle size non-uniformity) in transport and retention. Although, the density of the suspended particles is close to that of water, results reveal a slow transport of particles compared to the dissolved tracer whatever their size and flow velocity. The recovered particles in the mixture experiments show that the retention of large particles (10 and 16 \(\upmu \hbox {m}\)) enhances the retention of small ones (3 \(\upmu \hbox {m}\)). However, the straining of 10 and 16 \(\upmu \hbox {m}\) particles in “mixture experiments” is smaller than their straining in “monodisperse experiments”. A linear relationship summarizing the simultaneous effect of particle sizes and flow velocity on deposition kinetics coefficient is proposed.     

Modeling Chemotactic Waves in Saturated Porous Media using Adaptive Mesh Refinement

Bacterial transport is heavily influenced by chemical gradients and interfaces that exist in the subsurface. The main aim of this article is to describe a method of simulating the propagation of a traveling bacterial wave in a contaminated region and the resulting degradation of the contaminant. The presence of the chemotactic term and the relatively small bacterial diffusion means that the wave contains a very sharp wavefront. We, therefore, use an upwind conservative numerical scheme to obtain accurate and numerically stable solutions. The accuracy of the method is verified by comparisons with an exact one-dimensional solution of a simplified problem to give the same wavespeed. The method is then used to simulate the propagation of a realistic chemotactic wave in one dimension. We then use adaptive mesh refinement (AMR) to compute the propagation of chemotactic waves in two dimensions using the simplified model calibrated to give the same wavespeed as the full model.     

A Phenomenological Approach to Modeling Transport in Porous Media

The objective of this article is to make use of the phenomenological approach to construct models for the transport of extensive quantities, such as mass of a fluid phase, mass of a component of a fluid phase, momentum of a phase and energy, in porous medium domains. Special attention is devoted to express the fluxes of these extensive quantities, especially the non-advective ones, as functions of their relevant driving forces, obeying the principle of minimum entropy production. It is shown that for each extensive quantity, we have a linear diffusive flux term, a non-linear diffusive term, and a dispersive flux term. The latter is shown to be proportional to the velocity squared. In each case, the number of moduli that describe fluid and porous matrix properties is determined. The momentum balance equation for a porous medium domain, which is the “motion equation,” is analyzed and simplified for special cases, leading to Darcy’s law and to Brinkman’s equation.     

Thermodynamics, Kinetics, and Hydrodynamics of Mixed Salt Precipitation in Porous Media: Model Development and Parameter Estimation

Demands for hydrocarbon production have been increasing in recent years. Today many oilfields around the world are afflicted by the problem of scaling leading to severe formation damage and hampering of petroleum production from hydrocarbon reservoirs. In current study, a mathematical model for prediction of permeability reduction due to scale deposition is developed based on thermodynamics, kinetics, and hydrodynamics of mixed salt precipitation during flow through porous media. Model predictions are compared with sound experimental data for single deposition of barium sulfate and most importantly, for simultaneous precipitation of barium sulfate and strontium sulfate onto rock surface. Owing to high nonlinearity of the proposed model, kinetic parameters embedded in the mathematical model were tuned employing a new approach based on a hybrid algorithm consisting of particle swarm optimization (PSO) technique and pattern search (PS) algorithm. The average absolute deviations ranging from 1.03 to 9.3 % were observed between model forecasts and experimental data which corroborate the suitability and applicability of the model and also confirm the capability of PSO–PS hybrid algorithm as a highly efficient optimization tool. Estimated values for kinetic parameters are also in accordance with collision theory of chemical reactions.     

Predicting Vertical Flow Barriers Using Tracer Diffusion in Partially Saturated,Layered Porous Media

Sudden changes in isotopic tracer concentration in pore waters have been interpreted as indicating barriers to vertical advective flow through porous rocks in the subsurface, e.g. step changes in \(^{87}\hbox {Sr}/^{86}\) Sr ratio are often used in the oil and gas industry as a signature of reservoir compartmentalisation. This study shows that this is not necessarily the case. It can take millions of years for such step changes to equilibrate by diffusion if there is no flow resulting from pressure or density gradients even in high permeability, high porosity rocks, particularly if the water saturation is low. Changes in tracer concentration gradients can be good indicators of changes in porosity (or water saturation) between layers. In contrast changes in sorption without a change in porosity are almost impossible to identify. The time taken for concentration gradients to equilibrate is affected by the layer properties but can be quickly estimated from the harmonic average of the effective diffusion coefficient for each layer and a simple analytical expression for a homogeneous system. This was achieved by performing a sensitivity analysis on different layer properties (porosity contrast, saturation contrast, sorption contrast, thickness ratio) using existing analytical solutions for diffusion in layered systems.     

Laboratory Investigation of Variable-Density Flow and Solute Transport in Unsaturated–Saturated Porous Media

In many groundwater systems, fluid density and viscosity may vary in space and time as a function of changes in concentration and temperature of the fluid. When dense groundwater plumes interact with less dense ambient groundwater, these density variations can significantly affect flow and transport processes. Under certain conditions, gravitational instabilities in the form of lobe-shaped fingers can occur. This process is significant because it can lead to more rapid and spatially extensive solute transport. This paper presents new experiments carried out in a sand filled glass flow container under both fully saturated and variably-saturated conditions and focuses upon the processes that occur at the capillary fringe and below the water table, as affected by a dense contaminant plumes migration through the unsaturated zone. Source fluids stained with Rhodamine-WT were introduced at the upper boundary of the tank at a range of low and high densities. In addition to the fluid density gradients and porous medium permeability that determine the onset conditions for instabilities in fully saturated experiments, volumetric water content appears critical in the variably-saturated laboratory runs. Plume behaviour at the water table appears dependent upon the density of the fluid that accumulates there. For neutral and low density fluids, plumes accumulate at the water table and then spread laterally above it and the water table forms a barrier to further vertical flow as pore water velocities reduce with increasing water content. For medium and high density fluids, vertical movement continues as instabilities form at the capillary fringe and fingers begin to grow at the water table boundary and move downwards into the saturated zone. In these cases, lateral spreading of the plume is small. Despite these more qualitative observations, the exact nature of the relevant stability criteria for the onset and growth of instabilities in variably-saturated porous media presently remain unclear. All experimental results suggest, however, that the unsaturated zone and position of the water table must be considered in contaminant studies in order to predict the migration pathways, rates and ultimate fate of dense contaminant plumes. It is possible that the results of experiments presented in this paper could form a useful basis for the testing of variable-density (and variably-saturated) groundwater flow and solute transport numerical codes because they offer controlled physical laboratory analogs for comparison. They also provide a strong basis for the development of more rigorous mathematical formulations that are likely to be either developed or tested using numerical flow and solute transport simulators.     

Modeling Low Intensity Coupled Transport Processes in Fluid-Saturated Porous Media

Transport in Porous Media -     

Comparison of Iterative Methods for Improved Solutions of the Fluid Flow Equation in Partially Saturated Porous Media

Abstract. The Picard and modified Picard iteration schemes are often used to numerically solve the nonlinear Richards equation governing water flow in variably saturated porous media. While these methods are easy to implement, they are only linearly convergent. Another approach to solve the Richards equation is to use Newton's iterative method. This method, also known as Newton–Raphson iteration, is quadratically convergent and requires the computation of first derivatives. We implemented Newton's scheme into the mixed form of the Richards equation. As compared to the modified Picard scheme, Newton's scheme requires two additional matrices when the mixed form of the Richards equation is used and requires three additional matrices, when the pressure head-based form is used. The modified Picard scheme may actually be viewed as a simplified Newton scheme.Two examples are used to investigate the numerical performance of different forms of the 1D vertical Richards equation and the different iterative solution schemes. In the first example, we simulate infiltration in a homogeneous dry porous medium by solving both, the h based and mixed forms of Richards equation using the modified Picard and Newton schemes. Results shows that, very small time steps are required to obtain an accurate mass balance. These small times steps make the Newton method less attractive.In a second test problem, we simulate variable inflows and outflows in a heterogeneous dry porous medium by solving the mixed form of the Richards equation, using the modified Picard and Newton schemes. Analytical computation of the Jacobian required less CPU time than its computation by perturbation. A combination of the modified Picard and Newton scheme was found to be more efficient than the modified Picard or Newton scheme.     

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