Embedded discontinuity finite element modeling of fluid flow in fractured porous media

In many geomaterials, particularly rocks and clays, permeability is greatly enhanced by the presence of fractures. Fracture sets create an overall permeability that is anisotropic, enhanced in the directions of the fractures. In modeling the fractures via a finite element method, for example, meshing around these fractures can become quite difficult and result in computationally intensive systems. In this article, we develop a relatively simple method for including the fractures within the elements. Flow through the bulk medium is assumed to be governed by Darcy’s law, and the flow on the fracture by a generalization of the law. This model is embedded in a finite element framework, with the fractures passing through the elements. In this formulation, elements with fractures are given an enhanced permeability in the direction of the fractures. With these enhancements, the material essentially becomes anisotropically more permeable in the direction of fracture sets.     

A finite element method for modeling coupled flow and deformation in porous fractured media

Modeling the flow in highly fractured porous media by finite element method (FEM) has met two difficulties: mesh generation for fractured domains and a rigorous formulation of the flow problem accounting for fracture/matrix, fracture/fracture, and fracture/boundary fluid mass exchanges. Based on the recent theoretical progress for mass balance conditions in multifractured porous bodies, the governing equations for coupled flow and deformation in these bodies are first established in this paper. A weak formulation for this problem is then established allowing to build a FEM. Taking benefit from recent development of mesh‐generating tools for fractured media, this weak formulation has been implemented in a numerical code and applied to some typical problems of hydromechanical coupling in fractured porous media. It is shown that in this way, the FEM that has proved its efficiency to model hydromechanical phenomena in porous media is extended with all its performances (calculation time, couplings, and nonlinearities) to fractured porous media. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical modeling of the flow and transport of radionuclides in heterogeneous porous media

This paper is concerned with numerical methods for the modeling of flow and transport of contaminant in porous media. The numerical methods feature the mixed finite element method over triangles as a solver to the Darcy flow equation and a conservative finite volume scheme for the concentration equation. The convective term is approximated with a Godunov scheme over the dual finite volume mesh, whereas the diffusion–dispersion term is discretized by piecewise linear conforming triangular finite elements. It is shown that the scheme satisfies a discrete maximum principle. Numerical examples demonstrate the effectiveness of the methodology for a coupled system that includes an elliptic equation and a diffusion–convection–reaction equation arising when modeling flow and transport in heterogeneous porous media. The proposed scheme is robust, conservative, efficient, and stable, as confirmed by numerical simulations.      

Multibody approach for reactive transport modeling in discontinuous-heterogeneous porous media

Computational Geosciences - In the context of long-term degradation of porous media, the coupling between fracture mechanics and reactive transport is investigated. A reactive transport model in a...     

A fracture mapping and extended finite element scheme for coupled deformation and fluid flow in fractured porous media

This paper presents a fracture mapping (FM) approach combined with the extended finite element method (XFEM) to simulate coupled deformation and fluid flow in fractured porous media. Specifically, the method accurately represents the impact of discrete fractures on flow and deformation, although the individual fractures are not part of the finite element mesh. A key feature of FM‐XFEM is its ability to model discontinuities in the domain independently of the computational mesh. The proposed FM approach is a continuum‐based approach that is used to model the flow interaction between the porous matrix and existing fractures via a transfer function. Fracture geometry is defined using the level set method. Therefore, in contrast to the discrete fracture flow model, the fracture representation is not meshed along with the computational domain. Consequently, the method is able to determine the influence of fractures on fluid flow within a fractured domain without the complexity of meshing the fractures within the domain. The XFEM component of the scheme addresses the discontinuous displacement field within elements that are intersected by existing fractures. In XFEM, enrichment functions are added to the standard finite element approximation to adequately resolve discontinuous fields within the simulation domain. Numerical tests illustrate the ability of the method to adequately describe the displacement and fluid pressure fields within a fractured domain at significantly less computational expense than explicitly resolving the fracture within the finite element mesh. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical study of two‐phase fluid distributions in fractured porous media

Two‐phase fluid distributions in fractured porous media were studied using a single‐component multiphase (SCMP) lattice Boltzmann method (LBM), which was selected among three commonly used numerical approaches through a comparison against the available results of micro X‐ray computed tomography. The influence of the initial configuration and the periodic boundary conditions in the SCMP LBM for the fluid distribution analysis were investigated as well. It was revealed that regular porous media are sensitive to the initial distribution, whereas irregular porous media are insensitive. Moreover, to eliminate the influence of boundaries, the model's buffer size of an SCMP LBM simulation was suggested to be taken as approximately 12.5 times the average particle size. Then, the two‐phase fluid distribution of a porous medium was numerically studied using the SCMP LBM. Both detailed distribution patterns and macroscopic morphology parameters were reasonably well captured. Finally, the two‐phase fluid distributions in a fractured porous media were investigated. The influence of the degree of saturation, fracture length, and fracture width on the fluid distributions and migration was explored. Copyright © 2015 John Wiley & Sons, Ltd.     

模拟裂隙多孔介质中变饱和渗流的广义等效连续体方法

描述了一种计算裂隙多孔介质中变饱和渗流的广义等效连续体方法。这种方法忽略裂隙的毛细作用,设定一个与某孔隙饱和度相对应的综合饱和度极限值,并假定：（1）如果裂隙多孔介质的综合饱和度小于该极限值,水只在孔隙中存在并流动,而裂隙中则没有水的流动;（2）如果综合饱和度等于或大于该极限值,水将进入裂隙,并在裂隙内运动。分析比较了等效连续体模型的不同计算方法,并给出了一个模拟裂隙岩体中变饱和渗流与传热耦合问题的应用算例。结果表明,所述方法具有一般性,可以有效地模拟裂隙多孔介质中变饱和渗流的基本特征。     

Numerical simulation of liquefaction in porous media using nonlinear fluid flow law

It is well known that for a sufficiently high seepage velocity, the governing flow law of porous media is nonlinear (J. Computers & Fluids 2010; 39 : 2069–2077). However, this fact has not been considered in the studies of soil‐pore fluid interaction and in conventional soil mechanics. In the present paper, a fully explicit dynamic finite element method is developed for nonlinear Darcy law. The governing equations are expressed for saturated porous media based on the extension of the Biot (J. Appl. Phys. 1941; 12 : 155–164) formulation. The elastoplastic behavior of soil under earthquake loading is simulated using a generalized plasticity theory that is composed of a yield surface along with non‐associated flow rule. Numerical simulations of porous media subjected to horizontal and vertical components of ground motion excitations with different permeability coefficients are carried out; while computed maximum pore water pressure is specially taken into consideration to make the difference between Darcy and non‐Darcy flow regimes tangible. Finally, the effect of non‐Darcy flow on the evaluated liquefaction potential of sand in comparison to conventional Darcy law is examined. Copyright © 2014 John Wiley & Sons, Ltd.     

Boundary element analysis of contaminant transport in fractured porous media

A two-dimensional boundary integral method to analyse the flow of contaminant in fractured media having a two- or three-dimensional orthogonal fracture network is presented. The method assumes that the fractures provide the paths of least resistance for transport of contaminants while the matrix, because of its low permeability, acts as ‘storage blocks’ into which the contaminant diffuses. Laplace transform is used to eliminate the time variable in the governing equation in order to facilitate the formulation of a boundary integral equation in the Laplace transform space. Conventional boundary element techniques are applied to solve for the contaminant concentrations at specified locations in the spatial domain. The concentration in the time domain is then obtained by using an efficient inversion technique developed by Talbot. The method is able to analyse the behaviour of waste repositories which have diminishing concentration due to the mass transport of the contaminant into the surrounding fractured media.     

Numerical evaluation of multicomponent cation exchange reactive transport in physically and geochemically heterogeneous porous media

Experimental evidence and stochastic studies strongly show that the transport of reactive solutes in porous media is significantly influenced by heterogeneities in hydraulic conductivity, porosity, and sorption parameters. In this paper, we present Monte Carlo numerical simulations of multicomponent reactive transport involving competitive cation exchange reactions in a two-dimensional vertical physically and geochemically heterogeneous medium. Log hydraulic conductivity, log K, and log cation exchange capacity (log CEC) are assumed to be random Gaussian functions with spherical semivariograms. Random realizations of log K and log CEC are used as input data for the numerical simulation of multicomponent reactive transport with CORE^2D, a general purpose reactive transport code. Longitudinal features of the fronts of reactive and conservative species are computed from the temporal and spatial moments of depth-averaged concentrations. Monte Carlo simulations show that: (1) the displacement of reactive fronts increases with increasing variance of log K, while it decreases with the variance of log CEC; (2) second-order spatial moments increase with increasing variances of log K and log CEC; (3) uncertainties in the mean arrival time are largest (smallest) for negatively (positively) correlated log K and Log CEC; (4) cations undergoing competitive cation exchange exhibit different apparent velocities and retardation factors due to both physical and geochemical heterogeneities; and (5) the correlation between log K and log CEC affects significantly apparent cation retardation factors in heterogeneous aquifers.     

Transverse and longitudinal fluid flow modelling in fractured porous media with non‐matching meshes

A new discrete fracture model is introduced to simulate the steady‐state fluid flow in discontinuous porous media. The formulation uses a multi‐layered approach to capture the effect of both longitudinal and transverse permeability of the discontinuities in the pressure distribution. The formulation allows the independent discretisation of mesh and discontinuities, which do not need to conform. Given that the formulation is developed at the element level, no additional degrees of freedom or special integration procedures are required for coupling the non‐conforming meshes. The proposed model is shown to be reliable regardless of the permeability of the discontinuity being higher or lower than the surrounding domain. Four numerical examples of increasing complexity are solved to demonstrate the efficiency and accuracy of the new technique when compared with results available in the literature. Results show that the proposed method can simulate the fluid pressure distribution in fractured porous media. Furthermore, a sensitivity analysis demonstrated the stability regarding the condition number for wide range values of the coupling parameter.     

Smoothed particle hydrodynamics and its applications for multiphase flow and reactive transport in porous media

Smoothed particle hydrodynamics (SPH) is a Lagrangian method based on a meshless discretization of partial differential equations. In this review, we present SPH discretization of the Navier-Stokes and advection-diffusion-reaction equations, implementation of various boundary conditions, and time integration of the SPH equations, and we discuss applications of the SPH method for modeling pore-scale multiphase flows and reactive transport in porous and fractured media.     

Numerical modeling of two-phase fluid flow and oil slick transport in estuarine water

Oil spills is one of the most important hazards in the estuarine and coastal water. In recent decades, engineers try to predict the status of oil slick to manage the pollution spreading. The prediction of oil slick transport is carried out mainly by means of numerical models. In the current study, the development and application of a two-phase fluid flow model to simulate oil transport in the marine environment are presented. Different transport and fate processes are included in the developed model. The model consists of the Lagrangian method for the advection process, the Random Walk technique for horizontal diffusion process and the empirical equations for the fate processes. The major forces for driving oil particles are fluid current, wind speed and turbulent flow. Therefore, the multi-component hydrocarbon method has been included to the developed model in order to predict fate processes. As prediction of particle velocity components is of major importance for oil slick advection, therefore the binomial interpolation procedure has been chosen for the particle velocity components computations. In addition, shoreline boundary condition is included in the developed model to simulate shore response to oil slick transport near the beaches. The results of the model applications are compared with the analytical solutions, experimental measurements and other numerical models cited in literature. Comparisons of different sets of results represent the capability of developed model to predict the oil slick transport. In addition, the developed model is tested for two oil spill cases in the Persian Gulf.     

Embedded discontinuity approach for coupled hydromechanical analysis of fractured porous media

In this paper, a new continuum approach for the coupled hydromechanical analysis of fractured porous media is proposed. The methodology for describing the hydraulic characteristics invokes an enriched form of Darcy's law formulated in the presence of an embedded discontinuity. The constitutive relations governing the hydromechanical response are derived by averaging the fluid pressure gradient and the discontinuous displacement fields over a selected referential volume of the material, subject to some physical constraints. The framework incorporates an internal length scale which is explicitly embedded in the definition of gradient operators. The respective field equations are derived following the general form of balance equations in interacting continua. The conventional finite element method is then employed for the spatial discretization, and the generalized Newmark scheme is used for the temporal discretization. The proposed methodology is verified by some numerical examples dealing with a steady-state flow through fractured media as well as a time-dependent consolidation in the presence of a discontinuity.     

Numerical solutions for flow in porous media

A numerical approach is proposed to model the flow in porous media using homogenization theory. The proposed concept involves the analyses of micro‐true flow at pore‐level and macro‐seepage flow at macro‐level. Macro‐seepage and microscopic characteristic flow equations are first derived from the Navier–Stokes equation at low Reynolds number through a two‐scale homogenization method. This homogenization method adopts an asymptotic expansion of velocity and pressure through the micro‐structures of porous media. A slightly compressible condition is introduced to express the characteristic flow through only characteristic velocity. This characteristic flow is then numerically solved using a penalty FEM scheme. Reduced integration technique is introduced for the volumetric term to avoid mesh locking. Finally, the numerical model is examined using two sets of permeability test data on clay and one set of permeability test data on sand. The numerical predictions agree well with the experimental data if constraint water film is considered for clay and two‐dimensional cross‐connection effect is included for sand. Copyright © 2003 John Wiley & Sons, Ltd.     

Thermo-hydro-mechanical modeling of fracturing porous media with two-phase fluid flow using X-FEM technique

In this paper, a fully coupled thermo-hydro-mechanical model is presented for two-phase fluid flow and heat transfer in fractured/fracturing porous media using the extended finite element method. In the fractured porous medium, the traction, heat, and mass transfer between the fracture space and the surrounding media are coupled. The wetting and nonwetting fluid phases are water and gas, which are assumed to be immiscible, and no phase-change is considered. The system of coupled equations consists of the linear momentum balance of solid phase, wetting and nonwetting fluid continuities, and thermal energy conservation. The main variables used to solve the system of equations are solid phase displacement, wetting fluid pressure, capillary pressure, and temperature. The fracture is assumed to impose the strong discontinuity in the displacement field and weak discontinuities in the fluid pressure, capillary pressure, and temperature fields. The mode I fracture propagation is employed using a cohesive fracture model. Finally, several numerical examples are solved to illustrate the capability of the proposed computational algorithm. It is shown that the effect of thermal expansion on the effective stress can influence the rate of fracture propagation and the injection pressure in hydraulic fracturing process. Moreover, the effect of thermal loading is investigated properly on fracture opening and fluids flow in unsaturated porous media, and the convective heat transfer within the fracture is captured successfully. It is shown how the proposed computational model is capable of modeling the fully coupled thermal fracture propagation in unsaturated porous media.     

Numerical modelling of coupled flow and deformation in fractured rock specimens

A dual-porosity poroelastic model is extended to represent behaviour in cylindrical co-ordinates for the evaluation of flow-deformation effects in cylindrical laboratory samples incorporating a central wellbore or non-repeating axisymmetric injection on the periphery. Nine-node quadratic elements are used to represent mechanical deformation, while eight-node linear elements are used to interpolate the pressure fields, which offers significant advantages over the behaviour of other non-conforming elements. The model presented is validated against simplified analytical results, and extended to describe the behaviour of homogeneous and heterogeneous laboratory specimens subjected to controlled triaxial state of stress and injection tests. Apparent from the results is the significant influence of stress-deformation effects over system behaviour. Copyright © 1999 John Wiley & Sons, Ltd.     

Finite element simulation of fluid flow in fractured rock media

Fluid dynamics models are used by the petroleum industry to model single- and/or multi-phase flow within fractured rock formations, in order to facilitate extraction of fluids such as oil and natural gas, and in other areas of engineering to study groundwater flow, as well as to estimate contaminant seepage and transport. In this paper, the numerical modelling software Comsol is used to simulate air and water flow through a specimen of granite with a single vertical fracture subjected to triaxial loading conditions. The intent of the model is to simulate triaxial test findings on a rock specimen with a natural fracture. Fluid flow is simulated at various confining and inlet pressures using the cubic law. Model results were in good agreement with laboratory findings. Pressure distribution along the fracture and across the specimen are as expected with a near linear pressure distribution along the length of the fracture. A drawdown effect on pressure distribution across the specimen in the vicinity of the fracture is also observed. Pressure gradient was largely uniform; however, some localised zones of high gradient along the fracture are observed.     

含喉部裂隙介质CO_2反应迁移的格子Boltzmann模拟研究

CO_2地质封存是目前最经济、最可靠的CO_2减排技术之一,对CO_2反应迁移规律的研究具有重要的理论价值。采用Dardis多孔介质模型与已有的CO_2反应迁移模型耦合,对含喉部裂隙的复杂介质中反应迁移规律进行模拟研究。速度场的模拟结果表明,裂隙内的速度明显高于基质速度,在喉部中心线位置处速度达到最大值。溶解反应主要集中于入口段及裂隙的上、下边缘附近;受喉部的影响,喉部下游的裂隙边缘几乎不发生溶解反应。反应物H~+在裂隙及喉部中的浓度明显高于在基质中的浓度;而生成物Ca~(2+)的高浓度区则出现在下游基质区中。针对不同喉部位置的情况进行对比,分析其对反应率及组分浓度的影响规律。最后对含突扩孔的裂隙介质进行了模拟,发现其提高了裂隙内的速度及组分迁移,这与喉部裂隙介质内的规律正好相反。上述结果较好地说明了本模型具有模拟研究复杂裂隙介质内的CO_2反应迁移规律的能力。     

A mixed solution for two-dimensional unsteady flow in fractured porous media

A mixed finite element–boundary element solution for the analysis of two-dimensional flow in porous media composed of rock blocks and discrete fractures is described. The rock blocks are modelled implicitly by using boundary elements whereas finite elements are adopted to model the discrete fractures. The computational procedure has been implemented in a hybrid code which has been validated first by comparing the numerical results with the closed-form solution for flow in a porous aquifer intercepted by a vertical fracture only. Then, a more complex problem has been solved where a pervious, homogeneous and isotropic matrix containing a net of fractures is considered. The results obtained are shown to describe satisfactorily the main features of the flow problem under study. © 1997 by John Wiley & Sons, Ltd.