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.

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 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 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 eletrorresistividade 2-D a partir do potencial elétrico secundário

Um dos métodos clássicos da geofísica de exploração é o Método de Eletrorresistividade, estabelecido há um século pelos irmãos Schlumberger e desde então amplamente empregado em prospecção mineral, estudos ambientais e hidrogeologia e em pesquisa de fontes geotermais. Conceitualmente o método consiste de injeção de corrente elétrica na subsuperfície e de medida de diferença de potencial elétrico, resultante da interação da corrente com o meio. As localizações dessas fontes e receptores são determinadas pelo arranjo escolhido para o levantamento. Após o processamento, obtém-se pseudo-seções de resistividade aparente que indicam a distribuição de condutividade em subsuperfície. Devido à simplicidade dos fundamentos físicos de sua formulação, o método apresenta fácil implementação computacional quando comparado aos métodos eletromagnéticos de fonte controlada. Na literatura há inúmeros trabalhos de modelagem computacional, onde se calcula a resposta para problemas 2-D e 3-D. Nestes trabalhos, as pseudo-seções são obtidas a partir do cálculo do potencial elétrico total. Neste trabalho, apresentaremos a resposta da modelagem de eletroresistividade 2-D com o arranjo dipolo-dipolo, obtida a partir do potencial elétrico secundário. A solução é calculada através do método de elementos finitos usando malhas não estruturadas. Para efeito de validação, os resultados são comparados com a resposta 2-D obtida a partir dos potencias totais.

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.

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

Aplicação dos métodos geofísicos elétrico e eletromagnético na determinação de unidades sedimentares costeiras tropicais em Bragança, Nordeste do Pará

Esta dissertação mostra a aplicação de métodos geofísicos na determinação das unidades sedimentares da Planície Costeira Bragantina, que estão sendo estudados no Programa de Manejo e Dinâmica de Manguezais-MADAM (Mangrove Dynamics and Management). A finalidade deste trabalho é testar a metodologia geofísica no ambiente costeiro, visando auxiliar o mapeamento geológico e o entendimento da evolução estratigráfica holocênica da Planície Costeira Bragantina. Na realização deste trabalho, foram aplicados em diferentes unidades morfoestratigráficas da área de estudo: (a) o método convencional da eletrorresistividade através de sondagens elétricas verticais, usando o arranjo Schlumberger e (b) o método eletromagnético slingram (Horizontal Loop Eletromagnetic), com oito freqüências. Com as sondagens elétricas verticais determinaram-se os valores de resistividade em subsuperfície para os diferentes horizontes, correlacionando-os sempre que possível com as informações geológicas disponíveis. Em geral, os resultados obtidos com as sondagens elétricas não foram muito satisfatórios em termos das profundidades investigadas pelo fato do ambiente estudado ser muito condutivo, dificultando a penetração de corrente elétrica e contribuindo para que a profundidade de investigação fosse baixa. A aplicação do método eletromagnético mostrou-se mais simples e rápido do que as sondagens elétricas. A interpretação dos dados obtidos com esse método, realizada através da análise de perfis de medidas e de sondagens eletromagnéticas, permitiu inferir contatos laterais e estabelecer assinaturas geofísicas para os depósitos da planície costeira.

Comparação do desempenho dos métodos eletroresistivo e transiente eletromagnético num caso específico de prospecção de águas subterrâneas

O método eletroresistivo é um dos métodos de prospecção geofísica mais empregados na exploração de água subterrânea. Todavia, ele apresenta a desvantagem de ser de difícil uso em locais habitados. De fato, para realizar sondagens elétricas é necessário, muitas vezes, expandir cabos elétricos até um quilômetro ou mais, dependendo da profundidade que se deseja investigar, o que nem sempre é possível em áreas urbanas. O método transiente eletromagnético no domínio do tempo, é usado também para a exploração de água subterrânea e não apresenta as mesmas desvantagens, porque este método indutivo emprega uma bobina quadrada cujo tamanho não é um impedimento para ser utilizado em locais habitados. Por outro lado, o método eletromagnético é mais sensível a ruídos devido a instalações elétricas. Para comparar as vantagens e desvantagens dos dois métodos em áreas urbanas foram utilizados programas para a modelagem sintética das sondagens elétricas verticais e das sondagens transientes eletromagnéticas, ambos desenvolvidos por RIJO (1994ab, 1996). Para tal propósito foram considerados dois casos: no primeiro, foram definidos quatro tipos de modelos (H, K, A e Q) com base nos horizontes geoelétricos em Palestina do Pará. No segundo caso, se considerou os modelos obtidos pelas 21 sondagens elétricas verticais tipo Schlumberger da mesma área. Com base nestes modelos, se utilizou um programa de transientes electromagnéticas para gerar as curvas sintéticas das resistividades aparentes. Como resultado, os valores de resistividades aparentes dados em função do espaçamento de eletrodos (método eletroresistivo) e em função do tempo (método eletromagnético), permitiu obter nos dois casos mencionados uma relação linear em escala bilogaritmica entre o espaçamento (metros) e o tempo (milisegundos). Deste modo é possível inferir a profundidade de investigação desejada em áreas urbanas usando o método eletromagnético, quando a aplicação do método eletroresistivo ficar limitada pela acessibilidade.

Filtros otimizados para transformadas seno, co-seno e de Hankel j0, j1 e j2

Apresentamos dois algoritmos automáticos, os quais se utilizam do método dos mínimos quadrados de Wiener-Hopf, para o cálculo de filtros lineares digitais para as transformadas seno, co-seno e de Hankel J₀, J₁ e J₂. O primeiro, que otimiza os parâmetros: incremento das abscissas, abscissa inicial e o fator de deslocamento utilizados para os cálculos dos coeficientes dos filtros lineares digitais que são aferidos através de transformadas co-seno, seno e o segundo, que otimiza os parâmetros: incremento das abscissas e abscissa inicial utilizados para os cálculos dos coeficientes dos filtros lineares digitais que são aferidos através de transformadas de Hankel J₀, J₁ e J₂. Esses algoritmos levaram às propostas de novos filtros lineares digitais de 19, 30 e 40 pontos para as transformadas co-seno e seno e de novos filtros otimizados de 37 , 27 e 19 pontos para as transformadas J₀, J₁ e J₂, respectivamente. O desempenho dos novos filtros em relação aos filtros existentes na literatura geofísica é avaliado usando-se um modelo geofísico constituído por dois semi-espaços. Como fonte usou-se uma linha infinita de corrente entre os semi-espaços originando, desta forma, transformadas co-seno e seno. Verificou-se melhores desempenhos na maioria das simulações usando o novo filtro co-seno de 19 pontos em relação às simulações usando o filtro co-seno de 19 pontos existente na literatura. Verificou-se também a equivalência de desempenhos nas simulações usando o novo filtro seno de 19 pontos em relação às simulações usando o filtro seno de 20 pontos existente na literatura. Adicionalmente usou-se também como fonte um dipolo magnético vertical entre os semi-espaços originando desta forma, transformadas J₀ e J₁, verificando-se melhores desempenhos na maioria das simulações usando o novo filtro J₁ de 27 pontos em relação ao filtro J₁ de 47 pontos existente na literatura. Verificou-se também a equivalência de desempenhos na maioria das simulações usando o novo filtro J₀ de 37 pontos em relação ao filtro J₀ de 61 pontos existente na literatura. Usou-se também como fonte um dipolo magnético horizontal entre os semi-espaços, verificando-se um desempenho análogo ao que foi descrito anteriormente dos novos filtros de 37 e 27 pontos para as respectivas transformadas J₀ e J₁ em relação aos filtros de 61 e 47 pontos existentes na literatura, destas respectivas transformadas. Finalmente verificou-se a equivalência de desempenhos entre os novos filtros J₀ de 37 pontos e J₁ de 27 pontos em relação aos filtros de 61 e 47 pontos existentes na literatura destas transformadas, respectivamente, quando aplicados em modelos de sondagens elétricas verticais (Wenner e Schlumberger). A maioria dos nossos filtros contêm poucos coeficientes quando comparados àqueles geralmente usados na geofísica. Este aspecto é muito importante porque transformadas utilizando filtros lineares digitais são usadas maciçamente em problemas numéricos geofísicos.

Métodos sísmicos de refração e eletromagnético aplicados ao estudo de um paleocanal da Ilha do Marajó, Pará

Os métodos sísmico de refração e eletromagnético foram aplicados, em caráter experimental, à prospecção de água subterrânea. A eficiência destes métodos foi testada com cerca de 1.500m de perfis de refração rasa na delimitação geométrica em três dimensões de um paleocanal na região dos campos naturais da Ilha de Marajó e com um perfil eletromagnético (técnica Slingram) para a detecção de contatos laterais de água doce-água salobra. A interpretação dos dados sísmicos foi realizada com as equações de tempo de intersecção para n camadas inclinadas com o auxílio do controle geológico de 2 (dois) poços para contornar problemas de camadas não detectadas. Obteve-se uma baixa velocidade (0,81 km/s) para os sedimentos silte-arenosos saturados que preenchem o paleocanal. Em virtude da inversão de velocidades ter sido obtida como resultado para estes sedimentos, não se aconselha o mapeamento de paleocanais pela sísmica de refração, salvo com a utilização de técnicas especiais de campo (registro digital, 'up-hole survey', registro ondas S, reflexão simultânea, ...). Os dados eletromagnéticos foram analisados semi-quantitativamente com o auxílio de um modelo reduzido construído a partir de uma seção de resistividades aparentes. Um diagrama de Argand, traçado com os dados de laboratório, fornece subsídios para a determinação da condutância em ambientes de sedimentos não consolidados análogos ao de Marajó. Recomenda-se a utilização dos métodos eletromagnéticos para o reconhecimento lateral de zonas de água doce-água-salobra, uma vez que a componente quadratura da frequência alta (2640 Hz) define, com precisão suficiente, estas zonas. A comparação dos resultados dos dois métodos mostra que apenas parte do paleocanal se comporta como um aquífero.

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.

Modelagem unidimensional de dados do método MCSEM com variação na resistividade da água do mar e na orientação dos sensores

O Marine Controlled Source ElectroMagnetic (MCSEM) é um novo método geofísico que tem sido empregado para a detecção de camadas resistivas no sub-solo oceânico e tem se mostrado uma ferramenta promissora para a indústria do petróleo. Neste trabalho, é modelado levantamentos do método MCSEM incluindo variações nos valores de resistividade da água do mar. Tais variações podem resultar da influência de correntes marinhas, gradientes de temperatura ou de qualquer fator que possa afetar a salinidade da água. Observamos que os dados do método são afetados de maneira significativa por tais variações. Tanto a variação na resistividade quanto na espessura da camada de água pode afetar os resultados do MCSEM de maneira semelhante. É modelado também o posicionamento aleatório dos receptores, mostrando um decaimento significativo da amplitude do sinal quando os receptores se posicionam a uma certa distância da linha de medidas.

Modelamento eletromagnético analógico de corpos tabulares em contato e sem contato com o manto

Nas últimas três décadas os métodos eletromagnéticos vem se desenvolvendo satisfatoriamente em função da aplicabilidade na prospecção de corpos de sulfetos maciços. Em regiões de climas tropicais, normalmente o manto de intemperismo apresenta-se condutivo. E este, na maioria das vezes, não é levado em conta na prospecção eletromagnética, causando portanto erros consideráveis de interpretações. Neste trabalho consideramos o manto de intemperismo em contato e sem-contato com o corpo condutor. Com o objetivo de estudar os efeitos dos mantos sobre anomalias EM de corpos tabulares inclinados, foram feitos vários experimentos utilizando modelamento analógico em escala reduzida, admitindo-se diferentes parâmetros de resposta para o corpo e o manto. Para simular o corpo foram utilizadas placas de aço inoxidável com as dimensões suficientemente grande em relação ao espaçamento entre as bobinas, de tal modo que simulassem um semi-plano. Para simular o manto foi utilizado uma solução de sais, sendo que para o caso de manto-condutivo, o corpo foi colocado em contato galvânico com a solução. Para o manto-indutivo foi considerado sem-contato galvânico, de tal forma que o corpo e o manto fossem acoplados apenas indutivamente. Onde 1) o corpo foi colocado totalmente sem contato com o manto 2) o corpo foi revestido por uma película resistiva e colocado em contato com o manto. Com a presença de manto-indutivo, observamos que a amplitude dos perfis é levemente atenuada. Além disso, observamos na quadratura a reversão e o aparecimento de um pico-extra nas inclinações do corpo θ≤60º e nos valores de número de indução α_c≥78.16 e α_m≥0.5 respectivamente do corpo e do manto. E a rotação de fase se dava no sentido horário, sendo mais intensa para altos valores de número de indução do corpo. No manto-indutivo o corpo parece estar a uma profundidade maior que a verdadeira, e ser mais condutivo do que realmente é. Com manto-condutivo, observamos que as amplitudes dos perfis são ligeiramente acrescidos assim como, a rotação de fase se dava no sentido anti-horário sendo mais intensa para pequenos valores de número de indução do corpo. Os demais efeitos tais como reversão na quadratura e presença de pico-extra ocorrem de modo análogo ao ocorrido no caso de manto-indutivo. No manto-condutivo, o corpo parece estar a uma profundidade inferior à verdadeira e ser menos condutivo. As anomalias EM são ligeiramente modificadas em função da rotação de fase e atenuação de amplitudes que ocorrem nos campos primário e secundário quando atravessam o manto, e também em consequência da interação indutiva de corrente induzida entre corpo e manto. Além disso ocorre a redistribuição de corrente no manto devido à presença de corpo dentro do manto. No manto-condutivo as correntes são canalizadas dentro do corpo que está em contato galvânico com o manto, enquanto no manto-indutivo ocorre um desvio de corrente, devido à película resistiva que envolve o corpo.

Inversão de dados eletromagnéticos com o regularizador Variação Total e o uso da matriz de sensibilidade aproximada

Neste trabalho apresentamos um estudo da aplicação do regularizador “Variação Total” (VT) na inversão de dados geofísicos eletromagnéticos. O regularizador VT reforça a proximidade entre os parâmetros adjacentes, mas, quando a influência de uma descontinuidade é sentida nos dados, este permite mudanças abruptas sobre os parâmetros. Isso faz com que o método seja uma alternativa válida, quando os dados observados usados na inversão provém de um ambiente geológico com uma distribuição suave de condutividade, mas que pode apresentar descontinuidades em lugares como as interfaces entre as camadas geoelétricas, como na margem de uma zona de óleo ou de um corpo de sal, que podem ser zonas muito resistivas no interior de sedimentos condutivos. Quando, devido a baixa resolução nos dados, o método não tem informações o suficiente para identificar a interface, o regularizador variação total reforça a proximidade entre os parâmetros adjacentes fazendo um transição suave entre as condutividades camadas, da mesma forma que é apresentado pela suavidade global. O método de Variação Total permite que modelos menos suaves sejam alcançados porque na norma L₁ a medida de desajuste entre os pares de parâmetros adjacentes, dará o mesmo valor se a variação dos parâmetros é suave ou se a variação é abrupta, o que não é o caso se o mesmo desajuste é medido na norma L₂, pois em uma distribuição suave a medida do desajuste é menor, sendo assim favorecida pela minimização desta norma. O uso deste regularizador permite uma melhor estimativa do tamanho de um corpo, seja ele resistivo ou condutivo. O trabalho está apresentado na forma de três artigos, cada um descrevendo uma etapa no desenvolvimento do problema da inversão, seguindo uma sequência de complexidade crescente no problema direto. O primeiro artigo neste trabalho é intitulado “Inversão de dados do CSEM marinho 1D de meio estratificado anisotrópico com o regularizador Variação Total”. Este descreve o passo inicial no desenvolvimento do problema: a inversão de dados do CSEM marinho de modelos estratificados 1D com anisotropia na condutividade das camadas. Este problema se presta bem para este desenvolvimento, porque tem solução computacional muito mais rápida do que o 2D, e nele já estão presentes as características principais dos dados do método CSEM marinho, como a largura muito grande da faixa de amplitudes medidas em um levantamento, e a baixa resolução, inerente às baixas frequências empregadas. A anisotropia acrescenta uma dificuldade a mais no problema, por aumentar o nível de ambiguidade nos dados e demandar ainda mais informação do que no caso puramente isotrópico. Os resultados mostram que a aplicação dos vínculos de igualdade do método VT permite a melhor identificação de uma camada alvo resistiva do que a simples aplicação dos vínculos tradicionais de suavidade. Até onde podemos aferir, esta solução se mostra superior a qualquer outra já publicada para este problema. Além de ter sido muito importante para o desenvolvimento de códigos em paralelo. O segundo artigo apresentado aqui, “Inversão de dados Magnetotelúricos com o regularizador Variação Total e o uso da matriz de sensibilidade aproximada”, trata da inversão de dados do método Magnetotelúrico em ambientes 2D. Este problema demanda um esforço computacional muito maior do que o primeiro. Nele, estudamos a aplicação do método dos estados adjuntos para gerar uma boa aproximação para as derivadas necessárias para a construção da matriz de sensibilidade usada na inversão. A construção da matriz de sensibilidade é a etapa que demanda mais tempo no processo de inversão, e o uso do método de estados adjuntos foi capaz de reduzir muito este tempo, gerando derivadas com um bom nível de aproximação. Esta etapa da pesquisa foi fundamental pelo problema direto ser matematicamente e computacionalmente muito mais simples do que o do CSEM marinho 2D. Novamente em comparação com a aplicação do regularizador de suavidade global, o regularizador de Variação Total permitiu, neste problema, uma melhor delimitação das bordas de heterogeneidades bidimensionais. A terceira parte deste trabalho, apresentada no artigo “Inversão de dados do CSEM marinho 2.5D com o regularizador Variação Total e o uso da matriz de sensibilidade aproximada”, apresenta a apliação do método de Variação Total ao problema da inversão de dados CSEM marinho 2.5D. Usamos o método dos estados adjuntos para gerar uma boa aproximação para as derivadas necessárias para a construção da matriz de sensibilidade usada na inversão, acelerando assim o processo de inversão. Para deixar o processo de inversão ainda mais rápido, lançamos mão da programação em paralelo com o uso de topologia. A comparação entre a aplicação do regularizador de suavidade global, e o regularizador de Variação Total permitiu, assim como nos casos anteriores, uma melhor delimitação das bordas de heterogeneidades bidimensionais.