首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
为了及时有效地应对各种突发性环境污染事故,有必要开发一种简单实用、适于各类型污染物的场地污染数学模型。通过污染事故发生后污染物在包气带、饱和带迁移转化的概化,建立了污染物运移的自由入渗模型以及降雨入渗模型并给出各自相应的解析解。无降雨时,考虑污染物在重力作用下随包气带向下渗透的作用,建立一维垂直入渗模型。有降雨时,考虑污染场地（包气带）中污染物迁移和转化的对流作用、扩散作用及挥发、生物降解、吸附、根系吸收等作用,建立包气带剖面二维溶质运移模型和饱和带平面二维溶质运移数学模型。建模过程中,假定降雨量的平均分布及土壤质地、水力参数以及有机物成分、种类均相同,同时假定污染物与多孔介质间的作用为线性吸附,植物根系对污染物的吸附遵循一级动力学。基于模型的解析解,实现案例的模拟计算。模拟结果表明：该模型具有适用范围广、模拟高效快捷等优点,能够较准确预测污染发生后污染物在土壤中的动向、到达饱和带的时间以及饱和带中污染物的迁移情况。  相似文献   

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

3.
The principal aim of this study is to assess the scope of monitoring diesel plume migration in a scaled aquifer model with a miniaturised electrical resistivity array. Respectively 1000 and 500 ml of diesel were injected in both the unsaturated and water-saturated zones of a sand body overlying a clay aquitard, and diesel migration was monitored with a miniature electrode array and an off-the-shelf resistivity meter. Inverted time-lapse electrical resistivity tomography (ERT) data reflect downward and lateral spreading of the diesel plume away from the injection point in the unsaturated zone. Diesel was also imaged to spread upwards and laterally away from the injection point in the saturated zone, as controlled by capillary rise. In both cases later-time ERT images reflected preferential pooling of diesel on the water table, as well as vertical smearing of pooled diesel in response to simulated water-table fluctuations. Repeat fluid electrical conductivity (EC) and dissolved oxygen (DO) measurements validate the observed changes in bulk resistivity caused by both diesel injections. Artefacts introduced by 2D inversion of 3D contaminant transport were abound. Time-lapse ERT imaging of diesel transport is therefore inferred to be feasible and well-suited to complementing conventional techniques of intrusive site investigation, although time-lapse 3D or 4D ERT imaging is strongly advocated.  相似文献   

4.
Two geophysical techniques were integrated to map the groundwater aquifers on complex geological settings, in the crystalline basement terrain in northeast Nuba Mountains. The water flow is structurally controlled by the northwest–southeast extensional faults as one of several in-situ deformational patterns that are attributed to the collision of the Pan-African oceanic assemblage of the Nubian shield against the pre-Pan African continental crust to the west. The structural lineaments and drainage systems have been enhanced by the remote sensing technique. The geophysical techniques used are: vertical electrical soundings (VES) and electrical resistivity tomography (ERT), in addition to hydraulic conductivity measurements. These measurements were designed to overlap in order to improve the producibility of the geophysical data and to provide a better interpretation of the hydrogeological setting of the aquifer complex structure. Smooth and Block inversion schemes were attempted for the observed ERT data to study their reliability in mapping the different geometries in the complex subsurface. The VES data was conducted where ERT survey was not accessible, and inverted smoothly and merged with the ERT in the 3D resistivity grid. The hydraulic conductivity was measured for 42 water samples collected from the distributed dug wells in the study area; where extremely high saline zones were recorded and have been compared to the resistivity values in the 3D model.  相似文献   

5.
A series of geophysical parameters have been applied with geological perception to resolve the hydrogeological complexities over granitic terrain at Hyderabad, India. Frequent failure of borehole drillings and the thrust conditions of community have prompted a noninvasive suitable tool, applied at small scale for pinpointing potential well site. Geophysical scanning, viz. electrical resistivity tomography (ERT), spontaneous potential (SP), and electrical gradient profiling (GP) were employed within the restricted space of housing complex to obtain the true characteristics of the subsurface lithology, where anomalies by the underground utility structures have been nullified. Results showed, in ERT, the low order of resistivity range 123 to 200?Ωm showing a plume like weathered zone underlain by sudden slip of massive granite (>217?Ωm) was of great interest to proceed further in the process. Here, the switchover in SP value from +18 mV to ?17 mV and GP from mean value 10 to 90 mV/m was recorded. The anomalies in SP and GP were precisely coincided with the ERT where upcoming of subsurface massive granite next to the inferred fracture was noted. Drilling core logs satisfies the geophysical signatures ensuring the inferred saturated fracture with the total yield 1,302 gal/h.  相似文献   

6.
Geophysical methods have been applied to a wide range of hydrogeological problems. With improvement in geophysical inversion algorithms and measurement tools, significant achievements have been made in the characterization of subsurface architecture, time-lapse monitoring of hydrogeological process and contaminant plumes delineation. In this paper, we summarize the geophysical methods that are most widely used in hydrogeology including Electrical Resistivity Tomography (ERT), Induced Polarization (IP), Ground Penetrating Radar (GPR) and Electromagnetic Induction (EMI). Three examples including lab and field works are used to demonstrate current application of geophysical methods for characterizing subsurface architecture and contaminant plumes. Though great progress has been made in hydrogeohysics over the last few decades at home and abroad, challenges still remain in practical applications. More recently, hydrogeophysics continues to develop in the areas of establishment of hydrogeophysical models, large-scale architecture characterization, uncertainty analysis, biogeochemical process monitoring and ecosystem science.  相似文献   

7.
Monitoring transport of dissolved substances in soil deposits is particularly relevant where safety is concerned, as in the case of geo-environmental barriers. Geophysical methods are very appealing, since they cover a wide domain, localising possible preferential flow paths and providing reliable links between geophysical quantities and hydrological variables. This paper describes a 3D laboratory application of electrical resistivity tomography (ERT) used to monitor solute transport processes. Dissolution and transport tests on both homogeneous and heterogeneous samples were conducted in an instrumented oedometer cell. ERT was used to create maps of electrical conductivity of the monitored domain at different time intervals and to estimate concentration variations within the interstitial fluid. Comparisons with finite element simulations of the transport processes were performed to check the consistency of the results. Tests confirmed that the technique can monitor salt transport, infer the hydro-chemical behaviour of heterogeneous geomaterials and evaluate the performances of clay barriers.  相似文献   

8.
含水层非均质性的刻画是模拟地下水中污染物运移的关键。以渗透系数为研究对象,构建了综合集合卡尔曼滤波方法、有效电阻率模型与地下水运移模型的同化框架,通过融合地球物理观测数据与污染物浓度观测数据来推估渗透系数的空间分布。基于理想算例,验证了该同化框架刻画含水层非均质渗透系数场的有效性,并针对不同初始参数信息与观测类型对比了耦合与非耦合水文地球物理方法的适用性。研究结果表明:基于集合卡尔曼滤波方法同化多种类型的观测数据,可有效地推估非均质参数空间分布。当初始信息较准确时,耦合方法的参数推估精度更高;初始信息存在偏差时,非耦合方法有更好的同化效果。由于非耦合方法计算成本较低且对初始信息缺失时适用性更强,在实际应用中可先基于非耦合方法初步估计参数,再利用耦合方法进一步提高参数推估精度。融合多种类型观测数据可有效提高参数推估效果。  相似文献   

9.
Self potential (SP) and electrical resistivity tomography (ERT) methods are used together with the results of groundwater samples hydrogeochemical analysis to assess the impact of the water leak from the landfill garbage site at NamSon located in Northern Hanoi on causing pollution to the surrounding environment and affecting geological structure. Selected survey area covers an area of 180 × 300 m lying in the low land of the NamSon site with a slope ranging about 8 m in direction NW–SE. There are three geophysical measurements lines denoted as T1, T2 and T3. Processing 180 SP data points has allowed to draw maps of equipotential epoch in the two periods in 2015 and 2016. The maps show four zones of SP positive anomalies with maximum amplitudes of about +20 mV where the groundwater flow direction is downward and five zones of SP negative anomalies with minimum values in a range from ?180 to ?260 mV where the groundwater flow direction is upward. Resistivity values of the subsurface layers of soils and rocks have been aquired from 2D inverse model for measuring ERT in March 2015 and March 2016. The results of the ERT allowed to define the low resistivity in the range 15–20 Ωm related to leachate plume from NamSon landfill site. Results of the physico-chemical analysis of groundwater samples from the existing six boreholes show increases in concentration of the measured pollutant parameters indicating contamination of the groundwater as a result of solid waste leachate accumulation. This result is affirmative evidence for the survey results by geophysical technique. The rapid decrease in quality of groundwater over the last year is probably due to the influence of the leachate from the NamSon landfill site.  相似文献   

10.
Precise and efficient numerical simulation of transport processes in subsurface systems is a prerequisite for many site investigation or remediation studies. Random walk particle tracking (RWPT) methods have been introduced in the past to overcome numerical difficulties when simulating propagation processes in porous media such as advection-dominated mass transport. Crucial for the precision of RWPT methods is the accuracy of the numerically calculated ground water velocity field. In this paper, a global node-based method for velocity calculation is used, which was originally proposed by Yeh (Water Resour Res 7:1216–1225, 1981). This method is improved in three ways: (1) extension to unstructured grids, (2) significant enhancement of computational efficiency, and (3) extension to saturated (groundwater) as well as unsaturated systems (soil water). The novel RWPT method is tested with numerical benchmark examples from the literature and used in two field scale applications of contaminant transport in saturated and unsaturated ground water. To evaluate advective transport of the model, the accuracy of the velocity field is demonstrated by comparing several published results of particle pathlines or streamlines. Given the chosen test problem, the global node-based velocity estimation is found to be as accurate as the CK method (Cordes and Kinzelbach in Water Resour Res 28(11):2903–2911, 1992) but less accurate than the mixed or mixed-hybrid finite element methods for flow in highly heterogeneous media. To evaluate advective–diffusive transport, a transport problem studied by Hassan and Mohamed (J Hydrol 275(3–4):242–260, 2003) is investigated here and evaluated using different numbers of particles. The results indicate that the number of particles required for the given problem is decreased using the proposed method by about two orders of magnitude without losing accuracy of the concentration contours as compared to the published numbers.  相似文献   

11.
Polychlorinated biphenyls (PCBs) are the main constituents of clophen (the liquid of the electric transformers and capacitors) and have been characterized as potential human carcinogens. PCBs can be a hazardous contaminant of soil and groundwater. We used the mathematical model variably saturated 2D flow and transport (VS2DT model) to simulate the transport of PCBs from the soil surface to groundwater for a time period of 30 years. We also used a mathematical model to simulate the colloid-facilitated PCB transport, under saturated flow conditions. The results showed that PCBs dissolved in water cannot be transported to large depths in unsaturated soils, because of their strong sorption onto soil and low solubility in water. For soils with very low or no organic matter content, PCB transport is much faster and the probability of groundwater contamination is much higher. PCBs can partition to colloids originating from dissolved organic matter in groundwater. Colloid-facilitated PCB transport is faster compared to PCB transport in aqueous solution with no colloids present.  相似文献   

12.
Geoelectrical methods involving electrical resistivity tomography (ERT), self-potential (SP), frequency domain electromagnetic (FDEM), and very low frequency (VLF) methods have been used to provide valuable information in locating a known sulfide ore body in Soap Gulch, Montana. The study develops basis of comparison for the geophysical techniques employed. Ranges of resistivity along the area have been established using interpreted ERT which can help to understand the subsurface distribution of sulfides in the area. A sulfide body was delineated from the survey area corresponding to anomalously low resistivity values on the ERT section, negative SP, and high apparent current density zone in VLF. Depth to the localized ore zone ranges approximately from 10 to 20 m. FDEM data reflect the conductivity distribution of the shallow subsurface (less than 6 m deep); hence, the delineated sulfide zone had minimal contribution to FDEM measurements. The results of the study show that SP, VLF, and ERT methods provide significant information in localizing ore bodies. The survey revealed that the resistivity values obtained from ERT profile corroborate the FDEM, SP, and VLF from the area.  相似文献   

13.
This paper considers the hydrogeological simulation of groundwater movement in karstic regions using a hydrological modelling system (SHETRAN) which has been adapted for modelling flow in karstic aquifers. Flow and transport through karstic aquifers remains poorly understood, yet quantitative hydrogeological models are essential for developing and implementing groundwater protection policies. The new model has been developed and used within the STALAGMITE (Sustainable Management of Groundwater in Karstic Environments) project, funded by the European Commission. The SHETRAN model is physically based insofar as most of the parameters have some physical meaning. The SHETRAN model represents all of the key processes in the hydrological cycle, including subsurface flow in the saturated and unsaturated zones, surface flow over the ground surface and in channels, rainfall interception by vegetation canopies, evapotranspiration, snow-pack development and snowmelt. The modifications made to SHETRAN to simulate karstic aquifers are (1) the coupling of a pipe network model to a variably saturated, three-dimensional groundwater component (the VSS-NET component), to simulate flow under pressure in saturated conduits; (2) the coupling of surface water features (e.g. sinking streams or "ponors", and spring discharges) to the conduit system; (3) the addition of a preferential "bypass" flow mechanism to represent vertical infiltration through a high-conductivity epikarst zone. Lastly, a forward particle tracking routine has been developed to trace the path of hypothetical particles with matrix and pipe flow to springs or other discharge points. This component allows the definition of groundwater protection zones around a source for areas of the catchment (watershed) which are vulnerable to pollution from non-point sources (agriculture and forestry).  相似文献   

14.
Characterizing percolation patterns in unsaturated fractured rock has posed a greater challenge to modeling investigations than comparable saturated zone studies due to the heterogeneous nature of unsaturated media and the great number of variables impacting unsaturated flow. An integrated modeling methodology has been developed for quantitatively characterizing percolation patterns in the unsaturated zone of Yucca Mountain, Nevada (USA), a proposed underground repository site for storing high-level radioactive waste. The approach integrates moisture, pneumatic, thermal, and isotopic geochemical field data into a comprehensive three-dimensional numerical model for analyses. It takes into account the coupled processes of fluid and heat flow and chemical isotopic transport in Yucca Mountain’s highly heterogeneous, unsaturated fractured tuffs. Modeling results are examined against different types of field-measured data and then used to evaluate different hydrogeological conceptualizations through analyzing flow patterns in the unsaturated zone. In particular, this model provides clearer understanding of percolation patterns and flow behavior through the unsaturated zone, both crucial issues in assessing repository performance. The integrated approach for quantifying Yucca Mountain’s flow system is demonstrated to provide a practical modeling tool for characterizing flow and transport processes in complex subsurface systems.  相似文献   

15.
This paper proposes a two-stage geophysical approach to map the vertical cracking and the structural integrity of flood embankments made up of clay geomaterials susceptible to fissuring. The first stage is based on a ‘coarse-resolution’ investigation using conventional electrical resistivity tomography (ERT) equipment to identify the fissured zones in the embankment. This step is complemented by an additional geophysical technique, electromagnetic, to verify the ERT measurements. The second stage is based on a ‘high-resolution’ investigation using a miniature ERT system previously developed at the laboratory scale for detailed mapping of the fissure patterns. The ‘coarse-resolution’ stage is the major focus of this paper and was validated against two case studies in England and Scotland. Longitudinal ERT survey provided a tomographic picture of the upper desiccated zones of the embankments and fissured areas in 2-D, validating the range of resistivity results obtained previously on a fissured clay model in the laboratory. A transversal embankment resistivity tomography was also completed to show the positions of fissured zones in detail in the field. The electromagnetic technique as a fast screening tool allowed cross checking the ERT results and was also efficient in detecting high and low conductivity zones, indicating areas of potential weakness during flash floods and heavy rain. The southern embankment in England showed more fluctuations in the conductivity and resistivity than the north embankment in Scotland, likely to be due to the differences in climate, vegetation and location characteristics between the two sites. Conclusions were also drawn on the potential weaknesses for both embankments and the effect of vegetation on conductivity measurements.  相似文献   

16.
The movement of chemicals through soil to groundwater is a major cause of degradation of water resources. In many cases, serious human and stock health implications are associated with this form of pollution. The study of the effects of different factors involved in transport phenomena can provide valuable information to find the best remediation approaches. Numerical models are increasingly being used for predicting or analyzing solute transport processes in soils and groundwater. This article presents the development of a stochastic finite element model for the simulation of contaminant transport through soils with the main focus being on the incorporation of the effects of soil heterogeneity in the model. The governing equations of contaminant transport are presented. The mathematical framework and the numerical implementation of the model are described. The comparison of the results obtained from the developed stochastic model with those obtained from a deterministic method and some experimental results shows that the stochastic model is capable of predicting the transport of solutes in unsaturated soil with higher accuracy than deterministic one. The importance of the consideration of the effects of soil heterogeneity on contaminant fate is highlighted through a sensitivity analysis regarding the variance of saturated hydraulic conductivity as an index of soil heterogeneity. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

17.
The unsaturated fractured volcanic deposits at Yucca Mountain in Nevada, USA, have been intensively investigated as a possible repository site for storing high-level radioactive waste. Field studies at the site have revealed that there exist large variabilities in hydrological parameters over the spatial domain of the mountain. Systematic analyses of hydrological parameters using a site-scale three-dimensional unsaturated zone (UZ) flow model have been undertaken. The main objective of the sensitivity analyses was to evaluate the effects of uncertainties in hydrologic parameters on modeled UZ flow and contaminant transport results. Sensitivity analyses were carried out relative to fracture and matrix permeability and capillary strength (van Genuchten α) through variation of these parameter values by one standard deviation from the base-case values. The parameter variation resulted in eight parameter sets. Modeling results for the eight UZ flow sensitivity cases have been compared with field observed data and simulation results from the base-case model. The effects of parameter uncertainties on the flow fields were evaluated through comparison of results for flow and transport. In general, this study shows that uncertainties in matrix parameters cause larger uncertainty in simulated moisture flux than corresponding uncertainties in fracture properties for unsaturated flow through heterogeneous fractured rock.  相似文献   

18.
Subsurface geophysical surveys were carried out using a large range of methods in an unconfined sandstone aquifer in semiarid south-western Niger for improving both the conceptual model of water flow through the unsaturated zone and the parameterization of numerical a groundwater model of the aquifer. Methods included: electromagnetic mapping, electrical resistivity tomography (ERT), resistivity logging, time domain electromagnetic sounding (TDEM), and magnetic resonance sounding (MRS). Analyses of electrical conductivities, complemented by geochemical measurements, allowed us to identify preferential pathways for infiltration and drainage beneath gullies and alluvial fans. The mean water content estimated by MRS (13%) was used for computing the regional groundwater recharge from long-term change in the water table. The ranges in permeability and water content obtained with MRS allowed a reduction of the degree of freedom of aquifer parameters used in groundwater modelling.  相似文献   

19.
The impact of diffuse pollution, agricultural land use and climate change on the long-term response of subsurface–surface water quality is not well understood, but is a prerequisite for evaluation of water management options. The goal of this study is to model geochemical evolution of water chemistry from the infiltration through soil into the unsaturated zone, transport through bedrocks and granular aquifers to a river in order to identify zones of steep concentration gradients and high dynamics under transient flow conditions. A numerical model was constructed comprising a 2-D 1,500 m × 150 m vertical cross-section of typical sedimentary rock formations, a glacio-fluvial quaternary gravel aquifer in the valley and soil layers. The model coupled saturated/un-saturated flow and reactive transport under steady state and transient conditions. Geochemical interactions, include intra-aqueous kinetic reactions of oxygen with dissolved organic matter, as well as kinetics of carbonate dissolution/precipitation. This model section was chosen to provide insight in to the principal processes and time scales affecting water chemistry along different flow paths. The numerical simulator MIN3P was used, a finite volume program for variably saturated subsurface flow and multi-component reactive transport. The results show that subsurface water residence times range from approximately 2 to 2,000 years. Different zones are to be expected with respect to the development of mineral equilibria; namely, purely atmospherically influenced, as well as open and closed system carbonate dissolution. Short-term responses to daily averaged changes in precipitation, however, are only visible to some extent in the shallower and near-river parts of flow system and solute loads. This can most likely be explained by directional changes in flow paths, indicating that equilibrium geochemical condition predominate at the hillslope scale, i.e. water quality depends on transport pathways rather than on kinetic effects. The extent of reducing conditions is controlled by the presence of organic-rich layers (i.e. peat deposits), the dissolution kinetics of aquifer organic matter and the subsequent mixing with oxygenated water by hydrodynamic dispersion.  相似文献   

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

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

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