Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-18T16:11:18Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DispersaoOndasRayleigh.pdf: 20019750 bytes, checksum: afe92614b0dfea249cf425b1afcd0e2f (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:18:03Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-02T16:35:34Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DispersaoOndasRayleigh.pdf: 20019750 bytes, checksum: afe92614b0dfea249cf425b1afcd0e2f (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-04T16:43:52Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DispersaoOndasRayleigh.pdf: 20019750 bytes, checksum: afe92614b0dfea249cf425b1afcd0e2f (MD5) / Made available in DSpace on 2014-09-04T16:43:52Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DispersaoOndasRayleigh.pdf: 20019750 bytes, checksum: afe92614b0dfea249cf425b1afcd0e2f (MD5) Previous issue date: 1989 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Este trabalho representa um estudo de dispersão da componente vertical da onda de superfície de Rayleigh com trajetórias na plataforma Sulamericana. Os registros utilizados são provenientes das estações localizadas no território Brasileiro; sendo a do Rio de Janeiro (RDJ), a de Brasília (BDF), a de Caicó (CAI) e a de Belém (BEB), pois estas são as únicas estações sismológicas no Brasil que têm sensores de período longo e que servem para o estudo de dispersão no intervalo de 4 a 50 segundos, aqui realizados. Os terremotos utilizados estão localizados ao longo da parte leste da cadeia Andina e dentro da plataforma Sulamericana com trajetórias tipicamente continental. Foram selecionados 34 eventos com a utilização dos seguintes critérios práticos: a localização, a magnitude m_b e a profundidade, ocorridos durante o período de Janeiro de 1978 até Junho de 1987. O estudo de dispersão aqui abordado significa a determinação da velocidade de grupo e das amplitudes espectrais correspondentes aos harmônicos fundamental e primeiro superior. Normalmente os harmônicos de ordem segunda ou maior são raramente disponíveis na observação. Dois tipos de medidas foram feitas: (i) velocidade de grupo vs. período e (ii) amplitude vs. período. Os estudos de dispersão são fundamentais para determinação da estrutura da crosta e manto superior que estão diretamente relacionados com os fenômenos geológicos. Neste trabalho, regionalização é definida como a identificação das diferentes formas de curvas de dispersão, que estão relacionadas com as trajetórias epicentro-estação ao longo da plataforma Sulamericana e que venham ter uma correlação geológica como está descrito no item 4.3 deste trabalho. A distribuição dos epicentros se faz desde o extremo sul da Argentina até o extremo norte da Venezuela, objetivando iniciar com este trabalho uma sistemática voltada aos estudos de regionalização da plataforma Sulamericana na nossa instituição. Neste trabalho foram observados três tipos distintos de curvas em 27 trajetórias e agrupadas por famílias 1,2 e 3 respectivamente, onde procurou-se correlacionar suas diferentes formas com a geologia regional da plataforma Sulamericana. A obtenção da curva de dispersão foi feita através da técnica do filtro múltiplo (Dziewonski et al, 1969). Este filtro tem a propriedade de separar os harmônicos através das suas velocidades de grupo para cada frequência selecionada, e também de recuperar as amplitudes características dos harmônicos (Herrmann, 1973). O desenvolvimento teórico do filtro bem como suas limitações e forma de uso são tratados por Dziewonski et al (1972). Como parte do trabalho há a implantação, adaptações e o desenvolvimento de parte do fluxograma do filtro múltiplo, bem como a estruturação da digitalização dos dados para o processamento e interpretação não-automática dos resultados do processamento. / This work represents a dispersion study of the vertical component of the Rayleigh wave with trajectories on the South American platform. The records used are from stations located in the Brazilian territory such as the ones from Rio de Janeiro (RDJ), Caicó (CAI), Brasília (BDF) and Belém (BEB). These are the only seismological stations in Brazil with long period sensors that can be useful for the dispersion studies in the interval 4 to 50 seconds. The earthquakes used are localized along the eastern part of Andean mountain belt and inside the South American platform, with tipically continental trajectories. 34 events were selected using the following pratical criteria: the localization, the magnitude m_b, the depth, and that they occured in the period of January 1978 to June 1987. The dispersion study done here means the determination of group velocities and spectral amplitudes, corresponding to the fundamental and first higher modes. Normally the modes of second order or more are rarely observed. Two types of measurements were made: (i) group velocity versus period and (ii) spectral amplitude versus period. Dispersion studies are important for determining the structure of the crust and upper mantle which are directly related to gelogical phenomena. In this work, regionalization is defined as the identification of different forms of dispersion curves and their relation to the epicenter-station trajectory on the South American platform. We also correlate them geologicaly, as described in item 4.3 of this work. The distribution of epicenters occurs from the extreme South of Argentina to the extreme North of Venezuela, with the objective of iniciating with this work a sistematic study of the sub-Andean area in our institution. Three distinct types of curves were observed in 27 trajectories and grouped in three families (1, 2 and 3). Their forms were compared with the regional geological features of South American platform. The multiple filter technique was used to obtain the dispersion curves (Dziewonski et al 1969). This filter has the property of separating the modes through their group velocities for each frequency selected and also recovers the espectral amplitudes characteristic for each harmonic (Herrmann, 1973). The theorectical development of the filter properties as well as the limitations and aplicabilities are treated by Dziewonski et al (1972). As part of this work there was the implementation, adaptions and the layout of fluxograms of the multiple filter software, as well as the digitalization procedure for data processing and non-automatic interpretation of the results.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-07-09T12:29:50Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoEstabilidadeMecanica.pdf: 11477602 bytes, checksum: 6af6cc1a1e5c0885abc416d0e5bdc521 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:34:20Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-12T15:00:46Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoEstabilidadeMecanica.pdf: 11477602 bytes, checksum: 6af6cc1a1e5c0885abc416d0e5bdc521 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-19T13:25:48Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoEstabilidadeMecanica.pdf: 11477602 bytes, checksum: 6af6cc1a1e5c0885abc416d0e5bdc521 (MD5) / Made available in DSpace on 2014-09-19T13:25:48Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoEstabilidadeMecanica.pdf: 11477602 bytes, checksum: 6af6cc1a1e5c0885abc416d0e5bdc521 (MD5) Previous issue date: 1992 / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / FINEP - Financiadora de Estudos e Projetos / A análise da estabilidade mecânica de um poço pode ser feita a partir do cálculo de parâmetros elásticos da formação utilizando a densidade do meio e as velocidades de propagação das ondas compressional e cisalhante na formação rochosa, os quais podem ser obtidos de perfis geofísicos do poço. Em formações sedimentares pouco consolidadas as ferramentas de perfilagem sônica convencionais (monopolares) não conseguem registrar com acuidade a velocidade da onda cisalhante pois a primeira chegada dessa onda é camuflada pela chegada de outras ondas que podem ser mais rápidas que a onda cisalhante num poço perfurado neste tipo de formação. Medidas das velocidades sônicas são feitas em laboratório em amostras da formação, sob condições semelhantes às condições in situ, servindo como ajuste das velocidades registradas no poço pela ferramenta de perfilagem sônica. Para a análise de estabilidade da formação, perfis auxiliares são necessários como o perfil de porosidade, saturação de fluidos e perfis de composição mineralógica da formação rochosa. Exige-se ainda dados de testes de avaliação da formação e de condições do reservatório, mas que são comuns em poços de petróleo, como o teste de formação e os testes de pressurização do poço, tais como o teste de micro-fraturamento hidráulico ou o teste de absorção. A avaliação das tensões principais efetivas que atuam distante do poço, e que não são afetadas pela sua presença, é feita através da associação de um modelo de deformação elástica apropriado e o resultado do teste de pressurização disponível para o poço em estudo. Utilizando resultados clássicos da teoria da elasticidade geral pode-se calcular o campo de tensões modificado na vizinhança da parede do poço devido ao efeito da própria presença do poço ali perfurado e da diferença de pressão existente entre o interior do poço e a formação rochosa. A determinação das propriedades mecânicas da formação a partir das velocidades sônicas e a avaliação do campo de tensões assumindo um modelo elástico de deformação, supõem o meio rochoso no qual as ondas se propagam como elástico, homogêneo e isotrópico. Esta suposição representa a principal aproximação assumida pela metodologia descrita neste trabalho. De posse das propriedades mecânicas da formação e do campo de tensões que age na vizinhança do poço resta definir o critério segundo o qual a rocha sofre instabilidade mecânica quando submetida aquele campo de tensões. Isto permite determinar se, nas condições avaliadas do poço e da formação, haverá quebra da parede do poço por excesso de tensão e, se houver, qual a sua extensão. Assim o problema é como analisar o comportamento mecânico de um poço em uma formação pouco consolidada a partir de perfis geofísicos os quais podem ter problemas no registro das propriedades físicas do meio em formações deste tipo. A metodologia proposta é aplicada a dois intervalos de profundidade pertencentes a dois poços onde arenitos e folhelhos se intercalam e nos quais todos os dados necessários estão disponíveis. Os resultados obtidos mostram que, exceto quando outros mecanismos de quebra da parede do poço agem na formação, a metodologia proposta consegue com sucesso detectar zonas de ocorrência de instabilidade mecânica do poço provocadas por um campo de tensões que excede a resistência mecânica da formação. / The stability of boreholes can be evaluated from the rocks' mechanical properties, which are compared to the acting stress field using an appropriate failure criterion. The elastic parameters can be measured on rock samples in the laboratory (static parameters) or can be calculated from acoustic velocities and the formation bulk density (dynamic parameters). Dynamic elastic parameters can be obtained from wireline logs, but in poorly consolidated rocks the arrival of the shear wave is often masked by the Stoneley wave, which may travel faster and has a higher amplitude. The giant Campos basin offshore Brazil contains such rocks. We performed laboratory measurements of acoustic velocities on samples from this basin and we use them to correct the shear velocities obtained from waveform processing. Based on this we calculate elastic properties (e.g. Poisson's ratio, bulk modulus) of the rocks which we assume to be homogeneous and isotropic. Mechanical properties such as the shear, compressional and tensional strength are obtained from empirical relations with mineralogy and the elastic parameters. The stress field was evaluated using an uniaxial strain model and actual formation pressure data. From this, near-wellbore stresses were calculated using the modified Kirsch equations. The stress concentrations around the wellbore are illustrated for a number of representative cases. Wellbore failure occurs when these local stresses exceed the formation strength. Three different failure criteria were tested: The Mohr-Coulomb, the Griffith and the Fairhurst criterion. They respectively represent failure due to shear, tension and a combination of compression and tension. Actual wellbore instability was determined from the differential caliper measurement. Two sequences were analyzed, both consisting of intercalations of sands and shales. We found that the Fairhurst criterion offered the best match with the actually observed wellbore failure, and we propose it as a method to predict borehole instability for this basin.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-07-07T13:55:27Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EquipamentoEletromagneticoProspeccao.pdf: 10600477 bytes, checksum: 4318a450123d026ada0e49bfca7dbaf6 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Item sem assunto. on 2014-08-07T16:26:28Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-08T14:11:50Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EquipamentoEletromagneticoProspeccao.pdf: 10600477 bytes, checksum: 4318a450123d026ada0e49bfca7dbaf6 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-22T12:37:55Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EquipamentoEletromagneticoProspeccao.pdf: 10600477 bytes, checksum: 4318a450123d026ada0e49bfca7dbaf6 (MD5) / Made available in DSpace on 2014-09-22T12:37:55Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EquipamentoEletromagneticoProspeccao.pdf: 10600477 bytes, checksum: 4318a450123d026ada0e49bfca7dbaf6 (MD5) Previous issue date: 1979 / Nas últimas décadas, a prospecção por métodos eletromagnéticos vem-se constituindo numa técnica eficiente para prospecção mineral. O objetivo deste trabalho foi desenvolver um equipamento para prospecção eletromagnética quantitativa de corpos condutores, através do método dipolo-dipolo, podendo ainda ser usado em modelos reduzidos. Eletricamente, o sistema mede grandezas relacionadas ao acoplamento indutivo entre duas bobinas: transmissora e receptora. Elas são dispostas na superfície da terra, afastadas entre si, e a terra, desse modo, constitui o núcleo acoplador. Quando existem corpos condutores nas proximidades, estes são denunciados por alterações no comportamento do sinal induzido na bobina receptora. O equipamento compreende dois conjuntos: o transmissor e o receptor, além de acessórios. O transmissor gera um campo eletromagnético nas freqüências de 520 e 3.090 Hz, e um sinal de referência para o receptor, o qual é enviado através de um cabo. O receptor, inicialmente, separa o sinal induzido pelos campos secundários gerados por condutores, do campo normalmente recebido, quando a condutividade da subsuperfície é relativamente uniforme (campo primário). Em seguida, decompõe esse sinal em duas componentes ortogonais, uma em fase, e outra em quadratura com o campo primário. Através de duas escalas de precisão, as amplitudes dessas componentes são mostradas como percentagens do campo primário, com precisão de 1%. A sensibilidade do receptor é de 0,5 μV. O circuito eletrônico foi rigorosamente testado com preciso instrumental de laboratório. Em seguida, testou-se sua aplicação no Laboratório de Modelo Reduzido Eletromagnético do NCGG, refazendo-se experiências clássicas, encontradas na literatura especializada. No campo, foi experimentado próximo da cidade de Araci no Estado da Bahia, em áreas prospectadas pela "Rio Doce Geologia e Mineração S/A-DOCEGEO". Em ambos os casos, verificaram-se bons resultados. / In the last decade, electromagnetic methods constituted an efficient technique for prospecting of minerals. The objective of this work was to develop an equipment for quantitative electromagnetic prospecting of conducting bodies using dipole-dipole technique. The equipment can also be adopted for reduced scale model experiments. The equipment, in principle, measures the variation in the mutual electromagnetic coupling of two coils, namely the transmitting and the receiving coils. Ground will also contribute to this coupling when the coil system is used near it. This mutual impedance of the transmitter – receiver - ground system shall vary in the vicinity of an anomalous conductor which is manifested as a variation in the signal induced in the search coil. The equipment consists of two sets: a transmitter and a receiver alongwith accessories. The transmitter generates an electromagnetic field at frequencies of 520 Hz and 3090 Hz, and a reference signal for the receiver which is connected through a cable. The receiver first separates the anomalous secondary e.m. field from the uniform primary field. Then this secondary signal is decomposed into two orthogonal components, one inphase and other in quadrature with respect to the primary field. Amplitudes of these phase components can be measured directly in terms of the percentage of the primary field within a precision of 1%, using two amplitude scales. The sensibility of the receiver is 0.5 μV. Electronic circuits were rigorously tested by using precise laboratory instruments. The performance of the e.m. equipment was verified (a) by carrying out classical e.m. model experiments in the laboratory and comparing them with the published results; and (b) by taking some e. m. traverses in an area near Araci, in the State of Bahia, which has been prospected by "Rio Doce Geologia e Mineraçáo S/A-DOCEGEO". All these tests verify a highly good performance of the equipment in terms of its precision and usefulness for field work.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-04T13:28:01Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_FormulacaoImplementacaoPacote.pdf: 18521004 bytes, checksum: 5a1e56c851024f3a4b89fa031123237e (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Item sem assunto. on 2014-08-07T16:42:20Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-04T12:41:23Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_FormulacaoImplementacaoPacote.pdf: 18521004 bytes, checksum: 5a1e56c851024f3a4b89fa031123237e (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T14:50:06Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_FormulacaoImplementacaoPacote.pdf: 18521004 bytes, checksum: 5a1e56c851024f3a4b89fa031123237e (MD5) / Made available in DSpace on 2014-09-10T14:50:06Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_FormulacaoImplementacaoPacote.pdf: 18521004 bytes, checksum: 5a1e56c851024f3a4b89fa031123237e (MD5) Previous issue date: 1992 / PETROBRAS - Petróleo Brasileiro S.A. / UFPA - Universidade Federal do Pará / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / FINEP - Financiadora de Estudos e Projetos / O presente trabalho consiste na formulação de uma metodologia para interpretação automática de dados de campo magnético. Desta forma, a sua utilização tornará possível a determinação das fronteiras e magnetização de cada corpo. Na base desta metodologia foram utilizadas as características de variações abruptas de magnetização dos corpos. Estas variações laterais abruptas serão representadas por polinômios descontínuos conhecidos como polinômios de Walsh. Neste trabalho, muitos conceitos novos foram desenvolvidos na aplicação dos polinômios de Walsh para resolver problemas de inversão de dados aeromagnéticos. Dentre os novos aspectos considerados, podemos citar. (i) - O desenvolvimento de um algoritmo ótimo para gerar um jôgo dos polinômios "quase-ortogonais" baseados na distribuição de magnetização de Walsh. (ii) - O uso da metodologia damped least squares para estabilizar a solução inversa. (iii) - Uma investigação dos problemas da não-invariância, inerentes quando se usa os polinômios de Walsh. (iv) - Uma investigação da escolha da ordem dos polinômios, tomando-se em conta as limitações de resolução e o comportamento dos autovalores. Utilizando estas características dos corpos magnetizados é possível formular o problema direto, ou seja, a magnetização dos corpos obedece a distribuição de Walsh. É também possível formular o problema inverso, na qual a magnetização geradora do campo observado obedece a série de Walsh. Antes da utilização do método é necessária uma primeira estimativa da localização das fontes magnéticas. Foi escolhida uma metodologia desenvolvida por LOURES (1991), que tem como base a equação homogênea de Euler e cujas exigências necessárias à sua utilização é o conhecimento do campo magnético e suas derivadas. Para testar a metodologia com dados reais foi escolhida uma região localizada na bacia sedimentar do Alto Amazonas. Os dados foram obtidos a partir do levantamento aeromagnético realizado pela PETROBRÁS. / In this work we develop a method for the automatic interpretation of magnetic field data. This methodology determines the boundaries and magnetizations of geologic units. The often of abrupt variations of the magnetization of geological units are modeled in this methodology. These abrupt lateral variations will be represented by discontinuous polynomials known as Walsh polynomials. Several new concepts were developed in order to apply Walsh polynomials to solve aeromagnetic data inversion problems. Among these new concepts we can list. (i) - An optimum algorithm to generate a set of "quasi-orthogonal" polynomials using a Walsh magnetization distribution. (ii) - The use of the damped least squares method to stabilize the inverse solution. (iii) An investigation of positional invariance problems, which are inherent in the Walsh polynomials. (iv) - The choice of the order of the polynomials taking into account the limitations of resolution and the behaviour of the eigenvalues. Using these features of magnetized bodies it is possible to formulate the direct problem, i.e, the bodies magnetization obeys Walsh distribution. It is also possible to formulate the inverse problem, in which the magnetization that generates the observed field obeys the Walsh series. Prior to applying the method it is necessary to estimate the depths to the magnetic sources. We have chosen a method developed by LOURES (1991), which is based on the Euler homogenous equation which requires the knowledge of the magnetic field and its derivatives. We tested the method with real data from a region in the Alto Amazonas sedimentary basin. The aeromagnetic data were provided by PETROBRÁS.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-24T12:47:44Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoQuadraticaDados.pdf: 975991 bytes, checksum: 5c694f4239ad08ee42219b354ed909b9 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-30T13:41:17Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoQuadraticaDados.pdf: 975991 bytes, checksum: 5c694f4239ad08ee42219b354ed909b9 (MD5) / Made available in DSpace on 2014-09-30T13:41:17Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoQuadraticaDados.pdf: 975991 bytes, checksum: 5c694f4239ad08ee42219b354ed909b9 (MD5) Previous issue date: 2004 / Neste trabalho são apresentadas expressões exatas e aproximações quadráticas das conhecidas equações de Knott-Zöeppritz, que calculam as amplitudes dos coeficientes Rpp e Rsp em termos dos contrastes médios relativos, bem como é procedida a inversão dos parâmetros de dados de AVO a partir destas aproximações quadráticas. Nesta inversão é utilizado o algoritmo de Levenberg-Marquardt, e são considerados apenas os eventos refletidos Rpp e convertidos Rsp, não associados. Nos estudos dos parâmetros físicos dos meios contrastes de impedância (δz), módulo de cisalhamento (δμ) e velocidade da onda p (δα), verificou-se quais desses parâmetros podem ser invertidos. Os resultados obtidos mostram que o contraste de impedância (δz) é muito bem resolvido estando ele relacionado com o contraste de velocidade da onda p (δα) ou com o contraste de cisalhamento (δμ), no caso de eventos refletidos considerando modelos de alto, moderado e baixo contrastes. Por outro lado ao se fixar o contraste de impedância e relacionar os outros dois parâmetros em consideração, os resultados mostram que esses são mal resolvidos, ou seja, a região de ambiguidade torna-se muito grande e os parâmetros tornam-se ambíguos e instáveis. No caso do evento convertido e na combinação do evento refletido com o convertido, para os modelos de baixo e moderado contrastes, (δz) é muito bem resolvido, caso que não acontece para modelo de alto contraste. Diante desses resultados verifica-se que no procedimento de inversão quadrática de dados de AVO, fixado (δμ), a recuperação dos dois parâmetros variados é muito boa, no caso do evento refletido, e razoavelmente boa no caso do evento convertido, por esse motivo optou-se pela fixação do módulo de cisalhamento. / In this work, the exact and quadratic approach expression of Knott-Zöeprtiz are presented, as well as its coefficients, amplitudes, Rpp and Rsp, respectively, all of them as a function of the relative contrast, as well as the inversion of quadratic parameters of AVO data. In this inversion of the quadratic parameters, we have used the non-linear Levenberg-Marquadt algorithm, considering only converted and reflected events, either in the non-combined form. In the study of the physical parameters impedance contrast, shear modulus and P wave velocity, we have werified which of those parameters can be effectively inverted. The results have shown that the impedance contrast (δz) is very well solved when it is related with the P wave velocity contrast (δα) or with the shear contrast (δμ), considering reflected events in the case of high, moderate and low model contrasts. On the other hand, when the impedance constrast is fixed and the two other parameters in consideration are related, the results show that these parameters are bad solved, e.g., the ambiguity region becomes large and the parameters become ambiguous and unstable. In the case of a converted event and in the combination of converted and reflected events, considering low and moderate contrast models (δz) is well solved, which is not the case when we consider a high contrast model.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-10-07T12:27:18Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MedidasCoerenciaAnalise.pdf: 2796169 bytes, checksum: 3d1e7250bfa49952d03455487988a69a (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-10-09T13:15:11Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MedidasCoerenciaAnalise.pdf: 2796169 bytes, checksum: 3d1e7250bfa49952d03455487988a69a (MD5) / Made available in DSpace on 2014-10-09T13:15:11Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MedidasCoerenciaAnalise.pdf: 2796169 bytes, checksum: 3d1e7250bfa49952d03455487988a69a (MD5) Previous issue date: 2011 / ANP - Agência Nacional do Petróleo / PETROBRAS - Petróleo Brasileiro S.A. / FINEP - Financiadora de Estudos e Projetos / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Métodos iterativos para a análise de velocidade com migração dependem das funções objetivo medirem a horizontalização dos eventos de reflexão em famílias de ponto imagem comum (CIG). A migração em tempo é um método simples de imageamento para avaliar essas funções objetivo. Usando migração em tempo, estudamos a influência da função objetivo nos resultados da análise de velocidade com migração. Propomos duas novas funções objetivo para a análise velocidade com migração: Semblance Diferencial Estendido e o produto do Semblance Clássico com o Semblance Diferencial Estendido. Experimentos numéricos, utilizando os dados Marmousoft, mostram a eficácia das novas funções objetivo em estimar modelos velocidade produzindo eventos horizontalmente alinhados nas famílias de ponto imagem comum. / Iterative methods for migration velocity analysis depend on objective functions to measure the flatness of reflection events in common image gathers (CIG). Time migration is a simple imaging method to evaluate these objective functions. Using time migration we studied the influence of objective functions on the results of migration velocity analysis. We propose two new objective functions for migration velocity analysis: Extended Differential Semblance and the product of Classical Semblance times the Extended Differential Semblance. Numerical experiments using the Marmousoft data show the effectiveness of the new objective functions to estimate velocity models producing at events in common image gathers.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-10-20T13:51:09Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MetodosOtimizacaoGlobal.pdf: 2631790 bytes, checksum: d0343a988ffbcc231d9547a1fa5255c5 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-10-22T13:40:43Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MetodosOtimizacaoGlobal.pdf: 2631790 bytes, checksum: d0343a988ffbcc231d9547a1fa5255c5 (MD5) / Made available in DSpace on 2014-10-22T13:40:43Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MetodosOtimizacaoGlobal.pdf: 2631790 bytes, checksum: d0343a988ffbcc231d9547a1fa5255c5 (MD5) Previous issue date: 2004 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O método de empilhamento por Superfície de Reflexão Comum (SRC) produz seções simuladas de afastamento nulo (AN) por meio do somatório de eventos sísmicos dos dados de cobertura múltipla contidos nas superfícies de empilhamento. Este método não depende do modelo de velocidade do meio, apenas requer o conhecimento a priori da velocidade próxima a superfície. A simulação de seções AN por este método de empilhamento utiliza uma aproximação hiperbólica de segunda ordem do tempo de trânsito de raios paraxiais para definir a superfície de empilhamento ou operador de empilhamento SRC. Para meios 2D este operador depende de três atributos cinemáticos de duas ondas hipotéticas (ondas PIN e N), observados no ponto de emergência do raio central com incidência normal, que são: o ângulo de emergência do raio central com fonte-receptor nulo (β0) , o raio de curvatura da onda ponto de incidência normal (RPIN) e o raio de curvatura da onda normal (RN). Portanto, o problema de otimização no método SRC consiste na determinação, a partir dos dados sísmicos, dos três parâmetros (β0, RPIN, RN) ótimos associados a cada ponto de amostragem da seção AN a ser simulada. A determinação simultânea destes parâmetros pode ser realizada por meio de processos de busca global (ou otimização global) multidimensional, utilizando como função objetivo algum critério de coerência. O problema de otimização no método SRC é muito importante para o bom desempenho no que diz respeito a qualidade dos resultados e principalmente ao custo computacional, comparado com os métodos tradicionalmente utilizados na indústria sísmica. Existem várias estratégias de busca para determinar estes parâmetros baseados em buscas sistemáticas e usando algoritmos de otimização, podendo estimar apenas um parâmetro de cada vez, ou dois ou os três parâmetros simultaneamente. Levando em conta a estratégia de busca por meio da aplicação de otimização global, estes três parâmetros podem ser estimados através de dois procedimentos: no primeiro caso os três parâmetros podem ser estimados simultaneamente e no segundo caso inicialmente podem ser determinados simultaneamente dois parâmetros (β0, RPIN) e posteriormente o terceiro parâmetro (RN) usando os valores dos dois parâmetros já conhecidos. Neste trabalho apresenta-se a aplicação e comparação de quatro algoritmos de otimização global para encontrar os parâmetros SRC ótimos, estes são: Simulated Annealing (SA), Very Fast Simulated Annealing (VFSA), Differential Evolution (DE) e Controlled Rando Search - 2 (CRS2). Como resultados importantes são apresentados a aplicação de cada método de otimização e a comparação entre os métodos quanto a eficácia, eficiência e confiabilidade para determinar os melhores parâmetros SRC. Posteriormente, aplicando as estratégias de busca global para a determinação destes parâmetros, por meio do método de otimização VFSA que teve o melhor desempenho foi realizado o empilhamento SRC a partir dos dados Marmousi, isto é, foi realizado um empilhamento SRC usando dois parâmetros (β0, RPIN) estimados por busca global e outro empilhamento SRC usando os três parâmetros (β0, RPIN, RN) também estimados por busca global. / The Common Reflection Surface Stack (CRS) method provides the simulation of zerooffset (ZO) sections by means of the summing seismic events of the multicoverage data contained in the stacking surface. This method does not depend on the velocity macro-model of medium; it only requires a priori knowledge of the near-surface velocity. The simulation of ZO sections for this stacking method uses a hyperbolic second-order traveltime approximation of the paraxial rays to define the CRS stacking surface or CRS stack operator. For 2D media, this operator depends on three kinematic attributes of two hypothetical waves (NIP and N waves), observed in the point of emergency of the central ray with normal incidence, namely, the angle of emergency of the central ZO ray (β0), the radius of curvature of the Normal Incidence Point Wave (RNIP) and the radius of curvature of the Normal Wave (RN). Therefore, the optimization problem in the CRS method consists in the determination, from the seismic data, of the three optimal parameters (β0, RNIP, RN) associated to each sample point of ZO section to be simulated. The simultaneous determination of these parameters can be made by means of multidimensional global search process (or global optimization), using as objective function some coherence criterion. The optimization problem in CRS method is very important for the good performance with respect to quality of the results and mainly to computational cost, compared with the methods traditionally used in the seismic industry. There are several search strategies to determine these parameters, based on systematic searches and using optimization algorithms, where only one parameter at each time can be estimated, or the two or three parameters simultaneously. Taking in to account the search strategy by means of the application of global optimization, these three parameters can be estimated through of procedures: in the first case the three parameters can be simultaneously estimated and in second case initially two parameters can be determined simultaneously (β0, RNIP), and subsequently the third parameter (RN), using the values of the two parameters already known. In this work it is presented the application and comparison of four algorithms of global optimization to find the CRS optimal parameters: Simulated Annealing (SA), Very Fast Simulated Annealing (VFSA), Differential Evolution (DE) and Controlled Random Search – 2 (CRS2). As importants results of the application of each optimization method, as well as between the methods regarding the effectiveness, efficiency and reliability to determine the best CRS parameters are presented. Subsequently, applying the global search strategies for the determination of these parameters, by means of the optimization method VFSA that presented the best performance, the CRS stacking was applied to the Marmousi dataset, one stacking using two parameters (β0, RNIP), estimated by global search, and another CRS stacking using the three parameters (β0, RNIP, RN), also estimated by global search.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-10-20T15:48:49Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoAmplitudeVerdadeira.pdf: 9758204 bytes, checksum: 62046151ddeb3a3200cca375f8b96f32 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-10-22T13:46:36Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoAmplitudeVerdadeira.pdf: 9758204 bytes, checksum: 62046151ddeb3a3200cca375f8b96f32 (MD5) / Made available in DSpace on 2014-10-22T13:46:36Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoAmplitudeVerdadeira.pdf: 9758204 bytes, checksum: 62046151ddeb3a3200cca375f8b96f32 (MD5) Previous issue date: 1999 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Nos últimos anos tem-se verificado através de várias publicações um interesse crescente em métodos de migração com amplitude verdadeira, com o objetivo de obter mais informações sobre as propriedades de refletividade da subsuperfície da terra. A maior parte desses trabalhos tem tratado deste tema baseando-se na aproximação de Born, como em Bleistein (1987) e Bleistein et al. (1987), ou na aproximação do campo de ondas pela teoria do raio como Hubral et al. (1991), Schleicher et al. (1993) e Martins et al. (1997). Considerando configurações arbitrárias de fontes e receptores, as reflexões primárias compressionais podem ser imageadas em reflexões migradas no domínio do tempo ou profundidade de tal modo que as amplitudes do campo de ondas migrado são uma medida do coeficiente de reflexão dependente do ângulo de incidência. Para realizar esta tarefa, vários algoritmos têm sido propostos nos últimos anos baseados nas aproximações de Kirchhoff e Born. Essas duas abordagens utilizam um operador integral de empilhamento de difrações ponderado que é aplicado aos dados da seção sísmica de entrada. Como resultado obtém-se uma seção migrada onde, em cada ponto refletor, tem-se o pulso da fonte com amplitude proporcional ao coeficiente de reflexão naquele ponto. Baseando-se na aproximação de Kirchhoff e na aproximação da teoria do raio do campo de ondas, neste trabalho é obtida a função peso para modelos bidimensionais (2-D) e dois e meio dimensionais (2,5-D) que é aplicada a dados sintéticos com e sem ruído. O resultado mostra a precisão e estabilidade do método de migração em 2-D e 2,5-D como uma ferramenta para a obtenção de informações importantes da subsuperfície da terra, que é de grande interesse para a análise da variação da amplitude com o afastamento (ângulo). Em suma, este trabalho apresenta expressões para as funções peso 2-D e 2,5-D em função de parâmetros ao longo de cada ramo do raio. São mostrados exemplos da aplicação do algoritmo de migração em profundidade a dados sintéticos 2-D e 2,5-D obtidos por modelamento sísmico através da teoria do raio usando o pacote Seis88 (Cervený e Psencík, 1988) e os resultados confirmaram a remoção do espalhamento geométrico dos dados migrados mesmo na presença de ruído. Testes adicionais foram realizados para a análise do efeito de alongamento do pulso na migração em profundidade (Tygel et al., 1994) e a aplicação do empilhamento múltiplo (Tygel et al., 1993) para a estimativa de atributos dos pontos de reflexão - no caso o ângulo de reflexão e a posição do receptor. / In the recent past years we have seen through various published papers an increasing interest in true amplitude migration methods, in order to obtain more informations about the reflectivity properties of the earth subsurface. The most part of these works has treated of this thema either based on Born approximation as given by Beistein (1987) and Bleistein et al. (1987), or on ray theoretical wavefield approximation as given by Hubral et al. (1991), Schleicher et al. (1993) and Martins et al. (1997). By considering arbitrary source-receiver configurations the compressional primary reflections can be imaged into time or depth-migrated reflections so that the migrated wavefield amplitudes are a measure of angle-dependent reflection coeffients. In order to do this various migration algorithms were proposed in the recent past years based on Born or Kirchhoff approach. Both of them treats of a weighted diffraction stack integral operator that is applied to the input seismic data. As result we have a migrated seismic section where at each reflector point there is the source wavelet with the amplitude proportinal to the reflection coefficient at that point. Based on Kirchhoff approach, in this thesis we derive the weight function and the diffraction stack integral operator for the two dimensional (2-D) and for the two and one half (2.5-D) seismic model and apply it to a set of synthetic seismic data in noise environment. The result shows the accuracy and stability of the 2-D and 2.5-D migration methods as a tool for obtaining important information about the reflectivity properties of the earth subsurface, which is of great interest for the amplitude versus offset (angle) analysis. In summary, we present an expressions for the 2-D and 2.5-D weights as a function of parameters along each ray branch of the in-plane trajectory. Moreover, we show examples of application of the true-amplitude depth migration algorithm to synthetic seismic data obtained by ray theory seismic modeling using the Seis88 package (Cervený e Psencík, 1988), in order to make a numerical analysis and to verify the stability and accuracy of the algorithm. The results confirmed the removal of the geometrical spreading from migrated data, even in presence of noise. Additional tests were performed for pulse distortion analysis in depth rnigrated sections (Tygel et al., 1994) and to obtain reflection points attributes by multiple diffraction stack (Tygel et al., 1993).

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-10-20T14:02:56Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_Migracao3DTempo.pdf: 10471953 bytes, checksum: 415de3c21f90d60df1bc57b415e7a879 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-10-22T13:47:23Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_Migracao3DTempo.pdf: 10471953 bytes, checksum: 415de3c21f90d60df1bc57b415e7a879 (MD5) / Made available in DSpace on 2014-10-22T13:47:23Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_Migracao3DTempo.pdf: 10471953 bytes, checksum: 415de3c21f90d60df1bc57b415e7a879 (MD5) Previous issue date: 1998 / O trabalho aqui apresentado visa o estudo da migração 3D no tempo, em amplitudes verdadeiras, de dados da seção afastamento nulo ("zero-offset"), usando a aproximação de segunda ordem do tempo de trânsito de difração. O processo de migração no tempo usa a função difração determinada, corretamente, através de mesmos parâmetros que controlam a função tempo de trânsito de reflexão, além do auxílio dos raios paraxiais. Os parâmetros que controlam a função tempo de trânsito de reflexão são determinados através de pelo menos nove medidas de tempo de trânsito, obtidas a partir de dados pré-empilhados (tiro comum, ponto médio comum - CMP, etc.). Determinada a função difração, pode-se fazer a migração sem a necessidade de que se conheça o modelo de velocidades, pois a função tempo de trânsito de difração, como definido neste trabalho, não necessita do mesmo. O espalhamento geométrico é removido através do processo de migração e aplicação de uma função peso aos dados da seção afastamento nulo. Juntamente com a função peso, aplicou-se um filtro objetivando a recuperação da fase do sinal, ocorrida como consequência do processo de migração. Devido a limitação computacional, fez-se testes com dados sintéticos, em 2D, objetivando verificar a eficiência do programa. Os resultados obtidos foram satisfatórios, mostrando a eficiência e confiabilidade do processo. / This presentation aims at the 3-D time migration of zero offset data, in true amplitude. This method is based on paraxial ray theory and uses a diffraction time function which is directly and correctly determined by the measurement of pre-stack seismic data. It is not necessary to know a macro velocity model in order to apply the time migration. In order to obtain a true amplitude time migration the migration result must be multiplied by a scaling factor and convolved with a known function. Together with a scaling factor, a filter was applied in order to recover the signal phase altered during the migration process. Due to the computational limitation synthetic dada in 2-D was used aims to test the program efficiency. The result was satisfactory, showing the efficiency and robustness of process.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-10-31T16:47:53Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoFd3dProfundidade.pdf: 2169798 bytes, checksum: 22a97423796a43155cbe09c83257b278 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-11-03T12:25:54Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoFd3dProfundidade.pdf: 2169798 bytes, checksum: 22a97423796a43155cbe09c83257b278 (MD5) / Made available in DSpace on 2014-11-03T12:25:54Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoFd3dProfundidade.pdf: 2169798 bytes, checksum: 22a97423796a43155cbe09c83257b278 (MD5) Previous issue date: 2011 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / A implementação convencional do método de migração por diferenças finitas 3D, usa a técnica de splitting inline e crossline para melhorar a eficiência computacional deste algoritmo. Esta abordagem torna o algoritmo eficiente computacionalmente, porém cria anisotropia numérica. Esta anisotropia numérica por sua vez, pode levar a falsos posicionamentos de refletores inclinados, especialmente refletores com grandes ângulos de mergulho. Neste trabalho, como objetivo de evitar o surgimento da anisotropia numérica, implementamos o operador de extrapolação do campo de onda para baixo sem usar a técnica splitting inline e crossline no domínio frequência-espaço via método de diferenças finitas implícito, usando a aproximação de Padé complexa. Comparamos a performance do algoritmo iterativo Bi-gradiente conjugado estabilizado (Bi-CGSTAB) com o multifrontal massively parallel solver (MUMPS) para resolver o sistema linear oriundo do método de migração por diferenças finitas. Verifica-se que usando a expansão de Padé complexa ao invés da expansão de Padé real, o algoritmo iterativo Bi-CGSTAB fica mais eficientes computacionalmente, ou seja, a expansão de Padé complexa atua como um precondicionador para este algoritmo iterativo. Como consequência, o algoritmo iterativo Bi-CGSTAB é bem mais eficiente computacionalmente que o MUMPS para resolver o sistema linear quando usado apenas um termo da expansão de Padé complexa. Para aproximações de grandes ângulos, métodos diretos são necessários. Para validar e avaliar as propriedades desses algoritmos de migração, usamos o modelo de sal SEG/EAGE para calcular a sua resposta ao impulso. / Conventional implementations of 3D finite-difference (FD) migration use splitting techniques to accelerate performace and save computational cost. However, such techniques are plagued with numerical anisotropy that jeopardizes correct positioning of dipping reflectors in the directions not used for the splitted operators. We implement 3D downward continuation migration without splitting in the space coordinates using a complex Padé approximation and implicit finite differences. In this way, the numerical anisotropy is eliminated at the expense of a computationally more intensive solution of a large banded linear system. We compare the performance of the iterative stabilized biconjugate gradient (Bi-CGSTAB) and the multifrontal massively parallel direct solver (MUMPS). It turns out that the use of the complex Padé approximation provides an effective preconditioner for the Bi-CGSTAB, reducing the number of iterations relative to the real Padé expansion of the FD operator. As a consequence, the iterative Bi-CGSTABmethod ismore efficient than the directMUMPSmethodwhen solving for a single termin the Padé expansion. Forwide angle approximations direct methods are required. These algorithms are validated and the properties evaluated computing themigration impulse response in the SEG/EAGE salt model.