The presence of salt in dilatant normal faults may have a strong influence on fault mechanics in the Groningen field and on the related induced seismicity. At present, little is known of the structure of these fault zones. This study starts with the geological evolution of the Groningen area, where, during tectonic faulting, rock salt may have migrated downwards into dilatant faults. These fault zones therefore may contain inclusions of rock salt. Because of its rate-dependent mechanical properties, the presence of salt in a fault may introduce a loading-rate dependency into fault movement and affect the distribution of magnitudes of seismic events. We present a first-look study showing how these processes can be investigated using a combination of analogue and numerical modelling. Full scaling of the models and quantification of implications for induced seismicity in Groningen require further, more detailed studies: an understanding of fault zone structure in the Groningen field is required for improved predictions of induced seismicity. The analogue experiments are based on a simplified stratigraphy of the Groningen area, where it is generally thought that most of the Rotliegend faulting has taken place in the Jurassic, after deposition of the Zechstein. This suggests that, at the time of faulting, the sulphates were already transformed into brittle anhydrite. If these layers were sufficiently brittle to fault in a dilatant fashion, rock salt was able to flow downwards into the dilatant fractures. To test this hypothesis, we use sandbox experiments where we combine cohesive powder as analogue for brittle anhydrites and carbonates with viscous salt analogues to explore the developing fault geometry and the resulting distribution of salt in the faults. Using the observations from analogue models as input, numerical models investigate the stick-slip behaviour of fault zones containing ductile material qualitatively with the discrete element method (DEM). Results show that the DEM approach is suitable for modelling the seismicity of faults containing salt. The stick-slip motion of the fault becomes dependent on shear loading rate with a modification of the frequency–magnitude distribution of the generated seismic events.
