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

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

3.
Surface soil hydraulic properties are key factors controlling the partition of rainfall and snowmelt into runoff and soil water storage, and their knowledge is needed for sound land management. The objective of this study was to evaluate the effects of three land uses (native grass, brome grass and cultivated) on surface soil hydraulic properties under near‐saturated conditions at the St Denis National Wildlife Area, Saskatchewan, Canada. For each land use, water infiltration rates were measured using double‐ring and tension infiltrometers at ?0·3, ?0·7, ?1·5 and ?2·2 kPa pressure heads. Macroporosity and unsaturated hydraulic properties of the surface soil were estimated. Mean field‐saturated hydraulic conductivity (Kfs), unsaturated hydraulic conductivity at ?0·3 kPa pressure head, inverse capillary length scale (α) and water‐conducting macroporosity were compared for different land uses. These parameters of the native grass and brome grass sites were significantly (p < 0·1) higher than that of the cultivated sites. At the ?0·3 kPa pressure head, hydraulic conductivity of grasslands was two to three times greater than that of cultivated lands. Values of α were about two times and values of Kfs about four times greater in grasslands than in cultivated fields. Water‐conducting macroporosity of grasslands and cultivated fields were 0·04% and 0·01% of the total soil volume, respectively. Over 90% of the total water flux at ?0·06 kPa pressure head was transmitted through pores > 1·36 × 10?4 m in diameter in the three land uses. Land use modified near‐saturated hydraulic properties of surface soil and consequently may alter the water balance of the area by changing the amount of surface runoff and soil water storage. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

4.
Shah N  Nachabe M  Ross M 《Ground water》2007,45(3):329-338
In many landscapes, vegetation extracts water from both the unsaturated and the saturated zones. The partitioning of evapotranspiration (ET) into vadose zone evapotranspiration and ground water evapotranspiration (GWET) is complex because it depends on land cover and subsurface characteristics. Traditionally, the GWET fraction is assumed to decay with increasing depth to the water table (DTWT), attaining a value of 0 at what is termed the extinction depth. A simple assumption of linear decay with depth is often used but has never been rigorously examined using unsaturated-saturated flow simulations. Furthermore, it is not well understood how to relate extinction depths to characteristics of land cover and soil texture. In this work, variable saturation flow theory is used to simulate GWET for three land covers and a range of soil properties under drying soil conditions. For a water table within half a meter of the land surface, nearly all ET is extracted from ground water due to the close hydraulic connection between the unsaturated and the saturated zones. For deep-rooted vegetation, the decoupling of ground water and vadose zone was found to begin at water table depths between 30 and 100 cm, depending on the soil texture. The decline of ET with DTWT is better simulated by an exponential decay function than the commonly used linear decay. A comparison with field data is consistent with the findings of this study. Tables are provided to vary the extinction depth for heterogeneous landscapes with different vegetation cover and soil properties.  相似文献   

5.
Temperature and moisture content in the variably saturated subsurface are two of the most important physical parameters that govern a wide variety of geochemical and ecological processes. An understanding of thermal and hydraulic processes and properties of transient vadose zones is therefore fundamental in the evaluation of such processes. Here, an investigation of the thermal regime and subsurface properties of a tidally affected, variably saturated streambed is presented. Field and laboratory measurements, as well as a forward numerical model, are jointly employed in the investigation. Temperature, soil moisture, surface level, and water level data were recorded in a transect perpendicular to a tidally driven stream. Frequency‐domain analysis of the subsurface temperature measurements revealed the rapid decay of the tidal temperature driver within the top ～30 cm of sediment. Several techniques were used to evaluate subsurface thermal and hydraulic properties, including thermal conductivity and the soil water retention curve. These properties were used to constrain a forward numerical model that included coupled treatment of relevant variable saturation thermal and hydraulic physics. Even though the investigated vadose zone is intermittent and relatively shallow ( 20 cm), the results illustrate how error can be introduced into heat‐transport calculations if unsaturated conditions are not taken into account.  相似文献   

6.
Previous studies have shown that water retention curve (WRC) and the hydraulic conductivity vary because of changes of the void ratio or porosity of soil. However, limited documents pointed out the change of hydraulic properties of soil when compacted to different porosities while considering both of the drying and wetting processes of the WRC. This information is sometimes necessary for research like finger flow analysis or the occurrence of wetting and drying cycles as what would be seen in the field. Therefore, this study aims to examine the change of WRC characteristics with varied porosity considering both of the drying and wetting path in WRC by conducting a sand box experiment. Results show that the same type of sand compacted to various porosities have different hydraulic parameters. Hydraulic conductivities generally decrease with reduced porosities; shape parameter α of the van Genuchten equation (1980) linearly decreases with declining porosity and shape parameter n in a reversal manner for the sands of interest whether in the drying process or wetting process. The unsaturated properties of sand are further characterized by inspecting the variations of moisture content, matric suction and vertical displacement of soil body subject to periodic changes of the water level by another sand box experiment. The outcomes suggest that the saturated water content and residual water content are changing during the wetting–drying process, which can be an implication of the changed properties of WRC. The characteristics of volumetric deformation might be varied as well because of the observation of the dissimilar patterns of the changing vertical displacements among each wetting–drying process. Infiltration patterns of the sands also are identified through numerical modelling by introducing a constant infiltration flux from the surface followed by a no‐influx condition. Results indicate that less water accumulates in the sand near the surface for the sand compacted to higher porosity, but water can move deeper. Hydraulic conductivity is found as the prime factor dominating the evolvement of wetting fronts. However, shape parameters of water retention curves also affect the infiltration pattern to some extent. In addition, different sands with similar porosities can have quite different infiltrating characteristics. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

7.
In subsurface porous media, the soil water retention curve (WRC) and unsaturated hydraulic conductivity curve (UHC) are two important soil hydraulic property curves. Spatial heterogeneity is ubiquitous in nature, which may significantly affect soil hydraulic property curves. The main theme of this paper is to investigate how spatial heterogeneities, including their arrangements and amounts in soil flumes, affect soil hydraulic property curves. This paper uses a two‐dimensional variably saturated flow and solute transport finite element model to simulate variations of pressure and moisture content in soil flumes under a constant head boundary condition. To investigate the behavior of soil hydraulic property curves owing to variations of heterogeneities and their arrangements as well, cases with different proportions of heterogeneities are carried out. A quantitative evaluation of parameter variations in the van Genuchten model (VG model) resulting from heterogeneity is presented. Results show that the soil hydraulic properties are strongly affected by variations of heterogeneities and their arrangements. If the pressure head remains at a specific value, the soil moisture increases when heterogeneities increase in the soil flumes. On the other hand, the unsaturated hydraulic conductivity decreases when heterogeneities increase in the soil flumes under a constant pressure head. Moreover, results reveal that parameters estimated from both WRC and UHC also are affected by shapes of heterogeneity; this indicates that the parameters obtained from the WRC are not suitable for predicting the UHC of different shapes in heterogeneous media. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

8.
Simultaneous measurement of coupled water, heat, and solute transport in unsaturated porous media is made possible with the multi-functional heat pulse probe (MFHPP). The probe combines a heat pulse technique for estimating soil heat properties, water flux, and water content with a Wenner array measurement of bulk soil electrical conductivity (ECbulk). To evaluate the MFHPP, we conducted controlled steady-state flow experiments in a sand column for a wide range of water saturations, flow velocities, and solute concentrations. Flow and transport processes were monitored continuously using the MFHPP. Experimental data were analyzed by inverse modeling of simultaneous water, heat, and solute transport using an adapted HYDRUS-2D model. Various optimization scenarios yielded simultaneous estimation of thermal, solute, and hydraulic parameters and variables, including thermal conductivity, volumetric water content, water flux, and thermal and solute dispersivities. We conclude that the MFHPP holds great promise as an excellent instrument for the continuous monitoring and characterization of the vadose zone.  相似文献   

9.
The modeling of thick vadose zones is particularly challenging because of difficulties in collecting a variety of measured sediment properties, which are required for parameterizing the model. Some models rely on synthetic data, whereas others are simplified by running as homogeneous sediment domains and relying on a single set of sediment properties. Few studies have simulated flow processes through a thick vadose zone using real and comprehensive data sets comprising multiple measurements. Here, we develop a flow model for a 7-m-thick vadose zone. This model, combining the numerical codes CTRAN/W with SEEP/W, includes the measured sediment hydraulic properties of the investigated vadose zone and incorporates the actual climate and subsurface conditions of the study site (precipitations, water-table elevations, and stable isotope data). The model is calibrated by fitting the simulated and measured vertical profiles of water content. Our flow model calculates a transit time of 1 year for the travel of water through the 7-m vadose zone; this estimate matches stable isotope-based results obtained previously for this site. A homogeneous sediment domain flow model, which considers only a single set of sediment properties, produces a transit time that is approximately half the duration of that of the heterogeneous flow model. This difference highlights the importance of assuming heterogeneous material within models of thick vadose zones and testifies to the advantage gained when using real sediment hydraulic properties to parametrize a flow model.  相似文献   

10.
The “HYDRUS package for MODFLOW” is an existing MODFLOW package that allows MODFLOW to simultaneously evaluate transient water flow in both unsaturated and saturated zones. The package is based on incorporating parts of the HYDRUS-1D model (to simulate unsaturated water flow in the vadose zone) into MODFLOW (to simulate saturated groundwater flow). The coupled model is effective in addressing spatially variable saturated-unsaturated hydrological processes at the regional scale. However, one of the major limitations of this coupled model is that it does not have the capability to simulate solute transport along with water flow and therefore, the model cannot be employed for evaluating groundwater contamination. In this work, a modified unsaturated flow and transport package (modified HYDRUS package for MODFLOW and MT3DMS) has been developed and linked to the three-dimensional (3D) groundwater flow model MODFLOW and the 3D groundwater solute transport model MT3DMS. The new package can simulate, in addition to water flow in the vadose zone, also solute transport involving many biogeochemical processes and reactions, including first-order degradation, volatilization, linear or nonlinear sorption, one-site kinetic sorption, two-site sorption, and two-kinetic sites sorption. Due to complex interactions at the groundwater table, certain modifications of the pressure head (compared to the original coupling) and solute concentration profiles were incorporated into the modified HYDRUS package. The performance of the newly developed model is evaluated using HYDRUS (2D/3D), and the results indicate that the new model is effective in simulating the movement of water and contaminants in the saturated-unsaturated flow domains.  相似文献   

11.
Part of a small drainage basin on the Sevilleta National Wildlife Refuge (about 25 km north of Socorro, NM) was intensively instrumented with soil monitoring equipment to estimate natural ground-water recharge. Soil-moisture data were analysed with special attention to characterizing the influence of topography on the direction of vadose water flow paths in fine to medium aeolian sand. Moisture content data were obtained by the neutron scattering technique, and hydraulic head data were obtained using tensiometers. In addition, tracer experiments were conducted on a sandy hillslope to delineate the flow paths of vadose water. The results indicate that there is a strong lateral component to unsaturated flow on a hillslope, even in the absence of apparent sublayers of much lower permeability. Darcian calculations estimate the long-term, steady, deep flux beneath a concave location to be about 4 per cent of an assumed mean annual precipitation of 20 cm. The deep soil water flux downward varied by several orders of magnitude during the 17 month period of record.  相似文献   

12.
  相似文献   

13.
A borehole permeameter is well suited for testing saturated hydraulic conductivity (K(sat)) at specific depths in the vadose zone. Most applications of the method involve fine-grained soils that allow hand auguring of test holes and require a small water reservoir to maintain a constant head. In non-cohesive gravels, hand-dug test holes are difficult to excavate, holes are prone to collapse, and large volumes of water are necessary to maintain a constant head for the duration of the test. For coarse alluvial gravels, a direct-push steel permeameter was designed to place a slotted pipe at a specific sampling depth. Measurements can be made at successive depths at the same location. A 3790 L (1000 gallons) trailer-mounted water tank maintained a constant head in the permeameter. Head in the portable tank was measured with a pressure transducer and flow was calculated based on a volumetric rating curve. A U.S. Bureau of Reclamation analytical method was utilized to calculate K(sat). Measurements with the permeameter at a field site were similar to those reported from falling-head tests.  相似文献   

14.
After the initiation of gravity drainage, water is often assumed to be either (a) draining under unit gradient, or (b) at capillary/gravity equilibrium. Both of these simplifications can be useful, but the regimes of validity of each assumption must be delineated. Water pressures are measured versus time and distance as water drains out of a 1.6 m long sand column to determine the relative effects of capillary and gravitational forces during drainage. For medium sized sands (0.15–0.3 mm in diameter), the capillary pressure is constant in space in a large region of the column for over 12 days, and the water continues to flow under unit gradient for relatively long time scales. Similar results are seen for finer sands, but with a much faster approach to equilibrium. Numerical simulations and analytical estimates are presented and compare favorably to the measurements. Together, the experimental, theoretical and analytical results are used to calculate when capillary/gravity equilibrium is reached as a function of porous media properties and length of the unsaturated zone. The ratio of the length of the unsaturated zone to the bubbling pressure is a key parameter in determining the drainage regime, and that even for relatively short unsaturated zones the equilibrium time scale can be on the order of years.  相似文献   

15.
Migration of Petroleum Products Through Sandy Hydrogeologic Systems   总被引：5，自引：0，他引：5  
Laboratory column experiments were carried out to study the migration of petroleum products through a sandy porous medium. It was found that the oil pressure needed to displace water from the pores of the sand medium increased with depth below the top of the column. While oil under negative pressure displaced water through the variably unsaturated zone, a significant vertical column of oil was needed to displace water from the pores at the water table. These results indicate that oil penetration to and below the water table will occur only if the porous medium is highly conductive and the rate of oil leak is high. For small to moderate leak rates and hydraulic conductivities, oil would preferentially spread laterally through the zones above the water table. This process of spreading could serve as a natural barrier to severe aquifer contamination by petroleum products.
A simplified procedure was developed to give an order-of-magnitude estimate of the preferred subsurface migration pathways of leaked petroleum products. This procedure utilizes the main drainage and wetting curves for oil and water and the interfacial pressure between these two fluids. This approach could be useful in guiding exploratory investigations, reducing both the risk for further spreading of the contaminants and the investigative cost.  相似文献   

16.
A ground water monitoring program should include an investigation of all possible areas of concern. To be completely effective, the program should include soil sampling, soil analysis and water-quality examination of both the saturated and unsaturated zones. A well-tooled drill rig can take all the proper soil samples, perform all necessary tests and install a functional monitoring well. With the introduction of the fluoropolymer (Teflon(r)) sleeve lysimeter, a single monitoring well can be constructed to monitor both the saturated and unsaturated zones in one installation. The monitoring well screen and casing may also be completely constructed of fluoropolymer.
The sleeve lysimeter is designed with a threaded hollow inner diameter, allowing it to be attached between the joints of a casing string. This hollow I.D. acts as an extension of the casing; the lysimeter surrounds the casing. This creates an isolated vessel for sampling the vadose zone. Access to the screened monitoring well below is unaffected. Tests have shown that when properly installed, these porous fluoropolymer filter units can collect samples with no interaction between the filter and collected fluids.  相似文献   

17.
We present a geostatistically based inverse model for characterizing heterogeneity in parameters of unsaturated hydraulic conductivity for three-dimensional flow. Pressure and moisture content are related to perturbations in hydraulic parameters through cross-covariances, which are calculated to first-order. Sensitivities needed for covariance calculations are derived using the adjoint state sensitivity method. Approximations of the conditional mean parameter fields are then obtained from the cokriging estimator. Correlation between parameters and pressure – moisture content perturbations is seen to be strongly dependent on mean pressure or moisture content. High correlation between parameters and pressure data was obtained under saturated or near saturated flow conditions, providing accurate estimation of saturated hydraulic conductivity, while moisture content measurements provided accurate estimation of the pore size distribution parameter under unsaturated flow conditions.  相似文献   

18.
We present a geostatistically based inverse model for characterizing heterogeneity in parameters of unsaturated hydraulic conductivity for three-dimensional flow. Pressure and moisture content are related to perturbations in hydraulic parameters through cross-covariances, which are calculated to first-order. Sensitivities needed for covariance calculations are derived using the adjoint state sensitivity method. Approximations of the conditional mean parameter fields are then obtained from the cokriging estimator. Correlation between parameters and pressure – moisture content perturbations is seen to be strongly dependent on mean pressure or moisture content. High correlation between parameters and pressure data was obtained under saturated or near saturated flow conditions, providing accurate estimation of saturated hydraulic conductivity, while moisture content measurements provided accurate estimation of the pore size distribution parameter under unsaturated flow conditions.  相似文献   

19.
The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

20.
J. Lewis 《Ground water》2016,54(5):740-744
This technical note describes an effective and inexpensive field technique for measuring the saturated hydraulic conductivity of both undisturbed cores and repacked soil samples. The method requires no specialized equipment; everything that is required can be obtained in a hardware store. The method is a straightforward field implementation of the widely used falling‐head laboratory analysis directly derived from Darcy's law. As such, it sidesteps the need for empirical assumptions about soil texture and the relationship between saturated and unsaturated flow components which many permeameter‐based methods rely upon. The method is shown to produce results that are consistent with K values obtained elsewhere in the same homogeneous sand formation. Furthermore, the proposed method is useful for measuring hydraulic conductivity in drill cuttings obtained from direct push or auguring drill techniques, which cannot be done with any other field method. The range of hydraulic conductivity values that this test is appropriate for is on the order of 1E ? 7 m/s to 1E ? 3 m/s.  相似文献   

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

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