ThermoQuakes

ThermoQuakes – 3D thermo-mechanical modelling of seismicity related to the final disposal of spent nuclear fuel and radioactive waste in hard rock

Duration: 2014.11-2017.11

Funding: Swedish Radiation Safety Authority (SSM)

PI: Jeoung Seok Yoon, Arno Zang

Reference Group: Ove Stephansson (KTH, GFZ), Hilmar Bungum (NORSAR), Fabrice Cotton (GFZ)

To be able to assess the safety of a repository for radioactive waste and spent nuclear fuel, it is necessary to consider all possible threats that could impair the physical integrity of the barriers. One environmental threat is an earthquake occurring at nearby faults, and the effect of an earthquake relevant to the repository safety is the shear displacement of rock fractures induced by an earthquake event occurring at nearby large faults (Fig.1).

Fig.1. (a) Layout of the repository and the bounding faults (modified from SKB ; (b) Shear deformation of a rock fracture cutting across a deposition hole and damaging a canister containing spent nuclear fuels (Fälth et al. 2010).

Fig.1. (a) Layout of the repository and the bounding faults (modified from SKB ; (b) Shear deformation of a rock fracture cutting across a deposition hole and damaging a canister containing spent nuclear fuels (Fälth et al. 2010).


In this project, we use 3D discrete element based thermo-mechanical coupled dynamic modelling to investigate the safety of the repository at Forsmark Sweden for final disposal of spent nuclear fuels. Especially, we investigate the impacts of larger magnitude (M>6) earthquakes that could occur in the next few hundred years near the Forsmark site. Followings are the objectives of the project.

  • Generation of bonded-particle based 3D geological model of Forsmark repository site with heterogeneously represented geological discontinuities, e.g. faults and fractures
  • Development of a numerical workflow for simulation of fault dynamic rupture (Yoon et al. 2014, 2017)
  • Investigation of impacts of radioactive decaying heat to the repository fracture system (Yoon et al. 2016)
  • Investigation of impacts of earthquake events to the repository fracture system (Fig.2)
  • Verification and validation of the simulated earthquake against earthquake fault scaling relation (Fig.3)
Fig.2. Temporal change of the slip (in meter) distribution of a dynamically rupturing fault. The co-seismic slip concentrates at one point on the fault plane and propagates outward, which results in asymmetric and heterogeneous slip profiles. Black dots are the deposition holes.

 

Fig.2. Temporal change of the slip (in meter) distribution of a dynamically rupturing fault. The co-seismic slip concentrates at one point on the fault plane and propagates outward, which results in asymmetric and heterogeneous slip profiles. Black dots are the deposition holes.
(Click for animation)

Fig.3. Rupture area vs. moment magnitude relation of earthquake faults. The results of simulated earthquake faulting (black stars) are compared with the global tectonic earthquake fault data (dots).

 

Fig.3. Rupture area vs. moment magnitude relation of earthquake faults. The results of simulated earthquake faulting (black stars) are compared with the global tectonic earthquake fault data (dots).

Key references

  • Yoon JS, Stephansson O, Zang A, Min KB, Lanaro F (2017), Discrete bonded particle modelling of fault activation near a nuclear waste repository site and comparison to static rupture earthquake scaling laws, Int J Rock Mech Min Sci, 98, 1-9, http://dx.doi.org/10.1016/j.ijrmms.2017.07.008
  • Yoon JS, Stephansson O, Zang A, Min KB, Lanaro F (2016), Numerical modelling of earthquakes and induced seismicity under various in situ stress conditions at Forsmark, Sweden, the site for a final repository of spent nuclear fuel, RS2016 Symposium – 7th International Symposium on In-Situ Rock Stress, May 10-12, 2016, Tampere, Finland.

Kontakt

Profilfoto von  Dr.-Ing. Jeoung Seok Yoon

Dr.-Ing. Jeoung Seok Yoon
Erdbebengefährdung und Spannungsfeld

Helmholtzstraße 6/7
Haus H 6, Raum 207
14467 Potsdam
Tel. +49 331 288-1716

Kontakt

Profilfoto von  apl. Prof. Arno Zang

apl. Prof. Arno Zang
Erdbebengefährdung und Spannungsfeld

Helmholtzstraße 6/7
Haus H 6, Raum 210
14467 Potsdam
Tel. +49 331 288-1325