摘要

In seismic processing, the problems of low or irregular space sampling are solved by seismic interpolation techniques. In this paper, we propose an edge-preserving seismic interpolation method. For a 1D signal, many signal segments are extracted by a processing window sliding along the signal sequence. For one missing data sample, there are a lot of signal segments, including or closing to it, which can be used to recover it. These signal segments are approximated by polynomials, in order to preserve the edges, the one with the minimum fitting error is used to estimate the missing data sample. For 2D seismic data, a lot of 1D signals are extracted along a specified direction firstly, and then an interpolation result is obtained by applying above 1D edge-preserving interpolation method on all these 1D signals. For one possible direction, one interpolation result is achieved. For one missing data sample in this 2D data, among the above interpolation results, the one with the minimum fitting error is selected as the final output. Applications with synthetic and real data sets show that the proposed method is edge-preserving, anti-aliasing, and can be used to reconstruct irregular seismic data effectively.