Numerical simulation of stable isotope fractionation with density changes in regional-scale confined geothermal aquifers