Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-23T11:46:42Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoDadosSismicos.pdf: 9047067 bytes, checksum: 81c978f59e97e716df5d90a593c7766b (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-29T13:47:29Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoDadosSismicos.pdf: 9047067 bytes, checksum: 81c978f59e97e716df5d90a593c7766b (MD5) / Made available in DSpace on 2014-09-29T13:47:29Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoDadosSismicos.pdf: 9047067 bytes, checksum: 81c978f59e97e716df5d90a593c7766b (MD5) Previous issue date: 2007 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / A presente Dissertação de Mestrado tem como objetivo o estudo do problema de inversão sísmica baseada em refletores planos para arranjo fonte-comum (FC) e ponto-médiocomum (PMC). O modelo direto é descrito por camadas homogêneas, isotrópicas com interfaces plano-horizontais. O problema é relacionado ao empilhamento NMO baseado na otimização da função semblance, para seções PMC corrigidas de sobretempo normal (NMO). O estudo foi baseado em dois princípios. O primeiro princípio adotado foi de combinar dois grupos de métodos de inversão: um Método Global e um Método Local. O segundo princípio adotado foi o de cascata, segundo a teoria Wichert-Herglotz-Bateman, que estabelece que para conhecer uma camada inferior tem-se que conhecer primeiro a camada superior (dissecação). A aplicação do estudo é voltada à simulação sísmica de Bacia Sedimentar do Solimões e de Bacia Marinha para se obter uma distribuição local 1D de velocidades e espessuras para a subsuperfície em horizontes alvo. Sendo assim, limitamos a inversão entre 4 e 11 refletores, uma vez que na prática a indústria limita uma interpretação realizada apenas em número equivalente de 3 a 4 refletores principais. Ressalta-se que este modelo é aplicável como condição inicial ao imageamento de seções sísmicas em regiões geologicamente complexas com variação horizontal suave de velocidades. Os dados sintéticos foram gerados a partir dos modelos relacionados a informações geológicas, o que corresponde a uma forte informação a priori no modelo de inversão. Para a construção dos modelos relacionados aos projetos da Rede Risco Exploratório (FINEP) e de formação de recursos humanos da ANP em andamento, analisamos os seguintes assuntos relevantes: (1) Geologia de bacias sedimentares terrestre dos Solimões e ma rinha (estratigráfica, estrutural, tectônica e petrolífera); (2) Física da resolução vertical e horizontal; e (3) Discretização temporal-espacial no cubo de multi-cobertura. O processo de inversão é dependente do efeito da discretização tempo-espacial do campo de ondas, dos parâmetros físicos do levantamento sísmico, e da posterior reamostragem no cubo de cobertura múltipla. O modelo direto empregado corresponde ao caso do operador do empilhamento NMO (1D), considerando uma topografia de observação plana. O critério básico tomado como referência para a inversão e o ajuste de curvas é a norma 2 (quadrática). A inversão usando o presente modelo simples é computacionalmente atrativa por ser rápida, e conveniente por permitir que vários outros recursos possam ser incluídos com interpretação física lógica; por exemplo, a Zona de Fresnel Projetada (ZFP), cálculo direto da divergência esférica, inversão Dix, inversão linear por reparametrização, informações a priori, regularização. A ZFP mostra ser um conceito út il para estabelecer a abertura da janela espacial da inversão na seção tempo-distância, e representa a influência dos dados na resolução horizontal. A estimativa da ZFP indica uma abertura mínima com base num modelo adotado, e atualizável. A divergência esférica é uma função suave, e tem base física para ser usada na definição da matriz ponderação dos dados em métodos de inversão tomográfica. A necessidade de robustez na inversão pode ser analisada em seções sísmicas (FC, PMC) submetida a filtragens (freqüências de cantos: 5;15;75;85; banda-passante trapezoidal), onde se pode identificar, comparar e interpretar as informações contidas. A partir das seções, concluímos que os dados são contaminados com pontos isolados, o que propõe métodos na classe dos considerados robustos, tendo-se como referência a norma 2 (quadrados- mínimos) de ajuste de curvas. Os algoritmos foram desenvolvidos na linguagem de programação FORTRAN 90/95, usando o programa MATLAB para apresentação de resultados, e o sistema CWP/SU para modelagem sísmica sintética, marcação de eventos e apresentação de resultados. / The present Masters Thesis had for objective the study of the seismic inversion problem based on flat reflectors for common-source (SC) and common- mid-point (CMP) gathers. The forward model is described by homogeneous, isotropic, plane horizontal layers. The problem is related to the NMO stack based on the optimization of the semblance function, for CMP sections corrected for moveout time. The study was based on two principles. The first principle adopted was the combination of two groups of inversion methods: A global and a local method. The second principle adopted was stripping according to the Wichert-Herglotz-Bateman theory, that establishes that to know a lower layer it is necessary to know first the upper layer. The application of the study aims at the seismic simulation of the terrestrial Solimões and marine sedimentary basins to obtain a 1D distribution of velocity and layer thicknesses of the subsurface of target horizons. In this sense, we limited the inversion experiments to 4 to 11 reflectors, once in practice the industry limits the interpretations to be to about one same number of 3 to 4 main reflectors. Stands out that this model is applicable as initial condition to the imaging of seismic sections in geologically complex regions wit h slow lateral variation of velocities. The synthetic data was produced based on geological information that corresponds to strong a priori information for the inversion model. For the construction of models related to the projects in progress, we analyzed the following relevant subjects: (1) Geology of sedimentary basins terrestrial Solimões and marine (stratigraphy, structural, tectonics and petroliferous); (2) Physics of the vertical and horizontal seismic resolution; and (3) Temporal-spatial discretization of the multi-coverage cube. The inversion process is dependent on the discretization of the wave field in time-space, on the physical parameters of the seismic survey, and of further on the resampling in the multiple coverage cube. The direct model us ed corresponds to the case of the NMO (1D) stack operator, considering a flat observation topography. The basic criterion taken as reference for the inversion and curve fit is the norm 2 (quadratic). The inversion using the simple present model is computational attractive for being fast, and convenient for allowing several other techniques be included with a logical physical interpretation; e.g., the Fresnel projected zone (ZFP), the direct calculation of the spherical divergence, Dix inversion, linear inversion by reparametrization, a priori information, and regularization. The ZFP shows to be a useful concept to establish the aperture of the spatial inversion window in the time-distance section. The ZFP represents the influence of the data in the horizontal resolution. The estimative of the ZFP indicates a minimum aperture based on an adopted model. The spherical divergence is a smooth function, and it has physical basis to be used in the definition of a data weight matrix for tomographic inversion methods. The necessity of robustness in the inversion can be analyzed in seismic sections (CS and CMP) submitted to filtering (corners frequencies: 5; 15; 75; 85; pass-band trapezoidal), where one can identify, compare and interpret the information contained. From the sections, we conclude that the data are contaminated with isolated points, what proposes methods in the class considered as robust having as reference the norm 2 (least-square) of curve fitting. The development of the algorithms used the FORTRAN 90/95 programming language, the program MATLAB for presentation of results; and the package CWP/SU for synthetic seismic modeling, picking of events and presentation of results.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-23T13:59:49Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoGravimetricaRegularizacao.pdf: 2505565 bytes, checksum: 49a8f352bd30f36628c17ec58229301c (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-30T13:08:11Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoGravimetricaRegularizacao.pdf: 2505565 bytes, checksum: 49a8f352bd30f36628c17ec58229301c (MD5) / Made available in DSpace on 2014-09-30T13:08:11Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoGravimetricaRegularizacao.pdf: 2505565 bytes, checksum: 49a8f352bd30f36628c17ec58229301c (MD5) Previous issue date: 2005 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Apresentamos um novo método para inversão gravimétrica da distribuição espacial do contraste de densidade no plano horizontal, baseado na combinação da maximização da entropia de ordem zero com a minimização da entropia de ordem um. O topo e a base das fontes gravimétricas são presumidos ser planos e horizontais e o modelo interpretativo consiste de uma malha de prismas justapostos em ambas às direções horizontais, sendo os contrastes de densidade de cada prisma os parâmetros a serem estimados. A maximização da entropia de ordem zero é similar ao vínculo de suavidade global, enquanto a minimização da entropia de ordem um favorece descontinuidades na distribuição do contraste de densidade. Conseqüentemente a combinação judiciosa de ambas pode levar a soluções apresentando regiões com contrastes de densidade virtualmente constantes (no caso de corpos homogêneos), separadas por descontinuidades abruptas. O método foi aplicado a dados sintéticos simulando a presença de corpos intrusivos em sedimentos. A comparação dos resultados com aqueles obtidos através do método da suavidade global mostra que ambos os métodos localizam as fontes igualmente bem, mas o delineamento de seus contornos é efetuado com maior resolução pela regularização entrópica, mesmo no caso de fontes com 100 m de largura separadas entre si por uma distância de 50 m. No caso em que o topo da fonte causadora não é plano nem horizontal, tanto a regularização entrópica como a suavidade global produzem resultados semelhantes. A metodologia apresentada, bem como a suavidade global foram aplicadas a dois conjuntos de dados reais produzidos por intrusões em rochas metamórficas. O primeiro é proveniente da região de Matsitama, no nordeste de Botswana, centro sul da África. A aplicação das duas metodologias a estes dados produziu resultados similares, indicando que o topo das fontes não é plano nem horizontal. O segundo conjunto provém da região da Cornuália, Inglaterra e produziu uma distribuição estimada de contraste de densidade virtualmente constante para a regularização entrópica e oscilante para a suavidade global, indicando que a fonte gravimétrica apresenta topo aproximadamente plano e horizontal, o que é confirmado pela informação geológica disponível. / We present a new gravity inversion method, which produces an apparent density contrast mapping on the horizontal plane by combining the minimization of the first-order entropy with the maximization of the zero-order entropy of the estimated density contrasts. The interpretation model consists of a grid of vertical, juxtaposed prisms in both horizontal directions. The top and the bottom of the gravity sources are assumed to be flat and horizontal, and the parameters to be estimated are the prism density contrasts. The maximization of the zero-order entropy is similar to the global smoothness constraint whereas the minimization of the first-order entropy favors solutions presenting sharp borders, so a judicious combination of both constrains may lead to solutions characterized by regions where the estimated density contrasts are virtually constant (in the case of homogeneous bodies), separated by sharp discontinuities. The method has been applied to synthetic data simulating the presence of intrusive bodies in sediments. The comparison of the results with those obtained with the global smoothness constraint method shows that both methods produce good and equivalent locations of the source positions, but the entropic regularization delineates the contour of the bodies with greater resolution, even in the case of 100 m wide bodies separated by a distance as small as 50 m. In the case that the depth to the top of the causative sources is neither flat nor horizontal, both method produced similar results. Both the proposed and the global smoothness methods have been applied to two sets of real data produced by igneous intrusions into metamorphic rocks. The first one is from the Matsitama region, northeast of Botswana. The application of both methods to this data set lead to similar results, indicating that depth to the sources top is neither flat nor horizontal. The second set is from the Cornuália region, England. The results produced by both methods in this case are reasonably different, indicating that the gravity source present an approximately flat and horizontal top. This is confirmed by the available geological information.

A caracterização geológica e geotécnica dos locais de implantação de aterros sanitários é fundamental, tanto para fornecer subsídios adequados ao projeto de instalação das células de deposição dos resíduos sólidos urbanos (RSU), quanto para que se conheçam as características da área antes de ser modificada pela obra. Neste sentido, o presente trabalho visou mostrar a contribuição que a geofísica, com o uso da técnica da eletrorresistividade e do potencial espontâneo, associado a informações de investigação geotécnica direta pode dar para o estabelecimento de um background da área de instalação do aterro sanitário da cidade de São Carlos-SP. Foram realizados 7 caminhamentos elétricos (CE) com o arranjo dipolo-dipolo com eletrodos espaçados de 10 metros e três sondagens elétricas verticais (SEV), utilizando o arranjo Schlumberger, com profundidade de investigação de até 100 metros. Para processamento dos modelos bidimensionais dos CE utilizou-se o programa RES2DINV, e enquanto que para as SEV foi utilizado o modelo unidimensional. A interpretação dos modelos geoelétricos foi aperfeiçoada pela confrontação com informações geológico-geotécnicas prévias da área, obtidas nos estudos de impacto ambiental (EIA/RIMA) e resultados de sondagens de simples reconhecimento (SPT), realizadas na fase de estudo para implantação do aterro sanitário. Com isso foi possível definir com clareza a profundidade do topo rochoso e das zonas de maior umidade nos CE e SEV. O conjunto de seções de CE foi georeferenciada com o auxílio de GPS de precisão, o que permitiu o tratamento espacial das informações e o estabelecimento dos mapas da superfície do topo rochoso, das zonas de fluxo de água subterrânea e da espessura dos materiais inconsolidados. O conjunto de informações obtido permitiu identificar locais possíveis para a instalação futura de poços de monitoramento, assim como os dados geoelétricos das várias camadas servirá como padrão para identificar eventuais contaminações do solo ou da água subterrânea no futuro. A identificação da presença do topo rochoso em pequena profundidade, em algumas áreas do aterro, permite o melhor planejamento das escavações para instalação das células de deposição do RSU. / The geological and geotechnical characterization of landfill implantation areas is crucial; even to provide appropriate subsidies of urban solid residues (USR) installation project of deposition cells and for knowing the characteristics of the area before being modified by the work. In this sense, the present work aimed to show the contribution that the geophysics, with the use of electrical resistivity technique and the spontaneous potential, associated with direct geotechnical investigation information can give to the establishment of a background of landfill installation area from São Carlos-SP city. Were performed 7 electrical imaging (EI) with the dipole-dipole array with electrodes spaced by 10 meters and three vertical electrical sounding (VES), using the Schlumberger array with depth research up to 100 meters. For processing of two-dimensional models of the EI was used the RES2DINV program and while for the vertical electrical sounding (VES)-was used one-dimensional model. Interpretation of geoelectric models was enhanced by the confrontation with geological and geotechnical information of prior area, obtained in the environmental impact assessment (EIA) and results of a simple reconnaissance probes, performed in the study phase for the landfill implementation. It was possible to clearly define the depth of the bedrock and the higher humidity areas in the EI and VES. The set of sections electrical Imaging (EI) was georeferenced with the aid of GPS accuracy that allowed the spatial treatment of Information and the setting maps of the bedrock surface, the groundwater flow zones and the thickness of the unconsolidated materials. The obtained set of information enabling identification of possible locations for the wells future monitoring installation, as well as the geoelectrical data from several layers will be use as a standard to identify possible contamination of the soil or groundwater in the future. The identification of bedrock presence in small depth, in some areas of landfill, allows better planning of excavation for installation of MSW deposition cells.

Orientador: Emilson Pereira Leite / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Geociências / Made available in DSpace on 2018-08-27T18:29:04Z (GMT). No. of bitstreams: 1 Cataldo_RafaelAmaral_M.pdf: 3826932 bytes, checksum: 2389a2a5e78d7644118f034812b852e6 (MD5) Previous issue date: 2015 / Resumo: Esta pesquisa mostra a aplicação do método de inversão simultânea em dados sísmicos pré-empilhamento em um campo de petróleo localizado na Bacia de Campos, Brasil. Este método tem sido utilizado com sucesso em diversos reservatórios areníticos em várias regiões do mundo. O objetivo desta pesquisa foi o de avaliar e estudar a aplicação deste método em rochas carbonáticas albianas da Formação Quissamã. Com base em descobertas em campos próximos, esta Formação tem o potencial de possuir acumulações de hidrocarbonetos que podem ser revelados com o auxílio de modelos de inversão. Curvas de velocidades de ondas P (Vp), S (Vs) e densidade (?) ao longo de sete poços serviram como parâmetros de entrada para a construção de modelos iniciais tridimensionais. Os poços que não continham curvas de Vs e densidade foram calculados através de equações estabelecidas junto a literatura disponível. Assim, foram realizadas análises dos resultados da inversão ao longo dos poços. Com o intuito de criar os modelos iniciais de inversão, foram usados dados de impedâncias das ondas P (Zp), S (Zs) e ? ao longo dos poços. Os modelos tridimensionais finais foram calculados a partir dos desvios das relações lineares entre os logaritmos de Zs e Zp, assim como entre ? e Zp, os quais foram combinados com os modelos iniciais. Foram testados e analisados diversos parâmetros e filtros para a criação dos modelos iniciais e finais. As interpretações dos resultados da inversão foram conduzidas com base nos padrões encontrados em perfis de resistividade, raios gama, densidade, porosidade neutrão e sônico. Estas análises forneceram os critérios para selecionar os melhores modelos finais que são discutidos em detalhe. Os modelos finais mostram anomalias de baixa impedância consistentes com os estudos realizados junto aos perfis dos poços. Além disso, de acordo com as interpretações, um poço em particular possui alto potencial para conter hidrocarbonetos. O padrão de distribuição de impedâncias neste poço foi utilizado como referência para a identificação de outras regiões com alto potencial, dentro de todo o volume sísmico pré-empilhado disponível / Abstract: We have applied a method of simultaneous inversion of prestack seismic data of a hydrocarbon field located in the Campos Basin, Brazil. This method has been successfully applied to several sandstone reservoirs around the world. The purpose of this research was to test and study the application of this method to Albian carbonate rocks of the Quissamã Formation. Based on what is found in near fields, this Formation potentially contains hydrocarbon accumulations that could be revealed with the help of the inverted models. P-wave velocity (Vp), S-wave velocity (Vs) and density curves along seven wells were used as inputs to construct 3D initial models of acoustic impedance (Zp), shear impedance (Zs) and density (?). Wells without Vs and ? curves had those values calculated from well-established equations available in the literature. Analysis of the inversion results along wells were carried out. Final 3D models were calculated from deviations in linear relationships between the logarithms of Zs and Zp, as well as between density and Zp, which were merged with the initial models. We have tested and analyzed several parameters and filters to create initial and final models. Interpretations of the inversion results were conducted based on patterns found in logs such as resistivity, gamma ray, density, neutron porosity and sonic along the wells. These analyses provided criteria to select the best final models that are discussed in detail. The inverted models show low impedance anomalies that are consistent with previous studies performed with the well logs available. Furthermore, one well in particular was interpreted as having high potential to contain hydrocarbons. This well shows an impedance pattern that allowed us to highlight other areas with the same pattern throughout the entire prestack seismic volume / Mestrado / Geologia e Recursos Naturais / Mestre em Geociências

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-24T12:30:25Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoGravimetricaRelevoBacias.pdf: 5674544 bytes, checksum: c47719a6af47ec5688d59a7923b52522 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-24T15:37:25Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoGravimetricaRelevoBacias.pdf: 5674544 bytes, checksum: c47719a6af47ec5688d59a7923b52522 (MD5) / Made available in DSpace on 2014-09-24T15:37:25Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoGravimetricaRelevoBacias.pdf: 5674544 bytes, checksum: c47719a6af47ec5688d59a7923b52522 (MD5) Previous issue date: 2009 / Apresentamos um método de inversão de dados gravimétricos para a reconstrução do relevo descontínuo do embasamento de bacias sedimentares, nas quais o contraste de densidade entre o pacote sedimentar e o embasamento são conhecidos a priori podendo apresentar-se constante, ou decrescer monotonicamente com a profundidade. A solução é estabilizada usando o funcional variação total (VT), o qual não penaliza variações abruptas nas soluções. Comparamos o métodoproposto com os métodos da suavidade global (SG), suavidade ponderada (SP) e regularização entrópica (RE) usando dados sintéticos produzidos por bacias 2D e 3D apresentando relevos descontínuos do embasamento. As soluções obtidas com o método proposto foram melhores do que aquelas obtidas com a SG e similares às produzidas pela SP e RE. Por outro lado, diferentemente da SP, o método proposto não necessita do conhecimento a priori sobre a profundidade máxima do embasamento. Comparado com a RE, o método VT é operacionalmente mais simples e requer a especificação de apenas um parâmetro de regularização. Os métodos VT, SG e SP foram aplicados, também, às seguintes áreas: Ponte do Poema (UFPA), Steptoe Valley (Nevada, Estados Unidos), Graben de San Jacinto (Califórnia, Estados Unidos) e Büyük Menderes (Turquia). A maioria destas áreas são caracterizadas pela presença de falhas com alto ângulo. Em todos os casos, a VT produziu estimativas para a topografia do embasamento apresentando descontinuidades bruscas e com alto ângulo, em concordância com a configuração tectônica das áreas em questão. / We present a gravity inversion method to reconstruct the discontinuous basement relief of a sedimentary basin, whose density contrast between the basement and the sediments is constant, or decreases monotonically with depth, and known. The solution is stabilized using the total variation functional (VT), which does not impose smoothness on the solution. We compare the proposed method with the global smoothness (SG), weighted smoothness (SP) and entropic regularization (RE) methods using synthetic data produced by 2D and 3D basins presenting discontinuous relief. The solutions obtained with the proposed method were better than those obtained with SG and similar to the ones produced by SP and RE. Differently from SP, however, the proposed method does not require a priori knowledge about the basin maximum depth. As compared with RE, the present method is operationally simpler and requires the specification of just one regularization parameter. The VT, SG, and SP methods were also applied to the following areas: Ponte do Poema (UFPA), Steptoe Valley (Nevada, United States), San Jacinto Graben (California, United States) and Büyük Menderes (Turkey). Almost all those areas are characterized by the presence of highangle faults. In all cases, the VT produced basement relief estimates presenting sharp, high-angle discontinuities, in accordance with the tectonic setting of the areas.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-24T12:16:18Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoGravimetricaRelevo.pdf: 1116307 bytes, checksum: 24c74f0dccc22d0cfa831b6e7d373356 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-24T16:44:14Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoGravimetricaRelevo.pdf: 1116307 bytes, checksum: 24c74f0dccc22d0cfa831b6e7d373356 (MD5) / Made available in DSpace on 2014-09-24T16:44:14Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InversaoGravimetricaRelevo.pdf: 1116307 bytes, checksum: 24c74f0dccc22d0cfa831b6e7d373356 (MD5) Previous issue date: 2007 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Apresentamos um novo método de inversão gravimétrica para estimar o relevo descontínuo do embasamento de bacias sedimentares, cujos pacotes sedimentares podem ter contraste de densidade constante ou variando hiperbolicamente com a profundidade. O método combina a maximização da medida de entropia de ordem zero e a minimização da medida de entropia de primeira ordem do vetor de soluções (profundidades do embasamento). O modelointerpretativo consiste de um conjunto de primas retangulares verticais justapostos, com o contraste de densidade conhecido, cujas espessuras representam as profundidades do embasamento e são os parâmetros a serem determinados. A minimização da entropia de ordem um favorece soluções apresentando descontinuidades abruptas e a maximização da entropia de ordem zero é empregada apenas para evitar sua minimização excessiva. O método foi aplicado a dados sintéticos simulando: (i) bacias intracratônicas com o relevo do embasamento suave ou apresentando descontinuidades abruptas localizadas e (ii) bacias marginais falhadas. No caso de embasamento suave, a comparação dos resultados com aqueles obtidos através da suavidade global e da suavidade ponderada mostrou que ambos os métodos delinearam o embasamento. No caso de embasamentos falhados, seja em bacias intracratônicas ou marginais, as soluções obtidas com a regularização entrópica e a suavidade ponderada foram equivalentes, delineando o embasamento e as descontinuidades, o que não ocorreu com a suavidade global. A regularização entrópica, no entanto, não requereu informação a priori sobre a profundidade, diferentemente da suavidade ponderada. Tanto o método proposto como os métodos das suavidades global e ponderada foram aplicadas a quatro conjuntos de dados reais. O primeiro é proveniente da Ponte do Poema que está localizada no campus da Universidade Federal do Pará, Belém, Pará, permitiu certificar a aplicabilidade prática do método, uma vez que o relevo abaixo da ponte é acessível. O segundo, da porção norte do Steptoe Valley, Nevada, Estados Unidos. As soluções obtidas com a regularização entrópica e com a suavidade ponderada mostram um embasamento apresentando diversas descontinuidades verticais, o que não aconteceu com a suavidade global. O terceiro conjunto provém do graben de San Jacinto, situado na Califórnia, Estados Unidos, que levou a soluções em que o graben é assimétrico e apenas a suavidade ponderada apresentou solução com fundo achatado. O quarto conjunto de dados é oriundo do graben do Büyük Menderes, localizado na Turquia ocidental. As soluções do método proposto e da suavidade ponderada apresentaram grandes deslocamentos verticais, confirmados pela geologia, o que não aconteceu com a suavidade global. / We present a new gravity inversion method, which estimates the basement relief of a sedimentary basin, whose sediments may present a constant or a decreasing density contrast with depth relative to the basement. The method combines the minimization of the first-order entropy measure with the maximization of the zeroth-order entropy measure of the solution vector (depths to the basement). The interpretation model consists of a set of rectangular juxtaposed vertical prisms with known density contrasts and whose thicknesses represent the depths to the basement relief and are the parameters to be estimated. The minimization of the first-order entropy favors solutions presenting abrupt discontinuities, and the maximization of the zeroth-order entropy is employed just to prevent its excessive minimization. We applied our method to synthetic data simulating: (i) intracratonic basins with smooth or discontinuous relief, and (ii) faulted marginal basins. In the case of a smooth basement relief, it was well delineated by the entropic regularization, and by the global and weighted smoothness. In the case of a faulted basement, either in intracratonic or marginal basins, the entropic regularization and the weighted smoothness delineated the basement relief with good precision, differently from the global smoothness. The entropic regularization, however, did not require the knowledge of the basin’s maximum depth. We applied our method to four sets of real Bouguer anomalies. The first one comes from a profile across a bridge located in the campus of the Federal University of Pará, Belém. The solution delineated the known, discontinuous topography below the bridge. The second set comes from the northern portion of Steptoe Valley, Nevada. The solution delineated a discontinuous basement relief in accordance with the known geological setting of the area. The third set comes from the San Jacinto graben, situated in California, United States. The solution indicated an asymmetric graben, as already reported by other authors. The fourth data set comes from the Büyük Menderes graben, western Turkey. The solutions presented large vertical displacements on the northern border, which is confirmed by geological information.

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.