地下水中轻质有机污染物(LNAPL)透镜体研究

在二维砂槽模型中模拟了轻质油在均匀多孔介质地下水非饱和区中的运移过程。模拟结果表明,地下水毛细区是轻质油污染的重点区,除了ＬＮＡＰＬ的残留以外,进入地下水饱和的ＬＮＡＰＬ终将被地下水顶托回到毛细区中,毛细区以上的约大多民将在重力作用下进入毛细区中,试验中观察到达稳定状态时ＬＮＡＰＬ透镜体的上边缘略微高出毛细区。利用多孔介质毛细管模型,建立了利用界面张力、接触角、介质特征孔隙直径等物理量估算不同位置     

Source-zone characterization of a chlorinated-solvent contaminated Superfund site in Tucson, AZ

An extensive site-characterization project was conducted at a large chlorinated-solvent contaminated Superfund site in Tucson, AZ. The project consisted of several components, including traditional site-characterization activities, tracer tests, laboratory experiments conducted with core material collected from the site, and mathematical modeling. The primary focus of the work presented herein is the analysis of induced-gradient contaminant elution tests conducted in a source zone at the site, investigation of the potential occurrence of immiscible liquid in the saturated zone, characterization of the relationship between mass flux reduction and mass removal, and evaluation of the impact of source-zone management on site remediation. The results of the present study, along with those of prior work, indicate that immiscible liquid is likely present in the saturated zone at the site source zones. Extensive tailing and rebound was observed for the contaminant-elution tests, indicating nonideal transport and mass-transfer behavior. The elution data were analyzed with a source-zone-scale mathematical model, and the results indicated that nonideal immiscible-liquid dissolution was the primary cause of the observed behavior. The time-continuous relationship between mass flux reduction and mass removal associated with the plume-scale pump-and-treat operation exhibited an initial large drop in mass flux with minimal mass removed, followed by a period of minimal mass flux reduction and a second period of large reduction. This behavior reflects the impact of both source-zone and aqueous-plume mass removal dynamics. Ultimately, a greater than 90% reduction in mass flux was achieved for a mass removal of approximately 50%. The influence of source-zone management on site remediation was evaluated by conducting two predictive simulations, one for which the source zones were controlled and one for which they were not. A plume-scale model was used to simulate the composite contaminant concentrations associated with groundwater extracted with the pump-and-treat system, which were compared to measured data. The information generated from this study was used to enhance the site conceptual model, help optimize operation of the pump-and-treat system, and evaluate the utility of source-zone remediation.     

Effect of water-table fluctuation on dissolution and biodegradation of a multi-component, light nonaqueous-phase liquid

Light nonaqueous-phase liquids (LNAPLs) such as gasoline and diesel fuel are among the most common causes of soil and groundwater contamination. Dissolution and subsequent advective transport of LNAPL components can negatively impact water supplies, while biodegradation is thought to be an important sink for this class of contaminants. We present a laboratory investigation of the effect of a water-table fluctuation on dissolution and biodegradation of a multi-component LNAPL (85% hexadecane, 5% toluene, 5% ethylbenzene, and 5% 2-methylnapthalene on a molar basis) in a pair of similar model aquifers (80 cm x 50 cm x 3 cm), one of which was subjected to a water-table fluctuation. Water-table fluctuation resulted in LNAPL and air entrapment below the water table, an increase in the vertical extent of the LNAPL source zone (by factor 6.7), and an increase in the volume of water passing through the source zone (by factor ~18). Effluent concentrations of dissolved LNAPL components were substantially higher and those of dissolved nitrate lower in the model aquifer where a fluctuation had been induced. Thus, water-table fluctuation led to enhanced biodegradation activity (28.3 mmol of nitrate consumed compared to 16.3 mmol in the model without fluctuation) as well as enhanced dissolution of LNAPL components. Despite the increased biodegradation, fluctuation led to increased elution of dissolved LNAPL components from the system (by factors 10-20). Hence, water-table fluctuations in LNAPL-contaminated aquifers might be expected to result in increased exposure of downgradient receptors to LNAPL components. Accordingly, water-table fluctuations in contaminated aquifers are probably undesirable unless the LNAPL is of minimal solubility or the dissolved-phase plume is not expected to reach a receptor due to distance or the presence of some form of containment.     

Persistence of LNAPL sources: relationship between risk reduction and LNAPL recovery

Light nonaqueous phase liquids (LNAPLs), such as fuels, are the source of much soil and groundwater contamination. Though the mobility of LNAPLs is limited in many environments, dissolved-phase components, such as benzene, can produce groundwater plumes that are more mobile than the LNAPL source. In such a setting, it is commonly assumed that recovery of the LNAPL will result in a reduction in risk associated with the dissolved phase. This paper synthesizes several existing multiphase and chemical transport solutions into a single linked methodology that predicts concentrations of soluble constituents within and downgradient of LNAPL source zones from dissolution of those constituents into groundwater flowing through and below LNAPL sources. This approach has been applied to a variety of LNAPL spill conditions. For biodegradable compounds, these analyses show that the period of time where the dissolved-phase plume is expanding is very small compared to the duration of most LNAPL sources, and that the downgradient extent is generally less than about 100 m for BTEX compounds. Therefore, the risk to receptors, as measured by the maximum downgradient extent of dissolved-phase plume or the maximum concentration of these compounds at a downgradient receptor, is generally unrelated to the longevity of the LNAPL sources. The maximum downgradient extent of the dissolved-phase plume is determined almost entirely by the groundwater velocity and the biodegradation rate. These analyses further demonstrate that recovery of LNAPL by hydraulic methods is often ineffective at reducing risk. Except in coarse-grained soils or intermediate soils with significant LNAPL saturations, free-product recovery approaches do not result in significant reductions in the longevity of downgradient dissolved-phase contamination. Further, for biodegradable constituents, remediation does not result in a near-term decrease in the downgradient extent of contamination. Cleanup methods that act to change the composition of the LNAPL source are more effective at reducing the downgradient concentrations, particularly for fine-grained soils or when LNAPL saturations are low.     

A simple model which predicts some non-linear features between atmospheric sulphur and sulphur emissions

This paper develops a simple model and suggests a plausible chemico-physical mechanism for a non-linear response between atmospheric sulphur and sulphur emissions. It contains simplified representations of transport, deposition and conversion processes and uses a proxy in-cloud oxidant-limited reaction along a pathway connecting an emission source with a receptor site. Individual pathway responses to emissions show linear behaviour above a threshold. However, by averaging the values of SO2 at the receptor site from different pathways a continuous non-linear relationship is obtained. As emissions reduce, distant emission sources become less significant contributors of sulphur dioxide at a receptor site but their emissions are still counted in an emission inventory, leading to an apparent non-linearity. Sulphate is always found to contribute a signal to the receptor site total. This model goes someway to explaining a proposed 'crossover' between observed proportions of wet and dry deposited sulphur in the UK as emissions have been reduced.     

Assessing the impact of source-zone remediation efforts at the contaminant-plume scale through analysis of contaminant mass discharge

The long-term impact of source-zone remediation efforts was assessed for a large site contaminated by trichloroethene. The impact of the remediation efforts (soil vapor extraction and in-situ chemical oxidation) was assessed through analysis of plume-scale contaminant mass discharge, which was measured using a high-resolution data set obtained from 23 years of operation of a large pump-and-treat system. The initial contaminant mass discharge peaked at approximately 7kg/d, and then declined to approximately 2kg/d. This latter value was sustained for several years prior to the initiation of source-zone remediation efforts. The contaminant mass discharge in 2010, measured several years after completion of the two source-zone remediation actions, was approximately 0.2kg/d, which is ten times lower than the value prior to source-zone remediation. The time-continuous contaminant mass discharge data can be used to evaluate the impact of the source-zone remediation efforts on reducing the time required to operate the pump-and-treat system, and to estimate the cost savings associated with the decreased operational period. While significant reductions have been achieved, it is evident that the remediation efforts have not completely eliminated contaminant mass discharge and associated risk. Remaining contaminant mass contributing to the current mass discharge is hypothesized to comprise poorly accessible mass in the source zones, as well as aqueous (and sorbed) mass present in the extensive lower-permeability units located within and adjacent to the contaminant plume. The fate of these sources is an issue of critical import to the remediation of chlorinated-solvent contaminated sites, and development of methods to address these sources will be required to achieve successful long-term management of such sites and to ultimately transition them to closure.     

An inverse Gaussian plume approach for estimating atmospheric pollutant emissions from multiple point sources

A method is developed for estimating the emission rates of contaminants into the atmosphere from multiple point sources using measurements of particulate material deposited at ground level. The approach is based on a Gaussian plume type solution for the advection–diffusion equation with ground-level deposition and given emission sources. This solution to the forward problem is incorporated into an inverse algorithm for estimating the emission rates by means of a linear least squares approach. The results are validated using measured deposition and meteorological data from a large lead–zinc smelting operation in Trail, British Columbia. The algorithm is demonstrated to be robust and capable of generating reasonably accurate estimates of total contaminant emissions over the relatively short distances of interest in this study.     

Relationship between mass-flux reduction and source-zone mass removal: analysis of field data

The magnitude of contaminant mass-flux reduction associated with a specific amount of contaminant mass removed is a key consideration for evaluating the effectiveness of a source-zone remediation effort. Thus, there is great interest in characterizing, estimating, and predicting relationships between mass-flux reduction and mass removal. Published data collected for several field studies were examined to evaluate relationships between mass-flux reduction and source-zone mass removal. The studies analyzed herein represent a variety of source-zone architectures, immiscible-liquid compositions, and implemented remediation technologies. There are two general approaches to characterizing the mass-flux-reduction/mass-removal relationship, end-point analysis and time-continuous analysis. End-point analysis, based on comparing masses and mass fluxes measured before and after a source-zone remediation effort, was conducted for 21 remediation projects. Mass removals were greater than 60% for all but three of the studies. Mass-flux reductions ranging from slightly less than to slightly greater than one-to-one were observed for the majority of the sites. However, these single-snapshot characterizations are limited in that the antecedent behavior is indeterminate. Time-continuous analysis, based on continuous monitoring of mass removal and mass flux, was performed for two sites, both for which data were obtained under water-flushing conditions. The reductions in mass flux were significantly different for the two sites (90% vs. approximately 8%) for similar mass removals ( approximately 40%). These results illustrate the dependence of the mass-flux-reduction/mass-removal relationship on source-zone architecture and associated mass-transfer processes. Minimal mass-flux reduction was observed for a system wherein mass removal was relatively efficient (ideal mass-transfer and displacement). Conversely, a significant degree of mass-flux reduction was observed for a site wherein mass removal was inefficient (non-ideal mass-transfer and displacement). The mass-flux-reduction/mass-removal relationship for the latter site exhibited a multi-step behavior, which cannot be predicted using some of the available simple estimation functions.     

Toxic chemicals in an abandoned phenolic waste site

The nature and extent of pollution was determined at the site of a former pinetar manufacturer. Compound distributions at various areas about the site revealed that, in addition to groundwater leaching of soluble phenolics, insoluble contaminants were spread by a dike-breach incident and subsequent construction activities. Differences in the patterns of chemicals in various wells suggested that more than one source of pollution occurred. The distribution of compounds about the site indicated that a general clean-up would not be cost-effective. Placement of an interceptor to collect groundwater seepage that contaminates surface water is being considered as an alternative.     

Simulation of the effect of remediation on EDB and 1,2-DCA plumes at sites contaminated by leaded gasoline

An analytical model is used to simulate the effects of partial source removal and plume remediation on ethylene dibromide (EDB) and 1,2-dichloroethane (1,2-DCA) plumes at contaminated underground storage tank (UST) sites. The risk posed by EDB, 1,2-DCA, and commingled gasoline hydrocarbons varies throughout the plume over time. Dissolution from the light nonaqueous phase liquid (LNAPL) determines the concentration of each contaminant near the source, but biological decay in the plume has a greater influence as distance downgradient from the source increases. For this reason, compounds that exceed regulatory standards near the source may not in downgradient plume zones. At UST sites, partial removal of a residual LNAPL source mass may serve as a stand alone remedial technique if dissolved concentrations in the source zone are within several orders of magnitude of the applicable government or remedial standards. This may be the case with 1,2-DCA; however, EDB is likely to be found at concentrations that are orders of magnitude higher than its low Maximum Contaminant Level (MCL) of 0.05 μg/L (micrograms per liter). For sites with significant EDB contamination, even when plume remediation is combined with source depletion, significant timeframes may be required to mitigate the impact of this compound. Benzene and MTBE are commonly the focus of remedial efforts at UST sites, but simulations presented here suggest that EDB, and to a lesser extent 1,2-DCA, could be the critical contaminants to consider in the remediation design process at many sites.     

Influence of mass transfer characteristics for DNAPL source depletion and contaminant flux in a highly characterized glaciofluvial aquifer

The transfer of contaminant mass between the nonaqueous- and aqueous-phases is a process of central importance for the remediation of sites contaminated by dense nonaqueous-phase liquids (DNAPLs). This paper describes a comparison of the results obtained with various alternative DNAPL-aqueous-phase mass transfer models contained in the literature for predicting DNAPL source-zone depletion times in groundwater systems. These dissolution models were largely developed through laboratory column experiments. To gain insight into the implications of various representations of the local-scale kinetic as well as equilibrium DNAPL dissolution processes, aquifer heterogeneity and the complex architecture of a DNAPL source-zone, the aqueous-phase contaminant concentrations and mass fluxes arriving at a down-gradient compliance boundary are analyzed in a conditional stochastic framework. The hydrogeologic setting is a heterogeneous fluvial aquifer in Southwest Germany, referred to as the aquifer analog dataset, that was intensively characterized in three dimensions for hydrogeological parameters that include permeability, effective porosity, grain size, mineralogy and sorption coefficients. By embedding the various dissolution models into the compositional, multiphase flow model, CompFlow, the relative times predicted for complete depletion of a released DNAPL source due to natural dissolution are explored. Issues related to achieving environmental benefits through, for example, partial DNAPL-zone source removal via enhanced remedial technologies are also discussed. In this context, performance metrics in the form of peak aqueous-phase contaminant concentrations and mass fluxes arriving at a down-gradient compliance boundary are compared to each other. This is done for each of the alternative mass transfer models. A significant reduction in the fractional flux at a downstream location from the DNAPL source can be achieved by partial source-zone mass reduction; however, peak concentration levels at the same location remain much higher than the United States Environment Protection Agency (US-EPA) drinking water limits. Although groundwater quality was found to improve more rapidly for the equilibrium dissolution model, it is also shown that dissolution models that promote rapid DNAPL disappearance produce greater prediction uncertainty in the aqueous-phase flux reduction.     

Redistribution of contaminants by a fluctuating water table in a micro-porous, double-porosity aquifer: field observations and model simulations

Large seasonal fluctuations of the water table are characteristic of aquifers with a low specific yield, including those fractured, double-porosity aquifers that have significant matrix porosity containing virtually immobile porewater, such as the Chalk of northern Europe. Where these aquifers are contaminated, a strong relationship between water table elevation and contaminant concentration in groundwater is commonly observed, of significance to the assessment, monitoring, and remediation of contaminated groundwater. To examine the processes governing contaminant redistribution by a fluctuating water table within the 'seasonally unsaturated zone', or SUZ, profiles of porewater solute concentrations have been established at a contaminated site in southern England. These profiles document the contaminant distribution in porewater of the Chalk matrix over the SUZ at a greater level of detail than recorded previously. A novel double-porosity solute transport code has been developed to simulate the evolution of the SUZ matrix porewater contaminant profiles, given a fluctuating water table, when the groundwater is initially contaminated and the SUZ is initially free of contamination. The model is simply characterised by: the matrix-fracture porosity ratio, the matrix block geometry, and a characteristic diffusion time. De-saturation and re-saturation of fractures is handled by a new approximation method. Contaminant accumulates in the upper levels of the SUZ, where it is less accessible to mobile groundwater, and acts as a persistent secondary source of contamination once the original source of contamination has been removed or has become depleted. The 'SUZ process' first attenuates the progress of contaminants in groundwater, and subsequently controls the slow release of contamination back to the mobile groundwater, thus prolonging the duration of groundwater contamination by many years. The SUZ process should operate in any fractured, micro-porous lithology e.g. fractured clays and mudstones, making this approach widely applicable.     

Sequential anaerobic-aerobic treatment of an aquifer contaminated by halogenated organics: field results

In situ, sequential, anaerobic to aerobic treatment of groundwater removed perchloroethene (PCE, 1.1 microM) and benzene (0.8 microM) from a contaminated aquifer. Neither aerobic nor anaerobic treatment alone successfully degraded both the chlorinated and non-chlorinated organic contaminants in the aquifer. After the sequential treatment, PCE, trichloroethene (TCE), vinyl chloride (VC), chloroethane (CA), and benzene were not detectable in groundwater. Desorption of residual aquifer contaminants was tested by halting the groundwater recirculation and analyzing the groundwater after 3 and 7 weeks. No desorption of the chlorinated contaminants or daughter products was observed in the treated portion of the aquifer. Sequential anaerobic to aerobic treatment was successful in remediating the groundwater at this test site and may have broad applications at other contaminated sites. Over the 4-year course of the project, the predominant microbial environment of the test site varied from aerobic to sulfate-reducing, to methanogenic, and back to aerobic conditions. Metabolically active microbial populations developed under all conditions, demonstrating the diversity and robustness of natural microbial flora in the aquifer.     

Attenuation reactions in a multiple contaminated aquifer in Bitterfeld (Germany)

Large-scale contaminated sites with multiple contaminants in the groundwater present a challenge to risk assessment and remediation. Attenuation reactions take place in the subsurface and act to contain contaminants, but must be thoroughly investigated on a site-specific basis. Field data from monitoring wells at a contaminated industrial site in Bitterfeld, Germany, are presented and analyzed for evidence of the prevalent biodegradation reactions. The groundwater in the Tertiary aquifer is contaminated with large quantities of chlorinated aliphatic compounds, in addition to chlorobenzenes and BTEX. In this strictly anaerobic environment, geochemical indications for several microbial processes were found, including methanogenesis, sulfate and iron reduction as well as reductive dechlorination of the chlorinated hydrocarbons. Direct evidence for the latter degradation reaction was observed along the flowpath due to the appearance of intermediates and an increase in the degree of dechlorination.     

Evaluating LNAPL contamination using GPR signal attenuation analysis and dielectric property measurements: practical implications for hydrological studies

Groundwater and sub-surface contamination by Light Non-Aqueous Phase Liquids (LNAPLs) is one of the industrial world's most pressing environmental issues and a thorough understanding of the hydrological, physical and bio-chemical properties of the sub-surface is key to determining the spatial and temporal development of any particular contamination event. Non-invasive geophysical techniques (such as electrical resistivity, electromagnetic conductivity, Ground-Penetrating Radar, etc.) have proved to be successful sub-surface investigation and characterisation tools with Ground-Penetrating Radar (GPR) being particularly popular. Recent studies have shown that the spatial/temporal variation in GPR signal attenuation can provide important information on the electrical properties of the sub-surface materials that, in turn, can be used to assess the physical and hydrological nature of the pore fluids and associated contaminants. Unfortunately, a high percentage of current LNAPL-related GPR studies focus on contaminant mapping only, with little emphasis being placed on characterising the hydrological properties (e.g., determining contaminant saturation index, etc.). By comparing laboratory-based, dielectric measurements of LNAPL contaminated materials with the GPR signal attenuation observed in both contaminated and 'clean' areas of an LNAPL contaminated site, new insights have been gained into the nature of contaminant distribution/saturation and the likely signal attenuation mechanisms. The results show that, despite some practical limitations of the analysis technique, meaningful hydrological interpretations can be obtained on the contaminant properties, saturation index and bio-degradation processes. A generalised attenuation/saturation model has been developed that describes the physical and attenuation enhancement characteristics of the contaminated areas and reveals that the most significant attenuation is related to smeared zone surrounding the seasonally changing water table interface. It is envisaged that the model will provide a basis for the interpretation of GPR data from analogous LNAPL contaminated sites and provide investigators with an appreciation of the merits and limitations of GPR-based, attenuation analysis techniques for hydrological applications.     

Modeling dissolution and volatilization of LNAPL sources migrating on the groundwater table

A vertically averaged two-dimensional model was developed to describe areal spreading and migration of light nonaqueous-phase liquids (LNAPLs) introduced into the subsurface by spills or leaks from underground storage tanks. The NAPL transport model was coupled with two-dimensional contaminant transport models to predict contamination of soil gas and groundwater resulting from a LNAPL migrating on the water table. Numerical solutions were obtained by using the finite-difference method. Simulations and sensitivity analyses were conducted with a LNAPL of pure benzene to study LNAPL migration and groundwater contamination. The model was applied to subsurface contamination by jet fuel. Results indicated that LNAPL migration were affected mostly by volatilization. The generation and movement of the dissolved plume was affected by the geology of the site and the free-product plume. Most of the spilled mass remained as a free LNAPL phase 20 years after the spill. The migration of LNAPL for such a long period resulted in the contamination of both groundwater and a large volume of soil.     

Changes in contaminant mass discharge from DNAPL source mass depletion: evaluation at two field sites

Changes in contaminant fluxes resulting from aggressive remediation of dense nonaqueous phase liquid (DNAPL) source zone were investigated at two sites, one at Hill Air Force Base (AFB), Utah, and the other at Ft. Lewis Military Reservation, Washington. Passive Flux Meters (PFM) and a variation of the Integral Pumping Test (IPT) were used to measure fluxes in ten wells installed along a transect down-gradient of the trichloroethylene (TCE) source zone, and perpendicular to the mean groundwater flow direction. At both sites, groundwater and contaminant fluxes were measured before and after the source-zone treatment. The measured contaminant fluxes (J; ML(-2)T(-1)) were integrated across the well transect to estimate contaminant mass discharge (M(D); MT(-1)) from the source zone. Estimated M(D) before source treatment, based on both PFM and IPT methods, were approximately 76 g/day for TCE at the Hill AFB site; and approximately 640 g/day for TCE, and approximately 206 g/day for cis-dichloroethylene (DCE) at the Ft. Lewis site. TCE flux measurements made 1 year after source treatment at the Hill AFB site decreased to approximately 5 g/day. On the other hand, increased fluxes of DCE, a degradation byproduct of TCE, in tests subsequent to remediation at the Hill AFB site suggest enhanced microbial degradation after surfactant flooding. At the Ft. Lewis site, TCE mass discharge rates subsequent to remediation decreased to approximately 3 g/day for TCE and approximately 3 g/day for DCE approximately 1.8 years after remediation. At both field sites, PFM and IPT approaches provided comparable results for contaminant mass discharge rates, and show significant reductions (>90%) in TCE mass discharge as a result of DNAPL mass depletion from the source zone.     

A contaminated site investigation: comparison of information gained from geophysical measurements and hydrogeological modeling

The investigation of contaminated sites is usually a long and expensive process. It is therefore desirable to use a combination of methodologies in an integrative approach that can reduce redundant information gathering. The objective of this study was to examine the usefulness of 2 non-intrusive exploration techniques in a contaminated site investigation. Borehole positioning based on geophysical measurements was compared to positioning based on the Bayesian expert system for flow-field modeling. The goal set at the field site was the assessment of the type and load of contaminants transported from the landfill site to the adjacent aquifer and the extent of leachate plumes within the groundwater. The two methods made different demands on information gathering and were found to be complementary. The geophysical approach focused attention on the waste compartments at the site and on mixing plumes in the adjacent aquifer but could not, without prior information, provided information on the flow field. The Bayesian approach to flow-field modeling determined areas of greatest model uncertainty at the model boundaries. The model highlighted areas of greatest uncertainty that might otherwise have been overlooked and provided information on the most likely mean direction of the leachate plumes. It was concluded that both methods contribute to a site investigation and should be used before additional drilling is carried out.     

Quantification of groundwater contamination in an urban area using integral pumping tests

In this paper, the integral groundwater investigation method is used for the quantification of PCE and TCE mass flow rates at an industrialized urban area in Linz, Austria. In this approach, pumping wells positioned along control planes perpendicular to the groundwater flow direction are operated for a time period on the order of days and sampled for contaminants. The concentration time series of the contaminants measured during operation of the pumping wells are then used to determine contaminant mass flow rates, mean concentrations and the plume shapes and positions at the control planes. The three control planes used in Linz were positioned downstream of a number of potential source zones, which are distributed over the field site. By use of the integral investigation method, it was possible to identify active contaminant sources, quantify the individual source strength in terms of mass flow rates at the control planes and estimate the contaminant plume position relative to the control planes. The source zones emitting the highest PCE and TCE mass flow rates could be determined, representing the areas where additional investigation and remediation activities will be needed. Additionally, large parts of the area investigated could be excluded from further investigation and remediation activities.