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

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-05T13:45:51Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseComparativaTres.pdf: 6132554 bytes, checksum: 3f873423cc9d66ad912a70a0cf13c3b9 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:17:46Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-29T13:39:59Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseComparativaTres.pdf: 6132554 bytes, checksum: 3f873423cc9d66ad912a70a0cf13c3b9 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-09T14:38:20Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseComparativaTres.pdf: 6132554 bytes, checksum: 3f873423cc9d66ad912a70a0cf13c3b9 (MD5) / Made available in DSpace on 2014-09-09T14:38:20Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseComparativaTres.pdf: 6132554 bytes, checksum: 3f873423cc9d66ad912a70a0cf13c3b9 (MD5) Previous issue date: 1985 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / FINEP - Financiadora de Estudos e Projetos / Tradicionalmente, o método dos mínimos quadrados tem sido empregado na inversão não linear de dados de campo potencial. No caso em que as observações dos campos gravimétrico ou magnético contém apenas ruído Gaussiano. O método dos mínimos quadrados não apresenta problemas. Entretanto, quando as observações são perturbadas por ruído não Gaussiano, ou mesmo por ruído não aleatório, como é o caso de muitos ruídos geológicos, o método dos mínimos quadrados torna-se bastante ineficiente, e métodos alternativos devem ser empregados a fim de produzir interpretações realísticas. Neste trabalho, uma comparação é feita entre os métodos dos mínimos quadrados, dos mínimos absolutos e do ajuste-M, aplicados à inversão não linear de dados de campo potencial. A comparação é efetuada usando-se dados teóricos, onde diversas situações geológicas são simuladas. Os resultados mostram que na presença de ruído geológico, caracterizado por pequeno corpo raso acima do corpo principal, ou por corpo grande, adjacente ao corpo principal, o ajuste-M apresenta desempenho muito superior ao dos mínimos quadrados e dos mínimos absolutos. Na presença de ruído Gaussiano, entretanto, o ajuste-M tem um desempenho inferior aos outros dois métodos. Como o ruído Gaussiano é um ruído branco, parte dele pode ser removido por um filtro passa baixa adequado, sem muita perda do sinal, o que não ocorre com o ruído geológico que contém componentes importantes de baixo número de onda. Desse modo o ajuste-M se torna uma ferramenta importante na interpretação de áreas geologicamente complexas, onde é comum a contaminação das anomalias por ruído geológico. Os três métodos em estudo são aplicados a uma anomalia magnética real causada por uma intrusão de diabásio em forma de dique, em sedimentos arenosos da formação Piauí na Bacia do Parnaíba. Os três métodos apresentaram resultados semelhantes indicando que tanto o nível de ruído Gaussiano como geológico são baixos nesta anomalia. / Nonlinear inversion of potential field data has been traditionally accomplished by the least squares method. As far as the anomalous field is corrupted by Gaussian random noise. Least squares inversion has a good performance. However, when the data are contaminated by non Gaussian noise, which is the case of most geological noise, the least squares method presents an extremely poor performance. As a result, alternative methods must be employed in this case in order to produce realistic and meaningful interpretations. This paper presents a comparison among the least squares method, the minimum absolute error and M-fitting applied to non linear inversion of potential field data. The analysis is performed using theoretical data generated by synthetic models simulating several geological settings. The results show that in the presence of geological noise represented either by small shallow bodies above the main body, or by large interfering bodies adjacent to the main body, M-fitting presents a much batter performance as compared with the least squares or the minimun absolute error methods. In the presence of Gaussian random noise, however M-fitting has a poor performance. Since a Gaussian noise is a white noise, a law pass filter applied to the observed data would remove part of the Gaussian noise with a minimum loss of the low wavenumber signal. On the other hand, most geological noise have important low wavenumber spectral components so that this noise cannot be eliminated without a significant loss of signal. Therefore, the M-fitting method may become an important interpretation tool when applied to complex areas (where anomalies are usually contaminated by geological noise) provided the data have been previously filtered by a suitable law pass filter. All three methods analysed in the paper are applied to a real magnetic anomaly due to a dike of diabasic rock intruded in sandstones and shales from the Piauí Formation, in Parnaíba Basin, Brazil. All three methods yielded similar interpretations which are consistent with the available a priori geological information. The fact that all methods produced similar results indicates a low level of geological and Gaussian random noise in the data.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-10T15:42:02Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoAlgoritmosConversao.pdf: 3250632 bytes, checksum: a0e70466bb0b3ec4d4b51f6d03df6a9b (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Incluir palavras-chave on 2014-08-07T13:31:22Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-10T12:58:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoAlgoritmosConversao.pdf: 3250632 bytes, checksum: a0e70466bb0b3ec4d4b51f6d03df6a9b (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T14:49:27Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoAlgoritmosConversao.pdf: 3250632 bytes, checksum: a0e70466bb0b3ec4d4b51f6d03df6a9b (MD5) / Made available in DSpace on 2014-09-10T14:49:27Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoAlgoritmosConversao.pdf: 3250632 bytes, checksum: a0e70466bb0b3ec4d4b51f6d03df6a9b (MD5) Previous issue date: 2013 / Ainda hoje, a migração em tempo é o processo de imageamento substancialmente empregado na indústria do petróleo. Tal popularidade é devida ao seu alto grau de eficiência e robustez, além de sua habilidade em focalizar refletores nos mais variados ambientes geológicos. Entretanto, em áreas de alta complexidade geológica a migração em tempo falha de tal forma que a migração em profundidade e um campo de velocidade em profundidade são indispensáveis. Esse campo é geralmente obtido através de processos tomográficos partindo de um campo de velocidade inicial. A conversão de campos de velocidade de tempo para profundidade é uma forma rápida de se obter um campo inicial mais consistente geologicamente para tais processos. Alguns algoritmos de conversão tempo-profundidade recentemente desenvolvidos baseados no traçamento de raios-imagem são revistos e um algoritmo alternativo baseado na propagação da frente de onda-imagem é proposto. Os algoritmos são aplicados a dados sintéticos bidimensionais e avaliados de acordo com suas eficiência e acurácia, destacando suas vantagens, desvantagens e limitações na obtenção de campos de velocidade em profundidade. / Even today, time migration is the imaging process substantially employed in the oil industry. Such popularity is due to its high degree of efficiency and robustness, and its ability to focus reflectors in various geological environments. However, in areas of high geological complexity, time migration fails so that depth migration and a depth velocity field are indispensable. This field is usually obtained through tomographic iterative processes starting from an initial velocity field. The conversion of velocity fields from time to depth is a quick way of obtaining a geologically more consistent initial velocity field for such processes. Some algorithms for time-to-depth conversion based on newly developed image-ray tracing are reviewed and an alternative algorithm based on propagation of the image-wavefront is proposed. The algorithms are applied to two-dimensional synthetic data and evaluated according to their efficiency and accuracy, highlighting their advantages, disadvantages and limitations in obtaining depth velocity fields.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-12T13:56:10Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoFiltrosVelocidade.pdf: 4324452 bytes, checksum: 90dbfc8b340951c647cbf95c43f4bab3 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar as palavras-chave on 2014-08-07T13:38:08Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-10T13:51:07Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoFiltrosVelocidade.pdf: 4324452 bytes, checksum: 90dbfc8b340951c647cbf95c43f4bab3 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T14:49:50Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoFiltrosVelocidade.pdf: 4324452 bytes, checksum: 90dbfc8b340951c647cbf95c43f4bab3 (MD5) / Made available in DSpace on 2014-09-10T14:49:50Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoFiltrosVelocidade.pdf: 4324452 bytes, checksum: 90dbfc8b340951c647cbf95c43f4bab3 (MD5) Previous issue date: 2004 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / A motivação geológica deste trabalho reside no imageamento de estruturas de bacias sedimentares da região Amazônica, onde a geração e o acúmulo de hidrocarboneto estão relacionados com a presença de soleiras de diabásio. A motivação sísmica reside no fato de que essas rochas intrusivas possuem grandes contrastes de impedância com a rocha encaixante, o que resulta em múltiplas, externas e internas, com amplitudes semelhantes as das primárias. O sinal sísmico das múltiplas podem predominar sobre o sinal das reflexões primárias oriundas de interfaces mais profundas, o que pode dificultar o processamento, a interpretação e o imageamento da seção sísmica temporal. Neste trabalho, estudamos a atenuação de múltiplas em seções sintéticas fonte-comum (FC) através da comparação de dois métodos. O primeiro método resulta da combinação das técnicas Wiener-Hopf-Levinson de predição (WHLP) e o de empilhamento superfície-de-reflexão-comum (CRS), e denominando WHLP-CRS, onde o operador é desenhado exclusivamente no domínio do tempo-espaço. O segundo método utilizado é o filtro de velocidade (ω-k) aplicado após o empilhamento superfície-de-reflexão (CRS), onde o operador é desenhado exclusivamente no domínio bidimensional de freqüência temporal-espacial. A identificação das múltiplas é feita na seção de afastamento-nulo (AN) simulada com o empilhamento CRS, e utiliza o critério da periodicidade entre primária e suas múltiplas. Os atributos da frente de onda, obtidos através do empilhamento CRS, são utilizados na definição de janelas móveis no domínio tempo-espaço, que são usadas para calcular o operador WHLP-CRS. O cálculo do filtroω-k é realizado no domínio da freqüência temporal-espacial, onde os eventos são selecionados para corte ou passagem. O filtro (ω-k) é classificado como filtro de corte, com alteração de amplitude, mas não de fase, e limites práticos são impostos pela amostragem tempo-espaço. Em termos práticos, concluímos que, para o caso de múltiplas, os eventos separados no domínio x-t não necessariamente se separam no domínio ω-k, o que dificulta o desenho de um operador ω-k semelhante em performance ao operador x-t. / The geological motivation of this work is the imaging of sedimentary basin structures of the Amazon region, where the generation and accumulation of hydrocarbons is related to the existence of diabase sills. The seismic motivation is the fact that these intrusive rocks present a great impedance contrast with respect to the host rock, what gives rise to external and internal multiples, with primary-like amplitudes. The seismic signal of the multiples can predominate over the primary reflection signals from deeper interfaces, making difficult the processing, interpretation and imaging of seismic sections. In this work we study the attenuation of multiples in common shot (CS) sections by the comparison of two methods. The first one is the combination of the Wiener-Hopf-Levinson (WHLP) and the common reflection surface (CRS) stacking techniques, here called WHLPCRS, where the operator is exclusively designed in the space-time domain. The second method is a velocity filter (ω-k), applied after the CRS stacking, where the operator is exclusively designed in the frequency-wavenumber domain. The identification of the multiples is performed on the zero-offset (ZO) section simulated by the CRS stacking, using the periodicity between primaries and its multiples. The wavefront attributes, obtained through CRS stacking, are inserted on movable space-time windows, used to calculate the WHLP-CRS operator. The ω-k filter calculations are performed in the frequency-wavenumber domain, where the events of interest are selected for cutting or passage. The ω-k filter is classified as a cut-off filter, with amplitude alteration and preservation of phase, the limits of it are imposed by a space-time sampling. In practical aspects we conclude that for the case of multiples, separated events on the x-t domain are not necessarily separated on the ω-k domain, which raise difficulties in the designing of a ω-k operator with a similar performance when compared to the x-t operator.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-28T15:36:55Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Modelagem1d2.5dDados.pdf: 6036247 bytes, checksum: 09857d990695852fd78d1d492571ec9e (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:23:28Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-27T14:13:55Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Modelagem1d2.5dDados.pdf: 6036247 bytes, checksum: 09857d990695852fd78d1d492571ec9e (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T16:36:35Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Modelagem1d2.5dDados.pdf: 6036247 bytes, checksum: 09857d990695852fd78d1d492571ec9e (MD5) / Made available in DSpace on 2014-09-10T16:36:35Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Modelagem1d2.5dDados.pdf: 6036247 bytes, checksum: 09857d990695852fd78d1d492571ec9e (MD5) Previous issue date: 2014 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / INCT/GP - Instituto Nacional de Ciência e Tecnologia de Geofísica do Petróleo / Neste trabalho apresentamos a solução do campo eletromagnético gerado por um dipolo elétrico horizontal em meios transversalmente isotrópicos com eixo de simetria vertical (TIV) e com eixo de simetria inclinado (TII). Para modelos unidimensionais, o campo eletromagnético foi obtido por duas metodologias distintas: (1) solução semi-analítica das equações de Maxwell com auxílio de potenciais vetores no caso TIV e (2) em modelos com anisotropia transversal inclinada o campo eletromagnético foi separado em primário e secundário, e então, o campo secundário foi calculado pelo método de elementos finitos no domínio (kx, ky, z) da transformada de Fourier. Para estruturas bidimensionais, foi aplicada a mesma metodologia usado nos modelos TII unidimensionais, onde o campo secundário foi calculado pelo método de elementos finitos no domínio (x, ky, z), da transformada de Fourier, com a utilização de malhas não estruturadas para discretização dos modelos. Estas respostas foram usados para avaliar os efeitos da anisotropia elétrica nos dados CSEM marinho 1D e 2,5D. / In this work I present the solution to the electromagnetic field generated by a horizontal electric dipole in transversally isotropic media with vertical (TIV) as well as inclined (TII) symmetry axis. In one-dimensional models the electromagnetic field was obtained with two distinct methods: (1) For the TIV case, I have written a semi-analytical solution to the Maxwell’s equations, by using a vector potential formulation; (2) For the TII case, the field was represented as the composition of primary and secondary fields, where primary fields are those found in an underlying isotropic layered medium, and the secondary field is calculated numerically via the finite element method in the spatial Fourier transform domain (kx, ky, z). This last methodology was also used to calculate the fields in two-dimensional structures, including inclined anisotropy in any region of the models. In this 2,5D case, I have applied the finite element method in the (x, ky, z) do main. Here I have used unstructured meshes to discretise the media, and parallel programming to solve the linear systems of equations. The responses were used to study the effects of electrical anisotropy in marine CSEM data.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-11T13:07:43Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoInfluenciaEletrojato.pdf: 12249813 bytes, checksum: 25879a0f7f5bef64fad1033f047a5a8d (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:22:51Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-27T16:33:10Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoInfluenciaEletrojato.pdf: 12249813 bytes, checksum: 25879a0f7f5bef64fad1033f047a5a8d (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T16:36:54Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoInfluenciaEletrojato.pdf: 12249813 bytes, checksum: 25879a0f7f5bef64fad1033f047a5a8d (MD5) / Made available in DSpace on 2014-09-10T16:36:54Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoInfluenciaEletrojato.pdf: 12249813 bytes, checksum: 25879a0f7f5bef64fad1033f047a5a8d (MD5) Previous issue date: 1991 / O campo magnetotelúrico em regiões equatoriais viola a condição de ondas planas por causa de uma fonte fortemente concentrada na direção E-W na ionosfera, denominada eletrojato equatorial. No presente trabalho, procurou-se analisar a resposta magnetotelúrica de fontes que simulam o efeito do eletrojato equatorial. Foram considerados dois tipos de fontes para simular o eletrojato: uma linha infinita de corrente e uma distribuição gaussiana de densidade de corrente em relação a uma das coordenadas horizontais. A resistividade aparente foi obtida através da relação de Cagniard e comparada com os resultados de ondas planas. É mostrada também a comparação entre a fase da impedância na superfície, para os três tipos de fontes (ondas planas, eletrojato gaussiano e linha de corrente). O problema de meios com heterogeneidades laterais foi resolvido em termos de campos secundários, sendo as equações diferenciais solucionadas através da técnica de elementos finitos bidimensionais. Os resultados mostram que o eletrojato tem pouca influência nas respostas (resistividade aparente e fase) de estruturas geológicas rasas. Entretanto, a influência pode ser considerável nas estruturas profundas (maior que 5000 m), principalmente se suas resistividades são altas (maior que 100 Ω.m). Portanto, a influência do eletrojato equatorial deve ser considerada na interpretação de dados magnetotelúricos de bacias sedimentares profundas ou no estudo da crosta terrestre. / The magnetotelluric field in equatorial regions violates the plane wave condition due to a strongly concentrated electrical source of current in the ionosphere, called equatorial electrojet. In this thesis, it is analized the magnetotelluric reponses of electrical sources in order to simulates the effects of the equatorial electrojet. There were used two types of sources: an infinite current line and a gaussian distribuition of current density along one the horizontal cartisian coordenates. The apparent resistivities were computed using the Cagniard relationship and compared with the plane wave results. It was also carried out a comparison of the phase of the surface impedance for all three types of sources (plane wave, gaussian electrojet an line source). The differential equation of the models with lateral inhomogeneities were solved by the secundary field technique using the finite element method. The results show that the electrojet has minor influence in the MT responces (apparent resistivity and phase) of shallow geological structure. However, the effect can be very large in the case of deep strutures (more than 5000 m), mainly if their resistivities are larger than 100 Ωm. Therefore, the influence of the equatorial electrojet must be taken in consideration in the interpretation of the magnetotelluric data from sedimentary basins and from crustal studies.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-21T16:45:06Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InfluenciaEstruturasGeologicas.pdf: 5440337 bytes, checksum: dbc4ad289b2fcf0d65997eb2f07b3017 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:04:40Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-27T17:00:13Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InfluenciaEstruturasGeologicas.pdf: 5440337 bytes, checksum: dbc4ad289b2fcf0d65997eb2f07b3017 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T16:37:43Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InfluenciaEstruturasGeologicas.pdf: 5440337 bytes, checksum: dbc4ad289b2fcf0d65997eb2f07b3017 (MD5) / Made available in DSpace on 2014-09-10T16:37:43Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InfluenciaEstruturasGeologicas.pdf: 5440337 bytes, checksum: dbc4ad289b2fcf0d65997eb2f07b3017 (MD5) Previous issue date: 2005 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / A Terra atua como um grande magneto esférico, cujo campo assemelha-se àquele gerado por um dipolo magnético. Este campo apresenta mudanças de intensidade que variam com a localização e a hora local. A parte principal do campo geomagnético se origina no interior da Terra através de processos eletromagnéticos. Extensivos estudos mostraram ainda que existem contribuições de origem externa ao planeta, principalmente de origem solar. Dentre estas fontes há anomalias do campo magnético que surgem a partir de um aumento diurno da corrente elétrica em uma estreita faixa da ionosfera, de direção leste-oeste, centrada no equador magnético e denominada Eletrojato Equatorial (EEJ). Ocasionalmente estas correntes podem apresentar reversões de fluxo, sendo denominadas Contra-Eletrojato (CEJ). Vários autores têm estudado os efeitos do EEJ e CEJ sobre as observações geoeletromagnéticas. Eles estão interessados no efeito combinado do EEJ e estruturas geológicas condutivas 1-D e 2-D. Nestes trabalhos a estrutura 2-D sempre se apresentava paralela ao eletrojato, o que é uma hipótese bastante restritiva ao se modelar ambientes geológicos mais realistas, em que corpos bidimensionais podem ter qualquer strike em relação ao EEJ. Neste trabalho apresentamos a solução deste problema sem esta restrição. Assim, mostramos os campos geoeletromagnéticos devidos a estruturas bidimensionais que possuam strike oblíquo em relação ao EEJ, através de perfis dos campos elétrico e magnéticos calculados na superfície e formando direção arbitrária à heterogeneidade condutiva 2-D. Com esta resposta avaliamos ainda qual a influência que estruturas bidimensionais exercem sobre a resposta magnetotelúrica, sob influência do Eletrojato Equatorial. Durante o desenvolvimento deste trabalho, utilizamos o método de elementos finitos, tendo por fonte eletromagnética o EEJ e o CEJ, que por sua vez foram representados por uma combinação de distribuições gaussianas de densidade de corrente. Estas fontes foram decompostas nas direções paralela e perpendicular à estrutura 2-D, resultando nos modos de propagação TE1 e TE2 e TM acoplados, respectivamente. Resolvemos o modo acoplado aplicando uma Transformada de Fourier nas equações de Maxwell e uma Transformada Inversa de Fourier na solução encontrada. De acordo com os experimentos numéricos realizados em um modelo interpretativo da Anomalia Condutiva da Bacia do Parnaíba, formado por uma enorme estrutura de 3000 ohm-m dentro de um corpo externo condutivo (1 ohm-m), concluímos que a presença do CEJ causa uma inversão na anomalia, se compararmos com o resultado do EEJ. Concluímos também que para as frequências mais altas as componentes do campo elétrico apresentam menor influência da parte interna do corpo 2-D do que da parte externa. Já para frequências mais baixas este comportamento se observa com as componentes do campo magnético. Com relação à frequência, vimos os efeitos do “skin-depth”, principalmente nas respostas magnéticas. Além disso, quando a estrutura 2-D está paralela ao eletrojato, o campo elétrico é insensível à estrutura interna do modelo para todos os valores de frequência utilizados. Com respeito ao ângulo θh entre a heterogeneidade e a fonte, vimos que o modo TM se manifesta naturalmente quando θh é diferente de 0°. Neste caso, o modo TE é composto por uma parte devido à componente da fonte paralela à heterogeneidade e a outra devido à componente da fonte perpendicular, que é acoplada ao modo TM. Assim, os campos calculados têm relação direta com o valor de θh. Analisando a influência do ângulo entre a direção do perfil dos campos e o strike da heterogeneidade verificamos que, à medida que θh se aproxima de 90°, os campos primários tornam-se variáveis para valores de θp diferentes de 90°. Estas variações causam uma assimetria na anomalia e dão uma idéia da inclinação da direção do perfil em relação aos corpos. Finalmente, concluímos que uma das influências que a distância entre o centro do EEJ e o centro da estrutura 2-D, causa sobre as componentes dos campos está relacionado às correntes reversas do EEJ e CEJ, pois a 500 km do centro da fonte estas correntes têm máxima intensidade. No entanto, com o aumento da distância, as anomalias diminuem de intensidade. Nas sondagens MT, nós também usamos o EEJ e o CEJ como fonte primária e comparamos nossos resultados com a resposta da onda plana. Deste modo observamos que as componentes do campo geoeletromagnético, usadas para calcular a impedância, têm influência do fator de acoplamento entre os modos TE2 e TM. Além disso, esta influência se torna maior em meios resistivos e nas frequências mais baixas. No entanto, o fator de acoplamento não afeta os dados magnetotelúricos em frequências maiores de 10-2 Hz. Para frequências da ordem de 10-4 Hz os dados MT apresentam duas fontes de perturbação: a primeira e mais evidente é devido à presença fonte 2-D (EEJ e CEJ), que viola a hipótese da onda plana no método MT; e a segunda é causada pelo acoplamento entre os modos TE2 e TM, pois quando a estrutura bidimensional está obliqua à fonte 2-D temos correntes elétricas adicionais ao longo da heterogeneidade. Concluimos assim, que o strike de uma grande estrutura condutiva bidimensional relativamente à direção do EEJ ou CEJ tem de fato influência sobre o campo geomagnético. Por outro lado, para estudos magnetotelúricos rasos (frequências maiores de 10-3 Hz) o efeito do ângulo entre a estrutura geológica 2-D e a direção do EEJ não é tão importante. Contudo, em estudos de litosfera frequências menores de 10-3 Hz) o acoplamento entre os modos TE2 e TM não pode ser ignorado. / The Earth acts as a large magnet, whose field resembles one that is generated by a magnetic dipole. This field presents intensity changes that vary with observation location and the local time. The main part of the geomagnetic field is created within the Earth by electromagnetic processes. Extensive studies showed that there are also contributions from outside the Earth, mainly from solar origin. Among these outside sources there are anomalies of the magnetic field that arise from an diurnal increase of the electric current in a narrow strip located in the ionosphere, with east-west direction, centered above the magnetic equator and denominated Equatorial Electrojet (EEJ). Occasionally these currents present flow reversions, therefore denominated Counter-Electrojet (CEJ). Several authors have been studying the effects of the EEJ on the geomagnetic observations. They are interested in the combined effect of the equatorial electrojet and the 1-D e 2-D conductive geological structure underneath. In these works the 2-D structure is always considered parallel to the electrojet, which is a quite restrictive hypotheses in view to realistic geological situation, in that two-dimensional structures can have any direction in relation to the electrojet. We present in this work the solution of this problem without this restriction. Thus, here we present the geomagnetic fields due to a two-dimensional structure that possess oblique strike in relation to Equatorial Electrojet, through profiles of the electric and magnetic fields, calculated on the surface and forming arbitrary direction to the 2-D conductive heterogeneity. Further, we also evaluate the influence of an arbitrarily oriented two-dimensional structure would cause on the Magnetotelluric data, under the quatorial Electrojet. In the development of this work, we applied the method of finite elements with the EEJ and CEJ as electromagnetic source, that was represented by a sum of gaussians distributions of current density. This source was decomposed in the parallel and the perpendicular directions to the 2-D structure, resulting in the mode TE1 and the coupled mode TE2 and TM, respectively. We solved the coupled mode applying a Fourier Transform in the Maxwell equations and one Inverse Fourier Transform in the transformed-domain solution. According to the numerical experiments on a interpretative model of Parnaíba Basin Conductivity Anomaly, formed by a large 3000 ohm-m body inside a very conductive (1ohm-m) external structure, we conclude that the presence of the CEJ causes an inversion in the anomaly. We also conclude that at high frequencies the components of the electric field present smaller influence of the internal part of the 2-D body than the external part. Otherwise, we observed this behavior in the magnetic field at low frequencies. Varying the frequencies, we saw the effects of the “skin-depth” mainly in the magnetic field. Besides, there are situations where electric field is insensitive to the internal structure of the model for all values of the frequency used. With regard to the angle θh between the strike of the conductive heterogeneity and the EEJ direction, we observe the TM mode naturally when θh is different from 0°. In this case, the TE mode is composed of two components, one decomposition of the EEJ parallel to the heterogeneity and the other perpendicular to it. As consequence, the E and B fields have all their three components. When we analyzed the influence of the angle between the direction of the profile of fields and the strike of the 2-D heterogeneity, we conclude that its variation causes an asymmetry on the anomalies, which give an idea about the profile’s direction. Finally, we conclude that one of the influences that the distance between the center of the electrojet and center of the 2-D structure causes on the fields is related to the reverse currents, because at 500 km from the EEJ’s center, these currents have their maximum intensity. In the MT soundings, we also used the EEJ and CEJ as primary source and we compared our results with the plane wave response. We noted that the components of the geomagnetic field, used to evaluate the impedance, have an influence from the coupling factor between the TE2 and TM modes. Moreover, this influence become greater with decreasing of the frequency and for resistive host. However, the coupling factor do not affects the MT response at frequencies higher than 10-2 Hz. For lower frequencies, about 10-4 Hz, we detect two kinds of pertubations on the MT data with respect to the plane-wave one: the first is due the presence of the 2-D electromagnetic source (EEJ and CEJ) as primary field, which violates the plane wave hypothesis; and the second is caused by the coupled TE and TM modes because additional electric currents arise in the heterogeneity’s direction when it is oblique to EEJ. These efects increase with the resistivity of the environment. In conclusion, the strike of a large conductive 2-D structure relative to the direction of the EEJ or CEJ do have influence on the geomagnetic field. On the other hand, for shallow magnetotelluric studies (frequencies higher than 10-3 Hz) the effect of angle between the strike of the 2-D geological structure and the direction of the EEJ is not so important. However, for litospheric studies (frequencies lower than 10-3 Hz) the coupling between the two modes can not be ignored.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-17T12:28:14Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DelineamentoPeTalude.pdf: 13439786 bytes, checksum: 078ede70ecfa1f19ece5e14771d6405a (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Ausência de palavras-chave on 2014-08-07T12:41:09Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-12T12:23:14Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DelineamentoPeTalude.pdf: 13439786 bytes, checksum: 078ede70ecfa1f19ece5e14771d6405a (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-15T16:02:50Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DelineamentoPeTalude.pdf: 13439786 bytes, checksum: 078ede70ecfa1f19ece5e14771d6405a (MD5) / Made available in DSpace on 2014-09-15T16:02:50Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DelineamentoPeTalude.pdf: 13439786 bytes, checksum: 078ede70ecfa1f19ece5e14771d6405a (MD5) Previous issue date: 1992 / A metodologia de interpretação integrada dos dados geológicos e geofísicos observados em um perfil da margem continental do Ceará possibilita a identificação e a integração de características peculiares a cada tipo de dado. Dessa forma, é possível se definir a localização mais provável de feições estruturais importantes, tais como a fronteira entre as crostas continental e oceânica e o pé do talude, objeto do presente estudo. Segundo o Artigo 76 (parágrafo 4, item b) da Convenção das Nações Unidas para o Direito do Mar, o pé do talude é definido como o ponto de variação máxima do gradiente do talude na sua base. Entretanto, essa definição, apesar de simples no contexto fisiográfico, não é suficiente para propiciar a localização do pé do talude como preconiza a Convenção, razão pela qual se aplicam os métodos geofísicos. Dentro do contexto geofísico-geológico está implícita a interpretação quantitativa das anomalias gravimétricas ar-livre, que possibilita o delineamento do modelo geofísico representando a subsuperfície, cuja finalidade é subsidiar geologicamente a interpretação integrada dos dados aludidos. Um procedimento automático de ajuste de curvas combinando as técnicas de inversão de busca sistemática e a que utiliza derivadas foi usado com o propósito de gerar o modelo geofísico. A aplicação rigorosa de vínculos preliminarmente e a constante reavaliação desses vínculos através de um processo interativo entre a sísmica e a gravimetria, gerado durante a interpretação quantitativa das anomalias ar-livre, possibilitaram que o modelo geofísico final estivesse dentro dos padrões geológicos para área, notadamente quanto ao equilíbrio isostático (Teoria de Airy). O objetivo do presente trabalho é se estudar as características geológicas e geofísicas observadas ao longo de um perfil da margem continental do Ceará (LEPLAC III), notadamente quanto ao pé do talude, buscando estabelecer a aplicabilidade de uma metodologia de interpretação integrada desses dados, cuja finalidade é se definir de forma sistemática a localização mais provável para esta feição fisiográfica. A metodologia de interpretação integrada dos dados geológicos e geofísicos empregada mostrou-se eficiente para este objetivo. Foi possível se integrar: (i) a localização fisiográfica (distância da costa e profundidade) do pé do talude; (ii) a zona de instabilidade tectônica evidenciada pelos falhamentos, comuns nesta região; (iii) ao fim de uma zona magnética perturbada, associado a um ponto de mínimo na curva de anomalia magnética, e que possivelmente delimita o início de uma zona magnética quieta, denominada de anomalia E e (iv) um ponto de inflexão na curva de anomalia ar-livre, associado ao efeito gravimétrico do contraste de densidades entre as crostas continental, os sedimentos e a água do mar, evidenciado pela geometria do talude. Foi possível ainda se definir a localização mais provável para a fronteira entre as crostas continental e oceânica. Dada a rigorosidade na aplicação das técnicas de inversão e dos vínculos é provável que as correlações das características intrínsecas a cada tipo de dado efetuadas na conclusão desse trabalho tenham fundamento e possam ser confirmadas. A condição para isto é a aplicação da metodologia aqui estabelecida em um número maior de perfis. / The methodology of integrated interpretation applied to the geological and geophysical data observed along a profile at the continental margin of Ceará, enables the identification and integration of characteristics particular to each kind of datum. In this manner, it is possible to define the most probable location of important structural features, such as the boundary line between the continental and oceanic crust and the foot of the continental slope, which is the subject of this study. According to Article 76 (paragraph 4, item b) of the United Nations Covenant regarding the Law of the Sea, the foot of the continental slope is defined as the point of maximum variation of the gradient of the slope at its base. However this definition, despite its simplicity in a physiographic context, it is not sufficient to define the location of the foot of the slope, according to the Covenant, and for this reason geophysical methods are used. Within the geophysical-geological context, a quantitative interpretation of the free-air anomalies, can define the geophysical model, which represents the sub-surface, and can aid the integrated interpretation of the above mentioned data. An automatic procedure of curve adjustment, combining the inversion techniques of systematic search and gradients, was used to generate the geophysical model. The previous rigorous application of constraints and the constant re-evaluation of these constraints by means of an interactive process between seismics and gravimetrics, generated during the quantitative interpretation of free-air anomalies, constrained the final geophysical model to be within the geological framework of the area and within Airy's theory of isostatic equilibrium. The purpose of this research is to study the geological and geophysical characteristics observed on a profile at the continental margin of Ceará (LEPLAC III), especially at the foot of the continental slope, and to try to establish a methodology of integrated interpretation of these data, whose objective is to define in a systematic way, the most probable location for this physiographic feature. The methodology used turned out to be very efficient for the location of the foot of continental slope. In this sense, it was possible to integrate: (i) its physiographic location (distance from the coast and the water depth); (ii) a zone of tectonic instability inferred by the faulting, very common in the continental slope; (iii) the end of a disturbed magnetic anomaly zone, and that possibly delimits the beginning of magnetic quite zone, named anomaly E and (iv) a point of inflexion in the curve or free-air anomaly, associated to an gravimetric effect of contrast in the densities of the continental crusts, the sediment and the sea water. It was also possible to define the most probable location for the boundary line between the continental and oceanic crusts. Because of the rigorous application of the inversion techniques and constraints used, it is possible that the correlation of the intrinsic characteristics for each type of datum, performed at the conclusion of this research, have some foundation and can be confirmed if the methodology described is applied in a greater number of profiles.

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

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-10T13:10:33Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AproximacoesHiperbolicasTempos.pdf: 4094530 bytes, checksum: 2d1b416e88ea1617317477592ae645d5 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-06T12:35:22Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-10T12:19:23Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AproximacoesHiperbolicasTempos.pdf: 4094530 bytes, checksum: 2d1b416e88ea1617317477592ae645d5 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-18T11:33:11Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AproximacoesHiperbolicasTempos.pdf: 4094530 bytes, checksum: 2d1b416e88ea1617317477592ae645d5 (MD5) / Made available in DSpace on 2014-09-18T11:33:11Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AproximacoesHiperbolicasTempos.pdf: 4094530 bytes, checksum: 2d1b416e88ea1617317477592ae645d5 (MD5) Previous issue date: 2005 / A simulação de uma seção sísmica de afastamento nulo (AN) 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 melhor a razão sinal/ruído. Baseado na aproximação hiperbólica dos tempos de trânsito dependente de três parâmetros ou atributos cinemáticos de frentes de onda, recentemente, vem desenvolvendo-se um novo método para simular seções (AN) chamado método de empilhamento sísmico por Superfície de Reflexão Comum (ou empilhamento SRC). Também, seguindo este novo conceito de imageamento sísmico está surgindo um método para simular seções com afastamento comum (AC) a partir de dados de cobertura múltipla usando aproximações dos tempos de trânsito paraxiais na vizinhança de um raio central com afastamento finito. Esta nova aproximação dos tempos de trânsito depende de cinco atributos cinemáticos. Neste trabalho, a partir da aproximação dos tempos de trânsito paraxiais em relação a um raio central com afastamento finito, derivamos uma nova equação do tempo de trânsito usando a condição de um ponto difrator em profundidade, reduzindo a equação original para quatro parâmetros. Para ambas aproximações (reflexão e difração), mostramos a superfície de empilhamento SRC com afastamento finito (SRC-AF). Considerando um modelo sintético, realizamos um estudo comparativo das aproximações dos tempos de trânsito para as quatro configurações sísmicas (fonte comum (FC), receptor comum (RC), ponto-médio-comum (PMC) e afastamento comum (AC)). Para analisar o comportamento do operador SRC-AF, quando este é perturbado, discutimos sua sensibilidade em relação a cada um dos cinco parâmetros (K₁, K₂, K₃, β_S e β_G). Esta análise de sensibilidade é realizada em duas formas: Sensibilidade através da primeira derivada e Sensibilidade no Empilhamento SRC-AF. Após realizar a análise de sensibilidade utilizamos uma nova condição, K₂ = 0 e assim, obtemos uma nova aproximação, agora dependente de três parâmetros. Usando essas aproximações hiperbólicas (em função de cinco, quatro e três parâmetros), propomos um algoritmo para a simulação de seções AC a partir de dados de cobertura múltipla. Finalmente, é apresentado um estudo da zona de Fresnel, com o objetivo de determinar a delimitação da abertura da superfície de empilhamento SRC-AF. / The simulation of a zero-offset seismic (ZO) section from multicoverage, through of stack, is a very used method of seismic reflection imaging, that reduce the amount of data and is aimed to improve the signal/noise ratio. Based on hyperbolic traveltime approximation depending on three kinematic attributes, recently, it was developed a new method to simulate zero-offset (ZO) sections called common reflection surface stack (CRS) method. Also, following this new concept of seismic imaging it was introduzed a method to simulate common-offset (CO) sections from multicoverage data by using a paraxial traveltime approximation for paraxial rays in the vicinity of the finite-offset centra ray. This new traveltime approximation depend on five kinematic attributes. In this work, from the hyperbolic paraxial traveltime approximation with finite-offset central ray, we derive a new formula of traveltime approximation by using the diffraction point condition, reducing the original formula to four parameters. For both approximations, that´s for reflection and diffractions, we show the finite-offset FO-CRS stack operators. We, also obtain these formulas to the four seismic configurations (common-shot (CS), common-receiver (CR), common-midpoint (CMP) and common-offset (CO)). To analise behaviour FO-CRS stack, when this is perturbed, we debatable the sensibility with respect to each one the five parameters (K₁, K₂, K₃, β_S and β_G). This sensibility analysis was perfomed in two ways: Sensibility through the disturbance of each parameter visualized in the stacking surfaces FO-CRS- and FO-CDS, and the first derivative of the traveltimes FO-CRS and FO-CDS. After realize the sensibility analysis, we utility a new condition, K₂ = 0 and therefore, we obtained a new approximation, now depend on three parameters. Using that hyperbolic traveltime approximations (in function of five, four and three parameters), we propose an algorithm multicoverage seismic reflection data. Finally, is presented a studed of Fresnel zone, to delimit the aperture of FO-CRS stacking surface.