We are concerned with the acoustic scattering problem, at a frequency κ, by many small obstacles of arbitrary shapes with impedance boundary condition. These scatterers are assumed to be included in a bounded domain Ω in R 3 which is embedded in an acoustic background characterized by an eventually locally varying index of refraction. The collection of the scatterers Dm, m = 1, ..., M is modeled by four parameters: their number M , their maximum radius a, their minimum distance d and the surface impedances λm, m = 1, ..., M . We consider the parameters M, d and λm's having the following scaling properties:≈ a t and λm := λm(a) = λm,0a −β , as a → 0, with non negative constants s, t and β and complex numbers λm,0's with eventually negative imaginary parts.We derive the asymptotic expansion of the farfields with explicit error estimate in terms of a, as a → 0. The dominant term is the Foldy-Lax field corresponding to the scattering by the point-like scatterers located at the centers zm's of the scatterers Dm's with λm|∂Dm| as the related scattering coefficients. This asymptotic expansion is justified under the following conditionsand the error of the approximation is C a 3−2β−s , as a → 0, where the positive constants a0, λ−, λ+ and C depend only on the a priori uniform bounds of the Lipschitz characters of the obstacles Dm's and the ones of M (a)a s and d(a) a t . We do not assume the periodicity in distributing the small scatterers. In addition, the scatterers can be arbitrary close since t can be arbitrary large, i.e. we can handle the mesoscale regime. Finally, for spherical scatterers, we can also allow the limit case β = 1 with a slightly better error of the approximation.