Two-grid method on unstructured tetrahedra: Applying computational geometry to staggered solution of coupled flow and mechanics problems

Saumik Dana,Xiaoxi Zhao,Birendra Jha

DOI: https://doi.org/10.48550/arXiv.2102.04455

2021-02-08

Abstract:We develop a computational framework that leverages the features of sophisticated software tools and numerics to tackle some of the pressing issues in the realm of earth sciences. The algorithms to handle the physics of multiphase flow, concomitant geomechanics all the way to the surface of the earth and the complex geometries of field cases with surfaces of discontinuity are stacked on top of each other in a modular fashion which allows for easy use to the end user. The current focus of the framework is to provide the user with tools for assessing seismic risks associated with energy technologies as well as for use in generating forward simulations in inversion analysis from data obtained using GPS and InSAR. In this work, we focus on one critical aspect in the development of the framework: the use of computational geometry in a two-grid method for unstructured tetrahedral meshes

Computational Engineering, Finance, and Science

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop a computationally inexpensive framework that can simulate surface deformation and fault slip in deep overburden, which are related to multiphase flow and associated geomechanics. Specifically, the framework aims to evaluate the seismic risks associated with energy technologies such as carbon capture, enhanced geothermal systems, wastewater treatment, and enhanced oil recovery, and is used to generate forward simulations required for inversion analysis from data obtained from the Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar (InSAR).

The main challenge is that regardless of which numerical method is used to solve the coupled equations, the computational scale of the geomechanical grid is difficult to handle. Existing methods either directly impose overburden pressure on the reservoir, treating it as the coupled problem domain, which excludes the study of surface deformation, cannot truly reflect the influence of the overburden on pressure - sensitive reservoirs, and cannot capture induced earthquakes within faults far from the reservoir; or model the flow on a large domain and simulate the no - flow boundary conditions between the flowing and non - flowing regions through zero - permeability elements, but this is computationally infeasible for large - scale field problems due to memory requirements.

To address this challenge, the authors developed a two - grid - coupled multiphase - flow and geomechanical framework that allows for the spatial decoupling of the flow and geomechanical domains, in which the geomechanical sub - problem is solved on a separate grid with a larger spatial extent, extending all the way to the free surface. This method can not only study induced earthquakes within faults far from the reservoir, but also consider the mechanical behavior of the surrounding rocks and allow for the determination of surface deformation caused by underground pressure perturbations.