A Scalable Numerical Solution Method for Fine Scale Compositional Simulation

Summary Nonlinear iteration is a very important part of reservoir numerical simulation. If the iteration steps can be reduced, the speed of the iteration and the reservoir numerical simulation will be helped greatly. In this study, a self-equilibrium Newton method has been explored for nonlinear system of compositional simulator, which is composed of inexact Newton method, non-negative constraint line search method of molar concentration and automatic time step selection method. Furthermore, a multilevel preconditioning operator and domain decomposition method are adopted to realize parallel computation of the simulation. Comparison with traditional methods, numerical tests indicate that the self-equilibrium Newton method can improve the overall efficiency of the nonlinear system of the component model.