共查询到20条相似文献，搜索用时 31 毫秒
1.
The main purpose of this paper was to compare three approaches for predicting solute transport. The approaches include: (1) an effective parameter/macrodispersion approach (Gelhar and Axness 1983); (2) a heterogeneous approach using ordinary kriging based on core samples; and (3) a heterogeneous approach based on hydraulic tomography. We conducted our comparison in a heterogeneous sandbox aquifer. The aquifer was first characterized by taking 48 core samples to obtain local-scale hydraulic conductivity (K). The spatial statistics of these K values were then used to calculate the effective parameters. These K values and their statistics were also used for kriging to obtain a heterogeneous K field. In parallel, we performed a hydraulic tomography survey using hydraulic tests conducted in a dipole fashion with the drawdown data analyzed using the sequential successive linear estimator code (Yeh and Liu 2000) to obtain a K distribution (or K tomogram). The effective parameters and the heterogeneous K fields from kriging and hydraulic tomography were used in forward simulations of a dipole conservative tracer test. The simulated and observed breakthrough curves and their temporal moments were compared. Results show an improvement in predictions of drawdown behavior and tracer transport when the K tomogram from hydraulic tomography was used. This suggests that the high-resolution prediction of solute transport is possible without collecting a large number of small-scale samples to estimate flow and transport properties that are costly to obtain at the field scale. 相似文献
2.
Mark L. Kram Arturo A. Keller Joseph Rossabi Lome G. Everett 《Ground Water Monitoring & Remediation》2002,22(1):46-61
Contamination from the use of chlorinated solvents, often classified as dense nonaqueous phase liquids (DNAPLs) when in an undissolved state, pose environmental threats to ground water resources worldwide. DNAPL site characterization method performance comparisons are presented in a companion paper (Kram et al. 2001). This study compares the costs for implementing various characterization approaches using synthetic unit model scenarios (UMSs), each with particular physical characteristics. Unit costs and assumptions related to labor, equipment, and consumables are applied to determine costs associated with each approach for various UMSs. In general, the direct-push sensor systems provide cost-effective characterization information in soils that are penetrable with relatively shallow (less than 10 to 15 m) water tables. For sites with impenetrable lithology using direct-push techniques, the Ribbon NAPL Sampler Flexible Liner Underground Technologies Everting (FLUTe) membrane appears to be the most cost-effective approach. For all scenarios studied, partitioning interwell tracer tests (PITTs) are the most expensive approach due to the extensive pre-and post-PITT requirements. However, the PITT is capable of providing useful additional information, such as approximate DNAPL saturation, which is not generally available from any of the other approaches included in this comparison. 相似文献
3.
It has been known for many years that dispersivity increases with solute travel distance in a subsurface environment. The increase of dispersivity with solute travel distance results from the significant variation of hydraulic properties of heterogeneous media and was identified in the literature as scale-dependent dispersion. This study presents an analytical solution for describing two-dimensional non-axisymmetrical solute transport in a radially convergent flow tracer test with scale-dependent dispersion. The power series technique coupling with the Laplace and finite Fourier cosine transform has been applied to yield the analytical solution to the two-dimensional, scale-dependent advection–dispersion equation in cylindrical coordinates with variable-dependent coefficients. Comparison between the breakthrough curves of the power series solution and the numerical solutions shows excellent agreement at different observation points and for various ranges of scale-related transport parameters of interest. The developed power series solution facilitates fast prediction of the breakthrough curves at any observation point. 相似文献
4.
Hydraulic displacement is a mass removal technology suitable for stabilization of a dense, nonaqueous phase liquid (DNAPL) source zone, where stabilization is defined as reducing DNAPL saturations and reducing the risk of future pool mobilization. High resolution three-dimensional multiphase flow simulations incorporating a spatially correlated, heterogeneous porous medium illustrate that hydraulic displacement results in an increase in the amount of residual DNAPL present, which in turn results in increased solute concentrations in groundwater, an increase in the rate of DNAPL dissolution, and an increase in the solute mass flux. A higher percentage of DNAPL recovery is associated with higher initial DNAPL release volumes, lower density DNAPLs, more heterogeneous porous media, and increased drawdown of groundwater at extraction wells. The fact that higher rates of recovery are associated with more heterogeneous porous media stems from the fact that larger contrasts in permeability provide for a higher proportion of capillary barriers upon which DNAPL pooling and lateral migration can occur. Across all scenarios evaluated in this study, the ganglia-to-pool (GTP) ratio generally increased from approximately 0.1 to between approximately 0.3 and 0.7 depending on the type of DNAPL, the degree of heterogeneity, and the imposed hydraulic gradient. The volume of DNAPL recovered as a result of implementing hydraulic displacement ranged from between 9.4% and 45.2% of the initial release volume, with the largest percentage recovery associated with 1,1,1 trichloroethane, the least dense of the three DNAPLs considered. 相似文献
5.
6.
Modeling dispersion in homogeneous porous media with the convection–dispersion equation commonly requires computing effective transport coefficients. In this work, we investigate longitudinal and transverse dispersion coefficients arising from the method of volume averaging, for a variety of periodic, homogeneous porous media over a range of particle Péclet (Pep) numbers. Our objective is to validate the upscaled transverse dispersion coefficients and concentration profiles by comparison to experimental data reported in the literature, and to compare the upscaling approach to the more common approach of inverse modeling, which relies on fitting the dispersion coefficients to measured data. This work is unique in that the exact microscale geometry is available; thus, no simplifying assumptions regarding the geometry are required to predict the effective dispersion coefficients directly from theory. Transport of both an inert tracer and non-chemotactic bacteria is investigated for an experimental system that was designed to promote transverse dispersion. We highlight the occurrence of transverse dispersion coefficients that (1) depart from power-law behavior at relatively low Pep values and (2) are greater than their longitudinal counterparts for a specific range of Pep values. The upscaling theory provides values for the transverse dispersion coefficient that are within the 98% confidence interval of the values obtained from inverse modeling. The mean absolute error between experimental and upscaled concentration profiles was very similar to that between the experiments and inverse modeling. In all cases the mean absolute error did not exceed 12%. Overall, this work suggests that volume averaging can potentially be used as an alternative to inverse modeling for dispersion in homogeneous porous media. 相似文献
7.
When individual cells of a multiple-cell treatment wetland are hydraulically connected, the wetland has a cell-network structure. The hydraulic performance of treatment wetlands is often characterized using tracer residence time distributions (RTDs) measured between the wetland inlet and outlet, such that the wetland is considered as a single hydraulic unit, regardless of the extent of networking between individual internal cells. This work extends the single hydraulic unit approach to enable the specification of moments and RTD parameters for individual cells, or clusters of cells, within the cell-network based on inert tracer tests with injection only at the network inlet. Hydraulic performance is quantified in terms of hydraulic efficiency and travel time dimensionless variance using both the method of moments and RTD modeling. Cell-network analysis was applied to a case study from the Orlando Easterly Wetland (OEW), demonstrating the improvement in hydraulic performance of individual wetland cells following wetland restoration activities. Furthermore, cell-network analysis indicated that the location of water quality sampling station locations within the cell network can significantly affect the accuracy of pollutant removal effectiveness estimation when the individual sample station RTD does not represent the hydraulic unit RTD. At the OEW, it was determined that historical nutrient removal effectiveness estimation may be underestimated for one area and overestimated for another, and recommendations were provided for sample station locations to minimize future performance estimation errors. 相似文献
8.
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. 相似文献
9.
Ross [Ross PJ. Modeling soil water and solute transport – fast, simplified numerical solutions. Agron J 2003;95:1352–61] developed a fast, simplified method for solving Richards’ equation. This non-iterative 1D approach, using Brooks and Corey [Brooks RH, Corey AT. Hydraulic properties of porous media. Hydrol. papers, Colorado St. Univ., Fort Collins; 1964] hydraulic functions, allows a significant reduction in computing time while maintaining the accuracy of the results. The first aim of this work is to confirm these results in a more extensive set of problems, including those that would lead to serious numerical difficulties for the standard numerical method. The second aim is to validate a generalisation of the Ross method to other mathematical representations of hydraulic functions. 相似文献
10.
The success of transient storage (TS) modeling for natural streams depends, in part, on the ability to describe the dispersion process accurately. Evidence based on stream tracer data shows that solute transport processes often do not follow the classical second-order dispersion model (e.g., early breakthrough and faster than Fickian travel times were observed). While models based on space-fractional dispersion are a promising alternative, different definitions of fractional derivatives exist in the literature. Unlike integer-order derivatives, fractional derivatives represent convolutions of concentration with long-range spatial correlation and numerical approximations can produce dense matrices. Therefore issues of both accuracy and computational efficiency need to be examined to successfully identify model parameters for natural streams. In this paper, we first compare the performance of several numerical approaches for solving the space-fractional dispersion equation. We examine three different numerical approaches to approximate the space-fractional derivatives including: (a) a fully-implicit scheme based on the shifted Grünwald–Letnikov (GL) approximation (b) a three-point implicit representation based on the GL formula and (c) a three-point implicit scheme based on mass conservation and the Caputo definition of the fractional derivative. We then use an operator-splitting technique to evaluate a TS model based on space-fractional dispersion (the FSTS model) and test the model against analytical solutions and stream tracer data. A sequence acceleration method (Richardson extrapolation) significantly improves the performance of all schemes examined. Results indicate that the fully-implicit GL method with Richardson extrapolation produces the most accurate solutions while the three-point implicit GL scheme has a stringent time-step restriction to produce acceptable solutions. The three-point implicit scheme based on the Caputo derivative produces accurate solutions in a fraction of the time taken by the fully-implicit GL method and represents the best trade-off between accuracy and computational efficiency for practical applications. The scheme is suitable for parameter estimation and is used to successfully describe tracer data in a natural stream. 相似文献
11.
A tracer test in a carbonate aquifer is analyzed using the method of moments and two analytical advection-dispersion models (ADMs) as well as a numerical model. The numerical model is a coupled continuum-pipe flow and transport model that accounts for two different flow components in karstified carbonate aquifers, i.e., rapid and often turbulent conduit flow and Darcian flow in the fissured porous rock. All techniques employed provide reasonable fits to the tracer breakthrough curve (TBC) measured at a spring. The resulting parameter estimates are compared to investigate how each conceptual model of flow and transport processes that forms the basis of the analyses affects the interpretation of the tracer test. Numerical modeling results suggest that the method of moments and the analytical ADMs tend to overestimate the conduit volume because part of the water discharged at the spring is wrongly attributed to the conduit system if flow in the fissured porous rock is ignored. In addition, numerical modeling suggests that mixing of the two flow components accounts for part of the dispersion apparent in the measured TBC, while the remaining part can be attributed to Taylor dispersion. These processes, however, cannot reasonably explain the tail of the TBC. Instead, retention in immobile-fluid regions as included in a nonequilibrium ADM provides a possible explanation. 相似文献
12.
To the extent that sea surface temperature and colors can be considered passive tracers, their motions can be tracked to estimate the current velocities, or a conservation equation can be invoked to relate their temporal variations to the velocities. We investigate the latter, the so-called tracer inversion problem, with a particular focus on (1) the conditions under which the problem can be rendered over-determined for least squares solutions, (2) the possibility of using the tracer conservation equation within the “velocity projection” framework to estimate subsurface current profiles in shallow coastal waters, and (3) the accuracy of the tracer inversion calculation in terms of the data resolution and noise. The velocity projection framework refers to relating surface motion, either measured directly or made visible by tracers, to the subsurface current motion through the equations of motion. The accuracy of the tracer inversion calculation is quantified in terms of the spatial and temporal resolution of the tracer distribution. In the presence of irreducible tracer noise, the accuracy of the inversion rapidly degrades, and it is shown that the inversion with velocity projection can help improve accuracy. The tracer inversion method developed in this study is applied to the satellite sea surface temperature data, and the velocity result is compared to the velocity measurements made with the shore-based HF Coastal Current Radar. The potential of improving the velocity estimation with the present approach is indicated. 相似文献
13.
《Advances in water resources》2004,27(4):429-444
Partitioning interwell tracer tests (PITTs) are a relatively new technique for measuring the amount of nonaqueous phase liquid (NAPL) within saturated porous media. In this work we examined the influence of mass transfer limitations on the accuracy of measured NAPL from PITTs. Two mathematical models were used along with laboratory column experiments to explore the influence of tracer partition coefficient, tracer detection limit, and injected tracer mass on NAPL measurements. When dimensionless mass transfer coefficients were small, NAPL measurement errors decreased with decreasing tracer partition coefficient, decreasing tracer detection limit, and increasing injected tracer mass. Extrapolating breakthrough curves exponentially reduced but did not eliminate systematic errors in NAPL measurement. Although transport in a single stream tube was used in the mathematical models and laboratory experiments, the results from this simplified domain were supported by data taken from a three-dimensional computational experiment, where the NAPL resided as large pool. Based on these results, we suggest guidelines for interpreting tracer breakthrough data to ascertain the importance of mass transfer limitations on NAPL measurements. 相似文献
14.
Macropores are subsurface connected void spaces caused by processes such as fracture of soils, micro‐erosion, and fauna burrows. They are common near streams (e.g. hyporheic and riparian zones) and may act as preferential flow paths between surface and groundwaters, affecting hydrologic and biogeochemical processes. We tested the hydrologic function of macropores by constructing an artificial macropore within the saturated zone of a meander bend (open macropore, ‘OM’) and later filling its upstream end (partially filled macropore, ‘PFM’). For each treatment, we injected saline tracer at an upgradient monitoring well within the meander and monitored downgradient hydraulics and tracer transport. Pressure transducers in monitoring wells indicated hydraulic gradients within the meander were 32% higher perpendicular to and 6% higher parallel to the macropore for the OM than for the PFM. Additionally, hydraulic conductivities measured via falling head tests were 29 to 550 times higher along the macropore than in nearby sediment. We used electrical conductivity probes in wells and electrical resistivity imaging to track solute transport. Transport velocities through the meander were on average 9 and 21% higher (per temporal moment analysis and observed tracer peak, respectively) for the OM than for the PFM. Furthermore, temporal moments of tracer breakthrough analysis indicated downgradient longitudinal dispersion and breakthrough tracer curve tailing were on average 234% and 182% higher for the OM, respectively. This suggests the OM enabled solute transport at overall shorter timescales than the matrix but also increased tailing. Our results demonstrate the importance of macropores to meander bend hydrology and solute transport. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
15.
The inherent heterogeneity of geological media often results in anomalous dispersion for solute transport through them, and how to model it has been an interest over the past few decades. One promising approach that has been increasingly used to simulate the anomalous transport in surface and subsurface water is the fractional advection–dispersion equation (FADE), derived as a special case of the more general continuous time random walk or the stochastic continuum model. In FADE, the dispersion is not local and the solutes have appreciable probability to move long distances, and thus reach the boundary faster than predicted by the classical advection–dispersion equation (ADE). How to deal with different boundaries associated with FADE and their consequent impact is an issue that has not been thoroughly explored. In this paper we address this by taking one-dimensional solute movement in soil columns as an example. We show that the commonly used FADE with its fractional derivatives defined by the Riemann–Liouville definition is problematic and could result in unphysical results for solute transport in bounded domains; a modified method with the fractional dispersive flux defined by the Caputo derivatives is presented to overcome this problem. A finite volume approach is given to numerically solve the modified FADE and its associated boundaries. With the numerical model, we analyse the inlet-boundary treatment in displacement experiments in soil columns, and find that, as in ADE, treating the inlet as a prescribed concentration boundary gives rise to mass-balance errors and such errors could be more significant in FADE because of its non-local dispersion. We also discuss a less-documented but important issue in hydrology: how to treat the upstream boundary in analysing the lateral movement of tracer in an aquifer when the tracer is injected as a pulse. It is shown that the use of an infinite domain, as commonly assumed in literature, leads to unphysical backward dispersion, which has a significant impact on data interpretation. To avoid this, the upstream boundary should be flux-prescribed and located at the upstream edge of the injecting point. We apply the model to simulate the movement of Cl− in a tracer experiment conducted in a saturated hillslope, and analyse in details the significance of upstream-boundary treatments in parameter estimation. 相似文献
16.
The method of temporal moments is an efficient approach for analyzing breakthrough curves (BTCs). By matching the moments of the BTCs computed through parametric transfer-function models or one-dimensional transport models to those of the data, one can estimate the parameters characterizing the transfer function or apparent transport parameters. The classical method of moments presumes infinite duration. However, the measurement of BTCs is usually terminated prematurely, before the concentration has reached zero. Unless this truncation of the BTCs has been taken into account, the estimates of the parameters may be in error. Truncated measured BTCs are sometimes extrapolated assuming exponential decay. In this study, we use the concept of moments of the truncated impulse–response function [Jawitz JW. Moments of truncated continuous univariate distributions. Adv Water Res 2004;27:269–81] in the analysis of truncated BTCs corresponding to the commonly encountered step and step-pulse injection modes. The method is straightforward, based on the relation, which we derive, between truncated moments of the impulse–response function and the measured BTC. It is practical to apply and does not require the extrapolation of the measured BTC. The method is also accurate. In a numerical study we discuss how short a step-pulse injection may be so that we can approximate it as instantaneous. Finally, we apply the method to the analysis of a field-scale tracer test. 相似文献
17.
Natural Gradient Tracer Test to Evaluate Natural Attenuation of MTBE Under Anaerobic Conditions 总被引：1，自引：0，他引：1
A natural gradient tracer test using perdeuterated MTBE was conducted in an anaerobic aquifer to determine the relative importance of dispersion and degradation in reducing MTBE concentrations in ground water. Preliminary ground water chemistry and hydraulic conductivity data were used to place the tracer within an existing dissolved MTBE plume at Port Hueneme, California. Following one year of transport, the tracer plume was characterized in detail.
Longitudinal dispersion was identified as the dominant mechanism for lowering the perdeuterated MTBE concentrations. The method of moments was used to determine the longitudinal and lateral dispersion coefficients (0.85 m2 /day and 0.08 m2 /day, respectively). A mass-balance analysis, carried out after one year of transport, accounted for 110% of the injected mass and indicated that no significant mass loss occurred. The plume structure created by zones of higher and lower hydraulic conductivity at the site was complex, consisting of several localized areas of high tracer concentration in a lower concentration plume. This is important because the aquifer has generally been characterized as exhibiting fairly minor heterogeneity. In addition, the tracer plume followed a curved flowpath that deviated from the more macroscopic direction of ground water flow inferred from local ground water elevation measurements and the behavior of the existing plume. Understanding the mass balance, plume structure, curvature of the tracer plume, and consequently natural attenuation behavior required the detailed sampling approach employed in this study. These data imply that a detailed understanding of site hydrogeology and an extensive sampling network may be critical for the correct interpretation of monitored natural attenuation of MTBE. 相似文献
Longitudinal dispersion was identified as the dominant mechanism for lowering the perdeuterated MTBE concentrations. The method of moments was used to determine the longitudinal and lateral dispersion coefficients (0.85 m
18.
This paper describes the use of single-well tracer dilution techniques to resolve the rate of light nonaqueous phase liquid (LNAPL) flow through wells and the adjacent geologic formation. Laboratory studies are presented in which a fluorescing tracer is added to LNAPL in wells. An in-well mixer keeps the tracer well mixed in the LNAPL. Tracer concentrations in LNAPL are measured through time using a fiber optic cable and a spectrometer. Results indicate that the rate of tracer depletion is proportional to the rate of LNAPL flow through the well and the adjacent formation. Tracer dilution methods are demonstrated for vertically averaged LNAPL Darcy velocities of 0.00048 to 0.11 m/d and LNAPL thicknesses of 9 to 24 cm. Over the range of conditions studied, results agree closely with steady-state LNAPL flow rates imposed by pumping. A key parameter for estimating LNAPL flow rates in the formation is the flow convergence factor alpha. Measured convergence factors for 0.030-inch wire wrap, 0.030-inch-slotted polyvinyl chloride (PVC), and 0.010-inch-slotted PVC are 1.7, 0.91, and 0.79, respectively. In addition, methods for using tracer dilution data to determine formation transmissivity to LNAPL are presented. Results suggest that single-well tracer dilution techniques are a viable approach for measuring in situ LNAPL flow and formation transmissivity to LNAPL. 相似文献
19.
The internal tracer method for estimating contaminant degradation rates separates the attenuation effects not associated with degradation by using a codisposed recalcitrant internal tracer to normalize the degrading contaminant concentration. The remaining attenuation between the internal tracer and degrading contaminant is attributed to degradation and the degradation rate half-life is estimated from the first-order decay equation. An analytical solution of the advection- dispersion equation was used to evaluate flow-and-transport conditions that could result in incorrect estimates of contaminant degradation rate constants. Flow-and-transport characteristics that result in overestimating degradation rates were of particular interest because the internal tracer method often used to demonstrate natural attenuation can achieve remedial objectives. The analytical solution was also used to estimate the magnitude of error associated with using the internal tracer method at an example site and to explain different degradation rates estimated using tracers with different decay rate constants. 相似文献