Migração FFD 3D em profundidade usando aproximação de Padé complexa

Implementações dos métodos de migração diferença finita e Fourier (FFD) usam fatoração direcional para acelerar a performance e economizar custo computacional. Entretanto essa técnica introduz anisotropia numérica que podem erroneamente posicionar os refletores em mergulho ao longo das direções em que o não foi aplicado a fatoração no operador de migração. Implementamos a migração FFD 3D, sem usar a técnica do fatoração direcional, no domínio da frequência usando aproximação de Padé complexa. Essa aproximação elimina a anisotropia numérica ao preço de maior custo computacional buscando a solução do campo de onda para um sistema linear de banda larga. Experimentos numéricos, tanto no modelo homogêneo e heterogêneo, mostram que a técnica da fatoração direcional produz notáveis erros de posicionamento dos refletores em meios com forte variação lateral de velocidade. Comparamos a performance de resolução do algoritmo de FFD usando o método iterativo gradiente biconjugado estabilizado (BICGSTAB) e o multifrontal massively parallel direct solver (MUMPS). Mostrando que a aproximação de Padé complexa é um eficiente precondicionador para o BICGSTAB, reduzindo o número de iterações em relação a aproximação de Padé real. O método iterativo BICGSTAB é mais eficiente que o método direto MUMPS, quando usamos apenas um termo da expansão de Padé complexa. Para maior ângulo de abertura do operador, mais termos da série são requeridos no operador de migração, e neste caso, a performance do método direto é mais eficiente. A validação do algoritmo e as propriedades da evolução computacional foram avaliadas para a resposta ao impulso do modelo de sal SEG/EAGE.

Migração em profundidade usando a solução numérica da equação da eiconal

Nos últimos anos tem-se verificado um interesse crescente no desenvolvimento de algoritmos de imageamento sísmico com a finalidade de obter uma imagem da subsuperfície da terra. A migração pelo método de Kirchhoff, por exemplo, é um método de imageamento muito eficiente empregado na busca da localização de refletores na subsuperficie, quando dispomos do cálculo dos tempos de trânsito necessários para a etapa de empilhamento, sendo estes obtidos neste trabalho através da solução da equação eiconal. Primeiramente, é apresentada a teoria da migração de Kirchhoff em profundidade baseada na teoria do raio, sendo em seguida introduzida a equação eiconal, através da qual são obtidos os tempos de trânsitos empregados no empilhamento das curvas de difrações. Em seguida é desenvolvido um algoritmo de migração em profundidade fazendo uso dos tempos de trânsito obtidos através da equação eiconal. Finalmente, aplicamos este algoritmo a dados sintéticos contendo ruído aditivo e múltiplas e obtemos como resultado uma seção sísmica na profundidade. Através dos experimentos feitos neste trabalho observou-se que o algoritmo de migração desenvolvido mostrou-se bastante eficiente e eficaz na reconstrução da imagem dos refletores.

The use of instantaneous phase for improving sparse arrays images

Sparse arrays have pitch larger than half-wavelength (lambda/2) and there is a reduced number of elements in comparison with a full-populated array. Consequently, there is a reduction in cost, data acquisition and processing. However, conventional beamforming techniques result in large side and grating lobes, and consequently in image artifacts. In this work the instantaneous phase of the signals is used in a beamforming technique instead of the instantaneous amplitudes to improve images obtained from sparse arrays configurations. A threshold based on a statistical analysis and the number of signals used for imaging is applied to each pixel, in order to determine if that pixel is related to a defect or not. Three sets of data are used to evaluate the technique, considering medical and non-destructive testing: a simulated point spread function, a medical phantom and an aluminum plate with 2 lambda-, 7 lambda- and lambda-pitch, respectively. The conventional amplitude image is superposed by the image improved by the instantaneous phase, increasing the reflectors detectability and reducing artifacts for all cases, as well as dead zone for the tested plate.

Instantaneous Phase Threshold for Reflector Detection in Ultrasonic Images

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

Interpretação tectono-estratigráfica do reservatório Quissamã (Eo-meso-albiano) nos campos de Pampo, Bicudo e Bonito, Bacia de Campos

Hydrocarbon accumulations occur in marine carbonate reservoirs of Quissamã Formation (early to midlle Albian), southwest Campos Basin. We investigated Pampo, Bicudo and Bonito fields, in order to understand the tectono-structural framework of oncolite/bioclast coarse-grained calcarenite reservoir and the calcilutite-marl-shale sealing interval (Late Albian to Turonian). The database of eleven wells from those fields allowed to elaborate structural sections correlating the Macaé Group – both Quissamã reservoir and Outeiro seal, the latter corresponding to the tectonic deepening phase of basin evolution. Based on density and electric logs, it was prepared structural sections of the carbonate reservoirs with consequent identification of porous zones and oil-water contacts. An extensive 3-D seismic database (~300 Km2) allowed to map three reflectors which represent the limiting units of Macaé Group

Work conditions and occupational health of dentists in Brazilian public health system

Conselho Nacional de Desenvolvimento Ciêntífico e Tecnológico (CNPq)

Sistemas petrolíferos especulativos da bacia de Pelotas (offshore do Uruguai)

Pós-graduação em Geociências e Meio Ambiente - IGCE

Radiative transfer analyses of Titan's tropical atmosphere

Titan's optical and near-IR spectra result primarily from the scattering of sunlight by haze and its absorption by methane. With a column abundance of 92 km amagat (11 times that of Earth), Titan's atmosphere is optically thick and only similar to 10% of the incident solar radiation reaches the surface, compared to 57% on Earth. Such a formidable atmosphere obstructs investigations of the moon's lower troposphere and surface, which are highly sensitive to the radiative transfer treatment of methane absorption and haze scattering. The absorption and scattering characteristics of Titan's atmosphere have been constrained by the Huygens Probe Descent Imager/Spectral Radiometer (DISR) experiment for conditions at the probe landing site (Tomasko, M.G., Bezard, B., Doose, L., Engel, S., Karkoschka, E. 120084 Planet. Space Sci. 56, 624-247: Tomasko, M.G. et al. [2008b] Planet. Space Sci. 56, 669-707). Cassini's Visual and Infrared Mapping Spectrometer (VIMS) data indicate that the rest of the atmosphere (except for the polar regions) can be understood with small perturbations in the high haze structure determined at the landing site (Penteado, P.F., Griffith, CA., Tomasko, M.G., Engel, S., See, C., Doose, L, Baines, K.H., Brown, R.H., Buratti, B.J., Clark, R., Nicholson, P., Sotin, C. [2010]. Icarus 206, 352-365). However the in situ measurements were analyzed with a doubling and adding radiative transfer calculation that differs considerably from the discrete ordinates codes used to interpret remote data from Cassini and ground-based measurements. In addition, the calibration of the VIMS data with respect to the DISR data has not yet been tested. Here, VIMS data of the probe landing site are analyzed with the DISR radiative transfer method and the faster discrete ordinates radiative transfer calculation; both models are consistent (to within 0.3%) and reproduce the scattering and absorption characteristics derived from in situ measurements. Constraints on the atmospheric opacity at wavelengths outside those measured by DISR, that is from 1.6 to 5.0 mu m, are derived using clouds as diffuse reflectors in order to derive Titan's surface albedo to within a few percent error and cloud altitudes to within 5 km error. VIMS spectra of Titan at 2.6-3.2 mu m indicate not only spectral features due to CH4 and CH3D (Rannou, P., Cours, T., Le Mouelic, S., Rodriguez, S., Sotin, C., Drossart, P., Brown, R. [2010]. Icarus 208, 850-867), but also a fairly uniform absorption of unknown source, equivalent to the effects of a darkening of the haze to a single scattering albedo of 0.63 +/- 0.05. Titan's 4.8 mu m spectrum point to a haze optical depth of 0.2 at that wavelength. Cloud spectra at 2 mu m indicate that the far wings of the Voigt profile extend 460 cm(-1) from methane line centers. This paper releases the doubling and adding radiative transfer code developed by the DISR team, so that future studies of Titan's atmosphere and surface are consistent with the findings by the Huygens Probe. We derive the surface albedo at eight spectral regions of the 8 x 12 km(2) area surrounding the Huygens landing site. Within the 0.4-1.6 mu m spectral region our surface albedos match DISR measurements, indicating that DISR and VIMS measurements are consistently calibrated. These values together with albedos at longer 1.9-5.0 mu m wavelengths, not sampled by DISR, resemble a dark version of the spectrum of Ganymede's icy leading hemisphere. The eight surface albedos of the landing site are consistent with, but not deterministic of, exposed water ice with dark impurities. (C) 2011 Elsevier Inc. All rights reserved.

Efficient 1535 nm light emission from an all-Si-based optical micro-cavity containing Er3+ and Yb3+ ions

This work reports on the construction and spectroscopic analyses of optical micro-cavities (OMCs) that efficiently emit at ~1535 nm. The emission wavelength matches the third transmission window of commercial optical fibers and the OMCs were entirely based on silicon. The sputtering deposition method was adopted in the preparation of the OMCs, which comprised two Bragg reflectors and one spacer layer made of either Er- or ErYb-doped amorphous silicon nitride. The luminescence signal extracted from the OMCs originated from the 4I13/2→4I15/2 transition (due to Er3+ ions) and its intensity showed to be highly dependent on the presence of Yb3+ ions.According to the results, the Er3+-related light emission was improved by a factor of 48 when combined with Yb3+ ions and inserted in the spacer layer of the OMC. The results also showed the effectiveness of the present experimental approach in producing Si-based light-emitting structures in which the main characteristics are: (a) compatibility with the actual microelectronics industry, (b) the deposition of optical quality layers with accurate composition control, and (c) no need of uncommon elements-compounds nor extensive thermal treatments. Along with the fundamental characteristics of the OMCs, this work also discusses the impact of the Er3+-Yb3+ ion interaction on the emission intensity as well as the potential of the present findings.

Sea level measurement using single GNSS antenna SNR signals

This thesis presents a possible method to calculate sea level variation using geodetic-quality Global Navigate Satellite System (GNSS) receivers. Three antennas are used: two small antennas and a choke ring one, analyzing only Global Positioning System signals. The main goal of the thesis is to test a modified configuration for antenna set up. In particular, measurements obtained tilting one antenna to face the horizon are compared to measurements obtained from antennas looking upward. The location of the experiment is a coastal environment nearby the Onsala Space Observatory in Sweden. Sea level variations are obtained using periodogram analysis of the SNR signal and compared to synthetic gauge generated from two independent tide gauges. The choke ring antenna provides poor result, with an RMS around 6 cm and a correlation coefficients of 0.89. The smaller antennas provide correlation coefficients around 0.93. The antenna pointing upward present an RMS of 4.3 cm and the one pointing the horizon an RMS of 6.7 cm. Notable variation in the statistical parameters is found when modifying the length of the interval analyzed. In particular, doubts are risen on the reliability of certain scattered data. No relation is found between the accuracy of the method and weather conditions. Possible methods to enhance the available data are investigated, and correlation coefficient above 0.97 can be obtained with small antennas when sacrificing data points. Hence, the results provide evidence of the suitability of SNR signal analysis for sea level variation in coastal environment even in the case of adverse weather conditions. In particular, tilted configurations provides comparable result with upward looking geodetic antennas. A SNR signal simulator is also tested to investigate its performance and usability. Various configuration are analyzed in combination with the periodogram procedure used to calculate the height of reflectors. Consistency between the data calculated and those received is found, and the overall accuracy of the height calculation program is found to be around 5 mm for input height below 5 m. The procedure is thus found to be suitable to analyze the data provided by the GNSS antennas at Onsala.

Kolloide und ihre Überstrukturen als Bausteine zur Herstellung funktioneller Materialien

In dieser Arbeit wird die Synthese von Polymerkolloiden mit unterschiedlichen Formen und Funktionalitäten sowie deren Verwendung zur Herstellung kolloidaler Überstrukturen beschrieben. Über emulgatorfreie Emulsionspolymerisation (SFEP) erzeugte monodisperse sphärische Kolloide dienen als Bausteine von Polymeropalen, die durch die Selbstorganisation dieser Kolloide über vertikale Kristallisation (mit Hilfe einer Ziehmaschine) oder horizontale Kristallisation (durch Aufschleudern oder Aufpipettieren) entstehen. Durch die Kontrolle der Kugelgröße über die Parameter der Emulsionspolymerisation sowie die Einstellung der Schichtdicke der Kolloidkristalle über die Anpassung der Kristallisationsparameter ist die Erzeugung von qualitativ hochwertigen Opalen mit definierter Reflektionswellenlänge möglich. Darüber hinaus kann die chemische und thermische Beständigkeit der Opale durch den Einbau von Vernetzern oder vernetzbaren Gruppen in die Polymere erhöht werden. Die Opalfilme können als wellenlängenselektive Reflektoren in auf Fluoreszenzkonzentratoren basierenden Solarzellensystemen eingesetzt werden, um Lichtverluste in diesen Systemen zu reduzieren. Sie können auch als Template für die Herstellung invertierter Opale aus verschiedenen anorganischen Oxiden (TiO2, Al2O3, ZnO) dienen. Über einen CVD-Prozess erzeugte ZnO-Replika besitzen dabei den Vorteil, dass sie nicht nur eine hohe optische Qualität sondern auch eine elektrische Leitfähigkeit aufweisen. Dies ermöglicht sowohl deren Einsatz als Zwischenreflektor in Tandemsolarzellen als auch die Herstellung hierarchischer Strukturen über die Elektroabscheidung von Nanokristallen. In einem weiteren Teil der Arbeit wird die Herstellung funktioneller formanisotroper Partikel behandelt. Durch die Entmischung von mit Monomer gequollenen vernetzten Partikeln in einer Saatpolymerisation sind mehrere Mikrometer große Kolloide zugänglich, die aus zwei interpenetrierenden Halbkugeln aus gleichen oder verschiedenen Polymeren bestehen. Dadurch sind unter anderem Glycidyl-, Alkin- und Carbonsäuregruppen in die eine oder die andere Halbkugel integrierbar. Diese funktionellen Gruppen erlauben die Markierung bestimmter Partikelhälften mit Farbstoffen, die Beschichtung von Partikelbereichen mit anorganischen Oxiden wie SiO2 sowie die Erzeugung amphiphiler formanisotroper Partikel, die sich an Grenzflächen ausrichten lassen. Das Synthesekonzept kann - ausgehend von mittels SFEP erzeugten stark vernetzten PMMA-Partikeln - auch auf kleine Kolloide mit Größen von mehreren hundert Nanometern übertragen werden.

Geophysical Investigation of Subsurface Three-Dimensional Architecture of Icy Debris Fans With Ground Penetrating Radar in the Wrangell Mountains, Alaska

Icy debris fans have are newly-described landforms (Kochel and Trop, 2008 and 2012) as landforms developed immediately after deglaciation on Earth and similar features have been observed on Mars. Subsurface characteristics of Icy debris fans have not been previously investigated. Ground penetrating radar (GPR) was used to non-invasively investigate the subsurface characteristics of icy debris fans near McCarthy, Alaska, USA. The three fans investigated in Alaska are the East, West, and Middle fans (Kochel and Trop, 2008 and 2012) which below the Nabesna ice cap and on top of the McCarthy Creek Glacier. Icy debris fans in general are a largely unexplored suite of paraglacial landforms and processes in alpine regions. Recent field studies focused on direct observations and depositional processes. Their results showed that the fan's composition is primarily influenced by the type and frequency of depositional processes that supply the fan. Photographic studies show that the East Fan receives far more ice and snow avalanches whereas the Middle and West Fans receive fewer mass wasting events but more clastic debris is deposited on the Middle and West fan from rock falls and icy debris flows. GPR profiles and Wide-angle reflection and refraction (WARR) surveys consisting of both, common mid-point (CMP), and common shot-point (CSP) surveys investigated the subsurface geometry of the fans and the McCarthy Creek Glacier. All GPR surveys were collected in July of 2013 with 100MHz bi-static antennas. Four axial profiles and three cross-fan profiles were done on the West and Middle fans as well as the McCarthy Creek Glacier in order to investigate the relationship between the three features. GPR profiles yielded reflectors that were continuous for 10+ m and hyperbolic reflections in the subsurface. The depth to these reflections in the subsurface requires knowledge of the velocity of the subsurface. To find the velocity of the subsurface eight WARR surveys collected on the fans and on the McCarthy Creek glacier to provide information on variability of subsurface velocities. The profiles of the Middle and West fan have more reflections in their profiles compared to profiles done on the McCarthy Creek Glacier. Based on the WARR surveys, we interpret the lower energy return in the glacier to be caused by two reasons. 1) The increased attenuation due to wet ice versus drier ice and on the fan with GPR velocities >0.15m/ns. 2) Lack of interfaces in the glacier compared to those in the fans which are inferred to be produced by the alternating layers of stratified ice and lithic-rich layers. The GPR profiles on the West and Middle Fans show the shallow subsurface being dominated by lenticular reflections interpreted to be consistent with the shape of surficial deposits. The West Fan is distinguished from the Middle Fan by the nature of its reflections patterns and thicknesses of reflection packages that clearly shows the Middle fan with a greater thickness. The changes in subsurface reflections between the Middle and West Fans as well as the McCarthy Creek Glacier are thought to reflect the type and frequency of depositional processes and surrounding bedrock and talus slopes.

Multi-node TOA-DOA cooperative LOS-NLOS localization : enabling high accuracy and reliability

This dissertation investigates high performance cooperative localization in wireless environments based on multi-node time-of-arrival (TOA) and direction-of-arrival (DOA) estimations in line-of-sight (LOS) and non-LOS (NLOS) scenarios. Here, two categories of nodes are assumed: base nodes (BNs) and target nodes (TNs). BNs are equipped with antenna arrays and capable of estimating TOA (range) and DOA (angle). TNs are equipped with Omni-directional antennas and communicate with BNs to allow BNs to localize TNs; thus, the proposed localization is maintained by BNs and TNs cooperation. First, a LOS localization method is proposed, which is based on semi-distributed multi-node TOA-DOA fusion. The proposed technique is applicable to mobile ad-hoc networks (MANETs). We assume LOS is available between BNs and TNs. One BN is selected as the reference BN, and other nodes are localized in the coordinates of the reference BN. Each BN can localize TNs located in its coverage area independently. In addition, a TN might be localized by multiple BNs. High performance localization is attainable via multi-node TOA-DOA fusion. The complexity of the semi-distributed multi-node TOA-DOA fusion is low because the total computational load is distributed across all BNs. To evaluate the localization accuracy of the proposed method, we compare the proposed method with global positioning system (GPS) aided TOA (DOA) fusion, which are applicable to MANETs. The comparison criterion is the localization circular error probability (CEP). The results confirm that the proposed method is suitable for moderate scale MANETs, while GPS-aided TOA fusion is suitable for large scale MANETs. Usually, TOA and DOA of TNs are periodically estimated by BNs. Thus, Kalman filter (KF) is integrated with multi-node TOA-DOA fusion to further improve its performance. The integration of KF and multi-node TOA-DOA fusion is compared with extended-KF (EKF) when it is applied to multiple TOA-DOA estimations made by multiple BNs. The comparison depicts that it is stable (no divergence takes place) and its accuracy is slightly lower than that of the EKF, if the EKF converges. However, the EKF may diverge while the integration of KF and multi-node TOA-DOA fusion does not; thus, the reliability of the proposed method is higher. In addition, the computational complexity of the integration of KF and multi-node TOA-DOA fusion is much lower than that of EKF. In wireless environments, LOS might be obstructed. This degrades the localization reliability. Antenna arrays installed at each BN is incorporated to allow each BN to identify NLOS scenarios independently. Here, a single BN measures the phase difference across two antenna elements using a synchronized bi-receiver system, and maps it into wireless channel’s K-factor. The larger K is, the more likely the channel would be a LOS one. Next, the K-factor is incorporated to identify NLOS scenarios. The performance of this system is characterized in terms of probability of LOS and NLOS identification. The latency of the method is small. Finally, a multi-node NLOS identification and localization method is proposed to improve localization reliability. In this case, multiple BNs engage in the process of NLOS identification, shared reflectors determination and localization, and NLOS TN localization. In NLOS scenarios, when there are three or more shared reflectors, those reflectors are localized via DOA fusion, and then a TN is localized via TOA fusion based on the localization of shared reflectors.

Precise Orbit Determination of Low Earth Satellites at AIUB using GPS and SLR data

An ever increasing number of low Earth orbiting (LEO) satellites is, or will be, equipped with retro-reflectors for Satellite Laser Ranging (SLR) and on-board receivers to collect observations from Global Navigation Satellite Systems (GNSS) such as the Global Positioning Sys- tem (GPS) and the Russian GLONASS and the European Galileo systems in the future. At the Astronomical Insti- tute of the University of Bern (AIUB) LEO precise or- bit determination (POD) using either GPS or SLR data is performed for a wide range of applications for satellites at different altitudes. For this purpose the classical numeri- cal integration techniques, as also used for dynamic orbit determination of satellites at high altitudes, are extended by pseudo-stochastic orbit modeling techniques to effi- ciently cope with potential force model deficiencies for satellites at low altitudes. Accuracies of better than 2 cm may be achieved by pseudo-stochastic orbit modeling for satellites at very low altitudes such as for the GPS-based POD of the Gravity field and steady-state Ocean Circula- tion Explorer (GOCE).

In Situ Monitoring of Free-Phase Gas Accumulation and Release in Peatlands Using Ground Penetrating Radar (GPR)

We tested a set of surface common mid-point (CMP) ground penetrating radar (GPR) surveys combined with elevation rods ( to monitor surface deformation) and gas flux measurements to investigate in-situ biogenic gas dynamics and ebullition events in a northern peatland ( raised bog). The main findings are: ( 1) changes in the two-way travel time from the surface to prominent reflectors allow estimation of average gas contents and evolution of free-phase gas (FPG); ( 2) peat surface deformation and gas flux measurements are strongly consistent with GPR estimated changes in FPG content over time; ( 3) rapid decreases in atmospheric pressure are associated with increased gas flux; and ( 4) single ebullition events can induce releases of methane much larger ( up to 192 g/m(2)) than fluxes reported by others. These results indicate that GPR is a useful tool for assessing the spatial distribution, temporal variation, and volume of biogenic gas deposits in peatlands.