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.

Batólito Guaporeí: uma extensão do Complexo Granitoide Pensamiento em Mato Grosso, SW do Cráton Amazônico

O Batólito Guaporeí é um corpo de aproximadamente 240 km² alongado segundo a direção NW, localizado na região de Vila Bela da Santíssima Trindade, estado de Mato Grosso. Situa-se nos domínios da Província Rondoniana-San Ignácio, no Terreno Paraguá, na porção meridional do Cráton Amazônico. É formado por monzogranitos e, subordinadamente, granodioritos, quartzo-monzonitos e sienogranitos, caracterizados por granulação grossa e textura, em geral, porfirítica a porfiroclástica. Possui biotita como mineral máfico primário, por vezes, associada a anfibólio, e encontra-se metamorfizado na fácies xisto verde, exibindo estrutura milonítica, em estreitas zonas de cisalhamento. Evidências geoquímicas indicam que essas rochas derivam de um magma cálcio-alcalino de alto potássio a shoshonítico, metaluminoso a levemente peraluminoso evoluído por cristalização fracionada associada à assimilação crustal, possivelmente gerado em ambiente de arco continental. Duas fases de deformação relacionadas à Orogenia San Ignácio, caracterizadas pelo estiramento e alinhamento mineral evidenciadas pelas foliações S₁ e S₂, foram identificadas nestas rochas. Foi obtida pelo método de evaporação de Pb em zircão uma idade de 1.314 ± 3 Ma, interpretada como idade de cristalização do corpo granítico. Dados Sm-Nd em rocha total indicam idade modelo T_DM em torno de 1,7 Ga e valor negativo para _{εNd (t = 1,3)} (-14), corroborando a hipótese de envolvimento crustal na gênese do magma. Os resultados obtidos apontam semelhanças entre essas rochas e aquelas de região adjacente em território boliviano, sugerindo que o Granito Guaporeí representa uma extensão do Complexo Granitoide Pensamiento.

Geocronologia e evolução crustal da área do depósito de Cu-Au Gameleira, Província Mineral de Carajás (Pará), Brasil

O depósito de Cu-Au Gameleira está hospedado nas rochas do Grupo Igarapé Pojuca, pertencente ao Supergrupo Itacaiúnas, Província Mineral de Carajás, SE do Cráton Amazônico. Esse grupo está representado principalmente por rochas metavulcânicas máficas (RMV), anfibolitos, biotita xistos, formações ferríferas e/ou hidrotermalitos, cortadas por rochas intrusivas máficas (RIM), bem como por granitos arqueanos (2,56 Ga, Granito Deformado Itacaiúnas) e paleoproterozóicos (1,87 - 1,58 Ga, Granito Pojuca e Leucogranito do Gameleira). Cristais de zircão de um saprolito (2615 ± 10 Ma e 2683 ± 7 Ma) e de uma amostra de RIM (2705 ± 2 Ma), mostraram ser contemporâneos aos dos gabros do depósito Águas Claras. Datações Pb-Pb em rocha total e calcopirita de RMV indicaram idades de 2245 ± 29 Ma e 2419 ± 12 Ma, respectivamente, enquanto lixiviados de calcopirita indicaram idades de 2217 ± 19 Ma e 2180 ± 84 Ma. Essas idades são interpretadas como rejuvenescimento parcial provocado pelas intrusões graníticas proterozóicas (1,58 e 1,87 Ga) ou pelas reativações tectônicas associadas aos Sistemas Transcorrentes Carajás e Cinzento, ou total, provocada pelas últimas. As idades-modelo T_DM de 3,12 e 3,33 Ga para as RMV e RIM e os valores de _εNd (t) de -0,89 a -3,26 sugerem contribuição continental de rochas mais antigas e magmas gerados possivelmente em um ambiente de rifte continental ou de margem continental ativa.

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.

Tectonics and paleogeography of the Marajó Basin, northern Brazil

ABSTRACT: The Marajó Basin area presents geologic and geomorphologic features chiefly due to the Mesozoic extension and post-Miocene neotectonics. The extension event with an Early and a Late Cretaceous phases originated four sub-basins that constitutes the Marajó Basin, with a thick continental clastic sequence showing marine influence. NW and NNW normal faults and NE and ENE strike-slip faults controlled the basin geometry. The extension, related to the Equatorial Atlantic opening, propagated into the continent along crustal weakness zones of the Precambrian Tumucumaque, Amapá and Araguaia orogenic belts. The neotectonic event is a strike-slip regime which developed transtensional basins filled in by Upper Tertiary shallow marine (Pirabas Formation) and transitional sequences (Barreiras Group), followed by Quaternary fluvial deposits and transitional sequences derived from the Amazon and Tocantins rivers and the Marajoara estuary. The current landscape has a typical estuarine morphology. The coast morphology presents sea-cliffs on transitional Upper Tertiary sequences, while inwards dominate hills sustained by Mid-Pleistocene lateritic crust, with a flat erosive surface at 70 m. In the eastern Marajó Island several generations of paleochannels associated with fluvial-estuarine sequences are recognized, while a fluvial-marine plain is widespread on its western side.

Paleomagnetismo de rochas vulvânicas do Nordeste do Brasil e a época da abertura do Oceano Atlântico Sul

Na primeira parte deste trabalho foram desenvolvidos estudos de magnetismo de rochas e paleomagnetismo em amostras de rochas vulcânicas do Nordeste brasileiro. As idades das amostras compreende os períodos Jurássico e Cretáceo. Com este objetivo foram amostradas quatro áreas tendo sido estudado um total de 496 amostras em 55 sítios. Para a coleta foi utilizada uma perfuradora portátil que extrai amostras de 2.5 cm de diâmetro. A orientação das amostras foi feita por meio de uma bússola magnética e de um clinômetro. Os espécimes foram submetidos a desmagnetizações por campo magnético alternado e em alguns poucos casos foi empregada a desmagnetização térmica. Atribuindo-se peso unitário a cada sítio foi determinada a direção média da magnetização remanescente característica de cada uma das áreas estudadas. As rochas vulcânicas do período Jurássico, localizadas na borda oeste da Bacia do Maranhão (Porto Franco-Estreito), apresentaram uma direção media em que D= 3.9°, I= -17.9° com α₉₅= 9.3°, k= 17.9, N= 15 e todos os sítios apresentaram polaridade normal. Para esta área foi determinado o polo paleomagnético de coordenadas 85.3°N, 82.5°E (A₉₅= 6.9º) que se localiza próximo a outros polos paleomagnéticos conhecidos para esse período. As rochas da borda leste da Bacia do Maranhão (Teresina-Picos-Floriano) de idade cretácica inferior apresentaram uma direção média de magnetização remanescente característica tal que D= 174.7°, I= +6.0º com α₉₅= 2.8º, k= 122, N= 21 e todos os sítios apresentaram polaridade reversa. O polo paleomagnético associado a elas apresentou por coordenadas 83.6°N, 261.0°E (A₉₅=1.9°) e mostrou concordância com outros polos sul americanos de mesma idade. No Rio Grande do Norte foi estudado um enxame de diques toleíticos também de idade cretácica inferior, cuja direção média da magnetização remanescente característica encontrada foi D= 186.6º, I= +20.6º com α₉₅= 14.0° e k= 12.9, N= 10. Os sítios desta área apresentaram magnetizações com polaridades normal a reversa. O polo paleomagnético obtido se localiza em 80.6°N e 94.8°E com A₉₅= 9.5°. O estudo das rochas vulcânicas da província magnética do Cabo de Santo Agostinho indicou para a região um valor de D= 0.4º, I= -20.6º com α₉₅= 4.8° e k= 114, N= 9 para a magnetização remanescente característica. Todos os sítios apresentaram polaridade normal e o polo paleomagnético determinado apresentou as seguintes coordenadas: 87.6ºN, 135ºE com A₉₅= 4.5º. Foi discutida a eliminação da variação secular das direções obtidas, de forma que cada polo apresentado nesta dissertação é verdadeiramente um polo paleomagnético. A análise dos minerais magnéticos portadores da remanência, efetuada por curvas termomagnéticas ou por difração de Raio-X, indicou na maior parte das ocorrências, a presença de titanomagnetita pobre em titânio. A presença de maguemita e algumas vezes hematita, na maior parte das vezes resultado de intemperismo, não anulou a magnetização termoremanente associada à época de formação da rocha, que foi determinada após a aplicação de técnicas de desmagnetização aos espécimes. Pelas curvas termomagnéticas obteve-se, para a maioria das amostras, uma temperatura de Curie entre 500 e 600ºC. Os casos mais freqüentes indicaram a ocorrência de titanomagnetita exsolvida, em que foram observadas a presença de uma fase próxima à magnetita e outra fase rica em titânio, próxima à ilmenita, resultado de oxidação de alta temperatura. A segunda parte do trabalho diz respeito à determinação da época de abertura do oceano Atlântico Sul por meio de dados paleomagnéticos. Entretanto ao invés de se utilizar o procedimento comumente encontrado na literatura, e que se baseia nas curvas de deriva polar aparente de cada continente, foi aplicado um teste estatístico que avalia a probabilidade de determinada posição relativa entre os continentes ser válida ou não, para determinado período em estudo. Assim foi aplicado um teste F a polos paleomagnéticos da África e da América do Sul, dos períodos Triássico, Jurássico, Cretáceo Inferior e Cretáceo Médio-Superior, tendo sido estudadas situações que reconstituem a posição pré-deriva dos continentes e configurações que simulem um afastamento entre eles. Os resultados dos testes estatísticos indicaram, dentro de uma probabilidade de erro de menos de 5%, que a configuração pré-deriva de Martin et al (1981) é compatível com os dados paleomagnéticos do Triássico, mas apresenta uma diferença significativa para os paleopolos de Jurássico, Cretáceo Inferior, Cretáceo Médio-Superior. Outras reconstruções pré-deriva testadas apresentaram o mesmo resultado. A comparação entre os polos paleomagnéticos da América do Sul e da África, segundo uma reconstrução que admite uma pequena abertura entre os continentes, como a proposta por Sclater et al (1977) para 110 m.a. atrás, indicou que os dados do Triássico não são compatíveis com este afastamento. Por outro lado os paleopolos do Jurássico e do Cretáceo Inferior, embora mais antigos que a data sugerida pela reconstrução, são consistentes com esta separação dentro de uma probabilidade de erro de menos de 5%. Os dados do Cretáceo Médio-Superior se mostraram consistentes com a reconstrução sugerida para 80 m.a. atrás por Francheteau (1973) e que propõe uma separação maior entre os continentes. Com base na premissa de deslocamentos de blocos continentais rígidos a análise dos resultados obtidos indicou que América do Sul e África estavam unidas por suas margens continentais opostas no período Triássico e que uma pequena separação entre estes continentes, provavelmente devida a uma rutura inicial, ocorreu no Jurássico e se manteve, então, aproximadamente estacionária até o início do Cretáceo Inferior. Esta conclusão difere da maior parte dos trabalhos que discutem a abertura do oceano Atlântico Sul. Os dados do Cretáceo Médio-Superior são compatíveis com um afastamento rápido e significativo entre os continentes naquele período.

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.

Efeitos do manto de intemperismo sobre anomalias eletromagnéticas provocadas por corpos tabulares inclinados: um estudo através de modelamento reduzido

Electromagnetic methods have been extensively applied in the prospecting of sulphide bodies and other conducting materials. The interpretation of e. m. data is based on the results obtained either with reduced scale or analytical modelling. In most models, the host rocks, the overburden as well as the halo of disseminated sulphides are considered highly resistive although in nature they are often weakly conductive. Presence of a well-developed conducting overburden in tropical and sub-tropical regions and a saline rich crustal layer in semi-arid regions have been found to modify significantly the e. m. anomalies in practice. Therefore, the parameters of the target, determined on the basis of simple models where the presence of the conducting environment is neglected, are found to be in considerable error. The effects of the overburden on the e. m. anomalies of a dipping tabular ore body were studied with reduced scale models for varying response parameters of both the overburden and the orebody, and also for different depths and dips of the target. The overburden and the orebody were represented by metallic sheets of varying thickness in the scale model developed in accordance with the law of electromagnetic similitude. The results of these investigations show that the overburden affects the anomaly by causing: a) phase rotation; b) amplitude reduction; c) base level displacement; d) reversion of the quadrature; and e) appearance of an extra peak in the quadrature in the case of low dipping models. The last two effects complicate the quantification of the anomalies but, on the other hand, provide a qualitative indication of the response parameters of both the ore body and the overburden. The results were assembled in the form of Argand diagrams and, finally, an interpretation scheme is suggested for the e. m. field data on the basis of these diagrams.

The ceramic artifacts in archaeological black earth (terra preta) from Lower Amazon Region, Brazil: chemistry and geochemical evolution

Neste trabalho realizou-se a caracterização química de fragmentos de artefatos cerâmicos encontrados em sítios arqueológicos com terra preta no Baixo Amazonas (Cachoeira-Porteira, Pará, Brasil), representativos da cultura Konduri (de 900 a 400 anos AP). Esses fragmentos são constituídos de SiO₂, Al₂O₃, Fe₂O₃, Na₂O e P₂O₅, sendo que SiO₂ e Al₂O₃, juntos, perfazem mais de 80 % em peso. Os teores de P₂O,₅ são relativamente elevados (2,37 % em média) sob a forma de (Al,Fe)-fosfatos, incomuns em cerâmicas vermelhas primitivas, mas encontrados em algumas cerâmicas arqueológicas egípcias e romanas. As concentrações dos elementos traços são comparáveis ou mesmo inferiores ao nível crustal, embora a composição total seja próxima a mesma. A composição química (exceto P₂O₅) em conjunto com os dados mineralógicos e texturais indicam material saprolítico derivado de rochas ígneas félsicas ou rochas sedimentares como matéria-prima das cerâmicas. Os teores de K, Ca e Na mostram que os feldspatos e fragmentos de rochas foram adicionados ao material argiloso, como sugerido pela mineralogia. Os altos teores de sílica respondem pela presença de cauixi, cariapé e/ou areias quartzosas. Fósforo deve ter sido incorporadoà matriz argilosa da cerâmica, quando do cozimento de alimentos nos vasos cerâmicos, e ainda, em parte, durante a formação do perfil de solo tipo ABE sobre Latossolos Amarelos. A matéria prima e os temperos (cauixi, cariapé, rochas trituradas e fragmentos de vasos cerâmicos descartados) encontram-se disponíveis próximos aos sítios até a atualidade, e, portanto foram a área fonte dos mesmos para a confecção dos artefatos cerâmicos.

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 da estrutura geo-elétrica da Região do Juruá, AM, pelo método magnetotelúrico

Utilizando-se dados magnetotelúricos (MT), foi obtida uma imagem geo-elétrica nítida da região do Juruá, Bacia do Solimões, na forma de seções geo-elétricas. Os dados de campo foram registrados ao longo de três linhas de 15 km, espaçadas de 3.5 km, recobrindo uma área de 100 km². O espaçamento entre as 35 estações é irregular, variando de 400 m a 3500 m. A faixa de freqüências utilizada cobriu de 0.001 Hz até 300 Hz, o que permitiu investigar de 100 m até 60 km de profundidade. Os dados apresentam-se afetados pelo efeito de distorção estática. Para corrigir este efeito foi utilizada a mediana da resistividade do primeiro condutor, correspondente à Formação Solimões. Foi utilizado o invariante do tensor MT para interpretar a estrutura geo-elétrica do Juruá. As seções geo-elétricas foram obtidas a partir do agrupamento dos dados resultantes da transformação de Bostick e da inversão 1D de Occam, para cada estação. Foi identificada uma seqüência de camadas condutivas e resistivas, correspondentes ao pacote sedimentar, uma zona de falhas e o topo do embasamento geo-elétrico, caracterizando a Bacia do Solimões. Abaixo do embasamento geo-elétrico foram também identificados uma zona condutora, seguida por uma camada de baixa condutividade, a profundidades iguais ou superiores a 20 km. Esta camada é interpretada como sendo de composição de gabro, estando associada a processos de acreção vertical, intimamente ligados à estabilização crustal e espessamento da litosfera. Os resultados apresentam uma boa concordância com os perfis de resistividade de poços e dados sísmicos de superfície.

Estudo do método de empilhamento SRC e aplicação em dados sintéticos e reais

O empilhamento sísmico é um método de imageamento utilizado para simulação de seções sísmicas com afastamento nulo (AN) a partir de dados de cobertura múltipla, produzindo um aumento da razão sinal/ruído na seção simulada. Este método é realizado usando o processamento convencional denominado ponto médio comum (PMC), através das chamadas correções sobretempo normal e sobretempo oblíquo (STN/STO). Alternativamente, um novo método de empilhamento sísmico, denominado superfície de reflexão comum (SRC) pode ser usado para simular seções sísmicas em afastamento nulo (AN). O método SRC tem como principais características: 1) a construção de uma aproximação do tempo de trânsito de reflexões primárias na vizinhança de um raio de reflexão normal; 2) esta aproximação de empilhamento tem como parâmetros os ângulos de emergência do raio normal e os raios de curvaturas de duas ondas hipotéticas denominadas onda ponto de incidência normal (onda PIN) e onda normal (onda N). O método pressupõe ainda que a velocidade da onda longitudinal na parte superior do modelo próxima à superfície da Terra seja conhecida. Para determinar adequadamente os parâmetros supracitados, fundamentais para o método de imageamento SRC, são necessárias estratégias de busca que utilizem casos particulares da aproximação de segunda ordem hiperbólica dos tempos de trânsito. As estratégias de busca apresentadas são: SRC pragmático-extendito e SRC global-local. Neste trabalho é apresentada uma revisão extensa dos fundamentos teóricos de método de empilhamento SRC, como também das estratégias de busca dos parâmetros SRC, dos algoritmos de implementação, e da delimitação de aberturas do operador de empilhamento SRC usando conceitos de zonas de Fresnel e o fator de estiramento da correção STN. Como um exemplo de aplicação dos atributos ou parâmetros SRC é apresentado o resultado da correção STN, sem estiramento do sinal sísmico, obtido utilizando as velocidades estimadas dos parâmetros SRC. Também, são apresentados os resultados da aplicação dos três métodos de empilhamento (PMC (STN/STO), SRC pragmático-extendido, SRC global-local) nos dados sintéticos Marmousi e real terrestre da bacia paleozóica Amazônica. Finalmente apresentam-se as discussões sobre os resultados, quanto à resolução da imagem sísmica e custo computacional.

Interpretação gravimétrica-magnética da zona de fratura dupla Bode Verde

Em Janeiro de 1980, foi realizada a Operação Bode Verde I do Projeto CENTRATLAN, que constou da coleta de dados geofísicos e geológicos na área central da Cadeia Mesoatlântica (CMA), a sul da Ilha de Ascensão. Esta operação fez parte de um convênio firmado entre a Marinha Americana e a Marinha Brasileira, cujo objetivo era investigar o Oceano Atlântico Sul, especialmente a área onde se estende a província fisiográfica da Cadeia Mesoatlântica. Com o objetivo de estudar a Zona de Fratura Dupla Bode Verde, e assim conhecer melhor esta feição tão marcante no fundo oceânico e margem continental, foram tratados e interpretados os dados batimétricos, gravimétricos e magnéticos colhidos na operação anteriormente citada. A análise dos dados batimétricos permitiu uma definição para esta zona de fratura, com o seu caráter duplo marcado pela presença de duas calhas contínuas, paralelas entre si, e separadas por um alto que tem o seu próprio "rift-valley". As anomalias gravimétricas encontradas confirmam o caráter duplo, com os dois mínimos observados correspondendo às duas calhas da fratura. Estas calhas separam um bloco crustal de aproximadamente 40 km de largura. Com base nos dados gravimétricos, foram construídas nove seções crustais sobre a zona de fratura, que permitiram observar que a interface crosta-manto sofre um afinamento crustal embaixo das paredes da fratura. Esta interface, que encontra-se normalmente a uma profundidade de 8-9 km, passa a atingir profundidades de 5.5-6.0 km abaixo das paredes da fratura, com a crosta apresentando uma espessura de 2.0 km nestas partes. Este afinamento crustal é causado pela subida de material do manto. A extensiva alteração hidrotermal que ocorre na depressão central da fratura, pode ser a responsável pelo menor afinamento crustal observado nesta parte da zona de fratura. Para a interpretação dos dados magnéticos, foi usada uma escala de tempo de polaridade magnética entre o Cretáceo Inferior e o Cenozóico, e uma taxa de espalhamento oceânico de 2.0 cm/ano. Para a camada de basalto, que é responsável por parte das anomalias magnéticas observadas, foi usada uma espessura média de 0.5 km. Nas paredes da zona de fratura, ocorre uma diminuição na espessura desta camada, havendo uma interrupção da mesma na depressão central da fratura.

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.

Mapeamento de uma interface separando dois meios homogêneos através da reconstituição gravimétrica compacta: aplicação a dados do LEPLAC-IV

O problema do mapeamento de uma interface separando dois meios homogêneos tem sido tradicionalmente resolvido através da continuação para baixo da anomalia gravimétrica observada. Este procedimento requer a utilização de um filtro passa-baixa ou de um parâmetro de amortecimento, que tendem a diminuir a resolução do relevo estimado da interface. Além disso, o uso da continuação para baixo implica supor (desnecessariamente) que a interface é uma superfície harmônica. Desta forma, a utilização da continuação para baixo é restrita a interfaces que apresentam um relevo suave relativamente à sua profundidade média. As restrições impostas pela continuação para baixo do campo gravimétrico são reduzidas na técnica de inversão proposta para este estudo, cujo objetivo é mapear a interface crosta-manto baseada na incorporação de informações a priori sobre o mínimo momento de inércia da fonte causadora, em relação a um eixo de concentração de massa pré-estabelecido e coincidente com o topo da fonte anômala. Além da introdução de informações sobre o mínimo momento, foi incluída também a possibilidade de haver variação lateral e em profundidade de densidade em setores previamente estabelecidos, o que permite estender a aplicação desta técnica a situações geológicas mais complexas que envolvam a reconstituição de corpos justapostos, que sejam de diferentes naturezas, tais como aqueles presentes em zonas de transição crustal. O método de inversão apresentado foi avaliado em testes usando dados sintéticos e os resultados obtidos mostraram a eficiência desta técnica em recuperar corpos alongados horizontalmente, que apresentem como topo ou base, interfaces arbitrárias a serem mapeadas. Foi também constatada a eficácia do método em reconstituir corpos que comportem variação lateral de densidade que pudessem retratar situações em que se deseja investigar a natureza do embasamento em regiões de transição crustal. Com base nos resultados destes testes, a técnica foi então empregada na inversão gravimétrica de cinco perfis do LEPLAC IV que cruzam a região oceânica adjacente à margem continental sudeste brasileira. Os resultados desta aplicação permitiram inferir a natureza dos corpos que constituem o arcabouço estrutural da região onde encontram-se inseridas a Dorsal de São Paulo e a Zona de Fratura de Porto Alegre, assim como mapear a Moho subjacente.