Neste trabalho foi desenvolvido e testado um método de migração pré-empilhamento no domínio da freqüência para dados sísmicos ordenados em famílias de ponto de tiro comum. A migração é realizada simultaneamente para as fontes e receptores usando-se operadores do tipo "split-step". Os dados registrados nos receptores são depropagados em profundidade e a propagação da fonte é simulada também usando-se operadores "split-step". A imagem final em profundidade é obtida somando-se todas as freqüências de interesse, durante o processo de correlação dos campos propagados e depropagados, para cada nível de profundidade e pela soma de todos os tiros migrados. A técnica de migração "split-step" foi escolhida por ser de fácil implementação e também por ser bastante robusta em situações de forte contraste lateral de velocidade. Visando a diminuir o tempo computacional dessa técnica de migração, uma segunda alternativa foi implementada para o cálculo dos tempos das ondas diretas, através da solução por diferencias finitas da equação eiconal. Este segundo é referido como método híbrido.

Quanto ao custo computacional dos métodos implementados, o método híbrido foi em média 30% mais rápido do que o método "split-step" completo, mas, não apresentou resultados equivalentes quanto a qualidade de imagem.

Os algoritmos desenvolvidos para os métodos de migração foram testados em dados sintéticos com alta complexidade geológica, tais como os modelos EAGE/SEG e Marmousi. Também foram testados em dados reais de uma linha sísmica 2-D da Carolina do Norte. Os resultados obtidos neste trabalho com o método "split-step" completo foram bastante satisfatórios, especialmente com os dados Marmousi, que apresentou uma imagem bem melhor do que a conseguida com o método híbrido.

In this work it was developed and tested a prestack depth migration method in frequency domain for seismic data organized in common shot gathers. The migration is carried out simultaneously for the sources and receivers using split-step operators. The receiver data is depropagated in depth and the sources wavefields is propagated both using split-step operators. In this method the images is obtained by adding all the interest frequencies, during the process of correlation of the propagated and depropagated fields, for each depth level and also summing all the migrated shots. The split-step migration technique was chosen for being of easy implementation and also for being sufficiently robust in situations of strong lateral velocity contrast. In order to diminish the computational time of this migration technique a second alternative was used to calculate the direct wave traveltimes, through a finite difference solution of eiconal equation. We called this second method hybrid method.

Comparing the computational cost of the implemented methods, the hybrid method was in average 30% faster than the full split-step method, but on the other hand it did not present equivalent image quality.

The algorithms developed for the migration methods were tested in synthetic data with high geologic complexity, such as the models EAGE/SEG and Marmousi. Also they were tested in real data in a North Carolina 2D seismic line. The results with the full split-step method had been sufficiently satisfactory, especially with the Marmousi data, that presents a much better migrated section than the obtained with the hybrid method.