The access to real geometallurgical data is very limited in practice, making it difficult for practitioners, researchers and students to test methods, models and re- produce results in the field of geometallurgy. The aim of this work is to propose a methodology to simulate geometallurgical data with geostatistical tools preserving the coherent relationship among primary attributes, such as grades and geological attributes, with mineralogy and some response attributes, for example, grindability, throughput, kinetic flotation performance and recovery. The methodology is based in three main components: (i) definition of spatial relationship between Geomet- allurgical units, (ii) co-simulation of regionalized variables with geometallurgical coherence, and (iii) simulation of georeferenced drillholes based on geometrical and operational constraints. The simulated geometallurgical drillholes generated look very realistic, and they are consistent with the input statistics, coherent in terms of geology and mineralogy, and produce realistic processing metallurgical perfor- mance responses. These simulations can be used for several purposes, for example, benchmarking geometallurgical modelling methods and mine planning optimiza- tion solvers, or performing risk assessment under different blending schemes. 1