The invention discloses a three-dimensional VSP seismic data interpolation method based on weighted
matching pursuit, comprising the following steps: S1, respectively performing static correction,
deconvolution and VSP dynamic correction
processing on original seismic data collected in the field, eliminating the
time difference between traces of different offset distances, and obtaining the three-dimensional VSP seismic data before interpolation; S2, obtaining the spatial sampling position of the three-dimensional VSP seismic data according to the
shot point coordinates of the three-dimensional VSP seismic data before interpolation; S3, encrypting the three-dimensional VSP seismic data before interpolation in the
spatial direction with blank data, and obtaining a data body with spatial sampling rules; S4, performing F-K transformation on the data body with spatial sampling rules obtained in step S3, and obtaining an F-K spectrum; S5, performing linear
Radon transformation on the F-K spectrum, and obtaining
a weighting operator according to the
Radon transformation result; S6,
processing the F-K spectrum with the weighting operator obtained in step S5 to obtain a new F-K spectrum; S7, performing F-K inverse transformation on the new F-K spectrum to obtain the three-dimensional VSP seismic data after interpolation.