The secure and efficient storage of carbon dioxide (CO2) in deep saline aquifers is widely recognized as a critical component in the global strategy for mitigating anthropogenic climate change. Accurate and computationally efficient modeling of CO2
migration, phase behavior, and long-term trapping mechanisms in geologically heterogeneous formations remains an open and pressing challenge. This dissertation aims to advance modeling techniques for multiphase multicomponent flow in porous media,
with a focus on two complementary directions: physics-constrained deep learning and multiscale numerical solvers. Both approaches seek to provide computational efficiency while delivering accurate results on a desired level, thereby enabling predictive
simulations across a range of spatial and temporal scales.

This study introduces a multiscale simulation framework, termed Projection-based Embedded Discrete Fracture Modeling with Algebraic Dynamic Multilevel method (pEDFM-ADM), which integrates an embedded discrete fracture network representation with a fully algebraic, front-tracking-based mesh adaptation strategy. Incorporating a fully implicit scheme, compositional thermodynamics, and algebraic multilevel operators, the framework captures essential subsurface processes such as buoyancy-driven migration, convective dissolution, phase partitioning, and fracture-matrix interactions under geologically realistic conditions. The method constructs a hierarchy of multilevel grids and localized multiscale basis functions that introduce fine-scale heterogeneities at each coarse level. Adaptive mesh refinement and coarsening are driven by local variations in CO₂ mass fraction and executed through algebraic prolongation and restriction operators, enabling efficient projection between grid levels. The framework is systematically evaluated across a sequence of test cases with increasing complexity, including systems with low-permeability flow barriers, highly conductive fractures, striking a trade-off between computational resource and detailed simulation accuracy. Overall, the pEDFM-ADM framework provides a scalable, fully algebraic, and physically adaptive modeling tool for large-scale CO₂ storage simulations in fractured porous media, supporting predictive simulation and risk assessment for long-term carbon sequestration.

This work introduces a novel application of the Algebraic Dynamic Multilevel (ADM) method for simulating CO₂ storage in deep saline aquifers. By integrating a fully implicit coupling strategy, fully compositional thermodynamics, and adaptive mesh refinement, the ADM framework effectively models phenomena such as buoyancy-driven migration, convective dissolution, and phase partitioning under various subsurface conditions. The method starts with the construction of a hierarchy of multilevel grids and the generation of localized multiscale basis functions, which account for heterogeneities at each coarse level. During the simulation, the ADM method dynamically refines areas with significant overall CO₂ mass fraction gradients while coarsening smooth regions, thus optimizing computational resources without compromising the accuracy required to capture essential flow and transport characteristics. This dynamic grid adjustment is facilitated by algebraic prolongation and restriction operators, which map the fine-scale system onto a coarser grid suited to the evolving distribution of the CO₂ plume. This feature allows the ADM to navigate various coarsening thresholds efficiently, striking a trade-off between computational economy and detailed simulation accuracy. Even at relatively higher thresholds, key trapping mechanisms are captured with sufficient detail for quantification. These capabilities make the ADM framework well suited for long-term CO₂ sequestration in highly heterogeneous reservoirs, where large-scale models may otherwise become impractically expensive, offering a practical balance between the need for detailed simulations and manageable computational requirements. Overall, the ADM framework proves to be a robust tool for designing, monitoring, and analyzing large-scale CO₂ storage operations, supporting reliable and cost-effective solutions in carbon management.

CO₂ capture and storage is a viable solution in the effort to mitigate global climate change. Deep saline aquifers, in particular, have emerged as promising storage options, owing to their vast capacity and widespread distribution. However, the task of proficiently monitoring and simulating CO₂ behavior within these formations poses significant challenges. To address this, we introduce the physics-constraint neural network for CO₂ storage (CO₂PCNet), a model specifically designed for simulating and monitoring CO₂ storage in deep saline aquifers during injection and post-injection periods. Recognizing the significant challenges in accurately modeling the distribution and movement of CO₂ under varying permeability conditions, the CO₂PCNet integrates the principles of physics with the robustness of deep learning, serving as a powerful surrogate model. The architecture of CO₂PCNet starts with an encoder that adeptly processes spatial features from overall mole fraction (z_CO2) and pressure fields (P_l), capturing the complex dynamics of a CO₂ trajectory. By incorporating permeability information through a conditioning step, the network ensures a faithful representation of the influences on CO₂ behavior in subsurface conditions. A ConvLSTM module subsequently discerns temporal evolutions, reflecting the real-world progression of CO₂ plumes within the reservoir. Lastly, the decoder precisely reconstructs the predictive spatial profile of CO₂ distribution. CO₂PCNet, with its integration of convolutional layers, recurrent mechanisms, and physics-informed constraints, offers a refined approach to CO₂ storage simulation. This model offers the potential of utilizing advanced computational methods in advancing CCS practices.

CO₂ sequestration and storage in deep saline aquifers is a promising technology for mitigating the excessive concentration of the greenhouse gas in the atmosphere. However, accurately predicting the migration of CO₂ plumes requires complex multi-physics-based numerical simulation approaches, which are prohibitively expensive due to highly nonlinear coupled governing equations and uncertainties in heterogeneous spatial parameter distributions. To address this challenge, we developed an end-to-end deep learning workflow employing encoder–decoder architectures with residual network (ResNet) to efficiently predicts the spatial–temporal evolution of the solution CO₂-brine ratio (R_s) and gas saturation (S_g) – the two essential tasks for quantifying the amount of trapped CO₂ – given heterogeneous permeability fields as input. Specifically, we introduce a general multi-task learning with consistency (MTLC) framework to simultaneously predict R_s and S_g. The MTLC model leverages related tasks with less computational expensive labeled datasets to improve generalization ability. In our study, predictions for multiple tasks from the same permeability realization are not independent but expected to be consistent, as the proposed framework utilizes data-driven cross-task consistency constraints to augment learning of related tasks. Our deep learning model is trained based on physical trapping mechanisms, which play a dominant role in the CO₂ migration process. The results demonstrate that the MTLC model with joint learning yields more accurate predictions and improved generalization for predicting CO₂ migration in several test cases. Furthermore, our workflow is 10⁵ times faster than a high-fidelity physics-based numerical simulator, making it a viable alternative for field-scale applications.