the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Updating induced seismic hazard assessments during hydraulic stimulation experiments in underground laboratories: workflow and limitations
Abstract. Advancing technologies to harvest deep geothermal energy has seen backlashes related to unacceptable levels of induced seismic hazard during hydraulic stimulations. A thorough analysis of induced seismic hazard before these operations has recently become standard practice in the last decade. Additionally, more process understanding of the underlying causes of induced seismicity as well as novel approaches to develop geomechanical reservoirs are being explored in controlled underground laboratory experiments world-wide. Here, we present a probabilistic analysis of the seismic hazard induced by the ongoing hectometer scale stimulation experiments at the Bedretto Underground Laboratory for Geoenergies and Geosciences (BULGG). Our workflow allows for fast updates of the hazard computation as soon as new site-specific information on the seismogenic response (expressed primarily by the feedback afb-value and the Gutenberg Richter b-value) and ground motion models (GMM) become available. We present a sequence of hazard analyses corresponding to different project stages at the BULGG. These reveal the large uncertainty in a priori hazard estimations that only reduce once site-specific GMMs and information on the seismic response of specific stimulation stages are considered. The sources of uncertainty are 1) the large variability in the seismogenic response recorded across all stimulation case studies, as well as 2) uncertain GMMs on the underground laboratory scale. One implication for large-scale hydraulic stimulations is that hazard computation must be updated at different project stages. Additionally, stimulations have to be closely accompanied by a mitigation scheme, ideally in the form of an adaptive traffic light system (ATLS), which reassesses seismic hazard in near-real-time. Our study also shows that the observed seismogenic responses in underground laboratories differ from large-scale stimulations at greater depth in that the seismogenic response is substantially more variable and tends to be weaker. Reasons may be lower stress levels, but also smaller injected volumes accessing a more limited fracture network than large-scale stimulations. Exploring the physical reasons leading to the weaker seismogenic response may reveal ways for safer exploitation of geoenergy resources. Controlled underground laboratory experiments can readily contribute to this, and – as shown in the presented analysis – are likely to be safe in terms of induced seismic hazard.
- Preprint
(5323 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 09 Apr 2025)
-
RC1: 'Comment on egusphere-2024-3882', P. Martin Mai, 15 Feb 2025
reply
REVIEW FOR egusphere-2024-3882
Updating induced seismic hazard assessments during hydraulic stimulation experiments in underground laboratories: workflow and limitations
by Gischig et al
+++++++++++++++++++++++++++++
General Comments
—————————————
Induced seismicity in context of enhanced oil & gas exploitation, wastewater injection, and geothermal-energy harvesting is a recurring problem that operators, regulators and nearby communities have to deal with. In case of geothermal energy, hydraulic simulations have led to induced seismicity at a level that caused substantial shaking such that operations were stopped because the associated seismic hazard was not tolerable any more. The question then arises if this (time-dependent) seismic hazard can be quantified and “controlled” during hydraulic stimulation.
The study by Gischig and colleagues examines this question, with a focus on hectometer-scale stimulations in the Bedretto Underground Laboratory for Geoenergies and Geoscience (BULAGG) in Switzerland. Inspired from observations and lessons learned in several geothermal projects around the globe, the authors develop a work-flow to compute/update the seismic hazard at a location of interest given fluid-injection parameters and the overall boundary and initial conditions at the site in terms of geology, seismotectonics & regional stresses, whereby the seismic-hazard updates are based on the known injection history and measured seismicity parameters. Noting that hazard estimates may vary greatly depending on the state of information/data and can be better constrained with more data and refined seismicity parameters, the authors also stress that site-specific ground-motion data and a related ground-motion model (GMM) are critical to narrow down the hazard estimates to plausible ranges. The study concludes with proposing an adaptive traffic light system (aTLS) that capture the time-dependent seismic hazard changes in near-real time.
The manuscript is well written, with very accessible graphics and a well-composed structure that naturally navigates the reader through the rather comprehensive material in terms of previous studies, the site of interest, related observations, models developed in the past and applied for the chosen case study, hazard calculations and how these are eventually embedded into a traffic light system. That is, the paper is rich. It is dense. It contains a lot of information that the reader needs to digest. In my view, the authors did an excellent job in this regard, but, I do remark that for most 1st to 2nd year graduate students in this field and also the “general but interested reader”, this paper may not be an easy read.
From a technical point of view, I don’t have any major comments and concerns. The science is solid. The methods are well known (but not all explained in detail, hence readers need sufficient background knowledge), the data are exquisite, and the overall goal of the study is of importance scientifically and from a socio-economic point of view. Nevertheless, I have several remarks and questions related to the presentation, level of detail provided on certain aspects of the study (sometimes too much, and thus distracting from the “big picture”’; sometimes too sparse to be able to follow), and a few editorial remarks.
The one major point I would like to raise is the use of peak ground velocity (PGV) instead of peak ground acceleration (PGA) as ground-motion intensity measure. The earthquakes consider here are predominantly of small magnitudes and clearly dominated by high-frequency seismic radiation. The authors also state the most earthquakes studied radiate above 10 Hz. A rule-of-thumb is that PGV captures shaking intensities for waves around 1 Hz. Hence, the use of PGV is counter-intuitive and perhaps not physically justifiable. This aspects needs detailed consideration and explanation (see below for more on this).
Overall, I rate this paper as “publishable after minor revisions”. New calculations/analyses or substantial re-organization/rewriting are not needed, but I ask for clarifications and editorial improvements that should be straightforward to implement.
In summary, this is a very interesting and well-written manuscript that I think will be quite impactful.
Below, I provide a few technical comments follow by minor editorial suggestions.
Technical Comments
—————————-
+ Figure 2: It would be help to graphically show here the principal stress orientations (and magnitudes) discussed in Lines 139ff. That would help the readers to quickly grasp all tectonic details.
+ Figure 3: For the laymen readers, these 3D plots are hard-to-impossible to put in context. In essence, a detailed map / 3D graphic is needed that shows the locations of these boreholes within the Bedretto Lab. I am not sure if these locations can be easily added to Figure 2, or if another zoomed-in close up near the underground lab is needed. Please consider.
+ Section 3 (Instrumentation and Experiments) can be and perhaps should be deleted. In my opinion, these details are not needed to understand this study in terms of science, methods, scope and final results. On the other hand, this section distracts from the main “story line” and the main goals of this paper. If deemed important for this completeness purposes, I suggest to move this section into an Esupplement/Appendix.
+ Figure 4: Which magnitude scale is used here, Ml or Mw? Please indicate. In general, since this topic comes up later again, I suggest to explain already early on how magnitudes are estimated, if Ml or Mw is routinely/automatically determined, and with which uncertainties.
+ Line 275: The “simplifying assumption that the b-value remains constant during injection and after shut-in” is an interesting point to (re-)consider. First of all, is that assumption valid? Given the wealth of data and the experience of the team of author, this should be a very quick and easy point to check and verify. My suspicion is that this is not the case, looking at Figure 4. Perhaps time-dependent b-values, and the variations, over different time-window lengths could be computed to check if/when this assumption is correct. And if not, then we need to think about how this may be propagated into the later hazard calculation.
+ Figure 7: Panels a) and b) need some modifications. First, the y-axis range in both panels should be identical. Second, the grey-scale density plot in panel b) is too fuzzy and doesn’t allow being able to see details. I suggest to use a distinct colorbar with 6-10 visually clearly separable colors (say at 0.1, 0.2, 0.3 …) so that details can be seen.
+ Lines 326 - 335: Here, reference should be made to Galis et al (2017) (already in the reference list) and perhaps to Gabriel et al (2024, in Science) and Palgunadi et al (2024, JGR) on arrested and run-away ruptures in complex-geometry fault systems.
+ Lines 347 - 352: Using a simple constant stress-drop assumption to translate an estimated fault dimension to a possible event magnitude seems too simplistic, too approximate, and does not include any uncertainty. I strongly recommend to apply modern source-scaling relations (i.e. Thingbaijam et al, 2017), possibly also considering different faulting styles, to estimate potential event magnitude and its range.
+ Figure 8: I strongly recommend to plot the scaling relations van der Elst et al (2016) and Galis et al (2017) into this figure for completeness and reference. (This will also shorten the figure caption by two lines …).
+ Section Ground motion models: As someone who has experience in PSHA and GMM’s for “standard” regional/national seismic hazard assessment, I am puzzled that PGV is used as ground-motion intensity metric, instead of PGA. I realize this may be the engineering/operational practice in mining-seismicity studies, but this is very counter-intuitive, in particular because we are dealing with very small events that are dominated by high-frequency radiation, and hence PGV may not be an ideal shaking parameter to use. I suggest that the authors provide some clarification and rationale for their choice,.
+ Figure 9: The Cai-Kaiser (2018) model seems to be an almost exact replicate of McGarr & Fletcher (2005), just shifted downwards by “-1 log bias unit”. Is that the case? Perhaps an explanatory sentence.
+ Line 411: The wording “site-specific information” confuses me here, since it is not clear what the “site” is. In BULGG, there are numerous seismic sensors that each could be considered a “recording site”. On the other hand, the overall spatial foot-print of BULGG or any similar experimental facility is rather small and would be typically considered as a “single site” in any local/regional PSHA study. Please clarify.
+ Line 420: “induced earthquake … have frequencies higher than 10 Hz” —> this relates back to my comment above: Why is then PGV a useful ground-motion metric? And wouldn’t PGA make much more sense?
+ Figure 14: For the 10 panels shown on top, I suggest not to use a continuous colorbar-scale, but one with 8-12 clearly distinct color. Visually, the hues of red between, say 200 - 800 mm/s cannot be discriminated.
+ Line 537-539: The fact that PSHA estimate increase as more data are added is in fact a widely occurring but not well appreciated fact, in general; not only in the context of induced seismicity. I suggest to add corresponding references from the PSHA literature.
+ Sub-section Scale dependent seismogenic response: I would have expected at least a short discussion on whether there are dependencies of the b-value on the faulting-style of the earthquake. Schorlemmer at al (2005?) found a very compelling dependence of the b-value given the faulting-style, which in turn can be explained by the dominant acting stress regime. I suggest to add a few sentences on this here.
+ Line 606: somewhere close to the reference to Deichmann (2017) and in this section there should also be made reference to two papers by Bethmann et al (BSSA, 2011, and GJI, 2012) that examine Mw-Ml scaling relations and site/attenuation effects on Ml/Mw estimates in Switzerland.
Editorial Suggestions
——————————
+ The authors refer to “hazard” and “risk” numerous times in the paper, and I do understand that they want to clearly distinguish the two. However, in several instances this distinction is not clear and then things become confusing. Because there are no risk calculations included here and risk is only referred to in a general sense, I suggest the authors add a specific “item at risk” in corresponding statements, for example “risk for tunnel collapse”, or “risk to geothermal surface facilities” to better guide the readers what they in mind in each case.
+ Please carefully check the punctuation. I noticed many missing periods (“ .”) to conclude sentences, but even more so I found incorrect setting of commas (“ ,”) that lead to confusion in terms of meaning of the respective sentences.
Other points:
+ Line 100: move “Sweden” after Aspo (in Line 99)
+ Line 106: “intense” is not a good word here, as it cannot be quantified. Use something more specific: real-time; high-resolution (in space, time and frequency frequency) or something like that …
+ Line 109: “seismic risk” … see above …
+ Line 281: (and others) - the text refers to Table 1, but this does not exist; it is Table A1 in the Appendix. This may just be a formatting issue or problem with the latex-template, but please check such referencing carefully.
+ Line 314: The wording “moderate” seems unclear here. Moderate “magnitude”? But what magnitude would that be in the context of the event sizes shown here? Or perhaps better “more frequent events”?
+ Line 390: The “Table” mentioned here is given in the Appendix. Please correct.
+ Line 463: Figure caption to Figure 11: The “hazard” here should be clearly specified as “Hazard to exceed a certain earthquake magnitude”. Most readers associate “hazard” with “seismic” (i.e. “shaking hazard”) …
+ Line 487: The title to this sub-section should be “seismic” or ‘shaking’ hazard …
+ Line 505: See comment above to Line 463 / 483
+ Lines 540 - 544: This sentence seems garbled up; I cannot understand it. Also, change “became” to “become”.
+ Line 548: remove or quantify “somewhat”
+ Line 596: abbreviation “GSK” not defined
Final Note
——————
I did not check in detail the completeness and correctness of the references, nor the exact entries of the two tables in the Appendix.
Citation: https://doi.org/10.5194/egusphere-2024-3882-RC1
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
131 | 28 | 4 | 163 | 3 | 8 |
- HTML: 131
- PDF: 28
- XML: 4
- Total: 163
- BibTeX: 3
- EndNote: 8
Viewed (geographical distribution)
Country | # | Views | % |
---|---|---|---|
United States of America | 1 | 72 | 34 |
Switzerland | 2 | 34 | 16 |
Germany | 3 | 20 | 9 |
China | 4 | 10 | 4 |
France | 5 | 9 | 4 |
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
- 72