Optimized phase stability algorithm in three phase systems with an application to CO2 capture and sequestration in subsurface reservoirs

In optimizing CO2 capture and sequestration, accurate phase behavior prediction in subsurface reservoirs is critical. This is especially relevant for geologic CO2 storage in deep saline aquifers and depleted oil/gas reservoirs, and for CO2-EOR, where injected CO2 coexists with brine and resident hydrocarbons across wide pressure–temperature ranges (including near-critical CO2) and with non-negligible mutual solubilities. Our work introduces an efficient algorithm for phase stability testing, combined with multistage equilibrium multiphase negative flash calculations. We explore the stability of a three-phase system (V-L-Aq), using the criterion that the system is stable if the phase fractions of vapor, liquid, and water phases (βV,βL,βW) fall between 0 and 1. The stability testing begins with a negative flash calculation using diverse initial K-values to detect multiple phases. The same criterion is applied to two-phase systems (V-L, V-Aq, L-Aq). In case of instability, the algorithm shifts to the most probable two-phase system by systematically using negative flash and tangent plane distance calculations. It determines stable phases based on tangent plane solutions when the negative flash fails to converge, demonstrating adaptability across multiphase systems. The algorithm is validated through simulations of various mixtures, showing reliable convergence near critical points, and is applied to CO2 storage and EOR scenarios to predict different CO2 trapping mechanisms in the subsurface.