共查询到20条相似文献，搜索用时 31 毫秒
1.
A tracer test in a carbonate aquifer is analyzed using the method of moments and two analytical advection-dispersion models (ADMs) as well as a numerical model. The numerical model is a coupled continuum-pipe flow and transport model that accounts for two different flow components in karstified carbonate aquifers, i.e., rapid and often turbulent conduit flow and Darcian flow in the fissured porous rock. All techniques employed provide reasonable fits to the tracer breakthrough curve (TBC) measured at a spring. The resulting parameter estimates are compared to investigate how each conceptual model of flow and transport processes that forms the basis of the analyses affects the interpretation of the tracer test. Numerical modeling results suggest that the method of moments and the analytical ADMs tend to overestimate the conduit volume because part of the water discharged at the spring is wrongly attributed to the conduit system if flow in the fissured porous rock is ignored. In addition, numerical modeling suggests that mixing of the two flow components accounts for part of the dispersion apparent in the measured TBC, while the remaining part can be attributed to Taylor dispersion. These processes, however, cannot reasonably explain the tail of the TBC. Instead, retention in immobile-fluid regions as included in a nonequilibrium ADM provides a possible explanation. 相似文献
2.
Locating the zone of saline intrusion in a coastal karst aquifer using springflow data 总被引：2，自引：0，他引：2
Coastal fresh water aquifers are an increasingly desirable resource. In a karstic aquifer, sea water intrusion occurs as a salt water wedge, like in porous media. However, preferential flow conduits may alter the spatial and temporal distribution of the salt water. This is typically the case when the outlet of the aquifer is a brackish spring. This paper shows that salinity and flow rate variations at a spring, where salinity is inversely proportional to discharge, can help to understand the hydrodynamic functioning of the aquifer and to locate the fresh water-sea water mixing zone deep inside the aquifer. The volume of water-filled conduit between the sea water intrusion zone and the spring outlet is calculated by the integral over time of the flow rate during the time lag between the flow rate increase and the salinity decrease as measured at the spring. In the example of the spring at Almyros of Heraklio (Crete, Greece), this time lag is variable, depending on the discharge, but the volume of water-filled conduit appears to be constant, which shows that the processes of salt water intrusion and mixing in the conduit are constant throughout the year. The distance between the spring and the zone where sea water enters the conduit is estimated and provides an indication of the position where only fresh water is present in the conduit. 相似文献
3.
Studying the Flow Dynamics of a Karst Aquifer System with an Equivalent Porous Medium Model 总被引：2，自引：0，他引：2
The modeling of groundwater flow in karst aquifers is a challenge due to the extreme heterogeneity of its hydraulic parameters and the duality in their discharge behavior, that is, rapid response of highly conductive karst conduits and delayed drainage of the low‐permeability fractured matrix after recharge events. There are a number of different modeling approaches for the simulation of the karst groundwater dynamics, applicable to different aquifer as well as modeling problem types, ranging from continuum models to double continuum models to discrete and hybrid models. This study presents the application of an equivalent porous model approach (EPM, single continuum model) to construct a steady‐state numerical flow model for an important karst aquifer, that is, the Western Mountain Aquifer Basin (WMAB), shared by Israel and the West‐Bank, using MODFLOW2000. The WMAB was used as a catchment since it is a well‐constrained catchment with well‐defined recharge and discharge components and therefore allows a control on the modeling approach, a very rare opportunity for karst aquifer modeling. The model demonstrates the applicability of equivalent porous medium models for the simulation of karst systems, despite their large contrast in hydraulic conductivities. As long as the simulated saturated volume is large enough to average out the local influence of karst conduits and as long as transport velocities are not an issue, EPM models excellently simulate the observed head distribution. The model serves as a starting basis that will be used as a reference for developing a long‐term dynamic model for the WMAB, starting from the pre‐development period (i.e., 1940s) up to date. 相似文献
4.
by Guangquan Li 《Ground water》2009,47(5):714-722
The permeable conduit wall in a karst aquifer allows for water and solute to be exchanged between conduits and the limestone matrix. Contaminant sequestered in the limestone matrix is flushed into conduits following flood events. The contaminant released from the permeable wall will then mix with conduit water and will be transported downgradient in the conduit. A one-dimensional advection-dispersion equation is presented to describe this mixing-transport incorporating water flow and solute flux through the conduit wall. An analytical solution ignoring conduit dispersion is derived using the method of characteristics. Scale analysis is performed to provide a general guideline to estimate when conduit dispersion can be neglected. The solution also can be used to compute the distribution of solute in the matrix before flushing. 相似文献
5.
Flow and transport simulation in karst aquifers remains a significant challenge for the ground water modeling community. Darcy's law–based models cannot simulate the inertial flows characteristic of many karst aquifers. Eddies in these flows can strongly affect solute transport. The simple two-region conduit/matrix paradigm is inadequate for many purposes because it considers only a capacitance rather than a physical domain. Relatively new lattice Boltzmann methods (LBMs) are capable of solving inertial flows and associated solute transport in geometrically complex domains involving karst conduits and heterogeneous matrix rock. LBMs for flow and transport in heterogeneous porous media, which are needed to make the models applicable to large-scale problems, are still under development. Here we explore aspects of these future LBMs, present simple examples illustrating some of the processes that can be simulated, and compare the results with available analytical solutions. Simulations are contrived to mimic simple capacitance-based two-region models involving conduit (mobile) and matrix (immobile) regions and are compared against the analytical solution. There is a high correlation between LBM simulations and the analytical solution for two different mobile region fractions. In more realistic conduit/matrix simulation, the breakthrough curve showed classic features and the two-region model fit slightly better than the advection-dispersion equation (ADE). An LBM-based anisotropic dispersion solver is applied to simulate breakthrough curves from a heterogeneous porous medium, which fit the ADE solution. Finally, breakthrough from a karst-like system consisting of a conduit with inertial regime flow in a heterogeneous aquifer is compared with the advection-dispersion and two-region analytical solutions. 相似文献
6.
J.A. Barker 《Advances in water resources》1982,5(2):98-104
The main processes affecting the migration of a solute in a fissured aquifer will be advection and dispersion in the fissures, diffusion into the porous matrix; and adsorption. This paper considers solute transport in an idealized fissured aquifer consisting of slabs of saturated rock-matrix separated by equally spaced, planar fissures. The solution of the transport equations is developed as far as Laplace transforms of the solute concentrations in the fissure and matrix water. Numerical inversion of the transforms is used to investigate characteristic behaviour of the model for a number of special cases. 相似文献
7.
C. C. Smart 《水文研究》1988,2(1):31-41
Discharge exceedance probabilities are calculated for a simple model karst aquifer composed of a few multilevel conduits with recharge from a single sinking stream with an exponential flow exceedance distribution. It is assumed that outflow instantaneously matches inflow, so that the conduit volume is constant but discharge is governed by the head in a volumeless shaft at the top end of the system. It is shown that small single conduit aquifers will frequently overflow at the surface during floods and the exceedance probability of flow through the aquifer and over the surface can be defined as a function of the inflow distribution and the form of the aquifer. Systems with multiple conduits will overflow less frequently, but each conduit will exhibit a flow distribution characteristic of its form and position in the vertical hierarchy. Comparison of these findings with actual flow data from a conduit aquifer shows that the approach is valid, although imprecise. The model is unlikely to be applied directly, as it requires unusually detailed data. However, it provides for the first time an indication of typical flow durations for surface overflows and individual conduits in a karst aquifer. Contrasts in flow duration will have a profound influence on solutional and sedimentary processes in the karst system. 相似文献
8.
A solution conduit has a permeable wall allowing for water exchange and solute transfer between the conduit and its surrounding aquifer matrix. In this paper, we use Laplace Transform to solve a one‐dimensional equation constructed using the Euler approach to describe advective transport of solute in a conduit, a production‐value problem. Both nonuniform cross‐section of the conduit and nonuniform seepage at the conduit wall are considered in the solution. Physical analysis using the Lagrangian approach and a lumping method is performed to verify the solution. Two‐way transfer between conduit water and matrix water is also investigated by using the solution for the production‐value problem as a first‐order approximation. The approximate solution agrees well with the exact solution if dimensionless travel time in the conduit is an order of magnitude smaller than unity. Our analytical solution is based on the assumption that the spatial and/or temporal heterogeneity in the wall solute flux is the dominant factor in the spreading of spring‐breakthrough curves, and conduit dispersion is only a secondary mechanism. Such an approach can lead to the better understanding of water exchange and solute transfer between conduits and aquifer matrix. Highlights:
- 相似文献
9.
A structure model was used to analyse solute‐transport parameter estimates based on tracer breakthrough curves. In the model system, groundwater flow is envisioned to be organised in a complex conduit network providing a variety of short circuits with relative small carrying capacities along different erosion levels. The discharge through the fully filled conduits is limited owing to void geometries and turbulent flow; thus, a hierarchic overflow system evolves where conduits are (re‐)activated or dried up depending on the flow condition. Exemplified on the Lurbach–Tanneben karst aquifer, the applicability of the model approach was tested. Information derived from multi‐tracer experiments performed at different volumetric flow rates enabled to develop a structural model of the karst network, under constraint of the geomorphological and hydrological evolution of the site. Depending on the flow rate, groundwater is divided into up to eight flow paths. The spatial hierarchy of flow paths controls the sequence of flow path activation. Conduits of the topmost level are strongly influenced by reversible alteration processes. Sedimentation or blocking causes an overflow of water to the next higher conduit. Flow path specific dissolutional denudation rates were estimated using the temporal development of the partial discharge. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
10.
Christoph Butscher Adrian Auckenthaler Stefan Scheidler Peter Huggenberger 《Ground water》2011,49(1):66-76
Rapid changes in spring water quality in karst areas due to rapid recharge of bacterially contaminated water are a major concern for drinking water suppliers and users. The main objective of this study was to use field experiments with fecal indicators to verify the vulnerability of a karst spring to pathogens, as determined by using a numerical modeling approach. The groundwater modeling was based on linear storage models that can be used to simulate karst water flow. The vulnerability of the karst groundwater is estimated using such models to calculate criteria that influence the likelihood of spring water being affected by microbial contamination. Specifically, the temporal variation in the vulnerability, depending on rainfall events and overall recharge conditions, can be assessed and quantified using the dynamic vulnerability index (DVI). DVI corresponds to the ratio of conduit to diffuse flow contributions to spring discharge. To evaluate model performance with respect to predicted vulnerability, samples from a spring were analyzed for Escherichia coli, enterococci, Clostridium perfringens, and heterotrophic plate count bacteria during and after several rainfall events. DVI was shown to be an indication of the risk of fecal contamination of spring water with sufficient accuracy to be used in drinking water management. We conclude that numerical models are a useful tool for evaluating the vulnerability of karst systems to pathogens under varying recharge conditions 相似文献
11.
In this study, two conceptual models, the classic reservoir (CR) model and exchange reservoirs model embedded by dual porosity approach (DPR) are developed for simulation of karst aquifer functioning drained by multiple outlets. The performances of two developed models are demonstrated at a less developed karstic aquifer with three spring outlets located in Zagros Mountain in the south‐west of Iran using 22‐years of daily data. During the surface recharge, a production function based on water mass balance is implemented for computing the time series of surface recharge to the karst formations. The efficiency of both models has been assessed for simulation of daily spring discharge during the recession and also surface recharge periods. Results indicate that both CR and DPR models are capable of simulating the ordinates of spring hydrographs which drainage less developed karstic aquifer. However, the goodness of fit criteria indicates outperformance of DPR model for simulation of total hydrograph ordinates. In addition, the DPR model is capable of quantifying hydraulic properties of two hydrologically connected overlapping continua conduits network and fissure matrix which lays important foundations for the mining operation and water resource management whereas homogeneous model representations of the karstic subsurface (e.g., the CR) do not work accurately in the karstic environment. 相似文献
12.
Catchment‐scale conceptual modelling of water and solute transport in the dual flow system of the karst critical zone 
下载免费PDF全文
下载免费PDF全文 Hydrological and hydrochemical processes in the critical zone of karst environments are controlled by the fracture‐conduit network. Modelling hydrological and hydrochemical dynamics in such heterogeneous hydrogeological settings remains a research challenge. In this study, water and solute transport in the dual flow system of the karst critical zone were investigated in a 73.5‐km2 catchment in southwest China. We developed a dual reservoir conceptual run‐off model combined with an autoregressive and moving average model with algorithms to assess dissolution rates in the “fast flow” and “slow flow” systems. This model was applied to 3 catchments with typical karst critical zone architectures, to show how flow exchange between fracture and conduit networks changes in relation to catchment storage dynamics. The flux of bidirectional water and solute exchange between the fissure and conduit system increases from the headwaters to the outfall due to the large area of the developed conduits and low hydraulic gradient in the lower catchment. Rainfall amounts have a significant influence on partitioning the relative proportions of flow and solutes derived from different sources reaching the underground outlet. The effect of rainfall on catchment function is modulated by the structure of the karst critical zone (e.g., epikarst and sinkholes). Thin epikarst and well‐developed sinkholes in the headwaters divert more surface water (younger water) into the underground channel network, leading to a higher fraction of rainfall recharge into the fast flow system and total outflow. Also, the contribution of carbonate weathering to mass export is also higher in the headwaters due to the infiltration of younger water with low solute concentrations through sinkholes. 相似文献
13.
This paper presents a vertically averaged model for studying water and solute exchanges between a large river and its adjacent alluvial aquifer. The hydraulic model couples horizontal 2D Saint Venant equations for river flow and a 2D Dupuit equation for aquifer flow. The dynamic coupling between river and aquifer is provided by continuity of fluxes and water level elevation between the two domains. Equations are solved simultaneously by linking the two hydrological system matrices in a single global matrix in order to ensure the continuity conditions between river and aquifer and to accurately model two‐way coupling between these two domains. The model is applied to a large reach (about 36 km2) of the Garonne River (south‐western France) and its floodplain, including an instrumented site in a meander. Simulated hydraulic heads are compared with experimental measurements on the Garonne River and aquifer in the floodplain. Model verification includes comparisons for one point sampling date (27 piezometers, 30 March 2000) and for hydraulic heads variations measured continuously over 5 months (5 piezometers, 1 January to 1 June 2000). The model accurately reproduces the strong hydraulic connections between the Garonne River and its aquifer, which are confirmed by the simultaneous variation of the water level in the river and in piezometers located near the river bank. The simulations also confirmed that the model is able to reproduce groundwater flow dynamics during flood events. Given these results, the hydraulic model was coupled with a solute‐transport component, based on advection‐dispersion equations, to investigate the theoretical dynamics of a conservative tracer over 5 years throughout the 36 km2 reach studied. Meanders were shown to favour exchanges between river and aquifer, and although the tracer was diluted in the river, the contamination moved downstream from the injection plots and affected both river banks. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
14.
We infer information about the nature of groundwater flow within a karst aquifer from the physical and chemical response of a spring to storm events. The spring discharges from the Maynardville Limestone in Bear Creek Valley, Tennessee. Initially, spring discharge peaks approximately 1–2 h from the midpoint of summer storms. The initial peak is likely due to surface loading, which pressurizes the aquifer and results in water moving out of storage. All of the storms monitored exhibited recessions that follow a master recession curve very closely, indicating that storm response is fairly consistent and repeatable, independent of the time between storms and the configuration of the rain event itself. Electrical conductivity initially increases for 0.5–2.9 days (longer for smaller storms), the result of moving older water out of storage. This is followed by a 2.1–2.5 day decrease in conductivity, resulting from an increasing portion of low conductivity recharge water entering the spring. Stable carbon isotope data and the calcite saturation index of the spring water also support this conceptual model. Spring flow is likely controlled by displaced water from the aquifer rather than by direct recharge through the soil zone. 相似文献
15.
Flow from some springs in former glacial lakebeds of the Upper Midwest is extremely steady throughout the year and does not increase significantly after precipitation events or seasonal recharge. Analytical and simplified numerical models of spring systems were used to determine whether preferential ground water flow through high-permeability features in shallow sandstone aquifers could produce typical values of spring discharge and the unusually steady rates of spring flow. The analytical model is based on a one-dimensional solution for periodic ground water flow. Solutions to this model suggest that it is unlikely that a periodic forcing due to seasonal variations in areal recharge would propagate to springs in a setting where high-permeability features exist. The analytical model shows that the effective length of the aquifer, or the length of flowpaths to a spring, and the total transmissivity of the aquifer have the greatest potential to impact the nature of spring flow in this setting. The numerical models show that high-permeability features can influence the magnitude of spring flow and the results demonstrate that the lengths of ground water flowpaths increase when high-permeability features are explicitly modeled, thus decreasing the likelihood for temporal variations in spring flow. 相似文献
16.
A generalized discrete-continuum model is developed to simulate ground water flow in the karst aquifers of North China. The model is a hybrid numerical flow model, which takes into account both quick conduit flow and diffusive fissure flow. The conduit flow is represented by a discrete network model, and the fissure flow is modeled by a continuum approach. The developed model strongly emphasizes the function of the conduits in the flow fields. They control the general drainage pattern, as demonstrated in the simulation of a complex karst aquifer in North China. The model reproduces reasonably well the flow field in response to an unanticipated discharge of ground water from the karst aquifer into an underground mine based on the aquifer parameters that are manually calibrated from a multiple-well pumping test. Sensitivity of the model to the aquifer parameters was evaluated in the context of the case study. 相似文献
17.
Richard W. Healy Sosina S. Haile David L. Parkhurst Scott R. Charlton 《Ground water》2018,56(5):810-815
Variably saturated groundwater flow, heat transport, and solute transport are important processes in environmental phenomena, such as the natural evolution of water chemistry of aquifers and streams, the storage of radioactive waste in a geologic repository, the contamination of water resources from acid‐rock drainage, and the geologic sequestration of carbon dioxide. Up to now, our ability to simulate these processes simultaneously with fully coupled reactive transport models has been limited to complex and often difficult‐to‐use models. To address the need for a simple and easy‐to‐use model, the VS2DRTI software package has been developed for simulating water flow, heat transport, and reactive solute transport through variably saturated porous media. The underlying numerical model, VS2DRT, was created by coupling the flow and transport capabilities of the VS2DT and VS2DH models with the equilibrium and kinetic reaction capabilities of PhreeqcRM. Flow capabilities include two‐dimensional, constant‐density, variably saturated flow; transport capabilities include both heat and multicomponent solute transport; and the reaction capabilities are a complete implementation of geochemical reactions of PHREEQC. The graphical user interface includes a preprocessor for building simulations and a postprocessor for visual display of simulation results. To demonstrate the simulation of multiple processes, the model is applied to a hypothetical example of injection of heated waste water to an aquifer with temperature‐dependent cation exchange. VS2DRTI is freely available public domain software. 相似文献
18.
Sergei A. Fomin Vladimir A. ChugunovToshiyuki Hashida 《Advances in water resources》2011,34(2):205-214
The paper provides an introduction to fundamental concepts of mathematical modeling of mass transport in fractured porous heterogeneous rocks. Keeping aside many important factors that can affect mass transport in subsurface, our main concern is the multi-scale character of the rock formation, which is constituted by porous domains dissected by the network of fractures. Taking into account the well-documented fact that porous rocks can be considered as a fractal medium and assuming that sizes of pores vary significantly (i.e. have different characteristic scales), the fractional-order differential equations that model the anomalous diffusive mass transport in such type of domains are derived and justified analytically. Analytical solutions of some particular problems of anomalous diffusion in the fractal media of various geometries are obtained. Extending this approach to more complex situation when diffusion is accompanied by advection, solute transport in a fractured porous medium is modeled by the advection-dispersion equation with fractional time derivative. In the case of confined fractured porous aquifer, accounting for anomalous non-Fickian diffusion in the surrounding rock mass, the adopted approach leads to introduction of an additional fractional time derivative in the equation for solute transport. The closed-form solutions for concentrations in the aquifer and surrounding rocks are obtained for the arbitrary time-dependent source of contamination located in the inlet of the aquifer. Based on these solutions, different regimes of contamination of the aquifers with different physical properties can be readily modeled and analyzed. 相似文献
19.
Analytical solutions for the water table and lateral discharge in a heterogeneous unconfined aquifer with time-dependent source and fluctuating river stage were derived and compared with those in an equivalent homogeneous aquifer. The heterogeneous aquifer considered consists of a number of sections of different hydraulic conductivity values. The source term and river stage were assumed to be time-dependent but spatially uniform. The solutions derived is useful in studying various groundwater flow problems in a horizontally heterogeneous aquifer since the spatially piecewise-constant hydraulic conductivity and temporally piecewise-constant recharge and lateral discharge can be used to quantify variations in these processes commonly observed in reality. Applying the solutions derived to an aquifer of three sections of different hydraulic conductivity values shown that (1) the aquifer heterogeneity significantly increases the spatial variation of the water table and thus its gradient but it has little effect on lateral discharge in the case of temporally and spatially uniform recharge, (2) the time-dependent but spatially uniform recharge increases the temporal variation of groundwater table over the entire aquifer but its effect on lateral discharge is limited in the zone near the river, and (3) the effect of river stage fluctuation on the water table and lateral discharge is limited in the zone near the river and the effect of the heterogeneity is to increase lateral discharge to or recharge from the river. 相似文献
20.
《Advances in water resources》2002,25(7):733-746
Backward location and travel time probabilities, which provide information about the former location of contamination in an aquifer, can be used to identify unknown contamination sources. Backward location probability describes the possible upgradient positions of contamination at a known time in the past, and backward travel time probability describes the time required for contamination to travel from a known upgradient location to an observation point. These probabilities are related to adjoint states of resident concentration, and their governing equation is the adjoint of a forward contaminant transport model. Using adjoint theory to obtain the appropriate governing equation, we extend the backward probability model for conservative solutes to more general non-uniform and transient flow fields. In particular, we address three important extensions, spatially-varying porosity, transient flow and temporally-varying porosity, and internal distributed sources and sinks of solute and water. For the first time we learn that forward and backward location and travel time probabilities are not necessarily equivalent to adjoint states, but are related to them. The extensions are illustrated using a vertically-integrated groundwater model, creating transient flow by a step change in pumping and using areal recharge as an internal distributed source. Both the movement and spread of probabilities are affected. With internal sources of water, there are two interpretations of backward probability, depending on whether or not the source of water is also a source of solute. The results demonstrate how the backward probability model can be applied to other, perhaps more important, non-uniform and transient flow conditions, with time- and space-varying water storage, such as time-varying pumping or unsaturated (or saturated–unsaturated) flow and transport with spatially- and temporally-varying moisture content. 相似文献