A computer implemented method for the improved interpretation of seismic data, comprising the steps of: (a) obtaining a 3D seismic data set from a predetermined region; (b) selecting a zone of interest from said 3D seismic data set; (c) determining a frequency spectrum of said zone of interest; (d) performing a spectral decomposition into magnitude responses of a plurality of predetermined frequencies for said frequency spectrum; (e) generating a synthetic model representative of said zone of interest that is defined in at least a three dimensional space, wherein a first variable parameter is defined in a first dimension of said three dimensional space, and at least a second variable parameter is defined in at least a second and/or third dimension of said three dimensional space; (f) performing seismic modelling for said synthetic model on a frequency band that matches that of the 3D seismic data set, generating a synthetic seismic data set representing the seismic response of said synthetic model; (g) determining a frequency spectrum of said synthetic seismic data set and performing spectral decomposition of said frequency spectrum into magnitude responses of a plurality of predetermined frequencies; (h) calibrating said synthetic model utilising additional data indicative to physical properties of said zone of interest; (i) projecting said frequency spectrum of said synthetic seismic data set onto a horizon or slice of said synthetic model; and characterised in that the method further comprises the steps of: (j) generating parameter cross plot of the magnitude responses of said projected frequency spectrum with respect to the values of said first variable parameter and said at least a second variable parameter that corresponds to the projection onto a horizon or slice, and (k) identifying correlations between points of the cross plot and magnitude response values of the spectral decomposition of said zone of interest of said 3D seismic data set, and assigning, based on said identified correlations, values of said first and said at least second variable parameter to said zone of interest of said 3D seismic data set. (a) obtaining a 3D seismic data set from a predetermined region; (b) selecting a zone of interest from said 3D seismic data set; (c) determining a frequency spectrum of said zone of interest; (d) performing a spectral decomposition into magnitude responses of a plurality of predetermined frequencies for said frequency spectrum; (e) generating a synthetic model representative of said zone of interest that is defined in at least a three dimensional space, wherein a first variable parameter is defined in a first dimension of said three dimensional space, and at least a second variable parameter is defined in at least a second and/or third dimension of said three dimensional space; (f) performing seismic modelling for said synthetic model on a frequency band that matches that of the 3D seismic data set, generating a synthetic seismic data set representing the seismic response of said synthetic model; (g) determining a frequency spectrum of said synthetic seismic data set and performing spectral decomposition of said frequency spectrum into magnitude responses of a plurality of predetermined frequencies; (h) calibrating said synthetic model utilising additional data indicative to physical properties of said zone of interest; (i) projecting said frequency spectrum of said synthetic seismic data set onto a horizon or slice of said synthetic model; and (j) generating parameter cross plot of the magnitude responses of said projected frequency spectrum with respect to the values of said first variable parameter and said at least a second variable parameter that corresponds to the projection onto a horizon or slice, and (k) identifying correlations between points of the cross plot and magnitude response values of the spectral decomposition of said zone of interest of said 3D seismic data set, and assigning, based on said identified correlations, values of said first and said at least second variable parameter to said zone of interest of said 3D seismic data set.