Migração FD e FFD com aproximações de grande abertura angular

Métodos de migração baseados na equação da onda unidirecional apresentam limitações no imageamento de refletores com forte mergulho e no tratamento de ondas evanescentes. Utilizando a expansão de Padé complexa na solução da equação da onda unidirecional para imageamento em geofísica é possível tratar ondas evanescentes, estabilizando a migração e dispensando um tratamento especial para absorção nas fronteiras do domínio. Utilizando várias parcelas na expansão de Padé, refletores com forte mergulho podem ser corretamente migrados. Nesta dissertação aplicamos a equação da onda unidirecional com expansão de Padé complexa para implementar dois algoritmos de migração em profundidade pré e pós empilhamento: a migração por diferenças finitas (FD) e a migração por diferenças finitas e Fourier (FFD). O estudo das curvas de dispersão e da resposta impulsiva dos operadores de migração nos permitiu escolher o número adequado de parcelas na expansão de Padé e os coeficientes da expansão que garantem estabilidade à continuação do campo de onda para um ângulo de mergulho máximo prescrito. As implementações foram validadas nos dados Marmousi e no modelo de domo de sal da SEG/EAGE mostrando que refletores com forte mergulho foram corretamente migrados, mesmo na presença de forte variação lateral de velocidade. Esses resultados são comparados com outros métodos de migração baseados na equação da onda unidirecional ressaltando a qualidade da aproximação estudada neste trabalho.

Migração FFD 3D em profundidade usando aproximação de Padé complexa

Implementações dos métodos de migração diferença finita e Fourier (FFD) usam fatoração direcional para acelerar a performance e economizar custo computacional. Entretanto essa técnica introduz anisotropia numérica que podem erroneamente posicionar os refletores em mergulho ao longo das direções em que o não foi aplicado a fatoração no operador de migração. Implementamos a migração FFD 3D, sem usar a técnica do fatoração direcional, no domínio da frequência usando aproximação de Padé complexa. Essa aproximação elimina a anisotropia numérica ao preço de maior custo computacional buscando a solução do campo de onda para um sistema linear de banda larga. Experimentos numéricos, tanto no modelo homogêneo e heterogêneo, mostram que a técnica da fatoração direcional produz notáveis erros de posicionamento dos refletores em meios com forte variação lateral de velocidade. Comparamos a performance de resolução do algoritmo de FFD usando o método iterativo gradiente biconjugado estabilizado (BICGSTAB) e o multifrontal massively parallel direct solver (MUMPS). Mostrando que a aproximação de Padé complexa é um eficiente precondicionador para o BICGSTAB, reduzindo o número de iterações em relação a aproximação de Padé real. O método iterativo BICGSTAB é mais eficiente que o método direto MUMPS, quando usamos apenas um termo da expansão de Padé complexa. Para maior ângulo de abertura do operador, mais termos da série são requeridos no operador de migração, e neste caso, a performance do método direto é mais eficiente. A validação do algoritmo e as propriedades da evolução computacional foram avaliadas para a resposta ao impulso do modelo de sal SEG/EAGE.

Modelagem sísmica por diferenças finitas em meios bidimensionais com difratores

Este trabalho tem por objetivo a modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações. Como parte deste estudo foi feita a análise numérica visando as condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas. Para a validação do programa de diferenças finitas foi feita a comparação cinemática com a Teoria do Raio para um modelo simples. O motivo deste estudo é ter uma melhor compreensão e controle sobre os problemas de modelagem, visando contribuir para a solução de problemas na interpretação de dados sísmicos. Segundo vários autores na literatura geológica, Derby (1877), Evans (1906), Paiva (1929) e Moura (1938). A Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o recente, atingindo espessuras da ordem de 5 km. Os corpos de diabásio, presentes entre os sedimentos paleozóicos, estão dispostos na forma de soleiras, alcançando espessuras de centenas de metros, perfazendo um volume total de 90.000 quilômetros cúbicos. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica, impossibilitando uma melhor interpretação dos horizontes refletores que se encontram abaixo das soleiras. Para representar situações geológicas desse tipo foi usado um modelo acústico de velocidades. Para o cálculo dos sismogramas foi utilizado um programa de diferenças finitas com aproximação de quarta ordem da equação da onda acústica no espaço e no tempo. As aplicações dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas têm melhorado a compreensão sobre a propagação em meios onde existem heterogeneidades significativas, tendo como resultado boa resolução na interpretação dos eventos de reflexão sísmica em áreas de interesse. Como resultado dos experimentos numéricos realizados em meio de geologia complexa, foi observado a influência significativa das múltiplas devido a camada de alta velocidade, o que faz com que haja maior perda de energia dificultando a interpretação dos alvos. Por esta razão recomendo a integração de dados de superfície com os dados de poço, com o objetivo de se obter uma melhor imagem dos alvos abaixo das soleiras de diabásio.

Projeto de antenas multibandas integradas para terminais 4G

O presente trabalho apresenta um projeto e análise de uma antena monopolo planar multibanda operante nas frequências de 2000 MHz, 2500 MHz e 3500 MHz para aplicação em dispositivos móveis 3G e 4G. Dentre essas faixas estão contidos os espectros 3G (1895- 2165 MHz), 4G (2500-2690 MHz) e 4G (3400-3600 MHz) utilizadas pelo LTE e WIMAX 3,5. Para o desenvolvimento desse estudo foi tomado como referência uma antena dual-band de Hua- Ming Chen projetada para operar em sistemas ISM 1800 MHz, redes Wireless e Bluetooth 2400 MHz. O projeto proposto neste trabalho é uma modificação na estrutura da antena dual-band adaptada para tripla ressonância simultânea em frequências utilizadas em redes móveis 4G. Para alcançar o objetivo proposto utilizou-se para simulação o software CST® MICROWAVE STUDIO baseado no método Integração Finita no Domínio, as simulações foram feitas no domínio do tempo. Após obtenção satisfatória dos resultados simulados, um protótipo da antena foi confeccionado para medição e teste, os resultados apresentados são coerentes aos obtidos na simulação, o que valida o método proposto.

Estudo de técnicas de análise modal operacional em sistemas sujeitos a excitações aleatórias com a presença de componente harmônico

A metodologia tradicional de identificação de parâmetros na análise modal de estruturas é realizada a partir de sinais medidos de força de entrada e de movimento de resposta da estrutura em condições laboratoriais controladas. Entretanto, quando é necessária a obtenção dos parâmetros modais de estruturas de máquinas em operação, as condições para controlar e medir a excitação nestas situações impossibilita a realização da análise modal tradicional. Neste caso, o teste modal é realizado utilizando somente dados de resposta do sistema. A Análise Modal Operacional (AMO) é um método de extração modal em que nenhuma excitação artificial necessita ser aplicada ao sistema, utilizando-se a própria excitação operacional como entrada para medição da resposta do sistema. A técnica clássica de Análise Modal Operacional NExT considera, para isso, que a excitação operacional do sistema seja um ruído branco. Esta técnica faz a consideração de que as funções de correlação obtidas de estruturas podem ser consideradas como funções de resposta ao impulso e então métodos tradicionais de identificação modal no domínio do tempo podem ser empregados. Entretanto, caso a excitação operacional contenha componentes harmônicos que se sobressaiam, estes podem ser confundidos como modos naturais do sistema. Neste trabalho é demonstrada que através da função densidade de probabilidade da banda estreita contendo o pico de um modo, é possível identifica-lo como natural ou operacional (proveniente da excitação operacional da estrutura). É apresentada também uma modificação no método de identificação modal Exponencial Complexa Mínimos Quadrados (LSCE), passando a considerar sinais harmônicos de freqüências conhecidas presentes na excitação operacional, em um ensaio utilizando a técnica NExT. Para validação desses métodos, utiliza-se um modelo teórico de parâmetros modais conhecidos analiticamente e como estudo de caso experimental, um sistema formado por uma viga bi-apoiada suportando um motor elétrico com desbalanceamento de massa.

Regimes térmico e hídrico em solos sob ecossistemas naturais e área agrícola no Leste da Amazônia

A temperatura e umidade do solo são variáveis cujo conhecimento é fundamental para determinar os balanços de energia e água na biosfera. Os regimes térmico e hídrico dos solos sob cada ecossistema apresentam variações consideráveis, de acordo com sua mineralogia, o clima local e a vegetação. Nesse contexto, as temperaturas e umidades do solo foram medidas sob três ecossistemas existentes na região leste da Amazônia, a saber: floresta nativa (FLONA Caxiuanã, 01° 42' 30" S e 51° 31' 45" W), pastagem nativa (Soure, 00° 43' 25" S e 48° 30' 29" W) e área agrícola (!garapé-Açu, 01° 07' 59" S e 47° 36' 55" W). Os dados de campo na floresta e na pastagem foram coletados entre dezembro de 2001 e fevereiro de 2005; enquanto que na área agrícola, o monitoramento foi limitado de agosto de 2003 a fevereiro de 2005. Estas observações das variáveis físicas do solo foram analisadas levando em consideração as variáveis meteorológicas medidas simultaneamente tais como o fluxo de radiação solar incidente e a precipitação pluviométrica, que interferem diretamente nas variáveis do solo em cada sitio escolhido para estudo. As temperaturas do solo foram monitoradas por meio de sondas térmicas em profundidades de 0,05; 0,20 e 0,50 m. Fluxímetros de calor mediram esta variável em níveis de profundidade em 0,05 e 0,20 m. A umidade volumétrica do solo na camada superior de 0,30 m foi medida por sensor de sonda dupla por Reflectometria no Domínio do Tempo (TDR) em cada sitio. Foram feitas analises considerando as respostas do solo durante o período seco e chuvoso local, nestes três ecossistemas representativos do leste da Amazônia. Estimativas de difusividade térmica aparente do solo foram feitas pelos métodos da amplitude e da fase usando os dados de propagação do pulso diário de calor nesses solos. Os resultados mostraram valores bem diferentes, porém,no primeiro método pareceu mais confiável e adequado para o modelamento numérico. Como esperado, considerando a sua pouca cobertura vegetal, as temperaturas dos solos nos níveis superficiais, apresentaram grandes variações na pastagem e na área agrícola. Inesperadamente, as temperaturas na profundidade de 0,5 m abaixo da floresta mostraram maiores variações de amplitude que as profundidades de 0,20 e 0,05 m. O modelamento numérico das variações temporais da temperatura, em função da profundidade, para cada solo foi feito através do método harmônico Os resultados mostraram que o primeiro harmônico representou mais de 90% da variação total observada do pulso diário da temperatura da pastagem e área agrícola em 0,2 e 0,05 m de profundidade. Performance similar do modelamento foi observada na floresta nos níveis de 0,05 e 0,20 m. A magnitude dos fluxos de calor abaixo da pastagem e área agrícola atingiram valores seis vezes maiores que aqueles observados sob a floresta. Os resultados mostraram que, para a camada do solo superior de 0,30 m, a umidade volumétrica do solo sob a floresta é maior que sob os outros ecossistemas estudados neste trabalho. Este resultado é devido aparentemente; à proteção da floresta contra a evaporação da superfície do solo. Uma análise do comportamento sazonal e diário das temperaturas e umidade solos em resposta à radiação solar e precipitação é apresentada. Estudos de caso da taxa de perda da umidade do solo depois de significativa recarga de água por eventos de precipitação, também foram analisados. Algumas estimativas diárias de diminuição de água e recarga durante a noite e madrugada por subida de água de camadas subjacentes para a camada de 0.30 m foram feitas. Este trabalho analisou a maior serie temporal dos dados de temperatura e umidade dos solos coletados com alta freqüência de amostragem disponível até o momento, para o leste da Amazônia. Foi possível caracterizar as diferenças dos regimes destas variáveis físicas, abaixo de três ecossistemas importantes desta região. Estudos futuros dos minerais e materiais orgânicos nestes solos, bem como dos índices de área foliar e da biomassa das coberturas vegetais desses ecossistemas, melhoraria a compreensão dos regimes descritos neste trabalho.

Metodologia de discretização espacial para o método radial de interpolação por pontos (RPIM) aplicada para solução numérica das equações de Maxwell

Neste trabalho, foi desenvolvido e implementado um método de discretização espacial baseado na lei de Coulomb para geração de pontos que possam ser usados em métodos meshless para solução das equações de Maxwell. Tal método aplica a lei de Coulomb para gerar o equilíbrio espacial necessário para gerar alta qualidade de discretização espacial para um domínio de análise. Este método é denominado aqui de CLDM (Coulomb Law Discretization Method ) e é aplicado a problemas bidimensionais. Utiliza-se o método RPIM (Radial Point Interpolation Method) com truncagem por UPML (Uniaxial Perfectlly Matched Layers) para solução das equações de Maxwell no domínio do tempo (modo TMz).

Metodologia para estimação da topologia de linhas telefônicas através da análise de medições em uma e em duas portas

Nesta tese é proposta uma metodologia para identificação automática de topologias de linhas telefônicas compostas de uma ou mais seções de linha, as quais são utilizadas em sistemas de linhas digitais de assinante (DSL, do inglês digital subscriber line). Métodos com esta finalidade são fundamentais para a qualificação da linha com o objetivo de instalação do serviço DSL, em especial na atualização para serviços como o VDSL2 ou para adoção de "vetorização". Com o intuito de ser relevante para as operadoras na qualificação de uma rede inteira, composta de milhões de linhas, é fundamental que, além de precisão, métodos de qualificação tenham baixo custo computacional. Os métodos desenvolvidos são baseados nessa premissa e fazem análise da resposta ao impulso e da resposta à reflectometria no domínio do tempo de uma dada linha. Esses sinais são analisados utilizando-se um método de detecção de bordas, baseado em transformada wavelet, para identificar e extrair características de sinal que contenham informação sobre a topologia da linha. Dependendo das características disponíveis, é utilizado um dos três sistemas especialistas desenvolvidos para interpretação dessas informações e identificação da topologia. Estas metodologias são avaliadas através de um conjunto de teste de linhas reais medidas em laboratório. Seus resultados são comparados com os resultados de dois outros métodos implementados a partir da literatura. Os resultados obtidos mostram que os métodos propostos são eficientes na estimação de informações da topologia da linha e possuem reduzido custo computacional quando comparados às implementações das outras técnicas avaliadas.

Avaliação dinâmica da interligação Tucuruí-Manaus-Macapá frente a pequenas e grandes perturbações

Esta dissertação apresenta uma avaliação dinâmica da conexão Tucuruí-Manaus-Macapá ao Sistema Interligado Nacional frente a pequenas e grandes perturbações. Todos os dados utilizados foram obtidos da base de dados do Operador Nacional do Sistema Elétrico. No entanto, os parâmetros dos estabilizadores dos sistemas de potência das usinas de Macapá não estão sintonizados para a interconexão e as usinas de Manaus não apresentam dispositivos suplementares amortecedores. Sendo assim, o objetivo deste trabalho é projetar estabilizadores de sistemas de potência para as usinas hidrelétricas de Coaraci Nunes e Balbina, e também para as usinas térmicas de Mauá e Aparecida, através de técnicas de respostas em frequência baseada no critério de Nyquist. A utilização dos estabilizadores tem como finalidade o amortecimento dos modos eletromecânicos que surgem quando os sistemas isolados do Amapá e do Amazonas passam a ser integrados ao sistema interligado nacional. Os projetos são avaliados em diversos cenários operativos através de simulações no domínio do tempo.

Projeto e implementação em DSP de sistema de comunicação multiportadora baseado no padrão G.fast

As redes de acesso usando cabos de cobre atualmente utilizam bandas de frequência até 30 MHz, especificada no padrão VDSL2. À medida que arquiteturas híbridas de fibra e cobre se tornam mais proeminentes na indústria e academia, torna-se possível utilizar cabos metálicos mais curtos (i.e. até 250 metros) conectando o último ponto de distribuição aos usuários, de modo que frequências mais altas podem ser exploradas para se alcançar taxas de transmissão de dados de 500 Mbps ou mais, como é o caso do padrão G.fast atualmente em desenvolvimento no ITU-T. Nesse trabalho, um simulador no domínio do tempo foi desenvolvido para avaliar a capacidade do sistema G.fast com diferentes tamanhos de extensão cíclica e diferentes topologias de rede especificadas pelo ITU-T. Os resultados das simulações mostram que sistemas G.fast são robustos a bridged taps e capazes de atingir altas taxas de dados para todas as topologias simuladas, provendo suporte à próxima geração de serviços de banda larga. Além disso, esse trabalho descreve o progresso da implementação de um protótipo de modem baseado no padrão G.fast em um ambiente híbrido de DSP multicore e FPGA utilizando kits de avaliação adquiridos pela UFPA. Arquiteturas, protocolos de comunicação e benchmarks são apresentados e avaliados para se chegar à conclusão de que tal protótipo é factível e fornece suporte flexível a várias linhas de pesquisa em banda larga da próxima geração.

Treatment of geophysical data as a non-stationary process

ABSTRACT: The Kalman-Bucy method is here analized and applied to the solution of a specific filtering problem to increase the signal message/noise ratio. The method is a time domain treatment of a geophysical process classified as stochastic non-stationary. The derivation of the estimator is based on the relationship between the Kalman-Bucy and Wiener approaches for linear systems. In the present work we emphasize the criterion used, the model with apriori information, the algorithm, and the quality as related to the results. The examples are for the ideal well-log response, and the results indicate that this method can be used on a variety of geophysical data treatments, and its study clearly offers a proper insight into modeling and processing of geophysical problems.

Identificação das fontes alimentares de mosquitos transmissores da malária na Amazônia brasileira pela técnica de Bloodmeal ELISA

O ensaio imunoenzimático para identificação de repastos sangüíneos apresenta especificidade até nível de gênero, sensibilidade para identificação de repastos sangüíneos parciais e detecção de repastos múltiplos. Foram identificadas as fontes alimentares de 82% dos anofelinos coletados em campo. Foram testados seis anticorpos monoclonais (humano, suíno, cão, bovino, rato e galinha) e destes, apenas o anti- IgG bovino apresentou instabilidade. Obteve-se 55,7% (519/932) dos repastos positivos para sangue humano, o que demonstra, a preferência alimentar destes anofelinos por humanos. Dos 206 mosquitos que apresentaram repasto único, 27,6% foi em humanos. O Anopheles darlingi apresentou 41% dos repastos em humanos e o Índice de Sangue Humano (HBI) no intradomicílio foi de 0,71. An. marajoara apresentou 51,3% dos repastos em humanos, embora tenha sido encontrada em grande quantidade no ambiente extradomiciliar e apresentando HBI no intradomicílio de 0,76. O An. nuneztovari foi a espécie mais abundante, apresentando comportamento exofílico e antropofílico, com 53,8% de repastos em humano e HBI no intradomicílio de 0,65. O método ELISA Sanduíche está implantado, identificando a fonte alimentar das espécies de anofelinos coletadas em campo, exceto para bovino. É a primeira vez que esta técnica é utilizada para determinação de repastos sanguíneos em anofelinos na região Amazônica brasileira. É importante a determinação da fonte alimentar das espécies de anofelinos no sentido de caracterizar o comportamento antropofílico e assim associá-las ou não à transmissão de malária.

Influência da convecção termossolutal nos espaçamentos dendríticos primários durante a solidificação direcional horizontal de ligas Sn-Pb

Este trabalho apresenta um estudo teórico-experimental sobre os efeitos da convecção termossolutal sobre os espaçamentos dendríticos primários de ligas hipoeutéticas Sn-Pb. Assim, um dispositivo experimental de configuração horizontal refrigerado a água foi construído e amostras do referido sistema de ligas foram solidificadas direcionalmente sob condições transientes de extração de calor. Um método teórico-experimental é aplicado para determinar as velocidades de deslocamento da isoterma liquidus, os gradientes térmicos e as taxas de resfriamento. São também apresentados resultados para os respectivos coeficientes de transferência de calor na interface metal/molde os quais foram calculados a partir de uma análise comparativa realizada entre os perfis experimentais de temperatura e valores teóricos fornecidos por um método numérico baseado em volumes finitos. Alguns modelos teóricos para a previsão dos espaçamentos dendríticos primários são comparados com os resultados experimentais obtidos. Finalmente, um estudo comparativo é realizado entre resultados encontrados neste trabalho e aqueles apresentados na literatura para os espaçamentos dendríticos primários das ligas Sn-Pb investigadas quando solidificadas direcionalmente nos sistemas verticais ascendente e descendente, sob as mesmas condições. A análise das microestruturas indica que os espaçamentos dendríticos primários são bastante influenciados pela direção de crescimento do sólido.

Migração por equação de onda em meios anisotrópicos com correção de amplitude

Este trabalho discute dois aspectos da migração em profundidade através da continuação para baixo dos campos de onda: o tratamento de modos evanescentes e a correção da amplitude dos eventos migrados. Estes dois aspectos são discutidos em meios isotrópicos e para uma classe de meios anisotrópicos. Migrações por diferenças finitas (FD) e por diferenças finitas e Fourier (FFD) podem ser instáveis em meios com forte variação lateral de velocidade. Estes métodos utilizam aproximações de Padé reais para representar o operador que descreve a propagação de ondas descendentes. Estas abordagens não são capazes de tratar corretamente os modos evanescentes, o que pode levar à instabilidades numéricas em meios com forte variação lateral de velocidade. Uma solução possível para esse problema é utilizar aproximação de Padé complexa, que consegue melhor representar os modos evanescentes associados às reflexões pós-críticas, e neste trabalho esta aproximação é utilizada para obter algoritmos FD e híbrido FD/FFD estáveis para migração em meios transversalmente isotrópicos com eixo de simetria vertical (VTI), mesmo na presença de forte variação nas propriedades elásticas do meio. A estabilidade dos algoritmos propostos para meios VTI foi validada através da resposta ao impulso do operador de migração e pela sua aplicação na migração de dados sintéticos, em meios fortemente heterogêneos. Métodos de migração por equação de onda em meios heterogêneos não tratam corretamente a amplitude dos eventos durante a propagação. As equações de onda unidirecionais tradicionais descrevem corretamente apenas a parte cinemática da propagação do campo de onda. Assim, para uma descrição correta das amplitudes deve-se usar as equações de onda unidirecionais de amplitude verdadeira. Em meios verticalmente heterogêneos, as equações de onda unidirecionais de amplitude verdadeira podem ser resolvidas analiticamente. Em meios lateralmente heterogêneos, essas equações não possuem uma solução analítica. Mesmo soluções numéricas tendem a ser instáveis. Para melhorar a compensação de amplitude na migração, em meios com variação lateral de velocidade, é proposto uma aproximação estável para solução da equação de onda unidirecional de amplitude verdadeira. Esta nova aproximação é implementada nas migrações split-step e diferenças finitas e Fourier (FFD). O algoritmo split-step com correção de amplitude foi estendido para meios VTI. A migração pré e pós-empilhamento de dados sintéticos, em meios isotrópicos e anisotrópicos, confirmam o melhor tratamento das amplitudes e estabilidade dos algoritmos propostos.

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.