Dissertações em Geofísica (Mestrado) - CPGF/IG

URI Permanente para esta coleçãohttps://repositorio.ufpa.br/handle/2011/4993

O Mestrado Acadêmico pertente a o Programa de Pós-Graduação em Geofísica (CPGF) do Instituto de Geociências (IG) da Universidade Federal do Pará (UFPA).

Navegar

Agora exibindo 1 - 20 de 28

Acesso aberto (Open Access)
Uma análise de métodos de supressão de múltiplas de superfície livre aplicados a um dado real
(Universidade Federal do Pará, 2011-10-14) OLIVEIRA, Andrei Gomes de; GOMES, Ellen de Nazaré Souza; http://lattes.cnpq.br/1134403454849660
O registro em levantamentos sísmicos 2D de eventos lineares são comuns nos traços dos sismogramas. Estes eventos se sobrepõem às reflexões de interesse (reflexões primárias) complicando as etapas do processamento sísmico e a identificação dos possíveis horizontes geológicos da subsuperfície. Pode-se destacar as múltiplas de superfície que se distinguem das demais por apresentarem baixas frequências, altas amplitudes e baixa velocidade de propagação. A identifificação e remoção de reflexões múltiplas em dados sísmicos reais é uma etapa de grande importância no processamento sísmicos, pois a atenuação destas múltiplas pode evitar erros na interpretação geológica das imagens sísmicas. O objetivo do trabalho é a determinação de um fluxo de processamento, para a atenuação de múltiplas de superfície livre. Métodos de eliminação de múltiplas de superfície (SRME, filltragem Radon, deconvolução preditiva no domínio τ - p e filtragem f-k) aplicados a um dado real 2D, foram utilizados em diferentes combinações. Estas são analisadas e são apresentados os resultados das melhores combinações.
Acesso aberto (Open Access)
Análise de um esquema de inversão para estimativa de anisotropia local a partir de dados de VSP
(Universidade Federal do Pará, 2012-05-18) BARRETO, Adriano César Rodrigues; GOMES, Ellen de Nazaré Souza; http://lattes.cnpq.br/1134403454849660
É estudado o esquema de inversão linear para a determinação de anisotropia em torno de um receptor situado no interior do poço a partir de dados de VSP multiazimutal. O dado consiste da componente vertical do vetor de vagarosidade e do vetor de polarização de ondas P diretas e refletidas medidos no receptor. A inversão é feita a partir de uma aproximação de primeira ordem, em torno de um meio isotrópico de referência arbitrário, da relação entre os dados de vagarosidade e polarização de onda P e os parâmetros fracamente anisotrópicos do meio (parâmetros WA). Esse esquema de inversão vale para meios com anisotropia arbitrária, desde que o grau de anisotropia seja fraca e independe do tipo de poço (vertical, horizontal ou inclinado).
Acesso aberto (Open Access)
Aplicação de modelos de substituição de fluido em rochas sedimentares oriundas do nordeste brasileiro
(Universidade Federal do Pará, 2015-06-06) TROVÃO, Ana Alzira Fayal; FIGUEIREDO, José Jadsom Sampaio de; http://lattes.cnpq.br/1610827269025210
Reservatórios carbonáticos correspondem a cerca de 50% dos reservatórios de hidrocarbonetos do planeta. Esse tipo de litologia apresenta diferentes formas de heterogeneidade, que são as principais causas de erros em sua caracterização, podendo induzir resultados errôneos na estimativa dos módulos elásticos das rochas em estado saturado. O presente trabalho objetiva uma análise comparativa dos modelos de substituição de fluidos em reservatório não-convencional do tipo carbonato. Especificamente, analisou-se os processos de substituição de fluido em amostras de rochas provindas do Nordeste brasileiro, sob condições laboratoriais controladas (temperatura, pressão e graus de saturação) sob as perspectivas petrofísicas e ultrassônicas através de modelos de substituição de fluidos convencionais (Gassmann, Biot) e não-convencionais (Brown e Korringa, Muller e Sahay). Foram analisadas seis (6) amostras de rochas carbonáticas e uma (1) rocha arenítica obtendo-se valores de permeabilidade, porosidade efetiva, densidade de grão; além de medidas elásticas de velocidades compressionais (Vp) e cisalhantes (Vs1 e Vs2), na condição de 100% gás, sendo substituída em seguida por água. Previsões teóricas feitas pelos modelos convencionais se aproximaram dos valores experimentais no caso de rochas consideradas homogêneas. Enquanto que predições feitas pelo modelos não convencionais (por exemplo, Muller & Sahay) se ajustaram aos resultados experimentais da maioria das tufas assim como nos calcários.
Acesso aberto (Open Access)
Aplicação de Well-Tie no imageamento de uma área da bacia Tumbes-progresso
(Universidade Federal do Pará, 2014-08-29) CHAGAS, Tiago Paulo Câncio das; SILVA, Carolina Barros da; http://lattes.cnpq.br/5306784916926352; GOMES, Ellen de Nazaré Souza; http://lattes.cnpq.br/1134403454849660
Regiões geológicas que apresentam grandes variações laterais de velocidade, como em regiões com falhas e domos de sal, são áreas complexas que, em geral, o imageamento sísmico em tempo é limitado. Nestes casos, o imageamento sísmico em profundidade pode apresentar melhores resultados. Entretanto, é necessário um modelo de velocidade em profundidade com precisão, o que não é uma tarefa fácil. Uma alternativa para a estimativa do modelo de velocidade em profundidade é oferecida através da técnica do Well-Tie. Com a calibração do traço sintético, obtido a partir de dado de poço, com a seção sísmica, o mapa de velocidade em profundidade pode ser obtido através da curva tempo-profundidade. Neste trabalho, utilizamos o Well-Tie para melhorar o modelo de velocidade RMS obtido do processamento do dado sísmico. O modelo de velocidade atualizado foi convertido para profundidade e usado para gerar uma imagem sísmica migrada em profundidade de uma região da Bacia Tumbes-Progresso.
Acesso aberto (Open Access)
Aproximações não-hiperbólicas do tempo de trânsito utilizando aproximantes de Padé
(Universidade Federal do Pará, 2017-07-03) NEVES, Rodolfo André Cardoso; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
As aproximações de tempo de trânsito são ferramentas indispensáveis para as etapas de empilhamento e migração de dados sísmicos. Com o intuito de aumentar a acurácia das aproximações de tempo de trânsito, propomos a utilização dos aproximantes de Padé no desenvolvimento de novas aproximações para o configuração ponto médio comum e superfície de reflexão comum (SRC). As aproximações hiperbólicas, sobretempo normal e superfície de rexlexão comum, são aproximações em série de Taylor de segunda ordem do tempo de trânsito de reflexão. Os aproximantes de Padé surgem como alternativa à série de Taylor, pois como característica principal possuem raio de convergência maior, melhorando a acurácia da aproximação. Estes aproximantes são produzidos através da própria série de Taylor da função aproximada. Essa nova aproximação é obtida através da aproximação de Padé [2/2] da equação generalizada do sobretempo; e da aproximação de Padé [2/2] das expansões em série de Taylor de quarta ordem para a superfície de reflexão comum. A acurácia das aproximações de Padé é superior as aproximações convencionais da literatura: sobretempo normal, hipérbole deslocada e aproximação para o modelo transversalmente isotrópico com eixo de simetria vertical (TIV). As aproximações de Padé para a configuração ponto médio comum dependem apenas de um parâmetro a mais do que a equação do sobretempo normal e mantém a acurácia para longos afastamentos. As aproximações não hiperbólicas para a configuração superfície de reflexão comum: SRC não hiperbólico, SRC quarta ordem e SRC Padé, possuem acurácia superior a aproximação SRC hiperbólico comumente utilizada pela indústria, e aumentam a região de convergência da aproximação no domínio do afastamento e da separação entre os pontos médios. A aproximação quadrática do SRC quarta ordem consegue ser superior inclusive a aproximação não hiperbólica do SRC, produzindo erros de aproximação consideravelmente menores na inversão dos parâmetros ótimos do SRC através do método dos mínimos quadrados.
Acesso aberto (Open Access)
Caracterização de hidratos de gás a partir de dados geofísicos da região da Bacia Tumbes-Progresso (Peru)
(Universidade Federal do Pará, 2014-02-25) TEJADA SOTO, Zoraida Roxana; GOMES, Ellen de Nazaré Souza; http://lattes.cnpq.br/1134403454849660
O hidrato de gás é conhecido como uma fonte potencial de energia de hidrocarbonetos. Entretanto, a detecção e exploração desse tipo de recurso é ainda inexpressiva. Isto se deve a dois fatores: a dificuldade na detecção de áreas rentáveis em hidratos de gás e devido ao alto risco exploratório. Nesse trabalho, é apresentada uma análise para a detecção de prováveis áreas de ocorrência do hidrato de gás em uma região da Bacia Tumbes-Progresso (Peru), utilizando dados geofísicos: magnético, gravimétrico, sísmicos e de poço. Apresenta-se três abordagens para evidenciar as possíveis áreas de hidrato de gás: na primeira, as anomalias magnéticas e gravimétricas são analisadas com relação a geologia e a área de estudo que foi escolhida. Em seguida foi feito o imageamento sísmico da região de interesse e a identificação dos BSRs (Bottom Simulating Refletor), refletores que delimitam zonas de hidrato de gás e gás livre. Aos potenciais BSR foi aplicada a análise AVO (Amplitude Versus Offset). Por fim, na terceira abordagem, são avaliados os dados de poço, calculando-se a porosidade e saturação da rocha na zona de estabilidade de hidrato de gás. De acordo com a análise apresentada, verificou-se a ocorrência de hidrato de gás na área de estudo.
Acesso aberto (Open Access)
Conversão tempo-profundidade de seções sísmicas empilhadas por raio imagem e raio normal
(Universidade Federal do Pará, 2014-06-24) REZENDE, Diogo Pena; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
A conversão de seções sísmicas tempo em profundidade vem sendo cada vez mais realizada na indústria do petróleo pelo método do traçado de raio imagem. Este método converte do tempo para a profundidade, ponto a ponto, as amplitudes da seção sísmica migrada no tempo. Para cada ponto da seção migrada é traçado um raio, perpendicular a superfície. Ao fim deste, a amplitude do ponto da seção migrada é colocada em profundidade. O método de migração sísmica pré ou pós-empilhamento consiste em colocar os eventos sísmicos nas posições mais próximas das reais em seções em tempo ou em profundidade. Seções sísmicas em profundidade fornecem uma imagem aproximada da subsuperficie, de modo a facilitar a identificação de possíveis estruturas geológicas acumuladoras de petróleo. A conversão de seções do domínio do tempo para o domínio da profundidade é considerada uma etapa intermediária do processo de construção de imagens sísmicas em profundidade. Neste trabalho é desenvolvido e testado o método de conversão tempo para profundidade de seções sísmicas afastamento-nulo. A construção de seções em profundidade é feita pelo traçado do raio normal. Este método utiliza como dado de entrada seções afastamento-nulo no domínio (tempo de interseção versus vagarosidade). Cada ponto neste domínio fornece as condições iniciais para o traçamento do raio normal: a posição inicial de partida dos raios e os ângulos iniciais formados com a normal à superfície, ou seja, os parâmetros de vagarosidades iniciais. Diferente do método do raio imagem, vários raios são traçados para um mesmo tempo de trânsito e uma mesma posição inicial, definindo uma curva isócrona. A amplitude de cada ponto da seção afastamento-nulo é convertida para a profundidade a partir da distribuição destas ao longo de cada isócrona em profundidade. A conversão por raio imagem tem boa recuperação da profundidade dos refletores com curvatura suave além de baixo custo computacional, pois apenas um raio é traçado para cada ponto da seção em tempo. Porém a continuidade dos refletores pode ser prejudicada no caso de refletoras com curvaturas acentuadas. Por sua vez, a conversão por raios normais recupera a profundidade e continuidade dos refletores de modo satisfatório, já que um mesmo ponto em tempo é convertido varias vezes em profundidade. Entretanto, possui um custo computacional mais alto, pois vários raios devem ser traçados para um mesmo ponto em tempo além de converter os artefatos inerentes do Slant Stack.
Acesso aberto (Open Access)
Determinação de modelo de velocidades para a migração em profundidade de dados sísmicos 2-D da Bacia do Jequitinhonha
(Universidade Federal do Pará, 2013-12-13) OLIVEIRA, Lelis Araújo de; GOMES, Ellen de Nazaré Souza; http://lattes.cnpq.br/1134403454849660
A migração no domínio do tempo é bastante empregado pela indústria de hidrocarbonetos por sua baixa demanda computacional, entretanto em áreas onde o campo de velocidade tem grande variação lateral, em decorrência do meio em subsuperfície apresentar estruturas complexas, como falhas, domos de sal, etc., esta técnica torna-se ineficiente na caracterização das feições geológicas. Embora o custo computacional seja alto, a migração no domínio da profundidade, mostra-se mais efetiva que a migração em tempo ao lidar com essas estruturas. Todavia para obter uma seção em profundidade há a necessidade de um modelo de velocidade em profundidade com certo grau de precisão. Neste trabalho, são aplicadas duas metodologias de conversão tempo-profundidade para a obtenção do modelo de velocidades em profundidade para duas linhas marinhas reais 2-D da Bacia do Jequitinhonha. A precisão dos modelos gerados foram analisados a partir das seções sísmicas migradas obtidas por meios dessas velocidades.
Acesso aberto (Open Access)
Experimental verification of Hudson and Eshelby-Chen’s effective crack theory
(Universidade Federal do Pará, 2015-12-04) HENRIQUES, Jéssica Pena; FIGUEIREDO, José Jadsom Sampaio de; http://lattes.cnpq.br/1610827269025210
A modelagem física em escala reduzida de meios fissurados/fraturados em laboratório tem servido como uma ótima alternativa para compreender o comportamento de meios anisotrópicos. Neste trabalho foram realizadas medidas ultrassônicas em amostras com baixas densidades e diferentes razões de aspecto de fissura. O objetivo principal se deu em investigar o comportamento dos parâmetros elásticos (velocidade das ondas, parâmetros de Thomsen ε e γ) e dos coeficientes do tensor de rigidez elástica para meios transversalmente isotrópicos. Comparar os resultados obtidos com as previsões feitas pelos modelos de Hudson (1981) e Eshelby-Cheng (1993) também foi investigado neste trabalho. Foram confeccionadas doze amostras com dois tipos de densidade de fissuras, 5 e 8%. As fissuras que possuem três razões de aspecto diferentes (0,133, 0,1778 e 0,2667) eram formadas por inclusões de borrachas em uma matriz isotrópica homogênea de resina. Além disso, uma matriz puramente isotrópica foi construída apenas por resina epóxi. Dentre todas as amostras, seis (três para cada densidade) possuem apenas um tipo de razão de aspecto (puras), enquanto outras seis (três para cada densidade) possuem três tipos de razões de aspectos diferentes (mistas). Entre as previsões dos modelos, o de Eshelby-Cheng (1993) mostra um melhor ajuste em relação aos resultados experimentais para as amostras puras (para as duas densidades de inclusões). No entanto, nenhum dos modelos prevê com mínima precisão a tendência para as amostras mistas.
Acesso aberto (Open Access)
Implementações alternativas de tomografia do tempo de trânsito utilizando a equação da onda
(Universidade Federal do Pará, 2014-04-23) CARDOSO FILHO, Josafat Lopes; COSTA, Jessé Carvalho; http://lattes.cnpq.br/7294174204296739
A tomografia do tempo de trânsito utilizando a equação acústica é uma alternativa robusta para se estimar modelos de velocidade com fortes variações. Aplicaçoes desta metodologia a dados sísmico inter poços, auxiliam no monitoramento e caracterização de reservatórios. A escolha da função objetivo, estratégias de precondicionamento do gradiente e funcionais regularizadores, quando da implementação deste tipo de tomografia, influênciam a robustez, eficiência e qualidade das estimativas do modelo de velocidade. Estes três aspectos da implementação da tomografia através da equação de onda são investigados. Duas propostas de funções objetivo são utilizadas neste trabalho; a primeira é sensível a diferenças de fase entre os pulsos e a outra é proposta para ser menos sensível a fase do pulso fonte. Ambas não necessitam de marcações dos eventos e se mostram muito robustas em experimentos numéricos utilizando modelos com forte variação de velocidade. Uma estratégia de precondicionamento do gradiente da função objetivo, adaptada da literatura em processamento de imagens, permitiu acelerar a convergência do algoritmo ao eliminar eventos espúrios causados pela inevitável abertura limitada da geometria de aquisição dos dados, ruídos aleatórios e efeitos causados pelas fontes e receptores. A adição de funcionais regularizadores penalizando o desvio do modelo de velocidade de informação a priori a partir de perfis de poços suavizados, contribui adicionalmente para a estimativa de um modelo de velocidade mais consistente e com maior resolução.
Acesso aberto (Open Access)
Interpolação de eventos de reflexão em traços sísmicos de dados pré-empilhamento usando aproximação de tempo de trânsito SRC-AF
(Universidade Federal do Pará, 2014-06-24) FERNANDES, Alexandre Sodré; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
O método sísmico é uma técnica bem conhecida por revelar detalhes das estruturas geológicas em subsuperfície, pois os dados sísmicos são processados e produzem seções sísmicas, as quais são imagens das estruturas em subsuperfície. A qualidade desses dados sísmicos registrados depende de fatores como a complexidade da subsuperfície, nível de ruído, topografia da superfície de aquisição, heterogeneidades no manto de intemperismo, entre outros. Aquisição irregular, afastamentos curtos, baixa cobertura nos pontos comuns em profundidade levam a uma baixa qualidade do imageamento em subsuperfície e baixa resolução nas seções sísmicas. No imageamento superfície de reflexão comum, a aproximação hiperbólica dos tempos de trânsito para raios paraxiais na vizinhança de um raio central refletido com afastamento fonte-receptor finito, em função de cinco atributos cinemáticos do campo de onda, representa o campo temporal de reflexão melhor se comparado aos métodos convencionais. A Superfície de Reflexão Comum para Afastamento Finito (SRCAF) aproxima eventos de reflexão no espaço na vizinhança do traço com um afastamento fixo que se deseja interpolar, somando eventos correlacionados, onde a saída é definida como a média ponderada das amplitudes ao longo da aproximação do tempo de trânsito SRC-AF, atribuindo o resultado para o respectivo traço para todos os pontos de uma seção AC. O objetivo especifico deste trabalho é modelar a regularização de seções temporais preenchendo áreas onde faltam dados sísmicos, aumentando a razão sinal ruido através da interpolação de eventos de reflexão em traços sísmicos baseado na aproximação do tempo de trânsito Superfície de Reflexão Comum para Afastamento Finito (SRC-AF). O algorítimo foi aplicado nas configurações de aquisição sísmica Fonte Comum, Afastamento Comum e Ponto médio Comum para um conjunto de dados sintéticos 2D, modelados por traçamento de raios.
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)
Migração Kirchhoff paraxial pré-empilhamento em profundidade com amplitudes verdadeiras
(Universidade Federal do Pará, 2015-08-31) MOREIRA, Marcio Fernando de Andrade; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
A migração sísmica em profundidade é um das etapas mais importantes no processo de exploração de áreas geologicamente complexas de interesse da indústria do petróleo. Para se obter uma imagem da estrutura mais realística, necessita-se de um ótimo processo de recuperação da configuração dos refletores em subsuperfície e eficiente determinação dos coeficientes de reflexão para caracterizar os reservatórios de petróleo e gás. Uma das técnicas mais utilizada para se obter a imagem sísmica é o método de migração Kirchhoff, que pode ser muito melhorado pela utilização de pesos adequados, os quais quando aplicados às amplitudes durante o empilhamento de difrações possibilita uma estimativa dos coeficientes de reflexão. Neste trabalho, foi utilizado o método de migração Kirchhoff paraxial em profundidade 2D. Na etapa de cálculo dos tempos de trânsito faz-se uso da aproximação paraxial de segunda ordem, através da qual se pode obter uma boa aproximação dos tempos de trânsito na malha de migração. Os pesos utilizados durante a migração foram calculados usando-se as equações dinâmicas do raio e pela extrapolação para os pontos da malha de migração. Em comparação com o método Kirchhoff convencional, os resultados se mostraram bastantes satisfatórios no que diz respeito ao aumento da resolução da imagem dos refletores, bem como na obtenção dos coeficientes de reflexão. A eficiência da técnica proposta foi testada em dados sintéticos de estrutura geológica tipo anticlinal e no dados sísmico Marmousi.
Acesso aberto (Open Access)
Migração Kirchhoff pré-empilhamento em profundidade usando aproximação paraxial do tempo de trânsito
(Universidade Federal do Pará, 2014-08-28) CUNHA, Antonio Rizimar de Andrade; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
Para obter-se imagens mais realísticas possíveis da subsuperfície, necessita-se de técnicas eficientes de avaliar o campo de onda. A literatura geofísica apresenta vários métodos baseados na solução numérica da equação da onda sísmica. Dentre às várias técnicas de imageamento sísmico, a Migração Kirchhoff Pré-Empilhamento em Profundidade continua sendo largamente utilizada devido a sua praticidade no tratamento de dados organizados em diversas configurações sísmica e na solução de problemas relacionados a imagens. Devido a isto, buscase criar alternativas mais eficientes para o cálculo dos parâmetros envolvidos nos processos de migração. Considera-se como ponto crítico o cálculo dos tempos de trânsito, fundamentais na focalização das amplitudes em seus respectivos pontos em profundidade. O método baseia-se na Teoria Paraxial do Raio utilizando o método Runge-Kutta de Quarta Ordem. A extrapolação paraxial dos tempos de trânsito pertence ao sistema de traçamento dinâmico de raios onde, através deste, é possível determinar informações em regiões complexas do meio geológico que, comumente, causam o aparecimento de zonas de sombras. A aplicação desta técnica exige quantidades previamente dispostas em um raio central de referência que é obtido pelo traçamento cinemático de raios. È necessário um macro modelo de velocidade do meio para o cálculo dos tempos de trânsito. Para este fim utilizamos modelos exatos de velocidades determinados a partir do software MATLAB. Para efeito de comparação, consideramos duas maneiras distintas de calcular as tabelas dos tempos de trânsito: A primeira foi realizada por meio da rotina RAYT2D do pacote SEISMIC UNIX (SU) que é considerado um método robusto de avaliação; na segunda, utilizou-se o Método Paraxial05. A realização da migração em profundidade requer como dados de entradas o dado sísmico e as tabelas dos tempos de trânsito. As imagens obtidas foram determinadas com um algoritmo escrito em SHELL e uma terceira imagem foi obtida por meio da subtração entre os resultados iniciais.
Acesso aberto (Open Access)
Modelagem sísmica acústica e elástica por diferenças finitas e imageamento do depósito de minério de ferro N4WS no estado do Pará
(Universidade Federal do Pará, 2017-07-03) HOLANDA, Rafael Mansano; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
O método das diferenças finitas consiste na resolução numérica de equações diferenciais e uma de suas aplicações é na solução das equações da onda ou elastodinâmica, pela substituição das derivadas parciais no espaço e no tempo por suas aproximações de diferenças finitas. O presente trabalho constituiu-se na realização da modelagem sísmica empregando o método das diferenças finitas, tanto para o caso acústico como elástico. Em seguida, realizou-se o imageamento sísmico, para o caso acústico, em um modelo que simula a área do depósito de minério de ferro de N4WS, localizado na Província mineral do Carajás, no estado do Pará. Para isso, foi considerado um modelo de velocidade heterogêneo, isotrópico e bidimensional para os casos acústico e elástico. Para gerar o sismograma sintético, foram utilizados os programas FDSKALAR, para o caso acústico, e SOFI2D, para o caso elástico, empregando aproximação de quarta ordem da equação da onda acústica e elástica no espaço e no tempo. Para a validação dos programas de diferenças finitas e auxílio na interpretação dos eventos presentes nos sismogramas, realizou-se a comparação com os tempos de trânsito obtidos pela teoria do raio. Após a modelagem sísmica, efetuou-se o processamento dos dados gerados, através do pacote de código aberto de softwares Seismic Unix, com o objetivo de gerar uma imagem do refletor presente no modelo de mineração. Desta maneira, procurou-se um melhor entendimento e controle sobre os problemas da modelagem sísmica e imageamento, contribuindo para a interpretação de dados sísmicos e compreensão sobre a propagação das ondas sísmicas em meios onde existem heterogeneidades significativas. Além disso, buscou-se demonstrar a aplicação dos métodos sísmicos no estudo e delineamento de corpos minerais.
Acesso aberto (Open Access)
NIP-tomografia usando método CRS e dados sísmicos marinhos
(Universidade Federal do Pará, 2013) AFONSO, João Batista Rodrigues; LEITE, Lourenildo Williame Barbosa; http://lattes.cnpq.br/8588738536047617
Este trabalho constou da aplicação das técnicas de processamento, inversão e imageamento em profundidade de dado sintético Marmousoft e do real Jequitinhonha obtido na plataforma continental oriental atlântica brasileira do Estado da Bahia. Foram aplicados os métodos de empilhamento NMO convencional e CRS, e a técnica de inversão NIP-tomográfica aos dadoscitados. Com o empilhamento NMO foram obtidos os mapas de distribuição de velocidades RMS e intervalar através da análise de velocidade no mapa de coerência semblance. Com os resultados obtidos do empilhamento CRS de ambos os dados foi realizada a marcação dos eventos de reflexão para a obtenção dos parâmetros do campo de onda, através de análise de medida de coerência nos dados sísmicos de multicobertura, que serviram de entrada para a inversão NIP-tomográfica para se obter um modelo suavizado de velocidade. Posteriormente foi obtido uma seção migrada em profundidade pelo método Kirchhoff. Os resultados aqui apresentados foram utilizados para analizar criticamente as técnicas estudadas, comparando o empilhamento CRS com o empilhamento NMO convencional. Além disso, avaliar a evolução da qualidade visual das seções resultantes CRS e NMO através da continuidade dos eventos traço-a-traço e da relação sinal/ruído, analisar diferenças e melhoras no modelo de velocidade obtido pela NIP-tomografia A migração Kirchhoff em profundidade pré-empilhamento foi aplicada para realizar interpretações geológicas e promover melhores condições de processamento e imageamento.
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 de dados sísmicos reais da região amazônica
(Universidade Federal do Pará, 2006-05-17) GOMES, Anderson Batista; LEITE, Lourenildo Williame Barbosa; http://lattes.cnpq.br/8588738536047617
O tratamento de dados sísmicos divide-se basicamente em três partes: o pré-processamento, o processamento e o imageamento. Na presente dissertação discutimos as etapas de pré-processamento e dois importantes métodos de processamento voltados à simulação de seções afastamento-nulo (AN) a partir de dados de cobertura múltipla. O processamento Convencional (NMO/DMO) e o processamento Superfície de Reflexão Comum (CRS) foram aplicados a dados sísmicos de reflexão de algumas linhas sísmicas do grupo 204 do conjunto de dados do Graben do Tacutu (Brasil). Usamos o Sistema CWP/SU para realizar as etapas de pré-processamento e de processamento convencional (NMO/DMO). O processamento CRS foi realizado com o Sistema WIT/CRS. As etapas de pré-processamento consistiram basicamente de três partes: organização da geometria; cancelamento e silenciamento de traços ruidosos; e filtragem na freqüência temporal (filtro em f) e na velocidade (filtro f-k). A etapa de deconvolução foi realizada, porém devido a um resultado que não trouxe contribuições, os resultados não tiveram uso posterior. Também, a correção estática de elevação não foi realizada porque a topografia é muito suave (variação topográfica menor do que 20 m) no platô do Tacutu. A qualidade dos resultados do processamento NMO/DMO foi fortemente comprometida devido à dependência do método do modelo de velocidade, que neste caso não foi suficientemente preciso. Além disso, encontramos dificuldades na análise de velocidade (AV) devido à grande quantidade de ruído presente nos dados. Como conseqüência, a correção de sobre-tempo normal (¿Normal Moveout¿, NMO) e a migração não geraram resultados melhores. Fundamentado nos atributos estimados pelo método de empilhamento CRS, foi obtido um macro-modelo de velocidades através de inversão tomográfica de reflexão. Usando este macro-modelo, foram realizadas migrações em profundidade pré- e pós-empilhamento. Também, os atributos CRS são usados no método correção estática residual, e os resultados demonstraram uma melhora na resolução das seções empilhadas. As seções resultantes dos empilhamentos e das migrações foram interpretadas visando o delineamento de estruturas. A partir dos detalhes visuais dos painéis, interpretamos afinamentos, um anticlinal principal falhado e desconformidades, e conjuntos de horsts e grabens foram traçados. Em contrapartida, o embasamento não pode ser facilmente traçado.
Acesso aberto (Open Access)
Processamento e imageamento de dados sísmicos marinhos
(Universidade Federal do Pará, 2010) LIMA JÚNIOR, Hamilton Monteiro de; LEITE, Lourenildo Williame Barbosa; http://lattes.cnpq.br/8588738536047617
O presente trabalho constou do processamento e imageamento de dados sísmicos marinhos reais obtidos na plataforma continental oriental atlântica brasileira da Bacia do Jequitinhonha. A etapa de pré-processamento foi realizada com parceria com o Grupo de Pesquisa em Sísmica da Universidade Federal da Bahia que aplicou o método de atenuação de múltiplas denominado SRME ao lado. Este dado foi submetido em seguida à aplicação de várias outras técnicas utilizando o pacote livre da CWP/SU. A etapa de processamento WIT/CRS constou da aplicação de três fluxogramas que utilizam a tecnologia CRS para simulação de seções ZO. Estes fluxogramas diferiram segundo a inclusão sistemática dos processos denominados Correção Residual Estática e Densificação do Dado Pré-empilhamento. Os resultados dos três fluxogramas foram comparados entre si 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. Além de seção CRS empilhada, é também obtida uma seção de migração CRS que tem a intenção de ser usada para interpretação geológica, com o objetivo de promover melhores condições quanto a indicação de uma possível perfuração com sucesso.