Investigation of Thermal Stress of Cement Sheath for Geothermal Wells during Fracturing

Geothermal energy development has increasingly been studied in recently decades because of its renewable and sustainable features. It can be divided into two categories: traditional geothermal (hydrothermal) systems and enhanced geothermal systems (EGS) based on the type of exploitation. The hot dry rock (HDR) in the EGS incorporates about 80% of all thermal energy, and its value is about 100–1000 times that of fossil energy. It is pivotal for geothermal wells to improve the flow conductivity of the HDR mass, enhance the communication area of natural fractures, and constitute the fracture network between injection and production wells by hydraulic treatments. While the wellbore temperature significantly decreases because of fracturing, fluid injection will induce additional thermal stresses in the cement sheath, which will aggravate its failure. Considering the radial nonuniform temperature change, this paper proposes a new thermal stress model for a casing-cement sheath-formation combined system for geothermal wells during fracturing based on elastic mechanics and thermodynamics theory. This model is solved by the Gaussian main elimination method. Based on the analytical model, the thermal stresses of cement sheath have been analyzed. The effects of the main influencing parameters on thermal stresses have also been investigated. Results show that the radial and axial tensile thermal stresses are both obviously larger than tangential tensile thermal stress. The maximum radial and axial thermal stresses always occur at the casing interface while the location of the maximum tangential thermal stress varies. Generally, thermal stresses are more likely to induce radial and axial micro cracks in the cement sheath, and the cement sheath will fail more easily at the casing interface in fracturing geothermal wells. For integrity protection of the cement sheath, a proper decrease of casing wall thickness, casing linear thermal expansion coefficient, cement sheath elasticity modulus, and an increase of the fracturing fluid temperature has been suggested.