Modeling stress evolution around a rising salt diapir

We model the evolution of a salt diapir during sedimentation and study how deposition and salt movement affect stresses close to the diapir. We model the salt as a solid visco-plastic material and the sediments as a poro-elastoplastic material, using a generalized Modified Cam Clay model. The salt flows because ongoing sedimentation increases the average density within the overburden sediments, pressurizing the salt. Stresses rotate near a salt diapir, such that the maximum principal stress is perpendicular to the contact with the salt. The minimum principal stress is in the circumferential direction, and drops near the salt. The mean stress increases near the upper parts of the diapir, leading to a porosity that is lower than predicted for uniaxial burial at the same depth. We built this axisymmetric model within the large-strain finite-element program Elfen. Our results highlight the fact that forward modeling can provide a detailed understanding of the stress history of mudrocks close to salt diapirs; such an understanding is critical for predicting stress, porosity, and pore pressure in salt systems.     

Comparison of evolutionary and static modeling of stresses around a salt diapir

We compare an evolutionary with a static approach for modeling stress and deformation around a salt diapir; we show that the two approaches predict different stress histories and very different strains within adjacent wall rocks. Near the base of a rising salt diapir, significantly higher shear stresses develop when the evolutionary analysis is used. In addition, the static approach is not able to capture the decrease in the hoop stress caused by the circumferential diapir expansion, nor the increase in the horizontal stress caused by the rise of the diapir. Hence, only the evolutionary approach is able to predict a sudden decrease in the fracture gradient and identify areas of borehole instability near salt. Furthermore, the evolutionary model predicts strains an order of magnitude higher than the strains within the static model. More importantly, the evolutionary model shows significant shearing in the horizontal plane as a result of radial shortening accompanied by an almost-equivalent hoop extension. The evolutionary analysis is performed with ELFEN, and the static analysis with ABAQUS. We model the sediments using a poro-elastoplastic model. Overall, our results highlight the ability of forward evolutionary modeling to capture the stress history of mudrocks close to salt diapirs, which is essential for estimating the present strength and anisotropic characteristics of these sediments.