A numerical study of the thermohydraulics of an enhanced geothermal system project in Turkey is presented. The solid structures are modelled as porous media, using the numerically determined hydraulic fracturing data of other authors. The influence of several numerical modelling aspects such as the domain size, grid resolution, temporal resolution as well as the discretization scheme are investigated and assessed to obtain highly accurate numerical solutions under the applied modelling assumptions. Using the suggested mathematical and numerical model, different production scenarios are investigated.