Avaliação dos reservatórios do Campo de Ubarana-RN, com baixas resistividades

Desenvolveu-se neste trabalho um novo modelo de avaliação quantitativa para arenitos argilosos, baseado na bibliografia e no estudo da disposição dos argilominerais dentro dos poros das rochas reservatório. Este novo modelo leva em consideração a contribuição de duas condutividades extras, além da condutividade eletrolítica das rochas: uma se deve à contribuição efetiva da condutância superficial dos argilominerais, aumentada ou diminuida em função da capacidade de troca catiônica dos mesmos, e outra independente e em paralelo, devido a condutância que se desenvolve em uma rede microporosa contínua formada por argilominerais sobre o arcabouço das rochas. Para o estudo da validade da equação proposta, foram feitos vários gráficos comparativos entre as equações e/ou modelos existentes n” bibliografia especializada, deles resultando ótimas correlações, principalmente com a tradicional equação de Waxman & Smits. Foi feita uma comparação do novo modelo em três seções distintas do Campo de Ubarana, Rio Grande do Norte, duas com baixas resistividades supostamente com hidrocarbonetos e uma outra efetivamente produtora. Verificou-se que as baixas resistividades resultam do fato das seções estudadas serem subsaturadas, com menos de 50% de óleo e aproximadamente igual ou maior proporção de água relativamente salgada (85.000 ppm de NaCl equivalente). Esta mesma água, sem dúvida, muito contribui para as baixas resistividades, por formar a fase eletricamente condutiva das referidas seções. Sendo a equação de Waxman & Smits mundialmente reconhecida, pode-se afirmar que a equação proposta neste trabalho tem consistência teórica e prática e para o caso particular do campo de Ubarana, mostrou-se mais coerente com o histórico de produção dos poços estudados, do que as demais equações existentes e testadas.

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.

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 do embasamento e da variação da densidade com a profundidade em bacias sedimentares

O presente método postula uma variação hiperbólica para o contraste de densidade de uma bacia sedimentar em função da profundidade, e tem dois objetivos: (1) delinear o relevo do embasamento de uma bacia, conhecendo-se a anomalia gravimétrica, o contraste de densidade na superfície da bacia e o fator de decaimento do contraste de densidade com a profundidade; (2) estimar, além do relevo, o valor do contraste de densidade na superfície de uma bacia sedimentar e o fator de decaimento do contraste de densidade com a profundidade, sendo fornecida a anomalia gravimétrica e a profundidade do embasamento em alguns pontos da bacia. Nos dois casos o modelo interpretativo é um conjunto de prismas retangulares verticais justapostos, cujas espessuras, que são parâmetros a serem estimados, representam a profundidade da interface de separação entre os sedimentos e o embasamento. As soluções obtidas nos dois problemas acima formulados são estáveis devido à incorporação de informações adicionais sobre a suavidade do relevo estimado, e o conhecimento da profundidade do relevo do embasamento em alguns pontos, fornecido por furos de sondagem. O método foi testado em anomalias gravimétricas sintéticas produzidas pela simulação de bacias sedimentares com relevos suaves. Os resultados mostraram relevos com boa resolução e valores estimados do contraste de densidade na superfície da bacia e do fator de decaimento do contraste de densidade com a profundidade, próximos aos verdadeiros, indicando dessa maneira o potencial do método em interpretações gravimétricas de bacias sedimentares. O método foi aplicado à anomalia Bouguer da Bacia do Recôncavo, Brasil, delineando um relevo com um valor para a profundidade máxima de cerca de 6 km, semelhante ao estimado em interpretações sísmicas. As estimativas para o contraste de densidade na superfície da Bacia e o fator de decaimento com a profundidade foram -0,30 g/cm³ e 30 km, respectivamente, produzindo uma estimativa para a compactação máxima dos sedimentos da ordem de 4%.

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.

Inversão gravimétrica do relevo de bacias extencionais através da variação total

Apresentamos um método de inversão de dados gravimétricos para a reconstrução do relevo descontínuo do embasamento de bacias sedimentares, nas quais o contraste de densidade entre o pacote sedimentar e o embasamento são conhecidos a priori podendo apresentar-se constante, ou decrescer monotonicamente com a profundidade. A solução é estabilizada usando o funcional variação total (VT), o qual não penaliza variações abruptas nas soluções. Comparamos o métodoproposto com os métodos da suavidade global (SG), suavidade ponderada (SP) e regularização entrópica (RE) usando dados sintéticos produzidos por bacias 2D e 3D apresentando relevos descontínuos do embasamento. As soluções obtidas com o método proposto foram melhores do que aquelas obtidas com a SG e similares às produzidas pela SP e RE. Por outro lado, diferentemente da SP, o método proposto não necessita do conhecimento a priori sobre a profundidade máxima do embasamento. Comparado com a RE, o método VT é operacionalmente mais simples e requer a especificação de apenas um parâmetro de regularização. Os métodos VT, SG e SP foram aplicados, também, às seguintes áreas: Ponte do Poema (UFPA), Steptoe Valley (Nevada, Estados Unidos), Graben de San Jacinto (Califórnia, Estados Unidos) e Büyük Menderes (Turquia). A maioria destas áreas são caracterizadas pela presença de falhas com alto ângulo. Em todos os casos, a VT produziu estimativas para a topografia do embasamento apresentando descontinuidades bruscas e com alto ângulo, em concordância com a configuração tectônica das áreas em questão.

O Diagnóstico físico-conservacionista-DFC como subsídio à gestão ambiental da bacia hidrográfica do rio Quebra-Perna, Ponta Grossa-PR

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

Modelagem acústica por diferenças finitas e elementos finitos em 2-D e 2,5-D

A modelagem acústica fornece dados úteis para avaliação de metodologias de processamento e imageamento sísmico, em modelos com estrutura geológica complexa. Esquemas de diferenças finitas (DF) e elementos finitos (EF) foram implementados e avaliados em modelos homogêneos e heterogêneos. O algoritmo de diferenças finitas foi estendido para o caso 2,5-D em modelos com densidade variável. Foi apresentada a modelagem de alvos geológicos de interesse exploratório existentes na Bacia Paleozóica do Solimões na Amazônia. Reflexões múltiplas de longo período produzidas entre a superfície livre e a discordância Cretáceo-Paleozóica, a baixa resolução da onda sísmica nas proximidades do reservatório e as fracas reflexões na interface entre as rochas reservatório e as rochas selantes são as principais características dos dados sintéticos obtidos, os quais representam um grande desafio ao imageamento sísmico.

Morfometria de compartimentos hidrológicos e qualidade hídrica do município de Batatais-SP

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

Caracterização dos resíduos cauliníticos das indústrias de mineração de caulim da Amazônia como matéria-prima para produção de pozolanas de alta reatividade

As regiões do Rio Capim e do Rio Jari são os principais distritos caulinitícos da Região Amazônica, detentores das maiores reservas brasileiras de caulim de alta qualidade para aplicações como cobertura de papel. O caulim é lavrado e beneficiado por três grandes companhias que são responsáveis pela geração de aproximadamente 500 mil de toneladas anuais de um resíduo decorrente da etapa de centrifugação. Esse resíduo, na forma de polpa, é depositado em lagoas de sedimentação que ocupam grandes extensões de áreas. O objetivo da pesquisa foi investigar as características físicas, químicas e mineralógicas dos resíduos processados da Região do Jari e do Capim, de modo a avaliar se atendem aos requisitos como matéria-prima para a produção de uma pozolana de alta reatividade, o metacaulim, adição mineral que incorporada ao cimento Portland proporciona alto desempenho às misturas de concreto e argamassas. Os resíduos foram caracterizados por difração de raios X, análise térmica, espectroscopia de infravermelho, microscopia eletrônica de varredura, fluorescência de raios X e difração a laser. Ambos os resíduos são constituídos por no mínimo 92% de caulinita de baixa granulometria, cujas áreas superficiais específicas são superiores a 8 m²/g e os diâmetros médios de partículas inferiores a 1 µm. Os teores de sílica livre (quartzo) não foram superiores a 3%. O alto grau de concentração de caulinita destes resíduos dispensa os rígidos parâmetros de controle de remoção de impurezas, normalmente empregados na produção deste tipo de pozolana. O caulim do Rio Jari, com quantidade de defeitos na estrutura cristalina superior ao da caulinita do Rio Capim, proporcinou maior grau de desidroxilação a uma temperatura mais baixa quando calcinado, indicando a possibilidade de redução de gastos com energia para a produção da adição mineral. Os resultados das análises foram convergentes e apontam ambos os caulins estudados como excelentes matérias-primas para a produção do metacaulim de alta reatividade.

Caracterização mineralógica e geoquímica e estudo das transformações de fase do caulim duro da região do Rio Capim, Pará

A região do Rio Capim (Nordeste do Estado do Pará), destaca-se nacionalmente por suas grandes reservas de caulim para cobertura de papel. O minério extraído está localizado, em média, a 20 m de profundidade, recoberto por sedimentos argilo-arenosos da Formação Barreiras, e de um nível de caulim duro, também conhecido como flint ou semi-flint, considerado como estéril em função do teor de ferro elevado que inviabiliza sua aplicação para cobertura. Este trabalho objetiva a caracterização mineralógica e geoquímica do caulim duro além de acompanhar as transformações mineralógicas sofridas em função da temperatura de queima, com vistas a sua possível utilização como matéria-prima cerâmica. Foram utilizados dois tipos de caulim duro como material de partida, em função dos teores de ferro: o Caulim Duro Branco (CDB) e o Caulim Duro Ferruginoso (CDF), este último com teor de Fe2O3 de 10,36%. Os resultados indicam composição mineralógica dominada por caulinita, além de anatásio como acessório. No CDF ocorrem ainda goethita e hematita. Os estudo das transformações térmicas indicam que o início de formação de mullita se dá em temperaturas diferentes para as amostras estudadas.

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.