Método híbrido de migração pré-empilhamento em profundidade no domínio da freqüência em duas etapas com interpolação

Abstract

ABSTRACT. This paper shows the application of Phase Shift Plus Interpolation Split-Step (PSPI-SS) wavefield extrapolation operator coupled with the excitationtime imaging condition to perform prestack seismic depth migration of common shot profiles. The migration methods which are characterized by the use of different propagators for the source and receiver wavefields are called hybrid migration methods. Normally, the source wavefield is determined by the computation of the traveltimes from the solution of the eikonal equation. Here the traveltimes were calculated based on a finite difference solution of the eikonal equation (TGRID), a simple geometrical approach (FAST) and also by the paraxial ray tracing method (PARAX). For the extrapolation of the receiver data we applied the PSPI Split-Step operator, that is derived from the conventional Split-Step wave equation method using multiple references velocity. As the source traveltime computation is based on the ray assintotic theory the velocity field needs to be smoothed and the migration results obtained with the Kirchhof migration method, that is a very fast migration procedure, was used to choose the best level of smoothing of the velocity field. In this work the synthetic Marmousi data and a real dataset from Gulf of Mexico were used to test the migration methods proposed here. The migration results with the hybrid methods shown better images in quality then the kirchhoff's results and comparable to those obtained with the full PSPI-SS migration method, but a low computation cost.Keywords: Hybrid migration, Prestack migration, Split-Step Phase Shift Plus Interpolation. RESUMO.Neste trabalho apresentamos a aplicação do operador de extrapolação em profundidade em duas etapas com interpolação (PSPI-SS), operador que faz uso da correção Split-Step em cada passo de extrapolação, conjuntamente com a condição de imagem de tempo de excitação, para realizar a migração pré-empilhamento em profundidade de dados sísmicos ordenados em famílias de tiro comum. Os métodos de migração pré-empilhamento que realizam a extrapolação dos dados ao longo dos receptores e a continuação para baixo das fontes sísmicas com operadores diferentes, são chamados de métodos híbridos de migração. Nesses métodos o campo de onda de cada fonte sísmicaé determinado através do cálculo dos tempos de chegada da fonte a partir da aproximação de alta freqüência. Aqui utilizamos apenas os tempos das primeiras chegadas, que são calculados a partir da solução da equação iconal por diferença finitas (TGRID), por um método puramente geométrico (FAST) e também pelo método de traçamento de raios paraxiais (PARAX). Na implementação proposta, o campo de velocidades usado para o cálculo dos tempos de chegada precisa ser suavizado e o seu grau de suavizaçãoé determinado a partir dos resultados apresentados pela migração convencional Kirchhoff. A extrapolação dos dados dos receptoresé feita usando-se o campo de velocidades original sem nenhuma suavização, através do operador PSPI-SS. Para testar os novos métodos...