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.

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.

Mineralogy and chemistry of the green stone artifacts (muiraquitãs) of the museums of the Brazilian State of Pará

Os muiraquitãs foram considerados de proveniência asiática ou, também, como peças esculpidas pelas lendárias mulheres guerreiras, as Amazonas. São peças, hoje, muito raras, encontradas em alguns acervos de Museus. Estudos mineralógicos e químicos de 23 peças do acervo dos Museus de Gemas e do Encontro em Belém, Brasil, mostraram que os muiraquitãs podem ser constituídos, tanto de quartzo, como de albita, ou microclínio, pirofilita, variscita, anortita e tremolita (equivalente ao jade nefrítico), minerais frequentes em formações geológicas do Brasil. No entanto quatro peças são constituídas de jadeíta, ou seja, em jade jadeítico, raro e desconhecido na Amazônia e Brasil. A constatação da presença desse mineral reacende a discussão em torno da origem mineralógica dos muiraquitãs encontrados na Amazônia. Essa origem, antes da atual descoberta, era defendida como amazônica, devido à ausência de jade jadeítico nas peças pesquisadas e pelo fato de jadeíta não ter sido encontrada no Brasil, mas na América Central e na Ásia.

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.

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 anomalias eletromagnéticas de um condutor tabular vertical sob camadas parcialmente condutiva em multifrequência e multiseparação através de modelamento analógico

A resposta eletromagnética (EM) de um corpo condutivo envolvido por uma zona parcialmente condutiva, torna-se bastante diferente daquela de um corpo condutivo em um meio altamente resistivo. As zonas parcialmente condutivas, como por exemplo, rocha encaixante, halo de sulfetos disseminados ou manto de intemperismo, que envolvem o corpo condutivo, afetam a resposta EM de diferentes maneiras, dependendo de suas características físicas e geométricas e, em particular, do sistema de prospecção EM utilizado. Neste trabalho em modelamento analógico, foi feita uma análise de anomalias EM provocadas por corpos condutivos tabulares verticais sob manto de intemperismo, em levantamentos terrestres para diferentes sistemas de bobinas - horizontal coplanar, vertical coplanar e vertical coaxial - em oito frequências na faixa de 250 Hz a 35 kHz e separações entre as bobinas de 0,15; 0,20 e 0,25m. O manto de intemperismo foi simulado por folhas de aço finas dispostas horizontalmente e o corpo condutor principal por folhas de alumínio finas colocadas verticalmente. As dimensões das folhas foram determinadas de acordo com as condições de modelamento para o plano e o semi-plano. Foram utilizados três corpos e três mantos com diferentes espessuras e condutividades, simulando, deste modo, diversas situações geológicas. Os resultados mostraram que cada sistema de bobinas é afetado diferentemente pela presença do manto de intemperismo. Para a análise dos resultados foi plotado um conjunto de diagramas considerando os valores pico-a-pico das anomalias em fase e em quadratura. Um outro conjunto de diagramas mostra as amplitudes máximas em fase, que ocorrem quando a componente em quadratura se anula em uma frequência relativamente baixa para um conjunto de corpo-manto, e as amplitudes máximas em quadratura, que ocorrem quando a resposta em fase atinge um mínimo próximo de zero, em frequências relativamente altas. Com isto foi possível conhecer a faixa de frequências para cada sistema de bobinas, onde a resposta EM se encontra o mínimo afetada pela presença do manto de intemperismo. A maior amplitude na resposta é obtida no sistema horizontal coplanar e a menor no sistema vertical coplanar. Um aumento na separação entre as bobinas é acompanhado por um deslocamento da anomalia para baixas frequências. A faixa de frequências, onde a presença do manto tem pouca influência na resposta do corpo condutivo, e maior para o sistema vertical coaxial e menor para o sistema horizontal coplanar. Esses resultados dão uma luz para o conhecimento da posição e da largura da banda de frequências utilizável, assim como as melhores separações entre transmissor-receptor, para auxiliar no planejamento de sistemas de prospecção EM, de modo que a resposta fique o mais livre possível de sinais indesejáveis, tais como os causados pela presença do manto de intemperismo.

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 geofísico regional sobre águas subterrâneas na Ilha do Marajó-Pará-Brasil

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

Estudo em laboratório de anomalias de potencial espontâneo

Este trabalho é uma aplicação em laboratório do método de Potencial Espontâneo de Geofísica de Campo. Pares de diferentes rochas e solos residuais foram colocados em contato em laboratório para determinar se poderiam produzir diferença de potencial, como tem sido ocasionalmente observado durante investigações de campo. As amostras de rochas utilizadas foram calcário, basalto e riolito da área da Caldeira de Los Humeros em México. Os solos residuais são da área MM1 da Serra dos Carajás, Pará, Brasil. As medidas foram efetuadas usando sistema em seco (umidade relativa ambiente) e sistema úmido (com água adicionada). As mudanças nos potenciais medidos para as rochas e solos no laboratório mostram diferença no valor médio desde 5 mv até 50 mv entre os diferentes lados em contato. Assim o resultado desta investigação é positiva e indica que mudanças nos valores de Potencial Espontâneo podem originar-se desde mudanças nos tipos superficiais de solo ou pelo contato entre diferentes tipos de rochas.

Estudo da resposta de diferentes arranjos de bobinas na perfilagem de indução de poço: modelamento analógico

A sonda de indução tem sido largamente utilizada na perfilagem de poço para detectar reservatórios através do contraste de condutividade existente entre as zonas contendo hidrocarbonetos e as zonas contendo água conata. A sonda de indução, além de auxiliar na localização de formações que contêm hidrocarbonetos, possibilita o cálculo de reservas para análise da viabilidade de exploração do reservatório. Porém, apesar do grande avanço da sonda de indução, a qual utiliza arranjo coaxial, existem várias restrições ao seu uso. Camadas delgadas, poços preenchidos com lama condutiva e invasão profunda são algumas das situações onde não é adequada a utilização da sonda de indução. Considerando as limitações da sonda de indução é que desenvolvemos em laboratório um estudo sobre as diferentes disposições de bobinas para a perfilagem de indução. O estudo consiste em comparar a resposta da sonda de arranjo de bobinas coaxiais, com a resposta da sonda de arranjo de bobinas coplanares e verificar se as respostas encontradas resolvem alguns dos problemas apresentados pela sonda de bobinas coaxiais. As respostas das sondas de diferentes disposições de bobinas foram obtidas através da montagem de um sistema de modelamento analógico simulando um poço. Para o sistema modelo foi utilizado um fator de escala linear igual a 20. Um tanque contendo soluções salinas de diferentes condutividades simulava algumas situações geológicas controladas. O sistema permitia uma variação na frequência de 50 kHz a 100 kHz. Os resultados obtidos foram satisfatórios, em virtude do arranjo de bobinas coplanares ter apresentado melhor resolução vertical que o arranjo de bobinas coaxiais.

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.

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

Interpretação de anomalias magnéticas usando integrais finitas dos momentos das componentes

O objetivo do trabalho consiste em desenvolver uma metodologia para determinar os parâmetros de um corpo causador de uma anomalia magnética. Essa metodologia baseia-se na utilização de expressões analíticas, deduzidas para as integrais finitas dos momentos da anomalia de intensidade total e das componentes H_x, H_y e H_z ao longo dos três eixos de um sistema de coordenadas cartesianas. Por meio do ajuste entre essas expressões analíticas e as integrais numéricas finitas dos momentos das componentes magnéticas obtidas a partir de um levantamento geofísico, pode-se computar, através de um processo iterativo, os parâmetros de magnetização, posição, profundidade e dimensões de um corpo anômalo. No caso em que são conhecidas as medidas de somente uma componente magnética, ainda é possível aplicar o método pois as demais componentes podem ser obtidas por um esquema de filtragem matemática. A metodologia foi testada com sucesso para modelos dipolo pontual e linha de dipolos os quais são muito utilizados em Geofísica para interpretar anomalias magnéticas produzidas por corpos geológicos que possuem uma ou duas dimensões horizontais muito menores do que a profundidade.

Inversão conjunta de dados gravimétricos magnéticos e resistivo

Apesar das grandes vantagens decorrentes da interpretação conjunta de dados geofísicos, a aplicação da inversão simultânea destes dados tem sido pouco estudada, principalmente a nível de simulação. Este trabalho foi desenvolvido com intuito de cobrir parte desta deficiência. Neste trabalho, foi utilizado o método dos mínimos quadrados para comparar os resultados obtidos a partir de: a) inversão dos grupos de dados gravimétricos, magnéticos e resistivos tomados separadamente; b) inversão simultânea destes grupos de dados combinados dois a dois; c) inversão simultânea dos três grupos de dados. O trabalho é desenvolvido a partir de dados teóricos onde são simuladas diversas situações geológicas. A comparação dos resultados é efetuada a partir das estimativas dos parâmetros obtidos por cada inversão, pelos desvios padrões de cada parâmetro (inversão gravimetria-magnetometria para o ruído Gaussiano) e pela redução da ambiguidade, manifestada pela dependência das estimativas em relação a aproximação inicial dos parâmetros. Na maioria dos casos estudados as inversões conjuntas dos dados combinados dois a dois apresentam resultados bem superiores àqueles obtidos usando-se apenas um dos grupos de dados isoladamente, seja nas estimativas dos parâmetros, seja na redução da ambiguidade. Por sua vez, a inversão conjunta dos três grupos de dados apresentam resultados semelhantes às inversões dos grupos de dados combinados dois a dois, contudo em alguns casos a inversão conjunta dos três grupos de dados é a única totalmente independente da aproximação inicial. Dados gravimétricos e magnéticos de duas anomalias reais foram invertidos, produzindo sempre curvas estimadas bem ajustadas aos valores observados.

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.