Geostatistical simulations are the best practice for uncertainty quantification in non-linear systems. Simulations are currently undertaken on an ad hoc basis in the mining industry due to high requirements in computing environments and advanced geostatistics. In contrast to the estimation process, multivariate simulations can be performed using different methods, tools, and computational environments, similar to data science. For large companies with multiple commodities it is important to have a standard framework to test, build, and run workflows using different methods, applications and computational environments. This work presents an auditable, modular and scalable framework to create geostatistical simulation workflows where the workflow is abstracted from the implementation. Key aspects: automation – workflows should be able to run automatically given updated data (if hands on modelling is not required); auditability – workflows should document key considerations and parameters; modularisation – workflows should permit swapping out components, for example substituting domaining in Leapfrog software with automated domaining using machine learning scripts; scalability – easy, repeatable, portable deployments on diverse infrastructure (for example, experimenting on a laptop, then moving to an on premise cluster or to the Cloud); reproducibility – parameters should be archived to permit re-running and reproduction of workflows in future years; and sensitivity analysis – any component of a workflow should be able to run with multiple scenarios of parameters and data to enable global uncertainty quantification and sensitivity analysis. The presented framework can be used to create different workflows for geostatistical simulations in multiple assets to quantify uncertainty in the value chain for resource categorisation, reserves uncertainty, value of information and robust mine planning amongst others.