A thermo-poro-mechanical constitutive and numerical model for deformation in sedimentary basins

Comprehensive modeling of sedimentary basins should integrate the coupled nature of different phenomena such as sediment compaction, pore-fluid flow and heat transport. In this context, the thermal evolution of the basin plays an important role in the mechanical and chemo-mechanical deformation processes as heat modifies fluids viscosity and minerals physicochemical properties, thus affecting pore fluid expulsion and mineral stability. The present paper describes a three-dimensional constitutive and numerical model specifically devised for dealing with thermo-poromechanical deformation processes during diagenesis. The sedimentary basin is modeled as a fully saturated thermo-poro-elastoplastic-viscoplastic medium undergoing large strains. A key feature of the model is related to the evolution of the sediment material properties associated with temperature and large irreversible porosity changes. The computational model, which integrates both poroplastic and poroviscoplastic components of deformation at large strains, relies upon a parallel implementation of the finite element method with a shared memory multiprocessing interface. Numerical simulations of gravitational compaction during the sediment accretion phase as well as along geological period of pore pressure dissipation are performed in the context of oedometric setting. Special emphasis is given to temperature effects on the deformation history of the basin.