首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 156 毫秒
1.
This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater.  相似文献   

2.
Several laboratory experiments were conducted to identify the validity domain under which a Hele–Shaw cell may serve as a suitable analogue for variable-density flow in homogeneous porous media. These experiments are concerned with the injection into a Hele–Shaw cell of a salt solute at different concentrations and flow rates. The experimental data analysis highlighted two types of mixing zone shape: with and without ‘fingers’. A semi-empirical criterion based on the ratio between gravitational and injected velocities was used to forecast the change from one shape to another. The experimental data were then analysed using numerical solutions of the classical Hele–Shaw equations by taking into account an anisotropic dispersion tensor whose components depend on fluid density gradients. The good agreement between experimental and numerical results clearly shows that the validity of the concentration-dependent dispersion tensor strongly depends on the local Péclet number variation. For Péclet numbers lower 50, the Hele–Shaw cell can be considered as an analogous model of a homogeneous and isotropic 2D porous medium. It can be successfully used to study, at the laboratory scale, the gravitational instability effects induced by flow and transport phenomena into a porous medium.  相似文献   

3.
Optimal and sustainable extraction of groundwater in coastal aquifers   总被引：1，自引：0，他引：1  
Four examples are investigated for the optimal and sustainable extraction of groundwater from a coastal aquifer under the threat of seawater intrusion. The objectives and constraints of these management scenarios include maximizing the total volume of water pumped, maximizing the profit of selling water, minimizing the operational and water treatment costs, minimizing the salt concentration of the pumped water, and controlling the drawdown limits. The physical model is based on the density-dependent advective-dispersive solute transport model. Genetic algorithm is used as the optimization tool. The models are tested on a hypothetical confined aquifer with four pumping wells located at various depths. These solutions establish the feasibility of simulating various management scenarios under complex three-dimensional flow and transport processes in coastal aquifers for the optimal and sustainable use of groundwater.  相似文献   

4.
Analytical models have been exhaustively used to study simple seawater intrusion problems and the sustainable management of groundwater resources in coastal aquifers because of its simplicity, easy implementation, and low computational cost. Most of these models are based on the sharp‐interface approximation and the Ghyben–Herzberg relation, and their governing equations are expressed in terms of a single potential theory to calculate critical pumping rates in a coastal pumping scenario. The Ghyben–Herzberg approach neglects mixing of fresh water and seawater and implicitly assumes that salt water remains static. Therefore, the results of the analytical solutions may be inaccurate and unacceptable for some real‐complex case studies. This paper provides insight into the validity of sharp‐interface models to deal with seawater intrusion in coastal aquifers, i.e. when they can be applied to obtain accurate enough results. For that purpose, this work compares sharp‐interface solutions, based on the Ghyben–Herzberg approach, with numerical three‐dimensional variable‐density flow simulations for a set of heterogeneous groundwater flow and mass transport parameters, and different scenarios of spatially distributed recharge values and spatial wells placement. The numerical experiment has been carried out in a 3D unconfined synthetic aquifer using the finite difference numerical code SEAWAT for solving the coupled partial differential equations of flow and density‐dependent transport. This paper finds under which situations the sharp‐interface solution gives good predictions in terms of seawater penetration, transition zone width and critical pumping rates. Additionally, the simulation runs indicate to which parameters and scenarios the results are more sensitive. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

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

6.
  相似文献   

7.
Groundwater is a primary source of drinking water worldwide, but excess nutrients and emerging contaminants could compromise groundwater quality and limit its usage as a drinking water source. As such contaminants become increasingly prevalent in the biosphere, a fundamental understanding of their fate and transport in groundwater systems is necessary to implement successful remediation strategies. The dynamics of surface water-groundwater (hyporheic) exchange within a glacial, buried-valley aquifer system are examined in the context of their implications for the transport of nutrients and contaminants in riparian sediments. High conductivity facies act as preferential flow pathways which enhance nutrient and contaminant delivery, especially during storm events, but transport throughout the aquifer also depends on subsurface sedimentary architecture (e.g. interbedded high and low conductivity facies). Temperature and specific conductance measurements indicate extensive hyporheic mixing close to the river channel, but surface water influence was also observed far from the stream-aquifer interface. Measurements of river stage and hydraulic head indicate that significant flows during storms (i.e., hot moments) alter groundwater flow patterns, even between consecutive storm events, as riverbed conductivity and, more importantly, the hydraulic connectivity between the river and aquifer change. Given the similar mass transport characteristics among buried-valley aquifers, these findings are likely representative of glacial aquifer systems worldwide. Our results suggest that water resources management decisions based on average (base) flow conditions may inaccurately represent the system being evaluated, and could reduce the effectiveness of remediation strategies for nutrients and emerging contaminants.  相似文献   

8.
Patterns and Rates of Ground-Water Flow on Long Island, New York   总被引：3，自引：0，他引：3  
Increased ground-water contamination from human activities on Long Island has prompted studies to define the pattern and rate of ground-water movement. A two-dimensional, fine-mesh, finite-element model consisting of 11,969 nodes and 22,880 elements was constructed to represent ground-water flow along a north-south section through central Long Island. The model represents average hydrologic conditions within a corridor approximately 15 miles wide. The model solves discrete approximations of both the potential and stream functions. The resulting flownet depicts flow paths and defines the vertical distribution of flow within the section. Ground-water flow rates decrease with depth. Sixty-two percent of the water flows no deeper than the upper glacial (water-table) aquifer, 38 percent enters the underlying Magothy aquifer, and only 3.1 percent enters the Lloyd aquifer. The limiting streamlines for flow to the Magothy and Lloyd aquifers indicate that aquifer recharge areas are narrow east-west bands through the center of the island. The recharge area of the Magothy aquifer is only 5.4 miles wide; that of the Lloyd aquifer is less than 0.5 miles. The distribution of ground-water traveltime and a flownet are calculated from model results; both are useful in the investigation of contaminant transport or the chemical evolution of ground water within the flow system. A major discontinuity in traveltime occurs across the streamline which separates the flow subsystems of the two confined aquifers. Water that reaches the Lloyd aquifer attains traveltimes as high as 10,000 years, whereas water that has not penetrated deeper than the Magothy aquifer attains traveltimes of only 2,000 years. The finite-element approach used in this study is particularly suited to ground-water systems that have complex hydrostratigraphy and cross-sectional symmetry.  相似文献   

9.
Methods of estimation of the location of a sharp fresh water–salt water interface with hydraulic heads or pressures are relatively simple and are widely used. Progress has been made in the recent decade toward the mathematical relations describing the position of the sharp interface using hydraulic heads or pressures in coastal zones. This paper reviews several methods for estimation of the location of fresh water–salt water interface in coastal aquifers, including the classical Ghyben–Herzberg relation. The location of the fresh water–salt water interface in a coastal homogeneous, isotropic unconfined aquifer can be estimated based on piezometric heads at two points in the same vertical line tapping, respectively, the salt water zone (including the interface) and the fresh water zone (from the water table to the interface) when the groundwater flow system is in a steady state and satisfies the Dupuit assumption. If pressures are measured at two points in the fresh water and salt water zones in the same vertical line in the coastal aquifer under the same assumption, then the position of the interface can still be estimated with the pressure data. If the Dupuit assumption is not met in coastal aquifers and the vertical fresh water head gradients can be approximated with a straight line, the position of the interface can roughly be estimated by using the water level data in a partially penetrating well during drilling of the well.  相似文献   

10.
Volcanic aquifers supply a substantial portion of water resources in many parts of the world, including islands, and their productivity depends strongly on volcanic stratigraphy, which exhibits considerable heterogeneity. We investigated water inflow to lava tube caves formed from numerous basaltic lava flows in the northeastern coastal area of Jeju Island after storm events and monitored relative inflow rates monthly over 1 year to characterize groundwater flow processes in the upper parts of volcanic aquifers, and to evaluate the applicability of the previous hydrogeological models proposed for the island. Considerable water inflow arose shortly after storms from exposed palaeosol layers on the walls of the caves. The monthly monitoring results showed that wall inflow associated with these palaeosol layers is substantial. In both cases, discharge from ceiling drips was much less and more temporally variable compared to wall inflow discharge. Water flowing into the caves was rapidly drained through the floor at all monitoring sites. The lateral extent of the palaeosol layers was identified using drill core logs near the cave and outcrops in the coastal area. Based on these results, we inferred that multiple perched aquifers are formed by low-permeability palaeosol layers between lava flows, which are connected by vertical flows at discontinuities in the palaeosol layer, eventually reaching the basal aquifer. This study revealed the water inflow processes observed in lava tube caves constrained by palaeosol layers, and established a hydrogeological conceptual model incorporating multiple perched aquifers in both coastal and mountainous areas associated with extensive palaeosol layers formed during volcanic hiatuses. This finding would help elucidate recharge, groundwater flow, and contaminant transport processes in many volcanic aquifers that are not adequately represented by the previous models, and contribute to better management of groundwater in those areas.  相似文献   

11.
A two‐dimensional variable‐density groundwater flow and transport model was developed to provide a conceptual understanding of past and future conditions of nitrate (NO3) transport and estimate groundwater nitrate flux to the Gulf of Mexico. Simulation results show that contaminant discharge to the coast decreases as the extent of saltwater intrusion increases. Other natural and/or artificial surface waters such as navigation channels may serve as major sinks for contaminant loading and act to alter expected transport pathways discharging contaminants to other areas. Concentrations of NO3 in the saturated zone were estimated to range between 30 and 160 mg?L?1 as NO3. Relatively high hydraulic vertical gradients and mixing likely play a significant role in the transport processes, enhancing dilution and contaminant migration to depth. Residence times of NO3 in the deeper aquifers vary from 100 (locally) to about 300 years through the investigated aquifer system. NO3 mass fluxes from the shallow aquifers (0 to 5.7 × 104 mg?m?2?day?1) were primarily directed towards the navigation channel, which intersects and captures a portion of the shallow groundwater flow/discharge. Direct NO3 discharge to the sea (i.e. Gulf of Mexico) from the shallow aquifer was very low (0 to 9.0 × 101 mg · m?2?day?1) compared with discharge from the deeper aquifer system (0 to 8.2 × 103 mg?m?2?day?1). Both model‐calibrated and radiocarbon tracer‐determined contaminant flux estimates reveal similar discharge trends, validating the use of the model for density‐dependent flow conditions. The modelling approach shows promise to evaluate contaminant and nutrient loading for similar coastal regions worldwide. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

12.
The proper management of coastal aquifers commonly requires an understanding of regional mass flow and complete seawater–freshwater circulation. In this study, time series observations of seawater intrusion and refreshing were conducted using a column experiment based on natural flow conditions in coastal groundwater and a sampled medium from a coastal sandy aquifer without chemical treatment. Ranges of hydrodynamic and hydrochemical variables were tested and analysed. The results showed that the zeta potential of suspended colloids in aqueous solution in an aquifer polluted with 0.5 g/kg of heavy metals exhibited an isoelectric point for pH values ranging from 5.70 to 6.07 when freshwater or seawater completely occupied the aquifer pores, which is representative of natural hydrochemical conditions. In this scenario, a high background concentration of heavy metals induced colloidal immobilization. Otherwise, seawater–freshwater circulation enabled colloid mobilization due to ionic strength and pH fluctuations. The migration of multiple heavy metals occurred at a characteristic time of approximately 1 pore volume after each intrusion stage began and when the peak rate of colloid release was reached. At these times, the colloid behaviour determined the quantity and pathway of heavy metal transport. On the basis of the influences of seawater and freshwater interactions, the quantity of mobilized particles generally decreased and was uniformly distributed in each fraction due to particle loss and decreased porous connectivity. We speculate that the decrease in the total surface area of the migratory colloids may cause colloid‐associated heavy metal transport to decrease. The experimental results provide a useful basis for testing coastal groundwater flow and mass transport models because these phenomena require full characterization to precisely evaluate the associated fluxes from the field scale to the microscopic dimension.  相似文献   

13.
Karst aquifers are highly productive groundwater systems often associated with conduit flow. These systems can be highly vulnerable to contamination, resulting in a high potential for contaminant exposure to humans and ecosystems. This work develops statistical models to spatially characterize flow and transport patterns in karstified limestone and determines the effect of aquifer flow rates on these patterns. A laboratory‐scale Geo‐HydroBed model is used to simulate flow and transport processes in a karstic limestone unit. The model consists of stainless steel tanks containing a karstified limestone block collected from a karst aquifer formation in northern Puerto Rico. Experimental work involves making a series of flow and tracer injections, while monitoring hydraulic and tracer response spatially and temporally. Statistical mixed models (SMMs) are applied to hydraulic data to determine likely pathways of preferential flow in the limestone units. The models indicate a highly heterogeneous system with dominant, flow‐dependent preferential flow regions. Results indicate that regions of preferential flow tend to expand at higher groundwater flow rates, suggesting a greater volume of the system being flushed by flowing water at higher rates. Spatial and temporal distribution of tracer concentrations indicates the presence of conduit‐like and diffuse flow transport in the system, supporting the notion of both combined transport mechanisms in the limestone unit. The temporal response of tracer concentrations at different locations in the model coincide with, and confirms the preferential flow distribution generated with the SMMs used in the study.  相似文献   

14.
The Henry problem has played a key role in our understanding of seawater intrusion into coastal aquifers and in benchmarking density dependent flow codes. This paper seeks to modify Henry’s problem to ensure sensitivity to density variations and vertical salinity profiles that resemble field observations. In the proposed problem, the “dispersive Henry problem”, mixing is represented by means of the traditional Scheidegger dispersion tensor (dispersivity times water flux). Anisotropy in the hydraulic conductivity is acknowledged and Henry’s seaside boundary condition of prescribed salt concentration is replaced by a flux dependent boundary condition, which represents more realistically salt transport across the seaside boundary. This problem turns out to be very sensitive to density variations and its solution gets closer to reality. However, an improvement in the traditional Henry problem (gain in sensitivity and realism) can be also achieved if the value of the Peclet number is significantly reduced.Although the dispersive problem lacks an analytical solution, it can shed light on flow in coastal aquifers. It provides significant information about the factors controlling seawater penetration, width of the mixing zone and influx of seawater. The width of the mixing zone depends basically on dispersion with longitudinal and transverse dispersion controlling different parts of the mixing zone but displaying similar overall effects. Toe penetration is mainly controlled by the horizontal permeability and by the geometric mean of the dispersivities. Finally, transverse dispersivity and the geometric mean of the hydraulic conductivity are the leading parameters controlling the amount of saltwater that enters the aquifer.  相似文献   

15.
The high‐density Dead Sea water (1.235 g/cm3) forms a special interface configuration with the fresh groundwater resources of its surrounding aquifers. The fresh groundwater column beneath its surroundings is around one tenth of its length compared to oceanic water. This fact alone indicates the vulnerability of the fresh groundwater resources to the impacts of changes in the Dead Sea level and to saltwater migration. Ghyben‐Herzberg and Glover equations were used to calculate the volumes of water in coastal aquifers which were replaced by freshwater due to the interface seaward migration as a result of the drop in the level of the Dead Sea. For that purpose, the dynamic equation of Glover approach has been integrated to accommodate that type of interface readjustment. The calculated amounts of freshwater which substituted salt Dead Sea water due to the migration of interface are 3.21 · 1011 m3, from a Dead Sea level of –392 m to τ411 m below sea level. The average porosity of coastal aquifers was calculated to range from 2.8 to 2.94%. Geoelectric sounding measurements showed that areas underlying the coastal aquifers formerly occupied by the Dead Sea water are gradually becoming flushed and occupied by freshwater. The latter is becoming salinized due to the residuals of Dead Sea water in the aquifer matrix, the present salinity of which is lower than that of the Dead Sea water. At the same time salt dissolution from the Lisan Marl formation is causing collapses along the shorelines in the form of sinkholes, tens of meters in diameter and depth.  相似文献   

16.
Low-permeability layer (LPL), formed by natural deposit or artificial reclamation and commonly found below the intertidal zone of coastal groundwater system, can retard the ingress of seawater and contaminants, and shorten the travel time of the land-sourced contaminant to the marine environment compared with a homogenous sandy coastal aquifer. However, there is limited understanding on how an intertidal LPL, a condition occurred in a coastal aquifer at Moreton Bay, Australia, influences the groundwater and contaminant transport across the shallow beach aquifer system. We characterized the aquifer hydrological parameters, monitored the in situ groundwater heads, and constructed a 2-D numerical model to analyses the cross-shore hydrological processes in this stratified system. The calibrated model suggests that in the lower aquifer, the inland-source fresh groundwater flowed horizontally towards the sea, upwelled along the freshwater–saltwater interface, and exited the aquifer at the shore below the LPL. Whereas in the upper aquifer, the tidally driven seawater circulation formed a barrier that prevented fresh groundwater from horizontal transport and discharge to the beach above the LPL, thereby directing its leakage to the lower aquifer. A contaminant represented by a conservative tracer was ‘released’ the upper aquifer in the model and results showed that the spreading extent of the contaminant plume, the maximum rate of contaminant discharge to the ocean, and its plume length decreased compared with a simulation case in a homogenous sandy aquifer. Sensitivity analysis was also conducted to investigate the characteristics of the LPL, including its continuity and hydraulic conductivity, which were found to vary along the beach at Moreton Bay. The result shows that with a lower hydraulic conductivity and continuous layer of LPL reduced the groundwater exchange and contaminant transport between upper and lower aquifer. The findings from the combined field and modelling investigations on the impact of an intertidal LPL on coastal aquifer systems highlight its significant implications to alter the groundwater and mass transport across the land–ocean interface.  相似文献   

17.
Electrical and electromagnetic methods are well suited for coastal aquifer studies because of the large contrast in resistivity between fresh water-bearing and salt water-bearing formations. Interpretation models for these aquifers typically contain four layers: a highly resistive unsaturated zone; a surficial fresh water aquifer of intermediate resistivity; an underlying conductive, salt water saturated aquifer; and resistive substratum. Additional layers may be added to allow for variations in lithology within the fresh water and salt water layers. Two methods are evaluated: direct current resistivity and time domain electromagnetic soundings. Use of each method alone produces nonunique solutions for resistivities and/or thicknesses of the different layers. We show that joint inversion of vertical electric and time domain electromagnetic soundings produces a more tightly constrained interpretation model at three test sites than is produced by inversion methods applied to each data set independently.  相似文献   

18.
This review focuses on investigations of groundwater flow and solute transport in karst aquifers through laboratory scale models (LSMs). In particular, LSMs have been used to generate new data under different hydraulic and contaminant transport conditions, testing of new approaches for site characterization, and providing new insights into flow and transport processes through complex karst aquifers. Due to the increasing need for LSMs to investigate a wide range of issues, associated with flow and solute migration karst aquifers this review attempts to classify, and introduce a framework for constructing a karst aquifer physical model that is more representative of field conditions. The LSMs are categorized into four groups: sand box, rock block, pipe/fracture network, and pipe-matrix coupling. These groups are compared and their advantages and disadvantages highlighted. The capabilities of such models have been extensively improved by new developments in experimental methods and measurement devices. Newer technologies such as 3D printing, computed tomography scanning, X-rays, nuclear magnetic resonance, novel geophysical techniques, and use of nanomaterials allow for greater flexibilities in conducting experiments. In order for LSMs to be representative of karst aquifers, a few requirements are introduced: (1) the ability to simulate heterogeneous distributions of karst hydraulic parameters, (2) establish Darcian and non-Darcian flow regimes and exchange between the matrix and conduits, (3) placement of adequate sampling points and intervals, and (4) achieving some degree of geometric, kinematic, and dynamic similitude to represent field conditions.  相似文献   

19.
Numerical investigation of road salt impact on an urban wellfield   总被引：2，自引：0，他引：2  
The impact of road salt on a wellfield in a complex glacial moraine aquifer system is studied by numerical simulation. The moraine underlies an extensive urban and industrial landscape, which draws its water supply from >20 wellfields, several of which are approaching or have exceeded the drinking water limit for chloride. The study investigates the mechanisms of road salt infiltration, storage, and transport in the subsurface and assesses the effectiveness of mitigation measures designed to reduce the impact. The three-dimensional transport model accounts for increases in salt loading, as well as growth of the urbanized area and road network over the past 50 years. The simulations, which focus on one impacted wellfield, show chloride plumes originating mainly at arterial roads and migrating through aquitard windows into the water supply aquifers. The results suggest that the aquifer system contains a large and heterogeneously distributed mass of chloride and that concentrations in the aquifer can be substantially higher than the concentrations in the well water. Future impact scenarios indicate that although the system responds rapidly to reductions in salt loading, the residual chloride mass may take decades to flush out, even if road salting were discontinued. The implications with respect to urban wellfields in typical snow-belt areas are discussed.  相似文献   

20.
Groundwater flow and chemical transport in subterranean estuaries are poorly understood despite their potentially important implications for chemical fluxes from aquifers to coastal waters. Here, a numerical study of the dynamics in a subterranean estuary subject to tidal forcing is presented. Simulations show that salt transport associated with tidally driven seawater recirculation leads to the formation of an upper saline plume in the intertidal region. Computed transit times and flow velocities indicate that this plume represents a more active zone for mixing and reaction than the dispersion zone of the lower, classical salt wedge. Proper conceptualisation of this surficial mixing zone extends our understanding of processes within the subterranean estuary. Numerical tracer simulations reveal that tidal forcing may reduce the threat of a land-derived contaminant discharging to the marine environment by modifying the subsurface transport pathway and local geochemical conditions. Mixing and stratification in the subterranean estuary are strongly affected by both inland and tidal forcing. Based on the estuarine analogy we present a systematic classification of subterranean estuaries.  相似文献   

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

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