Hostname: page-component-77c89778f8-m8s7h Total loading time: 0 Render date: 2024-07-17T18:13:09.076Z Has data issue: false hasContentIssue false

Regional geothermal aquifer architecture of the fluvial Lower Cretaceous Nieuwerkerk Formation – a palynological analysis

Published online by Cambridge University Press:  15 January 2018

Cees J.L. Willems*
Affiliation:
Department of Geoscience and Engineering, Delft University of Technology, Delft, the Netherlands
Andrea Vondrak
Affiliation:
PanTerra Geoconsultants BV, Leiderdorp, the Netherlands
Dirk K. Munsterman
Affiliation:
TNO – Geological Survey of the Netherlands, P.O. Box 80015, 3508 TA Utrecht, the Netherlands
Marinus E. Donselaar
Affiliation:
Department of Geoscience and Engineering, Delft University of Technology, Delft, the Netherlands
Harmen F. Mijnlieff
Affiliation:
TNO – Geological Survey of the Netherlands, P.O. Box 80015, 3508 TA Utrecht, the Netherlands
*
*Corresponding author. Email: willems.cjl@gmail.com
Rights & Permissions [Opens in a new window]

Abstract

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.

Type
Review
Copyright
Copyright © Netherlands Journal of Geosciences Foundation 2018 

Introduction

In geothermal exploitation of sedimentary rocks, it is crucial to adequately predict the regional aquifer distribution. Often geological data are sparse and property extrapolation is required over large distances. This is especially challenging for fluvial aquifers, which are notorious for lateral variation in lithofacies and aquifer properties. The prediction of the regional sandstone distribution (henceforth fluvial architecture) from well logs in fluvial aquifers is often ambiguous because lithofacies distribution could be affected by both allogenic and autogenic processes (e.g. Hajek et al., Reference Hajek, Heller and Sheets2010; Donselaar et al., Reference Donselaar, Cuevas Gozalo and Moyano2013; Flood & Hampson, Reference Flood and Hampson2015; Van Toorenenburg et al., Reference Van Toorenenburg, Donselaar, Noordijk and Weltje2016). This is reflected by large aquifer thickness variations that are recognised in currently active geothermal doublet wells in the West Netherlands Basin (WNB). The fluvial, sandstone-rich successions that form the aquifer of the geothermal HON-GT doublet range in thickness from 50 to 150 m in approximately 1.5 km spaced wells (Fig. 1). In addition, the depth of this aquifer below the top of the marginally marine Rodenrijs Claystone Member (e.g. Van Adrichem Boogaert & Kouwe, Reference Van Adrichem Boogaert and Kouwe1993) ranges from almost 100 m to more than 200 m in different geothermal wells. Up to now, these variations cannot be explained and create uncertainty in the prediction of lifetime and drilling costs of future doublet systems in the WNB.

Fig. 1. Gamma-ray logs of two geothermal doublets HON-GT and PNA-GT. Fault interpretation is based on Duin et al. (Reference Duin, Doornenbal, Rijkers, Verbeek and Wong2006). Well-log correlation is based on ‘End-of-well reports’ (NLOG, 2017).

The aquifer in the geothermal wells of Figure 1 is interpreted as the Delft Sandstone Member which is part of the Lower Cretaceous Nieuwerkerk Formation (e.g. Van Adrichem Boogaert & Kouwe, Reference Van Adrichem Boogaert and Kouwe1993; Den Hartog Jager, Reference Den Hartog Jager, Rondeel, Batjes and Nieuwenhuijs1996; Herngreen & Wong, Reference Herngreen, Wong, Wong, Batjes and De Jager2007; Donselaar et al., Reference Donselaar, Groenenberg and Gilding2015). This member is characterised as a syn-rift, sandstone-rich interval ranging in age from Valanginian to Early Hauterivian, deposited in a meandering fluvial environment. Regional transgression and subsidence resulted in an increasingly marine character of the overlying sediments ranging from the restricted marine Rodenrijs Claystone Member to the marine Rijnland Group (Fig. 2). The interpretation of the Delft Sandstone Member as a single sandstone-rich interval in the upper section of the Nieuwerkerk Formation (Fig. 3A) is derived from lithostratigraphic regional well-log correlation from numerous hydrocarbon wells in the WNB (e.g. Racero-Baena & Drake, Reference Racero-Baena, Drake, Rondeel, Batjes and Nieuwenhuijs1996; Herngreen & Wong, Reference Herngreen, Wong, Wong, Batjes and De Jager2007). This model is commonly used in geothermal exploitation in the basin for doublet design and deployment. However, recent regional stratigraphic studies based on sequence stratigraphic principles did not acknowledge the Delft Sandstone Member (DeVault & Jeremiah, Reference DeVault and Jeremiah2002; Jeremiah et al., Reference Jeremiah, Duxbury and Rawson2010). DeVault & Jeremiah (Reference DeVault and Jeremiah2002) state that because of the syn-rift origin of the Nieuwerkerk Formation, clusters of amalgamated sandstone-rich zones can exist throughout the Nieuwerkerk Formation that do not necessarily form a single continuous sandstone-rich interval (Fig. 3B). The existence of two geological models that describe sandstone distribution in the Nieuwerkerk Formation creates uncertainty for geothermal exploitation because both models have a different impact on possible interference and aquifer thickness prediction. If the aquifer is formed by a single continuous sandstone-rich interval, pressure communication could affect injectivity and productivity of adjacent doublets (Fig. 3A). In contrast, pressure communication is less straightforward if different sandstone-rich zones occur with limited lateral extent. In the example of Figure 3B, claystone-dominated zones can form flow barriers or baffles between doublets 1 and 2. Furthermore, when the aquifer is not formed by a single sandstone-rich zone, the aquifer thickness depends on the lateral extent of the sandstone-rich zones that the doublets can encounter as is illustrated for doublet 2 in Figure 3B. Furthermore, the model in Figure 3B suggests that multiple aquifer targets can be present at deeper and hotter stratigraphic intervals affecting the geothermal potential in the region.

Fig. 2. Stratigraphic column for the Early Cretaceous section in the WNB indicating tectonic activity during deposition of the Rijnland Group, the Nieuwerkerk Formation and the main geothermal aquifers in the WNB: the Rijswijk Sandstone Member and the Delft Sandstone Member (Van Adrichem Boogaert and Kouwe, Reference Van Adrichem Boogaert and Kouwe1993).

Fig. 3. (A) Cartoons illustrating the difference in sandstone distribution in the Nieuwerkerk Formation in schematic strike sections on graben scale according to (A) the Delft Sandstone model and (B) the multiple sandstone-rich zones models.

The goal of this study is to place the fluvial aquifers in a chronostratigraphic framework. The results should decrease the uncertainty in the prediction of aquifer thickness for new doublet systems in the WNB and contribute to optimised doublet design. To reach this goal, palynological samples from drill cuttings are analysed in three geothermal wells, HON-GT-01, HON-GT-02 and PNA-GT-02, to define the chronostratigraphic position of fluvial intervals and identify regional maximum flooding surfaces (MFSs). The analyses are used to create a framework for a well-to-well correlation from which an explanation of aquifer thickness variations in the different doublets is proposed. This explanation is used to interpret regional aquifer architecture in different fault blocks.

Data and methods

Overview

This study was based on a combination of seismic interpretation, gamma-ray (GR) log correlation, and palynological analysis of cuttings. In the seismic interpretation, faults were identified in our study area in the WNB, which were active during deposition of the Nieuwerkerk Formation. In combination with a regional structural interpretation by Duin et al. (Reference Duin, Doornenbal, Rijkers, Verbeek and Wong2006) the lateral extent of these faults was identified. Secondly, by utilising palynological analyses of cuttings the chronostratigraphic position of each aquifer sandstone interval was identified and MFSs were interpreted. This formed the framework of improved geothermal GR well-log correlation. GR logs of 11 geothermal wells in our study area were used to compare fluvial architecture in different fault blocks. All results were finally combined in maps that predict the lateral extent of the sandstone-rich successions in the basin. An overview of the data used in our study is presented in Figure 4.

Fig. 4. Location of the geothermal wells for the GR well-log correlation, cuttings analysis, the outline of the seismic cross-section and the regional structural interpretation by Duin et al. (Reference Duin, Doornenbal, Rijkers, Verbeek and Wong2006).

Structural setting of the Nieuwerkerk Formation

On a seismic section perpendicular to the major fault trend, two seismic horizons were interpreted: the top and base of the Nieuwerkerk Formation. The basin-wide section was derived by merging ten 3D seismic sections (Fig. 4; Vondrak, Reference Vondrak2016). Using horizon flattening of the top of the formation, fault blocks were identified that experienced different tectonic movement affecting fluvial architecture of the Nieuwerkerk Formation. This is derived from thickness differences of the formation between the major faults. Using structural interpretation by Duin et al. (Reference Duin, Doornenbal, Rijkers, Verbeek and Wong2006), the regional outlines of the fault blocks that experienced different tectonic movement during deposition of the Nieuwerkerk Formation were mapped. This result was used as the basis for regional well-log correlations and generation of maps that describe the distribution of sandstone-rich successions.

Palynological analysis

A total of 42 cuttings samples from well PNA-GT-02, 40 samples from HON-GT-01 and 28 samples from HON-GT-02 were analysed. Two additional samples from well VDB-GT-04 (at depths 1890 m and 1910 m) complemented the palynological analysis of Munsterman (Reference Munsterman2012). The samples were processed at the TNO laboratory using the standard sample-processing procedures (e.g. Janssen & Dammers, Reference Janssen and Dammers2008), which involved HCl and HF treatment, and sieving over an 18µm mesh sieve. The well selection was based on the well location in different graben blocks and the total thickness of the Nieuwerkerk Formation that these wells encountered. Larger total thickness could potentially reveal more information from the fluvial interval. The cuttings descriptions and the GR logs in the ‘End-of-well-reports’ (NLOG, 2017), in combination with results from Munsterman (Reference Munsterman2012), provided a basis for the selection of sample depths. The location of the wells in different fault blocks allowed relating differences in fluvial architecture to the syn-tectonic origin of the interval. The palynological analysis consisted of age dating and identification of the Elegans MFS and the Paratollia MFS (e.g. Jeremiah et al., Reference Jeremiah, Duxbury and Rawson2010), which formed the framework of our regional correlation scheme.

Age dating

The age interpretation was based on the last occurrence datum (LOD) of palynomorphs, in particular dinoflagellate cysts (dinocysts), and pollen and spores (sporomorphs). Key-references concerning the palynostratigraphy of the Early Cretaceous from the North Sea region were Davey (Reference Davey1979, Reference Davey1982), Heilmann-Clausen (Reference Heilmann-Clausen1987), Costa & Davey (Reference Costa, Davey and Powell1992), Riding & Thomas (Reference Riding, Thomas and Powell1992), Partington et al. (Reference Partington, Copestake, Mitchener, Underhill and Parker1993), Abbink (Reference Abbink1998), Duxbury et al. (Reference Duxbury, Kadolsky, Johansen, Jones and Simmons1999), Herngreen et al. (Reference Herngreen, Kerstholt and Munsterman2000) and Duxbury (Reference Duxbury2001). The international geological timescale of Gradstein et al. (Reference Gradstein, Ogg, Smith and Ogg2012) was followed.

Sporomorph eco-grouping (SEG) method

The SEG method (Abbink, Reference Abbink1998; Abbink et al., Reference Abbink, Targarona, Brinkhuis and Visscher2001, Reference Abbink, Van Konijnenburg-Van Cittert, Van der Zwan and Visscher2004a,Reference Abbink, Van Konijnenburg-Van Cittert and Visscherb) was used to identify the Paratollia MFS in the fluvial aquifer interval of HON-GT-01, HON-GT-02 and PNA-GT-02. With this method, sporomorph types were related to vegetation eco-groups. Abbink et al. (Reference Abbink, Van Konijnenburg-Van Cittert, Van der Zwan and Visscher2004a, b) classified Jurassic to Lower Cretaceous sporomorphs into six eco-groups. (1) Upland vegetation grows on higher terrain well above ground water level, which is never submerged by water. (2) Lowland vegetation is found on plains with or without freshwater swamps. It is not influenced by salt water. When periodically submerged it is referred to as ‘Wet-Lowland’ otherwise ‘Dry-Lowland’. (3) River vegetation is found on riverbanks and could be periodically submerged. (4) Pioneering vegetation occupies recently developed eco-space that has been previously submerged by seawater. (5) Coastal vegetation is found along the coast. (6) Tidally influenced vegetation is daily influenced by tidal changes and regularly submerged in a salt-water regime. Quantitative analysis of sporomorphs indicated percentages of eco-groups that were represented in the cuttings samples. In the SEG method it is assumed that the lower coastal plain area is reduced during a transgression (Fig. 5A, B). Therefore the relative share of Lowland eco-group vegetation is minimal on the moment of maximum transgression, when a MFS is formed. Based on this assumption, trends in relative representation of eco-groups were related to sea-level fluctuation. MFSs were assigned to samples where the relative share of ‘Upland’ sporomorphs peaked with respect to the ‘Lowland’ eco-group while the marine-influenced eco-groups were poorly represented or absent. Cuttings samples with 10 m intervals were analysed in the 2560–2810 m (measured depth (MD)) interval in HON-GT-01, 2590–2860 m (MD) in HON-GT-02 and 2440–2850 m (MD) in PNA-GT-02. These intervals were selected based on their fluvial origin which was derived from the cuttings description in the ‘End-of-well-reports’ (NLOG, 2017).

Fig. 5. Schematic representation of the impact of (A) low sea level and (B) high sea level on the relative occurrence of eco-groups: 1=Upland, 2=Lowland, 3=River, 4=Pioneer, 5=Coastal, 6=Tidally influenced. Note the lower relative occurrence of the ‘Lowland’ eco group with high sea level in (B). Modified from Abbink et al. (Reference Abbink, Van Konijnenburg-Van Cittert, Van der Zwan and Visscher2004a,Reference Abbink, Van Konijnenburg-Van Cittert and Visscherb).

Results

Seismic interpretation

On the seismic cross-section three half-grabens and one horst block were recognised (Fig. 6A). The fault blocks were referred to as ‘Westland graben’, ‘Pijnacker graben’, ‘VDB graben’ and ‘Bergschenhoek horst’. The interpretation of the top and base of the Nieuwerkerk Formation indicated a lateral thickness variation of the Nieuwerkerk Formation in these grabens and horst, created by syn-depositional fault movement. Horizon flattening of the Top Nieuwerkerk Formation horizon was applied to highlight fault blocks where sedimentation might be affected by this tectonic movement (Fig. 6B). The associated faults that were active during deposition of the formation are highlighted in red. The regional extent of these faults was derived from the structural interpretation by Duin et al. (Reference Duin, Doornenbal, Rijkers, Verbeek and Wong2006). The three grabens are highlighted on the map in Figure 6. These results were used for the comparison of fluvial reservoir architecture in these three fault blocks.

Fig. 6. (A) Seismic section with interpretation of faults as well as top (yellow) and base (pink) of the Nieuwerkerk Formation horizons. (B) Horizon flattening on top Nieuwerkerk Formation. The outline of the seismic section and the outline of the interpreted faults are indicated on the map.

Palynological analysis

Palynological age dating formed the basis of the GR well-log correlation scheme. An overview of the results is presented in Figure 7 and Table 1. A detailed description of the analysis of all samples is presented in Appendix 1 (as Supplementary Material at https://doi.org/10.1017/njg.2017.23). Two MFSs were identified in HON-GT-01 and HON-GT-02 and four MFSs in PNA-GT-02. The MFS that is close to the Early Valanginian to Late Ryazanian boundary is associated with the Paratollia MFS in all wells. In well PNA-GT-02, this was based on the LOD of Stiphrosphaeridium dictyophorum (Sdi) at 2600 m depth and the LOD of Canningia compta (Cco) at 2620 m depth (Appendix 1). In HON-GT-01, the Paratollia MFS was interpreted at 2730 m depth based on the LODs of Canningia compta (and a morphologically closely related Escharisphaeridia spp. at 2730 m) and Perisseiasphaeridium insolitum at depth 2740 m MD (Costa & Davey, Reference Costa, Davey and Powell1992; Strauss et al., Reference Strauss, Elstner, du Chene, Mutterlose, Reiser and Brandt1993). In well HON-GT-02, this is based on the LOD of a single dinoflagellate cyst Stiphrosphaeridium dictyophorum.

Fig. 7. Combination GR logs, Palynological age dating of intervals and results of the SEG analysis. Age interpretation in VDB-GT-04 is based on Munsterman (2012). (A) Sandstone-rich zone of Valanginian age, (B) sandstone-rich zones of Early Valanginian/Late Ryazanian age.

Table 1. Overview of interval age dating.

The palynofacies and their relative occurrence in both HON-GT wells and PNA-GT-02 indicate that the Valanginian to Late Ryazanian interval was formed in a relatively humid, fluvial lowland environment, not directly positioned close to the coast. This last observation is derived from the absence or rare recognition of marine indicators. Both the sandstone content and the relative occurrence of sporomorphs associated to the ‘Lowland-dry’ eco-group is higher in HON-GT-02 compared to HON-GT-01, despite the relatively short distance of approximately 1.5 km between the wells in this doublet (Fig. 7).

In VDB-GT-04, the recognition of the Paratollia Ammonite Zone was based on the presence of Perisseiasphaeridium insolitum, Stiphrosphaeridium dictyophorum, Canningia compta, Hystrichosphaeridium scoriaceum and Oligosphaeridium diluculum in the samples at 1890 m and 1910 m MD (Costa & Davey, Reference Costa, Davey and Powell1992). A marine origin of the sample at 1890 m MD was recognised and therefore it may most likely be associated with the Paratollia MFS. Note that this is not based on SEG analysis in this well.

The results indicate that the aquifers in the four wells are not part of a single sandstone-rich succession. At least two sandstone-rich zones of Valanginian and Ryazanian age with limited lateral extent are encountered (Fig. 7). The Valanginian sandstone-rich zone A in HON-GT-01 relates to the upper section of the sandstone-rich zone in PNA-GT-02 with the same age. In contrast, the Valanginian succession in VDB-GT-04 is claystone-dominated. In this well, the aquifer is formed by a Ryazanian sandstone-rich zone B that relates to the Ryazanian sandstone-rich zone in PNA-GT-02. In HON-GT-01 the Ryazanian succession is claystone-dominated. Stacking of both sandstone-rich zones A and B accounts for the increased aquifer thickness in well PNA-GT-02. In contrast, lower aquifer thickness could be explained by the presence of a single sandstone-rich zone in HON-GT-01 and VDB-GT-04.

Regional well-log correlation

In the four wells with palynological analysis, the Paratollia MFS was recognised at approximately 300 m true vertical depth below the Elegans MFS (Fig. 8). In the wells without palynological cuttings analysis or in wells that did not reach sufficient depth, the top and base of the Valanginian succession were interpreted by extrapolation GR log patterns assuming constant thickness of the Valanginian interval. The resulting correlation scheme indicates that the Valanginian interval is sandstone-rich in both the Westland and Pijnacker grabens. In contrast, this interval has low sandstone content in the VDB fault block (Fig. 8). The Early Valanginian/Ryazanian interval has low sandstone content in the Westland fault block, but high sandstone content in the Pijnacker and VDB fault blocks. Our correlation scheme suggests that the prevailing position of the meander belts in which sand was deposited shifted from the east to the west side of the basin during the Ryazanian and Valanginian. In the Pijnacker fault block, both the Valanginian and the Early Valanginian/Ryazanian intervals have high sandstone content. This accounts for the largest interval with high sandstone content in all geothermal wells in the basin of approximately 250 m thickness in the PNA-GT-02 well. The other wells in this fault block have a limited total depth. Therefore they did not intersect the total Early Valanginian/Ryazanian interval. Similarly, the limited total depth of the VDB-GT-04 well only shows 70 m of the Early Valanginian/Ryazanian interval, of which approximately 50 m is sandstone-rich. In the Westland fault block, the Valanginian aquifer thickness ranges from 50 to 150 m. In the HON-GT-01 well, the lower half of the Valanginian interval is claystone-rich, unlike the other wells. The net-sandstone content (N/G) of the combined Valanginian and Early Valanginian intervals is calculated for each well. A specific GR log cut-off value is used in each well to take the differences in GR calibration into account. The N/G values range from 20% in the HON-GT-01 well where thick, non-aquifer intervals are included in the calculation, to 50% in the HAG-GT-01 well where the N/G calculation is based on the sandstone-dominated Valaginian interval. The arithmetic average N/G in all wells of the combined Valanginian and Early Valanginian / Late Ryazanian intervals is 35%.

Fig. 8. Well-log correlation of geothermal wells in three different fault blocks in the WNB. Solid lines indicate MFS interpretation based on cuttings analysis, dotted line is the projected MFS based on TVD, in wells without cuttings analysis.

Fig. 9. Series of maps indicating the location of sand-dominated sedimentation during the (A) Ryazanian / Early Valanginian and (B) Valanginian. Arrows indicate fluvial palaeoflow direction. (C) Cartoon illustrating facies distribution in the fluvial interval of the Nieuwerkerk Formation on a cross-section perpendicular to the fault trend.

Discussion

Ryazanian / Early Valanginian shift of sandstone-dominated facies

The seismic facies interpretation of a regional westward shift of sandstone-dominated facies during the Valanginian (Den Hartog Jager, Reference Den Hartog Jager, Rondeel, Batjes and Nieuwenhuijs1996) is corroborated by our palynology-based correlation (Fig. 9). Because of the unidirectional nature of the shift in sandstone-dominated facies in our study, we propose a tectonic origin of the shift. In the absence of such tectonic control, successive meander-belt avulsions and inherent compensational stacking would be the ruling processes, and a random spatial distribution of the fluvial sandstones the characteristic sedimentary architecture (e.g. Stouthamer & Berendsen, Reference Stouthamer and Berendsen2007; Hajek et al., Reference Hajek, Heller and Sheets2010; Donselaar et al., Reference Donselaar, Cuevas Gozalo and Moyano2013; Flood & Hampson, Reference Flood and Hampson2015; Van Toorenenburg et al., Reference Van Toorenenburg, Donselaar, Noordijk and Weltje2016). Based on this hypothesis, maps are generated that predict the lateral extent of the Ryazanian / Early Valanginian sandstone-rich zone (Fig. 9A) and the Valanginian sandstone-rich zone (Fig. 9B).

If tectonic movement had a strong impact on sedimentation, this might invalidate our assumption that the Valanginian interval has a constant thickness in our study area. In that case, it could be expected that the thickness of the Valanginian interval would increase in wells that are closer to the hanging wall of grabens. However, we expect that this would not have a significant impact on the trend in Figure 8 for two reasons. Firstly, the geothermal wells in the ‘Pijnacker Graben’ and ‘Westland Graben’ are drilled more or less parallel to the major fault trend and therefore the thickness correction would affect them equally. Secondly, in the wells with our palynological analysis an approximately constant thickness was observed, despite the fact that the four geothermal wells are located in different fault blocks. Additional palynological analysis in other WNB doublets could verify if the assumption is valid.

Regional aquifer architecture

In the present study, the correlation of the sandstone intervals in the Nieuwerkerk Formation is based on their chronostratigraphic position and the occurrence of two MFSs, whereas previous studies used a lithostratigraphic correlation. Van Adrichem Boogeart & Kouwe (Reference Van Adrichem Boogaert and Kouwe1993) identified the youngest fluvial sandstone-rich interval in the formation as the Delft Sandstone Member. Their regional aquifer architecture model (Fig. 2A) is used in current geothermal exploitation in the basin but does not adequately explain aquifer thickness variations like those shown in Figure 1. DeVault & Jeremiah (Reference DeVault and Jeremiah2002) described the regional aquifer architecture as a more random distribution of amalgamated sand complexes with limited lateral continuity that occur throughout the Nieuwerkerk Formation. These previous descriptions of the aquifer architecture were based on hydrocarbon wells on structural highs in the basin. Previously aquifer thickness prediction in the grabens was uncertain without well control, especially because these fault blocks might have experienced different tectonic movement. New well data from the graben fault blocks and the palynological analysis of our study suggest the aquifer architecture as sketched in Figure 9.

Because the number of geothermal wells in the grabens is still limited, the continuity of sandstone complexes is still uncertain. In the entire fluvial Valanginian to Late Ryazanian interval, N/G ranges from 20% to 50%, with an arithmetic average of 35% (Fig. 8). These percentages are an initial estimate of N/G, as no sensitivity study of GR cut-off value is included. Nevertheless this indicates that significant volumes of claystone are preserved. The Nieuwerkerk Formation is deposited by a relatively small meandering fluvial system with a palaeoflow depth of approximately 4 m (e.g. DeVault & Jeremiah, Reference DeVault and Jeremiah2002). The associated palaeochannel width and channel-belt width are therefore estimated to be approximately 40 m and 1–2 km, respectively (e.g. Bridge, Reference Bridge, Posamentier and Walker2006; Gibling, Reference Gibling2006). The maximum width of individual sandstone bodies is smaller than the channel-belt widths (e.g. Donselaar & Overeem, Reference Donselaar and Overeem2008; Donselaar et al., Reference Donselaar, Groenenberg and Gilding2015). Through amalgamation, sandstone-body width might extend further. However, claystone bodies are likely to form flow baffles or barriers perpendicular to the palaeoflow direction. This should be taken into account in doublet design and doublet placement as it will have an impact on possible interference between adjacent doublets and flowpath formation between injection- and production wells of individual doublets (e.g. Willems et al., Reference Willems, Nick, Donselaar, Weltje and Bruhn2017).

Our results have an impact on expected aquifer thickness in different fault blocks. Larger aquifer thickness could be expected in the Pijnacker fault block where sandstone-dominated zones in both successions overlap. Furthermore, our results affect expected aquifer depth and therefore drilling costs in different fault blocks. As shown in Figure 9C, the aquifer is found at greater depth in the VDB fault block than in the other two fault blocks. In addition, our results can be used for aquifer property extrapolation for new geothermal doublets. For example, the expected injectivity and productivity of future doublets should be based more on values which are measured in geothermal doublets in the same fault block. It is also possible that stratigraphically different sandstone successions have different properties. In current WNB doublets, productivity and injectivity vary considerably (Van Wees et al., Reference Van Wees, Degens, Zijp, De Boer, Obdam and Eyvazi2012). However, the variation could also be due to other factors such as scaling or skin formation. Van Wees et al. (Reference Van Wees, Degens, Zijp, De Boer, Obdam and Eyvazi2012) pointed out that unfortunately it is not possible to identify a single cause of this variability because of limited available data.

Palynological analyses and SEG method

Comparison of our results with those of Munsterman (Reference Munsterman2012) shows that the Valanginian interval in VDB-GT-04 has relatively more of a lower coastal plain character with respect to the HON-GT doublet and the PNA-GT-02 well. This could be due to a topographical difference between the fault blocks during the Valanginian. The SEG analysis indicated that in the HON-GT-02 well, a relatively higher fraction of ‘Lowland-dry’ type sporomorphs were recognised, which could point at a more inland location of the well compared to HON-GT-01 which is drilled more towards the palaeo-coastline in the north (e.g. Den Hartog Jager, Reference Den Hartog Jager, Rondeel, Batjes and Nieuwenhuijs1996). The Valanginian interval of HON-GT-01 also has a higher claystone content that could also be explained by a more near-coastal location of this well. In addition, the change in sandstone content could be explained by fault movement that directed the location of sand-rich meander belt deposits towards well HON-GT-02 while HON-GT-01 was located in the floodplain region. Due to the limited number of wells in our study, it is currently unclear how these observations are related to the fluvial architecture, sandstone distribution and the location of the palaeo-coastline.

The palynological age dating gives an indication of the age of interval and is not able to identify exact age boundaries or the exact location of the Paratollia MFS. The resolution of the age dating is limited by the sample spacing of 10 m and the risk of caving from higher sections. Uncertainty in age interpretation applies most to the Valanginian – Late Ryazanian boundary in our study area. Often, our interpretation of this boundary was based on recognition of a limited number of palynological indicators. In contrast, identification of the marine Elegans MFS has a lower degree of uncertainty because it was based on a combination of GR log interpretation and palynological analysis. This is because GR log signals in marine intervals are more often related to sea-level changes compared to GR log signals in fluvial intervals, like the Valanginian to Late Ryazanian interval in our study area. In addition, marine dinoflagellate cysts provide a higher-resolution dating than (long-ranging) terrestrial spores and pollen.

Our results underline the importance of palynological analysis for fluvial well-log correlation. These analyses enabled the identification of markers within fluvial claystone-dominated as well as sandstone-dominated successions, which would not have been possible based on GR log interpretation alone.

Conclusions

Based on the results of this study we can conclude that:

  • Current WNB geothermal doublets encounter sandstone-rich zones in at least two stratigraphic intervals of Valanginian age and of Early Valanginian/Ryazanian age.

  • Sandstone-rich zones in both intervals can overlap, which accounts for the large aquifer thickness in the PNA-GT-02 well.

  • Valanginian tectonic movement induced a shift of the deposition of sandstone-dominated facies from the east to the west of the basin.

  • This shift has an impact on expected aquifer thickness and aquifer depth in different fault blocks in the basin.

Acknowledgements

This study was carried out by Delft University of Technology in collaboration with Panterra Geoconsultants and TNO, in the context of the Research Agenda Geothermal Energy of the Ministry of Economic Affairs and LTO Glaskracht Nederland in the innovation programme ‘Kas als Energiebron’. We kindly thank the consortium of share- and stakeholders of the Delft Geothermal Project (DAP) for their support. We thank the associate editor Geert-Jan Vis and the reviewers Daan den Hartog Jager and Henk Kombrink for their constructive and helpful revisions of the original manuscript.

Supplementary material

Supplementary material is available online at https://doi.org/10.1017/njg.2017.23.

References

Abbink, O.A., 1998. Palynological investigations in the Jurassic of the North Sea region. PhD Thesis. Utrecht University (Utrecht).Google Scholar
Abbink, O.A., Targarona, J., Brinkhuis, H. & Visscher, H., 2001. Late Jurassic to earliest Cretaceous palaeoclimatic evolution of the southern North Sea. Global and Planetary Change 30: 231256.Google Scholar
Abbink, O.A., Van Konijnenburg-Van Cittert, J.H.A., Van der Zwan, C.J. & Visscher, H., 2004a. A sporomorph ecogroup model for the NW European Jurassic–Lower Cretaceous. Netherlands Journal of Geosciences 83 (2): 8192.Google Scholar
Abbink, O.A., Van Konijnenburg-Van Cittert, C.J. & Visscher, H., 2004b. A sporomorph ecogroup model for the NW European Jurassic–Lower Cretaceous: concepts and framework. Netherlands Journal of Geosciences 83 (1): 1738.Google Scholar
Bridge, J.S., 2006. Fluvial facies models: recent developments. In: Posamentier, H.W. & Walker, R.G. (eds.): Facies models revisited. SEPM Special Publication 84: 85170.Google Scholar
Costa, L.I. & Davey, R.J., 1992. Dinoflagellate cysts of the Cretaceous System. In: Powell, A.J. (ed.): A stratigraphic index of dinoflagellate cysts. Chapman & Hall (London/New York): 99154.Google Scholar
Davey, R.J., 1979. The stratigraphic distribution of dinocysts in the Portlandian (latest Jurassic) to Barremian (Early Cretaceous) of northwest Europe. AASP Contribution Series 5B: 4881.Google Scholar
Davey, R.J., 1982. Dinocyst stratigraphy of the latest Jurassic to Early Cretaceous of the Haldager No. 1 borehole, Denmark. Geological Survey of Denmark Series B, 6: 58 pp.Google Scholar
Den Hartog Jager, D.G., 1996. Fluviomarine sequences in the Lower Cretaceous of the West Netherlands Basin: correlation and seismic expression. In: Rondeel, H.E., Batjes, D.A.J. & Nieuwenhuijs, W.H. (eds): Geology of gas and oil under the Netherlands. Kluwer Academic Publishers (Dordrecht): 229241.Google Scholar
DeVault, B. & Jeremiah, J., 2002. Tectonostratigraphy of the Nieuwerkerk Formation (Delfland Subgroup), West Netherlands Basin. AAPG Bulletin 86 (10): 16791707.Google Scholar
Donselaar, M.E. & Overeem, I., 2008. Connectivity of fluvial point-bar deposits: an example from the Miocene Huesca fluvial fan, Ebro Basin, Spain. AAPG Bulletin 92 (9): 11091129.Google Scholar
Donselaar, M.E., Cuevas Gozalo, M.C. & Moyano, S., 2013. Avulsion processes at the terminus of low-gradient semi-arid fluvial systems: lessons from the Río Colorado, Altiplano endorheic basin, Bolivia. Sedimentary Geology 283: 114.Google Scholar
Donselaar, M.E., Groenenberg, R.M. & Gilding, D.T., 2015. Reservoir geology and geothermal potential of the Delft Sandstone Member in the West Netherlands Basin. World Geothermal Congress 2015, 19–25 April 2015, Melbourne, Australia: 9 pp. Conference proceedings.Google Scholar
Duin, E.J.T., Doornenbal, J.C., Rijkers, R.H.B., Verbeek, J.W. & Wong, Th.E., 2006. Subsurface structure of the Netherlands – results of recent onshore and offshore mapping. Netherlands Journal of Geosciences 85 (4): 245 276.Google Scholar
Duxbury, S., 2001. A palynological zonation scheme for the Lower Cretaceous United Kingdom Sector, Central North Sea. Neues Jahrbuch Geologisches Palaontologue Abhandlungen 219 (1): 95137.Google Scholar
Duxbury, S., Kadolsky, D. & Johansen, S., 1999. Sequence stratigraphic subdivision of the Humber Group in the Outer Moray Firth area (UKCS, North Sea). In: Jones, R.W. & Simmons, M.D. (eds): Biostratigraphy in production and development geology. Geological Society of London, Special Publication 152: 2354.Google Scholar
Flood, Y.S. & Hampson, G.J., 2015. Quantitative analysis of the dimensions and distribution of channelized fluvial sandbodies within a large outcrop dataset: Upper Cretaceous Blackhawk Formation, Wasatch Plateau, central Utah, U.S.A. Journal of Sedimentary Research 85 (4): 315336.Google Scholar
Gibling, M.R., 2006. Width and thickness of fluvial channel bodies and valley fills in the geological record: a literature compilation and classification. Journal of Sedimentary Research 76: 731770.Google Scholar
Gradstein, F., Ogg, J., Smith, A. & Ogg, G.M., 2012. A geologic time scale. Newsletters on Stratigraphy 45 (2): 171188.Google Scholar
Hajek, E.A., Heller, P.L. & Sheets, B.A., 2010. Significance of channel-belt clustering in alluvial basins. Geology 38 (6): 535538.Google Scholar
Heilmann-Clausen, C., 1987. Lower Cretaceous dinoflagellate biostratigraphy in the Danish Central Trough. Danmarks Geologische Undersogelse Serie A 17: 190.Google Scholar
Herngreen, G.F.W. & Wong, T.E., 2007. Cretaceous. In: Wong, Th., Batjes, D.A.J. & De Jager, J. (eds): Geology of the Netherlands. Royal Netherlands Academy of Arts and Sciences (Amsterdam): 127150.Google Scholar
Herngreen, G.F.W., Kerstholt, S.J. & Munsterman, D.K., 2000. Callovian – Ryazanian (‘Upper Jurassic’) palynostratigraphy of the Central North Sea Graben and Vlieland Basin, The Netherlands. Mededelingen Nederlands Instituut voor Toegepaste Geowetenschappen TNO 63.Google Scholar
Janssen, N.M.M. & Dammers, G., 2008. Sample processing for pre-Quaternary palynology. TNO report 2008-UR1190/A.Google Scholar
Jeremiah, J.M., Duxbury, S. & Rawson, P., 2010. Lower Cretaceous of the southern North Sea Basins: reservoir distribution within a sequence stratigraphic framework. Netherlands Journal of Geosciences 89 (3/4): 203237.Google Scholar
Munsterman, D.K., 2012. De resultaten van het palynologische onderzoek naar de ouderdom van de Onder Krijt successie in boring Van den Bosch-04 (VDB-04), interval 925–2006 m. TNO report TNO-060-UT-2011-02200/B.Google Scholar
NLOG, 2017. www.nlog.nl (Netherlands Oil and Gas portal).Google Scholar
Partington, M.A.P., Copestake, P., Mitchener, B.C. & Underhill, J.R., 1993. Biostratigraphic correlation of genetic stratigraphic sequences in the Jurassic – lowermost Cretaceous (Hettangian – Ryazanian) of the North Sea and adjacent areas. In: Parker, J. (ed.): Geological Society of London, Petroleum Geology Conference series 4: 371386.Google Scholar
Racero-Baena, A. & Drake, S., 1996. Structural style and reservoir development in the West Netherlands oil province. In: Rondeel, H., Batjes, D. & Nieuwenhuijs, W. (eds): Geology of gas and oil under the Netherlands. Kluwer Academic Publishers (Dordrecht): 211229.Google Scholar
Riding, J.B. & Thomas, J.E., 1992. Dinoflagellate cysts of the Jurassic System. In: Powell, A.J. (ed.): A stratigraphic index of dinoflagellate cysts. Kluwer Academic Publishers (Dordrecht): 798.Google Scholar
Stouthamer, E. & Berendsen, H.J.A., 2007. Avulsion, the relative roles of autogenic and allogenic processes. Sedimentary Geology 198 (3): 309325.Google Scholar
Strauss, C., Elstner, F., du Chene, R.J., Mutterlose, J., Reiser, H. & Brandt, K.-H., 1993. New micropaleontological and palynological evidence on the stratigraphic position of the ‘German Wealden’ in NW Germany. Zitteliana 20, 389401.Google Scholar
Van Adrichem Boogaert, H.A. & Kouwe, W.F.P., 1993. Stratigraphic nomenclature of the Netherlands, revision and update by Rijks Geologische Dienst (RGD) and Netherlands Oil and Gas Exploration and Production Association (NOGEPA). Mededelingen Rijks Geologische Dienst 50: 180 pp.Google Scholar
Van Toorenenburg, K.A., Donselaar, M.E., Noordijk, N.A. & Weltje, G.J., 2016. On the origin of crevasse-splay amalgamation in the Huesca fluvial fan (Ebro Basin, Spain): implications for connectivity in low net-to-gross fluvial deposits. Sedimentary Geology 343: 156164.Google Scholar
Van Wees, J.D.A.M., Degens, G., Zijp, M., De Boer, J., Obdam, A. & Eyvazi, F.J., 2012. BIA Geothermal – TNO umbrella report into the causes and solutions to poor well performance in Dutch geothermal projects. TNO report TNO 2012 R10719.Google Scholar
Vondrak, A., 2016. (Bio-) stratigraphic correlation of geothermal aquifers in the West Netherlands Basin. Panterra Geoconsultants report. Available at https://www.kasalsenergiebron.nl/duurzame-energie/aardwarmte.Google Scholar
Willems, C.J.L., Nick, H.M., Donselaar, M.E., Weltje, G.J. & Bruhn, D.F., 2017. On the connectivity anisotropy in fluvial Hot Sedimentary Aquifers and its influence on geothermal doublet performance. Geothermics 65: 222233.Google Scholar
Figure 0

Fig. 1. Gamma-ray logs of two geothermal doublets HON-GT and PNA-GT. Fault interpretation is based on Duin et al. (2006). Well-log correlation is based on ‘End-of-well reports’ (NLOG, 2017).

Figure 1

Fig. 2. Stratigraphic column for the Early Cretaceous section in the WNB indicating tectonic activity during deposition of the Rijnland Group, the Nieuwerkerk Formation and the main geothermal aquifers in the WNB: the Rijswijk Sandstone Member and the Delft Sandstone Member (Van Adrichem Boogaert and Kouwe, 1993).

Figure 2

Fig. 3. (A) Cartoons illustrating the difference in sandstone distribution in the Nieuwerkerk Formation in schematic strike sections on graben scale according to (A) the Delft Sandstone model and (B) the multiple sandstone-rich zones models.

Figure 3

Fig. 4. Location of the geothermal wells for the GR well-log correlation, cuttings analysis, the outline of the seismic cross-section and the regional structural interpretation by Duin et al. (2006).

Figure 4

Fig. 5. Schematic representation of the impact of (A) low sea level and (B) high sea level on the relative occurrence of eco-groups: 1=Upland, 2=Lowland, 3=River, 4=Pioneer, 5=Coastal, 6=Tidally influenced. Note the lower relative occurrence of the ‘Lowland’ eco group with high sea level in (B). Modified from Abbink et al. (2004a,b).

Figure 5

Fig. 6. (A) Seismic section with interpretation of faults as well as top (yellow) and base (pink) of the Nieuwerkerk Formation horizons. (B) Horizon flattening on top Nieuwerkerk Formation. The outline of the seismic section and the outline of the interpreted faults are indicated on the map.

Figure 6

Fig. 7. Combination GR logs, Palynological age dating of intervals and results of the SEG analysis. Age interpretation in VDB-GT-04 is based on Munsterman (2012). (A) Sandstone-rich zone of Valanginian age, (B) sandstone-rich zones of Early Valanginian/Late Ryazanian age.

Figure 7

Table 1. Overview of interval age dating.

Figure 8

Fig. 8. Well-log correlation of geothermal wells in three different fault blocks in the WNB. Solid lines indicate MFS interpretation based on cuttings analysis, dotted line is the projected MFS based on TVD, in wells without cuttings analysis.

Figure 9

Fig. 9. Series of maps indicating the location of sand-dominated sedimentation during the (A) Ryazanian / Early Valanginian and (B) Valanginian. Arrows indicate fluvial palaeoflow direction. (C) Cartoon illustrating facies distribution in the fluvial interval of the Nieuwerkerk Formation on a cross-section perpendicular to the fault trend.

Supplementary material: File

Willems et al. supplementary material 1

Willems et al. supplementary material

Download Willems et al. supplementary material 1(File)
File 27.2 KB