Geochemistry of silicate-rich rocks can curtail spreading of carbon dioxide in subsurface aquifers.

Pools of carbon dioxide are found in natural geological accumulations and in engineered storage in saline aquifers. It has been thought that once this CO2 dissolves in the formation water, making it denser, convection streams will transport it efficiently to depth, but this may not be so. Here, we assess theoretically and experimentally the impact of natural chemical reactions between the dissolved CO2 and the rock formation on the convection streams in the subsurface. We show that, while in carbonate rocks the streaming of dissolved carbon dioxide persists, the chemical interactions in silicate-rich rocks may curb this transport drastically and even inhibit it altogether. These results challenge our view of carbon sequestration and dissolution rates in the subsurface, suggesting that pooled carbon dioxide may remain in the shallower regions of the formation for hundreds to thousands of years. The deeper regions of the reservoir can remain virtually carbon free.

The behaviour of Sal'nikov's reaction, P --> A --> B, in a spherical batch reactor with the diffusion of heat and matter.

The behaviour of a simple chemical reaction, occurring with the release of heat in a closed batch reactor, is considered for the situation when matter and heat are transported only by diffusive processes; thus, the reacting fluid has negligible velocity, so that heat transfer is by thermal conduction. The reaction is Sal'nikov's, which consists of two, consecutive first-order steps, producing a product B, from a precursor P, via an active intermediate A, in P --> A --> B. The first of these steps is assumed to be thermoneutral, with zero activation energy, whilst the second is exothermic, with an appreciable activation energy. These features make Sal'nikov's reaction the simplest to display thermokinetic oscillations that characterise many, more complex schemes, e.g. cool flames in hydrocarbon combustion. This study involves identifying the regions of parameter space, in which these oscillations in the temperature and the concentration of the intermediate A occur, by means of numerical simulation. These regions are compared with previous analytical stability analyses in one-dimensional systems. It was found that oscillations occur over a much larger range of conditions in the case considered here, i.e. a reactor with spherical symmetry, than in the simple 1-D case, previously studied by Gray and Scott (P. Gray and S. K. Scott, Chemical Oscillations and Instabilities, Clarendon Press, Oxford, 1990, pp. 264-291). In addition, approximate analytical solutions for the temperature and concentration of A are presented for two limiting cases of non-oscillatory behaviour. These analytical solutions have been verified by comparison with full numerical solutions of the governing equations.