A implantação de obras civis e de saneamento na Bacia do Una, em Belém do Pará, e as condicionantes relacionadas às características geológicas e geotécnicas

Apresenta a implantação de obras de civis e de saneamento no Projeto de Recuperação das Baixadas de Belém -Bacia do Una, na cidade de Belém, no estado do Pará, que envolve obras de macrodrenagem, microdrenagem, redes de água, redes de esgoto e sistema viário. Também chamado de Projeto Una ou simplesmente de Projeto de Macrodrenagem, alcança uma área de 3.644 ha , correspondente a, aproximadamente, 21% da área continental e 7,20% da área total do município de Belém, desenvolve-se na Bacia do Una, situada em área de baixada, com grande parte de suas terras abaixo da cota de alagamento, isto é, até 3,60 m acima do nível do mar, apresenta uma camada espessa de argila orgânica, na maioria das vezes, superficial, atingindo em alguns lugares a grandes profundidades, inadequada, portanto para o suporte de fundações, estando as camadas de solos mais resistentes a profundidades variáveis, atingindo, em certas áreas, até 60 metros. Além das dificuldades de construção em solos dessa natureza, esta obra apresenta inúmeros condicionantes, como o remanejamento de famílias de áreas de risco, isto é, famílias que residem sobre o leito de canais, dificultando o escoamento das águas servidas e pluviais; as indenizações de residências abaladas pelas obras; a falta de espaço físico para a implantação das obras; a ruptura dos taludes das marginais de canais; as dificuldades de gerenciamento face a extensa abrangência da área; as relações com a comunidade, tumultuadas em alguns momentos e, principalmente, a priorização de algumas áreas, em detrimento de outras, para a execução das obras, face os recursos disponíveis. Fundamentado em relatórios técnicos, projetos, pesquisa de campo e experiências pessoais do autor, a dissertação faz uma descrição do projeto em seus diversos aspectos, bem como promove revisão critica do mesmo sob a ótica da engenharia civil, envolvendo aspectos de planejamento e execução de obras, com o objetivo de construir e divulgar ensinamentos que promovam a otimização, no desenvolvimento de obras semelhantes.

Influência de estruturas geológicas bidimensionais no campo geoeletromagnético na presença do eletrojato equatorial

A Terra atua como um grande magneto esférico, cujo campo assemelha-se àquele gerado por um dipolo magnético. Este campo apresenta mudanças de intensidade que variam com a localização e a hora local. A parte principal do campo geomagnético se origina no interior da Terra através de processos eletromagnéticos. Extensivos estudos mostraram ainda que existem contribuições de origem externa ao planeta, principalmente de origem solar. Dentre estas fontes há anomalias do campo magnético que surgem a partir de um aumento diurno da corrente elétrica em uma estreita faixa da ionosfera, de direção leste-oeste, centrada no equador magnético e denominada Eletrojato Equatorial (EEJ). Ocasionalmente estas correntes podem apresentar reversões de fluxo, sendo denominadas Contra-Eletrojato (CEJ). Vários autores têm estudado os efeitos do EEJ e CEJ sobre as observações geoeletromagnéticas. Eles estão interessados no efeito combinado do EEJ e estruturas geológicas condutivas 1-D e 2-D. Nestes trabalhos a estrutura 2-D sempre se apresentava paralela ao eletrojato, o que é uma hipótese bastante restritiva ao se modelar ambientes geológicos mais realistas, em que corpos bidimensionais podem ter qualquer strike em relação ao EEJ. Neste trabalho apresentamos a solução deste problema sem esta restrição. Assim, mostramos os campos geoeletromagnéticos devidos a estruturas bidimensionais que possuam strike oblíquo em relação ao EEJ, através de perfis dos campos elétrico e magnéticos calculados na superfície e formando direção arbitrária à heterogeneidade condutiva 2-D. Com esta resposta avaliamos ainda qual a influência que estruturas bidimensionais exercem sobre a resposta magnetotelúrica, sob influência do Eletrojato Equatorial. Durante o desenvolvimento deste trabalho, utilizamos o método de elementos finitos, tendo por fonte eletromagnética o EEJ e o CEJ, que por sua vez foram representados por uma combinação de distribuições gaussianas de densidade de corrente. Estas fontes foram decompostas nas direções paralela e perpendicular à estrutura 2-D, resultando nos modos de propagação TE₁ e TE₂ e TM acoplados, respectivamente. Resolvemos o modo acoplado aplicando uma Transformada de Fourier nas equações de Maxwell e uma Transformada Inversa de Fourier na solução encontrada. De acordo com os experimentos numéricos realizados em um modelo interpretativo da Anomalia Condutiva da Bacia do Parnaíba, formado por uma enorme estrutura de 3000 ohm-m dentro de um corpo externo condutivo (1 ohm-m), concluímos que a presença do CEJ causa uma inversão na anomalia, se compararmos com o resultado do EEJ. Concluímos também que para as frequências mais altas as componentes do campo elétrico apresentam menor influência da parte interna do corpo 2-D do que da parte externa. Já para frequências mais baixas este comportamento se observa com as componentes do campo magnético. Com relação à frequência, vimos os efeitos do “skin-depth”, principalmente nas respostas magnéticas. Além disso, quando a estrutura 2-D está paralela ao eletrojato, o campo elétrico é insensível à estrutura interna do modelo para todos os valores de frequência utilizados. Com respeito ao ângulo θ_h entre a heterogeneidade e a fonte, vimos que o modo TM se manifesta naturalmente quando θ_h é diferente de 0°. Neste caso, o modo TE é composto por uma parte devido à componente da fonte paralela à heterogeneidade e a outra devido à componente da fonte perpendicular, que é acoplada ao modo TM. Assim, os campos calculados têm relação direta com o valor de θ_h. Analisando a influência do ângulo entre a direção do perfil dos campos e o strike da heterogeneidade verificamos que, à medida que θ_h se aproxima de 90°, os campos primários tornam-se variáveis para valores de θ_p diferentes de 90°. Estas variações causam uma assimetria na anomalia e dão uma idéia da inclinação da direção do perfil em relação aos corpos. Finalmente, concluímos que uma das influências que a distância entre o centro do EEJ e o centro da estrutura 2-D, causa sobre as componentes dos campos está relacionado às correntes reversas do EEJ e CEJ, pois a 500 km do centro da fonte estas correntes têm máxima intensidade. No entanto, com o aumento da distância, as anomalias diminuem de intensidade. Nas sondagens MT, nós também usamos o EEJ e o CEJ como fonte primária e comparamos nossos resultados com a resposta da onda plana. Deste modo observamos que as componentes do campo geoeletromagnético, usadas para calcular a impedância, têm influência do fator de acoplamento entre os modos TE₂ e TM. Além disso, esta influência se torna maior em meios resistivos e nas frequências mais baixas. No entanto, o fator de acoplamento não afeta os dados magnetotelúricos em frequências maiores de 10^-2 Hz. Para frequências da ordem de 10^-4 Hz os dados MT apresentam duas fontes de perturbação: a primeira e mais evidente é devido à presença fonte 2-D (EEJ e CEJ), que viola a hipótese da onda plana no método MT; e a segunda é causada pelo acoplamento entre os modos TE₂ e TM, pois quando a estrutura bidimensional está obliqua à fonte 2-D temos correntes elétricas adicionais ao longo da heterogeneidade. Concluimos assim, que o strike de uma grande estrutura condutiva bidimensional relativamente à direção do EEJ ou CEJ tem de fato influência sobre o campo geomagnético. Por outro lado, para estudos magnetotelúricos rasos (frequências maiores de 10^-3 Hz) o efeito do ângulo entre a estrutura geológica 2-D e a direção do EEJ não é tão importante. Contudo, em estudos de litosfera frequências menores de 10^-3 Hz) o acoplamento entre os modos TE₂ e TM não pode ser ignorado.

Análise de feições geológicas na região de Tartarugalzinho, Amapá, através de métodos geofísicos

Dois métodos geofísicos foram aplicados para a detecção de estruturas geológicas na borda oriental do Escudo das Guianas. O método gravimétrico foi aplicado para encontrar falhas, discordâncias e descontinuidades no embasamento, associadas a cinturões de cisalhamento que possuem orientação NW-SE, e o posicionamento da profundidade com que o embasamento se encontra da superfície. Foram encontradas falhas no embasamento e também indicações de que a crosta apresenta um afinamento em certa região da borda. O método eletromagnético GPR foi utilizado para encontrar descontinuidades geológicas nas rochas sedimentares e averiguar se essas feições são o reflexo das falhas no embasamento e/ou reativação de algumas delas. O GPR localizou feições interpretadas como descontinuidades em rochas e paleoestruturas. O estudo teve o intuito de compreender como as descontinuidades do embasamento podem ter influenciado na formação de bacias sedimentares do Terciário e na formação de paleoestruturas rasas.

Efeitos do strike das estruturas geológicas de duas dimensões nas pseudos-seções de IP resistividade

Os levantamentos de IP-resistividade efetuados na Serra dos Carajás não foram executados ortogonalmente às estruturas geológicas, pois utilizaram linhas anteriormente abertas pelas equipes de geoquímica. Este fato motivou este estudo teórico da influência da direção das linhas de medidas de IP-resistividade em relação ao "strike" da estrutura. Foi usado o programa de elementos finitos de Rijo (1977), desenvolvido para levantamentos perpendiculares às estruturas com as adaptações necessárias. A modificação principal foi na rotina de transformação inversa de Fourier. Para o caso simples dos levantamentos perpendiculares, a transformada inversa é uma integral discreta com apenas sete pontos. No entanto, para as medidas obliquas, o integrando é oscilatório, e portanto, a integral a ser calculada é mais complexa. Foi adaptado um método apresentado por Ting e Luke (1981), usando dezoito pontos em cada integração. Foi constatado que o efeito da direção da linha em relação ao "strike" é desprezível para ângulos maiores que 60 graus. Para ângulos menores, o efeito consiste no alongamento da anomalia, com pequenas alterações em seu centro. Não há uma maneira simples de compensar este efeito com mudanças nos parâmetros do modelo.

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.

Hidroquímica do rio Solimões na região entre Manacapuru e Alvarães: Amazonas - Brasil

Este trabalho discute as características físico-químicas das águas dos rios Solimões, Purus e seus afluentes, coletadas em novembro de 2004 no Estado do Amazonas, entre as cidades de Manacapuru-Alvarães e Anamã-Pirarauara. Foram realizadas análises físico-químicas (temperatura, pH, condutividade elétrica, turbidez, Ca2+, Na+, K+, Mg2+, HCO3-, SO42-, Cl-), de elementos-traço (Li, B, Al, Sc, V, Cr, Mn, Fe, Co, Cu, Zn, As, Se, Rb, Sr, Mo, Cd, Sb, Cs, Ba, Pb, La, Ce e U) e isótopos de estrôncio. Os parâmetros analisados e a composição química mostram que as águas dos rios e igarapés da região central da Amazônia são quimicamente distintas entre si. As águas brancas do Solimões são cálcicas-bicarbonatadas e as do Purus bicarbonatadas, os respectivos afluentes são sódico-potássico-bicarbonatados e sódico-potássico-sulfatados. Isso acarreta águas brancas fracamente ácidas a neutras e mais condutivas, enquanto as pretas são menos mineralizadas, mais ácidas, especialmente as do Purus. O Ba, Sr, Cu, V e As mais elevados diferenciam as águas brancas do Solimões das do Purus, bem como os afluentes do primeiro em relação ao segundo. Esse conjunto de características indicam que tanto o Solimões, como o Purus e os respectivos afluentes, estão submetidos a condições geológicas/ambientais distintas. A influência do aporte de sedimentos dos Andes é diluída ao longo da bacia do Solimões e se reflete na formação das várzeas dos Solimões e Purus. Por outro lado as rochas crustais, representadas pelos escudos das Guianas e Brasileiro também contribuem, mas em menor proporção.

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.

Estoque de biomassa e carbono florestal em unidades de paisagem na Amazônia: uma análise a partir da abordagem metodológica ecologia da paisagem

O presente trabalho tem por objetivo central demonstrar a variabilidade existente na floresta no que tange aos estoques de biomassa e carbono florestal acima do solo, a partir da identificação e caracterização, com base em técnicas de sensoriamento remoto, de unidades de paisagem em uma área situada no município de Belterra, região oeste do Estado do Pará, a partir da matriz teórico-conceitual da abordagem Ecologia da Paisagem. Para o alcance de tal proposição, a metodologia empregada partiu da revisão da literatura sobre o tema, aquisição de dados cartográficos e orbitais, uso de técnicas de sensoriamento remoto, coleta de dados em campo, tratamento e análise estatística. O trabalho está dividido em quatro capítulos, seguidos pelas considerações gerais da obra. Partindo da matriz teórico-metodológica da Ecologia da Paisagem, analisa-se a dinâmica socioambiental do município de Belterra, que atualmente experimenta a expansão das atividades agrícolas, com destaque para a agricultura mecanizada da soja. A partir da análise multitemporal de imagens Landsat do município pôde-se avaliar a distribuição da cobertura florestal existente no mesmo, bem como o padrão espacial de distribuição das principais unidades de paisagem identificadas. Considerando esse recorte, realizou-se a coleta de dados em campo via inventário florestal em quatro tipologias florestais (floresta de alto platô, floresta de baixo platô, vegetação secundária e tensão ecológica) para obtenção de parâmetros morfométricos da vegetação e posterior quantificação dos estoques de biomassa e carbono contidos em cada unidade, bem como observar o comportamento estrutural da floresta nas mesmas. A adoção da paisagem como escala espacial de análise mostrou-se bastante satisfatória na quantificação dos estoques de biomassa e carbono florestal ao permitir considerar a influência da dinâmica socioeconômica na redução desses estoques. Além disso, possibilitou constatar que o reconhecimento da heterogeneidade da cobertura florestal é um elemento fundamental para a obtenção de estimativas de carbono de acordo com as características estruturais da vegetação, que varia de acordo com a topografia do terreno, com as espécies existentes e com as características geográficas, o que envolve a tipologia climática, as características geomorfológicas, pedológicas e geológicas da área.

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).

Determinação das velocidades intervalares usando a teoria paraxial do raio: aproximação de segunda ordem dos tempos de trânsito

Neste trabalho foi desenvolvido um método de solução ao problema inverso para modelos sísmicos compostos por camadas homogêneas e isotrópicas separadas por superfícies suaves, que determina as velocidades intervalares em profundidade e calcula a geometria das interfaces. O tempo de trânsito é expresso como uma função de parâmetros referidos a um sistema de coordenadas fixo no raio central, que é determinada numericamente na superfície superior do modelo. Essa função é posteriormente calculada na interface anterior que limita a camada não conhecida, através de um processo que determina a função característica em profundidade. A partir da função avaliada na interface anterior se calculam sua velocidade intervalar e a geometria da superfície posterior onde tem lugar a reflexão do raio. O procedimento se repete de uma forma recursiva nas camadas mais profundas obtendo assim a solução completa do modelo, não precisando em nenhum passo informação diferente à das camadas superiores. O método foi expresso num algoritmo e se desenvolveram programas de computador, os quais foram testados com dados sintéticos de modelos que representam feições estruturais comuns nas seções geológicas, fornecendo as velocidades em profundidade e permitindo a reconstrução das interfaces. Uma análise de sensibilidade sobre os programas mostrou que a determinação da função característica e a estimação das velocidades intervalares e geometria das interfaces são feitos por métodos considerados estáveis. O intervalo empírico de aplicabilidade das correções dinâmicas hiperbólicas foi tomado como uma estimativa da ordem de magnitude do intervalo válido para a aplicação do método.

Imageamento da porosidade através de perfis geofísicos de poço

O imageamento da porosidade é uma representação gráfica da distribuição lateral da porosidade da rocha, estimada a partir de dados de perfis geofísicos de poço. Apresenta-se aqui uma metodologia para produzir esta imagem geológica, totalmente independente da intervenção do intérprete, através de um algoritmo, dito, interpretativo baseado em dois tipos de redes neurais artificiais. A primeira parte do algoritmo baseia-se em uma rede neural com camada competitiva e é construído para realizar uma interpretação automática do clássico gráfico o P_b - Φ_N, produzindo um zoneamento do perfil e a estimativa da porosidade. A segunda parte baseia-se em uma rede neural com função de base radial, projetado para realizar uma integração espacial dos dados, a qual pode ser dividida em duas etapas. A primeira etapa refere-se à correlação de perfis de poço e a segunda à produção de uma estimativa da distribuição lateral da porosidade. Esta metodologia ajudará o intérprete na definição do modelo geológico do reservatório e, talvez o mais importante, o ajudará a desenvolver de um modo mais eficiente as estratégias para o desenvolvimento dos campos de óleo e gás. Os resultados ou as imagens da porosidade são bastante similares às seções geológicas convencionais, especialmente em um ambiente deposicional simples dominado por clásticos, onde um mapa de cores, escalonado em unidades de porosidade aparente para as argilas e efetiva para os arenitos, mostra a variação da porosidade e a disposição geométrica das camadas geológicas ao longo da seção. Esta metodologia é aplicada em dados reais da Formação Lagunillas, na Bacia do Lago Maracaibo, Venezuela.

Modelagem numérica da influência do eletrojato equatorial em dados magnetotelúricos produzidos por estruturas tridimensionais

A América do Sul apresenta várias peculiaridades geomagnéticas, uma delas, é a presença do Eletrojato Equatorial, o qual se estende de leste para oeste no Brasil ao longo de aproximadamente 3500 km. Considerando-se o fato de que a influência do Eletrojato Equatorial pode ser detectada a grandes distâncias do seu centro, isto suscita o interesse em se estudar os seus efeitos na exploração magnetotelúrica no Brasil. A influência do eletrojato equatorial na prospecção magnetotelúrica tem sido modelada para meios geológicos uni e bidimensionais valendo-se para isto de soluções analíticas fechadas e de técnicas numéricas tais como elementos finitos e diferenças finitas. Em relação aos meios geológicos tridimensionais, eles tem sido modelados na forma de "camadas finas", usando o algoritmo "thin sheet". As fontes indutoras utilizadas para simular o eletrojato equatorial nestes trabalhos, tem sido linhas de corrente, eletrojatos gaussianos e eletrojatos ondulantes. Por outro lado, o objetivo principal da nossa tese foi o modelamento dos efeitos que o eletrojato equatorial provoca em estruturas tridimensionais próprias da geofísica da prospecção. Com tal finalidade, utilizamos o esquema numérico da equação integral, com as fontes indutoras antes mencionadas. De maneira similar aos trabalhos anteriores, os nossos resultados mostram que a influência do eletrojato equatorial somente acontece em frequências menores que 10^-1 Hz. Este efeito decresce com a distância, mantendo-se até uns 3000 km do centro do eletrojato. Assim sendo, a presença de grandes picos nos perfis da resistividade aparente de um semi-espaço homogêneo, indica que a influência do eletrojato é notável neste tipo de meio. Estes picos se mostram com diferente magnitude para cada eletrojato simulado, sendo que a sua localização também muda de um eletrojato para outro. Entretanto, quando se utilizam modelos geo-elétricos unidimensionais mais de acordo com a realidade, tais como os meios estratificados, percebe-se que a resposta dos eletrojatos se amortece significativamente e não mostra muitas diferenças entre os diferentes tipos de eletrojato. Isto acontece por causa da dissipação da energia eletromagnética devido à presença da estratificação e de camadas condutivas. Dentro do intervalo de 3000 km, a resposta eletromagnética tridimensional pode ser deslocada para cima ou para baixo da resposta da onda plana, dependendo da localização do corpo, da frequência, do tipo de eletrojato e do meio geológico. Quando a resposta aparece deslocada para cima, existe um afastamento entre as sondagens uni e tridimensionais devidas ao eletrojato, assim como um alargamento da anomalia dos perfis que registra a presença da heterogeneidade tridimensional. Quando a resposta aparece deslocada para baixo, no entanto, há uma aproximação entre estes dois tipos de sondagens e um estreitamento da anomalia dos perfis. Por outro lado, a fase se mostra geralmente, de uma forma invertida em relação à resistividade aparente. Isto significa que quando uma sobe a outra desce, e vice-versa. Da mesma forma, comumente nas altas frequências as respostas uni e tridimensionais aparecem deslocadas, enquanto que nas baixas frequências se mostram com os mesmos valores, com exceção dos eletrojatos ondulantes com parâmetros de ondulação α = —2 e —3. Nossos resultados também mostram que características geométricas próprias das estruturas tridimensionais, tais como sua orientação em relação à direção do eletrojato e a dimensão da sua direção principal, afetam a resposta devido ao eletrojato em comparação com os resultados da onda plana. Desta forma, quando a estrutura tridimensional é rotacionada de 90°, em relação à direção do eletrojato e em torno do eixo z, existe uma troca de polarizações nas resistividades dos resultados, mas não existem mudanças nos valores da resistividade aparente no centro da estrutura. Ao redor da mesma, porém, se percebe facilmente alterações nos contornos dos mapas de resistividade aparente, ao serem comparadas com os mapas da estrutura na sua posição original. Isto se deve à persistência dos efeitos galvânicos no centro da estrutura e à presença de efeitos indutivos ao redor do corpo tridimensional. Ao alongar a direção principal da estrutura tridimensional, as sondagens magnetotelúricas vão se aproximando das sondagens das estruturas bidimensionais, principalmente na polarização XY. Mesmo assim, as respostas dos modelos testados estão muito longe de se considerar próximas das respostas de estruturas quase-bidimensionais. Porém, os efeitos do eletrojato em estruturas com direção principal alongada, são muito parecidos com aqueles presentes nas estruturas menores, considerando-se as diferenças entre as sondagens de ambos tipos de estruturas. Por outro lado, os mapas de resistividade aparente deste tipo de estrutura alongada, revelam um grande aumento nos extremos da estrutura, tanto para a onda plana como para o eletrojato. Este efeito é causado pelo acanalamento das correntes ao longo da direção principal da estrutura. O modelamento de estruturas geológicas da Bacia de Marajó confirma que os efeitos do eletrojato podem ser detetados em estruturas pequenas do tipo "horst" ou "graben", a grandes distâncias do centro do mesmo. Assim, os efeitos do eletrojato podem ser percebidos tanto nos meios estratificados como tridimensionais, em duas faixas de freqüência (nas proximidades de 10^-1 Hz e para freqüências menores que 10^-3 Hz), possivelmente influenciados pela presença do embasamento cristalino e a crosta inferior, respectivamente. Desta maneira, os resultados utilizando o eletrojato como fonte indutora, mostram que nas baixas freqüências as sondagens magnetotelúricas podem ser fortemente distorcidas, tanto pelos efeitos galvânicos da estrutura tridimensional como pela presença da influência do eletrojato. Conseqüêntemente, interpretações errôneas dos dados de campo podem ser cometidas, se não se corrigirem os efeitos do eletrojato equatorial ou, da mesma forma, não se utilisarem algoritmos tridimensionais para interpretar os dados, no lugar do usual modelo unidimensional de Tikhonov - Cagniard.

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.

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.

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.