Due to the existence of a constant sum of constraints, the geochemical data is presented as the compositional data that has a closed number system. A closed number system is a dataset that includes several variables. The summation value of variables is constant, being equal to one. By calculating the correlation coefficient of a closed number system and comparing it with an open number system, one can see an increase in the values of the closed number system, which is false. Such features of this data prevent the application of standard statistical techniques to process the data. Therefore, several methods have been proposed for transforming the data from closed to open number systems. There are various geostatistical methods consisting of estimation and simulation methods in order to model a deposit. Geostatistical simulations can produce various models for a deposit with different probability percentages. The most applicable geostatistical simulation method is the sequential Gaussian simulation technique, which is highly flexible. In this work, 392 Litho-geochemical data of the Baghqloom region of Kerman in Iran consisting of 20 elements were at first converted using an open number system. Afterwards, the elements that were helpful for exploring the area and were normally standard were simulated for 100 times. After the simulations, the valid output was chosen using geostatistical validation. The maps derived from the simulations revealed the enriched concentrations of mineralization elements in the central regions.