Time‐scale analysis in unsaturated porous media under external wave loads

Extreme waves caused by tsunamis and storm surges can lead to soil failures in the near‐shore region, which may have severe impact on coastal environments and communities. Multiphase flows in deformable porous media involve several coupled processes and multiple time scales, which are challenging for numerical simulations. The objective of this study is to investigate the roles of the various processes and their interactions in multiphase flows in unsaturated soils under external wave loading, via theoretical time‐scale analysis and numerical simulations. A coupled geomechanics–multiphase flow model based on conservation laws is used. Theoretical analysis based on coupled and decoupled models demonstrates that transient and steady‐state responses are governed by pore pressure diffusion and saturation front propagation, respectively, and that the two processes are essentially decoupled. Numerical simulations suggest that the compressibility of the pore fluids and the deformation of the soil skeleton are important when the transient responses of the media are of concern, while the steady‐state responses are not sensitive to these factors. The responses obtained from the fully coupled numerical simulations are explained by a simplified time‐scale analysis based on coupled and decoupled models. Copyright © 2010 John Wiley & Sons, Ltd.     

Two formulations for dynamic response of a cylindrical cavity in cross‐anisotropic porous media

Two formulations for calculating dynamic response of a cylindrical cavity in cross‐anisotropic porous media based on complex functions theory are presented. The basis of the method is the solution of Biot's consolidation equations in the complex plane. Employing two groups of potential functions for solid skeleton and pore fluid (each group includes three functions), the u–w formulation of Biot's equations are solved. Difference of these two solutions refers to use of two various potential functions. Equations for calculating stress, displacement and pore pressure fields of the medium are mentioned based on each two formulations. Copyright © 2009 John Wiley & Sons, Ltd.     

Two methods for pore network of porous media

Two approaches of generating pore networks of porous media are presented to capture the pore fabric. The first methodology extracted pore structure from a computer simulated packing of spheres. The modified Delaunay tessellation was used to describe the porous media, and modified Nelder–Mead method in conjunction with three pore‐merging algorithms was used to generate the pore size and coordination number distributions of the randomly packed spheres. The Biconical Abscissa Asymmetric CONcentric bond was used to describe the connection between two adjacent voids. This algorithm was validated by predicting pore structure of a cubic array of spheres of equal radius with known pore sizes, throat sizes and coordination number distributions. The predicted distributions of pore structure agreed well with the measured. Then, the algorithm was used to predict pore structure and permeability of randomly packed spherical particles, and predicted permeability values were compared with published experimental data. The results showed that the predicted permeability values were in good agreement with those measured, confirming the proposed algorithm can capture the main flow paths of packed beds. The second methodology generated an equivalent pore network of porous media, of which the centers of voids were located in a regular lattice with constant pore center distance. However, this network allowed for matching both main geometrical and topological characteristics of the porous media. A comparison of the two approaches suggested that the second approach can also be used as a predictive tool to quantitatively study the microscopic properties of flow through porous media. Copyright © 2012 John Wiley & Sons, Ltd.     

Variability of volume strain in bounded heterogeneous media

The deformation of the solid matrix affects the fluid pore pressure and flow by altering the pore volume. Such interaction in turn affects the storage of groundwater in the void space. Obviously, this subject is of interest in groundwater hydrology. This paper describes an investigation of the effect of aquifer heterogeneity on the variability of the fluid pressure head and solid's volume strain, where the assumption of a constant vertical total stress leads to a relatively simple relationship between changes in solid's volume strain and fluid pressure head. To solve the problem analytically, focus is placed on the one‐dimensional models. It is found from our closed‐form solutions that the variance and correlation length of the log hydraulic conductivity are important in increasing the variability of pressure head and solid's volume strain. It is hoped that our findings will provide a basic framework for understanding and quantifying field‐scale volume strain processes and be useful in stimulating further research in this area. Copyright © 2011 John Wiley & Sons, Ltd.     

河西地区重力动态变化的二维离散小波多尺度分解

采用二维离散小波的多尺度分解技术，对1994～1995年和2002～2003年的河西地区重力动态变化进行分解，提取该地区不同层次的重力动态变化异常，结果反映了不同深度在一定尺度意义上的介质密度的动态变化分布。     

Extension of the Griffith's fracture criteria to saturated clays

Inglis [1] has solved the problem of distribution of stress in an elastic plate around an elliptical hole. His works clarify the role of cracks in the failure of an elastic material. However, his solution cannot be applied to saturated clay because he considers only total stresses, while, in saturated clay, the criterion of rupture should be expressed in terms of effective and not total stresses. The solution of Atkinson and Craster [2] using Biot's poroelasticity theory, shows that there is no high pore pressure in the vicinity of the crack tips for saturated clay. The major difference between this approach and the Biot's theory of is that, in saturated clay, strain is a function of the variation of the effective stress [3], while, in poroelastic media, strain is only a function of the variation of the total stress [4, Equation 2.2]. Also in their solution there is continuity between the pore fluid and the inner fluid in the crack. Their solution is valid for poroelastic media involving a movement of the pore fluid. In our solution there is no movement of the pore fluid (Undrained condition). In this paper we have solved the same problem as Inglis [1], but for the particular case of saturated clay obeying elastic law. By solving this problem we obtained the expressions for pore pressure, effective stress, total stress and displacements. The results show that not only the total stress but also the pore pressure and the effective stress are also high in the vicinity of the crack tips. A new failure criterion, based on Griffith's strain energy principle [5] and maximum tensile stress [6], valid for saturated clay is developed in this paper. Copyright © 2003 John Wiley & Sons, Ltd.     

Effect of property variation and modelling on convection in a fluid overlying a porous layer

An accurate computational analysis is presented for the onset of thermal convection in a two‐layer system which is comprised of a saturated layer of porous material described by Darcy's law, over which lies a layer of the same saturating fluid. The two‐layer system is heated from below and the upper (fluid) surface is allowed to be fixed or stress free. The onset of convection may have a bi‐modal nature in which convection may be dominated by the porous medium or by the fluid depending on the depths of the relative layers, but this is strongly controlled by material parameters. The effect of variation of relevant fluid and porous material properties is investigated in detail, as is the effect of the interface boundary condition between the fluid and the porous medium. Copyright © 2001 John Wiley & Sons, Ltd.