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.

Inversão gravimétrica do relevo do embasamento usando regularização entrópica

Apresentamos um novo método de inversão gravimétrica para estimar o relevo descontínuo do embasamento de bacias sedimentares, cujos pacotes sedimentares podem ter contraste de densidade constante ou variando hiperbolicamente com a profundidade. O método combina a maximização da medida de entropia de ordem zero e a minimização da medida de entropia de primeira ordem do vetor de soluções (profundidades do embasamento). O modelointerpretativo consiste de um conjunto de primas retangulares verticais justapostos, com o contraste de densidade conhecido, cujas espessuras representam as profundidades do embasamento e são os parâmetros a serem determinados. A minimização da entropia de ordem um favorece soluções apresentando descontinuidades abruptas e a maximização da entropia de ordem zero é empregada apenas para evitar sua minimização excessiva. O método foi aplicado a dados sintéticos simulando: (i) bacias intracratônicas com o relevo do embasamento suave ou apresentando descontinuidades abruptas localizadas e (ii) bacias marginais falhadas. No caso de embasamento suave, a comparação dos resultados com aqueles obtidos através da suavidade global e da suavidade ponderada mostrou que ambos os métodos delinearam o embasamento. No caso de embasamentos falhados, seja em bacias intracratônicas ou marginais, as soluções obtidas com a regularização entrópica e a suavidade ponderada foram equivalentes, delineando o embasamento e as descontinuidades, o que não ocorreu com a suavidade global. A regularização entrópica, no entanto, não requereu informação a priori sobre a profundidade, diferentemente da suavidade ponderada. Tanto o método proposto como os métodos das suavidades global e ponderada foram aplicadas a quatro conjuntos de dados reais. O primeiro é proveniente da Ponte do Poema que está localizada no campus da Universidade Federal do Pará, Belém, Pará, permitiu certificar a aplicabilidade prática do método, uma vez que o relevo abaixo da ponte é acessível. O segundo, da porção norte do Steptoe Valley, Nevada, Estados Unidos. As soluções obtidas com a regularização entrópica e com a suavidade ponderada mostram um embasamento apresentando diversas descontinuidades verticais, o que não aconteceu com a suavidade global. O terceiro conjunto provém do graben de San Jacinto, situado na Califórnia, Estados Unidos, que levou a soluções em que o graben é assimétrico e apenas a suavidade ponderada apresentou solução com fundo achatado. O quarto conjunto de dados é oriundo do graben do Büyük Menderes, localizado na Turquia ocidental. As soluções do método proposto e da suavidade ponderada apresentaram grandes deslocamentos verticais, confirmados pela geologia, o que não aconteceu com a suavidade global.

Métodos geofísicos aplicados à prospecção de cobre na Serras de Carajás-Pará

Através de um convênio entre o Núcleo de Ciências Geofísicas e Geológicas da UFPa e a Companhia Rio Doce Geologia e Mineração S.A., um trabalho geofísico foi realizado na área denomina da ALVO 1, MM1, Serra dos Carajás, sul do Pará, Brasil. Foram aplicados os métodos magnético, polarização elétrica induzida e AFMAG, visando determinar a existência e a distribuição de sulfetos de cobre em subsuperfície. O trabalho foi realizado em uma malha de reconhecimento de 9 km² (ALVO 1) e em outra de detalhe com aproximadamente 0.63 km² (Corpo 2). O método AFMAG detectou vários condutores deslocados entre si por prováveis falhamentos. A orientação destes condutores, coincide com o "strike" dominante do xisto, o qual contém sulfetos de cobre associado a uma formação ferrífera. Pela estreita correlação que guardam com as anomalias magnéticas é possível supor-se que as anomalias de AFMAG são provocadas primeiramente pela formação ferrífera e secundariamente pelos sulfetos. O método magnético detectou fortes anomalias sobre a formação ferrífera. O deslocamento dos segmentos de anomalias implica em falhamento. Este método é muito importante na área, pela relação mantida entre a zona dos sulfetos e a formação ferrífera. O método polarização elétrica induzida/resistividade detectou anomalias com EPF variando de 6 a 12% aproximadamente. Várias destas anomalias foram confirmadas com furos de sondagem e trabalho de laboratório (Rocha, 1979). Verificou-se que essas anomalias provavelmente são provocadas por efeito combinado da formação ferrífera magnética e do sulfeto. A pesquisa geofísica realizada nesta área piloto tem grande importância porque os resultados obtidos poderão ser estendidos a uma área maior, em condições geológicas semelhantes.

Uso e gestão de águas subterrânea em Araguaína/TO

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

Mangrove sedimentary characteristics and implications for crab Ucides cordatus (Crustacea, Decapoda, Ucididae) distribution in an estuarine area of the Amazonian region

Na zona costeira, sedimentos, água e organismos interagem intensamente. Nas costas equatoriais dominadas por maré os manguezais são abundantes. Estas áreas são conhecidas por sua importância ecológica. No caso dos manguezais da costa atlântica da América do Sul o caranguejo-uçá Ucides cordatus (Linnaeus, 1763) se destaca por sua relevância ecológica e econômica, sendo que altas densidades deste organismo são encontradas na zona costeira amazônica. O presente estudo investiga a distribuição de sedimentos nos manguezais de Bragança (costa Amazônica, Brasil) e suas correlações com a vegetação e a distribuição do caranguejo-uçá. Quarenta e sete amostras de sedimento foram avaliadas, assim como caranguejos de 8 destas áreas foram coletados, onde o tipo dominante de vegetação foi também identificado. Os resultados demonstram que os sedimentos superficiais, assim como no extrato 0,8 a 1 m de profundidade, na área são principalmente compostos por silte (59%), incluindo em média 21% de areia e 20% de argila. O tipo de vegetação predominante foi significativamente correlacionado com a abundância e tamanho/peso dos caranguejos. As características sedimentares também foram substancialmente diferentes dependendo da vegetação. Áreas dominadas por Avicennia germinans tiveram mais areia e argila que as áreas dominadas por Rizophora mangle, onde a fração silte prevalece grandemente e os caranguejos eram significativamente maiores e mais abundantes. Os resultados demonstraram que sedimentos, invertebrados bentônicos e vegetação estão intimamente relacionados nos manguezais e devem ser estudados de maneira integrada.

Modelagem de dados MT 2D em multicamadas com anisotropia arbitrária

Neste trabalho apresentamos a modelagem direta bidimensional do método Magnetotelúrico (MT) com anisotropia arbitrária na condutividade elétrica para modelos de multicamadas estratificadas horizontalmente, separando a resposta do campo resultante em duas partes, primária e secundária. A parte primária refere-se ao campo 1D do modelo na ausência de heterogeneidades; a secundária, à contribuição no campo resultante gerada pela anomalia. Esta modelagem foi feita via técnica numérica dos Elementos Finitos (EF). Apresentamos também a modelagem direta do caso MT 1D para o mesmo modelo e anisotropia arbitrária, porém, com solução analítica pelo método da matriz de propagação; tal resultado é usado como fonte na formulação 2D. No estudo a respeito do comportamento dos campos 1D em profundidade na presença da anisotropia, verificamos nas componentes horizontais que suas curvas dentro de cada camada se assemelham a de um campo na presença de fontes, pois experimentam um aumento na amplitude; além de verificar o surgimento da componente normal E_z do campo, sem que haja uma densidade de corrente nesta direção. Com arbitrariedade na anisotropia queremos dizer que seus eixos principais podem assumir qualquer orientação em relação ao eixo de referência adotado para o problema; tal orientação é descrita sob três rotações consecutivas, chamadas de strike, dip e slant, que simulam inclinações bastante presentes na geologia dos ambientes terrestres. Verificamos que a anisotropia provoca alterações no comportamento das pseudo-seções de resistividades aparentes, interferindo na interpretação correta dos dados.

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.

Emissão do CO2 do solo em diferentes posições topográficas em área sob cultivo de cana-de-açúcar

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Variabilidade espacial da emissão de CO2 e sua relação com propriedades do solo em área de cana-de-açúcar no Sudeste do Brasil

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Distribuição espacial de Phytophthora nicotianae da cultura da cebola em amostras de solo e reação de cultivares ao patógeno

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

Distribuição espacial de Diaphorina citri Kuwayama (Hemiptera: Psyllidae) e sua relação com a expansão do greening em laranjeira Valência, utilizando geoestatística

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

Suscetibilidade magnética e sua associação com atributos do solo e do relevo em uma topossequência de Latossolos

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

A capa carbonática do sudoeste do cráton amazônico, estado de Rondônia: nova ocorrência e extensão dos eventos pós-glaciação marinoana (635 Ma)

Na porção oeste da Bacia dos Parecis, Estado de Rondônia, inserida no sudoeste do Cráton Amazônico, rochas carbonáticas expostas nas bordas dos grábens Pimenta Bueno e Colorado têm sido consideradas como parte do preenchimento eopaleozoico da bacia. A avaliação das fácies/microfácies e quimioestratigrafia dessas rochas nas regiões de Chupinguaia e Pimenta Bueno, confirmou a ocorrência de dolomitos rosados que sobrepõem, em contato direto, diamictitos glaciais previamente interpretados como depósitos de leques aluviais. Trabalhos prévios reportaram excursão negativa de δ¹³C, também confirmados neste trabalho, com variações entre -4.6 e -3,8‰VPDB em Chupinguaia e média de - 3,15‰VPDB em Pimenta Bueno. Esse padrão, de sedimentação e quimioestratigráfico, ausente nas rochas paleozoicas, é comumente encontrado nos depósitos carbonáticos anômalos do Neoproterozoico. No sul do Cráton Amazônico, Estado do Mato Grosso, rochas com essas mesmas características são descritas como capas carbonáticas relacionadas à glaciação marinoana (635 Ma). Neste trabalho, consideramos que os dolomitos rosados sobre diamictitos, em Rondônia, fazem parte do mesmo contexto das capas carbonáticas encontradas no Mato Grosso. Adicionalmente, destaca-se o contato brusco e deformado do dolomito sobre o diamictito, presente em ambas as ocorrências, configurando-se uma das feições típicas das capas carbonáticas do Cráton Amazônico. Essa relação paradoxal, entre diamictito e dolomito, tem sido interpretada como produto da mudança rápida das condições atmosféricas de icehouse para greenhouse, e a deformação da base foi gerada pelo rebound isostático. A capa carbonática de Rondônia compreende duas associações de fácies (AF2 e AF3) que recobrem depósitos glacio-marinhos compostos por paraconglomerados polimíticos (Pp), e arenito seixoso laminado (Asl), da AF1. A AF2 consiste em dolomudstone/dolopackstone peloidal com laminação plana a quasi-planar e com truncamentos de baixo-ângulo (fácies Dp), megamarcas onduladas (fácies Dm) e laminações truncadas por ondas (fácies Dt), interpretada como depósitos de plataforma rasa influenciada por ondas. Esta sucessão costeira é sucedida pela AF3, que compreende as fácies: dolomudstone/dolopackstone e dolomudstone/dolograinstone com partição de folhelho (Df) e siltito laminado (Sl). A fácies Df compreende um pacote de 6 metros de dolomito com partição de folhelho, apresentando lâminas de calcita fibrosa (pseudomorfos de evaporito) e dolomitos com laminações onduladas de corrente. Sobrejacente à fácies Df, ocorre a fácies Sl, apresentando 5 metros de siltito argiloso com laminação plana. Esta associação é interpretada como depósitos de plataforma rasa influenciada por maré, sendo sobreposta discordantemente, em contato angular, por depósitos glaciais do Eopaleozoico. Os valores isotópicos de C e O são negativos e refletem o sinal primário do C. No entanto, pode-se considerar uma leve influência da diagênese meteórica no sinal. As principais quebras nos sinais negativos podem estar associadas à influência meteórica, expressa pela substituição e preenchimento de poros por calcita e pela proximidade de superfícies estratigráficas, os quais refletem alguns padrões de alteração diagenética, representados nos sinais mais negativos. Diferentemente da capa carbonática do Mato Grosso, a capa de Rondônia possui níveis de pseudomorfos de evaporito e dolomitos com partição de folhelho (ritmito), em sucessão de fácies marinha rasa, onde os dolomitos de plataforma rasa influenciada por ondas passam para ritmitos e siltitos de plataforma rasa influenciada por maré (zona de inframaré), configurando uma sucessão retrogradante. Esta nova ocorrência de capa carbonática modifica a estratigrafia da base da Bacia dos Parecis, ao passo que exclui essas rochas carbonáticas da sequência eopaleozoica. Além disso, fornece informações que permitem reconstruir melhor a paleogeografia costeira da bacia neoproterozoica que acumulou os depósitos da plataforma carbonática do Grupo Araras, bem como estende os eventos pós-marinoanos ligados à hipótese do Snowball/Slushball Earth para o sudoeste do Cráton Amazônico, exposto no Estado de Rondônia.

Degradação ambiental da bacia hidrográfica do rio Uberaba - MG

Pós-graduação em Agronomia (Produção Vegetal) - FCAV

Desenvolvimento da vegetação e morfologia da foz do Amazonas-PA e rio Doce-ES durante o Quaternário tardio

Este trabalho compara as mudanças morfológicas e vegetacionais ocorridas ao longo da zona costeira da Ilha de Marajó, litoral amazônico, e da planície costeira do Rio Doce, sudeste do Brasil, durante o Holoceno e Pleistoceno tardio/Holoceno, respectivamente, com foco especificamente sobre a resposta dos manguezais para as flutuações do nível do mar e mudanças climáticas, já identificadas em vários estudos ao longo da costa brasileira. Esta abordagem integra datações por radiocarbono, descrição de características sedimentares, dados de pólen, e indicadores geoquímicos orgânicos (δ¹³C, δ¹₵N e C/N). Na planície costeira do Rio Doce entre ~47.500 e 29.400 anos cal AP, um sistema deltaico foi desenvolvido em resposta principalmente à diminuição do nível do mar. O aumento do nível do mar pós-glacial causou uma incursão marinha com invasão da zona costeira, favorecendo a evolução de um sistema estuarino/lagunar com planícies lamosas ocupadas por manguezais entre pelo menos ~7400 e ~5100 anos cal AP. Considerando a Ilha de Marajó durante o Holoceno inicial e médio (entre ~7500 e ~3200 anos cal AP) a área de manguezal aumentou nas planícies de maré lamosas com acúmulo de matéria orgânica estuarina/marinha. Provavelmente, isso foi resultado da incursão marinha causada pela elevação do nível do mar pós-glacial associada a uma subsidência tectônica da região. As condições de seca na região amazônica durante o Holoceneo inicial e médio provocou um aumento da salinidade no estuário, que contribuiu para a expansão do manguezal. Portanto, o efeito de subida do nível relativo do mar foi determinante para o estabelecimento dos manguezais na sua atual posição nas regiões norte e sudeste do Brasil. Entretanto, durante o Holoceno tardio (~3050-1880 anos cal AP) os manguezais em ambas as regiões retrairam para pequenas áreas, com algumas delas substituídas por vegetação de água doce. Isso foi causado pelo aumento da vazão dos rios associada a um período mais úmido registrado na região amazônica, enquanto que na planície costeira do Rio Doce, os manguezais encolheram em resposta a um aumento da entrada de sedimento fluvial associado a uma queda no nível relativo do mar.