The main objective of this study is to quantify interference effects between production and injection wells for geothermal projects that operate within the same reservoir. The secondary objective is to study the time it takes for the reservoir unit to thermally recharge after production has stopped. Two reservoir simulators are benchmarked to decide which one to use for reaching these objectives. The Ammerlaan, Duijvestijn and DAP doublets, all targeting the the Delft Sandstone Member in the West Netherlands Basin, serve as a case study. A literature review is presented to gain a better understanding of the geological history and the structural setting. A box-model of the study area is created and used to benchmark the two reservoir simulators. A static reservoir model is developed in Petrel through seismic interpretation, structural modelling, facies modelling and petrophysical modelling. Populating this static model with dynamic properties allows us to perform reservoir simulations in Eclipse100 and study the propagation of pressure and temperature. A discrete parameter analysis aims to capture the uncertainties that are associated with reservoir modelling and simulation. In this study, data is used from seismic surveys, wireline logs, core studies, a cuttings study and the monthly production volumes of the Ammerlaan and Duijvestijn doublets. It is shown that interference on temperature has a long term (20-30 years) effect on the Duijvestijn (positively) and the Ammerlaan and DAP doublets (negatively). The combined total energy production of the three doublets over 100 years is small with a decrease of 3% due to temperature interference, compared to when running the doublets in stand-alone configuration. Interference on pressure has a short term (days) effect on the achieved flow rate when the injection well cannot reach its target injection rate and is constraint on the maximum allowable pressure at which it is allowed to inject. This occurs under the low permeability scenarios and can be observed when the pressure wave arrives at the neighbouring injector. The Duijvestijn, Ammerlaan and DAP doublets all benefit from this through a combined increase in energy production of 8% over 100 years of production, averaged over all scenarios. During thermal recharge after production, the average reservoir temperature increases asymptotically. Temperature recharges to 96.1% to 97.4% of the initial reservoir temperature after 1000 years of recharging under different thickness and conductivity scenarios. The absence of a thermal influx at the bottom of the reservoir model is limiting the capacity of the reservoir to recharge to its initial average temperature. Interviews are conducted to investigate the implications that the findings of this research have on the policy measures for geothermal projects. Suggestions for changes in policy are made to optimize the recovery of heat in the subsurface.