Thermal decomposition and stability in a series of heterobimetallic carbonyl compounds of the type [Fe(CO)4(HgX)2] (X=Cl, Br, I)

Foram preparados compostos carbonílicos heterobimetálicos do tipo [Fe(CO)₄(HgX)₂] ( X= Cl, Br, I), contendo ligação metal-metal, objetivando investigar suas estabilidades térmicas em função do halogênio coordenado aos átomos de mercúrio. A caracterização destes complexos foi feita usando-se de técnicas espectroscópicas de infravermelho e ressonância magnética nuclear, além de análise elementar. O produto final das termodecomposições foi identificado através de espectroscopia no infravermelho e difratograma de raios-X, método de pó.

Potencial alelopático de duas neolignanas isoladas de folhas de Virola surinamensis (Myristicaceae)

Este trabalho teve por objetivos isolar, identificar e caracterizar a atividade alelopática de substâncias químicas presentes nas folhas de Virola surinamensis. O processo de isolamento e identificação das substâncias químicas envolveu o uso de solventes orgânicos e de Ressonância Magnética Nuclear (RMN ¹H, RMN ¹³C e RMN ¹³C-DEPT), espectro de COSY e de HETCOR. A avaliação da atividade alelopática foi realizada em bioensaios de germinação de sementes, em condições de 25 ºC de temperatura constante e fotoperíodo de 12 horas, e de desenvolvimento da radícula e do hipocótilo, com 25 ºC de temperatura constante e fotoperíodo de 24 horas, empregando-se concentrações variando de 1,0 a 8,0 mg L^-1. Como plantas receptoras, foram utilizadas as espécies daninhas Mimosa pudica, Senna obtusifolia e Senna occidentalis. Foram isoladas e identificadas duas neolignanas: a surinamensina e a virolina. A tendência geral observada nos resultados foi de aumento da intensidade dos efeitos alelopáticos inibitórios em função do aumento da concentração, com inibições máximas obtidas, sempre, na concentração de 8,0 mg L^-1. A surinamensina apresentou maior potencial para inibir a germinação e o desenvolvimento da radícula e do hipocótilo do que a virolina, independentemente da espécie receptora e do fator da planta analisado. Considerando-se as intensidades dos efeitos promovidos sobre os três fatores das plantas, o desenvolvimento da radícula e o do hipocótilo foram mais intensamente inibidos pelas duas substâncias do que a germinação das sementes. À exceção dos efeitos verificados sobre o desenvolvimento do hipocótilo, malícia foi a espécie de maior sensibilidade aos efeitos alelopáticos das duas neolignanas, enquanto mata-pasto foi aquela que evidenciou inibições de menor magnitude.

Variação temporal do fitoplâncton e dos parâmetros hidrológicos da zona de arrebentação da Ilha Canela (Bragança, Pará, Brasil)

As variações sazonal e nictemeral do fitoplâncton e dos parâmetros hidrológicos foram estudadas em uma estação fixa (00º46'37.2"S-046º43'24.5"W) localizada na Ilha Canela (Bragança-Pará) em setembro e dezembro de 2004 e março e junho de 2005. Amostras subsuperficiais de água foram coletadas para os estudos qualitativos, quantitativos e para a determinação das concentrações de clorofila a. Simultaneamente, foram medidos os parâmetros físico-químicos da água: salinidade, temperatura, pH, oxigênio dissolvido e percentual de saturação. Foram identificados 64 táxons entre Cyanophyta (um táxon), Bacillariophyta (54 táxons) e Dinophyta (nove táxons). A concentração de clorofila a variou de 4,67 mg m^-3 (período seco) a 5,44 mg m^-3 (período chuvoso) e acompanhou a densidade fitoplanctônica, que foi mais elevada durante o período chuvoso (média=1.870 x 10³ cél L^-1). Os fitoflagelados dominaram quantitativamente o fitoplâncton local, seguidos pelas diatomáceas. Dimeregramma minor e Skeletonema sp. constituíram espécies muito freqüentes e muito abundantes. A ressuspensão de sedimentos provocadas pelos intensos ventos e a arrebentação das ondas favoreceram a dominância de D. minor durante o período seco. No período chuvoso, a elevada pluviometria, os moderados ventos, bem como a influência das águas estuarinas do Taperaçu e Caeté, propiciaram a redução da salinidade e o desenvolvimento de outras espécies fitoplanctônicas.

Configurações magnéticas colineares e não-colineares em currais de Fe, Cr e Mn sobre a Pt (111)

Neste trabalho investigamos as propriedades magnéticas de currais de Fe, Cr e Mn adsorvidos sobre a superfície de Pt(111) utilizando o método RS-LMTO-ASA (Real Space Linear Muffin Tin Orbital - Atomic Sphere Approximation), o qual é um método de primeiros princípios baseado na Teoria do Funcional da Densidade (DFT-Density Functional Theory), que permite o cálculo de estruturas magnéticas não-colineares. Obtivemos que os átomos de Fe apresentam momentos magnéticos elevados, da ordem de 3.5µB /átomo, e têm uma interação de troca entre primeiros vizinhos forte e ferro-magnética. Isto leva a um arranjo magnético colinear no curral. Para os currais de Mn e Cr encontramos que estes possuem elevado momento magnético, da ordem de 4.51µB /átomo e 4.15µB /átomo, respectivamente, e interações de troca entre primeiros vizinhos antiferro-magnéticas. Isto conduz a arranjos magnéticos colineares em currais simples, assim como interessantes ordenamentos não-colineares, tais como estruturas tipo vértice (skyrmions), para os currais com uma geometria particular onde o antiferromagnetismo se apresenta frustado.

Investigação da transição isolante-metal do CDM via algoritmos computacionais e estudo de nano-dispositivos orgânicos

Apresentamos neste trabalho um estudo teórico sobre polímeros orgânicos conjugados. É conhecido que estes sistemas, em geral semicondutores ou isolantes, sob dopagem química podem vir a adquirir propriedades elétricas de material condutor. E ainda, sob ação de campo elétrico, pequenos oligômeros podem apresentar comportamento equivalente ao de dispositivos usuais, mas com inúmeras vantagens como, por exemplo, tamanho extremamente reduzido (alguns nanômetros). Dessa forma no primeiro capítulo faremos uma breve introdução sobre polímeros orgânicos conjugados mostrando alguns resultados experimentais obtidos para o polímero 4-dicianometileno-4,4-ciclopenta [2,1-b: 3,4b’] ditiofeno – CDM, que é o objeto central de estudo desta dissertação. O capítulo 2 trata dos métodos quânticos utilizados. Citaremos a Teoria de Hartre-Fock (HF) e suas derivações semi-empíricas. A técnica de Interação de configuração (CI) e a Teoria do Funcional da Densidade (DFT) também serão tratadas neste capítulo. O capítulo 3 é dedicado a descrever as características de alguns dispositivos usuais como diodos e transistores. Aqui o fundamental é entender a composição, o funcionamento e principalmente, como se comportam suas curvas características corrente versus voltagem (IxV). Citaremos ainda alguns dispositivos eletrônicos extremamente pequenos. No capítulo 4 começa nossos resultados e discussões referentes a análise da transição isolante-metal em CDM sob ação de dopagem. Primeiramente a nível semiempírico, obtivemos a caracterização ótica de oligômeros de CDM neutro e na presença de defeitos conformacionais do tipo bipólarons negativo e positivo. Partindo de geometrias otimizadas via métodos AM1 e PM3 obtivemos o espectro de absorção para sistemas com e sem carga. A nível Hartree-Fock calculamos a Estrutura de Bandas e a Densidade de Estados (DOS) para o PCDM no estado neutro e dopado. O cálculo da DOS e da Dispersão foram realizados através de programas computacionais desenvolvidos aqui no Grupo de Física de Materiais da Amazônia (GFMA). Apresentamos ainda neste capítulo o espectro de absorção teórico para oligômeros de CDM com diversas configurações com geometrias totalmente otimizadas pelo DFT. No capítulo 5 temos os resultados relativos à análise de nanodispositivos baseados em tetrâmeros de CDM com e sem carga. As curvas do deslocamento de carga versus voltagem apresentam características de curvas de dispositivos usuais. Analisamos também o espectro de absorção teórico dos nanodispositivos para valores de tensão nula e em pontos de saturação de corrente nas regiões direta e reversa.

Propriedades magnéticas de nanoestruturas de metais de transição 3d em superfícies de Pd

Motivados por estudos experimentais acerca de monocamadas de metais de transição 3d sobre superfícies de Pd, nesta dissertação investigamos o complexo magnetismo de nanoestruturas, embebidas ou adsorvidas, em superfícies metálicas através de cálculos de primeiros princípios. Utilizamos o método RS-LMTO-ASA (Real Space - Linear MuffinTin Orbital - Atomic Sphere Approximation), o qual é baseado na teoria do funcional da densidade (DFT - Density Functional Theory) e implementado para o cálculo de estruturas magnéticas não colineares. Com este propósito, investigamos nanoestruturas embebidas e ligas (2 x 2) de metais 3d (Cr, Mn, Fe, Co e Ni) na superfície Pd (110), além de nanoestruturas de Cr adsorvidas sobre a superfície de Pd (111). Primeiro, para as nanoestruturas embebidas na superfície Pd (110), analisamos a variação do momento magnético de spin orbital com relação ao número de vizinhos e de valência dos metais 3d. Também mostramos que estas estruturas têm ordenamento magnético colinear, exceto as de Cr e Mn, que apresentam magnetismo não colinear associado à frustração geométrica. Para o caso de nanofios de Cr adsorvidos sobre a superfície de Pd (111), verificamos uma configuração colinear antiferromagnética para cadeias com até 9 átomos. Para o nanofio com 10 átomos obtivemos uma configuração tipo antiferromagnética inclinada (canted). No caso de nanoestruturas de Cr bidimensionais, verificamos complexas configurações magnéticas não colineares com diferentes quiralidades.

Propriedades magnéticas de nanoestruturas adsorvidas em superfícies metálicas

Neste trabalho, utilizamos o método de primeiros princípios RS-LMTO-ASA (Real Space – Linear Muffin-Tin Orbital - Atomic Sphere Approximation) baseado na Teoria do Funcional da Densidade (DFT - Density Functional Theory) e implementado para o cálculo de estruturas magnéticas não-colineares, para investigar as propriedades magnéticas de nanoestruturas adsorvidas em superfícies metálicas. Consideramos aglomerados com diferentes geometrias e tamanhos como adátomos, dímeros, trímeros, nanofios e nanoestruturas de geometria triangular de Fe, Fe-Co e Fe-Pt adsorvidos sobre a superfície de Pt(111) e tratamos também nanoestruturas de Mn sobre a superfície de Ag(111). Mostramos que os nanofios de Fe-Co sobre a superfície de Pt(111) apresentam um ordenamento ferromagnético. Devido à redução do número de coordenação presente na superfície, os momentos de spin e orbital nos sítios de Fe e Co mostram-se elevados comparados com os respectivos valores dos momentos destes metais como bulk. Analisamos também como estes momentos variam em função da concentração destes elementos nos nanofios. Para os sistemas compostos por nanofios Fe-Pt adsorvidos em Pt(111), mostramos que é possível sintonizar as interações de troca entre os adátomos magnéticos Fe através da introdução de um diferente número de átomos Pt para ligá-los. Por exemplo, a interação de troca entre os adátomos de Fe pode ser consideravelmente aumentada pela introdução de cadeias de Pt que os conectem e tanto configurações ferromagnéticas, antiferromagnéticas ou não-colineares entre os adátomos de Fe podem ser estabilizadas, dependendo da espessura do espaçador Pt. Para os aglomerados Mn sobre a Ag(111) mostramos que a interação de troca entre os sítios de Mn depende não somente da distância entre os átomos, mas também do número de coordenação de cada sítio. Desta forma, verificamos um magnetismo não-colinear nestas nanoestruturas causado tanto por frustração geométrica, quanto pela competição de interações de curto e longo alcance. Nossos resultados estão em boa concordância com os resultados experimentais da literatura e com os resultados teóricos obtidos por outros métodos, quando existentes.

Propriedades magnéticas de nanoestruturas de metais de transição 3d adsorvidas na superfície de Pt(111)

Neste trabalho, utilizamos o método de primeiros princípios, RS-LMTO-ASA (“Real Space - Linear Muffin-Tin Orbital - Atomic Sphere Approximation”), baseado na Teoria do Funcional da Densidade (DFT) e implementado para o cálculo de estruturas magnéticas não-colineares, para investigar as propriedades magnéticas de nanoestruturas de metais de transição 3d (Cr, Mn, Fe, Co e Ni) adsorvidas na superfície de Pt(111). Diferentes geometrias como adátomos, dímeros, trímeros, fios lineares e zig-zag foram consideradas e, o tamanho dos aglomerados foi variado de 2 a 7 átomos. Mostramos que os aglomerados de Fe, Co e Ni sobre a superfície de Pt(111), para todas as geometrias simuladas, apresentam um ordenamento ferromagnético. Devido à redução do número de coordenação presente na superfície, os momentos de spin e orbital nos sítios de Fe, Co e Ni, para as diferentes geometrias, mostram-se elevados comparados com os respectivos valores dos momentos destes metais como bulk. Para os glomerados de Cr e Mn mostramos que a interação de troca antiferromagnética entre primeiros vizinhos leva a um ordenamento antiferromagnético colinear no caso de geometrias lineares. No entanto, se o antiferromagnetismo é frustrado por restrição geométrica imposta aos aglomerados pela superfície triangular do substrato, obtém-se um comportamento magnético não-colinear para aglomerados de Cr e Mn sobre a Pt(111). Nossos resultados estão em boa concordância com os resultados experimentais da literatura e com os resultados teóricos obtidos por outros métodos, quando existentes.

Interpolação de dados de campo potencial através da camada equivalente

O uso da técnica da camada equivalente na interpolação de dados de campo potencial permite levar em consideração que a anomalia, gravimétrica ou magnética, a ser interpolada é uma função harmônica. Entretanto, esta técnica tem aplicação computacional restrita aos levantamentos com pequeno número de dados, uma vez que ela exige a solução de um problema de mínimos quadrados com ordem igual a este número. Para viabilizar a aplicação da técnica da camada equivalente aos levantamentos com grande número de dados, nós desenvolvemos o conceito de observações equivalentes e o método EGTG, que, respectivamente, diminui a demanda em memória do computador e otimiza as avaliações dos produtos internos inerentes à solução dos problemas de mínimos quadrados. Basicamente, o conceito de observações equivalentes consiste em selecionar algumas observações, entre todas as observações originais, tais que o ajuste por mínimos quadrados, que ajusta as observações selecionadas, ajusta automaticamente (dentro de um critério de tolerância pré-estabelecido) todas as demais que não foram escolhidas. As observações selecionadas são denominadas observações equivalentes e as restantes são denominadas observações redundantes. Isto corresponde a partir o sistema linear original em dois sistemas lineares com ordens menores. O primeiro com apenas as observações equivalentes e o segundo apenas com as observações redundantes, de tal forma que a solução de mínimos quadrados, obtida a partir do primeiro sistema linear, é também a solução do segundo sistema. Este procedimento possibilita ajustar todos os dados amostrados usando apenas as observações equivalentes (e não todas as observações originais) o que reduz a quantidade de operações e a utilização de memória pelo computador. O método EGTG consiste, primeiramente, em identificar o produto interno como sendo uma integração discreta de uma integral analítica conhecida e, em seguida, em substituir a integração discreta pela avaliação do resultado da integral analítica. Este método deve ser aplicado quando a avaliação da integral analítica exigir menor quantidade de cálculos do que a exigida para computar a avaliação da integral discreta. Para determinar as observações equivalentes, nós desenvolvemos dois algoritmos iterativos denominados DOE e DOEg. O primeiro algoritmo identifica as observações equivalentes do sistema linear como um todo, enquanto que o segundo as identifica em subsistemas disjuntos do sistema linear original. Cada iteração do algoritmo DOEg consiste de uma aplicação do algoritmo DOE em uma partição do sistema linear original. Na interpolação, o algoritmo DOE fornece uma superfície interpoladora que ajusta todos os dados permitindo a interpolação na forma global. O algoritmo DOEg, por outro lado, otimiza a interpolação na forma local uma vez que ele emprega somente as observações equivalentes, em contraste com os algoritmos existentes para a interpolação local que empregam todas as observações. Os métodos de interpolação utilizando a técnica da camada equivalente e o método da mínima curvatura foram comparados quanto às suas capacidades de recuperar os valores verdadeiros da anomalia durante o processo de interpolação. Os testes utilizaram dados sintéticos (produzidos por modelos de fontes prismáticas) a partir dos quais os valores interpolados sobre a malha regular foram obtidos. Estes valores interpolados foram comparados com os valores teóricos, calculados a partir do modelo de fontes sobre a mesma malha, permitindo avaliar a eficiência do método de interpolação em recuperar os verdadeiros valores da anomalia. Em todos os testes realizados o método da camada equivalente recuperou mais fielmente o valor verdadeiro da anomalia do que o método da mínima curvatura. Particularmente em situações de sub-amostragem, o método da mínima curvatura se mostrou incapaz de recuperar o valor verdadeiro da anomalia nos lugares em que ela apresentou curvaturas mais pronunciadas. Para dados adquiridos em níveis diferentes o método da mínima curvatura apresentou o seu pior desempenho, ao contrário do método da camada equivalente que realizou, simultaneamente, a interpolação e o nivelamento. Utilizando o algoritmo DOE foi possível aplicar a técnica da camada equivalente na interpolação (na forma global) dos 3137 dados de anomalia ar-livre de parte do levantamento marinho Equant-2 e 4941 dados de anomalia magnética de campo total de parte do levantamento aeromagnético Carauari-Norte. Os números de observações equivalentes identificados em cada caso foram, respectivamente, iguais a 294 e 299. Utilizando o algoritmo DOEg nós otimizamos a interpolação (na forma local) da totalidade dos dados de ambos os levantamentos citados. Todas as interpolações realizadas não seriam possíveis sem a aplicação do conceito de observações equivalentes. A proporção entre o tempo de CPU (rodando os programas no mesmo espaço de memória) gasto pelo método da mínima curvatura e pela camada equivalente (interpolação global) foi de 1:31. Esta razão para a interpolação local foi praticamente de 1:1.

Processamento interpretativo de dados magnetométricos e inversão de dados gravimétricos aplicados à prospecção de hidrocarbonetos

Apresentamos dois métodos de interpretação de dados de campos potenciais, aplicados à prospecção de hidrocarbonetos. O primeiro emprega dados aeromagnéticos para estimar o limite, no plano horizontal, entre a crosta continental e a crosta oceânica. Este método baseia-se na existência de feições geológicas magnéticas exclusivas da crosta continental, de modo que as estimativas das extremidades destas feições são usadas como estimativas dos limites da crosta continental. Para tanto, o sinal da anomalia aeromagnética na região da plataforma, do talude e da elevação continental é amplificado através do operador de continuação analítica para baixo usando duas implementações: o princípio da camada equivalente e a condição de fronteira de Dirichlet. A maior carga computacional no cálculo do campo continuado para baixo reside na resolução de um sistema de equações lineares de grande porte. Este esforço computacional é minimizado através do processamento por janelas e do emprego do método do gradiente conjugado na resolução do sistema de equações. Como a operação de continuação para baixo é instável, estabilizamos a solução através do funcional estabilizador de primeira ordem de Tikhonov. Testes em dados aeromagnéticos sintéticos contaminados com ruído pseudo-aleatório Gaussiano mostraram a eficiência de ambas as implementações para realçar os finais das feições magnéticas exclusivas da crosta continental, permitindo o delineamento do limite desta com a crosta oceânica. Aplicamos a metodologia em suas duas implementações a dados aeromagnéticos reais de duas regiões da costa brasileira: Foz do Amazonas e Bacia do Jequitinhonha. O segundo método delineia, simultaneamente, a topografia do embasamento de uma bacia sedimentar e a geometria de estruturas salinas contidas no pacote sedimentar. Os modelos interpretativos consistem de um conjunto de prismas bidimensionais verticais justapostos, para o pacote sedimentar e de prismas bidimensionais com seções verticais poligonais para as estruturas salinas. Estabilizamos a solução, incorporando características geométricas do relevo do embasamento e das estruturas salinas compatíveis com o ambiente geológico através dos estabilizadores da suavidade global, suavidade ponderada e da concentração de massa ao longo de direções preferenciais, além de vínculos de desigualdade nos parâmetros. Aplicamos o método a dados gravimétricos sintéticos produzidos por fontes 2D simulando bacias sedimentares intracratônicas e marginais apresentando densidade do pacote sedimentar variando com a profundidade segundo uma lei hiperbólica e abrigando domos e almofadas salinas. Os resultados mostraram que o método apresenta potencial para delinear, simultaneamente, as geometrias tanto de almofadas e domos salinos, como de relevos descontínuos do embasamento. Aplicamos o método, também, a dados reais ao longo de dois perfis gravimétricos sobre as Bacias de Campos e do Jequitinhonha e obtivemos interpretações compatíveis com a geologia da área.

Paleomagnetismo de rochas vulvânicas do Nordeste do Brasil e a época da abertura do Oceano Atlântico Sul

Na primeira parte deste trabalho foram desenvolvidos estudos de magnetismo de rochas e paleomagnetismo em amostras de rochas vulcânicas do Nordeste brasileiro. As idades das amostras compreende os períodos Jurássico e Cretáceo. Com este objetivo foram amostradas quatro áreas tendo sido estudado um total de 496 amostras em 55 sítios. Para a coleta foi utilizada uma perfuradora portátil que extrai amostras de 2.5 cm de diâmetro. A orientação das amostras foi feita por meio de uma bússola magnética e de um clinômetro. Os espécimes foram submetidos a desmagnetizações por campo magnético alternado e em alguns poucos casos foi empregada a desmagnetização térmica. Atribuindo-se peso unitário a cada sítio foi determinada a direção média da magnetização remanescente característica de cada uma das áreas estudadas. As rochas vulcânicas do período Jurássico, localizadas na borda oeste da Bacia do Maranhão (Porto Franco-Estreito), apresentaram uma direção media em que D= 3.9°, I= -17.9° com α₉₅= 9.3°, k= 17.9, N= 15 e todos os sítios apresentaram polaridade normal. Para esta área foi determinado o polo paleomagnético de coordenadas 85.3°N, 82.5°E (A₉₅= 6.9º) que se localiza próximo a outros polos paleomagnéticos conhecidos para esse período. As rochas da borda leste da Bacia do Maranhão (Teresina-Picos-Floriano) de idade cretácica inferior apresentaram uma direção média de magnetização remanescente característica tal que D= 174.7°, I= +6.0º com α₉₅= 2.8º, k= 122, N= 21 e todos os sítios apresentaram polaridade reversa. O polo paleomagnético associado a elas apresentou por coordenadas 83.6°N, 261.0°E (A₉₅=1.9°) e mostrou concordância com outros polos sul americanos de mesma idade. No Rio Grande do Norte foi estudado um enxame de diques toleíticos também de idade cretácica inferior, cuja direção média da magnetização remanescente característica encontrada foi D= 186.6º, I= +20.6º com α₉₅= 14.0° e k= 12.9, N= 10. Os sítios desta área apresentaram magnetizações com polaridades normal a reversa. O polo paleomagnético obtido se localiza em 80.6°N e 94.8°E com A₉₅= 9.5°. O estudo das rochas vulcânicas da província magnética do Cabo de Santo Agostinho indicou para a região um valor de D= 0.4º, I= -20.6º com α₉₅= 4.8° e k= 114, N= 9 para a magnetização remanescente característica. Todos os sítios apresentaram polaridade normal e o polo paleomagnético determinado apresentou as seguintes coordenadas: 87.6ºN, 135ºE com A₉₅= 4.5º. Foi discutida a eliminação da variação secular das direções obtidas, de forma que cada polo apresentado nesta dissertação é verdadeiramente um polo paleomagnético. A análise dos minerais magnéticos portadores da remanência, efetuada por curvas termomagnéticas ou por difração de Raio-X, indicou na maior parte das ocorrências, a presença de titanomagnetita pobre em titânio. A presença de maguemita e algumas vezes hematita, na maior parte das vezes resultado de intemperismo, não anulou a magnetização termoremanente associada à época de formação da rocha, que foi determinada após a aplicação de técnicas de desmagnetização aos espécimes. Pelas curvas termomagnéticas obteve-se, para a maioria das amostras, uma temperatura de Curie entre 500 e 600ºC. Os casos mais freqüentes indicaram a ocorrência de titanomagnetita exsolvida, em que foram observadas a presença de uma fase próxima à magnetita e outra fase rica em titânio, próxima à ilmenita, resultado de oxidação de alta temperatura. A segunda parte do trabalho diz respeito à determinação da época de abertura do oceano Atlântico Sul por meio de dados paleomagnéticos. Entretanto ao invés de se utilizar o procedimento comumente encontrado na literatura, e que se baseia nas curvas de deriva polar aparente de cada continente, foi aplicado um teste estatístico que avalia a probabilidade de determinada posição relativa entre os continentes ser válida ou não, para determinado período em estudo. Assim foi aplicado um teste F a polos paleomagnéticos da África e da América do Sul, dos períodos Triássico, Jurássico, Cretáceo Inferior e Cretáceo Médio-Superior, tendo sido estudadas situações que reconstituem a posição pré-deriva dos continentes e configurações que simulem um afastamento entre eles. Os resultados dos testes estatísticos indicaram, dentro de uma probabilidade de erro de menos de 5%, que a configuração pré-deriva de Martin et al (1981) é compatível com os dados paleomagnéticos do Triássico, mas apresenta uma diferença significativa para os paleopolos de Jurássico, Cretáceo Inferior, Cretáceo Médio-Superior. Outras reconstruções pré-deriva testadas apresentaram o mesmo resultado. A comparação entre os polos paleomagnéticos da América do Sul e da África, segundo uma reconstrução que admite uma pequena abertura entre os continentes, como a proposta por Sclater et al (1977) para 110 m.a. atrás, indicou que os dados do Triássico não são compatíveis com este afastamento. Por outro lado os paleopolos do Jurássico e do Cretáceo Inferior, embora mais antigos que a data sugerida pela reconstrução, são consistentes com esta separação dentro de uma probabilidade de erro de menos de 5%. Os dados do Cretáceo Médio-Superior se mostraram consistentes com a reconstrução sugerida para 80 m.a. atrás por Francheteau (1973) e que propõe uma separação maior entre os continentes. Com base na premissa de deslocamentos de blocos continentais rígidos a análise dos resultados obtidos indicou que América do Sul e África estavam unidas por suas margens continentais opostas no período Triássico e que uma pequena separação entre estes continentes, provavelmente devida a uma rutura inicial, ocorreu no Jurássico e se manteve, então, aproximadamente estacionária até o início do Cretáceo Inferior. Esta conclusão difere da maior parte dos trabalhos que discutem a abertura do oceano Atlântico Sul. Os dados do Cretáceo Médio-Superior são compatíveis com um afastamento rápido e significativo entre os continentes naquele período.

Tomografia eletromagnética para caracterização de reservatórios de hidrocarbonetos

Na produção de petróleo é importante o monitoramento dos parâmetros do reservatório (permeabilidade, porosidade, saturação, pressão, etc) para o seu posterior gerenciamento. A variação dos parâmetros dinâmicos do reservatório induz variações na dinâmica do fluxo no reservatório, como por exemplo, perdas na pressão, dificultando o processo de extração do óleo. A injeção de fluidos aumenta a energia interna do reservatório e incrementa a pressão, estimulando o movimento do óleo em direção aos poços de extração. A tomografia eletromagnética poço-a-poço pode se tomar em uma técnica bastante eficaz no monitoramento dos processos de injeção, considerando-se o fato de ser altamente detectável a percolação de fluidos condutivos através das rochas. Esta tese apresenta o resultado de um algoritmo de tomografia eletromagnética bastante eficaz aplicado a dados sintéticos. O esquema de imageamento assume uma simetria cilíndrica em torno de uma fonte constituída por um dipolo magnético. Durante o processo de imageamento foram usados 21 transmissores e 21 receptores distribuídos em dois poços distanciados de 100 metros. O problema direto foi resolvido pelo método dos elementos finitos aplicado à equação de Helmhotz do campo elétrico secundário. O algoritmo resultante é válido para qualquer situação, não estando sujeito às restrições impostas aos algoritmos baseados nas aproximações de Born e Rytov. Por isso, pode ser aplicado eficientemente em qualquer situação, como em meios com contrastes de condutividade elétrica variando de 2 a 100, freqüências de 0.1 a 1000.0 kHz e heterogeneidades de qualquer dimensão. O problema inverso foi resolvido por intermédio do algoritmo de Marquardt estabilizado. A solução é obtida iterativamente. Os dados invertidos, com ruído Gaussiano aditivo, são as componentes em fase e em quadratura do campo magnético vertical. Sem o uso de vínculos o problema é totalmente instável, resultando em imagens completamente borradas. Duas categorias de vínculos foram usadas: vínculos relativos, do tipo suavidade, e vínculos absolutos. Os resultados obtidos mostram a eficiência desses dois tipos de vínculos através de imagens nítidas de alta resolução. Os tomogramas mostram que a resolução é melhor na direção vertical do que na horizontal e que é também função da freqüência. A posição e a atitude da heterogeneidade é bem recuperada. Ficou também demonstrado que a baixa resolução horizontal pode ser atenuada ou até mesmo eliminada por intermédio dos vínculos.

Análise comparativa de três métodos de inversão aplicados a dados gravimétricos e magnéticos em perfil

Tradicionalmente, o método dos mínimos quadrados tem sido empregado na inversão não linear de dados de campo potencial. No caso em que as observações dos campos gravimétrico ou magnético contém apenas ruído Gaussiano. O método dos mínimos quadrados não apresenta problemas. Entretanto, quando as observações são perturbadas por ruído não Gaussiano, ou mesmo por ruído não aleatório, como é o caso de muitos ruídos geológicos, o método dos mínimos quadrados torna-se bastante ineficiente, e métodos alternativos devem ser empregados a fim de produzir interpretações realísticas. Neste trabalho, uma comparação é feita entre os métodos dos mínimos quadrados, dos mínimos absolutos e do ajuste-M, aplicados à inversão não linear de dados de campo potencial. A comparação é efetuada usando-se dados teóricos, onde diversas situações geológicas são simuladas. Os resultados mostram que na presença de ruído geológico, caracterizado por pequeno corpo raso acima do corpo principal, ou por corpo grande, adjacente ao corpo principal, o ajuste-M apresenta desempenho muito superior ao dos mínimos quadrados e dos mínimos absolutos. Na presença de ruído Gaussiano, entretanto, o ajuste-M tem um desempenho inferior aos outros dois métodos. Como o ruído Gaussiano é um ruído branco, parte dele pode ser removido por um filtro passa baixa adequado, sem muita perda do sinal, o que não ocorre com o ruído geológico que contém componentes importantes de baixo número de onda. Desse modo o ajuste-M se torna uma ferramenta importante na interpretação de áreas geologicamente complexas, onde é comum a contaminação das anomalias por ruído geológico. Os três métodos em estudo são aplicados a uma anomalia magnética real causada por uma intrusão de diabásio em forma de dique, em sedimentos arenosos da formação Piauí na Bacia do Parnaíba. Os três métodos apresentaram resultados semelhantes indicando que tanto o nível de ruído Gaussiano como geológico são baixos nesta anomalia.

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.

Aplicação dos métodos gravimétrico e magnetométrico para a definição do comportamento estrutural da faixa de dobramentos Araguaia

O presente trabalho consiste na interpretação de informações gravimétricas e aeromagnetométricas (Projeto Geofísico Brasil-Canadá - PGBC) da região setentrional da Faixa de Dobramentos Araguaia, envolvendo uma área de aproximadamente 129.000km², compreendida entre os meridianos 47°50’W e 50°30’W e paralelos 4°50'S e 9°00'S. Abrange porções sudeste do Estado do Pará, noroeste de Goiás e oeste do Maranhão. Os trabalhos de campo constaram de levantamentos gravimétrico e altimétrico, ao longo de rodovias que constituem a rede viária regional, e os resultados obtidos, após correções e reduções, foram então organizados de forma a constituírem um mapa de anomalias Bouguer, apresentado em escala 1:500.000. O padrão gravimétrico da faixa de dobramentos Araguaia se caracteriza por mostrar feições predominantemente longitudinais, com curvas isoanômalas de direções submeridianas, concordante com o comportamento litológico-estrutural conhecido para a área. Na porção centro-meridional desse mapa aparece uma zona de anomalia negativa, alongada e intensa (menos de -100mgal), a qual sofre duas importantes inflexões para NW, sendo uma correlacionável ao Lineamento Carajás, e a outra na altura do paralelo 5°30'S. Apresenta zonas de altos gravimétricos, notadamente na porção noroeste (domínio cratônico), relacionada à influência das metavulcânicas do Grupo Grão Pará, e na porção nordeste, devida a massas densas introduzidas na crosta. Destaca-se também o relativo alto gravimétrico acompanhando o flanco oriental do eixo de ocorrência das braquidobras, podendo ser devido a massas excedentes colocadas sob essa região, através de esforços tectônicos, produzindo elevações do nível de base da crosta e consequentemente do embasamento. De forma genérica o flanco oriental da faixa mostra um gradiente mais intenso que o ocidental, devido principalmente à configuração geométrica das estruturas nesse setor. Da mesma forma merece destaque o gradiente regional ascendente de sul para norte. A análise das cartas de intensidade magnética do PGBC fornecem condições à individualização de diversos domínios magnéticos, perfeitamente correlacionáveis litológico e estruturalmente com as unidades geotectônicas que compõem o arcabouço regional, assim como suas subunidades. A execução de modelamento simples, quer para as informações gravimétricas como magnetométricas, objetivaram fornecer subsídios semi-quantitativos à interpretação, de forma a auxiliar no estabelecimento do padrão geométrico do embasamento da faixa de dobramentos Araguaia, assim como das estruturas a ela impostas. Assim, por meio do processamento das informações gravimétricas, verificou-se um padrão geométrico na forma de fatias imbricadas, estabelecidas sobre o substrato. A grande anomalia negativa da porção centro-meridional da faixa pode estar relacionada a uma deficiência de massa profunda (depressão na base da crosta) ou ainda a uma associação entre contraste de densidade e espessamento do pacote de supracrustais, nessa porção. O baixo gravimétrico da porção sudeste da área, entre Guarai e Tupirama, é devido a uma depressão do embasamento, da mesma forma como a anomalia negativa associada ao Lineamento Carajás está relacionada a uma estrutura sinformal conhecida.