首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 921 毫秒
1.
Although the importance to account for microrelief in the calculation of specific yields for shallow groundwater systems is well recognized, the microrelief influence is often treated very simplified, which can cause considerable errors. We provide a general one‐dimensional expression that correctly represents the effect of a microrelief on the total specific yield that is composed of the soil and surface specific yield. The one‐dimensional expression can be applied for different soil hydraulic parameterizations and soil surface elevation frequency distributions. Applying different van Genuchten parameters and a simple linear microrelief model, we demonstrate that the specific yield is influenced by the microrelief not only when surface storage directly contributes to specific yield by (partial) inundation but also when water levels are lower than the minimum surface elevation. Compared with a simplified representation of the soil specific yield, in which a mean soil surface is assumed for the calculation of soil specific yield, the correct representation can lead to lower as well as higher soil specific yields depending on the specific interaction of the soil water retention characteristics and the microrelief. The new equation can be used to obtain more accurate evapotranspiration estimates from water level fluctuations and to account for the effect of microtopographic subgrid variability on simulated water levels of spatially distributed hydrological models. Copyright © 2015 The Authors Hydrological Processes Published by John Wiley & Sons Ltd  相似文献   

2.
The ultimate specific yield, i.e. the specific yield under steady‐state condition, is generally dependent on the water table depth in the shallow groundwater systems. This paper presents a new expression to describe this relationship. The expression is based on the relationship between soil water content–soil depths below ground surface near the saturation zone, i.e. water content profile, which can be fitted by the van Genuchten model. Because this model is highly nonlinear, a Taylor series expansion was used to solve the integral equation for specific yield in our study. The new expression was verified by a drainage experiment, and the results demonstrated that the new expression is reliable for eolian sand when the range of water table fluctuation is less than 30 cm. The expression is easily applicable because knowledge of the water content profile near the saturation zone is the only requirement. Compared with the expression obtained from the Brooks–Corey model, the new expression can be used for a variety of aquifer and soil media. The improved accuracy of the specific yield provides a better estimate of discharge rates in shallow groundwater systems with water table fluctuations. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

3.
The analytic element method is well suited for the Gardner hydraulic conductivity function, but is limited in describing real soils. Therefore, parameter equivalence between the van Genuchten and Gardner hydraulic conductivity functions is explored for the case of steady vertical flow through a homogeneous medium with a single inclusion, i.e., a binary soil. The inclusion has different hydraulic parameters than the background medium. Equivalence is established using three methods: (1) effective capillary drive; (2) capillary length; (3) and a least-squares optimization method that aims to fit a Gardner function to a corresponding van Genuchten function by minimizing the difference in log conductivity over a specified pressure range. Comparisons between hydraulic models are made based on scatterplots of pressure head and the vertical Darcian flux obtained using a finite-element numerical solution with both constitutive relations. For applicability of an equivalent Gardner function over a broad range of pressure heads, the crossover pressure must be maintained between the two parametric functions. The crossover pressure is defined as the pressure in which the hydraulic conductivity of the inclusion is equal to the background. It can be shown that a hybrid methodology of preserving the crossover pressure exactly and using the effective capillary drive will result in hydraulic parameters that are easily obtained and provide good agreement between the conductivity functions of the GR model to the VG model.  相似文献   

4.
The Beerkan method based on in situ single‐ring water infiltration experiments along with the relevant specific Beerkan estimation of soil transfer parameters (BEST) algorithm is attractive for simple soil hydraulic characterization. However, the BEST algorithm may lead to erroneous or null values for the saturated hydraulic conductivity and sorptivity especially when there are only few infiltration data points under the transient flow state, either for sandy soil or soils in wet conditions. This study developed an alternative algorithm for analysis of the Beerkan infiltration experiment referred to as BEST‐generalized likelihood uncertainty estimation (GLUE). The proposed method estimates the scale parameters of van Genuchten water retention and Brooks–Corey hydraulic conductivity functions through the GLUE methodology. The GLUE method is a Bayesian Monte Carlo parameter estimation technique that makes use of a likelihood function to measure the goodness‐of‐fit between modelled and observed data. The results showed that using a combination of three different likelihood measurements based on observed transient flow, steady‐state flow and experimental steady‐state infiltration rate made the BEST‐GLUE procedure capable of performing an efficient inverse analysis of Beerkan infiltration experiments. Therefore, it is more applicable for a wider range of soils with contrasting texture, structure, and initial and saturated water content. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

5.

The present study proposes a novel and simplified methodology to assess the seismic bearing capacity (SBC) of a shallow strip footing by incorporating strength non-linearity arising due to partial saturation of a soil matrix. Furthermore, developed methodology incorporates the modal response analysis of soil layers to assess SBC. A constant matric suction distribution profile has been considered throughout the depth of the soil. The Van Genuchten equation and corresponding fitting parameters have been considered to quantify matric suction in the analysis. SBC has been obtained for three different geomaterials; viz. sand, fly ash and clay, based on their predominant grain size and diverse soil water characteristics curve (SWCC) attributes. Variation of SBC with different modes of vibration and damping ratio are reported for ranges of matric suction pertinent to the geomaterials considered in the study. The relative significance of matric suction on SBC has been reported for suction values within the transition zone of each geomaterial. It is observed that the SBC of sand is drastically reduced, with matric suction reaching beyond the residual suction value. The SBC of fly ash remains constant beyond the residual suction value, whereas the SBC of clay shows an increasing trend toward the practical range of matric suction values.

  相似文献   

6.
An estimation of the volume of light nonaqueous phase liquids (LNAPL) is often required during site assessment, remedial design, or litigation. LNAPL volume can be estimated by a strictly empirical approach whereby core samples, distributed throughout the vertical and lateral extent of LNAPL, are analyzed for LNAPL content, and these data are then integrated to compute a volume. Alternatively, if the LNAPL has obtained vertical equilibrium, the thickness of LNAPL in monitoring wells can be used to calculate of LNAPL in monitoring wells can be used to calculate LNAPL volume at the well locations if appropriate soil and LNAPL properties can be estimated.
A method is described for estimating key soil and LNAPL properties by nonlinear regression of vertical profiles of LNAPL saturation. The methods is relatively fast, cost effective, and amenable to quantitative analysis of uncertainty. Optionally, the method allows statistical determination of best-fit values for the Van Genuchten capillary parameters (n, αoil-water and αoil-air), residual water saturation and ANAPL density. The sensitivity of the method was investigated by fitting field LNAPL saturation profiles and then determining the variation in misfit (mean square residual) as a function of parameter value for each parameter. Using field data from a sandy aquifer, the fitting statistics were found to be highly sensitive to LNAPL density, αoil-water and αoil-air moderately sensitive to the Van Genuchten n value, and weakly sensitive to residual water saturation. The regression analysis also provides information that can be used to estimate uncertainty in the estimated parameters, which can then be used to estimate uncertainty in calculated values of specific volume.  相似文献   

7.
Maximum-likelihood estimators properly represent measurement error, thus provide a statistically sound basis for evaluating the adequacy of a model fit and for finding the multivariate parameter confidence region. We demonstrate the advantages of using maximum-likelihood estimators rather than simple least-squares estimators for the problem of finding unsaturated hydraulic parameters. Inversion of outflow data given independent retention data can be treated by an extension to a Bayesian estimator. As an example, we apply the methodology to retention and transient unsaturated outflow observations, both obtained on the same medium sand sample. We found the van Genuchten expression to be adequate for the retention data, as the best fit was within measurement error. The Cramer–Rao confidence bound described the true parameter uncertainty approximately. The Mualem–van Genuchten expression was, however, inadequate for our outflow observations, suggesting that the parameters (, n) may not always be equivalent in describing both retention and unsaturated conductivity.  相似文献   

8.
Erick Carlier 《水文研究》2007,21(21):2845-2849
The infiltration rate in the unsaturated zone is analysed from a probabilistic point of view. It is shown that the empirical formulas of Horton and Kostiakov, without apparent physical basis, are explained in a probabilistic approach. Horton's formula reflects a Markovian process contrary to Kostiakov's formula. This approach made it possible to explain why Kostiakov's formula is more powerful than that of Horton. A new equation of infiltration is proposed. The three formulas were compared, for four types of soil, with the model of van Genuchten based on the Richards equation. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

9.
We present herein clear field evidence for the persistence of a coarse surface layer in a gravel‐bed river during flows capable of transporting all grain sizes present on the channel bed. Detailed field measurements of channel topography and bed surface grain size were made in a gravel‐bed reach of the Colorado River prior to a flood in 2003. Runoff produced during the 2003 snowmelt was far above average, resulting in a sustained period of high flow with a peak discharge of 27 m3/s (170% of normal peak flow); all available grain sizes within the study reach were mobilized in this period of time. During the 2003 peak flow, the river avulsed immediately upstream of the study reach, thereby abandoning approximately one half kilometer of the former channel. The abandonment was rapid (probably within a few hours), leaving the bed texture essentially frozen in place at the peak of the flood. All locations sampled prior to the flood were resampled following the stream abandonment. In response to the high flow, the surface median grain size (D50s) coarsened slightly in the outer part of the bend while remaining nearly constant along the inner part of the bend, resulting in an overall increase from 18 to 21 mm for the study reach. Thus, the coarse bed surface texture persisted despite shear stresses throughout the bend that were well above the critical entrainment value. This may be explained because the response of the bed texture to increases in flow strength depends primarily upon the continued availability of the various grain size percentiles in the supply, which in this case was essentially unlimited for all sizes present in the channel. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

10.
  相似文献   

11.
Geological heterogeneities affect the dynamics of carbon dioxide (CO2) plumes in subsurface environments in important ways. Previously, we showed how the dynamics of CO2 plumes are influenced by the multiscaled sedimentary architecture in deep brine fluvial-type reservoirs. The results confirm that representing small-scale features and the corresponding heterogeneity in saturation functions, along with hysteresis in saturation functions, are all critical to understanding capillary trapping processes. Here, we show that when heterogeneity and hysteresis are represented, the two conventional approaches for defining saturation functions, Brooks-Corey and van Genuchten, represent fundamentally different physical systems. The Brooks-Corey approach represents heterogeneity in entry pressures, and leads to trapping by capillary pinning. The van Genuchten approach represents a network of pores transporting the nonwetting fluid, across rock types, with negligible capillary entry pressure, and leads to capillary retardation. These differences significantly affect the large-scale characteristics of CO2 plumes (i.e., their mass, shape, and position).  相似文献   

12.
Estimates of soil hydraulic properties using pedotransfer functions (PTF) are useful in many studies such as hydrochemical modelling and soil mapping. The objective of this study was to calibrate and test parametric PTFs that predict soil water retention and unsaturated hydraulic conductivity parameters. The PTFs are based on neural networks and the Bootstrap method using different sets of predictors and predict the van Genuchten/Mualem parameters. A Danish soil data set (152 horizons) dominated by sandy and sandy loamy soils was used in the development of PTFs to predict the Mualem hydraulic conductivity parameters. A larger data set (1618 horizons) with a broader textural range was used in the development of PTFs to predict the van Genuchten parameters. The PTFs using either three or seven textural classes combined with soil organic mater and bulk density gave the most reliable predictions of the hydraulic properties of the studied soils. We found that introducing measured water content as a predictor generally gave lower errors for water retention predictions and higher errors for conductivity predictions. The best of the developed PTFs for predicting hydraulic conductivity was tested against PTFs from the literature using a subdata set of the data used in the calibration. The test showed that the developed PTFs gave better predictions (lower errors) than the PTFs from the literature. This is not surprising since the developed PTFs are based mainly on hydraulic conductivity data near saturation and sandier soils than the PTFs from the literature. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

13.
Abstract. A simple closed-form expression relating saturated hydraulic conductivity to the van Genuchten capillary retention model parameters is derived. Application of this equation to an experimental data set shows reasonable agreement between measured and predicted saturated conductivity values. The proposed equation provides a consistent theoretical basis for estimating both saturated and unsaturated hydraulic conductivity from statistical pore structure models.  相似文献   

14.
Needle ice was grown in the laboratory to determine the limits of soil moisture and texture on the growth of needle ice. As the percentage of fines in the soil increases the lower limit of soil moisture required for growth decreases and may be defined by a linear boundary. The efficiency of conversion of soil water to needle ice in soil samples used, peaked at about 16 per cent soil fines. Contrary to expectations the finer textured soils produced significantly thicker ice crystals. Quantity of soil lifted by needle ice was related to ice height, the ice mass to atmospheric void ratio, soil texture and moisture and soil surface roughness. However, no definite, quantitative relationships could be determined. A test to determine if needle ice is selective in the grain sizes of particles lifted suggested no selectivity.  相似文献   

15.
The selective radius shift model was used to relate changes in mineral volume due to precipitation/dissolution reactions to changes in hydraulic properties affecting flow in porous media. The model accounts for (i) precipitation/dissolution taking place only in the water-filled part of the pore space and further that (ii) the amount of mineral precipitation/dissolution within a pore depends on the local pore volume. The pore bundle concept was used to connect pore-scale changes to macroscopic soil hydraulic properties. Precipitation/dissolution induces changes in the pore radii of water-filled pores and, consequently, in the effective porosity. In a time step of the numerical model, mineral reactions lead to a discontinuous pore-size distribution because only the water-filled pores are affected. The pore-size distribution is converted back to a soil moisture characteristic function to which a new water retention curve is fitted under physically plausible constraints. The model equations were derived for the commonly used van Genuchten/Mualem hydraulic properties. Together with a mixed-form solution of Richards’ equation for aqueous phase flow, the model was implemented into the geochemical modelling framework PHREEQC, thereby making available PHREEQC’s comprehensive geochemical reactions. Example applications include kinetic halite dissolution and calcite precipitation as a consequence of cation exchange. These applications showed marked changes in the soil’s hydraulic properties due to mineral precipitation/dissolution and the dependency of these changes on water contents. The simulations also revealed the strong influence of the degree of saturation on the development of the saturated hydraulic conductivity through its quadratic dependency on the van Genuchten parameter α. Furthermore, it was shown that the unsaturated hydraulic conductivity at fixed reduced water content can even increase during precipitation due to changes in the pore-size distribution.  相似文献   

16.
This study is an extension of the stochastic analysis of transient two-phase flow in randomly heterogeneous porous media (Chen et al. in Water Resour Res 42:W03425, 2006), by incorporating direct measurements of the random soil properties. The log-transformed intrinsic permeability, soil pore size distribution parameter, and van Genuchten fitting parameter are treated as stochastic variables that are normally distributed with a separable exponential covariance model. These three random variables conditioned on given measurements are decomposed via Karhunen–Loève decomposition. Combined with the conditional eigenvalues and eigenfunctions of random variables, we conduct a series of numerical simulations using stochastic transient water–oil flow model (Chen et al. in Water Resour Res 42:W03425, 2006) based on the KLME approach to investigate how the number and location of measurement points, different random soil properties, as well as the correlation length of the random soil properties, affect the stochastic behavior of water and oil flow in heterogeneous porous media.  相似文献   

17.
易朝路 《湖泊科学》2001,13(4):296-303
在江汉平原濒临长江北岸的洪湖及其周围的土壤中,分别取得沉积物柱状钻孔样品和土壤剖面样品,用偏光显微镜和电子显微镜观察了2.5kaBP以来形成的三种湖相粘土的微结构特征,并以此解释的其成因和沉积环境。0.9-2.5kaBP期间形成的青色粘土的典型微结构有：凝胶结构、细颗粒粒径和小孔隙、矿物颗粒低圆度和淡水中心冈硅藻等生物框架结构,主要是有机质胶体与粘土胶体相互作用形成的,此期间,河流带入湖泊的泥沙少,洪湖拥有一个开阔、稳定、浮游生物较多的淡水湖泊环境。0.45-1kaBP期间形成的黑色粘土的主要微结构类型有;凝胶结构、大的圆孔隙和植物纤维等生物框架结构,呈双峰分布的孔隙,它主要是由于维管束植物残体大量积累形成的,当时洪湖湖水变浅,大量挺水植物生长,湖泊已经沼泽化了。0.45kaBP以来形成的灰色粘土的微结构类型主要有：絮凝结构、较大的微孔隙和颗粒粒径、矿物颗粒较高的圆度等,物源主要来自河流带来的泥沙,当时河湖相通,河流入湖水沙量大。  相似文献   

18.
Understanding the response of partially saturated earth structures under various static and dynamic loads is important for the design and construction of economical and safe geotechnical engineering structures. In this study, the numerical approach is used to understand the dynamics of partially saturated soils. The mathematical equations governing the dynamics of partially saturated soils are derived based on the theory of mixtures and implemented within a finite element framework. The stress–strain behavior of the soil is represented by an elasto-plastic constitutive model for unsaturated soil based on bounding surface concept and the moisture-suction behavior is modeled using van Genuchten model. Fully coupled finite element simulations are performed to study the response of partially saturated soil embankment under earthquake loading and validated with centrifuge test results available in the literature. The predicted displacement responses are in good agreement with the measured responses. The pore water pressure, pore air pressure, matric suction, the degree of saturation in various elements and the response of the embankment under different initial moisture content are also discussed.  相似文献   

19.
Often the soil hydraulic parameters are obtained by the inversion of measured data (e.g. soil moisture, pressure head, and cumulative infiltration, etc.). However, the inverse problem in unsaturated zone is ill‐posed due to various reasons, and hence the parameters become non‐unique. The presence of multiple soil layers brings the additional complexities in the inverse modelling. The generalized likelihood uncertainty estimate (GLUE) is a useful approach to estimate the parameters and their uncertainty when dealing with soil moisture dynamics which is a highly non‐linear problem. Because the estimated parameters depend on the modelling scale, inverse modelling carried out on laboratory data and field data may provide independent estimates. The objective of this paper is to compare the parameters and their uncertainty estimated through experiments in the laboratory and in the field and to assess which of the soil hydraulic parameters are independent of the experiment. The first two layers in the field site are characterized by Loamy sand and Loamy. The mean soil moisture and pressure head at three depths are measured with an interval of half hour for a period of 1 week using the evaporation method for the laboratory experiment, whereas soil moisture at three different depths (60, 110, and 200 cm) is measured with an interval of 1 h for 2 years for the field experiment. A one‐dimensional soil moisture model on the basis of the finite difference method was used. The calibration and validation are approximately for 1 year each. The model performance was found to be good with root mean square error (RMSE) varying from 2 to 4 cm3 cm?3. It is found from the two experiments that mean and uncertainty in the saturated soil moisture (θs) and shape parameter (n) of van Genuchten equations are similar for both the soil types. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

20.
Recognizing that soil–structure interaction affects appreciably the earthquake response of highway overcrossings, this paper compares approximate analytical solutions and finite element results to conclude on a simple procedure that allows for the estimation of the kinematic response functions and dynamic stiffnesses of approach embankments. It is shown that the shear‐wedge model yields realistic estimates for the amplification functions of typical embankments and reveals the appropriate levels of dynamic strains which are subsequently used to estimate the stiffness and damping coefficients of embankments. The shear‐wedge model is extended to a two‐dimensional model in order to calculate the transverse static stiffness of an approach embankment loaded at one end. The formulation leads to a sound closed‐form expression for the critical length, Lc, that is the ratio of the transverse static stiffness of an approach embankment and the transverse static stiffness of a unit‐width wedge. It is shown through two case studies that the transverse dynamic stiffness (‘spring’ and ‘dashpot’) of the approach embankment can be estimated with confidence by multiplying the dynamic stiffness of the unit‐width wedge with the critical length, Lc. The paper concludes that the values obtained for the transverse kinematic response function and dynamic stiffness can also be used with confidence to represent the longitudinal kinematic response function and dynamic stiffness, respectively. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

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

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