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.

Aspectos fundamentais das medidas e interpretação de registros paleomagnéticos em rochas sedimentares da Formação Longá-Bacia do Parnaíba

Considerando que o número de dados paleomagnéticos para o Paleozóico Superior, principalmente Devoniano da América do Sul, é ainda insuficiente não permitindo que uma Curva de Deriva Polar confiável possa ser construída, foram amostrados 43 níveis estratigráficos da Formação Longá, os quais foram estudados paleomagneticamente com o propósito de contribuir para que uma Curva mais confiável possa ser construída. A amostragem foi feita segundo o método dos blocos ao longo dos perfis Teresina, Barras, Batalha, na rodovia PI-13 e Floriano, Nazaré do Piauí, Oeiras nas rodovias PI-24 e BR-230, no estado do Piauí. Os tratamentos foram iniciados no laboratório de Paleomagnetismo do IAG-USP e complementados no do NCGG-UFPa. Utilizou-se as técnicas de desmagnetização progressiva por campos alternados até 700 e/ou temperaturas até 670-700ºC. A interpretação dos resultados foi feita por meio dos diagramas vetoriais de Zijderveld, pelas curvas J-T/C de variação da intensidade magnética com os campos ou temperaturas, e pelos gráficos de variação na direção do vetor magnetização. Os cálculos de direção média e polos foram feitos segundo a Estatística de Fisher (1953). Foram identificadas 4 direções de magnetização remanente: 1. Uma secundária de origem química (CRM) e polaridade reversa, cujo mineral responsável é a hematita produzida provavelmente por alteração deutérica a partir da magnetita. Esta magnetização (identificada pela letra B) quando datada paleomagneticamente (coordenadas do polo 80ºS, 3°E, A₉₅ = 13.6°) indicou idade correspondente ao intervalo Carbonífero-Permiano. 2. Uma componente isotérmica (IRM) dura de espectro totalmente superposto à magnetização inicial que não foi afetada por nenhum dos tratamentos. Esta magnetização foi denominada D e apresentou direção muito estável em torno do ponto de declinação = 234.23º e inclinação = 41.94º. 3. Um grupo de direções de magnetização de origem viscosa (VRM) moles, identificados pela letra C, cuja direção média é dada por: declinação = 15º e inclinação = -20º. Foram removidas a temperatura entre 300 – 600ºC. 4. Finalmente uma magnetização principal, de polaridade normal, denominada A, provavelmente de origem detrítica (DRM) cujo correspondente polo paleomagnético (de coordenadas: 48ºS, 331.7ºE; A₉₅ = 9.9º) mostrou-se compatível com a idade da Formação (Devoniano Superior). Esta magnetização foi considerada inicial. Os polos paleomagnéticos correspondentes às magnetizações A e B, juntamente com outros da América do Sul, foram rotacionados para a África segundo a configuração pré-deriva de Smith e Hallam (1970) e comparados a polos Africanos e Australianos de mesma idade, mostrando-se coerentes. Suas polaridades também estão em acordo com as escalas magnetoestratigráficas publicadas por Irving e Pullaiah (1976) e Khramov e Rodionov (1981).

Estimativas da condutividade térmica dos minerais e rochas e influência de parâmetros térmicos e petrofísicos na resistividade aparente da formação

O presente estudo realiza estimativas da condutividade térmica dos principais minerais formadores de rochas, bem como estimativas da condutividade média da fase sólida de cinco litologias básicas (arenitos, calcários, dolomitos, anidritas e litologias argilosas). Alguns modelos térmicos foram comparados entre si, possibilitando a verificação daquele mais apropriado para representar o agregado de minerais e fluidos que compõem as rochas. Os resultados obtidos podem ser aplicados a modelamentos térmicos os mais variados. A metodologia empregada baseia-se em um algoritmo de regressão não-linear denominado de Busca Aleatória Controlada. O comportamento do algoritmo é avaliado para dados sintéticos antes de ser usado em dados reais. O modelo usado na regressão para obter a condutividade térmica dos minerais é o modelo geométrico médio. O método de regressão, usado em cada subconjunto litológico, forneceu os seguintes valores para a condutividade térmica média da fase sólida: arenitos 5,9 ± 1,33 W/mK, calcários 3.1 ± 0.12 W/mK, dolomitos 4.7 ± 0.56 W/mK, anidritas 6.3 ± 0.27 W/mK e para litologias argilosas 3.4 ± 0.48 W/mK. Na sequência, são fornecidas as bases para o estudo da difusão do calor em coordenadas cilíndricas, considerando o efeito de invasão do filtrado da lama na formação, através de uma adaptação da simulação de injeção de poços proveniente das teorias relativas à engenharia de reservatório. Com isto, estimam-se os erros relativos sobre a resistividade aparente assumindo como referência a temperatura original da formação. Nesta etapa do trabalho, faz-se uso do método de diferenças finitas para avaliar a distribuição de temperatura poço-formação. A simulação da invasão é realizada, em coordenadas cilíndricas, através da adaptação da equação de Buckley-Leverett em coordenadas cartesianas. Efeitos como o aparecimento do reboco de lama na parede do poço, gravidade e pressão capilar não são levados em consideração. A partir das distribuições de saturação e temperatura, obtém-se a distribuição radial de resistividade, a qual é convolvida com a resposta radial da ferramenta de indução (transmissor-receptor) resultando na resistividade aparente da formação. Admitindo como referência a temperatura original da formação, são obtidos os erros relativos da resistividade aparente. Através da variação de alguns parâmetros, verifica-se que a porosidade e a saturação original da formação podem ser responsáveis por enormes erros na obtenção da resistividade, principalmente se tais "leituras" forem realizadas logo após a perfuração (MWD). A diferença de temperatura entre poço e formação é a principal causadora de tais erros, indicando que em situações onde esta diferença de temperatura seja grande, perfilagens com ferramentas de indução devam ser realizadas de um a dois dias após a perfuração do poço.

Interpretação aeromagnética automática com uso da equação homogênea de Euler e sua aplicação na Bacia do Solimões

ABSTRACT: We present here a methodology for the rapid interpretation of aeromagnetic data in three dimensions. An estimation of the x, y and z coordinates of prismatic elements is obtained through the application of "Euler's Homogeneous equation" to the data. In this application, it is necessary to have only the total magnetic field and its derivatives. These components can be measured or calculated from the total field data. In the use of Euler's Homogeneous equation, the structural index, the coordinates of the corners of the prism and the depth to the top of the prism are unknown vectors. Inversion of the data by classical least-squares methods renders the problem ill-conditioned. However, the inverse problem can be stabilized by the introduction of both a priori information within the parameter vector together with a weighting matrix. The algorithm was tested with synthetic and real data in a low magnetic latitude region and the results were satisfactory. The applicability of the theorem and its ambiguity caused by the lack of information about the direction of total magnetization, inherent in all automatic methods, is also discussed. As an application, an area within the Solimões basin was chosen to test the method. Since 1977, the Solimões basin has become a center of exploration activity, motivated by the first discovery of gas bearing sandstones within the Monte Alegre formation. Since then, seismic investigations and drilling have been carried on in the region. A knowledge of basement structures is of great importance in the location of oil traps and understanding the tectonic history of this region. Through the application of this method a preliminary estimate of the areal distribution and depth of interbasement and sedimentary magnetic sources was obtained.

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.

Migração FD 3D em profundidade usando aproximação de Padé complexa

A implementação convencional do método de migração por diferenças finitas 3D, usa a técnica de splitting inline e crossline para melhorar a eficiência computacional deste algoritmo. Esta abordagem torna o algoritmo eficiente computacionalmente, porém cria anisotropia numérica. Esta anisotropia numérica por sua vez, pode levar a falsos posicionamentos de refletores inclinados, especialmente refletores com grandes ângulos de mergulho. Neste trabalho, como objetivo de evitar o surgimento da anisotropia numérica, implementamos o operador de extrapolação do campo de onda para baixo sem usar a técnica splitting inline e crossline no domínio frequência-espaço via método de diferenças finitas implícito, usando a aproximação de Padé complexa. Comparamos a performance do algoritmo iterativo Bi-gradiente conjugado estabilizado (Bi-CGSTAB) com o multifrontal massively parallel solver (MUMPS) para resolver o sistema linear oriundo do método de migração por diferenças finitas. Verifica-se que usando a expansão de Padé complexa ao invés da expansão de Padé real, o algoritmo iterativo Bi-CGSTAB fica mais eficientes computacionalmente, ou seja, a expansão de Padé complexa atua como um precondicionador para este algoritmo iterativo. Como consequência, o algoritmo iterativo Bi-CGSTAB é bem mais eficiente computacionalmente que o MUMPS para resolver o sistema linear quando usado apenas um termo da expansão de Padé complexa. Para aproximações de grandes ângulos, métodos diretos são necessários. Para validar e avaliar as propriedades desses algoritmos de migração, usamos o modelo de sal SEG/EAGE para calcular a sua resposta ao impulso.

Modelagem de dados MT 2D em multicamadas com anisotropia arbitrária

Neste trabalho apresentamos a modelagem direta bidimensional do método Magnetotelúrico (MT) com anisotropia arbitrária na condutividade elétrica para modelos de multicamadas estratificadas horizontalmente, separando a resposta do campo resultante em duas partes, primária e secundária. A parte primária refere-se ao campo 1D do modelo na ausência de heterogeneidades; a secundária, à contribuição no campo resultante gerada pela anomalia. Esta modelagem foi feita via técnica numérica dos Elementos Finitos (EF). Apresentamos também a modelagem direta do caso MT 1D para o mesmo modelo e anisotropia arbitrária, porém, com solução analítica pelo método da matriz de propagação; tal resultado é usado como fonte na formulação 2D. No estudo a respeito do comportamento dos campos 1D em profundidade na presença da anisotropia, verificamos nas componentes horizontais que suas curvas dentro de cada camada se assemelham a de um campo na presença de fontes, pois experimentam um aumento na amplitude; além de verificar o surgimento da componente normal E_z do campo, sem que haja uma densidade de corrente nesta direção. Com arbitrariedade na anisotropia queremos dizer que seus eixos principais podem assumir qualquer orientação em relação ao eixo de referência adotado para o problema; tal orientação é descrita sob três rotações consecutivas, chamadas de strike, dip e slant, que simulam inclinações bastante presentes na geologia dos ambientes terrestres. Verificamos que a anisotropia provoca alterações no comportamento das pseudo-seções de resistividades aparentes, interferindo na interpretação correta dos dados.

Modelagem direta bidimensional do método magnetotelúrico com o método dos elementos finitos de arestas

Fizemos a modelagem direta 2D do método magnetotelúrico (MT) com o método dos elementos finitos (MEF) de arestas em termos dos campos primários e secundários. Para usarmos modelos de maior complexidade e diminuirmos o custo computacional utilizamos malhas não estruturadas. Nas malhas utilizadas, introduzimos quatro nós em torno de cada estação MT, constituindo um quadrado alinhado nas direções dos eixos cartesianos