Deep geothermal energy in northern England: Insights from 3D finite difference temperature modelling

Abstract

Many of the most widely used deep geothermal resource maps for the UK are produced by contouring around sparsely distributed and often unreliable data points. We thus present a MATLAB-based 3D finite difference temperature modelling methodology, which provides a means for producing more resolute and geologically realistic versions of these maps. Our case study area in northern England represents an area where both sedimentary basins and radiothermal granite bodies comprise potential geothermal resources. We divide our 3D model into geological units, which are then assigned separate thermal properties. Assuming conductive heat transfer and steady-state and fixed boundary conditions, we calculate 3D regional subsurface temperature. Due to our averaging technique for thermal properties, the resolution of our geological model is scarcely compromised with respect to similar finite element methods. One predicted ‘hot spot’ at 1 km depth in the central part of our case study area corresponds with the granitic North Pennine Batholith. Other shallow hot spots correspond with thermally insulating sedimentary rock units and geological structures that incorporate these units. Predictive heat flow density maps highlight areas with accelerated surface heat flow associated with shallow conductive basement rock and heat producing granite bodies. Our predicted subsurface temperatures show broad similarities with measured equilibrium borehole temperatures. Inaccuracies may relate to convective heat transfer involving fault systems, or input variables relating to the geological model. Our predictive subsurface temperature and heat flow density maps are more resolute and geologically realistic relative to pre-existing contoured maps. The method presented here represents a useful tool for understanding controls on subsurface temperature distribution and geothermal potential.