El objetivo general del presente trabajo de título es definir e implementar una metodología para simulación de geoestadísticas de leyes considerando incertidumbre en el valor medio de la ley a simular atributo a estimar. En la actualidad, la simulación se condiciona a los datos por medio de un Kriging simple, el cual considera una hipótesis de estacionaridad estricta, que no se cumple a cabalidad en la mayoría de los casos.