The efficiency of mixing in density-driven natural-convection is largely governed by the aquifer permeability, which is heterogeneous in practice. The character (fingering, stable mixing or channeling) of flow-driven mixing processes depends primarily on the permeability heterogeneity character of the aquifer, i.e., on its degree of permeability variance (Dykstra-Parsons coefficient) and the correlation length. Here we follow the ideas of Waggoner et al. (1992) to identify different flow regimes of a density-driven natural convection flow by numerical simulation. Heterogeneous fields are generated with the spectral method of Shinozuka and Jan (1972), because the method allows the use of power-law variograms. We observe from our simulations that the rate of mass transfer of CO2 into water is higher for heterogeneous media.