首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
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.  相似文献   

2.
This article provides a partial answer to the question “What is the relation between excess hydraulic head and volume flux of water in a conduit within a porous matrix?”, focusing on the case that the forcing is steady. The conduit is modelled as a horizontal circular cylinder, imbedded within a porous matrix of rectangular cross section, having constant head prescribed on the sidewalls and being confined top and bottom. Laminar flow in the matrix is assumed to obey Darcy's law, while turbulent flow in the conduit is quantified using the Darcy–Weisbach equation. Analysis of the latter equation shows that the length scale of variations in the direction of the conduit is large compared with the scale of lateral and vertical variations. This permits separation of the full three-dimensional non-linear problem into a two-dimensional linear problem for head within the matrix and a one-dimensional non-linear problem for head within the conduit. Analytic solutions are obtained for the distribution of head in the matrix and in a conduit of either infinite or finite length. In both cases, the volume flux of water is proportional to the excess head to the 2/3 power, the conduit radius to the 5/3 power, the matrix permeability to the 1/3 power and gravity to the 1/3 power. The scale of variation of head along the conduit is proportional to the excess head to the ?1/3 power, the conduit radius to the 5/3 power, the matrix permeability to the ?2/3 power and gravity to the 1/3 power.  相似文献   

3.
Application of a Discrete-Continuum Model to Karst Aquifers in North China   总被引：7，自引：0，他引：7  
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.  相似文献   

4.
South China karst aquifer storm-scale hydrochemistry   总被引：7，自引：0，他引：7  
Liu Z  Groves C  Yuan D  Meiman J 《Ground water》2004,42(4):491-499
  相似文献   

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

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

7.
Karst aquifers represent dual flow systems consisting of a highly conductive conduit system embedded in a less permeable rock matrix. Hybrid models iteratively coupling both flow systems generally consume much time, especially because of the nonlinearity of turbulent conduit flow. To reduce calculation times compared to those of existing approaches, a new iterative equation solver for the conduit system is developed based on an approximated Newton–Raphson expression and a Gauß–Seidel or successive over‐relaxation scheme with a single iteration step at the innermost level. It is implemented and tested in the research code CAVE but should be easily adaptable to similar models such as the Conduit Flow Process for MODFLOW‐2005. It substantially reduces the computational effort as demonstrated by steady‐state benchmark scenarios as well as by transient karst genesis simulations. Water balance errors are found to be acceptable in most of the test cases. However, the performance and accuracy may deteriorate under unfavorable conditions such as sudden, strong changes of the flow field at some stages of the karst genesis simulations.  相似文献   

8.
Karst spring responses examined by process-based modeling   总被引：8，自引：0，他引：8  
Birk S  Liedl R  Sauter M 《Ground water》2006,44(6):832-836
Ground water in karst terrains is highly vulnerable to contamination due to the rapid transport of contaminants through the highly conductive conduit system. For contamination risk assessment purposes, information about hydraulic and geometric characteristics of the conduits and their hydraulic interaction with the fissured porous rock is an important prerequisite. The relationship between aquifer characteristics and short-term responses to recharge events of both spring discharge and physicochemical parameters of the discharged water was examined using a process-based flow and transport model. In the respective software, a pipe-network model, representing fast conduit flow, is coupled to MODFLOW, which simulates flow in the fissured porous rock. This hybrid flow model was extended to include modules simulating heat and reactive solute transport in conduits. The application of this modeling tool demonstrates that variations of physicochemical parameters, such as solute concentration and water temperature, depend to a large extent on the intensity and duration of recharge events and provide information about the structure and geometry of the conduit system as well as about the interaction between conduits and fissured porous rock. Moreover, the responses of solute concentration and temperature of spring discharge appear to reflect different processes, thus complementing each other in the aquifer characterization.  相似文献   

9.
Combination of a conduit flow law with the exceedance probability of stream discharge allows the estimation of the exceedance probability of water level in a simple conduit aquifer. The probability of water levels higher than available relief is interpreted as the probability of surface as opposed to underground runoff. Very high probability of surface runoff implies a fluvial environment, whereas very low probabilities define a mature karst or ‘holokarst’, Intermediate probabilities identify ‘fluviokarst’. Overflow probability depends on available relief, mean discharge, and especially conduit radius. Growth rate of the underground conduit depends on saturation deficit which thus controls the rate of evolution of the fluvial landscape, through fluviokarst to holokarst. However, variations of discharge and sedimentation through time can cause dramatic reversion of karst drainage into less mature states. Landscapes experiencing such periodic rejuvenation will have a confused morphology. A functional definition of landscape may be more objective and pertinent than arbitrary interpretation of form.  相似文献   

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

11.
Temporal age distributions are equivalent to probability density functions (PDFs) of transit time. The type and shape of a PDF provides important information related to ground-water mixing at the well or spring and the complex nature of flow networks in karst aquifers. Chlorofluorocarbon (CFC) concentrations measured for samples from 12 locations in the karstic Madison aquifer were used to evaluate the suitability of various PDF types for this aquifer. Parameters of PDFs could not be estimated within acceptable confidence intervals for any of the individual sites. Therefore, metrics derived from CFC-based apparent ages were used to evaluate results of PDF modeling in a more general approach. The ranges of these metrics were established as criteria against which families of PDFs could be evaluated for their applicability to different parts of the aquifer. Seven PDF types, including five unimodal and two bimodal models, were evaluated. Model results indicate that unimodal models may be applicable to areas close to conduits that have younger piston (i.e., apparent) ages and that bimodal models probably are applicable to areas farther from conduits that have older piston ages. The two components of a bimodal PDF are interpreted as representing conduit and diffuse flow, and transit times of as much as two decades may separate these PDF components. Areas near conduits may be dominated by conduit flow, whereas areas farther from conduits having bimodal distributions probably have good hydraulic connection to both diffuse and conduit flow.  相似文献   

12.
Ground water flow in karst terranes generally occurs in the solution channels of carbonate aquifers. A hydrogeologist may utilize borehole geophysical methods to identify these solution channels in aquifers. Two specific methods that are applicable in karst terrains are:
1. Natural gamma ray logging
2. Borehole caliper logging.
Gamma ray logging can detect the presence of inter-bedded strata in the main limestone unit, such as shale, which emit high levels of gamma radiation. Gamma ray logging can also detect clay deposits in solution channels that may act to restrict the flow of ground water. The areal extent of these rock strata or clay-filled solution channels can be determined when gamma ray logs are conducted at several borehole locations across the site of investigation.
Borehole caliper logging can be employed to determine the presences of solution channels within the aquifer when penetrated by a borehole. In addition, since shale layers and clay filling are less resistant than the surrounding limestone, the caliper log may detect both the presence and the thickness of shale or clay layers in the aquifer.
Gamma ray logs can be used in conjunction with caliper logs to provide data on the stratigraphic location and thickness of solution channels and clay and shale layers within a limestone aquifer. This information is valuable to the hydrogeologist performing investigations at sites located in limestone terranes because ground water flow preferentially occurs along solution channels.  相似文献   

13.
An analysis method for slug tests performed in a partially penetrating well within a vertical cutoff wall is presented. A steady‐state shape factor for evaluating hydraulic conductivity of the material within the wall was derived by applying the method of images to the previously developed analytical solution of Zlotnik et al. (2010) for an infinite aquifer. Two distinct boundary conditions were considered: constant‐head boundary for the case of direct contact between the wall and the aquifer, and no‐flux boundary representing an impermeable filter cake on the sides of the wall. The constant‐head and no‐flux boundary conditions yield significantly higher and lower shape factors, respectively, than those for the infinite aquifer. Consequently the conventional line‐fitting method for slug test analysis would yield an inaccurate estimate of the hydraulic conductivity of a vertical cutoff wall.  相似文献   

14.
In this paper, we present a conceptual‐numerical model that can be deduced from a calibrated finite difference groundwater‐flow model, which provides a parsimonious approach to simulate and analyze hydraulic heads and surface water body–aquifer interaction for linear aquifers (linear response of head to stresses). The solution of linear groundwater‐flow problems using eigenvalue techniques can be formulated with a simple explicit state equation whose structure shows that the surface water body–aquifer interaction phenomenon can be approached as the drainage of a number of independent linear reservoirs. The hydraulic head field could be also approached by the summation of the head fields, estimated for those reservoirs, defined over the same domain set by the aquifer limits, where the hydraulic head field in each reservoir is proportional to a specific surface (an eigenfunction of an eigenproblem, or an eigenvector in discrete cases). All the parameters and initial conditions of each linear reservoir can be mathematically defined in a univocal way from the calibrated finite difference model, preserving its characteristics (geometry, boundary conditions, hydrodynamic parameters (heterogeneity), and spatial distribution of the stresses). We also demonstrated that, in practical cases, an accurate solution can be obtained with a reduced number of linear reservoirs. The reduced computational cost of these solutions can help to integrate the groundwater component within conjunctive use management models. Conceptual approximation also facilitates understanding of the physical phenomenon and analysis of the factors that influence it. A simple synthetic aquifer has been employed to show how the conceptual model can be built for different spatial discretizations, the parameters required, and their influence on the simulation of hydraulic head fields and stream–aquifer flow exchange variables. A real‐world case was also solved to test the accuracy of the proposed approaches, by comparing its solution with that obtained using finite‐difference MODFLOW code. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

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

16.
As a result of rock dissolution processes, karst aquifers exhibit highly conductive features such as caves and conduits. Within these structures, groundwater flow can become turbulent and therefore be described by nonlinear gradient functions. Some numerical groundwater flow models explicitly account for pipe hydraulics by coupling the continuum model with a pipe network that represents the conduit system. In contrast, the Conduit Flow Process Mode 2 (CFPM2) for MODFLOW-2005 approximates turbulent flow by reducing the hydraulic conductivity within the existing linear head gradient of the MODFLOW continuum model. This approach reduces the practical as well as numerical efforts for simulating turbulence. The original formulation was for large pore aquifers where the onset of turbulence is at low Reynolds numbers (1 to 100) and not for conduits or pipes. In addition, the existing code requires multiple time steps for convergence due to iterative adjustment of the hydraulic conductivity. Modifications to the existing CFPM2 were made by implementing a generalized power function with a user-defined exponent. This allows for matching turbulence in porous media or pipes and eliminates the time steps required for iterative adjustment of hydraulic conductivity. The modified CFPM2 successfully replicated simple benchmark test problems.  相似文献   

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

18.
Complex aquifer systems are often modeled with quasi-three-dimensional models, which consider two-dimensional horizontal flow in the aquifers and one-dimensional vertical flow through aquitards. When the aquifer system consists of a phreatic aquifer and one or more semiconfined aquifers connected by aquitards, the discrete model consists of a nonlinear system of algebraic equations, because the transmissivity of the phreatic aquifer depends on the phreatic head. If the water extraction is very high, the phreatic aquifer can be depleted and the equations of the model must be modified accordingly. There are not simple and general criteria to state if the phreatic aquifer is depleted before solving the system of equations. Therefore, the iterative procedures (e.g., relaxation methods), used to find the solution to the forward problem, must handle these particular conditions and can suffer several problems of convergence. These problems can be caused by the choice of the initial head values or of the relaxation coefficient of the iterative algorithms; however, they can also be caused by the nonexistence or nonuniqueness of the solution to the system of nonlinear equations. The study of existence and uniqueness of the general problem is very difficult and, therefore, we consider a simplified problem, for which the discrete model can be handled analytically. The results of the numerical experiments show that the solution to the forward problem can be nonunique. Only for some cases it is possible to invoke physical arguments to eliminate tentative solutions.  相似文献   

19.
In variably confined carbonate platforms, impermeable confining units collect rainfall over large areas and deliver runoff to rivers or conduits in unconfined portions of platforms. Runoff can increase river stage or conduit heads in unconfined portions of platforms faster than local infiltration of rainfall can increase groundwater heads, causing hydraulic gradients between rivers, conduits and the aquifer to reverse. Gradient reversals cause flood waters to flow from rivers and conduits into the aquifer where they can dissolve limestone. Previous work on impacts of gradient reversals on dissolution has primarily emphasized individual caves and little research has been conducted at basin scales. To address this gap in knowledge, we used legacy data to assess how a gradient of aquifer confinement across the Suwannee River Basin, north‐central Florida affected locations, magnitudes and processes of dissolution during 2005–2007, a period with extreme ranges of discharge. During intense rain events, runoff from the confining unit increased river stage above groundwater heads in unconfined portions of the platform, hydraulically damming inputs of groundwater along a 200 km reach of river. Hydraulic damming allowed allogenic runoff with SICAL < ?4 to fill the entire river channel and flow into the aquifer via reversing springs. Storage of runoff in the aquifer decreased peak river discharges downstream and contributed to dissolution within the aquifer. Temporary storage of allogenic runoff in karst aquifers represents hyporheic exchange at a scale that is larger than found in streams flowing over non‐karst aquifers because conduits in karst aquifers extend the area available for exchange beyond river beds deep into aquifers. Post‐depositional porosity in variably confined carbonate platforms should thus be enhanced along rivers that originate on confining units. This distribution should be considered in models of porosity distribution used to manage water and hydrocarbon resources in carbonate rocks. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

20.
A confined aquifer may become unconfined near the pumping wells when the water level falls below the confining unit in the case where the pumping rate is great and the excess hydraulic head over the top of the aquifer is small. Girinskii's potential function is applied to analyze the steady ground water flow induced by pumping wells with a constant-head boundary in a mixed confined-unconfined aquifer. The solution of the single-well problem is derived, and the critical radial distance at which the flow changes from confined to unconfined condition is obtained. Using image wells and the superposition method, an analytic solution is presented to study steady ground water flow induced by a group of pumping wells in an aquifer bounded by a river with constant head. A dimensionless function is introduced to determine whether a water table condition exists or not near the pumping wells. An example with three pumping wells is used to demonstrate the patterns of potentiometric surface and development of water table around the wells.  相似文献   

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

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