Produção Científica



Dissertação de Mestrado
16/03/2022

Investigação de parâmetros para análise de velocidade de migração (MVA) usando dupla migração de integral de caminho (DPIM)
A migração de integral de caminho é um método para criar uma imagem migrada sem conhecimento prévio do modelo de velocidades por meio da soma ponderada de imagens migradas com um conjunto representativo de modelos de velocidade. Este conceito pode ser expandido para se extrair um modelo de velocidades, uma técnica denominada Análise de Velocidade de Migração via Dupla Migração de Integral de Caminho (Migration Velocity Analysis by Double Path Integral Migration – MVA via DPIM). Neste método, uma segunda imagem é criada, utilizando uma função peso ponderada pelo valor da velocidade empregada para a migração individual. Ao dividirmos esta segunda imagem pela primeira, extraimos o modelo de velocidades. Discutimos aqui diversos aspectos práticos da implementação de MVA by DPIM para migração em tempo, incluindo a parametrização da função peso, a estabilidade da divisão de imagens, a escolha da medida de coerência e a seleção do conjunto de velocidades relevantes. Através de testes realizados em uma ampla gama de modelos de velocidades, obtemos uma implementação robusta do método

Dissertação de Mestrado
16/03/2022

Redes Adversárias Generativas para Remoção de Ruídos em Dados Sísmicos
Dados sísmicos fornecem informações das características do ambiente geológico, que,após serem processados, são utilizados na interpretação para identificar e localizar as características do subsolo. O nível de complexidade desta tarefa é diretamente influenciado pela presença de ruídos, artefatos e/ou imprecisões em dados sísmicos. Logo, a remoção de ruídos e o aumento da qualidade destes dados são etapas essenciais no processamento sísmico. Para atenuar estes problemas, propomos abordagens baseadas em redes advers-rias generativas (generative adversarial networks, GANs) para melhorar os dados sísmicos através da remoção de ruídos e aumento de altas frequências nos dados. Nesse contexto,mapeamos o problema de remoção de ruído como um problema de tradução de imagem para imagem. Ou seja, utilizamos as GANs para aprender a mapear dados sísmicos com ruídos em dados sem ruídos, mantendo a integridade dos dados. Entretanto, o treinamento de redes neurais profundas (como as GANs) tem um problema inerente de escassez de dados rotulados. Para resolver esse problema, desenvolvemos um gerador de dados sísmicos sintéticos rotulados. Com esses dados, criamos nossos modelos baseados nas arquiteturas da Pix2Pix, Pix2PixHD e SPADE e utilizamos estes modelos para realizar a predição de dados sísmicos sintéticos e dados sísmicos reais. Nossos resultados, tanto em dados sísmicos sintéticos quanto em reais, sugerem que as abordagens propostas são promissoras, apresentando melhorias significativas na relação sinal-ruído e no ganho de altas frequências. A acurácia do nosso melhor modelo atingiu aproximadamente 87%, sendo este resultado considerado satisfatório para a utilização em cenários reais, de acordo com os especialistas. Com isso, é possível utilizar as GANs como ferramentas de apoio para os especialistas de sísmica, de modo a auxiliar as análises sísmicas e tomadas de decisões de forma mais rápida e assertiva

Dissertação de Mestrado
16/03/2022

Continuação da velocidade em família de ponto imagem comum em profundidade
Comparamos duas abordagens iterativas para análise de velocidade por migração em profundidade através da continuação da onda-imagem em famílias de ponto imagem comum (Common-Image Gathers - CIGs). Sendo elas, a Atualização Global, que visa atualizar todo o modelo de uma vez em cada iteração e o Layer Stripping, que atualiza o modelo uma camada por vez, de cima para baixo. Investigações anteriores mostraram que a atualização global das velocidades produz um modelo de velocidade razoável na primeira iteração, mas apresenta problemas de convergência nas iterações subsequentes. Por este motivo implementamos a abordagem Layer Stripping. Comparamos as duas abordagens iterativas numericamente, testando as implementações em dados sintéticos, a partir de um modelo bastante simples (mas representativo), que simula duas camadas sedimentares entre uma camada de água no topo e uma camada de sal no fundo. Nossos exemplos numéricos demonstram que a abordagem Layer Stripping melhora consideravelmente a convergência, tanto com, como sem o silenciamento dos afastamentos longos nos CIGs continuados. Sendo este um procedimento barato, que começa sem nenhum conhecimento a priori do meio, ele pode ser usado para construir modelos iniciais para técnicas de inversão mais sofisticadas.

Dissertação de Mestrado
16/03/2022

Inversão de Marchenko em uma dimensão.
A função de Green focada em um ponto arbitrário pode ser obtida por meio da solução da equação integral de Marchenko ou das equações acopladas de Marchenko, utlizando-se métodos iterativos. A equação integral de Marchenko é uma equação integral unidimensional que relaciona a amplitude de espalhamento e o coeficiente de reflexão que deu origem ao campo espalhado, sendo assim, um problema de espalhamento inverso. As equações acopladas de Marchenko são oriundas da relação entre o dado s??smico (função de Green) e uma solução fundamental de um problema de propagação (equação de Helmholtz homogênea com perturbação na velocidade). Este trabalho de dissertação visa reunir as deduções teóricas da equação integral de Marchenko e das equações acopladas de Marchenko, bem como apresentar o resultado da implementação dos dois esquemas iterativos e a comparação entre eles. Além disso, propõe uma aplicação para detecção de camadas utilizando as funções de Green unidirecionais redatumadas

Dissertação de Mestrado
16/03/2022

Uma comparação de métodos para transformação de seções sísmicas de afastamento comum em seções de afastamento nulo
Nesta dissertação fazemos uma comparação entre os primeiros métodos que surgiram com o objetivo de transformar seções sísmicas de afastamento comum em seções sísmicas de afastamento nulo (TZO). Os primeiro métodos que surgiram na literatura consideravam um meio isotrópico, velocidade e densidade constantes, e apenas tinham o objetivo de que a cinemática da transformação estivesse correta. Posteriormente com os conceitos de amplitude verdadeira, surgem métodos TZO de amplitude verdadeira. Com exceção do primeiro método apresentado (Yilmaz and Claerbout, 1980), os demais métodos fazem uma transformação cinemática perfeita. Para comparar a qualidade dos resultados obtidos por métodos de TZO de amplitude verdadeira, são feitos experimentos numéricos em três dos seis métodos apresentados. O primeiro é um método que leva em conta apenas a transformação cinemática, o segundo é um dos primeiros a trazer o conceito de amplitude verdadeira em TZO na literatura, mas considera que o refletor seja plano no cômputo do filtro de empilhamento do método, e o último método usa o mesmo conceito de amplitude verdadeira do segundo, entretanto considera o efeito da curvatura do refletor na amplitude, além de ser o único dos métodos apresentados que pode ser aplicado em um meio de velocidade não constante. Como era esperado, os experimentos numéricos mostram que cinematicamente os três métodos de TZO são idênticos, mas com relação à amplitude, os que utilizam do conceito de amplitude verdadeira são melhores. Curiosamente, mesmo que o segundo método não leve em conta a curvatura do refletor no efeito da amplitude, os resultados obtidos pelo mesmo e o terceiro método são praticamente iguais. Foi matematicamente provado que, de fato, o segundo e o terceiro métodos são assintoticamente equivalentes quando aplicados em um dado sísmico de conteúdo de alta frequência e em um meio de velocidade constante, inclusive, foi feito um experimento numérico mostrando que em situações que tanto o segundo quanto o terceiro método apresentam um maior erro de amplitude, o uso de um dado de mais alta frequência resulta em um erro menor ainda. Por último, foi feito um experimento numérico com o terceiro método em um modelo de velocidade com gradiente vertical constante. O experimento mostrou que a acurácia na amplitude atingida é tão boa quanto a obtida na aplicação do método em modelos de velocidade constante

Dissertação de Mestrado
16/03/2022

Redatumação utilizando retropropagação dos campos de onda
A redatumação consiste em mover fontes e receptores para um novo nível de profundidade. Esse objetivo pode ser alcançado por diferentes metodologias, como as que usam as técnicas interferométricas, ou por continuação dos campos de onda. Nesse trabalho, propomos usar a redatumação com base na retropropagação dos dados sísmicos. Para aplicar essa técnica, é necessário o conhecimento do modelo de velocidade da superfície de aquisição ao datum. Para validar a teoria, realizamos testes em dado sintéticos de diferentes complexidades e comparamos os resultados com testes usando a redatumação por correlação. Foi possível observar que a técninca de redatumação por retropropagação apresentou resultados bastante similar aos obtidos usando redatumação por correlação, tanto no posicionamento dos eventos, quanto na adição de artefatos ao resultado final. Também observamos que as técnicas apre- sentam a mesma sensibilidade quando se usa um modelo de velocidade inexato. Apesar das semelhanças, a técnica baseada na continuação dos campos de onda apresentou um tempo computacional muito menor, quando comparado a técnica interferométrica

Dissertação de Mestrado
16/03/2022

Caracterização de ruído sísmico ambiente em um campo de petróleo
Utilizamos um experimento de 5 horas de duração com 182 sensores verticais de 2 Hz implantados na superfície para caracterizar o ruído antes e durante um monitoramento de fraturamento hidráulico na Bacia Potiguar, Brasil. Observamos que o ruído sísmico é principalmente de induções eletromagnéticas e de vibração de máquinas perto da cabeça do poço, e de fontes localizadas a 2 km do centro da rede sísmica como as bombas de vareta de sucção, tubulações, estradas e instalações industriais. Também relatamos uma ressonância composta por ondas de corpo provenientes da área tratada, que só está presente quando a injeção ocorre. Interpretamos essa ressonância resultante de fraturas cheias de fluido na subsuperfície, mantendo ondas que reverberam nas fraturas. Diferentes estratégias foram empregadas para correlacionar e empilhar os dados para a Interferometria Sísmica de Ruído Ambiente: correlação cruzada normalizada geométrica (CCGN), correlação cruzada de fase (PCC), empilhamento linear e empilhamento ponderado de fase (PWS). PCC e PWS são baseados na coerência instantânea de fase dos sinais analíticos. Por causa da distribuição inadequada da fonte de ruído e da geometria da aquisição, artefatos surgem nos correlogramas. Nós propusemos uma metodologia simples para atenuar esses efeitos indesejados, que consiste em aplicar a correção estática (Linear Moveout - LMO), empilhar os dados no domínio do tiro e filtro f-k. As curvas de dispersão após esse processamento são aprimoradas e os resultados da análise de velocidade de fase são consistentes com os dados da literatura.

Dissertação de Mestrado
16/03/2022

Evidências geofísicas de domeamento na estruturação profunda da porção central da Faixa Seridó (NE-Brasil)
A Faixa Seridó (FS) é um cinturão de rochas metassedimentares localizado na Província Borborema (NE do Brasil). Ela foi deformada e metamorfizada no contexto da Orogênese Pan-Africana/Brasiliana. Neste evento tectônico, a crosta foi parcialmente fundida, particularmente na conexão da zona de cisalhamento Patos com a FS, alterando suas propriedades reológicas e, consequentemente, sua deformação. O domo anatético de Santa Luzia (~ 575 Ma) se destaca dentre as ocorrências de rochas associadas à fusão parcial. O mapa aeromagnético da FS mostra uma expressiva anomalia regional de longo comprimento de onda (> 25 km) em sua porção centro-sul. Esta anomalia apresenta forma sigmoidal e seu eixo central rotaciona de NNE-SSW para E-W próximo a zona de cisalhamento Patos. Anomalias de curto comprimento (< 25 km) e de grande amplitude (até 2.200 nT) se sobrepõem a anomalia regional. Os limites das anomalias magnéticas, frequentemente, coincidem com zonas de cisalhamento regionais, revelando que estas estruturas exercem importante controle tectono-estrutural sobre as fontes das anomalias. A comparação com dados gravimétricos evidencia que ambos os tipos de anomalias compartilham das mesmas fontes anômalas, sendo estas mais magnéticas e menos densas que a crosta circundante. Por outro lado, a correlação com a geologia de superfície mostra que as anomalias geofísicas muitas vezes ocorrem associadas com migmatitos e granitos ricos em magnetita. Nós realizamos duas abordagens de modelagem para os dados geofísicos: a primeira foi uma inversão automática 3D dos dados magnéticos e a segunda foi uma modelagem conjunta 2.5D, guiada pelo interprete, dos dados magnéticos e gravimétricos ao longo de três perfis. Como resultado da análise integrada dos dados geofísicos e geológicos e da modelagem, nós evidenciamos a existência de uma estrutura dômica regional na porção centro-sul da FS. Este domo gnáissico em grande escala é composto por quatro estruturas internas alongadas, geofisicamente modeladas como quatro corpos anômalos, além do material circundante. Em conjunto, eles formam uma estrutura crustal com eixo principal na direção NE-SW que rotaciona para E-W, seguindo a zona de cisalhamento Patos. O domo de Santa Luzia é apenas um domo local em pequena escala dessa estrutura. Os metassedimentos do Grupo Seridó ocorrem em sinformes em torno das estruturas que compõem o domo regional. Em particular, devido ao notável papel da zona de cisalhamento Patos na formação da grande estrutura crustal, é reforçada a concepção de que este lineamento foi um limite tectônico importante durante a Orogênese Pan-Africana/Brasiliana.

Dissertação de Mestrado
16/03/2022

Monitoramento de mudanças de velocidade no Arquipélago São Pedro São Paulo com interferometria de ruído sísmico
O Sistema Transformante São Paulo (STSP) está localizado na Dorsal Meso-Atlântica (MAR), Atlântico Equatorial e é formado por um complexo sistema de falhas, no qual suas interações provocam a ocorrência de eventos sísmicos com magnitudes pequenas a moderadas. A ausência de uma rede de estações sismográficas permanentes na região dificulta a caracterização da sismicidade neste local. Como resultado, os efeitos da sismicidade, de falhamento e de interação de fluidos na massa rochosa permanece pobremente amostrada. Nesta dissertação, nós quantificamos as perturbações de velocidade associadas à atividade sísmica no Arquipélago São Pedro São Paulo (ASPSP) usando uma estação sismográfica instalada na ilha Belmonte. Foram calculadas funções de autocorrelação (ACFs) de ruído a partir dos registros contínuos dessa estação ao longo do ano de 2012. As ACFs foram obtidas por duas abordagens diferentes: autocorrelação clássica normalizada geometricamente (ACGN) e autocorrelação por fase (PAC). Ambas as abordagens forneceram resultados similares e fomos capazes de estimar variações de velocidades devido a sismos com 3.0 ≤ ML ≤ 4.7. As mudanças no meio foram investigadas através da análise conjunta de curvas de decorrelação e de variações de velocidade, esta última sendo obtida pelo método de Moving Window Cross Spectral (MWCS). Nossas análises mostram que as variações de velocidade são possivelmente associadas ao acoplamento hidromecânico, no qual os aumentos de velocidade observados podem ser controlados por mudanças no esforço estático e interações rocha-fluido são responsáveis pelas quedas de velocidade observadas. Esses comportamentos são qualitativamente explicados pela natureza do complexo sistema de multi-falhas do STSP, onde interações rocha-fluido exercem um importante papel nas variações de velocidade observadas.

Dissertação de Mestrado
16/03/2022

Sismicidade do sistema transformante de São Paulo, no Atlântico Equatorial, de janeiro - julho de 2013
O Atlântico Equatorial é formado por diversas dorsais e falhas transformantes de escorregamento lento. Entre elas, o Sistema Transformante de São Paulo (STSP) é um complexo transformante, que se estende por 630 km, formado por múltiplas transformantes sendo quatro falhas e três segmentos intra-transformante. Na parte norte do STSP, a Zona da Serra do Atobá (ZSA) com 200 km de comprimento e 30 km de largura, é uma importante formação que atinge o nível do mar formando o Arquipélago de São Pedro e São Paulo (ASPSP). Foi realizado a localização hipocentral de 62 tremores do STSP. Eles ocorreram no ano de 2013, e foram registrados por um sismômetro instalado no ASPSP e três hidrofones implantados durante o cruzeiro COLMEIA. Usando os hipocentros, foi identificado uma zona sismogênica com profunda de transição frágil-dúctil à 18 km abaixo do oceano. Observamos que essa estrutura litosférica apresenta relação com a idade de deslocamento da transformante, na qual pode controlar as profundidades hipocentrais nas falhas transformantes do oceano. Além disso, os terremotos indicaram a existência de uma ampla extensão na profundidade da área de serpentinização, atingindo até 18 km abaixo da ZSA. Nós interpretamos isso como um efeito da percolação de águas em falhas e que possibilitam que elas atinjam até o manto abaixo do STSP, o que causa uma interação de rochas fluido-manto e o crescimento das falhas até o manto. Alguns hipocentros estavam localizados no segmento da Zona de Fratura Central (ZFC) e com profundidades que atingiam 8,8 km abaixo do fundo do mar. Nós interpretamos isso como a reativação de uma zona de fraqueza existente na ZFC, na qual ocorreu devido estresse induzido pela carga transpressiva da ZSA.

<<  <   1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36  37  38  39  40  41  42  43  44  45  46  47  48  49  50  51  52  53  54  55  56  57  58  59  60  61  62  63  64  65  66  67   >  >>