Seismic imaging is the final stage of the seismic processing allowing to reconstruct the internal subsurface structure. This procedure is one of the most time consuming and it requires huge computational resources to get high-quality amplitude-preserving images. In this paper, we present a parallel algorithm of seismic imaging, based on the solution of the one-way wave equation. The algorithm includes parallelization of the data flow, due to the multiple sources/receivers pairs processing. Wavefield extrapolation is performed by pseudo-spectral methods and applied by qFFT - each dataset is processed by a single MPI process. Common-offset vector images are constructed using all the solutions from all datasets thus by all-to-all MPI communications.