Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

Subsurface flow measurements using passive flux meters in variably-saturated cold-regions landscapes

To date, passive flux meters have predominantly been applied in temperate environments for tracking the movement of contaminants in groundwater. This study applies these instruments to reduce uncertainty in (typically instantaneous) flux measurements made in a low-gradient, wetland dominated, discontinuous permafrost environment. This method supports improved estimation of unsaturated and over-winter subsurface flows which are very difficult to quantify using hydraulic gradient-based approaches. Improved subsurface flow estimates can play a key role in understanding the water budget of this landscape.     

Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.     

Geostatistical features of streambed vertical hydraulic conductivities in Frenchman Creek Watershed in Western Nebraska

This study evaluated the spatial variability of streambed vertical hydraulic conductivity (K_v) in different stream morphologies in the Frenchman Creek Watershed, Western Nebraska, using different variogram models. Streambed K_v values were determined in situ using permeameter tests at 10 sites in Frenchman, Stinking Water and Spring Creeks during the dry season at baseflow conditions. Measurements were taken both in straight and meandering stream channels during a 5 day period at similar flow conditions. Each test site comprised of at least three transects and each transect comprised of at least three K_v measurements. Linear, Gaussian, exponential and spherical variogram models were used with Kriging gridding method for the 10 sites. As a goodness-of-fit statistic for the variogram models, cross-validation results showed differences in the median absolute deviation and the standard deviation of the cross-validation residuals. Results show that using the geometric means of the 10 sites for gridding performs better than using either all the K_v values from the 93 permeameter tests or 10 K_v values from the middle transects and centre permeameters. Incorporating both the spatial variability and the uncertainty involved in the measurement at a reach segment can yield more accurate grid results that can be useful in calibrating K_v at watershed or sub-watershed scales in distributed hydrological models.     

Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Simultaneous transport of water and solutes under transient unsaturated flow conditions — A case study

The imbalance between incoming and outgoing salt causes salinization of soils and sub-soils that result in increasing the salinity of stream-flows and agriculture land. This salinization is a serious environmental hazard particularly in semi-arid and arid lands. In order to estimate the magnitude of the hazard posed by salinity, it is important to understand and identify the processes that control salt movement from the soil surface through the root zone to the ground water and stream flows. In the present study, Malaprabha sub-basin (up to dam site) has been selected which has two distinct climatic zones, sub-humid (upstream of Khanapur) and semi-arid region (downstream of Khanapur). In the upstream, both surface and ground waters are used for irrigation, whereas in the downstream mostly groundwater is used. Both soils and ground waters are more saline in downstream parts of the study area. In this study we characterized the soil salinity and groundwater quality in both areas. An attempt is also made to model the distribution of potassium concentration in the soil profile in response to varying irrigation conditions using the SWIM (Soil-Water Infiltration and Movement) model. Fair agreement was obtained between predicted and measured results indicating the applicability of the model.     

煤岩对压裂裂缝长期导流能力影响的实验研究

煤层中一般都伴有煤层气,通过压裂可以将煤层气释放出来,一方面能够获得清洁能源,另一方面可以提高煤田的安全开采。煤岩层自身特点决定了其在压裂中应用的实验评价方法及结论应具有特殊性,本文主要针对煤层气井的煤岩层压裂裂缝长期导流能力进行实验评价研究。实验表明,与砂岩地层不同,煤岩的硬度较小,压裂中支撑剂嵌入情况较严重;煤岩易破碎,碎屑颗粒充填到支撑剂中;同时支撑剂存在破碎及吸附伤害等,多种作用导致支撑裂缝的导流能力降低。实验中,以煤岩层闭合压力为一个重要的参数,实验对比分析了在不同闭合压力下,铺砂浓度、支撑剂粒径、支撑剂种类及组合的选择、支撑剂颗粒的破碎与嵌入、压裂液残渣等因素对煤岩层导流能力的影响。本文的实验结论为,嵌入对煤岩层支撑裂缝的导流能力伤害很大,加大铺砂浓度、采用低破碎率支撑剂、降低压裂液伤害能在一定程度上提高煤层裂缝的导流能力,但影响程度随闭合压力的增加而逐渐降低。本文所得出的结论对今后煤层气的开发及设计施工具有一定的指导意义。     

The relationship between hydrodynamic properties and weathering of soils derived from volcanic rocks,Galapagos Islands (Ecuador)

The aim of this interdisciplinary study is to examine a component of the hydrological cycle in Galapagos by characterizing soil properties. Nine soil profiles were sampled on two islands. Their physical and hydrodynamic properties were analyzed, along with their mineralogical composition. Two groups of soils were identified, with major differences between them. The first group consists of soils located in the highlands (>350 m a.s.l.), characterized by low hydraulic conductivity (<10⁻⁵ m s⁻¹) and low porosity (<25%). These soils are thick (several meters) and homogeneous without coarse components. Their clay fraction is considerable and dominated by gibbsite. The second group includes soils located in the low parts of the islands (<300 m a.s.l.). These soils are characterized by high hydraulic conductivity (>10⁻³ m s⁻¹) and high porosity (>35%). The structure of these soils is heterogeneous and includes coarse materials. The physical properties of the soils are in good agreement with the variations of the rainfall according to the elevation, which appears as the main factor controlling the soil development. The clayey alteration products constrain soils physical and hydrodynamic properties by reducing the porosity and consequently the permeability and also by increasing water retention.     

Numerical modeling of rainstorm-induced shallow landslides in saturated and unsaturated soils

For the assessment of shallow landslides triggered by rainfall, the physically based model coupling the infinite slope stability analysis with the hydrological modeling in nearly saturated soil has commonly been used due to its simplicity. However, in that model the rainfall infiltration in unsaturated soil could not be reliably simulated because a linear diffusion-type Richards’ equation rather than the complete Richards’ equation was used. In addition, the effect of matric suction on the shear strength of soil was not actually considered. Therefore, except the shallow landslide in saturated soil due to groundwater table rise, the shallow landslide induced by the loss in unsaturated shear strength due to the dissipation of matric suction could not be reliably assessed. In this study, a physically based model capable of assessing shallow landslides in variably saturated soils is developed by adopting the complete Richards’ equation with the effect of slope angle in the rainfall infiltration modeling and using the extended Mohr–Coulomb failure criterion to describe the unsaturated shear strength in the soil failure modeling. The influence of rainfall intensity and duration on shallow landslide is investigated using the developed model. The result shows that the rainfall intensity and duration seem to have similar influence on shallow landslides respectively triggered by the increase of positive pore water pressure in saturated soil and induced by the dissipation of matric suction in unsaturated soil. The rainfall duration threshold decreases with the increase in rainfall intensity, but remains constant for large rainfall intensity.