Experiencias en innovación y mejora de la docencia de la ingenieria sísmica.

Resumen tomado de la publicación

Análise estratigráfica do Campo de Namorado (Bacia de Campos) com base na interpretação sísmica tridimensional

O presente trabalho analisa os reservatórios turbidíticos do Campo de Namorado, Bacia de Campos – RJ, com a apresentação de um novo modelo evolutivo para o intervalo entre o Albiano superior e Cenomaniano, na área do referido campo. As ferramentas utilizadas neste estudo consistiram da interpretação sísmica em ambiente tridimensional com o software VoxelGeo®, e da análise faciológica junto à perfilagem de poços do referido campo. A análise desenvolvida permitiu a individualização e a posterior visualização tridimensional de um paleocanal meandrante na base do intervalo estudado, feição esta até então não relatada em interpretações anteriores neste reservatório. Como resultado das análises sísmicas e faciológicas, foi possível elaborar um modelo deposicional, onde foram definidos quatro sistemas turbidíticos distintos, inclusos em duas seqüências de 3ª Ordem. Esses sistemas turbidíticos estariam, portanto, associados às seqüências de 4ª Ordem, que são interpretadas como parasseqüências, inseridas nos dois ciclos de 3ª Ordem. As seqüências de 3ª Ordem, que englobam os reservatórios do Campo de Namorado, representariam intervalos de alta freqüência no registro estratigráfico, dentro do contexto de afogamento (2ª Ordem) da Bacia de Campos. Pelas características da calha deposicional observada para o Campo de Namorado, é possível concluir que o sistema, como um todo, foi depositado em um complexo de canais, junto a sistemas de frentes deltaicas. Esses canais, provavelmente, foram esculpidos por fluxos hiperpicnais, formados a partir de inundações catastróficas. As informações provenientes deste estudo, proporcionaram uma melhor compreensão da gênese dos depósitos turbidíticos, acumuladores de hidrocarbonetos, no intervalo estudado, e cuja ocorrência está relacionada com etapas de rebaixamento relativo do nível do mar.

Critério constitutivo para o deslizamento com atrito ao longo da falha sísmica

A dinâmica de propagação da ruptura durante terremotos é um dos assuntos mais relevantes e complexos em Sismologia. Um critério constitutivo para a falha que descreva corretamente a relação das variáveis estáticas (tensões normais e tangenciais) com as variáveis cinéticas (deslocamentos e velocidades) na interface é necessário para efetuar uma análise dinâmica confiável. A fim de determinar um critério constitutivo para o deslizamento com atrito ao longo da falha sísmica, primeiramente apresentam-se métodos para caracterizar as superfícies deslizantes e discutem-se os principais critérios constitutivos que têm sido usados. Também são apresentados os resultados de um estudo experimental realizado, evidenciando, para sólidos metálicos em contato, uma lei constitutiva de variação do atrito com a velocidade. Um modelo numérico tridimensional baseado no Método dos Elementos Discretos (DEM) é usado para representar a região de rocha adjacente à falha. O modelo consiste de uma estrutura tridimensional periódica com massas concentradas nos nós, interconectadas por elementos visco-elásticos unidimensionais. Inicialmente, de acordo com modelos usuais em Sismologia, admite-se que o material é elástico, linear e homogêneo. Em uma segunda análise, a influência da não-homogeneidade do material é avaliada considerando que a massa específica, o módulo de Young e o coeficiente de atrito são campos aleatórios Gaussianos correlacionados. Na análise seguinte, o efeito da fratura na região de rocha adjacente à falha é também numericamente avaliado. Por fim, a influência de ruptura de micro-asperezas nas superfícies deslizantes é analisada. Através de simulação de Monte Carlo, as relações constitutivas macro (ou globais) para a falha são obtidas, admitindo como leis constitutivas micro (ou locais) os dois critérios mais usados em Sismologia: a lei de variação do atrito com a velocidade e a lei de variação do atrito com o deslizamento. Quando os blocos de rocha são admitidos serem elásticos e homogêneos não há um efeito de escala. Novamente, quando a rocha é considerada não-homogênea não há um efeito de escala significativo, apenas pequenas variações nos parâmetros das leis constitutivas macro em relação às leis micro são percebidas, exceto pela influência do campo aleatório do coeficiente de atrito, o qual apresenta um efeito de escala perceptível. A ocorrência de fratura nas proximidades da falha não causa modificações significativas nas leis constitutivas macro, apenas causa algumas oscilações em torno da lei sem fratura. Finalmente, obtém-se um critério constitutivo macro que leva em consideração a ruptura por cisalhamento das micro-asperezas nas superfícies deslizantes. O modelo é chamado lei modificada de variação do atrito com a velocidade e é considerado um critério constitutivo mais geral.

Ciclos deposionais evaporíticos da bacia de santos uma análise cicloestratigráfica a partir de dados de 2 poços e de traços de sísmica

A seqüência evaporítica Aptiana da Bacia de Santos constitui-se num dos maiores depósitos salíferos antigos conhecidos no mundo, tanto em área como em espessura. Esta última chega a atingir 2.500 m em alguns locais com baixa deformação. Devido à intensa halocinese ocorrida a partir do Albiano, sua estratificação deposicional ficou preservada em poucos locais. O Platô de São Paulo, feição fisiográfica que se estende até as águas profundas da Bacia de Santos, teve papel importante na preservação do acamadamento original. Dois poços e traços de uma seção sísmica, localizados nesta área, foram utilizados para o estudo cicloestratigráfico do pacote evaporítico. A correlação dos dois poços, a partir de seus perfis geofísicos, permitiu a identificação de um ciclo evaporítico principal com grande continuidade lateral e de espessura variável. Este ciclo deposicional foi classificado como do tipo brining upward/downward, e caracteriza-se pela precipitação rítmica de anidrita, halita, sais complexos, halita e, novamente, anidrita, no final de um ciclo e início do seguinte. Esta ordem de deposição foi interpretada como função da variação da espessura da lâmina de água e da concentração da salmoura. A análise cicloestratigráfica dos perfis geofísicos dos dois poços (raios gama, velocidade e densidade) e dos dois traços sísmicos resultou na correlação do principal ciclo identificado com o ciclo orbital de Milankovich de precessão longa, com 22 ka. Um outro ciclo secundário, com maior espessura, foi interpretado como resultado de variações climáticas em função do ciclo de obliqüidade, com um período de 39 ka. As taxas de deposição de cada ciclo são bastante variáveis em função da quantidade de cada tipo de evaporito. Através das análises cicloestratigráficas e do estudo das proporções em que cada sal ocorre, chegou-se a uma taxa de deposição da halita entre 10 e 11,5 m/ka, e de 0,61 a 0,64 m/ka para a anidrita. Esta diferença entre as taxas de deposição da halita e da anidrita explica a grande variação de espessura do ciclo principal, especialmente entre o depocentro e a borda da bacia. A menor concentração da salmoura nas proximidades da borda justifica uma menor espessura do ciclo e uma proporção bem maior de anidrita, quando se compara com uma deposição próxima ao depocentro. A cicloestratigrafia aplicada aos traços sísmicos apresentou resultados similares aos encontrados nos dois poços. A grande variação das propriedades físicas entre os diferentes tipos de evaporitos, especialmente entre a anidrita e os sais complexos, permite a intercalação de refletores com grande amplitude e continuidade. Foram contabilizados 22 ciclos nos dois traços analisados, os quais somados com o nível superior não analisado, mais rico em anidrita, totalizam 573 ka para a idade da seqüência evaporítica. Este resultado, mais os encontrados nos dois poços, permite uma previsão de tempo para a deposição de toda a seqüência entre 400 e 600 ka. Como um modelo mais provável para a deposição de toda a seqüência evaporítica, defende-se, neste trabalho, a hipótese de uma bacia profunda com lâmina de água também profunda. Os ciclos teriam se depositados a partir de mudanças de concentração da salmoura originadas pela variação de uma lâmina de água de cerca de 1000 m, na deposição da anidrita, para menos de 200 na deposição dos sais complexos, na parte mais profunda da bacia. Corroboram para este modelo, a grande continuidade, do depocentro até as bordas, das camadas de cada ciclo, a declividade, mesmo que bastante baixa, da bacia e a impossibilidade de valores de taxas de subsidência similares aos de deposição do ciclo evaporítico.

Interpretação integrada de sísmica de alta-resolução e da morfologia submarina da costa de fiordes na patagônia central - Chile

Este estudo apresenta o uso integrado da sísmica de alta resolução e da morfologia submarina para interpretar a evolução do ambiente glacimarinho da costa de fiordes da Patagônia Central, Chile. Foram analisados registros de perfilador de fundo e subfundo 3,5 kHz e modelos submarinos 3D de sete fiordes adjacentes ao campo de gelo Patagônico Sul (Eyre, Falcon, Penguin, Europa Peel, Calvo e Amalia) e parte do canal Icy. Os registros, com cerca de 300 km de levantamento acústico, foram obtidos pelo Servicio Hidrográfico y Oceanográfico de la Armada de Chile (SHOA), durante o cruzeiro de investigação científica Campo de Hielo Sur, realizado em 1995. Foram identificadas as principais fácies acústicas e geoformas sedimentares. A morfologia submarina e das bacias subaéreas adjacentes foi analisada com a elaboração de modelos tridimensionais subaéreos e submarinos e de perfis batimétricos longitudinais a partir de cartas náuticas, também do SHOA. Foram utilizadas imagens Landsat ETM+ na interpretação da geomorfologia glacial da área de entorno subaérea dos fiordes. O eixo longitudinal dos fiordes exibe morfologia irregular com bacias profundas e mostra fácies acústicas associadas ao sistema de depósitos de zonas de linha de encalhe (grounding line), às línguas de gelos flutuantes, e aos icebergs e ao gelo marinho. Refletores acústicos distinguem duas fácies principais, segundo sua configuração interna e geometria externa: caóticas e estratificadas. A geometria dos depocentros e as características dos refletores acústicos indicam a importante influência da batimetria e da topografia pré-existentes na dinâmica das geleiras e nos conseqüentes processos de sedimentação. Devido às grandes profundidades das bacias, as frentes das geleiras poderiam ser flutuantes ou aterradas ao fundo marinho. Em todo o caso, predomina um regime glacial onde o gelo está perto do ponto de fusão. É sugerido que os fiordes estudados resultam de um continuum de formas e de processos, que vão dos lineamentos e sistemas de falhas pré-existentes e controlados tectonicamente desde o Mioceno Inferior (ca. 25 milhões de anos AP), processos de denudação fluvial e de vertentes, vales fluviais modificados glacialmente e, finalmente, canais e fiordes erodidos pela ação do gelo. A tectônica originou uma topografia favorável para o desenvolvimento das geleiras, ocorreram processos de retroalimentação positiva entre aumento o da precipitação e expansão do campo de gelo no lado oeste da cordilheira. Provavelmente, as geleiras durante o Último Máximo Glacial não foram suficientemente espessas para aprofundar os fiordes, tampouco para dragar seus depósitos sedimentares para fora deles. Os processos de deposição de sedimentos caóticos acusticamente visíveis ocorreram durante o recuo das geleiras, já no Holoceno, onde elas alcançaram pontos de estabilidade, ajudadas pela morfologia dos fiordes, mesmo em águas profundas. Os depósitos estratificados localizados nas bacias intra-sills e em maiores profundidades, por sua vez, tiveram sua origem, provavelmente, das partes flutuantes do gelo, onde predominam os processos de desprendimento de icebergs (calving) e também do gelo marinho. Nessas mesmas bacias, depósitos mais profundos, não visíveis ao sistema de alta-resolução, poderiam estar preservados no fundo, pois não foram erodidos pelos sucessivos ciclos de avanço e recuo das geleiras, contendo assim informações sobre a evolução do campo de gelo ao longo do Quaternário.

Análise sísmica usando transformada de Curvelet

Oil prospecting is one of most complex and important features of oil industry Direct prospecting methods like drilling well logs are very expensive, in consequence indirect methods are preferred. Among the indirect prospecting techniques the seismic imaging is a relevant method. Seismic method is based on artificial seismic waves that are generated, go through the geologic medium suffering diffraction and reflexion and return to the surface where they are recorded and analyzed to construct seismograms. However, the seismogram contains not only actual geologic information, but also noise, and one of the main components of the noise is the ground roll. Noise attenuation is essential for a good geologic interpretation of the seismogram. It is common to study seismograms by using time-frequency transformations that map the seismic signal into a frequency space where it is easier to remove or attenuate noise. After that, data is reconstructed in the original space in such a way that geologic structures are shown in more detail. In addition, the curvelet transform is a new and effective spectral transformation that have been used in the analysis of complex data. In this work, we employ the curvelet transform to represent geologic data using basis functions that are directional in space. This particular basis can represent more effectively two dimensional objects with contours and lines. The curvelet analysis maps real space into frequencies scales and angular sectors in such way that we can distinguish in detail the sub-spaces where is the noise and remove the coefficients corresponding to the undesired data. In this work we develop and apply the denoising analysis to remove the ground roll of seismograms. We apply this technique to a artificial seismogram and to a real one. In both cases we obtain a good noise attenuation

Anisotropia sísmica crustal na região de Cascavel - CE

This dissertation presents a study on crustal seismic anisotropy in Cascavel - CE. The earthquake data employed here are from the Seismological Laboratory at Universidade Federal do Rio Grande do Norte (UFRN) and were colected from 29 September 1997 to 05 march 1998 using six three-component digital seismographic stations. In general, the cause of the observed seismic anisotropy in many regions of the world is interpreted in terms of fluid-filled stress aligned microcracks in the rockmass (EDA). In other words, the polarisation direction of the faster shear-wave splitting is parallel to SHmax. However, other researches on seismic anisotropy carried out in NE Brazil have shown a remarkable consistency of the faster shear-wave polarisation direction with the direction of the Precambrian fabric. The present work is another case study that is used to investigate this issue. In order to map the Precambrian fabric we used aeromagnetic data, since the study area is mostly covered with sediments (up to 50m thick) and in-situ field mapping would be very difficult to be carried out. According to the results from the present research, the observations of the faster shear-wave polarisation directions in two seismographic stations in Cascavel region are best explained in the framework of EDA. For the remaining two stations, the observed anisotropy may have two interpretions: (i) - 90_ flips of the direction of polarisation of the faster shear-wave, since that the event-to-station ray path would be through the fracture zone and hence would travel under a higher pore pressure and (ii) - the observed seismic anisotropy would agree with the direction on the ductile Precambrian fabric

Estudo da atividade sísmica em São Caetano-PE em 2007

In this dissertation we studied the seismic activity in the São Caetano county, Pernambuco State, Northeastern Brazil, located near the Pernambuco Lineament. The Pernambuco Lineament is a one of Neoproterozoic continental-scale shear zones that deforms the Borborema province. The seismicity estudied occurred in a NE trending branch of Pernambuco Lineament. The seismic activity in São Caetano started in 2006 and in May 20th, 2006 a 4,0 mb earthquake hit there. This was the largest earthquake ever reported in Pernambuco State. This dissertation is the result of a campaign done in the period from Februay 1th 2007 to July 31 th 2007. In this campaign up to nine three-component digital seismographic stations were deployed and the collected data was used to determine hypocenters and focal mechanism. A total of 214 earthquakes, recorded at least by three stations, were analyzed. To determine hypocenters and time origin the HYPO71 program was used assuming a half-space model with parameters : VP (P-wave velocity) equal to 5.90 km/s and the ratio VP/VS 1.70, where VS is the S-wave velocity. The earthquakes hypocentral distribution was approximately 4 km long and agrees with the NE-SW direction of the Pernambuco Lineamento branch. Hypocentres depth range from 2 to 8 km. The composed focal mechanism was made from a group of 14 selected earthquakes. We try firstly to find the fault plane solution matching the polarity distribution at stations, using the FPFIT program. The result was 43 deg ± 15 deg for strike, 59 deg ± 9 deg for dip and -142 deg ± 15 deg for rake. We also fitted a plane using the hypocentral distribution to obtain the dip and azimuth of the hypocentral distribution. The results obtained by this fit were 58 deg for the azimuth, 55 deg for the dip and -155 deg for rake. This result shows a mechanism of a strike-slip dextral fault with a normal component. This dissertation shows, once more, that there is a good correlation between the seismic activity and geological features in the region near the Pernambuco Lineament and its NE branches

Aquisição, processamento e análise de dados de sísmica de alta resolução na plataforma continental norte do Rio Grande do Norte: Vale inciso do rio açu

On the modern Continental Shelf to the north of Rio Grande do Norte state (NE Brazil) is located a paleo-valley, submerged during the last glacial sea-level lowstand, that marks continuation of the most important river of this area (Açu River). Despite the high level of exploration activity of oil industry, there is few information about shallow stratigraphy. Aiming to fill this gap, situated on the Neogene, was worked a marine seismic investigation, the development of a processing flow for high resolution data seismic, and the recognition of the main feature morphology of the study area: the incised valley of the River Açu. The acquisition of shallow seismic data was undertaken in conjunction with the laboratory of Marine Geology/Geophysics and Environmental Monitoring - GGEMMA of Federal University of Rio Grande do Norte UFRN, in SISPLAT project, where the geomorphological structure of the Rio paleovale Açu was the target of the investigation survey. The acquisition of geophysical data has been over the longitudinal and transverse sections, which were subsequently submitted to the processing, hitherto little-used and / or few addressed in the literature, which provided a much higher quality result with the raw data. Once proposed for the flow data was developed and applied to the data of X-Star (acoustic sensor), using available resources of the program ReflexW 4.5 A surface fluvial architecture has been constructed from the bathymetric data and remote sensing image fused and draped over Digital Elevation Models to create three-dimensional (3D) perspective views that are used to analyze the 3D geometry geological features and provide the mapping morphologically defined. The results are expressed in the analysis of seismic sections that extend over the region of the continental shelf and upper slope from mouth of the Açu River to the shelf edge, providing the identification / quantification of geometrical features such as depth, thickness, horizons and units seismic stratigraphyc area, with emphasis has been placed on the palaeoenvironmental interpretation of discordance limit and fill sediment of the incised valley, control by structural elements, and marked by the influence of changes in the sea level. The interpretation of the evolution of this river is worth can bring information to enable more precise descriptions and interpretations, which describes the palaeoenvironmental controls influencing incised valley evolution and preservation to provide a better comprehensive understanding of this reservoir analog system

Interpretação sísmica e análise de atributos - área do levantamento sísmico 3D de siririzinho (Bacia Sergipe-Alagoas, NE do Brasil)

The main objective of the present thesis was the seismic interpretation and seismic attribute analysis of the 3D seismic data from the Siririzinho high, located in the Sergipe Sub-basin (southern portion of Sergipe-Alagoas Basin). This study has enabled a better understanding of the stratigraphy and structure that the Siririzinho high experienced during its development. In a first analysis, we used two types of filters: the dip-steered median filter, was used to remove random noise and increase the lateral continuity of reflections, and fault-enhancement filter was applied to enhance the reflection discontinuities. After this filtering step similarity and curvature attributes were applied in order to identify and enhance the distribution of faults and fractures. The use of attributes and filtering greatly contributed to the identification and enhancement of continuity of faults. Besides the application of typical attributes (similarity and curvature) neural network and fingerprint techniques were also used, which generate meta-attributes, also aiming to highlight the faults; however, the results were not satisfactory. In a subsequent step, well log and seismic data analysis were performed, which allowed the understanding of the distribution and arrangement of sequences that occur in the Siririzinho high, as well as an understanding of how these units are affected by main structures in the region. The Siririzinho high comprises an elongated structure elongated in the NS direction, capped by four seismo-sequences (informally named, from bottom to top, the sequences I to IV, plus the top of the basement). It was possible to recognize the main NS-oriented faults, which especially affect the sequences I and II, and faults oriented NE-SW, that reach the younger sequences, III and IV. Finally, with the interpretation of seismic horizons corresponding to each of these sequences, it was possible to define a better understanding of geometry, deposition and structural relations in the area.

Análise da compacidade do solo na área urbana de Bauru/ Estado de São Paulo utilizando Ensaios de Penetração Padronizado e sísmica de refração

This study was done to determine compactness of soil strata in the urban area of Bauru, central São Paulo State using standard penetration test (SPT) and seismic refraction methods. Frequency distribution within boreholes was measured every meter to a depth of 30 meters. The geometric mean of standard penetration test values was computed every meter. A biplot graph showing the relationship of SPT with depth indicates soil stratification. In addition, five sections of seismic refraction values, which also indicate the compactness of soil strata, were obtained. A comparison of the two methods shows significant correlation between the results obtained from each.

Desenvolvimento de sistema para sísmica up-hole em conjunto com o ensaio SPT

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Modelagem sísmica via métodos das diferenças finitas: caso da bacia do Amazonas

Este trabalho tem por objetivo apresentar os resultados da modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações e múltiplas reflexões, tendo a Bacia do Amazonas como referência à modelagem. As condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas foram analisadas numericamente pelo método das diferenças finitas, visando melhor compreensão e controle da interpretação de dados sísmicos. A geologia da Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o Recente que atingem 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çam espessuras de centenas de metros e perfazem um volume total de aproximadamente 90000 Km³. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica o que impossibilita melhor interpretação dos horizontes refletores que se encontram abaixo destas soleiras. Para representar situações geológicas desse tipo foram usados um modelo (sintético) acústico de velocidades e um código computacional elaborado via método das diferenças finitas com aproximação de quarta ordem no espaço e no tempo da equação da onda. A aplicação dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas melhorou 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 observada a influência significativa das reflexões múltiplas devido à camada de alta velocidade, isto provocou maior perda de energia e dificultou a interpretação dos alvos. Por esta razão recomenda-se a integração de dados de superfície com os de poço, com o objetivo de obter melhor imagem dos alvos abaixo das soleiras de diabásio.

Atenuação de múltiplas e compressão do pulso fonte em dados de sísmica de reflexão utilizando o filtro Kalman-Bucy

O objetivo central deste trabalho é o estudo e a aplicação do método Kalman-Bucy no processo de deconvolução ao impulso e de deconvolução com predição, onde é considerado que os dados observados são classificados como não-estacionários. Os dados utilizados neste trabalho são sintéticos e, com isto, esta Tese tem características de um exercício numérico e investigativo. O operador de deconvolução ao impulso é obtido a partir da teoria de CRUMP (1974) fazendo uso das soluções das equações Wiener-Hopf apresentadas por KALMAN-BUCY (1961) nas formas contínuas e discretas considerando o processo como não estacionário. O operador de predição (KBCP) está baseado nas teorias de CRUMP (1974) e MENDEL ET AL (1979). Sua estrutura assemelha-se ao filtro Wiener-Hopf onde os coeficientes do operador (WHLP) são obtidos através da autocorrelação, e no caso (KBCP) são obtidos a partir da função b_i(k). o problema é definido em duas etapas: a primeira consta da geração do sinal, e a segunda da sua avaliação. A deconvolução realizada aqui é classificada como estatística, e é um modelo fortemente baseado nas propriedades do sinal registrado e de sua representação. Os métodos foram aplicados apenas em dados sintéticos de seção fonte-comum obtida a partir dos modelos com interfaces contínuas e camadas homogêneas.

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.