Inversão de dados de sísmica de refração profunda a partir da curva tempo-distância

O trabalho em pauta tem como objetivo o modelamento da crosta, através da inversão de dados de refração sísmica profunda, segundo camadas planas horizontais lateralmente homogêneas, sobre um semi-espaço. O modelo direto é dado pela expressão analítica da curva tempo-distância como uma função que depende da distância fonte-estação e do vetor de parâmetros velocidades e espessuras de cada camada, calculado segundo as trajetórias do raio sísmico, regidas pela Lei de Snell. O cálculo dos tempos de chegada por este procedimento, exige a utilização de um modelo cujas velocidades sejam crescentes com a profundidade, de modo que a ocorrência das camadas de baixa velocidade (CBV) é contornada pela reparametrização do modelo, levando-se em conta o fato de que o topo da CBV funciona apenas como um refletor do raio sísmico, e não como refrator. A metodologia de inversão utilizada tem em vista não só a determinação das soluções possíveis, mas também a realização de uma análise sobre as causas responsáveis pela ambiguidade do problema. A região de pesquisa das prováveis soluções é vinculada segundo limites superiores e inferiores para cada parâmetro procurado, e pelo estabelecimento de limites superiores para os valores de distâncias críticas, calculadas a partir do vetor de parâmetros. O processo de inversão é feito utilizando-se uma técnica de otimização do ajuste de curvas através da busca direta no espaço dos parâmetros, denominado COMPLEX. Esta técnica apresenta a vantagem de poder ser utilizada com qualquer função objeto, e ser bastante prática na obtenção de múltiplas soluções do problema. Devido a curva tempo-distância corresponder ao caso de uma multi-função, o algoritmo foi adaptado de modo a minimizar simultaneamente várias funções objetos, com vínculos nos parâmetros. A inversão é feita de modo a se obter um conjunto de soluções representativas do universo existente. Por sua vez, a análise da ambiguidade é realizada pela análise fatorial modo-Q, através da qual é possível se caracterizar as propriedades comuns existentes no elenco das soluções analisadas. Os testes com dados sintéticos e reais foram feitos tendo como aproximação inicial ao processo de inversão, os valores de velocidades e espessuras calculados diretamente da interpretação visual do sismograma. Para a realização dos primeiros, utilizou-se sismogramas calculados pelo método da refletividade, segundo diferentes modelos. Por sua vez, os testes com dados reais foram realizados utilizando-se dados extraídos de um dos sismogramas coletados pelo projeto Lithospheric Seismic Profile in Britain (LISPB), na região norte da Grã-Bretanha. Em todos os testes foi verificado que a geometria do modelo possui um maior peso na ambiguidade do problema, enquanto os parâmetros físicos apresentam apenas suaves variações, no conjunto das soluções obtidas.

A ictiofauna de poças rochosas de maré: padrões de distribuição associados ao espaço e ambiente

Em estudos com comunidades biológicas, a busca por padrões de ocorrência e distribuição das espécies vêm ganhando espaço em pesquisas realizadas nas últimas décadas. O conhecimento da estrutura das comunidades e suas relações, associados às interações desses organismos com o meio ambiente, fornecem subsídios necessários para a manutenção dessas comunidades, bem como informações relevantes em planos de conservação e manejo. Sendo assim, o objetivo desta pesquisa foi avaliar o efeito do ambiente e do espaço na estruturação da assembleia de peixes de poças rochosas de maré, descrevendo o padrão de ocupação das espécies nesse ambiente, testando a hipótese de que o espaço possui maior efeito na estruturação da ictiofauna. Foram amostradas 80 poças rochosas ao longo da Zona Costeira Amazônica, sendo 40 no período de maior precipitação e 40 no período de estiagem de 2011. Estas poças estão localizadas em cinco praias do litoral paraense e foram mensuradas quanto ao volume, a distância da margem e as variáveis físico-químicas da água, como pH, temperatura e salinidade. Foram coligidos 1.311 peixes, sendo 633 indivíduos na chuva e 648 na estiagem. Os indivíduos estão distribuídos em nove ordens, 14 famílias e 21 espécies, sendo a Ordem Perciformes a ordem mais abundante, com B. soporator e L. jocu sendo as espécies mais representativas. Os resultados obtidos na rotina BioEnv evidenciaram que as variáveis abióticas pH, temperatura e volume são responsáveis pela estruturação da comunidade. O cálculo da diversidade β evidenciou uma ampla variação de dados, abrangendo desde a substituição total das espécies até comunidades idênticas entre si, quanto à ocorrência e abundância. A DCA evidenciou que o período pluviométrico não é determinante na distribuição das espécies, uma vez que não há diferença na composição entre os dois períodos. O efeito do espaço e do ambiente sobre a comunidade foi pequeno, uma vez que a maior parte da variação não foi explicada utilizando-se todas as variáveis ambientais, como em uma segunda análise, em que só foram utilizadas as variáveis definidas pela rotina BioEnv. Seguindo o mesmo padrão, a pRDA apenas para as variáveis mais correlacionadas evidenciou que: 6% respondem ao ambiente, 4% ao espaço, 2% ao ambiente e espaço e 88% são outros fatores. Os resultados obtidos nessas análises comprovam que os fatores abióticos mensurados, bem como a distância entre as amostras, não são determinantes na composição e distribuição da ictiofauna em poças rochosas, refutando assim a hipótese de que o espaço possui efeito na estruturação da ictiofauna, uma vez que as condições ambientais sofrem constantes perturbações ocasionadas pelo movimento de marés. Sendo assim, as espécies que ocupam esse ambiente se mostram adaptadas, e por isso, não teriam sua distribuição afetada pela variação dos parâmetros ambientais nessa escala de estudo, respondendo somente ao efeito do espaço.

Propriedades magnéticas de nanoestruturas adsorvidas em superfícies metálicas

Neste trabalho, utilizamos o método de primeiros princípios RS-LMTO-ASA (Real Space – Linear Muffin-Tin Orbital - Atomic Sphere Approximation) baseado na Teoria do Funcional da Densidade (DFT - Density Functional Theory) e implementado para o cálculo de estruturas magnéticas não-colineares, para investigar as propriedades magnéticas de nanoestruturas adsorvidas em superfícies metálicas. Consideramos aglomerados com diferentes geometrias e tamanhos como adátomos, dímeros, trímeros, nanofios e nanoestruturas de geometria triangular de Fe, Fe-Co e Fe-Pt adsorvidos sobre a superfície de Pt(111) e tratamos também nanoestruturas de Mn sobre a superfície de Ag(111). Mostramos que os nanofios de Fe-Co sobre a superfície de Pt(111) apresentam um ordenamento ferromagnético. Devido à redução do número de coordenação presente na superfície, os momentos de spin e orbital nos sítios de Fe e Co mostram-se elevados comparados com os respectivos valores dos momentos destes metais como bulk. Analisamos também como estes momentos variam em função da concentração destes elementos nos nanofios. Para os sistemas compostos por nanofios Fe-Pt adsorvidos em Pt(111), mostramos que é possível sintonizar as interações de troca entre os adátomos magnéticos Fe através da introdução de um diferente número de átomos Pt para ligá-los. Por exemplo, a interação de troca entre os adátomos de Fe pode ser consideravelmente aumentada pela introdução de cadeias de Pt que os conectem e tanto configurações ferromagnéticas, antiferromagnéticas ou não-colineares entre os adátomos de Fe podem ser estabilizadas, dependendo da espessura do espaçador Pt. Para os aglomerados Mn sobre a Ag(111) mostramos que a interação de troca entre os sítios de Mn depende não somente da distância entre os átomos, mas também do número de coordenação de cada sítio. Desta forma, verificamos um magnetismo não-colinear nestas nanoestruturas causado tanto por frustração geométrica, quanto pela competição de interações de curto e longo alcance. Nossos resultados estão em boa concordância com os resultados experimentais da literatura e com os resultados teóricos obtidos por outros métodos, quando existentes.

Modelagem numérica de dados MCSEM 3D usando computação paralela

Desenvolvemos a modelagem numérica de dados sintéticos Marine Controlled Source Electromagnetic (MCSEM) usada na exploração de hidrocarbonetos para simples modelos tridimensionais usando computação paralela. Os modelos são constituidos de duas camadas estrati cadas: o mar e o sedimentos encaixantes de um delgado reservatório tridimensional, sobrepostas pelo semi-espaço correspondente ao ar. Neste Trabalho apresentamos uma abordagem tridimensional da técnica dos elementos nitos aplicada ao método MCSEM, usando a formulação da decomposição primária e secundária dos potenciais acoplados magnético e elétrico. Num pós-processamento, os campos eletromagnéticos são calculados a partir dos potenciais espalhados via diferenciação numérica. Exploramos o paralelismo dos dados MCSEM 3D em um levantamento multitransmissor, em que para cada posição do transmissor temos o mesmo processo de cálculos com dados diferentes. Para isso, usamos a biblioteca Message Passing Interface (MPI) e o modelo servidor cliente, onde o processador administrador envia os dados de entradas para os processadores clientes computar a modelagem. Os dados de entrada são formados pelos parâmetros da malha de elementos nitos, dos transmissores e do modelo geoelétrico do reservatório. Esse possui geometria prismática que representa lentes de reservatórios de hidrocarbonetos em águas profundas. Observamos que quando a largura e o comprimento horizontais desses reservatório têm a mesma ordem de grandeza, as resposta in-line são muito semelhantes e conseqüentemente o efeito tridimensional não é detectado. Por sua vez, quando a diferença nos tamanhos da largura e do comprimento do reservatório é signi cativa o efeito 3D é facilmente detectado em medidas in-line na maior dimensão horizontal do reservatório. Para medidas na menor dimensão esse efeito não é detectável, pois, nesse caso o modelo 3D se aproxima de um modelo bidimensional. O paralelismo dos dados é de rápida implementação e processamento. O tempo de execução para a modelagem multitransmissor em ambiente paralelo é equivalente ao tempo de processamento da modelagem para um único transmissor em uma máquina seqüêncial, com o acréscimo do tempo de latência na transmissão de dados entre os nós do cluster, o que justi ca o uso desta metodologia na modelagem e interpretação de dados MCSEM. Devido a reduzida memória (2 Gbytes) em cada processador do cluster do departamento de geofísica da UFPA, apenas modelos muito simples foram executados.

Inversão de dados sísmicos de reflexão a partir da curva do tempo de trânsito

A presente Dissertação de Mestrado tem como objetivo o estudo do problema de inversão sísmica baseada em refletores planos para arranjo fonte-comum (FC) e ponto-médiocomum (PMC). O modelo direto é descrito por camadas homogêneas, isotrópicas com interfaces plano-horizontais. O problema é relacionado ao empilhamento NMO baseado na otimização da função semblance, para seções PMC corrigidas de sobretempo normal (NMO). O estudo foi baseado em dois princípios. O primeiro princípio adotado foi de combinar dois grupos de métodos de inversão: um Método Global e um Método Local. O segundo princípio adotado foi o de cascata, segundo a teoria Wichert-Herglotz-Bateman, que estabelece que para conhecer uma camada inferior tem-se que conhecer primeiro a camada superior (dissecação). A aplicação do estudo é voltada à simulação sísmica de Bacia Sedimentar do Solimões e de Bacia Marinha para se obter uma distribuição local 1D de velocidades e espessuras para a subsuperfície em horizontes alvo. Sendo assim, limitamos a inversão entre 4 e 11 refletores, uma vez que na prática a indústria limita uma interpretação realizada apenas em número equivalente de 3 a 4 refletores principais. Ressalta-se que este modelo é aplicável como condição inicial ao imageamento de seções sísmicas em regiões geologicamente complexas com variação horizontal suave de velocidades. Os dados sintéticos foram gerados a partir dos modelos relacionados a informações geológicas, o que corresponde a uma forte informação a priori no modelo de inversão. Para a construção dos modelos relacionados aos projetos da Rede Risco Exploratório (FINEP) e de formação de recursos humanos da ANP em andamento, analisamos os seguintes assuntos relevantes: (1) Geologia de bacias sedimentares terrestre dos Solimões e ma rinha (estratigráfica, estrutural, tectônica e petrolífera); (2) Física da resolução vertical e horizontal; e (3) Discretização temporal-espacial no cubo de multi-cobertura. O processo de inversão é dependente do efeito da discretização tempo-espacial do campo de ondas, dos parâmetros físicos do levantamento sísmico, e da posterior reamostragem no cubo de cobertura múltipla. O modelo direto empregado corresponde ao caso do operador do empilhamento NMO (1D), considerando uma topografia de observação plana. O critério básico tomado como referência para a inversão e o ajuste de curvas é a norma 2 (quadrática). A inversão usando o presente modelo simples é computacionalmente atrativa por ser rápida, e conveniente por permitir que vários outros recursos possam ser incluídos com interpretação física lógica; por exemplo, a Zona de Fresnel Projetada (ZFP), cálculo direto da divergência esférica, inversão Dix, inversão linear por reparametrização, informações a priori, regularização. A ZFP mostra ser um conceito út il para estabelecer a abertura da janela espacial da inversão na seção tempo-distância, e representa a influência dos dados na resolução horizontal. A estimativa da ZFP indica uma abertura mínima com base num modelo adotado, e atualizável. A divergência esférica é uma função suave, e tem base física para ser usada na definição da matriz ponderação dos dados em métodos de inversão tomográfica. A necessidade de robustez na inversão pode ser analisada em seções sísmicas (FC, PMC) submetida a filtragens (freqüências de cantos: 5;15;75;85; banda-passante trapezoidal), onde se pode identificar, comparar e interpretar as informações contidas. A partir das seções, concluímos que os dados são contaminados com pontos isolados, o que propõe métodos na classe dos considerados robustos, tendo-se como referência a norma 2 (quadrados- mínimos) de ajuste de curvas. Os algoritmos foram desenvolvidos na linguagem de programação FORTRAN 90/95, usando o programa MATLAB para apresentação de resultados, e o sistema CWP/SU para modelagem sísmica sintética, marcação de eventos e apresentação de resultados.