The dispersive phase velocity of a wave propagating through a system is an important
parameter and carries valuable information in non-destructive tests related to multilayered systems such as a soil site. The dispersive phase velocity of a wave can be determined using the phase spectrum, which is easily evaluated through the cross power spectrum. However, the phase spectrum as determined using the cross power spectrum is sensitive to background noise which always exists in the field. This causes difficulties in the determination of the dispersive phase velocities. In this paper, a new method to evaluate the phase spectrum using the harmonic wavelet transform is proposed. The proposed method can successfully remove background noise effects. To evaluate the validity of the proposed method, numerical simulations of multi-layered systems were performed. Phase spectrums by the proposed method were found to be in good agreement with the actual phase spectrums under conditions characterized by heavy background noise. This shows the potential of the proposed method.