The invention discloses a
curvelet transform anti-
aliasing seismic
data reconstruction method based on a projection onto convex set
algorithm. The method is characterized in that aiming at
aliasing seismic data with rule lack or
aliasing seismic data with large sampling interval, firstly decomposing
curvelet transform into an f-k operator and a
curvelet splicing operator, starting from an f-k domain origin, obtaining an
energy distribution map by employing full-frequency-band information and an angle scanning policy, and picking a boundary
peak value of the distribution map, wherein the boundary represents an effective wave energy boundary; selecting a proper threshold, and establishing a corresponding
mask function according to the effective wave energy in the boundary peak; and finally introducing the
mask function to the conventional projection onto convex set
algorithm to obtain a new projection onto convex set
algorithm, and in an
iteration process, reconstructing the anti-aliasing seismic data by employing a hard threshold operator and an index threshold parameter formula. According to the method, the
signal-to-
noise ratio and the calculating speed of the reconstructed signalare increased, the weak effective wave
signal is protected, and reflecting wave events are more continuous and clear.