Avaliação de características dos regimes de umidade na flona de Caxiuanã-PA durante o experimento COBRA-PARÁ

Procura-se investigar a validade de um método de classificação de regimes de umidade, baseado na caracterização de diferentes "estados" da Camada Limite Atmosférica Tropical (CLAT), acima de uma área de floresta, de acordo com a metodologia proposta por Mahrt (1991). Para essas análises foram utilizados dados de radiossondagens e de uma torre micrometeorológica, coletados durante o período menos chuvoso da região, obtidos durante o experimento "COBRA-PARÁ" (realizado no período de 30/10 a 15/11 de 2006). A análise dos regimes de umidade consiste na representação em espaço de fase dos dados disponíveis da razão de Bowen (β), em função do parâmetro -h/L (onde h é a altura da camada de mistura turbulenta e L é o comprimento de Obukhov). Dependendo da localização dos dados nesse espaço foi possível caracterizar as seguintes classes: classe I - ar seco e instável; classe II - vento seco predominante; classe III - vento úmido; classe IV - condição úmida e instável; classe V - condensação de vapor d'água na superfície; classe VI - condição estável dominante; e classe VII -formação de orvalho induzido por radiação noturna resfriando a superfície. Das classes mencionadas, aquelas mais freqüentemente observadas em Caxiuanã, foram as III, IV e VI.

Caracterização de regimes de umidade em regiões tropicais: comparação entre floresta e savana

Este trabalho tem como objetivo investigar um método de classificação de regimes de umidade, baseado na caracterização de diferentes “estados” da Camada Limite Atmosférica Tropical (CLAT), tanto acima de uma área de floresta quanto acima de uma área de savana, de acordo com a metodologia proposta por Mahrt (1991). A partir dessa classificação é elaborado um aperfeiçoamento da mesma que incorpora tanto uma análise da estabilidade termodinâmica da CLAT para a área de floresta, quanto uma investigação sobre a contribuição da Energia Potencial Disponível para Convecção (CAPE) para a classificação dos regimes de umidade. Para essas análises foram utilizados dados de radiossondagens e de torres micrometeorológicas, coletados durante o período menos chuvoso de cada região. Esses dados foram obtidos durante experimentos de campo realizados nas áreas de estudo, sendo que para a área de floresta (Caxiuanã) se utilizaram os dados do experimento COBRA-PARÁ (realizado no período de 06 a 13/11 de 2006), enquanto que para a área de savana (Daly Waters) se utilizaram os dados do experimento “KOORIN” (realizado no período de 15/07 a 13/08 de 1974, na Austrália). A comparação entre os regimes de umidade de cada região sugere que, caso houvesse a savanização da Amazônia, a mesma apresentaria uma atmosfera seca, com a maior parte da energia utilizada para o aquecimento dessa atmosfera, com baixos índices de evapotranspiração, menores taxas de precipitação e inexistência de CAPE. Por outro lado, a análise da estabilidade da atmosfera para Caxiuanã mostrou que, contrariamente ao observado em experimentos na Amazônia Ocidental, na Amazônia Oriental, durante o experimento COBRA-PARÁ, os maiores valores de CAPE ocorreram às 18:00 HL, possivelmente, em decorrência da convergência de umidade que provém da baía de Caxiuanã através da circulação de brisa terrestre. Isso indica que nessa região os máximos de CAPE estiveram associados predominantemente aos campos de umidade e não aos de temperatura. Para essas condições de CLAT “perturbada” o espaço de fase proposto por Mahrt (1991) não caracteriza muito bem regimes de umidade associados a grandes valores da CAPE.

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.