Mineral resource evaluation requires defining geological rock-type domains. The traditional simulation methods have serious limitations for applications to large numbers of domains, which have complex contact relations. Plurigaussian simulation is an effective method which can be applied, in a simple way, to any number of domains, using both local and global geological information to infer the distributions of rock types. This work not only presents the application of the plurigaussian simulation method to the Balya lead-zinc deposit, but also assesses the spatially varying rock type proportions, and accounts for uncertainties between them. These parameters are extremely important for mining deposits, since the mineralizations of interest generally occur only in certain rock types. Furthermore, being able to model the different geological rock types is vital to good mine operations, production planning, and management. The results indicate that the plurigaussian method correctly reproduces the different orientations of the individual rock types, as seen in drill holes, and the proportion of each rock type, even if this varies in space.