The anomalies detected in phosphatic series of Sidi Chennane, one of phosphatic basins in Morocco, hinder the proper exploitation of phosphate levels and the assessing phosphate reserves seems incorrect. The purpose of this study was the evaluation of geostatistical methods for delimitation of these disturbances. To cover all the zones being able to be disturbed, we carried out, during the geophysical prospection in a parcel of 50 ha, 5151 resistivity measurements as horizontal profiling using the well-known Schlumberger array, in order to map the spatial distribution of the sterile hardpan inclusions. Geostatistical tools were used to quantify the spatial correlation between apparent resistivity data. Semivariograms were obtained using a classical Matheron semivariogram estimator and fit to the experimental semivariograms obtained. We have selected those with the best fit in terms of sum of squared residuals (SSR). Geostatistical analysis was performed using the software VESPER 1.63. Spatial distribution maps were made by ordinary kriging, the qualitative interpretation of these maps reflects that the exponential model is found to be the best model representing the spatial variability of our geoelectric data. The qualitative interpretation of the kriged resistivity maps allows defining resistivity contrast, consequently we have delimited the crossing dominate area from a ¿normal¿ into a ¿disturbed¿ area. Models of the geology were successfully obtained from geostatical method, which help mapping the phosphate deposit inclusions and the estimations of phosphate reserves were improved and better constrained.