Modelagem e migração em profundidade 2D em meios com simetria polar local

Abstract

ABSTRACT. This paper shows a technique based on the phase-shift method (PSM) to implement pre-stack depth migration on locally transverse isotropic media (LTI), with the symmetry axis direction varies continually along the layers. For seismic modeling, a generalization of the finite differences method for the solution of the elastic wave equation was used. With this procedure, it was possible to accommodate seismic modeling on LTI media defined by six parameters at each grid point, i.e. , density, P and S wave propagation velocities along the local symmetry axis, Thomsen parameters and the direction of the local symmetry axis itself. In order to separate from the seismograms the qP and qSV wavefields, an algorithm based on the Christoffel equation was implemented. The migration for each common shot gather is implemented solely by phase-shift based algorithms, which means that not only the depropagation of the registered wavefield, but also the generation of the time matrices involved in the imaging condition were obtained in this manner for each set of parameters at each depth level. The migration results using qP-qP and qP-qSV reflections show that the horizons were located precisely, and that the process is stable in relation to the symmetry axis variations. The proposed method is for multicomponent seismic acquisitions and might be applied to marine seismic data using streamers, or Ocean Bottom Cables or vertical cables. Since the proposed method uses phase-shift algorithms, its parallel implementation can be highly efficient.Keywords: migration, modeling, phase-shift, anisotropy, locally transverse isotropic media, finite differences. RESUMO.Este trabalho propõe uma técnica do tipo rotação de fase (PSM -Phase-Shift Method), para migração em profundidade de dados sísmicos para meios localmente transversalmente isotrópicos (LTI), nos quais a direção do eixo de simetria varia continuamente ao longo das camadas. Para a modelagem em meios LTI, onde cada ponto da malha tem suas características definidas pelas velocidades de fase P e SV, parâmetros de Thomsen, densidade e inclinação do eixo de simetria, desenvolveu-se um método baseado na solução da equação elástica da onda por diferenças finitas. Na separação dos modos de onda qP e qSV dos sismogramas, implementou-se um algoritmo baseado na solução da equação de Christoffel. A migração para cada família de tiro comumé realizada somente por PSM. Nos resultados das migrações usando reflexões dos tipos qP-qP e qP-qSV, os horizontes foram localizados precisamente e verificou-se que o processoé estável em relaçãoà variação do eixo de simetria. O método propostoé para aquisições sísmicas multicomponentes podendo ser aplicado em dados sísmicos marítimos convencionais, como também em dados provenientes de aquisições do tipo cabo de fundo e cabo vertical. Como o método proposto se baseia em algoritmos que utilizam PSM a sua implementação paralela pode ser altamente eficiente.Palavras-chave: migração, modelagem, rotação de fase, anisotropia, meios localmente transversalment...