A novel method to accelerate the three-phase flash calculation of the hydrocarbon-CO 2 system

To achieve high computational efficiency and simulation accuracy in the equation-of-state-based compositional simulation, the fast multiphase flash calculation approach is required for the hydrocarbon-CO 2 system, allowing for the simultaneous existence of vapour, CO 2-lean, and CO 2-rich phases. However, the conventional step-by-step method for multiphase flash calculation is associated with a high cost in identifying the three-phase equilibrium, thereby significantly impeding the computation efficiency of compositional simulation when the three-phase equilibrium appears in most grids. To address this challenge, a novel approach to multiphase flash calculation is proposed. In this study, all trial phases exhibiting negative tangent plane distance are employed to identify the three-phase equilibrium, eliminating the need for additional phase stability analysis and two-phase split calculations. Furthermore, this new method incurs minimal computational cost and requires only minor modifications to the existing stepwise procedure. Phase diagram calculations conducted on four cases demonstrate that the newly proposed method dramatically accelerates the identification of the three-phase equilibrium.