首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 375 毫秒
1.
Combining groundwater flow models with solute transport models represents a common challenge in groundwater resources assessments and contaminant transport modeling. Groundwater flow models are usually constructed at somewhat larger scales (involving a coarser discretization) to include natural boundary conditions. They are commonly calibrated using observed groundwater levels and flows (if available). The groundwater solute transport models may be constructed at a smaller scale with finer discretization than the flow models in order to accurately delineate the solute source and the modeled target, to capture any heterogeneity that may affect contaminant migration, and to minimize numerical dispersion while still maintaining a reasonable computing time. The solution that is explored here is based on defining a finer grid subdomain within a larger coarser domain. The local-grid refinement (LGR) implemented in the Modular 3D finite-difference ground-water flow model (MODFLOW) code has such a provision to simulate groundwater flow in two nested grids: a higher-resolution sub-grid within a coarse grid. Under the premise that the interface between both models was well defined, a comprehensive sensitivity and uncertainty analysis was performed whereby the effect of a parameter perturbation in a coarser-grid model on transport predictions using a higher-resolution grid was quantified. This approach was tested for a groundwater flow and solute transport analysis in support of a safety evaluation of the future Belgian near-surface radioactive waste disposal facility. Our reference coarse-grid groundwater flow model was coupled with a smaller fine sub-grid model in two different ways. While the reference flow model was calibrated using observed groundwater levels at a scale commensurate with that of the coarse-grid model, the fine sub-grid model was used to run a solute transport simulation quantifying concentrations in a hypothetical well nearby the disposal facility. When LGR coupling was compared to a one-way coupling, LGR was found to provide a smoother flow solution resulting in a more CPU-efficient transport solution. Parameter sensitivities performed with the groundwater flow model resulted in sensitivities at the head observation locations. These sensitivities identified the recharge as the most sensitive parameter, with the hydraulic conductivity of the upper aquifer as the second most sensitive parameter in regard to calculated groundwater heads. Based on one-percent sensitivity maps, the spatial distribution of the observations with the highest sensitivities is slightly different for the upper aquifer hydraulic conductivity than for recharge. Sensitivity analyses were further performed to assess the prediction scaled sensitivities for hypothetical contaminant concentrations using the combined groundwater flow and solute transport models. Including all pertinent parameters into the sensitivity analysis identified the hydraulic conductivity of the upper aquifer as the most sensitive parameter with regard to the prediction of contaminant concentrations.  相似文献   

2.
Flow and solute transport monitoring in the karst aquifer in SW Slovenia   总被引：1，自引：0，他引：1  
The role of the unsaturated zone in the karst aquifer hydraulic behaviour was brought into focus in these studies of the catchment of the Hubelj spring (SW Slovenia). The variations of natural tracers in precipitation and in groundwater during a summer storm event made it possible to trace local flow and solute transport in the observed aquifer. The results produced data on the aquifer recharge, storage and discharge processes, as well as on mechanisms that affected them, which reflects a karst groundwater dynamics also at a regional scale. They point out the significance of effects of the fast preferential flow—epiflow that is the main factor controlling solute/contaminant transport towards the aquifer saturated zone. Numerous arguments indicate that the karst aquifer flow and solute transport mechanisms depend on the hydraulic behaviour of the epikarst zone.  相似文献   

3.
潜水水流的动态随机模拟是一个复杂而难解决的问题.通过建立二维潜水非稳定流模拟的摄动随机有限元模型, 把控制方程的主要参数渗透系数和给水度随机变量、及源汇项和边界条件看作随机变量.在充分考虑4种随机因素的条件下, 推导出求解潜水二维非稳定流均值和方差的9个方程; 重点介绍了不同方程数值离散的特殊处理方法.通过设定理想例子对模拟结果进行了分析, 表明随机变量中边界条件值方差、渗透系数方差变化对水头方差变化的影响很小, 给水度方差的变化对水头方差的变化影响很大.本模型考虑因素全面, 对一般的潜水非稳定流随机模拟都可应用.本研究给出了边界、渗透系数、给水度的随机因素对潜水动态模拟的影响, 丰富和补充了地下水运动的随机理论.   相似文献   

4.
Layers of strong geologic contrast within the unsaturated zone can control recharge and contaminant transport to underlying aquifers. Slow diffuse flow in certain geologic layers, and rapid preferential flow in others, complicates the prediction of vertical and lateral fluxes. A simple model is presented, designed to use limited geological site information to predict these critical subsurface processes in response to a sustained infiltration source. The model is developed and tested using site-specific information from the Idaho National Laboratory in the Eastern Snake River Plain (ESRP), USA, where there are natural and anthropogenic sources of high-volume infiltration from floods, spills, leaks, wastewater disposal, retention ponds, and hydrologic field experiments. The thick unsaturated zone overlying the ESRP aquifer is a good example of a sharply stratified unsaturated zone. Sedimentary interbeds are interspersed between massive and fractured basalt units. The combination of surficial sediments, basalts, and interbeds determines the water fluxes through the variably saturated subsurface. Interbeds are generally less conductive, sometimes causing perched water to collect above them. The model successfully predicts the volume and extent of perching and approximates vertical travel times during events that generate high fluxes from the land surface. These developments are applicable to sites having a thick, geologically complex unsaturated zone of substantial thickness in which preferential and diffuse flow, and perching of percolated water, are important to contaminant transport or aquifer recharge.  相似文献   

5.
Sedimentological processes often result in complex three-dimensional subsurface heterogeneity of hydrogeological parameter values. Variogram-based stochastic approaches are often not able to describe heterogeneity in such complex geological environments. This work shows how multiple-point geostatistics can be applied in a realistic hydrogeological application to determine the impact of complex geological heterogeneity on groundwater flow and transport. The approach is applied to a real aquifer in Belgium that exhibits a complex sedimentary heterogeneity and anisotropy. A training image is constructed based on geological and hydrogeological field data. Multiple-point statistics are borrowed from this training image to simulate hydrofacies occurrence, while intrafacies permeability variability is simulated using conventional variogram-based geostatistical methods. The simulated hydraulic conductivity realizations are used as input to a groundwater flow and transport model to investigate the effect of small-scale sedimentary heterogeneity on contaminant plume migration. Results show that small-scale sedimentary heterogeneity has a significant effect on contaminant transport in the studied aquifer. The uncertainty on the spatial facies distribution and intrafacies hydraulic conductivity distribution results in a significant uncertainty on the calculated concentration distribution. Comparison with standard variogram-based techniques shows that multiple-point geostatistics allow better reproduction of irregularly shaped low-permeability clay drapes that influence solute transport.  相似文献   

6.
Removal of LNAPL (oil) from an aquifer is described using a multiphase flow model. At the well boundary seepage face conditions are imposed. These conditions are implemented in a numerical model and withdrawal in a twodimensional domain is simulated for two different geometries of the oil lens and for varied values of the physical parameters. Assuming vertical equilibrium, the oil flow equation is reduced by vertical integration. The well boundary condition is approximated by imposing zero oil lens thickness. Similarity solutions of the reduced equations for the two geometries show good agreement with the numerical results in most cases.  相似文献   

7.
The hydrological response of the Choshuishi alluvial fan to the 1999 Chi-Chi earthquake shows that the earthquake did impact the aquifer. The possible earthquake-induced changes in hydrogeological properties were investigated in this study. First, contour maps of the hydrologic anomaly, seismic factors, and vertical ground-surface displacement were compared qualitatively. Bulls eye patterns were found on the contour maps of hydraulic conductivity, coseismic groundwater-level change and vertical ground-surface displacement but did not occur with other seismic factors. The more permeable zones of the aquifer were found to coincide with the locations of greater vertical ground-surface displacement and coseismic groundwater-level change in the 1999 Chi-Chi earthquake. This indicates that the change of the hydrogeologic properties of Choshuishi alluvial fan due to the 1999 Chi-Chi earthquake may have mainly occurred in the highly permeable zones. Fractal, cross semivariogram and cross correlogram analyses were performed to quantitatively measure the persistency, variability and similarity, respectively, of spatial hydrologic response, seismic factors and hydraulic conductivity. The groundwater-level change, earthquake intensity, and vertical ground-surface displacement were found to show antipersistent tendencies while other factors showed the opposite. Higher correlations were found between hydraulic conductivity and groundwater-level change in aquifers 2–1 and 2–2, and between hydraulic conductivity and vertical ground-surface displacement in aquifer 3. Changes in porosities and hydraulic conductivity were evaluated in the main aquifers of the Choshuishi alluvial fan based on the data of hydrologic anomaly and the vertical ground-surface displacement. While both approaches show that the 1999 Chi-Chi earthquake has impacted the Choshuishi alluvial fan by reducing its porosity and hydraulic conductivity, these changes were not significant relative to natural variation in hydraulic conductivity.This revised version was published in May 2005 with correction to the rubric.  相似文献   

8.
The deep groundwater in the quaternary gravel sequence of the southern Upper Rhine Graben locally contains high chloride concentrations near the river Rhine between Fessenheim (France) in the South and Breisach (Germany) in the North. This historical pollution is mainly due to past infiltration from the former brine storage basins of the French potash mines on the Fessenheim Island and—to a lesser extent—from the leaching of the salt dumps of the German potash mines in Buggingen and Heitersheim. The spreading of the salt plume was investigated by means of a groundwater model. The aim of the model was to understand the brine movement, the present distribution of chloride as defined by recent hydrochemical investigations, and to select locations for new reconnaissance boreholes. The geological structure was reproduced by a three layer model, which was calibrated for steady state flow conditions. The hydraulic conductivity of the first layer was determined by comparing measured and calculated heads in the model area. The vertical resolution was refined to simulate the density-dependent salt transport processes. The transport of the salt plumes was simulated over a 40-year period, starting at the beginning of brine storage in the 1950s. The relevant transport parameters have been estimated in a sensitivity analysis, where the simulated breakthrough curves of chloride concentration have been compared with the measured data. The results of the groundwater model indicate that brines containing approximately 1 million tons of chloride are still present at the bottom of the aquifer. These highly concentrated salt brines mix with fresh water from the upper part of the aquifer. This dispersive process leads to the formation of a plume of chloride-rich water extending downstream, where pumping wells for several local water supplies are located.  相似文献   

9.
山东淄博市大武水源地裂隙岩溶水中污染物运移的数值研究   总被引：11，自引：0，他引：11  
朱学愚  刘建立 《地学前缘》2001,8(1):171-178
在分析研究淄博市大武水源地裂隙岩溶含水层的水力性质和污染物运移特征的基础上 ,对裂隙岩溶水的水头和污染物运移进行数值研究。目前国内外对裂隙岩溶水进行数值计算时 ,通常用等价多孔介质模型 ,但裂隙岩溶介质和多孔介质有很大不同。裂隙岩溶介质的储水和导水空间为裂隙网络 ,导水系数大 ,地下水的实际平均流速比孔隙水大得多 ,但给水度和贮水系数小。当用等价多孔介质模型进行模拟时应考虑这些特点。对于污染物运移的模拟 ,要同时求解水头方程和对流弥散方程 ,可采用MODFLOW和MT3D软件进行模拟。研究区裂隙岩溶水水头的数值计算表明 ,等效多孔介质模型水头的拟合误差能满足国标GB/T144 97- 93的要求。各时段地下水水量均衡计算的精度也满足要求。对流弥散方程的数值计算 ,由于Peclet数高达 95 .6 7,对流占绝对优势 ,可能存在数值弥散和数值振荡 ,因而采用多种方法进行了比较。对于同一问题 ,同时采用上游有限差分法 (UFDM) ,混合的欧拉拉格朗日方法 (特征线法MOC、改进特征线法MMOC和混合特征线法HMOC) ,总变异消减法(TVD)进行计算 ,并比较其结果。结果表明 ,混合特征线法 (HMOC)和总变异消减法 (TVD)比较适合于对流占优势的运移问题计算。由于渗透系数K和有效孔隙度θ对溶质运移结果的影响很大 ,?  相似文献   

10.
Hydrogeological research was performed to establish the origin of organic-rich water in a Quaternary aquifer (Wielkopolska Buried Valley aquifer, Poland). The spatial distribution of such water was determined by means of multilevel sampling with a packer. The nature, spatial distribution and chemistry of the organic-rich water suggest hydraulic connection between the Quaternary aquifer and the underlying Neogene aquifer. This connection may be enabled by an old abandoned, improperly plugged well, and must therefore be considered as an artificial hydraulic connection between hydrogeologically separated aquifers. For the verification of this hypothesis, a groundwater contaminant transport model was constructed. The results of contaminant transport modeling allow the contaminant transport parameters (transverse and longitudinal dispersivity) to be identified in the field. The numerical groundwater model was used for the prognosis of organic-rich water remediation, and for specification of water-resource management principles in the region where organic-rich water occurs. The results of this research enable the verification of assumptions regarding complex groundwater flow conditions in the boundary zone of a buried valley.  相似文献   

11.
Groundwater aquifer vulnerability has been assessed by incorporating the major geological and hydrogeological factors that affect and control the groundwater contamination using GIS-based DRASTIC model along with solute transport modeling. This work demonstrates the potential of GIS to derive a vulnerability map by overlying various spatially referenced digital data layers (i.e., depth to water, net recharge, aquifer media, soil media, topography, the impact of vadose zone and hydraulic conductivity) that portrays cumulative aquifer sensitivity ratings in Kishangarh, Rajasthan. It provides a relative indication of groundwater aquifer vulnerability to contamination. The soil moisture flow and solute transport regimes of the vadose zone associated with specific hydrogeological conditions play a crucial role in pollution risk assessment of the underlying groundwater resources. An effort has been made to map the vulnerability of shallow groundwater to surface pollutants of thestudy area, using soil moisture flow and contaminant transport modeling. The classical advection-dispersion equation coupled with Richard’s equation is numerically simulated at different point locations for assessing the intrinsic vulnerability of the valley. The role of soil type, slope, and the land-use cover is considered for estimating the transient flux at the top boundary from daily precipitation and evapotranspiration data of the study area. The time required by the solute peak to travel from the surface to the groundwater table at the bottom of the soil profile is considered as an indicator of avulnerability index. Results show a high vulnerability in the southern region, whereas low vulnerability is observed in the northeast and northern parts. The results have recognized four aquifer vulnerability zones based on DRASTIC vulnerability index (DVI), which ranged from 45 to 178. It has been deduced that approximately 18, 25, 34, and 23% of the area lies in negligible, low, medium and high vulnerability zones, respectively. The study may assist in decision making related to theplanning of industrial locations and the sustainable water resources development of the selected semi-arid area.  相似文献   

12.
This paper focuses on artificial groundwater recharge study in Ayyar basin, Tamil Nadu, India. The basin is covered by hard crystalline rock and overall has poor groundwater conditions. Hence, an artificial recharge study was carried out in this region through a project sponsored by Tamil Nadu State Council for Science and Technology. The Indian Remote Sensing satellite 1A Linear Imaging Self Scanning Sensor II (IRS 1A LISS II) satellite imagery, aerial photographs and geophysical resistivity data were used to prioritize suitable sites for artificial recharge and to estimate the volume of aquifer dimension available to recharge. The runoff water available for artificial recharge in the basin is estimated through Soil Conservation Service curve number method. The land use/land cover, hydrological soil group and storm rainfall data in different watershed areas were used to calculate the runoff in the watersheds. The weighted curve number for each watershed is obtained through spatial intersection of land use/land cover and hydrological soil group through GeoMedia 3.0 Professional GIS software. Artificial recharge planning was derived on the basis of availability of runoff, aquifer dimension, priority areas and water table conditions in different watersheds in the basin.  相似文献   

13.
Gradients controlling natural attenuation of ammonium   总被引：1，自引：0，他引：1  
Oxidation of reduced pollutants such as in groundwater often takes place at steep redox gradients where oxygenated water is being mixed into polluted water such as landfill leachate. In order to identify controlling parameters and quantify the influence of environmental factors for degradation, sensitivity analysis was performed by means of scenario specific numerical modelling. Geometrical factors such as aquifer thickness have been shown to be very influential on the capability of natural attenuation of pollutants in groundwater. The scenarios investigated here include biodegradation at redox gradients in groundwater, so called fringe processes, for (i) a partly contaminated aquifer with two reaction fronts, (ii) and a spatially variable aquifer thickness. In addition, (iii) the influence of groundwater recharge and (iv) restricted supply of O2 to contaminated water by slow dispersion and diffusion across the capillary fringe are investigated. Contaminated aquifer thickness, zones of enhanced mixing due to flow focussing and diffusion/dispersion coefficients in the capillary fringe are identified qualitatively as controlling factors for natural attenuation under complex conditions, whereas predictive functions will require further research.  相似文献   

14.
  相似文献   

15.
The sulfate pollution in an agriculturally used watershed has been investigated with respect to the transport in the saturated zone and the development of sulfate in the unsaturated zone. Besides of other sources such as acid wet and dry deposition or sulfate input by agricultural activities, most of the sulfate originates from oxidation of pyrite by either NO3 or O2. High sulfate concentrations coincided with high nitrate leaching caused by plowing of former grassland or by vegetable crop residues and with former wet lands that have become dry. By using soil water concentration data and maps showing the extension of former wetlands and grassland as well as agricultural land use, it was possible to delineate regions of high sulfate input. The transport of sulfate in the aquifer was analyzed with a modified version of the USGS MOC model, which takes into account the nonlinearity of the underlying equation describing unconfined groundwater flow. The calibration of the transport model showed good agreement between the estimated and modeled sulfate input rates. A prediction of future sulfate concentrations in the aquifer was feasible by using worst-case parameters.  相似文献   

16.
Starting with a stochastic differential equation with random coefficients describing steady-state flow, the effective hydraulic conductivity of 1-, 2-, and 3-dimensional aquifers is derived. The natural logarithm of hydraulic conductivity (lnK) is assumed to be heterogeneous, with a spatial trend, and isotropic. The effective conductivity relates the mean specific discharge in an aquifer to the mean hydraulic gradient, thus its importance in predicting Darcian discharge when field data represent mean or average values of conductivity or hydraulic head. Effective conductivity results are presented in exact form in terms of elementary functions after the introduction of special sets of coordinate transformations in two and three dimensions. It was determined that in one, two, and three dimensions, for the type of aquifer heterogeneity considered, the effective hydraulic conductivity depends on: (i) the angle between the gradient of the trend of lnK and the mean hydraulic gradient (which is zero in the one-dimensional situation); (2) (inversely) on the product of the magnitude of the trend gradient of lnK, b, and the correlation scale of lnK, and (3) (proportionally) on the variance of lnK, f 2 . The productb plays a central role in the stability of the results for effective hydraulic conductivity.  相似文献   

17.
The High Plains aquifer in the Southern High Plains (Texas and New Mexico), consisting of Tertiary, Cretaceous, and Triassic formations, has traditionally been considered to be recharged by its uppermost water-bearing unit, the Tertiary Ogallala aquifer. This article provides hydrologic, chemical, and isotopic evidence that in the Southern High Plains: (1) Cretaceous rocks actually contain independent recharge sources; (2) Triassic rocks cannot currently be recharged by the Ogallala aquifer in significant quantities; and (3) in places, both Cretaceous and Triassic aquifers recharge the overlying Ogallala aquifer. On the basis of chemical and isotopic data, playa lakes seem to act as the predominant recharge source of the Ogallala aquifer, suggesting recharge rates greater than 30 mm/yr, as opposed to the much lower rates reported by others. The Cretaceous aquifers are being recharged by cross-formational flow from the Ogallala aquifer but also from overlying Quaternary sands and the underlying Triassic aquifer in eastern New Mexico. Current recharge into the Triassic aquifer may be insignificant.  相似文献   

18.
Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.  相似文献   

19.
Groundwater is the major source of water supply for most uses in the rural settlements in Ghana. A groundwater flow model was calibrated for some aquifers of the southern Voltaian sedimentary system under steady-state conditions. The objective was to determine estimates of the hydraulic conductivities of the different hydrostratigraphic units of the southern Voltaian, and the distribution of recharge from precipitation. Data on the stable isotopes of oxygen and hydrogen from the study area suggest that groundwater recharge in the area is from rainfall. The calibrated steady-state model suggests that aquifer hydraulic conductivities in the study area range from 1.19 to 6.3 m/day. Hydrostratigraphic unit specific hydraulic conductivities are discussed. The calibrated recharge ranges from 3.81e−05 m/day to 6.0e−04 m/day, which represents 0.9–6% of the precipitation in the form of rainfall. Six distinct flowpaths have been defined using particle tracking. The particle tracking simulation suggests travel times in the range of 380 to 5,199 years from recharge areas to discharge areas along the flowpaths identified. A contaminant dropped at the recharge areas in the central sections of the model area would travel at these rates along the flowpaths, assuming that advection is the dominant transport process. Inverse geochemical modeling indicates the dissolution of albite, K-feldspars and anorthite, respectively, along flowpaths I and IV. The inverse modeling along flowpaths I and IV suggest the dissolution of albite, K-feldspar and anorthite, respectively, at 1.085e−06, 3.16e−08 and 3.168e−07 mmol/year.  相似文献   

20.
A modeling study of seawater intrusion in Alabama Gulf Coast,USA   总被引：4，自引：0，他引：4  
A numerical model of variable-density groundwater flow and miscible salt transport is developed to investigate the extent of seawater intrusion in the Gulf coast aquifers of Alabama, USA. The SEAWAT code is used to solve the density-dependent groundwater flow and solute transport governing equations. The numerical model is calibrated against the observed hydraulic heads measured in 1996 by adjusting the zonation and values of hydraulic conductivity and recharge rate. Using the calibrated model and assuming all the hydrogeologic conditions remain the same as those in 1996, a predictive 40-year simulation run indicates that further seawater intrusion into the coastal aquifers can occur in the study area. Moreover, the predicted intrusion may be more significant in the deeper aquifer than the shallower ones. As the population continues to grow and the demand for groundwater pumping intensifies beyond the 1996 level, it can be expected that the actual extent of seawater intrusion in the future would be more severe than the model prediction. Better strategies for groundwater development and management will be necessary to protect the freshwater aquifers from contamination by seawater intrusion.
Jin LinEmail:
  相似文献   

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

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