A step function density profile model for the convective stability of CO $$_2$$ 2 geological sequestration

Abstract

The convective stability associated with carbon sequestration is usually investigated by adopting an unsteady diffusive basic profile to account for the space and time development of the carbon-saturated boundary layer instability. The method of normal modes is not applicable due to the time dependence of the nonlinear base profile. Therefore, the instability is quantified either in terms of critical times at which the boundary layer instability sets in or in terms of long-time evolution of initial disturbances. This paper adopts an unstably stratified basic profile having a step function density with top heavy carbon-saturated layer (boundary layer) overlying a lighter carbon-free layer (ambient brine). The resulting configuration resembles that of the Rayleigh–Taylor problem with buoyancy diffusion at the interface separating the two layers. The discontinuous reference state satisfies the governing system of equations and boundary conditions and pertains to an unstably stratified motionless state. Our model accounts for anisotropy in both diffusion and permeability and chemical reaction between the carbon dioxide-rich brine and host mineralogy. We consider two cases for the boundary conditions, namely an impervious lower boundary with either a permeable (one-sided model) or poorly permeable upper boundary. These two cases possess neither steady nor unsteady unstably stratified equilibrium states. We proceed by supposing that the carbon dioxide that has accumulated below the top cap rock forms a layer of carbon-saturated brine of some thickness that overlies a carbon-free brine layer. The resulting stratification remains stable until the thickness, and by the same token, the density, of the carbon-saturated layer is sufficient to induce the fluid to overturn. The existence of a finite threshold value for the thickness is due to the stabilizing influence of buoyancy diffusion at the interface between the two layers. With this formulation for the reference state, the stability calculations will be in terms of critical boundary layer thickness instead of critical times, although the two formulations are homologous. This approach is tractable by the classical normal-mode analysis. Even though it yields only conservative threshold instability conditions, it offers the advantage for an analytically tractable study that puts forth expressions for the carbon concentration convective flux at the interface and explores the flow patterns through both linear and weakly nonlinear analyses.

http://ift.tt/2rsM6K9