Simulação de perfis nucleares de poço em formações complexas

A identificação e descrição dos caracteres litológicos de uma formação são indispensáveis à avaliação de formações complexas. Com este objetivo, tem sido sistematicamente usada a combinação de ferramentas nucleares em poços não-revestidos. Os perfis resultantes podem ser considerados como a interação entre duas fases distintas: • Fase de transporte da radiação desde a fonte até um ou mais detectores, através da formação. • Fase de detecção, que consiste na coleção da radiação, sua transformação em pulsos de corrente e, finalmente, na distribuição espectral destes pulsos. Visto que a presença do detector não afeta fortemente o resultado do transporte da radiação, cada fase pode ser simulada independentemente uma da outra, o que permite introduzir um novo tipo de modelamento que desacopla as duas fases. Neste trabalho, a resposta final é simulada combinando soluções numéricas do transporte com uma biblioteca de funções resposta do detector, para diferentes energias incidentes e para cada arranjo específico de fontes e detectores. O transporte da radiação é calculado através do algoritmo de elementos finitos (FEM), na forma de fluxo escalar 2½-D, proveniente da solução numérica da aproximação de difusão para multigrupos da equação de transporte de Boltzmann, no espaço de fase, dita aproximação P₁, onde a variável direção é expandida em termos dos polinômios ortogonais de Legendre. Isto determina a redução da dimensionalidade do problema, tornando-o mais compatível com o algoritmo FEM, onde o fluxo dependa exclusivamente da variável espacial e das propriedades físicas da formação. A função resposta do detector NaI(Tl) é obtida independentemente pelo método Monte Carlo (MC) em que a reconstrução da vida de uma partícula dentro do cristal cintilador é feita simulando, interação por interação, a posição, direção e energia das diferentes partículas, com a ajuda de números aleatórios aos quais estão associados leis de probabilidades adequadas. Os possíveis tipos de interação (Rayleigh, Efeito fotoelétrico, Compton e Produção de pares) são determinados similarmente. Completa-se a simulação quando as funções resposta do detector são convolvidas com o fluxo escalar, produzindo como resposta final, o espectro de altura de pulso do sistema modelado. Neste espectro serão selecionados conjuntos de canais denominados janelas de detecção. As taxas de contagens em cada janela apresentam dependências diferenciadas sobre a densidade eletrônica e a fitologia. Isto permite utilizar a combinação dessas janelas na determinação da densidade e do fator de absorção fotoelétrico das formações. De acordo com a metodologia desenvolvida, os perfis, tanto em modelos de camadas espessas quanto finas, puderam ser simulados. O desempenho do método foi testado em formações complexas, principalmente naquelas em que a presença de minerais de argila, feldspato e mica, produziram efeitos consideráveis capazes de perturbar a resposta final das ferramentas. Os resultados mostraram que as formações com densidade entre 1.8 e 4.0 g/cm³ e fatores de absorção fotoelétrico no intervalo de 1.5 a 5 barns/e^-, tiveram seus caracteres físicos e litológicos perfeitamente identificados. As concentrações de Potássio, Urânio e Tório, puderam ser obtidas com a introdução de um novo sistema de calibração, capaz de corrigir os efeitos devidos à influência de altas variâncias e de correlações negativas, observadas principalmente no cálculo das concentrações em massa de Urânio e Potássio. Na simulação da resposta da sonda CNL, utilizando o algoritmo de regressão polinomial de Tittle, foi verificado que, devido à resolução vertical limitada por ela apresentada, as camadas com espessuras inferiores ao espaçamento fonte - detector mais distante tiveram os valores de porosidade aparente medidos erroneamente. Isto deve-se ao fato do algoritmo de Tittle aplicar-se exclusivamente a camadas espessas. Em virtude desse erro, foi desenvolvido um método que leva em conta um fator de contribuição determinado pela área relativa de cada camada dentro da zona de máxima informação. Assim, a porosidade de cada ponto em subsuperfície pôde ser determinada convolvendo estes fatores com os índices de porosidade locais, porém supondo cada camada suficientemente espessa a fim de adequar-se ao algoritmo de Tittle. Por fim, as limitações adicionais impostas pela presença de minerais perturbadores, foram resolvidas supondo a formação como que composta por um mineral base totalmente saturada com água, sendo os componentes restantes considerados perturbações sobre este caso base. Estes resultados permitem calcular perfis sintéticos de poço, que poderão ser utilizados em esquemas de inversão com o objetivo de obter uma avaliação quantitativa mais detalhada de formações complexas.

Tomografia eletromagnética para caracterização de reservatórios de hidrocarbonetos

Na produção de petróleo é importante o monitoramento dos parâmetros do reservatório (permeabilidade, porosidade, saturação, pressão, etc) para o seu posterior gerenciamento. A variação dos parâmetros dinâmicos do reservatório induz variações na dinâmica do fluxo no reservatório, como por exemplo, perdas na pressão, dificultando o processo de extração do óleo. A injeção de fluidos aumenta a energia interna do reservatório e incrementa a pressão, estimulando o movimento do óleo em direção aos poços de extração. A tomografia eletromagnética poço-a-poço pode se tomar em uma técnica bastante eficaz no monitoramento dos processos de injeção, considerando-se o fato de ser altamente detectável a percolação de fluidos condutivos através das rochas. Esta tese apresenta o resultado de um algoritmo de tomografia eletromagnética bastante eficaz aplicado a dados sintéticos. O esquema de imageamento assume uma simetria cilíndrica em torno de uma fonte constituída por um dipolo magnético. Durante o processo de imageamento foram usados 21 transmissores e 21 receptores distribuídos em dois poços distanciados de 100 metros. O problema direto foi resolvido pelo método dos elementos finitos aplicado à equação de Helmhotz do campo elétrico secundário. O algoritmo resultante é válido para qualquer situação, não estando sujeito às restrições impostas aos algoritmos baseados nas aproximações de Born e Rytov. Por isso, pode ser aplicado eficientemente em qualquer situação, como em meios com contrastes de condutividade elétrica variando de 2 a 100, freqüências de 0.1 a 1000.0 kHz e heterogeneidades de qualquer dimensão. O problema inverso foi resolvido por intermédio do algoritmo de Marquardt estabilizado. A solução é obtida iterativamente. Os dados invertidos, com ruído Gaussiano aditivo, são as componentes em fase e em quadratura do campo magnético vertical. Sem o uso de vínculos o problema é totalmente instável, resultando em imagens completamente borradas. Duas categorias de vínculos foram usadas: vínculos relativos, do tipo suavidade, e vínculos absolutos. Os resultados obtidos mostram a eficiência desses dois tipos de vínculos através de imagens nítidas de alta resolução. Os tomogramas mostram que a resolução é melhor na direção vertical do que na horizontal e que é também função da freqüência. A posição e a atitude da heterogeneidade é bem recuperada. Ficou também demonstrado que a baixa resolução horizontal pode ser atenuada ou até mesmo eliminada por intermédio dos vínculos.

Análise de ambiguidade em métodos de campo potenciais através de análise fatorial, Q-Modal

Este trabalho apresenta uma metodologia para o estudo da ambiguidade na interpretação de dados geofísicos. Várias soluções alternativas, representativas da região de maior ambiguidade no espaço de parâmetros são obtidas, sendo posteriormente grupadas e ordenadas pela análise fatorial modo Q. Esta metodologia foi aplicada a dados sintéticos de campo potencial simulando-se causas de ambiguidade como discretização e truncamento da anomalia e a presença de ruídos aleatório e geológico. Um único prisma foi usado como modelo interpretativo, sendo a espessura a principal causa de ambiguidade tanto na gravimetria como na magnetometria. Segue-se a profundidade do topo sempre associada à espessura, quando o sinal da anomalia é alto. Quando a anomalia tem sinal baixo, a largura torna-se o segundo parâmetro mais importante, também associada à espessura. Ao contrário da presença de interferências geológicas, a presença de ruído aleatório nos campos, não é fator importante na ambiguidade. A aplicação da metodologia a dados reais ilustra o papel desta análise na caracterização de soluções alternativas e a importância da informação a priori na caracterização das causas de ambiguidade. A metodologia apresentada pode ser empregada em diversos estágios de um programa de prospecção fornecendo em cada estágio uma análise dos principais fatores causadores da ambiguidade, que poderá ser util no planejamento dos estágios seguintes. Comparada a outros métodos de análise de ambiguidade, como por exemplo regiões de confiança, a metodologia estudada destaca-se por não precisar satisfazer premissas estatísticas sobre a distribuição dos erros.

Aplicação dos métodos elétrico e eletricomagnético para o estudo hidrogeológico no município de Altamira-Pará

Este trabalho trata da aplicação dos métodos geofísicos eletroresistivo e eletromagnético no município de Altamira, região central do Estado do Pará, com a finalidade de fazer o mapeamento da espessura das camadas aqüíferas, bem como a ocorrência de intrusões de diabásio na área do município, permitindo uma melhor locação de pontos para perfuração com fins de captação à água subterrânea. Na área estudada, verificam-se afloramentos de rochas cristalinas do embasamento coberto por sedimentos com idade variando desde o Paleozóico ao Quaternário, com esperança de se obter água nos sedimentos do Devoniano Inferior, mais precisamente na formação Maecuru. Foram executadas 34 sondagens elétricas verticais com arranjo Schlumberger com AB/2 variando de 100 a 500 metros. 108 sondagens eletromagnéticas transientes também foram executadas com arranjo "loop coincidente", com 50 metros de lado. Por problemas eletrônicos adicionados aos ruídos de fenômenos supraparamagnéticos, os resultados das sondagens eletromagnéticas apresentaram-se muito distorcidos, não tendo sido possível extrair deles informação geológica de relevância. Por isso, foram criados modelos eletromagnéticos sintéticos representando a situação geológica da área. Os dados obtidos no campo em forma de curvas de resistividade aparente, foram plotados em papel bilogarítmo e interpretados com a ajuda de algoritmos direto e inverso os quais considera modelos de camadas horizontalmente estratificadas, homogêneas e isotrópicas. Os resultados da interpretação permitiram estabelecer modelos geoelétricos aproximados refletindo a situação geológica na subsuperfície. Os resultados da interpretação das sondagens elétricas, foram confrontados com os resultados da interpretação qualitativa do mapa de isogamas do campo magnético total, permitindo delimitar com maior precisão a zona de incidência de soleiras de diabásio, situada a norte da área. Dos quatro horizontes geoelétricos distinguidos, o terceiro é essencialmente constituído por arenitos, apresentando, por isso, melhores condições para armazenamento de água subterrânea.

Aplicação dos métodos geofísicos elétrico e eletromagnético na determinação de unidades sedimentares costeiras tropicais em Bragança, Nordeste do Pará

Esta dissertação mostra a aplicação de métodos geofísicos na determinação das unidades sedimentares da Planície Costeira Bragantina, que estão sendo estudados no Programa de Manejo e Dinâmica de Manguezais-MADAM (Mangrove Dynamics and Management). A finalidade deste trabalho é testar a metodologia geofísica no ambiente costeiro, visando auxiliar o mapeamento geológico e o entendimento da evolução estratigráfica holocênica da Planície Costeira Bragantina. Na realização deste trabalho, foram aplicados em diferentes unidades morfoestratigráficas da área de estudo: (a) o método convencional da eletrorresistividade através de sondagens elétricas verticais, usando o arranjo Schlumberger e (b) o método eletromagnético slingram (Horizontal Loop Eletromagnetic), com oito freqüências. Com as sondagens elétricas verticais determinaram-se os valores de resistividade em subsuperfície para os diferentes horizontes, correlacionando-os sempre que possível com as informações geológicas disponíveis. Em geral, os resultados obtidos com as sondagens elétricas não foram muito satisfatórios em termos das profundidades investigadas pelo fato do ambiente estudado ser muito condutivo, dificultando a penetração de corrente elétrica e contribuindo para que a profundidade de investigação fosse baixa. A aplicação do método eletromagnético mostrou-se mais simples e rápido do que as sondagens elétricas. A interpretação dos dados obtidos com esse método, realizada através da análise de perfis de medidas e de sondagens eletromagnéticas, permitiu inferir contatos laterais e estabelecer assinaturas geofísicas para os depósitos da planície costeira.

Aplicação dos métodos gravimétrico e magnetométrico para a definição do comportamento estrutural da faixa de dobramentos Araguaia

O presente trabalho consiste na interpretação de informações gravimétricas e aeromagnetométricas (Projeto Geofísico Brasil-Canadá - PGBC) da região setentrional da Faixa de Dobramentos Araguaia, envolvendo uma área de aproximadamente 129.000km², compreendida entre os meridianos 47°50’W e 50°30’W e paralelos 4°50'S e 9°00'S. Abrange porções sudeste do Estado do Pará, noroeste de Goiás e oeste do Maranhão. Os trabalhos de campo constaram de levantamentos gravimétrico e altimétrico, ao longo de rodovias que constituem a rede viária regional, e os resultados obtidos, após correções e reduções, foram então organizados de forma a constituírem um mapa de anomalias Bouguer, apresentado em escala 1:500.000. O padrão gravimétrico da faixa de dobramentos Araguaia se caracteriza por mostrar feições predominantemente longitudinais, com curvas isoanômalas de direções submeridianas, concordante com o comportamento litológico-estrutural conhecido para a área. Na porção centro-meridional desse mapa aparece uma zona de anomalia negativa, alongada e intensa (menos de -100mgal), a qual sofre duas importantes inflexões para NW, sendo uma correlacionável ao Lineamento Carajás, e a outra na altura do paralelo 5°30'S. Apresenta zonas de altos gravimétricos, notadamente na porção noroeste (domínio cratônico), relacionada à influência das metavulcânicas do Grupo Grão Pará, e na porção nordeste, devida a massas densas introduzidas na crosta. Destaca-se também o relativo alto gravimétrico acompanhando o flanco oriental do eixo de ocorrência das braquidobras, podendo ser devido a massas excedentes colocadas sob essa região, através de esforços tectônicos, produzindo elevações do nível de base da crosta e consequentemente do embasamento. De forma genérica o flanco oriental da faixa mostra um gradiente mais intenso que o ocidental, devido principalmente à configuração geométrica das estruturas nesse setor. Da mesma forma merece destaque o gradiente regional ascendente de sul para norte. A análise das cartas de intensidade magnética do PGBC fornecem condições à individualização de diversos domínios magnéticos, perfeitamente correlacionáveis litológico e estruturalmente com as unidades geotectônicas que compõem o arcabouço regional, assim como suas subunidades. A execução de modelamento simples, quer para as informações gravimétricas como magnetométricas, objetivaram fornecer subsídios semi-quantitativos à interpretação, de forma a auxiliar no estabelecimento do padrão geométrico do embasamento da faixa de dobramentos Araguaia, assim como das estruturas a ela impostas. Assim, por meio do processamento das informações gravimétricas, verificou-se um padrão geométrico na forma de fatias imbricadas, estabelecidas sobre o substrato. A grande anomalia negativa da porção centro-meridional da faixa pode estar relacionada a uma deficiência de massa profunda (depressão na base da crosta) ou ainda a uma associação entre contraste de densidade e espessamento do pacote de supracrustais, nessa porção. O baixo gravimétrico da porção sudeste da área, entre Guarai e Tupirama, é devido a uma depressão do embasamento, da mesma forma como a anomalia negativa associada ao Lineamento Carajás está relacionada a uma estrutura sinformal conhecida.

Avaliação de algoritmos para conversão de modelos de velocidade de tempo para profundidade

Ainda hoje, a migração em tempo é o processo de imageamento substancialmente empregado na indústria do petróleo. Tal popularidade é devida ao seu alto grau de eficiência e robustez, além de sua habilidade em focalizar refletores nos mais variados ambientes geológicos. Entretanto, em áreas de alta complexidade geológica a migração em tempo falha de tal forma que a migração em profundidade e um campo de velocidade em profundidade são indispensáveis. Esse campo é geralmente obtido através de processos tomográficos partindo de um campo de velocidade inicial. A conversão de campos de velocidade de tempo para profundidade é uma forma rápida de se obter um campo inicial mais consistente geologicamente para tais processos. Alguns algoritmos de conversão tempo-profundidade recentemente desenvolvidos baseados no traçamento de raios-imagem são revistos e um algoritmo alternativo baseado na propagação da frente de onda-imagem é proposto. Os algoritmos são aplicados a dados sintéticos bidimensionais e avaliados de acordo com suas eficiência e acurácia, destacando suas vantagens, desvantagens e limitações na obtenção de campos de velocidade em profundidade.

Construção de um modelo experimental simulando condições de perfilagem de indução

A ferramenta de indução tem sido, desde 1949, largamente utilizada na prospecção de petróleo pois auxilia na localização de formações que contém hidrocarbonetos. Ainda existem, porém, várias restrições quanto ao uso deste equipamento. Camadas excessivamente finas, poços preenchidos com lamas muito condutivas e grandes diâmetros de invasão são exemplos de situações onde seu desempenho não é bom. Para que se possa minimizar estas limitações é necessário o conhecimento da resposta da sonda diante de condições conhecidas. Isto pode ser feito via modelamento analógico e/ou através de modelamento numérico. O presente trabalho consiste na montagem de um sistema de modelamento analógico capaz de simular o perfil de duas bobinas coaxiais frente a meios infinitos e camadas espessas invadidas. O meio geológico é representado através de solução salina de diferentes condutividades. O fator de escala linear é igual a 20 e o equipamento executa medidas na faixa dos 50 kHz aos 100 kHz. Este sistema será utilizado em futuros estudos envolvendo disposições de bobinas diferentes das normalmente utilizadas e em estudos de outros problemas encontrados na perfilagem de indução.

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.

Os efeitos Delaware e Groningen: um estudo quantitativo por elementos finitos

Os efeitos Delaware e Groningen são dois tipos de anomalia que afetam ferramentas de eletrodos para perfilagem de resistividade. Ambos os efeitos ocorrem quando há uma camada muito resistiva, como anidrita ou halita, acima do(s) reservatório(s), produzindo um gradiente de resistividade muito similar ao produzido por um contato óleo-água. Os erros de interpretação produzidos têm ocasionado prejuízos consideráveis à indústria de petróleo. A PETROBRÁS, em particular, tem enfrentado problemas ocasionados pelo efeito Groningen sobre perfis obtidos em bacias paleozóicas da região norte do Brasil. Neste trabalho adaptamos, com avanços, uma metodologia desenvolvida por LOVELL (1990), baseada na equação de Helmholtz para H_Φ, para modelagem dos efeitos Delaware e Groningen. Solucionamos esta equação por elementos finitos triangulares e retangulares. O sistema linear gerado pelo método de elementos finitos é resolvido por gradiente bi-conjugado pré-condicionado, sendo este pré-condicionador obtido por decomposição LU (Low Up) da matriz de stiffness. As voltagens são calculadas por um algoritmo, mais preciso, recentemente desenvolvido. Os perfis são gerados por um novo algoritmo envolvendo uma sucessiva troca de resistividade de subdomínios. Este procedimento permite obter cada nova matriz de stiffness a partir da anterior pelo cálculo, muito mais rápido, da variação dessa matriz. Este método permite ainda, acelerar a solução iterativa pelo uso da solução na posição anterior da ferramenta. Finalmente geramos perfis sintéticos afetados por cada um dos efeitos para um modelo da ferramenta Dual Laterolog.

Equipamento eletromagnético para prospecção geofísica e modelamento reduzido

Nas últimas décadas, a prospecção por métodos eletromagnéticos vem-se constituindo numa técnica eficiente para prospecção mineral. O objetivo deste trabalho foi desenvolver um equipamento para prospecção eletromagnética quantitativa de corpos condutores, através do método dipolo-dipolo, podendo ainda ser usado em modelos reduzidos. Eletricamente, o sistema mede grandezas relacionadas ao acoplamento indutivo entre duas bobinas: transmissora e receptora. Elas são dispostas na superfície da terra, afastadas entre si, e a terra, desse modo, constitui o núcleo acoplador. Quando existem corpos condutores nas proximidades, estes são denunciados por alterações no comportamento do sinal induzido na bobina receptora. O equipamento compreende dois conjuntos: o transmissor e o receptor, além de acessórios. O transmissor gera um campo eletromagnético nas freqüências de 520 e 3.090 Hz, e um sinal de referência para o receptor, o qual é enviado através de um cabo. O receptor, inicialmente, separa o sinal induzido pelos campos secundários gerados por condutores, do campo normalmente recebido, quando a condutividade da subsuperfície é relativamente uniforme (campo primário). Em seguida, decompõe esse sinal em duas componentes ortogonais, uma em fase, e outra em quadratura com o campo primário. Através de duas escalas de precisão, as amplitudes dessas componentes são mostradas como percentagens do campo primário, com precisão de 1%. A sensibilidade do receptor é de 0,5 μV. O circuito eletrônico foi rigorosamente testado com preciso instrumental de laboratório. Em seguida, testou-se sua aplicação no Laboratório de Modelo Reduzido Eletromagnético do NCGG, refazendo-se experiências clássicas, encontradas na literatura especializada. No campo, foi experimentado próximo da cidade de Araci no Estado da Bahia, em áreas prospectadas pela "Rio Doce Geologia e Mineração S/A-DOCEGEO". Em ambos os casos, verificaram-se bons resultados.

Estudo de águas subterrâneas com métodos elétricos na região do semi-árido de Patos - Pb

O Nordeste do Brasil, em particular a região semi-árida do Estado da Paraíba, é muito castigada pelos períodos de grandes estiagens e escassez de recursos hídricos. A água de superfície não supre a demanda da região, portanto, a busca de água subterrânea e sua exploração é uma constante nos programas dos Governos Federal e Estadual. Nesta região, ocorrem geralmente apenas aquíferos aluviais rasos e aquíferos nas fraturas do embasamento cristalino. Tradicionalmente, as fraturas no cristalino são prospectados através de estudos aerofotogramétricos e geoestruturais, com média de 30% de furos secos. O dimensionamento dos aquíferos aluviais é feito através de perfurações sistemática do subsolo, que é uma metodologia muito demorosa e cara. O objetivo deste trabalho é mostrar que o emprego de métodos geofísicos, especialmente os métodos geoelétricos, são eficazes e econômicos para localizar e avaliar áreas promissoras de água subterrânea, tanto nos aquíferos aluviais como nos cristalinos. Neste trabalho, foram aplicados o convencional método da eletroresistividade e o método eletromagnético, ainda pouco usado na prospecção de água subterrânea, em três áreas selecionadas na região semi-árida de Patos - Pb, sendo duas representantes dos aquíferos aluviais e uma representante dos aquíferos cristalinos. Nos aquíferos aluviais, foi determinada a profundidade até o embasamento cristalino através de sondagens elétricas verticais (SEVs), e os respectivos limites laterais foram mapeados com perfilagens de eletroresistividade. Os resultados das SEVs são coerentes com as profundidades obtidas pelas sondagens mecânicas de confirmação mostrando, assim, a aplicabilidade desta metodologia na pesquisa deste tipo de aquífero. No aquífero cristalino, inicialmente determinou-se a orientação geral das zonas fraturadas do embasamento cristalino pela técnica modificada da sondagem elétrica vertical (SEV Radial) e, em seguida, a localização exata em planta das mesmas foram determinadas com perfilagens de eletroresistividade e eletromagnéticas no sistema horizontal coplanar nas frequências 880 Hz e 2640 Hz. Os resultados dos dois métodos utilizados foram satisfatórios, porém, a perfilagem eletromagnética é mais rápida, necessita de menos mão de obra, reduzindo assim os custos de prospecção de água subterrânea no embasamento cristalino.

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.

Geofísica aplicada à prospecção de água subterrânea na área do Rio Paracauari – Ilha de Marajó - Pará

A região do Rio Paracauari apresenta sérios problemas hídricos devido o excesso de chuvas no inverno e escassez de água no verão. A indústria agropecuária, principal atividade econômica da região, sofre tremendamente durante o período de secas. Foram realizados estudos geofísicos através do método de resistividade para localizar e mapear os aquíferos rasos de boa potencialidade. Foram feitas 53 sondagens elétricas verticais, do tipo Schlumberger, e 5 perfis de resistividade do tipo Wenner. Foi utilizada a técnica do ponto auxiliar e, depois, os diagramas coletivamente foram tratados e processados no computador, para obtenção de modelos aproximados. Depois da interpretação, foi possível separar duas áreas distintas. A primeira com resistividade alta até 950 Ω.m, com predominância de areia, é constituída de excelentes aquíferos e, provavelmente, trata-se de "Paleocanais". A segunda, com resistividade baixa até 0,5 Ω.m, com predominância de sedimentos argilosos, sílticos, com água salobra. Com os mapas de resistência transversal total e de resistividade foram definidas as áreas mais favoráveis à exploração de água potável. A área favorável é somente 30% da área total. Foi feita coleta das águas rasas e superficiais da área, e suas resistividades variam de 362 Ω.m a 1,1 Ω.m. Algumas perfurações foram feitas, e vieram constatar a existência destas áreas distintas. Com isto, foram solucionados parcialmente os problemas de água de algumas fazendas como São Lourenço, Gavinho e Conceição.

Métodos geofísicos aplicados ao mapeamento geológico numa área da Região Bragantina-Pará

Este trabalho trata da aplicação de métodos geofísicos, numa área piloto da região Bragantina, NE do Pará, com a finalidade de auxiliar no mapeamento geológico. A área acha-se coberta por Sedimentos Terciários espessos e apresenta ocorrências de rochas industriais como calcário e granito, representando uma superfície aflorante menor que 5% da área. Estudos geoquímicos, realizados por Goulart (1978), revelaram que no horizonte A, é impossível diferenciar os materiais de superfície quanto às suas origens. Isto torna difícil o mapeamento geológico de superfície, podendo levar a extrapolações exageradas dos contatos, o que ocorre no mapeamento pré-existente. Os métodos geofísicos aplicados foram: Eletroresistividade e Magnetometria. Também foram realizados alguns furos de sondagens mecânicas. Os resultados obtidos contribuíram para um melhor conhecimento das feições estruturais mapeadas, bem como tornou possível a descoberta de outras estruturas importantes. Definiram-se também duas áreas para subsequente estudo de detalhe, visando a separação do calcário, em melhor e pior qualidade para a fabricação do cimento. Quanto à região dos granitos, recomenda-se que seja feito o seu mapeamento sob forma de "ilhas de granito" em meio à espessas coberturas de solos e sedimentos.