The salt dome problem: a multilayered approach

Abstract

A two-dimensional finite element program has been developed to simulate the rise of individual diapirs through multi-layered overburdens exhibiting Newtonian-viscous rheologies. This program includes the following features: arbitrarily shaped interfaces between different fluid boundaries, layers of finite thickness and subsidence of the model within the model-frame as well as time- and depth-dependent material properties. This paper is concerned with single diapirs intruding upbuilt multilayered overburdens, the effect of a subsiding salt base due to pure downbuilding, and differential loading of gravitationally unstable distribution of mass. Variable density and viscosity of the overburden with depth control the shape and ascent rates of diapiric bodies. The subsidence history of "salt structures' dominates the development of maximum horizontal pressure gradients with time and thus the ascend histories of diapirs. -from Authors