Mineralogia e geoquímica de perfis de solo com Terra Preta Arqueológica de Bom Jesus do Tocantins, sudeste da Amazônia

A comparação de dados morfológicos, mineralógicos e químicos de solo com horizontes antrópicos - Terra Preta Arqueológica (TPA) com Argissolos adjacentes permitiu identificar os principais processos responsáveis pela formação da TPA em um sítio arqueológico no Município de Bom Jesus do Tocantins, sudeste do Estado do Pará. A similaridade entre os dados dos horizontes subsuperficiais do solo com TPA e solos adjacentes indica que o horizonte antrópico do solo TPA foi provavelmente desenvolvido a partir de um horizonte similar aos Argissolos adjacentes com posterior transformação pedogenética através da introdução de materiais orgânicos e inorgânicos por antigas colonizações humanas, resultando no espessamento do horizonte superficial e em concentrações maiores de CaO e P₂O₅ (teores totais), Zn (teor traço), P e Zn disponível (teores disponíveis), além de Ca e Mg trocáveis (teores trocáveis) em relação aos Argissolos adjacentes. Além disso, essa intervenção antrópica antiga também provocou modificações no horizonte subsuperficial do Argissolo com TPA, como concentrações altas de P₂O₅ e principalmente P disponível. O Soil Taxonomy e o Sistema Brasileiro de Classificação de Solos (SiBCS) são adequados para a identificação de solo com horizonte antrópico (p.exe. TPA), uma vez que priorizam nas ordens do solo os principais processos pedogenéticos atuantes na formação do solo, relacionados aos horizontes subsuperficiais, além das transformações pedogenéticas posteriores no horizonte superficial. Contudo, este trabalho recomenda o acréscimo de alguns atributos diagnósticos como quantidade de artefatos cerâmicos e líticos, P₂O₅, P e Zn disponíveis, C orgânico, Ca²++ Mg²⁺ (teores trocáveis), CTC e índice de saturação por bases no horizonte superficial para o agrupamento e distinção dos diversos tipos de solos antrópicos antigos da Amazônia.

Eletroresistividade na prospecção hidrogeológica no município de Ourilândia do Norte-PA

A prospecção geofísica para a exploração de água subterrânea em regiões cristalinas torna-se uma tarefa extremamente difícil pela presença de estruturas como falhas, fraturas e ocorrência de variação lateral de resistividade no manto intempérico. Para analisar a influência da variação bidimensional da propriedade física no comportamento das curvas de sondagem elétrica vertical Schlumberger, simulou-se modelos contendo tais estruturas encaixadas no manto de alteração com ângulo de mergulho de 30° e 90° e com a posição do centro de sondagem variando em relação a heterogeneidade. A modelagem numérica foi feita através do método dos elementos finitos com o programa EGSLIB/SEV2D. O resultado da simulação mostra que o emprego da técnica da sondagem elétrica vertical Schlumberger sobre meios bidimensionais não funciona eficientemente, como quando aplicada em regiões sedimentares e de aluvião, devido o fenômeno da ambigüidade que dificulta a interpretação das curvas. Dentre os métodos geofísicos para a investigação hidrogeológica, a aplicação dos métodos eletrorresistivos mostra-se mais simples, eficaz e econômica. Assim, o levantamento geofísico para orientar os trabalhos de perfurações de poços na cidade de Ourilândia do Norte contou de trinta e duas sondagens elétricas verticais tipo Schlumberger e um caminhamento elétrico dipolo-dipolo realizados em diversas ruas da cidade. As curvas de resistividade aparente, obtidas com as sondagens foram tratadas e processadas considerando um modelo de camadas planas, horizontais, homogêneas e isotropicas, inicialmente com o algoritmo EGSLIB/SEV1D e posteriormente com uso do programa EGSLIB/SEV1DINV, permitindo estabelecer modelos geoelétricos aproximados para a configuração da subsuperfície da cidade. Depois da interpretação quantitativa, foi possível separar duas zonas geoelétricas distintas. Esta diferenciação foi realizada com base nos valores de resistividade aparente e as informações geológicas de superfície. A primeira, com pouca representatividade na área da cidade, apresenta-se constituída por solo/aterro sobre saprolito de natureza argilo-arenosa recobrindo os granitóides tipo Rio Maria. A segunda composta pelo solo/aterro seguido de um horizonte preenchido por sedimentos arenosos superposto a camada argilo-arenosa e finalmente a rocha sã. Esta seqüência geoelétrica predomina na cidade. Os estratos arenoso e argilo-arenoso apresentam condições potenciais para armazenamento de água.

Estudo geofísico para a caracterização geológica e hidrogeológica da subsuperfície rasa em área de deposição de rejeitos sólidos no município de Barcarena-PA

Este estudo foi realizado na Área de Depósito de Rejeitos Sólidos (ADRS) da ALBRAS (Alumínio Brasileiro S.A.), localizada no município de Barcarena, Estado do Pará. Os resultados alcançados permitiram estabelecer um padrão de caracterização da subsuperfície, a partir de medidas geofísicas, fornecendo informações sobre as camadas geológicas e sobre parâmetros hidrogeológicos, tais como direção e velocidade do fluxo subterrâneo. Na realização do trabalho foram utilizados três métodos geofísicos: os métodos elétricos (Potencial Espontâneo e Eletrorresitividade) e eletromagnético (Slingram). A distribuição dos potenciais medidos através do Método do Potencial Espontâneo indicaram o sentido do fluxo subterrâneo local. Durante o estudo, foi montado um experimento de infiltração de solução salina no solo e seu efeito foi monitorado através de imageamento elétrico. Esse experimento permitiu que se estimasse a velocidade e o sentido do fluxo subterrâneo local ao longo da linha de imageamento. Os modelos interpretativos obtidos através da inversão de dados de resistividade aparente obtidos em sondagens elétricas verticais apresentaram uma boa correlação com perfis de raios gama corridos em poços tubulares da área de estudo. Estas Sondagens Elétricas Verticais permitiram que se detalhasse as unidades litológicas rasas através da estimativa de seus valores de resistividade e espessura. A interpretação dos dados eletromagnéticos através da análise de perfis de medidas, sondagens eletromagnéticas e mapas de contornos permitiu identificar a presença de zonas mais condutivas (material argiloso) e zonas mais resistivas (material arenoso). Além disso, observou-se que o mapa de contornos para a freqüência de 14080 Hz apresentou uma boa correlação com o mapa de contornos do Potencial Espontâneo referente a mesma área.

Estudo de pluma de contaminação por combustível com o método GPR e desenvolvimento de simulador geofísico para fins educativos

Este trabalho consistiu na investigação geofísica da possível contaminação de um posto de combustíveis em Belém-PA, bem como no desenvolvimento de um simulador geofísico para fins educativos. A investigação consistiu no levantamento de 13 perfis com o método geofísico eletromagnético Ground Penetrating Radar (GPR). Os radargramas obtidos, conjuntamente com os dados pré-existentes de sondagens mecânicas, Volatile Organics Compounds (VOC) e Benzeno-Tolueno-Etil-benzeno-Xileno (BTEX), permitiram identificar zonas de baixa reflexão do sinal eletromagnético (zonas atenuadas) em áreas com predominância de areia, o que indicaria contaminação. Como as medidas de GPR foram repetidas no período chuvoso, foi possível observar a redução da atenuação em perfil rico em areia, indicando a lixiviação dos ácidos graxos oriundos da biodegradação dos hidrocarbonetos, bem como o aumento da atenuação provocada pela argila, devido a sua elevada porosidade, e também, possivelmente, pelo aumento da biodegradação dos hidrocarbonetos que ficam retidos na argila. Este resultado constitui o primeiro caso inserido no simulador geofísico, experimento que permite mostrar o uso de diferentes métodos geofísicos. O simulador possui duas telas computacionais, a debaixo para simulação do modelo de subsuperfície escolhido e, a tela de cima, para representação da resposta geofísica obtida com o auxílio de uma unidade ótica que percorre o perfil representado na tela debaixo. Trata-se de uma modelagem analógica com recursos computacionais ainda não reportada na literatura, que permitirá mostrar diversos casos investigados constantes em dissertações e teses. O simulador poderá ser deslocado para empresas, universidades e secretarias, que fazem ou podem fazer uso da Geofísica, além de praças, escolas, etc., contribuindo de forma incisiva para a difusão da Geofísica como ferramenta na gestão de áreas contaminadas bem como para chamar a atenção da população para diferentes tipos de problemas, de modo a instruí-la sobre os mesmos.

Identificação de reflexões múltiplas, utilizando os parâmetros das frentes de onda PIN e normal

As reflexões múltiplas presentes nos sismogramas ocultam informações importantes sobre os refletores em subsuperfície e, podem até tornar completamente invisíveis as reflexões primárias, como no caso dos sismogramas marinhos, que muitas das vezes, exibem uma aparência anelar com fortes superposições das reflexões múltiplas, sobre as reflexões primárias. Problema este que tem sido alvo de importantes pesquisas, com o intuito de identificar, atenuar e/ou eliminá-las, através de vários métodos populares. O objetivo principal deste trabalho é a identificação das reflexões múltiplas. Com essa finalidade foi gerada, por modelamento direto, uma seção sísmica, com afastamento nulo (AN), contendo reflexões primárias e múltiplas simétricas de primeira ordem. Posteriormente, foi aplicada a migração cinemática do tipo Kirchhoff para obter o modelo em profundidade, apresentando uma boa recuperação dos refletores, bem como a presença de um refletor fictício, quando comparado com a seção anteriormente especificada. Foi obtida uma seção sísmica AN, do modelo migrado, na qual não é observado o segundo refletor, devido à ausência de contraste de impedância, entre a segunda e terceira camada, sendo este o primeiro indício de que o refletor fictício deste modelo é uma múltipla. Outro indício sobre a existência da múltipla foi a simetria encontrada entre as curvaturas do primeiro e terceiro refletor. Finalmente, foram calculados os parâmetros das frentes de ondas Hipotéticas Ponto de Incidência Normal (PIN) e Normal (N), bem como a velocidade Normal Moveout (NMO), tanto para os eventos de reflexões primárias como os eventos de reflexões múltiplas, para o modelo direto e para o modelo migrado. Em seguida, foram realizadas as comparações destes parâmetros, o que permitiu confirmar a veracidade dos indícios anteriores para a identificação das reflexões múltiplas.

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.

Investigação por eletroresistividade de ambientes cársticos no município de São João de Pirabas - PA

A presente tese procura investigar ambientes cársticos, usando os métodos eletroresistivos de sondagem elétrica vertical e caminhamento dipolo-dipolo. A área piloto para essa investigação está localizada na porção nordeste do Estado do Pará, no Município de São João de Pirabas. Nessa localidade, durante a perfuração de um poço para abastecimento de água para a cidade, a equipe de prospecção da Fundação Nacional de Saúde encontrou problemas devido a presença de cavidades nos calcários da Formação Pirabas. Para analisar esse problema foram realizadas 15 sondagens elétricas verticais, usando o arranjo Schlumberger, em diversos pontos da cidade. A abertura máxima dos eletrodos de corrente, para a realização das sondagens, foi em torno de 600 metros. As curvas de resistividade aparente das sondagens foram, inicialmente, interpretadas usando programas computacionais que calculam as resistividade do meio investigado com modelos unidimensionais, sem variações laterais de resistividade. Porém, as sondagens realizadas próximas ao poço, onde ocorrem cavernas, foram interpretadas usando modelos bidimensionais para representar as variações laterais de resistividade devido às cavernas. O algoritmo usado para essa interpretação se baseia na técnica dos elementos finitos. Os resultados das interpretações das sondagens elétricas, para os vários modelos de cavernas estudados, mostraram que as anomalias de resistividade devido a presença de cavernas são muito pequenas. Devido a este fato, o método de sondagem elétrica vertical não apresentou boa resolução para identificar com precisão a presença deste tipo de caverna em subsuperfície. Por outro lado, os resultados numéricos obtidos pelo caminhamento dipolo-dipolo mostraram que esse método é mais eficiente e menos ambíguo, para estudar cavernas, que o método de sondagem elétrica vertical. Contudo o método de caminhamento dipolo-dipolo é muito mais dispendioso do que o de sondagem elétrica verticais, tanto do ponto de vista de trabalho de campo como no processamento e interpretação de dados.

Métodos geofísicos aplicados a prospecção de água subterrânea no município de Ponta de Pedras - Ilha de Marajó - Pará

Este trabalho foi realizado numa área de 25 km² em Ponta de Pedras, Ilha de Marajó. Tem por objetivo o dimensionamento geométrico e seleção de zonas favoráveis à água subterrânea. Foram utilizados dois métodos geofísicos: sísmica de refração e eletroresistividade. O método sísmico foi empregado em caráter experimental, visando verificar a viabilidade de sua aplicação na área. Os resultados da interpretação confirmaram a suspeita inicial, da existência de horizontes geológicos pouco profundos, não possíveis de serem distinguidos pelo método sísmico. Com o método de eletroresistividade foram feitas sondagens verticais e perfis horizontais. Os resultados da interpretação indicaram áreas promissoras na parte central e numa faixa ao norte da área prospectada. A espessura do aquífero raso varia entre 0 e 50 m.

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.

Migração com amplitude verdadeira em meios bidimensionais (2-D) e introdução ao caso 2,5-D

Nos últimos anos tem-se verificado através de várias publicações um interesse crescente em métodos de migração com amplitude verdadeira, com o objetivo de obter mais informações sobre as propriedades de refletividade da subsuperfície da terra. A maior parte desses trabalhos tem tratado deste tema baseando-se na aproximação de Born, como em Bleistein (1987) e Bleistein et al. (1987), ou na aproximação do campo de ondas pela teoria do raio como Hubral et al. (1991), Schleicher et al. (1993) e Martins et al. (1997). Considerando configurações arbitrárias de fontes e receptores, as reflexões primárias compressionais podem ser imageadas em reflexões migradas no domínio do tempo ou profundidade de tal modo que as amplitudes do campo de ondas migrado são uma medida do coeficiente de reflexão dependente do ângulo de incidência. Para realizar esta tarefa, vários algoritmos têm sido propostos nos últimos anos baseados nas aproximações de Kirchhoff e Born. Essas duas abordagens utilizam um operador integral de empilhamento de difrações ponderado que é aplicado aos dados da seção sísmica de entrada. Como resultado obtém-se uma seção migrada onde, em cada ponto refletor, tem-se o pulso da fonte com amplitude proporcional ao coeficiente de reflexão naquele ponto. Baseando-se na aproximação de Kirchhoff e na aproximação da teoria do raio do campo de ondas, neste trabalho é obtida a função peso para modelos bidimensionais (2-D) e dois e meio dimensionais (2,5-D) que é aplicada a dados sintéticos com e sem ruído. O resultado mostra a precisão e estabilidade do método de migração em 2-D e 2,5-D como uma ferramenta para a obtenção de informações importantes da subsuperfície da terra, que é de grande interesse para a análise da variação da amplitude com o afastamento (ângulo). Em suma, este trabalho apresenta expressões para as funções peso 2-D e 2,5-D em função de parâmetros ao longo de cada ramo do raio. São mostrados exemplos da aplicação do algoritmo de migração em profundidade a dados sintéticos 2-D e 2,5-D obtidos por modelamento sísmico através da teoria do raio usando o pacote Seis88 (Cervený e Psencík, 1988) e os resultados confirmaram a remoção do espalhamento geométrico dos dados migrados mesmo na presença de ruído. Testes adicionais foram realizados para a análise do efeito de alongamento do pulso na migração em profundidade (Tygel et al., 1994) e a aplicação do empilhamento múltiplo (Tygel et al., 1993) para a estimativa de atributos dos pontos de reflexão - no caso o ângulo de reflexão e a posição do receptor.

Modelagem de eletrorresistividade 2-D a partir do potencial elétrico secundário

Um dos métodos clássicos da geofísica de exploração é o Método de Eletrorresistividade, estabelecido há um século pelos irmãos Schlumberger e desde então amplamente empregado em prospecção mineral, estudos ambientais e hidrogeologia e em pesquisa de fontes geotermais. Conceitualmente o método consiste de injeção de corrente elétrica na subsuperfície e de medida de diferença de potencial elétrico, resultante da interação da corrente com o meio. As localizações dessas fontes e receptores são determinadas pelo arranjo escolhido para o levantamento. Após o processamento, obtém-se pseudo-seções de resistividade aparente que indicam a distribuição de condutividade em subsuperfície. Devido à simplicidade dos fundamentos físicos de sua formulação, o método apresenta fácil implementação computacional quando comparado aos métodos eletromagnéticos de fonte controlada. Na literatura há inúmeros trabalhos de modelagem computacional, onde se calcula a resposta para problemas 2-D e 3-D. Nestes trabalhos, as pseudo-seções são obtidas a partir do cálculo do potencial elétrico total. Neste trabalho, apresentaremos a resposta da modelagem de eletroresistividade 2-D com o arranjo dipolo-dipolo, obtida a partir do potencial elétrico secundário. A solução é calculada através do método de elementos finitos usando malhas não estruturadas. Para efeito de validação, os resultados são comparados com a resposta 2-D obtida a partir dos potencias totais.

Modelagem do mCSEM no domínio do tempo usando transformada discreta de Fourier

A modelagem do mCSEM é feita normalmente no domínio da frequência, desde sua formulação teórica até a análise dos resultados, devido às simplificações nas equações de Maxwell, possibilitadas quando trabalhamos em um regime de baixa frequência. No entanto, a abordagem através do domínio do tempo pode em princípio fornecer informação equivalente sobre a geofísica da subsuperfície aos dados no domínio da frequência. Neste trabalho, modelamos o mCSEM no domínio da frequência em modelos unidimensionais, e usamos a transformada discreta de Fourier para obter os dados no domínio do tempo. Simulamos ambientes geológicos marinhos com e sem uma camada resistiva, que representa um reservatório de hidrocarbonetos. Verificamos que os dados no domínio do tempo apresentam diferenças quando calculados para os modelos com e sem hidrocarbonetos em praticamente todas as configurações de modelo. Calculamos os resultados considerando variações na profundidade do mar, na posição dos receptores e na resistividade da camada de hidrocarbonetos. Observamos a influência da airwave, presente mesmo em profundidades oceânicas com mais de 1000m, e apesar de não ser possível uma simples separação dessa influência nos dados, o domínio do tempo nos permitiu fazer uma análise de seus efeitos sobre o levantamento. Como parte da preparação para a modelagem em ambientes 2D e 3D, fazemos também um estudo sobre o ganho de desempenho pelo uso do paralelismo computacional em nossa tarefa.

Adubação nitrogenada em cultivos sucessivos de milho e soja sob irrigação

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Ocorrência de isoflavonas de soja no ambiente e correlação com atividades estrogênica : estudo de caso da região de Dourados (MS)

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

Implantação da cultura de cana-de-açúcar em argissolo vermelho do noroeste paulista: preparo do solo e uso de gesso

Pós-graduação em Agronomia - FEIS