Modeling internal deformation of salt structures targeted for radioactive waste disposal

Abstract: This thesis uses results of systematic numerical models to argue that externally inactive salt structures, which are potential targets for radioactive waste disposal, might be internally active due to the presence of dense layers or blocks within a salt layer.The three papers that support this thesis use the Gorleben salt diapir (NW Germany), which was targeted as a future final repository for high-grade radioactive waste, as a general guideline.The first two papers present systematic studies of the parameters that control the development of a salt diapir and how it entrains a dense anhydrite layer. Results from these numerical models show that the entrainment of a dense anhydrite layer within a salt diapir depends on four parameters: sedimentation rate, viscosity of salt, perturbation width and the stratigraphic location of the dense layer. The combined effect of these four parameters, which has a direct impact on the rate of salt supply (volume/area of the salt that is supplied to the diapir with time), shape a diapir and the mode of entrainment. Salt diapirs down-built with sedimentary units of high viscosity can potentially grow with an embedded anhydrite layer and deplete their source layer (salt supply ceases). However, when salt supply decreases dramatically or ceases entirely, the entrained anhydrite layer/segments start to sink within the diapir. In inactive diapirs, sinking of the entrained anhydrite layer is inevitable and strongly depends on the rheology of the salt, which is in direct contact with the anhydrite layer. During the post-depositional stage, if the effective viscosity of salt falls below the threshold value of around 1018-1019 Pa s, the mobility of anhydrite blocks might influence any repository within the diapir. However, the internal deformation of the salt diapir by the descending blocks decreases with increase in effective viscosity of salt.The results presented in this thesis suggest that it is highly likely that salt structures where dense and viscous layer/blocks are present undergo an internal deformation processes when these dense blocks start sinking within the diapir. Depending on size and orientation of these blocks, deformation pattern is significantly different within the diapir. Furthermore, model results applied to the Gorleben diapir show that the rate of descent of the entrained anhydrite blocks differs on different sides of the diapir. This suggests that if the anhydrite blocks descent within the Gorleben diapir, they initiate an asymmetric internal flow within it.