Diffusive transport has implications for the long-term status of underground storage of hydrogen (H_{2}) fuel and carbon dioxide (CO_{2}), technologies which are being pursued to mitigate climate change and advance the energy transition. Once injected underground, CO_{2} and H_{2} will exist in multiphase fluid-water-rock systems. The partially soluble injected fluids can flow through the porous rock in a connected plume, become disconnected and trapped as ganglia surrounded by groundwater within the storage rock pore space, and also dissolve and migrate through the aqueous phase once dissolved. Recent analyses have focused on the concentration gradients induced by differing capillary pressure between fluid ganglia which can drive diffusive transport ("Ostwald ripening"). However, studies have neglected or excessively simplified important factors, namely the nonideality of gases under geologic conditions, the opposing equilibrium state of dissolved CO_{2} and H_{2} driven by the partial molar density of dissolved solutes, and entropic and thermodiffusive effects resulting from geothermal gradients. We conduct an analysis from thermodynamic first principles and use this to provide numerical estimates for CO_{2} and H_{2} at conditions relevant to underground storage reservoirs. We show that while diffusive transport in isothermal systems is upwards for both gases, as indicated by previous analysis, entropic contributions to the free energy are so significant as to cause a reversal in the direction of diffusive transport in systems with geothermal gradients. For CO_{2}, even geothermal gradients less than 10^{∘}C/km (far less than typical gradients of 25^{∘}C/km) are sufficient to induce downwards diffusion at depths relevant to storage. Diffusive transport of H_{2} is less affected but still reverses direction under typical gradients, e.g., 30^{∘}C/km, at a depth of 1000m. This reversal occurs independent of the solute's thermophobicity or thermophilicity in aqueous solutions. The entropic contribution also modifies the magnitude of flux where geothermal gradients are present, with the largest diffusive fluxes estimated for CO_{2} with a 30^{∘}C/km gradient, despite the higher diffusion coefficient of H_{2}. We find a maximum flux on the order of 10^{-13} mol/(cm^{2}s) for CO_{2} in the 30^{∘}C/km scenario; significantly lower than literature estimates for maximum convective fluxes in moderate to high permeability formations. Contrary to previous studies, we find that in diffusion and convection will likely work in concert-both driving CO_{2} downwards, and both driving H_{2} upwards-for conditions representative of their respective storage reservoirs.
Read full abstract