Dissertação (mestrado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Civil e Ambiental, Pós-Graduação em Tecnologia Ambiental e Recursos Hídricos, 2013. / Submitted by Albânia Cézar de Melo (albania@bce.unb.br) on 2014-03-10T16:40:09Z No. of bitstreams: 1 2013_RafaelMachadoMello.pdf: 4535152 bytes, checksum: 8d61f1511b925a5bd2834fe5c45f58be (MD5) / Approved for entry into archive by Guimaraes Jacqueline(jacqueline.guimaraes@bce.unb.br) on 2014-03-19T14:27:31Z (GMT) No. of bitstreams: 1 2013_RafaelMachadoMello.pdf: 4535152 bytes, checksum: 8d61f1511b925a5bd2834fe5c45f58be (MD5) / Made available in DSpace on 2014-03-19T14:27:31Z (GMT). No. of bitstreams: 1 2013_RafaelMachadoMello.pdf: 4535152 bytes, checksum: 8d61f1511b925a5bd2834fe5c45f58be (MD5) / O presente trabalho consistiu na integração de diferentes ferramentas, modelagem matemática, levantamento geofísico, dados de monitoramento de qualidade e nível d´água para avaliar, por meio de um modelo computacional, o fluxo e a propagação dos contaminantes no aquífero freático de uma região do Distrito Federal, denominada Condomínio Vista Bela, onde esgoto doméstico bruto e de efluente industrial foram infiltrados no solo. O efluente industrial é derivado do processamento, lavagem, tratamento e embalagem de carnes de uma indústria charqueadora. O objetivo do trabalho foi simular o caso ocorrido em março de 2009, quando foi identificada a pluma de contaminação que se propagou até atingir poços e cisternas situadas à jusante do ponto de lançamento. Os poços da região onde ocorreu esse lançamento passaram a ser monitorados por meio de coleta mensal, ao longo de 3 anos. A análise de amostras da água subterrânea permitiu delinear a evolução da pluma em direção aos demais poços e a um rio a jusante do ocorrido. Aproveitando 14 cisternas existentes, utilizadas como única fonte de abastecimento de água da região, foi feito um monitoramento do nível piezométrico que subsidiou a modelagem matemática. Em janeiro de 2011, foi realizado um levantamento geofísico, utilizando a técnica de caminhamento elétrico com arranjo dipolo-dipolo, que possibilitou mapear a pluma de contaminação e também identificar contaminação em outro local a montante da empresa. Em poços da região também foram realizados testes de bombeamento para subsidiar a calibração dos parâmetros. Laudos de qualidade de água foram fornecidos pela a ADASA, órgão gestor de Recursos Hídricos do Distrito Federal. Utilizou-se o programa Visual Modflow para realizar as simulações. Cloreto foi adotado como contaminante de interesse por suas características conservativas, portanto excelente indicador para simulação do fluxo da água subterrânea e do transporte do contaminante. Os dados obtidos no monitoramento, o levantamento geofísico e a modelagem, indicam que a pluma alcançou 110m do ponto de lançamento, após 2 anos de injeção do contaminante, e que em 2 anos alcançará o reservatório no córrego Capão do Brejo, afluente do rio Descoberto. ______________________________________________________________________________ ABSTRACT / The present work is an integration of different tools, including mathematical modeling, geophysical surveying, water quality and level monitoring data to evaluate, by a computer model, the flow and contaminant propagation in the groundwater in a region of the Federal District called Bela Vista condominium, where raw domestic sewage and industrial effluents were infiltrated into the soil. The industrial effluent is derived from the processing, cleaning, treating and packaging processes in a meat processing industry. This work aimed to simulate the case occurred in March 2009, when the contamination plume was identified in the groundwater, spreading towards wells located downstream from the discharge point. The wells of the affected area were monitored monthly, during three years and the analysis of groundwater samples allowed to outline the plume’s evolution toward other wells and toward one river downstream the location. Using the existing 14 existing wells, the only water supply source in the area water level monitoring were carried out to support the mathematical modeling. In January 2011, a geophysical survey was performed, using the electrical profiling technique with dipole-dipole, which helped mapping the contamination plume and identification of another contamination upstream the dumping site. Pumping tests were carried out in some existing wells to provide data to subsidize the calibration. Water quality reports were provided by ADASA, the Federal District Water Resources and Sanitation Agency, into carry out the flow and contaminant propagation Visual Modflow program was used to perform the simulations. Chloride was adopted as the contaminant of interest in the study, because it is conservative and therefore it is an excellent indicator to simulate groundwater flow and transport of the contaminant. The study involving monitoring data geophysical survey and modeling indicated that the plume reached 110m from the point of release, after two years of injection of the contaminantand, in two years after the dumping it will reach the at the Capão do Brejo river dam, a Descoberto river tributary.

Orientadores: Ivan de Araujo Simões Filho, Martin Tygel / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Geociencias / Made available in DSpace on 2018-07-21T07:40:10Z (GMT). No. of bitstreams: 1 Rosalba_JoseFernando_M.pdf: 3884757 bytes, checksum: ebfa5591c5e70d83fb6997a0479cdb64 (MD5) Previous issue date: 1995 / Resumo: O interesse pelo estudo da propagação de ondas elásticas em meios anisotrópicos vem crescendo nos últimos anos em decorrência da necessidade de maior resolução nos dados sísmicos voltados para exploração de petróleo. Com o desenvolvimento de detectores para registro de múltiplos componentes do campo de ondas, de fontes sísmicas apropriadas, associadas a novas técnicas de aquisição de dados como walkway VSP e tomografia poço a poço, faz-se necessário a adoção de modelos elásticos anisotrópicos. A especificação e parametrização de um meio heterogêneo anisotrópico é bastante complexa. No caso mais geral, todos os 21 parâmetros elásticos variam com a posição independentemente. O conceito de meios anisotrópicos fatorados (HAF) aplicado às equações das perturbações dos tempos de trânsito, permite uma simplificação considerável da modelagem direta e inversa, trazendo, ainda, a vantagem de poder separar a contribuição das perturbações da anisotropia da contribuição das perturbações da heterogeneidade. Quando consideradas apenas as perturbações de primeira ordem e meio de referência isotrópico, a diferença entre os tempos de trânsito das duas ondas quase-cisalhantes (tempos de birreftingência) é independente da heterogeneidade do meio, sendo função apenas da anisotropia. Demonstra-se que o mesmo ocorre com a diferença entre os tempos da onda quase-compressional e a média das ondas quase-cisalhantes (tempos de P e S combinados). Neste trabalho, são investigados aspectos da parametrização na inversão dos tempos de birreftingência em meios fraturados e o conteúdo de informações que poderiam advir dos procedimentos de inversões adotados, como por exemplo alguns elementos de simetria, a inversão dos componentes elásticos, densidade de fraturas, etc.. É efetuada também uma análise da parametrização dos tempos de P e S combinados, mostrando-se o grande potencial do uso desta equação para a inversão anisotrópica dos 21 parâmetros elásticos, caso a inversão dos tempos fosse bem sucedida / Abstract: Interest in the study of the propagation of elastic waves in anisotropic media has been growing in recent years, as a result of the need for higher resolution in seismic data related to oil exploration. The development of detectors for recording multiple components of the waves field, from appropriate seismic sources, associated to new data acquisition techniques, such as the vertical seismic profile (VSP) walkway and well-to-well tomography, has resulted in the need to use anisotropic elastic models. The specification and parametrization of an anisotropic heterogeneous medium is very complexo In the most general case, alI the 21 elastic parameters vary according to position, independently. The concept of factorized anisotropic inhomogeneous media (FAI) applied to the equations of perturbations of transit times, makes direct and inverse modeling considerably simplified, and also offers the advantage of being able to separate the contribution of the anisotropy perturbations ftom the contribution of the heterogeneity perturbations. If on1y the first arder and isotropic reference medium perturbations are considered, the difference between the transit times of the two quasi-shearing waves (bireftingence times) does not depend on the heterogeneity of the medium, being a function on1y of the anisotropy. It is shown that the same occurs with the difference between the times of the quasi-compressional wave and the average of quasi-shearing waves (P and S times combined). This paper investigates aspects of parametrization in the inversion of the bireftingence times in ftactured media and the contents of information which could be forthcoming ftom the inversion procedures adopted, such as, for example, some elements of symmetry, the inversion of elastic components, density of ftactures, etc. It also analyzes the parametrization ofthe P and S times combined, showing this equation's great potential for anisotropic inversion of the 21 elastic parameters, if the inversion of the times is successful / Mestrado / Geoengenharia de Reservatorios / Mestre em Geociências

Submitted by Lúcia de Fátima Imbiriba de Sousa (lfis@ufpa.br) on 2018-09-17T19:29:55Z No. of bitstreams: 2 license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Tese_ModelagemEletromagnetica25D.pdf: 10641022 bytes, checksum: d953e3a4e0ec828b3262cfb6c21d10e8 (MD5) / Approved for entry into archive by Lúcia de Fátima Imbiriba de Sousa (lfis@ufpa.br) on 2018-09-17T19:31:51Z (GMT) No. of bitstreams: 2 license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Tese_ModelagemEletromagnetica25D.pdf: 10641022 bytes, checksum: d953e3a4e0ec828b3262cfb6c21d10e8 (MD5) / Made available in DSpace on 2018-09-17T19:31:51Z (GMT). No. of bitstreams: 2 license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Tese_ModelagemEletromagnetica25D.pdf: 10641022 bytes, checksum: d953e3a4e0ec828b3262cfb6c21d10e8 (MD5) Previous issue date: 2014-10-23 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / ANP - Agência Nacional do Petróleo / Apresentamos a formulação eletromagnética em geometria 2.5-D aplicada à modelagem do marine controlled-source electromagnetic (mCSEM) através do método de Diferenças Finitas. Utilizamos a separação dos sinais primário e secundário para evitar problemas de singularidade devido à característica pontual da fonte eletro-magnética, o dipolo elétrico. As componentes do campo eletromagnético são derivadas dos resultados obtidos para os potenciais vetor magnético e escalar elétrico, calculados em todo o domínio do problema, o qual deve ser completamente discretizado para o uso do método de Diferenças Finitas. A limitação imposta pelo uso de malhas estruturadas no delineamento das geometrias presentes nos modelos geológicos, serviu como motivação para introduzirmos o uso de malhas não-estruturadas em nossos problemas. Essas malhas são completamente adaptáveis aos modelos com que trabalhamos, promovendo um delineamento suave de suas estruturas, e podendo ser localmente refinadas apenas nas regiões de interesse. Apresentamos também o desenvolvimento do método RBF-DQ, que faz uso da técnica de aproximação de funções por meio de combinações lineares das funções de base radial (RBF) e da técnica de quadraturas diferenciais (DQ) para a aproximação das derivadas de nosso problema diferencial. Nossos resultados mostraram que o uso do método de Diferenças Finitas com malhas-não estruturadas pode ser aplicado nos problemas de modelagem geofísica, promovendo uma melhoria na qualidade dos dados modelados quando comparados com os resultados obtidos através das técnicas tradicionais de Diferenças Finitas. / We present a 2.5D electromagnetic formulation for modelling of the marine controlledsource electromagnetic (mCSEM) using a Finite Diference frequency domain (FDFD) method. The formulation is in terms of secondary fields thus removing the source point singularities. The components of the electromagnetic field are derived from the solution of the magnetic vector potential and electric scalar potential, evaluated in the entire problem domain that must be completely discretized for the use of the FDFD. Finite difference methods result in large sparse matrix equations that are efficiently solved by sparse matrix algebra preconditioned iterative methods. To overcome limitations imposed by structured grids in the traditional FDFD method, the new method is based upon unstructured grids allowing a better delineation of the geometries. These meshes are completely adaptable to the models we work with, promoting a smooth design of their structures, and may only be refined locally in regions of interest. We also present the development of RBF-DQ method, (radial basis function differential quadrature) which makes use of the technique of functions approximation by linear combinations of radial basis functions (RBF) and the technique of differential quadrature (DQ) for approximation of the derivatives. Our results show that the FDFD method with unstructured grids when applied to geophysical modeling problems, yield improved quality of modeled data in comparison with the results obtained by traditional techniques of FDFD method.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-26T12:54:11Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MapeamentoRelevoEmbasamento.pdf: 11738916 bytes, checksum: bd08738ac99694bdec77ea1cc8ca6a2c (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-06-27T12:30:47Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MapeamentoRelevoEmbasamento.pdf: 11738916 bytes, checksum: bd08738ac99694bdec77ea1cc8ca6a2c (MD5) / Made available in DSpace on 2014-06-27T12:30:47Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MapeamentoRelevoEmbasamento.pdf: 11738916 bytes, checksum: bd08738ac99694bdec77ea1cc8ca6a2c (MD5) Previous issue date: 1998 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Apresentamos três novos métodos estáveis de inversão gravimétrica para estimar o relevo de uma interface arbitrária separando dois meios. Para a garantia da estabilidade da solução, introduzimos informações a priori sobre a interface a ser mapeada, através da minimização de um (ou mais) funcional estabilizante. Portanto, estes três métodos se diferenciam pelos tipos de informação físico-geológica incorporados. No primeiro método, denominado suavidade global, as profundidades da interface são estimadas em pontos discretos, presumindo-se o conhecimento a priori sobre o contraste de densidade entre os meios. Para a estabilização do problema inverso introduzimos dois vínculos: (a) proximidade entre as profundidades estimadas e verdadeiras da interface em alguns pontos fornecidas por furos de sondagem; e (b) proximidade entre as profundidades estimadas em pontos adjacentes. A combinação destes dois vínculos impõe uma suavidade uniforme a toda interface estimada, minimizando, simultaneamente em alguns pontos, os desajustes entre as profundidades conhecidas pelas sondagens e as estimadas nos mesmos pontos. O segundo método, denominado suavidade ponderada, estima as profundidades da interface em pontos discretos, admitindo o conhecimento a priori do contraste de densidade. Neste método, incorpora-se a informação geológica que a interface é suave, exceto em regiões de descontinuidades produzidas por falhas, ou seja, a interface é predominantemente suave porém localmente descontínua. Para a incorporação desta informação, desenvolvemos um processo iterativo em que três tipos de vínculos são impostos aos parâmetros: (a) ponderação da proximidade entre as profundidades estimadas em pontos adjacentes; (b) limites inferior e superior para as profundidades; e (c) proximidade entre todas as profundidades estimadas e um valor numérico conhecido. Inicializando com a solução estimada pelo método da suavidade global, este segundo método, iterativamente, acentua as feições geométricas presentes na solução inicial; ou seja, regiões suaves da interface tendem a tornar-se mais suaves e regiões abruptas tendem a tornar-se mais abruptas. Para tanto, este método atribui diferentes pesos ao vínculo de proximidade entre as profundidades adjacentes. Estes pesos são automaticamente atualizados de modo a acentuar as descontinuidades sutilmente detectadas pela solução da suavidade global. Os vínculos (b) e (c) são usados para compensar a perda da estabilidade, devida à introdução de pesos próximos a zero em alguns dos vínculos de proximidade entre parâmetros adjacentes, e incorporar a informação a priori que a região mais profunda da interface apresenta-se plana e horizontal. O vínculo (b) impõe, de modo estrito, que qualquer profundidade estimada é não negativa e menor que o valor de máxima profundidade da interface conhecido a priori; o vínculo (c) impõe que todas as profundidades estimadas são próximas a um valor que deliberadamente viola a profundidade máxima da interface. O compromisso entre os vínculos conflitantes (b) e (c) resulta na tendenciosidade da solução final em acentuar descontinuidades verticais e apresentar uma estimativa suave e achatada da região mais profunda. O terceiro método, denominado mínimo momento de inércia, estima os contrastes de densidade de uma região da subsuperfície discretizada em volumes elementares prismáticos. Este método incorpora a informação geológica que a interface a ser mapeada delimita uma fonte anômala que apresenta dimensões horizontais maiores que sua maior dimensão vertical, com bordas mergulhando verticalmente ou em direção ao centro de massa e que toda a massa (ou deficiência de massa) anômala está concentrada, de modo compacto, em torno de um nível de referência. Conceitualmente, estas informações são introduzidas pela minimização do momento de inércia das fontes em relação ao nível de referência conhecido a priori. Esta minimização é efetuada em um subespaço de parâmetros consistindo de fontes compactas e apresentando bordas mergulhando verticalmente ou em direção ao centro de massa. Efetivamente, estas informações são introduzidas através de um processo iterativo inicializando com uma solução cujo momento de inércia é próximo a zero, acrescentando, em cada iteração, uma contribuição com mínimo momento de inércia em relação ao nível de referência, de modo que a nova estimativa obedeça a limites mínimo e máximo do contraste de densidade, e minimize, simultaneamente, os desajustes entre os dados gravimétricos observados e ajustados. Adicionalmente, o processo iterativo tende a "congelar" as estimativas em um dos limites (mínimo ou máximo). O resultado final é uma fonte anômala compactada em torno do nível de referência cuja distribuição de constraste de densidade tende ao limite superior (em valor absoluto) estabelecido a priori. Estes três métodos foram aplicados a dados sintéticos e reais produzidos pelo relevo do embasamento de bacias sedimentares. A suavidade global produziu uma boa reconstrução do arcabouço de bacias que violam a condição de suavidade, tanto em dados sintéticos como em dados da Bacia do Recôncavo. Este método, apresenta a menor resolução quando comparado com os outros dois métodos. A suavidade ponderada produziu uma melhoria na resolução de relevos de embasamentos que apresentam falhamentos com grandes rejeitos e altos ângulos de mergulho, indicando uma grande potencialidade na interpretação do arcabouço de bacias extensionais, como mostramos em testes com dados sintéticos e dados do Steptoe Valley, Nevada, EUA, e da Bacia do Recôncavo. No método do mínimo momento de inércia, tomou-se como nível de referência o nível médio do terreno. As aplicações a dados sintéticos e às anomalias Bouguer do Graben de San Jacinto, California, EUA, e da Bacia do Recôncavo mostraram que, em comparação com os métodos da suavidade global e ponderada, este método estima com excelente resolução falhamentos com pequenos rejeitos sem impor a restrição da interface apresentar poucas descontinuidades locais, como no método da suavidade ponderada. / We present three new stable gravity inversion methods to estimate the relief of an interface separating two media. Solution stability is attained by introducing a priori information about the interface, through the minimization of one (or more) stabilizing functional. These methods are, therefore, characterized by the physical and geological information incorporated to the problem. The first method, named global smoothness, estimates the depths to the interface at discrete points by assuming that the density contrast between the media is known. To stabilize the inverse problem, we introduce two different constraints: (a) proximity between the true and estimated interface depths at a few isolated points, and (b) proximity between the estimated depths at adjacent points. The combination of these two constraints impose a uniform degree of smoothness all over the estimated interface, minimizing, simultaneously, the misfit between the known and estimated depths at a few boreholes, for example. The second method, named weighted smoothness, estimates the interface depths at discrete points, assuming that the density contrast is known a priori. In this method, it is incorporated the information that the interface is smooth almost everywhere, but at a few fault discontinuities. To incorporate this attribute to the estimated relief, we developed an iterative process where three kinds of constraints are imposed on parameters: (a) weighted smoothness between values of adjacent parameters, (b) lower and upper bounds on the estimated depths, and (c) proximity between the values of the parameters and a known numerical value. Starting with an initial solution produced by the global smoothness method, this method enhances initially estimated geometric features of the interface; that is, flat areas will tend to become flatter and steep areas will tend to become steeper. This is accomplished by weighting the constraints which require proximity between adjacent parameters. The weights are updated at each iteration so as to enhance the discontinuities detected in a subtle way by the global smoothness method. Constraints (b) and (c) are used both to compensate for the decrease in solution stability due to the introduction of small weights, and to reinforce flatness at the basin bottom. Constraint (b) imposes that any depth be nonnegative and smaller than an a priori known maximum depth value whereas constraint (c) imposes that all depths be closest to a value deliberately violating the maximum depth. The trade-off between these conflicting constraints is attained with a final relief presenting fiat bottom and steep borders. The third method, named minimum moment of inertia, estimates the density contrasts of a subsurface region discretized into elementary prismatic cells. It incorporates the geological information that the interface to be mapped encompasses an anomalous source which besides presenting horizontal extents much larger than its largest vertical extent, exhibits bordes dipping either vertically or toward the center of mass, and that most of the anomalous mass (or mass deficiency) is concentrated, in a compact way, about a reference level. Conceptually, these information are introduced through the minimization of the moment of inertia of the anomalous sources with respect to a reference level coinciding with the mean topographic surface. This minimization is performed in a subspace of parameters consisting of compact sources and presenting bordes which dip either vertically or toward the ce4ter of mass. Effectivelly, these informations are introduced by means of an iterative process starting with a tentative solution dose to the null solution, and adds, at each iteration, a contribution which has minimum moment of inertia with respect to the reference level, in such a way that the estimate of the next iteration does not violate the bounds on the density contrast and minimizes, at the same time, the misfit between the observed and the fitted data. Additionally, the iterative process "freezes" a density estimate if it becomes very dose to either bound. The final solution at the end of the iterative process is an estimated solution exhibiting a compact mass distribution concentrated about the reference level, whose density contrast distribution is dose to the upper (in absolute value) bound established a priori. All three methods were applied to synthetic and field gravity data, produced, respectively, by simulated and real sedimentary basins. The global smoothness method produced a good reconstruction of the basin structural framework even when the true basements were not globally smooth, as was the case of the Recôncavo Basin, Brazil. This method presents, however, the lowest resolution as compared with the other two methods. The weighted smoothness method improved the resolution of basements presenting disontinuities produced by gravity faults with large vertical offsets. It is, therefore, potentially useful in interpreting the structural framework of extensional basins as illustrated both with synthetic data and data from the Steptoe Valley, Nevada, USA and from Recôncavo Basin, Brazil. The minimum moment of inertia method was also applied to synthetic data and data from Recôncavo Basin and from San Jacinto Graben, California, USA. The results showed that, as compared with the other two methods, this method produces excellent estimates of a basement relief consisting of several adjacent discontinuities with small vertical offsets. This is a remarkable advantage over the weighted smoothness method which requires that the interface present few, local discontinuities with large vertical offsets.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-13T12:21:24Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ConstrucaoModeloExperimental.pdf: 4841628 bytes, checksum: bbe77f18c2aa853d7c9afd3b138a6cf5 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:20:00Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-20T13:24:14Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ConstrucaoModeloExperimental.pdf: 4841628 bytes, checksum: bbe77f18c2aa853d7c9afd3b138a6cf5 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T16:56:03Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ConstrucaoModeloExperimental.pdf: 4841628 bytes, checksum: bbe77f18c2aa853d7c9afd3b138a6cf5 (MD5) / Made available in DSpace on 2014-08-22T16:56:03Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ConstrucaoModeloExperimental.pdf: 4841628 bytes, checksum: bbe77f18c2aa853d7c9afd3b138a6cf5 (MD5) Previous issue date: 1990 / A ferramenta de indução tem sido, desde 1949, largamente utilizada na prospecção de petróleo pois auxilia na localização de formações que contém hidrocarbonetos. Ainda existem, porém, várias restrições quanto ao uso deste equipamento. Camadas excessivamente finas, poços preenchidos com lamas muito condutivas e grandes diâmetros de invasão são exemplos de situações onde seu desempenho não é bom. Para que se possa minimizar estas limitações é necessário o conhecimento da resposta da sonda diante de condições conhecidas. Isto pode ser feito via modelamento analógico e/ou através de modelamento numérico. O presente trabalho consiste na montagem de um sistema de modelamento analógico capaz de simular o perfil de duas bobinas coaxiais frente a meios infinitos e camadas espessas invadidas. O meio geológico é representado através de solução salina de diferentes condutividades. O fator de escala linear é igual a 20 e o equipamento executa medidas na faixa dos 50 kHz aos 100 kHz. Este sistema será utilizado em futuros estudos envolvendo disposições de bobinas diferentes das normalmente utilizadas e em estudos de outros problemas encontrados na perfilagem de indução. / The induction well-logging tool has been widely used since 1949 in petroleum exploration because it helps in determining water saturation in potential hydrocarbon zones. There are many restrictions in using this tool. Excessively thin beds, boreholes filled with conductive muds, high resistivity formations, and large invasion diameters are situations where induction log response is poor. To minimize these limitations, it is necessary to know tool response in given conditions. This can be partially acomplished by scale modeling and/or digital modeling. The present thesis describes the construction of an experimental analog model system designed to simulate the two coaxial coils induction tool in three different geological situations. The geological medium is simulated by brine solution of diferent conductivities. The linear scale factor is twenty and the frequency range of the measuring system is 50 kHz to 100 kHz. This system can be utilized in future research work with different coil arrangements from the ones currently being used by the service companies and also to study other problems encountered in induction logging.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-05T12:05:08Z No. of bitstreams: 2 license_rdf: 24131 bytes, checksum: abbf6aca685b0252d838a9d6a83cd265 (MD5) Dissertacao_AnaliseAmbiguidadeMetodos.pdf: 9940360 bytes, checksum: 3150084e158f99d3ff58f006bbb33bbf (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:09:52Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-25T13:29:40Z No. of bitstreams: 2 license_rdf: 24131 bytes, checksum: abbf6aca685b0252d838a9d6a83cd265 (MD5) Dissertacao_AnaliseAmbiguidadeMetodos.pdf: 9940360 bytes, checksum: 3150084e158f99d3ff58f006bbb33bbf (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-28T13:46:38Z (GMT) No. of bitstreams: 2 license_rdf: 24131 bytes, checksum: abbf6aca685b0252d838a9d6a83cd265 (MD5) Dissertacao_AnaliseAmbiguidadeMetodos.pdf: 9940360 bytes, checksum: 3150084e158f99d3ff58f006bbb33bbf (MD5) / Made available in DSpace on 2014-08-28T13:46:38Z (GMT). No. of bitstreams: 2 license_rdf: 24131 bytes, checksum: abbf6aca685b0252d838a9d6a83cd265 (MD5) Dissertacao_AnaliseAmbiguidadeMetodos.pdf: 9940360 bytes, checksum: 3150084e158f99d3ff58f006bbb33bbf (MD5) Previous issue date: 1989 / Este trabalho apresenta uma metodologia para o estudo da ambiguidade na interpretação de dados geofísicos. Várias soluções alternativas, representativas da região de maior ambiguidade no espaço de parâmetros são obtidas, sendo posteriormente grupadas e ordenadas pela análise fatorial modo Q. Esta metodologia foi aplicada a dados sintéticos de campo potencial simulando-se causas de ambiguidade como discretização e truncamento da anomalia e a presença de ruídos aleatório e geológico. Um único prisma foi usado como modelo interpretativo, sendo a espessura a principal causa de ambiguidade tanto na gravimetria como na magnetometria. Segue-se a profundidade do topo sempre associada à espessura, quando o sinal da anomalia é alto. Quando a anomalia tem sinal baixo, a largura torna-se o segundo parâmetro mais importante, também associada à espessura. Ao contrário da presença de interferências geológicas, a presença de ruído aleatório nos campos, não é fator importante na ambiguidade. A aplicação da metodologia a dados reais ilustra o papel desta análise na caracterização de soluções alternativas e a importância da informação a priori na caracterização das causas de ambiguidade. A metodologia apresentada pode ser empregada em diversos estágios de um programa de prospecção fornecendo em cada estágio uma análise dos principais fatores causadores da ambiguidade, que poderá ser util no planejamento dos estágios seguintes. Comparada a outros métodos de análise de ambiguidade, como por exemplo regiões de confiança, a metodologia estudada destaca-se por não precisar satisfazer premissas estatísticas sobre a distribuição dos erros. / A method is presented to analyze ambiguity in geophysical interpretation. Initially, alternative solutions are obtained which are assumed to be representative of the region in parameter space where ambiguity is more pronounced. In a second stage, these solutions are grouped and ordered using Q-mode factor analysis. The method was applied to synthetic potential field data where some important causes of ambiguity are simulated, such as discretizing and truncation of the anomaly, and the presence of random and geologic noise. A single prism is employed as an interpretation model and in both gravity and magnetic cases the prism depth extent stands out as the main cause of ambiguity. The second cause is either the depth to the top when the anomaly presents a strong signal, or the width when the anomaly presents a weak signal. The presence of random noise is not an important factor leading to ambiguity. On the other hand, the presence of interfering sources is a strong fator in both gravity and magnetic interpretation. The application of the method to real anomalies illustrates its importance in defining alternative solutions, and the importance of the a priori information in the establishment of the causes of ambiguity. The presented method may be applied at any stage of an exploration project. At each stage it provides relevant information which may be useful in planning the next stage. Compared with other methods usually employed in ambiguity analysis such as confidence regions, for exemple, the proposed method has the advantage of requiring no statistical assumptions about the error distribuition.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-06T15:47:34Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnomaliasEletromagneticasCorpos.pdf: 9640358 bytes, checksum: 0ba56381c85eb0ebb5ef73165e7776b0 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-06T14:54:14Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-26T13:34:56Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnomaliasEletromagneticasCorpos.pdf: 9640358 bytes, checksum: 0ba56381c85eb0ebb5ef73165e7776b0 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T12:01:28Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnomaliasEletromagneticasCorpos.pdf: 9640358 bytes, checksum: 0ba56381c85eb0ebb5ef73165e7776b0 (MD5) / Made available in DSpace on 2014-09-10T12:01:28Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnomaliasEletromagneticasCorpos.pdf: 9640358 bytes, checksum: 0ba56381c85eb0ebb5ef73165e7776b0 (MD5) Previous issue date: 1997 / UFPA - Universidade Federal do Pará / PETROBRAS - Petróleo Brasileiro S.A. / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / A principal aplicação dos métodos eletromagnéticos é a prospecção de sulfetos maciços que pela própria natureza geológica são corpos longos e delgados, encrustados em rocha encaixante resistiva, permitindo a interpretação através de curvas de modelos reduzidos usando semi-planos em meio resistivo (ar-livre). No entanto, com a extensão do uso dos métodos eletromagnéticos em regiões que tem um manto de intemperismo de parcialmente condutivo a condutivo, como na região Amazônica e nas regiões semi-áridas ou de climas tropicais , esta técnica de interpretação tem levado a resultados bem diferentes da situação real. O objetivo deste trabalho é estudar a influência do manto de intemperismo em contato ôhmico com o corpo-alvo, utilizando os arranjos de bobinas horizontal coplanar (HCP), vertical coplanar (VCP), vertical coaxial (VCA) e mínimo (PERP), através do modelagem analógico. Para simular o corpo foram utilizadas placas de grafite, cujas dimensões satisfazem a condição de semi-plano. Para simular o manto foi usada uma solução de cloreto de amônia. De forma geral, os resultados obtidos com os diversos arranjos apresentaram as mesmas características, ressalvada as diferenças peculiares de cada arranjo. De forma resumida tem-se as seguintes alterações na anomalia atribuídas à um manto condutivo em contato ôhmico com o corpo-alvo: i) Há rotação de fase, que inicialmente é em sentido anti-horário, mas mais adiante, com o aumento do número de indução do manto , torna-se horária; ii) Alteração na forma padrão do perfil de quadratura devido ao surgimento do pico extra no lado a favor do mergulho; iii) Só há uniformidade de comportamento para o pico-a-pico contra da quadratura, que cresce com o aumento número de indução do manto. Maiores anomalias são obtidas com o sistema horizontal coplanar e as menores com o vertical coplanar. Neste estudo não foi observado formação de pico extra nas anomalias do arranjo vertical coplanar. Ainda, os efeitos de cobertura são mais acentuados no sistema horizontal coplanar comparados aos sistemas de bobina vertical e em anomalias de semi-planos menos inclinados e em menor profundidade. / Electromagnetic methods are mainly applied in the prospection of massive sulphide deposits. These deposits are frequently long in extension and thin, therefore the interpretation can be carried out considering a half-plane model in a resistive ambient. However, in Amazon region or in semi-arid and tropical regions where a partialy conductive to conductive overburden is encountered, interpretations using a resistive ambient model give erroneous results. The objective of this work is to study the influence of the conducting overburden on electromagnetic anomalies of inclined sheet-like bodies in ohmic contact with the overburden, using coil configurations - horizontal coplanar (HCP), vertical coplanar (VCP), vertical coaxial (VCA) and minium couple system (PERP). These studies are realized with the aid of analog modelling where the ore-body is simulated by a graphite sheet and the conducting overburden by ammonium choride solution. Generally, the diferente coil systems studied here show the similar characteristics of the overburden effects, reserving some peculiar effects of each coil system. The effects of a conductive overburden in ohmic contact with the orebody could be briefly described as follows: (i) There is an anti-clockwise phase rotation initially with the increase in the overburden induction number, but with further increase the phase rotates clockwise; (ii) The form of the quadrature anomaly profile is altered due to the appearence of an extra peak on the down-dip side of the ore-body, and (iii) The peak to peak quadrature anomaly on the up-dip side increases uniformly with increase in the overburden induction number. Highest anomalies are obtained with the horizontal coplanar system and lowest with vertical coplanar. No extra-peak in the quadrature profile is observed with the vertical coaxial system in this study. Also, the overbuden effects are more prominent in the horizontal coplanar system compared to vertical coil systems, and in the anomalies of the less inclined half-planes and at lesser depths.

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-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-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.