For the study of complex sound scattering phenomena, the Monte Carlo method was investigated to establish a novel particle size characterization model in concentrated particulate two-phase flow. Starting from the single particle scattering theory, the numerical simulation was carried out to predict the acoustic attenuation characteristics in the two-phase system of glass bead/water with different particle sizes, concentrations and ultrasonic frequencies. Through the comparison of ECAH model of single scattering theory, it was found that they met each other well at the concentration of 1%(φ), while clear divergence occurred with the increase of the concentration. At the same time, the simulation results were compared with experiment result from literature. It was shown that the attenuation in the particle radii of 35 and 90 μm in two-phase system predicted by Monte Carlo method agreed well with experiments. The proposed Monte Carlo method could be used for calculating the attenuation in high particle volume concentration of 50% and numerical modeling of particle size measurement in the complex scattering condition.