REVERSE TIME MIGRATION IN THE FREQUENCY DOMAIN BY THE RAPID EXPANSION METHOD

ABSTRACT. This paper proposes and describes the implementation of a new depth migration method in the frequency domain. The method, based in the reverse time migration (RTM) technique, extrapolates wavefields from the source and receivers to obtain migrated seismic images that are built directly into the frequency domain. In the proposed method, wavefields are propagated in the time domain and are then transformed into the frequency domain at each time extrapolation step through the discrete Fourier transform. Neither the forward nor backward wavefield is needed to be stored in memory or read from disk storage. To speed up the migration algorithm, the discrete Fourier transform kernel for each frequency is computed and salved before the time extrapolation procedure. At the imaging condition phase, both source and receiver wavefields are at the same frequency, so that, the construction of the image occurs by multiplying the forward source propagated wavefield with the backward propagated of the receivers wavefield for each frequency component. Subsequently, saving the source field at each step to later correlate it with the backpropagated receiver wave field, usually done in conventional RTM, becomes unnecessary. Nor is it necessary to invert a matrix for each frequency component, which is done in the migration technique that uses the Helmholtz equation solution in the frequency domain. Thus, the migration procedure in the frequency domain being proposed is more efficient from a computational point of view, and can also produce high quality migrated images as those produced by conventional RTM. The rapid expansion method (REM) is used for seismic forward modeling, which extrapolated data with good precision and free of numerical dispersion. Thus, with the transformed data at each step in the frequency domain, it is possible to construct high quality, in-depth seismic images at a lower computational cost. Moreover, this frequency domain migration with REM is an atractive strategy to design robust inverse algorithms, especially for 3D problems. To demonstrate the efficiency and applicability of the proposed method, two synthetic models were used and their results showed high quality images equivalent to those obtained by conventional RTM and thus proving the vality of the method. Keywords: wave equation migration, depth migration, imaging condition, frequency domain migration. RESUMO. Um metodo de migracao em profundidade no dominio da frequencia e proposto e implementado. O metodo consiste na extrapolacao dos campos de ondas da fonte e dos receptores e baseia-se na tecnica de migracao reversa no tempo (da sigla em ingles, RTM), obtendo imagens sismicas migradas, construidas diretamente no dominio da frequencia. No metodo que estamos propondo, os campos de ondas sao propagados no dominio do tempo e a cada passo de extrapolacao sao transformados para o dominio da frequencia, atraves da transformada de Fourier discreta (do ingles, on-the-fly transform). Para acelerar o algoritmo de migracao, o kernel da transformada de Fourier e calculado fora do loop do tempo. Alem disso, na etapa de condicao da imagem, os campos de onda, tanto da fonte como dos receptores, sao calculados no mesmo instante de tempo, ou seja, a construcao da imagem se da atraves da multiplicacao do campo de onda da fonte com o campo retropropagado dos receptores, para cada componente de frequencia. Portanto, nao precisamos salvar o campo da fonte a cada passo no tempo para posteriormente correlacionar com o campo de onda retropropagado dos receptores, como e usualmente feito na RTM convencional, nem e preciso inverter uma matriz para cada componente de frequencia, como e realizado normalmente pela tecnica de migracao no dominio da frequencia, utilizando a solucao da equacao de Helmholtz. Desta forma, o procedimento de migracao no dominio da frequencia que estamos propondo se torna mais eficiente do ponto de vista computacional, podendo produzir imagens migradas de alta qualidade, quando comparadas as obtidas atraves da RTM convencional no dominio do tempo. Para a extrapolacao dos campos de ondas no tempo foi empregado o metodo de expansao rapida (da sigla em ingles, REM), que permite a extrapolacao dos dados com boa precisao e livres de dispersao numerica. Desta forma, com os dados transformados para o dominio da frequencia, a cada passo no tempo, e possivel a construcao de imagens sismicas em profundidade de boa qualidade e a um menor custo computacional. Para demonstrar a eficiencia e aplicabilidade do metodo proposto, dois modelos sinteticos foram usados e seus resultados apresentaram imagens de alta qualidade equivalentes as obtidas pela RTM convencional. Palavras-chave: equacao de migracao da onda, migracao, condicao de imagem, migracao no dominio da frequencia.