Numerical simulations on the influence of matrix diffusion to carbon sequestration in double porosity fissured aquifers

The double porosity model for fissured rocks, such as limestones and dolomites, has some features that may be relevant for carbon sequestration. Numerical simulations were conducted to study the influence of matrix diffusion on the trapping mechanisms relevant for the long-term fate of CO2 injected in fissured rocks. The simulations show that, due to molecular diffusion of CO2 into the rock matrix, dissolution trapping and hydrodynamic trapping are more effective in double porosity aquifers than in an equivalent porous media. Mineral trapping, although assessed indirectly, is also probably more relevant in double porosity aquifers due to the larger contact surface and longer contact time between dissolved CO2 and rock minerals. However, stratigraphic/structural trapping is less efficient in double porosity media, because at short times CO2 is stored only in the fissures, requiring large aquifer volumes and increasing the risk associated to the occurrence of imperfections in the cap-rock through which leakage can occur. This increased risk is also a reality when considering storage in aquifers with a regional flow gradient, since the CO2 free-phase will move faster due to the higher flow velocities in fissured media and discharge zones may be reached sooner.