Produção Científica

**Artigo em Revista**

Multi-frequency electromagnetic method for inductive measurement of ground induced polarization and resistivityA geophysical electromagnetic method to inductively measure the ground electrical resistivity and induced polarization has recently been tested. Its basic characteristics involve three major differences from other methods: the two electrical ground parameters are obtained through measuring magnetic field. For this purpose, a transmitterâ€“receiver (T, R) electromagnetic system is used that operates in the frequency domain and consists of a horizontal loop as the transmitter for the perpendicular loops configuration on the ground surface; the measured function is the (T, R) inductive coupling main variation produced due to the presence of the earth, that is the magnetic field radial component; the measurements are conducted at a large number of frequencies (139 in the more advanced prototype), and the measured function is explored in the frequency interval 0.2 Hz to 1 kHz, a much broader frequency range of the induced polarization effect spectrum, than the one conventionally used in field exploration. Three major aspects are emphasized: (1) the existence of a small â€˜main zoneâ€™ interior to a half-space, which is responsible for most of the magnetic energy that the receiver measures on the half-space surface. This permits to substitute the entire half-space by the â€˜main zoneâ€™ and, in a second step, to substitute the â€˜main zoneâ€™ by an equiv- alent homogeneous half-space with the electrical characteristics of such â€˜main zone';(2) the existence of a closed solution for the fields that the (T, R) system generates on the surface of a homogeneous isotropic half-space, which provides exact functions with the two electrical parameters of interest as the variables (the apparent resistivity and relative polarization parameter); (3) the values of the electrical parameters so determined can be attributed to the central point of the â€˜main zoneâ€™. Three-horizontal layers half-space and a conductive sphere in the free-space are discussed as models. Four field surveys are analysed as examples and show a satisfactory performance of the method for detection of on-shore hydrocarbon reservoirs, description of induced reservoir variations and structural features mapping at depths up to 2.5 km. |

**Artigo em Revista**

Effects of torque produced by wake on the maneuverability of a flatfish autonomous underwater vehicleAutonomous underwater vehicles (AUV) are important resources to be used in the oil exploration industry in deep waters as well as a platform for scanning devices used in open sea regions of difficult human access. This work aims to analyze through computer simulations the influence of marine currents on the maneuverability of a flatfish shaped AUV. The 3D realistic scale simulations were performed on the Yemoja supercomputer located at SENAI-CIMATEC and describe the temporal evolution of the torques in the three rotational degrees of freedom - roll,yaw and pitch. The torques were calculated for two different inlet velocities and three angles (yaw) of attack showing a significant gain in the amplitude of these with increasing velocity and pitch being the component with the greatest amplitude of oscillation. |

**Artigo em Revista**

Cross-correlation in a turbulent flow: Analysis of the velocity field using the pDCCA coefficientThe stochastic process of a turbulent flow in a pipeline provides a time series of the velocity field at any point of the domain by solving numerically the Navier-Stokes equation. The turbulent flow was produced by obstacles near the inlet, injecting eddies into the current. Moving downstream, these eddies evolve to a fully turbulent flow. Many length and time scales are involved in this process. We explore the cross-correlations of the velocity field time series at different points and also at different time scales using the detrended cross-correlation coefficient, pDCCA, designed to analyze the cross-correlations in non-stationary time series. Thus, the results with DCCA allow interpreting how these eddies propagate downstream, and also quantify how adherent the velocity fields are with respect to the pipeline position. |

**Artigo em Revista**

Detection of the persistency of the blockages symmetry influence on the multi-scale cross-correlations of the velocity fields in internal turbulent flows in pipelinesIn this paper we analyze the influence of obstacles symmetry on the development of the turbulent flow of a fluid through a pipeline. The analysis is based on the numerical solutions of the Navierâ€“Stokes equations for the velocity field. The influence of the obstacle symmetry on the turbulence is detected by changing their shape while keeping the blockage ratio constant and calculating velocity field cross-correlations on the time series resulting from the simulation. The Detrended Cross-Correlation coefficient (ÏDCCA) is applied to obtain two-point correlations located at different regions of the channel: at mid-channel and near the walls (above and below). With this cross-correlation coefficient we quantify how far from the obstructions these coefficients become independent on the obstructions shape, establishing a scale for the obstruction symmetry memory loss. |

**Artigo em Revista**

A numerical viscoelastic model of ground response assimilating pore-water pressure measurementsWe consider a simple one-dimensional, viscoelastic model for shear-wave propagation on liquefiable soils. The soil is modelled as a layered medium parametrized by shear modulus and viscosity, which in turn depend on the excess pore-water pressure ratio. We numer ically solve the resulting wave propagation model with the spectral element method, and employ simulated annealing and weighted Gauss-Newton inversion algorithms to minimize the misfit of surface displacement, velocity, and acceleration. This procedure is validated us ing recorded ground motion and pore-water pressure data from the Imperial Valley Wildlife and the Garner Valley downhole arrays. Parameter inversion is also carried out with linear models with constant shear modulus and viscosity, and the proposed model provides better fitness in the presence of strong motion, especially in the 1987 Superstition Hills earthquake. Key words: viscoelastic wave equation, liquefaction, spectral element method, site response. |

**Artigo em Revista**

Complex Autoregressive Timeâ€“Frequency Analysis - Estimation of time-varying periodic signal componentsTimeâ€“frequency representations of nonstationary signals have a wide range of geophysical applications, including seismics, seismology, volcanology, and astrophysics. In this article, we estimate a complex autoregressive (AR) model from a short time window of the analytic signal. The local power spectrum is the inverse of the spectrum of this AR model. Since the coefficients are complex, the time window can be shorter than for the real AR model, which requires more coefficients. This results in higher timeâ€“frequency resolution, as seen in a synthetic data examples with different signal components. The new technique also gives good results when computing the instantaneous average frequency (IAF) of marine seismic data. Applied to digitized and downloaded data from the Laser Interferometer Gravitational-Wave Observatory (LIGO) in Hanford, Washington, the result clearly shows the linear chirp associated with the merger of two black holes. |

**Artigo em Revista**

APPLICATION OF TIME-FREQUENCY DECOMPOSITION METHOD IN THE STUDY OF GAS RESERVOIR IN THE SERGIPE-ALAGOAS BASINThe sedimentary basin of Sergipe-Alagoas, located on the Brazilian east bank, presents one of the most complete stratigraphic sections of the Brazilian continental margin. Hydrocarbon exploration activities began more than 50 years ago. The recent discoveries of hydrocarbons (gas and oil of high API grade) in turbiditic reservoirs of deep waters have further awakened the exploratory interest of the basin. Problems related to the processing and interpretation of seismic data have always received great attention from the scientific community. Currently, the use of time-frequency decomposition methods of the seismic signal is of great interest. Spectral decomposition has been widely used in reservoir characterization, such as determination of layer thickness, stratigraphic visualization with seismic attributes and identification of low frequency anomalies associated with the presence of gas. The mechanism causing these anomalies is not yet well known, but they are often attributed to the high attenuation of gas filled reservoirs. The approach used for spectral decomposition combines the maximum entropy method and the Wigner-Ville distribution, based on the idea of the Burg method that uses the prediction error operator to extend the Wigner-Ville kernel sequences by applying the Fourier transform to each extended sequence, thus allowing to obtain the Wigner-Ville distribution of maximum entropy. Keywords: Sergipe-Alagoas Basin,Wigner-Ville distribution, maximum entropy, spectral decomposition, seismic attributes, low frequency anomaly. |

**Artigo em Revista**

Processing of large offset data: experimental seismic line from Tenerife Field, ColombiaExploration seismology provides the main source of information about the Earthâ€™s subsurface, which in many cases can be presented as a simple model of horizontal or near-horizontal layers. After the seismic acquisition step, conventional seismic processing of reflection data provides an image of the subsurface by using information about the reflections of these layers. The traveltime from a source to different receivers is adjusted using a hyperbolic function. This expression is used in the case involving an isotropic medium, which is a simplification of nature, whereas geologically complex media are generally anisotropic. A subsurface model that more closely resembles reality is the vertical transverse isotropy, which defines two parameters that are required to correct the traveltimes: the NMO velocity and the anellipticity parameter. In this paper, we reviewed the literature and methodology for velocity analysis of seismic data acquired from anisotropic media. A model with horizontal layers and anisotropic behavior was developed and evaluated. The anisotropic velocity was compared to the isotropic velocity, and the results were analyzed. Finally, the methodology was applied to real seismic data, i.e. an experimental landline from Tenerife Field, Colombia. The results show the importance of the anellipticity parameter in models with anisotropic layers. |

**Artigo em Revista**

Filtering and frequency interpretations of Singular Spectrum AnalysisNew filtering and spectral interpretations of Singular Spectrum Analysis(SSA)are provided. It is shown that the variables reconstructed fromd iagonal averaging of reduced-rank approximations to the trajec-tory matrix can be obtained from a noncausal convolution filte rwith zero-phase characteristics. There-constructed variables are readily constructed using a two-pass filtering algorithm that is well known in the signal processing literature. When the number of rows in the trajectory matrix is much larger than number of columns, many results reported in the signal processing literature can be used to derive the properties of the resulting filters and their spectra. New features of there constructed series are revealed using these results. Two examples are used to illustrate the results derived in this paper. |

**Artigo em Revista**

OtimizacÌ§aÌƒo global para resolver problemas inversos em eletrorresistividade com flexibilidade na escolha dos viÌnculosInversion in DC-resistivity is an ill-posed inverse problem because different realizations of the same model might satisfy approximately the same data fitting criterium. It is therefore necessary to use constraints to obtain unique and / or stable solutions to small perturbations in the measurements. However, in general, the introduction of constraints has been restricted to cases of differentiable constraints, which can be treated with local optimization algorithms. 1D and 2D modeling in DC-resistivity is computationally inexpensive, allowing the use of global optimization methods (GOMs) to solve 1.5D and 2D inverse problems with flexibility in constraint incorporation. Changes in the cost function, either in the constraints or data fitting criteria, can be easily performed, since each term of the cost function is properly normalized to allow the approximate invariance of the Lagrange multipliers. GOMs have the potential to support a computational environment suitable for quantitative interpretation in which the comparison of solutions incorporating different constraints is one way of inferring characteristics of the actual distribution of the underground resistivity. In this work, we developed: (i) comparison of the performances of the Simulated Annealing (SA), Genetic Algorithm (GA) and Particle Swarm optimization (PSO) methods to solve the 1.5D inverse problem in DC resistivity using synthetic and field data; (ii) an inversion approach based on particle swarm optimization (PSO) to solve the 2D DC-resistivity inverse problem; (iii) exploration of several constraints in the variation of log-resistivity, including spatial continuity in both L1 andL2 norms, total variation and sparsity constraints using discrete cosine and Daubechies bases. In addition, we explore the minimum inertia constraint, including the case of using the Earthâ€™s surface as the target axis, to impose the concentration of resistive or conductive materials along target axes. The main results of the comparison for the 1.5D case are: a) all methods reproduce quite well the resistivity distribution of synthetic models, b) PSO and GA are very robust to changes in the cost function and SA is comparatively much more sensitive, c) PSO first and GA second present the best computational performances, requiring smaller number of forwarding modeling than SA, and d) GA shows the best performance with respect to the final attained value of the cost function and its standard deviation, whilst SA has the worst performance in this aspect. Equally important for both 1.5 and 2D cases, from the stopping criteria of the PSO algorithm results not only the best solution but also a cluster of suboptimal quasi-solutions from which uncertainty analyses can be performed. As a result, the interpreter has freedom to perform a quantitative interpretation process based on a feedback trial-and-error inversion approach, in a similar manner he/she has when using a friendly forward |

<< < 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 68 > >>