Marcação e atenuação de múltiplas de superfície livre, processamento e imageamento em dados sísmicos marinhos

O presente trabalho teve como objetivos a análise e atenuação de múltiplas de superfície livre, o processamento e o imageamento de dados sísmicos marinhos visando obter imagens migradas de utilidade à interpretação geológica voltada à exploração de petróleo. Foi dada uma atenção sistemática ao estudo de múltiplas de superfície livre do ponto de vista do filtro de predição baseado na teoria da comunicação com o objetivo de melhor aplicar o filtro WH de deconvolução preditiva, na etapa de processamento, posterior à correção NMO, embora outros métodos possam ser considerados mais competitivos. A identificação e a atenuação de reflexões múltiplas em dados sísmicos reais continua como um grande desafio no processamento de dados sísmicos, uma vez que elas são consideradas como ruído. No entanto, sendo este ruído classificado como coerente, várias técnicas foram desenvolvidas objetivando sua atenuação para evitar erros em cascata nas etapas posteriores como processamento, marcação de eventos, inversão tomográfica, imageamento, e finalmente na interpretação geológica das imagens obtidas. Outro aspecto do trabalho foi estabelecer um fluxograma de processamento-imageamento sendo a atenuação de múltiplas de superfície livre uma etapa central. Foram obtidas seções migradas em tempo e em profundidade onde se permite a interpretação que se deseja. O desenvolvimento deste trabalho foi realizado com os pacotes CWP/SU e MatLab.

Modelagem acústica por diferenças finitas e elementos finitos em 2-D e 2,5-D

A modelagem acústica fornece dados úteis para avaliação de metodologias de processamento e imageamento sísmico, em modelos com estrutura geológica complexa. Esquemas de diferenças finitas (DF) e elementos finitos (EF) foram implementados e avaliados em modelos homogêneos e heterogêneos. O algoritmo de diferenças finitas foi estendido para o caso 2,5-D em modelos com densidade variável. Foi apresentada a modelagem de alvos geológicos de interesse exploratório existentes na Bacia Paleozóica do Solimões na Amazônia. Reflexões múltiplas de longo período produzidas entre a superfície livre e a discordância Cretáceo-Paleozóica, a baixa resolução da onda sísmica nas proximidades do reservatório e as fracas reflexões na interface entre as rochas reservatório e as rochas selantes são as principais características dos dados sintéticos obtidos, os quais representam um grande desafio ao imageamento sísmico.

Refinamento da representação de raízes no modelo de biosfera SiB2 em área de floresta na Amazônia

O objetivo desta pesquisa foi estudar a sensibilidade dos fluxos de superfície e hidrologia do solo em relação à representação e distribuição de raízes no modelo de biosfera para uma floresta de terra firme na Amazônia. A finalidade foi avaliar o impacto na representatividade dos fluxos de energia considerando a sazonalidade da região amazônica, usando como suporte medidas intensivas realizadas em uma reserva biológica (Reserva Biológica do Cuieiras, em Manaus). Foram realizadas oito simulações com o modelo de biosfera SiB2 (“Simple Biosphere Model” – versão 2) , onde cada simulação representou um cenário diferente de distribuição de raízes em uma profundidade de 4 m de solo, dividido em três camadas: 0,5 m, 1,5 m e 2,0 m. As raízes foram distribuídas privilegiando a concentração de raízes na camada superficial, em seguida, na camada intermediária e, por fim, uma concentração de raízes abaixo de 2,0 m de profundidade. As simulações foram realizadas para o período de 2003 a 2006, enfatizando o ano de 2005 para avaliar o efeito da representação de raízes nos fluxos de energia (calor latente – LE e calor sensível – H) e de dióxido de carbono (CO₂). A partir da análise integrada dos fluxos simulados com dados observacionais medidos no sítio experimental foi possível perceber que uma redução na precipitação no ano de 2005, apesar de ter sido menor na parte central da Amazônia, implicou na diminuição da umidade do solo, mostrando que a floresta passou por um período de estresse hídrico maior do que os outros anos analisados. O modelo representou a energia disponível com valores muito próximos aos observados, variando sazonalmente em concordância com os dados medidos em 2005. No entanto, LE é superestimado durante a estação chuvosa, mas mostra juntamente com o fluxo de CO₂, a redução com a umidade do solo na estação seca, enquanto H é superestimado em até 20 W m^-2 durante todo o período simulado. Estes resultados mostram que, a consideração de raízes rasas é mais apropriada para regiões que possuem estação seca curta, conforme se caracteriza a área de estudo, e raízes profundas devem favorecer a modelagem dos processos para superfície de áreas com estação seca mais pronunciada. Com isto os resultados revelam que há necessidade de obter mais informações de propriedades físicas do solo, apropriadas às condições da região, para que outros refinamentos sejam efetivos na distinção do comportamento de florestas tropicais sob diferentes regimes de disponibilidade de água no solo.

Avaliação comparativa das perdas de água por evapotranspiração em mesocosmos colonizados por diferentes macrófitas aquáticas

Water loss increase through evapotranspiration is one the most important modes of interference of weeds in water bodies, especially in water storage reservoirs. Thus, this research was conducted to compare the water losses in mesocosms colonized with the aquatic weeds Myriophyllum aquaticum,Brachiaria subquadripara,Echinochloa polystachya,Typha latifolia and Pontederia lanceolata. The evaluations were carried out in July 2004, in three 4-day periods. The water losses were evaluated based on the amounts of water necessary to restore the existing water level before each evaluation period. The colonization of the aquatic weeds increased the water losses in the mesocosms, with the highest losses being observed in those colonized by T. latifolia, between 3.54 to 4.71 times the water surface without aquatic weeds. The losses in the mesocosms colonized by M. aquaticum,B. subquadripara, E. polystachya, and P. lanceolata were statistically similar and promoted increases between 1.54 to 2.21 times the free surface. The results showed that aquatic weed control is important to prevent water losses in reservoirs used for water storage.

Desenvolvimento de software utilizando a técnica sph (smoothed particle hydrodynamics) na geração de ondas de submersão

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

Redes neurais artificiais artmap-fuzzy aplicadas ao estudo de agitação marítima e ondas de lagos

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

A Numerical Study of the FENE-CR Model Applied to a Jet Flow Problem

The FENE-CR model is investigated through a numerical algorithm to simulate the time-dependent moving free surface flow produced by a jet impinging on a flat surface. The objective is to demonstrate that by increasing the extensibility parameter L, the numerical solutions converge to the solutions obtained with the Oldroyd-B model. The governing equations are solved by an established free surface flow solver based on the finite difference and marker-and-cell methods. Numerical predictions of the extensional viscosity obtained with several values of the parameter L are presented. The results show that if the extensibility parameter L is sufficiently large then the extensional viscosities obtained with the FENE-CR model approximate the corresponding Oldroyd-B viscosity. Moreover, the flow from a jet impinging on a flat surface is simulated with various values of the extensibility parameter L and the fluid flow visualizations display convergence to the Oldroyd-B jet flow results.

Investigation of the structural, optical and dielectric properties of highly (100)-oriented (Pb0.60Ca0.20Sr0.20)TiO3 thin films on LaNiO3 bottom electrode

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

Diodo planar orgânico Schottky: construção e caracterização

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

Effects of Ciprofloxacin-containing Scaffolds on Enterococcus faecalis Biofilms

Introduction: Antibiotic-containing polymer-based nanofibers (hereafter referred to as scaffolds) have demonstrated great potential for their use in regenerative endodontics from both an antimicrobial and cytocompatibility perspective. This study sought to evaluate in vitro the effects of ciprofloxacin (CIP)-containing polymer scaffolds against Enterococcus faecalis biofilms. Methods: Human mandibular incisors were longitudinally sectioned to prepare radicular dentin specimens. Sterile dentin specimens were distributed in 24-well plates and inoculated with E. faecalisfor biofilm formation. Infected dentin specimens were exposed to 3 groups of scaffolds, namely polydioxanone (PDS) (control), PDS + 5 wt% CIP, and PDS + 25 wt% CIP for 2 days. Colony-forming units (CFU/mL) (n = 10) and scanning electron microscopy (SEM) (n -= 2) were performed to quantitatively and qualitatively assess the antimicrobial effectiveness, respectively. Results: PDS scaffold containing CIP at 25 wt% showed maximum bacteria elimination with no microbial growth, differing statistically (P < .05) from the control (PDS) and from PDS scaffold containing CIP at 5 wt%. Statistical differences (P < .05) were also seen for the CFU/mL data between pure PDS (5.92-6.02 log CFU/mL) and the PDS scaffold containing CIP at 5 wt% (5.39 5.87 log CFU/mL). SEM images revealed a greater concentration of bacteria on the middle third of the dentin specimen. after 5 days of biofilm formation. On scaffold exposures, SEM images showed similar results when compared with the CFU/mL data. Dentin specimens exposed to PDS + 25 wt% CIP scaffolds displayed a practically bacteria-free surface. Conclusions: On the basis of the data presented, newly developed antibiotic-containing electrospun scaffolds hold promise as an intracanal medicament to eliminate biofilm/infection before regenerative procedures.

Magnetic-size effects in ultrathin layers

We study N-layer samples (N ≤ 10) for the Heisenberg model, with ferro- and antiferromagnetic exchange couplings, using a modified version of the Onsager reaction field approximation. The present scheme includes short-range spin-spin correlations, and allows for layer-dependent order parameters when free surface boundary conditions are imposed. The limits N = 1 (two dimensions) and N → ∞ (three dimensions) can be solved analytically, while systems with several layers have to be numerically calculated. We found no indication of a phase transition at finite temperature up to the sizes investigated (N = 10), the layered systems behaving essentially as two-dimensional. A phase transition is only obtained for the three-dimensional limit. © 1993.

NIR luminescent Er3+/Yb3+ co-doped SiO2-ZrO2 nanostructured planar and channel waveguides: Optical and structural properties

Optical and structural properties of planar and channel waveguides based on sol gel Er3+ and Yb3+ co-doped SiO2-ZrO2 are reported. Microstructured channels with high homogeneous surface profile were written onto the surface of multilayered densified films deposited on SiO2/Si substrates by a femtosecond laser etching technique. The densification of the planar waveguides was evaluated from changes in the refractive index and thickness, with full densification being achieved at 900 degrees C after annealing from 23 up to 500 min, depending on the ZrO2 content Crystal nucleation and growth took place together with densification, thereby producing transparent glass ceramic planar waveguides containing rare earth-doped ZrO2 nanocrystals dispersed in a silica-based glassy host Low roughness and crack-free surface as well as high confinement coefficient were achieved for all the compositions. Enhanced NIR luminescence of the Er3+ ions was observed for the Yb3+- codoped planar waveguides, denoting an efficient energy transfer from the Yb3+ to the Er3+ ion. (C) 2012 Elsevier B.V. All rights reserved.

Dependence of brittle-to-ductile transition on crystallographic direction in diamond turning of single-crystal silicon

The objective of this paper is to show the dependence relationship between the crystallographic orientations upon brittle-to-ductile transition during diamond turning of monocrystalline silicon. Cutting tests were performed using a -5 degrees rake angle round nose diamond tool at different machining scales. At the micrometre level, the feedrate was kept constant at 2.5 micrometres per revolution (mu m/r), and the depth of cut was varied from 1 to 5 mu m. At the submicrometre level, the depth of cut was kept constant at 500 nm and the feedrate varied from 5 to 10 mu m/r. At the micrometre level, the uncut shoulder generated with an interrupted cutting test procedure provided a quantitative measurement of the ductile-to-brittle transition. Results show that the critical chip thickness in silicon for ductile material removal reaches a maximum of 285 nm in the [100] direction and a minimum of 115 nm in the [110] direction, when the depth of cut was 5 mu m. It was found that when a submicrometre depth of cut was applied, microcracks were revealed in the [110] direction, which is the softer direction in silicon. Micro Raman spectroscopy was used to estimate surface residual stress after machining. Compressive residual stress in the range 142 MPa and smooth damage free surface finish was probed in the [100] direction for a depth of cut of 5 mu m, whereas residual stresses in the range 350 MPa and brittle damage was probed in the [110] direction for a depth of cut of 500 nm.

Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.

Two fluid numerical model for coastal applications

Wave breaking is an important coastal process, influencing hydro-morphodynamic processes such as turbulence generation and wave energy dissipation, run-up on the beach and overtopping of coastal defence structures. During breaking, waves are complex mixtures of air and water (“white water”) whose properties affect velocity and pressure fields in the vicinity of the free surface and, depending on the breaker characteristics, different mechanisms for air entrainment are usually observed. Several laboratory experiments have been performed to investigate the role of air bubbles in the wave breaking process (Chanson & Cummings, 1994, among others) and in wave loading on vertical wall (Oumeraci et al., 2001; Peregrine et al., 2006, among others), showing that the air phase is not negligible since the turbulent energy dissipation involves air-water mixture. The recent advancement of numerical models has given valuable insights in the knowledge of wave transformation and interaction with coastal structures. Among these models, some solve the RANS equations coupled with a free-surface tracking algorithm and describe velocity, pressure, turbulence and vorticity fields (Lara et al. 2006 a-b, Clementi et al., 2007). The single-phase numerical model, in which the constitutive equations are solved only for the liquid phase, neglects effects induced by air movement and trapped air bubbles in water. Numerical approximations at the free surface may induce errors in predicting breaking point and wave height and moreover, entrapped air bubbles and water splash in air are not properly represented. The aim of the present thesis is to develop a new two-phase model called COBRAS2 (stands for Cornell Breaking waves And Structures 2 phases), that is the enhancement of the single-phase code COBRAS0, originally developed at Cornell University (Lin & Liu, 1998). In the first part of the work, both fluids are considered as incompressible, while the second part will treat air compressibility modelling. The mathematical formulation and the numerical resolution of the governing equations of COBRAS2 are derived and some model-experiment comparisons are shown. In particular, validation tests are performed in order to prove model stability and accuracy. The simulation of the rising of a large air bubble in an otherwise quiescent water pool reveals the model capability to reproduce the process physics in a realistic way. Analytical solutions for stationary and internal waves are compared with corresponding numerical results, in order to test processes involving wide range of density difference. Waves induced by dam-break in different scenarios (on dry and wet beds, as well as on a ramp) are studied, focusing on the role of air as the medium in which the water wave propagates and on the numerical representation of bubble dynamics. Simulations of solitary and regular waves, characterized by both spilling and plunging breakers, are analyzed with comparisons with experimental data and other numerical model in order to investigate air influence on wave breaking mechanisms and underline model capability and accuracy. Finally, modelling of air compressibility is included in the new developed model and is validated, revealing an accurate reproduction of processes. Some preliminary tests on wave impact on vertical walls are performed: since air flow modelling allows to have a more realistic reproduction of breaking wave propagation, the dependence of wave breaker shapes and aeration characteristics on impact pressure values is studied and, on the basis of a qualitative comparison with experimental observations, the numerical simulations achieve good results.