Região do espaço que mais influencia em medidas eletromagnéticas no domínio da frequência: caso de uma linha de corrente sobre um semi-espaço condutor

Localizar em subsuperfície a região que mais influencia nas medidas obtidas na superfície da Terra é um problema de grande relevância em qualquer área da Geofísica. Neste trabalho, é feito um estudo sobre a localização dessa região, denominada aqui zona principal, para métodos eletromagnéticos no domínio da freqüência, utilizando-se como fonte uma linha de corrente na superfície de um semi-espaço condutor. No modelo estudado, tem-se, no interior desse semi-espaço, uma heterogeneidade na forma de camada infinita, ou de prisma com seção reta quadrada e comprimento infinito, na direção da linha de corrente. A diferença entre a medida obtida sobre o semi-espaço contendo a heterogeneidade e aquela obtida sobre o semi-espaço homogêneo, depende, entre outros parâmetros, da localização da heterogeneidade em relação ao sistema transmissor-receptor. Portanto, mantidos constantes os demais parâmetros, existirá uma posição da heterogeneidade em que sua influência é máxima nas medidas obtidas. Como esta posição é dependente do contraste de condutividade, das dimensões da heterogeneidade e da freqüência da corrente no transmissor, fica caracterizada uma região e não apenas uma única posição em que a heterogeneidade produzirá a máxima influência nas medidas. Esta região foi denominada zona principal. Identificada a zona principal, torna-se possível localizar com precisão os corpos que, em subsuperfície, provocam as anomalias observadas. Trata-se geralmente de corpos condutores de interesse para algum fim determinado. A localização desses corpos na prospecção, além de facilitar a exploração, reduz os custos de produção. Para localizar a zona principal, foi definida uma função Detetabilidade (∆), capaz de medir a influência da heterogeneidade nas medidas. A função ∆ foi calculada para amplitude e fase das componentes tangencial (H_x) e normal (H_z) à superfície terrestre do campo magnético medido no receptor. Estudando os extremos da função ∆ sob variações de condutividade, tamanho e profundidade da heterogeneidade, em modelos unidimensionais e bidimensionais, foram obtidas as dimensões da zona principal, tanto lateralmente como em profundidade. Os campos eletromagnéticos em modelos unidimensionais foram obtidos de uma forma híbrida, resolvendo numericamente as integrais obtidas da formulação analítica. Para modelos bidimensionais, a solução foi obtida através da técnica de elementos finitos. Os valores máximos da função ∆, calculada para amplitude de H_x, mostraram-se os mais indicados para localizar a zona principal. A localização feita através desta grandeza apresentou-se mais estável do que através das demais, sob variação das propriedades físicas e dimensões geométricas, tanto dos modelos unidimensionais como dos bidimensionais. No caso da heterogeneidade condutora ser uma camada horizontal infinita (caso 1D), a profundidade do plano central dessa camada vem dada pela relação p_o = 0,17 δ_o, onde p_o é essa profundidade e δ_o o "skin depth" da onda plana (em um meio homogêneo de condutividade igual à do meio encaixante (σ₁) e a freqüência dada pelo valor de w em que ocorre o máximo de ∆ calculada para a amplitude de H_x). No caso de uma heterogeneidade bidimensional (caso 2D), as coordenadas do eixo central da zona principal vem dadas por d_o = 0,77 r₀ (sendo d_o a distância horizontal do eixo à fonte transmissora) e p_o = 0,36 δ_o (sendo p_o a profundidade do eixo central da zona principal), onde r₀ é a distância transmissor-receptor e δ_o o "skin depth" da onda plana, nas mesmas condições já estipuladas no caso 1D. Conhecendo-se os valores de r₀ e δ_o para os quais ocorre o máximo de ∆, calculado para a amplitude de H_x, pode-se determinar (d_o, p_o). Para localizar a zona principal (ou, equivalentemente, uma zona condutora anômala em subsuperfície), sugere-se um método que consiste em associar cada valor da função ∆ da amplitude de H_x a um ponto (d, p), gerado através das relações d = 0,77 r e p = 0,36 δ, para cada w, em todo o espectro de freqüências das medidas, em um dado conjunto de configurações transmissor-receptor. São, então, traçadas curvas de contorno com os isovalores de ∆ que vão convergir, na medida em que o valor de ∆ se aproxima do máximo, sobre a localização e as dimensões geométricas aproximadas da heterogeneidade (zona principal).

Aspectos fundamentais das medidas e interpretação de registros paleomagnéticos em rochas sedimentares da Formação Longá-Bacia do Parnaíba

Considerando que o número de dados paleomagnéticos para o Paleozóico Superior, principalmente Devoniano da América do Sul, é ainda insuficiente não permitindo que uma Curva de Deriva Polar confiável possa ser construída, foram amostrados 43 níveis estratigráficos da Formação Longá, os quais foram estudados paleomagneticamente com o propósito de contribuir para que uma Curva mais confiável possa ser construída. A amostragem foi feita segundo o método dos blocos ao longo dos perfis Teresina, Barras, Batalha, na rodovia PI-13 e Floriano, Nazaré do Piauí, Oeiras nas rodovias PI-24 e BR-230, no estado do Piauí. Os tratamentos foram iniciados no laboratório de Paleomagnetismo do IAG-USP e complementados no do NCGG-UFPa. Utilizou-se as técnicas de desmagnetização progressiva por campos alternados até 700 e/ou temperaturas até 670-700ºC. A interpretação dos resultados foi feita por meio dos diagramas vetoriais de Zijderveld, pelas curvas J-T/C de variação da intensidade magnética com os campos ou temperaturas, e pelos gráficos de variação na direção do vetor magnetização. Os cálculos de direção média e polos foram feitos segundo a Estatística de Fisher (1953). Foram identificadas 4 direções de magnetização remanente: 1. Uma secundária de origem química (CRM) e polaridade reversa, cujo mineral responsável é a hematita produzida provavelmente por alteração deutérica a partir da magnetita. Esta magnetização (identificada pela letra B) quando datada paleomagneticamente (coordenadas do polo 80ºS, 3°E, A₉₅ = 13.6°) indicou idade correspondente ao intervalo Carbonífero-Permiano. 2. Uma componente isotérmica (IRM) dura de espectro totalmente superposto à magnetização inicial que não foi afetada por nenhum dos tratamentos. Esta magnetização foi denominada D e apresentou direção muito estável em torno do ponto de declinação = 234.23º e inclinação = 41.94º. 3. Um grupo de direções de magnetização de origem viscosa (VRM) moles, identificados pela letra C, cuja direção média é dada por: declinação = 15º e inclinação = -20º. Foram removidas a temperatura entre 300 – 600ºC. 4. Finalmente uma magnetização principal, de polaridade normal, denominada A, provavelmente de origem detrítica (DRM) cujo correspondente polo paleomagnético (de coordenadas: 48ºS, 331.7ºE; A₉₅ = 9.9º) mostrou-se compatível com a idade da Formação (Devoniano Superior). Esta magnetização foi considerada inicial. Os polos paleomagnéticos correspondentes às magnetizações A e B, juntamente com outros da América do Sul, foram rotacionados para a África segundo a configuração pré-deriva de Smith e Hallam (1970) e comparados a polos Africanos e Australianos de mesma idade, mostrando-se coerentes. Suas polaridades também estão em acordo com as escalas magnetoestratigráficas publicadas por Irving e Pullaiah (1976) e Khramov e Rodionov (1981).

Avaliação da influência do eletrojato equatorial na interpretação de dados magnetotelúrico bidimensionais

O campo magnetotelúrico em regiões equatoriais viola a condição de ondas planas por causa de uma fonte fortemente concentrada na direção E-W na ionosfera, denominada eletrojato equatorial. No presente trabalho, procurou-se analisar a resposta magnetotelúrica de fontes que simulam o efeito do eletrojato equatorial. Foram considerados dois tipos de fontes para simular o eletrojato: uma linha infinita de corrente e uma distribuição gaussiana de densidade de corrente em relação a uma das coordenadas horizontais. A resistividade aparente foi obtida através da relação de Cagniard e comparada com os resultados de ondas planas. É mostrada também a comparação entre a fase da impedância na superfície, para os três tipos de fontes (ondas planas, eletrojato gaussiano e linha de corrente). O problema de meios com heterogeneidades laterais foi resolvido em termos de campos secundários, sendo as equações diferenciais solucionadas através da técnica de elementos finitos bidimensionais. Os resultados mostram que o eletrojato tem pouca influência nas respostas (resistividade aparente e fase) de estruturas geológicas rasas. Entretanto, a influência pode ser considerável nas estruturas profundas (maior que 5000 m), principalmente se suas resistividades são altas (maior que 100 Ω.m). Portanto, a influência do eletrojato equatorial deve ser considerada na interpretação de dados magnetotelúricos de bacias sedimentares profundas ou no estudo da crosta terrestre.

Efeitos das variações sazonais do clima tropical úmido sobre as águas e sedimentos de manguezais do estuário do rio Marapanim, costa nordeste do Estado do Pará

ABSTRACT: To evaluate the short period climatic variations impact over Amazonic mangrove ecosystem, studies were carried out on the northeast coast of Para State. Sediments, surficial and interstitial waters were collected and examinated throughout salinity, pH and Eh (mV) measurements; mineralogical determination using X-ray diffraction and electronic microscopy. Chemical analysis of dissolved sulfides, sulfate and chloride, among others were made in samples collected seasonally, under spring and neap tides conditions. The seasonal chloride variations in the interstitial waters shows 20 g/l during the dry season and less of 10 g/l in the season; the concentration in surficial are higher at neap tides (rainy season) and spring tide (dry season). The dissolved sulfides were found only at 10 cm, indicating sediment exposure to the atmospheric oxygen advectives fluxes. The dissolved iron rates increases between 0-10 cm and the pH tends to neutrality. The saturation of interstitial waters at dry season is indicated by evaporitic minerals: gypsum and halite. The pluviometric variations are responsible by gradual changes in the nutrient and physical chemical properties of surficial and interstitial waters, into the saline equilibrium control at coastal waters, in the salinization and desalinization of sediments and the distribution of mangrove vegetation at the estuary. The prolonged exposure of sediments during the dry season and the morphological characteristics contribute to the total or partial oxidation of surficial sediments, modifying the mineralogy of sediments and the physical chemical characteristics of interstitial waters.

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.

Interpretação de anomalias magnéticas usando integrais finitas dos momentos das componentes

O objetivo do trabalho consiste em desenvolver uma metodologia para determinar os parâmetros de um corpo causador de uma anomalia magnética. Essa metodologia baseia-se na utilização de expressões analíticas, deduzidas para as integrais finitas dos momentos da anomalia de intensidade total e das componentes H_x, H_y e H_z ao longo dos três eixos de um sistema de coordenadas cartesianas. Por meio do ajuste entre essas expressões analíticas e as integrais numéricas finitas dos momentos das componentes magnéticas obtidas a partir de um levantamento geofísico, pode-se computar, através de um processo iterativo, os parâmetros de magnetização, posição, profundidade e dimensões de um corpo anômalo. No caso em que são conhecidas as medidas de somente uma componente magnética, ainda é possível aplicar o método pois as demais componentes podem ser obtidas por um esquema de filtragem matemática. A metodologia foi testada com sucesso para modelos dipolo pontual e linha de dipolos os quais são muito utilizados em Geofísica para interpretar anomalias magnéticas produzidas por corpos geológicos que possuem uma ou duas dimensões horizontais muito menores do que a profundidade.

Correlação entre parâmetros estimados pelos testes Colour Assessment and Diagnosis e Cambridge Colour Test na avaliação da discriminação de cores

Os testes Colour Assessment and Diagnosis (CAD) e Cambridge Colour Test (CCT) têm sido amplamente utilizados em pesquisas básicas e clínicas, devido à alta sensibilidade e especificidade de seus resultados. Estes testes utilizam diferentes paradigmas de estimulação para estimar os limiares de discriminação de cor. Pouco se sabe sobre a relação de cada paradigma na avaliação da discriminação de cor nesses testes. Sendo assim, este trabalho objetiva comparar os parâmetros de avaliação da discriminação de cor estimados pelos testes CAD e CCT em sujeitos tricromatas e com discromatopsia congênita. Foram avaliados 59 sujeitos tricromatas e 38 sujeitos discromatópsicos (16 protans, 22 deutans) com idade média de 26,32 ± 8,9 anos. Foram testados 66 sujeitos nos testes CAD e CCT, 29 sujeitos no teste CAD e 2 sujeitos no teste CCT. O fenótipo da visão de cores de todos os sujeitos foi determinado através de uma bateria de testes psicofísicos e a estimativa dos limiares de discriminação de cor foi avaliada pelos testes CAD e CCT. Os dados de limiares de discriminação de cor foram ajustados a funções de elipse. Os critérios analisados para cada sujeito foram: a área da elipse, o ângulo de rotação e tamanho dos vetores protan, deutan e tritan. Para cada um dos parâmetros foi realizada: estatística descritiva, análise da dispersão dos parâmetros entre os testes CAD e CCT e dos parâmetros em conjunto, razão entre os parâmetros, correlação dos parâmetros a três modelos matemáticos e análise de concordância. Os parâmetros de área e tamanho dos vetores deutan e tritan do subgrupo tricromata; área e tamanho do vetor tritan do subgrupo protan; e tamanho dos vetores protan e tritan do subgrupo deutan apresentaram equivalência entre os resultados de ambos os testes. Os parâmetros de área, ângulo de rotação e tamanho dos vetores protan e tritan apresentaram concordância de medidas entre os testes CAD e CCT. Fatores como as localizações distintas das coordenadas centrais dos testes CAD e CCT e a disposição espacial dos vetores no espaço de cor da CIE 1976 no teste CCT podem ter influenciado na determinação de limiares de discriminação cromática de ambos os testes. Apesar de utilizarem paradigmas distintos na configuração da estimulação, os testes CAD e CCT são equiparáveis.

Distribuição de raios, precipitação e casos de tempestades severas na Amazônia Oriental

O objetivo deste é analisar a relação entre raios e precipitação, sobre diferentes tipos de superfície terrestre e através de estudos de caso de tempestades severas, durante o período de dezembro de 2006 à dezembro de 2008. Os estudos foram feitos dentro de uma área definida pelas coordenadas geográficas 0° a 6° S e 50° a 47° W, abrangendo partes dos estados do Pará, Maranhão e extremo norte de Tocantins.Essa região de estudo foi escolhida por estar localizada em uma área central em relação à localização dos sensores da rede de detecção de raios do Sistema de Proteção Amazônia (RDR-SIPAM) e que contém a maior parte das linhas de transmissão de energia elétrica da região. Os dados de chuva foram coletados por 46 pluviômetros digitais das estações meteorológicas do INMET e HIDROWEB- ANA. Nas análises dos dias de tempestades severas, também foram utilizados dados de pluviômetro de alta resolução temporal operado pela CPRM em Belém. Para classificação dos tipos de superfície e em particular da cobertura vegetal, foram utilizadas as informações de banco de dados sobre Recursos Naturais do IBGE. Para um melhor detalhamento sobre a relação entre raios, precipitação e tipos de superfície, a área de estudo foi dividida em três sub-áreas que possuem características fisiográficas distintas. A Área 1 inclui a Ilha do Marajó e a grande massa d'água da foz do Rio Amazonas e oceano Atlântico próximo. A Área 2 possui 60% de cobertura por floresta. A Área 3 tem 62% de cobertura por áreas agrícolas e/ou pastagens. Os sistemas meteorológicos de larga escala produtores de chuva e raios foram analisados tanto para as 3 áreas consideradas, como nos casos de tempestades severas em torno de Belém. Nessas tempestades foi observado, em quase todos os casos, que os raios antecedem as chuvas. Em geral ficou evidenciada a influencia sazonal na densidade de ocorrência de raios sobre os diferentes tipos de superfície. Todavia, esse trabalho também mostrou que as condições meteorológicas determinam mais amplamente as magnitudes dessa variável. Isso pode ser observado com bastante clareza nos estudos de casos de tempestades severas, onde a ZCIT e as LI`s são fatores determinantes da precipitação e conseqüentemente dos raios ocorridos em Belém.