Aplicação de sondagem elétrica e GPR no mapeamento geotécnico no município de Belém-PA

Para a construção de edifícios o engenheiro necessita de informações geotécnicas do subsolo de rasa e média profundidades para projetar fundações. Com a implantação do projeto de Macrodrenagem da Bacia do Una de Belém, áreas que foram beneficiadas terão, num futuro próximo, um grande potencial para a ocupação e a construção de edifícios devido sua localização privilegiada. Por isso, essa área foi escolhida para investigação geofísica com o objetivo de mapeamento das camadas litológicas e consequentemente o topo da camada rija, onde geralmente se posicionam as fundações prediais profundas. Os métodos geofísicos utilizados foram de sondagem elétrica com arranjo Schlumberger, na determinação da resistividade aparente (p_α) em subsuperfície, e um eletromagnético, o radar de penetração no solo, através de um perfil, como uma ferramenta auxiliar das informações. As curvas de resistividade aparente obtidas com as sondagens Schlumberger foram processadas por inversão considerando um modelo de camadas planas, horizontais, homogêneas e isotrópicas e sua interpretação foi auxiliada por informações litológicas obtidas por sondagens a percussão com ensaio SPT. Assim, uma vez determinadas as variações de resitividades, foram confeccionadas três seções geoelétricas-litológicas e um mapa de contorno de profundidade do topo da camada rija na área estudada, com a finalidade de auxiliar o engenheiro em projetos de prédios nesta área.

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.

Determinação estatística dos contatos fluidos em perfis de poço

No interior dos reservatórios petrolíferos, o efeito da gravidade produz naturalmente a segregação dos fluidos presentes. Em razão basicamente das forças capilares não existe uma superfície definida de separação entre cada dois fluidos. No entanto, é comum na engenharia de petróleo admitir-se uma superfície convencional de separação dos fluidos. À essa superfície dá-se o nome de interface fluida ou contato fluido. A obtenção da localização dos contatos fluidos no interior dos reservatórios, como os contatos óleo-água, contatos gás-água e os contatos gás-óleo, sem a interpretação completa dos perfis geofísicos ou a realização dos procedimentos diretos dos testes de formação é um problema de difícil solução para a indústria do petróleo. O posicionamento em profundidade dos contatos pode ser utilizado para a definição dos procedimentos de completação do poço e para o posicionamento do canhoneio nos horizontes produtores dos poços verticais. Temos também a aplicação do mapeamento deste parâmetro, para os projetos e controle da perfuração no caso dos poços direcionais e horizontais, que utilizam tal informação para o posicionamento do eixo do poço no interior da zona de hidrocarboneto no reservatório, a fim de se evitar, principalmente, a produção de água. Este trabalho apresenta uma metodologia para a identificação e posicionamento dos contatos fluidos, através da aplicação da técnica da estatística multivariada, chamada de análise discriminante. Nesta aplicação, faz-se uso direto dos perfis de resistividade (R_T), de raios gama natural (RG), de densidade (ρ_b), de porosidade neutrônica (Φ_N), de caliper (Cal) e do perfil de volume de argila (Vcla), como dados de entrada no programa de análise discriminante. Para o caso de depósitos clásticos, com seqüências arenito-folhelho, a análise discriminante fornece ainda a indicação da litologia e as espessuras aparentes dos reservatórios. As curvas de perfis geofísicos aplicadas para a avaliação desta metodologia são provenientes de poços perfurados na bacia do lago Maracaibo, na Venezuela, onde, segundo estudos geológicos, existe a ocorrência de seqüências estratigráficas, com camadas de folhelhos e arenitos, sendo também confirmada a presença de hidrocarbonetos nos poços utilizados neste trabalho.

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.

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.

Espalhamento geométrico em modelos plano-estratificados

A medição de parâmetros físicos de reservatórios se constitui de grande importância para a detecção de hidrocarbonetos. A obtenção destes parâmetros é realizado através de análise de amplitude com a determinação dos coeficientes de reflexão. Para isto, faz-se necessário a aplicação de técnicas especiais de processamento capazes de corrigir efeitos de divergência esférica. Um problema pode ser estabelecido através da seguinte questão: Qual o efeito relativamente mais importante como responsável pela atenuação de amplitudes, o espalhamento geométrico ou a perda por transmissividade? A justificativa desta pergunta reside em que a correção dinâmica teórica aplicada a dados reais visa exclusivamente o espalhamento geométrico. No entanto, a análise física do problema por diferentes direções põe a resposta em condições de dúvida, o que é interessante e contraditório com a prática. Uma resposta embasada mais fisicamente pode dar melhor subsídio a outros trabalhos em andamento. O presente trabalho visa o cálculo da divergência esférica segundo a teoria Newman-Gutenberg e corrigir sismogramas sintéticos calculados pelo método da refletividade. O modelo-teste é crostal para que se possa ter eventos de refração crítica além das reflexões e para, com isto, melhor orientar quanto à janela de aplicação da correção de divergência esférica o que resulta em obter o então denominado “verdadeiras amplitudes”. O meio simulado é formado por camadas plano-horizontais, homogêneas e isotrópicas. O método da refletividade é uma forma de solução da equação de onda para o referido modelo, o que torna possível um entendimento do problema em estudo. Para se chegar aos resultados obtidos foram calculados sismogramas sintéticos através do programa P-SV-SH desenvolvido por Sandmeier (1998), e curvas do espalhamento geométrico em função do tempo para o modelo estudado como descrito por Newman (1973). Demonstramos como uma das conclusões que a partir dos dados do modelo (velocidades, espessuras, densidades e profundidades) uma equação para a correção de espalhamento geométrico visando às “verdadeiras amplitudes” não é de fácil obtenção. O objetivo maior então deveria ser obter um painel da função de divergência esférica para corrigir as verdadeiras amplitudes.

Estudo geofísico regional sobre águas subterrâneas na Ilha do Marajó-Pará-Brasil

Nos últimos dez anos foram realizadas na parte leste da Ilha de Marajó (região dos campos naturais) pelo IDESP e NCGG, mais de 800 SEVs para fins hidrogeológicos. Na época, grande parte dessas SEVs não foram totalmente interpretadas em forma quantitativa, devido à falta de recursos técnicos para fazê-lo de forma eficiente. Agora, usando meios mais modernos para interpretação automática de SEVs, voltou-se a interpretá-las com a finalidade de apresentar uma visão regional dos principais aquíferos da área, agrupar as SEVs em famílias características, testar até que ponto essa interpretação é confiável e propor o modelamento bidimensional como técnica alternativa para interpretar as SEVs realizadas em certos locais da área em questão. Como resultado dessa interpretação, com base na teoria convencional dos meios estratificados, foram definidos três tipos de sistemas de aquíferos. 1. O primeiro, denominado de aquífero profundo, situado a profundidades maiores que 50m, estende-se por toda a região prospectada, estando provavelmente associada às camadas superiores da Formação Marajó ou às litologias altamente resistivas das camadas mais profundas do Grupo Pará. 2. O segundo, denominado de aquífero raso e de média profundidade, localiza-se na parte sul e sudeste da região a profundidades compreendidas entre 10 a 50m, e está associado às lentes arenosas do Grupo Pará. 3. O terceiro, é constituído pelos paleocanais e estruturas similares, distribuídos aleatoriamente na região a pouca profundidade. A partir do estudo detalhado das SEVs, decidiu-se classificá-las em 3 famílias características com seus respectivos tipos e apresentar mapas de localização e da espessura dos aquíferos, bem como mapas de condutância longitudinal total e resistividade média da área. Estes últimos, permitem que se divida a região dos campos da Ilha de Marajó em três zonas principais: 1. Uma, altamente resistiva, situada ao sul e sudeste, a qual coincide com os terrenos aflorantes do Grupo Pará. 2. Outra, altamente condutiva, está localizada no centro e norte, onde se encontram aleatoriamente distribuídos os paleocanais e coincide com os terrenos topograficamente mais baixos, geralmente argilosos e embebidos de água salgada, que são procedentes da erosão dos terrenos circundantes topograficamente mais altos. 3. A última é medianamente resistiva e está relacionada com os terrenos vizinhos à cidade de Chaves (noroeste da região dos campos), os quais apresentam semelhanças com os do sul e sudeste da área. Usando-se a técnica de inversão na interpretação de uma SEV característica de cada família, testou-se, através do seu tratamento estatístico, até que ponto os modelos usados na interpretação dessas SEVs (teoria convencional dos meios estratificados) seriam confiáveis. Conclui-se, então, que a alta correlação existente entre os parâmetros dos modelos assumidos (camadas horizontais, isotrópicas e homogêneas) pode-se dever à utilização de modelos geofísicos muito simples para interpretar a complexa geologia de Marajó. Tendo-se verificado que nem sempre é possível aplicar a teoria das SEVs em meios horizontalmente estratificados para interpretar SEVs obtidas em certos locais de Marajó, os quais muitas vezes apresentam bruscas variações laterais de resistividade, passou-se a demonstrar que estas variações laterais afetam profundamente os dados das SEVs, utilizando-se para isto a técnica dos elementos finitos, a qual leva em conta essa variação bidimensional das propriedades físicas do meio. Foi também possível com esta técnica, modelar uma estrutura rasa, semelhante a um paleocanal, concluindo-se que estes resultados sugerem o emprego, duma forma mais profunda, deste tipo de tratamento para os dados obtidos na região dos campos da Ilha de Marajó.

Estudo da hidrogeologia no município de Palestina do Pará usando o método de eletroresistividade

Esse trabalho de geofísica aplicada à hidrogeologia, foi realizado numa área localizada na porção sudeste do estado do Pará, mas precisamente na cidade de Palestina à margem esquerda do rio Araguaia. Ironicamente existem problemas de abastecimento de água potável para a população da cidade. Esse problema deve-se ao fato que, no momento, é antieconômico o tratamento e transporte de água do rio para a cidade. Considerando que a extração de água subterrânea é economicamente mais viável para a solução desse problema, o Departamento de Geofísica e o curso de Pós-graduação de Geofísica da Universidade Federal do Pará (UFPa), por solicitação da prefeitura daquele município e da Fundação Nacional de Saúde (FNS), participaram no levantamento geofísico de eletroresistividade para investigar os possíveis locais de ocorrência de água subterrânea. Foram realizadas, na área, 21 Sondagens Elétricas Verticais (SEVs) medidas na superfície, utilizando o arranjo Schlumberger com abertura máxima dos eletrodos AB e MN de 500m e 50m, respectivamente. Para a sua interpretação foram usados os programas de modelagem direta (EGSLIB/SEV1D) e inversa (EGSLIB/SEV1DIN) de RIJO (1994), de modelos de n camadas horizontais, homogêneas e isotrópicas. Com base nos resultados da interpretação das 21 SEVs, e também das informações da geologia local foram confeccionadas três mapas de contorno e seis secções geoelétricas. Dos seis horizontes que aparecem nessas secções, o mais favorável para a exploração hidrogeológica é o quarto, com espessura variando de 10 a 79 m, estando a uma profundidade média de 50m e com valores de resistividade variando de 150 — 850 Ωm, sendo correlacionavel com a Formação Potí. A analise final dos mapas e das secções geoelétricas, indicou que o local mais apropriado para a perfuração de poços, é uma área próximo à atual área de captação.

“Eu voltei, agora pra ficar!”: um estudo sobre as estratégias empreendedoras adotadas por dekasseguis retornados no Estado do Pará

Consoante a importância da atividade empresarial para o desenrolar das migrações de retorno, sobretudo no que se refere a mitigação dos riscos e problemas tipicamente vivenciados pelos protagonistas de tais deslocamentos, esta tese objetivou compreender e explicar o processo de formação das estratégias empreendedoras adotadas por dekasseguis retornados no estado do Pará. Para tanto, desenvolveu-se um estudo de caso eminentemente qualitativo, utilizando-se um conjunto amplo de dados, com destaque para as narrativas obtidas através de entrevistas semiestruturadas, procedidas não somente com indivíduos que criaram negócios após regressarem do Japão, mas também com gestores de organizações envolvidas com o movimento dekassegui em um sentido mais amplo. Os achados obtidos permitiram elaborar um modelo de análise bastante amplo e de caráter multiescalar, que promove um entendimento mais denso sobre as trajetórias migratórias e as iniciativas empresariais observadas, capaz de captar uma série de elementos que se confluem no processo de criação e desenvolvimento das estratégias em vista. Através dele, foi possível defender a tese que: as estratégias empresariais adotadas pelos dekasseguis retornados foram formadas por combinações entre oportunidades e restrições de caráter econômico, sociocultural e político/institucional existentes em cada etapa de seus trajetos migratórios, bem como, pelo contexto espacial e temporal em que suas escolhas e ações ocorreram. O que reforça, por sua vez, a necessidade de se utilizar abordagens sociológicas e métodos complementares no estudo desses fenômenos. Além, é claro, de proporcionar um conjunto de considerações e recomendações interessantes aos estudiosos do tema e os praticantes do campo: aspirantes a dekasseguis, dekasseguis retornados que já tenham iniciado seus próprios negócios, líderes de associações nikkeys, gestores de órgãos federais e estaduais, entre outros.

Variação diária e caracterização morfológica das larvas de peixes do complexo estuarino do Rio Amazonas - PA

A grande importância dos recursos pesqueiros para a Amazônia, aliada à necessidade de ampliar os conhecimentos básicos sobre identificação das larvas de peixes (coletadas em ambiente natural), justifica o desenvolvimento deste trabalho, que tem como objetivo expandir as informações sobre o ictioplâncton, relacionando as com as tendências de variação diária e entre marés, do complexo estuarino do rio Amazonas – PA. As coletas foram realizadas durante o período diurno e noturno, no segundo semestre de 2007, pelo Projeto PIATAM mar II, sob ponto fixo na subárea 1 (estuário do rio Paracauari) e na subárea 2 (baía do Guajará) nas marés de sizígia e quadratura, em arrastos horizontais na sub-superfície da coluna d‟água com rede de plâncton cônico-cilíndrica e malha de 300μm. As amostras foram acondicionadas em recipientes contendo formalina a 4%. Os fatores hidrológicos foram obtidos in situ pelo Grupo de Oceanografia Química do Museu Paraense Emilio Goeldi. As amostras foram triadas e identificadas por meio de características morfológicas, morfométricas e merísticas, baseando-se na técnica de sequência regressiva de desenvolvimento e em bibliografias especializadas. As principais estruturas e características das fases iniciais dos peixes foram descritas e ilustradas, facilitando assim futuros estudos ictioplanctônicos para região. A temperatura superficial da água, potencial hidrogeniônico e oxigênio dissolvido não apresentaram diferenças significativas nas áreas estudadas. Os valores de salinidade não apresentaram diferença significativa entre as estações de coleta e marés, registrando apenas variação horizontal com aumento gradativo em direção à foz com valores máximos (12) e mínimos (0) para as subárea 1 e subárea 2, respectivamente. As maiores densidades de ovos foram registradas na subárea 1, em relação à subárea 2, com as maiores densidades para o período diurno (163,29 ovos/100m³) na subárea 1 e noturno (19,70 ovos/100m³) na subárea 2. As larvas foram distribuídas em 22 taxa representados por 13 famílias e 21 espécies, sendo os taxa dominantes: P. flavipinnis (46,29%), R. amazonica (19,75%), Engraulidae (10,70%), P. squamosissimus (7,55%), A. lineatus (5,19%), O. saurus (3,30%) e Gobiosoma sp. (2,15%), com elevada participação relativa dos Clupeiformes (76,75%). Quanto aos estágios de desenvolvimento, foi observada maior abundância de larvas em pré-flexão nas subáreas 1 e 2, sendo o estágio larval vitelino e pós-flexão os menos representativos. O período noturno apresentou as maiores densidade de larvas e número de taxa, evidenciando uma possível migração nictemeral do ictioplâncton. Apenas M. furnieri apresentou abundância significativamente maior nas amostras diurnas. A grande maioria dos taxa não apresentaram diferenças significativas entre as abundâncias diurnas e noturnas. Logo, a densidade de larvas e o número de taxa diferem entre o período diurno e noturno e entre maré. Portanto, as características morfológicas descritas no presente trabalho permitem uma adequada identificação das larvas, ampliando o conhecimento biológico das espécies estuarinas do litoral paraense, uma vez que as informações sobre larvas de peixes ainda são escassas, fazendo-se necessária uma intensificação nas pesquisas. Além disso, a compreensão da ecologia dos organismos, sobre tudo daqueles que apresentam seu ciclo de vida associado aos estuários, e as variações no transporte das larvas entre os períodos do dia e da noite e entre as marés são questões fundamentais para aprimorar o manejo e a conservação destes recursos renováveis.

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 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 quantitativa de anomalias gravimétricas ar-livre da margem continental norte brasileira, setor nordeste

Anomalias gravimétricas ar-livre de perfis perpendiculares a margem continental do tipo passiva apresentam uma configuração padrão. Esta configuração é, satisfatoriamente, explicada por um modelo geofísico formado por uma distribuição de descontinuidades horizontais bidimensionais. Um processo automático de busca aleatória é proposto para a interpretação quantitativa dos dados. Através do método de poliedros flexíves (Simplex), os parâmetros principais do modelo - o contraste de densidade, a profundidade, o rejeito e a localização de cada descontinuidade, puderam ser encontrados, admitindo uma relação número de pontos/número de parâmetros, a determinar, conveniente. Sobre a região do talude, as anomalias ar-livre da margem continental podem ser explicadas por uma única descontinuidade horizontal (degrau simples); e tendo que a resposta dos dados gravimétricos no domínio do número de onda contém informações sobre esta anomalia, foi proposto um procedimento gráfico iterativo para a análise espectral deste sinal. Aplicando a transformada de Fourier é possível determinar a profundidade e o rejeito da descontinuidade, e conhecendo estes parâmetros a densidade é calculada unicamente. O objetivo básico do uso destes procedimentos seria combinar os dois métodos de interpretação nos domínios do espaço e do número de onda, com a finalidade de obter soluções vinculadas mais plausíveis quanto ao contexto geológico esperado para a área estudada. Os dois procedimentos de interpretação foram aplicados nas anomalias gravimétricas ar-livre da margem continental norte brasileira, setor nordeste, abrangendo os estados do Maranhão ao Rio Grande do Norte. As respectivas capacidade de resolução de cada procedimento foram então analisadas. Demonstrou-se que a inversão realizada diretamente no domínio do espaço é mais favorável na interpretação das anomalias ar-livre, embora o tratamento espectral seja relativamente mais simples.

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

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

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

A presente Dissertação de Mestrado tem como objetivo o estudo do problema de inversão sísmica baseada em refletores planos para arranjo fonte-comum (FC) e ponto-médiocomum (PMC). O modelo direto é descrito por camadas homogêneas, isotrópicas com interfaces plano-horizontais. O problema é relacionado ao empilhamento NMO baseado na otimização da função semblance, para seções PMC corrigidas de sobretempo normal (NMO). O estudo foi baseado em dois princípios. O primeiro princípio adotado foi de combinar dois grupos de métodos de inversão: um Método Global e um Método Local. O segundo princípio adotado foi o de cascata, segundo a teoria Wichert-Herglotz-Bateman, que estabelece que para conhecer uma camada inferior tem-se que conhecer primeiro a camada superior (dissecação). A aplicação do estudo é voltada à simulação sísmica de Bacia Sedimentar do Solimões e de Bacia Marinha para se obter uma distribuição local 1D de velocidades e espessuras para a subsuperfície em horizontes alvo. Sendo assim, limitamos a inversão entre 4 e 11 refletores, uma vez que na prática a indústria limita uma interpretação realizada apenas em número equivalente de 3 a 4 refletores principais. Ressalta-se que este modelo é aplicável como condição inicial ao imageamento de seções sísmicas em regiões geologicamente complexas com variação horizontal suave de velocidades. Os dados sintéticos foram gerados a partir dos modelos relacionados a informações geológicas, o que corresponde a uma forte informação a priori no modelo de inversão. Para a construção dos modelos relacionados aos projetos da Rede Risco Exploratório (FINEP) e de formação de recursos humanos da ANP em andamento, analisamos os seguintes assuntos relevantes: (1) Geologia de bacias sedimentares terrestre dos Solimões e ma rinha (estratigráfica, estrutural, tectônica e petrolífera); (2) Física da resolução vertical e horizontal; e (3) Discretização temporal-espacial no cubo de multi-cobertura. O processo de inversão é dependente do efeito da discretização tempo-espacial do campo de ondas, dos parâmetros físicos do levantamento sísmico, e da posterior reamostragem no cubo de cobertura múltipla. O modelo direto empregado corresponde ao caso do operador do empilhamento NMO (1D), considerando uma topografia de observação plana. O critério básico tomado como referência para a inversão e o ajuste de curvas é a norma 2 (quadrática). A inversão usando o presente modelo simples é computacionalmente atrativa por ser rápida, e conveniente por permitir que vários outros recursos possam ser incluídos com interpretação física lógica; por exemplo, a Zona de Fresnel Projetada (ZFP), cálculo direto da divergência esférica, inversão Dix, inversão linear por reparametrização, informações a priori, regularização. A ZFP mostra ser um conceito út il para estabelecer a abertura da janela espacial da inversão na seção tempo-distância, e representa a influência dos dados na resolução horizontal. A estimativa da ZFP indica uma abertura mínima com base num modelo adotado, e atualizável. A divergência esférica é uma função suave, e tem base física para ser usada na definição da matriz ponderação dos dados em métodos de inversão tomográfica. A necessidade de robustez na inversão pode ser analisada em seções sísmicas (FC, PMC) submetida a filtragens (freqüências de cantos: 5;15;75;85; banda-passante trapezoidal), onde se pode identificar, comparar e interpretar as informações contidas. A partir das seções, concluímos que os dados são contaminados com pontos isolados, o que propõe métodos na classe dos considerados robustos, tendo-se como referência a norma 2 (quadrados- mínimos) de ajuste de curvas. Os algoritmos foram desenvolvidos na linguagem de programação FORTRAN 90/95, usando o programa MATLAB para apresentação de resultados, e o sistema CWP/SU para modelagem sísmica sintética, marcação de eventos e apresentação de resultados.