Multiscale Extended Finite Element Method for Deformable Fractured Porous Media

Deformable fractured porous media appear in many geoscience applications. While the extended finite element (XFEM) has been successfully developed within the computational mechanics community for accurate modeling of the deformation, its application in natural geoscientific applications is not straightforward. This is mainly due to the fact that subsurface formations are heterogeneous and span large length scales with many fractures at different scales. In this work, we propose a novel multiscale formulation for XFEM, based on locally computed basis functions. The local multiscale basis functions capture the heterogeneity and discontinuities introduced by fractures. Local boundary conditions are set to follow a reduced-dimensional system, in order to preserve the accuracy of the basis functions. Using these multiscale bases, a multiscale coarse-scale system is then governed algebraically and solved, in which no enrichment due to the fractures exist. Such formulation allows for significant computational cost reduction, at the same time, it preserves the accuracy of the discrete displacement vector space. The coarse-scale solution is finally interpolated back to the fine scale system, using the same multiscale basis functions. The proposed multiscale XFEM (MS-XFEM) is also integrated within a two-stage algebraic iterative solver, through which error reduction to any desired level can be achieved. Several proof-of-concept numerical tests are presented to assess the performance of the developed method. It is shown that the MS-XFEM is accurate, when compared with the fine-scale reference XFEM solutions. At the same time, it is significantly more efficient than the XFEM on fine-scale resolution. As such, it develops the first scalable XFEM method for large-scale heavily fractured porous media.

[1]  Yalchin Efendiev,et al.  Generalized Multiscale Finite Element Methods. Oversampling Strategies , 2013, 1304.4888.

[2]  Julian J. Bommer,et al.  Induced seismicity associated with Enhanced Geothermal Systems , 2007 .

[3]  Knut-Andreas Lie,et al.  Dynamic Coarsening and Local Reordered Nonlinear Solvers for Simulating Transport in Porous Media , 2020, SPE Journal.

[4]  Hong-wu Zhang,et al.  Extended multiscale finite element method for mechanical analysis of heterogeneous materials , 2010 .

[5]  M. Rashid The arbitrary local mesh replacement method: An alternative to remeshing for crack propagation analysis , 1998 .

[6]  Hadi Hajibeygi,et al.  Multiscale gradient computation for flow in heterogeneous porous media , 2017, J. Comput. Phys..

[7]  Iuliu Sorin Pop,et al.  Upscaling of reactive flows in domains with moving oscillating boundaries , 2012 .

[8]  Jiamin Jiang,et al.  Hybrid Coupled Discrete-Fracture/Matrix and Multicontinuum Models for Unconventional-Reservoir Simulation , 2016 .

[9]  Hamdi A. Tchelepi,et al.  Hybrid Multiscale Formulation for Coupled Flow and Geomechanics , 2016 .

[10]  A. Khoei,et al.  Fully coupled hydromechanical multiscale model with microdynamic effects , 2018 .

[11]  Hamdi A. Tchelepi,et al.  Multiscale two-stage solver for Biot’s poroelasticity equations in subsurface media , 2018, Computational Geosciences.

[12]  A macroelement stabilization for mixed finite element/finite volume discretizations of multiphase poromechanics , 2020, Computational Geosciences.

[13]  Hamdi A. Tchelepi,et al.  Unified thermo-compositional-mechanical framework for reservoir simulation , 2018, Computational Geosciences.

[14]  A. Khoei Extended Finite Element Method: Theory and Applications , 2015 .

[15]  Patrick Jenny,et al.  Iterative multiscale finite-volume method , 2008, J. Comput. Phys..

[16]  Ted Belytschko,et al.  A finite element method for crack growth without remeshing , 1999 .

[17]  Hui Zhou,et al.  Two-Stage Algebraic Multiscale Linear Solver for Highly Heterogeneous Reservoir Models , 2012, Annual Simulation Symposium.

[18]  Olav Møyner,et al.  The multiscale restriction smoothed basis method for fractured porous media (F-MsRSB) , 2016, J. Comput. Phys..

[19]  Robert Davis Cook,et al.  Finite Element Modeling for Stress Analysis , 1995 .

[20]  Hadi Hajibeygi,et al.  Algebraic multiscale method for flow in heterogeneous porous media with embedded discrete fractures (F-AMS) , 2016, J. Comput. Phys..

[21]  Cornelis Vuik,et al.  Algebraic dynamic multilevel method for embedded discrete fracture model (F-ADM) , 2018, J. Comput. Phys..

[22]  B. K. Mishra,et al.  A new multiscale XFEM for the elastic properties evaluation of heterogeneous materials , 2017 .

[23]  Olav Møyner,et al.  Use of Multiple Multiscale Operators To Accelerate Simulation of Complex Geomodels , 2017 .

[24]  Patrick Jenny,et al.  A hierarchical fracture model for the iterative multiscale finite volume method , 2011, J. Comput. Phys..

[25]  Yixuan Wang,et al.  Adaptive algebraic multiscale solver for compressible flow in heterogeneous porous media , 2014, J. Comput. Phys..

[26]  Amir R. Khoei,et al.  A mesh-independent finite element formulation for modeling crack growth in saturated porous media based on an enriched-FEM technique , 2014, International Journal of Fracture.

[27]  Jian-Ying Wu,et al.  An improved stable XFEM (Is-XFEM) with a novel enrichment function for the computational modeling of cohesive cracks , 2015 .

[28]  Anna Scotti,et al.  Unfitted FEM for Modelling the Interaction of Multiple Fractures in a Poroelastic Medium , 2017 .

[29]  M. Zoback,et al.  Fluid flow along potentially active faults in crystalline rock: Geology , 1995 .

[30]  Kes Heffer,et al.  Theory of linear poroelasticity with applications to geomechanics and hydrogeology , 2004 .

[31]  Yalchin Efendiev,et al.  An adaptive GMsFEM for high-contrast flow problems , 2013, J. Comput. Phys..

[32]  Jostein R. Natvig,et al.  Constrained pressure residual multiscale (CPR-MS) method for fully implicit simulation of multiphase flow in porous media , 2015, J. Comput. Phys..

[33]  Stein Krogstad,et al.  A Hierarchical Multiscale Method for Two-Phase Flow Based upon Mixed Finite Elements and Nonuniform Coarse Grids , 2006, Multiscale Model. Simul..

[34]  Kirill M. Terekhov,et al.  Cell-centered finite-volume method for heterogeneous anisotropic poromechanics problem , 2020, J. Comput. Appl. Math..

[35]  E Weinan,et al.  Finite difference heterogeneous multi-scale method for homogenization problems , 2003 .

[36]  A. Khoei,et al.  A bridge between dual porosity and multiscale models of heterogeneous deformable porous media , 2018, International Journal for Numerical and Analytical Methods in Geomechanics.

[37]  L. J. Sluys,et al.  A new method for modelling cohesive cracks using finite elements , 2001 .

[38]  L. J. Sluys,et al.  Three-dimensional embedded discontinuity model for brittle fracture , 2001 .

[39]  Hamdi A. Tchelepi,et al.  Algebraic multiscale solver for flow in heterogeneous porous media , 2014, J. Comput. Phys..

[40]  F. Gaspar,et al.  A finite difference analysis of Biot's consolidation model , 2003 .

[41]  Hadi Hajibeygi,et al.  Projection-based Embedded Discrete Fracture Model (pEDFM) , 2017 .

[42]  U. Hornung Homogenization and porous media , 1996 .

[43]  K. Bisdom,et al.  Calibrating discrete fracture-network models with a carbonate three-dimensional outcrop fracture network: Implications for naturally fractured reservoir modeling , 2014 .

[44]  T. Belytschko,et al.  Element‐free Galerkin methods , 1994 .

[45]  Yerlan Amanbek,et al.  Adaptive numerical homogenization for upscaling single phase flow and transport , 2019, J. Comput. Phys..

[46]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[47]  Y. Saad,et al.  Experimental study of ILU preconditioners for indefinite matrices , 1997 .

[48]  Paul A. Wawrzynek,et al.  Quasi-automatic simulation of crack propagation for 2D LEFM problems , 1996 .

[49]  Alessio Fumagalli,et al.  An Efficient XFEM Approximation of Darcy Flows in Arbitrarily Fractured Porous Media , 2014 .

[50]  G. Allaire Homogenization and two-scale convergence , 1992 .

[51]  P. Jenny,et al.  Finite volume–based modeling of flow‐induced shear failure along fracture manifolds , 2017 .

[52]  Patrick Jenny,et al.  Multiscale finite-volume method for parabolic problems arising from compressible multiphase flow in porous media , 2009, J. Comput. Phys..

[53]  L. J. Sluys,et al.  A phantom node formulation with mixed mode cohesive law for splitting in laminates , 2009 .

[54]  Ted Belytschko,et al.  Elastic crack growth in finite elements with minimal remeshing , 1999 .

[55]  H. Tchelepi,et al.  Multi-scale finite-volume method for elliptic problems in subsurface flow simulation , 2003 .

[56]  Guotong Ren,et al.  A Fully Coupled XFEM-EDFM Model for Multiphase Flow and Geomechanics in Fractured Tight Gas Reservoirs , 2016, ICCS.

[57]  Hamdi A. Tchelepi,et al.  Multiscale finite-element method for linear elastic geomechanics , 2017, J. Comput. Phys..

[58]  Hadi Hajibeygi,et al.  Multiscale finite volume method for finite-volume-based simulation of poroelasticity , 2019, J. Comput. Phys..

[59]  Cornelis Vuik,et al.  Adaptive dynamic multilevel simulation of fractured geothermal reservoirs , 2020, J. Comput. Phys. X.

[60]  Angelo Simone,et al.  The Discontinuity‐Enriched Finite Element Method , 2017 .

[61]  I. Babuska,et al.  The partition of unity finite element method: Basic theory and applications , 1996 .

[62]  A. Kozlova,et al.  A Real-Field Multiscale Black-Oil Reservoir Simulator , 2015, ANSS 2015.