Multiphase flow and transport through fractured heterogeneous porous media

The migration of Dense, Non-Aqueous Phase Liquid (DNAPL) and dissolved phase contamination through a fractured heterogeneous porous medium has been investigated through the use of a multiphase compositional model. The sensitivity of the timescales of migration and the distribution of contaminant in the subsurface to the mean permeability, the variance of the permeability, and the degree of fracturing of the domain were examined. It was found that increasing the mean permeability of the domain allowed the DNAPL to penetrate deeper into the subsurface, while decreasing the mean permeability caused the DNAPL to pool at shallower depths. The presence of fractures within the system was found to control the infiltration only in the most fractured domain. Moment analysis of the nonwetting phase showed that large-scale movement had ceased after approximately 9 years (maximum duration of the source-on condition was approximately 4.5 years). This tended to be due to a redistribution of the DNAPL towards a residual configuration, as was evidenced by the gradual trending of average nonwetting phase saturations within the domain to a static value. The dissolved phase plume was found to migrate at essentially the same rate as the nonwetting phase, due to the reduced relative permeability of lenses containing DNAPL, and due to diffusive losses of mass to the matrix of fractured clay and silty-clay lenses. Some exceptions to this were found when the DNAPL could not overcome the displacement pressure of a lens, and could not by-pass the lens due to the lack of available driving force after the source had been shut off.     

An experimental study of flow and transport in fractured slate

A series of laboratory tracer migration experiments in a single rock fracture have been performed, and the breakthrough curves have been interpreted using mathematical modelling. Discrepancies were observed between the experimental data and the predictions made using a simple advection-dispersion model. The potential reasons for these discrepancies have been investigated by applying more complex models: one model incorporates channelling of flow within the fracture, the other couples dispersion and advection in the fracture with rock-matrix diffusion. It is concluded that chanelling of flow can adequately explain the observed spreading behaviour; rock-matrix diffusion is not a significant mechanism influencing transport in these experiments.     

A triple-continuum approach for modeling flow and transport processes in fractured rock

This paper presents a triple-continuum conceptual model for simulating flow and transport processes in fractured rock. Field data collected from the unsaturated zone of Yucca Mountain, a repository site of high-level nuclear waste, show a large number of small-scale fractures. The effect of these small fractures has not been considered in previous modeling investigations within the context of a continuum approach. A new triple-continuum model (consisting of matrix, small-fracture, and large-fracture continua) has been developed to investigate the effect of these small fractures. This paper derives the model formulation and discusses the basic triple-continuum behavior of flow and transport processes under different conditions, using both analytical solutions and numerical approaches. The simulation results from the site-scale model of the unsaturated zone of Yucca Mountain indicate that these small fractures may have an important effect on radionuclide transport within the mountain.     

The impact of methanogenesis on flow and transport in coarse sand

The effects of biofilm growth and methane gas generation on water flow in porous media were investigated in an anaerobic two-dimensional sand-filled cell. Inoculation of the lower portion of the cell with a methanogenic culture and addition of methanol to the bottom of the cell led to biomass growth and formation of a gas phase. Biomass distributions in the water and on the sand in the cell were measured by protein analysis. The biofilm distribution on sand was observed by confocal laser scanning microscopy. The formation, migration, distribution and saturation of gases in the cell were visualized by the charge-coupled device (CCD) camera. The effects of biofilm and gas generation on water flow were separated by performing one tracer test in the presence of both biofilm and a gas phase and a second tracer test after removal of the gas phase through water flushing. The results of tracer tests demonstrated that flow and transport in the two-dimensional cell were significantly affected by both gas generation and biofilm growth. Gas generated at the bottom of the cell in the biologically active zone moved upwards in discrete fingers, so that gas phase saturations (gas-filled fraction of void space) in the biologically active zone at the bottom of the cell did not exceed 40-50%, while gas accumulation at the top of the cell produced gas phase saturations as high as 80%. The greatest reductions in water phase permeability, based on measurements of reductions in water phase saturations, occurred near the top of the box as a result of the gas accumulation. In contrast the greatest reductions in permeability due to biofilm growth, based on measurements of biofilm thickness, occurred in the most biologically active zone at the bottom of the cell, where gas phase saturations were approximately 40-50%, but permeability reductions due to biofilm growth were estimated to be 80-95%.     

Modeling of non-reactive solute transport in fractured clayey till during variable flow rate and time

Fractures and biopores can act as preferential flow paths in clay aquitards and may rapidly transmit contaminants into underlying aquifers. Reliable numerical models for assessment of groundwater contamination from such aquitards are needed for planning, regulatory and remediation purposes. In this investigation, high resolution preferential water-saturated flow and bromide transport data were used to evaluate the suitability of equivalent porous medium (EPM), dual porosity (DP) and discrete fracture/matrix diffusion (DFMD) numerical modeling approaches for assessment of flow and non-reactive solute transport in clayey till. The experimental data were obtained from four large undisturbed soil columns (taken from 1.5 to 3.5 m depth) in which biopores and channels along fractures controlled 96-99% of water-saturated flow. Simulating the transport data with the EPM effective porosity model (FRACTRAN in EPM mode) was not successful because calibrated effective porosity for the same column had to be varied up to 1 order of magnitude in order to simulate solute breakthrough for the applied flow rates between 11 and 49 mm/day. Attempts to simulate the same data with the DP models CXTFIT and MODFLOW/MT3D were also unsuccessful because fitted values for dispersion, mobile zone porosity, and mass transfer coefficient between mobile and immobile zones varied several orders of magnitude for the different flow rates, and because dispersion values were furthermore not physically realistic. Only the DFMD modeling approach (FRACTRAN in DFMD mode) was capable to simulate the observed changes in solute transport behavior during alternating flow rate without changing values of calibrated fracture spacing and fracture aperture to represent the macropores.     

Multiphase flow and transport caused by spontaneous gas phase growth in the presence of dense non-aqueous phase liquid

Disconnected bubbles or ganglia of trapped gas may occur below the top of the capillary fringe through a number of mechanisms. In the presence of dense non-aqueous phase liquid (DNAPL), the disconnected gas phase experiences mass transfer of dissolved gases, including volatile components from the DNAPL. The properties of the gas phase interface can also change. This work shows for the first time that when seed gas bubbles exist spontaneous gas phase growth can be expected to occur and can significantly affect water-gas-DNAPL distributions, fluid flow, and mass transfer. Source zone behaviour was observed in three different experiments performed in a 2-dimensional flow cell. In each case, a DNAPL pool was created in a zone of larger glass beads over smaller glass beads, which served as a capillary barrier. In one experiment effluent water samples were analyzed to determine the vertical concentration profile of the plume above the pool. The experiments effectively demonstrated a) a cycle of spontaneous gas phase expansion and vertical advective mobilization of gas bubbles and ganglia above the DNAPL source zone, b) DNAPL redistribution caused by gas phase growth and mobilization, and c) that these processes can significantly affect mass transport from a NAPL source zone.     

Flow and transport in fractured tuff at Yucca Mountain: numerical experiments on fast preferential flow mechanisms

Recent discovery of bomb-related ³⁶Cl at depth in fractured tuff in the unsaturated zone at the Yucca Mountain candidate high-level waste (HLW) repository site has called into question the usual modeling assumptions based on the equivalent continuum model (ECM). A dual continuum model (DCM) for simulating transient flow and transport at Yucca Mountain is developed. In order to ensure properly converged flow solutions, which are used in the transport simulation, a new flow solution convergence criteria is derived. An extensive series of simulation studies is presented which indicates that rapid movement of solute through the fractures will not occur unless there are intense episodic infiltration events. Movement of solute in the environs of the repository is enhanced if the properties of the tuff layer at the repository horizon are modified from current best-estimate values. Due to a large advective–dispersive coupling between the matrix and fractures, the matrix acts as a major buffer which inhibits rapid transport along the fractures. Consequently, fast movement of solutes through the fractures to the repository depth can only be explained if the matrix–fracture coupling term is significantly reduced from a value that would be calculated on the basis of data currently available.     

Particle transport in unsaturated fractured chalk under arid conditions

A series of field and laboratory experiments were conducted to study the mechanisms of particle detachment and transport from fractures in vadose chalk. Experiments of intermittent flow events along fracture surfaces were carried out in the laboratory. In the field, water was percolated from land surface via a discrete fracture into a compartmental sampler installed inside a horizontal corehole located I m below the surface. The mass, size distribution, and composition of the particles drained from the fracture voids were examined along with flow rates and salt dissolution. Two boreholes penetrating the underlying saturated zone were sampled and analyzed for colloidal concentration and composition. Most of the particle and solute release at the drained effluents occurred during the first several hours of flow, but erratic pulses of particles were still observed after long periods of time. Most of the detached particles had a mean diameter of >2 microm, while the mobile colloidal phase in the groundwater had a mean diameter of approximately 1 microm. Mineralogical composition of the groundwater colloids and the particles detached from the upper vadose fracture were similar. Laboratory observations demonstrated the importance of the existence of a coating layer, made of weathered particles and salts, on particle detachment. The results of this study suggest that: (1) particle detachment causes flow-rate variability in the unsaturated fracture; (2) the mechanisms of particle detachment and salt dissolution within the fracture are linked: and (3) although most of the detached particles are large and likely to accumulate inside fractures, some colloidal particles also eroded from the fracture void and are likely to be transported to the groundwater.     

Experimental determination of sorption in fractured flow systems

Fracture "skins" are alteration zones on fracture surfaces created by a variety of biological, chemical, and physical processes. Skins increase surface area, where sorption occurs, compared to the unaltered rock matrix. This study examines the sorption of organic solutes on altered fracture surfaces in an experimental fracture-flow apparatus. Fracture skins containing abundant metal oxides, clays, and organic material from the Breathitt Formation (Kentucky, USA) were collected in a manner such that skin surface integrity was maintained. The samples were reassembled in the lab in a flow-through apparatus that simulated approximately 2.7 m of a linear fracture "conduit." A dual-tracer injection scheme was utilized with the sorbing or reactive tracer compared to a non-reactive tracer (chloride) injected simultaneously. Sorption was assessed from the ratio of the first temporal moments of the breakthrough curves and from the loss of reactive tracer mass and evaluated as a function of flow velocity and solute type. The breakthrough curves suggest dual-flow regimes in the fracture with both sorbing and non-sorbing flow fields. Significant sorption occurs for the reactive components, and sorption increased with decreasing flow rate and decreasing compound solubility. Based on moment analysis, however, there was little retardation of the center of solute mass. These data suggest that non-equilibrium sorption processes dominate and that slow desorption and boundary layer diffusion cause extensive tailing in the breakthrough curves.     

Adsorption of a dye on clay and sand. Use of cyclodextrins as solubility-enhancement agents

Laboratory-scale studies were aimed at elucidating the physico-chemical aspects on the removal process of crystal violet (CV) from waters and solid substrates. The laponite clay (RD) and sand were chosen for the double aim at investigating them as CV adsorbents for water treatment and as substrates which mime the soil components. Sand is very effective in removing CV from waters. The cyclodextrins (CDs) were exploited as solubility-enhancement agents to remove CV from the solid substrates. They are powerful solvent media because they extract the CV from sand forming water-soluble CV/CD inclusion complexes and do not show affinity for sand. Optimum performance was shown by the modified CDs (i.e. hydroxypropyl-β-cyclodextrin and methyl-β-cyclodextrin). A linear correlation between the logarithm of the equilibrium constant for the CV/CD inclusion complexes formation (K_cpx) and the maximum amount of CV extracted from sand in the columns experiments at a flow rate of 1.5 ml min⁻¹ was drawn. This relationship predicts that CDs with K_cpx < 180 M⁻¹ are not suitable for CV removal from sand. CDs failed to displace CV from RD because they generate the formation of RD clusters where CV remains entrapped.     

Adsorption of a dye on clay and sand. Use of cyclodextrins as solubility-enhancement agents

Laboratory-scale studies were aimed at elucidating the physico-chemical aspects on the removal process of crystal violet (CV) from waters and solid substrates. The laponite clay (RD) and sand were chosen for the double aim at investigating them as CV adsorbents for water treatment and as substrates which mime the soil components. Sand is very effective in removing CV from waters. The cyclodextrins (CDs) were exploited as solubility-enhancement agents to remove CV from the solid substrates. They are powerful solvent media because they extract the CV from sand forming water-soluble CV/CD inclusion complexes and do not show affinity for sand. Optimum performance was shown by the modified CDs (i.e. hydroxypropyl-β-cyclodextrin and methyl-β-cyclodextrin). A linear correlation between the logarithm of the equilibrium constant for the CV/CD inclusion complexes formation (K_cpx) and the maximum amount of CV extracted from sand in the columns experiments at a flow rate of 1.5 ml min⁻¹ was drawn. This relationship predicts that CDs with K_cpx < 180 M⁻¹ are not suitable for CV removal from sand. CDs failed to displace CV from RD because they generate the formation of RD clusters where CV remains entrapped.     

Numerical analysis of water and solute transport in variably-saturated fractured clayey till

This study numerically investigates the influence of initial water content and rain intensities on the preferential migration of two fluorescent tracers, Acid Yellow 7 (AY7) and Sulforhodamine B (SB), through variably-saturated fractured clayey till. The simulations are based on the numerical model HydroGeoSphere, which solves 3D variably-saturated flow and solute transport in discretely-fractured porous media. Using detailed knowledge of the matrix, fracture, and biopore properties, the numerical model is calibrated and validated against experimental high-resolution tracer images/data collected under dry and wet soil conditions and for three different rain events. The model could reproduce reasonably well the observed preferential migration of AY7 and SB through the fractured till, although it did not capture the exact depth of migration and the negligible impact of the dead-end biopores in a near-saturated matrix. A sensitivity analysis suggests fast flow mechanisms and dynamic surface coating in the biopores, and the presence of a plough pan in the till.     

A comparison of hydraulic and transport parameters measured in low-permeability fractured media

Data from 90 tracer experiments performed in low-permeability fractured media have been studied to explore correlations among parameters controlling flow and transport. The original data had been interpreted by different authors using different models, which prevents direct comparison of their estimated parameters. In order to produce comparable parameters, the data have been reexamined using simple models (homogeneous domain, steady-state flow regime, single porosity). Specifically, hydraulic conductivity has been derived as the ratio of water flux to head gradient and apparent porosity as the ratio of water velocity to water flux; the former estimated from both first and peak arrival times. Hydraulic conductivity and porosity correlate along a straight line of slope 1:3 in log scale. While the regression is too noisy to be of predictive use, it lends some support to the use of a generalized cubic law. The fact that correlation for first arrival time porosity (0.77) is larger than for peak arrival porosity (0.62) suggests that first arrival is controlled by the same flow paths as hydraulic conductivity. Apparent porosity derived from peak arrival time is found to grow with travel time along a line of 0.55 slope (again log scale). The correlation coefficient ranges between 0.73 and 0.80 (depending on the data set) for hard rocks. The fact that this correlation is maintained when varying the flow rate at a given site leads us to suggest that it is caused by diffusion mechanisms. This conclusion is further supported by the increase of apparent porosity with the matrix porosity of the rock on which the experiments were performed.     

Chromium transport, oxidation, and adsorption in manganese-coated sand

We examine how the processes of advection, dispersion, oxidation-reduction, and adsorption combine to affect the transport of chromium through columns packed with pyrolusite (beta-MnO2)-coated sand. We find that beta-MnO2 effectively oxidizes Cr(III) to Cr(VI) and that the extent of oxidation is sensitive to changes in pH, pore water velocity, and influent concentrations of Cr(III). Cr(III) oxidation rates, although initially high, decline well before the supply of beta-MnO2 is depleted, suggesting that a reaction product inhibits the conversion of Cr(III) to Cr(VI). Rate-limited reactions govern the weak adsorption of each chromium species, with Cr(III) adsorption varying directly with pH and Cr(VI) adsorption varying inversely with pH. The breakthrough data on chromium transport can be matched closely by calculations of a simple model that accounts for (1) advective-dispersive transport of Cr(III), Cr(VI), and dissolved oxygen, (2) first-order kinetics adsorption of the reduced and oxidized chromium species, and (3) nonlinear rate-limited oxidation of Cr(III) to Cr(VI). Our work supplements the limited database on the transport of redox-sensitive metals in porous media and provides a means for quantifying the coupled processes that contribute to this transport.     

Analysis of field observations of tracer transport in a fractured till

We analyze a set of observations from a recently published, field-scale tracer test in a fractured till. These observations demonstrate a dominant, underlying non-Fickian behavior, which cannot be quantified using traditional modeling approaches. We use a continuous time random walk (CTRW) approach which thoroughly accounts for the measurements, and which is based on a physical picture of contaminant motion that is consistent with the geometric and hydraulic characterization of the fractured formation. We also incorporate convolution techniques in the CTRW theory, to consider transport between different regions containing distinct heterogeneity patterns. These results enhance the possibility that limitations in predicting non-Fickian modes of contaminant migration can be overcome.     

A remediation performance model for enhanced metabolic reductive dechlorination of chloroethenes in fractured clay till

A numerical model of metabolic reductive dechlorination is used to describe the performance of enhanced bioremediation in fractured clay till. The model is developed to simulate field observations of a full scale bioremediation scheme in a fractured clay till and thereby to assess remediation efficiency and timeframe. A relatively simple approach is used to link the fermentation of the electron donor soybean oil to the sequential dechlorination of trichloroethene (TCE) while considering redox conditions and the heterogeneous clay till system (clay till matrix, fractures and sand stringers). The model is tested on lab batch experiments and applied to describe sediment core samples from a TCE-contaminated site. Model simulations compare favorably to field observations and demonstrate that dechlorination may be limited to narrow bioactive zones in the clay matrix around fractures and sand stringers. Field scale simulations show that the injected donor is expected to be depleted after 5 years, and that without donor re-injection contaminant rebound will occur in the high permeability zones and the mass removal will stall at 18%. Long remediation timeframes, if dechlorination is limited to narrow bioactive zones, and the need for additional donor injections to maintain dechlorination activity may limit the efficiency of ERD in low-permeability media. Future work should address the dynamics of the bioactive zones, which is essential to understand for predictions of long term mass removal.     

Tylosin sorption to silty clay loam soils, swine manure, and sand

The objectives of this study were to assess sorption and desorption of tylosin, a macrolide antimicrobial chemical used in swine, cattle, and poultry production, in three silty clay loam soils of South Dakota and compare soil sorption to sand and manure sorption. The silty clay loam soils, from a toposequence in eastern South Dakota, standardized sand samples, and swine manure were used in 24-h batch sorption studies with tylosin concentrations ranging from 25 to 232 micro mole/L. Desorption from soil was conducted over a four-day period. Partition coefficients, based on the Freundlich isotherm (K(f)) or K(d) values, were calculated. K(f) values for the silty clay loams were similar, not influenced by landscape position, and averaged 1,350 with isotherm slopes ranging from 0.85 to 0.93. K(f) values for sand were dependent on solution/sand ratios and pH, ranging from 1.4 to 25.1. K(d) values of manure were dependent on the solution type and ranged from 840 L/kg with urine to about 175 L/kg when sorbed from water. Desorption of tylosin from each soil over the four-day period was < 0.2% of the amount added. The soils' high K(f) values and low desorption amounts suggest that once tylosin is in these soils, leaching to lower depths may not occur. However, this does not preclude runoff with soil eroded particles. If tylosin reaches a sand aquifer, through bypass flow or other mechanism(s), movement in the aquifer most likely would occur.     

Numerical modelling of nitrate transport in fractured porous media under non-isothermal conditions

Environmental Science and Pollution Research - Subsurface contamination is a frequent occurrence in fractured porous systems, posing a potential threat for the groundwater contamination. Tracking...     

Description of sorbing tracers transport in fractured media using the lattice model approach

The transport of contaminants in fractured media is a complex phenomenon with a great environmental impact. It has been described with several models, most of them based on complex partial differential equations, that are difficult to apply when equilibrium and nonequilibrium dynamics are considered in complex boundaries. With the aim of overcoming this limitation, a combination of two lattice Bathnagar, Gross and Krook (BGK) models, derived from the lattice Boltzmann model, is proposed in this paper. The fractured medium is assumed to be a single fissure in a porous rock matrix. The proposed approach permits us to deal with two processes with different length scales: advection–dispersion in the fissure and diffusion within the rock matrix. In addition to the mentioned phenomena, sorption reactions are also considered. The combined model has been tested using the experimental breakthrough curves obtained by Garnier et al. (Garnier, J.M., Crampon, N., Préaux, C., Porel, G., Vreulx, M., 1985. Traçage par ¹³C, ²H, I⁻ et uranine dans la nappe de la craie sénonienne en écoulement radial convergent (Béthune, France). J. Hidrol. 78, 379–392.) giving acceptable results. A study on the influence of the lattice BGK models parameters controlling sorption and matrix diffusion on the breakthrough curves shape is included.     

Explicit and implicit coupling during solute transport through clay membrane barriers

Simulations of salt (KCl) flux through a 1-m-thick clay membrane barrier (CMB) based on coupled solute transport theory are compared to simulated fluxes based on traditional advective-dispersive transport theory. The simulations are based on measured values for the effective salt-diffusion coefficient (Ds*) and chemico-osmotic efficiency coefficient (omega) for a bentonite-based barrier material subjected to KCl solutions. The results indicate that the exit salt flux is reduced due to both explicit coupling (hyperfiltration and chemico-osmotic counter-advection) and an implicit coupling effect resulting from the decrease in Ds* due to a decrease in the apparent tortuosity factor, tau a, with an increase in omega. Implicit coupling is shown to be more significant than explicit coupling for reducing and retarding salt flux through a CMB under diffusion-dominated conditions. Failure to account for the implicit coupling effect may result in unrealistic results, such as the existence of salt flux through a perfect (ideal) clay membrane (i.e., omega=1).