Estudo da perda de carga e calor no poço injetor no processo De drenagem gravitacional assistida por vapor (SAGD)

The oil companies in the area in general are looking for new technologies that can increase the recovery factor of oil contained in reservoirs. These investments are mainly aimed at reducing the costs of projects which are high. Steam injection is one of these special methods of recovery in which steam is injected into the reservoir in order to reduce the viscosity of the oil and make it more mobile. The process assisted gravity drainage steam (SAGD) using steam injection in its mechanism, as well as two parallel horizontal wells. In this process steam is injected through the horizontal injection well, then a vapor chamber is formed by heating the oil in the reservoir and, by the action of gravitational forces, this oil is drained down to where the production well. This study aims to analyze the influence of pressure drop and heat along the injection well in the SAGD process. Numerical simulations were performed using the thermal simulator STARS of CMG (Computer Modeling Group). The parameters studied were the thermal conductivity of the formation, the flow of steam injection, the inner diameter of the column, the steam quality and temperature. A factorial design was used to verify the influence of the parameters studied in the recovery factor. We also analyzed different injection flow rates for the model with pressure drop and no pressure drop, as well as different maximum flow rates of oil production. Finally, we performed an economic analysis of the two models in order to check the profitability of the projects studied. The results showed that the pressure drop in injection well have a significant influence on the SAGD process.

Efeito da perda de carga e calor no poço injetor no processo de drenagem gravitacional assistido com vapor e solvente

Nowadays, most of the hydrocarbon reserves in the world are in the form of heavy oil, ultra - heavy or bitumen. For the extraction and production of this resource is required to implement new technologies. One of the promising processes for the recovery of this oil is the Expanding Solvent Steam Assisted Gravity Drainage (ES-SAGD) which uses two parallel horizontal wells, where the injection well is situated vertically above the production well. The completion of the process occurs upon injection of a hydrocarbon additive at low concentration in conjunction with steam. The steam adds heat to reduce the viscosity of the oil and solvent aids in reducing the interfacial tension between oil/ solvent. The main force acting in this process is the gravitational and the heat transfer takes place by conduction, convection and latent heat of steam. In this study was used the discretized wellbore model, where the well is discretized in the same way that the reservoir and each section of the well treated as a block of grid, with interblock connection with the reservoir. This study aims to analyze the influence of the pressure drop and heat along the injection well in the ES-SAGD process. The model used for the study is a homogeneous reservoir, semi synthetic with characteristics of the Brazilian Northeast and numerical simulations were performed using the STARS thermal simulator from CMG (Computer Modelling Group). The operational parameters analyzed were: percentage of solvent injected, the flow of steam injection, vertical distance between the wells and steam quality. All of them were significant in oil recovery factor positively influencing this. The results showed that, for all cases analyzed, the model considers the pressure drop has cumulative production of oil below its respective model that disregards such loss. This difference is more pronounced the lower the value of the flow of steam injection

Modelagem e controle de um microveículo aéreo: uma aplicação de estabilidade robusta com a técnica backstepping em uma estrutura hexarrotor

In this Thesis, the development of the dynamic model of multirotor unmanned aerial vehicle with vertical takeoff and landing characteristics, considering input nonlinearities and a full state robust backstepping controller are presented. The dynamic model is expressed using the Newton-Euler laws, aiming to obtain a better mathematical representation of the mechanical system for system analysis and control design, not only when it is hovering, but also when it is taking-off, or landing, or flying to perform a task. The input nonlinearities are the deadzone and saturation, where the gravitational effect and the inherent physical constrains of the rotors are related and addressed. The experimental multirotor aerial vehicle is equipped with an inertial measurement unit and a sonar sensor, which appropriately provides measurements of attitude and altitude. A real-time attitude estimation scheme based on the extended Kalman filter using quaternions was developed. Then, for robustness analysis, sensors were modeled as the ideal value with addition of an unknown bias and unknown white noise. The bounded robust attitude/altitude controller were derived based on globally uniformly practically asymptotically stable for real systems, that remains globally uniformly asymptotically stable if and only if their solutions are globally uniformly bounded, dealing with convergence and stability into a ball of the state space with non-null radius, under some assumptions. The Lyapunov analysis technique was used to prove the stability of the closed-loop system, compute bounds on control gains and guaranteeing desired bounds on attitude dynamics tracking errors in the presence of measurement disturbances. The controller laws were tested in numerical simulations and in an experimental hexarotor, developed at the UFRN Robotics Laboratory

Utilização de redes neurais artificiais para detecção e diagnóstico de falhas

In a real process, all used resources, whether physical or developed in software, are subject to interruptions or operational commitments. However, in situations in which operate critical systems, any kind of problem may bring big consequences. Knowing this, this paper aims to develop a system capable to detect the presence and indicate the types of failures that may occur in a process. For implementing and testing the proposed methodology, a coupled tank system was used as a study model case. The system should be developed to generate a set of signals that notify the process operator and that may be post-processed, enabling changes in control strategy or control parameters. Due to the damage risks involved with sensors, actuators and amplifiers of the real plant, the data set of the faults will be computationally generated and the results collected from numerical simulations of the process model. The system will be composed by structures with Artificial Neural Networks, trained in offline mode using Matlab®

Pecolação em uma rede multifractal

In this work we present the principal fractals, their caracteristics, properties abd their classification, comparing them to Euclidean Geometry Elements. We show the importance of the Fractal Geometry in the analysis of several elements of our society. We emphasize the importance of an appropriate definition of dimension to these objects, because the definition we presently know doesn t see a satisfactory one. As an instrument to obtain these dimentions we present the Method to count boxes, of Hausdorff- Besicovich and the Scale Method. We also study the Percolation Process in the square lattice, comparing it to percolation in the multifractal subject Qmf, where we observe som differences between these two process. We analize the histogram grafic of the percolating lattices versus the site occupation probability p, and other numerical simulations. And finaly, we show that we can estimate the fractal dimension of the percolation cluster and that the percolatin in a multifractal suport is in the same universality class as standard percolation. We observe that the area of the blocks of Qmf is variable, pc is a function of p which is related to the anisotropy of Qmf

Dinâmica e estrutura de redes complexas no modelo de afinidade

In this work we elaborate and discuss a Complex Network model which presents connectivity scale free probability distribution (power-law degree distribution). In order to do that, we modify the rule of the preferential attachment of the Bianconi-Barabasi model, including a factor which represents the similarity of the sites. The term that corresponds to this similarity is called the affinity, and is obtained by the modulus of the difference between the fitness (or quality) of the sites. This variation in the preferential attachment generates very interesting results, by instance the time evolution of the connectivity, which follows a power-law distribution ki / ( t t0 )fi, where fi indicates the rate to the site gain connections. Certainly this depends on the affinity with other sites. Besides, we will show by numerical simulations results for the average path length and for the clustering coefficient

Efeitos de Confinamento de Vórtice Magnético em Nanodiscos de Fe e Py

We report two theoretical works, based in numerical simulations. The first study consists in the investigation of equilibrium phases and vortex formation in Ferro and Permalloy circular and square nanoelements.The another have the aim to investigate the magnetostatic interaction between pairs of nanodisks of Ferro and Permalloy and it`s impact in the vortex structure

Nucleação de vortices e paredes de domínio em nanoestruturas magnéticas

We report a theoretical investigation of the magnetic phases and hysteresis of exchange biased ferromagnetic (F) nanoelements for three di erent systems: exchange biased nanoparticles, exchange biased narrow ferromagnetic stripes and exchange biased thin ferromagnetic lms. In all cases the focus is on the new e ects produced by suitable patterns of the exchange energy coupling the ferromagnetic nanoelement with a large anisotropy antiferromagnetic (AF) substrate. We investigate the hysteresis of iron and permalloy nanoparticles with a square basis, with lateral dimensions between 45 nm and 120 nm and thickness between 12 nm and 21 nm. Interface bias is aimed at producing large domains in thin lms. Our results show that, contrary to intuition, the interface exchange coupling may generate vortex states along the hysteresis loop. Also, the threshold value of the interface eld strength for vortex nucleation is smaller for iron nanoelements. We investigate the nucleation and depinning of an array of domain walls pinned at interface defects of a vicinal stripe/AF bilayer. The interface exchange eld displays a periodic pattern corresponding to the topology of the AF vicinal substrate. The vicinal AF substrate consists of a sequence of terraces, each with spins from one AF subalattice, alternating one another. As a result the interface eld of neighboring terraces point in opposite direction, leading to the nucleation of a sequence of domain walls in the ferromagnetic stripe. We investigated iron an permalloy micrometric stripes, with width ranging from 100 nm and 300 nm and thickness of 5 nm. We focused in domain wall sequences with same chirality and alternate chirality. We have found that for 100nm terraces the same chiraility sequence is more stable, requiring a larger value of the external eld for depinning. The third system consists of an iron lm with a thickness of 5 nm, exchange coupled to an AF substrate with a periodic distribution of islands where the AF spins have the opposite direction of the spins in the background. This corresponds to a two-sublattice noncompensated AF plane (such as the surface of a (100) FeF2 lm), with monolayer-height islands containing spins of one sublattice on a surface containing spins of the opposite sublattice. The interface eld acting in the ferromagnetic spins over the islands points in the opposite direction of that in the spins over the background. This a model system for the investigation of interface roughness e ects. We have studied the coercicivity an exchange bias hysteresis shift as a function of the distance between the islands and the degree of interface roughness. We have found a relevant reduction of coercivity for nearly compensated interfaces. Also the e ective hysteresis shift is not proportional to the liquid moment of the AF plane. We also developed an analytical model which reproduces qualitatively the results of numerical simulations

Modelagem Matemática e Computacional de fenômenos eletrocinéticos em meios porosos carregados eletricamente

In this work we present a mathematical and computational modeling of electrokinetic phenomena in electrically charged porous medium. We consider the porous medium composed of three different scales (nanoscopic, microscopic and macroscopic). On the microscopic scale the domain is composed by a porous matrix and a solid phase. The pores are filled with an aqueous phase consisting of ionic solutes fully diluted, and the solid matrix consists of electrically charged particles. Initially we present the mathematical model that governs the electrical double layer in order to quantify the electric potential, electric charge density, ion adsorption and chemical adsorption in nanoscopic scale. Then, we derive the microscopic model, where the adsorption of ions due to the electric double layer and the reactions of protonation/ deprotanaç~ao and zeta potential obtained in modeling nanoscopic arise in microscopic scale through interface conditions in the problem of Stokes and Nerst-Planck equations respectively governing the movement of the aqueous solution and transport of ions. We developed the process of upscaling the problem nano/microscopic using the homogenization technique of periodic structures by deducing the macroscopic model with their respectives cell problems for effective parameters of the macroscopic equations. Considering a clayey porous medium consisting of kaolinite clay plates distributed parallel, we rewrite the macroscopic model in a one-dimensional version. Finally, using a sequential algorithm, we discretize the macroscopic model via the finite element method, along with the interactive method of Picard for the nonlinear terms. Numerical simulations on transient regime with variable pH in one-dimensional case are obtained, aiming computational modeling of the electroremediation process of clay soils contaminated

Parâmetro de regularização em problemas inversos: estudo numérico com a transformada de Radon

In general, an inverse problem corresponds to find a value of an element x in a suitable vector space, given a vector y measuring it, in some sense. When we discretize the problem, it usually boils down to solve an equation system f(x) = y, where f : U Rm ! Rn represents the step function in any domain U of the appropriate Rm. As a general rule, we arrive to an ill-posed problem. The resolution of inverse problems has been widely researched along the last decades, because many problems in science and industry consist in determining unknowns that we try to know, by observing its effects under certain indirect measures. Our general subject of this dissertation is the choice of Tykhonov´s regulaziration parameter of a poorly conditioned linear problem, as we are going to discuss on chapter 1 of this dissertation, focusing on the three most popular methods in nowadays literature of the area. Our more specific focus in this dissertation consists in the simulations reported on chapter 2, aiming to compare the performance of the three methods in the recuperation of images measured with the Radon transform, perturbed by the addition of gaussian i.i.d. noise. We choosed a difference operator as regularizer of the problem. The contribution we try to make, in this dissertation, mainly consists on the discussion of numerical simulations we execute, as is exposed in Chapter 2. We understand that the meaning of this dissertation lays much more on the questions which it raises than on saying something definitive about the subject. Partly, for beeing based on numerical experiments with no new mathematical results associated to it, partly for being about numerical experiments made with a single operator. On the other hand, we got some observations which seemed to us interesting on the simulations performed, considered the literature of the area. In special, we highlight observations we resume, at the conclusion of this work, about the different vocations of methods like GCV and L-curve and, also, about the optimal parameters tendency observed in the L-curve method of grouping themselves in a small gap, strongly correlated with the behavior of the generalized singular value decomposition curve of the involved operators, under reasonably broad regularity conditions in the images to be recovered

Desenvolvimento de dobras e falhas em ambiente distensional: Aplicação da modelagem física

The geological modeling allows, at laboratory scaling, the simulation of the geometric and kinematic evolution of geological structures. The importance of the knowledge of these structures grows when we consider their role in the creation of traps or conduits to oil and water. In the present work we simulated the formation of folds and faults in extensional environment, through physical and numerical modeling, using a sandbox apparatus and MOVE2010 software. The physical modeling of structures developed in the hangingwall of a listric fault, showed the formation of active and inactive axial zones. In consonance with the literature, we verified the formation of a rollover between these two axial zones. The crestal collapse of the anticline formed grabens, limited by secondary faults, perpendicular to the extension, with a curvilinear aspect. Adjacent to these faults we registered the formation of transversal folds, parallel to the extension, characterized by a syncline in the fault hangingwall. We also observed drag folds near the faults surfaces, these faults are parallel to the fault surface and presented an anticline in the footwall and a syncline hangingwall. To observe the influence of geometrical variations (dip and width) in the flat of a flat-ramp fault, we made two experimental series, being the first with the flat varying in dip and width and the second maintaining the flat variation in width but horizontal. These experiments developed secondary faults, perpendicular to the extension, that were grouped in three sets: i) antithetic faults with a curvilinear geometry and synthetic faults, with a more rectilinear geometry, both nucleated in the base of sedimentary pile. The normal antithetic faults can rotate, during the extension, presenting a pseudo-inverse kinematics. ii) Faults nucleated at the top of the sedimentary pile. The propagation of these faults is made through coalescence of segments, originating, sometimes, the formation of relay ramps. iii) Reverse faults, are nucleated in the flat-ramp interface. Comparing the two models we verified that the dip of the flat favors a differentiated nucleation of the faults at the two extremities of the mater fault. V These two flat-ramp models also generated an anticline-syncline pair, drag and transversal folds. The anticline was formed above the flat being sub-parallel to the master fault plane, while the syncline was formed in more distal areas of the fault. Due the geometrical variation of these two folds we can define three structural domains. Using the physical experiments as a template, we also made numerical modeling experiments, with flat-ramp faults presenting variation in the flat. Secondary antithetic, synthetic and reverse faults were generated in both models. The numerical modeling formed two folds, and anticline above the flat and a syncline further away of the master fault. The geometric variation of these two folds allowed the definition of three structural domains parallel to the extension. These data reinforce the physical models. The comparisons between natural data of a flat-ramp fault in the Potiguar basin with the data of physical and numerical simulations, showed that, in both cases, the variation of the geometry of the flat produces, variation in the hangingwall geometry