Estudo geofísico regional sobre águas subterrâneas na Ilha do Marajó-Pará-Brasil

Nos últimos dez anos foram realizadas na parte leste da Ilha de Marajó (região dos campos naturais) pelo IDESP e NCGG, mais de 800 SEVs para fins hidrogeológicos. Na época, grande parte dessas SEVs não foram totalmente interpretadas em forma quantitativa, devido à falta de recursos técnicos para fazê-lo de forma eficiente. Agora, usando meios mais modernos para interpretação automática de SEVs, voltou-se a interpretá-las com a finalidade de apresentar uma visão regional dos principais aquíferos da área, agrupar as SEVs em famílias características, testar até que ponto essa interpretação é confiável e propor o modelamento bidimensional como técnica alternativa para interpretar as SEVs realizadas em certos locais da área em questão. Como resultado dessa interpretação, com base na teoria convencional dos meios estratificados, foram definidos três tipos de sistemas de aquíferos. 1. O primeiro, denominado de aquífero profundo, situado a profundidades maiores que 50m, estende-se por toda a região prospectada, estando provavelmente associada às camadas superiores da Formação Marajó ou às litologias altamente resistivas das camadas mais profundas do Grupo Pará. 2. O segundo, denominado de aquífero raso e de média profundidade, localiza-se na parte sul e sudeste da região a profundidades compreendidas entre 10 a 50m, e está associado às lentes arenosas do Grupo Pará. 3. O terceiro, é constituído pelos paleocanais e estruturas similares, distribuídos aleatoriamente na região a pouca profundidade. A partir do estudo detalhado das SEVs, decidiu-se classificá-las em 3 famílias características com seus respectivos tipos e apresentar mapas de localização e da espessura dos aquíferos, bem como mapas de condutância longitudinal total e resistividade média da área. Estes últimos, permitem que se divida a região dos campos da Ilha de Marajó em três zonas principais: 1. Uma, altamente resistiva, situada ao sul e sudeste, a qual coincide com os terrenos aflorantes do Grupo Pará. 2. Outra, altamente condutiva, está localizada no centro e norte, onde se encontram aleatoriamente distribuídos os paleocanais e coincide com os terrenos topograficamente mais baixos, geralmente argilosos e embebidos de água salgada, que são procedentes da erosão dos terrenos circundantes topograficamente mais altos. 3. A última é medianamente resistiva e está relacionada com os terrenos vizinhos à cidade de Chaves (noroeste da região dos campos), os quais apresentam semelhanças com os do sul e sudeste da área. Usando-se a técnica de inversão na interpretação de uma SEV característica de cada família, testou-se, através do seu tratamento estatístico, até que ponto os modelos usados na interpretação dessas SEVs (teoria convencional dos meios estratificados) seriam confiáveis. Conclui-se, então, que a alta correlação existente entre os parâmetros dos modelos assumidos (camadas horizontais, isotrópicas e homogêneas) pode-se dever à utilização de modelos geofísicos muito simples para interpretar a complexa geologia de Marajó. Tendo-se verificado que nem sempre é possível aplicar a teoria das SEVs em meios horizontalmente estratificados para interpretar SEVs obtidas em certos locais de Marajó, os quais muitas vezes apresentam bruscas variações laterais de resistividade, passou-se a demonstrar que estas variações laterais afetam profundamente os dados das SEVs, utilizando-se para isto a técnica dos elementos finitos, a qual leva em conta essa variação bidimensional das propriedades físicas do meio. Foi também possível com esta técnica, modelar uma estrutura rasa, semelhante a um paleocanal, concluindo-se que estes resultados sugerem o emprego, duma forma mais profunda, deste tipo de tratamento para os dados obtidos na região dos campos da Ilha de Marajó.

Estudo da hidrogeologia no município de Palestina do Pará usando o método de eletroresistividade

Esse trabalho de geofísica aplicada à hidrogeologia, foi realizado numa área localizada na porção sudeste do estado do Pará, mas precisamente na cidade de Palestina à margem esquerda do rio Araguaia. Ironicamente existem problemas de abastecimento de água potável para a população da cidade. Esse problema deve-se ao fato que, no momento, é antieconômico o tratamento e transporte de água do rio para a cidade. Considerando que a extração de água subterrânea é economicamente mais viável para a solução desse problema, o Departamento de Geofísica e o curso de Pós-graduação de Geofísica da Universidade Federal do Pará (UFPa), por solicitação da prefeitura daquele município e da Fundação Nacional de Saúde (FNS), participaram no levantamento geofísico de eletroresistividade para investigar os possíveis locais de ocorrência de água subterrânea. Foram realizadas, na área, 21 Sondagens Elétricas Verticais (SEVs) medidas na superfície, utilizando o arranjo Schlumberger com abertura máxima dos eletrodos AB e MN de 500m e 50m, respectivamente. Para a sua interpretação foram usados os programas de modelagem direta (EGSLIB/SEV1D) e inversa (EGSLIB/SEV1DIN) de RIJO (1994), de modelos de n camadas horizontais, homogêneas e isotrópicas. Com base nos resultados da interpretação das 21 SEVs, e também das informações da geologia local foram confeccionadas três mapas de contorno e seis secções geoelétricas. Dos seis horizontes que aparecem nessas secções, o mais favorável para a exploração hidrogeológica é o quarto, com espessura variando de 10 a 79 m, estando a uma profundidade média de 50m e com valores de resistividade variando de 150 — 850 Ωm, sendo correlacionavel com a Formação Potí. A analise final dos mapas e das secções geoelétricas, indicou que o local mais apropriado para a perfuração de poços, é uma área próximo à atual área de captação.

Identificação de fácies em perfis de poço com algoritmo inteligente

A identificação de fácies em um poço não testemunhado é um dos problemas clássicos da avaliação de formação. Neste trabalho este problema é tratado em dois passos, no primeiro produz-se a codificação da informação geológica ou da descrição das fácies atravessadas em um poço testemunhado em termos das suas propriedades físicas registradas nos perfis geofísicos e traduzidas pelos parâmetros L e K, que são obtidos a partir dos perfis de porosidade (densidade, sônico e porosidade neutrônica) e pela argilosidade (Vsh) calculada pelo perfil de raio gama natural. Estes três parâmetros são convenientemente representados na forma do Gráfico Vsh-L-K. No segundo passo é realizada a interpretação computacional do Gráfico Vsh-L-K por um algoritmo inteligente construído com base na rede neural competitiva angular generalizada, que é especializada na classificação de padrões angulares ou agrupamento de pontos no espaço n-dimensional que possuem uma envoltória aproximadamente elipsoidal. Os parâmetros operacionais do algoritmo inteligente, como a arquitetura da rede neural e pesos sinápticos são obtidos em um Gráfico Vsh-L-K, construído e interpretado com as informações de um poço testemunhado. Assim, a aplicação deste algoritmo inteligente é capaz de identificar e classificar as camadas presentes em um poço não testemunhado, em termos das fácies identificadas no poço testemunhado ou em termos do mineral principal, quando ausentes no poço testemunhado. Esta metodologia é apresentada com dados sintéticos e com perfis de poços testemunhados do Campo de Namorado, na Bacia de Campos, localizada na plataforma continental do Rio de Janeiro, Brasil.

Fácies e proveniência de depósitos costeiros da Formação Raizama: evidências do registro Ediacarano-cambriano na faixa Paraguai, região de Nobres, Mato Grosso

Rochas siliciclásticas da Formação Raizama, unidade basal do Grupo Alto Paraguai de idade ediacarana-cambriana (635 – 541 Ma), ocorrem distribuídas descontinuamente ao longo da margem sul do Cráton Amazônico e segmento norte da Faixa Paraguai, centro-oeste do Brasil, estado do Mato Grosso. Estas recobrem discordantemente os depósitos de plataforma carbonática do Grupo Araras, onde foram registrados evidências do evento glacial Marinoano (635 Ma). O Grupo Alto Paraguai representa os estágios finais da colisão entre os blocos Amazônico e Paranapanema que culminaram no fechamento do Oceano Clymene (540-520 Ma). A Formação Raizama com espessura de aproximadamente 570 m é constituída por pelitos, arenitos finos a grossos, e arenitos com cimento dolomítico previamente interpretados como depósitos flúvio-costeiros distribuídos nos membros inferior (270 m) e superior (300 m). O estudo faciológico e estratigráfico desta unidade na região de Nobres, Estado do Mato Grosso, foi focado principalmente na seção aflorante de 600 m no leito do rio Serragem, que inclui a Cachoeira da Serra do Tombador. Foram definidas 17 fácies sedimentares, agrupadas em cinco associações de fácies (AF) representativas de uma sucessão costeira progradante iniciando por depósitos de shoreface inferior, os quais recobrem em conformidade correlativa os depósitos de plataforma carbonática da Formação Serra do Quilombo (Grupo Araras). A AF1 consiste em arenitos com laminação plano-paralela e laminação truncada por onda (microhummocky), individualizada por camadas de pelito laminado interpretadas como depósitos de shoreface inferior. Destaca-se na AF1 a primeira ocorrência de níveis centimétricos bioturbados por Skolithos em depósitos neoproterozoicos – cambrianos na Faixa Paraguai. A AF2 é formada por arenitos com estratificação cruzada swaley e estratificação plano-paralela interpretada como depósitos de shoreface superior. A AF3 é composta por arenitos com estratificações cruzadas tangenciais e acanaladas com recobrimentos de siltito/arenito muito fino representativos de depósitos de canal e barras de submaré. A AF4 é caracterizada por arenitos com estratificações cruzadas tangencial e sigmoidal, laminação plano-paralela a cruzadas de baixo-ângulo, ritmito arenito muito fino/siltito com acamamento flaser e gretas de contração, organizados em ciclos métricos de raseamento ascendente de planície de maré. A AF5 é constituída por arenito com estratificação cruzada acanalada marcada por lags residuais na base da associação, arenito com estratificações plano-paralela e cruzada de baixo-ângulo, interpretados como depósitos fluviais distais de rios entrelaçados, parcialmente retrabalhados por ondas. Grãos detríticos de zircão foram obtidos da AF3 e datados pelo método U-Pb, sendo a idade de 1001±9 Ma interpretada como a idade de máxima deposição da Formação Raizama. Aliado a tal análise, as paleocorrentes NE-SE mostram que estes grãos teriam como áreas fontes principais a Faixa Sunsás, SW do Cráton Amazônico, não sendo descartada contribuições oriundas da parte NW desse Cráton. A idade mesoproterozóica obtida serviu principalmente para desvendar a proveniência da Formação Raizama, enquanto que as datações da base do Grupo Araras, em torno de 627-622 Ma, associada à presença inequívoca do icnogênero Skolithos, tornam esta unidade muito mais próxima do limite com o Cambriano Inferior. Traços fósseis do Proterozoico são caracterizados quase que exclusivamente por traços horizontais, sendo que bioturbações verticais praticamente são ausentes ao longo do Neoproterozoico. Esta inferência vem de encontro com a idade máxima de 541 Ma obtida para a Formação Diamantino, a qual recobre a unidade estudada. Os dados radiométricos aliados com as interpretações paleoambientais, que incluem o registro das primeiras atividades de organismos perfurantes na Faixa Paraguai, abrem perspectivas de entender com maiores detalhes a sequência de eventos que tipificam os estratos do limite Ediacarano-Cambriano do Brasil, ainda pouco conhecidos.

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

Métodos geofísicos aplicados à prospecção de cobre na Serras de Carajás-Pará

Através de um convênio entre o Núcleo de Ciências Geofísicas e Geológicas da UFPa e a Companhia Rio Doce Geologia e Mineração S.A., um trabalho geofísico foi realizado na área denomina da ALVO 1, MM1, Serra dos Carajás, sul do Pará, Brasil. Foram aplicados os métodos magnético, polarização elétrica induzida e AFMAG, visando determinar a existência e a distribuição de sulfetos de cobre em subsuperfície. O trabalho foi realizado em uma malha de reconhecimento de 9 km² (ALVO 1) e em outra de detalhe com aproximadamente 0.63 km² (Corpo 2). O método AFMAG detectou vários condutores deslocados entre si por prováveis falhamentos. A orientação destes condutores, coincide com o "strike" dominante do xisto, o qual contém sulfetos de cobre associado a uma formação ferrífera. Pela estreita correlação que guardam com as anomalias magnéticas é possível supor-se que as anomalias de AFMAG são provocadas primeiramente pela formação ferrífera e secundariamente pelos sulfetos. O método magnético detectou fortes anomalias sobre a formação ferrífera. O deslocamento dos segmentos de anomalias implica em falhamento. Este método é muito importante na área, pela relação mantida entre a zona dos sulfetos e a formação ferrífera. O método polarização elétrica induzida/resistividade detectou anomalias com EPF variando de 6 a 12% aproximadamente. Várias destas anomalias foram confirmadas com furos de sondagem e trabalho de laboratório (Rocha, 1979). Verificou-se que essas anomalias provavelmente são provocadas por efeito combinado da formação ferrífera magnética e do sulfeto. A pesquisa geofísica realizada nesta área piloto tem grande importância porque os resultados obtidos poderão ser estendidos a uma área maior, em condições geológicas semelhantes.

Mangrove sedimentary characteristics and implications for crab Ucides cordatus (Crustacea, Decapoda, Ucididae) distribution in an estuarine area of the Amazonian region

Na zona costeira, sedimentos, água e organismos interagem intensamente. Nas costas equatoriais dominadas por maré os manguezais são abundantes. Estas áreas são conhecidas por sua importância ecológica. No caso dos manguezais da costa atlântica da América do Sul o caranguejo-uçá Ucides cordatus (Linnaeus, 1763) se destaca por sua relevância ecológica e econômica, sendo que altas densidades deste organismo são encontradas na zona costeira amazônica. O presente estudo investiga a distribuição de sedimentos nos manguezais de Bragança (costa Amazônica, Brasil) e suas correlações com a vegetação e a distribuição do caranguejo-uçá. Quarenta e sete amostras de sedimento foram avaliadas, assim como caranguejos de 8 destas áreas foram coletados, onde o tipo dominante de vegetação foi também identificado. Os resultados demonstram que os sedimentos superficiais, assim como no extrato 0,8 a 1 m de profundidade, na área são principalmente compostos por silte (59%), incluindo em média 21% de areia e 20% de argila. O tipo de vegetação predominante foi significativamente correlacionado com a abundância e tamanho/peso dos caranguejos. As características sedimentares também foram substancialmente diferentes dependendo da vegetação. Áreas dominadas por Avicennia germinans tiveram mais areia e argila que as áreas dominadas por Rizophora mangle, onde a fração silte prevalece grandemente e os caranguejos eram significativamente maiores e mais abundantes. Os resultados demonstraram que sedimentos, invertebrados bentônicos e vegetação estão intimamente relacionados nos manguezais e devem ser estudados de maneira integrada.

Modelagem de dados MT 2D em multicamadas com anisotropia arbitrária

Neste trabalho apresentamos a modelagem direta bidimensional do método Magnetotelúrico (MT) com anisotropia arbitrária na condutividade elétrica para modelos de multicamadas estratificadas horizontalmente, separando a resposta do campo resultante em duas partes, primária e secundária. A parte primária refere-se ao campo 1D do modelo na ausência de heterogeneidades; a secundária, à contribuição no campo resultante gerada pela anomalia. Esta modelagem foi feita via técnica numérica dos Elementos Finitos (EF). Apresentamos também a modelagem direta do caso MT 1D para o mesmo modelo e anisotropia arbitrária, porém, com solução analítica pelo método da matriz de propagação; tal resultado é usado como fonte na formulação 2D. No estudo a respeito do comportamento dos campos 1D em profundidade na presença da anisotropia, verificamos nas componentes horizontais que suas curvas dentro de cada camada se assemelham a de um campo na presença de fontes, pois experimentam um aumento na amplitude; além de verificar o surgimento da componente normal E_z do campo, sem que haja uma densidade de corrente nesta direção. Com arbitrariedade na anisotropia queremos dizer que seus eixos principais podem assumir qualquer orientação em relação ao eixo de referência adotado para o problema; tal orientação é descrita sob três rotações consecutivas, chamadas de strike, dip e slant, que simulam inclinações bastante presentes na geologia dos ambientes terrestres. Verificamos que a anisotropia provoca alterações no comportamento das pseudo-seções de resistividades aparentes, interferindo na interpretação correta dos dados.

Modelagem sísmica por diferenças finitas em meios bidimensionais com difratores

Este trabalho tem por objetivo a modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações. Como parte deste estudo foi feita a análise numérica visando as condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas. Para a validação do programa de diferenças finitas foi feita a comparação cinemática com a Teoria do Raio para um modelo simples. O motivo deste estudo é ter uma melhor compreensão e controle sobre os problemas de modelagem, visando contribuir para a solução de problemas na interpretação de dados sísmicos. Segundo vários autores na literatura geológica, Derby (1877), Evans (1906), Paiva (1929) e Moura (1938). A Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o recente, atingindo espessuras da ordem de 5 km. Os corpos de diabásio, presentes entre os sedimentos paleozóicos, estão dispostos na forma de soleiras, alcançando espessuras de centenas de metros, perfazendo um volume total de 90.000 quilômetros cúbicos. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica, impossibilitando uma melhor interpretação dos horizontes refletores que se encontram abaixo das soleiras. Para representar situações geológicas desse tipo foi usado um modelo acústico de velocidades. Para o cálculo dos sismogramas foi utilizado um programa de diferenças finitas com aproximação de quarta ordem da equação da onda acústica no espaço e no tempo. As aplicações dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas têm melhorado a compreensão sobre a propagação em meios onde existem heterogeneidades significativas, tendo como resultado boa resolução na interpretação dos eventos de reflexão sísmica em áreas de interesse. Como resultado dos experimentos numéricos realizados em meio de geologia complexa, foi observado a influência significativa das múltiplas devido a camada de alta velocidade, o que faz com que haja maior perda de energia dificultando a interpretação dos alvos. Por esta razão recomendo a integração de dados de superfície com os dados de poço, com o objetivo de se obter uma melhor imagem dos alvos abaixo das soleiras de diabásio.

A capa carbonática do sudoeste do cráton amazônico, estado de Rondônia: nova ocorrência e extensão dos eventos pós-glaciação marinoana (635 Ma)

Na porção oeste da Bacia dos Parecis, Estado de Rondônia, inserida no sudoeste do Cráton Amazônico, rochas carbonáticas expostas nas bordas dos grábens Pimenta Bueno e Colorado têm sido consideradas como parte do preenchimento eopaleozoico da bacia. A avaliação das fácies/microfácies e quimioestratigrafia dessas rochas nas regiões de Chupinguaia e Pimenta Bueno, confirmou a ocorrência de dolomitos rosados que sobrepõem, em contato direto, diamictitos glaciais previamente interpretados como depósitos de leques aluviais. Trabalhos prévios reportaram excursão negativa de δ¹³C, também confirmados neste trabalho, com variações entre -4.6 e -3,8‰VPDB em Chupinguaia e média de - 3,15‰VPDB em Pimenta Bueno. Esse padrão, de sedimentação e quimioestratigráfico, ausente nas rochas paleozoicas, é comumente encontrado nos depósitos carbonáticos anômalos do Neoproterozoico. No sul do Cráton Amazônico, Estado do Mato Grosso, rochas com essas mesmas características são descritas como capas carbonáticas relacionadas à glaciação marinoana (635 Ma). Neste trabalho, consideramos que os dolomitos rosados sobre diamictitos, em Rondônia, fazem parte do mesmo contexto das capas carbonáticas encontradas no Mato Grosso. Adicionalmente, destaca-se o contato brusco e deformado do dolomito sobre o diamictito, presente em ambas as ocorrências, configurando-se uma das feições típicas das capas carbonáticas do Cráton Amazônico. Essa relação paradoxal, entre diamictito e dolomito, tem sido interpretada como produto da mudança rápida das condições atmosféricas de icehouse para greenhouse, e a deformação da base foi gerada pelo rebound isostático. A capa carbonática de Rondônia compreende duas associações de fácies (AF2 e AF3) que recobrem depósitos glacio-marinhos compostos por paraconglomerados polimíticos (Pp), e arenito seixoso laminado (Asl), da AF1. A AF2 consiste em dolomudstone/dolopackstone peloidal com laminação plana a quasi-planar e com truncamentos de baixo-ângulo (fácies Dp), megamarcas onduladas (fácies Dm) e laminações truncadas por ondas (fácies Dt), interpretada como depósitos de plataforma rasa influenciada por ondas. Esta sucessão costeira é sucedida pela AF3, que compreende as fácies: dolomudstone/dolopackstone e dolomudstone/dolograinstone com partição de folhelho (Df) e siltito laminado (Sl). A fácies Df compreende um pacote de 6 metros de dolomito com partição de folhelho, apresentando lâminas de calcita fibrosa (pseudomorfos de evaporito) e dolomitos com laminações onduladas de corrente. Sobrejacente à fácies Df, ocorre a fácies Sl, apresentando 5 metros de siltito argiloso com laminação plana. Esta associação é interpretada como depósitos de plataforma rasa influenciada por maré, sendo sobreposta discordantemente, em contato angular, por depósitos glaciais do Eopaleozoico. Os valores isotópicos de C e O são negativos e refletem o sinal primário do C. No entanto, pode-se considerar uma leve influência da diagênese meteórica no sinal. As principais quebras nos sinais negativos podem estar associadas à influência meteórica, expressa pela substituição e preenchimento de poros por calcita e pela proximidade de superfícies estratigráficas, os quais refletem alguns padrões de alteração diagenética, representados nos sinais mais negativos. Diferentemente da capa carbonática do Mato Grosso, a capa de Rondônia possui níveis de pseudomorfos de evaporito e dolomitos com partição de folhelho (ritmito), em sucessão de fácies marinha rasa, onde os dolomitos de plataforma rasa influenciada por ondas passam para ritmitos e siltitos de plataforma rasa influenciada por maré (zona de inframaré), configurando uma sucessão retrogradante. Esta nova ocorrência de capa carbonática modifica a estratigrafia da base da Bacia dos Parecis, ao passo que exclui essas rochas carbonáticas da sequência eopaleozoica. Além disso, fornece informações que permitem reconstruir melhor a paleogeografia costeira da bacia neoproterozoica que acumulou os depósitos da plataforma carbonática do Grupo Araras, bem como estende os eventos pós-marinoanos ligados à hipótese do Snowball/Slushball Earth para o sudoeste do Cráton Amazônico, exposto no Estado de Rondônia.

Síntese e caracterização de materiais híbridos (PCHs e Organoargilas) com aplicabilidade na adsorção de fenol, benzeno e tolueno

Este estudo avalia o uso de argila do Estado do Acre de natureza esmectítica (argila original) como precursora na síntese de PCHs (Porous Clay Heterostructure) e na organofuncionalização com mercaptopropiltrimetoxisilano (MPTS-Arg) e aplicabilidade desses materiais como adsorventes de fenol, benzeno e tolueno em fase líquida. A síntese de PCHs envolve as seguintes etapas: preparação de organoargila com hexadeciltrimetilamônio (HDTMA-Argila); aumento do espaçamento basal (d₀₀₁) da organoargila usando o cosurfactante octilamina e adição concomitante de uma fonte de Si, como o tetraetilortosilicato (TEOS); lavagens com etanol em HCl, secagem a 60 o C (PCH-60 A e B); calcinação a 500 o C e 700 º C (PCH 500 A/B e PCH 700), com a remoção do surfatante/co-surfactante e consequente substituição destas substâncias orgânicas por prótons que permanecem nas paredes das galerias de pilares de SiO₂ até novo tratamento porventura realizado no PCH. O material MPTS-Arg foi obtido a partir da intercalação de mercaptopropiltrimetoxisilano na argila original sob agitação durante 1 h (Arg-MPTS 1) e na argila original previamente tratada com HCl 1 mol L^-1 e sob agitação durante 2 h (Arg-MPOS 2). Os materiais (argila original, HDTMA-Argila, PCH-60, PCH-500/700 e MPTS-Arg) foram caracterizados por difração de raios-X (DRX), microscopia eletrônica de varredura acoplado a espectroscopia de raios X por dispersão de energia (MEV/ EDS), espectroscopia de absorção molecular IV com transformada de Fourier (FTIR), técnicas de análises térmicas (DTA_TG) e análises texturais, como as medidas de área superficial específica (ASEBET, ASE_Langmuir), volume total de poros (VTP) e diâmetro médio de poros (DMP). Caracterização complementar por ressonância nuclear magnética no estado sólido acoplada a espectrometria de massa (RMN-M AS) foi realizada para o PCH-500 A. Nos experimentos de adsorção de fenol, benzeno e tolueno foram utilizadas suspensões aquosas dos adsorventes contendo separadamente os adsorvatos. As concentrações de equilíbrio de fenol, tolueno e benzeno foram medidas por espectroscopia de absorção molecular na região ultravioleta, conforme método aceito e aplicado em outros estudos de adsorção. Foram obtidos os seguintes resultados: identificação do argilomineral esmectita na argila natural com espaçamento basal de 1,25nm; intercalação do HDTMA e do MPTS com expansão do d(001) para 1,95nm (HDTMA-Arg 5) e 2,1nm (MPTS-Argila 1h). Após a modificação observou-se através da análise de EDS um decréscimo das concentrações dos elementos Na, Si e Al (%) nas posições interlamelares e aumento na concentração de C. O espectro ²⁹Si MAS NMR do PCH-500 A apresentou sinais químicos dos centro Q₃ Si(OSi)₃OH e Q₄ Si(OSi)₄ que indicam a formação de sílica em sua estrutura. Os PCHs (60 e 500 A) foram classificados como mesoporosos, apresentando ASE_BET entre 417 a 446 m² g^-1 e volume total de poros (cm³ g^-1) = 0,367 - 0,369. O aumento da temperatura de calcinação do PCH 500 para 700 o C (PCH 700) propiciou decréscimos em SBET (305 m² g^-1) e diâmetro médio de poros (0.976 nm). Foi verificado a partir da avaliação das medidas dos diâmetro de poros em PCH 700 que 46% das medidas correspondem a microporos. No estudo do equilíbrio de adsorção foram obtidos os valores de K_L (L mg^-1), K_F (L mg^-1) e q_Max (mg g^-1) nos seguintes processos: (i) adsorção de fenol (1.40; 2.26 e 3.80) e tolueno (145, 7.18 e 135) em PCH-500 A; (ii) adsorção de benzeno (0.33, 41.49 e 69.81) em PCH 500 B; (iii) na adsorção de tolueno em argila-MPOS 1 (1.06, 118.9 e 146) e argila-MPOS 2 (2.1, 92 e 97); (iv) na adsorção de fenol (0.07, 0.48 e 4.44) e benzeno (0.11, 4.95 e 35.8) em HDTA-Arg 5, respectivamente. Os dados obtidos sobre a caracterização dos materiais indicaram que a formação dos PCHs propiciou grande aumento nas concentrações de Si e decréscimo nas concentrações de Al, Fe e outros elementos, devido à delaminação- exfoliação da esmectita original. Foi verificado que os materiais são mais adequados para adsorção de tolueno e benzeno do que fenol.

Morfodinâmica costeira e o uso da orla oceânica de Salinópolis (Nordeste do Pará)

Os aspectos morfodinâmicos relacionados à erosão ou acresção da linha de costa são alguns dos assuntos analisados na gestão das zonas costeiras que vêem sendo tratada em todo mundo no sentido de monitorar e proteger essas zonas. Esta tese objetiva analisar o comportamento da morfodinâmica costeira de Salinópolis, relacionando-o ao uso da orla oceânica. A área de estudo foi compartimentada em três setores: Oeste (praias da Corvina e do Maçarico), Central (praia do Farol Velho) e Leste (praia do Atalaia). A metodologia consistiu na: (a) aquisição e tratamento de imagens multitemporais (1988-2001-2013) do satélite Landsat 5 TM, 7 ETM e 8 OLI; (b) aplicação de entrevistas/questionários com banhistas, (c) aquisição de dados de campo durante as estações chuvosa (26, 27, 28/04/2013) e menos chuvosa (04, 05, 06/10/2013); e (d) análise laboratorial para o tratamento dos dados adquiridos em campo (topografia das praias estudadas, amostragem de sedimentos superficiais das mesmas e com o uso de armadilhas, e medições oceanográficas de ondas, marés, correntes e turbidez). Foram feitas as representações gráficas dos perfis topográficos das praias, calculados os parâmetros estatísticos granulométricos de Folk & Ward (1957), as taxas do transporte sedimentar nas praias e os parâmetros morfométricos de Short & Hesp (1982), estes últimos foram calculados com o intuito de relacioná-los aos estados morfodinâmicos de praias propostos por Wright & Short (1984) e Masselink & Short (1993). Para a classificação da costa oceânica de Salinópolis em termos de uso e ocupação foi utilizado o decreto nº 5.300 de 7 de dezembro de 2004. A partir das pesquisas sobre a urbanização na costa e das obras situadas nos ambientes costeiros foi utilizada uma matriz proposta por Farinaccio & Tessler (2010) que lista uma série de impactos ambientais, e o quadro de geoindicadores do comportamento da linha de costa proposto por Bush et al. (1999), para a identificação de locais com vulnerabilidade à erosão ou acresção. Para as condições oceanográficas em cada praia e periculosidade ao banho nas mesmas, foram integralizados os dados de ondas, de correntes, de morfodinâmica praial e questionários aplicados com banhistas. Atualmente, a orla oceânica de Salinópolis possui diferentes características quanto à utilização e conservação, abrangendo desde a tipologia de orlas naturais (Classe A) até orlas com urbanização consolidada (Classe C). A primeira ocorre nos extremos da área de estudo e, a segunda, na região da sede municipal. Quatro tipos de praias foram identificados segundo a exposição marítima e o grau das condições oceanográficas: tipo 1 (Maçarico), tipo 2 (Corvina), tipo 3 (Farol Velho) e tipo 4 (Atalaia). O trecho de costa com maiores impactos ambientais e com elevada erosão costeira localiza-se na praia do Farol Velho. O grau de periculosidade ao banho foi de 4 (praia do Maçarico) a 7 (praia do Atalaia) – médio a alto grau de risco. As praias de Salinópolis apresentam declives suaves (< 1,5°), grandes variações na linha de costa entre as estações do ano (9,6 a 88, 4 m) e volume sedimentar variável dependendo do grau de exposição das praias ao oceano aberto. Predominou o estado morfodinâmico dissipativo (Ω>5,5) para estas praias, mas com ocorrência do estado de banco e calha longitudinais (4,7<Ω<5,5) no setor oeste. As macromarés na área de estudo apresentaram altura máxima de 5,3 m (Setor Central, durante a estação menos chuvosa) e mínima de 4 m no mesmo setor, durante a estação chuvosa. As correntes longitudinais foram mais intensas no setor leste (>0,45 m/s) durante as duas estaçoes do ano. As alturas de ondas foram também maiores no setor leste (máximo de 1,05 m durante a maré enchente na estação menos chuvosa) e os períodos de ondas foram mais curtos (<4,5 s) no setor oeste. A média granulométrica obtida dos sedimentos coletados na face praial apresentou escala mais freqüente entre 2,6 a 2,8 phi, indicando a predominância de areia fina. O grau de seleção predominante dos sedimentos foi de 0,2 a 0,5 phi (muito bem selecionados e bem selecionados), e da assimetria foi de positiva (0,10 a 0,30) e de aproximadamente simétrica (-0,10 a 0,10). O grau de curtose variou desde muito platicúrtica (<0,67) a muito leptocúrtica (1,50 a 3,00). Foram observados eventos de acresção sedimentar da estação chuvosa a menos chuvosa. De 22/07/1988 a 28/08/2013 (25 anos) também houve predomínio de acresção, onde o avanço médio linear da linha de costa foi de 190,26 m. O recuo médio linear obtido para toda área de estudo foi de -42,25 m. Áreas com maior erosão são pontuais: divisas das praias da Corvina e Maçarico, e Farol Velho e Atalaia. Os traps portáteis indicaram uma maior quantidade de sedimentos transportados longitudinalmente na estação menos chuvosa (Mín. 280 g/m³: enchente, setor oeste; Máx. 1098 g/m³: vazante, setor leste). Nos traps de espraiamento, o balanço entre a quantidade de sedimentos entrando e saindo nas praias foi menor no setor central (Mín. 80 g/m³: vazante, estação menos chuvosa; Máx. 690 g/m³: enchente, estação menos chuvosa). A circulação costeira sedimentar é proveniente, principalmente, do efeito das marés, com direção governada pela enchente e vazante dos rios que atravessam a costa. Os dados indicam o transporte longitudinal de sedimentos da ilha de Atalaia e rio Sampaio para o setor oeste e as margens das faixas praiais.

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.

Desenvolvimento metodológico e avaliação de contaminação por HPAs em sedimentos da baia de Guajará

Os Hidrocarbonetos Policíclicos Aromáticos (HPAs) são poluentes de efeito tóxico, prejudiciais ao meio ambiente e à saúde humana, fazem parte de um grupo de compostos poluentes orgânicos persistentes (POPs), que por suas características tem impactado o ambiente, sendo por esse motivo bastante estudados. Podem estar presentes nas formas particulada, dissolvida e/ou gasosa, estando presentes em diferentes ambientes; solo, sedimento, ar, água, material particulado na atmosfera, organismos e alimentos (Kennish, 2007). As fontes naturais de HPAs incluem atividades vulcânicas, queimadas naturais, exsudação de óleos, além de processos biogênicos. HPAs antrogênicos podem ocorrer pela combustão incompleta de óleos combustíveis (automotores e industriais), queima intencional de madeira e plantações, efluentes domésticos e/ou industriais, drenagens pluviais urbanas, derrames acidentais de óleos e derivados. Hidrofóbicos e lipofílicos, essas substâncias podem ser facilmente adsorvidas em sedimentos, sendo este compartimento um importante reservatório desses poluentes. Para avaliar a presença desses compostos no ambiente, utilizou-se nesse trabalho a Cromatografia Líquida de Alta Eficiência. As amostras utilizadas no desenvolvimento e otimização da metodologia foram coletadas na baia do Guajará (Belém – PA). O presente trabalho constituiu-se dessa forma em um desenvolvimento de um procedimento metodológico (com adaptações e otimizações) para quantificar 16 HPAs em 10 pontos na baia do Guajará, Belém – PA, em duas etapas de campo, totalizando 20 amostras analisadas. Na etapa de desenvolvimento do método analítico foram testados sistemas de eluição, polaridade do sistema e fluxo do eluente entre outros. Para validação do método foram avaliados os parâmetros fidelidade, linearidade, limite de detecção, limite de quantificação do método. Razões diagnósticas foram calculadas para identificação das fontes primárias do HPAs encontrados na baia. Foram identificadas, a partir de razões diagnósticas da ΣHPAs BMM/ΣHPAsAMM; Fen/Ant; Flt/Pir; Ant/Σ178; Flt/Σ202; B(a)P/Σ228 e Ind(123cd) pireno/Σ276 as fontes primárias dos 16 HPAs estudados no sedimento da baia. A somatória das concentrações dos HPAs leves na primeira etapa de campo, variou de 132,13 ng.g^-1 a 1704,14 ng.g^-1, a ΣHPAs dos pesados de 125,82 g.g^-1 a 1269,71 ng.g^-1 e ΣHPAs totais de 317,84 ng.g^-1 a 3117,06 ng.g^-1.. Na segunda etapa de campo, as concentrações dos HPAs leves variou de 76,12 ng.g^-1 a 1572,80 ng.g^-1 ; a ΣHPAs pesados variou entre 213,90 ng.g^-1 a 1423,03 ng.g^-1, e Σ HPAs totais teve concentrações de 290,02 ng.g^-1 a 2995,82 ng.g^-1. A partir dos resultados obtidos pode-se classificar a baia do Guajará como moderadamente impactada. A combustão constitui a fonte predominante de HPAs nos sedimentos da baia do Guajará, seguida da combustão de biomassa vegetal e aporte de petróleo e derivados. A maioria dos pontos estudados nesse trabalho, nas duas etapas de campo, apresentaram concentrações de HPAs individuais acima dos VGQS.