Uncertainties in geomechanical models – exhaustive vs. feasible approach
Moritz Ziegler
CORRESPONDING AUTHOR
Helmholtz-Zentrum Potsdam, Deutsches GeoForschungsZentrum
GFZ, 14473 Potsdam, Germany
Oliver Heidbach
Helmholtz-Zentrum Potsdam, Deutsches GeoForschungsZentrum
GFZ, 14473 Potsdam, Germany
Institut für Angewandte Geowissenschaften, TU Berlin, 10587 Berlin, Germany
Related authors
Denise Degen, Moritz Ziegler, Oliver Heidbach, Andreas Henk, Karsten Reiter, and Florian Wellmann
EGUsphere, https://doi.org/10.5194/egusphere-2024-2932, https://doi.org/10.5194/egusphere-2024-2932, 2024
Short summary
Short summary
Obtaining reliable estimates of the subsurface state distributions is essential to determine the location of e.g. potential nuclear waste disposal sites. However, providing these is challenging since it requires solving the problem numerous times yielding high computational cost. To overcome this, we use a physics-based machine learning method to construct surrogate models. We demonstrate how it produces physics-preserving predictions, which differentiates it from purely data-driven approaches.
Moritz O. Ziegler, Robin Seithel, Thomas Niederhuber, Oliver Heidbach, Thomas Kohl, Birgit Müller, Mojtaba Rajabi, Karsten Reiter, and Luisa Röckel
Solid Earth, 15, 1047–1063, https://doi.org/10.5194/se-15-1047-2024, https://doi.org/10.5194/se-15-1047-2024, 2024
Short summary
Short summary
The rotation of the principal stress axes in a fault structure because of a rock stiffness contrast has been investigated for the impact of the ratio of principal stresses, the angle between principal stress axes and fault strike, and the ratio of the rock stiffness contrast. A generic 2D geomechanical model is employed for the systematic investigation of the parameter space.
Karsten Reiter, Oliver Heidbach, and Moritz O. Ziegler
Solid Earth, 15, 305–327, https://doi.org/10.5194/se-15-305-2024, https://doi.org/10.5194/se-15-305-2024, 2024
Short summary
Short summary
It is generally assumed that faults have an influence on the stress state of the Earth’s crust. It is questionable whether this influence is still present far away from a fault. Simple numerical models were used to investigate the extent of the influence of faults on the stress state. Several models with different fault representations were investigated. The stress fluctuations further away from the fault (> 1 km) are very small.
Moritz O. Ziegler, Oliver Heidbach, and Mojtaba Rajabi
Saf. Nucl. Waste Disposal, 2, 79–80, https://doi.org/10.5194/sand-2-79-2023, https://doi.org/10.5194/sand-2-79-2023, 2023
Short summary
Short summary
The subsurface is subject to constant stress. With increasing depth, more rock overlies an area, thereby increasing the stress. There is also constant stress from the sides. Knowledge of this stress is fundamental to build lasting and safe underground structures. Very few data on the stress state are available; thus, computer models are used to predict this parameter. We present a method to improve the quality of the computer models, even if no direct data on the stress state are available.
Karsten Reiter, Oliver Heidbach, Moritz Ziegler, Silvio Giger, Rodney Garrard, and Jean Desroches
Saf. Nucl. Waste Disposal, 2, 71–72, https://doi.org/10.5194/sand-2-71-2023, https://doi.org/10.5194/sand-2-71-2023, 2023
Short summary
Short summary
Numerical methods can be used to estimate the stress state in the Earth’s upper crust. Measured stress data are needed for model calibration. High-quality stress data are available for the calibration of models for possible radioactive waste repositories in Switzerland. A best-fit model predicts the stress state for each point within the model volume. In this study, variable rock properties are used to predict the potential stress variations due to inhomogeneous rock properties.
Luisa Röckel, Steffen Ahlers, Sophia Morawietz, Birgit Müller, Tobias Hergert, Karsten Reiter, Andreas Henk, Moritz Ziegler, Oliver Heidbach, and Frank Schilling
Saf. Nucl. Waste Disposal, 2, 73–73, https://doi.org/10.5194/sand-2-73-2023, https://doi.org/10.5194/sand-2-73-2023, 2023
Short summary
Short summary
Stress data predicted by a geomechanical–numerical model are mapped onto 3D fault geometries. Then the slip tendency of these faults is calculated as a measure of their reactivation potential. Characteristics of the faults and the state of stress are identified that lead to a high fault reactivation potential. An overall high reactivation potential is observed in the Upper Rhine Graben area, whereas the reactivation potential is quite low in the Molasse Basin.
Tobias Hergert, Steffen Ahlers, Luisa Röckel, Sophia Morawietz, Karsten Reiter, Moritz Ziegler, Birgit Müller, Oliver Heidbach, Frank Schilling, and Andreas Henk
Saf. Nucl. Waste Disposal, 2, 65–65, https://doi.org/10.5194/sand-2-65-2023, https://doi.org/10.5194/sand-2-65-2023, 2023
Short summary
Short summary
In numerical geomechanical models, an initial stress state is established before displacement boundary conditions are applied in order to match calibration data. We present generic models to show that the choice of initial stress and boundary conditions affects the final state of stress in areas of the model domain where no stress data for calibration are available. These deviations are largest in the vicinity of lithological interfaces, and they can be reduced if more stress data exist.
Steffen Ahlers, Karsten Reiter, Tobias Hergert, Andreas Henk, Luisa Röckel, Sophia Morawietz, Oliver Heidbach, Moritz Ziegler, and Birgit Müller
Saf. Nucl. Waste Disposal, 2, 59–59, https://doi.org/10.5194/sand-2-59-2023, https://doi.org/10.5194/sand-2-59-2023, 2023
Short summary
Short summary
The recent crustal stress state is a crucial parameter in the search for a high-level nuclear waste repository. We present results of a 3D geomechanical numerical model that improves the state of knowledge by providing a continuum-mechanics-based prediction of the recent crustal stress field in Germany. The model results can be used, for example, for the calculation of fracture potential, for slip tendency analyses or as boundary conditions for smaller local models.
Sophia Morawietz, Moritz Ziegler, Karsten Reiter, and the SpannEnD Project Team
Saf. Nucl. Waste Disposal, 1, 71–72, https://doi.org/10.5194/sand-1-71-2021, https://doi.org/10.5194/sand-1-71-2021, 2021
Short summary
Short summary
Knowledge of the crustal stress state is important for the assessment of subsurface stability. In particular, stress magnitudes are essential for the calibration of geomechanical models that estimate a continuous description of the 3-D stress field from pointwise and incomplete stress data. We present the first comprehensive and open-access stress magnitude database for Germany, consisting of 568 data records. We introduce a quality ranking scheme for stress magnitude data for the first time.
Moritz O. Ziegler, Oliver Heidbach, John Reinecker, Anna M. Przybycin, and Magdalena Scheck-Wenderoth
Solid Earth, 7, 1365–1382, https://doi.org/10.5194/se-7-1365-2016, https://doi.org/10.5194/se-7-1365-2016, 2016
Short summary
Short summary
Subsurface engineering relies on sparsely distributed data points of the stress state of the earth's crust. 3D geomechanical--numerical modelling is applied to estimate the stress state in the entire volume of a large area. We present a multi-stage approach of differently sized models which provide the stress state in an area of interest derived from few and widely scattered data records. Furthermore we demonstrate the changes in reliability of the model depending on different input parameters.
Denise Degen, Moritz Ziegler, Oliver Heidbach, Andreas Henk, Karsten Reiter, and Florian Wellmann
EGUsphere, https://doi.org/10.5194/egusphere-2024-2932, https://doi.org/10.5194/egusphere-2024-2932, 2024
Short summary
Short summary
Obtaining reliable estimates of the subsurface state distributions is essential to determine the location of e.g. potential nuclear waste disposal sites. However, providing these is challenging since it requires solving the problem numerous times yielding high computational cost. To overcome this, we use a physics-based machine learning method to construct surrogate models. We demonstrate how it produces physics-preserving predictions, which differentiates it from purely data-driven approaches.
Moritz O. Ziegler, Robin Seithel, Thomas Niederhuber, Oliver Heidbach, Thomas Kohl, Birgit Müller, Mojtaba Rajabi, Karsten Reiter, and Luisa Röckel
Solid Earth, 15, 1047–1063, https://doi.org/10.5194/se-15-1047-2024, https://doi.org/10.5194/se-15-1047-2024, 2024
Short summary
Short summary
The rotation of the principal stress axes in a fault structure because of a rock stiffness contrast has been investigated for the impact of the ratio of principal stresses, the angle between principal stress axes and fault strike, and the ratio of the rock stiffness contrast. A generic 2D geomechanical model is employed for the systematic investigation of the parameter space.
Karsten Reiter, Oliver Heidbach, and Moritz O. Ziegler
Solid Earth, 15, 305–327, https://doi.org/10.5194/se-15-305-2024, https://doi.org/10.5194/se-15-305-2024, 2024
Short summary
Short summary
It is generally assumed that faults have an influence on the stress state of the Earth’s crust. It is questionable whether this influence is still present far away from a fault. Simple numerical models were used to investigate the extent of the influence of faults on the stress state. Several models with different fault representations were investigated. The stress fluctuations further away from the fault (> 1 km) are very small.
Oliver Heidbach, Karsten Reiter, Moritz O. Ziegler, and Birgit Müller
Saf. Nucl. Waste Disposal, 2, 185–185, https://doi.org/10.5194/sand-2-185-2023, https://doi.org/10.5194/sand-2-185-2023, 2023
Short summary
Short summary
When stresses yield a critical value, rock breaks and generate pathways for fluid migration. Thus, the contemporary undisturbed stress state is a key parameter for assessing the stability of deep geological repositories. In this workshop you can ask everything you always wanted to know about stress (but were afraid to ask), and this is divided into three parts. 1) How do we formally describe the stress field? 2) How do we to actually measure stress? 3) How do we go from points to 3D description?
Moritz O. Ziegler, Oliver Heidbach, and Mojtaba Rajabi
Saf. Nucl. Waste Disposal, 2, 79–80, https://doi.org/10.5194/sand-2-79-2023, https://doi.org/10.5194/sand-2-79-2023, 2023
Short summary
Short summary
The subsurface is subject to constant stress. With increasing depth, more rock overlies an area, thereby increasing the stress. There is also constant stress from the sides. Knowledge of this stress is fundamental to build lasting and safe underground structures. Very few data on the stress state are available; thus, computer models are used to predict this parameter. We present a method to improve the quality of the computer models, even if no direct data on the stress state are available.
Karsten Reiter, Oliver Heidbach, Moritz Ziegler, Silvio Giger, Rodney Garrard, and Jean Desroches
Saf. Nucl. Waste Disposal, 2, 71–72, https://doi.org/10.5194/sand-2-71-2023, https://doi.org/10.5194/sand-2-71-2023, 2023
Short summary
Short summary
Numerical methods can be used to estimate the stress state in the Earth’s upper crust. Measured stress data are needed for model calibration. High-quality stress data are available for the calibration of models for possible radioactive waste repositories in Switzerland. A best-fit model predicts the stress state for each point within the model volume. In this study, variable rock properties are used to predict the potential stress variations due to inhomogeneous rock properties.
Luisa Röckel, Steffen Ahlers, Sophia Morawietz, Birgit Müller, Tobias Hergert, Karsten Reiter, Andreas Henk, Moritz Ziegler, Oliver Heidbach, and Frank Schilling
Saf. Nucl. Waste Disposal, 2, 73–73, https://doi.org/10.5194/sand-2-73-2023, https://doi.org/10.5194/sand-2-73-2023, 2023
Short summary
Short summary
Stress data predicted by a geomechanical–numerical model are mapped onto 3D fault geometries. Then the slip tendency of these faults is calculated as a measure of their reactivation potential. Characteristics of the faults and the state of stress are identified that lead to a high fault reactivation potential. An overall high reactivation potential is observed in the Upper Rhine Graben area, whereas the reactivation potential is quite low in the Molasse Basin.
Tobias Hergert, Steffen Ahlers, Luisa Röckel, Sophia Morawietz, Karsten Reiter, Moritz Ziegler, Birgit Müller, Oliver Heidbach, Frank Schilling, and Andreas Henk
Saf. Nucl. Waste Disposal, 2, 65–65, https://doi.org/10.5194/sand-2-65-2023, https://doi.org/10.5194/sand-2-65-2023, 2023
Short summary
Short summary
In numerical geomechanical models, an initial stress state is established before displacement boundary conditions are applied in order to match calibration data. We present generic models to show that the choice of initial stress and boundary conditions affects the final state of stress in areas of the model domain where no stress data for calibration are available. These deviations are largest in the vicinity of lithological interfaces, and they can be reduced if more stress data exist.
Steffen Ahlers, Karsten Reiter, Tobias Hergert, Andreas Henk, Luisa Röckel, Sophia Morawietz, Oliver Heidbach, Moritz Ziegler, and Birgit Müller
Saf. Nucl. Waste Disposal, 2, 59–59, https://doi.org/10.5194/sand-2-59-2023, https://doi.org/10.5194/sand-2-59-2023, 2023
Short summary
Short summary
The recent crustal stress state is a crucial parameter in the search for a high-level nuclear waste repository. We present results of a 3D geomechanical numerical model that improves the state of knowledge by providing a continuum-mechanics-based prediction of the recent crustal stress field in Germany. The model results can be used, for example, for the calculation of fracture potential, for slip tendency analyses or as boundary conditions for smaller local models.
Michal Kruszewski, Alessandro Verdecchia, Oliver Heidbach, Rebecca M. Harrington, and David Healy
EGUsphere, https://doi.org/10.5194/egusphere-2023-1889, https://doi.org/10.5194/egusphere-2023-1889, 2023
Preprint archived
Short summary
Short summary
In this study, we investigate the evolution of fault reactivation potential in the greater Ruhr region (Germany) in respect to a future utilization of deep geothermal resources. We use analytical and numerical approaches to understand the initial stress conditions on faults as well as their evolution in space and time during geothermal fluid production. Using results from our analyses, we can localize areas more favorable for geothermal energy use based on fault reactivation potential.
Michal Kruszewski, Gerd Klee, Thomas Niederhuber, and Oliver Heidbach
Earth Syst. Sci. Data, 14, 5367–5385, https://doi.org/10.5194/essd-14-5367-2022, https://doi.org/10.5194/essd-14-5367-2022, 2022
Short summary
Short summary
The authors assemble an in situ stress magnitude and orientation database based on 429 hydrofracturing tests that were carried out in six coal mines and two coal bed methane boreholes between 1986 and 1995 within the greater Ruhr region (Germany). Our study summarises the results of the extensive in situ stress test campaign and assigns quality to each data record using the established quality ranking schemes of the World Stress Map project.
Luisa Röckel, Steffen Ahlers, Birgit Müller, Karsten Reiter, Oliver Heidbach, Andreas Henk, Tobias Hergert, and Frank Schilling
Solid Earth, 13, 1087–1105, https://doi.org/10.5194/se-13-1087-2022, https://doi.org/10.5194/se-13-1087-2022, 2022
Short summary
Short summary
Reactivation of tectonic faults can lead to earthquakes and jeopardize underground operations. The reactivation potential is linked to fault properties and the tectonic stress field. We create 3D geometries for major faults in Germany and use stress data from a 3D geomechanical–numerical model to calculate their reactivation potential and compare it to seismic events. The reactivation potential in general is highest for NNE–SSW- and NW–SE-striking faults and strongly depends on the fault dip.
Luisa Röckel, Steffen Ahlers, Sophia Morawietz, Birgit Müller, Karsten Reiter, Oliver Heidbach, Andreas Henk, Tobias Hergert, and Frank Schilling
Saf. Nucl. Waste Disposal, 1, 77–78, https://doi.org/10.5194/sand-1-77-2021, https://doi.org/10.5194/sand-1-77-2021, 2021
Karsten Reiter, Steffen Ahlers, Sophia Morawietz, Luisa Röckel, Tobias Hergert, Andreas Henk, Birgit Müller, and Oliver Heidbach
Saf. Nucl. Waste Disposal, 1, 75–76, https://doi.org/10.5194/sand-1-75-2021, https://doi.org/10.5194/sand-1-75-2021, 2021
Steffen Ahlers, Andreas Henk, Tobias Hergert, Karsten Reiter, Birgit Müller, Luisa Röckel, Oliver Heidbach, Sophia Morawietz, Magdalena Scheck-Wenderoth, and Denis Anikiev
Saf. Nucl. Waste Disposal, 1, 163–164, https://doi.org/10.5194/sand-1-163-2021, https://doi.org/10.5194/sand-1-163-2021, 2021
Sophia Morawietz, Moritz Ziegler, Karsten Reiter, and the SpannEnD Project Team
Saf. Nucl. Waste Disposal, 1, 71–72, https://doi.org/10.5194/sand-1-71-2021, https://doi.org/10.5194/sand-1-71-2021, 2021
Short summary
Short summary
Knowledge of the crustal stress state is important for the assessment of subsurface stability. In particular, stress magnitudes are essential for the calibration of geomechanical models that estimate a continuous description of the 3-D stress field from pointwise and incomplete stress data. We present the first comprehensive and open-access stress magnitude database for Germany, consisting of 568 data records. We introduce a quality ranking scheme for stress magnitude data for the first time.
Steffen Ahlers, Andreas Henk, Tobias Hergert, Karsten Reiter, Birgit Müller, Luisa Röckel, Oliver Heidbach, Sophia Morawietz, Magdalena Scheck-Wenderoth, and Denis Anikiev
Solid Earth, 12, 1777–1799, https://doi.org/10.5194/se-12-1777-2021, https://doi.org/10.5194/se-12-1777-2021, 2021
Short summary
Short summary
Knowledge about the stress state in the upper crust is of great importance for many economic and scientific questions. However, our knowledge in Germany is limited since available datasets only provide pointwise, incomplete and heterogeneous information. We present the first 3D geomechanical model that provides a continuous description of the contemporary crustal stress state for Germany. The model is calibrated by the orientation of the maximum horizontal stress and stress magnitudes.
Ershad Gholamrezaie, Magdalena Scheck-Wenderoth, Judith Bott, Oliver Heidbach, and Manfred R. Strecker
Solid Earth, 10, 785–807, https://doi.org/10.5194/se-10-785-2019, https://doi.org/10.5194/se-10-785-2019, 2019
Short summary
Short summary
Based on geophysical data integration and 3-D gravity modeling, we show that significant density heterogeneities are expressed as two large high-density bodies in the crust below the Sea of Marmara. The location of these bodies correlates spatially with the bends of the main Marmara fault, indicating that rheological contrasts in the crust may influence the fault kinematics. Our findings may have implications for seismic hazard and risk assessments in the Marmara region.
Moritz O. Ziegler, Oliver Heidbach, John Reinecker, Anna M. Przybycin, and Magdalena Scheck-Wenderoth
Solid Earth, 7, 1365–1382, https://doi.org/10.5194/se-7-1365-2016, https://doi.org/10.5194/se-7-1365-2016, 2016
Short summary
Short summary
Subsurface engineering relies on sparsely distributed data points of the stress state of the earth's crust. 3D geomechanical--numerical modelling is applied to estimate the stress state in the entire volume of a large area. We present a multi-stage approach of differently sized models which provide the stress state in an area of interest derived from few and widely scattered data records. Furthermore we demonstrate the changes in reliability of the model depending on different input parameters.
T. Hergert, O. Heidbach, K. Reiter, S. B. Giger, and P. Marschall
Solid Earth, 6, 533–552, https://doi.org/10.5194/se-6-533-2015, https://doi.org/10.5194/se-6-533-2015, 2015
Short summary
Short summary
A numerical model integrating the structure and mechanical properties of a sedimentary sequence in the Alpine foreland is presented to show that topography, tectonic faults and, most of all, spatialy variable rock properties affect the state of stress at depth. The tectonic forces acting on the sequence are primarily taken up by the stiff rock units leaving the weaker units in a stress shadow.
K. Reiter and O. Heidbach
Solid Earth, 5, 1123–1149, https://doi.org/10.5194/se-5-1123-2014, https://doi.org/10.5194/se-5-1123-2014, 2014
Short summary
The Earth's crust is subject to constant stress which is manifested by earthquakes at plate boundaries. This stress is not only at plate boundaries but everywhere in the crust. A profound knowledge of the magnitude and orientation of the stress is important to select and build a safe deep geological repository for nuclear waste. We demonstrate how to build computer models of the stress state and show how to deal with the associated uncertainties.
The Earth's crust is subject to constant stress which is manifested by earthquakes at plate...