In this work, we present a unified framework for the simulation of CO₂ sequestration problems at various time and space scales. The parametrization technique utilizes thermodynamic state-dependent operators expressing the governing equations for the thermal-compositional-reactive system to solve the nonlinear problem. This approach provides flexibility in the assembly of the Jacobian, which allows straightforward implementation of advanced thermodynamics. We validate our simulation framework through several simulation studies including complex physical phenomena relevant to CCUS. The proposed simulation framework is validated against a set of numerical and experimental benchmark tests, demonstrating the efficiency and accuracy of the modeling framework for CCUS-related subsurface applications. Important physical phenomena resulting from the complex thermodynamic interactions of CO₂ and impurities with reservoir fluids can be accurately captured now in detailed dynamic simulation. The investigated simulation scenarios include a reproduction of lab experiments at the core scale, investigation of macro-scale analog model and simulation of large-scale industrial application. The simulation time can also span from hours to years among various applications. Complex thermal-compositional-reactive phenomena can be addressed at each of these space and time scales. The unified thermodynamic description allows us to perform all these simulations for a reasonable CPU time due to advanced parametrization techniques and efficient GPU capabilities in our in-house reservoir simulator DARTS.

Several compositional reservoir simulators based on equations of state (EoS) have been designed. Yet, few of them can deal with thermal processes such as steam injection and no fully compositional thermal simulation of the steam injection process has been proposed with extra-heavy oils yet. Those simulations appear essential and will offer better tools to decide whether to carry out the exploitation of a heavy oil field or not. In those processes, the accurate modelling of the water/steam phase plays an important role and accurate multiphase equilibrium calculations are necessary. Thermodynamics generate highly nonlinear problems. Besides, in reservoir simulation a huge amount of phase equilibrium calculations is required, and a single failure may cause significant error propagations leading to false solutions. This study presents several improvements leading to a more robust and efficient phase equilibrium calculation (stability and flash) program. A general multiphase flash implementing all the developed algorithms is presented and tested against experimental and literature data. It can handle any number of phases; a numerical example consists in a four-phase simulation of CO₂ injection. Simulations of steam flooding are performed with highly heterogeneous reservoirs. Besides, a fully compositional simulation of the SAGD process for an extra heavy bitumen is performed, which appears to be the one of the first simulation of the kind in the literature.