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.

Empilhamento sísmico pela composição de ondas planas

No presente trabalho de tese é apresentada uma nova técnica de empilhamento de dados sísmicos para a obtenção da seção de incidência normal ou afastamento fonte-receptor nulo, aplicável em meios bidimensionais com variações laterais de velocidade. Esta nova técnica denominada Empilhamento Sísmico pela Composição de Ondas Planas (empilhamento PWC) foi desenvolvida tomando como base os conceitos físicos e matemáticos da decomposição do campo de ondas em ondas planas. Este trabalho pode ser dividido em três partes: Uma primeira parte, onde se apresenta uma revisão da técnica de empilhamento sísmico convencional e do processo de decomposição do campo de ondas produzido a partir de fontes pontuais em suas correspondentes ondas planas. Na segunda parte, é apresentada a formulação matemática e o procedimento de aplicação do método de empilhamento sísmico pela composição de ondas planas. Na terceira parte se apresenta a aplicação desta nova técnica de empilhamento na serie de dados Marmousi e uma analise sobre a atenuação de ruído. A formulação matemática desta nova técnica de empilhamento sísmico foi desenvolvida com base na teoria do espalhamento aplicado a ondas sísmicas sob a restrição do modelo de aproximação de Born. Nesse sentido, inicialmente se apresenta a determinação da solução da equação de onda caustica para a configuração com afastamento fonte-receptor finito, que posteriormente é reduzido para a configuração de afastamento fonte-receptor nulo. Por outra parte, com base nessas soluções, a expressão matemática deste novo processo de empilhamento sísmico é resolvida dentro do contexto do modelo de aproximação de Born. Verificou-se que as soluções encontradas por ambos procedimentos, isto é, por meio da solução da equação da onda e pelo processo de empilhamento proposto, são iguais, mostrando-se assim que o processo de empilhamento pela composição de ondas planas produz uma seção com afastamento fonte-receptor nulo. Esta nova técnica de empilhamento basicamente consiste na aplicação de uma dupla decomposição do campo de ondas em onda planas por meio da aplicação de dois empilhamentos oblíquos (slant stack), isto é um ao longo do arranjo das fontes e outro ao longo do arranjo dos detectores; seguido pelo processo de composição das ondas planas por meio do empilhamento obliquo inverso. Portanto, com base nestas operações e com a ajuda de um exemplo de aplicação nos dados gerados a partir de um modelo simples, são descritos os fundamentos e o procedimento de aplicação (ou algoritmo) desta nova técnica de obtenção da seção de afastamento nulo. Como exemplo de aplicação do empilhamento PWC em dados correspondentes a um meio com variações laterais de velocidade, foi aplicado nos dados Marmousi gerados segundo a técnica de cobertura múltipla a partir de um modelo que representa uma situação geológica real. Por comparação da seção resultante com a similar produzida pelo método de empilhamento convencional, observa-se que a seção de afastamento nulo desta nova técnica apresenta melhor definição e continuidade dos reflectores, como também uma melhor caracterização da ocorrência de difrações. Por último, da atenuação de ruído aleatório realizada nos mesmos dados, observa-se que esta técnica de empilhamento também produz uma atenuação do ruído presente no sinal, a qual implica um aumento na relação sinal ruído.

Estudo das feições lineares observadas em dados aeromagnéticos da Bacia do Solimões a partir de imagens digitais de sombreamento do relevo anômalo total

A Bacia do Solimões é caracterizada por um padrão magnético complexo, expresso ao longo de toda sua extensão, e produzido pela superposição de anomalias magnéticas que se apresentam sob a forma de feições lineares. Tais feições refletem as diferentes atividades tectônicas que atuaram na Região Amazônica através do Pré-Cambriano e Fanerozóico. Neste trabalho, foi aplicado nos dados aeromagnéticos da Bacia do Solimões um método de processamento de imagens digitais de sombreamento do relevo magnético anômalo total que, dadas suas características metodológicas, permitindo utilizá-lo como um filtro direcional, possibilitou a definição de aspectos relevantes no contexto das relações entre os lineamentos magnéticos E/W, NE/SW, e NW/SE. Nesse sentido, foram identificados padrões de lineamentos magnéticos que refletem a existência de zonas de cisalhamento transcorrentes dextrais, orientadas preferencialmente na direção E/W. A interação entre os diversos segmentos transcorrentes promoveu o desenvolvimento de um regime predominantemente transpressivo, representado por falhas reversas associadas aos lineamentos magnéticos E/W e NE/SW; e duplexes direcionais formando falhas em flor positivas e negativas associadas, respectivamente, aos lineamentos magnéticos orientados na direção N70-80E e N70-80W. A análise quantitativa permitiu explicar dois aspectos importantes em relação às feições lineares observadas nas imagens digitais. O primeiro mostra, através de modelamentos baseados na superposição de prismas bidimensionais, que estas feições lineares magnéticas podem ser explicadas pela superposição de fontes profundas intraembasamento altamente magnéticas, e fontes rasas de alta frequência, sendo estas associadas a falhas reversas ao longo dos níveis de diabásio, presentes em forma de intrusões nos sedimentos paleozóicos da Bacia do Solimões. O segundo aspecto, baseado na utilização dos conceitos da cross-covariância, constata a presença de 'offsets' dextrais, associados aos lineamentos magnéticos NE—SW, ao longo da direção E—W. Este fato mostra, quantitativamente, que o padrão magnético desta região pode ser explicado pela presença de zonas de cisalhamento transcorrestes dextrais, cujos processos tectônicos associados foram fortemente condicionados por zonas de fraquezas pré-existentes (Pré-Cambrianas, Paleozóicas) durante o Mesozóico e Cenozóico.

Estudo do método de empilhamento SRC e aplicação em dados sintéticos e reais

O empilhamento sísmico é um método de imageamento utilizado para simulação de seções sísmicas com afastamento nulo (AN) a partir de dados de cobertura múltipla, produzindo um aumento da razão sinal/ruído na seção simulada. Este método é realizado usando o processamento convencional denominado ponto médio comum (PMC), através das chamadas correções sobretempo normal e sobretempo oblíquo (STN/STO). Alternativamente, um novo método de empilhamento sísmico, denominado superfície de reflexão comum (SRC) pode ser usado para simular seções sísmicas em afastamento nulo (AN). O método SRC tem como principais características: 1) a construção de uma aproximação do tempo de trânsito de reflexões primárias na vizinhança de um raio de reflexão normal; 2) esta aproximação de empilhamento tem como parâmetros os ângulos de emergência do raio normal e os raios de curvaturas de duas ondas hipotéticas denominadas onda ponto de incidência normal (onda PIN) e onda normal (onda N). O método pressupõe ainda que a velocidade da onda longitudinal na parte superior do modelo próxima à superfície da Terra seja conhecida. Para determinar adequadamente os parâmetros supracitados, fundamentais para o método de imageamento SRC, são necessárias estratégias de busca que utilizem casos particulares da aproximação de segunda ordem hiperbólica dos tempos de trânsito. As estratégias de busca apresentadas são: SRC pragmático-extendito e SRC global-local. Neste trabalho é apresentada uma revisão extensa dos fundamentos teóricos de método de empilhamento SRC, como também das estratégias de busca dos parâmetros SRC, dos algoritmos de implementação, e da delimitação de aberturas do operador de empilhamento SRC usando conceitos de zonas de Fresnel e o fator de estiramento da correção STN. Como um exemplo de aplicação dos atributos ou parâmetros SRC é apresentado o resultado da correção STN, sem estiramento do sinal sísmico, obtido utilizando as velocidades estimadas dos parâmetros SRC. Também, são apresentados os resultados da aplicação dos três métodos de empilhamento (PMC (STN/STO), SRC pragmático-extendido, SRC global-local) nos dados sintéticos Marmousi e real terrestre da bacia paleozóica Amazônica. Finalmente apresentam-se as discussões sobre os resultados, quanto à resolução da imagem sísmica e custo computacional.

Imageamento multifoco de refletores sísmicos

A simulação de uma seção sísmica de afastamento-nulo (ZO) a partir de dados de cobertura múltipla para um meio 2-D, através do empilhamento, é um método de imageamento de reflexão sísmica muito utilizado, que permite reduzir a quantidade de dados e melhorar a relação sinal/ruído. Segundo Berkovitch et al. (1999) o método Multifoco está baseado na Teoria do Imageamento Homeomórfico e consiste em empilhar dados de cobertura múltipla com distribuição fonte-receptor arbitrária de acordo com uma nova correção de sobretempo, chamada Multifoco. Esta correção de sobretempo esta baseada numa aproximação esférica local da frente de onda focalizante na vizinhança da superfície da terra. Este método permite construir uma seção sísmica no domínio do tempo de afastamento nulo aumentando a relação sinal/ruído. A técnica Multifoco não necessita do conhecimento a priori de um macro-modelo de velocidades. Três parâmetros são usados para descrever a aproximação de tempo de trânsito, Multifoco, os quais são: 1) o ângulo de emergência do raio de afastamento nulo ou raio de reflexão normal (β₀), 2) a curvatura da frente de onda no Ponto de Incidência Normal (R_NIP) e 3) curvatura da frente de Onda Normal (R_N). Sendo também necessário a velocidade próximo a superfície da terra. Neste trabalho de tese aplico esta técnica de empilhamento Multifoco para dados de cobertura múltipla referidos a modelos de velocidade constante e modelo heterogêneos, com o objetivo de simular seções sísmicas afastamento-nulo. Neste caso, como se trata da solução de um problema direto, o macro-modelo de velocidades é considerado conhecido a priori. No contexto do problema inverso tem-se que os parâmetros R_NIP, R_N e β₀ podem ser determinados a partir da análise de coerência aplicada aos dados sísmicos de múltipla cobertura. Na solução deste problema a função objetivo, a ser otimizada, é definida pelo cálculo da máxima coerência existente entre os dados na superfície de empilhamento sísmico. Neste trabalho de tese nos discutimos a sensibilidade da aproximação do tempo de trânsito usado no empilhamento Multifoco, como uma função dos parâmetros R_NIP, R_N e β₀. Esta análise de sensibilidade é feita de três diferentes modos: 1) a primeira derivada da função objetivo, 2) a medida de coerência, denominada semelhança, e 3) a sensibilidade no Empilhamento Multifoco.

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.

Estratégia do planejamento e otimização de sistemas sem fio, considerando redes interferentes: abordagem baseada em cross-layer

Apesar do aumento significativo do uso de redes locais sem fio (WLAN) nos últimos anos, aspectos de projeto e planejamento de capacidade da rede são ainda sistematicamente negligenciados durante a implementação da rede. Tipicamente um projeto de rede local sem fio é feito e instalado por profissionais de rede. Esses profissionais são extremamente experientes com redes cabeadas, mas são ainda geralmente pouco experientes com redes sem fio. Deste modo, as instalações de redes locais sem fio são desvantajosas pela falta de um modelo de avaliação de desempenho e para determinar a localização do ponto de acesso (PA), além disso, fatores importantes do ambiente não são considerados no projeto. Esses fatores se tornam mais importante quando muitos pontos de acesso (PAs) são instalados para cobrir um único edifício, algumas vezes sem planejamento de freqüência. Falhas como essa podem causar interferência entre células geradas pelo mesmo PA. Por essa razão, a rede não obterá os padrões de qualidade de serviço (QoS) exigidos por cada serviço. O presente trabalho apresenta uma proposta para planejamento de redes sem fio levando em consideração a influência da interferência com o auxílio de inteligência computacional tais como a utilização de redes Bayesianas. Uma extensiva campanha de medição foi feita para avaliar o desempenho de dois pontos de acesso (PAs) sobre um cenário multiusuário, com e sem interferência. Os dados dessa campanha de medição foram usados como entrada das redes Bayesianas e confirmaram a influência da interferência nos parâmetros de QoS. Uma implementação de algoritmo genético foi utilizado permitindo uma abordagem híbrida para planejamento de redes sem fio. Como efeito de comparação para otimizar os parâmetros de QoS, de modo a encontrar a melhor distância do PA ao receptor garantindo as recomendações do International Telecomunication Union (ITU-T), a técnica de otimização por enxame de partículas foi aplicada.

Métodos de otimização global aplicados na busca dos parâmetros SRC

O método de empilhamento por Superfície de Reflexão Comum (SRC) produz seções simuladas de afastamento nulo (AN) por meio do somatório de eventos sísmicos dos dados de cobertura múltipla contidos nas superfícies de empilhamento. Este método não depende do modelo de velocidade do meio, apenas requer o conhecimento a priori da velocidade próxima a superfície. A simulação de seções AN por este método de empilhamento utiliza uma aproximação hiperbólica de segunda ordem do tempo de trânsito de raios paraxiais para definir a superfície de empilhamento ou operador de empilhamento SRC. Para meios 2D este operador depende de três atributos cinemáticos de duas ondas hipotéticas (ondas PIN e N), observados no ponto de emergência do raio central com incidência normal, que são: o ângulo de emergência do raio central com fonte-receptor nulo (β₀) , o raio de curvatura da onda ponto de incidência normal (R_PIN) e o raio de curvatura da onda normal (R_N). Portanto, o problema de otimização no método SRC consiste na determinação, a partir dos dados sísmicos, dos três parâmetros (β₀, R_PIN, R_N) ótimos associados a cada ponto de amostragem da seção AN a ser simulada. A determinação simultânea destes parâmetros pode ser realizada por meio de processos de busca global (ou otimização global) multidimensional, utilizando como função objetivo algum critério de coerência. O problema de otimização no método SRC é muito importante para o bom desempenho no que diz respeito a qualidade dos resultados e principalmente ao custo computacional, comparado com os métodos tradicionalmente utilizados na indústria sísmica. Existem várias estratégias de busca para determinar estes parâmetros baseados em buscas sistemáticas e usando algoritmos de otimização, podendo estimar apenas um parâmetro de cada vez, ou dois ou os três parâmetros simultaneamente. Levando em conta a estratégia de busca por meio da aplicação de otimização global, estes três parâmetros podem ser estimados através de dois procedimentos: no primeiro caso os três parâmetros podem ser estimados simultaneamente e no segundo caso inicialmente podem ser determinados simultaneamente dois parâmetros (β₀, R_PIN) e posteriormente o terceiro parâmetro (R_N) usando os valores dos dois parâmetros já conhecidos. Neste trabalho apresenta-se a aplicação e comparação de quatro algoritmos de otimização global para encontrar os parâmetros SRC ótimos, estes são: Simulated Annealing (SA), Very Fast Simulated Annealing (VFSA), Differential Evolution (DE) e Controlled Rando Search - 2 (CRS2). Como resultados importantes são apresentados a aplicação de cada método de otimização e a comparação entre os métodos quanto a eficácia, eficiência e confiabilidade para determinar os melhores parâmetros SRC. Posteriormente, aplicando as estratégias de busca global para a determinação destes parâmetros, por meio do método de otimização VFSA que teve o melhor desempenho foi realizado o empilhamento SRC a partir dos dados Marmousi, isto é, foi realizado um empilhamento SRC usando dois parâmetros (β₀, R_PIN) estimados por busca global e outro empilhamento SRC usando os três parâmetros (β₀, R_PIN, R_N) também estimados por busca global.

Migração com amplitude verdadeira em meios com gradiente constante de velocidade

Este trabalho tem por objetivo a aplicação de um método de migração com amplitudes verdadeiras, considerando-se um meio acústico onde a velocidade de propagação varia linearmente com a profundidade. O método de migração é baseado na teoria dos raios e na integral de migração de Kirchhoff, procurando posicionar de forma correta os refletores e recuperar os respetivos coeficientes de reflexão. No processo de recuperação dos coeficientes de reflexão, busca-se corrigir o fator de espalhamento geométrico de reflexões sísmicas primárias, sem o conhecimento a priori dos refletores procurados. Ao considerar-se configurações fonte-receptor arbitrárias, as reflexões primárias podem ser imageadas no tempo ou profundidade, sendo as amplitudes do campo de ondas migrado uma medida dos coeficientes de reflexão (função do ângulo de incidência). Anteriormente têm sido propostos alguns algoritmos baseados na aproximação de Born ou Kirchhoff. Todos são dados em forma de um operador integral de empilhamento de difrações, que são aplicados à entrada dos dados sísmicos. O resultado é uma seção sísmica migrada, onde cada ponto de reflexão é imageado com uma amplitude proporcional ao coeficiente de reflexão no ponto. No presente caso, o processo de migração faz uso de um modelo com velocidade que apresenta uma distribuição que varia linearmente com a profundidade, conhecido também como gradiente constante de velocidade. O esquema de migração corresponde a uma versão modificada da migração de empilhamento por difração e faz uso explícito da teoria do raio, por exemplo, na descrição de tempos de trânsito e amplitudes das reflexões primárias, com as quais a operação de empilhamento e suas propriedades podem ser entendidas geometricamente. Efeitos como o espalhamento geométrico devido à trajetória do raio levam a distorção das amplitudes. Estes efeitos têm que ser corregidos durante o processamento dos dados sísmicos. Baseados na integral de migração de Kirchhoff e na teoria paraxial dos raios, foi derivada a função peso e o operador da integral por empilhamento de difrações para um modelo sísmico 2,5-D, e aplicado a uma serie de dados sintéticos em ambientes com ruído e livre de ruído. O resultado mostra a precisão e estabilidade do método de migração em um meio 2,5-D como ferramenta para obter informação sobre as propriedades de refletividade da subsuperfície da terra. Neste método não são levados em consideração a existência de caústicas nem a atenuação devido a fricção interna.

Migração em profundidade usando a solução numérica da equação da eiconal

Nos últimos anos tem-se verificado um interesse crescente no desenvolvimento de algoritmos de imageamento sísmico com a finalidade de obter uma imagem da subsuperfície da terra. A migração pelo método de Kirchhoff, por exemplo, é um método de imageamento muito eficiente empregado na busca da localização de refletores na subsuperficie, quando dispomos do cálculo dos tempos de trânsito necessários para a etapa de empilhamento, sendo estes obtidos neste trabalho através da solução da equação eiconal. Primeiramente, é apresentada a teoria da migração de Kirchhoff em profundidade baseada na teoria do raio, sendo em seguida introduzida a equação eiconal, através da qual são obtidos os tempos de trânsitos empregados no empilhamento das curvas de difrações. Em seguida é desenvolvido um algoritmo de migração em profundidade fazendo uso dos tempos de trânsito obtidos através da equação eiconal. Finalmente, aplicamos este algoritmo a dados sintéticos contendo ruído aditivo e múltiplas e obtemos como resultado uma seção sísmica na profundidade. Através dos experimentos feitos neste trabalho observou-se que o algoritmo de migração desenvolvido mostrou-se bastante eficiente e eficaz na reconstrução da imagem dos refletores.

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.

Paleoambiente e diagênese da formação Itaituba, carbonífero da bacia do Amazonas, com base em testemunho de sondagem, região de Uruará, Pará

A Formação Itaituba de idade carbonífera representa a sedimentação carbonática de depósitos transgressivos do Grupo Tapajós da Bacia do Amazonas. A sucessão Itaituba é interpretada como depósitos de planície de maré mista, constituídos de calcários fossilíferos, dolomitos finos, arenitos finos a grossos e subordinadamente siltitos avermelhados, evaporitos e folhelhos negros. A análise de fácies e microfácies do testemunho de sondagem da região de Uruará, Estado do Pará, permitiu individualizar dezenove fácies agrupadas em cinco associações: planície de maré (AF1), canal de maré (AF2), laguna (AF3), barra bioclástica (AF4) e plataforma externa (AF5). AF1 é composta por arenito fino com rip-up clasts e gretas de contração, marga com grãos de quartzo e feldspato, dolomudstone laminado com grãos terrígenos e dolomito fino silicificado, com intercalação de argilito com grãos de quartzo disseminados, dolomitizado e localmente com sílica microcristalina. AF2 consiste em arenito médio a grosso com estratificação cruzada acanalada, recoberta por filmes pelíticos nos foresets, arenito muito fino a fino com acamamento wavy, siltito laminado com falhas sinsedimentares e acamamento convoluto. AF3 é constituída de siltito vermelho maciço, mudstone com fósseis, floatstone com braquiópodes e pirita disseminada e mudstone maciço com frequentes grãos de quartzo. AF4 e AF5 exibem abundantes bioclastos representados por espinhos e fragmentos de equinodermas, conchas, fragmentos e espinhos de braquiópodes, ostracodes, foraminíferos, algas vermelhas e conchas de bivalves. AF4 é formada por grainstone oolítico fossilífero e grainstone com terrígenos principalmente grãos de quartzo monocristalino e AF5 se compõe de wackestone fossilífero, wackestone com terrígenos e mudstone maciço com grãos de quartzo monocristalino. Subarcósios (AF1), arcósios (AF2) e arcósios líticos (AF2) são os tipos de arenitos da sucessão Itaituba e apresentam como principais constituintes grãos de quartzo monocristalino e policristalino, K-feldspato, plagioclásio, pirita, muscovita detrítica, fragmento de rocha pelítica, metamórfica e chert e raros bioclastos. O cimento é de calcita espática não ferrosa, óxido/hidróxido de ferro e sobrecrescimento de sílica. A porosidade é intergranular, móldica e às vezes alongada, sem permeabilidade perfazendo até 11% da rocha. Os processos diagenéticos dos arenitos são compactação física, sobrecrescimento de sílica, cimentação de calcita, formação de matriz diagenética, compactação química, substituição de grãos, autigênese de pirita, formação de óxido/hidróxido de ferro e alteração do plagioclásio. Os processos diagenéticos dos carbonatos são: micritização, neomorfismo, colomitização, fraturamento, compactação química, cimentação de calcita, dissolução secundária e autigênese de minerais. A sucessão da Formação Itaituba representa um sistema de laguna/planície de maré ligada a uma plataforma marinha carbonática. Planícies de maré desenvolveram-se nas margens das lagunas e eram periodicamente supridas por influxos de terrígenos finos (silte) que inibiam a precipitação carbonática. Barras bioclásticas eram cortadas por canais de maré (inlet) que conectavam a laguna com a plataforma rasa rica em organismos bentônicos.

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

Proveniência e idade do metamorfismo das rochas da Faixa Brasília, na Região de Tapira (SW de Minas Gerais)

A Faixa Brasília na região de Tapira, no SW do estado de Minas Gerais, exibe quatro domínios litotectônicos diferentes, imbricados através de falhas de empurrão, com vergência para o Cráton do São Francisco. Foram conduzidos estudos isotópicos através do método Sm/Nd, que revelam diferentes idades modelo e de metamorfismo para estas escamas. As rochas da escama inferior apresentam idade de metamorfismo de 543 Ma, enquanto a escama que a superpõem, apresentam idade de metamorfismo de 581 Ma. As idades modelo (TDM) são similares para as duas escamas, variando entre 1,7 e 2,2 Ga. Em função das características litológicas essas escamas são interpretadas como derivadas de rochas depositadas em plataforma continental distal, que tiveram como fonte principal rochas de idades paleoproterozóicas do Cráton do São Francisco. As rochas da escama superior apresentam idade de metamorfismo de 612 Ma e uma distribuição bimodal das idades modelo (TDM), 1,3 e 1,9 - 2,0 Ga. Em função das características litológicas e isotópicas as rochas dessas escamas são interpretadas como depositadas em ambiente de talude continental ou assoalho oceânico, tendo os sedimentos fonte mista, de idades paleo- e mesoproterozóicas do Cráton do São Francisco. Embora, com as limitações impostas tanto pelo pequeno número de amostras analisadas, quanto pelas incertezas do método Sm/Nd, interpreta-se que o metamorfismo não foi síncrono nas diferentes escamas. Isto é esperado em um sistema de cavalgamento, no qual as escamas mais metamórficas justapõem-se às escamas menos metamórficas.