摘要 |
The present invention relates to a method of filtering seismic data for noise attenuation. An embodiment of the present invention provides a method of processing seismic data in which the seismic data is transformed into an f-x domain (frequency space domain) using a discrete Fourier transform and is then filtered at each discrete frequency using an infinite impulse response (IIR) filter, such as a Butterworth or Chebychev filter. Initial conditions for the filter are set to minimize transients before applying the filter in both a forwards or backwards direction. Coefficients can be set so that the input and output functions of the filter are the same. Noise attenuation methods and methods for choosing initial step functions can be used when setting initial conditions. Preferably, extrapolated data is added before application of the filter. |