Geração de campo magnético direcional por combinação de dipolos magnéticos

A ferramenta de indução EM é utilizada na perfilagem de poço de prospecção com a finalidade de determinar a concentração de óleo na formação geológica. Formações com camadas que são relativamente de grande espessura já foram bem tratadas, entretanto zonas com camadas finas intercaladas com argila e areia com bom potencial de hidrocarbonetos precisam ser investigadas. Por isso, é necessário melhorar a resolução da sonda de indução para delimitar camadas finas com maior precisão. Teoricamente, se uma sonda de indução apresentar o feixe do campo magnético direcionado transversalmente ao eixo do poço, ela poderá ter um bom potencial para melhorar a resolução vertical das sondagens. Neste trabalho estudamos a geração do campo magnético direcional na forma de um feixe fino através de diferentes combinações de duas bobinas transmissoras (dipolos magnéticos) para determinar um arranjo que apresente uma melhor focalização na perfilagem de indução. Duas configurações com dois transmissores EM apresentam uma região em que só há contribuição da componente radial para a indução magnética: dois dipolos magnéticos paralelos transversais ao eixo do poço e dois dipolos magnéticos anticoaxiais ao eixo do poço. Para estas duas configurações foram analisados os possíveis espaçamentos entre os dois dipolos magnéticos transmissores, verificando-se que quanto mais próximos os dois dipolos magnéticos transmissores mais agrupadas ficam as linhas de fluxo magnético numa direção transversal ao eixo do poço. Conseqüentemente é gerado um feixe fino e direcional, concentrando energia em uma porção de volume estreita a partir dos transmissores, tendo um bom potencial para melhorar a resolução vertical da sonda de indução EM de poço.

Identificação automática das primeiras quebras em traços sísmicos por meio de uma rede neural direta

Apesar do avanço tecnológico ocorrido na prospecção sísmica, com a rotina dos levantamentos 2D e 3D, e o significativo aumento na quantidade de dados, a identificação dos tempos de chegada da onda sísmica direta (primeira quebra), que se propaga diretamente do ponto de tiro até a posição dos arranjos de geofones, permanece ainda dependente da avaliação visual do intérprete sísmico. O objetivo desta dissertação, insere-se no processamento sísmico com o intuito de buscar um método eficiente, tal que possibilite a simulação computacional do comportamento visual do intérprete sísmico, através da automação dos processos de tomada de decisão envolvidos na identificação das primeiras quebras em um traço sísmico. Visando, em última análise, preservar o conhecimento intuitivo do intérprete para os casos complexos, nos quais o seu conhecimento será, efetivamente, melhor aproveitado. Recentes descobertas na tecnologia neurocomputacional produziram técnicas que possibilitam a simulação dos aspectos qualitativos envolvidos nos processos visuais de identificação ou interpretação sísmica, com qualidade e aceitabilidade dos resultados. As redes neurais artificiais são uma implementação da tecnologia neurocomputacional e foram, inicialmente, desenvolvidas por neurobiologistas como modelos computacionais do sistema nervoso humano. Elas diferem das técnicas computacionais convencionais pela sua habilidade em adaptar-se ou aprender através de uma repetitiva exposição a exemplos, pela sua tolerância à falta de alguns dos componentes dos dados e pela sua robustez no tratamento com dados contaminados por ruído. O método aqui apresentado baseia-se na aplicação da técnica das redes neurais artificiais para a identificação das primeiras quebras nos traços sísmicos, a partir do estabelecimento de uma conveniente arquitetura para a rede neural artificial do tipo direta, treinada com o algoritmo da retro-propagação do erro. A rede neural artificial é entendida aqui como uma simulação computacional do processo intuitivo de tomada de decisão realizado pelo intérprete sísmico para a identificação das primeiras quebras nos traços sísmicos. A aplicabilidade, eficiência e limitações desta abordagem serão avaliadas em dados sintéticos obtidos a partir da teoria do raio.

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 em geofísica de poço: um estudo sobre ambiguidade

A ambiguidade na inversão de dados de geofísica de poço é estudada através da análise fatorial Q-modal. Este método é baseado na análise de um número finito de soluções aceitáveis, que são ordenadas, no espaço de soluções, segundo a direção de maior ambiguidade. A análise da variação dos parâmetros ao longo dessas soluções ordenadas permite caracterizar aqueles que são mais influentes na ambiguidade. Como a análise Q-modal é baseada na determinação de uma região de ambiguidade, obtida de modo empírico a partir de um número finito de soluções aceitáveis, é possível analisar a ambiguidade devida não só a erros nas observações, como também a pequenos erros no modelo interpretativo. Além disso, a análise pode ser aplicada mesmo quando os modelos interpretativos ou a relação entre os parâmetros não são lineares. A análise fatorial é feita utilizando-se dados sintéticos, e então comparada com a análise por decomposição em valores singulares, mostrando-se mais eficaz, uma vez que requer premissas menos restritivas, permitindo, desse modo, caracterizar a ambiguidade de modo mais realístico. A partir da determinação dos parâmetros com maior influência na ambiguidade do modelo é possível reparametrizá-lo, agrupando-os em um único parâmetro, redefinindo assim o modelo interpretativo. Apesar desta reparametrização incorrer na perda de resolução dos parâmetros agrupados, o novo modelo tem sua ambiguidade bastante reduzida.

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 geofísica ambiental e forense nos cemitérios do Benguí e do Tapanã (Belém - PA)

A Geofísica é uma ferramenta útil no estudo da contaminação de água provocada por cemitérios, assim como na investigação criminalística. A contaminação provocada por cemitérios é um tema ainda pouco estudado. Por outro lado, Geofísica forense é quase desconhecida no Brasil. Para estudar esses dois temas – cemitérios e alvos forenses – foram eleitos dois cemitérios públicos: o cemitério do Bengui, já desativado, e o cemitério do Tapanã. A investigação geofísica nos dois Cemitérios foi realizada por meio dos métodos eletromagnéticos ground penetrating radar (GPR) e Slingram de baixo número de indução (low induction number – LIN). No Cemitério do Bengui, testou-se o GPR na detecção do nível hidrostático, visto que, na implantação de um Cemitério, a informação mais importante é a profundidade do aquífero. No Cemitério do Tapanã, onde foi implantado o Campo de Testes controlados para Geofísica Forense, Ambiental e de Resgate, denominado FORAMB, utilizou-se GPR e Slingram em levantamentos mensais a fim de detectar o nível hidrostático, monitorar o desenvolvimento da contaminação e, especialmente, encontrar alvos forenses (um cadáver, uma caixa simulando um túnel e uma caixa com ferros simulando armas). O GPR mostrou-se eficaz em todos os estudos, especialmente na delineação dos alvos forenses. O Slingram não forneceu informações úteis.

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 aquíferos na cidade de Manaus (AM) utilizando perfilagem geofísica de poço e sondagem elétrica vertical

O sistema de águas subterrâneas vigentes na cidade de Manaus (AM) é pertencente ao Aqüífero Alter do Chão, o qual é caracterizado por fácies sedimentares arenosa, argilosa, areno-argilosa e “Arenito Manaus” depositadas em ambiente fluvial e flúvio-deltáico. Estima-se que 32.500 km³ de água possam ser explorados, fornecendo poços confinados a semiconfinados com vazão de até 300 m³/h em poços de aproximadamente 220 m de profundidade. Grandes distorções, contudo, têm sido notadas, Manaus, em foco neste trabalho, tem poços com baixo rendimento (de 78 m³/h), devido à má colocação de filtros em conseqüência da falta de informações sobre os aqüíferos. Através do uso integrado de perfilagem geofísica de poço (perfis de raio gama, de potencial espontâneo e de resistência elétrica), sondagem elétrica vertical e informações litológicas de amostragem de calha, mapearam-se camadas litológicas permo-porosas com elevado potencial aqüífero nas zonas norte e leste da cidade de Manaus. Foram identificadas duas zonas aqüíferas. A primeira zona é marcada pelos 50 m iniciais de profundidade. Nesta zona distinguem-se três tipos de litologias: argilosa, arenosa e areno-argilosa com as mais variadas dimensões. Estas duas últimas litologias apresentam elevado potencial aqüífero. Mesmo assim, nesta zona não há tendência a conter um volume de água considerável, devido à limitada continuidade lateral dos corpos arenosos e estar sujeita aos estorvos antrópicos superficiais. Portanto, não indicada para explotação de água em grandes quantidades. A segunda zona se estende a partir dos 50 m até aproximadamente 290 m de profundidade. É caracterizada por dois tipos de litologias: arenosa e areno-argilosa. Ao contrário da primeira, a esta profundidade há propensão do armazenamento e exploração de água. Apresenta litologias com elevado potencial aqüífero, cujos corpos arenosos são mais espessos e possuem uma maior continuidade lateral que os da zona 1; e valores da veza o de quase 300 m³/h. Serve, então, para abastecimento público por causa das boas vazões, garantia de água de boa qualidade protegida de ações antrópicas. A análise feita quanto à variação lateral, em mapas, de resistência transversal e condutância longitudinal, revelou que nas porções inferiores a 50 m de profundidade os setores NE e SW correspondem às zonas permeáveis. Portanto, possuem menos faixas argilosas, sendo considerados os setores mais transmissivos do sistema aqüífero. Os poços com profundidade maiores que 50 m recomenda-se sua locações nas faixas NNW, SW e extremo sul, em virtude dos resultados obtidos através da perfilagens geofísicas de poços, as quais delineiam camadas permo-porosas de elevado potencial aquífero para a área estudada.

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.

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 à 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.

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