We illustrate the performance of the generalized screen propagator on real seismic data for 3D zero-offset and prestack depth imaging. We use TotalFinaElf's L7D data set from the North Sea, a 3D marine seismic survey that contained limited azimuthal coverage. The subsurface shows significant tectonic deformation, including an intrusive salt body in sedimentary sequences. A transformation to common azimuth is applied prior to the 3D prestack depth imaging procedure. We compare the performance of the generalized screen propagator with that of a hybrid phase shift plus interpolation (PSPI)/split-step Fourier method. Three-dimensional prestack results confirm the generalized screen method handles multipathing more accurately. Comparisons are also made with Kirchhoff migration results. The results differ mainly in the fine-scale irregularities of the image and not in the wavefront set of the image. Using synthetic models of similar structure (the SEG/EAGE salt model), we further illustrate the importance of multipathing and multiple scattering. Overall, our results show that our wave-equation approach produces better images than the Kirchhoff approach to prestack depth migration; we attribute this mainly to a more complete handling of wave diffraction in the generalized screen expansion, which becomes important in strongly heterogeneous and irregular velocity models such as the ones containing salt bodies.