The invention discloses a three-dimensional seismic
data reconstruction method based on a weighted fast iterative shrinkage threshold
algorithm, and the method aims at solving a problem that a conventional three-dimensional seismic
data reconstruction method neglects the potential correlation between adjacent slices, and comprises the steps: firstly dividing three-dimensional seismic data into a series of time slices; and carrying out multi-scale multi-direction two-dimensional
curvelet transformation to obtain a
curvelet domain coefficient, then introducing a fast iterative shrinkage threshold
algorithm, and carrying out reconstruction by adopting an index threshold parameter formula. In the process, according to the characteristic that
curvelet coefficient support sets of two adjacent time slices have a large intersection,
a weighting operator is constructed in a curvelet domain by using a prior support set provided by a previous reconstructed time slice, and then the weighting operator is fused into a fast iterative shrinkage threshold
algorithm. According to the method, the
signal-to-
noise ratio of the reconstructed
signal is increased, the calculation speed is increased, the requirement for a
computer memory is reduced, weak effective wave signals are protected, and the reflected wave event is more continuous and clearer.