Improved phase-state identification bypass approach of the hydrocarbons-CO2-H2O system for compositional reservoir simulation

Abstract

CO2 injection is a highly effective technique to enhance oil recovery, achieved through continuous or alternative injection. However, the intricate interactions between different phases within porous media present significant challenges when predicting the performance of CO2 injection. To address this, it is crucial to employ compositional simulation, which accounts for the multiphase multicomponent transport. Nonetheless, conventional multiphase flash calculations can be computationally inefficient for large-scale reservoir simulations. Therefore, it is necessary to accelerate the Equation-of-State (EoS)-based compositional simulation, given the widespread use of CO2 enhanced oil recovery (CO2-EOR) in recent years. The phase-state identification bypass method has proven to be superior to other methods in terms of efficiency. However, this approach struggles with regions near phase boundaries, resulting in reduced computational efficiency in those areas.In this study, an enhanced phase-state identification bypass approach is developed to address this limitation. The first step involves discretising the pressure-temperature space using rectangular grids. Additionally, the tie-simplexes, which represent regions defined by the maximum number of phases formed by the fluid under consideration, are discretized in the phase-fraction space at the pressure and temperature of each discretization node. Subsequently, the discretization grid associated with the given point (the overall composition, pressure, and temperature) is located, and the phase states of the grid nodes are determined using the conventional multiphase flash method. If all nodes exhibit the same phase state, that phase state is assigned to the given point. However, if multiple phase states are obtained, a novel process is proposed to determine the phase state of the given point. To validate this improvement to the phase-state identification bypass method, phase diagram calculations and simulation cases are conducted, and the results demonstrate the robustness of the proposed method and its superior computational efficiency compared to the previous method.