Determinação automática da porosidade e zoneamento de perfis através da rede neural artificial competitiva

Duas das mais importantes atividades da interpretação de perfis para avaliação de reservatórios de hidrocarbonetos são o zoneamento do perfil (log zonation) e o cálculo da porosidade efetiva das rochas atravessadas pelo poço. O zoneamento é a interpretação visual do perfil para identificação das camadas reservatório e, consequentemente, dos seus limites verticais, ou seja, é a separação formal do perfil em rochas reservatório e rochas selante. Todo procedimento de zoneamento é realizado de forma manual, valendo-se do conhecimento geológico-geofísico e da experiência do intérprete, na avaliação visual dos padrões (características da curva do perfil representativa de um evento geológico) correspondentes a cada tipo litológico específico. O cálculo da porosidade efetiva combina tanto uma atividade visual, na identificação dos pontos representativos de uma particular rocha reservatório no perfil, como a escolha adequada da equação petrofísica que relaciona as propriedades físicas mensuradas da rocha com sua porosidade. A partir do conhecimento da porosidade, será estabelecido o volume eventualmente ocupado por hidrocarboneto. Esta atividade, essencial para a qualificação de reservatórios, requer muito do conhecimento e da experiência do intérprete de perfil para a efetiva avaliação da porosidade efetiva, ou seja, a porosidade da rocha reservatório, isenta do efeito da argila sobre a medida das propriedades físicas da mesma. Uma forma eficiente de automatizar estes procedimentos e auxiliar o geofísico de poço nestas atividades, que particularmente demandam grande dispêndio de tempo, é apresentado nesta dissertação, na forma de um novo perfil, derivado dos perfis tradicionais de porosidade, que apresenta diretamente o zoneamento. Pode-se destacar neste novo perfil as profundidades do topo e da base das rochas reservatório e das rochas selante, escalonado na forma de porosidade efetiva, denominado perfil de porosidade efetiva zoneado. A obtenção do perfil de porosidade efetiva zoneado é baseado no projeto e execução de várias arquiteturas de rede neural artificial, do tipo direta, com treinamento não supervisionado e contendo uma camada de neurônios artificiais, do tipo competitivo. Estas arquiteturas são projetadas de modo a simular o comportamento do intérprete de perfil, quando da utilização do gráfico densidade-neutrônico, para as situações de aplicabilidade do modelo arenito-folhelho. A aplicabilidade e limitações desta metodologia são avaliadas diretamente sobre dados reais, oriundos da bacia do Lago Maracaibo (Venezuela).

Efeitos do manto de intemperismo sobre anomalias “VLF” dos corpos tabulares inclinados

A utilização dos métodos indutivos de propagação E.M. na exploração mineral em regiões tropicais, apresenta grandes dificuldades devido a presença de uma camada superficial condutiva (manto de intemperismo) comumente encontrada nestas regiões. Na região Amazônica, o manto apresenta-se bastante desenvolvido e condutivo, e em regiões semi-áridas, pode-se formar uma fina crosta superficial de sal. Em conseqüência disto, a interpretação dos dados E.M. obtidos para modelos que não consideram uma cobertura condutiva levam a erros consideráveis. Objetivando-se estudar os efeitos do manto sobre anomalias VLF devidas a corpos tabulares inclinados em contato com o manto (manto ohmico), foi realizada uma série de experimentos através do modelamento analógico, considerando-se diferentes parâmetros de resposta para o manto e o corpo. O manto de intemperismo foi simulado por soluções de cloreto de amônia (NH₄Cl) dispostas horizontalmente e o corpo condutor por chapas de grafite colocadas em posições inclinadas verticalmente. Utilizou-se quatro corpos condutores e três mantos com diferentes espessuras e condutividades, simulando, desta forma, diversas situações geológicas. Os resultados são dados por simples situações dos corpos localizados em um meio não condutor (ar), onde os parâmetros variados são: profundidade do topo, condutividade e mergulho do corpo. Os efeitos da condutividade da cobertura são amplamente ilustrados e avaliados. Para a análise dos resultados, foi plotado um conjunto de curvas considerando-se os valores pico-a-pico das anomalias de "tilt angle" e de elipsidade. Os resultados foram sintetizados em um outro conjunto de curvas reunidas em diagramas de Argand. Estando ou não o manto presente, observou-se, tanto para o tilt angle quanto para a elipsidade, o efeito do aumento da profundidade é o de reduzir a magnitude pico-a-pico e a forma do pico da anomalia, fazendo com que este afaste-se do ponto de "cross-over". Para um condutor de mesma espessura, o aumento da condutividade causa um ligeiro aumento nas anomalias de tilt angle, e uma atenuação nas anomalias de elipsidade. O efeito geral na variação do mergulho do condutor é o de causar uma assimetria nos perfis de tilt angle e de elipsidade. O aumento da condutância do manto de intemperismo causa um acréscimo nas anomalias de elipsidade e uma ligeira diminuição nas anomalias de tilt angle; porém, a partir de um certo valor de condutância do manto (mantos mais condutivos) tanto as anomalias de tilt angle quanto as anomalias de elipsidade começam a atenuar. Há rotação de fase no sentido anti-horário, sendo mais intensa para grandes valores de número de indução do corpo. Na presença do manto, o corpo parece estar a uma profundidade inferior à verdadeira e a ser menos condutivo.

Eletroresistividade aplicada ao estudo de água subterrânea no município de Salvaterra-Marajó-Pará

A partir de 1974, quando da implantação do Núcleo de Ciências Geofísicas e Geológicas da UFPa., começou a ser desenvolvido um estudo de água subterrânea na Ilha do Marajó, sob o titulo de Projeto Marajó, financiado pela FINEP, que visava definir e delinear aquíferos em algumas áreas piloto da ilha. Este trabalho, realizado em uma dessas áreas piloto, no município de Salvaterra, Marajó, busca determinar áreas favoráveis para exploração de água subterrânea, através do uso do método geofísico de eletroresistividade. A área, de 200 Km², foi prospectada através de 65 sondagens elétricas verticais e 6 perfilagens horizontais, tendo sido utilizado o dispositivo Schlumberger nas sondagens e o dispositivo Wenner nas perfilagens. Os dados foram inicialmente interpretados através do método do Ponto Auxiliar, quando se obteve uma visão geral do comportamento resistivo dos terrenos sondados, sendo depois reinterpretados usando técnica computacional, indicando existir na área camadas com resistividades que variam desde 5 até 43000 ohms-metro. Observou-se a existência de um pacote formado por uma sucessão de camadas que apresentaram resistividades com valores de 500 a 4000 ohms-metro, interpretado como constituído de sedimentos arenosos com água doce (aquífero), tendo uma espessura que chegou a atingir mais de 60 metros em algumas áreas. Na área pesquisada, foram delimitadas três sub-áreas cobrindo cerca de 40% da área total e, nestas, as espessuras do aquífero vão de 20 a 80 metros. As profundidades do topo desse aquífero nas três sub-áreas variam de 1,3 a 6,6 metros. A interpretação geofísica foi comprovada através de quatro furos de sondagem até 20 metros. Dois destes furos, realizados nas áreas interpretadas como bons aquíferos, mostraram terrenos mais arenosos, e os outros dois, alotados fora daquelas áreas, apresentaram uma parcela considerável de sedimentos argilosos com algum material orgânico. Informações obtidas quando do levantamento dos poços escavados na região e das análises de amostras das águas destes poços, mostram que, em termos gerais, não existem problemas quanto à qualidade das águas desses aquíferos mais superficiais, com exceção de dois poços escavados muito próximos da borda da Baía do Marajó, nas localidades de Joanes e Monsarás, que apresentam um aumento na salinidade quando no pico da estação de estio. Pelo que foi visto, o método utilizado neste trabalho é uma técnica eficiente para localizar bons aquíferos em ambientes geológicos semelhantes ao da área pesquisada.

Estudo da estrutura geo-elétrica da Região do Juruá, AM, pelo método magnetotelúrico

Utilizando-se dados magnetotelúricos (MT), foi obtida uma imagem geo-elétrica nítida da região do Juruá, Bacia do Solimões, na forma de seções geo-elétricas. Os dados de campo foram registrados ao longo de três linhas de 15 km, espaçadas de 3.5 km, recobrindo uma área de 100 km². O espaçamento entre as 35 estações é irregular, variando de 400 m a 3500 m. A faixa de freqüências utilizada cobriu de 0.001 Hz até 300 Hz, o que permitiu investigar de 100 m até 60 km de profundidade. Os dados apresentam-se afetados pelo efeito de distorção estática. Para corrigir este efeito foi utilizada a mediana da resistividade do primeiro condutor, correspondente à Formação Solimões. Foi utilizado o invariante do tensor MT para interpretar a estrutura geo-elétrica do Juruá. As seções geo-elétricas foram obtidas a partir do agrupamento dos dados resultantes da transformação de Bostick e da inversão 1D de Occam, para cada estação. Foi identificada uma seqüência de camadas condutivas e resistivas, correspondentes ao pacote sedimentar, uma zona de falhas e o topo do embasamento geo-elétrico, caracterizando a Bacia do Solimões. Abaixo do embasamento geo-elétrico foram também identificados uma zona condutora, seguida por uma camada de baixa condutividade, a profundidades iguais ou superiores a 20 km. Esta camada é interpretada como sendo de composição de gabro, estando associada a processos de acreção vertical, intimamente ligados à estabilização crustal e espessamento da litosfera. Os resultados apresentam uma boa concordância com os perfis de resistividade de poços e dados sísmicos de superfície.

Estudos dos efeitos do manto de intemperismo sobre anomalias de VLF geradas por corpos inclinados: modelagem numérica através de programa eGs

As referências apresentadas e utilizadas através de métodos indutivos de propagação eletromagnética E.M., são utilizadas na prospecção de minerais, entre outros, haja vista que em regiões tropicais existem camadas superficiais condutivas (manto intemperismo) que levam a influenciar, evidentemente, na interpretação do resultado real. Em nossa região Amazônica, o manto apresenta-se bastante desenvolvido e condutivo. O objetivo deste trabalho é estudar os efeitos dos mantos sobre anomalias VLF, devidos a corpos tabulares inclinados em contato ou não com o manto (manto galvânico e manto indutivo). Deste modo, foram realizadas umas séries de experimentos em modelagens numéricas através da eGs, considerando-se diferentes parâmetros de resposta para o manto e para o corpo. Os trabalhos de modelagem numérica, bem como seus resultados foram apresentados por diversas situações, como: variação da profundidade do topo do condutor, variação do mergulho do condutor, variação da condutância do manto, entre outros. Para interpretação dos resultados, foi gerado um conjunto de curvas considerando-se valores pico-a-pico das anomalias de Tilt Angle e de Elipsidade. Os resultados foram sintetizados em um outro conjunto de curvas reunidas em diagramas de Argand. A medida em que se aumenta a profundidade, o pico-a-pico e a forma do pico da anomalia do Tilt Angle e da Elipsidade diminuem na presença do manto, ou na ausência do mesmo, fazendo com que este se afaste do ponto de “Cross-Over”. Para um condutor da mesma espessura, o aumento da condutância causa um aumento nas anomalias de Tilt Angle e da Elipsidade. O efeito geral na variação do mergulho do condutor é o de causar uma assimetria nos perfis de Tilt Angle e Elipsidade. O aumento da condutância do manto de intemperismo causa uma diminuição nas anomalias de Tilt Angle e de Elipsidade. Há rotação de fase no sentido anti-horário, sendo mais intensa para grandes valores de condutância do corpo. Na presença do manto, o corpo parece estar a uma profundidade inferior à verdadeira e, tendendo, a ser menos condutivo.

Hidrogeologia da área piloto Ponta de Pedras - Ilha do Marajó

Constatou-se um aquífero pouco profundo, substancial para o desenvolvimento da região em estudo: sistema aquífero Ponta de Pedras. Esse aquífero é livre em alguns locais e semi-confinado em sua maior parte. Sua profundidade de topo máxima encontrada foi de 14,0 m, porém, em várias sondagens no ocidente da área e na sondagem em Igarapé Vilar, seu topo não foi atingido, podendo estar a mais de 16,0 m de profundidade. As profundidades de sua base e as espessuras desse aquífero são maiores que 34,0 m e 17,0 m, respectivamente, na sondagem de Mangabeira. Sua alimentação é feita principalmente por águas meteóricas. Taxas de infiltração entre 10⁶ a 10⁷ m³ por dia em 1 Km², foram estimadas para o mês de fevereiro de 1977. Suas porosidades efetivas, estimadas entre 25% e 37%, permitiu calcular um volume de água subterrânea próximo de 250x10⁶m³. O coeficiente de Darcy (K) médio é de aproximadamente 200 litros por dia por centímetro quadrado do sistema aquífero. As águas subterrâneas estudadas têm as seguintes características físico-químicas: pH sempre ácido entre 2,4 e 6,7; condutividade elétrica entre 13 a 2.000 micromhos por centímetro, sendo que as mais condutivas são as da região da bacia do Rio Tijucaquara; sílica com teor médio de 10,4 mg/l; ferro total com teor máximo de 4,0 mg/l; cálcio e magnésio com teores bastante baixos implicando em águas moles na maioria das vezes; manganês com teor máximo de 0,15 mg/l; nitrogênio e fosforo com concentrações bastante baixas. O sódio e potássio são os elementos químicos que visualizam com facilidade a variação sazonal do quimismo dessas éguas. Nos finais dos períodos chuvosos essas águas subterrâneas são menos salinizadas, por outro lado, nos períodos de pequenas precipitações até início da estação chuvosa, o excesso de sais impede a potabilidade de várias dessas águas. O uso doméstico dessas águas está limitado normalmente pelas seguintes características: pH ácido; ferro total acima de 0,3 mg/l com máximo de 4,0 mg/l em várias delas; isentas de fluoretos; excesso de manganês em algumas delas (0,15 mg/l); excesso de cloreto em um piezômetro. O uso dessas águas subterrâneas na agricultura pode estar limitado em alguns locais da área onde: pH menor que 5,0; condutividades altas durante os períodos de menor precipitação; altas porcentagens de sódio (83% a 97%).

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.

Inversão gravimétrica usando regularização entrópica

Apresentamos um novo método para inversão gravimétrica da distribuição espacial do contraste de densidade no plano horizontal, baseado na combinação da maximização da entropia de ordem zero com a minimização da entropia de ordem um. O topo e a base das fontes gravimétricas são presumidos ser planos e horizontais e o modelo interpretativo consiste de uma malha de prismas justapostos em ambas às direções horizontais, sendo os contrastes de densidade de cada prisma os parâmetros a serem estimados. A maximização da entropia de ordem zero é similar ao vínculo de suavidade global, enquanto a minimização da entropia de ordem um favorece descontinuidades na distribuição do contraste de densidade. Conseqüentemente a combinação judiciosa de ambas pode levar a soluções apresentando regiões com contrastes de densidade virtualmente constantes (no caso de corpos homogêneos), separadas por descontinuidades abruptas. O método foi aplicado a dados sintéticos simulando a presença de corpos intrusivos em sedimentos. A comparação dos resultados com aqueles obtidos através do método da suavidade global mostra que ambos os métodos localizam as fontes igualmente bem, mas o delineamento de seus contornos é efetuado com maior resolução pela regularização entrópica, mesmo no caso de fontes com 100 m de largura separadas entre si por uma distância de 50 m. No caso em que o topo da fonte causadora não é plano nem horizontal, tanto a regularização entrópica como a suavidade global produzem resultados semelhantes. A metodologia apresentada, bem como a suavidade global foram aplicadas a dois conjuntos de dados reais produzidos por intrusões em rochas metamórficas. O primeiro é proveniente da região de Matsitama, no nordeste de Botswana, centro sul da África. A aplicação das duas metodologias a estes dados produziu resultados similares, indicando que o topo das fontes não é plano nem horizontal. O segundo conjunto provém da região da Cornuália, Inglaterra e produziu uma distribuição estimada de contraste de densidade virtualmente constante para a regularização entrópica e oscilante para a suavidade global, indicando que a fonte gravimétrica apresenta topo aproximadamente plano e horizontal, o que é confirmado pela informação geológica disponível.

Mapeamento de uma interface separando dois meios homogêneos através da reconstituição gravimétrica compacta: aplicação a dados do LEPLAC-IV

O problema do mapeamento de uma interface separando dois meios homogêneos tem sido tradicionalmente resolvido através da continuação para baixo da anomalia gravimétrica observada. Este procedimento requer a utilização de um filtro passa-baixa ou de um parâmetro de amortecimento, que tendem a diminuir a resolução do relevo estimado da interface. Além disso, o uso da continuação para baixo implica supor (desnecessariamente) que a interface é uma superfície harmônica. Desta forma, a utilização da continuação para baixo é restrita a interfaces que apresentam um relevo suave relativamente à sua profundidade média. As restrições impostas pela continuação para baixo do campo gravimétrico são reduzidas na técnica de inversão proposta para este estudo, cujo objetivo é mapear a interface crosta-manto baseada na incorporação de informações a priori sobre o mínimo momento de inércia da fonte causadora, em relação a um eixo de concentração de massa pré-estabelecido e coincidente com o topo da fonte anômala. Além da introdução de informações sobre o mínimo momento, foi incluída também a possibilidade de haver variação lateral e em profundidade de densidade em setores previamente estabelecidos, o que permite estender a aplicação desta técnica a situações geológicas mais complexas que envolvam a reconstituição de corpos justapostos, que sejam de diferentes naturezas, tais como aqueles presentes em zonas de transição crustal. O método de inversão apresentado foi avaliado em testes usando dados sintéticos e os resultados obtidos mostraram a eficiência desta técnica em recuperar corpos alongados horizontalmente, que apresentem como topo ou base, interfaces arbitrárias a serem mapeadas. Foi também constatada a eficácia do método em reconstituir corpos que comportem variação lateral de densidade que pudessem retratar situações em que se deseja investigar a natureza do embasamento em regiões de transição crustal. Com base nos resultados destes testes, a técnica foi então empregada na inversão gravimétrica de cinco perfis do LEPLAC IV que cruzam a região oceânica adjacente à margem continental sudeste brasileira. Os resultados desta aplicação permitiram inferir a natureza dos corpos que constituem o arcabouço estrutural da região onde encontram-se inseridas a Dorsal de São Paulo e a Zona de Fratura de Porto Alegre, assim como mapear a Moho subjacente.

Mapeamento de magnetização aparente usando regularização entrópica

Apresentamos um novo método de mapeamento de magnetização aparente no plano horizontal que combina a minimização da entropia de primeira ordem com a maximização da entropia de ordem zero dos contrastes de magnetização estimados. O modelo interpretativo é uma malha de prismas verticais justapostos ao longo de ambas as direções horizontais. Presumimos que o topo e a base das fontes magnéticas são planos e horizontais e estimamos os contrastes de magnetização dos prismas. A minimização da entropia de primeira ordem favorece soluções de bordas abruptas e a maximização da entropia de ordem zero evita a tendência de a fonte estimada ser um único prisma. Desta forma, uma combinação judiciosa de ambos os vínculos pode levar a soluções caracterizadas por regiões de contraste de magnetização virtualmente constantes separadas por descontinuidades abruptas. Aplicamos este método a dados sintéticos produzidos por intrusões simuladas em sedimentos e que apresentam topo e base planos e horizontais. Comparando nossos resultados com aqueles obtidos pelo vínculo da suavidade, mostramos que ambos os métodos produzem uma boa e equivalente localização do centro das fontes. Todavia, a regularização entrópica delineia as bordas do corpo com maior detalhe. Ambos os vínculos (suavidade e regularização entrópica) foram aplicados a uma anomalia real sobre um escarnito magnético em Butte Valley, Nevada, Estados Unidos. A regularização entrópica produziu uma estimativa da distribuição de magnetização com bordas mais abruptas, menor volume e maiores valores de magnetização aparente comparados àqueles produzidos pelo vinculo da suavidade.

Ecologia da ictiofauna e análise ecossistêmica das áreas de influência direta da UHE Coaracy Nunes, Ferreira Gomes - AP

O reservatório da UHE Coaracy Nunes no rio Araguari esta localizado entre os municípios de Ferreira Gomes e Porto Grande no estado do Amapá-Brasil, distando 200 km do Oceano Atlântico. A usina Coaracy Nunes foi a primeira hidrelétrica a ser construída na Amazônia brasileira, tendo suas obras iniciadas em 1967. O rio Araguari e o principal rio do estado do Amapá e representa fonte de geração de renda através da pesca, atividades agropecuárias em sua várzea, navegação, mineração, geração de energia e lazer. O presente estudo teve por objetivo avaliar as alterações impostas pela construção do reservatório da UHE Coaracy Nunes, através das assembleias de peixes de quatro áreas de influencia direta desta usina. Para isso, no período de maio de 2009 a julho de 2010, foram realizadas coletas bimensais, de peixes, com redes de malhas padronizadas variando de 1,0 a 10,0 cm entre nos adjacentes e outras técnicas auxiliares. A partir destas coletas, no capitulo 1 foi verificado a composição, abundancia (CPUEn) e biomassa (CPUEp) relativas da ictiofauna, eficiência amostral (curva do coletor, curva de rarefação e Jacknife) e descritores ecológicos de comunidades (riqueza, diversidade, equitabilidade e dominância) das assembleias das quatro áreas. Foram efetuadas analises de variância (ANOVA: bifatorial), Kruskal-Wallis, teste-T e Mann-Whitney para verificar se havia diferenças significativas dos descritores entre as áreas e períodos sazonais. Estas análises foram corroboradas por analises multivariadas de agrupamento (cluster), ordenamento (MDS), Anosim e Simper. No capitulo 2, os estados ecológicos das quatro áreas foram verificados utilizando como indicadores: curvas espécie abundancia, curvas K-dominância e curvas ABC, assim como modelos espécie-abundancia serie geométrica, log serie, log normal e broken stick, e modelo de regressão linear de espectros de tamanho. No capitulo 3, a estrutura trófica foi estimada a partir da categorização das espécies de cada área em 5 guildas: piscívora, onívora, detritívora, carnívora e herbívora. A abundancia, biomassa e índices ecológicos destas guildas foram estimados e verificados suas variações espaço-temporais, por analises de variância (ANOVA: bifatorial e Kruskal-Wallis) e teste t. No capitulo 4, a dieta das espécies mais abundantes das assembleias de cada área foi verificada e suas variações espaço-temporais detectadas por analise de variância (ANOVA: bifatorial e Kruskal-Wallis). Também foram estimados a amplitude e sobreposição de nicho das espécies mais abundantes, assim como a existência de competição entre as espécies através de modelagem nula. No capitulo 5 foi realizada a avaliação ecossistêmica das quatro áreas através de modelos de fluxo de biomassa na rede trófica do ecossistema, usando como instrumento de modelação o software Ecopath. Essas análises tinham por objetivo descrever as variações dos atributos ecológicos que quantificam as propriedades de maturidade, estabilidade e resiliência ecossistêmica que pudessem refletir os estados ecológicos dessas áreas. O modelo incluiu compartimentos funcionais desde produtores primários ate predadores de topo. No geral, todas as análises indicaram sensíveis alterações na ictiofauna atribuídas a implantação da UHE Coaracy Nunes, que se refletem nos três níveis de organização: ecossistema, comunidade (assembleia) e guilda. Os resultados indicaram a captura de 1.977 peixes distribuídos em 2 classes, 9 ordens, 23 famílias, 73 gêneros e 108 espécies. As curvas de acumulação de espécies e curvas de rarefação individualizadas demonstraram que houve suficiência amostral nas áreas Reservatório e Lacustre. Os resultados mostraram que a área Jusante foi mais rica, diversa e equitativa em relação as demais e que a sazonalidade não influenciou na variação destes índices. A abundancia relativa (CPUEn) foi superior nas áreas Reservatório e Lacustre e a biomassa relativa (CPUEb) foi superior na área Jusante, não havendo diferenças sazonais para esses descritores em todas as áreas. As analises de agrupamento (cluster) e ordenamento (MDS) da ictiofauna permitiram identificar a formação de três assembleias distintas: Jusante, Montante e uma assembleia que compreende as áreas Reservatório e Lacustre, ratificando a similaridade dessas duas áreas. Os resultados das curvas whitakeplot, ABC e K-dominância, assim como o ajuste satisfatório do modelo broken stick e os padrões das curvas de espectro de tamanho para a assembleia da área a jusante indicam que esta área foi a mais equilibrada em termos ecológicos. Nas áreas Lacustre e Reservatório, os resultados tanto do ajuste ao modelo serie geométrica, quanto os resultados das curvas whitake-plot, ABC e K-dominância e o espectro de tamanho, assim como os resultados das curvas e ajustes ao modelo série e menor espectro de tamanho para a assembleia da área Reservatório, refletem que os peixes destas áreas, em sua maioria, são indivíduos pequenos com elevada dominância e baixa equitabilidade, caracterizando comunidades típicas de áreas impactadas. A estrutura trófica das assembleias de peixes das áreas represadas (Reservatório e Lacustre) foram formatadas em função do barramento do rio, que isolou e fragmentou o ambiente, determinando sua modificação física, impondo o estabelecimento de uma ictiofauna de espécies pré-adaptadas as condições ambientais de represamento, diferente, em parte, da estrutura da ictiofauna fluvial pre-barramento, destacando as piscívoras, onívoras e detritívoras, que foram as mais ricas e abundantes em função da disponibilidade, nas duas áreas, dos recursos alimentares de sua preferencia. Os resultados demonstraram que as dietas das assembleias de todas as áreas foram similares quanto ao predomínio do consumo de peixes e detritos, seguidos de alimento vegetal aloctone, revelando um padrão com poucos nichos amplos e uma concentração maior de espécies com nichos mais estreitos. Contudo, o padrão de baixa amplitude trófica foi evidenciado pelo predomínio da guilda piscívora, somada as guildas detritívora e herbívora. A sazonalidade pouco influenciou na alimentação da maioria das espécies em todas as áreas. Os padrões comparativos da dieta entre as áreas Montante e Jusante com as áreas Reservatório e Lacustre indicam que a maioria das espécies das áreas represadas pertenciam as guildas piscívora, onívora e detritívora antes do barramento do rio, que colonizaram estes ambientes, influenciadas, principalmente, pela abundancia dos recursos alimentares de suas preferencias e das condições físicas ambientais favoráveis. Interações competitivas foram evidenciadas pelos modelos nulos, sugerindo que a competição também foi um fator importante na estruturação das assembleias. Ecossistemicamente, os quatro modelos de fluxo de biomassa representam ecossistemas com elevada produção primaria oriunda da floresta riparia e algas filamentosas, que são utilizadas parcialmente. A cadeia trófica baseada em detrito apresentou ser mais importante que a que tem como base a produção primaria nas áreas Reservatório e Lacustre. A maioria dos fluxos ocorre nos compartimentos de níveis tróficos baixos. As propriedades ecossistêmicas da área Jusante indicam que este ambiente se encontra mais desenvolvido e maduro em relação aos outros, caracterizado por resiliência e entropia altas. As áreas represadas (Reservatório e Lacustre) apresentaram atributos ecossistêmicos que lhe conferiram características de menos resiliente e menos maduro que as áreas de rio. A área Montante apresentou um padrão intermediário de resiliência, estabilidade e maturidade. Esses resultados evidenciam que apos quarenta anos da construção da barragem do reservatório de Coaracy Nunes, a fragmentação do ambiente proporcionou alterações ecossistêmicas negativas, refletidas nas assembleias de peixes das áreas acima do barramento e na analise ecossistêmica, evidenciando que a área jusante apresenta características de ambiente em bom estado ecológico, com baixa alteração de origem antrópica e capaz de suportar distúrbios.

Food categories reconstruction and feeding consumption estimates for the Sciaenid Macrodon ancylodon (Bloch & Schneider), and the congeneric fishes Stellifer rastrifer (Jordan) and Stellifer naso (Jordan) (Pisces, Perciformes) in the Caeté Estuary, Northern Coast of Brazil

Um método quantitativo para se estimar o consumo alimentar e o aporte energético das diferentes categorias alimentares é apresentado através da reconstrução das presas ingeridas com base em estruturas corporais não digeríveis. Para tal, o presente estudo estabelece, através do exame dos conteúdos estomacais de 1.086 exemplares dissecados de Macrodon ancylodon (Bloch & Schneider, 1801), Stellifer rastrifer (Jordan, 1889) e Stellifer naso (Jordan, 1889), as equações das relações funcionais entre o peso das presas e estruturas corporais. Com as categorias reconstruídas foi possível quantificar o alimento ingerido pelos espécimes. Os resultados indicaram que existe uma marcada diferença, tanto na composição das categorias alimentares, bem como no aporte energético acompanhando o desenvolvimento ontogênico do predador.

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.

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.