Application of Fourier Finite Differences and lowrank approximation method for seismic modeling and subsalt imaging