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

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-06T12:01:00Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseSensibilidadeEstereotomografia.pdf: 6412555 bytes, checksum: 4d766be627a360d62d4651c10cc03c7a (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:03:42Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-10T14:12:16Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseSensibilidadeEstereotomografia.pdf: 6412555 bytes, checksum: 4d766be627a360d62d4651c10cc03c7a (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-18T11:31:24Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseSensibilidadeEstereotomografia.pdf: 6412555 bytes, checksum: 4d766be627a360d62d4651c10cc03c7a (MD5) / Made available in DSpace on 2014-09-18T11:31:24Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseSensibilidadeEstereotomografia.pdf: 6412555 bytes, checksum: 4d766be627a360d62d4651c10cc03c7a (MD5) Previous issue date: 2005 / ANP - Agência Nacional do Petróleo / O método de estereotomografia é estendido para meios com anisotropia arbitrária e implementado para meios com anisotropia elíptica e anelíptica. Os modelos elípticos e anelípticos apresentam somente três parâmetros. Isto faz com que eles sejam menos sensíveis a ambiguidade, causada pela cobertura limitada dos raios em experimentos sísmicos de superfície e VSP, do que modelos transversalmente isotrópicos ou ortorrômbicos. As correspondentes aproximações para superfície de vagarosidade limita a validade desta implementação para eventos qP com anisotropia suave. Experimentos numéricos mostram o potencial e as limitações da estereotomografia para estimar macro modelos de velocidade adequados para o imageamento na presença de anisotropia e a importância dos eventos de transmissão de experimentos VSP de multiplo afastamento modelo para o sucesso desta abordagem. / Stereotomography is extended to general anisotropic models and implemented for elliptical and anelliptical anisotropy. The elliptical and anelliptical models present only three parameters. This makes them less sensitive to the ambiguity due to limited coverage of surface seismic experiments than transversaly isotropic or orthorhombic models. The corresponding approximations of the slowness surface restrict the validity of the present approach to qP events and mild anisotropy. Numerical experiments show the potential and the limitations of stereotomography in estimating macro-velocity models suitable for imaging in the presence of anisotropy as well as the importance of transmission events from multiple-offset VSP experiments for the success of the approach.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-21T13:00:33Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoSismicoSuperficie.pdf: 34321354 bytes, checksum: 98bc579e7260095e946f7fa7acc5a244 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:11:52Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-05T13:23:59Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoSismicoSuperficie.pdf: 34321354 bytes, checksum: 98bc579e7260095e946f7fa7acc5a244 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-18T11:35:24Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoSismicoSuperficie.pdf: 34321354 bytes, checksum: 98bc579e7260095e946f7fa7acc5a244 (MD5) / Made available in DSpace on 2014-09-18T11:35:24Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoSismicoSuperficie.pdf: 34321354 bytes, checksum: 98bc579e7260095e946f7fa7acc5a244 (MD5) Previous issue date: 2001 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / ANP - Agência Nacional do Petróleo / O método de empilhamento sísmico por Superfície de Reflexão Comum (ou empilhamento SRC) produz a simulação de seções com afastamento nulo (NA) a partir dos dados de cobertura múltipla. Para meios 2D, o operador de empilhamento SRC depende de três parâmetros que são: o ângulo de emergência do raio central com fonte-receptor nulo (β₀), o raio de curvatura da onda ponto de incidência normal (R_NIP) e o raio de curvatura da onda normal (R_N). O problema crucial para a implementação do método de empilhamento SRC consiste na determinação, a partir dos dados sísmicos, dos três parâmetros ótimos associados a cada ponto de amostragem da seção AN a ser simulada. No presente trabalho foi desenvolvido uma nova sequência de processamento para a simulação de seções AN por meio do método de empilhamento SRC. Neste novo algoritmo, a determinação dos três parâmetros ótimos que definem o operador de empilhamento SRC é realizada em três etapas: na primeira etapa são estimados dois parâmetros (β°₀ e R°_NIP) por meio de uma busca global bidimensional nos dados de cobertura múltipla. Na segunda etapa é usado o valor de β°₀ estimado para determinar-se o terceiro parâmetro (R°_N) através de uma busca global unidimensional na seção AN resultante da primeira etapa. Em ambas etapas as buscas globais são realizadas aplicando o método de otimização Simulated Annealing (SA). Na terceira etapa são determinados os três parâmetros finais (β₀, R_NIP e R_N) através uma busca local tridimensional aplicando o método de otimização Variable Metric (VM) nos dados de cobertura múltipla. Nesta última etapa é usado o trio de parâmetros (β°₀, R°_NIP, R°_N) estimado nas duas etapas anteriores como aproximação inicial. Com o propósito de simular corretamente os eventos com mergulhos conflitantes, este novo algoritmo prevê a determinação de dois trios de parâmetros associados a pontos de amostragem da seção AN onde há intersecção de eventos. Em outras palavras, nos pontos da seção AN onde dois eventos sísmicos se cruzam são determinados dois trios de parâmetros SRC, os quais serão usados conjuntamente na simulação dos eventos com mergulhos conflitantes. Para avaliar a precisão e eficiência do novo algoritmo, este foi aplicado em dados sintéticos de dois modelos: um com interfaces contínuas e outro com uma interface descontinua. As seções AN simuladas têm elevada razão sinal-ruído e mostram uma clara definição dos eventos refletidos e difratados. A comparação das seções AN simuladas com as suas similares obtidas por modelamento direto mostra uma correta simulação de reflexões e difrações. Além disso, a comparação dos valores dos três parâmetros otimizados com os seus correspondentes valores exatos calculados por modelamento direto revela também um alto grau de precisão. Usando a aproximação hiperbólica dos tempos de trânsito, porém sob a condição de R_NIP = R_N, foi desenvolvido um novo algoritmo para a simulação de seções AN contendo predominantemente campos de ondas difratados. De forma similar ao algoritmo de empilhamento SRC, este algoritmo denominado empilhamento por Superfícies de Difração Comum (SDC) também usa os métodos de otimização SA e VM para determinar a dupla de parâmetros ótimos (β₀, R_NIP) que definem o melhor operador de empilhamento SDC. Na primeira etapa utiliza-se o método de otimização SA para determinar os parâmetros iniciais β°₀ e R°_NIP usando o operador de empilhamento com grande abertura. Na segunda etapa, usando os valores estimados de β°₀ e R°_NIP, são melhorados as estimativas do parâmetro R_NIP por meio da aplicação do algoritmo VM na seção AN resultante da primeira etapa. Na terceira etapa são determinados os melhores valores de β°₀ e R°_NIP por meio da aplicação do algoritmo VM nos dados de cobertura múltipla. Vale salientar que a aparente repetição de processos tem como efeito a atenuação progressiva dos eventos refletidos. A aplicação do algoritmo de empilhamento SDC em dados sintéticos contendo campos de ondas refletidos e difratados, produz como resultado principal uma seção AN simulada contendo eventos difratados claramente definidos. Como uma aplicação direta deste resultado na interpretação de dados sísmicos, a migração pós-empilhamento em profundidade da seção AN simulada produz uma seção com a localização correta dos pontos difratores associados às descontinuidades do modelo. / By using an arbitrary source-receiver configuration and without knowledge of the velocity model the recently introduced seismic data stacking method called Common Reflection Surface (CRS) simulates a zero-offset (ZO) section from multi-coverage seismic reflection data. For 2-D acquisition, as by-products provides three normal ray parameters: 1) the emergence angle (β₀); 2) the radius of curvature of the Normal Incidence Point Wave (R_NIP); and 3) the radius of curvature of the Normal Wave (R_N). The CRS stack is based on the hyperbolic traveltime paraxial approximation depending on β₀, R_NIP and R_N. In this thesis is presented a new algorithm of the CRS stack based on two-parameters and one-parameter search strategy combining global and local optimization methods for determine the three parameters that define the stacking surface (or operator). This is performed in three steps: 1) two-parameters search by applying global optimization to determine β₀ and R_NIP; 2) one-parameter search by applying global optimization to determine R_N; and 3) three-parameters search by applying local optimization to determine three parameters, using as initial approximations the parameter triple of the earlier two steps. In the first two steps is used the Simulated Annealing (SA) algorithm and the Variable Metric algorithm is used in the third step. To simulate the conflicting dip events, for each ZO sample where there are interference of intersecting events is determined an additional parameter triple corresponding to a local minimum. The stacking along the respective operator for each particular event allows to simulate their interference in the simulated ZO section by means of their superposition. This new CRS stack algoritm was applied to synthetic data sets providing high-quality simulated ZO sections and high precision determination of the stack parameters in comparison with the forward modeling. Using the hyperbolic traveltime approximation for identical radii of curvature R_NIP = R_N, an algorithm called Common Diffraction Surface (CDS) stack was developed to simulate ZO sections for diffracted waves. In a similar way to the CRS stack procedure, this new algorithm also uses the SA and VM optimization methods to determine the optimal parameter couple (β₀, R_NIP) that define the best CDS operator. The main features of the algorithm are the data normalization, common-offset data, large aperture of the CDS operator and positive search space for R_NIP. The application of the CDS stack algorithm in a synthetic dataset containing reflected and diffracted wavefields provides as main result a simulated ZO section containing diffracted events clearly defined. The post-stack depth migration of the ZO section locates correctly the discontinuities of the second interface.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-05T11:44:12Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AdaptacaoMetodoBott.pdf: 2164475 bytes, checksum: baccb80d9f2f333abfebec4284c9394c (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:04:06Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-04T16:52:16Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AdaptacaoMetodoBott.pdf: 2164475 bytes, checksum: baccb80d9f2f333abfebec4284c9394c (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-22T11:21:30Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AdaptacaoMetodoBott.pdf: 2164475 bytes, checksum: baccb80d9f2f333abfebec4284c9394c (MD5) / Made available in DSpace on 2014-09-22T11:21:30Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AdaptacaoMetodoBott.pdf: 2164475 bytes, checksum: baccb80d9f2f333abfebec4284c9394c (MD5) Previous issue date: 2013 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Apresentamos um novo método para estimar o relevo do embasamento de bacias sedimentares através da extensão analítica da expressão da placa Bouguer para contraste de densidade entre o pacote sedimentar e o embasamento decrescendo monotonicamente com a profundidade de acordo com uma lei hiperbólica. O método requer ruído contido nos dados gravimétricos tenha desvio padrão inferior a 0,01 mGal. As estimativas do relevo do embasamento são obtidas nas posições espaciais localizadas diretamente abaixo de cada observação. A metodologia foi aplicada a dados sintéticos de bacias sedimentares simuladas apresentando relevo do embasamento suave. O método produziu relevo do embasamento estimado próximo do relevo simulado. O método foi também aplicado a três conjuntos de dados reais: Bacia do Recôncavo (Brasil), Graben do Büyük Menderes (Turquia) e Graben de San Jacinto (Estados Unidos). As soluções produzidas pelo método proposto apresentaram estimativas de profundidades máximas em acordo com as informações geológicas disponíveis: 6 km para Bacia do Recôncavo, 1,6 km para Vale do Büyük Menderes e 2,2 km para Graben de San Jacinto. / We present a new method to estimate the basement relief of a sedimentary basin through the analytic extension of the expression of the slab Bouguer for the case of a density contrast between the basement and the sediments decreasing monotonically with depth according to a hyperbolic law. The method requires that the noise in data have a standard deviation smaller than 0.01 mGal. The basement relief estimates are obtained at positions located directly below each observation. The methodology was applied on synthetic data from a simulated sedimentary basin presenting smooth basement relief. The method produced estimated basement topographies close to simulated basement topographies. The method was also applied to three sets of real Bouguer anomalies: Recôncavo Basin (Brazil), Büyük Menderes Graben (Turkey) and San Jacinto Graben (United States). The solutions produced by the proposed method attained maximum basement depth estimates in accordance with the available geological information: 6 km to Reconcavo Basin, 1.6 km to Büyük Menderes Valley, and 2.2 km to San Jacinto Graben.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-05T12:13:44Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseAproximacoesRpp.pdf: 4103377 bytes, checksum: 5e03cda58fdd1cf4609469596ddec5ad (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:07:51Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-05T13:52:00Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseAproximacoesRpp.pdf: 4103377 bytes, checksum: 5e03cda58fdd1cf4609469596ddec5ad (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-22T11:21:42Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseAproximacoesRpp.pdf: 4103377 bytes, checksum: 5e03cda58fdd1cf4609469596ddec5ad (MD5) / Made available in DSpace on 2014-09-22T11:21:42Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseAproximacoesRpp.pdf: 4103377 bytes, checksum: 5e03cda58fdd1cf4609469596ddec5ad (MD5) Previous issue date: 2002 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / ANP - Agência Nacional do Petróleo / FINEP - Financiadora de Estudos e Projetos / Neste trabalho são apresentadas aproximações lineares e quadráticas das equações de Zoeppritz para a obtenção dos coeficientes de reflexão e transmissão de eventos P-P e P-S em função dos ângulos de incidência e da média angular, bem como a análise de inversão linear de AVO, considerando os eventos de reflexão P-P e P-S dissociados e combinados. O uso das chamadas aproximações pseudo-quadráticas foi aplicado para obtenção de aproximações quadráticas apenas para eventos-PP, em torno dos contrastes médios das velocidades de ondas compressionais e cisalhantes e da razão Vs/Vp. Os resultados das aproximações desenvolvidas neste trabalho mostram que as aproximações quadráticas são mais precisas que as lineares, nas duas versões angulares. As comparações entre as aproximações em termos do ângulo de incidência e da média angular mostram que as aproximações quadráticas são equivalentes dentro do limite angular de [0º, 30º]. Por outro lado, as aproximações lineares em função do ângulo de incidência mostram-se mais precisas que as aproximações lineares em função da média angular. Na inversão linear, fez-se análises de sensibilidade e de ambigüidade e observou-se que, nos caso de eventos de reflexão P-P e P-S dissociados, apenas um parâmetro pode ser estimado e que a combinação destes eventos consegue estabilizar a inversão, permitindo a estimativa de dois dentre os parâmetros físicos dos meios (contrastes de impedância, de velocidade de onda P e de módulo de cisalhamento). / This work presents linear and quadratic approximations of the Zoeppritz equations for the derivation of reflection and transmission coefficients from P-P and P-S events as a function of incidence angles and angular average, as well as the linear inversion analysis, AVO, in respect to the disassociated and combined P-P and P-S reflection events. The use of the so-called pseudoquadratic approximations was applied for the derivation of quadratic approximations only for PPevents, around the average contrasts of compressional and shear waves velocities and Vs/Vp ratio. The results of approximations derived in this work show that quadratic approximations were more precise than the linear ones in the two angular versions. The comparisons between these approximations in terms incidence angle and angular average show that the quadratic approximations are equivalent within the angular limit of [0º to 30º]. In the other hand, the linear approximations as a function of incidence angle are more precise than the linear approximations as a function of the angular average. In the linear inversion, sensitivity and ambiguity analyses were carried out and one could see that in the case of disassociated reflection P-P and P-S events, just a parameter can be estimated and the combination of these events can stabilize the inversion permitting the estimation of two physical parameters for the media (impedance, P-wave velocity and shear bulk module contrasts).

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-18T11:54:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoContatoCamadas.pdf: 10720889 bytes, checksum: 3476f811330f5421d14e271a8e3c2428 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:18:59Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-12T12:35:41Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoContatoCamadas.pdf: 10720889 bytes, checksum: 3476f811330f5421d14e271a8e3c2428 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-22T13:15:32Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoContatoCamadas.pdf: 10720889 bytes, checksum: 3476f811330f5421d14e271a8e3c2428 (MD5) / Made available in DSpace on 2014-09-22T13:15:32Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoContatoCamadas.pdf: 10720889 bytes, checksum: 3476f811330f5421d14e271a8e3c2428 (MD5) Previous issue date: 1996 / PETROBRAS - Petróleo Brasileiro S.A. / Dois dos principais objetivos da interpretação petrofísica de perfis são a determinação dos limites entre as camadas geológicas e o contato entre fluidos. Para isto, o perfil de indução possui algumas importantes propriedades: É sensível ao tipo de fluido e a distribuição do mesmo no espaço poroso; e o seu registro pode ser modelado com precisão satisfatória como sendo uma convolução entre a condutividade da formação e a função resposta da ferramenta. A primeira propriedade assegura uma boa caracterização dos reservatórios e, ao mesmo tempo, evidencia os contatos entre fluidos, o que permite um zoneamento básico do perfil de poço. A segunda propriedade decorre da relação quasi-linear entre o perfil de indução e a condutividade da formação, o que torna possível o uso da teoria dos sistemas lineares e, particularmente, o desenho de filtros digitais adaptados à deconvolução do sinal original. A idéia neste trabalho é produzir um algoritmo capaz de identificar os contatos entre as camadas atravessadas pelo poço, a partir da condutividade aparente lida pelo perfil de indução. Para simplificar o problema, o modelo de formação assume uma distribuição plano-paralela de camadas homogêneas. Este modelo corresponde a um perfil retangular para condutividade da formação. Usando o perfil de entrada digitalizado, os pontos de inflexão são obtidos numericamente a partir dos extremos da primeira derivada. Isto gera uma primeira aproximação do perfil real da formação. Este perfil estimado é então convolvido com a função resposta da ferramenta gerando um perfil de condutividade aparente. Uma função custo de mínimos quadrados condicionada é definida em termos da diferença entre a condutividade aparente medida e a estimada. A minimização da função custo fornece a condutividade das camadas. O problema de otimização para encontrar o melhor perfil retangular para os dados de indução é linear nas amplitudes (condutividades das camadas), mas uma estimativa não linear para os contatos entre as camadas. Neste caso as amplitudes são estimadas de forma linear pelos mínimos quadrados mantendo-se fixos os contatos. Em um segundo passo mantem-se fixas as amplitudes e são calculadas pequenas mudanças nos limites entre as camadas usando uma aproximação linearizada. Este processo é interativo obtendo sucessivos refinamentos até que um critério de convergência seja satisfeito. O algoritmo é aplicado em dados sintéticos e reais demonstrando a robustez do método. / Two of the main objectives in petrophysical log interpretation are to determine geologic bed boundaries and fluid contacts. For these, the induction log has several important properties: it is sensitive to fluid type and distribution in the pore space, the induction measurement is accurately modeled as a convolution of formation conductivity and the tool response function. The first property assures good discrimination of the reserves and at the same time delineates oil-water contacts. This information permits a fundamental zoning of the well log. The second property follows from the quasi-linear relationship between the induction log and formation conductivity. Thus it is possible to use linear system theory, and in particular digital filter design to adaptively deconvolve the original signal. The idea in this work is to produce an algorithm capable of identifying contacts between beds traversed by a borehole, given the apparent conductivity read by an induction tool. To simplify the problem, the formation model is assumed to be a distribution of plane-parallel homogeneous beds. This model corresponds to a rectangular formation conductivity profile. Using the digitized input log, inflexion points are obtained numerically as extrema of first derivatives. This generates a first approximation of the real formation profile. This estimated profile is then convolved with the tool response function giving an estimated apparent conductivity log. A conditioned least-mean-square cost function is defined in terms of the difference between measured and estimated apparent conductivity. Minimizing the cost function yields the bed conductivities. The optimization problem of finding the best rectangular profile for induction data is linear for amplitudes (bed conductivities), but non-linear estimation for bed boundaries. In this case amplitudes are estimated by linear least-squares maintaining fixed contacts. A second pass maintains fixed amplitudes and computes small changes in bed boundaries using a linearized approximation. This processes is iterated to obtain successive refinement until a convergence criteria is satisfied. The algorithm is applied on synthetic and real data showing the robustness of the method.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-09T11:42:28Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoDeconvolucaoHomomorfica.pdf: 9295991 bytes, checksum: d8c9b4b1a029f7b919a238cb75a7235d (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar as palavras-chave on 2014-08-07T13:41:27Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-11T17:45:04Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoDeconvolucaoHomomorfica.pdf: 9295991 bytes, checksum: d8c9b4b1a029f7b919a238cb75a7235d (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-22T13:31:41Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoDeconvolucaoHomomorfica.pdf: 9295991 bytes, checksum: d8c9b4b1a029f7b919a238cb75a7235d (MD5) / Made available in DSpace on 2014-09-22T13:31:41Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoDeconvolucaoHomomorfica.pdf: 9295991 bytes, checksum: d8c9b4b1a029f7b919a238cb75a7235d (MD5) Previous issue date: 1998 / Um registro sísmico é frequentemente representado como a convolução de um pulso-fonte com a resposta do meio ao impulso, relacionada ao caminho da propagação. O processo de separação destes dois componentes da convolução é denominado deconvolução. Existe uma variedade de aproximações para o desenvolvimento de uma deconvolução. Uma das mais comuns é o uso da filtragem linear inversa, ou seja, o processamento do sinal composto, através de um filtro linear, cuja resposta de frequência é a recíproca da transformada de Fourier de um dos componentes do sinal. Obviamente, a fim de usarmos a filtragem inversa, tais componentes devem ser conhecidas ou estimadas. Neste trabalho, tratamos da aplicação a sinais sísmicos, de uma técnica de deconvolução não linear, proposta por Oppenheim (1965), a qual utiliza a teoria de uma classe de sistemas não lineares, que satisfazem um princípio generalizado de superposição, denominados de sistemas homomórficos. Tais sistemas são particularmente úteis na separação de sinais que estão combinados através da operação de convolução. O algoritmo da deconvolução homomórfica transforma o processo de convolução em uma superposição aditiva de seus componentes, com o resultado de que partes simples podem ser separadas mais facilmente. Esta classe de técnicas de filtragem representa uma generalização dos problemas de filtragem linear. O presente método oferece a considerável vantagem de que não é necessário fazer qualquer suposição prévia sobre a natureza do pulso sísmico fonte, ou da resposta do meio ao impulso, não requerendo assim, as considerações usuais de que o pulso seja de fase-mínima e que a distribuição dos impulsos seja aleatória, embora a qualidade dos resultados obtidos pela análise homomórfica seja muito sensível à razão sinal/ruído, como demonstrado. / A seismic record is often represented as the convolution of a wavelet with the impulse response relative to the reflection path. The process of separating these two components of the convolution is termed deconvolution. There are a number of approaches for carrying out a deconvolution. One of the most common is the use of linear inverse filtering, that is, processing the composite signal through a linear filter, whose frequency response is the reciprocal of the Fourier transform of one of the signal components. Obviously, in order to use inverse filtering, such components must be known or estimated. In this work, we deal with the application to seismic signals, of a nonlinear deconvolution technique, proposed by Oppenheim (1965), which uses the theory of a class of nonlinear systems, that satisfy a generalized principle of superposition, which are termed homomorphic systems. Such systems are particularly useful in separating signals which have been combined through the convolution operation. The homomorphic deconvolution algorithm transforms the convolutional process into an additive superposition of its components, with the result that the single parts can be separated more easily. This class of filtering techniques represent a generalization of linear filtering problems. This method offers the considerable advantage that no prior assumption about the nature of the seismic wavelet or the impulse response of the reflection path need be made, that is, it does not require the usual assumptions of a minimum-phase wavelet and a random distribution of impulses, although the quality of the results obtained by the homomorphic analysis is very sensible to the signal/noise ratio, as demonstrated.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-09T12:05:47Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoGravimetroPoco.pdf: 9662344 bytes, checksum: fa38b7f45df88e1e53f180b0c1508600 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:16:26Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-11T14:36:39Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoGravimetroPoco.pdf: 9662344 bytes, checksum: fa38b7f45df88e1e53f180b0c1508600 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-22T17:31:38Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoGravimetroPoco.pdf: 9662344 bytes, checksum: fa38b7f45df88e1e53f180b0c1508600 (MD5) / Made available in DSpace on 2014-09-22T17:31:38Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoGravimetroPoco.pdf: 9662344 bytes, checksum: fa38b7f45df88e1e53f180b0c1508600 (MD5) Previous issue date: 1992 / No período de outubro a novembro de 1990, foi realizado na Bacia Sergipe-Alagoas o primeiro levantamento gravimétrico de poço no Brasil, através de um programa com a participação da Universidade Federal do Pará, a PETROBRÁS e o U.S.G.S. (United States Geological Survey). Este levantamento teve como objetivos o teste do equipamento, a resposta dos perfis gravimétricos em relação aos problemas encontrados nas bacias sedimentares brasileiras e a comparação com os dados de densidade obtidos com o perfil de densidade compensada (CDL). Os levantamentos foram realizados em três poços e os dados obtidos passaram por um processo de redução, onde os valores foram transformados para miligals e corrigidos dos efeitos de maré, deriva e terreno, para obter as densidades preliminares. Porém, a região onde foram realizados os levantamentos apresenta além do embasamento muito raso, uma seqüência evaporítica com contraste de densidade grande em relação ao resto do pacote sedimentar, gerando gradientes verticais anômalos. Estes efeitos foram então corrigidos e obtidos os valores finais de densidade. Com base em uma caracterização litológica prévia, os valores de densidade foram então comparados com o perfil CDL. Verificou-se a boa qualidade dos dados, o que permitiu o cálculo de valores de porosidade, mostrando assim possíveis intervalos de interesse à exploração de hidrocarbonetos, tanto no embasamento, reservatório mais importante, quanto no intervalo sedimentar. / The first borehole gravity survey in Brazil was carried out during October and November of 1990 in Sergipe-Alagoas basin through a joint program between the Federal University of Pará, PETROBRAS and United States Geological Survey. The principal purpose of the survey was to test the applicability of this tool in Brazilian sedimentary basins. Another purpose was to compare the borehole gravimeter (BHGM) data with existing open-hole density log (CDL). The survey was conducted in three wells in Carmópolis field and the data were transformed to miligals, corrected for tide, drift and terrain effects and preliminary density values calculated. The region in which the surveys were carried out consists of an evaporite sequence overlying basement rock ocurring relatively near the surface. The large contrast between these units and the rest of sedimentary section results in vertical gravity anomalies. It is necessary to correct the BHGM density data for these anomalies. The resulting BHGM density data were then compared to that obtained from the CDL. Furthermore, lithology parameters defined from well log data were used to calculate porosities from both the BHGM and CDL density data, showing possible new oil intervals that can be tested.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-07-07T16:13:06Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoComparativoTecnicas.pdf: 927419 bytes, checksum: a5cdbec5e2bc326ca41fb9b7a4a7340b (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Item sem assunto. on 2014-08-07T16:27:49Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-19T13:23:47Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoComparativoTecnicas.pdf: 927419 bytes, checksum: a5cdbec5e2bc326ca41fb9b7a4a7340b (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-20T12:40:16Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoComparativoTecnicas.pdf: 927419 bytes, checksum: a5cdbec5e2bc326ca41fb9b7a4a7340b (MD5) / Made available in DSpace on 2014-08-20T12:40:16Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoComparativoTecnicas.pdf: 927419 bytes, checksum: a5cdbec5e2bc326ca41fb9b7a4a7340b (MD5) Previous issue date: 2007 / Os métodos numéricos de Elementos Finitos e Equação Integral são comumente utilizados para investigações eletromagnéticas na Geofísica, e, para essas modelagens é importante saber qual algoritmo é mais rápido num certo modelo geofísico. Neste trabalho são feitas comparações nos resultados de tempo computacional desses dois métodos em modelos bidimensionais com heterogeneidades condutivas num semiespaço resistivo energizados por uma linha infinita de corrente (com 1000Hz de freqüência) e situada na superfície paralelamente ao "strike" das heterogeneidades. Após a validação e otimização dos programas analisamos o comportamento dos tempos de processamento nos modelos de corpos retangulares variandose o tamanho, o número e a inclinação dos corpos. Além disso, investigamos nesses métodos as etapas que demandam maior custo computacional. Em nossos modelos, o método de Elementos Finitos foi mais vantajoso que o de Equação Integral, com exceção na situação de corpos com baixa condutividade ou com geometria inclinada. / The Finite Element and the Integral Equation numerical methods are commonly used in electromagnetic investigations in Geophysics. In those cases it is important to determine which algorithm is fastest for each geophysical model. In this work we compare computer times for two-dimensional modeling with both methods. The models studied here are formed by conductive rectangular bodies in a resistive half-space, energized by an infinite current line located on the surface and parallel to the strike of the inhomogeneities. The currents in all models operate in a frequency of 1000Hz. After validating and optimizing our codes, we analyze the processing times for several models, varying the number of bodies, their sizes and their tilt. We also investigated which parts of the programs demmand the greatest computer times. We found that the Finite Element method performs faster than the Integral Equation method in all models except those with low conductivity bodies or with inclined bodies.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-30T14:15:42Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EfeitoTopografiaSondagens.pdf: 7433935 bytes, checksum: 6339b45dae74232820291f9bd38da8fb (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-06T13:11:12Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-13T13:44:54Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EfeitoTopografiaSondagens.pdf: 7433935 bytes, checksum: 6339b45dae74232820291f9bd38da8fb (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-20T14:12:19Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EfeitoTopografiaSondagens.pdf: 7433935 bytes, checksum: 6339b45dae74232820291f9bd38da8fb (MD5) / Made available in DSpace on 2014-08-20T14:12:19Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EfeitoTopografiaSondagens.pdf: 7433935 bytes, checksum: 6339b45dae74232820291f9bd38da8fb (MD5) Previous issue date: 1990 / O efeito da topografia de um vale, uma colina, um declive e um aclive, em sondagens eletromagnéticas bidimensionais no domínio da freqüência, causa uma variação nos valores da amplitude e da fase da componente H_z devido a um meio homogêneo. A amplitude é menos afetada que a fase. A parede do vale mais próxima da linha de corrente causa uma forte diminuição dos valores da fase, enquanto que a parede do vale mais distante da linha causa um forte aumento. Os efeitos de uma colina são opostos aos do vale. Os efeitos do declive e do aclive num meio homogêneo, são similares, respectivamente, aos observados pelas paredes do vale e da colina mais próximas da linha de corrente. A resposta de um corpo condutivo retangular num meio homogêneo próximo a uma linha de corrente sofre pequenas variações devido a presença de um vale ou de uma colina situada longe da linha de corrente. Porém, se essas feições topográficas estiverem sobre o corpo, elas afetam fortemente a fase e a amplitude da componente H_z e apenas a amplitude, no caso da componente H_x. A resposta transiente não é afetada pela topografia para tempos muito baixos, pois nesse caso se investiga uma finíssima camada da superfície, assim como para tempos muito altos porque a profundidade de investigação é muito grande comparada com a dimensão da topografia. Para os modelos aqui estudados, a maior influência se dá para tempos intermediários, ao redor de 7 ms, ocasionando um retardo do ponto de "cross over" nas curvas de sondagens. / The topographic effects of a valley, a hill, a downward slope and an upward slope, on bidimensional frequency domain eletromagnetic soundings, cause a variation in the amplitude and the phase of H_z component due to a homogeneous medium. The amplitude is less affected than the phase. The wall of the valley nearer to the current line produces a reduction in the phase values while a farther wall causes an increase. The effects of a hill are opposite to that of a valley. The effects of a downward slope and of an upward slope in the homogeneous medium correspond to that of a single wall, nearer to the current line of the valley and of the hill respectively. The response of a conductive rectangular body placed in a homogeneous medium nearer to the current line sufferc small variations due to the presence of a valley or a hill farther away from the current source. However, if these topographic features lie above the body they affect strongly the phase and amplitude of H_z component and only the amplitude in the case of H_x component. The transient response is not affected by topography at early times because only a small portion of earth is sounded; as well as at later times because the depth of investigation is very large compared to the topographic dimensions. For the models studied here, the major influence is at an intermediate time, about 7 ms, causing a retardation in cross-over positions in the sounding curves.