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.

Processamento interpretativo de dados magnetométricos e inversão de dados gravimétricos aplicados à prospecção de hidrocarbonetos

Apresentamos dois métodos de interpretação de dados de campos potenciais, aplicados à prospecção de hidrocarbonetos. O primeiro emprega dados aeromagnéticos para estimar o limite, no plano horizontal, entre a crosta continental e a crosta oceânica. Este método baseia-se na existência de feições geológicas magnéticas exclusivas da crosta continental, de modo que as estimativas das extremidades destas feições são usadas como estimativas dos limites da crosta continental. Para tanto, o sinal da anomalia aeromagnética na região da plataforma, do talude e da elevação continental é amplificado através do operador de continuação analítica para baixo usando duas implementações: o princípio da camada equivalente e a condição de fronteira de Dirichlet. A maior carga computacional no cálculo do campo continuado para baixo reside na resolução de um sistema de equações lineares de grande porte. Este esforço computacional é minimizado através do processamento por janelas e do emprego do método do gradiente conjugado na resolução do sistema de equações. Como a operação de continuação para baixo é instável, estabilizamos a solução através do funcional estabilizador de primeira ordem de Tikhonov. Testes em dados aeromagnéticos sintéticos contaminados com ruído pseudo-aleatório Gaussiano mostraram a eficiência de ambas as implementações para realçar os finais das feições magnéticas exclusivas da crosta continental, permitindo o delineamento do limite desta com a crosta oceânica. Aplicamos a metodologia em suas duas implementações a dados aeromagnéticos reais de duas regiões da costa brasileira: Foz do Amazonas e Bacia do Jequitinhonha. O segundo método delineia, simultaneamente, a topografia do embasamento de uma bacia sedimentar e a geometria de estruturas salinas contidas no pacote sedimentar. Os modelos interpretativos consistem de um conjunto de prismas bidimensionais verticais justapostos, para o pacote sedimentar e de prismas bidimensionais com seções verticais poligonais para as estruturas salinas. Estabilizamos a solução, incorporando características geométricas do relevo do embasamento e das estruturas salinas compatíveis com o ambiente geológico através dos estabilizadores da suavidade global, suavidade ponderada e da concentração de massa ao longo de direções preferenciais, além de vínculos de desigualdade nos parâmetros. Aplicamos o método a dados gravimétricos sintéticos produzidos por fontes 2D simulando bacias sedimentares intracratônicas e marginais apresentando densidade do pacote sedimentar variando com a profundidade segundo uma lei hiperbólica e abrigando domos e almofadas salinas. Os resultados mostraram que o método apresenta potencial para delinear, simultaneamente, as geometrias tanto de almofadas e domos salinos, como de relevos descontínuos do embasamento. Aplicamos o método, também, a dados reais ao longo de dois perfis gravimétricos sobre as Bacias de Campos e do Jequitinhonha e obtivemos interpretações compatíveis com a geologia da área.

Aplicação dos métodos elétrico e eletricomagnético para o estudo hidrogeológico no município de Altamira-Pará

Este trabalho trata da aplicação dos métodos geofísicos eletroresistivo e eletromagnético no município de Altamira, região central do Estado do Pará, com a finalidade de fazer o mapeamento da espessura das camadas aqüíferas, bem como a ocorrência de intrusões de diabásio na área do município, permitindo uma melhor locação de pontos para perfuração com fins de captação à água subterrânea. Na área estudada, verificam-se afloramentos de rochas cristalinas do embasamento coberto por sedimentos com idade variando desde o Paleozóico ao Quaternário, com esperança de se obter água nos sedimentos do Devoniano Inferior, mais precisamente na formação Maecuru. Foram executadas 34 sondagens elétricas verticais com arranjo Schlumberger com AB/2 variando de 100 a 500 metros. 108 sondagens eletromagnéticas transientes também foram executadas com arranjo "loop coincidente", com 50 metros de lado. Por problemas eletrônicos adicionados aos ruídos de fenômenos supraparamagnéticos, os resultados das sondagens eletromagnéticas apresentaram-se muito distorcidos, não tendo sido possível extrair deles informação geológica de relevância. Por isso, foram criados modelos eletromagnéticos sintéticos representando a situação geológica da área. Os dados obtidos no campo em forma de curvas de resistividade aparente, foram plotados em papel bilogarítmo e interpretados com a ajuda de algoritmos direto e inverso os quais considera modelos de camadas horizontalmente estratificadas, homogêneas e isotrópicas. Os resultados da interpretação permitiram estabelecer modelos geoelétricos aproximados refletindo a situação geológica na subsuperfície. Os resultados da interpretação das sondagens elétricas, foram confrontados com os resultados da interpretação qualitativa do mapa de isogamas do campo magnético total, permitindo delimitar com maior precisão a zona de incidência de soleiras de diabásio, situada a norte da área. Dos quatro horizontes geoelétricos distinguidos, o terceiro é essencialmente constituído por arenitos, apresentando, por isso, melhores condições para armazenamento de água subterrânea.

Eletroresistividade na prospecção hidrogeológica no município de Ourilândia do Norte-PA

A prospecção geofísica para a exploração de água subterrânea em regiões cristalinas torna-se uma tarefa extremamente difícil pela presença de estruturas como falhas, fraturas e ocorrência de variação lateral de resistividade no manto intempérico. Para analisar a influência da variação bidimensional da propriedade física no comportamento das curvas de sondagem elétrica vertical Schlumberger, simulou-se modelos contendo tais estruturas encaixadas no manto de alteração com ângulo de mergulho de 30° e 90° e com a posição do centro de sondagem variando em relação a heterogeneidade. A modelagem numérica foi feita através do método dos elementos finitos com o programa EGSLIB/SEV2D. O resultado da simulação mostra que o emprego da técnica da sondagem elétrica vertical Schlumberger sobre meios bidimensionais não funciona eficientemente, como quando aplicada em regiões sedimentares e de aluvião, devido o fenômeno da ambigüidade que dificulta a interpretação das curvas. Dentre os métodos geofísicos para a investigação hidrogeológica, a aplicação dos métodos eletrorresistivos mostra-se mais simples, eficaz e econômica. Assim, o levantamento geofísico para orientar os trabalhos de perfurações de poços na cidade de Ourilândia do Norte contou de trinta e duas sondagens elétricas verticais tipo Schlumberger e um caminhamento elétrico dipolo-dipolo realizados em diversas ruas da cidade. As curvas de resistividade aparente, obtidas com as sondagens foram tratadas e processadas considerando um modelo de camadas planas, horizontais, homogêneas e isotropicas, inicialmente com o algoritmo EGSLIB/SEV1D e posteriormente com uso do programa EGSLIB/SEV1DINV, permitindo estabelecer modelos geoelétricos aproximados para a configuração da subsuperfície da cidade. Depois da interpretação quantitativa, foi possível separar duas zonas geoelétricas distintas. Esta diferenciação foi realizada com base nos valores de resistividade aparente e as informações geológicas de superfície. A primeira, com pouca representatividade na área da cidade, apresenta-se constituída por solo/aterro sobre saprolito de natureza argilo-arenosa recobrindo os granitóides tipo Rio Maria. A segunda composta pelo solo/aterro seguido de um horizonte preenchido por sedimentos arenosos superposto a camada argilo-arenosa e finalmente a rocha sã. Esta seqüência geoelétrica predomina na cidade. Os estratos arenoso e argilo-arenoso apresentam condições potenciais para armazenamento de água.

Interpretação gravimétrica-magnética da zona de fratura dupla Bode Verde

Em Janeiro de 1980, foi realizada a Operação Bode Verde I do Projeto CENTRATLAN, que constou da coleta de dados geofísicos e geológicos na área central da Cadeia Mesoatlântica (CMA), a sul da Ilha de Ascensão. Esta operação fez parte de um convênio firmado entre a Marinha Americana e a Marinha Brasileira, cujo objetivo era investigar o Oceano Atlântico Sul, especialmente a área onde se estende a província fisiográfica da Cadeia Mesoatlântica. Com o objetivo de estudar a Zona de Fratura Dupla Bode Verde, e assim conhecer melhor esta feição tão marcante no fundo oceânico e margem continental, foram tratados e interpretados os dados batimétricos, gravimétricos e magnéticos colhidos na operação anteriormente citada. A análise dos dados batimétricos permitiu uma definição para esta zona de fratura, com o seu caráter duplo marcado pela presença de duas calhas contínuas, paralelas entre si, e separadas por um alto que tem o seu próprio "rift-valley". As anomalias gravimétricas encontradas confirmam o caráter duplo, com os dois mínimos observados correspondendo às duas calhas da fratura. Estas calhas separam um bloco crustal de aproximadamente 40 km de largura. Com base nos dados gravimétricos, foram construídas nove seções crustais sobre a zona de fratura, que permitiram observar que a interface crosta-manto sofre um afinamento crustal embaixo das paredes da fratura. Esta interface, que encontra-se normalmente a uma profundidade de 8-9 km, passa a atingir profundidades de 5.5-6.0 km abaixo das paredes da fratura, com a crosta apresentando uma espessura de 2.0 km nestas partes. Este afinamento crustal é causado pela subida de material do manto. A extensiva alteração hidrotermal que ocorre na depressão central da fratura, pode ser a responsável pelo menor afinamento crustal observado nesta parte da zona de fratura. Para a interpretação dos dados magnéticos, foi usada uma escala de tempo de polaridade magnética entre o Cretáceo Inferior e o Cenozóico, e uma taxa de espalhamento oceânico de 2.0 cm/ano. Para a camada de basalto, que é responsável por parte das anomalias magnéticas observadas, foi usada uma espessura média de 0.5 km. Nas paredes da zona de fratura, ocorre uma diminuição na espessura desta camada, havendo uma interrupção da mesma na depressão central da fratura.

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.

Geology, petrology, U-Pb (shrimp) geochronology of the Morrinhos granite -Paraguá terrane, SW Amazonian craton: implications for the magmatic evolution of the San Ignácio orogeny

O Granito Morrinhos é um corpo batolítico levemente alongado segundo a direção NNW, de aproximadamente 1.140 km², localizado no município de Vila Bela da Santíssima Trindade, estado de Mato Grosso. Situa-se no Terreno Paraguá, Província Rondoniana-San Ignácio, na porção SW do Cráton Amazônico. Essa intrusão exibe uma variação composicional entre tonalito a monzogranito, textura inequigranular média a grossa, localmente, porfirítica, tendo biotita como máfico predominante em uma das fácies e hornblenda na outra, ambas metamorfizadas na fácies xisto verde. As rochas estudadas caracterizam uma sequência intermediária a ácida formada por um magmatismo subalcali no, do tipo álcali-cálcico, metaluminoso a levemente peraluminoso evoluído por meio de mecanismos de cristalização fracionada. Dados estruturais exibem registros de duas fases deformacionais, representadas pela foliação penetrativa (S₁) e dobras abertas (D₂) ambas, provavelmente, relacionadas à Orogenia San Ignácio. A investigação geocronológica (U-Pb SHRIMP) e geoquímica isotópica (Sm-Nd) dessas rochas indicaram, respectivamente, idade de cristalização 1.350 ± 12 Ma, T_DM em torno de 1,77 Ga e valor negativo para εNd(_1,35) de -2,57, sugerindo uma geração relacionada com processo de fusão parcial de uma crosta continental paleoproterozoica (estateriana). Os resultados aqui obtidos indicam que o Granito Morrinhos foi gerado em arco magmático continental, em estágio tardi a pós-orogênico, da Orogenia San Ignácio e permite reconhecê-lo como pertencente à Suíte Intrusiva Pensamiento.

Geology, petrology and geochronology (Pb-Pb) of the Serra da Bocaina Formation: evidence of an Orosirian Amoguijá Magmatic Arc in the Rio Apa Terrane, south of the Amazonian Craton

Neste trabalho apresentam-se os resultados do mapeamento geológico e caracterização petrológica da Formação Serra da Bocaina, pertencente ao Arco Magmático Amoguijá do Terreno Rio Apa, sul do Cráton Amazônico. A Formação Serra da Bocaina, na serra da homônima, consiste de rochas vulcânicas paleoproterozoicas de composição intermediária a predominantemente ácida, classificadas como andesito e riolitos, subdivididas em cinco fácies petrográficas sendo quatro piroclásticas e uma efusiva, que mantêm contato tectônico, a leste, com o Granito Carandá. Nas rochas estudadas estruturas tectônicas são formadas em duas fases deformacionais compressivas de natureza dúctil e dúctil-rúptil, respectivamente. A primeira fase, mais intensa, é observada ao longo de toda a área estudada e é responsável pela Zona de Cisalhamento Santa Rosa enquanto a segunda fase é mais discreta e localizada. O tratamento geoquímico indica que essas rochas foram geradas num ambiente de arco-vulcânico a partir de um magmatismo calcioalcalino de médio a alto-K, peraluminoso. Estas rochas retratam um evento magmático extrusivo, de natureza explosiva, relacionado à evolução do Arco Magmático Amoguijá, conforme resultado Pb-Pb em zircão de 1877,3 ± 3,9 Ma., interpretada como idade de cristalização destas rochas.

Geology, Geochemistry, and Geochronology (U-Pb) of the Rio Fortuna Gneiss - Serra do Baú Intrusive Suite - Paraguá Terrane - SW Amazonian Craton

O Gnaisse Rio Fortuna aflora na região da serra Santa Bárbara, nas imediações do Destacamento Militar Fortuna, na fronteira Brasil-Bolívia. Estes ortognaisses estão inseridos no Terreno Paraguá, em um setor afetado pela Orogenia Sunsás (1.0 a 0.9 Ga.). São classificados como ortognaisses de composição monzo a granodiorítica, com registros de, no mínimo, três fases de deformação. Idade U-Pb em zircão de 1.711 ± 13 Ma obtida por ablação a laser MC-ICP-MS, é considerada como correspondendo à idade de cristalização do protólito ígneo. Geoquimicamente, essas rochas constituem uma sequencia ácida formada por um magmatismo subalcalino, do tipo cálcio-alcalino de alto potássio, metaluminoso a peraluminoso.