# InterPore2018 New Orleans

14-17 May 2018
New Orleans
US/Central timezone

## Modeling the dissolution-driven convection as a Rayleigh-Benard problem

14 May 2018, 10:23
15m
New Orleans

#### New Orleans

Oral 20 Minutes MS 3.11: fundamental aspects of geological storage of CO2

### Speaker

Prof. Layachi Hadji (The University of Alabama)

### Description

We examine the linear and weakly nonlinear stability analyses of the dissolution-driven convection induced by the sequestration of carbon dioxide in a geological formation. The mathematical model consists of Darcy's equation, the conservation of mass and the conservation of solute equations. The model accounts for anisotropy in both carbon diffusion and permeability which is modeled by a decaying exponential function of depth. The presence of a first order reaction between the carbon-rich brine and host mineralogy is also included. We prescribe either Neumann or Dirichlet boundary condition for the concentration of carbon dioxide at the rigid upper and lower walls that bound a layer of infinite horizontal extent. We consider a Rayleigh-Taylor-like base state consisting of a carbon-rich heavy layer overlying a carbon-free lighter layer and seek the critical thickness at which this configuration becomes unstable. With this approach, standard mathematical methods that were successfully used in the study of Rayleigh-Benard convection can be applied to this problem. We quantify the influence of carbon diffusion anisotropy, permeability dependence on depth and the presence of the chemical reaction on the threshold instability conditions and associated flow patterns using the classical normal modes approach. The critical Rayleigh number and corresponding wavenumber are found to be independent of the depth of the formation. The weakly nonlinear analysis is performed using long wavelength asymptotics, the validity of which is limited to small Damk\"{o}hler numbers. We derive analytical expressions for the solute flux at the interface, the location of which corresponds to the minimum depth of the boundary layer at which instability sets in. We show that the interface acts as a sink leading to the formation of a self-organized exchange between descending carbon-rich brine and ascending carbon free brine. Plots of the high order perturbation terms for the concentration successfully reproduce the fingering pattern that is typically observed in experiments and full numerical simulations. Using the derived interface flux conditions, we put forth differential equations for the time evolution of the upward migration of the interface as the dissolution process progresses. We solve for the terminal time when the interface reaches the top boundary thereby quantifying the time it takes for an initial amount of injected super-critical Carbon dioxide to be completely dissolved. We also consider the case where the interface migration is accompanied by interface deformations that conform to the convection pattern.

### References

J. Ennis-King and L. Paterson, "Role of convective mixing in the long-term storage of carbon dioxide in deep saline
formations," SPE J. {\bf 10}, 349-356 (2005).

A. Riaz, M. Hesse, A. Tchelepi, and M. Orr, "Onset of convection in a gravitationally unstable diffusive boundary layer in porous media," J. Fluid Mech. {\bf 548}, 87-111 (2005).

J. Ennis-King,I. Preston, and L. Paterson, "Onset of convection in anisotropic porous media subject to rapid change in boundary conditions," Phys. Fluid, {\bf 17}, 84-107 (2005).

X. Xu, S. Chen S, and D. Zhang, "Convective stability analysis of the long term storage of carbon dioxide in deep saline aquifer," Advances in Water Resources {\bf 29}, 397-407 (2006).

C. Yang and Y. Gu, "Accelerated mass transfer of CO$_2$ in reservoir brine due to density-driven natural convection at high pressures and elevated temperatures," Ind. Eng. Chem. {\bf 45}, 2430-2436 (2006)

J.A. Neufeld, D. Vella, and H. E. Huppert, "The effect of a fissure on storage in a porous medium," J. Fluid Mech., {\bf 639} 239-259 (2007).

T.J. Kneafsey and K. Pruess, "Laboratory flow experiment for visualizing carbon dioxide-induced density driven brine convection," Transp Porous Med {\bf 82}, 123-139 (2010).

C. Slim and S. Ramakrishnan, "Onset and cessation of time-dependent, dissolution-driven convection in porous media," Phys. Fluids {\bf 22}, 124103 (2010).

J.A. Neufeld, A. Hesse, A. Riaz, A. Hallworth, A. Tchelepi, and H.E. Huppert, "Convective dissolution of carbon dioxide in saline aquifers," Geophys. Res. Lett. {\bf 37}, L22404 (2010).

D.R. Hewitt, J.A. Neufeld and J.R. Lister, "Convective shutdown in a porous medium at high Rayleigh number," J. Fluid Mech. {\bf 719}, 551â€“586 (2013).

A.C. Slim, M.M. Bandi, J.C. Miller, and L. Mahadevan, "Dissolution-driven convection in a Hele-Shaw cell", Phys. Fluids
{\bf 25}, 024101 (2013).

H.E. Huppert and J.A. Neufeld, " The fluid mechanics of carbon dioxide sequestration," J. Fluid. Mech. {\bf 46}, 255-272 (2014).

A.A. Hill and M.R. Morad, "Convective stability of carbon sequestration in anisotopic porous media," Proc. R. Soc. A {\bf 470}, 20140373, (2014)

A.C. Slim, "Solutal-Convection regimes in a two-dimeensional porous medium," J. Fluid Mech. {\bf 741} 461-491 (2014).

M. De Paoli, F., Zonta, and A. Soldati, "Influence of anisotropic permeability on convection in porous media: Implication for geological CO$_2$ sequestration," Phys. Fluids {\bf 28} 056601 (2016).

C.T. Wanstall and L. Hadji, "A step function density model for the convective stabilty of CO$_2$ geological sequestration," J. Eng. Math. (DOI 10.1007/s10665-017-9907-9) (2017).

L. Vo and L. Hadji, " Weakly nonlinear convection induced by the sequestration of $CO_2$ in a perfectly impervious geological formation, Physics of Fluids, (To appear, Dec. 2017)

K. Klinkenberg, "The permeability of porous media to liquids and gases," Drilling and production practice. American Petroleum Institute, (1941).

E.R. Lapwood, "Convection of a fluid in porous media," Proc. Phil. Soc. {\bf 44} 508-521 (1948).