Controle automático de adição de AlF3 em cubas de redução de alumínio, usando lógica fuzzy

Este trabalho descreve um sistema de análise de dados com a finalidade de gerar um sistema de controle utilizando técnica inteligente para adição de fluoreto de alumínio (AlF3) em fornos de redução de alumínio. O projeto baseia-se nos conceitos de lógica fuzzy, nos quais o conhecimento acumulado pelo especialista do processo é traduzido de maneira qualitativa em um conjunto de regras linguísticas do tipo SE ENTÃO. A utilização desta técnica inteligente para o controle de adição de fluoreto busca representar explicitamente um conhecimento qualitativo, detido pelos operadores de cubas eletrolíticas. Devido o sistema convencional não contemplar as variações dos fenômenos que envolvem a dinâmica do processo, um controlador fuzzy foi implmentado no sistema real para tomadas de decisões, utilizando o modelo mínimo de Mandani. Baseado neste modelo, as variáveis de processo para a entrada do sistema fuzzy, tais como temperatura de banho e percentual de fluoreto foram manipuladas para estimar a tendência de subida e descida, respectivamente, através do método mínimos quadrados(MMQ). O controlador fuzzy é aplicado para calcular a quantidade de fluoreto de alumínio (AlF3) a ser adicionado na cuba eletrolítica de forma automática sem a necessidade da intervenção do especialista do processo. A motivação para o uso de um sistema de controle fuzzy se deve ao fato de não se ter disponível um modelo dinâmico do processo de adição do fluoreto na cuba eletrolítica. Esta falta de modelagem se deve ao fato de grande complexidade dos fenômenos envolvidos em uma cuba que são processos termodinâmicos e eletromagnéticos acoplados.

Determinação de puberdade em cutias (Rodentia:Dasyproctidae) criadas em cativeiro, através da quantificação das células espermatogênicas

A cutia (Dasyprocta spp) é um roedor de médio porte, que apresenta o hábito de enterrar parte de seus alimentos, principalmente sementes, é considerado um dispersor em potencial, pois contribui com o reflorestamento natural. Por ser uma espécie muito procurada para consumo na área rural da Amazônia, o manejo em cativeiro e o conhecimento dos seus aspectos reprodutivos são formas alternativas para a criação de programas de produção e preservação dessa espécie com potencial econômico. O objetivo do presente trabalho, foi determinar o período que ocorre a puberdade, caracterizar os estádios do Ciclo do Epitélio Seminífero (CES), determinar a freqüência relativa dos estádios, calcular o rendimento geral da espermatogênese e o índice das células de sertoli. Foram utilizados 7 grupos com idade variando de 4 a 17 meses, os animais foram divididos em G1 (4 e 5 meses, n=4), G2 (6 e 7 meses, n=4), G3 (8 e 9 meses, n=4), G4 (10 e 11 meses, n=3), G5 (12 e 13 meses, n=4), G6 (14 e 15 meses, n=3) e G7 (16 e 17 meses, n=2). Ao atingir a idade programada os animais foram castrados, sob anestesia, e as amostras testiculares após biometria foram fixadas em ALFAC por 24 horas, submetidas ao processamento histológico de rotina, foram realizados cortes de 5 mm de espessura e os tecidos obtidos corados com HE. As fases de desenvolvimento reprodutivo, incluindo a puberdade, foram determinadas através da quantificação das células espermatogênicas de 10 túbulos seminíferos/animal com o contorno circular, os quais se encontravam no estádio I do CES, previamente caracterizado pelo método da morfologia tubular, que apresentou oito estádios do CES, e de 20 túbulos seminíferos/animal que não apresentavam espermatogênese completa, o número real das células foi obtido através da correção dos números brutos pelo diâmetro nuclear/nucleolar médio e espessura do corte histológico. Para determinar a freqüência relativa, 100 túbulos seminíferos de cada animal que já havia atingido a puberdade foram analisados. Os pesos corporal e testicular apresentaram correlações significativas com a idade e entre si. O peso corporal e biometria testicular aumentaram significativamente (P<0.05) até a maturidade sexual, onde na puberdade o PC e PT foram (1.903 + 0.55 Kg; 1.8 + 1.4 g, respectivamente) e na fase adulta foram (2.825 + 0.11Kg; 6.1 + 0.37g). Os grupos analisados foram classificados como G1 impúbere (13.21 + 1.6 mm); G2 pré-puberdade (7.09 + 1.7 mm); G3 puberdade (93.82 + 58.7 mm); G4 Pós-puberdade 1 (164.03 + 10.03 mm); G5 e G6 Pós-puberdade 2 de fase longa (173.8 + 10.5 mm; 185.9 + 1.5 mm) e G7 Adulto (238.8 + 72.6 mm). As células de sertoli, diminuiram significativamente (P<0.05) da fase impúbere (21.7 + 3.5 mm) até a puberdade (9.7 + 4.2 mm) onde iniciaram sua estabilização até a fase adulta (9.03 + 0.01mm). As células espermatogênicas apresentaram correlações altas e significativas com o peso testicular. As médias das seções transversais do diâmetro tubular aumentaram significativamente (P<0.05) entre os grupos analisados, G1 (109.5 + 4.6 mm); G2 (119.7 + 10 mm); G3 (174.6 + 24.1 mm); G4 (240.9 + 14.1 mm); G5 (219.6 + 9.8 mm); G6 (221.1 + 7.9 mm) e G7 (258.1 + 55.3 mm). Após a análise de 1600 túbulos seminíferos, os oito estádios caracterizados pelo método da morfologia tubular apresentaram as seguintes freqüências relativas: I (16.8 + 2.3%); II (18.8 + 3.4%); III (7.6 + 1.3%); IV (11.1 + 1.2%); V (21.2 + 4.2%); VI (10.4 + 3.4%); VII (7.9 + 2.4%); VIII (6.4 + 1.4%). A freqüência conjunta dos estádios foi pré-meiótica (43.5%), meiótica (10.7%) e pós-meiótica (45.8%). Portanto, a puberdade em cutias é alcançada no período de 8 e 9 meses, onde se observa maior produção das células espermatogênicas e estabelecimento das células de sertoli, a partir dessa fase inicia-se uma tentativa de estabilização das células espermatogênicas, o qual é acompanhado por uma tendência a estabilização do crescimento testicular.

Synthesis, X-ray crystal structure and theoretical calculations of antileishmanial neolignan analogues

A síntese e a estrutura cristalina por difração de raios-X de dois análogos de neolignanas, 2-(4-clorofenil)-1-feniletanona (20) e 2-[tio(4-clorofenil)]-1-(3,4-dimetoxifenil)propan-1-ona (12) são descritas. O composto 12 apresenta atividade intracelular contra Leishmania donovani e Leishmania amazonensis de amastigotas que causam a leishmaniose tegumentar e visceral. Além disso, a teoria do funcional de densidade (DFT) com o funcional híbrido B3LYP foi empregado para calcular um conjunto de descritores moleculares para dezenove análogos sintéticos de neolignanas com atividades antileishmaniose. Posteriormente, a análise discriminante stepwise foi realizada para investigar possíveis relações entre a estrutura molecular e atividades biológicas. Por meio dessa análise os compostos foram classificados em dois grupos ativos e inativos de acordo com seu grau de atividade biológica, e as propriedades mais importantes foram as cargas de alguns átomos, a afinidade eletrônica e o ClogP.

A study of neolignan compounds with biological activity against Paracoccidioides brasiliensis by using quantum chemical and chemometric methods

Métodos quimiométricos (estatísticos) são empregados para classificar um conjunto de compostos derivados de neolignanas com atividade biológica contra a Paracoccidioides brasiliensis. O método AM1 (Austin Model 1) foi utilizado para calcular um conjunto de descritores moleculares (propriedades) para os compostos em estudo. A seguir, os descritores foram analisados utilizando os seguintes métodos de reconhecimento de padrões: Análise de Componentes Principais (PCA), Análise Hierárquica de Agrupamentos (HCA) e o método de K-vizinhos mais próximos (KNN). Os métodos PCA e HCA mostraram-se bastante eficientes para classificação dos compostos estudados em dois grupos (ativos e inativos). Três descritores moleculares foram responsáveis pela separação entre os compostos ativos e inativos: energia do orbital molecular mais alto ocupado (EHOMO), ordem de ligação entre os átomos C1'-R7 (L14) e ordem de ligação entre os átomos C5'-R6 (L22). Como as variáveis responsáveis pela separação entre compostos ativos e inativos são descritores eletrônicos, conclui-se que efeitos eletrônicos podem desempenhar um importante papel na interação entre receptor biológico e compostos derivados de neolignanas com atividade contra a Paracoccidioides brasiliensis.

Structure-activity relationship of compounds which are anti-schistosomiasis active

O método do orbital molecular AM1 foi empregado para calcular um conjunto de descritores moleculares para vinte neolignanas sintéticas com atividade anti-esquistossomose. O método de reconhecimento de padrão (análise de componentes principais ACP, análise de conglomerados AC e análise de discriminante) foi utilizado para obter a relação entre a estrutura molecular e a atividade biológica. O conjunto de moléculas foi classificado em dois grupos de acordo com seus graus de atividade biológica. Estes resultados permitem que, projete-se racionalmente novos compostos, potenciais candidatos à síntese e à avaliação biológica.

Amplitudes e padrões de polarização de pulsos em meios anisotrópicos

Extrair informações litológicas da subsuperfície através de dados sísmicos constitui-se num grande desafio à prospecção sísmica, pois a hipótese de estratificações formadas por camadas isotrópicas se mostra insuficiente para representar o comportamento do campo elástico em levantamentos com grandes afastamentos entre fonte e receptor, geofones multicomponentes, medidas de VSP tridimensional, entre outros. Sob este panorama, a prospecção sísmica passa a considerar modelos anisotrópicos de subsuperfície para, por exemplo, caracterizar reservatórios. O objetivo deste texto é apresentar um formalismo para modelar o espalhamento de pulsos a partir de ondas planas incidentes em interfaces planas horizontais que separam meios anisotrópicos. Este espalhamento é obtido primeiramente, através da formulação explícita dos campos de deformação e tração como função das matrizes propagadoras, de polarização e de impedância do meio. Em seguidaeste formalismo é usado para a obtenção das matrizes dos coeficientes de reflexão e transmissão através de uma interface plana horizontal para posteriormente, ser generalizado para o espalhamento através de múltiplas camadas. Finalmente, inserem-se ao campo da onda incidente as amplitudes de um pulso analítico para calcular o espalhamento do pulso através de estratificações.

3D raytracing through homogeneous anisotropic media with smooth interfaces

O raio conectando dois pontos em um meio anisotrópico, homogêneo por partes e com variação lateral, é calculado utilizando-se técnicas de continuação em 3D. Se combinado com algoritmos para solução do problema de valor inicial, o método pode ser estendido para o cálculo de eventos qS₁ e qS₂. O algoritmo apresenta a mesma eficiência e robustez que implementações de técnicas de continuação em meios isotrópicos. Rotinas baseadas neste algoritmo têm várias aplicações de interesse. Primeiramente, na modelagem e inversão de parâmetros elásticos na presença de anisotropia. Em segundo lugar, as iterações de Newton-Raphson produzem atributos da frente de onda como vetor vagarosidade e a matrix hessiana do tempo de trânsito, quantidades que permitem determinar o espalhamento geométrico e aproximações de segunda ordem para o tempo de trânsito. Estes atributos permitem calcular as amplitudes ao longo do raio e investigar os efeitos da anisotropia no empilhamento CRS em modelos de velocidade simples.

Fenômeno de interferência na densidade de energia em cavidades com duas fronteiras oscilantes

O primeiro objetivo do presente trabalho é calcular a força quântica exata que atua sobre as fronteiras de uma cavidade, bem como o comportamento exato da densidade de energia numa cavidade não estática, onde ambas as fronteiras executam movimentos prescritos arbitrários. O modelo considerado é o do campo escalar não massivo em 1 + 1 dimensões, sendo que o campo obedece à condição de Dirichlet em cada uma das fronteiras. Considerando o vácuo como estado inicial do campo, nós mostramos que a densidade de energia em um dado ponto do espaço-tempo pode ser obtida através do traçado de uma sequência de linhas nulas, conectando o valor da densidade de energia nesse ponto a um certo valor conhecido da densidade de energia em um ponto das “zonas estáticas”. O segundo objetivo é mostrar que para movimentos específicos das fronteiras, particularmente para os quais ambas voltam às suas posições iniciais em instantes múltiplos do comprimento inicial da cavidade, o método exato por nós obtido permite encontrar soluções analíticas escritas como uma expansão em série na variável que controla as amplitudes de movimento das fronteiras. Os resultados analíticos por nós obtidos são aplicáveis a uma vasta classe de movimentos, a qual inclui a grande maioria dos casos ressonantes estudados na literatura. O terceiro objetivo do presente trabalho é investigar, através dos métodos de cálculos desenvolvidos aqui, o fenômeno da interferência na energie e na densidade de energia em cavidades com duas fronteiras móveis, obtendo fórmulas genéricas para os termos de interferência respectivos.

Potenciais pseudo-newtonianos e a radiação escalar emitida por uma fonte girando ao redor de um objeto estelar

Neste trabalho usamos os potenciais pseudo-newtonianos propostos por Paczynski e Wiita, Nowak e Wagoner e Artemova et al. para calcular a radiação escalar emitida por uma fonte em movimento circular e uniforme ao redor de um objeto estelar. Comparamos os resultados obtidos nessa abordagem com os resultados encontrados via teoria quântica de campos no espaço-tempo de Schwarzschild. Obtemos que, do infinito até a órbita circular marginalmente estável (R = 6M) o potencial que melhor reproduz os resultados de Schwarzschild é o de Nowak e Wagoner. Já entre esta órbita e a última órbita circular instável (R = 3M) nenhum dos potenciais pseudo-newtonianos produz resultados satisfatórios, e o potencial newtoniano mostra-se como a melhor aproximação. O potencial de Paczynski e Wiita, o mais utilizado na literatura para analisar discos de acre ção, gerou os resultados menos satisfatórios em nossa análise.

Propondo um problema de forrageio como meio de enriquecer o cativeiro: um estudo comparativo de duas espécies de primatas brasileiros (Callithrix penicillata e Saguinus imperator)

O presente estudo analisou a influência de um simples aparelho que simulou o forrageamento ativo por insetos, funcionando como um alimentador enigmático denominado puzzle, sobre o comportamento de duas espécies de calitriquíneos (Callitrichinae, Primates) mantidos em cativeiro no Centro Nacional de Primatas (CENP). Buscou-se comparar a reação das duas espécies frente ao aparelho proposto, e verificar a viabilidade deste artifício no enriquecimento ambiental das condições de cativeiro para estes primatas. Foram utilizados três casais de Saguinus imperator, conhecido por sagui imperador e três casais de Callithrix penicillata, conhecido por mico estrela, as comparações foram feitas entre as espécies e entre os sexos. Foram gravadas através de vídeo cassete seções de observação com duração de três horas, divididas em dois períodos, um controle (ausência do puzzle) e outro experimental (presença do puzzle), com um total de 36 horas de observação para cada casal. As sessões de controle foram usadas para calcular a linha base do orçamento de atividade para comparações com padrões de comportamento durante as sessões experimentais. As fitas foram transcritas, e todos os eventos de comportamento foram cronometrados, e registrados as medidas da freqüência e duração de eventos. Nas duas espécies a manipulação do puzzle não alcançou uma proporção de tempo muito elevada, porém demonstrou uma diferença clara entre as duas espécies, ocupando 3,96% do tempo total dos micos e 1,99% do tempo dos saguis. Os micos gastaram com repouso, durante as seções experimentais, 17% menos tempo em comparação com a situação controle. Considerando que a ociosidade reduziu em menos de 7% nos saguis, com a atividade geral aumentada em 10%. Entre os sexos, notou-se uma inversão, onde as fêmeas de sagui tiveram um aumento de 18% na ociosidade e queda de 14% na atividade, enquanto as fêmeas de mico tiveram uma diminuição nas duas categorias de pouco menos de 10%. Os machos responderam mais ao puzzle, os saguis aumentaram em 58% a atividade e diminuíram em 23% a ociosidade, os micos aumentaram a atividade em 4% e diminuíram a ociosidade em menos de 10%. Todos os animais aprenderam a manipular o equipamento e capturar os insetos, sendo que os micos tiveram um maior êxito no número de larvas capturadas, superando aos saguis em cerca de 54% de seu índice. As fêmeas tiveram mais êxitos em ambas as espécies, porém superficialmente nos micos. Em contraste, nos saguis as fêmeas foram 16% mais prósperas que os machos. O puzzle foi eficiente em enriquecer o ambiente cativo dos animais e estimular o comportamento manipulativo, sendo um instrumento importante na busca do bem estar para as duas espécies.

Atenuação de múltiplas pelo método WHLP-CRS

Nas bacias sedimentares da região Amazônica, a geração e o acúmulo de hidrocarboneto estão relacionados com a presença das soleiras de diabásio. Estas rochas magmáticas intrusivas possuem grandes contrastes de impedância com as rochas sedimentares encaixantes, resultando em múltiplas externas e internas, com amplitudes semelhantes às das reflexões sísmicas primárias. Estas múltiplas podem predominar sobre as informações oriundas de interfaces mais profundas, dificultando o processamento, a interpretação e o imageamento da seção de sísmica. O objetivo da presente tese é realizar a atenuação de múltiplas em seções sintéticas fontecomum (CS), através da combinação dos métodos Wiener-Hopf-Levinson de predição (WHLP) e o do empilhamento superfície-de-reflexão-comum (CRS), aqui denominando pela sigla WHLPCRS. O operador de deconvolução é calculado com as amplitudes reais do sinal sísmico e traço-a-traço, o que consideramos como uma melhor eficiência para a operação de atenuação. A identificação das múltiplas é feita na seção de afastamento-nulo (AN) simulada com o empilhamento CRS, utilizando o critério da periodicidade entre primária e suas múltiplas. Os atributos da frente de onda, obtidos através do empilhamento CRS, são utilizados na definição de janelas móveis no domínio tempo-espaço, e usados para calcular o operador WHLP-CRS. No desenvolvimento do presente trabalho, visamos evitar a inconveniência da seção processada ZO; desenhar e aplicar operadores na configuração CS; e estender o método WHL para camadas curvas.

Influência de estruturas geológicas bidimensionais no campo geoeletromagnético na presença do eletrojato equatorial

A Terra atua como um grande magneto esférico, cujo campo assemelha-se àquele gerado por um dipolo magnético. Este campo apresenta mudanças de intensidade que variam com a localização e a hora local. A parte principal do campo geomagnético se origina no interior da Terra através de processos eletromagnéticos. Extensivos estudos mostraram ainda que existem contribuições de origem externa ao planeta, principalmente de origem solar. Dentre estas fontes há anomalias do campo magnético que surgem a partir de um aumento diurno da corrente elétrica em uma estreita faixa da ionosfera, de direção leste-oeste, centrada no equador magnético e denominada Eletrojato Equatorial (EEJ). Ocasionalmente estas correntes podem apresentar reversões de fluxo, sendo denominadas Contra-Eletrojato (CEJ). Vários autores têm estudado os efeitos do EEJ e CEJ sobre as observações geoeletromagnéticas. Eles estão interessados no efeito combinado do EEJ e estruturas geológicas condutivas 1-D e 2-D. Nestes trabalhos a estrutura 2-D sempre se apresentava paralela ao eletrojato, o que é uma hipótese bastante restritiva ao se modelar ambientes geológicos mais realistas, em que corpos bidimensionais podem ter qualquer strike em relação ao EEJ. Neste trabalho apresentamos a solução deste problema sem esta restrição. Assim, mostramos os campos geoeletromagnéticos devidos a estruturas bidimensionais que possuam strike oblíquo em relação ao EEJ, através de perfis dos campos elétrico e magnéticos calculados na superfície e formando direção arbitrária à heterogeneidade condutiva 2-D. Com esta resposta avaliamos ainda qual a influência que estruturas bidimensionais exercem sobre a resposta magnetotelúrica, sob influência do Eletrojato Equatorial. Durante o desenvolvimento deste trabalho, utilizamos o método de elementos finitos, tendo por fonte eletromagnética o EEJ e o CEJ, que por sua vez foram representados por uma combinação de distribuições gaussianas de densidade de corrente. Estas fontes foram decompostas nas direções paralela e perpendicular à estrutura 2-D, resultando nos modos de propagação TE₁ e TE₂ e TM acoplados, respectivamente. Resolvemos o modo acoplado aplicando uma Transformada de Fourier nas equações de Maxwell e uma Transformada Inversa de Fourier na solução encontrada. De acordo com os experimentos numéricos realizados em um modelo interpretativo da Anomalia Condutiva da Bacia do Parnaíba, formado por uma enorme estrutura de 3000 ohm-m dentro de um corpo externo condutivo (1 ohm-m), concluímos que a presença do CEJ causa uma inversão na anomalia, se compararmos com o resultado do EEJ. Concluímos também que para as frequências mais altas as componentes do campo elétrico apresentam menor influência da parte interna do corpo 2-D do que da parte externa. Já para frequências mais baixas este comportamento se observa com as componentes do campo magnético. Com relação à frequência, vimos os efeitos do “skin-depth”, principalmente nas respostas magnéticas. Além disso, quando a estrutura 2-D está paralela ao eletrojato, o campo elétrico é insensível à estrutura interna do modelo para todos os valores de frequência utilizados. Com respeito ao ângulo θ_h entre a heterogeneidade e a fonte, vimos que o modo TM se manifesta naturalmente quando θ_h é diferente de 0°. Neste caso, o modo TE é composto por uma parte devido à componente da fonte paralela à heterogeneidade e a outra devido à componente da fonte perpendicular, que é acoplada ao modo TM. Assim, os campos calculados têm relação direta com o valor de θ_h. Analisando a influência do ângulo entre a direção do perfil dos campos e o strike da heterogeneidade verificamos que, à medida que θ_h se aproxima de 90°, os campos primários tornam-se variáveis para valores de θ_p diferentes de 90°. Estas variações causam uma assimetria na anomalia e dão uma idéia da inclinação da direção do perfil em relação aos corpos. Finalmente, concluímos que uma das influências que a distância entre o centro do EEJ e o centro da estrutura 2-D, causa sobre as componentes dos campos está relacionado às correntes reversas do EEJ e CEJ, pois a 500 km do centro da fonte estas correntes têm máxima intensidade. No entanto, com o aumento da distância, as anomalias diminuem de intensidade. Nas sondagens MT, nós também usamos o EEJ e o CEJ como fonte primária e comparamos nossos resultados com a resposta da onda plana. Deste modo observamos que as componentes do campo geoeletromagnético, usadas para calcular a impedância, têm influência do fator de acoplamento entre os modos TE₂ e TM. Além disso, esta influência se torna maior em meios resistivos e nas frequências mais baixas. No entanto, o fator de acoplamento não afeta os dados magnetotelúricos em frequências maiores de 10^-2 Hz. Para frequências da ordem de 10^-4 Hz os dados MT apresentam duas fontes de perturbação: a primeira e mais evidente é devido à presença fonte 2-D (EEJ e CEJ), que viola a hipótese da onda plana no método MT; e a segunda é causada pelo acoplamento entre os modos TE₂ e TM, pois quando a estrutura bidimensional está obliqua à fonte 2-D temos correntes elétricas adicionais ao longo da heterogeneidade. Concluimos assim, que o strike de uma grande estrutura condutiva bidimensional relativamente à direção do EEJ ou CEJ tem de fato influência sobre o campo geomagnético. Por outro lado, para estudos magnetotelúricos rasos (frequências maiores de 10^-3 Hz) o efeito do ângulo entre a estrutura geológica 2-D e a direção do EEJ não é tão importante. Contudo, em estudos de litosfera frequências menores de 10^-3 Hz) o acoplamento entre os modos TE₂ e TM não pode ser ignorado.

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.

Comparação dos filtros de velocidade e do operador WHLP-CRS na atenuação de múltiplas

A motivação geológica deste trabalho reside no imageamento de estruturas de bacias sedimentares da região Amazônica, onde a geração e o acúmulo de hidrocarboneto estão relacionados com a presença de soleiras de diabásio. A motivação sísmica reside no fato de que essas rochas intrusivas possuem grandes contrastes de impedância com a rocha encaixante, o que resulta em múltiplas, externas e internas, com amplitudes semelhantes as das primárias. O sinal sísmico das múltiplas podem predominar sobre o sinal das reflexões primárias oriundas de interfaces mais profundas, o que pode dificultar o processamento, a interpretação e o imageamento da seção sísmica temporal. Neste trabalho, estudamos a atenuação de múltiplas em seções sintéticas fonte-comum (FC) através da comparação de dois métodos. O primeiro método resulta da combinação das técnicas Wiener-Hopf-Levinson de predição (WHLP) e o de empilhamento superfície-de-reflexão-comum (CRS), e denominando WHLP-CRS, onde o operador é desenhado exclusivamente no domínio do tempo-espaço. O segundo método utilizado é o filtro de velocidade (ω-k) aplicado após o empilhamento superfície-de-reflexão (CRS), onde o operador é desenhado exclusivamente no domínio bidimensional de freqüência temporal-espacial. A identificação das múltiplas é feita na seção de afastamento-nulo (AN) simulada com o empilhamento CRS, e utiliza o critério da periodicidade entre primária e suas múltiplas. Os atributos da frente de onda, obtidos através do empilhamento CRS, são utilizados na definição de janelas móveis no domínio tempo-espaço, que são usadas para calcular o operador WHLP-CRS. O cálculo do filtroω-k é realizado no domínio da freqüência temporal-espacial, onde os eventos são selecionados para corte ou passagem. O filtro (ω-k) é classificado como filtro de corte, com alteração de amplitude, mas não de fase, e limites práticos são impostos pela amostragem tempo-espaço. Em termos práticos, concluímos que, para o caso de múltiplas, os eventos separados no domínio x-t não necessariamente se separam no domínio ω-k, o que dificulta o desenho de um operador ω-k semelhante em performance ao operador x-t.

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.