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 do gráfico de Hingle através de rede competitiva angular

A saturação de água é a principal propriedade petrofísica para a avaliação de reservatórios de hidrocarbonetos, pois através da análise dos seus valores é definida a destinação final do poço recém perfurado, como produtor ou poço seco. O cálculo da saturação de água para as formações limpas é, comumente, realizado a partir da equação de Archie, que envolve a determinação da resistividade da zona virgem, obtida a partir de um perfil de resistividade profunda e o cálculo de porosidade da rocha, obtida a partir dos perfis de porosidade. A equação de Archie envolve ainda, a determinação da resistividade da água de formação, que normalmente necessita de definição local e correção para a profundidade da formação e da adoção de valores convenientes para os coeficientes de Archie. Um dos métodos mais tradicionais da geofísica de poço para o cálculo da saturação de água é o método de Hingle, particularmente útil nas situações de desconhecimento da resistividade da água de formação. O método de Hingle estabelece uma forma linear para a equação de Archie, a partir dos perfis de resistividade e porosidade e a representa na forma gráfica, como a reta da água ou dos pontos, no gráfico de Hingle, com saturação de água unitária e o valor da resistividade da água de formação é obtido a partir da inclinação da reta da água. Independente do desenvolvimento tecnológico das ferramentas de perfilagem e dos computadores digitais, o geofísico, ainda hoje, se vê obrigado a realizar a interpretação de ábacos ou gráficos, sujeito a ocorrência de erros derivados da sua acuidade visual. Com o objetivo de mitigar a ocorrência deste tipo de erro e produzir uma primeira aproximação para a saturação de água em tempo real de perfilagem do poço, insere-se o trabalho apresentado nesta dissertação, com a utilização de uma conveniente arquitetura de rede neural artificial, a rede competitiva angular, capaz de identificar a localização da reta da água, a partir da identificação de padrões angulares presentes nos dados dos perfis de porosidade e resistividade representados no gráfico de Hingle. A avaliação desta metodologia é realizada sobre dados sintéticos, que satisfazem integralmente a equação de Archie, e sobre dados reais.

Interpretação gravimétrica-magnética da zona de fratura dupla Bode Verde

Em Janeiro de 1980, foi realizada a Operação Bode Verde I do Projeto CENTRATLAN, que constou da coleta de dados geofísicos e geológicos na área central da Cadeia Mesoatlântica (CMA), a sul da Ilha de Ascensão. Esta operação fez parte de um convênio firmado entre a Marinha Americana e a Marinha Brasileira, cujo objetivo era investigar o Oceano Atlântico Sul, especialmente a área onde se estende a província fisiográfica da Cadeia Mesoatlântica. Com o objetivo de estudar a Zona de Fratura Dupla Bode Verde, e assim conhecer melhor esta feição tão marcante no fundo oceânico e margem continental, foram tratados e interpretados os dados batimétricos, gravimétricos e magnéticos colhidos na operação anteriormente citada. A análise dos dados batimétricos permitiu uma definição para esta zona de fratura, com o seu caráter duplo marcado pela presença de duas calhas contínuas, paralelas entre si, e separadas por um alto que tem o seu próprio "rift-valley". As anomalias gravimétricas encontradas confirmam o caráter duplo, com os dois mínimos observados correspondendo às duas calhas da fratura. Estas calhas separam um bloco crustal de aproximadamente 40 km de largura. Com base nos dados gravimétricos, foram construídas nove seções crustais sobre a zona de fratura, que permitiram observar que a interface crosta-manto sofre um afinamento crustal embaixo das paredes da fratura. Esta interface, que encontra-se normalmente a uma profundidade de 8-9 km, passa a atingir profundidades de 5.5-6.0 km abaixo das paredes da fratura, com a crosta apresentando uma espessura de 2.0 km nestas partes. Este afinamento crustal é causado pela subida de material do manto. A extensiva alteração hidrotermal que ocorre na depressão central da fratura, pode ser a responsável pelo menor afinamento crustal observado nesta parte da zona de fratura. Para a interpretação dos dados magnéticos, foi usada uma escala de tempo de polaridade magnética entre o Cretáceo Inferior e o Cenozóico, e uma taxa de espalhamento oceânico de 2.0 cm/ano. Para a camada de basalto, que é responsável por parte das anomalias magnéticas observadas, foi usada uma espessura média de 0.5 km. Nas paredes da zona de fratura, ocorre uma diminuição na espessura desta camada, havendo uma interrupção da mesma na depressão central da fratura.

Interpretação de perfis dos carbonatos fraturados da Bacia do Pará-Maranhão

A finalidade deste trabalho é apresentar um estudo de caso dos carbonatos terciários da bacia do Pará-Maranhão, do ponto de vista da interpretação dos perfis registrados nessa área. Dois poços-chave, X e Y, foram escolhidos para o estudo. O tratamento dos dados foi realizado utilizando os recursos do programa LOGCALC, instalado no computador IBM-3090, da Petrobrás, e, também, o sistema DLPS, instalado no VAX-8600, da Universidade Federal de Pará. A avaliação da porosidade e, principalmente, das saturações, é dificultada pelas características não-convencionais dos carbonatos. A litologia é complexa, a salinidade da água de formação é relativamente baixa, em torno de 10.000 ppm de NaCl, a densidade das rochas é elevada e os carbonatos estão fraturados. Para diminuir o efeito da composição mineralógica, foi necessário discriminar os diferentes tipos litológicos. Três tipos de carbonatos foram identificados: calcário, calcário arenoso e calcário dolomitizado. A identificação litológica permite maior controle dos parâmetros da matriz e dos expoentes de porosidade das rochas, conduzindo a estimativas de porosidade e de saturações mais confiáveis. A presença de fraturas influencia, marcadamente, a resposta dos perfis, conforme pode ser notado no perfil de identificação de fraturas, nos perfis de resistividade, de densidade e na curva de raios-gama espectral do poço X. O expoente de porosidade, m, tomado do gráfico de Pickett, é frequentemente inferior a 1,5, valor considerado inerente a rochas fraturadas. Os modelos de Rasmus (1983) e de Porter et al.(1969), foram testados para calcular as saturações de água e óleo. A equação tradicional de Archie (1942), com a e m apropriados, também foi aplicada neste trabalho. O esquema para estimativa de saturações que conduz a resultados mais coerentes com os dados dos testes de formação, nos intervalos fraturados, é a saturação de Rasmus, calculada com a relação de Archie aplicada com o expoente de porosidade variável, tomado do modelo de Rasmus.

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 conjunta de dados gravimétricos magnéticos e resistivo

Apesar das grandes vantagens decorrentes da interpretação conjunta de dados geofísicos, a aplicação da inversão simultânea destes dados tem sido pouco estudada, principalmente a nível de simulação. Este trabalho foi desenvolvido com intuito de cobrir parte desta deficiência. Neste trabalho, foi utilizado o método dos mínimos quadrados para comparar os resultados obtidos a partir de: a) inversão dos grupos de dados gravimétricos, magnéticos e resistivos tomados separadamente; b) inversão simultânea destes grupos de dados combinados dois a dois; c) inversão simultânea dos três grupos de dados. O trabalho é desenvolvido a partir de dados teóricos onde são simuladas diversas situações geológicas. A comparação dos resultados é efetuada a partir das estimativas dos parâmetros obtidos por cada inversão, pelos desvios padrões de cada parâmetro (inversão gravimetria-magnetometria para o ruído Gaussiano) e pela redução da ambiguidade, manifestada pela dependência das estimativas em relação a aproximação inicial dos parâmetros. Na maioria dos casos estudados as inversões conjuntas dos dados combinados dois a dois apresentam resultados bem superiores àqueles obtidos usando-se apenas um dos grupos de dados isoladamente, seja nas estimativas dos parâmetros, seja na redução da ambiguidade. Por sua vez, a inversão conjunta dos três grupos de dados apresentam resultados semelhantes às inversões dos grupos de dados combinados dois a dois, contudo em alguns casos a inversão conjunta dos três grupos de dados é a única totalmente independente da aproximação inicial. Dados gravimétricos e magnéticos de duas anomalias reais foram invertidos, produzindo sempre curvas estimadas bem ajustadas aos valores observados.

Inversão de dados de sísmica de refração profunda a partir da curva tempo-distância

O trabalho em pauta tem como objetivo o modelamento da crosta, através da inversão de dados de refração sísmica profunda, segundo camadas planas horizontais lateralmente homogêneas, sobre um semi-espaço. O modelo direto é dado pela expressão analítica da curva tempo-distância como uma função que depende da distância fonte-estação e do vetor de parâmetros velocidades e espessuras de cada camada, calculado segundo as trajetórias do raio sísmico, regidas pela Lei de Snell. O cálculo dos tempos de chegada por este procedimento, exige a utilização de um modelo cujas velocidades sejam crescentes com a profundidade, de modo que a ocorrência das camadas de baixa velocidade (CBV) é contornada pela reparametrização do modelo, levando-se em conta o fato de que o topo da CBV funciona apenas como um refletor do raio sísmico, e não como refrator. A metodologia de inversão utilizada tem em vista não só a determinação das soluções possíveis, mas também a realização de uma análise sobre as causas responsáveis pela ambiguidade do problema. A região de pesquisa das prováveis soluções é vinculada segundo limites superiores e inferiores para cada parâmetro procurado, e pelo estabelecimento de limites superiores para os valores de distâncias críticas, calculadas a partir do vetor de parâmetros. O processo de inversão é feito utilizando-se uma técnica de otimização do ajuste de curvas através da busca direta no espaço dos parâmetros, denominado COMPLEX. Esta técnica apresenta a vantagem de poder ser utilizada com qualquer função objeto, e ser bastante prática na obtenção de múltiplas soluções do problema. Devido a curva tempo-distância corresponder ao caso de uma multi-função, o algoritmo foi adaptado de modo a minimizar simultaneamente várias funções objetos, com vínculos nos parâmetros. A inversão é feita de modo a se obter um conjunto de soluções representativas do universo existente. Por sua vez, a análise da ambiguidade é realizada pela análise fatorial modo-Q, através da qual é possível se caracterizar as propriedades comuns existentes no elenco das soluções analisadas. Os testes com dados sintéticos e reais foram feitos tendo como aproximação inicial ao processo de inversão, os valores de velocidades e espessuras calculados diretamente da interpretação visual do sismograma. Para a realização dos primeiros, utilizou-se sismogramas calculados pelo método da refletividade, segundo diferentes modelos. Por sua vez, os testes com dados reais foram realizados utilizando-se dados extraídos de um dos sismogramas coletados pelo projeto Lithospheric Seismic Profile in Britain (LISPB), na região norte da Grã-Bretanha. Em todos os testes foi verificado que a geometria do modelo possui um maior peso na ambiguidade do problema, enquanto os parâmetros físicos apresentam apenas suaves variações, no conjunto das soluções obtidas.

Inversão de dados sísmicos de reflexão a partir da curva do tempo de trânsito

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.

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.

Inversão gravimétrica usando regularização entrópica

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.

Inversão gravimétrica do relevo do embasamento usando regularização entrópica

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.

Inversão gravimétrica do relevo de bacias extencionais através da variação total

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.

Mapeamento de uma interface separando dois meios homogêneos através da reconstituição gravimétrica compacta: aplicação a dados do LEPLAC-IV

O problema do mapeamento de uma interface separando dois meios homogêneos tem sido tradicionalmente resolvido através da continuação para baixo da anomalia gravimétrica observada. Este procedimento requer a utilização de um filtro passa-baixa ou de um parâmetro de amortecimento, que tendem a diminuir a resolução do relevo estimado da interface. Além disso, o uso da continuação para baixo implica supor (desnecessariamente) que a interface é uma superfície harmônica. Desta forma, a utilização da continuação para baixo é restrita a interfaces que apresentam um relevo suave relativamente à sua profundidade média. As restrições impostas pela continuação para baixo do campo gravimétrico são reduzidas na técnica de inversão proposta para este estudo, cujo objetivo é mapear a interface crosta-manto baseada na incorporação de informações a priori sobre o mínimo momento de inércia da fonte causadora, em relação a um eixo de concentração de massa pré-estabelecido e coincidente com o topo da fonte anômala. Além da introdução de informações sobre o mínimo momento, foi incluída também a possibilidade de haver variação lateral e em profundidade de densidade em setores previamente estabelecidos, o que permite estender a aplicação desta técnica a situações geológicas mais complexas que envolvam a reconstituição de corpos justapostos, que sejam de diferentes naturezas, tais como aqueles presentes em zonas de transição crustal. O método de inversão apresentado foi avaliado em testes usando dados sintéticos e os resultados obtidos mostraram a eficiência desta técnica em recuperar corpos alongados horizontalmente, que apresentem como topo ou base, interfaces arbitrárias a serem mapeadas. Foi também constatada a eficácia do método em reconstituir corpos que comportem variação lateral de densidade que pudessem retratar situações em que se deseja investigar a natureza do embasamento em regiões de transição crustal. Com base nos resultados destes testes, a técnica foi então empregada na inversão gravimétrica de cinco perfis do LEPLAC IV que cruzam a região oceânica adjacente à margem continental sudeste brasileira. Os resultados desta aplicação permitiram inferir a natureza dos corpos que constituem o arcabouço estrutural da região onde encontram-se inseridas a Dorsal de São Paulo e a Zona de Fratura de Porto Alegre, assim como mapear a Moho subjacente.