Formulas for 3D migration of seismic data in frequency-wave number domain are derived. The software and procedure used on medium-sized computer and minicomputer are given。The rationality of this method is tested by theoretical model, The 3D theoretical model, 3D theoretical migration section, 3D migrated seismic section and their seiscrop sections are presented, The actual computation and the comparison between 2D migration and 3D migration show that 3D migration not only has higher accuracy and higher SjN ratio but also effectively suppresses side waves.