The quantitative analysis of Raman spectroscopic signals in biological tissue is generally difficult. Typical samples contain a multitude of molecular species and, in addition, measurements are altered by attenuation of the Raman signal. Realistic numerical modeling of the Raman process can help to facilitate the quantitative analysis of the Raman spectra, but approaches so far are scarce and often time-consuming. In this work, we report on two different and very efficient approaches for modeling of Raman scattering in turbid media irradiated by laser light. Both approaches utilize the Monte Carlo method to simulate the Raman scattering process. We compare the efficiency of both approaches and discuss possible future extensions and experimental validation.