1. A method of seismic imaging using reverse time migration with compensation for visco-acoustic effects, the method comprising: configuring a reverse time migration module with a visco-acoustic wave equation solved using a pseudo-analytical method; obtaining seismic data for a zone of interest; creating a model for the zone of interest, wherein the model comprises spatial variability in reference velocity and quality factor; performing reverse time migration using the reverse time migration module to obtain migrated seismic data for the zone of interest; and generating and displaying an image of the migrated seismic data.