Modelo fractal para resistividade complexa de rochas: interpretação petrofísica e aplicação à exploração geoelétrica

Rochas contendo metálicos disseminados ou partículas de argila em ambiente natural onde soluções eletrolíticas normalmente preenchem os poros das rochas, exibem um tipo de polarização em baixas freqüências conhecido como polarização induzida. Nesta tese foi desenvolvido um novo modelo para descrever o fenômeno de polarização das rochas, não apenas em baixas freqüências, mas compreendendo todo o espectro eletromagnético, possível de utilização na prospecção geoelétrica. Este novo modelo engloba a maioria dos modelos utilizados até o momento como casos especiais, além de superar as limitações dos mesmos. Seu circuito analógico inclui uma impedância não linear do tipo r (iwtf)^-n que simula o efeito das superfícies rugosas das interfaces entre os grãos bloqueadores (partículas metálicas e/ou de argilas) e o eletrólito. A impedância de Warburg generalizada está em série com a resistência dos grãos bloqueadores da passagem de corrente e em paralelo com a impedância da dupla camada associada a essas interfaces. Esta combinação está em série com a resistência do eletrólito nas passagens dos poros bloqueados. Os canais não bloqueados são representados por uma resistência que corresponde à resistividade normal CC da rocha. A combinação desta resistência com a capacitância "global" da rocha é finalmente conectada em paralelo ao resto do circuito mencionado acima. Os parâmetros deste modelo incluem a resistividade CC (p₀), a cargueabilidade (m), três tempos de relaxação (t, Tf and T2), um fator de resistividade de grãos (δ_r), e o expoente de freqüência (η). O tempo de relaxação fractal (Tf), e o expoente de frequencia (η) estão relacionados à geometria fractal das interfaces rugosas entre os minerais condutivos (grãos metálicos e/ou partículas de argila bloqueando os canais dos poros) e o eletrólito. O tempo de relaxação (T) é um resultado da relaxação em baixa freqüência das duplas camadas elétricas formadas nas interfaces eletrólito-cristais, enquanto (T₀) é o tempo de relaxação macroscópico da amostra como um todo. O fator de resistividade dos grãos (δ_r) relaciona a resistividade dos grãos condutivos com o valor de resistividade CC da rocha. A resistividade CC da rocha (p₀), e δ_r estão relacionados à porosidade, à condutividade do eletrólito e às relações mineralógicas entre a matriz e os grãos condutivos. O modelo foi testado sobre um intervalo largo de freqüências contra dados experimentais de amplitude e fase da resistividade bem como para dados de constante dielétrica complexa. Os dados utilizados neste trabalho foram obtidos a partir da digitalização de dados experimentais publicados, obtidos por diversos autores e englobando amostras de rochas sedimentares, ígneas e metam6rficas. É mostrado neste trabalho que os parâmetros deste modelo permitem identificar diferenças texturais e mineralógicas nas rochas. Bote modelo foi introduzido, primeiramente, como propriedade intrínseca de um semiespaço homogêneo sendo demonstrado, neste trabalho, que a resposta observada em superfície reflete as propriedades intrínsecas do meio polarizável, sendo o acoplamento eletromagnético desprezível em freqüências menores que 104 Hz. Em seguida, o meio polarizável foi embebido em um pacote de N camadas sendo demonstrado que os parâmetros fractais do meio polarizável podem ser obtidos do levantamento em superfície para diferentes espessuras dessa camada. Isto justifica a utilização pura e simples de modelos de polarização desenvolvidos para amostras em laboratório para ajustar dados de campo, o que vem sendo feito sem uma justificativa bem fundamentada. Estes resultados demonstram a importância para a prospecção geolétrica do modelo proposto nesta tese.

Uma nova abordagem para interpretação de anomalias gravimétricas regionais e residuais aplicada ao estudo da organização crustal: exemplo da Região Norte do Piauí e Noroeste do Ceará

A interpretação de anomalias gravimétricas é de grande importância no estudo de feições geológicas que ocorrem na crosta terrestre. Esta interpretação é, no entanto, dificultada pelo fato das anomalias gravimétricas serem resultantes da soma total dos efeitos produzidos por todos os contrastes de densidades de subsuperfície. Desse modo, com o intuito de separar efeitos de feições mais profundas de efeitos de feições mais rasas, bem como a caracterização da geometria desses dois conjuntos de feições, apresentamos um método de separação das componentes regional e residual do campo e a subsequente interpretação de cada componente. A separação regional-residual de dados gravimétricos é efetuada através da aproximação do campo regional por um polinômio ajustado ao campo observado por um método robusto. Este método é iterativo e usa como aproximação inicial a solução obtida através do ajuste polinomial pelo método dos mínimos quadrados. O método empregado minimiza a influência de observações contendo forte contribuição do campo residual no ajuste do campo regional. A componente regional obtida a partir da separação regional-residual é transformada em um mapa de distâncias verticais em relação a um nível de referência. Esta transformação compreende duas etapas. A primeira consiste na obtenção da continuação para baixo da componente regional, que é pressuposta ser causada por uma interface suave separando dois meios homogêneos, representando a interface crosta-manto, cujo contraste de densidade é supostamente conhecido. A segunda consiste na transformação do mapa de continuação para baixo em um mapa de distâncias verticais entre o nível de continuação (tomado como nível de referência) e a interface. Este método apresenta duas dificuldades. A primeira está ligada à instabilidade, havendo portanto a necessidade do emprego de um estabilizador o que acarreta a perda de resolução das feições que se desejam mapear. A segunda, inerente ao método gravimétrico, consiste na impossibilidade da determinação das profundidades absolutas da interface em cada ponto, bastando entretanto o conhecimento da profundidade absoluta em um ponto, através de informação independente, para que todas as outras profundidades absolutas sejam conhecidas. A componente residual obtida a partir da separação regional-residual é transformada em um mapa de contrastes de densidade aparente. Esta transformação consiste no cálculo do contraste de densidade de várias fontes prismáticas através de uma inversão linear pressupondo que as fontes reais estejam das a uma placa horizontal, com contrastes de densidade variando apenas nas direções horizontais. O desempenho do método de separação regional-residual apresentado foi avaliado, através de testes empregando dados sintéticos, fornecendo resultados superiores em relação aos métodos dos mínimos quadrados e da análise espectral. O método de interpretação da componente regional teve seu desempenho avaliado em testes com dados sintéticos onde foram produzidos mapeamentos de interfaces bem próximas das estruturas reais. O limite de resolução das feições que se desejam mapear depende não só do grau do polinômio ajustante, como também da própria limitação inerente ao método gravimétrico. Na interpretação da componente residual é necessário que se postule ou tenha informação a priori sobre a profundidade do topo e espessura da placa onde as fontes estão supostamente confinadas. No entanto, a aplicação do método em dados sintéticos, produziu estimativas razoáveis para os limites laterais das fontes, mesmo na presença de fontes interferentes, e pressupondo-se valores para profundidade do topo e espessura da placa, diferentes dos valores verdadeiros. A ambiguidade envolvendo profundidade do topo, espessura e densidade pode ser visualizada através de gráficos de valores de densidade aparente contra profundidade do topo presumida para a placa para vários valores postulados para a espessura da placa. Estes mesmos gráficos permitem, pelo aspecto das curvas, a elaboração de uma interpretação semi-quantitativa das profundidades das fontes reais. A seqüência dos três métodos desenvolvidos neste trabalho foi aplicada a dados gravimétricos da região norte do Piauí e noroeste do Ceará levando a um modelo de organização crustal que compreende espessamentos e adelgaçamentos crustais associados a um evento compressivo que possibilitou a colocação de rochas densas da base da crosta a profundidades rasas. Este modelo ê compatível com os dados geológicos de superfície. É ainda sugerida a continuidade, por mais 200 km em direção a sudoeste, do Cinturão de Cisalhamento Noroeste do Ceará por sob os sedimentos da Bacia do Parnaíba, com base nas evidências fornecidas pela interpretação da anomalia residual. Embora esta seqüência de métodos tenha sido desenvolvida com vistas ao estudo de feições crustais de porte continental, ela também pode ser aplicada ao estudo de feições mais localizadas como por exemplo no mapeamento do relevo do embasamento de/bacias sedimentares onde os sedimentos são cortados por rochas intrusivas mais densas.

Análise de resolução em tomografia interpoços em 3-D

A aplicação de métodos de inversão aos dados de múltiplos levantamentos sísmicos interpoços permite a reconstrução de modelos de vagarosidade em 3-D de alta resolução adequados para monitoramento de processos de recuperação avançada de petróleo e caracterização de reservatórios. Entretanto, a falta de cobertura volumétrica uniforme de raios de levantamentos interpoços exige informação adicional ao sistema tomográfico para obtenção de soluções estáveis. A discretização do modelo em uma malha 3-D com células prismáticas triangulares e a decomposição em valores singulares são utilizadas para avaliar a reconstrução tomográfica em 3-D. O ângulo da projeção de modelos-alvo no subespaço ortogonal ao espaço nulo efetivo da matriz tomográfica é um critério adequado para se otimizar a malha de discretização do modelo interpretativo e a geometria de aquisição dos dados de modo a melhorar o condicionamento da reconstrução tomográfica. Esta abordagem pode ser utilizada durante as iterações lineares para redefinir a malha ou avaliar a necessidade de informação a priori adicional ao sistema tomográfico.

Caracterização de reservatórios fraturados através de dados de ondas qP em levantamentos VSP Walkaway

Esta dissertação apresenta um método para estimar a orientação de fratura de um meio partindo do pressuposto que este se comporta efetivamente como um meio transversalmente isotrópico (TI). Desta forma, a orientação de fratura é obtida a partir da estimativa do eixo de simetria de um meio TI. Esta estimativa é obtida a partir dos vetores de vagarosidade e polarização de ondas qP medidos em experimentos VSP Walkaway. O processo de inversão é baseado em equações linearizadas dos vetores de vagarosidade e polarização de ondas qP e os parâmetros fracamente anisotrópicos do meio. São apresentados testes numéricos, em que a sensibilidade a fatores como grau de anisotropia, geometria do levantamento, tipo de onda utilizado e nível de ruído são analisados. São também apresentados os resultados para um conjunto de dados reais.

Delineamento do pé do talude na margem continental do Ceará através da integração de dados geológicos e geofísicos

A metodologia de interpretação integrada dos dados geológicos e geofísicos observados em um perfil da margem continental do Ceará possibilita a identificação e a integração de características peculiares a cada tipo de dado. Dessa forma, é possível se definir a localização mais provável de feições estruturais importantes, tais como a fronteira entre as crostas continental e oceânica e o pé do talude, objeto do presente estudo. Segundo o Artigo 76 (parágrafo 4, item b) da Convenção das Nações Unidas para o Direito do Mar, o pé do talude é definido como o ponto de variação máxima do gradiente do talude na sua base. Entretanto, essa definição, apesar de simples no contexto fisiográfico, não é suficiente para propiciar a localização do pé do talude como preconiza a Convenção, razão pela qual se aplicam os métodos geofísicos. Dentro do contexto geofísico-geológico está implícita a interpretação quantitativa das anomalias gravimétricas ar-livre, que possibilita o delineamento do modelo geofísico representando a subsuperfície, cuja finalidade é subsidiar geologicamente a interpretação integrada dos dados aludidos. Um procedimento automático de ajuste de curvas combinando as técnicas de inversão de busca sistemática e a que utiliza derivadas foi usado com o propósito de gerar o modelo geofísico. A aplicação rigorosa de vínculos preliminarmente e a constante reavaliação desses vínculos através de um processo interativo entre a sísmica e a gravimetria, gerado durante a interpretação quantitativa das anomalias ar-livre, possibilitaram que o modelo geofísico final estivesse dentro dos padrões geológicos para área, notadamente quanto ao equilíbrio isostático (Teoria de Airy). O objetivo do presente trabalho é se estudar as características geológicas e geofísicas observadas ao longo de um perfil da margem continental do Ceará (LEPLAC III), notadamente quanto ao pé do talude, buscando estabelecer a aplicabilidade de uma metodologia de interpretação integrada desses dados, cuja finalidade é se definir de forma sistemática a localização mais provável para esta feição fisiográfica. A metodologia de interpretação integrada dos dados geológicos e geofísicos empregada mostrou-se eficiente para este objetivo. Foi possível se integrar: (i) a localização fisiográfica (distância da costa e profundidade) do pé do talude; (ii) a zona de instabilidade tectônica evidenciada pelos falhamentos, comuns nesta região; (iii) ao fim de uma zona magnética perturbada, associado a um ponto de mínimo na curva de anomalia magnética, e que possivelmente delimita o início de uma zona magnética quieta, denominada de anomalia E e (iv) um ponto de inflexão na curva de anomalia ar-livre, associado ao efeito gravimétrico do contraste de densidades entre as crostas continental, os sedimentos e a água do mar, evidenciado pela geometria do talude. Foi possível ainda se definir a localização mais provável para a fronteira entre as crostas continental e oceânica. Dada a rigorosidade na aplicação das técnicas de inversão e dos vínculos é provável que as correlações das características intrínsecas a cada tipo de dado efetuadas na conclusão desse trabalho tenham fundamento e possam ser confirmadas. A condição para isto é a aplicação da metodologia aqui estabelecida em um número maior de perfis.

Estudo do transporte eletrônico em nanoestruturas baseadas em carotenoides e tétrades com fulereno C60

Neste Trabalho é apresentado um estudo teórico da estrutura eletrônica de uma molécula de fulereno C60 com junções em quatro terminais baseados em grupos doadores de elétron – etratiofulvaleno (TTF) – e grupos aceitadores de elétrons – fenilpropanodinila (FPP) e dispositivos moleculares baseados em derivados dos Carotenoides. O mecanismo de transporte investigado para os derivados dos Carotenoides foram utilizados para o melhor entendimento das curvas de Fowler- Nordheim (FN) e Millikan-Lauritsen (ML) para os sistemas baseados em fulereno C60. Em todos os casos foi possível confirmar que a análise empírica de Millikan-Lauritsen (ML) também é suficiente para descrever em todos os aspectos a espectroscopia de voltagem de transição (TVS). Para estudar os sistemas, foram feitas otimizações de geometria sistematicamente e observado uma transferência eletrônica calculada por métodos derivados de Hartree-Fock e Teoria do Funcional Densidade (DFT). Os resultados apresentados mostram um estudo detalhado do rearranjo de carga molecular para a estrutura, que sob a ação de um campo elétrico externo apontou que o transporte de carga está diretamente ligado ao tipo de junção que esse sistema é submetido de forma que a voltagem aplicada é intensa o bastante para criar um potencial de saturação nos sistemas em estudo: fulereno C60 com três terminais de tetratiofulvaleno e um terminal de fenil-propanodinila (C60-(TTF)₃-FPP); fulereno C60 com quatros terminais de fenil-propanodinila (C60-(FPP)₄). Os resultados mostram que se tem um retificador molecular que pode trabalhar corretamente como um retificador macroscópico.

Estudo comparativo das técnicas de elementos finitos e equação integral na modelagem eletromagnética bidimensional

Os métodos numéricos de Elementos Finitos e Equação Integral são comumente utilizados para investigações eletromagnéticas na Geofísica, e, para essas modelagens é importante saber qual algoritmo é mais rápido num certo modelo geofísico. Neste trabalho são feitas comparações nos resultados de tempo computacional desses dois métodos em modelos bidimensionais com heterogeneidades condutivas num semiespaço resistivo energizados por uma linha infinita de corrente (com 1000Hz de freqüência) e situada na superfície paralelamente ao "strike" das heterogeneidades. Após a validação e otimização dos programas analisamos o comportamento dos tempos de processamento nos modelos de corpos retangulares variandose o tamanho, o número e a inclinação dos corpos. Além disso, investigamos nesses métodos as etapas que demandam maior custo computacional. Em nossos modelos, o método de Elementos Finitos foi mais vantajoso que o de Equação Integral, com exceção na situação de corpos com baixa condutividade ou com geometria inclinada.

Interfaces gráficas de entrada e saída da eGs (electrical GEOPHYSICS suite)

Apresentamos uma interface gráfica para o programa eGs (electrical GEOPHYSICS suíte) usado na modelagem e interpretação de dados elétricos e eletromagnéticos na exploração geofísica. Esta dissertação corresponde a parcela do programa que trata das interfaces de entrada de dados de modelos geoelétricos bi e tridimensionais e das interfaces de apresentação dos resultados da modelagem como também dos dados observados. O programa contém interfaces para criação, modificação e visualização de modelos geoelétricos bi e tridimensionais, e para apresentação da resposta do modelo através de gráficos e imagens. Este trabalho resultou num programa com diversas opções para alterar a geometria e as propriedades físicas do modelo geoelétricos e recursos de visualização para facilitar a interpretação de dados geofísicos. A interface gráfica pode ser usada para criar modelos geoelétricos dos métodos elétrico ou eletromagnético mais usados e analisar as respostas dos modelos e de dados observados em campanha de campo. O programa é suficientemente flexível de modo que novas interfaces gráficas de outros métodos elétricos e eletromagnéticos podem ser facilmente implementadas.

Inversão de dados de sísmica de refração profunda a partir da curva tempo-distância

O trabalho em pauta tem como objetivo o modelamento da crosta, através da inversão de dados de refração sísmica profunda, segundo camadas planas horizontais lateralmente homogêneas, sobre um semi-espaço. O modelo direto é dado pela expressão analítica da curva tempo-distância como uma função que depende da distância fonte-estação e do vetor de parâmetros velocidades e espessuras de cada camada, calculado segundo as trajetórias do raio sísmico, regidas pela Lei de Snell. O cálculo dos tempos de chegada por este procedimento, exige a utilização de um modelo cujas velocidades sejam crescentes com a profundidade, de modo que a ocorrência das camadas de baixa velocidade (CBV) é contornada pela reparametrização do modelo, levando-se em conta o fato de que o topo da CBV funciona apenas como um refletor do raio sísmico, e não como refrator. A metodologia de inversão utilizada tem em vista não só a determinação das soluções possíveis, mas também a realização de uma análise sobre as causas responsáveis pela ambiguidade do problema. A região de pesquisa das prováveis soluções é vinculada segundo limites superiores e inferiores para cada parâmetro procurado, e pelo estabelecimento de limites superiores para os valores de distâncias críticas, calculadas a partir do vetor de parâmetros. O processo de inversão é feito utilizando-se uma técnica de otimização do ajuste de curvas através da busca direta no espaço dos parâmetros, denominado COMPLEX. Esta técnica apresenta a vantagem de poder ser utilizada com qualquer função objeto, e ser bastante prática na obtenção de múltiplas soluções do problema. Devido a curva tempo-distância corresponder ao caso de uma multi-função, o algoritmo foi adaptado de modo a minimizar simultaneamente várias funções objetos, com vínculos nos parâmetros. A inversão é feita de modo a se obter um conjunto de soluções representativas do universo existente. Por sua vez, a análise da ambiguidade é realizada pela análise fatorial modo-Q, através da qual é possível se caracterizar as propriedades comuns existentes no elenco das soluções analisadas. Os testes com dados sintéticos e reais foram feitos tendo como aproximação inicial ao processo de inversão, os valores de velocidades e espessuras calculados diretamente da interpretação visual do sismograma. Para a realização dos primeiros, utilizou-se sismogramas calculados pelo método da refletividade, segundo diferentes modelos. Por sua vez, os testes com dados reais foram realizados utilizando-se dados extraídos de um dos sismogramas coletados pelo projeto Lithospheric Seismic Profile in Britain (LISPB), na região norte da Grã-Bretanha. Em todos os testes foi verificado que a geometria do modelo possui um maior peso na ambiguidade do problema, enquanto os parâmetros físicos apresentam apenas suaves variações, no conjunto das soluções obtidas.

Modelamento e correção de descentralização das imagens de tempo de trânsito

As imagens da parede do poço obtidas através da ferramenta UBI (Ultrasonic Borehole Imager) são amplamente utilizadas por geólogos e engenheiros de petróleo para identificar eventos geológicos nas paredes de poços abertos e na inspeção dos tubos de revestimento, uma vez que este perfil praticamente fornece uma fotografia da parede do poço. As ferramentas de imageamento acústico produzem imagens do tempo de trânsito e da amplitude do pulso acústico gerado pela ferramenta e refletido na parede do poço. Entretanto, estas imagens podem ter uma interpretação não realista, uma vez que elas podem ter seu aspecto alterado em razão do movimento da ferramenta no interior do poço. Este trabalho apresenta o modelamento das imagens de tempo de trânsito a partir da aplicação do critério de Coulomb para a ruptura da parede do poço submetida a um estado plano de tensões, a qual fornecerá a seção do poço, que é a forma geométrica que será traçada pela ferramenta de imageamento acústico do poço. O deslocamento ascensional da ferramenta e as imperfeições da parede do poço, normalmente, são os responsáveis pelo deslocamento da posição do transdutor em relação ao eixo do poço. Este efeito pode ter grande responsabilidade nas imperfeições das imagens de tempo de trânsito. Assim, a correção dessas imagens, chamada de correção da descentralização, busca o reposicionamento do transdutor para a posição do eixo do poço. Apresenta-se, também, um método de correção do efeito da descentralização da ferramenta baseado neste modelamento. O método é proposto com base na geometria analítica plana e no método do raio para a definição do tempo de trânsito do pulso acústico, com o objetivo de reconstruir as imagens de tempo de trânsito obtidas com a ferramenta descentralizada, ou seja, corrigir estas imagens tornando-as como se fossem obtidas com a ferramenta centralizada em relação ao eixo do poço.

Modelamento da permissividade dielétrica de rochas saturadas de óleo e água e suas aplicações em perfilagem de poço

A ferramenta de propagação eletromagnética (EPT) fornece o tempo de propagação (T_pl) e a atenuação (A) de uma onda eletromagnética que se propaga num meio com perdas. Estas respostas da EPT são funções da permissividade dielétrica do meio. Existem vários modelos e fórmulas de misturas sobre a permissividade dielétrica de rochas reservatório que podem ser utilizados na interpretação da ferramenta de alta frequência. No entanto, as fórmulas de mistura não consideram a distribuição e a geometria do espaço poroso, e estes parâmetros são essenciais para que sejam obtidas respostas dielétricas mais próximas de uma rocha real. Foi selecionado um modelo baseado nos parâmetros descritos acima e este foi aplicado à dados dielétricos disponíveis na literatura. Foi obtida uma boa concordância entre as curvas teóricas e os dados experimentais, comprovando assim que a distribuição e a geometria dos poros têm que ser levadas em conta no desenvolvimento de um modelo realista. Foram conseguidas também funções de distribuição de razão de aspecto de poros, através das quais geramos várias curvas relacionando as respostas da EPT com diversas saturações de óleo/gás. Estas curvas foram aplicadas na análise de perfis. Como o modelo selecionado ajusta-se bem aos dados dielétricos disponíveis na literatura, torna-se atraente aplicá-lo à dados experimentais obtidos em rochas de campos brasileiros produtores de hidrocarbonetos para interpretação da EPT corrida em poços destes campos petrolíferos.

Pesagem em movimento e caracterização do tráfego ferroviário com uso da técnica B-WIM

Nesta dissertação é apresentado o desenvolvimento de algoritmos para aplicação do método Bridge-Weigh In Motion (B-WIM) para a pesagem em movimento de trens e para a caracterização do tráfego ferroviário, permitindo-se obter informações sobre a velocidade de passagem dos trens, número e espaçamento entre eixos. Os sistemas B-WIM a partir de uma simples instrumentação permitem determinar as cargas por eixo de veículos em movimento, eliminando o efeito dinâmico. Foram implementados os algoritmos para a determinação dos valores referentes a geometria do trem e das cargas, que foi validado a partir de um exemplo teórico, onde se simulou a passagem de um trem de características conhecidas sobre a ponte e as cargas por eixos foram determinadas com 100% de exatidão. Além disso, foi feito um exemplo numérico em elementos finitos, de um viaduto em concreto armado para aplicação do método, onde foi feita a determinação das cargas por eixo para diferentes velocidades de passagem do trem. A fim de reduzir o tempo de processamento nas análises do exemplo numérico, foi desenvolvido um algoritmo para a geração de cargas nodais no modelo numérico que reduziram o tempo de processamento em até 96% quando comparado com a análise de múltiplos passos (“Multi-Step”), que simula automaticamente a passagem do trem sobre a estrutura. Finalmente, o método foi testado em um caso real a partir de monitorações realizadas em um viaduto de concreto armado da Estrada de Ferro Carajás. Apesar de não ter sido possível a determinação das cargas por eixo da locomotiva, foi possível medir precisamente o peso bruto total da locomotiva quando se utilizou o modelo constitutivo de Collins & Mitchell (1991) para o concreto.

Paleoambiente da formação mosquito e a implantação do sistema desértico úmido da formação corda, jurássico superior, Centro-Oeste da Bacia do Parnaíba

O Mesozóico foi marcado por mudanças geológicas significativas, decorrentes de soerguimentos resultante da orogenia Gonduanide, que possibilitou a implantação de sistemas desérticos concomitantemente com expressivos eventos magmáticos. Na Bacia do Parnaíba, Nordeste do Brasil, estes eventos estão registrados nas unidades siliciclásticas do Triássico, os arenitos da Formação Sambaíba, representadas pelos derrames basálticos e arenitos fluviais e eólicos subordinados da Formação Mosquito e pelos arenitos flúvio-eólicos da Formação Corda. O estudo de fácies e estratigráfico realizado em afloramentos e testemunhos de sondagem na região entre Formosa da Serra Negra e Montes Altos, Estado do Maranhão, possibilitou reconstituir o paleoambiente do topo da Formação Mosquito e da Formação Corda, e inferir condições paleoclimáticas para a porção centro-oeste da Bacia do Parnaíba durante o Jurássico. Foram identificadas vinte fácies sedimentares agrupadas em cinco associações de fácies (AF) representativas de uma planície vulcânica com depósitos fluviais esporádicos e arenitos eólicos subordinados (AF1-Formação Mosquito), sucedida pela instalação de um sistema desértico úmido (AF2-AF5; Formação Corda). A planície vulcânica (AF1) constitui derrames basálticos intercalados com arenitos finos a grossos (arenitos intertrap) compostos por grãos arredondados a subangulosos de quartzo, feldspatos e fragmentos de vidro vulcânico. Os arenitos apresentam estratificações plano-paralela e cruzada de baixo ângulo, preenchendo geometria de canal ou em corpos tabulares. Depósitos de canal fluvial entrelaçado (AF2) consistem em conglomerados polimíticos, com grânulos e seixos subarredondados a angulosos de basalto, e arenitos grossos com estratificação cruzada acanalada e acamamento maciço. Os lençóis arenosos (AF3) foram divididos em dois elementos arquiteturais (EA), o primeiro (EA1) consistem em arenitos finos a muitos com geometria tabular e estruturas de deformação, o segundo (EA2) é composto por arenito fino a grosso com estratificação cruzada acanalada e laminação cruzada cavalgante, gutter cast de pequeno porte. O campo de dunas (AF4) foi subdividido em dois conjuntos de fácies (C), o primeiro (CI) é caracterizado por arenitos com estratificações cruzadas tabular e tangencial de pequeno a médio porte, estratificação planoparalela e laminação cruzada cavalgante transladante subcrítica. O segundo (CII) consiste de arenitos finos a médios, moderadamente selecionados, laminação ondulada e estruturas de adesão e gretas de contração com rip-up clast, curled mud flakes, forma ciclos de raseamento centimétricos, com topo marcado por horizontes mosqueados, ricos em óxido/hidróxido de ferro, bioturbações e gretas de contração, interpretados como depósitos de interdunas úmidas. Os lobos de suspensão (AF5) consistem em arenitos finos intercalados com pelitos e arenito/pelito com estratificação cruzada complexa. A abundância de esmectita na AF4 aponta para condições de clima semiárido. No Jurássico, a região centro-oeste da Bacia do Parnaíba, foi submetida a movimentos distensivos com recorrência de derrames básicos advindos de fissuras na crosta. Durante os intervalos de aquiescência sedimentos de rios efêmeros preenchiam depressões ou espraiavam-se na planície vulcânica. O final da atividade magmática foi sucedido pela implantação do desérto Corda com campo de dunas e canais fluviais efêmeros (wadi) que retrabalharam parte da planície vulcânica e esporadicamente invadiam os lençóis arenosos. Comparado aos ergs do Permo-Triássico (Formação Sambaíba), o deserto Jurássico da Formação Corda foi mais úmido e menos extenso precedendo os sistemas fluviais e costeiros de clima mais ameno do Cretáceo da Bacia do Parnaíba.

Geologia, mineralogia e geoquímica dos fosfatos de Sapucaia (Bonito - PA)

O depósito mineral de Sapucaia, situado no município de Bonito, região nordeste do Estado do Pará, é parte de um conjunto de ocorrências de fosfatos de alumínio lateríticos localizados predominantemente ao longo da zona costeira dos estados do Pará e Maranhão. Estes depósitos foram alvos de estudo desde o início do século passado, quando as primeiras descrições de “bauxitas fosforosas” foram mencionadas na região NW do Maranhão. Nas últimas décadas, com o crescimento acentuado da demanda por produtos fertilizantes pelo mercado agrícola mundial, diversos projetos de exploração mineral foram iniciados ou tiveram seus recursos ampliados no território brasileiro, dentre estes destaca-se a viabilização econômica de depósitos de fosfatos aluminosos, como o de Sapucaia, que vem a ser o primeiro projeto econômico mineral de produção e comercialização de termofosfatos do Brasil. Este trabalho teve como principal objetivo caracterizar a geologia, a constituição mineralógica e a geoquímica do perfil laterítico alumino-fosfático do morro Sapucaia. A macrorregião abrange terrenos dominados em sua maioria por rochas pré-cambrianas a paleozóicas, localmente definidas pela Formação Pirabas, Formação Barreiras, Latossolos e sedimentos recentes. A morfologia do depósito é caracterizada por um discreto morrote alongado que apresenta suaves e contínuos declives em suas bordas, e que tornam raras as exposições naturais dos horizontes do perfil laterítico. Desta forma, a metodologia aplicada para a caracterização do depósito tomou como base o programa de pesquisa geológica executada pela Fosfatar Mineração, até então detentora dos respectivos direitos minerais, onde foram disponibilizadas duas trincheiras e amostras de 8 testemunhos de sondagem. A amostragem limitou-se à extensão litológica do perfil laterítico, com a seleção de 44 amostras em intervalos médios de 1m, e que foram submetidas a uma rota de preparação e análise em laboratório. Em consonância com as demais ocorrências da região do Gurupi, os fosfatos de Sapucaia constituem um horizonte individualizado, de geometria predominantemente tabular, denominado simplesmente de horizonte de fosfatos de alumínio ou crosta aluminofosfática, que varia texturalmente de maciça a cavernosa, porosa a microporosa, que para o topo grada para uma crosta ferroalumino fosfática, tipo pele-de-onça, compacta a cavernosa, composta por nódulos de hematita e/ou goethita cimentados por fosfatos de alumínio, com características similares aos do horizonte de fosfatos subjacente. A crosta aluminofosfática, para a base do perfil, grada para um espesso horizonte argiloso caulinítico com níveis arenosos, que repousa sobre sedimentos heterolíticos intemperizados de granulação fina, aspecto argiloso, por vezes sericítico, intercalados por horizontes arenosos, e que não possuem correlação aparente com as demais rochas aflorantes da geologia na região. Aproximadamente 40% da superfície do morro é encoberta por colúvio composto por fragmentos mineralizados da crosta e por sedimentos arenosos da Formação Barreiras. Na crosta, os fosfatos de alumínio estão representados predominantemente pelo subgrupo da crandallita: i) série crandallita-goyazita (média de 57,3%); ii) woodhouseíta-svanbergita (média de 15,8%); e pela iii) wardita-millisita (média de 5,1%). Associados aos fosfatos encontram-se hematita, goethita, quartzo, caulinita, muscovita e anatásio, com volumes que variam segundo o horizonte laterítico correspondente. Como os minerais pesados em nível acessório a raro estão zircão, estaurolita, turmalina, anatásio, andalusita e silimanita. O horizonte de fosfatos, bem como a crosta ferroalumínio-fosfática, mostra-se claramente rica em P₂O₅, além de Fe₂O₃, CaO, Na₂O, SrO, SO₃, Th, Ta e em terras-raras leves como La e Ce em relação ao horizonte saprolítico. Os teores de SiO₂ são consideravelmente elevados, porém muito inferiores aqueles identificados no horizonte argiloso sotoposto. No perfil como um todo, observa-se uma correlação inversa entre SiO₂ e Al₂O₃; entre Al₂O₃ e Fe₂O₃, e positiva entre SiO₂ e Fe₂O₃, que ratificam a natureza laterítica do perfil. Diferente do que é esperado para lateritos bauxíticos, os teores de P₂O₅, CaO, Na₂O, SrO e SO₃ são fortemente elevados, concentrações consideradas típicas de depósitos de fosfatos de alumínio ricos em crandallita-goyazita e woodhouseítasvanbergita. A sucessão dos horizontes, sua composição mineralógica, e os padrões geoquímicos permitem correlacionar o presente depósito com os demais fosfatos de alumínio da região, mais especificamente Jandiá (Pará) e Trauíra (Maranhão), bem como outros situados além do território brasileiro, indicando portanto, que os fosfatos de alumínio de Sapucaia são produtos da gênese de um perfil laterítico maturo e completo, cuja rocha fonte pode estar relacionada a rochas mineralizadas em fósforo, tais como as observadas na Formação Pimenteiras, parcialmente aflorante na borda da Bacia do Parnaíba. Possivelmente, o atual corpo de minério integrou a paleocosta do mar de Pirabas, uma vez que furos de sondagem às proximidades do corpo deixaram claro a relação de contato lateral entre estas unidades.

Depósitos pleistocenos da formação Itaubal: paleoambiente e implicações na evolução da planície costeira do Amapá

A costa norte da América do Sul durante o Pleistoceno tardio esteve sujeita a oscilações do nível do mar, relacionadas a variações climáticas e influência da dinâmica de sedimentos carreados pelo Rio Amazonas, que moldaram a paisagem desta região. Terraços Pleistocenos da Formação Itaubal, anteriormente considerados como pertencentes à Formação Barreiras (Mioceno), constituem parte da Planície Costeira do Amapá e recobrem rochas do Escudo das Guianas. A integração das análises de fácies, estratigráfica e datações por Luminescência Opticamente Estimulada / regeneração de alíquota única e múltipla (LOE / SAR-MAR) entre 120.600 (± 12.000) e 23.150 (±6.800) anos AP permitiu o posicionamento da Formação Itaubal no Pleistoceno Superior. Estes depósitos siliciclásticos, de espessura máxima de 10 m, cor amarronzada a avermelhada com camadas de geometria tabular foram divididos em duas unidades separadas por inconformidade. A Unidade Inferior compreende as associações de fácies, de planície de inframaré (AF1) e de canal fluvial meandrante influenciado por maré (AF2), enquanto que a Unidade Superior, com maior concentração de argila que a Unidade Inferior, consiste em depósitos de planície de maré (AF3) e de canal fluvial entrelaçado (AF4). As duas unidades têm caracteristicas progradacionais dentro de trato de sistema de mar alto e regressivo, e foram depositadas diretamente sobre rochas do embasamento intensamente intemperizadas durante o Mioceno-Pleistoceno. Os depósitos da Formação Itaubal foram expostos durante o Último Máximo Glacial (22.000 - 18.000 anos AP) e posteriormente sobrepostos por depósitos finos do Rio Amazonas, que configuram a atual linha de costa da Costa norte da América do Sul. Pela primeira vez a Formação Itaubal define os eventos sedimentares do Pleistoceno na evolução da Planície Costeira do Amapá. A correlação de seus depósitos com os de Suriname e nordeste do Pará amplia a discussão sobre a configuração da linha de costa do norte da América do Sul desde o Pleistoceno.