Numerical simulation of fracture propagation and production performance in a fractured geothermal reservoir using a 2D FEM-based THMD coupling model