O objetivo principal deste trabalho é o desenvolvimento de modelos computacionais para analisar a resposta eletromagnética de ferramentas de perfilagem LWD/MWD em formações geofísicas arbitrárias. Essa modelagem envolve a determinação precisa de campos eletromagnéticos em regiões tridimensionais (3D) complexas e, conseqüentemente, a solução de sistemas lineares não-hermitianos de larga escala. A modelagem numérica é realizada através da aplicação do método dos volumes finitos (FVM) no domínio da freqüência. Desenvolvem-se dois modelos computacionais, o primeiro válido em regiões isotrópicas e o segundo considerando a presença de anisotropias no meio. As equações de Maxwell são resolvidas através de duas formulações distintas: formulação por campos e formulação por potenciais vetor e escalar. A discretização por volumes finitos utiliza um esquema de grades entrelaçadas em coordenadas cilíndricas para evitar erros de aproximação de escada da geometria da ferramenta. Os modelos desenvolvidos incorporam quatro técnicas numéricas para aumentar a eficiência computacional e a precisão do método. As formulações por campos e por potenciais vetor e escalar são comparadas em termos da taxa de convergência e do tempo de processamento em cenários tridimensionais. Os modelos foram validados e testados em cenários tridimensionais complexos, tais como: (i) poços horizontais ou direcionais; (ii) formações não homogêneas com invasões de fluído de perfuração; (iii) formações anisotrópicas e (iv) poços excêntricos. Motivado pela flexibilidade dos modelos e pelos resultados numéricos obtidos em diferentes cenários tridimensionais, estende-se a metodologia para analisar a resposta de ferramentas LWD que empregam antenas inclinadas em relação ao eixo da ferramenta. Tais ferramentas podem prover dados com sensibilidade azimutal, assim como estimativas da anisotropia da formação, auxiliando o geodirecionamento de poços direcionais e horizontais.