Summary The complex sedimentology and the irregular erosion process that affected the thickness of the N sandstone make it very difficult to generate a mathematical representation of the reservoir. Therefore the integration of all the available data is the key to represent the spatial variability of reservoir properties. Cohembi wells have been drilled in the best zones of the N sandstone just representing the sands and not showing the non-reservoir rocks. Thus all the produced variograms are biased by a dataset where reservoir heterogeneities cannot be represented. It is well known that baffles and barriers produce partial communication among some wells. Then geostatistics needed to be co-krigged with any variable that can be effectively used to estimate thickness changes along the reservoir. A correlation between amplitude and thickness was found. This is how a geobody was extracted from the 3D seismic. A cross plot of the geobody amplitude values vs thickness was obtained and the resulting equation from the linear regression was used to generate a “Pseudo Thickness Distribution Map” The purpose of this map is to allow a 2D trend to be used when co-krigging with the variograms and to create a 3D grid that accurately represent sand body distribution.