Abstract: The paper presents a numerical approach to modelling the linear and nonlinear scattering at a layer of bubbles. Numerical studies of noncoherent scattering using a Monte Carlo method that can cope with various effects of propagation inside of a bubbles layer such as boundary reflections, scattering anisotropy, sound dispersion, sound attenuation and time dependence are performed. Based on Rayleigh-Plesset bubble’s oscillation equation, the propagation and backscattering of signals at bubble population of different size spectra is analyzed. Propagated through the layer and backscattering changes in the signal spectra are studied and modeled according to the properties of bubble populations similar as in the sea.