Controle de vibração em reatores elétricos através de absorvedores dinâmicos viscoelásticos

Reatores elétricos utilizados em subestações de energia têm funções específicas na operacionalização e manutenção da funcionalidade destas instalações, sendo equipamentos aplicados ao controle de reativos, de correntes de curto circuito e surtos de manobra em sistemas de transmissão de energia elétrica. Detectou-se que alguns reatores em funcionamento apresentaram elevados níveis de vibração estrutural. Dois destes reatores operam na subestação da cidade de Rurópolis, situada a oeste do estado do Pará. Outro problema decorrente do elevado nível de vibração destas estruturas é a geração de pontos quentes no núcleo eletromagnético do reator, podendo gerar bolhas de gases altamente inflamáveis. A existência destes problemas operacionais pode ocasionar ao Sistema de Transmissão do Norte do Brasil interrupções na transmissão de energia elétrica. Isto geraria na região Norte e em outras interligadas (Nordeste, Sudeste e Centro-Oeste) elevados prejuízos sociais e financeiros. A aplicação de ferramentas para a minimização do problema vibro-acústico (redução dos níveis vibração e, por conseguinte, ruído) pode evitar danos estruturais e operacionais que resultem na sua indisponibilidade operativa reduzindo, conseqüentemente, perdas e a necessidade do desenvolvimento de equipamentos de custo elevado com baixos índices de vibração. No intuito de se reduzir tais níveis, propõe-se à aplicação de absorvedores dinâmicos na estrutura externa dos reatores. Uma vez dimensionados em quantidade, posição e parâmetros definidores (massa, rigidez e amortecimento) espera-se criar um sistema, que sintonizado com a freqüência da excitação de natureza eletromagnética, absorverá parte da energia vibratória do reator, minimizando-se, os elevados índices detectados e, por conseqüência, todos os problemas subseqüentes. São indicados compostos viscoelásticos nos absorvedores dinâmicos (com propriedades dinâmicas dependentes tanto da temperatura como da freqüência de trabalho), em virtude de estes materiais aumentarem o campo de ação dos absorvedores, adequando-se à elevada densidade modal verificada nas chaparias metálicas do reator.

Determinação estatística dos contatos fluidos em perfis de poço

No interior dos reservatórios petrolíferos, o efeito da gravidade produz naturalmente a segregação dos fluidos presentes. Em razão basicamente das forças capilares não existe uma superfície definida de separação entre cada dois fluidos. No entanto, é comum na engenharia de petróleo admitir-se uma superfície convencional de separação dos fluidos. À essa superfície dá-se o nome de interface fluida ou contato fluido. A obtenção da localização dos contatos fluidos no interior dos reservatórios, como os contatos óleo-água, contatos gás-água e os contatos gás-óleo, sem a interpretação completa dos perfis geofísicos ou a realização dos procedimentos diretos dos testes de formação é um problema de difícil solução para a indústria do petróleo. O posicionamento em profundidade dos contatos pode ser utilizado para a definição dos procedimentos de completação do poço e para o posicionamento do canhoneio nos horizontes produtores dos poços verticais. Temos também a aplicação do mapeamento deste parâmetro, para os projetos e controle da perfuração no caso dos poços direcionais e horizontais, que utilizam tal informação para o posicionamento do eixo do poço no interior da zona de hidrocarboneto no reservatório, a fim de se evitar, principalmente, a produção de água. Este trabalho apresenta uma metodologia para a identificação e posicionamento dos contatos fluidos, através da aplicação da técnica da estatística multivariada, chamada de análise discriminante. Nesta aplicação, faz-se uso direto dos perfis de resistividade (R_T), de raios gama natural (RG), de densidade (ρ_b), de porosidade neutrônica (Φ_N), de caliper (Cal) e do perfil de volume de argila (Vcla), como dados de entrada no programa de análise discriminante. Para o caso de depósitos clásticos, com seqüências arenito-folhelho, a análise discriminante fornece ainda a indicação da litologia e as espessuras aparentes dos reservatórios. As curvas de perfis geofísicos aplicadas para a avaliação desta metodologia são provenientes de poços perfurados na bacia do lago Maracaibo, na Venezuela, onde, segundo estudos geológicos, existe a ocorrência de seqüências estratigráficas, com camadas de folhelhos e arenitos, sendo também confirmada a presença de hidrocarbonetos nos poços utilizados neste trabalho.

Inclusões fluidas crepitadas, fluidos hipersalinos e aquo-carbônicos em quartzo associado a rochas micáceas no Granito Xinguara - Terreno Granito-Greenstone de Rio Maria, PA

As rochas micáceas encontradas no Granito Xinguara, terreno Granito-Greenstone de Rio Maria, Pará, são compostas por muscovita e clorita com níveis de quartzo intercalados, que formam uma xistosidade bem desenvolvida. Essa xistosidade é cortada por veios de quartzo. Ambas as gerações de quartzo apresentam os mesmos tipos de inclusões fluidas em halos ou trilhas secundárias de composições variadas entre aquosas, aquo-carbônicas e saturadas em torno de grandes inclusões primárias crepitadas ou em trilhas transgranulares secundárias. A grande variação de temperaturas de homogeneização, a alta salinidade, as evidências de estrangulamento e a existência das inclusões crepitadas permitem supor forte influência de alterações pós-formacionais e reequilíbrio relacionados à intrusão do granito. Essas rochas foliadas são, portanto, enclaves metassedimentares afetados por fluidos graníticos hipersalinos aquo-carbônicos.

Evolução dos fluidos mineralizantes e modelo genético dos veios de quartzo auríferos em zona de cisalhamento do depósito Ouro Roxo, Província Tapajós, Jacareacanga (PA), Brasil

O depósito aurífero Ouro Roxo, localizado no município de Jacareacanga, Província Aurífera do Tapajós, sudoeste do Pará, formou-se em um sistema hidrotermal que gerou veios de quartzo sulfetados, em zona de cisalhamento N-S, dúctil-rúptil, oblíqua, denominada Ouro Roxo-Canta Galo, cortando granitoides calcioalcalinos da Suíte Intrusiva Tropas, de idade paleoproterozoica e hospedeira da mineralização, em rochas localmente milonitizadas. Três tipos de fluidos foram caracterizados como geradores do depósito: 1) fluido aquoso H₂O-NaCl-MgCl₂-FeCl₂ de salinidade baixa a moderada, com temperatura de homogeneização total (Th) = 180-280°C; 2) salmoura H₂O-NaCl-CaCl₂ com Th = 270-400°C, provavelmente portadoras de Cu e Bi, relacionadas geneticamente a um evento magmático contemporâneo ao cisalhamento que sofreu diluição pela mistura com água meteórica, baixando sua salinidade e temperatura (Th = 120-380°C); 3) fluido aquocarbônico de média salinidade, com Th = 230-430°C, que foi interpretado como o fluido mineralizante mais primitivo, provavelmente aurífero, relacionado com o cisalhamento. As condições de temperatura e pressão (T-P) de formação do minério, estimadas conjuntamente pelo geotermômetro da clorita e as isócoras das inclusões fluidas, situam-se entre 315 e 388°C e 2 a 4,1kb. Dois mecanismos simultâneos provocaram a deposição do minério em sítios de transtensão da zona de cisalhamento: 1) mistura de fluido aquocarbônico com salmoura magmática com aumento de fO2 e redução de pH; 2) interação entre os fluidos e os feldspatos e minerais ferromagnesianos do granitoide hospedeiro, com reações de hidrólise e sulfetação, provocaram redução de fO2 e fS2, com precipitação de sulfetos de Fe juntamente com ouro. O ambiente orogênico, o estilo filoneano do depósito, o controle estrutural pela zona de cisalhamento, a alteração hidrotermal (propilítica + fílica + carbonatação), a associação metálica (Au + Cu + Bi), o fluido mineralizante aquocarbônico associado com salmoura magmática na deposição do minério são compatíveis com um modelo orogênico com participação magmática para a gênese do depósito Ouro Roxo.

Estudo da influência da granulometria no comportamento reológico de polpa de bauxita e no fator de atrito

O presente trabalho apresenta estudos de caracterização reológica e determinações de fatores de atrito em três polpas de bauxita, originárias da Mina de Miltônia, Paragominas – Pará, resultantes de condições operacionais com graus diferentes de moagem. Deste modo, a diferenciação básica entre as mesmas reside na distribuição granulométrica. O objetivo é fazer uma proposta de revisão da especificação granulométrica do produto (polpa de bauxita para o mineroduto), permitindo uma comparação com as polpas do projeto (planta piloto) e da operação atual (usina de beneficiamento) quanto a resposta a alterações, principalmente, da viscosidade e fatores de atrito. A polpa proposta incorpora um percentual maior de finos (< 10 microns) em relação as demais polpas. Os benefícios possibilitarão ganhos de recuperação com a incorporação desta parcela de finos, a qual atualmente é descartada para a bacia de rejeitos. O material foi caracterizado por análises granulométricas da série Tyler, físicas e físico-químicas, microscopia eletrônica de varredura (MEV) e espectrometria de raios-x dispersiva em energia (EDS). Utilizando-se de um viscosímetro de cilindros coaxiais, tipo Searle, abordo use os aspectos reológicos destas polpas, em função da sua distribuição granulométrica e da variação da concentração de sólidos. As propriedades reológicas avaliadas foram a viscosidade, a tensão de escoamento e também os comportamentos que mostram a dependência da viscosidade com o tempo de aplicação de uma taxa de cisalhamento constante, verificando-se por meio de curva de histerese, a existência ou não da tixotropia e da reopexia. Os resultados experimentais mostraram que para as polpas ensaiadas, a granulometria e a concentração de sólidos são parâmetros relevantes na determinação da viscosidade. Os testes com estas três polpas de bauxita, de granulometrias diferentes, foram realizados a concentrações pré-determinadas, de modo a determinar o modelo reológico cujos parâmetros apresentaram os melhores coeficientes de correlação (R²), sendo que os melhores ajustes encontrados foram segundo o modelo de Herschel-Bulkley. Os fatores de atrito foram obtidos por formulações que envolvem uma relação logarítmica, utilizando-se da Lei da Parede como método de obtenção das expressões aplicáveis a estes fluidos não- Newtonianos viscoplásticos, sendo comparados com os calculados pelas correlações de Dodge e Metzner e verificou-se que de uma maneira geral há a concordância entre o modelo estudado nas três polpas. A análise dos fatores de atrito na concentração de 50% de sólidos, concentração de bombeio por mineroduto, mostra que o comportamento das polpas estudadas é semelhante, com os melhores ajustes para a faixa de 10.000 a 100.000 Reynolds (faixa média de trabalho para o escoamento turbulento de polpas de minério, homogêneas e heterogêneas).

Avaliação do teste de aglutinação com partículas de látex sensibilizadas com exoantígeno bruto de Paracoccidioides brasiliensis no sorodiagnóstico da Paracoccidioidomicose

Paracoccidioidomicose (PCM), uma micose sistêmica granulomatosa causada pelo fungo dimórfico Paracoccidioides brasiliensis, é endêmica na América do Sul. Conídios provavelmente agem como propágulos infectantes e são inalados para os pulmões, onde ocorre a transformação à forma leveduriforme patogênica. Duas principais formas clínicas são consideradas: a forma aguda ou subaguda (tipo juvenil) e a forma crônica (tipo adulto). O diagnóstico definitivo da PCM inclui a observação direta da levedura multibrotante característica em fluidos biológicos e secções teciduais ou isolamento do fungo de materiais clínicos. Na PCM, testes sorológicos, além de auxílio diagnóstico, têm a função de acompanhamento durante e pós-tratamento. Portanto, a técnica utilizada precisa aliar sensibilidade e especificidade, para que o valor preditivo seja máximo e reprodutível. O propósito deste estudo foi avaliar um teste de aglutinação com látex (LA) para detectar anticorpos anti-P.brasiliensis contra antígeno bruto do fungo. Cinqüenta e uma (51) amostras de soro de pacientes com PCM foram testadas. Positividade foi observada em 84,31% (43/51), cujos padrões de aglutinação variaram de 1+ a 4+. Reatividade dessas amostras foi verificada em títulos variando entre 1:2 e 1:64. Reatividade cruzada foi observada com outras doenças fúngicas (aspergilose e histoplasmose), e com doenças não fúngicas. Amostras de soro humano normal não foram reativas. A sensibilidade, especificidade e valores preditivos, positivo e negativo, produzidos pelo teste LA foram 84,31%, 81,05%, 70,49% e 90,59%, respectivamente. Em conclusão, estes resultados mostram que o teste LA é instrumento útil no sorodiagnóstico da PCM, além de vantagens como baixo custo e rápida execução, a despeito de outros testes, tais como ID e Western blotting.

Interação CFD-DEM em fluidização: aplicação para o setor mineral

A fluidização de partículas é amplamente utilizada na indústria, principalmente devido às altas taxas de transferência de calor e massa entre as fases. O acoplamento entre a Dinâmica dos Fluidos Computacional (CFD – Computational Fluid Dynamics) e o Método dos Elementos Discretos (DEM – Discrete Element Method) tem se tornado atrativo para a simulação de fluidização, já que nesse caso o movimento das partículas é analisado de forma mais direta do que em outros tipos de abordagens. O grande problema do acoplamento CFD-DEM é a alta exigência computacional para rastrear todas as partículas do sistema, o que leva ao uso de estratégias de redução do tempo de simulação que em caso de utilização incorreta podem comprometer os resultados. O presente trabalho trata da aplicação do acoplamento CFD-DEM na análise de fluidização de alumina, que é um problema importante para o setor mineral. Foram analisados diversos parâmetros capazes de influenciar os resultados e o tempo de simulação como os passos de tempo, os modelos de arrasto, a distribuição granulométrica das partículas, a constante de rigidez, a utilização de partículas representativas com tamanho maior que o das partículas reais, etc. O modelo de força de interação DEM utilizado foi o modelo de mola e amortecedor lineares (LSD – Linear Spring Dashpot). Todas as simulações foram realizadas com o software ANSYS FLUENT 14.5 e os resultados obtidos foram comparados com dados experimentais e da literatura. Tais resultados permitiram comprovar a capacidade do modelo linear LSD em predizer o comportamento global de leitos de alumina e reduzir o tempo de simulação, desde que os parâmetros do modelo sejam definidos de forma adequada.

Modos quasinormais e pólos de Regge para os buracos acústicos canônicos

Usando o formalismo relativístico no estudo da propagação de perturbações lineares em fluidos ideais, obtêm-se fortes analogias com os resultados encontrados na Teoria da Relatividade Geral. Neste contexto, de acordo com Unruh [W. Unruh, Phys. Rev. Letters 46, 1351 (1981)], é possível simular um espaço-tempo dotado de uma métrica efetiva em um fluído ideal barotrópico, irrotacional e perturbado por ondas acústicas. Esse espaço-tempo efetivo é chamado de espaço-tempo acústico e satisfaz as propriedades geométricas e cinemáticas de um espaço-tempo curvo. Neste trabalho estudamos os modos quasinormais (QNs) e os pólos de Regge (PRs) para um espaço-tempo acústico conhecido como buraco acústico canônico (BAC). No nosso estudo, usamos o método de expansão assintótica proposto por Dolan e Ottewill [S. R. Dolan e A. C. Ottewill, Class. Quantum Gravity 26, 225003 (2009)] para calcularmos, em termos arbitrários do número de overtone n, as frequências QNs e os momentos angulares para os PRs, bem como suas respectivas funções de onda. As frequências e as funções de onda dos modos QNs são expandidas em termos de potências inversas de L = l + 1/2 , onde l é o momento angular, enquanto que os momentos angulares e funções de onda dos PRs são expandidos em termos do inverso das frequências de oscilação do buraco acústico canônico. Comparamos os nossos resultados com os já existentes na literatura, que usam a aproximação de Wentzel-Kramers-Brillouin (WKB) como método de determinação dos modos QNs e dos PRs, e obtemos uma excelente concordância dentro do limite da aproximação eikonal (l ≥ 2 e l > n).

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 tratamento estatístico multivariante em dados de perfis de poços da Bacia de Sergipe-Alagoas

Utilizando-se dados do campo de Camorim (Bacia de Sergipe-Alagoas), foi testado e aplicado um conjunto de técnicas estatísticas multivariantes (análises de agrupamentos, de componentes principais e discriminante) no intuito de identificar as fácies, previamente definidas em poços testemunhados, por meio dos perfis, viabilizando o reconhecimento das mesmas nos demais poços não testemunhados da área. A segunda etapa do processo de determinação das fácies consistiu no emprego de métodos auxiliares (análises composicional e de seqüência de fácies), que, combinados com as técnicas multivariantes, propiciaram melhores resultados na calibração rocha-perfil. A determinação das fácies, uma vez estabelecida, possibilitou o refinamento do processo de avaliação de formações ao viabilizar o exame de cada fácies-reservatório, isoladamente. Assim, esse procedimento tornou possível a escolha, para cada litologia, dos parâmetros utilizados na interpretação dos perfis ao mesmo tempo em que permitiu a totalização em separado dos valores de espessura, porosidade e saturação dos fluidos, bem como a adoção de diferentes valores de corte (cut-offs) para cada grupo considerado. Outras aplicações incluíram a melhoria na estimativa da porosidade e da permeabilidade, a adaptação de algoritmos para o cálculo preliminar de porosidade, a confecção de mapas de fácies e a geração automática de seções estratigráficas. Finalmente, foram destacadas a perspectiva de integração desse estudo com sistemas estatísticos de descrição de reservatórios, outras técnicas de determinação de fácies em desenvolvimento e a retomada da utilização de métodos estatísticos multivariantes em dados de perfis, como ferramenta de exploração.

Detecção do contato entre camadas através do perfil de indução

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.

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.

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

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.

Um estudo da estabilidade mecânica de poços a partir de perfis geofísicos em formações pouco consolidadas

A análise da estabilidade mecânica de um poço pode ser feita a partir do cálculo de parâmetros elásticos da formação utilizando a densidade do meio e as velocidades de propagação das ondas compressional e cisalhante na formação rochosa, os quais podem ser obtidos de perfis geofísicos do poço. Em formações sedimentares pouco consolidadas as ferramentas de perfilagem sônica convencionais (monopolares) não conseguem registrar com acuidade a velocidade da onda cisalhante pois a primeira chegada dessa onda é camuflada pela chegada de outras ondas que podem ser mais rápidas que a onda cisalhante num poço perfurado neste tipo de formação. Medidas das velocidades sônicas são feitas em laboratório em amostras da formação, sob condições semelhantes às condições in situ, servindo como ajuste das velocidades registradas no poço pela ferramenta de perfilagem sônica. Para a análise de estabilidade da formação, perfis auxiliares são necessários como o perfil de porosidade, saturação de fluidos e perfis de composição mineralógica da formação rochosa. Exige-se ainda dados de testes de avaliação da formação e de condições do reservatório, mas que são comuns em poços de petróleo, como o teste de formação e os testes de pressurização do poço, tais como o teste de micro-fraturamento hidráulico ou o teste de absorção. A avaliação das tensões principais efetivas que atuam distante do poço, e que não são afetadas pela sua presença, é feita através da associação de um modelo de deformação elástica apropriado e o resultado do teste de pressurização disponível para o poço em estudo. Utilizando resultados clássicos da teoria da elasticidade geral pode-se calcular o campo de tensões modificado na vizinhança da parede do poço devido ao efeito da própria presença do poço ali perfurado e da diferença de pressão existente entre o interior do poço e a formação rochosa. A determinação das propriedades mecânicas da formação a partir das velocidades sônicas e a avaliação do campo de tensões assumindo um modelo elástico de deformação, supõem o meio rochoso no qual as ondas se propagam como elástico, homogêneo e isotrópico. Esta suposição representa a principal aproximação assumida pela metodologia descrita neste trabalho. De posse das propriedades mecânicas da formação e do campo de tensões que age na vizinhança do poço resta definir o critério segundo o qual a rocha sofre instabilidade mecânica quando submetida aquele campo de tensões. Isto permite determinar se, nas condições avaliadas do poço e da formação, haverá quebra da parede do poço por excesso de tensão e, se houver, qual a sua extensão. Assim o problema é como analisar o comportamento mecânico de um poço em uma formação pouco consolidada a partir de perfis geofísicos os quais podem ter problemas no registro das propriedades físicas do meio em formações deste tipo. A metodologia proposta é aplicada a dois intervalos de profundidade pertencentes a dois poços onde arenitos e folhelhos se intercalam e nos quais todos os dados necessários estão disponíveis. Os resultados obtidos mostram que, exceto quando outros mecanismos de quebra da parede do poço agem na formação, a metodologia proposta consegue com sucesso detectar zonas de ocorrência de instabilidade mecânica do poço provocadas por um campo de tensões que excede a resistência mecânica da formação.

Estudo da resposta de perfilagem de indução de camadas finas com diferentes arranjos de bobinas: modelamento analógico

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.