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.

Interface gráfica para pré e pós-processamentos de dados elétricos e eletromagnéticos em Geofísica

O presente trabalho consiste no desenvolvimento de uma interface gráfica para pré e pós-processamentos de dados elétricos e eletromagnéticos em geofísica utilizando-a na transformação da EGSLIB – Electrical Geophysics Software Library, a qual é constituída por um acervo de programas FORTRAN para DOS, em um software integrado e que possua uma eficiente interface gráfica com o usuário. O resultado deste esforço foi a surgimento de um produto que foi denominado EGS 2000 – Electrical Geophysics Suite - 2000 version. Construído em ambiente DELPHI da Borland, através de um processo de desenvolvimento que, além de orientado a objetos, foi centrado em uma arquitetura que permitiu um alto grau de aproveitamento de programas FORTRAN já desenvolvidos e testados. Tal produto, constitui-se em uma poderosa ferramenta de auxílio à modelagem e interpretação de dados geofísicos totalmente executável no sistema operacional MS-Windows. Com isso, alcançou-se uma forma economicamente viável para dar uma destinação prática aos resultados teóricos obtidos nas universidades e institutos de pesquisa. Os exemplos apresentados, mostram-nos que, auxiliadas por computador, modelagens e interpretações de dados geofísicos alcançam um nível de produtividade e versatilidade jamais obtidos pelas ferramentas de outrora, uma vez que diversos modelos podem ser armazenados e recuperados e seus resultados comparados com dados observados sistematicamente colecionados, conferindo ao processo um alto grau de confiabilidade e rapidez.

Inversão quadrática de dados de AVO usando algoritmo de Marquardt

Neste trabalho são apresentadas expressões exatas e aproximações quadráticas das conhecidas equações de Knott-Zöeppritz, que calculam as amplitudes dos coeficientes Rpp e Rsp em termos dos contrastes médios relativos, bem como é procedida a inversão dos parâmetros de dados de AVO a partir destas aproximações quadráticas. Nesta inversão é utilizado o algoritmo de Levenberg-Marquardt, e são considerados apenas os eventos refletidos Rpp e convertidos Rsp, não associados. Nos estudos dos parâmetros físicos dos meios contrastes de impedância (δz), módulo de cisalhamento (δμ) e velocidade da onda p (δα), verificou-se quais desses parâmetros podem ser invertidos. Os resultados obtidos mostram que o contraste de impedância (δz) é muito bem resolvido estando ele relacionado com o contraste de velocidade da onda p (δα) ou com o contraste de cisalhamento (δμ), no caso de eventos refletidos considerando modelos de alto, moderado e baixo contrastes. Por outro lado ao se fixar o contraste de impedância e relacionar os outros dois parâmetros em consideração, os resultados mostram que esses são mal resolvidos, ou seja, a região de ambiguidade torna-se muito grande e os parâmetros tornam-se ambíguos e instáveis. No caso do evento convertido e na combinação do evento refletido com o convertido, para os modelos de baixo e moderado contrastes, (δz) é muito bem resolvido, caso que não acontece para modelo de alto contraste. Diante desses resultados verifica-se que no procedimento de inversão quadrática de dados de AVO, fixado (δμ), a recuperação dos dois parâmetros variados é muito boa, no caso do evento refletido, e razoavelmente boa no caso do evento convertido, por esse motivo optou-se pela fixação do módulo de cisalhamento.

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.

Investigação com os métodos GPR e magnético de urnas funerárias do Sítio Nossa Senhora do Perpétuo Socorro (Pará)

Este trabalho apresenta os resultados da aplicação de métodos geofísicos na eleição de áreas potenciais para salvamento arqueológico no Sítio Nossa Senhora do Perpétuo Socorro no estado do Pará. O Sítio possui urnas funerárias indígenas aflorantes e sub-aflorantes de bastante interesse arqueológico. Os métodos utilizados foram Magnetometria e Radar de Penetração no Solo (GPR). As medidas magnéticas foram utilizadas para indicar os locais mais apropriados para escavação reduzindo bastante as áreas potenciais. O GPR foi usado para confirmar as anomalias detectadas pelas medidas magnéticas e proporcionar uma melhor avaliação espacial tanto horizontal quanto verticalmente, diminuindo os erros cometidos quando se identificam anomalias magnéticas que não são causadas por feições arqueológicas. A metodologia de estudo baseou-se na comparação de anomalias obtidas em terrenos desconhecidos com terrenos onde comprovadamente havia a existência de elementos mapeáveis como urnas aflorantes e raízes, caracterizando e agrupando essas anomalias.

Mapeamento de jazidas de bauxitas no Distrito Paragominas – Açailândia (PA-MA) pelo método da sísmica de refração

Nos últimos 30 anos no Brasil, o mapeamento de camada bauxítica tem sido sempre efetuado através de perfilagens de poços. Neste estudo, aplicou-se o método da sísmica de refração rasa pela primeira vez, em cerca de 3.000m de perfis, em caráter experimental em algumas localidades ao longo da rodovia BR-10, assim como no Platô Esperança (região de Jabuti-PA) e perto da cidade de Açailândia (MA); com objetivo de delimitar a camada de bauxita. Embora se tenha utilizado, em algumas ocasiões, durante a interpretação, o método de distância crítica no cálculo de profundidade, foi possível, mediante o emprego do Método Recíproco Generalizado (GRM), o processamento e interpretação dos dados, para se obter a delimitação contínua da camada investigada. O GRM apresenta vantagens sobre os métodos tradicionais, tais como poder delimitar camadas onduladas, detetar camadas escondidas, além da sensibilidade para camadas com mergulho até 20°. A partir dós gráficos da análise de velocidade, foi calculada para a camada de bauxita uma velocidade média de 840 m/s, e mediante o emprego dos tempos de profundidade chegou-se a sua delimitação: ela possui um comportamento irregular, tendo uma espessura média de 5m. A pesquisa geofísica realizada nas três áreas-piloto tem grande importância porque o método utilizado poderá ser aplicado a outras regiões, em condições geológicas semelhantes.

O método GPR aplicado à localização de tubulações utilizadas no abastecimento de água na região urbana do município de Belém - Pará

Este trabalho consistiu em avaliar a viabilidade do método GPR (Ground Penetrating Radar) na localização de tubulações de aço, cimento-amianto, PVC e ferro fundido, de diversos diâmetros, utilizadas pela Companhia de Saneamento do Pará (COSANPA), no abastecimento de água à população da região urbana do município de Belém-Pará. Para o alcance deste objetivo foram realizados 14 perfis em 6 locais, distintos da área urbana de Belém, com o intuito de identificar 13 alvos, utilizando-se antenas de 200 MHz e 400 MHz, geralmente usadas neste tipo de levantamento. Os radargramas obtidos com o GPR permitiram identificar quase todos os alvos pesquisados. Apenas 3 dos 13 alvos, não foram plenamente identificados, ficando de fora, apenas os alvos de PVC, de pequena bitola (50 mm e 75 mm). Após o estudo concluímos que o método GPR é muito eficaz neste tipo de aplicação e que sua utilização tem grande valia na localização das redes de abastecimento de água. O que deve facilitar enormemente as obras de expansão da rede, assim como a sua manutenção.

Metodologia geofísica aplicada à pesquisa de salvamento arqueológico nos sítios Bittencourt, Jambuaçu e Jaburu, Estado do Pará

Este trabalho apresenta os resultados da aplicação de métodos geofísicos no salvamento arqueológico em três sítios no estado do Pará: Bittencourt, localizado em Abaetetuba; Jambuaçu, localizado em Moju; e Jaburu, localizado em Almeirim. Os métodos utilizados foram Magnetometria, Cintilometria e Radar de Penetração no Solo (GPR). As medidas magnéticas foram utilizadas para indicar os locais para escavação. O GPR foi usado para confirmar as anomalias detectadas pelas medidas magnéticas e proporcionou a diminuição de erros cometidos quando se identificam anomalias magnéticas que não são causadas por feições arqueológicas. As medidas de cintilometria foram usadas para tentar mapear os solos de Terra Preta Arqueológica (TPA) através da radiação gama natural do solo. Os locais indicados para escavação pela magnetometria e GPR possibilitaram a descoberta de materiais arqueológicos como cerâmica, louça, lítico e ferro. Os resultados obtidos com a cintilometria foram satisfatórios, mostrando haver relação entre as medidas mais baixas de cintilometria e as áreas com maior espessura de TPA. A metodologia geofísica utilizada demonstrou eficiência e rapidez sem causar impactos destrutivos nos sítios.

Métodos geofísicos aplicados à arqueologia no Estado do Pará

Este trabalho constitui uma primeira tentativa de aplicação de métodos geofísicos à pesquisa de sítios arqueológicos no Brasil. Esses métodos foram usados com a finalidade de detetar feições arqueológicas no subsolo, em dois sítios arqueológicos no estado do Pará: o sítio PA-SA-4, da fase arqueológica Areão e o sítio Ilha dos Bichos da fase Marajoara. No PA-SA-4, foi empregado o método magnético, enquanto que na Ilha dos Bichos, além do método magnético, foi usado o eletro-resistivo. O método eletro-resistivo possibilitou não só a determinação da camada de ocupação, pois esta camada se comporta geoeletricamente diferente das camadas que compõem as formações geológicas da região, como também permitiu detetar camadas compactas de argila, planas e horizontais, de origem arqueológica, situadas no lado oriental da Ilha dos Bichos. O método magnético levou à deteção de antigos locais de queima, devido a magnetização termo-remanente adquirida pelas argilas no processo de aquecimento/resfriamento, permitindo a descoberta de antigos fornos, bem como de uma urna funerária sob espessa camada de argila queimada.

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.

Métodos geofísicos aplicados ao ultramafito da Serra de Quatipuru - Estado do Pará

O objetivo deste trabalho foi a aplicação de métodos geofísicos para a avaliação econômica de ocorrência de níquel, cromo e amianto, em uma área de 6,25 km² - Malha 2 - Serra de Quatipuru, localizada a cerca de 50 km a SW de Conceição do Araguaia, no sul do Estado do Pará. Inicialmente foi aplicado o método magnético em caráter de reconhecimento, usando magnetômetros de campo total e componente vertical visando localizar zonas anômalas de amianto e cromita, uma vez que estes minerais se encontram associados com a magnetita. Em seguida foi usado o método eletromagnético de "dip angle" com o objetivo de definir zonas condutoras associadas a possíveis depósitos minerais ou características estruturais. Com base na interpretação destes dois métodos, definiu-se os locais de aplicação do método de polarização elétrica induzida-resistividade, que teve caráter de detalhe visando como alvo principal depósitos de sulfetos disseminados. As possíveis ocorrências minerais na área foram correlacionadas com anomalias geofísicas e confirmadas através de testemunhos de sondagem concluindo-se não serem de importância econômica. Contudo a aplicação destes métodos geofísicos foi útil, pois possibilitou o levantamento detalhado da área, que de outra forma só poderia ser realizado por meio dos métodos mais convencionais, tais como sondagem mecânica através de perfurações, os quais são muito caros e demorados.

Métodos sísmicos de refração e eletromagnético aplicados ao estudo de um paleocanal da Ilha do Marajó, Pará

Os métodos sísmico de refração e eletromagnético foram aplicados, em caráter experimental, à prospecção de água subterrânea. A eficiência destes métodos foi testada com cerca de 1.500m de perfis de refração rasa na delimitação geométrica em três dimensões de um paleocanal na região dos campos naturais da Ilha de Marajó e com um perfil eletromagnético (técnica Slingram) para a detecção de contatos laterais de água doce-água salobra. A interpretação dos dados sísmicos foi realizada com as equações de tempo de intersecção para n camadas inclinadas com o auxílio do controle geológico de 2 (dois) poços para contornar problemas de camadas não detectadas. Obteve-se uma baixa velocidade (0,81 km/s) para os sedimentos silte-arenosos saturados que preenchem o paleocanal. Em virtude da inversão de velocidades ter sido obtida como resultado para estes sedimentos, não se aconselha o mapeamento de paleocanais pela sísmica de refração, salvo com a utilização de técnicas especiais de campo (registro digital, 'up-hole survey', registro ondas S, reflexão simultânea, ...). Os dados eletromagnéticos foram analisados semi-quantitativamente com o auxílio de um modelo reduzido construído a partir de uma seção de resistividades aparentes. Um diagrama de Argand, traçado com os dados de laboratório, fornece subsídios para a determinação da condutância em ambientes de sedimentos não consolidados análogos ao de Marajó. Recomenda-se a utilização dos métodos eletromagnéticos para o reconhecimento lateral de zonas de água doce-água-salobra, uma vez que a componente quadratura da frequência alta (2640 Hz) define, com precisão suficiente, estas zonas. A comparação dos resultados dos dois métodos mostra que apenas parte do paleocanal se comporta como um aquífero.

Modelagem direta bidimensional do método magnetotelúrico com o método dos elementos finitos de arestas

Fizemos a modelagem direta 2D do método magnetotelúrico (MT) com o método dos elementos finitos (MEF) de arestas em termos dos campos primários e secundários. Para usarmos modelos de maior complexidade e diminuirmos o custo computacional utilizamos malhas não estruturadas. Nas malhas utilizadas, introduzimos quatro nós em torno de cada estação MT, constituindo um quadrado alinhado nas direções dos eixos cartesianos

Modelagem direta de dados de eletrorresistividade 3-D

Dentre os métodos geofísicos, o da Eletrorresistividade é um dos mais tradicionais, com o seu desenvolvimento ocorrido há mais de 80 anos. Durante esse tempo o seu uso acompanhou o avanço no poder de processamento numérico e mais recentemente, a modelagem e inversão tridimensional tornou-se uma possibilidade para o geofísico. Apresentamos, neste trabalho, a técnica de elementos finitos aplicada ao método da eletrorresistividade 3-D, através do cálculo do potencial secundário. Para o desenvolvimento da metodologia, simulamos o levantamento do método da eletrorresistividade 3-D com os arranjos Dipolo-Dipolo e Schlumberger, visando medir as variações laterais e verticais da resistividade aparente do solo. Estes arranjos consistem na injeção de corrente elétrica na superfície e de medidas de diferenças de potencial elétrico, resultante da interação da corrente elétrica com o solo. Sendo que, as fontes e receptores são localizados de acordo com os arranjos escolhidos para o levantamento. Neste trabalho, as curvas de sondagem e as pseudo-secções de resistividade aparente, são obtidas através da modelagem de eletrorresistividade 3-D, usando malha de elementos finitos regular. Para efeito de validação, os resultados são comparados com a resposta 3-D obtida a partir dos potenciais totais.

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

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