Modelagem numérica de dados MCSEM 3D usando computação paralela

Desenvolvemos a modelagem numérica de dados sintéticos Marine Controlled Source Electromagnetic (MCSEM) usada na exploração de hidrocarbonetos para simples modelos tridimensionais usando computação paralela. Os modelos são constituidos de duas camadas estrati cadas: o mar e o sedimentos encaixantes de um delgado reservatório tridimensional, sobrepostas pelo semi-espaço correspondente ao ar. Neste Trabalho apresentamos uma abordagem tridimensional da técnica dos elementos nitos aplicada ao método MCSEM, usando a formulação da decomposição primária e secundária dos potenciais acoplados magnético e elétrico. Num pós-processamento, os campos eletromagnéticos são calculados a partir dos potenciais espalhados via diferenciação numérica. Exploramos o paralelismo dos dados MCSEM 3D em um levantamento multitransmissor, em que para cada posição do transmissor temos o mesmo processo de cálculos com dados diferentes. Para isso, usamos a biblioteca Message Passing Interface (MPI) e o modelo servidor cliente, onde o processador administrador envia os dados de entradas para os processadores clientes computar a modelagem. Os dados de entrada são formados pelos parâmetros da malha de elementos nitos, dos transmissores e do modelo geoelétrico do reservatório. Esse possui geometria prismática que representa lentes de reservatórios de hidrocarbonetos em águas profundas. Observamos que quando a largura e o comprimento horizontais desses reservatório têm a mesma ordem de grandeza, as resposta in-line são muito semelhantes e conseqüentemente o efeito tridimensional não é detectado. Por sua vez, quando a diferença nos tamanhos da largura e do comprimento do reservatório é signi cativa o efeito 3D é facilmente detectado em medidas in-line na maior dimensão horizontal do reservatório. Para medidas na menor dimensão esse efeito não é detectável, pois, nesse caso o modelo 3D se aproxima de um modelo bidimensional. O paralelismo dos dados é de rápida implementação e processamento. O tempo de execução para a modelagem multitransmissor em ambiente paralelo é equivalente ao tempo de processamento da modelagem para um único transmissor em uma máquina seqüêncial, com o acréscimo do tempo de latência na transmissão de dados entre os nós do cluster, o que justi ca o uso desta metodologia na modelagem e interpretação de dados MCSEM. Devido a reduzida memória (2 Gbytes) em cada processador do cluster do departamento de geofísica da UFPA, apenas modelos muito simples foram executados.

Uma nova abordagem para interpretação de anomalias gravimétricas regionais e residuais aplicada ao estudo da organização crustal: exemplo da Região Norte do Piauí e Noroeste do Ceará

A interpretação de anomalias gravimétricas é de grande importância no estudo de feições geológicas que ocorrem na crosta terrestre. Esta interpretação é, no entanto, dificultada pelo fato das anomalias gravimétricas serem resultantes da soma total dos efeitos produzidos por todos os contrastes de densidades de subsuperfície. Desse modo, com o intuito de separar efeitos de feições mais profundas de efeitos de feições mais rasas, bem como a caracterização da geometria desses dois conjuntos de feições, apresentamos um método de separação das componentes regional e residual do campo e a subsequente interpretação de cada componente. A separação regional-residual de dados gravimétricos é efetuada através da aproximação do campo regional por um polinômio ajustado ao campo observado por um método robusto. Este método é iterativo e usa como aproximação inicial a solução obtida através do ajuste polinomial pelo método dos mínimos quadrados. O método empregado minimiza a influência de observações contendo forte contribuição do campo residual no ajuste do campo regional. A componente regional obtida a partir da separação regional-residual é transformada em um mapa de distâncias verticais em relação a um nível de referência. Esta transformação compreende duas etapas. A primeira consiste na obtenção da continuação para baixo da componente regional, que é pressuposta ser causada por uma interface suave separando dois meios homogêneos, representando a interface crosta-manto, cujo contraste de densidade é supostamente conhecido. A segunda consiste na transformação do mapa de continuação para baixo em um mapa de distâncias verticais entre o nível de continuação (tomado como nível de referência) e a interface. Este método apresenta duas dificuldades. A primeira está ligada à instabilidade, havendo portanto a necessidade do emprego de um estabilizador o que acarreta a perda de resolução das feições que se desejam mapear. A segunda, inerente ao método gravimétrico, consiste na impossibilidade da determinação das profundidades absolutas da interface em cada ponto, bastando entretanto o conhecimento da profundidade absoluta em um ponto, através de informação independente, para que todas as outras profundidades absolutas sejam conhecidas. A componente residual obtida a partir da separação regional-residual é transformada em um mapa de contrastes de densidade aparente. Esta transformação consiste no cálculo do contraste de densidade de várias fontes prismáticas através de uma inversão linear pressupondo que as fontes reais estejam das a uma placa horizontal, com contrastes de densidade variando apenas nas direções horizontais. O desempenho do método de separação regional-residual apresentado foi avaliado, através de testes empregando dados sintéticos, fornecendo resultados superiores em relação aos métodos dos mínimos quadrados e da análise espectral. O método de interpretação da componente regional teve seu desempenho avaliado em testes com dados sintéticos onde foram produzidos mapeamentos de interfaces bem próximas das estruturas reais. O limite de resolução das feições que se desejam mapear depende não só do grau do polinômio ajustante, como também da própria limitação inerente ao método gravimétrico. Na interpretação da componente residual é necessário que se postule ou tenha informação a priori sobre a profundidade do topo e espessura da placa onde as fontes estão supostamente confinadas. No entanto, a aplicação do método em dados sintéticos, produziu estimativas razoáveis para os limites laterais das fontes, mesmo na presença de fontes interferentes, e pressupondo-se valores para profundidade do topo e espessura da placa, diferentes dos valores verdadeiros. A ambiguidade envolvendo profundidade do topo, espessura e densidade pode ser visualizada através de gráficos de valores de densidade aparente contra profundidade do topo presumida para a placa para vários valores postulados para a espessura da placa. Estes mesmos gráficos permitem, pelo aspecto das curvas, a elaboração de uma interpretação semi-quantitativa das profundidades das fontes reais. A seqüência dos três métodos desenvolvidos neste trabalho foi aplicada a dados gravimétricos da região norte do Piauí e noroeste do Ceará levando a um modelo de organização crustal que compreende espessamentos e adelgaçamentos crustais associados a um evento compressivo que possibilitou a colocação de rochas densas da base da crosta a profundidades rasas. Este modelo ê compatível com os dados geológicos de superfície. É ainda sugerida a continuidade, por mais 200 km em direção a sudoeste, do Cinturão de Cisalhamento Noroeste do Ceará por sob os sedimentos da Bacia do Parnaíba, com base nas evidências fornecidas pela interpretação da anomalia residual. Embora esta seqüência de métodos tenha sido desenvolvida com vistas ao estudo de feições crustais de porte continental, ela também pode ser aplicada ao estudo de feições mais localizadas como por exemplo no mapeamento do relevo do embasamento de/bacias sedimentares onde os sedimentos são cortados por rochas intrusivas mais densas.

Análise comparativa de três métodos de inversão aplicados a dados gravimétricos e magnéticos em perfil

Tradicionalmente, o método dos mínimos quadrados tem sido empregado na inversão não linear de dados de campo potencial. No caso em que as observações dos campos gravimétrico ou magnético contém apenas ruído Gaussiano. O método dos mínimos quadrados não apresenta problemas. Entretanto, quando as observações são perturbadas por ruído não Gaussiano, ou mesmo por ruído não aleatório, como é o caso de muitos ruídos geológicos, o método dos mínimos quadrados torna-se bastante ineficiente, e métodos alternativos devem ser empregados a fim de produzir interpretações realísticas. Neste trabalho, uma comparação é feita entre os métodos dos mínimos quadrados, dos mínimos absolutos e do ajuste-M, aplicados à inversão não linear de dados de campo potencial. A comparação é efetuada usando-se dados teóricos, onde diversas situações geológicas são simuladas. Os resultados mostram que na presença de ruído geológico, caracterizado por pequeno corpo raso acima do corpo principal, ou por corpo grande, adjacente ao corpo principal, o ajuste-M apresenta desempenho muito superior ao dos mínimos quadrados e dos mínimos absolutos. Na presença de ruído Gaussiano, entretanto, o ajuste-M tem um desempenho inferior aos outros dois métodos. Como o ruído Gaussiano é um ruído branco, parte dele pode ser removido por um filtro passa baixa adequado, sem muita perda do sinal, o que não ocorre com o ruído geológico que contém componentes importantes de baixo número de onda. Desse modo o ajuste-M se torna uma ferramenta importante na interpretação de áreas geologicamente complexas, onde é comum a contaminação das anomalias por ruído geológico. Os três métodos em estudo são aplicados a uma anomalia magnética real causada por uma intrusão de diabásio em forma de dique, em sedimentos arenosos da formação Piauí na Bacia do Parnaíba. Os três métodos apresentaram resultados semelhantes indicando que tanto o nível de ruído Gaussiano como geológico são baixos nesta anomalia.

Aplicação dos métodos elétrico e eletricomagnético para o estudo hidrogeológico no município de Altamira-Pará

Este trabalho trata da aplicação dos métodos geofísicos eletroresistivo e eletromagnético no município de Altamira, região central do Estado do Pará, com a finalidade de fazer o mapeamento da espessura das camadas aqüíferas, bem como a ocorrência de intrusões de diabásio na área do município, permitindo uma melhor locação de pontos para perfuração com fins de captação à água subterrânea. Na área estudada, verificam-se afloramentos de rochas cristalinas do embasamento coberto por sedimentos com idade variando desde o Paleozóico ao Quaternário, com esperança de se obter água nos sedimentos do Devoniano Inferior, mais precisamente na formação Maecuru. Foram executadas 34 sondagens elétricas verticais com arranjo Schlumberger com AB/2 variando de 100 a 500 metros. 108 sondagens eletromagnéticas transientes também foram executadas com arranjo "loop coincidente", com 50 metros de lado. Por problemas eletrônicos adicionados aos ruídos de fenômenos supraparamagnéticos, os resultados das sondagens eletromagnéticas apresentaram-se muito distorcidos, não tendo sido possível extrair deles informação geológica de relevância. Por isso, foram criados modelos eletromagnéticos sintéticos representando a situação geológica da área. Os dados obtidos no campo em forma de curvas de resistividade aparente, foram plotados em papel bilogarítmo e interpretados com a ajuda de algoritmos direto e inverso os quais considera modelos de camadas horizontalmente estratificadas, homogêneas e isotrópicas. Os resultados da interpretação permitiram estabelecer modelos geoelétricos aproximados refletindo a situação geológica na subsuperfície. Os resultados da interpretação das sondagens elétricas, foram confrontados com os resultados da interpretação qualitativa do mapa de isogamas do campo magnético total, permitindo delimitar com maior precisão a zona de incidência de soleiras de diabásio, situada a norte da área. Dos quatro horizontes geoelétricos distinguidos, o terceiro é essencialmente constituído por arenitos, apresentando, por isso, melhores condições para armazenamento de água subterrânea.

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

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

Eletroresistividade aplicada ao estudo de água subterrânea no município de Salvaterra-Marajó-Pará

A partir de 1974, quando da implantação do Núcleo de Ciências Geofísicas e Geológicas da UFPa., começou a ser desenvolvido um estudo de água subterrânea na Ilha do Marajó, sob o titulo de Projeto Marajó, financiado pela FINEP, que visava definir e delinear aquíferos em algumas áreas piloto da ilha. Este trabalho, realizado em uma dessas áreas piloto, no município de Salvaterra, Marajó, busca determinar áreas favoráveis para exploração de água subterrânea, através do uso do método geofísico de eletroresistividade. A área, de 200 Km², foi prospectada através de 65 sondagens elétricas verticais e 6 perfilagens horizontais, tendo sido utilizado o dispositivo Schlumberger nas sondagens e o dispositivo Wenner nas perfilagens. Os dados foram inicialmente interpretados através do método do Ponto Auxiliar, quando se obteve uma visão geral do comportamento resistivo dos terrenos sondados, sendo depois reinterpretados usando técnica computacional, indicando existir na área camadas com resistividades que variam desde 5 até 43000 ohms-metro. Observou-se a existência de um pacote formado por uma sucessão de camadas que apresentaram resistividades com valores de 500 a 4000 ohms-metro, interpretado como constituído de sedimentos arenosos com água doce (aquífero), tendo uma espessura que chegou a atingir mais de 60 metros em algumas áreas. Na área pesquisada, foram delimitadas três sub-áreas cobrindo cerca de 40% da área total e, nestas, as espessuras do aquífero vão de 20 a 80 metros. As profundidades do topo desse aquífero nas três sub-áreas variam de 1,3 a 6,6 metros. A interpretação geofísica foi comprovada através de quatro furos de sondagem até 20 metros. Dois destes furos, realizados nas áreas interpretadas como bons aquíferos, mostraram terrenos mais arenosos, e os outros dois, alotados fora daquelas áreas, apresentaram uma parcela considerável de sedimentos argilosos com algum material orgânico. Informações obtidas quando do levantamento dos poços escavados na região e das análises de amostras das águas destes poços, mostram que, em termos gerais, não existem problemas quanto à qualidade das águas desses aquíferos mais superficiais, com exceção de dois poços escavados muito próximos da borda da Baía do Marajó, nas localidades de Joanes e Monsarás, que apresentam um aumento na salinidade quando no pico da estação de estio. Pelo que foi visto, o método utilizado neste trabalho é uma técnica eficiente para localizar bons aquíferos em ambientes geológicos semelhantes ao da área pesquisada.

Morfoestratigrafia do noroeste da planície costeira de São João de Pirabas (Nordeste do Pará)

Na Planície Costeira de São João de Pirabas (NE do Pará), ocorrem principalmente depósitos terciários e quaternários, cuja distribuição e espessura foram influenciadas por movimentos neotectônicos e oscilações do nível do mar atuantes desde o Mioceno. A análise estratigráfica com base em testemunhos (com até 6 m de comprimento), permitiu a identificação de um padrão de sedimentação, visualizado através de quatro sucessões marinhas: sucessão marinha retrogradacional basal - SB (lamas de intermaré, areias de antigos cordões praiais e areias de canais de maré); sucessão marinha retrogradacional - S1 (sedimentos predominantemente arenosos de face praial); sucessão progradacional - S2 (ambiente de planície de maré e "chêniers" associados) e; sucessão retrogradacional atual - S3 (cordões duna-praia, barras arenosas longitudinais e de baías, que migram sobre os manguezais). A evolução desta planície está relacionada às oscilações do nível do mar que, inicialmente, deram origem a sucessão Retrogradacional Basal (SB), durante uma fase transgressiva relacionada ao Pleistoceno Terminal (?), enquanto as sucessões S1, S2 e S3 teriam evoluído a partir da atuação de ciclos transgressivos e regressivos, desde aproximadamente 5.100 anos A.P até os dias atuais.

Mineralogical and geochemical influences on sediment color of Amazon wetlands analyzed by visible spectrophotometry

Com base em dados sedimentológicos e geoquímicos, este trabalho relaciona medições espectrofotométricas com a composição do sedimento, e sua aplicação em estudos paleoecológicos das áreas alagáveis da Amazônia. Os dados CIELAB estão diretamente relacionados à composição mineralógica e química dos sedimentos, especialmente quartzo, oxihidróxidos e sulfetos de ferro, e carbono orgânico total. Conteúdos de carbono orgânico total entre 0,4-1%, 1-2%, 3-5% e 15-40% foram relacionados a dados de L* (luminosidade) de 27, 26-15, 7-10 e 7 ou menos, respectivamente. Os valores CIELAB de um depósito com turfa em Marabá, Pará, foram proporcionais a variações no conteúdo de quartzo e carbono orgânico total, mas mudanças nas zonas de cores similares, principalmente nos valores de +a* (vermelho) e +b* (amarelo), ao longo de outros depósitos em Calçoene, Amapá e Soure, Pará, indicam uma relação muito próxima entre os conteúdos de carbono orgânico total, oxihidróxidos e sulfetos de ferro. Além disso, o diagrama Q7/4 (razão entre valores percentuais de refletância em 700 nm e 400 nm, juntamente com dados de L*) indicou sedimentos ricos em ferro para a fácies lama bioturbada no depósito do Amapá, fácies lama bioturbada e areia bioturbada do depósito de Soure, e das fácies areia com laminação cruzada e areia maciça do depósito de Marabá. Ainda, sedimentos ricos em carbono orgânico foram encontrados na lama bioturbada no depósito do Amapá, fácies heterolítica lenticular e lama bioturbada do depósito de Soure, e das fácies lama laminado e turfa do depósito de Marabá. Na área de Marabá, os dados sugerem uma influência autóctone com formação de turfa. As áreas de zonas úmidas costeiras no Marajó e Amapá representam o desenvolvimento típico de planícies de maré com formação de sulfetos e oxihidróxidos de ferro durante alternâncias entre inundação e exposição.

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.

Geoquímica das águas do médio e baixo rio Madeira e seus principais tributários - Amazonas - Brasil

Este trabalho teve como objetivo estudar as águas do rio Madeira e seus principais tributários entre a cidade de Humaitá e sua foz no rio Amazonas. Foram analisados pH, condutividade, turbidez, íons maiores, elementos traço e isótopos de Sr nos períodos de seca, cheia e transição para a seca entre 2009 e 2010. As águas do Madeira, classificadas com brancas, são bicarbonatadas-cálcicas, têm pH entre 5 e 6 e são mais concentradas que as dos tributários. Estes têm águas de cor preta, mais ácidas e quimicamente heterogêneas, os da margem esquerda são quimicamente mais semelhantes as do Madeira, enquanto os da margem direita têm alta concentração em SiO2. Os cátions, Cl- e NO3- são mais concentrados na cheia o que sugere influência do solo, da vegetação e da composição da água da chuva (Cl-), enquanto HCO3-, SO42-, Al, Br e P, com maiores concentrações na seca, devem estar relacionados com a química das rochas. A SiO2 e os elementos terras raras (ETR) com concentrações elevadas na seca e na cheia, estão associados tanto a vegetação e ao solo como as rochas. A interação desses fatores é a causa da heterogeneidade química das águas. Contudo, a semelhança entre as águas dos tributários da margem esquerda e as do Madeira são consequência das rochas dos Andes serem a fonte dos sedimentos cenozóicos percolados por elas, enquanto a química das águas dos tributários da margem direita retrata a estabilidade tectônica, o intenso intemperismo e a baixa taxa de erosão das rochas do cráton Amazônico.

Jarina: o marfim das biojóias da Amazônia

A palmeira jarina (Phytelephas macrocarpa) é endêmica da Amazônia, onde se desenvolve sobre antigas planícies de inundação, cujos sedimentos são constituídos por quartzo, minerais de argila 2:1 e feldspatos, constituindo solos férteis e pouco ácidos a neutros. As sementes dessa palmeira são incluídas entre as gemas orgânicas raras. Devido a sua cor e brilho, as sementes são comparadas ao marfim animal, apesar da baixa dureza e baixa densidade, sendo empregadas na manufatura de biojóias e artefatos. Esses produtos são bem aceitos comercialmente devido às sementes serem susceptíveis a mudança de coloração e outros melhoramentos. Infelizmente, as jóias não apresentam vida longa, pois as sementes podem sofrer ataque de microorganismos entre 5 e 10 anos. Se houver uma política adequada para cadeia produtiva das sementes de jarina, a mesma poderá se tornar de grande importância para o desenvolvimento da região Amazônica, ao criar novas oportunidades de trabalho e agregação de valor aos produtos. No entanto faz-se necessário um especial cuidado para evitar exploração inadequada das sementes para assegurar a preservação da espécie.

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.

Estudo das feições lineares observadas em dados aeromagnéticos da Bacia do Solimões a partir de imagens digitais de sombreamento do relevo anômalo total

A Bacia do Solimões é caracterizada por um padrão magnético complexo, expresso ao longo de toda sua extensão, e produzido pela superposição de anomalias magnéticas que se apresentam sob a forma de feições lineares. Tais feições refletem as diferentes atividades tectônicas que atuaram na Região Amazônica através do Pré-Cambriano e Fanerozóico. Neste trabalho, foi aplicado nos dados aeromagnéticos da Bacia do Solimões um método de processamento de imagens digitais de sombreamento do relevo magnético anômalo total que, dadas suas características metodológicas, permitindo utilizá-lo como um filtro direcional, possibilitou a definição de aspectos relevantes no contexto das relações entre os lineamentos magnéticos E/W, NE/SW, e NW/SE. Nesse sentido, foram identificados padrões de lineamentos magnéticos que refletem a existência de zonas de cisalhamento transcorrentes dextrais, orientadas preferencialmente na direção E/W. A interação entre os diversos segmentos transcorrentes promoveu o desenvolvimento de um regime predominantemente transpressivo, representado por falhas reversas associadas aos lineamentos magnéticos E/W e NE/SW; e duplexes direcionais formando falhas em flor positivas e negativas associadas, respectivamente, aos lineamentos magnéticos orientados na direção N70-80E e N70-80W. A análise quantitativa permitiu explicar dois aspectos importantes em relação às feições lineares observadas nas imagens digitais. O primeiro mostra, através de modelamentos baseados na superposição de prismas bidimensionais, que estas feições lineares magnéticas podem ser explicadas pela superposição de fontes profundas intraembasamento altamente magnéticas, e fontes rasas de alta frequência, sendo estas associadas a falhas reversas ao longo dos níveis de diabásio, presentes em forma de intrusões nos sedimentos paleozóicos da Bacia do Solimões. O segundo aspecto, baseado na utilização dos conceitos da cross-covariância, constata a presença de 'offsets' dextrais, associados aos lineamentos magnéticos NE—SW, ao longo da direção E—W. Este fato mostra, quantitativamente, que o padrão magnético desta região pode ser explicado pela presença de zonas de cisalhamento transcorrestes dextrais, cujos processos tectônicos associados foram fortemente condicionados por zonas de fraquezas pré-existentes (Pré-Cambrianas, Paleozóicas) durante o Mesozóico e Cenozóico.

Geofísica aplicada ao mapeamento de aquíferos da Região Metropolitana de Belém

Hoje, cerca de 46% da água potável consumida na região de Belém provém dos mananciais subterrâneos, o que mostra a grande importância desta fonte, mesmo em uma região rica em água superficial (rios, igarapés, furos e até mesmo uma baía). Isso decorre da falta de execução do planejamento previsto, além de que os custos de obras de captação de águas superficiais são muito elevados. A alternativa para Belém tem sido a água subterrânea, como um bom, saudável e econômico paliativo. A existência de uma enorme quantidade desses recursos nos sedimentos Barreiras e na Formação Pirabas contribui para esta alternativa. Os sedimentos Barreiras foram depositados em ambiente fluvial e estuarino, enquanto que a Formação Pirabas possui fácies de origem estearina, de mangue e marinha. No contexto apresentado, mostra-se então a necessidade de investigar tais recursos com o intuito de fazer ótimo proveito destes. Para tanto, foram usadas duas metodologias geofísicas de forma integrada: a Sondagem Elétrica Vertical (SEV) e a Perfilagem Geofísica de Poço, usando raios gama, potencial espontâneo e resistência elétrica. Pode-se discernir dois ambientes aquíferos. O primeiro até cerca de 160 m de profundidade, em que prevalecem corpos arenosos em formas de canal e lentes entrecortados ou não, imersos em uma matriz argilosa. O segundo é caracterizado por pacotes arenosos bastante extensos com lentes de argila, cuja profundidade é desde 160 m até 285 m, aproximadamente. O consumo de água da porção superior é adequado às residências e condomínios, visto o custo elevado de obras de captação a grandes profundidades. Já o abastecimento público deve utilizar a porção inferior, devido às vazões bastante elevadas (maiores que 100 m³/h) e devido à garantia de água de boa qualidade, protegida de ações antrópicas, como vem sendo feito.

Geofísica aplicada à prospecção de água subterrânea na área do Rio Paracauari – Ilha de Marajó - Pará

A região do Rio Paracauari apresenta sérios problemas hídricos devido o excesso de chuvas no inverno e escassez de água no verão. A indústria agropecuária, principal atividade econômica da região, sofre tremendamente durante o período de secas. Foram realizados estudos geofísicos através do método de resistividade para localizar e mapear os aquíferos rasos de boa potencialidade. Foram feitas 53 sondagens elétricas verticais, do tipo Schlumberger, e 5 perfis de resistividade do tipo Wenner. Foi utilizada a técnica do ponto auxiliar e, depois, os diagramas coletivamente foram tratados e processados no computador, para obtenção de modelos aproximados. Depois da interpretação, foi possível separar duas áreas distintas. A primeira com resistividade alta até 950 Ω.m, com predominância de areia, é constituída de excelentes aquíferos e, provavelmente, trata-se de "Paleocanais". A segunda, com resistividade baixa até 0,5 Ω.m, com predominância de sedimentos argilosos, sílticos, com água salobra. Com os mapas de resistência transversal total e de resistividade foram definidas as áreas mais favoráveis à exploração de água potável. A área favorável é somente 30% da área total. Foi feita coleta das águas rasas e superficiais da área, e suas resistividades variam de 362 Ω.m a 1,1 Ω.m. Algumas perfurações foram feitas, e vieram constatar a existência destas áreas distintas. Com isto, foram solucionados parcialmente os problemas de água de algumas fazendas como São Lourenço, Gavinho e Conceição.