In the past 10 years the mature hydrocarbon province the West Netherlands Basin has hosted rapidly expanding geothermal development. Upper Jurassic to Lower Cretaceous strata from which gas and oil had been produced since the 1950s became targets for geothermal exploitation. The extensive publicly available subsurface data including seismic surveys, several cores and logs from hundreds of hydrocarbon wells, combined with understanding of the geology after decades of hydrocarbon exploitation, facilitated the offtake of geothermal exploitation. Whilst the first geothermal projects proved the suitability of the permeable Upper Jurassic to Lower Cretaceous sandstones for geothermal heat production, they also made clear that much detail of the aquifer geology is not yet fully understood. The aquifer architecture varies significantly across the basin because of the syn-tectonic sedimentation. The graben fault blocks that contain the geothermal targets experienced a different tectonic history compared to the horst and pop-up structures that host the hydrocarbon fields from which most subsurface data are derived. Accurate prediction of the continuity and thickness of aquifers is a prerequisite for efficient geothermal well deployment that aims at increasing heat recovery while avoiding the risk of early cold-water breakthrough. The potential recoverable heat and the current challenges to enhance further expansion of heat exploitation from this basin are evident. This paper presents an overview of the current understanding and uncertainties of the aquifer geology of the Upper Jurassic to Lower Cretaceous strata and discusses new sequence-stratigraphic updates of the regional sedimentary aquifer architecture.

The primary challenge for efficient geothermal doublet design and deployment is the adequate prediction of the size, shape, lateral extent and thickness (or aquifer architecture) of aquifers. In the West Netherlands Basin, fluvial Lower Cretaceous sandstone-rich successions form the main aquifers for geothermal heat exploitation. Large variations in the thickness of these successions are recognised in currently active doublet systems that cannot be explained. This creates an uncertainty in aquifer thickness prediction, which increases the uncertainty in doublet lifetime prediction as it has an impact on net aquifer volume. The goal of this study was to improve our understanding of the thickness variations and regional aquifer architecture of the Nieuwerkerk Formation geothermal aquifers. For this purpose, new palynological data were evaluated to correlate aquifers in currently active doublet systems based on their chronostratigraphic position and regional Maximum Flooding Surfaces. Based on the palynological cuttings analysis, the fluvial interval of the Nieuwerkerk Formation was subdivided into two successions: a Late Ryazanian to Early Valanginian succession and a Valanginian succession. Within these successions trends were identified in sandstone content. In combination with seismic interpretation, maps were constructed that predict aquifer thickness and their lateral extent in the basin. The study emphasises the value of palynological analyses to reduce the uncertainty of fluvial hot sedimentary aquifer exploitation.