Gaussian Packet describes the propagation of seismic wave in the high-frequency asymptotic theory whose
envelopes at given time are nearly Gaussian functions. We prove that the wave field generated by a given
Gabor function can be calculated by integrating Gaussian packets over take-off angles. The seismic source
function can be decomposed into a series of Gabor functions, and we presented a new method which can
simulate the wave field generated by any source function using Gaussian packets. Also, this method can be
used in reverse time migration. Numerical experiments of prestack depth migration demonstrate the validity
of Gaussian Packets modeling.