Calculating Groundwater Response Times for Flow in Heterogeneous Porous Media

Predicting the amount of time required for a transient groundwater response to take place is a practical question that is of interest in many situations. This time scale is often called the response time. In the groundwater hydrology literature, there are two main methods used to calculate the response time: (1) both the transient and steady‐state groundwater flow equations are solved, and the response time is taken to be amount of time required for the transient solution to approach the steady solution within some tolerance; and (2) simple scaling arguments are adopted. Certain limitations restrict both of these approaches. In this study, we outline a third method, based on the theory of mean action time. We derive the governing boundary value problem for both the mean and variance of action time for confined flow in two‐dimensional heterogeneous porous media. Importantly, we show that these boundary value problems can be solved using widely available software. Applying these methods to a test case reveals the advantages of the theory of mean action time relative to standard methods.     

Flow Characteristics of Microbubble Suspensions in Porous Media as an Oxygen Carrier

Microbubble suspensions were generated as an oxygen carrier for aerobic biodegradation, and their flow characteristics in porous media were investigated. Commercial surfactants including sodium dodecyl sulfate (SDS), and dodecylethyldimethylammonium bromide (DEDAB), saponin (a natural surfactant), and collagen (a protein hydrolysate) were examined as base materials for generating microbubble suspensions. Among them, 2×CMC (critical micellar concentration) of SDS, DEDAB, and saponin developed microbubble suspensions with the highest gas hold‐up and half‐drainage time. Visualization of the flow patterns in sand showed that the microbubble suspensions were separated into a liquid and gas phase directly after injection, showing much faster movement of liquid phase flow. The gas front of the microbubble suspensions flowed in a plug‐flow manner, particularly in cases of SDS and DEDAB. The experimental results from both homogeneous and heterogeneous cells confirmed that the microbubble flow could overcome the heterogeneity in porous media. However, the plug‐flow characteristics and flow propagation of the microbubble suspensions to the low‐permeability zone was accompanied by a large pressure drop, which needs to be considered for future field application.     

An Experimenting Field Approach for the Numerical Solution of Multiphase Flow in Porous Media

In this work, we apply the experimenting pressure field technique to the problem of the flow of two or more immiscible phases in porous media. In this technique, a set of predefined pressure fields are introduced to the governing partial differential equations. This implies that the velocity vector field and the divergence at each cell of the solution mesh can be determined. However, since none of these fields is the true pressure field entailed by the boundary conditions and/or the source terms, the divergence at each cell will not be the correct one. Rather the residue which is the difference between the true divergence and the calculated one is obtained. These fields are designed such that these residuals are used to construct the matrix of coefficients of the pressure equation and the right‐hand side. The experimenting pressure fields are generated in the solver routine and are fed to the different routines, which may be called physics routines, which return to the solver the elements of the matrix of coefficients. Therefore, this methodology separates the solver routines from the physics routines and therefore results in simpler, easy to construct, maintain, and update algorithms.     

Application of the Theory of Extreme Events to Problems of Approximating Probability Distributions of Water Flow Peaks

Statistically treated data of long-term observational series (>30 years) for 66 rivers are used to study the character of distribution of water flow maximums during rain-induced floods in Maritime Territory, Russia. Two probability models are discussed: generalized Pareto distribution (for the domain of very large values) and generalized distribution of extremums (for the rest of the range). The issue of optimal conjunction point of these distributions is discussed. The problem of increasing the accuracy of distribution parameter estimates through data grouping is considered.     

Two-Dimensional Flow and Mass Transport through Variably Saturated Porous Media from a Surface Waste Storage

Water Resources - Numerical calculations in DHI FEFLOW software environment were used to assess the effect of variably saturated zone parameters on capillary barrier formation. The sensitivity of...     

Specification of Matrix Cleanup Goals in Fractured Porous Media

Semianalytical transient solutions have been developed to evaluate what level of fractured porous media (e.g., bedrock or clay) matrix cleanup must be achieved in order to achieve compliance of fracture pore water concentrations within a specified time at specified locations of interest. The developed mathematical solutions account for forward and backward diffusion in a fractured porous medium where the initial condition comprises a spatially uniform, nonzero matrix concentration throughout the domain. Illustrative simulations incorporating the properties of mudstone fractured bedrock demonstrate that the time required to reach a desired fracture pore water concentration is a function of the distance between the point of compliance and the upgradient face of the domain where clean groundwater is inflowing. Shorter distances correspond to reduced times required to reach compliance, implying that shorter treatment zones will respond more favorably to remediation than longer treatment zones in which back‐diffusion dominates the fracture pore water response. For a specified matrix cleanup goal, compliance of fracture pore water concentrations will be reached sooner for decreased fracture spacing, increased fracture aperture, higher matrix fraction organic carbon, lower matrix porosity, shorter aqueous phase decay half‐life, and a higher hydraulic gradient. The parameters dominating the response of the system can be measured using standard field and laboratory techniques.     

Simulation of Tidal Effects on Contaminant Transport in Porous Media

A one-dimensional numerical model is developed with oscillating velocities and dispersions to simulate the migration process of a contaminant plume within tidally influenced aquifers. Model simulations demonstrate that a major effect the tidal fluctuation has on the migration process of a contaminant plume is the exit concentration discharging to the tidal estuary. Tidal fluctuation causes the exit concentration levels to be significantly diluted by the surface-water body of the estuary. Sensitivity analyses demonstrate that tidal fluctuation hastens the rate of plume migration near the bank of the estuary because of the relatively high advective and dispersive fluxes induced by tides. However, tides affect the migration process only over a short distance from the tidal-water interface (about 40 ft for the parameters used in this study). If the contaminant plume is located far beyond the interface, tidal fluctuations will not affect the rate of plume migration until an existing regional ground-water flow velocity brings the plume to the tidally active zone. With or without tides, the rate of contaminant migration increases with higher regional hydraulic gradient. Furthermore, the effects of tidal fluctuations on the transport process become insignificant with higher regional hydraulic gradients.     

Experiment and Simulation of Non-Reactive Solute Transport in Porous Media

To more accurately predict the migration behavior of pollutants in porous media, we conduct laboratory scale experiments and model simulation. Aniline (AN) is used in one-dimensional soil column experiments designed under various media and hydrodynamic conditions. The advection-dispersion equation (ADE) and the continuous-time random walk (CTRW) were used to simulate the breakthrough curves (BTCs) of the solute transport. The results show that the media and hydrodynamic conditions are two important factors affecting solute transport and are related to the degree of non-Fickian transport. The simulation results show that CTRW can more effectively describe the non-Fickian phenomenon in the solute transport process than ADE. The sensitive parameter in the CTRW simulation process is , which can reflect the degree of non-Fickian diffusion in the solute transport. Understanding the relationship of with velocity and media particle size is conducive to improving the reactive solute transport model. The results of this study provide a theoretical basis for better prediction of pollutant transport in groundwater.     

Innovative Techniques for Measuring the Oil Content of Oil-Contaminated Porous Media

The oil content of oil-contaminated porous media is an important parameter for the assessment and remediation of oil pollution in soil and groundwater. However, an accurate measurement method is not available. In this study, we propose a new equation to calculate the oil content of water-bearing media based on traditional extraction–ultraviolet spectroscopy. Further, an improved experimental method was developed. The results indicate that the pure solid weight and oil content of oil-contaminated media can be accurately determined by introducing the oil drying loss coefficient (γ). The average relative errors of the improved method range from −0.14% to −0.96%. They are much smaller than errors of −4.49% to −10.97% of the original methods, indicating that the accuracy of measured oil content is enhanced. In addition, the accuracy of the new method does not depend on oil volatility and oil content.