Modelagem sísmica via métodos das diferenças finitas: caso da bacia do Amazonas

Este trabalho tem por objetivo apresentar os resultados da modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações e múltiplas reflexões, tendo a Bacia do Amazonas como referência à modelagem. As condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas foram analisadas numericamente pelo método das diferenças finitas, visando melhor compreensão e controle da interpretação de dados sísmicos. A geologia da Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o Recente que atingem espessuras da ordem de 5 km. Os corpos de diabásio, presentes entre os sedimentos paleozóicos, estão dispostos na forma de soleiras, alcançam espessuras de centenas de metros e perfazem um volume total de aproximadamente 90000 Km³. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica o que impossibilita melhor interpretação dos horizontes refletores que se encontram abaixo destas soleiras. Para representar situações geológicas desse tipo foram usados um modelo (sintético) acústico de velocidades e um código computacional elaborado via método das diferenças finitas com aproximação de quarta ordem no espaço e no tempo da equação da onda. A aplicação dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas melhorou a compreensão sobre a propagação em meios onde existem heterogeneidades significativas, tendo como resultado boa resolução na interpretação dos eventos de reflexão sísmica em áreas de interesse. Como resultado dos experimentos numéricos realizados em meio de geologia complexa, foi observada a influência significativa das reflexões múltiplas devido à camada de alta velocidade, isto provocou maior perda de energia e dificultou a interpretação dos alvos. Por esta razão recomenda-se a integração de dados de superfície com os de poço, com o objetivo de obter melhor imagem dos alvos abaixo das soleiras de diabásio.

Estudo da geometria e cinemática das rochas sedimentares arqueanas da mina do Igarapé do Azul – Carajás-PA

Mina de Manganês do Igarapé Azul posiciona-se geologicamente no interior do feixe da Falha Carajás, na porção central do Sistema Transcorrente de Carajás. O depósito do Manganês do Azul relaciona-se a rochas sedimentares pelíticas do Membro Azul, na base da Formação Águas Claras (Arqueano), em contato discordante, acima do Grupo Grão Pará (Nogueira et al, 1995). Três frentes de lavra a céu aberto estão atualmente em andamento na área: (1) Mina Principal (Mina 1), (2) Mina 2 e (3) Mina 3. Nestes locais encontram-se excelentes afloramentos de siltitos intercalados com argilitos e arenitos finos, intercalados com níveis manganesíferos. Essas rochas estão organizadas em conjuntos de dobras e falhas normais e inversas sob deformação heterogênea, particionada em diferentes escalas. As seções geológicas realizadas nas frentes de lavra mostram a predominância de siltitos intercalados com argilitos em contato com rochas pelíticas manganesíferas e minério (bióxido de Mn). Nessas rochas são comuns estruturas primárias tipo hummocky, estratificações cruzadas, e laminações plano-paralelas. O acamamento centimétrico a métrico (em média de 30 a 50 cm ) representa a principal estrutura primária, usada como marcador de deformação, observada nas rochas. A Mina do Igarapé Azul encontra-se dividida em dois blocos, separados por falha normal com rejeito de até dezenas de metros, com o bloco norte alto em relação ao bloco sul. O bloco sul encontra-se pouco deformado, apresentado uma regularidade no acamamento que mergulha com ângulos suaves para sul, colocando a camada de minério sucessivamente em níveis mais profundos na direção S. No bloco norte o acamamento apresenta um comportamento heterogêneo. A deformação é mais expressiva nessa região, estando o nível de minério deformado por dobras e falhas inversas. Além da cinemática vertical, as falhas apresentam deslocamento conjugado destral dando a essas feições um caráter oblíquo. Essa região pode ser definida como um corredor de deformação. O corredor observado no bloco norte, de acordo com os domínios principais separados pelas falhas anteriormente descritas, possui orientação NW-SE, com aproximadamente um quilômetro de extensão, sendo caracterizado por dobras assimétricas curvilineares com eixos de mergulhos suaves (10° a 25°) para NW e SE. Essas dobras são seccionadas por falhas normais sinuosas NW-SE e/ou E-W, com baixo ângulo de mergulho (em torno de 10° a 30°), subordinadas a transcorrências destrais, gerando em escala de detalhes feições como drag folds. Observam-se ainda falhas inversas retas e/ou sinuosas NW-SE e zonas de falhas sub-verticalizadas WNW -ESE. As dobras individuais nesta área são estruturas do tipo reversas, flexurais e com geometria en echelon com orientação semelhante às dobras curvilineares: eixos com baixo ângulo de mergulho (10° a 25°) e caimento para SE. O conjunto de feições anteriormente descrito desenha, em escala quilométrica, um antiforme aberto, provavelmente resultante da acomodação do acamamento em resposta a deformação dessas falhas. O paralelismo entre feições observadas na área da Mina do Igarapé Azul e os lineamentos maiores que desenham a Falha Carajás em planta sugere uma relação com dois importantes episódios deformacionais ocorridos durante a história tectônica da Falha Carajás. As falhas normais associadas a componente direcional destral, de maior expressividade da área da mina, estariam relacionadas ao episódio de transtensão destral responsável pela instalação da Falha Carajás anterior a 2.6 Ga (Pinheiro, 1997). As dobras, as falhas de cavalgamento e as zonas de falhas sub-verticalizadas estariam relacionados a deformações sob regime de transpressão sinistral, um segundo evento atuante na região, responsável pela reativação e inversão da maioria das estruturas próximas à zona da Falha Carajás (Pinheiro, 1997; Pinheiro & Holdsworth, 2000; Lima, 2002).

Processamento interpretativo de dados magnetométricos e inversão de dados gravimétricos aplicados à prospecção de hidrocarbonetos

Apresentamos dois métodos de interpretação de dados de campos potenciais, aplicados à prospecção de hidrocarbonetos. O primeiro emprega dados aeromagnéticos para estimar o limite, no plano horizontal, entre a crosta continental e a crosta oceânica. Este método baseia-se na existência de feições geológicas magnéticas exclusivas da crosta continental, de modo que as estimativas das extremidades destas feições são usadas como estimativas dos limites da crosta continental. Para tanto, o sinal da anomalia aeromagnética na região da plataforma, do talude e da elevação continental é amplificado através do operador de continuação analítica para baixo usando duas implementações: o princípio da camada equivalente e a condição de fronteira de Dirichlet. A maior carga computacional no cálculo do campo continuado para baixo reside na resolução de um sistema de equações lineares de grande porte. Este esforço computacional é minimizado através do processamento por janelas e do emprego do método do gradiente conjugado na resolução do sistema de equações. Como a operação de continuação para baixo é instável, estabilizamos a solução através do funcional estabilizador de primeira ordem de Tikhonov. Testes em dados aeromagnéticos sintéticos contaminados com ruído pseudo-aleatório Gaussiano mostraram a eficiência de ambas as implementações para realçar os finais das feições magnéticas exclusivas da crosta continental, permitindo o delineamento do limite desta com a crosta oceânica. Aplicamos a metodologia em suas duas implementações a dados aeromagnéticos reais de duas regiões da costa brasileira: Foz do Amazonas e Bacia do Jequitinhonha. O segundo método delineia, simultaneamente, a topografia do embasamento de uma bacia sedimentar e a geometria de estruturas salinas contidas no pacote sedimentar. Os modelos interpretativos consistem de um conjunto de prismas bidimensionais verticais justapostos, para o pacote sedimentar e de prismas bidimensionais com seções verticais poligonais para as estruturas salinas. Estabilizamos a solução, incorporando características geométricas do relevo do embasamento e das estruturas salinas compatíveis com o ambiente geológico através dos estabilizadores da suavidade global, suavidade ponderada e da concentração de massa ao longo de direções preferenciais, além de vínculos de desigualdade nos parâmetros. Aplicamos o método a dados gravimétricos sintéticos produzidos por fontes 2D simulando bacias sedimentares intracratônicas e marginais apresentando densidade do pacote sedimentar variando com a profundidade segundo uma lei hiperbólica e abrigando domos e almofadas salinas. Os resultados mostraram que o método apresenta potencial para delinear, simultaneamente, as geometrias tanto de almofadas e domos salinos, como de relevos descontínuos do embasamento. Aplicamos o método, também, a dados reais ao longo de dois perfis gravimétricos sobre as Bacias de Campos e do Jequitinhonha e obtivemos interpretações compatíveis com a geologia da área.

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.

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

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

Modelagem sísmica por diferenças finitas em meios bidimensionais com difratores

Este trabalho tem por objetivo a modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações. Como parte deste estudo foi feita a análise numérica visando as condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas. Para a validação do programa de diferenças finitas foi feita a comparação cinemática com a Teoria do Raio para um modelo simples. O motivo deste estudo é ter uma melhor compreensão e controle sobre os problemas de modelagem, visando contribuir para a solução de problemas na interpretação de dados sísmicos. Segundo vários autores na literatura geológica, Derby (1877), Evans (1906), Paiva (1929) e Moura (1938). A Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o recente, atingindo espessuras da ordem de 5 km. Os corpos de diabásio, presentes entre os sedimentos paleozóicos, estão dispostos na forma de soleiras, alcançando espessuras de centenas de metros, perfazendo um volume total de 90.000 quilômetros cúbicos. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica, impossibilitando uma melhor interpretação dos horizontes refletores que se encontram abaixo das soleiras. Para representar situações geológicas desse tipo foi usado um modelo acústico de velocidades. Para o cálculo dos sismogramas foi utilizado um programa de diferenças finitas com aproximação de quarta ordem da equação da onda acústica no espaço e no tempo. As aplicações dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas têm melhorado a compreensão sobre a propagação em meios onde existem heterogeneidades significativas, tendo como resultado boa resolução na interpretação dos eventos de reflexão sísmica em áreas de interesse. Como resultado dos experimentos numéricos realizados em meio de geologia complexa, foi observado a influência significativa das múltiplas devido a camada de alta velocidade, o que faz com que haja maior perda de energia dificultando a interpretação dos alvos. Por esta razão recomendo a integração de dados de superfície com os dados de poço, com o objetivo de se obter uma melhor imagem dos alvos abaixo das soleiras de diabásio.

Avaliação preliminar do emprego de arenito zeolítico da região nordeste do Brasil como material pozolânico para cimento Portland

As zeolitas possuem atividade pozolânica normalmente sem a necessidade de ativação térmica, por isto têm sido empregadas na produção de cimento e concreto hidráulicos desde a época do império romano. Hoje em dia são utilizadas na fabricação do cimento Portland através da substituição do clinquer em percentuais que variam entre 5 e 20%, dependendo da reatividade e da finura da zeólita. Em razão disto, são muito importantes do ponto de vista econômico e ambiental, principalmente quando não necessitam de tratamento térmico para adquirirem caráter pozolânico satisfatório, porque reduzem significativamente a energia de produção do clinquer e a liberação de CO₂ proveniente tanto da descarbonatação da calcita como da combustão de combustíveis fósseis. Contudo, dados sobre reservas de zeólitas naturais são escassos e imprecisos. No Brasil, não existe conhecimento sobre depósitos naturais de zeólitas que possam ser explorados comercialmente. No nordeste do Brasil existe a ocorrência de zeolitas sedimentares relacionadas a arenitos descoberta nos anos 2000. Estes arenitos são constituídos de quartzo, argilominerais e zeolitas naturais (estilbita). O objetivo geral desse trabalho foi avaliar se esta zeólita natural presente no arenito possui atividade pozolânica satisfatória para ser empregada como adição mineral em cimentos Portland. No programa experimental o arenito zeolítico passou por beneficiamento através da remoção, por peneiramento, do quartzo e outros minerais inertes, de modo a concentrar a zeólita estilbita e com isto verificar as propriedades pozolânicas deste mineral. No estudo experimental foram empregadas as técnicas de difração de raios X, calorimetria, ensaios químicos e de determinação da atividade pozolânica em argamassas de cal hidratada e cimento Portland. Os resultados mostraram que o arenito zeolítico acelerou a hidratação do cimento Portland devido a extrema finura do material. O arenito apresentou atividade pozolânica, sendo a estilbita responsável por este comportamento. Entretanto, a reatividade foi ligeiramente inferior ao mínimo exigido para ser empregado em escala industrial como pozolana. Estudos complementares são necessários para averiguar se o tratamento térmico entre 300 °C e 500 °C pode aumentar a atividade pozolânica do arenito devido a destruição da estrutura cristalina tanto da estilbita quanto da esmectita presente no arenito.

Estudos isotópicos (Pb, O, H, S) em zonas alteradas e mineralizadas do depósito cupro-aurífero Visconde, Província Mineral de Carajás

O depósito cupro-aurífero Visconde está localizado na Província Mineral de Carajás, a cerca de 15 km a leste do depósito congênere de classe mundial Sossego. Encontra-se em uma zona de cisalhamento de direção WNW-ESE, que marca o contato das rochas metavulcanossedimentares da Bacia Carajás com o embasamento. Nessa zona ocorrem outros depósitos hidrotermais cupro-auríferos com características similares (Alvo 118, Cristalino, Jatobá, Bacaba, Bacuri, Castanha), que têm sido enquadrados na classe IOCG (Iron Oxide Copper-Gold), embora muitas dúvidas ainda existam quanto a sua gênese, principalmente no que diz respeito à idade da mineralização e fontes dos fluidos, ligantes e metais. O depósito Visconde está hospedado em rochas arqueanas variavelmente cisalhadas e alteradas hidrotermalmente, as principais sendo metavulcânicas félsicas (2968 ± 15 Ma), o Granito Serra Dourada (2860 ± 22 Ma) e gabros/dioritos. Elas registram diversos tipos de alteração hidrotermal com forte controle estrutural, destacando-se as alterações sódica (albita + escapolita) e sódico-cálcica (albita + actinolita ± turmalina ± quartzo ± magnetita ± escapolita), mais precoces, que promoveram a substituição ubíqua de minerais primários das rochas e a disseminação de calcopirita, pirita, molibdenita e pentlandita. Dados isotópicos de oxigênio e hidrogênio de minerais representativos desses tipos de alteração mostram que os fluidos hidrotermais foram quentes (410 – 355°C) e ricos em ¹⁸O (δ¹⁸O_H2O= +4,2 a 9,4‰). Sobreveio a alteração potássica, caracterizada pela intensa biotitização das rochas, a qual ocorreu concomitantemente ao desenvolvimento de foliação milonítica, notavelmente desenhada pela orientação de palhetas de biotita, que precipitaram de fluidos com assinatura isotópica de oxigênio similar à dos estágios anteriores (δ¹⁸O_H2O entre +4,8 e +7,2‰, a 355°C). Microclina e alanita são outras fases características desse estágio, além da calcopirita precipitada nos planos da foliação. A temperaturas mais baixas (230 ± 11°C), fluidos empobrecidos em ¹⁸O (δ¹⁸O_H2O = -1,3 a +3,7‰) geraram associações de minerais cálcico-magnesianos (albita + epidoto + clorita ± calcita ± actinolita) que são contemporâneas à mineralização. Valores de δ¹⁸D_H2O e δO_H2O indicam que os fluidos hidrotermais foram inicialmente formados por águas metamórficas e formacionais, a que se misturou alguma água de fonte magmática. Nos estágios tardios, houve considerável influxo de águas superficiais. Diluição e queda da temperatura provocaram a precipitação de abundantes sulfetos (calcopirita ± bornita ± calcocita ± digenita), os quais se concentraram principalmente em brechas tectônicas - os principais corpos de minério - que chegam a conter até cerca de 60% de sulfetos. Veios constituídos por minerais sódico-cálcicos também apresentam comumente sulfetos. A associação de minerais de minério e ganga indica uma assinatura de Cu-Au- Fe-Ni-ETRL-B-P para a mineralização. Os valores de δ³⁴S (-1,2 a +3,4‰) de sulfetos sugerem enxofre de origem magmática (proveniente da exsolução de magmas ou da dissolução de sulfetos das rochas ígneas pré-existentes) e precipitação em condições levemente oxidantes. Datação do minério por lixiviação e dissolução total de Pb em calcopirita forneceu idades de 2736 ± 100 Ma e 2729 ± 150 Ma, que indicam ser a mineralização neoarqueana e, a despeito dos altos erros, permite descartar um evento mineralizador paleoproterozoico. A idade de 2746 ± 7 Ma (MSDW=4,9; evaporação de Pb em zircão), obtida em um corpo granítico não mineralizado (correlacionado à Suíte Planalto) que ocorre na área do depósito, foi interpretada como a idade mínima da mineralização. Assim, a formação do depósito Visconde teria relação com o evento transpressivo ocorrido entre 2,76 e 2,74 Ga, reponsável pela inversão da Bacia Carajás e pela geração de magmatismo granítico nos domínios Carajás e de Transição. Esse evento teria desencadeado reações de devolatilização em rochas do Supergrupo Itacaiúnas, ou mesmo, provocado a expulsão de fluidos conatos salinos aprisionados em seus intertícios. Esses fluidos teriam migrado pelas zonas de cisalhamento e reagido com as rochas (da bacia e do embasamento) pelas quais se movimentaram durante a fase dúctil. As concentrações subeconômicas do depósito Visconde devem ser resultado da ausência de grandes estruturas que teriam favorecido maior influxo de fluidos superficiais, tal como ocorreu na formação dos depósitos Sossego e Alvo 118.

Aspectos texturais do magmatismo e tramas da tectônica impostas ao Granito Chaval na Zona de Cisalhamento Santa Rosa, extremo Noroeste da Província Borborema

A região Noroeste da Província Borborema apresenta uma diversidade de corpos graníticos de natureza e evolução tectônica diversificadas, do Paleoproterozoico ao Paleozoico, com maior incidência relacionada ao Neoproterozoico e alojamento em diferentes fases da orogenia Brasiliana. Um desses exemplos é o Granito Chaval, que representa um batólito aflorante próximo à costa Atlântica do Ceará e Piauí, intrusivo em ortognaisses do Complexo Granja e supracrustais do Grupo Martinópole. Ele é, em parte, coberto por depósitos cenozoicos costeiros e rochas sedimentares paleozoicas da Bacia do Parnaíba. O Granito Chaval tem como característica marcante a textura porfirítica, destacando-se megacristais de microclina, em sienogranitos e monzogranitos, e outras feições texturais/estruturais de origem magmática, Essas permitiram interpretar sua evolução como de alojamento relativamente raso do plúton, conduzido por processos de cristalização fracionada, mistura de magmas com fluxo magmático e ação gravitacional em função da diferença de densidade do magma, levando à flutuação e ascensão de megacristais de microclina no magma residual, com alojamento de leucogranitos e pegmatitos nos estágios finais da evolução deste plutonismo. Por outro lado, em toda a metade Leste do plúton, encontra-se um rico acervo de estruturas tectógenas de cisalhamento, relacionada à implantação da Zona de Cisalhamento Transcorrente Santa Rosa, que levou a transformações tectonometamórficas superpostas às feições magmáticas, as quais atingiram condições metamórficas máximas na fácies anfibolito baixo. Cartograficamente, foram individualizados três domínios estruturais em que estão presentes uma gama de variações petroestruturais do Granito Chaval, sejam feições texturais/estruturais ígneas e tectônicas. As rochas plutônicas foram deformadas e modificadas progressivamente à medida que se dirige para Leste, no qual as rochas mudam-se para tonalidades mais escuras do cinza e os processos de cominuição e recristalização dinâmica reduzem, progressivamente, a granulação grossa desses granitos bem como o tamanho dos fenocristais para dimensões mais finas, mantendo-se suas características porfiroides. Desse modo, a trama milonítica se torna evidente, acentuando-se ao atingir a porção principal da Zona de Cisalhamento Transcorrente Santa Rosa. Como principais feições estruturais, destacam-se extinção ondulante forte; encurvamento e segmentação de cristais; geminação de deformação; rotação de cristais; microbudinagem; foliação anastomosada, inclusive S-C; lineação de estiramento; formas amendoadas de porfiroclastos, fitas e folhas de quartzo e recristalização. Os produtos desses processos de cisalhamento resultam na formação de protomilonitos, milonitos e ultramilonitos. Essas faixas miloníticas representam os locais de maior concentração da deformação, por isso é possível acompanhar progressivamente suas modificações texturais e mineralógicas, configurando uma sequência clássica de deformação progressiva heterogênea, por cisalhamento simples, em condições frágil-dúctil e dúctil. O alojamento do Granito Chaval aconteceu no final do Criogeniano (aproximadamente 630 Ma) e pode ser interpretado como magmatismo sin a tardi-tectônico em relação ao evento Brasiliano. O processo de cisalhamento que gerou a Zona de Cisalhamento Transcorrente Santa Rosa se formou nos incrementos finais da deformação de uma colisão continental em um sistema de cavalgamento oblíquo, em que se edificou o Cinturão de Cisalhamento Noroeste do Ceará, devido ao extravasamento lateral de massas crustais em fluxo dúctil acontecido no final da orogenia Brasiliana no Noroeste da Província Borborema.

DesG: uma nova metodologia para intervenção e coleta de dados magnetotelúricos

Esta tese apresenta duas contribuições para a interpretação geofísica, voltadas ao Método Magnetotelúrico. A primeira trata de uma nova abordagem para a interpretação MT, denominada MÉTODO DESCRITIVO-GEOLÓGICO (DesG), em alusão à incorporação explícita de informação a priori de correlação fácil com a descrição geológica tradicional. O intérprete define por meio de elementos geométricos (pontos e linhas) o arcabouço de feições geológicas bem como fornece valores de resistividade aos corpos geológicos presumidos. O método estima a distribuição de resistividade subsuperficial em termos de fontes anômalas próximas aos elementos geométricos, ajustando as respostas produzidas por estes corpos às medidas de campo. A solução obtida fornece então informações que podem auxiliar na modificação de algumas informações a priori imprecisas, permitindo que sucessivas inversões sejam realizadas até que a solução ajuste os dados e faça sentido geológico. Entre as características relevantes do método destacam-se: (i) os corpos podem apresentar resistividade maior ou menor do que a resistividade do meio encaixante, (ii) vários meios encaixantes contendo ou não corpos anômalos podem ser cortados pelo perfil e (iii) o contraste de resistividade entre corpo e encaixante pode ser abrupto ou gradativo. A aplicação do método a dados sintéticos evidencia, entre outras vantagens, a sua potencialidade para estimar o mergulho de falhas com inclinação variável, que merece especial interesse em Tectônica, e delinear soleiras de diabásio em bacias sedimentares, um sério problema para a prospecção de petróleo. O método permite ainda a interpretação conjunta das causas do efeito estático e das fontes de interesse. A aplicação a dados reais é ilustrada tomando-se como exemplo dados do COPROD2, cuja inversão produziu soluções compatíveis com o conhecimento geológico sobre a área. A segunda contribuição refere-se a desenho de experimento geofísico. Por meio de indicadores diversos, em especial a matriz densidade de informação, é mostrado que a resolução teórica dos dados pode ser estudada, que guia o planejamento da prospecção. A otimização do levantamento permite determinar os períodos e as posições das estações de medida mais adequados ao delineamento mais preciso de corpos cujas localizações são conhecidas aproximadamente.

Fácies sedimentares da Formação Prosperança, Neoproterozóico do Sul do Escudo das Guianas, estado do Amazonas

Rochas neoproterozóicas da Formação Prosperança, cobertura sedimentar da porção sul do Escudo das Guianas, são pobremente expostas quando comparadas com o registro sedimentar fanerozóico das bacias da Amazônia Ocidental. A Formação Prosperança consiste em conglomerados, arenitos arcosianos e pelitos que preenchem grábens no embasamento no Escudo das Guianas orientados segundo NW-SE. Esta unidade é sotoposta por rochas carbonáticas da Formação Acarí que ocorre apenas em subsuperfície, compondo o Grupo Purus, embasamento sedimentar das bacias paleozóicas produtoras de óleo da Amazônia. A análise de fácies e estratigráfica realizada na região do baixo rio Negro, Estado do Amazonas, indicou que a Formação Prosperança consiste em quatro associações de fácies representativas de um sistema flúvio-deltáico: 1) prodelta/lacustre, 2) frente deltaica, 3) foreshore/shoreface e 4) planície braided distal. Camadas tabulares de pelitos de prodelta extensos por quilômetros sugerem uma bacia sedimentar ampla de provável origem lacustre ou mar restrito. Lobos de frente deltaica complexos foram alimentados por distributários braided, formados por dunas subaquosas que migravam para NW do Bloco Amazônico durante o Neoproterozóico. Arenitos gerados sob condições de fluxo oscilatório/combinado são compatíveis com depósitos de praia e de face litorânea. A Bacia Prosperança foi fracamente invertida durante o intervalo Cambriano-Ordoviciano. Posteriormente estes depósitos são recobertos erosivamente por conglomerados e arenitos com estratificações cruzadas e planares interpretadas como depósitos fluviais do tipo braided, com migração para SE, pertencentes ao Grupo Trombetas do Eo-paleozóico da Bacia do Amazonas.

Estratigrafia da sucessão sedimentar Pós-Barreiras (Zona Bragantina, Pará) com base em radar de penetração no solo

Investigação pioneira aplicando Radar de Penetração no Solo (GPR) na área da Praia do Atalaia, norte do Brasil, nos permitiu caracterizar, pela primeira vez, fácies e estratigrafia dos depósitos conhecidos informalmente como Sedimentos Pós-Barreiras (Plioceno e mais recente). Esta sucessão sedimentar recobre discordantemente o embasamento miocênico, representado pelas formações Pirabas/Barreiras. Três unidades estratigráficas foram reconhecidas. A Unidade 1, inferior, consiste em um intervalo com 6 m de espessura média dominado por reflexões pobremente definidas e de baixa amplitude, e que intergradam com reflexões de média escala dos tipos tangencial, obliquo e hummocky. A Unidade 2, intermediária, possui cerca de 9 m de espessura e inclui, principalmente, reflexões oblíquas de larga escala, cujas configurações variam de paralela, tangencial, sigmoidal a sigmoidal-complexa. A Unidade 3, superior, corresponde a um intervalo entre 3,5 e 9 m de espessura, sendo dominada por reflexões hummocky, seguidas por reflexões de média escala dos tipos oblíquo, paralelo a sub-paralelo, e em corte-e-preenchimento. A análise da configuração das reflexões internas e geometria das reflexões nos leva a propor que a unidade correspondente aos Sedimentos Pós-Barreiras é mais variável faciologicamente que inicialmente imaginado, incluindo depósitos eólicos (dunas costeiras), bem como depósitos de cordão litorâneo, planície de maré, canal e mangue. Além disto, o mapeamento das três unidades descritas acima é importante para desvendar a complexidade de sedimentação versus erosão durante o Neógeno tardio no norte do Brasil.

Atenuação de múltiplas pelo método WHLP-CRS

Nas bacias sedimentares da região Amazônica, a geração e o acúmulo de hidrocarboneto estão relacionados com a presença das soleiras de diabásio. Estas rochas magmáticas intrusivas possuem grandes contrastes de impedância com as rochas sedimentares encaixantes, resultando em múltiplas externas e internas, com amplitudes semelhantes às das reflexões sísmicas primárias. Estas múltiplas podem predominar sobre as informações oriundas de interfaces mais profundas, dificultando o processamento, a interpretação e o imageamento da seção de sísmica. O objetivo da presente tese é realizar a atenuação de múltiplas em seções sintéticas fontecomum (CS), através da combinação dos métodos Wiener-Hopf-Levinson de predição (WHLP) e o do empilhamento superfície-de-reflexão-comum (CRS), aqui denominando pela sigla WHLPCRS. O operador de deconvolução é calculado com as amplitudes reais do sinal sísmico e traço-a-traço, o que consideramos como uma melhor eficiência para a operação de atenuação. A identificação das múltiplas é feita na seção de afastamento-nulo (AN) simulada com o empilhamento CRS, utilizando o critério da periodicidade entre primária e suas múltiplas. Os atributos da frente de onda, obtidos através do empilhamento CRS, são utilizados na definição de janelas móveis no domínio tempo-espaço, e usados para calcular o operador WHLP-CRS. No desenvolvimento do presente trabalho, visamos evitar a inconveniência da seção processada ZO; desenhar e aplicar operadores na configuração CS; e estender o método WHL para camadas curvas.

Modelo fractal para resistividade complexa de rochas: interpretação petrofísica e aplicação à exploração geoelétrica

Rochas contendo metálicos disseminados ou partículas de argila em ambiente natural onde soluções eletrolíticas normalmente preenchem os poros das rochas, exibem um tipo de polarização em baixas freqüências conhecido como polarização induzida. Nesta tese foi desenvolvido um novo modelo para descrever o fenômeno de polarização das rochas, não apenas em baixas freqüências, mas compreendendo todo o espectro eletromagnético, possível de utilização na prospecção geoelétrica. Este novo modelo engloba a maioria dos modelos utilizados até o momento como casos especiais, além de superar as limitações dos mesmos. Seu circuito analógico inclui uma impedância não linear do tipo r (iwtf)^-n que simula o efeito das superfícies rugosas das interfaces entre os grãos bloqueadores (partículas metálicas e/ou de argilas) e o eletrólito. A impedância de Warburg generalizada está em série com a resistência dos grãos bloqueadores da passagem de corrente e em paralelo com a impedância da dupla camada associada a essas interfaces. Esta combinação está em série com a resistência do eletrólito nas passagens dos poros bloqueados. Os canais não bloqueados são representados por uma resistência que corresponde à resistividade normal CC da rocha. A combinação desta resistência com a capacitância "global" da rocha é finalmente conectada em paralelo ao resto do circuito mencionado acima. Os parâmetros deste modelo incluem a resistividade CC (p₀), a cargueabilidade (m), três tempos de relaxação (t, Tf and T2), um fator de resistividade de grãos (δ_r), e o expoente de freqüência (η). O tempo de relaxação fractal (Tf), e o expoente de frequencia (η) estão relacionados à geometria fractal das interfaces rugosas entre os minerais condutivos (grãos metálicos e/ou partículas de argila bloqueando os canais dos poros) e o eletrólito. O tempo de relaxação (T) é um resultado da relaxação em baixa freqüência das duplas camadas elétricas formadas nas interfaces eletrólito-cristais, enquanto (T₀) é o tempo de relaxação macroscópico da amostra como um todo. O fator de resistividade dos grãos (δ_r) relaciona a resistividade dos grãos condutivos com o valor de resistividade CC da rocha. A resistividade CC da rocha (p₀), e δ_r estão relacionados à porosidade, à condutividade do eletrólito e às relações mineralógicas entre a matriz e os grãos condutivos. O modelo foi testado sobre um intervalo largo de freqüências contra dados experimentais de amplitude e fase da resistividade bem como para dados de constante dielétrica complexa. Os dados utilizados neste trabalho foram obtidos a partir da digitalização de dados experimentais publicados, obtidos por diversos autores e englobando amostras de rochas sedimentares, ígneas e metam6rficas. É mostrado neste trabalho que os parâmetros deste modelo permitem identificar diferenças texturais e mineralógicas nas rochas. Bote modelo foi introduzido, primeiramente, como propriedade intrínseca de um semiespaço homogêneo sendo demonstrado, neste trabalho, que a resposta observada em superfície reflete as propriedades intrínsecas do meio polarizável, sendo o acoplamento eletromagnético desprezível em freqüências menores que 104 Hz. Em seguida, o meio polarizável foi embebido em um pacote de N camadas sendo demonstrado que os parâmetros fractais do meio polarizável podem ser obtidos do levantamento em superfície para diferentes espessuras dessa camada. Isto justifica a utilização pura e simples de modelos de polarização desenvolvidos para amostras em laboratório para ajustar dados de campo, o que vem sendo feito sem uma justificativa bem fundamentada. Estes resultados demonstram a importância para a prospecção geolétrica do modelo proposto nesta tese.

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.