Hydro‐mechanical modeling of cohesive crack propagation in multiphase porous media using the extended finite element method

In this paper, a numerical model is developed for the fully coupled hydro‐mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non‐wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two‐phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid skeleton is accounted. The cohesive crack model is integrated into the numerical modeling by means of which the nonlinear fracture processes occurring along the fracture process zone are simulated. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three‐phase formulation. The other variables are incorporated into the model via the experimentally determined functions, which specify the relationship between the hydraulic properties of the fracturing porous medium, that is saturation, permeability and capillary pressure. The spatial discretization is implemented by employing the extended finite element method, and the time domain discretization is performed using the generalized Newmark scheme to derive the final system of fully coupled nonlinear equations of the hydro‐mechanical problem. It is illustrated that by allowing for the interaction between various processes, that is the solid skeleton deformation, the wetting and the non‐wetting pore fluid flow and the cohesive crack propagation, the effect of the presence of the geomechanical discontinuity can be completely captured. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical modeling of multiphase fluid flow in deforming porous media: A comparison between two- and three-phase models for seismic analysis of earth and rockfill dams

In this paper, a fully coupled numerical model is presented for the finite element analysis of the deforming porous medium interacting with the flow of two immiscible compressible wetting and non-wetting pore fluids. The governing equations involving coupled fluid flow and deformation processes in unsaturated soils are derived within the framework of the generalized Biot theory. The displacements of the solid phase, the pressure of the wetting phase and the capillary pressure are taken as the primary unknowns of the present formulation. The other variables are incorporated into the model using the experimentally determined functions that define the relationship between the hydraulic properties of the porous medium, i.e. saturation, relative permeability and capillary pressure. It is worth mentioning that the imposition of various boundary conditions is feasible notwithstanding the choice of the primary variables. The modified Pastor–Zienkiewicz generalized constitutive model is introduced into the mathematical formulation to simulate the mechanical behavior of the unsaturated soil. The accuracy of the proposed mathematical model for analyzing coupled fluid flows in porous media is verified by the resolution of several numerical examples for which previous solutions are known. Finally, the performance of the computational algorithm in modeling of large-scale porous media problems including the large elasto-plastic deformations is demonstrated through the fully coupled analysis of the failure of two earth and rockfill dams. Furthermore, the three-phase model is compared to its simplified one which simulates the unsaturated porous medium as a two-phase one with static air phase. The paper illustrates the shortcomings of the commonly used simplified approach in the context of seismic analysis of two earth and rockfill dams. It is shown that accounting the pore air as an independent phase significantly influences the unsaturated soil behavior.     

Irreducible water distribution in sandstone rock: Two phase flow simulations in CT-based pore network

A simple cutoff approach based on the capillary bundle model has become an industrial standard method to quickly obtain the irreducible water saturation from NMR T₂ distribution. However this approach is not always valid. To overcome the shortcoming of the capillary bundle model that ignores pore-pore connectivity, we have conducted a two-phase flow simulation in a CT-based pore network. The CT -based pore network is a representation of a real rock pore structure that is described by a binary X-ray tomographic data set. Simulation on such network mimics a process of the porous plate measurement. The generated capillary curve is quite reasonable. The oil accession distributions at different water saturations plotted as a function of pore size provide an insight for the immiscible displacement process in the real rock pore structure. Water is trapped not only in dead ends but also in the well-connected pores due to a pore level by-pass mechanism. At the capillary end point pressure, a plot of the trapped water distribution as a function of pore size has the same lognormal distribution as the pore size distribution, which is much different from what the simple capillary bundle model suggests.     

碳酸盐岩储层参数对微观渗流的影响

针对碳酸盐岩油藏孔隙大小的双峰分布特征,首先利用计算机模拟建立了描述不同孔隙特征的大孔隙网络模型和微孔隙网络模型,在此基础上提出一种耦合算法构建出的同时描述大孔隙和微孔隙特征的碳酸盐岩双孔隙网络模型;然后,基于侵入-逾渗理论,模拟双孔隙网络模型中油水两相流体的一次驱替和二次吸吮过程,并建立了毛细管压力和相对渗透率的求解模型;最后,通过调整双孔隙网络结构参数,模拟水湿油藏条件下碳酸盐岩储层参数对相对渗透率曲线的影响. 结果表明,随着微孔隙比例因子和平均配位数的增加,油相相对渗透率曲线升高;随着双孔隙半径比的增加,油相和水相相对渗透率曲线下降,这对碳酸盐岩油藏渗流机理研究有着重要的指导意义.      

Pore to pore validation of pore network modelling against micromodel experiment results

Pore network modelling (PNM) has been widely used to study the multiphase flow and transport in porous media. Although a number of recent papers discussed the PNM validation on core-scale parameters such as permeability, relative permeability and capillary pressure; quantitative predictive potential of PNM on pore by pore basis has rarely been studied. The aim of this paper is to present a direct comparison between PNM simulations and corresponding micro-model experiments at the same scale and the same geometry. A number of well-defined and constrained two-phase flow in porous medium experimental scenarios were utilized to validate the physics solving part in PNM (filling rules, capillary and viscous pressure). This work validates that a dynamic pore network flow solver can predict two-phase flow displacements for these experiments for drainage situations at both pore and plug scales. A glass-etched micro-model is used to quantify the accuracy of a dynamic PNM solver on pore and core levels. Two-phase drainage micro fluidic experiments at different flow conditions are performed on micro-models. PNM simulations are performed on the same pattern and flow conditions as used in micro-model experiments. The two-phase distribution extracted from experiment images is registered onto rsults of PNM simulations for direct pore to pore comparison. Pore-scale matching level is found at around 75 % for all three test cases. The matching level of core-scale parameters such as S _{w c} and oil-phase permeability varies from case to case; the relative error to micro-model experiment measurements varies from 15 to 60 %. Possible reasons leading to discrepancies on core-scale parameters are discussed: missing considerations during validation of the combination of uncertainty in both simulator input parameters and experiments are seen as the principal factors.     

A dynamic network model for two‐phase immiscible flow

A dynamic pore‐scale network model is formulated for two‐phase immiscible flow. Interfaces are tracked through the pore throats using a modified Poiseuille equation, whereas special displacement rules are used at the pore bodies. The model allows interfaces to move over several pore‐lengths within a time step. Initial computational results are presented for a drainage experiment to demonstrate some of the features of the model. This revised version was published online in July 2006 with corrections to the Cover Date.     

Modeling of dynamic cohesive fracture propagation in porous saturated media

In this paper, a mathematical model is presented for the analysis of dynamic fracture propagation in the saturated porous media. The solid behavior incorporates a discrete cohesive fracture model, coupled with the flow in porous media through the fracture network. The double‐nodded zero‐thickness cohesive interface element is employed for the mixed mode fracture behavior in tension and contact behavior in compression. The crack is automatically detected and propagated perpendicular to the maximum effective stress. The spatial discretization is continuously updated during the crack propagation. Numerical examples from the hydraulic fracturing test and the concrete gravity dam show the capability of the model to simulate dynamic fracture propagation. The comparison is performed between the quasi‐static and fully dynamic solutions, and the performance of two analyses is investigated on the values of crack length and crack mouth opening. Copyright © 2010 John Wiley & Sons, Ltd.     

库水位下降时的岸坡非稳定渗流问题研究

水位下降时岸坡的渗流是涉及土体由饱和向非饱和状态过渡的水-气二相流过程,目前相关研究成果大都假设孔隙气压力为0,忽略孔隙气的影响。根据水、空气的质量守恒定律和达西定律,结合多相流理论建立水-气二相流模型,采用高效的积分有限差分法求解,通过变换主要变量,实现饱和（单相）与非饱和（二相）的相互转变,并给出各种边界条件下合理的数学处理方法。通过Muskat渗流问题,验证了上述模型的正确性;并对某土质岸坡水位下降时的非稳定渗流问题进行分析,结果表明,岸坡的基质吸力小于浸润线以上的负孔隙水压力,在浸润线以上的很大区域为毛细管水饱和带,其土体饱和且基质吸力为0,这对边坡稳定十分不利,精确分析水位下降的边坡稳定问题时,孔隙气压力变化的影响值得研究。     

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.     

考虑气相影响的降雨入渗过程分析研究

降雨入渗过程是水在下渗的过程中驱替空气的水-气二相流过程,对这一过程的精确模拟一直是渗流计算的难点,目前的处理方法通常是忽略孔隙气压力变化的影响。根据多相流理论,结合质量守恒定律和达西定律,建立了水-气二相流模型,模型的求解采用积分有限差分法和Newton-Raphson迭代方法,通过变换主要变量来表达相态的变化,实现了水相、气相边界条件及降雨入渗边界的精确模拟。利用上述模型对一土柱试验进行模拟,从而验证了模型的正确性,研究了一均质土层的降雨入渗过程,得到了孔隙水压力、孔隙气压力和毛细压力及含水率的变化过程。根据入渗率与地表孔隙气压力的变化关系,验证了孔隙气压力的增大对入渗水流产生阻滞作用。在求解非稳定渗流问题中,考虑空气压力变化的影响是值得研究的。     

On computation of strain‐dependent permeability of rocks and rock‐like porous media

Determination of transport properties of geomaterials is an important issue in many fields of engineering analysis and design. For example, in petroleum engineering, in situ permeability of an oil reservoir may be crucial in establishing its viability for exploitation, whilst prevention of leakage from underground storage facilities for oil and gas, nuclear waste as well as viability of CO₂ sequestration projects crucially depends on its long‐term values. Permeability is indirectly related to the porosity, pore‐size distribution and pore architecture of the porous media. These parameters evolve when a strain field is imposed. Physical measurement of permeability under a strain field in laboratory conditions is difficult, expensive and prone to a number of uncertainties. In the past, pore network models have been used to compute permeability of materials under stress/strain‐free conditions. In this paper, we propose an enhanced pore network model to compute permeability of rocks and rock‐like porous media under a stress/strain field. Data of pore‐size distribution obtained from mercury intrusion porosimetry are used to compute permeability of rock samples from various unspecified oilfields in the world. It is shown that the two permeabilities can be predicted from the model with sufficient accuracy. A hypothesis for change in porosity, pore‐size distribution and pore architecture as a result of imposed mechanical strains is then proposed. Based on this, permeability is computed again for one of the rock samples under uniaxial and triaxial compressive and tensile strain fields. It is shown that depending on the state of strain field imposed, permeability evolves in an anisotropic manner. Permeability under tensile strain field increases dramatically compared with the reduction that takes place under compressive strain field of the same magnitude. Copyright © 2014 John Wiley & Sons, Ltd.     

A micromechanics model for molecular diffusion in materials with complex pore structure

Molecular diffusion in fully saturated porous materials is strongly influenced by the pore space, which, in general, is characterized by a complex topological structure. Hence, information on macroscopic diffusion properties requires up‐scaling of transport processes within nano‐pores and micro‐pores over several spatial scales. A new model in the framework of continuum micromechanics is proposed for predicting the effective molecular diffusivity in porous materials. Considering a representative volume element, characterizing a porous material without any information about the pore space microstructure complexity, the uniform flux is perturbed by recursively embedding shape information hierarchically in the form of the ESHELBY matrix‐inclusion morphology to obtain the effective diffusivity as a function of the recurrence level and the porosity. The model predicts a threshold value for the porosity, below which no molecular diffusion can occur because of the presence of isolated pore clusters that are not connected and unavailable for transport. The maximum porosity, below which no molecular transport is possible, is predicted as one‐third for spherical inclusions. The model allows for extensions to more complex morphologies of the inclusions. We also identify, that the effects of the micro‐structure on molecular transport are characterized by porosity dependent long‐range and short‐range interactions. The developed framework is extended to incorporate realistic pore size distributions across several spatial scales by means of a distribution function within the hierarchical homogenization scheme. Available experimental results assert the model predictions. Copyright © 2015 John Wiley & Sons, Ltd.     

The effect of capillary pressure on the saturation equation of two‐phase flow in porous media

A complete and accurate simulation of two‐phase flow in porous media requires knowledge of all the controlling physics (and values of physical parameters) that play a relevant role and an understanding of the effects of each one on the solution. Of particular concern here is the effect of capillary pressure and the length scale over which it is relevant. The goal of this paper is to provide guidance onto when to include the effects of capillary pressure in the model, and onto what are the resulting length scale restrictions if those effects are to be included. Copyright © 2011 John Wiley & Sons, Ltd.     

An elasto‐plastic three phase model for partially saturated soil for the finite element simulation of compressed air support in tunnelling

This paper presents a fully coupled finite element formulation for partially saturated soil as a triphasic porous material, which has been developed for the simulation of shield tunnelling with heading face support using compressed air. While for many numerical simulations in geotechnics use of a two‐phase soil model is sufficient, the simulation of compressed air support demands the use of a three‐phase model with the consideration of air as a separate phase. A multiphase model for soft soils is developed, in which the individual constituents of the soil—the soil skeleton, the fluid and the gaseous phase—and their interactions are considered. The triphasic model is formulated within the framework of the theory of porous media, based upon balance equations and constitutive relations for the soil constituents and their mixture. An elasto‐plastic, cam–clay type model is extended to partially saturated soil conditions by incorporating capillary pressure according to the Barcelona basic model. The hydraulic properties of the soil are described via DARCY 's law and the soil–water characteristic curve after VAN GENUCHTEN . Water is modelled as an incompressible and air as a compressible phase. The model is validated by means of selected benchmark problems. The applicability of the model to geotechnical problems is demonstrated by results from the simulation of a compressed air intervention in shield tunnelling. Copyright © 2009 John Wiley & Sons, Ltd.     

Frictionless contact formulation for dynamic analysis of nonlinear saturated porous media based on the mortar method

A finite element algorithm for frictionless contact problems in a two‐phase saturated porous medium, considering finite deformation and inertia effects, has been formulated and implemented in a finite element programme. The mechanical behaviour of the saturated porous medium is predicted using mixture theory, which models the dynamic advection of fluids through a fully saturated porous solid matrix. The resulting mixed formulation predicts all field variables including the solid displacement, pore fluid pressure and Darcy velocity of the pore fluid. The contact constraints arising from the requirement for continuity of the contact traction, as well as the fluid flow across the contact interface, are enforced using a penalty approach that is regularised with an augmented Lagrangian method. The contact formulation is based on a mortar segment‐to‐segment scheme that allows the interpolation functions of the contact elements to be of order N. The main thrust of this paper is therefore how to deal with contact interfaces in problems that involve both dynamics and consolidation and possibly large deformations of porous media. The numerical algorithm is first verified using several illustrative examples. This algorithm is then employed to solve a pipe‐seabed interaction problem, involving large deformations and dynamic effects, and the results of the analysis are also compared with those obtained using a node‐to‐segment contact algorithm. The results of this study indicate that the proposed method is able to solve the highly nonlinear problem of dynamic soil–structure interaction when coupled with pore water pressures and Darcy velocity. Copyright © 2015 John Wiley & Sons, Ltd.     

A simulation study of capillary transport,preferential retention and distribution of salts in historic sandstone buildings

Salt crystallisation is a major problem of deterioration in historic stone buildings, monuments and sculptures. The capillary rise of soil water is one of the primary sources of salts in stone structures, which evaporates leaving the salts behind. It has been noted that the spatial distribution profile of different species of salts crystallised in historic stone buildings is not homogeneous, i.e. different salts crystallise at different locations. The capillary transport and inhomogeneous spatial distribution of different salts in the porous building materials has been considered to be a result of solubility-dependent crystallisation; however, the factors responsible for this phenomenon are not clearly known. This paper aims to investigate the factors influencing the differential distribution of salts during capillary rise of soil water. In this study, the capillary transport of salts was simulated on two different sandstones—Locharbriggs, a Permo–Triassic, red sandstone and Stoke Hall, a Carboniferous, buff sandstone. The experiments were carried out under controlled environmental conditions to eliminate the possibility of evaporation-driven crystallisation of salts depending on their solubilities. The results indicate that fractionation or differential distribution of salts takes place even in the absence of evaporation and crystallisation. The sandstones exhibit properties like an ion exchange column, and ionic species present in the salt solution show differential distribution within the porous network of sandstone.     

A pore-throat model based on grain-size distribution to quantify gravity-dominated DNAPL instabilities in a water-saturated homogeneous porous medium

A pore-scale network model based on spherical pore bodies and cylindrical pore throats was developed to describe the displacement of water by DNAPL. The pore body size and the pore throat size were given by statistical distributions with user-specified values for the minimum, mean and maximum sizes. The numerical model was applied to a laboratory experiment conducted on a sand-filled glass column. The parameters relative to pore body size and pore throat size that were used in the construction of the equivalent network were derived from the discrete grain-size distribution of the real porous medium. The calculated arrival times of the DNAPL front were compared with those measured using optic fibre sensors placed at different points on the control section of the experimental device. Furthermore, the model simulated DNAPL pressure measured at the entrance section of the system. In general, the numerical results obtained with the model were in good agreement with the actual measurements.     

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

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.     

Advances in Air Sparging Technology of Saturated Zone

In situ air sparging involves injecting atmospheric air, under pressure, into the saturated zone to remove those volatile and semi-volatile organic groundwater contaminants and to promote their biodegradation by increasing subsurface oxygen concentrations. Due to the advantages of low cost, high efficiency and in situ constructability, groundwater Air Sparging (AS) technology has been quickly developed in the world recently. Based on the explanation of its remediation principle, literature review is done on the research advancement of air sparging technology mainly from three aspects. First, various methods for determination of the zone of influence and visualization techniques of air flow forms during air sparging are summarized. Then the influence of environmental geological conditions and construction technology parameters on the remediation effect of air sparging is systematically analyzed. Thereafter, two main types of air sparging theoretical models including lumped-parameter model and multiphase fluid flow model are discussed respectively in detail. Finally, based on the problems and difficulties existing in present research and engineering practice, several future tasks such as the enhancement remediation techniques in complex geological sites, microscopic intrinsic mechanisms, and establishment of related design and construction standards which require to be done are briefly analyzed.