Pore Scale Numerical Modelling of Geological Carbon Storage Through Mineral Trapping Using True Pore Geometries

Mineral trapping (MT)is the most secure method of sequestering carbon for geologically significant periods of time. The processes behind MT fundamentally occur at the pore scale, therefore understanding which factors control MT at this scale is crucial. We present a finite elements advection–diffusion–reaction numerical model which uses true pore geometry model domains generated from μCT imaging. Using this model, we investigate the impact of pore geometry features such as branching, tortuosity and throat radii on the distribution and occurrence of carbonate precipitation in different pore networks over 2000 year simulated periods. We find evidence that a greater tortuosity, greater degree of branching of a pore network and narrower pore throats are detrimental to MT and contribute to the risk of clogging and reduction of connected porosity. We suggest that a tortuosity of less than 2 is critical in promoting greater precipitation per unit volume and should be considered alongside porosity and permeability when assessing reservoirs for geological carbon storage (GCS). We also show that the dominant influence on precipitated mass is the Damköhler number, or reaction rate, rather than the availability of reactive minerals, suggesting that this should be the focus when engineering effective subsurface carbon storage reservoirs for long term security.