A Blended Transient/Quasistatic Lagrangian Framework For Salt Tectonics Simulations With Stabilized Tetrahedral Finite Elements

We propose a Lagrangian solid mechanics framework for the simulation of salt tectonics and other large-deformation geomechanics problems at the basin scale. Our approach relies on general elastic-viscoplastic constitutive models to characterize the deformation of geologic strata, in contrast with the majority of published works on the subject, which utilize nonlinear Stokes flow models. By means of multiscale asymptotics, we also show that the inertia term in the momentum balance equation can be safely neglected, if the goal is to track the Earth's crust deformation over long periods of time. Our time integration strategy is a blended transient/quasistatic approach, in that it consists of a constitutive stress update, subject to the constraint that the stresses must satisfy static equilibrium. In addition, we use stabilized finite element methods specifically built for triangular and tetrahedral grids, which can also perform well under incompressibility constraints. Our approach offers computational geologists the following advantages: (1) improved flexibility in the choice of subsurface constitutive models with respect to the nonlinear Stokes flow; (2) improved efficiency over transient dynamics algorithms used in this context in the past, which are forced to resolve seismic events over geologic time scales; and (3) improved robustness in large strain computations over quadrilateral/hexahedral finite elements. We demonstrate the performance of the proposed approach with simulations of passive diapirism.