Modelagem sísmica e inversão na presença de anisotropia

A necessidade da adoção de modelos elásticos anisotrópicos, no contexto da sísmica de exploração, vem crescendo com o advento de novas técnicas de aquisição de dados como VSP, walkway VSP, tomografia poço a poço e levantamentos sísmicos com grande afastamento. Meios anisotrópicos, no contexto da sísmica de exploração, são modelos efetivos para explicar a propagação de ondas através de meios que apresentam padrões de heterogeneidade em escala muito menor que o comprimento de onda das ondas sísmicas. Particularmente, estes modelos são muito úteis para explicar o dado sísmico mais robusto que são as medidas de tempo de trânsito. Neste trabalho, são investigados aspectos da propagação de ondas, traçado de raios e inversão de tempos de trânsito em meios anisotrópicos. É estudada a propagação de ondas SH em meios anisotrópicos estratificados na situação mais geral onde estas ondas podem ocorrer, ou seja, em meios monoclínicos com um plano vertical de simetria especular. É mostrado que o campo de ondas SH refletido a partir de um semi-espaço estratificado, não apresenta qualquer informação sobre a possível presença de anisotropia em subsuperfície. São apresentados métodos simples e eficientes para o traçado de raios em 3D através de meios anisotrópicos estratificados, baseados no princípio de Fermat. Estes métodos constituem o primeiro passo para o desenvolvimento de algoritmos de inversão de tempos de trânsito para meios anisotrópicos em 3D, a partir de dados de VSP e walkaway VSP. Esta abordagem é promissora para determinação de modelos de velocidade, que são necessários para migração de dados sísmicos 3D na presença de anisotropia. É efetuada a análise da inversão tomográfica não linear, para meios estratificados transversalmente isotrópicos com um eixo de simetria vertical(TIV). As limitações dos dados de tempo de trânsito de eventos qP para determinação das constantes elásticas, são estabelecidas e caracterizados os efeitos da falta de cobertura angular completa na inversão tomográfica. Um algoritmo de inversão foi desenvolvido e avaliado em dados sintéticos. A aplicação do algoritmo a dados reais demonstra a consistência de meios TIV. Esta abordagem é útil para casos onde há informação a priori sobre a estratificação quase plana das formações e onde os próprios dados do levantamento poço a poço apresentam um alto grau de simetria especular em relação a um plano vertical. Também pode ser útil em interpretações preliminares, onde a estimativa de um meio estratificado, serve como modelo de fundo para se efetuar análises mais detalhadas, por exemplo, como um modelo de velocidades anisotrópico para migração, ou como um modelo de calibração para análises de AVO.

Análise das aproximações RPP e RSP para meios isotrópicos

Neste trabalho são apresentadas aproximações lineares e quadráticas das equações de Zoeppritz para a obtenção dos coeficientes de reflexão e transmissão de eventos P-P e P-S em função dos ângulos de incidência e da média angular, bem como a análise de inversão linear de AVO, considerando os eventos de reflexão P-P e P-S dissociados e combinados. O uso das chamadas aproximações pseudo-quadráticas foi aplicado para obtenção de aproximações quadráticas apenas para eventos-PP, em torno dos contrastes médios das velocidades de ondas compressionais e cisalhantes e da razão Vs/Vp. Os resultados das aproximações desenvolvidas neste trabalho mostram que as aproximações quadráticas são mais precisas que as lineares, nas duas versões angulares. As comparações entre as aproximações em termos do ângulo de incidência e da média angular mostram que as aproximações quadráticas são equivalentes dentro do limite angular de [0º, 30º]. Por outro lado, as aproximações lineares em função do ângulo de incidência mostram-se mais precisas que as aproximações lineares em função da média angular. Na inversão linear, fez-se análises de sensibilidade e de ambigüidade e observou-se que, nos caso de eventos de reflexão P-P e P-S dissociados, apenas um parâmetro pode ser estimado e que a combinação destes eventos consegue estabilizar a inversão, permitindo a estimativa de dois dentre os parâmetros físicos dos meios (contrastes de impedância, de velocidade de onda P e de módulo de cisalhamento).

Análise do efeito da discretização do modelo de velocidades nas migrações Kirchhoff e Kirchhoff-Gaussian- Beam 2D pré-empilhamento em profundidade

O Feixe Gaussiano (FG) é uma solução assintótica da equação da elastodinâmica na vizinhança paraxial de um raio central, a qual se aproxima melhor do campo de ondas do que a aproximação de ordem zero da Teoria do Raio. A regularidade do FG na descrição do campo de ondas, assim como a sua elevada precisão em algumas regiões singulares do meio de propagação, proporciona uma forte alternativa no imageamento sísmicos. Nesta dissertação, apresenta-se um novo procedimento de migração sísmica pré-empilhamento em profundidade com amplitudes verdadeiras, que combina a flexibilidade da migração tipo Kirchhoff e a robustez da migração baseada na utilização de Feixes Gaussianos para a representação do campo de ondas. O algoritmo de migração proposto é constituído por dois processos de empilhamento: o primeiro é o empilhamento de feixes (“beam stack”) aplicado a subconjuntos de dados sísmicos multiplicados por uma função peso definida de modo que o operador de empilhamento tenha a mesma forma da integral de superposição de Feixes Gaussianos; o segundo empilhamento corresponde à migração Kirchhoff tendo como entrada os dados resultantes do primeiro empilhamento. Pelo exposto justifica-se a denominação migração Kirchhoff-Gaussian-Beam (KGB).Afim de comparar os métodos Kirchhoff e KGB com respeito à sensibilidade em relação ao comprimento da discretização, aplicamos no conjunto de dados conhecido como Marmousi 2-D quatro grids de velocidade, ou seja, 60m, 80m 100m e 150m. Como resultado, temos que ambos os métodos apresentam uma imagem muito melhor para o menor intervalo de discretização da malha de velocidade. O espectro de amplitude das seções migradas nos fornece o conteúdo de frequência espacial das seções das imagens obtidas.

Análise de parâmetros da função característica de Hamilton

A análise de velocidades é um processo fundamental na sísmica de reflexão, onde as velocidades de empilhamento bem como o tempo de trânsito de afastamento nulo (onde supõe-se que a fonte e o detetor ocupam a mesma posição) são parâmetros suficientes na determinação do modelo geológico para meios com camadas horizontais. Por outro lado quando expresso através do formalismo estabelecido por Hamilton, os mesmos parâmetros (velocidade e tempo de afastamento nulo) são suficientes para determinar a função característica para este mesmo tipo de meio. Para o caso de um modelo geológico heterogêneo com interfaces arbitrariamente curvas, a função característica de Hamilton é dada a partir da estimativa de nove parâmetros, onde os mesmos parâmetros são necessários na determinação do modelo geológico em 3D. Este trabalho tem por objetivo estimar os parâmetros que determinam a função característica de Hamilton para meios 3D e estudar a influência de cada parâmetro na função, através de cortes horizontais nas seções de tempos de trânsito (conhecidos como time slices), nas configurações de ponto médio comum e afastamento nulo. Dentro desta abordagem é dado um exemplo a partir de um modelo sintético onde, aqueles resultados obtidos com estudo da influência de cada parâmetro na função característica, são aplicados como um critério de ajuste entre a função característica calculada e a função de tempos de trânsito obtida no levantamento de dados.

Análise de resolução em tomografia interpoços em 3-D

A aplicação de métodos de inversão aos dados de múltiplos levantamentos sísmicos interpoços permite a reconstrução de modelos de vagarosidade em 3-D de alta resolução adequados para monitoramento de processos de recuperação avançada de petróleo e caracterização de reservatórios. Entretanto, a falta de cobertura volumétrica uniforme de raios de levantamentos interpoços exige informação adicional ao sistema tomográfico para obtenção de soluções estáveis. A discretização do modelo em uma malha 3-D com células prismáticas triangulares e a decomposição em valores singulares são utilizadas para avaliar a reconstrução tomográfica em 3-D. O ângulo da projeção de modelos-alvo no subespaço ortogonal ao espaço nulo efetivo da matriz tomográfica é um critério adequado para se otimizar a malha de discretização do modelo interpretativo e a geometria de aquisição dos dados de modo a melhorar o condicionamento da reconstrução tomográfica. Esta abordagem pode ser utilizada durante as iterações lineares para redefinir a malha ou avaliar a necessidade de informação a priori adicional ao sistema tomográfico.

Análise de sensibilidade para estereotomografia em meios elípticos e anelípticos

O método de estereotomografia é estendido para meios com anisotropia arbitrária e implementado para meios com anisotropia elíptica e anelíptica. Os modelos elípticos e anelípticos apresentam somente três parâmetros. Isto faz com que eles sejam menos sensíveis a ambiguidade, causada pela cobertura limitada dos raios em experimentos sísmicos de superfície e VSP, do que modelos transversalmente isotrópicos ou ortorrômbicos. As correspondentes aproximações para superfície de vagarosidade limita a validade desta implementação para eventos qP com anisotropia suave. Experimentos numéricos mostram o potencial e as limitações da estereotomografia para estimar macro modelos de velocidade adequados para o imageamento na presença de anisotropia e a importância dos eventos de transmissão de experimentos VSP de multiplo afastamento modelo para o sucesso desta abordagem.

Inversão de dados de sísmica de refração profunda a partir da curva tempo-distância

O trabalho em pauta tem como objetivo o modelamento da crosta, através da inversão de dados de refração sísmica profunda, segundo camadas planas horizontais lateralmente homogêneas, sobre um semi-espaço. O modelo direto é dado pela expressão analítica da curva tempo-distância como uma função que depende da distância fonte-estação e do vetor de parâmetros velocidades e espessuras de cada camada, calculado segundo as trajetórias do raio sísmico, regidas pela Lei de Snell. O cálculo dos tempos de chegada por este procedimento, exige a utilização de um modelo cujas velocidades sejam crescentes com a profundidade, de modo que a ocorrência das camadas de baixa velocidade (CBV) é contornada pela reparametrização do modelo, levando-se em conta o fato de que o topo da CBV funciona apenas como um refletor do raio sísmico, e não como refrator. A metodologia de inversão utilizada tem em vista não só a determinação das soluções possíveis, mas também a realização de uma análise sobre as causas responsáveis pela ambiguidade do problema. A região de pesquisa das prováveis soluções é vinculada segundo limites superiores e inferiores para cada parâmetro procurado, e pelo estabelecimento de limites superiores para os valores de distâncias críticas, calculadas a partir do vetor de parâmetros. O processo de inversão é feito utilizando-se uma técnica de otimização do ajuste de curvas através da busca direta no espaço dos parâmetros, denominado COMPLEX. Esta técnica apresenta a vantagem de poder ser utilizada com qualquer função objeto, e ser bastante prática na obtenção de múltiplas soluções do problema. Devido a curva tempo-distância corresponder ao caso de uma multi-função, o algoritmo foi adaptado de modo a minimizar simultaneamente várias funções objetos, com vínculos nos parâmetros. A inversão é feita de modo a se obter um conjunto de soluções representativas do universo existente. Por sua vez, a análise da ambiguidade é realizada pela análise fatorial modo-Q, através da qual é possível se caracterizar as propriedades comuns existentes no elenco das soluções analisadas. Os testes com dados sintéticos e reais foram feitos tendo como aproximação inicial ao processo de inversão, os valores de velocidades e espessuras calculados diretamente da interpretação visual do sismograma. Para a realização dos primeiros, utilizou-se sismogramas calculados pelo método da refletividade, segundo diferentes modelos. Por sua vez, os testes com dados reais foram realizados utilizando-se dados extraídos de um dos sismogramas coletados pelo projeto Lithospheric Seismic Profile in Britain (LISPB), na região norte da Grã-Bretanha. Em todos os testes foi verificado que a geometria do modelo possui um maior peso na ambiguidade do problema, enquanto os parâmetros físicos apresentam apenas suaves variações, no conjunto das soluções obtidas.

Mapeamento de jazidas de bauxitas no Distrito Paragominas – Açailândia (PA-MA) pelo método da sísmica de refração

Nos últimos 30 anos no Brasil, o mapeamento de camada bauxítica tem sido sempre efetuado através de perfilagens de poços. Neste estudo, aplicou-se o método da sísmica de refração rasa pela primeira vez, em cerca de 3.000m de perfis, em caráter experimental em algumas localidades ao longo da rodovia BR-10, assim como no Platô Esperança (região de Jabuti-PA) e perto da cidade de Açailândia (MA); com objetivo de delimitar a camada de bauxita. Embora se tenha utilizado, em algumas ocasiões, durante a interpretação, o método de distância crítica no cálculo de profundidade, foi possível, mediante o emprego do Método Recíproco Generalizado (GRM), o processamento e interpretação dos dados, para se obter a delimitação contínua da camada investigada. O GRM apresenta vantagens sobre os métodos tradicionais, tais como poder delimitar camadas onduladas, detetar camadas escondidas, além da sensibilidade para camadas com mergulho até 20°. A partir dós gráficos da análise de velocidade, foi calculada para a camada de bauxita uma velocidade média de 840 m/s, e mediante o emprego dos tempos de profundidade chegou-se a sua delimitação: ela possui um comportamento irregular, tendo uma espessura média de 5m. A pesquisa geofísica realizada nas três áreas-piloto tem grande importância porque o método utilizado poderá ser aplicado a outras regiões, em condições geológicas semelhantes.

Modelagem sísmica por diferenças finitas em meios bidimensionais com difratores

Este trabalho tem por objetivo a modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações. Como parte deste estudo foi feita a análise numérica visando as condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas. Para a validação do programa de diferenças finitas foi feita a comparação cinemática com a Teoria do Raio para um modelo simples. O motivo deste estudo é ter uma melhor compreensão e controle sobre os problemas de modelagem, visando contribuir para a solução de problemas na interpretação de dados sísmicos. Segundo vários autores na literatura geológica, Derby (1877), Evans (1906), Paiva (1929) e Moura (1938). A Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o recente, atingindo espessuras da ordem de 5 km. Os corpos de diabásio, presentes entre os sedimentos paleozóicos, estão dispostos na forma de soleiras, alcançando espessuras de centenas de metros, perfazendo um volume total de 90.000 quilômetros cúbicos. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica, impossibilitando uma melhor interpretação dos horizontes refletores que se encontram abaixo das soleiras. Para representar situações geológicas desse tipo foi usado um modelo acústico de velocidades. Para o cálculo dos sismogramas foi utilizado um programa de diferenças finitas com aproximação de quarta ordem da equação da onda acústica no espaço e no tempo. As aplicações dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas têm melhorado a compreensão sobre a propagação em meios onde existem heterogeneidades significativas, tendo como resultado boa resolução na interpretação dos eventos de reflexão sísmica em áreas de interesse. Como resultado dos experimentos numéricos realizados em meio de geologia complexa, foi observado a influência significativa das múltiplas devido a camada de alta velocidade, o que faz com que haja maior perda de energia dificultando a interpretação dos alvos. Por esta razão recomendo a integração de dados de superfície com os dados de poço, com o objetivo de se obter uma melhor imagem dos alvos abaixo das soleiras de diabásio.

Caracterização dos reservatórios do grupo Macaé no Campo de Garoupa - Bacia de Campos - por meio de interpretação sísmica e perfis de poço

The study area is the Garoupa Field, located 75 km from the coast of Rio de Janeiro in bathymetric elevation of 125 m and it is the first important offshore discovery of the Campos Basin with hydrocarbon production in the Albian carbonates. The high significance of carbonate rock reservoirs for the global oil industry and its degree of complexity, combined with the fact that the Garoupa Field is a mature field, make this research project interesting to develop a re-exploratory vision of the study area. The objective of this paper is to characterize the physical properties of the Macae Group reservoirs in the Garoupa Field through well correlation techniques, petrophysical analysis and seismic interpretation (3D), in order to evaluate qualitatively and quantitatively the response of geophysical data to the presence of carbonate rock reservoirs. From the analysis of well profiles (gamma ray, resistivity, density and sonic) it was possible to identify six reservoir levels in the Quissamã Formation. As a result, the reservoirs have good correlation between the wells in the Garoupa Field and show lateral facies variation. In the seismic interpretation, structural contour maps were generated in time (ms) from the top of the Macae Group and two horizons that correspond to the reservoir levels. These maps show that the oil producing area of this field is a structural high located at the west of the interpreted seismic cube, and it was classified as a rollover limited by lateral faults forming a horst. The seismic attribute maps show a good response to the reservoirs of the Garoupa Field because they are associated with amplitude anomalies, verified by the correlation between the physical properties of the reservoir rock and seismic data through cross plots, with emphasis in a linear correlation between the resistivity profile and the values of Maximum Absolute Amplitude and RMS Amplitude attributes

Sismicidade e esforços tectônicos na zona sísmica Acaraú, Nordeste do Brasil

Inside the Borborema Province the Northwestern Ceará (NC) is one of the most seismic active regions. There are reports of an earthquake occurred in 1810 in the Granja town. On January, 2008 the seismic activity in NC has increased and it was deployed a seismographic network with 11 digital stations. In 2009, another earthquake sequence began and it was deployed another seismographic network in the Santana do Acaraú town with 6 stations. This thesis presents the results obtained by analyzing the data recorded in these two networks. The epicentral areas are located near the northeastern part of the Transbrasiliano Lineament, a shear zone with NE-SW-trending that cuts the study area. The hypocenters are located between 1km and 8km. The strike-slip focal mechanisms were found, which is predominant in the Borborema Province. An integration of seismological, geological and geophysical data was performed and it show that the seismogenic faults found are oriented in the same direction to the local brittle structures observed in field and magnetic lineaments. The SHmax (maximum compressional stress) direction in NC was estimated using an inversion of seven focal mechanisms. The horizontal maximum compression stress (σ1 = 300°) with orientation NW-SE and extension (σ3 = 210°) with NE-SW and σ2 vertical. These results are consistent with results of previous studies. The seismic activity recorded in NC is not related to a possible reactivation of the Transbrasiliano Lineament, by now.

Análise estratigráfica de subsuperfície do Devoniano Inferior da Bacia do Rio do Peixe, Nordeste do Brasil

Subsurface stratigraphic analysis of Devonian strata from the Rio do Peixe Basin, newly recognized by palynological studies, has resulted in the identification of two new lithostratigraphic units assembled in the Santa Helena Group. The Pilões Formation, the lower unit, is composed mainly of dark mudstones and medium-tovery fine-grained sandstones, with minor conglomerates and breccias. The Triunfo Formation, the upper unit, comprises whitish grey, kaolinitic, coarse-grained to conglomeratic, cross stratified sandstones and conglomerates, with interbedded mudstones and fine-grained sandstones. These units were characterized using cores, sidewall and cuttings samples, conventional logs and image log, from three wells drilled by PETROBRAS, and 3D seismic data. The Pilões Formation is interpreted as prodeltaic facies, with lesser associated subaqueous talus, debrite and sandy turbidite lobe facies, distal part of fandelta and braided fluviodeltaic facies of Triunfo Formation. The Santa Helena Group corresponds to the Lower Devonian tectono-sequence deposited in a NW-SE-trending graben during a transgressiveregressive cycle. With 343 meters of thickness (isochore) in well 1-PIL-1-PB (Pilões 1), this sequence has a non-conformity at the lower boundary and its upper boundary is an unconformity with the Lower Cretaceous tectono-sequence (Rio do Peixe Group), that represents a hiatus of about 265 million years. Ignimbrites and coignimbrite breccias (Poço da Jurema volcanic breccia), related to an unknown pyroclastic volcanic event, were recognized at the northern margin of the Sousa halfgraben. Evidence from well data suggests that this event is coeval with the Devonian graben filling. The present study indicates a polyhistorical tectono-volcanosedimentary evolution of the basin. This lithostratigraphic update brings new perspectives for geological research in the Rio do Peixe Basin, as well as in other inland basins of the Northeastern of Brazil. The results of the research also contribute to the kwnoledge of the Borborema Province and western Gondwana paleogeography during the Early Devonian.

Análise estratigráfica da sequência mesodevoniana-eocarbonífera da Bacia do Parnaíba, nordeste do Brasil

The Middle Devonian-Early Carboniferous sequence of the Parnaíba Basin, lithostratigraphically defined as Canindé Group, has been reinterpreted using the basic model of sequence stratigraphy. Therefore, lithology and gamma ray well-logs and seismic lines of central portion of the basin were analyzed, producing up from there diagrams 1D, isochore maps and stratigraphic sections. As results of this study, were defined two depositional cycles of second order, referred as Depositional Sequence 1 (SEQ1) and the Depositional Sequence 2 (SEQ2). The SEQ1, with interval about 37 Ma, is limited below by Early Devonian Unconformity and is equivalent to the formations Itaim, Pimenteiras and Cabeças. The SEQ2, which follows, comprises a range of about 15 Ma and is equivalent to the Longá Formation The SEQ1 starts with the lowstand systems tract, consisting of progradational parasequence set in the basal part, predominantly pelitic, deposited on a prodelta under influence of storms and the upper part consists in sandstones of deltaic front, with the maximum regressive surface on the upper limit. The transgressive systems tract, deposited above, is characterized by retrogradacional parasequence set composed of shallow shelf mudstones, deposited under storm conditions. The maximum flooding surface, upper limit of this tract, is positioned in a shale level whose radioactivity in gammaray well-log is close to 150 API. The highstand systems tract presents progradational parasequence set, comprising mudstones and sandstones deposited in shelf, fluvial-estuarine or deltaic and periglacial environments, with the upper limit the Late Devonian Unconformity. The SEQ2 was deposited in shelf environment, starting with the lowstand systems tract, that is characterized by a progradational parasequence set, followed by the transgressive systems tract, with retrogradational character. The upper limit of the tract corresponding to the fusion between maximum flooding surface with the upper limit of this sequence, which is the Early Carboniferous Unconformity, where the overlapping section was eroded. This section, which corresponds the highstand systems tract is restricted to portions at which the erosion that generate the Early-Carboniferous Unconformity was less effective, preserving the records of this unit.

Análise estratigráfica da sequência mesodevoniana-eocarbonífera da Bacia do Parnaíba, nordeste do Brasil

The Middle Devonian-Early Carboniferous sequence of the Parnaíba Basin, lithostratigraphically defined as Canindé Group, has been reinterpreted using the basic model of sequence stratigraphy. Therefore, lithology and gamma ray well-logs and seismic lines of central portion of the basin were analyzed, producing up from there diagrams 1D, isochore maps and stratigraphic sections. As results of this study, were defined two depositional cycles of second order, referred as Depositional Sequence 1 (SEQ1) and the Depositional Sequence 2 (SEQ2). The SEQ1, with interval about 37 Ma, is limited below by Early Devonian Unconformity and is equivalent to the formations Itaim, Pimenteiras and Cabeças. The SEQ2, which follows, comprises a range of about 15 Ma and is equivalent to the Longá Formation The SEQ1 starts with the lowstand systems tract, consisting of progradational parasequence set in the basal part, predominantly pelitic, deposited on a prodelta under influence of storms and the upper part consists in sandstones of deltaic front, with the maximum regressive surface on the upper limit. The transgressive systems tract, deposited above, is characterized by retrogradacional parasequence set composed of shallow shelf mudstones, deposited under storm conditions. The maximum flooding surface, upper limit of this tract, is positioned in a shale level whose radioactivity in gammaray well-log is close to 150 API. The highstand systems tract presents progradational parasequence set, comprising mudstones and sandstones deposited in shelf, fluvial-estuarine or deltaic and periglacial environments, with the upper limit the Late Devonian Unconformity. The SEQ2 was deposited in shelf environment, starting with the lowstand systems tract, that is characterized by a progradational parasequence set, followed by the transgressive systems tract, with retrogradational character. The upper limit of the tract corresponding to the fusion between maximum flooding surface with the upper limit of this sequence, which is the Early Carboniferous Unconformity, where the overlapping section was eroded. This section, which corresponds the highstand systems tract is restricted to portions at which the erosion that generate the Early-Carboniferous Unconformity was less effective, preserving the records of this unit.

Análise da trajetória de poços horizontais em reservatório de óleo pesado

The application of thermal methods, to increase the recovery of heavy oil in mature fields through drainage with multilateral and horizontal wells, has been thoroughly studied, theorically, experimentally, testing new tools and methods. The continuous injection of steam, through a steam injector well and a horizontal producer well in order to improve horizontal sweep of the fluid reservoir, it is an efficient method. Starting from an heterogeneous model, geologically characterized, modeling geostatistics, set history and identification of the best path of permeability, with seismic 3D, has been dubbed a studying model. It was studied horizontal wells in various directions in relation to the steam and the channel of higher permeability, in eight different depths. Into in the same area were studied, the sensitivity of the trajectories of horizontal wells, according to the depth of navigation. With the purpose of obtaining the highest output of oil to a particular flow, quality, temperature and time for the injection of steam. The wells studied showed a significant improvement in the cumulative oil recovery in one of the paths by promoting an alternative to application in mature fields or under development fields with heavy oil