Towards the application of Stokes flow equations to structural restoration simulations

Abstract

Structural restoration is commonly used to assess the deformation of geological structures and to reconstruct past basin geometries. For this, geomechanical restoration considers faults as frictionless contact surfaces. To bring more physical behavior and better handle large deformations, we build on a reverse-time Stokes-based method, previously applied to restore salt structures with negative time step advection. We test the applicability of the method to structures including sediments of variable viscosity, faults and non-flat topography. We present a simulation code that uses a combination of arbitrary Lagrangian-Eulerian methods and particle-in-cell methods, and is coupled with adaptive mesh refinement. It is used to apply the reverse-time Stokes-based method on simple two-dimensional geological cross-sections and shows that reasonable restored geometries can be obtained.