Hostname: page-component-7c8c6479df-nwzlb Total loading time: 0 Render date: 2024-03-28T05:40:19.872Z Has data issue: false hasContentIssue false

Onset of convection in a gravitationally unstable diffusive boundary layer in porous media

Published online by Cambridge University Press:  01 February 2006

A. RIAZ
Affiliation:
Department of Petroleum Engineering, Stanford University, Stanford, CA 94305, USA
M. HESSE
Affiliation:
Department of Petroleum Engineering, Stanford University, Stanford, CA 94305, USA
H. A. TCHELEPI
Affiliation:
Department of Petroleum Engineering, Stanford University, Stanford, CA 94305, USA
F. M. ORR
Affiliation:
Department of Petroleum Engineering, Stanford University, Stanford, CA 94305, USA

Abstract

We present a linear stability analysis of density-driven miscible flow in porous media in the context of carbon dioxide sequestration in saline aquifers. Carbon dioxide dissolution into the underlying brine leads to a local density increase that results in a gravitational instability. The physical phenomenon is analogous to the thermal convective instability in a semi-infinite domain, owing to a step change in temperature at the boundary. The critical time for the onset of convection in such problems has not been determined accurately by previous studies. We present a solution, based on the dominant mode of the self-similar diffusion operator, which can accurately predict the critical time and the associated unstable wavenumber. This approach is used to explain the instability mechanisms of the critical time and the long-wave cutoff in a semi-infinite domain. The dominant mode solution, however, is valid only for a small parameter range. We extend the analysis by employing the quasi-steady-state approximation (QSSA) which provides accurate solutions in the self-similar coordinate system. For large times, both the maximum growth rate and the most dangerous mode decay as $t^{1/4}$. The long-wave and the short-wave cutoff modes scale as $t^{1/5}$ and $t^{4/5}$, respectively. The instability problem is also analysed in the nonlinear regime by high-accuracy direct numerical simulations. The nonlinear simulations at short times show good agreement with the linear stability predictions. At later times, macroscopic fingers display intense nonlinear interactions that significantly influence both the front propagation speed and the overall mixing rate. A dimensional analysis for typical aquifers shows that for a permeability variation of 1—3000 mD, the critical time can vary from 2000 yrs to about 10 days while the critical wavelength can be between 200 m and 0.3 m.

Type
Papers
Copyright
© 2006 Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)