Avaliação geofísica e geológica de um porção de quebra de talude da Bacia do Jequitinhonha

A utilização dos métodos de reflexão sísmica na exploração e desenvolvimento de reservatórios de hidrocarbonetos ocorre devido à sua vasta e densa amostragem, tanto em área quanto em profundidade, aliada ao refinamento de técnicas para o tratamento dos dados de reflexão sísmica, a partir destes dados, são geradas seções sísmicas, que após a aplicação de tratamento adequado, são utilizadas na interpretação dos estratos e/ou estruturas geológicas da subsuperfície. Neste trabalho é feita uma análise Geofísica Geológica de duas linhas sísmicas reais 2D marinhas da porção de quebra de talude da Bacia do Jequitinhonha. Para tanto, foi realizado um conjunto de processamento sísmico com objetivo de atenuar as reflexões múltiplas comuns em dados marinhos, além disso, foram estimados os modelos de velocidade em profundidade, utilizados para determinação das seções sísmicas migradas em profundidade. Nestas foram identificadas as superfícies refletoras. Através da análise dessas superfícies foram feitas as marcações de sismofácies, com base nos conceitos iniciais da sismoestratigrafia, com a finalidade de avaliar a qualidade do produto derivado do processamento sísmico, empregado neste estudo, para uma interpretação sismoestratigráfica, a qual está fundamentada na análise dos padrões de terminações dos refletores e padrão interno das sismofácies.

Caracterização geoelétrica da área de curtume localizado no Distrito Industrial de Icoaraci, Belém-Pará

Foi desenvolvido um estudo da subsuperfície rasa para investigar a possibilidade de contaminação do subsolo por líquidos usados no tratamento de couros em área industrial. O estudo foi realizado na área da indústria Couro do Norte, Bairro de Maracacuera, no Distrito Industrial de Icoaraci (município de Belém-Pará). Durante o estudo foram feitas medidas para a determinação do fluxo de água do primeiro aqüífero e aplicada metodologia geofísica com a finalidade de detectar a possível contaminação do subsolo por poluentes liberados pelo curtume. As metodologias geofísicas utilizadas foram: eletrorresistividade (na forma de imageamento e través de sondagens verticais), potencial espontâneo, eletromagnético (Slingram) e radar de penetração no solo (GPR). A interpretação dos dados eletromagnéticos permitiu identificar a presença de zonas mais condutivas (provável pluma de contaminação) e zonas mais resistivas (blocos lateríticos). A distribuição dos potenciais medidos através do Método do Potencial Espontâneo indicou o sentido do fluxo subterrâneo local. Com as Sondagens Elétricas Verticais determinaram-se os valores de resistividade em subsuperfície para os diferentes horizontes correlacionando-os sempre com as informações obtidas no Caminhamento Elétrico. Através dos dados de GPR identificou-se a presença de blocos lateríticos, tubulações de metal e zonas com elevada atenuação (zona de sombra) que indicam a presença de contaminação.

Construção de banco de dados petrofísicos para determinação automática de litologias

Neste trabalho desenvolveu-se uma solução para o clássico problema da determinação automática de litologias. Construiu-se um banco de dados utilizando-se a técnica de análise discriminante aplicada sobre trechos testemunhados de poços perfilados de um campo situado em lâmina d'água profunda localizado na bacia de Campos. A técnica de análise de agrupamento, utilizada por um grande número de pesquisadores, foi testada para o mesmo conjunto de dados porém não forneceu resultados satisfatórios em comparação com os obtidos com a análise discriminante. Os dados de perfis, após sofrerem correções ambientais, são colocados em profundidade com os dados de testemunhagem e determina-se uma função discriminante que classifica as respostas dos perfis de acordo com as litofácies descritas nos testemunhos. Dessa forma obtém-se um banco de dados, a partir do qual se faz a discriminação de litologias em poços, testemunhados ou não, da mesma área. O banco de dados assim construído permite sua atualização através da incorporação de trechos testemunhados de novos poços. Mostra-se a necessidade de agrupar as litologias em quatro litofácies principais para este campo de maneira a padronizar os resultados, diminuir as indeterminações e melhorar a apresentação final. Como produto final obtém-se uma curva discriminada de eletrofácies que reproduz, com algum grau de certeza, as litofácies que reproduzem as litologias esperadas em subsuperfície. Como exemplo de aplicação deste método, utilizou-se as curvas discriminadas para correlacionar três poços, utilizando para isto um algoritmo de correlação automática.

Deconvolução de perfis de poço através de rede neural recorrente

Para a indústria do petróleo, a interpretação dos perfis de poço é a principal fonte de informação sobre a presença e quantificação de hidrocarbonetos em subsuperfície. Entretanto, em duas situações as novas tecnologias, tanto em termos do processo construtivo das ferramentas, quanto da transmissão dos dados não têm justificativa econômica, ensejando a utilização de um conjunto de perfis convencionais: reavaliações de campos maduros e avaliações de campos marginais. Os procedimentos de aquisição dos perfis convencionais podem alterar o valor da propriedade física bem como a localização dos limites verticais de uma camada rochosa. Este é um antigo problema na geofísica de poço – o paradoxo entre a resolução vertical e a profundidade de investigação de uma ferramenta de perfilagem. Hoje em dia, isto é contornado através da alta tecnologia na construção das novas ferramentas, entretanto, este problema ainda persiste no caso das ferramentas convencionais como, a ferramenta de raio gama natural (GR). Apresenta-se, neste trabalho, um novo método para atenuar as alterações induzidas no perfil pela ferramenta, através da integração do clássico modelo convolucional do perfil com as redes neurais recorrentes. Assume-se que um perfil de poço pode ser representado através da operação de convolução em profundidade entre a variação da propriedade física da rocha (perfil ideal) e uma função que representa a alteração produzida sobre a propriedade física, chamada como resposta vertical da ferramenta. Assim, desenvolve-se um processamento iterativo dos perfis, o qual atua na forma da operação de deconvolução, composto por três redes neurais recorrentes. A primeira visa estimar a resposta vertical da ferramenta; a segunda procura definir os limites verticais de cada camada rochosa e a última é construída para estimar o valor real da propriedade física. Este processamento é iniciado com uma estimativa externa tanto para o perfil ideal, quanto para a resposta vertical da ferramenta. Finalmente, mostram-se as melhorias na resolução vertical e na avaliação da propriedade física produzida por esta metodologia em perfis sintéticos e em perfis reais da formação Lagunillas, bacia do Lago Maracaibo, Venezuela.

Desenvolvimento de um conjunto de modelagem de alta resolução para os métodos eletromagnéticos indutivos

A Modelagem Eletromagnética Analógica tem sido uma poderosa ferramenta no estudo das respostas eletromagnéticas de uma grande variedade das técnicas de levantamentos e de modelos geológicos. Nas aplicações de métodos eletromagnéticos nas investigações de contaminação no subsolo, recentemente em grande demanda, são obtidas anomalias muito baixas (da ordem de 0,01% do campo primário). Com o objetivo de estudar em laboratório os problemas de levantamentos eletromagnéticos para alvos de baixo número de indução, um conjunto de modelagem em escala reduzida de alta sensibilidade, foi projetado e construído. Este conjunto opera na faixa de freqüência de 1 kHz até 1 MHz com excelente linearidade. O transdutor de transmissão gera suficiente densidade de fluxo magnético com corrente de até 3 App (pico a pico) circulando nele. O sistema de recepção, isto é, o conjunto de bobinas receptoras mais o pré-amplificador acoplado apresenta sensibilidade 3 V/A/m na freqüência de 10 kHz com inclinação de 20 dB/dec de freqüência. A alta sensibilidade do sistema de recepção combinada com a boa capacidade de geração do sistema transmissor, permite alcançar a resolução de 0,025% do campo primário. O conjunto dispõe ainda, de um sistema mecânico que permite posicionar o sistema de transdutores (transmissor e receptor) em posições arbitrárias com resolução de 1 mm e precisão de 0,05 mm. A precisão no deslocamento linear do sistema de sonda numa linha de 1715 mm de comprimento é alcançada por meio de um motor de passo controlado por um microcontrolador, que é rigidamente acoplado a um fuso que converte o movimento angular do motor em movimento linear. Com este conjunto, é possível simular uma grande variedade de técnicas de levantamentos eletromagnéticos no formato Slingram, isto é, dipolo-dipolo. Além da boa resolução mecânica e boa resolução de medida eletromagnética do conjunto, todo o processo de controle e aquisição de dados é operacionalizado por meio de um computador mestre.

Eletroresistividade aplicada à caracterização hidrogeológica para o abastecimento da cidade de São Domingos do Araguaia-Pará

O uso de métodos geofísicos (elétricos, eletromagnéticos e perfilagem de poço) na prospecção de água subterrânea em muitas localidades no Estado do Pará, tornou-se possível graças a um convênio firmado em 1988, entre a Fundação Nacional de Saúde e a Universidade Federal do Pará, através do Departamento de Geofísica e do Curso de Pós-Graduação em Geofísica, sendo estes responsáveis pelo estudo técnico das áreas prospectadas. Nosso objetivo nesta tese é contribuir com critérios geoelétricos, com base na utilização de dados de resistividade aparente (p_a), medidos em superfície, visando melhorar a caracterização do quadro geológico de subsuperfície e por conseguinte prover informações mais confiáveis quanto aos recursos dos mananciais subterrâneos. Inicialmente, analisamos a influência que algumas estruturas 2D, em subsuperfície, tem sobre os dados de eletroresistividade, medidos em Sondagens Elétricas Verticais (SEVs), com o arranjo Schlumberger. Este estudo foi realizado através de simulações numéricas utilizando o programa computacional SEV2D desenvolvido no Departamento de Geofísica, baseado na técnica dos elementos finitos. Resultados mostram ambiguidades geradas pelo uso de métodos de interpretação 1D de SEVs, executados em ambientes geológicos de características bidimensionais. A parte prática da tese se constituiu no tratamento interpretativo de dados de eletroresistividade coletados na sede do município de São Domingos do Araguaia, cidade localizada na região sudeste do Estado do Pará. A aquisição deste dados se deu através de SEVs, utilizando o arranjo Schlumberger. Correlacionando os resultados das interpretações geofísicas com informações geológicas foi possível definir um quadro geológico para a área que serve como referência para a prospecção de água subterrânea na referida cidade.

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.