Geological carbon storage represents a new and substantial challenge for the subsurface geosciences. To increase understanding and make good engineering decisions, containment processes and large-scale storage operations must be simulated in a thousand year perspective. A hierarchy of models of increasing computational complexity for analysis and simulation of large-scale CO₂ storage has been implemented as a separate module of the open-source Matlab Reservoir Simulation Toolbox (MRST). This paper describes a general family of two-scale models available in this module. The models consist of two-dimensional flow equations formulated in terms of effective quantities obtained from hydrostatic reconstructions of vertical pressure and saturation distributions. The corresponding formulation is fully implicit and is the first to give a mass-conservative treatment and include general (non-linearized) CO₂ properties. In particular, the models account for compressibility, dissolution, and hysteresis effects in the fine-scale capillary and relative permeability functions and can be used to accurately and efficiently study the combined large-scale and long-term effects of structural, residual, and solubility trapping.

Modeling geological carbon storage represents a new and substantial challenge for the subsurface geosciences. To increase understanding and make good engineering decisions, containment processes and large-scale storage operations must be simulated in a thousand-year perspective. Large differences in spatial and temporal scales make it prohibitively expensive to compute the fate of injected CO₂ using traditional 3D simulators. Instead, accurate forecast can be computed using simplified models that are adapted to the specific setting of the bouyancy-driven migration of the light fluid phase. This paper presents a family of vertically integrated models for studying the combined large-scale and long-term effects of structural, residual, and solubility trapping of CO₂. The models are based on an assumption of a sharp interface separating CO₂ and brine and can provide a detailed inventory of the injected CO₂ volumes over periods of thousands of years within reasonable computational time. To be compatible with simulation tools used in industry, the models are formulated in a black-oil framework. The models are implemented in MRST-co2lab, which is an open community software developed especially to study and optimize large-scale, long-term geological storage of CO₂. The resulting simulators are fully implicit and handle input from standard geomodeling tools.