Abstract: Data of vertical seismic profile in the geological formation are recorded and at least one estimation of a geological formation local model dip is performed, also at least one dynamic dip estimation is performed for at least one wavetype excited by sources and registered by receivers using directionally-based vectors. A difference between the local model dip estimation and the dynamic dip estimation is calculated. Based on the calculated difference between the local model dip estimation and the dynamic dip estimation at least one weighting coefficient is calculated. The calculated at least one weighting coefficient is used for determining at least one formulae for imaging with the use of elastic waves at reverse time migration conditions taking into account a dip angle of the geological formation. The determined at least one formulae is applied to the obtained vertical seismic profile data and at least one image of the geological formation is produced.