Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-05T14:17:45Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseEsquemaLinear.pdf: 958488 bytes, checksum: d78a32e1815ce14cc9e903e1c610c694 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Incluir palavras-chave on 2014-08-07T14:17:45Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-11T12:42:20Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseEsquemaLinear.pdf: 958488 bytes, checksum: d78a32e1815ce14cc9e903e1c610c694 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-15T16:01:43Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseEsquemaLinear.pdf: 958488 bytes, checksum: d78a32e1815ce14cc9e903e1c610c694 (MD5) / Made available in DSpace on 2014-09-15T16:01:43Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseEsquemaLinear.pdf: 958488 bytes, checksum: d78a32e1815ce14cc9e903e1c610c694 (MD5) Previous issue date: 2013 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Neste trabalho é apresentada uma análise do esquema de inversão linear para a estimativa de anisotropia na vizinhança de um receptor situado em um poço a partir de da componente vertical do vetor de vagarosidade e do vetor de polarização de ondas P medidops em experimentos de VSP walkaway multiazimutal. Independente do meio acima do geofone (homogêneo ou heterogêneo) e da forma do poço (pode ser direcional ou curvado, vertical e inclinado), a inversão é feita a partir de uma aproximação de primeira ordem em torno de um meio isotrópico de referência. O esquma da inversão é analisado considerando fatores como: o nível de ruído nos dados, o tipo de onda P, o grau de anisotropia do meio, a escolha dos parâmetros no meio isotrópico de referência e grau de heterogeneidade do meio. Os resultados são apresentados. / This study presents an analysis of the linear inversion scheme for estimating anisotropy in the neighborhood of a receiver placed in a well using the vertical components of the slowness and polarization vectors of P-waves measured in multi-azimuth walkaway vertical seismic profile (VSP) surveys. Independently of the medium above the geophone (homogeneous or heterogeneous) and the shape of the well (directional or curved, vertical or sloped), an inversion is performed from a first-order approximation around a reference isotropic medium. The analysis of the inversion scheme considers factors such as the noise level of the data, the type of P-wave, the degree of the anisotropy of the medium, the choice of parameters in the reference isotropic medium and the degree of heterogeneity of the medium.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-12T11:26:19Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_CaracterizacaoReservatoriosFraturados.pdf: 1892695 bytes, checksum: c9b2c5a8a38811f366576b899c2b9603 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:23:47Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-04T13:34:51Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_CaracterizacaoReservatoriosFraturados.pdf: 1892695 bytes, checksum: c9b2c5a8a38811f366576b899c2b9603 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-15T16:04:02Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_CaracterizacaoReservatoriosFraturados.pdf: 1892695 bytes, checksum: c9b2c5a8a38811f366576b899c2b9603 (MD5) / Made available in DSpace on 2014-09-15T16:04:02Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_CaracterizacaoReservatoriosFraturados.pdf: 1892695 bytes, checksum: c9b2c5a8a38811f366576b899c2b9603 (MD5) Previous issue date: 2008 / ANP - Agência Nacional do Petróleo / Esta dissertação apresenta um método para estimar a orientação de fratura de um meio partindo do pressuposto que este se comporta efetivamente como um meio transversalmente isotrópico (TI). Desta forma, a orientação de fratura é obtida a partir da estimativa do eixo de simetria de um meio TI. Esta estimativa é obtida a partir dos vetores de vagarosidade e polarização de ondas qP medidos em experimentos VSP Walkaway. O processo de inversão é baseado em equações linearizadas dos vetores de vagarosidade e polarização de ondas qP e os parâmetros fracamente anisotrópicos do meio. São apresentados testes numéricos, em que a sensibilidade a fatores como grau de anisotropia, geometria do levantamento, tipo de onda utilizado e nível de ruído são analisados. São também apresentados os resultados para um conjunto de dados reais. / This dissertation presents a method to estimate the orientation of the axis of symmetry of a medium assuming it behaves effectively as a transversely isotropic medium (TI). The fracture orientation is then achieved from the estimation of the axis of symmetry of a TI medium. This estimative is done with the slowness and polarization vectors from qP waves, measured in VSP Walkaway experiments. The inversion process is based on linearized equations of the slowness and polarization vectors from qP waves and the weakly anisotropic parameters of the medium. Numerical tests are presented, on which the sensibility to factors like the strength of anisotropy, survey geometry, type of wave utilized and noise level are analyzed. Test results for a set of real data are also shown.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-11T11:56:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoDadosAltimetria.pdf: 11101007 bytes, checksum: 7977990454c0f6748cccf8baf41fdf95 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Ausência de palavras-chave on 2014-08-07T12:35:54Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-16T12:25:14Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoDadosAltimetria.pdf: 11101007 bytes, checksum: 7977990454c0f6748cccf8baf41fdf95 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-16T13:20:55Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoDadosAltimetria.pdf: 11101007 bytes, checksum: 7977990454c0f6748cccf8baf41fdf95 (MD5) / Made available in DSpace on 2014-09-16T13:20:55Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoDadosAltimetria.pdf: 11101007 bytes, checksum: 7977990454c0f6748cccf8baf41fdf95 (MD5) Previous issue date: 1992 / O objetivo principal deste trabalho é avaliar o uso dos dados de altimetria de satélites para mapear a superfície do potencial gravitacional (geóide) no mar. Esta avaliação se faz por comparações da resolução e precisão entre os dados de altimetria processados numa superfície equipotencial (o mar) e dados obtidos a partir de levantamentos convencionais. Uma vez processada a superfície equipotencial, quantidades tais como a anomalia "ar livre" juntamente com o desvio vertical podem ser calculadas. Os dados altimétricos ("altura do mar") utilizados neste trabalho foram coletados pelo satélite GEOSAT. Este satélite rastreou diversas áreas oceânicas do globo processando 44 ciclos em dois anos. Alguns pesquisadores utilizaram os valores médios da "altura do mar" deste satélite para melhoramentos em precisão e resolução dos registros. Estes valores tratados estão disponíveis em NOAA (National Oceanic and Atmospheric Administration) sendo deste modo repassados à UFPa para utilização nesta tese. Os dados de gravimetria marinha utilizados neste trabalho são aqueles obtidos do levantamento "Equatorial Atlantic" (EQUANT I e II) resultantes de uma pesquisa conjunta entre várias instituições com objetivos científicos de conhecer o comportamento da margem equatorial Brasileira. Para comparação e integração entre os dois tipos de dados obtidos através de fontes distintas (medidas de satélite e do navio), poder-se-ia obter a aceleração vertical numa superfície equipotencial partindo-se de um tratamento algébrico dos dados coletados por rastreamento altimétrico do satélite Geosat ou alternativamente poder-se-ia processar transformações dos dados de gravimetria marinha em uma superfície equipotencial equivalente. Em decorrência de diferenças no espaçamento entre as linhas dos levantamentos, ou seja, as linhas das trajetórias do satélite estão largamente espaçadas em comparação com aquelas do levantamento marinho, optou-se por transformar os dados gravimétricos do navio em uma superfície equipotencial. Neste tipo de transformação, consideraram-se vários fatores tais como efeitos "aliasing", nível de ruídos nos levantamentos do navio, redução ao geóide (correção "ar livre"), bem como erros computacionais durante as transformações. Com a supressão parcial desses efeitos (enfatizando o "aliasing") encontrou-se forte correlação entre os dois conjuntos de dados verificando-se um nível de coerência satisfatório a partir do comprimento de onda de 11 km. Comparando-se este resultado com o nível de resolução do satélite Geosat largamente estudado por outros pesquisadores enfatiza-se que de fato a resolução dos valores médios (2 anos) do satélite Geosat aproxima-se da resolução dos levantamentos do Equant I e Equant II. / The principal objective of this thesis is to evaluate the use of satellite altimetry data to map the potential surface of the earth's gravity field at sea. This evaluation is done by comparing the resolution and accuracy of the altimeter data over an equipotential surface (the sea) with that computed from conventional sea survey data. From the potential surface, quantities such as the free air anomaly or the vertical deflection can be calculated. The altimetry data utilized in this thesis were collected by the Geosat satellite and each data point was averaged over a two year (44 orbital cycles) period. This averaging was done to improve the accuracy and resolution of the sea height over ocean areas. These data are available from NOAA and were supplied to UFPa for use in this thesis. The ship survey data from EQUANT 1 and EQUANT 2 were generated from a joint research project among several institutions whose scientific objective was to study the behavior of the Brazilian Equatorial Continental Margin. In order to integrate and compare the data between these two different sources (satellite and ship measurements), either the satellite data must be reduced to a vertical acceleration on the equipotential surface or alternately the ship survey data must be transformed into an equivalent potential surface. Taking into account the wide separation of the satellite tracks relative to the ship survey track, it was decided to transform the ship survey data into a potential surface, and the correlation and comparisons were carried out along satellite tracks. In this type of transformation, various factors must be considered such as the effects of aliasing, ship survey noise levels, reduction to the geoid, computational errors in the transformation, etc. In spite of these effects (particularly the aliasing problem) the results reflect a strong correlation between the two sets of data verifying that a satisfactory level of coherence exists for wavelengths greater than 11 km. The resolution for the 2 year cycle average of the satellite is, in fact, almost as good as the ship survey data.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-30T16:18:05Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_RegiaoEspacoInfluencia.pdf: 18156456 bytes, checksum: bb858a842e65e92f7d3ec93381566520 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:32:09Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-16T14:13:13Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_RegiaoEspacoInfluencia.pdf: 18156456 bytes, checksum: bb858a842e65e92f7d3ec93381566520 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-16T14:35:30Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_RegiaoEspacoInfluencia.pdf: 18156456 bytes, checksum: bb858a842e65e92f7d3ec93381566520 (MD5) / Made available in DSpace on 2014-09-16T14:35:30Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_RegiaoEspacoInfluencia.pdf: 18156456 bytes, checksum: bb858a842e65e92f7d3ec93381566520 (MD5) Previous issue date: 1994 / UFPA - Universidade Federal do Pará / PETROBRAS - Petróleo Brasileiro S.A. / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / FINEP - Financiadora de Estudos e Projetos / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Localizar em subsuperfície a região que mais influencia nas medidas obtidas na superfície da Terra é um problema de grande relevância em qualquer área da Geofísica. Neste trabalho, é feito um estudo sobre a localização dessa região, denominada aqui zona principal, para métodos eletromagnéticos no domínio da freqüência, utilizando-se como fonte uma linha de corrente na superfície de um semi-espaço condutor. No modelo estudado, tem-se, no interior desse semi-espaço, uma heterogeneidade na forma de camada infinita, ou de prisma com seção reta quadrada e comprimento infinito, na direção da linha de corrente. A diferença entre a medida obtida sobre o semi-espaço contendo a heterogeneidade e aquela obtida sobre o semi-espaço homogêneo, depende, entre outros parâmetros, da localização da heterogeneidade em relação ao sistema transmissor-receptor. Portanto, mantidos constantes os demais parâmetros, existirá uma posição da heterogeneidade em que sua influência é máxima nas medidas obtidas. Como esta posição é dependente do contraste de condutividade, das dimensões da heterogeneidade e da freqüência da corrente no transmissor, fica caracterizada uma região e não apenas uma única posição em que a heterogeneidade produzirá a máxima influência nas medidas. Esta região foi denominada zona principal. Identificada a zona principal, torna-se possível localizar com precisão os corpos que, em subsuperfície, provocam as anomalias observadas. Trata-se geralmente de corpos condutores de interesse para algum fim determinado. A localização desses corpos na prospecção, além de facilitar a exploração, reduz os custos de produção. Para localizar a zona principal, foi definida uma função Detetabilidade (∆), capaz de medir a influência da heterogeneidade nas medidas. A função ∆ foi calculada para amplitude e fase das componentes tangencial (Hx) e normal (Hz) à superfície terrestre do campo magnético medido no receptor. Estudando os extremos da função ∆ sob variações de condutividade, tamanho e profundidade da heterogeneidade, em modelos unidimensionais e bidimensionais, foram obtidas as dimensões da zona principal, tanto lateralmente como em profundidade. Os campos eletromagnéticos em modelos unidimensionais foram obtidos de uma forma híbrida, resolvendo numericamente as integrais obtidas da formulação analítica. Para modelos bidimensionais, a solução foi obtida através da técnica de elementos finitos. Os valores máximos da função ∆, calculada para amplitude de Hx, mostraram-se os mais indicados para localizar a zona principal. A localização feita através desta grandeza apresentou-se mais estável do que através das demais, sob variação das propriedades físicas e dimensões geométricas, tanto dos modelos unidimensionais como dos bidimensionais. No caso da heterogeneidade condutora ser uma camada horizontal infinita (caso 1D), a profundidade do plano central dessa camada vem dada pela relação po = 0,17 δo, onde po é essa profundidade e δo o "skin depth" da onda plana (em um meio homogêneo de condutividade igual à do meio encaixante (σ1) e a freqüência dada pelo valor de w em que ocorre o máximo de ∆ calculada para a amplitude de Hx). No caso de uma heterogeneidade bidimensional (caso 2D), as coordenadas do eixo central da zona principal vem dadas por do = 0,77 r0 (sendo do a distância horizontal do eixo à fonte transmissora) e po = 0,36 δo (sendo po a profundidade do eixo central da zona principal), onde r0 é a distância transmissor-receptor e δo o "skin depth" da onda plana, nas mesmas condições já estipuladas no caso 1D. Conhecendo-se os valores de r0 e δo para os quais ocorre o máximo de ∆, calculado para a amplitude de Hx, pode-se determinar (do, po). Para localizar a zona principal (ou, equivalentemente, uma zona condutora anômala em subsuperfície), sugere-se um método que consiste em associar cada valor da função ∆ da amplitude de Hx a um ponto (d, p), gerado através das relações d = 0,77 r e p = 0,36 δ, para cada w, em todo o espectro de freqüências das medidas, em um dado conjunto de configurações transmissor-receptor. São, então, traçadas curvas de contorno com os isovalores de ∆ que vão convergir, na medida em que o valor de ∆ se aproxima do máximo, sobre a localização e as dimensões geométricas aproximadas da heterogeneidade (zona principal). / One of the major interpretation problems in geophysics is to determine the region in the subsurface which generates the main part of the signal. In this thesis, the position and size of this region, hereinafter called the main zone, have been found by modelling an electromagnetic system in which the source is an infinite line of electric current, extended over a conductive half-space. The earth has been modelled as a conductive half-space with an inhomogeneity in it as being an infinite layer or a prism of infinite length in the direction of the source line. The signal in the receiver of an electromagnetic system over a conductive homogeneous half-space is different from the one taken over the half space including an inhomogeneity. This difference is a function of the position of the inhomogeneity in relation to the transmiter-receiver system, besides other parameters. Therefore, with the other parameters fixed, there will be a specific position where this difference will maximize. Since this position depends on conductivity contrast, inhomogeneity dimensions and on source frequency, instead of a single position one will have a region where the inhomogeneity will give the maximum contribution to the measured signal. This region is called the main zone. Once the main zone is identified, the targets in the subsurface can be more precisely located. Usually they are conductive parts of the earth with some specific interest. One can facilitate the exploration and reduce production costs if these conductors are well identified during prospecting. A detectability function (∆) has been defined to measure the contribution to the signal due to the inhomogeneity. The ∆ function has been computed using amplitude and phase of the magnetic field components: Hx and Hz which are, respectively, the tangential and the normal to the earth's surface. The size and position of the main zone has been identified using the extremais of the ∆ function, which change with conductivity contrast, and the inhomogeneities' size and depth. Electromagnetic fields for one-dimensional models were calculated using a hybrid form, numerically solving the integrals that were obtained analytically. Two-dimensional models were computed numerically, by the finite elements technique. The maximum values of ∆ function, computed with amplitude of Hx, have been chosen to locate the main zone. This shows more stable results than other amplitude and phase components, both for one and two-dimensional models, when physical properties and geometric dimensions are changed. For the one-dimensional model, where the inhomogeneity is an infinitely extended horizontal layer, the depth of its central plane was found to be po = 0.17 δo, where po is the depth of this central plane and δo is the skin depth for the plane wave (in an homogeneous half-space having a conductivity σ1 equal to that of the backgound, and the frequency w corresponding to the maximum value of ∆ calculatede for the amplitude of Hx). For two-dimensional inhomogeneities, the co-ordinates of the main zone central axis was found to be do = 0,77 r0 (where do is the horizontal distance from this axis to the source) and po = 0,36 δo (where po is the depth of this central axis), with r0 being the source-receiver separation and δo the skin depth in the same conditions as in the one-dimensional case. If the values of r0 and δo are known, it is possible to determine (do, po). Associating each value of ∆ function (calculated using the amplitude of Hx) with the values of d = 0,77 r and p = 0,36 δ for each r and w used to generate ∆, a method to locate the main zone is sugested. The isovalue curves of ∆ are plotted to construct sections of ∆. These sections indicate the conductors position and provide some helpful insight into their geometric forms when the values of ∆ get dose to the maximum.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-06T12:01:00Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseSensibilidadeEstereotomografia.pdf: 6412555 bytes, checksum: 4d766be627a360d62d4651c10cc03c7a (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:03:42Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-10T14:12:16Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseSensibilidadeEstereotomografia.pdf: 6412555 bytes, checksum: 4d766be627a360d62d4651c10cc03c7a (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-18T11:31:24Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseSensibilidadeEstereotomografia.pdf: 6412555 bytes, checksum: 4d766be627a360d62d4651c10cc03c7a (MD5) / Made available in DSpace on 2014-09-18T11:31:24Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseSensibilidadeEstereotomografia.pdf: 6412555 bytes, checksum: 4d766be627a360d62d4651c10cc03c7a (MD5) Previous issue date: 2005 / ANP - Agência Nacional do Petróleo / O método de estereotomografia é estendido para meios com anisotropia arbitrária e implementado para meios com anisotropia elíptica e anelíptica. Os modelos elípticos e anelípticos apresentam somente três parâmetros. Isto faz com que eles sejam menos sensíveis a ambiguidade, causada pela cobertura limitada dos raios em experimentos sísmicos de superfície e VSP, do que modelos transversalmente isotrópicos ou ortorrômbicos. As correspondentes aproximações para superfície de vagarosidade limita a validade desta implementação para eventos qP com anisotropia suave. Experimentos numéricos mostram o potencial e as limitações da estereotomografia para estimar macro modelos de velocidade adequados para o imageamento na presença de anisotropia e a importância dos eventos de transmissão de experimentos VSP de multiplo afastamento modelo para o sucesso desta abordagem. / Stereotomography is extended to general anisotropic models and implemented for elliptical and anelliptical anisotropy. The elliptical and anelliptical models present only three parameters. This makes them less sensitive to the ambiguity due to limited coverage of surface seismic experiments than transversaly isotropic or orthorhombic models. The corresponding approximations of the slowness surface restrict the validity of the present approach to qP events and mild anisotropy. Numerical experiments show the potential and the limitations of stereotomography in estimating macro-velocity models suitable for imaging in the presence of anisotropy as well as the importance of transmission events from multiple-offset VSP experiments for the success of the approach.

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

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-05T11:44:12Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AdaptacaoMetodoBott.pdf: 2164475 bytes, checksum: baccb80d9f2f333abfebec4284c9394c (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:04:06Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-04T16:52:16Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AdaptacaoMetodoBott.pdf: 2164475 bytes, checksum: baccb80d9f2f333abfebec4284c9394c (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-22T11:21:30Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AdaptacaoMetodoBott.pdf: 2164475 bytes, checksum: baccb80d9f2f333abfebec4284c9394c (MD5) / Made available in DSpace on 2014-09-22T11:21:30Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AdaptacaoMetodoBott.pdf: 2164475 bytes, checksum: baccb80d9f2f333abfebec4284c9394c (MD5) Previous issue date: 2013 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Apresentamos um novo método para estimar o relevo do embasamento de bacias sedimentares através da extensão analítica da expressão da placa Bouguer para contraste de densidade entre o pacote sedimentar e o embasamento decrescendo monotonicamente com a profundidade de acordo com uma lei hiperbólica. O método requer ruído contido nos dados gravimétricos tenha desvio padrão inferior a 0,01 mGal. As estimativas do relevo do embasamento são obtidas nas posições espaciais localizadas diretamente abaixo de cada observação. A metodologia foi aplicada a dados sintéticos de bacias sedimentares simuladas apresentando relevo do embasamento suave. O método produziu relevo do embasamento estimado próximo do relevo simulado. O método foi também aplicado a três conjuntos de dados reais: Bacia do Recôncavo (Brasil), Graben do Büyük Menderes (Turquia) e Graben de San Jacinto (Estados Unidos). As soluções produzidas pelo método proposto apresentaram estimativas de profundidades máximas em acordo com as informações geológicas disponíveis: 6 km para Bacia do Recôncavo, 1,6 km para Vale do Büyük Menderes e 2,2 km para Graben de San Jacinto. / We present a new method to estimate the basement relief of a sedimentary basin through the analytic extension of the expression of the slab Bouguer for the case of a density contrast between the basement and the sediments decreasing monotonically with depth according to a hyperbolic law. The method requires that the noise in data have a standard deviation smaller than 0.01 mGal. The basement relief estimates are obtained at positions located directly below each observation. The methodology was applied on synthetic data from a simulated sedimentary basin presenting smooth basement relief. The method produced estimated basement topographies close to simulated basement topographies. The method was also applied to three sets of real Bouguer anomalies: Recôncavo Basin (Brazil), Büyük Menderes Graben (Turkey) and San Jacinto Graben (United States). The solutions produced by the proposed method attained maximum basement depth estimates in accordance with the available geological information: 6 km to Reconcavo Basin, 1.6 km to Büyük Menderes Valley, and 2.2 km to San Jacinto Graben.

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

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-18T11:54:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoContatoCamadas.pdf: 10720889 bytes, checksum: 3476f811330f5421d14e271a8e3c2428 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:18:59Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-12T12:35:41Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoContatoCamadas.pdf: 10720889 bytes, checksum: 3476f811330f5421d14e271a8e3c2428 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-22T13:15:32Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoContatoCamadas.pdf: 10720889 bytes, checksum: 3476f811330f5421d14e271a8e3c2428 (MD5) / Made available in DSpace on 2014-09-22T13:15:32Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoContatoCamadas.pdf: 10720889 bytes, checksum: 3476f811330f5421d14e271a8e3c2428 (MD5) Previous issue date: 1996 / PETROBRAS - Petróleo Brasileiro S.A. / Dois dos principais objetivos da interpretação petrofísica de perfis são a determinação dos limites entre as camadas geológicas e o contato entre fluidos. Para isto, o perfil de indução possui algumas importantes propriedades: É sensível ao tipo de fluido e a distribuição do mesmo no espaço poroso; e o seu registro pode ser modelado com precisão satisfatória como sendo uma convolução entre a condutividade da formação e a função resposta da ferramenta. A primeira propriedade assegura uma boa caracterização dos reservatórios e, ao mesmo tempo, evidencia os contatos entre fluidos, o que permite um zoneamento básico do perfil de poço. A segunda propriedade decorre da relação quasi-linear entre o perfil de indução e a condutividade da formação, o que torna possível o uso da teoria dos sistemas lineares e, particularmente, o desenho de filtros digitais adaptados à deconvolução do sinal original. A idéia neste trabalho é produzir um algoritmo capaz de identificar os contatos entre as camadas atravessadas pelo poço, a partir da condutividade aparente lida pelo perfil de indução. Para simplificar o problema, o modelo de formação assume uma distribuição plano-paralela de camadas homogêneas. Este modelo corresponde a um perfil retangular para condutividade da formação. Usando o perfil de entrada digitalizado, os pontos de inflexão são obtidos numericamente a partir dos extremos da primeira derivada. Isto gera uma primeira aproximação do perfil real da formação. Este perfil estimado é então convolvido com a função resposta da ferramenta gerando um perfil de condutividade aparente. Uma função custo de mínimos quadrados condicionada é definida em termos da diferença entre a condutividade aparente medida e a estimada. A minimização da função custo fornece a condutividade das camadas. O problema de otimização para encontrar o melhor perfil retangular para os dados de indução é linear nas amplitudes (condutividades das camadas), mas uma estimativa não linear para os contatos entre as camadas. Neste caso as amplitudes são estimadas de forma linear pelos mínimos quadrados mantendo-se fixos os contatos. Em um segundo passo mantem-se fixas as amplitudes e são calculadas pequenas mudanças nos limites entre as camadas usando uma aproximação linearizada. Este processo é interativo obtendo sucessivos refinamentos até que um critério de convergência seja satisfeito. O algoritmo é aplicado em dados sintéticos e reais demonstrando a robustez do método. / Two of the main objectives in petrophysical log interpretation are to determine geologic bed boundaries and fluid contacts. For these, the induction log has several important properties: it is sensitive to fluid type and distribution in the pore space, the induction measurement is accurately modeled as a convolution of formation conductivity and the tool response function. The first property assures good discrimination of the reserves and at the same time delineates oil-water contacts. This information permits a fundamental zoning of the well log. The second property follows from the quasi-linear relationship between the induction log and formation conductivity. Thus it is possible to use linear system theory, and in particular digital filter design to adaptively deconvolve the original signal. The idea in this work is to produce an algorithm capable of identifying contacts between beds traversed by a borehole, given the apparent conductivity read by an induction tool. To simplify the problem, the formation model is assumed to be a distribution of plane-parallel homogeneous beds. This model corresponds to a rectangular formation conductivity profile. Using the digitized input log, inflexion points are obtained numerically as extrema of first derivatives. This generates a first approximation of the real formation profile. This estimated profile is then convolved with the tool response function giving an estimated apparent conductivity log. A conditioned least-mean-square cost function is defined in terms of the difference between measured and estimated apparent conductivity. Minimizing the cost function yields the bed conductivities. The optimization problem of finding the best rectangular profile for induction data is linear for amplitudes (bed conductivities), but non-linear estimation for bed boundaries. In this case amplitudes are estimated by linear least-squares maintaining fixed contacts. A second pass maintains fixed amplitudes and computes small changes in bed boundaries using a linearized approximation. This processes is iterated to obtain successive refinement until a convergence criteria is satisfied. The algorithm is applied on synthetic and real data showing the robustness of the 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.