Two-grid finite element method for the dual-permeability-Stokes fluid flow model

In this paper, two-grid finite element method for the steady dual-permeability-Stokes fluid flow model is proposed and analyzed. Dual-permeability-Stokes interface system has vast applications in many areas such as hydrocarbon recovery process, especially in hydraulically fractured tight/shale oil/gas reservoirs. Two-grid method is popular and convenient to solve a large multiphysics interface system by decoupling the coupled problem into several subproblems. Herein, the two-grid approach is used to reduce the coding task substantially, which provides computational flexibility without losing the approximate accuracy. Firstly, we solve a global problem through standard Pk − Pk− 1 − Pk − Pk finite elements on the coarse grid. After that, a coarse grid solution is applied for the decoupling between the interface terms and the mass exchange terms to solve three independent subproblems on the fine grid. The three independent parallel subproblems are the Stokes equations, the microfracture equations, and the matrix equations, respectively. Four numerical tests are presented to validate the numerical methods and illustrate the features of the dual-permeability-Stokes model.

[1]  Zhangxin Chen,et al.  Transient flow of slightly compressible fluids through double-porosity, double-permeability systems — A state-of-the-art review , 1989 .

[2]  Xiao-Chuan Cai,et al.  A scalable nonlinear fluid-structure interaction solver based on a Schwarz preconditioner with isogeometric unstructured coarse spaces in 3D , 2017, J. Comput. Phys..

[3]  Xiaoming He,et al.  A Data Assimilation Enabled Model for Coupling Dual Porosity Flow with Free Flow , 2018, 2018 17th International Symposium on Distributed Computing and Applications for Business Engineering and Science (DCABES).

[4]  Wei Liu,et al.  Two-grid finite volume element methods for semilinear parabolic problems , 2010 .

[5]  Jinchao Xu.  Two-grid Discretization Techniques for Linear and Nonlinear PDEs , 1996 .

[6]  Olivier Gosselin,et al.  Matrix-Fracture Transfer Function in Dual-Media Flow Simulation: Review, Comparison and Validation , 2008 .

[7]  Bin Jiang.  A parallel domain decomposition method for coupling of surface and groundwater flows , 2009 .

[8]  E. Miglio,et al.  Mathematical and numerical models for coupling surface and groundwater flows , 2002 .

[9]  Li Shan,et al.  Partitioned Time Stepping Method for a Dual-Porosity-Stokes Model , 2018, J. Sci. Comput..

[10]  Bin Yuan,et al.  A composite dual-porosity fractal model for channel-fractured horizontal wells , 2018 .

[11]  Jinchao Xu,et al.  Numerical Solution to a Mixed Navier-Stokes/Darcy Model by the Two-Grid Approach , 2009, SIAM J. Numer. Anal..

[12]  Xiaoming He,et al.  A domain decomposition method for the time-dependent Navier-Stokes-Darcy model with Beavers-Joseph interface condition and defective boundary condition , 2020, J. Comput. Phys..

[13]  Weidong Zhao,et al.  Finite Element Approximations for Stokes–darcy Flow with Beavers–joseph Interface Conditions * , 2022 .

[14]  Tong Zhang,et al.  Two novel decoupling algorithms for the steady Stokes-Darcy model based on two-grid discretizations , 2014 .

[15]  Jinchao Xu,et al.  A Two-Grid Method of a Mixed Stokes-Darcy Model for Coupling Fluid Flow with Porous Media Flow , 2007, SIAM J. Numer. Anal..

[16]  Yanren Hou,et al.  A decoupling two‐grid algorithm for the mixed Stokes‐Darcy model with the Beavers‐Joseph interface condition , 2014 .

[17]  Jirka Simunek,et al.  Single-porosity and dual-porosity modeling of water flow and solute transport in subsurface-drained fields using effective field-scale parameters , 2005 .

[18]  Md. Abdullah Al Mahbub,et al.  Two-Level Finite Element Approximation for Oseen Viscoelastic Fluid Flow , 2018 .

[19]  Haibiao Zheng.  Two-grid finite element method for the stabilization of mixed Stokes-Darcy model , 2016, Discrete & Continuous Dynamical Systems - B.

[20]  G. I. Barenblatt,et al.  Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata] , 1960 .

[21]  Júlio Carneiro,et al.  Numerical simulations on the influence of matrix diffusion to carbon sequestration in double porosity fissured aquifers , 2009 .

[22]  Todd Arbogast,et al.  Derivation of the double porosity model of single phase flow via homogenization theory , 1990 .

[23]  Li Shan,et al.  Partitioned Time Stepping Method for Fully Evolutionary Stokes-Darcy Flow with Beavers-Joseph Interface Conditions , 2013, SIAM J. Numer. Anal..

[24]  Xiaoming He,et al.  An artificial compressibility ensemble algorithm for a stochastic Stokes‐Darcy model with random hydraulic conductivity and interface conditions , 2019, International Journal for Numerical Methods in Engineering.

[25]  Xiaoming He,et al.  Parallel, non-iterative, multi-physics domain decomposition methods for time-dependent Stokes-Darcy systems , 2014, Math. Comput..

[26]  M. Moraveji,et al.  Experimental and simulation study of gas diffusion effect during gas injection into naturally fractured reservoirs , 2016 .

[27]  Francisco-Javier Sayas,et al.  A Decoupled Preconditioning Technique for a Mixed Stokes–Darcy Model , 2013, J. Sci. Comput..

[28]  Md. Abdullah Al Mahbub,et al.  Mixed stabilized finite element method for the stationary Stokes-dual-permeability fluid flow model , 2020 .

[29]  Sheik S. Rahman,et al.  Estimation of permeability of naturally fractured reservoirs by pressure transient analysis: An innovative reservoir characterization and flow simulation , 2016 .

[30]  Zhangxin Chen,et al.  A stabilized finite element method based on two local Gauss integrations for a coupled Stokes-Darcy problem , 2016, J. Comput. Appl. Math..

[31]  Roberto Aguilera,et al.  Modeling Fractured Horizontal Wells As Dual Porosity Composite Reservoirs - Application To Tight Gas, Shale Gas And Tight Oil Cases , 2011 .

[32]  Xiaoming He,et al.  A stabilized finite volume element method for a coupled Stokes–Darcy problem , 2017, Applied Numerical Mathematics.

[33]  Jinchao Xu,et al.  A Novel Two-Grid Method for Semilinear Elliptic Equations , 1994, SIAM J. Sci. Comput..

[34]  Alfio Quarteroni,et al.  Robin-Robin Domain Decomposition Methods for the Stokes-Darcy Coupling , 2007, SIAM J. Numer. Anal..

[35]  Dongyang Shi,et al.  Unconditional optimal error estimates of a two-grid method for semilinear parabolic equation , 2017, Appl. Math. Comput..

[36]  Xiaoming He,et al.  Improved Numerical Simulation for Shale Gas Reservoirs , 2014 .

[37]  Yanping Chen,et al.  Two-grid finite element methods combined with Crank-Nicolson scheme for nonlinear Sobolev equations , 2018, Adv. Comput. Math..

[38]  Deepjyoti Goswami,et al.  A two-level finite element method for time-dependent incompressible Navier-Stokes equations with non-smooth initial data , 2012, 1211.3342.

[39]  J. E. Warren,et al.  The Behavior of Naturally Fractured Reservoirs , 1963 .

[40]  M. Abbasi,et al.  Fluid flow in fractured reservoirs: Exact analytical solution for transient dual porosity model with variable rock matrix block size , 2018 .

[41]  R. Raghavan,et al.  Pressure transient analysis methods for bounded naturally fractured reservoirs , 1985 .

[42]  J. Galvis,et al.  NON-MATCHING MORTAR DISCRETIZATION ANALYSIS FOR THE COUPLING STOKES-DARCY EQUATIONS , 2007 .

[43]  A. de Swaan O..  Analytic Solutions for Determining Naturally Fractured Reservoir Properties by Well Testing , 2017 .

[44]  Shuyu Sun,et al.  Coupling nonlinear Stokes and Darcy flow using mortar finite elements , 2011 .

[45]  Horst H. Gerke,et al.  Evaluation of a first-order water transfer term for variably saturated dual-porosity flow models , 1993 .

[46]  Xiaohong Zhu,et al.  Decoupled schemes for a non-stationary mixed Stokes-Darcy model , 2009, Math. Comput..

[47]  Xiaoming He,et al.  On Stokes-Ritz Projection and Multistep Backward Differentiation Schemes in Decoupling the Stokes-Darcy Model , 2016, SIAM J. Numer. Anal..

[48]  Xiaoming He,et al.  A Dual-Porosity-Stokes Model and Finite Element Method for Coupling Dual-Porosity Flow and Free Flow , 2016, SIAM J. Sci. Comput..

[49]  P. Saffman.  On the Boundary Condition at the Surface of a Porous Medium , 1971 .

[50]  Guangzhi Du,et al.  A parallel two-grid linearized method for the coupled Navier-Stokes-Darcy problem , 2017, Numerical Algorithms.

[51]  W. Layton,et al.  A decoupling method with different subdomain time steps for the nonstationary stokes–darcy model , 2013 .

[52]  O. Cicek.  Compositional and Non-Isothermal Simulation of CO2 Sequestration in Naturally Fractured Reservoirs/Coalbeds: Development and Verification of the Model , 2003 .

[53]  Zhangxin Chen,et al.  Discontinuous Finite Volume Element Method for a Coupled Non-stationary Stokes–Darcy Problem , 2018, J. Sci. Comput..

[54]  Chenchen Wang,et al.  Discrete Fracture-Vug Network Model for Modeling Fluid Flow in Fractured Vuggy Porous Media , 2010 .

[55]  Béatrice Rivière,et al.  A Finite Element Method with Strong Mass Conservation for Biot’s Linear Consolidation Model , 2017, J. Sci. Comput..

[56]  Rocky Allen Seale,et al.  An Effective Openhole Horizontal Completion System for Multistage Fracturing and Stimulation , 2008 .

[57]  D. Roehl,et al.  Integrated discrete fracture and dual porosity - Dual permeability models for fluid flow in deformable fractured media , 2019, Journal of Petroleum Science and Engineering.

[58]  Hassan Hassanzadeh,et al.  Matrix-fracture transfer shape factor for modeling flow of a compressible fluid in dual-porosity media , 2011 .

[59]  Nanfeng Sun,et al.  Effects of transmural pressure and wall shear stress on LDL accumulation in the arterial wall: a numerical study using a multilayered model. , 2007, American journal of physiology. Heart and circulatory physiology.

[60]  Xiao-Chuan Cai,et al.  A Highly Scalable Multilevel Schwarz Method with Boundary Geometry Preserving Coarse Spaces for 3D Elasticity Problems on Domains with Complex Geometry , 2016, SIAM J. Scientific Computing.

[61]  Ivan Yotov,et al.  Coupling Fluid Flow with Porous Media Flow , 2002, SIAM J. Numer. Anal..

[62]  Guangzhi Du,et al.  A two‐grid method with backtracking for the mixed Navier–Stokes/Darcy model , 2020, Numerical Methods for Partial Differential Equations.

[63]  Vartan Kurtcuoglu,et al.  Computational modeling of coupled blood-wall mass transport of LDL: effects of local wall shear stress. , 2008, American journal of physiology. Heart and circulatory physiology.

[64]  VIVETTE GIRAULT,et al.  DG Approximation of Coupled Navier-Stokes and Darcy Equations by Beaver-Joseph-Saffman Interface Condition , 2009, SIAM J. Numer. Anal..

[65]  M. Gunzburger,et al.  Coupled Stokes-Darcy model with Beavers-Joseph interface boundary condition , 2010 .

[66]  Eric S. Carlson,et al.  Devonian shale gas production: mechanisms and simple models , 1991 .

[67]  D. Joseph,et al.  Boundary conditions at a naturally permeable wall , 1967, Journal of Fluid Mechanics.

[68]  Y. Rong,et al.  A second-order artificial compression method for the evolutionary Stokes-Darcy system , 2019, Numerical Algorithms.

[69]  Béatrice Rivière,et al.  A strongly conservative finite element method for the coupling of Stokes and Darcy flow , 2010, J. Comput. Phys..

[70]  Ivan Yotov,et al.  Discontinuous Galerkin and mimetic finite difference methods for coupled Stokes–Darcy flows on polygonal and polyhedral grids , 2013, Numerische Mathematik.

[71]  K Perktold,et al.  Mathematical and numerical models for transfer of low-density lipoproteins through the arterial walls: a new methodology for the model set up with applications to the study of disturbed lumenal flow. , 2005, Journal of biomechanics.

[72]  Xiaoming He,et al.  A Domain Decomposition Method for the Steady-State Navier-Stokes-Darcy Model with Beavers-Joseph Interface Condition , 2015, SIAM J. Sci. Comput..

[73]  Lawrence W. Teufel,et al.  Coupling Fluid-Flow and Geomechanics in Dual-Porosity Modeling of Naturally Fractured Reservoirs - Model Description and Comparison , 2000 .

[74]  Yanren Hou,et al.  Optimal error estimates of a decoupled scheme based on two-grid finite element for mixed Stokes-Darcy model , 2015, Appl. Math. Lett..

[75]  Sheik S. Rahman,et al.  Numerical simulation of Fluid-Rock coupling heat transfer in naturally fractured geothermal system , 2011 .

[76]  G. Gatica,et al.  A conforming mixed finite-element method for the coupling of fluid flow with porous media flow , 2008 .

[77]  Vahid Nassehi,et al.  Modelling of combined Navier–Stokes and Darcy flows in crossflow membrane filtration , 1998 .

[78]  Zhangxin Chen,et al.  Effect of fracture pressure depletion regimes on the dual-porosity shape factor for flow of compressible fluids in fractured porous media , 2011 .

[79]  Li Shan,et al.  Optimal error estimates of both coupled and two-grid decoupled methods for a mixed Stokes–Stokes model , 2018, Applied Numerical Mathematics.

[80]  Yunqing Huang,et al.  A modified two-grid decoupling method for the mixed Navier-Stokes/Darcy Model , 2016, Comput. Math. Appl..

[81]  Xiaoming He,et al.  Multi-Stage Fractured Horizontal Well Numerical Simulation and its Application in Tight Shale Reservoirs , 2015 .

[82]  Hoang Tran,et al.  Analysis of Long Time Stability and Errors of Two Partitioned Methods for Uncoupling Evolutionary Groundwater-Surface Water Flows , 2013, SIAM J. Numer. Anal..

[83]  Iraj Ershaghi,et al.  Pressure Transient Analysis of Dually Fractured Reservoirs , 1996 .

[84]  Yunqing Huang,et al.  Two-grid methods of expanded mixed finite-element solutions for nonlinear parabolic problems , 2019, Applied Numerical Mathematics.

[85]  Roberto Aguilera,et al.  Naturally Fractured Reservoirs , 1980 .

[86]  Gabriel N. Gatica,et al.  A Residual-Based A Posteriori Error Estimator for the Stokes-Darcy Coupled Problem , 2010, SIAM J. Numer. Anal..

[87]  Svetlana Tlupova,et al.  Domain Decomposition Methods for Solving Stokes-Darcy Problems with Boundary Integrals , 2013, SIAM J. Sci. Comput..

[88]  Benoit Noetinger,et al.  Scaling Up Matrix-Fracture Transfers in Dual-Porosity Models: Theory and Application , 1999 .

[89]  Mingchao Cai,et al.  A multilevel decoupled method for a mixed Stokes/Darcy model , 2012, J. Comput. Appl. Math..

[90]  K. Aziz,et al.  Matrix-fracture transfer shape factors for dual-porosity simulators , 1995 .

[91]  Xiaoming He,et al.  A Coupled Multiphysics Model and a Decoupled Stabilized Finite Element Method for the Closed-Loop Geothermal System , 2020, SIAM J. Sci. Comput..

[92]  Robert A. Wattenbarger,et al.  Rate Transient Analysis in Naturally Fractured Shale Gas Reservoirs , 2008 .

[93]  Xiaoming He,et al.  Coupled and decoupled stabilized mixed finite element methods for nonstationary dual‐porosity‐Stokes fluid flow model , 2019, International Journal for Numerical Methods in Engineering.