首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
Wang F  Bright J 《Ground water》2004,42(5):760-766
The influence on solute transport of the small-scale spatial variation of aquifer hydraulic conductivity (K) was analyzed by comparing results from fine-grid (2 m by 2 m) simulations of a synthetic heterogeneous aquifer to those from coarse-grid (8 m by 4 m) simulations of an equivalent homogeneous aquifer. Realizations of the K field of the heterogeneous aquifer were generated, using the Monte Carlo approach, from a lognormal distribution with mean log K of 2 (K in m/d) and three levels of log K variance of 0.1, 0.5, and 1.0. Numerical simulation results show that the average standard deviation of point concentrations increased from 1.21 to 5.78 when the value of log K variance was increased from 0.1 to 1.0. The average discrepancy between modeled concentrations (obtained from a coarse-grid deterministic numerical simulation) and the actual mean point concentrations (obtained from fine-grid Monte Carlo numerical simulations) increased from 0.91 to 4.23 with the increase in log K variance. The results from this study illustrate the uncertainty in predictions from contaminant transport models due to their inability to simulate the effects of heterogeneities at scales smaller than the model grid.  相似文献   

2.
Monte Carlo simulations are conducted to evaluate microbial-mediated contaminant reactions in an aquifer comprised of spatially variable microbial biomass concentrations, aquifer hydraulic conductivities, and initial electron donor/acceptor concentrations. A finite element simulation model is used that incorporates advection, dispersion, and Monod kinetic expressions to describe biological processes. Comparisons between Monte Carlo simulations of heterogeneous systems and simulations using homogeneous formulation of the same two-dimensional transport problem are presented. For the assumed set of parameters, physical aquifer heterogeneity is found to have a minor effect on the mass of contaminant biodegraded/transformed when compared to a homogeneous system; however, it noticeably changes the dispersion, skewness, and peakness of contaminant concentration distributions. Similarly, for low microbial growth rate, given favorable microbial growth characteristics, biological heterogeneity has minor effect on the mass of contaminant biodegraded/transformed when compared to a homogeneous system. On the other hand, when higher effective growth rates are assumed, biological heterogeneity and spatial heterogeneities in essential electron donor/acceptors reduce the efficiency of biotic contaminant reactions; consequently, model simulations derived from heterogeneous biomass distributions predict remediation time scales that are longer than those simulated for homogeneous systems. When correlations between physical aquifer and biological heterogeneities are considered, the assumed correlation affects predicted mean and variance of contaminant concentration and biomass distributions. For example, an assumed negative correlation between hydraulic conductivity and the initial biomass distribution produces a plume where less efficient biotic contaminant reactions occur at the leading edge of the plume; this is consistent with less degradation/transformation occurring over regions of higher groundwater velocities. However, the presence and absence of these correlations do not appear to affect the efficiency of microbial-mediated contaminant attenuation.  相似文献   

3.
Flow and transport models in heterogeneous geological formations are usually large-scale with excessive computational complexity and uncertain characteristics. Uncertainty quantification for predicting subsurface flow and transport often entails utilizing a numerical Monte Carlo framework, which repeatedly simulates the model according to a random field parameter representing hydrogeological characteristics of the aquifer. The physical resolution (e.g. spatial grid resolution) for the simulation is customarily chosen based on recommendations in the literature, independent of the number of Monte Carlo realizations. This practice may lead to either excessive computational burden or inaccurate solutions. We develop an optimization-based methodology that considers the trade-off between the following conflicting objectives: time associated with computational costs, statistical convergence of the model prediction and physical errors corresponding to numerical grid resolution. Computational resources are allocated by considering the overall error based on a joint statistical–numerical analysis and optimizing the error model subject to a given computational constraint. The derived expression for the overall error explicitly takes into account the joint dependence between the discretization error of the physical space and the statistical error associated with Monte Carlo realizations. The performance of the framework is tested against computationally extensive simulations of flow and transport in spatially heterogeneous aquifers. Results show that modelers can achieve optimum physical and statistical resolutions while keeping a minimum error for a given computational time. The physical and statistical resolutions obtained through our analysis yield lower computational costs when compared to the results obtained with prevalent recommendations in the literature. Lastly, we highlight the significance of the geometrical characteristics of the contaminant source zone on the optimum physical and statistical resolutions.  相似文献   

4.
Today, in different countries, there exist sites with contaminated groundwater formed as a result of inappropriate handling or disposal of hazardous materials or wastes. Numerical modeling of such sites is an important tool for a correct prediction of contamination plume spreading and an assessment of environmental risks associated with the site. Many uncertainties are associated with a part of the parameters and the initial conditions of such environmental numerical models. Statistical techniques are useful to deal with these uncertainties. This paper describes the methods of uncertainty propagation and global sensitivity analysis that are applied to a numerical model of radionuclide migration in a sandy aquifer in the area of the RRC “Kurchatov Institute” radwaste disposal site in Moscow, Russia. We consider 20 uncertain input parameters of the model and 20 output variables (contaminant concentration in the observation wells predicted by the model for the end of 2010). Monte Carlo simulations allow calculating uncertainty in the output values and analyzing the linearity and the monotony of the relations between input and output variables. For the non monotonic relations, sensitivity analyses are classically done with the Sobol sensitivity indices. The originality of this study is the use of modern surrogate models (called response surfaces), the boosting regression trees, constructed for each output variable, to calculate the Sobol indices by the Monte Carlo method. It is thus shown that the most influential parameters of the model are distribution coefficients and infiltration rate in the zone of strong pipe leaks on the site. Improvement of these parameters would considerably reduce the model prediction uncertainty.  相似文献   

5.
This study investigates stochastic optimization of dense nonaqueous phase liquid (DNAPL) remediation design at Dover Air Force Base Area 5 using emulsified vegetable oil (EVO) injection. The Stochastic Cost Optimization Toolkit (SCOToolkit) is used for the study, which couples semianalytical DNAPL source depletion and transport models with parameter estimation, error propagation, and stochastic optimization modules that can consider multiple sources and remediation strategies. Model parameters are calibrated to field data conditions on prior estimates of parameters and their uncertainty. Monte Carlo simulations are then performed to identify optimal remediation decisions that minimize the expected net present value (NPV) cleanup cost while maintaining concentrations at compliance wells under the maximum contaminant level (MCL). The results show that annual operating costs could be reduced by approximately 50% by implementing the identified optimal remediation strategy. We also show that recalibration and reoptimization after 50 years using additional monitoring data could lead to a further 60% reduction in annual operating cost increases the reliability of the proposed remediation actions.  相似文献   

6.
We perform global sensitivity analysis (GSA) through polynomial chaos expansion (PCE) on a contaminant transport model for the assessment of radionuclide concentration at a given control location in a heterogeneous aquifer, following a release from a near surface repository of radioactive waste. The aquifer hydraulic conductivity is modeled as a stationary stochastic process in space. We examine the uncertainty in the first two (ensemble) moments of the peak concentration, as a consequence of incomplete knowledge of (a) the parameters characterizing the variogram of hydraulic conductivity, (b) the partition coefficient associated with the migrating radionuclide, and (c) dispersivity parameters at the scale of interest. These quantities are treated as random variables and a variance-based GSA is performed in a numerical Monte Carlo framework. This entails solving groundwater flow and transport processes within an ensemble of hydraulic conductivity realizations generated upon sampling the space of the considered random variables. The Sobol indices are adopted as sensitivity measures to provide an estimate of the role of uncertain parameters on the (ensemble) target moments. Calculation of the indices is performed by employing PCE as a surrogate model of the migration process to reduce the computational burden. We show that the proposed methodology (a) allows identifying the influence of uncertain parameters on key statistical moments of the peak concentration (b) enables extending the number of Monte Carlo iterations to attain convergence of the (ensemble) target moments, and (c) leads to considerable saving of computational time while keeping acceptable accuracy.  相似文献   

7.
A groundwater plume containing high concentrations of pharmaceutical compounds, mainly sulfonamides, barbiturates, and ethyl urethane, in addition to chlorinated ethenes and benzene was investigated. The contamination originating from a former pharmaceutical industry discharges into a multilayered aquifer system and a downgradient stream. In this study, geological and hydrogeological data were integrated into a numerical flow model to examine identified trends using statistical approaches, including principal component analysis and hierarchal cluster analysis. A joint interpretation of the groundwater flow paths and contaminant concentrations in the different compartments (i.e., groundwater and hyporheic zone) provided insight on the transport processes of the different contaminant plumes to the stream. The analysis of historical groundwater concentrations of pharmaceutical compounds at the site suggested these compounds are slowly degrading. The pharmaceutical compounds migrate in both a deep semiconfined aquifer, as well as in the shallow unconfined aquifer, and enter the stream along a 2-km stretch. This contrasted with the chlorinated ethenes, which mainly discharge to the stream as a focused plume from the unconfined aquifer. The integrated approach developed here, combining groundwater flow modeling and statistical analyses of the contaminant concentration data collected in groundwater and the hyporheic zone, lead to an improved understanding of the observed distribution of contaminants in the unconfined and semiconfined aquifers, and thus to their discharge to the stream. This approach is particularly relevant for large and long-lasting contaminant sources and plumes, such as abandoned landfills and industrial production sites, where field investigations may be very expensive.  相似文献   

8.
A finite element 2D Monte Carlo approach is used to evaluate the sensitivity of groundwater contaminant discharges to a Damkohler number ω and spatial variability in aquifer hydraulic conductivity, initial microbial biomass concentrations, and electron acceptor/donor concentrations. Bioattenuation is most sensitive to spatial variations in incipient biomass and critical electron donors/acceptors for ω ≥ 1 (i.e., when pore-water residence times are high compared to the time needed for microbial growth or contaminant attenuation). Under these conditions, critical reaction processes can become substrate-limited at multiple locations throughout the aquifer; which in turn increases expected contaminant discharges and their uncertainties at monitored transects. For ω ≤ 0.2, contaminant discharge is not sensitive to incipient biomass variations. Physical heterogeneities expedite plume arrival and delay departure at transects and in turn attenuate peak discharges but do not affect cumulative contaminant discharges. Physical heterogeneities do, however, induce transect mass discharge variances that are bimodal functions of time; the first peak beings consistently higher. A simple stream tube model is invoked to explain the occurrence of peaks in contaminant discharge variance.  相似文献   

9.
Purge and pump samples from screened wells reflect concentration averaging and contaminant redistribution by wellbore flow. These issues were assessed in a screened well at the Hanford Site by investigating the vertical profile of a technetium-99 plume in a conventional well under static and pumped conditions. Specific conductance and technetium-99 concentrations were well correlated, and this enabled measurement of specific conductance to be used as a surrogate for technetium-99 concentration. Time-series measurements were collected during purging from three specific conductance probes installed in the well at 1.2, 3.1, and 4.9 m below the static water level in a 7.7-m-deep screened well. The vertical contaminant profile adjacent to the well in the aquifer was calculated using the concentration profile in the well during pumping, the pumping flow rate, and a wellbore flow and mixing model. The plume was found to be stratified in the aquifer—the highest concentrations occurred adjacent to the upper part of the screened interval. The purge and pump sample concentrations were 41% to 58% of the calculated peak concentration in the aquifer. Plume stratification in the aquifer adjacent to the well screen became more pronounced as pumping continued. Extended pumping may have partially reversed the effect of contaminant redistribution in the aquifer by wellbore flow and allowed the stratification of the plume to be more observable. It was also found that the vertical profile of contamination in the well under static (i.e., nonpumping conditions) was not representative of the profile in the aquifer. Thus, passive or micropurge sampling techniques, which sample the wellbore water at different depths, would not yield results representative of the aquifer in this well.  相似文献   

10.
High concentrations of ammonia in a river can cause fish kills and harms to other aquatic organisms. A simple water quality model is needed to predict the probability of ammonia concentration violations as compared to the US Environmental Protection Agencys ammonia criteria. A spreadsheet with Random Monte Carlo (RMC) simulations to model ammonia concentrations at the mixing point (between a river and the effluent of a wastewater treatment plant) was developed with the use of Microsoft Excel and Crystal Ball add-in software. The model uses effluent and river ammonia, alkalinity, and total carbonate data to determine the probability density functions (PDFs) for the Monte Carlo simulations. Normal, lognormal, exponential and uniform probability distributions were tested using the Chi-square method and p-value associated with it to choose the best fit to the random data selected from the East Burlington wastewater treatment plant in North Carolina and the Clinch River in Tennessee. It is suggested that different options be tested with a minimum of three classes and a maximum of n/5 classes (n = number of data points) and the highest probability (p-value) for the PDF being tested be chosen. The results indicted that six violations to the EPA criterion for maximum concentration (CMC) were predicted when using 2000 RMC simulations and PDFs fitted to the available data, which violate the current criterion of no more than one violation over 3 years. All violations occur when the pH of the blend ranges from 8.0 to 9.0. No violations were found to the criteria of chronic concentration (CCC) using RMC.  相似文献   

11.
In risk analysis, a complete characterization of the concentration distribution is necessary to determine the probability of exceeding a threshold value. The most popular method for predicting concentration distribution is Monte Carlo simulation, which samples the cumulative distribution function with a large number of repeated operations. In this paper, we first review three most commonly used Monte Carlo (MC) techniques: the standard Monte Carlo, Latin Hypercube sampling, and Quasi Monte Carlo. The performance of these three MC approaches is investigated. We then apply stochastic collocation method (SCM) to risk assessment. Unlike the MC simulations, the SCM does not require a large number of simulations of flow and solute equations. In particular, the sparse grid collocation method and probabilistic collocation method are employed to represent the concentration in terms of polynomials and unknown coefficients. The sparse grid collocation method takes advantage of Lagrange interpolation polynomials while the probabilistic collocation method relies on polynomials chaos expansions. In both methods, the stochastic equations are reduced to a system of decoupled equations, which can be solved with existing solvers and whose results are used to obtain the expansion coefficients. Then the cumulative distribution function is obtained by sampling the approximate polynomials. Our synthetic examples show that among the MC methods, the Quasi Monte Carlo gives the smallest variance for the predicted threshold probability due to its superior convergence property and that the stochastic collocation method is an accurate and efficient alternative to MC simulations.  相似文献   

12.
Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.  相似文献   

13.
Aquifer vulnerability to pesticide migration through till aquitards   总被引：1，自引：0，他引：1  
This study investigates the influence of key factors-mainly recharge rate and degradation half-life--on downward migration of the widely used pesticide mecoprop (MCPP) through a typical clayey till aquitard. The study uses the numerical model FRAC3Dvs, which is a three-dimensional discrete fracture/matrix diffusion (DFMD) numerical transport model. The model was calibrated with laboratory and field data from a site near Havdrup, Denmark, but the overall findings are expected to be relevant to many other sites in similar settings. Fracture flow and MCPP transport parameters for the model were obtained through calibration using well-characterized laboratory experiments with large (0.5 m diameter by 0.5 m high) undisturbed columns of the fractured till and a field experiment. A second level of upscaling and sensitivity analysis was then carried out using data on hydraulic head, fracture spacing, and water budget from the field site. The simulations of downward migration of MCPP show that MCPP concentration and mass flux into the underlying aquifer, and hence the aquifer vulnerability to this pesticide compound, is mainly dependent on the degradation rate of the pesticide, the overall aquitard water budget, and the ground water recharge rate into the aquifer. The influence of flow rate, matrix diffusion, and degradation rate are intertwined. This results in one to four orders of magnitude higher MCPP flux into the aquifer from aquifer recharge rates of 20 and 120 mm/yr, respectively, for no degradation and MCPP half-life of 0.5 yr. From a sensitivity analysis it was found that the range of MCPP flux into the aquifer varied less than one order of magnitude due to (1) changing fracture spacing from 1 to 10 m, or (2) preferential flow along inclined thin sand layers, which represent common conditions for the current and other settings of clayey till in Denmark and other glaciated areas in Europe and North America. The results indicate that for aquifers overlain by fractured clayey tills, the vulnerability to contamination with pesticides (pesticide flux into the aquifer) and other widespread agricultural contaminants is going to vary strongly in the watershed as a function of the distribution of aquitard water budget (flow rate) and aquitard redox environment (controlling contaminant degradation rates), even if the thickness of the till is relatively constant. DFMD modeling of cause-effect relationships within such systems has great potential to support decisions in planning, regulation, and contaminant remediation.  相似文献   

14.
A solute breakthrough curve measured during a two-well tracer test was successfully predicted in 1986 using specialized contaminant transport models. Water was injected into a confined, unconsolidated sand aquifer and pumped out 125 feet (38.3 m) away at the same steady rate. The injected water was spiked with bromide for over three days; the outflow concentration was monitored for a month. Based on previous tests, the horizontal hydraulic conductivity of the thick aquifer varied by a factor of seven among 12 layers. Assuming stratified flow with small dispersivities, two research groups accurately predicted breakthrough with three-dimensional (12-layer) models using curvilinear elements following the arc-shaped flowlines in this test.
Can contaminant transport models commonly used in industry, that use rectangular blocks, also reproduce this breakthrough curve? The two-well test was simulated with four MODFLOW-based models, MT3D (FD and HMOC options), MODFLOWT, MOC3D, and MODFLOW-SURFACT.
Using the same 12 layers and small dispersivity used in the successful 1986 simulations, these models fit almost as accurately as the models using curvilinear blocks. Subtle variations in the curves illustrate differences among the codes. Sensitivities of the results to number and size of grid blocks, number of layers, boundary conditions, and values of dispersivity and porosity are briefly presented. The fit between calculated and measured breakthrough curves degenerated as the number of layers and/or grid blocks decreased, reflecting a loss of model predictive power as the level of characterization lessened. Therefore, the breakthrough curve for most field sites can be predicted only qualitatively due to limited characterization of the hydrogeology and contaminant source strength.  相似文献   

15.
A transient axisymmetric saturated-unsaturated numerical flow model was coupled with a particle tracking model to investigate the movement of contaminants when a shallow unconfined aquifer is pumped at a constant rate. The particle tracking model keeps track of locations and masses of solutes in the aquifer, and the time of capture by the well. At the end of each time-step the flow model solves the Richard's equation for the hydraulic head distribution from which elemental velocities are calculated. Solutes are then displaced for a period equivalent to the time-step using both the magnitude and direction of the elemental velocities. Numerical experiments were performed to investigate effluent concentrations in wells with screens of different length and in different positions relative to zones of stratified contamination. At early times of pumping the effluent concentrations were similar to the concentrations adjacent to the well screen, but at late times, the concentrations approached the vertically averaged concentration in the aquifer. Time to attain the vertically averaged concentration was determined by the well geometry, initial location of the contaminant plume in relation to the well screen, and hydraulic properties of the aquifer. The results are consistent with the hydraulics of flow to a pumping well and of particular importance, they demonstrate that short-term pump tests could give erroneous design concentrations for pump-and-treat systems. The model provides a means of quantifying arrival times and mixing ratios. It could therefore provide a useful means of designing production wells in aquifers with stratified contamination and more efficient recovery systems for aquifer remediation.  相似文献   

16.
Bayer P  Finkel M 《Ground water》2006,44(2):234-243
We investigate the performance of vertical hydraulic barriers in combination with extraction wells for the partial hydraulic isolation of contaminated aquifer areas. The potential advantage of such combinations compared to a conventional pump-and-treat system has already been demonstrated in a previous study. Here we extend the scope of the performance analysis to the impact of uncertainty in the regional flow direction as well as to highly heterogeneous aquifer transmissivity distributions. In addition, two new well-barrier scenarios are proposed and analyzed. The hydraulic efficiency of the scenarios is rated based on the expected (mean) reduction of the pumping rate that is required to achieve downgradient contaminant capture. The uncertain spatial distribution of aquifer transmissivity is considered by means of unconditioned Monte Carlo simulations. The significance of uncertain background flow conditions is incorporated by computing minimized pumping rates for deviations of the regional flow direction up to 30 degrees from a normative base case. The results give an answer on how pumping rates have to be changed for each barrier-well combination in order to achieve robust systems. It is exposed that in comparison to installing exclusively wells, the barrier-supported approach generally yields savings in the (average) pumping rate. The particular efficiency is shown to be highly dependent on the interaction of variance and integral scale of transmissivity distribution, well and barrier position, as well as direction of background flow.  相似文献   

17.
 Stochastic techniques, such as Monte Carlo experiments, are more and more frequently used for the study of flow and transport in heterogeneous aquifers. When the aquifer is composed of distinct hydrofacies, a common way to model heterogeneity is to first generate equally-possible hydrofacies fields, and then convert these hydrofacies fields into hydraulic conductivity (K) fields by assigning a single K value to each facies. This technique assumes relative homogeneity of K within each facies but may not be appropriate for the most conductive facies that often exhibits substantial variability. In this paper, we assessed the impacts of assigning multiple random K, rather than a uniform K value, to the highly conductive facies on the results of a flow and transport model. A set of fifty stochastic hydrofacies maps depicting an environment similar to the Snake River Plain aquifer (SRPA) in south-east Idaho were generated. Simulations demonstrated that a uniform K value, if carefully chosen, can reasonably reproduce the specific discharges and early particle arrival times produced by multiple K values. Yet, the results obtained with a uniform K value are dramatically less variable than those obtained with multiple K values. It is therefore concluded that stochastic simulations with uniform K assigned to the most conductive and variable facies do not necessarily portray the entire uncertainty in the analyses.  相似文献   

18.
At a study site in the midwestern United States, multiple-completion wells demonstrated that a vertical hydraulic gradient was responsible for the contamination pattern exhibited by chlorinated solvent plumes. The typical pattern consisted of little or no contamination in the upper portion of the aquifer with concentrations increasing with depth. When ground water contamination was discovered in an unexpected portion of the site, water level elevations and contaminant distribution data obtained from multiple-completion wells resulted in identification of the source location. The well eventually determined to be located in the source area displayed contaminant levels much higher in the upper zone of the aquifer — the opposite contamination pattern of other on-site wells. Such results indicated that the spill had occurred near this location and that solvent residing along the capillary fringe was continuing to contaminate the aquifer.  相似文献   

19.
Most established methods to characterize aquifer structure and hydraulic conductivities of hydrostratigraphical units are not capable of delivering sufficient information in the spatial resolution that is desired for sophisticated numerical contaminant transport modeling and adapted remediation design. With hydraulic investigation methods based on the direct-push (DP) technology such as DP slug tests, DP injection logging, and the hydraulic profiling tool, it is possible to rapidly delineate hydrogeological structures and estimate their hydraulic conductivity in shallow unconsolidated aquifers without the need for wells. A combined application of these tools was used for the investigation of a contaminated German refinery site and for the setup of hydraulic aquifer models. The quality of DP investigation and the models was evaluated by comparisons of tracer transport simulations using these models and measured breakthroughs of two natural gradient tracer tests. Model scenarios considering the information of all tools together showed good reproduction of the measured breakthroughs, indicating the suitability of the approach and a minor impact of potential technical limitations. Using the DP slug tests alone yielded significantly higher deviations for the determined hydraulic conductivities compared to considering two or three of the tools. Realistic aquifer models developed on basis of such combined DP investigation approaches can help optimize remediation concepts or identify flow regimes for aquifers with a complex structure.  相似文献   

20.
A numerical study was conducted to investigate the influence of tides on the fate of terrestrially derived BTEX discharging through an unconfined aquifer to coastal waters. Previous studies have revealed that tide-induced seawater circulations create an active salt–freshwater mixing zone in the near-shore aquifer and alter the specific subsurface pathway for contaminants discharging to the coastal environment. Here the coupled density-dependent flow and multi-species reactive transport code PHWAT was used to examine the impact of these tidal effects on the aerobic biodegradation of BTEX released in a coastal aquifer and its subsequent loading to coastal waters. Simulations indicated that tides significantly enhance BTEX attenuation in the near-shore aquifer. They also reduce the rate of chemical transfer from the aquifer to the ocean and exit concentrations at the beach face. For the base case consisting of toluene transport and biodegradation, 79% of toluene initially released in the aquifer was attenuated prior to discharge with tides present, compared to only 1.8% for the non-tidal case. The magnitude of tidal forcing relative to the fresh groundwater flow rate was shown to influence significantly the extent of biodegradation as it controls the intensity of salt–freshwater mixing, period of exposure of the contaminant to the mixing zone and rate of oxygen delivery to the aquifer. The oxygen available for biodegradation also depends on the rate at which oxygen is consumed by natural processes such as organic matter decomposition. While simulations conducted with heterogeneous conductivity fields highlighted the uncertainties associated with predicting contaminant loadings, the study revealed overall that BTEX may undergo significant attenuation in tidally influenced aquifers prior to discharge.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号