Permeability and P-wave velocity change in granitic rocks under freeze–thaw cycles

An extensive experimental investigation of microstructural changes in granites under freeze–thaw cycles using permeability and P-wave velocity measurements is described. Two types of natural granite rocks are considered and tested under dry and saturated conditions. The specimens were subjected to 200 heating–cooling cycles (??20°C/?+?20°C); each cycle had a duration of 24 h. The results indicate that the ageing process decreases the permeability and the P-wave velocity. The reduction in P-wave velocity is likely to be due to microcrack development (material damage). In dry samples, the microcracks result from the repeated differential contraction–dilatation of the mineral components. In water-saturated samples, there is an additional effect of freezing and thawing of water in the porous network. The decrease in permeability in the dry samples is due to partial closure of existing microcracks. In water-saturated samples, there was no increase in the permeability. A physically acceptable explanation is that new microcracks are not necessarily connected with those that already exist. Therefore the physicochemical process resulting from water–rock interactions also affects the permeability. This phenomenon reduces fluid flow in the material.     

An expression for the permeability of anisotropic granular media

There are many expressions proposed for the permeability of isotropic media based on flow channel and pore size distribution concepts, but there are no such expressions for anisotropic media. In this paper an expression for the permeability of an anisotropic medium is proposed, which has been verified in the laboratory. The mechanism behind fluid flow through soil was investigated using microscopic computer simulations to propose an expression for macroscopic permeability. The soil was assumed to be a spatially periodic porous medium, and the Navier-Stokes equation was solved using the FEM with appropriate boundary conditions for several different arrangements of the porous medium. The basic variables influencing flow through soil at the microscopic level were identified as specific surface area, void ratio, particle shape, material heterogeneity and the arrangement of particles in a porous medium. A sensitivity analysis was carried out to obtain an expression for the permeability in terms of the above variables. The corresponding macroscopic variables for the above microscopic variables are average specific surface area, average void ratio, anisotropy, tortuosity due to material heterogeneity, and the arrangement of particles respectively. An expression for the directional permeability is proposed in terms of these variables for the most common occurrence of particles in a porous medium. For the verification of the proposed equation, the permeability values of a fine-grained sand were measured at different void ratios and were compared with those predicted by the proposed equation. The results show that the predicted permeability values from the proposed equation are very close to the measured values.     

Micromechanical Modeling of Anisotropic Damage-Induced Permeability Variation in Crystalline Rocks

This paper presents a study on the initiation and progress of anisotropic damage and its impact on the permeability variation of crystalline rocks of low porosity. This work was based on an existing micromechanical model considering the frictional sliding and dilatancy behaviors of microcracks and the recovery of degraded stiffness when the microcracks are closed. By virtue of an analytical ellipsoidal inclusion solution, lower bound estimates were formulated through a rigorous homogenization procedure for the damage-induced effective permeability of the microcracks-matrix system, and their predictive limitations were discussed with superconducting penny-shaped microcracks, in which the greatest lower bounds were obtained for each homogenization scheme. On this basis, an empirical upper bound estimation model was suggested to account for the influences of anisotropic damage growth, connectivity, frictional sliding, dilatancy, and normal stiffness recovery of closed microcracks, as well as tensile stress-induced microcrack opening on the permeability variation, with a small number of material parameters. The developed model was calibrated and validated by a series of existing laboratory triaxial compression tests with permeability measurements on crystalline rocks, and applied for characterizing the excavation-induced damage zone and permeability variation in the surrounding granitic rock of the TSX tunnel at the Atomic Energy of Canada Limited’s (AECL) Underground Research Laboratory (URL) in Canada, with an acceptable agreement between the predicted and measured data.     

Modelling of cement suspension flow in granular porous media

A theoretical model of cement suspensions flow in granular porous media considering particle filtration is presented in this paper. Two phenomenological laws have been retained for the filtration rate and the intrinsic permeability evolution. A linear evolution with respect to the volume fraction of cement in the grout has been retained for the filtration rate. The intrinsic permeability of the porous medium is looked for in the form of a hyperbolic function of the porosity change. The model depends on two phenomenological parameters only. The equations of this model are solved analytically in the one‐dimensional case. Besides, a numerical resolution based on the finite element method is also presented. It could be implemented easily in situations where no analytical solution is available. Finally, the predictions of the model are compared to the results of a grout injection test on a long column of sand. Copyright © 2005 John Wiley & Sons, Ltd.     

Assessment of acceleration modelling for fluid‐filled porous media subjected to dynamic loading

The purpose of this paper is to examine the importance of different possible simplifying approximations when performing numerical simulations of fluid‐filled porous media subjected to dynamic loading. In particular, the relative importance of the various acceleration terms for both the solid and the fluid, especially the convective contribution, is assessed. The porous medium is modelled as a binary mixture of a solid phase, in the sense of a porous skeleton, and a fluid phase that represents both liquid and air in the pores. The solid particles are assumed to be intrinsically incompressible, whereas the fluid is assigned a finite intrinsic compressibility. Finite element (FE) simulations are carried out while assuming material properties and loading conditions representative for a road structure. The results show that, for the range of the material data used in the simulations, omitting the relative acceleration gives differences in the solution of the seepage velocity field, whereas omitting only the convective term does not lead to significant differences. Copyright © 2007 John Wiley & Sons, Ltd.     

Coupling between anisotropic damage and permeability variation in brittle rocks

In this paper, a coupled constitutive model is proposed for anisotropic damage and permeability variation in brittle rocks under deviatoric compressive stresses. The formulation of the model is based on experimental evidences and main physical mechanisms involved in the scale of microcracks are taken into account. The proposed model is expressed in the macroscopic framework and can be easily implemented for engineering application. The macroscopic free enthalpy of cracked solid is first determined by approximating crack distribution by a second‐order damage tensor. The effective elastic properties of damaged material are then derived from the free enthalpy function. The damage evolution is related to the crack growth in multiple orientations. A pragmatic approach inspired from fracture mechanics is used for the formulation of the crack propagation criterion. Compressive stress induced crack opening is taken into account and leads to macroscopic volumetric dilatancy and permeability variation. The overall permeability tensor of cracked material is determined using a micro–macro averaging procedure. Darcy's law is used for fluid flow at the macroscopic scale whereas laminar flow is assumed at the microcrack scale. Hydraulic connectivity of cracks increases with crack growth. The proposed model is applied to the Lac du Bonnet granite. Generally, good agreement is observed between numerical simulations and experimental data. Copyright © 2005 John Wiley & Sons, Ltd.     

Modified Markov Chain Monte Carlo Method for Dynamic Data Integration Using Streamline Approach

In this paper, the Markov Chain Monte Carlo (MCMC) approach is used for sampling of the permeability field conditioned on the dynamic data. The novelty of the approach consists of using an approximation of the dynamic data based on streamline computations. The simulations using the streamline approach allows us to obtain analytical approximations in the small neighborhood of the previously computed dynamic data. Using this approximation, we employ a two-stage MCMC approach. In the first stage, the approximation of the dynamic data is used to modify the instrumental proposal distribution. The obtained chain correctly samples from the posterior distribution; the modified Markov chain converges to a steady state corresponding to the posterior distribution. Moreover, this approximation increases the acceptance rate, and reduces the computational time required for MCMC sampling. Numerical results are presented.     

An analytical solution for the nonlinear distribution of effective and total stresses in vertical backfilled stopes

The increasing use of backfill in underground mines requires a proper evaluation of the stress state in and around the filled openings. This is, however, a relatively complex issue due, in part, to the large contrast in strength and stiffness between the backfill material and surrounding rock mass. In recent years, it has been shown that arching theory, based on limit equilibrium analysis, can be used to estimate the stress distribution in backfilled stopes. Nonetheless, many simplifications are involved in such analytical solutions and this affects the precision and significance of the calculated results. In this paper, a previously developed solution is enhanced by introducing the combined effects of non-uniform vertical stress distribution and positive pore water pressure. This leads to a more representative analytical solution, as demonstrated by successful comparisons with numerical simulations. The results indicate that the proposed solution can be used to estimate the effective and total stress state in submerged or partially submerged backfilled stopes with a simple geometry.     

Hydrogen penetration in water through porous medium: application to a radioactive waste storage site

Hydrogen penetration in water through porous medium was analyzed in the paper. A two-phase compositional model approach was considered. The first part of the work deals with the thermodynamic analysis of the hydrogen–water system. The thermodynamic model was calibrated using the experimental data of hydrogen solubility in water. The phase densities, viscosities and phase concentrations were presented in an analytical form. Moreover, the domain of validity of analytical laws—such as Henry’s, Raoult’s and Kelvin’s laws—for the estimation of phase properties was presented for the analyzed system. The second part deals with two-phase hydrodynamic behaviors. An analytical solution for the non-compressible flow was constructed. In general case, the influence of relative permeabilities on the flow regimes was analyzed numerically. The notion pseudo-saturation was introduced to define phase appearance. Actually, mobile gas created a time displaced front relatively slower than mobile gas flow. Diffusion becomes really important for low mobile gas case as the penetration accelerates for the large range of saturation. In contrast, the mass exchange phenomena have a small influence on the flow type. Thus, the regimes of hydrogen penetration in liquid were shown really sensitive to the relative permeability form.     

Effective medium methods and a computational approach for estimating geomaterial properties of porous materials with randomly oriented ellipsoidal pores

The role of effective medium approaches and the differential scheme in estimating the overall elastic moduli of a porous medium with randomly oriented pores is examined. The analytical estimates of the elastic moduli are compared with approximations available in the literature that are valid for small pore aspect ratios. Accuracy of the analytical estimates is further established by performing finite element simulations. Finite element estimates are obtained for a model of a porous medium containing three families of spheroidal pores arranged in a mutually orthogonal configuration. This model is regarded as a close approximation to a porous medium with randomly oriented pores. Also, experimental data available for several sandstones and granites was used to analytically examine the influence of the aspect ratio of pores on overall properties. This information is also used to provide an estimate for the permeability of the porous medium.     

Fully coupled poromechanical back analysis of the pulse test by inverse method

In this paper, we investigate the pulse test, which is usually considered as efficient for measuring the permeability of weakly permeable porous media. The pulse is first analyzed and we show that it is a fully poromechanical coupled problem. Owing to those couplings, the problem is 2D‐axisymmetrical, rather than 1D as is usually considered to be the case. As a consequence, the 1D solutions, for example under constant mean stress hypothesis, although giving good approximates of permeability and storage coefficient, are not rigorous and an enhanced back analysis of the test requires 2D solutions. Therefore, no analytical solution exists, and the problem has to be solved using 2D‐axisymmetrical numerical modelings of the pulse test. The finite element method is considered in this paper. We then proceed to formulate the pulse test back analysis as a parameter identification problem, and we focus on intrinsic permeability, Biot coefficient, drained Young's modulus and reservoir compressibility levels. This parameter identification problem is solved by an inverse method consisting of the minimization of a cost‐functional, through a gradient‐based algorithm. This new method of interpretation of the pulse test is finally applied to laboratory tests on Meuse/Haute–Marne argillite and is shown to give encouraging results. Copyright © 2010 John Wiley & Sons, Ltd.     

Efficient block preconditioners for the coupled equations of pressure and deformation in highly discontinuous media

Large‐scale simulations of flow in deformable porous media require efficient iterative methods for solving the involved systems of linear algebraic equations. Construction of efficient iterative methods is particularly challenging in problems with large jumps in material properties, which is often the case in geological applications, such as basin evolution at regional scales. The success of iterative methods for this type of problems depends strongly on finding effective preconditioners. This paper investigates how the block‐structured matrix system arising from single‐phase flow in elastic porous media should be preconditioned, in particular for highly discontinuous permeability and significant jumps in elastic properties. The most promising preconditioner combines algebraic multigrid with a Schur complement‐based exact block decomposition. The paper compares numerous block preconditioners with the aim of providing guidelines on how to formulate efficient preconditioners. Copyright © 2010 John Wiley & Sons, Ltd.     

Stochastic representation for anisotropic permeability tensor random fields

In this paper, we introduce a novel stochastic model for the permeability tensor associated with stationary random porous media. In the light of recent works on mesoscale modeling of permeability, we first discuss the physical interpretation of the permeability tensor randomness. Subsequently, we propose a nonparametric prior probabilistic model for non‐Gaussian permeability tensor random fields, making use of the information theory and a maximum entropy procedure, and provide a physical interpretation of the model parameters. Finally, we demonstrate the capability of the considered class of random fields to generate higher levels of statistical fluctuations for selected stochastic principal permeabilities. This unique flexibility offered by the parameterization of the model opens up many new possibilities for both forward simulations (e.g. for uncertainty propagation in predictive simulations) and stochastic inverse problem solving. Copyright © 2011 John Wiley & Sons, Ltd.     

Analytical solution of the harmonic waves diffraction by a cylindrical lined cavity in poroelastic saturated medium

This paper presents a model for the analysis of plane waves diffraction at a cavity in an infinite homogeneous poroelastic saturated medium, lined by a lining composed of four equal segments. An elastic boundary layer is placed between the cavity lining and the infinite porous medium. The boundary layer is simulated by ‘elastic boundary conditions’ in which the bulk matrix stress is proportional to the relative displacement between the lining and the surrounding medium matrix boundary. In addition, fluid impermeability through the intermediate layer is assumed. For the frequencies, that differ from the pseudoresonanse frequencies, the problem was reduced to the problem of an ideal elastic medium. A closed‐form analytical solution of the problem was obtained using Fourier–Bessel series, the convergence of which was proven. It was shown that the number of series terms required to obtain a desired level of accuracy can be determined in advance. The influence of the medium porosity on the medium dynamic stress concentration was studied. Copyright © 2006 John Wiley & Sons, Ltd.     

Adaptive Multiscale Permeability Estimation

With multiscale permeability estimation one does not select parameterization prior to the estimation. Instead, one performs a hierarchical search for the right parameterization while solving a sequence of estimation problems with an increasing parameterization dimension. In some previous works on the subject, the same refinement is applied all over the porous medium. This may lead to over-parameterization, and subsequently, to unrealistic permeability estimates and excessive computational work. With adaptive multiscale permeability estimation, the new parameterization at an arbitrary stage in the estimation sequence is such that new degrees of freedom are not necessarily introduced all over the porous medium. The aim is to introduce new degrees of freedom only where it is warranted by the data. In this paper, we introduce a novel adaptive multiscale estimation. The approach is used to estimate absolute permeability from two-phase pressure data in several numerical examples.     

多孔介质中非水相流体运移的数值模拟

针对多孔介质中水、气和非水相流体(NAPLs)的多相流动特点,建立了非水相流体(NAPLs)污染物迁移模型,分析了非水相流体在土壤非饱和区和地下水系统中的运移规律。通过有限元数值解对轻非水相流体和重非水相流体在土壤系统中的迁移过程进行模拟,得到了污染物的时空分布特征和污染范围。计算结果表明,数值模拟方法能够合理地描述非水相流体的运移过程和污染特征。土体渗透性和污染物残余饱和度是其重要影响因素。     

Responses of chemically active and naturally fractured shale under time‐dependent mechanical loading and ionic solution exposure

In this paper, the analytical dual‐porosity dual‐permeability poromechanics solution for saturated cylinders is extended to account for electrokinetic effects and material transverse isotropy, which simulate the responses of chemically active naturally fractured shale under time‐dependent mechanical loading and ionic solution exposure. The solution addresses the stresses, fracture pore pressure, matrix pore pressure, fluid fluxes, ion concentration evolution, and displacements due to the applied stress, pore pressure, and solute concentration difference between the sample and the circulation fluid. The presented solution will not only help validate numerical simulations but also assist in calibrating and interpreting laboratory results on dual‐porosity dual‐permeability shale. It is recommended that the analytical solutions of radial and axial displacements be used to match the corresponding laboratory‐recorded data to determine shale dual permeability and chemo‐electrical parameters including membrane coefficient, ions diffusion coefficients, and electro‐osmotic permeability.     

Modelling of steady‐state fluid flow in 3D fractured isotropic porous media: application to effective permeability calculation

In this paper, 3D steady‐state fluid flow in a porous medium with a large number of intersecting fractures is derived numerically by using collocation method. Fluid flow in the matrix and fractures is described by Darcy's law and Poiseuille's law, respectively. The recent theoretical development presented a general potential solution to model the steady‐state flow in fractured porous media under a far‐field condition. This solution is a hypersingular integral equation with pressure field in the fracture surfaces as the main unknown. The numerical procedure can resolve the problem for any form of fractures and also takes into account the interactions and the intersection between fractures. Once the pressure field and then the flux field in fractures have been determined, the pressure field in the porous matrix is computed completely. The basic problem of a single fracture is investigated, and a semi‐analytical solution is presented. Using the solution obtained for a single fracture, Mori‐Tanaka and self‐consistent schemes are employed for upscaling the effective permeability of 3D fractured porous media. Copyright © 2012 John Wiley & Sons, Ltd.     

Simulation of damage–permeability coupling for mortar under dynamic loads

The results reported in this paper deal with the simulation of damage in cohesive geomaterials such as rocks or concrete subjected to dynamic loads. The practical objective is to stimulate the production of tight gas reservoirs with a technique that is an alternative to hydraulic fracturing. The principle is that when subjected to dynamic loads, cohesive materials such as concrete, rocks or ceramics exhibit distributed micro‐cracking as opposed to localised cracking observed under static loads. Hence, a low permeability rock containing gas trapped into occluded pores can be fragmented with the help of dynamic loads, and gas can be extracted in a much more efficient way compared with hydraulic fracturing, where only large macro cracks are formed with very few connections between occluded pores. At the stage of laboratory development of this technique, compressive underwater shock waves have been used to increase the intrinsic permeability of concrete specimens. In a previous study, pressure waves generated by pulsed arc electrohydraulic discharges in water were used in order to induce micro‐cracking and an increase of average permeability of concrete hollow cylinders subjected to confinement stresses (equivalent to geostatic stresses). We discuss here a 3‐D anisotropic constitutive model aimed at describing the dynamic response of these specimens. It is based on rate‐dependent continuum damage constitutive relations. Crack closure effects and damage‐induced anisotropy are included in the model. The directional growth of damage is related to the directional growth of material intrinsic permeability. Numerical simulations of damage induced by shock waves show good agreement with the experiments for various confinement levels of the specimens. Copyright © 2013 John Wiley & Sons, Ltd.