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.