Depth remigration by means of the image wave equation

Este trabalho aborda a questão de como resolver a equação da onda imagem para o problema de remigração na profundidade através de métodos numéricos. O objetivo deste problema é a reconstrução de uma imagem das camadas geológicas do subsolo a partir de uma imagem previamente migrada com um modelo de velocidade, geralmente, incorreto. Nosso principal objetivo neste trabalho é a investigação de possíveis métodos que possam resolver os problemas que surgiram ao usarmos esquemas explícitos do método de diferenças _nitas na solução da equação da onda imagem em trabalhos anteriores, como, por exemplo, a dispersão numérica. Para isso, estudamos aqui o método de volumes _nitos, assim como esquemas implícitos do método de diferenças _nitas. O método de volumes _nitos possui como característica principal propagar as médias das células da malha ao invés de simplesmente os dados pontuais como é feito no método de diferenças _nitas. As outras tentativas para solucionar o problema da dispersão foram dois tipos de implementação de esquemas implícitos do método de diferenças _nitas, isto é, implementações implícitas de esquemas convencionais avaliados em pontos da malha e um esquema avaliado nos centros das células. A qualidade dos algoritmos estudados foi testada numericamente. Estes testes numéricos mostram que o método de volumes _nitos não é adequado para resolver o problema da dispersão, uma vez que a média calculada a cada passo aumenta o estiramento do pulso. Além disso, as implementações implícitas dos esquemas convencionais mostram o mesmo comportamento de dispersão que as implementações explícitas. Unicamente o esquema centrado foi capaz de melhorar a dispersão numérica em comparação com as implementações anteriores,porém somente para dados contendo exclusivamente baixas freqüênciasThis work approaches the question of how to solve the image-wave equation for depth remigration by numerical methods. The objective is the reconstruction of an image of the geologic layers of the subsoil from a previously migrated image with a different velocity model. Our main objective in this work is the investigation of possible methods that can solve the problems that appeared when using explicit _nite-difference schemes for the solution of the image-wave equation in previous works, particularly numerical dispersion. For this purpose, we study the method of _nite volumes, as well as implicit _nite-difference schemes. The main characteristic of the _nite-volume method is to simply propagate the averages in the cells of the mesh instead of the discretized data themselves as it is done in the _nitedifference method. As another attempt to solve the problem of the dispersion, we study two types of implementation of implicit _nite-difference schemes, that is, implicit implementations of conventional schemes evaluated out the edge of the cell and a scheme evaluated in the center of the cell. The quality of the studied algoritms has been tested numerically. These numerical tests show that the method of _nite volumes is not adequate to solve the problem of dispersion, for the average calculated in each step additionally increases the pulse stretch. Moreover, the implicit implementations of the conventional schemes show the same dispersion behavior as the explicit implementations. Solely the centered scheme was capable to improve the numerical dispersion in comparison with the previous implementations, however only for data containin