Processamento interpretativo de dados magnetométricos e inversão de dados gravimétricos aplicados à prospecção de hidrocarbonetos

Apresentamos dois métodos de interpretação de dados de campos potenciais, aplicados à prospecção de hidrocarbonetos. O primeiro emprega dados aeromagnéticos para estimar o limite, no plano horizontal, entre a crosta continental e a crosta oceânica. Este método baseia-se na existência de feições geológicas magnéticas exclusivas da crosta continental, de modo que as estimativas das extremidades destas feições são usadas como estimativas dos limites da crosta continental. Para tanto, o sinal da anomalia aeromagnética na região da plataforma, do talude e da elevação continental é amplificado através do operador de continuação analítica para baixo usando duas implementações: o princípio da camada equivalente e a condição de fronteira de Dirichlet. A maior carga computacional no cálculo do campo continuado para baixo reside na resolução de um sistema de equações lineares de grande porte. Este esforço computacional é minimizado através do processamento por janelas e do emprego do método do gradiente conjugado na resolução do sistema de equações. Como a operação de continuação para baixo é instável, estabilizamos a solução através do funcional estabilizador de primeira ordem de Tikhonov. Testes em dados aeromagnéticos sintéticos contaminados com ruído pseudo-aleatório Gaussiano mostraram a eficiência de ambas as implementações para realçar os finais das feições magnéticas exclusivas da crosta continental, permitindo o delineamento do limite desta com a crosta oceânica. Aplicamos a metodologia em suas duas implementações a dados aeromagnéticos reais de duas regiões da costa brasileira: Foz do Amazonas e Bacia do Jequitinhonha. O segundo método delineia, simultaneamente, a topografia do embasamento de uma bacia sedimentar e a geometria de estruturas salinas contidas no pacote sedimentar. Os modelos interpretativos consistem de um conjunto de prismas bidimensionais verticais justapostos, para o pacote sedimentar e de prismas bidimensionais com seções verticais poligonais para as estruturas salinas. Estabilizamos a solução, incorporando características geométricas do relevo do embasamento e das estruturas salinas compatíveis com o ambiente geológico através dos estabilizadores da suavidade global, suavidade ponderada e da concentração de massa ao longo de direções preferenciais, além de vínculos de desigualdade nos parâmetros. Aplicamos o método a dados gravimétricos sintéticos produzidos por fontes 2D simulando bacias sedimentares intracratônicas e marginais apresentando densidade do pacote sedimentar variando com a profundidade segundo uma lei hiperbólica e abrigando domos e almofadas salinas. Os resultados mostraram que o método apresenta potencial para delinear, simultaneamente, as geometrias tanto de almofadas e domos salinos, como de relevos descontínuos do embasamento. Aplicamos o método, também, a dados reais ao longo de dois perfis gravimétricos sobre as Bacias de Campos e do Jequitinhonha e obtivemos interpretações compatíveis com a geologia da área.

Modelo fractal para resistividade complexa de rochas: interpretação petrofísica e aplicação à exploração geoelétrica

Rochas contendo metálicos disseminados ou partículas de argila em ambiente natural onde soluções eletrolíticas normalmente preenchem os poros das rochas, exibem um tipo de polarização em baixas freqüências conhecido como polarização induzida. Nesta tese foi desenvolvido um novo modelo para descrever o fenômeno de polarização das rochas, não apenas em baixas freqüências, mas compreendendo todo o espectro eletromagnético, possível de utilização na prospecção geoelétrica. Este novo modelo engloba a maioria dos modelos utilizados até o momento como casos especiais, além de superar as limitações dos mesmos. Seu circuito analógico inclui uma impedância não linear do tipo r (iwtf)^-n que simula o efeito das superfícies rugosas das interfaces entre os grãos bloqueadores (partículas metálicas e/ou de argilas) e o eletrólito. A impedância de Warburg generalizada está em série com a resistência dos grãos bloqueadores da passagem de corrente e em paralelo com a impedância da dupla camada associada a essas interfaces. Esta combinação está em série com a resistência do eletrólito nas passagens dos poros bloqueados. Os canais não bloqueados são representados por uma resistência que corresponde à resistividade normal CC da rocha. A combinação desta resistência com a capacitância "global" da rocha é finalmente conectada em paralelo ao resto do circuito mencionado acima. Os parâmetros deste modelo incluem a resistividade CC (p₀), a cargueabilidade (m), três tempos de relaxação (t, Tf and T2), um fator de resistividade de grãos (δ_r), e o expoente de freqüência (η). O tempo de relaxação fractal (Tf), e o expoente de frequencia (η) estão relacionados à geometria fractal das interfaces rugosas entre os minerais condutivos (grãos metálicos e/ou partículas de argila bloqueando os canais dos poros) e o eletrólito. O tempo de relaxação (T) é um resultado da relaxação em baixa freqüência das duplas camadas elétricas formadas nas interfaces eletrólito-cristais, enquanto (T₀) é o tempo de relaxação macroscópico da amostra como um todo. O fator de resistividade dos grãos (δ_r) relaciona a resistividade dos grãos condutivos com o valor de resistividade CC da rocha. A resistividade CC da rocha (p₀), e δ_r estão relacionados à porosidade, à condutividade do eletrólito e às relações mineralógicas entre a matriz e os grãos condutivos. O modelo foi testado sobre um intervalo largo de freqüências contra dados experimentais de amplitude e fase da resistividade bem como para dados de constante dielétrica complexa. Os dados utilizados neste trabalho foram obtidos a partir da digitalização de dados experimentais publicados, obtidos por diversos autores e englobando amostras de rochas sedimentares, ígneas e metam6rficas. É mostrado neste trabalho que os parâmetros deste modelo permitem identificar diferenças texturais e mineralógicas nas rochas. Bote modelo foi introduzido, primeiramente, como propriedade intrínseca de um semiespaço homogêneo sendo demonstrado, neste trabalho, que a resposta observada em superfície reflete as propriedades intrínsecas do meio polarizável, sendo o acoplamento eletromagnético desprezível em freqüências menores que 104 Hz. Em seguida, o meio polarizável foi embebido em um pacote de N camadas sendo demonstrado que os parâmetros fractais do meio polarizável podem ser obtidos do levantamento em superfície para diferentes espessuras dessa camada. Isto justifica a utilização pura e simples de modelos de polarização desenvolvidos para amostras em laboratório para ajustar dados de campo, o que vem sendo feito sem uma justificativa bem fundamentada. Estes resultados demonstram a importância para a prospecção geolétrica do modelo proposto nesta tese.

Detecção de fraturas no poço usando ferramenta de indução coplanar: modelamento analógico

A aplicabilidade de uma ferramenta de indução coplanar de duas bobinas na detecção de zonas fraturadas é estudado através de modelamento reduzido em laboratório. O fator de escala escolhido foi 10 e as condutividades de campo foram representadas por soluções de Cloreto de Amônia. O modelo geológico é composto por uma rocha homogênea resistiva com porções fraturadas e não-fraturada. As resistividades de campo são 4 Ω-m e 200 Ω-m para a zona fraturada e zona não-fraturada, respectivamente. As zonas fraturadas estudadas eram verticais, a várias distâncias da parede do poço, e cruzando a perfuração com mergulhos de 0°, 45° e 60°. A maioria dos perfis foram corridos em uma freqüência de laboratório de 100 kHz. Os resultados apresentados sob a forma de perfis mostraram que a ferramenta coplanar fornece informações sobre o mergulho, espessura e distância da parede do poço, dependendo da posição relativa das bobinas e da zona fraturada. Como a ferramenta coplanar não fornece informação de toda a formação ao redor da mesma, foi sugerido a construção de um sistema de duplo acoplamento semi-nulo (DASN) e alguns perfis são apresentados. Além do maior recobrimento, essa ferramenta fornece uma "técnica de olhada rápida" (quick-look technique) para estimar o mergulho da zona fraturada qualitativamente.

Inversão de dados de sísmica de refração profunda a partir da curva tempo-distância

O trabalho em pauta tem como objetivo o modelamento da crosta, através da inversão de dados de refração sísmica profunda, segundo camadas planas horizontais lateralmente homogêneas, sobre um semi-espaço. O modelo direto é dado pela expressão analítica da curva tempo-distância como uma função que depende da distância fonte-estação e do vetor de parâmetros velocidades e espessuras de cada camada, calculado segundo as trajetórias do raio sísmico, regidas pela Lei de Snell. O cálculo dos tempos de chegada por este procedimento, exige a utilização de um modelo cujas velocidades sejam crescentes com a profundidade, de modo que a ocorrência das camadas de baixa velocidade (CBV) é contornada pela reparametrização do modelo, levando-se em conta o fato de que o topo da CBV funciona apenas como um refletor do raio sísmico, e não como refrator. A metodologia de inversão utilizada tem em vista não só a determinação das soluções possíveis, mas também a realização de uma análise sobre as causas responsáveis pela ambiguidade do problema. A região de pesquisa das prováveis soluções é vinculada segundo limites superiores e inferiores para cada parâmetro procurado, e pelo estabelecimento de limites superiores para os valores de distâncias críticas, calculadas a partir do vetor de parâmetros. O processo de inversão é feito utilizando-se uma técnica de otimização do ajuste de curvas através da busca direta no espaço dos parâmetros, denominado COMPLEX. Esta técnica apresenta a vantagem de poder ser utilizada com qualquer função objeto, e ser bastante prática na obtenção de múltiplas soluções do problema. Devido a curva tempo-distância corresponder ao caso de uma multi-função, o algoritmo foi adaptado de modo a minimizar simultaneamente várias funções objetos, com vínculos nos parâmetros. A inversão é feita de modo a se obter um conjunto de soluções representativas do universo existente. Por sua vez, a análise da ambiguidade é realizada pela análise fatorial modo-Q, através da qual é possível se caracterizar as propriedades comuns existentes no elenco das soluções analisadas. Os testes com dados sintéticos e reais foram feitos tendo como aproximação inicial ao processo de inversão, os valores de velocidades e espessuras calculados diretamente da interpretação visual do sismograma. Para a realização dos primeiros, utilizou-se sismogramas calculados pelo método da refletividade, segundo diferentes modelos. Por sua vez, os testes com dados reais foram realizados utilizando-se dados extraídos de um dos sismogramas coletados pelo projeto Lithospheric Seismic Profile in Britain (LISPB), na região norte da Grã-Bretanha. Em todos os testes foi verificado que a geometria do modelo possui um maior peso na ambiguidade do problema, enquanto os parâmetros físicos apresentam apenas suaves variações, no conjunto das soluções obtidas.

Inversão em geofísica de poço: um estudo sobre ambiguidade

A ambiguidade na inversão de dados de geofísica de poço é estudada através da análise fatorial Q-modal. Este método é baseado na análise de um número finito de soluções aceitáveis, que são ordenadas, no espaço de soluções, segundo a direção de maior ambiguidade. A análise da variação dos parâmetros ao longo dessas soluções ordenadas permite caracterizar aqueles que são mais influentes na ambiguidade. Como a análise Q-modal é baseada na determinação de uma região de ambiguidade, obtida de modo empírico a partir de um número finito de soluções aceitáveis, é possível analisar a ambiguidade devida não só a erros nas observações, como também a pequenos erros no modelo interpretativo. Além disso, a análise pode ser aplicada mesmo quando os modelos interpretativos ou a relação entre os parâmetros não são lineares. A análise fatorial é feita utilizando-se dados sintéticos, e então comparada com a análise por decomposição em valores singulares, mostrando-se mais eficaz, uma vez que requer premissas menos restritivas, permitindo, desse modo, caracterizar a ambiguidade de modo mais realístico. A partir da determinação dos parâmetros com maior influência na ambiguidade do modelo é possível reparametrizá-lo, agrupando-os em um único parâmetro, redefinindo assim o modelo interpretativo. Apesar desta reparametrização incorrer na perda de resolução dos parâmetros agrupados, o novo modelo tem sua ambiguidade bastante reduzida.

Métodos geofísicos aplicados a prospecção de água subterrânea no município de Ponta de Pedras - Ilha de Marajó - Pará

Este trabalho foi realizado numa área de 25 km² em Ponta de Pedras, Ilha de Marajó. Tem por objetivo o dimensionamento geométrico e seleção de zonas favoráveis à água subterrânea. Foram utilizados dois métodos geofísicos: sísmica de refração e eletroresistividade. O método sísmico foi empregado em caráter experimental, visando verificar a viabilidade de sua aplicação na área. Os resultados da interpretação confirmaram a suspeita inicial, da existência de horizontes geológicos pouco profundos, não possíveis de serem distinguidos pelo método sísmico. Com o método de eletroresistividade foram feitas sondagens verticais e perfis horizontais. Os resultados da interpretação indicaram áreas promissoras na parte central e numa faixa ao norte da área prospectada. A espessura do aquífero raso varia entre 0 e 50 m.

Modelagem de efeitos plasmônicos em um metamaterial nanoestruturado para aplicação em sensoriamento óptico

Metamateriais são meios artificiais estruturados para a obtenção de propriedades físicas não prontamente encontradas na natureza. Estes meios, em geral, são organizados periodicamente com uma célula unitária que possui dimensões menores que o comprimento de onda da radiação eletromagnética incidente. Esta característica permite que estes meios possam ser tratados como materiais homogêneos e, por conseguinte, descritos em termos de parâmetros constitutivos efetivos. Este trabalho realiza a modelagem analítica e numérica de um metamaterial nanoestruturado com inclusões metálicas. É desenvolvido um formalismo fundamentado no modelo de Maxwell-Garnett para descrever o índice de refração efetivo do metamaterial levando em conta a forma geométrica das nanopartículas e a natureza dos modos locais excitados. É evidenciado que estes modos correspodem a excitações plasmônicas superfíciais. Por meio da caracterização do metamaterial no domínio óptico, é investigado o potencial de tais meios para o desenvolvimento de sensores de índice de refração mais compactos e com valores de sensibilidade superiores aos encontrados em sistemas baseados unicamente em fibras ópticas.

Modelamento da permissividade dielétrica de rochas saturadas de óleo e água e suas aplicações em perfilagem de poço

A ferramenta de propagação eletromagnética (EPT) fornece o tempo de propagação (T_pl) e a atenuação (A) de uma onda eletromagnética que se propaga num meio com perdas. Estas respostas da EPT são funções da permissividade dielétrica do meio. Existem vários modelos e fórmulas de misturas sobre a permissividade dielétrica de rochas reservatório que podem ser utilizados na interpretação da ferramenta de alta frequência. No entanto, as fórmulas de mistura não consideram a distribuição e a geometria do espaço poroso, e estes parâmetros são essenciais para que sejam obtidas respostas dielétricas mais próximas de uma rocha real. Foi selecionado um modelo baseado nos parâmetros descritos acima e este foi aplicado à dados dielétricos disponíveis na literatura. Foi obtida uma boa concordância entre as curvas teóricas e os dados experimentais, comprovando assim que a distribuição e a geometria dos poros têm que ser levadas em conta no desenvolvimento de um modelo realista. Foram conseguidas também funções de distribuição de razão de aspecto de poros, através das quais geramos várias curvas relacionando as respostas da EPT com diversas saturações de óleo/gás. Estas curvas foram aplicadas na análise de perfis. Como o modelo selecionado ajusta-se bem aos dados dielétricos disponíveis na literatura, torna-se atraente aplicá-lo à dados experimentais obtidos em rochas de campos brasileiros produtores de hidrocarbonetos para interpretação da EPT corrida em poços destes campos petrolíferos.

Pesagem em movimento e caracterização do tráfego ferroviário com uso da técnica B-WIM

Nesta dissertação é apresentado o desenvolvimento de algoritmos para aplicação do método Bridge-Weigh In Motion (B-WIM) para a pesagem em movimento de trens e para a caracterização do tráfego ferroviário, permitindo-se obter informações sobre a velocidade de passagem dos trens, número e espaçamento entre eixos. Os sistemas B-WIM a partir de uma simples instrumentação permitem determinar as cargas por eixo de veículos em movimento, eliminando o efeito dinâmico. Foram implementados os algoritmos para a determinação dos valores referentes a geometria do trem e das cargas, que foi validado a partir de um exemplo teórico, onde se simulou a passagem de um trem de características conhecidas sobre a ponte e as cargas por eixos foram determinadas com 100% de exatidão. Além disso, foi feito um exemplo numérico em elementos finitos, de um viaduto em concreto armado para aplicação do método, onde foi feita a determinação das cargas por eixo para diferentes velocidades de passagem do trem. A fim de reduzir o tempo de processamento nas análises do exemplo numérico, foi desenvolvido um algoritmo para a geração de cargas nodais no modelo numérico que reduziram o tempo de processamento em até 96% quando comparado com a análise de múltiplos passos (“Multi-Step”), que simula automaticamente a passagem do trem sobre a estrutura. Finalmente, o método foi testado em um caso real a partir de monitorações realizadas em um viaduto de concreto armado da Estrada de Ferro Carajás. Apesar de não ter sido possível a determinação das cargas por eixo da locomotiva, foi possível medir precisamente o peso bruto total da locomotiva quando se utilizou o modelo constitutivo de Collins & Mitchell (1991) para o concreto.

Deformação das rochas siliciclásticas paleoproterozoicas do Grupo Araí como exemplo das reativações de falhas do embasamento, Serra do Tepequém, Roraima, norte do Brasil

As rochas siliciclásticas da Serra do Tepequém são correlacionadas aos Grupos Araí e Suapi, pertencentes ao Supergrupo Roraima de idade paleoproterozoica. A análise estrutural destas rochas revelou que o acamamento exibe arranjos com mergulhos preferencialmente para SE e NW, individualizados em domínios limitados por zonas de falhas oblíquas sinistrais com rejeitos normais e inversos, com direção NE-SW. Essa estruturação é formada por dobras forçadas quilométricas do tipo kink bands e chevrons. O novo arcabouço geométrico observado é característico de um ambiente de deformação de nível crustal raso a médio. A história tectônica é controlada por reativações dos planos de fraqueza preexistentes nas tramas antigas, dúcteis, do embasamento. Esse modelo diverge dos regionais prévios para a região, os quais consideram as dobras existentes como produtos de ambiente dúctil sob tectônica colisional. Os resultados evidenciam a importância da presença de estruturas antigas do embasamento, relacionado ao Escudo das Guianas, como controladoras da geometria das rochas da Serra do Tepequém, em ambiente rúptil.

Sedimentação na bacia do Pantanal Mato-Grossense, Centro-Oeste do Brasil

The Pantanal is a Quaternary sedimentary basin located at the left margin of the Upper Paraguay River, west-central Brazil. Basin infilling was mainly by siliciclastic sediments and the stratigraphic succession exhibits an overall finingupward pattern. The depositional system tract is composed by a large meandering fluvial plain and several marginal alluvial fans, being the Taquari megafan the most striking feature. The present landscape is a complex tropical wetland, with geomorphic features derived from the present conditions and other inherited from successive Pleistocene and Holocene climates. During the Pleistocene, the sedimentary environment was dominated by braided alluvial fans, the original geometry of which is preserved as relict forms, permitting remarkable patterns of distributary paleochannels to be easily recognized in satellite images. Eolian processes were active in some abandoned lobes, contemporaneously with sedimentation in active fan lobes. Closed ponds bordered by lunette sand dunes, originally salt pans produced by eolian deflation, are relict eolian landforms in the Pantanal landscape. Eolian processes were probably more effective at the glacial maximum. Landscape has been changing in the Pantanal area since the end of the Pleistocene in adaptation to a more humid and warmer environment prevailing during Holocene. Initiation of the modern wetland has occurred during the Pleistocene / Holocene transition, with the change to a more humid climate and the individualization of lacustrine systems. The modern Pantanal wetland is a vast expanse of poorly drained lowlands that experiences annual flooding from summer to fall months. Although climatic fluctuations have occurred during all the Holocene, the alluvial fans have remained active depositional systems and lobes were formed by progradation and abandonment. Abandoned lobes were subjected... (Complete abstract click electronic address below)

Síntese, caracterização e estudo da atividade antitumoral de complexos de paládio(II) com ligantes sulfurados e trifenilfosfina

Pós-graduação em Química - IQ

Análise de robustez de estações da rede GNSS-SP

Currently, the need of reliable coordinates has been one of the main objectives of the scientific and practice community. Thus, the robustness analysis of a geodetic network, aims, at analyzing if the network is robust or not, based on the maximum undetectable errors. The network will be robust if the influence of these errors is small, otherwise it is weak, or not robust. This analysis is performed with the merger of two techniques, one which deals with the statistical analysis of reliability and the other one with the geometric strength analysis. The reliability analysis will provide the maximum error that cannot be detected by tests, after the adjustment. After finding these errors, the geometric strength analysis will determine the potential strain that the network will have, based on these errors. It is emphasized that the robustness analysis doesn't depend of the datum, reflecting only the geometry of the network and the accuracy of the observations (VANÌCEK et al., 2001). Therefore, this work aims at contributing to the scientific research on geodetic networks, checking the same, based on their geometry and observations.

Influência de parâmetros físicos e do processamento digital utilizado nos sinais oriundos do ensaio charpy instrumentado

Pós-graduação em Engenharia Civil - FEIS

Unsupervised manifold learning using Reciprocal kNN Graphs in image re-ranking and rank aggregation tasks

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)