Several non-destructive evaluation techniques to characterize concrete structures are based on ultrasonic wave propagation. The interpretation of the results is often limited by the scattering phenomena between the ultrasonic wave and the high concentration aggregates contained in the cement matrix. Numerical simulations allow for further insights. This study aims to build a two-dimensional numerical model in order to reproduce and interpret ultrasonic wave propagations in concrete. The model is built in a spectral-element software package called SPECFEM2D. The validation of the numerical tool is based on the use of resin samples containing different amount of aluminum rods from low (5%) to high concentration (40%), the last one being representative of aggregate concentration in concrete. These samples are characterized using an ultrasonic testing bench (ultrasonic water tank) from 150kHz to 370kHz. The measured results are analyzed in terms of phase velocity and attenuation which are the main parameters of coherent waves. As homogenization models such as the Waterman-Truell or Conoir-Norris models are usually used to model coherent waves in two-phase systems, we also compare the experimental and numerical results against them. We confirm that the use of these homogenization models is limited to low concentration of scattering phase, which is not adapted to applications to concrete. Finally, we use our numerical tool to carry out a parametric study on scatterer concentration, shape, orientation and size distribution of aggregates in concrete. We show that aggregate orientation has an influence on coherent wave parameters, but aggregate shape has not.