Nonlinear solver based on trust region approximation for <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" altimg="si129.svg" display="inline" id="d1e590"><mml:msub><mml:mrow><mml:mi mathvariant="normal">CO</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:math> utilization and storage in subsurface reservoir

Abstract

Simulation of CO 2 utilization and storage (CCUS) in subsurface reservoirs with complex heterogeneous structures requires a model that captures multiphase compositional flow and transport. Accurate simulation of these processes necessitates the use of stable numerical methods that are based on an implicit treatment of the flux term in the conservation equation. Due to the complicated thermodynamic phase behavior, including the appearance and disappearance of multiple phases, the discrete approximation of the governing equations is highly nonlinear. Consequently, robust and efficient techniques are needed to solve the resulting nonlinear system of algebraic equations. In this study, we present a powerful nonlinear solver based on a generalization of the trust-region technique for compositional multiphase flows. The approach is designed to embed a newly introduced Operator-Based Linearization technique and is grounded on the analysis of multi-dimensional tables related to parameterized convection operators. We split the parameter space of the nonlinear problem into a set of trust regions where the convection operators preserve the second-order behavior (i.e., they remain positive or negative definite). We approximate these trust regions in the solution process by detecting the boundary of convex regions via analysis of the directional derivative. This analysis is performed adaptively while tracking the nonlinear update trajectory in the parameter space. The proposed nonlinear solver locally constrains the update of the overall compositions across the boundaries of convex regions. We tested the performance of the proposed nonlinear solver for various scenarios. In many cases, our approach yields an improved behavior of the nonlinear solution in comparison to state-of-the-art solvers. • We study the nonlinearity of compositional flow and transport applied for CCUS simulation. • We utilize operator-based linearization to facilitate nonlinear analysis of the equation. • There are inflection/kink lines in highly nonlinear convective operators that can lead to nonlinear convergence issues. • Our proposed algorithm approximates TR boundaries accurately, reducing detection cost compared to Hessian assembly. • Our TR solver enables larger timesteps and achieves better convergence compared to existing methods.

Topics & Concepts