Modelagem sísmica por diferenças finitas em meios bidimensionais com difratores

Este trabalho tem por objetivo a modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações. Como parte deste estudo foi feita a análise numérica visando as condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas. Para a validação do programa de diferenças finitas foi feita a comparação cinemática com a Teoria do Raio para um modelo simples. O motivo deste estudo é ter uma melhor compreensão e controle sobre os problemas de modelagem, visando contribuir para a solução de problemas na interpretação de dados sísmicos. Segundo vários autores na literatura geológica, Derby (1877), Evans (1906), Paiva (1929) e Moura (1938). A Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o recente, atingindo 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çando espessuras de centenas de metros, perfazendo um volume total de 90.000 quilômetros cúbicos. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica, impossibilitando uma melhor interpretação dos horizontes refletores que se encontram abaixo das soleiras. Para representar situações geológicas desse tipo foi usado um modelo acústico de velocidades. Para o cálculo dos sismogramas foi utilizado um programa de diferenças finitas com aproximação de quarta ordem da equação da onda acústica no espaço e no tempo. As aplicações dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas têm melhorado 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 observado a influência significativa das múltiplas devido a camada de alta velocidade, o que faz com que haja maior perda de energia dificultando a interpretação dos alvos. Por esta razão recomendo a integração de dados de superfície com os dados de poço, com o objetivo de se obter uma melhor imagem dos alvos abaixo das soleiras de diabásio.

Inspeção de materiais compósitos utilizando ondas acústicas guiadas

Pós-graduação em Engenharia Elétrica - FEIS

Interpretação racional do ensaio DPL

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

Um algoritmo paralelo eficiente de migração reversa no tempo (rtm) 3d com granularidade fina

The reverse time migration algorithm (RTM) has been widely used in the seismic industry to generate images of the underground and thus reduce the risk of oil and gas exploration. Its widespread use is due to its high quality in underground imaging. The RTM is also known for its high computational cost. Therefore, parallel computing techniques have been used in their implementations. In general, parallel approaches for RTM use a coarse granularity by distributing the processing of a subset of seismic shots among nodes of distributed systems. Parallel approaches with coarse granularity for RTM have been shown to be very efficient since the processing of each seismic shot can be performed independently. For this reason, RTM algorithm performance can be considerably improved by using a parallel approach with finer granularity for the processing assigned to each node. This work presents an efficient parallel algorithm for 3D reverse time migration with fine granularity using OpenMP. The propagation algorithm of 3D acoustic wave makes up much of the RTM. Different load balancing were analyzed in order to minimize possible losses parallel performance at this stage. The results served as a basis for the implementation of other phases RTM: backpropagation and imaging condition. The proposed algorithm was tested with synthetic data representing some of the possible underground structures. Metrics such as speedup and efficiency were used to analyze its parallel performance. The migrated sections show that the algorithm obtained satisfactory performance in identifying subsurface structures. As for the parallel performance, the analysis clearly demonstrate the scalability of the algorithm achieving a speedup of 22.46 for the propagation of the wave and 16.95 for the RTM, both with 24 threads.

O estudo de caos na propagação do raio de som em um guia de onda no oceano

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

Análise do risco associado à rotura da barragem de Fernandilho. Modelação hidrodinâmica da onda de cheia

As cheias provocadas por roturas de barragens constituem eventos catastróficos que podem conduzir a perdas consideráveis de vidas humanas e de bens materiais. Para melhorar a segurança, em relação a este risco potencial, torna-se necessário avaliar a segurança da população e das estruturas, existentes ou a construir, e estimar a extensão dos danos em caso de acidente. Neste estudo, após a indicação da legislação em vigor, apresenta-se uma descrição sumária do modelo DamBreak. Descrevem-se as equações básicas que definem o modelo de formação da brecha e que caracterizam o processo de formação da cheia. Também se referem as equações de Saint-Venant, as quais caracterizam o processo de propagação da cheia. De seguida é feita uma descrição da barragem de Fernandilho, do leito do Barranco d

Análise de incerteza associada à determinação da velocidade de onda sônica em depósito de carvão obtida por perfilagem geofísica

A modelagem de um depósito mineral é realizada por meio da combinação de diversas fontes de informações. Dentre estas fontes pode-se citar a sísmica de reflexão. A sísmica de reflexão fornece dados de tempos de propagação de ondas sísmicas até essas serem refletidas pelas estruturas dos depósitos minerais. As profundidades dessas estruturas podem ser obtidas multiplicando-se os tempos pelas velocidades de propagação das ondas sísmicas. Normalmente, a velocidade de uma onda sísmica é determinada indiretamente por meio do processamento dos próprios dados sísmicos, o que pode gerar erros na interpretação de seções geológicas. A perfilagem geofísica é uma alternativa na determinação dessa velocidade, uma vez que a velocidade de onda acústica é obtida ao longo do furo perfilado e a velocidade de onda acústica pode ser relacionada com a velocidade de onda sísmica. As estimativas de valores de velocidade na região entre os furos perfilados permite as estimativas de valores de profundidade nessa região. Neste estudo, foram analisadas possibilidades de se estimar, em um grid, valores de velocidade e incertezas associadas a esses valores, utilizando-se ferramentas geoestatísticas. A simulação seqüencial Gaussiana, dentre as ferramentas analisadas, pareceu ser a mais adequada para obtenção de velocidades e incerteza dos valores estimados em cada nó do grid considerado. Para o caso abordado, alguns valores de profundidade da estrutura de interesse apresentaram variações significativas em função da incerteza de velocidade. Essas variações são muito importantes na execução de certos métodos de lavra subterrânea, o que enfatiza a importância da determinação da incerteza das velocidades estimadas. A metodologia é apresentada e ilustrada em um importante depósito de carvão em Queensland, Austrália.

Rompimento de barreiras: análise experimental e numérica na previsão de velocidades de propagação de frentes de material hiperconcentrado

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Modelo de propagação para redes sem fio fixas na banda de 5,8 GHZ em cidades típicas da região amazônica

O estudo da perda de propagação, nas cidades da região amazônica, envolve ambiente caracterizado pelo clima tropical e, suburbano densamente arborizado. Levando consideração à importância da faixa ISM 5,8 GHz, esta dissertação apresenta um modelo propagação para a faixa de frequência em questão, agregando as características da atenuação experimentada pela onda de rádio quando se propaga em ambientes de cidades típicas região amazônica. Para tanto, medidas de potência recebida foram coletadas em 335 clientes fixos, distribuídos em 12 cidades na região norte do Brasil, sendo estes atendidos pelo programa de inclusão digital do estado do Pará, Navega Pará. Também foram realizadas medidas com mobilidade no campus da Universidade Federal do Pará (UFPA). Apresenta ainda o desempenho do modelo proposto sobre outros modelos (Modelo SUI e COST231-Hata) descritos na literatura, para redes sem fio fixas e com mobilidade. As métricas desempenho utilizadas foram o erro RMS e o desvio padrão com relação aos dados medidos. O ajuste dos parâmetros do modelo proposto é realizado através do método de mínimos quadrados lineares, aplicado em duas etapas para diminuir a incerteza sobre os parâmetros ajustados. O modelo proposto alcançou um erro RMS de 3,8 dB e desvio padrão de 2,3 dB, superando os demais modelos que obtiveram erros RMS acima de 10 dB e desvios padrão acima de 5 dB. Os resultados obtidos mostram a sua eficiência sobre outros modelos para predição de perdas na faixa de 5,8 GHz em sistemas fixos e móveis.

Descargas parciais em barras estatóricas de hidrogeradores: modelagem por guia de onda e análise de um acoplador

Neste trabalho, tem-se como objetivo modelar as barras estatóricas de hidrogeradores por guias de onda multicamadas e estudar a aplicabilidade de um acoplador (sensor) para monitoramento das ocorrências de descargas parciais no interior dessas barras (mais precisamente, na sua camada isolante). Para isso, inicia-se por uma abordagem conceitual do fenômeno de descargas parciais, na qual é mostrada a importância do monitoramento das mesmas para identificar os sintomas de fragilidade dielétrica dos materiais isolantes dessas barras. Posteriormente, realiza-se uma apresentação sobre os métodos de monitoramento das descargas parciais, analisam-se os sinais de pulsos elétricos produzidos por elas e seus parâmetros típicos. Em seguida, estudam-se os modos de propagação de ondas eletromagnéticas produzidas pelas descargas parciais nas barras estatóricas de hidrogeradores, modelando-as por meio de guias de onda multicamadas. Para finalizar, apresenta-se um modelo de acoplador eletromagnético ou sensor, formado por uma linha de microfita, enfocando a sua aplicabilidade no sistema de monitoramento dessas ondas eletromagnéticas.

Migração por equação de onda em meios anisotrópicos com correção de amplitude

Este trabalho discute dois aspectos da migração em profundidade através da continuação para baixo dos campos de onda: o tratamento de modos evanescentes e a correção da amplitude dos eventos migrados. Estes dois aspectos são discutidos em meios isotrópicos e para uma classe de meios anisotrópicos. Migrações por diferenças finitas (FD) e por diferenças finitas e Fourier (FFD) podem ser instáveis em meios com forte variação lateral de velocidade. Estes métodos utilizam aproximações de Padé reais para representar o operador que descreve a propagação de ondas descendentes. Estas abordagens não são capazes de tratar corretamente os modos evanescentes, o que pode levar à instabilidades numéricas em meios com forte variação lateral de velocidade. Uma solução possível para esse problema é utilizar aproximação de Padé complexa, que consegue melhor representar os modos evanescentes associados às reflexões pós-críticas, e neste trabalho esta aproximação é utilizada para obter algoritmos FD e híbrido FD/FFD estáveis para migração em meios transversalmente isotrópicos com eixo de simetria vertical (VTI), mesmo na presença de forte variação nas propriedades elásticas do meio. A estabilidade dos algoritmos propostos para meios VTI foi validada através da resposta ao impulso do operador de migração e pela sua aplicação na migração de dados sintéticos, em meios fortemente heterogêneos. Métodos de migração por equação de onda em meios heterogêneos não tratam corretamente a amplitude dos eventos durante a propagação. As equações de onda unidirecionais tradicionais descrevem corretamente apenas a parte cinemática da propagação do campo de onda. Assim, para uma descrição correta das amplitudes deve-se usar as equações de onda unidirecionais de amplitude verdadeira. Em meios verticalmente heterogêneos, as equações de onda unidirecionais de amplitude verdadeira podem ser resolvidas analiticamente. Em meios lateralmente heterogêneos, essas equações não possuem uma solução analítica. Mesmo soluções numéricas tendem a ser instáveis. Para melhorar a compensação de amplitude na migração, em meios com variação lateral de velocidade, é proposto uma aproximação estável para solução da equação de onda unidirecional de amplitude verdadeira. Esta nova aproximação é implementada nas migrações split-step e diferenças finitas e Fourier (FFD). O algoritmo split-step com correção de amplitude foi estendido para meios VTI. A migração pré e pós-empilhamento de dados sintéticos, em meios isotrópicos e anisotrópicos, confirmam o melhor tratamento das amplitudes e estabilidade dos algoritmos propostos.

Estudo de propagação de ondas em planície de inundação para elaboração de plano de ação emergencial de barragens: UHE Três irmãos estudo de caso

Pós-graduação em Engenharia Civil - FEIS

Influência da espessura do osso cortical sobre a velocidade de propagação do ultrassom

OBJETIVO: Avaliar a influência da espessura do osso cortical sobre a velocidade de propagação do ultrassom (in vitro). MÉTODO: Foram utilizadas 60 lâminas ósseas confeccionadas a partir do fêmur de bovinos, com diferentes espessuras, variando de 1 a 6mm (10 de cada). As medidas da velocidade do ultrassom foram realizadas por aparelho projetado para este fim, utilizando técnica subaquática e por contato direto com auxílio de gel de acoplamento. Os transdutores foram posicionados de duas maneiras diferentes; opostos entre si, com o osso entre eles, sendo a medida chamada de transversal; e, paralelos na mesma superfície cortical, sendo a medida chamada de axial. RESULTADOS: Com o modo de transmissão axial, a velocidade de propagação do ultrassom aumenta conforme a espessura do osso cortical aumenta, independente da distância entre os transdutores, até a espessura de 5mm, mantendo-se constante após. Não houve alteração da velocidade quando o modo de transmissão foi transversal. CONCLUSÃO: A velocidade de propagação do ultrassom aumenta com o aumento da espessura da cortical óssea, no modo de transmissão axial, até o momento em que a espessura supera o comprimento da onda, mantendo a velocidade constante a partir de então. Nível de Evidência: Estudo Experimental.

Produção e caracterização de guias de onda dopados com terras-raras contendo nanopartículas semicondutoras.

O presente trabalho tem como objetivo estudar a produção e caracterização de filmes finos do tipo GeO2-Bi2O3 (BGO) produzidos por sputtering-RF com e sem nanopartículas (NPs) semicondutoras, dopados e codopados com íons de Er3+ ou Er3+/Yb3+ para a produção de amplificadores ópticos. A produção de guias de onda do tipo pedestal baseados nos filmes BGO foi realizada a partir de litografia óptica seguida por processo de corrosão por plasma e deposição física a vapor. A incorporação dos íons de terras-raras (TRs) foi verificada a partir dos espectros de emissão. Análises de espectroscopia e microscopia foram indispensáveis para otimizar os parâmetros dos processos para a construção dos guias de onda. Foi observado aumento significativo da luminescência do Er3+ (região do visível e do infravermelho), em filmes finos codopados com Er3+/Yb3+ na presença de nanopartículas de Si. As perdas por propagação mínimas observadas foram de ~1,75 dB/cm para os guias pedestal em 1068 nm. Para os guias dopados com Er3+ foi observado aumento significativo do ganho na presença de NPs de silício (1,8 dB/cm). O ganho óptico nos guias de onda amplificadores codopados com Er3+/Yb3+ e dopados com Er3+ com e sem NPs de silício também foi medido. Ganho de ~8dB/cm em 1542 nm, sob excitação em 980 nm, foi observado para os guias pedestal codopados com Er3+/Yb3+ (Er = 4,64.1019 átomos/cm3, Yb = 3,60.1020 átomos/cm3) com largura de 80 µm; para os guias codopados com concentração superior de Er3+/Yb3+ (Er = 1,34.1021 átomos/cm3, Yb = 3,90.1021 átomos/cm3) e com NPs de Si, foi observado aumento do ganho óptico de 50% para guia com largura de 100 µm. Os resultados apresentados demonstram que guias de onda baseados em germanatos, com ou sem NPs semicondutoras, são promissores para aplicações em dispositivos fotônicos.

Metodologias de determinação de alturas de onda para dimensionamento de obras marítimas

O clima marítimo ao longo da costa da Região Autónoma da Madeira (RAM) apresenta algumas diferenças, sendo para isso fundamental efetuar uma análise das características oceanográficas e geotécnicas do Arquipélago da Madeira. A referida análise tem como pressuposto retratar, o mais fiável possível, as diferenças, incidindo especialmente nos Portos do Funchal, Caniçal, Porto Moniz e Porto Santo. A caracterização e simulação da propagação da agitação marítima, aplicada a estruturas portuárias da RAM, representam o cerne desta dissertação, sendo os seus objetivos primordiais: a definição da onda de projeto de uma obra marítima; a descrição da metodologia de determinação de estimativas de longo prazo; a análise comparativa; a definição de metodologias de cálculo e a aplicação a registos da costa da RAM. Numa primeira fase procedeu-se à recolha de informação bibliográfica e de todos os registos informáticos disponibilizados pelas identidades responsáveis, nomeadamente a Administração dos Portos da Região Autónoma da Madeira (APRAM). Posteriormente efetuou-se um registo fotográfico das estruturas portuárias de modo a poder enquadra-las geograficamente na região. Com recurso aos dados fornecidos pela APRAM e a folhas de cálculo desenvolvidas no programa Microsoft Excel foi possível efetuar uma análise das alturas de onda significativa, do período de pico e da direção de ocorrência, por forma a poder caracterizar a agitação marítima que ocorre em cada um dos locais alvo deste estudo, sendo esta complementada pela análise das variáveis oceanográficas através do programa Carol. O passo seguinte incidiu sobre a análise da propagação da agitação marítima em cada um dos locais analisados, com recurso ao programa Mopla e à caracterização supracitada, por forma a conseguir definir a altura de onda que ocorre na estrutura portuária em causa. Como tarefa final são expostas as principais conclusões, algumas considerações finais e alguns objetivos que deveriam ser alcançados num futuro próximo.