A conservative sequential fully implicit method for compositional reservoir simulation

Lee, S.H.; Ţene, M.; Du, Shichao; Wen, Xian‐Huan; Efendiev, Yalchin

doi:10.1016/j.jcp.2020.109961

Cited by 9 publications

(1 citation statement)

References 36 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The Sequential Fully-Implicit (SFI) scheme, developed in the Multiscale Finite Volume Method (MSFV) framework [20], is an option to deal with this issue. See [32,33,34,35] for extensions to compositional flow simulations, and [36] for developments on nonlinear acceleration techniques. We also study in this work the MRCM combined with the SFI scheme to approximate a challenging two-phase flow problem, with strong fingering instabilities [37].…”

Section: Introductionmentioning

confidence: 99%

A multiscale Robin-coupled implicit method for two-phase flows in high-contrast formations

Rocha¹,

Sousa²,

Ausas³

et al. 2021

Preprint

View full text Add to dashboard Cite

In the presence of strong heterogeneities, it is well known that the use of explicit schemes for the transport of species in a porous medium suffers from severe restrictions on the time step. This has led to the development of implicit schemes that are increasingly favoured by practitioners for their computational efficiency. The transport equation requires knowledge of the velocity field, which results from an elliptic problem (Darcy problem) that is the most expensive part of the computation. When considering large reservoirs, a cost-effective way of approximating the Darcy problems is using multiscale domain decomposition (MDD) methods. They allow for the pressure and velocity fields to be computed on coarse meshes (large scale), while detailed basis functions are defined locally, usually in parallel, in a much finer grid (small scale). In this work we adopt the Multiscale Robin Coupled Method (MRCM, [Guiraldello, et al., J. Comput. Phys., 355 (2018) pp. 1-21], [Rocha, et al., J. Comput. Phys., (2020) 109316]), which is a generalization of previous MDD methods that allows for great flexibility in the choice of interface spaces. In this article we investigate the combination of the MRCM with implicit transport schemes. A sequentially implicit strategy is proposed, with different trust-region algorithms ensuring the convergence of the transport solver. The method is assessed on several very stringent 2D two-phase problems, demonstrating its stability even for large time steps. It is also shown that the best accuracy is achieved by considering recently introduced non-polynomial interface spaces, since polynomial spaces are not optimal for high-contrast channelized permeability fields.

show abstract