The sand-prone fluvial Delft Sandstone Member of the Nieuwerkerk Formation is the main target for geothermal exploitation in the West Netherlands Basin. Uncertainty reduction in the modelling of permeability heterogeneity and connectivity of the sandstone reservoir units is of prime importance in the low-profit geothermal energy market. Lithostratigraphic correlation of the Delft Sandstone Member dates back to the 1970s, and considers the fluvial deposits as Valanginian (140-134 Ma) in age and synchronously deposited throughout the basin. Analysis of newly-acquired palynological data, in combination with well-log correlation and seismic interpretation, enabled the construction of a sequence-stratigraphic framework for the Nieuwerkerk Formation, in which the fluvial sandstones were deposited in a much wider age range than previously believed, from Late Ryazanian to Valanginian (142-134 Ma). The depocentre of fluvial deposition shifted in time from NE to SW across the width of the basin. This diachronous development resulted in a discontinuous spatial arrangement of fluvial sandstones. The new reservoir architecture model will aid in the well placement of geothermal doublets and the assessment of interference risk of adjacent geothermal projects.

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.