A computationally efficient numerical model to understand potential CO2 leakage risk within gigatonne scale geologic storage

Abstract

<p>The majority of available climate change mitigation pathways, targeting net-zero CO<sub>2</sub> emissions by 2050, rely heavily on the permanent storage of CO<sub>2</sub> in deep geologic formations at the gigatonne scale. The spatial and temporal scales of interest to geologic carbon storage (GCS) raise concerns about CO<sub>2</sub> leakage to shallow sediments or back into the atmosphere. The assessment of CO<sub>2</sub> storage performance is subject to huge computational costs of numerically simulating CO<sub>2</sub> migration across geologic layers at the basin scale and is therefore restricted in practice to multi-century periods. Here, we present a computationally affordable and yet physically sound model to understand the likelihood of CO<sub>2</sub> leakage over geologic time scales (millions of years) (Kivi et al., 2022). The model accounts for vertical two-phase flow and transport of CO<sub>2</sub> and brine in a multi-layered system, comprising a sequence of aquifers and sealing rocks from the crystalline basement up to the surface (a total thickness of 1600 m), representative of a sedimentary basin. We argue that the model is capable of capturing the dynamics of CO<sub>2</sub> leakage during basin-wide storage because the lateral advancement of CO<sub>2</sub> plume injected from a dense grid of wellbores transforms into buoyant vertical rise within a short period after shut-in. A critical step in the proposed model is its initialization, which should reproduce the average CO<sub>2</sub> saturation column and pressure profiles. We initialize the model by injecting CO<sub>2</sub> at a constant overpressure into an upper lateral portion of the target aquifer while the bottom boundary is permeable to brine, resembling brine displacement by CO<sub>2</sub> plume or leakage at basin margins. The optimum model setting can be achieved by adjusting the brine leakage parameter through calibration. We solve the governing equations using the finite element code CODE_BRIGHT. Discretizing the model with 7,100 quadrilateral elements and using an adaptive time-stepping scheme, the CPU time for the simulation of CO<sub>2</sub> containment in the subsurface for 1 million years is around 140 hours on a Xeon CPU of speed 2.5 GHz. The obtained results suggest that the upward CO<sub>2</sub> flow in free phase is strongly hindered by the sequence of caprocks even if they are pervasively fractured. CO<sub>2</sub> leakage towards the surface is governed by the intrinsically slow molecular diffusion process, featuring aqueous CO<sub>2</sub> transport rates as low as 1 meter per several thousand years. The model shows that GCS in multi-layered geologic settings is extremely unlikely to be associated with leakage, implying that GCS is a secure carbon removal technology. </p> <p> </p> <p><strong>Reference</strong></p> <p>Kivi, I. R., Makhnenko, R. Y., Oldenburg, C. M., Rutqvist, J., & Vilarrasa, V. (2022). Multi-layered systems for permanent geologic storage of CO<sub>2</sub> at the gigatonne scale. <em>Geophysical Research Letters</em>, 49, e2022GL100443. https://doi.org/10.1029/2022GL100443</p>