首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 958 毫秒
1.
As pumping proceeds in an unconfined aquifer the water-table is lowered and slowly drains out of the initially saturated material giving a delayed yield from storage. Although considerable work has been carried out on delayed yield, a definitive solution is still lacking. The unsaturated zone was treated here as a significant factor in delayed yield analysis. The one-dimensional unsaturated flow equation was solved numerically to simulate the flow in the zone located between successive water levels in the unconfined aquifer for the case with zero flux at the upper boundary. The moisture content distribution curve was integrated numerically to determine the amount of water drained with time. The study shows that the numerical results are in close agreement with the monitored field data, which implies the important role of the unsaturated zone and the reliability of the numerical simulation in describing this physical phenomenon. The results show that the relationship between the flux and time is sensibly in exponential form, which agrees with the general applicability of Boulton's assumption.  相似文献   

2.
Equations are derived for the flow to a pumped well in an aquifer having uniform anisotropy and overlain by a low-permeability aquitard. The water-table is assumed to be located in the aquitard. Drainage from the capillary zone above the water-table is taken into account.The differential equation for the flow in the aquifer is identical with that derived in a previous paper. The formation constants may therefore be evaluated by using type curves as described in that paper.A well-known pumping test is reanalysed, using the equations in the present paper. It is shown that the time-drawdown curves can be explained only by the existence of a low-permeability stratum in the vicinity of the water-table. In this example the slow draining of the unsaturated zone above the water-table seems to be a significant factor in determining the shape of the time-drawdown curves.  相似文献   

3.
In this second of two papers, analytical step-response functions, developed in the companion paper for several cases of transient hydraulic interaction between a fully penetrating stream and a confined, leaky, or water-table aquifer, are used in the convolution integral to calculate aquifer heads, streambank seepage rates, and bank storage that occur in response to stream-stage fluctuations and basinwide recharge or evapotranspiration. Two computer programs developed on the basis of these step-response functions and the convolution integral are applied to the analysis of hydraulic interaction of two alluvial stream–aquifer systems in the northeastern and central United States. These applications demonstrate the utility of the analytical functions and computer programs for estimating aquifer and streambank hydraulic properties, recharge rates, streambank seepage rates, and bank storage. Analysis of the water-table aquifer adjacent to the Blackstone River in Massachusetts suggests that the very shallow depth of water table and associated thin unsaturated zone at the site cause the aquifer to behave like a confined aquifer (negligible specific yield). This finding is consistent with previous studies that have shown that the effective specific yield of an unconfined aquifer approaches zero when the capillary fringe, where sediment pores are saturated by tension, extends to land surface. Under this condition, the aquifer's response is determined by elastic storage only. Estimates of horizontal and vertical hydraulic conductivity, specific yield, specific storage, and recharge for a water-table aquifer adjacent to the Cedar River in eastern Iowa, determined by the use of analytical methods, are in close agreement with those estimated by use of a more complex, multilayer numerical model of the aquifer. Streambank leakance of the semipervious streambank materials also was estimated for the site. The streambank-leakance parameter may be considered to be a general (or lumped) parameter that accounts not only for the resistance of flow at the river–aquifer boundary, but also for the effects of partial penetration of the river and other near-stream flow phenomena not included in the theoretical development of the step-response functions.  相似文献   

4.
The one-dimensional pressure head profile above a fixed water-table was studied for different steady infiltration rates. As shown in previous studies, when the infiltration rate (qin) is less than the soil saturated hydraulic conductivity (Ks), this profile can be divided into two parts: (1) from the water-table surface (z0) to an elevation zγ, the pressure head varies from 0 to a value hγ such as K(hγ)=−qin; (2) above the elevation zγ the pressure head is constant and equal to hγ. Above the water-table the zone where the pressure head is variable has been called ‘transition zone’. Its height is shown to be highly variable and to depend on soil properties as well as on the infiltration rate. This transition zone is not the ‘capillary fringe’ as defined by Gillham (Gillham R.W., 1984. The capillary fringe and its effect on water-table response. J. Hydrol. 67, 307–324). Numerical experiments performed with HYDRUS-2D® for the case of a water-table drained by parallel ditches have shown that the height of the transition zone is similar in the one-dimensional profile and in the two-dimensional system as long as the local slope of the water-table is small. This result is important since in a two-dimensional system, the transition zone is the place where all the horizontal component of the unsaturated flow occurs. The ratio of the horizontal component of the unsaturated flow vs the total horizontal flow in both the unsaturated and saturated zones has been computed. For a given soil, this ratio decreases as the infiltration rate increases; for a given infiltration rate, the soil with the thinnest transition zone transfers the largest amount of water above the water-table.  相似文献   

5.
  相似文献   

6.
Mexico City is situated in the Valley of Mexico on the extensive lacustrine clays that overlay highly productive aquifers of both volcanic and sedimentary origin. The Valley is closed by volcanic mountains. The natural hydraulic boundary conditions associated withe mountain ranges and their relationship to the important aquifers were studied using a two-dimensional, steady-state finite-element model in cross section. Four cross sections were analysed under hydrologic conditions existing prior to the large scale pumping of the aquifers. Factors such as bulk hydraulic conductivities and regional infiltration rates were obtained from field observations and the literature to assess location of the associated groundwater divides, and the water-table in the mountains. The modeled flow patterns are consistent with the historical hydrologic records piezometric characteristics and observed surface features of the groundwater in the Basin of Mexico. From the modeling results, the groundwater recharge in the mountains is 30–50% of the mean average precipitation. Higher and lower rates result in a flow regime that is not compatible with field observations. In general the location of the divides in the mountains is displaced towards the Valley of Mexico, which influences the groundwater budget of the Valley. The water table in places is several hundred metres below ground surface, in accordance with field observations of a very thick unsaturated zone. Before major aquifer exploitation began about 50 years ago, 40–50% of the total discharge into the Valley was by upward flow through the lacustrine deposits. The best results were obtained using a subsurface distribution of hydrostratigraphic units based on recently published geological interpretations.  相似文献   

7.
This paper describes a stochastic analysis of steady state flow in a bounded, partially saturated heterogeneous porous medium subject to distributed infiltration. The presence of boundary conditions leads to non-uniformity in the mean unsaturated flow, which in turn causes non-stationarity in the statistics of velocity fields. Motivated by this, our aim is to investigate the impact of boundary conditions on the behavior of field-scale unsaturated flow. Within the framework of spectral theory based on Fourier–Stieltjes representations for the perturbed quantities, the general expressions for the pressure head variance, variance of log unsaturated hydraulic conductivity and variance of the specific discharge are presented in the wave number domain. Closed-form expressions are developed for the simplified case of statistical isotropy of the log hydraulic conductivity field with a constant soil pore-size distribution parameter. These expressions allow us to investigate the impact of the boundary conditions, namely the vertical infiltration from the soil surface and a prescribed pressure head at a certain depth below the soil surface. It is found that the boundary conditions are critical in predicting uncertainty in bounded unsaturated flow. Our analytical expression for the pressure head variance in a one-dimensional, heterogeneous flow domain, developed using a nonstationary spectral representation approach [Li S-G, McLaughlin D. A nonstationary spectral method for solving stochastic groundwater problems: unconditional analysis. Water Resour Res 1991;27(7):1589–605; Li S-G, McLaughlin D. Using the nonstationary spectral method to analyze flow through heterogeneous trending media. Water Resour Res 1995; 31(3):541–51], is precisely equivalent to the published result of Lu et al. [Lu Z, Zhang D. Analytical solutions to steady state unsaturated flow in layered, randomly heterogeneous soils via Kirchhoff transformation. Adv Water Resour 2004;27:775–84].  相似文献   

8.
The unsteady water-table movement caused by intermittent rainfall and varying evaporation in flat lands intersected by a network of ditches is modelled using land-drainage theory. The unsteady water tables are assumed to behave as a continuous succession of steady states with the flux through the water table given by the sum of components due to rainfall and evaporation through the soil surface and due to water released or taken up by the unsaturated soil above the water table. A simple steady-state drainage equation is used for the relationship between water-table height and flux, and the specific yield is assumed to have a constant value. The simulated seasonal water table using estimated hydraulic soil properties and meteorological records for a field site agreed with available dip-well observations. The water table was much lower than the ditch-water level during the summer months. The sensitivity of simulated results to model parameters is demonstrated.  相似文献   

9.
Data are presented of vertical profiles of soil water potential and soil temperature from a field lysimeter experiment with automatic water-table control at Silwood Park, Ascot, UK. The observed water-table fluxes confirm the analysis of profile response, which shows a diurnal pattern of upward water movement as a result of surface freezing, followed by reversed flow when soil thawing occurs. The results provide a clear, and relatively rare, illustration of the phenomenon of freezing-induced redistribution.  相似文献   

10.
A numerical algorithm is described for solving the free-surface groundwater flow equations in 3-D large-scale unconfined aquifers with strongly heterogeneous conductivity and surface recharge. The algorithm uses a moving mesh to track the water-table as it evolves according to kinematic and seepage face boundary conditions. Both steady-state and transient algorithms are implemented in the SECO-Flow 3-D code and demonstrated on stratigraphy based on the Delaware Basin of south-eastern New Mexico.  相似文献   

11.
The present paper describes an approach to modelling the unsaturated soil-moisture zone in the framework of an integrated physically-based hydrologic response model. It is supposed that the subsurface flow regime may be viewed as two separate entities — a saturated flow system which may be modelled by standard two-dimensional regional techniques, and a single overlying unsaturated zone in which the flow is essentially vertical. Coupling takes place via the definition of saturation at the lower boundary of the unsaturated zone, and via a conservative water balance. Attention is focused on the computational procedure for the unsaturated zone as a self-contained module. The major difficulties are the definition of the interface between the saturated and unsaturated zones, the nonlinear character of the equation used to describe unsaturated flow, the inclusion of realistic atmospheric boundary conditions, and, the interaction between water uptake by plants and available soil-moisture. Each of these points is discussed, in turn, with the emphasis on mathematically formulating the problem in such a way that the most important physical features are reproduced with a minimal amount of computational effort. The text concludes with a few illustrative examples.  相似文献   

12.
This work presents a rigorous numerical validation of analytical stochastic models of steady state unsaturated flow in heterogeneous porous media. It also provides a crucial link between stochastic theory based on simplifying assumptions and empirical field and simulation evidence of variably saturated flow in actual or realistic hypothetical heterogeneous porous media. Statistical properties of unsaturated hydraulic conductivity, soil water tension, and soil water flux in heterogeneous soils are investigated through high resolution Monte Carlo simulations of a wide range of steady state flow problems in a quasi-unbounded domain. In agreement with assumptions in analytical stochastic models of unsaturated flow, hydraulic conductivity and soil water tension are found to be lognormally and normally distributed, respectively. In contrast, simulations indicate that in moderate to strong variable conductivity fields, longitudinal flux is highly skewed. Transverse flux distributions are leptokurtic. the moments of the probability distributions obtained from Monte Carlo simulations are compared to modified first-order analytical models. Under moderate to strong heterogeneous soil flux conditions (σ2y≥1), analytical solutions overestimate variability in soil water tension by up to 40% as soil heterogeneity increases, and underestimate variability of both flux components by up to a factor 5. Theoretically predicted model (cross-)covariance agree well with the numerical sample (cross-)covarianaces. Statistical moments are shown to be consistent with observed physical characteristics of unsaturated flow in heterogeneous soils.©1998 Elsevier Science Limited. All rights reserved  相似文献   

13.
Groundwater flows and travel times were analyzed under a landfill site characterized by a thick unsaturated zone, in an arid environment, using two-dimensional models. Transport in the unsaturated zone is important in arid environments with thick vadose zones and usually determines the travel time of groundwater from the landfill site to the regulatory compliance surface. The combined use of a parallel flow (Dupuit-Forchheimer) areal model and vertical cross-sectional model is explored at a site typical in the utility industry. Flow paths and travel times to the compliance zone are obtained for use in landfill site design or environmental impact assessments. Sensitivity analysis shows that boundary conditions are important factors influencing the entire analysis and thus need to be confirmed with field measurements.  相似文献   

14.
Groundwater elevation fluctuation has been recognized as one mechanism causing temporal indoor air volatile organic chemical (VOC) impacts in vapor intrusion risk assessment guidance. For dissolved VOC sources, groundwater table fluctuation shortens/lengthens the transport pathway, and delivers dissolved contaminants to soils that are alternating between water saturated and variably saturated conditions, thereby enhancing volatilization potential. To date, this mechanism has not been assessed with field data, but enhanced VOC emission flux has been observed in lab-scale and modeling studies. This work evaluates the impact of groundwater elevation changes on VOC emission flux from a dissolved VOC plume into a house, supplemented with modeling results for cyclic groundwater elevation changes. Indoor air concentrations, air exchange rates, and depth to groundwater (DTW) were collected at the study house during an 86-d constant building underpressurization test. These data were used to calculate changes in trichloroethylene (TCE) emission flux to indoor air, during a period when DTW varied daily and seasonally from about 3.1 to 3.4 m below the building foundation (BF). Overall, TCE flux to indoor air varied by about 50% of the average, without any clear correlation to changes in DTW or its change rate. To complement the field study, TCE surface emission fluxes were simulated using a one-dimensional model (HYDRUS 1D) for conditions similar to the field site. Simulation results showed time-averaged surface TCE fluxes for cyclic water-table elevations were greater than for stationary water-table conditions at an equivalent time-averaged water-table position. The magnitudes of temporal TCE emission flux changes were generally less than 50% of the time-averaged flux, consistent with the field site observations. Simulation results also suggested that TCE emission flux changes due to groundwater fluctuation are likely to be significant at sites with shallow groundwater (e.g., < 0.5 m BF) and permeable soil types (e.g., sand).  相似文献   

15.
A field investigation of unsaturated flow through a lithophysal unit of fractured welded tuff containing lithophysal cavities has been initiated. To characterize flow in this spatially heterogeneous medium, a systematic approach has been developed to perform tests in boreholes drilled at regular intervals in an underground tunnel (drift). The purpose of the testing is to quantify the amounts of water seeping into the drift versus the amount of water moving around the drift when released into boreholes at many equidistant locations along the drift. In this paper, we describe the test equipment system that has been built for this purpose. Because the field-scale measurements--of liquid flow in the unsaturated, fractured rocks--require continuous testing for periods of days to weeks, the control of test equipment has been fully automated, allowing operation with no human presence at the field site. Preliminary results from the first set of tests indicate that, while the effects of evaporation on characterization of hydrological properties of the rock can be significant, these effects can be controlled and quantified. These tests give insight into the role of the cavities as potential storage during the initial transient flow prior to the breakthrough of water at the drift crown, as well as the role of connected fractures that provide the subsequent quasi-steady flow. In addition to the stated purpose of realizing the flow partitioning, the results yield values for the effective porosity in the pathways for liquid flow in the regions tested thus far.  相似文献   

16.
Abstract

A resistivity survey in Wenner arrangement has been conducted in the water-table aquifer in the eastern part of the Coastal Area of Belgium. Bore holes have been drilled to test the geoelectrical data and also to collect water samples. The chemical characteristics and the distribution of different water types explain the hydrochemistry of the water-table aquifer. The fresh-/brackish-water boundary as determined by the resistivity survey has been expressed in terms of total dissolved-solids content.  相似文献   

17.
A second order stochastic differential equation is used for modeling of water-table elevation. The data were sampled at the Borden Aquifer as a part of a tracer experiment. The purpose of the water-table data collection was to determine presence of a water flow. We argue that the water-table surface is a simple plane oscillating up and down in time according to an equation for a stochastic oscillator. We derive the model, estimate its parameters and provide arguments for goodness-of-fit of the model.  相似文献   

18.
A second order stochastic differential equation is used for modeling of water-table elevation. The data were sampled at the Borden Aquifer as a part of a tracer experiment. The purpose of the water-table data collection was to determine presence of a water flow. We argue that the water-table surface is a simple plane oscillating up and down in time according to an equation for a stochastic oscillator. We derive the model, estimate its parameters and provide arguments for goodness-of-fit of the model.  相似文献   

19.
The dynamics (from rotation to rotation) of the absolute values of the large-scale open solar magnetic field fluxes in the four-sector field structure has been considered for the first time, using CRs 2032–2035 in July–October 2005 as examples. An important role of the ratio of the fluxes at the eastern and western sector boundaries (Φ E /Φ W ) is confirmed. As in the cases of the two-sector structure, Φ E /Φ W > 1 is typical of active rigidly corotating boundaries with intense sunspot formation, flares, and interplanetary and geomagnetic disturbances. A remarkable property of the considered structure was the presence of a rapidly increasing flux in an initially narrow sector and the flux interaction with a stable rigidly corotating sector in the zone of the main active longitudes, which caused an unexpectedly strong geoeffective long-range action of flares near the corresponding active boundary.  相似文献   

20.
Non-unique solutions of inverse problems arise from a lack of information that satisfies necessary conditions for the problem to be well defined. This paper investigates these conditions for inverse modeling of water flow through multi-dimensional variably saturated porous media. It shows that in order to obtain a unique estimate of hydraulic parameters, along each streamline of the flow field (1) spatial and temporal head observations must be given; (2) the number of spatial and temporal head observations required should be greater or equal to the number of unknown parameters; (3) the flux boundary condition or the pumping rate of a well must be specified for the homogeneous case and both boundary flux and pumping rate are a must for the heterogeneous case; (4) head observations must encompass both saturated and unsaturated conditions, and the functional relationships for unsaturated hydraulic conductivity/pressure head and for the moisture retention should be given, and (5) the residual water content value also need to be specified a priori or water content measurements are needed for the estimation of the saturated water content.For field problems, these necessary conditions can be collected or estimated but likely involve uncertainty. While the problems become well defined and have unique solutions, the solutions likely will be uncertain. Because of this uncertainty, stochastic approaches are deemed to be appropriate for inverse problems as they are for forward problems to address uncertainty. Nevertheless, knowledge of these necessary conditions is critical to reduce uncertainty in both characterization of the vadose zone and the aquifer, and prediction of water flow and solute migration in the subsurface.  相似文献   

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

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