Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-11T12:34:33Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoEfeitoJanela.pdf: 12983544 bytes, checksum: d8af62d9b160f107ec5594c633121e76 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:24:38Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-04T13:08:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoEfeitoJanela.pdf: 12983544 bytes, checksum: d8af62d9b160f107ec5594c633121e76 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T16:37:25Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoEfeitoJanela.pdf: 12983544 bytes, checksum: d8af62d9b160f107ec5594c633121e76 (MD5) / Made available in DSpace on 2014-09-10T16:37:25Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoEfeitoJanela.pdf: 12983544 bytes, checksum: d8af62d9b160f107ec5594c633121e76 (MD5) Previous issue date: 1999 / A presente dissertação consta de estudos sobre deconvolução sísmica, onde buscamos otimizar desempenhos na operação de suavização, na resolução da estimativa da distribuição dos coeficientes de reflexão e na recuperação do pulso-fonte. Os filtros estudados são monocanais, e as formulações consideram o sismograma como o resultado de um processo estocástico estacionário, e onde demonstramos os efeitos de janelas e de descoloração. O principio aplicado é o da minimização da variância dos desvios entre o valor obtido e o desejado, resultando no sistema de equações normais Wiener-Hopf cuja solução é o vetor dos coeficientes do filtro para ser aplicado numa convolução. O filtro de deconvolução ao impulso é desenhado considerando a distribuição dos coeficientes de reflexão como uma série branca. O operador comprime bem os eventos sísmicos a impulsos, e o seu inverso é uma boa aproximação do pulso-fonte. O janelamento e a descoloração melhoram o resultado deste filtro. O filtro de deconvolução aos impulsos é desenhado utilizando a distribuição dos coeficientes de reflexão. As propriedades estatísticas da distribuição dos coeficientes de reflexão tem efeito no operador e em seu desempenho. Janela na autocorrelação degrada a saída, e a melhora é obtida quando ela é aplicada no operador deconvolucional. A transformada de Hilbert não segue o princípio dos mínimos-quadrados, e produz bons resultados na recuperação do pulso-fonte sob a premissa de fase-mínima. O inverso do pulso-fonte recuperado comprime bem os eventos sísmicos a impulsos. Quando o traço contém ruído aditivo, os resultados obtidos com auxilio da transformada de Hilbert são melhores do que os obtidos com o filtro de deconvolução ao impulso. O filtro de suavização suprime ruído presente no traço sísmico em função da magnitude do parâmetro de descoloração utilizado. A utilização dos traços suavizados melhora o desempenho da deconvolução ao impulso. A descoloração dupla gera melhores resultados do que a descoloração simples. O filtro casado é obtido através da maximização de uma função sinal/ruído. Os resultados obtidos na estimativa da distribuição dos coeficientes de reflexão com o filtro casado possuem melhor resolução do que o filtro de suavização. / The present masters dissertation consists of studies on seismic deconvolution where we look for otimizing the operations of smoothing, of resolution for the estimation of distribution of reflection coefficients, and of recovery of source-pulse. The studied filters are single channel, and the formulations consider the seismic trace as the result of a stationary stochastic process, and we demonstrate the effects of taper windows and of prewhitening in resolution. The applied principle is the minimization of the difference's variance between real and desired outputs, resulting on a system Wiener-Hopf normal equations whose solution is the vector of filter coefficients to be applied in a convolution. Spike deconvolution is designed considering the distribution of reflection coefficients as a white series. The operator compresses the seismic events to impulses, and its inverse is a good approximation to the source-pulse. The application of taper windows and of prewhitening improve the output of this filter. Spike-series filters are designed using the distribution of reflection coefficients. The statistical properties of the reflection coefficients distribution affect the operator and its performance. Taper windows on autocorrelation degrades the output, and improvement is achieved when it is applied to the deconvolutional operator. The Hilbert transform produces good results in the recovery of source-pulse, under the premise of minimum-phase. The inverse of the recovered source-pulse compresses the seismic events to impulses. When the seismic trace contains additive noise, the results obtained with Hilbert transform are better than with spike deconvolution. Smoothing filter suppresses noise in the seismic trace as a function of a prewhitening parameter. The use of smoothed traces improves the spike deconvolution. Double prewhitening generates better results than the simple prewhitening. The matched filter operator is obtained from maximization of a signal/noise ratio function. Deconvolving the output of the matched filter for the estimation of the distribution of reflection coefficients possess better resolution than when using a smoothing filter.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-21T13:00:33Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoSismicoSuperficie.pdf: 34321354 bytes, checksum: 98bc579e7260095e946f7fa7acc5a244 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:11:52Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-05T13:23:59Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoSismicoSuperficie.pdf: 34321354 bytes, checksum: 98bc579e7260095e946f7fa7acc5a244 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-18T11:35:24Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoSismicoSuperficie.pdf: 34321354 bytes, checksum: 98bc579e7260095e946f7fa7acc5a244 (MD5) / Made available in DSpace on 2014-09-18T11:35:24Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoSismicoSuperficie.pdf: 34321354 bytes, checksum: 98bc579e7260095e946f7fa7acc5a244 (MD5) Previous issue date: 2001 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / ANP - Agência Nacional do Petróleo / O método de empilhamento sísmico por Superfície de Reflexão Comum (ou empilhamento SRC) produz a simulação de seções com afastamento nulo (NA) a partir dos dados de cobertura múltipla. Para meios 2D, o operador de empilhamento SRC depende de três parâmetros que são: o ângulo de emergência do raio central com fonte-receptor nulo (β₀), o raio de curvatura da onda ponto de incidência normal (R_NIP) e o raio de curvatura da onda normal (R_N). O problema crucial para a implementação do método de empilhamento SRC consiste na determinação, a partir dos dados sísmicos, dos três parâmetros ótimos associados a cada ponto de amostragem da seção AN a ser simulada. No presente trabalho foi desenvolvido uma nova sequência de processamento para a simulação de seções AN por meio do método de empilhamento SRC. Neste novo algoritmo, a determinação dos três parâmetros ótimos que definem o operador de empilhamento SRC é realizada em três etapas: na primeira etapa são estimados dois parâmetros (β°₀ e R°_NIP) por meio de uma busca global bidimensional nos dados de cobertura múltipla. Na segunda etapa é usado o valor de β°₀ estimado para determinar-se o terceiro parâmetro (R°_N) através de uma busca global unidimensional na seção AN resultante da primeira etapa. Em ambas etapas as buscas globais são realizadas aplicando o método de otimização Simulated Annealing (SA). Na terceira etapa são determinados os três parâmetros finais (β₀, R_NIP e R_N) através uma busca local tridimensional aplicando o método de otimização Variable Metric (VM) nos dados de cobertura múltipla. Nesta última etapa é usado o trio de parâmetros (β°₀, R°_NIP, R°_N) estimado nas duas etapas anteriores como aproximação inicial. Com o propósito de simular corretamente os eventos com mergulhos conflitantes, este novo algoritmo prevê a determinação de dois trios de parâmetros associados a pontos de amostragem da seção AN onde há intersecção de eventos. Em outras palavras, nos pontos da seção AN onde dois eventos sísmicos se cruzam são determinados dois trios de parâmetros SRC, os quais serão usados conjuntamente na simulação dos eventos com mergulhos conflitantes. Para avaliar a precisão e eficiência do novo algoritmo, este foi aplicado em dados sintéticos de dois modelos: um com interfaces contínuas e outro com uma interface descontinua. As seções AN simuladas têm elevada razão sinal-ruído e mostram uma clara definição dos eventos refletidos e difratados. A comparação das seções AN simuladas com as suas similares obtidas por modelamento direto mostra uma correta simulação de reflexões e difrações. Além disso, a comparação dos valores dos três parâmetros otimizados com os seus correspondentes valores exatos calculados por modelamento direto revela também um alto grau de precisão. Usando a aproximação hiperbólica dos tempos de trânsito, porém sob a condição de R_NIP = R_N, foi desenvolvido um novo algoritmo para a simulação de seções AN contendo predominantemente campos de ondas difratados. De forma similar ao algoritmo de empilhamento SRC, este algoritmo denominado empilhamento por Superfícies de Difração Comum (SDC) também usa os métodos de otimização SA e VM para determinar a dupla de parâmetros ótimos (β₀, R_NIP) que definem o melhor operador de empilhamento SDC. Na primeira etapa utiliza-se o método de otimização SA para determinar os parâmetros iniciais β°₀ e R°_NIP usando o operador de empilhamento com grande abertura. Na segunda etapa, usando os valores estimados de β°₀ e R°_NIP, são melhorados as estimativas do parâmetro R_NIP por meio da aplicação do algoritmo VM na seção AN resultante da primeira etapa. Na terceira etapa são determinados os melhores valores de β°₀ e R°_NIP por meio da aplicação do algoritmo VM nos dados de cobertura múltipla. Vale salientar que a aparente repetição de processos tem como efeito a atenuação progressiva dos eventos refletidos. A aplicação do algoritmo de empilhamento SDC em dados sintéticos contendo campos de ondas refletidos e difratados, produz como resultado principal uma seção AN simulada contendo eventos difratados claramente definidos. Como uma aplicação direta deste resultado na interpretação de dados sísmicos, a migração pós-empilhamento em profundidade da seção AN simulada produz uma seção com a localização correta dos pontos difratores associados às descontinuidades do modelo. / By using an arbitrary source-receiver configuration and without knowledge of the velocity model the recently introduced seismic data stacking method called Common Reflection Surface (CRS) simulates a zero-offset (ZO) section from multi-coverage seismic reflection data. For 2-D acquisition, as by-products provides three normal ray parameters: 1) the emergence angle (β₀); 2) the radius of curvature of the Normal Incidence Point Wave (R_NIP); and 3) the radius of curvature of the Normal Wave (R_N). The CRS stack is based on the hyperbolic traveltime paraxial approximation depending on β₀, R_NIP and R_N. In this thesis is presented a new algorithm of the CRS stack based on two-parameters and one-parameter search strategy combining global and local optimization methods for determine the three parameters that define the stacking surface (or operator). This is performed in three steps: 1) two-parameters search by applying global optimization to determine β₀ and R_NIP; 2) one-parameter search by applying global optimization to determine R_N; and 3) three-parameters search by applying local optimization to determine three parameters, using as initial approximations the parameter triple of the earlier two steps. In the first two steps is used the Simulated Annealing (SA) algorithm and the Variable Metric algorithm is used in the third step. To simulate the conflicting dip events, for each ZO sample where there are interference of intersecting events is determined an additional parameter triple corresponding to a local minimum. The stacking along the respective operator for each particular event allows to simulate their interference in the simulated ZO section by means of their superposition. This new CRS stack algoritm was applied to synthetic data sets providing high-quality simulated ZO sections and high precision determination of the stack parameters in comparison with the forward modeling. Using the hyperbolic traveltime approximation for identical radii of curvature R_NIP = R_N, an algorithm called Common Diffraction Surface (CDS) stack was developed to simulate ZO sections for diffracted waves. In a similar way to the CRS stack procedure, this new algorithm also uses the SA and VM optimization methods to determine the optimal parameter couple (β₀, R_NIP) that define the best CDS operator. The main features of the algorithm are the data normalization, common-offset data, large aperture of the CDS operator and positive search space for R_NIP. The application of the CDS stack algorithm in a synthetic dataset containing reflected and diffracted wavefields provides as main result a simulated ZO section containing diffracted events clearly defined. The post-stack depth migration of the ZO section locates correctly the discontinuities of the second interface.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-05T12:13:44Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseAproximacoesRpp.pdf: 4103377 bytes, checksum: 5e03cda58fdd1cf4609469596ddec5ad (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:07:51Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-05T13:52:00Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseAproximacoesRpp.pdf: 4103377 bytes, checksum: 5e03cda58fdd1cf4609469596ddec5ad (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-22T11:21:42Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseAproximacoesRpp.pdf: 4103377 bytes, checksum: 5e03cda58fdd1cf4609469596ddec5ad (MD5) / Made available in DSpace on 2014-09-22T11:21:42Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseAproximacoesRpp.pdf: 4103377 bytes, checksum: 5e03cda58fdd1cf4609469596ddec5ad (MD5) Previous issue date: 2002 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / ANP - Agência Nacional do Petróleo / FINEP - Financiadora de Estudos e Projetos / Neste trabalho são apresentadas aproximações lineares e quadráticas das equações de Zoeppritz para a obtenção dos coeficientes de reflexão e transmissão de eventos P-P e P-S em função dos ângulos de incidência e da média angular, bem como a análise de inversão linear de AVO, considerando os eventos de reflexão P-P e P-S dissociados e combinados. O uso das chamadas aproximações pseudo-quadráticas foi aplicado para obtenção de aproximações quadráticas apenas para eventos-PP, em torno dos contrastes médios das velocidades de ondas compressionais e cisalhantes e da razão Vs/Vp. Os resultados das aproximações desenvolvidas neste trabalho mostram que as aproximações quadráticas são mais precisas que as lineares, nas duas versões angulares. As comparações entre as aproximações em termos do ângulo de incidência e da média angular mostram que as aproximações quadráticas são equivalentes dentro do limite angular de [0º, 30º]. Por outro lado, as aproximações lineares em função do ângulo de incidência mostram-se mais precisas que as aproximações lineares em função da média angular. Na inversão linear, fez-se análises de sensibilidade e de ambigüidade e observou-se que, nos caso de eventos de reflexão P-P e P-S dissociados, apenas um parâmetro pode ser estimado e que a combinação destes eventos consegue estabilizar a inversão, permitindo a estimativa de dois dentre os parâmetros físicos dos meios (contrastes de impedância, de velocidade de onda P e de módulo de cisalhamento). / This work presents linear and quadratic approximations of the Zoeppritz equations for the derivation of reflection and transmission coefficients from P-P and P-S events as a function of incidence angles and angular average, as well as the linear inversion analysis, AVO, in respect to the disassociated and combined P-P and P-S reflection events. The use of the so-called pseudoquadratic approximations was applied for the derivation of quadratic approximations only for PPevents, around the average contrasts of compressional and shear waves velocities and Vs/Vp ratio. The results of approximations derived in this work show that quadratic approximations were more precise than the linear ones in the two angular versions. The comparisons between these approximations in terms incidence angle and angular average show that the quadratic approximations are equivalent within the angular limit of [0º to 30º]. In the other hand, the linear approximations as a function of incidence angle are more precise than the linear approximations as a function of the angular average. In the linear inversion, sensitivity and ambiguity analyses were carried out and one could see that in the case of disassociated reflection P-P and P-S events, just a parameter can be estimated and the combination of these events can stabilize the inversion permitting the estimation of two physical parameters for the media (impedance, P-wave velocity and shear bulk module contrasts).

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-07-09T16:11:24Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoMetodoEmpilhamento.pdf: 12404346 bytes, checksum: 35af7af06fdd3fc800f4767bf97f530e (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar as palavras-chave on 2014-08-07T13:38:31Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-21T13:06:48Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoMetodoEmpilhamento.pdf: 12404346 bytes, checksum: 35af7af06fdd3fc800f4767bf97f530e (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T13:27:33Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoMetodoEmpilhamento.pdf: 12404346 bytes, checksum: 35af7af06fdd3fc800f4767bf97f530e (MD5) / Made available in DSpace on 2014-08-22T13:27:33Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoMetodoEmpilhamento.pdf: 12404346 bytes, checksum: 35af7af06fdd3fc800f4767bf97f530e (MD5) Previous issue date: 2005 / O empilhamento sísmico é um método de imageamento utilizado para simulação de seções sísmicas com afastamento nulo (AN) a partir de dados de cobertura múltipla, produzindo um aumento da razão sinal/ruído na seção simulada. Este método é realizado usando o processamento convencional denominado ponto médio comum (PMC), através das chamadas correções sobretempo normal e sobretempo oblíquo (STN/STO). Alternativamente, um novo método de empilhamento sísmico, denominado superfície de reflexão comum (SRC) pode ser usado para simular seções sísmicas em afastamento nulo (AN). O método SRC tem como principais características: 1) a construção de uma aproximação do tempo de trânsito de reflexões primárias na vizinhança de um raio de reflexão normal; 2) esta aproximação de empilhamento tem como parâmetros os ângulos de emergência do raio normal e os raios de curvaturas de duas ondas hipotéticas denominadas onda ponto de incidência normal (onda PIN) e onda normal (onda N). O método pressupõe ainda que a velocidade da onda longitudinal na parte superior do modelo próxima à superfície da Terra seja conhecida. Para determinar adequadamente os parâmetros supracitados, fundamentais para o método de imageamento SRC, são necessárias estratégias de busca que utilizem casos particulares da aproximação de segunda ordem hiperbólica dos tempos de trânsito. As estratégias de busca apresentadas são: SRC pragmático-extendito e SRC global-local. Neste trabalho é apresentada uma revisão extensa dos fundamentos teóricos de método de empilhamento SRC, como também das estratégias de busca dos parâmetros SRC, dos algoritmos de implementação, e da delimitação de aberturas do operador de empilhamento SRC usando conceitos de zonas de Fresnel e o fator de estiramento da correção STN. Como um exemplo de aplicação dos atributos ou parâmetros SRC é apresentado o resultado da correção STN, sem estiramento do sinal sísmico, obtido utilizando as velocidades estimadas dos parâmetros SRC. Também, são apresentados os resultados da aplicação dos três métodos de empilhamento (PMC (STN/STO), SRC pragmático-extendido, SRC global-local) nos dados sintéticos Marmousi e real terrestre da bacia paleozóica Amazônica. Finalmente apresentam-se as discussões sobre os resultados, quanto à resolução da imagem sísmica e custo computacional. / Seismic stacking is a method designed to simulate zero-offset (ZO) seismic sections from multi-coverage seismic data. The resulting simulated ZO section gives rise to a significant increase of the signal to noise ratio. This method is done by means of seismic processing socalled the common mid point (CMP) method, through makes use of the so-called normal moveout and dip moveout (NMO/DMO) corrections. In this work, a new seismic stacking method called common reflection surface (CRS) is used to simulate ZO seismic sections. The main characteristic of the CRS method is: 1) the use of a more general, multi-parametric stacking operator that approximates the traveltimes of the primary reflections for sources and receiver pairs arbitrarily located on the vicinity of the normal ray. The parameters or attributes of the stacking operator are the radii of curvatures of two hypothetical waves called normal incidence point (NIP) wave and normal (N) wave; as well as the emergency angle of the normal ray. The CRS method assumes that the near-surface velocity is a priori known. To determine the abovementioned parameters appropriately, main for the CRS imaging method, is necessary search strategies that use cases special of the approach of second hyperbolic order of the traveltimes. The presented search strategies are: extended-pragmatic CRS and global-local CRS. To show the efficiency of these strategies they are applied in the synthetic Marmousi and real land datasets of the Amazon palaeozoic basin. As result sections ZO simulated by three different stacking methods (CMP, CRS extended- pragmatic and global-local CRS), to compare the efficiency of these, with relationship I cost her computational and the resolution of the seismic image.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-07-07T14:11:35Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EspalhamentoGeometricoModelos.pdf: 2612660 bytes, checksum: d364c2cff234f3468ec78391f90f4839 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Item sem assunto. on 2014-08-07T16:26:58Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-20T13:05:28Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EspalhamentoGeometricoModelos.pdf: 2612660 bytes, checksum: d364c2cff234f3468ec78391f90f4839 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T16:55:01Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EspalhamentoGeometricoModelos.pdf: 2612660 bytes, checksum: d364c2cff234f3468ec78391f90f4839 (MD5) / Made available in DSpace on 2014-08-22T16:55:01Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EspalhamentoGeometricoModelos.pdf: 2612660 bytes, checksum: d364c2cff234f3468ec78391f90f4839 (MD5) Previous issue date: 2004 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / A medição de parâmetros físicos de reservatórios se constitui de grande importância para a detecção de hidrocarbonetos. A obtenção destes parâmetros é realizado através de análise de amplitude com a determinação dos coeficientes de reflexão. Para isto, faz-se necessário a aplicação de técnicas especiais de processamento capazes de corrigir efeitos de divergência esférica. Um problema pode ser estabelecido através da seguinte questão: Qual o efeito relativamente mais importante como responsável pela atenuação de amplitudes, o espalhamento geométrico ou a perda por transmissividade? A justificativa desta pergunta reside em que a correção dinâmica teórica aplicada a dados reais visa exclusivamente o espalhamento geométrico. No entanto, a análise física do problema por diferentes direções põe a resposta em condições de dúvida, o que é interessante e contraditório com a prática. Uma resposta embasada mais fisicamente pode dar melhor subsídio a outros trabalhos em andamento. O presente trabalho visa o cálculo da divergência esférica segundo a teoria Newman-Gutenberg e corrigir sismogramas sintéticos calculados pelo método da refletividade. O modelo-teste é crostal para que se possa ter eventos de refração crítica além das reflexões e para, com isto, melhor orientar quanto à janela de aplicação da correção de divergência esférica o que resulta em obter o então denominado “verdadeiras amplitudes”. O meio simulado é formado por camadas plano-horizontais, homogêneas e isotrópicas. O método da refletividade é uma forma de solução da equação de onda para o referido modelo, o que torna possível um entendimento do problema em estudo. Para se chegar aos resultados obtidos foram calculados sismogramas sintéticos através do programa P-SV-SH desenvolvido por Sandmeier (1998), e curvas do espalhamento geométrico em função do tempo para o modelo estudado como descrito por Newman (1973). Demonstramos como uma das conclusões que a partir dos dados do modelo (velocidades, espessuras, densidades e profundidades) uma equação para a correção de espalhamento geométrico visando às “verdadeiras amplitudes” não é de fácil obtenção. O objetivo maior então deveria ser obter um painel da função de divergência esférica para corrigir as verdadeiras amplitudes. / The measurement of physical parameters of reservoirs is of great importance to the detection of hydrocarbons. To obtain these parameters, an amplitude analysis is performed with the determination of the reflection coefficients. For this, it is necessary the application of special processing techniques able to correct the spherical divergence effects on seismic time sections. A problem can be established through the following question: What is the relatively more important effect responsible for the amplitude attenuation: geometrical spreading or the loss by transmissivity? A justification for this question resides in that the theoretical dynamic correction applied to real data aims exclusively to the geometrical spreading. On the other side, a physical analysis of the problem by different directions places the answer in conditions of doubt, what is interesting and contradictory with the practice. A more physically based answer to this question can give better grounds to other works in progress. The present work aims at the calculus of the spherical divergence according to the Newman-Gutemberg theory, and to correct synthetic seismograms calculated by the reflectivity method. The test model considered is crustal in order to have critical refraction events besides reflection events, and to better position with respect to the time window for application of the spherical divergence correction, which results in obtaining the denoted “true amplitudes”. The simulated medium is formed by plane-horizontal, homogeneous and isotropic layers. The reflectivity method is a form of solution of the elastic wave equation for this reference model, what makes possible an understanding of the structured problem. To arrive at the obtained results, synthetic seismogram were calculated by using the fortran program P-SV-SH written and supplied by Sandmeier (1998), and reflection geometrical spreading curves as function of time were calculated as described by Newman (1973). As a conclusion, we have demonstrated that from the model information (velocities, thicknesses, densities and depths) it is not simple to obtain an equation for geometrical spreading correction aiming at the true amplitudes. The major aim would then be to obtain a panel of the spherical divergence function to correct for true amplitudes.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-07-07T16:04:37Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoComparativoEstereotomografia.pdf: 2594656 bytes, checksum: 282190512a1d4a327115bef3b5492988 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Item sem assunto. on 2014-08-07T16:27:22Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-08T12:51:48Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoComparativoEstereotomografia.pdf: 2594656 bytes, checksum: 282190512a1d4a327115bef3b5492988 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-18T11:35:38Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoComparativoEstereotomografia.pdf: 2594656 bytes, checksum: 282190512a1d4a327115bef3b5492988 (MD5) / Made available in DSpace on 2014-09-18T11:35:38Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoComparativoEstereotomografia.pdf: 2594656 bytes, checksum: 282190512a1d4a327115bef3b5492988 (MD5) Previous issue date: 2010 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / A determinação de um acurado modelo de velocidades é um requisito fundamental para a realização do imageamento sísmico. Métodos novos como a Estereotomografia préempilhamento e a Tomografia da onda NIP são ferramentas poderosas e bastante sugestivas para este propósito. Basicamente, a Estereotomografia pré-empilhamento se baseia no conceito de eventos localmente coerentes interpretados como reflexões primárias e associados com pares de segmentos de raios, que tem um mesmo ponto de reflexão em profundidade. Na Tomografia da onda NIP um evento sísmico é representado por uma onda hipotética NIP, que está relacionada a um ponto de reflexão em profundidade. Os atributos da onda NIP são determinados no decorrer do Empilhamento de Superfície de Reflexão Comum (empilhamento CRS). Este trabalho tem como objetivo, fazer um estudo comparativo de ambos os métodos de determinação do modelo de velocidades em profundidade. Assim, é realizada uma revisão dos fundamentos teóricos de ambos os métodos tomográficos, destacando as principais diferenças e aplicando estas aproximações em um dado sintético e um dado real marinho (linha sísmica 214-2660 da Bacia do Jequitinhonha). Para avaliar os modelos de velocidades encontrados pelas aproximações, foi utilizada a migração pré-empilhamento em profundidade do tipo Kirchhoff e também as famílias de imagem comum (CIG). Os resultados mostraram que ambos os métodos tomográficos fornecem modelos de velocidades representativos. Contudo, constatou-se que a estereotomografia tem melhor desempenho em meios com variações laterais de velocidades, porém, aplicável somente em dados pré-empilhados com alta razão sinal/ruído. / The determination of an accurate velocity model is a fundamental requirement for the seismic imaging. New methods, such as prestack stereotomography and poststack NIP wave tomography, are powerful and very suggestive tools for this task. The prestack stereotomography is basically based on the concept of locally coherent events interpreted as primary reflections and that are associated to ray segments that are linked through the same reflection point in depth. In NIP wave tomography a seismic event is represented by a hypothetic NIP wave that is associated to a reflection point in depth. The NIP wave attributes are determined during Common Reflection Surface (CRS) procedure. The objective of this work is to compare both methods of velocity model determination in depth. Then a review of the theoretical foundations of both tomographic methods are made, considering its main differences, and then applied to a synthetic data and a real marine dataset (seismic line 214-2660 of the Jequitinhonha Basin, Brazil). In order to evaluate the velocity models determined by these two approximations, the data were prestack depth migrated using the Kirchhoff algorithm and also generated Common Image Gathers (CIG). The results have shown that both tomographic methods yield representative velocity models. However, it was noticed that the velocity model estimated by stereotomography behaved better in laterally varying media, but only applied in prestack data with a high signal-to-noise ratio.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-05T16:07:43Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_IdentificacaoReflexoesMultiplas.pdf: 1217629 bytes, checksum: 141af446a98307cf3a885b2c7606e5f7 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-06T13:21:46Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-19T14:38:30Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_IdentificacaoReflexoesMultiplas.pdf: 1217629 bytes, checksum: 141af446a98307cf3a885b2c7606e5f7 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-20T13:11:07Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_IdentificacaoReflexoesMultiplas.pdf: 1217629 bytes, checksum: 141af446a98307cf3a885b2c7606e5f7 (MD5) / Made available in DSpace on 2014-08-20T13:11:07Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_IdentificacaoReflexoesMultiplas.pdf: 1217629 bytes, checksum: 141af446a98307cf3a885b2c7606e5f7 (MD5) Previous issue date: 2005 / As reflexões múltiplas presentes nos sismogramas ocultam informações importantes sobre os refletores em subsuperfície e, podem até tornar completamente invisíveis as reflexões primárias, como no caso dos sismogramas marinhos, que muitas das vezes, exibem uma aparência anelar com fortes superposições das reflexões múltiplas, sobre as reflexões primárias. Problema este que tem sido alvo de importantes pesquisas, com o intuito de identificar, atenuar e/ou eliminá-las, através de vários métodos populares. O objetivo principal deste trabalho é a identificação das reflexões múltiplas. Com essa finalidade foi gerada, por modelamento direto, uma seção sísmica, com afastamento nulo (AN), contendo reflexões primárias e múltiplas simétricas de primeira ordem. Posteriormente, foi aplicada a migração cinemática do tipo Kirchhoff para obter o modelo em profundidade, apresentando uma boa recuperação dos refletores, bem como a presença de um refletor fictício, quando comparado com a seção anteriormente especificada. Foi obtida uma seção sísmica AN, do modelo migrado, na qual não é observado o segundo refletor, devido à ausência de contraste de impedância, entre a segunda e terceira camada, sendo este o primeiro indício de que o refletor fictício deste modelo é uma múltipla. Outro indício sobre a existência da múltipla foi a simetria encontrada entre as curvaturas do primeiro e terceiro refletor. Finalmente, foram calculados os parâmetros das frentes de ondas Hipotéticas Ponto de Incidência Normal (PIN) e Normal (N), bem como a velocidade Normal Moveout (NMO), tanto para os eventos de reflexões primárias como os eventos de reflexões múltiplas, para o modelo direto e para o modelo migrado. Em seguida, foram realizadas as comparações destes parâmetros, o que permitiu confirmar a veracidade dos indícios anteriores para a identificação das reflexões múltiplas. / The multiple reflections present in seismograms hide important information of reflectors in subsurface,sometimes turning primary reflections invisible. This is the case of marine seismograms, in which sometimes they present a ring-like appearance, with strong superposition of multiple reflections, together with primary reflections. This problem has been the subject of important researches, whose principal objective is the identification, attenuation and elimination of multiples, using in this respect some widespread methods. The objective of this work is the identification of multiple reflections. With this purpose, we have forward modeled zero-offset (ZO) seismic sections containing primary reflections, together with first-order multiple reflections. Later on,we have performed a kinematic Kirchhoff depth migration in order to obtain a simulated geological model in depth, with true and fictitious reflectors at the same time. From this depth migrated section, we have modeled a ZO section, in which the fictitious reflector is not recorded, due to the impedance contrast in the second and third layers of the model. This is the first clue that the referred fictitious reflector is, in fact, the result of imaging a multiple reflection. Another clue for this assertion is the semblance of the curvatures of the first and third reflectors. Finally, the wavefront parameters of two hypothetical waves Normal Incidence Point (NIP) and Normal (N) were computed, as well as the NMO velocity for the primary and multiple events, in the forward and migrated models. We have finally compared these parameters in order to verify the later clues with respect to the identification of multiple reflections.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-18T16:54:07Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ImageamentoMultifocoRefletores.pdf: 15149810 bytes, checksum: b0d10801d9480960ffdf7f59a15484e1 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-18T17:14:34Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ImageamentoMultifocoRefletores.pdf: 15149810 bytes, checksum: b0d10801d9480960ffdf7f59a15484e1 (MD5) / Made available in DSpace on 2014-09-18T17:14:34Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ImageamentoMultifocoRefletores.pdf: 15149810 bytes, checksum: b0d10801d9480960ffdf7f59a15484e1 (MD5) Previous issue date: 2000 / DAAD - Serviço Alemão de Intercâmbio Acadêmico / Deutscher Akademischer Austauschdiens / A simulação de uma seção sísmica de afastamento-nulo (ZO) a partir de dados de cobertura múltipla para um meio 2-D, através do empilhamento, é um método de imageamento de reflexão sísmica muito utilizado, que permite reduzir a quantidade de dados e melhorar a relação sinal/ruído. Segundo Berkovitch et al. (1999) o método Multifoco está baseado na Teoria do Imageamento Homeomórfico e consiste em empilhar dados de cobertura múltipla com distribuição fonte-receptor arbitrária de acordo com uma nova correção de sobretempo, chamada Multifoco. Esta correção de sobretempo esta baseada numa aproximação esférica local da frente de onda focalizante na vizinhança da superfície da terra. Este método permite construir uma seção sísmica no domínio do tempo de afastamento nulo aumentando a relação sinal/ruído. A técnica Multifoco não necessita do conhecimento a priori de um macro-modelo de velocidades. Três parâmetros são usados para descrever a aproximação de tempo de trânsito, Multifoco, os quais são: 1) o ângulo de emergência do raio de afastamento nulo ou raio de reflexão normal (β0), 2) a curvatura da frente de onda no Ponto de Incidência Normal (RNIP) e 3) curvatura da frente de Onda Normal (RN). Sendo também necessário a velocidade próximo a superfície da terra. Neste trabalho de tese aplico esta técnica de empilhamento Multifoco para dados de cobertura múltipla referidos a modelos de velocidade constante e modelo heterogêneos, com o objetivo de simular seções sísmicas afastamento-nulo. Neste caso, como se trata da solução de um problema direto, o macro-modelo de velocidades é considerado conhecido a priori. No contexto do problema inverso tem-se que os parâmetros RNIP, RN e β0 podem ser determinados a partir da análise de coerência aplicada aos dados sísmicos de múltipla cobertura. Na solução deste problema a função objetivo, a ser otimizada, é definida pelo cálculo da máxima coerência existente entre os dados na superfície de empilhamento sísmico. Neste trabalho de tese nos discutimos a sensibilidade da aproximação do tempo de trânsito usado no empilhamento Multifoco, como uma função dos parâmetros RNIP, RN e β0. Esta análise de sensibilidade é feita de três diferentes modos: 1) a primeira derivada da função objetivo, 2) a medida de coerência, denominada semelhança, e 3) a sensibilidade no Empilhamento Multifoco. / The simulation of a zero-offset section (ZO) from multi-coverage seismic reflection data for a 2-D media, through the stack, is a widely used seismic reflection imaging method, that allows to reduce the amount of data and to improve the signal-to-noise ratio. According to Berkovitch et al. (1999) the Multifocusing method is based on Theory of Homeomorphic Imaging and consists of stacking multi-coverage seismic reflection data with arbitrary distribution source-receiver according to a new local moveout correction, called Multifocusing. This moveout correction is based on a local spherical approximation of the focusing wavefront in the vicinity of the surface of the earth. This method allows to build a seismic section in the domain of the time of a zero-offset increasing the signal-to-noise ratio. The Multifocusing technique does not need any knowledge a priori of the macro-model velocity. Three parameters are used to describe the Multifocusing moveout correction, which are: 1) emergence angle of the zero-offset ray or normal reflection ray (β0), 2) the wavefront curvature at the Point of Normal Incidence (RNIP and 3) the wavefront curvature of Normal Wave (RN). Being also necessary the near-surface velocity. In this thesis work I apply this Multifocusing Stack technique for multi-coverage seismic reflection data covering the models of constant velocity and heterogeneous model, with the objective of simulating of zero-offset sections. In this case as it is the solution of forward problem, this macro velociy model is considered know apriori. In the context of the inverse problem it is had the parameters, RNIP, RN and β0 and can be determining through the analysis of applied coherence to the multi-coverage seismic reflection data. In the solution of this problem the objective function, to be optimized, is defined by the calculation of the maximum existent coherence among the data in the surface of seismic stack. In this thesis we discuss the sensibility of the travetime, used in Multifocusing Stack, as a function of the parameters RNIP, RN and β0. This sensibility analysis is done of three different manners: 1) the first derivate of the objetive function, 2) the coherence measure, denominated semblance, and 3) the sensibility in the Multifocusing Stack.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-10-20T13:51:09Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MetodosOtimizacaoGlobal.pdf: 2631790 bytes, checksum: d0343a988ffbcc231d9547a1fa5255c5 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-10-22T13:40:43Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MetodosOtimizacaoGlobal.pdf: 2631790 bytes, checksum: d0343a988ffbcc231d9547a1fa5255c5 (MD5) / Made available in DSpace on 2014-10-22T13:40:43Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MetodosOtimizacaoGlobal.pdf: 2631790 bytes, checksum: d0343a988ffbcc231d9547a1fa5255c5 (MD5) Previous issue date: 2004 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / 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. / The Common Reflection Surface Stack (CRS) method provides the simulation of zerooffset (ZO) sections by means of the summing seismic events of the multicoverage data contained in the stacking surface. This method does not depend on the velocity macro-model of medium; it only requires a priori knowledge of the near-surface velocity. The simulation of ZO sections for this stacking method uses a hyperbolic second-order traveltime approximation of the paraxial rays to define the CRS stacking surface or CRS stack operator. For 2D media, this operator depends on three kinematic attributes of two hypothetical waves (NIP and N waves), observed in the point of emergency of the central ray with normal incidence, namely, the angle of emergency of the central ZO ray (β0), the radius of curvature of the Normal Incidence Point Wave (RNIP) and the radius of curvature of the Normal Wave (RN). Therefore, the optimization problem in the CRS method consists in the determination, from the seismic data, of the three optimal parameters (β0, RNIP, RN) associated to each sample point of ZO section to be simulated. The simultaneous determination of these parameters can be made by means of multidimensional global search process (or global optimization), using as objective function some coherence criterion. The optimization problem in CRS method is very important for the good performance with respect to quality of the results and mainly to computational cost, compared with the methods traditionally used in the seismic industry. There are several search strategies to determine these parameters, based on systematic searches and using optimization algorithms, where only one parameter at each time can be estimated, or the two or three parameters simultaneously. Taking in to account the search strategy by means of the application of global optimization, these three parameters can be estimated through of procedures: in the first case the three parameters can be simultaneously estimated and in second case initially two parameters can be determined simultaneously (β0, RNIP), and subsequently the third parameter (RN), using the values of the two parameters already known. In this work it is presented the application and comparison of four algorithms of global optimization to find the CRS optimal parameters: Simulated Annealing (SA), Very Fast Simulated Annealing (VFSA), Differential Evolution (DE) and Controlled Random Search – 2 (CRS2). As importants results of the application of each optimization method, as well as between the methods regarding the effectiveness, efficiency and reliability to determine the best CRS parameters are presented. Subsequently, applying the global search strategies for the determination of these parameters, by means of the optimization method VFSA that presented the best performance, the CRS stacking was applied to the Marmousi dataset, one stacking using two parameters (β0, RNIP), estimated by global search, and another CRS stacking using the three parameters (β0, RNIP, RN), also estimated by global search.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-10-20T15:48:49Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoAmplitudeVerdadeira.pdf: 9758204 bytes, checksum: 62046151ddeb3a3200cca375f8b96f32 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-10-22T13:46:36Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoAmplitudeVerdadeira.pdf: 9758204 bytes, checksum: 62046151ddeb3a3200cca375f8b96f32 (MD5) / Made available in DSpace on 2014-10-22T13:46:36Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoAmplitudeVerdadeira.pdf: 9758204 bytes, checksum: 62046151ddeb3a3200cca375f8b96f32 (MD5) Previous issue date: 1999 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / 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. / In the recent past years we have seen through various published papers an increasing interest in true amplitude migration methods, in order to obtain more informations about the reflectivity properties of the earth subsurface. The most part of these works has treated of this thema either based on Born approximation as given by Beistein (1987) and Bleistein et al. (1987), or on ray theoretical wavefield approximation as given by Hubral et al. (1991), Schleicher et al. (1993) and Martins et al. (1997). By considering arbitrary source-receiver configurations the compressional primary reflections can be imaged into time or depth-migrated reflections so that the migrated wavefield amplitudes are a measure of angle-dependent reflection coeffients. In order to do this various migration algorithms were proposed in the recent past years based on Born or Kirchhoff approach. Both of them treats of a weighted diffraction stack integral operator that is applied to the input seismic data. As result we have a migrated seismic section where at each reflector point there is the source wavelet with the amplitude proportinal to the reflection coefficient at that point. Based on Kirchhoff approach, in this thesis we derive the weight function and the diffraction stack integral operator for the two dimensional (2-D) and for the two and one half (2.5-D) seismic model and apply it to a set of synthetic seismic data in noise environment. The result shows the accuracy and stability of the 2-D and 2.5-D migration methods as a tool for obtaining important information about the reflectivity properties of the earth subsurface, which is of great interest for the amplitude versus offset (angle) analysis. In summary, we present an expressions for the 2-D and 2.5-D weights as a function of parameters along each ray branch of the in-plane trajectory. Moreover, we show examples of application of the true-amplitude depth migration algorithm to synthetic seismic data obtained by ray theory seismic modeling using the Seis88 package (Cervený e Psencík, 1988), in order to make a numerical analysis and to verify the stability and accuracy of the algorithm. The results confirmed the removal of the geometrical spreading from migrated data, even in presence of noise. Additional tests were performed for pulse distortion analysis in depth rnigrated sections (Tygel et al., 1994) and to obtain reflection points attributes by multiple diffraction stack (Tygel et al., 1993).