Nós discutimos e aplicamos um método para predizer e atenuar múltiplas em dados sísmicos não empilhados. Para meios estratificados horizontalmente as múltiplas são periódicas em relação ao tempo de interseção no domínio tau-p. Assim, elas podem ser estimadas e removidas em cada traço. Mas, este não é o caso em meios 2D. Para meios 2D implementamos o método bidimensional, no domínio das ondas planas, onde a variação lateral da estrutura é considerada pelo acoplamento de diferentes parâmetros de raio. As múltiplas são estimadas a partir das primárias, da parte rasa dos dados, pela convolução destas com o dado, seguido pelo filtro de mínimos-quadrados para melhor corrigir as múltiplas presentes nos dados. A aplicação do método é ilustrada em dados sintéticos e reais e, em ambos os casos, para efeito de comparação, aplicamos também a técnica 1D de atenuação de múltiplas. Computacionalmente o método 2D de atenuação de múltiplas é mais pesado pois, as múltiplas são estimadas no domínio do parâmetro de raio descendente (pd,po,tau) e para isso é necessário o uso de duas transformadas tau-p diretas; uma ao longo do afastamento (xo --> po) e outra ao longo da fonte (xs --> ps), e a partir desse domínio, um mapeamento (ps --> pd) e, no final do processamento, uma tansformada tau-p inversa para retornar ao domínio inical. O método 1D trabalha com os dados no domínio (xs,po,tau), fazendo-se necessário somente uma transformada tau-p direta (xo --> po) tornando-se com isso computacionalmente mais eficiente. Os resultados obtidos demonstram o ótimo desempenho do presente método de atenuação principalmente quando comparado ao método 1D, especialmente quando o modelo apresenta camada com mergulho, onde a periodicidade entre as múltiplas e as primárias é violada em cada um dos diversos parâmetros de raio. Nesse caso, o método 1D de atenuação não consegue estimar corretamente as múltiplas e como consequência afeta a remoção destas do dado original. Nos dados reais ambos os métodos mostraram-se eficientes conseguindo o método 2D de atenuação de múltiplas melhores resultados principalmente no dado do Golfo do México que possui uma superfície irregular e várias múltiplas associadas à almofada de sal.

In this work we discuss and apply a wave based multiple attenuation method in plane wave domain. In the tau-p domain the multiples are periodic related to the intercerpt time only for horizontally layered medium and can be estimate and attenuated in the same ray parameter trace. For 2D media this periodicity is not verified then a bidimension method is used and the lateral variation of the structure is handled by coupling different upgoing and downgoing plane waves. The method belong to the prediction and subtraction approach. The predicted multiples are estimated through the convolution of the shallow primaries with the data and the corrected multiples are made by least-square filtering. We illustrate the aplication of this method for both sinthetic and real seismic data. Such results are compared with 1D atenuation multiples method. The estimation of multiples in 2D method are made on (pd,po,tau) domain, wich requires two direct tau-p transforms, one mapping of pd from ps and one inverse tau-p transform. These transforms and calculations are computationaly expensive compared to 1D method. In the latter method only one tau-p transform is necessary, leading to (xs,po,tau) domain. The final results show the optimum performance of the 2D atenuation method, compared to 1D method. For 2D medium, the 1D atenuation method fails because the multiples are not correctly estimated. In real data both methods are efficient, however the results of 2D atenuation method are more sensitive to lateral variations, like salt domes, and irregular surfaces.