The discrete equivalent spherical scatterers aggregate model is developed to investigate the photon′s transport characteristics in highly scattering media. Based on Monte-Carlo method, we propose a new ellipse algorithm to determine the scattering point distributions. According to Mie theory and the radiative transfer equation, combined with the obtained scattering point distributions, path-length dependent scattering order distribution and scattering intensity in different scattering orders with different particle sizes are numerically calculated. The results show that the photon's transport in highly media is essentially a multiple scattering process. The numerical algorithm can be used effectively to analyze dynamic scattering properties of scattered light from single-scattering to multiple-scattering regime.