Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-22T14:44:50Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InversaoVinculadaDados.pdf: 34199418 bytes, checksum: fccd8f42ebfa3799899d9531b88a0cb8 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:13:16Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-29T17:16:12Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InversaoVinculadaDados.pdf: 34199418 bytes, checksum: fccd8f42ebfa3799899d9531b88a0cb8 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T12:55:01Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InversaoVinculadaDados.pdf: 34199418 bytes, checksum: fccd8f42ebfa3799899d9531b88a0cb8 (MD5) / Made available in DSpace on 2014-09-10T12:55:01Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InversaoVinculadaDados.pdf: 34199418 bytes, checksum: fccd8f42ebfa3799899d9531b88a0cb8 (MD5) Previous issue date: 2000 / Uma técnica para a inversão de dados magnetotelúricos é apresentada neste trabalho. Dois tipos de dados são tratados aqui, dados gerados por modelos unidimensionais com anisotropia na condutividade das camadas e dados bi-dimensionais de levantamentos do método EMAP (ElectroMagnetic Array Profiling). Em ambos os casos fazemos a inversão usando vínculos aproximados de igualdade para estabilizar as soluções. Mostramos as vantagens e as limitações do uso destes vínculos nos processos de inversão. Mesmo vinculada a inversão ainda pode se tornar instável. Para inverter os dados 2-D do EMAP, apresentamos um processo que consiste de três partes: 1 – A construção de um modelo interpretativo e da aproximação inicial para a inversão a partir dos dados de seções de resistividade aparente filtradas pelo processo de filtragem do EMAP; 2 – a inclusão de uma camada de corpos pequenos aflorantes, chamada de camada destatic shift, aos modelos interpretativos para resolver as fontes de distorções estáticas que contaminam os dados; 3 – o uso dos vínculos aproximados de igualdade absoluta para estabilizar as soluções. Os dois primeiros passos nos permitem extrair o máximo de informação possível dos dados, enquanto que o uso dos vínculos de igualdade nos permite incluir informação a priori que possua significado físico e geológico. Com estes passos, obtemos uma solução estável e significativa. Estudaremos o método em dados sintéticos de modelos bi-dimensionais e em dados reais de uma linha de EMAP feita na Bacia do Paraná. / We present a technique to invert magnetotelluric data. Two kinds of MT data are treated: data generated by anisotropic one-dimensional models and two-dimensional data from EMAP (ElectroMagnetic Array Profiling) surveys. In both cases we perform the inversions using approximate equality constraints to stabilize the solutions. We show the advantages and the limitations of those constraints in the inversion process. Even when the inversion is constrained, it can still be unstable. To invert the two-dimensional EMAP data we present a process that consists of three steps: 1 – The interpretation model and the initial guess for each case are derived from the filtered apparent resistivity sections which are the result of the filtering process of the EMAP method; 2 – One layer formed by small outcropping bodies (termed static shift layer) is added to the interpretation model in order to resolve the sources of static distortions that contaminate the data; 3 – The inversion is constrained by approximate equality constraints. The first two steps enable us to extract the largest possible amount of information from the data, while the use of the approximate equality constraints allow us to incorporate a priori information which has true physical or geological meaning. With these steps we obtain solutions that are stable and meaningful. We study the method in two-dimensional models as well as in real data from an EMAP line surveyed at the Paraná Basin.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-22T13:48:05Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InversaoDadosOndas.pdf: 140994516 bytes, checksum: 401ef5e70ef9a780409ce40a7c36d018 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:22:02Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-01T16:35:48Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InversaoDadosOndas.pdf: 140994516 bytes, checksum: 401ef5e70ef9a780409ce40a7c36d018 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T12:55:18Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InversaoDadosOndas.pdf: 140994516 bytes, checksum: 401ef5e70ef9a780409ce40a7c36d018 (MD5) / Made available in DSpace on 2014-09-10T12:55:18Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InversaoDadosOndas.pdf: 140994516 bytes, checksum: 401ef5e70ef9a780409ce40a7c36d018 (MD5) Previous issue date: 1996 / Neste trabalho compilamos informações sobre um grande número de medidas de velocidade de grupo para ondas Rayleigh do modo fundamental, com período até 100 segundos. Tais dados consistiram de informações retiradas da literatura geofísica e cobriram toda a Terra. Parte dos dados foi organizada em trabalhos anteriores e uma segunda parte foi apresentada aqui de forma inédita. Para a América do Sul, selecionamos os principais conjuntos de dados de tais ondas e elaboramos diversos perfis onde a distribuição de velocidade de ondas cisalhantes foi obtida a partir da inversão das curvas de dispersão de velocidade de grupo. Tais perfis serviram para termos uma ideia inicial da estrutura interna da Terra em nosso continente. Com o conjunto global de dados de velocidade de grupo foi possível obtermos os mapas de distribuição lateral de valores de velocidade para cada período referencial entre 20 e 100 segundos. Tais mapas foram produzidos da mesma forma que os mapas de velocidade de fase de ROSA (1986), onde a amostragem for para realizada para blocas medindo 10x10 graus, englobando toda a Terra, em projeção mercator. O valor de velocidade de grupo em cada bloco, para cada período, foi obtido a partir da inversão estocástica dos dados de anomalia de velocidade em relação aos modelos regionalizados de JORDAN (1981) com os valores de velocidade de grupo de ROSA et al. (1992). Os mapas de velocidade de grupo obtidos aqui foram então empregados, na América do Sul, com os valores de velocidade de fase dos mapas obtidos por ROSA (1986). Assim, foi possível determinarmos, em profundidade, os mapas de variação de velocidade de onda cisalhante e os mapas de distribuição de valores de densidade. Com isto, pudemos construir o primeiro mapa de profundidade do Moho (todo do Manto Superior) da América do Sul. / We have compiled a large set of group velocity measurements for fundamental mode Rayleigh waves with period smaller than 100 seconds. These data are composed of worldwide information gathered from the geophysical literature. Part of the data was compiled in previous works, and a second part of these data was measured in this thesis. For the South American plate, we have selected the main sets of measurements for such waves and gathered them in several profiles, for which the depth distribution of the shear wave velocity was determined from the inversion of the group velocity dispersion curves. These depth profiles were useful to have an idea of the internal structure of the Earth underneath the South American plate. Using the global group velocity data set, it was possible to determine the maps showing the lateral variation of group velocity for reference period values ranging from 20 to 100 seconds. These maps were produced in the same way that ROSA (1986) used to generate the corresponding phase velocity maps, using blocks measuring 10x10 degrees around the earth and considering the mercator projection. The group velocity value at each block was obtained, for each reference period value, from the stochastic inversion of the travel time delays from the JORDAN (1981) each model, considering the average group velocity values determined by ROSA et al. (1992) for this model. The group velocity maps obtained here were then jointly used with the phase velocity values determined by ROSA (1986) for the determination of the earth’s internal distribution of shear wave velocity as well as the density distribution, from the depth inversion of both phase and group velocity dispersion curves. The results enabled us to construct the first map showing the possible depth of the Moho discontinuity (upper mantle depth) for South America.

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

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

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

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

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-30T14:40:43Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_NovaAbordagemInterpretacao.pdf: 15232651 bytes, checksum: 8838ce5885d35952a159a082516bf81f (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:32:59Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-12T13:53:09Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_NovaAbordagemInterpretacao.pdf: 15232651 bytes, checksum: 8838ce5885d35952a159a082516bf81f (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-15T16:03:50Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_NovaAbordagemInterpretacao.pdf: 15232651 bytes, checksum: 8838ce5885d35952a159a082516bf81f (MD5) / Made available in DSpace on 2014-09-15T16:03:50Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_NovaAbordagemInterpretacao.pdf: 15232651 bytes, checksum: 8838ce5885d35952a159a082516bf81f (MD5) Previous issue date: 1989 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / CNRS - Centre National de la Recherche Scientifique / FINEP - Financiadora de Estudos e Projetos / A interpretação de anomalias gravimétricas é de grande importância no estudo de feições geológicas que ocorrem na crosta terrestre. Esta interpretação é, no entanto, dificultada pelo fato das anomalias gravimétricas serem resultantes da soma total dos efeitos produzidos por todos os contrastes de densidades de subsuperfície. Desse modo, com o intuito de separar efeitos de feições mais profundas de efeitos de feições mais rasas, bem como a caracterização da geometria desses dois conjuntos de feições, apresentamos um método de separação das componentes regional e residual do campo e a subsequente interpretação de cada componente. A separação regional-residual de dados gravimétricos é efetuada através da aproximação do campo regional por um polinômio ajustado ao campo observado por um método robusto. Este método é iterativo e usa como aproximação inicial a solução obtida através do ajuste polinomial pelo método dos mínimos quadrados. O método empregado minimiza a influência de observações contendo forte contribuição do campo residual no ajuste do campo regional. A componente regional obtida a partir da separação regional-residual é transformada em um mapa de distâncias verticais em relação a um nível de referência. Esta transformação compreende duas etapas. A primeira consiste na obtenção da continuação para baixo da componente regional, que é pressuposta ser causada por uma interface suave separando dois meios homogêneos, representando a interface crosta-manto, cujo contraste de densidade é supostamente conhecido. A segunda consiste na transformação do mapa de continuação para baixo em um mapa de distâncias verticais entre o nível de continuação (tomado como nível de referência) e a interface. Este método apresenta duas dificuldades. A primeira está ligada à instabilidade, havendo portanto a necessidade do emprego de um estabilizador o que acarreta a perda de resolução das feições que se desejam mapear. A segunda, inerente ao método gravimétrico, consiste na impossibilidade da determinação das profundidades absolutas da interface em cada ponto, bastando entretanto o conhecimento da profundidade absoluta em um ponto, através de informação independente, para que todas as outras profundidades absolutas sejam conhecidas. A componente residual obtida a partir da separação regional-residual é transformada em um mapa de contrastes de densidade aparente. Esta transformação consiste no cálculo do contraste de densidade de várias fontes prismáticas através de uma inversão linear pressupondo que as fontes reais estejam das a uma placa horizontal, com contrastes de densidade variando apenas nas direções horizontais. O desempenho do método de separação regional-residual apresentado foi avaliado, através de testes empregando dados sintéticos, fornecendo resultados superiores em relação aos métodos dos mínimos quadrados e da análise espectral. O método de interpretação da componente regional teve seu desempenho avaliado em testes com dados sintéticos onde foram produzidos mapeamentos de interfaces bem próximas das estruturas reais. O limite de resolução das feições que se desejam mapear depende não só do grau do polinômio ajustante, como também da própria limitação inerente ao método gravimétrico. Na interpretação da componente residual é necessário que se postule ou tenha informação a priori sobre a profundidade do topo e espessura da placa onde as fontes estão supostamente confinadas. No entanto, a aplicação do método em dados sintéticos, produziu estimativas razoáveis para os limites laterais das fontes, mesmo na presença de fontes interferentes, e pressupondo-se valores para profundidade do topo e espessura da placa, diferentes dos valores verdadeiros. A ambiguidade envolvendo profundidade do topo, espessura e densidade pode ser visualizada através de gráficos de valores de densidade aparente contra profundidade do topo presumida para a placa para vários valores postulados para a espessura da placa. Estes mesmos gráficos permitem, pelo aspecto das curvas, a elaboração de uma interpretação semi-quantitativa das profundidades das fontes reais. A seqüência dos três métodos desenvolvidos neste trabalho foi aplicada a dados gravimétricos da região norte do Piauí e noroeste do Ceará levando a um modelo de organização crustal que compreende espessamentos e adelgaçamentos crustais associados a um evento compressivo que possibilitou a colocação de rochas densas da base da crosta a profundidades rasas. Este modelo ê compatível com os dados geológicos de superfície. É ainda sugerida a continuidade, por mais 200 km em direção a sudoeste, do Cinturão de Cisalhamento Noroeste do Ceará por sob os sedimentos da Bacia do Parnaíba, com base nas evidências fornecidas pela interpretação da anomalia residual. Embora esta seqüência de métodos tenha sido desenvolvida com vistas ao estudo de feições crustais de porte continental, ela também pode ser aplicada ao estudo de feições mais localizadas como por exemplo no mapeamento do relevo do embasamento de/bacias sedimentares onde os sedimentos são cortados por rochas intrusivas mais densas. / Despite its great importance to the study of global geologic structures, interpreting gravity anomalies is not a trivial task because the observed gravity field is the resultant of every gravity effect produced by every elementary density contrast. Therefore, in order to isolate the effects produced by shallow sources from those produced by deep sources, I present a new method for regional-residual separation and methods for interpreting each isolated component. The regional-residual separation is perfomed by approximating the regional field by a polynomial fitted to the observed field by a robust method. This method is iterative and its starting value is the least-squares fitting. Also, the influence of observations containing substantial contributions of the residual field in the regional field fitting is minimized. The computed regional field is transformed into a map of vertical distances relative to a given datum. This transformation consists of two stages. The first one is the downward continuation of the regional field which is assumed to be produced by a smooth interface separating two homogeneous media: the crust and the mantle. The density contrast between the media is presumably known. The second stage consists in transforming the downward continued field into a map of vertical distances relative to a given datum by means of simple operations. This method presents two difficulties. The first one is related to the instability inherent to the downward continuation operation. The use of a stabilizer is therefore mandatory, leading to an inevitable loss of resolution of the features being mapped. The second difficulty, inherent to the gravity method, is the impossibility of determining the interface absolute depths. However, the knowledge of the absolute depth at one single point of the interface by independent means allows the computation of all absolute depths. The computed residual component is transformed into an apparent density map. This transformation consists in calculating the intensity of several prismatic sources by linear inversion, assuming that the real sources are confined to a horizontal slab and have density contrasts varying only along the horizontal directions. The performance of the regional-residual separation method was assessed in tests using synthetic data, always producing better results as compared either with polynomial fitting by least-squares or with the spectral analysis method. The method for interpreting the regional component was applied to synthetic data producing interfaces very close to the true ones. The limit of resolution of the features being mapped depend not only on the degree of the fitting polynomial, but also on the limitation imposed by the gravity method itself. In interpreting the residual component, a priori information is needed about the depth and thickness of the slab confining the true sources. However, results of tests using synthetic data showed that reasonable estimates for the h6rizontal limits of the sources can be obtained, even when the depth and thickness of the slab are not known. The ambiguity involving depth to the top, thickness and the apparent density can be visualized by means of curves of apparent density as a function of the presumed depth to the top of the slab, each curve corresponding to a particular assumed value for the slab thickness. An analysis of the configuration of the curves allows a semi-quantitative interpretation of the real sources depths. The sequence of all three methods described above was applied to gravity data from northern Piauí and northwestern Ceará state. As a result, a crustal organization model was obtained consisting of crustal thickenings and thinnings related to a compressive event which caused the raise of dense, lower crust rocks to shallower depths. This model is consistent with surface geological information. Also, the .gravity interpretation suggests the continuity of the Northwestern Ceará Shear Belt for more than 200 km under the Parnaíba Basin sedimentary cover. Although the sequence of methods presented here has been developed for the study of large scale crustal structures, it could also be applied to the interpretation of smaller structures, as, for example, the basement relief of a sedimentary basin where the sediments have been intruded by mafic rocks.

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

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

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