Reservoir compaction, surface subsidence and induced seismicity are often associated with prolonged hydrocarbon production. Recent experiments conducted on the Groningen gas field's Slochteren sandstone reservoir rock, at in-situ conditions, have shown that compaction involves both poro-elastic strain and time-independent, permanent strain, caused by consolidation and shear of clay films coating the sandstone grains, with grain failure occurring at higher stresses. To model compaction of the reservoir in space and time, numerical approaches, such as the Discrete Element Method (DEM), populated with realistic grain-scale mechanisms are needed. We developed a new particle-interaction law (contact model) for classic DEM to explicitly account for the experimentally observed mechanisms of non-linear elasticity, intergranular clay film deformation, and grain breakage. It was calibrated against both hydrostatic and conventional triaxial
Accepted ArticleThis article has been accepted for publication and undergone full peer review but has not been through the copyediting, typesetting, pagination and proofreading process, which may lead to differences between this version and the Version of Record. Please cite this article as