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.

Análise de contaminantes de fase líquida não aquosa (NAPLs) por aplicação do método eletromagnético indutivo (EM)

O descarte indevido de resíduos industriais no solo tem causado danos ao meio ambiente e à saúde humana. A identificação e o diagnóstico de áreas contaminadas é um procedimento necessário afim de que seja possível tomar medidas de contenção da fonte poluidora e posterior remediação. As técnicas geofísicas podem auxiliar na caracterização do meio geológico, definição das dimensões da pluma de contaminação, seleção de pontos para locação de poços de monitoramento e de bombeamento. Este trabalho apresenta os resultados da aplicação do método Eletromagnético Indutivo em uma área industrial contaminada por benzeno, tolueno, xileno (BTX) e dicloroetano. Os resultados indicam dois padrões de migração dos contaminantes em subsuperfície. Uma tendência de migração no plano horizontal até o nível de 15 m de profundidade e uma tendência de migração no plano vertical entre 15 e 60 m de profundidade. As tendências de migração divergentes estão condicionadas a densidade dos compostos que constituem a pluma de contaminação, ou seja, um intervalo LNAPL sobrenadante (BTX), sobreposto a um intervalo DNAPL (1,2 dicloroetano).

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

Aplicação do método eletromagnético indutivo (EM) no monitoramento de contaminantes em subsuperfície

O uso combinado de poços de amostragem e técnicas geofísicas constitui a maneira mais adequada para identificação e o monitoramento de áreas contaminadas. Este procedimento permite a locação de poços de amostragem e de monitoramento em pontos estratégicos, otimizando resultados e minimizando gastos. O método geofísico Eletromagnético Indutivo (EM) apresenta amplas possibilidades de aplicação em estudos ambientais devido à facilidade na aquisição de dados, versatilidade do equipamento em campo e a possibilidade de varredura de grandes áreas num curto espaço de tempo. Este trabalho realiza uma análise comparativa de dados em obtidos no ano de 1992 e em 2003. O local estudado é uma área industrial que apresenta contaminação do solo e água subterrânea por Benzeno, Tolueno, Xileno, 1,2 dicloroetano, Sódio e Cloreto, produto da infiltração de efluentes químicos diretamente do solo. Os resultados indicam uma acentuada atenuação da pluma contaminante, com provável redução no grau de contaminação. Os compostos de fase leve (LNAPLs) apresentaram maior tendência de migração horizontal, concomitante ao movimento da água subterrânea. Os compostos de fase densa (DNAPLs) apresentaram uma tendência mais acentuada de migração no sentido vertical, possivelmente devido à ausência de superfícies impermeáveis. Os compostos inorgânicos acompanharam o fluxo dos compostos de fase líquida não aquosa (NAPLs), o que por sua vez permitiu a caracterização da pluma de contaminação como um corpo condutivo em relação às áreas adjacentes.

Aplicação do método eletromagnético transiente em área contaminada sob atenuação natural

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

Emprego do método eletromagnético (em) na investigação e mapeamento de contaminações geradas por cemitérios: o caso de Monte Alto/SP

The use of geophysical methods in the analysis and mapping of likely contamination plumes has shown a huge value in the development of researches related to environmental issues. Among the main geophysical methods applied to this purpose, the geoelectrical methods stand out, being based on the use of electric and electromagnetic fields. The electromagnetic method (EM) uses physical environment responses created by electromagnetic field propagation, consisting of an alternated electric output and a magnetization force, used to obtain the electric conductivity (σ) of the subsurface materials. It has a broad application in the mineral exploration and prospection, and its use has been also efficient in environmental researches. By the electromagnetic method, this work examined and mapped possible contaminations created by cemetery leachate in the Saudade's Cemetery, located in the city of Monte Alto, Sao Paulo, where sedimentary rocks of the Bauru Group are exposed. Cemetery leachate is a putrefactive liquid, generated from the decomposition of organic matter, which contains substances that can contaminate soil and groundwater. Results achieved in this work show that the use of the electromagnetic method is an important indirect tool in the investigation of environmental contaminations in areas with highly contrasting physical properties

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.

Métodos espectronodais para cálculos de transporte de partículas neutras com fonte fixa na formulação de ordenadas discretas e multigrupo de energia

Um método espectronodal é desenvolvido para problemas de transporte de partículas neutras de fonte fixa, multigrupo de energia em geometria cartesiana na formulação de ordenadas discretas (SN). Para geometria unidimensional o método espectronodal multigrupo denomina-se método spectral Greens function (SGF) com o esquema de inversão nodal (NBI) que converge solução numérica para problemas SN multigrupo em geometria unidimensional, que são completamente livre de erros de truncamento espacial para ordem L de anisotropia de espalhamento desde que L < N. Para geometria X; Y o método espectronodal multigrupo baseia-se em integrações transversais das equações SN no interior dos nodos de discretização espacial, separadamente nas direções coordenadas x e y. Já que os termos de fuga transversal são aproximados por constantes, o método nodal resultante denomina-se SGF-constant nodal (SGF-CN), que é aplicado a problemas SN multigrupo de fonte fixa em geometria X; Y com espalhamento isotrópico. Resultados numéricos são apresentados para ilustrar a eficiência dos códigos SGF e SGF-CN e a precisão das soluções numéricas convergidas em cálculos de malha grossa.

Compressão de dados sísmicos com perda controlada

Neste trabalho apresentamos um novo método de compressão, com perda controlada de dados, que tem a vantagem de ter uma taxa significativa de compressão sem introduzir nenhuma perda superior a um parâmetro escolhido pelo usuário. Esta abordagem é uma abordagem mista, pois usa técnicas de compactação de dados tanto com perda quanto sem perda. Isto quer dizer que conseguimos um método que alia as vantagens da alta compressão, sem introduzir distorções indesejáveis nos dados. Mostramos como a massa de dados utilizada nos nossos estudos é obtida e a sua importância na prospecção de depósitos de hidrocarbonetos. É apresentado um levantamento bibliográfico com técnicas de compressão aplicadas a dados sísmicos tipicamente utilizadas em aplicações comerciais. Por fim, apresentamos os resultados da compressão utilizando o método em conjuntos de dados sísmicos reais. Para 1% de erro, os arquivos de dados sísmicos compactados passaram a ter algo próximo a 25% de seus tamanhos originais, o que representa um fator de compressão de aproximadamente 4

Insensibilização de frangos de corte em atmosfera controlada na promoção do bem estar animal: Diego Augusto Bitencourt. -

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

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.