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.