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.