Method of obtaining an attenuation model estimate for a vertical seismic profile(VSP), comprising: receiving a vertical seismic profile dataset, the VSP dataset having been generated by recording a wavefield at a plurality of depth levels;building an estimate of a plurality of changes in an attribute of the wavefield sensitive to attenuation between respective pairs of depth levels based on the VSP dataset;producing a plurality of corrected changes in the attribute between the respective pairs of depth levels by modeling local interference of the wavefield from interfaces near each of the depth levels using information comprising measured well log and/or borehole information and correcting the estimated changes in the attribute for this interference;choosing and fitting an attenuation law to the corrected transfer functions;outputting an attenuation model.