Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-06T11:33:01Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseResolucaoTomografia.pdf: 4620366 bytes, checksum: 020dce7bff43549e65896f9f48082777 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-06-26T14:58:18Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseResolucaoTomografia.pdf: 4620366 bytes, checksum: 020dce7bff43549e65896f9f48082777 (MD5) / Made available in DSpace on 2014-06-26T14:58:18Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseResolucaoTomografia.pdf: 4620366 bytes, checksum: 020dce7bff43549e65896f9f48082777 (MD5) Previous issue date: 2000 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / A aplicação de métodos de inversão aos dados de múltiplos levantamentos sísmicos interpoços permite a reconstrução de modelos de vagarosidade em 3-D de alta resolução adequados para monitoramento de processos de recuperação avançada de petróleo e caracterização de reservatórios. Entretanto, a falta de cobertura volumétrica uniforme de raios de levantamentos interpoços exige informação adicional ao sistema tomográfico para obtenção de soluções estáveis. A discretização do modelo em uma malha 3-D com células prismáticas triangulares e a decomposição em valores singulares são utilizadas para avaliar a reconstrução tomográfica em 3-D. O ângulo da projeção de modelos-alvo no subespaço ortogonal ao espaço nulo efetivo da matriz tomográfica é um critério adequado para se otimizar a malha de discretização do modelo interpretativo e a geometria de aquisição dos dados de modo a melhorar o condicionamento da reconstrução tomográfica. Esta abordagem pode ser utilizada durante as iterações lineares para redefinir a malha ou avaliar a necessidade de informação a priori adicional ao sistema tomográfico. / The joint inversion of multiple crosswell surveys traveltime data allows the reconstruction of high resolution 3-D slowness models suitable for monitoring enhanced oil recovery processes and detailed reservoir characterization. However, the non-uniform volumetric ray coverage unavoidable in crosswell surveys requires additional information to the tomographic system in order to obtain stable solutions. The 3-D tomographic problem is analyzed discretizing the models in 3-D prismatic meshes and using the singular value decomposition. The projection angle of target models in the subspace orthogonal to the effective null space of the tomographic matrix is an useful criteria for designing mesh and data acquisition geometry in order to better conditioning the tomographic reconstruction. The approach can also be used during the linear iterations to redefine the mesh or evaluate the requirement of additional prior information to the tomographic system.

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-06T12:37:09Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseVelocidadeOtimizacao.pdf: 2382926 bytes, checksum: 413bbca447f7bac1271afc2aae31c08d (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-06T12:35:46Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-07T12:54:14Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseVelocidadeOtimizacao.pdf: 2382926 bytes, checksum: 413bbca447f7bac1271afc2aae31c08d (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-18T14:09:19Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseVelocidadeOtimizacao.pdf: 2382926 bytes, checksum: 413bbca447f7bac1271afc2aae31c08d (MD5) / Made available in DSpace on 2014-08-18T14:43:01Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseVelocidadeOtimizacao.pdf: 2382926 bytes, checksum: 413bbca447f7bac1271afc2aae31c08d (MD5) Previous issue date: 2011 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / FINEP - Financiadora de Estudos e Projetos / PETROBRAS - Petróleo Brasileiro S.A. / Este trabalho teve como objetivo geral desenvolver uma metodologia sistemática para a inversão de dados de reflexão sísmica em arranjo ponto-médio-comum (PMC), partindo do caso 1D de variação vertical de velocidade e espessura que permite a obtenção de modelos de velocidades intervalares, vint,n, as espessuras intervalares, zn, e as velocidades média-quadrática, vRMS,n, em seções PMC individualizadas. Uma consequência disso é a transformação direta destes valores do tempo para profundidade. Como contribuição a análise de velocidade, foram desenvolvidos dois métodos para atacar o problema baseado na estimativa de velocidade intervalar. O primeiro método foi baseado na marcação manual em seções PMC, e inversão por ajuste de curvas no sentido dos quadrados-mínimos. O segundo método foi baseado na otimização da função semblance para se obter uma marcação automática. A metodologia combinou dois tipos de otimização: um Método Global (Método Price ou Simplex), e um Método Local (Gradiente de Segunda Ordem ou Conjugado), submetidos a informação à priori e vínculos. A marcação de eventos na seção tempo-distância faz parte dos processos de inversão, e os pontos marcados constituem os dados de entrada juntamente com as informações à priori do modelo a ser ajustado. A marcação deve, por princípio, evitar eventos que representem múltiplas, difrações e interseções, e numa seção pode ser feita mais de 50 marcações de eventos, enquanto que num mapa semblance não se consegue marcar mais de 10 eventos de reflexão. A aplicação deste trabalho é voltada a dados sísmicos de bacias sedimentares em ambientes marinhos para se obter uma distribuição de velocidades para a subsuperfície, onde o modelo plano-horizontal é aplicado em seções PMC individualizadas, e cuja solução pode ser usada como um modelo inicial em processos posteriores. Os dados reais da Bacia Marinha usados neste trabalho foram levantados pela PETROBRAS em 1985, e a linha sísmica selecionada foi a de número L5519 da Bacia do Camamu, e o PMC apresentado é a de número 237. A linha é composta de 1098 pontos de tiro, com arranjo unilateraldireito. O intervalo de amostragem é 4 ms. O espaçamento entre os geofones é 13,34 m com o primeiro geofone localizado a 300 m da fonte. O espaçamento entre as fontes é de 26,68 m. Como conclusão geral, o método de estimativa de velocidade intervalar apresentada neste trabalho fica como suporte alternativo ao processo de análise de velocidades, onde se faz necessário um controle sobre a sequência de inversão dos PMCs ao longo da linha sísmica para que a solução possa ser usada como modelo inicial ao imageamento, e posterior inversão tomográfica. Como etapas futuras, podemos propor trabalhos voltados direto e especificamente a análise de velocidade sísmica estendendo o caso 2D de otimização do semblance ao caso 3D, estender o presente estudo para o caso baseado na teoria do raio imagem com a finalidade de produzir um mapa continuo de velocidades para toda a seção sísmica de forma automática. / This work had as a general aim to develop a systematic methodology for the inversion of seismic reflection data organized in common-midpoint gathers (CMP), starting from 1D vertical variation of velocity and thickness that allows to obtain interval velocity, vint,n, in time, the correspondent interval thickness, zn, and the correspondent mean-square velocity, vRMS,n in individualized CMP gathers. A direct consequence of this work the transformation of these values from time to depth. Two methods were developed to attack the problem defined as velocity analysis based on the estimation of interval velocity. The first method was based on manual picking of reflection events on CMP gathers, and inversion by curve fitting in the least-square sense. The second method was based on the otimization of the semblance function to obtain an automatic picking. The methodology combined two types of optimization: a Global Method (Price or Simplex), and Local Method (second order gradient or cojugate), subject to a priori information and constraints. The picking of events in time-distance section is of fundamental importance in the process of inversion, and the picked points are the input data along with a priori information of the model to be adjusted. The picking must, in principle, avoid events that represent multiples, diffractions and intersections, and in a section over 50 pickings can be made, while in a semblance map not more than 10 events could usually be picked by eye. The application of this work is focused on seismic data of marine sedimentary basins to obtain a distribution of velocities for the subsurface, where a plane-horizontal model is applied for individual CMP sections, and that the solution can be used as an initial model in subsequent processes. The real data used in this study were collected by Petrobras in 1985, and the selected seismic line was of number L5519 of the Camamu Basin, and the CMP presented is of number 237. The line consists of 1098 shot points with right-lateral arrangement. The sampling interval is 4 ms. The spacing between the geophones is 13.34 m with the first geophone located at 300 m from the source. The spacing between the sources is 26.68 m. As a general conclusion, the method for estimating interval velocity in this work stands as an alternative support to velocity analysis, where it is necessary a control over the sequential inversion of CMP gathers along the seismic line such that the solution can be used as an initial model for imaging, and further tomographic inversion. As future work, we can be propose studies directely and specifically related to seismic velocity analysis by extending the 2D semblance optimization method to 3D, and extending the present studies to the method based on the image ray, aiming at producing a continuous velocity map for the entire section in an automatic way.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-28T12:03:20Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InversaoMomentosFonte.pdf: 26856422 bytes, checksum: 4c2d4fd4c3ccfc0072a0e1986811b2a7 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:14:12Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-19T12:49:34Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InversaoMomentosFonte.pdf: 26856422 bytes, checksum: 4c2d4fd4c3ccfc0072a0e1986811b2a7 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-20T11:56:40Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InversaoMomentosFonte.pdf: 26856422 bytes, checksum: 4c2d4fd4c3ccfc0072a0e1986811b2a7 (MD5) / Made available in DSpace on 2014-08-20T11:56:41Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_InversaoMomentosFonte.pdf: 26856422 bytes, checksum: 4c2d4fd4c3ccfc0072a0e1986811b2a7 (MD5) Previous issue date: 1993 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / UFPA - Universidade Federal do Pará / A inversão de momentos de fonte gravimétrica tridimensional é analisada em duas situações. Na primeira se admite conhecer apenas a anomalia. Na segunda se admite conhecer, além da anomalia, informação a priori sobre o corpo anômalo. Sem usar informação a priori, mostramos que é possível determinar univocamente todo momento, ou combinação linear de momentos, cujo núcleo polinomial seja função apenas das coordenadas Cartesianas que definem o plano de medida e que tenha Laplaciano nulo. Além disso, mostramos que nenhum momento cujo núcleo polinomial tenha Laplaciano não nulo pode ser determinado. Por outro lado, informação a priori é implicitamente introduzida se o método de inversão de momentos se baseia na aproximação da anomalia pela série truncada obtida de sua expansão em multipolos. Dado um centro de expansão qualquer, o truncamento da série impõe uma condição de regularização sobre as superfícies equipotenciais do corpo anômalo, que permite estimar univocamente os momentos e combinações lineares de momentos que são os coeficientes das funções-bases da expansão em multipolos. Assim, uma distribuição de massa equivalente à real é postulada, sendo o critério de equivalência especificado pela condição de ajuste entre os campos observado e calculado com a série truncada em momentos de uma ordem máxima pré-estabelecida. Os momentos da distribuição equivalente de massa foram identificados como a solução estacionária de um sistema de equações diferenciais lineares de 1a. ordem, para a qual se asseguram unicidade e estabilidade assintótica. Para a série retendo momentos até 2a. ordem, é implicitamente admitido que o corpo anômalo seja convexo e tenha volume finito, que ele esteja suficientemente distante do plano de medida e que a sua distribuição espacial de massa apresente três planos ortogonais de simetria. O método de inversão de momentos baseado na série truncada (IMT) é adaptado para o caso magnético. Para este caso, mostramos que, para assegurar unicidade e estabilidade assintótica, é suficiente pressupor, além da condição de regularização, a condição de que a magnetização total tenha direção e sentido constantes, embora desconhecidos. O método IMT baseado na série de 2a. ordem (IMT2) é aplicado a anomalias gravimétricas e magnéticas tridimensionais sintéticas. Mostramos que se a fonte satisfaz as condições exigidas, boas estimativas da sua massa ou vetor momento de dipolo anômalo total, da posição de seu centro de massa ou de momento de dipolo e das direções de seus três eixos principais são obtidas de maneira estável. O método IMT2 pode falhar parcialmente quando a fonte está próxima do plano de medida ou quando a anomalia tem efeitos localizados e fortes de um corpo pequeno e raso e se tenta estimar os parâmetros de um corpo grande e profundo. Definimos por falha parcial a situação em que algumas das estimativas obtidas podem não ser boas aproximações dos valores verdadeiros. Nas duas situações acima descritas, a profundidade do centro da fonte (maior) e as direções de seus eixos principais podem ser erroneamente estimadas, embora que a massa ou vetor momento de dipolo anômalo total e a projeção do centro desta fonte no plano de medida ainda sejam bem estimados. Se a direção de magnetização total não for constante, o método IMT2 pode fornecer estimativas erradas das direções dos eixos principais (mesmo se a fonte estiver distante do plano de medida), embora que os demais parâmetros sejam bem estimados. O método IMT2 pode falhar completamente se a fonte não tiver volume finito. Definimos por falha completa a situação em que qualquer estimativa obtida pode não ser boa aproximação do valor verdadeiro. O método IMT2 é aplicado a dados reais gravimétricos e magnéticos. No caso gravimétrico, utilizamos uma anomalia situada no estado da Bahia, que se supõe ser causada por um batólito de granito. Com base nos resultados, sugerimos que as massas graníticas geradoras desta anomalia tenham sido estiradas na direção NNW e adelgaçadas na direção vertical durante o evento compressivo que causou a orogênese do Sistema de Dobramentos do Espinhaço. Além disso, estimamos que a profundidade do centro de massa da fonte geradora é cerca de 20 km. No caso magnético, utilizamos a anomalia de um monte submarino situado no Golfo da Guiné. Com base nos resultados, estimamos que o paleopolo magnético do monte submarino tem latitude 50°48'S e longitude 74°54'E e sugerimos que não exista contraste de magnetização expressivo abaixo da base do monte submarino. / The inversion of three-dimensional gravity source moments is analyzed in two situations. In the first one only the anomalous field is assumed to be known. In the second situation a priori information about the anomalous body is assumed to be known besides the field data. Without using a priori information, we show that it is possible to determine uniquely any moment, or linear combination of moments, whose polynomial kernel: (a) is not a function of the Cartesian coordinate which is orthogonal to the measuring plane and (b) has null Laplacian. Besides, we show that it is impossible to determine any moment whose polynomial kernel has non-null Laplacian. On the other hand, we show that a priori information is implicitly introduced if the source moment inversion method is based on the approximation of the anomalous field by the truncated series obtained from its multipole expansion. Given any center of expansion, the series truncation impores a regularization condition on the equipotential surfaces of the anomalous body that allows the moments and linear combination of moments (which are the coefficients of the multipole expansion basis function) to be uniquely estimated. So, a mass distribution equivalent to the real mass distribution is postulated, being the equivalence criterion specified by the fitting conditions between the observed anomaly and the anomaly calculated with the truncated multipole expansion series. The highest order for the retained terms in the truncated series is specified by the previously defined maximum order for the moments. The moments of the equivalent mass distribution were identified as the stationary solution of a system of first order linear differential equations, for which uniqueness and assymptotic stability are assured. For the series having moments up to 2nd order, it is implicitly assumed that the anomalous body: (1) has finite volume, (2) that it is sufficiently far from the measuring plane and (3) that its spatial naass distribution is convex and presents three orthogonal planes of symmetry. The source moment inversion method based on the approximation of the anomalous field by a truncated series (MIT) is adapted to the magnetic case. In this case, we show that in order to guarantee uniqueness and assymptotic stability it is sufficient to assume, besides the regularization condition, that the total magnetization has constant but unknown direction. The MIT method based on the 2nd order series (MIT2) is applied to three-dimensional synthetic gravity and magnetic anomalies. If the source satisfies all imposed conditions, we show that it is possible to obtain in a stable way good estimates of the total anomalous mass or dipole moment vector, of the position of center of mass or dipole moment and of the directions of all three principal axes. A partia' failure of MIT2 method may occur either if the source is dose to the measuring plane or if the anomaly presents a localized but strong effect due to a shallow and small body and an attempt is made to estimate the moments of a large and deep body. By partial failure we mean the situation when some of the estimates may be poor aproximations of the true values. In these two cases we show that the estimates of the depth and the directions of the principal axes of the (main) source may be poor but the estimates of the total anomalous mass or dipole moment vector and the projection on the measuring plane of the center of mass or dipole moment of the source are good. If the total magnetization direction is not constant, MIT2 method may produce poor estimates of the directions of the principal axes (even if the source is far from the measuring plane) but good estimates are obtained for the other parameters. A complete failure of MIT2 method may occur if the source does not have finite volume. By complete failure we mean the situation when any obtained estimate may be a poor aproximation of the true value. MIT2 method is applied to real gravity and magnetic data. In the gravimetric case we used an anomaly located in Bahia state, Brazil, which is assumed to be produced by the presence of a large granitic body. Based on the inversion results, we propose that the grafite was deformed into an oblate ellipsoid during the compressive event that generated the Middle Proterozoic Espinhaço orogeny. The center of mass estimated for this body is about 20 km. In the magnetic case, we used an anomaly produced by a seamount located in the Gulf of Guinea. Based on the inversion results, we estimate a magnetic palaeopole for the seamount at 50°48'S and 74°54'E and we suggest that no important magnetization contrast exists below the bottom of the seamount.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-06T16:39:28Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoAnaliseSeries.pdf: 9563203 bytes, checksum: 19126d425ed2de3d0aeb62c84c46a64e (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar as palavras-chave on 2014-08-07T13:42:22Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-20T12:25:05Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoAnaliseSeries.pdf: 9563203 bytes, checksum: 19126d425ed2de3d0aeb62c84c46a64e (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-20T13:20:01Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoAnaliseSeries.pdf: 9563203 bytes, checksum: 19126d425ed2de3d0aeb62c84c46a64e (MD5) / Made available in DSpace on 2014-08-20T13:20:01Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoAnaliseSeries.pdf: 9563203 bytes, checksum: 19126d425ed2de3d0aeb62c84c46a64e (MD5) Previous issue date: 1993 / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / FINEP - Financiadora de Estudos e Projetos / Esta tese descreve a aplicação de análise de séries temporais em perfis de poço. Através desta técnica é possível avaliar-se a repetição e a resolução vertical de perfis, e determinar-se o intervalo de amostragem e a velocidade de perfilagem ideais para diferentes perfis. A comparação entre três poços é também feita, baseada num mesmo tipo de perfil. Para tanto, na seqüência utilizada, procurou-se manter num mesmo domínio os dados cuja quantidade total de amostras (N) por perfil não ultrapassou 2048. Desses dados, foram inicialmente retirados o valor médio das amostras e o alinhamento polinomial algébrico eventualmente nelas embutido. Em seguida, foi efetuada a aplicação do ponderador cossenoidal, do filtro passa-alta, da janela Hanning, do cálculo da função coerência, do espectro de fase, da razão sinal-ruído e dos espectros de potência do sinal e do ruído, nesta ordem. Para a função coerência, fez-se necessário o cálculo dos níveis de confiança de 50%, 90% e 95%. O cálculo do primeiro nível teve por base a necessidade de se determinar a resolução vertical de alguns perfis, e dos demais, a fim de que fosse obtida uma informação referente à localização daqueles níveis para a coerência calculada. Em relação ao espectro de fase, seu cálculo surgiu da necessidade de se obter uma informação adicional a respeito dos perfis manipulados, ou seja, o conhecimento da ocorrência ou não de deslocamento relativo de profundidade entre a seção principal e a seção repetida. A razão sinal-ruído foi calculada no sentido de possibilitar a comparação, como elemento avaliador dos diversos tipos de perfis, com a coerência e o cálculo dos espectros de potência. Os espectros de potência do sinal e do ruído foram calculados para se ter mais um parâmetro de avaliação da seção repetida, já que em tese, os espectros de potência do sinal e do ruído da seção repetida devem ser iguais aos respectivos espectros da seção principal. Os dados utilizados na aplicação da metodologia proposta foram fornecidos pela PETROBRÁS e oriundos de quatro poços da Bacia Potiguar emersa. Por questões de sigilo empresarial, os poços foram identificados como poços A, B, C e D. A avaliação da repetição entre diferentes tipos de perfis indica que, para o poço A, o perfil micro-esférico (MSFL) tem repetição melhor do que o perfil de porosidade neutrônica (CNL), o qual tem, por sua vez, repetição melhor do que o perfil de raios gama normal (GR). Para os perfis do poço D, uma diminuição da velocidade de perfilagem de 550 m/h para 275 m/h é vantajosa apenas para o perfil de porosidade neutrônica. Já a velocidade de perfilagem de 920 m/h, utilizada, na obtenção dos perfis do poço C, é totalmente inadequada para os perfis de resistividade (MSFL, ILD e ILM). A diminuição do intervalo de amostragem de 0,20 m para 0,0508 m, nos perfis de raios gama e de porosidade neutrônica, e 0,0254 m para o perfil de densidade, apresenta bons resultados quando aplicada no poço D. O cálculo da resolução vertical indica, para o perfil de porosidade neutrônica, uma superioridade qualitativa em relação ao perfil de raios gama normal, ambos pertencentes ao poço A. Para o poço C, o perfil micro-esférico apresenta uma resolução vertical na mesma ordem de grandeza da resolução do perfil de raios gama do poço B, o que evidencia ainda mais a inconveniência da velocidade de perfilagem utilizada no poço C. Já para o poço D, o cálculo da resolução vertical indica uma superioridade qualitativa do perfil de densidade de alta resolução em relação ao perfil de raios gama de alta resolução. A comparação entre os poços A, B e D, levada a efeito através dos respectivos perfis de porosidade neutrônica normais, comprova que a presença de ruído aleatório, em geral, está diretamente ligada à porosidade da formação - uma maior porosidade indica uma presença maior de ruído e, por conseguinte, uma queda qualitativa no perfil obtido. A análise do espectro de fase de cada perfil indica um deslocamento em profundidade, existente entre as seções principal e repetida de todos os perfis do poço C. E isto pôde ser confirmado com a posterior superposição das seções. / This thesis describes the application of time seris analysis to wireline logs. Through this technique it is possible to evaluate both their repeatability and vertical resolution, and determine the optimum sampling interval and acquisition speed for different logs. A comparison between three wells is also made, based on the same type of log. The sequence used, is to obtain data, in the same domain, for which the number of samples (N) does not exceed 2048. For these data the mean sample value and the algebraic polynomial are determined. The following were then applied, in order: cosine taper, high pass filter, Hanning window, the calculation of the coherence function, the phase spectra, and the signal to noise ratio of the power spectra for both the signal and the noise. For the coherence function, it was necessary to calculate the level of confidence for 50 necessary to determine the vertical resolution of some logs. The others were calculated to provide information concerning the position of the levels of coherence calculated. In relation to the phase spectra, its calculation suggests it is necessary to obtain additional information in respect of the processed logs, or in other words knowledge of any relative depth shifts made between the principal and repeat sections. The signal to noise ratio was calculated to investigate the possibility of evaluating the different types of logs by making a comparison with the coherence and the calculated power spectra. The power spectra of the signal and the noise were calculated to provide one additional parameter to evaluate the repeat section. In theory the power spectra of the signal and the noise of the repeat section should be equal to their respective spectra for the principal sections. The data used in this work were provided by PETROBRAS and originated in four boles of the Potiguar Basin. These are referred to as Holes A, B, C, and D. The evaluation of the repeatability between different types of logs indicates that, for Hole A, the microspherical log (MSFL) has better repeatability than the porosity log (CNL), and that this repeats better than the gamma ray log (GR). For the logs run in Hole D, a decrease in the speed at which the logs are run, from 550 m/hr to 275 m/hr, is advantageous only for the neutron porosity log. The velocity of 920 m/hr used in Hole C is totally inadequate for the resistivity logs (MSFL, ILD, and ILM). A decrease in the sampling interval from 0.2m to 0.0508m, for the gamma ray and neutron porosity logs, and to 0.0254m for the density log, gives good results when applied to Hole D. The calculation of vertical resolution indicates, that the neutron porosity is of superior quality when compared with the natural gamma ray in Hole A. In hole C the microspherical log (MSFL) has a vertical resolution of the same order as the natural gamma ray log (GR) in Hole B. This confirms the inappropriate speed for the logs acquired in Hole C. For Hole D, the calculation of vertical resolution indicates a superior quality for the high resolution density log when compared with the high resolution gamma ray log. The comparison between Holes A, B and D, uses the respective neutron porosity logs, confirming that the presence of noise, in general, is directly linked with the porosity of the formation - a higher porosity indicates more noise and consequently, a lower quality log is obtained. Analysis of phase spectra of each log indicates a depth shift exists between the principal and repeat sections for all the logs of Hole C. This could be confirmed through superposition of the sections.

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-18T12:26:26Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoFraturasPoco.pdf: 6764040 bytes, checksum: fb2713019e47094663a20f610f82b95d (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Incluir palavras-chave on 2014-08-07T16:02:57Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-20T16:48:10Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoFraturasPoco.pdf: 6764040 bytes, checksum: fb2713019e47094663a20f610f82b95d (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T16:57:12Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoFraturasPoco.pdf: 6764040 bytes, checksum: fb2713019e47094663a20f610f82b95d (MD5) / Made available in DSpace on 2014-08-22T16:57:12Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoFraturasPoco.pdf: 6764040 bytes, checksum: fb2713019e47094663a20f610f82b95d (MD5) Previous issue date: 1992 / A aplicabilidade de uma ferramenta de indução coplanar de duas bobinas na detecção de zonas fraturadas é estudado através de modelamento reduzido em laboratório. O fator de escala escolhido foi 10 e as condutividades de campo foram representadas por soluções de Cloreto de Amônia. O modelo geológico é composto por uma rocha homogênea resistiva com porções fraturadas e não-fraturada. As resistividades de campo são 4 Ω-m e 200 Ω-m para a zona fraturada e zona não-fraturada, respectivamente. As zonas fraturadas estudadas eram verticais, a várias distâncias da parede do poço, e cruzando a perfuração com mergulhos de 0°, 45° e 60°. A maioria dos perfis foram corridos em uma freqüência de laboratório de 100 kHz. Os resultados apresentados sob a forma de perfis mostraram que a ferramenta coplanar fornece informações sobre o mergulho, espessura e distância da parede do poço, dependendo da posição relativa das bobinas e da zona fraturada. Como a ferramenta coplanar não fornece informação de toda a formação ao redor da mesma, foi sugerido a construção de um sistema de duplo acoplamento semi-nulo (DASN) e alguns perfis são apresentados. Além do maior recobrimento, essa ferramenta fornece uma "técnica de olhada rápida" (quick-look technique) para estimar o mergulho da zona fraturada qualitativamente. / The applicability of a coplanar induction two-coil system in the detection of a fractured zone is studied through reduced scale model experiments in the laboratory. The scale factor chosen was 10 and the field conductivities were represented by Amonium Chloride solutions. The geologic model was compounded by an homogeneous resistive rock with fractured and non-fractured parts. The field resistivities are 4 Ω-m and 200 Ω-m of the fractured zone and of the non-fractured zone, respectively. The fractured zones studied are vertical, at varying distances from the borehole wall and those crossing the well are dipping at 0°, 45° and 60°. The majority of the logs were run in a laboratory frequency of 100 kHz. The results presented in the form of profiles (S/P%) show that the coplanar tool gives information about the dip, thickness and distance of the borehole wall depending on the relative position of the fractured zone and axes of the coils. Some blind zones were found to exist depending upon the geometry of fractures in relation with the coplanar system. To resolve this problem, we suggested a Semi—Null Double Coupling (SNDC) System and some logs are presented. The SNDC System provides a quick-look technique to estimate qualitatively the dip of the fractured zone.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-10T13:43:12Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ArranjosCoplanarCoaxial.pdf: 1708656 bytes, checksum: d6198612e2f0ea877bfc56118d676869 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:16:28Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-20T14:17:29Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ArranjosCoplanarCoaxial.pdf: 1708656 bytes, checksum: d6198612e2f0ea877bfc56118d676869 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T17:01:06Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ArranjosCoplanarCoaxial.pdf: 1708656 bytes, checksum: d6198612e2f0ea877bfc56118d676869 (MD5) / Made available in DSpace on 2014-08-22T17:01:06Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ArranjosCoplanarCoaxial.pdf: 1708656 bytes, checksum: d6198612e2f0ea877bfc56118d676869 (MD5) Previous issue date: 2007 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / A ferramenta de indução triaxial ou multicomponente foi desenvolvida recentemente para solucionar os problemas de medição da convencional perfilagem de indução em formações anisotrópicas. Neste trabalho são apresentadas as respostas da sonda de indução triaxial em reservatórios que possuem seqüências finamente laminadas de arenito e folhelho. Reservatórios como estes apresentam a conhecida anisotropia macroscópica, portanto, os resultados demonstrarão a existência de uma resistividade horizontal e outra vertical. Devido a sonda triaxial de indução apresentar dois pares de bobinas coplanares, os quais em geral fornecem os valores de resistividade vertical, este trabalho apresenta inicialmente um estudo analítico e qualitativo das respostas obtidas pelo arranjo coplanar onde o transmissor é tratado como um dipolo. Também será apresentada a solução analítica em que a bobina transmissora é considerada como um loop de dimensão finita. / The triaxial or multicomponent induction tool was recently developed to solve the problems of conventional induction logging measurements in anisotropic formations. In this work are presented the triaxial induction sonde answers in reservoirs with thinly laminated sand-shale sequences. Reservoirs as these present the known macrocospic anisotropy, therefore, the results will demonstrate to the existence of both horizontal and vertical resistivity. Due the triaxial induction sonde to present two pairs of coplanar coils, which in general provide the values of vertical resistivity, this work presents initially an analytical and qualitative study of the answers obtained for coplanar array where the transmitter is considered a dipole. Also will be presented the analytical solution where the transmitter coil is evaluated as a loop of finite dimension.

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

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