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).

Made available in DSpace on 2011-03-23T21:19:32Z (GMT). No. of bitstreams: 0 / Item created via OAI harvest from source: http://www.bdtd.ufpa.br/tde_oai/oai2.php on 2011-03-23T21:19:32Z (GMT). Item's OAI Record identifier: oai:bdtd.ufpa.br:233 / ANP - Agência Nacional do Petróleo / The Common-Reflection-Surface (CRS) stack is a new seismic processing method for simulating zero-offset (ZO) and common-offset (CO) sections. It is based on a second-order hyperbolic paraxial approximation of reflection traveltimes in the vicinity of a central ray. For ZO section simulation the central ray is a normal ray, while for CO section simulation the central ray is a finite-offset ray. In addition to the ZO section, the CRS stack method also provides estimates of wavefield kinematic attributes useful for solving interval velocity inversion, geometrical spreading calculation, Fresnel zone estimate, and also diffraction events simulation. In this work, Its proposed a new strategy to do a pre-stack depth migration by using the CRS derived wavefield kinematic attributes, so-called CRS based pre-stack depth migration (CRS-PSDM) method. The CRS-PSDM method uses the CRS results (ZO section and kinematic attributes) to construct an optimized stack traveltime surface along which the amplitudes of the multi-coverage seismic data are to be summed and the result is put in a point of the migration target zone in depth. In the same sense as in Kirchhoff type pre-stack depth migration (K-PSDM), the CRSPSDM method needs a migration velocity model. Unlike the K-PSDM method, the CRS-PSDM needs only to calculate the zero-offset traveltimes, i.e, along only ray conecting the considered point in depth to a given coincident position of source-receptor at surface. The final result is a zero-offset time-to-depth converted seismic image of reflectors from pre-stack seismic data. / O empilhamento por superfície de reflexão comum (ou empilhamento SRC), conhecido como empilhamento CRS, do inglês Commom reflection surface, constitui-se em um novo método para o processamento sísmico na simulação de seções afastamento nulo (AN) e afastamento comum (AC). Este método é baseado em uma aproximação paraxial hiperbólica de segunda ordem dos tempos de trânsito de reflexão na vizinhança de um raio central. Para a simulação de seção AN, o raio central é um raio normal, enquanto que para a simulação de uma seção AC o raio central é um raio de afastamento finito. Em adição à seção AN, o método de empilhamento SRC também fornece estimativas dos atributos cinemáticos do campo de onda, sendo aplicados, por exemplo, na determinação (por um processo de inversão) da velocidade intervalar, no cálculo do espalhamento geométrico, na estimativa da zona de Fresnel, e também na simulação de eventos de tempos de difrações, este último tendo uma grande importância para a migração pré-empilhamento. Neste trabalho é proposta uma nova estratégia para fazer uma migração em profundidade pré-empilhamento, que usa os atributos cinemáticos do campo de onda derivados do empilhamento SRC, conhecido por método CRS-PSDM, do inglês CRS based pre-stack depth migration. O método CRS-PSDM usa os resultados obtidos do método SRC, isto é, as seções dos atributos cinemáticos do campo de onda, para construir uma superfície de tempos de trânsito de empilhamento, ao longo da qual as amplitudes do dado sísmico de múltipla cobertura são somadas, sendo o resultado da soma atribuído a um dado ponto em profundidade, na zona alvo de migração que é definida por uma malha regular. Similarmente ao método convencional de migração tipo Kirchhoff (K-PSDM), o método CRS-PSDM precisa de um modelo de velocidade de migração. Contrário ao método K-PSDM, o método CRS-PSDM necessita apenas computar os tempos de trânsito afastamento nulo, ao seja, ao longo de um único raio ligando o ponto considerado em profundidade a uma dada posição de fonte e receptor coincidentes na superfície. O resultado final deste procedimento é uma imagem sísmica em profundidade dos refletores a partir do dado de múltipla cobertura.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-20T12:07:56Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AtenuacaoMultiplasCompressao.pdf: 24452115 bytes, checksum: 8631331284b78acc0b168236fbe6c6ea (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Ausência de palavras-chave on 2014-08-07T12:42:09Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-10T15:30:42Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AtenuacaoMultiplasCompressao.pdf: 24452115 bytes, checksum: 8631331284b78acc0b168236fbe6c6ea (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-18T11:32:11Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AtenuacaoMultiplasCompressao.pdf: 24452115 bytes, checksum: 8631331284b78acc0b168236fbe6c6ea (MD5) / Made available in DSpace on 2014-09-18T11:32:11Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AtenuacaoMultiplasCompressao.pdf: 24452115 bytes, checksum: 8631331284b78acc0b168236fbe6c6ea (MD5) Previous issue date: 2003 / O objetivo central deste trabalho é o estudo e a aplicação do método Kalman-Bucy no processo de deconvolução ao impulso e de deconvolução com predição, onde é considerado que os dados observados são classificados como não-estacionários. Os dados utilizados neste trabalho são sintéticos e, com isto, esta Tese tem características de um exercício numérico e investigativo. O operador de deconvolução ao impulso é obtido a partir da teoria de CRUMP (1974) fazendo uso das soluções das equações Wiener-Hopf apresentadas por KALMAN-BUCY (1961) nas formas contínuas e discretas considerando o processo como não estacionário. O operador de predição (KBCP) está baseado nas teorias de CRUMP (1974) e MENDEL ET AL (1979). Sua estrutura assemelha-se ao filtro Wiener-Hopf onde os coeficientes do operador (WHLP) são obtidos através da autocorrelação, e no caso (KBCP) são obtidos a partir da função bi(k). o problema é definido em duas etapas: a primeira consta da geração do sinal, e a segunda da sua avaliação. A deconvolução realizada aqui é classificada como estatística, e é um modelo fortemente baseado nas propriedades do sinal registrado e de sua representação. Os métodos foram aplicados apenas em dados sintéticos de seção fonte-comum obtida a partir dos modelos com interfaces contínuas e camadas homogêneas. / The main objective of this work is the study and the application of the Kalman-Bucy method in the processo f deconvolution to the impulse and deconvolution with prediction, considering the observed data as no stationary. The data used in this work are synthetic and, with this, this Thesis has characteristics of a numerical and search. The operator of deconvolution to the impulse is obtained from the Crump theory (1974), doing use of the solution of equation of Wiener-Holp presented by Kalman-Bucy in the continuoun and discrete forms considering the stacionary process. The prediction operator (KBCP) is based the Crump (1974) and Mendel et al (1979) theorics. Its structure resembles the Wiener-Hopf filter, where the coefficients of the operator are obtained through the autocorrelation, in the case (KBCP) are obtained from the function bi(k). The problem is defined in two steps: the first consists of the generation of the signal, and second of its evaluation. The deconvolution performed is classified as statistics, and is a model based in the properties of the registered signal and its representation. The method were applied only in synthetic data with common-shot section obtained from models with continuous interfaces and homogeneous layers.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-11T12:53:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoGeofisicaGeologica.pdf: 9834744 bytes, checksum: a2405c7a3ab5bd40721dc8d10090d688 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-07-03T14:41:41Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoGeofisicaGeologica.pdf: 9834744 bytes, checksum: a2405c7a3ab5bd40721dc8d10090d688 (MD5) / Made available in DSpace on 2014-07-03T14:41:41Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoGeofisicaGeologica.pdf: 9834744 bytes, checksum: a2405c7a3ab5bd40721dc8d10090d688 (MD5) Previous issue date: 2013 / UFPA - Universidade Federal do Pará / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / INCT - Instituto Nacional de Ciência e Tecnologia / PETROBRAS - Petróleo Brasileiro S.A. / A utilização dos métodos de reflexão sísmica na exploração e desenvolvimento de reservatórios de hidrocarbonetos ocorre devido à sua vasta e densa amostragem, tanto em área quanto em profundidade, aliada ao refinamento de técnicas para o tratamento dos dados de reflexão sísmica, a partir destes dados, são geradas seções sísmicas, que após a aplicação de tratamento adequado, são utilizadas na interpretação dos estratos e/ou estruturas geológicas da subsuperfície. Neste trabalho é feita uma análise Geofísica Geológica de duas linhas sísmicas reais 2D marinhas da porção de quebra de talude da Bacia do Jequitinhonha. Para tanto, foi realizado um conjunto de processamento sísmico com objetivo de atenuar as reflexões múltiplas comuns em dados marinhos, além disso, foram estimados os modelos de velocidade em profundidade, utilizados para determinação das seções sísmicas migradas em profundidade. Nestas foram identificadas as superfícies refletoras. Através da análise dessas superfícies foram feitas as marcações de sismofácies, com base nos conceitos iniciais da sismoestratigrafia, com a finalidade de avaliar a qualidade do produto derivado do processamento sísmico, empregado neste estudo, para uma interpretação sismoestratigráfica, a qual está fundamentada na análise dos padrões de terminações dos refletores e padrão interno das sismofácies. / The seismic reflection methods use at the exploration and development hydrocarbon reservoirs is due to its vast and dense sampling in both area and depth, combined with techniques refinement to processing seismic reflection data, from these data, are generated seismic sections, which after applying appropriate treatment, are used in the interpretation of the strata and/or the subsurface geologic structures. At this paper have made one Geophysics Geological analysis at two real marine seismic lines 2D of break portion embankment Jequitinhonha Basin. Therefore, we performed a set of seismic processing in order to mitigate the multiple reflections common marine data, moreover, were estimated velocity models in depth, used to determine the depth migrated seismic sections. At these were identified reflective surfaces. Through analysis at these surfaces markings of sismofacies. have made based on the initial concepts seismic stratigraphy in order to evaluate the product quality derived from the seismic processing used at this study for seismic stratigraphic interpretation, which is based on the analysis termination patterns of reflectors and internal standard of sismofacies.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-20T14:53:12Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoMetodoSuperficie.pdf: 6623577 bytes, checksum: 67b656d135334821f2f993c805d88214 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Ausência de palavras-chave on 2014-08-07T12:41:40Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-13T14:27:13Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoMetodoSuperficie.pdf: 6623577 bytes, checksum: 67b656d135334821f2f993c805d88214 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-20T14:39:54Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoMetodoSuperficie.pdf: 6623577 bytes, checksum: 67b656d135334821f2f993c805d88214 (MD5) / Made available in DSpace on 2014-08-20T14:39:54Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoMetodoSuperficie.pdf: 6623577 bytes, checksum: 67b656d135334821f2f993c805d88214 (MD5) Previous issue date: 2003 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / CTPETRO - Plano Nacional de Ciência e Tecnologia do Setor Petróleo e Gás Natural / FINEP - Financiadora de Estudos e Projetos / O método de empilhamento sísmico CRS simula seções sísmicas ZO a partir de dados de cobertura múltipla, independente do macro-modelo de velocidades. Para meios 2-D, a função tempo de trânsito de empilhamento depende de três parâmetros, a saber: do ângulo de emergência do raio de reflexão normal (em relação à normal da superfície) e das curvaturas das frentes de onda relacionadas às ondas hipotéticas, denominadas NIP e Normal. O empilhamento CRS consiste na soma das amplitudes dos traços sísmicos em dados de múltipla cobertura, ao longo da superfície definida pela função tempo de trânsito do empilhamento CRS, que melhor se ajusta aos dados. O resultado do empilhamento CRS é assinalado a pontos de uma malha pré-definida na seção ZO. Como resultado tem-se a simulação de uma seção sísmica ZO. Isto significa que para cada ponto da seção ZO deve-se estimar o trio de parâmetros ótimos que produz a máxima coerência entre os eventos de reflexão sísmica. Nesta Tese apresenta-se fórmulas para o método CRS 2-D e para a velocidade NMO, que consideram a topografia da superfície de medição. O algoritmo é baseado na estratégia de otimização dos parâmetros de fórmula CRS através de um processo em três etapas: 1) Busca dos parâmetros, o ângulo de emergência e a curvatura da onda NIP, aplicando uma otimização global, 2) busca de um parâmetro, a curvatura da onda N, aplicando uma otimização global, e 3) busca de três parâmetros aplicando uma otimização local para refinar os parâmetros estimados nas etapas anteriores. Na primeira e segunda etapas é usado o algoritmo Simulated Annealing (SA) e na terceira etapa é usado o algoritmo Variable Metric (VM). Para o caso de uma superfície de medição com variações topográficas suaves, foi considerada a curvatura desta superfície no algoritmo do método de empilhamento CRS 2-D, com aplicação a dados sintéticos. O resultado foi uma seção ZO simulada, de alta qualidade ao ser comparada com a seção ZO obtida por modelamento direto, com uma alta razão sinal-ruído, além da estimativa do trio de parâmetros da função tempo de trânsito. Foi realizada uma nálise de sensibilidade para a nova função de tempo de trânsito CRS em relação à curvatura da superfície de medição. Os resultados demonstraram que a função tempo de trânsito CRS é mais sensível nos pontos-médios afastados do ponto central e para grandes afastamentos. As expressões da velocidade NMO apresentadas foram aplicadas para estimar as velocidades e as profundidades dos refletores para um modelo 2-D com topografia suave. Para a inversão destas velocidades e profundidades dos refletores, foi considerado o algoritmo de inversão tipo Dix. A velocidade NMO para uma superfície de medição curva, permite estimar muito melhor estas velocidades e profundidades dos refletores, que as velocidades NMO referidas as superfícies planas. Também apresenta-se uma abordagem do empilhamento CRS no caso 3-D. neste caso a função tempo de trânsito depende de oito parâmetros. São abordadas cinco estratégias de busca destes parâmetros. A combinação de duas destas estratégias (estratégias das três aproximações dos tempos de trânsito e a estratégia das configurações e curvaturas arbitrárias) foi aplicada exitosamente no empilhamento CRS 3-D de dados sintéticos e reais. / The CRS stacking method simulates ZO seismic sections from multi-coverage data and does not dependente on a macro-velocity model. For 2-D medium the stacking traveltime depends on three parameters: the emergence angle of the normal ray (with respect to the measurement surface normal) and the wavefront curvatures of two hypothetical waves, called Normal-Incidence-Point (NIP) wave and Normal (N) wave. The CRS method consists of summing the amplitudes of the seismic traces in the multicoverage data along the surface defined by CRS stacking traveltime which that fits best the data set. The result of the CRS stack is assigned to points of a grid pre-defined in the ZO section. As the result obtain a simulated ZO section. This means that for each point of the ZO section must be estimated the three optimal parameters that yield the maximum coherence between the events of seismic reflection. In this Thesis I present formulae for the 2-D CRS method and for the NMO velocity that consider the topography of the measurement surface. The algorithm is based on the optimization strategy divided into three steps: 1) To search for the emergence angle and the curvature of the NIP wave, by applying a global optimization, 2) to search for the curvature of the N wave, by applying global optimization, and 3) to refine the initial parameters estimated in first two steps by applying local optimization. In the first two steps is used the Simulated Annealing (SA) algorithm and in the third step the Variable Metric (VM) algorithm. For the case of a measurement surface with smooth topography the curvature of this surface is included in the 2-D CRS stack formalism. This CRS algorithm implemented was applied to synthetic data set. The result is a simulated ZO section of high quality, with a high signal-to-noise ratio, and the estimative of the parameter triplet. It is performed a sensibility analysis for the new CRS stacking traveltime with respect to the curvature in several points of the curved measurement surface. This study showed that the CRS traveltime is more sensitive for fast midpoints of the central points and larger offsets. The expressions for the NMO velocities presented here is applied to estimate the interval velocities and the depth of the reflectors for 2-D model with a smooth topography. For the inversion of the velocities and the depth of the reflectors is considered the Dix-type inversion algorithm. The NMO velocity for a curved measurement surface deserves to best estimate the velocities and the depths of the reflectors than NMO velocities referred to planar surfaces. Also, I present an introduction to 3-D stack. In this case, the stacking traveltime depends on eight parameters. These parameters can be obtained by using some parameter-search strategies that I have showed in this Thesis. The combination of the strategy of the Traveltime Approximations and the strategy of Arbitrary Curvatures is used to apply 3-D CRS stack successful in synthetic and real data sets, respectively.