Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-23T13:38:42Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoGravimetricaRapida.pdf: 2097359 bytes, checksum: 644d6d953cc244368b99024e23a9dbad (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-29T11:50:51Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoGravimetricaRapida.pdf: 2097359 bytes, checksum: 644d6d953cc244368b99024e23a9dbad (MD5) / Made available in DSpace on 2014-09-29T11:50:51Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoGravimetricaRapida.pdf: 2097359 bytes, checksum: 644d6d953cc244368b99024e23a9dbad (MD5) Previous issue date: 2013 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Apresentamos um novo método de inversão linear bidimensional de dados gravimétricos produzidos por bacias sedimentares com relevo do embasamento descontínuo. O método desenvolvido utiliza um modelo interpretativo formado por um conjunto de fitas horizontais bidimensionais justapostas cujas espessuras são os parâmetros a serem estimados. O contraste de densidade entre o embasamento e os sedimentos é presumido constante e conhecido. As estimativas das espessuras foram estabilizadas com o funcional da Variação Total (VT) que permite soluções apresentando descontinuidades locais no relevo do embasamento. As estimativas do relevo são obtidas através da resolução de um sistema de equações lineares, resolvido na norma L1. Como métodos lineares subestimam as estimativas de profundidade do embasamento de bacias maiores que cerca de 500 m, amplificamos as estimativas de profundidade através da modificação da matriz associada ao modelo interpretativo de fitas. As estimativas obtidas através deste procedimento são em geral ligeiramente superestimadas. Desse modo, elas são corrigidas através de uma correção definida pela expressão da placa Bouguer. Testes em dados sintéticos e reais produziram resultados comparáveis aos produzidos pelo método não linear, mas exigiram menor tempo computacional. A razão R entre os tempos exigidos pelo método não linear e o método proposto cresce com o número de observações e parâmetros. Por exemplo, para 60 observações e 60 parâmetros, R é igual a 4, enquanto para 2500 observações e 2500 parâmetros R cresce para 16,8. O método proposto e o método de inversão não linear foram aplicados também em dados reais do Steptoe Valley, Nevada, Estados Unidos, e da ponte do POEMA, no Campus do Guamá em Belém, produzindo soluções similares às obtidas com o método não linear exigindo menor tempo computacional. / We present a new 2D linear inversion gravity method to estimate the discontinuous basement relief of a sedimentary basin. The method uses an interpretation model consisting of a set of 2D juxtaposed horizontal ribbons whose thicknesses are the parameters to be estimated. The density contrast between the basement and the sediments is assumed constant and known. The estimates of the thicknesses were stabilized with the Total Variation functional (TV) which permits solutions presenting local discontinuities in the basement relief. The estimates of the basement relief are produced by solving a single linear system of equations in the L1 norm. Because linear inversion methods underestimate the depth of the basement relief for sedimentary basins thicker than 500 m, we amplify the depth estimates through by modifying the matrix associated with the interpretation model. This procedure in general leads to overestimated depths which are corrected through the expression of the anomaly of a Bouguer slab. Tests on synthetic and real data produced similar results as compared with the nonlinear method, but required less computational time. The ratio R between the computational times required by the nonlinear and the proposed methods to produce similar solutions increases with the number of observations and parameters. For example, for 60 observations and 60 parameters R is equal to 4 whereas for 2500 observations and 2500 R increases to 16.8. The proposed method and the nonlinear inversion method were also applied to real gravity data from Steptoe Valley, Nevada, U.S.A. and from Ponte do POEMA at the Campus of Guamá, in Belém, producing solutions which are similar to the ones produced by the nonlinear method but required a smaller computer time.