Dissipation-Based Continuation Method for Multiphase Flow in Heterogeneous Porous Media

TL;DR

This paper introduces a dissipation-based continuation method to improve the convergence of nonlinear solvers in reservoir simulations involving multiphase flow in heterogeneous porous media, demonstrating enhanced efficiency over standard methods.

Contribution

A novel dissipation-based continuation approach is developed to enhance nonlinear solver convergence in complex multiphase flow simulations, with adaptive dissipation control and demonstrated effectiveness.

Findings

Improved convergence in 1D and 2D heterogeneous flow problems.

Adaptive dissipation strategy enhances solver robustness.

Outperforms standard damped-Newton methods in tests.

Abstract

In reservoir simulation, solution of the coupled systems of nonlinear algebraic equations that are associated with fully-implicit (backward Euler) discretization is challenging. Having a robust and efficient nonlinear solver is necessary in order for reservoir simulation to serve as the primary tool for managing the recovery processes of large-scale reservoirs. Here, we develop a continuation method based on the use of a dissipation operator. We focus on nonlinear two-phase flow and transport in heterogeneous formations in the presence of viscous, gravitational, and capillary forces. The homotopy is constructed by adding numerical dissipation to the coupled discrete conservation equations. A continuation parameter is introduced to control the amount of dissipation. Numerical evidence of multi-dimensional models and detailed analysis of single-cell problems are used to explain how the dissipation operator improves the nonlinear convergence of the coupled system of equations. An adaptive strategy to determine the dissipation coefficient is proposed. The dissipation level is computed locally for each cell interface. We demonstrate the efficiency of the dissipation-based continuation (DBC) nonlinear solver using several examples, including 1D scalar transport and 2D heterogeneous problems with fully-coupled flow and transport. The DBC solver has better convergence properties compared with the standard damped-Newton solvers used in reservoir simulation.