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 in executable code in memory of a computer system with a visco-acoustic wave equation solved using a pseudo-analytical method, wherein the visco-acoustic wave equation for media with heterogeneous attenuation is time discretized by a second-order finite difference approximation that incurs errors; obtaining seismic data for a zone of interest, wherein the seismic data is acquired using a seismic energy source and a plurality of towed seismic sensors; creating a model for the zone of interest by the reverse time migration module in the executable code in the memory of the computer system, wherein the model comprises spatial variability in reference velocity and quality factor; performing reverse time migration using the reverse time migration module in the executable code in the memory of the computer system to obtain migrated seismic data for the zone of interest; and generating and displaying an image of the migrated seismic data by the computer system.