El conocimiento de la arquitectura estratigráfica y las propiedades geotécnicas de los sedimentos superficiales del suelo es esencial para la evaluación del riesgo geotécnico. En el área de estudio de Saguenay, los depósitos cuaternarios consisten en una capa de till basal y depósitos postglaciares heterogéneos. Considerando el contexto estratigráfico y la heterogeneidad del tipo de suelo, se desarrolla una metodología estocástica de múltiples pasos para la modelización geológica en 3D y la cuantificación de las incertidumbres asociadas. Esta metodología se adopta para estudios regionales e involucra métodos de interpolación y simulación geoestadística. Se aplica el kriging bayesiano empírico (EBK) para generar el mapa de topografía de la roca madre y determinar el grosor de los sedimentos de till y sus incertidumbres. La media y varianza localmente variables del método EBK permiten tener en cuenta la complejidad de los datos y la no estacionaridad moderada. Luego se realiza una simulación de indicadores secuenciales para determinar la probabilidad de ocurrencia de los sedimentos postglaciares discontinuos (arcilla, arena y grava) sobre la capa de till basal. Los mapas de grosor individuales de las capas de suelo discontinuas y sus incertidumbres se generan de manera probabilística. El marco estocástico propuesto es adecuado para depósitos de suelo heterogéneos caracterizados por conjuntos de datos superficiales y subsuperficiales complejos.