Seismic wave propagation in attenuative media can be efficiently modeled with ray-based methods. Here, we present a methodology to generate reflection data from attenuative anisotropic media using the Kirchhoff scattering integral and summation of Gaussian beams. Green's functions are computed in the reference elastic model by Gaussian-beam summation, and the influence of attenuation is incorporated as a perturbation along the central ray. The reflected wavefield is obtained by substituting the approximate Green's functions into the Kirchhoff scattering integral. Numerical examples for P-waves in transversely isotropic (TI) media with a horizontal reflector demonstrate the accuracy of the method.