Tradicionalmente, para visualização de campos vetoriais em modelos discretos de reservatórios naturais de petróleo, traça-se linhas de fluxo resolvendo a sua equação diferencial ordinária célula-a-célula, seja através de soluções analíticas ou numéricas, considerando o campo de velocidade local de cada célula. Essa estratégia tem como desvantagem traçar a linha considerando um campo de velocidade discreto e portanto descontínuo. Além disso, para modelos massivos, resolver a equação célula-a-célula pode tornar o método ineficiente. Neste trabalho, exploramos uma estratégia diferente: ao traçar as linhas de fluxo considera-se um campo de velocidade contínuo, representado pelo modelo discreto do reservatório. Para tanto, propõe-se: (i) o uso de uma estrutura espacial para acelerar a localização de um ponto no modelo de reservatório; (ii) o uso de interpolação esférica para avaliação do campo de velocidade a partir do modelo discreto; (iii) o uso de um método numérico adaptativo para controlar o erro numérico da integração. Os resultados obtidos em modelos de reservatórios reais demonstram que o método proposto atende aos requisitos de precisão, mantendo um bom desempenho.