Los suelos residuales y saprolitos son frecuentes en yacimientos lateríticos de níquel, se caracterizan por la heterogeneidad de su configuración y la variabilidad de las propiedades mecánicas, mayormente debido al proceso de formación de la roca madre y a las etapas muy erráticas de meteorización. Este trabajo propone una metodología que permite el ingreso de la cohesión como una variable aleatoria en zonas con estratigrafía similar. Se desarrolló un código de automatización de la herramienta Plaxis para el ingreso de los datos, ejecución repetitiva de la rutina de cálculo y la gestión de los datos de salida para detectar la superficie de falla más probable. Se presenta un caso de aplicación al yacimiento Fenix que la Compañía Guatemalteca de Níquel explota en Izabal, Guatemala. Los resultados se integraron en un análisis de riesgo de la operación minera que permitió la toma de decisiones y la optimización de las operaciones.
