Migração Kirchhoff pré-empilhamento em profundidade modificada usando o operador de feixes gaussianos

A teoria dos feixes gaussianos foi introduzida na literatura sísmica no início dos anos 80 por pesquisadores russos e tchecos, e foi originalmente utilizada no cálculo do campo de ondas eletromagnéticas, baseado na teoria escalar da difração. Na teoria dos feixes gaussianos, o campo de ondas sísmicas é obtido por uma integral, cujo o integrando é constituído de duas partes, a saber: (1) as amplitudes dos campos das ondas na vizinhança do ponto de observação e (2) a função fase de cada um desses campos de ondas, que neste caso é representada por um tempo de trânsito paraxial complexo. Como ferramenta de imageamento, mais precisamente como operador de migração, os primeiros trabalhos usando feixes gaussianos datam do final da década de 80 e início dos anos 90. A regularidade dos campos de ondas descritos pelos feixes gaussianos, além de sua alta precisão em regiões singulares do modelo de velocidades, tornaram o uso de feixes gaussianos como uma alternativa híbrida viável para a migração. Nesse trabalho, unimos a flexibilidade da migração tipo Kirchhoff em profundidade em verdadeira amplitude com a regularidade da descrição do campo de ondas, representado pela sobreposição de feixes gaussianos. Como forma de controlar de forma estável quantidades usadas na construção de feixes gaussianos, utilizamos informações advindas do volume de Fresnel, mais precisamente a zona de Fresnel ao redor do ponto de reflexão e a zona de Fresnel projetada, localizada ao redor do ponto de registro do sismograma e cuja a informação se encontra nas curvas de reflexão de dados sísmico. Nosso processo de migração pode ser chamado como uma migração Kirchhoff em verdadeira amplitude usando um operador de feixes gaussianos.

Migração pré-empilhamento Kirchhoff feixes gaussianos 2,5D nos domínios afastamento comum e ângulo-comum

O método de migração do tipo Kirchhoff se apresenta na literatura como uma das ferramentas mais importantes de todo o processamento sísmico, servindo de base para a resolução de outros problemas de imageamento, devido ao um menor custo computacional em relação aos métodos que tem por base a solução numérica da equação da onda. No caso da aplicação em três dimensões (3D), mesmo a migração do tipo Kirchhoff torna-se dispendiosa, no que se refere aos requisitos computacionais e até mesmo numéricos para sua efetiva aplicação. Desta maneira, no presente trabalho, objetivando produzir resultados com uma razão sinal/ruído maior e um menor esforço computacional, foi utilizado uma simplificação do meio denominado 2.5D, baseado nos fundamentos teóricos da propagação de feixes gaussianos. Assim, tendo como base o operador integral com feixes gaussianos desenvolvido por Ferreira e Cruz (2009), foi derivado um novo operador integral de superposição de campos paraxiais (feixes gaussianos), o mesmo foi inserido no núcleo do operador integral de migração Kirchhoff convencional em verdadeira amplitude, para a situação 2,5D, definindo desta maneira um novo operador de migração do tipo Kirchhoff para a classe pré-empilhamento em verdadeira amplitude 2.5D (KGB,do inglês Kirchhoff-Gausian-Beam). Posteriormente, tal operador foi particularizado para as configurações de medida afastamento comum (CO, do inglês common offset) e ângulo de reflexão comum (CA, do inglês common angle), ressaltando ainda, que na presente Tese foi também idealizada uma espécie de flexibilização do operador integral de superposição de feixes gaussianos, no que concerne a sua aplicação em mais de um domínio, quais sejam, afastamento comum e fonte comum. Nesta Tese são feitas aplicações de dados sintéticos originados a partir de um modelo anticlinal.

Migração 3-D Kirchhoff-Gaussian-Beam (KGB) pré-empilhamento no domínio da profundidade

O Feixe Gaussiano (FG) é uma solução assintótica da equação da elastodinâmica na vizinhança paraxial de um raio central, a qual se aproxima melhor do campo de ondas do que a aproximação de ordem zero da Teoria do Raio. A regularidade do FG na descrição do campo de ondas, assim como a sua elevada precisão em algumas regiões singulares do meio de propagação, proporciona uma forte alternativa na solução de problemas de modelagem e imageamento sísmicos. Nesta Tese, apresenta-se um novo procedimento de migração sísmica pré-empilhamento em profundidade com amplitudes verdadeiras, que combina a flexibilidade da migração tipo Kirchhoff e a robustez da migração baseada na utilização de Feixes Gaussianos para a representação do campo de ondas. O algoritmo de migração proposto é constituído por dois processos de empilhamento: o primeiro é o empilhamento de feixes (“beam stack”) aplicado a subconjuntos de dados sísmicos multiplicados por uma função peso definida de modo que o operador de empilhamento tenha a mesma forma da integral de superposição de Feixes Gaussianos; o segundo empilhamento corresponde à migração Kirchhoff tendo como entrada os dados resultantes do primeiro empilhamento. Pelo exposto justifica-se a denominação migração Kirchhoff-Gaussian-Beam (KGB). As principais características que diferenciam a migração KGB, durante a realização do primeiro empilhamento, de outros métodos de migração que também utilizam a teoria dos Feixes Gaussianos, são o uso da primeira zona de Fresnel projetada para limitar a largura do feixe e a utilização, no empilhamento do feixe, de uma aproximação de segunda ordem do tempo de trânsito de reflexão. Como exemplos são apresentadas aplicações a dados sintéticos para modelos bidimensionais (2-D) e tridimensionais (3-D), correspondentes aos modelos Marmousi e domo de sal da SEG/EAGE, respectivamente.

Análise do efeito da discretização do modelo de velocidades nas migrações Kirchhoff e Kirchhoff-Gaussian- Beam 2D pré-empilhamento em profundidade

O Feixe Gaussiano (FG) é uma solução assintótica da equação da elastodinâmica na vizinhança paraxial de um raio central, a qual se aproxima melhor do campo de ondas do que a aproximação de ordem zero da Teoria do Raio. A regularidade do FG na descrição do campo de ondas, assim como a sua elevada precisão em algumas regiões singulares do meio de propagação, proporciona uma forte alternativa no imageamento sísmicos. Nesta dissertação, apresenta-se um novo procedimento de migração sísmica pré-empilhamento em profundidade com amplitudes verdadeiras, que combina a flexibilidade da migração tipo Kirchhoff e a robustez da migração baseada na utilização de Feixes Gaussianos para a representação do campo de ondas. O algoritmo de migração proposto é constituído por dois processos de empilhamento: o primeiro é o empilhamento de feixes (“beam stack”) aplicado a subconjuntos de dados sísmicos multiplicados por uma função peso definida de modo que o operador de empilhamento tenha a mesma forma da integral de superposição de Feixes Gaussianos; o segundo empilhamento corresponde à migração Kirchhoff tendo como entrada os dados resultantes do primeiro empilhamento. Pelo exposto justifica-se a denominação migração Kirchhoff-Gaussian-Beam (KGB).Afim de comparar os métodos Kirchhoff e KGB com respeito à sensibilidade em relação ao comprimento da discretização, aplicamos no conjunto de dados conhecido como Marmousi 2-D quatro grids de velocidade, ou seja, 60m, 80m 100m e 150m. Como resultado, temos que ambos os métodos apresentam uma imagem muito melhor para o menor intervalo de discretização da malha de velocidade. O espectro de amplitude das seções migradas nos fornece o conteúdo de frequência espacial das seções das imagens obtidas.

Antenas cornetas para o satélite CosmoGal

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Electrónica e Telecomunicações

Análise por feixes de íons de filmes finos dielétricos depositados por sputtering reativo e crescidos termicamente

Esta tese trata de filmes finos de materiais dielétricos depositados por sputtering reativo e crescidos termicamente que são analisados por métodos de feixes de íons. Como introdução aos nossos trabalhos, são abordados os princípios da deposição por sputtering com corrente contínua, com rádio-freqüências e em atmosferas reativas. Além disso, são apresentados os princípios de crescimento térmico em fornos clássicos e em fornos de tratamento térmico rápido. Também são descritos em detalhe os métodos análiticos de espectropia de retroespalhamento Rutheford e reações nucleares ressonantes e não-ressonantes para medir as quantidades totais e os perfis de concentração em função da profundidade dos vários isótopos utilizados, bem como o método de difração de raios-X. No que concerne à deposição por sputtering reativo de filmes finos de nitreto de silicio, nitreto de alumínio e nitreto duplo de titânio e alumínio, são estudados os processos de deposição, são caracterizados os filmes desses materiais e são estabelecidas correlações entre características dos filmes e os parâmetros de deposição. Quanto ao crescimento térmico em atmosferas reativas de filmes muito finos de óxido, nitreto e oxinitreto de sílicio, são apresentados modelos para o crescimento dos filmes de óxido de sílicio, é estudada a influência da limpeza do substrato de sílicio nos mecanismos de crescimento dos filmes de óxido de sílicio e são elucidados aspectos dos mecanismos de nitretação térmica do sílicio e de filmes de óxido de sílicio.

Estudo da instabilidade de feixes de condutores em linhas aéreas de transmissão

A demanda de energia elétrica tem aumentado consideravelmente com o observado crescimento da atividade industrial, exigindo o desenvolvimento de novas tecnologias para aumentar a capacidade de transporte de energia por linha. Neste sentido, as linhas de transmissão de alta voltagem, anteriormente formadas por condutores isolados, podem ter a sua capacidade expressivamente incrementada dispondo-se os condutores em feixes, ou seja, com mais de um condutor por fase. A ação do vento nos condutores em forma de feixes originam problemas vibratórios diferentes daqueles comumente encontrados em condutores isolados. Além de oscilações por galope e desprendimento de vórtices, feixes são suscetíveis a oscilações de baixa freqüência causada pelo efeito de interferência das esteiras entre os condutores. Estas oscilações são de dois tipos: (i) movimento de sub-vão, que é característico do condutor a sotavento, o qual pode-se deslocar independentemente dos outros condutores vizinhos, e (ii) movimento de vão completo, no qual movimentos do feixe ocorrem no vão compreendido entre as torres de suporte. Alguns trabalhos na literatura consideram os movimentos de vão completo em condutores coberto por gelo, mas pouca informação se tem sobre o estudo de oscilações em feixes sem a presença de gelo. Devido a complexidade destes fenômenos, não há ainda critérios claros quanto à estabilidade do feixe para as configurações usuais de linhas deste tipo. Este trabalho tem então como objetivo inicial apresentar as principais características dinâmicas de condutores isolados e em feixes. Além disso, apresentar uma análise dos parâmetros que influenciam os movimentos, dentre os quais são citados: características do terreno, velocidade e turbulência do vento, número e arranjo dos condutores, rugosidade do condutor, espaçamento entre condutores, inclinação do feixe, sistemas de espaçadores e de suspensão, tração do condutor, freqüências naturais e efeito de flecha. Ensaios em túnel de vento são normalmente realizados com o propósito de obter-se dados de tipos específicos de condutores e configurações de feixes de condutores. Posteriormente, estes dados podem ser utilizados na análise teórica dos movimentos dos condutores a fim de prever-se a estabilidade dos mesmos. Neste trabalho foram projetados experimentos do tipo estático, que permitiram a obtenção dos coeficientes aerodinâmicos e suas derivadas em relação aos ângulos de incidência do vento em feixes de cabos. Para isto utilizaram-se duas células de carga com extensômetros, onde em uma delas mediu-se diretamente as forças de arrasto e de sustentação e na outra os momentos de torção. Os experimentos foram conduzidos no túnel de vento Prof. Joaquim Blessmann, da UFRGS. Os modelos eram compostos por condutores rígidos de pequeno comprimento que estavam dispostos isoladamente, em feixe de 02 cabos dispostos lado-a-lado ou em feixes de 04 cabos dispostos na configuração de um quadrado. Os ensaios foram realizados para diferentes espaçamentos entre os cabos (10, 14, 18 e 22 diâmetros), para diversos ângulos de inclinação do feixe (0o a 45o) e ainda para diferentes velocidades e condições de escoamento (suave ou turbulento). Os ensaios foram realizados com cabos lisos e cabos do tipo Rook ACSR 24/7, a fim de analisar-se a influência da rugosidade no comportamento dos modelos. Adicionalmente, foi desenvolvida uma aplicação da análise do problema de instabilidade dinâmica através da utilização de equações linearizadas do movimento de condutores em feixe. A partir destas equações, determinam-se as regiões de instabilidade das oscilações de vão completo, em feixes de dois e quatro condutores. Os coeficientes aerodinâmicos utilizados nestas equações foram aqueles determinados nos ensaios no túnel de vento Prof. Blessmann. Finalmente, o conhecimento da influência dos parâmetros vinculados as características do feixe, vão e vento incidente nos coeficientes aerodinâmicos e na conseqüente estabilidade do feixe, possibilita a determinação de alguns critérios de projeto que garantam maior estabilidade dos feixes de condutores.

Estudo do crescimento de filmes de carbono sobre silício devido à irradiação com feixes de H e He

É bem conhecido que técnicas experimentais que fazem uso de feixes iônicos induzem sobre a superfície dos alvos irradiados o depósito de diversos tipos de impurezas. Este depósito é causado pelas interações entre o feixe de íons, as moléculas de gás residual no interior das câmaras de irradiação e a superfície da amostra. Apesar do fato deste processo poder alterar significativamente os parâmetros de irradiações, bem como afetar a análise de materiais tratados, os parâmetros experimentais que influenciam na deposição das impurezas ainda não são bem conhecidos e nem o depósito se encontra suficientemente quantificado. No presente trabalho relatamos um estudo quantitativo da deposição de carbono sobre amostras de Si (100) irradiadas com feixes de He e H. A deposição de carbono foi determinada em função da fluência de irradiação, variando diversos parâmetros experimentais, tais como: pressão na câmara de irradiação, temperatura do alvo, densidade de corrente do feixe, energia do feixe e o estado da carga do íon. Em todos os casos a análise das amostras irradiadas foi feita pela técnica de Retroespalhamento de Rutherford em direção Canalizada (RBS/C) e através da reação nuclear ressonante 12C(a, a´)12C na energia de 4265 keV. Os resultados experimentais mostram que se consegue minimizar a deposição de C através: a) da redução do tempo de irradiação, b) da redução da pressão na câmara de irradiação, c) do aumento da temperatura do alvo durante a irradiação e d) minimização do poder de freamento do íon no alvo.

Determinação da resposta dinâmica de feixes de condutores em linhas aéreas de transmissão

Até os dias de hoje, a abordagem estática tem sido a usualmente empregada para a avaliação da resposta de cabos de Linhas de Transmissão (LTs), apesar da resposta dinâmica ser, em muitos casos, reconhecidamente importante na avaliação do desempenho dos cabos. Para uma análise completa, não se pode desconhecer a natureza dinâmica da maioria dos fenômenos a que as LTs estão submetidas, sendo exemplos típicos as excitações mecânicas causadas pela ação do vento (carregamentos de baixa e alta freqüência) e a ruptura de cabos. Cabe salientar que mesmo a análise estática empregada tem sido simplificada, considerando normalmente apenas casos contemplados com soluções analíticas. A justificativa para emprego da análise estática baseia-se no menor esforço numérico exigido. Entretanto, essa justificativa já não se sustenta dado os grandes avanços na área computacional, possibilitando o estudo do desempenho de LTs na ocorrência de fenômenos que provocam carregamentos dinâmicos. Adicionalmente, há uma crescente demanda para que se obtenha um melhor entendimento de muitas questões relativas ao comportamento dinâmico das LTs, que possam explicar desempenhos observados e/ou sustentar o desenvolvimento de novas alternativas mais arrojadas de projeto e construção Entre os diversos métodos utilizados na engenharia estrutural, a integração direta das equações do movimento, através de métodos numéricos como as diferenças finitas centrais, se constitui numa poderosa ferramenta de cálculo. Esta ferramenta possibilita o tratamento de problemas envolvendo não linearidades geométricas e do material, como são os casos onde se avalia a resposta de cabos suspensos submetidos a carregamentos variáveis no tempo. Esta pesquisa objetiva a aplicação do método da integração direta das equações do movimento na análise de feixes de cabos de LT, quando submetidos à ação de carregamentos mecânicos variáveis no tempo, principalmente à excitação de ventos oriundos de fenômenos com natureza complexa (tormentas elétricas, por exemplo). É apresentado um método para determinar a resposta dinâmica de feixes que considera a interação entre o vento incidente e o movimento do condutor. A solução é obtida por integração numérica, no domínio do tempo, das equações de movimento de um modelo tridimensional não-linear discreto de feixe, as quais definem as forças nos espaçadores e cabos através de coeficientes aerodinâmicos obtidos experimentalmente. Também são apresentados modelos dos fenômenos meteorológicos mais comuns em nosso país (Brasil): tormenta extratropical (EPS) e tormenta elétrica (TS) Como ilustração, são apresentados exemplos de modelagem de vãos de LTs com condutor singelo e com feixes de condutores. Os exemplos demonstram a capacidade de avaliação de Estados Limites relacionados à distância relativa entre subcondutores, à estabilidade do feixe, à representação das suas propriedades e ao comportamento dinâmico, bem como aos carregamentos transmitidos às estruturas. A análise emprega conhecidas relações constitutivas para representar o comportamento tensão-deformação dos cabos. O enfoque utilizado possibilita a avaliação mais precisa de casos reais que ainda não podiam ser convenientemente tratados, além de permitir a extensão para estudos bem mais complexos, tais como feixes com disposições assimétricas.

Efeitos não-gaussianos em astrofísica e cosmologia

The recent astronomical observations indicate that the universe has null spatial curvature, is accelerating and its matter-energy content is composed by circa 30% of matter (baryons + dark matter) and 70% of dark energy, a relativistic component with negative pressure. However, in order to built more realistic models it is necessary to consider the evolution of small density perturbations for explaining the richness of observed structures in the scale of galaxies and clusters of galaxies. The structure formation process was pioneering described by Press and Schechter (PS) in 1974, by means of the galaxy cluster mass function. The PS formalism establishes a Gaussian distribution for the primordial density perturbation field. Besides a serious normalization problem, such an approach does not explain the recent cluster X-ray data, and it is also in disagreement with the most up-to-date computational simulations. In this thesis, we discuss several applications of the nonextensive q-statistics (non-Gaussian), proposed in 1988 by C. Tsallis, with special emphasis in the cosmological process of the large structure formation. Initially, we investigate the statistics of the primordial fluctuation field of the density contrast, since the most recent data from the Wilkinson Microwave Anisotropy Probe (WMAP) indicates a deviation from gaussianity. We assume that such deviations may be described by the nonextensive statistics, because it reduces to the Gaussian distribution in the limit of the free parameter q = 1, thereby allowing a direct comparison with the standard theory. We study its application for a galaxy cluster catalog based on the ROSAT All-Sky Survey (hereafter HIFLUGCS). We conclude that the standard Gaussian model applied to HIFLUGCS does not agree with the most recent data independently obtained by WMAP. Using the nonextensive statistics, we obtain values much more aligned with WMAP results. We also demonstrate that the Burr distribution corrects the normalization problem. The cluster mass function formalism was also investigated in the presence of the dark energy. In this case, constraints over several cosmic parameters was also obtained. The nonextensive statistics was implemented yet in 2 distinct problems: (i) the plasma probe and (ii) in the Bremsstrahlung radiation description (the primary radiation from X-ray clusters); a problem of considerable interest in astrophysics. In another line of development, by using supernova data and the gas mass fraction from galaxy clusters, we discuss a redshift variation of the equation of state parameter, by considering two distinct expansions. An interesting aspect of this work is that the results do not need a prior in the mass parameter, as usually occurs in analyzes involving only supernovae data.Finally, we obtain a new estimate of the Hubble parameter, through a joint analysis involving the Sunyaev-Zeldovich effect (SZE), the X-ray data from galaxy clusters and the baryon acoustic oscillations. We show that the degeneracy of the observational data with respect to the mass parameter is broken when the signature of the baryon acoustic oscillations as given by the Sloan Digital Sky Survey (SDSS) catalog is considered. Our analysis, based on the SZE/X-ray data for a sample of 25 galaxy clusters with triaxial morphology, yields a Hubble parameter in good agreement with the independent studies, provided by the Hubble Space Telescope project and the recent estimates of the WMAP

Avaliação por meio da tomografia computadorizada de feixes cônicos e radiografia cefalométrica lateral do espaço aéreo faríngeo em indivíduos submetidos à cirurgia ortognática

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

Estudo da interação de prótons com alvos não homogêneos, aplicados a tomografia com feixes de prótons

Pós-graduação em Biologia Geral e Aplicada - IBB

Análise da viabilidade de simulações com feixes de prótons de 150MeV em alvos heterogêneos

A protonterapia é um tipo de tratamento promissor quanto a diminuição de dose nos tecidos adjacentes ao tumor, devido às suas característica de deposição de energia ao interagir com a matéria. Para aumentar a acurácia do tratamento é necessário que as simulações sejam feitas a partir de um sistema de formação de imagem que utilize também um feixe de prótons, como a tomografia por feixe de prótons (pCT). Essa área atualmente está em desenvolvimento e se defronta com inúmeros problemas na formação de uma imagem precisa. Este trabalho tem o intuito de analisar a viabilidade de simulações computacionais para este tipo de feixe. O código computacional SRIM é um conjunto de programas que calcula o freamento e o alcance de íons ao interagir com a matéria. Esse cálculo baseia-se em ferramentas oriundas da mecânica quântica e também no método de Monte Carlo, para determinar, de modo probabilístico, trajetórias, colisões e ângulos. Foram realizadas simulações no programa SRIM para um alvo heterogêneo envolvendo água e osso cortical. Para analisar o comportamento do feixe transmitido, foram feitas alterações na posição e no tamanho da heterogeneidade. A partir dos dados obtidos pode-se obter valores médios, histogramas e gaussianas dos mesmos para cada simulação e tirar conclusões sobre a viabilidade e as dependências desse tipo de simulação

Uso do método de Monte Carlo para simulação e detecção do espectro de raios-X gerado com feixes de elétrons em alvo de Molibdênio, utilizando o Geant4

The huge demand for procedures involving ionizing radiation promotes the need for safe methods of experimentation considering the danger of their biological e ects with consequent risk to humans. Brazilian's legislation prohibits experiments involving this type of radiation in humans through Decree 453 of Ministry of Health with determines that such procedures comply with the principles of justi cation, optimization and dose limitation. In this line, concurrently with the advancement of available computer processing power, computing simulations have become relevant in those situations where experimental procedures are too cost or impractical. The Monte Carlo method, created along the Manhattan Project duringWorldWar II, is a powerful strategy to simulations in computational physics. In medical physics, this technique has been extensively used with applications in diagnostics and cancer treatment. The objective of this work is to simulate the production and detection of X-rays for the energy range of diagnostic radiology, for molybdenum target, using the Geant4 toolkit. X-ray tubes with this kind of target material are used in diagnostic radiology, speci cally in mammography, one of the most used techniques for screening of breast cancer in women. During the simulations, we used di erent models for bremsstrahlung available in physical models for low energy, in situations already covered by the literature in earlier versions of Geant4. Our results show that although the physical situations seems qualitatively adequate, quantitative comparisons to available analytical data shows aws in the code of Geant4 Low Energy source

Sistema de cálculo computacional da unidade monitora para feixes de elétrons em radioterapia

In radiation theraphy with electron beam, the electrons are produced in linear accelerators, and energy the most used have between 4MeV and 20MeV. Generally, the treatments are done for superficial injuries, because the low penetration of these particles. In this work a system for calculation of monitor units (U.M.) for cases of treatments with electron beam was developed. The Excel program of Microsoft was used and is easily found in the operational system of the personal microcomputers. In the Excel has been inserted the pertinent data of the linear accelerator of Varian, model 2100C, used in the Service of radiation theraphy of the Hospital of the Clinics of the College of Medicine of the UNESP of Botucatu. For some values of the physical parameters, such as: factors field and factors calibration, not supplied in the tests of acceptance of the machine, still proceeded calculations from interpolation and extrapolation. The mathematical formulas for automatic search of these and others factors used in the calculations of the determination of the U.M had been developed in agreement available routines in Excel. For this the functions had been used the function IF (that it imposes search condition) and the PROCH (that looks a value in a column from determined line), beyond the basic functions of addition, multiplication and division. It is intended to optimize the routine of the Services of radiation theraphy that perform through eletrontheraphy procedures, speeding the calculations and minimizing the occurrence of errors and uncertainties deriving of the maken a mistake manipulation of the parameters gotten in tables of data of electron beams