Wave equation depth migration using complex Padé approximation

Propomos um novo método de migração em profundidade baseado na solução da equação da onda com densidade constante no domínio da freqüência. Uma aproximação de Padé complexa é usada para aproximar o operador de evolução aplicado na extrapolação do campo de ondas. Esse método reduz as imprecisões e instabilidades devido às ondas evanescentes e produz imagens com menos ruídos numéricos que aquelas obtidas usando-se a aproximação de Padé real para o operador exponencial, principalmente em meios com fortes variações de velocidades. Testes em dados de afastamento nulo do modelo de sal SEG/EAGE e nos dados de tiro comum 2-D Marmousi foram realizados. Os resultados obtidos mostram que o método de migração proposto consegue lidar com fortes variações laterais e também tem uma boa resposta para refletores com mergulhos íngremes. Os resultados foram comparados àqueles resultados obtidos com os métodos split-step Fourier (SSF), phase shift plus interpolarion (PSPI) e Fourier diferenças-finitas (FFD).

Modelagem da propagação sonora em dutos : abordagem matemática no domínio da frequência e do tempo utilizando transformada de laplace

Dissertação (mestrado)—Universidade de Brasília, Faculdade UnB Gama, Faculdade de Tecnologia, Programa de Pós-graduação em Integridade de Materiais da Engenharia, 2016.

Migração usando equação da onda unidirecional de abertura arbitrária

Neste trabalho avaliamos uma classe de operadores de continuação de campos de onda, baseados em equações unidirecionais e com aplicação direta à migração sísmica. O método de representação de equações de onda unidirecionais, desenvolvido neste trabalho, é válido para abertura angular arbitrária, baseia-se no conceito de rigidez de um semiespaço, na transformação Dirichlet-Neumann e em sua discretização por elementos finitos. O método de construção dos operadores de continuação requer a introdução de variáveis auxiliares cujo número cresce em função da maior abertura angular desejada para o operador. Efetuamos a implementação no domínio do espaço e da frequência o que permite sua imediata paralelização. Baseados em experimentos numéricos, que avaliam a relação de dispersão e a resposta ao impulso do operador, propomos prescrições que permitem especificar o número de variáveis auxiliares e o passo de continuação para o operador de migração. A aplicação do algoritmo nos dados do modelo de domo salino da SEG-EAGE demonstra a capacidade do algoritmo em migrar refletores com forte mergulho em meios com forte variação lateral de velocidade.

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 uma equação de onda não-linear

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

Controle exato para a equação de onda: o método HUM

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

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.

Soluções da equação de Fokker-Planck para um potencial isoespectral ao potencial de Morse

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

Solução da equação de Fokker-Planck para o potencial do oscilador harmônico duplo

Pós-graduação em Biofísica Molecular - IBILCE

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

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

Comportamento assintótico e controlabilidade exata para a equação de Klein-Gordon

Pós-graduação em Matemática - IBILCE

Equação de Fokker-Planck, supersimetria e enovelamento de proteína

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Soluções analíticas da equação de Duffin-Kemmer-Petiau

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

Migração FD e FFD com aproximações de grande abertura angular

Métodos de migração baseados na equação da onda unidirecional apresentam limitações no imageamento de refletores com forte mergulho e no tratamento de ondas evanescentes. Utilizando a expansão de Padé complexa na solução da equação da onda unidirecional para imageamento em geofísica é possível tratar ondas evanescentes, estabilizando a migração e dispensando um tratamento especial para absorção nas fronteiras do domínio. Utilizando várias parcelas na expansão de Padé, refletores com forte mergulho podem ser corretamente migrados. Nesta dissertação aplicamos a equação da onda unidirecional com expansão de Padé complexa para implementar dois algoritmos de migração em profundidade pré e pós empilhamento: a migração por diferenças finitas (FD) e a migração por diferenças finitas e Fourier (FFD). O estudo das curvas de dispersão e da resposta impulsiva dos operadores de migração nos permitiu escolher o número adequado de parcelas na expansão de Padé e os coeficientes da expansão que garantem estabilidade à continuação do campo de onda para um ângulo de mergulho máximo prescrito. As implementações foram validadas nos dados Marmousi e no modelo de domo de sal da SEG/EAGE mostrando que refletores com forte mergulho foram corretamente migrados, mesmo na presença de forte variação lateral de velocidade. Esses resultados são comparados com outros métodos de migração baseados na equação da onda unidirecional ressaltando a qualidade da aproximação estudada neste trabalho.