Aplicação do gravímetro de poço em problemas de avaliação de formação no Brasil

No período de outubro a novembro de 1990, foi realizado na Bacia Sergipe-Alagoas o primeiro levantamento gravimétrico de poço no Brasil, através de um programa com a participação da Universidade Federal do Pará, a PETROBRÁS e o U.S.G.S. (United States Geological Survey). Este levantamento teve como objetivos o teste do equipamento, a resposta dos perfis gravimétricos em relação aos problemas encontrados nas bacias sedimentares brasileiras e a comparação com os dados de densidade obtidos com o perfil de densidade compensada (CDL). Os levantamentos foram realizados em três poços e os dados obtidos passaram por um processo de redução, onde os valores foram transformados para miligals e corrigidos dos efeitos de maré, deriva e terreno, para obter as densidades preliminares. Porém, a região onde foram realizados os levantamentos apresenta além do embasamento muito raso, uma seqüência evaporítica com contraste de densidade grande em relação ao resto do pacote sedimentar, gerando gradientes verticais anômalos. Estes efeitos foram então corrigidos e obtidos os valores finais de densidade. Com base em uma caracterização litológica prévia, os valores de densidade foram então comparados com o perfil CDL. Verificou-se a boa qualidade dos dados, o que permitiu o cálculo de valores de porosidade, mostrando assim possíveis intervalos de interesse à exploração de hidrocarbonetos, tanto no embasamento, reservatório mais importante, quanto no intervalo sedimentar.

Avaliação da influência do eletrojato equatorial na interpretação de dados magnetotelúrico bidimensionais

O campo magnetotelúrico em regiões equatoriais viola a condição de ondas planas por causa de uma fonte fortemente concentrada na direção E-W na ionosfera, denominada eletrojato equatorial. No presente trabalho, procurou-se analisar a resposta magnetotelúrica de fontes que simulam o efeito do eletrojato equatorial. Foram considerados dois tipos de fontes para simular o eletrojato: uma linha infinita de corrente e uma distribuição gaussiana de densidade de corrente em relação a uma das coordenadas horizontais. A resistividade aparente foi obtida através da relação de Cagniard e comparada com os resultados de ondas planas. É mostrada também a comparação entre a fase da impedância na superfície, para os três tipos de fontes (ondas planas, eletrojato gaussiano e linha de corrente). O problema de meios com heterogeneidades laterais foi resolvido em termos de campos secundários, sendo as equações diferenciais solucionadas através da técnica de elementos finitos bidimensionais. Os resultados mostram que o eletrojato tem pouca influência nas respostas (resistividade aparente e fase) de estruturas geológicas rasas. Entretanto, a influência pode ser considerável nas estruturas profundas (maior que 5000 m), principalmente se suas resistividades são altas (maior que 100 Ω.m). Portanto, a influência do eletrojato equatorial deve ser considerada na interpretação de dados magnetotelúricos de bacias sedimentares profundas ou no estudo da crosta terrestre.

Comparação dos filtros de velocidade e do operador WHLP-CRS na atenuação de múltiplas

A motivação geológica deste trabalho reside no imageamento de estruturas de bacias sedimentares da região Amazônica, onde a geração e o acúmulo de hidrocarboneto estão relacionados com a presença de soleiras de diabásio. A motivação sísmica reside no fato de que essas rochas intrusivas possuem grandes contrastes de impedância com a rocha encaixante, o que resulta em múltiplas, externas e internas, com amplitudes semelhantes as das primárias. O sinal sísmico das múltiplas podem predominar sobre o sinal das reflexões primárias oriundas de interfaces mais profundas, o que pode dificultar o processamento, a interpretação e o imageamento da seção sísmica temporal. Neste trabalho, estudamos a atenuação de múltiplas em seções sintéticas fonte-comum (FC) através da comparação de dois métodos. O primeiro método resulta da combinação das técnicas Wiener-Hopf-Levinson de predição (WHLP) e o de empilhamento superfície-de-reflexão-comum (CRS), e denominando WHLP-CRS, onde o operador é desenhado exclusivamente no domínio do tempo-espaço. O segundo método utilizado é o filtro de velocidade (ω-k) aplicado após o empilhamento superfície-de-reflexão (CRS), onde o operador é desenhado exclusivamente no domínio bidimensional de freqüência temporal-espacial. A identificação das múltiplas é feita na seção de afastamento-nulo (AN) simulada com o empilhamento CRS, e utiliza o critério da periodicidade entre primária e suas múltiplas. Os atributos da frente de onda, obtidos através do empilhamento CRS, são utilizados na definição de janelas móveis no domínio tempo-espaço, que são usadas para calcular o operador WHLP-CRS. O cálculo do filtroω-k é realizado no domínio da freqüência temporal-espacial, onde os eventos são selecionados para corte ou passagem. O filtro (ω-k) é classificado como filtro de corte, com alteração de amplitude, mas não de fase, e limites práticos são impostos pela amostragem tempo-espaço. Em termos práticos, concluímos que, para o caso de múltiplas, os eventos separados no domínio x-t não necessariamente se separam no domínio ω-k, o que dificulta o desenho de um operador ω-k semelhante em performance ao operador x-t.

Model of wetland development of the Amapá coast during the late Holocene

Os tipos de vegetação atual, sequências sedimentares, dados de pólen e datações por radiocarbono obtidas em três testemunhos de sedimento da planície costeira de Calçoene foram utilizados para estabelecer uma história paleoecológica durante o Holoceno superior das zonas úmidas costeiras do Amapá conforme as mudanças no regime de inundação, nível do mar e clima. Baseado nestes três registros, quatro fases de desenvolvimento da vegetação são apresentadas e provavelmente refletem a interação entre o fluxo de energia na acumulação do sedimento e a influência das águas salobras e doces na vegetação. Este trabalho sugere alternâncias entre períodos caracterizados por influências marinha e fluvial. O perfil longitudinal não revelou a ocorrência de manguezais nos sedimentos depositados por volta de 2100 anos A.P. Durante a segunda fase, a lama preencheu progressivamente as depressões e canais de maré. Provavelmente, os manguezais iniciaram seu desenvolvimento nas margens dos canais, e os campos herbáceos nos setores elevados. A terceira fase é caracterizada por uma interrupção no desenvolvimento dos manguezais e a expansão da vegetação de várzea devido a uma diminuição na influência das águas marinhas. A última fase é representada pela expansão de manguezais e várzeas. A correlação entre os padrões atuais de distribuição das unidades geobotânicas e a paleovegetação indica que os manguezais e as florestas de várzea estão migrando sobre os campos herbáceos nos setores topograficamente mais elevados do litoral em estudo, o que pode estar relacionado a um aumento do nível relativo do mar.

The ceramic artifacts in archaeological black earth (terra preta) from lower Amazon region, Brazil: Mineralogy

Sítios arqueológicos com Terra Preta, denominados de Terra Preta de Índio ou ainda Terra Preta Arqueológica (TPA) são muito freqüentes na Amazônia. As TPA geralmente contém fragmentos de vasos cerâmicos, por vezes abundantes, além de líticos, que são materiais de grande importância para os estudos arqueológicos. Para consubstanciar esses estudos, realizou-se pesquisas mineralógicas e químicas em fragmentos cerâmicos provenientes de dois sítios arqueológicos da região de Cachoeira-Porteira, Estado do Pará. Os fragmentos foram classificados segundo seus principais temperos em: cauixi, cariapé, areia+feldspatos e caco de vaso cerâmico. Mineralogicamente são compostos de quartzo, minerais de argila calcinados (especialmente caulinita), feldspatos (albita e microclínio), hematita, goethita, maghemita, variscita-estrengita, fosfatos amorfos, anatásio, e raramente apatita, rhabdophana e óxidos de Mn e Ba. Cauixi e cariapé são componentes orgânicos silicosos e amorfos a DRX. A composição mineralógica e a morfologia dos seus grãos indicam saprólito (material argiloso rico em quartzo) derivado de rochas ígneas félsicas de granulação fina ou rochas sedimentares ricas em argilominerais como matéria-prima dos vasos cerâmicos. Neste material argiloso cauixi, cariapé e/ou areias, ricas em sílica, foram intencionalmente adicionados. O elevado conteúdo de fosfatos de Al-Fe, amorfos ou como de baixa cristalinidade, originou-se a partir do contato entre a matriz argilosa da parede do vaso cerâmico com a solução aquosa quente durante o cozimento diário de alimentos de origem animal (principal fonte de fósforo). A cristalização dos fosfatos deve ter prosseguida mesmo depois que os vasos foram descartados, e juntos com os restos de matéria orgânica vegetal e animal incorporaram-se aos solos residuais. Participaram desta forma na formação dos solos tipo TPA.

Eletroresistividade na prospecção hidrogeológica no município de Ourilândia do Norte-PA

A prospecção geofísica para a exploração de água subterrânea em regiões cristalinas torna-se uma tarefa extremamente difícil pela presença de estruturas como falhas, fraturas e ocorrência de variação lateral de resistividade no manto intempérico. Para analisar a influência da variação bidimensional da propriedade física no comportamento das curvas de sondagem elétrica vertical Schlumberger, simulou-se modelos contendo tais estruturas encaixadas no manto de alteração com ângulo de mergulho de 30° e 90° e com a posição do centro de sondagem variando em relação a heterogeneidade. A modelagem numérica foi feita através do método dos elementos finitos com o programa EGSLIB/SEV2D. O resultado da simulação mostra que o emprego da técnica da sondagem elétrica vertical Schlumberger sobre meios bidimensionais não funciona eficientemente, como quando aplicada em regiões sedimentares e de aluvião, devido o fenômeno da ambigüidade que dificulta a interpretação das curvas. Dentre os métodos geofísicos para a investigação hidrogeológica, a aplicação dos métodos eletrorresistivos mostra-se mais simples, eficaz e econômica. Assim, o levantamento geofísico para orientar os trabalhos de perfurações de poços na cidade de Ourilândia do Norte contou de trinta e duas sondagens elétricas verticais tipo Schlumberger e um caminhamento elétrico dipolo-dipolo realizados em diversas ruas da cidade. As curvas de resistividade aparente, obtidas com as sondagens foram tratadas e processadas considerando um modelo de camadas planas, horizontais, homogêneas e isotropicas, inicialmente com o algoritmo EGSLIB/SEV1D e posteriormente com uso do programa EGSLIB/SEV1DINV, permitindo estabelecer modelos geoelétricos aproximados para a configuração da subsuperfície da cidade. Depois da interpretação quantitativa, foi possível separar duas zonas geoelétricas distintas. Esta diferenciação foi realizada com base nos valores de resistividade aparente e as informações geológicas de superfície. A primeira, com pouca representatividade na área da cidade, apresenta-se constituída por solo/aterro sobre saprolito de natureza argilo-arenosa recobrindo os granitóides tipo Rio Maria. A segunda composta pelo solo/aterro seguido de um horizonte preenchido por sedimentos arenosos superposto a camada argilo-arenosa e finalmente a rocha sã. Esta seqüência geoelétrica predomina na cidade. Os estratos arenoso e argilo-arenoso apresentam condições potenciais para armazenamento de água.

Um estudo da estabilidade mecânica de poços a partir de perfis geofísicos em formações pouco consolidadas

A análise da estabilidade mecânica de um poço pode ser feita a partir do cálculo de parâmetros elásticos da formação utilizando a densidade do meio e as velocidades de propagação das ondas compressional e cisalhante na formação rochosa, os quais podem ser obtidos de perfis geofísicos do poço. Em formações sedimentares pouco consolidadas as ferramentas de perfilagem sônica convencionais (monopolares) não conseguem registrar com acuidade a velocidade da onda cisalhante pois a primeira chegada dessa onda é camuflada pela chegada de outras ondas que podem ser mais rápidas que a onda cisalhante num poço perfurado neste tipo de formação. Medidas das velocidades sônicas são feitas em laboratório em amostras da formação, sob condições semelhantes às condições in situ, servindo como ajuste das velocidades registradas no poço pela ferramenta de perfilagem sônica. Para a análise de estabilidade da formação, perfis auxiliares são necessários como o perfil de porosidade, saturação de fluidos e perfis de composição mineralógica da formação rochosa. Exige-se ainda dados de testes de avaliação da formação e de condições do reservatório, mas que são comuns em poços de petróleo, como o teste de formação e os testes de pressurização do poço, tais como o teste de micro-fraturamento hidráulico ou o teste de absorção. A avaliação das tensões principais efetivas que atuam distante do poço, e que não são afetadas pela sua presença, é feita através da associação de um modelo de deformação elástica apropriado e o resultado do teste de pressurização disponível para o poço em estudo. Utilizando resultados clássicos da teoria da elasticidade geral pode-se calcular o campo de tensões modificado na vizinhança da parede do poço devido ao efeito da própria presença do poço ali perfurado e da diferença de pressão existente entre o interior do poço e a formação rochosa. A determinação das propriedades mecânicas da formação a partir das velocidades sônicas e a avaliação do campo de tensões assumindo um modelo elástico de deformação, supõem o meio rochoso no qual as ondas se propagam como elástico, homogêneo e isotrópico. Esta suposição representa a principal aproximação assumida pela metodologia descrita neste trabalho. De posse das propriedades mecânicas da formação e do campo de tensões que age na vizinhança do poço resta definir o critério segundo o qual a rocha sofre instabilidade mecânica quando submetida aquele campo de tensões. Isto permite determinar se, nas condições avaliadas do poço e da formação, haverá quebra da parede do poço por excesso de tensão e, se houver, qual a sua extensão. Assim o problema é como analisar o comportamento mecânico de um poço em uma formação pouco consolidada a partir de perfis geofísicos os quais podem ter problemas no registro das propriedades físicas do meio em formações deste tipo. A metodologia proposta é aplicada a dois intervalos de profundidade pertencentes a dois poços onde arenitos e folhelhos se intercalam e nos quais todos os dados necessários estão disponíveis. Os resultados obtidos mostram que, exceto quando outros mecanismos de quebra da parede do poço agem na formação, a metodologia proposta consegue com sucesso detectar zonas de ocorrência de instabilidade mecânica do poço provocadas por um campo de tensões que excede a resistência mecânica da formação.

Fácies e proveniência de depósitos costeiros da Formação Raizama: evidências do registro Ediacarano-cambriano na faixa Paraguai, região de Nobres, Mato Grosso

Rochas siliciclásticas da Formação Raizama, unidade basal do Grupo Alto Paraguai de idade ediacarana-cambriana (635 – 541 Ma), ocorrem distribuídas descontinuamente ao longo da margem sul do Cráton Amazônico e segmento norte da Faixa Paraguai, centro-oeste do Brasil, estado do Mato Grosso. Estas recobrem discordantemente os depósitos de plataforma carbonática do Grupo Araras, onde foram registrados evidências do evento glacial Marinoano (635 Ma). O Grupo Alto Paraguai representa os estágios finais da colisão entre os blocos Amazônico e Paranapanema que culminaram no fechamento do Oceano Clymene (540-520 Ma). A Formação Raizama com espessura de aproximadamente 570 m é constituída por pelitos, arenitos finos a grossos, e arenitos com cimento dolomítico previamente interpretados como depósitos flúvio-costeiros distribuídos nos membros inferior (270 m) e superior (300 m). O estudo faciológico e estratigráfico desta unidade na região de Nobres, Estado do Mato Grosso, foi focado principalmente na seção aflorante de 600 m no leito do rio Serragem, que inclui a Cachoeira da Serra do Tombador. Foram definidas 17 fácies sedimentares, agrupadas em cinco associações de fácies (AF) representativas de uma sucessão costeira progradante iniciando por depósitos de shoreface inferior, os quais recobrem em conformidade correlativa os depósitos de plataforma carbonática da Formação Serra do Quilombo (Grupo Araras). A AF1 consiste em arenitos com laminação plano-paralela e laminação truncada por onda (microhummocky), individualizada por camadas de pelito laminado interpretadas como depósitos de shoreface inferior. Destaca-se na AF1 a primeira ocorrência de níveis centimétricos bioturbados por Skolithos em depósitos neoproterozoicos – cambrianos na Faixa Paraguai. A AF2 é formada por arenitos com estratificação cruzada swaley e estratificação plano-paralela interpretada como depósitos de shoreface superior. A AF3 é composta por arenitos com estratificações cruzadas tangenciais e acanaladas com recobrimentos de siltito/arenito muito fino representativos de depósitos de canal e barras de submaré. A AF4 é caracterizada por arenitos com estratificações cruzadas tangencial e sigmoidal, laminação plano-paralela a cruzadas de baixo-ângulo, ritmito arenito muito fino/siltito com acamamento flaser e gretas de contração, organizados em ciclos métricos de raseamento ascendente de planície de maré. A AF5 é constituída por arenito com estratificação cruzada acanalada marcada por lags residuais na base da associação, arenito com estratificações plano-paralela e cruzada de baixo-ângulo, interpretados como depósitos fluviais distais de rios entrelaçados, parcialmente retrabalhados por ondas. Grãos detríticos de zircão foram obtidos da AF3 e datados pelo método U-Pb, sendo a idade de 1001±9 Ma interpretada como a idade de máxima deposição da Formação Raizama. Aliado a tal análise, as paleocorrentes NE-SE mostram que estes grãos teriam como áreas fontes principais a Faixa Sunsás, SW do Cráton Amazônico, não sendo descartada contribuições oriundas da parte NW desse Cráton. A idade mesoproterozóica obtida serviu principalmente para desvendar a proveniência da Formação Raizama, enquanto que as datações da base do Grupo Araras, em torno de 627-622 Ma, associada à presença inequívoca do icnogênero Skolithos, tornam esta unidade muito mais próxima do limite com o Cambriano Inferior. Traços fósseis do Proterozoico são caracterizados quase que exclusivamente por traços horizontais, sendo que bioturbações verticais praticamente são ausentes ao longo do Neoproterozoico. Esta inferência vem de encontro com a idade máxima de 541 Ma obtida para a Formação Diamantino, a qual recobre a unidade estudada. Os dados radiométricos aliados com as interpretações paleoambientais, que incluem o registro das primeiras atividades de organismos perfurantes na Faixa Paraguai, abrem perspectivas de entender com maiores detalhes a sequência de eventos que tipificam os estratos do limite Ediacarano-Cambriano do Brasil, ainda pouco conhecidos.

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.

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

Investigação por eletroresistividade de ambientes cársticos no município de São João de Pirabas - PA

A presente tese procura investigar ambientes cársticos, usando os métodos eletroresistivos de sondagem elétrica vertical e caminhamento dipolo-dipolo. A área piloto para essa investigação está localizada na porção nordeste do Estado do Pará, no Município de São João de Pirabas. Nessa localidade, durante a perfuração de um poço para abastecimento de água para a cidade, a equipe de prospecção da Fundação Nacional de Saúde encontrou problemas devido a presença de cavidades nos calcários da Formação Pirabas. Para analisar esse problema foram realizadas 15 sondagens elétricas verticais, usando o arranjo Schlumberger, em diversos pontos da cidade. A abertura máxima dos eletrodos de corrente, para a realização das sondagens, foi em torno de 600 metros. As curvas de resistividade aparente das sondagens foram, inicialmente, interpretadas usando programas computacionais que calculam as resistividade do meio investigado com modelos unidimensionais, sem variações laterais de resistividade. Porém, as sondagens realizadas próximas ao poço, onde ocorrem cavernas, foram interpretadas usando modelos bidimensionais para representar as variações laterais de resistividade devido às cavernas. O algoritmo usado para essa interpretação se baseia na técnica dos elementos finitos. Os resultados das interpretações das sondagens elétricas, para os vários modelos de cavernas estudados, mostraram que as anomalias de resistividade devido a presença de cavernas são muito pequenas. Devido a este fato, o método de sondagem elétrica vertical não apresentou boa resolução para identificar com precisão a presença deste tipo de caverna em subsuperfície. Por outro lado, os resultados numéricos obtidos pelo caminhamento dipolo-dipolo mostraram que esse método é mais eficiente e menos ambíguo, para estudar cavernas, que o método de sondagem elétrica vertical. Contudo o método de caminhamento dipolo-dipolo é muito mais dispendioso do que o de sondagem elétrica verticais, tanto do ponto de vista de trabalho de campo como no processamento e interpretação de dados.