(8) |

Often with seismic data, the amplitude varies across offset, midpoint and time. These amplitude variations can be troublesome with least-squares inversion because they tend to bias the final result Claerbout (1992). Therefore, it is important to make sure that these amplitude variations do not affect our processing. One solution is to apply a weight to the data like Automatic Gain Control (AGC) or a geometrical spreading correction. However, a better way is to incorporate this weight inside the inversion by scaling the residual Guitton (2003a). Introducing a weighting function in the PEFs estimation leads to:

(9) |

(10) |

(11) |

Now, prior to the signal estimation in equation (5), and need to be computed from a signal and noise model, respectively. The multiple model is often (but not necessarily) derived by auto-convolving the recorded data Verschhur et al. (1992), thus obtaining a prestack model of the multiples later used to estimate a bank of non-stationary PEFs .

At this stage, a key assumption is that the relative amplitude of all order of multiples is preserved. In theory, an accurate surface-related multiple model can be derived if (1) the source wavelet is known, (2) the surface coverage is large enough, and (3) all the terms of the Taylor series that model different orders of multiples are incorporated Verschhur et al. (1992). In practice, however, a single convolution is performed (first term of the Taylor series) which gives a multiple model with erroneous relative amplitude for high-order multiples Guitton et al. (2001); Wang and Levin (1994). In addition, the surface coverage might not be sufficient. This can generate wrong amplitudes for short offset traces and complex structures. Because PEFs estimate patterns, wrong relative amplitude can affect the noise estimation. However, as we shall see later, 3D filters seem to better cope with noise modeling inadequacies.

The signal PEFs are more difficult to estimate since the signal is usually unknown. However, Spitz (1999) shows that for uncorrelated signal and noise, the data PEFs can be approximated with

(12) |

(13) |

(14) |

(15) |

Thanks to the Helix Claerbout (1998); Mersereau and Dudgeon (1974), the PEFs can have any dimension. In this paper, I use 2D and 3D filters and demonstrate that 3D filters lead to the best noise attenuation result. When 2D filters are used, the multiple attenuation is performed on one shot gather at a time. When 3D filters are used, the multiple attenuation is performed on one macro-gather at a time. A macro-gather is a cube made of fifty consecutive shots with all the offsets and time samples. When the multiple attenuation is done, the macro-gathers are reassembled to form the final result. There is an overlap of five shots between successive macro-gathers. In the next section, I show a prestack multiple attenuation example with the synthetic Sigsbee2B dataset.

Figure 1

Figure 2

5/23/2004