Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-18T16:54:07Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ImageamentoMultifocoRefletores.pdf: 15149810 bytes, checksum: b0d10801d9480960ffdf7f59a15484e1 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-18T17:14:34Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ImageamentoMultifocoRefletores.pdf: 15149810 bytes, checksum: b0d10801d9480960ffdf7f59a15484e1 (MD5) / Made available in DSpace on 2014-09-18T17:14:34Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ImageamentoMultifocoRefletores.pdf: 15149810 bytes, checksum: b0d10801d9480960ffdf7f59a15484e1 (MD5) Previous issue date: 2000 / DAAD - Serviço Alemão de Intercâmbio Acadêmico / Deutscher Akademischer Austauschdiens / A simulação de uma seção sísmica de afastamento-nulo (ZO) a partir de dados de cobertura múltipla para um meio 2-D, através do empilhamento, é um método de imageamento de reflexão sísmica muito utilizado, que permite reduzir a quantidade de dados e melhorar a relação sinal/ruído. Segundo Berkovitch et al. (1999) o método Multifoco está baseado na Teoria do Imageamento Homeomórfico e consiste em empilhar dados de cobertura múltipla com distribuição fonte-receptor arbitrária de acordo com uma nova correção de sobretempo, chamada Multifoco. Esta correção de sobretempo esta baseada numa aproximação esférica local da frente de onda focalizante na vizinhança da superfície da terra. Este método permite construir uma seção sísmica no domínio do tempo de afastamento nulo aumentando a relação sinal/ruído. A técnica Multifoco não necessita do conhecimento a priori de um macro-modelo de velocidades. Três parâmetros são usados para descrever a aproximação de tempo de trânsito, Multifoco, os quais são: 1) o ângulo de emergência do raio de afastamento nulo ou raio de reflexão normal (β0), 2) a curvatura da frente de onda no Ponto de Incidência Normal (RNIP) e 3) curvatura da frente de Onda Normal (RN). Sendo também necessário a velocidade próximo a superfície da terra. Neste trabalho de tese aplico esta técnica de empilhamento Multifoco para dados de cobertura múltipla referidos a modelos de velocidade constante e modelo heterogêneos, com o objetivo de simular seções sísmicas afastamento-nulo. Neste caso, como se trata da solução de um problema direto, o macro-modelo de velocidades é considerado conhecido a priori. No contexto do problema inverso tem-se que os parâmetros RNIP, RN e β0 podem ser determinados a partir da análise de coerência aplicada aos dados sísmicos de múltipla cobertura. Na solução deste problema a função objetivo, a ser otimizada, é definida pelo cálculo da máxima coerência existente entre os dados na superfície de empilhamento sísmico. Neste trabalho de tese nos discutimos a sensibilidade da aproximação do tempo de trânsito usado no empilhamento Multifoco, como uma função dos parâmetros RNIP, RN e β0. Esta análise de sensibilidade é feita de três diferentes modos: 1) a primeira derivada da função objetivo, 2) a medida de coerência, denominada semelhança, e 3) a sensibilidade no Empilhamento Multifoco. / The simulation of a zero-offset section (ZO) from multi-coverage seismic reflection data for a 2-D media, through the stack, is a widely used seismic reflection imaging method, that allows to reduce the amount of data and to improve the signal-to-noise ratio. According to Berkovitch et al. (1999) the Multifocusing method is based on Theory of Homeomorphic Imaging and consists of stacking multi-coverage seismic reflection data with arbitrary distribution source-receiver according to a new local moveout correction, called Multifocusing. This moveout correction is based on a local spherical approximation of the focusing wavefront in the vicinity of the surface of the earth. This method allows to build a seismic section in the domain of the time of a zero-offset increasing the signal-to-noise ratio. The Multifocusing technique does not need any knowledge a priori of the macro-model velocity. Three parameters are used to describe the Multifocusing moveout correction, which are: 1) emergence angle of the zero-offset ray or normal reflection ray (β0), 2) the wavefront curvature at the Point of Normal Incidence (RNIP and 3) the wavefront curvature of Normal Wave (RN). Being also necessary the near-surface velocity. In this thesis work I apply this Multifocusing Stack technique for multi-coverage seismic reflection data covering the models of constant velocity and heterogeneous model, with the objective of simulating of zero-offset sections. In this case as it is the solution of forward problem, this macro velociy model is considered know apriori. In the context of the inverse problem it is had the parameters, RNIP, RN and β0 and can be determining through the analysis of applied coherence to the multi-coverage seismic reflection data. In the solution of this problem the objective function, to be optimized, is defined by the calculation of the maximum existent coherence among the data in the surface of seismic stack. In this thesis we discuss the sensibility of the travetime, used in Multifocusing Stack, as a function of the parameters RNIP, RN and β0. This sensibility analysis is done of three different manners: 1) the first derivate of the objetive function, 2) the coherence measure, denominated semblance, and 3) the sensibility in the Multifocusing Stack.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-10-03T14:05:59Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MetodosGeofisicosAplicadosProspeccaoMineral.pdf: 16126636 bytes, checksum: 323e2f8be50081544aa5c2e91ed305a1 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-10-08T15:24:58Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MetodosGeofisicosAplicadosProspeccaoMineral.pdf: 16126636 bytes, checksum: 323e2f8be50081544aa5c2e91ed305a1 (MD5) / Made available in DSpace on 2014-10-08T15:24:58Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MetodosGeofisicosAplicadosProspeccaoMineral.pdf: 16126636 bytes, checksum: 323e2f8be50081544aa5c2e91ed305a1 (MD5) Previous issue date: 1982 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Universidade Nacional de Panamá / O presente trabalho de prospecção geofísica abrange o levantamento de uma faixa de direção NW-SE, desde o local denomina do Chatão, à margem direita do rio Gurupi, até 5 km além de Chega Tudo, no estado de Maranhão. Esse levantamento foi feito em travessas perpendiculares àquela faixa, a qual constitui uma extensa zona de falha. Os trabalhos de geofísica de campo, envolvendo os métodos magnético (campo total), cintilómétrico, eletro-resistivo (Schlumberger, Wenner e Half-Schlumberger) e E.M. (Slingram), foram realizados em duas etapas: a) levantamento de reconhecimento em uma área de 22 km2, onde foram aplicados magnetometria e cintilometria e b) levantamento de detalhe em uma área de 3,50 km2, onde foram aplicados eletro-resistividade e E.M. Os objetivos principais deste trabalho foram determinar e localizar as principais feições estruturais e geo-elétricas associadas às zonas mineralizadas, e selecionar, entre os métodos geofísicos empregados, qual o de melhor poder resolutivo na indicação de veios de quartzo em sub-superfície, possivelmente associados à mineração de ouro. O método magnético indicou uma anomalia principal no lado NE da área de reconhecimento, enquanto que o cintilométrico apresentou somente o valor do "background" da região. Por isso o método cintilométrico não será discutido nesse trabalho. Com os métodos eletro-resistivo e E.M. foram observadas várias anomalias ao longo de toda a área de detalhe. A interpretação dos dados de magnetometria e de eletro-resistividade foi realizada com técnicas computacionais, enquanto que para a dos dados eletromagnéticos foram empregados diagramas de fase. Os resultados apresentam uma estreita correlação das anomalias magnéticas e eletro-resistivas com a zona mineralizada. O método E.M. não apresentou nenhum poder resolutivo, uma vez que detectou anomalias semelhantes dentro e fora das áreas mais promissoras. Todas as anomalias foram confirmadas com furos de sondagem. Considerando-se que a área prospectada apresenta, devido ao difícil acesso, o mesmo grau de dificuldade para realização da pesquisa por métodos diretos (perfuração, poços ou trincheiras), como qualquer outra área similar na região amazônica, os resultados aqui apresentados confirmaram que os métodos magnético e eletro-resistivo podem ser utilizados com sucesso como ferramenta indireta na localização de zonas promissoras à exploração de ouro dentro da zona de falha de Chega Tudo ou em outros ambientes geológicos de semelhantes características. / This work of geophysical prospecting was carried out on a NW-SE geological lineament zone from Chatão - on the right side of Gurupi river up to 5 km beyond Chega Tudo in the State of Maranhão. The field work with the application of magnetic (total field), radiometric, electroresistivity (Schlumberger, Wenner and Half-Schlumbergere) and EM (Sligram) methods was carried out in two stages: a) reconnaissance survey of a 22 km2 area with magnetic and radiometric methods, and b) detailed survey of a 3.5 km2 area with electroresistivity and EM methods. The objectives were to study the main geoelectrical structures associated with mineralized zones, and to select among the used geophysical techniques the one with greater resolution to locate sub-surface quartz veins with gold mineralization. The magnetic method shows a large anomaly in the NE region of the area of reconnaissance, whereas the radiometric method indicates only the background value everywhere. Therefore the latter will be not discussed. The electroresistivity and EM surveys indicate many important anomalies on the detailed area. The electroresistivity and magnetic interpretation were carried out by computacional techniques whereas the EM data were interpreted by phase diagrams. The results show a good correlation between the magnetic and electroresistivity anomalies and the mineralized zone. The EM method was not found effective because it could not discriminate the anomalies inside and outside the mineralized zone. All anomalies were tested with well logs. Since the studied area has the same degree of difficulties for geological prospecting as any other similar area in the Amazon region, the results presented here show that the magnetic and electroresistivity methods, which turned out to be very effective as exploration tools for prospecting gold at Chega Tudo region, could be applied to any other environment with similar geological characteristics.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-10-07T13:08:10Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MapeamentoInterfaceSeparando.pdf: 20588420 bytes, checksum: 3a6fa78f0f81a9732415bdbf453cbc4a (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-10-09T14:18:57Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MapeamentoInterfaceSeparando.pdf: 20588420 bytes, checksum: 3a6fa78f0f81a9732415bdbf453cbc4a (MD5) / Made available in DSpace on 2014-10-09T14:18:57Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MapeamentoInterfaceSeparando.pdf: 20588420 bytes, checksum: 3a6fa78f0f81a9732415bdbf453cbc4a (MD5) Previous issue date: 1995 / O problema do mapeamento de uma interface separando dois meios homogêneos tem sido tradicionalmente resolvido através da continuação para baixo da anomalia gravimétrica observada. Este procedimento requer a utilização de um filtro passa-baixa ou de um parâmetro de amortecimento, que tendem a diminuir a resolução do relevo estimado da interface. Além disso, o uso da continuação para baixo implica supor (desnecessariamente) que a interface é uma superfície harmônica. Desta forma, a utilização da continuação para baixo é restrita a interfaces que apresentam um relevo suave relativamente à sua profundidade média. As restrições impostas pela continuação para baixo do campo gravimétrico são reduzidas na técnica de inversão proposta para este estudo, cujo objetivo é mapear a interface crosta-manto baseada na incorporação de informações a priori sobre o mínimo momento de inércia da fonte causadora, em relação a um eixo de concentração de massa pré-estabelecido e coincidente com o topo da fonte anômala. Além da introdução de informações sobre o mínimo momento, foi incluída também a possibilidade de haver variação lateral e em profundidade de densidade em setores previamente estabelecidos, o que permite estender a aplicação desta técnica a situações geológicas mais complexas que envolvam a reconstituição de corpos justapostos, que sejam de diferentes naturezas, tais como aqueles presentes em zonas de transição crustal. O método de inversão apresentado foi avaliado em testes usando dados sintéticos e os resultados obtidos mostraram a eficiência desta técnica em recuperar corpos alongados horizontalmente, que apresentem como topo ou base, interfaces arbitrárias a serem mapeadas. Foi também constatada a eficácia do método em reconstituir corpos que comportem variação lateral de densidade que pudessem retratar situações em que se deseja investigar a natureza do embasamento em regiões de transição crustal. Com base nos resultados destes testes, a técnica foi então empregada na inversão gravimétrica de cinco perfis do LEPLAC IV que cruzam a região oceânica adjacente à margem continental sudeste brasileira. Os resultados desta aplicação permitiram inferir a natureza dos corpos que constituem o arcabouço estrutural da região onde encontram-se inseridas a Dorsal de São Paulo e a Zona de Fratura de Porto Alegre, assim como mapear a Moho subjacente. / Mapping an arbitrary interface separating two homogeneous media, traditionally uses the downward continuation of the observed gravity anomaly. This procedure requires either a low-pass filtering or the introduction of a damping parameter. In both cases there is a loss in the resolution of the estimated interface relief. Moreover, the use of the downward continuation operator implies the unnecessary assumption that the interface is a harmonic surface. As a result, the use of the downward continuation operator constrains the relief mapping to interfaces which are smooth relative to its average depth. This difficulty is reduced in the inversion procedure presented in this thesis, whose aim is to map the crust-mantle interface using the minimum moment of inertia method by assuming a single horizontal axis of mass concentration, coinciding with the earth's surface. Horizontal and vertical density variations may be taken into account so that the method may be extended to complex geological settings involving the interpretation of adjacent bodies having different origins such as those present in crustal transitional zones. The perfomance of the inversion method was evaluated using sinthetic data and the results confirmed its efficiency in recovering the shape of horizontally elongated sources whose top or bottom may be represented by an arbitrary surface. The results also confirmed the method's efficiency in mapping the shape of bodies presenting lateral density variation, reflecting, therefore, a geological situation where the nature of the basement of the crustal transitional zone is being investigated. Based on these tests, the method was applied to the gravity inversion of five profiles from LEPLAC-IV, offshore the southeastern brazilian continental margin. The results permitted to infer the type of bodies making up the structural framework of the region comprising the São Paulo Ridge and the Porto Alegre Fracture Zone; in addition, it was possible to map the Moho in this region.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-10-31T16:46:15Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoEquacaoOnda.pdf: 1000885 bytes, checksum: 3e58b46bd6f6e3d79fcdf93d365719a0 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-11-03T12:15:28Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoEquacaoOnda.pdf: 1000885 bytes, checksum: 3e58b46bd6f6e3d79fcdf93d365719a0 (MD5) / Made available in DSpace on 2014-11-03T12:15:28Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoEquacaoOnda.pdf: 1000885 bytes, checksum: 3e58b46bd6f6e3d79fcdf93d365719a0 (MD5) Previous issue date: 2007 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Neste trabalho avaliamos uma classe de operadores de continuação de campos de onda, baseados em equações unidirecionais e com aplicação direta à migração sísmica. O método de representação de equações de onda unidirecionais, desenvolvido neste trabalho, é válido para abertura angular arbitrária, baseia-se no conceito de rigidez de um semiespaço, na transformação Dirichlet-Neumann e em sua discretização por elementos finitos. O método de construção dos operadores de continuação requer a introdução de variáveis auxiliares cujo número cresce em função da maior abertura angular desejada para o operador. Efetuamos a implementação no domínio do espaço e da frequência o que permite sua imediata paralelização. Baseados em experimentos numéricos, que avaliam a relação de dispersão e a resposta ao impulso do operador, propomos prescrições que permitem especificar o número de variáveis auxiliares e o passo de continuação para o operador de migração. A aplicação do algoritmo nos dados do modelo de domo salino da SEG-EAGE demonstra a capacidade do algoritmo em migrar refletores com forte mergulho em meios com forte variação lateral de velocidade. / The arbitrarily wide angle wave equations are applied to post-stack migration. This class of wavefield continuation operators are based on the concept of semi-space stiffness, Dirichlet to Neumann transformations and its discretization using finite elements. The specification of the continuation operator requires auxiliar variables. The number of these variable increase with the maximum angular aperture. The algorithm is implemented in space-frequency domain which allows straightforward parallelization. Based on numerical experiments using the dispersion relation and the impulse response of these operators we prescribe the number of auxiliar variables and downward continuation step. The migration of the zero offset SEG-EAGE data using the wide angle continuation operators shows their efectiveness for imaging steep dips in a complex velocity model.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-10-31T16:48:42Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoFdFfd.pdf: 1152510 bytes, checksum: 3aee1c0b06480886813dd15b93fe803d (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-11-03T12:34:41Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoFdFfd.pdf: 1152510 bytes, checksum: 3aee1c0b06480886813dd15b93fe803d (MD5) / Made available in DSpace on 2014-11-03T12:34:41Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoFdFfd.pdf: 1152510 bytes, checksum: 3aee1c0b06480886813dd15b93fe803d (MD5) Previous issue date: 2007 / Métodos de migração baseados na equação da onda unidirecional apresentam limitações no imageamento de refletores com forte mergulho e no tratamento de ondas evanescentes. Utilizando a expansão de Padé complexa na solução da equação da onda unidirecional para imageamento em geofísica é possível tratar ondas evanescentes, estabilizando a migração e dispensando um tratamento especial para absorção nas fronteiras do domínio. Utilizando várias parcelas na expansão de Padé, refletores com forte mergulho podem ser corretamente migrados. Nesta dissertação aplicamos a equação da onda unidirecional com expansão de Padé complexa para implementar dois algoritmos de migração em profundidade pré e pós empilhamento: a migração por diferenças finitas (FD) e a migração por diferenças finitas e Fourier (FFD). O estudo das curvas de dispersão e da resposta impulsiva dos operadores de migração nos permitiu escolher o número adequado de parcelas na expansão de Padé e os coeficientes da expansão que garantem estabilidade à continuação do campo de onda para um ângulo de mergulho máximo prescrito. As implementações foram validadas nos dados Marmousi e no modelo de domo de sal da SEG/EAGE mostrando que refletores com forte mergulho foram corretamente migrados, mesmo na presença de forte variação lateral de velocidade. Esses resultados são comparados com outros métodos de migração baseados na equação da onda unidirecional ressaltando a qualidade da aproximação estudada neste trabalho. / Seismic Migration by downward continuation using the unidirectional wave equation approximations has two shortcomings: imaging steep dip reflectors and handling evanescent waves. Complex Padé approximations allow a better treatment of evanescent modes stabilizing the finite difference migration, and does not require special treatment for domain boundaries. Imaging of steep dip reflectors can be improved using several terms in the Padé expansion. This dissertation discuss the implementation and evaluation of complex Padé approximation for finite difference migration and Fourier finite difference migration. The study of the dispersion relation and impulsive response associated to the migration operator is used to select the number of terms and coefficients in the Padé expansion which assures stability for a prescribed maximum reflector dip. The implementations are validated in the Marmousi and SEG/EAGE salt model datasets, and compared to other wave equation migration methods. The results of FD and FFD complex Padé migrations can handle steeper dips, and present a much lower signal to noise ratio than their real valued counterparts.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-10-31T16:49:13Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoProfundidadeSolucao.pdf: 2614276 bytes, checksum: 4c6495ea7cb216062cd9c3623481dc55 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-11-03T13:08:31Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoProfundidadeSolucao.pdf: 2614276 bytes, checksum: 4c6495ea7cb216062cd9c3623481dc55 (MD5) / Made available in DSpace on 2014-11-03T13:08:31Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_MigracaoProfundidadeSolucao.pdf: 2614276 bytes, checksum: 4c6495ea7cb216062cd9c3623481dc55 (MD5) Previous issue date: 2001 / ANP - Agência Nacional do Petróleo / Nos últimos anos tem-se verificado um interesse crescente no desenvolvimento de algoritmos de imageamento sísmico com a finalidade de obter uma imagem da subsuperfície da terra. A migração pelo método de Kirchhoff, por exemplo, é um método de imageamento muito eficiente empregado na busca da localização de refletores na subsuperficie, quando dispomos do cálculo dos tempos de trânsito necessários para a etapa de empilhamento, sendo estes obtidos neste trabalho através da solução da equação eiconal. Primeiramente, é apresentada a teoria da migração de Kirchhoff em profundidade baseada na teoria do raio, sendo em seguida introduzida a equação eiconal, através da qual são obtidos os tempos de trânsitos empregados no empilhamento das curvas de difrações. Em seguida é desenvolvido um algoritmo de migração em profundidade fazendo uso dos tempos de trânsito obtidos através da equação eiconal. Finalmente, aplicamos este algoritmo a dados sintéticos contendo ruído aditivo e múltiplas e obtemos como resultado uma seção sísmica na profundidade. Através dos experimentos feitos neste trabalho observou-se que o algoritmo de migração desenvolvido mostrou-se bastante eficiente e eficaz na reconstrução da imagem dos refletores. / In the last years we have seen an increasing interest in seismic imaging algorithms in order to obtain better informations about the earth interior. The Kirchhoff migration method is very useful for determining the position of seismic reflectors, if is known the seismic wave velocity model and the traveltimes are well determined through the earth model. The traveltime calculation is a necessary step for stacking the seismic data by means of the Kirchhoff migration operator. In this work the traveltimes are obtained by solving the eiconal equation of the ray theory. At first, the theory of Kirchhoff migration is reviewed, by considering depth migration in heterogeneous media with arbitrary curved reflectors. Secondly, the numerical solution of the eiconal equation is presented for transmited, diffracted and head waves. There offer, the depth migration algorithm is presented, must makes use of traveltimes obtained by the eiconal equation. Finally, the developed migration algorithm is applied to synthetic models, providing a very good image resolution in comparison with the conventional ray tracing migration methods, even in the presence of random or coherent (multiple reflections) noise.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-10-31T16:50:17Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_Modelagem25DElementos.pdf: 1680198 bytes, checksum: 53aafa3ebd7bc308d0665ed0a3d2d22a (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-11-03T13:19:12Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_Modelagem25DElementos.pdf: 1680198 bytes, checksum: 53aafa3ebd7bc308d0665ed0a3d2d22a (MD5) / Made available in DSpace on 2014-11-03T13:19:12Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_Modelagem25DElementos.pdf: 1680198 bytes, checksum: 53aafa3ebd7bc308d0665ed0a3d2d22a (MD5) Previous issue date: 2014 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Este trabalho consiste em realizar a modelagem, via elementos finitos (EF) 2,5D, do efeito da topografia do terreno sobre dados obtidos com o método eletromagnético a multi-frequência (EMMF). Este método usa como fonte uma grande espira quadrada de corrente elétrica com centenas de metros de lado, e como receptores, bobinas posicionadas na horizontal em alinhamento com o transmissor. A subsuperfície é representada por heterogeneidades bidimensionais imersas em um meio horizontalmente estratificado. A formulação, partindo das equações de Maxwell, é desenvolvida a partir da separação do campo eletromagnético em primário (campos no hospedeiro multi-estratificado) e secundário (diferença entre o campo total e o primário). O domínio discretizado é descrito por uma malha não estruturada, com elementos triangulares. Para calcular as componentes derivadas da solução de elementos finitos, em um determinado nó da malha, foi usada a média aritmética das derivadas das funções bases de EF em torno daquele nó. O código de modelagem construído permite quantificar e analisar como os gradientes topográficos influenciam as medidas dos campos eletromagnéticos gerados. A aplicação é a avaliação dessas influências sobre a componente radial do campo da espira na superfície terrestre, que é a componente empregada no método eletromagnético a multi-frequência (EMMF). / This dissertation shows the numerical modeling results of the effect of the terrain topography on the multi-frequency electromagnetic method (EMMF). In this method the source is a very large square current loop, with a side of a few hundreds meters. The used numerical method is the Finite Element (FE) method in a 2,5D formulation. The geo-electrical models are composed of two-dimensional structures inside layered one-dimensional media. The formulation uses the separation of the EM field in two parts called the primary field (the field of the same source in the layered medium only) and the secondary field (the difference between the complete field in the 2D model and the primary field), to avoid the difficulty of representing the singular source in the discrete representation of the media in the form of the FE mesh. The discretized domain is build with a non-structured triangular mesh. In order to calculate the field components that are derived from those given by the FE solution in a given node, it is necessary to calculate numerical derivatives on the discrete solutions. Here we have used the average value of the derivatives of the basis functions around the node, which is a simple and effective way of estimating the derivative. The resulting modeling code allow us to quantitatively analyse the influence of topographic gradients on the measured electromagnetic fields from this method.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-11-14T13:29:37Z No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelagemMcsemDominio.pdf: 2587982 bytes, checksum: 128281a90d0470ea7eb6ca5d50968966 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-11-14T17:15:28Z (GMT) No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelagemMcsemDominio.pdf: 2587982 bytes, checksum: 128281a90d0470ea7eb6ca5d50968966 (MD5) / Made available in DSpace on 2014-11-14T17:15:28Z (GMT). No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelagemMcsemDominio.pdf: 2587982 bytes, checksum: 128281a90d0470ea7eb6ca5d50968966 (MD5) Previous issue date: 2009 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / A modelagem do mCSEM é feita normalmente no domínio da frequência, desde sua formulação teórica até a análise dos resultados, devido às simplificações nas equações de Maxwell, possibilitadas quando trabalhamos em um regime de baixa frequência. No entanto, a abordagem através do domínio do tempo pode em princípio fornecer informação equivalente sobre a geofísica da subsuperfície aos dados no domínio da frequência. Neste trabalho, modelamos o mCSEM no domínio da frequência em modelos unidimensionais, e usamos a transformada discreta de Fourier para obter os dados no domínio do tempo. Simulamos ambientes geológicos marinhos com e sem uma camada resistiva, que representa um reservatório de hidrocarbonetos. Verificamos que os dados no domínio do tempo apresentam diferenças quando calculados para os modelos com e sem hidrocarbonetos em praticamente todas as configurações de modelo. Calculamos os resultados considerando variações na profundidade do mar, na posição dos receptores e na resistividade da camada de hidrocarbonetos. Observamos a influência da airwave, presente mesmo em profundidades oceânicas com mais de 1000m, e apesar de não ser possível uma simples separação dessa influência nos dados, o domínio do tempo nos permitiu fazer uma análise de seus efeitos sobre o levantamento. Como parte da preparação para a modelagem em ambientes 2D e 3D, fazemos também um estudo sobre o ganho de desempenho pelo uso do paralelismo computacional em nossa tarefa. / The mCSEM modelling is usually done in the frequency domain, from its theoretical formulation to the analysis of the results. However, the time domain approach is, in principle, capable of providing equivalent information about the geo-electric structure of the subsurface. In this work, we model frequency domain mCSEM data in 1-D environments, then we perform the discrete Fourier transform to obtain time domain results. We simulated marine geological environments with and without the resistive layer that represents the hydrocarbon reservoir. We verified that the time domain data are significantly different when calculated for models with and without hydrocarbons in almost all model configurations. We calculated the results considering variations in the sea depth, in the position of the receivers and in the resistivity of the hydrocarbon layer. We observed the influence of the airwave, even at sea depths greater than 1000m, and although a simple separation of this influence on data is not possible, the time domain allowed us to do an analysis of its effects on the survey. As part of the preparation for the 2-D and 3-D modelling, we also have studied the gain in performance from the use of parallel processing in our task.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-11-14T13:30:11Z No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelagemFiltragemDistorcoes.pdf: 7630304 bytes, checksum: f4663ae7412a96d4ebf62dbb11a7a564 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-11-17T12:55:55Z (GMT) No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelagemFiltragemDistorcoes.pdf: 7630304 bytes, checksum: f4663ae7412a96d4ebf62dbb11a7a564 (MD5) / Made available in DSpace on 2014-11-17T12:55:55Z (GMT). No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelagemFiltragemDistorcoes.pdf: 7630304 bytes, checksum: f4663ae7412a96d4ebf62dbb11a7a564 (MD5) Previous issue date: 1994 / PETROBRAS - Petróleo Brasileiro S.A. / As distorções estáticas são um grave problema que afeta o método magnetotelúrico e muitas tentativas têm sido feitas para eliminar ou minimizar os seus efeitos. Este trabalho trata de uma técnica nova que trata deste problema, o método EMAP, que é uma adaptação do método magnetotelúrico na qual as medidas de campo elétrico são feitas em uma linha contínua de dipolos conectados entre si, e os dados coletados desta maneira são tratados com um filtro espacial passa-baixa dependente da frequência. Este trabalho é composto de duas partes principais, a simulação numérica dos dados contaminados com as distorções estáticas e a filtragem destes dados com o filtro espacial passa-baixa. Na primeira parte, aplicamos o método dos elementos finitos para simular a resposta dos dipolos elétricos, os quais fazem as medidas do campo elétrico. Na segunda parte aplicamos a janela de Hanning como filtro passa-baixa apropriado para tratar os dados. / Static distortions are a serius problem to the magnetotelluric method and many atempts to eliminate or to minimize them in the MT data have been tried with more or less success. This work is an overview on a new technique that deals with this problem, the EMAP method, witch is a modification of the traditional MT method in witch the measures of the electric field are made in a continuous line of dipoles equally spaced, and the data so obtained are filtered with a spatial low-pass, frequency dependent digital filter. This work is composed of two main parts, the numerical simulation of electric field data distorted by the static shift and the aplication of the spatial low-pass filter to minimize the static shift. In the first part, we use the finite elements method to simulate the response of the electric dipole witch actually makes the electric field measurements; in the second part, we aply the Hanning window as the spatial filter to treat the data.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-11-17T12:23:15Z No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoCorrecaoDescentralizacao.pdf: 1000871 bytes, checksum: 9e13030d31de4a03b7be9f569ee09e76 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-11-17T13:30:40Z (GMT) No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoCorrecaoDescentralizacao.pdf: 1000871 bytes, checksum: 9e13030d31de4a03b7be9f569ee09e76 (MD5) / Made available in DSpace on 2014-11-17T13:30:40Z (GMT). No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoCorrecaoDescentralizacao.pdf: 1000871 bytes, checksum: 9e13030d31de4a03b7be9f569ee09e76 (MD5) Previous issue date: 2003 / ANP - Agência Nacional do Petróleo / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / As imagens da parede do poço obtidas através da ferramenta UBI (Ultrasonic Borehole Imager) são amplamente utilizadas por geólogos e engenheiros de petróleo para identificar eventos geológicos nas paredes de poços abertos e na inspeção dos tubos de revestimento, uma vez que este perfil praticamente fornece uma fotografia da parede do poço. As ferramentas de imageamento acústico produzem imagens do tempo de trânsito e da amplitude do pulso acústico gerado pela ferramenta e refletido na parede do poço. Entretanto, estas imagens podem ter uma interpretação não realista, uma vez que elas podem ter seu aspecto alterado em razão do movimento da ferramenta no interior do poço. Este trabalho apresenta o modelamento das imagens de tempo de trânsito a partir da aplicação do critério de Coulomb para a ruptura da parede do poço submetida a um estado plano de tensões, a qual fornecerá a seção do poço, que é a forma geométrica que será traçada pela ferramenta de imageamento acústico do poço. O deslocamento ascensional da ferramenta e as imperfeições da parede do poço, normalmente, são os responsáveis pelo deslocamento da posição do transdutor em relação ao eixo do poço. Este efeito pode ter grande responsabilidade nas imperfeições das imagens de tempo de trânsito. Assim, a correção dessas imagens, chamada de correção da descentralização, busca o reposicionamento do transdutor para a posição do eixo do poço. Apresenta-se, também, um método de correção do efeito da descentralização da ferramenta baseado neste modelamento. O método é proposto com base na geometria analítica plana e no método do raio para a definição do tempo de trânsito do pulso acústico, com o objetivo de reconstruir as imagens de tempo de trânsito obtidas com a ferramenta descentralizada, ou seja, corrigir estas imagens tornando-as como se fossem obtidas com a ferramenta centralizada em relação ao eixo do poço. / The imaging tools used to borehole wall features identification have been largely utilized by geologist and oil engineers to identify geological events in a open hole and inspect the casing tube. The acoustic borehole imaging tools generate a transit time image and an acoustic amplitude image that are used to this proposes. However those logs may have a non-realist interpretation, since some tools effect can negatively influence in the images appearance. This paper presents a transit time image model starting from the application of the Coulomb’s approach to the borehole wall rupture submitted to a plane state of tensions which will supply the borehole section that is the geometric form that will be mapped by the acoustic borehole imaging toll. The tool up displacement and the borehole wall imperfections are usually the responsible for the transducer displacement in relation to the borehole axis. This effect may have important responsibility in the acoustic images imperfections. Thus, a computational process of transducer repositioning in the borehole axis position obtains the correction of those images, called decentralization correction. A method of tool decentralization effect correction is presented too based on this model which is proposed based on the plane analytic geometry and in the ray method to the definition of the transit time of the acoustic pulse, with the objective of reconstruct the transit time images achieved by the decentralized tool, that is to say, correct these images becoming as they were achieved by the centralized tool in relation to the borehole axis.