A three-dimensional (3-D) prestack reverse-time migration algorithm for common-source P-wave data from anisotropic media is developed and illustrated by application to synthetic data. Both extrapolation of the data and computation of the excitation-time imaging condition are implemented using a second-order finite-difference solution of the 3-D anisotropic scalar-wave equation. Poorly focused, distorted images are obtained if data from anisotropic media are migrated using isotropic extrapolation; well focused, clear images are obtained using anisotropic extrapolation. A priori estimation of the 3-D anisotropic velocity distribution is required. Zones of anomalous, directionally dependent reflectivity associated with anisotropic fracture zones are detectable in both the 3-D common-source data and the corresponding migrated images.