Coupled CO₂ injection well flow model to assess thermal stresses under geomechanical uncertainty

We develop a two-phase transient non-isothermal wellbore flow model, augmented with a radial heat conduction in the annulus, casing, and the reservoir. Using the available data for a saline aquifer in Denmark, we build a one-dimensional geomechanical well model and assess the stresses at the wellbore wall using the analytical Kirsch formula. Using the temperature at the wellbore wall, we calculate the corresponding thermal stresses. Furthermore, we assess the impact of the uncertainty in thermal expansion coefficients on the magnitudes of thermal stresses. For the cases considered, the magnitude of the changes in the critical pressure and in the fracture pressure with and without thermal stresses does not exceed 3%.