The invention brings forward a frequency domain three-dimensional irregular earthquake data reconstruction method. The method is characterized in that first of all, three-dimensional earthquake data in a time domain is converted to a frequency domain by use of Fourier transform, and then, a projection onto convex set (POCS) algorithm is employed and curvelet transform capable of describing localized features of the earthquake data is introduced; and in an iteration process, a new threshold parameter attenuating according to an index rule is brought forward, and each frequency slice is individually reconstructed by use of a soft threshold operator, such that the iteration frequency is reduced, the reconstruction precision is improved, and the purpose of reconstructing the three-dimensional earthquake data is realized. According to the invention, the new threshold parameter attenuating according to the index rule is brought forward and each frequency slice is reconstructed in individually by use of the soft threshold operator, such that the disadvantage of quite slow convergence speed of a conventional threshold parameter is overcome, the calculation complexity of an algorithm is reduced, the calculation efficiency is substantially improved, and the operation time is reduced.