Geothermal energy extraction changes reservoir pressure and temperature due to colder fluid injection, which impact the stress field. In naturally fractured or faulted reservoirs, stress changes introduce risks of fault reactivation and induced seismicity. In this study, the capabilities of the Delft Advanced Research Terra Simulator (DARTS) framework have been extended to account for the thermo-poroelastic response relevant to geothermal applications. In the proposed modeling approach, a low-fidelity model is used for geomechanics where stresses are calculated from pressure and temperature changes. The fault slippage is calculated using Mohr-Coulomb friction criteria. The developed model has been applied to heterogeneous geothermal field to evaluate how energy production and thermal front can affect the geomechanical state of the field and trigger induced seismicity.