Estudo comparativo entre estereotomografia e da tomografia da onda NIP: aplicação em dados sintéticos e reais

A determinação de um acurado modelo de velocidades é um requisito fundamental para a realização do imageamento sísmico. Métodos novos como a Estereotomografia préempilhamento e a Tomografia da onda NIP são ferramentas poderosas e bastante sugestivas para este propósito. Basicamente, a Estereotomografia pré-empilhamento se baseia no conceito de eventos localmente coerentes interpretados como reflexões primárias e associados com pares de segmentos de raios, que tem um mesmo ponto de reflexão em profundidade. Na Tomografia da onda NIP um evento sísmico é representado por uma onda hipotética NIP, que está relacionada a um ponto de reflexão em profundidade. Os atributos da onda NIP são determinados no decorrer do Empilhamento de Superfície de Reflexão Comum (empilhamento CRS). Este trabalho tem como objetivo, fazer um estudo comparativo de ambos os métodos de determinação do modelo de velocidades em profundidade. Assim, é realizada uma revisão dos fundamentos teóricos de ambos os métodos tomográficos, destacando as principais diferenças e aplicando estas aproximações em um dado sintético e um dado real marinho (linha sísmica 214-2660 da Bacia do Jequitinhonha). Para avaliar os modelos de velocidades encontrados pelas aproximações, foi utilizada a migração pré-empilhamento em profundidade do tipo Kirchhoff e também as famílias de imagem comum (CIG). Os resultados mostraram que ambos os métodos tomográficos fornecem modelos de velocidades representativos. Contudo, constatou-se que a estereotomografia tem melhor desempenho em meios com variações laterais de velocidades, porém, aplicável somente em dados pré-empilhados com alta razão sinal/ruído.

Vibrações naturais em um sistema de interação viga-água incluindo o efeito de onda de superfície

Este trabalho tem como objetivo o estudo modal de um sistema estruturafluído, modelado pela equação de Euler-Bernoulli, para uma viga elástica, sujeita a pressão da água, e por ondas de superfície livre. O sistema acoplado possui condições para o domínio sólido (viga fixa-livre), para o domínio fluído (impermeabilidade e rigidez inferior), com ondas de superfície, e de interfacec fluido-estrutura (condições de continunidade na deflexão, ângulo de rotação, força de corte interno e momento de curvatura). Para deteerminar as vibrações livres e deslocamento no seco e no molhado, utiliza-se o método espectral para eliminar a dependência oscilatória temporal e concentrar-se na determinação dos modos através do estudo de problemas de contornos espaciais. Os modos podem ser calculados com o uso da base clássica de Euler ou da base dinâmica gerada pela resposta impulso. Foram feitas simulações para um material específico, e apresentados os resultados obtidos.

Vibrações naturais em um sistema de interação viga-água sem o efeito de onda de superfície livre

Neste trabalho é apresentado um estudo sobre o comportamento dinâmico de sistemas flexíveis com interação viga-água. Especificamente, é escolhida uma barragem freqüentemente encontrada em problemas de engenharia de irrigação e represas. O modelo matemático utilizado para descrever tal fenômeno tem como variáveis principais a pressão hidrodinâmica e os deslocamentos horizontais da barragem. Uma vez formuladas as equações governantes de tal sistema, que resultam em uma equação de onda para a pressão e uma equação de Euler-Bernoulli para os deslocamentos, é utilizada a técnica de separação de variáveis para proceder à sua solução. Por simplicidade, assumiu-se que no modelo não havia o efeito de onda da superfície livre. Primeiro, foi resolvida a equação para a pressão, sendo calculadas as suas componentes temporais e espaciais. Depois, foi resolvida a equação de Euler-Bernoulli junto com a condição de interface. Foram simulações para o material específico, e apresentados os resultados conseguidos. No final do trabalho são enunciadas as conclusões pertinentes.

Estudo da estrutura subsuperficial da província Borborema com correlação de ruído sísmico

Ambient seismic noise has traditionally been considered as an unwanted perturbation in seismic data acquisition that "contaminates" the clean recording of earthquakes. Over the last decade, however, it has been demonstrated that consistent information about the subsurface structure can be extracted from cross-correlation of ambient seismic noise. In this context, the rules are reversed: the ambient seismic noise becomes the desired seismic signal, while earthquakes become the unwanted perturbation that needs to be removed. At periods lower than 30 s, the spectrum of ambient seismic noise is dominated by microseism, which originates from distant atmospheric perturbations over the oceans. The microsseism is the most continuous seismic signal and can be classified as primary – when observed in the range 10-20 s – and secondary – when observed in the range 5-10 s. The Green‘s function of the propagating medium between two receivers (seismic stations) can be reconstructed by cross-correlating seismic noise simultaneously recorded at the receivers. The reconstruction of the Green‘s function is generally proportional to the surface-wave portion of the seismic wavefield, as microsseismic energy travels mostly as surface-waves. In this work, 194 Green‘s functions obtained from stacking of one month of daily cross-correlations of ambient seismic noise recorded in the vertical component of several pairs of broadband seismic stations in Northeast Brazil are presented. The daily cross-correlations were stacked using a timefrequency, phase-weighted scheme that enhances weak coherent signals by reducing incoherent noise. The cross-correlations show that, as expected, the emerged signal is dominated by Rayleigh waves, with dispersion velocities being reliably measured for periods ranging between 5 and 20 s. Both permanent stations from a monitoring seismic network and temporary stations from past passive experiments in the region are considered, resulting in a combined network of 33 stations separated by distances between 60 and 1311 km, approximately. The Rayleigh-wave, dispersion velocity measurements are then used to develop tomographic images of group velocity variation for the Borborema Province of Northeast Brazil. The tomographic maps allow to satisfactorily map buried structural features in the region. At short periods (~5 s) the images reflect shallow crustal structure, clearly delineating intra-continental and marginal sedimentary basins, as well as portions of important shear zones traversing the Borborema Province. At longer periods (10 – 20 s) the images are sensitive to deeper structure in the upper crust, and most of the shallower anomalies fade away. Interestingly, some of them do persist. The deep anomalies do not correlate with either the location of Cenozoic volcanism and uplift - which marked the evolution of the Borborema Province in the Cenozoic - or available maps of surface heat-flow, and the origin of the deep anomalies remains enigmatic.

Dispersão das ondas de Rayleigh na Plataforma Sulamericana

Este trabalho representa um estudo de dispersão da componente vertical da onda de superfície de Rayleigh com trajetórias na plataforma Sulamericana. Os registros utilizados são provenientes das estações localizadas no território Brasileiro; sendo a do Rio de Janeiro (RDJ), a de Brasília (BDF), a de Caicó (CAI) e a de Belém (BEB), pois estas são as únicas estações sismológicas no Brasil que têm sensores de período longo e que servem para o estudo de dispersão no intervalo de 4 a 50 segundos, aqui realizados. Os terremotos utilizados estão localizados ao longo da parte leste da cadeia Andina e dentro da plataforma Sulamericana com trajetórias tipicamente continental. Foram selecionados 34 eventos com a utilização dos seguintes critérios práticos: a localização, a magnitude m_b e a profundidade, ocorridos durante o período de Janeiro de 1978 até Junho de 1987. O estudo de dispersão aqui abordado significa a determinação da velocidade de grupo e das amplitudes espectrais correspondentes aos harmônicos fundamental e primeiro superior. Normalmente os harmônicos de ordem segunda ou maior são raramente disponíveis na observação. Dois tipos de medidas foram feitas: (i) velocidade de grupo vs. período e (ii) amplitude vs. período. Os estudos de dispersão são fundamentais para determinação da estrutura da crosta e manto superior que estão diretamente relacionados com os fenômenos geológicos. Neste trabalho, regionalização é definida como a identificação das diferentes formas de curvas de dispersão, que estão relacionadas com as trajetórias epicentro-estação ao longo da plataforma Sulamericana e que venham ter uma correlação geológica como está descrito no item 4.3 deste trabalho. A distribuição dos epicentros se faz desde o extremo sul da Argentina até o extremo norte da Venezuela, objetivando iniciar com este trabalho uma sistemática voltada aos estudos de regionalização da plataforma Sulamericana na nossa instituição. Neste trabalho foram observados três tipos distintos de curvas em 27 trajetórias e agrupadas por famílias 1,2 e 3 respectivamente, onde procurou-se correlacionar suas diferentes formas com a geologia regional da plataforma Sulamericana. A obtenção da curva de dispersão foi feita através da técnica do filtro múltiplo (Dziewonski et al, 1969). Este filtro tem a propriedade de separar os harmônicos através das suas velocidades de grupo para cada frequência selecionada, e também de recuperar as amplitudes características dos harmônicos (Herrmann, 1973). O desenvolvimento teórico do filtro bem como suas limitações e forma de uso são tratados por Dziewonski et al (1972). Como parte do trabalho há a implantação, adaptações e o desenvolvimento de parte do fluxograma do filtro múltiplo, bem como a estruturação da digitalização dos dados para o processamento e interpretação não-automática dos resultados do processamento.

Análise de onda completa de superfície seletiva em freqüência do tipo anteparo duplo

This work presents a theoretical analysis and numerical and experimental results of the scattering characteristics of frequency selective surfaces, using elements of type patch perfectly conductor. The structures are composed of two frequency selective surfaces on isotropic dielectric substrates cascaded, separated by a layer of air. The analysis is performed using the method of equivalent transmission line in combination with the Galerkin method, to determine the transmission and reflection characteristics of the structures analyzed. Specifically, the analysis uses the impedance method, which models the structure by an equivalent circuit, and applies the theory of transmission lines to determine the dyadic Green's function for the cascade structure. This function relates the incident field and surface current densities. These fields are determined algebraically by means of potential incidents and the imposition of the continuity of the fields in the dielectric interfaces. The Galerkin method is applied to the numerical determination of the unknown weight coefficients and hence the unknown densities of surface currents, which are expanded in terms of known basis functions multiplied by these weight coefficients. From the determination of these functions, it becomes possible to obtain numerical scattered fields at the top and bottom of the structures and characteristics of transmission and reflection of these structures. At work, we present numerical and experimental results for the characteristics of transmission and reflection. Comparisons were made with other results presented in literature, and it was observed a good agreement in the cases presented suggestions continuity of the work are presented

Estudo de acoplamento de luz a sistemas multicamadas : plasmons de superfícies e guias de onda

O acoplamento de radiação óptica em sistemas multicamadas tem sido objeto de diversas pesquisas, não somente acadêmicas, mas também para aplicações industriais, tanto na área de sensores ópticos, como na de processadores ópticos e na de comunicações ópticas. Existe uma variedade de técnicas que são exploradas nestes estudos. Nesta tese, focalizamos nossa atenção no acoplamento de radiação laser ou mesmo de luz branca incoerente a um filme tipo multicamadas e estudamos os mecanismos físicos que governam as reflexões e as perdas por absorção nas multicamadas, efeitos que normalmente não parecem quando a incidência ocorre diretamente do ar para o filme. A técnica que exploramos é conhecida como reflexão interna total atenuada, ATR. A presença de um filme fino metálico permite o estudo experimental de plasmons de superfície e do acoplamento a modos guiados do sistema multicamadas. Além dos estudos experimentais, apresentamos um cálculo teórico iterativo para a refletividade dos filmes do tipo multicamadas para um número qualquer de camadas, que apresenta vantagens computacionais no ajuste dos dados experimentais ou em simulações. Esta contribuição não necessita as aproximações encontradas em um grande número de trabalhos que envolvem sistemas mais simples. Apresentamos também o cálculo do fluxo de energia dentro de cada camada individual do sistema, o que nos permite determinar o tipo de modo acoplado e a sua localização. O método foi aplicado a diversos sistemas, quando uma das camadas foi modificada. Estes estudos foram realizados como função do ângulo de incidência e do comprimento de onda da radiação incidente para uma variedade de sistemas multicamadas. Nossa simulação teórica se mostra bastante adequada e útil para a projeção de sistemas multicamadas complexos com camadas metálicas e dielétricas para sensores óticos.

Controle de qualidade e dosimetria em equipamentos de tomografia computadorizada

OBJETIVO: Avaliação de condições dos equipamentos e dosimetria em setores de tomografia computadorizada utilizando protocolos de cabeça, abdome e coluna lombar em pacientes adultos (em três equipamentos distintos) e pediátricos com até um ano e meio de vida (em um dos equipamentos avaliados). MATERIAIS E MÉTODOS: Foram estimados o índice de dose em tomografia computadorizada e a dose média em cortes múltiplos, em exames com pacientes adultos, em três distintos equipamentos. Ainda foram estimadas as doses na superfície de entrada e as doses absorvidas em exame de cabeça para pacientes adultos e pediátricos em um dos equipamentos avaliados. RESULTADOS: Foram realizados testes de controle de qualidade, mecânicos, demonstrando que os equipamentos satisfazem as especificações de uso estabelecidas pelas normas vigentes. Os resultados da dosimetria mostraram que valores de dose média em cortes múltiplos excederam em até 109,0% os valores de níveis de referência, apresentando consideráveis variações entre os equipamentos avaliados neste estudo. As doses absorvidas obtidas com protocolos pediátricos são inferiores aos de pacientes adultos, apresentando redução de até 51,0% na tireoide. CONCLUSÃO: Neste estudo foram avaliadas as condições de operação de três equipamentos tomográficos, estabelecendo quais parâmetros devem ser trabalhados para a implantação de um programa de controle de qualidade nas instituições onde esta pesquisa foi desenvolvida.

Abordagem matemática de roll waves em escoamentos hiperconcentrados com superfície livre

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Geoprocessamento na identificação de ilhas de calor e influência do uso e ocupação do solo na temperatura aparente da superfície no município de Botucatu/SP

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Equação reduzida para ondas curtas na superfície da água

Pós-graduação em Física - IFT

Construção de um interferômetro de Bath para análise do comportamento da superfície de vidros usados na confecção de espelhos côncavos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Migração em profundidade pré-empilhamento utilizando os atributos cinemáticos do empilhamento por superfície de reflexão comum

O empilhamento por superfície de reflexão comum (ou empilhamento SRC), conhecido como empilhamento CRS, do inglês Commom reflection surface, constitui-se em um novo método para o processamento sísmico na simulação de seções afastamento nulo (AN) e afastamento comum (AC). Este método é baseado em uma aproximação paraxial hiperbólica de segunda ordem dos tempos de trânsito de reflexão na vizinhança de um raio central. Para a simulação de seção AN, o raio central é um raio normal, enquanto que para a simulação de uma seção AC o raio central é um raio de afastamento finito. Em adição à seção AN, o método de empilhamento SRC também fornece estimativas dos atributos cinemáticos do campo de onda, sendo aplicados, por exemplo, na determinação (por um processo de inversão) da velocidade intervalar, no cálculo do espalhamento geométrico, na estimativa da zona de Fresnel, e também na simulação de eventos de tempos de difrações, este último tendo uma grande importância para a migração pré-empilhamento. Neste trabalho é proposta uma nova estratégia para fazer uma migração em profundidade pré-empilhamento, que usa os atributos cinemáticos do campo de onda derivados do empilhamento SRC, conhecido por método CRS-PSDM, do inglês CRS based pre-stack depth migration. O método CRS-PSDM usa os resultados obtidos do método SRC, isto é, as seções dos atributos cinemáticos do campo de onda, para construir uma superfície de tempos de trânsito de empilhamento, ao longo da qual as amplitudes do dado sísmico de múltipla cobertura são somadas, sendo o resultado da soma atribuído a um dado ponto em profundidade, na zona alvo de migração que é definida por uma malha regular. Similarmente ao método convencional de migração tipo Kirchhoff (K-PSDM), o método CRS-PSDM precisa de um modelo de velocidade de migração. Contrário ao método K-PSDM, o método CRS-PSDM necessita apenas computar os tempos de trânsito afastamento nulo, ao seja, ao longo de um único raio ligando o ponto considerado em profundidade a uma dada posição de fonte e receptor coincidentes na superfície. O resultado final deste procedimento é uma imagem sísmica em profundidade dos refletores a partir do dado de múltipla cobertura.

Estimativa da radiação de onda longa atmosférica em áreas de floresta e de pastagem no sudoeste da Amazônia

A radiação de onda longa proveniente da atmosfera (Lin) é a componente do balanço de radiação mais difícil de ser medida. Na Amazônia praticamente não existem medidas regulares dessa componente, mesmo sendo uma importante variável no cálculo do balanço de radiação à superfície e muito usada para alimentar modelos climáticos. Tendo em vista a necessidade desses dados, o objetivo do presente trabalho é avaliar o desempenho de sete equações na estimativa da Lin para dias de céu claro em áreas de floresta (Reserva Biológica do Jaru, 10º4'48''S; 61º55'48''W) e de pastagem (Fazenda Nossa Senhora, 10º45'S; 62º21'W) no sudoeste da Amazônia. Medidas de radiação de onda longa atmosférica realizadas no período de junho de 2005 a maio de 2006 foram comparadas com as estimativas. As equações testadas tiveram desempenho satisfatório apenas durante a estação seca. As condições de alta nebulosidade, dominantes na estação chuvosa, restringiram a quantidade de dados utilizados na avaliação das equações. As equações que utilizam informações de temperatura do ar e pressão de vapor d'água para a estimativa da Lin tiveram melhor desempenho em relação às que utilizam apenas a temperatura do ar. As equações de Brutsaert (1975), Idso (1981) e Prata (1996) foram as que apresentaram melhor desempenho, apresentando os maiores índices de concordância, e sendo, portanto, as equações mais indicadas para a estimativa da radiação de onda longa atmosférica no sudoeste da Amazônia.

Empilhamento pelo método superfície de reflexão comum 2-D com topografia e introdução ao caso 3-D

O método de empilhamento sísmico CRS simula seções sísmicas ZO a partir de dados de cobertura múltipla, independente do macro-modelo de velocidades. Para meios 2-D, a função tempo de trânsito de empilhamento depende de três parâmetros, a saber: do ângulo de emergência do raio de reflexão normal (em relação à normal da superfície) e das curvaturas das frentes de onda relacionadas às ondas hipotéticas, denominadas NIP e Normal. O empilhamento CRS consiste na soma das amplitudes dos traços sísmicos em dados de múltipla cobertura, ao longo da superfície definida pela função tempo de trânsito do empilhamento CRS, que melhor se ajusta aos dados. O resultado do empilhamento CRS é assinalado a pontos de uma malha pré-definida na seção ZO. Como resultado tem-se a simulação de uma seção sísmica ZO. Isto significa que para cada ponto da seção ZO deve-se estimar o trio de parâmetros ótimos que produz a máxima coerência entre os eventos de reflexão sísmica. Nesta Tese apresenta-se fórmulas para o método CRS 2-D e para a velocidade NMO, que consideram a topografia da superfície de medição. O algoritmo é baseado na estratégia de otimização dos parâmetros de fórmula CRS através de um processo em três etapas: 1) Busca dos parâmetros, o ângulo de emergência e a curvatura da onda NIP, aplicando uma otimização global, 2) busca de um parâmetro, a curvatura da onda N, aplicando uma otimização global, e 3) busca de três parâmetros aplicando uma otimização local para refinar os parâmetros estimados nas etapas anteriores. Na primeira e segunda etapas é usado o algoritmo Simulated Annealing (SA) e na terceira etapa é usado o algoritmo Variable Metric (VM). Para o caso de uma superfície de medição com variações topográficas suaves, foi considerada a curvatura desta superfície no algoritmo do método de empilhamento CRS 2-D, com aplicação a dados sintéticos. O resultado foi uma seção ZO simulada, de alta qualidade ao ser comparada com a seção ZO obtida por modelamento direto, com uma alta razão sinal-ruído, além da estimativa do trio de parâmetros da função tempo de trânsito. Foi realizada uma nálise de sensibilidade para a nova função de tempo de trânsito CRS em relação à curvatura da superfície de medição. Os resultados demonstraram que a função tempo de trânsito CRS é mais sensível nos pontos-médios afastados do ponto central e para grandes afastamentos. As expressões da velocidade NMO apresentadas foram aplicadas para estimar as velocidades e as profundidades dos refletores para um modelo 2-D com topografia suave. Para a inversão destas velocidades e profundidades dos refletores, foi considerado o algoritmo de inversão tipo Dix. A velocidade NMO para uma superfície de medição curva, permite estimar muito melhor estas velocidades e profundidades dos refletores, que as velocidades NMO referidas as superfícies planas. Também apresenta-se uma abordagem do empilhamento CRS no caso 3-D. neste caso a função tempo de trânsito depende de oito parâmetros. São abordadas cinco estratégias de busca destes parâmetros. A combinação de duas destas estratégias (estratégias das três aproximações dos tempos de trânsito e a estratégia das configurações e curvaturas arbitrárias) foi aplicada exitosamente no empilhamento CRS 3-D de dados sintéticos e reais.