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.

Tide distortion and attenuation in an Amazonian tidal river

O presente estudo teve como objetivo investigar a propagação da maré no sistema fluvial Guamá-Capim, na região amazônica, considerando aspectos hidrológicos e geomorfológicos. Os métodos empregados incluem dados históricos de vazões fluviais e níveis d'água, além de medições próprias de maré em diferentes locais e períodos ao longo do sistema. Os principais pontos defendidos no presente trabalho incluem a vazão fluvial como principal fator para a distorção da maré e consequente formação de pororoca no sistema, assim como o baixo relevo da área seria responsável por incremento na incursão da maré para o continente. Os resultados revelam uma penetração da maré de mais de 200 km, ocorrendo também uma forte deformação da maré, se intensificando gradualmente a montante, resultando em uma vazante até 5 horas mais longa que a enchente 161 km a montante, apresentando também velocidades de enchente levemente superiores, incluindo o decaimento contínuo da altura da maré, intensificado a partir desse ponto. Sazonalmente, a vazão fluvial aumenta 10 vezes no rio Guamá e 4 vezes no rio Capim. Durante períodos de alta descarga, combinados com marés equinociais (e.g. março-abril), o fenômeno da pororoca ocorre no sistema, em associação ao baixo relevo da área. Desta forma as principais conclusões são de que a forte deformação da maré relacionada a altas descargas fluviais e o baixo relevo da área são os fatores preponderantes na propagação da maré no sistema e formação de pororoca. Mais além, o sistema foi classificado como tidal river, onde a enorme descarga de água doce na região costeira amazônica resulta na prática ausência de salinidade no sistema estudado. Assim, seu estudo vem contribuir na definição e classificação de sistemas estuarinos.

Circulation and salt intrusion in the Piaçaguera Channel, Santos (SP)

A análise de dados termohalinos e correntes medidos em uma estação fixa no Canal de Piaçaguera (Estuário de Santos) no inverno foi feita em termos de condições cíclicas da maré (quadratura e sizígia) e quase-estacionária, com o objetivo de caracterizar a estratificação da massa de água estuarina, sua circulação e transporte de sal forçados pela modulação quinzenal da maré. Foram utilizados métodos clássicos de análise de dados observacionais horários e quase sinóticos e de simulações analíticas de perfis estacionários de salinidade e do componente longitudinal da velocidade. Durante o ciclo de maré de quadratura as velocidades de enchente (v<0) e vazante (v>0) variaram de -0.20 m/s a 0.30 m/s, associadas à pequena variação de salinidade entre a superfície e o fundo (26.4 psu a 30.7 psu). No ciclo de sizígia a velocidade aumentou de -0.40 m/s a 0.45 m/s, mas a estratificação de salinidade permaneceu praticamente a mesma. Os perfis estacionários teóricos de salinidade e de velocidade apresentaram boa concordância (Skill próximo a 1,0) quando comparados aos perfis observacionais. Durante a modulação quinzenal da maré não houve alteração na classificação do canal estuarino (tipo 2a-parcialmente misturado e fracamente estratificado), pois a taxa de aumento da energia potencial não foi suficiente para ocasionar a erosão da haloclina. Esses resultados, associados à alta estabilidade vertical (RiL >20) e ao número de Richardson estuarino (1,6), permitem as seguintes conclusões: i) o mecanismo que forçou a circulação e os processos de mistura foi principalmente o balanço da descarga fluvial com a maré, associado ao componente baroclínico da força de gradiente de pressão; ii) não houve variações nas principais características termohalinas e da circulação devido à modulação quinzenal da maré; e iii) os perfis quase estacionários de salinidade e da velocidade foram adequadamente simulados com um modelo analítico clássico.