In this work, the propagation of a thermal-hydraulic front caused by cold fluid injection is investigated, using a finite element (FE) model. The objective is to simulate the propagation of temperature changes in the subsurface and to quantify the resulting thermally induced stresses. The model is based on the setup of a geothermal powerplant in the greater Munich area. In the vicinity of the injection well of this powerplant seismic events occurred five years after commissioning. The powerplant uses the Malm aquifer as geothermal reservoir and reinjects cold fluid in the vicinity of a fault zone. In order to approach to the complex reservoirs’ hydrogeology and the fault zone characteristics, a mesh generator was created using MATLAB. The generator creates a model based on regular cube-shaped elements, which can be separately parameterized and thus grouped to geological units. Furthermore, the cube's surfaces can be integrated as discrete, two-dimensional elements, which can also be separately parameterized. Thus, a 3D model containing a 3D fault zone was created with a reservoir containing a regularly distributed 2D fracture network and wells as 1D line element. ... mehrA scenario analysis was used to investigate the influence of the fault and the fractures on the thermal-hydraulic front propagation. The results of the simulations show that especially the fracture network can significantly influence the range of the area influenced thermally by fluid injection. The fluid flow within the fault zone can lead to temperature changes in several hundred meters underneath the transition zone between the reservoir and the crystalline basement. To quantify the resulting thermo-hydraulic stresses analytical approaches were used. The results show that the thermal-induced stresses exceed the hydraulic induced stresses clearly and can have influences on a faults' integrity.