A methodology for regionalizing 3‐D effective porosity at watershed scale in crystalline aquifers          下载免费PDF全文

Benoît Dewandel  Yvan Caballero  Jérôme Perrin  Alexandre Boisson  Fabrice Dazin  Sylvain Ferrant  Subash Chandra  Jean‐Christophe Maréchal 《水文研究》2017,31(12):2277-2295

An innovative approach for regionalizing the 3‐D effective porosity field is presented and applied to two large, overexploited, and deeply weathered crystalline aquifers located in southern India. The method derives from earlier work on regionalizing a 2‐D effective porosity field in that part of an aquifer where the water table fluctuates, which is now extended over the entire aquifer using a 3‐D approach. A method based on geological and geophysical surveys has also been developed for mapping the weathering profile layers (saprolite and fractured layers). The method for regionalizing 3‐D effective porosity combines water table fluctuation and groundwater budget techniques at various cell sizes with the use of satellite‐based data (for groundwater abstraction), the structure of the weathering profile, and geostatistical techniques. The approach is presented in detail for the Kudaliar watershed (983 km²) and tested on the 730 km² Anantapur watershed. At watershed scale, the effective porosity of the aquifer ranges from 0.5% to 2% in Kudaliar and between 0.3% and 1% in Anantapur, which agrees with earlier works. Results show that (a) depending on the geology and on the structure of the weathering profile, the vertical distribution of effective porosity can be very different and that the fractured layers in crystalline aquifers are not necessarily characterized by a rapid decrease in effective porosity and (b) that the lateral variations in effective porosity can be larger than the vertical ones. These variations suggest that within a same weathering profile, the density of open fractures and/or degree of weathering in the fractured zone may significantly vary from a place to another. The proposed method provides information on the spatial distribution of effective porosity that is of prime interest in terms of flux and contaminant transport in crystalline aquifers. Implications for mapping groundwater storage and scarcity are also discussed, which should help in improving groundwater resource management strategies.  相似文献   

The effective stress concept and the evaluation of changes in pore air pressure under jacketed isotropic compression tests for dry sand based on the 2‐phase porous theory          下载免费PDF全文

Sayuri Kimoto  Fusao Oka  Yasuhiro Morimoto 《国际地质力学数值与分析法杂志》2017,41(18):1894-1907

The effective stress concept for solid‐fluid 2‐phase media was revisited in this work. In particular, the effects of the compressibility of both the pore fluid and the soil particles were studied under 3 different conditions, i.e., undrained, drained, and unjacketed conditions based on a Biot‐type theory for 2‐phase porous media. It was confirmed that Terzaghi effective stress holds at the moment when soil grains are assumed to be incompressible and when the compressibility of the pore fluid is small enough compared to that of the soil skeleton. Then, isotropic compression tests for dry sand under undrained conditions were conducted within the triaxial apparatus in which the changes in the pore air pressure could be measured. The ratio of the increment in the cell pressure to the increment in the pore air pressure, m, corresponds to the inverse of the B value by Bishop and was obtained during the step loading of the cell pressure. In addition, the m values were evaluated by comparing them with theoretically obtained values based on the solid‐fluid 2‐phase mixture theory. The experimental m values were close to the theoretical values, as they were in the range of approximately 40 to 185, depending on the cell pressure. Finally, it was found that the soil material with a highly compressible pore fluid, such as air, must be analyzed with the multi‐phase porous mixture theory. However, Terzaghi effective stress is practically applicable when the compressibilities of both the soil particles and the pore fluid are small enough compared to that of the soil skeleton.  相似文献   

Elastic lateral dynamic impedance functions for a rigid cylindrical shell type foundation          下载免费PDF全文

R. He  R.Y.S. Pak  L.Z. Wang 《国际地质力学数值与分析法杂志》2017,41(4):508-526

Elastic lateral dynamic impedance functions are defined as the ratio of the lateral dynamic force/moment to the corresponding lateral displacement/rotation at the top ending of a foundation at very small strains. Elastic lateral dynamic impedance functions have a defining influence on the natural frequencies of offshore wind turbines supported on cylindrical shell type foundations, such as suction caissons, bucket foundations, and monopiles. This paper considers the coupled horizontal and rocking vibration of a cylindrical shell type foundation embedded in a fully saturated poroelastic seabed in contact with a seawater half‐space. The formulation of the coupled seawater–shell–seabed vibration problem is simplified by treating the shell as a rigid one. The rigid shell vibration problem is approached by the integral equation method using ring‐load Green's functions for a layered seawater‐seabed half‐space. By considering the boundary conditions at the shell–soil interface, the shell vibration problem is reduced to Fredholm integral equations. Through an analysis of the corresponding Cauchy singular equations, the intrinsic singular characteristics of the problem are rendered explicit. With the singularities incorporated into the solution representation, an effective numerical method involving Gauss–Chebyshev method is developed for the governing Fredholm equations. Selected numerical results for the dynamic contact load distributions, displacements of the shell, and lateral dynamic impedance functions are examined for different shell length–radius ratio, poroelastic materials, and frequencies of excitation. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

On the dependence of the saturated hydraulic conductivity upon the effective porosity through a power law model at different scales          下载免费PDF全文

Carmine Fallico  Samuele De Bartolo  Massimo Veltri  Gerardo Severino 《水文研究》2016,30(13):2366-2372

We study the scale dependence of the saturated hydraulic conductivity K_s through the effective porosity n_e by means of a newly developed power‐law model (PLM) which allows to use simultaneously measurements at different scales. The model is expressed as product between a single PLM (capturing the impact of the dominating scale) and a characteristic function κ^? accounting for the correction because of the other scale(s). The simple (closed form) expression of the κ^?‐function enables one to easily identify the scales which are relevant for K_s. The proposed model is then applied to a set of real data taken at the experimental site of Montalto Uffugo (Italy), and we show that in this case two (i.e. laboratory and field) scales appear to be the main ones. The implications toward an important application (solute transport) in Hydrology are finally discussed. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Inverse modelling using PS‐InSAR data for improved land subsidence simulation in Las Vegas Valley,Nevada          下载免费PDF全文

Meijing Zhang  Thomas J. Burbey 《水文研究》2016,30(24):4494-4516

Las Vegas Valley has had a long history of groundwater development and subsequent surface deformation. InSAR interferograms have revealed detailed and complex spatial patterns of subsidence in the Las Vegas Valley area that do not coincide with major pumping regions. This research represents the first effort to use high spatial and temporal resolution subsidence observations from InSAR and hydraulic head data to inversely calibrate transmissivities (T), elastic and inelastic skeletal storage coefficients (S_ke and S_kv) of the developed‐zone aquifer and conductance (CR) of the basin‐fill faults for the entire Las Vegas basin. The results indicate that the subsidence observations from InSAR are extremely beneficial for accurately quantifying hydraulic parameters, and the model calibration results are far more accurate than when using only groundwater levels as observations, and just a limited number of subsidence observations. The discrepancy between distributions of pumping and greatest levels of subsidence is found to be attributed to spatial variations in clay thickness. The Eglington fault separates thicker interbeds to the northwest from thinner interbeds to the southeast and the fault may act as a groundwater‐flow barrier and/or subsidence boundary, although the influence of the groundwater barrier to this area is found to be insignificant. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Displacement analytical prediction of shallow tunnel based on unified displacement function under slope boundary     

Fanchao Kong  Dechun Lu  Xiuli Du  Chenpeng Shen 《国际地质力学数值与分析法杂志》2019,43(1):183-211

For slope condition of ground surface, the asymmetrical deformation about the vertical center line and the horizontal center line of the tunnel cross section can be formed. A unified displacement function expressed by the Fourier series is presented to express the asymmetrical deformation of the tunnel cross section. Five basic deformation modes corresponding to the expansion order 2 are a complete deformation mode to reflect deformation behaviors of the tunnel cross section under slope boundary. Such this complete displacement mode is implemented into the complex variable solution for analytically predicting tunneling-induced ground deformation under slope boundary. All of these analytical solutions are verified by good agreements of the comparison between the analytical solutions and finite element method results. A parameter study is carried out to investigate the influence of deformation modes of the tunnel cross section, geometrical conditions of the tunnel and the slope angle, and “Buoyancy effect” on the displacement field. Finally, the proposed method is consistent with measured data of the Hejie tunnel in China qualitatively. The presented solution can provide a simplified indication for evaluating the ground deformation under slope condition of ground surface.  相似文献   

A multidirectional semi‐analytical method for analysis of laterally loaded pile groups in multi‐layered elastic strata          下载免费PDF全文

F.S. Tehrani  M. Prezzi  R. Salgado 《国际地质力学数值与分析法杂志》2016,40(12):1730-1757

A semi‐analytical method for calculating the response of single piles and pile groups subjected to lateral loading is developed in this paper. Displacements anywhere in the soil domain are tied to the displacements of the piles through decay functions. The principle of virtual work and the calculus of variations are used to derive the governing differential equations that describe the response of the piles and soil. The eigenvalue method and the finite difference technique are used to solve the system of coupled differential equations for the piles and soil, respectively. The proposed method takes into account the soil surface displacement along and perpendicular to the loading direction and produces displacement fields that are very close to those produced by the finite element method but at lower computational effort. Compared with the previous method that considered only the soil displacement along the loading direction, accounting for the multi‐directional soil displacement field produces responses for the piles and soil that are closer to those approximated by the finite element method. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

A smoothed cap model for variably saturated soils and its robust numerical implementation          下载免费PDF全文

P. Gamnitzer  G. Hofstetter 《国际地质力学数值与分析法杂志》2015,39(12):1276-1303

This paper deals with the numerical implementation of a cap model for unsaturated soils. It provides a brief review of existing cap model approaches, based on which an improved model formulated in terms of generalised effective stress and matric suction is derived and described in detail. Although the proposed model is a multisurface plasticity model, it can efficiently be implemented using only single‐surface projections because of the smoothness of the model, which is obtained by construction. Numerical algorithms are provided for these single‐surface stress projections, using a single‐equation approach whenever possible. The robustness of the utilised single‐equation approaches is enhanced by proposing problem‐fitted start‐up procedures based on investigations of the nonlinear projection equations. A comparison of the model response with extensive material test data is used to validate the model and to demonstrate the robust application of the approach to silty sands and low to medium plasticity clays. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Predicting the dependency of a degree of saturation on void ratio and suction using effective stress principle for unsaturated soils     

David Mašín 《国际地质力学数值与分析法杂志》2010,34(1):73-90

The paper presents an approach to predicting variation of a degree of saturation in unsaturated soils with void ratio and suction. The approach is based on the effective stress principle for unsaturated soils and several underlying assumptions. It focuses on the main drying and wetting processes and does not incorporate the effects of hydraulic hysteresis. It leads to the dependency of water retention curve (WRC) on void ratio, which does not require any material parameters apart from the parameters specifying WRC for the reference void ratio. Its validity is demonstrated by comparing predictions with the experimental data on four different soils taken over from the literature. Good correlation between the measured and predicted behaviour indirectly supports applicability of the effective stress principle for unsaturated soils. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates     

J. Mertens  H. Madsen  M. Kristensen  D. Jacques  J. Feyen 《水文研究》2005,19(8):1611-1633

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献