Teses em Geofísica (Doutorado) - CPGF/IG

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

O Doutorado 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 - 18 de 18

Acesso aberto (Open Access)
Análise de processos oceanográficos no estuário do rio Pará
(Universidade Federal do Pará, 2016-11-04) ROSÁRIO, Renan Peixoto; ROLLNIC, Marcelo; http://lattes.cnpq.br/6585442266149471
Esta pesquisa de doutorado investigou processos oceanográficos físicos no estuário do Rio Pará, com foco no processo de intrusão salina e hidrodinâmica. A escolha desse tema surgiu a partir da necessidade de se consolidar o entendimento dos aspectos hidrodinâmicos e hidrográficos no estuário do Rio Pará, já que esta região da Zona Costeira Amazônica ainda se apresenta como um desafio à pesquisa. Um dos desafios consistiu em definir métodos e parâmetros para resolver diferentes escalas de espaço e tempo. Neste contexto, observações diretas no ambiente estuarino foram realizadas através de medições de intensidade e direção de correntes, perfis verticais e longitudinais de salinidade e temperatura, durante um período considerado de baixa descarga fluvial e outro de alta descarga fluvial. Além disso, de forma inédita, foi realizado durante um ano e dez meses o monitoramento da salinidade e nível de água (maré) em pontos estratégicos do estuário. As principais conclusões que esta pesquisa obteve a partir desse conjunto de dados foi a identificação da intrusão salina no estuário do Rio Pará, adentrando cerca de 100 km da foz. A sensibilidade da frente de salinidade está sujeita a variabilidade sazonal, devido a descarga fluvial, e variabilidade diária, devido à grande energia das marés na região. O transporte de Stokes, gerado pela propagação da onda de maré no estuário foi a parcela responsável pelo transporte de sal estuário acima, intensificando essa intrusão salina. A porção mais interna do estuário (mais de 60 km da foz) não existe circulação gravitacional e o transporte de sal estuário acima é realizado totalmente por difusão turbulenta; e na porção externa o fluxo resultante reverte com a profundidade e os processos advectivo e difusivos são importantes para contribuir para o transporte de sal no estuário.
Acesso aberto (Open Access)
Determinação das velocidades intervalares usando a teoria paraxial do raio: aproximação de segunda ordem dos tempos de trânsito
(Universidade Federal do Pará, 1998) MONTES VIDES, Luis Alfredo; SÖLLNER, Walter Franz
Neste trabalho foi desenvolvido um método de solução ao problema inverso para modelos sísmicos compostos por camadas homogêneas e isotrópicas separadas por superfícies suaves, que determina as velocidades intervalares em profundidade e calcula a geometria das interfaces. O tempo de trânsito é expresso como uma função de parâmetros referidos a um sistema de coordenadas fixo no raio central, que é determinada numericamente na superfície superior do modelo. Essa função é posteriormente calculada na interface anterior que limita a camada não conhecida, através de um processo que determina a função característica em profundidade. A partir da função avaliada na interface anterior se calculam sua velocidade intervalar e a geometria da superfície posterior onde tem lugar a reflexão do raio. O procedimento se repete de uma forma recursiva nas camadas mais profundas obtendo assim a solução completa do modelo, não precisando em nenhum passo informação diferente à das camadas superiores. O método foi expresso num algoritmo e se desenvolveram programas de computador, os quais foram testados com dados sintéticos de modelos que representam feições estruturais comuns nas seções geológicas, fornecendo as velocidades em profundidade e permitindo a reconstrução das interfaces. Uma análise de sensibilidade sobre os programas mostrou que a determinação da função característica e a estimação das velocidades intervalares e geometria das interfaces são feitos por métodos considerados estáveis. O intervalo empírico de aplicabilidade das correções dinâmicas hiperbólicas foi tomado como uma estimativa da ordem de magnitude do intervalo válido para a aplicação do método.
Acesso aberto (Open Access)
Espalhamento elástico em meios anisotrópicos estratificados
(Universidade Federal do Pará, 1994-06-09) PROTÁZIO, João dos Santos; SCHOENBERG, Michael
A análise de AVO constitui-se, atualmente, numa importante ferramenta para a extração de informações litológicas a partir de dados sísmicos, através do uso dos contrastes de impedância acústica nas interfaces que separam diferentes litologias. A hipótese usual de isotropia deixa de valer, em muitos casos, após o advento de arranjos de grande afastamento e geofones com multi-superfície. Para a interpretação destes dados, a análise de AVO deve incluir anisotropia. Este trabalho apresenta uma teoria de AVO e resultados numéricos para um meio anisotrópico estratificado. Esta tese contém três contribuições. Inicialmente, é apresentada uma nova abordagem para o estudo da reflexão-transmissão através de interface plana que separam dois meios anisotrópicos com pelo menos um plano horizontal de simetria especular. As equações de Zoeppritz são generalizadas para incluir anisotropia, através da introdução das chamadas matrizes de impedância, o que simplifica bastante o formalismo anterior. Posteriormente, é descrito o estudo da reflexão de ondas P através de interface entre um meio isotrópico e outro transversalmente isotrópico (TI). É mostrado que a reflexão de ondas P, neste tipo de experimento, não fornece informações sobre a presença de anisotropia do semi-espaço TI, pelo menos em incidência pré-crítica. Finalmente, é discutido o comportamento da reflexão e transmissão de pulsos, em incidência pós-crítica, através de meios anisotrópicos estratificados. Observa-se que o comportamento pós-crítico dos pulsos espalhados carregam valiosa informações sobre a anisotropia dos meios atravessados por eles.
Acesso aberto (Open Access)
Imageamento homeomórfico de refletores sísmicos
(Universidade Federal do Pará, 1994-10-06) CRUZ, João Carlos Ribeiro; HUBRAL, Peter; http://lattes.cnpq.br/7703430139551941
Neste trabalho é apresentada uma nova técnica para a realização do empilhamento sísmico, aplicada ao problema do imageamento de refletores fixos em um meio bidimensional, suavemente heterogêneo, isotrópico, a partir de dados de reflexão. Esta nova técnica chamada de imageamento homeomórfico tem como base a aproximação geométrica do raio e propriedades topológicas dos refletores. São utilizados, portanto, os conceitos de frente de onda, ângulo de incidência, raio de curvatura da frente de onda, cáustica e definição da trajetória do raio; de tal modo que a imagem obtida mantém relações de homeomorfismo com o objeto que se deseja imagear. O empilhamento sísmico é feito, nesta nova técnica de imageamento, aplicando-se uma correção local do tempo, ∆ t, ao tempo de trânsito, t, do raio que parte da fonte sísmica localizada em xo, reflete-se em um ponto de reflexão, Co, sendo registrado como uma reflexão primária em um geofone localizado em xg, em relação ao tempo de referência to no sismograma, correspondente ao tempo de trânsito de um raio central. A fórmula utilizada nesta correção temporal tem como parâmetros o raio de curvatura Ro, o ângulo de emergência βo da frente de onda, no instante em que a mesma atinge a superfície de observação, e a velocidade vo considerada constante nas proximidades da linha sísmica. Considerando-se uma aproximação geométrica seguido um círculo para a frente de onda, pode-se estabelecer diferentes métodos de imageamento homeomórfico dependendo da configuração de processamento. Sendo assim tem-se: 1) Método Elemento de Fonte (Receptor) Comum (EF(R)C). Utiliza-se uma configuração onde se tem um conjunto de sismogramas relacionado com uma única fonte (receptor), e considera-se uma frente de onda real (de reflexão); 2) Método Elemento de Reflexão Comum (ERC). Utiliza-se uma configuração onde um conjunto de sismogramas é relacionado com um único ponto de reflexão, e considera-se uma frente de onda hipoteticamente originada neste ponto; 3) Método Elemento de Evoluta Comum (EEC). Utiliza-se uma configuração onde cada sismograma está relacionado com um par de fonte e geofone coincidentemente posicionados na linha sísmica, e considera-se uma frente de onda hipoteticamente originada no centro de curvatura do refletor. Em cada um desses métodos tem-se como resultados uma seção sísmica empilhada, u(xo, to); e outras duas seções denominadas de radiusgrama, Ro (xo, to), e angulograma, βo(xo, to), onde estão os valores de raios de curvatura e ângulos de emergência da frente de onda considerada no instante em que a mesma atinge a superfície de observação, respectivamente. No caso do método denominado elemento refletor comum (ERC), a seção sísmica resultante do empilhamento corresponde a seção afastamento nulo. Pode-se mostrar que o sinal sísmico não sofre efeitos de alongamento como consequência da correção temporal, nem tão pouco apresenta problemas de dispersão de pontos de reflexão como consequência da inclinação do refletor, ao contrário do que acontece com as técnicas de empilhamento que tem por base a correção NMO. Além disto, por não necessitar de um macro modelo de velocidades a técnica de imageamento homeomórfico, de um modo geral, pode também ser aplicada a modelos heterogêneos, sem perder o rigor em sua formulação. Aqui também são apresentados exemplos de aplicação dos métodos elemento de fonte comum (EFC) (KEYDAR, 1993), e elemento refletor comum (ERC) (STEENTOFT, 1993), ambos os casos com dados sintéticos. No primeiro caso, (EFC), onde o empilhamento é feito tendo como referência um raio central arbitrário, pode-se observar um alto nível de exatidão no imageamento obtido, além do que é dada uma interpretação para as seções de radiusgrama e angulograma, de modo a se caracterizar aspectos geométricos do model geofísico em questão. No segundo caso, (ERC), o método é aplicado a série de dados Marmousi, gerados pelo método das diferenças finitas, e o resultado é comparado com aquele obtido por métodos convecionais (NMO/DMO) aplicados aos mesmos dados. Como consequência, observa-se que através do método ERC pode-se melhor detectar a continuidade de refletores, enquanto que através dos métodos convencionais caracterizam-se melhor a ocorrência de difrações. Por sua vez, as seções de radiusgrama e angulograma, no método (ERC), apresentam um baixo poder de resolução nas regiões do modelo onde se tem um alto grau de complexidade das estruturas. Finalmente, apresenta-se uma formulação unificada que abrange os diferentes métodos de imageamento homeomórfico citados anteriormente, e também situações mais gerais onde a frente de onda não se aproxima a um círculo, mas a uma curva quadrática qualquer.
Acesso aberto (Open Access)
Influência de estruturas geológicas bidimensionais no campo geoeletromagnético na presença do eletrojato equatorial
(Universidade Federal do Pará, 2005) SILVA, Marcos Welby Correa; RIJO, Luiz; http://lattes.cnpq.br/3148365912720676
A Terra atua como um grande magneto esférico, cujo campo assemelha-se àquele gerado por um dipolo magnético. Este campo apresenta mudanças de intensidade que variam com a localização e a hora local. A parte principal do campo geomagnético se origina no interior da Terra através de processos eletromagnéticos. Extensivos estudos mostraram ainda que existem contribuições de origem externa ao planeta, principalmente de origem solar. Dentre estas fontes há anomalias do campo magnético que surgem a partir de um aumento diurno da corrente elétrica em uma estreita faixa da ionosfera, de direção leste-oeste, centrada no equador magnético e denominada Eletrojato Equatorial (EEJ). Ocasionalmente estas correntes podem apresentar reversões de fluxo, sendo denominadas Contra-Eletrojato (CEJ). Vários autores têm estudado os efeitos do EEJ e CEJ sobre as observações geoeletromagnéticas. Eles estão interessados no efeito combinado do EEJ e estruturas geológicas condutivas 1-D e 2-D. Nestes trabalhos a estrutura 2-D sempre se apresentava paralela ao eletrojato, o que é uma hipótese bastante restritiva ao se modelar ambientes geológicos mais realistas, em que corpos bidimensionais podem ter qualquer strike em relação ao EEJ. Neste trabalho apresentamos a solução deste problema sem esta restrição. Assim, mostramos os campos geoeletromagnéticos devidos a estruturas bidimensionais que possuam strike oblíquo em relação ao EEJ, através de perfis dos campos elétrico e magnéticos calculados na superfície e formando direção arbitrária à heterogeneidade condutiva 2-D. Com esta resposta avaliamos ainda qual a influência que estruturas bidimensionais exercem sobre a resposta magnetotelúrica, sob influência do Eletrojato Equatorial. Durante o desenvolvimento deste trabalho, utilizamos o método de elementos finitos, tendo por fonte eletromagnética o EEJ e o CEJ, que por sua vez foram representados por uma combinação de distribuições gaussianas de densidade de corrente. Estas fontes foram decompostas nas direções paralela e perpendicular à estrutura 2-D, resultando nos modos de propagação TE1 e TE2 e TM acoplados, respectivamente. Resolvemos o modo acoplado aplicando uma Transformada de Fourier nas equações de Maxwell e uma Transformada Inversa de Fourier na solução encontrada. De acordo com os experimentos numéricos realizados em um modelo interpretativo da Anomalia Condutiva da Bacia do Parnaíba, formado por uma enorme estrutura de 3000 ohm-m dentro de um corpo externo condutivo (1 ohm-m), concluímos que a presença do CEJ causa uma inversão na anomalia, se compararmos com o resultado do EEJ. Concluímos também que para as frequências mais altas as componentes do campo elétrico apresentam menor influência da parte interna do corpo 2-D do que da parte externa. Já para frequências mais baixas este comportamento se observa com as componentes do campo magnético. Com relação à frequência, vimos os efeitos do “skin-depth”, principalmente nas respostas magnéticas. Além disso, quando a estrutura 2-D está paralela ao eletrojato, o campo elétrico é insensível à estrutura interna do modelo para todos os valores de frequência utilizados. Com respeito ao ângulo θh entre a heterogeneidade e a fonte, vimos que o modo TM se manifesta naturalmente quando θh é diferente de 0°. Neste caso, o modo TE é composto por uma parte devido à componente da fonte paralela à heterogeneidade e a outra devido à componente da fonte perpendicular, que é acoplada ao modo TM. Assim, os campos calculados têm relação direta com o valor de θh. Analisando a influência do ângulo entre a direção do perfil dos campos e o strike da heterogeneidade verificamos que, à medida que θh se aproxima de 90°, os campos primários tornam-se variáveis para valores de θp diferentes de 90°. Estas variações causam uma assimetria na anomalia e dão uma idéia da inclinação da direção do perfil em relação aos corpos. Finalmente, concluímos que uma das influências que a distância entre o centro do EEJ e o centro da estrutura 2-D, causa sobre as componentes dos campos está relacionado às correntes reversas do EEJ e CEJ, pois a 500 km do centro da fonte estas correntes têm máxima intensidade. No entanto, com o aumento da distância, as anomalias diminuem de intensidade. Nas sondagens MT, nós também usamos o EEJ e o CEJ como fonte primária e comparamos nossos resultados com a resposta da onda plana. Deste modo observamos que as componentes do campo geoeletromagnético, usadas para calcular a impedância, têm influência do fator de acoplamento entre os modos TE2 e TM. Além disso, esta influência se torna maior em meios resistivos e nas frequências mais baixas. No entanto, o fator de acoplamento não afeta os dados magnetotelúricos em frequências maiores de 10-2 Hz. Para frequências da ordem de 10-4 Hz os dados MT apresentam duas fontes de perturbação: a primeira e mais evidente é devido à presença fonte 2-D (EEJ e CEJ), que viola a hipótese da onda plana no método MT; e a segunda é causada pelo acoplamento entre os modos TE2 e TM, pois quando a estrutura bidimensional está obliqua à fonte 2-D temos correntes elétricas adicionais ao longo da heterogeneidade. Concluimos assim, que o strike de uma grande estrutura condutiva bidimensional relativamente à direção do EEJ ou CEJ tem de fato influência sobre o campo geomagnético. Por outro lado, para estudos magnetotelúricos rasos (frequências maiores de 10-3 Hz) o efeito do ângulo entre a estrutura geológica 2-D e a direção do EEJ não é tão importante. Contudo, em estudos de litosfera frequências menores de 10-3 Hz) o acoplamento entre os modos TE2 e TM não pode ser ignorado.
Acesso aberto (Open Access)
Inversão da forma de onda orientada ao alvo
(Universidade Federal do Pará, 2016-09-16) COSTA, Carlos Alexandre Nascimento da; COSTA, Jessé Carvalho; http://lattes.cnpq.br/7294174204296739
Propomos uma nova metodologia de inversão da forma de onda orientada ao alvo para estimar os parâmetros físicos de uma área alvo em subsuperfície para dados sísmicos adquiridos com aquisição VSP-desviado ou com aquisição com fontes e receptores localizados na superfície. Além disso, investigamos a importância de eventos de múltiplos espalhamentos no conjunto de dados usados como inputs para estimar as respostas ao impulso da área alvo em subsuperfície iterativamente através de um esquema de inversão esparso para as mesma geometria de aquisição anteriormente citadas. Essas metodologias são baseadas no ajuste entre os campos de onda ascendente observado e modelado em um específico nível em profundidade próximo da área alvo, onde o campo de onda ascendente modelado é estimado através da representação tipo-convolução para função de Green. A principal característica de nossa metodologia de inversão da forma de onda orientada ao alvo é usar como inputs os campos de onda descendente e ascendente para estimar os parâmetros físicos locais sem necessitar que estes campos de onda sejam locais, dessa maneira evitamos adotar qualquer metodologia de redatumação para estimar as respostas ao impulso locais. Mostramos através de exemplos numéricos que podemos relaxar a necessidade de uma densa amostragem de fontes e receptores, característicos de esquemas de redatumação baseados em interferometria, para estimar as respostas ao impulso da área alvo através do esquema de inversão esparso usando como inputs os campos de onda com eventos de múltiplos espalhamentos. Estas metodologia são atraentes para dados sísmicos de aquisição VSP-desviado, pois para este tipo de dado não é necessário conhecer a região localizada acima da área alvo para estimar os campos de onda ascendente e descendente usados como inputs em ambas as metodologias. No entanto, para dados sísmicos adquiridos com fontes e receptores localizados na superfície, em princípio é necessário conhecer um modelo de velocidade a priori para estimar os campos de onda ascendente e descendente próximo da área alvo. Para ambas as metodologias investigadas, usamos a inversão da forma de onda baseada em migração para estimar os campos de onda ascendente e descendente próximo da área alvo a partir de dados sísmicos adquiridos com fontes e receptores na superfície.
Acesso aberto (Open Access)
Inversão de dados eletromagnéticos com o regularizador Variação Total e o uso da matriz de sensibilidade aproximada
(Universidade Federal do Pará, 2012-12-20) LUZ, Edelson da Cruz; RÉGIS, Cícero Roberto Teixeira; http://lattes.cnpq.br/7340569532034401
Neste trabalho apresentamos um estudo da aplicação do regularizador “Variação Total” (VT) na inversão de dados geofísicos eletromagnéticos. O regularizador VT reforça a proximidade entre os parâmetros adjacentes, mas, quando a influência de uma descontinuidade é sentida nos dados, este permite mudanças abruptas sobre os parâmetros. Isso faz com que o método seja uma alternativa válida, quando os dados observados usados na inversão provém de um ambiente geológico com uma distribuição suave de condutividade, mas que pode apresentar descontinuidades em lugares como as interfaces entre as camadas geoelétricas, como na margem de uma zona de óleo ou de um corpo de sal, que podem ser zonas muito resistivas no interior de sedimentos condutivos. Quando, devido a baixa resolução nos dados, o método não tem informações o suficiente para identificar a interface, o regularizador variação total reforça a proximidade entre os parâmetros adjacentes fazendo um transição suave entre as condutividades camadas, da mesma forma que é apresentado pela suavidade global. O método de Variação Total permite que modelos menos suaves sejam alcançados porque na norma L1 a medida de desajuste entre os pares de parâmetros adjacentes, dará o mesmo valor se a variação dos parâmetros é suave ou se a variação é abrupta, o que não é o caso se o mesmo desajuste é medido na norma L2, pois em uma distribuição suave a medida do desajuste é menor, sendo assim favorecida pela minimização desta norma. O uso deste regularizador permite uma melhor estimativa do tamanho de um corpo, seja ele resistivo ou condutivo. O trabalho está apresentado na forma de três artigos, cada um descrevendo uma etapa no desenvolvimento do problema da inversão, seguindo uma sequência de complexidade crescente no problema direto. O primeiro artigo neste trabalho é intitulado “Inversão de dados do CSEM marinho 1D de meio estratificado anisotrópico com o regularizador Variação Total”. Este descreve o passo inicial no desenvolvimento do problema: a inversão de dados do CSEM marinho de modelos estratificados 1D com anisotropia na condutividade das camadas. Este problema se presta bem para este desenvolvimento, porque tem solução computacional muito mais rápida do que o 2D, e nele já estão presentes as características principais dos dados do método CSEM marinho, como a largura muito grande da faixa de amplitudes medidas em um levantamento, e a baixa resolução, inerente às baixas frequências empregadas. A anisotropia acrescenta uma dificuldade a mais no problema, por aumentar o nível de ambiguidade nos dados e demandar ainda mais informação do que no caso puramente isotrópico. Os resultados mostram que a aplicação dos vínculos de igualdade do método VT permite a melhor identificação de uma camada alvo resistiva do que a simples aplicação dos vínculos tradicionais de suavidade. Até onde podemos aferir, esta solução se mostra superior a qualquer outra já publicada para este problema. Além de ter sido muito importante para o desenvolvimento de códigos em paralelo. O segundo artigo apresentado aqui, “Inversão de dados Magnetotelúricos com o regularizador Variação Total e o uso da matriz de sensibilidade aproximada”, trata da inversão de dados do método Magnetotelúrico em ambientes 2D. Este problema demanda um esforço computacional muito maior do que o primeiro. Nele, estudamos a aplicação do método dos estados adjuntos para gerar uma boa aproximação para as derivadas necessárias para a construção da matriz de sensibilidade usada na inversão. A construção da matriz de sensibilidade é a etapa que demanda mais tempo no processo de inversão, e o uso do método de estados adjuntos foi capaz de reduzir muito este tempo, gerando derivadas com um bom nível de aproximação. Esta etapa da pesquisa foi fundamental pelo problema direto ser matematicamente e computacionalmente muito mais simples do que o do CSEM marinho 2D. Novamente em comparação com a aplicação do regularizador de suavidade global, o regularizador de Variação Total permitiu, neste problema, uma melhor delimitação das bordas de heterogeneidades bidimensionais. A terceira parte deste trabalho, apresentada no artigo “Inversão de dados do CSEM marinho 2.5D com o regularizador Variação Total e o uso da matriz de sensibilidade aproximada”, apresenta a apliação do método de Variação Total ao problema da inversão de dados CSEM marinho 2.5D. Usamos o método dos estados adjuntos para gerar uma boa aproximação para as derivadas necessárias para a construção da matriz de sensibilidade usada na inversão, acelerando assim o processo de inversão. Para deixar o processo de inversão ainda mais rápido, lançamos mão da programação em paralelo com o uso de topologia. A comparação entre a aplicação do regularizador de suavidade global, e o regularizador de Variação Total permitiu, assim como nos casos anteriores, uma melhor delimitação das bordas de heterogeneidades bidimensionais.
Acesso aberto (Open Access)
Inversão de velocidades por otimização global usando a aproximação superfície de reflexão comum com afastamento finito
(Universidade Federal do Pará, 2016-08-25) MESQUITA, Marcelo Jorge Luz; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
A literatura geofísica recente tem mostrado que a construção de um modelo inicial mais realístico possível é a forma mais apropriada de se reduzir complicações do problema mal posto da inversão da forma da onda completa, e de fornecer as condições necessárias de convergência da função objetivo em direção ao mínimo global. Modelos otimizados são úteis como estimativas iniciais para métodos mais sofisticados de inversão e migração. Desenvolvo um método de inversão de velocidade da onda P usando dados sísmicos pré-empilhados para grandes afastamentos e baseado em medidas de coerência. A estratégia de inversão proposta é totalmente automática, baseada no cálculo do semblance e regida pela aproximação de tempo de trânsito paraxial, o chamado método da Superfície de Reflexão Comum com Afastamento Finito. Ela é realizada em dois passos, a princípio, usando raios imagens e o conhecimento de um modelo de velocidades conhecido a priori, determino as interfaces refletoras em profundidade a partir de seção migrada em tempo. A seguir, o modelo em profundidade gerado é usado como entrada na parametrização do modelo de velocidades, o qual é feito camada por camada. A estratégia de inversão é baseada na análise dos semblances calculados em cada ponto médio comum pelo método da Superfície de Reflexão Comum com Afastamento Finito. Para iniciar a inversão no segundo passo, a aproximação paraxial é feita pelo traçamento de raios no modelo de velocidades em profundidade obtido no primeiro passo. Usando a média aritmética dos semblances calculados em todos os pontos médios comuns como função objetivo, camada após camada, o algoritmo de otimização global Very Fast Simulated Annealing é aplicado para obter a convergência da função objetivo em direção ao máximo global. Ao aplicar a estratégia de inversão em dados sintéticos e reais, mostro a robustez do algoritmo de inversão proposto, produzindo modelos de velocidades da onda P otimizados a partir de dados pré-empilhados.
Acesso aberto (Open Access)
Modelagem e imageamento 2.5D no domínio do tempo através de diferenças finitas
(Universidade Federal do Pará, 2010) SILVA NETO, Francisco de Assis; SCHLEICHER, Maria Amélia Novais; http://lattes.cnpq.br/4767998352165705; COSTA, Jessé Carvalho; http://lattes.cnpq.br/7294174204296739
A modelagem 2.5D consiste em simular a propagação do campo de ondas em 3D em meios com simetria de translação em uma direção. Nesta tese esta abordagem é formulada para meios elásticos e anisotrópicos com classe de simetria arbitrária e a geometria de aquisição não precisa coincidir com um plano de simetria do meio. A migração por reversão no tempo do campo de ondas é formulada e implementada através de diferenças finitas 2.5D. Para reduzir os efeitos de retro-espalhamento e melhorar a recuperação da amplitude dos eventos migrados, propomos uma nova condição de imagem para migração reversa no tempo baseada na análise assintótica da condição de imagem clássica por correlação cruzada. Experimentos numéricos indicam que a migração reversa no tempo 2.5D com a nova condição de imagem proposta, melhora a resolução da imagem em relação à migração reversa no tempo 2D e reduz acentuadamente os ruídos causados por retro-espalhamento.
Acesso aberto (Open Access)
Modelagem eletromagnética 2.5-D de dados geofísicos através do método de diferenças finitas com malhas não-estruturadas
(Universidade Federal do Pará, 2014-10-23) MIRANDA, Diego da Costa; RÉGIS, Cícero Roberto Teixeira; http://lattes.cnpq.br/7340569532034401; HOWARD JUNIOR, Allen Quentin; http://lattes.cnpq.br/6447166738854045
Apresentamos a formulação eletromagnética em geometria 2.5-D aplicada à modelagem do marine controlled-source electromagnetic (mCSEM) através do método de Diferenças Finitas. Utilizamos a separação dos sinais primário e secundário para evitar problemas de singularidade devido à característica pontual da fonte eletro-magnética, o dipolo elétrico. As componentes do campo eletromagnético são derivadas dos resultados obtidos para os potenciais vetor magnético e escalar elétrico, calculados em todo o domínio do problema, o qual deve ser completamente discretizado para o uso do método de Diferenças Finitas. A limitação imposta pelo uso de malhas estruturadas no delineamento das geometrias presentes nos modelos geológicos, serviu como motivação para introduzirmos o uso de malhas não-estruturadas em nossos problemas. Essas malhas são completamente adaptáveis aos modelos com que trabalhamos, promovendo um delineamento suave de suas estruturas, e podendo ser localmente refinadas apenas nas regiões de interesse. Apresentamos também o desenvolvimento do método RBF-DQ, que faz uso da técnica de aproximação de funções por meio de combinações lineares das funções de base radial (RBF) e da técnica de quadraturas diferenciais (DQ) para a aproximação das derivadas de nosso problema diferencial. Nossos resultados mostraram que o uso do método de Diferenças Finitas com malhas-não estruturadas pode ser aplicado nos problemas de modelagem geofísica, promovendo uma melhoria na qualidade dos dados modelados quando comparados com os resultados obtidos através das técnicas tradicionais de Diferenças Finitas.
Acesso aberto (Open Access)
Modelagem numérica da influência do eletrojato equatorial em dados magnetotelúricos produzidos por estruturas tridimensionais
(Universidade Federal do Pará, 1993-02-05) GONZÁLEZ CARRASQUILLA, Antonio Abel; RIJO, Luiz; http://lattes.cnpq.br/3148365912720676
A América do Sul apresenta várias peculiaridades geomagnéticas, uma delas, é a presença do Eletrojato Equatorial, o qual se estende de leste para oeste no Brasil ao longo de aproximadamente 3500 km. Considerando-se o fato de que a influência do Eletrojato Equatorial pode ser detectada a grandes distâncias do seu centro, isto suscita o interesse em se estudar os seus efeitos na exploração magnetotelúrica no Brasil. A influência do eletrojato equatorial na prospecção magnetotelúrica tem sido modelada para meios geológicos uni e bidimensionais valendo-se para isto de soluções analíticas fechadas e de técnicas numéricas tais como elementos finitos e diferenças finitas. Em relação aos meios geológicos tridimensionais, eles tem sido modelados na forma de "camadas finas", usando o algoritmo "thin sheet". As fontes indutoras utilizadas para simular o eletrojato equatorial nestes trabalhos, tem sido linhas de corrente, eletrojatos gaussianos e eletrojatos ondulantes. Por outro lado, o objetivo principal da nossa tese foi o modelamento dos efeitos que o eletrojato equatorial provoca em estruturas tridimensionais próprias da geofísica da prospecção. Com tal finalidade, utilizamos o esquema numérico da equação integral, com as fontes indutoras antes mencionadas. De maneira similar aos trabalhos anteriores, os nossos resultados mostram que a influência do eletrojato equatorial somente acontece em frequências menores que 10-1 Hz. Este efeito decresce com a distância, mantendo-se até uns 3000 km do centro do eletrojato. Assim sendo, a presença de grandes picos nos perfis da resistividade aparente de um semi-espaço homogêneo, indica que a influência do eletrojato é notável neste tipo de meio. Estes picos se mostram com diferente magnitude para cada eletrojato simulado, sendo que a sua localização também muda de um eletrojato para outro. Entretanto, quando se utilizam modelos geo-elétricos unidimensionais mais de acordo com a realidade, tais como os meios estratificados, percebe-se que a resposta dos eletrojatos se amortece significativamente e não mostra muitas diferenças entre os diferentes tipos de eletrojato. Isto acontece por causa da dissipação da energia eletromagnética devido à presença da estratificação e de camadas condutivas. Dentro do intervalo de 3000 km, a resposta eletromagnética tridimensional pode ser deslocada para cima ou para baixo da resposta da onda plana, dependendo da localização do corpo, da frequência, do tipo de eletrojato e do meio geológico. Quando a resposta aparece deslocada para cima, existe um afastamento entre as sondagens uni e tridimensionais devidas ao eletrojato, assim como um alargamento da anomalia dos perfis que registra a presença da heterogeneidade tridimensional. Quando a resposta aparece deslocada para baixo, no entanto, há uma aproximação entre estes dois tipos de sondagens e um estreitamento da anomalia dos perfis. Por outro lado, a fase se mostra geralmente, de uma forma invertida em relação à resistividade aparente. Isto significa que quando uma sobe a outra desce, e vice-versa. Da mesma forma, comumente nas altas frequências as respostas uni e tridimensionais aparecem deslocadas, enquanto que nas baixas frequências se mostram com os mesmos valores, com exceção dos eletrojatos ondulantes com parâmetros de ondulação α = —2 e —3. Nossos resultados também mostram que características geométricas próprias das estruturas tridimensionais, tais como sua orientação em relação à direção do eletrojato e a dimensão da sua direção principal, afetam a resposta devido ao eletrojato em comparação com os resultados da onda plana. Desta forma, quando a estrutura tridimensional é rotacionada de 90°, em relação à direção do eletrojato e em torno do eixo z, existe uma troca de polarizações nas resistividades dos resultados, mas não existem mudanças nos valores da resistividade aparente no centro da estrutura. Ao redor da mesma, porém, se percebe facilmente alterações nos contornos dos mapas de resistividade aparente, ao serem comparadas com os mapas da estrutura na sua posição original. Isto se deve à persistência dos efeitos galvânicos no centro da estrutura e à presença de efeitos indutivos ao redor do corpo tridimensional. Ao alongar a direção principal da estrutura tridimensional, as sondagens magnetotelúricas vão se aproximando das sondagens das estruturas bidimensionais, principalmente na polarização XY. Mesmo assim, as respostas dos modelos testados estão muito longe de se considerar próximas das respostas de estruturas quase-bidimensionais. Porém, os efeitos do eletrojato em estruturas com direção principal alongada, são muito parecidos com aqueles presentes nas estruturas menores, considerando-se as diferenças entre as sondagens de ambos tipos de estruturas. Por outro lado, os mapas de resistividade aparente deste tipo de estrutura alongada, revelam um grande aumento nos extremos da estrutura, tanto para a onda plana como para o eletrojato. Este efeito é causado pelo acanalamento das correntes ao longo da direção principal da estrutura. O modelamento de estruturas geológicas da Bacia de Marajó confirma que os efeitos do eletrojato podem ser detetados em estruturas pequenas do tipo "horst" ou "graben", a grandes distâncias do centro do mesmo. Assim, os efeitos do eletrojato podem ser percebidos tanto nos meios estratificados como tridimensionais, em duas faixas de freqüência (nas proximidades de 10-1 Hz e para freqüências menores que 10-3 Hz), possivelmente influenciados pela presença do embasamento cristalino e a crosta inferior, respectivamente. Desta maneira, os resultados utilizando o eletrojato como fonte indutora, mostram que nas baixas freqüências as sondagens magnetotelúricas podem ser fortemente distorcidas, tanto pelos efeitos galvânicos da estrutura tridimensional como pela presença da influência do eletrojato. Conseqüêntemente, interpretações errôneas dos dados de campo podem ser cometidas, se não se corrigirem os efeitos do eletrojato equatorial ou, da mesma forma, não se utilisarem algoritmos tridimensionais para interpretar os dados, no lugar do usual modelo unidimensional de Tikhonov - Cagniard.
Acesso aberto (Open Access)
Modelagem sísmica e inversão na presença de anisotropia
(Universidade Federal do Pará, 1993-11-04) COSTA, Jessé Carvalho; SCHOENBERG, Michael
A necessidade da adoção de modelos elásticos anisotrópicos, no contexto da sísmica de exploração, vem crescendo com o advento de novas técnicas de aquisição de dados como VSP, walkway VSP, tomografia poço a poço e levantamentos sísmicos com grande afastamento. Meios anisotrópicos, no contexto da sísmica de exploração, são modelos efetivos para explicar a propagação de ondas através de meios que apresentam padrões de heterogeneidade em escala muito menor que o comprimento de onda das ondas sísmicas. Particularmente, estes modelos são muito úteis para explicar o dado sísmico mais robusto que são as medidas de tempo de trânsito. Neste trabalho, são investigados aspectos da propagação de ondas, traçado de raios e inversão de tempos de trânsito em meios anisotrópicos. É estudada a propagação de ondas SH em meios anisotrópicos estratificados na situação mais geral onde estas ondas podem ocorrer, ou seja, em meios monoclínicos com um plano vertical de simetria especular. É mostrado que o campo de ondas SH refletido a partir de um semi-espaço estratificado, não apresenta qualquer informação sobre a possível presença de anisotropia em subsuperfície. São apresentados métodos simples e eficientes para o traçado de raios em 3D através de meios anisotrópicos estratificados, baseados no princípio de Fermat. Estes métodos constituem o primeiro passo para o desenvolvimento de algoritmos de inversão de tempos de trânsito para meios anisotrópicos em 3D, a partir de dados de VSP e walkaway VSP. Esta abordagem é promissora para determinação de modelos de velocidade, que são necessários para migração de dados sísmicos 3D na presença de anisotropia. É efetuada a análise da inversão tomográfica não linear, para meios estratificados transversalmente isotrópicos com um eixo de simetria vertical(TIV). As limitações dos dados de tempo de trânsito de eventos qP para determinação das constantes elásticas, são estabelecidas e caracterizados os efeitos da falta de cobertura angular completa na inversão tomográfica. Um algoritmo de inversão foi desenvolvido e avaliado em dados sintéticos. A aplicação do algoritmo a dados reais demonstra a consistência de meios TIV. Esta abordagem é útil para casos onde há informação a priori sobre a estratificação quase plana das formações e onde os próprios dados do levantamento poço a poço apresentam um alto grau de simetria especular em relação a um plano vertical. Também pode ser útil em interpretações preliminares, onde a estimativa de um meio estratificado, serve como modelo de fundo para se efetuar análises mais detalhadas, por exemplo, como um modelo de velocidades anisotrópico para migração, ou como um modelo de calibração para análises de AVO.
Acesso aberto (Open Access)
Post-imaging analysis of pressure prediction in productive sedimentary basins for oil and gas exploration
(Universidade Federal do Pará, 2015-05-26) VIEIRA, Wildney Wallacy da Silva; LEITE, Lourenildo Williame Barbosa; http://lattes.cnpq.br/8588738536047617
Esta tese tem vários aspectos relacionados à modelagem de bacia sedimentar na exploração de óleo e gás, e com duas divisões gerais: estimativa de parâmetros, e predição de pressão. Para a estrutura do presente trabalho, o primeiro tópico está relacionada com a análise de velocidade e meios efetivos, onde se estima uma distribuição para a velocidade da onda P no tempo, seguido da transformação para a profundidade, e usar um modelo efetivo para a densidade e para a distribuição de velocidades da onda S. A razão para esta focalização inicialmente destas estimativas é porque eles representam a principal informação de base que se pode ter a partir do domínio sísmico, de onde os outros parâmetros sísmicos podem ser calculados, e que serve de base para a segunda parte deste trabalho. O segundo tópico está relacionado à cálculo de tensão, deformação e pressão na subsuperfície utilizando os dados das velocidades das ondas P e S e os modelos de densidade, com a finalidade de localizar áreas de altas e baixas pressões que atuam como bombas de sucção naturais para a mecânica da acumulação de óleo e gás em zonas produtivas e camadas reservatórios. Destacamos na segunda parte para a apresentação, chamar atenção para a sensibilidade do mapeamento de pressão em função da variação de velocidade e densidade. Classificamos a primeira divisão como dedicado ao processamento e imageamento sísmico convencional, e nomeamos a segunda divisão como predição de tensão-deformação-pressão pós-imageamento. Como o objetivo final da geofísica é obter imagens da subsuperfície sob diferentes propriedades, o cálculo de tensão só faz total sentido para o caso de dados reais, e isto faz com que os dados adquiridos seja obrigatoriamente em três componentes. Uma conclusão importante dos experimentos numéricos, mostramos que a pressão não tem um comportamento trivial, uma vez que pode diminuir com a profundidade e criar bombas naturais responsáveis pelo acúmulo de fluidos. A teoria de meios porosos baseia-se integralmente em geometria diferencial, porque esta disciplina matemática lida com propriedades geométricas coletivos para reservatórios reais. Mostrouse que tais propriedades coletivas são, nomeadamente, a porosidade, a área da superfície específica, a curvatura média e a curvatura Gaussiana. Por exemplo, meios fraturados tem, como regra, uma pequena porosidade, mas área da superfície específica muito grande, o que cria a razão 𝛾 = 𝑣𝑆/𝑣𝑃 anômala e alta, e isto significa um coeficiente de Poisson, 𝜎, negativo. Outra conclusão é relacionado ao cálculo da descontinuidade de pressão entre sólido e líquido, o que depende da estrutura de poros.
Acesso aberto (Open Access)
Região do espaço que mais influencia em medidas eletromagnéticas no domínio da frequência: caso de uma linha de corrente sobre um semi-espaço condutor
(Universidade Federal do Pará, 1994-07-28) BRITO, Licurgo Peixoto de; DIAS, Carlos Alberto; http://lattes.cnpq.br/9204009150155131
Localizar em subsuperfície a região que mais influencia nas medidas obtidas na superfície da Terra é um problema de grande relevância em qualquer área da Geofísica. Neste trabalho, é feito um estudo sobre a localização dessa região, denominada aqui zona principal, para métodos eletromagnéticos no domínio da freqüência, utilizando-se como fonte uma linha de corrente na superfície de um semi-espaço condutor. No modelo estudado, tem-se, no interior desse semi-espaço, uma heterogeneidade na forma de camada infinita, ou de prisma com seção reta quadrada e comprimento infinito, na direção da linha de corrente. A diferença entre a medida obtida sobre o semi-espaço contendo a heterogeneidade e aquela obtida sobre o semi-espaço homogêneo, depende, entre outros parâmetros, da localização da heterogeneidade em relação ao sistema transmissor-receptor. Portanto, mantidos constantes os demais parâmetros, existirá uma posição da heterogeneidade em que sua influência é máxima nas medidas obtidas. Como esta posição é dependente do contraste de condutividade, das dimensões da heterogeneidade e da freqüência da corrente no transmissor, fica caracterizada uma região e não apenas uma única posição em que a heterogeneidade produzirá a máxima influência nas medidas. Esta região foi denominada zona principal. Identificada a zona principal, torna-se possível localizar com precisão os corpos que, em subsuperfície, provocam as anomalias observadas. Trata-se geralmente de corpos condutores de interesse para algum fim determinado. A localização desses corpos na prospecção, além de facilitar a exploração, reduz os custos de produção. Para localizar a zona principal, foi definida uma função Detetabilidade (∆), capaz de medir a influência da heterogeneidade nas medidas. A função ∆ foi calculada para amplitude e fase das componentes tangencial (Hx) e normal (Hz) à superfície terrestre do campo magnético medido no receptor. Estudando os extremos da função ∆ sob variações de condutividade, tamanho e profundidade da heterogeneidade, em modelos unidimensionais e bidimensionais, foram obtidas as dimensões da zona principal, tanto lateralmente como em profundidade. Os campos eletromagnéticos em modelos unidimensionais foram obtidos de uma forma híbrida, resolvendo numericamente as integrais obtidas da formulação analítica. Para modelos bidimensionais, a solução foi obtida através da técnica de elementos finitos. Os valores máximos da função ∆, calculada para amplitude de Hx, mostraram-se os mais indicados para localizar a zona principal. A localização feita através desta grandeza apresentou-se mais estável do que através das demais, sob variação das propriedades físicas e dimensões geométricas, tanto dos modelos unidimensionais como dos bidimensionais. No caso da heterogeneidade condutora ser uma camada horizontal infinita (caso 1D), a profundidade do plano central dessa camada vem dada pela relação po = 0,17 δo, onde po é essa profundidade e δo o "skin depth" da onda plana (em um meio homogêneo de condutividade igual à do meio encaixante (σ1) e a freqüência dada pelo valor de w em que ocorre o máximo de ∆ calculada para a amplitude de Hx). No caso de uma heterogeneidade bidimensional (caso 2D), as coordenadas do eixo central da zona principal vem dadas por do = 0,77 r0 (sendo do a distância horizontal do eixo à fonte transmissora) e po = 0,36 δo (sendo po a profundidade do eixo central da zona principal), onde r0 é a distância transmissor-receptor e δo o "skin depth" da onda plana, nas mesmas condições já estipuladas no caso 1D. Conhecendo-se os valores de r0 e δo para os quais ocorre o máximo de ∆, calculado para a amplitude de Hx, pode-se determinar (do, po). Para localizar a zona principal (ou, equivalentemente, uma zona condutora anômala em subsuperfície), sugere-se um método que consiste em associar cada valor da função ∆ da amplitude de Hx a um ponto (d, p), gerado através das relações d = 0,77 r e p = 0,36 δ, para cada w, em todo o espectro de freqüências das medidas, em um dado conjunto de configurações transmissor-receptor. São, então, traçadas curvas de contorno com os isovalores de ∆ que vão convergir, na medida em que o valor de ∆ se aproxima do máximo, sobre a localização e as dimensões geométricas aproximadas da heterogeneidade (zona principal).
Acesso aberto (Open Access)
Respostas eletromagnéticas dos arranjos coplanar e coaxial em poço
(Universidade Federal do Pará, 2000-09-29) CARVALHO, Paulo Roberto de; VERMA, Om Prakash; http://lattes.cnpq.br/2723609019309173
Foi estudado a viabilidade de aplicação do arranjo coplanar de bobinas nas sondas de perfilagem em poço por indução eletromagnética. Paralelamente foram geradas as respostas do convencional arranjo coaxial, que é o amplamente utilizado nas sondas comerciais, com o propósito de elaborar uma análise comparativa. Através da solução analítica (meios homogêneos) e semi-analítica (meios heterogêneos) foram geradas inicialmente as respostas para modelos mais simples, tais como os do (1) meio homogêneo, isotrópico e ilimitado; (2) uma casca cilíndrica simulando a frente de invasão; (3) duas cascas cilíndricas para simular o efeito annulus; (4) uma interface plana e dois semi-espaços simulando o contato entre duas camadas espessas e (5) uma camada plano-horizontal e dois semi-espaços iguais. Apesar da simplicidade destes modelos, eles permitem uma análise detalhada dos efeitos que alguns parâmetros geoelétricos têm sobre as respostas. Aí então, aplicando ainda as condições de contorno nas fronteiras (Sommerfeld Boundary Value Problem), obtivemos as soluções semi-analíticas que nos permitiram simular as respostas em modelos relativamente mais complexos, tais como (1) zonas de transição gradacional nas frentes de invasão; (2) seqüências de camadas plano-paralelas horizontais e inclinadas; (3) seqüências laminadas que permitem simular meios anisotrópicos e (4) passagem gradacional entre duas camadas espessas. Concluimos que o arranjo coplanar de bobinas pode ser uma ferramenta auxiliar na (1) demarcação das interfaces de camadas espessas; (2) posicionamento dos reservatórios de pequenas espessuras; (3) avaliação de perfis de invasão e (4) localizar variações de condutividade azimutalmente.
Acesso aberto (Open Access)
Seismic amplitude analysis and quality factor estimation based on redatuming
(Universidade Federal do Pará, 2015-04-25) OLIVEIRA, Francisco de Souza; FIGUEIREDO, José Jadsom Sampaio de; http://lattes.cnpq.br/1610827269025210
A correção de amplitude é uma tarefa importante para corrigir a dissipação de energia sísmica por espalhamento geometrico ou atenuação durante a propagação da onda acústica / elástica em sólidos. Neste trabalho, propomos uma forma de estimar o fator de qualidade dos dados de reflexão sísmica, com uma metodologia baseada na combinação do método de deslocamento da frequência de pico (PFS) e do operador de redatumação. A contribuição deste trabalho está em corrigir os tempos de trânsito quando o meio é formado por muitas camadas. Em outras palavras, a correção da tabela de tempo de trânsito utilizada no método PFS é realizada utilizando um operador de redatumação. A operação proposta, é realizada de forma iterativa, com isto, permitindo estimar o fator de qualidade Q, camada por camada de um modo mais preciso. A operação de redatumação é usada para simular a aquisição de dados em novos níveis, evitando distorções produzidas por irregularidade próximas da superfície relacionadas com a geometria ou com as propriedades de heterogeneidade do meio. Propomos uma aplicação do operador de redatumação Kirchhoff em verdadeira ampilitude (TAKR) em meios homogêneos e comparamos com o operador de redatumação Kirchhoff convencional (KR) restrito ao caso de afastamento nulo. Nossa metodologia é baseada na combinação do método de deslocamento da frequência de pico e o operador de redatumação (TAKR com peso igual a 1). Aplicação em dados sintéticos e em dados reais sísmico (Viking Graben) e GPR(Siple Dome) demonstra a viabilidade de nossa análise.
Acesso aberto (Open Access)
Solução da equação de Archie com algoritmos inteligentes
(Universidade Federal do Pará, 2011) SILVA, Carolina Barros da; ANDRADE, André José Neves; http://lattes.cnpq.br/8388930487104926
A equação de Archie é um marco histórico da Avaliação de Formação por ser a primeira relação envolvendo as propriedades físicas das rochas e as suas propriedades petrofísicas possibilitando a identificação e a quantificação de hidrocarbonetos em subsuperfície. A saturação de água é a solução da equação de Archie obtida a partir da medida da resistividade e da estimativa da porosidade da formação. No entanto, a solução da equação de Archie é não trivial, na dependência do conhecimento prévio da resistividade da água de formação e dos expoentes de Archie (cimentação e saturação). Esta tese apresenta um conjunto de algoritmos inteligentes inéditos, que possibilitam a solução da equação de Archie. Uma variação da rede neural competitiva, denominada como rede neural bicompetitiva realiza o zoneamento do poço, delimitando as camadas reservatório. Para cada camada reservatório, um novo algoritmo genético, com uma estratégia evolutiva baseada na reprodução de fungos produz estimativas para os parâmetros de porosidade da matriz (densidade, tempo de trânsito e porosidade neutrônica), que aliados a um novo modelo de rocha produzem estimativas realistas da porosidade, considerando os efeitos da argilosidade. Uma nova rede neural competitiva denominada como rede competitiva angular realiza a interpretação do Gráfico de Pickett fornecendo as informações da resistividade da água de formação e do expoente de cimentação. Todos os resultados da metodologia aqui apresentada são obtidos com dados sintéticos e perfis convencionais.
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 Teses em Geofísica (Doutorado) - CPGF/IG por Agência de fomento "CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior"