Análise eletromagnética de cabos OPGW utilizando métodos numérico e analítico

Este trabalho apresenta o estudo eletromagnético de cabos OPGW (Optical Ground Wire) os quais têm dupla função: de pára-raios para linhas de transmissão de alta tensão e de canal de comunicação através de fibras ópticas embutidas na estrutura do cabo. Descargas atmosféricas ou curtos-circuitos podem comprometer a integridade do cabo, devido ao aquecimento nas regiões onde há maior concentração de corrente. Para a análise deste problema foram feitos cálculos eletromagnéticos relacionando-os aos efeitos térmicos no cabo. Nesta análise foram consideradas três diferentes geometrias: o modelo de cabo real, o modelo de cabo com camadas homogêneas e o modelo de cabo com uma camada modificada; esta modificação está relacionada à forma geométrica dos fios da armação do cabo. As ferramentas utilizadas em tal estudo foram o software comercial FEMLAB Multiphysics, baseado no método dos elementos finitos, e um método analítico desenvolvido a partir das equações de Maxwell no domínio da freqüência, que foi implementado utilizando o software MATLAB. Os principais resultados deste trabalho são gráficos de distribuição de densidade de corrente na seção reta do cabo para diferentes freqüências, estudo do efeito pelicular e do efeito de proximidade entre os condutores do cabo.

Análise teórica de novos circuladores de 3-portas em cristais fotônicos bidimensionais

Este trabalho apresenta um estudo teórico sobre novos circuladores compactos com 3-portas tipos W e Y, baseados em cristais fotônicos bidimensionais. No circulador tipo Y, os guias de onda que o compõem formam ângulos de 120° entre si (com formato assemelhado ao da letra Y). O circulador tipo W é uma modificação do tipo Y, obtido a partir do reposicionamento de uma das portas entre as outras duas com um ângulo de 60° entre os guias de onda (com formato assemelhado ao da letra W). Os parâmetros geométricos dos cristais foram obtidos dos diagramas de bandas proibidas. O circulador de três portas tipo Y, projetado para operar em frequências de micro-ondas, foi investigado com o objetivo de gerar um protótipo inédito, enquanto que o tipo W, para frequências ópticas, foi investigado para demonstrar a possibilidade de desenvolver um circulador mais compacto em comparação com o tipo Y conhecido. O tipo W pode ser também uma alternativa geométrica mais adequada no design de circuitos integrados. Os modelos são bons no sentido em que possuem elevada isolação (maior que -20 dB em ambos os circuladores) e baixa perda de inserção (maior que -0,5 dB no caso do circulador tipo Y). O circulador tipo W apresenta uma largura de banda de operação em torno de 100 GHz para um nível de -20 dB de isolação, centrado no comprimento de onda de 1,5um. As simulações foram feitas utilizando-se o software comercial COMSOL Multiphysics, o qual se baseia no método dos elementos finitos.

Segurança estrutural de uma ponte ferroviária em concreto armado sobre o rio Mãe Maria

O estudo do comportamento de pontes ferroviárias vem tomando grande destaque na área experimental, devido muitas das vezes haver a necessidade de um aumento dessas cargas móveis. A companhia VALE a partir de um convênio estabelecido com a UFPA (Universidade Federal do Pará), buscou realizar estudos das capacidades resistentes das pontes ferroviárias, que compõem a Estrada de Ferro Carajás (EFC), para futuros casos de carregamentos que passarão de 32 tf para 40 tf, o que provoca mudanças consideráveis no comportamento da estrutura. A falta de manutenção ou uma manutenção irregular das pontes de concreto armado podem resultar em manifestações patológicas que podem comprometer a segurança e ao risco de intervenção do escoamento. Dentre as OAE’s (Obras de Arte Especial) já analisadas, encontra-se a ponte sobre o Rio Mãe Maria localizada no município de Marabá, estado do Pará, transpondo o rio Mãe Maria, identificada como OAE 50A, apresenta um traçado retilíneo executada em concreto armado moldada “in loco”, é constituída por dois vãos hiperestáticos, totalizando 64,20 m, transpondo o leito do rio em dois vãos de 20 m e com seção transversal oferecendo a largura de 5,85 m. Com isso tornou-se indispensável para avaliação de sua integridade física a monitoração com extensômetros elétricos de resistência (EER) dos elementos (laje, longarina, blocos, encontros e tubulões) que compões a estrutura durante a passagem dos trens carregados de minério. Além da monitoração, foram realizados ensaios não-destrutivos como: dureza superficial através da esclerometria e a estimativa do cobrimento da armadura e localização da mesma por meio da pacometria, ensaios destrutivos como: extração de corpos-de-prova para determinar a compressão diametral e o módulo de elasticidade do concreto e a resistência a penetração através do teste de carbonatação inerentes a questão de durabilidade do mesmo, que serão comparados com a memória de cálculo e com os modelos computacionais tridimensionais desenvolvidos via “software” comercialSAP2000® (COMPUTERS AND STRUCTURES), que utiliza o método dos elementos finitos (MEF). O MEF foi utilizado para obtenção dos esforços permitindo as verificações de segurança de acordo com as recomendações normativas, tais análises comparativas são apresentadas e discutidas para a conclusão do comportamento das longarinas. Conhecendo com precisão a capacidade estrutural da ponte e a sua vida útil, pode-se assegurar o tráfego de composições com maior segurança.

Projeto de um monopolo Hi para aplicações no sistema brasileiro de TV digital

A presente pesquisa trata o projeto e análise de uma antena monopolo planar com geometria modificada visando sua utilização para recepção do sinal de TV digital operante no Brasil na faixa de 470 MHz a 806 MHz. Faixa essa contida no espectro de UHF – Ultra High Frequency (300 MHz a 3 GHz). Para desenvolvimento desse trabalho foi tomado como referência à antena denominada “The Hi Monopole”. Que originalmente foi apresentada para operar em sistemas UWB (Ultra Wide Band) em 3,1 a 10,6 GHz. Para o desenvolvimento do trabalho proposto, diferentes técnicas de adequação da antena podem ser utilizadas para operação em banda larga, tais como: modificação na estrutura da antena, carregamento resistivo, chaveamento, utilização de elementos parasitas e estruturas de casamento. O projeto de antenas banda larga pode ser realizado a partir de três abordagens diferentes: domínio do tempo, domínio da frequência e método de expansão por singularidades. O método no domínio da frequência foi empregado neste trabalho para o projeto da antena proposta, algumas das técnicas supracitadas foram analisadas almejando o aumento da largura de banda, sendo confeccionado um protótipo da antena para validar os conceitos empregados. A antena foi então projetada para a faixa de 470 MHz a 890 MHz. O protótipo construído para essa mesma faixa apresentou bons resultados, o que valida à técnica empregada. Aspectos positivos e negativos do uso desta técnica são discutidos ao longo do trabalho. O programa computacional comercial CST® MICROWAVE STUDIO, baseado na Técnica da Integração Finita (FIT), foi usado para simulações no domínio da frequência.

Influência do tipo de conexão protética, do platform-switching e do designe do implante no ambiente biomecânico de implantes imediatos com carga imediata

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Avaliação de protocolos para colagem de bráquetes cerâmicos e remoção do adesivo remanescente com vistas á preservação da integridade de esmalte dentário

Pós-graduação em Ciências Odontológicas - FOAR

Construção de um modelo de dispersão de poluentes na Baía do Guajará

Os ambientes marinho e costeiro do Brasil vêm sofrendo nos últimos anos um considerável processo de degradação ambiental, gerado pela crescente pressão sobre os recursos naturais marinhos e continentais e pela capacidade limitada desses ecossistemas absorverem os impactos resultantes. No caso específico da Baía do Guajará, situada no Estuário Guajará, na região delimitada pela Baía do Marajó e pelos rios Pará e Guamá, estes riscos estão associados principalmente à poluição oriunda da cidade de Belém e às operações de transporte de derivados de petróleo. No caso da Baía do Guajará, sabe-se que ocorre o despejo de esgoto periodicamente, de acordo com a maré. Desta forma, pouco tem sido feito para amenizar os efeitos do problema. Os modelos numéricos, por considerarem a complexidade dos fenômenos físicos que governam o funcionamento dos ambientes costeiros, vêm sendo utilizados como ferramentas em simulação hidrodinâmica, com reconhecida importância no gerenciamento costeiro. Sua aplicação pode dar suporte a estudos diversos, como no entendimento dos processos dispersivos de poluentes, implantação de sistemas de monitoramento de qualidade de água, planejamento de ações em casos de derrame de óleo, e avaliar as conseqüências de alterações na geometria de canais e baías e estuários. Este trabalho tem como objetivo esclarecer essas questões e descreve as etapas necessárias para a construção de um modelo numérico de dispersão na Baía do Guajará usando as ferramentas de Sistemas de Informação Geográfica e o Método de Elementos Finitos.

Avaliação de métodos eletromagnéticos aplicando campos polarizados e focalizados

Em investigações geofísicas rasas que empregam os métodos eletromagnéticos indutivos mais avançados, alvos com baixo número de indução (Low Induction Number – LIN) produzem anomalias eletromagnéticas muito baixas e de difícil interpretação. Para suprir esta deficiência, neste trabalho são estudados a aplicabilidade de campos eletromagnéticos polarizados e focalizados – POLFOCEM como fonte primária de indução. Os campos E.M. focalizados e polarizados, vertical e horizontalmente, são obtidos pelas combinações vetoriais de pares de dipolos transmissores e, ocorrem na região central entre eles. A focalização é observada nesta região na profundidade de 0,25 do espaçamento entre esses transmissores – L. Portanto, máximos acoplamentos podem ser obtidos através da seleção da polarização de acordo com a geometria do alvo, ocorrendo um aumento na densidade de fluxo magnético sobre ele e, máximas anomalias produzidas. É utilizada uma metodologia numérica para o cômputo dessas anomalias por meio da técnica dos elementos finitos para solução do problema 2,5-D. Em todos os experimentos numéricos são realizadas comparações qualitativas e quantitativas entre as respostas obtidas pelos sistemas POLFOCEM e convencional, o qual emprega um único dipolo como transmissor (dipolo-dipolo). As anomalias produzidas pelo sistema POLFOCEM, em que os dipolos transmissores são acionados simultaneamente, correspondem à soma das anomalias produzidas por cada um desses dipolos independentes, caracterizando, desta forma, a linearidade dos campos eletromagnéticos. Os experimentos numéricos são realizados para alvos prismáticos bidimensionais com três diferentes inclinações, inseridos num semi-espaço resistivo, e para as freqüências das fontes na faixa das ondas de rádio. As anomalias assimétricas no sistema convencional, que se tornam simétricas no sistema POLFOCEM, apresentam valores menores em amplitude. Contudo, aquelas anomalias tanto assimétricas quanto simétricas que se tornam anti-simétricas apresentam valores maiores. Em decorrência dessas diminuições e aumentos nas amplitudes ocorrem rotações nos diagramas de Argand, no sentido horário e anti-horário para alvos com baixos valores de condutividade, respectivamente. Em experimentos de identificação de presença de dois alvos próximos, o sistema convencional é capaz de identificá-los primeiramente, prevalecendo o seu uso.

Modelagem numérica de dados MCSEM 3D usando computação paralela

Desenvolvemos a modelagem numérica de dados sintéticos Marine Controlled Source Electromagnetic (MCSEM) usada na exploração de hidrocarbonetos para simples modelos tridimensionais usando computação paralela. Os modelos são constituidos de duas camadas estrati cadas: o mar e o sedimentos encaixantes de um delgado reservatório tridimensional, sobrepostas pelo semi-espaço correspondente ao ar. Neste Trabalho apresentamos uma abordagem tridimensional da técnica dos elementos nitos aplicada ao método MCSEM, usando a formulação da decomposição primária e secundária dos potenciais acoplados magnético e elétrico. Num pós-processamento, os campos eletromagnéticos são calculados a partir dos potenciais espalhados via diferenciação numérica. Exploramos o paralelismo dos dados MCSEM 3D em um levantamento multitransmissor, em que para cada posição do transmissor temos o mesmo processo de cálculos com dados diferentes. Para isso, usamos a biblioteca Message Passing Interface (MPI) e o modelo servidor cliente, onde o processador administrador envia os dados de entradas para os processadores clientes computar a modelagem. Os dados de entrada são formados pelos parâmetros da malha de elementos nitos, dos transmissores e do modelo geoelétrico do reservatório. Esse possui geometria prismática que representa lentes de reservatórios de hidrocarbonetos em águas profundas. Observamos que quando a largura e o comprimento horizontais desses reservatório têm a mesma ordem de grandeza, as resposta in-line são muito semelhantes e conseqüentemente o efeito tridimensional não é detectado. Por sua vez, quando a diferença nos tamanhos da largura e do comprimento do reservatório é signi cativa o efeito 3D é facilmente detectado em medidas in-line na maior dimensão horizontal do reservatório. Para medidas na menor dimensão esse efeito não é detectável, pois, nesse caso o modelo 3D se aproxima de um modelo bidimensional. O paralelismo dos dados é de rápida implementação e processamento. O tempo de execução para a modelagem multitransmissor em ambiente paralelo é equivalente ao tempo de processamento da modelagem para um único transmissor em uma máquina seqüêncial, com o acréscimo do tempo de latência na transmissão de dados entre os nós do cluster, o que justi ca o uso desta metodologia na modelagem e interpretação de dados MCSEM. Devido a reduzida memória (2 Gbytes) em cada processador do cluster do departamento de geofísica da UFPA, apenas modelos muito simples foram executados.

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).

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.

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 teórica do problema de Weaver da falha infinita, modo TE magnetotelúrico

É apresentada uma solução totalmente analítica do modelo da falha infinita para o modo TE magnetotelúrico, levando em conta a presença do ar, utilizando como base o trabalho de Sampaio apresentado em 1985, que apresenta uma solução parcialmente analítica e parcialmente numérica – solução híbrida. Naquela solução foram aplicadas oito condições de contorno, sendo que em quatro delas foram encontradas inconsistências matemáticas que foram dirimidas com alterações adequadas nas soluções propostas por Sampaio. Tais alterações propiciaram a chegarse à solução totalmente analítica aqui apresentada. A solução obtida foi comparada com a solução de Weaver, com a de Sampaio e com o resultado do método numérico dos elementos finitos para contrastes de resistividade iguais a 2, 10 e 50. A comparação da solução analítica, para o campo elétrico normalizado, com a solução de elementos finitos mostra que a solução analítica proporcionou resultados mais próximos, em comparação aos fornecidos por Weaver e por Sampaio. Este é um problema muito difícil, aberto para uma solução analítica definitiva. A solução apresentada aqui é, nesta direção, um grande passo.

Efeito da topografia em sondagens eletromagnéticas bidimensionais nos domínios da frequência e do tempo

O efeito da topografia de um vale, uma colina, um declive e um aclive, em sondagens eletromagnéticas bidimensionais no domínio da freqüência, causa uma variação nos valores da amplitude e da fase da componente H_z devido a um meio homogêneo. A amplitude é menos afetada que a fase. A parede do vale mais próxima da linha de corrente causa uma forte diminuição dos valores da fase, enquanto que a parede do vale mais distante da linha causa um forte aumento. Os efeitos de uma colina são opostos aos do vale. Os efeitos do declive e do aclive num meio homogêneo, são similares, respectivamente, aos observados pelas paredes do vale e da colina mais próximas da linha de corrente. A resposta de um corpo condutivo retangular num meio homogêneo próximo a uma linha de corrente sofre pequenas variações devido a presença de um vale ou de uma colina situada longe da linha de corrente. Porém, se essas feições topográficas estiverem sobre o corpo, elas afetam fortemente a fase e a amplitude da componente H_z e apenas a amplitude, no caso da componente H_x. A resposta transiente não é afetada pela topografia para tempos muito baixos, pois nesse caso se investiga uma finíssima camada da superfície, assim como para tempos muito altos porque a profundidade de investigação é muito grande comparada com a dimensão da topografia. Para os modelos aqui estudados, a maior influência se dá para tempos intermediários, ao redor de 7 ms, ocasionando um retardo do ponto de "cross over" nas curvas de sondagens.

Modelagens direta e inversa de dados mCSEM 3D

Marine Controlled Source Electromagnetic - mCSEM é um método geofísico eletromagnético que nos últimos dez anos vem sendo usado na prospecção de hidrocarbonetos com bastante êxito. Este método consiste em um dipolo elétrico horizontal (DEH) localizado um pouco acima do assoalho marinho, operando em baixa frequência (0,1-1,0 Hz) e receptores regularmente distribuídos no fundo do mar que captam os campos eletromagnéticos provenientes da difusão de energia gerada pelo dipolo transmissor. Neste trabalho vamos apresentar o problema direto do método mCSEM 3D, propondo soluções numéricas, através do método dos elementos finitos tridimensionais, para modelos geoelétricos mCSEM 3D. Para fins de análise de coerência, os resultados obtidos são comparados com soluções disponíveis na literatura. Em seguida, apresentaremos a inversão de um de seus modelos segundo uma proposta de metodologia de inversão juntamente com a proposta de solução direta para o mCSEM 3D, acima mencionada, realizando assim a inversão de um modelo geoelétrico do mCSEM 3D para duas frequências.