Produção Científica
Artigo em Revista
| A Numerical Comparison Between Quasi-MonteCarlo and Sparse Grid Stochastic Collocation Methods Quasi-Monte Carlo methods and stochastic collocation methods based on sparse grids have become popular with solving stochastic partial differential equations.These methods use deterministic points for multi-dimensional integration or interpolation without suffering from the curse of dimensionality. It is not evident which method is best, specially on random models of physical phenomena. We numerically study the error of quasi-Monte Carlo and sparse gridmethods in the context of groundwater flow in heterogeneous media. In particular, we consider the dependence of the variance error on the stochastic dimension and the number of samples/collocation points for steady flow problems in which the hydraulic conductivity is a lognormal process. The suitability of each technique is identified in terms of computational cost and error tolerance. |
Artigo em Revista
| Effect of Element Distortion on the Numerical Dispersion of Spectral Element Methods Spectral element methods are well established in the field of wave propagation,in particular because they inherit the flexibility of finite element methods and have low numerical dispersion error. The latter is experimentally acknowledged, but has been theoretically shown only in limited cases, such as Cartesian meshes. It is well known that a finite element mesh can contain distorted elements that generate numerical errors for very large distortions. In the present work, we study the effect of element distortion on the numerical dispersion error and determine the distortion range in which an accurate solution is obtained for a given error tolerance. We also discuss a double-grid calculation of the spectral element matrices that preserves accuracy in deformed geometries. |
Dissertação de Mestrado
| Filtragem adaptativa SVD de volumes sÃsmicos 3D para realçar refletores e estruturas geológicas. Washington Oliveira Martins. Filtragem adaptativa SVD de volumes sÃsmicos 3D para realçar refletores e estruturas geológicas. 2012. Dissertação (Mestrado em GeofÃsica) - Universidade Federal da Bahia, . Orientador: Milton José Porsani. |
Tese de Doutorado
| Dorian Caraballo Ledesma. Deconvolução de dados sÃsmicos de reflexão utilizando mudança de fase do filtro de Wiener-Levinson. 2011. Dorian Caraballo Ledesma. Deconvolução de dados sÃsmicos de reflexão utilizando mudança de fase do filtro de Wiener-Levinson. 2011. Tese (Doutorado em GeofÃsica) - Universidade Federal da Bahia, Coordenação de Aperfeiçoamento de Pessoal de NÃvel Superior. Orientador: Milton José Porsani. |
Dissertação de Mestrado
| Processamento de Dados SÃsmicos com Grandes Afastamentos: Dados Sintéticos e Linha SÃsmica do Campo de Tenerife, Colômbia. Francisco Ortega Gamboa. Processamento de Dados SÃsmicos com Grandes Afastamentos: Dados Sintéticos e Linha SÃsmica do Campo de Tenerife, Colômbia. 2012. Dissertação (Mestrado em GeofÃsica) - Universidade Federal da Bahia, Coordenação de Aperfeiçoamento de Pessoal de NÃvel Superior. Orientador: Amin Bassrei. |
Apresentação
| GêBR: a free seismic processing interface There are many programs for processing seismic data that are freely available and widespread, for example Seismic Un*x, Madagascar, FreeUSP, and SEPlib, among others. All these packages consist of packages of command-line-oriented programs that are designed to be used in sequence; the data conceptually flow in a pipeline through one program after another. Each program is generally controlled by its own set of command-line options. To take full advantage of such a toolkit, the user must have considerable knowledge beyond general geophysical expertise: shell scripting, process submission and management, and batch queue processing, to name a few. While these skills are useful, they should not be a requirement for seismic data processing. A suitable graphical user interface could take care of these computational details, allowing the user to focus on the central problem of processing seismic data. This is particularly important during training courses, where the limited duration of the does not leave time for learning skills that are not essential to the material being taught. A graphical user interface may also boost the uptake of a new program, by making it more accessible to users and allowing its easy integration with other programs available within the same interface. These principles have guided the development of GêBR, a graphical user interface to control commandline programs for seismic processing. It permits users to build complex processing flows from predefined modules known as menus. Menus describing new programs can be easily added to the interface, extending its capabilities. GêBR is also designed to be simple, in the sense that a couple of hours is enough to introduce the core features of the interface, to allow the user to start working with the seismic data. |
Apresentação
| Minimum-delay seismic trace decomposition and SVD filtering for seismic reflection enhancement Spiking deconvolution corrects for the effect of the seismic wavelet, assumed to be minimum delay, by applying an inverse filter to the seismic trace to get an estimate of reflectivity. In order to compensate for propagation and absorption effects one may use time-varying deconvolution where a different inverse filter is computed and applied for each output sample position. We modify this procedure by estimating a minimum-delay wavelet for each time-sample position of the seismic trace. This gives a decomposition of the seismic trace as a sum of minimum-delay wavelets, each multiplied by a reflectivity coefficient. The reflectivity estimation is a single-trace process which is sensitive to non-white noise, and it does not take into account lateral continuity of reflections. We therefore have developed a new data processing strategy by combining it with adaptive SVD filtering. The SVD filtering process is applied to the data in two steps. First, in a sliding spatial window on NMO-corrected CMP or common shot gathers. Next, after local dip estimation and correction, on local patches in common-offset panels. After the SVD filtering, we apply the new reflectivity estimation procedure. The SVD filtering removes noise and improves lateral continuity while the reflectivity estimation increases the high-frequency content in the data and improves vertical resolution. The new data processing strategy was successfully applied to land seismic data from North-east in Brazil. Improvements in data quality are evident in prestack data panels, velocity analysis and the stacked section. |
Apresentação
| Migration velocity analysis using residual diffraction moveout in the pre/post-stack depth domain Diffraction events contain direct information on the medium velocity. In this work, we develop a method for migration velocity improvement and diffraction localization based on a moveout analysis of over or undermigrated diffraction events in the depth domain. The method uses an arbitrary initial velocity model as input. It provides an update to the velocity model and diffraction locations in the depth domain as a result. The algorithm is based on the focusing of remigration velocity rays from incorrectly migrated diffraction curves. These velocity rays are constructed from a ray-tracing like approach applied to the image-wave equation for velocity continuation. After picking the diffraction events in the migrated domain, the method has a very low computational cost, and the diffraction points are located automatically. We demonstrate the feasibility of our methods using a synthetic data example. |
Apresentação
| Experimental relations between stress and fracture properties on synthetic anisotropic media Elastic anisotropy due to aligned cracks has been the subject of many seismic physical modeling experiments. In earlier investigations, different experimental approaches have taken into account the size, shape and density of cracks. In this paper we present a physical study of the aspect ratio as a function of applied uniaxial stress. We carried out pulse transmission measurements of P- and S-wave velocities in a reference model without inclusions and in a model with penny-shaped neoprene inclusions. The reference model is an anisotropic matrix that consists of stacked plexiglass plates. Rubber discs were used as inclusions in that anisotropic matrix leading to secondary anisotropy. We recorded ultrasonic seismic data using P-wave transducers with central frequency 120 kHz and S-wave transducers with 90 kHz. We compressed the physical models using pressures ranging from 3 to 15.8 MPa. Full crack closure occurs at stress 14.6 MPa normal to model faces. Our analysis indicates three different regimens for the behavior of the inclusions. These results suggest a different dependence of the crack aspect ratio on uniaxial stress at the low state of stress than usually described in the literature. Though our results are not extensive, they show that simple experimental approaches might provide valuable insight into the behavior of cracked rocks at reservoir stress levels. |
Apresentação
| Estimation of fracture orientation through elastic ultrasonic waves Estimate the preferential fracture orientation based on an analysis of cross-correlated S-wave seismograms and Thomsen parameter g is the main goal of this work. For this purpose, we analyzed ultrasonic measurements of elastic S-waves in a physical-modeling experiment with an artificially anisotropic cracked model. The solid matrix of the model consisted of epoxy-resin; small rubber strips simulate cracks with a compliant fill. The anisotropic cracked model has three regions each with a different fracture orientation. We used a rotation of the S-wave polarizations for a cross-correlation analysis of the orientations, and S-wave measurements to evaluate the weak anisotropic parameter g . The shear-wave source had a dominant frequency of 90 kHz. This frequency corresponds to long wavelengths compared to the crack aperture, ensuring effective-media behavior. Integrating the results from crosscorrelation with anisotropic parameter analysis, we were able to estimate the fracture orientation in our anisotropic cracked physical model. The g parameter has shown good agreement with the cross-correlation analysis and, beyond that, provided additional information about the crack orientation that crosscorrelation alone did not fully resolve. Moreover, our results show that the shear waves are strongly influenced by crack orientation. |
<< < 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 > >>
