We model P and S-wave mean square envelopes in 2-D random media by solving the Radiative Transfer Theory (RTT) equation using a Monte Carlo (MC) particle based approach. Using a source frequency in the order of kilohertz, typical of cross-well seismic sources, we find good agreement with finite differences in both effective and forward scattering regimes. We show that MC simulations is a much faster way of modeling scattering attenuation, compared to finite differences, and therefore is a suitable tool to characterize the medium heterogeneities associated, for example, with highly fracture reservoirs.