Upscaling of CO2 vertical migration through a periodic layered porous medium: The capillary-free and capillary-dominant cases

• Published in: Advances in water resources 2010-09, Vol.33 (9), p.1164-1175
• Main Authors: MOUCHE, Emmanuel, HAYEK, Mohamed, MÜGLER, Claude
• Format: Article
• Language: English
• Subjects: Buoyancy; Carbon dioxide; Continental interfaces, environment; Continuity; Earth sciences; Earth, ocean, space; Exact sciences and technology; Flux; Fluxes; Hydrology. Hydrogeology; Mathematical analysis; Mathematical models; Migration; Ocean, Atmosphere; Saturation; Sciences of the Universe
• ISSN: 0309-1708; 1872-9657
• DOI: 10.1016/j.advwatres.2010.07.005

We present an upscaled model for the vertical migration of a CO2 plume through a vertical column filled with a periodic layered porous medium. This model may describe the vertical migration of a CO2 plume in a perfectly layered horizontal aquifer. Capillarity and buoyancy are taken into account and semi-explicit upscaled flux functions are proposed in the two following cases: (i) capillarity is the main driving force and (ii) buoyancy is the only driving force. In both cases, we show that the upscaled buoyant flux is a bell-shaped function of the saturation, as in the case of a homogeneous porous medium. In the capillary-dominant case, we show that the upscaled buoyant flux is the harmonic mean of the buoyant fluxes in each layer. The upscaled saturation is governed by the continuity of the capillary pressure at the interface between layers. In the capillary-free case, the upscaled buoyant flux and upscaled saturation are determined by the flux continuity condition at the interface. As the flux is not continuous over the entire range of saturation, the upscaled saturation is only defined where continuity is verified, i.e. in two saturation domains. As a consequence, the upscaled buoyant flux is described by a piecewise continuous function. Two analytical approximations of this flux are proposed and this capillary-free upscaled model is validated for two cases of heterogeneity. Upscaled and cell averaged saturations are in good agreement. Furthermore, the proposed analytical upscaled fluxes provide satisfactory approximations as long as the saturation set at the inlet of the column is in a range where analytical and numerical upscaled fluxes are close.