Inversão de dados sísmicos de reflexão a partir da curva do tempo de trânsito

A presente Dissertação de Mestrado tem como objetivo o estudo do problema de inversão sísmica baseada em refletores planos para arranjo fonte-comum (FC) e ponto-médiocomum (PMC). O modelo direto é descrito por camadas homogêneas, isotrópicas com interfaces plano-horizontais. O problema é relacionado ao empilhamento NMO baseado na otimização da função semblance, para seções PMC corrigidas de sobretempo normal (NMO). O estudo foi baseado em dois princípios. O primeiro princípio adotado foi de combinar dois grupos de métodos de inversão: um Método Global e um Método Local. O segundo princípio adotado foi o de cascata, segundo a teoria Wichert-Herglotz-Bateman, que estabelece que para conhecer uma camada inferior tem-se que conhecer primeiro a camada superior (dissecação). A aplicação do estudo é voltada à simulação sísmica de Bacia Sedimentar do Solimões e de Bacia Marinha para se obter uma distribuição local 1D de velocidades e espessuras para a subsuperfície em horizontes alvo. Sendo assim, limitamos a inversão entre 4 e 11 refletores, uma vez que na prática a indústria limita uma interpretação realizada apenas em número equivalente de 3 a 4 refletores principais. Ressalta-se que este modelo é aplicável como condição inicial ao imageamento de seções sísmicas em regiões geologicamente complexas com variação horizontal suave de velocidades. Os dados sintéticos foram gerados a partir dos modelos relacionados a informações geológicas, o que corresponde a uma forte informação a priori no modelo de inversão. Para a construção dos modelos relacionados aos projetos da Rede Risco Exploratório (FINEP) e de formação de recursos humanos da ANP em andamento, analisamos os seguintes assuntos relevantes: (1) Geologia de bacias sedimentares terrestre dos Solimões e ma rinha (estratigráfica, estrutural, tectônica e petrolífera); (2) Física da resolução vertical e horizontal; e (3) Discretização temporal-espacial no cubo de multi-cobertura. O processo de inversão é dependente do efeito da discretização tempo-espacial do campo de ondas, dos parâmetros físicos do levantamento sísmico, e da posterior reamostragem no cubo de cobertura múltipla. O modelo direto empregado corresponde ao caso do operador do empilhamento NMO (1D), considerando uma topografia de observação plana. O critério básico tomado como referência para a inversão e o ajuste de curvas é a norma 2 (quadrática). A inversão usando o presente modelo simples é computacionalmente atrativa por ser rápida, e conveniente por permitir que vários outros recursos possam ser incluídos com interpretação física lógica; por exemplo, a Zona de Fresnel Projetada (ZFP), cálculo direto da divergência esférica, inversão Dix, inversão linear por reparametrização, informações a priori, regularização. A ZFP mostra ser um conceito út il para estabelecer a abertura da janela espacial da inversão na seção tempo-distância, e representa a influência dos dados na resolução horizontal. A estimativa da ZFP indica uma abertura mínima com base num modelo adotado, e atualizável. A divergência esférica é uma função suave, e tem base física para ser usada na definição da matriz ponderação dos dados em métodos de inversão tomográfica. A necessidade de robustez na inversão pode ser analisada em seções sísmicas (FC, PMC) submetida a filtragens (freqüências de cantos: 5;15;75;85; banda-passante trapezoidal), onde se pode identificar, comparar e interpretar as informações contidas. A partir das seções, concluímos que os dados são contaminados com pontos isolados, o que propõe métodos na classe dos considerados robustos, tendo-se como referência a norma 2 (quadrados- mínimos) de ajuste de curvas. Os algoritmos foram desenvolvidos na linguagem de programação FORTRAN 90/95, usando o programa MATLAB para apresentação de resultados, e o sistema CWP/SU para modelagem sísmica sintética, marcação de eventos e apresentação de resultados.

Modelamento e correção de descentralização das imagens de tempo de trânsito

As imagens da parede do poço obtidas através da ferramenta UBI (Ultrasonic Borehole Imager) são amplamente utilizadas por geólogos e engenheiros de petróleo para identificar eventos geológicos nas paredes de poços abertos e na inspeção dos tubos de revestimento, uma vez que este perfil praticamente fornece uma fotografia da parede do poço. As ferramentas de imageamento acústico produzem imagens do tempo de trânsito e da amplitude do pulso acústico gerado pela ferramenta e refletido na parede do poço. Entretanto, estas imagens podem ter uma interpretação não realista, uma vez que elas podem ter seu aspecto alterado em razão do movimento da ferramenta no interior do poço. Este trabalho apresenta o modelamento das imagens de tempo de trânsito a partir da aplicação do critério de Coulomb para a ruptura da parede do poço submetida a um estado plano de tensões, a qual fornecerá a seção do poço, que é a forma geométrica que será traçada pela ferramenta de imageamento acústico do poço. O deslocamento ascensional da ferramenta e as imperfeições da parede do poço, normalmente, são os responsáveis pelo deslocamento da posição do transdutor em relação ao eixo do poço. Este efeito pode ter grande responsabilidade nas imperfeições das imagens de tempo de trânsito. Assim, a correção dessas imagens, chamada de correção da descentralização, busca o reposicionamento do transdutor para a posição do eixo do poço. Apresenta-se, também, um método de correção do efeito da descentralização da ferramenta baseado neste modelamento. O método é proposto com base na geometria analítica plana e no método do raio para a definição do tempo de trânsito do pulso acústico, com o objetivo de reconstruir as imagens de tempo de trânsito obtidas com a ferramenta descentralizada, ou seja, corrigir estas imagens tornando-as como se fossem obtidas com a ferramenta centralizada em relação ao eixo do poço.