Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-10T13:35:01Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AproximacoesHiperbolicasTemposTransito.pdf: 2945184 bytes, checksum: cf9d6cdca9479875cda71b0d7c22cfda (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar as palavras-chave on 2014-08-07T13:37:42Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-10T12:08:20Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AproximacoesHiperbolicasTemposTransito.pdf: 2945184 bytes, checksum: cf9d6cdca9479875cda71b0d7c22cfda (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T12:54:16Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AproximacoesHiperbolicasTemposTransito.pdf: 2945184 bytes, checksum: cf9d6cdca9479875cda71b0d7c22cfda (MD5) / Made available in DSpace on 2014-09-10T12:54:16Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AproximacoesHiperbolicasTemposTransito.pdf: 2945184 bytes, checksum: cf9d6cdca9479875cda71b0d7c22cfda (MD5) Previous issue date: 2005 / A simulação de uma seção sísmica de afastamento-nulo (AN) a partir de dados de cobertura múltipla pode ser realizada através do empilhamento sísmico, o qual é um método de imageamento de reflexão sísmica muito utilizado na indústria. O processo de empilhamento sísmico permite reduzir a quantidade de dados e melhorar a razão sinal/ruído. Baseado na aproximação hiperbólica dos tempos de trânsito dependente de três parâmetros ou atributos cinemáticos de frentes de onda. Recentemente, vem desenvolvendo-se um novo método para simular seções (AN) chamado método de empilhamento sísmico por Superfície de Reflexão Comum (ou empilhamento SRC). Este novo formalismo pode ser estendido para construir seções de afastamento-nulo (AN) a partir de dados de cobertura múltipla, usando aproximações dos tempos de trânsito paraxiais na vizinhança de um raio central com afastamento-nulo (AN), para o caso de uma linha de medição com topografia suave e rugosa. Essas duas aproximações de tempos de trânsito também dependem de três atributos cinemáticos de frentes de ondas. Nesta dissertação, apresenta-se uma revisão teórica da teoria paraxial do raio para a obtenção das aproximações dos tempos de trânsito paraxiais considerando uma linha de medição com topografia rugosa e suave. A partir das aproximações dos tempos de trânsito paraxiais em relação a um raio central com afastamento-nulo (AN), foram obtidas duas novas aproximações de tempos de trânsito usando a condição de um ponto difrator em profundidade, reduzindo as equações originais para dois parâmetros. Também foram obtidas as aproximações para o caso de raios paraxiais com afastamento-nulo (AN). Para as aproximações de reflexão e difração utilizando um mesmo modelo sintético, foram comparadas através da representação gráfica as superfícies de empilhamento das aproximações dos tempos de trânsito para a topografia suave e rugosa. Em seguida, para analisar o comportamento dos operadores associados a reflexão e a difração, quando estes são perturbados, discutimos suas sensibilidades em relação a cada um dos três parâmetros (β₀, K_PIN, K_N). Esta análise de sensibilidade foi realizada em duas formas: Sensibilidade através da perturbação de cada parâmetro visualizado nas superfícies de empilhamento SRC-TR e SDC-TR e da primeira derivada dos tempos de trânsito SRC-TR e SDC-TR. Finalmente, usando essas aproximações hiperbólicas em função de três e dois parâmetros e com base nos resultados da análise de sensibilidade, foi proposto um algoritmo para simular seções AN a partir de dados de cobertura múltipla. / The simulation of a zero-offset seismic section (ZO) from multicoverage data can be carried through the seismic stacking, which is a very used method of seismic reflection imaging in the oil industry. Seismic stacking process allows to reduce the amount of data and is aimed to improve the signal/noise ratio. Based in the hyperbolic approach of the traveltime, dependent on three parameters or kinematic attributes of wavefronts, recently it was developed a new method to simulate zero-offset sections (ZO), called seismic stack for Common Reflection Surface method (or SRC stack). This new formalism can be extended to construct ZO sections from multicoverage data, using the approach of the paraxial traveltimes in the neighborhood of a central ray with zero offset, for the case of a line of measurement with soft and rugged topography. These two approaches of traveltime also depend on three kinematic attributes of wavefronts. In this work, a theoretical revision of the paraxial theory of the ray for the attainment of the approaches of the paraxial traveltimes is presented, considering a line of measurement with rugged and soft topography. From the approaches of the paraxial traveltimes with relation to a central ray with zero-offset (ZO), two new approaches had been derived from traveltimes using the condition of a diffraction point in depth, reducing the original equations for two parameters. For the approaches of reflection and diffraction, using the same synthetic model, a graphical representation of their respective stacking surfaces were compared, using the traveltime approaches for the soft and rugged topography. After that, we analyze the behavior of the operators associated with the reflection and the diffraction surfaces when these are disturbed, where we investigate the sensibilities in relation to each one of the three parameters (β₀, K_NIP, K_M). This sensibility analysis was perfomed in two ways: Sensibility through the disturbance of each parameter visualized in the stacking surfaces SRC-TR and SDC-TR, and the first derivative of the traveltimes SRC-TR and SDC-TR. Finally, using these hyperbolic approaches as a function of three and two parameters, and based on the results of the sensibility analysis, were considered an algorithm to simulate ZO sections from data of multi coverage.

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

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

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-12T13:56:10Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoFiltrosVelocidade.pdf: 4324452 bytes, checksum: 90dbfc8b340951c647cbf95c43f4bab3 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar as palavras-chave on 2014-08-07T13:38:08Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-10T13:51:07Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoFiltrosVelocidade.pdf: 4324452 bytes, checksum: 90dbfc8b340951c647cbf95c43f4bab3 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T14:49:50Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoFiltrosVelocidade.pdf: 4324452 bytes, checksum: 90dbfc8b340951c647cbf95c43f4bab3 (MD5) / Made available in DSpace on 2014-09-10T14:49:50Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoFiltrosVelocidade.pdf: 4324452 bytes, checksum: 90dbfc8b340951c647cbf95c43f4bab3 (MD5) Previous issue date: 2004 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / A motivação geológica deste trabalho reside no imageamento de estruturas de bacias sedimentares da região Amazônica, onde a geração e o acúmulo de hidrocarboneto estão relacionados com a presença de soleiras de diabásio. A motivação sísmica reside no fato de que essas rochas intrusivas possuem grandes contrastes de impedância com a rocha encaixante, o que resulta em múltiplas, externas e internas, com amplitudes semelhantes as das primárias. O sinal sísmico das múltiplas podem predominar sobre o sinal das reflexões primárias oriundas de interfaces mais profundas, o que pode dificultar o processamento, a interpretação e o imageamento da seção sísmica temporal. Neste trabalho, estudamos a atenuação de múltiplas em seções sintéticas fonte-comum (FC) através da comparação de dois métodos. O primeiro método resulta da combinação das técnicas Wiener-Hopf-Levinson de predição (WHLP) e o de empilhamento superfície-de-reflexão-comum (CRS), e denominando WHLP-CRS, onde o operador é desenhado exclusivamente no domínio do tempo-espaço. O segundo método utilizado é o filtro de velocidade (ω-k) aplicado após o empilhamento superfície-de-reflexão (CRS), onde o operador é desenhado exclusivamente no domínio bidimensional de freqüência temporal-espacial. A identificação das múltiplas é feita na seção de afastamento-nulo (AN) simulada com o empilhamento CRS, e utiliza o critério da periodicidade entre primária e suas múltiplas. Os atributos da frente de onda, obtidos através do empilhamento CRS, são utilizados na definição de janelas móveis no domínio tempo-espaço, que são usadas para calcular o operador WHLP-CRS. O cálculo do filtroω-k é realizado no domínio da freqüência temporal-espacial, onde os eventos são selecionados para corte ou passagem. O filtro (ω-k) é classificado como filtro de corte, com alteração de amplitude, mas não de fase, e limites práticos são impostos pela amostragem tempo-espaço. Em termos práticos, concluímos que, para o caso de múltiplas, os eventos separados no domínio x-t não necessariamente se separam no domínio ω-k, o que dificulta o desenho de um operador ω-k semelhante em performance ao operador x-t. / The geological motivation of this work is the imaging of sedimentary basin structures of the Amazon region, where the generation and accumulation of hydrocarbons is related to the existence of diabase sills. The seismic motivation is the fact that these intrusive rocks present a great impedance contrast with respect to the host rock, what gives rise to external and internal multiples, with primary-like amplitudes. The seismic signal of the multiples can predominate over the primary reflection signals from deeper interfaces, making difficult the processing, interpretation and imaging of seismic sections. In this work we study the attenuation of multiples in common shot (CS) sections by the comparison of two methods. The first one is the combination of the Wiener-Hopf-Levinson (WHLP) and the common reflection surface (CRS) stacking techniques, here called WHLPCRS, where the operator is exclusively designed in the space-time domain. The second method is a velocity filter (ω-k), applied after the CRS stacking, where the operator is exclusively designed in the frequency-wavenumber domain. The identification of the multiples is performed on the zero-offset (ZO) section simulated by the CRS stacking, using the periodicity between primaries and its multiples. The wavefront attributes, obtained through CRS stacking, are inserted on movable space-time windows, used to calculate the WHLP-CRS operator. The ω-k filter calculations are performed in the frequency-wavenumber domain, where the events of interest are selected for cutting or passage. The ω-k filter is classified as a cut-off filter, with amplitude alteration and preservation of phase, the limits of it are imposed by a space-time sampling. In practical aspects we conclude that for the case of multiples, separated events on the x-t domain are not necessarily separated on the ω-k domain, which raise difficulties in the designing of a ω-k operator with a similar performance when compared to the x-t operator.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-10T13:10:33Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AproximacoesHiperbolicasTempos.pdf: 4094530 bytes, checksum: 2d1b416e88ea1617317477592ae645d5 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-06T12:35:22Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-10T12:19:23Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AproximacoesHiperbolicasTempos.pdf: 4094530 bytes, checksum: 2d1b416e88ea1617317477592ae645d5 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-18T11:33:11Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AproximacoesHiperbolicasTempos.pdf: 4094530 bytes, checksum: 2d1b416e88ea1617317477592ae645d5 (MD5) / Made available in DSpace on 2014-09-18T11:33:11Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AproximacoesHiperbolicasTempos.pdf: 4094530 bytes, checksum: 2d1b416e88ea1617317477592ae645d5 (MD5) Previous issue date: 2005 / A simulação de uma seção sísmica de afastamento nulo (AN) 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 melhor a razão sinal/ruído. Baseado na aproximação hiperbólica dos tempos de trânsito dependente de três parâmetros ou atributos cinemáticos de frentes de onda, recentemente, vem desenvolvendo-se um novo método para simular seções (AN) chamado método de empilhamento sísmico por Superfície de Reflexão Comum (ou empilhamento SRC). Também, seguindo este novo conceito de imageamento sísmico está surgindo um método para simular seções com afastamento comum (AC) a partir de dados de cobertura múltipla usando aproximações dos tempos de trânsito paraxiais na vizinhança de um raio central com afastamento finito. Esta nova aproximação dos tempos de trânsito depende de cinco atributos cinemáticos. Neste trabalho, a partir da aproximação dos tempos de trânsito paraxiais em relação a um raio central com afastamento finito, derivamos uma nova equação do tempo de trânsito usando a condição de um ponto difrator em profundidade, reduzindo a equação original para quatro parâmetros. Para ambas aproximações (reflexão e difração), mostramos a superfície de empilhamento SRC com afastamento finito (SRC-AF). Considerando um modelo sintético, realizamos um estudo comparativo das aproximações dos tempos de trânsito para as quatro configurações sísmicas (fonte comum (FC), receptor comum (RC), ponto-médio-comum (PMC) e afastamento comum (AC)). Para analisar o comportamento do operador SRC-AF, quando este é perturbado, discutimos sua sensibilidade em relação a cada um dos cinco parâmetros (K₁, K₂, K₃, β_S e β_G). Esta análise de sensibilidade é realizada em duas formas: Sensibilidade através da primeira derivada e Sensibilidade no Empilhamento SRC-AF. Após realizar a análise de sensibilidade utilizamos uma nova condição, K₂ = 0 e assim, obtemos uma nova aproximação, agora dependente de três parâmetros. Usando essas aproximações hiperbólicas (em função de cinco, quatro e três parâmetros), propomos um algoritmo para a simulação de seções AC a partir de dados de cobertura múltipla. Finalmente, é apresentado um estudo da zona de Fresnel, com o objetivo de determinar a delimitação da abertura da superfície de empilhamento SRC-AF. / The simulation of a zero-offset seismic (ZO) section from multicoverage, through of stack, is a very used method of seismic reflection imaging, that reduce the amount of data and is aimed to improve the signal/noise ratio. Based on hyperbolic traveltime approximation depending on three kinematic attributes, recently, it was developed a new method to simulate zero-offset (ZO) sections called common reflection surface stack (CRS) method. Also, following this new concept of seismic imaging it was introduzed a method to simulate common-offset (CO) sections from multicoverage data by using a paraxial traveltime approximation for paraxial rays in the vicinity of the finite-offset centra ray. This new traveltime approximation depend on five kinematic attributes. In this work, from the hyperbolic paraxial traveltime approximation with finite-offset central ray, we derive a new formula of traveltime approximation by using the diffraction point condition, reducing the original formula to four parameters. For both approximations, that´s for reflection and diffractions, we show the finite-offset FO-CRS stack operators. We, also obtain these formulas to the four seismic configurations (common-shot (CS), common-receiver (CR), common-midpoint (CMP) and common-offset (CO)). To analise behaviour FO-CRS stack, when this is perturbed, we debatable the sensibility with respect to each one the five parameters (K₁, K₂, K₃, β_S and β_G). This sensibility analysis was perfomed in two ways: Sensibility through the disturbance of each parameter visualized in the stacking surfaces FO-CRS- and FO-CDS, and the first derivative of the traveltimes FO-CRS and FO-CDS. After realize the sensibility analysis, we utility a new condition, K₂ = 0 and therefore, we obtained a new approximation, now depend on three parameters. Using that hyperbolic traveltime approximations (in function of five, four and three parameters), we propose an algorithm multicoverage seismic reflection data. Finally, is presented a studed of Fresnel zone, to delimit the aperture of FO-CRS stacking surface.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-07-09T16:19:06Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoMetodoInversao.pdf: 6449611 bytes, checksum: fb08c6d600df18e726054ead8c8a3047 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:05:28Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-11T12:17:31Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoMetodoInversao.pdf: 6449611 bytes, checksum: fb08c6d600df18e726054ead8c8a3047 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-15T16:01:19Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoMetodoInversao.pdf: 6449611 bytes, checksum: fb08c6d600df18e726054ead8c8a3047 (MD5) / Made available in DSpace on 2014-09-15T16:01:19Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoMetodoInversao.pdf: 6449611 bytes, checksum: fb08c6d600df18e726054ead8c8a3047 (MD5) Previous issue date: 2007 / Neste trabalho, estuda-se um novo método de inversão tomográfica de reflexão para a determinação de um modelo isotrópico e suave de velocidade por meio da aplicação, em dados sintéticos e reais, do programa Niptomo que é uma implementação do método de inversão tomográfica dos atributos cinemáticos da onda hipotética do ponto de incidência normal (PIN). Os dados de entrada para a inversão tomográfica, isto é, o tempo de trânsito e os atributos da onda PIN (raio de curvatura da frente de onda emergente e ângulo de emergência), são retirados de uma série de pontos escolhidos na seção afastamento nulo (AN) simulada, obtida pelo método de empilhamento por superfícies de reflexão comum (SRC). Normalmente, a escolha destes pontos na seção AN é realizada utilizando-se programas de picking automático, que identificam eventos localmente coerentes na seção sísmica com base nos parâmetros fornecidos pelo usuário. O picking é um dos processos mais críticos dos métodos de inversão tomográfica, pois a inclusão de dados de eventos que não sejam de reflexões primárias podem ser incluídos neste processo, prejudicando assim o modelo de velocidades a ser obtido pela inversão tomográfica. Este trabalho tem por objetivo de construir um programa de picking interativo para fornecer ao usuário o controle da escolha dos pontos de reflexões sísmicas primárias, cujos dados serão utilizados na inversão tomográfica. Os processos de picking e inversão tomográfica são aplicados nos dados sintéticos Marmousi e nos dados da linha sísmica 50-RL-90 da Bacia do Tacutu. Os resultados obtidos mostraram que o picking interativo para a escolha de pontos sobre eventos de reflexões primárias favorece na obtenção de um modelo de velocidade mais preciso. / In this work, a new tomographic reflection inversion method is studied, with the objective of determining smooth and isotropic velocity models for real and synthetic seismic data. The velocity model is obtained by the program NIPTOMO, which is an implementation of a tomographic inversion method of the kinematic attributes of the hypothetical waves NIP (normal incidence point). The input data for the tomographic inversion, e. g., traveltime tables and NIP attributes (wavefront curvature radius and emergence angle) are obtained through a set of picked points on a simulated zero-offset (ZO) section, previously obtained by the common-reflection surface (CRS) method. The choice of picked points in a ZO section is normally made using automatic picking programs that identify locally coherent events in a seismic section, based in user-defined parameters. The picking process is one of the most critical steps in this new method, since the choice of non primary reflection events, such as diffractions and multiples, may lead to errors in the velocity model obtained by the tomographic inversion. This work has as main objective develop an interative picking program to give to the user the control of the choice of primary reflection points that are used as input for the tomographic inversion. The present picking program was tested on the synthetic Marmousi dataset and on a real data from the Tacutu Basin, Brazil (seismic line 50-RL-90). The results obtained show that the interative picking process for the choice of primary reflection events points favors the determination of a precise velocity model.