Effect of soil moisture dynamics on dense nonaqueous phase liquid (DNAPL) spill zone architecture in heterogeneous porous media

The amount, location, and form of NAPL in contaminated vadose zones are controlled by the spatial distribution of water saturation and soil permeability, the NAPL spill scenario, water infiltration events, and vapor transport. To evaluate the effects of these processes, we used the three-phase flow simulator STOMP, which includes a new permeability-liquid saturation-capillary pressure (k-S-P) constitutive model. This new constitutive model considers three NAPL forms: free, residual, and trapped. A 2-D vertical cross-section with five stratigraphic layers was assumed, and simulations were performed for seven cases. The conceptual model of the soil heterogeneity was based upon the stratigraphy at the Hanford carbon tetrachloride (CT) spill site. Some cases considered co-disposal of NAPL with large volumes of wastewater, as also occurred at the Hanford CT site. In these cases, the form and location of NAPL were most strongly influenced by high water discharge rates and NAPL evaporation to the atmosphere. In order to investigate the impact of heterogeneity, the hydraulic conductivity within the lower permeability layer was modeled as a realization of a random field having three different classes. For six extreme cases of 100 realizations, the CT mass that reached the water table varied by a factor of two, and was primarily controlled by the degree of lateral connectivity of the low conductivity class within the lowest permeability layer. The grid size at the top boundary had a dramatic impact on NAPL diffusive flux just after the spill event when the NAPL was present near the ground surface. NAPL evaporation with a fine grid spacing at the top boundary decreased CT mass that reached the water table by 74%, compared to the case with a coarse grid spacing, while barometric pumping had a marginal effect for the case of a continuous NAPL spill scenario considered in this work. For low water infiltration rate scenarios, the distribution of water content prior to a NAPL spill event decreased CT mass that reached the water table by 98% and had a significant impact on the formation of trapped NAPL. For all cases simulated, use of the new constitutive model that allows the formation of residual NAPL increased the amount of NAPL retained in the vadose zone. Density-driven advective gas flow from the ground surface controlled vapor migration in strongly anisotropic layers, causing NAPL mass flux to the lower layer to be reduced. These simulations indicate that consideration of the formation of residual and trapped NAPLs and dynamic boundary conditions (e.g., areas, rates, and periods of different NAPL and water discharge and fluctuations of atmospheric pressure) in the context of full three-phase flow are needed, especially for NAPL spill events at the ground surface. In addition, NAPL evaporation, density-driven gas advection, and NAPL vertical movement enhanced by water flow must be considered in order to predict NAPL distribution and migration in the vadose zone.     

Impact of nonaqueous phase liquid (NAPL) source zone architecture on mass removal mechanisms in strongly layered heterogeneous porous media during soil vapor extraction

An existing multiphase flow simulator was modified in order to determine the effects of four mechanisms on NAPL mass removal in a strongly layered heterogeneous vadose zone during soil vapor extraction (SVE): a) NAPL flow, b) diffusion and dispersion from low permeability zones, c) slow desorption from sediment grains, and d) rate-limited dissolution of trapped NAPL. The impacts of water and NAPL saturation distribution, NAPL-type (i.e., free, residual, or trapped) distribution, and spatial heterogeneity of the permeability field on these mechanisms were evaluated. Two different initial source zone architectures (one with and one without trapped NAPL) were considered and these architectures were used to evaluate seven different SVE scenarios. For all runs, slow diffusion from low permeability zones that gas flow bypassed was a dominant factor for diminished SVE effectiveness at later times. This effect was more significant at high water saturation due to the decrease of gas-phase relative permeability. Transverse dispersion contributed to fast NAPL mass removal from the low permeability layer in both source zone architectures, but longitudinal dispersion did not affect overall mass removal time. Both slow desorption from sediment grains and rate-limited mass transfer from trapped NAPL only marginally affected removal times. However, mass transfer from trapped NAPL did affect mass removal at later time, as well as the NAPL distribution. NAPL flow from low to high permeability zones contributed to faster mass removal from the low permeability layer, and this effect increased when water infiltration was eliminated. These simulations indicate that if trapped NAPL exists in heterogeneous porous media, mass transfer can be improved by delivering gas directly to zones with trapped NAPL and by lowering the water content, which increases the gas relative permeability and changes trapped NAPL to free NAPL.     

Impact of varying soil structure on transport processes in different diagnostic horizons of three soil types

When soil structure varies in different soil types and the horizons of these soil types, it has a significant impact on water flow and contaminant transport in soils. This paper focuses on the effect of soil structure variations on the transport of pesticides in the soil above the water table. Transport of a pesticide (chlorotoluron) initially applied on soil columns taken from various horizons of three different soil types (Haplic Luvisol, Greyic Phaeozem and Haplic Cambisol) was studied using two scenarios of ponding infiltration. The highest infiltration rate and pesticide mobility were observed for the Bt₁ horizon of Haplic Luvisol that exhibited a well-developed prismatic structure. The lowest infiltration rate was measured for the Bw horizon of Haplic Cambisol, which had a poorly developed soil structure and a low fraction of large capillary pores and gravitational pores. Water infiltration rates were reduced during the experiments by a soil structure breakdown, swelling of clay and/or air entrapped in soil samples. The largest soil structure breakdown and infiltration decrease was observed for the Ap horizon of Haplic Luvisol due to the low aggregate stability of the initially well-aggregated soil. Single-porosity and dual-permeability (with matrix and macropore domains) flow models in HYDRUS-1D were used to estimate soil hydraulic parameters via numerical inversion using data from the first infiltration experiment. A fraction of the macropore domain in the dual-permeability model was estimated using the micro-morphological images. Final soil hydraulic parameters determined using the single-porosity and dual-permeability models were subsequently used to optimize solute transport parameters. To improve numerical inversion results, the two-site sorption model was also applied. Although structural changes observed during the experiment affected water flow and solute transport, the dual-permeability model together with the two-site sorption model proved to be able to approximate experimental data.     

土壤柴油污染修复的抽气提取去除实验研究

为得到土壤气相抽提(SVE)去除柴油的优化条件，进行了一维土柱抽气提取去除柴油污染物的实验研究，研究不同初始含水率、不同抽气量对污染土壤中柴油去除率的影响及不同深度残留柴油的变化规律。结果表明：在本实验模拟的范围内，抽气量越大，SVE处理效果越好；初始含水率越低，处理效果越好；此外，不同深度去除率变化的规律基本上是随深度的增大而减小。实验结果可为土壤轻油污染实际治理提供实验数据基础。     

Conceptual and numerical model of uranium(VI) reductive immobilization in fractured subsurface sediments

A conceptual model and numerical simulations of bacterial U(VI) reduction in fractured subsurface sediments were developed to assess the potential feasibility of biomineralization at the fracture/matrix interface as a mechanism for immobilization of uranium in structured subsurface media. The model envisions flow of anaerobic groundwater, with or without acetate as an electron donor for stimulation of U(VI) reduction by dissimilatory metal-reducing bacteria (DMRB), within mobile macropores along a one-dimensional flow path. As the groundwater moves along the flow path, U(VI) trapped in the immobile mesopore and micropore domains (the sediment matrix) becomes desorbed and transferred to the mobile macropores (fractures) via a first-order exchange mechanism. By allowing bacterial U(VI) reduction to occur in the mesopore domain (assumed to account for 12% of total sediment pore volume) according to experimentally-determined kinetic parameters and an assumed DMRB abundance of 10(7) cells per cm3 bulk sediment (equivalent to 4 mg of cells per dm3 bulk sediment), the concentration of U(VI) in the macropore domain was reduced ca. 10-fold compared to that predicted in the absence of mesopore DMRB activity after a 6-month simulation period. The results suggest that input of soluble electron donors over a period of years could lead to a major redistribution of uranium in fractured subsurface sediments, converting potentially mobile sorbed U(VI) to an insoluble reduced phase (i.e. uraninite) in the mesopore domain that is expected to be permanently immobile under sustained anaerobic conditions.     

有机污染土壤通风去污技术研究进展

土壤气体抽排净化技术是一种安全、高效、经济的有机污染土壤治理技术。本文系统分析了土壤气体抽排净化技术的特点、现状及发展趋势 ;总结了估算污染物净化时间的数学模型 ;概括了影响土壤气体抽排净化技术净化时间的主要因素 ,包括土壤透气率、含水率、有机污染物气体饱和蒸气压、抽排流速、环境温度 ;本文还对目前土壤气体抽排净化技术需要研究的问题及发展前景进行了展望。     

Petroleum mass removal from low permeability sediment using air sparging/soil vapor extraction: impact of continuous or pulsed operation

Air sparging and soil vapor extraction (AS/SVE) are innovative remediation techniques that utilize volatilization and microbial degradation to remediate petroleum spills from soils and groundwater. This in situ study investigated the use of AS/SVE to remediate a gasoline spill from a leaking underground storage tank (UST) in the low permeability, clayey soil of the Appalachian Piedmont. The objectives of this study were to evaluate AS/SVE in low permeability soils by quantifying petroleum mass removal rates, monitoring vadose zone contaminant levels, and comparing the mass extraction rates of continuous AS/SVE to 8 and 24 h pulsed operation. The objectives were met by collecting AS/SVE exhaust gas samples and vadose zone air from multi-depth soil vapor probes. Samples were analyzed for O₂, CO₂, BTEX (benzene, toluene, ethylbenzene, xylene), and total combustible hydrocarbon (TCH) concentrations using portable hand meters and gas chromatography. Continuous AS/SVE was effective in removing 608 kg of petroleum hydrocarbons from low permeability soil in 44 days (14.3 kg day⁻¹). Mass removal rates ranged from 2.6 times higher to 5.1 times lower than other AS/SVE studies performed in sandy sediments. BTEX levels in the vadose zone were reduced from about 5 ppm to 1 ppm. Ten pulsed AS/SVE tests removed 78 kg in 23 days and the mean mass removal rate (17.6 kg day⁻¹) was significantly higher than the last 15 days of continuous extraction. Pulsed operation may be preferable to continuous operation because of increased mass removal and decreased energy consumption.     

Effect of soil layering on NAPL removal behavior in soil-heated vapor extraction

Soil heating has been proposed as a method to enhance the vapor extraction of NAPLs from contaminated soils. Three-dimensional fluid flow and heat transfer simulations have been performed for soil-heated vapor extraction to determine the transient system performance for a hypothetical configuration. Soil layering has been considered in evaluation of the initial non-aqueous phase liquid (NAPL) distribution and in evaporation and transport to the vapor extraction location. Results from this layered model are compared with results for a homogeneous system with an initially uniform NAPL, indicating the influence of layering, the initial NAPL distribution, the type of NAPL, and the possibility of enhanced vapor diffusion. Not only is the NAPL removal time reduced significantly with the addition of heat, but the uncertainty in the removal time owing to a number of difficult to characterize in situ factors, such as layering and the initial NAPL distribution, is much less than for standard soil vapor extraction without heating, owing to the rise in temperature and increase in NAPL vapor pressure with time.     

A multi-flowpath model for the interpretation of immiscible displacement experiments in heterogeneous soil columns

This work focuses on the phenomenon of the immiscible two-phase flow of water and oil in saturated heterogeneous soil columns. The goal is to develop a fast and reliable method for quantifying soil heterogeneities for incorporation into the relevant capillary pressure and relative permeability functions. Such data are commonly used as input data in simulators of contaminant transport in the subsurface. Rate-controlled drainage experiments are performed on undisturbed soil columns and the transient response of the axial distribution of water saturation is determined from electrical measurements. The transient responses of the axial distribution of water saturation and total pressure drop are fitted with the multi-flowpath model (MFPM) where the pore space is regarded as a system of parallel paths of different permeability. The MFPM enables us to quantify soil heterogeneity at two scales: the micro-scale parameters describe on average the effects of pore network heterogeneities on the two-phase flow pattern; the macro-scale parameters indicate the variability of permeability at the scale of interconnected pore networks. The capillary pressure curve is consistent with that measured with mercury intrusion porosimetry over the low pressure range. The oil relative permeability increases sharply at a very low oil saturation (< 10^− 3) and tends to a high end value. The water relative permeability decreases abruptly at a low oil saturation (~ 0.1), whereas the irreducible wetting phase saturation is quite high. The foregoing characteristics of the two-phase flow properties are associated with critical (preferential) flowpaths that comprise a very small percentage of the total pore volume, control the overall hydraulic conductivity, and are consistent with the very broad range of pore-length scales usually probed in soil porous matrix.     

One-at-a-time sensitivity analysis of pollutant loadings to subsurface properties for the assessment of soil and groundwater pollution potential

Chemical leak was numerically simulated for four chemical substances: benzene (light non-aqueous phase liquid (NAPL)), tetrachloroethylene (dense NAPL), phenol (soluble in water), and pentachlorophenol (white crystalline solid) in a hypothetical subsurface leak situation using a multiphase compositional transport model. One metric ton of chemical substances was assumed to leak at a point 3.51 m above the water table in a homogeneous unconfined aquifer which had the depth to water table of 7.135 m, the hydraulic gradient of 0.00097, the recharge rate of 0.7 mm/day, and the permeability of 2.92?×?10^?10 m². For comparison, surface spill scenarios, which had a long pathway from source to the water table, were simulated. Using the model results, point-source pollutant loadings to soil and groundwater were calculated by multiplying mass, impact area, and duration above and below the water table respectively. Their sensitivity to subsurface properties (depth to water table, recharge rate, porosity, organic carbon content, decay rate, hydraulic gradient, capillary pressure, relative permeability, permeability) was analyzed, with changing each parameter within acceptable ranges. The study result showed that the pollutant loading to groundwater was more sensitive to the subsurface properties than the pollutant loading to soil. Decay rate, groundwater depth, hydraulic gradient and porosity were influential to pollutant loadings. The impact of influential parameters on pollutant loadings was nonlinear. The dominant subsurface properties of pollution loadings (e.g., decay rate, groundwater depth, hydraulic gradient, and porosity for groundwater) also affect the vulnerability, and the subsurface pollutant loadings defined in this study are dependent on chemical properties as well, which indicates that the influential hydrogeological and physicochemical parameters to pollutant loadings can be used for pollution potential assessment. The contribution of this work is the suggestion that the sensitivity of pollutant loadings can be used for pollution potential assessment. Soil and groundwater pollution potential of chemicals are discussed altogether for leak scenarios. A physics-based model is used to understand the impact of subsurface properties on the fate and transport of chemicals above and below the water table, and consequently their impact on the pollutant loading to soil and groundwater.

     

Modeling of air sparging of VOC-contaminated soil columns

Air sparging is a remediation technology currently being applied for the restoration of sites contaminated with volatile organic compounds (VOCs). Attempts have been made by various researchers to model the fate of VOCs in the gas and liquid phase during air sparging. In this study, a radial diffusion model with an air–water mass transfer boundary condition was developed and applied for the prediction of VOC volatilization from air sparging of contaminated soil columns. The approach taken was to use various parameters such as mass transfer coefficients and tortuosity factors determined previously in separate experiments using a single air channel apparatus and applying these parameters to a complex system with many air channels. Incorporated in the model, is the concept of mass transfer zone (MTZ) where diffusion of VOCs in this zone was impacted by the volatilization of VOCs at the air–water interface but with negligible impact outside the zone. The model predicted fairly well the change in the VOC concentrations in the exhaust air, the final average aqueous VOC concentration, and the total mass removed. The predicted mass removal was within 1% to 20% of the actual experimental mass removed. The results of the model seemed to suggest that air-sparged soil columns may be modeled as a composite of individual air channels surrounded by a MTZ. For a given air flow rate and air saturation, the VOC removal was found to be inversely proportional to the radius of the air channel. The approach taken provided conceptual insights on mass transfer processes during air sparging operations.     

Dual-domain solute transfer and transport processes: evaluation in batch and transport experiments

Soil macropore networks establish a dual-domain transport scenario in which water and solutes are preferentially channeled through soil macropores while slowly diffusing into and out of the bulk soil matrix. The influence of macropore networks on intra-ped solute diffusion and preferential transport in a soil typical of subsurface-drained croplands in the Midwestern United States was studied in batch- and column-scale experiments. In the batch diffusion studies with soil aggregates, the estimated diffusion radius (length) of the soil aggregates corresponded to the half-spacing of the aggregate fissures, suggesting that the intra-ped fissures reduced the diffusion impedance and preferentially allowed solutes to diffuse into the soil matrix. In the column-scale solute transport experiments, the average diffusion radius (estimated from HYDRUS-2D simulations and a first-order diffusive transfer term) was nearly double that of the batch-scale study. This increase may be attributed to a loss of pore continuity and a compounding of the small diffusion impedance through macropores at the larger scale. The column-scale solute transport experiments also suggest that two preferential networks exist in the soil. At and near soil saturation, a primary network of large macropores (possibly root channels and earthworm burrows) dominate advective transport, causing a high degree of physical and sorption nonequilibrium and simultaneous breakthrough of a nonreactive (bromide) and a reactive (alachlor) solute. As the saturation level decreases, the primary network drains, while transport through smaller macropores (possibly intra-ped features) continues, resulting in a reduced degree of nonequilibrium and separation in the breakthrough curves of bromide and alachlor.     

Effect of pH control at the anode for the electrokinetic removal of phenanthrene from kaolin soil

Polycyclic aromatic hydrocarbon (PAH)-contaminated soils exist at numerous sites, and these sites may threaten public health and the environment because many PAH compounds are toxic, mutagenic, and/or carcinogenic. PAHs are also hydrophobic and persistent, so conventional remediation methods are often costly or inefficient, especially when the contaminants are present in low permeability and/or organic soils. An innovative technique, electrokinetically enhanced in situ flushing, has the potential to increase soil-solution-contaminant interaction and PAH removal efficiency for low permeability soils; however, the electrolysis reaction at the anode may adversely affect the remediation of low acid buffering capacity soils, such as kaolin. Therefore, the objective of this study was to improve the remediation of low acid buffering soils by controlling the pH at the anode to counteract the electrolysis reaction. Six bench-scale electrokinetic experiments were conducted, where each test employed one of three different flushing solutions, deionized water, a surfactant, or a cosolvent. For each of these solutions, tests were performed with and without a 0.01 M NaOH solution at the anode to control the pH. The test using deionized water with pH control generated a higher electroosmotic flow than the equivalent test performed without pH control, but the electroosmotic flow difference between the surfactant and cosolvent tests with and without pH control was minor compared to that observed with the deionized water tests. Controlling the pH was beneficial for increasing contaminant solubilization and migration from the soil region adjacent to the anode, but the high contaminant concentrations that resulted in the middle or cathode soil regions indicates that subsequent changes in the soil and/or solution chemistry caused contaminant deposition and low overall contaminant removal efficiency.     

Modeling radium and radon transport through soil and vegetation

A one-dimensional flow and transport model was developed to describe the movement of two fluid phases, gas and water, within a porous medium and the transport of 226Ra and 222Rn within and between these two phases. Included in this model is the vegetative uptake of water and aqueous 226Ra and 222Rn that can be extracted from the soil via the transpiration stream. The mathematical model is formulated through a set of phase balance equations and a set of species balance equations. Mass exchange, sink terms and the dependence of physical properties upon phase composition couple the two sets of equations. Numerical solution of each set, with iteration between the sets, is carried out leading to a set-iterative compositional model. The Petrov-Galerkin finite element approach is used to allow for upstream weighting if required for a given simulation. Mass lumping improves solution convergence and stability behavior. The resulting numerical model was applied to four problems and was found to produce accurate, mass conservative solutions when compared to published experimental and numerical results and theoretical column experiments. Preliminary results suggest that the model can be used as an investigative tool to determine the feasibility of phytoremediating radium and radon-contaminated soil.     

Water and solute movement in soil as influenced by macropore characteristics: 2. Macropore tortuosity

The paper describes the results of a laboratory study on the effects of macropore tortuosity on breakthrough curves BTCs and solute distribution in a Forman loam (fine loamy-mixed Udic Haploborolls) soil. BTC were obtained using 2-D columns (slab) containing artificial macropores of five different tortuosity levels. The BTCs were run under a constant hydraulic head of 0.08 m over an initially air dry soil. The input solutions contained 1190 mg l⁻¹ of potassium bromide, 10 mg l⁻¹ of Rhodamine WT, and 100 mg l⁻¹ of FD&C Blue #1. A soil column without macropores served as a control. The displacement of a non-adsorbed tracer was not affected by the tortuosity level. An increase in macropore tortuosity progressively increased the breakthrough time, increased the apparent retardation coefficient (R′), decreased the depth to the center of mass of a given adsorbed tracer, and increased the anisotropy in tracer distribution profile. The relative importance of macropore tortuosity increased with an increase in the adsorption coefficient of the tracer. Compared to macropore continuity, the macropore tortuosity had greater impact on solute distribution profile than in its leaching.     

Kinematic wave approximation to solute transport along preferred flow paths in soils

A routing procedure is introduced which accounts for the loss of a conservative solute tracer from preferred paths during macropore flow. Water flow is treated as a series of kinematic waves from which the tracer is lost due to mixing previously stored soil water, and an expression for solute loss is added to a previously developed model. The model parameters are estimated through experiments at three different input rates applied to a column of a macroporous forest soil.The results of seven experimental runs indicate that solute losses are consistently highest at the early stages of infiltration and drainage flow. An empirical relationship is proposed which links the frequency distribution of the flow parameter with that for solute loss from the preferred path during transient water flow and solute transport.