Avaliação da influência do eletrojato equatorial na interpretação de dados magnetotelúrico bidimensionais

O campo magnetotelúrico em regiões equatoriais viola a condição de ondas planas por causa de uma fonte fortemente concentrada na direção E-W na ionosfera, denominada eletrojato equatorial. No presente trabalho, procurou-se analisar a resposta magnetotelúrica de fontes que simulam o efeito do eletrojato equatorial. Foram considerados dois tipos de fontes para simular o eletrojato: uma linha infinita de corrente e uma distribuição gaussiana de densidade de corrente em relação a uma das coordenadas horizontais. A resistividade aparente foi obtida através da relação de Cagniard e comparada com os resultados de ondas planas. É mostrada também a comparação entre a fase da impedância na superfície, para os três tipos de fontes (ondas planas, eletrojato gaussiano e linha de corrente). O problema de meios com heterogeneidades laterais foi resolvido em termos de campos secundários, sendo as equações diferenciais solucionadas através da técnica de elementos finitos bidimensionais. Os resultados mostram que o eletrojato tem pouca influência nas respostas (resistividade aparente e fase) de estruturas geológicas rasas. Entretanto, a influência pode ser considerável nas estruturas profundas (maior que 5000 m), principalmente se suas resistividades são altas (maior que 100 Ω.m). Portanto, a influência do eletrojato equatorial deve ser considerada na interpretação de dados magnetotelúricos de bacias sedimentares profundas ou no estudo da crosta terrestre.

Comparação dos filtros de velocidade e do operador WHLP-CRS na atenuação de múltiplas

A motivação geológica deste trabalho reside no imageamento de estruturas de bacias sedimentares da região Amazônica, onde a geração e o acúmulo de hidrocarboneto estão relacionados com a presença de soleiras de diabásio. A motivação sísmica reside no fato de que essas rochas intrusivas possuem grandes contrastes de impedância com a rocha encaixante, o que resulta em múltiplas, externas e internas, com amplitudes semelhantes as das primárias. O sinal sísmico das múltiplas podem predominar sobre o sinal das reflexões primárias oriundas de interfaces mais profundas, o que pode dificultar o processamento, a interpretação e o imageamento da seção sísmica temporal. Neste trabalho, estudamos a atenuação de múltiplas em seções sintéticas fonte-comum (FC) através da comparação de dois métodos. O primeiro método resulta da combinação das técnicas Wiener-Hopf-Levinson de predição (WHLP) e o de empilhamento superfície-de-reflexão-comum (CRS), e denominando WHLP-CRS, onde o operador é desenhado exclusivamente no domínio do tempo-espaço. O segundo método utilizado é o filtro de velocidade (ω-k) aplicado após o empilhamento superfície-de-reflexão (CRS), onde o operador é desenhado exclusivamente no domínio bidimensional de freqüência temporal-espacial. A identificação das múltiplas é feita na seção de afastamento-nulo (AN) simulada com o empilhamento CRS, e utiliza o critério da periodicidade entre primária e suas múltiplas. Os atributos da frente de onda, obtidos através do empilhamento CRS, são utilizados na definição de janelas móveis no domínio tempo-espaço, que são usadas para calcular o operador WHLP-CRS. O cálculo do filtroω-k é realizado no domínio da freqüência temporal-espacial, onde os eventos são selecionados para corte ou passagem. O filtro (ω-k) é classificado como filtro de corte, com alteração de amplitude, mas não de fase, e limites práticos são impostos pela amostragem tempo-espaço. Em termos práticos, concluímos que, para o caso de múltiplas, os eventos separados no domínio x-t não necessariamente se separam no domínio ω-k, o que dificulta o desenho de um operador ω-k semelhante em performance ao operador x-t.

Os efeitos Delaware e Groningen: um estudo quantitativo por elementos finitos

Os efeitos Delaware e Groningen são dois tipos de anomalia que afetam ferramentas de eletrodos para perfilagem de resistividade. Ambos os efeitos ocorrem quando há uma camada muito resistiva, como anidrita ou halita, acima do(s) reservatório(s), produzindo um gradiente de resistividade muito similar ao produzido por um contato óleo-água. Os erros de interpretação produzidos têm ocasionado prejuízos consideráveis à indústria de petróleo. A PETROBRÁS, em particular, tem enfrentado problemas ocasionados pelo efeito Groningen sobre perfis obtidos em bacias paleozóicas da região norte do Brasil. Neste trabalho adaptamos, com avanços, uma metodologia desenvolvida por LOVELL (1990), baseada na equação de Helmholtz para H_Φ, para modelagem dos efeitos Delaware e Groningen. Solucionamos esta equação por elementos finitos triangulares e retangulares. O sistema linear gerado pelo método de elementos finitos é resolvido por gradiente bi-conjugado pré-condicionado, sendo este pré-condicionador obtido por decomposição LU (Low Up) da matriz de stiffness. As voltagens são calculadas por um algoritmo, mais preciso, recentemente desenvolvido. Os perfis são gerados por um novo algoritmo envolvendo uma sucessiva troca de resistividade de subdomínios. Este procedimento permite obter cada nova matriz de stiffness a partir da anterior pelo cálculo, muito mais rápido, da variação dessa matriz. Este método permite ainda, acelerar a solução iterativa pelo uso da solução na posição anterior da ferramenta. Finalmente geramos perfis sintéticos afetados por cada um dos efeitos para um modelo da ferramenta Dual Laterolog.

Interpretação de perfis dos carbonatos fraturados da Bacia do Pará-Maranhão

A finalidade deste trabalho é apresentar um estudo de caso dos carbonatos terciários da bacia do Pará-Maranhão, do ponto de vista da interpretação dos perfis registrados nessa área. Dois poços-chave, X e Y, foram escolhidos para o estudo. O tratamento dos dados foi realizado utilizando os recursos do programa LOGCALC, instalado no computador IBM-3090, da Petrobrás, e, também, o sistema DLPS, instalado no VAX-8600, da Universidade Federal de Pará. A avaliação da porosidade e, principalmente, das saturações, é dificultada pelas características não-convencionais dos carbonatos. A litologia é complexa, a salinidade da água de formação é relativamente baixa, em torno de 10.000 ppm de NaCl, a densidade das rochas é elevada e os carbonatos estão fraturados. Para diminuir o efeito da composição mineralógica, foi necessário discriminar os diferentes tipos litológicos. Três tipos de carbonatos foram identificados: calcário, calcário arenoso e calcário dolomitizado. A identificação litológica permite maior controle dos parâmetros da matriz e dos expoentes de porosidade das rochas, conduzindo a estimativas de porosidade e de saturações mais confiáveis. A presença de fraturas influencia, marcadamente, a resposta dos perfis, conforme pode ser notado no perfil de identificação de fraturas, nos perfis de resistividade, de densidade e na curva de raios-gama espectral do poço X. O expoente de porosidade, m, tomado do gráfico de Pickett, é frequentemente inferior a 1,5, valor considerado inerente a rochas fraturadas. Os modelos de Rasmus (1983) e de Porter et al.(1969), foram testados para calcular as saturações de água e óleo. A equação tradicional de Archie (1942), com a e m apropriados, também foi aplicada neste trabalho. O esquema para estimativa de saturações que conduz a resultados mais coerentes com os dados dos testes de formação, nos intervalos fraturados, é a saturação de Rasmus, calculada com a relação de Archie aplicada com o expoente de porosidade variável, tomado do modelo de Rasmus.

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 estável do relevo da base e da variação da densidade com a profundidade em aterros sanitários

O presente trabalho apresenta um método que utiliza uma função hiperbólica para descrever o decaimento do contraste de densidade com a profundidade em aterros sanitários e lixões, e consiste de duas abordagens: (1) o mapeamento da base de um lixão ou de um aterro sanitário conhecendo-se a anomalia gravimétrica, o contraste de densidade na superfície e o fator de decaimento do contraste de densidade com a profundidade; (2) estimação do contraste de densidade na superfície e o fator de decaimento do contraste de densidade com a profundidade a partir do conhecimento da anomalia gravimétrica e do relevo da base do lixão ou do aterro sanitário; uma variante desta última abordagem utiliza os valores conhecidos de profundidade desses ambientes para gerar a anomalia ajustada e estimar o valor do contraste de densidade na superfície e o fator de decaimento do contraste com a profundidade através da minimização da norma Euclidiana do vetor de resíduos entre as observações e o ajuste. Em ambas abordagens o modelo interpretativo é representado por um conjunto de prismas retangulares justapostos, cujas espessuras são os parâmetros a serem determinados. As soluções encontradas através das duas abordagens são estáveis devido à incorporação de informação sobre a suavidade do relevo estimado. O método foi aplicado a dados sintéticos produzidos pela simulação de aterros sanitários e lixões com até 15 m de profundidades e de relevos suaves. Os resultados das interpretações gravimétricas foram consistentes e revelaram a eficácia do método tanto na abordagem (1) como na abordagem (2). O método também foi aplicado à anomalia Bouguer residual do aterro sanitário da Fazenda Thomas, Indiana, USA, indicando que o contraste de densidade do aterro pode ou ser constante ou apresentar variação máxima de 0,09 g/cm³.

Inversão gravimétrica rápida do relevo do embasamento aplicando o funcional da variação total

Apresentamos um novo método de inversão linear bidimensional de dados gravimétricos produzidos por bacias sedimentares com relevo do embasamento descontínuo. O método desenvolvido utiliza um modelo interpretativo formado por um conjunto de fitas horizontais bidimensionais justapostas cujas espessuras são os parâmetros a serem estimados. O contraste de densidade entre o embasamento e os sedimentos é presumido constante e conhecido. As estimativas das espessuras foram estabilizadas com o funcional da Variação Total (VT) que permite soluções apresentando descontinuidades locais no relevo do embasamento. As estimativas do relevo são obtidas através da resolução de um sistema de equações lineares, resolvido na norma L₁. Como métodos lineares subestimam as estimativas de profundidade do embasamento de bacias maiores que cerca de 500 m, amplificamos as estimativas de profundidade através da modificação da matriz associada ao modelo interpretativo de fitas. As estimativas obtidas através deste procedimento são em geral ligeiramente superestimadas. Desse modo, elas são corrigidas através de uma correção definida pela expressão da placa Bouguer. Testes em dados sintéticos e reais produziram resultados comparáveis aos produzidos pelo método não linear, mas exigiram menor tempo computacional. A razão R entre os tempos exigidos pelo método não linear e o método proposto cresce com o número de observações e parâmetros. Por exemplo, para 60 observações e 60 parâmetros, R é igual a 4, enquanto para 2500 observações e 2500 parâmetros R cresce para 16,8. O método proposto e o método de inversão não linear foram aplicados também em dados reais do Steptoe Valley, Nevada, Estados Unidos, e da ponte do POEMA, no Campus do Guamá em Belém, produzindo soluções similares às obtidas com o método não linear exigindo menor tempo computacional.

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.

Investigação por eletroresistividade de ambientes cársticos no município de São João de Pirabas - PA

A presente tese procura investigar ambientes cársticos, usando os métodos eletroresistivos de sondagem elétrica vertical e caminhamento dipolo-dipolo. A área piloto para essa investigação está localizada na porção nordeste do Estado do Pará, no Município de São João de Pirabas. Nessa localidade, durante a perfuração de um poço para abastecimento de água para a cidade, a equipe de prospecção da Fundação Nacional de Saúde encontrou problemas devido a presença de cavidades nos calcários da Formação Pirabas. Para analisar esse problema foram realizadas 15 sondagens elétricas verticais, usando o arranjo Schlumberger, em diversos pontos da cidade. A abertura máxima dos eletrodos de corrente, para a realização das sondagens, foi em torno de 600 metros. As curvas de resistividade aparente das sondagens foram, inicialmente, interpretadas usando programas computacionais que calculam as resistividade do meio investigado com modelos unidimensionais, sem variações laterais de resistividade. Porém, as sondagens realizadas próximas ao poço, onde ocorrem cavernas, foram interpretadas usando modelos bidimensionais para representar as variações laterais de resistividade devido às cavernas. O algoritmo usado para essa interpretação se baseia na técnica dos elementos finitos. Os resultados das interpretações das sondagens elétricas, para os vários modelos de cavernas estudados, mostraram que as anomalias de resistividade devido a presença de cavernas são muito pequenas. Devido a este fato, o método de sondagem elétrica vertical não apresentou boa resolução para identificar com precisão a presença deste tipo de caverna em subsuperfície. Por outro lado, os resultados numéricos obtidos pelo caminhamento dipolo-dipolo mostraram que esse método é mais eficiente e menos ambíguo, para estudar cavernas, que o método de sondagem elétrica vertical. Contudo o método de caminhamento dipolo-dipolo é muito mais dispendioso do que o de sondagem elétrica verticais, tanto do ponto de vista de trabalho de campo como no processamento e interpretação de dados.

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).

Avaliação quantitativa da dinâmica espaço-temporal da precipitação na região hidrográfica Tocantins-Araguaia

A análise do comportamento da precipitação em uma bacia hidrográfica é fundamental para a engenharia e gerenciamento dos recursos hídricos. A Região Hidrográfica Tocantins-Araguaia (RHTA) pela sua ocupação recente e potencialidades econômicas, ganha destaque no cenário nacional. Este trabalho avalia quantitativamente a dinâmica espaço-temporal da precipitação anual nesta região durante um período de 30 anos de dados. A dinâmica da precipitação pode ser analisada pelo cálculo da precipitação média em uma dada área, compondo assim mapas de isoietas de precipitação anual. No entanto, a confecção destes mapas requer um método de interpolação que melhor represente as características pluviométricas em locais não amostrados para posterior análise quantitativa do comportamento da precipitação. Para tanto, foram realizados análises exploratórias descritivas amostral e espacial como requisito de estacionaridade do método de interpolação geoestatístico, ajuste e validação do modelo teórico que se adéque ao variograma de precipitação anual. Após a confecção do mapa de isoietas pelo método de Krigagem Ordinária (sem tendência) e Krigagem Universal (com tendência) foi realizado o cálculo do volume precipitado na região hidrográfica pelo método dos contornos. A dinâmica espacial da precipitação foi realizada com base na análise da estatística descritiva, mapa de isoietas, mapa hipsométrico, Índice de Irregularidade Meteorológica (IMM) e Coeficiente de Variação. A dinâmica temporal foi analisada pela distribuição dos totais anuais de precipitação volumétrica para cada sub-bacia da RHTA, pelo Índice de Anomalia Padronizada, na variação interanual de precipitação e teste de tendência e magnitude representados respectivamente pelo Teste de Mann Kendall e Sen’s. Os resultados correlacionados com as anomalias meteorológicas do Oceano Atlântico (Dipolo) e Pacífico (ENOS) indicam o comportamento da precipitação bastante heterogêneo e com grande variabilidade temporal principalmente na sub-bacia Tocantins-Alto (TOA) (14%). Diminuição da amplitude pluviométrica, em anos de anomalia meteorológica intensa ocasionando um incremento de precipitação ao sul das sub-bacias TOA e ARA e diminuição da precipitação na sub-bacia TOB, em eventos de El Niño. Não se pode comprovar pelo teste de Mann Kendall que há uma tendência estatisticamente significativa no volume precipitado na RHTA, mas o estimador Sen’s dá indícios de queda na precipitação na sub-bacia TOA (-1,24 Km^³/ano) e Araguaia (ARA) (-1,13 Km^³/ano) e aumento da precipitação na sub-bacia do Tocantins Baixo (TOB) (0,53 Km^³/ano) e para a RHTA (-1,5 Km^³/ano). Assim a variabilidade espacial e temporal nas sub-bacias está intimamente relacionada aos eventos de anomalia meteorológica, na qual, a sua ação ocorre de maneira irregular ao longo da área de estudo e pode influenciar as diversas atividades sócio-econômicas na RHTA de acordo com sua magnitude e área de ocorrência.

Paleoambiente e proveniência da formação cabeças da bacia do Parnaíba: evidências da glaciação famenniana e implicações na potencialidade do reservatório

O histórico de prospecção de hidrocarbonetos da Bacia Paleozoica do Parnaíba, situada no norte-nordeste do Brasil, sempre foi considerado desfavorável quando comparado aos super-reservatórios estimados do Pré-Sal das bacias da Margem Atlântica e até mesmo interiores, como a Bacia do Solimões. No entanto, a descoberta de gás natural em depósitos da superseqüência mesodevoniana-eocarbonífera do Grupo Canindé, que incluem as formações Pimenteiras, Cabeças e Longá, impulsionou novas pesquisas no intuito de refinar a caracterização paleoambiental, paleogeográfica, bem como, entender o sistema petrolífero, os possíveis plays e a potencialidade do reservatório Cabeças. A avaliação faciológica e estratigráfica com ênfase no registro da tectônica glacial, em combinação com a geocronologia de zircão detrítico permitiu interpretar o paleoambiente e a proveniência do reservatório Cabeças. Seis associações de fácies agrupadas em sucessões aflorantes, com espessura máxima de até 60m registram a evolução de um sistema deltaico Devoniano influenciado por processos glaciais principalmente no topo da unidade. 1) frente deltaica distal, composta por argilito maciço, conglomerado maciço, arenito com acamamento maciço, laminação plana e estratificação cruzada sigmoidal 2) frente deltaica proximal, representada pelas fácies arenito maciço, arenito com laminação plana, arenito com estratificação cruzada sigmoidal e conglomerado maciço; 3) planície deltaica, representada pelas fácies argilito laminado, arenito maciço, arenito com estratificação cruzada acanalada e conglomerado maciço; 4) shoreface glacial, composta pelas fácies arenito com marcas onduladas e arenito com estratificação cruzada hummocky; 5) depósitos subglaciais, que englobam as fácies diamictito maciço, diamictito com pods de arenito e brecha intraformacional; e 6) frente deltaica de degelo, constituída pelas fácies arenito maciço, arenito deformado, arenito com laminação plana, arenito com laminação cruzada cavalgante e arenito com estratificação cruzada sigmoidal. Durante o Fammeniano (374-359 Ma) uma frente deltaica dominada por processos fluviais progradava para NW (borda leste) e para NE (borda oeste) sobre uma plataforma influenciada por ondas de tempestade (Formação Pimenteiras). Na borda leste da bacia, o padrão de paleocorrente e o espectro de idades U-Pb em zircão detrítico indicam que o delta Cabeças foi alimentado por áreas fonte situadas a sudeste da Bacia do Parnaíba, provavelmente da Província Borborema. Grãos de zircão com idade mesoproterozóica (~ 1.039 – 1.009 Ma) e neoproterozóica (~ 654 Ma) são os mais populosos ao contrário dos grãos com idade arqueana (~ 2.508 – 2.678 Ma) e paleoproterozóica (~ 2.054 – 1.992 Ma). O grão de zircão concordante mais novo forneceu idade ²⁰⁶Pb/²³⁸U de 501,20 ± 6,35 Ma (95% concordante) indicando idades de áreas-fonte cambrianas. As principais fontes de sedimentos do delta Cabeças na borda leste são produto de rochas do Domínio Zona Transversal e de plútons Brasilianos encontrados no embasamento a sudeste da Bacia do Parnaíba, com pequena contribuição de sedimentos oriundos de rochas do Domínio Ceará Central e da porção ocidental do Domínio Rio Grande do Norte. No Famenniano, a movimentação do supercontinente Gondwana para o polo sul culminou na implantação de condições glaciais concomitantemente com o rebaixamento do nível do mar e exposição da região costeira. O avanço das geleiras sobre o embasamento e depósitos deltaicos gerou erosão, deposição de diamictons com clastos exóticos e facetados, além de estruturas glaciotectônicas tais como plano de descolamento, foliação, boudins, dobras, duplex, falhas e fraturas que refletem um cisalhamento tangencial em regime rúptil-dúctil. O substrato apresentava-se inconsolidado e saturados em água com temperatura levemente abaixo do ponto de fusão do gelo (permafrost quente). Corpos podiformes de arenito imersos em corpos lenticulares de diamicton foram formados pela ruptura de camadas pelo cisalhamento subglacial. Lentes de conglomerados esporádicas (dump structures) nos depósitos de shoreface sugere queda de detritos ligados a icebergs em fases de recuo da geleira. A elevação da temperatura no final do Famenniano reflete a rotação destral do Gondwana e migração do polo sul da porção ocidental da América do Sul e para o oeste da África. Esta nova configuração paleogeográfica posicionou a Bacia do Parnaíba em regiões subtropicais iniciando o recuo de geleiras e a influência do rebound isostático. O alívio de pressão é indicado pela geração de sills e diques clásticos, estruturas ball-and-pillow, rompimento de camadas e brechas. Falhas de cavalgamento associadas à diamictitos com foliação na borda oeste da bacia sugerem que as geleiras migravam para NNE. O contínuo aumento do nível do mar relativo propiciou a instalação de sedimentação deltaica durante o degelo e posteriormente a implantação de uma plataforma transgressiva (Formação Longá). Diamictitos interdigitados com depósitos de frente deltaica na porção superior da Formação Cabeças correspondem a intervalos com baixo volume de poros e podem representar trapas estratigráficas secundárias no reservatório. As anisotropias primárias subglaciais do topo da sucessão Cabeças, em ambas as bordas da Bacia do Parnaíba, estende a influência glacial e abre uma nova perspectiva sobre a potencialidade efetiva do reservatório Cabeças do sistema petrolífero Mesodevoniano-Eocarbonífero da referida bacia.

A zona de contato entre as formações Motuca e Sambaíba, Permo-Triássico da bacia do Parnaíba, regiões de Filadélfia (TO), Riachão (MA) e Loreto (MA)

O intervalo que compreende o final do Paleozóico e início do Mesozóico foi marcado por mudanças globais paleogeográficas e paleoclimáticas, em parte atribuídas a eventos catastróficos. A intensa continentalização do supercontinente Pangéia, com a implantação de extensos desertos, sucedeu os ambientes costeiros-plataformais do início do Permiano. Os registros desses eventos no norte do Brasil são encontrados nas bacias intracratônicas, particularmente na Bacia do Parnaíba, representados pela zona de contato entre as formações Motuca e Sambaíba. A Formação Motuca é constituída predominantemente por pelitos vermelhos laminados com lentes de gipsita, calcita e marga. Na porção leste da Bacia do Parnaíba, as fácies da Formação Motuca tornam-se mais arenosas com a ocorrência expressiva de arenitos com estratificação cruzada sigmoidal. A Formação Sambaíba consiste em arenitos de coloração creme alaranjada com estratificação plano-paralela e estratificação cruzada de médio a grande porte. Em geral, o contato entre as unidades é brusco, representado pela passagem de arenitos finos com laminação cruzada cavalgante e acamamento flaser/wavy da Formação Motuca para arenitos médios com falhas/microfalhas sinsedimentares e laminações convolutas da Formação Sambaíba. Foram individualizadas 14 fácies sedimentares, agrupadas em quatro associações: AF1 – Lacustre raso / Planície de lama (mudflat), AF2 – “Panela” salina (saline pan), AF3 – Lençol de areia e AF4 – Campo de dunas. A AF1 foi depositada dominantemente por processos de decantação em um extenso ambiente lacustre raso de baixa energia, influenciado por influxos esporádicos de areias oriundos de rios efêmeros. Este sistema lacustre foi, provavelmente, influenciado por períodos de contração e expansão, devido às variações das condições climáticas predominantemente áridas. Os mais expressivos períodos de contração ocorreram na porção oeste da Bacia do Parnaíba, representados pelo desenvolvimento de planícies de lama (mudflats) associadas a lagoas efêmeras saturadas em carbonatos e a “panelas” salinas (saline pans- AF2). Os lençóis de areia (AF3) são planícies arenosas extensas, localmente com área úmidas, intensamente retrabalhadas por processos eólicos. A AF4 é interpretada como parte de um erg composto por dunas/draas em zona saturada em areia, com interdunas secas subordinadas. Intervalos deformados lateralmente contínuos por centenas de quilômetros ocorrem na zona de contato entre as formações Motuca e Sambaíba. Pelitos com camadas contorcidas e brechadas (Formação Motuca) e arenitos com falhas/microfalhas sinsedimentares, laminação convoluta e diques de injeção preenchidos por argilitos (Formação Sambaíba) são interpretados como sismitos induzidos por terremotos de alta magnitude (>8 na escala Ritcher). Anomalias geoquímicas de elementos traços como Mn, Cr, Co, Cu e Ni na zona de contato entre as formações, juntamente com a presença de micropartículas de composição metálica na matriz argilosa dos sismitos, corroboram com impactos de meteoritos no limite c, possivelmente do astroblema Riachão.

Evolução geológica das sequências do embasamento do Cinturão Araguaia na região de Paraíso do Tocantins (TO), Brasil

Estudos isotópicos baseados nas metodologias Pb-Pb em zircão e Sm/Nd (rocha total), permitiram avanços no entendimento do quadro geológico evolutivo e litoestratigráfico do embasamento do segmento sul do Cinturão Araguaia. Os processos geológicos identificados aconteceram a partir do Arqueano (2,6 Ga e T_DM 2,78 – 3,25 Ga) e se estenderam até o Neoproterozoico. Os ortognaisses do Complexo Rio_dos Mangues posicionam-se no Paleoproterozoico (2,05 – 2,08 Ga) e T_DM 2,35 – 2,21 Ga. Um forte encurtamento crustal e fusão parcial de compartimentos isolados e espessados, gerou corpos ígneos de 1,85 e 1,82 Ga e o Granito Serrote (1,86 Ga), que provêm de fontes entre 2,50 e 2,43 Ga. No final do Mesoproterozoico a região foi marcada por processos tafrogenéticos, evidenciados por magmatismo máfico e alcalino (1,05 Ga) e bacias deposicionais, como a que acolheu os sedimentos que originaram as supracrustais do Cinturão Araguaia. No Neoproterozoico, através da inversão nas condições geodinâmicas, ocorreu novo processo de encurtamento/espessamento crustal com fusões que originaram expressivas massas batolíticas (Granitos Matança e Santa Luzia). O Cinturão Araguaia foi edificado a partir dessa movimentação tectônica. O transporte de massas tectônicas no sentido do Cráton Amazônico teria ocorrido, resultando na atual arquitetura em que se encontram as várias unidades lito-estratigráficas, organizadas sob a forma de lascas imbricadas.

Reconstituição paleoambiental das formações Motuca e Sambaíba, Permo-Triássico da Bacia do Parnaíba no sudoeste do Estado do Maranhão, Brasil

O intervalo que compreende o final do Paleozoico e início do Mesozoico foi marcado por mudanças globais paleogeográficas e paleoclimáticas, em parte atribuídas a eventos catastróficos. A intensa continentalização do supercontinente Pangeia, com a implantação de extensos desertos, sucedeu os ambientes costeiro-plataformais do início do Permiano. Os registros desses eventos no norte do Brasil são encontrados nas bacias intracratônicas, particularmente na sucessão Permotriássica da Bacia do Parnaíba. A análise de fácies e estratigráfica de afloramentos desta sucessão permitiu a individualização de 14 fácies sedimentares agrupadas em 4 associações de fácies (AF): AF1 e AF2, relacionadas aos depósitos da Formação Motuca, e AF3 e AF4, representativas da base da Formação Sambaíba. A AF1 - Lacustre raso/Mudflat consiste em pelitos vermelhos laminados com lentes de gipsita, calcita e marga, além de lobos de arenitos sigmoidais. A AF2 - Saline pan é constituída por corpos lenticulares de gipso laminado, gipso nodular e gipsarenito, sobrepostos por pelitos esverdeados com nódulos de dolomita e palygorskita. A AF3 - Lençol de areia e AF4 - Campo de dunas são formadas, respectivamente, por arenitos de coloração creme alaranjada com estratificação plano-paralela e estratificação cruzada de médio a grande porte. Destaca-se o registro de intervalos deformados lateralmente contínuos por centenas de quilômetros na zona de contato entre as formações Motuca e Sambaíba. Nestes, ocorrem pelitos com camadas contorcidas e brechadas (Formação Motuca) e arenitos com falhas/microfalhas sinsedimentares, laminação convoluta e diques de injeção preenchidos por argilitos (Formação Sambaíba), interpretados como sismitos induzidos por terremotos de alta magnitude (> 8 na escala Richter).