This paper proposes an advanced stochastic workflow to jointly model sedimentary facies and diagenesis. The formation of interest is the Early Eocene Alveolina Limestone Formation, which outcrops in the Serraduy area (Graus–Tremp Basin, NE Spain). Ten sedimentary lithotypes representing facies or facies associations of a mixed siliciclastic–carbonate ramp system were identified within the succession. A 3D model describing the depositional architecture is also proposed. The results from the diagenetic study evidenced the occurrence of several successive calcite cements, which were grouped into five diagenetic imprints for modeling. These imprints were then quantified to ease their integration into numerical models. The following step consisted in building a 3D gridded model with seven different modeling units. They were populated using a bi-plurigaussian simulation approach that reproduced both the sedimentary organization and the observed diagenetic imprint distributions. Last, the simulation results were validated referring to paleogeographic and diagenetic conceptual maps.