Patent ID: 8855934

Claim:
A method of processing seismic data comprising: receiving the seismic data; and using a processor to process the seismic data, wherein the processor comprises a non-transitory computer-readable medium comprising instructions that when executed by the processor cause the processor to perform that following method: transforming the seismic data into an f-x domain using a discrete Fourier transform; at each discrete frequency, applying an infinite impulse response (IIR) filter; applying each IIR filter spatially in both a forward direction and a backward direction, wherein initial conditions are set for the IIR filter at the beginning of the data to minimize transients before applying the IIR filter in the forward direction; assuming an input step function and setting filter coefficients such that the input and output functions of the filter are the same, wherein at least one of a first sample is used as the step size of the step function for setting the initial conditions at the beginning of the data and a last sample is used as the step size of the step function for setting the initial conditions at the end of the data, and wherein: a noise attenuation method is applied to the first sample; the noise attenuation is only applied to the first sample; the attenuated first sample is used as the step size of the step function for setting the initial conditions at the beginning of the data; and the step size used for setting the initial conditions at the end of the data is the last sample.