In this research, the spatial pattern of seismic activity in the Zagros region has been investigated using Variogram analysis method. For this purpose, based on the analysis of the seismicity data of the region, the spatial variations of the density of seismic energy released, have been investigated and the Variograms representing these spatial variations in different directions, have been provided and interpreted. The results obtained show that across the directions close to NW-SE (nearly parallel to the general trend of the Zagros belt), the least spatial variations and in contrast in directions nearly perpendicular to the axis, the highest degree of spatial variations and spatial heterogeneities are observed. The results also indicate that geostatistics methods, and in particular Variogram modeling, can be very useful in evaluating spatial variations associated with seismic activity of regions.