Identificação de ambientes de sedimentação na Área Metropolitana de Belém a partir de perfis de poço

Esta dissertação apresenta um estudo realizado com base em padrões associados à forma das curvas de perfis geofísicos de poços perfurados para captação de água subterrânea na Região Metropolitana de Belém. O objetivo do estudo foi o reconhecimento de depósitos cenozóicos, bem como a identificação de seus ambientes de sedimentação. Para isto, foram utilizados os perfis de Raios-Gama, Resistência Elétrica e Potencial Espontâneo corridos em 21 poços, cobrindo a profundidade máxima de 300 m. O estudo permitiu identificar dois pacotes de sedimentos: um até a profundidade de 160 m, contendo bastante argila e finas camadas de areia, e o outro abaixo de 160 m, com espessas camadas arenosas. Esses pacotes foram correlacionados com seqüências deposicionais do Cenozóico Inferior descritas por Rossetti (2000), as quais correspondem às formações Pirabas, Barreiras e aos sedimentos Pós-Barreiras e Holocênicos. No pacote superior, a forma das curvas dos perfis mostra padrões que podem ser associados a seqüências transgressivas e regressivas, enquanto no pacote inferior a forma das curvas sugere um ambiente de centro de canal. O estudo demonstra que a metodologia empregada para o reconhecimento dos ambientes de sedimentação pode ser considerada satisfatória quando não há nenhuma outra informação geológica adicional disponível.

Identificação de fácies em perfis de poço com algoritmo inteligente

A identificação de fácies em um poço não testemunhado é um dos problemas clássicos da avaliação de formação. Neste trabalho este problema é tratado em dois passos, no primeiro produz-se a codificação da informação geológica ou da descrição das fácies atravessadas em um poço testemunhado em termos das suas propriedades físicas registradas nos perfis geofísicos e traduzidas pelos parâmetros L e K, que são obtidos a partir dos perfis de porosidade (densidade, sônico e porosidade neutrônica) e pela argilosidade (Vsh) calculada pelo perfil de raio gama natural. Estes três parâmetros são convenientemente representados na forma do Gráfico Vsh-L-K. No segundo passo é realizada a interpretação computacional do Gráfico Vsh-L-K por um algoritmo inteligente construído com base na rede neural competitiva angular generalizada, que é especializada na classificação de padrões angulares ou agrupamento de pontos no espaço n-dimensional que possuem uma envoltória aproximadamente elipsoidal. Os parâmetros operacionais do algoritmo inteligente, como a arquitetura da rede neural e pesos sinápticos são obtidos em um Gráfico Vsh-L-K, construído e interpretado com as informações de um poço testemunhado. Assim, a aplicação deste algoritmo inteligente é capaz de identificar e classificar as camadas presentes em um poço não testemunhado, em termos das fácies identificadas no poço testemunhado ou em termos do mineral principal, quando ausentes no poço testemunhado. Esta metodologia é apresentada com dados sintéticos e com perfis de poços testemunhados do Campo de Namorado, na Bacia de Campos, localizada na plataforma continental do Rio de Janeiro, Brasil.

Interpretação do gráfico de Hingle através de rede competitiva angular

A saturação de água é a principal propriedade petrofísica para a avaliação de reservatórios de hidrocarbonetos, pois através da análise dos seus valores é definida a destinação final do poço recém perfurado, como produtor ou poço seco. O cálculo da saturação de água para as formações limpas é, comumente, realizado a partir da equação de Archie, que envolve a determinação da resistividade da zona virgem, obtida a partir de um perfil de resistividade profunda e o cálculo de porosidade da rocha, obtida a partir dos perfis de porosidade. A equação de Archie envolve ainda, a determinação da resistividade da água de formação, que normalmente necessita de definição local e correção para a profundidade da formação e da adoção de valores convenientes para os coeficientes de Archie. Um dos métodos mais tradicionais da geofísica de poço para o cálculo da saturação de água é o método de Hingle, particularmente útil nas situações de desconhecimento da resistividade da água de formação. O método de Hingle estabelece uma forma linear para a equação de Archie, a partir dos perfis de resistividade e porosidade e a representa na forma gráfica, como a reta da água ou dos pontos, no gráfico de Hingle, com saturação de água unitária e o valor da resistividade da água de formação é obtido a partir da inclinação da reta da água. Independente do desenvolvimento tecnológico das ferramentas de perfilagem e dos computadores digitais, o geofísico, ainda hoje, se vê obrigado a realizar a interpretação de ábacos ou gráficos, sujeito a ocorrência de erros derivados da sua acuidade visual. Com o objetivo de mitigar a ocorrência deste tipo de erro e produzir uma primeira aproximação para a saturação de água em tempo real de perfilagem do poço, insere-se o trabalho apresentado nesta dissertação, com a utilização de uma conveniente arquitetura de rede neural artificial, a rede competitiva angular, capaz de identificar a localização da reta da água, a partir da identificação de padrões angulares presentes nos dados dos perfis de porosidade e resistividade representados no gráfico de Hingle. A avaliação desta metodologia é realizada sobre dados sintéticos, que satisfazem integralmente a equação de Archie, e sobre dados reais.

Interpretação de perfis dos carbonatos fraturados da Bacia do Pará-Maranhão

A finalidade deste trabalho é apresentar um estudo de caso dos carbonatos terciários da bacia do Pará-Maranhão, do ponto de vista da interpretação dos perfis registrados nessa área. Dois poços-chave, X e Y, foram escolhidos para o estudo. O tratamento dos dados foi realizado utilizando os recursos do programa LOGCALC, instalado no computador IBM-3090, da Petrobrás, e, também, o sistema DLPS, instalado no VAX-8600, da Universidade Federal de Pará. A avaliação da porosidade e, principalmente, das saturações, é dificultada pelas características não-convencionais dos carbonatos. A litologia é complexa, a salinidade da água de formação é relativamente baixa, em torno de 10.000 ppm de NaCl, a densidade das rochas é elevada e os carbonatos estão fraturados. Para diminuir o efeito da composição mineralógica, foi necessário discriminar os diferentes tipos litológicos. Três tipos de carbonatos foram identificados: calcário, calcário arenoso e calcário dolomitizado. A identificação litológica permite maior controle dos parâmetros da matriz e dos expoentes de porosidade das rochas, conduzindo a estimativas de porosidade e de saturações mais confiáveis. A presença de fraturas influencia, marcadamente, a resposta dos perfis, conforme pode ser notado no perfil de identificação de fraturas, nos perfis de resistividade, de densidade e na curva de raios-gama espectral do poço X. O expoente de porosidade, m, tomado do gráfico de Pickett, é frequentemente inferior a 1,5, valor considerado inerente a rochas fraturadas. Os modelos de Rasmus (1983) e de Porter et al.(1969), foram testados para calcular as saturações de água e óleo. A equação tradicional de Archie (1942), com a e m apropriados, também foi aplicada neste trabalho. O esquema para estimativa de saturações que conduz a resultados mais coerentes com os dados dos testes de formação, nos intervalos fraturados, é a saturação de Rasmus, calculada com a relação de Archie aplicada com o expoente de porosidade variável, tomado do modelo de Rasmus.

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.

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.

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.

Melhoramento da resolução para camadas delgadas de perfilagens existentes em áreas específicas de produção-exploração

No campo da perfilagem de poço existe um grande interesse voltado para as técnicas que melhorem a resolução vertical das respostas das ferramentas. Neste trabalho, optamos por desenvolver uma técnica de aumento de resolução vertical dos perfis obtidos com a ferramenta de indução denominada na bibliografia de 6FF40, através de um algoritmo que utiliza as técnicas de estimativa dos mínimos quadrados. Este método torna possível o processamento dos dados registrados de uma maneira computacionalmente eficiente. O algoritmo apresentado necessita apenas dos dados registrados, dos valores adotados como coeficientes da ferramenta, e de uma estimativa dos ruídos existentes. Como ilustração foram utilizados trechos de um perfil de uma área que reconhecidamente apresenta problemas relacionados à resolução vertical da ferramenta 6FF40. Com o objetivo de verificar a eficiência do algoritmo utilizado, os perfis processados foram correlacionados com os perfis de raio gama e com os perfis esféricos focalizados do poço 7-LOR-18-RN, perfis estes que possuem resolução vertical maior do que os perfis convencionais de indução, comprovando a individualização das camadas delgadas que antes do processamento não eram facilmente reconhecidas no perfil original. O algoritmo foi também testado com dados sintéticos demonstrando sua eficiência na recuperação de valores mais representativos para R_t.

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.

Contaminação por herbicidas em corpos hídricos da microbacia do Córrego Rico (SP) e aspectos toxicológicos de atrazine a juvenis de Piaractus mesopotamicus

Pós-graduação em Agronomia (Produção Vegetal) - FCAV

Modelamento e correção de descentralização das imagens de tempo de trânsito

As imagens da parede do poço obtidas através da ferramenta UBI (Ultrasonic Borehole Imager) são amplamente utilizadas por geólogos e engenheiros de petróleo para identificar eventos geológicos nas paredes de poços abertos e na inspeção dos tubos de revestimento, uma vez que este perfil praticamente fornece uma fotografia da parede do poço. As ferramentas de imageamento acústico produzem imagens do tempo de trânsito e da amplitude do pulso acústico gerado pela ferramenta e refletido na parede do poço. Entretanto, estas imagens podem ter uma interpretação não realista, uma vez que elas podem ter seu aspecto alterado em razão do movimento da ferramenta no interior do poço. Este trabalho apresenta o modelamento das imagens de tempo de trânsito a partir da aplicação do critério de Coulomb para a ruptura da parede do poço submetida a um estado plano de tensões, a qual fornecerá a seção do poço, que é a forma geométrica que será traçada pela ferramenta de imageamento acústico do poço. O deslocamento ascensional da ferramenta e as imperfeições da parede do poço, normalmente, são os responsáveis pelo deslocamento da posição do transdutor em relação ao eixo do poço. Este efeito pode ter grande responsabilidade nas imperfeições das imagens de tempo de trânsito. Assim, a correção dessas imagens, chamada de correção da descentralização, busca o reposicionamento do transdutor para a posição do eixo do poço. Apresenta-se, também, um método de correção do efeito da descentralização da ferramenta baseado neste modelamento. O método é proposto com base na geometria analítica plana e no método do raio para a definição do tempo de trânsito do pulso acústico, com o objetivo de reconstruir as imagens de tempo de trânsito obtidas com a ferramenta descentralizada, ou seja, corrigir estas imagens tornando-as como se fossem obtidas com a ferramenta centralizada em relação ao eixo do poço.

Modelamento de perfis de indução

Neste trabalho foi implementado um algoritmo que permite o modelamento das respostas das sondas de indução (2C40-6FF40-ILD) em meios estratificados. O algoritmo é baseado na solução exata do problema de Green para a obtenção do campo eletromagnético em meios estratificados com a fonte no interior dos estratos. O código computacional desenvolvido neste trabalho para meios horizontalmente estratificados tem o objetivo de auxiliar o desenvolvimento de códigos mais avançados, uma vez que é difícil o acesso para a pesquisa dos códigos citados na literatura. Os códigos disponíveis são comerciais, mas estes têm o objetivo único de produção, além de serem, bastante dispendiosos. O código também foi utilizado para a obtenção da resposta da sonda ILD em reservatórios do tipo laminado, onde a sonda tem baixa resolução. E verificamos, a exemplo de outros autores, que a condutividade obtida pela sonda pode ser expressa pela soma das condutividades dos constituintes litológicos do reservatório ponderada pelos seus percentuais em volume. Esta relação foi definida para reservatórios ainda não analisados na literatura, tipo arenito moderadamente silicificado/arenito/folhelho, onde a lâmina de interesse tem resistividade superior a resistividade de uma das lâminas, e inferior a resistividade da outra. Foi analisado também reservatórios tipo arenito moderadamente silicificado/arenito, onde a lâmina de interesse tem a menor resistividade. Com esta relação pode-se obter para estes reservatórios laminados estimativas mais precisas de Rt e consequentemente determinar as saturações água/óleo e produtibilidade, também de forma mais precisa.

Modelamento da permissividade dielétrica de rochas saturadas de óleo e água e suas aplicações em perfilagem de poço

A ferramenta de propagação eletromagnética (EPT) fornece o tempo de propagação (T_pl) e a atenuação (A) de uma onda eletromagnética que se propaga num meio com perdas. Estas respostas da EPT são funções da permissividade dielétrica do meio. Existem vários modelos e fórmulas de misturas sobre a permissividade dielétrica de rochas reservatório que podem ser utilizados na interpretação da ferramenta de alta frequência. No entanto, as fórmulas de mistura não consideram a distribuição e a geometria do espaço poroso, e estes parâmetros são essenciais para que sejam obtidas respostas dielétricas mais próximas de uma rocha real. Foi selecionado um modelo baseado nos parâmetros descritos acima e este foi aplicado à dados dielétricos disponíveis na literatura. Foi obtida uma boa concordância entre as curvas teóricas e os dados experimentais, comprovando assim que a distribuição e a geometria dos poros têm que ser levadas em conta no desenvolvimento de um modelo realista. Foram conseguidas também funções de distribuição de razão de aspecto de poros, através das quais geramos várias curvas relacionando as respostas da EPT com diversas saturações de óleo/gás. Estas curvas foram aplicadas na análise de perfis. Como o modelo selecionado ajusta-se bem aos dados dielétricos disponíveis na literatura, torna-se atraente aplicá-lo à dados experimentais obtidos em rochas de campos brasileiros produtores de hidrocarbonetos para interpretação da EPT corrida em poços destes campos petrolíferos.

Resposta da bananeira Prata-Anã, à aplicação de zinco e boro no rizoma

Pós-graduação em Agronomia (Produção Vegetal) - FCAV