3D raytracing through homogeneous anisotropic media with smooth interfaces

O raio conectando dois pontos em um meio anisotrópico, homogêneo por partes e com variação lateral, é calculado utilizando-se técnicas de continuação em 3D. Se combinado com algoritmos para solução do problema de valor inicial, o método pode ser estendido para o cálculo de eventos qS₁ e qS₂. O algoritmo apresenta a mesma eficiência e robustez que implementações de técnicas de continuação em meios isotrópicos. Rotinas baseadas neste algoritmo têm várias aplicações de interesse. Primeiramente, na modelagem e inversão de parâmetros elásticos na presença de anisotropia. Em segundo lugar, as iterações de Newton-Raphson produzem atributos da frente de onda como vetor vagarosidade e a matrix hessiana do tempo de trânsito, quantidades que permitem determinar o espalhamento geométrico e aproximações de segunda ordem para o tempo de trânsito. Estes atributos permitem calcular as amplitudes ao longo do raio e investigar os efeitos da anisotropia no empilhamento CRS em modelos de velocidade simples.

Estudo clínico do tempo de fechamento de espaço e do movimento dentário durante a retratação de caninos entre dois tipos de braquetes

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

Análise experimental sobre velocidade de captura em sistemas de transporte pneumático

A velocidade de transporte é um dos mais importantes parâmetros no transporte pneumático de sólidos. O êxito no transporte de materiais particulados depende da predição ou da determinação da velocidade mínima de transporte. Uma velocidade acima daquela necessária ao transporte estável das partículas sólidas conduz a um grande consumo de energia devido ao aumento na perda de pressão do sistema, degradação dos sólidos e abrasão dos sólidos e da tubulação. Já uma velocidade abaixo desse valor limite certamente resultará na deposição das partículas sólidas para o fundo da tubulação e conseqüentemente o entupimento desta. Neste trabalho, uma técnica para medir a velocidade mínima de captura de partículas sólidas em uma tubulação na direção horizontal em um sistema de Transporte Pneumático é desenvolvida. Ela baseia-se em observações visuais do comportamento das partículas ao ocorrer a captura, em medidas da velocidade de operação do gás e da massa das partículas capturadas. É realizada ainda a análise qualitativa de alguns parâmetros que influenciam na velocidade de captura das partículas, permitindo uma maior compreensão dos fenômenos envolvidos.

Análise de velocidade por otimização do semblance na reflexão sísmica

Este trabalho teve como objetivo geral desenvolver uma metodologia sistemática para a inversão de dados de reflexão sísmica em arranjo ponto-médio-comum (PMC), partindo do caso 1D de variação vertical de velocidade e espessura que permite a obtenção de modelos de velocidades intervalares, v_int,n, as espessuras intervalares, z_n, e as velocidades média-quadrática, v_RMS,n, em seções PMC individualizadas. Uma consequência disso é a transformação direta destes valores do tempo para profundidade. Como contribuição a análise de velocidade, foram desenvolvidos dois métodos para atacar o problema baseado na estimativa de velocidade intervalar. O primeiro método foi baseado na marcação manual em seções PMC, e inversão por ajuste de curvas no sentido dos quadrados-mínimos. O segundo método foi baseado na otimização da função semblance para se obter uma marcação automática. A metodologia combinou dois tipos de otimização: um Método Global (Método Price ou Simplex), e um Método Local (Gradiente de Segunda Ordem ou Conjugado), submetidos a informação à priori e vínculos. A marcação de eventos na seção tempo-distância faz parte dos processos de inversão, e os pontos marcados constituem os dados de entrada juntamente com as informações à priori do modelo a ser ajustado. A marcação deve, por princípio, evitar eventos que representem múltiplas, difrações e interseções, e numa seção pode ser feita mais de 50 marcações de eventos, enquanto que num mapa semblance não se consegue marcar mais de 10 eventos de reflexão. A aplicação deste trabalho é voltada a dados sísmicos de bacias sedimentares em ambientes marinhos para se obter uma distribuição de velocidades para a subsuperfície, onde o modelo plano-horizontal é aplicado em seções PMC individualizadas, e cuja solução pode ser usada como um modelo inicial em processos posteriores. Os dados reais da Bacia Marinha usados neste trabalho foram levantados pela PETROBRAS em 1985, e a linha sísmica selecionada foi a de número L5519 da Bacia do Camamu, e o PMC apresentado é a de número 237. A linha é composta de 1098 pontos de tiro, com arranjo unilateraldireito. O intervalo de amostragem é 4 ms. O espaçamento entre os geofones é 13,34 m com o primeiro geofone localizado a 300 m da fonte. O espaçamento entre as fontes é de 26,68 m. Como conclusão geral, o método de estimativa de velocidade intervalar apresentada neste trabalho fica como suporte alternativo ao processo de análise de velocidades, onde se faz necessário um controle sobre a sequência de inversão dos PMCs ao longo da linha sísmica para que a solução possa ser usada como modelo inicial ao imageamento, e posterior inversão tomográfica. Como etapas futuras, podemos propor trabalhos voltados direto e especificamente a análise de velocidade sísmica estendendo o caso 2D de otimização do semblance ao caso 3D, estender o presente estudo para o caso baseado na teoria do raio imagem com a finalidade de produzir um mapa continuo de velocidades para toda a seção sísmica de forma automática.

Aplicação de análise de séries temporais em perfis de poço

Esta tese descreve a aplicação de análise de séries temporais em perfis de poço. Através desta técnica é possível avaliar-se a repetição e a resolução vertical de perfis, e determinar-se o intervalo de amostragem e a velocidade de perfilagem ideais para diferentes perfis. A comparação entre três poços é também feita, baseada num mesmo tipo de perfil. Para tanto, na seqüência utilizada, procurou-se manter num mesmo domínio os dados cuja quantidade total de amostras (N) por perfil não ultrapassou 2048. Desses dados, foram inicialmente retirados o valor médio das amostras e o alinhamento polinomial algébrico eventualmente nelas embutido. Em seguida, foi efetuada a aplicação do ponderador cossenoidal, do filtro passa-alta, da janela Hanning, do cálculo da função coerência, do espectro de fase, da razão sinal-ruído e dos espectros de potência do sinal e do ruído, nesta ordem. Para a função coerência, fez-se necessário o cálculo dos níveis de confiança de 50%, 90% e 95%. O cálculo do primeiro nível teve por base a necessidade de se determinar a resolução vertical de alguns perfis, e dos demais, a fim de que fosse obtida uma informação referente à localização daqueles níveis para a coerência calculada. Em relação ao espectro de fase, seu cálculo surgiu da necessidade de se obter uma informação adicional a respeito dos perfis manipulados, ou seja, o conhecimento da ocorrência ou não de deslocamento relativo de profundidade entre a seção principal e a seção repetida. A razão sinal-ruído foi calculada no sentido de possibilitar a comparação, como elemento avaliador dos diversos tipos de perfis, com a coerência e o cálculo dos espectros de potência. Os espectros de potência do sinal e do ruído foram calculados para se ter mais um parâmetro de avaliação da seção repetida, já que em tese, os espectros de potência do sinal e do ruído da seção repetida devem ser iguais aos respectivos espectros da seção principal. Os dados utilizados na aplicação da metodologia proposta foram fornecidos pela PETROBRÁS e oriundos de quatro poços da Bacia Potiguar emersa. Por questões de sigilo empresarial, os poços foram identificados como poços A, B, C e D. A avaliação da repetição entre diferentes tipos de perfis indica que, para o poço A, o perfil micro-esférico (MSFL) tem repetição melhor do que o perfil de porosidade neutrônica (CNL), o qual tem, por sua vez, repetição melhor do que o perfil de raios gama normal (GR). Para os perfis do poço D, uma diminuição da velocidade de perfilagem de 550 m/h para 275 m/h é vantajosa apenas para o perfil de porosidade neutrônica. Já a velocidade de perfilagem de 920 m/h, utilizada, na obtenção dos perfis do poço C, é totalmente inadequada para os perfis de resistividade (MSFL, ILD e ILM). A diminuição do intervalo de amostragem de 0,20 m para 0,0508 m, nos perfis de raios gama e de porosidade neutrônica, e 0,0254 m para o perfil de densidade, apresenta bons resultados quando aplicada no poço D. O cálculo da resolução vertical indica, para o perfil de porosidade neutrônica, uma superioridade qualitativa em relação ao perfil de raios gama normal, ambos pertencentes ao poço A. Para o poço C, o perfil micro-esférico apresenta uma resolução vertical na mesma ordem de grandeza da resolução do perfil de raios gama do poço B, o que evidencia ainda mais a inconveniência da velocidade de perfilagem utilizada no poço C. Já para o poço D, o cálculo da resolução vertical indica uma superioridade qualitativa do perfil de densidade de alta resolução em relação ao perfil de raios gama de alta resolução. A comparação entre os poços A, B e D, levada a efeito através dos respectivos perfis de porosidade neutrônica normais, comprova que a presença de ruído aleatório, em geral, está diretamente ligada à porosidade da formação - uma maior porosidade indica uma presença maior de ruído e, por conseguinte, uma queda qualitativa no perfil obtido. A análise do espectro de fase de cada perfil indica um deslocamento em profundidade, existente entre as seções principal e repetida de todos os perfis do poço C. E isto pôde ser confirmado com a posterior superposição das seções.

Avaliação de algoritmos para conversão de modelos de velocidade de tempo para profundidade

Ainda hoje, a migração em tempo é o processo de imageamento substancialmente empregado na indústria do petróleo. Tal popularidade é devida ao seu alto grau de eficiência e robustez, além de sua habilidade em focalizar refletores nos mais variados ambientes geológicos. Entretanto, em áreas de alta complexidade geológica a migração em tempo falha de tal forma que a migração em profundidade e um campo de velocidade em profundidade são indispensáveis. Esse campo é geralmente obtido através de processos tomográficos partindo de um campo de velocidade inicial. A conversão de campos de velocidade de tempo para profundidade é uma forma rápida de se obter um campo inicial mais consistente geologicamente para tais processos. Alguns algoritmos de conversão tempo-profundidade recentemente desenvolvidos baseados no traçamento de raios-imagem são revistos e um algoritmo alternativo baseado na propagação da frente de onda-imagem é proposto. Os algoritmos são aplicados a dados sintéticos bidimensionais e avaliados de acordo com suas eficiência e acurácia, destacando suas vantagens, desvantagens e limitações na obtenção de campos de velocidade em profundidade.

Influência do cálculo subgengival nas medidas de sondagem periodontal com sonda manual e eletrônica de força controlada: estudo em humanos

Pós-graduação em Odontologia - FOAR

Um estudo da estabilidade mecânica de poços a partir de perfis geofísicos em formações pouco consolidadas

A análise da estabilidade mecânica de um poço pode ser feita a partir do cálculo de parâmetros elásticos da formação utilizando a densidade do meio e as velocidades de propagação das ondas compressional e cisalhante na formação rochosa, os quais podem ser obtidos de perfis geofísicos do poço. Em formações sedimentares pouco consolidadas as ferramentas de perfilagem sônica convencionais (monopolares) não conseguem registrar com acuidade a velocidade da onda cisalhante pois a primeira chegada dessa onda é camuflada pela chegada de outras ondas que podem ser mais rápidas que a onda cisalhante num poço perfurado neste tipo de formação. Medidas das velocidades sônicas são feitas em laboratório em amostras da formação, sob condições semelhantes às condições in situ, servindo como ajuste das velocidades registradas no poço pela ferramenta de perfilagem sônica. Para a análise de estabilidade da formação, perfis auxiliares são necessários como o perfil de porosidade, saturação de fluidos e perfis de composição mineralógica da formação rochosa. Exige-se ainda dados de testes de avaliação da formação e de condições do reservatório, mas que são comuns em poços de petróleo, como o teste de formação e os testes de pressurização do poço, tais como o teste de micro-fraturamento hidráulico ou o teste de absorção. A avaliação das tensões principais efetivas que atuam distante do poço, e que não são afetadas pela sua presença, é feita através da associação de um modelo de deformação elástica apropriado e o resultado do teste de pressurização disponível para o poço em estudo. Utilizando resultados clássicos da teoria da elasticidade geral pode-se calcular o campo de tensões modificado na vizinhança da parede do poço devido ao efeito da própria presença do poço ali perfurado e da diferença de pressão existente entre o interior do poço e a formação rochosa. A determinação das propriedades mecânicas da formação a partir das velocidades sônicas e a avaliação do campo de tensões assumindo um modelo elástico de deformação, supõem o meio rochoso no qual as ondas se propagam como elástico, homogêneo e isotrópico. Esta suposição representa a principal aproximação assumida pela metodologia descrita neste trabalho. De posse das propriedades mecânicas da formação e do campo de tensões que age na vizinhança do poço resta definir o critério segundo o qual a rocha sofre instabilidade mecânica quando submetida aquele campo de tensões. Isto permite determinar se, nas condições avaliadas do poço e da formação, haverá quebra da parede do poço por excesso de tensão e, se houver, qual a sua extensão. Assim o problema é como analisar o comportamento mecânico de um poço em uma formação pouco consolidada a partir de perfis geofísicos os quais podem ter problemas no registro das propriedades físicas do meio em formações deste tipo. A metodologia proposta é aplicada a dois intervalos de profundidade pertencentes a dois poços onde arenitos e folhelhos se intercalam e nos quais todos os dados necessários estão disponíveis. Os resultados obtidos mostram que, exceto quando outros mecanismos de quebra da parede do poço agem na formação, a metodologia proposta consegue com sucesso detectar zonas de ocorrência de instabilidade mecânica do poço provocadas por um campo de tensões que excede a resistência mecânica da formação.

Imageamento multifoco de refletores sísmicos

A simulação de uma seção sísmica de afastamento-nulo (ZO) a partir de dados de cobertura múltipla para um meio 2-D, através do empilhamento, é um método de imageamento de reflexão sísmica muito utilizado, que permite reduzir a quantidade de dados e melhorar a relação sinal/ruído. Segundo Berkovitch et al. (1999) o método Multifoco está baseado na Teoria do Imageamento Homeomórfico e consiste em empilhar dados de cobertura múltipla com distribuição fonte-receptor arbitrária de acordo com uma nova correção de sobretempo, chamada Multifoco. Esta correção de sobretempo esta baseada numa aproximação esférica local da frente de onda focalizante na vizinhança da superfície da terra. Este método permite construir uma seção sísmica no domínio do tempo de afastamento nulo aumentando a relação sinal/ruído. A técnica Multifoco não necessita do conhecimento a priori de um macro-modelo de velocidades. Três parâmetros são usados para descrever a aproximação de tempo de trânsito, Multifoco, os quais são: 1) o ângulo de emergência do raio de afastamento nulo ou raio de reflexão normal (β₀), 2) a curvatura da frente de onda no Ponto de Incidência Normal (R_NIP) e 3) curvatura da frente de Onda Normal (R_N). Sendo também necessário a velocidade próximo a superfície da terra. Neste trabalho de tese aplico esta técnica de empilhamento Multifoco para dados de cobertura múltipla referidos a modelos de velocidade constante e modelo heterogêneos, com o objetivo de simular seções sísmicas afastamento-nulo. Neste caso, como se trata da solução de um problema direto, o macro-modelo de velocidades é considerado conhecido a priori. No contexto do problema inverso tem-se que os parâmetros R_NIP, R_N e β₀ podem ser determinados a partir da análise de coerência aplicada aos dados sísmicos de múltipla cobertura. Na solução deste problema a função objetivo, a ser otimizada, é definida pelo cálculo da máxima coerência existente entre os dados na superfície de empilhamento sísmico. Neste trabalho de tese nos discutimos a sensibilidade da aproximação do tempo de trânsito usado no empilhamento Multifoco, como uma função dos parâmetros R_NIP, R_N e β₀. Esta análise de sensibilidade é feita de três diferentes modos: 1) a primeira derivada da função objetivo, 2) a medida de coerência, denominada semelhança, e 3) a sensibilidade no Empilhamento Multifoco.

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.

Medidas de coerência para análise de velocidade na migração em tempo

Métodos iterativos para a análise de velocidade com migração dependem das funções objetivo medirem a horizontalização dos eventos de reflexão em famílias de ponto imagem comum (CIG). A migração em tempo é um método simples de imageamento para avaliar essas funções objetivo. Usando migração em tempo, estudamos a influência da função objetivo nos resultados da análise de velocidade com migração. Propomos duas novas funções objetivo para a análise velocidade com migração: Semblance Diferencial Estendido e o produto do Semblance Clássico com o Semblance Diferencial Estendido. Experimentos numéricos, utilizando os dados Marmousoft, mostram a eficácia das novas funções objetivo em estimar modelos velocidade produzindo eventos horizontalmente alinhados nas famílias de ponto imagem comum.

Migração com amplitude verdadeira em meios com gradiente constante de velocidade

Este trabalho tem por objetivo a aplicação de um método de migração com amplitudes verdadeiras, considerando-se um meio acústico onde a velocidade de propagação varia linearmente com a profundidade. O método de migração é baseado na teoria dos raios e na integral de migração de Kirchhoff, procurando posicionar de forma correta os refletores e recuperar os respetivos coeficientes de reflexão. No processo de recuperação dos coeficientes de reflexão, busca-se corrigir o fator de espalhamento geométrico de reflexões sísmicas primárias, sem o conhecimento a priori dos refletores procurados. Ao considerar-se configurações fonte-receptor arbitrárias, as reflexões primárias podem ser imageadas no tempo ou profundidade, sendo as amplitudes do campo de ondas migrado uma medida dos coeficientes de reflexão (função do ângulo de incidência). Anteriormente têm sido propostos alguns algoritmos baseados na aproximação de Born ou Kirchhoff. Todos são dados em forma de um operador integral de empilhamento de difrações, que são aplicados à entrada dos dados sísmicos. O resultado é uma seção sísmica migrada, onde cada ponto de reflexão é imageado com uma amplitude proporcional ao coeficiente de reflexão no ponto. No presente caso, o processo de migração faz uso de um modelo com velocidade que apresenta uma distribuição que varia linearmente com a profundidade, conhecido também como gradiente constante de velocidade. O esquema de migração corresponde a uma versão modificada da migração de empilhamento por difração e faz uso explícito da teoria do raio, por exemplo, na descrição de tempos de trânsito e amplitudes das reflexões primárias, com as quais a operação de empilhamento e suas propriedades podem ser entendidas geometricamente. Efeitos como o espalhamento geométrico devido à trajetória do raio levam a distorção das amplitudes. Estes efeitos têm que ser corregidos durante o processamento dos dados sísmicos. Baseados na integral de migração de Kirchhoff e na teoria paraxial dos raios, foi derivada a função peso e o operador da integral por empilhamento de difrações para um modelo sísmico 2,5-D, e aplicado a uma serie de dados sintéticos em ambientes com ruído e livre de ruído. O resultado mostra a precisão e estabilidade do método de migração em um meio 2,5-D como ferramenta para obter informação sobre as propriedades de refletividade da subsuperfície da terra. Neste método não são levados em consideração a existência de caústicas nem a atenuação devido a fricção interna.

Análise do limite da capacidade de escoamento de esgoto sanitário em tubulação interceptora assentada em área urbana

O objetivo deste trabalho foi analisar o limite da capacidade de escoamento de esgoto sanitário em tubulação interceptora assentada na área urbana do município de Belém – PA. Para isso, foram considerados a Bacia de Esgotamento 5 (BE 5), as projeções populacionais e de produção de esgoto (2010-2030) da área oficial com coleta e transporte de esgoto da BE 5 (Área 1) e da área total dessa Bacia (Área 2), os 100m (cem) do interceptor assentado na BE 5 (chegada da Estação Elevatória de Esgoto do Una), bem como as informações identificadas dessa tubulação e a alteração de declividade proposta neste trabalho. Foi executado o dimensionamento hidráulico e verificou-se que na Área 1, para todas as declividades escolhidas, o interceptor terá capacidade de escoamento no período de 2010-2030, chegando a lâmina líquida máxima de 0,225 em 2030. No entanto, para a Área 2 verificou-se que a capacidade de escoamento terá comprometimento, quando utilizamos três das sete declividades estimadas, chegando a lâmina líquida máxima de 0,950 em 2030, indo de encontro ao recomentado na NBR 9649/1986 Projeto de Rede Coletora de Esgoto Sanitário (0,750). Recomenda-se, neste caso, a reestruturação/ampliação do sistema de coleta e transporte de esgoto da BE 5, em virtude do aumento populacional e de produção de esgoto, resultante da nova divisão de Bacias de Esgotamento realizado no PDSES (2010).

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.