AM

A.W. Martinius

info

Please Note

42 records found

Journal article (2026) - Lucie Alain, Christian Haug Eide, Henrik Nygaard Hansen, Agustin Argüello Scotti, Monika Shukla, Bjørn Terje Oftedal, Ian Sharp, Allard Martinius
Chlorite-coated intervals occurring in deeply buried sandstone reservoirs preserve porosity and permeability by inhibiting quartz cementation. However, the geometry of such chlorite-coated intervals and the controls on their distribution remain unclear. The Lower Jurassic Tilje Formation of the Halten Terrace, Norwegian Continental Shelf, shows sparsely distributed chlorite coats preserving high permeability at great depth. Changes between chlorite-coated and quartz-cemented intervals can occur at a centimetre scale and are thus hard to predict. In this study, detailed sedimentological investigations of 12 wells in the Smørbukk Field, integrated with biostratigraphic data and combined with permeability analysis of more than 5000 core plug samples, facilitate the interpretation of depositional elements and the tracing of chlorite-coated intervals across the field. Ten depositional environments are recognised, representing deposition as a large tide-dominated delta that prograded in an embayment for most of the Tilje Formation. A shift in depositional processes occurs towards the top of the formation, with the development of a latest Pliensbachian candidate relative sea-level fall and incision rapidly followed by the onset of a latest Pliensbachian–earliest Toarcian transgressive phase. The transgressive phase is associated with a change to an open-marine seaway and the deposition of marine tidal dunes. Two types of chlorite-coated intervals occur: (i) thick (20 to 30 m), stratigraphically bounded intervals and (ii) thin (<3 m), sparsely distributed intervals. Chlorite-coated intervals occur mostly in fluvial-tidal channels but also occasionally in other depositional environments, suggesting that sedimentological controls do not fully explain their distribution. Instead, regional controls such as the composition of the hinterland, climate, relative sea-level change, tectonics and geochemical conditions in the primary depositional environments are considered important factors. This study illustrates the various geometries of chlorite-coated intervals in the subsurface, having important implications for building predictive models of their distribution and for conventional and CO2 reservoir characterisation. ...
Conference paper (2026) - Q. Zhang, S. Geiger, J.E. Storms, M.D. Jackson, C. Jacquemyn, G.J. Hampson, A.W. Martinius
Stratigraphic and diagenetic heterogeneities exert a strong control on CO₂ plume dynamics and long-term storage performance in shallow-marine reservoirs. Using multiphase flow simulations conditioned to a geologically realistic reservoir model, this study demonstrates that repeated cemented barriers and stratigraphic baffles significantly slow vertical plume ascent and promote lateral plume spreading. This behavior increases overall CO₂ storage through a combined action of stratigraphic trapping, capillary pinning, residual trapping, and dissolution. Plume dispersion and local spreading were found to strongly enhance dissolution trapping. By subdividing a single buoyant plume into smaller, vertically confined accumulations, stratigraphic barriers increase CO₂–brine interfacial area and prolong residence times, resulting in substantially greater solubility trapping than would occur in homogeneous reservoir models. Structural configuration further modulates this behavior, with reduced buoyant drive leading to higher dissolved fractions. The strong sensitivity of plume migration and trapping behavior to geological heterogeneity highlights the importance of fitness-for-purpose modeling in CCS site assessment. While simplified models may be appropriate for early-stage screening, they may fail to capture critical plume dispersion and trapping processes if applied beyond their intended scope. These findings emphasize the need for detailed geological characterization and advanced modeling approaches tailored specifically to CCS applications. ...
Conference paper (2026) - C. Jacquemyn, M.D. Jackson, G.J. Hampson, D. Petrosvkyy, Q. Zhang, J. Storms, S. Geiger, A.W. Martinius
Assessment of CO2 storage capacity in shallow marine reservoirs must be underpinned by understanding of the impact of geological heterogeneity (Zhang et al., 2025). Heterogeneity is present at multiple scales from laminae (cm’s) to facies (m’s), to facies associations (10’s m), to parasequences (100’s m), to parasequence set (km’s) scale. Using a multiscale representative elementary volume (REV) approach, we investigate how heterogeneity at the facies and facies association scales impact properties such as effective permeability and effective relative permeability that are key inputs for simulating storage at reservoir scale. Facies-scale (m’s) heterogeneity clearly impacts effective permeability values and effective relative permeability curves, and their anisotropy. Models that capture the heterogeneity architecture of facies in 3D form a robust basis to link plug-scale (cm’s) measurements to larger scales, providing values that are potentially directly usable in flow simulation assessing CO2 flow dynamics and trapping for CCS. Upscaled relative permeability curves are anisotropic and can vary outside the bounds of the input curves as a direct result of 3D heterogeneity. This finding demonstrates that sedimentological heterogeneity affects CO2 flow and storage and illustrates how core-plug-scale measurements can be used to make predictions that incorporate the effects of facies-scale heterogeneity at larger scales. ...
Reducing the uncertainty of reservoir characterization requires to better identify the small-scale structures of the subsurface from the available data. Studying the seismic response of meter-scale, stratigraphic heterogeneities typically relies on the generation of reservoir models based on outcrop examples and their forward seismic modelling. To bridge geological information and seismic modelling, these methods allocate values of acoustic properties, such as mass-density and P-wave velocity, according to discretized properties like layer-type lithology or facies units. This strategy matches the current workflow in seismic data inversion in industry, where modelling workflows are based on lithofacies distributions. However, from stratigraphic modelling, we know that meter-scale heterogeneities occur within certain facies and lithologies. Here, we evaluate the difference on the seismic response between allocating acoustic properties in a grain size–based, semi-continuous manner versus discretized manners based on lithology and facies classifications. To do so, we generate a reference geological simulation that we populate with acoustic properties, mass-density and P-wave velocity, using three different strategies: (1) based on grain size distribution; (2) based on facies distribution; and (3) based on lithology. The method we propose includes the generation of realistic geological simulations based on stratigraphic modelling and the transformation of its output into acoustic properties, honouring the intra-lithology and intra-facies, small-scale structures. We, then, generate seismic data by applying a forward seismic modelling workflow. The synthetic data show that the grain size–based simulation allows the identification of small-scale, stratigraphic heterogeneities, such as beds with strong density and velocity contrasts. These stratigraphic structures are smoothened or may completely disappear in the facies and lithology discretized simulations and, therefore, are not (well) represented in the synthetic seismic data. Recognizing meter-scale, stratigraphic heterogeneities is relevant for the characterization of the fluid flow in the reservoir. However, current discrete and lithology-based strategies in seismic inversion are not able to resolve such heterogeneities because real subsurface properties are not discrete properties but continuous, unless there are stratigraphic discontinuities such as erosional surfaces or faults. This research works towards a better understanding of the relationship between changes in these continuous properties and the observed seismic data by introducing greater complexity into the discretized geological simulations. Here, we use synthetic seismic images with the goal of eventually aiding in fine-tuning seismic inversion methodologies applied to real seismic data. One pathway is to foster the development of inversion approaches that can leverage stratigraphic modelling to get stronger geological priors and replace the standard but inadequate multi-Gaussian prior. ...

Mechanisms, terminology and State-of-the-Art

Review (2025) - Qin Zhang, Sebastian Geiger, Joep E.A. Storms, Denis V. Voskov, Matthew D. Jackson, Gary J. Hampson, Carl Jacquemyn, Allard W. Martinius
Capillary pinning refers to the immobilization of CO₂ at capillary barriers when the uprising CO2 pressure is lower than the capillary entry pressure of the overlaying pore throats. Also known as local capillary trapping, it has been proposed as a fifth geologic CO₂ storage mechanism, alongside structural, solubility, residual, and mineral trapping. Despite extensive research, the fragmented terminology surrounding capillary pinning has led to confusion, making it challenging to synthesize findings effectively. Often conflated with mechanisms such as residual and hysteresis trapping, capillary pinning is commonly underestimated or completely overlooked in reservoir-scale models. Furthermore, difficulties in characterizing and upscaling small-scale geologic heterogeneities that influence capillary pinning contribute to significant uncertainties, with estimates of CO₂ trapped via this mechanism ranging from 3 % to 100 % of total CO₂ trapped via capillary actions. This review explores the fundamental mechanisms, experimental findings, and modeling approaches for assessing CO₂ capillary pinning in carbon capture and storage (CCS). It seeks to bridge the gap between the reservoir engineering community, with its extensive expertise in hydrocarbon recovery but that needs adjustments for CCS applications, and the subsurface storage community, which stands to benefit from this knowledge but often lacks access to relevant literature. Additionally, the study identifies key research opportunities to advance the understanding of capillary pinning in sedimentary rocks, ultimately enhancing the efficacy and reliability of CCS operations. ...
Journal article (2025) - Thomas E. Nichols, Richard H. Worden, James E. Houghton, Joshua Griffiths, Christian Brostrøm, Allard W. Martinius
We have developed a generalisable machine learning framework for reservoir quality prediction in deeply buried clastic systems. Applied to the Lower Jurassic deltaic sandstones of the Tilje Formation (Halten Terrace, North Sea), the approach integrates sedimentological facies modelling with mineralogical and petrophysical prediction in a single workflow. Using supervised Extreme Gradient Boosting (XGBoost) models, we classify reservoir facies, predict permeability directly from standard wireline log parameters and estimate the abundance of porosity-preserving grain coating chlorite (gamma ray, neutron porosity, caliper, photoelectric effect, bulk density, compressional and shear sonic, and deep resistivity). Model development and evaluation employed stratified K-fold cross-validation to preserve facies proportions and mineralogical variability across folds, supporting robust performance assessment and testing generalisability across a geologically heterogeneous dataset. Core description, point count petrography, and core plug analyses were used for ground truthing. The models distinguish chlorite-associated facies with up to 80% accuracy and estimate permeability with a mean absolute error of 0.782 log(mD), improving substantially on conventional regression-based approaches. The models also enable prediction, for the first time using wireline logs, grain-coating chlorite abundance with a mean absolute error of 1.79% (range 0–16%). The framework takes advantage of diagnostic petrophysical responses associated with chlorite and high porosity, yielding geologically consistent and interpretable results. It addresses persistent challenges in characterising thinly bedded, heterogeneous intervals beyond the resolution of traditional methods and is transferable to other clastic reservoirs, including those considered for carbon storage and geothermal applications. The workflow supports cost-effective, high-confidence subsurface characterisation and contributes a flexible methodology for future work at the interface of geoscience and machine learning. ...

Storage Potential and Impacts of Heterogeneity in Pressure Front

Conference paper (2025) - F. Silva Lira, M. Erdtmann, G. Gantois, F. Da Costa, H. Zerfass, L. Cassino, P. Walter, A. Martinius, S. Geiger, L. Menezes, S. Bortolini, A. Guirro, G. Waisman, F. Bulhoes, G. Freitas, G. Vieira, N. Lima
Brazil’s industrial emissions are 180 million tons of CO2 per year, and approximately 60% of these emissions are coming from industrial clusters located in the southeast. Tthe development of new offshore storage locations in this region is hence of strategic importance for future Carbon Capture and Storage (CCS) projects in Brazil. This study presents an evaluation of CO2 storage in the deep saline aquifers of the Jureia-Ponta Aguda formation, a gigaton-scale storage resource located in the shallow waters of Santos Basin that has the potential to support the development of at least three larger CCS hubs, each with a target injection of 1Gt of CO2. We show how Lorenz Coefficient map allows us to screen the pressure influence areas for each hub by linking reservoir heterogeneity to the spatio-temporal evolution of the pressure front, thereby identifying potential risks of pressure interference between neighboring CCS hubs. ...
Conference paper (2025) - Q. Zhang, S. Geiger, J. Storms, H. Hajibeygi, M. Jackson, G. Hampson, C. Jacquemyn, S. Krevor, A. Martinius
The North Sea’s potential as a Green Energy Hub depends on large-scale CO2 storage in shallow-marine sandstones, but the effects of geologic heterogeneity, such as permeability barriers and capillary entry pressure contrasts, remain underexplored. This study uses multiphase flow simulations on geologically realistic, surface-based reservoir models informed by outcrop analogue data from wave-dominated shoreface sandstones. We investigate how sedimentological heterogeneity influences CO2 plume migration, pressure evolution, and storage capacity.

Preliminary results show that capillary barriers tied to facies architecture and early cementation, conditioned to clinoform geometries, significantly control plume movement. These barriers promote lateral spreading and residual trapping, representing a potential upper limit on long-term CO2 storage when stable. Clinoform-related heterogeneity also induces flow compartmentalization, limiting pressure dissipation and enhancing anisotropy, which may reduce injectivity and cause spatially variable pressure buildup.

Comparisons with waterflood simulations reveal contrasting dynamics: water advances more uniformly, while CO2 migration is more sensitive to fine-scale architecture due to its lower interfacial tension and capillary entry pressures. These findings underscore the need to incorporate realistic sedimentological heterogeneity in dynamic models to avoid misestimating injectivity, pressure behavior, and storage security. This approach offers a robust framework for early-stage screening and risk assessment in complex storage settings. ...
Journal article (2025) - E. Cecchetti, A.W. Martinius, M. Felder, M.E. Donselaar, H.A. Abels
The Main Buntsandstein Subgroup in the Roer Valley Graben in the southern Dutch subsurface is a sand-prone sedimentary interval deposited in a fluvial-aeolian environment, and is currently investigated for its suitability as target for low-entalpy geothermal exploration. The current depositional models in the Roer Valley Graben do not fully address the facies heterogeneities within and between Buntsandstein sedimentary units and their impact on the prediction of reservoir architecture. A detailed analysis of the Main Buntsandstein sedimentary facies heterogeneities to de-risk future sustainable energy operations is therefore crucial. In the present study, the sedimentology and lithostratigraphy of the Buntsandstein are assessed in a multidisciplinary analysis by use of a subsurface dataset composed of well cores, gamma-ray logs, and thin section data. The deposition of the Main Buntsandstein sediments in the Roer Valley Graben is dominated by different fluvial processes, with minor aeolian reworking. River planform style evolved through geological time from highly mobile and ephemeral to more perennial in nature. These changes in river style seem to be dictated by a decrease in climatic aridity along with a decrease in tectonic activity. The depositional processes resulted in the development of six lithofacies associations, developing three different types of reservoir architectures with their own set of heterogeneities at different spatial scales. Amalgamated, stacked sandstones have the highest net-to-gross (N/G) with a high degree of lateral and vertical connectivity, and the highest average porosity and permeability. Compensational-stacked sandstone reservoirs have a lower N/G and are the most heterogeneous due to the frequent occurrence of cemented intervals as well as mud drapes in the sandstone bodies. Marginal isolated sandstones show a well preserved relationship between reservoir properties and depositional facies, while more data are needed to resolve the spatial connectivity and lateral continuity of these sandstone bodies. The results of this study enhance the understanding of Lower Triassic reservoir architecture and sedimentary heterogeneities in the Roer Valley Graben that can be applied well beyond the area and provide a solid basis for future investigation of the relationship between sedimentary facies, diagenesis, and reservoir quality. ...
Journal article (2024) - Youwei Wang, Timothy F. Baars, Joep E.A. Storms, Allard W. Martinius, Philip D. Gingerich, Magda Chmielewska, Simon J. Buckley, Hemmo A. Abels
Sedimentation on river floodplains is a complex process that involves overbank flooding, crevasse splaying, and river avulsion. The resulting floodplain stratigraphy often exhibits floodplain aggradation cycles with alternating fine-grained overbank flooding deposits that underwent significant petrogenesis, and coarser-grained, avulsion-belt deposits largely devoid of pedogenic impact. These cycles are linked to lateral migration and avulsion of channels driven by internal dynamics, external factors, or a combination of both. To better understand the spatial and vertical variability of such floodplain aggradation cycles, we map these in three dimensions using a photogrammetric model of the lower Eocene Willwood Formation in the northern Bighorn Basin, Wyoming, USA. This allows identifying 44 floodplain aggradation cycles in ∼300 m of strata with an average thickness of 6.8 m and a standard deviation of 2.0 m. All the cycles are traceable over the entire model, pointing to their spatial consistency over the 10 km2 study area. At the same time, rapid lateral thickness changes of the floodplain aggradation cycles occur with changes up to 4 m over a lateral distance of 400 m. Variogram analyses of both field and numerical-model results reveal stronger consistency of floodplain aggradation cycle thicknesses along the paleoflow direction compared to perpendicular to paleoflow. Strong compensational stacking occurs at the vertical scale of 2–3 floodplain aggradation cycles (14–20 m), while full compensational stacking occurs at larger scales of more than six floodplain aggradation cycles (>41 m). The lateral and vertical thickness variability of the floodplain aggradation cycles, as well as their compensational stacking behavior, are interpreted to be dominantly driven by autogenic processes such as crevasse splaying and avulsing that preferentially fill topographic lows. External climate forcing may have interacted with these autogenic processes, producing the laterally persistent and vertically repetitive floodplain aggradation cycles. The spatial variability of floodplain aggradation cycles demonstrated in this study highlights again the need for three-dimensional data collection in alluvial floodplain settings rather than depending on one-dimensional records. ...
Journal article (2024) - Youwei Wang, Timothy F. Baars, Joep E.A. Storms, Allard W. Martinius, Philip D. Gingerich, Hemmo A. Abels
Alluvial stratigraphy builds up over geologic time under the complex interplay of external climatic and tectonic forces and internal stochastic processes. This complexity makes it challenging to attribute alluvial stratigraphic changes to specific factors. Geological records indicate pronounced and persistent climatic changes during the Phanerozoic, while the effects of these changes on alluvial stratigraphy remain insufficiently documented. We provide evidence for 405 k.y. long-eccentricity climate forcing of alluvial stratigraphy in the lower Eocene Willwood Formation of the Bighorn Basin, Wyoming (USA). Two ∼90-m-thick intervals, characterized by a relative paucity of sand, dominance of sinuous-river channels, and floodplain sediments with better-developed paleosols, coincide with eccentricity maxima as determined through integrated stratigraphic methods. These intervals are interspersed with three contrasting intervals, marked by relatively high sand content, prevalent braided-river channels, and less-developed paleosols, corresponding to eccentricity minima. A comprehensive genetic model that integrates climate, source-to-sink system, and alluvial dynamics to explain these findings remains to be elucidated. Given the consistent presence of the 405 k.y. eccentricity cycle throughout Earth’s history, it is plausible to infer that its influence may be discernible across a wide array of alluvial stratigraphic records. ...
Conference paper (2024) - R. van der Kooij, A. Martinius, S. Geiger, G. Hampson, C. Jacquemyn, M. Jackson, D. Petrovskyy, K. Baird
In this contribution we compare results between fast screening methods based on flow diagnostics with results from full-physics simulations for CO2-brine displacement. Specifically, we analyse how either method identifies the key geological features that control CO2 plume migration across an ensemble of reservoir models for the shallow marine system present in the Johansen and Cook formations of the Northern Lights project. Our results indicate that many of the heterogeneities that appear to control CO2 plume migration based on the fast flow diagnostics also control CO2 plume migration when considering the complex fluid flow physics inherent to CO2-brine displacement. However, the magnitude of impact differs between both methods, conforming the hypothesis that more complex fluid flow processes increase the impact of geological heterogeneity on reservoir flow. Our results suggest that fast-screening methods like the one proposed by Jackson et al. (2022) are a reliable approach to analyse which heterogeneities at a given storage site are most likely to control CO2 plume migration, helping geoscientist and reservoir engineers to design more meaningful reservoir models that contain the key geological heterogeneities and allow us to predict CO2 plume migration more reliably. ...
Journal article (2024) - Domenico Chiarella, Sergio G. Longhitano, Luigi Spalluto, Allard W. Martinius
The hierarchies of the stratigraphic discontinuity surfaces observed in ancient tidalites are qualitatively assessed, aiming to evaluate their role as possible preferential conduits for fluid migration. Three outcrop examples are presented from microtidal settings of southern Italy: (i) siliciclastic tidalites consisting of quartz-rich cross-stratified sandstones generated by strong two-directional tidal currents flowing along a tidal strait; (ii) carbonate tidalites, which accumulated in a Cretaceous lagoon and tidal flat where peritidal cycles formed vertically-stacked sequences of biopeloidal and fenestral packstones, wackestones and bindstones during repeated phases of Milankovitch-scale sea-level changes; (iii) mixed, siliciclastic-bioclastic tidalites, deposited in a bay and recording offshore-transition, to shoreface wave-dominated and tide-influenced environments. Observations made during this study suggest that fluid movement can be controlled by the presence of main bounding surfaces that occur at different dimensions, from large (hectometre)-scale, to medium (decametre)-scale, to smaller (metre)-scales. These surfaces produced either by depositional or erosional processes, are characterised by different features and geometries in siliciclastic, carbonate and mixed siliciclastic-bioclastic tidalites arguably revealing complex internal pathways for fluid flows. These results suggest that fluids propagating along the main discontinuities follow a dominant sub-horizontal direction of propagation, associated with minor sub-vertical movements, due to local internal surface geometries and interconnections and a general lack of fractures. This surface-based approach to the study of fluid-flow transmission within stratified rocks represents a conceptual attempt to predict fluid mobility and reservoir potential in tidalite-bearing siliciclastic, carbonate and mixed reservoir rocks. ...
Journal article (2024) - Emilio Cecchetti, Allard W. Martinius, Pierre-Olivier Bruna, Annelies Bender, Hemmo A. Abels
The lower Triassic Main Buntsandstein Subgroup represents a promising, but high-risk geothermal play in the Netherlands. Although the gross thickness in boreholes locally exceeds 200 m, the spatial distribution, geometries and preservation of these sedimentary units remained uncertain due to the lack of seismic data with sufficient resolution and the sparse borehole network. This creates uncertainty in the quantification of the aquifer dimensions that is essential for the planning of geothermal operations.

In this study, seismic interpretation and 2D palinspastic restoration of new and reprocessed seismic data were conducted and combined with borehole data to assess the tectonic evolution of the Roer Valley Graben in the southeastern Netherlands and its control on the spatial distribution of the Main Buntsandstein Subgroup sediments. Our results show that the central and southern parts of the Roer Valley Graben were active depocenters in the Early to Middle Triassic times dominated by fluvial sandstone deposition, providing important play elements for prospective leads on geothermal exploration. The northern part of the basin was a more marginal area where mostly fine-grained sediments were deposited. To the northwest, differential subsidence resulted in the development of areas where the Buntsandstein thickness is reduced to ∼150 m.

After deposition, the Main Buntsandstein sediments were compartmentalised by faulting related to post-depositional tectonic activity, locally reducing the lateral extent of the geothermal target areas down to 1–2 km in a ∼NE–SW direction. On the platform areas adjacent to the Roer Valley Graben and to the southeast, Jurassic sediments are largely absent and the Main Buntsandstein sediments are present at depths shallower than 2 km. These platforms are promising targets for further investigation, as the relatively shallow burial depths, compared to the central part of the Graben, may have contributed to the preservation of more favourable reservoir properties. ...
Journal article (2024) - Jan H. Van Den Berg, Allard W. Martinius
In the past decades it has become evident that the often cyclic cross-shore migration of longshore bars is significantly influenced by wave climate. This study demonstrates that this cyclic migration, whether landward or seaward, leads to the formation of low-angle seaward-inclined stratification (SIS) spanning the beach and upper shoreface in ground-penetrating radar (GPR) data from prograding sandy coastal barriers including strandplains. Previously documented radargrams of these systems are reinterpreted using this knowledge of cross-shore dynamics of longshore bars. Five distinct wave-climate-related radar architectures are identified. A notable observation is the prevalence of SIS as the dominant upper-shoreface to beach structure in most radargrams, despite its infrequently described occurrence from outcrop observations. To address this disparity, this paper also focuses on recognizing SIS in outcrop and core data. Compelling evidence of SIS was discovered in a late Messinian coastal barrier parasequence in SE Spain, and the associated sedimentary architecture is described in detail. It is hypothesized that this example holds generic significance, illustrated using core data from Middle Jurassic strata of the Norwegian continental shelf, and suggesting that the sedimentary architectures of numerous other ancient prograding barrier deposits may exhibit variations within a similar SIS framework. ...
Abstract (2024) - Emilio Cecchetti, Allard W. Martinius, Marita Felder, Marinus E. Donselaar, Annelies Bender, Giovanni Bertotti, Pierre-Olivier Bruna, Hemmo A. Abels
Sandstones from the Main Buntsandstein Subgroup represent a promising deep geothermal target in the subsurface of the Netherlands considering their widespread distribution and temperatures locally reaching 140-150 °C at depths of ~ 3 to 5 km. The Main Buntsandstein Subgroup is a sand-prone interval, but the reservoir quality of these sandstones is known to be heterogeneous as result of an interplay between depositional and diagenetic processes. This makes the Buntsandstein sediments an uncertain and risky geothermal play.

In this project, we assess the syn- and post-depositional history of these sediments. The aim is to define structural, sedimentary, and diagenetic heterogeneities within the Main Buntsandstein sediments and assess their impact on reservoir quality. This will help reduce uncertainties for geothermal operations in the Triassic in the southern Netherlands and beyond.

The structural analysis of the study area using seismic and well data reveals that the Main Buntsandstein sediments represent an early syn-rift sequence and that their present-day distribution is strongly controlled by faulting. In parallel, the study of the sedimentology and stratigraphy conducted on core and wireline data indicates that the depositional environment evolves through the Buntsandstein stratigraphy, resulting in the development of different reservoir architectures. Diagenesis has largely altered the primary relationship between sedimentary facies and porosity and permeability. Overall cementation seems to have a larger impact on reducing reservoir quality than compaction, with quartz, dolomite, and illite representing the most abundant types of cement. The analysis of fractures using core and image logs suggests that the fracture density is driven by the lithological variability within the Main Buntsandstein and that fracture joints and stylolites locally may contribute to enhancing the system permeability.

The integrated assessment of the results allows the development of prospect play maps for the Buntsandstein in the southern Netherlands, addressing uncertainties and providing future recommendations for further exploration and optimizing geothermal operations in the Triassic.
...
Conference paper (2023) - G. Hampson, A. Martinius, M. Jackson, S. Krevor, J. Storms, D. Voskov, H. Hajibeygi, S. Geiger
This poster outlines a hierarchical, multiscale modelling approach that is adapted from proven hydrocarbon reservoir characterization workflows to determine which 3D sedimentological and stratigraphic heterogeneity types at which temporal and spatial scales and in which configurations are most important for successful long-term CO2 storage. The approach is particularly in saline aquifers that are data-poor but have the potential to store large CO2 volumes. It uses the Representative Element Volume (REV) concept and associated upscaling methodology to characterise sedimentological heterogeneity types, and it leverages novel modelling tools that facilitate rapid model construction and prototyping, geometrically accurate representation of key geological heterogeneities in models, and computationally efficient simulation of all CO2 trapping mechanisms over relevant spatial and temporal scales. ...
Many stratigraphic features occur at a scale that is at the edge or below vertical seismic resolution. Thus, they cannot be directly observed in the seismic data, while still having an important effect on the fluid flow within the system. The better understanding of these sub-seismic scale features or heterogeneities can help decrease subsurface uncertainty. Here we present a novel method that integrates forward stratigraphic modelling, petrophysics, and geophysics to decipher the seismic imprint of heterogeneities in wave-dominated, shallow marine environments. The proposed three-stepped method starts with defining geology-related input parameters for BarSim, a stratigraphic forward modelling software that produces models that include stratigraphic architecture, grain size distribution, and facies distribution. Then, the geological data is translated, cell by cell, into petrophysical data (density, Vp, and Vs) using emphirical relationships. Finally, the forward seismic modelling is performed by combining a finite difference approach strategy and angle-dependent full wavefield migration to retrieve the angle gathers This method also allows the generation of large amounts of field-independent data suitable for machine learning applications. ...
Journal article (2023) - Timothy F. Baars, Richard Huis in ‘t Veld, Linzhi Zhang, Maaike Koopmans, Duncan McLean, Allard W. Martinius, Hemmo A. Abels
Orbital driven climate control on sedimentation produces regional, stratigraphically repetitive characters and so cyclostratigraphic correlation can improve correlation and identify stratigraphic trends in borehole sections. This concept is commonly used to correlate marine and lacustrine strata. However, in the alluvial domain, its use is more challenging because internal, local dynamics controlling sedimentation may interfere with the expression of cyclic climate forcing. Intervals of low net-to-gross may be important for successful application in this domain as they tend to better document regional changes. This study applies climate-based stratigraphic correlation concepts to improve well correlations, characterise vertical sand distribution, and identify potential reservoir targets in a generally low net-to-gross interval. Coarsening upward sedimentary repetitions (cyclothems) are identified and correlated with high certainty in nineteen well sections in the upper Carboniferous Westoe and Cleaver formations of the Silverpit Basin. Local sedimentary dynamics provide variability in the character of the cyclothems and several types of cyclothem are classified. Correlation of sections using cyclothems recognised on wireline logs is done twice: once manually and once semi-automatically. The semi-automated correlation is based on calculation of deviation curves which depict stratigraphic changes that are less dependent on absolute wireline values and follow vertical trends more clearly. The correlations provide composite stratigraphies that are analysed using vertical proportions curves. Both approaches yield similar results in terms of stratigraphic trends. However, for detailed correlation of wells, the manual correlation is better at accounting for any local variability within the system. The same two zones of higher net-to-gross ratios are found using both correlation methods. These are linked to palaeoclimatic changes driven by long eccentricity and the proposed climate stratigraphic model has predictive value for identifying sandstone occurrence. The climate-based stratigraphic correlation improves the assessment reservoir distribution and properties on small (10–20 m thickness) and large (100–200 m thickness) stratigraphical scales. ...

Investigating the Role of Sedimentary Facies

Abstract (2023) - E. Cecchetti, A.W. Martinius, M Felder, M.E. Donselaar, H.A. Abels
The Lower Triassic Main Buntsandstein Subgroup represents one of the most promising deep geothermal plays in the Netherlands. Reservoir zones consist of sandstone units deposited in an arid to semi-arid alluvial plain reaching total thicknesses of over 200 m in the Roer Valley Graben. However, the lack of a comprehensive reservoir model integrating stratigraphy, sedimentology, and petrography makes the Buntsandstein a high-risk target. Here, we present a reservoir architecture model of the spatial and temporal variation of sedimentary facies and inherent permeability barriers and baffles, based on the integration of data from 33 wells from different stratigraphic levels and different areas of the Roer Valley Graben.

Core samples and thin sections analysis revealed two major reservoir facies that were interpreted as the products of transport and depositional processes in braided and sinuous river settings. The identified facies were then coupled to wireline logs to assess the spatial and temporal variation in sedimentary architecture. The lower part of the stratigraphy is dominated by braided river reservoir facies with a high degree of connectivity, where regional lacustrine-playa lake sediments represent the main potential permeability barriers. By contrast, the upper part of the stratigraphy is characterized by an increase in the proportion of sinuous river complexes. These latter yield a lower degree of connectivity with different types of baffles such as intercalated fine-grained overbank sediments, abandonment plugs, bar-draping fines, and cemented dolocrete scour fills. These are much more localized compared to the braided complexes-related barriers, making the prediction of the upper stratigraphy architecture uncertain.
...