首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
Subsurface drainage systems have been widely used to deal with soil salinization and waterlogging problems around the world. In this paper, a mathematical model was introduced to quantify the transient behavior of the groundwater table and the seepage from a subsurface drainage system. Based on the assumption of a hydrostatic pressure distribution, the model considered the pore-water flow in both the phreatic and vadose soil zones. An approximate analytical solution for the model was derived to quantify the drainage of soils which were initially water-saturated. The analytical solution was validated against laboratory experiments and a 2-D Richards equation-based model, and found to predict well the transient water seepage from the subsurface drainage system. A saturated flow-based model was also tested and found to over-predict the time required for drainage and the total water seepage by nearly one order of magnitude, in comparison with the experimental results and the present analytical solution. During drainage, a vadose zone with a significant water storage capacity developed above the phreatic surface. A considerable amount of water still remained in the vadose zone at the steady state with the water table situated at the drain bottom. Sensitivity analyses demonstrated that effects of the vadose zone were intensified with an increased thickness of capillary fringe, capillary rise and/or burying depth of drains, in terms of the required drainage time and total water seepage. The analytical solution provides guidance for assessing the capillary effects on the effectiveness and efficiency of subsurface drainage systems for combating soil salinization and waterlogging problems.  相似文献   

2.
Analytical solutions for the water flow and solute transport equations in the unsaturated zone are presented. We use the Broadbridge and White nonlinear model to solve the Richards’ equation for vertical flow under a constant infiltration rate. Then we extend the water flow solution and develop an exact parametric solution for the advection-dispersion equation. The method of characteristics is adopted to determine the location of a solute front in the unsaturated zone. The dispersion component is incorporated into the final solution using a singular perturbation method. The formulation of the analytical solutions is simple, and a complete solution is generated without resorting to computationally demanding numerical schemes. Indeed, the simple analytical solutions can be used as tools to verify the accuracy of numerical models of water flow and solute transport. Comparison with a finite-element numerical solution indicates that a good match for the predicted water content is achieved when the mesh grid is one-fourth the capillary length scale of the porous medium. However, when numerically solving the solute transport equation at this level of discretization, numerical dispersion and spatial oscillations were significant.  相似文献   

3.
Changes in groundwater elevation may cause a change in the net normal stress and matric potential within the soil mass, which results in volume changes of unsaturated soil. This research investigated the relationship between the drawdown of groundwater and the characteristics of volumetric compressibility of unsaturated soil. Sand column experiments were designed and conducted to measure the volume changes of coarse and fine sands under different types of drainage conditions at fast and slow drainage rates. The finite element program FEMWATER was calibrated and used to simulate the distributions of stress, tension and moisture content within the sands. Finally, based on the changes of net normal stress and matric potential and the observed volume change of the sands, a least‐square method was applied to determine the volumetric consolidation parameters of the unsaturated soils. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

4.
We derive an analytical model of soil-gas contamination sparged into an imlined unsaturated zone. A nonaqueous phase liquid (NAPL) source lies in the capillary fringe, with an exponential sparge constant within the radius of influence and a constant ambient evaporation rate beyond. Advection, diffusion, and dispersion govern the conservative soil-gas response, expressed as a quasi-steady series solution with radial Bessel and hyperbolic vertical dependence. Simulations suggest that sparged contamination initially spreads beyond the radius of influence down a negative gradient. This gradient eventually reverses, leading to a subsequent influx of ambient contamination. Soil-gas concentrations accordingly reflect slowly varying source conditions as well as slowly varying diffusive transport through the radius of influence. The two time scales are independent: One depends on NAPL, airflow, and capillary fringe characteristics, the other on soil moisture, gaseous diffusivity, and unsaturated zone thickness. The influx of ambient contamination generates an asymptotic soil-gas concentration much less than the initial source concentration. The simple model is applied to a pilot-scale sparging study at Plattsburgh Air Force Base in upstate New York, with physically plausible results.  相似文献   

5.
Estimation of groundwater recharge to an unconfined aquifer is studied using analytical and numerical techniques and results are compared with field observations. There is an acute need for such estimation in water balance studies in arid climates, and the case study in this paper is for such a region. The wetting front movement in the unsaturated zone depends on antecedent soil moisture, the ponded water depth and its duration, and on the position of the water table and the hydraulic properties of the unsaturated zone. A hydraulic connection between the recharge basin and the aquifer is not immediately established because the wetting front is unsaturated. A numerical model is applied to estimate recharge in an arid-zone wadi, and its validity is tested by comparing it with an analytical solution of the equations. The calculated recharge values matched the piezometric levels observed at a well site at the edge of the wadi channel. The total recharge depths found by integration in the time domain provided a good estimate of the transmitted volume of water per unit length of wadi channel. The findings were confirmed by runoff volume measurements at gauging stations located in the basin. © 1997 John Wiley & Sons, Ltd.  相似文献   

6.
A periodic vertical movement of the groundwater table results in a subsequent cyclic response of the water content and pressure profiles in the vadose zone. The sequence of periodic wetting and drying processes can be affected by hysteresis effects in this zone. A one-dimensional saturated/unsaturated flow model based on Richards’ equation and the Mualem (Soil Sci. 137 (1984) 283) hysteresis model is formulated which can take into account multi-cycle hysteresis effects in the relation between capillary pressure and water content. The numerical integration of the unsaturated flow equation is based on a Galerkin-type finite element method. The flow domain is discretised by finite elements with linear shape functions. Simulations start with static water content and pressure profiles, which correspond to either a boundary drying or wetting retention curve. To facilitate the numerical solution of the hysteretic case an implicit non-iterative procedure was chosen for the solution of the nonlinear differential equation. Laboratory experiments were performed with a vertical sand column by imposing a high frequency periodic pressure head at the lower end of the column. The total water volume in the column, and the periodic water content profile averaged over time were measured. The boundary drying and wetting curves of the relation between water content and capillary pressure were determined by independent experiments. The simulations of the experimental conditions show a clear effect of the hysteresis phenomenon on the water content profile. The simulations with hysteresis agree well with the measurements. Computed dimensionless water content profiles are presented for different oscillation frequencies with and without consideration of hysteresis.  相似文献   

7.
The biogeochemical composition of stream water and the surrounding riparian water is mainly defined by the exchange of water and solutes between the stream and the riparian zone. Short-term fluctuations in near stream hydraulic head gradients (e.g., during stream flow events) can significantly influence the extent and rate of exchange processes. In this study, we simulate exchanges between streams and their riparian zone driven by stream stage fluctuations during single stream discharge events of varying peak height and duration. Simulated results show that strong stream flow events can trigger solute mobilization in riparian soils and subsequent export to the stream. The timing and amount of solute export is linked to the shape of the discharge event. Higher peaks and increased durations significantly enhance solute export, however, peak height is found to be the dominant control for overall mass export. Mobilized solutes are transported to the stream in two stages (1) by return flow of stream water that was stored in the riparian zone during the event and (2) by vertical movement to the groundwater under gravity drainage from the unsaturated parts of the riparian zone, which lasts for significantly longer time (> 400 days) resulting in long tailing of bank outflows and solute mass outfluxes. We conclude that strong stream discharge events can mobilize and transport solutes from near stream riparian soils into the stream. The impact of short-term stream discharge variations on solute exchange may last for long times after the flow event.  相似文献   

8.
Volume reduction and lowering of capillary pressure within a large DNAPL pool are utilized as objectives in the design of a large-scale dual phase recovery system at a chemical manufacturing facility in the United States. By reducing DNAPL pool height through mass removal, capillary pressure is lowered, resulting in a reduced potential for future vertical and horizontal mobilization of the chlorinated solvent DNAPL pool. The DNAPL pool extends over an approximately 200 m by 275 m area in low permeability fill deposits overlying a clay aquitard. A three-dimensional multiphase flow model was employed to arrive at a final design incorporating nine horizontal drains (total length 664 m) and a pulsed pumping system. The numerical model was calibrated to the results of a 42-day field pilot-test involving the removal of approximately 25,000 L of DNAPL from a single, 55 m long horizontal drain. Numerical simulation revealed that gravity drainage, as opposed to hydraulic gradients in the water phase, is the dominant recovery mechanism at this site. This stems from the relatively high density and the viscosity of the DNAPL, and the relatively low permeability of the formation deposits. The use of pulsed pumping is shown to reduce the volume of contaminated ground water recovered from the 9-drain system, without significant reduction of the total volume of DNAPL recovered.  相似文献   

9.
An analytical solution is presented for the slug tests conducted in a partially penetrating well in an unconfined aquifer affected from above by an unsaturated zone. The solution considers the effects of wellbore skin and oscillatory responses on underdamped slug tests. The flow in the saturated zone is described by a two‐dimensional, axially symmetric governing equation, and the flow in the unsaturated zone above the water table by a linearized one‐dimensional Richards' equation. The unsaturated medium properties are represented by the exponential constitutive relationships. A Laplace domain solution is derived using the Laplace and finite Fourier transform and the solution in the real‐time domain is evaluated using the numerical inverse Laplace transform method. The solution derived in this study is more general and reduces to the most commonly used solutions for slug tests in their specified conditions. It is found that the unsaturated flow has a significant impact on the slug test conducted in an unconfined aquifer. The impact of unsaturated flow on such a slug test is enhanced with a larger anisotropy ratio, a shorter well screen length, a shorter distance between the well screen and the water table, or a larger well screen radius. The impact of unsaturated flow on slug tests decreases as the degree of penetration (the length of well screen) increases. For a fixed well screen length, the impact of unsaturated flow on slug tests decreases as the distance between the centre of screen and the water table increases. A large dimensionless well screen radius (>0.01) leads to significant effects of unsaturated flow on slug tests. The unsaturated flow reduces the oscillatory responses to underdamped slug tests. The unsaturated zone has significant impact on slug test under high‐permeability wellbore skin.  相似文献   

10.
Evaluation of volatilization as a natural attenuation pathway for MTBE   总被引：2，自引：0，他引：2  
Lahvis MA  Baehr AL  Baker RJ 《Ground water》2004,42(2):258-267
Volatilization and diffusion through the unsaturated zone can be an important pathway for natural attenuation remediation of methyl tert-butyl ether (MTBE) at gasoline spill sites. The significance of this pathway depends primarily on the distribution of immiscible product within the unsaturated zone and the relative magnitude of aqueous-phase advection (ground water recharge) to gaseous-phase diffusion. At a gasoline spill site in Laurel Bay, South Carolina, rates of MTBE volatilization from ground water downgradient from the source are estimated by analyzing the distribution of MTBE in the unsaturated zone above a solute plume. Volatilization rates of MTBE from ground water determined by transport modeling ranged from 0.0020 to 0.0042 g m(-2)/year, depending on the assumed rate of ground water recharge. Although diffusive conditions at the Laurel Bay site are favorable for volatilization, mass loss of MTBE is insignificant over the length (230 m) of the solute plume. Based on this analysis, significant volatilization of MTBE from ground water downgradient from source areas at other sites is not likely. In contrast, model results indicate that volatilization coupled with diffusion to the atmosphere could be a significant mass loss pathway for MTBE in source areas where residual product resides above the capillary zone. Although not documented, mass loss of MTBE at the Laurel Bay site due to volatilization and diffusion to the atmosphere are predicted to be two to three times greater than mass loading of MTBE to ground water due to dissolution and recharge. This result would imply that volatilization in the source zone may be the critical natural attenuation pathway for MTBE at gasoline spill sites, especially when considering capillary zone limitations on volatilization of MTBE from ground water and the relative recalcitrance of MTBE to biodegradation.  相似文献   

11.
Dunn AM  Silliman SE 《Ground water》2003,41(6):729-734
A laboratory tank was used to study entrapment of water in coarse sand lenses above the water table and of air in coarse sand lenses below the water table. Monitoring of these experiments involved a combination of visual inspection, measurement of moisture content, and measurement of air/water pressure. The medium consisted of coarse sand lenses with various degrees of vertical connectivity embedded within a fine sand matrix. Experiments were performed under conditions of both drainage (from a fully saturated medium) and imbibition. Observations during drainage included: (1) water was trapped in the coarse sand zones above the water table at heights significantly greater than anticipated from consideration of capillary rise in the coarse sand; (2) rapid drainage of these same coarse zones occurred when air penetrated into these zones through the surrounding fine sands; and (3) prior to the time of penetration of the coarse sand by air, water pressure in the coarse zone dropped significantly below atmospheric pressure. Observations during imbibition included: (1) entrapment of air within coarse sands below the water table, (2) the pore fluids in these zones varied spatially from predominantly air to predominantly water, and (3) pressure in the trapped air phase was significantly greater than pressure in the water phase in the surrounding fine sand. Overall, these results demonstrated significant sensitivity to the geometry of the coarse sand inclusions, particularly the vertical connectivity of the coarse sand lens.  相似文献   

12.
For certain porous media and initial conditions, constant flux infiltrations show a saturation profile which exhibits overshoot. This overshoot is the cause of gravity driven fingering, cannot be described by standard models of unsaturated flow, and is likely controlled by the exact nature of the pore filling at the initial front. Here we report synchrotron X-ray micro-tomography measurements of the porous medium and measure which pores are filled by water and air at the initial wetting front as a function of flux. We find that at high fluxes all the pores are filled with water; for intermediate fluxes, the pores along the edge of the column remain unsaturated; and for low fluxes the pores in the bulk of the experimental column remain unsaturated. This suggests that the unsaturated overshoot conditions observed at higher fluxes are primarily an edge effect of the column. The results can help delineate the correct continuum models that can capture overshoot and gravity driven fingering.  相似文献   

13.
Systematic experimental deviations from theoretical predictions derived for water retention characteristics of fractal porous media have previously been interpreted in terms of continuum percolation theory (at low moisture contents, below the critical volume fraction of water, αc capillary flow ceases). In other work, continuum percolation theory was applied to find the hydraulic conductivity as a function of saturation for saturations high enough to guarantee percolation of capillary flow. Now these two problems are further linked, using percolation theory to estimate non-equilibrium water retention at matric potential values such that the equilibrium water content is too low for percolation of capillary flow paths. In particular, a procedure for developing a time-dependent moisture content is developed for experimental time scales long enough that film flow can provide an alternate mechanism for equilibrating when continuous capillary flow is not possible. The time scales are defined in terms of moisture-dependent length scales and film flow and capillary flow hydraulic conductivities. Imbibition is treated in the extreme case of no film-flow contribution to equilibration. In another application at higher matric potentials, recursive relations are derived for the water content of porous media during drying when external pressures are changed at rates too rapid for equilibrium to be attained by capillary flow.  相似文献   

14.
  相似文献   

15.
Monitoring of the vadose zone is a potentially complex, time-consuming, and expensive problem. The location of monitoring points and selection of monitoring instruments can be optimized by using computer models. Numerical models developed for this purpose in the past have often been expensive and difficult to use. This paper describes a fast, three-dimensional, approximate analytical solution to the moisture content in the unsaturated zone. An analytical solution is available for steady-state drainage, whereas an approximate analytical solution is available for the transient case. The model will handle an arbitrary distribution of fluid sources, as well as vertical and horizontal impermeable boundaries.
The model may be applied to predict the incursion of fluid from accidental leakage or infiltration over large areas from unlined ponds and land treatment sites. The model is quite useful as an aid in designing monitoring or premonitoring programs near hazardous waste sites. Examples are presented to demonstrate the model's utility in estimating the maximum spread of a contaminant, the extent to which the fluid may spread with depth, the regions of high and low capillary pressure, and the non-linear behavior of the saturation when drainage from several sources in considered. These results are useful for the placement of monitoring locations and the selection of appropriate instruments, and as a tool in working with regulatory agencies to design monitoring programs. A glimpse of the future is necessary for today's planning. Computer models are some of the most useful crystal balls we have available.  相似文献   

16.
A simple process‐based approach to predict regional‐scale loading of nitrate at the water table was implemented in a GIS for Great Britain. This links a nitrate input function, unsaturated zone thickness, and lithologically dependent rate of nitrate unsaturated zone travel to estimate arrival time of nitrate at the water table. The nitrate input function is the loading at the base of the soil and has been validated using unsaturated zone porewater profiles. The unsaturated zone thickness uses groundwater levels based on regional‐scale observations infilled by interpolated river base levels. Estimates of the rate of unsaturated zone travel are attributed from regional‐scale hydrogeological mapping. The results indicate that peak nitrate loading may have already arrived at the water table for many aquifers, but that it has not where the unsaturated zone is relatively thick There are contrasting outcomes for the two main aquifers which have similar unsaturated zone velocities, the predominantly low relief Permo‐Triassic sandstones, and the Chalk, which forms significant topographic features. For about 60% of the Chalk, the peak input has not yet reached the water table and will continue to arrive over the next 60 years. The methodology is readily transferable and provides a robust method for estimating peak arrival time for any diffuse conservative pollutant where an input function can be defined at a regional scale and requires only depth to groundwater and a hydrogeological classification. The methodology is extendable in that if additional information is available this can easily be incorporated into the model scheme. British Geology Survey © NERC 2011. Hydrological Process © 2011 John Wiley & Sons, Ltd  相似文献   

17.
An important quantity in groundwater protection is the residence time of water in an aquifer. It relates to both the travel time of a pollutant to arrive at a well and the time span required for self-purification of a polluted aquifer after removal of pollutant inputs. Time scales for aquifers can be gained from artificial tracer experiments or from environmental tracer data, the latter offering the only realistic alternative if time scales of years or decades have to be taken into account.

Different tracers show different time scales due to their different transport mechanisms especially in the unsaturated zone. While solute tracers are moved advectively with the seepage water, gas tracers pass the unsaturated zone diffusively through the air phase. Depending on the properties of the unsaturated zone (hydraulic properties, thickness) this difference in behavior can be used to separate the subsurface transport process into the unsaturated and the saturated parts.

In a field study in Germany, SF6 and 3H were used as environmental tracers. Both have a relatively well-known input function. Interpretation of data from observation wells by a box model approach led to spatially and temporally varying residence times. This was an indication that the influence of the unsaturated zone could not be neglected. While the gas tracer SF6 shows only residence times in the saturated zone, the tracer 3H reflects the whole travel time of water including both the unsaturated and saturated zones. Using a one-dimensional plug-flow model for the unsaturated zone combined with a detailed two-dimensional flow and transport model for the saturated zone leads to a holistic and consistent interpretation of the measured tracer concentrations. The observed pattern of old water under thick loess cover and younger water under areas where the fractured basalt aquifer crops out is reproduced after adjusting only two parameters: the effective porosity of the saturated aquifer and the product of field capacity and thickness of the unsaturated zone. While the effective porosity of the saturated zone is adjusted by means of the SF6 data, the field capacity of the loess layer is adjusted by means of the 3H observations. The thickness of the unsaturated zone is deduced from geological and pedological maps. All flow data are obtained from a calibrated flow model, which is based on geological data, observed heads and pumping tests only.

The transport model for the saturated zone was calibrated by fitting the porosity by means of gaseous tracer concentrations (SF6). The combined saturated–unsaturated zone model was then calibrated by fitting the field capacity of the unsaturated zone by means of 3H concentrations. With this model it was possible to verify the observed NO3 concentrations at the drinking water wells and to develop predictions for their future development under various scenarios of fertilizer input reduction in specific areas.  相似文献   


18.
The hydraulic gradient between aquifers and rivers is one of the most variable properties in a river/aquifer system. Detailed process understanding of bank storage under hydraulic gradients is obtained from a two‐dimensional numerical model of a variably saturated aquifer slice perpendicular to a river. Exchange between the river and the aquifer occurs first at the interface with the unsaturated zone. The proportion of total water exchanged through the river bank compared to the river bed is a function of aquifer hydraulic conductivity, partial penetration, and hydraulic gradient. Total exchange may be estimated to within 50% using existing analytical solutions provided that unsaturated zone processes do not strongly influence exchange. Model‐calculated bank storage is at a maximum when no hydraulic gradient is present and increases as the hydraulic conductivity increases. However, in the presence of a hydraulic gradient, the largest exchange flux or distance of penetration does not necessarily correspond to the highest hydraulic conductivity, as high hydraulic conductivity increases the components of exchange both into and out of an aquifer. Flood wave characteristics do not influence ambient groundwater discharge, and so in large floods, hydraulic gradients must be high to reduce the volume of bank storage. Practical measurement of bank storage metrics is problematic due to the limitations of available measurement technologies and the nested processes of exchange that occur at the river‐aquifer interface. Proxies, such as time series concentration data in rivers and groundwater, require further development to be representative and quantitative.  相似文献   

19.
We performed quasi‐two‐dimensional flow through laboratory experiments to study the effect of a coarse‐material inclusion, located in the proximity of the water table, on flow and oxygen transfer in the capillary fringe. The experiments investigate different phases of mass transfer from the unsaturated zone to anoxic groundwater under both steady‐state and transient flow conditions, the latter obtained by fluctuating the water table. Monitoring of flow and transport in the different experimental phases was performed by visual inspection of the complex flow field using a dye tracer solution, measurement of oxygen profiles across the capillary fringe, and determination of oxygen fluxes in the effluent of the flow‐through chamber. Our results show significant effects of the coarse‐material inclusion on oxygen transfer during the different phases of the experiments. At steady state, the oxygen flux across the unsaturated/saturated interface was considerably enhanced due to flow focusing in the fully water‐saturated coarse‐material inclusion. During drainage, a zone of higher water saturation formed in the fine material overlying the coarse lens. The entrapped oxygen‐rich aqueous phase contributed to the total amount of oxygen supplied to the system when the water table was raised back to its initial level. In case of imbibition, pronounced air entrapment occurred in the coarse lens, causing oxygen to partition between the aqueous and gaseous phases. The oxygen mass supplied to the anoxic groundwater following the imbibition event was found to be remarkably higher (approximately seven times) in the heterogeneous system compared with a similar experiment performed in a homogeneous porous medium.  相似文献   

20.
The finite-element method based on a Galerkin technique was used to formulate the problem of simulating the two-dimensional (cross-sectional) transient movement of water and solute in saturated or partially saturated nonuniform porous media. The numerical model utilizes linear triangular elements. Nonreactive, as well as reactive solutes whose behaviour can be described by a distribution coefficient or first-order reaction term were considered. The flow portion of the model was tested by comparison of the model results with experimental and finite-difference results for transient flow in an unsaturated sand column and the solute transport portion of the model was tested by comparison with analytical solution results. The model was applied to a hypothetical case involving movement of water and solutes in tile-drained soils. The simulation results showed the development of distinct solute leaching patterns in the soil as drainage proceeded. Although applied to a tile drainage problem in this study, the model should be equally useful in the study of a wide range of two-dimensional water and solute migration problems.  相似文献   

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

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