Modelagem 2,5D por elementos finitos dos efeitos da topografia do terreno sobre dados obtidos com o método eletromagnético a multifrequência

Este trabalho consiste em realizar a modelagem, via elementos finitos (EF) 2,5D, do efeito da topografia do terreno sobre dados obtidos com o método eletromagnético a multi-frequência (EMMF). Este método usa como fonte uma grande espira quadrada de corrente elétrica com centenas de metros de lado, e como receptores, bobinas posicionadas na horizontal em alinhamento com o transmissor. A subsuperfície é representada por heterogeneidades bidimensionais imersas em um meio horizontalmente estratificado. A formulação, partindo das equações de Maxwell, é desenvolvida a partir da separação do campo eletromagnético em primário (campos no hospedeiro multi-estratificado) e secundário (diferença entre o campo total e o primário). O domínio discretizado é descrito por uma malha não estruturada, com elementos triangulares. Para calcular as componentes derivadas da solução de elementos finitos, em um determinado nó da malha, foi usada a média aritmética das derivadas das funções bases de EF em torno daquele nó. O código de modelagem construído permite quantificar e analisar como os gradientes topográficos influenciam as medidas dos campos eletromagnéticos gerados. A aplicação é a avaliação dessas influências sobre a componente radial do campo da espira na superfície terrestre, que é a componente empregada no método eletromagnético a multi-frequência (EMMF).

Modelagem de transformadores de núcleo amorfo usando método de elementos finitos

Os materiais amorfos vêm ganhando grande espaço na indústria de transformadores, devido às baixas perdas no núcleo, visto que estes possuem um ciclo de histerese mais estreito quando comparados com os núcleos tradicionais de aço silício. Entretanto, seu custo ainda tem sido um grande fator para pouca inserção deste tipo de equipamento nos sistemas elétricos de potência. Estudos sobre o custo/benefício em termos técnicos de desempenho e robustez devem ser considerados quando se deseja projetar transformadores que utilizam esse tipo de material. A grande contribuição deste trabalho encontra-se na análise do desempenho do transformador de núcleo amorfo diante de curtos-circuitos, visto que tais falhas neste equipamento ocasionam redução da receita, não apenas por gastos com manutenção, mas também porque a concessionária deixa de vender seu produto, energia elétrica, além de poder estar sujeita a penalidades por parte dos órgãos de regulação do setor elétrico. Quando em condições de curto-circuito os enrolamentos dos transformadores ficam submetidos a esforços mecânicos, produzidos por forças de Lorentz, essas forças surgem como resultado do fluxo produzido pelos próprios condutores em paralelo que transportam corrente na mesma direção. Diante disso, estudar o comportamento eletromagnético do transformador é fundamental para obtenção de tais forças. Para o desenvolvimento deste trabalho, utilizou-se o software Finite Element Method Magnetics (FEMM). Esta ferramenta se baseia no método de elementos finitos para realizar os cálculos das magnitudes eletromecânicas e, consequentemente, o cálculo das forças atuando nas espiras, as quais permitem realizar os cálculos dos esforços mecânicos. Por fim, este trabalho aborda a aplicação da ferramenta FEMM para o cálculo de esforços mecânicos e simulação do comportamento eletromagnético de um transformador de distribuição.

Estudo comparativo das técnicas de elementos finitos e equação integral na modelagem eletromagnética bidimensional

Os métodos numéricos de Elementos Finitos e Equação Integral são comumente utilizados para investigações eletromagnéticas na Geofísica, e, para essas modelagens é importante saber qual algoritmo é mais rápido num certo modelo geofísico. Neste trabalho são feitas comparações nos resultados de tempo computacional desses dois métodos em modelos bidimensionais com heterogeneidades condutivas num semiespaço resistivo energizados por uma linha infinita de corrente (com 1000Hz de freqüência) e situada na superfície paralelamente ao "strike" das heterogeneidades. Após a validação e otimização dos programas analisamos o comportamento dos tempos de processamento nos modelos de corpos retangulares variandose o tamanho, o número e a inclinação dos corpos. Além disso, investigamos nesses métodos as etapas que demandam maior custo computacional. Em nossos modelos, o método de Elementos Finitos foi mais vantajoso que o de Equação Integral, com exceção na situação de corpos com baixa condutividade ou com geometria inclinada.

Modelagem acústica por diferenças finitas e elementos finitos em 2-D e 2,5-D

A modelagem acústica fornece dados úteis para avaliação de metodologias de processamento e imageamento sísmico, em modelos com estrutura geológica complexa. Esquemas de diferenças finitas (DF) e elementos finitos (EF) foram implementados e avaliados em modelos homogêneos e heterogêneos. O algoritmo de diferenças finitas foi estendido para o caso 2,5-D em modelos com densidade variável. Foi apresentada a modelagem de alvos geológicos de interesse exploratório existentes na Bacia Paleozóica do Solimões na Amazônia. Reflexões múltiplas de longo período produzidas entre a superfície livre e a discordância Cretáceo-Paleozóica, a baixa resolução da onda sísmica nas proximidades do reservatório e as fracas reflexões na interface entre as rochas reservatório e as rochas selantes são as principais características dos dados sintéticos obtidos, os quais representam um grande desafio ao imageamento sísmico.

Modelagem direta bidimensional do método magnetotelúrico com o método dos elementos finitos de arestas

Fizemos a modelagem direta 2D do método magnetotelúrico (MT) com o método dos elementos finitos (MEF) de arestas em termos dos campos primários e secundários. Para usarmos modelos de maior complexidade e diminuirmos o custo computacional utilizamos malhas não estruturadas. Nas malhas utilizadas, introduzimos quatro nós em torno de cada estação MT, constituindo um quadrado alinhado nas direções dos eixos cartesianos

Os efeitos Delaware e Groningen: um estudo quantitativo por elementos finitos

Os efeitos Delaware e Groningen são dois tipos de anomalia que afetam ferramentas de eletrodos para perfilagem de resistividade. Ambos os efeitos ocorrem quando há uma camada muito resistiva, como anidrita ou halita, acima do(s) reservatório(s), produzindo um gradiente de resistividade muito similar ao produzido por um contato óleo-água. Os erros de interpretação produzidos têm ocasionado prejuízos consideráveis à indústria de petróleo. A PETROBRÁS, em particular, tem enfrentado problemas ocasionados pelo efeito Groningen sobre perfis obtidos em bacias paleozóicas da região norte do Brasil. Neste trabalho adaptamos, com avanços, uma metodologia desenvolvida por LOVELL (1990), baseada na equação de Helmholtz para H_Φ, para modelagem dos efeitos Delaware e Groningen. Solucionamos esta equação por elementos finitos triangulares e retangulares. O sistema linear gerado pelo método de elementos finitos é resolvido por gradiente bi-conjugado pré-condicionado, sendo este pré-condicionador obtido por decomposição LU (Low Up) da matriz de stiffness. As voltagens são calculadas por um algoritmo, mais preciso, recentemente desenvolvido. Os perfis são gerados por um novo algoritmo envolvendo uma sucessiva troca de resistividade de subdomínios. Este procedimento permite obter cada nova matriz de stiffness a partir da anterior pelo cálculo, muito mais rápido, da variação dessa matriz. Este método permite ainda, acelerar a solução iterativa pelo uso da solução na posição anterior da ferramenta. Finalmente geramos perfis sintéticos afetados por cada um dos efeitos para um modelo da ferramenta Dual Laterolog.

Análise termoelástica experimental e computacional por elementos finitos da parede externa em concreto armado de um forno de cozimento de anodo

Nesta dissertação foi desenvolvido um estudo sobre o comportamento de uma estrutura em concreto armado submetida a elevadas variações térmicas. A estrutura analisada consiste na parede externa de um forno de cozimento de anodo utilizado na produção de Alumínio, e está submetida a variações térmicas de cerca de 125C devido à própria operação do forno. As motivações principais para o estudo foram a ocorrência de grandes deformações e o surgimento de fissuras na estrutura, o que poderia provocar a inutilização do forno. O objetivo do trabalho foi a investigação das causas destas patologias, assim como o estabelecimento de propostas para reforço da estrutura. Para isso, foram realizados estudos experimentais e computacionais do comportamento da estrutura. Inicialmente, a estrutura foi monitorada utilizando-se transdutores de deslocamento e sensores de temperatura (termopares), conectados a um sistema de aquisição de dados para obtenção e armazenamento automática das amostras ao longo do tempo. Em seguida, foram desenvolvidos modelos computacionais em Elementos Finitos com auxilio do programa computacional Algor, para determinação da distribuição de temperatura e as correspondentes tensões e deformações de origem térmica na estrutura. Nestes estudos, foram realizadas análises estacionárias e transientes de condução de calor, seguidas de análises de tensões de origem térmica. Como conclusão do estudo, tem-se que a metodologia proposta para a solução do problema foi bastante satisfatória, solucionando o problema de forma precisa, porém econômica.

Modelagem 1D e 2,5D de dados do método CSEM marinho em meios com anisotropia transversal inclinada

Neste trabalho apresentamos a solução do campo eletromagnético gerado por um dipolo elétrico horizontal em meios transversalmente isotrópicos com eixo de simetria vertical (TIV) e com eixo de simetria inclinado (TII). Para modelos unidimensionais, o campo eletromagnético foi obtido por duas metodologias distintas: (1) solução semi-analítica das equações de Maxwell com auxílio de potenciais vetores no caso TIV e (2) em modelos com anisotropia transversal inclinada o campo eletromagnético foi separado em primário e secundário, e então, o campo secundário foi calculado pelo método de elementos finitos no domínio (k_x, k_y, z) da transformada de Fourier. Para estruturas bidimensionais, foi aplicada a mesma metodologia usado nos modelos TII unidimensionais, onde o campo secundário foi calculado pelo método de elementos finitos no domínio (x, k_y, z), da transformada de Fourier, com a utilização de malhas não estruturadas para discretização dos modelos. Estas respostas foram usados para avaliar os efeitos da anisotropia elétrica nos dados CSEM marinho 1D e 2,5D.

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.

Modelagem de dados MT 2D em multicamadas com anisotropia arbitrária

Neste trabalho apresentamos a modelagem direta bidimensional do método Magnetotelúrico (MT) com anisotropia arbitrária na condutividade elétrica para modelos de multicamadas estratificadas horizontalmente, separando a resposta do campo resultante em duas partes, primária e secundária. A parte primária refere-se ao campo 1D do modelo na ausência de heterogeneidades; a secundária, à contribuição no campo resultante gerada pela anomalia. Esta modelagem foi feita via técnica numérica dos Elementos Finitos (EF). Apresentamos também a modelagem direta do caso MT 1D para o mesmo modelo e anisotropia arbitrária, porém, com solução analítica pelo método da matriz de propagação; tal resultado é usado como fonte na formulação 2D. No estudo a respeito do comportamento dos campos 1D em profundidade na presença da anisotropia, verificamos nas componentes horizontais que suas curvas dentro de cada camada se assemelham a de um campo na presença de fontes, pois experimentam um aumento na amplitude; além de verificar o surgimento da componente normal E_z do campo, sem que haja uma densidade de corrente nesta direção. Com arbitrariedade na anisotropia queremos dizer que seus eixos principais podem assumir qualquer orientação em relação ao eixo de referência adotado para o problema; tal orientação é descrita sob três rotações consecutivas, chamadas de strike, dip e slant, que simulam inclinações bastante presentes na geologia dos ambientes terrestres. Verificamos que a anisotropia provoca alterações no comportamento das pseudo-seções de resistividades aparentes, interferindo na interpretação correta dos dados.

Modelagem 2,5D dos campos usados no Método Eletromagnético a Multi-Frequência - EMMF

Esta tese mostra a modelagem 2,5D de dados sintéticos do Método Eletromagnético a Multi-frequência (EMMF). O trabalho é apresentado em duas partes: a primeira apresenta os detalhes dos métodos usados nos cálculos dos campos gerados por uma bobina horizontal de corrente colocada sobre a superfície de modelos bidimensionais; e a segunda, usa os resultados obtidos para simular os dados medidos no método EMMF, que são as partes real e imaginária da componente radial do campo magnético gerado pela bobina. Nesta segunda parte, observamos o comportamento do campo calculado em diversos modelos, incluindo variações nas propriedades físicas e na geometria dos mesmos, com o intuito de verificar a sensibilidade do campo observado com relação às estruturas presentes em uma bacia sedimentar. Com esta modelagem, podemos observar as características dos dados e como as duas partes, real e imaginária, contribuem com informações distintas e complementares. Os resultados mostram que os dados da componente radial do campo magnético apresentam muito boa resolução lateral, mesmo estando a fonte fixa em uma única posição. A capacidade desses dados em distinguir e resolver estruturas alvo será fundamental para o trabalho futuro de inversão, bem como para a construção de seções de resistividade aparente.

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.

Ensaios sobre a contribuição do ar nos resultados da modelagem do modo transverso-elétrico do método magnetotelúrico

Apresentamos a solução numérica do modelo da falha infinita para o modo TE utilizando o método dos elementos finitos para investigar o efeito do ar na solução TE. A comparação de nossa solução para o campo elétrico normalizado com a solução de Weaver evidenciou uma discrepância esperada devido ao efeito do ar, que foi intencionalmente negligenciado por Weaver, com o intento de facilitar a solução do problema. O problema analítico, levando em conta a presença do ar, é substancialmente mais difícil. Uma solução do problema real, parcialmente analítica e parcialmente numérica, foi apresentada por Sampaio. Devido a uma discrepância inesperada entre a nossa solução por elementos finitos e a solução exata de Sampaio, recalculamos a sua solução. Iterando a parte numérica da solução de Sampaio até ordens mais altas pudemos verificar que esta apresenta comportamento divergente. Para examinar o efeito do ar na resistividade aparente, comparamos três perfis de resistividade aparente obtidos a partir de nossa solução por elementos finitos com os mesmos perfis obtidos pela solução de Weaver. A discrepância encontrada entre estas soluções foi muito pequena, situando-se na faixa de erro instrumental de campo.

Estudos dos efeitos do manto de intemperismo sobre anomalias de VLF geradas por corpos inclinados: modelagem numérica através de programa eGs

As referências apresentadas e utilizadas através de métodos indutivos de propagação eletromagnética E.M., são utilizadas na prospecção de minerais, entre outros, haja vista que em regiões tropicais existem camadas superficiais condutivas (manto intemperismo) que levam a influenciar, evidentemente, na interpretação do resultado real. Em nossa região Amazônica, o manto apresenta-se bastante desenvolvido e condutivo. O objetivo deste trabalho é estudar os efeitos dos mantos sobre anomalias VLF, devidos a corpos tabulares inclinados em contato ou não com o manto (manto galvânico e manto indutivo). Deste modo, foram realizadas umas séries de experimentos em modelagens numéricas através da eGs, considerando-se diferentes parâmetros de resposta para o manto e para o corpo. Os trabalhos de modelagem numérica, bem como seus resultados foram apresentados por diversas situações, como: variação da profundidade do topo do condutor, variação do mergulho do condutor, variação da condutância do manto, entre outros. Para interpretação dos resultados, foi gerado um conjunto de curvas considerando-se valores pico-a-pico das anomalias de Tilt Angle e de Elipsidade. Os resultados foram sintetizados em um outro conjunto de curvas reunidas em diagramas de Argand. A medida em que se aumenta a profundidade, o pico-a-pico e a forma do pico da anomalia do Tilt Angle e da Elipsidade diminuem na presença do manto, ou na ausência do mesmo, fazendo com que este se afaste do ponto de “Cross-Over”. Para um condutor da mesma espessura, o aumento da condutância causa um aumento nas anomalias de Tilt Angle e da Elipsidade. O efeito geral na variação do mergulho do condutor é o de causar uma assimetria nos perfis de Tilt Angle e Elipsidade. O aumento da condutância do manto de intemperismo causa uma diminuição nas anomalias de Tilt Angle e de Elipsidade. Há rotação de fase no sentido anti-horário, sendo mais intensa para grandes valores de condutância do corpo. Na presença do manto, o corpo parece estar a uma profundidade inferior à verdadeira e, tendendo, a ser menos condutivo.

Modelagem computacional de dados magnetotelúricos marinhos 2-D

O presente trabalho tem como objetivo investigar o espalhamento das ondas planas causados pela variação lateral das propriedades físicas das rochas. Essa é uma das questões mais importante para o sucesso da exploração geofísica. O método geofísico usado nesta dissertação foi método magnetotelúrico (MT) e método Magnetotelúrico marinho (mMT). Umas das ferramentas utilizadas foram o método dos elementos finitos, este método é notável pela flexibilidade em resolver numericamente as equações diferenciais para campos eletromagnéticos secundários, gerados por estruturas geológicas com geometrias complexas, sendo utilizado para obter a solução numérica do espalhamento do campo eletromagnético da onda plana por um dique vertical (modelos bidimensionais). Os algoritmos desenvolvidos foram aplicados em diversos modelos de dique vertical, com diferentes parâmetros geoelétricos.