Amplitudes e padrões de polarização de pulsos em meios anisotrópicos

Extrair informações litológicas da subsuperfície através de dados sísmicos constitui-se num grande desafio à prospecção sísmica, pois a hipótese de estratificações formadas por camadas isotrópicas se mostra insuficiente para representar o comportamento do campo elástico em levantamentos com grandes afastamentos entre fonte e receptor, geofones multicomponentes, medidas de VSP tridimensional, entre outros. Sob este panorama, a prospecção sísmica passa a considerar modelos anisotrópicos de subsuperfície para, por exemplo, caracterizar reservatórios. O objetivo deste texto é apresentar um formalismo para modelar o espalhamento de pulsos a partir de ondas planas incidentes em interfaces planas horizontais que separam meios anisotrópicos. Este espalhamento é obtido primeiramente, através da formulação explícita dos campos de deformação e tração como função das matrizes propagadoras, de polarização e de impedância do meio. Em seguidaeste formalismo é usado para a obtenção das matrizes dos coeficientes de reflexão e transmissão através de uma interface plana horizontal para posteriormente, ser generalizado para o espalhamento através de múltiplas camadas. Finalmente, inserem-se ao campo da onda incidente as amplitudes de um pulso analítico para calcular o espalhamento do pulso através de estratificações.

Deconvolução de Euler: passado, presente e futuro – um tutorial

Neste tutorial apresentamos uma revisão da deconvolução de Euler que consiste de três partes. Na primeira parte, recordamos o papel da clássica formulação da deconvolução de Euler 2D e 3D como um método para localizar automaticamente fontes de campos potenciais anômalas e apontamos as dificuldades desta formulação: a presença de uma indesejável nuvem de soluções, o critério empírico usado para determinar o índice estrutural (um parâmetro relacionado com a natureza da fonte anômala), a exeqüibilidade da aplicação da deconvolução de Euler a levantamentos magnéticos terrestres, e a determinação do mergulho e do contraste de susceptibilidade magnética de contatos geológicos (ou o produto do contraste de susceptibilidade e a espessura quando aplicado a dique fino). Na segunda parte, apresentamos as recentes melhorias objetivando minimizar algumas dificuldades apresentadas na primeira parte deste tutorial. Entre estas melhorias incluem-se: i) a seleção das soluções essencialmente associadas com observações apresentando alta razão sinal-ruído; ii) o uso da correlação entre a estimativa do nível de base da anomalia e a própria anomalia observada ou a combinação da deconvolução de Euler com o sinal analítico para determinação do índice estrutural; iii) a combinação dos resultados de (i) e (ii), permitindo estimar o índice estrutural independentemente do número de soluções; desta forma, um menor número de observações (tal como em levantamentos terrestres) pode ser usado; iv) a introdução de equações adicionais independentes da equação de Euler que permitem estimar o mergulho e o contraste de susceptibilidade das fontes magnéticas 2D. Na terceira parte apresentaremos um prognóstico sobre futuros desenvolvimentos a curto e médio prazo envolvendo a deconvolução de Euler. As principais perspectivas são: i) novos ataques aos problemas selecionados na segunda parte deste tutorial; ii) desenvolvimento de métodos que permitam considerar interferências de fontes localizadas ao lado ou acima da fonte principal, e iii) uso das estimativas de localização da fonte anômala produzidas pela deconvolução de Euler como vínculos em métodos de inversão para obter a delineação das fontes em um ambiente computacional amigável.

Composição, riqueza e abundância de espécies de anfíbios na região do médio rio Xingú

Este estudo teve como objetivo inventariar a fauna de anfíbios nos diferentes hábitats da região do Médio Rio Xingu. Para tanto, foram selecionadas seis localidades, sendo as amostragens realizadas em três incursões, entre novembro de 2007 e março de 2008, totalizando 48 dias. Foram identificados seis tipos de hábitats na região amostrada: floresta de terra-firme, floresta sazonalmente alagável (várzea), margem do rio, “pedral”, lagoa e área alterada. Para a coleta de dados foram utilizadas duas metodologias: armadilhas de interceptação e queda e procura ativa. As armadilhas foram instaladas apenas nas áreas de floresta de terra-firme, enquanto a procura ativa foi empregada em todos os tipos de hábitats identificados. As amostragens resultaram no registro de 56 espécies de anfíbios e outras oito espécies foram identificadas e registradas em estudos anteriores na área. A diversidade e riqueza de espécies foram maiores na localidade Caracol, onde predominou a floresta de terra-firme, e menor na localidade Ilha Grande, dominada por floresta sazonalmente alagável. Foram registrados nove modos reprodutivos na área de estudo, sendo todos observados na floresta de terra-firme e apenas três nos “pedrais”, o que parece refletir a baixa heterogeneidade ambiental dessa área. Para comparar as localidades estudadas foram realizadas duas análises de similaridade, uma para cada metodologia de coleta utilizada. A análise de similaridade dos dados de procura ativa apontou maior semelhança na composição das espécies entre as áreas de floresta de terra-firme do que entre os outros hábitats. A análise de agrupamento realizada entre a composição de espécies obtida nesse estudo e outros levantamentos realizados na Amazônia agrupou esta área com outra também localizada no Médio Xingu.

Identificação molecular de espécies do complexo Anopheles albitarsis (Diptera, Culicidae, Anophelinae), coletadas em dois municípios da Amazônia brasileira, com análise de suscetibilidade natural a plasmódios humanos

Anofelinos membros de complexos de espécies crípticas podem exibir diferenças comportamentais, de susceptibilidade a infecção malárica, e resistência a inseticidas. Assim, a identificação de espécies vetoras tem relevância epidemiológica, o que nem sempre é possível por critérios morfológicos. Métodos alternativos têm sido empregados para tal, como os que analisam regiões altamente conservadas do DNA ribossômico, variável entre as espécies, conhecidas como espaçadoras internas transcritas (ITS). Considera-se atualmente que o complexo Anopheles c seja composto por seis espécies: An. albitarsis s.s., An. oryzalimnetes, An. albitarsis F, An. marajoara, An. deaneorum, e An. janconnae. Destas, pelo menos as três últimas são incriminadas como vetores de malária na Amazônia brasileira. O objetivo deste estudo foi realizar identificação molecular de espécies do complexo An. albitarsis, por análise da seqüências do ITS2 do rDNA, com vistas a analisar sua importância na transmissão de malária nos municípios de Macapá, Amapá e Peixe-Boi, Pará, inclusive investigando pela primeira vez a ocorrência do An. albitarsis F nestas duas áreas epidemiologicamente distintas: a primeira com histórico de alto risco de transmissão de malária e a segunda não. O estudo foi realizado entre janeiro de 2009 e abril de 2010, e consistiu de capturas de anofelinos de 12 horas de duração (ecostofase) no peridomicílio. Todas as fêmeas coletadas foram morfologicamente identificadas e apenas os An. albitarsis s.l. tiveram cabeça e tórax separadas para análise da infecção natural por ELISA; ovários para análise de paridade e patas, asas e carcaça para identificação molecular. Em Macapá foram realizadas seis coletas, obtendo-se um total de 584 anofelinos, sendo 366 An. albitarsis s.l. (62,7%), 167 An. darlingi (28,6%), 33 An. triannulatus s.l (5,6%), 15 An. braziliensis (2,6%) e 3 An. nuneztovari (0,5%). Pela PCR foi possível visualizar a banda específica de An. marajoara em 320 espécimes dos An. albitarsis s.l testados. Do restante, 33 foram negativos e 13 amplificaram um fragmento de ~490 pb nos iniciadores empregados, não permitindo chegar ao diagnóstico específico. O An. marajoara apresentou características biológicas e comportamentais que ratificam sua importância epidemiológica na transmissão de malária em Macapá, tais como: ser a espécie mais prevalente, com maior proporção de fêmeas paridas (73,0%), e portanto com maiores chances de se infectarem com o plasmódio, ocorrer tanto na estação menos quanto na mais chuvosa, e apresentar atividade hematofágica durante toda a ecostofase, alem disso, foi encontrado naturalmente infectado por P. vivax e P. falciparum (taxa de infecção natural de 3,1%). Em Peixe-Boi, foram capturados 43 anofelinos: An. triannulatus s.l (20 espécimes, 46,5 %), An. albitarsis s.l. (13: 30,2 %), An. darlingi (8: 18,6%), e An. nuneztovari (2: 4,7%). Todos os An. albitarsis s.l. coletados foram identificados pela ITS2 como An. oryzalimnetes. Nenhum deles foi encontrado infectado pelos plasmódios testados, e a maioria das fêmeas era parida (84,6%). São necessários levantamentos entomológicos sistemáticos que analisem a importância deste anofelino na transmissão de malária na cidade. O An. albitarsis F não foi encontrado nas duas áreas estudadas. Nossos resultados contribuem para o entendimento da epidemiologia da malária na região Amazônica brasileira.

Ordem e desordem do território turístico: a chegada do estranho e os conflitos de territorialidades na orla oeste de Mosqueiro, Belém/PA

A Ilha de Mosqueiro, localizada ao norte da sede municipal de Belém, consiste no principal espaço turístico para a população belenense e num tradicional espaço do turismo de segunda-residência, principalmente as praias de sua orla oeste. Diante da realidade construída neste espaço a partir do aumento da acessibilidade e da conseqüente intensificação e diversificação das práticas turísticas, o objetivo geral da pesquisa consiste em analisar o conflito entre as territorialidades constituídas pelas práticas turísticas de segunda residência e excursionista nas praias da orla oeste de Mosqueiro e seus reflexos no reordenamento territorial dessa localidade. Para viabilizar o alcance desse objetivo foram realizados levantamentos bibliográficos e documentais, entrevistas com gestores e especialistas ligados à questão e entrevistas com turistas de segunda residência e excursionistas. Os resultados da pesquisa demonstram que a dinâmica de ordem/desordem do território turístico da orla oeste de Mosqueiro é fruto do conflito de territorialidades entre as práticas turísticas com vistas ao maior controle do espaço e à implementação de um ordenamento territorial mais adequado à natureza de sua respectiva prática. Os discursos adquirem grande importância na manifestação desse conflito, especialmente no que se refere aos turistas de segunda residência, que evidenciam forte intencionalidade nas ações. Os principais reordenamentos territoriais inerentes a esse conflito se materializam por meio de três planos principais: a territorialização da prática turística excursionista na orla oeste de Mosqueiro; a desterritorialização da prática turística de segunda residência; e o rearranjo territorial da prática turística de segunda residência.

Rede turística e organização espacial: uma análise da Ilha de Mosqueiro, Belém/PA

A Ilha do Mosqueiro, distrito administrativo do município de Belém localizado ao norte da sede deste, consiste num destacado espaço turístico, cuja tradicional importância regional ligada ao turismo de praia é acrescentada a emergência de práticas turísticas delineadas a partir da exuberância natural tipicamente amazônica, como o turismo de natureza e o ecoturismo; e pelo aumento da importância de fluxos nacionais e internacionais. Tomando por área de estudo, mais precisamente, a faixa de orla da Ilha do Mosqueiro que se estende do extremo oeste ao extremo norte dessa localidade, objetivou-se analisar e compreender a configuração da rede turística da Ilha de Mosqueiro, seus condicionantes sociais e políticos, níveis de solidariedade transescalar e mobilidade sócio-espacial, e como isso pode contribuir ao entendimento da organização espacial do turismo no momento contemporâneo.Os procedimentos metodológicos utilizados com vista ao alcance do objetivo exposto acima consistiram em levantamentos bibliográficos, levantamentos documentais, trabalho de campo de caráter exploratório a partir da técnica de observação dirigida; entrevistas semi-estruturadas com representante das Secretarias Municipais de Belém, da Secretaria Estadual de Turismo, da Agência Distrital de Mosqueiro, da Polícia Militar paraense e de Associações locais; entrevistas estruturadas com turistas, agentes e operadores turísticos, proprietários de restaurantes, bares e outros comerciantes locais. A rede turística, que pode ser definida a partir da configuração das redes de suporte, tem seu estabelecimento a partir da articulação entre áreas emissoras e receptoras de fluxos turísticos, nas quais de configuram nódulos privilegiados conectados por meios de alimentação, meios de hospedagem, meios de entretenimento, empresas, atrativos e demanda turísticos nas mais variadas escalas geográficas (COSTA, RIBEIRO e TAVARES, 2004; COSTA e TAVARES, 2005). Em Mosqueiro a estruturação espacial das redes técnicas deu suporte a configuração espacial da rede turística, a qual abarca três níveis de solidariedade: local, nacional e internacional; que delineiam as escalas das ações abarcadas pela prática turística desenvolvida nessa localidade. O resultado dessa dinâmica da rede turística pode ser percebido através do estabelecimento de uma organização espacial estruturada a partir de três nodosidades: orlas das praias do Murubira, Chapéu Virado e Farol; orla da praia do Paraíso e orla de várzea das trilhas ecológicas; as quais, respectivamente, estão relacionadas a eventos turísticos de sol e praia, de turismo na natureza e de turismo ecológico de excursão, os quais são efetuados por movimentos de segunda residência, de hotelaria e de excursionismo.

Regionalização e gestão no espaço turístico: o processo de roteirização e de gestão participativa do Pólo Marajó, Pará

O presente trabalho analisa no contexto da regionalização e gestão no espaço turístico, dois processos inerentes ao Programa de Regionalização do Turismo (PRT), a roteirização e a gestão participativa, tendo como lócus de pesquisa o pólo Marajó. O texto procura evidenciar os fundamentos teóricos do espaço turístico, a partir do pensamento de BOULLÓN (1985), associando esta abordagem com as formulações contidas na teoria dos pólos de crescimento (PERROUX, 1967), no aménagement du territoire (ANDRADE, 1971, 1987), na teoria da inovação (PORTER, 1998) e na desenvolvimento local (Almeida, 2002, Barqueiro, 2001, Buarque, 2002, Rodrigues, 1977). Analisa-se as políticas públicas nacionais recentes, o Programa Nacional de Municipalização do Turismo (PNMT), o PRT, além do Plano de Desenvolvimento Turístico do Estado do Pará, identificando a orientação desses documentos técnicos quanto ao processo de desenvolvimento do turismo. Para que fossem alcançados os objetivos propostos foram formuladas questões orientadoras da análise do objeto de estudo, quanto a Roteirização e Gestão Participativa. Na seqüência, cumpriu-se uma etapa de pesquisa, de caráter bibliográfico e documental, que permitiu a construção do referencial teórico. Complementarmente, realizou-se o trabalho de campo com entrevistas, observações, levantamentos de informações na área em estudo. O resultado desse processo foi sistematizado, a partir da interpretação das informações obtidas no trabalho de gabinete e de campo. Consolidadas estas informações, as mesmas foram organizadas em forma de texto analítico, com o suporte de mapas, fotos, quadro e tabelas. Considerando os dois principais objetivos propostos (análise da roteirização e da gestão participativa no pólo Marajó), observou-se que o trabalho de roteirização na vila do Pesqueiro concentrou-se, basicamente, em dois dos quatro componentes do espaço turístico: a matéria prima e a superestrutura. Os demais componentes (planta turística e infra-estrutura) não receberam atenção adequada. Tais fatos, conjugados a fatores externos ao universo composto pela comunidade local, como a questão do transporte para o Marajó, e a atuação dos agentes de receptivo de Belém - que têm grande influência na comercialização dos roteiros - dificultaram o processo de venda do produto. Quanto a Gestão Participativa , observou-se que foram instituídos dois comitês gestores, um em Soure e outro em Salvaterra, dissociados do processo de constituição do Fórum Regional de Turismo do Pólo Marajó. Além disso, a estrutura do Fórum não obedeceu ao princípio da equidade, em termos de representação dos atores, e a concepção da governança regional do turismo ocorreu após o processo de roteirização do Amazônia do Marajó. Esses elementos indicam que ainda não se pratica, concretamente, uma gestão participativa no pólo Marajó.

Interpolação de dados de campo potencial através da camada equivalente

O uso da técnica da camada equivalente na interpolação de dados de campo potencial permite levar em consideração que a anomalia, gravimétrica ou magnética, a ser interpolada é uma função harmônica. Entretanto, esta técnica tem aplicação computacional restrita aos levantamentos com pequeno número de dados, uma vez que ela exige a solução de um problema de mínimos quadrados com ordem igual a este número. Para viabilizar a aplicação da técnica da camada equivalente aos levantamentos com grande número de dados, nós desenvolvemos o conceito de observações equivalentes e o método EGTG, que, respectivamente, diminui a demanda em memória do computador e otimiza as avaliações dos produtos internos inerentes à solução dos problemas de mínimos quadrados. Basicamente, o conceito de observações equivalentes consiste em selecionar algumas observações, entre todas as observações originais, tais que o ajuste por mínimos quadrados, que ajusta as observações selecionadas, ajusta automaticamente (dentro de um critério de tolerância pré-estabelecido) todas as demais que não foram escolhidas. As observações selecionadas são denominadas observações equivalentes e as restantes são denominadas observações redundantes. Isto corresponde a partir o sistema linear original em dois sistemas lineares com ordens menores. O primeiro com apenas as observações equivalentes e o segundo apenas com as observações redundantes, de tal forma que a solução de mínimos quadrados, obtida a partir do primeiro sistema linear, é também a solução do segundo sistema. Este procedimento possibilita ajustar todos os dados amostrados usando apenas as observações equivalentes (e não todas as observações originais) o que reduz a quantidade de operações e a utilização de memória pelo computador. O método EGTG consiste, primeiramente, em identificar o produto interno como sendo uma integração discreta de uma integral analítica conhecida e, em seguida, em substituir a integração discreta pela avaliação do resultado da integral analítica. Este método deve ser aplicado quando a avaliação da integral analítica exigir menor quantidade de cálculos do que a exigida para computar a avaliação da integral discreta. Para determinar as observações equivalentes, nós desenvolvemos dois algoritmos iterativos denominados DOE e DOEg. O primeiro algoritmo identifica as observações equivalentes do sistema linear como um todo, enquanto que o segundo as identifica em subsistemas disjuntos do sistema linear original. Cada iteração do algoritmo DOEg consiste de uma aplicação do algoritmo DOE em uma partição do sistema linear original. Na interpolação, o algoritmo DOE fornece uma superfície interpoladora que ajusta todos os dados permitindo a interpolação na forma global. O algoritmo DOEg, por outro lado, otimiza a interpolação na forma local uma vez que ele emprega somente as observações equivalentes, em contraste com os algoritmos existentes para a interpolação local que empregam todas as observações. Os métodos de interpolação utilizando a técnica da camada equivalente e o método da mínima curvatura foram comparados quanto às suas capacidades de recuperar os valores verdadeiros da anomalia durante o processo de interpolação. Os testes utilizaram dados sintéticos (produzidos por modelos de fontes prismáticas) a partir dos quais os valores interpolados sobre a malha regular foram obtidos. Estes valores interpolados foram comparados com os valores teóricos, calculados a partir do modelo de fontes sobre a mesma malha, permitindo avaliar a eficiência do método de interpolação em recuperar os verdadeiros valores da anomalia. Em todos os testes realizados o método da camada equivalente recuperou mais fielmente o valor verdadeiro da anomalia do que o método da mínima curvatura. Particularmente em situações de sub-amostragem, o método da mínima curvatura se mostrou incapaz de recuperar o valor verdadeiro da anomalia nos lugares em que ela apresentou curvaturas mais pronunciadas. Para dados adquiridos em níveis diferentes o método da mínima curvatura apresentou o seu pior desempenho, ao contrário do método da camada equivalente que realizou, simultaneamente, a interpolação e o nivelamento. Utilizando o algoritmo DOE foi possível aplicar a técnica da camada equivalente na interpolação (na forma global) dos 3137 dados de anomalia ar-livre de parte do levantamento marinho Equant-2 e 4941 dados de anomalia magnética de campo total de parte do levantamento aeromagnético Carauari-Norte. Os números de observações equivalentes identificados em cada caso foram, respectivamente, iguais a 294 e 299. Utilizando o algoritmo DOEg nós otimizamos a interpolação (na forma local) da totalidade dos dados de ambos os levantamentos citados. Todas as interpolações realizadas não seriam possíveis sem a aplicação do conceito de observações equivalentes. A proporção entre o tempo de CPU (rodando os programas no mesmo espaço de memória) gasto pelo método da mínima curvatura e pela camada equivalente (interpolação global) foi de 1:31. Esta razão para a interpolação local foi praticamente de 1:1.

Inversão vinculada de dados de EMAP e MT-1D anisotrópico

Uma técnica para a inversão de dados magnetotelúricos é apresentada neste trabalho. Dois tipos de dados são tratados aqui, dados gerados por modelos unidimensionais com anisotropia na condutividade das camadas e dados bi-dimensionais de levantamentos do método EMAP (ElectroMagnetic Array Profiling). Em ambos os casos fazemos a inversão usando vínculos aproximados de igualdade para estabilizar as soluções. Mostramos as vantagens e as limitações do uso destes vínculos nos processos de inversão. Mesmo vinculada a inversão ainda pode se tornar instável. Para inverter os dados 2-D do EMAP, apresentamos um processo que consiste de três partes: 1 – A construção de um modelo interpretativo e da aproximação inicial para a inversão a partir dos dados de seções de resistividade aparente filtradas pelo processo de filtragem do EMAP; 2 – a inclusão de uma camada de corpos pequenos aflorantes, chamada de camada destatic shift, aos modelos interpretativos para resolver as fontes de distorções estáticas que contaminam os dados; 3 – o uso dos vínculos aproximados de igualdade absoluta para estabilizar as soluções. Os dois primeiros passos nos permitem extrair o máximo de informação possível dos dados, enquanto que o uso dos vínculos de igualdade nos permite incluir informação a priori que possua significado físico e geológico. Com estes passos, obtemos uma solução estável e significativa. Estudaremos o método em dados sintéticos de modelos bi-dimensionais e em dados reais de uma linha de EMAP feita na Bacia do Paraná.

Modelagem sísmica e inversão na presença de anisotropia

A necessidade da adoção de modelos elásticos anisotrópicos, no contexto da sísmica de exploração, vem crescendo com o advento de novas técnicas de aquisição de dados como VSP, walkway VSP, tomografia poço a poço e levantamentos sísmicos com grande afastamento. Meios anisotrópicos, no contexto da sísmica de exploração, são modelos efetivos para explicar a propagação de ondas através de meios que apresentam padrões de heterogeneidade em escala muito menor que o comprimento de onda das ondas sísmicas. Particularmente, estes modelos são muito úteis para explicar o dado sísmico mais robusto que são as medidas de tempo de trânsito. Neste trabalho, são investigados aspectos da propagação de ondas, traçado de raios e inversão de tempos de trânsito em meios anisotrópicos. É estudada a propagação de ondas SH em meios anisotrópicos estratificados na situação mais geral onde estas ondas podem ocorrer, ou seja, em meios monoclínicos com um plano vertical de simetria especular. É mostrado que o campo de ondas SH refletido a partir de um semi-espaço estratificado, não apresenta qualquer informação sobre a possível presença de anisotropia em subsuperfície. São apresentados métodos simples e eficientes para o traçado de raios em 3D através de meios anisotrópicos estratificados, baseados no princípio de Fermat. Estes métodos constituem o primeiro passo para o desenvolvimento de algoritmos de inversão de tempos de trânsito para meios anisotrópicos em 3D, a partir de dados de VSP e walkaway VSP. Esta abordagem é promissora para determinação de modelos de velocidade, que são necessários para migração de dados sísmicos 3D na presença de anisotropia. É efetuada a análise da inversão tomográfica não linear, para meios estratificados transversalmente isotrópicos com um eixo de simetria vertical(TIV). As limitações dos dados de tempo de trânsito de eventos qP para determinação das constantes elásticas, são estabelecidas e caracterizados os efeitos da falta de cobertura angular completa na inversão tomográfica. Um algoritmo de inversão foi desenvolvido e avaliado em dados sintéticos. A aplicação do algoritmo a dados reais demonstra a consistência de meios TIV. Esta abordagem é útil para casos onde há informação a priori sobre a estratificação quase plana das formações e onde os próprios dados do levantamento poço a poço apresentam um alto grau de simetria especular em relação a um plano vertical. Também pode ser útil em interpretações preliminares, onde a estimativa de um meio estratificado, serve como modelo de fundo para se efetuar análises mais detalhadas, por exemplo, como um modelo de velocidades anisotrópico para migração, ou como um modelo de calibração para análises de AVO.

Análise de resolução em tomografia interpoços em 3-D

A aplicação de métodos de inversão aos dados de múltiplos levantamentos sísmicos interpoços permite a reconstrução de modelos de vagarosidade em 3-D de alta resolução adequados para monitoramento de processos de recuperação avançada de petróleo e caracterização de reservatórios. Entretanto, a falta de cobertura volumétrica uniforme de raios de levantamentos interpoços exige informação adicional ao sistema tomográfico para obtenção de soluções estáveis. A discretização do modelo em uma malha 3-D com células prismáticas triangulares e a decomposição em valores singulares são utilizadas para avaliar a reconstrução tomográfica em 3-D. O ângulo da projeção de modelos-alvo no subespaço ortogonal ao espaço nulo efetivo da matriz tomográfica é um critério adequado para se otimizar a malha de discretização do modelo interpretativo e a geometria de aquisição dos dados de modo a melhorar o condicionamento da reconstrução tomográfica. Esta abordagem pode ser utilizada durante as iterações lineares para redefinir a malha ou avaliar a necessidade de informação a priori adicional ao sistema tomográfico.

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.

Aplicação dos métodos gravimétrico e magnetométrico para a definição do comportamento estrutural da faixa de dobramentos Araguaia

O presente trabalho consiste na interpretação de informações gravimétricas e aeromagnetométricas (Projeto Geofísico Brasil-Canadá - PGBC) da região setentrional da Faixa de Dobramentos Araguaia, envolvendo uma área de aproximadamente 129.000km², compreendida entre os meridianos 47°50’W e 50°30’W e paralelos 4°50'S e 9°00'S. Abrange porções sudeste do Estado do Pará, noroeste de Goiás e oeste do Maranhão. Os trabalhos de campo constaram de levantamentos gravimétrico e altimétrico, ao longo de rodovias que constituem a rede viária regional, e os resultados obtidos, após correções e reduções, foram então organizados de forma a constituírem um mapa de anomalias Bouguer, apresentado em escala 1:500.000. O padrão gravimétrico da faixa de dobramentos Araguaia se caracteriza por mostrar feições predominantemente longitudinais, com curvas isoanômalas de direções submeridianas, concordante com o comportamento litológico-estrutural conhecido para a área. Na porção centro-meridional desse mapa aparece uma zona de anomalia negativa, alongada e intensa (menos de -100mgal), a qual sofre duas importantes inflexões para NW, sendo uma correlacionável ao Lineamento Carajás, e a outra na altura do paralelo 5°30'S. Apresenta zonas de altos gravimétricos, notadamente na porção noroeste (domínio cratônico), relacionada à influência das metavulcânicas do Grupo Grão Pará, e na porção nordeste, devida a massas densas introduzidas na crosta. Destaca-se também o relativo alto gravimétrico acompanhando o flanco oriental do eixo de ocorrência das braquidobras, podendo ser devido a massas excedentes colocadas sob essa região, através de esforços tectônicos, produzindo elevações do nível de base da crosta e consequentemente do embasamento. De forma genérica o flanco oriental da faixa mostra um gradiente mais intenso que o ocidental, devido principalmente à configuração geométrica das estruturas nesse setor. Da mesma forma merece destaque o gradiente regional ascendente de sul para norte. A análise das cartas de intensidade magnética do PGBC fornecem condições à individualização de diversos domínios magnéticos, perfeitamente correlacionáveis litológico e estruturalmente com as unidades geotectônicas que compõem o arcabouço regional, assim como suas subunidades. A execução de modelamento simples, quer para as informações gravimétricas como magnetométricas, objetivaram fornecer subsídios semi-quantitativos à interpretação, de forma a auxiliar no estabelecimento do padrão geométrico do embasamento da faixa de dobramentos Araguaia, assim como das estruturas a ela impostas. Assim, por meio do processamento das informações gravimétricas, verificou-se um padrão geométrico na forma de fatias imbricadas, estabelecidas sobre o substrato. A grande anomalia negativa da porção centro-meridional da faixa pode estar relacionada a uma deficiência de massa profunda (depressão na base da crosta) ou ainda a uma associação entre contraste de densidade e espessamento do pacote de supracrustais, nessa porção. O baixo gravimétrico da porção sudeste da área, entre Guarai e Tupirama, é devido a uma depressão do embasamento, da mesma forma como a anomalia negativa associada ao Lineamento Carajás está relacionada a uma estrutura sinformal conhecida.

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.

Desenvolvimento de um conjunto de modelagem de alta resolução para os métodos eletromagnéticos indutivos

A Modelagem Eletromagnética Analógica tem sido uma poderosa ferramenta no estudo das respostas eletromagnéticas de uma grande variedade das técnicas de levantamentos e de modelos geológicos. Nas aplicações de métodos eletromagnéticos nas investigações de contaminação no subsolo, recentemente em grande demanda, são obtidas anomalias muito baixas (da ordem de 0,01% do campo primário). Com o objetivo de estudar em laboratório os problemas de levantamentos eletromagnéticos para alvos de baixo número de indução, um conjunto de modelagem em escala reduzida de alta sensibilidade, foi projetado e construído. Este conjunto opera na faixa de freqüência de 1 kHz até 1 MHz com excelente linearidade. O transdutor de transmissão gera suficiente densidade de fluxo magnético com corrente de até 3 App (pico a pico) circulando nele. O sistema de recepção, isto é, o conjunto de bobinas receptoras mais o pré-amplificador acoplado apresenta sensibilidade 3 V/A/m na freqüência de 10 kHz com inclinação de 20 dB/dec de freqüência. A alta sensibilidade do sistema de recepção combinada com a boa capacidade de geração do sistema transmissor, permite alcançar a resolução de 0,025% do campo primário. O conjunto dispõe ainda, de um sistema mecânico que permite posicionar o sistema de transdutores (transmissor e receptor) em posições arbitrárias com resolução de 1 mm e precisão de 0,05 mm. A precisão no deslocamento linear do sistema de sonda numa linha de 1715 mm de comprimento é alcançada por meio de um motor de passo controlado por um microcontrolador, que é rigidamente acoplado a um fuso que converte o movimento angular do motor em movimento linear. Com este conjunto, é possível simular uma grande variedade de técnicas de levantamentos eletromagnéticos no formato Slingram, isto é, dipolo-dipolo. Além da boa resolução mecânica e boa resolução de medida eletromagnética do conjunto, todo o processo de controle e aquisição de dados é operacionalizado por meio de um computador mestre.