• Energy Research
• Embargo close
• NL close
• Delft University of Technology close

Periglacial alluvial fans are common in northwestern and central Europe and their pre-Holocene stratigraphic records typically date back to late Middle Pleniglacial and Late Pleniglacial (late MIS3 and 2). Preserved stratigraphic records that include an entire interglacial-glacial cycle have, so far, not been described and it is thus unknown how periglacial alluvial fans responded during a full cycle of interglacial-glacial climate changes. In this paper, we reconstruct the evolution of the Eerbeek periglacial alluvial fan in the Netherlands which was deposited during the late Saalian (MIS 6) to late Weichselian (MIS 2) period, including the entire last interglacial–glacial cycle (MIS 5-2). Our reconstruction is based on 48, up-to 45-m deep borehole and Cone Penetration Test (CPT) logs that allowed the construction of an 8-km long longitudinal and a 7-km long transverse cross section over the Eerbeek periglacial alluvial fan. Age control was provided by means of 17, previously published, Optically Stimulated Luminescence ages of two boreholes on the fan, and 14 14C ages from three boreholes and a nearby, now abandoned, quarry.Overlying a thick, late Saalian (MIS 6) alluvial fan record, is a 4- to 18-m thick alternation of distinct organic (mainly peat and humic clays), siliciclastic alluvial fan (coarse- and medium-grained sands), Rhine (coarse- and medium grained sands), and aeolian (mainly medium-grained sands) stratigraphic units. Organic levels indicate fan stability during the Eemian interglacial (MIS 5e), and Brørup (MIS 5c), Odderade–Ognon interstadial complex (MIS 5a), and Middle Pleniglacial (MIS 3) interstadials 14, 13, 12 and 11 as well as late MIS 2 interstadial 1a. Clastic sediments indicate alluvial fan activity during the Herning (MIS 5d), Rederstall (MIS 5b), Ognon stadial complex (late MIS 5a), Early Pleniglacial (MIS 4) and upper Middle Pleniglacial (upper MIS 3) stadials 13, 12 and 11. Sediments from the coldest and driest period of the Last Glacial (late MIS 3 and MIS 2) are absent and following a phase of aeolian activity, the fan was only reactivated at the MIS 2 to MIS 1 transition (stadial 1). We attribute the absence of fan activity during the coldest period of the last interglacial-glacial cycle to the eastward orientation of the fan making it less sensitive to permafrost melt.The colder MIS substages and stadials in which the Eerbeek fan was active coincided with the presence of permafrost and/or a seasonal, deeply frozen soil, and a relatively humid climate during which vegetation was largely absent. The presence of channels that dissect the underlying organic units suggests that the Eerbeek fan initially responded to the changes from interstadials to stadials by means of erosion. As climate cooled and permafrost/deep frost developed, the fan switched to alluvial aggradation. The consistent presence of coarsening-fining upward sequences suggests a relation with cycles of increased overland flow due to increasingly more frozen subsoil conditions. The fan stratigraphy therefore shows the direct coupling between warmer-colder MIS substages and interstadial-stadial climate cyclicity and alluvial fan response over the entire last interglacial-glacial cycle.

Periglacial alluvial fans are common in northwestern and central Europe and their pre-Holocene stratigraphic records typically date back to late Middle Pleniglacial and Late Pleniglacial (late MIS3 and 2). Preserved stratigraphic records that include an entire interglacial-glacial cycle have, so far, not been described and it is thus unknown how periglacial alluvial fans responded during a full cycle of interglacial-glacial climate changes. In this paper, we reconstruct the evolution of the Eerbeek periglacial alluvial fan in the Netherlands which was deposited during the late Saalian (MIS 6) to late Weichselian (MIS 2) period, including the entire last interglacial–glacial cycle (MIS 5-2). Our reconstruction is based on 48, up-to 45-m deep borehole and Cone Penetration Test (CPT) logs that allowed the construction of an 8-km long longitudinal and a 7-km long transverse cross section over the Eerbeek periglacial alluvial fan. Age control was provided by means of 17, previously published, Optically Stimulated Luminescence ages of two boreholes on the fan, and 14 14C ages from three boreholes and a nearby, now abandoned, quarry.Overlying a thick, late Saalian (MIS 6) alluvial fan record, is a 4- to 18-m thick alternation of distinct organic (mainly peat and humic clays), siliciclastic alluvial fan (coarse- and medium-grained sands), Rhine (coarse- and medium grained sands), and aeolian (mainly medium-grained sands) stratigraphic units. Organic levels indicate fan stability during the Eemian interglacial (MIS 5e), and Brørup (MIS 5c), Odderade–Ognon interstadial complex (MIS 5a), and Middle Pleniglacial (MIS 3) interstadials 14, 13, 12 and 11 as well as late MIS 2 interstadial 1a. Clastic sediments indicate alluvial fan activity during the Herning (MIS 5d), Rederstall (MIS 5b), Ognon stadial complex (late MIS 5a), Early Pleniglacial (MIS 4) and upper Middle Pleniglacial (upper MIS 3) stadials 13, 12 and 11. Sediments from the coldest and driest period of the Last Glacial (late MIS 3 and MIS 2) are absent and following a phase of aeolian activity, the fan was only reactivated at the MIS 2 to MIS 1 transition (stadial 1). We attribute the absence of fan activity during the coldest period of the last interglacial-glacial cycle to the eastward orientation of the fan making it less sensitive to permafrost melt.The colder MIS substages and stadials in which the Eerbeek fan was active coincided with the presence of permafrost and/or a seasonal, deeply frozen soil, and a relatively humid climate during which vegetation was largely absent. The presence of channels that dissect the underlying organic units suggests that the Eerbeek fan initially responded to the changes from interstadials to stadials by means of erosion. As climate cooled and permafrost/deep frost developed, the fan switched to alluvial aggradation. The consistent presence of coarsening-fining upward sequences suggests a relation with cycles of increased overland flow due to increasingly more frozen subsoil conditions. The fan stratigraphy therefore shows the direct coupling between warmer-colder MIS substages and interstadial-stadial climate cyclicity and alluvial fan response over the entire last interglacial-glacial cycle.

The aim of this study was to analyze whether, and if so, in what way risks would influence the design,costs and routing of CO2pipelines. This article assesses locational and societal risks of CO2pipelinetransport and analyses whether rerouting or implementing additional risk mitigation measures is themost cost-effective option. The models EFFECTS and RISKCURVES are used to estimate the dispersion andrisk, respectively. The pipeline routes are optimized by using the least cost path function in ArcGIS.This article evaluates three case studies in the Netherlands. The results show that pipelines transportingdense phase CO2(8–17 MPa) with a minimal amount of risk mitigation measures already meet the 10−6locational risk required in the Netherlands. 10−6locational risks of 135 m are calculated for intermediatepumping stations, handling 450 kg CO2/s (about 14 Mt CO2/year). In all the cases, pumping stations couldbe located along the pipeline route without any problem.For the cases studied transporting gaseous CO2(1.5–3 MPa) leads to larger 10−6locational risk distancesthan transporting dense phase CO2. This is caused by the large momentum behind a dense phase CO2release, leading to smaller but higher jet and to a higher mixing rate with the surrounding air than for agaseous CO2release.Based on our analysis, it can be concluded that dense phase CO2transport is safe if it is well organized.The risks are manageable and widely accepted under current legislation. In addition, risk mitigationmeasures, like marker tape and increased surveillance, are available which reduce the risk significantlyand increase the costs only slightly. Pipeline routing for gaseous CO2transport appears more challengingin densely populated areas, because larger safety zones are attached to it.

The aim of this study was to analyze whether, and if so, in what way risks would influence the design,costs and routing of CO2pipelines. This article assesses locational and societal risks of CO2pipelinetransport and analyses whether rerouting or implementing additional risk mitigation measures is themost cost-effective option. The models EFFECTS and RISKCURVES are used to estimate the dispersion andrisk, respectively. The pipeline routes are optimized by using the least cost path function in ArcGIS.This article evaluates three case studies in the Netherlands. The results show that pipelines transportingdense phase CO2(8–17 MPa) with a minimal amount of risk mitigation measures already meet the 10−6locational risk required in the Netherlands. 10−6locational risks of 135 m are calculated for intermediatepumping stations, handling 450 kg CO2/s (about 14 Mt CO2/year). In all the cases, pumping stations couldbe located along the pipeline route without any problem.For the cases studied transporting gaseous CO2(1.5–3 MPa) leads to larger 10−6locational risk distancesthan transporting dense phase CO2. This is caused by the large momentum behind a dense phase CO2release, leading to smaller but higher jet and to a higher mixing rate with the surrounding air than for agaseous CO2release.Based on our analysis, it can be concluded that dense phase CO2transport is safe if it is well organized.The risks are manageable and widely accepted under current legislation. In addition, risk mitigationmeasures, like marker tape and increased surveillance, are available which reduce the risk significantlyand increase the costs only slightly. Pipeline routing for gaseous CO2transport appears more challengingin densely populated areas, because larger safety zones are attached to it.

In this paper we present a review of controlling geological, tectonic and engineering factors for induced seismicity associated to gas depletion in the Netherlands and we place experiences from extensive Dutch geomechanical studies in the past decade in the context of generic models for induced seismicity. Netherlands is in a mature gas production phase, marked by excellent subsurface structural and stratigraphic characterization. Over 190 gas fields of varying size have been exploited. No more than 15% of these fields show seismicity. Geomechanical studies show that, similar to the EGS stimulation phase, largest seismicity is localized on pre-existing fault structures. However, the prime cause for seismicity in gas depletion is differential compaction, whereas in EGS stimulation related pressure build-up and fluid pressure diffusion along the faults form the prime mechanism. On the other hand, our study has a close theoretical analogy to reservoirs where the fluid volumes extracted are significantly larger than the re-injected volumes, and which can result in (differential) reservoir compaction.The observed onset of induced seismicity in the Netherlands occurs after a considerable pressure drop in the gas fields. Geomechanical models show that both the delay in the onset of induced seismicity as well as the non-linear increase in seismic moment observed in the induced seismicity, can be explained using a model of differential compaction, if the faults involved in induced seismicity are not critically stressed at the onset of depletion. The presented model serves to highlight key aspects of the interaction of initial stress and differential compaction in the framework of induced seismicity in Dutch gas fields. It is not intended as predictive model for induced seismicity in a particular field. To this end, a much more detailed field specific study, taking into account the full complexity of reservoir geometry, depletion history, mechanical properties and initial stress field conditions is required.

In this paper we present a review of controlling geological, tectonic and engineering factors for induced seismicity associated to gas depletion in the Netherlands and we place experiences from extensive Dutch geomechanical studies in the past decade in the context of generic models for induced seismicity. Netherlands is in a mature gas production phase, marked by excellent subsurface structural and stratigraphic characterization. Over 190 gas fields of varying size have been exploited. No more than 15% of these fields show seismicity. Geomechanical studies show that, similar to the EGS stimulation phase, largest seismicity is localized on pre-existing fault structures. However, the prime cause for seismicity in gas depletion is differential compaction, whereas in EGS stimulation related pressure build-up and fluid pressure diffusion along the faults form the prime mechanism. On the other hand, our study has a close theoretical analogy to reservoirs where the fluid volumes extracted are significantly larger than the re-injected volumes, and which can result in (differential) reservoir compaction.The observed onset of induced seismicity in the Netherlands occurs after a considerable pressure drop in the gas fields. Geomechanical models show that both the delay in the onset of induced seismicity as well as the non-linear increase in seismic moment observed in the induced seismicity, can be explained using a model of differential compaction, if the faults involved in induced seismicity are not critically stressed at the onset of depletion. The presented model serves to highlight key aspects of the interaction of initial stress and differential compaction in the framework of induced seismicity in Dutch gas fields. It is not intended as predictive model for induced seismicity in a particular field. To this end, a much more detailed field specific study, taking into account the full complexity of reservoir geometry, depletion history, mechanical properties and initial stress field conditions is required.

We developed a coupled code to obtain a better understanding of the role of pore pressure changes in causing fracture reactivation and seismicity during EGS. We implemented constitutive models for fractures in a continuum approach, which is advantageous because of the ease of integration in existing geomechanical codes (FLAC3D), the speed of the calculations and the flexibility of the fracture representation. For modeling the mechanical behavior of the fracture zone the softening ubiquitous joints model was used with random strength properties for the fracture zone. We implemented a hyperbolic deformation with effective stress for the reversible tensile fracture opening, and a linear relationship between plastic shear strain and ir eversible fracture opening. The effective permeability of the fracture network was described by a cubic dependence on the fracture aperture. We created a model inspired on the Soultz-la-Forêts GPK3 injection well, intersected by a dominant fracture zone. Due to injection of water, the model reproduced reactivation of the fracture zone and we observed the growth of a zone with large directional permeability. The model was used to perform a sensitivity analysis on parameters like injection rates, in-situ stress regime, fracture strength, and frictional weakening. This allowed us to evaluate the trends of their impact on fracture reactivation, including reactivated area, seismic moment and moment magnitudes.