首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 24 毫秒
1.
The process of attempting to model ground-water systems requires a good understanding of the spatial variation of aquifer hydraulic properties. The capabilities of the more recent innovative flowmeters such as the electromagnetic and heat pulse flowmeters provide the sensitivity to measure ambient flows and pump-induced flows. These flowmeters provide the measurements of pump-induced vertical flows which are analyzed to obtain vertical variations in horizontal hydraulic conductivity, K(z). With discrete areal K-values, K(x, y), and vertical profiles of K, provided by multiwell testing, the essential elements are present to produce a three-dimensional hydraulic conductivity field. The advent of these new flow measuring devices has contributed much to the motivation behind this paper. This paper presents the results of applying deterministic and stochastic methodology to the three-dimensional interpolation of hydraulic properties, specifically, hydraulic conductivity, K. Three of the approaches applied in this paper are deterministic in nature, inverse-distance weighting, inverse-distance-squared weighting, and ordinary kriging, while the fourth is a stochastic approach based on self-affine fractals. All of the methods are applied to measured data collected from 14 wells at a site in the United States near Mobile, Alabama. The three-dimensional K-distributions generated by each of the methods are used as inputs to an advective based transport model with the resulting model output compared to a two-well tracer study run previously at the same site.  相似文献   

2.
A new tracer experiment (referred to as MADE‐5) was conducted at the well‐known Macrodispersion Experiment (MADE) site to investigate the influence of small‐scale mass‐transfer and dispersion processes on well‐to‐well transport. The test was performed under dipole forced‐gradient flow conditions and concentrations were monitored in an extraction well and in two multilevel sampler (MLS) wells located at 6, 1.5, and 3.75 m from the source, respectively. The shape of the breakthrough curve (BTC) measured at the extraction well is strongly asymmetric showing a rapidly arriving peak and an extensive late‐time tail. The BTCs measured at seven different depths in the two MLSs are radically different from one another in terms of shape, arrival times, and magnitude of the concentration peaks. All of these characteristics indicate the presence of a complex network of preferential flow pathways controlling solute transport at the test site. Field‐experimental data were also used to evaluate two transport models: a stochastic advection‐dispersion model (ADM) based on conditional multivariate Gaussian realizations of the hydraulic conductivity field and a dual‐domain single‐rate (DDSR) mass‐transfer model based on a deterministic reconstruction of the aquifer heterogeneity. Unlike the stochastic ADM realizations, the DDSR accurately predicted the magnitude of the concentration peak and its arrival time (within a 1.5% error). For the multilevel BTCs between the injection and extraction wells, neither model reproduced the observed values, indicating that a high‐resolution characterization of the aquifer heterogeneity at the subdecimeter scale would be needed to fully capture 3D transport details.  相似文献   

3.
A tracer test was used to evaluate whether cross contamination exists along a monitoring well completed through a shallow ground water system in fractured clay and screened in a sand and gravel aquifer. The fractured clay is separated from the sand and gravel deposit by a layer of highly plastic unfractured clay. A natural vertical downward hydraulic gradient of approximately 0.5 exists between the shallow system and the sand and gravel aquifer. Ground water contamination was detected in an adjacent monitoring well screened in the fractured clay and in the monitoring well screened in the sand and gravel deposit. No ground water contamination was apparent in an intermediate well screened in the unfractured clay layer. A tracer of sodium bromide was injected into a shallow boring near the monitoring wells. The tracer was detected in the monitoring well in the sand and gravel aquifer after three to seven days. The bromide concentration continued to increase in this well with time while the concentration in the shallow boring declined. This trend of tracer concentration indicates the tracer has in fact migrated downward and possibly traveled along the well column.  相似文献   

4.
Illman WA  Berg SJ  Yeh TC 《Ground water》2012,50(3):421-431
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.  相似文献   

5.
We investigate the importance of selecting two different methodologies for the determination of hydraulic conductivity from available grain-size distributions on the stochastic modeling of the depth-averaged breakthrough curve observed during a forced-gradient tracer test experiment. The latter was performed in the Lauswiesen alluvial aquifer, located near the city of Tübingen, Germany, by injecting NaBr into a well at a distance of about 50 m from a pumping well. We also examine the joint effect of the choice of the transport model adopted to describe solute transport at the site and the way the spatial distribution of porosity is assessed. In the absence of direct measurements of porosity, we consider: (a) the model used by Riva et al. (J Contam Hydrol 88:92–118, 2006; J Contam Hydrol 101:1–13, 2008), which relates the natural logarithms of effective porosity and conductivity through an empirical, experimentally-based, linear relationship derived for a nearby experimental site; and (b) a model based on a commonly used relationship linking the total porosity to the coefficient of uniformity of grain size distributions. Transport is described in terms of a purely advective process and/or by including mass exchange processes between mobile and immobile regions. Modeling of flow and transport is performed within a Monte Carlo framework, upon conceptualizing the aquifer as a random composite medium. Our results indicate that the model adopted to describe the correlation between conductivity and porosity and the way grain-sieve information are incorporated to depict the heterogeneous distribution of hydraulic conductivity can have relevant effects in the interpretation of the data at the site. All the conceptual models employed to describe the structural heterogeneity of the system and transport features can reasonably reproduce the global characteristics of the experimental depth-averaged breakthrough curve. Specific details, such as the peak concentration and the time of first arrival, can be better reproduced by a double porosity transport model when a correlation between conductivity and porosity based on grain size information at the site is considered. The best prediction of the late-time behavior of the measured breakthrough curves, in terms of the observed heavy tailing, is offered by directly linking porosity distribution to the spatial variability of particle size information.  相似文献   

6.
Accurate quantification of in situ heterogeneity and flow processes through fractured geologic media remains elusive for hydrogeologists due to the complexity in fracture characterization and its multiscale behavior. In this research, we demonstrated the efficacy of tracer-electrical resistivity tomography (ERT) experiments combined with numerical simulations to characterize heterogeneity and delineate preferential flow paths in a fractured granite aquifer. A series of natural gradient saline tracer experiments were conducted from a depth window of 18 to 22 m in an injection well (IW) located inside the Indian Institute of Technology Hyderabad campus. Tracer migration was monitored in a time-lapse mode using two cross-sectional surface ERT profiles placed in the direction of flow gradient. ERT data quality was improved by considering stacking, reciprocal measurements, resolution indicators, and geophysical logs. Dynamic changes in subsurface electrical properties inferred via resistivity anomalies were used to highlight preferential flow paths of the study area. Temporal changes in electrical resistivity and tracer concentration were monitored along the vertical in an observation well located at 48 m to the east of the IW. ERT-derived tracer breakthrough curves were in agreement with geochemical sample measurements. Fracture geometry and hydraulic properties derived from ERT and pumping tests were further used to evaluate two mathematical conceptualizations that are relevant to fractured aquifers. Results of numerical analysis conclude that dual continuum model that combines matrix and fracture systems through a flow exchange term has outperformed equivalent continuum model in reproducing tracer concentrations at the monitoring wells (evident by a decrease in RMSE from 199 to 65 mg/L). A sensitivity analysis on model simulations conclude that spatial variability in hydraulic conductivity, local-scale dispersion, and flow exchange at fracture-matrix interface have a profound effect on model simulations.  相似文献   

7.
Traditional numerical methods for the delineation of wellhead protection areas span deterministic and probabilistic approaches. They provide time-related capture zones. However, none of the existing approaches identifies the groundwater contribution areas related to each source or sink. In this work, the worthiness of the so-called double delineation approach was extended. This task was achieved by simple postprocessing of its dual outputs leading to a highly efficient screening tool. In the particular context of geothermal resources management through the well doublets of the Dogger aquifer in the Paris Basin (France), the approach was extended to forecast the compositional heat breakthrough at production wells. Hence, cold-water breakthrough and temperature decline in production wells are timely assessed in low-enthalpy geothermal reservoirs. The method quantifies how groundwater volumes are moving through space and time between any couple of source and sink. It provides unprecedented tools advancing the enhanced understanding of water resources systems functioning. It is highly recommended to implement the presented concepts in the current and future generations of community groundwater models.  相似文献   

8.
Water from the San Joaquin Delta, having chloride concentrations up to 3590 mg/L, has intruded fresh water aquifers underlying Stockton, California. Changes in chloride concentrations at depth within these aquifers were evaluated using sequential electromagnetic (EM) induction logs collected during 2004 through 2007 at seven multiple‐well sites as deep as 268 m. Sequential EM logging is useful for identifying changes in groundwater quality through polyvinyl chloride‐cased wells in intervals not screened by wells. These unscreened intervals represent more than 90% of the aquifer at the sites studied. Sequential EM logging suggested degrading groundwater quality in numerous thin intervals, typically between 1 and 7 m in thickness, especially in the northern part of the study area. Some of these intervals were unscreened by wells, and would not have been identified by traditional groundwater sample collection. Sequential logging also identified intervals with improving water quality—possibly due to groundwater management practices that have limited pumping and promoted artificial recharge. EM resistivity was correlated with chloride concentrations in sampled wells and in water from core material. Natural gamma log data were used to account for the effect of aquifer lithology on EM resistivity. Results of this study show that a sequential EM logging is useful for identifying and monitoring the movement of high‐chloride water, having lower salinities and chloride concentrations than sea water, in aquifer intervals not screened by wells, and that increases in chloride in water from wells in the area are consistent with high‐chloride water originating from the San Joaquin Delta rather than from the underlying saline aquifer.  相似文献   

9.
Heterogeneity in the physical properties of an aquifer can significantly affect the viability of aquifer storage and recovery (ASR) by reducing the recoverable proportion of low-salinity water where the ambient ground water is brackish or saline. This study investigated the relationship between knowledge of heterogeneity and predictions of solute transport and recovery efficiency by combining permeability and ASR-based tracer testing with modeling. Multiscale permeability testing of a sandy limestone aquifer at an ASR trial site showed that small-scale core data give lower-bound estimates of aquifer hydraulic conductivity (K), intermediate-scale downhole flowmeter data offer valuable information on variations in K with depth, and large-scale pumping test data provide an integrated measure of the effective K that is useful to constrain ground water models. Chloride breakthrough and thermal profiling data measured during two cycles of ASR showed that the movement of injected water is predominantly within two stratigraphic layers identified from the flowmeter data. The behavior of the injectant was reasonably well simulated with a four-layer numerical model that required minimal calibration. Verification in the second cycle achieved acceptable results given the model's simplicity. Without accounting for the aquifer's layered structure, high precision could be achieved on either piezometer breakthrough or recovered water quality, but not both. This study demonstrates the merit of an integrated approach to characterizing aquifers targeted for ASR.  相似文献   

10.
The comparison between two series of optimal remediation designs using deterministic and stochastic approaches showed a number of converging features. Limited sampling measurements in a supposed contaminated aquifer formed the hydraulic conductivity field and the initial concentration distribution used in the optimization process. The deterministic and stochastic approaches employed a single simulation–optimization method and a multiple realization approach, respectively. For both approaches, the optimization model made use of a genetic algorithm. In the deterministic approach, the total cost, extraction rate, and the number of wells used increase when the design must satisfy the intensified concentration constraint. Growing the stack size in the stochastic approach also brings about same effects. In particular, the change in the selection frequency of the used extraction wells, with increasing stack size, for the stochastic approach can indicate the locations of required additional wells in the deterministic approach due to the intensified constraints. These converging features between the two approaches reveal that a deterministic optimization approach with controlled constraints is achievable enough to design reliable remediation strategies, and the results of a stochastic optimization approach are readily available to real contaminated sites.  相似文献   

11.
Two-well tracer tests are often conducted to investigate subsurface solute transport in the field. Analyzing breakthrough curves in extraction and monitoring wells using numerical methods is nontrivial due to highly nonuniform flow conditions. We extended approximate analytical solutions for the advection-dispersion equation for an injection-extraction well doublet in a homogeneous confined aquifer under steady-state flow conditions for equal injection and extraction rates with no transverse dispersion and negligible ambient flow, and implemented the solutions in Microsoft Excel using Visual Basic for Application (VBA). Functions were implemented to calculate concentrations in extraction and monitoring wells at any location due to a step or pulse injection. Type curves for a step injection were compared with those calculated by numerically integrating the solution for a pulse injection. The results from the two approaches are similar when the dispersivity is small. As the dispersivity increases, the latter was found to be more accurate but requires more computing time. The code was verified by comparing the results with published-type curves and applied to analyze data from the literature. The method can be used as a first approximation for two-well tracer test design and data analysis, and to check accuracy of numerical solutions. The code and example files are publicly available.  相似文献   

12.
To better understand the groundwater resources of southern Nye County, Nevada, a multipart distributed thermal perturbation sensing (DTPS) test was performed on a complex of three wells. These wells penetrate an alluvial aquifer that drains the Nevada National Security Site, and characterizing the hydraulic properties and flow paths of the regional groundwater flow system has proven very difficult. The well complex comprised one pumping well and two observation wells, both located 18 m from the pumping well. Using fiber‐optic cables and line heaters, DTPS tests were performed under both stressed and unstressed conditions. Each test injects heat into the water column over a period of one to two days, and observes the rising temperature during heat injection and falling temperatures after heating ceases. Aquifer thermal properties are inferred from temperature patterns in the cased section of the wells, and fluxes through the 30‐m screened section are estimated based on a model that incorporates conductive and advective heat fluxes. Vertical variations in flux are examined on a scale of tens of cm. The actively flowing zones of the aquifer change between the stressed and unstressed test, and anisotropy in the aquifer permeability is apparent from the changing fluxes between tests. The fluxes inferred from the DTPS tests are compared to solute tracer tests previously performed on the same site. The DTPS‐based fluxes are consistent with the fastest solute transport observed in the tracer test, but appear to overestimate the mean flux through the system.  相似文献   

13.
The geologic architecture in aquifer systems affects the behavior of fluid flow and the dispersion of mass. The spatial distribution and connectivity of higher-permeability facies play an important role. Models that represent this geologic structure have reduced entropy in the spatial distribution of permeability relative to models without structure. The literature shows that the stochastic model with the greatest variance in the distribution of predictions (i.e., the most conservative model) will not simply be the model representing maximum disorder in the permeability field. This principle is further explored using the Shannon entropy as a single metric to quantify and compare model parametric spatial disorder to the temporal distribution of mass residence times in model predictions. The principle is most pronounced when geologic structure manifests as preferential-flow pathways through the system via connected high-permeability sediments. As per percolation theory, at certain volume fractions the full connectivity of the high-permeability sediments will not be represented unless the model is three-dimensional. At these volume fractions, two-dimensional models can profoundly underrepresent the entropy in the real, three-dimensional, aquifer system. Thus to be conservative, stochastic models must be three-dimensional and include geologic structure.  相似文献   

14.
Wells designated as groundwater under the direct influence (GUDI) of surface water have caused an ongoing boil-water advisory afflicting the island of Tutuila, American Samoa for almost a decade. Regulatory testing at these wells found turbidity and indicator bacteria spikes correlated with heavy rainfall events. However, the mechanism of this contamination has, until now, remained unknown. Surface water may reach wells through improperly sealed well casings, or through the aquifer matrix itself. In this study, three independent surface water tracers, turbidity, indicator bacteria, and water isotopes were used to assess recharge timing and determine contamination mechanisms. Results from each method were reasonably consistent, revealing average GUDI well breakthrough times of 37 ± 21 h for turbidity, 18 to 63 h for bacteria, and 1 to 5 days for water isotopes. These times match well with estimated subsurface flow rates through highly permeable aquifer materials. In contrast, where one well casing was found to be compromised, turbidity breakthrough was observed at 3 to 4 h. These results support local management decisions and show repairing or replacing wells will likely result in continued GUDI contamination. Additionally, differences in observed rainfall response for each tracer provide insight into the recharge dynamics and subsurface flow characteristics of this and other highly conductive young-basaltic aquifers.  相似文献   

15.
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.  相似文献   

16.
Discrete-fracture and dual-porosity models are infrequently used to simulate solute transport through fractured unconsolidated deposits, despite their more common application in fractured rock where distinct flow regimes are hypothesized. In this study, we apply four fracture transport models--the mobile-immobile model (MIM), parallel-plate discrete-fracture model (PDFM), and stochastic and deterministic discrete-fracture models (DFMs)--to demonstrate their utility for simulating solute transport through fractured till. Model results were compared to breakthrough curves (BTCs) for the conservative tracers potassium bromide (KBr), pentafluorobenzoic acid (PFBA), and 1,4-piperazinediethanesulfonic acid (PIPES) in a large-diameter column of fractured till. Input parameters were determined from independent field and laboratory methods. Predictions of Br BTCs were not significantly different among models; however, the stochastic and deterministic DFMs were more accurate than the MIM or PDFM when predicting PFBA and PIPES BTCs. DFMs may be more applicable than the MIM for tracers with small effective diffusion coefficients (De) or for short timescales due to differences in how these models simulate diffusion or incorporate heterogeneities by their fracture networks. At large scales of investigation, the more computationally efficient MIM and PDFM may be more practical to implement than the three-dimensional DFMs, or a combination of model approaches could be employed. Regardless of the modeling approach used, fractures should be incorporated routinely into solute transport models in glaciated terrain.  相似文献   

17.
Characterization of a multilayer aquifer using open well dilution tests   总被引:1,自引:0,他引:1  
West LJ  Odling NE 《Ground water》2007,45(1):74-84
An approach to characterization of multilayer aquifer systems using open well borehole dilution is described. The approach involves measuring observation well flow velocities while a nearby extraction well is pumped by introducing a saline tracer into observation wells and collecting dilution vs. depth profiles. Inspection of tracer profile evolution allows discrete permeable layers within the aquifer to be identified. Dilution profiles for well sections between permeable layers are then converted into vertical borehole flow velocities and their evolution, using an analytic solution to the advection-dispersion equation applied to borehole flow. The dilution approach is potentially able to measure much smaller flow velocities that would be detectable using flowmeters. Vertical flow velocity data from the observation wells are then matched to those generated using a hydraulic model of the aquifer system, "shorted" by the observation wells, to yield the hydraulic properties of the constituent layers. Observation well flow monitoring of pumping tests represents a cost-effective alternative or preliminary approach to pump testing each layer of a multilayer aquifer system separately using straddle packers or screened wells and requires no prior knowledge of permeable layer depths and thicknesses. The modification described here, of using tracer dilution rather than flowmeter logging to obtain well flow velocities, allows the approach to be extended to greater well separations, thus characterizing a larger volume of the aquifer. An example of the application of this approach to a multilayer Chalk Aquifer in Yorkshire, Northeast England, is presented.  相似文献   

18.
19.
This paper reviews different borehole flowmeter analysis methods and evaluates their applicability to a test site composed of fluvial deposits. Results from tracer and aquifer tests indicate that the aquifer is highly heterogeneous and that low-K skin effects exist at the wells. Borehole flowmeter tests were performed at 37 wells. An appropriate method for calculating borehole flowmeter K values was developed based on results from multiwell pumping tests, single-well pumping tests, and slug tests. The flowmeter data produced 881 K values. The trends and the magnitude of the K values are consistent with results from geologic investigations, recirculating tracer tests, and large-scale multiwell pumping tests. The field tests illustrate that high-K deposits can significantly affect ground-water flows in some heterogeneous fluvial aquifers.  相似文献   

20.
A tracer test was conducted to characterize the flow of groundwater across a permeable reactive barrier constructed with plant mulch (a biowall) at the OU‐1 site on Altus Air Force Base, Oklahoma. This biowall is intended to intercept and treat groundwater contaminated by trichloroethylene (TCE) in a shallow aquifer. The biowall is 139‐m long, 7.3‐m deep, and 0.5‐m wide. Bromide was injected from an upgradient well into the groundwater as a conservative tracer, and was subsequently observed breaking through in monitoring wells within and downgradient of the biowall. The bromide breakthrough data demonstrate that groundwater entering the biowall migrated across it, following the slope of the local groundwater surface. The average seepage velocity of groundwater was approximately 0.06 m/d. On the basis of the Darcy velocity of groundwater and geometry of the biowall, the average residence time of groundwater in the biowall was estimated at 10 d. Assuming all TCE removal occurred in the biowall, the reduction in TCE concentrations in groundwater across the biowall corresponds to a first‐order attenuation rate constant in the range of 0.38 to 0.15 per d. As an independent estimate of the degradation rate constant, STANMOD software was used to fit curves through data on the breakthrough of bromide and TCE in selected wells downgradient of the injection wells. Best fits to the data required a first‐order degradation rate constant for TCE removal in the range of 0.13 to 0.17 per d. The approach used in this study provides an objective evaluation of the remedial performance of the biowall that can provide a basis for design of other biowalls that are intended to remediate TCE‐contaminated groundwater.  相似文献   

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

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