FOG: função de objetivo genérico para roteamento e alocação de comprimentos de onda em redes ópticas com tráfego de demanda estática

Esta tese de doutorado propõe uma formulação matemática para simulação de roteamento e alocação de comprimentos de onda em redes ópticas, sem a inclusão de restrições que não são inerentes ao problema básico e com o objetivo de ser aplicável a qualquer tipo de rede óptica com tráfego de demanda estática. O estabelecimento de uma rota seguida da seleção de um comprimento de onda é um dos pontos chave para o bom funcionamento de uma rede óptica, pois influencia na forma como os recursos da rede serão gerenciados. Assim, o processo de roteamento e alocação de comprimentos de onda em redes ópticas, conhecido como RWA (Routing and Wavelength Assignment), necessita de soluções que busquem a sua otimização. Entretanto, a despeito dos inúmeros estudos com o objetivo de otimizar o processo RWA, observa-se que não há, a priori, nenhuma solução que possa levar a uma padronização do referido processo. Considerando que a padronização é desejável na consolidação do uso de qualquer tecnologia, a Tese descrita neste trabalho é uma Função de Objetivo Genérico (FOG) que trata do processo de roteamento e alocação de comprimentos de onda, visando estabelecer uma base a partir da qual seja possível desenvolver um padrão ou vários padrões para redes ópticas. A FOG foi testada, via simulação, no processo de alocação de comprimentos de onda do inglês, Wavelength Assignment e no processo RWA como um todo. Em ambos os casos, os testes foram realizados considerando redes opacas, trazendo resultados surpreendentes, considerando a simplicidade da solução para um problema não trivial.

A Criança na fotografia: o retrato da infância na primeira metade do século XX em Belém do Pará (1900 a 1950)

O presente estudo tem por objetivo analisar a infância na fotografia de crianças na cidade de Belém do Pará da primeira metade do século XX. Portanto, a pesquisa versa sobre a importância da fotografia enquanto documento histórico, que revela múltiplos significados; importando, sobretudo, o seu potencial enquanto linguagem representativa da infância e ou de uma concepção de infância construída no cenário da cidade de Belém do Pará, no período que compreende a primeira metade do século XX. A fotografia é uma fonte histórica, é uma imagem/documento. No caso de fotos para álbum de família, tem-se por registro da história familiar e sócio familiar, não necessariamente devem estar dentro de um álbum, mas apenas, que apresentem tais características, pois literalmente revelam os modos e as circunstancias em que a família ou parte dela é fotografada. Já as fotografias de crianças em revistas, são deslocamentos das fotos que geralmente eram produzidas para interesse da própria família para colocá-las em seus álbuns, entretanto, também eram utilizadas como forma de promoção familiar dentro de publicações avulsas nas revistas ou periódicos das primeiras décadas do século XX. Geralmente acompanhados de uma mensagem discursiva elevando o “nome” da família. A fotografia permite que quase toda pessoa (não só as mais abastadas) possa se transformar em objeto-imagem, ou numa série de imagens que retratam momentos de suas vidas. No caso, as imagens fotográficas encontradas nos álbuns dos festejos escolares são a prova disso: fotografias exclusivamente de crianças e para com fins comemorativos de alguns eventos anuais presentes nos lazeres infantis educacionais. A fotografia poderia ser tomada como um equivalente da memória individual e coletiva, com a imagem fixada de um tempo que, aparentemente, foi “recortado”. A fotografia de crianças em periódicos e álbum de família é um caminho para se compreender e conhecer a história social, cultural e educacional da infância na Amazônia, pois as imagens fotográficas podem revelar muitos índices a respeito de uma família e, por extensão, da sociedade: seus ritos, modos de vida, afetividade e ideais. O corpus da pesquisa é composto de fotografias de crianças arquivadas na Biblioteca Artur Viana. Para aprofundamento teórico sobre fotografia e imagem utilizaremos Dubois, Kossoy, Benjamin, Fabris, Leite, Mauad, Burke e Fischman. Quanto à leitura sobre infância destacamos Rizzini, Del Priori, Freitas, e Kulmman. No contexto da realidade da Amazônia privilegiamos os seguintes autores: Bezzera Neto e Sarges.

Caracterização morfológica de Liophis reginae semilineatus (Wagler, 1824) e Liophis reginae macrosomus (Amaral, 1935), e o status taxonômico de Liophis oligolepis Boulenger, 1905

Foi analisado um total de 249 espécimes das subespécies de Liophis reginae que ocorrem no Brasil (Liophis reginae macrosomus e Liophis reginae semilineatus), com o objetivo de caracterizar morfologicamente os dois táxons e verificar a validade de Liophis oligolepis, táxon considerado sinônimo de Liophis reginae semilineatus. Foram considerados 22 dados merísticos e 18 morfométricos. Os complexos osteológicos e cefálicos, assim como os hemipênis foram analisados comparativamente entre os táxons. Foram identificadas, utilizando MANOVA, diferenças sexuais significativas em variáveis mensuradas da cabeça e do corpo. Uma análise da função discriminante (AFD) foi utilizada em sexos separados para maximizar a separação, num espaço rnultivariado, dos três táxons definidos a priori: Liophis regime semilineatus, Liophis reginae macrosomus e Liophis oligolepis. Os escores dos espécimes machos separaram claramente Liophis reginae semilineatus de Liophis oligolepis no eixo da primeira função discriminante e Liophis reginae semilineatus de Liophis reginae macrosomus na segunda função. Em fêmeas, a primeira função discriminante separou Liophis oligolepis de Liophis reginae semilineatus e Liophis reginae macrosomus. Na AFD de Liophis oligolepis e Liophis reginae semilineatus, considerando os dois sexos juntos, houve uma diferenciação entre os dois táxons na primeira função discriminante. A mesma análise foi feita para Liophis reginae semilineatus e Liophis reginae rnacrosomus, sendo observado uma separação destes táxons na primeira função discriminante. Verificou-se alguns exemplares de Liophis reginae semilineatus e Liophis reginae macrosomus além das áreas de distribuição registradas para essas subespécies. Para verificar se estes e os outros exemplares, ocorrentes nas áreas de distribuição citadas em bibliografia, formavam o mesmo grupo, foi feita uma análise da função discriminante. Os resultados indicam que as novas ocorrências pertencem ao mesmo grupo de espécimes previamente estabelecidos. A morfologia hemipeniana não diferiu entre Liophis reginae semilineatus, Liophis oligolepis e Liophis reginae macrosomus. Os três táxons apresentam crânios semelhantes no aspecto geral, porém existem diferenças, no osso parietal de Liophis reginae semilineatus e Liophis oligolepis. O status taxonômico de Liophis oligolepis foi definido na categoria específica com base nos caracteres meristicos, morfométricos e morfológicos. As subespécies Liophis reginae semilineatus e Liophis reginae macrosomus deverão ser reavaliadas no futuro, para a verificação de seus status. Com a análise de mais exemplares e a inclusão de Liophis reginae regime e Liophis reginae zweiftli nas análises.

Quando a paisagem torna-se obra: uma abordagem ecológica das composições do tipo paisagem sonora

Pós-graduação em Música - IA

Deconvolução de processo sísmico não-estacionário

O presente trabalho trata da aplicação do filtro Kalman-Bucy (FKB), organizado como uma deconvolução (FKBD), para extração da função refletividade a partir de dados sísmicos. Isto significa que o processo é descrito como estocástico não-estacionário, e corresponde a uma generalização da teoria de Wiener-Kolmogorov. A descrição matemática do FKB conserva a relação com a do filtro Wiener-Hopf (FWH) que trata da contra-parte com um processo estocástico estacionário. A estratégia de ataque ao problema é estruturada em partes: (a) Critério de otimização; (b) Conhecimento a priori; (c) Algoritmo; e (d) Qualidade. O conhecimento a priori inclui o modelo convolucional, e estabelece estatísticas para as suas componentes do modelo (pulso-fonte efetivo, função refletividade, ruídos geológico e local). Para demostrar a versatilidade, a aplicabilidade e limitações do método, elaboramos experimentos sistemáticos de deconvolução sob várias situações de nível de ruídos aditivos e de pulso-fonte efetivo. Demonstramos, em primeiro lugar, a necessidade de filtros equalizadores e, em segundo lugar, que o fator de coerência espectral é uma boa medida numérica da qualidade do processo. Justificamos também o presente estudo para a aplicação em dados reais, como exemplificado.

Contribuição da dinâmica morfoestrutural e morfoescultural na estruturação da bacia do Ribeirão Entupido, complexo alcalino do Passa Quatro - SP

Pós-graduação em Geografia - IGCE

Interrogando práticas do UNICEF para os adolescentes no Brasil

Este trabalho investigou e interrogou as práticas discursivas do UNICEF direcionadas aos “adolescentes” brasileiros. Utilizou-se o método histórico-genealógico foucaultiano para interrogar o relatório “Situação da Adolescência Brasileira” (2002), que se constituiu como fonte privilegiada desta pesquisa. Desse modo, os questionamentos que moveram o estudo foram: que práticas do UNICEF incidem sobre os corpos de adolescentes brasileiros, no século XX e início do século XXI? Que subjetividades essas práticas produzem? Como objetivam a adolescência? Que relações de poder acionam frente a esses corpos? Que efeitos elas produzem? Tais problematizações não tiveram por finalidade, fazer a história do falso ou do verdadeiro, pois isso não tem importância política, mas problematizar a produção dos regimes de verdades a respeito destes sujeitos e os efeitos destes na atualidade. Dessa forma, marcar a singularidade dos acontecimentos que forjaram este objeto como um problema para as ciências humanas, e como uma questão para o UNICEF e para o Sistema de Garantia de Direitos. O objetivo do estudo foi analisar as práticas discursivas de poder e subjetivação que objetivam e subjetivam a adolescência brasileira. De posse da ferramenta foucaultiana, desmontamos o documento, cortamos as séries que o compõem, desarticulamos as pretensas continuidades, reescrevemos e reinventamos o objeto adolescência, deixando em suspenso as certezas e verdades que o atravessam e que pretendem constituí-lo como objeto natural, imersos em essencialismos e homogeneizações. Como resultados, identificamos dicotomias no documento, como: potencialidade/risco, fase positiva/negativa, por exemplo, que tentam naturalizar o sujeito como algo dado a priori, portador de uma essência objetivado e subjetivado por uma perspectiva linear do desenvolvimento humano, como: adaptação/desadaptação, normal/anormal, maturidade/imaturidade e uma sequência linear de fases, que atende também a concepções econômicas desenvolvimentistas e neoliberais preocupadas com a equação custo-benefício.Foi com um olhar atento às ninharias do poder, que buscamos destruir certezas e evidências, atentando não para as intencionalidades dos jogos de forças, mas, ao acaso das lutas.

Estados emocionais de estresse e cortisol salivar na escalada em rocha

Pós-graduação em Ciências da Motricidade - IBRC

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.

Geotermia rasa em Belém

Um estudo detalhado de geotermia rasa foi realizado no intervalo de profundidade de 0,02 a 210,0 m, na região metropolitana de Belém. As medidas de temperatura foram efetuadas com o uso de termômetros de mercúrio e termistor, enquanto as medidas de condutividade térmica foram feitas em testemunhos de sondagem, utilizando-se o aparato tipo agulha. O fluxo de radiação solar incidente foi registrado com um actinógrafo do tipo Robitzech. O gradiente geotérmico médio obtido na região metropolitana de Belém foi de 0,0254 ± 0,0007 °C.m^-1. O valor médio de condutividade térmica dos testemunhos de sondagem coletados nesta região foi de 1,66 ± 0,52 W.m^-1.°C^-1. Por sua vez, o fluxo geotermal médio na região estudada é de 42,16 ± 1,14 mW.m^-2. A 1,0 m de profundidade ocorrem mudanças temporais nos valores de temperatura, as quais não podem ser consideradas como desprezíveis. Estas, mudanças estão diretamente relacionadas com as variações do fluxo de radiação solar incidente na superfície. O maior incremento deste fluxo, observado de um dia para o outro, foi de cerca de 30 W.m^-2, o que correspondeu a um aumento de temperatura a 1,0 m de profundidade da ordem de 0,22° C. Os perfis de temperatura obtidos durante este trabalho apresentaram deriva em sentidos alternados, durante o ciclo de período de um ano. Os perfis geotérmicos rasos são caracterizados por uma zona de fluxo de calor nulo, denominada "cotovelo", a partir da qual os valores de temperatura crescem com o aumento da profundidade. A profundidade do cotovelo é fundamentalmente influenciada pelo fluxo de calor gerado pela radiação solar incidente na superfície da área em estudo. Elaborou-se um modelo de evolução temporal o qual representa a estrutura térmica e suas variações, da zona compreendida de 0,02 a 10,0 m de profundidade. Os resultados oriundos deste modelo de transferência de calor por condução foram comparados com os obtidos em trabalhos de campo. Observa-se uma boa concordância entre esses resultados. Porém, o ajuste numérico é mais evidente para o período de setembro a fevereiro. Este modelamento poderá ser utilizado também para previsões de deriva dos perfis geotérmicos, desde que sejam conhecidos, a priori, os valores correspondentes às temperaturas médias mensais na superfície, o gradiente geotérmico regional e um dado perfil geotérmico. Este trabalho comprova que o fluxo de radiação solar incidente na superfície é a principal fonte de influência sobre os perfis geotérmicos rasos. A profundidade máxima dessa influência depende principalmente da magnitude desse fluxo, do grau de proteção superficial à incidência direta da radiação e da litologia do local em estudo.

Inversão vinculada de dados de EMAP e MT-1D anisotrópico

Uma técnica para a inversão de dados magnetotelúricos é apresentada neste trabalho. Dois tipos de dados são tratados aqui, dados gerados por modelos unidimensionais com anisotropia na condutividade das camadas e dados bi-dimensionais de levantamentos do método EMAP (ElectroMagnetic Array Profiling). Em ambos os casos fazemos a inversão usando vínculos aproximados de igualdade para estabilizar as soluções. Mostramos as vantagens e as limitações do uso destes vínculos nos processos de inversão. Mesmo vinculada a inversão ainda pode se tornar instável. Para inverter os dados 2-D do EMAP, apresentamos um processo que consiste de três partes: 1 – A construção de um modelo interpretativo e da aproximação inicial para a inversão a partir dos dados de seções de resistividade aparente filtradas pelo processo de filtragem do EMAP; 2 – a inclusão de uma camada de corpos pequenos aflorantes, chamada de camada destatic shift, aos modelos interpretativos para resolver as fontes de distorções estáticas que contaminam os dados; 3 – o uso dos vínculos aproximados de igualdade absoluta para estabilizar as soluções. Os dois primeiros passos nos permitem extrair o máximo de informação possível dos dados, enquanto que o uso dos vínculos de igualdade nos permite incluir informação a priori que possua significado físico e geológico. Com estes passos, obtemos uma solução estável e significativa. Estudaremos o método em dados sintéticos de modelos bi-dimensionais e em dados reais de uma linha de EMAP feita na Bacia do Paraná.

Mapeamento do relevo do embasamento de bacias sedimentares através da inversão gravimétrica vinculada

Apresentamos três novos métodos estáveis de inversão gravimétrica para estimar o relevo de uma interface arbitrária separando dois meios. Para a garantia da estabilidade da solução, introduzimos informações a priori sobre a interface a ser mapeada, através da minimização de um (ou mais) funcional estabilizante. Portanto, estes três métodos se diferenciam pelos tipos de informação físico-geológica incorporados. No primeiro método, denominado suavidade global, as profundidades da interface são estimadas em pontos discretos, presumindo-se o conhecimento a priori sobre o contraste de densidade entre os meios. Para a estabilização do problema inverso introduzimos dois vínculos: (a) proximidade entre as profundidades estimadas e verdadeiras da interface em alguns pontos fornecidas por furos de sondagem; e (b) proximidade entre as profundidades estimadas em pontos adjacentes. A combinação destes dois vínculos impõe uma suavidade uniforme a toda interface estimada, minimizando, simultaneamente em alguns pontos, os desajustes entre as profundidades conhecidas pelas sondagens e as estimadas nos mesmos pontos. O segundo método, denominado suavidade ponderada, estima as profundidades da interface em pontos discretos, admitindo o conhecimento a priori do contraste de densidade. Neste método, incorpora-se a informação geológica que a interface é suave, exceto em regiões de descontinuidades produzidas por falhas, ou seja, a interface é predominantemente suave porém localmente descontínua. Para a incorporação desta informação, desenvolvemos um processo iterativo em que três tipos de vínculos são impostos aos parâmetros: (a) ponderação da proximidade entre as profundidades estimadas em pontos adjacentes; (b) limites inferior e superior para as profundidades; e (c) proximidade entre todas as profundidades estimadas e um valor numérico conhecido. Inicializando com a solução estimada pelo método da suavidade global, este segundo método, iterativamente, acentua as feições geométricas presentes na solução inicial; ou seja, regiões suaves da interface tendem a tornar-se mais suaves e regiões abruptas tendem a tornar-se mais abruptas. Para tanto, este método atribui diferentes pesos ao vínculo de proximidade entre as profundidades adjacentes. Estes pesos são automaticamente atualizados de modo a acentuar as descontinuidades sutilmente detectadas pela solução da suavidade global. Os vínculos (b) e (c) são usados para compensar a perda da estabilidade, devida à introdução de pesos próximos a zero em alguns dos vínculos de proximidade entre parâmetros adjacentes, e incorporar a informação a priori que a região mais profunda da interface apresenta-se plana e horizontal. O vínculo (b) impõe, de modo estrito, que qualquer profundidade estimada é não negativa e menor que o valor de máxima profundidade da interface conhecido a priori; o vínculo (c) impõe que todas as profundidades estimadas são próximas a um valor que deliberadamente viola a profundidade máxima da interface. O compromisso entre os vínculos conflitantes (b) e (c) resulta na tendenciosidade da solução final em acentuar descontinuidades verticais e apresentar uma estimativa suave e achatada da região mais profunda. O terceiro método, denominado mínimo momento de inércia, estima os contrastes de densidade de uma região da subsuperfície discretizada em volumes elementares prismáticos. Este método incorpora a informação geológica que a interface a ser mapeada delimita uma fonte anômala que apresenta dimensões horizontais maiores que sua maior dimensão vertical, com bordas mergulhando verticalmente ou em direção ao centro de massa e que toda a massa (ou deficiência de massa) anômala está concentrada, de modo compacto, em torno de um nível de referência. Conceitualmente, estas informações são introduzidas pela minimização do momento de inércia das fontes em relação ao nível de referência conhecido a priori. Esta minimização é efetuada em um subespaço de parâmetros consistindo de fontes compactas e apresentando bordas mergulhando verticalmente ou em direção ao centro de massa. Efetivamente, estas informações são introduzidas através de um processo iterativo inicializando com uma solução cujo momento de inércia é próximo a zero, acrescentando, em cada iteração, uma contribuição com mínimo momento de inércia em relação ao nível de referência, de modo que a nova estimativa obedeça a limites mínimo e máximo do contraste de densidade, e minimize, simultaneamente, os desajustes entre os dados gravimétricos observados e ajustados. Adicionalmente, o processo iterativo tende a "congelar" as estimativas em um dos limites (mínimo ou máximo). O resultado final é uma fonte anômala compactada em torno do nível de referência cuja distribuição de constraste de densidade tende ao limite superior (em valor absoluto) estabelecido a priori. Estes três métodos foram aplicados a dados sintéticos e reais produzidos pelo relevo do embasamento de bacias sedimentares. A suavidade global produziu uma boa reconstrução do arcabouço de bacias que violam a condição de suavidade, tanto em dados sintéticos como em dados da Bacia do Recôncavo. Este método, apresenta a menor resolução quando comparado com os outros dois métodos. A suavidade ponderada produziu uma melhoria na resolução de relevos de embasamentos que apresentam falhamentos com grandes rejeitos e altos ângulos de mergulho, indicando uma grande potencialidade na interpretação do arcabouço de bacias extensionais, como mostramos em testes com dados sintéticos e dados do Steptoe Valley, Nevada, EUA, e da Bacia do Recôncavo. No método do mínimo momento de inércia, tomou-se como nível de referência o nível médio do terreno. As aplicações a dados sintéticos e às anomalias Bouguer do Graben de San Jacinto, California, EUA, e da Bacia do Recôncavo mostraram que, em comparação com os métodos da suavidade global e ponderada, este método estima com excelente resolução falhamentos com pequenos rejeitos sem impor a restrição da interface apresentar poucas descontinuidades locais, como no método da suavidade ponderada.

Inversão de momentos de fonte em métodos potenciais

A inversão de momentos de fonte gravimétrica tridimensional é analisada em duas situações. Na primeira se admite conhecer apenas a anomalia. Na segunda se admite conhecer, além da anomalia, informação a priori sobre o corpo anômalo. Sem usar informação a priori, mostramos que é possível determinar univocamente todo momento, ou combinação linear de momentos, cujo núcleo polinomial seja função apenas das coordenadas Cartesianas que definem o plano de medida e que tenha Laplaciano nulo. Além disso, mostramos que nenhum momento cujo núcleo polinomial tenha Laplaciano não nulo pode ser determinado. Por outro lado, informação a priori é implicitamente introduzida se o método de inversão de momentos se baseia na aproximação da anomalia pela série truncada obtida de sua expansão em multipolos. Dado um centro de expansão qualquer, o truncamento da série impõe uma condição de regularização sobre as superfícies equipotenciais do corpo anômalo, que permite estimar univocamente os momentos e combinações lineares de momentos que são os coeficientes das funções-bases da expansão em multipolos. Assim, uma distribuição de massa equivalente à real é postulada, sendo o critério de equivalência especificado pela condição de ajuste entre os campos observado e calculado com a série truncada em momentos de uma ordem máxima pré-estabelecida. Os momentos da distribuição equivalente de massa foram identificados como a solução estacionária de um sistema de equações diferenciais lineares de 1a. ordem, para a qual se asseguram unicidade e estabilidade assintótica. Para a série retendo momentos até 2a. ordem, é implicitamente admitido que o corpo anômalo seja convexo e tenha volume finito, que ele esteja suficientemente distante do plano de medida e que a sua distribuição espacial de massa apresente três planos ortogonais de simetria. O método de inversão de momentos baseado na série truncada (IMT) é adaptado para o caso magnético. Para este caso, mostramos que, para assegurar unicidade e estabilidade assintótica, é suficiente pressupor, além da condição de regularização, a condição de que a magnetização total tenha direção e sentido constantes, embora desconhecidos. O método IMT baseado na série de 2a. ordem (IMT2) é aplicado a anomalias gravimétricas e magnéticas tridimensionais sintéticas. Mostramos que se a fonte satisfaz as condições exigidas, boas estimativas da sua massa ou vetor momento de dipolo anômalo total, da posição de seu centro de massa ou de momento de dipolo e das direções de seus três eixos principais são obtidas de maneira estável. O método IMT2 pode falhar parcialmente quando a fonte está próxima do plano de medida ou quando a anomalia tem efeitos localizados e fortes de um corpo pequeno e raso e se tenta estimar os parâmetros de um corpo grande e profundo. Definimos por falha parcial a situação em que algumas das estimativas obtidas podem não ser boas aproximações dos valores verdadeiros. Nas duas situações acima descritas, a profundidade do centro da fonte (maior) e as direções de seus eixos principais podem ser erroneamente estimadas, embora que a massa ou vetor momento de dipolo anômalo total e a projeção do centro desta fonte no plano de medida ainda sejam bem estimados. Se a direção de magnetização total não for constante, o método IMT2 pode fornecer estimativas erradas das direções dos eixos principais (mesmo se a fonte estiver distante do plano de medida), embora que os demais parâmetros sejam bem estimados. O método IMT2 pode falhar completamente se a fonte não tiver volume finito. Definimos por falha completa a situação em que qualquer estimativa obtida pode não ser boa aproximação do valor verdadeiro. O método IMT2 é aplicado a dados reais gravimétricos e magnéticos. No caso gravimétrico, utilizamos uma anomalia situada no estado da Bahia, que se supõe ser causada por um batólito de granito. Com base nos resultados, sugerimos que as massas graníticas geradoras desta anomalia tenham sido estiradas na direção NNW e adelgaçadas na direção vertical durante o evento compressivo que causou a orogênese do Sistema de Dobramentos do Espinhaço. Além disso, estimamos que a profundidade do centro de massa da fonte geradora é cerca de 20 km. No caso magnético, utilizamos a anomalia de um monte submarino situado no Golfo da Guiné. Com base nos resultados, estimamos que o paleopolo magnético do monte submarino tem latitude 50°48'S e longitude 74°54'E e sugerimos que não exista contraste de magnetização expressivo abaixo da base do monte submarino.

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.

Modelagem sísmica e inversão na presença de anisotropia

A necessidade da adoção de modelos elásticos anisotrópicos, no contexto da sísmica de exploração, vem crescendo com o advento de novas técnicas de aquisição de dados como VSP, walkway VSP, tomografia poço a poço e levantamentos sísmicos com grande afastamento. Meios anisotrópicos, no contexto da sísmica de exploração, são modelos efetivos para explicar a propagação de ondas através de meios que apresentam padrões de heterogeneidade em escala muito menor que o comprimento de onda das ondas sísmicas. Particularmente, estes modelos são muito úteis para explicar o dado sísmico mais robusto que são as medidas de tempo de trânsito. Neste trabalho, são investigados aspectos da propagação de ondas, traçado de raios e inversão de tempos de trânsito em meios anisotrópicos. É estudada a propagação de ondas SH em meios anisotrópicos estratificados na situação mais geral onde estas ondas podem ocorrer, ou seja, em meios monoclínicos com um plano vertical de simetria especular. É mostrado que o campo de ondas SH refletido a partir de um semi-espaço estratificado, não apresenta qualquer informação sobre a possível presença de anisotropia em subsuperfície. São apresentados métodos simples e eficientes para o traçado de raios em 3D através de meios anisotrópicos estratificados, baseados no princípio de Fermat. Estes métodos constituem o primeiro passo para o desenvolvimento de algoritmos de inversão de tempos de trânsito para meios anisotrópicos em 3D, a partir de dados de VSP e walkaway VSP. Esta abordagem é promissora para determinação de modelos de velocidade, que são necessários para migração de dados sísmicos 3D na presença de anisotropia. É efetuada a análise da inversão tomográfica não linear, para meios estratificados transversalmente isotrópicos com um eixo de simetria vertical(TIV). As limitações dos dados de tempo de trânsito de eventos qP para determinação das constantes elásticas, são estabelecidas e caracterizados os efeitos da falta de cobertura angular completa na inversão tomográfica. Um algoritmo de inversão foi desenvolvido e avaliado em dados sintéticos. A aplicação do algoritmo a dados reais demonstra a consistência de meios TIV. Esta abordagem é útil para casos onde há informação a priori sobre a estratificação quase plana das formações e onde os próprios dados do levantamento poço a poço apresentam um alto grau de simetria especular em relação a um plano vertical. Também pode ser útil em interpretações preliminares, onde a estimativa de um meio estratificado, serve como modelo de fundo para se efetuar análises mais detalhadas, por exemplo, como um modelo de velocidades anisotrópico para migração, ou como um modelo de calibração para análises de AVO.