Programa de Pós-Graduação em Geofísica - CPGF/IG

URI Permanente desta comunidadehttps://repositorio.ufpa.br/handle/2011/2355

O Programa de Pós-Graduação em Geofísica da UFPA (CPGF) do Instituto de Geociências (IG) da Universidade Federal do Pará (UFPA). Foi o segundo no Brasil a formar recursos humanos em Geofísica em nível de pós-graduação stricto sensu. Criado em 1972, funcionou até 1992 junto com os Cursos de Pós-Graduação em Geoquímica e Geologia.

Navegar

Agora exibindo 1 - 20 de 49

Acesso aberto (Open Access)
Análise da resolução do modelo de Hughes em sondagens magnetotelúricas
(Universidade Federal do Pará, 1996-02-13) SANTOS, Elinei Pinto dos; RIJO, Luiz; http://lattes.cnpq.br/3148365912720676
Na maioria dos métodos de exploração geofísica, a interpretação é feita assumindo-se um modelo da Terra uniformemente estratificado. Todos os métodos de inversão, inclusive o de dados eletromagnéticos, exigem técnica de modelamento teórico de modo a auxiliar a interpretação. Na literatura os dados são geralmente interpretados em termos de uma estrutura condutiva unidimensional; comumente a Terra é assumida ser horizontalmente uniforme de modo que a condutividade é função somente da profundidade. Neste trabalho uma técnica semi-analítica de modelagem desenvolvida por Hughes (1973) foi usada para modelar a resposta magnética de duas camadas na qual a interface separando as camadas pode ser representada por uma expansão em série de Fourier. A técnica envolve um método de perturbação para encontrar o efeito de um contorno senoidal com pequenas ondulações. Como a perturbação é de primeira ordem a solução obtida é linear, podemos então usar o princípio da superposição e combinar soluções para várias senoides de forma a obter a solução para qualquer dupla camada expandida em série de Fourier. Da comparação com a técnica de elementos finitos, as seguintes conclusões podem ser tiradas: • Para um modelo de dupla camada da Terra, as camadas separadas por uma interface cuja profundidade varia senoidalmente em uma direção, as respostas eletromagnética são muito mais fortes quando a espessura da primeira camada é da ordem do skin depth da onda eletromagnética no meio, e será tanto maior quanto maior for o contraste de condutividade entre as camadas; • Por outro lado, a resistividade aparente para este modelo não é afetada pela mudança na frequência espacial (v) do contorno; • Em caso do uso da solução geral para qualquer dupla camada na Terra cuja interface possa ser desenvolvida em série de Fourier, esta técnica produziu bons resultados quando comparado com a técnica de elementos finitos. A linerização restringe a aplicação da técnica para pequenas estruturas, apesar disso, uma grande quantidade de estruturas pode ser modelada de modo simples e com tempo computacional bastante rápido; • Quando a dimensão da primeira camada possui a mesma ordem de grandeza da estrutura, esta técnica não é recomendada, porque para algumas posições de sondagem, as curvas de resistividade aparente obtidas mostram um pequeno deslocamento quando comparadas com as curvas obtidas por elementos finitos.
Acesso aberto (Open Access)
Uma análise de métodos de supressão de múltiplas de superfície livre aplicados a um dado real
(Universidade Federal do Pará, 2011-10-14) OLIVEIRA, Andrei Gomes de; GOMES, Ellen de Nazaré Souza; http://lattes.cnpq.br/1134403454849660
O registro em levantamentos sísmicos 2D de eventos lineares são comuns nos traços dos sismogramas. Estes eventos se sobrepõem às reflexões de interesse (reflexões primárias) complicando as etapas do processamento sísmico e a identificação dos possíveis horizontes geológicos da subsuperfície. Pode-se destacar as múltiplas de superfície que se distinguem das demais por apresentarem baixas frequências, altas amplitudes e baixa velocidade de propagação. A identifificação e remoção de reflexões múltiplas em dados sísmicos reais é uma etapa de grande importância no processamento sísmicos, pois a atenuação destas múltiplas pode evitar erros na interpretação geológica das imagens sísmicas. O objetivo do trabalho é a determinação de um fluxo de processamento, para a atenuação de múltiplas de superfície livre. Métodos de eliminação de múltiplas de superfície (SRME, filltragem Radon, deconvolução preditiva no domínio τ - p e filtragem f-k) aplicados a um dado real 2D, foram utilizados em diferentes combinações. Estas são analisadas e são apresentados os resultados das melhores combinações.
Acesso aberto (Open Access)
Análise de velocidade por otimização do semblance na reflexão sísmica
(Universidade Federal do Pará, 2010) VIEIRA, Wildney Wallacy da Silva; LEITE, Lourenildo Williame Barbosa; http://lattes.cnpq.br/8588738536047617
Este trabalho teve como objetivo geral desenvolver uma metodologia sistemática para a inversão de dados de reflexão sísmica em arranjo ponto-médio-comum (PMC), partindo do caso 1D de variação vertical de velocidade e espessura que permite a obtenção de modelos de velocidades intervalares, vint,n, as espessuras intervalares, zn, e as velocidades média-quadrática, vRMS,n, em seções PMC individualizadas. Uma consequência disso é a transformação direta destes valores do tempo para profundidade. Como contribuição a análise de velocidade, foram desenvolvidos dois métodos para atacar o problema baseado na estimativa de velocidade intervalar. O primeiro método foi baseado na marcação manual em seções PMC, e inversão por ajuste de curvas no sentido dos quadrados-mínimos. O segundo método foi baseado na otimização da função semblance para se obter uma marcação automática. A metodologia combinou dois tipos de otimização: um Método Global (Método Price ou Simplex), e um Método Local (Gradiente de Segunda Ordem ou Conjugado), submetidos a informação à priori e vínculos. A marcação de eventos na seção tempo-distância faz parte dos processos de inversão, e os pontos marcados constituem os dados de entrada juntamente com as informações à priori do modelo a ser ajustado. A marcação deve, por princípio, evitar eventos que representem múltiplas, difrações e interseções, e numa seção pode ser feita mais de 50 marcações de eventos, enquanto que num mapa semblance não se consegue marcar mais de 10 eventos de reflexão. A aplicação deste trabalho é voltada a dados sísmicos de bacias sedimentares em ambientes marinhos para se obter uma distribuição de velocidades para a subsuperfície, onde o modelo plano-horizontal é aplicado em seções PMC individualizadas, e cuja solução pode ser usada como um modelo inicial em processos posteriores. Os dados reais da Bacia Marinha usados neste trabalho foram levantados pela PETROBRAS em 1985, e a linha sísmica selecionada foi a de número L5519 da Bacia do Camamu, e o PMC apresentado é a de número 237. A linha é composta de 1098 pontos de tiro, com arranjo unilateraldireito. O intervalo de amostragem é 4 ms. O espaçamento entre os geofones é 13,34 m com o primeiro geofone localizado a 300 m da fonte. O espaçamento entre as fontes é de 26,68 m. Como conclusão geral, o método de estimativa de velocidade intervalar apresentada neste trabalho fica como suporte alternativo ao processo de análise de velocidades, onde se faz necessário um controle sobre a sequência de inversão dos PMCs ao longo da linha sísmica para que a solução possa ser usada como modelo inicial ao imageamento, e posterior inversão tomográfica. Como etapas futuras, podemos propor trabalhos voltados direto e especificamente a análise de velocidade sísmica estendendo o caso 2D de otimização do semblance ao caso 3D, estender o presente estudo para o caso baseado na teoria do raio imagem com a finalidade de produzir um mapa continuo de velocidades para toda a seção sísmica de forma automática.
Acesso aberto (Open Access)
Anomalias eletromagnéticas de corpos tubulares inclinados em contato com o manto de intemperismo usando diferentes arranjos de bobinas-modelamento analógico
(Universidade Federal do Pará, 1997-03-16) SILVA, Roberto Cesar Alves; VERMA, Om Prakash; http://lattes.cnpq.br/2723609019309173
A principal aplicação dos métodos eletromagnéticos é a prospecção de sulfetos maciços que pela própria natureza geológica são corpos longos e delgados, encrustados em rocha encaixante resistiva, permitindo a interpretação através de curvas de modelos reduzidos usando semi-planos em meio resistivo (ar-livre). No entanto, com a extensão do uso dos métodos eletromagnéticos em regiões que tem um manto de intemperismo de parcialmente condutivo a condutivo, como na região Amazônica e nas regiões semi-áridas ou de climas tropicais , esta técnica de interpretação tem levado a resultados bem diferentes da situação real. O objetivo deste trabalho é estudar a influência do manto de intemperismo em contato ôhmico com o corpo-alvo, utilizando os arranjos de bobinas horizontal coplanar (HCP), vertical coplanar (VCP), vertical coaxial (VCA) e mínimo (PERP), através do modelagem analógico. Para simular o corpo foram utilizadas placas de grafite, cujas dimensões satisfazem a condição de semi-plano. Para simular o manto foi usada uma solução de cloreto de amônia. De forma geral, os resultados obtidos com os diversos arranjos apresentaram as mesmas características, ressalvada as diferenças peculiares de cada arranjo. De forma resumida tem-se as seguintes alterações na anomalia atribuídas à um manto condutivo em contato ôhmico com o corpo-alvo: i) Há rotação de fase, que inicialmente é em sentido anti-horário, mas mais adiante, com o aumento do número de indução do manto , torna-se horária; ii) Alteração na forma padrão do perfil de quadratura devido ao surgimento do pico extra no lado a favor do mergulho; iii) Só há uniformidade de comportamento para o pico-a-pico contra da quadratura, que cresce com o aumento número de indução do manto. Maiores anomalias são obtidas com o sistema horizontal coplanar e as menores com o vertical coplanar. Neste estudo não foi observado formação de pico extra nas anomalias do arranjo vertical coplanar. Ainda, os efeitos de cobertura são mais acentuados no sistema horizontal coplanar comparados aos sistemas de bobina vertical e em anomalias de semi-planos menos inclinados e em menor profundidade.
Acesso aberto (Open Access)
Aplicação dos métodos eletromagnéticos VLF e HLEM na prospecção hidrogeológica no município de São Domingos do Araguaia-PA
(Universidade Federal do Pará, 1996-04-25) FARIAS, Valcir João da Cunha; RIJO, Luiz; http://lattes.cnpq.br/3148365912720676
Inserido no convênio de cooperação firmado 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, este trabalho tem como objetivo o mapeamento geológico com a aplicação dos métodos eletromagnéticos VLF e HLEM, visando o estudo hidrogeológico do Município de São Domingos do Araguaia no Sudeste do Estado do Pará. Também, estudamos o desempenho quantitativo do método VLF, no que diz respeito a prospecção de água subterrânea. Primeiramente, analisamos a resposta VLF, através dos parâmetros da elipse de polarização (inclinação e elipsidade), de diferentes situações geológicas possíveis de serem encontradas no local de estudo e voltadas a hidrogeologia. Também, avaliamos a influência que os parâmetro físicos do modelo exercem na anomalia, utilizando para estas finalidades modelagem numérica por elemento finito. A análise quantitativa foi realizada através de ajuste, por tentativa e erro, com curvas teóricas obtidas da simulação de modelos idealizados na interpretação qualitativa, a qual fizemos com o auxilio da técnica de filtragem de Fraser e com perfís do método HLEM. Os resultados mostraram a existência de feições estruturais de "trends" N-S e NE-SW, que correlacionadas à geologia regional são falhas normais e de transferências, respectivamente. Com base nesses resultados, foram selecionados alguns locais com maiores probabilidades de sucesso na captação de água subterrânea através de poços.
Acesso aberto (Open Access)
Avaliação geofísica e geológica de um porção de quebra de talude da Bacia do Jequitinhonha
(Universidade Federal do Pará, 2013) SANTOS, Renata Vieira dos; SARGES, Roseane Ribeiro; http://lattes.cnpq.br/9737447758598999; GOMES, Ellen de Nazaré Souza; http://lattes.cnpq.br/1134403454849660
A utilização dos métodos de reflexão sísmica na exploração e desenvolvimento de reservatórios de hidrocarbonetos ocorre devido à sua vasta e densa amostragem, tanto em área quanto em profundidade, aliada ao refinamento de técnicas para o tratamento dos dados de reflexão sísmica, a partir destes dados, são geradas seções sísmicas, que após a aplicação de tratamento adequado, são utilizadas na interpretação dos estratos e/ou estruturas geológicas da subsuperfície. Neste trabalho é feita uma análise Geofísica Geológica de duas linhas sísmicas reais 2D marinhas da porção de quebra de talude da Bacia do Jequitinhonha. Para tanto, foi realizado um conjunto de processamento sísmico com objetivo de atenuar as reflexões múltiplas comuns em dados marinhos, além disso, foram estimados os modelos de velocidade em profundidade, utilizados para determinação das seções sísmicas migradas em profundidade. Nestas foram identificadas as superfícies refletoras. Através da análise dessas superfícies foram feitas as marcações de sismofácies, com base nos conceitos iniciais da sismoestratigrafia, com a finalidade de avaliar a qualidade do produto derivado do processamento sísmico, empregado neste estudo, para uma interpretação sismoestratigráfica, a qual está fundamentada na análise dos padrões de terminações dos refletores e padrão interno das sismofácies.
Acesso aberto (Open Access)
Comparação dos efeitos do eletrojato equatorial nos dados magnetotelúricos bi e tridimensionais
(Universidade Federal do Pará, 1995-11-01) PONTES, Altem Nascimento; RIJO, Luiz; http://lattes.cnpq.br/3148365912720676
Este trabalho trata da comparação dos efeitos do eletrojato equatorial nos dados magnetotelúricos bi e tridimensionais. Ele compõe-se fundamentalmente de duas partes: Na primeira, que consiste na comparação dos efeitos dos modelos bi e tridimensionais (linha de corrente e eletrojato gaussiano), determinamos as dimensões que deve apresentar uma estrutura tridimensional para que a resposta magnetotelúrica proveniente dessa estrutura, nas polarizações YX e XY, possa ser substituída pela resposta magnetotelúrica proveniente de uma estrutura bidimensional nos modos TE e TM respectivamente, para a resistividade aparente e fase. Os resultados mostram, para um embasamento condutivo ou resistivo, que é necessário aumentar a dimensão da estrutura tridimensional na direção principal acima de dezesseis vezes, em relação à dimensão inicial, para que a resposta magnetotelúrica proveniente dessa estrutura na polarização YX, possa ser substituída pela resposta magnetotelúrica proveniente de uma estrutura bidimensional no modo TE; no caso da polarização ser XY, essa substituição ocorrerá, no modo TM bidimensional, se aumentarmos a dimensão da estrutura tridimensional na direção principal acima de vinte e seis vezes em relação à dimensão inicial. Na segunda parte, que consiste na determinação dos efeitos do eletrojato equatorial Onwumechilliano, posicionamos a estrutura tridimensional a 10 km, 100 km, 300 km, 500 km, 1000 km e 1500 km, respectivamente, de distância do eletrojato. Dessa forma, determinamos a distância na qual os efeitos do eletrojato equatorial Onwumechilliano sobre as sondagens magnetotelúricas são praticamente desprezíveis. Verificamos que, posicionando a estrutura tridimensional a 1500 km de distância da fonte, para embasamento condutivo, o eletrojato não afeta as sondagens magnetotelúricas. Quando o embasamento é resistivo, o eletrojato afeta as sondagens magnetotelúricas, mesmo posicionando-se a estrutura tridimensional a 1500 km de distância da fonte. De um modo geral, constatamos que o eletrojato equatorial afeta as sondagens magnetotelúricas para freqüências inferiores a 10-2 Hz, tanto na resistividade aparente quanto na fase e em ambas as polarizações, principalmente, quando o embasamento é resistivo. Observamos também que o efeito galvânico é mais pronunciado na polarização YX que na XY.
Acesso aberto (Open Access)
Construção de banco de dados petrofísicos para determinação automática de litologias
(Universidade Federal do Pará, 1994-04-14) COUTO, Nuno André Pomar; LUTHI, Stefan Moritz
Neste trabalho desenvolveu-se uma solução para o clássico problema da determinação automática de litologias. Construiu-se um banco de dados utilizando-se a técnica de análise discriminante aplicada sobre trechos testemunhados de poços perfilados de um campo situado em lâmina d'água profunda localizado na bacia de Campos. A técnica de análise de agrupamento, utilizada por um grande número de pesquisadores, foi testada para o mesmo conjunto de dados porém não forneceu resultados satisfatórios em comparação com os obtidos com a análise discriminante. Os dados de perfis, após sofrerem correções ambientais, são colocados em profundidade com os dados de testemunhagem e determina-se uma função discriminante que classifica as respostas dos perfis de acordo com as litofácies descritas nos testemunhos. Dessa forma obtém-se um banco de dados, a partir do qual se faz a discriminação de litologias em poços, testemunhados ou não, da mesma área. O banco de dados assim construído permite sua atualização através da incorporação de trechos testemunhados de novos poços. Mostra-se a necessidade de agrupar as litologias em quatro litofácies principais para este campo de maneira a padronizar os resultados, diminuir as indeterminações e melhorar a apresentação final. Como produto final obtém-se uma curva discriminada de eletrofácies que reproduz, com algum grau de certeza, as litofácies que reproduzem as litologias esperadas em subsuperfície. Como exemplo de aplicação deste método, utilizou-se as curvas discriminadas para correlacionar três poços, utilizando para isto um algoritmo de correlação automática.
Acesso aberto (Open Access)
Decomposição em valores singulares aplicada a dados de campo magnético
(Universidade Federal do Pará, 1992-12-15) MOURA, Helyelson Paredes; O'BRIEN, Douglas Patrick
Neste trabalho, a decomposição em valores singulares (DVS) de uma matriz A, n x m, que representa a anomalia magnética, é vista como um método de filtragem bidimensional de coerência que separa informações correlacionáveis e não correlacionáveis contidas na matriz de dados magnéticos A. O filtro DVS é definido através da expansão da matriz A em autoimagens e valores singulares. Cada autoimagem é dada pelo produto escalar dos vetores de base, autovetores, associados aos problemas de autovalor e autovetor das matrizes de covariância ATA e AAT. Este método de filtragem se baseia no fato de que as autoimagens associadas a grandes valores singulares concentram a maior parte da informação correlacionável presente nos dados, enquanto que a parte não correlacionada, presumidamente constituída de ruídos causados por fontes magnéticas externas, ruídos introduzidos pelo processo de medida, estão concentrados nas autoimagens restantes. Utilizamos este método em diferentes exemplos de dados magnéticos sintéticos. Posteriormente, o método foi aplicado a dados do aerolevantamento feito pela PETROBRÁS no Projeto Carauari-Norte (Bacia do Solimões), para analisarmos a potencialidade deste na identificação, eliminação ou atenuação de ruídos e como um possível método de realçar feições particulares da anomalia geradas por fontes profundas e rasas. Este trabalho apresenta também a possibilidade de introduzir um deslocamento estático ou dinâmico nos perfis magnéticos, com a finalidade de aumentar a correlação (coerência) entre eles, permitindo assim concentrar o máximo possível do sinal correlacionável nas poucas primeiras autoimagens. Outro aspecto muito importante desta expansão da matriz de dados em autoimagens e valores singulares foi o de mostrar, sob o ponto de vista computacional, que a armazenagem dos dados contidos na matriz, que exige uma quantidade n x m de endereços de memória, pode ser diminuída consideravelmente utilizando p autoimagens. Assim o número de endereços de memória cai para p x (n + m + 1), sem alterar a anomalia, na reprodução praticamente perfeita. Dessa forma, concluímos que uma escolha apropriada do número e dos índices das autoimagens usadas na decomposição mostra potencialidade do método no processamento de dados magnéticos.
Acesso aberto (Open Access)
Deconvolução de perfis de poços através do ajuste de energia
(Universidade Federal do Pará, 1994-04-14) GUERRA, Carlos Eduardo; LOVELL, Michael Anthony
As medidas de resistividade são de fundamental importância para o cálculo da saturação de óleo em reservatórios potencialmente produtores. A combinação das medidas de resistividade rasa e profunda permite a obtenção dos parâmetros Rt, Rxo e di. Mas, em reservatórios complexos existem dificuldades em se obter leituras confiáveis de Rt, devido à baixa resolução vertical das ferramentas de investigação profunda. Em reservatórios laminados, por exemplo, as leituras obtidas pela ferramenta de indução profunda (ILD) podem levar a uma interpretação errônea das mesmas, levando a acreditar que as medidas obtidas do perfil referem-se a uma única camada. Este problema pode ser em parte resolvido através de uma metodologia que melhore a resolução vertical dos perfis de investigação profunda, valendo-se do uso de informações obtidas de um perfil de alta resolução vertical, i.e; a curva de resistividade rasa. Uma abordagem neste sentido seria usar um perfil de alta resolução que apresente uma boa correlação com o perfil de investigação profunda. Esta correlação pode ser melhor avaliada se aplicarmos um filtro no perfil de alta resolução, de tal forma que o perfil resultante tenha teoricamente a mesma resolução vertical do perfil de baixa resolução. A obtenção deste filtro, porém, recai na premissa de que as funções respostas verticais para as ferramentas de alta e baixa resolução são disponíveis, o que não ocorre na prática. Este trabalho se propõe mostrar uma nova abordagem onde o filtro pode ser obtido a partir de um tratamento no domínio da freqüência. Tal tratamento visa igualar a energia espectral do perfil de alta resolução à energia espectral do perfil de baixa resolução tendo como base o Teorema de Parseval. Será mostrado que a resolução vertical depende fundamentalmente da energia espectral do perfil em questão. A seguir, uma regressão linear é aplicada sobre os perfis de alta resolução filtrado e de baixa resolução. Para cada ponto amostrado dos perfis, uma rotina de minimização é aplicada visando escolher o melhor intervalo de correlação entre os perfis. Finalmente, um fator de correção é aplicado sobre cada ponto do perfil de baixa resolução. Os resultados obtidos com os perfis de indução são promissores, demonstrando a eficácia da abordagem e mesmo quando aplicada em perfis com diferentes propriedades petrofísicas, a metodologia funcionou satisfatoriamente, sem danificar os perfis originais.
Acesso aberto (Open Access)
Detecção do contato entre camadas através do perfil de indução
(Universidade Federal do Pará, 1996-08-21) MARQUES JÚNIOR, José; HOWARD JUNIOR, Allen Quentin; http://lattes.cnpq.br/6447166738854045
Dois dos principais objetivos da interpretação petrofísica de perfis são a determinação dos limites entre as camadas geológicas e o contato entre fluidos. Para isto, o perfil de indução possui algumas importantes propriedades: É sensível ao tipo de fluido e a distribuição do mesmo no espaço poroso; e o seu registro pode ser modelado com precisão satisfatória como sendo uma convolução entre a condutividade da formação e a função resposta da ferramenta. A primeira propriedade assegura uma boa caracterização dos reservatórios e, ao mesmo tempo, evidencia os contatos entre fluidos, o que permite um zoneamento básico do perfil de poço. A segunda propriedade decorre da relação quasi-linear entre o perfil de indução e a condutividade da formação, o que torna possível o uso da teoria dos sistemas lineares e, particularmente, o desenho de filtros digitais adaptados à deconvolução do sinal original. A idéia neste trabalho é produzir um algoritmo capaz de identificar os contatos entre as camadas atravessadas pelo poço, a partir da condutividade aparente lida pelo perfil de indução. Para simplificar o problema, o modelo de formação assume uma distribuição plano-paralela de camadas homogêneas. Este modelo corresponde a um perfil retangular para condutividade da formação. Usando o perfil de entrada digitalizado, os pontos de inflexão são obtidos numericamente a partir dos extremos da primeira derivada. Isto gera uma primeira aproximação do perfil real da formação. Este perfil estimado é então convolvido com a função resposta da ferramenta gerando um perfil de condutividade aparente. Uma função custo de mínimos quadrados condicionada é definida em termos da diferença entre a condutividade aparente medida e a estimada. A minimização da função custo fornece a condutividade das camadas. O problema de otimização para encontrar o melhor perfil retangular para os dados de indução é linear nas amplitudes (condutividades das camadas), mas uma estimativa não linear para os contatos entre as camadas. Neste caso as amplitudes são estimadas de forma linear pelos mínimos quadrados mantendo-se fixos os contatos. Em um segundo passo mantem-se fixas as amplitudes e são calculadas pequenas mudanças nos limites entre as camadas usando uma aproximação linearizada. Este processo é interativo obtendo sucessivos refinamentos até que um critério de convergência seja satisfeito. O algoritmo é aplicado em dados sintéticos e reais demonstrando a robustez do método.
Acesso aberto (Open Access)
Efeitos do manto de intemperismo sobre anomalias “VLF” dos corpos tabulares inclinados
(Universidade Federal do Pará, 1997-03-16) ROZAL, Edilberto Oliveira; VERMA, Om Prakash; http://lattes.cnpq.br/2723609019309173
A utilização dos métodos indutivos de propagação E.M. na exploração mineral em regiões tropicais, apresenta grandes dificuldades devido a presença de uma camada superficial condutiva (manto de intemperismo) comumente encontrada nestas regiões. Na região Amazônica, o manto apresenta-se bastante desenvolvido e condutivo, e em regiões semi-áridas, pode-se formar uma fina crosta superficial de sal. Em conseqüência disto, a interpretação dos dados E.M. obtidos para modelos que não consideram uma cobertura condutiva levam a erros consideráveis. Objetivando-se estudar os efeitos do manto sobre anomalias VLF devidas a corpos tabulares inclinados em contato com o manto (manto ohmico), foi realizada uma série de experimentos através do modelamento analógico, considerando-se diferentes parâmetros de resposta para o manto e o corpo. O manto de intemperismo foi simulado por soluções de cloreto de amônia (NH4Cl) dispostas horizontalmente e o corpo condutor por chapas de grafite colocadas em posições inclinadas verticalmente. Utilizou-se quatro corpos condutores e três mantos com diferentes espessuras e condutividades, simulando, desta forma, diversas situações geológicas. Os resultados são dados por simples situações dos corpos localizados em um meio não condutor (ar), onde os parâmetros variados são: profundidade do topo, condutividade e mergulho do corpo. Os efeitos da condutividade da cobertura são amplamente ilustrados e avaliados. Para a análise dos resultados, foi plotado um conjunto de curvas considerando-se os valores pico-a-pico das anomalias de "tilt angle" e de elipsidade. Os resultados foram sintetizados em um outro conjunto de curvas reunidas em diagramas de Argand. Estando ou não o manto presente, observou-se, tanto para o tilt angle quanto para a elipsidade, o efeito do aumento da profundidade é o de reduzir a magnitude pico-a-pico e a forma do pico da anomalia, fazendo com que este afaste-se do ponto de "cross-over". Para um condutor de mesma espessura, o aumento da condutividade causa um ligeiro aumento nas anomalias de tilt angle, e uma atenuação nas anomalias de elipsidade. O efeito geral na variação do mergulho do condutor é o de causar uma assimetria nos perfis de tilt angle e de elipsidade. O aumento da condutância do manto de intemperismo causa um acréscimo nas anomalias de elipsidade e uma ligeira diminuição nas anomalias de tilt angle; porém, a partir de um certo valor de condutância do manto (mantos mais condutivos) tanto as anomalias de tilt angle quanto as anomalias de elipsidade começam a atenuar. Há rotação de fase no sentido anti-horário, sendo mais intensa para grandes valores de número de indução do corpo. Na presença do manto, o corpo parece estar a uma profundidade inferior à verdadeira e a ser menos condutivo.
Acesso aberto (Open Access)
Eletroresistividade aplicada à caracterização hidrogeológica para o abastecimento da cidade de São Domingos do Araguaia-Pará
(Universidade Federal do Pará, 1996-02-13) SILVA, Arnaldo de Queiroz da; RIJO, Luiz; http://lattes.cnpq.br/3148365912720676
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 (pa), 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.
Acesso aberto (Open Access)
Eletroresistividade na prospecção hidrogeológica no município de Ourilândia do Norte-PA
(Universidade Federal do Pará, 1996-05-22) SOUSA, José Antonio Ferreira de; RIJO, Luiz; http://lattes.cnpq.br/3148365912720676
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.
Acesso aberto (Open Access)
Espalhamento elástico em meios anisotrópicos estratificados
(Universidade Federal do Pará, 1994-06-09) PROTÁZIO, João dos Santos; SCHOENBERG, Michael
A análise de AVO constitui-se, atualmente, numa importante ferramenta para a extração de informações litológicas a partir de dados sísmicos, através do uso dos contrastes de impedância acústica nas interfaces que separam diferentes litologias. A hipótese usual de isotropia deixa de valer, em muitos casos, após o advento de arranjos de grande afastamento e geofones com multi-superfície. Para a interpretação destes dados, a análise de AVO deve incluir anisotropia. Este trabalho apresenta uma teoria de AVO e resultados numéricos para um meio anisotrópico estratificado. Esta tese contém três contribuições. Inicialmente, é apresentada uma nova abordagem para o estudo da reflexão-transmissão através de interface plana que separam dois meios anisotrópicos com pelo menos um plano horizontal de simetria especular. As equações de Zoeppritz são generalizadas para incluir anisotropia, através da introdução das chamadas matrizes de impedância, o que simplifica bastante o formalismo anterior. Posteriormente, é descrito o estudo da reflexão de ondas P através de interface entre um meio isotrópico e outro transversalmente isotrópico (TI). É mostrado que a reflexão de ondas P, neste tipo de experimento, não fornece informações sobre a presença de anisotropia do semi-espaço TI, pelo menos em incidência pré-crítica. Finalmente, é discutido o comportamento da reflexão e transmissão de pulsos, em incidência pós-crítica, através de meios anisotrópicos estratificados. Observa-se que o comportamento pós-crítico dos pulsos espalhados carregam valiosa informações sobre a anisotropia dos meios atravessados por eles.
Acesso aberto (Open Access)
Estimativa de parâmetros elásticos em meios anisotrópicos
(Universidade Federal do Pará, 2003-06-20) GOMES, Ellen de Nazaré Souza; PROTÁZIO, João dos Santos; http://lattes.cnpq.br/4210442535067685
As medidas de amplitude, polarização e vagarosidade contem informações sobre o meio onde a propagação de onda ocorre. Esta tese investiga esses dados com objetivo de estimar as propriedades elásticas deste meio. Coeficientes de reflexão podem ser estimados das amplitudes dos dados e dependem de forma não linear dos contrastes dos parâmetros elásticos e do contraste de densidade entre os meios separados por uma interface. Quando o contraste de impedância é fraco, as aproximações lineares para a refletividade qP são mais convenientes para inversão da densidade e dos parâmetros elásticos usando as análises de amplitude versus ângulo de incidência (AVO) e amplitude versus a direção do plano de incidência (AVD). Escrevendo as equações de Zoepprittz de forma separada nos permite escrever uma solução destas equações em termos das matrizes de impedância e polarização. Usando esta solução são determinadas aproximações lineares para a refletividade da onda qP considerando fraco contraste de impedância, fraca anisotropia mas com classe de simetria de arbitrária. As linearizações são avaliadas para diferentes geometrias de aquisição e várias escolhas do meio de referência. Estas aproximações apresentam bom desempenho comparado com o valor exato do coeficiente de reflexão da onda qP e de suas ondas convertidas para incidências de até 30° e meios que obedecem à hipótese de fraca anisotropia. Um conjunto de fraturas orientado é representado efetivamente por um meio transversalmente isotrópico (TI), as aproximações lineares da refletividade da onda qP podem ser usadas para estimar a orientação de fratura. Partindo deste pressuposto este problema consiste em estimar a orientação do eixo de simetria a partir de dados de refletividade de onda qP. Este trabalho mostra que são necessários múltiplos azimutes e múltiplas incidências para se obter uma estimativa estável. Também é mostrado que apenas os coeficientes das ondas qS e qT são sensíveis ao mergulho da fratura. Foi investigada a estimativa da anisotropia local através de dados de VSP multiazimutal dos vetores de polarização e vagarosidade. Foram usadas medidas da componente vertical do vetor de vagarosidade e o vetor de polarização de ondas qP diretas e refletidas. O esquema de inversão é validado através de exemplos sintéticos considerando diferentes escolhas do vetor normal à frente de onda no meio de referência, meios de referências e geometria de aquisição. Esta análise mostra que somente um subgrupo dos parâmetros elástico pode ser estimado. Uma importante aplicação desta metodologia é o seu potencial para a determinação de classes de anisotropia. A aplicação desta metodologia aos dados do mar de Java mostra que os modelos isotrópicos e TIV são inadequados para o ajuste desses dados.
Acesso aberto (Open Access)
Estimativa de porosidade em lâminas petrográficas através da morfologia matemática binária
(Universidade Federal do Pará, 2013-08-02) CASTELO, Fernando Walleson Louzada; ANDRADE, André José Neves; http://lattes.cnpq.br/8388930487104926
Na exploração de petróleo em bacias submersas, ocorre a necessidade da perfuração de poços com alto ângulo de desvio e de poços horizontais, que visam otimizar o número de alvos exploratórios atingidos por uma única plataforma. Nestes casos, é tecnicamente impossível realização das operações de testemunhagem, o que inviabiliza a análise de testemunho para a medida direta da porosidade. Em algumas situações da avaliação de formação em que existe baixa confiança nos valores de porosidade, pode-se apoiar no conhecimento geológico da área. Para as bacias semi-submersas, a continuidade lateral das camadas pode permitir a coleta de amostras de mão em afloramentos, na parte emersa da bacia. No caso das bacias submersas, pode-se colher amostras de mão em afloramentos de formações análogas. Um problema relativamente comum em um laboratório de análises petrofísicas é a extração de amostras (plugs) adaptadas aos equipamentos de medida, diretamente das amostras de rocha recolhidas em afloramentos. Por outro lado, para este tipo de amostra de rocha é trivial a obtenção de laminas petrográficas. Esta dissertação busca produzir uma estimativa da porosidade diretamente sobre as imagens de laminas petrográficas com a utilização do método de processamento de imagens, conhecido como morfologia matemática, que objetiva descrever quantitativamente as estruturas (formas) geométricas presente nas imagens.
Acesso aberto (Open Access)
Estimativas da condutividade térmica dos minerais e rochas e influência de parâmetros térmicos e petrofísicos na resistividade aparente da formação
(Universidade Federal do Pará, 1995-08-09) COZZOLINO, Klaus; HOWARD JUNIOR, Allen Quentin; http://lattes.cnpq.br/6447166738854045
O presente estudo realiza estimativas da condutividade térmica dos principais minerais formadores de rochas, bem como estimativas da condutividade média da fase sólida de cinco litologias básicas (arenitos, calcários, dolomitos, anidritas e litologias argilosas). Alguns modelos térmicos foram comparados entre si, possibilitando a verificação daquele mais apropriado para representar o agregado de minerais e fluidos que compõem as rochas. Os resultados obtidos podem ser aplicados a modelamentos térmicos os mais variados. A metodologia empregada baseia-se em um algoritmo de regressão não-linear denominado de Busca Aleatória Controlada. O comportamento do algoritmo é avaliado para dados sintéticos antes de ser usado em dados reais. O modelo usado na regressão para obter a condutividade térmica dos minerais é o modelo geométrico médio. O método de regressão, usado em cada subconjunto litológico, forneceu os seguintes valores para a condutividade térmica média da fase sólida: arenitos 5,9 ± 1,33 W/mK, calcários 3.1 ± 0.12 W/mK, dolomitos 4.7 ± 0.56 W/mK, anidritas 6.3 ± 0.27 W/mK e para litologias argilosas 3.4 ± 0.48 W/mK. Na sequência, são fornecidas as bases para o estudo da difusão do calor em coordenadas cilíndricas, considerando o efeito de invasão do filtrado da lama na formação, através de uma adaptação da simulação de injeção de poços proveniente das teorias relativas à engenharia de reservatório. Com isto, estimam-se os erros relativos sobre a resistividade aparente assumindo como referência a temperatura original da formação. Nesta etapa do trabalho, faz-se uso do método de diferenças finitas para avaliar a distribuição de temperatura poço-formação. A simulação da invasão é realizada, em coordenadas cilíndricas, através da adaptação da equação de Buckley-Leverett em coordenadas cartesianas. Efeitos como o aparecimento do reboco de lama na parede do poço, gravidade e pressão capilar não são levados em consideração. A partir das distribuições de saturação e temperatura, obtém-se a distribuição radial de resistividade, a qual é convolvida com a resposta radial da ferramenta de indução (transmissor-receptor) resultando na resistividade aparente da formação. Admitindo como referência a temperatura original da formação, são obtidos os erros relativos da resistividade aparente. Através da variação de alguns parâmetros, verifica-se que a porosidade e a saturação original da formação podem ser responsáveis por enormes erros na obtenção da resistividade, principalmente se tais "leituras" forem realizadas logo após a perfuração (MWD). A diferença de temperatura entre poço e formação é a principal causadora de tais erros, indicando que em situações onde esta diferença de temperatura seja grande, perfilagens com ferramentas de indução devam ser realizadas de um a dois dias após a perfuração do poço.
Acesso aberto (Open Access)
Estudo da resposta de perfilagem de indução de camadas finas com diferentes arranjos de bobinas: modelamento analógico
(Universidade Federal do Pará, 1993-12-16) CARVALHO, Paulo Roberto de; VERMA, Om Prakash; http://lattes.cnpq.br/2723609019309173
Objetivando contribuir para a melhoria da resolução vertical das sondas de indução, utilizando arranjos de bobinas não-convencionais, fizemos um estudo comparativo das respostas obtidas com os arranjos coaxial e coplanar, através do modelamento analógico em escala reduzida. Construímos sondas de indução com um par de bobinas, bem como modelos geoelétricos que simulam seqüências litoestratigráficas formadas por camadas tanto espessas como delgadas, com ou sem invasão de fluidos, utilizando um fator de redução de escala igual a 20. O sistema de instrumentação nos permitiu medições da razão entre o campo secundário com relação ao primário na ordem de 0,01 %. Analisando os perfis obtidos com ambos os arranjos, coaxial e coplanar, chegamos a conclusão que: • quando se refere a camadas delgadas de condutividade elétrica relativamente elevadas, como é o caso de níveis argilíticos num pacote arenítico contendo hidrocarbonetos, o arranjo de bobinas coaxial é visivelmente superior ao coplanar, no que se refere ao posicionamento e estimativa das espessuras destas finas camadas; • por outro lado, quando se trata de camadas delgadas de condutividade relativamente baixa, como é o caso de lentes areníticas saturadas em hidrocarbonetos num pacote de folhelho, verificamos que o arranjo coplanar apresenta uma resolução vertical sensivelmente melhor, tanto para camadas finas quanto para as de maior espessura; • o efeito de camadas adjacentes (shoulder effect) se apresenta bem mais acentuado nos perfis obtidos com o arranjo coaxial; • o arranjo coplanar apresenta uma melhor definição de bordas para as camadas espessas. Entretanto, em camadas de menor espessura, o arranjo coplanar perde aquela ligeira oscilação do sinal que posiciona as interfaces de contato entre camadas.
Acesso aberto (Open Access)
Estudo de geração de campo uniforme para modelagem eletromagnética
(Universidade Federal do Pará, 1996-10-25) CHAVES NETO, Antonio Maia de Jesus; VERMA, Om Prakash; http://lattes.cnpq.br/2723609019309173
Para estudar os problemas de prospecção geofísica eletromagnética através de modelagem analógica, as condições em escala natural são representadas, no laboratório, em escala reduzida de acordo com a teoria da similitude. Portanto, para investigar os problemas de técnicas VLF, AFMAG e MT, frequentemente é necessário criar campo uniforme no arranjo experimental. Vários sistemas físicos para geração de campos uniformes são analisados teoricamente nesta tese. Os sistemas estudados aqui são a bobina circular, bobina de Helmholtz, solenóide, um plano de corrente e dois planos paralelos de correntes. As equações analíticas foram obtidas para campo magnético num ponto do espaço e subsequentemente as condições de campo uniforme. Nos casos em que as condições para o campo uniforme não puderam ser obtidas analiticamente, a porcentagem de desvio do campo em relação a um ponto pré-selecionado foi calculada. Contudo, os mapas de campo magnético, assim como o mapa de porcentagem de desvio, estão presentes para todos os sistemas estudados aqui. Também, foram calculados as áreas e os volumes espaciais de vários desvios de porcentagem do campo uniforme. Um estudo comparativo desses sistemas mostra que o solenóide é a maneira mais eficiente para criar um campo uniforme, seguido pelo sistema de bobinas de Helmholtz. Porém, o campo criado em um solenóide está em um espaço fechado onde é difícil colocar modelos e substituí-los para executar experimentos. Portanto, recomenda-se o uso de bobinas de Helmholtz para criar um campo uniforme. Este último sistema fornece campo uniforme com espaço aberto suficiente, o que facilita o experimento.