Groundwater heat pump systems cause thermal anomalies in the aquifer that can impact upon downstream systems and reduce their efficiency. Therefore, it is important to optimally position the extraction and injection wells of such systems to avoid negative interactions and maximize the thermal potential of the aquifer. This paper presents a new method to determine optimal well layouts of groundwater heat pumps using the adjoint approach, which is an efficient way to solve the underlying PDE-constrained optimization problem. An integral part of the method is the numerical groundwater simulation, which here is based on the finite element method. In addition, a multi-start initialization strategy is introduced in an attempt to
better reach the global optimum. The method is applied to a real case study with 10 groundwater heat pumps, i.e. 20 wells, and two optimization scenarios with different natural groundwater temperatures. In both scenarios, the optimization method successfully determines a well layout that maximizes groundwater temperatures at all extraction wells. Comparing the results from these scenarios demonstrates that hydro-geological conditions can have a significant impact on the optimal well layout. The proposed method is equally applicable to systems with multiple extraction and injection wells and can be extended to various other shallow geothermal applications, such as combined heating and cooling systems.