Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-21T13:47:09Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EstimativaParametrosElasticos.pdf: 4171380 bytes, checksum: fd321286a8d47fdb48182e29b3e07b48 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:11:27Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-08T16:40:47Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EstimativaParametrosElasticos.pdf: 4171380 bytes, checksum: fd321286a8d47fdb48182e29b3e07b48 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T12:55:33Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EstimativaParametrosElasticos.pdf: 4171380 bytes, checksum: fd321286a8d47fdb48182e29b3e07b48 (MD5) / Made available in DSpace on 2014-09-10T12:55:33Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EstimativaParametrosElasticos.pdf: 4171380 bytes, checksum: fd321286a8d47fdb48182e29b3e07b48 (MD5) Previous issue date: 2003 / FINEP - Financiadora de Estudos e Projetos / PETROBRAS - Petróleo Brasileiro S.A. / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / UFPA - Universidade Federal do Pará / CTPETRO - Plano Nacional de Ciência e Tecnologia do Setor Petróleo e Gás Natural / As medidas de amplitude, polarização e vagarosidade contem informações sobre o meio onde a propagação de onda ocorre. Esta tese investiga esses dados com objetivo de estimar as propriedades elásticas deste meio. Coeficientes de reflexão podem ser estimados das amplitudes dos dados e dependem de forma não linear dos contrastes dos parâmetros elásticos e do contraste de densidade entre os meios separados por uma interface. Quando o contraste de impedância é fraco, as aproximações lineares para a refletividade qP são mais convenientes para inversão da densidade e dos parâmetros elásticos usando as análises de amplitude versus ângulo de incidência (AVO) e amplitude versus a direção do plano de incidência (AVD). Escrevendo as equações de Zoepprittz de forma separada nos permite escrever uma solução destas equações em termos das matrizes de impedância e polarização. Usando esta solução são determinadas aproximações lineares para a refletividade da onda qP considerando fraco contraste de impedância, fraca anisotropia mas com classe de simetria de arbitrária. As linearizações são avaliadas para diferentes geometrias de aquisição e várias escolhas do meio de referência. Estas aproximações apresentam bom desempenho comparado com o valor exato do coeficiente de reflexão da onda qP e de suas ondas convertidas para incidências de até 30° e meios que obedecem à hipótese de fraca anisotropia. Um conjunto de fraturas orientado é representado efetivamente por um meio transversalmente isotrópico (TI), as aproximações lineares da refletividade da onda qP podem ser usadas para estimar a orientação de fratura. Partindo deste pressuposto este problema consiste em estimar a orientação do eixo de simetria a partir de dados de refletividade de onda qP. Este trabalho mostra que são necessários múltiplos azimutes e múltiplas incidências para se obter uma estimativa estável. Também é mostrado que apenas os coeficientes das ondas qS e qT são sensíveis ao mergulho da fratura. Foi investigada a estimativa da anisotropia local através de dados de VSP multiazimutal dos vetores de polarização e vagarosidade. Foram usadas medidas da componente vertical do vetor de vagarosidade e o vetor de polarização de ondas qP diretas e refletidas. O esquema de inversão é validado através de exemplos sintéticos considerando diferentes escolhas do vetor normal à frente de onda no meio de referência, meios de referências e geometria de aquisição. Esta análise mostra que somente um subgrupo dos parâmetros elástico pode ser estimado. Uma importante aplicação desta metodologia é o seu potencial para a determinação de classes de anisotropia. A aplicação desta metodologia aos dados do mar de Java mostra que os modelos isotrópicos e TIV são inadequados para o ajuste desses dados. / Amplitude, polarization and the slowness vector measurements carry information about the medium where wave propagation occurs. This thesis investigates these data aiming at the recovery of elastic properties in anisotropic media. Reflection coefficients can be estimated from amplitude data and depend nonlinearly on elastic and density contrasts across an interface. When the impedance contrast is weak, the linear approximations for the qP reflectivity are more convenient for inversion of density and elastic parameters using analysis of amplitude versus the angle of incidence (AVO) and amplitude versus the direction of the incidence plane (AVD). Partitioning the linear system defined by Zoepprittz equations allows one to write the solution of these equations in terms of impedance and polarization matrices. Using this solution, linear approximations for the qP reflectivity are derived for weak impedance contrasts and arbitrary symmetry classes of anisotropy. The linear approximations are evaluated for different acquisition geometries and choice of the reference medium. The approximations for the reflection coefficients of the reflected qP and the converted waves are in good agreement with the exact solution for incidence angles up to 30° for media that satisfy the weak impedance assumption. If a single oriented set of fractures is represented by a transversely isotropic effective medium, the linear approximations for qP reflectivity can be used to estimate the fractures orientation. Under these assumptions this problem is reframed as the estimation of the symmetry axis orientation from qP reflectivity data. This work shows the requirement of multiple components and multiple azimuthal data and quantifies the minimum amount of data for stable estimation. Also it is shown that the reflection coefficients of converted waves qS and qT only are sensitive to fractures dip. The inversion of polarization and slowness from multiazimutal VSP data are investigated for the estimation of local anisotropy. We use measurements of the vertical component of the slowness vector and the qP polarization data of direct and reflected waves. The inversion algorithm is validated in synthetic data sets for different choices of the wave front normal, reference medium and acquisition geometries. This analysis shows that only a subset of elastic parameters is recovered. An important application of this approach is its potential to determine the class of anisotropy. The application of this methodology to the Java Sea data set shows that isotropy and transversely isotropic models are inadequate to fit the data.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-12T13:42:26Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoEfeitosEletrojato.pdf: 19652715 bytes, checksum: 6ae97b0886be3a40c3f5c2510c0ee42a (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:23:18Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-27T16:03:30Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoEfeitosEletrojato.pdf: 19652715 bytes, checksum: 6ae97b0886be3a40c3f5c2510c0ee42a (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T16:37:09Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoEfeitosEletrojato.pdf: 19652715 bytes, checksum: 6ae97b0886be3a40c3f5c2510c0ee42a (MD5) / Made available in DSpace on 2014-09-10T16:37:09Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoEfeitosEletrojato.pdf: 19652715 bytes, checksum: 6ae97b0886be3a40c3f5c2510c0ee42a (MD5) Previous issue date: 1995 / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / PETROBRAS - Petróleo Brasileiro S.A. / Este trabalho trata da comparação dos efeitos do eletrojato equatorial nos dados magnetotelúricos bi e tridimensionais. Ele compõe-se fundamentalmente de duas partes: Na primeira, que consiste na comparação dos efeitos dos modelos bi e tridimensionais (linha de corrente e eletrojato gaussiano), determinamos as dimensões que deve apresentar uma estrutura tridimensional para que a resposta magnetotelúrica proveniente dessa estrutura, nas polarizações YX e XY, possa ser substituída pela resposta magnetotelúrica proveniente de uma estrutura bidimensional nos modos TE e TM respectivamente, para a resistividade aparente e fase. Os resultados mostram, para um embasamento condutivo ou resistivo, que é necessário aumentar a dimensão da estrutura tridimensional na direção principal acima de dezesseis vezes, em relação à dimensão inicial, para que a resposta magnetotelúrica proveniente dessa estrutura na polarização YX, possa ser substituída pela resposta magnetotelúrica proveniente de uma estrutura bidimensional no modo TE; no caso da polarização ser XY, essa substituição ocorrerá, no modo TM bidimensional, se aumentarmos a dimensão da estrutura tridimensional na direção principal acima de vinte e seis vezes em relação à dimensão inicial. Na segunda parte, que consiste na determinação dos efeitos do eletrojato equatorial Onwumechilliano, posicionamos a estrutura tridimensional a 10 km, 100 km, 300 km, 500 km, 1000 km e 1500 km, respectivamente, de distância do eletrojato. Dessa forma, determinamos a distância na qual os efeitos do eletrojato equatorial Onwumechilliano sobre as sondagens magnetotelúricas são praticamente desprezíveis. Verificamos que, posicionando a estrutura tridimensional a 1500 km de distância da fonte, para embasamento condutivo, o eletrojato não afeta as sondagens magnetotelúricas. Quando o embasamento é resistivo, o eletrojato afeta as sondagens magnetotelúricas, mesmo posicionando-se a estrutura tridimensional a 1500 km de distância da fonte. De um modo geral, constatamos que o eletrojato equatorial afeta as sondagens magnetotelúricas para freqüências inferiores a 10-2 Hz, tanto na resistividade aparente quanto na fase e em ambas as polarizações, principalmente, quando o embasamento é resistivo. Observamos também que o efeito galvânico é mais pronunciado na polarização YX que na XY. / This work deals with the comparison between the effects of the equatorial electrojet in two and three-dimensional magnetotelluric data. It is composed of two parts: in the first part, which consists in comparing the effects of two-dimensional with three-dimensional models (current line and gaussian electrojet), we determine the dimensions that a three-dimensional structure should have so that the MT response of that structure, in YX and XY polarizations, can be replaced by the MT response of a 2D structure in TE and TM modes, for apparent resistivity and phase. The results show that the dimension of the 3D structure in the main direction must be over 16 times its length in the initial dimension, so that the MT response of that structure in YX polarization can be replaced by the MT response at a 2D structure in TE mode. If the polarization is XY the main dimension must be 26 times the initial dimension. In the second part, which consists in determining the effects of the Onwumechillian equatorial electrojet, we placed the 3D structure at 10 km, 100 km, 300 km, 500 km, 1000 km and 1500 km, from the electrojet. In this way we determine the distance at which the effects of the Onwumechillian equatorial electrojet on the MT soundings are practically negligible. We observe that, for a conductive basement, the equatorial electrojet does not affect the MT soundings if the 3D structure is placed at 1500 km from the source. When the basement is resistive, the electrojet affects the MT soundings even if the structure is at 1500 km away from the source. We verified that the equatorial electrojet affects the MT soundings in frequencies under 10-2Hz, in both resistivity and phase and in both polarizations. We also observed that the galvanic effect is stronger on the YX than on the XY polarization.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-11T12:34:33Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoEfeitoJanela.pdf: 12983544 bytes, checksum: d8af62d9b160f107ec5594c633121e76 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:24:38Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-04T13:08:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoEfeitoJanela.pdf: 12983544 bytes, checksum: d8af62d9b160f107ec5594c633121e76 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T16:37:25Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoEfeitoJanela.pdf: 12983544 bytes, checksum: d8af62d9b160f107ec5594c633121e76 (MD5) / Made available in DSpace on 2014-09-10T16:37:25Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoEfeitoJanela.pdf: 12983544 bytes, checksum: d8af62d9b160f107ec5594c633121e76 (MD5) Previous issue date: 1999 / A presente dissertação consta de estudos sobre deconvolução sísmica, onde buscamos otimizar desempenhos na operação de suavização, na resolução da estimativa da distribuição dos coeficientes de reflexão e na recuperação do pulso-fonte. Os filtros estudados são monocanais, e as formulações consideram o sismograma como o resultado de um processo estocástico estacionário, e onde demonstramos os efeitos de janelas e de descoloração. O principio aplicado é o da minimização da variância dos desvios entre o valor obtido e o desejado, resultando no sistema de equações normais Wiener-Hopf cuja solução é o vetor dos coeficientes do filtro para ser aplicado numa convolução. O filtro de deconvolução ao impulso é desenhado considerando a distribuição dos coeficientes de reflexão como uma série branca. O operador comprime bem os eventos sísmicos a impulsos, e o seu inverso é uma boa aproximação do pulso-fonte. O janelamento e a descoloração melhoram o resultado deste filtro. O filtro de deconvolução aos impulsos é desenhado utilizando a distribuição dos coeficientes de reflexão. As propriedades estatísticas da distribuição dos coeficientes de reflexão tem efeito no operador e em seu desempenho. Janela na autocorrelação degrada a saída, e a melhora é obtida quando ela é aplicada no operador deconvolucional. A transformada de Hilbert não segue o princípio dos mínimos-quadrados, e produz bons resultados na recuperação do pulso-fonte sob a premissa de fase-mínima. O inverso do pulso-fonte recuperado comprime bem os eventos sísmicos a impulsos. Quando o traço contém ruído aditivo, os resultados obtidos com auxilio da transformada de Hilbert são melhores do que os obtidos com o filtro de deconvolução ao impulso. O filtro de suavização suprime ruído presente no traço sísmico em função da magnitude do parâmetro de descoloração utilizado. A utilização dos traços suavizados melhora o desempenho da deconvolução ao impulso. A descoloração dupla gera melhores resultados do que a descoloração simples. O filtro casado é obtido através da maximização de uma função sinal/ruído. Os resultados obtidos na estimativa da distribuição dos coeficientes de reflexão com o filtro casado possuem melhor resolução do que o filtro de suavização. / The present masters dissertation consists of studies on seismic deconvolution where we look for otimizing the operations of smoothing, of resolution for the estimation of distribution of reflection coefficients, and of recovery of source-pulse. The studied filters are single channel, and the formulations consider the seismic trace as the result of a stationary stochastic process, and we demonstrate the effects of taper windows and of prewhitening in resolution. The applied principle is the minimization of the difference's variance between real and desired outputs, resulting on a system Wiener-Hopf normal equations whose solution is the vector of filter coefficients to be applied in a convolution. Spike deconvolution is designed considering the distribution of reflection coefficients as a white series. The operator compresses the seismic events to impulses, and its inverse is a good approximation to the source-pulse. The application of taper windows and of prewhitening improve the output of this filter. Spike-series filters are designed using the distribution of reflection coefficients. The statistical properties of the reflection coefficients distribution affect the operator and its performance. Taper windows on autocorrelation degrades the output, and improvement is achieved when it is applied to the deconvolutional operator. The Hilbert transform produces good results in the recovery of source-pulse, under the premise of minimum-phase. The inverse of the recovered source-pulse compresses the seismic events to impulses. When the seismic trace contains additive noise, the results obtained with Hilbert transform are better than with spike deconvolution. Smoothing filter suppresses noise in the seismic trace as a function of a prewhitening parameter. The use of smoothed traces improves the spike deconvolution. Double prewhitening generates better results than the simple prewhitening. The matched filter operator is obtained from maximization of a signal/noise ratio function. Deconvolving the output of the matched filter for the estimation of the distribution of reflection coefficients possess better resolution than when using a smoothing filter.

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.