Digital simulation of stochastic wind velocity field is prerequisite to flutter or buffeting analysis of long-span bridges in time-domain. Wavelet method was applied to the simulation of the stochastic wind field, performing the effective simulation of spatial stochastic wind field. According to the idea of multi-resolution analysis and orthonormal wavelets, relationship between wavelet coefficients and PSD(power spectral density)function is deduced , then the wavelet coefficients on each and every scale are obtained from the given PSD function, the intermittency was introduced into wavelet coefficients while at the same time preserves the target spectral characters. As a result, the inverse wavelet transform is applied to generate stationary wind velocity history through the given wavelet coefficients. A numerical example to illustrate the application of the proposed method for the simulation of a bridge wind velocity ﬁeld is provided.