Produção Científica
Artigo em Revista
Interval mineral and fluid densities estimation from welllogs: Application to the Norne Field dataset Formation evaluation techniques are the key to understand subsurface rocks properties from welllogs, especially those drilled in hydrocarbon exploration wells. Knowledge of the parameters related to different types of rocks is traditionally used in forward determination of lithology, porosity and water saturation, which can be refined by calibrating the input models. In this work, we perform welllog interval linear inversion with respect to formation density to investigate mineral and fluid properties in a real dataset. The method is based on an overdetermined problem, which supposes a homogeneous distribution of petrophysical parameters through stratigraphic layers and is applied in conventional reservoir rocks from the Norne Field (offshore Norway). Bulk density, gammaray and neutron porosity logs are employed in a workflow that relies on layerbylayer leastsquares regressions to estimate matrix, shale and fluid apparent densities. In this process, shale volume and the total density porosity are calculated depth by depth from empirical equations feed by the input logs. Therewith, the introduced inversion scheme stands as an alternative approach for well data interpretation that focuses on computing the densities of the rock constituents instead of fixing these parameters to invert fractional volumes. Furthermore, the application in two wellbores resulted in geological consistent individual densities in most intervals, except for a gasbearing zone observed in one of the boreholes, where porosity uncertainty caused anomalous variationin grain and fluid densities. 

Artigo em Revista
Seismic Stratigraphy of Camamu Basin, Northeastern Brazil â€”The Camamu Basin is located at the northeastern Brazilian coastline and has significant hydrocarbon potential in both shallow and deep water settings. However, despite an already operating productive gas field, the basin is not well known. Herein, a regional stratigraphic interpretation of the rift, the transitional and the drift megasequences is reported, based on a data set of 152 poststacked 2D seismic lines, 1 3D seismic cube and 34 wells with eletrologs. The study revealed that the rift megasequence is much more complex than previously thought. In the southern region of the Camamu Basin, five rift sequences were mapped, while in the northern region, in the area of the regional tectonic lineament known as the Salvador Transcurrent Zone (STZ), seven rift sequences were recognized. This difference suggests a tectonic control during the rifting process, because the northern region is intensely affected by shear stress induced by transcurrent tectonics of the STZ during the crustal breakup. For the postrift or â€˜â€˜transitionalâ€™â€™ phase associated with thermal subsidence installed after the rift, as well as for the drift succession, the tectonic control exerted by Salvadorâ€™s Transcurrent Zone is not detectable and the tectonosedimentary evolution of the basin follows the general pattern of the Brazilian marginal basins. 
Artigo em Revista
NonHyperbolic Velocity Analysis of Seismic Data from Jequitinhonha Basin, Northeastern Brazil Normal moveout (NMO) velocity is used in seismic data processing to correct the data from the moveout effect. This velocity depends on the medium above the reflector and it is estimated from the adjustment of a hyperbolic function that approximates the reflection time. This approximation is reasonable for media formed by isotropic layers. For deeper exploration targets, which effectively behave as anisotropic media,the NMO velocity estimate from the hyperbolic approximation becomes imprecise. One possibility is the use of nonhyperbolic approximations for the reflection time and deeming the medium to be anisotropic. However, these approximations make the NMO velocity estimation a more complex problem, since the anisotropic parameters are unknown. In this study the NMO velocities for a vertical transverse isotropy medium are estimated using two nonhyperbolic reflection time approaches. For comparing the two methodologies that estimate NMO velocity, a 2D dataset from Jequitinhonha Basin is used and it presents anisotropic behavior. The results show that this approach produces more consistent results than the conventional approach, which ignores the anisotropy of the medium. 

Artigo em Revista
Inversion of satellite gravimetric data from RecÃ´ncavoTucanoJatobÃ¡ Basin System Density differences among subsurface rocks cause variations in the gravitational field of Earth, which is known as gravity anomaly. Interpretation of these gravity anomalies allows assessment of the probable depth and shape of the causative body. For several decades, gravity data were acquired on the surface, but after the scientific and technological advances of the last decades, geopotential models were developed, including gravitational observations on a global scale through space satellite missions. This paper investigated the Moho structure in the region of RecÃ´ncavoTucanoJatobÃ¡ riftbasin system based on the information of the terrestrial gravity field from the EIGEN6C4 geopotential model. The frequency domain inversion technique was applied, which is known as the ParkerOldenburg iterative method. Bouguer anomaly data were used in the inversion procedure to determine the thickness and geometry of the crust in the region. Data inversion considered a twolayer model with constant density contrast, in which the entire signal was related to Moho topography. In addition, data inversion was carried out to determine the basement depths. The program proved to be efficient and able to manage large data sets. The results, both of the crust thickness and the sedimentary package, validated the geodynamic evolution understanding of the basin system. 

Artigo em Revista
EVOLUTIONARY PROBLEMS OF NONLINEAR MAGNETOELASTICITY We consider mixed problems for nonlinear equations of magnetoelasticity. Our main result in the threedimensional case is the proof of an existence and uniqueness theorem; uniqueness is established under some extra restrictions on the smoothness of solutions. We also manage to prove the existence and uniqueness of a weak solution to the problem in the twodimensional case; uniqueness is established without any additional a priori assumptions on the smoothness of solutions. 

Artigo em Revista
Petrofacies classification using machine learning algorithms Carbonate reservoirs represent a large portion of the worldâ€™s oil and gas reserves, exhibiting specific characteristics that pose complex challenges to the reservoirsâ€™ characterization, production, and management. Therefore, the evaluation of the relationships between the key parameters, such as porosity,permeability, water saturation, and pore size distribution, is a complex task considering only welllog data, due to the geologic heterogeneity. Hence, the petrophysical parameters are the key to assess the original composition and postsedimentological aspects of the carbonate reservoirs. The concept of reservoir petrofacies was proposed as a tool for the characterization and prediction of the reservoir quality as it combines primary textural analysis with laboratory measurements of porosity, permeability, capillary pressure, photomicrograph descriptions, and other techniques, which contributes to understanding the postdiagenetic events. We have adopted a workflow to petrofacies classification of a carbonate reservoir from the Campos Basin in southeastern Brazil, using the following machine learning methods: decision tree, random forest, gradient boosting, Knearest neighbors, and naÃ¯ve Bayes. The data set comprised1477 wireline data from two wells (A3 and A10) that had petrofacies classes already assigned based on core descriptions. It was divided into two subsets, one for training and one for testing the capability of the trained models to assign petrofacies. The supervisedlearning models have used labeled training data to learn the relationships between the input measurements and the petrofacies to be assigned. Additionally, we have developed a comparison of the modelsâ€™ performance using the testing set according to accuracy, precision, recall, and F1score evaluation metrics. Our approach has proved to be a valuable ally in petrofacies classification, especially for analyzing a welllogging database with no prior petrophysical information. 

Artigo em Revista
Mathematical Model of Water Alternated Polymer Injection Chemical enhanced oil recovery (EOR) methods include the injection of aqueous polymer solutions slugs driven by water. Polymer solutions increase water viscosity, decreasing the water phase mobility and improving oil recovery through better sweep efciency. In this paper, we present the water alternated polymer EOR technique, which is based on the injection of successive polymer slugs alternated by water slugs. The mathematical problem is composed by two conservation equations: one of them is related to the water volume and the other one to the polymer mass. We assume that the polymer may be adsorbed by the rock, and the relation between the concentration in the aqueous solution and the solid is governed by a Langmuir type adsorption isotherm. The water viscosity is a function of the polymer concentration in water. The 2Ã—2 system of hyperbolic equations was decoupled by introducing a potential function instead of time as an independent variable. The water alternated polymer injection is represented by a varying boundary condition. The analytical solution presents interactions between waves of diferent families. It is shown that the polymer slugs always catch up each other along the porous media generating a single slug. As a consequence, the water slugs will disappear. This solution is new and was compared to numerical results with close agreement. It also can be used for the selection of the most suitable enhanced oil recovery technique for a particular oil field. 

Artigo em Revista
The effectiveness of spectral decompositionbased layer thickness estimation: A seismic physical modeling example We have constructed a channel complex model at a scale of 1:10,000 by stacking 3Dprinted polylactide layers with negative relief meandering channels. This model was subjected to an ultrasonic commonoffset acquisition in a water tank (with the water filling the channels), and the result was treated as a zerooffset 3D acoustic reflection seismogram, receiving a deterministic deconvolution and a poststack migration as data treatment. We then developed an algorithm to yield volumes of estimated twoway time layer thickness from multiplefrequency volumes obtained through the shorttime Fourier transform. The estimated thicknesses were compared with the measurements of the physical model obtained through Xray computed tomography. Despite the strong signal attenuation and imaging issues, the results were rather satisfactory, increasing the confidence in using spectral decomposition for quantitative seismic analysis. 

Artigo em Revista
Core plug and 2D/3Dimage integrated analysis for improving permeability estimation based on the differences between micro and macroporosity in Middle East carbonate rocks Carbonate rocks are porous systems, with pores and pore throats of varying morphologies that result from depositional and diagenetic processes. Such heterogeneity produces a complex arrangement between the grains and pores that affects the petrophysical properties while limiting the utility of measurement techniques. Petrophysical properties are generally acquired by conventional laboratory methods, although to provide accurate results, core plugs need to be recovered intact. Twodimensional digital image analysis (2D DIA) enables the processing of any core cut and requires minimal data manipulation and computation when compared to threedimensional approaches. In DIA, permeability is calculated using models described in the literature that often do not provide good predictions for carbonate rocks. Often, the permeability, which is controlled by the size and shape of the pores and pore throats, is related to porosity values; however, the porosity of a rock varies from the micro to macro level, resulting in enormous uncertainty in estimating permeability. In this article, we present a new strategy to improve the prediction of permeability by using poreshape parameters from 2D DIA, which provides data related to the macropores resulting from the optical resolution. This gas technique measures the absolute permeability, which is used as a calibration parameter, and the total porosity, which is used to calculate the microporosity. The test samples used are from Oman outcrops of the Huqf Supergroup and Salalah Formation, which are analogous to the carbonates of the giant reservoirs in the Middle East. Microporosity was characteristic of all the samples due to the calcite mud matrix, recrystallized calcite cement, microcracks and the crystalline texture caused by dolomitization. The poreshape parameters from the 2D DIA improved the permeability prediction and were found to relate to the pore types that make up the rock, whereas the 3D technique did not provide a good result. The R2 of the 2D data was 0.96, demonstrating the efficiency of the procedures applied in mitigating the uncertainties of the models for the set of samples studied. 

Artigo em Revista
Dynamics of the nonequilibrium flow in a duct with obstruction This work aims to numerically simulate the dynamics of a channel flow with an obstruction since the moment we inject a fluid with an homogeneous velocity profile. The simulations uses the open source tools of the OpenFOAM platform, the pisoFoam and the LES turbulence model, describing in detail the velocity profiles of laminar and turbulent flows. We also perform a boundary layer mapping in the presence of an obstacle. We used three different domains to follow the evolution of the velocity profile while the fluid progresses downstream and passes the obstruction. The results reproduce the wellknown results of laminar flow in a channel, as well as the average velocity profile in the turbulent regime and the occurrence of attachments by the obstruction. These preliminary results are used to validate the solvers and the mesh used. Next,an analysis of the velocity profile dynamics resulted in determining an exponential decay of the root mean square deviations of the homogeneous to the parabolic, and to the turbulent regime in the channel. 
