首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 297 毫秒
1.
The interpretation of groundwater ages in coastal aquifers requires an improved understanding of relationships between age distributions and the processes accompanying dispersive, density-dependent flow and transport. This study uses numerical modeling to examine the influence of mixing and a selection of other hydrogeological factors on steady-state age distributions in coastal aquifers. Three methods of age estimation are compared: the piston flow age, the direct age, and the tracer-based age. These are applied to various forms of the Henry problem, as well as to three variants of a larger, hypothetical coastal aquifer. Circulation of water within the seawater wedge results in markedly higher ages in the transition zone than in the underlying saltwater or overlying freshwater. Piston flow ages show a sharp increase where the freshwater and saltwater systems meet, whereas direct- and tracer-based simulations result in a smoother age distribution, as expected. Greater degrees of mixing result in larger differences between piston flow and direct or tracer-based ages, and bring about lower ages in the saltwater wedge. Tracer-based ages are preferred over direct- and piston flow ages for comparison with field data, especially in cases with wide transition zones. Despite the relatively simple conditions used for the simulations, complex age distributions with depth were obtained. Hence, the assessment of ages in field cases will be difficult, particularly where accurate ages in the transition zone are sought.  相似文献   

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

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

4.
A benchmark test for flow in karstic aquifers is presented in the form of an exact solution of the harmonic variations of water flux and head within a karst conduit that is imbedded within a three-dimensional porous matrix having a free surface. The variations are driven by a prescribed variation of head applied at one end of the conduit. The benchmark consists of expressions for the spring discharge as a function of time and the conduit head and flux as functions of distance along the conduit and time. These expressions contain three dimensionless parameters, permitting development of a wide range of specific benchmark tests. The expressions are particularly simple in the case of an infinitely deep aquifer. This limiting solution should provide the most severe test for two-dimensional models of karst aquifer flow. Another limiting case of interest is that in which the conduit diameter is equal to the water depth. This limiting solution should provide the easiest test for two-dimensional models.  相似文献   

5.
In many fields of study, and certainly in hydrogeology, uncertainty propagation is a recurring subject. Usually, parametrized probability density functions (PDFs) are used to represent data uncertainty, which limits their use to particular distributions. Often, this problem is solved by Monte Carlo simulation, with the disadvantage that one needs a large number of calculations to achieve reliable results. In this paper, a method is proposed based on a piecewise linear approximation of PDFs. The uncertainty propagation with these discretized PDFs is distribution independent. The method is applied to the upscaling of transmissivity data, and carried out in two steps: the vertical upscaling of conductivity values from borehole data to aquifer scale, and the spatial interpolation of the transmissivities. The results of this first step are complete PDFs of the transmissivities at borehole locations reflecting the uncertainties of the conductivities and the layer thicknesses. The second step results in a spatially distributed transmissivity field with a complete PDF at every grid cell. We argue that the proposed method is applicable to a wide range of uncertainty propagation problems.  相似文献   

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

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

8.
In the last few years, the scientific community has developed several hydrological models aimed at the simulation of hydrological processes acting at the basin scale. In this context, the portion of peak runoff contributing areas represents a critical variable for a correct estimate of surface runoff. Such areas are strongly influenced by the saturated portion of a river basin (influenced by antecedent conditions) but may also evolve during a specific rainfall event. In the recent years, we have developed 2 theoretically derived probability distributions that attempt to interpret these 2 processes adopting daily runoff and flood‐peak time series. The probability density functions (PDFs) obtained by these 2 schematisations were compared for humid river basins in southern Italy. Results highlighted that the PDFs of the peak runoff contributing areas can be interpreted by a gamma distribution and that the PDF of the relative saturated area provides a good interpretation of such process that can be used for flood prediction.  相似文献   

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

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

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

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

14.
We apply geospatial analysis to borehole imagery in an effort to develop new techniques to evaluate the spatial distribution and internal structure of karst conduits. Remote sensing software is used to classify a high resolution, digital borehole image of limestone bedrock from the Biscayne aquifer (South Florida, USA) into a binary image divided into cells of rock matrix and pores. Within a GIS, 2D porosity is calculated for a series of rectangular sampling windows placed over the binary image and then plotted as a function of depth. Potential conduits that intersect the borehole are identified as peaks of high porosity. A second GIS technique identifies a conduit as a continuous object that spans the entire borehole width. According to these criteria, geospatial analysis reveals ∼10 discrete conduits along the ∼15 m borehole image. Continuous sampling of the geologic medium intersected by the borehole provides insight into the internal structure of karst aquifers and the evolution of karst features. Most importantly, this pilot study demonstrates that GIS-based techniques are capable of quantifying the depths, dimensions, shapes, apertures and connectivity of potential conduits, physical attributes that impact flow in karst aquifers.  相似文献   

15.
The evolution of bubble size distributions in volcanic eruptions   总被引:1,自引:0,他引:1  
We review observations of bubble size distributions (BSDs) generated during explosive volcanic eruptions and laboratory explosions, as inferred from vesicle size distributions found in the end products. Unimodal, polymodal, exponential and power law BSDs are common, even in the absence of coalescence, and both power law and exponential distributions have been generated in the same eruption. To date theoretical models have proposed incompatible mechanisms for producing the various distributions. We here present a unifying mechanism. Data from our laboratory analogue experiments suggest that power law distributions are associated with highly non-equilibrium degassing. A numerical model is developed in which bubbles nucleate repeatedly and grow in the spaces between those of previous generations, where, in a non-equilibrium degassing scenario, the volatile concentration remains high. This process causes the BSD to evolve from unimodal, through exponential, into a power law. The exponent of the power law is a measure of the number of nucleation events, or the duration of the nucleation period compared with the timescale of bubble growth. The mathematical inevitability of the evolution from unimodal (Poissonian) to power law is discussed. The findings may resolve the apparent contradiction between the equilibrium degassing conduit flow models and the non-equilibrium degassing conditions derived from bubble growth models of explosive volcanic eruptions. The process of ongoing nucleation is the mechanism whereby the volcanic system maintains near-equilibrium in the case of rapid depressurisation and slow volatile diffusion.  相似文献   

16.
Hydrological and biogeochemical processes in karst environments are strongly controlled by heterogeneous fracture-conduit networks. Quantifying the spatio-temporal variability of water transit time and young water fractions in such heterogeneous hydrogeological systems is fundamental to linking discharge and water quality dynamics in the karst critical zone. We used a tracer-aided conceptual hydrological model to track the fate of each hour of rain input individually. Using this approach, the variability of transit time distributions and young water fraction were estimated in the main landscape units in a karst catchment of Chenqi in Guizhou Province, Southwest China. The model predicted that the mean young water (i.e., <~2 months old) fraction of ground conduit flow is 0.31. Marked seasonal variabilities in water storage and hydrological connectivity between the conduit network and fractured matrix, as well as between hillslopes and topographic depression, drive the dynamics of young water fraction and travel time distributions in each landscape unit. Especially, the strong hydrological connectivity between the land surface and underground conduits caused by the direct infiltration through large fractures and sinkholes, leads the drastic increasement in young water fraction of runoff after heavy rain. Even though the contribution of young water to runoff is greater, the strong mixing and drainage of small fractures accelerate the old water release during high flows during the wet season. It is notable that the young water may sometimes be the most contaminated component contributing to the underground conduit network in karst catchments, because of the direct transfer of contaminants from the ground surface with rain water via large fractures and sinkholes.  相似文献   

17.
Many geothermal anomalies are intersected by vertical fault zones (narrow zones of fractured material with large effective permeability). These conduits are probably responsible for much of the upwelling of hot water from depth. This paper considers a shallow aquifer intersected by a vertical fault. The fluid flow in the aquifer is numerically modeled as a two-dimensional problem. It is observed that the temperature distribution in the aquifer is governed primarily by lateral flow of hot water supplied from the intersecting vertical fault and only secondarily by conduction. The numerical results also provide a possible explanation for the local temperature maxima and inversions occasionally observed in borehole measurements. The present model is an alternative to that based on mushroom-shaped isotherm distributions found in high Rayleigh number large-scale circulation cell calculations.  相似文献   

18.
Geochemical data indicate that the Springfield Plateau aquifer, a carbonate aquifer of the Ozark Plateaus Province in central USA, has two distinct hydrochemical zones. Within each hydrochemical zone, water from springs is geochemically and isotopically different than water from wells. Geochemical data indicate that spring water generally interacts less with the surrounding rock and has a shorter residence time, probably as a result of flowing along discrete fractures and solution openings, than water from wells. Water type throughout most of the aquifer was calcium bicarbonate, indicating that carbonate‐rock dissolution is the primary geochemical process occurring in the aquifer. Concentrations of calcium, bicarbonate, dissolved oxygen and tritium indicate that most ground water in the aquifer recharged rapidly and is relatively young (less than 40 years). In general, field‐measured properties, concentrations of many chemical constituents, and calcite saturation indices were greater in samples from the northern part of the aquifer (hydrochemical zone A) than in samples from the southern part of the aquifer (hydrochemical zone B). Factors affecting differences in the geochemical composition of ground water between the two zones are difficult to identify, but could be related to differences in chert content and possibly primary porosity, solubility of the limestone, and amount and type of cementation between zone A than in zone B. In addition, specific conductance, pH, alkalinity, concentrations of many chemical constituents and calcite saturation indices were greater in samples from wells than in samples from springs in each hydrochemical zone. In contrast, concentrations of dissolved oxygen, nitrite plus nitrate, and chloride generally were greater in samples from springs than in samples from wells. Water from springs generally flows rapidly through large conduits with minimum water–rock interactions. Water from wells flow through small fractures, which restrict flow and increase water–rock interactions. As a result, springs tend to be more susceptible to surface contamination than wells. The results of this study have important implications for the geochemical and hydrogeological processes of similar carbonate aquifers in other geographical locations. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

19.
Based on Helios measurements, seven quantities of normalized PDF (Probability Distribution Function) associated with magnetic field and its disturbances are utilized to characterize the intermittency in the fast solar wind using Castaing distributions and the idea of “Flatness”. The magnetic field fluctuations are found to be more intermittent at farther distances from the sun. The “Flatness” decreases with increasing time scales, with the corresponding PDF eventually approaching Gaussian distributions. Such a transition occurs at a relatively small time scale for the perpendicular component of perturbed field. The increase in “Flatness” with decreasing time scale is more apparent farther from the sun. By examining how the relative energy density of magnetic disturbances at various time scales changes with the mean field, our study supports the idea that the perturbed fields in the fast solar wind in the frequency range considered are consistent with cross-scale redistribution of wave energy favoring larger scales.  相似文献   

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

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

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