This work presents a simulation of the response of packets of microbubbles in an ultrasonic pulse-echo scan line. Rayleigh-Plesset equation has been used to predict the echo from numerically obtained radial dynamics of microbubbles. Varying the number of scattering microbubbles on the pulse wave form has been discussed. To improve microbubble-specific imaging at high frequencies, the subharmonic and second harmonic signals from individual microbubbles as well as microbubbles packets were simulated as a function of size and pressure. Two different modes of harmonic generation have been distinguished. The strength and bandwidth of the subharmonic component in the scattering spectrum of microbubbles is greater than that of the second harmonic. The pressure spectra provide quantitative and detailed information on the dynamic behaviour of ultrasound contrast agent microbubbles packet.