Influência de variáveis meteorológicas na produção de liteira na Estação Científica Ferreira Penna, Caxiuanã, Pará

A avaliação da influência da exclusão da precipitação sobre a variação na produção de liteira foi realizada na Reserva Florestal de Caxiuanã – Estação Científica Ferreira Penna (1º 42' S, 51º 31' W) durante o período de março de 2001 a fevereiro de 2003, com o objetivo de identificar as principais variáveis meteorológicas e do balanço hídrico mensal que afetaram a produção e decomposição de liteira. Este trabalho foi parte do subprojeto EXPERIMENTO DE SECA NA FLORESTA (ESECAFLOR), que tinha o objetivo de estudar o impacto de seca prolongada na floresta nos fluxos de água, energia e dióxido de carbono na floresta amazônica. A sazonalidade na produção de liteira e de seus componentes (folhas, gravetos e partes reprodutivas) foi bem estabelecida, com a ocorrência da maior produção de liteira nos meses com menor precipitação. A produção total mensal de liteira variou durante o período experimental, de 294,78 kg.ha-1 a 1758,69 kg.ha-1, com um valor médio de 777,70 kg.ha-1. A produção total de liteira foi distribuída em folhas (61,40 %), gravetos (18,45 %) e partes reprodutivas (20,14 %). Os resultados obtidos na parcela sob condições naturais apresentaram uma maior produção de liteira, aproximadamente 25 %, comparada com a produção de liteira na parcela submetida à exclusão da água de chuva. As variáveis mais fortemente correlacionadas com a produção de liteira e suas componentes foram a velocidade do vento, a densidade de fluxo de radiação solar global, a densidade de fluxo de radiação fotossintéticamente

Variabilidade sazonal da condutância estomática em um ecossistema de manguezal amazônico e suas relações com variáveis meteorológicas

No presente trabalho foram estudadas as variações da condutância estomática (g_s) para o período chuvoso (março) e seco (agosto) do ano de 2003, e suas relações de dependência com algumas variáveis meteorológicas medidas em um ecossistema de manguezal amazônico. As informações utilizadas foram do projeto ECOBIOMA, parte integrante do Experimento de Grande Escala da Biosfera-Atmosfera da Amazônia (LBA). A g_s acompanha a tendência de variação do balanço de radiação, atingindo valores máximos durante o dia e mínimos durante a noite. A condutância apresentou maiores flutuações no período chuvoso, com valor médio de g_s = 0,015 m s^-1, porém com magnitudes inferiores as do período seco. Durante a época seca apresentou um valor médio de g_s = 0,027 m s^-1, com menor amplitude, variando de 0,010 < g_s < 0,042 m s^-1. As variáveis meteorológicas utilizadas para o estabelecimento de relações de dependência com a variabilidade diária de g_s foram déficit de umidade específica (δq), déficit de pressão de vapor (DPV), saldo de radiação (Rn) e velocidade do vento (Vv). O DPV apresentou as melhores correlações com a g_s sendo o R² = 0,99 em ambos os períodos. Apesar de também ser importante nas trocas gasosas entre a vegetação e a atmosfera, a Vv apresentou a menor influência na variação média da g_s, com um R² = 0,44 para época chuvosa e R² =0,51 para o período seco.

Metodologia geofísica aplicada ao estudo arqueológico dos sítios Bittencourt e Jambuaçu, Estado do Pará

Medidas magnéticas, de radiação gama e eletromagnéticas com o Radar de Penetração no Solo (GPR) foram realizadas em dois sítios arqueológicos localizados nas proximidades de uma área destinada a construção de um mineroduto. O objetivo das medidas foi identificar locais para escavação arqueológica visando a descoberta de remanescentes de habitantes que povoaram a região em tempos pretéritos, auxiliando nos estudos para recompor a sua história. As medidas magnéticas evidenciaram anomalias que levaram a descoberta de uma grande quantidade de fragmentos cerâmicos e de outros artefatos. As medidas de radiação gama permitiram a discriminação entre solo arenoso e argiloso, o que representa um resultado importante, pois o solo arenoso geralmente é correlacionado a camada de ocupação arqueológica onde são encontrados os remanescentes. O Radar de Penetração no Solo confirmou as anomalias magnéticas e mostrou outras feições anômalas não relacionadas aos remanescentes arqueológicos.

Metodologia de planejamento para inserção de geradores fotovoltaicos em redes elétricas isoladas e supridas por geradores a diesel

Fontes renováveis de energia baseadas na geração fotovoltaica (GFV) são alternativas energéticas promissoras para a complementação da geração elétrica convencional e centralizada, como usinas térmicas a diesel que suprem potência para redes elétricas isoladas em cidades e localidades remotas na Amazônia. A alocação e o dimensionamento de geradores para aplicação como geração distribuída (GD) é um problema desafiador, com implicações técnicas e econômicas, relacionadas ao planejamento, projeto e operação da rede e, particularmente, a GFV em função das condições ambientais, principalmente radiação solar e temperatura ambiente. Esta tese apresenta uma metodologia analítica para alocar e dimensionar a potência ativa de unidades de geração fotovoltaica, composta pelo gerador FV integrado ao inversor CC/CA (GDFV) para integração, de forma concentrada ou dispersa em redes isoladas de média tensão, e contempla o atendimento de múltiplos objetivos, melhoria do perfil de tensão da rede, redução das perdas ativas e redução da participação da geração a diesel, proporcionando redução no consumo de óleo diesel e, consequentemente, redução da poluição ambiental. A solução global do método proposto constitui um compromisso em relação a esses objetivos, apresentando ponderações diferenciadas para os mesmos, de acordo com prioridades estabelecidas no planejamento do sistema elétrico sob estudo. Para validação da metodologia proposta, foram modeladas e simuladas as redes de 33 e 69 barras do IEEE e um sistema elétrico isolado, cuja usina térmica a diesel supre alimentação para a cidade de Aveiro-PA, Região Amazônica, obtendo-se como resultados dessas simulações melhoria significativa no perfil de tensão, redução nas perdas ativas e na potência de geração a diesel, de acordo com indicadores técnicos que permitem a avaliação quantitativa da integração da GDFV na rede elétrica.

Difração de policristais e difração múltipla de raios X para o estudo da influência dos íons Mn+2, Mg+2 e Cu+2 nas estruturas cristalinas da L asparagina monohidratada e do sulfato de níquel hexahidratado

Neste trabalho realizamos um estudo sobre a influência dos dopantes Mn⁺², Mg⁺² e Cu⁺² nas estruturas cristalinas de cristais de Sulfato de Níquel hexahidratado (NSH) e L Asparagina Monohidratada (LAM). A introdução de dopantes em uma rede cristalina pode alterar suas propriedades físicas ou seu hábito de crescimento. Estas alterações podem favorecer as aplicações tecnológicas destes cristas em diversas áreas como medicina, agricultura, óptica e eletrônica. Os cristais de NSH foram crescidos pelo método da evaporação lenta do solvente e dopados com íons de Mn⁺² e Mg⁺², resultando em cristais de boa qualidade. Realizamos medidas de Difração de raios X de policristais nos cristais puros e dopados e a partir dos resultados obtidos fizemos refinamentos, usando o método de Rietiveld, onde foi observado que os cristais dopados apresentavam a mesma estrutura tetragonal e grupo espacial que o cristal puro, havendo uma pequena mudança em seus parâmetros de rede e volume de suas células unitárias. Observamos que a introdução de dopantes causou alterações nos comprimentos das ligações e nos ângulos entre os átomos de níquel e oxigênio, isso pode explicar porque as temperaturas de desidratação dos cristais de NSH:Mg e NSH:Mn são maiores que a do NSH puro. Usamos a técnica de Difração Mútipla de raios X com radiação síncroton em diferentes energias na estação de trabalho XRD1, do Laboratorio Nacional de Luz Síncroton (LNLS) a fim de identificarmos possíveis mudanças nas estruturas dos cristais dopados de Sulfato de Níquel e de L Asparagina. Os diagramas Renninger mostram mudanças na intensidade, perfil e posições dos picos secundários dos cristais dopados causadas pela introdução dos dopantes. Os cristais de L Asparagina Monohidratada foram crescidos pelo método da evaporação lenta do solvente, sendo dopados com íons de Cu⁺². As medidas de difração múltipla mostram que o cristal dopado possui a mesma estrutura ortorrômbica que o cristal puro. Foram detectadas mudanças nas intensidades, assim como, nas posições e perfil de picos secundários no diagramas Renninger para o cristal dopado. Nossos resultados indicam que o mecanismo de incorporação dos íons de Cu⁺² na rede cristalina da L Asparagina Monohidratada ocorre de forma intersticial.

Fonte escalar acoplada ao campo de Klein-Gordon orbitando um objeto estelar

Neste trabalho determinamos, utilizando Teoria Quântica de Campos em nível de árvore, a radiação escalar emitida por uma fonte em movimento circular uniforme no espaço-tempo plano de Minkowski, assumindo Gravitação Newtoniana, e no espaço-tempo curvo de um buraco negro sem carga e com momento angular nulo, assumindo Relatividade Geral. Efetuamos este cálculo analiticamente para o caso de Minkowski e numericamente no âmbito do espaço-tempo de Schwarzschild, sendo que neste espaço-tempo curvo obtivemos a forma analítica e a normalização dos modos nas regiões assintóticas. Verificamos que, para as órbitas circulares estáveis de acordo com a Relatividade Geral, a potência irradiada no caso de um buraco negro de Schwarzschild é menor do que a obtida no espaço-tempo de Minkowski assumindo a Gravitação Newtoniana. Obtemos também que apenas uma pequena parcela da radiação emitida é absorvida pelo buraco negro. Verificamos que a diferença entre as potências irradiadas em Schwarzschild e Minkowski diminui na medida em que aumentamos o valor da massa do campo. Em Schwarzschild, uma parcela cada vez maior da radiação emitida é absorvida pelo buraco negro na medida em que aumentamos o valor da massa do campo.

Quantização canônica do campo de Proca no espaço-tempo de Rindler e interação de uma fonte uniformemente acelerada com o banho térmico de Unruh

Fazemos a quantização canônica do campo vetorial massivo, primeiro com relação a observadores inerciais e depois com relação a observadores acelerados. Investigamos como uma fonte uniformemente acelerada em Minkowski interage com o campo vetorial massivo no vácuo inercial, através do cálculo da taxa de resposta total. Esta taxa de resposta é calculada em dois referenciais diferentes, um inercial e outro co-acelerado com a fonte. De acordo com o efeito Unruh, no referencial acelerado, o vácuo inercial corresponde a um banho térmico de partículas. Levando em conta este efeito, mostramos, explicitamente, que estas taxas de resposta são idênticas. Este resultado pode ser usado para descrever a interação de elétrons estáticos com partículas Z⁰ presentes na radiação Hawking, desde que os elétrons estejam muito próximos do horizonte de eventos de um buraco negro.

Geotermia rasa em Belém

Um estudo detalhado de geotermia rasa foi realizado no intervalo de profundidade de 0,02 a 210,0 m, na região metropolitana de Belém. As medidas de temperatura foram efetuadas com o uso de termômetros de mercúrio e termistor, enquanto as medidas de condutividade térmica foram feitas em testemunhos de sondagem, utilizando-se o aparato tipo agulha. O fluxo de radiação solar incidente foi registrado com um actinógrafo do tipo Robitzech. O gradiente geotérmico médio obtido na região metropolitana de Belém foi de 0,0254 ± 0,0007 °C.m^-1. O valor médio de condutividade térmica dos testemunhos de sondagem coletados nesta região foi de 1,66 ± 0,52 W.m^-1.°C^-1. Por sua vez, o fluxo geotermal médio na região estudada é de 42,16 ± 1,14 mW.m^-2. A 1,0 m de profundidade ocorrem mudanças temporais nos valores de temperatura, as quais não podem ser consideradas como desprezíveis. Estas, mudanças estão diretamente relacionadas com as variações do fluxo de radiação solar incidente na superfície. O maior incremento deste fluxo, observado de um dia para o outro, foi de cerca de 30 W.m^-2, o que correspondeu a um aumento de temperatura a 1,0 m de profundidade da ordem de 0,22° C. Os perfis de temperatura obtidos durante este trabalho apresentaram deriva em sentidos alternados, durante o ciclo de período de um ano. Os perfis geotérmicos rasos são caracterizados por uma zona de fluxo de calor nulo, denominada "cotovelo", a partir da qual os valores de temperatura crescem com o aumento da profundidade. A profundidade do cotovelo é fundamentalmente influenciada pelo fluxo de calor gerado pela radiação solar incidente na superfície da área em estudo. Elaborou-se um modelo de evolução temporal o qual representa a estrutura térmica e suas variações, da zona compreendida de 0,02 a 10,0 m de profundidade. Os resultados oriundos deste modelo de transferência de calor por condução foram comparados com os obtidos em trabalhos de campo. Observa-se uma boa concordância entre esses resultados. Porém, o ajuste numérico é mais evidente para o período de setembro a fevereiro. Este modelamento poderá ser utilizado também para previsões de deriva dos perfis geotérmicos, desde que sejam conhecidos, a priori, os valores correspondentes às temperaturas médias mensais na superfície, o gradiente geotérmico regional e um dado perfil geotérmico. Este trabalho comprova que o fluxo de radiação solar incidente na superfície é a principal fonte de influência sobre os perfis geotérmicos rasos. A profundidade máxima dessa influência depende principalmente da magnitude desse fluxo, do grau de proteção superficial à incidência direta da radiação e da litologia do local em estudo.

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.

Modelagem numérica da influência do eletrojato equatorial em dados magnetotelúricos produzidos por estruturas tridimensionais

A América do Sul apresenta várias peculiaridades geomagnéticas, uma delas, é a presença do Eletrojato Equatorial, o qual se estende de leste para oeste no Brasil ao longo de aproximadamente 3500 km. Considerando-se o fato de que a influência do Eletrojato Equatorial pode ser detectada a grandes distâncias do seu centro, isto suscita o interesse em se estudar os seus efeitos na exploração magnetotelúrica no Brasil. A influência do eletrojato equatorial na prospecção magnetotelúrica tem sido modelada para meios geológicos uni e bidimensionais valendo-se para isto de soluções analíticas fechadas e de técnicas numéricas tais como elementos finitos e diferenças finitas. Em relação aos meios geológicos tridimensionais, eles tem sido modelados na forma de "camadas finas", usando o algoritmo "thin sheet". As fontes indutoras utilizadas para simular o eletrojato equatorial nestes trabalhos, tem sido linhas de corrente, eletrojatos gaussianos e eletrojatos ondulantes. Por outro lado, o objetivo principal da nossa tese foi o modelamento dos efeitos que o eletrojato equatorial provoca em estruturas tridimensionais próprias da geofísica da prospecção. Com tal finalidade, utilizamos o esquema numérico da equação integral, com as fontes indutoras antes mencionadas. De maneira similar aos trabalhos anteriores, os nossos resultados mostram que a influência do eletrojato equatorial somente acontece em frequências menores que 10^-1 Hz. Este efeito decresce com a distância, mantendo-se até uns 3000 km do centro do eletrojato. Assim sendo, a presença de grandes picos nos perfis da resistividade aparente de um semi-espaço homogêneo, indica que a influência do eletrojato é notável neste tipo de meio. Estes picos se mostram com diferente magnitude para cada eletrojato simulado, sendo que a sua localização também muda de um eletrojato para outro. Entretanto, quando se utilizam modelos geo-elétricos unidimensionais mais de acordo com a realidade, tais como os meios estratificados, percebe-se que a resposta dos eletrojatos se amortece significativamente e não mostra muitas diferenças entre os diferentes tipos de eletrojato. Isto acontece por causa da dissipação da energia eletromagnética devido à presença da estratificação e de camadas condutivas. Dentro do intervalo de 3000 km, a resposta eletromagnética tridimensional pode ser deslocada para cima ou para baixo da resposta da onda plana, dependendo da localização do corpo, da frequência, do tipo de eletrojato e do meio geológico. Quando a resposta aparece deslocada para cima, existe um afastamento entre as sondagens uni e tridimensionais devidas ao eletrojato, assim como um alargamento da anomalia dos perfis que registra a presença da heterogeneidade tridimensional. Quando a resposta aparece deslocada para baixo, no entanto, há uma aproximação entre estes dois tipos de sondagens e um estreitamento da anomalia dos perfis. Por outro lado, a fase se mostra geralmente, de uma forma invertida em relação à resistividade aparente. Isto significa que quando uma sobe a outra desce, e vice-versa. Da mesma forma, comumente nas altas frequências as respostas uni e tridimensionais aparecem deslocadas, enquanto que nas baixas frequências se mostram com os mesmos valores, com exceção dos eletrojatos ondulantes com parâmetros de ondulação α = —2 e —3. Nossos resultados também mostram que características geométricas próprias das estruturas tridimensionais, tais como sua orientação em relação à direção do eletrojato e a dimensão da sua direção principal, afetam a resposta devido ao eletrojato em comparação com os resultados da onda plana. Desta forma, quando a estrutura tridimensional é rotacionada de 90°, em relação à direção do eletrojato e em torno do eixo z, existe uma troca de polarizações nas resistividades dos resultados, mas não existem mudanças nos valores da resistividade aparente no centro da estrutura. Ao redor da mesma, porém, se percebe facilmente alterações nos contornos dos mapas de resistividade aparente, ao serem comparadas com os mapas da estrutura na sua posição original. Isto se deve à persistência dos efeitos galvânicos no centro da estrutura e à presença de efeitos indutivos ao redor do corpo tridimensional. Ao alongar a direção principal da estrutura tridimensional, as sondagens magnetotelúricas vão se aproximando das sondagens das estruturas bidimensionais, principalmente na polarização XY. Mesmo assim, as respostas dos modelos testados estão muito longe de se considerar próximas das respostas de estruturas quase-bidimensionais. Porém, os efeitos do eletrojato em estruturas com direção principal alongada, são muito parecidos com aqueles presentes nas estruturas menores, considerando-se as diferenças entre as sondagens de ambos tipos de estruturas. Por outro lado, os mapas de resistividade aparente deste tipo de estrutura alongada, revelam um grande aumento nos extremos da estrutura, tanto para a onda plana como para o eletrojato. Este efeito é causado pelo acanalamento das correntes ao longo da direção principal da estrutura. O modelamento de estruturas geológicas da Bacia de Marajó confirma que os efeitos do eletrojato podem ser detetados em estruturas pequenas do tipo "horst" ou "graben", a grandes distâncias do centro do mesmo. Assim, os efeitos do eletrojato podem ser percebidos tanto nos meios estratificados como tridimensionais, em duas faixas de freqüência (nas proximidades de 10^-1 Hz e para freqüências menores que 10^-3 Hz), possivelmente influenciados pela presença do embasamento cristalino e a crosta inferior, respectivamente. Desta maneira, os resultados utilizando o eletrojato como fonte indutora, mostram que nas baixas freqüências as sondagens magnetotelúricas podem ser fortemente distorcidas, tanto pelos efeitos galvânicos da estrutura tridimensional como pela presença da influência do eletrojato. Conseqüêntemente, interpretações errôneas dos dados de campo podem ser cometidas, se não se corrigirem os efeitos do eletrojato equatorial ou, da mesma forma, não se utilisarem algoritmos tridimensionais para interpretar os dados, no lugar do usual modelo unidimensional de Tikhonov - Cagniard.

Respostas eletromagnéticas dos arranjos coplanar e coaxial em poço

Foi estudado a viabilidade de aplicação do arranjo coplanar de bobinas nas sondas de perfilagem em poço por indução eletromagnética. Paralelamente foram geradas as respostas do convencional arranjo coaxial, que é o amplamente utilizado nas sondas comerciais, com o propósito de elaborar uma análise comparativa. Através da solução analítica (meios homogêneos) e semi-analítica (meios heterogêneos) foram geradas inicialmente as respostas para modelos mais simples, tais como os do (1) meio homogêneo, isotrópico e ilimitado; (2) uma casca cilíndrica simulando a frente de invasão; (3) duas cascas cilíndricas para simular o efeito annulus; (4) uma interface plana e dois semi-espaços simulando o contato entre duas camadas espessas e (5) uma camada plano-horizontal e dois semi-espaços iguais. Apesar da simplicidade destes modelos, eles permitem uma análise detalhada dos efeitos que alguns parâmetros geoelétricos têm sobre as respostas. Aí então, aplicando ainda as condições de contorno nas fronteiras (Sommerfeld Boundary Value Problem), obtivemos as soluções semi-analíticas que nos permitiram simular as respostas em modelos relativamente mais complexos, tais como (1) zonas de transição gradacional nas frentes de invasão; (2) seqüências de camadas plano-paralelas horizontais e inclinadas; (3) seqüências laminadas que permitem simular meios anisotrópicos e (4) passagem gradacional entre duas camadas espessas. Concluimos que o arranjo coplanar de bobinas pode ser uma ferramenta auxiliar na (1) demarcação das interfaces de camadas espessas; (2) posicionamento dos reservatórios de pequenas espessuras; (3) avaliação de perfis de invasão e (4) localizar variações de condutividade azimutalmente.

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.

Análise da resolução do modelo de Hughes em sondagens magnetotelúricas

Na maioria dos métodos de exploração geofísica, a interpretação é feita assumindo-se um modelo da Terra uniformemente estratificado. Todos os métodos de inversão, inclusive o de dados eletromagnéticos, exigem técnica de modelamento teórico de modo a auxiliar a interpretação. Na literatura os dados são geralmente interpretados em termos de uma estrutura condutiva unidimensional; comumente a Terra é assumida ser horizontalmente uniforme de modo que a condutividade é função somente da profundidade. Neste trabalho uma técnica semi-analítica de modelagem desenvolvida por Hughes (1973) foi usada para modelar a resposta magnética de duas camadas na qual a interface separando as camadas pode ser representada por uma expansão em série de Fourier. A técnica envolve um método de perturbação para encontrar o efeito de um contorno senoidal com pequenas ondulações. Como a perturbação é de primeira ordem a solução obtida é linear, podemos então usar o princípio da superposição e combinar soluções para várias senoides de forma a obter a solução para qualquer dupla camada expandida em série de Fourier. Da comparação com a técnica de elementos finitos, as seguintes conclusões podem ser tiradas: • Para um modelo de dupla camada da Terra, as camadas separadas por uma interface cuja profundidade varia senoidalmente em uma direção, as respostas eletromagnética são muito mais fortes quando a espessura da primeira camada é da ordem do skin depth da onda eletromagnética no meio, e será tanto maior quanto maior for o contraste de condutividade entre as camadas; • Por outro lado, a resistividade aparente para este modelo não é afetada pela mudança na frequência espacial (v) do contorno; • Em caso do uso da solução geral para qualquer dupla camada na Terra cuja interface possa ser desenvolvida em série de Fourier, esta técnica produziu bons resultados quando comparado com a técnica de elementos finitos. A linerização restringe a aplicação da técnica para pequenas estruturas, apesar disso, uma grande quantidade de estruturas pode ser modelada de modo simples e com tempo computacional bastante rápido; • Quando a dimensão da primeira camada possui a mesma ordem de grandeza da estrutura, esta técnica não é recomendada, porque para algumas posições de sondagem, as curvas de resistividade aparente obtidas mostram um pequeno deslocamento quando comparadas com as curvas obtidas por elementos finitos.

Deconvolução de perfis de poço

A maioria dos perfis de poço utilizados nas avaliações petrofísicas de reservatórios possuem uma resolução vertical na ordem de um metro. Isto cria um problema quando as espessuras típicas das camadas são inferiores a um metro, uma vez que não há correção das leituras. Os perfis de alta resolução vertical como da ferramenta de propagação eletromagnética (EPT, Schlumberger), o dipmeter (SHDT, Schlumberger) ou das ferramentas de varredura acústica ou elétrica possuem uma resolução vertical da ordem de centimetros, mas apresentam uma limitada aplicação para as avaliações petrofísicas. Nós apresentamos um método para a deconvolução de um perfil de baixa resolução vertical que utiliza informações de um perfil de alta resolução vertical para identificar uma nítida interface entre camadas que apresentam valores da propriedade petrofísica contrastante, mas localmente constante em ambos os lados. A partir desse intervalo de controle, nós determinamos a função resposta vertical da ferramenta sob as condições atuais do poço com base no teorema da convolução. Utilizamos várias interfaces de modo a obter valores mais representativos da resposta da ferramenta. O perfil de baixa resolução é então deconvoluido utilizando a transformada discreta de Fourier (FFT) sobre todo o intervalo de interesse. É importante destacar que a invasão do filtrado da lama e a presença do bolo de lama não produzem efeitos danosos sobre o método, que foi aplicado a perfis sintéticos e a dados de campo, onde a aplicação de filtros com um correto ajuste de profundidade, bem como a própria escolha do intervalo de controle, antes da deconvolução, são de extrema importância para o sucesso do método.

Desenvolvimento de um conjunto de modelagem de alta resolução para os métodos eletromagnéticos indutivos

A Modelagem Eletromagnética Analógica tem sido uma poderosa ferramenta no estudo das respostas eletromagnéticas de uma grande variedade das técnicas de levantamentos e de modelos geológicos. Nas aplicações de métodos eletromagnéticos nas investigações de contaminação no subsolo, recentemente em grande demanda, são obtidas anomalias muito baixas (da ordem de 0,01% do campo primário). Com o objetivo de estudar em laboratório os problemas de levantamentos eletromagnéticos para alvos de baixo número de indução, um conjunto de modelagem em escala reduzida de alta sensibilidade, foi projetado e construído. Este conjunto opera na faixa de freqüência de 1 kHz até 1 MHz com excelente linearidade. O transdutor de transmissão gera suficiente densidade de fluxo magnético com corrente de até 3 App (pico a pico) circulando nele. O sistema de recepção, isto é, o conjunto de bobinas receptoras mais o pré-amplificador acoplado apresenta sensibilidade 3 V/A/m na freqüência de 10 kHz com inclinação de 20 dB/dec de freqüência. A alta sensibilidade do sistema de recepção combinada com a boa capacidade de geração do sistema transmissor, permite alcançar a resolução de 0,025% do campo primário. O conjunto dispõe ainda, de um sistema mecânico que permite posicionar o sistema de transdutores (transmissor e receptor) em posições arbitrárias com resolução de 1 mm e precisão de 0,05 mm. A precisão no deslocamento linear do sistema de sonda numa linha de 1715 mm de comprimento é alcançada por meio de um motor de passo controlado por um microcontrolador, que é rigidamente acoplado a um fuso que converte o movimento angular do motor em movimento linear. Com este conjunto, é possível simular uma grande variedade de técnicas de levantamentos eletromagnéticos no formato Slingram, isto é, dipolo-dipolo. Além da boa resolução mecânica e boa resolução de medida eletromagnética do conjunto, todo o processo de controle e aquisição de dados é operacionalizado por meio de um computador mestre.