A fully-coupled geomechanical approach to capture the true compaction occurring in an evolving 3D stress regime is adopted. 2D models of clay columns accounting for sedimentation and tectonic activity are developed. A high deviatoric stress is generated in compressional tectonic basins resulting in a substantial decrease in porosity. The more the deviation from the uniaxial burial trend, the more inaccurate the EDM overpressure predictions.