Operator splitting technique using streamline projection for two-phase flow in highly heterogeneous and anisotropic porous media
This work addresses the need for efficient simulation of two-phase flow in highly heterogeneous porous media, which is relevant for applications like oil reservoir engineering.
The authors developed a fast streamline-based numerical method for two-phase flow in heterogeneous and anisotropic porous media, combining a discontinuous Galerkin method for pressure with a one-dimensional front tracking method for transport. The method was validated against standard fully-implicit finite volume methods in 2D and 3D tests, showing accurate flow front propagation.
In this paper, we present a fast streamline-based numerical method for the two-phase flow equations in high-rate flooding scenarios for incompressible fluids in heterogeneous and anisotropic porous media. A fractional flow formulation is adopted and a discontinuous Galerkin method (DG) is employed to solve the pressure equation. Capillary effects can be neglected in high-rate flooding scenarios. This allows us to present an improved streamline approach in combination with the one-dimensional front tracking method to solve the transport equation. To handle the high computational costs of the DG approximation, domain decomposition is applied combined with an algebraic multigrid preconditioner to solve the linear system. Special care at the interior interfaces is required and the streamline tracer has to include a dynamic communication strategy. The method is validated in various two- and three-dimensional tests, where comparisons of the solutions in terms of approximation of flow front propagation with standard fully-implicit finite volume methods are provided.