Stochastic long-term mine planning has evolved to account for different sources of uncertainty. Typically, the uncertainty and local variability of boundaries in geological domains have been overlooked by experts through their deterministic interpretation of available data. Categorical attributes are used to model geological domains, and their stochastic simulation accounts for the mentioned issues. The ability of two-points simulation methods to reproduce complex patterns or the requirement of a training image in multiple-points simulation methods has limited their implementation in mining environments. The high-order simulation of categorical attributes presents a mathematically consistent framework that overcomes these limitations by using high-order spatial statistics from sample data. The case study at a gold mining complex shows two stochastic mine plans based on two sets of geological realisations: geological domains in the first set are modelled using conventional wireframes, while, in the second, they are simulated through the high-order method. The resulting mine plans are substantially different; while both plans present a similar quantity of metal recovered and lifespan, risk profiles are up to 40% wider, and the expected NPV is 20% higher for the case of simulated geological domains, given the decrease of waste handling costs and the corresponding reduction in environmental footprint.