Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-20T13:16:05Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_DeterminacaoVelocidadesIntervalares.pdf: 7481021 bytes, checksum: 44368dd9cbbfc6d84ad2252e68b0db37 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:12:24Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-09T11:40:45Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_DeterminacaoVelocidadesIntervalares.pdf: 7481021 bytes, checksum: 44368dd9cbbfc6d84ad2252e68b0db37 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T12:00:44Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_DeterminacaoVelocidadesIntervalares.pdf: 7481021 bytes, checksum: 44368dd9cbbfc6d84ad2252e68b0db37 (MD5) / Made available in DSpace on 2014-09-10T12:00:44Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_DeterminacaoVelocidadesIntervalares.pdf: 7481021 bytes, checksum: 44368dd9cbbfc6d84ad2252e68b0db37 (MD5) Previous issue date: 1998 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Neste trabalho foi desenvolvido um método de solução ao problema inverso para modelos sísmicos compostos por camadas homogêneas e isotrópicas separadas por superfícies suaves, que determina as velocidades intervalares em profundidade e calcula a geometria das interfaces. O tempo de trânsito é expresso como uma função de parâmetros referidos a um sistema de coordenadas fixo no raio central, que é determinada numericamente na superfície superior do modelo. Essa função é posteriormente calculada na interface anterior que limita a camada não conhecida, através de um processo que determina a função característica em profundidade. A partir da função avaliada na interface anterior se calculam sua velocidade intervalar e a geometria da superfície posterior onde tem lugar a reflexão do raio. O procedimento se repete de uma forma recursiva nas camadas mais profundas obtendo assim a solução completa do modelo, não precisando em nenhum passo informação diferente à das camadas superiores. O método foi expresso num algoritmo e se desenvolveram programas de computador, os quais foram testados com dados sintéticos de modelos que representam feições estruturais comuns nas seções geológicas, fornecendo as velocidades em profundidade e permitindo a reconstrução das interfaces. Uma análise de sensibilidade sobre os programas mostrou que a determinação da função característica e a estimação das velocidades intervalares e geometria das interfaces são feitos por métodos considerados estáveis. O intervalo empírico de aplicabilidade das correções dinâmicas hiperbólicas foi tomado como uma estimativa da ordem de magnitude do intervalo válido para a aplicação do método. / In this work a method was developed to solve the inverse seismic problem in models consisting of isotropic and homogeneous layers separated by smooth interfaces, which determines the interval velocities in depth and calculates the geometry of the interfaces. The traveltime is expressed by a function with parameters referred to a coordinated system fixed at the central ray, and numerically estimated at the superior surface of the model in the vicinity of the normal ray. The function is later recalculated at the anterior interface limiting the unknown layer, through a process which determines the characteristic function in depth. The characteristic function of the traveltimes evaluated at the anterior interface allows to know the interval velocity of the layer and the geometry of the posterior interface where the normal reflection takes place. The procedure is repeated recursively at deeper layers getting the complete solution without a priori knowledge but the upper determined layers. Computer’s programs expressing the algorithm of the method were developed and tested with synthetic seismic data, generated through models with structural factions very common in geological sections, obtaining the interval velocities in depth with considered acceptable errors and reconstructing the interfaces. A sensibility analysis was done in order to verify the stability of the two methods. The empirical range of applicability of hyperbolic dynamic corrections was taken for the range of applicability of the developed method.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-09T11:42:28Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoDeconvolucaoHomomorfica.pdf: 9295991 bytes, checksum: d8c9b4b1a029f7b919a238cb75a7235d (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar as palavras-chave on 2014-08-07T13:41:27Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-11T17:45:04Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoDeconvolucaoHomomorfica.pdf: 9295991 bytes, checksum: d8c9b4b1a029f7b919a238cb75a7235d (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-22T13:31:41Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoDeconvolucaoHomomorfica.pdf: 9295991 bytes, checksum: d8c9b4b1a029f7b919a238cb75a7235d (MD5) / Made available in DSpace on 2014-09-22T13:31:41Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoDeconvolucaoHomomorfica.pdf: 9295991 bytes, checksum: d8c9b4b1a029f7b919a238cb75a7235d (MD5) Previous issue date: 1998 / Um registro sísmico é frequentemente representado como a convolução de um pulso-fonte com a resposta do meio ao impulso, relacionada ao caminho da propagação. O processo de separação destes dois componentes da convolução é denominado deconvolução. Existe uma variedade de aproximações para o desenvolvimento de uma deconvolução. Uma das mais comuns é o uso da filtragem linear inversa, ou seja, o processamento do sinal composto, através de um filtro linear, cuja resposta de frequência é a recíproca da transformada de Fourier de um dos componentes do sinal. Obviamente, a fim de usarmos a filtragem inversa, tais componentes devem ser conhecidas ou estimadas. Neste trabalho, tratamos da aplicação a sinais sísmicos, de uma técnica de deconvolução não linear, proposta por Oppenheim (1965), a qual utiliza a teoria de uma classe de sistemas não lineares, que satisfazem um princípio generalizado de superposição, denominados de sistemas homomórficos. Tais sistemas são particularmente úteis na separação de sinais que estão combinados através da operação de convolução. O algoritmo da deconvolução homomórfica transforma o processo de convolução em uma superposição aditiva de seus componentes, com o resultado de que partes simples podem ser separadas mais facilmente. Esta classe de técnicas de filtragem representa uma generalização dos problemas de filtragem linear. O presente método oferece a considerável vantagem de que não é necessário fazer qualquer suposição prévia sobre a natureza do pulso sísmico fonte, ou da resposta do meio ao impulso, não requerendo assim, as considerações usuais de que o pulso seja de fase-mínima e que a distribuição dos impulsos seja aleatória, embora a qualidade dos resultados obtidos pela análise homomórfica seja muito sensível à razão sinal/ruído, como demonstrado. / A seismic record is often represented as the convolution of a wavelet with the impulse response relative to the reflection path. The process of separating these two components of the convolution is termed deconvolution. There are a number of approaches for carrying out a deconvolution. One of the most common is the use of linear inverse filtering, that is, processing the composite signal through a linear filter, whose frequency response is the reciprocal of the Fourier transform of one of the signal components. Obviously, in order to use inverse filtering, such components must be known or estimated. In this work, we deal with the application to seismic signals, of a nonlinear deconvolution technique, proposed by Oppenheim (1965), which uses the theory of a class of nonlinear systems, that satisfy a generalized principle of superposition, which are termed homomorphic systems. Such systems are particularly useful in separating signals which have been combined through the convolution operation. The homomorphic deconvolution algorithm transforms the convolutional process into an additive superposition of its components, with the result that the single parts can be separated more easily. This class of filtering techniques represent a generalization of linear filtering problems. This method offers the considerable advantage that no prior assumption about the nature of the seismic wavelet or the impulse response of the reflection path need be made, that is, it does not require the usual assumptions of a minimum-phase wavelet and a random distribution of impulses, although the quality of the results obtained by the homomorphic analysis is very sensible to the signal/noise ratio, as demonstrated.

O delta submarino do Amazonas é a principal feição geológica da Plataforma Continental Amazônica e é influenciado por diversos fatores oceanográficos e meteorológicos como a forte amplitude das marés, a Corrente Norte Brasileira e os ventos alíseos, e, principalmente, pela descarga do rio Amazonas. A construção do delta e a distribuição sedimentar tem sido influenciada por estes fatores ao longo do tempo e para investigar estes aspectos, foram utilizados registros de sísmica rasa de alta resolução e testemunhos de sedimento. O principal enfoque da pesquisa consistiu no processamento e na interpretação de dados de sísmica Parasound, para o qual ainda não há uma classificação específica para os tipos de eco. Os dados sísmicos foram processados utilizando-se dois softwares, o SENT da Universidade de Bremen (Alemanha) e o SEISMIC UNIX (SU) do Center for Waver Phenomena da Colorado School of Mines (Estados Unidos). Em uma primeira etapa foi feita uma filtragem de freqüências aplicando-se filtros passa-banda para se eliminar os ruídos e melhorar a visualização do sinal sísmico. Foram confeccionados espectros de amplitude do sinal para se checar a atenuação do sinal em áreas com presença de gás disperso no sedimento. Após a etapa do processamento foi realizada a interpretação dos perfis com a utilização de dados de análise granulométrica do sedimento e dos espectros de amplitude. Foram interpretados dez perfis, sendo três perpendiculares a linha de costa e o restante paralelos. A interpretação destes perfis permitiu o mapeamento das zonas com presença de gás, caracterização das feições geológicas e distinção dos diferentes tipos de caráter de eco. Os dados sísmicos foram processados utilizando-se dois softwares, o SENT da Universidade de Bremen (Alemanha) e o SEISMIC UNIX (SU) do Center for Waver Phenomena da Colorado School of Mines (Estados Unidos). Em uma primeira etapa foi feita uma filtragem de freqüências aplicando-se filtros passa-banda para se eliminar os ruídos e melhorar a visualização do sinal sísmico. Foram confeccionados espectros de amplitude do sinal para se checar a atenuação do sinal em áreas com presença de gás disperso no sedimento. Após a etapa do processamento foi realizada a interpretação dos perfis com a utilização de dados de análise granulométrica do sedimento e dos espectros de amplitude. Foram interpretados dez perfis, sendo três perpendiculares a linha de costa e o restante paralelos. A interpretação destes perfis permitiu o mapeamento das zonas com presença de gás, caracterização das feições geológicas e distinção dos diferentes tipos de caráter de eco. / Amazon submarine delta is the main feature of the Amazon shelf. This area is influenced by several energetic forces, such as the large tidal amplitudes, the North Brazilian Current, the Trade Winds and specially by the enormous water discharge from the Amazon River. The delta construction and the sedimentary distribution have been influenced by this factors during the geologic time. In order to understand these aspects, high resolution shallow seismic and sediment cores were used.The main focus of this research was seismic processing and interpretation of Parasound data to which no echo-character classification was yet established. Data was processed using two softwares, SENT of Bremen University and Seismic Unix (SU) of Center for Wave Phenomena of Colorado School of Mines. At the first stage “pass-band” filters were applied in order to reduce noise and improve the quality of seismic signal. Amplitude spectra were constructed to check the signal attenuation when gas is present in sediment. Profile nterpretation was based on sediment grain-size analyses and spectra amplitude. Ten profiles were interpreted, being three of them perpendicular to the coast and others parallel. Interpretation allowed location of gas charged sediment,characterization of geologic features and distinction of echo character.

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-05-28T13:29:48Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MigracaoKirchhoffPreEmpilhamento.pdf: 7520606 bytes, checksum: 04380685465a2a6e5ef4c422ca54216a (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir assuntos on 2014-08-06T16:21:01Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-08T14:11:40Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MigracaoKirchhoffPreEmpilhamento.pdf: 7520606 bytes, checksum: 04380685465a2a6e5ef4c422ca54216a (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-09T12:08:59Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MigracaoKirchhoffPreEmpilhamento.pdf: 7520606 bytes, checksum: 04380685465a2a6e5ef4c422ca54216a (MD5) / Made available in DSpace on 2014-09-09T12:08:59Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MigracaoKirchhoffPreEmpilhamento.pdf: 7520606 bytes, checksum: 04380685465a2a6e5ef4c422ca54216a (MD5) Previous issue date: 2007 / 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 / A teoria dos feixes gaussianos foi introduzida na literatura sísmica no início dos anos 80 por pesquisadores russos e tchecos, e foi originalmente utilizada no cálculo do campo de ondas eletromagnéticas, baseado na teoria escalar da difração. Na teoria dos feixes gaussianos, o campo de ondas sísmicas é obtido por uma integral, cujo o integrando é constituído de duas partes, a saber: (1) as amplitudes dos campos das ondas na vizinhança do ponto de observação e (2) a função fase de cada um desses campos de ondas, que neste caso é representada por um tempo de trânsito paraxial complexo. Como ferramenta de imageamento, mais precisamente como operador de migração, os primeiros trabalhos usando feixes gaussianos datam do final da década de 80 e início dos anos 90. A regularidade dos campos de ondas descritos pelos feixes gaussianos, além de sua alta precisão em regiões singulares do modelo de velocidades, tornaram o uso de feixes gaussianos como uma alternativa híbrida viável para a migração. Nesse trabalho, unimos a flexibilidade da migração tipo Kirchhoff em profundidade em verdadeira amplitude com a regularidade da descrição do campo de ondas, representado pela sobreposição de feixes gaussianos. Como forma de controlar de forma estável quantidades usadas na construção de feixes gaussianos, utilizamos informações advindas do volume de Fresnel, mais precisamente a zona de Fresnel ao redor do ponto de reflexão e a zona de Fresnel projetada, localizada ao redor do ponto de registro do sismograma e cuja a informação se encontra nas curvas de reflexão de dados sísmico. Nosso processo de migração pode ser chamado como uma migração Kirchhoff em verdadeira amplitude usando um operador de feixes gaussianos. / The Gaussian Beam (GB) concept was introduced in the seismic literature by Russian and Czech researchers in the begining of the 80’s. This theory, which by its turn was based on the scalar electromagnetic diffraction theory, is in fact a (zero order) complex paraxial ray theory, designed to satisfactorilly describe the seismic wavefield propagation beyond the standard zero order ray theory, up to then the only theory used to describe the high frequency seismic wavefield propagation in smoothed velocity models. As an imaging tool, the first works to deal with GB’s were published in the end of the 80’s and in the begining of the 90’s. The regularity in the description of the wavefield by GB’ s, as well as its high accuracy in some singular regions of the velocity model, transformed the use of GB’s into a viable hybrid alternative in the migration theory. In this work, we unite the flexibility in imaging of the true amplitude prestack Kirchhoff depth migration with the regularity in the description of the wavefield by a superposition of GB’s. As a way of controlling in a very stable way some quantities used in the construction of the beams, we have made use of some informations based on the Fresnel volume elements, more especifically speaking the Fresnel zone radius around the reflection point in depth and its counterpart, projected towards the acquisition surface. This information is centred around the recording point of the seismogram and is also present in the seismic data reflection traveltime curves. Our migration process can be named a true amplitude prestack Kirchhoff depth migration using GB’s as Green function, namely KGB-PSDM.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-26T12:49:54Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ImageamentoHomeomorficoRefletores.pdf: 21596438 bytes, checksum: f55e0802039a072693fa618437a54300 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:14:41Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-04T16:32:56Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ImageamentoHomeomorficoRefletores.pdf: 21596438 bytes, checksum: f55e0802039a072693fa618437a54300 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-16T14:44:09Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ImageamentoHomeomorficoRefletores.pdf: 21596438 bytes, checksum: f55e0802039a072693fa618437a54300 (MD5) / Made available in DSpace on 2014-09-16T14:44:09Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ImageamentoHomeomorficoRefletores.pdf: 21596438 bytes, checksum: f55e0802039a072693fa618437a54300 (MD5) Previous issue date: 1994 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / DAAD - Serviço Alemão de Intercâmbio Acadêmico / Deutscher Akademischer Austauschdiens / Neste trabalho é apresentada uma nova técnica para a realização do empilhamento sísmico, aplicada ao problema do imageamento de refletores fixos em um meio bidimensional, suavemente heterogêneo, isotrópico, a partir de dados de reflexão. Esta nova técnica chamada de imageamento homeomórfico tem como base a aproximação geométrica do raio e propriedades topológicas dos refletores. São utilizados, portanto, os conceitos de frente de onda, ângulo de incidência, raio de curvatura da frente de onda, cáustica e definição da trajetória do raio; de tal modo que a imagem obtida mantém relações de homeomorfismo com o objeto que se deseja imagear. O empilhamento sísmico é feito, nesta nova técnica de imageamento, aplicando-se uma correção local do tempo, ∆ t, ao tempo de trânsito, t, do raio que parte da fonte sísmica localizada em xo, reflete-se em um ponto de reflexão, Co, sendo registrado como uma reflexão primária em um geofone localizado em xg, em relação ao tempo de referência to no sismograma, correspondente ao tempo de trânsito de um raio central. A fórmula utilizada nesta correção temporal tem como parâmetros o raio de curvatura Ro, o ângulo de emergência βo da frente de onda, no instante em que a mesma atinge a superfície de observação, e a velocidade vo considerada constante nas proximidades da linha sísmica. Considerando-se uma aproximação geométrica seguido um círculo para a frente de onda, pode-se estabelecer diferentes métodos de imageamento homeomórfico dependendo da configuração de processamento. Sendo assim tem-se: 1) Método Elemento de Fonte (Receptor) Comum (EF(R)C). Utiliza-se uma configuração onde se tem um conjunto de sismogramas relacionado com uma única fonte (receptor), e considera-se uma frente de onda real (de reflexão); 2) Método Elemento de Reflexão Comum (ERC). Utiliza-se uma configuração onde um conjunto de sismogramas é relacionado com um único ponto de reflexão, e considera-se uma frente de onda hipoteticamente originada neste ponto; 3) Método Elemento de Evoluta Comum (EEC). Utiliza-se uma configuração onde cada sismograma está relacionado com um par de fonte e geofone coincidentemente posicionados na linha sísmica, e considera-se uma frente de onda hipoteticamente originada no centro de curvatura do refletor. Em cada um desses métodos tem-se como resultados uma seção sísmica empilhada, u(xo, to); e outras duas seções denominadas de radiusgrama, Ro (xo, to), e angulograma, βo(xo, to), onde estão os valores de raios de curvatura e ângulos de emergência da frente de onda considerada no instante em que a mesma atinge a superfície de observação, respectivamente. No caso do método denominado elemento refletor comum (ERC), a seção sísmica resultante do empilhamento corresponde a seção afastamento nulo. Pode-se mostrar que o sinal sísmico não sofre efeitos de alongamento como consequência da correção temporal, nem tão pouco apresenta problemas de dispersão de pontos de reflexão como consequência da inclinação do refletor, ao contrário do que acontece com as técnicas de empilhamento que tem por base a correção NMO. Além disto, por não necessitar de um macro modelo de velocidades a técnica de imageamento homeomórfico, de um modo geral, pode também ser aplicada a modelos heterogêneos, sem perder o rigor em sua formulação. Aqui também são apresentados exemplos de aplicação dos métodos elemento de fonte comum (EFC) (KEYDAR, 1993), e elemento refletor comum (ERC) (STEENTOFT, 1993), ambos os casos com dados sintéticos. No primeiro caso, (EFC), onde o empilhamento é feito tendo como referência um raio central arbitrário, pode-se observar um alto nível de exatidão no imageamento obtido, além do que é dada uma interpretação para as seções de radiusgrama e angulograma, de modo a se caracterizar aspectos geométricos do model geofísico em questão. No segundo caso, (ERC), o método é aplicado a série de dados Marmousi, gerados pelo método das diferenças finitas, e o resultado é comparado com aquele obtido por métodos convecionais (NMO/DMO) aplicados aos mesmos dados. Como consequência, observa-se que através do método ERC pode-se melhor detectar a continuidade de refletores, enquanto que através dos métodos convencionais caracterizam-se melhor a ocorrência de difrações. Por sua vez, as seções de radiusgrama e angulograma, no método (ERC), apresentam um baixo poder de resolução nas regiões do modelo onde se tem um alto grau de complexidade das estruturas. Finalmente, apresenta-se uma formulação unificada que abrange os diferentes métodos de imageamento homeomórfico citados anteriormente, e também situações mais gerais onde a frente de onda não se aproxima a um círculo, mas a uma curva quadrática qualquer. / This thesis presents a new technique for seismic stacking called homeomorphic imaging, which is applicable to the imaging of seismic reflectors in a bidimensional, inhomogeneous and isotropic medium. This new technique is based on ray geometrical approximation and topological properties of reflection surfaces. For this purpose the concepts of wavefront, incidence angle, radius and caustic of wavefront and ray trajetory are used. Considering a circle as the geometrical approximation of the wavefront in propagation, it is possible to define diferent homeomorphic imaging methods, depending on processing configuration. In this way, the following methods are possible: 1) Common Source (Receiver) Element (CS(R)E), which relate to a set of seismograms with a single source (receiver) and a real reflected wavefront is considered; 2) Common-Reflecting-Element (CRE), which relate to a set of seismograms with a single reflection point and a wavefront hipotetically generated in the same reflection point is considered; 3) Common Evolute Element (CEE), which relate to a set of seismograms with each pair of source and geophone located in the same point on the seismic line and a wavefront hipothetically generated in the curvature center of the reflector is considered. In the first method is obtained a stacked seismic section using arbitrary central rays. In the last two methods the result is a zero-offset seismic section. These methods give also other two sections called radiusgram and anglegram, the latter being emergence angles and the former radii of wavefront in the moment that it reaches the observational surface. The seismic stacking is made using a local correction-time applied to the travel time of a ray that leaves the source, and after reflection, is registered as a primary reflection at a geophone, in relation to the reference time which is the travel time of the central ray. The formula used for the temporal correction depends on the radius, the emergence angle of the wavefront and the velocity which is considered constant near the seismic line. It is possible to show that in this new technique the registered signal is not submitted to stretch effects as a consequence of the temporal correction, furthermore there is no problem with reflector point dispersal as a consequence of dip reflectors, in contrast with the techniques that are based on NMO/DMO. In addition, considering that no a prori knowledge of a macromodel is necessary but the velocity near the seismic line, the homeomorphic imaging can be applied to inhomogeneous models without losing the strictness of the formulation.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-09T12:19:45Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoMetodoKalman.pdf: 5931360 bytes, checksum: 9d5bc90e28509e1df9337a46bafe80c4 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-06T14:42:05Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-11T13:45:39Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoMetodoKalman.pdf: 5931360 bytes, checksum: 9d5bc90e28509e1df9337a46bafe80c4 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-19T17:03:06Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoMetodoKalman.pdf: 5931360 bytes, checksum: 9d5bc90e28509e1df9337a46bafe80c4 (MD5) / Made available in DSpace on 2014-09-19T17:03:06Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoMetodoKalman.pdf: 5931360 bytes, checksum: 9d5bc90e28509e1df9337a46bafe80c4 (MD5) Previous issue date: 1998 / O filtro de Kalman é aplicado para filtragem inversa ou problema de deconvolução. Nesta dissertação aplicamos o método de Kalman, considerado como uma outra visão de processamento no domínio do tempo, para separar sinal-ruído em perfil sônico admitido como uma realização de um processo estocástico não estacionário. Em um trabalho futuro estudaremos o problema da deconvolução. A dedução do filtro de Kalman destaca a relação entre o filtro de Kalman e o de Wiener. Estas deduções são baseadas na representação do sistema por variáveis de estado e modelos de processos aleatórios, com a entrada do sistema linear acrescentado com ruído branco. Os resultados ilustrados indicam a aplicabilidade dessa técnica para uma variedade de problemas de processamento de dados geofísicos, por exemplo, ideal para well log. O filtro de Kalman oferece aos geofísicos de exploração informações adicionais para o processamento, problemas de modelamento e a sua solução. / The Kalman filter is applied to the inverse filtefing or deconvolution problem. In this dissertation we applied the Kalman method, it is considered like a processament vition on time domain, to separet signal-noise within sonic perfil which is admited like no stationary stochastic process. In next work will survey deconvolution problem. The derivation given of the Kalman filter emphasizes the relationship between the Kalman and Wiener filter. This derivation is based on the modeling of randon processes as the output of linear systems excited by white noise. Ilustrative results indicate the applicability of these tchniques to a variety of geophysical data processing problems, for example the ideal well log teated here. The Kalman filter offters exploration geophysicists addtition insight into processing problem modeling and solution.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-07-07T13:24:07Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EmpilhamentoSismicoComposicao.pdf: 15271732 bytes, checksum: 07aae8678ad182a0825a6c915a960864 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:25:09Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-08T14:38:21Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EmpilhamentoSismicoComposicao.pdf: 15271732 bytes, checksum: 07aae8678ad182a0825a6c915a960864 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-19T14:43:56Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EmpilhamentoSismicoComposicao.pdf: 15271732 bytes, checksum: 07aae8678ad182a0825a6c915a960864 (MD5) / Made available in DSpace on 2014-09-19T14:43:56Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EmpilhamentoSismicoComposicao.pdf: 15271732 bytes, checksum: 07aae8678ad182a0825a6c915a960864 (MD5) Previous issue date: 1997 / DAAD - Serviço Alemão de Intercâmbio Acadêmico / Deutscher Akademischer Austauschdiens / No presente trabalho de tese é apresentada uma nova técnica de empilhamento de dados sísmicos para a obtenção da seção de incidência normal ou afastamento fonte-receptor nulo, aplicável em meios bidimensionais com variações laterais de velocidade. Esta nova técnica denominada Empilhamento Sísmico pela Composição de Ondas Planas (empilhamento PWC) foi desenvolvida tomando como base os conceitos físicos e matemáticos da decomposição do campo de ondas em ondas planas. Este trabalho pode ser dividido em três partes: Uma primeira parte, onde se apresenta uma revisão da técnica de empilhamento sísmico convencional e do processo de decomposição do campo de ondas produzido a partir de fontes pontuais em suas correspondentes ondas planas. Na segunda parte, é apresentada a formulação matemática e o procedimento de aplicação do método de empilhamento sísmico pela composição de ondas planas. Na terceira parte se apresenta a aplicação desta nova técnica de empilhamento na serie de dados Marmousi e uma analise sobre a atenuação de ruído. A formulação matemática desta nova técnica de empilhamento sísmico foi desenvolvida com base na teoria do espalhamento aplicado a ondas sísmicas sob a restrição do modelo de aproximação de Born. Nesse sentido, inicialmente se apresenta a determinação da solução da equação de onda caustica para a configuração com afastamento fonte-receptor finito, que posteriormente é reduzido para a configuração de afastamento fonte-receptor nulo. Por outra parte, com base nessas soluções, a expressão matemática deste novo processo de empilhamento sísmico é resolvida dentro do contexto do modelo de aproximação de Born. Verificou-se que as soluções encontradas por ambos procedimentos, isto é, por meio da solução da equação da onda e pelo processo de empilhamento proposto, são iguais, mostrando-se assim que o processo de empilhamento pela composição de ondas planas produz uma seção com afastamento fonte-receptor nulo. Esta nova técnica de empilhamento basicamente consiste na aplicação de uma dupla decomposição do campo de ondas em onda planas por meio da aplicação de dois empilhamentos oblíquos (slant stack), isto é um ao longo do arranjo das fontes e outro ao longo do arranjo dos detectores; seguido pelo processo de composição das ondas planas por meio do empilhamento obliquo inverso. Portanto, com base nestas operações e com a ajuda de um exemplo de aplicação nos dados gerados a partir de um modelo simples, são descritos os fundamentos e o procedimento de aplicação (ou algoritmo) desta nova técnica de obtenção da seção de afastamento nulo. Como exemplo de aplicação do empilhamento PWC em dados correspondentes a um meio com variações laterais de velocidade, foi aplicado nos dados Marmousi gerados segundo a técnica de cobertura múltipla a partir de um modelo que representa uma situação geológica real. Por comparação da seção resultante com a similar produzida pelo método de empilhamento convencional, observa-se que a seção de afastamento nulo desta nova técnica apresenta melhor definição e continuidade dos reflectores, como também uma melhor caracterização da ocorrência de difrações. Por último, da atenuação de ruído aleatório realizada nos mesmos dados, observa-se que esta técnica de empilhamento também produz uma atenuação do ruído presente no sinal, a qual implica um aumento na relação sinal ruído. / In this thesis we present a new seismic data stacking method called Plane Wave Composition (PWC). This method, applicable in a bidimensional medium with lateral velocity gradients, is developed on the basis of physical and mathematical concepts on the plane wave decomposition of spherical wave fields. In the initial part of this work, we present a review on the conventional stacking method and on plane wave decomposition of the point-source seismograms. The stacking by plane wave composition is a method which produce a normal incidence (or zero offset) section by the application of the following main processes: A double plane wave decomposition, achieved by a slant stack along the shot array and another slant stack along the receiver array, followed by a plane wave composition achieved by an inverse slant stack. The PWC stacking method is theoretically formulated here on the basis of the scattering theory applied to seismic waves, within the constraint of the Born approximation. Initially, starting with the acoustic wave equation, for a finite source-receiver configuration, a solution for the direct single scatter (Born) model is derived. That result is reduced for the coincident source-receiver (zero offset) configuration. Afterward, the mathematical expression of PWC stacking method is solved replacing the observed data function by the scattered field obtained by the Born approximation. For most clarity, the algorithm to obtain the zero offset seismic section by the PWC stacking method is described by applying it to the data corresponding to a simple model. A successful application is performed using the Marmousi seismic data set, corresponding to a geological complex model. Finally, in the same data set, a noise analysis shows that this method increase the signal-noise ratio in the seismic trace. Thus, it has been showed that the PWC stacking method is an efficient alternative to process seismic data of complex models.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-05T13:44:13Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_IdentificacaoAutomaticaPrimeiras.pdf: 11576524 bytes, checksum: ad81c2874a0a83c3c5bd73b07a0c091c (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-06T13:10:44Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-11T13:31:43Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_IdentificacaoAutomaticaPrimeiras.pdf: 11576524 bytes, checksum: ad81c2874a0a83c3c5bd73b07a0c091c (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-22T11:21:58Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_IdentificacaoAutomaticaPrimeiras.pdf: 11576524 bytes, checksum: ad81c2874a0a83c3c5bd73b07a0c091c (MD5) / Made available in DSpace on 2014-09-22T11:21:58Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_IdentificacaoAutomaticaPrimeiras.pdf: 11576524 bytes, checksum: ad81c2874a0a83c3c5bd73b07a0c091c (MD5) Previous issue date: 2000 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Apesar do avanço tecnológico ocorrido na prospecção sísmica, com a rotina dos levantamentos 2D e 3D, e o significativo aumento na quantidade de dados, a identificação dos tempos de chegada da onda sísmica direta (primeira quebra), que se propaga diretamente do ponto de tiro até a posição dos arranjos de geofones, permanece ainda dependente da avaliação visual do intérprete sísmico. O objetivo desta dissertação, insere-se no processamento sísmico com o intuito de buscar um método eficiente, tal que possibilite a simulação computacional do comportamento visual do intérprete sísmico, através da automação dos processos de tomada de decisão envolvidos na identificação das primeiras quebras em um traço sísmico. Visando, em última análise, preservar o conhecimento intuitivo do intérprete para os casos complexos, nos quais o seu conhecimento será, efetivamente, melhor aproveitado. Recentes descobertas na tecnologia neurocomputacional produziram técnicas que possibilitam a simulação dos aspectos qualitativos envolvidos nos processos visuais de identificação ou interpretação sísmica, com qualidade e aceitabilidade dos resultados. As redes neurais artificiais são uma implementação da tecnologia neurocomputacional e foram, inicialmente, desenvolvidas por neurobiologistas como modelos computacionais do sistema nervoso humano. Elas diferem das técnicas computacionais convencionais pela sua habilidade em adaptar-se ou aprender através de uma repetitiva exposição a exemplos, pela sua tolerância à falta de alguns dos componentes dos dados e pela sua robustez no tratamento com dados contaminados por ruído. O método aqui apresentado baseia-se na aplicação da técnica das redes neurais artificiais para a identificação das primeiras quebras nos traços sísmicos, a partir do estabelecimento de uma conveniente arquitetura para a rede neural artificial do tipo direta, treinada com o algoritmo da retro-propagação do erro. A rede neural artificial é entendida aqui como uma simulação computacional do processo intuitivo de tomada de decisão realizado pelo intérprete sísmico para a identificação das primeiras quebras nos traços sísmicos. A aplicabilidade, eficiência e limitações desta abordagem serão avaliadas em dados sintéticos obtidos a partir da teoria do raio. / In spite of the technologic development happened at seismic prospection, and the significative amount of data with seismic two-dimensional (2D) and three-dimensional (3D) surveys, some process in the seismic interpretation task like the first break picking, remains in a manual version, that still needs an intuitive human intervention. This dissertation purpose, fill in the seismic processing with the intention to look for an efficient method to enable the computational simulation of the human visual system behavior, through decision process automation involved in first break picking in a seismic trace; looking at to preserv the interpreter intuitive knowledgement to more complex tasks, where your knowledgement will be better profitable. Neural networks, the most important implementation of neurocomputing systems, were initially developed by neurobiologists as computer models of the neural system in the brain. They differ from conventional computation techniques in their ability to adaptively discriminate or learn through repeated exposure to examples, their tolerance to data component failure and their robustness in the presence of high noise levels. This computing technology provide some techniques that can reduce the labor intensive aspects of the first break picking, maintaining the quality and reliability of the results. The method here presented is an application of an artificial neural network computational process, known as feedforward multilayer perceptron trained with the error back-propagation algorithm; from the establishment of a convenient neural network architecture and learning set that make possible its application over seismic data. This method is a computational simulation of seismic interpreter decision intuitive process for first break picking in seismic traces. The applicability, efficiency and limitations of this approach will be appraised in synthetic data obtained starting out the ray theoretical method.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-28T14:33:29Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MigracaoPreEmpilhamento.pdf: 3640959 bytes, checksum: 76a4c529f9749947cfa0045c3df999c5 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:19:11Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-08T13:13:54Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MigracaoPreEmpilhamento.pdf: 3640959 bytes, checksum: 76a4c529f9749947cfa0045c3df999c5 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-08T15:57:26Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MigracaoPreEmpilhamento.pdf: 3640959 bytes, checksum: 76a4c529f9749947cfa0045c3df999c5 (MD5) / Made available in DSpace on 2014-09-08T15:57:26Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MigracaoPreEmpilhamento.pdf: 3640959 bytes, checksum: 76a4c529f9749947cfa0045c3df999c5 (MD5) Previous issue date: 2012 / FAPESPA - Fundação Amazônia de Amparo a Estudos e Pesquisas / O método de migração do tipo Kirchhoff se apresenta na literatura como uma das ferramentas mais importantes de todo o processamento sísmico, servindo de base para a resolução de outros problemas de imageamento, devido ao um menor custo computacional em relação aos métodos que tem por base a solução numérica da equação da onda. No caso da aplicação em três dimensões (3D), mesmo a migração do tipo Kirchhoff torna-se dispendiosa, no que se refere aos requisitos computacionais e até mesmo numéricos para sua efetiva aplicação. Desta maneira, no presente trabalho, objetivando produzir resultados com uma razão sinal/ruído maior e um menor esforço computacional, foi utilizado uma simplificação do meio denominado 2.5D, baseado nos fundamentos teóricos da propagação de feixes gaussianos. Assim, tendo como base o operador integral com feixes gaussianos desenvolvido por Ferreira e Cruz (2009), foi derivado um novo operador integral de superposição de campos paraxiais (feixes gaussianos), o mesmo foi inserido no núcleo do operador integral de migração Kirchhoff convencional em verdadeira amplitude, para a situação 2,5D, definindo desta maneira um novo operador de migração do tipo Kirchhoff para a classe pré-empilhamento em verdadeira amplitude 2.5D (KGB,do inglês Kirchhoff-Gausian-Beam). Posteriormente, tal operador foi particularizado para as configurações de medida afastamento comum (CO, do inglês common offset) e ângulo de reflexão comum (CA, do inglês common angle), ressaltando ainda, que na presente Tese foi também idealizada uma espécie de flexibilização do operador integral de superposição de feixes gaussianos, no que concerne a sua aplicação em mais de um domínio, quais sejam, afastamento comum e fonte comum. Nesta Tese são feitas aplicações de dados sintéticos originados a partir de um modelo anticlinal. / A Kirchhoff-type migration is considered in the geophysics literature as one of the most fundamental tools in seismic data processing, the base for solution of several imaging problems. In this respect, it must be considered its wide use and its successful history for the oil and gas industry, associated with its low computational cost and flexibility to deal with non-wavefield datasets when compared to other methods. However in 3D, even when compared to other existing and most effective methods, its computational cost and implementation is still considered high, due to several reasons: new acquisition technologies, data storage and burden, azimuth richness, etc. Thus the main objective of the present work is to implement and simulate migration results (i.e., images) with high signal-to-noise ratios and with a less computer burdens in 2.5D media, using the theoretical framework of Gaussian Beams (GBs). By considering one implementation of a superposition of GBs integral operator studied by Ferreira and Cruz (2009) and by the use of the stationary phase method (Bleistein, 2000), a new integral superposition migration operator using paraxial fields (i.e., GBs) was implemented and studied. Theoretically speaking, the present migration operator was inserted in the kernel of a conventional, 2.5D, true-amplitude, prestack Kirchhoff migration integral operator, thus defining a 2.5D prestack Kirchhoff-Gaussian Beam (KGB) migration operator. The present migration operator was later configured to hold commonoffset (CO) and common-angle (CA) seismic acquisition configurations. I remark that in the present thesis one flexibility of the GB migration operator was idealized in order to handle its effective application in more the one sorting configuration, i.e., common-offset and commonsource.