Programa de Pós-Graduação em Geofísica - CPGF/IG

URI Permanente desta comunidadehttps://repositorio.ufpa.br/handle/2011/2355

O Programa de Pós-Graduação em Geofísica da UFPA (CPGF) do Instituto de Geociências (IG) da Universidade Federal do Pará (UFPA). Foi o segundo no Brasil a formar recursos humanos em Geofísica em nível de pós-graduação stricto sensu. Criado em 1972, funcionou até 1992 junto com os Cursos de Pós-Graduação em Geoquímica e Geologia.

Navegar

Agora exibindo 1 - 16 de 16

Acesso aberto (Open Access)
Análise do efeito da discretização do modelo de velocidades nas migrações Kirchhoff e Kirchhoff-Gaussian- Beam 2D pré-empilhamento em profundidade
(Universidade Federal do Pará, 2014-02-28) PAIXÃO, Marcelo Tavares; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
O Feixe Gaussiano (FG) é uma solução assintótica da equação da elastodinâmica na vizinhança paraxial de um raio central, a qual se aproxima melhor do campo de ondas do que a aproximação de ordem zero da Teoria do Raio. A regularidade do FG na descrição do campo de ondas, assim como a sua elevada precisão em algumas regiões singulares do meio de propagação, proporciona uma forte alternativa no imageamento sísmicos. Nesta dissertação, apresenta-se um novo procedimento de migração sísmica pré-empilhamento em profundidade com amplitudes verdadeiras, que combina a flexibilidade da migração tipo Kirchhoff e a robustez da migração baseada na utilização de Feixes Gaussianos para a representação do campo de ondas. O algoritmo de migração proposto é constituído por dois processos de empilhamento: o primeiro é o empilhamento de feixes (“beam stack”) aplicado a subconjuntos de dados sísmicos multiplicados por uma função peso definida de modo que o operador de empilhamento tenha a mesma forma da integral de superposição de Feixes Gaussianos; o segundo empilhamento corresponde à migração Kirchhoff tendo como entrada os dados resultantes do primeiro empilhamento. Pelo exposto justifica-se a denominação migração Kirchhoff-Gaussian-Beam (KGB).Afim de comparar os métodos Kirchhoff e KGB com respeito à sensibilidade em relação ao comprimento da discretização, aplicamos no conjunto de dados conhecido como Marmousi 2-D quatro grids de velocidade, ou seja, 60m, 80m 100m e 150m. Como resultado, temos que ambos os métodos apresentam uma imagem muito melhor para o menor intervalo de discretização da malha de velocidade. O espectro de amplitude das seções migradas nos fornece o conteúdo de frequência espacial das seções das imagens obtidas.
Acesso aberto (Open Access)
Estudo do método de inversão tomográfica 2D baseada nos atributos cinemáticos da onda NIP e construção de interface gráfica para o processo de picking
(Universidade Federal do Pará, 2007) OLIVEIRA, Iran Gadelha de; CALLAPINO, German Garabito; http://lattes.cnpq.br/6064981270181319
Neste trabalho, estuda-se um novo método de inversão tomográfica de reflexão para a determinação de um modelo isotrópico e suave de velocidade por meio da aplicação, em dados sintéticos e reais, do programa Niptomo que é uma implementação do método de inversão tomográfica dos atributos cinemáticos da onda hipotética do ponto de incidência normal (PIN). Os dados de entrada para a inversão tomográfica, isto é, o tempo de trânsito e os atributos da onda PIN (raio de curvatura da frente de onda emergente e ângulo de emergência), são retirados de uma série de pontos escolhidos na seção afastamento nulo (AN) simulada, obtida pelo método de empilhamento por superfícies de reflexão comum (SRC). Normalmente, a escolha destes pontos na seção AN é realizada utilizando-se programas de picking automático, que identificam eventos localmente coerentes na seção sísmica com base nos parâmetros fornecidos pelo usuário. O picking é um dos processos mais críticos dos métodos de inversão tomográfica, pois a inclusão de dados de eventos que não sejam de reflexões primárias podem ser incluídos neste processo, prejudicando assim o modelo de velocidades a ser obtido pela inversão tomográfica. Este trabalho tem por objetivo de construir um programa de picking interativo para fornecer ao usuário o controle da escolha dos pontos de reflexões sísmicas primárias, cujos dados serão utilizados na inversão tomográfica. Os processos de picking e inversão tomográfica são aplicados nos dados sintéticos Marmousi e nos dados da linha sísmica 50-RL-90 da Bacia do Tacutu. Os resultados obtidos mostraram que o picking interativo para a escolha de pontos sobre eventos de reflexões primárias favorece na obtenção de um modelo de velocidade mais preciso.
Acesso aberto (Open Access)
Inversão quadrática de dados de AVO usando algoritmo de Marquardt
(Universidade Federal do Pará, 2004-11-12) MENDES, Ana Cristina Marques; PROTÁZIO, João dos Santos; http://lattes.cnpq.br/4210442535067685
Neste trabalho são apresentadas expressões exatas e aproximações quadráticas das conhecidas equações de Knott-Zöeppritz, que calculam as amplitudes dos coeficientes Rpp e Rsp em termos dos contrastes médios relativos, bem como é procedida a inversão dos parâmetros de dados de AVO a partir destas aproximações quadráticas. Nesta inversão é utilizado o algoritmo de Levenberg-Marquardt, e são considerados apenas os eventos refletidos Rpp e convertidos Rsp, não associados. Nos estudos dos parâmetros físicos dos meios contrastes de impedância (δz), módulo de cisalhamento (δμ) e velocidade da onda p (δα), verificou-se quais desses parâmetros podem ser invertidos. Os resultados obtidos mostram que o contraste de impedância (δz) é muito bem resolvido estando ele relacionado com o contraste de velocidade da onda p (δα) ou com o contraste de cisalhamento (δμ), no caso de eventos refletidos considerando modelos de alto, moderado e baixo contrastes. Por outro lado ao se fixar o contraste de impedância e relacionar os outros dois parâmetros em consideração, os resultados mostram que esses são mal resolvidos, ou seja, a região de ambiguidade torna-se muito grande e os parâmetros tornam-se ambíguos e instáveis. No caso do evento convertido e na combinação do evento refletido com o convertido, para os modelos de baixo e moderado contrastes, (δz) é muito bem resolvido, caso que não acontece para modelo de alto contraste. Diante desses resultados verifica-se que no procedimento de inversão quadrática de dados de AVO, fixado (δμ), a recuperação dos dois parâmetros variados é muito boa, no caso do evento refletido, e razoavelmente boa no caso do evento convertido, por esse motivo optou-se pela fixação do módulo de cisalhamento.
Acesso aberto (Open Access)
Joining diffraction filter and residual diffraction moveout to constructo a velocity model in the depth and time domains: application to a Viking Graben data set
(Universidade Federal do Pará, 2014-08-14) COLLAZOS GONZALEZ, Jaime Andres; FIGUEIREDO, José Jadsom Sampaio de; http://lattes.cnpq.br/1610827269025210
Ondas sísmicas difratadas são geradas por descontinuidades na subsuperfície da Terra com o tamanho da ordem do comprimento de onda sísmico. Uma vez que o campo de onda incidente pode ser significativamente afetado por essas descontinuidades, muitas propriedades importantes destes eventos podem ser usadas para melhorar a prática de imageamento sísmico. Neste trabalho propomos uma abordagem prática para construir modelos de velocidade no domínio do tempo e profundidade usando difrações. Esta metodologia consiste na aplicação do filtro destrutor de onda plana (plane wave destruction - PWD) juntamente com método residual difraction moveout (RDM), de modo a construir modelos de velocidade nos domínios do tempo e da profundidade. Nosso método depende apenas de difrações (identificadas) filtradas a partir de eventos de reflexão e um modelo de velocidade inicial arbitrário de entrada. As imagens migradas no domínio pós-empilhado (nos domínios do tempo e da profundidade) são comparados com imagens migradas derivadas do processamento sísmico convencional. Nestes domínios, usamos a migracão Kirchhoff pós-empilhamento. Desconsiderando a necessidade de identificar e escolher os eventos de difração na migração pós-empilhamento no domínio da profundidade, o método apresenta um custo computacional muito baixo. Para alcançar um modelo de velocidade aceitável o tempo de processamento comparado ao método convencional foi menor. A viabilidade de nossa metodologia é testada num dado sísmico real do Viking Graben.
Acesso aberto (Open Access)
Medidas de coerência para análise de velocidade na migração em tempo
(Universidade Federal do Pará, 2011) MACIEL, Jonathas da Silva; COSTA, Jessé Carvalho; http://lattes.cnpq.br/7294174204296739
Métodos iterativos para a análise de velocidade com migração dependem das funções objetivo medirem a horizontalização dos eventos de reflexão em famílias de ponto imagem comum (CIG). A migração em tempo é um método simples de imageamento para avaliar essas funções objetivo. Usando migração em tempo, estudamos a influência da função objetivo nos resultados da análise de velocidade com migração. Propomos duas novas funções objetivo para a análise velocidade com migração: Semblance Diferencial Estendido e o produto do Semblance Clássico com o Semblance Diferencial Estendido. Experimentos numéricos, utilizando os dados Marmousoft, mostram a eficácia das novas funções objetivo em estimar modelos velocidade produzindo eventos horizontalmente alinhados nas famílias de ponto imagem comum.
Acesso aberto (Open Access)
Métodos de otimização global aplicados na busca dos parâmetros SRC
(Universidade Federal do Pará, 2004-06-21) CARMO, Klaus Mairan Laurido do; GARABITO CALLAPINO, German; http://lattes.cnpq.br/6064981270181319
O método de empilhamento por Superfície de Reflexão Comum (SRC) produz seções simuladas de afastamento nulo (AN) por meio do somatório de eventos sísmicos dos dados de cobertura múltipla contidos nas superfícies de empilhamento. Este método não depende do modelo de velocidade do meio, apenas requer o conhecimento a priori da velocidade próxima a superfície. A simulação de seções AN por este método de empilhamento utiliza uma aproximação hiperbólica de segunda ordem do tempo de trânsito de raios paraxiais para definir a superfície de empilhamento ou operador de empilhamento SRC. Para meios 2D este operador depende de três atributos cinemáticos de duas ondas hipotéticas (ondas PIN e N), observados no ponto de emergência do raio central com incidência normal, que são: o ângulo de emergência do raio central com fonte-receptor nulo (β0) , o raio de curvatura da onda ponto de incidência normal (RPIN) e o raio de curvatura da onda normal (RN). Portanto, o problema de otimização no método SRC consiste na determinação, a partir dos dados sísmicos, dos três parâmetros (β0, RPIN, RN) ótimos associados a cada ponto de amostragem da seção AN a ser simulada. A determinação simultânea destes parâmetros pode ser realizada por meio de processos de busca global (ou otimização global) multidimensional, utilizando como função objetivo algum critério de coerência. O problema de otimização no método SRC é muito importante para o bom desempenho no que diz respeito a qualidade dos resultados e principalmente ao custo computacional, comparado com os métodos tradicionalmente utilizados na indústria sísmica. Existem várias estratégias de busca para determinar estes parâmetros baseados em buscas sistemáticas e usando algoritmos de otimização, podendo estimar apenas um parâmetro de cada vez, ou dois ou os três parâmetros simultaneamente. Levando em conta a estratégia de busca por meio da aplicação de otimização global, estes três parâmetros podem ser estimados através de dois procedimentos: no primeiro caso os três parâmetros podem ser estimados simultaneamente e no segundo caso inicialmente podem ser determinados simultaneamente dois parâmetros (β0, RPIN) e posteriormente o terceiro parâmetro (RN) usando os valores dos dois parâmetros já conhecidos. Neste trabalho apresenta-se a aplicação e comparação de quatro algoritmos de otimização global para encontrar os parâmetros SRC ótimos, estes são: Simulated Annealing (SA), Very Fast Simulated Annealing (VFSA), Differential Evolution (DE) e Controlled Rando Search - 2 (CRS2). Como resultados importantes são apresentados a aplicação de cada método de otimização e a comparação entre os métodos quanto a eficácia, eficiência e confiabilidade para determinar os melhores parâmetros SRC. Posteriormente, aplicando as estratégias de busca global para a determinação destes parâmetros, por meio do método de otimização VFSA que teve o melhor desempenho foi realizado o empilhamento SRC a partir dos dados Marmousi, isto é, foi realizado um empilhamento SRC usando dois parâmetros (β0, RPIN) estimados por busca global e outro empilhamento SRC usando os três parâmetros (β0, RPIN, RN) também estimados por busca global.
Acesso aberto (Open Access)
Migração 3D no tempo usando a teoria dos raios paraxiais
(Universidade Federal do Pará, 1998-09) QUEIROZ, Norcirio Pantoja; SÖLLNER, Walter Franz
O trabalho aqui apresentado visa o estudo da migração 3D no tempo, em amplitudes verdadeiras, de dados da seção afastamento nulo ("zero-offset"), usando a aproximação de segunda ordem do tempo de trânsito de difração. O processo de migração no tempo usa a função difração determinada, corretamente, através de mesmos parâmetros que controlam a função tempo de trânsito de reflexão, além do auxílio dos raios paraxiais. Os parâmetros que controlam a função tempo de trânsito de reflexão são determinados através de pelo menos nove medidas de tempo de trânsito, obtidas a partir de dados pré-empilhados (tiro comum, ponto médio comum - CMP, etc.). Determinada a função difração, pode-se fazer a migração sem a necessidade de que se conheça o modelo de velocidades, pois a função tempo de trânsito de difração, como definido neste trabalho, não necessita do mesmo. O espalhamento geométrico é removido através do processo de migração e aplicação de uma função peso aos dados da seção afastamento nulo. Juntamente com a função peso, aplicou-se um filtro objetivando a recuperação da fase do sinal, ocorrida como consequência do processo de migração. Devido a limitação computacional, fez-se testes com dados sintéticos, em 2D, objetivando verificar a eficiência do programa. Os resultados obtidos foram satisfatórios, mostrando a eficiência e confiabilidade do processo.
Acesso aberto (Open Access)
Migração com amplitude verdadeira em meios bidimensionais (2-D) e introdução ao caso 2,5-D
(Universidade Federal do Pará, 1999) URBAN, Jaime Antonio; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
Nos últimos anos tem-se verificado através de várias publicações um interesse crescente em métodos de migração com amplitude verdadeira, com o objetivo de obter mais informações sobre as propriedades de refletividade da subsuperfície da terra. A maior parte desses trabalhos tem tratado deste tema baseando-se na aproximação de Born, como em Bleistein (1987) e Bleistein et al. (1987), ou na aproximação do campo de ondas pela teoria do raio como Hubral et al. (1991), Schleicher et al. (1993) e Martins et al. (1997). Considerando configurações arbitrárias de fontes e receptores, as reflexões primárias compressionais podem ser imageadas em reflexões migradas no domínio do tempo ou profundidade de tal modo que as amplitudes do campo de ondas migrado são uma medida do coeficiente de reflexão dependente do ângulo de incidência. Para realizar esta tarefa, vários algoritmos têm sido propostos nos últimos anos baseados nas aproximações de Kirchhoff e Born. Essas duas abordagens utilizam um operador integral de empilhamento de difrações ponderado que é aplicado aos dados da seção sísmica de entrada. Como resultado obtém-se uma seção migrada onde, em cada ponto refletor, tem-se o pulso da fonte com amplitude proporcional ao coeficiente de reflexão naquele ponto. Baseando-se na aproximação de Kirchhoff e na aproximação da teoria do raio do campo de ondas, neste trabalho é obtida a função peso para modelos bidimensionais (2-D) e dois e meio dimensionais (2,5-D) que é aplicada a dados sintéticos com e sem ruído. O resultado mostra a precisão e estabilidade do método de migração em 2-D e 2,5-D como uma ferramenta para a obtenção de informações importantes da subsuperfície da terra, que é de grande interesse para a análise da variação da amplitude com o afastamento (ângulo). Em suma, este trabalho apresenta expressões para as funções peso 2-D e 2,5-D em função de parâmetros ao longo de cada ramo do raio. São mostrados exemplos da aplicação do algoritmo de migração em profundidade a dados sintéticos 2-D e 2,5-D obtidos por modelamento sísmico através da teoria do raio usando o pacote Seis88 (Cervený e Psencík, 1988) e os resultados confirmaram a remoção do espalhamento geométrico dos dados migrados mesmo na presença de ruído. Testes adicionais foram realizados para a análise do efeito de alongamento do pulso na migração em profundidade (Tygel et al., 1994) e a aplicação do empilhamento múltiplo (Tygel et al., 1993) para a estimativa de atributos dos pontos de reflexão - no caso o ângulo de reflexão e a posição do receptor.
Acesso aberto (Open Access)
Migração FD 3D em profundidade usando aproximação de Padé complexa
(Universidade Federal do Pará, 2011-06-21) COSTA, Carlos Alexandre Nascimento da; COSTA, Jessé Carvalho; http://lattes.cnpq.br/7294174204296739
A implementação convencional do método de migração por diferenças finitas 3D, usa a técnica de splitting inline e crossline para melhorar a eficiência computacional deste algoritmo. Esta abordagem torna o algoritmo eficiente computacionalmente, porém cria anisotropia numérica. Esta anisotropia numérica por sua vez, pode levar a falsos posicionamentos de refletores inclinados, especialmente refletores com grandes ângulos de mergulho. Neste trabalho, como objetivo de evitar o surgimento da anisotropia numérica, implementamos o operador de extrapolação do campo de onda para baixo sem usar a técnica splitting inline e crossline no domínio frequência-espaço via método de diferenças finitas implícito, usando a aproximação de Padé complexa. Comparamos a performance do algoritmo iterativo Bi-gradiente conjugado estabilizado (Bi-CGSTAB) com o multifrontal massively parallel solver (MUMPS) para resolver o sistema linear oriundo do método de migração por diferenças finitas. Verifica-se que usando a expansão de Padé complexa ao invés da expansão de Padé real, o algoritmo iterativo Bi-CGSTAB fica mais eficientes computacionalmente, ou seja, a expansão de Padé complexa atua como um precondicionador para este algoritmo iterativo. Como consequência, o algoritmo iterativo Bi-CGSTAB é bem mais eficiente computacionalmente que o MUMPS para resolver o sistema linear quando usado apenas um termo da expansão de Padé complexa. Para aproximações de grandes ângulos, métodos diretos são necessários. Para validar e avaliar as propriedades desses algoritmos de migração, usamos o modelo de sal SEG/EAGE para calcular a sua resposta ao impulso.
Acesso aberto (Open Access)
Migração FFD 3D em profundidade usando aproximação de Padé complexa
(Universidade Federal do Pará, 2011) CAMPOS, Itamara do Socorro da Silveira; COSTA, Jessé Carvalho; http://lattes.cnpq.br/7294174204296739
Implementações dos métodos de migração diferença finita e Fourier (FFD) usam fatoração direcional para acelerar a performance e economizar custo computacional. Entretanto essa técnica introduz anisotropia numérica que podem erroneamente posicionar os refletores em mergulho ao longo das direções em que o não foi aplicado a fatoração no operador de migração. Implementamos a migração FFD 3D, sem usar a técnica do fatoração direcional, no domínio da frequência usando aproximação de Padé complexa. Essa aproximação elimina a anisotropia numérica ao preço de maior custo computacional buscando a solução do campo de onda para um sistema linear de banda larga. Experimentos numéricos, tanto no modelo homogêneo e heterogêneo, mostram que a técnica da fatoração direcional produz notáveis erros de posicionamento dos refletores em meios com forte variação lateral de velocidade. Comparamos a performance de resolução do algoritmo de FFD usando o método iterativo gradiente biconjugado estabilizado (BICGSTAB) e o multifrontal massively parallel direct solver (MUMPS). Mostrando que a aproximação de Padé complexa é um eficiente precondicionador para o BICGSTAB, reduzindo o número de iterações em relação a aproximação de Padé real. O método iterativo BICGSTAB é mais eficiente que o método direto MUMPS, quando usamos apenas um termo da expansão de Padé complexa. Para maior ângulo de abertura do operador, mais termos da série são requeridos no operador de migração, e neste caso, a performance do método direto é mais eficiente. A validação do algoritmo e as propriedades da evolução computacional foram avaliadas para a resposta ao impulso do modelo de sal SEG/EAGE.
Acesso aberto (Open Access)
Modelagem acústica por diferenças finitas e elementos finitos em 2-D e 2,5-D
(Universidade Federal do Pará, 2004-10-26) SILVA NETO, Francisco de Assis; COSTA, Jessé Carvalho; http://lattes.cnpq.br/7294174204296739
A modelagem acústica fornece dados úteis para avaliação de metodologias de processamento e imageamento sísmico, em modelos com estrutura geológica complexa. Esquemas de diferenças finitas (DF) e elementos finitos (EF) foram implementados e avaliados em modelos homogêneos e heterogêneos. O algoritmo de diferenças finitas foi estendido para o caso 2,5-D em modelos com densidade variável. Foi apresentada a modelagem de alvos geológicos de interesse exploratório existentes na Bacia Paleozóica do Solimões na Amazônia. Reflexões múltiplas de longo período produzidas entre a superfície livre e a discordância Cretáceo-Paleozóica, a baixa resolução da onda sísmica nas proximidades do reservatório e as fracas reflexões na interface entre as rochas reservatório e as rochas selantes são as principais características dos dados sintéticos obtidos, os quais representam um grande desafio ao imageamento sísmico.
Acesso aberto (Open Access)
Modelagem sísmica por diferenças finitas em meios bidimensionais com difratores
(Universidade Federal do Pará, 2000-09-14) FERNANDES, Lindemberg Lima; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
Este trabalho tem por objetivo a modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações. Como parte deste estudo foi feita a análise numérica visando as condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas. Para a validação do programa de diferenças finitas foi feita a comparação cinemática com a Teoria do Raio para um modelo simples. O motivo deste estudo é ter uma melhor compreensão e controle sobre os problemas de modelagem, visando contribuir para a solução de problemas na interpretação de dados sísmicos. Segundo vários autores na literatura geológica, Derby (1877), Evans (1906), Paiva (1929) e Moura (1938). A Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o recente, atingindo espessuras da ordem de 5 km. Os corpos de diabásio, presentes entre os sedimentos paleozóicos, estão dispostos na forma de soleiras, alcançando espessuras de centenas de metros, perfazendo um volume total de 90.000 quilômetros cúbicos. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica, impossibilitando uma melhor interpretação dos horizontes refletores que se encontram abaixo das soleiras. Para representar situações geológicas desse tipo foi usado um modelo acústico de velocidades. Para o cálculo dos sismogramas foi utilizado um programa de diferenças finitas com aproximação de quarta ordem da equação da onda acústica no espaço e no tempo. As aplicações dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas têm melhorado a compreensão sobre a propagação em meios onde existem heterogeneidades significativas, tendo como resultado boa resolução na interpretação dos eventos de reflexão sísmica em áreas de interesse. Como resultado dos experimentos numéricos realizados em meio de geologia complexa, foi observado a influência significativa das múltiplas devido a camada de alta velocidade, o que faz com que haja maior perda de energia dificultando a interpretação dos alvos. Por esta razão recomendo a integração de dados de superfície com os dados de poço, com o objetivo de se obter uma melhor imagem dos alvos abaixo das soleiras de diabásio.
Acesso aberto (Open Access)
Um novo algorítimo para implementação do método de empilhamento SRC 3-D e modelagem sísmica da Bacia do Amazonas
(Universidade Federal do Pará, 2006-05-23) LIMA, Charles Cristiano de Assis; CALLAPINO, German Garabito; http://lattes.cnpq.br/6064981270181319
O método de empilhamento Superfície de Reflexão Comum (SRC) é uma técnica de imageamento de reflexão sísmica que não depende do modelo de velocidades para simular uma seção ou um volume Afastamento-Nulo (AN) a partir dos dados de cobertura múltipla. O método SRC tem sido usado como uma alternativa atraente para o processamento de dados de reflexão e muitos trabalhos vem demostrando que este método simula seções ou volumes AN com melhor razão sinal/ruído e melhor resolução lateral que os obtidos pelos métodos convencionais, e.g. empilhamento NMO/DMO. O empilhamento SRC fornece importantes atributos ou parâmetros do campo de onda que podem ser usados em diversas aplicações como, por exemplo, na determinação do modelo de velocidades, análises AVO e A V A, determinação de fator de espalhamento geométrico, zonas de fresnel projetadas, entre outros. O método SRC utiliza como operador uma aproximação de segunda ordem dos tempos de trânsito hiperbólicos na vizinhança de um raio central. No caso 3-D, para um raio central normal ou AN, o operador SRC depende de oito parâmetros, que podem ser determinados através de análises de coerência a partir dos dados sísmicos. Neste trabalho é examinado o operador SRC-AN 3D para eventos de reflexão e difração com relação aos tempos de trânsitos obtidos por modelagem direta. Os resultados desta comparação demonstram que o operador SRC-AN 3D apresenta um bom ajuste com a superfície de tempos de trânsito obtida por modelagem direta. Com base nestes resultados numéricos e nas particularizações e simplificações da aproximação hiperbólica dos tempos de trânsito, é apresentado neste trabalho um novo algoritmo de empilhamento SRC-3D para simular um volume de dados com afastamento nulo. Nesta dissertação também é apresentada uma modelagem sísmica 3-D completa da Bacia do Amazonas com os pacotes FORMS, MODEL e CCRT, do consórcio Seismic Waves in Complex 3-D Structures - SW3D, de domínio público. O modelo foi construído com o pacote MODEL e foi inteiramente baseado em dados reais de sísmica e poço da bacia. O modelo é constituído por três soleiras de diabásio e um delgado reservatório, situação geológica típica da região. Com o pacote CRT foi realizado traçamento de raio 3-D e a geração de sismogramas sintéticos. Todos os modelos foram vizualizados com VRML e GOCAD. Pretende-se, em trabalhos futuros, utilzar este modelo para gerar dados sísmicos com cobertura múltipla e validar o método SRC-AN 3D, assim como outras técnicas de imageamento sísmico.
Acesso aberto (Open Access)
Processamento de dados sísmicos da bacia do Camamu
(Universidade Federal do Pará, 2009) MATOS, Ligia Naia de; GOMES, Ellen de Nazaré Souza; http://lattes.cnpq.br/1134403454849660
O imageamento sísmico é uma maneira de visualizar meios em superfície. Para se chegar a uma imagem sísmica confiável, em um longo caminho deve ser percorrido. A primeira etapa consiste na aquisição dos dados sísmicos. uma aquisição mal planejada ou com problemas na sua realização compromete todas as etapas de tratamento dos dados até o imageamento. Em seguida os dados sísmicos foram processados. O processamento sísmicos é a alteração dos dados sismicos para a supressão de ruídos, realçando o sinal e migrando os eventos sísmicos para as suas posições corretas. As etapas do processamento típicamente incluem a análise de velocidade e frequências, correção estática, deconvolução, correção de sobretempo, empilhamento e migração, esta pode ser realizada antes ou depois do empilhamento. Neste trabalho, nosso objetivo é a obtenção de uma seção migrada em tempo de dados 2D da bacia do Camamu. Estes dados obtidos na região de quebra do talude apresenta complexidades na sua interpretação devido a presença de eventos di múltiplas de superfície livre, múltiplas de curto período além de forte variação de velocidade. Para o imageamento, nós utilizamos a técnica do ponto médio comum e migração em tempo pós-empilhamento utilizando o software de processamento sísmico Promax (Landmark Inc). Para a atenuação das múltiplas utilizamos as técnicas de deconvolução preditiva, filtro radom e filtro f-k. Os resultados foram satisfatórios quando comparados com os resultados obtidos sem a aplicação dessas técnicas.
Acesso aberto (Open Access)
Processamento, inversão e imageamento de dados sísmicos marinhos
(Universidade Federal do Pará, 2012) SILVA, Douglas Augusto Barbosa da; LEITE, Lourenildo Williame Barbosa; http://lattes.cnpq.br/8588738536047617
Este trabalho constou do estudo e do processamento, inversão e imageamento em tempo e em profundidade de dados sísmicos Marmousi, e do real Jequitinhonha obtido na plataforma continental atlântica brasileira do estado da Bahia. Foram aplicados os métodos de empilhamento NMO e CRS, e a técnica de inversão niptomográfica dos atributos cinemáticos do campo de onda. Com o empilhamento NMO foi obtido um mapa de distribuição de velocidades através da análise de velocidade no mapa de coerência semblance e, posteriormente, as seções empilhada e migrada no domínio do tempo e da profundidade. O método de empilhamento da superfície de reflexão comum (CRS) foi aplicado com o programa crsstack-511 para obter as seções empilhada e migrada no domínio do tempo, e extrair os parâmetros do campo de onda, através de análise de coerência e da redundância contida nos dados sísmicos de multicobertura. O método de inversão niptomográfica de reflexão foi realizado por meio da aplicação do programa niptomo, que é uma implementação do método de inversão dos atributos cinemáticos da onda hipotética NIP, extraídos diretamente do empilhamento CRS, para se obter um modelo suavizado de velocidade e, posteriormente, uma seção migrada em profundidade. As seções migradas são do tipo Kirchhoff. As técnicas utilizadas seguem um fluxograma pré-estabelecido e realizado seguindo um arquivo Makefile, que funciona como um organizador de etapas. Estas etapas foram realizadas no ambiente de trabalho Linux e no sistema SeismicUn*x da Center for WavePhenomena (CWP) da Colorado Schoolof Mines. Os resultados das três técnicas foram comparados com o objetivo de mostrar a evolução da qualidade visual das seções resultantes através da continuidade dos eventos traço-a-traço e da relação sinal/ruído, analisar diferenças e melhoras nas seções migradas visando uma melhor interpretação geológica e promover melhores condições de processamento e imageamento, buscando auxiliar possíveis perfurações com sucesso
Acesso aberto (Open Access)
Structural constraints for image-based inversion methods
(Universidade Federal do Pará, 2016-04-22) MACIEL, Jonathas da Silva; COSTA, Jessé Carvalho; http://lattes.cnpq.br/7294174204296739
Esta tese apresenta duas metodologias de regularização estrutural para os métodos de análise de velocidade com migração e inversão conjunta com migração: regularização gradiente cruzado e filtragem com operadores morfológicos. Na análise de velocidade com migração, a regularização de gradiente cruzado tem como objetivo vincular os contrates de velocidade com o mapa de refletividade, através da paralelização dos vetores gradiente de velocidade com os vetores gradiente da imagem. Propõe-se uma versão com gradiente cruzado das funções objeto de minimização: Differential Semblance, Stack Power e Partial Stack Power. Combina-se a função Partial Stack Power com sua versão de gradiente cruzados, com o objetivo de aumentar gradativamente a resolução do modelo de velocidade, sem comprometer o ajuste das componentes de longo comprimento de onda do modelo de velocidade. Na inversão conjunta com migração propõe-se aplicar os operadores morfológicos de erosão e dilatação, no pré-condicionamento do modelo de velocidade em cada iteração. Os operadores usam o mapa de refletividade para delimitar as regiões com mesmo valor de propriedade física. Eles homogenizam a camada geológica e acentuam o contraste de velocidade nas bordas. Os vínculos estruturais não apenas irão reduzir a ambiguidade na estimativa do modelo de velocidade, mas tornará os métodos de inversão com migração mais estáveis, reduzindo artefatos, delineando soluções geologicamente plausíveis e acelerando a convergência da função objeto de minimização.

Navegar

Navegando Programa de Pós-Graduação em Geofísica - CPGF/IG por CNPq "CNPQ::CIENCIAS EXATAS E DA TERRA::GEOCIENCIAS::GEOFISICA::SISMOLOGIA"