This paper explores the possibility of a development project with a geothermal well doublet in the Pisa plain, Italy. The performance of the system has been evaluated with a 3-dimensional field-scale numerical model that simulates the evolution of temperature and pressure conditions in the aquifer, under different exploitation scenarios. Coupled groundwater flow and thermal transport processes in the reservoir are considered together with non-Darcy fluid flow in the wellbores, and heat exchange between boreholes and surrounding rock formations. Calculations are performed with a parallelized version of the wellbore-reservoir simulator T2Well. This code allows for the efficient modeling of coupled hydraulic-thermal processes over a domain about 40 km(2) wide and 1.5 km thick. Simulation results indicate that the energy of the reservoir is sufficient for the designed extraction rate (between 80 and 150 m(3)/h), but also suitable for much larger rates, up to 250 m3/h. Although aimed at assessing the long-term performance of a specific system, this modeling approach could be profitably applied for the design of similar projects elsewhere. (C) 2017 Elsevier Ltd. All rights reserved.