Single scattering events recorded in surface seismic data do not fully illuminate the subsurface structure, especially if it is complicated. In such cases, multiple internal scatterings (internal multiples) can help improve the illumination. We devise a generalized internal multiple imaging (GIMI) procedure that maps internal multiple energy to their true location with a relatively mild addition to the computational cost. GIMI theory relies heavily on seismic interferometry, which often involves cumbersome algebra, especially when one is dealing with high-order terms in the perturbation series. To make the derivations, and inference of the results easier, we introduce Feynman-like diagrams to represent different terms of the perturbation series (solution to the Lippman–Schwinger equation). The rules we define for the diagrams allow operations like convolution and cross-correlation in the series to be compressed in diagram form. The application of the theory to a double scattering example demonstrates the power of the method.