Orientador: Martin Tygel / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica e Instituto de Geociências / Made available in DSpace on 2018-08-23T04:02:48Z (GMT). No. of bitstreams: 1 RuedaSerrano_Dany_M.pdf: 36981147 bytes, checksum: 4818cdd1a4fd5bda74072b4e6d58a7ba (MD5) Previous issue date: 2013 / Resumo: A razão sinal-ruído, geralmente abreviada por razão S/N (do Inglês signal-to-noise ratio) é bastante afetada por ruídos aleatórios, os quais degradam a continuidade e identificação de refletores, com prejuízos para interpretação geológica. Com vistas à superação dessas dificuldades, apresentamos neste trabalho uma técnica baseado na suavização de parâmetros obtidos pela aplicação do método de empilhamento Common-Reection-Surface (CRS). A suavização é realizada através da utilização de estatísticas locais em pequenas janelas alinhadas com os eventos de reflexão. Primeiramente, os parâmetros CRS são estimados através de uma aplicação convencional do método CRS. Em seguida, os parâmetros CRS são suavizados de modo a eliminar flutuações e valores espúrios. Finalmente, os dados originais são submetidos a um novo empilhamento CRS, com utilização dos parâmetros suavizados recém obtidos. O processo pode ser aplicado de maneira iterativa para otimização de resultados. O esquema proposto foi aplicado em um dois conjuntos de dados marítimos 2D, um sintético e um de campo. No dado sintético, obtivemos uma efetiva atenuação do ruído aleatório com significativa ênfase dos eventos de reflexão. A aplicação nos dados de campo resultou num significativo increment da razão S/N com consequente ênfase e aumento de continuidade dos refletores. A interpretação sísmica dos elementos estratigráficos tornou-se mais precisa na seção pós-migrada em tempo (PostSTM) correspondente ao método CRS. A seção migrada foi obtida com a utilização dos parâmetros suavizados em conjunto com um guia de velocidade / Abstract: Signal-to-noise (S/N) ratio in seismic data is affected by random noise, influencing the continuity and identification of reflectors. In this work, I present a method to overcome this problem based on smoothing Common-Reflection-Surface (CRS) parameters through application of local statistics in small windows aligned with reflection events. First, the CRS parameters are obtained by a standard application of the CRS stack method. Subsequently, the CRS parameters are smoothed so as to eliminate fluctuations and outliers. Finally, a CRS stack is performed with the new, smoothed parameters. The process may be iterative, to achieve optimal results. The proposed scheme has been applied on a 2D synthetic data set and marine field data. The synthetic data application showed effective random noise attenuation plus highlighting of the reflection events. Application to the real marine data resulted in an increase of S/N ratio with consequent highlight and greater continuity of the reflections. The seismic interpretation of stratigraphic elements has been more precise in the PostSTM section related with the CRS stack method. It is performed using smoothed attributes and a velocity model as a guide / Mestrado / Reservatórios e Gestão / Mestre em Ciências e Engenharia de Petróleo

Orientadores: Martin Tygel, Eduardo Filpo Ferreira da Silva / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-25T22:08:42Z (GMT). No. of bitstreams: 1 Cardoso_ClaudioGuerra_M.pdf: 5757657 bytes, checksum: 5b4ea3b5608c625fce78a747dffdf8e0 (MD5) Previous issue date: 1999 / Resumo: Atualmente, cerca de 90% das reservas brasileiras de óleo localizam-se nas porções submersas das bacias sedimentares costeiras. Nessas bacias, a sísmica de reflexão é uma das principais ferramentas de que dispõe o profissional de exploração e de desenvolvimento de campos de petróleo. Como a interpretação dos dados sísmicos considera somente as reflexões primárias, é necessário que estes estejam livres de ruídos, particularmente, da reflexão múltipla do fundo do mar, que é o ruído mais danoso em dados marítimos. A presente dissertação descreve um método novo de eliminação de múltipla do fundo do mar de 1ª ordem em configuração de afastamento fonte - receptor comum, por empilhamento Kirchhoff, mostrando exemplos numéricos de aplicação do mesmo em dados sintéticos. Num desses exemplos, no qual concorrem a múltipla e uma reflexão primária - que poderia advir de um refletor de interesse - após a eliminação da múltipla, as diferenças nas amplitudes da primária, quando comparadas com as amplitudes sem a interferência da múltipla, não ultrapassam 5%. Dessa forma, o método, por preservar as amplitudes das primárias após a eliminação da múltipla, se configura como muito adequado à geofísica de reservatórios / Abstract: At present, about 90% of Brazilian oil reserves are located in offshore sedimentary basins. In these basins, reflection seismic is one of the most important tools in exploration and field development. As seismic data interpretation only deals with primary reflections, is important to have noise-free, amplitude preserved primaries. The most harmful noise in marine data is the first-order, sea-bottom multiple reflection. This dissertation describes a new method for first order, sea-bottom, multiple elimination by means of a Kirchhoff - type stacking, in common offset configuration, showing its application to synthetic seismic data. In one of the examples, in which the multiple and a primary reflection interfere, the differences between amplitudes of primary reflection after application of the method and primary computed with no such interference are smaller than 5%. Due to the primaries amplitude-preserving property, this method provides a very reliable image to reservoir geophysics, such as amplitude versus offset analysis (AVO) and seismic attribute determination / Mestrado / Mestre em Engenharia de Petróleo

Orientador: Martin Tygel / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica e Instituto de Geociencias / Made available in DSpace on 2018-08-04T20:27:31Z (GMT). No. of bitstreams: 1 Santiago_TaniaMariaGodinho_M.pdf: 6170921 bytes, checksum: fb1caa1410e2409aa2455f02c3c86e58 (MD5) Previous issue date: 2004 / Resumo: Atualmente, a indústria do petróleo, tanto na área de caracterização de reservatórios quanto na de exploração, tem feito cada vez mais uso do estudo da variação da amplitude sísmica com o afastamento fonte-receptor, ou com o ângulo de reflexão, como ferramenta auxiliar na previsão da existência de hidrocarbonetos nos reservatórios. Como o método, rotineiramente utilizado, de transformação de afastamento para ângulo de reflexão, em dados migrados na configuração afastamento comum, envolve erros que podem diminuir a confiabilidade deste estudo, toma-se necessário uma nova maneira para a obtenção de dados que produzam curvas mais acuradas da variação da amplitude com o ângulo de reflexão (AVA). Neste sentido, apresenta-se nesta dissertação a migração Kirchhoff no domínio de ângulo comum 2,5D por empilhamento em tempo e em amplitude verdadeira, e os resultados de sua aplicação em dados sísmicos sintéticos de três modelos geológicos. Discute-se a influência, nos resultados, dos parâmetros de aquisição dos dados sísmicos e a parametrização da migração, especificamente a abertura e intervalo entre os ângulos de mergulho da migração. São feitas comparações das curvas de amplitude versus o ângulo de reflexão oriundas dos conjuntos de dados migrados no domínio de ângulo comum e no domínio de afastamento comum, e mostra-se a maior eficiência da primeira migração em obter curvas de AVA que se aproximam mais das curvas teóricas / Abstract: The analysis of amplitude variation with offset or reflection angle is of great importance in the oil industry to predict hydrocarbon presence in the reservoir. This analysis is applied to both reservoir and exploration areas. The usual transformation method from offset to reflection angle, for common-offset migrated data, involves errors that may produce unreliable results. Thus, it becomes necessary to find altemative procedures to obtain appropriate data to provide more reliable curves of amplitude versus reflection angle (AVA). This work considers 2.5D true amplitude Kirchhoff-time migration in the common-angle domain, and its application to synthetic seismic data. It also examines the influence of seismic acquisition and migration parameters, aperture and migration dip-angle increment, on the migration results. Comparison between AVA curves obtained from common-angle and common-offset migrated data confirms that the former is a more reliable procedure / Mestrado / Mestre em Ciências e Engenharia de Petróleo

Orientadores: Martin Tygel, Lucio Tunes dos Santos / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica e Instituto de Geociencias / Made available in DSpace on 2018-08-12T15:47:44Z (GMT). No. of bitstreams: 1 Nunes_RodrigoSilva_M.pdf: 3222542 bytes, checksum: a51b6fa6537055a2ce60aa8066fc2be2 (MD5) Previous issue date: 2005 / Resumo: O método CRS (do Inglês "Common Reflection Surface"), introduzido por Peter Hubral em 1997, é uma técnica alternativa ao convencional empilhamento CMP. O empilhamento CRS é capaz de fornecer não apenas uma melhoria nas seções simuladas de afastamento nulo (ZO), mas também um trio de parâmetros que extraem mais informações sobre a propagação do meio do que o único parâmetro do clássico método CMP, a velocidade NMO. O objetivo deste trabalho é analisar o potencial do método CRS em produzir modelos de velocidades sísmicas, analogamente para a computação da velocidade NMO no método CMP. Nós descrevemos os métodos CMP e CRS e comparamos suas performances sobre ilustrações sintéticas e exemplos de dados reais. / Abstract: The Common Reflection Surface (CRS) method, introduced by Peter Hubral in 1997, is an alternative technique to the conventional CMP stacking. The CRS stack is able to provide not only improved simulated zero-offset sections, but also a set of parameters that convey more information of the propagating medium than the single NMO-velocity parameter of the classical CMP method. The objective of this work is to analyze the potential of the CRS method to produce seismic velocity models, analogously to the computation of the NMO velocities in the CMP method. We review the CMP and CRS methods, and compare their performances on illustrative synthetic and real data examples. / Mestrado / Reservatórios e Gestão / Mestre em Ciências e Engenharia de Petróleo

Orientador: Jorg Dietrich Wilhelm Schleicher / Tese (doutorado) - Universidade Estadual de Campinas, Instituto de Matematica, Estatistica e Computação Cientifica / Made available in DSpace on 2018-08-14T00:48:52Z (GMT). No. of bitstreams: 1 Carvalho_RafaelAleixode_D.pdf: 9585016 bytes, checksum: 31508e558a4b01742760ced71847289a (MD5) Previous issue date: 2009 / Resumo: Como os alvos de exploração tornaram-se mais profundos, os comprimentos dos cabos têm aumentado em conformidade, fazendo a aproximação hiperbólica convencional produzir tempos de trânsito cada vez mais imprecisos. Em outras palavras, para as modernas geometrias de aquisição para grandes afastamentos, a aproximação hiperbólica já não é suficiente para horizontalizar a família CMP por causa da não homogeneidade ou anisotropia dos meios. Para resolver este problema, muitas fórmulas para o tempo de trânsito foram propostas na literatura que fornecem aproximações de qualidade diferente. Demonstrou-se que para meios com isotropia transversal vertical (meios VTI), apenas dois parâmetros do tempo de trânsito são suficientes para a realização de todo o processamento temporal, sendo a velocidade NMO e um parâmetro de anisotropia. Por isso, nesta tese, nos concentramos, na dedução de aproximações simples para o tempo de trânsito que dependem de um único parâmetro de anisotropia. Começamos por dar uma visão geral de uma coleção de tais aproximações para o tempo de trânsito encontradas na literatura e comparar suas qualidades. Em seguida, deduzimos um conjunto de novas aproximações para o tempo de trânsito que dependem de um parâmetro baseado em aproximações encontradas na literatura. A principal vantagem das nossas aproximações é que algumas delas são expressões analíticas bastante simples que as tornam fáceis de serem utilizadas, ao mesmo tempo que têm a mesma qualidade ou maior que as fórmulas já estabelecidas. Utilizamos estas aproximações para o tempo de trânsito para uma inversão dos parâmetros de anisotropia. Utilizando uma estimativa da velocidade NMO a partir da análise de velocidades hiperbólica, pode-se estimar o parâmetro anisotrópico a partir de uma aproximação para o tempo de trânsito mais geral. Estendemos o procedimento em dois passos utilizando um termo não hiperbólico mais preciso na aproximação para o tempo de trânsito. As aproximações para o tempo de trânsito deduzidas permitem predizer o viés na estimativa da velocidade NMO, proporcionando assim um meio de corrigir, tanto a estimativa a velocidade NMO, quanto o conseqüente valor do parâmetro de anisotropia. Por meio de um exemplo numérico, demonstramos que a estimativa dos parâmetros do tempo de trânsito, usando este processo iterativo, apresenta considerável melhora. Palavras-chave: Aproximação para tempos de trânsito, meios VTI, análise de velocidade e geofísica / Abstract: As exploration targets have become deeper, cable lengths have increased accordingly, making the conventional two term hyperbolic traveltime approximation produce increasingly erroneous traveltimes. In other words, for modern long-offset acquisition geometries, a hyperbolic traveltime approximation is no longer sufficient to flatten the CMP gather because of medium inhomogeneity or anisotropy. To overcome this problem, many traveltime formulas were proposed in the literature that provide approximations of different quality. It has been demonstrated that for transversly isotropic media with a vertical symmetry axis (VTI media), just two traveltime parameters are sufficient to perform all time-related processing, being the NMO velocity and one anisotropy parameter. Therefore, we concentrate in this thesis,on simple traveltime approximations that depend on a single anisotropy parameter. We start by giving an overview of a collection of such traveltime approximations found in the literature and compare their quality. Next, we derive a set of new single-parameter traveltime approximations based on the ones found in the literature. The main advantage of our approximations is that some of them are rather simple analytic expressions that make them easy to use, while achieving the same quality as the better of the established formulas. We then use these traveltime aproximations for an inversion of the anisotropy parameters. Using an estimate of the NMO velocity from a hyperbolic velocity analysis, one can estimate the anisotropic parameter from a more general traveltime approximation. We extend this two-step procedure using a more accurate nonhyperbolicity term in the traveltime approximation. The used traveltime approximations allow to predict the bias in the NMO velocity estimate, thus providing a means of correcting both the estimated NMO velocity and the resulting anisotropy parameter value. By means of a numerical example, we demonstrate that the estimation of the traveltime parameters, using this iterative procedure, is improved considerably. Keywords: Traveltime approximations, VTI media, velocity analysis and geophysics / Doutorado / Geofisica Computacional / Doutor em Matemática Aplicada

Orientador: Armando Zaupa Remacre / Tese (doutorado) - Universidade Estadual de Campinas, Instituto de Geociencias / Made available in DSpace on 2018-08-09T03:37:21Z (GMT). No. of bitstreams: 1 Sancevero_SergioSacani_D.pdf: 13993022 bytes, checksum: beb2507aee5ca130897dad57c706b2a9 (MD5) Previous issue date: 2007 / Resumo: O processo de caracterização de reservatórios é atualmente uma das etapas mais importantes na exploração, desenvolvimento e produção de um campo. Porém, para que esse processo seja realizado da melhor forma possível é preciso se ter o conhecimento de determinados métodos, que integram as diferentes informações disponíveis. Desse modo, o objetivo principal dessa tese é estudar de forma criteriosa e quantitativa o processo de caracterização de reservatórios do ponto de vista dos dados sísmicos, avaliando antigos e novos métodos, e definindo novas metodologias que possam ser aplicadas de maneira decisiva neste processo. Para que esses métodos pudessem ser avaliados de maneira conclusiva foi utilizado nesta tese um modelo sintético que reproduzisse minimamente algumas características cruciais de determinados reservatórios como a complexa distribuição dos corpos de areia e a presença de corpos com espessura subsísmica que levassem ao limite as técnicas de modelagem tradicionais, proporcionando avaliar novos métodos. Assim, para caracterizar essas complexas feições, foram utilizados nesta tese dois meios principais de interpretação, primeiro a inversão sísmica dando um caráter preditivo ao dado sísmico e por fim a análise multiatributos, dando um caráter classificatório. No caso da inversão sísmica foram utilizados três métodos de obtenção da impedância acústica. Entre eles, foi a inversão geoestatística que demonstrou ser a mais eficiente das técnicas no que diz respeito à caracterização de reservatórios com espessura subsísmica e complexa distribuição dos corpos de areia. No caso dos atributos, pôde-se demonstrar que é necessário que sejam tratados com uma abordagem multivariada para que seja aproveitada a correlação entre eles e que por meio de técnicas de classificação e modelagem possa se decidir os mais relevantes para o processo. Neste caso 3 métodos de análise multivariada foram apresentados e testados, sendo que dois deles (ICA e MAF) de maneira inédita e que produziram resultados superiores àqueles alcançados quando a tradicional técnica de PCA é aplicada. Assim, com o que foi apresentado, pode-se concluir que o processo de caracterização é um estágio crucial para o desenvolvimento dos campos, mas não é fácil de ser realizado, a menos que os métodos e as técnicas envolvidas sejam conhecidas de maneira profunda. Só assim, é possível extrair o máximo de informações do dado sísmico, caracterizando o reservatório de forma quantitativa e integrada, otimizando sua produção e reduzindo os riscos e custos com a sua explotação / Abstract: The reservoir characterization process can be considered curretly the most important stage in the exploration, development and production of the oil field. However, this process is only carried out in the best way if the geologists, geophisicist and engineering has the knowledge of some definitive methods and techniques that integrated all information available about the field. Thus, the aim of this thesis is to study in a criterious and quantitative way the reservoirscharacterization process, analyzing the seismic data,by the evaluation of classic and novel methods, to defining new methodlogies that can be applied in decisive way into this process. So, for these methods could be evaluated in a conclusive way, were used in this thesis asynthetic reference model that reproduced some critical features of determined reservoirs, as the complex distribution of sand bodies and the subseismic thickness. These characteristics pushing to the limit the traditional modelling techniques. In this thesis to characterizze the complex features present in the reference model we used two interpretation techniques, first we analyze the seismic inversion that give a preditive character to the seismic data and after we study the multiattribute analysis that give a classificatory caracter to the seismic interpretation. For the seismic inversion, the stochastic or geostatistical inversion, that demonstrated to be the most efficient technique to characterized the complex and the subseismic features present in the model. About the seismic attributes it could be demonstrated that even so in some cases they represent the features of the model, are necessary that they are dealt with a multivariate approach, to used the advantage of the correlation between them. For the seismic attribute analysis, 3 methods of multivariate statistics analysis were used, two of them (ICA and MAF) for the first time in the reservoir characterization processo With the results we can proved that these 2 new methods improved the process of multiattribute anlysis prducing superior results when compare with the results obtained by the application of traditional PCA technique. With it was presented, can be concluded that the reservoir characterization process is a crucial stage and have some difficults to be accomplishment, unless the methods and the involved techniques are known deeply. Thus it is possible to extract the maximum informations from the datasets, characterizing the reservoir in a quantitative and integrated environmental, optimizing its production and reducing the risks and the costs with its explotation / Doutorado / Administração e Politica de Recursos Minerais / Doutor em Ciências

Orientadores: Rodrigo de Souza Portugal, Carlos Alves da Cunha Filho / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica, Instituto de Geociencias / Made available in DSpace on 2018-08-12T03:43:20Z (GMT). No. of bitstreams: 1 Moraes_DioneCherpinsky_M.pdf: 4463877 bytes, checksum: 6418a90fbf3389f9233ef9eb721da563 (MD5) Previous issue date: 2004 / Resumo: Os levantamentos sísmicos são parametrizados para que os dados sejam adquiridos segundo uma malha regular. Tal regularidade quase nunca é possível, pois durante a aquisição dos dados ocorrem obstáculos operacionais como cidades, estradas, plataformas e áreas de preservação ambiental, dentre muitos outros. Em dados marítimos sempre ocorre outro tipo de irregularidade, que é a deriva do cabo de hidrofones devido a correntes oceânicas. Tenta-se então regularizar esses dados no início do processamento sísmico, para que processos cruciais como análise de velocidades e migração tenham melhores desempenhos. Neste trabalho, a interpolação e regularização dos dados são feitas com auxílio da transformada T - p. Os algoritmos desenvolvidos utilizam a técnica do empilhamento oblíquo. Para os casos 2D e 3D, os parâmetros ideais são discutidos para que o dado retome do domínio T - P com a menor quantidade de artefatos possível. A regularização dos dados é realizada quando há deriva de cabos de hidrofones. A forma e a amplitude do sinal são preservadas quando realizam-se as transformadas T - P direta e invesa. Para o caso 2D, também são estudadas a interpolação de dados e a regularização quando ocorre um longo trecho sem informação sísmica. / Abstract: Seismic surveys are designed on the purpose that all samples collected during the seismic experiment fall on a specific regular grid. Nevertheless, this data regularity is almost impossible to achieve due to different obstacles during seismic acquisition such as constructions (cities, pipelines or other facilities), roads, platforms, preservation areas and so on. A very important non-cultural irregularity which occurs during marine seismic surveys and shall be part of our main concern is hydrophone cable drift caused by ocean currents (cable feathering). These irregularities shall be treated in the first steps of seismic data processing and data regularization can be the right tool to be used. Regularization may improve the overall performance of important steps in seismic processing like velocity analysis and migration. Data interpolation and regularization are performed using the T - P transform, with slant stack 2D and 3D algorithms. We discuss different issues in 2D and 3D data regularization using T - p transforms such as ideal parameterization to avoid artifacts and the 2D experiments related to interpolation and regularization of gaps in seismic information and cable feathering. Requirements for amplitude and phase preservation when the pair of T - P transforms is performed are also discussed. / Mestrado / Mestre em Ciências e Engenharia de Petróleo

Orientador: Martin Tygel / Tese (doutorado) - Universidade Estadual de Campinas, Facudade de Engenharia Mecanica e Instituto de Geociencias / Made available in DSpace on 2018-08-12T05:48:28Z (GMT). No. of bitstreams: 1 UniversidadeEstadualdeCampinas.FaculdadedeEngenhariaMecanica.ProgramadePos-GraduacaoemCienciaseEngenhariadePetroleo_D.pdf: 5963808 bytes, checksum: 11070db7b13faf0139f14fda5d6a4fb6 (MD5) Previous issue date: 2007 / Resumo: O método Common Reflection Surface (CRS) tem por objetivo a obtenção de seções empilhadas de qualidade e extração de parâmetros do campo de onda que caracterizam o meio geológico, através de análises de coerência e com máxima utilização da redundância contida nos dados sísmicos de multicobertura. Para a realização dessas finalidades, o método CRS utiliza "famílias" de pares de fontes e receptores arbitrariamente dispostos em torno de um ponto central e empilhamentos através de tempos de trânsito multi-paramétricos. O CRS representa um significativo aprimoramento do clássico método Common midpoint (CMP), o qual utiliza famílias de pares fonte e receptor simetricamente dispostos em torno do ponto central (denominadas famílias CMP). Além disso, o número maior de parâmetros CRS permite melhores informações sobre o meio geológico. Neste trabalho, discutimos os fundamentos e principais resultados do método CRS na situação 2D, comparando-o ao método CMP. Dentre os vários estudos e aperfeiçoamentos realizados no método CRS, destacam-se (a) nova implementação do programa CRS, com utilização de parâmetros diretamente advindos dos dados; (b) implementação conjunta dos métodos CRS e CMP visando melhor contribuição de cada um deles para o imageamento sísmico; (c) introdução de nova metodologia para eliminação ou tenuação de ruído alinhado nas seções CRS e (d) investigação qualitativa e quantitativa das resoluções vertical e lateral do método CRS. / Abstract: The Common Reflection Surface (CRS) method is designed to produce stacked sections of superior quality, as well as a number of useful wavefield attributes that caracterize the geological subsurface. These attributes are estimated by means of coherency analysis methods, which optimally exploit the redundancy contained in the multicoverge seismic data. To take advantage of this redundancy, the CRS method employs "gathers" of source and receiver pairs, arbitrarily located around a fixed central point, with the help of multi-parametric traveltime moveouts. Waving the restriction of symmetric source-receiver pairs around the central point allows full use of the multicoverage data. As a result, CRS stacked sections are seen to be cleaner, with better signal-tonoise ratio and continuity of reflection events. We discuss the basics and main results of the CRS method, in the 2D situation, in close comparison with the classical CMP method. A number of improvements on the method itself, as well as on its practical applications, are presented. These include (a) a new implementation of CRS, with the use of parameters most directly related to the input data; (b) joint use of CRS and CMP, with the aim of obtaining the best possible processing and imaging results; (c) introduction of a new methodology to eliminate or attenuate aligned noise and (d) qualitative and quantitative investigation on vertical and lateral resolution in CRS sections. / Doutorado / Reservatórios e Gestão / Doutor em Ciências e Engenharia de Petróleo

Orientadores: Maria Amélia Novais Schleicher, Lúcio Tunes dos Santos / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Matemática Estatística e Computação Científica / Made available in DSpace on 2018-08-25T10:12:39Z (GMT). No. of bitstreams: 1 Kamioka_DanielaMidori_M.pdf: 9312784 bytes, checksum: cd5e7b440a4b44168749d7cc069f1ac0 (MD5) Previous issue date: 2014 / Resumo: Uma das principais etapas do processamento sísmico é a análise de velocidade. Essa é uma das técnicas que utilizam dados sísmicos organizados em famílias de ponto médio comum (CMP ¿ Common Midpoint). Encontrar curvas (por exemplo, hiperbólicas) pré-determinadas que se ajustam aos tempos de trânsito de reflexão da melhor maneira possível é de grande importância no processo CMP. Além disso, é necessário determinar corretamente os parâmetros que definem essas curvas de melhores ajustes, pois estes parâmetros transmitem informações importantes a serem extraídos dos dados sísmicos. Portanto, é essencial uma medida que estabeleça se certa curva se ajusta aos tempos de trânsito. A função semblance mede o grau de coerência (ou alinhamento) dos traços sísmicos ao longo das curvas testadas. A semblance convencional é uma medida robusta e fácil de ser calculada, e por esses motivos tem sido a mais utilizada dentre outras medidas. Além da função convencional, aplicamos outros dois tipos de funções: a Semblance com Peso e a Semblance AB. Os exemplos numéricos desse trabalho mostram que a medida de coerência com peso exibe uma melhor resolução no sentido de que as seções apresentam picos mais destacados com relação à semblance convencional, mas pelo fato dos valores obtidos pelas duas funções serem semelhantes e o custo computacional do cálculo da medida com peso ser bem superior à convencional, o esforço não vale a pena. Já a semblance AB tenta ser superior à convencional uma vez que para dados que apresentam um trend de amplitude (reversão de polaridade), a função AB apresenta uma pequeno aumento nos seus valores na região de reversão em questão, mas nada significativo, pelo menos nos exemplos de seção CMP utilizados nessa dissertação. Apesar de todas as tentativas de se obter uma medida de coerência que apresente melhores resultados com um baixo custo computacional, ainda assim a semblance convencional continua a apresentar resultados semelhantes aos outros métodos e com o menor tempo computacional para o seu cálculo / Abstract: One of the main steps of seismic processing is the velocity analysis. This is one of the methods that uses seismic data arranged in common midpoint (CMP) gathers. Finding predetermined curves (such as hyperbolic ones) that fit the reflection traveltimes the best possible way is of great importance in the CMP process. Moreover, it is necessary to correctly determine the parameters that define those best-fitting curves, as these parameters convey important information to be extracted from seismic data. Therefore, it is essential to establish a measure that determines if any curve fits the travel time. The semblance function is such a measure that determines the degree of coherence (or alignment) of the seismic traces along trial curves. Conventional semblance is a robust and easy-to-calculate measure, and for these reasons it has been the most used among other measures. Recently, two variations of semblance have been proposed in the literature: Weighted Semblance and AB Semblance. In this thesis, we show by means of numerical examples that Weighted Semblance provides better resolution in the sense that the semblance sections present more pronounced peaks than conventional semblance. However, the values obtained by the two functions are rather similar and the computational cost of calculating the weighted measure is much higher. Therefore, the effort is not worth it. AB semblance tries to be superior to the conventional function for data presenting an amplitude trend (for example, a polarity reversal). In our numerical tests, the AB function indeed provides a small increase in its values in the region of a polarity reversal, but nothing significant, at least for the CMP section examples used in this dissertation. Despite all attempts to obtain a coherence measure that offers better results with a low computational cost, conventional semblance still maintains similar results to other methods and the lowest computational time for its calculation / Mestrado / Matematica Aplicada / Mestra em Matemática Aplicada

Orientadores: Joerg Dietrich Wilhelm Schleicher, Maria Amélia Novais Schleicher / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Matemática, Estatística e Computação Científica / Made available in DSpace on 2018-08-19T06:03:21Z (GMT). No. of bitstreams: 1 Amaro_BrunoDias_M.pdf: 4962738 bytes, checksum: a4af95a3c26d6474e07216773530662e (MD5) Previous issue date: 2011 / Resumo: Em uma migração sísmica, após o processo de propagação da onda, um importante passo para se obter uma imagem migrada é a aplicação de uma condição de imagem, que distingue entre pontos refletores e não refletores. A condição de imagem teórica desenvolvida por Claerbout usa uma deconvolução entre os campos de ondas descendente e ascendente (ou uma divisão no domínio da frequência). Uma vez que não conhecemos a posição real do refletor, esta divisão precisa ser executada em todos os pontos da imagem, causando uma certa instabilidade no processo, já que o campo de onda descendente é igual a zero em determinados pontos do meio. Pensando em tentar solucionar esse problema, diferentes técnicas de estabilização foram propostas, sendo a primeira uma correlação ao invés da deconvolução teoricamente exigida. Entretanto este processo destrói as informações de amplitude contidas nos dados. Como fica cada vez mais importante fornecer compensação de iluminação e preservação de amplitudes, particularmente, se o objetivo é recuperar com sucesso uma refletividade do meio, uma condição de imagem que destrói amplitudes é inaceitável. Por esta razão, várias formas alternativas de condições de imagem têm surgido nos últimos anos, que são avaliadas pela qualidade das amplitudes e artefatos produzidos. Neste trabalho estudamos uma série de condições de imagem com compensação de iluminação. Apresentamos também novas condições de imagem por quadrados mínimos estabilizados e as comparamos com as condições de imagem anteriormente propostas. Nossos experimentos numéricos em um simples modelo horizontal que utiliza um modelo de velocidade vertical não homogêneo e no dado Marmousi mostram que tais condições produzem resultados satisfatórios. Uma observação geral sobre todas as comparações feitas é que a condição de imagem por quadrados mínimos totais estabilizados que zera os pontos da imagem onde o denominador é de baixo valor produziu a melhor imagem com o mínimo de artefatos de migração e o mínimo de amplitudes de migração afetados. A qualidade desta imagem chega muito próximo da simples condição de imagem por correlação / Abstract: In a seismic migration, after the wave propagation process, an important step to obtain a migrated image is the application of an image condition, that distinguishes between reflector and non reflector points. The theoretical imaging condition addressed by Claerbout uses a deconvolution between downgoing and upcoming wavefields (or a division in the frequency domain). Since we do not know the actual position of the reflector, this division must be performed at all points of the image, causing instability in the process, since the downgoing wavefield is zero at certain points. Trying to remedy this problem, different stabilization techniques have been proposed, being the first a crosscorrelation rather than the theoretically required deconvolution. However, this process destroys the information contained in the data amplitudes. As it becomes increasingly important to provide compensation for illumination and amplitude recovery, particularly, if the goal is to successfully recover a reflectivity of the medium, a condition that destroys the image amplitude is unacceptable. For this reason, several alternative forms of imaging conditions have emerged in the recent past that are evaluated by the quality of the output amplitudes and artifacts produced. In this work we study a set of imaging conditions with illumination compensation. We also present new stabilized least-squares image conditions and compare them to previously proposed forms. Our numerical experiments on a simple horizontal interface model using a vertically inhomogeneous velocity model and on the Marmousi data set show that they produce satisfactory results. The general observation from the overall comparison is that the stabilized total least-squares imaging condition produced the best image with the least migration artifacts and the least affected migration amplitudes. Its image quality comes very close to the one of the simple crosscorrelation imaging condition, however with correctly preserved amplitudes / Mestrado / Matematica Aplicada / Mestre em Matemática Aplicada