No abstract

No abstract

Numerical models of geologic carbon sequestration (GCS) in saline aquifers use multiphase fluid flow-characteristic curves (relative permeability and capillary pressure) to represent the interactions of the non-wetting CO2 and the wetting brine. Relative permeability data for many sedimentary formations is very scarce, resulting in the utilisation of mathematical correlations to generate the fluid flow characteristics in these formations. The flow models are essential for the prediction of CO2 storage capacity and trapping mechanisms in the geological media. The observation of pressure dissipation across the storage and sealing formations is relevant for storage capacity and geomechanical analysis during CO2 injection. This paper evaluates the relevance of representing relative permeability variations in the sealing formation when modelling geological CO2 sequestration processes. Here we concentrate on gradational changes in the lower part of the caprock, particularly how they affect pressure evolution within the entire sealing formation when duly represented by relative permeability functions. The results demonstrate the importance of accounting for pore size variations in the mathematical model adopted to generate the characteristic curves for GCS analysis. Gradational changes at the base of the caprock influence the magnitude of pressure that propagates vertically into the caprock from the aquifer, especially at the critical zone (i.e. the region overlying the CO2 plume accumulating at the reservoir-seal interface). A higher degree of overpressure and CO2 storage capacity was observed at the base of caprocks that showed gradation. These results illustrate the need to obtain reliable relative permeability functions for GCS, beyond just permeability and porosity data. The study provides a formative principle for geomechanical simulations that study the possibility of pressure-induced caprock failure during CO2 sequestration.

This document is the author's post-print version, incorporating any revisions agreed during the peer-review process. Some differences between the published version and this version may remain and you are advised to consult the published version if you wish to cite from it.

In reservoir engineering, the predictive analyses of CO2 sequestration in subsurface formations commonly employ numerical models of subsurface formations. A significant number of work have utilised numerical modelling techniques to predict the impact of the reservoir's boundary conditions and interlayer communication on CO2 storage capacity in aquifers. To the best of our knowledge, no study on the impact of boundary conditions on CO2 storage efficiency has focused on the combined effect of this factor in the reservoir and saturation functions in the caprock. To this end, this study examined the effect of integrating both processes on pressure evolution in the caprock during the numerical simulation of CO2 injection into a deep saline aquifer. Utilising the Sleipner benchmark model, we also showed how varying saturation functions in the caprock can affect the storage efficiency in the reservoir formation.