Avaliação de dados de altimetria do satélite Geosat a partir de comparação com dados do levantamento gravimétrico marinho Equant

O objetivo principal deste trabalho é avaliar o uso dos dados de altimetria de satélites para mapear a superfície do potencial gravitacional (geóide) no mar. Esta avaliação se faz por comparações da resolução e precisão entre os dados de altimetria processados numa superfície equipotencial (o mar) e dados obtidos a partir de levantamentos convencionais. Uma vez processada a superfície equipotencial, quantidades tais como a anomalia "ar livre" juntamente com o desvio vertical podem ser calculadas. Os dados altimétricos ("altura do mar") utilizados neste trabalho foram coletados pelo satélite GEOSAT. Este satélite rastreou diversas áreas oceânicas do globo processando 44 ciclos em dois anos. Alguns pesquisadores utilizaram os valores médios da "altura do mar" deste satélite para melhoramentos em precisão e resolução dos registros. Estes valores tratados estão disponíveis em NOAA (National Oceanic and Atmospheric Administration) sendo deste modo repassados à UFPa para utilização nesta tese. Os dados de gravimetria marinha utilizados neste trabalho são aqueles obtidos do levantamento "Equatorial Atlantic" (EQUANT I e II) resultantes de uma pesquisa conjunta entre várias instituições com objetivos científicos de conhecer o comportamento da margem equatorial Brasileira. Para comparação e integração entre os dois tipos de dados obtidos através de fontes distintas (medidas de satélite e do navio), poder-se-ia obter a aceleração vertical numa superfície equipotencial partindo-se de um tratamento algébrico dos dados coletados por rastreamento altimétrico do satélite Geosat ou alternativamente poder-se-ia processar transformações dos dados de gravimetria marinha em uma superfície equipotencial equivalente. Em decorrência de diferenças no espaçamento entre as linhas dos levantamentos, ou seja, as linhas das trajetórias do satélite estão largamente espaçadas em comparação com aquelas do levantamento marinho, optou-se por transformar os dados gravimétricos do navio em uma superfície equipotencial. Neste tipo de transformação, consideraram-se vários fatores tais como efeitos "aliasing", nível de ruídos nos levantamentos do navio, redução ao geóide (correção "ar livre"), bem como erros computacionais durante as transformações. Com a supressão parcial desses efeitos (enfatizando o "aliasing") encontrou-se forte correlação entre os dois conjuntos de dados verificando-se um nível de coerência satisfatório a partir do comprimento de onda de 11 km. Comparando-se este resultado com o nível de resolução do satélite Geosat largamente estudado por outros pesquisadores enfatiza-se que de fato a resolução dos valores médios (2 anos) do satélite Geosat aproxima-se da resolução dos levantamentos do Equant I e Equant II.

Decomposição em valores singulares aplicada a dados de campo magnético

Neste trabalho, a decomposição em valores singulares (DVS) de uma matriz A, n x m, que representa a anomalia magnética, é vista como um método de filtragem bidimensional de coerência que separa informações correlacionáveis e não correlacionáveis contidas na matriz de dados magnéticos A. O filtro DVS é definido através da expansão da matriz A em autoimagens e valores singulares. Cada autoimagem é dada pelo produto escalar dos vetores de base, autovetores, associados aos problemas de autovalor e autovetor das matrizes de covariância A^TA e AA^T. Este método de filtragem se baseia no fato de que as autoimagens associadas a grandes valores singulares concentram a maior parte da informação correlacionável presente nos dados, enquanto que a parte não correlacionada, presumidamente constituída de ruídos causados por fontes magnéticas externas, ruídos introduzidos pelo processo de medida, estão concentrados nas autoimagens restantes. Utilizamos este método em diferentes exemplos de dados magnéticos sintéticos. Posteriormente, o método foi aplicado a dados do aerolevantamento feito pela PETROBRÁS no Projeto Carauari-Norte (Bacia do Solimões), para analisarmos a potencialidade deste na identificação, eliminação ou atenuação de ruídos e como um possível método de realçar feições particulares da anomalia geradas por fontes profundas e rasas. Este trabalho apresenta também a possibilidade de introduzir um deslocamento estático ou dinâmico nos perfis magnéticos, com a finalidade de aumentar a correlação (coerência) entre eles, permitindo assim concentrar o máximo possível do sinal correlacionável nas poucas primeiras autoimagens. Outro aspecto muito importante desta expansão da matriz de dados em autoimagens e valores singulares foi o de mostrar, sob o ponto de vista computacional, que a armazenagem dos dados contidos na matriz, que exige uma quantidade n x m de endereços de memória, pode ser diminuída consideravelmente utilizando p autoimagens. Assim o número de endereços de memória cai para p x (n + m + 1), sem alterar a anomalia, na reprodução praticamente perfeita. Dessa forma, concluímos que uma escolha apropriada do número e dos índices das autoimagens usadas na decomposição mostra potencialidade do método no processamento de dados magnéticos.

Delineamento do pé do talude na margem continental do Ceará através da integração de dados geológicos e geofísicos

A metodologia de interpretação integrada dos dados geológicos e geofísicos observados em um perfil da margem continental do Ceará possibilita a identificação e a integração de características peculiares a cada tipo de dado. Dessa forma, é possível se definir a localização mais provável de feições estruturais importantes, tais como a fronteira entre as crostas continental e oceânica e o pé do talude, objeto do presente estudo. Segundo o Artigo 76 (parágrafo 4, item b) da Convenção das Nações Unidas para o Direito do Mar, o pé do talude é definido como o ponto de variação máxima do gradiente do talude na sua base. Entretanto, essa definição, apesar de simples no contexto fisiográfico, não é suficiente para propiciar a localização do pé do talude como preconiza a Convenção, razão pela qual se aplicam os métodos geofísicos. Dentro do contexto geofísico-geológico está implícita a interpretação quantitativa das anomalias gravimétricas ar-livre, que possibilita o delineamento do modelo geofísico representando a subsuperfície, cuja finalidade é subsidiar geologicamente a interpretação integrada dos dados aludidos. Um procedimento automático de ajuste de curvas combinando as técnicas de inversão de busca sistemática e a que utiliza derivadas foi usado com o propósito de gerar o modelo geofísico. A aplicação rigorosa de vínculos preliminarmente e a constante reavaliação desses vínculos através de um processo interativo entre a sísmica e a gravimetria, gerado durante a interpretação quantitativa das anomalias ar-livre, possibilitaram que o modelo geofísico final estivesse dentro dos padrões geológicos para área, notadamente quanto ao equilíbrio isostático (Teoria de Airy). O objetivo do presente trabalho é se estudar as características geológicas e geofísicas observadas ao longo de um perfil da margem continental do Ceará (LEPLAC III), notadamente quanto ao pé do talude, buscando estabelecer a aplicabilidade de uma metodologia de interpretação integrada desses dados, cuja finalidade é se definir de forma sistemática a localização mais provável para esta feição fisiográfica. A metodologia de interpretação integrada dos dados geológicos e geofísicos empregada mostrou-se eficiente para este objetivo. Foi possível se integrar: (i) a localização fisiográfica (distância da costa e profundidade) do pé do talude; (ii) a zona de instabilidade tectônica evidenciada pelos falhamentos, comuns nesta região; (iii) ao fim de uma zona magnética perturbada, associado a um ponto de mínimo na curva de anomalia magnética, e que possivelmente delimita o início de uma zona magnética quieta, denominada de anomalia E e (iv) um ponto de inflexão na curva de anomalia ar-livre, associado ao efeito gravimétrico do contraste de densidades entre as crostas continental, os sedimentos e a água do mar, evidenciado pela geometria do talude. Foi possível ainda se definir a localização mais provável para a fronteira entre as crostas continental e oceânica. Dada a rigorosidade na aplicação das técnicas de inversão e dos vínculos é provável que as correlações das características intrínsecas a cada tipo de dado efetuadas na conclusão desse trabalho tenham fundamento e possam ser confirmadas. A condição para isto é a aplicação da metodologia aqui estabelecida em um número maior de perfis.

Os efeitos Delaware e Groningen: um estudo quantitativo por elementos finitos

Os efeitos Delaware e Groningen são dois tipos de anomalia que afetam ferramentas de eletrodos para perfilagem de resistividade. Ambos os efeitos ocorrem quando há uma camada muito resistiva, como anidrita ou halita, acima do(s) reservatório(s), produzindo um gradiente de resistividade muito similar ao produzido por um contato óleo-água. Os erros de interpretação produzidos têm ocasionado prejuízos consideráveis à indústria de petróleo. A PETROBRÁS, em particular, tem enfrentado problemas ocasionados pelo efeito Groningen sobre perfis obtidos em bacias paleozóicas da região norte do Brasil. Neste trabalho adaptamos, com avanços, uma metodologia desenvolvida por LOVELL (1990), baseada na equação de Helmholtz para H_Φ, para modelagem dos efeitos Delaware e Groningen. Solucionamos esta equação por elementos finitos triangulares e retangulares. O sistema linear gerado pelo método de elementos finitos é resolvido por gradiente bi-conjugado pré-condicionado, sendo este pré-condicionador obtido por decomposição LU (Low Up) da matriz de stiffness. As voltagens são calculadas por um algoritmo, mais preciso, recentemente desenvolvido. Os perfis são gerados por um novo algoritmo envolvendo uma sucessiva troca de resistividade de subdomínios. Este procedimento permite obter cada nova matriz de stiffness a partir da anterior pelo cálculo, muito mais rápido, da variação dessa matriz. Este método permite ainda, acelerar a solução iterativa pelo uso da solução na posição anterior da ferramenta. Finalmente geramos perfis sintéticos afetados por cada um dos efeitos para um modelo da ferramenta Dual Laterolog.

Efeitos do manto de intemperismo sobre anomalias eletromagnéticas provocadas por corpos tabulares inclinados: um estudo através de modelamento reduzido

Electromagnetic methods have been extensively applied in the prospecting of sulphide bodies and other conducting materials. The interpretation of e. m. data is based on the results obtained either with reduced scale or analytical modelling. In most models, the host rocks, the overburden as well as the halo of disseminated sulphides are considered highly resistive although in nature they are often weakly conductive. Presence of a well-developed conducting overburden in tropical and sub-tropical regions and a saline rich crustal layer in semi-arid regions have been found to modify significantly the e. m. anomalies in practice. Therefore, the parameters of the target, determined on the basis of simple models where the presence of the conducting environment is neglected, are found to be in considerable error. The effects of the overburden on the e. m. anomalies of a dipping tabular ore body were studied with reduced scale models for varying response parameters of both the overburden and the orebody, and also for different depths and dips of the target. The overburden and the orebody were represented by metallic sheets of varying thickness in the scale model developed in accordance with the law of electromagnetic similitude. The results of these investigations show that the overburden affects the anomaly by causing: a) phase rotation; b) amplitude reduction; c) base level displacement; d) reversion of the quadrature; and e) appearance of an extra peak in the quadrature in the case of low dipping models. The last two effects complicate the quantification of the anomalies but, on the other hand, provide a qualitative indication of the response parameters of both the ore body and the overburden. The results were assembled in the form of Argand diagrams and, finally, an interpretation scheme is suggested for the e. m. field data on the basis of these diagrams.

Efeitos do manto de intemperismo sobre anomalias “VLF” dos corpos tabulares inclinados

A utilização dos métodos indutivos de propagação E.M. na exploração mineral em regiões tropicais, apresenta grandes dificuldades devido a presença de uma camada superficial condutiva (manto de intemperismo) comumente encontrada nestas regiões. Na região Amazônica, o manto apresenta-se bastante desenvolvido e condutivo, e em regiões semi-áridas, pode-se formar uma fina crosta superficial de sal. Em conseqüência disto, a interpretação dos dados E.M. obtidos para modelos que não consideram uma cobertura condutiva levam a erros consideráveis. Objetivando-se estudar os efeitos do manto sobre anomalias VLF devidas a corpos tabulares inclinados em contato com o manto (manto ohmico), foi realizada uma série de experimentos através do modelamento analógico, considerando-se diferentes parâmetros de resposta para o manto e o corpo. O manto de intemperismo foi simulado por soluções de cloreto de amônia (NH₄Cl) dispostas horizontalmente e o corpo condutor por chapas de grafite colocadas em posições inclinadas verticalmente. Utilizou-se quatro corpos condutores e três mantos com diferentes espessuras e condutividades, simulando, desta forma, diversas situações geológicas. Os resultados são dados por simples situações dos corpos localizados em um meio não condutor (ar), onde os parâmetros variados são: profundidade do topo, condutividade e mergulho do corpo. Os efeitos da condutividade da cobertura são amplamente ilustrados e avaliados. Para a análise dos resultados, foi plotado um conjunto de curvas considerando-se os valores pico-a-pico das anomalias de "tilt angle" e de elipsidade. Os resultados foram sintetizados em um outro conjunto de curvas reunidas em diagramas de Argand. Estando ou não o manto presente, observou-se, tanto para o tilt angle quanto para a elipsidade, o efeito do aumento da profundidade é o de reduzir a magnitude pico-a-pico e a forma do pico da anomalia, fazendo com que este afaste-se do ponto de "cross-over". Para um condutor de mesma espessura, o aumento da condutividade causa um ligeiro aumento nas anomalias de tilt angle, e uma atenuação nas anomalias de elipsidade. O efeito geral na variação do mergulho do condutor é o de causar uma assimetria nos perfis de tilt angle e de elipsidade. O aumento da condutância do manto de intemperismo causa um acréscimo nas anomalias de elipsidade e uma ligeira diminuição nas anomalias de tilt angle; porém, a partir de um certo valor de condutância do manto (mantos mais condutivos) tanto as anomalias de tilt angle quanto as anomalias de elipsidade começam a atenuar. Há rotação de fase no sentido anti-horário, sendo mais intensa para grandes valores de número de indução do corpo. Na presença do manto, o corpo parece estar a uma profundidade inferior à verdadeira e a ser menos condutivo.

Efeitos do strike das estruturas geológicas de duas dimensões nas pseudos-seções de IP resistividade

Os levantamentos de IP-resistividade efetuados na Serra dos Carajás não foram executados ortogonalmente às estruturas geológicas, pois utilizaram linhas anteriormente abertas pelas equipes de geoquímica. Este fato motivou este estudo teórico da influência da direção das linhas de medidas de IP-resistividade em relação ao "strike" da estrutura. Foi usado o programa de elementos finitos de Rijo (1977), desenvolvido para levantamentos perpendiculares às estruturas com as adaptações necessárias. A modificação principal foi na rotina de transformação inversa de Fourier. Para o caso simples dos levantamentos perpendiculares, a transformada inversa é uma integral discreta com apenas sete pontos. No entanto, para as medidas obliquas, o integrando é oscilatório, e portanto, a integral a ser calculada é mais complexa. Foi adaptado um método apresentado por Ting e Luke (1981), usando dezoito pontos em cada integração. Foi constatado que o efeito da direção da linha em relação ao "strike" é desprezível para ângulos maiores que 60 graus. Para ângulos menores, o efeito consiste no alongamento da anomalia, com pequenas alterações em seu centro. Não há uma maneira simples de compensar este efeito com mudanças nos parâmetros do modelo.

Estudo de anomalias eletromagnéticas de um condutor tabular vertical sob camadas parcialmente condutiva em multifrequência e multiseparação através de modelamento analógico

A resposta eletromagnética (EM) de um corpo condutivo envolvido por uma zona parcialmente condutiva, torna-se bastante diferente daquela de um corpo condutivo em um meio altamente resistivo. As zonas parcialmente condutivas, como por exemplo, rocha encaixante, halo de sulfetos disseminados ou manto de intemperismo, que envolvem o corpo condutivo, afetam a resposta EM de diferentes maneiras, dependendo de suas características físicas e geométricas e, em particular, do sistema de prospecção EM utilizado. Neste trabalho em modelamento analógico, foi feita uma análise de anomalias EM provocadas por corpos condutivos tabulares verticais sob manto de intemperismo, em levantamentos terrestres para diferentes sistemas de bobinas - horizontal coplanar, vertical coplanar e vertical coaxial - em oito frequências na faixa de 250 Hz a 35 kHz e separações entre as bobinas de 0,15; 0,20 e 0,25m. O manto de intemperismo foi simulado por folhas de aço finas dispostas horizontalmente e o corpo condutor principal por folhas de alumínio finas colocadas verticalmente. As dimensões das folhas foram determinadas de acordo com as condições de modelamento para o plano e o semi-plano. Foram utilizados três corpos e três mantos com diferentes espessuras e condutividades, simulando, deste modo, diversas situações geológicas. Os resultados mostraram que cada sistema de bobinas é afetado diferentemente pela presença do manto de intemperismo. Para a análise dos resultados foi plotado um conjunto de diagramas considerando os valores pico-a-pico das anomalias em fase e em quadratura. Um outro conjunto de diagramas mostra as amplitudes máximas em fase, que ocorrem quando a componente em quadratura se anula em uma frequência relativamente baixa para um conjunto de corpo-manto, e as amplitudes máximas em quadratura, que ocorrem quando a resposta em fase atinge um mínimo próximo de zero, em frequências relativamente altas. Com isto foi possível conhecer a faixa de frequências para cada sistema de bobinas, onde a resposta EM se encontra o mínimo afetada pela presença do manto de intemperismo. A maior amplitude na resposta é obtida no sistema horizontal coplanar e a menor no sistema vertical coplanar. Um aumento na separação entre as bobinas é acompanhado por um deslocamento da anomalia para baixas frequências. A faixa de frequências, onde a presença do manto tem pouca influência na resposta do corpo condutivo, e maior para o sistema vertical coaxial e menor para o sistema horizontal coplanar. Esses resultados dão uma luz para o conhecimento da posição e da largura da banda de frequências utilizável, assim como as melhores separações entre transmissor-receptor, para auxiliar no planejamento de sistemas de prospecção EM, de modo que a resposta fique o mais livre possível de sinais indesejáveis, tais como os causados pela presença do manto de intemperismo.

Estudo do comportamento das feições lineares observado em dados aeromagnéticos da Bacia do Paraná

O processamento dos dados aeromagnéticos em imagens revelou-se como um bom método de interpretação, que aqui denominamos de método das imagens sombreadas. A aplicação deste método foi feita em dados aeromagnéticos de uma área da Bacia do Paraná, com o objetivo de realçar os lineamentos que não são evidentes nos mapas de contorno. Para subsidiar os nossos estudos, utilizamos imagens de LANDSAT porque muitos lineamentos presentes na bacia não possuem resposta ao levantamento aeromagnético. Com a integração destas duas imagens definimos duas direções dominantes para os lineamentos nos seguintes intervalos: N40° — 60°W e N40° — 65°E. Com a definição deste padrão, obtivemos informações do arcabouço estrutural e, consequentemente, uma idéia da evolução tectônica da bacia. A outra parte desenvolvida neste trabalho foi o cálculo do vetor de magnetização total nos principais lineamentos da bacia, e eles apresentaram um total de 5 inclinações de magnetização diferentes. Essas inclinações foram explicadas pela superposição do efeito da magnetização induzida, efeito de reversão do campo geomagnético, e o efeito da magnetização anisotrópica causada pela desmagnetização nos diques.

Estudos dos efeitos do manto de intemperismo sobre anomalias de VLF geradas por corpos inclinados: modelagem numérica através de programa eGs

As referências apresentadas e utilizadas através de métodos indutivos de propagação eletromagnética E.M., são utilizadas na prospecção de minerais, entre outros, haja vista que em regiões tropicais existem camadas superficiais condutivas (manto intemperismo) que levam a influenciar, evidentemente, na interpretação do resultado real. Em nossa região Amazônica, o manto apresenta-se bastante desenvolvido e condutivo. O objetivo deste trabalho é estudar os efeitos dos mantos sobre anomalias VLF, devidos a corpos tabulares inclinados em contato ou não com o manto (manto galvânico e manto indutivo). Deste modo, foram realizadas umas séries de experimentos em modelagens numéricas através da eGs, considerando-se diferentes parâmetros de resposta para o manto e para o corpo. Os trabalhos de modelagem numérica, bem como seus resultados foram apresentados por diversas situações, como: variação da profundidade do topo do condutor, variação do mergulho do condutor, variação da condutância do manto, entre outros. Para interpretação dos resultados, foi gerado um conjunto de curvas considerando-se valores pico-a-pico das anomalias de Tilt Angle e de Elipsidade. Os resultados foram sintetizados em um outro conjunto de curvas reunidas em diagramas de Argand. A medida em que se aumenta a profundidade, o pico-a-pico e a forma do pico da anomalia do Tilt Angle e da Elipsidade diminuem na presença do manto, ou na ausência do mesmo, fazendo com que este se afaste do ponto de “Cross-Over”. Para um condutor da mesma espessura, o aumento da condutância causa um aumento nas anomalias de Tilt Angle e da Elipsidade. O efeito geral na variação do mergulho do condutor é o de causar uma assimetria nos perfis de Tilt Angle e Elipsidade. O aumento da condutância do manto de intemperismo causa uma diminuição nas anomalias de Tilt Angle e de Elipsidade. Há rotação de fase no sentido anti-horário, sendo mais intensa para grandes valores de condutância do corpo. Na presença do manto, o corpo parece estar a uma profundidade inferior à verdadeira e, tendendo, a ser menos condutivo.

Interpretação de anomalias magnéticas usando integrais finitas dos momentos das componentes

O objetivo do trabalho consiste em desenvolver uma metodologia para determinar os parâmetros de um corpo causador de uma anomalia magnética. Essa metodologia baseia-se na utilização de expressões analíticas, deduzidas para as integrais finitas dos momentos da anomalia de intensidade total e das componentes H_x, H_y e H_z ao longo dos três eixos de um sistema de coordenadas cartesianas. Por meio do ajuste entre essas expressões analíticas e as integrais numéricas finitas dos momentos das componentes magnéticas obtidas a partir de um levantamento geofísico, pode-se computar, através de um processo iterativo, os parâmetros de magnetização, posição, profundidade e dimensões de um corpo anômalo. No caso em que são conhecidas as medidas de somente uma componente magnética, ainda é possível aplicar o método pois as demais componentes podem ser obtidas por um esquema de filtragem matemática. A metodologia foi testada com sucesso para modelos dipolo pontual e linha de dipolos os quais são muito utilizados em Geofísica para interpretar anomalias magnéticas produzidas por corpos geológicos que possuem uma ou duas dimensões horizontais muito menores do que a profundidade.

Interpretação quantitativa de anomalias gravimétricas ar-livre da margem continental norte brasileira, setor nordeste

Anomalias gravimétricas ar-livre de perfis perpendiculares a margem continental do tipo passiva apresentam uma configuração padrão. Esta configuração é, satisfatoriamente, explicada por um modelo geofísico formado por uma distribuição de descontinuidades horizontais bidimensionais. Um processo automático de busca aleatória é proposto para a interpretação quantitativa dos dados. Através do método de poliedros flexíves (Simplex), os parâmetros principais do modelo - o contraste de densidade, a profundidade, o rejeito e a localização de cada descontinuidade, puderam ser encontrados, admitindo uma relação número de pontos/número de parâmetros, a determinar, conveniente. Sobre a região do talude, as anomalias ar-livre da margem continental podem ser explicadas por uma única descontinuidade horizontal (degrau simples); e tendo que a resposta dos dados gravimétricos no domínio do número de onda contém informações sobre esta anomalia, foi proposto um procedimento gráfico iterativo para a análise espectral deste sinal. Aplicando a transformada de Fourier é possível determinar a profundidade e o rejeito da descontinuidade, e conhecendo estes parâmetros a densidade é calculada unicamente. O objetivo básico do uso destes procedimentos seria combinar os dois métodos de interpretação nos domínios do espaço e do número de onda, com a finalidade de obter soluções vinculadas mais plausíveis quanto ao contexto geológico esperado para a área estudada. Os dois procedimentos de interpretação foram aplicados nas anomalias gravimétricas ar-livre da margem continental norte brasileira, setor nordeste, abrangendo os estados do Maranhão ao Rio Grande do Norte. As respectivas capacidade de resolução de cada procedimento foram então analisadas. Demonstrou-se que a inversão realizada diretamente no domínio do espaço é mais favorável na interpretação das anomalias ar-livre, embora o tratamento espectral seja relativamente mais simples.

Inversão gravimétrica estável do relevo do embasamento e da variação da densidade com a profundidade em bacias sedimentares

O presente método postula uma variação hiperbólica para o contraste de densidade de uma bacia sedimentar em função da profundidade, e tem dois objetivos: (1) delinear o relevo do embasamento de uma bacia, conhecendo-se a anomalia gravimétrica, o contraste de densidade na superfície da bacia e o fator de decaimento do contraste de densidade com a profundidade; (2) estimar, além do relevo, o valor do contraste de densidade na superfície de uma bacia sedimentar e o fator de decaimento do contraste de densidade com a profundidade, sendo fornecida a anomalia gravimétrica e a profundidade do embasamento em alguns pontos da bacia. Nos dois casos o modelo interpretativo é um conjunto de prismas retangulares verticais justapostos, cujas espessuras, que são parâmetros a serem estimados, representam a profundidade da interface de separação entre os sedimentos e o embasamento. As soluções obtidas nos dois problemas acima formulados são estáveis devido à incorporação de informações adicionais sobre a suavidade do relevo estimado, e o conhecimento da profundidade do relevo do embasamento em alguns pontos, fornecido por furos de sondagem. O método foi testado em anomalias gravimétricas sintéticas produzidas pela simulação de bacias sedimentares com relevos suaves. Os resultados mostraram relevos com boa resolução e valores estimados do contraste de densidade na superfície da bacia e do fator de decaimento do contraste de densidade com a profundidade, próximos aos verdadeiros, indicando dessa maneira o potencial do método em interpretações gravimétricas de bacias sedimentares. O método foi aplicado à anomalia Bouguer da Bacia do Recôncavo, Brasil, delineando um relevo com um valor para a profundidade máxima de cerca de 6 km, semelhante ao estimado em interpretações sísmicas. As estimativas para o contraste de densidade na superfície da Bacia e o fator de decaimento com a profundidade foram -0,30 g/cm³ e 30 km, respectivamente, produzindo uma estimativa para a compactação máxima dos sedimentos da ordem de 4%.

Inversão gravimétrica estável do relevo da base e da variação da densidade com a profundidade em aterros sanitários

O presente trabalho apresenta um método que utiliza uma função hiperbólica para descrever o decaimento do contraste de densidade com a profundidade em aterros sanitários e lixões, e consiste de duas abordagens: (1) o mapeamento da base de um lixão ou de um aterro sanitário conhecendo-se a anomalia gravimétrica, o contraste de densidade na superfície e o fator de decaimento do contraste de densidade com a profundidade; (2) estimação do contraste de densidade na superfície e o fator de decaimento do contraste de densidade com a profundidade a partir do conhecimento da anomalia gravimétrica e do relevo da base do lixão ou do aterro sanitário; uma variante desta última abordagem utiliza os valores conhecidos de profundidade desses ambientes para gerar a anomalia ajustada e estimar o valor do contraste de densidade na superfície e o fator de decaimento do contraste com a profundidade através da minimização da norma Euclidiana do vetor de resíduos entre as observações e o ajuste. Em ambas abordagens o modelo interpretativo é representado por um conjunto de prismas retangulares justapostos, cujas espessuras são os parâmetros a serem determinados. As soluções encontradas através das duas abordagens são estáveis devido à incorporação de informação sobre a suavidade do relevo estimado. O método foi aplicado a dados sintéticos produzidos pela simulação de aterros sanitários e lixões com até 15 m de profundidades e de relevos suaves. Os resultados das interpretações gravimétricas foram consistentes e revelaram a eficácia do método tanto na abordagem (1) como na abordagem (2). O método também foi aplicado à anomalia Bouguer residual do aterro sanitário da Fazenda Thomas, Indiana, USA, indicando que o contraste de densidade do aterro pode ou ser constante ou apresentar variação máxima de 0,09 g/cm³.

Inversão gravimétrica rápida do relevo do embasamento aplicando o funcional da variação total

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