TY - JOUR
T1 - Multiscale finite volume method for finite-volume-based simulation of poroelasticity
AU - Sokolova, Irina
AU - Bastisya, Muhammad Gusti
AU - Hajibeygi, Hadi
PY - 2019
Y1 - 2019
N2 - We propose a multiscale finite volume method (MSFV) for simulation of coupled flow-deformation in heterogeneous porous media under elastic deformation (i.e., poroelastic model). The fine-scale fully resolved system of equations is obtained based on a conservative finite-volume method in which the displacement and pore pressure unknowns are located in a staggered configuration. The coupling is treated through a fully-coupled fully-implicit formulation. On this fully-coupled finite-volume system, coarse-scale grids for flow and deformation are imposed. Local basis functions for scalar pore pressure and vectorial displacement unknowns are then solved over their respective local domains at the beginning of the simulation, and reused for the rest of the time-dependent simulations. These local basis functions are then clustered to form the prolongation operator. As for the finite-volume nature of the proposed multiscale system, finite-volume restriction operators for poroelastic systems are utilised. Once the coarse-scale system is solved, its solution is prolonged back to the original fine-scale resolution, providing approximate fine-scale solution. The finite-volume multiscale formulation provides conservative stress and mass flux both at fine and coarse scale. Several numerical test cases are provided first to validate the fine-scale finite-volume discrete fully-implicit simulation, and then to investigate the accuracy of the proposed multiscale formulation. Moreover we compare our fully implicit MSFV method with hybrid multiscale Finite Element-Finite Volume (h-MSFE-FV). Our multiscale method allows for quantification of the elastic geomechanical behaviour with using only a fraction of the fine-scale grid cells, even for highly heterogeneous time-dependent models. As such, it casts a promising approach for field-scale quantification of the mechanical deformation and stress field due to injection and production in a subsurface formation.
AB - We propose a multiscale finite volume method (MSFV) for simulation of coupled flow-deformation in heterogeneous porous media under elastic deformation (i.e., poroelastic model). The fine-scale fully resolved system of equations is obtained based on a conservative finite-volume method in which the displacement and pore pressure unknowns are located in a staggered configuration. The coupling is treated through a fully-coupled fully-implicit formulation. On this fully-coupled finite-volume system, coarse-scale grids for flow and deformation are imposed. Local basis functions for scalar pore pressure and vectorial displacement unknowns are then solved over their respective local domains at the beginning of the simulation, and reused for the rest of the time-dependent simulations. These local basis functions are then clustered to form the prolongation operator. As for the finite-volume nature of the proposed multiscale system, finite-volume restriction operators for poroelastic systems are utilised. Once the coarse-scale system is solved, its solution is prolonged back to the original fine-scale resolution, providing approximate fine-scale solution. The finite-volume multiscale formulation provides conservative stress and mass flux both at fine and coarse scale. Several numerical test cases are provided first to validate the fine-scale finite-volume discrete fully-implicit simulation, and then to investigate the accuracy of the proposed multiscale formulation. Moreover we compare our fully implicit MSFV method with hybrid multiscale Finite Element-Finite Volume (h-MSFE-FV). Our multiscale method allows for quantification of the elastic geomechanical behaviour with using only a fraction of the fine-scale grid cells, even for highly heterogeneous time-dependent models. As such, it casts a promising approach for field-scale quantification of the mechanical deformation and stress field due to injection and production in a subsurface formation.
KW - Algebraic multiscale solver
KW - Finite volume method
KW - Geomechanics
KW - Multiscale finite-volume method
KW - Poromechanics
KW - Porous media
UR - http://www.scopus.com/inward/record.url?scp=85058932508&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2018.11.039
DO - 10.1016/j.jcp.2018.11.039
M3 - Article
AN - SCOPUS:85058932508
SN - 0021-9991
VL - 379
SP - 309
EP - 324
JO - Journal of Computational Physics
JF - Journal of Computational Physics
ER -