Systems and methods for imaging seismic data

A non-transitory program carrier device tangibly carrying computer executable instructions for imaging seismic data, the instructions being executable to implement:

• approximating TTI coefficients using a Pade approximation and a dispersion relation equation;
• applying hybrid TTI-WEM to a velocity model and anisotropic parameters for a pre-stack shot gather using the approximated TTI coefficients and a computer system to determine a source side wavefield propagation value and a receiver side wavefield propagation value in a frequency-space domain;
• converting the source side wavefield propagation value and the receiver side wavefield propagation value from the frequency-space domain to a time-space domain;
• and
• applying a zero-lag cross-correlation image condition equation to form a partial output image using the converted source side wavefield propagation value and the converted reservoir side wavefield propagation value.