An adaptive two-grid preconditioner and linearly implicit scheme for shale gas transport in fractured porous media

TL;DR

This paper develops an adaptive two-grid preconditioner combined with a linearly implicit scheme to efficiently simulate nonlinear gas transport in fractured shale formations, capturing complex fracture-matrix interactions.

Contribution

It introduces an adaptive multiscale space based on local spectral problems for improved two-grid preconditioning in fractured media simulations.

Findings

Iterative convergence independent of fracture-matrix permeability contrast

Effective capture of local fracture features through spectral space

Demonstrated scalability on models with up to 160 fractures

Abstract

We consider a nonlinear mixed-dimensional model for simulating gas transport in shale formation. The mathematical model consists of a coupled system of nonlinear equations, where flow within fractures is represented using a lower-dimensional representation. For the numerical solution of the coupled transport problem, we construct an unstructured mesh that resolves lower dimensional fractures on the grid level and use the finite element approximation to build a discrete system. To construct an efficient scheme for the resulting nonlinear problem, we use an explicit-implicit method for time integration, where we carefully choose an additive partition of the nonlinear operators to separate the stiff linear component and integrate it implicitly to ensure the stability of the time integration. Next, we invert the linear partition of the operator by constructing an efficient two-grid preconditioner for shale gas transport in fractured porous media. We use a local pointwise smoother on the fine grid and carefully design an adaptive multiscale space for coarse grid approximation based on local generalized eigenvalue problems. We utilize an adaptive thresholding to automatically identify local dominant modes which correspond to the very small eigenvalues in local domains. We remark that such spatial features are automatically captured through our local spectral problems, and connect these to fracture information in the global formulation of the problem. Approximation properties of the local spectral space with convergence of the proposed two-grid algorithm are given. Numerical results are presented for two fracture distributions with 30 and 160 fractures, demonstrating iterative convergence independent of the contrast of fracture and porous matrix permeability.