(1) |

where A, B and C are parameters related to the average and contrasts of the elastic properties of the limiting layers. Therefore, if it is possible to determine these parameters, one can simulate the amplitudes of primary events at every depth or time step. I refer to the simulation of the primary amplitudes as AVA modeling. The aim of this work is not to invert the AVA parameters for rock and fluid properties, but just to obtain a reasonable estimate of the primaries to be further used in a residual-multiple-attenuation sequence.

The AVA modeling problem can be addressed by solving locally, for every different depth (or time) step, the normal equations ,where are the estimated AVA parameters, has the column vector form , , and is the ADCIG (or a CMP gather after NMO). The simulated primaries are obtained by using the estimated AVA parameters in equation (1). One problem that arises from this local solution is that the AVA parameters are not constrained to be smooth along depth. This can lead to anomalous amplitudes in the simulated primaries for a certain depth.

In the present approach, the AVA modeling problem consists of determining the three parameters A, B and C of equation (1) by solving the following data fitting problem:

(2) |

(3) |

Generally, in ADCIGs, residual multiples are more persistent at near angles, because of insufficient moveout difference between them and the flattened primaries. Additionally, at the farthest angles, stretch occurs. To avoid these imperfections in the input data, the fitting goal is to minimize the residual,

(4) |

If unrealistic variations of the AVA parameters with depth are an issue, the following regularization goal can be introduced:

(5) |

To accelerate the solution, the problem can be solved with preconditioning Claerbout and Fomel (2001), using the transformation , where and is the preconditioned variable. Finally, the fitting goals reduce to

(6) |

After the determination of the three AVA parameters for every depth step, primaries are simulated by computing the reflection coefficient for all reflection angles using equation (1). Of course, the method relies on the flatness of the reflectors in the CIG to correctly extract the 3 parameters of the AVA curve.

To get the residual-multiple-attenuated data, the adaptive subtraction must be applied in two steps. The first one aims to obtain an estimate of the residual multiples. This is done by subtracting the adjusted version of the simulated primaries from the original data. Amplitude and phase adjustments are achieved by the convolution of the simulated primaries with a prediction-error Wiener filter.

As the estimated residual multiples may contain some primary information (mainly at small reflection angles), I use the strategy proposed by Guitton et al. (2001), the so-called `subtraction method,' in which two different prediction-error filters (PEFs), which model primaries and multiples, are computed. Their method allows regularization, to decrease the crosstalk between multiples and primaries. The corresponding fitting goals are

subjected to

5/6/2007