The invention relates to a seismic wave attenuation compensation method in a curvelet domain. The method comprises the following steps of (1) collecting a seismic signal and carrying out preprocessing to obtain a time domain signal x(t) of an earthquake, (2) carrying out curvelet transform on the time domain signal x(t) to obtain a curvelet domain signal psi j, l (omega j, Tk), (3) carrying out smoothing process on the curvelet domain signal psi j, l (omega j, Tk) to obtain a curvelet coefficient Phi j, l (omega j, Tk), (4) carrying out threshold iterative processing on the curvelet coefficient Phi j, l (omega j, Tk) to obtain a denoised signal Psi j, l (omega j, Tk), (5) carrying out pointwise recursion to obtain the deep and shallow layer energy ratio delta j, ll (omega j, Tk) of the denoised signal Psi j, l (omega j, Tk) in each frequency band compensation angle, (6) calculating all frequency band compensation factor coefficient matrices Dj, ll (omega j, Tk) in a compensation angle, (7) weighing the reciprocal values of the compensation factor coefficient matrices Dj, ll (omega j, Tk) in all frequency band compensation angles to obtain a signal Cj, l (omega j, Tk) after curvelet wave weighing, (8) carrying out inverse curvelet transform on the signal Cj, l (omega j, Tk) after the curvelet wave weighing to obtain a time domain signal xx(t) after frequency division directional compensation. The method can be widely applied to the seismic data processing.