共查询到20条相似文献，搜索用时 31 毫秒
1.
The time required for dense nonaqueous phase liquid (DNAPL) to cease migrating following release to the subsurface is a valuable component of a site conceptual model. This study uses numerical simulation to investigate the migration of six different DNAPLs in sandy aquifers. The most influential parameters governing migration cessation time are the density and viscosity of the DNAPL and the mean hydraulic conductivity of the aquifer. Releases of between 1 and 40 drums of chlorinated solvent DNAPLs, characterized by relatively high density and low viscosity, require on the order of months to a few years to cease migrating in a heterogeneous medium sand aquifer having an average hydraulic conductivity of 7.4 x 10(-3) cm/s. In contrast to this, the release of 20 drums of coal tar (rho(D)= 1061 kg/m(3), micro(D)= 0.161 Pa.s) requires more than 100 years to cease migrating in the same aquifer. Altering the mean hydraulic conductivity of the aquifer results in a proportional change in cessation times. Parameters that exhibit relatively little influence on migration time scales are the DNAPL-water interfacial tension, release volume, source capillary pressure, mean aquifer porosity, and ambient ground water hydraulic gradient. This study also demonstrates that low-density DNAPLs (e.g., coal tar) give rise to greater amounts of lateral spreading and greater amounts of pooling on capillary barriers than high-density DNAPLs such as trichloroethylene or tetrachloroethylene. 相似文献
2.
Mohamed I.M. Darwish John E. McCray Peter K. Currie Pacelli L.J. Zitha 《Ground Water Monitoring & Remediation》2003,23(2):92-101
Remediation of dense nonaqueous phase liquids (DNAPLs) is recognized as one of the most difficult problems associated with ground water pollution. The pump-and-treat technique, usually consisting of a continuous operation of extraction-injection wells, is widely used for ground water remediation. In a stratified or otherwise heterogeneous aquifer, however, this technique suffers from tailing and rebound problems, which limit its cleanup efficiency and result in higher operation costs. The tailing and rebound is usually due to slow diffusion of contaminants out of lower- permeability heterogeneities into the flow regime of the higher-permeability zone. In this study, we conduct bench-scale experiments to investigate a novel polymer system and injection method to improve the pump-and-treat technique for DNAPL trapped in a layer of porous media that has a relatively low permeability compared to the surrounding media. This technique might be useful, for example, to remove DNAPL from these low-permeability zones after removal of DNAPL from the higher-permeability zones by a more traditional remediation method. The polymer system consists of a mixture of anionic and cationic polyacrylamides in solution and the injection method is based on flow-induced polymer adsorption, called bridging adsorption. The study includes single and parallel-column experiments. The measured polymer penetration depths were compared with values predicted from a numerical simulation, which was developed previously by the authors of this paper. The experiments and simulations show that the polymer injection leads to a modification of the permeability contrast that favors a more efficient pump-and-treat process. These results suggest that additional research to upscale the technology to pilot scales is warranted. 相似文献
3.
《Advances in water resources》1998,22(2):117-132
The migration and entrapment of dense nonaqueous phase liquids (DNAPLs) in aquifer formations is typically believed to be controlled by physical heterogeneities. This belief is based upon the assumption that permeability and capillary properties are determined by the soil texture. Capillarity and relative permeability, however, will also depend on porous medium wettability characteristics. This wettability may vary spatially in a formation due to variations in aqueous phase chemistry, contaminant aging, and/or variations in mineralogy and organic matter distributions. In this work, a two-dimensional multiphase flow simulator is modified to simulate coupled physical and chemical formation heterogeneity. To model physical heterogeneity, a spatially correlated permeability field is generated, and then related to the capillary pressure-saturation function according to Leverett scaling. Spatial variability of porous medium wettability is assumed to be correlated with the natural logarithm of the intrinsic permeability. The influence of wettability on the hysteretic hydraulic property relations is also modeled. The simulator is then employed to investigate the potential influence of coupled physical and chemical heterogeneity on DNAPL flow and entrapment. For reasonable ranges of wettability characteristics, simulations demonstrate that spatial variations in wettability can have a dramatic impact on DNAPL distributions. Higher organic saturations, increased lateral spreading, and decreased depth of infiltration were predicted when the contact angle was varied spatially. When chemical heterogeneity was defined by spatial variation of organic-wet solid fractions (fractional wettability porous media), however, the resultant organic saturation distributions were more similar to those for perfectly water-wet media, due to saturation dependent wettability effects on the hydraulic property relations. 相似文献
4.
The recovery of coal tar or creosote, which are dense nonaqueous phase liquids (DNAPLs), from the subsurface has been used as a means of site remediation at several former manufactured-gas plant sites and wood-treating facilities. Surface-active components in these complex DNAPL can have acid/base and surfactant characteristics that can significantly affect the wettability in these systems, with a reversal to DNAPL wetting at low pHs. These changes in wettability as a function of pH were employed to evaluate their significance on waterflooding efficiency in one-dimensional soil columns. The use of pH as a means of controlling wettability resulted in identical density and viscosity properties between the water-and DNAPL-wet conditions. At some of the higher pHs, the interracial tension changed as well as the wettability. Maintaining a constant dimensionless capillary number was used to minimize the effect of this variable. DNAPL saturations remaining in the DNAPL-wet systems after waterflooding ranged from 38% to 47%. Significantly lower final DNAPL saturations were achieved for water-wetting systems (15% to 30%). 相似文献
5.
E. T. Vogler C. V. Chrysikopoulos 《Stochastic Environmental Research and Risk Assessment (SERRA)》2001,15(1):33-46
A two-dimensional numerical transport model is developed to determine the effect of aquifer anisotropy and heterogeneity
on mass transfer from a dense nonaqueous phase liquid (DNAPL) pool. The appropriate steady state groundwater flow equation
is solved implicitly whereas the equation describing the transport of a sorbing contaminant in a confined aquifer is solved
by the alternating direction implicit method. Statistical anisotropy in the aquifer is introduced by two-dimensional, random
log-normal hydraulic conductivity field realizations with different directional correlation lengths. Model simulations indicate
that DNAPL pool dissolution is enhanced by increasing the mean log-transformed hydraulic conductivity, groundwater flow velocity,
and/or anisotropy ratio. The variance of the log-transformed hydraulic conductivity distribution is shown to be inversely
proportional to the average mass transfer coefficient. 相似文献
6.
《Advances in water resources》2001,24(7):803-819
Alcohol flooding, consisting of injection of a mixture of alcohol and water, is one source removal technology for dense non-aqueous phase liquids (DNAPLs) currently under investigation. An existing compositional multiphase flow simulator (UTCHEM) was adapted to accurately represent the equilibrium phase behavior of ternary and quaternary alcohol/DNAPL systems. Simulator predictions were compared to laboratory column experiments and the results are presented here. It was found that several experiments involved unstable displacements of the NAPL bank by the alcohol flood or of the alcohol flood by the following water flood. Unstable displacement led to additional mixing compared to ideal displacement. This mixing was approximated by a large dispersion in one-dimensional simulations and or by including permeability heterogeneities on a very small scale in three-dimensional simulations. Three-dimensional simulations provided the best match. Simulations of unstable displacements require either high-resolution grids, or need to consider the mixing of fluids in a different manner to capture the resulting effects on NAPL recovery. 相似文献
7.
Chunlin Huang Bill X. Hu Xin Li Ming Ye 《Stochastic Environmental Research and Risk Assessment (SERRA)》2009,23(8):1155-1167
Hydraulic conductivity distribution and plume initial source condition are two important factors affecting solute transport
in heterogeneous media. Since hydraulic conductivity can only be measured at limited locations in a field, its spatial distribution
in a complex heterogeneous medium is generally uncertain. In many groundwater contamination sites, transport initial conditions
are generally unknown, as plume distributions are available only after the contaminations occurred. In this study, a data
assimilation method is developed for calibrating a hydraulic conductivity field and improving solute transport prediction
with unknown initial solute source condition. Ensemble Kalman filter (EnKF) is used to update the model parameter (i.e., hydraulic
conductivity) and state variables (hydraulic head and solute concentration), when data are available. Two-dimensional numerical
experiments are designed to assess the performance of the EnKF method on data assimilation for solute transport prediction.
The study results indicate that the EnKF method can significantly improve the estimation of the hydraulic conductivity distribution
and solute transport prediction by assimilating hydraulic head measurements with a known solute initial condition. When solute
source is unknown, solute prediction by assimilating continuous measurements of solute concentration at a few points in the
plume well captures the plume evolution downstream of the measurement points. 相似文献
8.
Richard T. Gill Steven F. Thornton Michael J. Harbottle Jonathan W.N. Smith 《Ground Water Monitoring & Remediation》2015,35(3):46-56
This study investigates and quantifies the influence of physical heterogeneity in granular porous media, represented by materials with different hydraulic conductivity, on the migration of nitrate, used as an amendment to enhance bioremediation, under an electric field. Laboratory experiments were conducted in a bench‐scale test cell under a low applied direct current using glass bead and clay mixes and synthetic groundwater to represent ideal conditions. The experiments included bromide tracer tests in homogeneous settings to deduce controls on electrokinetic transport of inorganic solutes in the different materials, and comparison of nitrate migration under homogeneous and heterogeneous scenarios. The results indicate that physical heterogeneity of subsurface materials, represented by a contrast between a higher‐hydraulic conductivity and lower‐hydraulic conductivity material normal to the direction of the applied electric field exerts the following controls on nitrate migration: (1) a spatial change in nitrate migration rate due to changes in effective ionic mobility and subsequent accumulation of nitrate at the interface between these materials; and (2) a spatial change in the voltage gradient distribution across the hydraulic conductivity contrast, due to the inverse relationship with effective ionic mobility. These factors will contribute to higher mass transport of nitrate through low hydraulic conductivity zones in heterogeneous porous media, relative to homogeneous host materials. Overall electrokinetic migration of amendments such as nitrate can be increased in heterogeneous granular porous media to enhance the in situ bioremediation of organic contaminants present in low hydraulic conductivity zones. 相似文献
9.
Hydraulic/partitioning tracer tomography (HPTT) was recently developed by Yeh and Zhu [Yeh T-CJ, Zhu J. Hydraulic/partitioning tracer tomography for characterization of dense nonaqueous phase liquid source zones, Water Resour Res 2007;43:W06435. doi:10.1029/2006WR004877.] for estimating spatial distribution of dense nonaqueous phase liquids (DNAPLs) in the subsurface. Since discrete tracer concentration data are directly utilized for the estimation of DNAPLs, this approach solves the hyperbolic convection–dispersion equation. Solution to the convection–dispersion equation however demands fine temporal and spatial discretization, resulting in high computational cost for an HPTT analysis. In this work, we use temporal moments of tracer breakthrough curves instead of discrete concentration data to estimate DNAPL distribution. This approach solves time independent partial differential equations of the temporal moments, and therefore avoids solving the convection–dispersion equation using a time marching scheme, resulting in a dramatic reduction of computational cost. To reduce numerical oscillations associated with convection dominated transport problems such as in inter-well tracer tests, the approach uses a finite element solver adopting the streamline upwind Petrov–Galerkin method to calculate moments and sensitivities. We test the temporal moment approach through numerical simulations. Comparing the computational costs between utilizing moments and discrete concentrations, we find that temporal moments significantly reduce the computation time. We also find that tracer moment data collected through a tomographic survey alone are able to yield reasonable estimates of hydraulic conductivity, as indicated by a correlation of 0.588 between estimated and true hydraulic conductivity fields in the synthetic case study. 相似文献
10.
This pore-scale modeling study in saturated porous media shows that compound-specific effects are important not only at steady-state and for the lateral displacement of solutes with different diffusivities but also for transient transport and solute breakthrough. We performed flow and transport simulations in two-dimensional pore-scale domains with different arrangement of the solid grains leading to distinct characteristics of flow variability and connectivity, representing mildly and highly heterogeneous porous media, respectively. The results obtained for a range of average velocities representative of groundwater flow (0.1–10 m/day), show significant effects of aqueous diffusion on solute breakthrough curves. However, the magnitude of such effects can be masked by the flux-averaging approach used to measure solute breakthrough and can hinder the correct interpretation of the true dilution of different solutes. We propose, as a metric of mixing, a transient flux-related dilution index that allows quantifying the evolution of solute dilution at a given position along the main flow direction. For the different solute transport scenarios we obtained dilution breakthrough curves that complement and add important information to traditional solute breakthrough curves. Such dilution breakthrough curves allow capturing the compound-specific mixing of the different solutes and provide useful insights on the interplay between advective and diffusive processes, mass transfer limitations, and incomplete mixing in the heterogeneous pore-scale domains. The quantification of dilution for conservative solutes is in good agreement with the outcomes of mixing-controlled reactive transport simulations, in which the mass and concentration breakthrough curves of the product of an instantaneous transformation of two initially segregated reactants were used as measures of reactive mixing. 相似文献
11.
Non-local stochastic moment equations are used successfully to analyze groundwater flow in randomly heterogeneous media. Here we present a moment equations-based approach to quantify the uncertainty associated with the estimation of well catchments. Our approach is based on the development of a complete second order formalism which allows obtaining the first statistical moments of the trajectories of conservative solute particles advected in a generally non-uniform groundwater flow. Approximate equations of moments of particles’ trajectories are then derived on the basis of a second order expansion in terms of the standard deviation of the aquifer log hydraulic conductivity. Analytical expressions are then obtained for the predictors of locations of mean stagnation points, together with their associated uncertainties. We implement our approach on heterogeneous media in bounded two-dimensional domains, with and without including the effect of conditioning on hydraulic conductivity information. The impact of domain size, boundary conditions, heterogeneity and non-stationarity of hydraulic conductivity on the prediction of a well catchment is explored. The results are compared against Monte Carlo simulations and semi-analytical solutions available in the literature. The methodology is applicable to both infinite and bounded domains and is free of distributional assumptions (and so applies to both Gaussian and non-Gaussian log hydraulic conductivity fields) and formally includes the effect of conditioning on available information. 相似文献
12.
Leanna Woods Pan Robert L. Siegrist Michelle Crimi 《Ground Water Monitoring & Remediation》2012,32(2):96-105
In situ remediation technologies have the potential to alter subsurface properties such as natural organic matter (NOM) content or character, which could affect the organic carbon‐water partitioning behavior of chlorinated organic solvents, including dense nonaqueous phase liquids (DNAPLs). Laboratory experiments were completed to determine the nature and extent of changes in the partitioning behavior of trichloroethene (TCE) caused by in situ chemical oxidation or in situ surfactant flushing. Sandy porous media were obtained from the subsurface at a site in Orlando, Florida. Experiments were run using soil slurries in zero‐headspace reactors (ZHRs) following a factorial design to study the effects of porous media properties (sand vs. loamy sand with different total organic carbon [TOC] contents), TCE concentration (DNAPL presence or absence), and remediation agent type (potassium permanganate vs. activated sodium persulfate, Dowfax 8390 vs. Tween 80). Results revealed that the fraction of organic carbon (foc) of porous media after treatment by oxidants or surfactants was higher or lower relative to that in the untreated media controls. Isotherm experiments were run using the treated and control media to measure the distribution coefficient (Kd) of TCE. Organic carbon‐water partitioning coefficient values (Koc) calculated from the experimental data revealed that Koc values for TCE in the porous media were altered via treatment using oxidants and surfactants. This alteration can affect the validity of estimates of contaminant mass remaining after remediation. Thus, potential changes in partitioning behavior should be considered to help avoid decision errors when judging the effectiveness of an in situ remediation technology. 相似文献
13.
F.A. Radu N. Suciu J. HoffmannA. Vogel O. Kolditz C.-H. ParkS. Attinger 《Advances in water resources》2011,34(1):47-61
This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater. 相似文献
14.
DNAPL to LNAPL Transitions During Horizontal Cosolvent Flooding 总被引：1，自引：0，他引：1
Eberhard Roeder Ronald W. Falta Jr. Cindy M. Lee John T. Coates 《Ground Water Monitoring & Remediation》2001,21(1):77-88
Cosolvent flooding is a technology with the potential to remove nonaqueous phase liquid (NAPL) sources from the subsurface. It can be used to initiate separate phase mobilization, which allows removal of NAPL within very few pore volumes. Mobilization may result in a sinking DNAPL bank during horizontal flooding of NAPLs denser than water. Reversal of phase density difference between aqueous and DNAPL phases could potentially avoid this downward migration of mobilized DNAPLs. We achieved phase density difference reversal and made DNAPLs float using two components in the cosolvent flooding solution. A low-density cosolvent partitions preferentially into the DNAPL and swells it, which causes a reduction in density of the DNAPL and reversal of the density difference between the NAPL and aqueous phases. A highdensity additive that remains in the aqueous phase allows the cosolvent flooding solution overall to have a density greater than that of water and permits control of the flooding instability. This study focused on tert-butanol as the swelling cosolvent and tetrachloroethylene as the contaminant. In batch tests with sucrose and glycerol as dense additives, phase density difference reversal occurred. To investigate the applicability of phase density difference reversal as a remediation technology, horizontal column and sandbox experiments were performed. These experiments demonstrated the occurrence of phase density difference reversal and effective remediation in horizontal cosolvent floods. 相似文献
15.
《Ground Water Monitoring & Remediation》2001,21(2):71-88
Volume reduction and lowering of capillary pressure within a large DNAPL pool are utilized as objectives in the design of a large-scale dual phase recovery system at a chemical manufacturing facility in the United States. By reducing DNAPL pool height through mass removal, capillary pressure is lowered, resulting in a reduced potential for future vertical and horizontal mobilization of the chlorinated solvent DNAPL pool. The DNAPL pool extends over an approximately 200 m by 275 m area in low permeability fill deposits overlying a clay aquitard. A three-dimensional multiphase flow model was employed to arrive at a final design incorporating nine horizontal drains (total length 664 m) and a pulsed pumping system. The numerical model was calibrated to the results of a 42-day field pilot-test involving the removal of approximately 25,000 L of DNAPL from a single, 55 m long horizontal drain. Numerical simulation revealed that gravity drainage, as opposed to hydraulic gradients in the water phase, is the dominant recovery mechanism at this site. This stems from the relatively high density and the viscosity of the DNAPL, and the relatively low permeability of the formation deposits. The use of pulsed pumping is shown to reduce the volume of contaminated ground water recovered from the 9-drain system, without significant reduction of the total volume of DNAPL recovered. 相似文献
16.
17.
Haley A. Schneider W. Andrew Jackson Ph.D. P.E. Paul B. Hatzinger Ph.D Charles E. Schaefer Ph.D 《Ground Water Monitoring & Remediation》2020,40(4):27-43
A direct-drive high-resolution passive profiler (HRPP) was developed to quantify and delineate concentrations of chlorinated volatile organic compounds (CVOCs), geochemical indicators and CVOC-degrading microorganisms/genes, as well as to perform compound-specific stable isotope analysis (CSIA) of CVOCs and estimate interstitial velocity at <30-cm resolution. The profilers can be coupled together to provide a continuous sample interval and advanced to depths up to approximately 9 m below-ground surface (bgs) within saturated media where direct-push techniques are feasible. The HRPP was field tested in a previous dense nonaqueous phase liquid (DNAPL) source zone at the former Naval Air Station in Alameda, CA. HRPP data sets were compared to the following traditional groundwater data sets: CVOC and anion concentrations in standard and multilevel monitoring well water samples, CVOC concentrations in soil core samples, qualitative contaminant profiles delineated with a membrane interface probe (MIP), microbial community and CSIA profiles from Bio-Traps® deployed in wells, groundwater velocity from passive flux meters (PFMs), lithologic profiles correlated with MIP electrical conductivity (EC), and velocity estimates based on permeability profiles measured with a Geoprobe hydraulic profiling tool (HPT). In some cases, the HRPP data were equivalent to traditional techniques and, in other cases, the HRPP data were more representative of local variability rather than bulk aquifer conditions. Overall the results support the use of the HRPP to provide high-resolution data on concentrations, velocity, and microbial activity in temporary direct-push deployments without well installation, providing a new tool to better assess source zones and contaminated groundwater plumes, even in low permeability media, and to increase the fidelity of site transport models. 相似文献
18.
Numerical simulations of variable-density flow and solute transport have been conducted to investigate dense plume migration for various configurations of 2D fracture networks. For orthogonal fractures, simulations demonstrate that dispersive mixing in fractures with small aperture does not stabilize vertical plume migration in fractures with large aperture. Simulations in non-orthogonal 2D fracture networks indicate that convection cells form and that they overlap both the porous matrix and fractures. Thus, transport rates in convection cells depend on matrix and fracture flow properties. A series of simulations in statistically equivalent networks of fractures with irregular orientation show that the migration of a dense plume is highly sensitive to the geometry of the network. If fractures in a random network are connected equidistantly to the solute source, few equidistantly distributed fractures favor density-driven transport. On the other hand, numerous fractures have a stabilizing effect, especially if diffusive transport rates are high. A sensitivity analysis for a network with few equidistantly distributed fractures shows that low fracture aperture, low matrix permeability and high matrix porosity impede density-driven transport because these parameters reduce groundwater flow velocities in both the matrix and the fractures. Enhanced molecular diffusion slows down density-driven transport because it favors solute diffusion from the fractures into the low-permeability porous matrix where groundwater velocities are smaller. For the configurations tested, variable-density flow and solute transport are most sensitive to the permeability and porosity of the matrix, which are properties that can be determined more accurately than the geometry and hydraulic properties of the fracture network, which have a smaller impact on density-driven transport. 相似文献
19.
The spatial distribution of hydraulic properties in the subsurface controls groundwater flow and solute transport. However, many approaches to modeling these distributions do not produce geologically realistic results and/or do not model the anisotropy of hydraulic conductivity caused by bedding structures in sedimentary deposits. We have developed a flexible object-based package for simulating hydraulic properties in the subsurface—the Hydrogeological Virtual Realities (HyVR) simulation package. This implements a hierarchical modeling framework that takes into account geological rules about stratigraphic bounding surfaces and the geometry of specific sedimentary structures to generate realistic aquifer models, including full hydraulic-conductivity tensors. The HyVR simulation package can create outputs suitable for standard groundwater modeling tools (e.g., MODFLOW), is written in Python, an open-source programming language, and is openly available at an online repository. This paper presents an overview of the underlying modeling principles and computational methods, as well as an example simulation based on the Macrodispersion Experiment site in Columbus, Mississippi. Our simulation package can currently simulate porous media that mimic geological conceptual models in fluvial depositional environments, and that include fine-scale heterogeneity in distributed hydraulic parameter fields. The simulation results allow qualitative geological conceptual models to be converted into digital subsurface models that can be used in quantitative numerical flow-and-transport simulations, with the aim of improving our understanding of the influence of geological realism on groundwater flow and solute transport. 相似文献
20.
本文以Biot提出的流体饱和多孔介质波动理论为基础, 建立了成层地基模型, 把地下水位以下的饱和土层用水饱和多孔介质模拟, 地下水位以上土层用气饱和多孔介质模拟. 通过研究入射平面简谐波在成层地基中的传播, 分析了地下水位变化对地震地面运动的影响. 结果表明: P波入射下, 当土体骨架相对刚度较小时, 地下水位变化对地表位移尤其是竖向地表位移幅值影响较大, 随着地下水位的下降, 竖向位移逐渐增加, 相对应的峰值频率逐渐减小; 当土体骨架相对刚度较大时, 地下水位变化对地面运动影响不大. 相似文献