Imageamento da porosidade através de perfis geofísicos de poço

O imageamento da porosidade é uma representação gráfica da distribuição lateral da porosidade da rocha, estimada a partir de dados de perfis geofísicos de poço. Apresenta-se aqui uma metodologia para produzir esta imagem geológica, totalmente independente da intervenção do intérprete, através de um algoritmo, dito, interpretativo baseado em dois tipos de redes neurais artificiais. A primeira parte do algoritmo baseia-se em uma rede neural com camada competitiva e é construído para realizar uma interpretação automática do clássico gráfico o P_b - Φ_N, produzindo um zoneamento do perfil e a estimativa da porosidade. A segunda parte baseia-se em uma rede neural com função de base radial, projetado para realizar uma integração espacial dos dados, a qual pode ser dividida em duas etapas. A primeira etapa refere-se à correlação de perfis de poço e a segunda à produção de uma estimativa da distribuição lateral da porosidade. Esta metodologia ajudará o intérprete na definição do modelo geológico do reservatório e, talvez o mais importante, o ajudará a desenvolver de um modo mais eficiente as estratégias para o desenvolvimento dos campos de óleo e gás. Os resultados ou as imagens da porosidade são bastante similares às seções geológicas convencionais, especialmente em um ambiente deposicional simples dominado por clásticos, onde um mapa de cores, escalonado em unidades de porosidade aparente para as argilas e efetiva para os arenitos, mostra a variação da porosidade e a disposição geométrica das camadas geológicas ao longo da seção. Esta metodologia é aplicada em dados reais da Formação Lagunillas, na Bacia do Lago Maracaibo, Venezuela.

Determinação automática da porosidade e zoneamento de perfis através da rede neural artificial competitiva

Duas das mais importantes atividades da interpretação de perfis para avaliação de reservatórios de hidrocarbonetos são o zoneamento do perfil (log zonation) e o cálculo da porosidade efetiva das rochas atravessadas pelo poço. O zoneamento é a interpretação visual do perfil para identificação das camadas reservatório e, consequentemente, dos seus limites verticais, ou seja, é a separação formal do perfil em rochas reservatório e rochas selante. Todo procedimento de zoneamento é realizado de forma manual, valendo-se do conhecimento geológico-geofísico e da experiência do intérprete, na avaliação visual dos padrões (características da curva do perfil representativa de um evento geológico) correspondentes a cada tipo litológico específico. O cálculo da porosidade efetiva combina tanto uma atividade visual, na identificação dos pontos representativos de uma particular rocha reservatório no perfil, como a escolha adequada da equação petrofísica que relaciona as propriedades físicas mensuradas da rocha com sua porosidade. A partir do conhecimento da porosidade, será estabelecido o volume eventualmente ocupado por hidrocarboneto. Esta atividade, essencial para a qualificação de reservatórios, requer muito do conhecimento e da experiência do intérprete de perfil para a efetiva avaliação da porosidade efetiva, ou seja, a porosidade da rocha reservatório, isenta do efeito da argila sobre a medida das propriedades físicas da mesma. Uma forma eficiente de automatizar estes procedimentos e auxiliar o geofísico de poço nestas atividades, que particularmente demandam grande dispêndio de tempo, é apresentado nesta dissertação, na forma de um novo perfil, derivado dos perfis tradicionais de porosidade, que apresenta diretamente o zoneamento. Pode-se destacar neste novo perfil as profundidades do topo e da base das rochas reservatório e das rochas selante, escalonado na forma de porosidade efetiva, denominado perfil de porosidade efetiva zoneado. A obtenção do perfil de porosidade efetiva zoneado é baseado no projeto e execução de várias arquiteturas de rede neural artificial, do tipo direta, com treinamento não supervisionado e contendo uma camada de neurônios artificiais, do tipo competitivo. Estas arquiteturas são projetadas de modo a simular o comportamento do intérprete de perfil, quando da utilização do gráfico densidade-neutrônico, para as situações de aplicabilidade do modelo arenito-folhelho. A aplicabilidade e limitações desta metodologia são avaliadas diretamente sobre dados reais, oriundos da bacia do Lago Maracaibo (Venezuela).

Estimativa dos perfis de permeabilidade e de porosidade utilizando rede neural artificial

A permeabilidade e a porosidade são duas das mais importantes propriedades petrofísicas para a qualificação dos reservatórios de óleo e gás. A porosidade está relacionada à capacidade de armazenamento de fluidos e a permeabilidade, com a capacidade de produção destes fluidos. Suas medidas são, normalmente, realizadas em laboratório, através de testemunhos da rocha. Esses processos têm custos elevados e nem todos os poços são testemunhados. As estimativas da permeabilidade e da porosidade são de fundamental importância para os engenheiros de reservatório e geofísicos, uma vez que seus valores podem definir a completação ou não de um poço petrolífero. O perfil de porosidade e sua relação com o perfil de densidade, é bem conhecida na geofísica de poço. No entanto, existem poucas relações quantitativas e/ou qualitativas entre a porosidade e a permeabilidade, como por exemplo as relações de Kozeny. Sendo assim, este trabalho busca o estabelecimento do perfil de permeabilidade e do perfil de porosidade, a partir de informações do perfil de densidade. Para tanto, buscamos a relação entre a propriedade física da rocha (densidade) e as propriedades petrofísicas: permeabilidade e porosidade, utilizando como metodologia à técnica de redes neurais artificiais, como a rede neural artificial com função de base radial. A obtenção da permeabilidade e da porosidade a partir da rede neural artificial, que possui como entrada a informação da densidade possibilita um menor custo para a aquisição dessas importantes informações petrofísicas, permite ao intérprete de perfis de poço optar ou não pela exploração de uma unidade estudada, além de uma visão mais completa do reservatório. Os procedimentos para a estimativa da permeabilidade e da porosidade estão direcionados para uma única formação, mas os intérpretes de perfis poderão aplicar a diretriz apresentada no programa de rede neural artificial com função de base radial, utilizando a estimativa dessas propriedades petrofísicas para outras formações, inclusive de outros campos petrolíferos. Portanto, recomenda-se a utilização de um conjunto de dados completo, com quantidade de dados suficientes de um mesmo poço, a fim de viabilizar corretamente a melhor interpretação.

Reciclagem de resíduo a partir de cinzas de carvão mineral: produção de argamassas

Este trabalho tem como objetivo avaliar a possibilidade de aproveitamento do resíduo obtido a partir da combustão do carvão mineral em caldeiras de leito fluidizado, de uma refinaria de alumina no estado do Pará. Neste contexto, foi avaliada a incorporação de cinzas volantes, como pozolana, em substituição parcial do cimento na produção de argamassas. Para tanto, foram elaborados corpos de prova utilizando-se o cimento do tipo Portland CPII-E-32. As misturas foram definidas na proporção 4:1, ou seja, 4 partes de agregado (sílica) e 1 parte de aglomerante (cimento e cinza), com a inserção de cinza nos teores de 10, 20, 30, 40 e 50 em porcentagem com relação ao cimento, além da argamassa sem adição. Foram definidas duas relações água/aglomerante, 0,4 e 0,8. Após tempo de cura de 7 e 28 dias, foram feitos testes de resistência à compressão para análise de comportamento. Além disso, as argamassas com a relação água/aglomerante 0,8 foram sujeitas a análises complementares de Difração de Raios-X, Microscopia Eletrônica de Varredura, absorção de água, porosidade aparente e massa específica aparente, cujos resultados obtidos mostraram-se compatíveis quando comparados com os dados da literatura, demonstrando ser viável a aplicação das cinzas estudadas na indústria da construção civil.

Utilização do lodo de estação de tratamento de efluentes da indústria de injeção plástica como matéria-prima para indústria cerâmica

A indústria de injeção plástica é de fundamental importância para a economia nacional, pois é responsável pela produção de insumos primários para várias indústrias de eletroeletrônicos. Estas empresas geram efluentes procedentes das unidades de pintura por jateamento e aspersão. Este resíduo necessita de tratamento especial para sua separação, armazenagem e destino final. Este trabalho estuda uma forma alternativa de utilização deste resíduo na indústria de cerâmica vermelha, avaliando parâmetros que influenciam em propriedades importantes nos artefatos fabricados com argila vermelha como: densidade, porosidade e resistência à flexão, avaliados para misturas de 2%, 5% e 10% em peso de resíduo, em relação a argila utilizada, os resultados mostram que há influência deste material nas reações de estado sólido podendo haver melhora nas propriedades citadas.

Estudo e avaliação do uso de resíduos do processo Bayer como matéria-prima na produção de agregados sintéticos para a construção civil

A grande quantidade de resíduos sólidos, oriundos da indústria de mineração e metalurgia, constitui um sério problema sócio-ambiental. As características físicoquímicas destes resíduos despertam grande interesse para outras indústrias. A indústria cerâmica pode ser uma interessante consumidora da maioria deste material, para suprir a grande escassez das reservas de matérias-primas atual. Neste contexto, este trabalho mostra os estudos realizados para a reciclagem da lama vermelha, como matéria-prima na produção de agregado sintético, visando à construção civil. A lama vermelha, principal rejeito industrial da fabricação de alumina, mostrou-se um insumo de grande interesse na fabricação de diferentes tipos de agregado para ser utilizado na produção de concreto, para a construção civil. Pelas suas características físico-químicas e a grande quantidade que é produzida anualmente (cerca de 10.000.000 t em duas fábricas, só no Norte do Brasil). Estudos realizados neste trabalho mostram a possibilidade de fabricação de agregados, com diferentes propriedades e possibilidades de aplicação, na indústria da construção civil. Estas propriedades dependem do controle de parâmetros, como o teor de sílica livre e argila, a granulometria e a temperatura de sinterização. Tais variáveis permitem controlar perfeitamente a formação de fase vítrea que é a responsável pelas propriedades dos agregados como: porosidade, resistência mecânica e densidade. Este material pode ser utilizado em concreto convencional ou em concretos especiais, para atender a demanda da construção civil.

Produção de carvão ativado a partir da biomassa residual da castanha do Brasil (Bertholletia excelsa L.) para adsorção de cobre (II)

Produziu-se carvão ativado a partir da casca da castanha-do-Brasil (Bertholletia excelsa L.) para ser utilizado na remoção de cobre (II), no processo de adsorção em sistema de batelada. A casca é resíduo do beneficiamento da castanha, que foi coletada, selecionada, lavada em água corrente e depois secada em estufa a 150ºC por 24 h. Os carvões foram carbonizados a 400 ºC por 3 h e ativados termicamente a 800ºC em tempos de 1, 2 e 3 h, quando receberam as respectivas codificaçõesCA1, CA2 e CA3. Depois foram caracterizadas quanto: à área superficial específica, ao volume e tamanho de poros, à difração de raios-X, à microscopia eletrônica de varredura (MEV) acoplada ao EDS(sistema de energia dispersiva por raios-X) e à espectroscopia de infravermelho com transformada de Fourier (FTIR). Ensaios preliminares foram realizados para avaliar a eficiência dos carvões quanto à remoção de cobre (II) em solução sintética de concentração inicial 50 mg L^-1. Como os resultados foram satisfatórios para CA1, CA2 e CA3 (93,43, 97,23 e 96, 92 % para os respectivos carvões), decidiu-se pelo que apresentou maior percentual de remoção. O CA2 foi produzido e caracterizado quanto: às densidades reais e aparentes, à porosidade em leito fixo, ao pH, à umidade (em base úmida), às cinzas, ao carbono fixo e aos grupos funcionais pelo método de Boehm. Realizaram-se ensaios para determinar a eficiência de remoção de cobre (II) quanto à influência do diâmetro da partícula do carvão, ao pH da solução, à influência do tempo de contato e à variação da concentração inicial. Os resultados de maior percentual de remoção foram para o diâmetro 0,595≤D≤1,19mm, pH 5,09, tempo de 5 min e concentrações de 50, 100 e 150mg L^-1.O modelo cinético de adsorção que melhor se ajustou aos dados foi o de pseudo2^a ordem. Os dados experimentais apresentaram bom ajuste aos modelos matemáticos de isotermas de Langmuir e Frendlich. Sendo assim, obteve-se carvão ativado de baixo custo a partir da casca da castanha-do-Brasil a qual apresentou boa eficiência na remoção de cobre (II) (acima de 90 % para maioria das análises) possibilitando também a utilização no tratamento de efluentes.

Obtenção e caracterização de carvão ativado de caroço de buriti (Mauritia flexuosa L. f.) para avaliar o processo de adsorção de uma solução de Cu (II)

Produziu-se carvão de caroço de buriti (CCB) a partir do rejeito da produção artesanal do fruto, destinado à extração de óleo, à temperatura de 400ºC. O CCB então foi ativado a temperaturas de 800ºC e 900ºC. Ensaios de adsorção foram executados para se avaliar o desempenho dessas temperaturas de ativação na adsorção de uma solução de cobre (II) de concentração inicial conhecida de 50 mg/L. Após a análise dos resultados, decidiu-se pela ativação do carvão a 900ºC. Caracterizou-se o carvão ativado do caroço de buriti (CACB) a 900ºC de acordo com as propriedades comerciais, tais como: área específica, porosimetria, densidades aparente e real, porosidade de um leito fixo, microscopia eletrônica de varredura, conteúdo de cinzas, pH, umidade, carbono fixo e grupos funcionais de superfície ácida presentes no CACB. Realizaram-se ensaios de equilíbrio de adsorção para se determinar a influência do diâmetro das partículas de CACB (D < 0,595 mm; 0,595 < D < 1,19 mm e D > 1,19 mm); a influência do tempo de contato adsorvente/adsorvato (15, 30, 60, 120, 180, 240 e 300 minutos); a influência do pH (3,00; 4,01; 5,18 e 6,00) e a influência da concentração inicial da solução de cobre (II) (5, 10, 30, 50, 80, 100 e 200 mg/L) para se avaliar a eficiência de remoção. Os resultados mostraram uma maior eficiência de remoção de cobre (II) para o diâmetro D < 0,595 mm; para o tempo de contato de 300 minutos; para o pH de 4,01 e para as concentrações iniciais de cobre (II) de 50 e 80 mg/L. Os modelos matemáticos de Langmuir e Freundlich foram aplicados para os dados de equilíbrio de adsorção. O modelo matemático de Langmuir foi o que melhor se ajustou aos dados de equilíbrio. De acordo com os dados da cinética de equilíbrio, observa-se que a partir do tempo de contato de 15 minutos todas as concentrações de equilíbrio ficaram abaixo do máximo permitido de 1,0 mg/L previsto pela legislação CONAMA nº 357/2005 para lançamento de efluentes em corpos receptores. Os resultados experimentais encontrados são indicativos de que é possível a remoção de cobre (II) de efluentes industriais utilizando CACB ativado fisicamente a 900ºC por um período de 60 minutos.

Projeto de um controlador Fuzzy Takagi-Sugeno para uma estação de tratamento de esgoto por lodos ativados

Este trabalho propõe o projeto de um controlador Fuzzy do tipo Takagi-Sugeno em uma estação de tratamento de esgoto por lodos ativados. Este tipo de tratamento ocorre na presença de oxigênio, pois microorganismos aeróbios presentes no licor misto irão proporcionar a remoção tanto da matéria carbonácea quanto nutrientes formados por compostos à base de nitrogênio. O controlador atua via mecanismo de aeração e foi projetado para interpolar os ganhos proporcionais e integrativos de três controladores fixos locais que, por sua vez, foram projetados a partir da linearização de um modelo contínuo de balanço de massa de oxigênio. O controlador auxiliará manter a concentração de oxigênio dissolvido desejável na faixa de operação do processo para que ocorra a nitrificação, reação química de oxiredução que transforma amônio em nitrito, que é fundamental para o sucesso do processo, e poderá também economizar energia elétrica utilizada pelo mecanismo de aeração. Os ensaios foram realizados via simulação computacional em quatro cenários idealizados comparando o desempenho do controlador fuzzy e o desempenho de um controlador fixo projetado em um ponto de operação diferente do qual se utilizou como valor de referência na planta. Ao aplicar um degrau no valor de referência estabelecido no cenário de análise, observou-se o tempo de subida, o tempo de acomodação e o erro em regime. Após os ensaios, observou-se que o desempenho do controlador fuzzy nos quesitos tempo de subida e tempo de acomodação em relação ao controlador fixo foi melhor, enquanto nos quesitos sobressinal e erro em regime foi semelhante. Após as análises, concluiu-se que a estratégia de controle escolhida neste trabalho é viável, pois de acordo com o valor auferido de oxigênio dissolvido na entrada da planta, o controlador fuzzy irá interpolar os ganhos proporcionais e integrativos de um controlador fixo projetado na vizinhança deste valor e assim, atuar de forma bastante satisfatória.

Estudo da migração internacional na Amazônia brasileira a partir da análise dos dados de espacialidade e seletividade dos censos de 2000 e 2010

A migração internacional tornou-se um tema de grande importância no cenário mundial tanto para os países de destino quanto para os países de origem. Na Amazônia ainda existe certa falta de interesse com a questão. No entanto, já existem alguns trabalhos específicos de cada país amazônico sobre o assunto. Porém as informações tratam estudos de caso sem levar em consideração a análise de toda a região. Isto se dá entre outros motivos pela falta de um banco de dados sobre migrações para a Amazônia. Entretanto os censos demográficos tornam-se importantes instrumentos de análise. Com base nos censos demográficos brasileiros de 2000 e 2010 que se busca questionar a dinâmica migratória recente para a Amazônia brasileira a fim de dar prosseguimento a uma série de discussões sobre as mudanças de origem, a distribuição espacial e o perfil desse novo migrante. Em um primeiro momento cria-se uma breve revisão teórica sobre a migração internacional, posteriormente discute-se a história da migração internacional na Amazônia brasileira e logo após analisam-se os dados dos dois últimos censos demográficos brasileiros sobre a migração internacional Num segundo momento analisa-se a distribuição espacial da migração de brasileiros retornados a Amazônia brasileira com intuito de fazer uma comparação entre os municípios que evidenciam essa dinâmica tanto para estrangeiros quanto para brasileiros. Essa espacialização retrata o processo histórico de ocupação e urbanização da Amazônia e também a porosidade da fronteira. Assim como o perfil migratório dá indícios de ser resultante do crescimento econômico da exploração de recursos naturais na Amazônia. Contudo, diante da complexidade do tema migração, é necessário integrar conhecimentos para a análise do processo migratório que ultrapassem os limites dos paradigmas clássicos.

Confecção de discos de desbaste a partir de AI2O3 e pó de vidro

Este trabalho descreve o desenvolvimento de discos abrasivos e de corte confeccionados com Al²O³ e pó de vidro reciclado para aplicação em ferramenta de corte. Os discos são utilizados para micromizar termoplásticos pós-consumo rígidos e levá-los a diversos tipos de processamento, como por exemplo, a rotomoldagem. Este disco abrasivo utiliza pó de vidro como um aditivo ligante e fundente, e pó de poli(tereftalato de etileno) pós-consumo como plastificante, facilitando um melhor preenchimento do molde. As proporções de pó de vidro reciclado foram analisadas a partir de estudos anteriores com outros minerais, como o diopsídio e o feldspato, já que os estudos com o pó de vidro reciclado substituindo os minerais naturais são recentes. O desenvolvimento dos discos acontece em duas etapas. A primeira é de homogeneização dos componentes com os percentuais calculados. Posteriormente é prensado e levado a mufla até uma temperatura de 900 ^oC para obter a percolação da resina e a consequente degradação do poli (tereftalato de etileno), pré-sinterização e amolecimento dos silicatos . Após a obtenção destes compostos verdes acontece a segunda e última etapa, que consiste na sinterização dos corpos de prova para a diminuição da porosidade e consequente aumento da resistência mecânica, seguida de acabamento superficial e testes de corte.

Modelo fractal para resistividade complexa de rochas: interpretação petrofísica e aplicação à exploração geoelétrica

Rochas contendo metálicos disseminados ou partículas de argila em ambiente natural onde soluções eletrolíticas normalmente preenchem os poros das rochas, exibem um tipo de polarização em baixas freqüências conhecido como polarização induzida. Nesta tese foi desenvolvido um novo modelo para descrever o fenômeno de polarização das rochas, não apenas em baixas freqüências, mas compreendendo todo o espectro eletromagnético, possível de utilização na prospecção geoelétrica. Este novo modelo engloba a maioria dos modelos utilizados até o momento como casos especiais, além de superar as limitações dos mesmos. Seu circuito analógico inclui uma impedância não linear do tipo r (iwtf)^-n que simula o efeito das superfícies rugosas das interfaces entre os grãos bloqueadores (partículas metálicas e/ou de argilas) e o eletrólito. A impedância de Warburg generalizada está em série com a resistência dos grãos bloqueadores da passagem de corrente e em paralelo com a impedância da dupla camada associada a essas interfaces. Esta combinação está em série com a resistência do eletrólito nas passagens dos poros bloqueados. Os canais não bloqueados são representados por uma resistência que corresponde à resistividade normal CC da rocha. A combinação desta resistência com a capacitância "global" da rocha é finalmente conectada em paralelo ao resto do circuito mencionado acima. Os parâmetros deste modelo incluem a resistividade CC (p₀), a cargueabilidade (m), três tempos de relaxação (t, Tf and T2), um fator de resistividade de grãos (δ_r), e o expoente de freqüência (η). O tempo de relaxação fractal (Tf), e o expoente de frequencia (η) estão relacionados à geometria fractal das interfaces rugosas entre os minerais condutivos (grãos metálicos e/ou partículas de argila bloqueando os canais dos poros) e o eletrólito. O tempo de relaxação (T) é um resultado da relaxação em baixa freqüência das duplas camadas elétricas formadas nas interfaces eletrólito-cristais, enquanto (T₀) é o tempo de relaxação macroscópico da amostra como um todo. O fator de resistividade dos grãos (δ_r) relaciona a resistividade dos grãos condutivos com o valor de resistividade CC da rocha. A resistividade CC da rocha (p₀), e δ_r estão relacionados à porosidade, à condutividade do eletrólito e às relações mineralógicas entre a matriz e os grãos condutivos. O modelo foi testado sobre um intervalo largo de freqüências contra dados experimentais de amplitude e fase da resistividade bem como para dados de constante dielétrica complexa. Os dados utilizados neste trabalho foram obtidos a partir da digitalização de dados experimentais publicados, obtidos por diversos autores e englobando amostras de rochas sedimentares, ígneas e metam6rficas. É mostrado neste trabalho que os parâmetros deste modelo permitem identificar diferenças texturais e mineralógicas nas rochas. Bote modelo foi introduzido, primeiramente, como propriedade intrínseca de um semiespaço homogêneo sendo demonstrado, neste trabalho, que a resposta observada em superfície reflete as propriedades intrínsecas do meio polarizável, sendo o acoplamento eletromagnético desprezível em freqüências menores que 104 Hz. Em seguida, o meio polarizável foi embebido em um pacote de N camadas sendo demonstrado que os parâmetros fractais do meio polarizável podem ser obtidos do levantamento em superfície para diferentes espessuras dessa camada. Isto justifica a utilização pura e simples de modelos de polarização desenvolvidos para amostras em laboratório para ajustar dados de campo, o que vem sendo feito sem uma justificativa bem fundamentada. Estes resultados demonstram a importância para a prospecção geolétrica do modelo proposto nesta tese.

Simulação de perfis nucleares de poço em formações complexas

A identificação e descrição dos caracteres litológicos de uma formação são indispensáveis à avaliação de formações complexas. Com este objetivo, tem sido sistematicamente usada a combinação de ferramentas nucleares em poços não-revestidos. Os perfis resultantes podem ser considerados como a interação entre duas fases distintas: • Fase de transporte da radiação desde a fonte até um ou mais detectores, através da formação. • Fase de detecção, que consiste na coleção da radiação, sua transformação em pulsos de corrente e, finalmente, na distribuição espectral destes pulsos. Visto que a presença do detector não afeta fortemente o resultado do transporte da radiação, cada fase pode ser simulada independentemente uma da outra, o que permite introduzir um novo tipo de modelamento que desacopla as duas fases. Neste trabalho, a resposta final é simulada combinando soluções numéricas do transporte com uma biblioteca de funções resposta do detector, para diferentes energias incidentes e para cada arranjo específico de fontes e detectores. O transporte da radiação é calculado através do algoritmo de elementos finitos (FEM), na forma de fluxo escalar 2½-D, proveniente da solução numérica da aproximação de difusão para multigrupos da equação de transporte de Boltzmann, no espaço de fase, dita aproximação P₁, onde a variável direção é expandida em termos dos polinômios ortogonais de Legendre. Isto determina a redução da dimensionalidade do problema, tornando-o mais compatível com o algoritmo FEM, onde o fluxo dependa exclusivamente da variável espacial e das propriedades físicas da formação. A função resposta do detector NaI(Tl) é obtida independentemente pelo método Monte Carlo (MC) em que a reconstrução da vida de uma partícula dentro do cristal cintilador é feita simulando, interação por interação, a posição, direção e energia das diferentes partículas, com a ajuda de números aleatórios aos quais estão associados leis de probabilidades adequadas. Os possíveis tipos de interação (Rayleigh, Efeito fotoelétrico, Compton e Produção de pares) são determinados similarmente. Completa-se a simulação quando as funções resposta do detector são convolvidas com o fluxo escalar, produzindo como resposta final, o espectro de altura de pulso do sistema modelado. Neste espectro serão selecionados conjuntos de canais denominados janelas de detecção. As taxas de contagens em cada janela apresentam dependências diferenciadas sobre a densidade eletrônica e a fitologia. Isto permite utilizar a combinação dessas janelas na determinação da densidade e do fator de absorção fotoelétrico das formações. De acordo com a metodologia desenvolvida, os perfis, tanto em modelos de camadas espessas quanto finas, puderam ser simulados. O desempenho do método foi testado em formações complexas, principalmente naquelas em que a presença de minerais de argila, feldspato e mica, produziram efeitos consideráveis capazes de perturbar a resposta final das ferramentas. Os resultados mostraram que as formações com densidade entre 1.8 e 4.0 g/cm³ e fatores de absorção fotoelétrico no intervalo de 1.5 a 5 barns/e^-, tiveram seus caracteres físicos e litológicos perfeitamente identificados. As concentrações de Potássio, Urânio e Tório, puderam ser obtidas com a introdução de um novo sistema de calibração, capaz de corrigir os efeitos devidos à influência de altas variâncias e de correlações negativas, observadas principalmente no cálculo das concentrações em massa de Urânio e Potássio. Na simulação da resposta da sonda CNL, utilizando o algoritmo de regressão polinomial de Tittle, foi verificado que, devido à resolução vertical limitada por ela apresentada, as camadas com espessuras inferiores ao espaçamento fonte - detector mais distante tiveram os valores de porosidade aparente medidos erroneamente. Isto deve-se ao fato do algoritmo de Tittle aplicar-se exclusivamente a camadas espessas. Em virtude desse erro, foi desenvolvido um método que leva em conta um fator de contribuição determinado pela área relativa de cada camada dentro da zona de máxima informação. Assim, a porosidade de cada ponto em subsuperfície pôde ser determinada convolvendo estes fatores com os índices de porosidade locais, porém supondo cada camada suficientemente espessa a fim de adequar-se ao algoritmo de Tittle. Por fim, as limitações adicionais impostas pela presença de minerais perturbadores, foram resolvidas supondo a formação como que composta por um mineral base totalmente saturada com água, sendo os componentes restantes considerados perturbações sobre este caso base. Estes resultados permitem calcular perfis sintéticos de poço, que poderão ser utilizados em esquemas de inversão com o objetivo de obter uma avaliação quantitativa mais detalhada de formações complexas.

Tomografia eletromagnética para caracterização de reservatórios de hidrocarbonetos

Na produção de petróleo é importante o monitoramento dos parâmetros do reservatório (permeabilidade, porosidade, saturação, pressão, etc) para o seu posterior gerenciamento. A variação dos parâmetros dinâmicos do reservatório induz variações na dinâmica do fluxo no reservatório, como por exemplo, perdas na pressão, dificultando o processo de extração do óleo. A injeção de fluidos aumenta a energia interna do reservatório e incrementa a pressão, estimulando o movimento do óleo em direção aos poços de extração. A tomografia eletromagnética poço-a-poço pode se tomar em uma técnica bastante eficaz no monitoramento dos processos de injeção, considerando-se o fato de ser altamente detectável a percolação de fluidos condutivos através das rochas. Esta tese apresenta o resultado de um algoritmo de tomografia eletromagnética bastante eficaz aplicado a dados sintéticos. O esquema de imageamento assume uma simetria cilíndrica em torno de uma fonte constituída por um dipolo magnético. Durante o processo de imageamento foram usados 21 transmissores e 21 receptores distribuídos em dois poços distanciados de 100 metros. O problema direto foi resolvido pelo método dos elementos finitos aplicado à equação de Helmhotz do campo elétrico secundário. O algoritmo resultante é válido para qualquer situação, não estando sujeito às restrições impostas aos algoritmos baseados nas aproximações de Born e Rytov. Por isso, pode ser aplicado eficientemente em qualquer situação, como em meios com contrastes de condutividade elétrica variando de 2 a 100, freqüências de 0.1 a 1000.0 kHz e heterogeneidades de qualquer dimensão. O problema inverso foi resolvido por intermédio do algoritmo de Marquardt estabilizado. A solução é obtida iterativamente. Os dados invertidos, com ruído Gaussiano aditivo, são as componentes em fase e em quadratura do campo magnético vertical. Sem o uso de vínculos o problema é totalmente instável, resultando em imagens completamente borradas. Duas categorias de vínculos foram usadas: vínculos relativos, do tipo suavidade, e vínculos absolutos. Os resultados obtidos mostram a eficiência desses dois tipos de vínculos através de imagens nítidas de alta resolução. Os tomogramas mostram que a resolução é melhor na direção vertical do que na horizontal e que é também função da freqüência. A posição e a atitude da heterogeneidade é bem recuperada. Ficou também demonstrado que a baixa resolução horizontal pode ser atenuada ou até mesmo eliminada por intermédio dos vínculos.

Aplicação de análise de séries temporais em perfis de poço

Esta tese descreve a aplicação de análise de séries temporais em perfis de poço. Através desta técnica é possível avaliar-se a repetição e a resolução vertical de perfis, e determinar-se o intervalo de amostragem e a velocidade de perfilagem ideais para diferentes perfis. A comparação entre três poços é também feita, baseada num mesmo tipo de perfil. Para tanto, na seqüência utilizada, procurou-se manter num mesmo domínio os dados cuja quantidade total de amostras (N) por perfil não ultrapassou 2048. Desses dados, foram inicialmente retirados o valor médio das amostras e o alinhamento polinomial algébrico eventualmente nelas embutido. Em seguida, foi efetuada a aplicação do ponderador cossenoidal, do filtro passa-alta, da janela Hanning, do cálculo da função coerência, do espectro de fase, da razão sinal-ruído e dos espectros de potência do sinal e do ruído, nesta ordem. Para a função coerência, fez-se necessário o cálculo dos níveis de confiança de 50%, 90% e 95%. O cálculo do primeiro nível teve por base a necessidade de se determinar a resolução vertical de alguns perfis, e dos demais, a fim de que fosse obtida uma informação referente à localização daqueles níveis para a coerência calculada. Em relação ao espectro de fase, seu cálculo surgiu da necessidade de se obter uma informação adicional a respeito dos perfis manipulados, ou seja, o conhecimento da ocorrência ou não de deslocamento relativo de profundidade entre a seção principal e a seção repetida. A razão sinal-ruído foi calculada no sentido de possibilitar a comparação, como elemento avaliador dos diversos tipos de perfis, com a coerência e o cálculo dos espectros de potência. Os espectros de potência do sinal e do ruído foram calculados para se ter mais um parâmetro de avaliação da seção repetida, já que em tese, os espectros de potência do sinal e do ruído da seção repetida devem ser iguais aos respectivos espectros da seção principal. Os dados utilizados na aplicação da metodologia proposta foram fornecidos pela PETROBRÁS e oriundos de quatro poços da Bacia Potiguar emersa. Por questões de sigilo empresarial, os poços foram identificados como poços A, B, C e D. A avaliação da repetição entre diferentes tipos de perfis indica que, para o poço A, o perfil micro-esférico (MSFL) tem repetição melhor do que o perfil de porosidade neutrônica (CNL), o qual tem, por sua vez, repetição melhor do que o perfil de raios gama normal (GR). Para os perfis do poço D, uma diminuição da velocidade de perfilagem de 550 m/h para 275 m/h é vantajosa apenas para o perfil de porosidade neutrônica. Já a velocidade de perfilagem de 920 m/h, utilizada, na obtenção dos perfis do poço C, é totalmente inadequada para os perfis de resistividade (MSFL, ILD e ILM). A diminuição do intervalo de amostragem de 0,20 m para 0,0508 m, nos perfis de raios gama e de porosidade neutrônica, e 0,0254 m para o perfil de densidade, apresenta bons resultados quando aplicada no poço D. O cálculo da resolução vertical indica, para o perfil de porosidade neutrônica, uma superioridade qualitativa em relação ao perfil de raios gama normal, ambos pertencentes ao poço A. Para o poço C, o perfil micro-esférico apresenta uma resolução vertical na mesma ordem de grandeza da resolução do perfil de raios gama do poço B, o que evidencia ainda mais a inconveniência da velocidade de perfilagem utilizada no poço C. Já para o poço D, o cálculo da resolução vertical indica uma superioridade qualitativa do perfil de densidade de alta resolução em relação ao perfil de raios gama de alta resolução. A comparação entre os poços A, B e D, levada a efeito através dos respectivos perfis de porosidade neutrônica normais, comprova que a presença de ruído aleatório, em geral, está diretamente ligada à porosidade da formação - uma maior porosidade indica uma presença maior de ruído e, por conseguinte, uma queda qualitativa no perfil obtido. A análise do espectro de fase de cada perfil indica um deslocamento em profundidade, existente entre as seções principal e repetida de todos os perfis do poço C. E isto pôde ser confirmado com a posterior superposição das seções.