the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Evolution of fault reactivation potential in deep geothermal systems. Insights from the greater Ruhr region, Germany
Michal Kruszewski
Alessandro Verdecchia
Oliver Heidbach
Rebecca M. Harrington
David Healy
Abstract. The success of deep geothermal systems depends on the presence of fault zones in the subsurface. Faults play a vital role in the Earth's plumbing system by facilitating fluid flow when they dilate, but are simultaneously known to enhance the hazard of the system once slipping in shear mode. As dilation of a fault enhances its permeability significantly, shear failure can lead to loss of boreholes or seismic events of economic concern. In this study, we present the evolution of reactivation potential of major faults during 25-year production period in deep generic geothermal systems in the greater Ruhr region in western Germany. To determine the pre-operational in situ stress state we use a recently published comprehensive dataset of stress magnitude data from the greater Ruhr region in an analytical-probabilistic model accounting for uncertainties of in situ stress, fault geometry, and frictional properties for a prospective reservoir in the Devonian Massenkalk formations. The resulting cumulative distribution functions of dilation and slip tendency of given fault sets suggests that more than half of the combined length of NW-SE-striking faults have a high reactivation probability, whereas the NE-SW-striking faults remain not optimally-oriented in the regional stress field. Using the relationship between dilation and slip tendency, we propose fault segments suitable for geothermal development that exhibit high hydraulic conductivity, i.e. high dilation tendency, and lower potential for shear failure, i.e. low slip tendency. In the second step, we employ generic thermo-hydro-mechanical models to quantify induced spatio-temporal stress changes on selected fault planes due to long-term geothermal production. We find that after 25 years thermal stress changes contribute significantly to the change of the reactivation potential which should be accounted for while planning deep geothermal systems.
- Preprint
(34542 KB) - Metadata XML
- BibTeX
- EndNote
Michal Kruszewski et al.
Status: open (until 10 Oct 2023)
-
RC1: 'Comment on egusphere-2023-1889', Anonymous Referee #1, 27 Sep 2023
reply
General comments
This manuscript intends to discuss the fault reactivation potential in deep geothermal systems that would be developed in the greater Ruhr region, Germany. To do so, two main approaches are considered. The first one is based on the well-known static slip and dilation tendencies complemented by the critical pore pressure. The second one is based on a thermo-hydro-mechanical numerical modelling of a synthetic geothermal doublet.
This manuscript presents approaches and results very similar to previous studies without relevant increase of information. The concepts and ideas are already well established. Along the manuscript, necessary details and quality controls are missing, especially regarding the numerical reservoir modelling, which prevents from being confident in the quantified results. Moreover, the deterministic THM modelling remains very generic and therefore delivers common knowledge and general results. The scientific significance is missing over the manuscript. Besides, the manuscript is wordy and extrapolates the results. Consequently, I recommend that this manuscript be rejected for publication in Solid Earth and any other journal.
Specific comments
The two main criticisms are associated with both approaches developed in the manuscript.
Slip tendency, dilation tendency and critical pore pressure
1. The methodology and the results presented have been discussed in a very similar way and for the same area by the same main author in:
- Kruszewski, M., G. Montegrossi, T. Backers, and E. H. Saenger, 2021, In Situ Stress State of the Ruhr Region (Germany) and Its Implications for Permeability Anisotropy: Rock Mechanics and Rock Engineering, 54, 6649–6663.
- Kruszewski, M., G. Klee, T. Niederhuber, and O. Heidbach, 2022a, In situ stress database of the greater Ruhr region (Germany) derived from hydrofracturing tests and borehole logs: Earth System Science Data, 14, 5367–5385.
- Kruszewski, M., G. Montegrossi, M. Balcewicz, G. de Los Angeles Gonzalez de Lucio, O. A. Igbokwe, T. Backers, and E. H. Saenger, 2022b, 3D in situ stress state modelling and fault reactivation risk exemplified in the Ruhr region (Germany): Geomechanics for Energy and the Environment, 32, 100386.
Despite the probabilistic approach was applied in this manuscript, the added-value is questionable.
2. Regarding the probabilistic approach, why the procedure presented by Seithel et al. (2019) (a paper your refer to) is not used although it was developed for the same purpose?
3. The slip tendency results (Fig. 4) highlight fault patches that have Ts higher than 1. How can it be?
4. With the presented results, many faults should be already critically stressed. How do you explain that no natural seismicity is observed in the area? A chapter discussing the natural seismicity of the area is missing.
5. For the range of pressure found in the Ruhr area at that depth (<2 kbar), Byerlee (1978) observed friction coefficients of 0.85. However, the limit assumed in the manuscript is 0.6 (L132, L187, Table 2) and not 0.85, using the same reference, why? L210, however, refereeing to Byerlee (1978) again, the 0.85 friction is written to be a possible value!
6. For Sf, in Eq. 3, Co is accounted for, but it does not appear for the slip tendency although both parameters (Ts and Sf) have the same theoretical background (Mohr-Coulomb failure criterion). Why is it so?
7. In section 5.2, L378-379: “Scalar values used for fault stability evaluation based on the contribution of fluid pressure only, such as Sf, will not provide a full picture of the fault stability in situ”. This is also true for the slip tendency, so mention it as well.
THM modeling
8. The last sentence of the conclusion is not surprising and does not need any result of the numerical simulation that was described in the manuscript. Below is a (non-exhaustive) list of papers that are already discussing the importance of thermally induced stress changes on a long term basis in geothermal contexts:
- De Simone, S., V. Vilarrasa, J. Carrera, A. Alcolea, and P. Meier, 2013, Thermal coupling may control mechanical stability of geothermal reservoirs during cold water injection: Physics and Chemistry of the Earth, Parts A/B/C, 64, 117–126.
- Egert, R., Gaucher, E., Savvatis, A., Goblirsch, P., Kohl, T., 2022. Numerical determination of long-term alterations of THM characteristics of a Malm geothermal reservoir during continuous exploitation. Presented at the European Geothermal Congress 2022, Berlin, Germany.
- Jeanne, P., J. Rutqvist, and P. F. Dobson, 2017, Influence of injection-induced cooling on deviatoric stress and shear reactivation of preexisting fractures in Enhanced Geothermal Systems: Geothermics, 70, 367–375.
- Jeanne, P., J. Rutqvist, P. F. Dobson, J. Garcia, M. Walters, C. Hartline, and A. Borgia, 2015, Geomechanical simulation of the stress tensor rotation caused by injection of cold water in a deep geothermal reservoir: Journal of Geophysical Research: Solid Earth, 120, 8422–8438.
- Kivi, I. R., E. Pujades, J. Rutqvist, and V. Vilarrasa, 2022, Cooling-induced reactivation of distant faults during long-term geothermal energy production in hot sedimentary aquifers: Scientific Reports, 12, 2065.
- Koh, J., H. Roshan, and S. S. Rahman, 2011, A numerical study on the long term thermo-poroelastic effects of cold water injection into naturally fractured geothermal reservoirs: Computers and Geotechnics, 38, 669–682.
- Wassing, B. B. T., T. Candela, S. Osinga, E. Peters, L. Buijze, P. A. Fokker, and J. D. Van Wees, 2021, Time-dependent Seismic Footprint of Thermal Loading for Geothermal Activities in Fractured Carbonate Reservoirs: Frontiers in Earth Science, 9.
Many references regarding THM modelling in similar contexts should be given but they are missing. They could have been used as inspiration source.
9. When presenting THM numerical modelling, it is necessary to develop much more what is actually done to give confidence in the results. So far, it is not the case, and a lot of information is missing, e.g. what are the physical processes activated (equations)? The above-mentioned papers could help to do so.
10. Was a mesh sensitivity study carried out? This is questionable when looking at the discontinuous curves of Fig. 7a and b.
11. In the THM results, it would be most important to see space and time distribution of, at least, the pore-pressure field and the temperature field before jumping directly to the shortest distance between wells and fault.
12. Section 3.2, L165-168: “Effects such as fault permeability enhancement due to the dilation, change of rock properties due to Pp or temperature, T, the influence of fluid chemistry on rock mass and fault properties, mechanisms of earthquake interactions, and the Kaiser effect are not considered in the simulation” This looks like COMSOL could account for all of these. I am not aware that COMSOL can simulate earthquakes.
Additional comments
13. First sentence of abstract: This is wrong as underlined e.g. by the deep geothermal exploitation in the Paris basin for many decades.
14. Nothing in the manuscript supports the simulation of seismicity or aseismic slip or seismic hazard. Consequently, these aspects should be mentioned with care.
Technical corrections
15. Second sentence of abstract: what is the Earth’s “plumbing” system? I have never read such wording in a geothermal context. Do you mean “circulation”?
16. L18: […] a complex “web” of faults […]? I have never read such wording in a geothermal context. Do you mean “network”? It is found L401 as well.
17. Avoid using “the distance to failure” (e.g. L31), you mean in meters (?), prefer the “reactivation potential”.
18. The Appendix does not correspond at all to what is announced in the main part of the manuscript.
There would be a lot more to say, but that was the gist.
Citation: https://doi.org/10.5194/egusphere-2023-1889-RC1
Michal Kruszewski et al.
Michal Kruszewski et al.
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
301 | 123 | 11 | 435 | 6 | 6 |
- HTML: 301
- PDF: 123
- XML: 11
- Total: 435
- BibTeX: 6
- EndNote: 6
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1