[pt] Este trabalho propõe uma técnica do tipo rotação de fase para migração em profundidade de dados sísmicos para meios com simetria polar local (localmente transversalmente isotrópicos, LTI), nos quais a direção do eixo de simetria varia continuamente ao longo das camadas. São testadas, através de simulações numéricas de levantamentos sísmicos, a precisão e a estabilidade do método, em função da variação do eixo de simetria. Para a realização das simulações, desenvolveu-se um método a partir da solução da equação elástica da onda usando-se a técnica das diferenças finitas, que possibilita 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. 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 para determinar os operadores de separação. A migração para cada família de tiro comum é realizada por soluções da equação da onda usando somente técnicas de rotações de fase. De fato, tanto a depropagação do campo registrado quanto a geração das matrizes de tempo utilizadas na condição de imageamento, são realizadas por soluções que envolvem rotações de fase para cada conjunto de parâmetros, em cada nível de profundidade. 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. Vale ressaltar que o método não está restrito a 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 OBC (Ocean Bottom Cable) e com cabo vertical. Como o método proposto se baseia em algoritmos que utilizam técnicas de rotação de fase, a sua implementação conta com o beneficio de ser altamente paralelizável. / [en] This work proposes a technique based on the phase-shift method to implement pre-stack depth migration on locally transverse isotropic media (LTI), in which the direction of the symmetry axis varies continually along the layers. Through numerical seismic data simulations the methods robustness and stability were tested in relation to the axis symmetry variations. 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 not restricted to multicomponent seismic acquisitions, but it can be applied to marine seismic data using streamers, or Ocean Bottom Cables or vertical cables. Since the proposed method uses phaseshift algorithms, its parallel implementation can be highly efficient.