Transport in Porous Media, Vol.33, No.1, 29-63, 1998
Layered thermohaline convection in hypersaline geothermal systems
Thermohaline convection occurs in hypersaline geothermal systems due to thermal and salinity effects on liquid density. Because of its importance in oceanography, thermohaline convection in viscous liquids has received more attention than thermohaline convection in porous media. The fingered and layered convection patterns observed in viscous liquid thermohaline convection have been hypothesized to occur also in porous media. However, the extension of convective dynamics from viscous liquid systems to porous media systems is complicated by the presence of the solid matrix in porous media. The solid grains cause thermal retardation, hydrodynamic dispersion, and permeability effects. We present simulations of thermohaline convection in model systems based on the Salton Sea Geothermal System, California, that serve to point out the general dynamics of porous media thermohaline convection in the diffusive regime, and the effects of porosity and permeability, in particular. We use the TOUGH2 simulator with residual formulation and fully coupled solution technique for solving the strongly coupled equations governing thermohaline convection in porous media. We incorporate a model for brine density that takes into account the effects of NaCl and CaCl2. Simulations show that in forced convection, the increased pore velocity and thermal retardation in low-porosity regions enhances brine transport relative to heat transport. In thermohaline convection, the heat and brine transport are strongly coupled and enhanced transport of brine over heat cannot occur because buoyancy caused by heat and brine together drive the flow. Random permeability heterogeneity has a limited effect if the scale of flow is much larger than the scale of permeability heterogeneity. For the system studied here, layered thermohaline convection persists for more than one million years for a variety of initial conditions. Our simulations suggest that layered thermohaline convection is possible in hypersaline geothermal systems provided the vertical permeability is smaller than the horizontal permeability, as is likely in sedimentary basins such as the Salton Trough. Layered thermohaline convection can explain many of the observations made at the Salton Sea Geothermal System over the years.